Elongation rate of RNA polymerase II affects pausing patterns across 3' UTRs.

Yeast mRNAs are polyadenylated at multiple sites in their 3' untranslated regions (3' UTRs), and poly(A) site usage is regulated by the rate of transcriptional elongation by RNA polymerase II (Pol II). Slow Pol II derivatives favor upstream poly(A) sites, and fast Pol II derivatives favor downstream poly(A) sites. Transcriptional elongation and polyadenylation are linked at the nucleotide level, presumably reflecting Pol II dwell time at each residue that influences the level of polyadenylation. Here, we investigate the effect of Pol II elongation rate on pausing patterns and the relationship between Pol II pause sites and poly(A) sites within 3' UTRs. Mutations that affect Pol II elongation rate alter sequence preferences at pause sites within 3' UTRs, and pausing preferences differ between 3' UTRs and coding regions. In addition, sequences immediately flanking the pause sites show preferences that are largely independent of Pol II speed. In wild-type cells, poly(A) sites are preferentially located < 50 nucleotides upstream from Pol II pause sites, but this spatial relationship is diminished in cells harboring Pol II speed mutants. Based on a random forest classifier, Pol II pause sites are modestly predicted by the distance to poly(A) sites but are better predicted by the chromatin landscape in Pol II speed derivatives. Transcriptional regulatory proteins can influence the relationship between Pol II pausing and polyadenylation but in a manner distinct from Pol II elongation rate derivatives. These results indicate a complex relationship between Pol II pausing and polyadenylation.

Transcription elongation is finely tuned by dozens of regulatory factors.

Understanding the complex network that regulates transcription elongation requires the quantitative analysis of RNA polymerase II (Pol II) activity in a wide variety of regulatory environments. We performed native elongating transcript sequencing (NET-seq) in 41 strains of Saccharomyces cerevisiae lacking known elongation regulators, including RNA processing factors, transcription elongation factors, chromatin modifiers, and remodelers. We found that the opposing effects of these factors balance transcription elongation and antisense transcription. Different sets of factors tightly regulate Pol II progression across gene bodies so that Pol II density peaks at key points of RNA processing. These regulators control where Pol II pauses with each obscuring large numbers of potential pause sites that are primarily determined by DNA sequence and shape. Antisense transcription varies highly across the regulatory landscapes analyzed, but antisense transcription in itself does not affect sense transcription at the same locus. Our findings collectively show that a diverse array of factors regulate transcription elongation by precisely balancing Pol II activity.

Simultaneous quantification of mRNA and protein in single cells reveals post-transcriptional effects of genetic variation.

Trans-acting DNA variants may specifically affect mRNA or protein levels of genes located throughout the genome. However, prior work compared trans-acting loci mapped in separate studies, many of which had limited statistical power. Here, we developed a CRISPR-based system for simultaneous quantification of mRNA and protein of a given gene via dual fluorescent reporters in single, live cells of the yeast Saccharomyces cerevisiae. In large populations of recombinant cells from a cross between two genetically divergent strains, we mapped 86 trans-acting loci affecting the expression of ten genes. Less than 20% of these loci had concordant effects on mRNA and protein of the same gene. Most loci influenced protein but not mRNA of a given gene. One locus harbored a premature stop variant in the YAK1 kinase gene that had specific effects on protein or mRNA of dozens of genes. These results demonstrate complex, post-transcriptional genetic effects on gene expression.

Pervasive Phenotypic Impact of a Large Nonrecombining Introgressed Region in Yeast.

To explore the origin of the diversity observed in natural populations, many studies have investigated the relationship between genotype and phenotype. In yeast species, especially in Saccharomyces cerevisiae, these studies are mainly conducted using recombinant offspring derived from two genetically diverse isolates, allowing to define the phenotypic effect of genetic variants. However, large genomic variants such as interspecies introgressions are usually overlooked even if they are known to modify the genotype-phenotype relationship. To have a better insight into the overall phenotypic impact of introgressions, we took advantage of the presence of a 1-Mb introgressed region, which lacks recombination and contains the mating-type determinant in the Lachancea kluyveri budding yeast. By performing linkage mapping analyses in this species, we identified a total of 89 loci affecting growth fitness in a large number of conditions and 2,187 loci affecting gene expression mostly grouped into two major hotspots, one being the introgressed region carrying the mating-type locus. Because of the absence of recombination, our results highlight the presence of a sexual dimorphism in a budding yeast for the first time. Overall, by describing the phenotype-genotype relationship in the Lachancea kluyveri species, we expanded our knowledge on how genetic characteristics of large introgression events can affect the phenotypic landscape.

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.

RNAi as a Tool to Study Virulence in the Pathogenic Yeast Candida glabrata.

The yeast Candida glabrata is a major opportunistic pathogen causing mucosal and systemic infections in humans. Systemic infections caused by this yeast have high mortality rates and are difficult to treat due to this yeast's intrinsic and frequently adapting antifungal resistance. To understand and treat C. glabrata infections, it is essential to investigate the molecular basis of C. glabrata virulence and resistance. We established an RNA interference (RNAi) system in C. glabrata by expressing the Dicer and Argonaute genes from Saccharomyces castellii (a budding yeast with natural RNAi). Our experiments with reporter genes and putative virulence genes showed that the introduction of RNAi resulted in 30 and 70% gene-knockdown for the construct-types antisense and hairpin, respectively. The resulting C. glabrata RNAi strain was used for the screening of a gene library for new virulence-related genes. Phenotypic profiling with a high-resolution quantification of growth identified genes involved in the maintenance of cell integrity, antifungal drugs, and ROS resistance. The genes identified by this approach are promising targets for the treatment of C. glabrata infections.

Variation of the meiotic recombination landscape and properties over a broad evolutionary distance in yeasts.

Meiotic recombination is a major factor of genome evolution, deeply characterized in only a few model species, notably the yeast Saccharomyces cerevisiae. Consequently, little is known about variations of its properties across species. In this respect, we explored the recombination landscape of Lachancea kluyveri, a protoploid yeast species that diverged from the Saccharomyces genus more than 100 million years ago and we found striking differences with S. cerevisiae. These variations include a lower recombination rate, a higher frequency of chromosomes segregating without any crossover and the absence of recombination on the chromosome arm containing the sex locus. In addition, although well conserved within the Saccharomyces clade, the S. cerevisiae recombination hotspots are not conserved over a broader evolutionary distance. Finally and strikingly, we found evidence of frequent reversal of commitment to meiosis, resulting in return to mitotic growth after allele shuffling. Identification of this major but underestimated evolutionary phenomenon illustrates the relevance of exploring non-model species.

Dissection of quantitative traits by bulk segregant mapping in a protoploid yeast species.

Since more than a decade ago, Saccharomyces cerevisiae has been used as a model to dissect complex traits, revealing the genetic basis of a large number of traits in fine detail. However, to have a more global view of the genetic architecture of traits across species, the examination of the molecular basis of phenotypes within non-conventional species would undoubtedly be valuable. In this respect, the Saccharomycotina yeasts represent ideal and potential non-model organisms. Here we sought to assess the feasibility of genetic mapping by bulk segregant analysis in the protoploid Lachancea kluyveri (formerly S. kluyveri) yeast species, a distantly related species to S. cerevisiae For this purpose, we designed a fluorescent mating-type marker, compatible with any mating-competent strains representative of this species, to rapidly create a large population of haploid segregants (>10(5) cells). Quantitative trait loci can be mapped by selecting and sequencing an enriched pool of progeny with extreme phenotypic values. As a test bed, we applied this strategy and mapped the causal loci underlying halotolerance phenotypes in L. kluyveri Overall, this study demonstrates that bulk segregant mapping is a powerful way for investigating the genetic basis of natural variations in non-model yeast organisms and more precisely in L. kluyveri.

Differences in environmental stress response among yeasts is consistent with species-specific lifestyles.

Defining how organisms respond to environmental change has always been an important step toward understanding their adaptive capacity and physiology. Variation in transcription during stress has been widely described in model species, especially in the yeast Saccharomyces cerevisiae, which helped to shape general rules regarding how cells cope with environmental constraints, as well as to decipher the functions of many genes. Comparison of the environmental stress response (ESR) across species is essential to obtaining better insight into the common and species-specific features of stress defense. In this context, we explored the transcriptional landscape of the yeast Lachancea kluyveri (formerly Saccharomyces kluyveri) in response to diverse stresses, using RNA sequencing. We investigated variation in gene expression and observed a link between genetic plasticity and environmental sensitivity. We identified the ESR genes in this species and compared them to those already found in S. cerevisiae We observed common features between the two species, as well as divergence in the regulatory networks involved. Of interest, some changes were related to differences in species lifestyle. Thus we were able to decipher how adaptation to stress has evolved among different yeast species. Finally, by analyzing patterns of coexpression, we were able to propose potential biological functions for 42% of genes and also annotate 301 genes for which no function could be assigned by homology. This large data set allowed for the characterization of the evolution of gene regulation and provides an efficient tool for assessing gene function.

Evolution of intraspecific transcriptomic landscapes in yeasts.

Variations in gene expression have been widely explored in order to obtain an accurate overview of the changes in regulatory networks that underlie phenotypic diversity. Numerous studies have characterized differences in genomic expression between large numbers of individuals of model organisms such as Saccharomyces cerevisiae. To more broadly survey the evolution of the transcriptomic landscape across species, we measured whole-genome expression in a large collection of another yeast species: Lachancea kluyveri (formerly Saccharomyces kluyveri), using RNAseq. Interestingly, this species diverged from the S. cerevisiae lineage prior to its ancestral whole genome duplication. Moreover, L. kluyveri harbors a chromosome-scale compositional heterogeneity due to a 1-Mb ancestral introgressed region as well as a large set of unique unannotated genes. In this context, our comparative transcriptomic analysis clearly showed a link between gene evolutionary history and expression behavior. Indeed, genes that have been recently acquired or under function relaxation tend to be less transcribed show a higher intraspecific variation (plasticity) and are less involved in network (connectivity). Moreover, utilizing this approach in L. kluyveri also highlighted specific regulatory network signatures in aerobic respiration, amino-acid biosynthesis and glycosylation, presumably due to its different lifestyle. Our data set sheds an important light on the evolution of intraspecific transcriptomic variation across distant species.

Deciphering regulatory variation of THI genes in alcoholic fermentation indicate an impact of Thi3p on PDC1 expression.

BACKGROUND: Thiamine availability is involved in glycolytic flux and fermentation efficiency. A deficiency of this vitamin may be responsible for sluggish fermentations in wine making. Therefore, both thiamine uptake and de novo synthesis could have key roles in fermentation processes. Thiamine biosynthesis is regulated in response to thiamine availability and is coordinated by the thiamine sensor Thi3p, which activates Pdc2p and Thi2p. We used a genetic approach to identify quantitative trait loci (QTLs) in wine yeast and we discovered that a set of thiamine genes displayed expression-QTL on a common locus, which contains the thiamine regulator THI3. RESULTS: We deciphered here the source of these regulatory variations of the THI and PDC genes. We showed that alteration of THI3 results in reduced expression of the genes involved in thiamine biosynthesis (THI11/12/13 and THI74) and increased expression of the pyruvate decarboxylase gene PDC1. Functional analysis of the allelic effect of THI3 confirmed the control of the THI and PDC1 genes. We observed, however, only a small effect of the THI3 on fermentation kinetics. We demonstrated that the expression levels of several THI genes are correlated with fermentation rate, suggesting that decarboxylation activity could drive gene expression through a modulation of thiamine content. Our data also reveals a new role of Thi3p in the regulation of the main pyruvate decarboxylase gene, PDC1. CONCLUSIONS: This highlights a switch from PDC1 to PDC5 gene expression during thiamine deficiency, which may improve the thiamine affinity or conservation during the enzymatic reaction. In addition, we observed that the lab allele of THI3 and of the thiamin transporter THI7 have diverged from the original alleles, consistent with an adaptation of lab strains to rich media containing an excess of thiamine.

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.

QTL mapping of the production of wine aroma compounds by yeast.

BACKGROUND: Wine aroma results from the combination of numerous volatile compounds, some produced by yeast and others produced in the grapes and further metabolized by yeast. However, little is known about the consequences of the genetic variation of yeast on the production of these volatile metabolites, or on the metabolic pathways involved in the metabolism of grape compounds. As a tool to decipher how wine aroma develops, we analyzed, under two experimental conditions, the production of 44 compounds by a population of 30 segregants from a cross between a laboratory strain and an industrial strain genotyped at high density. RESULTS: We detected eight genomic regions explaining the diversity concerning 15 compounds, some produced de novo by yeast, such as nerolidol, ethyl esters and phenyl ethanol, and others derived from grape compounds such as citronellol, and cis-rose oxide. In three of these eight regions, we identified genes involved in the phenotype. Hemizygote comparison allowed the attribution of differences in the production of nerolidol and 2-phenyl ethanol to the PDR8 and ABZ1 genes, respectively. Deletion of a PLB2 gene confirmed its involvement in the production of ethyl esters. A comparison of allelic variants of PDR8 and ABZ1 in a set of available sequences revealed that both genes present a higher than expected number of non-synonymous mutations indicating possible balancing selection. CONCLUSIONS: This study illustrates the value of QTL analysis for the analysis of metabolic traits, and in particular the production of wine aromas. It also identifies the particular role of the PDR8 gene in the production of farnesyldiphosphate derivatives, of ABZ1 in the production of numerous compounds and of PLB2 in ethyl ester synthesis. This work also provides a basis for elucidating the metabolism of various grape compounds, such as citronellol and cis-rose oxide.

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.

FSY1, a horizontally transferred gene in the Saccharomyces cerevisiae EC1118 wine yeast strain, encodes a high-affinity fructose/H+ symporter.

Transport of glucose and fructose in the yeast Saccharomyces cerevisiae plays a crucial role in controlling the rate of wine fermentation. In S. cerevisiae, hexoses are transported by facilitated diffusion via hexose carriers (Hxt), which prefer glucose to fructose. However, utilization of fructose by wine yeast is critically important at the end of fermentation. Here, we report the characterization of a fructose transporter recently identified by sequencing the genome of the commercial wine yeast strain EC1118 and found in many other wine yeasts. This transporter is designated Fsy1p because of its homology with the Saccharomyces pastorianus fructose/H(+) symporter Fsy1p. A strain obtained by transformation of the V5 hxt1-7Δ mutant with FSY1 grew well on fructose, but to a much lesser extent on glucose as the sole carbon source. Sugar uptake and symport experiments showed that FSY1 encodes a proton-coupled symporter with high affinity for fructose (K(m) 0.24±0.04mM). Using real-time RT-PCR, we also investigated the expression pattern of FSY1 in EC1118 growing on various carbon sources. FSY1 was repressed by high concentrations of glucose or fructose and was highly expressed on ethanol as the sole carbon source. The characteristics of this transporter indicate that its acquisition could confer a significant advantage to S. cerevisiae during the wine fermentation process. This transporter is a good example of acquisition of a new function in yeast by horizontal gene transfer.

Teolenn: an efficient and customizable workflow to design high-quality probes for microarray experiments.

Despite the development of new high-throughput sequencing techniques, microarrays are still attractive tools to study small genome organisms, thanks to sample multiplexing and high-feature densities. However, the oligonucleotide design remains a delicate step for most users. A vast array of software is available to deal with this problem, but each program is developed with its own strategy, which makes the choice of the best solution difficult. Here we describe Teolenn, a universal probe design workflow developed with a flexible and customizable module organization allowing fixed or variable length oligonucleotide generation. In addition, our software is able to supply quality scores for each of the designed probes. In order to assess the relevance of these scores, we performed a real hybridization using a tiling array designed against the Trichoderma reesei fungus genome. We show that our scoring pipeline correlates with signal quality for 97.2% of all the designed probes, allowing for a posteriori comparisons between quality scores and signal intensities. This result is useful in discarding any bad scoring probes during the design step in order to get high-quality microarrays. Teolenn is available at http://transcriptome.ens.fr/teolenn/.