Spontaneous Mutation Rate in the Smallest Photosynthetic Eukaryotes

Low Mutation Rate and Atypical Mutation Spectrum in Prasinoderma coloniale: Insights From an Early Diverging Green Lineage

Mutations are the ultimate source of genetic diversity on which natural selection and genetic drift act, playing a crucial role in evolution and long-term adaptation. At the molecular level, the spontaneous mutation rate (µ), defined as the number of mutations per base per generation, thus determines the adaptive potential of a species. Through a mutation accumulation experiment, we estimate the mutation rate and spectrum in Prasinoderma coloniale, a phytoplankton species from an early-branching lineage within the Archaeplastida, characterized by an unusually high genomic guanine-cytosine (GC) content (69.8%). We find that P. coloniale has a very low total mutation rate of µ = 2.00 × 10-10. The insertion-deletion mutation rate is almost 5 times lesser than the single nucleotide mutation rate with µID = 3.40 × 10-11 and µSNM = 1.62 × 10-10. Prasinoderma coloniale also exhibits an atypical mutational spectrum: While essentially all other eukaryotes show a bias toward GC to AT mutations, no evidence of this AT-bias is observed in P. coloniale. Since cytosine methylation is known to be mutagenic, we hypothesized that this may result from an absence of C-methylation. Surprisingly, we found high levels of C-methylation (14% in 5mC, 25% in 5mCG contexts). Methylated cytosines did not show increased mutation rates compared with unmethylated ones, not supporting the prevailing notion that C-methylation universally leads to higher mutation rates. Overall, P. coloniale combines a GC-rich genome with a low mutation rate and original mutation spectrum, suggesting the almost universal AT-bias may not have been present in the ancestor of the green lineage.

Flipping the switch on some of the slowest mutating genomes: Direct measurements of plant mitochondrial and plastid mutation rates in msh1 mutants

Plant mitochondrial and plastid genomes have exceptionally slow rates of sequence evolution, and recent work has identified an unusual member of the MutS gene family ("plant MSH1") as being instrumental in preventing point mutations in these genomes. However, the effects of disrupting MSH1-mediated DNA repair on "germline" mutation rates have not been quantified. Here, we used Arabidopsis thaliana mutation accumulation (MA) lines to measure mutation rates in msh1 mutants and matched wild type (WT) controls. We detected 124 single nucleotide variants (SNVs: 49 mitochondrial and 75 plastid) and 668 small insertions and deletions (indels: 258 mitochondrial and 410 plastid) in msh1 MA lines. In striking contrast, we did not find any organelle mutations in the WT MA lines, and reanalysis of data from a much larger WT MA experiment also failed to detect any variants. The observed number of SNVs in the msh1 MA lines corresponds to estimated mutation rates of 6.1×10-7 and 3.2 ×10-6 per bp per generation in mitochondrial and plastid genomes, respectively. These rates exceed those of species known to have very high mitochondrial mutation rates (e.g., nematodes and fruit flies) by an order of magnitude or more and are on par with estimated rates in humans despite the generation times of A. thaliana being nearly 100-fold shorter. Therefore, disruption of a single plant-specific genetic factor in A. thaliana is sufficient to erase or even reverse the enormous difference in organelle mutation rates between plants and animals.

Protists and protistology in the Anthropocene: challenges for a climate and ecological crisis

Eukaryotic microorganisms, or "protists," while often inconspicuous, play fundamental roles in the Earth ecosystem, ranging from primary production and nutrient cycling to interactions with human health and society. In the backdrop of accelerating climate dysregulation, alongside anthropogenic disruption of natural ecosystems, understanding changes to protist functional and ecological diversity is of critical importance. In this review, we outline why protists matter to our understanding of the global ecosystem and challenges of predicting protist species resilience and fragility to climate change. Finally, we reflect on how protistology may adapt and evolve in a present and future characterized by rapid ecological change.

Genome-wide patterns of selection-drift variation strongly associate with organismal traits across the green plant lineage

There are many gaps in our knowledge of how life cycle variation and organismal body architecture associate with molecular evolution. Using the diverse range of green algal body architectures and life cycle types as a test case, we hypothesize that increases in cytomorphological complexity are likely to be associated with a decrease in the effective population size, because larger-bodied organisms typically have smaller populations, resulting in increased drift. For life cycles, we expect haploid-dominant lineages to evolve under stronger selection intensity relative to diploid-dominant life cycles owing to masking of deleterious alleles in heterozygotes. We use a genome-scale data set spanning the phylogenetic diversity of green algae and phylogenetic comparative approaches to measure the relative selection intensity across different trait categories. We show stronger signatures of drift in lineages with more complex body architectures compared with unicellular lineages, which we consider to be a consequence of smaller effective population sizes of the more complex algae. Significantly higher rates of synonymous as well as nonsynonymous substitutions relative to other algal body architectures highlight that siphonous and siphonocladous body architectures, characteristic of many green seaweeds, form an interesting test case to study the potential impacts of genome redundancy on molecular evolution. Contrary to expectations, we show that levels of selection efficacy do not show a strong association with life cycle types in green algae. Taken together, our results underline the prominent impact of body architecture on the molecular evolution of green algal genomes.

H3K4me1 recruits DNA repair proteins in plants

DNA repair proteins can be recruited by their histone reader domains to specific epigenomic features, with consequences on intragenomic mutation rate variation. Here, we investigated H3K4me1-associated hypomutation in plants. We first examined 2 proteins which, in plants, contain Tudor histone reader domains: PRECOCIOUS DISSOCIATION OF SISTERS 5 (PDS5C), involved in homology-directed repair, and MUTS HOMOLOG 6 (MSH6), a mismatch repair protein. The MSH6 Tudor domain of Arabidopsis (Arabidopsis thaliana) binds to H3K4me1 as previously demonstrated for PDS5C, which localizes to H3K4me1-rich gene bodies and essential genes. Mutations revealed by ultradeep sequencing of wild-type and msh6 knockout lines in Arabidopsis show that functional MSH6 is critical for the reduced rate of single-base substitution (SBS) mutations in gene bodies and H3K4me1-rich regions. We explored the breadth of these mechanisms among plants by examining a large rice (Oryza sativa) mutation data set. H3K4me1-associated hypomutation is conserved in rice as are the H3K4me1-binding residues of MSH6 and PDS5C Tudor domains. Recruitment of DNA repair proteins by H3K4me1 in plants reveals convergent, but distinct, epigenome-recruited DNA repair mechanisms from those well described in humans. The emergent model of H3K4me1-recruited repair in plants is consistent with evolutionary theory regarding mutation modifier systems and offers mechanistic insight into intragenomic mutation rate variation in plants.

Prognostic analysis of mutated genes and insight into effects of DNA damage and repair on mutational strand asymmetries in gastric cancer

Gastric cancer (GACA) is a complex and multifaceted disease influenced by a variety of environmental and genetic factors. Somatic mutations play a major role in its development, and their characteristics, including the asymmetry between two DNA strands, are of great interest and appear as a signal of information and guidance, revealing mechanisms of DNA damage and repair. Here, we analyzed the impact of High-frequency mutated genes on patient prognosis and found that the effect of expression levels of tumor protein p53 (TP53) and lysine methyltransferase 2C (KMT2C) genes remained high throughout the development of GACA, with similar expression patterns. After investigating mutation asymmetry across mutagenic processes, we found that transcriptional asymmetry was dominated by T > G mutations under the influence of transcription couples repair and damage. The apolipoprotein B mRNA editing enzyme catalytic polypeptide like (APOBEC) enzyme that induces mutations during DNA replication has been identified here and we identified a replicative asymmetry, which was dominated by C > A mutations in left-replicating. Strand bias in different mutation classes at transcription factor binding sites and enhancer regions were also confirmed, which implies the important role of non-coding regulatory elements in the occurrence of mutations. This work systematically describes mutational strand asymmetries in specific genomic regions, shedding light on the DNA damage and repair mechanisms underlying somatic mutations in cohorts of GACA patients with gastric cancer.

Enhancing cold resistance in Banana (Musa spp.) through EMS-induced mutagenesis, L-Hyp pressure selection: phenotypic alterations, biomass composition, and transcriptomic insights

BACKGROUND

The cultivation of bananas encounters substantial obstacles, particularly due to the detrimental effects of cold stress on their growth and productivity. A potential remedy that has gained attention is the utilization of ethyl mesylate (EMS)-induced mutagenesis technology, which enables the creation of a genetically varied group of banana mutants. This complex procedure entails subjecting the mutants to further stress screening utilizing L-Hyp in order to identify those exhibiting improved resistance to cold. This study conducted a comprehensive optimization of the screening conditions for EMS mutagenesis and L-Hyp, resulting in the identification of the mutant cm784, which exhibited remarkable cold resistance. Subsequent investigations further elucidated the physiological and transcriptomic responses of cm784 to low-temperature stress.

RESULTS

EMS mutagenesis had a substantial effect on banana seedlings, resulting in modifications in shoot and root traits, wherein a majority of seedlings exhibited delayed differentiation and limited elongation. Notably, mutant leaves displayed altered biomass composition, with starch content exhibiting the most pronounced variation. The application of L-Hyp pressure selection aided in the identification of cold-resistant mutants among seedling-lethal phenotypes. The mutant cm784 demonstrated enhanced cold resistance, as evidenced by improved survival rates and reduced symptoms of chilling injury. Physiological analyses demonstrated heightened activities of antioxidant enzymes and increased proline production in cm784 when subjected to cold stress. Transcriptome analysis unveiled 946 genes that were differentially expressed in cm784, with a notable enrichment in categories related to 'Carbohydrate transport and metabolism' and 'Secondary metabolites biosynthesis, transport, and catabolism'.

CONCLUSION

The present findings provide insights into the molecular mechanisms that contribute to the heightened cold resistance observed in banana mutants. These mechanisms encompass enhanced carbohydrate metabolism and secondary metabolite biosynthesis, thereby emphasizing the adaptive strategies employed to mitigate the detrimental effects induced by cold stress.

Low Spontaneous Mutation Rate in Complex Multicellular Eukaryotes with a Haploid-Diploid Life Cycle

The spontaneous mutation rate µ is a crucial parameter to understand evolution and biodiversity. Mutation rates are highly variable across species, suggesting that µ is susceptible to selection and drift and that species life cycle and life history may impact its evolution. In particular, asexual reproduction and haploid selection are expected to affect the mutation rate, but very little empirical data are available to test this expectation. Here, we sequence 30 genomes of a parent-offspring pedigree in the model brown alga Ectocarpus sp.7, and 137 genomes of an interspecific cross of the closely related brown alga Scytosiphon to have access to the spontaneous mutation rate of representative organisms of a complex multicellular eukaryotic lineage outside animals and plants, and to evaluate the potential impact of life cycle on the mutation rate. Brown algae alternate between a haploid and a diploid stage, both multicellular and free living, and utilize both sexual and asexual reproduction. They are, therefore, excellent models to empirically test expectations of the effect of asexual reproduction and haploid selection on mutation rate evolution. We estimate that Ectocarpus has a base substitution rate of µbs = 4.07 × 10-10 per site per generation, whereas the Scytosiphon interspecific cross had µbs = 1.22 × 10-9. Overall, our estimations suggest that these brown algae, despite being multicellular complex eukaryotes, have unusually low mutation rates. In Ectocarpus, effective population size (Ne) could not entirely explain the low µbs. We propose that the haploid-diploid life cycle, combined with extensive asexual reproduction, may be additional key drivers of the mutation rate in these organisms.

Similar mutation rates but different mutation spectra in moderate and extremely halophilic archaea

Archaea are a major part of Earth's microbiota and extremely diverse. Yet, we know very little about the process of mutation that drives such diversification. To expand beyond previous work with the moderate halophilic archaeal species Haloferax volcanii, we performed a mutation-accumulation experiment followed by whole-genome sequencing in the extremely halophilic archaeon Halobacterium salinarum. Although Hfx. volcanii and Hbt. salinarum have different salt requirements, both species have highly polyploid genomes and similar GC content. We accumulated mutations for an average of 1250 generations in 67 mutation accumulation lines of Hbt. salinarum, and revealed 84 single-base substitutions and 10 insertion-deletion mutations. The estimated base-substitution mutation rate of 3.99 × 10-10 per site per generation or 1.0 × 10-3 per genome per generation in Hbt. salinarum is similar to that reported for Hfx. volcanii (1.2 × 10-3 per genome per generation), but the genome-wide insertion-deletion rate and spectrum of mutations are somewhat dissimilar in these archaeal species. The spectra of spontaneous mutations were AT biased in both archaea, but they differed in significant ways that may be related to differences in the fidelity of DNA replication/repair mechanisms or a simple result of the different salt concentrations.

Very Low Rates of Spontaneous Gene Deletions and Gene Duplications in Dictyostelium discoideum

The study of spontaneous mutation rates has revealed a wide range of heritable point mutation rates across species, but there are comparatively few estimates for large-scale deletion and duplication rates. The handful of studies that have directly calculated spontaneous rates of deletion and duplication using mutation accumulation lines have estimated that genes are duplicated and deleted at orders of magnitude greater rates than the spontaneous point mutation rate. In our study, we tested whether spontaneous gene deletion and gene duplication rates are also high in Dictyostelium discoideum, a eukaryote with among the lowest point mutation rates (2.5 × 10-11 per site per generation) and an AT-rich genome (GC content of 22%). We calculated mutation rates of gene deletions and duplications using whole-genome sequencing data originating from a mutation accumulation experiment and determined the association between the copy number mutations and GC content. Overall, we estimated an average of 3.93 × 10-8 gene deletions and 1.18 × 10-8 gene duplications per gene per generation. While orders of magnitude greater than their point mutation rate, these rates are much lower compared to gene deletion and duplication rates estimated from mutation accumulation lines in other organisms (that are on the order of ~ 10-6 per gene/generation). The deletions and duplications were enriched in regions that were AT-rich even compared to the genomic background, in contrast to our expectations if low GC content was contributing to low mutation rates. The low deletion and duplication mutation rates in D. discoideum compared to other eukaryotes mirror their low point mutation rates, supporting previous work suggesting that this organism has high replication fidelity and effective molecular machinery to avoid the accumulation of mutations in their genome.

Rates and spectra of de novo structural mutations in Chlamydomonas reinhardtii

Genetic variation originates from several types of spontaneous mutation, including single-nucleotide substitutions, short insertions and deletions (indels), and larger structural changes. Structural mutations (SMs) drive genome evolution and are thought to play major roles in evolutionary adaptation, speciation, and genetic disease, including cancers. Sequencing of mutation accumulation (MA) lines has provided estimates of rates and spectra of single-nucleotide and indel mutations in many species, yet the rate of new SMs is largely unknown. Here, we use long-read sequencing to determine the full mutation spectrum in MA lines derived from two strains (CC-1952 and CC-2931) of the green alga Chlamydomonas reinhardtii The SM rate is highly variable between strains and between MA lines, and SMs represent a substantial proportion of all mutations in both strains (CC-1952 6%; CC-2931 12%). The SM spectra differ considerably between the two strains, with almost all inversions and translocations occurring in CC-2931 MA lines. This variation is associated with heterogeneity in the number and type of active transposable elements (TEs), which comprise major proportions of SMs in both strains (CC-1952 22%; CC-2931 38%). In CC-2931, a Crypton and a previously undescribed type of DNA element have caused 71% of chromosomal rearrangements, whereas in CC-1952, a Dualen LINE is associated with 87% of duplications. Other SMs, notably large duplications in CC-2931, are likely products of various double-strand break repair pathways. Our results show that diverse types of SMs occur at substantial rates, and support prominent roles for SMs and TEs in evolution.

Advances in engineering algae for biofuel production

While algae demonstrate potential as a sustainable fuel source, low productivities limit the economic realization of algal biofuels. High-throughput strain engineering, omics-informed genome-scale modeling, and microbiome engineering are key technologies for enabling algal biofuels. High-throughput strain engineering efforts generate improved traits, including high biomass productivity and lipid content, in diverse algal species. Genome-scale models, constructed with the aid of omics data, provide insight into metabolic limitations and guide rational algal strain engineering efforts. As outdoor cultivation systems introduce exogenous organisms, microbiome engineering seeks to eliminate harmful organisms and introduce beneficial species. Optimizing algal biomass production and lipid content using these technologies may overcome the productivity barrier for the commercialization of algal biofuels.

Salt stress alters the spectrum of de novo mutation available to selection during experimental adaptation of Chlamydomonas reinhardtii

The genetic basis of adaptation is driven by both selection and the spectrum of available mutations. Given that the rate of mutation is not uniformly distributed across the genome and varies depending on the environment, understanding the signatures of selection across the genome is aided by first establishing what the expectations of genetic change are from mutation. To determine the interaction between salt stress, selection, and mutation across the genome, we compared mutations observed in a selection experiment for salt tolerance in Chlamydomonas reinhardtii to those observed in mutation accumulation (MA) experiments with and without salt exposure. MA lines evolved under salt stress had a single-nucleotide mutation rate of 1.1 × 10 - 9 $1.1 \times 10^{-9}$ , similar to that of MA lines under standard conditions ( 9.6 × 10 - 10 $9.6 \times 10^{-10}$ ). However, we found that salt stress led to an increased rate of indel mutations, but that many of these mutations were removed under selection. Finally, lines adapted to salt also showed excess clustering of mutations in the genome and the co-expression network, suggesting a role for positive selection in retaining mutations in particular compartments of the genome during the evolution of salt tolerance. Our study shows that characterizing mutation rates and spectra expected under stress helps disentangle the effects of environment and selection during adaptation.

Long-Term Stability of Bacterial Associations in a Microcosm of Ostreococcus tauri (Chlorophyta, Mamiellophyceae)

Phytoplankton-bacteria interactions rule over carbon fixation in the sunlit ocean, yet only a handful of phytoplanktonic-bacteria interactions have been experimentally characterized. In this study, we investigated the effect of three bacterial strains isolated from a long-term microcosm experiment with one Ostreococcus strain (Chlorophyta, Mamiellophyceae). We provided evidence that two Roseovarius strains (Alphaproteobacteria) had a beneficial effect on the long-term survival of the microalgae whereas one Winogradskyella strain (Flavobacteriia) led to the collapse of the microalga culture. Co-cultivation of the beneficial and the antagonistic strains also led to the loss of the microalga cells. Metagenomic analysis of the microcosm is consistent with vitamin B12 synthesis by the Roseovarius strains and unveiled two additional species affiliated to Balneola (Balneolia) and Muricauda (Flavobacteriia), which represent less than 4% of the reads, whereas Roseovarius and Winogradskyella recruit 57 and 39% of the reads, respectively. These results suggest that the low-frequency bacterial species may antagonize the algicidal effect of Winogradskyella in the microbiome of Ostreococcus tauri and thus stabilize the microalga persistence in the microcosm. Altogether, these results open novel perspectives into long-term stability of phytoplankton cultures.

Diversity and Evolution of Mamiellophyceae: Early-Diverging Phytoplanktonic Green Algae Containing Many Cosmopolitan Species

The genomic revolution has bridged a gap in our knowledge about the diversity, biology and evolution of unicellular photosynthetic eukaryotes, which bear very few discriminating morphological features among species from the same genus. The high-quality genome resources available in the class Mamiellophyceae (Chlorophyta) have been paramount to estimate species diversity and screen available metagenomic data to assess the biogeography and ecological niches of different species on a global scale. Here we review the current knowledge about the diversity, ecology and evolution of the Mamiellophyceae and the large double-stranded DNA prasinoviruses infecting them, brought by the combination of genomic and metagenomic analyses, including 26 metabarcoding environmental studies, as well as the pan-oceanic GOS and the Tara Oceans expeditions.

Codon Usage Bias in Phytoplankton

Non-random usage of synonymous codons, known as “codon bias”, has been described in many organisms, from bacteria to Drosophila, but little is known about it in phytoplankton. This phenomenon is thought to be driven by selection for translational efficiency. As the efficacy of selection is proportional to the effective population size, species with large population sizes, such as phytoplankton, are expected to have strong codon bias. To test this, we measured codon bias in 215 strains from Haptophyta, Chlorophyta, Ochrophyta (except diatoms that were studied previously), Dinophyta, Cryptophyta, Ciliophora, unicellular Rhodophyta and Chlorarachniophyta. Codon bias is modest in most groups, despite the astronomically large population sizes of marine phytoplankton. The strength of the codon bias, measured with the effective number of codons, is the strongest in Haptophyta and the weakest in Chlorarachniophyta. The optimal codons are GC-ending in most cases, but several shifts to AT-ending codons were observed (mainly in Ochrophyta and Ciliophora). As it takes a long time to reach a new equilibrium after such shifts, species having AT-ending codons show a lower frequency of optimal codons compared to other species. Genetic diversity, calculated for species with more than three strains sequenced, is modest, indicating that the effective population sizes are many orders of magnitude lower than the astronomically large census population sizes, which helps to explain the modest codon bias in marine phytoplankton. This study represents the first comparative analysis of codon bias across multiple major phytoplankton groups.

The genome-wide rate and spectrum of EMS-induced heritable mutations in the microcrustacean Daphnia: on the prospect of forward genetics

Forward genetic screening using the alkylating mutagen ethyl methanesulfonate (EMS) is an effective method for identifying phenotypic mutants of interest, which can be further genetically dissected to pinpoint the causal genetic mutations. An accurate estimate of the rate of EMS-induced heritable mutations is fundamental for determining the mutant sample size of a screening experiment that aims to saturate all the genes in a genome with mutations. This study examines the genome-wide EMS-induced heritable base-substitutions in three species of the freshwater microcrustacean Daphnia to help guide screening experiments. Our results show that the 10 mM EMS treatment induces base substitutions at an average rate of 1.17 × 10^-6/site/generation across the three species, whereas a significantly higher average mutation rate of 1.75 × 10^-6 occurs at 25 mM. The mutation spectrum of EMS-induced base substitutions at both concentration is dominated by G:C to A:T transitions. Furthermore, we find that female Daphnia exposed to EMS (F₀ individuals) can asexually produce unique mutant offspring (F₁) for at least 3 consecutive broods, suggestive of multiple broods as F₁ mutants. Lastly, we estimate that about 750 F₁s are needed for all genes in the Daphnia genome to be mutated at least once with a 95% probability. We also recommend 4-5 F₂s should be collected from each F₁ mutant through sibling crossing so that all induced mutations could appear in the homozygous state in the F₂ population at 70-80% probability.

Epigenetic modifications affect the rate of spontaneous mutations in a pathogenic fungus

Mutations are the source of genetic variation and the substrate for evolution. Genome-wide mutation rates appear to be affected by selection and are probably adaptive. Mutation rates are also known to vary along genomes, possibly in response to epigenetic modifications, but causality is only assumed. In this study we determine the direct impact of epigenetic modifications and temperature stress on mitotic mutation rates in a fungal pathogen using a mutation accumulation approach. Deletion mutants lacking epigenetic modifications confirm that histone mark H3K27me3 increases whereas H3K9me3 decreases the mutation rate. Furthermore, cytosine methylation in transposable elements (TE) increases the mutation rate 15-fold resulting in significantly less TE mobilization. Also accessory chromosomes have significantly higher mutation rates. Finally, we find that temperature stress substantially elevates the mutation rate. Taken together, we find that epigenetic modifications and environmental conditions modify the rate and the location of spontaneous mutations in the genome and alter its evolutionary trajectory.

Dictyostelium discoideum as a Model to Assess Genome Stability Through DNA Repair

Preserving genome integrity through repair of DNA damage is critical for human health and defects in these pathways lead to a variety of pathologies, most notably cancer. The social amoeba Dictyostelium discoideum is remarkably resistant to DNA damaging agents and genome analysis reveals it contains orthologs of several DNA repair pathway components otherwise limited to vertebrates. These include the Fanconi Anemia DNA inter-strand crosslink and DNA strand break repair pathways. Loss of function of these not only results in malignancy, but also neurodegeneration, immune-deficiencies and congenital abnormalities. Additionally, D. discoideum displays remarkable conservations of DNA repair factors that are targets in cancer and other therapies, including poly(ADP-ribose) polymerases that are targeted to treat breast and ovarian cancers. This, taken together with the genetic tractability of D. discoideum, make it an attractive model to assess the mechanistic basis of DNA repair to provide novel insights into how these pathways can be targeted to treat a variety of pathologies. Here we describe progress in understanding the mechanisms of DNA repair in D. discoideum, and how these impact on genome stability with implications for understanding development of malignancy.

Phytoplankton biodiversity and the inverted paradox

Earth's aquatic food webs are overwhelmingly supported by planktonic microalgae that live in the sunlit water column where only a minimum number of physical niches are readily identifiable. Despite this paucity of environmental differentiation, these "phytoplankton" populations exhibit a rich biodiversity, an observation not easily reconciled with broadly accepted rules of resource-based competitive exclusion. This conundrum is referred to as the "Paradox of the Plankton". Consideration of physical distancing between nutrient depletion zones around individual phytoplankton, however, suggests a competition-neutral resource landscape. Application of neutral theory to the sheer number of phytoplankton in physically-mixed water masses yields a prediction of astronomical biodiversity, suggesting the inverted paradox: Why are there so few phytoplankton species? Here, we introduce a trophic constraint on phytoplankton that, when combined with stochastic principals of ecological drift, predicts only modest levels of diversity in an otherwise competition-neutral landscape. Our "trophic exclusion" principle predicts diversity to be independent of population size and yields a species richness across cell-size classes that is consistent with broad oceanographic survey observations.

Evolutionary Genomics of Sex-Related Chromosomes at the Base of the Green Lineage

Although sex is now accepted as a ubiquitous and ancestral feature of eukaryotes, direct observation of sex is still lacking in most unicellular eukaryotic lineages. Evidence of sex is frequently indirect and inferred from the identification of genes involved in meiosis from whole genome data and/or the detection of recombination signatures from genetic diversity in natural populations. In haploid unicellular eukaryotes, sex-related chromosomes are named mating-type (MTs) chromosomes and generally carry large genomic regions where recombination is suppressed. These regions have been characterized in Fungi and Chlorophyta and determine gamete compatibility and fusion. Two candidate MT+ and MT− alleles, spanning 450–650 kb, have recently been described in Ostreococcus tauri, a marine phytoplanktonic alga from the Mamiellophyceae class, an early diverging branch in the green lineage. Here, we investigate the architecture and evolution of these candidate MT+ and MT− alleles. We analyzed the phylogenetic profile and GC content of MT gene families in eight different genomes whose divergence has been previously estimated at up to 640 Myr, and found evidence that the divergence of the two MT alleles predates speciation in the Ostreococcus genus. Phylogenetic profiles of MT trans-specific polymorphisms in gametologs disclosed candidate MTs in two additional species, and possibly a third. These Mamiellales MT candidates are likely to be the oldest mating-type loci described to date, which makes them fascinating models to investigate the evolutionary mechanisms of haploid sex determination in eukaryotes.

A neutral model for the loss of recombination on sex chromosomes

The loss of recombination between sex chromosomes has occurred repeatedly throughout nature, with important implications for their subsequent evolution. Explanations for this remarkable convergence have generally invoked only adaptive processes (e.g. sexually antagonistic selection); however, there is still little evidence for these hypotheses. Here we propose a model in which recombination on sex chromosomes is lost due to the neutral accumulation of sequence divergence adjacent to (and thus, in linkage disequilibrium with) the sex determiner. Importantly, we include in our model the fact that sequence divergence, in any form, reduces the probability of recombination between any two sequences. Using simulations, we show that, under certain conditions, a region of suppressed recombination arises and expands outwards from the sex-determining locus, under purely neutral processes. Further, we show that the rate and pattern of recombination loss are sensitive to the pre-existing recombination landscape of the genome and to sex differences in recombination rates, with patterns consistent with evolutionary strata emerging under some conditions. We discuss the applicability of these results to natural systems. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

De Novo Mutation Rate Variation and Its Determinants in Chlamydomonas

De novo mutations are central for evolution, since they provide the raw material for natural selection by regenerating genetic variation. However, studying de novo mutations is challenging and is generally restricted to model species, so we have a limited understanding of the evolution of the mutation rate and spectrum between closely related species. Here, we present a mutation accumulation (MA) experiment to study de novo mutation in the unicellular green alga Chlamydomonas incerta and perform comparative analyses with its closest known relative, Chlamydomonas reinhardtii. Using whole-genome sequencing data, we estimate that the median single nucleotide mutation (SNM) rate in C. incerta is μ = 7.6 × 10-10, and is highly variable between MA lines, ranging from μ = 0.35 × 10-10 to μ = 131.7 × 10-10. The SNM rate is strongly positively correlated with the mutation rate for insertions and deletions between lines (r > 0.97). We infer that the genomic factors associated with variation in the mutation rate are similar to those in C. reinhardtii, allowing for cross-prediction between species. Among these genomic factors, sequence context and complexity are more important than GC content. With the exception of a remarkably high C→T bias, the SNM spectrum differs markedly between the two Chlamydomonas species. Our results suggest that similar genomic and biological characteristics may result in a similar mutation rate in the two species, whereas the SNM spectrum has more freedom to diverge.

The spontaneous mutation rate of Drosophila pseudoobscura

The spontaneous mutation rate is a very variable trait that is subject to drift, selection and is sometimes highly plastic. Consequently, its variation between close species, or even between populations from the same species, can be very large. Here, I estimated the spontaneous mutation rate of Drosophila pseudoobscura and Drosophila persimilis crosses to explore the mutation rate variation within the Drosophila genus. All mutation rate estimations in Drosophila varied fourfold, probably explained by the sensitivity of the mutation rate to environmental and experimental conditions. Moreover, I found a very high mutation rate in the hybrid cross between D. pseudoobscura and D. persimilis, in agreement with known elevated mutation rate in hybrids. This mutation rate increase can be explained by heterozygosity and fitness decrease effects in hybrids.

Evolution of Mutation Rate in Astronomically Large Phytoplankton Populations

Genetic diversity is expected to be proportional to population size, yet, there is a well-known, but unexplained lack of genetic diversity in large populations-the "Lewontin's paradox." Larger populations are expected to evolve lower mutation rates, which may help to explain this paradox. Here, we test this conjecture by measuring the spontaneous mutation rate in a ubiquitous unicellular marine phytoplankton species Emiliania huxleyi (Haptophyta) that has modest genetic diversity despite an astronomically large population size. Genome sequencing of E. huxleyi mutation accumulation lines revealed 455 mutations, with an unusual GC-biased mutation spectrum. This yielded an estimate of the per site mutation rate µ = 5.55×10-10 (CI 95%: 5.05×10-10 - 6.09×10-10), which corresponds to an effective population size Ne ∼ 2.7×106. Such a modest Ne is surprising for a ubiquitous and abundant species that accounts for up to 10% of global primary productivity in the oceans. Our results indicate that even exceptionally large populations do not evolve mutation rates lower than ∼10-10 per nucleotide per cell division. Consequently, the extreme disparity between modest genetic diversity and astronomically large population size in the plankton species cannot be explained by an unusually low mutation rate.

Stability across the Whole Nuclear Genome in the Presence and Absence of DNA Mismatch Repair

We describe the contribution of DNA mismatch repair (MMR) to the stability of the eukaryotic nuclear genome as determined by whole-genome sequencing. To date, wild-type nuclear genome mutation rates are known for over 40 eukaryotic species, while measurements in mismatch repair-defective organisms are fewer in number and are concentrated on Saccharomyces cerevisiae and human tumors. Well-studied organisms include Drosophila melanogaster and Mus musculus, while less genetically tractable species include great apes and long-lived trees. A variety of techniques have been developed to gather mutation rates, either per generation or per cell division. Generational rates are described through whole-organism mutation accumulation experiments and through offspring–parent sequencing, or they have been identified by descent. Rates per somatic cell division have been estimated from cell line mutation accumulation experiments, from systemic variant allele frequencies, and from widely spaced samples with known cell divisions per unit of tissue growth. The latter methods are also used to estimate generational mutation rates for large organisms that lack dedicated germlines, such as trees and hyphal fungi. Mechanistic studies involving genetic manipulation of MMR genes prior to mutation rate determination are thus far confined to yeast, Arabidopsis thaliana, Caenorhabditis elegans, and one chicken cell line. A great deal of work in wild-type organisms has begun to establish a sound baseline, but far more work is needed to uncover the variety of MMR across eukaryotes. Nonetheless, the few MMR studies reported to date indicate that MMR contributes 100-fold or more to genome stability, and they have uncovered insights that would have been impossible to obtain using reporter gene assays.

Fitness effects of spontaneous mutations in a warming world

Spontaneous mutations fuel evolutionary processes and differ in consequence, but the consequences depend on the environment. Biophysical considerations of protein thermostability predict that warm temperatures may systematically increase the deleteriousness of mutation. We sought to test whether mutation reduced fitness more when measured in an environment that reflected climate change projections for temperature. We investigated the effects of spontaneous mutations on life history, size, and fitness in 21 mutation accumulation lines and 12 control lines of Daphnia pulex at standard and elevated (+4℃) temperatures. Warmer temperature accelerated life history and reduced body length and clutch sizes. Mutation led to reduced mean clutch sizes and fitness estimates at both temperatures. We found no evidence of a systematic temperature–mutation interaction on trait means, although some lines showed evidence of beneficial mutation at one temperature and deleterious mutation at the other. However, trait variances are also influenced by mutation, and we observed increased variances due to mutation for most traits. For variance of the intrinsic rate of increase and some reproductive traits, we found significant temperature–mutation interactions, with a larger increase due to mutation in the warmer environment. This suggests that selection on new mutations will be more efficient at elevated temperatures.

Streptomycin and nalidixic acid elevate the spontaneous genome-wide mutation rate in Escherichia coli

Since antibiotic resistance is a growing public health problem worldwide, it is important to understand how antibiotics and spontaneous mutations cooperate and shape the genome-wide mutation rate and spectrum. Here, we quantitatively evaluate genome-wide mutational profiles of Escherichia coli after long-term subinhibitory exposure to a broad-spectrum (streptomycin) and a narrow-spectrum antibiotic (nalidixic acid), using a mutation accumulation design combined with whole-genome resequencing of replicate lines as a mutagenicity test. We determined that, while the genome-wide mutation rate is slightly higher in the streptomycin-treated lines compared to the control lines, there is a significant increase in the nalidixic acid-treated lines. Our findings suggest that both broad and narrow-spectrum antibiotics may elevate the mutation rates in E. coli, but mechanisms of action may affect the consequence, thus contribute to accelerating the rate of adaptation and conferring antibiotic resistance.

Estimation of the SNP Mutation Rate in Two Vegetatively Propagating Species of Duckweed

Mutation rate estimates for vegetatively reproducing organisms are rare, despite their frequent occurrence across the tree of life. Here we report mutation rate estimates in two vegetatively reproducing duckweed species, Lemna minor and Spirodela polyrhiza We use a modified approach to estimating mutation rates by taking into account the reduction in mutation detection power that occurs when new individuals are produced from multiple cell lineages. We estimate an extremely low per generation mutation rate in both species of duckweed and note that allelic coverage at de novo mutation sites is very skewed. We also find no substantial difference in mutation rate between mutation accumulation lines propagated under benign conditions and those grown under salt stress. Finally, we discuss the implications of interpreting mutation rate estimates in vegetatively propagating organisms.

A model of autowave self-organization as a hierarchy of active media in the biological evolution

Within the framework of the active media concept, we develop a biophysical model of autowave self-organization which is treated as a hierarchy of active media in the evolution of the biosphere. We also propose a mathematical model of the autowave process of speciation in a flow of mutations for the three main taxonometric groups (prokaryotes, unicellular and multicellular eukaryotes) with a naturally determined lower boundary of living matter (the appearance of prokaryotes) and an open upper boundary for the formation of new species. It is shown that the fluctuation-bifurcation description of the evolution for the formation of new taxonometric groups as a trajectory of transformation of small fluctuations into giant ones adequately reflects the process of self-organization during the formation of taxa. The major concepts of biological evolution, conditions of hierarchy formation as a fundamental manifestation of self-organization and complexity in the evolution of biological systems are considered.

Low Base-Substitution Mutation Rate but High Rate of Slippage Mutations in the Sequence Repeat-Rich Genome of Dictyostelium discoideum

We describe the rate and spectrum of spontaneous mutations for the social amoeba Dictyostelium discoideum, a key model organism in molecular, cellular, evolutionary and developmental biology. Whole-genome sequencing of 37 mutation accumulation lines of D. discoideum after an average of 1,500 cell divisions yields a base-substitution mutation rate of 2.47 × 10⁻¹¹ per site per generation, substantially lower than that of most eukaryotic and prokaryotic organisms, and of the same order of magnitude as in the ciliates Paramecium tetraurelia and Tetrahymena thermophila. Known for its high genomic AT content and abundance of simple sequence repeats, we observe that base-substitution mutations in D. discoideum are highly A/T biased. This bias likely contributes both to the high genomic AT content and to the formation of simple sequence repeats in the AT-rich genome of Dictyostelium discoideum. In contrast to the situation in other surveyed unicellular eukaryotes, indel rates far exceed the base-substitution mutation rate in this organism with a high proportion of 3n indels, particularly in regions without simple sequence repeats. Like ciliates, D. discoideum has a large effective population size, reducing the power of random genetic drift, magnifying the effect of selection on replication fidelity, in principle allowing D. discoideum to evolve an extremely low base-substitution mutation rate.

Single Nucleotide Polymorphism Charting of P. patens Reveals Accumulation of Somatic Mutations During in vitro Culture on the Scale of Natural Variation by Selfing

Introduction: Physcomitrium patens (Hedw.) Mitten (previously known as Physcomitrella patens) was collected by H.L.K. Whitehouse in Gransden Wood (Huntingdonshire, United Kingdom) in 1962 and distributed across the globe starting in 1974. Hence, the Gransden accession has been cultured in vitro in laboratories for half a century. Today, there are more than 13 different pedigrees derived from the original accession. Additionally, accessions from other sites worldwide were collected during the last decades. Methods and Results: In this study, 250 high throughput RNA sequencing (RNA-seq) samples and 25 gDNA samples were used to detect single nucleotide polymorphisms (SNPs). Analyses were performed using five different P. patens accessions and 13 different Gransden pedigrees. SNPs were overlaid with metadata and known phenotypic variations. Unique SNPs defining Gransden pedigrees and accessions were identified and experimentally confirmed. They can be successfully employed for PCR-based identification. Conclusion: We show independent mutations in different Gransden laboratory pedigrees, demonstrating that somatic mutations occur and accumulate during in vitro culture. The frequency of such mutations is similar to those observed in naturally occurring populations. We present evidence that vegetative propagation leads to accumulation of deleterious mutations, and that sexual reproduction purges those. Unique SNP sets for five different P. patens accessions were isolated and can be used to determine individual accessions as well as Gransden pedigrees. Based on that, laboratory methods to easily determine P. patens accessions and Gransden pedigrees are presented.

Synergism between the Black Queen effect and the proteomic constraint on genome size reduction in the photosynthetic picoeukaryotes

The photosynthetic picoeukaryotes (PPEs) comprise a rare example of free-living eukaryotes that have undergone genome reduction. Here, we examine a duality in the process; the proposed driver of genome reduction (the Black Queen hypothesis, BQH), and the resultant impact of genome information loss (the Proteomic Constraint hypothesis, PCH). The BQH predicts that some metabolites may be shared in the open ocean, thus driving loss of redundant metabolic pathways in individual genomes. In contrast, the PCH predicts that as the information content of a genome is reduced, the total mutation load is also reduced, leading to loss of DNA repair genes due to the resulting reduction in selective constraint. Consistent with the BQH, we observe that biosynthetic pathways involved with soluble metabolites such as amino acids and carotenoids are preferentially lost from the PPEs, in contrast to biosynthetic pathways involved with insoluble metabolites, such as lipids, which are retained. Consistent with the PCH, a correlation between proteome size and the number of DNA repair genes, and numerous other informational categories, is observed. While elevated mutation rates resulting from the loss of DNA repair genes have been linked to reduced effective population sizes in intracellular bacteria, this remains to be established. This study shows that in microbial species with large population sizes, an underlying factor in modulating their DNA repair capacity appears to be information content.

Virus-host coexistence in phytoplankton through the genomic lens

Virus-microbe interactions in the ocean are commonly described by "boom and bust" dynamics, whereby a numerically dominant microorganism is lysed and replaced by a virus-resistant one. Here, we isolated a microalga strain and its infective dsDNA virus whose dynamics are characterized instead by parallel growth of both the microalga and the virus. Experimental evolution of clonal lines revealed that this viral production originates from the lysis of a minority of virus-susceptible cells, which are regenerated from resistant cells. Whole-genome sequencing demonstrated that this resistant-susceptible switch involved a large deletion on one chromosome. Mathematical modeling explained how the switch maintains stable microalga-virus population dynamics consistent with their observed growth pattern. Comparative genomics confirmed an ancient origin of this "accordion" chromosome despite a lack of sequence conservation. Together, our results show how dynamic genomic rearrangements may account for a previously overlooked coexistence mechanism in microalgae-virus interactions.

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.

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.

Investigating Evolutionary Rate Variation in Bacteria

Rates of molecular evolution are known to vary between species and across all kingdoms of life. Here, we explore variation in the rate at which bacteria accumulate mutations (accumulation rates) in their natural environments over short periods of time. We have compiled estimates of the accumulation rate for over 34 species of bacteria, the majority of which are pathogens evolving either within an individual host or during outbreaks. Across species, we find that accumulation rates vary by over 3700-fold. We investigate whether accumulation rates are associated to a number potential correlates including genome size, GC content, measures of the natural selection and the time frame over which the accumulation rates were estimated. After controlling for phylogenetic non-independence, we find that the accumulation rate is not significantly correlated to any factor. Furthermore, contrary to previous results, we find that it is not impacted by the time frame of which the estimate was made. However, our study, with only 34 species, is likely to lack power to detect anything but large effects. We suggest that much of the rate variation may be explained by differences between species in the generation time in the wild.

Evolution of Codon Usage Bias in Diatoms

Codon usage bias (CUB)-preferential use of one of the synonymous codons, has been described in a wide range of organisms from bacteria to mammals, but it has not yet been studied in marine phytoplankton. CUB is thought to be caused by weak selection for translational accuracy and efficiency. Weak selection can overpower genetic drift only in species with large effective population sizes, such as Drosophila that has relatively strong CUB, while organisms with smaller population sizes (e.g., mammals) have weak CUB. Marine plankton species tend to have extremely large populations, suggesting that CUB should be very strong. Here we test this prediction and describe the patterns of codon usage in a wide range of diatom species belonging to 35 genera from 4 classes. We report that most of the diatom species studied have surprisingly modest CUB (mean Effective Number of Codons, ENC = 56), with some exceptions showing stronger codon bias (ENC = 44). Modest codon bias in most studied diatom species may reflect extreme disparity between astronomically large census and modest effective population size (N_e), with fluctuations in population size and linked selection limiting long-term N_e and rendering selection for optimal codons less efficient. For example, genetic diversity (pi ~0.02 at silent sites) in Skeletonema marinoi corresponds to N_e of about 10 million individuals, which is likely many orders of magnitude lower than its census size. Still, N_e ~10⁷ should be large enough to make selection for optimal codons efficient. Thus, we propose that an alternative process-frequent changes of preferred codons, may be a more plausible reason for low CUB despite highly efficient selection for preferred codons in diatom populations. The shifts in the set of optimal codons should result in the changes of the direction of selection for codon usage, so the actual codon usage never catches up with the moving target of the optimal set of codons and the species never develop strong CUB. Indeed, we detected strong shifts in preferential codon usage within some diatom genera, with switches between preferentially GC-rich and AT-rich 3^rd codon positions (GC3). For example, GC3 ranges from 0.6 to 1 in most Chaetoceros species, while for Chaetoceros dichaeta GC3 = 0.1. Both variation in selection intensity and mutation spectrum may drive such shifts in codon usage and limit the observed CUB. Our study represents the first genome-wide analysis of CUB in diatoms and the first such analysis for a major phytoplankton group.

Considering the Role of Adaptive Evolution in Models of the Ocean and Climate System

Numerical models have been highly successful in simulating global carbon and nutrient cycles in today's ocean, together with observed spatial and temporal patterns of chlorophyll and plankton biomass at the surface. With this success has come some confidence in projecting the century-scale response to continuing anthropogenic warming. There is also increasing interest in using such models to understand the role of plankton ecosystems in past oceans. However, today's marine environment is the product of billions of years of continual evolution-a process that continues today. In this paper, we address the questions of whether an assumption of species invariance is sufficient, and if not, under what circumstances current model projections might break down. To do this, we first identify the key timescales and questions asked of models. We then review how current marine ecosystem models work and what alternative approaches are available to account for evolution. We argue that for timescales of climate change overlapping with evolutionary timescales, accounting for evolution may to lead to very different projected outcomes regarding the timescales of ecosystem response and associated global biogeochemical cycling. This is particularly the case for past extinction events but may also be true in the future, depending on the eventual degree of anthropogenic disruption. The discipline of building new numerical models that incorporate evolution is also hugely beneficial in itself, as it forces us to question what we know about adaptive evolution, irrespective of its quantitative role in any specific event or environmental changes.

Mutational Landscape of Spontaneous Base Substitutions and Small Indels in Experimental Caenorhabditis elegans Populations of Differing Size

Experimental investigations into the rates and fitness effects of spontaneous mutations are fundamental to our understanding of the evolutionary process. To gain insights into the molecular and fitness consequences of spontaneous mutations, we conducted a mutation accumulation (MA) experiment at varying population sizes in the nematode Caenorhabditis elegans, evolving 35 lines in parallel for 409 generations at three population sizes (N = 1, 10, and 100 individuals). Here, we focus on nuclear SNPs and small insertion/deletions (indels) under minimal influence of selection, as well as their accrual rates in larger populations under greater selection efficacy. The spontaneous rates of base substitutions and small indels are 1.84 (95% C.I. ± 0.14) × 10-9 substitutions and 6.84 (95% C.I. ± 0.97) × 10-10 changes/site/generation, respectively. Small indels exhibit a deletion bias with deletions exceeding insertions by threefold. Notably, there was no correlation between the frequency of base substitutions, nonsynonymous substitutions, or small indels with population size. These results contrast with our previous analysis of mitochondrial DNA mutations and nuclear copy-number changes in these MA lines, and suggest that nuclear base substitutions and small indels are under less stringent purifying selection compared to the former mutational classes. A transition bias was observed in exons as was a near universal base substitution bias toward A/T. Strongly context-dependent base substitutions, where 5'-Ts and 3'-As increase the frequency of A/T → T/A transversions, especially at the boundaries of A or T homopolymeric runs, manifest as higher mutation rates in (i) introns and intergenic regions relative to exons, (ii) chromosomal cores vs. arms and tips, and (iii) germline-expressed genes.

Extreme Lewontin's Paradox in Ubiquitous Marine Phytoplankton Species

Larger populations are expected to have larger genetic diversity. However, as pointed out by Lewontin in 1974, the range of population sizes exceeds the range of genetic diversity by many orders of magnitude (a.k.a. "Lewontin's paradox," LP). The reasons for LP remain obscure. Here, This paper reports an extreme case of LP in astronomically large populations of the ubiquitous unicellular marine phytoplankton species Emiliania huxleyi (Haptophyta)-the species that accounts for 10-20% of primary productivity in the oceans and its blooms are so extensive that they are visible from space. This study demonstrates that despite the wide distribution and enormous population size, the world-wide sample of E. huxleyi strains with sequenced genomes represents a single cohesive species and contains surprisingly limited genetic diversity (π ∼ 0.006 per silent site). The patterns of polymorphism reveal even larger populations in the past, and frequent recombination (ρ ∼ 0.006) throughout the genome, ruling out demographic history and asexual reproduction as possible causes of low polymorphism in E. huxleyi. Natural selection wiping out genetic diversity at linked sites (a.k.a. "genetic draft") must be strong and frequent to account for low polymorphism in E. huxleyi. This study sheds the first light on poorly understood evolutionary genetic processes in astronomically large populations of marine microplankton.

Somatic mutations substantially increase the per-generation mutation rate in the conifer Picea sitchensis

The rates and biological significance of somatic mutations have long been a subject of debate. Somatic mutations in plants are expected to accumulate with vegetative growth and time, yet rates of somatic mutations are unknown for conifers, which can reach exceptional sizes and ages. We investigated somatic mutation rates in the conifer Sitka spruce (Picea sitchensis (Bong.) Carr.) by analyzing a total of 276 Gb of nuclear DNA from the tops and bottoms of 20 old‐growth trees averaging 76 m in height. We estimate a somatic base substitution rate of 2.7 × 10⁻⁸ per base pair within a generation. To date, this is one of the highest estimated per‐generation rates of mutation among eukaryotes, indicating that somatic mutations contribute substantially to the total per‐generation mutation rate in conifers. Nevertheless, as the sampled trees are centuries old, the per‐year rate is low in comparison with nontree taxa. We argue that although somatic mutations raise genetic load in conifers, they generate important genetic variation and enable selection both among cell lineages within individual trees and among offspring.

Old Trade, New Tricks: Insights into the Spontaneous Mutation Process from the Partnering of Classical Mutation Accumulation Experiments with High-Throughput Genomic Approaches

Mutations spawn genetic variation which, in turn, fuels evolution. Hence, experimental investigations into the rate and fitness effects of spontaneous mutations are central to the study of evolution. Mutation accumulation (MA) experiments have served as a cornerstone for furthering our understanding of spontaneous mutations for four decades. In the pregenomic era, phenotypic measurements of fitness-related traits in MA lines were used to indirectly estimate key mutational parameters, such as the genomic mutation rate, new mutational variance per generation, and the average fitness effect of mutations. Rapidly emerging next-generating sequencing technology has supplanted this phenotype-dependent approach, enabling direct empirical estimates of the mutation rate and a more nuanced understanding of the relative contributions of different classes of mutations to the standing genetic variation. Whole-genome sequencing of MA lines bears immense potential to provide a unified account of the evolutionary process at multiple levels-the genetic basis of variation, and the evolutionary dynamics of mutations under the forces of selection and drift. In this review, we have attempted to synthesize key insights into the spontaneous mutation process that are rapidly emerging from the partnering of classical MA experiments with high-throughput sequencing, with particular emphasis on the spontaneous rates and molecular properties of different mutational classes in nuclear and mitochondrial genomes of diverse taxa, the contribution of mutations to the evolution of gene expression, and the rate and stability of transgenerational epigenetic modifications. Future advances in sequencing technologies will enable greater species representation to further refine our understanding of mutational parameters and their functional consequences.

Comparative demography elucidates the longevity of parasitic and symbiotic relationships

Parasitic and symbiotic relationships govern vast nutrient and energy flows, yet controversy surrounds their longevity. Enduring relationships may engender parallel phylogenies among hosts and parasites, but so may ephemeral relationships when parasites colonize related hosts. An understanding of whether symbiont and host populations have grown and contracted in concert would be useful when considering the temporal durability of these relationships. Here, we devised methods to compare demographic histories derived from genomic data. We compared the historical growth of the agent of severe human malaria, Plasmodium falciparum, and its mosquito vector, Anopheles gambiae, to human and primate histories, thereby discerning long-term parallels and anthropogenic population explosions. The growth history of Trichinella spiralis, a zoonotic parasite disseminated by swine, proved regionally specific, paralleling distinctive growth histories for wild boar in Asia and Europe. Parallel histories were inferred for an anemone and its algal symbiont (Exaiptasia pallida and Symbiodinium minutum). Concerted growth in potatoes and the agent of potato blight (Solanum tuberosum and Phytophthora infestans) did not commence until the age of potato domestication. Through these examples, we illustrate the utility of comparative historical demography as a new exploratory tool by which to interrogate the origins and durability of myriad ecological relationships. To facilitate future use of this approach, we introduce a tool called C-PSMC to align and evaluate the similarity of demographic history curves.

Rapidity of Genomic Adaptations to Prasinovirus Infection in a Marine Microalga

Prasinoviruses are large dsDNA viruses commonly found in aquatic systems worldwide, where they can infect and lyse unicellular prasinophyte algae such as Ostreococcus. Host susceptibility is virus strain-specific, but resistance of susceptible Ostreococcus tauri strains to a virulent virus arises frequently. In clonal resistant lines that re-grow, viruses are usually present for many generations, and genes clustered on chromosome 19 show physical rearrangements and differential expression. Here, we investigated changes occurring during the first two weeks after inoculation of the prasinovirus OtV5. By serial dilutions of cultures at the time of inoculation, we estimated the frequency of resistant cells arising in virus-challenged O. tauri cultures to be 10^-3⁻10^-4 of the inoculated population. Re-growing resistant cells were detectable by flow cytometry 3 days post-inoculation (dpi), visible re-greening of cultures occurred by 6 dpi, and karyotypic changes were visually detectable at 8 dpi. Resistant cell lines showed a modified spectrum of host-virus specificities and much lower levels of OtV5 adsorption.

The Mutation Rate and the Age of the Sex Chromosomes in Silene latifolia

Many aspects of sex chromosome evolution are common to both plants and animals [1], but the process of Y chromosome degeneration, where genes on the Y become non-functional over time, may be much slower in plants due to purifying selection against deleterious mutations in the haploid gametophyte [2, 3]. Testing for differences in Y degeneration between the kingdoms has been hindered by the absence of accurate age estimates for plant sex chromosomes. Here, we used genome resequencing to estimate the spontaneous mutation rate and the age of the sex chromosomes in white campion (Silene latifolia). Screening of single nucleotide polymorphisms (SNPs) in parents and 10 F₁ progeny identified 39 de novo mutations and yielded a rate of 7.31 × 10^-9 (95% confidence interval: 5.20 × 10^-9 - 8.00 × 10^-9) mutations per site per haploid genome per generation. Applying this mutation rate to the synonymous divergence between homologous X- and Y-linked genes (gametologs) gave age estimates of 11.00 and 6.32 million years for the old and young strata, respectively. Based on SNP segregation patterns, we inferred which genes were Y-linked and found that at least 47% are already dysfunctional. Applying our new estimates for the age of the sex chromosomes indicates that the rate of Y degeneration in S. latifolia is nearly 2-fold slower when compared to animal sex chromosomes of a similar age. Our revised estimates support Y degeneration taking place more slowly in plants, a discrepancy that may be explained by differences in the life cycles of animals and plants.

Muver, a computational framework for accurately calling accumulated mutations

BACKGROUND

Identification of mutations from next-generation sequencing data typically requires a balance between sensitivity and accuracy. This is particularly true of DNA insertions and deletions (indels), that can impart significant phenotypic consequences on cells but are harder to call than substitution mutations from whole genome mutation accumulation experiments. To overcome these difficulties, we present muver, a computational framework that integrates established bioinformatics tools with novel analytical methods to generate mutation calls with the extremely low false positive rates and high sensitivity required for accurate mutation rate determination and comparison.

RESULTS

Muver uses statistical comparison of ancestral and descendant allelic frequencies to identify variant loci and assigns genotypes with models that include per-sample assessments of sequencing errors by mutation type and repeat context. Muver identifies maximally parsimonious mutation pathways that connect these genotypes, differentiating potential allelic conversion events and delineating ambiguities in mutation location, type, and size. Benchmarking with a human gold standard father-son pair demonstrates muver's sensitivity and low false positive rates. In DNA mismatch repair (MMR) deficient Saccharomyces cerevisiae, muver detects multi-base deletions in homopolymers longer than the replicative polymerase footprint at rates greater than predicted for sequential single-base deletions, implying a novel multi-repeat-unit slippage mechanism.

CONCLUSIONS

Benchmarking results demonstrate the high accuracy and sensitivity achieved with muver, particularly for indels, relative to available tools. Applied to an MMR-deficient Saccharomyces cerevisiae system, muver mutation calls facilitate mechanistic insights into DNA replication fidelity.

Population genomics of picophytoplankton unveils novel chromosome hypervariability

Tiny photosynthetic microorganisms that form the picoplankton (between 0.3 and 3 μm in diameter) are at the base of the food web in many marine ecosystems, and their adaptability to environmental change hinges on standing genetic variation. Although the genomic and phenotypic diversity of the bacterial component of the oceans has been intensively studied, little is known about the genomic and phenotypic diversity within each of the diverse eukaryotic species present. We report the level of genomic diversity in a natural population of Ostreococcus tauri (Chlorophyta, Mamiellophyceae), the smallest photosynthetic eukaryote. Contrary to the expectations of clonal evolution or cryptic species, the spectrum of genomic polymorphism observed suggests a large panmictic population (an effective population size of 1.2 × 10⁷) with pervasive evidence of sexual reproduction. De novo assemblies of low-coverage chromosomes reveal two large candidate mating-type loci with suppressed recombination, whose origin may pre-date the speciation events in the class Mamiellophyceae. This high genetic diversity is associated with large phenotypic differences between strains. Strikingly, resistance of isolates to large double-stranded DNA viruses, which abound in their natural environment, is positively correlated with the size of a single hypervariable chromosome, which contains 44 to 156 kb of strain-specific sequences. Our findings highlight the role of viruses in shaping genome diversity in marine picoeukaryotes.