Hybridization and adaptive evolution of diverse Saccharomyces species for cellulosic biofuel production

Synthetic hybrids have intermediate traits or display heterosis

To determine whether interspecies hybridization could introduce traits from other species relevant to biofuel production, we generated eight synthetic hybrids. Specifically, we crossed haploid strains of S. cerevisiae previously engineered and experimentally evolved for xylose fermentation (Y128 and Y101, a heterothallic haploid derivative of Y73) [24, 25] with previously generated haploid derivatives of the reference strains of S. mikatae, S. kudriavzevii, S. uvarum [28, 50], and our new heterothallic haploid derivative of S. eubayanus (Table 1; Fig. 1), four of the most divergent members of the genus.

Strain	Synonym	Parent	Species	Population	Genotype	Comment	Publication
Wild strains
yHDPN14	FM1117		Saccharomyces paradoxus	A	MAT a/MAT α?; ρ +  (C2_H69)		This study
yHAB336	SDPGP3		Saccharomyces mikatae		MAT a/MAT α?; ρ +  (C2_H162)		This study
yHAB407	XXYS16L-5		Saccharomyces kudriavzevii		MAT a/MAT α?; ρ +  (C2_H146)		This study
yHAB413	HZZt19L.1		Saccharomyces arboricola		MAT a/MAT α?; ρ +  (C2_H158)		This study
yHCT77	FM1277		Saccharomyces uvarum	HOL-NA	MAT a/MAT α?; ρ +  (C2_H77)		[60]
yHCT69	yHDPN9		Saccharomyces eubayanus	PB-HOL	MAT a/MAT α?; ρ +  (C2_H97)		This study
Engineered strains
GLBRCY73	yHDPN35		Saccharomyces cerevisiae		MAT a/MAT α; hoΔ::XYL123-KanMX; ρ +  (C2_H13)	Evolved in YPX, aerobically	[24]
GLBRCY101	yHDPN47		Saccharomyces cerevisiae		MAT a hoΔ::XYL123-KanMX; ρ +  (C2_H13)	Dissected from GLBRCY73	This study
GLBRCY128	yHDPN36		Saccharomyces cerevisiae		hoΔ::ScTAL1-CpXylA-SchXYL3-KanMX; MAT a; ρ +  (C2_H13)	Evolved in YPX, anaerobically	[25]
yHSSS217		IFO1815T	Saccharomyces mikatae	Asia-A	MAT α hoΔ::NatMX; ρ +  (C2_H74)		[28]
FM1097		IFO1802T	Saccharomyces kudriavzevii	Asia-A	MAT α hoΔ::NatMX; ρ +  (C2_H83)		[28]
yHSSS101		CBS7001	Saccharomyces uvarum	HOL-EU	MAT α hoΔ::NatMX; ρ +  (C2_H77)		[28]
yHEB69		FM1318	Saccharomyces eubayanus	PB-HOL	MAT α hoΔ::TKMX; ρ +  (C2_H103)		This study
Synthetic hybrids
yHDPN1	yHWA103		Y101 × S. mikatae yHSSS217		MAT a/MAT α; hoΔ::XYL123-KanMX/hoΔ::NatMX; ρ +  (Recombinant)		This study
yHDPN5	yHWA93		Y101 × S. kudriavzevii FM1097		MAT a/MAT α; hoΔ::XYL123-KanMX/hoΔ::NatMX; ρ +  (S. cerevisiae)		This study
yHDPN379			Y101 × S. kudriavzevii FM1097		MAT a/MAT α; hoΔ::XYL123-KanMX/hoΔ::NatMX; ρ +  (S. cerevisiae)	yHDPN5 evolved in ACSH, microaerobically, 50 G	This study
yHDPN399			Y101 × S. mikatae yHSSS217		MAT a/MAT α; hoΔ::XYL123-KanMX/hoΔ::NatMX; ρ +  (Recombinant)	yHDPN1 evolved in ACSH, microaerobically, 50 G	This study

Geographical, ecological, and genotypic information of selected Saccharomyces species for generating hybrids, or for conducting the three replicate fermentations in four media conditions at 24 °C, for wild isolates, or 30 °C for synthetic hybrids

C2_HXXX: Indicates the COX2 haplotype based on [100]

Sch, Scheffersomyces (Pichia) stipitis; Cp, Clostridium phytofermentans; Sc, Saccharomyces cerevisiae; G, generations

XYL123, abbreviation for the SchXYL1-SchXYL2-SchXYL3 cassette

In microaerobic conditions (microtiter plates), Y128 grew more slowly than Y73 (t test, p value 2.379 × 10⁻⁰⁶) (Fig. 2; Additional file 1), likely due to the disruption of key signaling pathways during the adaptation of Y128 to anaerobic growth in YPX media [25, 26] (Fig. 2). In general, hybrids isolated from fermented beverages grew more slowly than Y73 (t test, p value 1.551e−06) and the synthetic hybrids made with Y73 (t test, p value = 9.741e−05) (Fig. 2a). Thus, hybrid status itself does not enable better growth in ACSH; instead, the specific parents of the hybrids likely confer the relevant properties.

During growth in ACSH microtiter plates, the synthetic hybrids generally had intermediate kinetic traits relative to both parents (Fig. 2; Additional file 1). For example, Y73 × S. mikatae, Y128 × S. mikatae, and Y128 × S. uvarum grew faster than their S. cerevisiae parents (t test, p value <0.0010), while Y73 × S. kudriavzevii and Y73 × S. eubayanus grew faster than their non-S. cerevisiae parents (t test, p value <0.0300). One hybrid, Y128 × S. uvarum, showed significant heterosis by growing faster than both parents (t test, p value <0.0057). These data suggest that ACSH kinetic traits are heritable and, in many cases, at least partially dominant.

The synthetic hybrids retained the critical xylose consumption trait of their engineered S. cerevisiae parents, except for Y73 × S. kudriavzevii, whose extracellular xylose concentration remained unchanged after 14 days in culture at 24 °C (Additional file 2). Although the concentration of extracellular xylose decreased in most of the synthetic hybrids, the final concentration was higher than the S. cerevisiae parents (t test, p value = 0.0309). These results suggest that interspecies hybridization can generate phenotypic diversity, while often retaining important parental kinetic and fermentative traits.

Adaptive evolution improved the fermentative traits of synthetic hybrids

To determine whether synthetic hybrids could be improved through adaptive evolution, six synthetic hybrids (those crossing Y73 or Y128 to either S. mikatae, S. kudriavzevii, or S. uvarum) were evolved for 50 generations in triplicate (Fig. 1). Each round consisted of fermentations in ACSH at 30 °C for 14 days under microaerobic conditions. Since glucose was exhausted by the second day, this stringent condition forced the cells toward starvation if they did not consume xylose, while simultaneously forcing them to retain tolerance to HTs and ethanol. Most of the synthetic hybrids did not survive these harsh conditions, and only two replicates of Y73 × S. mikatae and one replicate of Y73 × S. kudriavzevii reached the complete 50 generations of adaptive evolution. At the end of the adaptive evolution experiment, 10 colonies of each surviving hybrid were grown in Minimal Medium (MM) containing 2% xylose, and the fastest-growing variant was selected for further analysis.

Next, we tested the growth rate of the ancestral synthetic hybrids and evolved hybrids in microtiter cultures we also examined metabolites of the glycolytic and pentose phosphate pathways (Fig. 3b; Additional file 4) during fermentations of the evolved hybrids, focusing on ACSH, YPDX (a rich medium containing glucose and xylose that matches their concentrations in ACSH), YPDX plus the HTs found in ACSH, and YPDX plus feruloyl amide (FA), a product of AFEX pretreatment previously shown to have inhibitory effects in E. coli [9, 12]. Despite lack of growth improvement in ACSH (Additional file 3), extracellular xylose concentrations in microaerobic ACSH fermentations decreased significantly faster in cultures of evolved hybrids than ancestral hybrids (t test, p value <0.0368) (Figs. 3, 4). The consumption of xylose translated into a slight but significant increase of ethanol during ACSH fermentation in the evolved hybrid Y73 × S. mikatae, as well as significant increases in other conditions for both evolved hybrids (Additional file 5). These results suggest that the evolved hybrids adapted to improve xylose fermentation at 30 °C in the presence of HTs, although respiration in the microaerobic conditions also likely occurred at later time points (Fig. 4). Other relevant traits included a significant reduction in extracellular concentrations of acetate and glycerol (t test, p value <0.0121) (Fig. 3).

To test the performance of our synthetic hybrids in more industrially relevant conditions, we also performed ACSH fermentations at 30 °C in anaerobic conditions for 7 days (Fig. 3; Additional files 6A, F, 7A, F). Xylose consumption and ethanol production in ACSH under anaerobic conditions remained qualitatively similar compared to microaerobic conditions. Interestingly, the extracellular xylose concentrations decreased faster in both evolved hybrids relative to ancestral hybrids (t test, p value <0.04886).

Adaptive evolution drives genome rearrangements

Previous studies have shown that the genomes of interspecies hybrids are unstable [51–53]. To test whether our evolved hybrids experienced genome instability during their adaptation to ACSH, we sequenced their genomes, measured their relative DNA content by flow cytometry, and compared these data with that from their ancestors. As expected, ancestral hybrids were diploid, and the evolved Y73 × S. mikatae hybrid remained approximately diploid (Additional file 8). In contrast, flow cytometry and read coverage data indicated that the evolved hybrid Y73 × S. kudriavzevii was approximately triploid, suggesting an increase of the DNA content per cell (Additional file 8).

Genome-wide ploidy levels were estimated using a reference-based coverage count of sequencing reads (Fig. 5). Chromosome copy numbers for both subgenomes in the ancestral hybrids followed the expected 1:1 ratio (Fig. 5a, c). However, both evolved hybrids contained gross chromosomal rearrangements, as suggested by increased read coverage for one of the parents and decreased read coverage for the other parent. The evolved hybrid of Y73 × S. mikatae contained interspecies translocations and copy number changes involving at least five chromosomes: IV, VII, VIII, XI, and XII. For example, the amplified left arm of Y73 chromosome IV in the evolved hybrid Y73 × S. mikatae contained the Sch. stipitis xylose utilization genes engineered into Y73 to confer xylose utilization, suggesting that these genes were amplified to increase their expression, while the homologous region of the S. mikatae genome was lost (Fig. 5b). The evolved hybrid of Y73 × S. kudriavzevii was mostly triploid, but several chromosomes were tetrasomic (2:2), including chromosomes II, III, V, VI, IX, XI, and XV. In both evolved hybrids, several non-S. cerevisiae genomic regions were replaced by their S. cerevisiae homologous counterparts (Fig. 5b, d). A translocation in chromosome XII was shared among both evolved hybrids, occurring ~50 Kb to the left of the cluster of ribosomal DNA (rDNA). In both cases, the cluster of rDNA was inherited from S. cerevisiae subgenome. Thus, our stringent adaptive evolution experiment generated considerable genomic diversity, some of which might be involved in the fermentation improvements.

Heritable hydrolysate tolerance of S. mikatae and other Saccharomyces strains

Having shown that interspecies hybridization and adaptive evolution can generate diverse strains relevant to biofuel production, we interrogated whether genetically diverse strains of the various species of Saccharomyces might harbor additional potential (Fig. 1) for future research. We screened 507 wild strains of Saccharomyces from species other than S. cerevisiae and found a wide range of characteristics (Fig. 6A; Additional file 9). Notably, the kinetic parameters of ACSH microtiter cultures of the non-S. cerevisiae strains chosen as parents to construct the synthetic hybrids above were not necessarily the best strains for their respective species. These results suggest that additional hybridization experiments with diverse strains of Saccharomyces might unlock additional biofuel potential.

The Saccharomyces species differed significantly in their responses to ACSH (Kruskal–Wallis test p value <2.2 × 10⁻¹⁶) (Fig. 6A). A principal component analysis (PCA) demonstrated a strong correlation between maximum OD and growth rate among strains (Fig. 6C, Spearman’s correlation ρ = 0.67, p value = 2.2 × 10⁻¹⁶), while the length of the lag phase constituted a second principal component. S. mikatae and S. paradoxus produced more biomass and had faster maximum growth rates in ACSH than the other species (Fig. 6A, B). Although none of the species as a whole grew significantly faster than Y73 (Dunn’s test, corrected p value >0.05), whose genetic background was chosen specifically for ACSH tolerance [22], both species included several strains that grew faster (Fig. 6A). For example, S. mikatae IFO1815 grew faster than Y73 (t test, p value <9.7 × 10⁻⁴, Fig. 2a), a capability shared with many other S. mikatae strains. Indeed, S. mikatae grew significantly faster than the other Saccharomyces species (Dunn’s test, corrected p value <0.0455) (Fig. 6A), which suggests an innate tolerance to ACSH.

Several of our Saccharomyces strains have been classified into different populations or phylogenetic lineages (Additional file 9) [24, 25, 28, 32, 33, 43, 44, 54–63]. When we examined the kinetic parameters by splitting the strains by their lineage designation (Additional file 10), we found significant differences in growth rates between lineages for S. paradoxus and S. eubayanus (Kruskal–Wallis test, p value <0.05). For example, S. paradoxus strains from the America C (Québec) lineage [62] grew significantly faster than other lineages of this species (Dunn’s test corrected p value <0.0354). Similarly, strains from the Patagonia A and B lineages grew significantly faster than a set of closely related admixture strains (Dunn’s test corrected p value <0.0157), suggesting that this specific combination of alleles from these two lineages [33, 44] is detrimental in this particular condition. These associations further suggest that ACSH kinetic traits are heritable.

Yeast diversity translates into fermentation trait diversity

To test the growth kinetics (microtiter cultures) in several inhibitory conditions, we selected the fastest- and slowest-growing strains of each species from the above screen (Additional file 9), as well as the fastest from each lineage (Fig. 1). Specifically, we tested ACSH, YPDX, and YPDX with various individual HTs (Fig. 7A; Additional file 11). In general, the Saccharomyces species had different responses to each condition (Kruskal–Wallis test p value <2.2 × 10⁻¹⁶). ACSH was the most inhibitory condition (Dunn’s test corrected p value <0.05), while YPDX plus a cocktail of HTs was the second most inhibitory (Fig. 7B). Although the latter was designed to mimic the effects of inhibitors found in ACSH, strains still grew an average of 47% more slowly in ACSH, further highlighting the chemical complexity of ACSH [11]. Strains growing in YPDX plus a single individual HT grew, on average, 9% more slowly than in YPDX without HTs (Fig. 7B). The best-growing strains of S. paradoxus, S. mikatae, S. kudriavzevii, and S. arboricola grew faster than other strains of their respective species for most individual HTs (Fig. 7A), suggesting that these strains bear novel genes or alleles that encode greater HT tolerance.

Individual HTs affected each species differently (Additional file 11). For example, compared to other Saccharomyces species, S. mikatae grew faster in most of the conditions testing individual HTs, except for in YPDX plus sodium acetate and YPDX plus the cocktail of HTs, which includes sodium acetate (Additional file 11). Thus, sodium acetate was the most inhibitory HT for our collection of S. mikatae strains (Dunn’s test, corrected p value = 0.0093). The kinetic parameters of strains were correlated among media, but the correlation values of the various media with ACSH and YPDX plus the cocktail of HTs were the lowest (Additional file 12).

Next, we examined the fermentation traits of the selected strains (Additional file 9) in ACSH microaerobic fermentations at 24 °C (to ensure that even psychrophilic (cold-tolerant) strains grew well) for 14 days, measuring several key compounds, such as extracellular sugar and ethanol concentrations (Additional file 2). To determine how HTs specifically affected these fermentations, we then selected one wild strain of each non-S. cerevisiae species, together with Y73 and Y128 (Table 1), for additional fermentations (triplicate assays) in ACSH, YPDX plus a cocktail of HTs, YPDX plus feruloyl amide, and YPDX. In these triplicate fermentations, glucose was generally consumed within the first 48 h (Fig. 7; Additional files 4, 13). After glucose was consumed, many strains utilized xylose to some degree, mostly in non-ACSH media. Xylitol accumulation was highly correlated with decreased extracellular xylose (Spearman’s correlation ρ = −1, which is a rank-based test, p value = 0.0167) (Additional file 14A, B). Glycerol accumulated during the first 46 h (Additional file 14C), while acetate accumulated between 46 and 166 h (Fig. 8). In our microaerobic conditions, ethanol levels decreased for most of the strains after the second day (Additional file 14D), suggesting that biomass continued to be produced by the consumption of accumulated ethanol and glycerol (Additional file 14E).

As expected, the S. cerevisiae engineered strain Y73 was able to utilize xylose better than other Saccharomyces strains, except for in YPDX where Y73 and Y128 consumed similar levels (Wilcoxon rank sum test p value >0.05). Although xylose utilization was highly inhibited in ACSH compared to YPDX (Wilcoxon rank sum test, p value = 2.2 × 10⁻¹⁶), we were able to detect xylose utilization by other non-S. cerevisiae strains in other media (Fig. 8; Additional file 13). In particular, S. uvarum consumed significantly more xylose (Additional file 13G) than other non-S. cerevisiae strains in YPDX plus the cocktail of HTs and in YPDX (t test, corrected p value <0.0240), except when compared to S. paradoxus in YPDX medium (Additional file 13C), where they consumed xylose similarly.

Interspecies hybridization and adaptive evolution combine and improve traits

Lignocellulosic biomass hydrolysates, such as ACSH, AFEX switchgrass hydrolysates, and wheat straw hydrolysates, are complex media [11, 64, 65] that pose particular challenges for efficient fermentations [9]. One method to begin to generate new strains with higher yields or tolerance is by selecting a chassis strain with useful properties and introducing biological parts from other organisms by genetic engineering [66]. However, this scheme requires specific knowledge of which biological parts are responsible for the traits of interest. Quantitative trait locus (QTL) mapping [67], comparative genomics [21, 68], chemical genomics [69], and several other methods have been successfully used to identify such parts, but each approach is labor intensive and costly.

Here we followed a rationale being applied in the fermented beverage industry, where several types of industrial hybrids are known, and researchers have recently constructed synthetic hybrids deliberately [70–72]. By crossing strains of S. cerevisiae engineered for xylose fermentation [21, 24] with different species of Saccharomyces, we efficiently combined their biological parts and, by adaptive evolution, established a selective regime where the hybrid strains adapted to ACSH while retaining hydrolysate tolerance and improving fermentative traits relative to the ancestral synthetic hybrids. Although none of these strains outperformed the Y73, repeating the process for additional generations may yield further improvements.

We detected several genomic rearrangements and copy number changes in these evolved hybrids, events frequently observed in Saccharomyces hybrids used in the fermented beverage industry [33, 40, 42, 43, 47, 73, 74]. At least some of these rearrangements and copy number changes probably confer adaptive properties, such as the increase in the number of copies of the Sch. stipitis xylose utilization (XYL) genes. Interestingly, the Y73 × S. kudriavzevii hybrid significantly improved its xylose consumption without amplifying these XYL genes, suggesting that its improved xylose utilization was driven by a distinct mechanism. Other rearrangements might be important for genome stability, such as the fixation of the S. cerevisiae cluster of rDNA, which is known to undergo concerted evolution in several domesticated yeasts [43, 75], as well as in interspecies hybrids in other types of organisms [76].

The genome content of hybrids from the fermented beverage industry ranges from approximately diploid to approximately tetraploid, and most strains are aneuploid [40, 47, 74]. Hybrids that are approximately triploid or tetraploid can be generated in two ways. The first model proposes haploid–haploid hybridization, followed by increases of ploidy, such as through endoreduplication [77]. A second model proposes so-called “rare mating” events where a diploid MAT a /MAT α cell is converted to MAT a /MAT a or MAT α /MAT α by gene conversion and then mates with a haploid or another compatible diploid [46, 78]. Here one of our diploid synthetic hybrids increased its genome content to approximately triploid with several tetrasomic chromosomes. Although we cannot rule out a rare mating event in liquid culture, endoreplication followed by the losses of one of the four copies of most chromosomes seems to be the simplest explanation.

Saccharomyces diversity as a source of biological innovation

Saccharomyces species offer a great opportunity to use interspecies hybridization to explore the vast genetic diversity present in the genus [30]. The likelihood that any natural strain would be optimal for all traits required for a particular industrial process is infinitesimal. Indeed, we found that none of the engineered strains, nor wild strains from the various Saccharomyces species were optimal for HT tolerance, xylose consumption, or ethanol production. This result highlights the importance of determining the genetic bases of industrially desirable traits so that they can be used to further engineer chassis strains for biofuel production.

Specifically, here we demonstrated that many strains of S. paradoxus and S. mikatae tolerate ACSH better than a strain of S. cerevisiae specifically chosen for its ACSH tolerance [24]. S. paradoxus was previously shown to tolerate furfural and hydroxymethylfurfural (HMF) [79, 80], results that our data replicate on a much broader strain collection. We also found that the under-characterized species S. mikatae might have even better tolerance traits than S. paradoxus for HMF and several other HTs. Interestingly, S. paradoxus, S. mikatae, and especially S. uvarum can also consume xylose to some extent in ACSH. In general, Saccharomyces have been found to not utilize xylose [81], but a handful of studies in S. cerevisiae [82, 83], together with this study on other Saccharomyces, have demonstrated that some strains do consume xylose, even though they do so slowly and mainly accumulate xylitol. Variation among strains in HT tolerance and xylose utilization could both be ecologically relevant. In nature, yeasts are likely to be differentially exposed to xylose and HTs, which can be derived from the degradation of lignocellulosic biomass by other fungi [84], deployed as defensive compounds by plants [85], or encountered in the guts of insects that feed on plants [86]. Thus, understanding the genetic basis of variation in these traits among yeast species and populations might simultaneously illuminate their ecologies and provide biotechnological benefits.

Exploiting the vast genetic diversity of Saccharomyces

In this study, we chose Saccharomyces strains for hybridization prior to our comprehensive screen for traits of interest to biofuel production, in part because of the considerable labor required to generate stable haploid strains for crosses. To maximize the exploration of genetic traits from Saccharomyces species and allow crosses with asexual industrial strains, we recently developed the Hybrid Production (HyPr) system [72]. The utility of HyPr resides in its ability to generate allopolyploids by efficiently inducing rare mating between diploids and selecting for hybrids with reciprocal markers. With this methodology in hand, we can now easily generate allotetraploids and expose them to adaptive evolution to generate and recombine genomic diversity. Tetraploids may be of particular interest for these types of experiments because they have been shown to adapt faster than diploids and haploids [87]. Coupled with the identification of top performers for each Saccharomyces species, this new technology will allow us to explore a broader genomic landscape than we have done here with haploid–haploid synthetic hybrids, with the goal of ultimately determining the genetic basis of variation in traits relevant to biofuel production.

Yeast strains and maintenance

Our strain survey included five hundred and seven wild Saccharomyces strains, thirteen Saccharomyces hybrids isolated from fermented beverages, and eight synthetic hybrids. Strain information can be found in Additional files 1 and 9. As a standard reference, we used a strain of S. cerevisiae engineered and aerobically evolved for the consumption of xylose, Y73 [21, 24]. Yeast strains were stored in cryotubes with YPD (1% yeast extract, 2% peptone and 2% glucose) and 15% glycerol at −80 °C. Routine cultures were maintained in YPD plus 2% agar plates at 24 °C.

Haploid-haploid crosses for hybrid production and confirmation

Haploid derivatives of S. mikatae IFO1815, S. kudriavzevii IFO1802, S. uvarum CBS7001 were previously generated [28]. Heterothallic haploid derivatives of Y73 (GLBRCY101) and S. eubayanus FM1318 (yHEB67) were generated in this study. For the generation of haploid derivative strains, the HO locus was replaced by a drug-resistance or TK [88] marker, and transformants were selected (Table 1). Sporulation was induced on sporulation plates (1% Potassium Acetate, 0.1% yeast extract, 0.05% glucose, 2% agar) or YPD at 24 °C, and tetrads were microdissected and micromanipulated using the SporePlay Dissector (Singer Instruments, UK). The strain of S. cerevisiae engineered and anaerobically evolved for the consumption of xylose, Y128 [25], is a haploid strain, so sporulation was not required. Mating type was determined by crossing the hoΔ::DominantMarker haploid strains with diagnostic auxotrophic haploid strains that were MAT a or MAT α. Successful crosses were able to grow in minimal medium (MM: 0.67% Yeast Nitrogen Base, 0.1% monosodium glutamate, 2% glucose) plus the conditions required to select for the dominant marker.

Heterothallic haploid Y101 or Y128 strains were each crossed with four haploid, MAT-compatible, non-S. cerevisiae strains. After pre-culturing strains at 24 °C overnight in YPD liquid, MAT-compatible cells were cultured in a 1.5-mL Eppendorf tube and patched in YPD agar plate to allow mating. After 12–36 h of incubation, a small amount of the patch was transferred to a YPD plate with the two corresponding drugs or other selective cocktails. DNA from successful growing colonies was extracted by the NaOH method [63, 89]. The hybrid nature of the nuclear genome was tested by a PCR-based random fragment length polymorphism (RFLP) method [42], and parental mitochondrial contribution was confirmed by PCR and Sanger-sequencing the mitochondrial COX2 gene using primers previously described [90]. Eight synthetic hybrids were isolated and stored in cryotubes with YPD plus 15% glycerol at −80 °C.

Microtiter plate growth curves

Our complete collection of strains was pre-cultured in deep 96-well plates with 500 µl of synthetic complete (SC: 0.67% Yeast Nitrogen Base without amino acids and without ammonium sulfate, 0.079% Drop-out Mix Complete, 0.5% ammonium sulfate, 2% glucose) until saturation at 24 °C. After pre-culture, 10 µl of saturated culture was inoculated into a 96-well plate (Nunc, Roskilde, Denmark) containing 240 µl of ACSH. Outer wells contained 250 µl of sterile water and were not inoculated. Each plate contained the reference strain Y73. Growth curve analysis and kinetic parameter estimation are detailed below.

Based on the results of this initial screen, the fastest-growing strain of each species or lineage and the slowest-growing strain of each species were selected, resulting in a total of forty-five Saccharomyces strains that also included the type strains of each species, the reference biofuel chassis strains, and eight synthetic hybrids (Fig. 7). The strain of S. mikatae chosen as its slowest-growing representative was discovered to not be a strain of Saccharomyces during the study, so it was removed from further analysis; all other strains were confirmed to be correctly identified. These representative strains are stored in a unique deep 96-well plate with YPD plus 15% glycerol. We calculated growth kinetic parameters in YPDX (1% yeast extract, 2% peptone, 6% glucose and 3% xylose), YPDX plus a cocktail of HTs (80 mM acetamide, 5.5 mM coumaroyl amide, 0.71 mM ferulic acid, 5.5 mM feruloyl amide, 1.1 mM hydroxy-methyl-furfural, 2.1 mM p-coumaric acid, and 32 mM sodium acetate), YPDX plus each individual HT (at concentrations indicated in the recipe for YPDX plus the cocktail of HTs). All media were adjusted to pH 5.2 with HCl. The concentrations of HTs are based their estimated concentrations in ACSH, previously published as Synthetic Hydrolysate v2 (SynH2) [11]. Feruloyl and coumaroyl amide were synthesized as described previously [11]. Each plate contained five copies of the reference strain Y73.

To monitor the growth of strains in the different conditions, inoculated 96-well plates were placed in Tecan F500 (Tecan Trading, Männedorf, Switzerland) maintaining an interior temperature set to a minimum of 28 °C. Selected synthetic ancestral and evolved hybrids (see below) and GLRBCY73 were tested at 30 °C (Fig. 1). Absorbance at 595 nm was monitored every 15 min for 5 days. Background absorbance was subtracted from the average of three negative controls (media without cells). Kinetic parameters for each condition were calculated in GCAT [91]. Average and standard deviations from two independent biological replicates were calculated in R [92]. Kinetic parameter values from the first ACSH screen were normalized against values from the reference strain, Y73.

Fermentations and extracellular compound measurement

ACSH was prepared following the previously described method [93] with autoclaving prior to enzymatic hydrolysis, and its approximate composition has been determined [11]. Exploratory ACSH fermentations of selected strains (Fig. 1; Additional file 9), S. paradoxus strain CBS432, and eight synthetic hybrids were performed in 50 mL of ACSH at 24 °C in microaerobic conditions (250-mL flasks with airlocks) to minimize the sensitivity to higher temperatures of psychrophilic Saccharomyces species, such as S. eubayanus, S. kudriavzevii, and S. uvarum [94]. The starting OD at λ = 600 nm was 0.2, and shaking was performed at 110 rpm. A 600-µl sample from each individual flask was collected from fourteen time points: 0, 22.5, 46.5, 70.5, 94.5, 118.5, 142.5, 166.5, 190.5, 214.5, 238.5, 262.5, 286.5, and 334.5 h. Absorbance at 600 nm was measured using a DiluPhotometer spectrophotometer (Implen, Munich, Germany). Samples were centrifuged at 15,000 rcf for 1 min to pellet cells. Supernatant was analyzed in high-performance liquid chromatography (HPLC, see below) to measure extracellular compounds.

Six wild non-S. cerevisiae strains, one from each species, were selected based on their kinetics and potential xylose consumption (Table 1). These promising strains, along with the two engineered S. cerevisiae strains, were selected for more detailed fermentations in ACSH, YPDX, YPDX plus a cocktail of HTs (see above), and YPDX plus 1.2 mM feruloyl amide at 24 °C for 14 days. Three biological replicates were performed in independent 150-mL flask (see above) experiments a week or more apart. Similar fermentations were performed for S. cerevisiae Y73, two synthetic hybrids, and two evolved hybrids (see below), except that they were performed at 30 °C. Samples from five time points (0, 46.5, 166.5, 262.5, and 334.5 h) were processed as above for absorbance and extracellular compounds measurements.

To assess performance of key strains in a more industrially relevant condition, anaerobic fermentations were performed in ACSH at 30 °C for Y73, two synthetic hybrids, and two evolved hybrids. Media and cultures were degassed in an anaerobic chamber where O₂ levels remained below 30 ppm. In this case, we sampled each day during the 7 days of fermentation.

Extracellular acetate, d-glucose, d-xylose, ethanol, glycerol, succinate, and xylitol in ACSH, YPDX and derivatives were determined by HPLC (Agilent 1260 infinity) using a quaternary pump, chilled (4 °C) autosampler, vacuum degasser, and refractive index detector (Agilent Technologies, Inc., Palo Alto, CA). The HPLC column consisted of an Aminex HPX-87H column (Bio-Rad) operating at 50 °C, a mobile phase of 0.02N H₂SO₄, and a flow rate of 0.5 ml/min.

Adaptive evolution

Six of the synthetic hybrids (those generated between crosses of GLRBCY101 (heterothallic haploid derivative of Y73) or Y128 with haploid derivatives of S. mikatae, S. kudriavzevii, and S. uvarum) were evolved in triplicate for 50 generations at 30 °C in ACSH in microaerobic conditions (shake flasks with airlocks). 14 days of fermentation were performed to allow cells to consume xylose, and an aliquot of each replicate was transferred to a fresh ACSH until it reached 50 generations. The absorbance of the starting cultures and pitching density was 0.2 OD at λ = 600 nm, and OD was monitored at the end of each fermentation to calculate the number of generations per round of fermentation. When necessary, we performed sample dilutions to stay in the linear range of the spectrophotometer. From each round and replicate, three samples of 1 ml of the population culture were stored at −80 °C. From fermentation round 5 to the end of adaptive evolution experiment, a sample from each individual replicate was cultured on YPD plates to isolate 10 colonies. Colonies were stored at −80 °C. When the experiments reached 50 generations, the 10 colonies were pre-cultured overnight in SC plus 0.2% glucose, and 10 µl was transferred into a 96-well plate with 240 µl of MM plus 2% xylose (Table 1) to select the fastest-growing hybrid variant for further analyses.

DNA content and genome sequencing

DNA contents of ancestral and evolved hybrids were determined by flow cytometry as previously described [72]. The diploid Y73 (2n) and the tetraploid W34/70 (4n) were used as references.

Paired-end Illumina HiSeq2000 platform 2 × 100 reads were generated (Additional file 15). Reads were demultiplexed, and adapters were removed using Trimmomatic [95], with parameters: 2:30:10 TRAILING:3 MINLEN:25. Species genome contributions to the ancestral and evolved synthetic hybrids were estimated as described previously [33] using genome assemblies of S. cerevisiae, S. mikatae, and S. kudriavzevii downloaded from www.saccharomycessensustricto.org [28]. Illumina reads were deposited in the SRA public database of NCBI under BioProject accession number SRP090190. Mitochondrial inheritance was assessed by assembling the synthetic hybrid genomes using SPAdes [96] as implemented in the wrapper iWGS v1 [97]. Scaffolds with a GC content lower than 30% were selected and mapped against the S. cerevisiae ancestor of the GLRBCY strains, Y22-3 [23], S. mikatae (accession number KX657747), and S. kudriavzevii (accession number KX657746) mitochondrial genomes, using Geneious R6 [98].

Statistical analysis

A non-parametric Kruskal–Wallis test was applied to infer significant differences among data by using a Dunn’s test, which also conducts an ad hoc test for pairwise comparisons. P values were corrected for multiple comparisons to control the false discovery rate (FDR), by applying the Benjamini-Hochberg correction. Values were considered significant when the corrected/adjusted p value was below the FDR threshold of 5%. Correlation tests were conducted using the non-parametric Spearman’s rank correlation test. The reported t test was approached as an unpaired t test because independent samples were assayed for significance. Statistics, Principal Component Analysis (PCA), boxplots, and other graphs were conducted in R (R Development Core Team 2010).