The mitochondrial genome of the fission yeast Schizosaccharomyces pombe. 2. Localization of genes by interspecific hybridization in strain ade7-50h- and cloning of the genome in small fragments

A series of 18 small overlapping restriction fragments has been cloned, covering the complete mitochondrial genome of Schizosaccharomyces pombe. By hybridizing mitochondrial gene probes from Saccharomyces cerevisiae and Neurospora crassa with restriction fragments of Schizosaccharomyces pombe mitochondrial DNA, the following homologous genes were localized on the mitochondrial genome of S. pombe: cob, cox1, cox2 and cox3, ATPase subunit 6 and 9 genes, the large rRNA gene and both types of open reading frames occurring in mitochondrial introns of various ascomycetes. The region of the genome, hybridizing with cob exon probes is separated by an intervening sequence of about 2500 bp, which is homologous with the first two introns of the cox1 gene in Saccharomyces cerevisiae (class II introns according to Michel et al. 1982). Similarly, in the cox1 homologous region, which covers about 4000 bp, two regions were detected hybridizing with class I intron probes, suggesting the existence of two cox1 introns in Schizosaccharomyces pombe. Hybridization with several specific exon probes with a determined order has revealed that cob, cox1, cox3 and the large rRNA gene are all transcribed from the same DNA strand. The low intensities of hybridization signals suggest a large evolutionary distance between Schizosaccharomyces pombe and Saccharomyces cerevisiae or Neurospora crassa mitochondrial genes. Considering the length of the mitochondrial DNA of Schizosaccharomyces pombe (about 19.4 kbp) and the expected length of the localized genes and intron sequences there is enough space left for encoding the expected set of tRNAs and the small rRNA gene. The existence of leader-, trailer-, ori- and spacer sequences or further unassigned reading frames is then restricted to a total length of about 3000 bp only.     

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.     

小长蝽线粒体基因组全序列测定与长蝽总科系统发育分析

为了解小长蝽Nysius ericae(Schilling)线粒体基因组结构及长蝽总科的分子系统发育关系。本试验采用Illumina MiSeq测序平台对小长蝽线粒体基因进行测序,对基因组序列进行拼装、注释和特征分析,利用最大似然法和贝叶斯法构建基于12种长蝽总科昆虫线粒体全基因组核苷酸序列的系统发育树。小长蝽线粒体基因组全长为16 330 bp(GenBank登录号：MW465654),基因组包括13个蛋白编码基因(PCGs),22个tRNA基因,2个rRNA基因和1段非编码控制区。11个蛋白质编码基因的起始密码子为典型的ATN;cox1,nad4l的起始密码子为TTG。cob的终止密码子为TAG,其余蛋白编码基因的终止密码子为TAA。只有trnS1缺少DHU臂,其余tRNA基因均能形成典型的三叶草结构。12种长蝽总科昆虫线粒体全基因组序列构建的昆虫系统发育树结果显示,小长蝽与Nysius plebeius具有更近的亲缘关系,且与传统形态学分类基本一致。小长蝽线粒体基因组符合长蝽总科线粒体基因组的一般特征。结果表明小长蝽与N.plebeius的亲缘关系更近。     

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.     

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

The Agrocybe aegerita mitochondrial genome contains two inverted repeats of the nad4 gene arisen by duplication on both sides of a linear plasmid integration site.

The Agrocybe aegerita mitochondrial genome possesses two polB genes with linear plasmid origin. The cloning and sequencing of the regions flanking Aa-polB P1 revealed two large inverted repeats (higher than 2421 nt) separated by a single copy region of 5834 nt. Both repeats contain identical copies of the nad4 gene. The single copy region contains two disrupted genes with plasmid origin Aa-polB P1 and a small ORF homologous to a small gene described in two basidiomycete linear plasmids. The phylogenetic analyses argue in favor of a same plasmid origin for both genes but, surprisingly, these genes were separated by a mitochondrial tRNA-Met. Both strands of the complete region containing the two nad4 inverted copies and the tRNA-Met appear to be transcribed on large polycistronic mRNAs. A model summarizing the events that would have occurred is proposed: (1) capture of the tRNA by the plasmid before its integration in the mtDNA or acquisition of the tRNA gene by recombination after the plasmid integration, (2) integration of the plasmid in the mtDNA, accompanied by a large duplication containing the nad4 gene and (3) erosion of the plasmid sequences by large deletions and mutations.     

The Agrocybe aegerita mitochondrial genome contains two inverted repeats of the nad4 gene arisen by duplication on both sides of a linear plasmid integration site

The Agrocybe aegerita mitochondrial genome possesses two polB genes with linear plasmid origin. The cloning and sequencing of the regions flanking Aa-polB P1 revealed two large inverted repeats (higher than 2421 nt) separated by a single copy region of 5834 nt. Both repeats contain identical copies of the nad4 gene. The single copy region contains two disrupted genes with plasmid origin Aa-polB P1 and a small ORF homologous to a small gene described in two basidiomycete linear plasmids. The phylogenetic analyses argue in favor of a same plasmid origin for both genes but, surprisingly, these genes were separated by a mitochondrial tRNA-Met. Both strands of the complete region containing the two nad4 inverted copies and the tRNA-Met appear to be transcribed on large polycistronic mRNAs. A model summarizing the events that would have occurred is proposed: (1) capture of the tRNA by the plasmid before its integration in the mtDNA or acquisition of the tRNA gene by recombination after the plasmid integration, (2) integration of the plasmid in the mtDNA, accompanied by a large duplication containing the nad4 gene and (3) erosion of the plasmid sequences by large deletions and mutations.     

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.     

Characterization of the Complete Chloroplast Genome of Apple (Malus × domestica,Rosaceae)*

Apple (Malus × domestica) is one of the most important temperate fruits. To better understand the molecular basis of this species, we characterized the complete chloroplast (cp) genome sequence downloaded from Genome Database for Rosaceae. The cp genome of apple is a circular molecule of 160068bp in length with a typical quadripartite structure of two inverted repeats (IRs) of 26352bp, separated by a small single copy region of 19180bp (SSC) and a large single copy region (LSC) of 88184bp. A total of 135 predicted genes (115 unique genes, and another 20 genes were duplicated in the IR) were identified, including 81 protein coding genes, four rRNA genes and 30 tRNA genes. Three genes of ycf15, ycf68 and infA contain several internal stop codons, which were interpreted as pseudogenes. The genome structure, gene order, GC content and codon usage of apple are similar to the typical angiosperm cp genomes. Thirty repeat regions (≥30bp) were detected, twenty one of which are tandem, six are forward and three are inverted repeats. Two hundred thirty seven simple sequence repeat (SSR) loci were revealed and most of them are composed of A or T, contributing to a distinct bias in base composition. Additionally, average 10000bp non coding region contains 24 SSR sites, while protein coding region contains five SSR sites, indicating an uneven distribution of SSRs. The complete cp genome sequence of apple reported in this paper will facilitate the future studies of its population genetics, phylogenetics and chloroplast genetic engineering.     

Separate origins of group I introns in two mitochondrial genes of the katablepharid Leucocryptos marina

Mitochondria are descendants of the endosymbiotic α-proteobacterium most likely engulfed by the ancestral eukaryotic cells, and the proto-mitochondrial genome should have been severely streamlined in terms of both genome size and gene repertoire. In addition, mitochondrial (mt) sequence data indicated that frequent intron gain/loss events contributed to shaping the modern mt genome organizations, resulting in the homologous introns being shared between two distantly related mt genomes. Unfortunately, the bulk of mt sequence data currently available are of phylogenetically restricted lineages, i.e., metazoans, fungi, and land plants, and are insufficient to elucidate the entire picture of intron evolution in mt genomes. In this work, we sequenced a 12 kbp-fragment of the mt genome of the katablepharid Leucocryptos marina. Among nine protein-coding genes included in the mt genome fragment, the genes encoding cytochrome b and cytochrome c oxidase subunit I (cob and cox1) were interrupted by group I introns. We further identified that the cob and cox1 introns host open reading frames for homing endonucleases (HEs) belonging to distantly related superfamilies. Phylogenetic analyses recovered an affinity between the HE in the Leucocryptos cob intron and two green algal HEs, and that between the HE in the Leucocryptos cox1 intron and a fungal HE, suggesting that the Leucocryptos cob and cox1 introns possess distinct evolutionary origins. Although the current intron (and intronic HE) data are insufficient to infer how the homologous introns were distributed to distantly related mt genomes, the results presented here successfully expanded the evolutionary dynamism of group I introns in mt genomes.     

Complete mitochondrial genome of the boreal digging frog Kaloula borealis (Anura, Microhylidae)

We sequenced the complete mitochondrial genome from the boreal digging frog Kaloula borealis. The genome sequence was 17,173 bp in size, and the gene order and contents were identical to those of previously reported amphibian mitochondrial genomes. Of 13 protein-coding genes (PCGs), 5 genes (CO2, ATPase 6, CO3, ND3, and ND4) had incomplete stop codons. Also ND1 gene used GTG as a start codon, while CO1 and ND5 genes used AGG as a stop codon. The base composition of K. borealis mitogenome showed a strong anti-G bias (6.11%) on the 3rd position of PCGs.     

Mitochondrial genome of a tertiary endosymbiont retains genes for electron transport proteins

Mitochondria and plastids originated through endosymbiosis, and subsequently became reduced and integrated with the host in similar ways. Plastids spread between lineages through further secondary or even tertiary endosymbioses, but mitochondria appear to have originated once and have not spread between lineages. Mitochondria are also generally lost in secondary and tertiary endosymbionts, with the single exception of the diatom tertiary endosymbiont of dinoflagellates like Kryptoperidinium foliaceum, where both host and endosymbiont are reported to contain mitochondria. Here we describe the first mitochondrial genes from this system: cytochrome c oxidase 1 (cox1), cytochrome oxidase 3 (cox3), and cytochrome b (cob). Phylogenetic analyses demonstrated that all characterized genes were derived from the pennate diatom endosymbiont, and not the host. We also demonstrated that all three genes are expressed, that cox1 contains spliced group II introns, and that cob and cox3 form an operon, all like their diatom relatives. The endosymbiont mitochondria not only retain a genome, but also express their genes, and are therefore likely involved in electron transport. Ultrastructural examination confirmed the endosymbiont mitochondria retain normal tubular cristae. Overall, these data suggest the endosymbiont mitochondria have not reduced at the genomic or functional level.     

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

Evolutionary rates, divergence dates, and the performance of mitochondrial genes in Bayesian phylogenetic analysis

The mitochondrial genome is one of the most frequently used loci in phylogenetic and phylogeographic analyses, and it is becoming increasingly possible to sequence and analyze this genome in its entirety from diverse taxa. However, sequencing the entire genome is not always desirable or feasible. Which genes should be selected to best infer the evolutionary history of the mitochondria within a group of organisms, and what properties of a gene determine its phylogenetic performance? The current study addresses these questions in a Bayesian phylogenetic framework with reference to a phylogeny of plethodontid and related salamanders derived from 27 complete mitochondrial genomes; this topology is corroborated by nuclear DNA and morphological data. Evolutionary rates for each mitochondrial gene and divergence dates for all nodes in the plethodontid mitochondrial genome phylogeny were estimated in both Bayesian and maximum likelihood frameworks using multiple fossil calibrations, multiple data partitions, and a clock-independent approach. Bayesian analyses of individual genes were performed, and the resulting trees compared against the reference topology. Ordinal logistic regression analysis of molecular evolution rate, gene length, and the G-shape parameter a demonstrated that slower rate of evolution and longer gene length both increased the probability that a gene would perform well phylogenetically. Estimated rates of molecular evolution vary 84-fold among different mitochondrial genes and different salamander lineages, and mean rates among genes vary 15-fold. Despite having conserved amino acid sequences, cox1, cox2, cox3, and cob have the fastest mean rates of nucleotide substitution, and the greatest variation in rates, whereas rrnS and rrnL have the slowest rates. Reasons underlying this rate variation are discussed, as is the extensive rate variation in cox1 in light of its proposed role in DNA barcoding.     

The mitochondrial genome of Strongyloides stercoralis (Nematoda) - idiosyncratic gene order and evolutionary implications

The complete mitochondrial genome sequence of the parasitic nematode Strongyloides stercoralis was determined, and its organisation and structure compared with other nematodes for which complete mitochondrial sequence data were available. The mitochondrial genome of S. stercoralis is 13,758 bp in size and contains 36 genes (all transcribed in the clockwise direction) but lacks the atp8 gene. This genome has a high T content (55.9%) and a low C content (8.3%). Corresponding to this T content, there are 16 (poly-T) tracts of >/=12 Ts distributed across the genome. In protein-coding genes, the T bias is greatest (76.4%) at the third codon position compared with the first and second codon positions. Also, the C content is higher at the first (9.3%) and second (13.4%) codon positions than at the third (2%) position. These nucleotide biases have a significant effect on predicted codon usage patterns and, hence, on amino acid compositions of the mitochondrial proteins. Interestingly, six of the 12 protein-coding genes are predicted to employ a unique initiation codon (TTT), which has not yet been reported for any other animal mitochondrial genome. The secondary structures predicted for the 22 transfer RNA (trn) genes and the two ribosomal RNA (rrn) genes are similar to those of other nematodes. In contrast, the gene arrangement in the mitochondrial genome of S. stercoralis is different from all other nematodes studied to date, revealing only a limited number of shared gene boundaries (atp6-nad2 and cox2-rrnL). Evolutionary analyses of mitochondrial nucleotide and amino acid sequence data sets for S. stercoralis and seven other nematodes demonstrate that the mitochondrial genome provides a rich source of phylogenetically informative characters. In conclusion, the S. stercoralis mitochondrial genome, with its unique gene order and characteristics, should provide a resource for comparative mitochondrial genomics and systematics studies of parasitic nematodes.