Mitochondrial genomes: anything goes

Recent advances in understanding mitochondrial genome diversity

Ever since its discovery, the double-stranded DNA contained in the mitochondria of eukaryotes has fascinated researchers because of its bacterial endosymbiotic origin, crucial role in encoding subunits of the respiratory complexes, compact nature, and specific inheritance mechanisms. In the last few years, high-throughput sequencing techniques have accelerated the sequencing of mitochondrial genomes (mitogenomes) and uncovered the great diversity of organizations, gene contents, and modes of replication and transcription found in living eukaryotes. Some early divergent lineages of unicellular eukaryotes retain certain synteny and gene content resembling those observed in the genomes of alphaproteobacteria (the inferred closest living group of mitochondria), whereas others adapted to anaerobic environments have drastically reduced or even lost the mitogenome. In the three main multicellular lineages of eukaryotes, mitogenomes have pursued diverse evolutionary trajectories in which different types of molecules (circular versus linear and single versus multipartite), gene structures (with or without self-splicing introns), gene contents, gene orders, genetic codes, and transfer RNA editing mechanisms have been selected. Whereas animals have evolved a rather compact mitochondrial genome between 11 and 50 Kb in length with a highly conserved gene content in bilaterians, plants exhibit large mitochondrial genomes of 66 Kb to 11.3 Mb with large intergenic repetitions prone to recombination, and fungal mitogenomes have intermediate sizes of 12 to 236 Kb.

The draft mitochondrial genome of Magnolia biondii and mitochondrial phylogenomics of angiosperms

The mitochondrial genomes of flowering plants are well known for their large size, variable coding-gene set and fluid genome structure. The available mitochondrial genomes of the early angiosperms show extreme genetic diversity in genome size, structure, and sequences, such as rampant HGTs in Amborella mt genome, numerous repeated sequences in Nymphaea mt genome, and conserved gene evolution in Liriodendron mt genome. However, currently available early angiosperm mt genomes are still limited, hampering us from obtaining an overall picture of the mitogenomic evolution in angiosperms. Here we sequenced and assembled the draft mitochondrial genome of Magnolia biondii Pamp. from Magnoliaceae (magnoliids) using Oxford Nanopore sequencing technology. We recovered a single linear mitochondrial contig of 967,100 bp with an average read coverage of 122 × and a GC content of 46.6%. This draft mitochondrial genome contains a rich 64-gene set, similar to those of Liriodendron and Nymphaea, including 41 protein-coding genes, 20 tRNAs, and 3 rRNAs. Twenty cis-spliced and five trans-spliced introns break ten protein-coding genes in the Magnolia mt genome. Repeated sequences account for 27% of the draft genome, with 17 out of the 1,145 repeats showing recombination evidence. Although partially assembled, the approximately 1-Mb mt genome of Magnolia is still among the largest in angiosperms, which is possibly due to the expansion of repeated sequences, retention of ancestral mtDNAs, and the incorporation of nuclear genome sequences. Mitochondrial phylogenomic analysis of the concatenated datasets of 38 conserved protein-coding genes from 91 representatives of angiosperm species supports the sister relationship of magnoliids with monocots and eudicots, which is congruent with plastid evidence.

Evidence of Extensive Intraspecific Noncoding Reshuffling in a 169-kb Mitochondrial Genome of a Basidiomycetous Fungus

Comparative genomics of fungal mitochondrial genomes (mitogenomes) have revealed a remarkable pattern of rearrangement between and within major phyla owing to horizontal gene transfer and recombination. The role of recombination was exemplified at a finer evolutionary time scale in basidiomycetes group of fungi as they display a diversity of mitochondrial DNA inheritance patterns. Here, we assembled mitogenomes of six species from the Hymenochaetales order of basidiomycetes and examined 59 mitogenomes from 2 genetic lineages of Phellinus noxius. Gene order is largely collinear, while intergene regions are major determinants of mitogenome size variation. Substantial sequence divergence was found in shared introns consistent with high horizontal gene transfer frequency observed in yeasts, but we also identified a rare case where an intron was retained in five species since speciation. In contrast to the hyperdiversity observed in nuclear genomes of Phellinus noxius, mitogenomes’ intraspecific polymorphisms at protein-coding sequences are extremely low. Phylogeny network based on introns revealed turnover as well as exchange of introns between two lineages. Strikingly, some strains harbor a mosaic origin of introns from both lineages. Analysis of intergenic sequence indicated substantial differences between and within lineages, and an expansion may be ongoing as a result of exchange between distal intergenes. These findings suggest that the evolution in mitochondrial DNAs is usually lineage specific but chimeric mitotypes are frequently observed, thus capturing the possible evolutionary processes shaping mitogenomes in a basidiomycete. The large mitogenome sizes reported in various basidiomycetes appear to be a result of interspecific reshuffling of intergenes.

Phylogenetic relationship between Australian Fusarium oxysporum isolates and resolving the species complex using the multispecies coalescent model

BACKGROUND

The Fusarium oxysporum species complex (FOSC) is a ubiquitous group of fungal species readily isolated from agroecosystem and natural ecosystem soils which includes important plant and human pathogens. Genetic relatedness within the complex has been studied by sequencing either the genes or the barcoding gene regions within those genes. Phylogenetic analyses have demonstrated a great deal of diversity which is reflected in the differing number of clades identified: three, five and eight. Genetic limitation within the species in the complex has been studied through Genealogical Concordance Phylogenetic Species Recognition (GCPSR) analyses with varying number of phylogenetic 'species' identified ranging from two to 21. Such differing views have continued to confuse users of these taxonomies.

RESULTS

The phylogenetic relationships between Australian F. oxysporum isolates from both natural and agricultural ecosystems were determined using three datasets: whole genome, nuclear genes, and mitochondrial genome sequences. The phylogenies were concordant except for three isolates. There were three concordant clades from all the phylogenies suggesting similar evolutionary history for mitochondrial genome and nuclear genes for the isolates in these three clades. Applying a multispecies coalescent (MSC) model on the eight single copy nuclear protein coding genes from the nuclear gene dataset concluded that the three concordant clades correspond to three phylogenetic species within the FOSC. There was 100% posterior probability support for the formation of three species within the FOSC. This is the first report of using the MSC model to estimate species within the F. oxysporum species complex. The findings from this study were compared with previously published phylogenetics and species delimitation studies.

CONCLUSION

Phylogenetic analyses using three different gene datasets from Australian F. oxysporum isolates have all supported the formation of three major clades which delineated into three species. Species 2 (Clade 3) may be called F. oxysporum as it contains the neotype for F. oxysporum.

Co-evolution in the Jungle: From Leafcutter Ant Colonies to Chromosomal Ends

Biological entities are multicomponent systems where each part is directly or indirectly dependent on the others. In effect, a change in a single component might have a consequence on the functioning of its partners, thus affecting the fitness of the entire system. In this article, we provide a few examples of such complex biological systems, ranging from ant colonies to a population of amino acids within a single-polypeptide chain. Based on these examples, we discuss one of the central and still challenging questions in biology: how do such multicomponent consortia co-evolve? More specifically, we ask how telomeres, nucleo-protein complexes protecting the integrity of linear DNA chromosomes, originated from the ancestral organisms having circular genomes and thus not dealing with end-replication and end-protection problems. Using the examples of rapidly evolving topologies of mitochondrial genomes in eukaryotic microorganisms, we show what means of co-evolution were employed to accommodate various types of telomere-maintenance mechanisms in mitochondria. We also describe an unprecedented runaway evolution of telomeric repeats in nuclei of ascomycetous yeasts accompanied by co-evolution of telomere-associated proteins. We propose several scenarios derived from research on telomeres and supported by other studies from various fields of biology, while emphasizing that the relevant answers are still not in sight. It is this uncertainty and a lack of a detailed roadmap that makes the journey through the jungle of biological systems still exciting and worth undertaking.

Direct RNA sequencing approach to compare non-model mitochondrial transcriptomes: An application to a cephalopod host and its mesozoan parasite

To identify non-protein coding as well as truncated or premature RNA sequences expressed and obtain more complete transcriptome information, we combined the MinION direct RNA-sequencing of a conventional poly(A) RNA purification method with poly(A)-tagging of the non-coding RNA (ncRNA) fraction. This approach was applied to transcriptome sequencing of the dicyemid mesozoan, Dicyema misakiense, which has minicircular mitochondrial DNA molecules where each molecule encodes a single gene, as well as the host. Using informatics analysis, we distinguished dicyemid RNAs from those of the host squid. The poly(A) RNAs were assigned to host mitochondrial genes, host nuclear protein-coding genes, Dicyema nuclear protein-coding genes, and Dicyema mitochondrial genes in the decreasing order. Our poly(A)-tailing method recovered significantly more ncRNAs from the host compared with the sequencing of poly(A) RNAs. Furthermore, our method captured various lengths of squid mitochondrial DNA (mtDNA) transcripts at different steps of maturation including a read of 3,500 bp, which covers 21% of the squid mitochondrial genome, possibly a premature host RNA product. In contrast, shorter and less abundant reads were recovered from the dicyemid mitochondrial RNAs (mtRNAs). Even the longest read was 307 bp covering only a part of a minicircle. This study revealed significantly different modes of the mitochondrial transcription between a mesozoan and the host. Our approach to perform direct RNA-sequencing combined with the poly(A)-tailing reaction can be an effective method to fully capture non-poly(A) transcripts in a wide range of organisms.

Divergent Mitochondrial Ribosomes in Unicellular Parasitic Protozoans

The mitochondrion in parasitic protozoans is a clinically proven drug target. A specialized ribosome (mitoribosome) is required to translate genes encoded on the mitochondrial (mt) DNA. Despite the significance, little is known about mitoribosomes in many medically and economically important unicellular protozoans.

Pathogenic Mitochondria DNA Mutations: Current Detection Tools and Interventions

Mitochondria are best known for their role in energy production, and they are the only mammalian organelles that contain their own genomes. The mitochondrial genome mutation rate is reported to be 10–17 times higher compared to nuclear genomes as a result of oxidative damage caused by reactive oxygen species during oxidative phosphorylation. Pathogenic mitochondrial DNA mutations result in mitochondrial DNA disorders, which are among the most common inherited human diseases. Interventions of mitochondrial DNA disorders involve either the transfer of viable isolated mitochondria to recipient cells or genetically modifying the mitochondrial genome to improve therapeutic outcome. This review outlines the common mitochondrial DNA disorders and the key advances in the past decade necessary to improve the current knowledge on mitochondrial disease intervention. Although it is now 31 years since the first description of patients with pathogenic mitochondrial DNA was reported, the treatment for mitochondrial disease is often inadequate and mostly palliative. Advancements in diagnostic technology improved the molecular diagnosis of previously unresolved cases, and they provide new insight into the pathogenesis and genetic changes in mitochondrial DNA diseases.

Mitochondrial genome and diverse inheritance patterns in Pleurotus pulmonarius

Pleurotus pulmonarius, a member of the Pleurotaceae family in Basidiomycota, is an edible, economically important mushroom in most Asian countries. In this study, the complete mitochondrial genomes (mtDNA) of three P. pulmonarius strains - two monokaryotic commercial (J1-13 and ZA3) and one wild (X1-15) - were sequenced and analyzed. In ZA3 and X1-15, the mtDNA molecule was found to be a single circle of 68,305 bp and 73,435 bp, respectively. Both strains contain 14 core protein-coding genes and two ribosomal RNA (rRNA) subunit genes. The ZA3 strain has 22 transfer RNA (tRNA) genes and nine introns: eight in cytochrome c oxidase subunit 1 (coxl), and one in the rRNA large subunit (rnl). Monokaryotic J1-13 and ZA3 mtDNAs were found to be similar in their structure. However, the wild strain X1-15 contains 25 tRNA genes and only seven introns in coxl. Open reading frames (ORFs) of ZA3/J1-13 and X1-15 encode LAGLIDADG, ribosomal protein S3, and DNA polymerase II. In addition, mtDNA inheritance in J1-13, ZA3, and X1-15 was also studied. Results showed that the mtDNA inheritance pattern was uniparental and closely related to dikaryotic hyphal location with respect to the parent. Results also show that mtDNA inheritance is influenced by both the parental nuclear genome and mitogenome in the zone of contact between two compatible parents. In summary, this analysis provides valuable information and a basis for further studies to improve our understanding of the inheritance of fungal mtDNA.

Subchromosome-Scale Nuclear and Complete Mitochondrial Genome Characteristics of Morchella crassipes

Morchella crassipes (Vent.) Pers., a typical yellow morel species with high economic value, is mainly distributed in the low altitude plains of Eurasia. However, rare research has been performed on its genomics and polarity, thus limiting its research and development. Here, we reported a fine physical map of the nuclear genome at the subchromosomal-scale and the complete mitochondrial genome of M. crassipes. The complete size of the nuclear genome was 56.7 Mb, and 23 scaffolds were assembled, with eight of them being complete chromosomes. A total of 11,565 encoding proteins were predicted. The divergence time analysis showed that M. crassipes representing yellow morels differentiated with black morels at ~33.98 Mya (million years), with 150 gene families contracted and expanded in M. crassipes versus the two black morels (M. snyderi and M. importuna). Furthermore, 409 CAZYme genes were annotated in M. crassipes, containing almost all plant cell wall degrading enzymes compared with the mycorrhizal fungi (truffles). Genomic annotation of mating type loci and amplification of the mating genes in the monospore population was conducted, the results indicated that M. crassipes is a heterothallic fungus. Additionally, a complete circular mitochondrial genome of M. crassipes was assembled, the size reached as large as 531,195 bp. It can be observed that the strikingly large size was the biggest up till now, coupled with 14 core conserved mitochondrial protein-coding genes, two rRNAs, 31 tRNAs, 51 introns, and 412 ncORFs. The total length of intron sequences accounted for 53.67% of the mitochondrial genome, with 19 introns having a length over 5 kb. Particularly, 221 of 412 ncORFs were distributed within 51 introns, and the total length of the ncORFs sequence accounted for 40.83% of the mitochondrial genome, and 297 ncORFs had expression activity in the mycelium stage, suggesting their potential functions in M. crassipes. Meanwhile, there was a high degree of repetition (51.31%) in the mitochondria of M. crassipes. Thus, the large number of introns, ncORFs and internal repeat sequences may contribute jointly to the largest fungal mitochondrial genome to date. The fine physical maps of nuclear genome and mitochondrial genome obtained in this study will open a new door for better understanding of the mysterious species of M. crassipes.

The inflated mitochondrial genomes of siphonous green algae reflect processes driving expansion of noncoding DNA and proliferation of introns

Within the siphonous green algal order Bryopsidales, the size and gene arrangement of chloroplast genomes has been examined extensively, while mitochondrial genomes have been mostly overlooked. The recently published mitochondrial genome of Caulerpa lentillifera is large with expanded noncoding DNA, but it remains unclear if this is characteristic of the entire order. Our study aims to evaluate the evolutionary forces shaping organelle genome dynamics in the Bryopsidales based on the C. lentillifera and Ostreobium quekettii mitochondrial genomes. In this study, the mitochondrial genome of O. quekettii was characterised using a combination of long and short read sequencing, and bioinformatic tools for annotation and sequence analyses. We compared the mitochondrial and chloroplast genomes of O. quekettii and C. lentillifera to examine hypotheses related to genome evolution. The O. quekettii mitochondrial genome is the largest green algal mitochondrial genome sequenced (241,739 bp), considerably larger than its chloroplast genome. As with the mtDNA of C. lentillifera, most of this excess size is from the expansion of intergenic DNA and proliferation of introns. Inflated mitochondrial genomes in the Bryopsidales suggest effective population size, recombination and/or mutation rate, influenced by nuclear-encoded proteins, differ between the genomes of mitochondria and chloroplasts, reducing the strength of selection to influence evolution of their mitochondrial genomes.

Two Pentatricopeptide Repeat Proteins Are Required for the Splicing of nad5 Introns in Maize

Mitochondrial genes in flowering plants contain predominantly group II introns that require precise splicing before translation into functional proteins. Splicing of these introns is facilitated by various nucleus-encoded splicing factors. Due to lethality of mutants, functions of many splicing factors have not been revealed. Here, we report the function of two P-type PPR proteins PPR101 and PPR231, and their role in maize seed development. PPR101 and PPR231 are targeted to mitochondria. Null mutation of PPR101 and PPR231 arrests embryo and endosperm development, generating empty pericarp and small kernel phenotype, respectively, in maize. Loss-of-function in PPR101 abolishes the splicing of nad5 intron 2, and reduces the splicing of nad5 intron 1. Loss-of-function in PPR231 reduces the splicing of nad5 introns 1, 2, 3 and nad2 intron 3. The absence of Nad5 protein eliminates assembly of complex I, and activates the expression of alternative oxidase AOX2. These results indicate that both PPR101 and PPR231 are required for mitochondrial nad5 introns 1 and 2 splicing, while PPR231 is also required for nad5 intron 3 and nad2 intron 3. Both genes are essential to complex I assembly, mitochondrial function, and maize seed development. This work reveals that the splicing of a single intron involves multiple PPRs.

Lose-of-Function of a Rice Nucleolus-Localized Pentatricopeptide Repeat Protein Is Responsible for the floury endosperm14 Mutant Phenotypes

BACKGROUND

The endosperm of rice (Oryza sativa) has been usually used for the study of starch synthesis. Although several related factors have been revealed, other unknown members remain to be identified, given that starch synthesis is a complicated and sophisticated process.

RESULTS

Here, we identified and characterized a new rice seed mutant, floury endosperm14 (flo14), which showed chalked endosperm and seed-lethal phenotypes. Map-based cloning indicated FLO14 encodes a novel P-family PPR protein which contains ten PPR motifs. Afterwards the gene was named OsNPPR3. Subcellular localization showed OsNPPR3 was targeted to nucleolus. Quantitative RT-PCR analysis demonstrated that OsNPPR3 was universally expressed in various tissues, with pronounced levels during rice endosperm development. Molecular analysis further suggested that OsNPPR3 was involved in the regulation of expression levels and splicing of a few genes in mitochondria.

CONCLUSION

The study demonstrates that the nucleolus-localized PPR protein is responsible for the flo14 mutant phenotypes through affecting nuclear and mitochondrial gene expression and splicing.

DNA Replication in Human Mitochondria

DNA replication in human mitochondria has been studied for several decades; however, its mechanism still remains unclear. During the last 15 years, many new experimental data on the mitochondrial replication have appeared, although extremely contradictory. Two asynchronous (strand displacement and RITOLS) and one synchronous (strand-coupled) replication models have been proposed. In the asynchronous models, replication from the origin in the H-chain starts earlier, so that the replication of the two chains ends at different times. The synchronous model is more traditional and implies two replication forks with leading and lagging strands initiated at the same origin. For each of the three models, both confirming and contradicting experimental data exist. Most likely, there is no single model of mitochondrial replication. It is possible that the unique mitochondrial replication machinery that has originated as a results of endosymbiosis has an unexpected variety of replication strategies to maintain the mitochondrial genome. An unusual combination of enzymes of different origin (phage, bacterial, eukaryotic) and unique features of the mitochondrial genome (existance of heavy and light chains, insertions of ribonucleotides, a variety of origins) can allow replication through different mechanisms. In human mitochondria, asynchronous replication seems to dominate; however, synchronous replication is also possible under certain conditions. In the human heart mitochondria, circular mitochondrial DNA (mtDNA) molecules can rearrange in a network of rapidly replicating linear genomes, thereby suggesting possible existence of a wide range of replication mechanisms in the mitochondria. The review describes the main stages of mtDNA replication and enzymes involved in this process, as well as discusses the prospects of mitochondrial replication studies.

The mitochondrial genome of Morchella importuna (272.2 kb) is the largest among fungi and contains numerous introns, mitochondrial non-conserved open reading frames and repetitive sequences

The complete mitochondrial genome of Morchella importuna, the famous edible and medicinal mushroom, was assembled as a 272,238 bp single circular dsDNA. As the largest mitogenome among fungi, it exhibits several distinct characteristics. The mitogenome of M. importuna encoded 14 core conserved mitochondrial protein-coding genes and 151 mitochondrial non-conserved open reading frames (ncORFs) were predicted, of which 61 were annotated as homing endonuclease genes, and 108 were confirmed to be expressed during the vegetative growth stages of M. importuna. In addition, 34 introns were identified in seven core genes (cob, cox1, cox2, cox3, nad1, nad4 and nad5) and two rRNA genes (rrnS and rrnL) with a length from 383 bp to 7453 bp, and eight large introns with a length range of 2340 bp to 7453 bp contained multiple intronic mtORFs. Moreover, 34 group I (IA, IB, IC1, IC2, ID and derived group I introns) and four group II intron domains were identified for the 34 introns, including five hybrid ones. Furthermore, the M. importuna mitogenome showed the presence of about 18.7% mitogenomic interspersed repeats. These and the aforementioned ncORFs and introns, contributed to the enlarged size of the mitogenome.

Mapping and editing animal mitochondrial genomes: can we overcome the challenges?

The animal mitochondrial genome, although small, can have a big impact on health and disease. Non-pathogenic sequence variation among mitochondrial DNA (mtDNA) haplotypes influences traits including fertility, healthspan and lifespan, whereas pathogenic mutations are linked to incurable mitochondrial diseases and other complex conditions like ageing, diabetes, cancer and neurodegeneration. However, we know very little about how mtDNA genetic variation contributes to phenotypic differences. Infrequent recombination, the multicopy nature and nucleic acid-impenetrable membranes present significant challenges that hamper our ability to precisely map mtDNA variants responsible for traits, and to genetically modify mtDNA so that we can isolate specific mutants and characterize their biochemical and physiological consequences. Here, we summarize the past struggles and efforts in developing systems to map and edit mtDNA. We also assess the future of performing forward and reverse genetic studies on animal mitochondrial genomes. This article is part of the theme issue 'Linking the mitochondrial genotype to phenotype: a complex endeavour'.

Plant Organelle Genome Replication

Mitochondria and chloroplasts perform essential functions in respiration, ATP production, and photosynthesis, and both organelles contain genomes that encode only some of the proteins that are required for these functions. The proteins and mechanisms for organelle DNA replication are very similar to bacterial or phage systems. The minimal replisome may consist of DNA polymerase, a primase/helicase, and a single-stranded DNA binding protein (SSB), similar to that found in bacteriophage T7. In Arabidopsis, there are two genes for organellar DNA polymerases and multiple potential genes for SSB, but there is only one known primase/helicase protein to date. Genome copy number varies widely between type and age of plant tissues. Replication mechanisms are only poorly understood at present, and may involve multiple processes, including recombination-dependent replication (RDR) in plant mitochondria and perhaps also in chloroplasts. There are still important questions remaining as to how the genomes are maintained in new organelles, and how genome copy number is determined. This review summarizes our current understanding of these processes.

Comparative analysis of Chlorosarcinopsis eremi mitochondrial genome with some Chlamydomonadales algae

Chlorosarcinopsis eremi is a member of Chlamydomonadales algae which is isolated from terrestrial environments. In this study, the mitochondrial genome of C. eremi isolated from desert region of Iran, was represented for the first time. Following sequencing, assembly and annotation, comparative analyses of C. eremi and other available Chlamydomonadales algae complete mitochondrial genomes were performed. The mitochondrial genome of C. eremi was circular, had a low number of genes coding in the same strand with a minor amount of repeated sequences; same as other non-Reinhardtinia species of Chlamydomonadales algae. GC content of C. eremi mitochondrial genome was in normal range when compared with non-Chlamydomonadales organisms, but among Chlamydomonadales algae, C. eremi had a low GC content mitochondrial genome. C. eremi had the highest percent of non-coding sequences in comparison with other available Chlamydomonadales mitochondrial genomes which was related to intergenic regions. Identity analysis of protein-coding sequences of Chlamydomonadales mitochondrial genomes showed more divergences and may be related to the high mutation rate of mitochondrial genome reported in microbial eukaryotes.

Ancient Mitochondrial Gene Transfer between Fungi and the Orchids

The mitochondrial genomes (mitogenomes) of plants are known to incorporate and accumulate DNA from intra- and extracellular donors. Despite the intimate relationships formed between flowing plants (angiosperms) and fungi, lengthy fungal-like sequence has not been identified in angiosperm mitogenomes to date. Here, we present multiple lines of evidence documenting horizontal gene transfer (HGT) between the mitogenomes of fungi and the ancestors of the orchids, plants that are obligate parasites of fungi during their early development. We show that the ancestor of the orchids acquired an ∼270-bp fungal mitogenomic region containing three transfer RNA genes. We propose that the short HGT was later replaced by a second HGT event transferring >8 kb and 14 genes from a fungal mitogenome to that of the ancestor of the largest orchid subfamily, Epidendroideae. Our results represent the first evidence of genomic-scale HGT between fungal and angiosperm mitogenomes and demonstrate that the length intergenic spacer regions of angiosperm mitogenomes can effectively fossilize the genomic remains of ancient, nonplant organisms.

Maize pentatricopeptide repeat protein DEK41 affects cis-splicing of mitochondrial nad4 intron 3 and is required for normal seed development

The splicing of organelle-encoded mRNA in plants requires proteins encoded in the nucleus. The mechanism of splicing and the factors involved are not well understood. Pentatricopeptide repeat (PPR) proteins are known to participate in such RNA-protein interactions. Maize defective kernel 41 (dek41) is a seedling-lethal mutant that causes developmental defects. In this study, the Dek41 gene was cloned by Mutator tag isolation and allelic confirmation, and was found to encode a P-type PPR protein that targets mitochondria. Analysis of the mitochondrial RNA transcript profile revealed that dek41 mutations cause reduced splicing efficiency of mitochondrial nad4 intron 3. Immature dek41 kernels exhibited severe reductions in complex I assembly and NADH dehydrogenase activity. Up-regulated expression of alternative oxidase genes and deformed inner cristae of mitochondria in dek41, as revealed by TEM, indicated that proper splicing of nad4 is essential for correct mitochondrial functioning and morphology. Consistent with this finding, differentially expressed genes in the dek41 endosperm included those related to mitochondrial function and activity. Our results indicate that DEK41 is a PPR protein that affects cis-splicing of mitochondrial nad4 intron 3 and is required for correct mitochondrial functioning and maize kernel development.

The translational landscape of Arabidopsis mitochondria

Messenger RNA translation is a complex process that is still poorly understood in eukaryotic organelles like mitochondria. Growing evidence indicates though that mitochondrial translation differs from its bacterial counterpart in many key aspects. In this analysis, we have used ribosome profiling technology to generate a genome-wide snapshot view of mitochondrial translation in Arabidopsis. We show that, unlike in humans, most Arabidopsis mitochondrial ribosome footprints measure 27 and 28 bases. We also reveal that respiratory subunits encoding mRNAs show much higher ribosome association than other mitochondrial mRNAs, implying that they are translated at higher levels. Homogenous ribosome densities were generally detected within each respiratory complex except for complex V, where higher ribosome coverage corroborated with higher requirements for specific subunits. In complex I respiratory mutants, a reorganization of mitochondrial mRNAs ribosome association was detected involving increased ribosome densities for certain ribosomal protein encoding transcripts and a reduction in translation of a few complex V mRNAs. Taken together, our observations reveal that plant mitochondrial translation is a dynamic process and that translational control is important for gene expression in plant mitochondria. This study paves the way for future advances in the understanding translation in higher plant mitochondria.

Key within-membrane residues and precursor dosage impact the allotopic expression of yeast subunit II of cytochrome c oxidase

Experimentally relocating mitochondrial genes to the nucleus for functional expression (allotopic expression) is a challenging process. The high hydrophobicity of mitochondria-encoded proteins seems to be one of the main factors preventing this allotopic expression. We focused on subunit II of cytochrome c oxidase (Cox2) to study which modifications may enable or improve its allotopic expression in yeast. Cox2 can be imported from the cytosol into mitochondria in the presence of the W56R substitution, which decreases the protein hydrophobicity and allows partial respiratory rescue of a cox2-null strain. We show that the inclusion of a positive charge is more favorable than substitutions that only decrease the hydrophobicity. We also searched for other determinants enabling allotopic expression in yeast by examining the COX2 gene in organisms where it was transferred to the nucleus during evolution. We found that naturally occurring variations at within-membrane residues in the legume Glycine max Cox2 could enable yeast COX2 allotopic expression. We also evidence that directing high doses of allotopically synthesized Cox2 to mitochondria seems to be counterproductive because the subunit aggregates at the mitochondrial surface. Our findings are relevant to the design of allotopic expression strategies and contribute to the understanding of gene retention in organellar genomes.

Enrichment of G4DNA and a Large Inverted Repeat Coincide in the Mitochondrial Genomes of Termitomyces

Mitochondria retain their own genome, a hallmark of their bacterial ancestry. Mitochondrial genomes (mtDNA) are highly diverse in size, shape, and structure, despite their conserved function across most eukaryotes. Exploring extreme cases of mtDNA architecture can yield important information on fundamental aspects of genome biology. We discovered that the mitochondrial genomes of a basidiomycete fungus (Termitomyces spp.) contain an inverted repeat (IR), a duplicated region half the size of the complete genome. In addition, we found an abundance of sequences capable of forming G-quadruplexes (G4DNA); structures that can disrupt the double helical formation of DNA. G4DNA is implicated in replication fork stalling, double-stranded breaks, altered gene expression, recombination, and other effects. To determine whether this occurrence of IR and G4DNA was correlated within the genus Termitomyces, we reconstructed the mitochondrial genomes of 11 additional species including representatives of several closely related genera. We show that the mtDNA of all sampled species of Termitomyces and its sister group, represented by the species Tephrocybe rancida and Blastosporella zonata, are characterized by a large IR and enrichment of G4DNA. To determine whether high mitochondrial G4DNA content is common in fungi, we conducted the first broad survey of G4DNA content in fungal mtDNA, revealing it to be a highly variable trait. The results of this study provide important direction for future research on the function and evolution of G4DNA and organellar IRs.

Evolution: In Chloroplast Genomes, Anything Goes

A new study shows that Cladophorales green algae have the most unconventional chloroplast DNAs ever observed, whereby genes are located on small linear single-stranded palindromic elements. This puzzling architecture has parallels with mini-circular chloroplast genomes of dinoflagellates and raises many questions about how it arose and is maintained.

Mobile genetic elements explain size variation in the mitochondrial genomes of four closely-related Armillaria species

Background

Species in the genus Armillaria (fungi, basidiomycota) are well-known as saprophytes and pathogens on plants. Many of them cause white-rot root disease in diverse woody plants worldwide. Mitochondrial genomes (mitogenomes) are widely used in evolutionary and population studies, but despite the importance and wide distribution of Armillaria, the complete mitogenomes have not previously been reported for this genus. Meanwhile, the well-supported phylogeny of Armillaria species provides an excellent framework in which to study variation in mitogenomes and how they have evolved over time.

Results

Here we completely sequenced, assembled, and annotated the circular mitogenomes of four species: A. borealis, A. gallica, A. sinapina, and A. solidipes (116,443, 98,896, 103,563, and 122,167 bp, respectively). The variation in mitogenome size can be explained by variable numbers of mobile genetic elements, introns, and plasmid-related sequences. Most Armillaria introns contained open reading frames (ORFs) that are related to homing endonucleases of the LAGLIDADG and GIY-YIG families. Insertions of mobile elements were also evident as fragments of plasmid-related sequences in Armillaria mitogenomes. We also found several truncated gene duplications in all four mitogenomes.

Conclusions

Our study showed that fungal mitogenomes have a high degree of variation in size, gene content, and genomic organization even among closely related species of Armillara. We suggest that mobile genetic elements invading introns and intergenic sequences in the Armillaria mitogenomes have played a significant role in shaping their genome structure. The mitogenome changes we describe here are consistent with widely accepted phylogenetic relationships among the four species.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5732-z) contains supplementary material, which is available to authorized users.

Mitochondrial mRNA translation initiation contributes to oxidative metabolism in the myocardia of aged, obese mice

Being both advanced in age and obese each contribute to cardiac hypertrophy in a unique manner. Electron transport complexes I and IV are implicated in deficient electron transport during cardiomyopathies and contain the majority of protein subunits that are transcribed and translated by machinery localized within the mitochondria.

PURPOSE

To assess myocardial mt-mRNA translation factors in relation to mitochondrial content and mtDNA-encoded protein using a mouse model of aged obesity and to test the relationship of mt-mRNA translation initiation factor 2 (mtIF2) to oxidative capacity and the cellular oxidation-reduction (redox) state in cardiomyocytes.

METHODS

Male C56BL/6 J mice fed lean or high fat diet were aged to either ~3 months or ~22 months, the heart was excised and analyzed using immunoblot and qPCR to assess differences in mitochondrial mRNA translation machinery. Using H9c2 cardiomyocytes, mtIF2 was knocked-down and oxidative metabolic characteristics assessed including oxidation/reduction state, bioenergetic flux, and hypoxic resistance was tested.

RESULTS

Aged, obese mouse hearts were ~40% larger than young, lean controls and contained ~50% less mtIF2 protein alongside ~25-50% lower content of Cytb, a protein encoded by mtDNA. Reducing the level of mtIF2 by shRNA is associated with ~15-20% lower content of OXPHOS complex I and IV, ~30% lower optical redox ratio, ~40% oxygen reserve capacity, and ~20% less cell survival following hypoxia.

CONCLUSION

We present evidence of altered mt-mRNA translation during cardiac hypertrophy in aged obesity. We build on these results by demonstrating the necessity of mtIF2 in maintaining oxidative characteristics of cardiac muscle cells.

Mitochondrial presequence import: Multiple regulatory knobs fine-tune mitochondrial biogenesis and homeostasis

Mitochondria are pivotal organelles for cellular signaling and metabolism, and their dysfunction leads to severe cellular stress. About 60-70% of the mitochondrial proteome consists of preproteins synthesized in the cytosol with an amino-terminal cleavable presequence targeting signal. The TIM23 complex transports presequence signals towards the mitochondrial matrix. Ultimately, the mature protein segments are either transported into the matrix or sorted to the inner membrane. To ensure accurate preprotein import into distinct mitochondrial sub-compartments, the TIM23 machinery adopts specific functional conformations and interacts with different partner complexes. Regulatory subunits modulate the translocase dynamics, tailoring the import reaction to the incoming preprotein. The mitochondrial membrane potential and the ATP generated via oxidative phosphorylation are key energy sources in driving the presequence import pathway. Thus, mitochondrial dysfunctions have rapid repercussions on biogenesis. Cellular mechanisms exploit the presequence import pathway to monitor mitochondrial dysfunctions and mount transcriptional and proteostatic responses to restore functionality.

[The underestimated coding potential of mitochondrial DNA]

Mitochondria are ancient organelles that emerged from the endosymbiosis of free-living proto-bacteria. They still retain a semi-autonomous genetic system with a small genome. Mitochondrial DNA (mtDNA) codes for 13 essential proteins for the production of ATP, the sequences of which are relatively conserved across Metazoans. The discovery of additional mitochondria-derived peptides (MDPs) indicates an underestimated coding potential. Humanin, an anti-apoptotic peptide, is likely independently transcribed from within the 16S rRNA gene, as are recently described SHLPs. MOTS-c, discovered in silico, has been demonstrated to be involved in metabolism and insulin sensitivity. Gau, is a positionally conserved open reading frame (ORF) sequence found in the antisense strand of the COX1 gene and its corresponding peptide is strictly colocalized with mitochondrial markers. In bivalves with doubly uniparental inheritance of mtDNA, male and female mtDNAs each carry a separate additional gene possibly involved in sex determination. Other MDPs likely exist and their investigation will shed light on the underestimated functional repertoire of mitochondria.

Release of Mitochondrial and Nuclear DNA During On-Pump Heart Surgery: Kinetics and Relation to Extracellular Vesicles

During heart surgery with cardiopulmonary bypass (CPB), the release of mitochondrial (mtDNA) and nuclear DNA (nDNA) and their association to extracellular vesicles were investigated. In patients undergoing elective coronary artery bypass grafting (CABG, n = 12), blood was sampled before, during, and after surgery from peripheral artery, pulmonary artery, and the coronary sinus. Plasma was separated in three fractions: microvesicles, exosomes, and supernatant. mtDNA and nDNA were measured by qPCR. mtDNA and nDNA levels increased after start of surgery, but before CPB, and increased further during CPB. mtDNA copy number was about 1000-fold higher than nDNA. mtDNA was predominantly localized to the vesicular fractions in plasma, whereas nDNA was predominantly in the supernatant. The amount of free mtDNA increased after surgery. There was no net release or disappearance of DNAs across the pulmonary, systemic, or coronary circulation. Extracellular DNAs, in particular mtDNA, may be important contributors to the whole-body inflammation during CPB.

Complete mitochondrial genome of Benthodytes marianensis (Holothuroidea: Elasipodida: Psychropotidae): Insight into deep sea adaptation in the sea cucumber

Complete mitochondrial genomes play important roles in studying genome evolution, phylogenetic relationships, and species identification. Sea cucumbers (Holothuroidea) are ecologically important and diverse members, living from the shallow waters to the hadal trench. In this study, we present the mitochondrial genome sequence of the sea cucumber Benthodytes marianensis collected from the Mariana Trench. To our knowledge, this is the first reported mitochondrial genome from the genus Benthodytes. This complete mitochondrial genome is 17567 bp in length and consists of 13 protein-coding genes, two ribosomal RNA genes and 22 transfer RNA genes (duplication of two tRNAs: trnL and trnS). Most of these genes are coded on the positive strand except for one protein-coding gene (nad6) and five tRNA genes which are coded on the negative strand. Two putative control regions (CRs) have been found in the B. marianensis mitogenome. We compared the order of genes from the 10 available holothurian mitogenomes and found a novel gene arrangement in B. marianensis. Phylogenetic analysis revealed that B. marianensis clustered with Peniagone sp. YYH-2013, forming the deep-sea Elasipodida clade. Positive selection analysis showed that eleven residues (24 S, 45 S, 185 S, 201 G, 211 F and 313 N in nad2; 108 S, 114 S, 322 C, 400 T and 427 S in nad4) were positively selected sites with high posterior probabilities. We predict that nad2 and nad4 may be the important candidate genes for the further investigation of the adaptation of B. marianensis to the deep-sea environment.

Genetic Drift and Indel Mutation in the Evolution of Yeast Mitochondrial Genome Size

Mitochondrial genomes (mitogenomes) are remarkably diverse in genome size and organization, but the origins of dynamic mitogenome architectures are still poorly understood. For instance, the mutational burden hypothesis postulates that the drastic difference between large plant mitogenomes and streamlined animal mitogenomes can be driven by their different mutation rates. However, inconsistent trends between mitogenome sizes and mutation rates have been documented in several lineages. These conflicting results highlight the need of systematic and sophisticated investigations on the evolution and diversity of mitogenome architecture. This study took advantage of the strikingly variable mitogenome size among different yeast species and also among intraspecific strains, examined sequence dynamics of introns, GC-clusters, tandem repeats, mononucleotide repeats (homopolymers) and evaluated their contributions to genome size variation. The contributions of these sequence features to mitogenomic variation are dependent on the timescale, over which extant genomes evolved from their last common ancestor, perhaps due to a combination of different turnover rates of mobile sequences, variable insertion spaces, and functional constraints. We observed a positive correlation between mitogenome size and the level of genetic drift, suggesting that mitogenome expansion in yeast is likely driven by multiple types of sequence insertions in a primarily nonadaptive manner. Although these cannot be explained directly by the mutational burden hypothesis, our results support an important role of genetic drift in the evolution of yeast mitogenomes.

Annotation and analysis of the mitochondrial genome of Coniothyrium glycines, causal agent of red leaf blotch of soybean, reveals an abundance of homing endonucleases

Coniothyrium glycines, the causal agent of soybean red leaf blotch, is a USDA APHIS-listed Plant Pathogen Select Agent and potential threat to US agriculture. Sequencing of the C. glycines mt genome revealed a circular 98,533-bp molecule with a mean GC content of 29.01%. It contains twelve of the mitochondrial genes typically involved in oxidative phosphorylation (atp6, cob, cox1-3, nad1-6, and nad4L), one for a ribosomal protein (rps3), four for hypothetical proteins, one for each of the small and large subunit ribosomal RNAs (rns and rnl) and a set of 30 tRNAs. Genes were encoded on both DNA strands with cox1 and cox2 occurring as adjacent genes having no intergenic spacers. Likewise, nad2 and nad3 are adjacent with no intergenic spacers and nad5 is immediately followed by nad4L with an overlap of one base. Thirty-two introns, comprising 54.1% of the total mt genome, were identified within eight protein-coding genes and the rnl. Eighteen of the introns contained putative intronic ORFs with either LAGLIDADG or GIY-YIG homing endonuclease motifs, and an additional eleven introns showed evidence of truncated or degenerate endonuclease motifs. One intron possessed a degenerate N-acetyl-transferase domain. C. glycines shares some conservation of gene order with other members of the Pleosporales, most notably nad6-rnl-atp6 and associated conserved tRNA clusters. Phylogenetic analysis of the twelve shared protein coding genes agrees with commonly accepted fungal taxonomy. C. glycines represents the second largest mt genome from a member of the Pleosporales sequenced to date. This research provides the first genomic information on C. glycines, which may provide targets for rapid diagnostic assays and population studies.

Mitochondrial and Plastid Genomes from Coralline Red Algae Provide Insights into the Incongruent Evolutionary Histories of Organelles

Mitochondria and plastids are generally uniparentally inherited and have a conserved gene content over hundreds of millions of years, which makes them potentially useful phylogenetic markers. Organelle single gene-based trees have long been the basis for elucidating interspecies relationships that inform taxonomy. More recently, high-throughput genome sequencing has enabled the construction of massive organelle genome databases from diverse eukaryotes, and these have been used to infer species relationships in deep evolutionary time. Here, we test the idea that despite their expected utility, conflicting phylogenetic signal may exist in mitochondrial and plastid genomes from the anciently diverged coralline red algae (Rhodophyta). We generated complete organelle genome data from five coralline red algae (Lithothamnion sp., Neogoniolithon spectabile, Renouxia sp., Rhodogorgon sp., and Synarthrophyton chejuensis) for comparative analysis with existing organelle genome data from two other species (Calliarthron tuberculosum and Sporolithon durum). We find strong evidence for incongruent phylogenetic signal from both organelle genomes that may be explained by incomplete lineage sorting that has maintained anciently derived gene copies or other molecular evolutionary processes such as hybridization or gene flow during the evolutionary history of coralline red algae.

Insertion Defects of Mitochondrially Encoded Proteins Burden the Mitochondrial Quality Control System

The mitochondrial proteome contains proteins from two different genetic systems. Proteins are either synthesized in the cytosol and imported into the different compartments of the organelle or directly produced in the mitochondrial matrix. To ensure proteostasis, proteins are monitored by the mitochondrial quality control system, which will degrade non-native polypeptides. Defective mitochondrial membrane proteins are degraded by membrane-bound AAA-proteases. These proteases are regulated by factors promoting protein turnover or preventing their degradation. Here we determined genetic interactions between the mitoribosome receptors Mrx15 and Mba1 with the quality control system. We show that simultaneous absence of Mrx15 and the regulators of the i-AAA protease Mgr1 and Mgr3 provokes respiratory deficiency. Surprisingly, mutants lacking Mrx15 were more tolerant against proteotoxic stress. Furthermore, yeast cells became hypersensitive against proteotoxic stress upon deletion of MBA1. Contrary to Mrx15, Mba1 cooperates with the regulators of the m-AAA and i-AAA proteases. Taken together, these results suggest that membrane protein insertion and mitochondrial AAA-proteases are functionally coupled, possibly reflecting an early quality control step during mitochondrial protein synthesis.

Perfection of eccentricity: Mitochondrial genomes of diplonemids

Mitochondria are the sandbox of evolution as exemplified most particularly by the diplonemids, a group of marine microeukaryotes. These protists are uniquely characterized by their highly multipartite mitochondrial genome and systematically fragmented genes whose pieces are spread out over several dozens of chromosomes. The type species Diplonema papillatum was the first member of this group in which the expression of fragmented mitochondrial genes was investigated experimentally. We now know that gene expression involves separate transcription of gene pieces (modules), RNA editing of module transcripts, and module joining to mature mRNAs and rRNAs. The mechanism of cognate module recognition and ligation is distinct from known intron splicing and remains to be uncovered. Here, we review the current status of research on mitochondrial genome architecture, as well as gene complement, structure, and expression modes in diplonemids. Further, we discuss the potential molecular mechanisms of posttranscriptional processing, and finally reflect on the evolutionary trajectories and trends of mtDNA evolution as seen in this protist group. © 2018 IUBMB Life, 70(12):1197-1206, 2018.

Massive mitochondrial DNA content in diplonemid and kinetoplastid protists

The mitochondrial DNA of diplonemid and kinetoplastid protists is known for its suite of bizarre features, including the presence of concatenated circular molecules, extensive trans‐splicing and various forms of RNA editing. Here we report on the existence of another remarkable characteristic: hyper‐inflated DNA content. We estimated the total amount of mitochondrial DNA in four kinetoplastid species (Trypanosoma brucei, Trypanoplasma borreli, Cryptobia helicis, and Perkinsela sp.) and the diplonemid Diplonema papillatum. Staining with 4′,6‐diamidino‐2‐phenylindole and RedDot1 followed by color deconvolution and quantification revealed massive inflation in the total amount of DNA in their organelles. This was further confirmed by electron microscopy. The most extreme case is the ∼260 Mbp of DNA in the mitochondrion of Diplonema, which greatly exceeds that in its nucleus; this is, to our knowledge, the largest amount of DNA described in any organelle. Perkinsela sp. has a total mitochondrial DNA content ~6.6× greater than its nuclear genome. This mass of DNA occupies most of the volume of the Perkinsela cell, despite the fact that it contains only six protein‐coding genes. Why so much DNA? We propose that these bloated mitochondrial DNAs accumulated by a ratchet‐like process. Despite their excessive nature, the synthesis and maintenance of these mtDNAs must incur a relatively low cost, considering that diplonemids are one of the most ubiquitous and speciose protist groups in the ocean. © 2018 The Authors. IUBMB Life published by Wiley Periodicals, Inc. on behalf of International Union of Biochemistry and Molecular Biology., 70(12):1267–1274, 2018

The Multiple Levels of Mitonuclear Coregulation

Together, the nuclear and mitochondrial genomes encode the oxidative phosphorylation (OXPHOS) complexes that reside in the mitochondrial inner membrane and enable aerobic life. Mitochondria maintain their own genome that is expressed and regulated by factors distinct from their nuclear counterparts. For optimal function, the cell must ensure proper stoichiometric production of OXPHOS subunits by coordinating two physically separated and evolutionarily distinct gene expression systems. Here, we review our current understanding of mitonuclear coregulation primarily at the levels of transcription and translation. Additionally, we discuss other levels of coregulation that may exist but remain largely unexplored, including mRNA modification and stability and posttranslational protein degradation.

The first complete mitochondrial genome from the family Hygrophoraceae (Hygrophorus russula) by next-generation sequencing and phylogenetic implications

Hygrophorus russula (Schaeff.) Kauffman is an edible ectomycorrhizal fungus that is widely distributed in the world. In this study, the mitogenome of H. russula was sequenced and assembled. The mitogenome of H. russula is composed of circular DNA molecules, with a total size of 55,769 bp. Further analysis indicated that the frequent use of A and T in codons contributes to the high AT content (80.87%) in the H. russula mitogenome. Comparative analysis indicated that the length and base composition of the core protein-encoding genes, and the number of tRNA genes in the H. russula mitogenome varied from that of other Agaricales mitogenomes. Gene arrangement analysis revealed a novel gene order in the H. russula mitogenome. In addition, the expansion of the mitogenome in Agaricales was found to be closely related to the increase in the number of introns. Phylogenetic analysis of the combined mitochondrial gene set showed strong support for tree topologies, and H. russula was determined to be relatively distant from other Agaricales species. This study is the first report on the mitogenome of a member of genus Hygrophorus as well as family Hygrophoraceae, which improves our understanding of mitochondrial differentiation and evolution in the important ectomycorrhizal fungi Hygrophorus species.

Fungal mitochondrial genomes and genetic polymorphisms

Mitochondria are the powerhouses of eukaryotic cells, responsible for ATP generation and playing a role in a diversity of cellular and organismal functions. Different from the majority of other intracellular membrane structures, mitochondria contain their own genetic materials that are capable of independent replication and inheritance. In this mini-review, we provide brief summaries of fungal mitochondrial genome structure, size, gene content, inheritance, and genetic variation. We pay special attention to the relative genetic polymorphisms of the mitochondrial vs nuclear genomes at the population level within individual fungal species. Among the 20 species/groups of species reviewed here, there is a range of variation among genes and species in the relative nuclear and mitochondrial genetic polymorphisms. Interestingly, most (15/20) showed a greater genetic diversity for nuclear genes and genomes than for mitochondrial genes and genomes, with the remaining five showing similar or slower nuclear genome genetic variations. This fungal pattern is different from the dominant pattern in animals, but more similar to that in plants. At present, the mechanisms for the variations among fungal species and the overall low level of mitochondrial sequence polymorphisms are not known. The increasing availability of population genomic data should help us reveal the potential genetic and ecological factors responsible for the observed variations.

The linear mitochondrial genome of the quarantine chytrid Synchytrium endobioticum; insights into the evolution and recent history of an obligate biotrophic plant pathogen

BACKGROUND

Chytridiomycota species (chytrids) belong to a basal lineage in the fungal kingdom. Inhabiting terrestrial and aquatic environments, most are free-living saprophytes but several species cause important diseases: e.g. Batrachochytrium dendrobatidis, responsible for worldwide amphibian decline; and Synchytrium endobioticum, causing potato wart disease. S. endobioticum has an obligate biotrophic lifestyle and isolates can be further characterized as pathotypes based on their virulence on a differential set of potato cultivars. Quarantine measures have been implemented globally to control the disease and prevent its spread. We used a comparative approach using chytrid mitogenomes to determine taxonomical relationships and to gain insights into the evolution and recent history of introductions of this plant pathogen.

RESULTS

We assembled and annotated the complete mitochondrial genome of 30 S. endobioticum isolates and generated mitochondrial genomes for five additional chytrid species. The mitochondrial genome of S. endobioticum is linear with terminal inverted repeats which was validated by tailing and PCR amplifying the telomeric ends. Surprisingly, no conservation in organisation and orientation of mitochondrial genes was observed among the Chytridiomycota except for S. endobioticum and its sister species Synchytrium microbalum. However, the mitochondrial genome of S. microbalum is circular and comprises only a third of the 72.9 Kbp found for S. endobioticum suggesting recent linearization and expansion. Four mitochondrial lineages were identified in the S. endobioticum mitochondrial genomes. Several pathotypes occur in different lineages, suggesting that these have emerged independently. In addition, variations for polymorphic sites in the mitochondrial genome of individual isolates were observed demonstrating that S. endobioticum isolates represent a community of different genotypes. Such communities were shown to be complex and stable over time, but we also demonstrate that the use of semi-resistant potato cultivars triggers a rapid shift in the mitochondrial haplotype associated with increased virulence.

CONCLUSIONS

Mitochondrial genomic variation shows that S. endobioticum has been introduced into Europe multiple times, that several pathotypes emerged multiple times, and that isolates represent communities of different genotypes. Our study represents the most comprehensive dataset of chytrid mitogenomes, which provides new insights into the extraordinary dynamics and evolution of mitochondrial genomes involving linearization, expansion and reshuffling.

The complete mitochondrial genome of the early flowering plant Nymphaea colorata is highly repetitive with low recombination

BACKGROUND

Mitochondrial genomes of flowering plants (angiosperms) are highly dynamic in genome structure. The mitogenome of the earliest angiosperm Amborella is remarkable in carrying rampant foreign DNAs, in contrast to Liriodendron, the other only known early angiosperm mitogenome that is described as 'fossilized'. The distinctive features observed in the two early flowering plant mitogenomes add to the current confusions of what early flowering plants look like. Expanded sampling would provide more details in understanding the mitogenomic evolution of early angiosperms. Here we report the complete mitochondrial genome of water lily Nymphaea colorata from Nymphaeales, one of the three orders of the earliest angiosperms.

RESULTS

Assembly of data from Pac-Bio long-read sequencing yielded a circular mitochondria chromosome of 617,195 bp with an average depth of 601×. The genome encoded 41 protein coding genes, 20 tRNA and three rRNA genes with 25 group II introns disrupting 10 protein coding genes. Nearly half of the genome is composed of repeated sequences, which contributed substantially to the intron size expansion, making the gross intron length of the Nymphaea mitochondrial genome one of the longest among angiosperms, including an 11.4-Kb intron in cox2, which is the longest organellar intron reported to date in plants. Nevertheless, repeat mediated homologous recombination is unexpectedly low in Nymphaea evidenced by 74 recombined reads detected from ten recombinationally active repeat pairs among 886,982 repeat pairs examined. Extensive gene order changes were detected in the three early angiosperm mitogenomes, i.e. 38 or 44 events of inversions and translocations are needed to reconcile the mitogenome of Nymphaea with Amborella or Liriodendron, respectively. In contrast to Amborella with six genome equivalents of foreign mitochondrial DNA, not a single horizontal gene transfer event was observed in the Nymphaea mitogenome.

CONCLUSIONS

The Nymphaea mitogenome resembles the other available early angiosperm mitogenomes by a similarly rich 64-coding gene set, and many conserved gene clusters, whereas stands out by its highly repetitive nature and resultant remarkable intron expansions. The low recombination level in Nymphaea provides evidence for the predominant master conformation in vivo with a highly substoichiometric set of rearranged molecules.

Pervasive Transcription of Mitochondrial, Plastid, and Nucleomorph Genomes across Diverse Plastid-Bearing Species

Organelle genomes exhibit remarkable diversity in content, structure, and size, and in their modes of gene expression, which are governed by both organelle- and nuclear-encoded machinery. Next generation sequencing (NGS) has generated unprecedented amounts of genomic and transcriptomic data, which can be used to investigate organelle genome transcription. However, most of the available eukaryotic RNA-sequencing (RNA-seq) data are used to study nuclear transcription only, even though large numbers of organelle-derived reads can typically be mined from these experiments. Here, we use publicly available RNA-seq data to assess organelle genome transcription in 59 diverse plastid-bearing species. Our RNA mapping analyses unraveled pervasive (full or near-full) transcription of mitochondrial, plastid, and nucleomorph genomes. In all cases, 85% or more of the organelle genome was recovered from the RNA data, including noncoding (intergenic and intronic) regions. These results reinforce the idea that organelles transcribe all or nearly all of their genomic material and are dependent on post-transcriptional processing of polycistronic transcripts. We explore the possibility that transcribed intergenic regions are producing functional noncoding RNAs, and that organelle genome noncoding content might provide raw material for generating regulatory RNAs.