Sylvain Glémin

ADAPTATION AND MALADAPTATION IN SELFING AND OUTCROSSING SPECIES: NEW MUTATIONS VERSUS STANDING VARIATION

Evolution of selfing from outcrossing recurrently occurred in many lineages, especially in flowering plants. Evolution of selfing induces dramatic changes in the population genetics functioning but its consequences on the dynamics of adaptation have been overlooked. We studied a simple one‐locus model of adaptation where a population experiences an environmental change at a given time. We first determined the effect of the mating system on the genetic bases and the speed of adaptation, focusing on the dominance of beneficial mutations and the respective part of standing variation and new mutations. Then, we assumed that the environmental change is associated with population decline to determine the effect of the mating system on the probability of population extinction. Extending previous results, we found that adaptation is more efficient and extinction less likely in outcrossers when beneficial mutations are dominant and codominant and when standing variation plays a significant role in adaptation. However, given adaptation does occur, it is usually more rapid in selfers than in outcrossers. Our results bear implications for the evolution of the selfing syndrome, the dynamics of the domestication process, and the dead‐end hypothesis that posits that selfing lineages are doomed to extinction on the long run.

Patterns and Evolution of Nucleotide Landscapes in Seed Plants

Analysis of EST sequence data in more than 200 species reveals unexpected continuous variations in GC content in seed plant genomes, with several independent enrichment episodes from GC-poor and homogeneous genomes to more derived GC-rich and highly heterogeneous ones. It reveals the possible role of GC-biased gene conversion, a recombination-associated process favoring G and C bases. Nucleotide landscapes, which are the way base composition is distributed along a genome, strongly vary among species. The underlying causes of these variations have been much debated. Though mutational bias and selection were initially invoked, GC-biased gene conversion (gBGC), a recombination-associated process favoring the G and C over A and T bases, is increasingly recognized as a major factor. As opposed to vertebrates, evolution of GC content is less well known in plants. Most studies have focused on the GC-poor and homogeneous Arabidopsis thaliana genome and the much more GC-rich and heterogeneous rice (Oryza sativa) genome and have often been generalized as a dicot/monocot dichotomy. This vision is clearly phylogenetically biased and does not allow understanding the mechanisms involved in GC content evolution in plants. To tackle these issues, we used EST data from more than 200 species and provided the most comprehensive description of gene GC content across the seed plant phylogeny so far available. As opposed to the classically assumed dicot/monocot dichotomy, we found continuous variations in GC content from the probably ancestral GC-poor and homogeneous genomes to the more derived GC-rich and highly heterogeneous ones, with several independent enrichment episodes. Our results suggest that gBGC could play a significant role in the evolution of GC content in plant genomes.

Male‐specific DNA markers provide genetic evidence of an XY chromosome system, a recombination arrest and allow the tracing of paternal lineages in date palm

Whether sex chromosomes are differentiated is an important aspect of our knowledge of dioecious plants, such as date palm (Phoenix dactylifera). In this crop plant, the female individuals produce dates, and are thus the more valuable sex. However, there is no way to identify the sex of date palm plants before reproductive age, and the sex-determining mechanism is still unclear. To identify sex-linked microsatellite markers, we surveyed a set of 52 male and 55 female genotypes representing the geographical diversity of the species. We found three genetically linked loci that are heterozygous only in males. Male-specific alleles allowed us to identify the gender in 100% of individuals. These results confirm the existence of an XY chromosomal system with a nonrecombining XY-like region in the date palm genome. The distribution of Y haplotypes in western and eastern haplogroups allowed us to trace two male ancestral paternal lineages that account for all known Y diversity in date palm. The very low diversity associated with Y haplotypes is consistent with clonal paternal transmission of a nonrecombining male-determining region. Our results establish the date palm as a biological model with one of the most ancient sex chromosomes in flowering plants.

Transcriptome population genomics reveals severe bottleneck and domestication cost in the African rice (#Oryza glaberrima#)

The African cultivated rice (Oryza glaberrima) was domesticated in West Africa 3000 years ago. Although less cultivated than the Asian rice (O. sativa), O. glaberrima landraces often display interesting adaptation to rustic environment (e.g. drought). Here, using RNA‐seq technology, we were able to compare more than 12 000 transcripts between 9 O. glaberrima, 10 wild O. barthii and one O. meridionalis individuals. With a synonymous nucleotide diversity πs = 0.0006 per site, O. glaberrima appears as the least genetically diverse crop grass ever documented. Using approximate Bayesian computation, we estimated that O. glaberrima experienced a severe bottleneck during domestication. This demographic scenario almost fully accounts for the pattern of genetic diversity across O. glaberrima genome as we detected very few outliers regions where positive selection may have further impacted genetic diversity. Moreover, the large excess of derived nonsynonymous substitution that we detected suggests that the O. glaberrima population suffered from the ‘cost of domestication’. In addition, we used this genome‐scale data set to demonstrate that (i) O. barthii genetic diversity is positively correlated with recombination rate and negatively with gene density, (ii) expression level is negatively correlated with evolutionary constraint, and (iii) one region on chromosome 5 (position 4–6 Mb) exhibits a clear signature of introgression with a yet unidentified Oryza species. This work represents the first genome‐wide survey of the African rice genetic diversity and paves the way for further comparison between the African and the Asian rice, notably regarding the genetics underlying domestication traits.

Quantification of GC-biased gene conversion in the human genome

Much evidence indicates that GC-biased gene conversion (gBGC) has a major impact on the evolution of mammalian genomes. However, a detailed quantification of the process is still lacking. The strength of gBGC can be measured from the analysis of derived allele frequency spectra (DAF), but this approach is sensitive to a number of confounding factors. In particular, we show by simulations that the inference is pervasively affected by polymorphism polarization errors and by spatial heterogeneity in gBGC strength. We propose a new general method to quantify gBGC from DAF spectra, incorporating polarization errors, taking spatial heterogeneity into account, and jointly estimating mutation bias. Applying it to human polymorphism data from the 1000 Genomes Project, we show that the strength of gBGC does not differ between hypermutable CpG sites and non-CpG sites, suggesting that in humans gBGC is not caused by the base-excision repair machinery. Genome-wide, the intensity of gBGC is in the nearly neutral area. However, given that recombination occurs primarily within recombination hotspots, 1%-2% of the human genome is subject to strong gBGC. On average, gBGC is stronger in African than in non-African populations, reflecting differences in effective population sizes. However, due to more heterogeneous recombination landscapes, the fraction of the genome affected by strong gBGC is larger in non-African than in African populations. Given that the location of recombination hotspots evolves very rapidly, our analysis predicts that, in the long term, a large fraction of the genome is affected by short episodes of strong gBGC.

Extreme recombination frequencies shape genome variation and evolution in the honeybee, Apis mellifera.

Meiotic recombination is a fundamental cellular process, with important consequences for evolution and genome integrity. However, we know little about how recombination rates vary across the genomes of most species and the molecular and evolutionary determinants of this variation. The honeybee, Apis mellifera, has extremely high rates of meiotic recombination, although the evolutionary causes and consequences of this are unclear. Here we use patterns of linkage disequilibrium in whole genome resequencing data from 30 diploid honeybees to construct a fine-scale map of rates of crossing over in the genome. We find that, in contrast to vertebrate genomes, the recombination landscape is not strongly punctate. Crossover rates strongly correlate with levels of genetic variation, but not divergence, which indicates a pervasive impact of selection on the genome. Germ-line methylated genes have reduced crossover rate, which could indicate a role of methylation in suppressing recombination. Controlling for the effects of methylation, we do not infer a strong association between gene expression patterns and recombination. The site frequency spectrum is strongly skewed from neutral expectations in honeybees: rare variants are dominated by AT-biased mutations, whereas GC-biased mutations are found at higher frequencies, indicative of a major influence of GC-biased gene conversion (gBGC), which we infer to generate an allele fixation bias 5 – 50 times the genomic average estimated in humans. We uncover further evidence that this repair bias specifically affects transitions and favours fixation of CpG sites. Recombination, via gBGC, therefore appears to have profound consequences on genome evolution in honeybees and interferes with the process of natural selection. These findings have important implications for our understanding of the forces driving molecular evolution.

GC-Biased Gene Conversion Impacts Ribosomal DNA Evolution in Vertebrates, Angiosperms and Other Eukaryotes

Ribosomal DNA (rDNA) is one of the most conserved genes in eukaryotes. The multiples copies of rDNA in the genome evolve in a concerted manner, through unequal crossing over and/or gene conversion, two mechanisms related to homologous recombination. Recombination increases local GC content in several organisms through a process known as GC-biased gene conversion (gBGC). gBGC has been well characterized in mammals, birds, and grasses, but its phylogenetic distribution across the tree of life is poorly understood. Here, we test the hypothesis that recombination affects the evolution of base composition in 18S rDNA and examine the reliability of this thoroughly studied molecule as a marker of gBGC in eukaryotes. Phylogenetic analyses of 18S rDNA in vertebrates and angiosperms reveal significant heterogeneity in the evolution of base composition across both groups. Mammals, birds, and grasses experience increases in the GC content of the 18S rDNA, consistent with previous genome-wide analyses. In addition, we observe increased GC contents in Ostariophysi ray-finned fishes and commelinid monocots (i.e., the clade including grasses), suggesting that the genomes of these two groups have been affected by gBGC. Polymorphism analyses in rDNA confirm that gBGC, not mutation bias, is the most plausible explanation for these patterns. We also find that helix and loop sites of the secondary structure of ribosomal RNA do not evolve at the same pace: loops evolve faster than helices, whereas helices are GC richer than loops. We extend analyses to major lineages of eukaryotes and suggest that gBGC might have also affected base composition in Giardia (Diplomonadina), nudibranch gastropods (Mollusca), and Asterozoa (Echinodermata).

GC content evolution in coding regions of angiosperm genomes: a unifying hypothesis

In angiosperms (as in other species), GC content varies along and between genes, within a genome, and between genomes of different species, but the reason for this distribution is still an open question. Grass genomes are particularly intriguing because they exhibit a strong bimodal distribution of genic GC content and a sharp 5-3 decreasing GC content gradient along most genes. Here, we propose a unifying model to explain the main patterns of GC content variation at the gene and genome scale. We argue that GC content patterns could be mainly determined by the interactions between gene structure, recombination patterns, and GC-biased gene conversion. Recent studies on fine-scale recombination maps in angiosperms support this hypothesis and previous results also fit this model. We propose that our model could be used as a null hypothesis to search for additional forces that affect GC content in angiosperms.

Hitchhiking of Deleterious Alleles and the Cost of Adaptation in Partially Selfing Species

Self-fertilization is generally seen to be disadvantageous in the long term. It increases genetic drift, which subsequently reduces polymorphism and the efficiency of selection, which also challenges adaptation. However, high selfing rates can increase the fixation probability of recessive beneficial mutations, but existing theory has generally not accounted for the effect of linked sites. Here, we analyze a model for the fixation probability of deleterious mutants that hitchhike with selective sweeps in diploid, partially selfing populations. Approximate analytical solutions show that, conditional on the sweep not being lost by drift, higher inbreeding rates increase the fixation probability of the deleterious allele, due to the resulting reduction in polymorphism and effective recombination. When extending the analysis to consider a distribution of deleterious alleles, as well as the average fitness increase after a sweep, we find that beneficial alleles generally need to be more recessive than the previously assumed dominance threshold (h < 1/2) for selfing to be beneficial from one-locus theory. Our results highlight that recombination aiding the efficiency of selection on multiple loci amplifies the fitness benefits of outcrossing over selfing, compared to results obtained from one-locus theory. This effect additionally increases the parameter range under which obligate outcrossing is beneficial over partial selfing.

Marker-based investigation of inbreeding depression in the endangered species Brassica insularis

Various methods have been proposed to estimate inbreeding depression and to assess its consequences for natural populations. As an alternative to controlled crosses, the use of molecular markers has allowed direct investigation of inbreeding depression in natural populations, but usually suffers from low statistical power. Here, we investigated the effect of inbreeding depression on survival in two populations of the rare species Brassica insularis, using both controlled crosses and a marker-based approach. We compare the respective merits of the two approaches for studying inbreeding depression. We also use information from the molecular markers to dissect in detail patterns of inbreeding depression in this species. A posteriori, we find that combining the approaches was not necessary to obtain simple point estimates of inbreeding depression. However, using molecular markers may give insight into the genetic basis of inbreeding depression, such as the occurrence of epistatic interactions among deleterious alleles or purging.