Structure and evolution of teleost mitochondrial control regions

We amplified and sequenced the mitochondrial control region from 23 species representing six families of teleost fish. The length of this segment is highly variable among even closely related species due to the presence of tandemly repeated sequences and large insertions. The position of the repetitive sequences suggests that they arise during replication both near the origin of replication and at the site of termination of the D-loop strand. Many of the conserved sequence blocks (CSBs) observed in mammals are also found among fish. In particular, the mammalian CSB-D is present in all of the fish species studied. Study of potential secondary structures of RNAs from the conserved regions provides little insight into the functional constraints on these regions. The variable structure of these control regions suggests that particular care should be taken to identify the most appropriate segment for studies of intraspecific variation. Correspondence to: T.D. Kocher     

大林姬鼠两亚种线粒体DNA 细胞色素b 基因和控制区的序列多样性

为了研究大林姬鼠两亚种(韩国的指名亚种及中国东北和内蒙古地区的东北亚种)线粒体DNA 的变异程度并确定朝鲜亚种的分类地位,我们分别将来自韩国和中国东北长白山地区的两亚种的线粒体DNA 的细胞色素b 基因和控制区进行了测序分析。我们将测序所得到细胞色素b 基因序列与来自基因库的大林姬鼠5 个亚种的相应的单倍型进行了分析,结果显示大林姬鼠可分为4 个类群[类群1:韩国大林姬鼠指名亚种;类群2:中国长白山和内蒙古地区的东北亚种、俄罗斯外加贝尔的majuculus 亚种;类群3:中国长春的东北亚种、俄罗斯Primorye(俄罗斯远东地区) rufulus 亚种、俄罗斯库页岛(俄罗斯远东地区)和日本北海道地区的giliacus 亚种;类群4:中国黑龙江海林地区的东北亚种]。线粒体的控制区序列分析显示韩国指名亚种也不同于中国东北地区的东北亚种。本研究的类群1,2 和3 与Serizawa et al. (2002)的研究的K、S 和R 的分支相对应。这表明韩国指名亚种(类群1 和分支K)的线粒体DNA 与其他类群不同。另外,我们还发现在细胞色素b 基因构建的系统树中,东北亚种可以与类群2 (分支S)及类群3 (分支R 的不同亚种聚合在一起。我们认为线粒体DNA 的母性遗传与两个相邻亚种的个体之间的种内杂交造成了基于细胞色素b 序列对东北亚种的聚类分析结果与基于形态学特征的分类结果的不一致。因此,我们提出对这些显示出核苷酸序列多样性的东北亚种不能只用细胞色素b的数据进行亚种分类,还应该结合形态学和核DNA 特征进行进一步分析。最后,我们还发现韩国的指名亚种的细胞色素b 序列在平均距离16. 93% 的基础上不同于来自基因库的A. speciosus。Jones and Johnson (1965)指出了韩国的大林姬鼠在形态上的区别,所以我们认为韩国的大林姬鼠指名亚种A. p. peninsulae 是一种具有形态和遗传特异性的地方亚种。     

Morphological, mitochondrial DNA and allozyme evolution in representative amphibians and reptiles inhabiting each side of the Strait of Gibraltar

Patterns of differentiation in morphology, mitochondrial DNA and allozymes in amphibians and reptiles inhabiting northern and southern shores of the Strait of Gibraltar are not concordant, suggesting that each taxon was affected differently by events preceding or following the formation of the Strait of Gibraltar. Mitochondrial DNA and allozyme differentiation between Discoglossus jeanneae and Discoglossus scovazzi (Anura, Discoglossidae), Rana perezi and Rana saharica (Anura, Ranidae), and Blanus cinereus and Blanus tingitanus (Squamata, Amphisbaenia, Amphisbaenidae) is substantial, whereas morphological differentiation is moderate in Rana and Blanus , but is substantial in Discoglossus . Differentiation in mitochondrial DNA and morphology between Timon ( Lacerta ) lepidus and Timon ( Lacerta ) tangitanus (Squamata, Lacertoidea, Lacertidae) is considerable, but allozyme differentiation is low. In members of type-I and -II Podarcis vaucheri (Squamata, Lacertoidea, Lacertidae), morphology and mitochondrial DNA are moderately differentiated, but allozyme differentiation is low. Spanish and Moroccan populations of Hyla meridionalis (Anura, Hylidae), Mauremys leprosa (Testudines, Geoemydidae), and Macroprotodon brevis (Squamata, Serpentes, Colubridae) demonstrate little allozyme and mitochondrial DNA differentiation, but whereas morphological differentiation between Mauremys and Macroprotodon populations is moderate, Hyla demonstrate substantial morphological differentiation between continental populations. These data suggest that sex-limited mitochondrial markers are reflective of ancient phylogenetic history, whereas biparentally inherited allozyme markers and morphological characteristics reflect more recent population structure and movement.  © 2008 The Linnean Society of London, Biological Journal of the Linnean Society , 2008, 94 , 445–461.     

基于3种线粒体标记探讨中日沿海角木叶鲽遗传多样性差异

角木叶鲽(Pleuronichthys cornutus)是东亚沿海重要的鲽形目经济鱼类, 为更好地保护和开发利用其种质资源, 有必要全面了解其遗传背景。本研究测定了中国和日本沿海7个群体200尾角木叶鲽线粒体控制区(CR) 5'端、细胞色素b (Cytb)和NADH脱氢酶第二亚基(ND2)基因序列, 比较不同标记在解析遗传多样性和种群结构上的可行性与有效性, 阐明中日沿海角木叶鲽群体间出现遗传分化的分子机制。CR序列分析发现中日沿海7个角木叶鲽群体遗传多样性表现出较高的单倍型多样性(H_d = 0.9699)和较低的核苷酸多样性(π = 0.0061); 各群体间无显著的遗传分化(F_ST = -0.0197-0.0184, P > 0.05); 单倍型网络未显示出明显的地理聚群和谱系结构; 分子方差分析(AMOVA)表明变异主要发生在群体内部(> 99.17%)。进一步通过Cytb和ND2基因分别与CR序列对比分析, 结果表明群体遗传多样性均表现为高H_d (0.9683-0.9829)低π (0.0050-0.0063)模式, 仅有ND2基因分析F_ST值(F_ST = 0.0302, P < 0.05)显示了中国碣石(GDJS)和日本明石(JAP)群体间显著的低水平遗传分化现象。CR、Cytb和ND2的单倍型网络图均无明显的地理聚类和谱系结构, AMOVA分析也显示变异主要来源于群体内(> 98.39%)。种群历史动态分析结果显示, 角木叶鲽可能在第四纪中更新世晚期经历了群体扩张事件, 扩张时间分别为31.93-9.58万年前(CR)、27.53-22.02万年前(Cytb)和26.99-18.75万年前(ND2)。综上所述, 中日沿海的角木叶鲽具有较高遗传多样性, GDJS和JAP群体间存在低度分化; ND2基因比CR和Cytb序列更适于分析角木叶鲽种群遗传结构, 选择多个遗传标记可有效弥补单一标记分析遗传多样性的局限性; 推测冰期两大独立避难所的形成及GDJS和JAP群体距离相隔较远是其发生遗传分化的主要原因。研究结果为中日沿海角木叶鲽渔业资源的种质保护与可持续利用提供了理论依据。     

翘嘴鳜线粒体DNA细胞色素b基因序列分析

测定了翘嘴鳜（Siniperca chuatsi）的线粒体细胞色素b基因编码区全序列1140个核苷酸序列,并且由此推导出对应的氨基酸序列,A/T所占比例为53.4％,G C约占46.6％。该基因中密码子第1位核苷酸中4种碱基组成较为均衡,第2位核苷酸中T的使用率较高,G的使用率较低;而密码子第3位C/A的使用率高,而G的使用率仅为3.7％;在氨基酸序列中Leu所占比例16.36％,远高于其它氨基酸;而Cys使用比例仅为0.79％。本研究为鳜类的系统发育,资源保护和利用提供了核苷酸序列方面的资料,至于翘嘴鳜线粒体细胞色素b基本序列中密码子不同位置碱基使用比率以及氨基酸组成是翘嘴鳜特有还是鳜类共有特征,尚有待进一步研究。     

Tempo and mode of sequence evolution in mitochondrial DNA of HawaiianDrosophila

Summary Sequence comparisons were made for up to 667 bp of DNA cloned from 14 kinds of HawaiianDrosophila and five other dipteran species. These sequences include parts of the genes for NADH dehydrogenase (subunits 1, 2, and 5) and rRNA (from the large ribosomal subunit). Because the times of divergence among these species are known approximately, the sequence comparisons give insight into the evolutionary dynamics of this molecule. Transitions account for nearly all of the differences between sequences that have diverged by less than 2%; for these sequences the mean rate of divergence appears to be about 2%/Myr. In comparisons involving greater divergence times and greater sequence divergence, relatively more of the sequence differences are due to transversions. Specifically, the fraction of these differences that are counted as transversions rises from an initial value of less than 0.1 to a plateau value of nearly 0.6. The time required to reach half of the plateau value, about 10 Myr, is similar to that for mammalian mtDNA. The mtDNAs of flies and mammals are also alike in the shape of the curve relating the percentage of positions at which there are differences in protein-coding regions to the time of divergence. For both groups of animals, the curve has a steep initial slope ascribable to fast accumulation of synonymous substitutions and a shallow final slope resulting from the slow accumulation of substitutions causing amino acid replacements. However, the percentage of all sites that can experience a high rate of substitution appears to be only about 8% for fly mtDNA compared to about 20% for mammalian mtDNA. The low percentage of hypervariable sites may be a consequence of a functional constraint associated with the low content of guanine and cytosine in fly mtDNA.     

Sequence length and variation in the mitochondrial control region of two freshwater gobiid fishes belonging to Rhinogobius (Teleostei: Gobioidei)

The control region (D-loop) of mitochondrial DNA (mtDNA) was amplified and sequenced for eight samples of the rhinogobies Rhinogobius maculafasciatus and R. giurinus from Taiwan and southern China. The control regions of both species are of 841–842 bp; the length of these sequences being the most compact among all known sequences in teleost fishes. Three conserved sequence blocks (CSB) were observed. The full D-loop and tRNA Phe gene sequences were determined and compared with other fishes. The interspecific sequence divergence between the two species is 11.3–11.7%; and the intraspecific variation in R. guirinus 0.8–1.8%. Results suggest that the control region of Rhinogobius is informative for phylogenetic reconstruction at both intraspecific and interspecific levels in this gobiid genus.     

大趾鼠耳蝠线粒体DNA控制区结构及变异

采用PCR产物直接测序法首次测定大趾鼠耳蝠(Myotis macrodactylus)10个个体的线粒体DNA(mtDNA)控制区全序列,并进行了结构和变异分析。结果表明,大趾鼠耳蝠的控制区结构与其他哺乳动物相似,可分为一个中央保守区(包括F、E、D、C、B元件)和两个外围结构域:延伸的终止结合序列区(包括ETAS1和ETAS2元件)和保守序列区(包括CSB1、CSB2和CSB3元件),其中最为保守的是中央保守区(核苷酸变异度为1.8%)。大趾鼠耳蝠控制区核苷酸全序列具有丰富的长度多态性(1778～2048bp),主要是由在碱基组成、重复数目和排列方式上异质的串联重复序列造成的。在ETAS内发现了TACAT及其反向互补序列ATGTA,支持滑移错配模式(slipped mispairing model)。本研究为该物种的进一步研究和保护提供基础遗传数据。     

Diversity and evolution of T-cell receptor variable region genes in mammals and birds

 The receptor of a T lymphocyte (TCR) recognizes nonself antigens in the company of major histocompatibility complex (MHC) molecules presented to it by the antigen-presenting cell. The variable region of TCR is encoded by either a concatenation of variable region (TCR-V), diversity region (TCR-D), and joining region (TCR-J) genes, or a concatenation of TCR-V and TCR-J genes. The TCR-V genes exist as a multigene family in vertebrate species. Here we study the evolutionary relationships of TCR-V genes from humans, sheep, cattle, rabbits, mice, and chicken. These six species can be classified into two groups according to the frequency of γδ T-cells in their peripheral T-cell populations. The "γδ low" group of species includes humans and mice, in which γδ T-cells constitute very limited portion of the T-cell population. The "γδ high" group includes sheep, cattle, rabbits, and chicken, in which γδ T-cells comprise up to 60% of the T-cell population. Here, we compiled TCR-V sequences from the six species and conducted a phylogenetic analysis. We identified various TCR-V gene subgroups based on the analysis. We found that humans and mice have representatives from nearly all of the subgroups identified, while other species have lost subgroups to different extent. Therefore, the γδ low species have a high degree of diversity of TCR-V genes, while γδ high species all have limited diversity of TCR-V genes. This pattern is similar to that found for immunoglobulin variable region (IGV) genes. Received: 20 May 1999 / Revised: 13 July 1999     

中华攀雀线粒体基因组全序列测定与分析

该研究使用长PCR扩增和引物步移法测定了中华攀雀(Remiz consobrinus)线粒体基因组全序列,在对序列进行拼接和注释的基础上,分析了其结构、序列组成及蛋白编码基因密码子使用情况等,并对22个tRNA和2个rRNA的二级结构以及控制区结构进行了预测及系统发育分析,为雀形目鸟类的系统发育研究提供了新信息。中华攀雀线粒体基因组全长16737bp,GenBank登录号KC463856,碱基A、T、C、G的含量分别为27.8%、21.5%、35.4%及15.3%,37个基因排列顺序与已报道的其他鸟类基本一致,包含13个蛋白编码基因、22个tRNA基因、2个rRNA基因及1个非编码的控制区(D-loop),有18对基因间共存在77bp的间隔,7对基因间共存在30bp的重叠。除ND3基因的起始密码子为ATT外,其余均为标准的ATG,11个蛋白编码基因的终止密码子为TAA、TAG、AGA或AGG,2个为不完全终止密码子T(COⅢ、ND4)。除tRNASer-AGNDHU臂缺失外,其余21个tRNA均可形成典型的三叶草结构,在出现的27处碱基错配中有19处为常见的G-U错配。SrRNA和LrRNA二级结构分别包含3个结构域47个茎环结构和6个结构域60个茎环结构,与所发表的其他鸟类rRNA二级结构大体一致。中华攀雀控制区发现了同样存在于其他鸟类控制区的保守框F-box、D-box、C-box、B-box、Bird similarity-box和CSB1-box。该研究支持将攀雀科作为独立的科,同时,支持莺总科与攀雀科的单系性。     

安徽汉族群体mtDNA高变区I序列多态性分析

以Anderson标准序列作为对照,用GeneDOC软件确定42个安徽汉族无关个体的mtDNA高变区I序列在线粒体基因组中的位置,通过序列比对软件clustalX分析安徽汉族群体mt DNA高变区I序列多态性,共检测到38种单倍型和57个变异位点.在mtDNA高变区I序列中14个bp的高变结构域中,安徽汉人16183位点变异率高达38%,在16187位点的变异率为4.8%.同时发现,安徽汉人与成都汉人在mtDNA高变区I 16183和16189位点的变异率接近,明显高于广东汉人.     

Molecular characterization and mitochondrial sequence variation in two sympatric species of Proechimys (Rodentia: Echimyidae) in French Guiana

Variations in mitochondrial DNA characters were used to characterize two morphologically similar and sympatric species of Neotropical terrestrial rodents of the genus Proechimys (Mammalia: Echimyidae). We sequenced both cytochrome b (1140 pb) and part of the control region (445 pb) from four individuals of P. cuvieri and five of P. cayennensis from French Guiana, which allowed us to depict intra- and inter-specific patterns of variation. The phylogenetic relationships between the nine sequences evidence the monophyly of each species, and illustrate that more polymorphism might exist within P. cuvieri than within P. cayennensis. By developing species-specific primers to amplify a fragment of the cytochrome b gene, we were able to identify 50 individuals of Proechimys spp. caught in two localities of French Guiana. In both sites of primary rainforests, we showed that the two species live in syntopy, and this observation emphasizes the need to document ecological differences which should exist in order to diminish inter-specific competition.     

Specific point mutations may not accumulate with aging in the mouse mitochondrial DNA control region

Increasing evidence suggests that mitochondrial function declines during aging in various tissues and in a wide range of organisms. This correlates with an age-dependent large accumulation of specific point mutations in the mtDNA control region that was reported recently in human fibroblast and skeletal muscle. However, evaluations of aging-related mtDNA mutations in other model animal systems. In this study, we analyzed mtDNA control regions of brain, skeletal muscle, heart, and other tissues from aged mice, in search of specific point mutations. A 948-bp fragment covering the entire mtDNA control region from various tissues of mice at the age of 25–26 months was sequenced. The sequence analysis was accomplished with a newly developed program Mutation Quantifier, which was able to accurately detect mutations with frequencies as low as 3%. Probably due to the relative shorter life-span, unlike what has been reported in human mtDNA, our results indicated there might be no significant accumulation of specific mutations in mouse mtDNA control region during aging.     

The complete mitochondrial genome of the leafminer Liriomyza sativae (Diptera: Agromyzidae): great difference in the A+T-rich region compared to Liriomyza trifolii

The complete mitochondrial genome sequence of Liriomyza sativae Blanchard (15,551 bp) was determined and analyzed in this study. The circular genome contained 37 genes including 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and an A + T-rich region. The initiation codons of COI and ND1 were ‘ATCA’ and ‘GTG’, respectively. ND2 gene used the truncated termination codon ‘T’. All the tRNA genes had the typical cloverleaf secondary structures except for tRNA^Ser(AGN) gene, which was found with the absence of a DHU arm. In addition, a tRNA-like secondary structure (tRNA^Met) was found in the A + T-rich region. The great difference was that the length of L. sativae A + T-rich region was 597 bp shorter than that of Liriomyza trifolii (Burgess). Meanwhile, some minor differences such as ‘TATA’ block were also observed in L. sativae in contrast to ‘TACA’ block in L. trifolii. There were also some essential structure elements such as ‘TATA’ block, ‘G(A)_nT’ block, poly-T stretch and stem-and-loop structure in the A + T-rich region of L. sativae mitochondrial genome.     

Evolutionary Relationships in the Sand-Dwelling Cichlid Lineage of Lake Tanganyika Suggest Multiple Colonization of Rocky Habitats and Convergent Origin of Biparental Mouthbrooding

The cichlid species flock of Lake Tanganyika is comprised of seven seeding lineages that evolved in step with changes of the lake environment. One seeding lineage diversified into at least six lineages within a short period of time. Our study focuses on the diversification of one of these lineages, the Ectodini, comprising highly specialized, sand- and rock-dwelling species. They display two distinct breeding styles: maternal and biparental mouthbrooding. By analyzing three mtDNA gene segments in 30 species representing all 13 described genera, we show that the Ectodini rapidly diversified into four clades at the onset of their radiation. The monotypic genus Grammatotria is likely to represent the most ancestral split, followed by the almost contemporary origin of three additional clades, the first comprising the benthic genus Callochromis, the second comprising the benthic genera Asprotilapia, Xenotilapia, Enantiopus, and Microdontochromis, and the third comprising the semi-pelagic genera Ophthalmotilapia, Cardiopharynx, Cyathopharynx, Ectodus, Aulonocranus, Lestradea, and Cunningtonia. Our study confirms the benthic and sand-dwelling life-style as ancestral. Rocky habitats were colonized independently in the Xenotilapia- and Ophthalmotilapia-clade. The Xenotilapia-clade comprises both maternal and biparental mouthbrooders. Their mode of breeding appears to be highly plastic: biparental mouthbrooding either evolved once in the common ancestor of the clade, to be reverted at least three times, or evolved at least five times independently from a maternally mouthbrooding ancestor. Furthermore, the genera Xenotilapia, Microdontochromis, Lestradea, and Ophthalmotilapia appeared paraphyletic in our analyses, suggesting the need of taxonomic revision.     

Sequence variation and evolution of the mitochondrial DNA control region in the musk shrew, Suncus murinus

The complete mitochondrial DNA (mtDNA) control region was cloned and sequenced in the musk shrew, Suncus murinus, Insectivora. The general aspect was similar to that found in other mammals. We have found in two locations of this region the presence of arrays of tandem repeats like those in other shrew species. One array was located in the left domain containing the termination-associated sequences (TAS) and the length of a copy was 77 bp. The other repeats were situated upstream from the recognition site for the end of H-strand replication in the right domain and were 20 bp long. The left halves of the control region containing the former repeats were sequenced and compared in several laboratory lines and wild animals from different localities, variations in copy number of repeated sequences were found both among individuals and within an individual. A comparative study of repeated sequences provides useful indication for the origin and evolution of tandem repeated sequences. Strand slippage and mispairing during replication of mtDNA with concerted manner is currently regarded as a dominant theory to account molecular mechanism for tandemly repeated sequences, and the pattern of sequence and length variation in our study supports this theory. Our results, however, suggest that the evolution of the repeated sequences containing the TAS in the musk shrew might go through the process of two steps; at the first step one complete repeated and several incomplete repeated sequences had reproduced in common ancestor of the shrew, and the second stage step-up of complete repeated sequences occurred with concerted evolution after differentiation into continental and insular groups.     

Ancient Divergence in Bathypelagic Lake Tanganyika Deepwater Cichlids: Mitochondrial Phylogeny of the Tribe Bathybatini

The cichlid species flock of Lake Tanganyika represents a polyphyletic assemblage of eight ancestral lineages, which colonized the emerging lake independently. Our study is focused on one of these lineages, the Bathybatini, a tribe of specialized piscivorous cichlids of the deep pelagic zone. By analyzing three mtDNA gene segments of all eight species of the tribe and two species of the closely related Trematocarini, we propose on the basis of a linearized tree analysis that the Bathybatini comprise two distinct lineages, the genera Hemibates and Bathybates, that seeded the primary lacustrine Tanganyika radiation independently. The genus Hemibates is likely to represent a distinct lineage that emerged simultaneously with the tribe Trematocarini and the genus Bathybates and should be therefore treated as a distinct tribe. Within the genus Bathybates, B. minor clearly represents the most ancestral split and is likely to have diverged from the remaining species in the course of the primary lacustrine Tanganyika radiation during which also the radiations of the Lamprologini and the H-lineage took place. The remaining large Bathybates species also diversified almost simultaneously and in step with the diversification of other Tanganyikan lineages—the Limnochromini and Cyprichromini—with B. graueri occupying the most ancestral branch, suggesting that these were induced by the same environmental changes. The lack of geographic color morphs suggests that competition and resource partitioning, rather than allopatric speciation, promoted speciation within the genus Bathybates.Reviewing Editor: Dr. Axel Meyer     

Structure and patterns of sequence variation in the mitochondrial DNA control region of the great cats

Mitochondrial DNA control region structure and variation were determined in the five species of the genus Panthera. Comparative analyses revealed two hypervariable segments, a central conserved region, and the occurrence of size and sequence heteroplasmy. As observed in the domestic cat, but not commonly seen in other animals, two repetitive sequence arrays (RS-2 with an 80-bp motif and RS-3 with a 6-10-bp motif) were identified. The 3' ends of RS-2 and RS-3 were highly conserved among species, suggesting that these motifs have different functional constraints. Control region sequences provided improved phylogenetic resolution grouping the sister taxa lion (Panthera leo) and leopard (Panthera pardus), with the jaguar (Panthera onca).     

Parallel origins of duplications and the formation of pseudogenes in mitochondrial DNA from parthenogenetic lizards (Heteronotia binoei; Gekkonidae)

Summary Analysis of mitochondrial DNAs (mtDNAs) from parthenogenetic lizards of theHeteronotia binoei complex with restriction enzymes revealed an 5-kb addition present in all 77 individuals. Cleavage site mapping suggested the presence of a direct tandem duplication spanning the 16S and 12S rRNA genes, the control region and most, if not all, of the gene for the subunit 1 of NADH dehydrogenase (ND1). The location of the duplication was confirmed by Southern hybridization. A restriction enzyme survey provided evidence for modifications to each copy of the duplicated sequence, including four large deletions. Each gene affected by a deletion was complemented by an intact version in the other copy of the sequence, although for one gene the functional copy was heteroplasmic for another deletion. Sequencing of a fragment from one copy of the duplication which encompassed the tRNA^leu(UUR) and parts of the 16S rRNA and ND1 genes, revealed mutations expected to disrupt function. Thus, evolution subsequent to the duplication event has resulted in mitochondrial pseudogenes. The presence of duplications in all of these parthenogens, but not among representatives of their maternal sexual ancestors, suggests that the duplications arose in the parthenogenetic form. This provides the second instance inH. binoei of mtDNA duplication associated with the transition from sexual to parthenogenetic reproduction. The increased incidence of duplications in parthenogenetic lizards may be caused by errors in mtDNA replication due to either polyploidy or hybridity of their nuclear genomes.