Similar Relative Mutation Rates in the Three Genetic Compartments of Mesostigma and Chlamydomonas

A complex interplay of evolutionary forces continues to shape ancient co-occurring symbiont genomes

Many insects depend on ancient associations with intracellular bacteria for essential nutrition. The genomes of these bacteria are often highly reduced. Although drift is a major driver of symbiont evolution, other evolutionary forces continue to influence them. To understand how ongoing molecular evolution and gene loss shape symbiont genomes, we sequenced two of the most ancient symbionts known, Sulcia and Nasuia, from 20 Hawaiian Nesophrosyne leafhoppers. We leveraged the parallel divergence of Nesophrosyne lineages throughout Hawaii as a natural experimental framework. Sulcia and Nasuia experience ongoing-but divergent-gene loss, often in a convergent fashion. Although some genes are under relaxed selection, purifying and positive selection are also important drivers of genome evolution, particularly in maintaining certain nutritional and cellular functions. Our results further demonstrate that symbionts experience dramatically different evolutionary environments, as evidenced by the finding that Sulcia and Nasuia have one of the slowest and fastest rates of molecular evolution known.

Plastid genomes and phylogenomics of liverworts (Marchantiophyta): Conserved genome structure but highest relative plastid substitution rate in land plants

With some 7300 species of small nonvascular spore-producing plants, liverworts represent one of the major lineages of land plants. Although multi-locus molecular phylogenetic studies have elucidated relationships of liverworts at different taxonomic categories, the backbone phylogeny of liverworts is still to be fully resolved, especially for the placement of Ptilidiales and the relationships within Jungermanniales and Marchantiales. Here, we provided phylogenomic inferences of liverworts based on 42 newly sequenced and 24 published liverwort plastid genomes representing all but two orders of liverworts, and characterized the evolution of the plastome in liverworts. The structure of the plastid genome is overall conserved across the phylogeny of liverworts, with only two structural variants detected from simple thalloids, besides 18 out of 43 liverwort genera showing intron variations in their plastomes. Complex thalloid liverworts maintain the most plastid genes, and seem to undergo fewer gene deletions and pseudogenization events than other liverworts. Plastid phylogenetic inferences yielded mostly robustly supported relationships, and consistently resolved Ptilidiales as the sister to Porellales. The relative ratio of silent substitutions across the three genetic compartments (i.e., 1:15:10, for mitochondrial:plastid:nuclear) suggests that liverwort plastid genes have the potential to evolve faster than their nuclear counterparts, unlike in any other major land plant lineages where the mutation rate of nuclear genes overwhelm those of their plastid and mitochondrial counterparts.

Evolutionary Rates are Correlated Between Buchnera Endosymbionts and the Mitochondrial Genomes of Their Aphid Hosts

The evolution of bacterial endosymbiont genomes is strongly influenced by host-driven selection. Factors affecting host genome evolution will potentially affect endosymbiont genomes in similar ways. One potential outcome is correlations in molecular rates between the genomes of the symbiotic partners. Recently, we presented the first evidence of such correlations between the mitochondrial genomes of cockroaches and the genomes of their endosymbiont (Blattabacterium cuenoti). Here we investigate whether similar patterns are found in additional host-symbiont partners. We use partial genome data from multiple strains of the bacterial endosymbionts Buchnera aphidicola and Sulcia muelleri, and the mitochondrial genomes of their sap-feeding insect hosts. Both endosymbionts show phylogenetic congruence with the mitochondria of their hosts, a result that is expected due to their identical mode of inheritance. We compared root-to-tip distances and branch lengths of phylogenetically independent species pairs. Both analyses showed a highly significant correlation of molecular rates between the genomes of Buchnera and the mitochondrial genomes of their hosts. A similar correlation was detected between Sulcia and their hosts, but was not statistically significant. Our results indicate that evolutionary rate correlations between hosts and long-term symbionts may be a widespread phenomenon.

Distinctive evolutionary pattern of organelle genomes linked to the nuclear genome in Selaginellaceae

Plastids and mitochondria are endosymbiotic organelles that store genetic information. The genomes of these organelles generally exhibit contrasting patterns regarding genome architecture and genetic content. However, they have similar genetic features in Selaginellaceae, and little is known about what causes parallel evolution. Here, we document the multipartite plastid genomes (plastomes) and the highly divergent mitochondrial genomes (mitogenomes) from spikemoss obtained by combining short- and long-reads. The 188-kb multipartite plastome has three ribosomal operon copies in the master genomic conformation, creating the alternative subgenomic conformation composed of 110- and 78-kb subgenomes. The long-read data indicated that the two different genomic conformations were present in almost equal proportions in the plastomes of Selaginella nipponica. The mitogenome of S. nipponica was assembled into 27 contigs with a total size of 110 kb. All contigs contained directly arranged repeats at both ends, which introduced multiple conformations. Our results showed that plastomes and mitogenomes share high tRNA losses, GC-biased nucleotides, elevated substitution rates and complicated organization. The exploration of nuclear-encoded organelle DNA replication, recombination and repair proteins indicated that, several single-targeted proteins, particularly plastid-targeted recombinase A1, have been lost in Selaginellaceae; conversely, the dual-targeted proteins remain intact. According to the reported function of recombinase A1, we propose that the plastomes of spikemoss often fail to pair homologous sequences during recombination, and the dual-targeted proteins play a key role in the convergent genetic features of plastomes and mitogenomes. Our results provide a distinctive evolutionary pattern of the organelle genomes in Selaginellaceae and evidence of their convergent evolution.

First Estimation of the Spontaneous Mutation Rate in Diatoms

Mutations are the origin of genetic diversity, and the mutation rate is a fundamental parameter to understand all aspects of molecular evolution. The combination of mutation-accumulation experiments and high-throughput sequencing enabled the estimation of mutation rates in most model organisms, but several major eukaryotic lineages remain unexplored. Here, we report the first estimation of the spontaneous mutation rate in a model unicellular eukaryote from the Stramenopile kingdom, the diatom Phaeodactylum tricornutum (strain RCC2967). We sequenced 36 mutation accumulation lines for an average of 181 generations per line and identified 156 de novo mutations. The base substitution mutation rate per site per generation is μbs = 4.77 × 10-10 and the insertion-deletion mutation rate is μid = 1.58 × 10-11. The mutation rate varies as a function of the nucleotide context and is biased toward an excess of mutations from GC to AT, consistent with previous observations in other species. Interestingly, the mutation rates between the genomes of organelles and the nucleus differ, with a significantly higher mutation rate in the mitochondria. This confirms previous claims based on indirect estimations of the mutation rate in mitochondria of photosynthetic eukaryotes that acquired their plastid through a secondary endosymbiosis. This novel estimate enables us to infer the effective population size of P. tricornutum to be Ne∼8.72 × 106.

Evolutionary rates are correlated between cockroach symbionts and mitochondrial genomes

Bacterial endosymbionts evolve under strong host-driven selection. Factors influencing host evolution might affect symbionts in similar ways, potentially leading to correlations between the molecular evolutionary rates of hosts and symbionts. Although there is evidence of rate correlations between mitochondrial and nuclear genes, similar investigations of hosts and symbionts are lacking. Here, we demonstrate a correlation in molecular rates between the genomes of an endosymbiont (Blattabacterium cuenoti) and the mitochondrial genomes of their hosts (cockroaches). We used partial genome data for multiple strains of B. cuenoti to compare phylogenetic relationships and evolutionary rates for 55 cockroach/symbiont pairs. The phylogenies inferred for B. cuenoti and the mitochondrial genomes of their hosts were largely congruent, as expected from their identical maternal and cytoplasmic mode of inheritance. We found a correlation between evolutionary rates of the two genomes, based on comparisons of root-to-tip distances and on comparisons of the branch lengths of phylogenetically independent species pairs. Our results underscore the profound effects that long-term symbiosis can have on the biology of each symbiotic partner.

Relative Mutation Rates in Nucleomorph-Bearing Algae

Chlorarachniophyte and cryptophyte algae are unique among plastid-containing species in that they have a nucleomorph genome: a compact, highly reduced nuclear genome from a photosynthetic eukaryotic endosymbiont. Despite their independent origins, the nucleomorph genomes of these two lineages have similar genomic architectures, but little is known about the evolutionary pressures impacting nucleomorph DNA, particularly how their rates of evolution compare to those of the neighboring genetic compartments (the mitochondrion, plastid, and nucleus). Here, we use synonymous substitution rates to estimate relative mutation rates in the four genomes of nucleomorph-bearing algae. We show that the relative mutation rates of the host versus endosymbiont nuclear genomes are similar in both chlorarachniophytes and cryptophytes, despite the fact that nucleomorph gene sequences are notoriously highly divergent. There is some evidence, however, for slightly elevated mutation rates in the nucleomorph DNA of chlorarachniophytes-a feature not observed in that of cryptophytes. For both lineages, relative mutation rates in the plastid appear to be lower than those in the nucleus and nucleomorph (and, in one case, the mitochondrion), which is consistent with studies of other plastid-bearing protists. Given the divergent nature of nucleomorph genes, our finding of relatively low evolutionary rates in these genomes suggests that for both lineages a burst of evolutionary change and/or decreased selection pressures likely occurred early in the integration of the secondary endosymbiont.

The mutational hazard hypothesis of organelle genome evolution: 10 years on

Why is there such a large variation in size and noncoding DNA content among organelle genomes? One explanation is that this genomic variation results from differences in the rates of organelle mutation and random genetic drift, as opposed to being the direct product of natural selection. Along these lines, the mutational hazard hypothesis (MHH) holds that 'excess' DNA is a mutational liability (because it increases the potential for harmful mutations) and, thus, has a greater tendency to accumulate in an organelle system with a low mutation rate as opposed to one with a high rate of mutation. Various studies have explored this hypothesis and, more generally, the relationship between organelle genome architecture and the mode and efficiency of organelle DNA repair. Although some of these investigations are in agreement with the MHH, others have contradicted it; nevertheless, they support a central role of mutation, DNA maintenance pathways and random genetic drift in fashioning organelle chromosomes. Arguably, one of the most important contributions of the MHH is that it has sparked crucial, widespread discussions about the importance of nonadaptive processes in genome evolution.

Rickettsial endosymbiont in the "early-diverging" streptophyte green alga Mesostigma viride

A bacterial endosymbiont was unexpectedly found in the "axenic" culture strain of the streptophyte green alga Mesostigma viride (NIES-995). Phylogenetic analyses based on 16S rRNA gene sequences showed that the symbiont belongs to the order Rickettsiales, specifically to the recently designated clade "Candidatus Megaira," which is closely related to the well-known Rickettsia clade. Rickettsiales bacteria of the "Ca. Megaira" clade are found in a taxonomically diverse array of eukaryotic hosts, including chlorophycean green algae, several ciliate species, and invertebrates such as Hydra. Transmission electron microscopy, fluorescence in situ hybridi-zation, and SYBR Green I staining experiments revealed that the endosymbiont of M. viride NIES-995 is rod shaped, typically occurs in clusters, and is surrounded by a halo-like structure, presumably formed by secretory substances from the bacterium. Two additional M. viride strains (NIES-296 and NIES-475), but not SAG50-1, were found to house the rickettsial endosymbiont. Analyses of strain NIES-995 transcriptome data indicated the presence of at least 91 transcriptionally active genes of symbiont origins. These include genes for surface proteins (e.g., rOmpB) that are known to play key roles in bacterial attachment onto host eukaryotes in related Rickettsia species. The assembled M. viride transcriptome includes transcripts that code for a suite of predicted algal-derived proteins, such as Ku70, WASH, SCAR, and CDC42, which may be important in the formation of the algal-rickettsial association.

Auxenochlorella protothecoides and Prototheca wickerhamii plastid genome sequences give insight into the origins of non-photosynthetic algae

The forfeiting of photosynthetic capabilities has occurred independently many times throughout eukaryotic evolution. But almost all non-photosynthetic plants and algae still retain a colorless plastid and an associated genome, which performs fundamental processes apart from photosynthesis. Unfortunately, little is known about the forces leading to photosynthetic loss; this is largely because there is a lack of data from transitional species. Here, we compare the plastid genomes of two “transitional” green algae: the photosynthetic, mixotrophic Auxenochlorella protothecoides and the non-photosynthetic, obligate heterotroph Prototheca wickerhamii. Remarkably, the plastid genome of A. protothecoides is only slightly larger than that of P. wickerhamii, making it among the smallest plastid genomes yet observed from photosynthetic green algae. Even more surprising, both algae have almost identical plastid genomic architectures and gene compositions (with the exception of genes involved in photosynthesis), implying that they are closely related. This close relationship was further supported by phylogenetic and substitution rate analyses, which suggest that the lineages giving rise to A. protothecoides and P. wickerhamii diverged from one another around six million years ago.

Guanchochroma wildpretii gen. et spec. nov. (Ochrophyta) Provides New Insights into the Diversification and Evolution of the Algal Class Synchromophyceae

A new relative of the chrysophyte genus Chrysopodocystis was found in Tenerife and termed Guanchochroma wildpretii. This unicellular alga was most noticeably discernible from Chrysopodocystis socialis (the only species of this genus) by the presence of a cyst-like stage with a multilayered lorica, which also functions as a dispersal unit and shows secondary wall growth. Secondary expansion of loricae (cell casings not involved in cell division, usually with a more or less pronounced opening) has never been observed previously and marks a unique feature of the new taxon. Plastids are non-randomly distributed within cells of G. wildpretii. 18S rRNA gene analyses identified the two species as sister lineages and placed them in a monophyletic group with the Synchromophyceae, a heterokont algal (Ochrophyta) class characterized by the presence of chloroplast complexes. Yet, neither Chrysopodocystis nor Guanchochroma showed this feature in ultrastructure analyses. Additionally, their 18S rRNA genes possessed distinct inserts, the highest GC-content known for Ochrophyta and exceptionally long branches on the Ochrophyta 18S rDNA phylogenetic tree, suggesting substantially increased substitution rates along their branch compared to Synchromophyceae. Plastid marker data (rbcL) recovered a monophyletic clade of Chrysopodocystis, Guanchochroma and Synchromophyceae as well, yet with lower supports for internal split order due to limited resolution of the marker. Evidence for the sequence of events leading to the formation of the plastid complex of Synchromophyceae still remains ambiguous because of the apparently short timeframe in which they occurred.

Mutation rates in plastid genomes: they are lower than you might think

Within plastid-bearing species, the mutation rate of the plastid genome is often assumed to be greater than that of the mitochondrial genome. This assumption is based on early, pioneering studies of land plant molecular evolution, which uncovered higher rates of synonymous substitution in plastid versus mitochondrial DNAs. However, much of the plastid-containing eukaryotic diversity falls outside of land plants, and the patterns of plastid DNA evolution for embryophytes do not necessarily reflect those of other groups. Recent analyses of plastid and mitochondrial substitution rates in diverse lineages have uncovered very different trends than those recorded for land plants. Here, I explore these new data and argue that for many protists the plastid mutation rate is lower than that of the mitochondrion, including groups with primary or secondary plastids as well as nonphotosynthetic algae. These findings have far-reaching implications for how we view plastid genomes and how their sequences are used for evolutionary analyses, and might ultimately reflect a general tendency toward more efficient DNA repair mechanisms in plastids than in mitochondria.

Mitochondrial and plastid genome architecture: Reoccurring themes, but significant differences at the extremes

Mitochondrial and plastid genomes show a wide array of architectures, varying immensely in size, structure, and content. Some organelle DNAs have even developed elaborate eccentricities, such as scrambled coding regions, nonstandard genetic codes, and convoluted modes of posttranscriptional modification and editing. Here, we compare and contrast the breadth of genomic complexity between mitochondrial and plastid chromosomes. Both organelle genomes have independently evolved many of the same features and taken on similar genomic embellishments, often within the same species or lineage. This trend is most likely because the nuclear-encoded proteins mediating these processes eventually leak from one organelle into the other, leading to a high likelihood of processes appearing in both compartments in parallel. However, the complexity and intensity of genomic embellishments are consistently more pronounced for mitochondria than for plastids, even when they are found in both compartments. We explore the evolutionary forces responsible for these patterns and argue that organelle DNA repair processes, mutation rates, and population genetic landscapes are all important factors leading to the observed convergence and divergence in organelle genome architecture.

Massive and widespread organelle genomic expansion in the green algal genus Dunaliella

The mitochondrial genomes of chlamydomonadalean green algae are renowned for their highly reduced and conserved gene repertoires, which are almost fixed at 12 genes across the entire lineage. The sizes of these genomes, however, are much more variable, with some species having small, compact mitochondrial DNAs (mtDNAs) and others having expanded ones. Earlier work demonstrated that the halophilic genus Dunaliella contains extremely inflated organelle genomes, but to date the mtDNA of only one isolate has been explored. Here, by surveying mtDNA architecture across the Chlamydomonadales, we show that various Dunaliella species have undergone massive levels of mitochondrial genomic expansion, harboring the most inflated, intron-dense mtDNAs available from chlorophyte green algae. The same also appears to be true for their plastid genomes, which are potentially among the largest of all plastid-containing eukaryotes. Genetic divergence data are used to investigate the underlying causes of such extreme organelle genomic architectures, and ultimately reveal order-of-magnitude differences in mitochondrial versus plastid mutation rates within Dunaliella.

Nucleotide substitution analyses of the glaucophyte Cyanophora suggest an ancestrally lower mutation rate in plastid vs mitochondrial DNA for the Archaeplastida

A lot is known about the evolution and architecture of plastid, mitochondrial, and nuclear genomes, but surprisingly little is known about their relative rates of mutation. Most available relative-rate data come from seed plants, which, with few exceptions, have a mitochondrial mutation rate that is lower than those of the plastid and nucleus. But new findings from diverse plastid-bearing lineages have shown that for some eukaryotes the mitochondrial mutation rate is an order of magnitude greater than those of the plastid and nucleus. Here, we explore for the first time relative rates of mutation within the Glaucophyta-one of three main lineages that make up the Archaeplastida (or Plantae sensu lato). Nucleotide substitution analyses from distinct isolates of the unicellular glaucophyte Cyanophora paradoxa reveal 4-5-fold lower rates of mutation in the plastid and nucleus than the mitochondrion, which is similar to the mutational pattern observed in red algae and haptophytes, but opposite to that of seed plants. These data, together with data from previous reports, suggest that for much of the known photosynthetic eukaryotic diversity, plastid DNA mutations occur less frequently than those in mitochondrial DNA.

Organellar inheritance in the green lineage: insights from Ostreococcus tauri

Along the green lineage (Chlorophyta and Streptophyta), mitochondria and chloroplast are mainly uniparentally transmitted and their evolution is thus clonal. The mode of organellar inheritance in their ancestor is less certain. The inability to make clear phylogenetic inference is partly due to a lack of information for deep branching organisms in this lineage. Here, we investigate organellar evolution in the early branching green alga Ostreococcus tauri using population genomics data from the complete mitochondrial and chloroplast genomes. The haplotype structure is consistent with clonal evolution in mitochondria, while we find evidence for recombination in the chloroplast genome. The number of recombination events in the genealogy of the chloroplast suggests that recombination, and thus biparental inheritance, is not rare. Consistent with the evidence of recombination, we find that the ratio of the number of nonsynonymous to the synonymous polymorphisms per site is lower in chloroplast than in the mitochondria genome. We also find evidence for the segregation of two selfish genetic elements in the chloroplast. These results shed light on the role of recombination and the evolutionary history of organellar inheritance in the green lineage.

Relative rates of evolution among the three genetic compartments of the red alga Porphyra differ from those of green plants and do not correlate with genome architecture

In photosynthetic eukaryotes, relative silent-site nucleotide substitution rates (which can be used to approximate relative mutation rates) among mitochondrial, plastid, and nuclear genomes (mtDNAs, ptDNAs, and nucDNAs) are estimated to be 1:3:10 respectively for seed plants and roughly equal for green algae. These estimates correlate with certain genome characteristics, such as size and coding density, and have therefore been taken to support a relationship between mutation rate and genome architecture. Plants and green algae, however, represent a small fraction of the major eukaryotic plastid-bearing lineages. Here, we investigate relative rates of mutation within the model red algal genus Porphyra. In contrast to plants, we find that the levels of silent-site divergence between the Porphyra purpurea and Porphyra umbilicalis mtDNAs are three times that of their ptDNAs and five times that of their nucDNAs. Moreover, relative mutation rates do not correlate with genome architecture: despite an estimated three-fold difference in their mutation rate, the mitochondrial and plastid genome coding densities are equivalent - an observation that extends to organisms with secondary red algal plastids. These findings are supported by within-species silent-site polymorphism data from P. purpurea.