Cascading transcriptional effects of a naturally occurring frameshift mutation in Saccharomyces cerevisiae

Using wasps as a tool to restore a functioning vine grape mycobiota and preserve the mycobial "terroir"

In the last one-hundred years, the exponential expansion of wine making has artificialized the agricultural landscape as well as its microbial diversity, spreading human selected Saccharomyces cerevisiae strains. Evidence showed that social wasps can harbor a significant fraction of the yeast phenotypic diversity of a given area of wine production, allowing different strains to overwinter and mate in their gut. The integrity of the wasp-yeast ecological interaction is of paramount importance to maintain the resilience of microbial populations associated to wine aromatic profiles. In a field experiment, we verified whether Polistes dominula wasps, reared in laboratory and fed with a traceable S. cerevisiae strain, could be a useful tool to drive the controlled yeast dispersion directly on grapes. The demonstration of the biotechnological potential of social insects in organic wine farming lays the foundations for multiple applications including maintenance of microbial biodiversity and rewilding vineyards through the introduction of wasp associated microbiomes.

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.

Molecular and evolutionary processes generating variation in gene expression

Heritable variation in gene expression is common within and between species. This variation arises from mutations that alter the form or function of molecular gene regulatory networks that are then filtered by natural selection. High-throughput methods for introducing mutations and characterizing their cis- and trans-regulatory effects on gene expression (particularly, transcription) are revealing how different molecular mechanisms generate regulatory variation, and studies comparing these mutational effects with variation seen in the wild are teasing apart the role of neutral and non-neutral evolutionary processes. This integration of molecular and evolutionary biology allows us to understand how the variation in gene expression we see today came to be and to predict how it is most likely to evolve in the future.

Mechanism of Long-Range Chromosome Motion Triggered by Gene Activation

Movement of chromosome sites within interphase cells is critical for numerous pathways including RNA transcription and genome organization. Yet, a mechanism for reorganizing chromatin in response to these events had not been reported. Here, we delineate a molecular chaperone-dependent pathway for relocating activated gene loci in yeast. Our presented data support a model in which a two-authentication system mobilizes a gene promoter through a dynamic network of polymeric nuclear actin. Transcription factor-dependent nucleation of a myosin motor propels the gene locus through the actin matrix, and fidelity of the actin association was ensured by ARP-containing chromatin remodelers. Motor activity of nuclear myosin was dependent on the Hsp90 chaperone. Hsp90 further contributed by biasing the remodeler-actin interaction toward nucleosomes with the non-canonical histone H2A.Z, thereby focusing the pathway on select sites such as transcriptionally active genes. Together, the system provides a rapid and effective means to broadly yet selectively mobilize chromatin sites.

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.

Genetics of trans-regulatory variation in gene expression

Heritable variation in gene expression forms a crucial bridge between genomic variation and the biology of many traits. However, most expression quantitative trait loci (eQTLs) remain unidentified. We mapped eQTLs by transcriptome sequencing in 1012 yeast segregants. The resulting eQTLs accounted for over 70% of the heritability of mRNA levels, allowing comprehensive dissection of regulatory variation. Most genes had multiple eQTLs. Most expression variation arose from trans-acting eQTLs distant from their target genes. Nearly all trans-eQTLs clustered at 102 hotspot locations, some of which influenced the expression of thousands of genes. Fine-mapped hotspot regions were enriched for transcription factor genes. While most genes had a local eQTL, most of these had no detectable effects on the expression of other genes in trans. Hundreds of non-additive genetic interactions accounted for small fractions of expression variation. These results reveal the complexity of genetic influences on transcriptome variation in unprecedented depth and detail.

Non-invasive real time monitoring of yeast volatilome by PTR-ToF-MS

INTRODUCTION

Producing a wide range of volatile secondary metabolites Saccharomyces cerevisiae influences wine, beer, and bread sensory quality and hence selection of strains based on their volatilome becomes pivotal. A rapid on-line method for volatilome assessing of strains growing on standard solid media is still missing.

OBJECTIVES

Methodologically, the aim of this study was to demonstrate the automatic, real-time, direct, and non-invasive monitoring of yeast volatilome in order to rapidly produce a robust large data set encompassing measurements relative to many strains, replicates and time points. The fundamental scope was to differentiate volatilomes of genetically similar strains of oenological relevance during the whole growing process.

METHOD

Six different S. cerevisiae strains (four meiotic segregants of a natural strain and two laboratory strains) inoculated onto a solid medium have been monitored on-line by Proton Transfer Reaction-Time-of-Flight-Mass Spectrometry for 11 days every 4 h (3540 time points). FastGC PTR-ToF-MS was performed during the stationary phase on the 5th day.

RESULTS

More than 300 peaks have been extracted from the average spectra associated to each time point, 70 have been tentatively identified. Univariate and multivariate analyses have been performed on the data matrix (3640 measurements × 70 peaks) highlighting the volatilome evolution and strain-specific features. Laboratory strains with opposite mating type, and meiotic segregants of the same natural strain showed significantly different profiles.

CONCLUSIONS

The described set-up allows the on-line high-throughput screening of yeast volatilome of S. cerevisiae strains and the identification of strain specific features and new metabolic pathways, discriminating also genetically similar strains, thus revealing a novel method for strain phenotyping, identification, and quality control.

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.

Gene Expression Evolves under a House-of-Cards Model of Stabilizing Selection

Divergence in gene regulation is hypothesized to underlie much of phenotypic evolution, but the role of natural selection in shaping the molecular phenotype of gene expression continues to be debated. To resolve the mode of gene expression, evolution requires accessible theoretical predictions for the effect of selection over long timescales. Evolutionary quantitative genetic models of phenotypic evolution can provide such predictions, yet those predictions depend on the underlying hypotheses about the distributions of mutational and selective effects that are notoriously difficult to disentangle. Here, we draw on diverse genomic data sets including expression profiles of natural genetic variation and mutation accumulation lines, empirical estimates of genomic mutation rates, and inferences of genetic architecture to differentiate contrasting hypotheses for the roles of stabilizing selection and mutation in shaping natural expression variation. Our analysis suggests that gene expression evolves in a domain of phenotype space well fit by the House-of-Cards (HC) model. Although the strength of selection inferred is sensitive to the number of loci controlling gene expression, the model is not. The consistency of these results across evolutionary time from budding yeast through fruit fly implies that this model is general and that mutational effects on gene expression are relatively large. Empirical estimates of the genetic architecture of gene expression traits imply that selection provides modest constraints on gene expression levels for most genes, but that the potential for regulatory evolution is high. Our prediction using data from laboratory environments should encourage the collection of additional data sets allowing for more nuanced parameterizations of HC models for gene expression.

Unravelling glutathione conjugate catabolism in Saccharomyces cerevisiae: the role of glutathione/dipeptide transporters and vacuolar function in the release of volatile sulfur compounds 3-mercaptohexan-1-ol and 4-mercapto-4-methylpentan-2-one

Sulfur-containing aroma compounds are key contributors to the flavour of a diverse range of foods and beverages, such as wine. The tropical fruit characters of Sauvignon Blanc wines are attributed to the presence of the aromatic thiols 3-mercaptohexan-1-ol (3-MH), its acetate ester 3-mercaptohexyl acetate (3-MHA), and 4-mercapto-4-methylpentan-2-one (4-MMP). These aromatic thiols are not detectable in grape juice to any significant extent but are released by yeast during alcoholic fermentation. While the processes involved in the release of 3-MH and 4-MMP from their cysteinylated precursors have been studied extensively, degradation pathways for glutathione S-conjugates (GSH-3-MH and GSH-4-MMP) have not. In this study, a candidate gene approach was taken, focusing on genes known to play a role in glutathione and glutathione-S-conjugate turnover in Saccharomyces cerevisiae. Our results confirm the role of Opt1p as the major transporter responsible for uptake of GSH-3-MH and GSH-4-MMP, and identify vacuolar Ecm38p as a key determinant of 3-MH release from GSH-3-MH. ECM38 was unimportant, on the other hand, for release of 4-MMP, and abolition of vacuolar biogenesis caused an increase in the amount of 4-MMP released. The alternative cytosolic glutathione degradation pathway was not involved in release of either thiol from their glutathionylated precursors. Finally, cycling of GSH-3-MH and/or its breakdown intermediates between the cytosol and the vacuole or extracellular space was implicated in modulation of 3-MH formation. Together, these results provide new targets for development of yeast strains that optimize release of these potent volatile sulfur compounds, and further our understanding of the processes involved in glutathione-S-conjugate turnover.

Reduced Ssy1-Ptr3-Ssy5 (SPS) signaling extends replicative life span by enhancing NAD+ homeostasis in Saccharomyces cerevisiae

Attenuated nutrient signaling extends the life span in yeast and higher eukaryotes; however, the mechanisms are not completely understood. Here we identify the Ssy1-Ptr3-Ssy5 (SPS) amino acid sensing pathway as a novel longevity factor. A null mutation of SSY5 (ssy5Δ) increases replicative life span (RLS) by ∼50%. Our results demonstrate that several NAD(+) homeostasis factors play key roles in this life span extension. First, expression of the putative malate-pyruvate NADH shuttle increases in ssy5Δ cells, and deleting components of this shuttle, MAE1 and OAC1, largely abolishes RLS extension. Next, we show that Stp1, a transcription factor of the SPS pathway, directly binds to the promoter of MAE1 and OAC1 to regulate their expression. Additionally, deletion of SSY5 increases nicotinamide riboside (NR) levels and phosphate-responsive (PHO) signaling activity, suggesting that ssy5Δ increases NR salvaging. This increase contributes to NAD(+) homeostasis, partially ameliorating the NAD(+) deficiency and rescuing the short life span of the npt1Δ mutant. Moreover, we observed that vacuolar phosphatase, Pho8, is partially required for ssy5Δ-mediated NR increase and RLS extension. Together, our studies present evidence that supports SPS signaling is a novel NAD(+) homeostasis factor and ssy5Δ-mediated life span extension is likely due to concomitantly increased mitochondrial and vacuolar function. Our findings may contribute to understanding the molecular basis of NAD(+) metabolism, cellular life span, and diseases associated with NAD(+) deficiency and aging.

Experimental selection of long-term intracellular mycobacteria

Some intracellular bacteria are known to cause long-term infections that last decades without compromising the viability of the host. Although of critical importance, the adaptations that intracellular bacteria undergo during this long process of residence in a host cell environment remain obscure. Here, we report a novel experimental approach to study the adaptations of mycobacteria imposed by a long-term intracellular lifestyle. Selected Mycobacterium bovis BCG through continuous culture in macrophages underwent an adaptation process leading to impaired phenolic glycolipids (PGL) synthesis, improved usage of glucose as a carbon source and accumulation of neutral lipids. These changes correlated with increased survival of mycobacteria in macrophages and mice during re-infection and also with the specific expression of stress- and survival-related genes. Our findings identify bacterial traits implicated in the establishment of long-term cellular infections and represent a tool for understanding the physiological states and the environment that bacteria face living in fluctuating intracellular environments.

Proteome variance differences within populations of European whitefish (Coregonus lavaretus) originating from contrasting salinity environments

Variation in gene expression is an important component of the phenotypic differences observed in nature. Gene expression variance across biological groups and environmental conditions has been studied extensively and has revealed specific genes and molecular mechanisms of interest. However, little is known regarding the importance of within-population gene expression variation to environmental adaptation. To address this issue, we quantified the proteomes of individuals of European whitefish (Coregonus lavaretus) from populations that have previously been shown to have adapted during early development to freshwater and brackishwater salinity environments. Using MS-based label-free proteomics, we studied 955 proteins in eight hatch-stage fish embryos from each population that had been reared in either freshwater or brackishwater salinity conditions. By comparing the levels of within-population protein expression variance over individuals and per protein between populations, we found that fish embryos from the population less affected by salinity level had also markedly higher levels of expression variance. Gene Ontologies and molecular pathways associated with osmoregulation showed the most significant difference of within-population proteome variance between populations. Several new candidate genes for salinity adaptation were identified, emphasising the added value of combining assessments of within-population gene expression variation with standard gene expression analysis practices for better understanding the mechanisms of environmental adaptation.

BIOLOGICAL SIGNIFICANCE

We demonstrate the benefits of studying within-population gene expression variance together with more typical methods of gene expression profiling. Proteome variance differences within European whitefish populations originating from different salinity environments allowed us to identify several new candidate genes for salinity adaptation in teleost fish and generate many further hypotheses to be tested. This article is part of a Special Issue entitled: Proteomics of non-model organisms.

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.

Differential adaptation to multi-stressed conditions of wine fermentation revealed by variations in yeast regulatory networks

BACKGROUND

Variation of gene expression can lead to phenotypic variation and have therefore been assumed to contribute the diversity of wine yeast (Saccharomyces cerevisiae) properties. However, the molecular bases of this variation of gene expression are unknown. We addressed these questions by carrying out an integrated genetical-genomic study in fermentation conditions. We report here quantitative trait loci (QTL) mapping based on expression profiling in a segregating population generated by a cross between a derivative of the popular wine strain EC1118 and the laboratory strain S288c.

RESULTS

Most of the fermentation traits studied appeared to be under multi-allelic control. We mapped five phenotypic QTLs and 1465 expression QTLs. Several expression QTLs overlapped in hotspots. Among the linkages unraveled here, several were associated with metabolic processes essential for wine fermentation such as glucose sensing or nitrogen and vitamin metabolism. Variations affecting the regulation of drug detoxification and export (TPO1, PDR12 or QDR2) were linked to variation in four genes encoding transcription factors (PDR8, WAR1, YRR1 and HAP1). We demonstrated that the allelic variation of WAR1 and TPO1 affected sorbic and octanoic acid resistance, respectively. Moreover, analysis of the transcription factors phylogeny suggests they evolved with a specific adaptation of the strains to wine fermentation conditions. Unexpectedly, we found that the variation of fermentation rates was associated with a partial disomy of chromosome 16. This disomy resulted from the well known 8-16 translocation.

CONCLUSIONS

This large data set made it possible to decipher the effects of genetic variation on gene expression during fermentation and certain wine fermentation properties. Our findings shed a new light on the adaptation mechanisms required by yeast to cope with the multiple stresses generated by wine fermentation. In this context, the detoxification and export systems appear to be of particular importance, probably due to nitrogen starvation. Furthermore, we show that the well characterized 8-16 translocation located in SSU1, which is associated with sulfite resistance, can lead to a partial chromosomic amplification in the progeny of strains that carry it, greatly improving fermentation kinetics. This amplification has been detected among other wine yeasts.

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 computational pipeline to discover highly phylogenetically informative genes in sequenced genomes: application to Saccharomyces cerevisiae natural strains

The quest for genes representing genetic relationships of strains or individuals within populations and their evolutionary history is acquiring a novel dimension of complexity with the advancement of next-generation sequencing (NGS) technologies. In fact, sequencing an entire genome uncovers genetic variation in coding and non-coding regions and offers the possibility of studying Saccharomyces cerevisiae populations at the strain level. Nevertheless, the disadvantageous cost-benefit ratio (the amount of details disclosed by NGS against the time-expensive and expertise-demanding data assembly process) still precludes the application of these techniques to the routinely assignment of yeast strains, making the selection of the most reliable molecular markers greatly desirable. In this work we propose an original computational approach to discover genes that can be used as a descriptor of the population structure. We found 13 genes whose variability can be used to recapitulate the phylogeny obtained from genome-wide sequences. The same approach that we prove to be successful in yeasts can be generalized to any other population of individuals given the availability of high-quality genomic sequences and of a clear population structure to be targeted.

The genotype-phenotype maps of systems biology and quantitative genetics: distinct and complementary

The processes by which genetic variation in complex traits is generated and maintained in populations has for a long time been treated in abstract and statistical terms. As a consequence, quantitative genetics has provided limited insights into our understanding of the molecular bases of quantitative trait variation. With the developing technological and conceptual tools of systems biology, cellular and molecular processes are being described in greater detail. While we have a good description of how signaling and other molecular networks are organized in the cell, we still do not know how genetic variation affects these pathways, because systems and molecular biology usually ignore the type and extent of genetic variation found in natural populations. Here we discuss the quantitative genetics and systems biology approaches for the study of complex trait architecture and discuss why these two disciplines would synergize with each other to answer questions that neither of the two could answer alone.

Deciphering the molecular basis of wine yeast fermentation traits using a combined genetic and genomic approach

The genetic basis of the phenotypic diversity of yeast is still poorly understood. Wine yeast strains have specific abilities to grow and ferment under stressful conditions compared with other strains, but the genetic basis underlying these traits is unknown. Understanding how sequence variation influences such phenotypes is a major challenge to address adaptation mechanisms of wine yeast. We aimed to identify the genetic basis of fermentation traits and gain insight into their relationships with variations in gene expression among yeast strains. We combined fermentation trait QTL mapping and expression profiling of fermenting cells in a segregating population from a cross between a wine yeast derivative and a laboratory strain. We report the identification of QTL for various fermentation traits (fermentation rates, nitrogen utilization, metabolites production) as well as expression QTL (eQTL). We found that many transcripts mapped to several eQTL hotspots and that two of them overlapped with QTL for fermentation traits. A QTL controlling the maximal fermentation rate and nitrogen utilization overlapping with an eQTL hotspot was dissected. We functionally demonstrated that an allele of the ABZ1 gene, localized in the hotspot and involved in p-aminobenzoate biosynthesis, controls the fermentation rate through modulation of nitrogen utilization. Our data suggest that the laboratory strain harbors a defective ABZ1 allele, which triggers strong metabolic and physiological alterations responsible for the generation of the eQTL hotspot. They also suggest that a number of gene expression differences result from some alleles that trigger major physiological disturbances.

Expression variability of co-regulated genes differentiates Saccharomyces cerevisiae strains

BACKGROUND

Saccharomyces cerevisiae (Baker's yeast) is found in diverse ecological niches and is characterized by high adaptive potential under challenging environments. In spite of recent advances on the study of yeast genome diversity, little is known about the underlying gene expression plasticity. In order to shed new light onto this biological question, we have compared transcriptome profiles of five environmental isolates, clinical and laboratorial strains at different time points of fermentation in synthetic must medium, during exponential and stationary growth phases.

RESULTS

Our data unveiled diversity in both intensity and timing of gene expression. Genes involved in glucose metabolism and in the stress response elicited during fermentation were among the most variable. This gene expression diversity increased at the onset of stationary phase (diauxic shift). Environmental isolates showed lower average transcript abundance of genes involved in the stress response, assimilation of nitrogen and vitamins, and sulphur metabolism, than other strains. Nitrogen metabolism genes showed significant variation in expression among the environmental isolates.

CONCLUSIONS

Wild type yeast strains respond differentially to the stress imposed by nutrient depletion, ethanol accumulation and cell density increase, during fermentation of glucose in synthetic must medium. Our results support previous data showing that gene expression variability is a source of phenotypic diversity among closely related organisms.

Gene network architecture as a canvas for the interpretation of ecological genomics investigations

New technologies promise to revolutionize the field of molecular ecology. This technological progress comes with its own set of challenges. Among the most important ones is the analysis and interpretation of the data in a way that tells us about the molecular causes of the phenotype of interest and its consequences. In this issue, Whitehead et al. (2010) reveal part of the mechanistic basis of evolved pollution tolerance by studying the developmental and transcriptional response of tolerant and sensitive fish embryos to polychlorinated biphenyls (PCBs), a pollutant commonly found in coastal waters of the United States. By integrating their gene expression profiling data with phenotypic data on individuals along with what is known about pathways by which this pollutant acts in zebrafish and mammals, they are able to suggest detailed mechanisms that have evolved to allow a fish population to adapt to a very damaging pollutant and develop normally.

Evolution of transcriptional regulatory networks in yeast populations

Saccharomyces cerevisiae is the most thoroughly studied eukaryote at the cellular, molecular, and genetic level. Recent boost in whole-genome sequencing, array-based allelic variation mapping, and genome-wide transcriptional profiling have unprecedentedly advanced knowledge on cell biology and evolution of this organism. It is now possible to investigate how evolution shapes the functional architecture of yeast genomes and how this architecture relates to the evolution of the regulatory networks controlling the expression of genes that make up an organism. A survey of the information on genetic and whole-genome expression variations in yeast populations shows that a significant score of gene expression variation is dependent on genotype-by-environment interaction. In some cases, large trans effects are the result of mutations in the promoters of key master regulator genes. Yet trans-variation in environmental sensor proteins appears to explain the majority of the expression patterns differentiating strains in natural populations. The challenge is now to use this information to model how individual genetic polymorphisms interact in a condition-dependent fashion to produce phenotypic change. In this study, we show how fruitful application of systems biology to the progress of science and medicine requires the use of evolution as a lens to reconstruct the hierarchical structure of regulation of biological systems. The lessons learned in yeast can be of paramount importance in advancing the application of genomics and systems biology to emerging fields including personalized medicine.

Natural variation in gene expression in the early development of dauer larvae of Caenorhabditis elegans

BACKGROUND

The free-living nematode Caenorhabditis elegans makes a developmental decision based on environmental conditions: larvae either arrest as dauer larva, or continue development into reproductive adults. There is natural variation among C. elegans lines in the sensitivity of this decision to environmental conditions; that is, there is variation in the phenotypic plasticity of dauer larva development. We hypothesised that these differences may be transcriptionally controlled in early stage larvae. We investigated this by microarray analysis of different C. elegans lines under different environmental conditions, specifically the presence and absence of dauer larva-inducing pheromone.

RESULTS

There were substantial transcriptional differences between four C. elegans lines under the same environmental conditions. The expression of approximately 2,000 genes differed between genetically different lines, with each line showing a largely line-specific transcriptional profile. The expression of genes that are markers of larval moulting suggested that the lines may be developing at different rates. The expression of a total of 89 genes was putatively affected by dauer larva or non-dauer larva-inducing conditions. Among the upstream regions of these genes there was an over-representation of DAF-16-binding motifs.

CONCLUSION

Under the same environmental conditions genetically different lines of C. elegans had substantial transcriptional differences. This variation may be due to differences in the developmental rates of the lines. Different environmental conditions had a rather smaller effect on transcription. The preponderance of DAF-16-binding motifs upstream of these genes was consistent with these genes playing a key role in the decision between development into dauer or into non-dauer larvae. There was little overlap between the genes whose expression was affected by environmental conditions and previously identified loci involved in the plasticity of dauer larva development.

Temporal constraints on the incorporation of regulatory mutants in evolutionary pathways

Understanding the molecular details of the sequence of events in multistep evolutionary pathways can reveal the extent to which natural selection exploits regulatory mutations affecting expression, amino acid replacements affecting the active site, amino acid replacements affecting protein folding or stability, or variations affecting gene copy number. In experimentally exploring the adaptive landscape of the evolution of resistance to beta-lactam antibiotics in enteric bacteria, we noted that a regulatory mutation that increases beta-lactamase expression by about 2-fold has a very strong tendency to be fixed at or near the end of the evolutionary pathway. This pattern contrasts with previous experiments selecting for the utilization of novel substrates, in which regulatory mutations that increase expression are often fixed early in the process. To understand the basis of the difference, we carried out experiments in which the expression of beta-lactamase was under the control of a tunable arabinose promoter. We find that the fitness effect of an increase in gene expression is highly dependent on the catalytic activity of the coding sequence. An increase in expression of an inefficient enzyme has a negligible effect on drug resistance; however, the effect of an increase in expression of an efficient enzyme is very large. The contrast in the temporal incorporation of regulatory mutants between antibiotic resistance and the utilization of novel substrates is related to the nature of the function that relates enzyme activity to fitness. A mathematical model of beta-lactam resistance is examined in detail and shown to be consistent with the observed results.

Dissecting the pleiotropic consequences of a quantitative trait nucleotide

The downstream consequences of a single quantitative trait polymorphism can provide important insight into the molecular basis of a trait. However, the molecular consequences of a polymorphism may be complex and only a subset of these may influence the trait of interest. In natural isolates of Saccharomyces cerevisiae, a nonsynonymous polymorphism in cystathione beta-synthase (CYS4) causes a deficiency in both cysteine and glutathione that results in rust-colored colonies and drug-dependent growth defects. Using a single-nucleotide allele replacement, we characterized the effects of this polymorphism on gene expression levels across the genome. To determine whether any of the differentially expressed genes are necessary for the production of rust-colored colonies, we screened the yeast deletion collection for genes that enhance or suppress rust coloration. We found that genes in the sulfur assimilation pathway are required for the production of rust color but not the drug-sensitivity phenotype. Our results show that a single quantitative trait polymorphism can generate a complex set of downstream changes, providing a molecular basis for pleiotropy.

Variations in stress sensitivity and genomic expression in diverse S. cerevisiae isolates

Interactions between an organism and its environment can significantly influence phenotypic evolution. A first step toward understanding this process is to characterize phenotypic diversity within and between populations. We explored the phenotypic variation in stress sensitivity and genomic expression in a large panel of Saccharomyces strains collected from diverse environments. We measured the sensitivity of 52 strains to 14 environmental conditions, compared genomic expression in 18 strains, and identified gene copy-number variations in six of these isolates. Our results demonstrate a large degree of phenotypic variation in stress sensitivity and gene expression. Analysis of these datasets reveals relationships between strains from similar niches, suggests common and unique features of yeast habitats, and implicates genes whose variable expression is linked to stress resistance. Using a simple metric to suggest cases of selection, we found that strains collected from oak exudates are phenotypically more similar than expected based on their genetic diversity, while sake and vineyard isolates display more diverse phenotypes than expected under a neutral model. We also show that the laboratory strain S288c is phenotypically distinct from all of the other strains studied here, in terms of stress sensitivity, gene expression, Ty copy number, mitochondrial content, and gene-dosage control. These results highlight the value of understanding the genetic basis of phenotypic variation and raise caution about using laboratory strains for comparative genomics.

Much attention has been given to the ways in which organisms evolve new phenotypes and the influence of the environment on this process. A major focus of study is defining the genetic basis for phenotypes important for organismal fitness. As a first step toward this goal, we surveyed phenotypic variation in diverse yeast strains collected from different environments by characterizing variations in stress resistance and genomic expression. We uncovered many phenotypic differences across yeast strains, both in stress tolerance and gene expression. The similarities and differences of the strains analyzed uncovered phenotypes shared by strains that live in similar environments, suggesting common features of yeast niches as well as mechanisms that different strains use to thrive in those conditions. We provide evidence that some characteristics of strains isolated from oak tree soil have been selected for, perhaps because of the shared selective pressures imposed by their environment. One theme emerging from our studies is that the laboratory strain of yeast, long used as a model for yeast physiology and basic biology, is aberrant compared to all other strains. This result raises caution about making general conclusions about yeast biology based on a single strain with a specific genetic makeup.

Dominance and the evolutionary accumulation of cis- and trans-effects on gene expression

Gene expression levels appear to be under pervasive stabilizing selection. Yet the genetic architecture underlying abundant gene expression diversity within and between populations remains elusive. Here, we investigated the role of dominance in the segregation of cis- and trans-regulation within and between populations. We used chromosome substitution lines of Drosophila melanogaster to show that (i) >70% of the genes that are differentially expressed between two homozygous lines are masked in the heterozygous, suggesting that one of the substituted chromosomes contains a recessive allele; (ii) such large masking is already obtained with heterozygous chromosomes originating from the same population, with the time of divergence between chromosomes in heterozygous lines making only a small but significant contribution to the masking of variation observed in homozygous lines; (iii) variation in gene expression due to trans-regulation is biased toward greater deviations from additivity because of recessive and dominant alleles, whereas variation due to cis-regulation shows higher additivity; and (iv) genetic divergence between second chromosomes is associated with increased cis-regulation, whereas the level of trans-regulation shows little increase over the time scale studied. Our results indicate that cis-acting alleles may be preferentially fixed by positive natural selection because of their higher additivity, and that the disruption of gene expression by recessive variation with pervasive trans-effects may be important for understanding gene expression variation within populations. We suggest that widespread regulatory effects of recessive low-frequency homozygous variation may provide a general mechanism mediating disease phenotypes and the genetic load of natural populations.

Sleuthing the difference a nucleotide can make

Identifying the genetic basis of natural variation in genome-wide gene expression is a daunting task, yet central to characterizing the potential evolutionary impact of gene expression polymorphism. In this issue, Brown et al. reveal the genetic variation that underlies a significant portion of the gene expression polymorphism that segregates in a natural isolate of the wine and bread yeast Saccharomyces cerevisiae. That polymorphism is shown to be associated with a single nucleotide insertion. Interestingly, the major expression polymorphism is not caused by an insertion in a cis-regulatory region or in a transcription factor, but in a homonucleotide repeat within the coding sequence of SSY1, which encodes an amino acid sensor. The expression polymorphism, perhaps unsurprisingly, is also environment-dependent and leads to differential fitness across environments. The result calls for a deeper consideration within molecular ecological genomics of the gene-by-environment interactions that lead to differential gene expression. The finding also calls for further data, clarifying the kinds of genetic variation that constitute gene expression polymorphism between individuals in natural populations.