Genome streamlining via complete loss of introns has occurred multiple times in lichenized fungal mitochondria

Reductions in genome size and complexity are a hallmark of obligate symbioses. The mitochondrial genome displays clear examples of these reductions, with the ancestral alpha‐proteobacterial genome size and gene number having been reduced by orders of magnitude in most descendent modern mitochondrial genomes. Here, we examine patterns of mitochondrial evolution specifically looking at intron size, number, and position across 58 species from 21 genera of lichenized Ascomycete fungi, representing a broad range of fungal diversity and niches. Our results show that the cox1gene always contained the highest number of introns out of all the mitochondrial protein‐coding genes, that high intron sequence similarity (>90%) can be maintained between different genera, and that lichens have undergone at least two instances of complete, genome‐wide intron loss consistent with evidence for genome streamlining via loss of parasitic, noncoding DNA, in Phlyctis boliviensisand Graphis lineola. Notably, however, lichenized fungi have not only undergone intron loss but in some instances have expanded considerably in size due to intron proliferation (e.g., Alectoria fallacina and Parmotrema neotropicum), even between closely related sister species (e.g., Cladonia). These results shed light on the highly dynamic mitochondrial evolution that is occurring in lichens and suggest that these obligate symbiotic organisms are in some cases undergoing recent, broad‐scale genome streamlining via loss of protein‐coding genes as well as noncoding, parasitic DNA elements.     

Footprints of unconventional mitochondrial inheritance in bivalve phylogeny: Signatures of positive selection on clades with doubly uniparental inheritance

The doubly uniparental inheritance (DUI) of some bivalve mollusks is the major exception to the common maternal inheritance of mitochondria in animals. DUI involves two mitochondrial lineages with paternal and maternal transmission routes, and it appears as a complex phenomenon requiring both nuclear and mitochondrial adaptations. DUI distribution seems to be scattered among the Bivalvia, and there are several clues for its multiple origins. In this paper, we investigate whether the incipient DUI systems had left possible selective signatures on mitochondrial genomes. Alongside the outstanding divergence of amino acid sequences, we confirmed strong purifying selection to act on mitochondrial genes. However, we found evidence that distinct episodes of intense directional pressure are associated with the origins of different DUI systems: We interpret these signals as footprints of the coevolution with the nuclear genome that ought to take place at the base of a DUI clade. Six genes (atp6, cox1, cox2, cox3, nad4L, and nad6) seem to be more commonly linked to the appearance of DUI. We also identified few putative DUI‐specific mutations, thus extending support to the hypothesis of multiple independent origins of this complex phenomenon.     

Characterization of the complete mitochondrial genome of Diaphania pyloalis (Lepidoptera: Pyralididae)

The complete mitochondrial genome (mitogenome) of Diaphania pyloalis (Lepidoptera: Pyralididae) was determined to be 15,298 bp and has the typical gene organization of mitogenomes from lepidopteran insects. It consists of 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and an A + T-rich region. The A + T content of this mitogenome is 80.83% and the AT skew is slightly positive. All PCGs are initiated by ATN codons, except for cytochrome c oxidase subunit 1 (cox1) gene which is initiated by CGA. Only the cox2 gene has an incomplete stop codon consisting of just a T. All the tRNA genes display a typical clover-leaf structure of mitochondrial tRNA. The A + T-rich region of the mitogenome is 332 bp in length, including several common features found in lepidopteran mitogenomes. Phylogenetic analysis showed that the D. pyloalis is close to Pyralididae.     

The complete mitochondrial genome of the verongid sponge Aplysina cauliformis: implications for DNA barcoding in demosponges

DNA “barcoding,” the determination of taxon-specific genetic variation typically within a fragment of the mitochondrial cytochrome oxidase 1 (cox1) gene, has emerged as a useful complement to morphological studies, and is routinely used by expert taxonomists to identify cryptic species and by non-experts to better identify samples collected during field surveys. The rate of molecular evolution in the mitochondrial genomes (mtDNA) of nonbilaterian animals (sponges, cnidarians, and placozoans) is much slower than in bilaterian animals for which DNA barcoding strategies were developed. If sequence divergence among nonbilaterian mtDNA and specifically cox1 is too slow to generate diagnostic variation, alternative genes for DNA barcoding and species-level phylogenies should be considered. Previous study across the Aplysinidae (Demospongiae, Verongida) family of sponges demonstrated no nucleotide substitutions in the traditional cox1 barcoding fragment among the Caribbean species of Aplysina. As the mitochondrial genome of Aplysina fulva has previously been sequenced, we are now able to make the first comparisons between complete mtDNA of congeneric demosponges to assess whether potentially informative variation exists in genes other than cox1. In this article, we present the complete mitochondrial genome of Aplysina cauliformis, a circular molecule 19620 bp in size. The mitochondrial genome of A. cauliformis is the same length as is A. fulva and shows six confirmed nucleotide differences and an additional 11 potential SNPs. Of the six confirmed SNPs, NADH dehydrogenase subunit 5 (nad5) and nad2 each contain two, and in nad2 both yield amino acid substitutions, suggesting balancing selection may act on this gene. Thus, while the low nucleotide diversity in Caribbean aplysinid cox1 extends to the entire mitochondrial genome, some genes do display variation. If these represent interspecific differences, then they may be useful alternative markers for studies in recently diverged sponge clades.     

肺囊康定产生菌Glarea lozoyensis线粒体基因组序列的再分析

[目的] Glarea lozoyensis是抗真菌药物卡泊芬净的产生菌,其突变菌株ATCC 74030的线粒体基因组已被报道。我们此前的研究发现诱变剂能引起该菌某些细胞核基因的突变,但诱变剂是否也能引起线粒体DNA序列的改变并不清楚。[方法] 组装野生型菌株ATCC 20868的线粒体基因组,并与发表的突变型菌株ATCC 74030的线粒体基因组进行比较。通过PCR验证野生和突变菌株线粒体基因组间表现差异之处,并利用正确的线粒体基因组序列进行新的分析。[结果] 我们成功组装出野生型菌株ATCC 20868的线粒体基因组,通过比较其与发表的ATCC 74030的线粒体基因组序列,发现存在6处单核苷酸变异位点和2处具有长度差异的区域。然而,随后的PCR验证和序列比较并没有发现2个菌株间存在这些差异。最初观察到的碱基差异是因为发表的ATCC 74030线粒体基因组存在序列错误。有趣的是,在Glarea lozoyensis的线粒体基因组中,我们发现存在3个具有内含子的tRNA基因和1个rnpB基因。同时,该菌线粒体基因组中存在多种重复序列,在其线粒体和细胞核基因组间也存在明显的DNA片段重复事件。[结论] 诱变剂没有引起G. lozoyensis线粒体DNA的任何改变;发表的ATCC 74030的线粒体基因组存在序列错误。我们报道G. lozoyensis正确的线粒体基因组序列,并且发现该菌线粒体和细胞核基因组间频繁的基因交流。     

Complete mitochondrial exploration of Echinococcus multilocularis from French alveolar echinococcosis patients

Alveolar echinococcosis (AE) is a parasitosis that is expanding worldwide, including in Europe. The development of genotypic markers is essential to follow its spatiotemporal evolution. Sequencing of the commonly used mitochondrial genes cob, cox1, and nad2 shows low discriminatory power, and analysis of the microsatellite marker EmsB does not allow nucleotide sequence analysis. We aimed to develop a new method for the genotyping of Echinococcus multilocularis based on whole mitochondrial genome (mitogenome) sequencing, to determine the genetic diversity among 30 human visceral samples from French patients, and compare this method with those currently in use. Sequencing of the whole mitochondrial genome was carried out after amplification by PCR, using one uniplex and two multiplex reactions to cover the 13,738 bp of the mitogenome, combined with Illumina technology. Thirty complete mitogenome sequences were obtained from AE lesions. One showed strong identity with Asian genotypes (99.98% identity) in a patient who had travelled to China. The other 29 mitogenomes could be differentiated into 13 haplotypes, showing higher haplotype and nucleotide diversity than when using the cob, cox1, and nad2 gene sequences alone. The mitochondrial genotyping data and EmsB profiles did not overlap, probably because one method uses the mitochondrial genome and the other the nuclear genome. The pairwise fixation index (F_st) value between individuals living inside and those living outside the endemic area was high (F_st = 0.222, P = 0.002). This is consistent with the hypothesis of an expansion from historical endemic areas to peripheral regions.     

The complete nucleotide sequence of the mitochondrial genome of <Emphasis Type="Italic">Phthonandria atrilineata</Emphasis> (Lepidoptera: Geometridae)

Using long-polymerase chain reaction (Long-PCR) method, we determined the complete nucleotide sequence of the mitochondrial genome (mitogenome) of Phthonandria atrilineata. The complete mtDNA from P. atrilineata was 15,499 base pairs in length and contained 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes, and a control region. The P. atrilineata genes were in the same order and orientation as the completely sequenced mitogenomes of other lepidopteran species. The nucleotide composition of P. atrilineata mitogenome was biased toward A + T nucleotides (81.02%), and the 13 PCGs show different A + T contents that range from 73.25% (cox1) to 92.12% (atp8). Phthonandria had the canonical set of 22 tRNA genes, that fold in the typical cloverleaf structure described for metazoan mt tRNAs, with the unique exception of trnS(AGN). The phylogenetic relationships were reconstructed with the concatenated sequences of the 13 PCGs of the mitochondrial genome, which confirmed that P. atrilineata is most closely related to the superfamily Bombycoidea.     

Identification of the Mitochondrial Genome in the Chrysophyte Alga Ochromonas danica

ABSTRACT. Analysis of total DNA isolated from the Chrysophyte alga Ochromonas danica revealed, in addition to nuclear DNA, two genomes present as numerous copies per cell. The larger genome (?120 kilobase pairs or kbp) is the plastid DNA, which is identified by its hybridization to plasmids containing sequences for the photosynthesis genes rbcL, psbA, and psbC. The smaller genome (40 kbp) is the mitochondrial genome as identified by its hybridization with plasmids containing gene sequences of plant cytochrome oxidase subunits I and II. Both the 120- and 40-kbp genomes contain genes for the small and large subunits of rDNA. The mitochondrial genome is linear with terminal inverted repeats of about 1.6 kbp. Two other morphologically similar species were examined, Ochromonas minuta and Poteriochromonas malhamensis. All three species have linear mitochondrial DNA of 40 kbp. Comparisons of endonuclease restriction-fragment patterns of the mitochondrial and chloroplast DNAs as well as those of their nuclear rDNA repeats failed to reveal any fragment shared by any two of the species. Likewise, no common fragment size was detected by hybridization with plasmids containing heterologous DNA or with total mitochondrial DNA of O. danica; these observations support the taxonomic assignment of these three organisms to different species. The Ochromonas mitochondrial genomes are the first identified in the chlorophyll a/c group of algae. Combining these results with electron microscopic observations of putative mitochondrial genomes reported for other chromophytes and published molecular studies of other algal groups suggests that all classes of eukaryote algae may have mitochondrial genomes < 100 kbp in size, more like other protistans than land plants.     

Rapid structural evolution of Dendrobium mitogenomes and mito-nuclear phylogeny discordances in Dendrobium (Orchidaceae)

Reconstructing mitochondrial genomes of angiosperms is extremely intricate due to frequent recombinations which give rise to varied sized in Dendrobium mitogenomes and their structural variations, even in most orchid species. In this study, we first sequenced two complete and five draft mitochondrial genomes of Dendrobium using next-generation and third-generation sequencing technologies. The mitochondrial genomes were 420 538–689 048 bp long, showing multipartite (multichromosomal) structures that consisted of variably sized circular or linear-mapping isoforms (chromosomes). The comparison of mitochondrial genomes showed frequent gene losses in Dendrobium species. To explore structure variations of mitochondrial genomes in vivo, we quantified copy numbers of five mitochondrial genes and DNA contents per mitochondrion. The gene copy numbers and the DNA contents showed extreme differences during Dendrobium development, suggesting dynamic structures of mitochondrial genomes. Furthermore, phylogenetic relationships of 97 accessions from 39 Dendrobium species were constructed based on 12 nuclear single-copy genes and 15 mitochondrial genes. We discovered obvious discordance between the nuclear and mitochondrial trees. Reticulate evolution was inferred from the species network analysis in Dendrobium. Our findings revealed the rapid structural evolution of Dendrobium mitochondrial genomes and the existence of hybridization events in Dendrobium species, which provided new insights into in vivo structural variations of plant mitochondrial genomes and the strong potential of mitochondrial genes in deciphering plant evolution history.     

The complete mitochondrial genome of fall armyworm <Emphasis Type="Italic">Spodoptera frugiperda</Emphasis> (Lepidoptera:Noctuidae)

The mitochondrial genome (mitogenome) can provide important information for understanding molecular evolution and phylogenetic analyses. The complete mitogenome of Spodoptera frugiperda (Lepidoptera:Noctuidae) was determined to be 15,365 bp in length and has the typical gene order found in Noctuidae mitogenomes, it includes 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a A+T-rich region. The nucleotide composition was biased toward A+T nucleotides (81.09 %) and the AT skew of this mitogenome was slightly positive (0.004). All PCGs were initiated by ATN codons, except for cytochrome c oxidase subunit 1 (cox1) gene which was initiated by CGA. Eight of the 13 PCGs have the incomplete termination codon, T or TA. All the tRNA genes displayed the typical clover-leaf structure of mitochondrial tRNAs, with the exception of trnS1 (AGN). The A+T-rich region was 328 bp in length and consisted of several features common to the Noctuidae insects. Phylogenetic analysis showed that the S. frugiperda was within the Noctuidae.     

Characterization of the complete mitochondrial genome of Cryptotermes domesticus (Blattodea: Kalotermitidae): Genome description and phylogenetic implications

The complete mitochondrial genome of Cryptotermes domesticus (Haviland) was sequenced and annotated to study its characteristics and the phylogenetic relationship of C. domesticus to other termite species. The mitogenome of C. domesticus is a circular, close, and double-stranded molecule with a length of 15,655 bp. The sequenced mitogenome contains 37 typical genes, which are highly conserved in gene size, organization, and codon usage. Transfer RNA genes (tRNAs) also have typical secondary structures. All of the 13 protein-coding genes (PCGs) start with an ATN codon, except for nad4, which starts with GTG and terminates with the terminal codon TAA and TAG or the incomplete form T-- (cox2 and nad5). Most tRNAs have a typical cloverleaf structure, except for trnS1, in which this form is replaced by a simple loop and lacks the dihydrouridine (DHU) arm. The nucleotide diversity (Pi) and nonsynonymous (Ka)/synonymous (Ks) mutation rate ratios indicate that nad1, cox1, and cox3 are the most conserved genes, and that cox1 has the lowest rate of evolution. In addition, an 89 bp repeated sequence was found in the A + T-rich region. Phylogenetic analysis was performed using Bayesian inference (BI) and maximum likelihood (ML) methods based on 13 PCGs, and the monophyly of Kalotermitidae was supported.     

The complete mitogenome of Bombyx mori strain Dazao (Lepidoptera: Bombycidae) and comparison with other lepidopteran insects

The complete mitochondrial genome (mitogenome) of Bombyx mori strain Dazao (Lepidoptera: Bombycidae) was determined to be 15,653 bp, including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a A + T-rich region. It has the typical gene organization and order of mitogenomes from lepidopteran insects. The AT skew of this mitogenome was slightly positive and the nucleotide composition was also biased toward A + T nucleotides (81.31%). All PCGs were initiated by ATN codons, except for cytochrome c oxidase subunit 1 (cox1) gene which was initiated by CGA. The cox1 and cox2 genes had incomplete stop codons consisting of just a T. All the tRNA genes displayed a typical clover-leaf structure of mitochondrial tRNA. The A + T-rich region of the mitogenome was 495 bp in length and consisted of several features common to the lepidopteras. Phylogenetic analysis showed that the B. mori Dazao was close to Bombycidae.     

Regular Spliceosomal Introns Are Invasive in Chlamydomonas reinhardtii: 15 Introns in the Recently Relocated Mitochondrial cox2 and cox3 Genes

In the unicellular green alga, Chlamydomonas reinhardtii, cytochrome oxidase subunit 2 (cox2) and 3 (cox3) genes are missing from the mitochondrial genome. We isolated and sequenced a BAC clone that carries the whole cox3 gene and its corresponding cDNA. Almost the entire cox2 gene and its cDNA were also determined. Comparison of the genomic and the corresponding cDNA sequences revealed that the cox3 gene contains as many as nine spliceosomal introns and that cox2 bears six introns. Putative mitochondria targeting signals were predicted at each N terminal of the cox genes. These spliceosomal introns were typical GT–AG-type introns, which are very common not only in Chlamydomonas nuclear genes but also in diverse eukaryotic taxa. We found no particular distinguishing features in the cox introns. Comparative analysis of these genes with the various mitochondrial genes showed that 8 of the 15 introns were interrupting the conserved mature protein coding segments, while the other 7 introns were located in the N-terminal target peptide regions. Phylogenetic analysis of the evolutionary position of C. reinhardtii in Chlorophyta was carried out and the existence of the cox2 and cox3 genes in the mitochondrial genome was superimposed in the tree. This analysis clearly shows that these cox genes were relocated during the evolution of Chlorophyceae. It is apparent that long before the estimated period of relocation of these mitochondrial genes, the cytosol had lost the splicing ability for group II introns. Therefore, at least eight introns located in the mature protein coding region cannot be the direct descendant of group II introns. Here, we conclude that the presence of these introns is due to the invasion of spliceosomal introns, which occurred during the evolution of Chlorophyceae. This finding provides concrete evidence supporting the ``intron-late' model, which rests largely on the mobility of spliceosomal introns. Received: 22 August 2000 / Accepted: 28 February 2001     

球孢白僵菌种内比较线粒体基因组学分析

【目的】明确球孢白僵菌种内线粒体基因组的分化程度。【方法】从GenBank下载已知的球孢白僵菌6个菌株线粒体基因组序列,详细分析基因组的组成结构,比较外显子区、内含子区和基因间区的碱基变异情况,分析菌株间的系统发育关系。【结果】球孢白僵菌不同菌株的线粒体基因组大小为28.8–32.3 kb,都有14个常见的核心蛋白编码基因、2个rRNA基因和25个tRNA基因,具有很强的共线性关系。但是,不同菌株含有的线粒体内含子数目存在差异(2–5个/菌株),在cox1、cox2和nad1基因中表现出内含子插入/缺失多态性,这是导致线粒体基因组大小变化的主要因素。对外显子、内含子和基因间区的碱基变异情况进行分析,发现内含子和基因间区相对变异较大,而外显子区相对变异较小。系统发育分析发现,这些球孢白僵菌菌株以很高的支持度聚在一起,具有相同内含子分布规律的菌株也具有较近的聚类关系。【结论】本研究首次报道球孢白僵菌因内含子数目不同、插入缺失突变和单核苷酸变异等在线粒体基因组上表现出较大程度的遗传分化,为认识真菌种内线粒体基因组分化提供了新的证据。     

Promiscuous mitochondrial group II intron sequences in plant nuclear genomes

Gene translocations from the organelles to the nucleus are postulated by the endosymbiont hypothesis. We here report evidence for sequence insertions in the nuclear genomes of plants that are derived from noncoding regions of the mitochondrial genome. Fragments of mitochondrial group II introns are identified in the nuclear genomes of tobacco and a bean species. The duplicated intron sequences of 75–140 bp are derived from cis- and trans-splicing introns of genes encoding subunits 1 and 5 of the NADH dehydrogenase. The mitochondrial sequences are inserted in the vicinities of a lectin gene, different glucanase genes and a gene encoding a subunit of photosystem II. Sequence similarities between the nuclear and mitochondrial copies are in the range of 80 to 97%, suggesting recent transfer events that occurred in the basic glucanase genes before and in the lectin gene after the gene duplications in the evolution of the nuclear gene families. Overlapping regions of the same introns are in two instances also involved in intramitochondrial sequence duplications. Correspondence to: V. Knoop     

Jumping and gliding rodents: Mitogenomic affinities of Pedetidae and Anomaluridae deduced from an RNA-Seq approach

An RNA-Seq strategy was used to obtain the complete set of protein-coding mitochondrial genes from two rodent taxa. Thanks to the next generation sequencing (NGS) 454 approach, we determined the complete mitochondrial DNA genome from Graphiurus kelleni (Mammalia: Rodentia: Gliridae) and partial mitogenome from Pedetes capensis (Pedetidae), and compared them with published rodent and outgroup mitogenomes. We finished the mitogenome sequencing by a series of amplicons using conserved PCR primers to fill the gaps corresponding to tRNA, rRNA and control regions. Phylogenetic analyses of the mitogenomes suggest a well-supported rodent phylogeny in agreement with nuclear gene trees. Pedetes groups with Anomalurus into the clade Anomaluromorpha, while Graphiurus branches within the squirrel-related clade. Moreover, Pedetes + Anomalurus branch with Castor into the mouse-related clade. Our study demonstrates the utility of NGS for obtaining new mitochondrial genomes as well as the importance of choosing adequate models of sequence evolution to infer the phylogeny of rodents.     

Complete mitochondrial genome of the darkling beetle Gonocephalum outreyi (Coleoptera: Tenebrionidae) with phylogenetic implications

The complete mitochondrial genome (mitogenome) of Gonocephalum outreyi was determined by using next-generation sequencing approach. The full length of this mitogenome is 15,836?bp, which consists of 37 typical metazoan mitochondrial genes with an identical genome organization to ancestral insects. The majority of the protein-coding genes begin with the codon ATN, except for cox1 and cox2 with AAT and AAA, respectively. To elucidate the phylogenetic position of G. outreyi, we used various sequence coding schemes for protein-coding genes and the combined nucleotide sequences of all mitochondrial genes for tree building under the Bayesian and Maximum Likelihood inferences. The phylogenetic results consistently supported G. outreyi as a member of the family Tenebrionidae. The monophyly of both Tenebrionoidea and Tenebrionidae were strongly supported. The Scraptiidae and Melandryidae were recovered to be non-monophyletic in regards to the Osphya. Within Tenebrionidae, the subfamilies Diaperinae and Tenebrioninae were found to be non-monophyletic.     

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.