Recombination Rate Variation and Infrequent Sex Influence Genetic Diversity in Chlamydomonas reinhardtii

Male-biased recombination at chromosome ends in a songbird revealed by precisely mapping crossover positions

Recombination plays a crucial role in evolution by generating novel haplotypes and disrupting linkage between genes, thereby enhancing the efficiency of selection. Here, we analyze the genomes of 12 great reed warblers (Acrocephalus arundinaceus) in a 3-generation pedigree to identify precise crossover positions along the chromosomes. We located more than 200 crossovers and found that these were highly concentrated toward the telomeric ends of the chromosomes. Apart from this major pattern in the recombination landscape, we found significantly higher frequencies of crossovers in genic compared with intergenic regions, and in exons compared with introns. Moreover, while the number of recombination events was similar between the sexes, the crossovers were located significantly closer to the ends of paternal compared with maternal chromosomes. In conclusion, our study of the great reed warbler revealed substantial variation in crossover frequencies within chromosomes, with a distinct bias toward the sub-telomeric regions, particularly on the paternal side. These findings emphasize the importance of thoroughly screening the entire length of chromosomes to characterize the recombination landscape and uncover potential sex-biases in recombination.

Population bottlenecks and sexual recombination shape diatom microevolution

Diatoms are single-celled organisms that contribute approximately 20% of the global primary production and play a crucial role in biogeochemical cycles and trophic chains. Despite their ecological importance, our knowledge of microevolution is limited. We developed a model using the SLiM evolutionary framework to address this knowledge gap. As a reference, we used the diatom Pseudo-nitzschia multistriata, which has been extensively studied in the Gulf of Naples. Our model recapitulates what we observe in natural populations, with microevolutionary processes that occur annually during a three-stage bloom phase. Interestingly, we found that non-bloom phases allow the population to maintain sex-generated diversity produced during blooms. This finding suggests that non-bloom phases are critical to counteract bloom-related pressures and mitigate genetic divergence at the species level. Moreover, our model showed that despite the consistent genetic differentiation during bloom phases, the population tends to return to pre-bloom states. While our model is limited to neutral dynamics, our study provides valuable insights into diatoms' microevolution, paving the way to explore the ecological implications of the life history dynamics of these organisms.

Multiannual patterns of genetic structure and mating type ratios highlight the complex bloom dynamics of a marine planktonic diatom

Understanding the genetic structure of populations and the processes responsible for its spatial and temporal dynamics is vital for assessing species' adaptability and survival in changing environments. We investigate the genetic fingerprinting of blooming populations of the marine diatom Pseudo-nitzschia multistriata in the Gulf of Naples (Mediterranean Sea) from 2008 to 2020. Strains were genotyped using microsatellite fingerprinting and natural samples were also analysed with Microsatellite Pool-seq Barcoding based on Illumina sequencing of microsatellite loci. Both approaches revealed a clonal expansion event in 2013 and a more stable genetic structure during 2017-2020 compared to previous years. The identification of a mating type (MT) determination gene allowed to assign MT to strains isolated over the years. MTs were generally at equilibrium with two notable exceptions, including the clonal bloom of 2013. The populations exhibited linkage equilibrium in most blooms, indicating that sexual reproduction leads to genetic homogenization. Our findings show that P. multistriata blooms exhibit a dynamic genetic and demographic composition over time, most probably determined by deeper-layer cell inocula. Occasional clonal expansions and MT imbalances can potentially affect the persistence and ecological success of planktonic diatoms.

Population Genomics of the Facultatively Sexual Liverwort Marchantia polymorpha

The population genomics of facultatively sexual organisms are understudied compared with their abundance across the tree of life. We explore patterns of genetic diversity in two subspecies of the facultatively sexual liverwort Marchantia polymorpha using samples from across Southern Ontario, Canada. Despite the ease with which M. polymorpha should be able to propagate asexually, we find no evidence of strictly clonal descent among our samples and little to no signal of isolation by distance. Patterns of identity-by-descent tract sharing further showed evidence of recent recombination and close relatedness between geographically distant isolates, suggesting long distance gene flow and at least a modest frequency of sexual reproduction. However, the M. polymorpha genome contains overall very low levels of nucleotide diversity and signs of inefficient selection evidenced by a relatively high fraction of segregating deleterious variants. We interpret these patterns as possible evidence of the action of linked selection and a small effective population size due to past generations of asexual propagation. Overall, the M. polymorpha genome harbors signals of a complex history of both sexual and asexual reproduction.

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.

Variation in fine-scale recombination rate in temperature-evolved Drosophila melanogaster populations in response to selection

Meiotic recombination plays a critical evolutionary role in maintaining fitness in response to selective pressures due to changing environments. Variation in recombination rate has been observed amongst and between species and populations and within genomes across numerous taxa. Studies have demonstrated a link between changes in recombination rate and selection, but the extent to which fine-scale recombination rate varies between evolved populations during the evolutionary period in response to selection is under active research. Here, we utilize a set of 3 temperature-evolved Drosophila melanogaster populations that were shown to have diverged in several phenotypes, including recombination rate, based on the temperature regime in which they evolved. Using whole-genome sequencing data from these populations, we generated linkage disequilibrium-based fine-scale recombination maps for each population. With these maps, we compare recombination rates and patterns among the 3 populations and show that they have diverged at fine scales but are conserved at broader scales. We further demonstrate a correlation between recombination rates and genomic variation in the 3 populations. Lastly, we show variation in localized regions of enhanced recombination rates, termed warm spots, between the populations with these warm spots and associated genes overlapping areas previously shown to have diverged in the 3 populations due to selection. These data support the existence of recombination modifiers in these populations which are subject to selection during evolutionary change.

Recombination-aware phylogeographic inference using the structured coalescent with ancestral recombination

Movement of individuals between populations or demes is often restricted, especially between geographically isolated populations. The structured coalescent provides an elegant theoretical framework for describing how movement between populations shapes the genealogical history of sampled individuals and thereby structures genetic variation within and between populations. However, in the presence of recombination an individual may inherit different regions of their genome from different parents, resulting in a mosaic of genealogical histories across the genome, which can be represented by an Ancestral Recombination Graph (ARG). In this case, different genomic regions may have different ancestral histories and so different histories of movement between populations. Recombination therefore poses an additional challenge to phylogeographic methods that aim to reconstruct the movement of individuals from genealogies, although also a potential benefit in that different loci may contain additional information about movement. Here, we introduce the Structured Coalescent with Ancestral Recombination (SCAR) model, which builds on recent approximations to the structured coalescent by incorporating recombination into the ancestry of sampled individuals. The SCAR model allows us to infer how the migration history of sampled individuals varies across the genome from ARGs, and improves estimation of key population genetic parameters such as population sizes, recombination rates and migration rates. Using the SCAR model, we explore the potential and limitations of phylogeographic inference using full ARGs. We then apply the SCAR to lineages of the recombining fungus Aspergillus flavus sampled across the United States to explore patterns of recombination and migration across the genome.

How Can We Resolve Lewontin's Paradox?

We discuss the genetic, demographic, and selective forces that are likely to be at play in restricting observed levels of DNA sequence variation in natural populations to a much smaller range of values than would be expected from the distribution of census population sizes alone-Lewontin's Paradox. While several processes that have previously been strongly emphasized must be involved, including the effects of direct selection and genetic hitchhiking, it seems unlikely that they are sufficient to explain this observation without contributions from other factors. We highlight a potentially important role for the less-appreciated contribution of population size change; specifically, the likelihood that many species and populations may be quite far from reaching the relatively high equilibrium diversity values that would be expected given their current census sizes.

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.

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.

Comparative genomics of Chlamydomonas

Despite its role as a reference organism in the plant sciences, the green alga Chlamydomonas reinhardtii entirely lacks genomic resources from closely related species. We present highly contiguous and well-annotated genome assemblies for three unicellular C. reinhardtii relatives: Chlamydomonas incerta, Chlamydomonas schloesseri, and the more distantly related Edaphochlamys debaryana. The three Chlamydomonas genomes are highly syntenous with similar gene contents, although the 129.2 Mb C. incerta and 130.2 Mb C. schloesseri assemblies are more repeat-rich than the 111.1 Mb C. reinhardtii genome. We identify the major centromeric repeat in C. reinhardtii as a LINE transposable element homologous to Zepp (the centromeric repeat in Coccomyxa subellipsoidea) and infer that centromere locations and structure are likely conserved in C. incerta and C. schloesseri. We report extensive rearrangements, but limited gene turnover, between the minus mating type loci of these Chlamydomonas species. We produce an eight-species core-Reinhardtinia whole-genome alignment, which we use to identify several hundred false positive and missing genes in the C. reinhardtii annotation and >260,000 evolutionarily conserved elements in the C. reinhardtii genome. In summary, these resources will enable comparative genomics analyses for C. reinhardtii, significantly extending the analytical toolkit for this emerging model system.

Switching environments, synchronous sex, and the evolution of mating types

While facultative sex is common in sexually reproducing species, for reasons of tractability most mathematical models assume that such sex is asynchronous in the population. In this paper, we develop a model of switching environments to instead capture the effect of an entire population transitioning synchronously between sexual and asexual modes of reproduction. We use this model to investigate the evolution of the number of self-incompatible mating types in finite populations, which empirically can range from two to thousands. When environmental switching is fast, we recover the results of earlier studies that implicitly assumed populations were engaged in asynchronous sexual reproduction. However when the environment switches slowly, we see deviations from previous asynchronous theory, including a lower number of mating types at equilibrium and bimodality in the stationary distribution of mating types. We provide analytic approximations for both the fast and slow switching regimes, as well as a numerical scheme based on the Kolmogorov equations for the system to quickly evaluate the model dynamics at intermediate parameters. Our approach exploits properties of integer partitions in number theory. We also demonstrate how additional biological processes such as selective sweeps can be accounted for in this switching environment framework, showing that beneficial mutations can further erode mating type diversity in synchronous facultatively sexual populations.