Complex Genetic Interactions in a Quantitative Trait Locus

Whether in natural populations or between two unrelated members of a species, most phenotypic variation is quantitative. To analyze such quantitative traits, one must first map the underlying quantitative trait loci. Next, and far more difficult, one must identify the quantitative trait genes (QTGs), characterize QTG interactions, and identify the phenotypically relevant polymorphisms to determine how QTGs contribute to phenotype. In this work, we analyzed three Saccharomyces cerevisiae high-temperature growth (Htg) QTGs (MKT1, END3, and RHO2). We observed a high level of genetic interactions among QTGs and strain background. Interestingly, while the MKT1 and END3 coding polymorphisms contribute to phenotype, it is the RHO2 3′UTR polymorphisms that are phenotypically relevant. Reciprocal hemizygosity analysis of the Htg QTGs in hybrids between S288c and ten unrelated S. cerevisiae strains reveals that the contributions of the Htg QTGs are not conserved in nine other hybrids, which has implications for QTG identification by marker-trait association. Our findings demonstrate the variety and complexity of QTG contributions to phenotype, the impact of genetic background, and the value of quantitative genetic studies in S. cerevisiae.


Introduction
Identification of quantitative trait genes (QTGs) has been hampered by multiple factors, including variable quantitative trait locus (QTL) contributions, complex architectures [1,2], epistasis [3,4], gene-environment interactions, and genetic heterogeneity [5][6][7]. With the availability of whole genome sequences, several strategies have been proposed to facilitate the identification of QTGs [7][8][9][10][11][12], including marker-trait association, systematic functional screening, and candidate gene prediction by expression profiling or phenotype screening [13]. However, the major gap in QTL analysis is the lack of detailed molecular information for any one quantitative trait regarding QTL structure, genetic interactions, and the contributions (or relevance) of QTGs identified in one genetic background or population to other genetic backgrounds or populations of the same species.
We previously mapped a Saccharomyces cerevisiae hightemperature growth (Htg) QTL and identified MKT1, END3, and RHO2 as Htg QTGs [1]. The QTL complexity was surprising not only because there were three QTGs in the 32-kb interval, but also because the contributing alleles derived from both the Htg þ YJM145 (alleles designated as MKT1-145 and RHO2-145) and Htg À S288c (END3-288) backgrounds. In the YJM145, S288c, and 11 other Htg þ and Htg À genetic backgrounds, there were multiple coding and noncoding QTG polymorphisms; however, by marker-trait association the phenotypically relevant polymorphisms were not identifiable [1].
In this study, we assessed genes encoding genetic or physical interactors of MKT1, RHO2, and END3 as novel QTGs and found that none of these genes are Htg QTGs; we also characterized the three QTGs in isolation and in different combinations using deletion and reciprocal hemizygosity analysis (RHA; see Figure 1) [1]. We performed expression and homologous allele replacement analysis to identify phenotypically relevant polymorphisms and used RHA to determine allele contributions to Htg in ten different genetic backgrounds. Our study sheds light on the complexity underlying quantitative traits and emphasizes the utility of different approaches to dissect and analyze quantitative traits.

Known Genetic and Physical Interactors of the Htg QTGs Are Not Htg QTGs Themselves
One potential strategy to identify novel QTGs is to test genes that interact genetically (e.g., synthetic lethality) or physically (e.g., in complexes) with known QTGs. No genetic interactions have been reported for mkt1 or rho2 mutations [14,15]. However, end3 has been shown to be synthetically lethal with apm3 [16], arp2 [17], and pan1 [18] mutations and has a synthetic slow-growth phenotype with rvs161D and rvs167D deletions [19]. End3p also physically interacts with Pan1p [18] and, in addition, Mkt1p physically interacts with Pbp1p [20]. Therefore, we performed RHA in YJM145/S288c hybrids to determine if these genes, the gene products of which interact genetically or physically with the gene products of known Htg QTGs, are Htg QTGs themselves.
However, RHA showed that none of these interacting genes make allele-specific contributions to the Htg phenotype; that is, APM3, ARP2, PAN1, RVS161, RVS167, and PBP1 are not Htg QTGs.

Deletion Analysis of Htg QTGs
To determine the deletion phenotypes of Htg QTGs in the S288c and YJM145 backgrounds, we constructed single (qtg1D/ qtg1D) and double (qtg1D/qtg1D qtg2D/qtg2D) diploid deletion strains and tested their Htg phenotypes. Consistent with previous work [21,22], end3D S288c and YJM145 background strains were temperature-sensitive at 37 8C. Therefore, End3p is required for growth at even modestly elevated temperatures in both backgrounds.
In the S288c background at 39 8C, the highest temperature tested in this Htg À background, rho2D strains grew better than wild-type strains. However, in the Htg þ YJM145 background, rho2D strains grew less than wild type at 40 8C and 41 8C (Table 1). In the S288c background, mkt1D and wild-type strains grew equivalently at 39 8C. However, in the YJM145 background, mkt1D strains grew less than wild type at 40 8C, while at 41 8C mkt1D strains did not grow. Therefore, while end3D mutants in both parental backgrounds were unable to grow at 37 8C, the rho2D and mkt1D phenotypes were different in the two genetic backgrounds, which suggests that phenotype screening of null mutants may have limited predictive value as a QTG identification strategy.
In both the S288c and YJM145 backgrounds, we constructed double mutants to test for synthetic interactions between individual Htg QTGs. No synthetic lethal interactions were noted for any pair of deletions or for deletion of all three Htg QTGs in either background. All end3Dcontaining strains (end3D mkt1D, end3D rho2D, and end3D mkt1D rho2D) grew similarly to end3D strains in both backgrounds. However, in contrast to rho2D and mkt1D strains, mkt1D rho2D strains in the S288c background grew less than wild type at 39 8C, and thus exhibit a synthetic Htg À phenotype (Table 1). When tested in the YJM145 background at 40 8C, however, the growth defect of mkt1D rho2D strains was the sum of the mkt1D and rho2D growth defects, consistent with Mkt1-145p and Rho2-145p contributing additively to Htg (Table 1).

Genetic Interactions between Htg QTG Alleles
As the first step in assessing the genetic interactions between Htg QTG alleles, we determined the Htg phenotypes of a heteroallelic YJM145/S288c hybrid versus deletion heterozygotes; for example, END3-145/END3-288 versus END3-145/D, and END3-145/END3-288 versus D/END3-288. Consistent with the Htg þ alleles of MKT1, END3, and RHO2 being dominant, deletion heterozygotes containing only the Htg þ QTG alleles had Htg phenotypes similar to a heteroallelic hybrid containing both Htg þ and Htg À alleles ( Table 2). As predicted from our previous work (Table S1) [1], deletion heterozygotes containing only an Htg À allele had a highly significant growth disadvantage at high temperature relative to a heteroallelic hybrid.
Having determined the Htg phenotypes of a heteroallelic hybrid versus all single-deletion heterozygotes, we tested for interactions between Htg QTG alleles by determining the relative Htg phenotypes of a heteroallelic hybrid versus double-deletion heterozygotes. In competitions of a heteroallelic hybrid versus hybrids with the Htg À alleles of RHO2 and END3 deleted (RHO2-145/D D/END3-288), there was no

Synopsis
Most of the differences in phenotype between unrelated members of a species are polygenic in nature. Because of their ubiquity and importance, these polygenic (or quantitative) traits have been intensively studied, and a variety of techniques have been proposed to identify and characterize quantitative trait genes (QTGs). Indeed, the main application of the recently published human HapMap project is to identify the genes responsible for diseases that are quantitative in nature. Using a well-defined Saccharomyces cerevisiae quantitative trait locus containing three QTGs (MKT1, END3, and RHO2), the authors used deletions to analyze the contributions of each gene to phenotype, singly and in combination, and found a variety of interactions. Expression analysis showed no difference in steady-state mRNA levels between alleles of the three genes. Homologous allele replacement identified the phenotypically relevant differences between alleles of each gene, which were single coding polymorphisms for two genes (MKT1 and END3) and the 39 untranslated region of one gene (RHO2). Finally, analysis of multiple genetic backgrounds showed that the phenotypes conferred by these genetic variants were not conserved. The results show that the techniques proposed to identify QTGs, such as expression analysis and marker-trait association, have profound limitations, and that unbiased genome-wide approaches are needed to dissect quantitative traits. The results also demonstrate the complexity of the genetic interactions that affect quantitative traits and the value of the S. cerevisiae system in studying these traits.
difference between the hybrids. In contrast, in competitions of the heteroallelic hybrid versus hybrids with the Htg þ alleles of RHO2 and END3 deleted (D/RHO2-288 END3-145/D), the deletion heterozygotes were Htg À relative to the heteroallelic hybrid. However, the approximately two-fold growth differential in these competitions was much less than that observed when single Htg þ alleles of RHO2 and END3 were deleted ( Table 2), suggesting a complex genetic interaction between alleles of RHO2 and END3.
In competitions of a heteroallelic hybrid versus hybrids with Htg þ alleles of MKT1 and END3 deleted (D/MKT1-288 END3-145/D), the deletion heterozygotes had Htg phenotypes significantly less than the heteroallelic hybrid. The Htg phenotypic differential in these competitions was approximately equal to the sum of the differences observed in the single-deletion heterozygote analyses (Table 2), indicating MKT1-145 and END3-288 additivity. Similarly, in competitions of a heteroallelic hybrid versus hybrids with Htg þ alleles of MKT1 and RHO2 deleted (D/MKT1-288 D/RHO2-288), the deletion heterozygotes had Htg phenotypes significantly less than the heteroallelic hybrid, consistent with the earlier observed MKT1-145 and RHO2-145 additivity.

Expression Analysis of the Htg QTGs
To further characterize the Htg QTGs and gain insight into the mechanism of the quantitative trait and the location(s) of the phenotypically relevant polymorphism(s), we carefully examined the steady state mRNA levels of the three Htg QTGs as a function of both temperature and allele. Since isogenic reciprocal hemizygotes contain only one allele of a QTG, in a comparison between reciprocally hemizygous strains, any allele-specific differences in steady state mRNA levels must be due to cis-acting polymorphisms. Therefore, we performed Northern analysis of two pairs of reciprocal hemizygotes for each QTG. Temperature (25 8C versus 37 8C) had no effect on RHO2 and END3 steady state mRNA levels. However, MKT1 steady state mRNA levels increased in an allele-nonspecific manner at 37 8C (37 8C/25 8C ¼ 2.0). No allele-specific fold differences larger than two-fold were observed in steady state mRNA at either temperature for any of the genes.

Homologous Allele Replacement of the Htg QTGs
As a first test of the hypothesis that the phenotypically relevant differences between Htg þ and Htg À QTG alleles were contained within the genes, and most likely from expression results within the coding regions, we ectopically integrated alleles of RHO2 and alleles of MKT1 (separately) at the HO locus. However, we found that the resulting strains had insertspecific, Htg QTG allele-nonspecific, dominant-negative growth defects (unpublished data).
Since these ectopic integrations were problematic, we then constructed isogenic homologous allele replacement strains that differed (see Table S2) solely with respect to the S288c versus YJM145 nonsynonymous coding region polymorphisms and one 59 and one 39 UTR polymorphism in MKT1 (aÀ104g D30G K453R aþ54t), or nonsynonymous coding polymorphisms in END3 (S258N D268N) or RHO2 (F91C) (see Table S3). We then compared the Htg phenotypes of isogenic strains that differed solely in these MKT1, END3, or RHO2 polymorphisms.
In the S288c background, there was no significant difference in growth between strains with different MKT1 alleles (Tables 3 and S2). However, in the YJM145 background, Htg þ MKT1-145 allele-containing strains had highly significant Htg advantages over Htg À MKT1-288 allele-containing strains. Similarly, in the YJM145/S288c hybrid background, Htg þ MKT1-145 allele-containing strains had highly significant Htg advantages over Htg À MKT1-288 allele-containing strains. Therefore, these END3 coding (S258N D268N) and MKT1 coding/noncoding (aÀ104g D30G K453R aþ54t) polymorphisms are responsible for much of the END3 and MKT1 allele-specific contributions to the Htg phenotype, and their effects are highly influenced by genetic background.

Identification of the MKT1 and END3 Quantitative Trait Nucleotide
To determine which of the coding single-nucleotide polymorphisms (SNPs) in the two major Htg QTGs, MKT1 and END3, are relevant quantitative trait nucleotides (QTNs) for the Htg phenotype, we used site-directed mutagenesis to replace the two coding SNPs in MKT1 and END3 individually. The isogenic strains were created such that they differed solely in the coding SNP under investigation. RHA competitions were done to determine the effect of each SNP replacement in each gene on the Htg phenotype.
For MKT1, RHA analysis of SNP replacement strains showed that the A1358G SNP (where the S288c basenucleotide numbers from the start-YJM145 base), resulting in a conservative amino acid change from lysine (K) in S288c to arginine (R) in YJM145 (K453R), had no effect on Htg. However, the A89G SNP, resulting in a nonconservative amino acid change from aspartic acid (D) in S288c to glycine (G) in YJM145 (D30G), affected Htg with 89G (30G)-containing strains growing better at high temperature (Table 5).
For END3, RHA analysis of SNP replacement strains showed that the C802T SNP, resulting in an amino acid change from aspartic acid (D) in S288c to asparagine (N) in YJM145 (D268N), had no effect on Htg. However, the C773T SNP, resulting in a nonconservative amino acid change from serine (S) in S288c to asparagine (N) in YJM145 (S258N),  affected Htg with 773C (258S)-containing strains growing better at high temperature (Table 5).

Phenotypic Heterogeneity among the Htg QTGs in Diverse Genetic Backgrounds
Marker-trait association studies assess the association between a sequence variant and a trait. Marker-trait association, which may be affected by genetic background [2,3] and is one of the applications of the data coming from the HapMap project, aims to aid QTG discovery [23]. As a first step in determining the impact of genetic background on phenotype and on marker-trait association, we determined the Htg phenotype of the S288c background and of ten additional unrelated genetic backgrounds. The Htg phenotypes of the 11 unrelated genetic background strains varied over a greater than 10 4 -fold range, demonstrating a high degree of Htg phenotypic heterogeneity, and there was no association between the specific MKT1, RHO2, and END3 polymorphisms shown to be relevant in YJM145/S288c and the Htg phenotype (Table 6).
To further assess the impact of genetic background on phenotype and on marker-trait association, we crossed an S288c background strain with strains from each of the ten unrelated Htg þ and Htg À genetic backgrounds and determined the Htg phenotypes of the resulting hybrids. All hybrids, except YJM627/S288c, showed positive heterosis ( Table 6); that is, the Htg phenotype of the hybrids was greater than either parent, consistent with a general benefit of heterozygosity at many loci. Therefore, although the S288c genome constitutes one-half of the genome content of all hybrids, the genetic composition of the entire genome affects the Htg phenotype. Again, there was no association between the specific MKT1, RHO2, and END3 polymorphisms shown to be relevant in YJM145/S288c and Htg phenotype.
As a final step in determining the impact of genetic background on phenotype and on marker-trait association, we performed RHA of the Htg QTGs in the ten hybrid genetic backgrounds (Table 6) [1]. RHA of MKT1, END3, and RHO2 in hybrids between S288c and the ten unrelated genetic backgrounds further demonstrated the complex impact of genetic background on QTG contributions to phenotype. For example, in MKT1, all of the non-S288c backgrounds have the same coding polymorphism (30G) that was found to be Htg þ in the YJM145/S288c hybrid. However, in the hybrids between S288c (30D) and these other genetic backgrounds, the 30G-containing MKT1 alleles present in all of the hybrids could be phenotypically Htg þ , neutral, or Htg À and, conversely, the MKT1-288 (30D) allele present in all of the hybrids could be Htg À , neutral, or Htg þ . There are similar cases for the END3 Htg þ coding polymorphism 258S. Our analysis shows therefore that QTGs identified in one background do not necessarily make the same contribution to phenotype in other genetic backgrounds, and that QTL architecture (that is, the   identity and arrangement of QTGs) is therefore not conserved and instead varies with genetic background.

Discussion
Here we greatly extended our previous analysis [1] by determining the genetic interactions between the Htg QTGs and identifying the QTNs in the Htg genes. We determined the conservation of these QTNs in various strain backgrounds and tested various methods of QTG identification and analysis.
Our results show that not only is the architecture of this QTL complex, but that genetic interactions further add to the complexity. We found both additive interactions, such as between MKT1 and RHO2, and more complex genetic interactions, between END3 and RHO2. In addition, we showed that alleles have different effects in different genetic backgrounds. For example, END3-288 and END3-145 are phenotypically equivalent in the parental backgrounds, but differ phenotypically in the hybrid background. That is, alleles of genes from both parental backgrounds are needed to produce the END3-288 versus END3-145 phenotypic differential.
We find that none of the known interacting genes of the Htg QTGs are themselves Htg QTGs. We also find that none of the three Htg QTGs differ in their steady state mRNA levels. Indeed, our results show that RHO2 39UTR polymorphisms that have no effect on steady state mRNA levels but may instead affect mRNA translation [24,25] or localization [26] contribute to Htg; as with QTGs containing coding QTNs, such QTGs are not identifiable by expression analysis. Since the above functional and expression techniques form the basis for candidate gene approaches to identify QTGs, our results demonstrate the need for unbiased genome-wide approaches to identify QTGs.
These three Htg QTGs have no obvious functional relationship, which may have implications for other quantitative traits. END3, which has a temperature-sensitive deletion phenotype [22] and a phenotypically relevant coding polymorphism, is involved in endocytosis. RHO2 encodes a nonessential small GTPase [27]; we found that no allelespecific differences exist in RHO2 steady state mRNA levels and that the RHO2 39UTR is phenotypically relevant. MKT1 is involved in the maintenance of the M2 dsRNA satellite of the L-A virus; the S288c-allele is thought to be a loss-of-function mutation [28]. MKT1 also plays a role in the posttranscriptional regulation of HO by binding to the 39UTR of HO mRNA, without affecting steady state mRNA levels [20]. We find that the MKT1 D30G polymorphism affects Htg and others have shown that this same MKT1 polymorphism also plays a role in sporulation efficiency [29]. Since none of the mkt1 mutant phenotypes appear to be related, this represents an intriguing case of pleiotropy. This mkt1 mutant pleiotropy may be related to the function of Mkt1p in binding to the 39UTR of HO mRNA and controlling translation (without affecting steady state mRNA levels): in addition to HO, Mkt1p may regulate the translation of other mRNAs, such as RHO2.
Our analysis of Htg allele contributions in multiple genetic backgrounds demonstrates the impact of genetic background on allelic phenotype. Although our sample size was small, the lack of association of the sequence variants to the contribution of an allele has implications for marker-trait association. A variety of theoretical and experimental studies of QTLs in higher organisms [5,7,8,30] have suggested that epistasis and/or genetic background is a complicating factor in QTG identi-fication; our results provide a graphic demonstration of these effects. To fully understand a phenotype, one must understand the background effects on phenotype, which requires the identification of the genetic variants in each background that are causative. One way to identify the causative genetic variants is to ''homozygose'' known QTGs, and then map novel QTLs; another is to backcross F1 segregants and map novel Htg QTLs. These strategies will assist in identifying additional Htg QTLs that contribute to phenotype in the presence and absence of the identified Htg QTGs.
Our results show the importance of yeast as a model system for QTL analysis and highlight important issues related to QTL identification and analysis. Our results and conclusions from quantitative studies in yeast show that the study of these phenotypes in higher organisms remains a difficult challenge. Improved understanding of quantitative traits in simpler eukaryotes such as yeast will facilitate the analysis of complex traits in higher eukaryotes.

Materials and Methods
Yeast strains, growth media, and conditions. Yeast strains were grown at 30 8C unless otherwise stated. The S. cerevisiae strains used in this work (Table S4) are derivatives of S288c (a laboratory genetic background, designated ''S''), YJM145 (a clinically-derived genetic background, designated ''YBN,'' ''YJM,'' or ''YHS'' [31]), and other clinically and nonclinically derived genetic backgrounds (designated YHS, YRP, and YJM). To make hybrids (designated as ''XHS'') of test and S288c background strains, spores of lys2 and lys5 diploid strains were allowed to mate on YPD plates that were then replica plated to SDC-LYS plates to select Lys þ diploid hybrids.
To induce sporulation, cultures were first grown overnight in liquid YPD at 30 8C, after which a 1-ml aliquot was removed, washed in sterile water twice, and resuspended in 3 ml of sterile 0.5% (w/v) potassium acetate medium [36]. Sporulating cultures were incubated at 30 8C for 48 h with aeration in a roller drum and examined microscopically for the presence of tetrads.
High-temperature growth competitions between two or three strains were performed using a technique similar to RHA, as previously described [1].
PCR and yeast transformation. The PCR to make gene deletion constructs containing dominant drug-resistance markers was done as described previously [35]. All PCR products were purified with Qiagen's PCR Purification kit (Qiagen, Valencia, California, United States) before introduction into yeast. DNA was introduced into yeast by the lithium acetate transformation protocol [37]. Homologous integration of dominant drug-resistance cassettes was verified by colony PCR [38]. In cases where colony PCR did not yield a product, genomic DNA was isolated [39] and PCR was performed using colony PCR conditions. The primers used in the study (IDT, Corwallis, Illinois, United States) are listed in Table S5.
Reciprocal hemizygosity analysis. RHA, a technique for identifying allele-specific contributions to a phenotype, was performed as previously described (see Figure 1) [1]. The resulting data were used to calculate means and standard errors and to perform tests of significance [40].
Reciprocal hemizygosity expression analysis. To assay the expression levels of each allele of the Htg quantitative trait genes, we performed Northern analysis of isogenic reciprocal hemizygotes. For each Htg quantitative trait gene, a pair of reciprocal hemizygotes were grown (separately) at 25 8C and 37 8C in 50 ml of YPD with aeration to an OD 600 of 1.0. Total RNA was extracted as previously described, with the following modifications. After hot phenol-chloroform extraction, samples were frozen overnight at À80 8C and then centrifuged for 15 min at 13,000 rpm; the aqueous phase was processed as previously described. Cells were then washed and resuspended in 400 ll of extraction buffer (50 mM sodium acetate [pH 5.3], 1 mM EDTA) and 40 ll of SDS (10% w/v) followed by extraction with 440 ll of hot buffered phenol-chloroform (65 8C) for 30 min with vortexing every 5 min. After being frozen overnight at À80 8C, samples were centrifuged for 15 min at 13,000 rpm. After transfer of the aqueous phase into a new tube, the samples were phenol-chloroform extracted once. The RNA was precipitated with two volumes of 100% ethanol and 0.1 volume of 3 M sodium acetate (pH 5.3) for 30 min at À20 8C. After centrifugation, the RNA was resuspended in 100 ll of DEPC-treated water. For each sample, 10 lg of RNA was size-fractionated on a MOPS/formaldehyde gel and Northern blotting was performed. The membranes were probed with 32 P-labeled gene-specific PCR probe in ULTRAhyb solution (Ambion, Austin, Texas, United States) and were incubated overnight at 42 8C. The blots were washed twice in 23 SSC and 0.1% SDS for 5 min and twice again with 0.13 SSC and 0.1% SDS for 15 min, exposed to a phosphor-imager screen that was read with a Storm imaging system and quantified using ImageQuant software (Amersham Biosciences, Piscataway, New Jersey, United States). The mRNA levels of the Htg quantitative trait genes were compared against actin mRNA as a loading control.
Cloning alleles of quantitative trait genes. The Htg þ and Htg À alleles of MKT1-END3 (pHS07 and pHS08; see Figure S1) and RHO2 (pHS09 and pHS10; see Figure S2) were cloned using a modified gap repair [41] technique. To gap repair cloned alleles of each gene, primers with 59 tails homologous to the pRS316 [42] multicloning site were used to amplify from genomic DNA 0.5-kb regions of the open reading frames of the genes upstream and downstream of the Htg QTG. The PCR products were then purified using a Qiagen PCR Purification kit and introduced into ura3 strains (S19 and YJM791, Table S4) along with SmaI-digested pRS316 plasmid. Circularization of the SmaI-digested plasmid occurred via in vivo recombination, with the upstream and downstream PCR products providing homology between the multicloning site and the chromosomal regions flanking the Htg gene. Plasmids from Ura þ transformants, which were selected on SDC-URA plates, were introduced into the E. coli strain DH5a followed by selection on LB (containing 100 lg/ml ampicillin and 40 lg/ml X-gal) [43,44]. Plasmids were isolated from white, ampicillinresistant E. coli transformants using the Qiagen Miniprep kit and were checked by restriction digestion for the presence of the desired inserts. The inserts were then sequenced to confirm their sequence and the presence of the allele-specific polymorphisms.
Homologous coding allele replacements of Htg QTGs. Homologous allele replacements were done in haploid strains in three different ways for the three Htg genes. For MKT1 coding and noncoding allele replacements, a LYS5 counterselection method described by Ito-Harashima and McCusker [45] was used. A region of the MKT1 gene from 69 to 1,379 bp was first deleted in haploid lys5 strains (YHS731 and S1288) using the LYS5MX4 cassette to generate mkt1D::LYS5MX4 strains. Into these mkt1D::LYS5MX4 strains were introduced (i) a 3,344-bp BcgI-BsrGI MKT1 allele fragment from either pHS07 (CEN URA3 MKT1-288 END3-288) or pHS08 (CEN URA3 MKT1-145 END3-145) (see Figure S1) and (ii) pSA41 (CEN NatMX4) (see Figure S3); cells were plated on YPDþNat plates. Replacements of the mkt1D::LYS5MX4 by the cloned MKT1 DNA were selected by replica plating the colonies growing on YPDþNat onto aaminoadipate plates; replica plating to a-aminoadipate plates was repeated to reduce background growth. Lys À colonies growing on aaminoadipate plates were tested by PCR for the desired allele replacement, e.g., mkt1D::LYS5MX4 lys5 ! MKT1-145 lys5. DNA sequencing was performed to confirm the presence of the desired allelic polymorphism (see Table S2).
RHO2 and MKT1 coding allele replacements were performed in the manner similar to the MKT1 coding and noncoding allele replacements except that a Hyg/T-urf13MX4 cassette (in pBN49, see Figure S4) was used instead of LYS5MX4 for selection and counterselection. The Hyg/T-urf13MX4 cassette is a counterselectable drug marker cassette in which hygromycin resistance is used to select transformants, after which the cassette can be counterselected on methomyl-containing medium; T-urf13 is a maize mitochondrial gene that confers sensitivity to methomyl [46]; a similar pop-in/pop-out gene replacement strategy has been described in Pichia pastoris [47]. For MKT1, the polymorphic regions of 69-104 bp, containing the coding polymorphism at 89 bp, and 1,339-1,379, containing the coding polymorphism at 1,358 bp, were deleted separately using the Hyg/T-urf13MX4 cassette. For RHO2, the polymorphic region of 216-324 bp containing the coding polymorphism at 271 bp was deleted in haploid strains using the Hyg/T-urf13MX4 cassette. Using a technique similar to that described by Storici et al. [48], we introduced a 120-bp dsDNA containing the appropriate allele-specific polymorphism into these hygromycin-resistant methomyl-sensitive rho2D::Hyg/T-urf13MX4 transformants. We then selected for loss of Hyg/T-urf13MX4 cassette on YPEG þ 1 mM methomyl plates. The putative allele replacements were checked for hygromycin sensitivity and by PCR. DNA sequencing was performed to confirm the presence of the desired polymorphisms.
END3 allele replacements were based on rescuing the temperaturesensitive phenotype of end3D strains. First, the polymorphic region of END3 (773-802 bp) was deleted using the KanMX4 cassette to generate end3D::KanMX4 strains. Into the end3D::KanMX4 strains was introduced (i) a 120-bp dsDNA containing allele-specific polymorphisms using a technique similar to that described by Storici et al. [48] and (ii) pAG26 [35] (CEN HygMX4); cells were plated onto YPDþHyg and incubated at 30 8C. The hygromycin-resistant transformants were then replica plated onto YPD plates, which were incubated at 39 8C to select for END3 þ transformants-that is, transformants in which end3D::KanMX4 has been replaced by the dsDNA repairing the deletion and generating END3 þ . The putative END3 þ transformants were checked for Kan sensitivity and temperature resistance, and were confirmed by PCR and DNA sequencing.