Analysis of the evolutionary forces shaping mitochondrial genomes of a Neotropical malaria vector complex

Many vectors of human malaria belong to complexes of morphologically indistinguishable cryptic species. Here we report the analysis of the newly sequenced complete mitochondrial DNA molecules from six recognized or putative species of one such group, the Neotropical Anopheles albitarsis complex. The molecular evolution of these genomes had been driven by purifying selection, particularly strongly acting on the RNA genes. Directional mutation pressure associated with the strand-asynchronous asymmetric mtDNA replication mechanism may have shaped a pronounced DNA strand asymmetry in the nucleotide composition in these and other Anopheles species. The distribution of sequence polymorphism, coupled with the conflicting phylogenetic trees inferred from the mitochondrial DNA and from the published white gene fragment sequences, indicates that the evolution of the complex may have involved ancient mtDNA introgression. Six protein coding genes (nad5, nad4, cox3, atp6, cox1 and nad2) have high levels of sequence divergence and are likely informative for population genetics studies. Finally, the extent of the mitochondrial DNA variation within the complex supports the notion that the complex consists of a larger number of species than until recently believed.     

The complete DNA sequence of the mitochondrial genome of Physarum polycephalum

The complete sequence of the mitochondrial DNA (mtDNA) of the true slime mold Physarun polycephalum has been determined. The mtDNA is a circular 62,862-bp molecule with an A+T content of 74.1%. A search with the program BLAST X identified the protein-coding regions. The mitochondrial genome of P. polycephalum was predicted to contain genes coding for 12 known proteins [for three cytochrome c oxidase subunits, apocytochrome b, two F1Fo-ATPase subunits, five NADH dehydrogenase (nad) subunits, and one ribosomal protein], two rRNA genes, and five tRNA genes. However, the predicted ORFs are not all in the same frame, because mitochondrial RNA in P. polycephalum undergoes RNA editing to produce functional RNAs. The nucleotide sequence of an nad7 cDNA showed that 51 nucleotides were inserted at 46 sites in the mRNA. No guide RNA-like sequences were observed in the mtDNA of P. polycephalum. Comparison with reported Physarum mtDNA sequences suggested that sites of RNA editing vary among strains. In the Physarum mtDNA, 20 ORFs of over 300 nucleotides were found and ORFs 14 19 are transcribed.     

The analysis of the complete mitochondrial genome of Lecanicillium muscarium (synonym Verticillium lecanii) suggests a minimum common gene organization in mtDNAs of Sordariomycetes: phylogenetic implications

The mitochondrial genome (mtDNA) of the entomopathogenic fungus Lecanicillium muscarium (synonym Verticillium lecanii) with a total size of 24,499-bp has been analyzed. So far, it is the smallest known mitochondrial genome among Pezizomycotina, with an extremely compact gene organization and only one group-I intron in its large ribosomal RNA (rnl) gene. It contains the 14 typical genes coding for proteins related to oxidative phosphorylation, the two rRNA genes, one intronic ORF coding for a possible ribosomal protein (rps), and a set of 25 tRNA genes which recognize codons for all amino acids, except alanine and cysteine. All genes are transcribed from the same DNA strand. Gene order comparison with all available complete fungal mtDNAs-representatives of all four Phyla are included-revealed some characteristic common features like uninterrupted gene pairs, overlapping genes, and extremely variable intergenic regions, that can all be exploited for the study of fungal mitochondrial genomes. Moreover, a minimum common mtDNA gene order could be detected, in two units, for all known Sordariomycetes namely nad1-nad4-atp8-atp6 and rns-cox3-rnl, which can be extended in Hypocreales, to nad4L-nad5-cob-cox1-nad1-nad4-atp8-atp6 and rns-cox3-rnl nad2-nad3, respectively. Phylogenetic analysis of all fungal mtDNA essential protein-coding genes as one unit, clearly demonstrated the superiority of small genome (mtDNA) over single gene comparisons.     

The complete mitochondrial genome of the black coral Chrysopathes formosa (Cnidaria:Anthozoa:Antipatharia) supports classification of antipatharians within the subclass Hexacorallia

Black corals comprise a globally distributed shallow- and deep-water taxon whose phylogenetic position within the Anthozoa has been debated. We sequenced the complete mitochondrial genome of the antipatharian Chrysopathes formosa to further evaluate its phylogenetic relationships. The circular mitochondrial genome (18,398 bp) consists of 13 energy pathway protein-coding genes and two ribosomal RNAs, but only two transfer RNA genes (trnM and trnW), as well as a group I intron within the nad5 gene that contains the only copies of nad1 and nad3. No novel genes were found in the antipatharian mitochondrial genome. Gene order and genome content are most similar to those of the sea anemone Metridium senile (subclass Hexacorallia), with differences being the relative location of three contiguous genes (cox2-nad4-nad6) and absence (from the antipatharian) of a group I intron within the cox1 gene. Phylogenetic analyses of multiple protein-coding genes support classifying the Antipatharia within the subclass Hexacorallia and not the subclass Ceriantipatharia; however, the sister-taxon relationships of black corals within Hexacorallia remain inconclusive.     

Extensive variation in nuclear mitochondrial DNA content between the genomes of Phytophthora sojae and Phytophthora ramorum

Fragments of mitochondrial DNA (mtDNA) transferred to the nuclear genome are called nuclear mitochondrial DNAs (NUMTs). We report here a comparison of NUMT content between genomes from two species of the same genus. Analysis of the genomes of Phytophthora sojae and P. ramorum revealed large differences in the NUMT content of the two genomes: 16.27 x 10(-3) and 2.28 x 10(-3)% of each genome, respectively. Substantial differences also exist between the two species in the sizes of the NUMTs found in each genome, with ranges of 20 to 405 bp for P. sojae and 19 to 137 bp for P. ramorum. Furthermore, in P. sojae, fragments from the mitochondrial genes rns, rnl, coxl, and nad (various subunits) are found most frequently, whereas P. ramorum NUMTs most often originate from the cox3, rpsl4, nad4, and nad5 genes. The large differences in the presumptive mtDNA insertions suggest that the insertions occurred subsequent to the divergence of the two species, and this is supported by sequence comparisons among the NUMTs and the mtDNA sequences of the two species. P. sojae mtDNA sequences inserted in the nuclear genome appear to have been altered as a result of insertions, deletions, inversions, and translocations and provide insights into active mechanisms of sequence divergence in this plant pathogen. No clear examples were found of NUMTs forming functional nuclear genes or of NUMTs inserted into exons or introns of any nuclear gene.     

Mitochondrial genome of the colorless green alga Polytomella parva: two linear DNA molecules with homologous inverted repeat Termini

Most of the well-characterized mitochondrial genomes from diverse green algal lineages are circular mapping DNA molecules; however, Chlamydomonas reinhardtii has a linear 15.8 kb unit mitochondrial genome with 580 or 581 bp inverted repeat ends. In mitochondrial-enriched fractions prepared from Polytomella parva (=P. agilis), a colorless, naturally wall-less relative of C. reinhardtii, we have detected two linear mitochondrial DNA (mtDNA) components with sizes of 13.5 and 3.5 kb. Sequences spanning 97% and 86% of the 13.5- and 3.5-kb mtDNAs, respectively, reveal that these molecules contain long, at least 1.3 kb, homologous inverted repeat sequences at their termini. The 3.5-kb mtDNA has only one coding region (nad6), the functionality of which is supported by both the relative rate at which it has accumulated nonsynonymous nucleotide substitutions and its absence from the 13.5-kb mtDNA which encodes nine genes (i.e., large and small subunit rRNA [LSU and SSU rRNA] genes, one tRNA gene, and six protein-coding genes). On the basis of DNA sequence data, we propose that a variant start codon, GTG, is utilized by the P. parva 13.5-kb mtDNA-encoded gene, nad5. Using the relative rate test with Chlamydomonas moewusii (=C. eugametos) as the outgroup, we conclude that the nonsynonymous nucleotide substitution rate in the mitochondrial protein-coding genes of P. parva is on an average about 3.3 times that of the C. reinhardtii counterparts.     

Rapid evolution of the compact and unusual mitochondrial genome in the ctenophore, Pleurobrachia bachei

Ctenophores are one of the most basally branching lineages of metazoans with the largest mitochondrial organelles in the animal kingdom. We sequenced the mitochondrial (mtDNA) genome from the Pacific cidipid ctenophore, Pleurobrachia bachei. The circular mitochondrial genome is 11,016 nts, with only 12 genes, and one of the smallest metazoan mtDNA genomes recorded. The protein coding genes are intronless cox1-3, cob, nad1, 3, 4, 4L and 5. The nad2 and 6 genes are represented as short fragments whereas the atp6 gene was found in the nuclear genome. Only the large ribosomal RNA subunit and two tRNAs were present with possibly the small subunit unidentifiable due to extensive fragmentation. The observed unique features of this mitochondrial genome suggest that nuclear and mitochondrial genomes have evolved at very different rates. This reduced mtDNA genome sharply contrasts with the very large sizes of mtDNA found in other basal metazoans including Porifera (sponges), and Placozoa (Trichoplax).     

Complete sequence of the mitochondrial genome of Tetrahymena pyriformis and comparison with Paramecium aurelia mitochondrial DNA

We report the complete nucleotide sequence of the Tetrahymena pyriformis mitochondrial genome and a comparison of its gene content and organization with that of Paramecium aurelia mtDNA. T. pyriformis mtDNA is a linear molecule of 47,172 bp (78.7 % A+T) excluding telomeric sequences (identical tandem repeats of 31 bp at each end of the genome). In addition to genes encoding the previously described bipartite small and large subunit rRNAs, the T. pyriformis mitochondrial genome contains 21 protein-coding genes that are clearly homologous to genes of defined function in other mtDNAs, including one (yejR) that specifies a component of a cytochrome c biogenesis pathway. As well, T. pyriformis mtDNA contains 22 open reading frames of unknown function larger than 60 codons, potentially specifying proteins ranging in size from 74 to 1386 amino acid residues. A total of 13 of these open reading frames ("ciliate-specific") are found in P. aurelia mtDNA, whereas the remaining nine appear to be unique to T. pyriformis; however, of the latter, five are positionally equivalent and of similar size in the two ciliate mitochondrial genomes, suggesting they may also be homologous, even though this is not evident from sequence comparisons. Only eight tRNA genes encoding seven distinct tRNAs are found in T. pyriformis mtDNA, formally confirming a long-standing proposal that most T. pyriformis mitochondrial tRNAs are nucleus-encoded species imported from the cytosol. Atypical features of mitochondrial gene organization and expression in T. pyriformis mtDNA include split and rearranged large subunit rRNA genes, as well as a split nad1 gene (encoding subunit 1 of NADH dehydrogenase of respiratory complex I) whose two segments are located on and transcribed from opposite strands, as is also the case in P. aurelia. Gene content and arrangement are very similar in T. pyriformis and P. aurelia mtDNAs, the two differing by a limited number of duplication, inversion and rearrangement events. Phylogenetic analyses using concatenated sequences of several mtDNA-encoded proteins provide high bootstrap support for the monophyly of alveolates (ciliates, dinoflagellates and apicomplexans) and slime molds.     

The complete mitochondrial genome sequence of Oncicola luehei (Acanthocephala: Archiacanthocephala) and its phylogenetic position within Syndermata

In the present study, we determined the complete mitochondrial genome sequence of Oncicola luehei (14,281bp), the first archiacanthocephalan representative and the second complete sequence from the phylum Acanthocephala. The complete genome contains 36 genes including 12 protein coding genes, 22 transfer RNA (tRNA) genes and 2 ribosomal RNA genes (rrnL and rrnS) as reported for other syndermatan species. All genes are encoded on the same strand. The overall nucleotide composition of O. luehei mtDNA is 37.7% T, 29.6% G, 22.5% A, and 10.2% C. The overall A+T content (60.2%) is much lower, compared to other syndermatan species reported so far, due to the high frequency (18.3%) of valine encoded by GTN in its protein-coding genes. Results from phylogenetic analyses of amino acid sequences for 10 protein-coding genes from 41 representatives of major metazoan groups including O. luehei supported monophyly of the phylum Acanthocephala and of the clade Syndermata (Acanthocephala+Rotifera), and the paraphyly of the clade Eurotatoria (classes Bdelloidea+Monogononta from phylum Rotifera). Considering the position of the acanthocephalan species within Syndermata, it is inferred that obligatory parasitism characteristic of acanthocephalans was acquired after the common ancestor of acanthocephalans diverged from its sister group, Bdelloidea. Additional comparison of complete mtDNA sequences from unsampled acanthocephalan lineages, especially classes Polyacanthocephala and Eoacanthocephala, is required to test if mtDNA provides reliable information for the evolutionary relationships and pattern of life history diversification found in the syndermatan groups.     

Numerous gene rearrangements in the mitochondrial genome of the wallaby louse, Heterodoxus macropus (Phthiraptera)

The complete arrangement of genes in the mitochondrial (mt) genome is known for 12 species of insects, and part of the gene arrangement in the mt genome is known for over 300 other species of insects. The arrangement of genes in the mt genome is very conserved in insects studied, since all of the protein-coding and rRNA genes and most of the tRNA genes are arranged in the same way. We sequenced the entire mt genome of the wallaby louse, Heterodoxus macropus, which is 14,670 bp long and has the 37 genes typical of animals and some noncoding regions. The largest noncoding region is 73 bp long (93% A+T), and the second largest is 47 bp long (92% A+T). Both of these noncoding regions seem to be able to form stem-loop structures. The arrangement of genes in the mt genome of this louse is unlike that of any other animal studied. All tRNA genes have moved and/or inverted relative to the ancestral gene arrangement of insects, which is present in the fruit fly Drosophila yakuba. At least nine protein-coding genes (atp6, atp8, cox2, cob, nad1-nad3, nad5, and nad6) have moved; moreover, four of these genes (atp6, atp8, nad1, and nad3) have inverted. The large number of gene rearrangements in the mt genome of H. macropus is unprecedented for an arthropod.     

Complete sequence of the mitochondrial genome of Taenia saginata: comparison with T. solium and T. asiatica

The complete sequence of the Taenia saginata mitochondrial genome was determined, and its organization and structure were compared to other human-tropic Taenia tapeworms for which complete mitochondrial sequence data were available. The mitochondrial genome was 13,670 bp long, contained 12 protein-coding genes, two ribosomal RNAs (rRNAs, a small and a large subunit), and 22 transfer RNAs (tRNAs). It did not encode the atp8 gene. Overlapping regions were found between nad4L and nad4, nad1 and trnN, and cox1 and trnT. The ATG initiation codon was used for 10 protein-coding genes, and the GTG initiation codon was used for the remaining 2 genes (nad4 and atp6). The size of the protein-coding genes of the three human Taenia tapeworms did not vary, except for Taenia solium nad1 (891 aa) and nad4 (1212 aa) and Taenia asiatica cox2 (576 aa). The tRNA genes were 57-75 bp long, and the predicted secondary structures of 18 of these genes had typical clover-leaf shapes with paired dihydrouridine (DHU) arms. The genes in all human Taenia tapeworms for the two mitochondrial rRNA subunits rrnL and rrnS are separated by trnC. The putative T. saginata rrnL and rrnS are 972 and 732 bp long, respectively. The non-coding regions of the mt genome of T. saginata consisted of 2 regions: a short non-coding region (SNR, 66 nucleotides) and a long non-coding region (LNR, 159 nucleotides). The overall sequence difference in the full mitochondrial genome between T. saginata and T. asiatica was 4.6%, while T. solium differed by 11%. In conclusion, the complete sequence of the T. saginata mitochondrial genome will serve as a resource for comparative mitochondrial genomics and systematic studies of the parasitic cestodes.     

The complete mitochondrial genome of Flustrellidra hispida and the phylogenetic position of Bryozoa among the Metazoa

The complete mitochondrial genome of Flustrellidra hispida (Bryozoa, Ctenostomata, Flustrellidridae) was sequenced using a transposon-mediated approach. All but one of the 36 genes were identified (trnS2). The genome is 13,026 bp long, being one of the smallest metazoan mitochondrial genomes sequenced to date with a unique gene order when compared to other Metazoa. The genome has an overall AT richness of 59.4%. We found seven regions of overlaps between tRNAs and protein-coding genes ranging from 2 to 11 nt, and seven regions of overlap between tRNAs, ranging from 1 to 8 nt, resulting in a total number of 46 overlapping nucleotides. Genes nad4, cox2, atp8, and nad3 are terminated by the abbreviated stop codon T and cytb is suggested to terminate on (ACT)AA; we postulate that mRNA editing is required to remove AC for TAA to be functional in terminating translation. Phylogenetic analysis of nucleotide and amino acid data place Flustrellidra in the Lophotrochozoa. DNA for this study originated from two populations resulting in a contig consisting of multiple haplotypes. Twenty-seven SNP sites were detected, the majority occurring in cox1 and nad5. With cox1 already established as a marker in bryozoan studies, we advocate the further testing of nad5.     

The complete mitochondrial genome of the nudibranch Roboastra europaea (Mollusca: Gastropoda) supports the monophyly of opisthobranchs

The complete nucleotide sequence (14,472 bp) of the mitochondrial genome of the nudibranch Roboastra europaea (Gastropoda: Opisthobranchia) was determined. This highly compact mitochondrial genome is nearly identical in gene organization to that found in opisthobranchs and pulmonates (Euthyneura) but not to that in prosobranchs (a paraphyletic group including the most basal lineages of gastropods). The newly determined mitochondrial genome differs only in the relative position of the trnC gene when compared with the mitochondrial genome of Pupa strigosa, the only opisthobranch mitochondrial genome sequenced so far. Pupa and Roboastra represent the most basal and derived lineages of opisthobranchs, respectively, and their mitochondrial genomes are more similar in sequence when compared with those of pulmonates. All phylogenetic analyses (maximum parsimony, minimum evolution, maximum likelihood, and Bayesian) based on the deduced amino acid sequences of all mitochondrial protein-coding genes supported the monophyly of opisthobranchs. These results are in agreement with the classical view that recognizes Opisthobranchia as a natural group and contradict recent phylogenetic studies of the group based on shorter sequence data sets. The monophyly of opisthobranchs was further confirmed when a fragment of 2,500 nucleotides including the mitochondrial cox1, rrnL, nad6, and nad5 genes was analyzed in several species representing five different orders of opisthobranchs with all common methods of phylogenetic inference. Within opisthobranchs, the polyphyly of cephalaspideans and the monophyly of nudibranchs were recovered. The evolution of mitochondrial tRNA rearrangements was analyzed using the cox1+rrnL+nad6+nad5 gene phylogeny. The relative position of the trnP gene between the trnA and nad6 genes was found to be a synapomorphy of opisthobranchs that supports their monophyly.     

残锷线蛱蝶线粒体基因组全序列及其系统学意义(英文)

对残锷线蛱蝶(Parathymasulpitia)(鳞翅目:蛱蝶科)线粒体基因组全序列进行了测定。结果表明:残锷线蛱蝶线粒体基因组全序列全长为15268bp,除了在trnS1(AGN)和trnE基因之间有一段121bp长的基因间隔外,其基因的排列顺序及排列方向与大多数已测鳞翅目物种基本一致。在蛋白质编码基因中,除cox1以CGA作为其起始密码子之外,其余12个蛋白质编码基因都以标准的ATN作为起始密码子。此外,除nad4基因以单独的T为终止密码子,其余12个蛋白质编码基因都以TAA结尾。除trnS1(AGN)缺少DHU臂之外,22个tRNA基因都显示典型的三叶草形二级结构。除A+T富集区外的非编码序列中,线粒体基因组共含有11个基因间隔区。其中,最长的一个121bp的基因间隔区位于trnS1(AGN)和trnE之间,其A+T含量高达100%。另外,和其他鳞翅目物种一样,在其A+T富集区的3’端有一段长达18bp的poly-T结构。A+T富集区内部没有明显的小卫星样多拷贝重复序列,而含有一些微卫星样的重复结构。本研究基于13种蛋白编码基因序列的组合数据,用最大似然法和贝叶斯法对蛱蝶科几个主要亚科间共9个代表物种间的系统发生关系进行了分析。结果表明,本研究的结果与前人的分子系统学研究结论基本吻合(其中,线蛱蝶亚科和釉蛱蝶亚科互为姐妹群),而与形态学的研究结论不一致。     

The complete mitochondrial genome of the helmet catfish Cranoglanis bouderius (Siluriformes: Cranoglanididae) and the phylogeny of otophysan fishes

The complete sequence of the 16,539 nucleotide mitochondrial genome from the single species of the catfish family Cranoglanididae, the helmet catfish Cranoglanis bouderius, was determined using the long and accurate polymerase chain reaction (LA PCR) method. The nucleotide sequences of C. bouderius mitochondrial DNA have been compared with those of three other catfish species in the same order. The contents of the C. bouderius mitochondrial genome are 13 protein-coding genes, two ribosomal RNA and 22 transfer RNA genes, and a non-coding control region, the gene order of which is identical to that observed in most other vertebrates. Phylogenetic analyses for 13 otophysan fishes were performed using Bayesian method based on the concatenated mtDNA protein-coding gene sequence and the individual protein-coding gene sequence data set. The competing otophysan topologies were then tested by using the approximately unbiased test, the Kishino-Hasegawa test, and the Shimodaira-Hasegawa test. The results show that the grouping ((((Characiformes, Gymnotiformes), Siluriformes), Cypriniformes), outgroup) is the most likely but there is no significant difference between this one and the other alternative hypotheses. In addition, the phylogenetic placement of the family Cranoglanididae among siluriform families was also discussed.     

The complete mitochondrial genome of the grand jackknife clam, Solen grandis (Bivalvia: Solenidae): a novel gene order and unusual non-coding region

Molluscs in general, and bivalves in particular, exhibit an extraordinary degree of mitochondrial gene order variation when compared with other metazoans. The complete mitochondrial genome of Solen grandis (Bivalvia: Solenidae) was determined using long-PCR and genome walking techniques. The entire mitochondrial genome sequence of S. grandis is 16,784 bp in length, and contains 36 genes including 12 protein-coding genes (atp8 is absent), 2 ribosomal RNAs, and 22 tRNAs. All genes are encoded on the same strand. Compared with other species, it bears a novel gene order. Besides these, we find a peculiar non-coding region of 435 bp with a microsatellite-like (TA)₁₂ element, poly-structures and many hairpin structures. In contrast to the available heterodont mitochondrial genomes from GenBank, the complete mtDNA of S. grandis has the shortest cox3 gene, and the longest atp6, nad4, nad5 genes.     

Characterization of mitochondrial genome of Chinese wild mulberry silkworm, Bomyx mandarina (Lepidoptera: Bombycidae)

The complete mitochondrial genome of Chinese Bombyx mandarina(ChBm) was determined.The cir-cular genome is 15682 bp long, and contains a typical gene complement, order, and arrangement iden-tical to that of Bombyx mori(B.mori) and Japanese Bombyx mandarina(JaBm) except for two addi-tional tRNA-like structures:tRNASer(TGA)-like and tRNAIle(TAT)-like.All protein-coding sequences are initi-ated with a typical ATN codon except for the COI gene, which has a 4-bp TTAG putative initiator codon.Eleven of 13 protein-coding genes(PCGs) have a complete termination codon(all TAA), but the re-maining two genes terminate with incomplete codons.All tRNAs have the typical clover-leaf structures of mitochondrial tRNAs, with the exception of tRNASer(TGA)-like, with a four stem-and-loop structure.The length of the A T-rich region of ChBm is 484 bp, shorter than those of JaBm(747 bp) and B.mori(494―499 bp).Phylogenetic analysis among B.mori, ChBm, JaBm, and Antheraea pernyi(Anpe) showed that B.mori is more closely related to ChBm than JaBm.The earliest divergence time estimate for B.mori-ChBm and B.mori-JaBm is about 1.08±0.18―1.41±0.24 and 1.53±0.20―2.01±0.26 Mya, respec-tively.ChBm and JaBm diverged around 1.11±0.16―1.45±0.21 Mya.