Aymé Spor #

Daniel J Kvitek #

Thibault Nidelet

Juliette Martin

Judith Legrand

Christine Dillmann

Aurélie Bourgais

Dominique De Vienne

Gavin Sherlock

Delphine Sicard
email: sicard@moulon.inra.fr

Phenotypic and genotypic convergences are influenced by historical contingency and environment in yeast

Keywords: adaptive landscape, pleiotropy, Saccharomyces cerevisiae, 14-3-3 protein, BMH1 Evolution

Different organisms have independently and recurrently evolved similar phenotypic traits at different points throughout history. This phenotypic convergence may be caused by genotypic convergence and constrained by historical contingency. To investigate how convergence may be driven by selection in a particular environment and constrained by history, we analyzed nine lifehistory traits and four metabolic traits during an experimental evolution of six yeast strains in four different environments. In each of the environments, the population converged towards a different life-history strategy. However, phenotypic convergence was partly associated with the selection of mutations in genes involved in the same pathway. In a fifth of our evolution experiments, mutations in the same gene, BMH1, were selected, in three out of the six ancestral genotypes. Two types of BMH1 mutation with opposite phenotypic effects on several traits were found. The evolution of most traits, as well as the occurrence of BMH1 mutations, was significantly influenced by the ancestral strain. However, this effect could not be easily predicted from ancestors' phylogeny or past-selection. All together, our data demonstrate that phenotypic and its underlying genotypic convergence depends on a complex interplay between the evolutionary environment, pleiotropy and the ancestor genetic background but are not straightforwardly predicable.

Introduction

Throughout the tree of life, evolutionarily divergent lineages have recurrently and independently evolved similar phenotypic traits. Such convergent evolution has been observed in numerous plant and animal species [START_REF] Arendt | Convergence and parallelism reconsidered: what have we learned about the genetics of adaptation?[END_REF]Manceau et al. 2010), but evidence for trait convergence is sparse in microorganisms, likely due to the paucity of visually observable phenotypes. Examples of convergent evolution in microbes include the independent evolution of fruiting bodies and thallus type in microscopic fungi [START_REF] Rivas Plata | Parallel evolution and phenotypic divergence in lichenized fungi: A case study in the lichen-forming fungal family Graphidaceae (Ascomycota: Lecanoromycetes: Ostropales)[END_REF], and fermentation in yeast species [START_REF] Ostrowski | The genetic basis of parallel and divergent phenotypic responses in evolving populations of Escherichia coli[END_REF].

Phenotypic convergence can be caused either by natural selection or by chance. Physical or biological constraints can restrict available phenotypes to a subset of phenotypic space [START_REF] Depristo | Missense meanderings in sequence space: a biophysical view of protein evolution[END_REF][START_REF] Weinreich | Darwinian Evolution Can Follow Only Very Few Mutational Paths to Fitter Proteins[END_REF][START_REF] Gompel | The causes of repeated genetic evolution[END_REF][START_REF] Chevin | RESTRICTED PLEIOTROPY, HETEROGENOUS MUTATION, AND PARALLEL EVOLUTION[END_REF]Feldman et al. 2012), and in this restricted space, chance alone can cause phenotypic convergence. However, if convergent evolution occurred only by chance, it would not be driven by the ecological niche. By contrast, if natural selection drives convergent evolution, convergence would be expected in environments that are similar in some aspects. However, the extent to which the environment drives convergent evolution is unknown, as most studies have focused on either natural populations in uncontrolled environments [START_REF] Arendt | Convergence and parallelism reconsidered: what have we learned about the genetics of adaptation?[END_REF][START_REF] Gompel | The causes of repeated genetic evolution[END_REF][START_REF] Elmer | Adaptation in the age of ecological genomics: insights from parallelism and convergence[END_REF] or laboratoryevolved populations in a single well-controlled environment [START_REF] Wichman | Different Trajectories of Parallel Evolution During Viral Adaptation[END_REF]). In addition, most studies have focused on a limited number of traits, and it is unclear how convergent evolution for a single phenotypic trait is constrained by changes in correlated traits. Co-variation of life-history traits (i.e. traits involved in the life-cycle of an organism) and/or morphological traits has been well described (Roff 2002), with the course of evolution shaped by natural selection [START_REF] Schluter | Adaptive radiation along genetic lines of least resistance[END_REF][START_REF] Arnold | The adaptive landscape as a conceptual bridge between micro-and macroevolution[END_REF]. However, how the convergence of one trait impacts the evolution of other traits remains to be studied empirically [START_REF] Kolbe | CONVERGENT EVOLUTION OF PHENOTYPIC INTEGRATION AND ITS ALIGNMENT WITH MORPHOLOGICAL DIVERSIFICATION IN CARIBBEAN ANOLIS ECOMORPHS[END_REF].

Phenotypic convergence can occur through mutations in different sets of genes that cause similar phenotypes in distinct lineages. Alternatively, phenotypic convergence may be caused by convergent or parallel evolution at the genotypic level [START_REF] Wichman | Different Trajectories of Parallel Evolution During Viral Adaptation[END_REF][START_REF] Arendt | Convergence and parallelism reconsidered: what have we learned about the genetics of adaptation?[END_REF][START_REF] Remold | Evolutionary Genomics of Host Adaptation in Vesicular Stomatitis Virus[END_REF][START_REF] Gompel | The causes of repeated genetic evolution[END_REF][START_REF] Elmer | Adaptation in the age of ecological genomics: insights from parallelism and convergence[END_REF]Feldman et al. 2012;[START_REF] Tenaillon | The Molecular Diversity of Adaptive Convergence[END_REF]. Parallel genotypic evolution arises when mutations occur in independent lineages that start from the same genotype, while convergence refers to mutations produced from different ancestral genotypes [START_REF] Zhang | Detection of convergent and parallel evolution at the amino acid sequence level[END_REF]. Parallel and convergent genotypic evolution may arise at several levels: the same nucleotide mutating independently several times [START_REF] Wichman | Different Trajectories of Parallel Evolution During Viral Adaptation[END_REF]Rozpedowska et al. 2011), different mutations in the same gene (Rosenblum et al. 2010) or in a multigene family [START_REF] Christin | C4 Photosynthesis Evolved in Grasses via Parallel Adaptive Genetic Changes[END_REF][START_REF] Srithayakumar | Evidence for evolutionary convergence at MHC in two broadly distributed mesocarnivores[END_REF], through mutations in different genes sharing the same function (Elias and Tawfik 2012) or in the same network [START_REF] Lozovsky | Stepwise acquisition of pyrimethamine resistance in the malaria parasite[END_REF].

Phenotypic and/or genotypic evolution can be constrained by historical factors, which can produce different phenotypic and/or genotypic outcomes despite similar environmental conditions. This has been defined as historical contingency (Travisano et al 1995[START_REF] Blount | Historical contingency and the evolution of a key innovation in an experimental population of Escherichia coli[END_REF]. Historical contingency arises because the effects of mutations are contingent on the alleles that have been retained from history through epistasis. This constrains mutation roads and as a consequence the evolution of phenotypes. Selection in similar environment is supposed to eliminate the effect of historical contingency. Indeed, selection is supposed to lead to the same phenotypic solution regardless of the genotypic background, as illustrated by the numerous examples of convergence. This has lead to the idea that the effect of historical contingency should be detected at the genomic level, but may be less frequently detected at the phenotypic level, especially for traits correlated to fitness (Teotonio et al. 2009, Joshi et al. 2033, Nguyen et al. 2011[START_REF] Bedhomme | Genotypic but not phenotypic historical contingency revealed by viral experimental evolution[END_REF]. However, an extensive study including multiple traits in multiple environments is lacking to test for this hypothesis.

In this study, we analyzed the evolution of 13 metabolic and life-history traits across multiple environments and genotypes using the budding yeast Saccharomyces cerevisiae as a model system. Yeast is one of the few microorganisms for which life-history traits have been studied [START_REF] Spor | Ant" and "Grasshopper" Life-History Strategies in Saccharomyces cerevisiae[END_REF][START_REF] Spor | Niche-driven evolution of metabolic and life-history strategies in natural and domesticated populations of Saccharomyces cerevisiae[END_REF][START_REF] Granek | Pleiotropic signaling pathways orchestrate yeast development[END_REF][START_REF] Magwene | Outcrossing, mitotic recombination, and life-history trade-offs shape genome evolution in Saccharomyces cerevisiae[END_REF]Wang et al. 2011), and yeast populations display different life-history strategies depending on their ecological niche of origin [START_REF] Spor | Niche-driven evolution of metabolic and life-history strategies in natural and domesticated populations of Saccharomyces cerevisiae[END_REF]). Using experimental evolution, we ask how selection and historical contingency interplay on phenotypic and genotypic evolution. We found evidence of genotypic convergence underlying multiple trait convergence in specific environments, suggesting that mutational paths on the adaptive landscape are restricted by selection. We also found that the evolution of most traits, including fitness components, is constrained by history. By further investigating the convergence in one gene, we demonstrate that genotypic convergence underlying multi-trait convergence depends partly on the environment and on the ancestors' genetic background, highlighting the role of pleiotropy and history in shaping rugged fitness landscape.

Results

We chose six diverse S. cerevisiae strains for which there exists whole-genome sequence, with the goal of maximizing habitat, genetic and phenotypic variability [START_REF] Spor | Ant" and "Grasshopper" Life-History Strategies in Saccharomyces cerevisiae[END_REF][START_REF] Liti | Population genomics of domestic and wild yeasts[END_REF][START_REF] Spor | Niche-driven evolution of metabolic and life-history strategies in natural and domesticated populations of Saccharomyces cerevisiae[END_REF]) (Table 1). Three replicates of each strain were propagated in serial batch cultures under four selection regimes (1%_48h, 1%_96h, 15%_48h and 15%_96h) differing by the glucose content (1% or 15%) and by the cycle length (48 hours or 96 hours). The 48 h and 96 h cycle length experiments were stopped after at least 325 or 165 generations, respectively. At the end of the experiment, we studied nine life-history traits (growth rates during fermentation (R ferm ) and respiration (R resp ), the time to diauxic shift (T shift ), population sizes at T shift and the serial transfer time points (K ferm, PopSize 48h , PopSize 96h ) and cell sizes at T shift and serial transfer time points (S ferm, S 48h and S 96h )), and four metabolic traits (the specific glucose consumption rate (J spec ), the biomass yield from fermentation (Y ferm ) and the quantity of ethanol produced and released in the medium at the two serial transfer time points (Eth 48h and Eth 96h )).

Phenotypic responses to selection

We analyzed the evolution of each trait relative to its ancestral population using mixedeffects analysis of variance (see Material and Methods). The selection regime effect, which measures the differences in phenotypic evolution between the four selection regimes, across all strains, was significant for all traits except R resp , indicating that the strains evolved towards different phenotypes depending on the selection regime (Table 2). We summarize the response to selection in each selection regime using the average trait evolution across all strains (Table 3). A value of 0 indicates no detectable response to the selection on average, while a positive (negative) value indicates an increase (decrease) of the average value of the trait in the evolved populations compared to the ancestral strains. The strongest response to selection occurred in the selection regime with the strongest glucose starvation, i.e. the 1%_96h selection regime. The 1%_96h regime selected for a more extreme strategy (decreased cell size S and increased population size, PopSize), which tends to a previously described life-history strategy found among strains originating from soil and oak bark [14]. Populations selected under this regime also displayed an increase of R ferm and J spec . The 1%_48h regime selected for an intermediate life-history strategy that is close to the average life-history strategy of the ancestral populations at 48h (PopSize 48h and S 48h are near zero), however the characteristics of the evolved population are different regarding the traits in fermentation (increased R ferm , K ferm , J spec , and decreased S ferm, Table 2 and3). Globally, both 1% regimes selected for a decreased cell size compared to the 15% regimes. The difference in cycle lengths (48h versus 96h) mainly led to the evolution of metabolic changes -both 96h cycle length regimes selected for populations that left more residual ethanol compared to the 48h cycle length regimes. In 15% glucose regimes, the 96h regime also selected for an increase in the amount of time spent in fermentation (T shift ) with a decrease in the glucose consumption rate (J spec ).

Multiple traits convergence within selection regimes

To determine which traits converged most within a selection regime and diverged most between selection regimes, we carried out one linear discriminant analysis (LDA) on all traits for each evaluation medium (1 % and 15 % glucose, respectively Figure 1A and1B) setting the four selection regimes as four a priori categories. The LDA axes represent the combination of life-history trait values towards which directional selection has driven the populations in each selection regime, i.e. the multivariate phenotypic convergence landscape. The first two axes explained 89.1 % and 95.9 % of the variance in 1% glucose and 15% glucose evaluation media respectively. The probability of a population being correctly assigned to its selection regime was 0.99 in the 15% glucose evaluation medium, indicating convergence within selection regime and diversifying selection among selection regimes for several traits (see Figure 1B). In the 1% glucose evaluation medium, this probability was 0.88 in the 1% glucose evaluation medium as we were not able to discriminate between the 1%_48h regime and the 15%_48h regime (see Figure 1A). Among the 13 traits, cell size (S 96h , S 48h , S ferm ), traits related to fitness during the respiration phase (PopSize 96h , PopSize 48h , R resp ) as well as metabolic traits (J spec , Eth 48h , Eth 96h , T shift ) were significantly correlated to at least one of the LDA axes indicating that they were the most convergent traits within selection regime and divergent traits between selection regimes. By contrast, the traits related to fitness during fermentation (R ferm , K ferm , Y ferm ) did not significantly explain any of the axes of the phenotypic convergence landscape. Indeed, these traits evolved towards the same points in the different selection regimes.

Effect of historical contingency on phenotypic convergence

The effect of historical contingency was tested by a MANOVA performed on all traits of the evolved populations evaluated either in the 1% glucose or in the 15% glucose media. The ancestor effects were significant in both evaluation media (Pillai=1,5, P<0.001 in 1% glucose media; Pillai=1.9, P<0.001 in 15% glucose media) indicating that the genetic background constrained the phenotypic landscape. An ancestor effect was found separately for each trait except for T shift and Eth 48h in the 1% medium of evaluation and for R ferm. R resp . Eth 48h in the 15% medium of evaluation (Table S1). In 14 cases out of the 26 evaluated traits (13 traits in two evaluation media), an interaction effect between selection regime and ancestor effect was found indicating that the ancestor effect may depend on the environment (Table S1).

We then asked whether ancestral strains that were phenotypically more similar evolved more similar phenotypes. We measured the distances between pairwise ancestral strains and the distances between the derived pairwise evolved populations in each selection regime on the phenotypic landscape. On the LDA axis, we did not find any significant correlation between the phenotypic distances of evolved populations and their ancestors ones. When analyzing each trait separately, only five of the 52 correlation tests (13 traits × 4 selection regimes) were significant indicating that phenotypic evolution cannot be predicted easily from the ancestors' phenotype. Significant correlations were found for S ferm in the 15%_96h selection regime (r=0.7, p<0.001) as well as K ferm (r=0.76, p<0.001), Popsize96h (r=0.62, p=0.002), Y ferm (r=0.59, p=0.003) and Eth 96h (r=0.69, p<0.001) in the 1%_96H selection regime. We did not find any significant correlations between the phenotypic distances of evolved populations and the phylogenetic relationships of their ancestral strains.

Genotypic response to selection

To determine if genotypic convergence was underlying the observed phenotypic convergence, we sequenced the genomes of the ancestral strains Y55, YPS128, UWOPS83-787.3 and YJM981, and of clones that were evolved in the 1%_48h and 15%_48h conditions, for a total of eight ancestral-evolved pairs. We identified 27 SNPs and indels across the eight evolved clones relative to their respective ancestor (Table 4).

The vast majority of SNPs were non-synonymous changes (21/27). Of the remaining mutations, three were synonymous changes in protein sequences, and three were located in intergenic regions. Of the 21 non-synonymous variants, ten were nonsense mutations and eleven were missense mutations. Unlike nonsense mutations, which typically result in loss of protein or domain function, missense mutations may cause loss of function, gain of function, or no functional change. Of the eleven missense mutations, SIFT [START_REF] Kumar | Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm[END_REF] predicted that six would disrupt protein function, four would be tolerated, while one had insufficient homology data for an accurate prediction. Assuming that protein-modifying mutations are most likely to influence phenotype, these data suggest that 59% of all found variants (16/27) most likely cause the evolved phenotypes, with at least one proteinmodifying mutation per evolved clone.

We found evidence of seven copy number variants (CNVs) in the evolved clones relative to their ancestral strain (Figure 2), with four in YPS128 evolved in 1%_48h strain alone. Three of these four increase copy number, most likely resulting from a heterozygous duplication event, while the other CNV is likely a heterozygous deletion. All four CNVs are in telomeric or subtelomeric regions. The three remaining CNVs each occur in one strain apiece. Lastly, chromosome 15 appears to be a homozygously duplicated, creating two extra copies, in Y55 evolved in 15%_48h.

Convergence on genes involved in same function or pathway-We determined Gene Ontology (GO) enrichments (Boyle et al. 2004) for the genes mutated across all strains (Table 5). GO enrichments shared by both the 1%_48h and 15%_48h conditions include general signal transduction as well as Ras protein signal transduction indicating genetic convergence on genes involved in these pathway.

The 1%_48h condition was enriched for the GO terms signal transduction, MAPKKK cascade and cellular response to stimulus (Table 5). The drivers of these enrichments were BMH1, PKH1, SSK2 and SOG2. Of these, only BMH1 obviously functions in nutrient signaling [START_REF] Burbelo | 14-3-3 proteins: hot numbers in signal transduction[END_REF][START_REF] Bertram | The 14-3-3 proteins positively regulate rapamycin-sensitive signaling[END_REF]. The others are involved in cell wall integrity and osmolarity sensing through the HOG pathway. The 15%_48h condition did not have significant GO enrichments when all mutated genes were considered, but when only genes with protein-modifying mutations were used, "mitotic cell cycle" was enriched (Table 5). This enrichment was driven by TPD3, BMH1, TVP38 and CDC25. Tpd3p is a regulatory subunit of protein phosphatase 2A (PP2A) and is required for transcription by RNA pol III [START_REF] Van Zyl | Inactivation of the protein phosphatase 2A regulatory subunit A results in morphological and transcriptional defects in Saccharomyces cerevisiae[END_REF], which when mutated results in larger cell size [START_REF] Jorgensen | Systematic Identification of Pathways That Couple Cell Growth and Division in Yeast[END_REF], as well as a general decrease in resistance to stress [START_REF] Auesukaree | Genome-wide identification of genes involved in tolerance to various environmental stresses inSaccharomyces cerevisiae[END_REF]. We observed a nonsense mutation in the clone evolved from Y55 in the 15%_48h condition, which may cause the larger cell size characteristic of the 15%_48h phenotype. Cdc25 activates Ras by phosphorylation of GDP, leading to increased signaling through the cAMP/PKA pathway, which stimulates progression through the cell cycle [START_REF] Broek | The S. cerevisiae< i> CDC25</i> gene product regulates the< i> RAS</i>/adenylate cyclase pathway[END_REF]. The CDC25 mutation in the evolved clone of YJM981 under the 15%_48h condition is a nonsense mutation, likely resulting in reduced Ras activity and decreased signaling through the cAMP/PKA pathway.

Genotypic convergence at the BMH1 locus-Because four out of the eight sequenced clones had BMH1 mutations (Table 4), we determined whether there were BMH1 mutations in the other evolved populations, and also sequenced BMH2, which encodes a paralogous protein of the 14-3-3 protein family. No mutations were detected in BMH2, but 12 out of the 60 clones had BMH1 mutations (Figure 3A), with eight being heterozygous for the mutant allele and four being homozygous. The 12 mutations were distributed across nine sites and can be grouped in two classes: six mutations are STOP mutations that occurred within the 642-690 bp positions from the ATG and six mutations are non synonymous mutations that mostly occurred at 522-534 bp from the ATG (Figure 3A). At one site, the exact same mutation occurred twice independently.

Using a mathematical model, we checked whether the number of BMH1 mutations observed at the end of the experimental evolution could be found under the hypotheses that there is no selection. By simulating 10 7 experimental evolutions of M 1 =32 populations (respectively M 2 =28 populations) with 34 bottlenecks (96h selection regime) or 69 bottlenecks (48h selection regime), we estimated that the probabilities of observing 4 mutants after evolution in the 96h selection regime or 8 mutants after evolution in the 48h selection regime under the null hypotheses that there is no selection were below 10 -7 .

BMH1 mutations were found in Y55, YJM981, UWOPS83-787.3 but not in the other ancestor strains leading to a significant ancestor effect on convergence (Figure 3A and 3B, Fisher's Exact Test, P=0.005). The occurrence of mutations in BMH1 was not related to the ancestral BMH1 allele (Table 6) and does not appear to be related to the phylogenetic distance between the ancestral strains (Figure 3B).

Phenotypic effect of BMH1 mutations-The distribution of the two classes of BMH1 mutations (truncation versus amino-acid changes) was different between the two media although not statistically significant (Fisher's Exact Test, P=0.06). STOP mutations, resulting in truncated proteins, emerged only in the 1 % glucose medium (Table 6). From the predicted 3D protein structure, these truncated proteins will lack the final helix of the Cterminus as well as a glutamine repeat known to be involved in protein/protein interactions (Figure 3C). Non synonymous mutations lead to amino-acid changes at positions 55, 101, 174, 178 of the protein (Table 6) and were found in strains evolved in both glucose concentrations. Amino-acid changes at 55, 174, and 178 are located inside the groove where Bmh1p is expected to interact with other proteins as predicted by docking with arbitrary peptides and with two known Bmh1 protein partners (Serine/threonine-protein phosphatase PP1-2 and Heat shock protein Ssb1).

We tested the effect of BMH1 mutations on the value of each of the 13 traits separately. We distinguished five classes of Bmh1 protein: the ancestral sequences of strains S288c, Y55, YJM981, the ancestral sequences of strains of YPS128, NCYC110, UWOP83-787.3, the truncated proteins by stop mutations, the proteins having amino-acid substitutions pointing toward the groove, and the protein having the unique amino-acid substitution pointing toward the back inside. The last three classes of BMH1 mutations changed the phenotypic value of a similar set of traits: Rferm, Tshift, Sferm and Rresp (Figure 4, Table 7). Interestingly, knowing the class of BMH1 mutation provides as good or better prediction of the variation of these traits than knowing the ancestral strain (ANOVA's R 2 comparison, Table 7). Globally, the mutation effects were clearer in 15 % evaluation medium than in 1 %. In 15 % glucose medium, substitutions G55D, G174D and N178S, which are located in the Bmh1p groove of strain UWOP83-787.3, decreased average R ferm (P=0.026), increased T shift (P=0.01) and R resp (P=0.047), when compared to the ancestral protein sequence W2. Interestingly, the substitution D101N located in the back of the protein in strain Y55 had also an effect on T shift and R ferm (decreased R ferm P=0.05, increased T shift, P=0.019). The effects of stop mutations were opposite to the effects of amino-acid substitutions: they increased R ferm (P=0.013 for populations evolved in 1%_48h, P=0.168 for populations evolved in1%96h), decreased T shift (P=0.02 for populations evolved in 1%_48h, P=0.056 for populations evolved in 1%_96h) and increased S ferm (P=0.086 for populations evolved in 1%_48h, P=0.047 for populations evolved in 1%_96h). Because stop mutations only occurred in the 1 % selection regimes, we also analyzed whether the effect of stop mutations changed with the media of evaluation. We found that stop mutations had antagonistic effects in 15 % and 1 % glucose media: while stop mutations increased R ferm and decreases T shift in 15 % glucose media, it decreased R ferm and increased T shift in 1 % glucose media. Note that stop mutations were selected in the media where their effect on the reproduction rate was negative, demonstrating that reproduction rate is not always the best proxy for fitness.

Discussion

By varying both the environment and genetic background during experimental evolution, we found direct evidence of multiple traits divergence between experimental conditions, as well as extensive convergent evolution of diverse genetic backgrounds to common environments at the phenotypic and genetic levels.

Divergent evolution between environments

Each environment selected for a different combination of phenotypic values, indicating that diversifying selection had occurred between the four selection regimes, and that divergence was due to various life-history and metabolic traits rather than a single fitness component. Results from numerous experiments have shown that environmental variation in space and time indeed leads to ecologically divergent populations (Rainey et al. 2000;Jessup et al. 2004;Kassen and Rainey 2004;MacLean 2005;Habets et al. 2006, Kolbe et al. 2012). In yeast, it is known that populations from natural habitats such as oak and soil have a small cell size, large population and high growth rate, while yeast populations from domesticated habitats (beer, wine, bread) have a bigger cell size, a higher survival rate but a lower growth rate and population size [START_REF] Spor | Ant" and "Grasshopper" Life-History Strategies in Saccharomyces cerevisiae[END_REF][START_REF] Orr | Fitness and its role in evolutionary genetics[END_REF][START_REF] Albertin | Population Size Drives Industrial Saccharomyces cerevisiae Alcoholic Fermentation and Is under Genetic Control[END_REF]. Here, we show in the lab that decreasing the amount of glucose resource and increasing the time of glucose starvation (1% _96h selection regime) indeed select for yeast populations of small cell size but high growth rate and high population size in comparison to environments with high sugar availability. By analyzing evolution in action under differing controlled environments, we provide direct evidence that differences in sugar availability, as found between forest and domesticated environments, can explain the difference of life-history strategies in the yeast ecological niches.

Convergent evolution in phenotype and genotype

Phenotypic convergent evolution is a well-documented phenomenon that shapes the adaptation of diverged organisms evolving in a common niche. Most studies investigating phenotypic convergence have focused on a single trait and/or a single environment. Here, we provide direct experimental support for multi-trait convergence in response to selection in four different environments. We chose ancestral strains that were widely distributed on the S. cerevisiae phylogenomic tree [START_REF] Liti | Population genomics of domestic and wild yeasts[END_REF]) that were also phenotypically divergent for their life-history and metabolic strategies [START_REF] Spor | Niche-driven evolution of metabolic and life-history strategies in natural and domesticated populations of Saccharomyces cerevisiae[END_REF]. Despite this diversity, strains evolved towards a common phenotype in each selection regime. In Escherichia coli, Fong et al demonstrated parallel genotypic evolution and convergent growth phenotypes after evolving nine replicates of a single strain in two environmental conditions [START_REF] Fong | Parallel adaptive evolution cultures of Escherichia coli lead to convergent growth phenotypes with different gene expression states[END_REF]. Our study extends these findings to show that phenotypic convergence occurred between different genetic backgrounds.

After sequencing the genomes in a subset of our evolved strains, we find evidence that the same pathways, genes and even nucleotides are recurrently mutated. The Ras/cAMP/PKA pathway and protein phosphatase 2A (PP2A) complex are repeated targets, with mutations in the genes BMH1, CDC25, CYR1 and TPD3, PPM1, respectively. The gene BMH1 is mutated in at least 20% of the populations, and four sites were mutated independently twice. These repeated mutational events strongly suggest that phenotypic convergence can be caused by convergence at the genotypic and pathway level. Mutations in the Ras/cAMP pathway had already been found in many other evolution experiments performed in chemostats [START_REF] Kao | Molecular characterization of clonal interference during adaptive evolution in asexual populations of Saccharomyces cerevisiae[END_REF][START_REF] Kvitek | Reciprocal Sign Epistasis between Frequently Experimentally Evolved Adaptive Mutations Causes a Rugged Fitness Landscape[END_REF][START_REF] Wenger | Hunger Artists: Yeast Adapted to Carbon Limitation Show Trade-Offs under Carbon Sufficiency[END_REF], indicating that this pathway is a frequent target of adaptive mutation both in continuous and batch culture.

Finding recurrent mutations in BMH1 was unexpected, and to our knowledge no previous experimental evolution studies have discovered adaptive mutations in this gene. BMH1 is highly pleiotropic, and have been shown to bind at least 271 proteins [START_REF] Kakiuchi | Proteomic Analysis of in Vivo 14-3-3 Interactions in the Yeast Saccharomyces cerevisiae †[END_REF]). Molecular analyses have implicated it in many functions such as chaperone-like protein, protein tethers, cell signaling, cell cycle control, transcription regulation, post-transcriptional regulation, life-span and apoptosis [START_REF] Aitken | 14-3-3 proteins: A historic overview[END_REF][START_REF] Van Heusden | Yeast 14-3-3 proteins[END_REF][START_REF] Bruckmann | Post-Transcriptional Control of the Saccharomyces cerevisiae Proteome by 14-3-3 Proteins[END_REF]Paul and van Heusden 2009;[START_REF] Wang | Deleting the 14-3-3 Protein Bmh1 Extends Life Span in Saccharomyces cerevisiae by Increasing Stress Response[END_REF][START_REF] Clapp | 14-3-3 Protects against stress-induced apoptosis[END_REF][START_REF] Veisova | Role of individual phosphorylation sites for the 14-3-3-protein-dependent activation of yeast neutral trehalase Nth1[END_REF]). Here, we showed that Bmh1 has an impact on the reproduction rate in fermentation and respiration, the length of fermentation (until the diauxic shift) and on cell size. Which traits controlled by Bmh1p have been selected for in our experiment remains to be studied.

We also see evidence for phenotypic convergence caused by mutations in disparate pathways in the 15%_48h condition. Mutations in genes within the PP2A complex (TPD3, PPM1), Ras pathway/network (CDC25, BMH1), glucose sensing pathway (MTH1), and membrane potential/cytoplasmic pH regulation (HRK1) are all likely causative mutations. However, each mutation affects a different pathway, with multiple pathways often being mutated in individual clones, suggesting that disparate genetic changes can converge on a common phenotype. [START_REF] Fong | Parallel adaptive evolution cultures of Escherichia coli lead to convergent growth phenotypes with different gene expression states[END_REF] demonstrated that the transcriptional states of their evolved populations were very different from each other, despite similarity in endpoint growth phenotypes. They also showed that the evolutionary response involved an initial widespread expression shift followed by a large number of compensatory gene expression changes. Recently, other studies have also highlighted that phenotypic convergence could be achieved through diverse genetic mechanisms (Steiner et al. 2009;[START_REF] Chou | Optimization of Gene Expression through Divergent Mutational Paths[END_REF].

The cause of convergence-It is a common debate whether convergent evolution occurs because of natural selection toward a common adaptive phenotype, or for reasons unrelated to adaptation such as random evolutionary changes or constraints inherent to the biological system [START_REF] Depristo | Missense meanderings in sequence space: a biophysical view of protein evolution[END_REF][START_REF] Weinreich | Darwinian Evolution Can Follow Only Very Few Mutational Paths to Fitter Proteins[END_REF][START_REF] Gompel | The causes of repeated genetic evolution[END_REF]Feldman et al. 2012). Here, we provide strong evidence of convergent phenotypic and genotypic evolution occurring because of adaptation.

We show that mutations in a single gene converge specifically in one environment and have antagonistic effects in other environments, which is strong evidence that convergence occurred due to selection. Other mutations in the same gene, even at the same site, occurred and were selected for in different genetic backgrounds and in different environments. In addition, we show that positions of the mutations in BMH1 appear to be nonrandom. None of the changes are located in the N-terminus of the protein, which is involved in dimerization. Instead, observed missense and nonsense mutations are predicted to be located in regions where protein-protein interactions occur and may change the capacity for interaction of Bmh1p with its potential partners. Finally, by simulating our experimental evolution under the neutral hypothesis, we show that such genetic convergence is unlikely without selection.

Genotypic convergence is thought to occur because of variation in adaptive potential between loci: variation in mutation rate, variation in the magnitude of mutation effect, variation in the number and type of traits controlled directly or indirectly by the locus (pleiotropy, epistasy). Theoretical studies [START_REF] Chevin | RESTRICTED PLEIOTROPY, HETEROGENOUS MUTATION, AND PARALLEL EVOLUTION[END_REF]) have shown that the probability of parallel evolution is increased when pleiotropy increased. Here we found that among 8 genomes, the only gene that has been mutated independently several times is highly pleiotropic.

The effect of historical contingency

We showed that phenotypic and genotypic evolution depended on the ancestor. This was true for all but two life-history traits (R resp and R ferm in 15% evaluation medium), for all but two metabolic traits (T shift in 1% evaluation medium and Eth 48h , Table S1) as well as for the occurrence of selected mutations in BMH1 (Table 6). The influence of the initial genotypes remains apparent even if adaptive convergence within a selection regime occurred; in fact, the phenotypic evolution was even better explained by initial genotype than by selection regime in most traits (Table S1). However, predicting the phenotypic evolution based on phenotypic or genetic distances between ancestors appear not to be feasible. Results from previous studies have led to the idea that the influence of historical contingency is higher for traits that have less impact on fitness (Travisano et al. 2005[START_REF] Joshi | The contribution of ancestry, chance, and past and ongoing selection to adaptive evolution[END_REF]. Here, we present the most complete study thus far for testing this idea by including the analysis of many different cellular levels: genomic, metabolic, morphological and reproduction traits. By contrast to previous studies, we found that the ancestors effect persists for all traits' and has not disappear for traits that are strongly selected for. As illustrated by convergence in BMH1, we show that this can be partly explained by the fact that genetic convergence that is contingent on an ancestor's genetic background occur in genes that are highly pleiotropic and that determine variation in many different types of traits (metabolic, morphological, reproduction traits). Finally, we found that the effect of historical contingency depends on the environment indicating that trait plasticity should also be considered when analyzing the effect of history. All together our data highlight the role of epistasis, pleiotropy and environment into life-history evolution.

Fitness landscape in microbes

The fitness landscape topology found here correspond to a rugged landscape with various connected peaks organized around a convergent point. Whether the top of the peaks is reached in our experiment remains to be studied. Ruggedness has been related to a high level of epistasis [START_REF] Kvitek | Reciprocal Sign Epistasis between Frequently Experimentally Evolved Adaptive Mutations Causes a Rugged Fitness Landscape[END_REF]. We show here that it is also related to constraints emerging from pleiotropy. By and large, researchers have used the noncompetitive reproduction rate (also called growth rate) as a proxy for evolutionary fitness in microbes [START_REF] Orr | Fitness and its role in evolutionary genetics[END_REF]. It remains debatable how accurate of a proxy this is, since adaptive strategies can impact traits other than reproduction rate, such as the combination of traits representing fitness, often referred to as life-history traits (Roff 2002). For example, it has been shown that a simple increase in reproduction rate does not have to reflect the advantage of diploids over haploids in experimental evolution [START_REF] Gerstein | Cryptic Fitness Advantage: Diploids Invade Haploid Populations Despite Lacking Any Apparent Advantage as Measured by Standard Fitness Assays[END_REF]. Similarly, in our case, the 15%_96h selection regime did not select for a higher reproduction rate during fermentation (R ferm ). However, it selected for a decreased population size, increased cell size and increase in specific glucose consumption rate, phenotypes that suggest the cells are not increasing reproduction rate but instead are metabolizing the abundant glucose to create biomass. Thus, phenotypic adaptation of multiple traits should be considered when investigating the fitness landscape of evolved microbial populations.

In conclusion, by analyzing the evolution in action and its genotypic cause, we found direct evidence of both convergent evolution in a particular environment and divergent selection among contrasted environments. In addition we highlight the need to analyze multiple fitness components for a better understanding of adaptation. We demonstrated examples of distinct genetic changes from disparate genetic backgrounds converging on a common phenotype, as well as these disparate genetic backgrounds finding convergent genotypic solutions. Finally we show that genotypic convergence underlying multiple traits evolution is constrained by historical contingency and depends on the environment.

Materials & Methods

Yeast Strains

Six strains were chosen from the Saccharomyces Genome Re-sequencing Project (Table 1) that are distributed on the S. cerevisiae phylogenomic tree [START_REF] Liti | Population genomics of domestic and wild yeasts[END_REF], and differ for several life-history and metabolic traits [START_REF] Spor | Ant" and "Grasshopper" Life-History Strategies in Saccharomyces cerevisiae[END_REF][START_REF] Spor | Niche-driven evolution of metabolic and life-history strategies in natural and domesticated populations of Saccharomyces cerevisiae[END_REF]. They are all homothallic autodiploids, except S288c, which is haploid. For each phenotypic evaluation, a single new colony was isolated from the -80°C stock.

Experimental Evolution

Three replicates of each of the six strains were serially propagated (2 mL into 40 mL of fresh liquid medium) under four selection regimes over 5 months. Two different glucose conditions (1 % and 15 % glucose in 3 % YNB with amino acids) and two different "winter" lengths (48 h and 96 h) were chosen. The cultures were incubated at 30 °C at 200 rpm. In 1% glucose, the yeast populations first grow fermentatively and then switch to respirative growth after the diauxic shift. Cycles of 48 h have a very short period of respiration, unlike the cycles of 96 h. In 15 % glucose, the yeast populations grow only fermentatively but are subject to osmotic stress. Cycles of 48 h are composed of a long period of growth followed by a short period of stationary phase whereas cycles of 96 h have a long stationary phase. The four selection regimes are denoted 1%_48h, 1%_96h, 15%_48h and 15%_96h. All six strains were propagated in these four conditions in triplicate. Over 5 months, the populations transferred every 48 hours were transferred 74 times (about 325 generations), while the populations transferred every 96h were transferred 37 times. At each transfer, the populations were diluted 20 fold leading to an effective population size ranging from 10 6 to 10 7 depending on the strain and the selection regime. Out of the 72 evolved populations, 12 were discarded due to cross-contamination, as detected using microsatellite analysis. The 60 remaining evolved population included a duplicate of evolution for all ancestors strains in all four selection regimes, except for strain S288c (the haploid strain) in the 1_48h and 15_48h selection regimes and for strain YJM981 in the 1_96h selection regime.

Phenotyping Ancestral Strains and Evolved Populations

Phenotyping of two clones from each of the six ancestral strains and a single end-clone from each of the 60 evolved populations was performed in both 1 % and 15 % glucose as described in [START_REF] Spor | Niche-driven evolution of metabolic and life-history strategies in natural and domesticated populations of Saccharomyces cerevisiae[END_REF], leading to a total of 144 sets of phenotyping kinetics. For each kinetic, nine life-history (reproduction rate in fermentation R ferm , time to diauxic shift T shift , reproduction rate in respiration R resp , carrying capacity K ferm and cell size S ferm at the diauxic shift time point, population sizes PopSize 48h , PopSize 96h , and cell sizes S 48h ,S 96h measured at the serial transfer time points) and four metabolic traits (specific glucose consumption rate denoted J spec (g × min -1 × cell -1 ), ethanol amount denoted Eth 48h and Eth 96h expressed in g × L -1 and yield Y ferm , measured as the ratio of the biomass to the quantity of glucose consumed) were quantified.

Statistical Analysis

Response to selection-The 144 kinetics were randomly distributed into five blocks of experiments. For each trait Z, the phenotypic evolution was calculated as follows:

where the value of the ancestral strain is averaged over the replicates.

Variation of each variable ΔZ among media and strains was first analyzed using the following mixed model of analysis of variance (ANOVA):

where ΔZ is the variable, Block is the random block effect (experimental repetition, m = 1, 2, 3, 4, 5), eval medium is the evaluation condition effect (l = 1, 2), strain is the diploid strain effect (i = 1, 2, 3, 4, 5), selection regime is the selection regime effect (j = 1, 2, 3, 4), eval medium*strain (fixed), eval medium*selection regime (fixed) and strain*selection regime (fixed) are interaction effects and ε is the residual error. The interaction effects were chosen to study the interactions between the historical contingency (ancestors' effect) and selection effects. To account for multiple testing, the significance of the different effects was assessed using the false discovery rate method (ref).

For each trait, normality and homogeneity of residual distributions were checked. Least square means and REML variance estimates were obtained using the JMP ® Software. The significance of differences between selection regime means was assessed using the Tukey HSD method.

Identifying convergence and divergence-To identify whether the selection regime and the genetic background of the ancestor had an impact on the traits we carried out a MANOVA on the 13 traits. Those analyses were conducted on the phenotypic values of each trait for each replicate of evolution from each diploid ancestor in each selection regime. However, prior to those analyses, we corrected the data for a block effect estimated by a simple regression model performed trait by trait and for each evaluation medium separately. Missing values were replaced by the mean over replicates of evolution. The PopSize 48h ,PopSize 96h , and K ferm were log10 transformed. We analyzed separately the data obtained with the two evaluation medium. We discarded the S288c strain from those analyses as it was the single haploid strain.

To study further the phenotypic landscape of our evolved strains, we performed linear discriminant analyses (LDA) of those data for each evaluation medium. The four selection regimes defined the four categories used a priori to compute the discriminant functions. The a posteriori assignment probabilities to a priori categories indicate if the objects can be properly assigned to an a priori group. After using the data of diploid evolved strains to build the axes of the LDA, the ancestral strains as well as the evolved strains (including the discared xx strain) were projected on those axes. The output of the analyses in each evaluation medium (l =1, 2) were the coordinates of each strain evolved in each selection regime (j = 1, 2, 3, 4) as well as the coordinates of each ancestor strain. The haploid S288c evolved populations were excluded from the LDA analysis but projected afterwards on the phenotypic landscape in the cases of evolved populations with no missing data.

The axes of the LDA represent a linear combination of the traits that discriminate the most selection regime. Hence, the traits that correlate with the axes are those for which convergence was observed within a selection regime and divergence was observed between selection regimes. To account for multiple testing, the significance of correlations was assessed using the false discovery rate method.

Determining the genotype of evolved clones

Whole-genome sequencing and SNP/indel detection-Genomic DNA was extracted from final clones either by spooling (Treco 1987), or using Qiagen G/100 Genomic Tips. Single-end Illumina sequencing libraries were generated using the Illumina Genomic DNA Sample Prep Kit starting with 5 μg of genomic DNA, and sequencing flow cells were prepared using the Illumina Standard Cluster generation Kit. Samples were sequenced on the Illumina Genome Analyzer II, generating 36 bp reads. The data were mapped to the S288c reference (from SGD, Dec 3, 2008) and to additional contigs generated from other S. cerevisiae strains [START_REF] Dunn | Analysis of the Saccharomyces cerevisiae pangenome reveals a pool of copy number variants distributed in diverse yeast strains from differing industrial environments[END_REF]) using BWA v0.5.8 [START_REF] Li | Genome Project Data Processing Subgroup. The Sequence Alignment/Map format and SAMtools[END_REF] with default parameters.

SNP calling was performed on uniquely mapping reads only using the Genome Analysis Toolkit v1.0.4905 (McKenna et al. 2011) on all strains simultaneously with ad hoc quality filtering (see Supp. Material and Methods). SNPs were validated using Sanger sequencing.

Indels were called using Samtools v0.1.7 [START_REF] Li | Genome Project Data Processing Subgroup. The Sequence Alignment/Map format and SAMtools[END_REF]) with default settings; the false positive rate was decreased using an in-house developed method (see Supp. Material and methods). Indels were validated using Sanger sequencing.

Copy-number variation (CNV) detection-Sequence coverage was used to identify CNVs between each evolved and ancestral strain pair as in [START_REF] Araya | Whole-genome sequencing of a laboratory-evolved yeast strain[END_REF]. Briefly, raw sequencing coverage was averaged over 25 bp segments across the genomes of each evolved clone and ancestor and log2(evolved/ancestor) ratios were calculated. These log2 ratios were then adjusted by either the genome-wide log2 mean to identify whole-chromosome CNVs, or each chromosome's log2 mean to identify intrachromosomal CNVs. Genome segments were identified using the R package DNAcopy v1.22.1 [START_REF] Venkatraman | A faster circular binary segmentation algorithm for the analysis of array CGH data[END_REF] with default parameters except: smooth.CNA(smooth.region=5); segment(min.width=5, undo.splits="sdundo", undo.SD=4). Results of DNAcopy were plotted and regions of varying copy number were identified visually. Putative CNVs occupying partial chromosomes were tested using real-time quantitative PCR of genomic DNA as described [START_REF] Hoebeeck | Real-time quantitative PCR as an alternative to Southern blot or fluorescence in situ hybridization for detection of gene copy number changes[END_REF], except that each putative CNV was normalized to a single nonvarying locus located on the same chromosome as the putative CNV. The wholechromosome 15 CNV in Y55:15%_48h was normalized to a single non-varying locus on a different chromosome.

Gene Ontology (GO) enrichment analysis-GO biological process enrichments were determined using GO::TermFinder (Boyle et al. 2004) at SGD (http:// www.yeastgenome.org) using default parameters, except dubious ORFs were omitted from the analysis, with a False Discovery Rate cutoff of 0.05. All genes containing at least one mutation within its ORF were included.

Convergence on BMH1

Testing the presence of selection on BMH1-We built a theoretical model (see Supp. Material and methods) of our experimental evolution and computed the probability of observing recurrent BMH1 mutations among M evolved clones at the end of our experimental evolution under the hypotheses that there is no selection.

Evaluating the impact of BMH1 mutations-We tested for the effect of BMH1 mutations on Z the value of a single trait in an evolved population. We used the following linear model:

where Y ijlb is Z, eval medium is the effect of the culture medium (l = 1,2), selection regime (j = 1,2,3,4) is the effect of the selection regime, Bmh1p the effect of mutations in Bmh1p (b=1,2,3,4,5), (Bmh1p*eval medium) lb , (Bmh1p*selection regime) jb , (eval medium*selection regime) jl , (Bmh1p*eval medium*selection regime] ljb interaction effects, and ε ijlb the residual error. The effects of the Bmh1p protein changes (Bmh1p effect) were tested from the ANOVA using contrast. FDR corrections were carried out using a value of 15%. For three traits (Y ferm , Eth_48h, Eth_96h), the ANOVA was carried out separately for each evaluation medium because of residuals heterocedasticity.

This ANOVA model using BMH1 as sole genetic information was compared to the following model using the genetic differences between ancestral strains:

Similar R 2 between the two models would indicate that genetic variation in the gene where genetic convergence occurred explained as much traits variation than genome differences between evolved strains.

BMH1: structure and docking predictions

We generated a structural model of yeast dimeric Bmh1p to locate the amino-acid changes and truncation on the 3D protein structure. We used its homology with the human protein Characteristics of the mutations discovered in the evolved populations

Distribution of the 13 Bmh1p mutations

There are 2 different ancestral sequences of Bmh1p (W1 and W2), one in S288C, Y55, YJM981 strains, the other in YPS128, NCYC110, UWOP83-787.3 strains.

Table 7

Testing the effect of Bmh1p mutations Two ANOVA models were carried out, the first one with ancestor strains as genetic information, the other one with Bmh1p class of sequences as genetic information. The first column presents the R 2 of the first model. The other columns show the R 2 of the second model as well as the P-values associated with the Bmh1p main and interactions effects. Significant effects after FDR correction are highlighted in grey.

Figure 1 .

 1 Figure 1. Canonical plots of the linear discriminant analysisA. Phenotypic evaluation of ancestral strains and evolved populations in the 1 % glucose medium. B. Phenotypic evaluation of ancestral strains and evolved populations in the 15 % glucose medium. The correlations between traits and how they explain the linear discriminant axis are presented on the right. The colors of the symbols correspond to the different selection regimes (empty blue: 1%_48h, filled blue: 1%_96h, empty red: 15%_48h and filled red: 15%_96h). The black symbols correspond to the ancestral strains. The shapes of the symbols correspond to the evolved and ancestral state of a given strain (stars for S288c, down triangle for Y55, up triangle for YPS128, diamond for NCYC110, square for UWOPS83-787.3, and circle for YJM981).

Figure 2 .

 2 Figure 2. Copy Number Variants (CNVs) identified by sequencing coverage of the evolved strain versus the ancestral strain A) Zoomed-in coverage plots of each sub-chromosomal CNV. Red-dotted lines delimit the regions identified as CNVs by DNAcopy. Blue line is a running median of log2 ratios. B) Whole-genome coverage plot for the 97: 15%_48h evolved strain showing duplication of chromosome 15. Alternating grey/black colors are chromosomes. Red line is a running median of log2 ratios.

Figure 3 .

 3 Figure 3. Polymorphism in the BMH1 gene and encoded protein A) Each grey horizontal line represents the allele sequence of a clone isolated from a population after evolution. The ancestral strain is shown on the left. Polymorphic sites are shown with vertical lines. Nucleotides are indicated with color (T in red, A in green, C in blue, G in black). Bars represent the ancestral allele and stars show mutations that had occurred along the evolutionary course. Corresponding changes in the protein are indicated on top of the nucleotide mutations. Crosses are stop mutations. Amino-acid substitutions are indicated by their position in the protein. Mutations that have occurred in 1% selection regimes are indicated in blue, mutations that have occurred in 15% selection regimes in red. B) Neighbor-joining tree based on the proportion of different nucleotides between genomes of the six ancestral strains. The BMH1 allele is indicated by the color of the strain's name (orange or purple). The number of mutations that has occurred in BMH1 during the course of evolution is indicated below the name of each strain. C) Model of dimeric Bmh1p (region 1 to 235). A: front view, B: side view, C: back view, top row: cartoon representation, bottom row : surface representation. The two chains composing the dimer are shown in different colors (cyan and green). The positions affected by the mutation are highlighted in colors: position 55 in yellow, position 101 in orange, position 174 in red and position 178 in purple.

Figure 4 .

 4 Figure 4. Effect of Bmh1p mutations in each selection regime where the mutations had occurred as evaluated in 15% glucose medium W indicates evolved populations that have kept their ancestral protein sequences (W1 or W2), Gr indicates evolved populations that have had a substitution pointing on the groove of the protein, STOP indicates evolved populations having a truncated protein.

  Author manuscript; available in PMC 2015 May 20. Author manuscript; available in PMC 2015 May 20.

  Author manuscript; available in PMC 2015 May 20.

  Author manuscript; available in PMC 2015 May 20.

Table 1 Saccharomyces cerevisiae strains used in this study

 1

	Strains Lab Genetic Geographical Habitat Environmental	ID ID Group Origin Source	NCYC110 247 WA West Africa Industrial Ginger Beer from Z. officinale	YJM981 304 W/E Bergamo, Italy Clinical Isolate Vagina	YPS128 104 NA Pennsylvania, USA Forest Oak Exudate	UWOPS83-787.3 270 Mosaic Bahamas Fruit Opuntia stricta	S288C 96 Mosaic California, USA Laboratory Rotting fig	Y55 97 Mosaic / WA WA France Laboratory Wine	Genetic groups were obtained from inference of S. cerevisiae population structure using the program Structure (Liti et al. 2009). WA (West Africa), W/E (Wine/European), NA (North America), Mosaic /	WA (Mosaic, mostly West African).

Table 2 Analysis of variance of the relative response to selection of each life-history and metabolic trait

 2 Source is the tested effect, df the degrees of freedom. MS the mean square and F Fisher's F. The significant p-values after FDR correction (global level 0.05) are indicated in grey.

	K

ferm PopSize 48h PopSize 96h S ferm S 48h S 96h

	P	0.42	0.26	2.10 -08	6.10 -10	3.10 -08	0.1	0.017
	F	1	1.3	13.7	20.9	13.2	2.2	2.3	0.002
	P	0.024	6.10 -17	1.10 -15	1.10 -15	4.10 -14	2.10 -07	0.002
	F	2.8	116.4	32.6	40.2	27.8	14.4	2.95	0.001
	P	0.35	6.10 -18	5.10 -06	4.10 -05	4.10 -17	0.01	10 -04
	F	1.1	128.7	9.1	8.9	37.2	4	4	0.002
	P	6.10 -07	0.95	7.10 -05	7.10 -06	8.10 -13	1.10 -04	4.10 -04
	F	9.4	0.003	7.1	10.6	24.1	7.7	3.5	0.046
	P	8. 10 -04	3.10 -05	7.10 -11	9.10 -05	6.10 -04	0.216	0.05
	F	4.7	19.9	19.2	8.1	5.5	1.5	1.9	0.027
	P	5.10 -04	7.10 -26	1.10 -04	0.017	4.10 -08	0.48	0.52
	F	5	259.6	6.6	3.6	13.1	0.84	0.94	0.1	R
	df	5	1	4	3	4	3	12
	Source	Block	eval medium	strain	selection regime	eval medium*strain	eval medium*selection regime	strain*selection regime	Residuals means square

ferm R resp T shift J spec Y ferm Eth 48h Eth 96h

	P	0.300	4.10 -08	7.10 -08	6.10 -07	0.012	1.10 -06	0.104
	F	1.24	37.27	12.53	12.98	3.45	12.14	1.62	0.05
	P	6.10 -04	2.10 -21	8.10 -23	0.004	3.10 -22	0.007	0.014
	F	4.92	177.65	60.83	4.87	58.33	4.40	2.31	1.92
	P	0.004	2.10 -24	4.10 -05	0.027	2.10 -06	0.79	0.88
	F	3.8	236	7.5	3.2	10	0.35	0.55	0.09
	F P	0.77 0.552	2.23 0.141	13.4 8.10 -08	3.8 0.014	13.3 9.10 -08	3.25 0.028	1.88°0.056	3.59
	F P F P	3.01 0.016 13.67 2.10 -09	0.75 0.39 56.28 1.10 -10	4.11 0.005 1.96 0.11	0.51 0.68 5.51 0.002	1.87 0.13 1.79 0.14	2.29 0.085 1.7 0.18	0.48 0.92 1.81°0.063	0.05 0.01
	df F P	5 15.2 3.10 -10	1 18.4 5.10 -05	4 9.6 2.10 -06	3 3.4 0.02	4 6.8 1.10 -04	3 1.28 0.28	12 1.80°0.063	0.02
	Source	Block	eval medium	ancestor	selection regime	eval mediumstrain	eval medium*selection regime	ancestor*selection regime	Residuals

Evolution. Author manuscript; available in PMC 2015 May 20.

Table 3 Main effect of the selection regime on the relative response to selection of life-history and metabolic traits

 3 Mean relative response to selection is given for each trait in each selection regime. For each trait, significance of the differences was assessed by HSD Tukey tests. Levels that are not connected by a common letter (A, B or C) are significantly different.

	K
	Selection	regime

ferm PopSize 48h PopSize 96h S ferm S 48h S 96h R ferm R resp T shift J spec Eth 48h Eth 96h Y

	ferm	0.67 A	0.713 A	0.514 A	0.554
		-0.381 C	-0.208 B	-0.239 C	-0.004 A
		0.73 A	1.72 A	1.61 A	1.61 A
		1.467 A	0.4 AB	1.61 AB	0.155 B
		-0.030 B	0.013 B	-0.028 B	0.070 A
		NS	0.113 A	0.085 A	0.069 A
		0.083 A	0.060 AB	0.021 AB	-0.03 B
		-0.061 B	-0.098 C	-0.018 A	-0.008 A
		-0.070 B	-0.115 C	-0.003 A	-0.012 A
		-0.137 BC	-0.160 C	-0.085 A	-0.116 AB
		0.569 B	0.832 A	0.499 B	0.461 B
		0.088 B	0.247 A	-0.018 B	0.013 B
		1.1 A	1.08 AB	0.87 AB	0.82 B
		1%_48h	1%_96h	15%_48h	15%_96h

A

Evolution. Author manuscript; available in PMC 2015 May 20.

Table 4

 4

Table 6

 6

Evolution. Author manuscript; available in PMC 2015 May 20.

Acknowledgements

The authors would like to thank Gianni Liti for kindly providing strains. We are also grateful to Mickael Brockhurst, the GQF team and Frank Rosenzweig for fruitful discussions and critical reading of the manuscript.

Read Archive under accession SRA029322.1.

Supplementary Material

Refer to Web version on PubMed Central for supplementary material.

Author Manuscript