Diversity matters in wheat mixtures: A genomic survey of the impact of genetic diversity on the performance of 12 way durum wheat mixtures grown in two contrasted and controlled environments.

In ecology, an increase in genetic diversity within a community in natural ecosystems increases its productivity, while in evolutionary biology, kinship selection predicts that relatedness on social traits improves fitness. Varietal mixtures, where different genotypes are grown together, show contrasting results, especially for grain yield where both positive and negative effects of mixtures have been reported. To understand the effect of diversity on field performance, we grew 96 independent mixtures each composed with 12 durum wheat (Triticum turgidum ssp. durum Thell.) inbred lines, under two contrasting environmental conditions for water availability. Using dense genotyping, we imputed allelic frequencies and a genetic diversity index on more than 96000 loci for each mixture. We then analyzed the effect of genetic diversity on agronomic performance using a genome-wide approach. We explored the stress gradient hypothesis, which proposes that the greater the unfavourable conditions, the more beneficial the effect of diversity on mixture performance. We found that diversity on average had a negative effect on yield and its components while it was beneficial on grain weight. There was little support for the stress gradient theory. We discuss how to use genomic data to improve the assembly of varietal mixtures.

Evolution of flowering time in a selfing annual plant: Roles of adaptation and genetic drift.

Resurrection studies are a useful tool to measure how phenotypic traits have changed in populations through time. If these trait modifications correlate with the environmental changes that occurred during the time period, it suggests that the phenotypic changes could be a response to selection. Selfing, through its reduction of effective size, could challenge the ability of a population to adapt to environmental changes. Here, we used a resurrection study to test for adaptation in a selfing population of Medicago truncatula, by comparing the genetic composition and flowering times across 22 generations. We found evidence for evolution toward earlier flowering times by about two days and a peculiar genetic structure, typical of highly selfing populations, where some multilocus genotypes (MLGs) are persistent through time. We used the change in frequency of the MLGs through time as a multilocus fitness measure and built a selection gradient that suggests evolution toward earlier flowering times. Yet, a simulation model revealed that the observed change in flowering time could be explained by drift alone, provided the effective size of the population is small enough (<150). These analyses suffer from the difficulty to estimate the effective size in a highly selfing population, where effective recombination is severely reduced.

Fitness consequences of hybridization in a predominantly selfing species: insights into the role of dominance and epistatic incompatibilities.

Studying the consequences of hybridization on plant performance is insightful to understand the adaptive potential of populations, notably at local scales. Due to reduced effective recombination, predominantly selfing species are organized in highly homozygous multi-locus-genotypes (or lines) that accumulate genetic differentiation both among- and within-populations. This high level of homozygosity facilitates the dissection of the genetic basis of hybrid performance in highly selfing species, which gives insights into the mechanisms of reproductive isolation between lines. Here, we explored the fitness consequences of hybridization events between natural inbred lines of the predominantly selfing species Medicago truncatula, at both within- and among-populations scales. We found that hybridization has opposite effects pending on studied fitness proxies, with dry mass showing heterosis, and seed production showing outbreeding depression. Although we found significant patterns of heterosis and outbreeding depression, they did not differ significantly for within- compared to among-population crosses. Family-based analyses allowed us to determine that hybrid differentiation was mostly due to dominance and epistasis. Dominance and/or dominant epistatic interactions increased dry mass, while decreasing seed production, and recessive epistatic interactions mostly had a positive effect on both fitness proxies. Our results illustrate how genetic incompatibilities can accumulate at a very local scale among multi-locus-genotypes, and how non-additive genetic effects contribute to heterosis and outbreeding depression.

How and When Does Outcrossing Occur in the Predominantly Selfing Species Medicago truncatula?

Empirical studies on natural populations of Medicago truncatula revealed selfing rates higher than 80%, but never up to 100%. Similarly, several studies of predominantly selfing species show variability in the level of residual outcrossing between populations and also between temporal samples of the same population. However, these studies measure global selfing rates at the scale of the population and we do not know whether there is intra-population variation and how outcrossing events are distributed, between genotypes, plants, flowers, or seeds. Theoretical studies predict the maintenance of residual outcrossing in highly selfing species due to environmental (e.g., pollen biology) and/or genetic determinants and decompositions of the variation in outcrossing rate using experimental data can be very informative to test these hypotheses. Here, we focus on one natural population of M. truncatula in order to describe precisely its mating system. In particular, we investigated the determinants of the selfing rate by testing for seasonal variations (environmental determinism) and variations between genotypes (genetic determinism). We measured selfing rates in maternal progenies from plants collected widely across a natural population. For each plant, we collected pods from flowers produced at the beginning and at the end of the flowering season to test for a seasonal variation in the outcrossing rate. For each collected offspring, we also estimated the likelihood that it was issued from a self-fertilization event and assessed the genetic component of variation of this mating system measure. We found a significant, albeit small, increase in outcrossing rate in progenies collected at the end [t m = 0.137 (SD = 0.025)] compared to those collected at the beginning [t m = 0.083 (0.016)] of the flowering season. A significant between genotypes variation in selfing rate was also detected, resulting in a heritability of 9% for the rate of residual outcrossing. Altogether, our work shows that despite a predominantly selfing reproductive mode, M. truncatula displays variation in residual outcrossing rate, and that this trait is likely under a complex determinism combining environmental and genetic factors. We discuss the evolutionary implications of our results for the population.

How does selfing affect the genetic variance of quantitative traits? An updated meta-analysis on empirical results in angiosperm species.

Most theoretical works predict that selfing should reduce the level of additive genetic variance available for quantitative traits within natural populations. Despite a growing number of quantitative genetic studies undertaken during the last two decades, this prediction is still not well supported empirically. To resolve this issue and confirm or reject theoretical predictions, we reviewed quantitative trait heritability estimates from natural plant populations with different rates of self-fertilization and carried out a meta-analysis. In accordance with models of polygenic traits under stabilizing selection, we found that the fraction of additive genetic variance is negatively correlated with the selfing rate. Although the mating system explains a moderate fraction of the variance, the mean reduction of narrow-sense heritability values between strictly allogamous and predominantly selfing populations is strong, around 60%. Because some nonadditive components of genetic variance become selectable under inbreeding, we determine whether self-fertilization affects the relative contribution of these components to genetic variance by comparing narrow-sense heritability estimates from outcrossing populations with broad-sense heritability estimated in autogamous populations. Results suggest that these nonadditive components of variance may restore some genetic variance in predominantly selfing populations; it remains, however, uncertain how these nonadditive components will contribute to adaptation.

Structure of multilocus genetic diversity in predominantly selfing populations.

Predominantly selfing populations are expected to have reduced effective population sizes due to nonrandom sampling of gametes, demographic stochasticity (bottlenecks or extinction-recolonization), and large scale hitchhiking (reduced effective recombination). Thus, they are expected to display low genetic diversity, which was confirmed by empirical studies. The structure of genetic diversity in predominantly selfing species is dramatically different from outcrossing ones, with populations often dominated by one or a few multilocus genotypes (MLGs) coexisting with several rare genotypes. Therefore, multilocus diversity indices are relevant to describe diversity in selfing populations. Here, we use simulations to provide analytical expectations for multilocus indices and examine whether selfing alone can be responsible for the high-frequency MLGs persistent through time in the absence of selection. We then examine how combining single and multilocus indices of diversity may be insightful to distinguish the effects of selfing, population size, and more complex demographic events (bottlenecks, migration, admixture, or extinction-recolonization). Finally, we examine how temporal changes in MLG frequencies can be insightful to understand the evolutionary trajectory of a given population. We show that combinations of selfing and small demographic sizes can result in high-frequency MLGs, as observed in natural populations. We also show how different demographic scenarios can be distinguished by the parallel analysis of single and multilocus indices of diversity, and we emphasize the importance of temporal data for the study of predominantly selfing populations. Finally, the comparison of our simulations with empirical data on populations of Medicago truncatula confirms the pertinence of our simulation framework.

Discontinuities in quinoa biodiversity in the dry Andes: An 18-century perspective based on allelic genotyping.

History and environment shape crop biodiversity, particularly in areas with vulnerable human communities and ecosystems. Tracing crop biodiversity over time helps understand how rural societies cope with anthropogenic or climatic changes. Exceptionally well preserved ancient DNA of quinoa (Chenopodium quinoa Willd.) from the cold and arid Andes of Argentina has allowed us to track changes and continuities in quinoa diversity over 18 centuries, by coupling genotyping of 157 ancient and modern seeds by 24 SSR markers with cluster and coalescence analyses. Cluster analyses revealed clear population patterns separating modern and ancient quinoas. Coalescence-based analyses revealed that genetic drift within a single population cannot explain genetic differentiation among ancient and modern quinoas. The hypothesis of a genetic bottleneck related to the Spanish Conquest also does not seem to apply at a local scale. Instead, the most likely scenario is the replacement of preexisting quinoa gene pools with new ones of lower genetic diversity. This process occurred at least twice in the last 18 centuries: first, between the 6th and 12th centuries-a time of agricultural intensification well before the Inka and Spanish conquests-and then between the 13th century and today-a period marked by farming marginalization in the late 19th century likely due to a severe multidecadal drought. While these processes of local gene pool replacement do not imply losses of genetic diversity at the metapopulation scale, they support the view that gene pool replacement linked to social and environmental changes can result from opposite agricultural trajectories.

Author Correction: Intermediate degrees of synergistic pleiotropy drive adaptive evolution in ecological time.

In the version of this Article previously published, there was a typographical error ('4' instead of '2') in the equations relating F ST and effective population size (N e) in the Methods section 'Genome-wide scan for selection based on temporal differentiation'. The correct equations are given below.[Formula: see text] [Formula: see text].

How to optimize the precision of allele and haplotype frequency estimates using pooled-sequencing data.

Sequencing pools of individuals rather than individuals separately reduces the costs of estimating allele frequencies at many loci in many populations. Theoretical and empirical studies show that sequencing pools comprising a limited number of individuals (typically fewer than 50) provides reliable allele frequency estimates, provided that the DNA pooling and DNA sequencing steps are carefully controlled. Unequal contributions of different individuals to the DNA pool and the mean and variance in sequencing depth both can affect the standard error of allele frequency estimates. To our knowledge, no study separately investigated the effect of these two factors on allele frequency estimates; so that there is currently no method to a priori estimate the relative importance of unequal individual DNA contributions independently of sequencing depth. We develop a new analytical model for allele frequency estimation that explicitly distinguishes these two effects. Our model shows that the DNA pooling variance in a pooled sequencing experiment depends solely on two factors: the number of individuals within the pool and the coefficient of variation of individual DNA contributions to the pool. We present a new method to experimentally estimate this coefficient of variation when planning a pooled sequencing design where samples are either pooled before or after DNA extraction. Using this analytical and experimental framework, we provide guidelines to optimize the design of pooled sequencing experiments. Finally, we sequence replicated pools of inbred lines of the plant Medicago truncatula and show that the predictions from our model generally hold true when estimating the frequency of known multilocus haplotypes using pooled sequencing.

Intermediate degrees of synergistic pleiotropy drive adaptive evolution in ecological time.

Rapid phenotypic evolution of quantitative traits can occur within years, but its underlying genetic architecture remains uncharacterized. Here we test the theoretical prediction that genes with intermediate pleiotropy drive adaptive evolution in nature. Through a resurrection experiment, we grew Arabidopsis thaliana accessions collected across an 8-year period in six micro-habitats representative of that local population. We then used genome-wide association mapping to identify the single-nucleotide polymorphisms (SNPs) associated with evolved and unevolved traits in each micro-habitat. Finally, we performed a selection scan by testing for temporal differentiation in these SNPs. Phenotypic evolution was consistent across micro-habitats, but its associated genetic bases were largely distinct. Adaptive evolutionary change was most strongly driven by a small number of quantitative trait loci (QTLs) with intermediate degrees of pleiotropy; this pleiotropy was synergistic with the per-trait effect size of the SNPs, increasing with the degree of pleiotropy. In addition, weak selection was detected for frequent micro-habitat-specific QTLs that shape single traits. In this population, A. thaliana probably responded to local warming and increased competition, in part mediated by central regulators of flowering time. This genetic architecture, which includes both synergistic pleiotropic QTLs and distinct QTLs within particular micro-habitats, enables rapid phenotypic evolution while still maintaining genetic variation in wild populations.

Adaptation to climate through flowering phenology: a case study in Medicago truncatula.

Local climatic conditions likely constitute an important selective pressure on genes underlying important fitness-related traits such as flowering time, and in many species, flowering phenology and climatic gradients strongly covary. To test whether climate shapes the genetic variation on flowering time genes and to identify candidate flowering genes involved in the adaptation to environmental heterogeneity, we used a large Medicago truncatula core collection to examine the association between nucleotide polymorphisms at 224 candidate genes and both climate variables and flowering phenotypes. Unlike genome-wide studies, candidate gene approaches are expected to enrich for the number of meaningful trait associations because they specifically target genes that are known to affect the trait of interest. We found that flowering time mediates adaptation to climatic conditions mainly by variation at genes located upstream in the flowering pathways, close to the environmental stimuli. Variables related to the annual precipitation regime reflected selective constraints on flowering time genes better than the other variables tested (temperature, altitude, latitude or longitude). By comparing phenotype and climate associations, we identified 12 flowering genes as the most promising candidates responsible for phenological adaptation to climate. Four of these genes were located in the known flowering time QTL region on chromosome 7. However, climate and flowering associations also highlighted largely distinct gene sets, suggesting different genetic architectures for adaptation to climate and flowering onset.

Spatio-temporal changes in the structure of an Australian frog hybrid zone: a 40-year perspective.

Spatio-temporal studies of hybrid zones provide an opportunity to test evolutionary hypotheses of hybrid zone maintenance and movement. We conducted a landscape genetics study on a classic hybrid zone of the south-eastern Australian frogs, Litoria ewingii and Litoria paraewingi. This hybrid zone has been comprehensively studied since the 1960s, providing the unique opportunity to directly assess changes in hybrid zone structure across time. We compared both mtDNA and male advertisement call data from two time periods (present and 1960s). Clinal analysis of the coincidence (same center) and concordance (same width) of these traits indicated that the center of the hybrid zone has shifted 1 km south over the last 40 years, although the width of the zone and the rate of introgression remained unchanged. The low frequency of hybrids, the strong concordance of clines within a time period, and the small but significant movement across the study period despite significant anthropogenic changes through the region, suggest the hybrid zone is a tension zone located within a low-density trough. Hybrid zone movement has not been considered common in the past but our findings highlight that it should be considered a crucial component to our understanding of evolution.

Pedigree-free estimates of heritability in the wild: promising prospects for selfing populations.

Estimating the genetic variance available for traits informs us about a population's ability to evolve in response to novel selective challenges. In selfing species, theory predicts a loss of genetic diversity that could lead to an evolutionary dead-end, but empirical support remains scarce. Genetic variability in a trait is estimated by correlating the phenotypic resemblance with the proportion of the genome that two relatives share identical by descent ('realized relatedness'). The latter is traditionally predicted from pedigrees (Φ A : expected value) but can also be estimated using molecular markers (average number of alleles shared). Nevertheless, evolutionary biologists, unlike animal breeders, remain cautious about using marker-based relatedness coefficients to study complex phenotypic traits in populations. In this paper, we review published results comparing five different pedigree-free methods and use simulations to test individual-based models (hereafter called animal models) using marker-based relatedness coefficients, with a special focus on the influence of mating systems. Our literature review confirms that Ritland's regression method is unreliable, but suggests that animal models with marker-based estimates of relatedness and genomic selection are promising and that more testing is required. Our simulations show that using molecular markers instead of pedigrees in animal models seriously worsens the estimation of heritability in outcrossing populations, unless a very large number of loci is available. In selfing populations the results are less biased. More generally, populations with high identity disequilibrium (consanguineous or bottlenecked populations) could be propitious for using marker-based animal models, but are also more likely to deviate from the standard assumptions of quantitative genetics models (non-additive variance).

Detecting selection along environmental gradients: analysis of eight methods and their effectiveness for outbreeding and selfing populations.

Thanks to genome-scale diversity data, present-day studies can provide a detailed view of how natural and cultivated species adapt to their environment and particularly to environmental gradients. However, due to their sensitivity, up-to-date studies might be more sensitive to undocumented demographic effects such as the pattern of migration and the reproduction regime. In this study, we provide guidelines for the use of popular or recently developed statistical methods to detect footprints of selection. We simulated 100 populations along a selective gradient and explored different migration models, sampling schemes and rates of self-fertilization. We investigated the power and robustness of eight methods to detect loci potentially under selection: three designed to detect genotype-environment correlations and five designed to detect adaptive differentiation (based on F(ST) or similar measures). We show that genotype-environment correlation methods have substantially more power to detect selection than differentiation-based methods but that they generally suffer from high rates of false positives. This effect is exacerbated whenever allele frequencies are correlated, either between populations or within populations. Our results suggest that, when the underlying genetic structure of the data is unknown, a number of robust methods are preferable. Moreover, in the simulated scenario we used, sampling many populations led to better results than sampling many individuals per population. Finally, care should be taken when using methods to identify genotype-environment correlations without correcting for allele frequency autocorrelation because of the risk of spurious signals due to allele frequency correlations between populations.

The evolution of harm--effect of sexual conflicts and population size.

Conflicts of interest between mates can promote the evolution of male traits that reduce female fitness and that drive coevolution between the sexes. The rate of adaptation depends on the intensity of selection and its efficiency, which depends on drift and genetic variability. This leads to the largely untested prediction that coevolutionary adaptations such as those driven by sexual conflict should evolve faster in large populations. We tested this using the bruchid beetle Callosobruchus maculatus, a species where harm inflicted by males is well documented. Although most experimental evolution studies remove sexual conflict, we reintroduced it in populations in which it had been experimentally removed. Both population size and standing genetic variability were manipulated in a factorial experimental design. After 90 generations of relaxed conflict (monogamy), the reintroduction of sexual conflicts for 30 generations favored males that harmed females and females that were more resistant to the genital damage inflicted by males. Males evolved to become more harmful when population size was large rather than when initial genetic variation was enriched. Our study shows that sexual selection can create conditions in which males can benefit from harming females and that selection may tend to be more intense and effective in larger populations.

Genetic architecture of metabolic rate: environment specific epistasis between mitochondrial and nuclear genes in an insect.

The extent to which mitochondrial DNA (mtDNA) variation is involved in adaptive evolutionary change is currently being reevaluated. In particular, emerging evidence suggests that mtDNA genes coevolve with the nuclear genes with which they interact to form the energy producing enzyme complexes in the mitochondria. This suggests that intergenomic epistasis between mitochondrial and nuclear genes may affect whole-organism metabolic phenotypes. Here, we use crossed combinations of mitochondrial and nuclear lineages of the seed beetle Callosobruchus maculatus and assay metabolic rate under two different temperature regimes. Metabolic rate was affected by an interaction between the mitochondrial and nuclear lineages and the temperature regime. Sequence data suggests that mitochondrial genetic variation has a role in determining the outcome of this interaction. Our genetic dissection of metabolic rate reveals a high level of complexity, encompassing genetic interactions over two genomes, and genotype × genotype × environment interactions. The evolutionary implications of these results are twofold. First, because metabolic rate is at the root of life histories, our results provide insights into the complexity of life-history evolution in general, and thermal adaptation in particular. Second, our results suggest a mechanism that could contribute to the maintenance of nonneutral mtDNA polymorphism.

Genetic, morphological, and feather isotope variation of migratory willow warblers show gradual divergence in a ring.

The circular distribution of the willow warbler Phylloscopus trochilus around the Baltic Sea shares many features with the classic examples of ring species; however, the system is much younger. It has previously been shown that a secondary contact zone is located in central Scandinavia, where there are narrow clines for several morphological traits coincident with a migratory divide. Here we analyse multiple traits and genes from > 1700 males captured on breeding territories at 77 sites spread around the Baltic Sea to test the following hypothesis. If the secondary contact zone in Scandinavia is a result of divergence in two allopatric refuge populations during the last glaciation, we expect to find a similar secondary contact zone somewhere else around the circular distribution. Our results show that the trait clines were wider and displaced from each other along the eastern side of the Baltic Sea. Analyses of 12 microsatellite loci confirmed that the genome is very similar between the terminal forms (F(ST) = 0). Two AFLP-derived markers filtered out from a genomic scan instead appear to be maintained by selection. These markers exhibited steep clines at the secondary contact zone in Scandinavia, but as for the phenotypic traits, had vastly different cline centres east of the Baltic Sea. The trait clines along the ring distribution outside the Scandinavian secondary contact zone thus seem to have been shaped by independent action of selection or drift during the process of postglacial colonization.

Comparing clines on molecular and phenotypic traits in hybrid zones: a window on tension zone models.

The study of zones of secondary contact provides insight into the maintenance of reproductive isolation. Tension zone theory supplies powerful tools for assessing how dispersal and selection shape hybrid zones. We present a multimodal analysis of phenotypic clines in conjunction with clines at molecular markers in a hybrid zone between Larus glaucescens and Larus occidentalis. We developed a new method to analyze simultaneously clines of quantitative traits and molecular data. Low linkage disequilibrium and the lack of coincidence between clines at six microsatellites, a mitochondrial DNA region, and two phenotypic traits indicated introgression. However, the hypothesis of neutral diffusion was rejected based on evidence that all of the clines were concordant and narrower than expected for neutral clines, indicating some indirect selection. The analysis of phenotypic variance gave evidence of restricted phenotypic introgression and together with the bimodal distribution of phenotypes suggested that disruptive selection is acting across the hybrid zone, especially on the coloration of bare parts. Multimodal analysis of phenotypic clines also highlighted a shift between the peak of intermediates and the cline center, left behind by hybrid zone motion. High-resolution analysis of phenotypes distribution thus proved useful for detecting hybrid zone movement even without temporal data.