The distribution of fitness effects of spontaneous mutations in Chlamydomonas reinhardtii inferred using frequency changes under experimental evolution.

The distribution of fitness effects (DFE) for new mutations is fundamental for many aspects of population and quantitative genetics. In this study, we have inferred the DFE in the single-celled alga Chlamydomonas reinhardtii by estimating changes in the frequencies of 254 spontaneous mutations under experimental evolution and equating the frequency changes of linked mutations with their selection coefficients. We generated seven populations of recombinant haplotypes by crossing seven independently derived mutation accumulation lines carrying an average of 36 mutations in the haploid state to a mutation-free strain of the same genotype. We then allowed the populations to evolve under natural selection in the laboratory by serial transfer in liquid culture. We observed substantial and repeatable changes in the frequencies of many groups of linked mutations, and, surprisingly, as many mutations were observed to increase as decrease in frequency. Mutation frequencies were highly repeatable among replicates, suggesting that selection was the cause of the observed allele frequency changes. We developed a Bayesian Monte Carlo Markov Chain method to infer the DFE. This computes the likelihood of the observed distribution of changes of frequency, and obtains the posterior distribution of the selective effects of individual mutations, while assuming a two-sided gamma distribution of effects. We infer that the DFE is a highly leptokurtic distribution, and that approximately equal proportions of mutations have positive and negative effects on fitness. This result is consistent with what we have observed in previous work on a different C. reinhardtii strain, and suggests that a high fraction of new spontaneously arisen mutations are advantageous in a simple laboratory environment.

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.

Patterns of population structure and complex haplotype sharing among field isolates of the green alga Chlamydomonas reinhardtii.

The nature of population structure in microbial eukaryotes has long been debated. Competing models have argued that microbial species are either ubiquitous, with high dispersal and low rates of speciation, or that for many species gene flow between populations is limited, resulting in evolutionary histories similar to those of macroorganisms. However, population genomic approaches have seldom been applied to this question. Here, we analyse whole-genome resequencing data for all 36 confirmed field isolates of the green alga Chlamydomonas reinhardtii. At a continental scale, we report evidence for putative allopatric divergence, between both North American and Japanese isolates, and two highly differentiated lineages within N. America. Conversely, at a local scale within the most densely sampled lineage, we find little evidence for either spatial or temporal structure. Taken together with evidence for ongoing admixture between the two N. American lineages, this lack of structure supports a role for substantial dispersal in C. reinhardtii and implies that between-lineage differentiation may be maintained by reproductive isolation and/or local adaptation. Our results therefore support a role for allopatric divergence in microbial eukaryotes, while also indicating that species may be ubiquitous at local scales. Despite the high genetic diversity observed within the most well-sampled lineage, we find that pairs of isolates share on average ~9% of their genomes in long haplotypes, even when isolates were sampled decades apart and from different locations. This proportion is several orders of magnitude higher than the Wright-Fisher expectation, raising many further questions concerning the evolutionary genetics of C. reinhardtii and microbial eukaryotes generally.

Inferring the distribution of fitness effects of spontaneous mutations in Chlamydomonas reinhardtii.

Spontaneous mutations are the source of new genetic variation and are thus central to the evolutionary process. In molecular evolution and quantitative genetics, the nature of genetic variation depends critically on the distribution of effects of mutations on fitness and other quantitative traits. Spontaneous mutation accumulation (MA) experiments have been the principal approach for investigating the overall rate of occurrence and cumulative effect of mutations but have not allowed the phenotypic effects of individual mutations to be studied directly. Here, we crossed MA lines of the green alga Chlamydomonas reinhardtii with its unmutated ancestral strain to create haploid recombinant lines, each carrying an average of 50% of the accumulated mutations in a large number of combinations. With the aid of the genome sequences of the MA lines, we inferred the genotypes of the mutations, assayed their growth rate as a measure of fitness, and inferred the distribution of fitness effects (DFE) using a Bayesian mixture model. We infer that the DFE is highly leptokurtic (L-shaped). Of mutations with absolute fitness effects exceeding 1%, about one-sixth increase fitness in the laboratory environment. The inferred distribution of effects for deleterious mutations is consistent with a strong role for nearly neutral evolution. Specifically, such a distribution predicts that nucleotide variation and genetic variation for quantitative traits will be insensitive to change in the effective population size.

Signatures of natural selection in abiotic stress-responsive genes of Solanum chilense.

Environmental conditions are strong selective forces, which may influence adaptation and speciation. The wild tomato species Solanum chilense, native to South America, is exposed to a range of abiotic stress factors. To identify signatures of natural selection and local adaptation, we analysed 16 genes involved in the abiotic stress response and compared the results to a set of reference genes in 23 populations across the entire species range. The abiotic stress-responsive genes are characterized by elevated nonsynonymous nucleotide diversity and divergence. We detected signatures of positive selection in several abiotic stress-responsive genes on both the population and species levels. Local adaptation to abiotic stresses is particularly apparent at the boundary of the species distribution in populations from coastal low-altitude and mountainous high-altitude regions.

Fitness change in relation to mutation number in spontaneous mutation accumulation lines of Chlamydomonas reinhardtii.

Although all genetic variation ultimately stems from mutations, their properties are difficult to study directly. Here, we used multiple mutation accumulation (MA) lines derived from five genetic backgrounds of the green algae Chlamydomonas reinhardtii that have been previously subjected to whole genome sequencing to investigate the relationship between the number of spontaneous mutations and change in fitness from a nonevolved ancestor. MA lines were on average less fit than their ancestors and we detected a significantly negative correlation between the change in fitness and the total number of accumulated mutations in the genome. Likewise, the number of mutations located within coding regions significantly and negatively impacted MA line fitness. We used the fitness data to parameterize a maximum likelihood model to estimate discrete categories of mutational effects, and found that models containing one to two mutational effect categories (one neutral and one deleterious category) fitted the data best. However, the best-fitting mutational effects models were highly dependent on the genetic background of the ancestral strain.

Adaptation to low temperatures in the wild tomato species Solanum chilense.

Molecular adaptation to abiotic stresses in plants is a complex process based mainly on the modifications of gene transcriptional activity and the alteration of protein-protein interactions. We used a combination of population genetic, comparative transcriptomic and plant physiology approaches to investigate the mechanisms of adaptation to low temperatures in Solanum chilense populations distributed along Andean altitudinal gradients. We found that plants from all populations have high chilling tolerance, which does not correlate with temperatures in their native habitats. In contrast, tolerance to freezing shows a significant association with altitude and temperature variables. We also observed the differences in expression patterns of cold-response genes between plants from high- and low-altitude populations. These results suggest that genetic adaptations to low temperatures evolved in high-altitude populations of S. chilense. At the transcriptional level, these adaptations may include high levels of constitutive expression of the genes encoding ICE1, the key transcription factor of the cold signalling pathway, and chloroplast ω-3 fatty acid desaturase FAD7. At the sequence level, a signature of selection associated with the adaptation to high altitudes was detected at the C-terminal part of ICE1 encoding the ACT regulatory domain.

North-South Colonization Associated with Local Adaptation of the Wild Tomato Species Solanum chilense.

After colonization population sizes may vary across the species range depending on environmental conditions and following colonizations. An interesting question is whether local adaptation occurs more frequently in large ancestral populations or in small derived populations. A higher number of new mutations and a lower effect of genetic drift should favor selection in large populations, whereas small derived populations may require an initial local adaptation event to facilitate the colonization of new habitats. Wild tomatoes are native to a broad range of different habitats characterized by variable abiotic conditions in South America, and represent an ideal system to study this interplay between demography and natural selection. Population genetic analyses and statistical inference of past demography were conducted on pooled-sequencing data from 30 genes (8,080 single nucleotide polymorphisms) from an extensive sampling of 23 Solanum chilense populations over Chile and Peru. We reveal first a north-south colonization associated with relaxed purifying selection in the south as shown by a decrease of genetic variation and an increasing proportion of nonsynonymous polymorphism from north to south, and population substructure with at least four genetic groups. Second, we uncover a dual picture of adaptation consisting of 1) a decreasing proportion of adaptive amino acid substitutions from north to south suggesting that adaptation is favored in large populations, whereas 2) signatures of local adaptation predominantly occur in the smaller populations from the marginal ranges in the south.