Dissecting the pleiotropic consequences of a quantitative trait nucleotide

Role of genetic architecture in phenotypic plasticity

Phenotypic plasticity, the ability of an organism to display different phenotypes across environments, is widespread in nature. Plasticity aids survival in novel environments. Herein, we review studies from yeast that allow us to start uncovering the genetic architecture of phenotypic plasticity. Genetic variants and their interactions impact the phenotype in different environments, and distinct environments modulate the impact of genetic variants and their interactions on the phenotype. Because of this, certain hidden genetic variation is expressed in specific genetic and environmental backgrounds. A better understanding of the genetic mechanisms of phenotypic plasticity will help to determine short- and long-term responses to selection and how wide variation in disease manifestation occurs in human populations.

Multiple epistatic DNA variants in a single gene affect gene expression in trans

DNA variants that alter gene expression in trans are important sources of phenotypic variation. Nevertheless, the identity of trans-acting variants remains poorly understood. Single causal variants in several genes have been reported to affect the expression of numerous distant genes in trans. Whether these simple molecular architectures are representative of trans-acting variation is unknown. Here, we studied the large RAS signaling regulator gene IRA2, which contains variants with extensive trans-acting effects on gene expression in the yeast Saccharomyces cerevisiae. We used systematic CRISPR-based genome engineering and a sensitive phenotyping strategy to dissect causal variants to the nucleotide level. In contrast to the simple molecular architectures known so far, IRA2 contained at least seven causal nonsynonymous variants. The effects of these variants were modulated by nonadditive, epistatic interactions. Two variants at the 5'-end affected gene expression and growth only when combined with a third variant that also had no effect in isolation. Our findings indicate that the molecular basis of trans-acting genetic variation may be considerably more complex than previously appreciated.

A neglected conceptual problem regarding phenotypic plasticity's role in adaptive evolution: The importance of genetic covariance and social drive

There is tantalizing evidence that phenotypic plasticity can buffer novel, adaptive genetic variants long enough to permit their evolutionary spread, and this process is often invoked in explanations for rapid adaptive evolution. However, the strength and generality of evidence for it is controversial. We identify a conceptual problem affecting this debate: recombination, segregation, and independent assortment are expected to quickly sever associations between genes controlling novel adaptations and genes contributing to trait plasticity that facilitates the novel adaptations by reducing their indirect fitness costs. To make clearer predictions about this role of plasticity in facilitating genetic adaptation, we describe a testable genetic mechanism that resolves the problem: genetic covariance between new adaptive variants and trait plasticity that facilitates their persistence within populations. We identify genetic architectures that might lead to such a covariance, including genetic coupling via physical linkage and pleiotropy, and illustrate the consequences for adaptation rates using numerical simulations. Such genetic covariances may also arise from the social environment, and we suggest the indirect genetic effects that result could further accentuate the process of adaptation. We call the latter mechanism of adaptation social drive, and identify methods to test it. We suggest that genetic coupling of plasticity and adaptations could promote unusually rapid ‘runaway’ evolution of novel adaptations. The resultant dynamics could facilitate evolutionary rescue, adaptive radiations, the origin of novelties, and other commonly studied processes.

Sulfate assimilation regulates hydrogen sulfide production independent of lifespan and reactive oxygen species under methionine restriction condition in yeast

Endogenously produced hydrogen sulfide was proposed to be an underlying mechanism of lifespan extension via methionine restriction. However, hydrogen sulfide regulation and its beneficial effects via methionine restriction remain elusive. Here, we identified the genes required to increase hydrogen sulfide production under methionine restriction condition using genome-wide high-throughput screening in yeast strains with single-gene deletions. Sulfate assimilation-related genes, such as MET1, MET3, MET5, and MET10, were found to be particularly crucial for hydrogen sulfide production. Interestingly, methionine restriction failed to increase hydrogen sulfide production in mutant strains; however, it successfully extended chronological lifespan and reduced reactive oxygen species levels. Altogether, our observations suggested that increased hydrogen sulfide production via methionine restriction is not the mechanism underlying extended yeast lifespan, even though increased hydrogen sulfide production occurred simultaneously with yeast lifespan extension under methionine restriction condition.

DNA variants affecting the expression of numerous genes in trans have diverse mechanisms of action and evolutionary histories

DNA variants that alter gene expression contribute to variation in many phenotypic traits. In particular, trans-acting variants, which are often located on different chromosomes from the genes they affect, are an important source of heritable gene expression variation. However, our knowledge about the identity and mechanism of causal trans-acting variants remains limited. Here, we developed a fine-mapping strategy called CRISPR-Swap and dissected three expression quantitative trait locus (eQTL) hotspots known to alter the expression of numerous genes in trans in the yeast Saccharomyces cerevisiae. Causal variants were identified by engineering recombinant alleles and quantifying the effects of these alleles on the expression of a green fluorescent protein-tagged gene affected by the given locus in trans. We validated the effect of each variant on the expression of multiple genes by RNA-sequencing. The three variants differed in their molecular mechanism, the type of genes they reside in, and their distribution in natural populations. While a missense leucine-to-serine variant at position 63 in the transcription factor Oaf1 (L63S) was almost exclusively present in the reference laboratory strain, the two other variants were frequent among S. cerevisiae isolates. A causal missense variant in the glucose receptor Rgt2 (V539I) occurred at a poorly conserved amino acid residue and its effect was strongly dependent on the concentration of glucose in the culture medium. A noncoding variant in the conserved fatty acid regulated (FAR) element of the OLE1 promoter influenced the expression of the fatty acid desaturase Ole1 in cis and, by modulating the level of this essential enzyme, other genes in trans. The OAF1 and OLE1 variants showed a non-additive genetic interaction, and affected cellular lipid metabolism. These results demonstrate that the molecular basis of trans-regulatory variation is diverse, highlighting the challenges in predicting which natural genetic variants affect gene expression.

Differences in the DNA sequence of individual genomes contribute to differences in many traits, such as appearance, physiology, and the risk for common diseases. An important group of these DNA variants influences how individual genes across the genome are turned on or off. In this paper, we describe a strategy for identifying such “trans-acting” variants in different strains of baker’s yeast. We used this strategy to reveal three single DNA base changes that each influences the expression of dozens of genes. These three DNA variants were very different from each other. Two of them changed the protein sequence, one in a transcription factor and the other in a sugar sensor. The third changed the expression of an enzyme, a change that in turn caused other genes to alter their expression. One variant existed in only a few yeast isolates, while the other two existed in many isolates collected from around the world. This diversity of DNA variants that influence the expression of many other genes illustrates how difficult it is to predict which DNA variants in an individual’s genome will have effects on the organism.

Quantitative Trait Nucleotides Impacting the Technological Performances of Industrial Saccharomyces cerevisiae Strains

The budding yeast Saccharomyces cerevisiae is certainly the prime industrial microorganism and is related to many biotechnological applications including food fermentations, biofuel production, green chemistry, and drug production. A noteworthy characteristic of this species is the existence of subgroups well adapted to specific processes with some individuals showing optimal technological traits. In the last 20 years, many studies have established a link between quantitative traits and single-nucleotide polymorphisms found in hundreds of genes. These natural variations constitute a pool of QTNs (quantitative trait nucleotides) that modulate yeast traits of economic interest for industry. By selecting a subset of genes functionally validated, a total of 284 QTNs were inventoried. Their distribution across pan and core genome and their frequency within the 1,011 Saccharomyces cerevisiae genomes were analyzed. We found that 150 of the 284 QTNs have a frequency lower than 5%, meaning that these variants would be undetectable by genome-wide association studies (GWAS). This analysis also suggests that most of the functional variants are private to a subpopulation, possibly due to their adaptive role to specific industrial environment. In this review, we provide a literature survey of their phenotypic impact and discuss the opportunities and the limits of their use for industrial strain selection.

Novel biosynthetic pathway for sulfur amino acids in Cryptococcus neoformans

We elucidated a unique feature of sulfur metabolism in Cryptococcus neoformans. C. neoformans produces cysteine solely by the O-acetylserine pathway that consists of serine-O-acetyl transferase and cysteine synthase. We designated the gene encoding the former enzyme CYS2 (locus tag CNE02740) and the latter enzyme CYS1 (locus tag CNL05880). The cys1Δmutant strain was found to be avirulent in a murine infection model. Methionine practically does not support growth of the cys1Δ strain, and cysteine does not serve as a methionine source, indicating that the transsulfuration pathway does not contribute to sulfur amino acid synthesis in C. neoformans. Among the genes encoding enzymes catalyzing the reactions from homoserine to methionine, the gene corresponding to the Saccharomyces cerevisiae MET17 encoding O-acetylhomoserine sulfhydrylase (Met17p) had remained to be identified in C. neoformans. By genetic analysis of Met^- mutants obtained by Agrobacterium tumefaciens-mediated mutagenesis, we concluded that Cnc01220, most similar to Str2p (36% identity), cystathionine-γ-synthase, in the Saccharomyces genome, is the C. neoformans version of O-acetylhomoserine sulfhydrylase. We designated CNC01220 as MET17. The C. neoformans met3Δ mutant defective in the first step of the sulfate assimilation pathway, sulfate adenylyltransferase, barely uses methionine as a sulfur source, whereas it uses cysteine efficiently. The poor utilization of methionine by the met3Δ mutant is most probably due to the absence of the transsulfuration pathway, causing an incapability of C. neoformans to produce cysteine and hydrogen sulfide from methionine. When cysteine is used as a sulfur source, methionine is likely produced de novo by using hydrogen sulfide derived from cysteine via an unidentified pathway. Altogether, the unique features of sulfur amino acid metabolism in C. neoformans will make this fungus a valuable experimental system to develop anti-fungal agents and to investigate physiology of hydrogen sulfide.

Population perspectives on functional genomic variation in yeast

Advances in high-throughput sequencing have facilitated large-scale surveys of genomic variation in the budding yeast,Saccharomyces cerevisiae These surveys have revealed extensive sequence variation between yeast strains. However, much less is known about how such variation influences the amount and nature of variation for functional genomic traits within and between yeast lineages. We review population-level studies of functional genomic variation, with a particular focus on how population functional genomic approaches can provide insights into both genome function and the evolutionary process. Although variation in functional genomics phenotypes is pervasive, our understanding of the consequences of this variation, either in physiological or evolutionary terms, is still rudimentary and thus motivates increased attention to appropriate null models. To date, much of the focus of population functional genomic studies has been on gene expression variation, but other functional genomic data types are just as likely to reveal important insights at the population level, suggesting a pressing need for more studies that go beyond transcription. Finally, we discuss how a population functional genomic perspective can be a powerful approach for developing a mechanistic understanding of the processes that link genomic variation to organismal phenotypes through gene networks.

Temporal expression profiling identifies pathways mediating effect of causal variant on phenotype

Even with identification of multiple causal genetic variants for common human diseases, understanding the molecular processes mediating the causal variants’ effect on the disease remains a challenge. This understanding is crucial for the development of therapeutic strategies to prevent and treat disease. While static profiling of gene expression is primarily used to get insights into the biological bases of diseases, it makes differentiating the causative from the correlative effects difficult, as the dynamics of the underlying biological processes are not monitored. Using yeast as a model, we studied genome-wide gene expression dynamics in the presence of a causal variant as the sole genetic determinant, and performed allele-specific functional validation to delineate the causal effects of the genetic variant on the phenotype. Here, we characterized the precise genetic effects of a functional MKT1 allelic variant in sporulation efficiency variation. A mathematical model describing meiotic landmark events and conditional activation of MKT1 expression during sporulation specified an early meiotic role of this variant. By analyzing the early meiotic genome-wide transcriptional response, we demonstrate an MKT1-dependent role of novel modulators, namely, RTG1/3, regulators of mitochondrial retrograde signaling, and DAL82, regulator of nitrogen starvation, in additively effecting sporulation efficiency. In the presence of functional MKT1 allele, better respiration during early sporulation was observed, which was dependent on the mitochondrial retrograde regulator, RTG3. Furthermore, our approach showed that MKT1 contributes to sporulation independent of Puf3, an RNA-binding protein that steady-state transcription profiling studies have suggested to mediate MKT1-pleiotropic effects during mitotic growth. These results uncover interesting regulatory links between meiosis and mitochondrial retrograde signaling. In this study, we highlight the advantage of analyzing allele-specific transcriptional dynamics of mediating genes. Applications in higher eukaryotes can be valuable for inferring causal molecular pathways underlying complex dynamic processes, such as development, physiology and disease progression.

The causal path from a genetic variant to a complex phenotype such as disease progression is often not known. Studying gene expression variation is one approach to identify the mediating genes, however, it is difficult to distinguish causative from correlative genes. This becomes a challenge especially when studying developmental and physiological traits, since they involve dynamic processes contributing to the variation and only single static expression profiling is performed. As a proof of concept, we addressed this challenge here in yeast, by studying genome-wide gene expression in the presence of the causative polymorphism of MKT1 as the sole genetic variant, during the time phase when it contributes to sporulation efficiency variation. Our analysis during early sporulation identified mitochondrial retrograde signaling and nitrogen starvation as novel regulators, acting additively to regulate sporulation efficiency. Furthermore, we showed that PUF3, a known interactor of MKT1 had an independent role in sporulation. Our results highlight the role of differential mitochondrial signaling for efficient meiosis, providing insights into the factors regulating infertility. In addition, our study has implications for characterizing the molecular effects of causal genetic variants on dynamic biological processes during development and disease progression.

Phytochelatin synthase is required for tolerating metal toxicity in a basidiomycete yeast and is a conserved factor involved in metal homeostasis in fungi

Background

Phytochelatin synthase (PCS) is an enzyme that catalyzes the biosynthesis of phytochelatin from glutathione. Phytochelatins protect cells against the toxic effects of non-essential heavy metals, such as cadmium, and hence growth is restricted in the presence of these metals in mutants in PCS-encoding genes. PCS genes from fungi have been characterized in only two species in the Ascomycota, and these genes are considered sparsely distributed in the fungal kingdom.

Results

A gene encoding a putative PCS was identified in Sporobolomyces sp. strain IAM 13481, a fungus that is a member of the Pucciniomycotina subphylum of the Basidiomycota. The function of this PCS1 gene was assessed by heterologous expression in the Ascomycota yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe, and by mutating the gene in Sporobolomyces. The gene is required for tolerance to toxic concentrations of non-essential cadmium as well as the essential metal copper. Pcs1 homologs in fungi and other eukaryotes have putative targeting sequences for mitochondrial localization: the S. pombe homolog was fused to green fluorescent protein and it co-localized with a mitochondrial dye. Evaluation of the presence or absence of PCS and PCS-like homologs in the genome sequences of fungi indicates that they have a wide distribution, and the absence in most Ascomycota and Basidiomycota (the Dikarya) species can be explained by a small number of gene losses.

Conclusions

The ecology of the species within the fungi carrying putative PCS genes, the phenotypes of phytochelatin synthase mutants in two major fungal lineages, and the presence of homologs in many non-Dikarya lineages parallel what is seen in the plant and animal kingdoms. That is, PCS is a protein present early during the evolution of the fungi and whose role is not solely dedicated to combating toxic concentrations of non-essential metals.

Electronic supplementary material

The online version of this article (doi:10.1186/s40694-015-0013-3) contains supplementary material, which is available to authorized users.

Formation of hydrogen sulfide from cysteine in Saccharomyces cerevisiae BY4742: genome wide screen reveals a central role of the vacuole

Discoveries on the toxic effects of cysteine accumulation and, particularly, recent findings on the many physiological roles of one of the products of cysteine catabolism, hydrogen sulfide (H₂S), are highlighting the importance of this amino acid and sulfur metabolism in a range of cellular activities. It is also highlighting how little we know about this critical part of cellular metabolism. In the work described here, a genome-wide screen using a deletion collection of Saccharomyces cerevisiae revealed a surprising set of genes associated with this process. In addition, the yeast vacuole, not previously associated with cysteine catabolism, emerged as an important compartment for cysteine degradation. Most prominent among the vacuole-related mutants were those involved in vacuole acidification; we identified each of the eight subunits of a vacuole acidification sub-complex (V₁ of the yeast V-ATPase) as essential for cysteine degradation. Other functions identified included translation, RNA processing, folate-derived one-carbon metabolism, and mitochondrial iron-sulfur homeostasis. This work identified for the first time cellular factors affecting the fundamental process of cysteine catabolism. Results obtained significantly contribute to the understanding of this process and may provide insight into the underlying cause of cysteine accumulation and H₂S generation in eukaryotes.

Genetic architecture of ethanol-responsive transcriptome variation in Saccharomyces cerevisiae strains

Natural variation in gene expression is pervasive within and between species, and it likely explains a significant fraction of phenotypic variation between individuals. Phenotypic variation in acute systemic responses can also be leveraged to reveal physiological differences in how individuals perceive and respond to environmental perturbations. We previously found extensive variation in the transcriptomic response to acute ethanol exposure in two wild isolates and a common laboratory strain of Saccharomyces cerevisiae. Many expression differences persisted across several modules of coregulated genes, implicating trans-acting systemic differences in ethanol sensing and/or response. Here, we conducted expression QTL mapping of the ethanol response in two strain crosses to identify the genetic basis for these differences. To understand systemic differences, we focused on "hotspot" loci that affect many transcripts in trans. Candidate causal regulators contained within hotspots implicate upstream regulators as well as downstream effectors of the ethanol response. Overlap in hotspot targets revealed additive genetic effects of trans-acting loci as well as "epi-hotspots," in which epistatic interactions between two loci affected the same suites of downstream targets. One epi-hotspot implicated interactions between Mkt1p and proteins linked to translational regulation, prompting us to show that Mkt1p localizes to P bodies upon ethanol stress in a strain-specific manner. Our results provide a glimpse into the genetic architecture underlying natural variation in a stress response and present new details on how yeast respond to ethanol stress.

MET2 affects production of hydrogen sulfide during wine fermentation

The production of hydrogen sulfide (H2S) during yeast fermentation contributes negatively to wine aroma. We have mapped naturally occurring mutations in commercial wine strains that affect production of H2S. A dominant R310G mutant allele of MET2, which encodes homoserine O-acetyltransferase, is present in several wine yeast strains as well as in the main lab strain S288c. Reciprocal hemizygosity and allele swap experiments demonstrated that the MET2 R310G allele confers reduced H2S production. Mutations were also identified in genes encoding the two subunits of sulfite reductase, MET5 and MET10, which were associated with reduced H2S production. The most severe of these, an allele of MET10, showed five additional phenotypes: reduced growth rate on sulfate, elevated secretion of sulfite, and reduced production in wine of three volatile sulfur compounds: methionol, carbon disulfide and methylthioacetate. Alleles of MET5 and MET10, but not MET2, affected H2S production measured by colour assays on BiGGY indicator agar, but MET2 effects were seen when bismuth was added to agar plates made with Sauvignon blanc grape juice. Collectively, the data are consistent with the hypothesis that H2S production during wine fermentation results predominantly from enzyme activity in the sulfur assimilation pathway. Lower H2S production results from mutations that reduce the activity of sulfite reductase, the enzyme that produces H2S, or that increase the activity of L-homoserine-O-acetyltransferase, which produces substrate for the next step in the sulfur assimilation pathway.

The evolution of phenotypic correlations and "developmental memory"

Development introduces structured correlations among traits that may constrain or bias the distribution of phenotypes produced. Moreover, when suitable heritable variation exists, natural selection may alter such constraints and correlations, affecting the phenotypic variation available to subsequent selection. However, exactly how the distribution of phenotypes produced by complex developmental systems can be shaped by past selective environments is poorly understood. Here we investigate the evolution of a network of recurrent nonlinear ontogenetic interactions, such as a gene regulation network, in various selective scenarios. We find that evolved networks of this type can exhibit several phenomena that are familiar in cognitive learning systems. These include formation of a distributed associative memory that can "store" and "recall" multiple phenotypes that have been selected in the past, recreate complete adult phenotypic patterns accurately from partial or corrupted embryonic phenotypes, and "generalize" (by exploiting evolved developmental modules) to produce new combinations of phenotypic features. We show that these surprising behaviors follow from an equivalence between the action of natural selection on phenotypic correlations and associative learning, well-understood in the context of neural networks. This helps to explain how development facilitates the evolution of high-fitness phenotypes and how this ability changes over evolutionary time.

The molecular basis of phenotypic variation in yeast

The power of yeast genetics has now been extensively applied to phenotypic variation among strains of Saccharomyces cerevisiae. As a result, over 100 genes and numerous sequence variants have been identified, providing us with a general characterization of mutations underlying quantitative trait variation. Most quantitative trait alleles exert considerable phenotypic effects and alter conserved amino acid positions within protein coding sequences. When examined, quantitative trait alleles influence the expression of numerous genes, most of which are unrelated to an allele's phenotypic effect. The profile of quantitative trait alleles has proven useful to reverse quantitative genetics approaches and supports the use of systems genetics approaches to synthesize the molecular basis of trait variation across multiple strains.

A noncomplementation screen for quantitative trait alleles in saccharomyces cerevisiae

Both linkage and linkage disequilibrium mapping provide well-defined approaches to mapping quantitative trait alleles. However, alleles of small effect are particularly difficult to refine to individual genes and causative mutations. Quantitative noncomplementation provides a means of directly testing individual genes for quantitative trait alleles in a fixed genetic background. Here, we implement a genome-wide noncomplementation screen for quantitative trait alleles that affect colony color or size by using the yeast deletion collection. As proof of principle, we find a previously known allele of CYS4 that affects colony color and a novel allele of CTT1 that affects resistance to hydrogen peroxide. To screen nearly 4700 genes in nine diverse yeast strains, we developed a high-throughput robotic plating assay to quantify colony color and size. Although we found hundreds of candidate alleles, reciprocal hemizygosity analysis of a select subset revealed that many of the candidates were false positives, in part the result of background-dependent haploinsufficiency or second-site mutations within the yeast deletion collection. Our results highlight the difficulty of identifying small-effect alleles but support the use of noncomplementation as a rapid means of identifying quantitative trait alleles of large effect.

Divergence of the yeast transcription factor FZF1 affects sulfite resistance

Changes in gene expression are commonly observed during evolution. However, the phenotypic consequences of expression divergence are frequently unknown and difficult to measure. Transcriptional regulators provide a mechanism by which phenotypic divergence can occur through multiple, coordinated changes in gene expression during development or in response to environmental changes. Yet, some changes in transcriptional regulators may be constrained by their pleiotropic effects on gene expression. Here, we use a genome-wide screen for promoters that are likely to have diverged in function and identify a yeast transcription factor, FZF1, that has evolved substantial differences in its ability to confer resistance to sulfites. Chimeric alleles from four Saccharomyces species show that divergence in FZF1 activity is due to changes in both its coding and upstream noncoding sequence. Between the two closest species, noncoding changes affect the expression of FZF1, whereas coding changes affect the expression of SSU1, a sulfite efflux pump activated by FZF1. Both coding and noncoding changes also affect the expression of many other genes. Our results show how divergence in the coding and promoter region of a transcription factor alters the response to an environmental stress.

Changes in gene regulation are thought to play an important role in evolution. While variation in gene expression between species is common, it is hard to identify the phenotypic consequences of this variation since many changes in gene expression may have subtle or no phenotypic effects. In this study, we investigate changes in sulfite resistance and gene expression caused by the transcription factor, FZF1, that has evolved rapidly during the divergence of related yeast species. We find that divergence in the ability of FZF1 to confer sulfite resistance is mediated by changes in its expression as well as changes in its protein structure, both of which cause changes in the expression of other genes. Our results show how the combination of multiple changes within a transcription factor can produce substantial changes in phenotype and the expression of many genes.

A model of developmental evolution: selection, pleiotropy and compensation

Development and physiology translate genetic variation into phenotypic variation and determine the genotype-phenotype map, such as which gene affects which character (pleiotropy). Any genetic change in this mapping reflects a change in development. Here, we discuss evidence for variation in pleiotropy and propose the selection, pleiotropy and compensation model (SPC) for adaptive evolution. It predicts that adaptive change in one character is associated with deleterious pleiotropy in others and subsequent selection to compensate for these pleiotropic effects. The SPC model provides a unifying perspective for a variety of puzzling phenomena, including developmental systems drift and character homogenization. The model suggests that most adaptive signatures detected in genome scans could be the result of compensatory changes, rather than of progressive character adaptations.

Phenotypic landscape of Saccharomyces cerevisiae during wine fermentation: evidence for origin-dependent metabolic traits

The species Saccharomyces cerevisiae includes natural strains, clinical isolates, and a large number of strains used in human activities. The aim of this work was to investigate how the adaptation to a broad range of ecological niches may have selectively shaped the yeast metabolic network to generate specific phenotypes. Using 72 S. cerevisiae strains collected from various sources, we provide, for the first time, a population-scale picture of the fermentative metabolic traits found in the S. cerevisiae species under wine making conditions. Considerable phenotypic variation was found suggesting that this yeast employs diverse metabolic strategies to face environmental constraints. Several groups of strains can be distinguished from the entire population on the basis of specific traits. Strains accustomed to growing in the presence of high sugar concentrations, such as wine yeasts and strains obtained from fruits, were able to achieve fermentation, whereas natural yeasts isolated from “poor-sugar” environments, such as oak trees or plants, were not. Commercial wine yeasts clearly appeared as a subset of vineyard isolates, and were mainly differentiated by their fermentative performances as well as their low acetate production. Overall, the emergence of the origin-dependent properties of the strains provides evidence for a phenotypic evolution driven by environmental constraints and/or human selection within S. cerevisiae.

Evolution of adaptive phenotypic variation patterns by direct selection for evolvability

A basic assumption of the Darwinian theory of evolution is that heritable variation arises randomly. In this context, randomness means that mutations arise irrespective of the current adaptive needs imposed by the environment. It is broadly accepted, however, that phenotypic variation is not uniformly distributed among phenotypic traits, some traits tend to covary, while others vary independently, and again others barely vary at all. Furthermore, it is well established that patterns of trait variation differ among species. Specifically, traits that serve different functions tend to be less correlated, as for instance forelimbs and hind limbs in bats and humans, compared with the limbs of quadrupedal mammals. Recently, a novel class of genetic elements has been identified in mouse gene-mapping studies that modify correlations among quantitative traits. These loci are called relationship loci, or relationship Quantitative Trait Loci (rQTL), and affect trait correlations by changing the expression of the existing genetic variation through gene interaction. Here, we present a population genetic model of how natural selection acts on rQTL. Contrary to the usual neo-Darwinian theory, in this model, new heritable phenotypic variation is produced along the selected dimension in response to directional selection. The results predict that selection on rQTL leads to higher correlations among traits that are simultaneously under directional selection. On the other hand, traits that are not simultaneously under directional selection are predicted to evolve lower correlations. These results and the previously demonstrated existence of rQTL variation, show a mechanism by which natural selection can directly enhance the evolvability of complex organisms along lines of adaptive change.

Evolution of pleiotropy: epistatic interaction pattern supports a mechanistic model underlying variation in genotype-phenotype map

The genotype-phenotype (GP) map consists of developmental and physiological mechanisms mapping genetic onto phenotypic variation. It determines the distribution of heritable phenotypic variance on which selection can act. Comparative studies of morphology as well as of gene regulatory networks show that the GP map itself evolves, yet little is known about the actual evolutionary mechanisms involved. The study of such mechanisms requires exploring the variation in GP maps at the population level, which presently is easier to quantify by statistical genetic methods rather than by regulatory network structures. We focus on the evolution of pleiotropy, a major structural aspect of the GP map. Pleiotropic genes affect multiple traits and underlie genetic covariance between traits, often causing evolutionary constraints. Previous quantitative genetic studies have demonstrated population-level variation in pleiotropy in the form of loci, at which genotypes differ in the genetic covariation between traits. This variation can potentially fuel evolution of the GP map under selection and/or drift. Here, we propose a developmental mechanism underlying population genetic variation in covariance and test its predictions. Specifically, the mechanism predicts that the loci identified as responsible for genetic variation in pleiotropy are involved in trait-specific epistatic interactions. We test this prediction for loci affecting allometric relationships between traits in an advanced intercross between inbred mouse strains. The results consistently support the prediction. We further find a high degree of sign epistasis in these interactions, which we interpret as an indication of adaptive gene complexes within the diverged parental lines.

Sulfate assimilation mediates tellurite reduction and toxicity in Saccharomyces cerevisiae

Despite a century of research and increasing environmental and human health concerns, the mechanistic basis of the toxicity of derivatives of the metalloid tellurium, Te, in particular the oxyanion tellurite, Te(IV), remains unsolved. Here, we provide an unbiased view of the mechanisms of tellurium metabolism in the yeast Saccharomyces cerevisiae by measuring deviations in Te-related traits of a complete collection of gene knockout mutants. Reduction of Te(IV) and intracellular accumulation as metallic tellurium strongly correlated with loss of cellular fitness, suggesting that Te(IV) reduction and toxicity are causally linked. The sulfate assimilation pathway upstream of Met17, in particular, the sulfite reductase and its cofactor siroheme, was shown to be central to tellurite toxicity and its reduction to elemental tellurium. Gene knockout mutants with altered Te(IV) tolerance also showed a similar deviation in tolerance to both selenite and, interestingly, selenomethionine, suggesting that the toxicity of these agents stems from a common mechanism. We also show that Te(IV) reduction and toxicity in yeast is partially mediated via a mitochondrial respiratory mechanism that does not encompass the generation of substantial oxidative stress. The results reported here represent a robust base from which to attack the mechanistic details of Te(IV) toxicity and reduction in a eukaryotic organism.

A combined-cross analysis reveals genes with drug-specific and background-dependent effects on drug sensitivity in Saccharomyces cerevisiae

Effective pharmacological therapy is often inhibited by variable drug responses and adverse drug reactions. Dissecting the molecular basis of different drug responses is difficult due to complex interactions involving multiple genes, pathways, and cellular processes. We previously found a single nucleotide polymorphism within cystathionine beta-synthase (CYS4) that causes multi-drug sensitivity in a vineyard strain of Saccharomyces cerevisiae. However, not all variation was accounted for by CYS4. To identify additional genes influencing drug sensitivity, we used CYS4 as a covariate and conducted both single- and combined-cross linkage mapping. After eliminating numerous false-positive associations, we identified 16 drug-sensitivity loci, only 3 of which had been previously identified. Of 4 drug-sensitivity loci selected for validation, 2 showed replicated associations in independent crosses, and two quantitative trait genes within these regions, AQY1 and MKT1, were found to have drug-specific and background-dependent effects. Our results suggest that drug response may often depend on interactions between genes with multi-drug and drug-specific effects.