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.     

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.     

The mitochondrial genome of Priapulus caudatus Lamarck (Priapulida: Priapulidae)

We sequenced and annotated the complete mitochondrial (mt) genome of the priapulid Priapulus caudatus in order to provide a source of phylogenetic characters including an assessment of gene order arrangement. The genome was 14,919 bp in its entirety with few, short non-coding regions. A number of protein-coding and tRNA genes overlapped, making the genome relatively compact. The gene order was: cox1, cox2, trnK, trnD, atp8, atp6, cox3, trnG, nad3, trnA, trnR, trnN, rrnS, trnV, rrnL, trnL(yaa), trnL(nag), nad1, -trnS(nga), -cob, -nad6, trnP, -trnT, nad4L, nad4, trnH, nad5, trnF, -trnE, -trnS(nct), trnI, -trnQ, trnM, nad2, trnW, -trnC, -trnY; where '-' indicates genes transcribed on the opposite strand. The gene order, although unique amongst Metazoa, shared the greatest number of gene boundaries and the longest contiguous fragments with the chelicerate Limulus polyphemus. The mt genomes of these taxa differed only by a single inversion of 18 contiguous genes bounded by rrnS and trnS(nct). Other arthropods and nematodes shared fewer gene boundaries but considerably more than the most similar non-ecdysozoan.     

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 mitochondrial genome sequence of Whitmania pigra (Annelida, Hirudinea): the first representative from the class Hirudinea

The mitochondrial genome is a significant tool for investigating the evolutionary history of metazoan animals. The currently available mitochondrial genome data in GenBank is limited to understand the detail evolutionary relationship among the metazoan animals, especially in the phylum Annelida. Here we present the mitochondrial gene organization, gene order and codon usage of the leech Whitmania pigra (Annelida), which is the first representative from the class Hirudinea. It is a circular molecule of 14,426bp, and encodes 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. All 37 genes of W. pigra mitochondrial genome are transcribed from the same strand, which is identical to studied annelids, two echiurans, two sipunculans and many other lophotrochozoans. Five conserved gene clusters can be found in mitochondrial genomes of nine studied annelids, including (1) cox1-N-cox2; (2) cox3-Q-nad6-cob-W-atp6; (3) H-nad5-F-E-P-T-nad4L-nad4; (4) srRNA-V-lrRNA; and (5) nad3-S(1)-nad2. Compared with that of other studied annelids, translocations of transfer RNAs were found in the gene arrangement of W. pigra mitochondrial genome. Phylogenetic analysis strongly support that the species from Hirudinina and Oligochaeta form a monophyletic group Clitellata (BPM=100, BPP=100), which is consistent with previous research based on morphological and other molecular data. Both gene order data and amino acid sequences reveal that echiurans are derived annelids and sipunculans should be clustered with annelids and echiurans.     

Evidence for the existence of some dissociation in an otherwise strong linkage disequilibrium between mitochondrial and chloroplastic genomes in Cyclobalanopsis glauca

Variations in mitochondrial DNA in Cyclobalanopsis glauca (Thunb. ex Murray) Oerst. were studied in 140 trees from 32 populations collected from within the tree's natural range. By sequencing two mitochondrial DNA intron fragments (nad4/3-nad4/4r and nad7/2-nad7/3r), we revealed a total of 1788 bp and five polymorphic sites which allowed us to distinguish six mitotypes. The mitochondrial DNA markers provided replicated data to support population phylogeographical scenarios suggested previously using chloroplastic DNA markers. The gene genealogical tree of mitochondrial DNA was partially congruent with the chloroplastic DNA tree owing to the slower mutation rate and different mutational direction. Significant linkage disequilibrium existed between the two organellar genomes. Further paring analyses between fragments synthesized using different primers, accompanied by exclusion of polymorphic sites, showed that the random association could be attributed specifically to one of the polymorphic sites of the petG-trnP fragment of the chloroplastic genome, and the three polymorphic sites of the nad4/3-nad4/4r fragment of the mitochondrial genome. The former was inferred to derive from paternal leakage, and the latter from recurrent mutation. These polymorphic sites were also responsible for uncoupling of the combined gene tree of mitotype and chlorotype. In conclusion, specific fragments found in this study contribute to the incomplete congruence of the two organellar lineages that otherwise associate well phylogeographically.     

The complete mitochondrial genomes of Schistosoma haematobium and Schistosoma spindale and the evolutionary history of mitochondrial genome changes among parasitic flatworms

Complete mitochondrial genome sequences for the schistosomes Schistosoma haematobium and Schistosoma. spindale have been characterized. S. haematobium is the causative agent of urinary schistosomiasis in humans and S. spindale uses ruminants as its definitive host; both are transmitted by freshwater snail intermediate hosts. Results confirm a major gene order rearrangement among schistosomes in all traditional Schistosoma species groups other than Schistosoma japonicum; i.e., species groups S. mansoni, S. haematobium, and S. indicum. These data lend support to the 'out of Asia' (East and Southeast Asia) hypothesis for Schistosoma. The gene order change involves translocation of atp6-nad2-trnA and a rearrangement of nad3-nad1 relative to other parasitic flatworm mt genomes so far sequenced. Gene order and tRNA secondary structure changes (loss and acquisition of the DHU and/or TPsiC arms of trnC, trnF, and trnR) between mitochondrial genomes of these and other (digenean and cestode) flatworms were inferred by character mapping onto a phylogeny estimated from nuclear small subunit rRNA gene sequences of these same species, in order to find additional rare genomic changes suitable as synapomorphies. Denser and wider taxon sampling of mt genomes across the Platyhelminthes will validate these putative characters.     

The mitochondrial genome of the thermal dimorphic fungus Penicillium marneffei is more closely related to those of molds than yeasts

We report the complete sequence of the mitochondrial genome of Penicillium marneffei, the first complete mitochondrial DNA sequence of a thermal dimorphic fungus. This 35 kb mitochondrial genome contains the genes encoding ATP synthase subunits 6, 8, and 9 (atp6, atp8, and atp9), cytochrome oxidase subunits I, II, and III (cox1, cox2, and cox3), apocytochrome b (cob), reduced nicotinamide adenine dinucleotide ubiquinone oxireductase subunits (nad1, nad2, nad3, nad4, nad4L, nad5, and nad6), ribosomal protein of the small ribosomal subunit (rps), 28 tRNAs, and small and large ribosomal RNAs. Analysis of gene contents, gene orders, and gene sequences revealed that the mitochondrial genome of P. marneffei is more closely related to those of molds than yeasts.     

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.     

甜果螨全线粒体基因组测序与分析

【目的】测定和分析甜果螨Carpoglyphus lactis线粒体基因组全序列,并在线粒体基因组水平探讨其在真螨总目(Acariformes)中的系统发育地位,为真螨总目分类及果螨科线粒体基因组研究提供科学依据。【方法】挑取实验室饲养的甜果螨成螨,用传统的酚氯仿抽提法和试剂盒提取法提取甜果螨基因组DNA。然后采用节肢动物或螨类线粒体基因的通用引物PCR扩增出甜果螨线粒体基因cox1,cob,rrnS和nad4-nad5的部分序列;再设计种特异性引物进行Long-PCR扩增和步移法测序,测出甜果螨线粒体基因组全序列。应用SeqMan, SEQUIN 9.0和tRNAscan等生物信息学软件,对甜果螨线粒体基因组的基因结构等进行生物信息学分析。最后基于17种真螨总目螨类的蛋白质编码基因,采用最大似然法构建系统发育树。【结果】甜果螨线粒体全基因组总长为14 060 bp(GenBank登录号:MN073839),为典型的闭合双链DNA分子,共由37个基因组成,包括13个蛋白质编码基因(PCGs)、22个tRNA基因和2个rRNA基因;甜果螨线粒体基因组还包括1个大的非编码区(large non-coding region, LNR)。系统发育分析结果显示,甜果螨Carpoglyphus lactis属于无气门亚目粉螨总科(Acaroidae),与椭圆食粉螨Aleuroglyphus ovatus构成一支。粉螨总科(Acaroidae)和薄口螨总科(Histiostomatoidae)聚成一簇,与痒螨股(Psoroptidia)构成姐妹群。【结论】本研究首次获得并分析了甜果螨线粒体基因组全序列。甜果螨与椭圆食粉螨的亲缘关系较近。     

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).     

Molecular phylogeny of euthyneura (mollusca: gastropoda)

A new phylogenetic hypothesis for Euthyneura is proposed based on the analysis of primary sequence data (mitochondrial cox1, trnV, rrnL, trnL(cun), trnA, trnP, nad6, and nad5 genes) and the phylogenetic utility of two rare genomic changes (the relative position of the mitochondrial trnP gene, and an insertion/deletion event in a conserved region of the mitochondrial Cox1 protein) is addressed. Both sources of phylogenetic information clearly rejected the monophyly of pulmonates, a group of gastropods well supported so far by morphological evidence. The marine basommatophoran pulmonate Siphonaria was placed within opisthobranchs and shared with them the insertion of a Glycine in the Cox 1 protein. The marine systellommatophoran pulmonate Onchidella was recovered at the base of the opisthobranch + Siphonaria clade. Opisthobranchs, Siphonaria, and Onchidella shared the relative position of the mitochondrial trnP gene between the mitochondrial trnA and nad6 genes. The land snails and slugs (stylommatophoran pulmonates) were recovered as an early split in the phylogeny of advanced gastropods. The monophyly of the Euthyneura (Opisthobranchia + Pulmonata) was rejected by the inclusion of the heterostrophan Pyramidella.     

The Mitochondrial Genome of Chlorogonium elongatum Inferred from the Complete Sequence

The 22,704-bp circular mitochondrial DNA (mtDNA) of the chlamydomonad alga Chlorogonium elongatum was completely cloned and sequenced. The genome encodes seven proteins of the respiratory electron transport chain, subunit 1 of the cytochrome oxidase complex (cox1), apocytochrome b (cob), five subunits of the NADH dehydrogenase complex (nad1, nad2, nad4, nad5, and nad6), a set of three tRNAs (Q, W, M), and the large (LSU)- and small (SSU)-subunit ribosomal RNAs. Six group-I introns were found, two each in the cox1, cob, and nad5 genes. In each intron an open reading frame (ORF) related to maturases or endonucleases was identified. Both the LSU and the SSU rRNA genes are split into fragments intermingled with each other and with other genes. Although the average A + T content is 62.2%, GC-rich clusters were detected in intergenic regions, in variable domains of the rRNA genes, and in introns and intron-encoded ORFs. A comparison of the genome maps reveals that C. elongatum and Chlamydomonas eugametos mtDNAs are more closely related to one another than either is to Chlamydomonas reinhardtii mtDNA. Received: 3 November 1997 / Accepted: 12 January 1998     

Variation in sequence and RNA editing within core domains of mitochondrial group II introns among plants

The 3' regions of several group II introns within the mitochondrial genes nad1 and nad7 show unexpected sequence divergence among flowering plants, and the core domains 5 and 6 are predicted to have weaker helical structure than those in self-splicing group II introns. To assess whether RNA editing improves helical stability by the conversion of A-C mispairs to A-U pairs, we sequenced RT-PCR amplification products derived from excised intron RNAs or partially spliced precursors. Only in some cases was editing observed to strengthen the predicted helices. Moreover, the editing status within nad1 intron 1 and nad7 intron 4 was seen to differ among plant species, so that homologous intron sequences shared lower similarity at the RNA level than at the DNA level. Plant-specific variation was also seen in the length of the linker joining domains 5 and 6 of nad7 intron 3; it ranged from 4 nt in wheat to 11 nt in soybean, in contrast to the 2-4 nt length typical of classical group II introns. However, this intron is excised as a lariat structure with a domain 6 branchpoint adenosine. Our observations suggest that the core structures and sequences of these plant mitochondrial introns are subject to less stringent evolutionary constraints than conventional group II introns.     

萧氏松茎象线粒体基因组全序列测定与分析

象甲是鞘翅目中物种最丰富的类群, 目前关于其线粒体基因组全序列的研究还未见报道。本研究利用长距PCR和引物步移法对萧氏松茎象Hylobitelus xiaoi Zhang线粒体基因组全序列进行了测定。结果显示： 萧氏松茎象线粒体基因组序列全长16 123 bp（GenBank登录号为JX847496）, 共编码37个基因和1个非编码的控制区, 基因次序与典型的六足动物线粒体基因排列一致, 未发现基因重排现象。在基因组中两个值得注意的发现分别是： 1）N链上存在1个额外的trnV-like序列, 反密码子为GAC, 长度为69 bp, 其中65 bp与J链上的trnD重叠; 2）trnS^UCN和nad1之间存在1个长度为232 bp的基因间隔区。全部13个蛋白质编码基因的起始密码子均为ATN, 9个蛋白质编码基因的终止密码子为TAA, 其余4个蛋白质编码基因中, nad1和cox2的终止密码子为TAG, nad4和nad5则以不完整的终止密码子T作为终止信号。除trnS^AGN外, 其余的tRNAs均可形成典型的三叶草结构。而trnS^AGN的反密码子由TCT替代GCT, 反密码子臂延长形成9 bp（中间含1个碱基突起）, TΨC臂由正常的5 bp变为6 bp, DHU臂缩短仅1 bp, 各个臂之间没有连接碱基。线粒体控制区中包括10处长度不少于5 bp的poly-T（最长poly-T长度为14 bp）和2处微卫星样重复序列 (TA)₆和(TA)₉。本研究结果为探讨象甲总科在鞘翅目中的系统学地位及其与其他总科间的系统发生关系等问题提供了重要的分子生物学数据。     

Characterization and Phylogenetic Analysis of the Mitochondrial Genome of Shiraia bambusicola Reveals Special Features in the Order of Pleosporales

Shiraia bambusicola P. Henn. is a pathogenic fungus of bamboo, and its fruiting bodies are regarded as folk medicine. We determined and analyzed its complete mitochondrial DNA sequence (circular DNA molecule of 39,030 bp, G + C content of 25.19%). It contains the typical genes encoding proteins involved in electron transport and coupled oxidative phosphorylation (nad1-6 and nad4L, cob and cox1-3), one ATP synthase subunit (atp6), 4 hypothetical proteins, and two genes for large and small rRNAs (rnl and rns). There is a set of 32 tRNA genes comprising all 20 amino acids, and these genes are evenly distributed on the two strands. Phylogenetic analyses based on concatenated mitochondrial proteins indicated that S. bambusicola clustered with members of the order Pleosporales, which is in agreement with previous results. The gene arrangements of Dothideomycetes species contained three regions of gene orders partitioned in their mitochondrial genomes, including block 1 (nad6-atp6), block 2 (nad1-cox3) and block 3 (genes around rns). S. bambusicola displayed unique special features that differed from the other Pleosporales species, especially in the coding regions around rns (trnR-trnY). Moreover, a comparison of gene orders in mitochondrial genomes from Pezizomycotina revealed that although all encoded regions are located on the same strand in most Pezizomycotina mtDNAs, genes from Dothideomycetes species had different orientations, as well as diverse positions and colocalization of genes (such as cox3, cox1-cox2 and nad2–nad3); these distinctions were regarded as class-specific features. Interestingly, two incomplete copies of the atp6 gene were found on different strands of the mitogenomic DNA, a finding that has not been observed in the other analyzed fungal species. In our study, mitochondrial genomes from Dothideomycetes species were comprehensively analyzed for the first time, including many species that have not appeared in previous reports.