三斑海马线粒体基因组核苷酸全序列结构与分析

通过PCR扩增并测序获得了三斑海马(Hippocampus trimaculatus)线粒体DNA(mt DNA)全序列。三斑海马线粒体基因组全序列长度为16 534 bp(Gen Bank登录号为KJ956892),编码37个基因,包括13个蛋白编码基因、22个t RNA基因和2个r RNA基因。非编码区域包括1个控制区(D-loop)及一个轻链复制起始区域。大部分基因由H-链编码,包括14个t RNA基因、2个r RNA基因、12个蛋白编码基因;只有ND6和8个t RNA基因是在L-链编码。预测的22个t RNA基因的二级结构均为典型的三叶草状。基因间隔一般1~14 bp不等。此外,还存在7处碱基重叠,其中,4处是鱼类和脊椎动物典型的基因重叠位点。总的碱基含量分别为,A 32.7%,C 23.4%,G 14.6%,T 29.3%,A+T含量为62.0%。其线粒体基因组序列的结构与脊椎动物的典型结构近似。邻接法和贝叶斯法构建的三斑海马系统进化树的拓扑结构相似,这与现有的三斑海马的系统演化地位一致。本研究为海马的进化研究以及保护工作提供了基础数据。     

Nucleotide sequence of nine protein-coding genes and 22 tRNAs in the mitochondrial DNA of the sea starPisaster ochraceus

Summary We have cloned and sequenced over 9 kb of the mitochondrial genome from the sea starPisaster ochraceus. Within a continuous 8.0-kb fragment are located the genes for NADH dehydrogenase subunits 1, 2, 3, and 4L (ND1, ND2, ND3, and ND4L), cytochrome oxidase subunits I, II, and III (COI, COII, and COIII), and adenosine triphosphatase subunits 6 and 8 (ATPase 6 and ATPase 8). This large fragment also contains a cluster of 13 tRNA genes between ND1 and COI as well as the genes for isoleucine tRNA between ND1 and ND2, arginine tRNA between COI and ND4L, lysine tRNA between COII and ATPase 8, and the serine (UCN) tRNA between COIII and ND3. The genes for the other five tRNAs lie outside this fragment. The gene for phenylalanine tRNA is located between cytochrome b and the 12S ribosomal genes. The genes for tRNA^glu and tRNA^thr are 3 to the 12S ribosomal gene. The tRNAs for histidine and serine (AGN) are adjacent to each other and lie between ND4 and ND5. These data confirm the novel gene order in mitochondrial DNA (mtDNA) of sea stars and delineate additional distinctions between the sea star and other mtDNA molecules.     

赤麂线粒体全基因组的序列和结构

提取赤麂细胞株总DNA,参照我们实验室已测定的同属动物小麂线粒体全基因组序列设计引物,PCR扩增、测序、拼接,获得赤麂线粒体全基因组序列并进行生物信息学分析。赤麂线粒体全基因组序列全长16354bp。定位了22个tRNA基因、2个rRNA基因、13个蛋白编码基因和1个D-loop区。赤麂与小麂及其它哺乳动物线粒体的基因组结构相同,它们的序列同源性都较高。     

Genetic characterization of the mitochondrial DNA from Lepeophtheirus salmonis (Crustacea; Copepoda). A new gene organization revealed

The mitochondrial DNA (mtDNA) from the salmon louse, Lepeophtheirus salmonis, is 15445 bp. It includes the genes coding for cytochrome B (Cyt B), ATPase subunit 6 and 8 (A6 and A8), NADH dehydrogenase subunits 1-6 and 4L (ND1, ND2, ND3, ND4, ND4L, ND5 and ND6), cytochrome c oxidase subunits I-III (COI, COII and COIII), two rRNA genes (12S rRNA and 16S rRNA) and 22 tRNAs. Two copies of tRNA-Lys are present in the mtDNA of L. salmonis, while tRNA-Cys was not identified. Both DNA strands contain coding regions in the salmon louse, in contrast to the other copepod characterized Tigriopus japonicus, but only a few genes overlap. In vertebrates, ND4 and ND4L are transcribed as one bicistronic mRNA, and are therefore localized together. The same organization is also found in crustaceans, with the exceptions of T. japonicus, Neocalanus cristatus and L. salmonis that deviate from this pattern. Another exception of the L. salmonis mtDNA is that A6 and A8 do not overlap, but are separated by several genes. The protein-coding genes have a bias towards AT-rich codons. The mitochondrial gene order in L. salmonis differs significantly from the copepods T. japonicus, Eucalanus bungii, N. cristatus and the other 13 crustaceans previously characterized. Furthermore, the mitochondrial rRNA genes are encoded on opposite strands in L. salmonis. This has not been found in any other arthropods, but has been reported in two starfish species. In a phylogenetic analysis, using an alignment of mitochondrial protein sequences, L. salmonis groups together with T. japonicus, being distant relatives to the other crustaceans.     

Differences in mtDNA whole sequence between Tibetan and Han populations suggesting adaptive selection to high altitude

We performed a mitochondrial whole-genome comparison study in 40 Tibetan and 50 Han Chinese. All subjects could be classified into 13 haplogroups pertained to the Macrohaplogroup M and N that pitched different quadrants by principal component analysis. We observed a difference in the M9 haplogroup and identified 18 significant variants by comparing whole sequences between Tibetan and Han populations. Variants in ND2, COX2, tRNA alanine and 12S rRNA were predicted to confer increased protein stability in Tibetans. We compared the base substitutions of nonsynonymous (NS) versus synonymous (S) of 13 protein-encoding genes and found the NS/S values of the ATP6, ATP8, and Cyt b genes were larger (>1) in Tibetans than that in Han population. Our findings provide clues for the existence of adaptive selection for the ATP6, ATP8, Cyt b, ND2, COX2, tRNA alanine and 12S rRNA genes in Tibetans which likely contributed to adaptation to their specific geographic environment, such as high altitude.     

黑尾地鸦线粒体基因组序列测定与分析

采用长PCR扩增的线粒体DNA和引物步移法, 测定并注释了中国特有鸟类-黑尾地鸦(Podoces hendersoni)的线粒体基因组全序列。黑尾地鸦的mtDNA序列全长16 867 bp, GenBank登录号GU592504。基因含量和排列次序与原鸡的一致, 包含13个蛋白编码基因、22个tRNA、2个rRNA和1个控制区(D-loop)。除COI基因以GTG作为起始密码子外, 其余12个蛋白质编码基因均以典型ATG密码子起始。11个蛋白编码基因以完全终止密码子TAA、AGG或AGA终止, COIII和ND4基因终止密码子为不完整的T。tRNA^Ser(AGY)的DHU臂缺失, tRNA^Leu(CUN)的反密码子环由9个碱基构成, 而不是标准的7个碱基。其余的20个tRNA基因的二级结构均属典型的三叶草结构。预测了rRNA的二级结构, 其中, 12S rRNA二级结构包含4个结构域, 43个茎环结构; 16S rRNA的二级结构包含6个结构域, 55个茎环结构。此外, 在其他鸟类控制区中所发现的F-box、D-box、C-box、B-box、Bird similarity-box和CSB1-box也同样存在于黑尾地鸦中。     

黄毛纺蚋线粒体基因组全序列测定与分析

利用PCR步移法对黄毛纺蚋的线粒体基因组全序列进行了测定和分析。黄毛纺蚋线粒体基因组全长15904 bp(Gen Bank序列号KP793690),包括13个蛋白编码基因、22个tRNA基因、2个rRNA基因以及长度为939 bp的非编码区。A、T、C、G碱基含量分别为39.1%、35.8%、10.4%、14.7%。9个蛋白编码基因和14个tRNA基因在J链编码,其余4个蛋白编码基因和8个tRNA基因在N链编码,基因排列顺序与其它已知双翅目昆虫相同。13个蛋白编码基因中除COI以TTG作为起始密码外,其余蛋白质基因均以ATN作为起始密码子,终止密码子多数为典型的TAA、TAG,只有COI和ND4L以单独的T作为终止密码子。在所测得的22个tRNA基因中,除tRNASer(AGN)缺少DHU臂外,其余tRNA均能形成典型的三叶草结构。     

Phylogenetic analyses of the harvest mouse,Micromys minutus (Rodentia: Muridae) based on the complete mitogenome sequences

Harvest Mouse (Micromys minutus) has a very wide range of distribution in Asia and Europe. However, the phylogenetic relationship of M. minutus is still uncertain. In this study, we determined the complete mitochondrial (mt) genome sequences of M. minutus, and used the complete mitochondrial genome sequences constructed the phylogenetic tree of Muroidea. The size of the genome is 16,232 bp in length and has a base composition of 33.6% A, 29.1% T, 24.8% C, and 12.5% G. The mitogenome structure was similar to that of typical vertebrate and other rodents' mitochondrial genomes, includes 13 protein-coding genes, 2 rRNA genes (12S rRNA and 16S rRNA), 22 tRNA genes, and 1 control region. We suggested a new initiation codon for ND5 (NADH dehydrogenase subunit), which has been never reported in the mitochondrial genome of vertebrate. The ML and BI phylogenetic trees, which based on the combination of the 12 protein-coding genes, supported strongly that the genus Micromys was represent an early offshoot within the Muridae with high support values (BI = 1.00, ML = 100).     

Mitochondrial DNA in the sea urchin Arbacia lixula: evolutionary inferences from nucleotide sequence analysis.

From the stirodont Arbacia lixula we determined the sequence of 5,127 nucleotides of mitochondrial DNA (mtDNA) encompassing 18 tRNAs, two complete coding genes, parts of three other coding genes, and part of the 12S ribosomal RNA (rRNA). The sequence confirms that the organization of mtDNA is conserved within echinoids. Furthermore, it underlines the following peculiar features of sea urchin mtDNA: the clustering of tRNAs, the short noncoding regulatory sequence, and the separation by the ND1 and ND2 genes of the two rRNA genes. Comparison with the orthologous sequences from the camarodont species Paracentrotus lividus and Strongylocentrotus purpuratus revealed that (1) echinoids have an extra piece on the amino terminus of the ND5 gene that is probably the remnant of an old leucine tRNA gene; (2) third-position codon nucleotide usage has diverged between A. lixula and the camarodont species to a significant extent, implying different directional mutational pressures; and (3) the stirodont-camarodont divergence occurred twice as long ago as did the P. lividus-S. purpuratus divergence.     

Mitochondrial DNA Genomes Organization and Phylogenetic Relationships Analysis of Eight Anemonefishes (Pomacentridae: Amphiprioninae)

Anemonefishes (Pomacentridae Amphiprioninae) are a group of 30 valid coral reef fish species with their phylogenetic relationships still under debate. The eight available mitogenomes of anemonefishes were used to reconstruct the molecular phylogenetic tree; six were obtained from this study (Amphiprion clarkii, A. frenatus, A. percula, A. perideraion, A. polymnus and Premnas biaculeatus) and two from GenBank (A. bicinctus and A. ocellaris). The seven Amphiprion species represent all four subgenera and P. biaculeatus is the only species from Premnas. The eight mitogenomes of anemonefishes encoded 13 protein-coding genes, two rRNA genes, 22 tRNA genes and two main non-coding regions, with the gene arrangement and translation direction basically identical to other typical vertebrate mitogenomes. Among the 13 protein-coding genes, A. ocellaris (AP006017) and A. percula (KJ174497) had the same length in ND5 with 1,866 bp, which were three nucleotides less than the other six anemonefishes. Both structures of ND5, however, could translate to amino acid successfully. Only four mitogenomes had the tandem repeats in D-loop; the tandem repeats were located in downstream after Conserved Sequence Block rather than the upstream and repeated in a simply way. The phylogenetic utility was tested with Bayesian and Maximum Likelihood methods using all 13 protein-coding genes. The results strongly supported that the subfamily Amphiprioninae was monophyletic and P. biaculeatus should be assigned to the genus Amphiprion. Premnas biaculeatus with the percula complex were revealed to be the ancient anemonefish species. The tree forms of ND1, COIII, ND4, Cytb, Cytb+12S rRNA, Cytb+COI and Cytb+COI+12S rRNA were similar to that 13 protein-coding genes, therefore, we suggested that the suitable single mitochondrial gene for phylogenetic analysis of anemonefishes maybe Cytb. Additional mitogenomes of anemonefishes with a combination of nuclear markers will be useful to substantiate these conclusions in future studies.     

Three divergent mitochondrial genomes from California populations of the copepod Tigriopus californicus

Previous work on the harpacticoid copepod Tigriopus californicus has focused on the extensive population differentiation in three mtDNA protein coding genes (COXI, COXII, Cytb). In order to get a more complete understanding of mtDNA evolution in this species, we sequenced three complete mitochondrial genomes (one from each of three California populations) and compared them to two published mtDNA genomes from an Asian congener, Tigriopus japonicus. Several features of the mtDNA genome appear to be conserved within the genus: 1) the unique order of the protein coding genes, rRNA genes and most of the tRNA genes, 2) the genome is compact, varying between 14.3 and 14.6 kb, and 3) all genes are encoded on the same strand of the mtDNA. Within T. californicus, extremely high levels of nucleotide divergence (>20%) are observed across much of the mitochondrial genome. Inferred amino acid sequences of the proteins encoded in the mtDNAs also show high levels of divergence; at the extreme, the three ND3 variants in T. californicus showed >25% amino acid substitutions, compared with <3% amino acid divergence at the previously studied COXI locus. Unusual secondary structures make functional assignments of some tRNAs difficult. The only apparent tRNA(trp) in these genomes completely overlaps the 5' end of the 16S rRNA in all three T. californicus mtDNAs. Although not previously noted, this feature is also conserved in T. japonicus mtDNAs; whether this sequence is processed into a functional tRNA has not been determined. The putative control region contains a duplicated segment of different length (from 88 to 155 bp) in each of the T. californicus sequences. In each case, the duplicated segments are not tandem repeats; despite their different lengths, the distance between the start of the first and the start of the second repeat is conserved (520 bp). The functional significance, if any, of this repeat structure remains unknown.     

Nucleotide sequence and evolution of coding and noncoding regions of a quail mitochondrial genome

Summary Segments of the Japanese quail mito-chondrial genome encompassing many tRNA and protein genes, the small and part of the large rRNA genes, and the control region have been cloned and sequenced. Analysis of the relative position of these genes confirmed that the tRNA^Glu and ND6 genes in galliform mitochondrial DNA are located immediately adjacent to the control region of the molecule instead of between the cytochrome b and ND5 genes as in other vertebrates. Japanese quail and chicken display another distinctive characteristic, that is, they both lack an equivalent to the light-strand replication origin found between the tRNA^Cys and tRNA^Asn genes in all vertebrate mitochondrial genomes sequenced thus far. Comparison of the protein-encoding genes revealed that a great proportion of the substitutions are silent and involve mainly transitions. This bias toward transitions also occurs in the tRNA and rRNA genes but is not observed in the control region where transversions account for many of the substitutions. Sequence alignment indicated that the two avian control regions evolve mainly through base substitutions but are also characterized by the occurrence of a 57-bp deletion/addition event at their 5′ end. The overall sequence divergence between the two gallinaceous birds suggests that avian mitochondrial genomes evolve at a similar rate to other vertebrate mitochondrial DNAs.     

Maxicircle (mitochondrial) genome sequence (partial) of Leishmania major: gene content, arrangement and composition compared with Leishmania tarentolae

We report 8420 bp of DNA sequence data from the maxicircle (mitochondrial) genome of Leishmania major (MHOM/SU/73/5ASKH), a much larger portion of this genome than has been reported previously from any Leishmania species infecting humans. This region contains 10 partial and complete genes: 5 protein-encoding genes (COII, COIII, ND1, ND7 and Cyt b); two ribosomal RNA subunits (12S and 9S) and three unidentified open reading frames (MURF1, MURF4 (ATPase6) and MURF5), as in the lizard-infecting species L. tarentolae. The genes from L. major exhibit 85-87% identity with those of L. tarentolae at the nucleotide level and 71-94% identity at the amino acid level. Most differences between sequences from the two species are transversions. The gene order and arrangement within the maxicircle of L. major are similar to those in L. tarentolae, but base composition and codon usage differ between the species. Codons assigned for initiation for protein-coding genes available for comparison are similar in five genes in the two species. Pre-editing was identified in some of the protein-coding genes. Short intergenic non-coding regions are also present in L. major as they are in L. tarentolae. Intergenic regions between 9S rRNA and MURF5, MURF1 and ND1 genes are G+C rich and considered to be extensive RNA editing regions. The RNA editing process is likely to be conserved in similar pattern in L. major as in L. tarentolae.     

Intraspecific phylogeny of the New Zealand short-tailed bat Mystacina tuberculata inferred from multiple mitochondrial gene sequences

An intraspecific phylogeny was established for the New Zealand short-tailed bat Mystacina tuberculata using a 2,878-bp sequence alignment from multiple mitochondrial genes (control region, ND2, 12S ribosomal RNA [rRNA], 16S rRNA, and tRNA). The inferred phylogeny comprises six lineages, with estimated divergences extending back between 0.93 and 0.68 million years to the middle Pleistocene. The lineages do not correspond to the existing subspecific taxonomy. Although multiple lineages occur sympatrically in many populations, the lineages are geographically structured. This structure has persisted despite repeated cycles of range expansion and contraction in response to climatic oscillations and catastrophic volcanic eruptions. The distribution of lineages among populations in central North Island indicates that a hybrid zone was formed by simultaneous colonization from single-lineage source populations inhabiting remote forest refugia. The observed pattern is not typical of microbats, which because of their high mobility generally exhibit low levels of genetic differentiation and geographic structure over continental ranges. Although lineages of M. tuberculata occur sympatrically in many populations, genetic distances between them are sufficiently large to suggest that they may be considered evolutionary significant units or taxonomic subspecies.     

The Complete Mitochondrial Genome of the Stalk-Eyed Bug Chauliops fallax Scott,and the Monophyly of Malcidae (Hemiptera: Heteroptera)

Chauliops fallax Scott, 1874 (Hemiptera: Heteroptera: Malcidae: Chauliopinae) is one of the most destructive insect pests of soybean and rice fields in Asia. Here we sequenced the complete mitochondrial genome of this pest. This genome is 15,739 bp long, with an A+T content of 73.7%, containing 37 typical animal mitochondrial genes and a control region. All genes were arranged in the same order as most of other Heteroptera. A remarkable strand bias was found for all nine protein coding genes (PCGs) encoded by the majority strand were positive AT-skew and negative GC-skew, whereas the reverse were found in the remaining four PCGs encoded by the minority strand and two rRNA genes. The models of secondary structures for the two rRNA genes of sequenced true bugs and Lygaeoidea were predicted. 16S rRNA consisted of six domains (domain III is absent as in other known arthropod mitochondrial genomes) and 45 helices, while three domains and 27 helices for 12S rRNA. The control region consists of five subregions: a microsatellite-like region, a tandem repeats region and other three motifs. The unusual intergenic spacer between tRNA-H and ND4 only found in the species of Lygaeoidea, not in other heteropteran species, may be the synapomorphy of this superfamily. Phylogenetic analyses were carried out based on all the 13 PCGs showed that Chauliopinae was the sister group of Malcinae and the monophyly of Lygaeoidea.     

Neoplastic transformation is associated with coordinate induction of nuclear and cytoplasmic oxidative phosphorylation genes

Neoplastic transformation was found to have a marked effect on the expression of nuclear DNA (nDNA)- and mitochondrial DNA (mtDNA)-encoded oxidative phosphorylation (OXPHOS) genes. Examining three pairs of human diploid fibroblasts and their SV 40-transformed counterparts revealed that mRNAs for the nuclear-encoded ATP synthase beta and the adenine nucleotide translocator (ANT) isoform 1 and 2 genes were markedly induced, whereas the mRNA for the ANT isoform 3 gene remained unchanged. The mRNA levels for the mtDNA-encoded 12 S rRNA, ND2, ATPase6+8, COIII, ND5+6, and Cytb genes were also increased, whereas the mtDNA number declined. Similar analysis of a cervical carcinoma (HeLa), fibrosarcoma (HT1080), and an Epstein-Barr virus (EBV)-transformed lymphoblastoid line (EBV-L) revealed that all three ANT isoforms were also expressed in these cells. Hence, changes in the expression of OXPHOS genes may be a common feature of transformed cells.     

Use of mitogenomic information in teleostean molecular phylogenetics: a tree-based exploration under the maximum-parsimony optimality criterion

We explored the phylogenetic utility and limits of the individual and concatenated mitochondrial genes for reconstructing the higher-level relationships of teleosts, using the complete (or nearly complete) mitochondrial DNA sequences of eight teleosts (including three newly determined sequences), whose relative phylogenetic positions were noncontroversial. Maximum-parsimony analyses of the nucleotide and amino acid sequences of 13 protein-coding genes from the above eight teleosts, plus two outgroups (bichir and shark), indicated that all of the individual protein-coding genes, with the exception of ND5, failed to recover the expected phylogeny, although unambiguously aligned sequences from 22 concatenated transfer RNA (tRNA) genes (stem regions only) recovered the expected phylogeny successfully with moderate statistical support. The phylogenetic performance of the 13 protein-coding genes in recovering the expected phylogeny was roughly classified into five groups, viz. very good (ND5, ND4, COIII, COI), good (COII, cyt b), medium (ND3, ND2), poor (ND1, ATPase 6), and very poor (ND4L, ND6, ATPase 8). Although the universality of this observation was unclear, analysis of successive concatenation of the 13 protein-coding genes in the same ranking order revealed that the combined data sets comprising nucleotide sequences from the several top-ranked protein-coding genes (no 3rd codon positions) plus the 22 concatenated tRNA genes (stem regions only) best recovered the expected phylogeny, with all internal branches being supported by bootstrap values >90%. We conclude that judicious choice of mitochondrial genes and appropriate data weighting, in conjunction with purposeful taxonomic sampling, are prerequisites for resolving higher-level relationships in teleosts under the maximum-parsimony optimality criterion.