# Comprehensive survey of condition-specific reproductive isolation reveals genetic incompatibility in yeast

Jing Hou, Anne Friedrich, Jean-Sebastien Gounot, Joseph Schacherer

► To cite this version:

Jing Hou, Anne Friedrich, Jean-Sebastien Gounot, Joseph Schacherer. Comprehensive survey of condition-specific reproductive isolation reveals genetic incompatibility in yeast. Nature Communications, 2015, 6 (1), pp.7214. 10.1038/ncomms8214. hal-05022660

HAL Id: hal-05022660

https://hal.science/hal-05022660v1

Submitted on 7 Apr 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

# Comprehensive survey of condition-specific reproductive isolation reveals genetic incompatibility in yeast

Jing Hou<fnref n="1" />, Anne Friedrich<fnref n="1" />, Jean-Sebastien Gounot<fnref n="1" /> & Joseph Schacherer<fnref n="1" />

Genetic variation within a species could cause negative epistasis leading to reduced hybrid fitness and post-zygotic reproductive isolation. Recent studies in yeasts revealed chromosomal rearrangements as a major mechanism dampening intraspecific hybrid fertility on rich media. Here, by analysing a large number of *Saccharomyces cerevisiae* crosses on different culture conditions, we show environment-specific genetic incompatibility segregates readily within yeast and contributes to reproductive isolation. Over 24% (117 out of 481) of cases tested show potential epistasis, among which 6.7% (32 out of 481) are severe, with at least 20% of progeny loss on tested conditions. Based on the segregation patterns, we further characterize a two-locus Dobzhansky-Müller incompatibility case leading to offspring respiratory deficiency caused by nonsense mutation in a nuclear-encoding mitochondrial gene and tRNA suppressor. We provide evidence that this precise configuration could be adaptive in fluctuating environments, highlighting the role of ecological selection in the onset of genetic incompatibility and reproductive isolation in yeast.

<footnote n="1">¹Department of Genetics, Genomics and Microbiology, University of Strasbourg/CNRS, UMR7156, 28 rue Goethe, 67083 Strasbourg, France. Correspondence and requests for materials should be addressed to J.S.(email: schacherer@unistra.fr).</footnote>
Genetic variation accumulated in different natural popula-
tions could occasionally cause deleterious epistatic effects
leading to reduced hybrid fitness and intrinsic post-
zygotic reproductive isolation. The Dobzhansky–Müller model
described such negative epistasis as genetic incompatibility,
where independently fixed mutations in allopatric populations
could not properly function together when combined in
hybrids¹. Although most prominently studied between closely
related species²⁻⁹, recent efforts have been made to characterize
genetic incompatibilities within the same species, and it has
been shown that deleterious epistasis segregates readily at an
intraspecific scale¹⁰⁻¹⁶. The current focus on intraspecific genetic
incompatibilities underscores the importance of the ongoing
phenotypic consequences of genetic diversity within an inter-
mating population, and captures the evolutionary origin of the
early onset of reproductive isolation and speciation.

In different model organisms, multiple molecular mechanisms
have been identified to explain the observed genetic incompat-
ibility within populations. For example, in *Caenorhabditis elegans*,
the specific allelic combination of a paternal-delivered toxin
*peel-1* and a zygotic-acting antidote *zeel-1* was incompatible,
which led to F2 offspring inviability in a cross between Hawaii
and Bristol isolates¹⁰,¹⁷. Another example in *Arabidopsis thaliana*
emphasized the role of genetic drift, where reciprocal inactivation
of a duplicated essential gene histidinol-phosphate amino-
transferase (HPA) in different parental genomes led to F2
segregation distortion¹². In addition, several studies highlighted
the effect of the molecular arms race for genes involved in the
plant pathogen defence immune system leading to hybrid
necrosis¹⁸⁻²⁰. More recently, a systematic analysis across the
species *A. thaliana* further confirmed the essential role of plant
autoimmunity and identified several hot spots leading to hybrid
necrosis¹⁶. However, regardless of the molecular causes, the
evolutionary origin of incompatibilities depends largely on the life
history context of a species. The interplay between selection and
drift could play a prominent role in shaping genetic variation and
their phenotypic outcome in nature. Systematic exploration of
intraspecific reproductive isolation in additional species is
therefore valuable to obtain a global picture of the evolutionary
and molecular mechanisms that could be involved in various
species.

Yeasts are free-living unicellular eukaryotes, which can be
isolated from various ecological (for example, tree exudate, wine,
different fermentations and immuno-compromised patients) and
geographical (Europe, Asia, Africa and America) niches²¹.
Genetic variation accumulated over their evolution in these
niches constitutes a rich repertoire of potential variants leading to
reproductive isolation. Recent efforts of species-wide surveys of
post-zygotic reproductive isolation within yeasts showed frequent
occurrence of reduced hybrid fertility when crossing genetically
divergent parental isolates, where the observed cases were shown
to be mostly due to the presence of large-scale chromosomal
rearrangements that led to the unbalanced distribution of gene
sets in the offspring¹⁴,¹⁵. By contrast, no evident case of classic
genetic incompatibility has been found so far in natural
populations of yeast. Nevertheless, as most studies were only
focused on offspring viabilities in permissive laboratory
conditions, that is, rich media yeast extract peptone dextrose
(YPD) at 30°C, the extent to which different environmental
conditions has a part in the onset of reproductive isolation was
largely overlooked.

Here, we perform a global evaluation of the onset of genetic
incompatibility across various environments using the yeast
*Saccharomyces cerevisiae*. We measure the offspring viability of 27
crosses on a large number of conditions (different carbon sources,
chemicals and temperatures). We show condition-specific

incompatibilities are frequent: over one-fourth of the tested
cases show loss of offspring viability to varied degrees (from
1 to 62%), possibly indicating the presence of negative
epistasis. We analyse the segregation patterns of identified
cases and focus on one case that demonstrated a potential
recessive two loci Dobzhansky-Müller incompatibility related
to the loss of offspring viability on media containing non-
fermentable carbon sources (for example, glycerol and ethanol).
By further analysing the genes and mutations involved, we show
that the incompatibility is due to the presence of a nonsense
mutation in a nuclear-encoding mitochondrial gene and a
tRNA nonsense suppressor. This is the first identified case of
classic two loci Dobzhansky-Müller genetic incompatibility
within a yeast species. To identify the evolutionary origin of the
onset and maintenance of this particular configuration, we
evaluate the effect of the suppressor in different genetic
backgrounds and show that the presence of this suppressor
confers a surprising level of phenotypic variation, which could
potentially offer selective advantages across various environmen-
tal stresses. Interestingly, despite frequent observation of
tRNA suppressors in laboratory selections²²⁻²⁶, the occurrence
of such suppressors in natural isolates of yeast is extremely rare,
possibly suggesting a transient role of suppressors in adaptation
in *S. cerevisiae*.

Results

**Survey of environment-dependent reproductive isolation.** In total, 27 natural isolates previously shown to be compatible with the reference strain S288c (offspring viability >90% on YPD) were selected¹⁵. All isolates were crossed with S288c and offspring viability was scored and confirmed on YPD (Supplementary Table 1). For each cross, 20 full tetrads (containing only viable spores) were chosen to be tested on 20 conditions, including different temperatures, carbon sources and various chemical compounds (Fig. 1 and Supplementary Table 2). This summed up to a total of 540 instances spanning 27 crosses on 20 conditions (Fig. 1 and Supplementary Fig. 1). Among all 540 instances assessed, 59 involved at least one parental strain being non-viable on the condition tested and were excluded for further analysis (Fig. 1). Overall, 24.3% of all instances (117 out of 481) showed signs of negative epistasis with different degrees of loss of offspring viability ranging from 1 to 62% (Fig. 1). Among these cases, 6.7% (32 out of 481) showed moderate to severe incompatibility, with at least 20% of the segregants being non-viable on the condition tested (Fig. 1).

**Potential case of two loci Dobzhansky-Müller incompatibility.**

To assess the genetic complexity of the observed cases, we ana-
lysed the segregation patterns of the lethal phenotype
(Supplementary Fig. 1). Most of the cases were consistent with
complex epistasis and were not characterized in this study (~103
out of 117, Supplementary Fig. 1). We focused on one cross
between a clinical isolate YJM421 (ref. 27) and S288c, which
showed a clear pattern of recessive two loci Dobzhansky-Müller
incompatibility in several conditions related to respiration effi-
ciency (2% YP sorbitol, 2% YP glycerol, 2% YP ethanol and 2%
YP galactose). In this scenario, the lethal allelic combination
should follow Mendelian segregation, which leads to one-fourth
in the loss of viability in the offspring, resulting in a ratio of 1:4:1
for tetrads containing 4, 3 or 2 viable spores, respectively,
assuming the interacting loci are independent. For this cross,
additional tetrads were tested on 2% YP glycerol and the segre-
gation pattern was confirmed (Fig. 2a). Approximately 25% of the
offspring were respiration deficient and were unable to grow on
2% YP glycerol.
<figure><img src="image_3.png" /><figcaption>**Figure 1 | Offspring viability of 27 natural isolates crossed with S288c on 20 conditions.** Offspring viabilities estimated based on 20 full tetrads are colour coded with the vertical axis representing isolates crossed and horizontal axis representing the 20 conditions tested. All isolates were previously shown to produce high offspring viability (> 90%) on YPD when crossed with S288c¹⁵. Conditions where either one or both parental strains were non-viable are coloured in grey. The case followed is circled in black.</figcaption></figure>

**Mapping of the loci involved using bulk segregant analysis.** To map the loci involved, we used a bulk segregant analysis strategy followed by whole-genome sequencing (BSA-seq)¹⁵. Briefly, 80 segregants that were non-viable on 2% YP glycerol from independent tetrads were pooled and sequenced. The sequences obtained were aligned to the genome of S288c and the allele frequency of S288c was scored at each polymorphic position. For most genomic regions, the expected allele frequency for both parental strains was ~0.5, whereas the loci involved in the incompatibility would have deviated allele frequencies. Using this strategy, we mapped two regions with significant allele frequency deviation (see details in Methods), one located on the right arm of chromosome V (position 413,107 to 458,959) and the second on the left arm of chromosome X (position 331,633 to 364,022), spanning approximately 46 and 33 kb, respectively (Fig. 2b).

**Identification and functional validation of candidate genes.** To identify the causative genes for the observed respiratory deficiency, we closely examined the mapped regions for potential candidates involved in respiration. In total, five genes in these regions were potentially involved in respiration according to the SGD annotations (http://www.yeastgenome.org/), among which three were found in the region on chromosome V (EMP65, COX15 and FTR1) and two in the region on chromosome X (TIM54 and AIM22). We examined the DNA sequences of these genes in the YJM421 background and found a nonsense mutation at the position + 115 in the open reading frame of COX15 (CAA to TAA; position 453574 on chromosome V). COX15 encodes an inner membrane cargo protein in the mitochondria, the function of which is essential for respiration²⁸. The observation of a nonsense mutation in this gene was surprising as the presence of such mutation would likely abolish the function of COX15 and

lead to respiratory deficiency; whereas the strain YJM421 carrying the mutation was respiratory competent.

To verify if the YJM421 allele of COX15 (*cox15stop*) was functional, we deleted the *cox15stop* in YJM421, and the resulting mutant was unable to grow on media containing non-fermentable carbon sources (Supplementary Fig. 2a). Moreover, allele replacement of *cox15stop* with the wild-type COX15 from S288c in the YJM421 background resulted in total rescue of the genetic incompatibility observed: cross between YJM421 *cox15stop*: COX15 and S288c led to 98.6% offspring viability on YP glycerol 2% (400 segregants tested; Fig. 2c). These results confirmed that *cox15stop* was functional in the YJM421 background and was involved in the incompatibility between YJM421 and S288c. The fact that *cox15stop* was functionally active strongly suggests the presence of a genetic element at the interacting loci on chromosome X that compensates the effect of the nonsense mutation in YJM421.

Indeed, when examining the DNA sequence of YJM421 in the mapped region on chromosome X, we found a mutation at the anticodon position of a tyrosine tRNA tY(GUA)J1 (GTA to TTA, position 354280 on chromosome X), which in turn transformed this tRNA into a TAA nonsense suppressor (*SUP7*). The presence of this suppressor would effectively read-through the premature stop codon in *cox15stop*, which leads to a functional protein product in YJM421. However, this configuration of *cox15stop1* *SUP7* in YJM421 renders the strain incompatible when crossed with S288c, as one-fourth of the segregants would inherent only the non-functional *cox15stop* allele but not the suppressor, leading to respiratory deficiency. To confirm this hypothesis, we transformed segregants that are non-viable on 2% YP glycerol with a yeast centromeric plasmid containing the suppressor *SUP7* (*CEN SUP7*), and confirmed that the presence of this suppressor restored their respiration capacity (Supplementary Fig. 2b). These results demonstrated the first identified pair of
<figure><img src="image_4.png" /><figcaption>**Figure 2** | Molecular characterization of the incompatible case related to respiratory deficiency between YJM421 and S288c. (a) Phenotypic segregation of the cross. The frequency of tetrads containing 4, 3, 2, 1 or 0 viable segregants was presented for YPD (upper plot) and YP glycerol 2% (lower plot). The number of tetrads tested is as indicated. (b) Genomic regions with skewed allele frequencies identified using bulk segregant analysis. Eighty non-viable segregants on 2% YP glycerol were pooled and sequenced and two candidate regions were identified. The horizontal axis represents the coordinates of chromosome V and X. The vertical axis corresponds to the allele frequencies of S288c. Identified regions are coloured in red. (c) Comparison of the phenotype distributions of the segregants obtained before and after allele replacement of COX15. The normalized growth ratio of 400 segregants from the cross between YJM421 (cox15<sup>stop</sup>/SUP7; left panel) or YJM421* (COX15/SUP7; right panel) and S288c are presented as colour coded frequency distributions. Shaded areas indicate the fractions of segregants that are respiratory deficient.</figcaption></figure>

Dobzhansky–Müller incompatibility genes within a yeast species. Nevertheless, the evolutionary and physiological implications of this specific combination of *cox15<sup>stop</sup>/SUP7* are still unclear.

**Differential fitness effect of SUP7 in diverse isolates.** In fact, tRNA suppressors are well known to effectively suppress non-sense mutations by stop codon read-through, although the presence of such suppressors is likely detrimental because of the perturbation of cellular translational fidelity. As the incompatible strain YJM421 did not show any apparent growth defect, we sought to evaluate the effect of the suppressor *SUP7* on growth in different genetic backgrounds. We transformed 23 diverse natural isolates with a plasmid containing *SUP7* (CEN_SUP7) or an empty control plasmid (CEN_Ctrl) and measured their growth rate in a non-stressful condition (YPD at 30 °C) using micro-cultures (Supplementary Table 1). A mean reduction of the growth rate across all strains tested was observed (mean reduction

22.8%, N = 138, two-sided *t*-test *P*-value << 0.005), with the most severe case of 2.53-folds lower growth in the presence of the suppressor compared with the strain carrying the control plasmid (Fig. 3a; strain Y9 with two-sided *t*-test *P*-value < 0.05, N = 6). Interestingly, despite an overall deleterious effect of *SUP7*, several isolates, including YJM421, YJM320 and T7, showed similar or higher growth rates in the presence of the suppressor (Fig. 3a). These results suggest that the effect of *SUP7* on growth is background dependent and different levels of genetic assimilation could be observed, such as the case for YJM421, thus allowing for the persistence of *SUP7* in this strain.

**Impacts of SUP7 in stress conditions across natural isolates.** To further investigate the phenotypic consequences of *SUP7*, we evaluated the fitness of the same set of 23 isolates carrying the plasmid with *SUP7* (CEN_SUP7) on solid media for various stress conditions (membrane stability, proteome perturbation,
<figure><img src="image_5.png" /><figcaption>**Figure 3** | Phenotypic consequences of *SUP7* in various isolates. (a) Growth variation in the presence of *SUP7* in non-stressful conditions (YPD, 30 °C). Growth rates of 23 strains measured in liquid YPD at 30 °C in the presence of a plasmid containing *SUP7* (blue dots) or empty control vector (red dots) are presented with error bars (mean ± s.d.; *N* = 6). Percentage of growth variation is calculated and presented as bars (*N* = 6; two-sided t-test \*P value < 0.05; \*\*P value < 0.005). (b) Suppressor induced phenotypic variation in stress conditions. Significant variation of the normalized growth ratio (>10%) due to the presence of *SUP7* of 23 strains is presented for 11 stress conditions using a heatmap with blue indicating loss of growth and red for gain of growth compared with strains carrying the empty control vector. (c) Significant increase of phenotypic variance in the presence of *SUP7*. Distribution of the normalized growth ratio in stress conditions was compared for strains carrying the suppressor *SUP7* (SUP) or control (Ctrl) in seven conditions. Median values for each condition are indicated with a bar. Statistical significance is as shown on the plot (*N* = 288, two-sided F-test \*P value < 0.05, \*\*P value < 0.005, \*\*\*P value < 0.0005; Levene test †P value < 0.05).</figcaption></figure>

osmotic stress, different carbon sources and high temperatures; Fig. 3b and Supplementary Table 2). The normalized growth ratio was calculated by comparing the colony size on tested conditions versus YPD to eliminate the effects of growth differences on YPD and pinning density on solid plates. We then calculated the percentage of fitness variation for each isolate in the presence of *SUP7* compared with the same isolate carrying the control plasmid on each condition. Significant variation due to the presence of *SUP7* was observed in most of the conditions tested, with nearly half of the cases showing a gain of fitness higher than 10% (Fig. 3b). These variations appeared to be strain and condition specific, with exceptions for some conditions (YPD 37 °C and YPD ethanol 15%) where all strains grew better in the presence of *SUP7*, and some strains (CLIB294, YJM269 and CLIB272) with an overall gain of fitness across all conditions. In addition, for most of the conditions tested (7 out of 11), significantly increased phenotypic variance was observed in the presence of *SUP7* compared with the controls across all strains (Fig. 3c and

Supplementary Fig. 3). These results suggest that the suppressor *SUP7* contributes to marked phenotypic variation across different genetic backgrounds in stress conditions, and carrying the suppressor might, in turn, offer some selective advantages in the presence of environmental challenges.

**Frequency of nonsense mutation and tRNA suppressor in yeast.**
To explore the prevalence of nonsense mutations and tRNA suppressors in natural populations of yeast, we surveyed 100 genomes of *S. cerevisiae* that are publically available²⁹,³⁰. Compared with common lab strains, nonsense mutations in natural isolates were quite frequent, with an average of ~10 nonsense mutations in each stop codon class per strain (Fig. 4). The frequency of nonsense mutations globally followed a normal distribution, with a maximum frequency of nonsense mutations per strain of ~16 mutations for each anticodon class (Fig. 4). Genes bearing any class of nonsense mutations were functionally
<figure><img src="image_6.png" /><figcaption>**Figure 4 |** Nonsense mutations in natural populations. Distribution and number of nonsense mutations in verified ORFs across 100 sequenced natural isolates<sup>29,30</sup>. Mutations in different stop codon classes are colour coded.</figcaption></figure>

enriched for stress-related activities, such as transmembrane transporter activity, detoxification and transcription regulation (Supplementary Table 3). More than 40% (215 out of 500) of the detected genes with nonsense mutations were shared by at least two isolates.

In addition to nonsense mutations, we also looked for the presence of potential tRNA suppressors in these genomes. No tRNA suppressor of any anticodon class was found in this rather large data set. In contrast to the prevalence of nonsense mutations in natural populations, the frequency of tRNA suppressors is extremely rare, possibly suggesting a transient role of the suppressors in adaptation in *S. cerevisiae*.

## Discussion

We performed a species-wide survey of environment-dependent reproductive isolation and identified the first Dobzhansky–Müller incompatibility gene pair related to offspring respiratory deficiency within the yeast species *Saccharomyces cerevisiae*. We showed the incompatibility was due to a combination of a nonsense mutation in *COX15* and a tRNA suppressor *SUP7* in a single isolate, which leads to one-fourth of the offspring having only a non-functional copy of *COX15* upon crossing with the reference strain S288c. We also provided evidence that the persistence of this particular allelic combination might potentially be related to increased evolutionary potential when facing fluctuating environmental conditions in nature. Our study highlights the importance of understanding the ecological context to hybrid fitness and extends the overview of possible mechanisms involved in the onset of intraspecific post-zygotic reproductive isolation in yeast.

In yeasts, multiple mechanisms, such as chromosomal rearrangements, anti-recombination, cyto-nuclear incompatibility and meiotic drive elements, have been identified to explain the observed loss of hybrid fertility between different species<sup>7,8,31–35</sup>. However, the relative role of Dobzhansky–Müller genetic incompatibility to the onset of reproductive isolation in yeasts has long been a subject of debate, primarily due to lack of

empirical support<sup>6,36–38</sup>. At the intraspecific level, large-scale chromosomal rearrangements such as reciprocal translocations were considered to be the major mechanism leading to reduced offspring viability when crossing natural populations<sup>14,15</sup>, whereas cases of deleterious genic interactions were found to be rare, with the only example demonstrated in *S. cerevisiae* related to interactions between genes in the mismatch repair system leading to sporadic progeny loss<sup>39,40</sup>. In addition to these mechanisms, here we showed that classic two loci Dobzhansky–Müller incompatibilities do exist in natural isolates of yeast and could readily lead to reproductive isolation in different environmental conditions. Although the conditions investigated in the current study do not represent the true ecological contexts encountered, it is evident that the overall picture of molecular mechanisms affecting reproductive traits in nature is far more complex than previously envisioned within a yeast species.

Nevertheless, the evolutionary forces driving the onset and maintenance of isolating mechanisms are still in question. In the identified case of respiration-related genetic incompatibility, both mutations were found in the genome of the incompatible strain YJM421, resulting in a ‘derived-ancestral’ type of interaction. This observation confirmed that the onset of genetic incompatibility does not necessarily require independent fixation of causative mutations in allopatry, as was initially proposed by the Dobzhansky–Müller model<sup>1</sup>. However, although the origin of the onset and maintenance of this particular configuration is unclear, several possibilities could be envisioned. On one hand, as the YJM421 strain showed phenotypic tolerance of the suppressor *SUP7* (Fig. 3a), it is possible that this suppressor was acquired in conditions where possessing a suppressor was beneficial. In this scenario, the loss-of-function allele of *cox15stop* might arise due to random genetic drift, the effect of which was buffered by the pre-existence of *SUP7*. Alternatively, the fixation of the *cox15stop* could arise before the apparition of *SUP7*. In fact, the incompatible isolate YJM421 was of clinical origin and it has been shown that deletion of *COX15* could confer higher levels of resistance to antifungal drugs and biofilm formation, two traits that are particularly advantageous for clinical propagation<sup>41</sup>. In this scenario, it is possible that the loss-of-function allele *cox15stop* originally arose due to selection pressure in clinical conditions. When the strain was replaced in favourable conditions and the original selective pressure was removed, the suppressor could arise to rescue the loss of respiratory capacity because of the adaptation to a new environment. This strain could then become integrated in subsequent genetic and phenotypic assimilation, allowing the particular allelic combination to persist. Even though these scenarios remain only conjectures, it is likely that environmental fluctuation and selection might at least partly contribute to the onset and maintenance of this particular case.

tRNA suppressors perturb the translational fidelity by stop codon read-through, the effect of which resembles the yeast prion [PSI +]<sup>42,43</sup>. In *S. cerevisiae*, tRNA suppressors have been frequently selected in numerous genetic screens<sup>22,25,26</sup>. However, how frequently such suppressors occur in natural isolates was unknown. We surveyed over 100 publically available genomes of *S. cerevisiae* natural isolates and found no tRNA suppressor of any known family in these genomes except for that of YJM421. The scarcity of tRNA suppressors in natural isolates in contrast with their relatively high frequency of occurrence in the presence of strong selection suggests that tRNA suppressors might represent a transient survival mechanism, which could subsequently be lost in the absence of selection<sup>44</sup>. Nevertheless, the evolutionary fate of such suppressors probably depends on the specific genomic and environmental context of the strain in question. In fact, we demonstrated that the presence of the suppressor *SUP7* conferred an increased phenotypic capacity
across multiple stress conditions, possibly fuelled by the relatively high number of naturally occurring nonsense mutations in natural isolates. Therefore, much like the prion [PSI +]⁴⁵, tRNA suppressors could offer context-dependent selective advantages in yeast. However, as opposed to prions, tRNA suppressors are stably transmitted in a Mendelian manner, which in turn could drive the fixation of allelic combinations leading to the onset of negative epistasis, as is evident in the identified case of genetic incompatibility.

By taking into account environmental factors in the onset of reproductive isolation across a large number of crosses in S. cerevisiae, we revealed that context-dependent negative epistasis readily segregates in this species and the frequency of which might be more common than previously thought. Nevertheless, cases of the classic two loci Dobzhansky–Müller incompatibility appeared to be rare, with most cases of identified potential negative epistasis apparently reaching a higher complexity, even at an intraspecific scale. The origin of such epistasis was potentially due to the combinatorial effect of selection and drift. Further understanding of the onset of intraspecific genetic incompatibilities will extend our perspectives regarding the ongoing phenotypic consequences of genetic diversity within a species, as well as the underlying evolutionary forces that shape the patterns of such variation.

## Methods

**Strains.** A collection of 27 strains isolated from diverse ecological (tree exudate, wine, different fermentations and clinical) and geographical (Europe, Asia, Africa and America) origins was used in this study for crossing, phenotyping and survey of nonsense mutations (Supplementary Table 1). Laboratory strains isogenic to S288c, FY4 (*MAT*a) and FY5 (*MAT*a) were also used. A uracil auxotrophic mutant of YJM421 (YJM421 *ura3Δ0*) was generated by deleting the *URA3* gene using 5-FOA selection.

**Media and culture conditions.** Standard media conditions were used (YPD 1% yeast extract, 2% peptone and 2% glucose at 30 °C) for growth and maintenance of the strains. A final concentration of 200 µg ml⁻¹ hygromycin (Euromedex) was supplemented to maintain the plasmid CEN_SUP7 or CEN_Ctrl carrying a resistance marker gene hygroMX. Sporulation was induced on potassium acetate plates (1% potassium acetate, 2% agar). The 5-FOA selection plate was made by supplementing 25 µg ml⁻¹ uracil and 1 mg ml⁻¹ 5-FOA in a synthetic complete uracil dropout media (SC-URA). Detailed media composition for the screening of offspring viability and stress conditions are listed in Supplementary Table 2. All procedures are performed at 30 °C unless otherwise indicated.

**Screen of offspring viability on various conditions.** The 27 strains are all *MAT*a and were previously shown to be compatible when crossed with FY4 (*MAT*a, isogenic to S288c)¹⁵. High offspring viability (>90%) of all crosses was confirmed using microdissection on YPD using a micromanipulator (Singer MSM-400) after asci digestion by zymolase (MT ImmunO 20 T). For each cross, 20 full tetrad (containing 4 viable spores) were suspended in liquid YPD and pinned onto corresponding condition plates as well as a YPD control plate using a frogger. Growth was scored by eye after 48 h where viable segregants formed a patch. The offspring viability for each cross and condition corresponds to the ratio between the number of viable segregants and the total number of segregants viable on YPD. See Supplementary Table 1 and Supplementary Fig. 1 for detailed offspring viability and phenotype segregation for each tested case.

**Bulk segregant analysis strategy.** For the cross between YJM421 and S288c, 25% of the segregants were respiratory deficient and not viable on 2% YP glycerol. In total, 80 independent respiratory-deficient segregants were separately cultured then pooled by equal OD readings at 600 nm and the DNA was sequenced. Causative genes for the observed incompatibility are presumably united in these segregants and genomic regions involved are mapped by analysing the allele frequency variation.

**DNA extraction, sequencing, SNP calling and data analysis.** Genomic DNA was extracted using the Qiagen Genomic-tip kit. DNA of the pooled segregants was sequenced using Illumina HiSeq 2000. Detailed procedures of DNA preparation were described previously⁴⁶. We used paired-end libraries, 101 bp per read, the coverage was ~50 ×. Reads were aligned to the S288c genome using BWA with '-n 5 -o 2' options. Single-nucleotide polymorphism (SNP) calling was done using SAMtools⁴⁷. The allele frequency of S288c was scored at each polymorphic

position. As indicated by the segregation pattern, two genomic regions enriched for each parental strain are expected. We defined the region involved in the YJM421 parent as the genomic region comprising most of the mapped SNPs with allele frequency of S288c <0.1, as the S288c alleles were supposed to be absent (chromosome V). As for the region involved in the S288c genome, the YJM421 alleles should have been absent, which would result in a region with no mapped SNP markers (chromosome X). Inconsecutive SNPs that showed high allelic enrichment are considered to be sequencing noise and were not considered in the analysis.

**Plasmid construction and transformation.** Empty centromeric Gateway plasmid with hygroMX as a resistance marker was kindly provided by Dr S. Treusch (CEN_Ctrl)⁴⁸. The fragment containing *SUP7* was amplified with its original regulatory regions from the genomic DNA of YJM421, tagged with attB1/attB2 sites compatible with the Gateway cloning system (Invitrogen). Cloning of the fragment was performed according to instruction. The resulting plasmid (CEN_SUP7) was verified using restriction enzymes as well as Sanger sequencing. Transformation of yeast strains was performed using EZ transformation kit (MP biomedicals). Single colonies of transformants were maintained on YPD media containing 200 µg ml⁻¹ hygromycin.

**Growth measurements using microculture.** To measure the growth of strains in the presence or absence of the suppressor in non-stressful condition, strains transformed with the corresponding plasmids (CEN_Ctrl or CEN_SUP) were pregrown for 24 h in YPD containing 200 µg ml⁻¹ hygromycin and then transferred into 150 µl of the same media in flat bottom 96-well plates (NucIon, ThermoFischer) using a long pin replicator. Cultures were grown in six replicates for 48 h and the absorbance of each well was read at 595 nm in 10 min intervals (four spatial positions and three flashes) using a TECAN plate reader (Infinite series) with horizontal and orbital shaking. Growth curves were retrieved using the programme GATHODE⁴⁹ and the growth rate was calculated using an exponential curve fit.

**Quantitative phenotyping using solid media.** To measure growth variation due to the presence of *SUP7* in stress conditions, 23 isolates carrying different plasmids (CEN_Ctrl or CEN_SUP) were pregrown in liquid YPD with 200 µg ml⁻¹ hygromycin and pinned onto a solid YPD plate with hygromycin to a 384 format using a replicating robot RoToR (Singer instruments). The resulting matrix plate was allowed to grow for 48 h at 30 °C, and was subsequently replicated on 12 different conditions including YPD as a pinning control. All media conditions tested were supplemented with 200 µg ml⁻¹ hygromycin to maintain the plasmids, detailed compositions are listed in Supplementary Table 2. Plates were scanned for colony size at three time points (24, 40 and 48 h) and analysed using the R package Gitter⁵⁰. Each strain was tested in 12 replicates and growth was measured by analysing colony size. Normalized growth ratios were calculated first by normalizing the colony size by the growth on YPD for pinning effect, then by the value of the strain carrying the suppressor (CEN_SUP) as a control (CEN_Ctrl) for each condition. All calculations and statistical analyses were performed using R.

**Annotation of nonsense mutations and functional enrichment.** A total of 100 recently sequenced isolates of S. cerevisiae (NCBI BioProject PRJNA189847 to PRJNA189936, PRJNA189300, PRJNA188959, PRJEB2299 and refs 29,30) were used for nonsense mutation detection. The reads were retrieved and cleaned using cutAdapt, and aligned to the reference genome S288c. Read alignments and SNP calling were performed as before and SNPs were annotated using the EMBL annotation using a customized Python script. Gene Ontology (GO) term enrichment was performed using FunSpec⁵¹.

**Genome assemblies and detection of tRNA suppressor.** For the same set of 100 strains, cleaned reads were assembled using SOAPdenovo2 (ref. 52), version 2.0.4, with a k-mer size of 75. The assemblies were then surveyed for potential tRNA suppressors using tRNAscan-SE (version 1.3.1) and no suppressor of any codon families were found in this data set.

## References

1. Coyne, J. A. & Orr, H. A. *Speciation*, xiii, 545, 2 p. of plates (Sinauer Associates, 2004).

2. Presgraves, D. C. The molecular evolutionary basis of species formation. *Nat. Rev. Genet.* **11**, 175–180 (2010).

3. Chen, C. et al. A two-locus interaction causes interspecific hybrid weakness in rice. *Nat. Commun.* **5**, 3357 (2014).

4. Maheshwari, S. & Barbash, D. A. The genetics of hybrid incompatibilities. *Annu. Rev. Genet.* **45**, 331–355 (2011).

5. Ouyang, Y. & Zhang, Q. Understanding reproductive isolation based on the rice model. *Annu. Rev. Plant Biol.* **64**, 111–135 (2013).

6. Greig, D. Reproductive isolation in Saccharomyces. *Heredity (Edinb)* **102**, 39–44 (2009).
7. Chou, J. Y. & Leu, J. Y. Speciation through cytonuclear incompatibility: insights from yeast and implications for higher eukaryotes. *Bioessays* **32**, 401–411 (2010).

8. Zanders, S. *et al.* Genome rearrangements and pervasive meiotic drive cause hybrid infertility in fission yeast. *Elife* **3**, e02630 (2014).

9. Bomblies, K. & Weigel, D. Arabidopsis: a model genus for speciation. *Curr. Opin. Genet. Dev.* **17**, 500–504 (2007).

10. Seidel, H. S., Rockman, M. V. & Kruglyak, L. Widespread genetic incompatibility in *C. elegans* maintained by balancing selection. *Science* **319**, 589–594 (2008).

11. Paliwal, S., Fiumera, A. C. & Fiumera, H. L. Mitochondrial-nuclear epistasis contributes to phenotypic variation and coevolution in natural isolates of *Saccharomyces cerevisiae*. *Genetics* **198**, 1251–1265 (2014).

12. Bikard, D. *et al.* Divergent evolution of duplicate genes leads to genetic incompatibilities within *A. thaliana*. *Science* **323**, 623–626 (2009).

13. Corbett-Detig, R. B., Zhou, J., Clark, A. G., Hartl, D. L. & Ayroles, J. F. Genetic incompatibilities are widespread within species. *Nature* **504**, 135–137 (2013).

14. Charron, G., Leducq, J. B. & Landry, C. R. Chromosomal variation segregates within incipient species and correlates with reproductive isolation. *Mol. Ecol.* **23**, 4362–4372 (2014).

15. Hou, J., Friedrich, A., de Montigny, J. & Schacherer, J. Chromosomal rearrangements as a major mechanism in the onset of reproductive isolation in *Saccharomyces cerevisiae*. *Curr. Biol.* **24**, 1153–1159 (2014).

16. Chae, E. *et al.* Species-wide genetic incompatibility analysis identifies immune genes as hot spots of deleterious epistasis. *Cell* **159**, 1341–1351 (2014).

17. Seidel, H. S. *et al.* A novel sperm-delivered toxin causes late-stage embryo lethality and transmission ratio distortion in *C elegans*. *PLoS Biol.* **9**, e1001115 (2011).

18. Alcazar, R., Garcia, A. V., Parker, J. E. & Reymond, M. Incremental steps toward incompatibility revealed by Arabidopsis epistatic interactions modulating salicylic acid pathway activation. *Proc. Natl Acad. Sci. USA* **106**, 334–339 (2009).

19. Bomblies, K. & Weigel, D. Hybrid necrosis: autoimmunity as a potential gene-flow barrier in plant species. *Nat. Rev. Genet.* **8**, 382–393 (2007).

20. Alcazar, R. & Parker, J. E. The impact of temperature on balancing immune responsiveness and growth in Arabidopsis. *Trends Plant Sci.* **16**, 666–675 (2011).

21. Schacherer, J., Shapiro, J. A., Ruderfer, D. M. & Kruglyak, L. Comprehensive polymorphism survey elucidates population structure of *Saccharomyces cerevisiae*. *Nature* **458**, 342–345 (2009).

22. Adams, A. E. & Botstein, D. Dominant suppressors of yeast actin mutations that are reciprocally suppressed. *Genetics* **121**, 675–683 (1989).

23. Broach, J. R., Friedman, L. & Sherman, F. Correspondence of yeast UAA suppressors to cloned tRNAserUCA genes. *J. Mol. Biol.* **150**, 375–387 (1981).

24. Gesteland, R. F. *et al.* Yeast suppressors of UAA and UAG nonsense codons work efficiently *in vitro* via tRNA. *Cell* **7**, 381–390 (1976).

25. Novick, P., Osmond, B. C. & Botstein, D. Suppressors of yeast actin mutations. *Genetics* **121**, 659–674 (1989).

26. Ono, B. I., Tanaka, M., Kominami, M., Ishino, Y. & Shinoda, S. Recessive UAA suppressors of the yeast *Saccharomyces cerevisiae*. *Genetics* **102**, 653–664 (1982).

27. McCusker, J. H., Clemons, K. V., Stevens, D. A. & Davis, R. W. Genetic characterization of pathogenic *Saccharomyces cerevisiae* isolates. *Genetics* **136**, 1261–1269 (1994).

28. Merz, S. & Westermann, B. Genome-wide deletion mutant analysis reveals genes required for respiratory growth, mitochondrial genome maintenance and mitochondrial protein synthesis in *Saccharomyces cerevisiae*. *Genome Biol.* **10**, R95 (2009).

29. Skelly, D. A. *et al.* Integrative phenomics reveals insight into the structure of phenotypic diversity in budding yeast. *Genome Res.* **23**, 1496–1504 (2013).

30. Strope, P. K. *et al.* The 100-genomes strains, an *S. cerevisiae* resource that illuminates its natural phenotypic and genotypic variation and emergence as an opportunistic pathogen. *Genome. Res.* **25**, 762–774 (2015).

31. Fischer, G., James, S. A., Roberts, I. N., Oliver, S. G. & Louis, E. J. Chromosomal evolution in *Saccharomyces*. *Nature* **405**, 451–454 (2000).

32. Hunter, N., Chambers, S. R., Louis, E. J. & Borts, R. H. The mismatch repair system contributes to meiotic sterility in an interspecific yeast hybrid. *EMBO J.* **15**, 1726–1733 (1996).

33. Delneri, D. *et al.* Engineering evolution to study speciation in yeasts. *Nature* **422**, 68–72 (2003).

34. Chou, J. Y., Hung, Y. S., Lin, K. H., Lee, H. Y. & Leu, J. Y. Multiple molecular mechanisms cause reproductive isolation between three yeast species. *PLoS Biol.* **8**, e1000432 (2010).

35. Lee, H. Y. *et al.* Incompatibility of nuclear and mitochondrial genomes causes hybrid sterility between two yeast species. *Cell* **135**, 1065–1073 (2008).

36. Li, C., Wang, Z. & Zhang, J. Toward genome-wide identification of Bateson-Dobzhansky-Muller incompatibilities in yeast: a simulation study. *Genome Biol. Evol.* **5**, 1261–1272 (2013).

37. Greig, D. A screen for recessive speciation genes expressed in the gametes of F1 hybrid yeast. *PLoS Genet.* **3**, e21 (2007).

38. Greig, D., Borts, R. H., Louis, E. J. & Travisano, M. Epistasis and hybrid sterility in *Saccharomyces*. *Proc. Biol. Sci.* **269**, 1167–1171 (2002).

39. Demogines, A., Wong, A., Aquadro, C. & Alani, E. Incompatibilities involving yeast mismatch repair genes: a role for genetic modifiers and implications for disease penetrance and variation in genomic mutation rates. *PLoS Genet.* **4**, e1000103 (2008).

40. Heck, J. *et al.* Negative epistasis between natural variants of the *Saccharomyces cerevisiae* MLH1 and PMS1 genes results in a defect in mismatch repair. *Proc. Natl Acad. Sci. USA* **103**, 3256–3261 (2006).

41. Vandenbosch, D. *et al.* Genomewide screening for genes involved in biofilm formation and miconazole susceptibility in *Saccharomyces cerevisiae*. *FEMS Yeast Res.* **13**, 720–730 (2013).

42. Shorter, J. & Lindquist, S. Prions as adaptive conduits of memory and inheritance. *Nat. Rev. Genet.* **6**, 435–450 (2005).

43. Torabi, N. & Kruglyak, L. Genetic basis of hidden phenotypic variation revealed by increased translational readthrough in yeast. *PLoS Genet.* **8**, e1002546 (2012).

44. Yu, W. & Spreitzer, R. J. Chloroplast heteroplasmicity is stabilized by an amber-suppressor tryptophan tRNA(CUA). *Proc. Natl Acad. Sci. USA* **89**, 3904–3907 (1992).

45. Halfmann, R. *et al.* Prions are a common mechanism for phenotypic inheritance in wild yeasts. *Nature* **482**, 363–368 (2012).

46. Friedrich, A., Jung, P., Reisser, C., Fischer, G. & Schacherer, J. Population genomics reveals chromosome-scale heterogeneous evolution in a protoploid yeast. *Mol. Biol. Evol.* **32**, 184–192 (2015).

47. Li, H. *et al.* The Sequence Alignment/Map format and SAMtools. *Bioinformatics* **25**, 2078–2079 (2009).

48. Treusch, S., Albert, F. W., Bloom, J. S., Kotenko, I. E. & Kruglyak, L. Genetic mapping of MAPK-mediated complex traits across *S. cerevisiae*. *PLoS Genet.* **11**, e1004913 (2015).

49. Jung, P. P., Christian, N., Kay, D. P., Skupin, A. & Linster, C. L. Protocols and programs for high-throughput growth and aging phenotyping in yeast. *PLoS One* **10**, e0119807.

50. Wagih, O. & Parts, L. jitter: a robust and accurate method for quantification of colony sizes from plate images. *G3 (Bethesda)* **4**, 547–552 (2014).

51. Robinson, M. D., Grigull, J., Mohammad, N. & Hughes, T. R. FunSpec: a web-based cluster interpreter for yeast. *BMC Bioinformatics* **3**, 35 (2002).

52. Luo, R. *et al.* SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. *Gigascience* **1**, 18 (2012).

## Acknowledgements

We sincerely thank Bernard Dujon for helpful discussions and Kelle Freel for critical reading of the manuscript. We also thank Sebastian Treusch for kindly offering the Gateway compatible plasmid and Paul Jung for providing a demo version of the GATHODE programme for growth rate extraction. This work was supported by an Agence nationale de la recherche grant (2011-JSV6-004-01) and a National Institutes of Health (NIH) grant R01 GM101091-01 (J.S.). We also thank the Université de Strasbourg (IdEx 2012 Attraction) for their financial support. J.H. was supported by a grant from the French 'Ministère de l'Enseignement Supérieur et de la Recherche'.

## Author contributions

J.S. and J.H. conceived and designed the experiments; J.H. performed the experiments and data analyses; A.F. gathered sequence data and performed *de novo* assemblies and tRNA scans; J.-S.G. performed SNP annotations for nonsense mutations; J.H. and J.S. analysed the data and wrote the paper.

## Additional information

**Accession codes:** The short-read sequences of the pool of segregants generated in this study have been deposited in the NCBI Sequence Read Archive (SRA) database under the accession codes PRJNA280047 and SRP056771.

**Supplementary Information** accompanies this paper at http://www.nature.com/naturecommunications

**Competing financial interests:** The authors declare no competing financial interests.

**Reprints and permission information** is available online at http://npg.nature.com/reprintsandpermissions/

**How to cite this article:** Hou, J. *et al.* Comprehensive survey of condition-specific reproductive isolation reveals genetic incompatibility in yeast. *Nat. Commun.* **6**:7214 (2015).

<figure><img src="image_7.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/