Mitochondrial genomes: anything goes

Mt-rps3 is an ancient gene which provides insight into the evolution of fungal mitochondrial genomes

The nuclear ribosomal protein S3 (Rps3) is implicated in the assembly of the ribosomal small subunit. Fungi and plants present a gene copy in their mitochondrial (mt) genomes. An analysis of 303 complete fungal mt genomes showed that, when rps3 is found, it is either a free-standing gene or an anchored gene within the omega intron of the rnl gene. Early divergent fungi, Basidiomycota and all yeasts but the CTG group belong to the first case, and Pezizomycotina to the second. Its position, size and genetic code employed are conserved within species of the same Order. Size variability is attributed to different number of repeats. These repeats consist of AT-rich sequences. MtRps3 proteins lack the KH domain, necessary for binding to rRNA, in their N-terminal region. Their C-terminal region is conserved in all Domains of life. Phylogenetic analysis showed that nuclear and mtRps3 proteins are descendants of archaeal and a-proteobacterial homologues, respectively. Thus, fungal mt-rps3 gene is an ancient gene which evolved within the endosymbiotic model and presents different evolutionary routes: (a) coming from a-proteobacteria, it was relocated to another region of the mt genome, (b) via its insertion to the omega intron, it was transferred to the nucleus and/or got lost, and (c) it was re-routed to the mt genome again. Today, Basidiomycota and Saccharomycetales seem to follow the first evolutionary route and almost all Pezizomycotina support the second scenario with their exceptions being the result of the third scenario, i.e., the gene's re-entry to the mt genome.

The bipartite mitochondrial genome of Ruizia karukerae (Rhigonematomorpha, Nematoda)

Mitochondrial genes and whole mitochondrial genome sequences are widely used as molecular markers in studying population genetics and resolving both deep and shallow nodes in phylogenetics. In animals the mitochondrial genome is generally composed of a single chromosome, but mystifying exceptions sometimes occur. We determined the complete mitochondrial genome of the millipede-parasitic nematode Ruizia karukerae and found its mitochondrial genome consists of two circular chromosomes, which is highly unusual in bilateral animals. Chromosome I is 7,659 bp and includes six protein-coding genes, two rRNA genes and nine tRNA genes. Chromosome II comprises 7,647 bp, with seven protein-coding genes and 16 tRNA genes. Interestingly, both chromosomes share a 1,010 bp sequence containing duplicate copies of cox2 and three tRNA genes (trnD, trnG and trnH), and the nucleotide sequences between the duplicated homologous gene copies are nearly identical, suggesting a possible recent genesis for this bipartite mitochondrial genome. Given that little is known about the formation, maintenance or evolution of abnormal mitochondrial genome structures, R. karukerae mtDNA may provide an important early glimpse into this process.

Biogenesis of the bc1 Complex of the Mitochondrial Respiratory Chain

The oxidative phosphorylation system contains four respiratory chain complexes that connect the transport of electrons to oxygen with the establishment of an electrochemical gradient over the inner membrane for ATP synthesis. Due to the dual genetic source of the respiratory chain subunits, its assembly requires a tight coordination between nuclear and mitochondrial gene expression machineries. In addition, dedicated assembly factors support the step-by-step addition of catalytic and accessory subunits as well as the acquisition of redox cofactors. Studies in yeast have revealed the basic principles underlying the assembly pathways. In this review, we summarize work on the biogenesis of the bc₁ complex or complex III, a central component of the mitochondrial energy conversion system.

Comparative mitochondrial genomics of cryptophyte algae: gene shuffling and dynamic mobile genetic elements

BACKGROUND

Cryptophytes are an ecologically important group of algae comprised of phototrophic, heterotrophic and osmotrophic species. This lineage is of great interest to evolutionary biologists because their plastids are of red algal secondary endosymbiotic origin. Cryptophytes have a clear phylogenetic affinity to heterotrophic eukaryotes and possess four genomes: host-derived nuclear and mitochondrial genomes, and plastid and nucleomorph genomes of endosymbiotic origin.

RESULTS

To gain insight into cryptophyte mitochondrial genome evolution, we sequenced the mitochondrial DNAs of five species and performed a comparative analysis of seven genomes from the following cryptophyte genera: Chroomonas, Cryptomonas, Hemiselmis, Proteomonas, Rhodomonas, Storeatula and Teleaulax. The mitochondrial genomes were similar in terms of their general architecture, gene content and presence of a large repeat region. However, gene order was poorly conserved. Characteristic features of cryptophyte mtDNAs included large syntenic clusters resembling α-proteobacterial operons that encode bacteria-like rRNAs, tRNAs, and ribosomal protein genes. The cryptophyte mitochondrial genomes retain almost all genes found in many other eukaryotes including the nad, sdh, cox, cob, and atp genes, with the exception of sdh2 and atp3. In addition, gene cluster analysis showed that cryptophytes possess a gene order closely resembling the jakobid flagellates Jakoba and Reclinomonas. Interestingly, the cox1 gene of R. salina, T. amphioxeia, and Storeatula species was found to contain group II introns encoding a reverse transcriptase protein, as did the cob gene of Storeatula species CCMP1868.

CONCLUSIONS

These newly sequenced genomes increase the breadth of data available from algae and will aid in the identification of general trends in mitochondrial genome evolution. While most of the genomes were highly conserved, extensive gene arrangements have shuffled gene order, perhaps due to genome rearrangements associated with hairpin-containing mobile genetic elements, tRNAs with palindromic sequences, and tandem repeat sequences. The cox1 and cob gene sequences suggest that introns have recently been acquired during cryptophyte evolution. Comparison of phylogenetic trees based on plastid and mitochondrial genome data sets underscore the different evolutionary histories of the host and endosymbiont components of present-day cryptophytes.

Mitochondrial versus nuclear gene expression and membrane protein assembly: the case of subunit 2 of yeast cytochrome c oxidase

Deletion of the yeast mitochondrial gene COX2, encoding subunit 2 (mtCox2) of cytochrome c oxidase (CcO), results in a respiratory-incompetent Δcox2 strain. For a cytosol-synthesized Cox2 to restore respiratory growth, it must carry the W56R mutation (cCox2W56R). Nevertheless, only a fraction of cCox2W56R is matured in mitochondria, allowing ∼60% steady-state accumulation of CcO. This can be attributed either to the point mutation or to an inefficient biogenesis of cCox2W56R. We generated a strain expressing the mutant protein mtCox2W56R inside mitochondria which should follow the canonical biogenesis of mitochondria-encoded Cox2. This strain exhibited growth rates, CcO steady-state levels, and CcO activity similar to those of the wild type; therefore, the efficiency of Cox2 biogenesis is the limiting step for successful allotopic expression. Upon coexpression of cCox2W56R and mtCox2, each protein assembled into CcO independently from its genetic origin, resulting in a mixed population of CcO with most complexes containing the mtCox2 version. Notably, the presence of the mtCox2 enhances cCox2W56R incorporation. We provide proof of principle that an allotopically expressed Cox2 may complement a phenotype due to a mutant mitochondrial COX2 gene. These results are relevant to developing a rational design of genes for allotopic expression intended to treat human mitochondrial diseases.

Comparative analysis of the mitochondrial genome of the fungus Colletotrichum lindemuthianum, the causal agent of anthracnose in common beans

Fungi of the genus Colletotrichum are economically important and are used as models in plant-pathogen interaction studies. In this study, the complete mitochondrial genomes of two Colletotrichum lindemuthianum isolates were sequenced and compared with the mitochondrial genomes of seven species of Colletotrichum. The mitochondrial genome of C. lindemuthianum is a typical circular molecule 37,446 bp (isolate 89 A₂ 2-3) and 37,440 bp (isolate 83.501) in length. The difference of six nucleotides between the two genomes is the result of a deletion in the ribosomal protein S3 (rps3) gene in the 83.501 isolate. In addition, substitution of adenine for guanine within the rps3 gene in the mitochondrial genome of the 83.501 isolate was observed. Compared to the previously sequenced C. lindemuthianum mitochondrial genome, an exon no annotated in the cytochrome c oxidase I (cox1) gene and a non-conserved open reading frame (ncORF) were observed. The size of the mitochondrial genomes of the seven species of Colletotrichum was highly variable, being attributed mainly to the ncORF, ranging from one to 10 and also from introns ranging from one to 11 and which encode a total of up to nine homing endonucleases. This paper reports for the first time by means of transcriptome that then ncORFs are transcribed in Colletotrichum spp. Phylogeny data revealed that core mitochondrial genes could be used as an alternative in phylogenetic relationship studies in Colletotrichum spp. This work contributes to the genetic and biological knowledge of Colletotrichum spp., which is of great economic and scientific importance.

Replication stress in mitochondria

Mitochondrial DNA (mtDNA), which is essential for mitochondrial and cell function, is replicated and transcribed in the organelle by proteins that are entirely coded in the nucleus. Replication of mtDNA is challenged not only by threats related to the replication machinery and orchestration of DNA synthesis, but also by factors linked to the peculiarity of this genome. Indeed the architecture, organization, copy number, and location of mtDNA, which are markedly distinct from the nuclear genome, require ad hoc and complex regulation to ensure coordinated replication. As a consequence sub-optimal mtDNA replication, which results from compromised regulation of these factors, is generally associated with mitochondrial dysfunction and disease. Mitochondrial DNA replication should be considered in the context of the organelle and the whole cell, and not just a single genome or a single replication event. Major threats to mtDNA replication are linked to its dependence on both mitochondrial and nuclear factors, which require exquisite coordination of these crucial subcellular compartments. Moreover, regulation of replication events deals with a dynamic population of multiple mtDNA molecules rather than with a fixed number of genome copies, as it is the case for nuclear DNA. Importantly, the mechanistic aspects of mtDNA replication are still debated. We describe here major challenges for human mtDNA replication, the mechanistic aspects of the process that are to a large extent original, and their consequences on disease.

Polymorphic duplicate genes and persistent non-coding sequences reveal heterogeneous patterns of mitochondrial DNA loss in salamanders

Background

Mitochondria are the site of the citric acid cycle and oxidative phosphorylation (OXPHOS). In metazoans, the mitochondrial genome is a small, circular molecule averaging 16.5 kb in length. Despite evolutionarily conserved gene content, metazoan mitochondrial genomes show a diversity of gene orders most commonly explained by the duplication-random loss (DRL) model. In the DRL model, (1) a sequence of genes is duplicated in tandem, (2) one paralog sustains a loss-of-function mutation, resulting in selection to retain the other copy, and (3) the non-functional paralog is eventually deleted from the genome. Despite its apparent role in generating mitochondrial gene order diversity, little is known about the tempo and mode of random gene loss after duplication events. Here, we determine mitochondrial gene order across the salamander genus Aneides, which was previously shown to include at least two DRL-mediated rearrangement events. We then analyze these gene orders in a phylogenetic context to reveal patterns of DNA loss after mitochondrial gene duplication.

Results

We identified two separate duplication events that resulted in mitochondrial gene rearrangements in Aneides; one occurred at the base of the clade tens of millions of years ago, while the other occurred much more recently (i.e. within a single species), resulting in gene order polymorphism and paralogs that are readily identifiable. We demonstrate that near-complete removal of duplicate rRNA genes has occurred since the recent duplication event, whereas duplicate protein-coding genes persist as pseudogenes and duplicate tRNAs persist as functionally intact paralogs. In addition, we show that non-coding DNA duplicated at the base of the clade has persisted across species for tens of millions of years.

Conclusions

The evolutionary history of the mitochondrial genome, from its inception as a bacterial endosymbiont, includes massive genomic reduction. Consistent with this overall trend, selection for efficiency of mitochondrial replication and transcription has been hypothesized to favor elimination of extra sequence. Our results, however, suggest that there may be no strong disadvantage to extraneous sequences in salamander mitochondrial genomes, although duplicate rRNA genes may be deleterious.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4358-2) contains supplementary material, which is available to authorized users.

Phylogenetic association of Schizothorax plagiostomus with other schizothoracine fishes based on mitochondrial cytochrome B gene and control region

Cytochrome B (Cyt B) gene and control region of mitochondrial DNA are considered important for evaluating phylogenetic association of a species. In this study, we sequenced Cyt B and control region of Schizothorax plagiostomus and constructed phylogenetic association tree of S. plagiostomus with 23 schizothoracine fishes. We found S. plagiostomus to be closely associated with S. labiatus, S. richardsonii, S. progastus, and S. esocinus, with high-bootstraps values. Several conserved sequence blocks were identified within D-loop sequences. These are highly conserved within genus Schizothorax compared to other. This study reports the phylogenetic position of the S. plagiostomus among schizothoracines fishes and organization of D-loop region in S. plagiostomus from Pakistan.

Mitochondrial Genome Evolution and a Novel RNA Editing System in Deep-Branching Heteroloboseids

Discoba (Excavata) is an evolutionarily important group of eukaryotes that includes Jakobida, with the most bacterial-like mitochondrial genomes known, and Euglenozoa, many of which have extensively fragmented mitochondrial genomes. However, little is known about the mitochondrial genomes of Heterolobosea, the third main group of Discoba. Here, we studied two heteroloboseids—an undescribed amoeba “BB2” and Pharyngomonas kirbyi. Phylogenomic analysis revealed that they form a clade that is a sister group to all other Heterolobosea. We characterized the mitochondrial genomes of BB2 and P. kirbyi, which encoded 44 and 48 putative protein-coding genes respectively. Their gene contents were similar to that of Naegleria. In BB2, mitochondrially encoded RNAs were heavily edited, with ∼500 mononucleotide insertion events, mostly guanosines. These insertions always have the same identity as an adjacent nucleotide. Editing occurs in all ribosomal RNAs and protein-coding transcripts except one, and half of the transfer RNAs. Analysis of Illumina deep-sequencing data suggested that this RNA editing is very accurate and efficient, and most likely co-transcriptional. The dissimilarity of this editing process to other RNA editing phenomena in discobids, as well as its apparent absence in P. kirbyi, suggest that this remarkably extensive system of insertional editing evolved independently in the BB2 lineage, after its divergence from the P. kirbyi lineage.

Multiple reversals of strand asymmetry in molluscs mitochondrial genomes, and consequences for phylogenetic inferences

Strand asymmetry in nucleotide composition is a remarkable feature of animal mitochondrial genomes. The strand-specific bias in the nucleotide composition of the mtDNA has been known to be highly problematic for phylogenetic analyses. Here, the strand asymmetry was compared across 140 mollusc species and analyzed for a mtDNA fragment including twelve protein-coding genes. The analyses show that almost all species in Gastropoda (except Heterobranchia) and all species in Bivalvia present reversals of strand bias. The skew values on individual genes for all codon positions (P₁₂₃), third codon positions (P₃), and fourfold redundant third codon positions (P_4FD) indicated that CG skews are the best indicators of strand asymmetry. The differences in the patterns of strand asymmetry significantly influenced the amino acid composition of the encoded proteins. These biases are most striking for the amino acids Valine, Cysteine, Asparagine and Threonines, which appear to have evolved asymmetrical exchanges in response to shifts in nucleotide composition. Molluscs with strong variability of genome architectures (ARs) are usually characterized by a reversal of the usual strand bias. Phylogenetic analyses show that reversals of asymmetric mutational constraints have consequences on the phylogenetic inferences, as taxa characterized by reverse strand bias (Heterobranchia and Bivalvia) tend to group together due to long-branch attraction (LBA) artifacts. Neutral Transitions Excluded (NTE) model did not overcome the problem of heterogeneous biases present in molluscs mt genomes, suggested it may not be appropriate for molluscs mt genome data. Further refinement phylogenetic models may help us better understand internal relationships among these diverse organisms.

Lynn Margulis and the endosymbiont hypothesis: 50 years later

The 1967 article "On the Origin of Mitosing Cells" in the Journal of Theoretical Biology by Lynn Margulis (then Lynn Sagan) is widely regarded as stimulating renewed interest in the long-dormant endosymbiont hypothesis of organelle origins. In her article, not only did Margulis champion an endosymbiotic origin of mitochondria and plastids from bacterial ancestors, but she also posited that the eukaryotic flagellum (undulipodium in her usage) and mitotic apparatus originated from an endosymbiotic, spirochete-like organism. In essence, she presented a comprehensive symbiotic view of eukaryotic cell evolution (eukaryogenesis). Not all of the ideas in her article have been accepted, for want of compelling evidence, but her vigorous promotion of the role of symbiosis in cell evolution unquestionably had a major influence on how subsequent investigators have viewed the origin and evolution of mitochondria and plastids and the eukaryotic cell per se.

The complete mitochondrial genome of the edible Basidiomycete mushroom Phlebopus Portentosus

The complete mitochondrial genome of Phlebopus portentosus was determined using Illumina sequencing. The circular genome is 42,963 bp in length with GC content of 21.37%. It contains 14 putative protein-coding genes, the ribosomal RNA subunits and 23 tRNAs, all located on the same strand. The evolutionary relationships between P. portentosus and other 23 representative basidiomycete species were revealed based on sequences at the 14 concatenated mitochondrial protein-coding genes.

The mitochondrial genome of the plant-pathogenic fungus Stemphylium lycopersici uncovers a dynamic structure due to repetitive and mobile elements

Stemphylium lycopersici (Pleosporales) is a plant-pathogenic fungus that has been associated with a broad range of plant-hosts worldwide. It is one of the causative agents of gray leaf spot disease in tomato and pepper. The aim of this work was to characterize the mitochondrial genome of S. lycopersici CIDEFI-216, to use it to trace taxonomic relationships with other fungal taxa and to get insights into the evolutionary history of this phytopathogen. The complete mitochondrial genome was assembled into a circular double-stranded DNA molecule of 75,911 bp that harbors a set of 37 protein-coding genes, 2 rRNA genes (rns and rnl) and 28 tRNA genes, which are transcribed from both sense and antisense strands. Remarkably, its gene repertoire lacks both atp8 and atp9, contains a free-standing gene for the ribosomal protein S3 (rps3) and includes 13 genes with homing endonuclease domains that are mostly located within its 15 group I introns. Strikingly, subunits 1 and 2 of cytochrome oxidase are encoded by a single continuous open reading frame (ORF). A comparative mitogenomic analysis revealed the large extent of structural rearrangements among representatives of Pleosporales, showing the plasticity of their mitochondrial genomes. Finally, an exhaustive phylogenetic analysis of the subphylum Pezizomycotina based on mitochondrial data reconstructed their relationships in concordance with several studies based on nuclear data. This is the first report of a mitochondrial genome belonging to a representative of the family Pleosporaceae.

The pentatricopeptide repeat protein MTSF2 stabilizes a nad1 precursor transcript and defines the 3΄ end of its 5΄-half intron

RNA expression in plant mitochondria implies a large number of post-transcriptional events in which transcript processing and stabilization are essential. In this study, we analyzed the function of the Arabidopsis mitochondrial stability factor 2 gene (MTSF2) and show that the encoded pentatricopeptide repeat protein is essential for the accumulation of stable nad1 mRNA. The production of mature nad1 requires the assembly of three independent RNA precursors via two trans-splicing reactions. Genetic analyses revealed that the lack of nad1 in mtsf2 mutants results from the specific destabilization of the nad1 exons 2-3 precursor transcript. We further demonstrated that MTSF2 binds to its 3΄ extremity with high affinity, suggesting a protective action by blocking exoribonuclease progression. By defining the 3΄ end of nad1 exons 2-3 precursor, MTSF2 concomitantly determines the 3΄ extremity of the first half of the trans-intron found at the end of the transcript. Therefore, binding of the MTSF2 protein to nad1 exons 2-3 precursor evolved both to stabilize the transcript and to define a 3΄ extremity compatible with the trans-splicing reaction needed to reconstitute mature nad1. We thus reveal that the range of transcripts stabilized by association with protective protein on their 3΄ end concerns also mitochondrial precursor transcripts.

Mitigating Mitochondrial Genome Erosion Without Recombination

Mitochondria are ATP-producing organelles of bacterial ancestry that played a key role in the origin and early evolution of complex eukaryotic cells. Most modern eukaryotes transmit mitochondrial genes uniparentally, often without recombination among genetically divergent organelles. While this asymmetric inheritance maintains the efficacy of purifying selection at the level of the cell, the absence of recombination could also make the genome susceptible to Muller's ratchet. How mitochondria escape this irreversible defect accumulation is a fundamental unsolved question. Occasional paternal leakage could in principle promote recombination, but it would also compromise the purifying selection benefits of uniparental inheritance. We assess this tradeoff using a stochastic population-genetic model. In the absence of recombination, uniparental inheritance of freely-segregating genomes mitigates mutational erosion, while paternal leakage exacerbates the ratchet effect. Mitochondrial fusion-fission cycles ensure independent genome segregation, improving purifying selection. Paternal leakage provides opportunity for recombination to slow down the mutation accumulation, but always at a cost of increased steady-state mutation load. Our findings indicate that random segregation of mitochondrial genomes under uniparental inheritance can effectively combat the mutational meltdown, and that homologous recombination under paternal leakage might not be needed.

Quantification of mitochondrial DNA copy number in suspected cancer patients by a well optimized ddPCR method

Changes in mitochondrial DNA (mtDNA) content is a useful clinical biomarker for various diseases, however results are controversial as several analytical factors can affect measurement of mtDNA. MtDNA is often quantified by taking ratio between a target mitochondrial gene and a reference nuclear gene (mtDNA/nDNA) using quantitative real time PCR often on two separate experiments. It measures relative levels by using external calibrator which may not be comparable across laboratories. We have developed and optimized a droplet digital PCR (ddPCR) based method for quantification of absolute copy number of both mtDNA and nDNA gene in whole blood. Finally, the role of mtDNA in suspected cancer patients referred to a cancer diagnostic center was investigated.

Analytical factors which can result in false quantification of mtDNA have been optimized and both target and reference have been quantified simultaneously with intra- and inter-assay coefficient variances as 3.1% and 4.2% respectively. Quantification of mtDNA show that compared to controls, solid tumors (but not hematologic malignancies) and other diseases had significantly lower copy number of mtDNA. Higher mtDNA (highest quartile) was associated with a significantly lower risk of both solid tumors and other diseases, independent of age and sex. Receiver operating curve demonstrated that mtDNA levels could differentiate controls from patients with solid tumors and other diseases.

Quantification of mtDNA by a well optimized ddPCR method showed that its depletion may be a hallmark of general illness and can be used to stratify healthy individuals from patients diagnosed with cancer and other chronic diseases.

Biological action in Read-Write genome evolution

Many of the most important evolutionary variations that generated phenotypic adaptations and originated novel taxa resulted from complex cellular activities affecting genome content and expression. These activities included (i) the symbiogenetic cell merger that produced the mitochondrion-bearing ancestor of all extant eukaryotes, (ii) symbiogenetic cell mergers that produced chloroplast-bearing ancestors of photosynthetic eukaryotes, and (iii) interspecific hybridizations and genome doublings that generated new species and adaptive radiations of higher plants and animals. Adaptive variations also involved horizontal DNA transfers and natural genetic engineering by mobile DNA elements to rewire regulatory networks, such as those essential to viviparous reproduction in mammals. In the most highly evolved multicellular organisms, biological complexity scales with 'non-coding' DNA content rather than with protein-coding capacity in the genome. Coincidentally, 'non-coding' RNAs rich in repetitive mobile DNA sequences function as key regulators of complex adaptive phenotypes, such as stem cell pluripotency. The intersections of cell fusion activities, horizontal DNA transfers and natural genetic engineering of Read-Write genomes provide a rich molecular and biological foundation for understanding how ecological disruptions can stimulate productive, often abrupt, evolutionary transformations.

Heteroplasmy due to coexistence of mtCOI haplotypes from different lineages of the Thrips tabaci cryptic species group

Heteroplasmy is the existence of multiple mitochondrial DNA haplotypes within the cell. Although the number of reports of heteroplasmy is increasing for arthropods, the occurrence, number of variants, and origins are not well studied. In this research, the occurrence of heteroplasmy was investigated in Thrips tabaci, a putative species complex whose lineages can be distinguished by their mitochondrial DNA haplotypes. The results from this study showed that heteroplasmy was due to the occurrence of mitochondrial cytochrome oxydase I (mtCOI) haplotypes from two different T. tabaci lineages. An assay using flow cytometry and quantitative real-time PCR was then used to quantify the per cell copy number of the two mtCOI haplotypes present in individuals exhibiting heteroplasmy from nine geographically distant populations in India. All of the T. tabaci individuals in this study were found to exhibit heteroplasmy, and in every individual the per cell copy number of mtCOI from lineage 3 comprised 75-98% of the haplotypes detected and was variable among individuals tested. There was no evidence to suggest that the presense of lineage-specific haplotypes was due to nuclear introgression; however, further studies are needed to investigate nuclear introgression and paternal leakage during rare interbreeding between individuals from lineages 2 and 3.

SMRT Sequencing Revealed Mitogenome Characteristics and Mitogenome-Wide DNA Modification Pattern in Ophiocordyceps sinensis

Single molecule, real-time (SMRT) sequencing was used to characterize mitochondrial (mt) genome of Ophiocordyceps sinensis and to analyze the mt genome-wide pattern of epigenetic DNA modification. The complete mt genome of O. sinensis, with a size of 157,539 bp, is the fourth largest Ascomycota mt genome sequenced to date. It contained 14 conserved protein-coding genes (PCGs), 1 intronic protein rps3, 27 tRNAs and 2 rRNA subunits, which are common characteristics of the known mt genomes in Hypocreales. A phylogenetic tree inferred from 14 PCGs in Pezizomycotina fungi supports O. sinensis as most closely related to Hirsutella rhossiliensis in Ophiocordycipitaceae. A total of 36 sequence sites in rps3 were under positive selection, with dN/dS >1 in the 20 compared fungi. Among them, 16 sites were statistically significant. In addition, the mt genome-wide base modification pattern of O. sinensis was determined in this study, especially DNA methylation. The methylations were located in coding and uncoding regions of mt PCGs in O. sinensis, and might be closely related to the expression of PCGs or the binding affinity of transcription factor A to mtDNA. Consequently, these methylations may affect the enzymatic activity of oxidative phosphorylation and then the mt respiratory rate; or they may influence mt biogenesis. Therefore, methylations in the mitogenome of O. sinensis might be a genetic feature to adapt to the cold and low PO2 environment at high altitude, where O. sinensis is endemic. This is the first report on epigenetic modifications in a fungal mt genome.

Complete mitochondrial genomes of prasinophyte algae Pyramimonas parkeae and Cymbomonas tetramitiformis

Mitochondria are archetypal eukaryotic organelles that were acquired by endosymbiosis of an ancient species of alpha-proteobacteria by the last eukaryotic common ancestor. The genetic information contained within the mitochondrial genome has been an important source of information for resolving relationships among eukaryotic taxa. In this study, we utilized mitochondrial and chloroplast genomes to explore relationships among prasinophytes. Prasinophytes are represented by diverse early-diverging green algae whose physical structures and genomes have the potential to elucidate the traits of the last common ancestor of the Viridiplantae (or Chloroplastida). We constructed de novo mitochondrial genomes for two prasinophyte algal species, Pyramimonas parkeae and Cymbomonas tetramitiformis, representing the prasinophyte clade. Comparisons of genome structure and gene order between these species and to those of other prasinophytes revealed that the mitochondrial genomes of P. parkeae and C. tetramitiformis are more similar to each other than to other prasinophytes, consistent with other molecular inferences of the close relationship between these two species. Phylogenetic analyses using the inferred amino acid sequences of mitochondrial and chloroplast protein-coding genes resolved a clade consisting of P. parkeae and C. tetramitiformis; and this group (representing the prasinophyte clade I) branched with the clade II, consistent with previous studies based on the use of nuclear gene markers.

Unparalleled replacement of native mitochondrial genes by foreign homologs in a holoparasitic plant

Horizontal gene transfer (HGT) among flowering plant mitochondria occurs frequently and, in most cases, leads to nonfunctional transgenes in the recipient genome. Parasitic plants are particularly prone to this phenomenon, but their mitochondrial genomes (mtDNA) have been largely unexplored. We undertook a large-scale mitochondrial genomic study of the holoparasitic plant Lophophytum mirabile (Balanophoraceae). Comprehensive phylogenetic analyses were performed to address the frequency, origin, and impact of HGT. The sequencing of the complete mtDNA of L. mirabile revealed the unprecedented acquisition of host-derived mitochondrial genes, representing 80% of the protein-coding gene content. All but two of these foreign genes replaced the native homologs and are probably functional in energy metabolism. The genome consists of 54 circular-mapping chromosomes, 25 of which carry no intact genes. The likely functional replacement of up to 26 genes in L. mirabile represents a stunning example of the potential effect of rampant HGT on plant mitochondria. The use of host-derived genes may have a positive effect on the host-parasite relationship, but could also be the result of nonadaptive forces.

Dek35 Encodes a PPR Protein that Affects cis-Splicing of Mitochondrial nad4 Intron 1 and Seed Development in Maize

In higher plants, the splicing of organelle-encoded mRNA involves a complex collaboration with nuclear-encoded proteins. Pentatricopeptide repeat (PPR) proteins have been implicated in these RNA-protein interactions. In this study, we performed the cloning and functional characterization of maize Defective kernel 35 (Dek35). The dek35-ref mutant is a lethal-seed mutant with developmental deficiency. Dek35 was cloned through Mutator tag isolation and further confirmed by four additional independent mutant alleles. Dek35 encodes an P-type PPR protein that targets the mitochondria. The dek35 mutation causes significant reduction in the accumulation of DEK35 proteins and reduced splicing efficiency of mitochondrial nad4 intron 1. Analysis of mitochondrial complex in dek35 immature seeds indicated severe deficiency in the complex I assembly and NADH dehydrogenase activity. Transcriptome analysis of dek35 endosperm revealed enhanced expression of genes involved in the alternative respiratory pathway and extensive differentially expressed genes related to mitochondrial function and activity. Collectively, these results indicate that Dek35 encodes an PPR protein that affects the cis-splicing of mitochondrial nad4 intron 1 and is required for mitochondrial function and seed development.

The first complete mitochondrial genomes of subterranean dytiscid diving beetles (Limbodessus and Paroster) from calcrete aquifers of Western Australia

Comparative analyses of mitochondrial (mt) genomes may provide insights into the genetic changes, associated with metabolism, that occur when surface species adapt to living in underground habitats. Such analyses require comparisons among multiple independently evolved subterranean species, with the dytiscid beetle fauna from the calcrete archipelago of central Western Australia providing an outstanding model system to do this. Here, we present the first whole mt genomes from four subterranean dytiscid beetle species of the genera Limbodessus (L. palmulaoides) and Paroster (P. macrosturtensis, P. mesosturtensis and P. microsturtensis) and compare genome sequences with those from surface dytiscid species. The mt genomes were sequenced using a next-generation sequencing approach employing the Illumina Miseq system and assembled de novo. All four mt genomes are circular, ranging in size from 16 504 to 16 868 bp, and encode 37 genes and a control region. The overall structure (gene number, orientation and order) of the mt genomes is the same as that found in eight sequenced surface species, but with genome size variation resulting from length variation of intergenic regions and the control region . Our results provide a basis for future investigations of adaptive evolutionary changes that may occur in mt genes when species move underground.

The Role of Mitochondria in AMD: Current Knowledge and Future Applications

Mitochondria are organelles which comprise the main respiratory machinery in the eukaryotic cells. In addition to their crucial role in energy production, they have profound effects on apoptosis and retrograde signaling to nucleus. Mitochondria have their own DNA, which codes for different proteins mostly involved in oxidative phosphorylation. Significant changes in the mitochondria of retinal pigment epithelium have been reported in age-related macular degeneration (AMD), which is correlated with the severity of the disease. Cybrid cell lines that have identical nuclei but mitochondria from different individuals can provide a unique means for studying the relationship between mitochondria and AMD. Different approaches for protection of mitochondria have been introduced which can be considered as potential future treatments for AMD and other age- related disorders.

Great migration: epigenetic reprogramming and germ cell-oocyte metamorphosis determine individual ovarian reserve

Emigration is defined as a synchronized movement of germ cells between the yolk sack and genital ridges. The miraculous migration of germ cells resembles the remigration of salmon traveling from one habitat to other. This migration of germ cells is indispensible for the development of new generations. It is not, however, clear why germ cells differentiate during migration but not at the place of origin. In order to escape harmful somatic signals which might disturb the proper establishment of germ cells forced germ cell migration may be necessary. Another reason may be to benefit from the opportunities of new habitats. Therefore, emigration may have powerful effects on the population dynamics of the immigrant germ cells. While some of these cells do reach their target, some others die or reach to wrong targets. Only germ cell precursors with genetically, and structurally powerful can reach their target. Likewise, epigenetic reprogramming in both migratory and post-migratory germ cells is essential for the establishment of totipotency. During this journey some germ cells may sacrifice themselves for the goodness of the others. The number and quality of germ cells reaching the genital ridge may vary depending on the problems encountered during migration. If the aim in germ cell specification is to provide an optimal ovarian reserve for the continuity of the generation, then this cascade of events cannot be only accomplished at the same level for every one but also are manifested by several outcomes. This is significant evidence supporting the possibility of unique individual ovarian reserve.

Monogonont Rotifer, Brachionus calyciflorus, Possesses Exceptionally Large, Fragmented Mitogenome

In contrast to the highly conserved mitogenomic structure and organisation in most animals (including rotifers), the two previously sequenced monogonont rotifer mitogenomes were fragmented into two chromosomes similar in size, each of which possessed one major non-coding region (mNCR) of about 4–5 Kbp. To further explore this phenomenon, we have sequenced and analysed the mitogenome of one of the most studied monogonont rotifers, Brachionus calyciflorus. It is also composed of two circular chromosomes, but the chromosome-I is extremely large (27 535 bp; 3 mNCRs), whereas the chromosome-II is relatively small (9 833 bp; 1 mNCR). With the total size of 37 368 bp, it is one of the largest metazoan mitogenomes ever reported. In comparison to other monogononts, gene distribution between the two chromosomes and gene order are different and the number of mNCRs is doubled. Atp8 was not found (common in rotifers), and Cytb was present in two copies (the first report in rotifers). A high number (99) of SNPs indicates fast evolution of the Cytb-1 copy. The four mNCRs (5.3–5.5 Kb) were relatively similar. Publication of this sequence shall contribute to the understanding of the evolutionary history of the unique mitogenomic organisation in this group of rotifers.

Mitochondrial-nuclear co-evolution leads to hybrid incompatibility through pentatricopeptide repeat proteins

Mitochondrial-nuclear incompatibility has a major role in reproductive isolation between species. However, the underlying mechanism and driving force of mitochondrial-nuclear incompatibility remain elusive. Here, we report a pentatricopeptide repeat-containing (PPR) protein, Ccm1, and its interacting partner, 15S rRNA, to be involved in hybrid incompatibility between two yeast species, Saccharomyces cerevisiae and Saccharomyces bayanus S. bayanus-Ccm1 has reduced binding affinity for S. cerevisiae-15S rRNA, leading to respiratory defects in hybrid cells. This incompatibility can be rescued by single mutations on several individual PPR motifs, demonstrating the highly evolvable nature of PPR proteins. When we examined other PPR proteins in the closely related Saccharomyces sensu stricto yeasts, about two-thirds of them showed detectable incompatibility. Our results suggest that fast co-evolution between flexible PPR proteins and their mitochondrial RNA substrates may be a common driving force in the development of mitochondrial-nuclear hybrid incompatibility.

The Complete Mitochondrial Genome of Aleurocanthus camelliae: Insights into Gene Arrangement and Genome Organization within the Family Aleyrodidae

There are numerous gene rearrangements and transfer RNA gene absences existing in mitochondrial (mt) genomes of Aleyrodidae species. To understand how mt genomes evolved in the family Aleyrodidae, we have sequenced the complete mt genome of Aleurocanthus camelliae and comparatively analyzed all reported whitefly mt genomes. The mt genome of A. camelliae is 15,188 bp long, and consists of 13 protein-coding genes, two rRNA genes, 21 tRNA genes and a putative control region (GenBank: KU761949). The tRNA gene, trnI, has not been observed in this genome. The mt genome has a unique gene order and shares most gene boundaries with Tetraleurodes acaciae. Nineteen of 21 tRNA genes have the conventional cloverleaf shaped secondary structure and two (trnS₁ and trnS₂) lack the dihydrouridine (DHU) arm. Using ARWEN and homologous sequence alignment, we have identified five tRNA genes and revised the annotation for three whitefly mt genomes. This result suggests that most absent genes exist in the genomes and have not been identified, due to be lack of technology and inference sequence. The phylogenetic relationships among 11 whiteflies and Drosophila melanogaster were inferred by maximum likelihood and Bayesian inference methods. Aleurocanthus camelliae and T. acaciae form a sister group, and all three Bemisia tabaci and two Bemisia afer strains gather together. These results are identical to the relationships inferred from gene order. We inferred that gene rearrangement plays an important role in the mt genome evolved from whiteflies.

Moramonas marocensis gen. nov., sp. nov.: a jakobid flagellate isolated from desert soil with a bacteria-like, but bloated mitochondrial genome

A new jakobid genus has been isolated from Moroccan desert soil. The cyst-forming protist Moramonas marocensis gen. nov., sp. nov. has two anteriorly inserted flagella of which one points to the posterior cell pole accompanying the ventral feeding groove and is equipped with a dorsal vane-a feature typical for the Jakobida. It further shows a flagellar root system consisting of singlet microtubular root, left root (R1), right root (R2) and typical fibres associated with R1 and R2. The affiliation of M. marocensis to the Jakobida was confirmed by molecular phylogenetic analyses of the SSU rRNA gene, five nuclear genes and 66 mitochondrial protein-coding genes. The mitochondrial genome has the high number of genes typical for jakobids, and bacterial features, such as the four-subunit RNA polymerase and Shine-Dalgarno sequences upstream of the coding regions of several genes. The M. marocensis mitochondrial genome encodes a similar number of genes as other jakobids, but is unique in its very large genome size (greater than 264 kbp), which is three to four times higher than that of any other jakobid species investigated yet. This increase seems to be due to a massive expansion in non-coding DNA, creating a bloated genome like those of plant mitochondria.

Mitochondrial Function and Maize Kernel Development Requires Dek2, a Pentatricopeptide Repeat Protein Involved in nad1 mRNA Splicing

In flowering plants, many respiration-related proteins are encoded by the mitochondrial genome and the splicing of mitochondrion-encoded messenger RNA (mRNA) involves a complex collaboration with nuclear-encoded proteins. Pentatricopeptide repeat (PPR) proteins have been implicated in these RNA-protein interactions. Maize defective kernel 2 (dek2) is a classic mutant with small kernels and delayed development. Through positional cloning and allelic confirmation, we found Dek2 encodes a novel P-type PPR protein that targets mitochondria. Mitochondrial transcript analysis indicated that dek2 mutation causes reduced splicing efficiency of mitochondrial nad1 intron 1. Mitochondrial complex analysis in dek2 immature kernels showed severe deficiency of complex I assembly. Dramatically up-regulated expression of alternative oxidases (AOXs), transcriptome data, and TEM analysis results revealed that proper splicing of nad1 is critical for mitochondrial functions and inner cristaes morphology. This study indicated that Dek2 is a new PPR protein that affects the splicing of mitochondrial nad1 intron 1 and is required for mitochondrial function and kernel development.

Mitochondrial translation factors reflect coordination between organelles and cytoplasmic translation via mTOR signaling: Implication in disease

Mitochondria are semi-autonomous organelle possessing their own translation machinery to biosynthesize mitochondrial DNA (mtDNA)-encoded polypeptides, which are the core subunits of oxidative phosphorylation (OXPHOS) complexes. Mitochondrial translation elongation factor 4 (mtEF4) is a key quality control factor in mitochondrial translation (mt-translation) that regulates mitochondrial tRNA translocation and modulates cellular responses by influencing cytoplasmic translation (ct-translation). In addition to mtEF4, mt-translational activators, mitochondrial microRNAs (mitomiRs), and MITRAC have been reported recently as crucial mt-translation regulators. Here, we focus on the novel ways how these factors regulate mt-translation, discuss the main cellular response of mammalian target of rapamycin (mTOR) signalling upon mt-translation defects, and summarize the related human diseases.

Analysis of mitochondrial genetic diversity of Ustilago maydis in Mexico

The current understanding of the genetic diversity of the phytopathogenic fungus Ustilago maydis is limited. To determine the genetic diversity and structure of U. maydis, 48 fungal isolates were analyzed using mitochondrial simple sequence repeats (SSRs). Tumours (corn smut or 'huitlacoche') were collected from different Mexican states with diverse environmental conditions. Using bioinformatic tools, five microsatellites were identified within intergenic regions of the U. maydis mitochondrial genome. SSRMUM4 was the most polymorphic marker. The most common repeats were hexanucleotides. A total of 12 allelic variants were identified, with a mean of 2.4 alleles per locus. An estimate of the genetic diversity using analysis of molecular variance (AMOVA) revealed that the highest variance component is within states (84%), with moderate genetic differentiation between states (16%) (F_ST= 0.158). A dendrogram generated using the unweighted paired-grouping method with arithmetic averages (UPGMA) and the Bayesian analysis of population structure grouped the U. maydis isolates into two subgroups (K = 2) based on their shared SSRs.

Complete mitochondrial genomes of Trisidos kiyoni and Potiarca pilula: Varied mitochondrial genome size and highly rearranged gene order in Arcidae

We present the complete mitochondrial genomes (mitogenomes) of Trisidos kiyoni and Potiarca pilula, both important species from the family Arcidae (Arcoida: Arcacea). Typical bivalve mtDNA features were described, such as the relatively conserved gene number (36 and 37), a high A + T content (62.73% and 61.16%), the preference for A + T-rich codons, and the evidence of non-optimal codon usage. The mitogenomes of Arcidae species are exceptional for their extraordinarily large and variable sizes and substantial gene rearrangements. The mitogenome of T. kiyoni (19,614 bp) and P. pilula (28,470 bp) are the two smallest Arcidae mitogenomes. The compact mitogenomes are weakly associated with gene number and primarily reflect shrinkage of the non-coding regions. The varied size in Arcidae mitogenomes reflect a dynamic history of expansion. A significant positive correlation is observed between mitogenome size and the combined length of cox1-3, the lengths of Cytb, and the combined length of rRNAs (rrnS and rrnL) (P < 0.001). Both protein coding genes (PCGs) and tRNA rearrangements is observed in P. pilula and T. kiyoni mitogenomes. This analysis imply that the complicated gene rearrangement in mitochondrial genome could be considered as one of key characters in inferring higher-level phylogenetic relationship of Arcidae.

Genes in Hiding

Unrecognizable genes are an unsettling problem in genomics. Here, we survey the various types of cryptic genes and the corresponding deciphering strategies employed by cells. Encryption that renders genes substantially different from homologs in other species includes sequence substitution, insertion, deletion, fragmentation plus scrambling, and invasion by mobile genetic elements. Cells decode cryptic genes at the DNA, RNA or protein level. We will focus on a recently discovered case of unparalleled encryption involving massive gene fragmentation and nucleotide deletions and substitutions, occurring in the mitochondrial genome of a poorly understood protist group, the diplonemids. This example illustrates that comprehensive gene detection requires not only auxiliary sequence information - transcriptome and proteome data - but also knowledge about a cell's deciphering arsenal.

Complete Chloroplast and Mitochondrial Genome Sequences of the Hydrocarbon Oil-Producing Green Microalga Botryococcus braunii Race B (Showa)

The green alga Botryococcus braunii is capable of the production and excretion of high quantities of long-chain hydrocarbons and exopolysaccharides. In this study, we present the complete plastid and mitochondrial genomes of the hydrocarbon-producing microalga Botryococcus braunii race B (Showa), with a total length of 156,498 and 129,356 bp, respectively.

Nothing in Evolution Makes Sense Except in the Light of Genomics: Read-Write Genome Evolution as an Active Biological Process

The 21st century genomics-based analysis of evolutionary variation reveals a number of novel features impossible to predict when Dobzhansky and other evolutionary biologists formulated the neo-Darwinian Modern Synthesis in the middle of the last century. These include three distinct realms of cell evolution; symbiogenetic fusions forming eukaryotic cells with multiple genome compartments; horizontal organelle, virus and DNA transfers; functional organization of proteins as systems of interacting domains subject to rapid evolution by exon shuffling and exonization; distributed genome networks integrated by mobile repetitive regulatory signals; and regulation of multicellular development by non-coding lncRNAs containing repetitive sequence components. Rather than single gene traits, all phenotypes involve coordinated activity by multiple interacting cell molecules. Genomes contain abundant and functional repetitive components in addition to the unique coding sequences envisaged in the early days of molecular biology. Combinatorial coding, plus the biochemical abilities cells possess to rearrange DNA molecules, constitute a powerful toolbox for adaptive genome rewriting. That is, cells possess "Read-Write Genomes" they alter by numerous biochemical processes capable of rapidly restructuring cellular DNA molecules. Rather than viewing genome evolution as a series of accidental modifications, we can now study it as a complex biological process of active self-modification.

Intron Derived Size Polymorphism in the Mitochondrial Genomes of Closely Related Chrysoporthe Species

In this study, the complete mitochondrial (mt) genomes of Chrysoporthe austroafricana (190,834 bp), C. cubensis (89,084 bp) and C. deuterocubensis (124,412 bp) were determined. Additionally, the mitochondrial genome of another member of the Cryphonectriaceae, namely Cryphonectria parasitica (158,902 bp), was retrieved and annotated for comparative purposes. These genomes showed high levels of synteny, especially in regions including genes involved in oxidative phosphorylation and electron transfer, unique open reading frames (uORFs), ribosomal RNAs (rRNAs) and transfer RNAs (tRNAs), as well as intron positions. Comparative analyses revealed signatures of duplication events, intron number and length variation, and varying intronic ORFs which highlighted the genetic diversity of mt genomes among the Cryphonectriaceae. These mt genomes showed remarkable size polymorphism. The size polymorphism in the mt genomes of these closely related Chrysoporthe species was attributed to the varying number and length of introns, coding sequences and to a lesser extent, intergenic sequences. Compared to publicly available fungal mt genomes, the C. austroafricana mt genome is the second largest in the Ascomycetes thus far.

From microbiology to cell biology: when an intracellular bacterium becomes part of its host cell

Mitochondria and chloroplasts are now called organelles, but they used to be bacteria. As they transitioned from endosymbionts to organelles, they became more and more integrated into the biochemistry and cell biology of their hosts. Work over the last 15 years has shown that other symbioses show striking similarities to mitochondria and chloroplasts. In particular, many sap-feeding insects house intracellular bacteria that have genomes that overlap mitochondria and chloroplasts in terms of size and coding capacity. The massive levels of gene loss in some of these bacteria suggest that they, too, are becoming highly integrated with their host cells. Understanding these bacteria will require inspiration from eukaryotic cell biology, because a traditional microbiological framework is insufficient for understanding how they work.

Distinct role of Arabidopsis mitochondrial P-type pentatricopeptide repeat protein-modulating editing protein, PPME, in nad1 RNA editing

The mitochondrion is an important power generator in most eukaryotic cells. To preserve its function, many essential nuclear-encoded factors play specific roles in mitochondrial RNA metabolic processes, including RNA editing. RNA editing consists of post-transcriptional deamination, which alters specific nucleotides in transcripts to mediate gene expression. In plant cells, many pentatricopeptide repeat proteins (PPRs) participate in diverse organellar RNA metabolic processes, but only PLS-type PPRs are involved in RNA editing. Here, we report a P-type PPR protein from Arabidopsis thaliana, P-type PPR-Modulating Editing (PPME), which has a distinct role in mitochondrial nad1 RNA editing via RNA binding activity. In the homozygous ppme mutant, cytosine (C)-to-uracil (U) conversions at both the nad1-898 and 937 sites were abolished, disrupting Arg³⁰⁰-to-Trp³⁰⁰ and Pro³¹³-to-Ser³¹³ amino acid changes in the mitochondrial NAD1 protein. NAD1 is a critical component of mitochondrial respiration complex I; its activity is severely reduced in the homozygous ppme mutant, resulting in significantly altered growth and development. Both abolished RNA editing and defective complex I activity were completely rescued by CaMV 35S promoter- and PPME native promoter-driven PPME genomic fragments tagged with GFP in a homozygous ppme background. Our experimental results demonstrate a distinct role of a P-type PPR protein, PPME, in RNA editing in plant organelles.

Towards understanding the evolution and functional diversification of DNA-containing plant organelles

Plastids and mitochondria derive from prokaryotic symbionts that lost most of their genes after the establishment of endosymbiosis. In consequence, relatively few of the thousands of different proteins in these organelles are actually encoded there. Most are now specified by nuclear genes. The most direct way to reconstruct the evolutionary history of plastids and mitochondria is to sequence and analyze their relatively small genomes. However, understanding the functional diversification of these organelles requires the identification of their complete protein repertoires – which is the ultimate goal of organellar proteomics. In the meantime, judicious combination of proteomics-based data with analyses of nuclear genes that include interspecies comparisons and/or predictions of subcellular location is the method of choice. Such genome-wide approaches can now make use of the entire sequences of plant nuclear genomes that have emerged since 2000. Here I review the results of these attempts to reconstruct the evolution and functions of plant DNA-containing organelles, focusing in particular on data from nuclear genomes. In addition, I discuss proteomic approaches to the direct identification of organellar proteins and briefly refer to ongoing research on non-coding nuclear DNAs of organellar origin (specifically, nuclear mitochondrial DNA and nuclear plastid DNA).