The complete mitochondrial genome sequence of the pathogenic yeast Candida (Torulopsis) glabrata

We report here the complete sequence of the mitochondrial (mt) genome of the pathogenic yeast Candida glabrata. This 20 kb mt genome is the smallest among sequenced hemiascomycetous yeasts. Despite its compaction, the mt genome contains the genes encoding the apocytochrome b (COB), three subunits of ATP synthetase (ATP6, 8 and 9), three subunits of cytochrome oxidase (COX1, 2 and 3), the ribosomal protein VAR1, 23 tRNAs, small and large ribosomal RNAs and the RNA subunit of RNase P. Three group I introns each with an intronic open reading frame are present in the COX1 gene. This sequence is available under accession number AJ511533.     

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.     

Structure and organization of the mitochondrial genome of the unicellular red alga Cyanidioschyzon merolae deduced from the complete nucleotide sequence.

The complete nucleotide sequence of the mitochondrial genome of a very primitive unicellular red alga, Cyanidioschyzon merolae , has been determined. The mitochondrial genome of C.merolae contains 34 genes for proteins including unidentified open reading frames (ORFs) (three subunits of cytochrome c oxidase, apocytochrome b protein, three subunits of F1F0-ATPase, seven subunits of NADH ubiquinone oxidoreductase, three subunits of succinate dehydrogenase, four proteins implicated in c-type cytochrome biogenesis, 11 ribosomal subunits and two unidentified open reading frames), three genes for rRNAs and 25 genes for tRNAs. The G+C content of this mitochondrial genome is 27.2%. The genes are encoded on both strands. The genome size is comparatively small for a plant mitochondrial genome (32 211 bp). The mitochondrial genome resembles those of plants in its gene content because it contains several ribosomal protein genes and ORFs shared by other plant mitochondrial genomes. In contrast, it resembles those of animals in the genome organization, because it has very short intergenic regions and no introns. The gene set in this mitochondrial genome is a subset of that of Reclinomonas americana , an amoeboid protozoan. The results suggest that plant mitochondria originate from the same ancestor as other mitochondria and that most genes were lost from the mitochondrial genome at a fairly early stage of the evolution of the plants.     

The complete mitochondrial genome of yarrowia lipolytica

We here report the complete nucleotide sequence of the 47.9 kb mitochondrial (mt) genome from the obligate aerobic yeast Yarrowia lipolytica. It encodes, all on the same strand, seven subunits of NADH: ubiquinone oxidoreductase (ND1-6, ND4L), apocytochrome b (COB), three subunits of cytochrome oxidase (COX1, 2, 3), three subunits of ATP synthetase (ATP6, 8 and 9), small and large ribosomal RNAs and an incomplete set of tRNAs. The Y. lipolytica mt genome is very similar to the Hansenula wingei mt genome, as judged from blocks of conserved gene order and from sequence homology. The extra DNA in the Y. lipolytica mt genome consists of 17 group 1 introns and stretches of A+Trich sequence, interspersed with potentially transposable GC clusters. The usual mould mt genetic code is used. Interestingly, there is no tRNA able to read CGN (arginine) codons. CGN codons could not be found in exonic open reading frames, whereas they do occur in intronic open reading frames. However, several of the intronic open reading frames have accumulated mutations and must be regarded as pseudogenes. We propose that this may have been triggered by the presence of untranslatable CGN codons. This sequence is available under EMBL Accession No. AJ307410.     

Two circular chromosomes of unequal copy number make up the mitochondrial genome of the rotifer Brachionus plicatilis

The monogonont rotifer Brachionus plicatilis is an emerging model system for a diverse array of questions in limnological ecosystem dynamics, the evolution of sexual recombination, cryptic speciation, and the phylogeny of basal metazoans. We sequenced the complete mitochondrial genome of B. plicatilis sensu strictu NH1L and found that it is composed of 2 circular chromosomes, designated mtDNA-I (11,153 bp) and mtDNA-II (12,672 bp). Hybridization to DNA isolated from mitochondria demonstrated that mtDNA-I is present at 4 times the copy number of mtDNA-II. The only nucleotide similarity between the 2 chromosomes is a 4.9-kbp region of 99.5% identity including a transfer RNA (tRNA) gene and an extensive noncoding region that contains putative D-loop and control sequence. The mtDNA-I chromosome encodes 4 proteins (ATP6, COB, NAD1, and NAD2), 13 tRNAs, and the large and small subunit ribosomal RNAs; mtDNA-II encodes 8 proteins (COX1-3, NAD3-6, and NAD4L) and 9 tRNAs. Gene order is not conserved between B. plicatilis and its closest relative with a sequenced mitochondrial genome, the acanthocephalan Leptorhynchoides thecatus, or other sequenced mitochondrial genomes. Polymerase chain reaction assays and Southern hybridization to DNA from 18 strains of Brachionus suggest that the 2-chromosome structure has been stable for millions of years. The novel organization of the B. plicatilis mitochondrial genome into 2 nearly equal chromosomes of 4-fold different copy number may provide insight into the evolution of metazoan mitochondria and the phylogenetics of rotifers and other basal animal phyla.     

A nuclear gene encoding the beta subunit of the mitochondrial ATP synthase in Nicotiana plumbaginifolia.

The beta subunit of the mitochondrial ATP synthase in Nicotiana plumbaginifolia is encoded by two nuclear genes, atp2-1 and atp2-2, which are both expressed. The complete nucleotide sequence of atp2-1 has been determined. It contains eight introns ranging from 88 to 1453 bp. The last intron contains a putative insertion element (Inp), of 812 bp bordered by 35-bp inverted repeats which share an 11-bp homology with Agrobacterium tumefaciens T-DNA borders. Sequences homologous to Inp are present in multiple copies in the N. plumbaginifolia and the N. tabacum genome but not in more distant species. The atp2-1 encoded polypeptide is highly homologous to beta subunits from other ATP synthases but it contains an extension at the N terminus which is probably involved in mitochondrial targeting. A sequence homology between exon 4 of atp2-1 and exon 1 of the human ras genes suggests a common ancestral origin for these exons.     

大卫绢蛱蝶线粒体基因组全序列测定和分析

目前有关蝶类线粒体基因组全序列及其分子进化的研究报道还不多见。本文利用long PCR和引物步移法得到大卫绢蛱蝶Calinaga davidis的线粒体基因组全序列, 同时就其基因组成和结构特点作了初步分析。结果显示： 其基因组全长为15 267 bp (GenBank登录号为HQ658143), 包括13个蛋白质编码基因(ATP6, ATP8, COI-III, ND1-6, ND4L, Cytb)、22个tRNA基因、2个rRNA基因(16S和12S)以及非编码的控制区。与其他鳞翅目昆虫相一致, 其基因组未出现基因重排现象。基因组共包含11个基因间隔区,总长度为130 bp, 间隔长度1~46 bp, 最大间隔在tRNA^Gln与ND2基因之间; 基因间共存在13处重叠, 总长度为66 bp, 重叠碱基数1~35 bp, 最长的重叠区位于COII与tRNA^Lys基因。lrRNA和srRNA基因长度分别为1 337 bp和773 bp; 除tRNA^Ser(AGN)缺少二氢尿嘧啶臂(DHU stem), 在相应的位置上只形成一个简单环外, 其余的tRNA基因都能形成典型的三叶草结构。13个蛋白编码基因总长度为11 247 bp, 共有3 737个密码子, 它们的碱基组成和密码子的使用具有明显的偏倚性; 除COI外(起始密码子TTG), 其余的12个蛋白质编码基因都以标准的ATN作为起始密码子; COI基因终止密码子为不完全T, ND4基因终止密码子为不完全TA, 其余基因都以TAA为终止密码子。A+T丰富区全长为389 bp, A+T含量高达92.0%, 其中存在2段类似微卫星的重复序列(TA)₆和(AAT)₄。本文的研究结果为探讨绢蛱蝶亚科在蛱蝶科中的系统学地位及其与其他亚科间的系统发生关系等问题提供了重要的分子生物学数据。     

池蝶蚌线粒体基因组全序列分析

采用普通PCR扩增、SHOT-GUN测序、软件拼接首次获得了池蝶蚌（Hyriopsis schlegelii）线粒体基因组全序列。线粒体基因组全长为15939 bp,由13个蛋白质编码基因、22个tRNA基因、2个SrRNA基因和28个长度为1393 bp的非编码区组成;除ND3-ND5、ND4L、ATP6、ATP8、COX1-COX3、tRNA-D、tRNA-H之外,其他大多数基因在L链编码。池蝶蚌线粒体全基因组序列、蛋白编码基因、tRNA基因、rRNA基因及非编码区的A+T含量分别为60.36%、59.84%、61.7%、60.23%及62.5%,与其他淡水蚌类一致,均表现出A+T偏好性,淡水蚌类线粒体基因组长度的差异主要表现在非编码区长度的差异。池蝶蚌mtDNA的COX2-12SrRNA区域基因排列存在差异,是ND3、tRNAHis、tRNAAla、tRNASer1、tRNASer2、tRNAGlu、ND2、tRNAMet 8个基因发生重组造成。22个tRNA基因都具有典型的三叶草二级结构,tRNA-E与 tRNA-W间的非编码区含有一个ORF区,而控制区并未发现。从GenBank上下载的14种双壳纲贝类的mtDNA序列构建的系统进化树,显示池蝶蚌与三角帆蚌亲缘关系最近。研究结果为淡水珍珠蚌线粒体基因重排及进化特征提供理论依据。       

Mitogenome phylogenetics: the impact of using single regions and partitioning schemes on topology, substitution rate and divergence time estimation

The availability of mitochondrial genome sequences is growing as a result of recent technological advances in molecular biology. In phylogenetic analyses, the complete mitogenome is increasingly becoming the marker of choice, usually providing better phylogenetic resolution and precision relative to traditional markers such as cytochrome b (CYTB) and the control region (CR). In some cases, the differences in phylogenetic estimates between mitogenomic and single-gene markers have yielded incongruent conclusions. By comparing phylogenetic estimates made from different genes, we identified the most informative mitochondrial regions and evaluated the minimum amount of data necessary to reproduce the same results as the mitogenome. We compared results among individual genes and the mitogenome for recently published complete mitogenome datasets of selected delphinids (Delphinidae) and killer whales (genus Orcinus). Using Bayesian phylogenetic methods, we investigated differences in estimation of topologies, divergence dates, and clock-like behavior among genes for both datasets. Although the most informative regions were not the same for each taxonomic group (COX1, CYTB, ND3 and ATP6 for Orcinus, and ND1, COX1 and ND4 for Delphinidae), in both cases they were equivalent to less than a quarter of the complete mitogenome. This suggests that gene information content can vary among groups, but can be adequately represented by a portion of the complete sequence. Although our results indicate that complete mitogenomes provide the highest phylogenetic resolution and most precise date estimates, a minimum amount of data can be selected using our approach when the complete sequence is unavailable. Studies based on single genes can benefit from the addition of a few more mitochondrial markers, producing topologies and date estimates similar to those obtained using the entire mitogenome.     

The complete mitochondrial genome sequence and gene organization of the mud crab (Scylla paramamosain) with phylogenetic consideration

The complete mitochondrial genome is of great importance for better understanding the genome-level characteristics and phylogenetic relationships among related species. In the present study, we determined the complete mitochondrial genome DNA sequence of the mud crab (Scylla paramamosain) by 454 deep sequencing and Sanger sequencing approaches. The complete genome DNA was 15,824 bp in length and contained a typical set of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and a putative control region (CR). Of 37 genes, twenty-three were encoded by the heavy strand (H-strand), while the other ones were encoded by light strand (L-strand). The gene order in the mitochondrial genome was largely identical to those obtained in most arthropods, although the relative position of gene tRNA^His differed from other arthropods. Among 13 protein-coding genes, three (ATPase subunit 6 (ATP6), NADH dehydrogenase subunits 1 (ND1) and ND3) started with a rare start codon ATT, whereas, one gene cytochrome c oxidase subunit I (COI) ended with the incomplete stop codon TA. All 22 tRNAs could fold into a typical clover-leaf secondary structure, with the gene sizes ranging from 63 to 73 bp. The phylogenetic analysis based on 12 concatenated protein-coding genes showed that the molecular genetic relationship of 19 species of 11 genera was identical to the traditional taxonomy.     

绿眼赛茧蜂线粒体基因组全序列测定和分析

【目的】测定绿眼赛茧蜂Zele chlorophthalmus线粒体基因组全序列,分析其基因组结构及茧蜂科(Braconidae)部分类群的系统发育关系。【方法】利用Illumina MiSeq二代测序技术对绿眼赛茧蜂的线粒体基因组进行测序,对基因组序列进行拼装、注释,分析其结构特点和碱基组成;基于22种茧蜂科昆虫的COX1蛋白编码基因序列,应用最大似然法(ML)和邻接法(NJ)构建系统发育树,分析绿眼赛茧蜂与茧蜂科其他昆虫的系统发育关系。【结果】绿眼赛茧蜂线粒体基因组全长16 661 bp(GenBank登录号: MG822749),包含13个蛋白质编码基因、22个tRNA基因和2个rRNA基因,共37个基因,以及1个控制区。线粒体基因组有明显的核苷酸组成的偏倚,AT偏正,GC偏负,其A+T含量为82.83%。基因排列顺序与推测的昆虫祖先的序列不完全一致,tRNA基因7处发生重排。13个蛋白质编码基因均以ATN为起始密码子,以TAA为终止密码子。在22个tRNA基因二级结构中,除tRNA^His(H)缺失TΨC环和tRNA^Cys(C)仅剩二氢尿嘧啶(DHU)臂和反密码子臂外,其余tRNA基因均能形成典型的三叶草结构。基于COX1蛋白编码序列的系统发育分析结果显示,与绿眼赛茧蜂亲缘关系最近的是同属于赛茧蜂属的雪跗赛茧蜂Z. niveitarsis。【结论】本研究首次获得绿眼赛茧蜂线粒体基因组全序列。结果表明绿眼赛茧蜂隶属于优茧蜂亚科(Euphorinae)赛茧蜂属,并支持赛茧蜂属的单系性。     

A Novel Mitochondrial Genome Organization for the Blue Mussel, Mytilus Edulis

The sequence of 13.9 kilobases (kb) of the 17.1-kb mitochondrial genome of Mytilus edulis has been determined, and the arrangement of all genes has been deduced. Mytilus mitochondrial DNA (mtDNA) contains 37 genes, all of which are transcribed from the same DNA strand. The gene content of Mytilus is typically metazoan in that it includes genes for large and small ribosomal RNAs, for a complete set of transfer RNAs and for 12 proteins. The protein genes encode the cytochrome b apoenzyme, cytochrome c oxidase (CO) subunits I-III, NADH dehydrogenase (ND) subunits 1-6 and 4L, and ATP synthetase (ATPase) subunit 6. No gene for ATPase subunit 8 could be found. The reading frames for the ND1, COI, and COIII genes contain long extensions relative to those genes in other metazoan mtDNAs. There are 23 tRNA genes, one more than previously found in any metazoan mtDNA. The additional tRNA appears to specify methionine, making Mytilus mtDNA unique in having two tRNA(Met) genes. Five lengthy unassigned intergenic sequences are present, four of which vary in length from 79 to 119 nucleotides and the largest of which is 1.2 kb. The base compositions of these are unremarkable and do not differ significantly from that of the remainder of the mtDNA. The arrangement of genes in Mytilus mtDNA is remarkably unlike that found in any other known metazoan mtDNA.     

Molecular evolution of cytochrome c oxidase: rate variation among subunit VIa isoforms

Cytochrome c oxidase (COX) consists of 13 subunits, 3 encoded in the mitochondrial genome and 10 in the nucleus. Little is known of the role of the nuclear-encoded subunits, some of which exhibit tissue-specific isoforms. Subunit VIa is unique in having tissue-specific isoforms in all mammalian species examined. We examined relative evolutionary rates for the COX6A heart (H) and liver (L) isoform genes along the length of the molecule, specifically in relation to the tissue-specific function(s) of the two isoforms. Nonsynonymous (amino acid replacement) substitutions in the COX6AH gene occurred more frequently than in the ubiquitously expressed COX6AL gene. Maximum-parsimony analysis and sequence divergences from reconstructed ancestral sequences revealed that after the ancestral COX6A gene duplicated to yield the genes for the H and L isoforms, the sequences encoding the mitochondrial matrix region of the COX VIa protein experienced an elevated rate of nonsynonymous substitutions relative to synonymous substitutions. This is expected for relaxed selective constraints after gene duplication followed by purifying selection to preserve the replacements with tissue-specific functions.      

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.     

Genomic characterization of POS5, the Saccharomyces cerevisiae mitochondrial NADH kinase

Disruption of the Saccharomyces cerevisiae mitochondrial NADH kinase POS5 increases the mitochondrial mutation rate 50-fold. Whereas most multicellular eukaryotic genomes have one NADH kinase gene, the yeast genome contains three distinct genes encoding NAD/H kinase activity. To determine if all three genes are essential for viability we constructed combinations of gene knockouts. We show that only the pos5Deltautr1Delta combination is synthetically lethal, demonstrating an essential overlapping function, and showing that NAD/H kinase activity is essential for eukaryotic viability. The single human NAD/H kinase gene can rescue the lethality of the double knockout in yeast, demonstrating that the single human gene can fill the various functions provided by the three yeast genes. The human NAD/H kinase gene harbors very common sequence variants, but all of these equally complement the synthetic lethality in yeast, illustrating that each of these are functionally wild-type. To understand the molecular mechanism of the mitochondrial genome instability of pos5 mutation we performed gene expression analysis on the pos5Delta. The pos5Delta resulted in an increase in expression of most of the iron transport genes including key genes involved in iron-sulfur cluster assembly. Decreased expression occurred in many genes involved in the electron transport chain. We show that the pos5Delta expression pattern is similar to the frataxin homolog knockout (yfh1Delta), the yeast model for Friedreich's ataxia. These combined data show that the POS5 NAD/H kinase is an important protein required for a variety of essential cellular pathways and that deficient iron-sulfur cluster assembly may play a critical role in the mitochondrial mutator phenotype observed in the pos5Delta.