Genes and junk in plant mitochondria-repair mechanisms and selection

Evolution and maintenance of mtDNA gene content across eukaryotes

Across eukaryotes, most genes required for mitochondrial function have been transferred to, or otherwise acquired by, the nucleus. Encoding genes in the nucleus has many advantages. So why do mitochondria retain any genes at all? Why does the set of mtDNA genes vary so much across different species? And how do species maintain functionality in the mtDNA genes they do retain? In this review, we will discuss some possible answers to these questions, attempting a broad perspective across eukaryotes. We hope to cover some interesting features which may be less familiar from the perspective of particular species, including the ubiquity of recombination outside bilaterian animals, encrypted chainmail-like mtDNA, single genes split over multiple mtDNA chromosomes, triparental inheritance, gene transfer by grafting, gain of mtDNA recombination factors, social networks of mitochondria, and the role of mtDNA dysfunction in feeding the world. We will discuss a unifying picture where organismal ecology and gene-specific features together influence whether organism X retains mtDNA gene Y, and where ecology and development together determine which strategies, importantly including recombination, are used to maintain the mtDNA genes that are retained.

Expansion of the MutS Gene Family in Plants

The MutS gene family is distributed across the tree of life and is involved in recombination, DNA repair, and protein translation. Multiple evolutionary processes have expanded the set of MutS genes in plants relative to other eukaryotes. Here, we investigate the origins and functions of these plant-specific genes. Land plants, green algae, red algae, and glaucophytes share cyanobacterial-like MutS1 and MutS2 genes that presumably were gained via plastid endosymbiotic gene transfer. MutS1 was subsequently lost in some taxa, including seed plants, whereas MutS2 was duplicated in Viridiplantae (i.e., land plants and green algae) with widespread retention of both resulting paralogs. Viridiplantae also have two anciently duplicated copies of the eukaryotic MSH6 gene (i.e., MSH6 and MSH7) and acquired MSH1 via horizontal gene transfer - potentially from a nucleocytovirus. Despite sharing the same name, "plant MSH1" is not directly related to the gene known as MSH1 in some fungi and animals, which may be an ancestral eukaryotic gene acquired via mitochondrial endosymbiosis and subsequently lost in most eukaryotic lineages. There has been substantial progress in understanding the functions of MSH1 and MSH6/MSH7 in plants, but the roles of the cyanobacterial-like MutS1 and MutS2 genes remain uncharacterized. Known functions of bacterial homologs and predicted protein structures, including fusions to diverse nuclease domains, provide hypotheses about potential molecular mechanisms. Because most plant-specific MutS proteins are targeted to the mitochondria and/or plastids, the expansion of this family appears to have played a large role in shaping plant organelle genetics.

Disruption of recombination machinery alters the mutational landscape in plant organellar genomes

Land plant organellar genomes have extremely low rates of point mutation yet also experience high rates of recombination and genome instability. Characterizing the molecular machinery responsible for these patterns is critical for understanding the evolution of these genomes. While much progress has been made towards understanding recombination activity in land plant organellar genomes, the relationship between recombination pathways and point mutation rates remains uncertain. The organellar targeted mutS homolog MSH1 has previously been shown to suppress point mutations as well as non-allelic recombination between short repeats in Arabidopsis thaliana. We therefore implemented high-fidelity Duplex Sequencing to test if other genes that function in recombination and maintenance of genome stability also affect point mutation rates. We found small to moderate increases in the frequency of single nucleotide variants (SNVs) and indels in mitochondrial and/or plastid genomes of A. thaliana mutant lines lacking radA, recA1, or recA3. In contrast, osb2 and why2 mutants did not exhibit an increase in point mutations compared to wild type (WT) controls. In addition, we analyzed the distribution of SNVs in previously generated Duplex Sequencing data from A. thaliana organellar genomes and found unexpected strand asymmetries and large effects of flanking nucleotides on mutation rates in WT plants and msh1 mutants. Finally, using long-read Oxford Nanopore sequencing, we characterized structural variants in organellar genomes of the mutant lines and show that different short repeat sequences become recombinationally active in different mutant backgrounds. Together, these complementary sequencing approaches shed light on how recombination may impact the extraordinarily low point mutation rates in plant organellar genomes.

Plant Organellar MSH1 Is a Displacement Loop-Specific Endonuclease

MutS HOMOLOG 1 (MSH1) is an organellar-targeted protein that obstructs ectopic recombination and the accumulation of mutations in plant organellar genomes. MSH1 also modulates the epigenetic status of nuclear DNA, and its absence induces a variety of phenotypic responses. MSH1 is a member of the MutS family of DNA mismatch repair proteins but harbors an additional GIY-YIG nuclease domain that distinguishes it from the rest of this family. How MSH1 hampers recombination and promotes fidelity in organellar DNA inheritance is unknown. Here, we elucidate its enzymatic activities by recombinantly expressing and purifying full-length MSH1 from Arabidopsis thaliana (AtMSH1). AtMSH1 is a metalloenzyme that shows a strong binding affinity for displacement loops (D-loops). The DNA-binding abilities of AtMSH1 reside in its MutS domain and not in its GIY-YIG domain, which is the ancillary nickase of AtMSH1. In the presence of divalent metal ions, AtMSH1 selectively executes multiple incisions at D-loops, but not other DNA structures including Holliday junctions or dsDNA, regardless of the presence or absence of mismatches. The selectivity of AtMSH1 to dismantle D-loops supports the role of this enzyme in preventing recombination between short repeats. Our results suggest that plant organelles have evolved novel DNA repair routes centered around the anti-recombinogenic activity of MSH1.

Collective mitochondrial dynamics resolve conflicting cellular tensions: From plants to general principles

Mitochondria play diverse and essential roles in eukaryotic cells, and plants are no exception. Plant mitochondria have several differences from their metazoan and fungal cousins: they often exist in a fragmented state, move rapidly on actin rather than microtubules, have many plant-specific metabolic features and roles, and usually contain only a subset of the complete mtDNA genome, which itself undergoes frequent recombination. This arrangement means that exchange and complementation is essential for plant mitochondria, and recent work has begun to reveal how their collective dynamics and resultant "social networks" of encounters support this exchange, connecting plant mitochondria in time rather than in space. This review will argue that this social network perspective can be extended to a "societal network", where mitochondrial dynamics are an essential part of the interacting cellular society of organelles and biomolecules. Evidence is emerging that mitochondrial dynamics allow optimal resolutions to competing cellular priorities; we will survey this evidence and review potential future research directions, highlighting that plant mitochondria can help reveal and test principles that apply across other kingdoms of life. In parallel with this fundamental cell biology, we also highlight the translational "One Health" importance of plant mitochondrial behaviour - which is exploited in the production of a vast amount of crops consumed worldwide - and the potential for multi-objective optimisation to understand and rationally re-engineer the evolved resolutions to these tensions.

Genome copy number predicts extreme evolutionary rate variation in plant mitochondrial DNA

Nuclear and organellar genomes can evolve at vastly different rates despite occupying the same cell. In most bilaterian animals, mitochondrial DNA (mtDNA) evolves faster than nuclear DNA, whereas this trend is generally reversed in plants. However, in some exceptional angiosperm clades, mtDNA substitution rates have increased up to 5,000-fold compared with closely related lineages. The mechanisms responsible for this acceleration are generally unknown. Because plants rely on homologous recombination to repair mtDNA damage, we hypothesized that mtDNA copy numbers may predict evolutionary rates, as lower copy numbers may provide fewer templates for such repair mechanisms. In support of this hypothesis, we found that copy number explains 47% of the variation in synonymous substitution rates of mtDNA across 60 diverse seed plant species representing ~300 million years of evolution. Copy number was also negatively correlated with mitogenome size, which may be a cause or consequence of mutation rate variation. Both relationships were unique to mtDNA and not observed in plastid DNA. These results suggest that homologous recombinational repair plays a role in driving mtDNA substitution rates in plants and may explain variation in mtDNA evolution more broadly across eukaryotes. Our findings also contribute to broader questions about the relationships between mutation rates, genome size, selection efficiency, and the drift-barrier hypothesis.

Stochastic organelle genome segregation through Arabidopsis development and reproduction

Organelle DNA (oDNA) in mitochondria and plastids is vital for plant (and eukaryotic) life. Selection against damaged oDNA is mediated in part by segregation - sorting different oDNA types into different cells in the germline. Plants segregate oDNA very rapidly, with oDNA recombination protein MSH1 a key driver of this segregation, but we have limited knowledge of the dynamics of this segregation within plants and between generations. Here, we reveal how oDNA evolves through Arabidopsis thaliana development and reproduction. We combine stochastic modelling, Bayesian inference, and model selection with new and existing tissue-specific oDNA measurements from heteroplasmic Arabidopsis plant lines through development and between generations. Segregation proceeds gradually but continually during plant development, with a more rapid increase between inflorescence formation and the next generation. When MSH1 is compromised, the majority of observed segregation can be achieved through partitioning at cell divisions. When MSH1 is functional, mtDNA segregation is far more rapid; we show that increased oDNA gene conversion is a plausible mechanism quantitatively explaining this acceleration. These findings reveal the quantitative, time-dependent details of oDNA segregation in Arabidopsis. We also discuss the support for different models of the plant germline provided by these observations.

Morphological characteristics and genetic diversity of floating and attached Ulva prolifera--A case study in the Yellow Sea, China

Green tides in the Yellow Sea have occurred periodically since 2007, impacting the ecological environment and green algal communities along the coasts of Jiangsu and Shandong provinces in China. To investigate the morphological characteristics and genetic diversity of Ulva prolifera, we conducted surveys and comparative analyses of both its floating and attached forms along the coastal areas of Jiangsu and Shandong. The results revealed that the external morphology of the floating U. prolifera was multibranched. The attached U. prolifera displayed significant morphological variation among individuals. Based on the analysis of the amplified characteristic bands of the chloroplast gene, it was shown that both floating and attached U. prolifera could hybridize with Ulva linza. The genetic diversity of U. prolifera was studied using mitochondrial and chloroplast genome fragments. All floating U. prolifera and three strains of attached U. prolifera belonged to the same haplotype. The genetic diversity of floating U. prolifera was low, and there were some genetic differences with attached U. prolifera. The attached U. prolifera displayed a higher level of genetic diversity with abundant sites of variation and haplotypes.

Nuclear DNA segments homologous to mitochondrial DNA are obstacles for detecting heteroplasmy in sugar beet (Beta vulgaris L.)

Heteroplasmy, the coexistence of multiple mitochondrial DNA (mtDNA) sequences in a cell, is well documented in plants. Next-generation sequencing technology (NGS) has made it feasible to sequence entire genomes. Thus, NGS has the potential to detect heteroplasmy; however, the methods and pitfalls in heteroplasmy detection have not been fully investigated and identified. One obstacle for heteroplasmy detection is the sequence homology between mitochondrial-, plastid-, and nuclear DNA, of which the influence of nuclear DNA segments homologous to mtDNA (numt) need to be minimized. To detect heteroplasmy, we first excluded nuclear DNA sequences of sugar beet (Beta vulgaris) line EL10 from the sugar beet mtDNA sequence. NGS reads were obtained from single plants of sugar beet lines NK-195BRmm-O and NK-291BRmm-O and mapped to the unexcluded mtDNA regions. More than 1000 sites exhibited intra-individual polymorphism as detected by genome browsing analysis. We focused on a 309-bp region where 12 intra-individual polymorphic sites were closely linked to each other. Although the existence of DNA molecules having variant alleles at the 12 sites was confirmed by PCR amplification from NK-195BRmm-O and NK-291BRmm-O, these variants were not always called by six variant-calling programs, suggesting that these programs are inappropriate for intra-individual polymorphism detection. When we changed the nuclear DNA reference, a numt absent from EL10 was found to include the 309-bp region. Genetic segregation of an F2 population from NK-195BRmm-O x NK-291BRmm-O supported the numt origin of the variant alleles. Using four references, we found that numt detection exhibited reference dependency, and extreme polymorphism of numts exists among sugar beet lines. One of the identified numts absent from EL10 is also associated with another intra-individual polymorphic site in NK-195mm-O. Our data suggest that polymorphism among numts is unexpectedly high within sugar beets, leading to confusion about the true degree of heteroplasmy.

Fragaria mitogenomes evolve rapidly in structure but slowly in sequence and incur frequent multinucleotide mutations mediated by microinversions

Plant mitochondrial DNA has been described as evolving rapidly in structure but slowly in sequence. However, many of the noncoding portions of plant mitogenomes are not homologous among species, raising questions about the rate and spectrum of mutations in noncoding regions. Recent studies have suggested that the lack of homology in noncoding regions could be due to increased sequence divergence. We compared 30 kb of coding and 200 kb of noncoding DNA from 13 sequenced Fragaria mitogenomes, followed by analysis of the rate of sequence divergence, microinversion events and structural variations. Substitution rates in synonymous sites and nongenic sites are nearly identical, suggesting that the genome-wide point mutation rate is generally consistent. A surprisingly high number of large multinucleotide substitutions were detected in Fragaria mitogenomes, which may have resulted from microinversion events and could affect phylogenetic signal and local rate estimates. Fragaria mitogenomes preferentially accumulate deletions relative to insertions and substantial genomic arrangements, whereas mutation rates could positively associate with these sequence and structural changes among species. Together, these observations suggest that plant mitogenomes exhibit low point mutations genome-wide but exceptionally high structural variations, and our results favour a gain-and-loss model for the rapid loss of homology among plant mitogenomes.

Altered collective mitochondrial dynamics in the Arabidopsis msh1 mutant compromising organelle DNA maintenance

Mitochondria form highly dynamic populations in the cells of plants (and almost all eukaryotes). The characteristics and benefits of this collective behaviour, and how it is influenced by nuclear features, remain to be fully elucidated. Here, we use a recently developed quantitative approach to reveal and analyse the physical and collective 'social' dynamics of mitochondria in an Arabidopsis msh1 mutant where the organelle DNA maintenance machinery is compromised. We use a newly created line combining the msh1 mutant with mitochondrially targeted green fluorescent protein (GFP), and characterize mitochondrial dynamics with a combination of single-cell time-lapse microscopy, computational tracking, and network analysis. The collective physical behaviour of msh1 mitochondria is altered from that of the wild type in several ways: mitochondria become less evenly spread, and networks of inter-mitochondrial encounters become more connected, with greater potential efficiency for inter-organelle exchange-reflecting a potential compensatory mechanism for the genetic challenge to the mitochondrial DNA population, supporting more inter-organelle exchange. We find that these changes are similar to those observed in friendly, where mitochondrial dynamics are altered by a physical perturbation, suggesting that this shift to higher connectivity may reflect a general response to mitochondrial challenges, where physical dynamics of mitochondria may be altered to control the genetic structure of the mtDNA population.

Social networks in the single cell

This article comments on: Chustecki JM, Etherington RD, Gibbs DJ, Johnston IG. 2022. Altered collective mitochondrial dynamics in the Arabidopsis msh1 mutant compromising organelle DNA maintenance. Journal of Experimental Botany 73, 5428-5439. Plant mitochondrial DNA (mtDNA) can become damaged in many ways. A major repair mechanism is homologous recombination, which requires an undamaged DNA template. Presumably, this template comes from a different mitochondrion in the same cell. Plant mitochondria undergo fission and fusion to form transient networks which could allow the exchange of genetic information. To test this hypothesis, Chustecki et al. (2022) used msh1 mutants with defective DNA repair, and showed that mitochondrial interactions increased, revealing a link between the physical and genetic behavior of mitochondria.

Sorting of mitochondrial and plastid heteroplasmy in Arabidopsis is extremely rapid and depends on MSH1 activity

The fate of new mitochondrial and plastid mutations depends on their ability to persist and spread among the numerous organellar genome copies within a cell (heteroplasmy). The extent to which heteroplasmies are transmitted across generations or eliminated through genetic bottlenecks is not well understood in plants, in part because their low mutation rates make these variants so infrequent. Disruption of MutS Homolog 1 (MSH1), a gene involved in plant organellar DNA repair, results in numerous de novo point mutations, which we used to quantitatively track the inheritance of single nucleotide variants in mitochondrial and plastid genomes in Arabidopsis. We found that heteroplasmic sorting (the fixation or loss of a variant) was rapid for both organelles, greatly exceeding rates observed in animals. In msh1 mutants, plastid variants sorted faster than those in mitochondria and were typically fixed or lost within a single generation. Effective transmission bottleneck sizes (N) for plastids and mitochondria were N ∼ 1 and 4, respectively. Restoring MSH1 function further increased the rate of heteroplasmic sorting in mitochondria (N ∼ 1.3), potentially because of its hypothesized role in promoting gene conversion as a mechanism of DNA repair, which is expected to homogenize genome copies within a cell. Heteroplasmic sorting also favored GC base pairs. Therefore, recombinational repair and gene conversion in plant organellar genomes can potentially accelerate the elimination of heteroplasmies and bias the outcome of this sorting process.

The Roles of Mutation and Selection Acting on Mitochondrial Genomes Inferred from Intraspecific Variation in Seed Plants

There is a paradox in the plant mitochondrial genome, that is, the genic region evolves slowly while the intergenic region evolves rapidly. Thus, the intergenic regions of the plant mitochondrial genome are difficult to align across different species, even in closely related species. Here, to character the mechanism of this paradox, we identified interspecific variations in the Ginkgo biloba, Oryza sativa, and Arabidopsis thaliana mitochondrial and plastid genome at a genome-wide level. The substitution rate of synonymous sites in genic regions was similar to the substitution rate of intergenic regions, while the substitution rate of nonsynonymous sites in genic regions was lower than that in intergenic regions, suggesting the mutation inputs were the same among different categories within the organelle genome, but the selection pressure varied. The substitution rate of single-copy regions was higher than that of IR (inverted repeats) in the plastid genome at an intraspecific level. The substitution rate of single-copy regions was higher than that of repeats in the G. biloba and A. thaliana mitochondrial genomes, but lower in that of O. sativa. This difference may be related to the length and distribution of repeats. Copy number variations that existed in the G. biloba and O. sativa mitochondrial genomes were confirmed. This study reveals the intraspecific variation pattern of organelle genomes at a genome-wide level, and that copy number variations were common in plant mitochondrial genomes.

Did doubly uniparental inheritance (DUI) of mtDNA originate as a cytoplasmic male sterility (CMS) system?

Animal and plant species exhibit an astonishing diversity of sexual systems, including environmental and genetic determinants of sex, with the latter including genetic material in the mitochondrial genome. In several hermaphroditic plants for example, sex is determined by an interaction between mitochondrial cytoplasmic male sterility (CMS) genes and nuclear restorer genes. Specifically, CMS involves aberrant mitochondrial genes that prevent pollen development and specific nuclear genes that restore it, leading to a mixture of female (male-sterile) and hermaphroditic individuals in the population (gynodioecy). Such a mitochondrial-nuclear sex determination system is thought to be rare outside plants. Here, we present one possible case of CMS in animals. We hypothesize that the only exception to the strict maternal mtDNA inheritance in animals, the doubly uniparental inheritance (DUI) system in bivalves, might have originated as a mitochondrial-nuclear sex-determination system. We document and explore similarities that exist between DUI and CMS, and we propose various ways to test our hypothesis.

Unprecedented organelle genomic variations in morning glories reveal independent evolutionary scenarios of parasitic plants and the diversification of plant mitochondrial complexes

Background

The morning glories (Convolvulaceae) are distributed worldwide and produce economically important crops, medicinal herbs, and ornamentals. Members of this family are diverse in morphological characteristics and trophic modes, including the leafless parasitic Cuscuta (dodders). Organelle genomes were generally used for studying plant phylogeny and genomic variations. Notably, plastomes in parasitic plants always show non-canonical features, such as reduced size and accelerated rates. However, few organelle genomes of this group have been sequenced, hindering our understanding of their evolution, and dodder mitogenome in particular.

Results

We assembled 22 new mitogenomes and 12 new plastomes in Convolvulaceae. Alongside previously known ones, we totally analyzed organelle genomes of 23 species in the family. Our sampling includes 16 leafy autotrophic species and 7 leafless parasitic dodders, covering 8 of the 12 tribes. Both the plastid and mitochondrial genomes of these plants have encountered variations that were rarely observed in other angiosperms. All of the plastomes possessed atypical IR boundaries. Besides the gene and IR losses in dodders, some leafy species also showed gene and intron losses, duplications, structural variations, and insertions of foreign DNAs. The phylogeny reconstructed by plastid protein coding sequences confirmed the previous relationship of the tribes. However, the monophyly of ‘Merremieae’ and the sister group of Cuscuta remained uncertain. The mitogenome was significantly inflated in Cuscuta japonica, which has exceeded over 800 kb and integrated massive DNAs from other species. In other dodders, mitogenomes were maintained in small size, revealing divergent evolutionary strategies. Mutations unique to plants were detected in the mitochondrial gene ccmFc, which has broken into three fragments through gene fission and splicing shift. The unusual changes likely initially happened to the common ancestor of the family and were caused by a foreign insertion from rosids followed by double-strand breaks and imprecise DNA repairs. The coding regions of ccmFc expanded at both sides after the fission, which may have altered the protein structure.

Conclusions

Our family-scale analyses uncovered unusual scenarios for both organelle genomes in Convolvulaceae, especially in parasitic plants. The data provided valuable genetic resources for studying the evolution of Convolvulaceae and plant parasitism.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01250-1.

Chloroplast Genome Traits Correlate With Organismal Complexity and Ecological Traits in Chlorophyta

A positive relationship between cell size and chloroplast genome size within chloroplast-bearing protists has been hypothesized in the past and shown in some case studies, but other factors influencing chloroplast genome size during the evolution of chlorophyte algae have been less studied. We study chloroplast genome size and GC content as a function of habitats and cell size of chlorophyte algae. The chloroplast genome size of green algae in freshwater, marine and terrestrial habitats was differed significantly, with terrestrial algae having larger chloroplast genome sizes in general. The most important contributor to these enlarged genomes in terrestrial species was the length of intergenic regions. There was no clear difference in the GC content of chloroplast genomes from the three habitats categories. Functional morphological categories also showed differences in chloroplast genome size, with filamentous algae having substantially larger genomes than other forms of algae, and foliose algae had lower GC content than other groups. Chloroplast genome size showed no significant differences among the classes Ulvophyceae, Trebouxiophyceae, and Chlorophyceae, but the GC content of Chlorophyceae chloroplast genomes was significantly lower than that of Ulvophyceae and Trebouxiophyceae. There was a certain positive relationship between chloroplast genome size and cell size for the Chlorophyta as a whole and within each of three major classes. Our data also confirmed previous reports that ancestral quadripartite architecture had been lost many times independently in Chlorophyta. Finally, the comparison of the phenotype of chlorophytes algae harboring plastids uncovered that most of the investigated Chlorophyta algae housed a single plastid per cell.

Organelle Genome Variation in the Red Algal Genus Ahnfeltia (Florideophyceae)

The agarophyte Ahnfeltia (Ahnfeltiales, Rhodophyta) is a globally widespread genus with 11 accepted species names. Two of the most widespread species in this genus, A. plicata and A. fastigiata, may have diverged genetically due to past geographic changes and subsequent geographic isolation. To investigate this genomic and genetic diversity, we generated new plastid (ptDNAs) and mitochondrial genomes (mtDNAs) of these Ahnfeltia species from four different regions (A. plicata - Chile and UK and A. fastigiata - Korea and Oregon). Two architecture variations were found in the Ahnfeltia genomes: in ptDNA of A. fastigiata Oregon, the hypothetical pseudogene region was translocated, likely due to recombination with palindromic repeats or a gene transfer from a red algal plasmid. In mtDNA of A. fastigiata Korea, the composition of the group II intronic ORFs was distinct from others suggesting different scenarios of gain and loss of group II intronic ORFs. These features resulted in genome size differences between the two species. Overall gene contents of organelle genomes of Ahnfeltia were conserved. Phylogenetic analysis using concatenated genes from ptDNAs and mtDNAs supported the monophyly of the Ahnfeltiophycidae. The most probable individual gene trees showed that the Ahnfeltia populations were genetically diversified. These trees, the cox1 haplotype network, and a dN/dS analysis all supported the theory that these Ahnfeltia populations have diversified genetically in accordance with geographic distribution.

Microhomologies Are Associated with Tandem Duplications and Structural Variation in Plant Mitochondrial Genomes

Short tandem repeats (STRs) contribute to structural variation in plant mitochondrial genomes, but the mechanisms underlying their formation and expansion are unclear. In this study, we detected high polymorphism in the nad7-1 region of the Pinus tabuliformis mitogenome caused by the rapid accumulation of STRs and rearrangements over a few million years ago. The STRs in nad7-1 have a 7-bp microhomology (TAG7) flanking the repeat array. We then scanned the mitogenomes of 136 seed plants to understand the role of microhomology in the formation of STR and mitogenome evolution. A total of 13,170 STRs were identified, and almost half of them were associated with microhomologies. A substantial amount (1,197) of microhomologies was long enough to mediate structural variation, and the length of microhomology is positively correlated with the length of tandem repeat unit. These results suggest that microhomology may be involved in the formation of tandem repeat via microhomology-mediated pathway, and the formation of longer duplicates required greater length of microhomology. We examined the abundance of these 1,197 microhomologies, and found 75% of them were enriched in the plant mitogenomes. Further analyses of the 400 prevalent microhomologies revealed that 175 of them showed differential enrichment between angiosperms and gymnosperms and 186 differed between angiosperms and conifers, indicating lineage-specific usage and expansion of microhomologies. Our study sheds light on the sources of structural variation in plant mitochondrial genomes and highlights the importance of microhomology in mitochondrial genome evolution.

Slow crabs - fast genomes: Locomotory capacity predicts skew magnitude in crustacean mitogenomes

Base composition skews (G-C/G+C) of mitochondrial genomes are believed to be primarily driven by mutational pressure, which is positively correlated with metabolic rate. In marine animals, metabolic rate is also positively correlated with locomotory capacity. Given the central role of mitochondria in energy metabolism, we hypothesised that selection for locomotory capacity should be positively correlated with the strength of purifying selection (dN/dS), and thus be negatively correlated with the skew magnitude. Therefore, these two models assume diametrically opposite associations between the metabolic rate and skew magnitude: positive correlation in the prevailing paradigm, and negative in our working hypothesis. We examined correlations between the skew magnitude, metabolic rate, locomotory capacity, and several other variables previously associated with mitochondrial evolution on 287 crustacean mitogenomes. Weakly locomotory taxa had higher skew magnitude and ω (dN/dS) values, but not the gene order rearrangement rate. Skew and ω magnitudes were correlated. Multilevel regression analyses indicated that three competing variables, body size, gene order rearrangement rate, and effective population size, had negligible impacts on the skew magnitude. In most crustacean lineages selection for locomotory capacity appears to be the primary factor determining the skew magnitude. Contrary to the prevailing paradigm, this implies that adaptive selection outweighs nonadaptive selection (mutation pressure) in crustaceans. However, we found indications that effective population size (nonadaptive factor) may outweigh the impact of locomotory capacity in sessile crustaceans (Thecostraca). In conclusion, skew magnitude is a product of the interplay between adaptive and nonadaptive factors, the balance of which varies among lineages.

In and out: Evolution of viral sequences in the mitochondrial genomes of legumes (Fabaceae)

Plant specific mitoviruses in the 'genus' Mitovirus (Narnaviridae) and their integrated sequences (non-retroviral endogenous RNA viral elements or NERVEs) have been recently identified in various plant lineages. However, the sparse phylogenetic coverage of complete plant mitochondrial genome (mitogenome) sequences and the non-conserved nature of mitochondrial intergenic regions have hindered comparative studies on mitovirus NERVEs in plants. In this study, 10 new mitogenomes were sequenced from legumes (Fabaceae). Based on comparative genomic analysis of 27 total mitogenomes, we identified mitovirus NERVEs and transposable elements across the family. All legume mitogenomes included NERVEs and total NERVE length varied from ca. 2 kb in the papilionoid Trifolium to 35 kb in the mimosoid Acacia. Most of the NERVE integration sites were in highly variable intergenic regions, however, some were positioned in six cis-spliced mitochondrial introns. In the Acacia mitogenome, there were L1-like transposon sequences including an almost full-length copy with target site duplications (TSDs). The integration sites of NERVEs in four introns showed evidence of L1-like retrotransposition events. Phylogenetic analysis revealed that there were multiple instances of precise deletion of NERVEs between TSDs. This study provides clear evidence that a L1-like retrotransposition mechanism has a long history of contributing to the integration of viral RNA into plant mitogenomes while microhomology-mediated deletion can restore the integration site.

Detecting de novo mitochondrial mutations in angiosperms with highly divergent evolutionary rates

Although plant mitochondrial genomes typically show low rates of sequence evolution, levels of divergence in certain angiosperm lineages suggest anomalously high mitochondrial mutation rates. However, de novo mutations have never been directly analyzed in such lineages. Recent advances in high-fidelity DNA sequencing technologies have enabled detection of mitochondrial mutations when still present at low heteroplasmic frequencies. To date, these approaches have only been performed on a single plant species (Arabidopsis thaliana). Here, we apply a high-fidelity technique (Duplex Sequencing) to multiple angiosperms from the genus Silene, which exhibits extreme heterogeneity in rates of mitochondrial sequence evolution among close relatives. Consistent with phylogenetic evidence, we found that Silene latifolia maintains low mitochondrial variant frequencies that are comparable with previous measurements in Arabidopsis. Silene noctiflora also exhibited low variant frequencies despite high levels of historical sequence divergence, which supports other lines of evidence that this species has reverted to lower mitochondrial mutation rates after a past episode of acceleration. In contrast, S. conica showed much higher variant frequencies in mitochondrial (but not in plastid) DNA, consistent with an ongoing bout of elevated mitochondrial mutation rates. Moreover, we found an altered mutational spectrum in S. conica heavily biased towards AT→GC transitions. We also observed an unusually low number of mitochondrial genome copies per cell in S. conica, potentially pointing to reduced opportunities for homologous recombination to accurately repair mismatches in this species. Overall, these results suggest that historical fluctuations in mutation rates are driving extreme variation in rates of plant mitochondrial sequence evolution.

Plant Mitochondria are a Riddle Wrapped in a Mystery Inside an Enigma

A fundamental paradox motivates the study of plant mitochondrial genomics: the mutation rate is very low (lower than in the nucleus) but the rearrangement rate is high. A landmark paper published in Journal of Molecular Evolution in 1988 established these facts and revealed the paradox. Jeffrey Palmer and Laura Herbon did a prodigious amount of work in the pre-genome sequencing era to identify both the high frequency of rearrangements between closely related species, and the low frequency of mutations, observations that have now been confirmed many times by sequencing. This paper was also the first to use molecular data on rearrangements as a phylogenetic trait to build a parsimonious tree. The work was a technical tour-de-force, its findings are still at the heart of plant mitochondrial genomics, and the underlying molecular mechanisms that produce this paradox are still not completely understood.

Six Newly Sequenced Chloroplast Genomes From Trentepohliales: The Inflated Genomes, Alternative Genetic Code and Dynamic Evolution

Cephaleuros is often known as an algal pathogen with 19 taxonomically valid species, some of which are responsible for red rust and algal spot diseases in vascular plants. No chloroplast genomes have yet been reported in this genus, and the limited genetic information is an obstacle to understanding the evolution of this genus. In this study, we sequenced six new Trentepohliales chloroplast genomes, including four Cephaleuros and two Trentepohlia. The chloroplast genomes of Trentepohliales are large compared to most green algae, ranging from 216 to 408 kbp. They encode between 93 and 98 genes and have a GC content of 26-36%. All new chloroplast genomes were circular-mapping and lacked a quadripartite structure, in contrast to the previously sequenced Trentepohlia odorata, which does have an inverted repeat. The duplicated trnD -GTC, petD, and atpA genes in C. karstenii may be remnants of the IR region and shed light on its reduction. Chloroplast genes of Trentepohliales show elevated rates of evolution, strong rearrangement dynamics and several genes display an alternative genetic code with reassignment of the UGA/UAG codon presumably coding for arginine. Our results present the first whole chloroplast genome of the genus Cephaleuros and enrich the chloroplast genome resources of Trentepohliales.

The development of a novel diagnostic PCR for Madurella mycetomatis using a comparative genome approach

Background

Eumycetoma is a neglected tropical disease most commonly caused by the fungus Madurella mycetomatis. Identification of eumycetoma causative agents can only be reliably performed by molecular identification, most commonly by species-specific PCR. The current M. mycetomatis specific PCR primers were recently discovered to cross-react with Madurella pseudomycetomatis. Here, we used a comparative genome approach to develop a new M. mycetomatis specific PCR for species identification.

Methodology

Predicted-protein coding sequences unique to M. mycetomatis were first identified in BLASTCLUST based on E-value, size and presence of orthologues. Primers were then developed for 16 unique sequences and evaluated against 60 M. mycetomatis isolates and other eumycetoma causing agents including the Madurella sibling species. Out of the 16, only one was found to be specific to M. mycetomatis.

Conclusion

We have discovered a predicted-protein coding sequence unique to M. mycetomatis and have developed a new species-specific PCR to be used as a novel diagnostic marker for M. mycetomatis.

Mycetoma is a neglected tropical disease characterised by tumorous swellings and grain formation. This disease can be caused by more than 70 different micro-organisms and is categorised into actinomycetoma (caused by bacteria) and eumycetoma (caused by fungi). The most common causative agent of mycetoma is the fungus Madurella mycetomatis. Diagnosis of eumycetoma is often only done clinically or by histopathological examination and culturing of the grains. Unfortunately, that often leads to misidentifications. Molecular identification is currently the most reliable method to identify the causative agents. However, we have recently discovered that the only M. mycetomatis species-specific PCR primers cross-reacts to Madurella pseudomycetomatis. Since all Madurella species cause eumycetoma and have different susceptibilities to antifungal agents, it is important to be able to accurately identify them to the species level. Here we have used a comparative genome approach to identify and design new M. mycetomatis species-specific PCR primers. These primers can be used to identify M. mycetomatis directly from grains and do not cross-react with any of the other eumycetoma causative agents tested. We, therefore, recommended the use of these primers in reference centres and local laboratories to identify M. mycetomatis to the species level.

Can Green Algal Plastid Genome Size Be Explained by DNA Repair Mechanisms?

A major finding in organelle biology over the past decade is that land plant mitochondrial genomes, which are the largest among eukaryotes, can have a “Jekyll and Hyde” mutational pattern: low for synonymous sites, high for intergenic ones. This has led to the theory that double-strand breaks (DSBs) in the intergenic DNA of plant mitogenomes are repaired by inaccurate mechanisms, such as break-induced replication, which can result in large insertions and, thus, could explain why these genomes are so prone to expansion. But how universal is this theory? Can it apply to other giant organelle DNAs, such as the massive plastid DNAs (ptDNAs) of chlamydomonadalean green algae? Indeed, it can. Analysis of the expanded plastomes from two distinct isolates of the unicellular chlamydomonadalean Chlorosarcinopsis eremi uncovered exceptionally low rates of synonymous substitution in the coding regions but high substitution rates, including frequent indels, in the noncoding ptDNA, mirroring the trend from land plant mitogenomes. Remarkably, nearly all of the substitutions and indels identified in the noncoding ptDNA of C. eremi occur adjacent to or within short inverted palindromic repeats, suggesting that these elements are mutational hotspots. Building upon earlier studies, I propose that these palindromic repeats are predisposed to DSBs and that error-prone repair of these breaks is contributing to genomic expansion. Short palindromic repeats are a common theme among bloated plastomes, including the largest one on record, meaning that these data could have wide-reaching implications for our understanding of ptDNA expansion.

MSH1 is required for maintenance of the low mutation rates in plant mitochondrial and plastid genomes

Mitochondrial and plastid genomes in land plants exhibit some of the slowest rates of sequence evolution observed in any eukaryotic genome, suggesting an exceptional ability to prevent or correct mutations. However, the mechanisms responsible for this extreme fidelity remain unclear. We tested seven candidate genes involved in cytoplasmic DNA replication, recombination, and repair (POLIA, POLIB, MSH1, RECA3, UNG, FPG, and OGG1) for effects on mutation rates in the model angiosperm Arabidopsis thaliana by applying a highly accurate DNA sequencing technique (duplex sequencing) that can detect newly arisen mitochondrial and plastid mutations even at low heteroplasmic frequencies. We find that disrupting MSH1 (but not the other candidate genes) leads to massive increases in the frequency of point mutations and small indels and changes to the mutation spectrum in mitochondrial and plastid DNA. We also used droplet digital PCR to show transmission of de novo heteroplasmies across generations in msh1 mutants, confirming a contribution to heritable mutation rates. This dual-targeted gene is part of an enigmatic lineage within the mutS mismatch repair family that we find is also present outside of green plants in multiple eukaryotic groups (stramenopiles, alveolates, haptophytes, and cryptomonads), as well as certain bacteria and viruses. MSH1 has previously been shown to limit ectopic recombination in plant cytoplasmic genomes. Our results point to a broader role in recognition and correction of errors in plant mitochondrial and plastid DNA sequence, leading to greatly suppressed mutation rates perhaps via initiation of double-stranded breaks and repair pathways based on faithful homologous recombination.

MSH1 is required for maintenance of the low mutation rates in plant mitochondrial and plastid genomes

Mitochondrial and plastid genomes in land plants exhibit some of the slowest rates of sequence evolution observed in any eukaryotic genome, suggesting an exceptional ability to prevent or correct mutations. However, the mechanisms responsible for this extreme fidelity remain unclear. We tested seven candidate genes involved in cytoplasmic DNA replication, recombination, and repair (POLIA, POLIB, MSH1, RECA3, UNG, FPG, and OGG1) for effects on mutation rates in the model angiosperm Arabidopsis thaliana by applying a highly accurate DNA sequencing technique (duplex sequencing) that can detect newly arisen mitochondrial and plastid mutations even at low heteroplasmic frequencies. We find that disrupting MSH1 (but not the other candidate genes) leads to massive increases in the frequency of point mutations and small indels and changes to the mutation spectrum in mitochondrial and plastid DNA. We also used droplet digital PCR to show transmission of de novo heteroplasmies across generations in msh1 mutants, confirming a contribution to heritable mutation rates. This dual-targeted gene is part of an enigmatic lineage within the mutS mismatch repair family that we find is also present outside of green plants in multiple eukaryotic groups (stramenopiles, alveolates, haptophytes, and cryptomonads), as well as certain bacteria and viruses. MSH1 has previously been shown to limit ectopic recombination in plant cytoplasmic genomes. Our results point to a broader role in recognition and correction of errors in plant mitochondrial and plastid DNA sequence, leading to greatly suppressed mutation rates perhaps via initiation of double-stranded breaks and repair pathways based on faithful homologous recombination.

The Mitochondrial Genome of Eleusine indica and Characterization of Gene Content within Poaceae

Plant mitochondrial (mt) genome assembly provides baseline data on size, structure, and gene content, but resolving the sequence of these large and complex organelle genomes remains challenging due to fragmentation, frequent recombination, and transfers of DNA from neighboring plastids. The mt genome for Eleusine indica (Poaceae: goosegrass) is comprehensibly analyzed here, providing key reference data for an economically significant invasive species that is also the maternal parent of the allotetraploid crop Finger millet (Eleusine coracana). The assembled E. indica genome contains 33 protein coding genes, 6 rRNA subunits, 24 tRNA, 8 large repetitive regions 15 kb of transposable elements across a total of 520,691 bp. Evidence of RNA editing and loss of rpl2, rpl5, rps14, rps11, sdh4, and sdh3 genes is evaluated in the context of an updated survey of mt genomic gene content across the grasses through an analysis of publicly available data. Hypothesized patterns of Poaceae mt gene loss are examined in a phylogenetic context to clarify timing, showing that rpl2 was transferred to the nucleus from the mitochondrion prior to the origin of the PACMAD clade.

Comparative Analysis of the Mitochondrial Genomes of Nicotiana tabacum: Hints Toward the Key Factors Closely Related to the Cytoplasmic Male Sterility Mechanism

Background

Cytoplasmic male sterility (CMS) is a complex phenomenon of plant sterility that can produce non-functional pollen. It is caused by mutation, rearrangement or recombination in the mitochondrial genome. So far, the systematic structural characteristics of the changes in the mitochondrial genome from the maintainer lines to the CMS lines have not been reported in tobacco.

Results

The mitochondrial genomes of the flower buds from both CMS lines and maintainer lines of two Nicotiana tabacum cultivars (YY85, sYY85, ZY90, and sZY90) were sequenced using the PacBio and Illumina Hiseq technology, and several findings were produced by comparative analysis based on the de novo sequencing. (1) The genomes of the CMS lines were larger, and the different areas were mostly non-coding regions. (2) A large number of rearrangement regions were detected in the CMS lines, with many translocation regions. (3) Thirteen gene clusters were shared by the four mitochondrial genomes, among which two of the gene clusters, nad2-sdh3 and nad6-rps4, were far from each other in the CMS lines. (4) Thirty-three protein-coding genes were conserved in four mitochondrial genomes. However, nad3 was detected one additional copy in the maintainer lines, and sequence differences were revealed in the four candidate genes (atp6, cox2, nad2, and sdh3). Importantly, the evolutionary tree based on the four genes could be used to distinguish the CMS lines and the maintainer lines well for the sequenced mitochondrial genomes of the tobacco. (5) Sixteen CMS-specific open reading frames (ORFs) were found, three of which (orf91, orf115b, and orf100) were previously reported. (6) The differences in intensity of the protein–protein (PPI) interaction in ATP6 were further verified using the yeast two-hybrid analysis.

Conclusion

Although the majority of the sequences, genes and gene clusters were shared by the mitochondrial genomes of the maintainer and the CMS lines in tobacco, extensive structural variations identified with comprehensive analysis based on the mitochondrial genomes, including rearrangement, gene order, the mitochondrial genome expansion and shrinkage events, might be related to CMS. Additionally, the candidate protein-coding genes and CMS-specific ORFs were closely associated with the CMS mechanism. Verification experiments of one of the candidate genes were performed, and the validity of our research results was supported.

Comparative Mitogenome Analysis of the Genus Trifolium Reveals Independent Gene Fission of ccmFn and Intracellular Gene Transfers in Fabaceae

The genus Trifolium is the largest of the tribe Trifolieae in the subfamily Papilionoideae (Fabaceae). The paucity of mitochondrial genome (mitogenome) sequences has hindered comparative analyses among the three genomic compartments of the plant cell (nucleus, mitochondrion and plastid). We assembled four mitogenomes from the two subgenera (Chronosemium and Trifolium) of the genus. The four Trifolium mitogenomes were compact (294,911-348,724 bp in length) and contained limited repetitive (6.6-8.6%) DNA. Comparison of organelle repeat content highlighted the distinct evolutionary trajectory of plastid genomes in a subset of Trifolium species. Intracellular gene transfer (IGT) was analyzed among the three genomic compartments revealing functional transfer of mitochondrial rps1 to nuclear genome along with other IGT events. Phylogenetic analysis based on mitochondrial and nuclear rps1 sequences revealed that the functional transfer in Trifolieae was independent from the event that occurred in robinioid clade that includes genus Lotus. A novel, independent fission event of ccmFn in Trifolium was identified, caused by a 59 bp deletion. Fissions of this gene reported previously in land plants were reassessed and compared with Trifolium.

The Tempo and Mode of Angiosperm Mitochondrial Genome Divergence Inferred from Intraspecific Variation in Arabidopsis thaliana

The mechanisms of sequence divergence in angiosperm mitochondrial genomes have long been enigmatic. In particular, it is difficult to reconcile the rapid divergence of intergenic regions that can make non-coding sequences almost unrecognizable even among close relatives with the unusually high levels of sequence conservation found in genic regions. It has been hypothesized that different mutation and repair mechanisms act on genic and intergenic sequences or alternatively that mutational input is relatively constant but that selection has strikingly different effects on these respective regions. To test these alternative possibilities, we analyzed mtDNA divergence within Arabidopsis thaliana, including variants from the 1001 Genomes Project and changes accrued in published mutation accumulation (MA) lines. We found that base-substitution frequencies are relatively similar for intergenic regions and synonymous sites in coding regions, whereas indel and nonsynonymous substitutions rates are greatly depressed in coding regions, supporting a conventional model in which mutation/repair mechanisms are consistent throughout the genome but differentially filtered by selection. Most types of sequence and structural changes were undetectable in 10-generation MA lines, but we found significant shifts in relative copy number across mtDNA regions for lines grown under stressed vs. benign conditions. We confirmed quantitative variation in copy number across the A. thaliana mitogenome using both whole-genome sequencing and droplet digital PCR, further undermining the classic but oversimplified model of a circular angiosperm mtDNA structure. Our results suggest that copy number variation is one of the most fluid features of angiosperm mitochondrial genomes.

Evolving mtDNA populations within cells

Mitochondrial DNA (mtDNA) encodes vital respiratory machinery. Populations of mtDNA molecules exist in most eukaryotic cells, subject to replication, degradation, mutation, and other population processes. These processes affect the genetic makeup of cellular mtDNA populations, changing cell-to-cell distributions, means, and variances of mutant mtDNA load over time. As mtDNA mutant load has nonlinear effects on cell functionality, and cell functionality has nonlinear effects on tissue performance, these statistics of cellular mtDNA populations play vital roles in health, disease, and inheritance. This mini review will describe some of the better-known ways in which these populations change over time in different organisms, highlighting the importance of quantitatively understanding both mutant load mean and variance. Due to length constraints, we cannot attempt to be comprehensive but hope to provide useful links to some of the many excellent studies on these topics.

Mitochondrial DNA Repair in an Arabidopsis thaliana Uracil N-Glycosylase Mutant

Substitution rates in plant mitochondrial genes are extremely low, indicating strong selective pressure as well as efficient repair. Plant mitochondria possess base excision repair pathways; however, many repair pathways such as nucleotide excision repair and mismatch repair appear to be absent. In the absence of these pathways, many DNA lesions must be repaired by a different mechanism. To test the hypothesis that double-strand break repair (DSBR) is that mechanism, we maintained independent self-crossing lineages of plants deficient in uracil-N-glycosylase (UNG) for 11 generations to determine the repair outcomes when that pathway is missing. Surprisingly, no single nucleotide polymorphisms (SNPs) were fixed in any line in generation 11. The pattern of heteroplasmic SNPs was also unaltered through 11 generations. When the rate of cytosine deamination was increased by mitochondrial expression of the cytosine deaminase APOBEC3G, there was an increase in heteroplasmic SNPs but only in mature leaves. Clearly, DNA maintenance in reproductive meristem mitochondria is very effective in the absence of UNG while mitochondrial genomes in differentiated tissue are maintained through a different mechanism or not at all. Several genes involved in DSBR are upregulated in the absence of UNG, indicating that double-strand break repair is a general system of repair in plant mitochondria. It is important to note that the developmental stage of tissues is critically important for these types of experiments.

The complete mitochondrial genome of Taxus cuspidata (Taxaceae): eight protein-coding genes have transferred to the nuclear genome

BACKGROUND

Gymnosperms represent five of the six lineages of seed plants. However, most sequenced plant mitochondrial genomes (mitogenomes) have been generated for angiosperms, whereas mitogenomic sequences have been generated for only six gymnosperms. In particular, complete mitogenomes are available for all major seed plant lineages except Conifer II (non-Pinaceae conifers or Cupressophyta), an important lineage including six families, which impedes a comprehensive understanding of the mitogenomic diversity and evolution in gymnosperms.

RESULTS

Here, we report the complete mitogenome of Taxus cuspidata in Conifer II. In comparison with previously released gymnosperm mitogenomes, we found that the mitogenomes of Taxus and Welwitschia have lost many genes individually, whereas all genes were identified in the mitogenomes of Cycas, Ginkgo and Pinaceae. Multiple tRNA genes and introns also have been lost in some lineages of gymnosperms, similar to the pattern observed in angiosperms. In general, gene clusters could be less conserved in gymnosperms than in angiosperms. Moreover, fewer RNA editing sites were identified in the Taxus and Welwitschia mitogenomes than in other mitogenomes, which could be correlated with fewer introns and frequent gene losses in these two species.

CONCLUSIONS

We have sequenced the Taxus cuspidata mitogenome, and compared it with mitogenomes from the other four gymnosperm lineages. The results revealed the diversity in size, structure, gene and intron contents, foreign sequences, and mutation rates of gymnosperm mitogenomes, which are different from angiosperm mitogenomes.

DNA Repair and the Stability of the Plant Mitochondrial Genome

The mitochondrion stands at the center of cell energy metabolism. It contains its own genome, the mtDNA, that is a relic of its prokaryotic symbiotic ancestor. In plants, the mitochondrial genetic information influences important agronomic traits including fertility, plant vigor, chloroplast function, and cross-compatibility. Plant mtDNA has remarkable characteristics: It is much larger than the mtDNA of other eukaryotes and evolves very rapidly in structure. This is because of recombination activities that generate alternative mtDNA configurations, an important reservoir of genetic diversity that promotes rapid mtDNA evolution. On the other hand, the high incidence of ectopic recombination leads to mtDNA instability and the expression of gene chimeras, with potential deleterious effects. In contrast to the structural plasticity of the genome, in most plant species the mtDNA coding sequences evolve very slowly, even if the organization of the genome is highly variable. Repair mechanisms are probably responsible for such low mutation rates, in particular repair by homologous recombination. Herein we review some of the characteristics of plant organellar genomes and of the repair pathways found in plant mitochondria. We further discuss how homologous recombination is involved in the evolution of the plant mtDNA.

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.

Structure-Function Analysis Reveals the Singularity of Plant Mitochondrial DNA Replication Components: A Mosaic and Redundant System

Plants are sessile organisms, and their DNA is particularly exposed to damaging agents. The integrity of plant mitochondrial and plastid genomes is necessary for cell survival. During evolution, plants have evolved mechanisms to replicate their mitochondrial genomes while minimizing the effects of DNA damaging agents. The recombinogenic character of plant mitochondrial DNA, absence of defined origins of replication, and its linear structure suggest that mitochondrial DNA replication is achieved by a recombination-dependent replication mechanism. Here, I review the mitochondrial proteins possibly involved in mitochondrial DNA replication from a structural point of view. A revision of these proteins supports the idea that mitochondrial DNA replication could be replicated by several processes. The analysis indicates that DNA replication in plant mitochondria could be achieved by a recombination-dependent replication mechanism, but also by a replisome in which primers are synthesized by three different enzymes: Mitochondrial RNA polymerase, Primase-Helicase, and Primase-Polymerase. The recombination-dependent replication model and primers synthesized by the Primase-Polymerase may be responsible for the presence of genomic rearrangements in plant mitochondria.

Fluctuations in Fabaceae mitochondrial genome size and content are both ancient and recent

BACKGROUND

Organelle genome studies of Fabaceae, an economically and ecologically important plant family, have been biased towards the plastid genome (plastome). Thus far, less than 15 mitochondrial genome (mitogenome) sequences of Fabaceae have been published, all but four of which belong to the subfamily Papilionoideae, limiting the understanding of size variation and content across the family. To address this, four mitogenomes were sequenced and assembled from three different subfamilies (Cercidoideae, Detarioideae and Caesalpinioideae).

RESULTS

Phylogenetic analysis based on shared mitochondrial protein coding regions produced a fully resolved and well-supported phylogeny that was completely congruent with the plastome tree. Comparative analyses suggest that two kinds of mitogenome expansions have occurred in Fabaceae. Size expansion of four genera (Tamarindus, Libidibia, Haematoxylum, and Leucaena) in two subfamilies (Detarioideae and Caesalpinioideae) occurred in relatively deep nodes, and was mainly caused by intercellular gene transfer and/or interspecific horizontal gene transfer (HGT). The second, more recent expansion occurred in the Papilionoideae as a result of duplication of native mitochondrial sequences. Family-wide gene content analysis revealed 11 gene losses, four (rps2, 7, 11 and 13) of which occurred in the ancestor of Fabaceae. Losses of the remaining seven genes (cox2, rpl2, rpl10, rps1, rps19, sdh3, sdh4) were restricted to specific lineages or occurred independently in different clades. Introns of three genes (cox2, ccmFc and rps10) showed extensive lineage-specific length variation due to large sequence insertions and deletions. Shared DNA analysis among Fabaceae mitogenomes demonstrated a substantial decay of intergenic spacers and provided further insight into HGT between the mimosoid clade of Caesalpinioideae and the holoparasitic Lophophytum (Balanophoraceae).

CONCLUSION

This study represents the most exhaustive analysis of Fabaceae mitogenomes so far, and extends the understanding the dynamic variation in size and gene/intron content. The four newly sequenced mitogenomes reported here expands the phylogenetic coverage to four subfamilies. The family has experienced multiple mitogenome size fluctuations in both ancient and recent times. The causes of these size variations are distinct in different lineages. Fabaceae mitogenomes experienced extensive size fluctuation by recruitment of exogenous DNA and duplication of native mitochondrial DNA.

Plant organellar DNA polymerases repair double-stranded breaks by microhomology-mediated end-joining

Double-stranded breaks (DSBs) in plant organelles are repaired via genomic rearrangements characterized by microhomologous repeats. These microhomologous signatures predict the existence of an unidentified enzymatic machinery capable of repairing of DSBs via microhomology-mediated end-joining (MMEJ) in plant organelles. Here, we show that organellar DNA polymerases from Arabidopsis thaliana (AtPolIA and AtPolIB) perform MMEJ using microhomologous sequences as short as six nucleotides. AtPolIs execute MMEJ by virtue of two specialized amino acid insertions located in their thumb subdomains. Single-stranded binding proteins (SSBs) unique to plants, AtWhirly2 and organellar single-stranded binding proteins (AtOSBs), hinder MMEJ, whereas canonical mitochondrial SSBs (AtmtSSB1 and AtmtSSB2) do not interfere with MMEJ. Our data predict that organellar DNA rearrangements by MMEJ are a consequence of a competition for the 3'-OH of a DSBs. If AtWhirlies or AtOSBs gain access to the single-stranded DNA (ssDNA) region of a DSB, the reaction will shift towards high-fidelity routes like homologous recombination. Conversely MMEJ would be favored if AtPolIs or AtmtSSBs interact with the DSB. AtPolIs are not phylogenetically related to metazoan mitochondrial DNA polymerases, and the ability of AtPolIs to execute MMEJ may explain the abundance of DNA rearrangements in plant organelles in comparison to animal mitochondria.

Tracing the Evolution of the Plastome and Mitogenome in the Chloropicophyceae Uncovered Convergent tRNA Gene Losses and a Variant Plastid Genetic Code

The tiny green algae belonging to the Chloropicophyceae play a key role in marine phytoplankton communities; this newly erected class of prasinophytes comprises two genera (Chloropicon and Chloroparvula) containing each several species. We sequenced the plastomes and mitogenomes of eight Chloropicon and five Chloroparvula species to better delineate the phylogenetic affinities of these taxa and to infer the suite of changes that their organelle genomes sustained during evolution. The relationships resolved in organelle-based phylogenomic trees were essentially congruent with previously reported rRNA trees, and similar evolutionary trends but distinct dynamics were identified for the plastome and mitogenome. Although the plastome sustained considerable changes in gene content and order at the time the two genera split, subsequently it remained stable and maintained a very small size. The mitogenome, however, was remodeled more gradually and showed more fluctuation in size, mainly as a result of expansions/contractions of intergenic regions. Remarkably, the plastome and mitogenome lost a common set of three tRNA genes, with the trnI(cau) and trnL(uaa) losses being accompanied with important variations in codon usage. Unexpectedly, despite the disappearance of trnI(cau) from the plastome in the Chloroparvula lineage, AUA codons (the codons recognized by this gene product) were detected in certain plastid genes. By comparing the sequences of plastid protein-coding genes from chloropicophycean and phylogenetically diverse chlorophyte algae with those of the corresponding predicted proteins, we discovered that the AUA codon was reassigned from isoleucine to methionine in Chloroparvula. This noncanonical genetic code has not previously been uncovered in plastids.

High and Variable Rates of Repeat-Mediated Mitochondrial Genome Rearrangement in a Genus of Plants

For 30 years, it has been clear that angiosperm mitochondrial genomes evolve rapidly in sequence arrangement (i.e., synteny), yet absolute rates of rearrangement have not been measured in any plant group, nor is it known how much these rates vary. To investigate these issues, we sequenced and reconstructed the rearrangement history of seven mitochondrial genomes in Monsonia (Geraniaceae). We show that rearrangements (occurring mostly as inversions) not only take place at generally high rates in these genomes but also uncover significant variation in rearrangement rates. For example, the hyperactive mitochondrial genome of Monsonia ciliata has accumulated at least 30 rearrangements over the last million years, whereas the branch leading to M. ciliata and its sister species has sustained rearrangement at a rate that is at least ten times lower. Furthermore, our analysis of published data shows that rates of mitochondrial genome rearrangement in seed plants vary by at least 600-fold. We find that sites of rearrangement are highly preferentially located in very close proximity to repeated sequences in Monsonia. This provides strong support for the hypothesis that rearrangement in angiosperm mitochondrial genomes occurs largely through repeat-mediated recombination. Because there is little variation in the amount of repeat sequence among Monsonia genomes, the variable rates of rearrangement in Monsonia probably reflect variable rates of mitochondrial recombination itself. Finally, we show that mitochondrial synonymous substitutions occur in a clock-like manner in Monsonia; rates of mitochondrial substitutions and rearrangements are therefore highly uncoupled in this group.

Tension and Resolution: Dynamic, Evolving Populations of Organelle Genomes within Plant Cells

Mitochondria and plastids form dynamic, evolving populations physically embedded in the fluctuating environment of the plant cell. Their evolutionary heritage has shaped how the cell controls the genetic structure and the physical behavior of its organelle populations. While the specific genes involved in these processes are gradually being revealed, the governing principles underlying this controlled behavior remain poorly understood. As the genetic and physical dynamics of these organelles are central to bioenergetic performance and plant physiology, this challenges both fundamental biology and strategies to engineer better-performing plants. This article reviews current knowledge of the physical and genetic behavior of mitochondria and chloroplasts in plant cells. An overarching hypothesis is proposed whereby organelles face a tension between genetic robustness and individual control and responsiveness, and different species resolve this tension in different ways. As plants are immobile and thus subject to fluctuating environments, their organelles are proposed to favor individual responsiveness, sacrificing genetic robustness. Several notable features of plant organelles, including large genomes, mtDNA recombination, fragmented organelles, and plastid/mitochondrial differences may potentially be explained by this hypothesis. Finally, the ways that quantitative and systems biology can help shed light on the plethora of open questions in this field are highlighted.

Repeats of Unusual Size in Plant Mitochondrial Genomes: Identification, Incidence and Evolution

Plant mitochondrial genomes have excessive size relative to coding capacity, a low mutation rate in genes and a high rearrangement rate. They also have abundant non-tandem repeats often including pairs of large repeats which cause isomerization of the genome by recombination, and numerous repeats of up to several hundred base pairs that recombine only when the genome is stressed by DNA damaging agents or mutations in DNA repair pathway genes. Early work on mitochondrial genomes led to the suggestion that repeats in the size range from several hundred to a few thousand base pair are underrepresented. The repeats themselves are not well-conserved between species, and are not always annotated in mitochondrial sequence assemblies. We systematically identified and compared these repeats, which are important clues to mechanisms of DNA maintenance in mitochondria. We developed a tool to find and curate non-tandem repeats larger than 50bp and analyzed the complete mitochondrial sequences from 157 plant species. We observed an interesting difference between taxa: the repeats are larger and more frequent in the vascular plants. Analysis of closely related species also shows that plant mitochondrial genomes evolve in dramatic bursts of breakage and rejoining, complete with DNA sequence gain and loss. We suggest an adaptive explanation for the existence of the repeats and their evolution.

Homologous recombination changes the context of Cytochrome b transcription in the mitochondrial genome of Silene vulgaris KRA

Background

Silene vulgaris (bladder campion) is a gynodioecious species existing as two genders – male-sterile females and hermaphrodites. Cytoplasmic male sterility (CMS) is generally encoded by mitochondrial genes, which interact with nuclear fertility restorer genes. Mitochondrial genomes of this species vary in DNA sequence, gene order and gene content. Multiple CMS genes are expected to exist in S. vulgaris, but little is known about their molecular identity.

Results

We assembled the complete mitochondrial genome from the haplotype KRA of S. vulgaris. It consists of five chromosomes, two of which recombine with each other. Two small non-recombining chromosomes exist in linear, supercoiled and relaxed circle forms. We compared the mitochondrial transcriptomes from females and hermaphrodites and confirmed the differentially expressed chimeric gene bobt as the strongest CMS candidate gene in S. vulgaris KRA. The chimeric gene bobt is co-transcribed with the Cytochrome b (cob) gene in some genomic configurations. The co-transcription of a CMS factor with an essential gene may constrain transcription inhibition as a mechanism for fertility restoration because of the need to maintain appropriate production of the necessary protein. Homologous recombination places the gene cob outside the control of bobt, which allows for the suppression of the CMS gene by the fertility restorer genes. We found the loss of three editing sites in the KRA mitochondrial genome and identified four sites with highly distinct editing rates between KRA and another S. vulgaris haplotypes (KOV). Three of these highly differentially edited sites were located in the transport membrane protein B (mttB) gene. They resulted in differences in MttB protein sequences between haplotypes.

Conclusions

Frequent homologous recombination events that are widespread in plant mitochondrial genomes may change chromosomal configurations and also the control of gene transcription including CMS gene expression. Posttranscriptional processes, e.g. RNA editing shall be evaluated in evolutionary and co-evolutionary studies of mitochondrial genes, because they may change protein composition despite the sequence identity of the respective genes. The investigation of natural populations of wild species such as S. vulgaris are necessary to reveal important aspects of CMS missed in domesticated crops, the traditional focus of the CMS studies.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5254-0) contains supplementary material, which is available to authorized users.

Plant organellar DNA polymerases paralogs exhibit dissimilar nucleotide incorporation fidelity

The coding sequences of plant mitochondrial and chloroplast genomes present a lower mutation rate than the coding sequences of animal mitochondria. However, plant mitochondrial genomes frequently rearrange and present high mutation rates in their noncoding sequences. DNA replication in plant organelles is carried out by two DNA polymerases (DNAP) paralogs. In Arabidopsis thaliana at least one DNAP paralog (AtPolIA or AtPolIB) is necessary for plant viability, suggesting that both genes are partially redundant. To understand how AtPolIs replicate genomes that present low and high mutation rates, we measured their nucleotide incorporation for all 16-base pair combinations in vitro. AtPolIA presents an error rate of 7.26 × 10^-5 , whereas AtPolIB has an error rate of 5.45 × 10^-4 . Thus, AtPolIA and AtPolIB are 3.5 and 26-times less accurate than human mitochondrial DNAP γ. The 8-fold difference in fidelity between both AtPolIs results from a higher catalytic efficiency in AtPolIA. Both AtPolIs extend from mismatches and the fidelity of AtPolIs ranks between high fidelity and lesion bypass DNAPs. The different nucleotide incorporation fidelity between AtPolIs predicts a prevalent role of AtPolIA in DNA replication and AtPolIB in DNA repair. We hypothesize that in plant organelles, DNA mismatches generated during DNA replication are repaired via recombination-mediated or DNA mismatch repair mechanisms that selectively target the coding region and that the mismatches generated by AtPolIs may result in the frequent expansion and rearrangements present in plant mitochondrial genomes.

Causes and Consequences of Rapidly Evolving mtDNA in a Plant Lineage

Understanding mechanisms of coevolution between nuclear and mitochondrial (mt) genomes is a defining challenge in eukaryotic genetics. The angiosperm genus Silene is a natural system to investigate the causes and consequences of mt mutation rate variation because closely related species have highly divergent rates. In Silene species with fast-evolving mtDNA, nuclear genes that encode mitochondrially targeted proteins (N-mt genes) are also fast-evolving. This correlation could indicate positive selection to compensate for mt mutations, but might also result from a recent relaxation of selection. To differentiate between these interpretations, we used phylogenetic and population-genetic methods to test for positive and relaxed selection in three classes of N-mt genes (oxidative phosphorylation genes, ribosomal genes, and “RRR” genes involved in mtDNA recombination, replication, and repair). In all three classes, we found that species with fast-evolving mtDNA had: 1) elevated d_N/d_S, 2) an excess of nonsynonymous divergence relative to levels of intraspecific polymorphism, which is a signature of positive selection, and 3) no clear signals of relaxed selection. “Control” genes exhibited comparatively few signs of positive selection. These results suggest that high mt mutation rates can create selection on N-mt genes and that relaxed selection is an unlikely cause of recent accelerations in the evolution of N-mt genes. Because mt-RRR genes were found to be under positive selection, it is unlikely that elevated mt mutation rates in Silene were caused by inactivation of these mt-RRR genes. Therefore, the causes of extreme increases in angiosperm mt mutation rates remain uncertain.

Positive Selection Driving Cytoplasmic Genome Evolution of the Medicinally Important Ginseng Plant Genus Panax

Panax L. (the ginseng genus) is a shade-demanding group within the family Araliaceae and all of its species are of crucial significance in traditional Chinese medicine. Phylogenetic and biogeographic analyses demonstrated that two rounds of whole genome duplications accompanying with geographic and ecological isolations promoted the diversification of Panax species. However, contributions of the cytoplasmic genomes to the adaptive evolution of Panax species remained largely uninvestigated. In this study, we sequenced the chloroplast and mitochondrial genomes of 11 accessions belonging to seven Panax species. Our results show that heterogeneity in nucleotide substitution rate is abundant in both of the two cytoplasmic genomes, with the mitochondrial genome possessing more variants at the total level but the chloroplast showing higher sequence polymorphisms at the genic regions. Genome-wide scanning of positive selection identified five and 12 genes from the chloroplast and mitochondrial genomes, respectively. Functional analyses further revealed that these selected genes play important roles in plant development, cellular metabolism and adaptation. We therefore conclude that positive selection might be one of the potential evolutionary forces that shaped nucleotide variation pattern of these Panax species. In particular, the mitochondrial genes evolved under stronger selective pressure compared to the chloroplast genes.

Rapid evolutionary divergence of diploid and allotetraploid Gossypium mitochondrial genomes

Background

Cotton (Gossypium spp.) is commonly grouped into eight diploid genomic groups and an allotetraploid genomic group, AD. The mitochondrial genomes supply new information to understand both the evolution process and the mechanism of cytoplasmic male sterility. Based on previously released mitochondrial genomes of G. hirsutum (AD₁), G. barbadense (AD₂), G. raimondii (D₅) and G. arboreum (A₂), together with data of six other mitochondrial genomes, to elucidate the evolution and diversity of mitochondrial genomes within Gossypium.

Results

Six Gossypium mitochondrial genomes, including three diploid species from D and three allotetraploid species from AD genome groups (G. thurberi D₁, G. davidsonii D_3-d and G. trilobum D₈; G. tomentosum AD₃, G. mustelinum AD₄ and G. darwinii AD₅), were assembled as the single circular molecules of lengths about 644 kb in diploid species and 677 kb in allotetraploid species, respectively. The genomic structures of mitochondrial in D group species were identical but differed from the mitogenome of G. arboreum (A₂), as well as from the mitogenomes of five species of the AD group. There mainly existed four or six large repeats in the mitogenomes of the A + AD or D group species, respectively. These variations in repeat sequences caused the major inversions and translocations within the mitochondrial genome. The mitochondrial genome complexity in Gossypium presented eight unique segments in D group species, three specific fragments in A + AD group species and a large segment (more than 11 kb) in diploid species. These insertions or deletions were most probably generated from crossovers between repetitive or homologous regions. Unlike the highly variable genome structure, evolutionary distance of mitochondrial genes was 1/6th the frequency of that in chloroplast genes of Gossypium. RNA editing events were conserved in cotton mitochondrial genes. We confirmed two near full length of the integration of the mitochondrial genome into chromosome 1 of G. raimondii and chromosome A03 of G. hirsutum, respectively, with insertion time less than 1.03 MYA.

Conclusion

Ten Gossypium mitochondrial sequences highlight the insights to the evolution of cotton mitogenomes.

Electronic supplementary material

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

Mitochondrial genome evolution in Alismatales: Size reduction and extensive loss of ribosomal protein genes

The order Alismatales is a hotspot for evolution of plant mitochondrial genomes characterized by remarkable differences in genome size, substitution rates, RNA editing, retrotranscription, gene loss and intron loss. Here we have sequenced the complete mitogenomes of Zostera marina and Stratiotes aloides, which together with previously sequenced mitogenomes from Butomus and Spirodela, provide new evolutionary evidence of genome size reduction, gene loss and transfer to the nucleus. The Zostera mitogenome includes a large portion of DNA transferred from the plastome, yet it is the smallest known mitogenome from a non-parasitic plant. Using a broad sample of the Alismatales, the evolutionary history of ribosomal protein gene loss is analyzed. In Zostera almost all ribosomal protein genes are lost from the mitogenome, but only some can be found in the nucleus.