Natural variation in non-coding regions underlying phenotypic diversity in budding yeast

Trans-eQTL hotspots shape complex traits by modulating cellular states

Regulatory genetic variation shapes gene expression, providing an important mechanism connecting DNA variation and complex traits. The causal relationships between gene expression and complex traits remain poorly understood. Here, we integrated transcriptomes and 46 genetically complex growth traits in a large cross between two strains of the yeast Saccharomyces cerevisiae. We discovered thousands of genetic correlations between gene expression and growth, suggesting potential functional connections. Local regulatory variation was a minor source of these genetic correlations. Instead, genetic correlations tended to arise from multiple independent trans-acting regulatory loci. Trans-acting hotspots that affect the expression of numerous genes accounted for particularly large fractions of genetic growth variation and of genetic correlations between gene expression and growth. Genes with genetic correlations were enriched for similar biological processes across traits but with heterogeneous direction of effect. Our results reveal how trans-acting regulatory hotspots shape complex traits by altering cellular states.

Phenotyping of a new yeast mapping population reveals differences in the activation of the TORC1 signalling pathway between wild and domesticated yeast strains

Domestication can be understood as a symbiotic relationship that benefits both domesticator and domesticated species, involving multiple genetic changes that configure the phenotype of the domesticated species. One of the most important domesticated species is the yeast Saccharomyces cerevisiae, with both domesticated strains used for different fermentations processes for thousands of years and wild strains existing only in environments without human intervention; however, little is known about the phenotypic effects associated with its domestication. In the present work, we studied the effect of domestication on yeast TORC1 activation, a pleiotropic signalling pathway conserved across the eukaryotic domain. To achieve this goal, we improved a previously generated methodology to assess TORC1 activation, which turned out to be as effective as the original one but also presents several practical advantages for its application (such as facilitating confirmation of transformants and putting the Luc reporter gene under the control of the same PRPL26A promoter for each transformed strain). We then generated a mapping population, the so-called TOMAN-G population, derived from the "1002 Yeast Genomes Project" population, the most comprehensive catalogue of the genetic variation in yeasts. Finally, strains belonging to the TOMAN-G population were phenotyped for TORC1 activation, and then we compared the results obtained between yeast strains with different ecological origins, finding differences in TORC1 activation between wild and domesticated strains, particularly wine strains. These results are indicative of the effect of domestication on TORC1 activation, specifically that the different evolutionary trajectories of wild and domesticated strains have in fact caused differences in the activation of this pathway; furthermore, the phenotypic data obtained in this work could be used to continue underlying the genetic bases of TORC1 activation, a process that is still not fully understood, using techniques such as GWAS to search for specific genetic variants underlying the observed phenotypic variability and phylogenetic tree inferences to gain insight into the evolutionary relationships between these genetic variants.

Signatures of selection in Angus and Hanwoo beef cattle using imputed whole genome sequence data

In this study, we detected signatures of selection in Hanwoo and Angus beef cattle using allele frequency and haplotype-based methods based on imputed whole genome sequence variants. Our dataset included 13,202 Angus animals with 10,057,633 imputed SNPs and 10,437 Hanwoo animals with 13,241,550 imputed SNPs. The dataset was subset down to 6,873,624 SNPs in common between the two populations to identify within population (runs of homozygosity, extended haplotype homozygosity) and between population signals of selection (allele fixation index, extended haplotype homozygosity). Assuming these selection signals were complementary to each other, they were combined into a decorrelated composite of multiple signals to identify regions under selection for each of the breeds. 27 genomic regions spanning 25.15 Mb and harboring 360 genes were identified in Angus on chromosomes 1,3, 4, 5, 6, 7, 8, 12, 13, 14, 16, 20, 21 and 28. Similarly, in Hanwoo, 59 genes and 17 genomic regions spanning 5.21 Mb on chromosomes 2, 4, 5, 6, 7, 8, 9, 10, 13, 17, 20 and 24 were identified. Apart from a small region on chromosome 13, there was no major overlap of selection signals between the two breeds reflecting their largely different selection histories, environmental challenges, breeding objectives and breed characteristics. Positional candidate genes identified in selected genomic regions in Angus have been previously associated with growth, immunity, reproductive development, feed efficiency and adaptation to environment while the candidate genes identified in Hanwoo included important genes regulating meat quality, fat deposition, cholesterol metabolism, lipid synthesis, neuronal development, and olfactory reception.

Inter-subspecies mouse F1 hybrid embryonic stem cell lines newly established for studies of allelic imbalance in gene expression

Allele-specific monoallelic gene expression is a unique phenomenon and a great resource for analyzing gene regulation. To study this phenomenon, we established new embryonic stem (ES) cell lines derived from F1 hybrid blastocysts from crosses between four mouse subspecies (Mus musculus domesticus, C57BL/6; M. musculus molossinus, MSM/Ms; M. musculus musculus, PWK; M. musculus castaneus, HMI/Ms) and analyzed the expression levels of undifferentiated pluripotent stem cell markers and karyotypes of each line. To demonstrate the utility of our cell lines, we analyzed the allele-specific expression pattern of the Inpp5d gene as an example. The allelic expression depended on the parental alleles; this dependence could be a consequence of differences in compatibility between cis- and trans-elements of the Inpp5d gene from different subspecies. The use of parental mice from four subspecies greatly enhanced genetic polymorphism. The F1 hybrid ES cells retained this polymorphism not only in the Inpp5d gene, but also at a genome-wide level. As we demonstrated for the Inpp5d gene, the established cell lines can contribute to the analysis of allelic expression imbalance based on the incompatibility between cis- and trans-elements and of phenotypes related to this incompatibility.

Trans-eQTL hotspots shape complex traits by modulating cellular states

Regulatory genetic variation shapes gene expression, providing an important mechanism connecting DNA variation and complex traits. The causal relationships between gene expression and complex traits remain poorly understood. Here, we integrated transcriptomes and 46 genetically complex growth traits in a large cross between two strains of the yeast Saccharomyces cerevisiae. We discovered thousands of genetic correlations between gene expression and growth, suggesting potential functional connections. Local regulatory variation was a minor source of these genetic correlations. Instead, genetic correlations tended to arise from multiple independent trans-acting regulatory loci. Trans-acting hotspots that affect the expression of numerous genes accounted for particularly large fractions of genetic growth variation and of genetic correlations between gene expression and growth. Genes with genetic correlations were enriched for similar biological processes across traits, but with heterogeneous direction of effect. Our results reveal how trans-acting regulatory hotspots shape complex traits by altering cellular states.

Genome-wide expression QTL mapping reveals the highly dynamic regulatory landscape of a major wheat pathogen

BACKGROUND

In agricultural ecosystems, outbreaks of diseases are frequent and pose a significant threat to food security. A successful pathogen undergoes a complex and well-timed sequence of regulatory changes to avoid detection by the host immune system; hence, well-tuned gene regulation is essential for survival. However, the extent to which the regulatory polymorphisms in a pathogen population provide an adaptive advantage is poorly understood.

RESULTS

We used Zymoseptoria tritici, one of the most important pathogens of wheat, to generate a genome-wide map of regulatory polymorphism governing gene expression. We investigated genome-wide transcription levels of 146 strains grown under nutrient starvation and performed expression quantitative trait loci (eQTL) mapping. We identified cis-eQTLs for 65.3% of all genes and the majority of all eQTL loci are within 2kb upstream and downstream of the transcription start site (TSS). We also show that polymorphism in different gene elements contributes disproportionally to gene expression variation. Investigating regulatory polymorphism in gene categories, we found an enrichment of regulatory variants for genes predicted to be important for fungal pathogenesis but with comparatively low effect size, suggesting a separate layer of gene regulation involving epigenetics. We also show that previously reported trait-associated SNPs in pathogen populations are frequently cis-regulatory variants of neighboring genes with implications for the trait architecture.

CONCLUSIONS

Overall, our study provides extensive evidence that single populations segregate large-scale regulatory variation and are likely to fuel rapid adaptation to resistant hosts and environmental change.

Single nucleotide polymorphisms associated with wine fermentation and adaptation to nitrogen limitation in wild and domesticated yeast strains

For more than 20 years, Saccharomyces cerevisiae has served as a model organism for genetic studies and molecular biology, as well as a platform for biotechnology (e.g., wine production). One of the important ecological niches of this yeast that has been extensively studied is wine fermentation, a complex microbiological process in which S. cerevisiae faces various stresses such as limited availability of nitrogen. Nitrogen deficiencies in grape juice impair fermentation rate and yeast biomass production, leading to sluggish or stuck fermentations, resulting in considerable economic losses for the wine industry. In the present work, we took advantage of the "1002 Yeast Genomes Project" population, the most complete catalogue of the genetic variation in the species and a powerful resource for genotype-phenotype correlations, to study the adaptation to nitrogen limitation in wild and domesticated yeast strains in the context of wine fermentation. We found that wild and domesticated yeast strains have different adaptations to nitrogen limitation, corroborating their different evolutionary trajectories. Using a combination of state-of-the-art bioinformatic (GWAS) and molecular biology (CRISPR-Cas9) methodologies, we validated that PNP1, RRT5 and PDR12 are implicated in wine fermentation, where RRT5 and PDR12 are also involved in yeast adaptation to nitrogen limitation. In addition, we validated SNPs in these genes leading to differences in fermentative capacities and adaptation to nitrogen limitation. Altogether, the mapped genetic variants have potential applications for the genetic improvement of industrial yeast strains.

Impact of the acquired subgenome on the transcriptional landscape in Brettanomyces bruxellensis allopolyploids

Gene expression variation can provide an overview of the changes in regulatory networks that underlie phenotypic diversity. Certain evolutionary trajectories such as polyploidization events can have an impact on the transcriptional landscape. Interestingly, the evolution of the yeast species Brettanomyces bruxellensis has been punctuated by diverse allopolyploidization events leading to the coexistence of a primary diploid genome associated with various haploid acquired genomes. To assess the impact of these events on gene expression, we generated and compared the transcriptomes of a set of 87 B. bruxellensis isolates, selected as being representative of the genomic diversity of this species. Our analysis revealed that acquired subgenomes strongly impact the transcriptional patterns and allow discrimination of allopolyploid populations. In addition, clear transcriptional signatures related to specific populations have been revealed. The transcriptional variations observed are related to some specific biological processes such as transmembrane transport and amino acids metabolism. Moreover, we also found that the acquired subgenome causes the overexpression of some genes involved in the production of flavor-impacting secondary metabolites, especially in isolates of the beer population.

A Saccharomyces eubayanus haploid resource for research studies

Since its identification, Saccharomyces eubayanus has been recognized as the missing parent of the lager hybrid, S. pastorianus. This wild yeast has never been isolated from fermentation environments, thus representing an interesting candidate for evolutionary, ecological and genetic studies. However, it is imperative to develop additional molecular genetics tools to ease manipulation and thus facilitate future studies. With this in mind, we generated a collection of stable haploid strains representative of three main lineages described in S. eubayanus (PB-1, PB-2 and PB-3), by deleting the HO gene using CRISPR-Cas9 and tetrad micromanipulation. Phenotypic characterization under different conditions demonstrated that the haploid derivates were extremely similar to their parental strains. Genomic analysis in three strains highlighted a likely low frequency of off-targets, and sequencing of a single tetrad evidenced no structural variants in any of the haploid spores. Finally, we demonstrate the utilization of the haploid set by challenging the strains under mass-mating conditions. In this way, we found that S. eubayanus under liquid conditions has a preference to remain in a haploid state, unlike S. cerevisiae that mates rapidly. This haploid resource is a novel set of strains for future yeast molecular genetics studies.

Comprehensive Genomic Discovery of Non-Coding Transcriptional Enhancers in the African Malaria Vector Anopheles coluzzii

Almost all regulation of gene expression in eukaryotic genomes is mediated by the action of distant non-coding transcriptional enhancers upon proximal gene promoters. Enhancer locations cannot be accurately predicted bioinformatically because of the absence of a defined sequence code, and thus functional assays are required for their direct detection. Here we used a massively parallel reporter assay, Self-Transcribing Active Regulatory Region sequencing (STARR-seq), to generate the first comprehensive genome-wide map of enhancers in Anopheles coluzzii, a major African malaria vector in the Gambiae species complex. The screen was carried out by transfecting reporter libraries created from the genomic DNA of 60 wild A. coluzzii from Burkina Faso into A. coluzzii 4a3A cells, in order to functionally query enhancer activity of the natural population within the homologous cellular context. We report a catalog of 3,288 active genomic enhancers that were significant across three biological replicates, 74% of them located in intergenic and intronic regions. The STARR-seq enhancer screen is chromatin-free and thus detects inherent activity of a comprehensive catalog of enhancers that may be restricted in vivo to specific cell types or developmental stages. Testing of a validation panel of enhancer candidates using manual luciferase assays confirmed enhancer function in 26 of 28 (93%) of the candidates over a wide dynamic range of activity from two to at least 16-fold activity above baseline. The enhancers occupy only 0.7% of the genome, and display distinct composition features. The enhancer compartment is significantly enriched for 15 transcription factor binding site signatures, and displays divergence for specific dinucleotide repeats, as compared to matched non-enhancer genomic controls. The genome-wide catalog of A. coluzzii enhancers is publicly available in a simple searchable graphic format. This enhancer catalogue will be valuable in linking genetic and phenotypic variation, in identifying regulatory elements that could be employed in vector manipulation, and in better targeting of chromosome editing to minimize extraneous regulation influences on the introduced sequences.

Importance: Understanding the role of the non-coding regulatory genome in complex disease phenotypes is essential, but even in well-characterized model organisms, identification of regulatory regions within the vast non-coding genome remains a challenge. We used a large-scale assay to generate a genome wide map of transcriptional enhancers. Such a catalogue for the important malaria vector, Anopheles coluzzii, will be an important research tool as the role of non-coding regulatory variation in differential susceptibility to malaria infection is explored and as a public resource for research on this important insect vector of disease.

Uncovering Divergence in Gene Expression Regulation in the Adaptation of Yeast to Nitrogen Scarcity

Saccharomyces cerevisiae rewires its transcriptional output to survive stressful environments, such as nitrogen scarcity under fermentative conditions. Although divergence in nitrogen metabolism among natural yeast populations has been reported, the impact of regulatory genetic variants modulating gene expression and nitrogen consumption remains to be investigated. Here, we employed an F1 hybrid from two contrasting S. cerevisiae strains, providing a controlled genetic environment to map cis factors involved in the divergence of gene expression regulation in response to nitrogen scarcity. We used a dual approach to obtain genome-wide allele-specific profiles of chromatin accessibility, transcription factor binding, and gene expression through ATAC-seq (assay for transposase accessible chromatin) and RNA-seq (transcriptome sequencing). We observed large variability in allele-specific expression and accessibility between the two genetic backgrounds, with a third of these differences specific to a deficient nitrogen environment. Furthermore, we discovered events of allelic bias in gene expression correlating with allelic bias in transcription factor binding solely under nitrogen scarcity, where the majority of these transcription factors orchestrates the nitrogen catabolite repression regulatory pathway and demonstrates a cis × environment-specific response. Our approach allowed us to find cis variants modulating gene expression, chromatin accessibility, and allelic differences in transcription factor binding in response to low nitrogen culture conditions. IMPORTANCE Historically, coding variants were prioritized when searching for causal mechanisms driving adaptation of natural populations to stressful environments. However, the recent focus on noncoding variants demonstrated their ubiquitous role in adaptation. Here, we performed genome-wide regulatory variation profiles between two divergent yeast strains when facing nitrogen nutritional stress. The open chromatin availability of several regulatory regions changes in response to nitrogen scarcity. Importantly, we describe regulatory events that deviate between strains. Our results demonstrate a widespread variation in gene expression regulation between naturally occurring populations in response to stressful environments.

The mode of expression divergence in Drosophila fat body is infection-specific

Transcription is controlled by interactions of cis-acting DNA elements with diffusible trans-acting factors. Changes in cis or trans factors can drive expression divergence within and between species, and their relative prevalence can reveal the evolutionary history and pressures that drive expression variation. Previous work delineating the mode of expression divergence in animals has largely used whole-body expression measurements in one condition. Because cis-acting elements often drive expression in a subset of cell types or conditions, these measurements may not capture the complete contribution of cis-acting changes. Here, we quantify the mode of expression divergence in the Drosophila fat body, the primary immune organ, in several conditions, using two geographically distinct lines of D. melanogaster and their F1 hybrids. We measured expression in the absence of infection and in infections with Gram-negative S. marcescens or Gram-positive E. faecalis bacteria, which trigger the two primary signaling pathways in the Drosophila innate immune response. The mode of expression divergence strongly depends on the condition, with trans-acting effects dominating in response to Gram-negative infection and cis-acting effects dominating in Gram-positive and preinfection conditions. Expression divergence in several receptor proteins may underlie the infection-specific trans effects. Before infection, when the fat body has a metabolic role, there are many compensatory effects, changes in cis and trans that counteract each other to maintain expression levels. This work shows that within a single tissue, the mode of expression divergence varies between conditions and suggests that these differences reflect the diverse evolutionary histories of host-pathogen interactions.

QTL analysis of natural Saccharomyces cerevisiae isolates reveals unique alleles involved in lignocellulosic inhibitor tolerance

Decoding the genetic basis of lignocellulosic inhibitor tolerance in Saccharomyces cerevisiae is crucial for rational engineering of bioethanol strains with enhanced robustness. The genetic diversity of natural strains present an invaluable resource for the exploration of complex traits of industrial importance from a pan-genomic perspective to complement the limited range of specialised, tolerant industrial strains. Natural S. cerevisiae isolates have lately garnered interest as a promising toolbox for engineering novel, genetically encoded tolerance phenotypes into commercial strains. To this end, we investigated the genetic basis for lignocellulosic inhibitor tolerance of natural S. cerevisiae isolates. A total of 12 quantitative trait loci underpinning tolerance were identified by next-generation sequencing linked bulk-segregant analysis of superior interbred pools. Our findings corroborate the current perspective of lignocellulosic inhibitor tolerance as a multigenic, complex trait. Apart from a core set of genetic variants required for inhibitor tolerance, an additional genetic background-specific response was observed. Functional analyses of the identified genetic loci revealed the uncharacterised ORF, YGL176C and the bud-site selection XRN1/BUD13 as potentially beneficial alleles contributing to tolerance to a complex lignocellulosic inhibitor mixture. We present evidence for the consideration of both regulatory and coding sequence variants for strain improvement.

Differential Gene Expression and Allele Frequency Changes Favour Adaptation of a Heterogeneous Yeast Population to Nitrogen-Limited Fermentations

Alcoholic fermentation is fundamentally an adaptation process, in which the yeast Saccharomyces cerevisiae outperforms its competitors and takes over the fermentation process itself. Although wine yeast strains appear to be adapted to the stressful conditions of alcoholic fermentation, nitrogen limitations in grape must cause stuck or slow fermentations, generating significant economic losses for the wine industry. One way to discover the genetic bases that promote yeast adaptation to nitrogen-deficient environments are selection experiments, where a yeast population undergoes selection under conditions of nitrogen restriction for a number of generations, to then identify by sequencing the molecular characteristics that promote this adaptation. In this work, we carried out selection experiments in bioreactors imitating wine fermentation under nitrogen-limited fermentation conditions (SM60), using the heterogeneous SGRP-4X yeast population, to then sequence the transcriptome and the genome of the population at different time points of the selection process. The transcriptomic results showed an overexpression of genes from the NA strain (North American/YPS128), a wild, non-domesticated isolate. In addition, genome sequencing and allele frequency results allowed several QTLs to be mapped for adaptation to nitrogen-limited fermentation. Finally, we validated the ECM38 allele of NA strain as responsible for higher growth efficiency under nitrogen-limited conditions. Taken together, our results revealed a complex pattern of molecular signatures favouring adaptation of the yeast population to nitrogen-limited fermentations, including differential gene expression, allele frequency changes and loss of the mitochondrial genome. Finally, the results suggest that wild alleles from a non-domesticated isolate (NA) may have a relevant role in the adaptation to the assayed fermentation conditions, with the consequent potential of these alleles for the genetic improvement of wine yeast strains.

GTR1 Affects Nitrogen Consumption and TORC1 Activity in Saccharomyces cerevisiae Under Fermentation Conditions

The TORC1 pathway coordinates cell growth in response to nitrogen availability present in the medium, regulating genes related to nitrogen transport and metabolism. Therefore, the adaptation of Saccharomyces cerevisiae to changes in nitrogen availability implies variations in the activity of this signaling pathway. In this sense, variations in nitrogen detection and signaling pathway are one of the main causes of differences in nitrogen assimilation during alcoholic fermentation. Previously, we demonstrated that allelic variants in the GTR1 gene underlying differences in ammonium and amino acids consumption between Wine/European (WE) and West African (WA) strains impact the expression of nitrogen transporters. The GTR1 gene encodes a GTPase that participates in the EGO complex responsible for TORC1 activation in response to amino acids availability. In this work, we assessed the role of the GTR1 gene on nitrogen consumption under fermentation conditions, using a high sugar concentration medium with nitrogen limitation and in the context of the WE and WA genetic backgrounds. The gtr1Δ mutant presented a reduced TORC1 activity and increased expression levels of nitrogen transporters, which in turn favored ammonium consumption, but decreased amino acid assimilation. Furthermore, to identify the SNPs responsible for differences in nitrogen consumption during alcoholic fermentation, we studied the polymorphisms present in the GTR1 gene. We carried out swapping experiments for the promoter and coding regions of GTR1 between the WE and WA strains. We observed that polymorphisms in the coding region of the WA GTR1 gene are relevant for TORC1 activity. Altogether, our results highlight the role of the GTR1 gene on nitrogen consumption in S. cerevisiae under fermentation conditions.

The Transcriptional Aftermath in Two Independently Formed Hybrids of the Opportunistic Pathogen Candida orthopsilosis

How new pathogens emerge is an important question that remains largely unanswered. Some emerging yeast pathogens are hybrids originated through the crossing of two different species, but how hybridization contributes to higher virulence is unclear. Here, we show that hybrids selectively retain gene regulation plasticity inherited from the two parents and that this plasticity affects genes involved in virulence.

Interspecific hybridization can drive evolutionary adaptation to novel environments. The Saccharomycotina clade of budding yeasts includes many hybrid lineages, and hybridization has been proposed as a source for new pathogenic species. Candida orthopsilosis is an emerging opportunistic pathogen for which most clinical isolates are hybrids, each derived from one of at least four independent crosses between the same two parental lineages. To gain insight into the transcriptomic aftermath of hybridization in these pathogens, we analyzed allele-specific gene expression in two independently formed hybrid strains and in a homozygous strain representative of one parental lineage. Our results show that the effect of hybridization on overall gene expression is rather limited, affecting ∼4% of the genes studied. However, we identified a larger effect in terms of imbalanced allelic expression, affecting ∼9.5% of the heterozygous genes in the hybrids. This effect was larger in the hybrid with more extensive loss of heterozygosity, which may indicate a tendency to avoid loss of heterozygosity in these genes. Consistently, the number of shared genes with allele-specific expression in the two independently formed hybrids was higher than random expectation, suggesting selective retention. Some of the imbalanced genes have functions related to pathogenicity, including zinc transport and superoxide dismutase activities. While it remains unclear whether the observed imbalanced genes play a role in virulence, our results suggest that differences in allele-specific expression may add an additional layer of phenotypic plasticity to traits related to virulence in C. orthopsilosis hybrids.

IMPORTANCE How new pathogens emerge is an important question that remains largely unanswered. Some emerging yeast pathogens are hybrids originated through the crossing of two different species, but how hybridization contributes to higher virulence is unclear. Here, we show that hybrids selectively retain gene regulation plasticity inherited from the two parents and that this plasticity affects genes involved in virulence.

Transcriptional Activity and Protein Levels of Horizontally Acquired Genes in Yeast Reveal Hallmarks of Adaptation to Fermentative Environments

In the past decade, the sequencing of large cohorts of Saccharomyces cerevisiae strains has revealed a landscape of genomic regions acquired by Horizontal Gene Transfer (HGT). The genes acquired by HGT play important roles in yeast adaptation to the fermentation process, improving nitrogen and carbon source utilization. However, the functional characterization of these genes at the molecular level has been poorly attended. In this work, we carried out a systematic analysis of the promoter activity and protein level of 30 genes contained in three horizontally acquired regions commonly known as regions A, B, and C. In three strains (one for each region), we used the luciferase reporter gene and the mCherry fluorescent protein to quantify the transcriptional and translational activity of these genes, respectively. We assayed the strains generated in four different culture conditions; all showed low levels of transcriptional and translational activity across these environments. However, we observed an increase in protein levels under low nitrogen culture conditions, suggesting a possible role of the horizontally acquired genes in the adaptation to nitrogen-limited environments. Furthermore, since the strains carrying the luciferase reporter gene are null mutants for the horizontally acquired genes, we assayed growth parameters (latency time, growth rate, and efficiency) and the fermentation kinetics in this set of deletion strains. The results showed that single deletion of 20 horizontally acquired genes modified the growth parameters, whereas the deletion of five of them altered the maximal CO₂ production rate (Vmax). Interestingly, we observed a correlation between growth parameters and Vmax for an ORF within region A, encoding an ortholog to a thiamine (vitamin B1) transporter whose deletion decreased the growth rate, growth efficiency, and CO₂ production. Altogether, our results provided molecular and phenotypic evidence highlighting the importance of horizontally acquired genes in yeast adaptation to fermentative environments.

Disentangling the genetic bases of Saccharomyces cerevisiae nitrogen consumption and adaptation to low nitrogen environments in wine fermentation

The budding yeast Saccharomyces cerevisiae has been considered for more than 20 years as a premier model organism for biological sciences, also being the main microorganism used in wide industrial applications, like alcoholic fermentation in the winemaking process. Grape juice is a challenging environment for S. cerevisiae, with nitrogen deficiencies impairing fermentation rate and yeast biomass production, causing stuck or sluggish fermentations, thus generating sizeable economic losses for wine industry. In the present review, we summarize some recent efforts in the search of causative genes that account for yeast adaptation to low nitrogen environments, specially focused in wine fermentation conditions. We start presenting a brief perspective of yeast nitrogen utilization under wine fermentative conditions, highlighting yeast preference for some nitrogen sources above others. Then, we give an outlook of S. cerevisiae genetic diversity studies, paying special attention to efforts in genome sequencing for population structure determination and presenting QTL mapping as a powerful tool for phenotype-genotype correlations. Finally, we do a recapitulation of S. cerevisiae natural diversity related to low nitrogen adaptation, specially showing how different studies have left in evidence the central role of the TORC1 signalling pathway in nitrogen utilization and positioned wild S. cerevisiae strains as a reservoir of beneficial alleles with potential industrial applications (e.g. improvement of industrial yeasts for wine production). More studies focused in disentangling the genetic bases of S. cerevisiae adaptation in wine fermentation will be key to determine the domestication effects over low nitrogen adaptation, as well as to definitely proof that wild S. cerevisiae strains have potential genetic determinants for better adaptation to low nitrogen conditions.

Comparison of Phylogenetic Tree Topologies for Nitrogen Associated Genes Partially Reconstruct the Evolutionary History of Saccharomyces cerevisiae

Massive sequencing projects executed in Saccharomyces cerevisiae have revealed in detail its population structure. The recent "1002 yeast genomes project" has become the most complete catalogue of yeast genetic diversity and a powerful resource to analyse the evolutionary history of genes affecting specific phenotypes. In this work, we selected 22 nitrogen associated genes and analysed the sequence information from the 1011 strains of the "1002 yeast genomes project". We constructed a total evidence (TE) phylogenetic tree using concatenated information, which showed a 27% topology similarity with the reference (REF) tree of the "1002 yeast genomes project". We also generated individual phylogenetic trees for each gene and compared their topologies, identifying genes with similar topologies (suggesting a shared evolutionary history). Furthermore, we pruned the constructed phylogenetic trees to compare the REF tree topology versus the TE tree and the individual genes trees, considering each phylogenetic cluster/subcluster within the population, observing genes with cluster/subcluster topologies of high similarity to the REF tree. Finally, we used the pruned versions of the phylogenetic trees to compare four strains considered as representatives of S. cerevisiae clean lineages, observing for 15 genes that its cluster topologies match 100% the REF tree, supporting that these strains represent main lineages of yeast population. Altogether, our results showed the potential of tree topologies comparison for exploring the evolutionary history of a specific group of genes.

Divergent Allometric Trajectories in Gene Expression and Coexpression Produce Species Differences in Sympatrically Speciating Midas Cichlid Fish

The mechanisms of speciation without geographic isolation (i.e., sympatric speciation) remain debated. This is due in part to the fact that the genomic landscape that could promote or hinder species divergence in the presence of gene flow is still largely unknown. However, intensive research is now centered on understanding the genetic architecture of adaptive traits associated with this process as well as how gene expression might affect these traits. Here, using RNA-Seq data, we investigated gene expression of sympatrically speciating benthic and limnetic Neotropical cichlid fishes at two developmental stages. First, we identified groups of coexpressed genes (modules) at each stage. Although there are a few large and well-preserved modules, most of the other modules are not preserved across life stages. Second, we show that later in development more and larger coexpression modules are associated with divergence between benthic and limnetic fish compared with the earlier life stage. This divergence between benthic and limnetic fish in coexpression mirrors divergence in overall expression between benthic and limnetic fish, which is more pronounced later in life. Our results reveal that already at 1-day posthatch benthic and limnetic fish diverge in (co)expression, and that this divergence becomes more substantial when fish are free-swimming but still unlikely to have divergent swimming and feeding habits. More importantly, our study describes how the coexpression of several genes through development, as opposed to individual genes, is associated with benthic–limnetic species differences, and how two morphogenetic trajectories diverge as fish grow older.

Compensatory trans-regulatory alleles minimizing variation in TDH3 expression are common within Saccharomyces cerevisiae

Heritable variation in gene expression is common within species. Much of this variation is due to genetic differences outside of the gene with altered expression and is trans-acting. This trans-regulatory variation is often polygenic, with individual variants typically having small effects, making the genetic architecture and evolution of trans-regulatory variation challenging to study. Consequently, key questions about trans-regulatory variation remain, including the variability of trans-regulatory variation within a species, how selection affects trans-regulatory variation, and how trans-regulatory variants are distributed throughout the genome and within a species. To address these questions, we isolated and measured trans-regulatory differences affecting TDH3 promoter activity among 56 strains of Saccharomyces cerevisiae, finding that trans-regulatory backgrounds varied approximately twofold in their effects on TDH3 promoter activity. Comparing this variation to neutral models of trans-regulatory evolution based on empirical measures of mutational effects revealed that despite this variability in the effects of trans-regulatory backgrounds, stabilizing selection has constrained trans-regulatory differences within this species. Using a powerful quantitative trait locus mapping method, we identified ∼100 trans-acting expression quantitative trait locus in each of three crosses to a common reference strain, indicating that regulatory variation is more polygenic than previous studies have suggested. Loci altering expression were located throughout the genome, and many loci were strain specific. This distribution and prevalence of alleles is consistent with recent theories about the genetic architecture of complex traits. In all mapping experiments, the nonreference strain alleles increased and decreased TDH3 promoter activity with similar frequencies, suggesting that stabilizing selection maintained many trans-acting variants with opposing effects. This variation may provide the raw material for compensatory evolution and larger scale regulatory rewiring observed in developmental systems drift among species.

KAE1 Allelic Variants Affect TORC1 Activation and Fermentation Kinetics in Saccharomyces cerevisiae

The eukaryotic domain-conserved TORC1 signalling pathway connects growth with nutrient sufficiency, promoting anabolic processes such as ribosomal biogenesis and protein synthesis. In Saccharomyces cerevisiae, TORC1 is activated mainly by the nitrogen sources. Recently, this pathway has gotten renewed attention but now in the context of the alcoholic fermentation, due to its key role in nitrogen metabolism regulation. Although the distal and proximal effectors downstream TORC1 are well characterised in yeast, the mechanism by which TORC1 is activated by nitrogen sources is not fully understood. In this work, we took advantage of a previously developed microculture-based methodology, which indirectly evaluates TORC1 activation in a nitrogen upshift experiment, to identify genetic variants affecting the activation of this pathway. We used this method to phenotype a recombinant population derived from two strains (SA and WE) with different geographic origins, which show opposite phenotypes for TORC1 activation by glutamine. Using this phenotypic information, we performed a QTL mapping that allowed us to identify several QTLs for TORC1 activation. Using a reciprocal hemizygous analysis, we validated GUS1, KAE1, PIB2, and UTH1 as genes responsible for the natural variation in the TORC1 activation. We observed that reciprocal hemizygous strains for KAE1 (ATPase required for t⁶A tRNA modification) gene showed the greatest phenotypic differences for TORC1 activation, with the hemizygous strain carrying the SA allele (KAE1^SA) showing the higher TORC1 activation. In addition, we evaluated the fermentative capacities of the hemizygous strains under low nitrogen conditions, observing an antagonistic effect for KAE1^SA allele, where the hemizygous strain containing this allele presented the lower fermentation rate. Altogether, these results highlight the importance of the tRNA processing in TORC1 activation and connects this pathway with the yeasts fermentation kinetics under nitrogen-limited conditions.

Genetic variants of TORC1 signaling pathway affect nitrogen consumption in Saccharomyces cerevisiae during alcoholic fermentation

In the alcoholic fermentation process, Saccharomyces cerevisiae strains present differences in their nitrogen consumption profiles, these phenotypic outcomes have complex genetic and molecular architectures. In this sense, variations in nitrogen signaling pathways regulated by TORC1 represent one of the main sources of phenotypic diversity in nitrogen consumption. This emphasizes the possible roles that allelic variants from the TORC1 pathway have in the nitrogen consumption differences observed in yeast during the alcoholic fermentation. Here, we studied the allelic diversity in the TORC1 pathway across four yeast strains and determined how these polymorphisms directly impact nitrogen consumption during alcoholic fermentation. Using a reciprocal hemizygosity approach combined with phenotyping under fermentative conditions, we found that allelic variants of GTR1, TOR2, SIT4, SAP185, EAP1, NPR1 and SCH9 underlie differences in the ammonium and amino acids consumption phenotypes. Among these, GTR1 alleles from the Wine/European and West African genetic backgrounds showed the greatest effects on ammonium and amino acid consumption, respectively. Furthermore, we identified allelic variants of SAP185, TOR2, SCH9 and NPR1 from an oak isolate that increased the amino acid consumption preference over ammonium; representing putative candidates coming from a non-domesticated strain that could be used for genetic improvement programs. In conclusion, our results demonstrated that a large number of allelic variants within the TORC1 pathway significantly impacts on regulatory mechanisms of nitrogen assimilation during alcoholic fermentation.

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.

Indirect monitoring of TORC1 signalling pathway reveals molecular diversity among different yeast strains

Saccharomyces cerevisiae is the main species responsible for the alcoholic fermentation in wine production. One of the main problems in this process is the deficiency of nitrogen sources in the grape must, which can lead to stuck or sluggish fermentations. Currently, yeast nitrogen consumption and metabolism are under active inquiry, with emphasis on the study of the TORC1 signalling pathway, given its central role responding to nitrogen availability and influencing growth and cell metabolism. However, the mechanism by which different nitrogen sources activates TORC1 is not completely understood. Existing methods to evaluate TORC1 activation by nitrogen sources are time-consuming, making difficult the analyses of large numbers of strains. In this work, a new indirect method for monitoring TORC1 pathway was developed on the basis of the luciferase reporter gene controlled by the promoter region of RPL26A gene, a gene known to be expressed upon TORC1 activation. The method was tested in strains representative of the clean lineages described so far in S. cerevisiae. The activation of the TORC1 pathway by a proline-to-glutamine upshift was indirectly evaluated using our system and the traditional direct methods based on immunoblot (Sch9 and Rps6 phosphorylation). Regardless of the different molecular readouts obtained with both methodologies, the general results showed a wide phenotypic variation between the representative strains analysed. Altogether, this easy-to-use assay opens the possibility to study the molecular basis for the differential TORC1 pathway activation, allowing to interrogate a larger number of strains in the context of nitrogen metabolism phenotypic differences.

Fungal Light-Oxygen-Voltage Domains for Optogenetic Control of Gene Expression and Flocculation in Yeast

Optogenetic switches permit accurate control of gene expression upon light stimulation. These synthetic switches have become a powerful tool for gene regulation, allowing modulation of customized phenotypes, overcoming the obstacles of chemical inducers, and replacing their use by an inexpensive resource: light. In this work, we implemented FUN-LOV, an optogenetic switch based on the photon-regulated interaction of WC-1 and VVD, two LOV (light-oxygen-voltage) blue-light photoreceptors from the fungus Neurospora crassa. When tested in yeast, FUN-LOV yields light-controlled gene expression with exquisite temporal resolution and a broad dynamic range of over 1,300-fold, as measured by a luciferase reporter. We also tested the FUN-LOV switch for heterologous protein expression in Saccharomyces cerevisiae, where Western blot analysis confirmed strong induction upon light stimulation, surpassing by 2.5 times the levels achieved with a classic GAL4/galactose chemical-inducible system. Additionally, we utilized FUN-LOV to control the ability of yeast cells to flocculate. Light-controlled expression of the flocculin-encoding gene FLO1, by the FUN-LOV switch, yielded flocculation in light (FIL), whereas the light-controlled expression of the corepressor TUP1 provided flocculation in darkness (FID). Altogether, the results reveal the potential of the FUN-LOV optogenetic switch to control two biotechnologically relevant phenotypes such as heterologous protein expression and flocculation, paving the road for the engineering of new yeast strains for industrial applications. Importantly, FUN-LOV’s ability to accurately manipulate gene expression, with a high temporal dynamic range, can be exploited in the analysis of diverse biological processes in various organisms.

Optogenetic switches are molecular devices which allow the control of different cellular processes by light, such as gene expression, providing a versatile alternative to chemical inducers. Here, we report a novel optogenetic switch (FUN-LOV) based on the LOV domain interaction of two blue-light photoreceptors (WC-1 and VVD) from the fungus N. crassa. In yeast cells, FUN-LOV allowed tight regulation of gene expression, with low background in darkness and a highly dynamic and potent control by light. We used FUN-LOV to optogenetically manipulate, in yeast, two biotechnologically relevant phenotypes, heterologous protein expression and flocculation, resulting in strains with potential industrial applications. Importantly, FUN-LOV can be implemented in diverse biological platforms to orthogonally control a multitude of cellular processes.

Transposable element insertions shape gene regulation and melanin production in a fungal pathogen of wheat

BACKGROUND

Fungal plant pathogens pose major threats to crop yield and sustainable food production if they are highly adapted to their host and the local environment. Variation in gene expression contributes to phenotypic diversity within fungal species and affects adaptation. However, very few cases of adaptive regulatory changes have been reported in fungi and the underlying mechanisms remain largely unexplored. Fungal pathogen genomes are highly plastic and harbor numerous insertions of transposable elements, which can potentially contribute to gene expression regulation. In this work, we elucidated how transposable elements contribute to variation in melanin accumulation, a quantitative trait in fungi that affects survival under stressful conditions.

RESULTS

We demonstrated that differential transcriptional regulation of the gene encoding the transcription factor Zmr1, which controls expression of the genes in the melanin biosynthetic gene cluster, is responsible for variation in melanin accumulation in the fungal plant pathogen Zymoseptoria tritici. We show that differences in melanin levels between two strains of Z. tritici are due to two levels of transcriptional regulation: (1) variation in the promoter sequence of Zmr1 and (2) an insertion of transposable elements upstream of the Zmr1 promoter. Remarkably, independent insertions of transposable elements upstream of Zmr1 occurred in 9% of Z. tritici strains from around the world and negatively regulated Zmr1 expression, contributing to variation in melanin accumulation.

CONCLUSIONS

Our studies identified two levels of transcriptional control that regulate the synthesis of melanin. We propose that these regulatory mechanisms evolved to balance the fitness costs associated with melanin production against its positive contribution to survival in stressful environments.

GPD1 and ADH3 Natural Variants Underlie Glycerol Yield Differences in Wine Fermentation

Glycerol is one of the most important by-products of alcohol fermentation, and depending on its concentration it can contribute to wine flavor intensity and aroma volatility. Here, we evaluated the potential of utilizing the natural genetic variation of non-coding regions in budding yeast to identify allelic variants that could modulate glycerol phenotype during wine fermentation. For this we utilized four Saccharomyces cerevisiae strains (WE - Wine/European, SA – Sake, NA – North American, and WA – West African), which were previously profiled for genome-wide Allele Specific Expression (ASE) levels. The glycerol yields under Synthetic Wine Must (SWM) fermentations differed significantly between strains; WA produced the highest glycerol yields while SA produced the lowest yields. Subsequently, from our ASE database, we identified two candidate genes involved in alcoholic fermentation pathways, ADH3 and GPD1, exhibiting significant expression differences between strains. A reciprocal hemizygosity assay demonstrated that hemizygotes expressing GPD1^WA, GPD1^SA, ADH3^WA and ADH3^SA alleles had significantly greater glycerol yields compared to GPD1^WE and ADH3^WE. We further analyzed the gene expression profiles for each GPD1 variant under SWM, demonstrating that the expression of GPD1^WE occurred earlier and was greater compared to the other alleles. This result indicates that the level, timing, and condition of expression differ between regulatory regions in the various genetic backgrounds. Furthermore, promoter allele swapping demonstrated that these allele expression patterns were transposable across genetic backgrounds; however, glycerol yields did not differ between wild type and modified strains, suggesting a strong trans effect on GPD1 gene expression. In this line, Gpd1 protein levels in parental strains, particularly Gpd1p^WE, did not necessarily correlate with gene expression differences, but rather with glycerol yield where low Gpd1p^WE levels were detected. This suggests that GPD1^WE is influenced by recessive negative post-transcriptional regulation which is absent in the other genetic backgrounds. This dissection of regulatory mechanisms in GPD1 allelic variants demonstrates the potential to exploit natural alleles to improve glycerol production in wine fermentation and highlights the difficulties of trait improvement due to alternative trans-regulation and gene-gene interactions in the different genetic background.

Evolutionary divergence of 3' UTRs in cichlid fishes

Background

Post-transcriptional regulation is crucial for the control of eukaryotic gene expression and might contribute to adaptive divergence. The three prime untranslated regions (3’ UTRs), that are located downstream of protein-coding sequences, play important roles in post-transcriptional regulation. These regions contain functional elements that influence the fate of mRNAs and could be exceptionally important in groups such as rapidly evolving cichlid fishes.

Results

To examine cichlid 3’ UTR evolution, we 1) identified gene features in nine teleost genomes and 2) performed comparative analyses to assess evolutionary variation in length, functional motifs, and evolutionary rates of 3’ UTRs. In all nine teleost genomes, we found a smaller proportion of repetitive elements in 3’ UTRs than in the whole genome. We found that the 3’ UTRs in cichlids tend to be longer than those in non-cichlids, and this was associated, on average, with one more miRNA target per gene in cichlids. Moreover, we provided evidence that 3’ UTRs on average have evolved faster in cichlids than in non-cichlids. Finally, analyses of gene function suggested that both the top 5% longest and 5% most rapidly evolving 3’ UTRs in cichlids tended to be involved in ribosome-associated pathways and translation.

Conclusions

Our results reveal novel patterns of evolution in the 3’ UTRs of teleosts in general and cichlids in particular. The data suggest that 3’ UTRs might serve as important meta-regulators, regulators of other mechanisms governing post-transcriptional regulation, especially in groups like cichlids that have undergone extremely fast rates of phenotypic diversification and speciation.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4821-8) contains supplementary material, which is available to authorized users.

Gene-gene and gene-environment interactions in complex traits in yeast

One of the fundamental questions in biology is how the genotype regulates the phenotype. An increasing number of studies indicate that, in most cases, the effect of a genetic locus on the phenotype is context-dependent, i.e. it is influenced by the genetic background and the environment in which the phenotype is measured. Still, the majority of the studies, in both model organisms and humans, that map the genetic regulation of phenotypic variation in complex traits primarily identify additive loci with independent effects. This does not reflect an absence of the contribution of genetic interactions to phenotypic variation, but instead is a consequence of the technical limitations in mapping gene-gene interactions (GGI) and gene-environment interactions (GEI). Yeast, with its detailed molecular understanding, diverse population genomics and ease of genetic manipulation, is a unique and powerful resource to study the contributions of GGI and GEI in the regulation of phenotypic variation. Here we review studies in yeast that have identified GGI and GEI that regulate phenotypic variation, and discuss the contribution of these findings in explaining missing heritability of complex traits, and how observations from these GGI and GEI studies enhance our understanding of the mechanisms underlying genetic robustness and adaptability that shape the architecture of the genotype-phenotype map.

Bayesian Inference of Allele-Specific Gene Expression Indicates Abundant Cis-Regulatory Variation in Natural Flycatcher Populations

Polymorphism in cis-regulatory sequences can lead to different levels of expression for the two alleles of a gene, providing a starting point for the evolution of gene expression. Little is known about the genome-wide abundance of genetic variation in gene regulation in natural populations but analysis of allele-specific expression (ASE) provides a means for investigating such variation. We performed RNA-seq of multiple tissues from population samples of two closely related flycatcher species and developed a Bayesian algorithm that maximizes data usage by borrowing information from the whole data set and combines several SNPs per transcript to detect ASE. Of 2,576 transcripts analyzed in collared flycatcher, ASE was detected in 185 (7.2%) and a similar frequency was seen in the pied flycatcher. Transcripts with statistically significant ASE commonly showed the major allele in >90% of the reads, reflecting that power was highest when expression was heavily biased toward one of the alleles. This would suggest that the observed frequencies of ASE likely are underestimates. The proportion of ASE transcripts varied among tissues, being lowest in testis and highest in muscle. Individuals often showed ASE of particular transcripts in more than one tissue (73.4%), consistent with a genetic basis for regulation of gene expression. The results suggest that genetic variation in regulatory sequences commonly affects gene expression in natural populations and that it provides a seedbed for phenotypic evolution via divergence in gene expression.

Mating of natural Saccharomyces cerevisiae strains for improved glucose fermentation and lignocellulosic inhibitor tolerance

Natural Saccharomyces cerevisiae isolates from vineyards in the Western Cape, South Africa were evaluated for ethanol production in industrial conditions associated with the production of second-generation biofuels. The strains displayed high phenotypic diversity including the ability to grow at 45 °C and in the presence of 20% (v/v) ethanol, strain YI13. Strains HR4 and YI30 were inhibitor-tolerant under aerobic and oxygen-limited conditions, respectively. Spore-to-spore hybridization generated progeny that displayed heterosis, including increased ethanol productivity and improved growth in the presence of a synthetic inhibitor cocktail. Hybrid strains HR4/YI30#6 and V3/YI30#6 were able to grow at a high salt concentration (2 mol/L NaCl) with V3/YI30#6 also able to grow at a high temperature (45 °C). Strains HR4/YI30#1 and #3 were inhibitor-tolerant, with strain HR4/YI30#3 having similar productivity (0.36 ± 0.0036 g/L per h) as the superior parental strain, YI30 (0.35 ± 0.0058 g/L per h). This study indicates that natural S. cerevisiae strains display phenotypic variation and heterosis can be achieved through spore-to-spore hybridization. Several of the phenotypes (temperature-, osmo-, and inhibitor tolerance) displayed by both the natural strains and the generated progeny were at the maximum conditions reported for S. cerevisiae strains.

Genetic basis of mycotoxin susceptibility differences between budding yeast isolates

Micophenolic acid (MPA) is an immunosuppressant mycotoxin which impairs yeast cell growth to variable degrees depending on the genetic background. Such variation could have emerged from several phenomena, including MPA gene resistance mutations and variations in copy number and localisation of resistance genes. To test this, we evaluated MPA susceptibility in four S. cerevisiae isolates and genetically dissected variation through the identification of Quantitative Trait Loci. Via linkage analysis we identified six QTLs, majority of which were located within subtelomeres and co-localised with IMD2, an inosine monophosphate dehydrogenase previously identified underlying MPA drug resistance in yeast cells. From chromosome end disruption and bioinformatics analysis, it was found that the subtelomere localisation of IMD2 within chromosome ends is variable depending on the strain, demonstrating the influence of IMD2 on the natural variation in yeast MPA susceptibility. Furthermore, GxE gene expression analysis of strains exhibiting opposite phenotypes indicated that ribosome biogenesis, RNA transport, and purine biosynthesis were impaired in strains most susceptible to MPA toxicity. Our results demonstrate that natural variation can be exploited to better understand the molecular mechanisms underlying mycotoxin susceptibility in eukaryote cells and demonstrate the role of subtelomeric regions in mediating interactions with the environment.

Identification of Nitrogen Consumption Genetic Variants in Yeast Through QTL Mapping and Bulk Segregant RNA-Seq Analyses

Saccharomyces cerevisiae is responsible for wine must fermentation. In this process, nitrogen represents a limiting nutrient and its scarcity results in important economic losses for the wine industry. Yeast isolates use different strategies to grow in poor nitrogen environments and their genomic plasticity enables adaptation to multiple habitats through improvements in nitrogen consumption. Here, we used a highly recombinant S. cerevisiae multi-parent population (SGRP-4X) derived from the intercross of four parental strains of different origins to identify new genetic variants responsible for nitrogen consumption differences during wine fermentation. Analysis of 165 fully sequenced F12 segregants allowed us to map 26 QTL in narrow intervals for 14 amino acid sources and ammonium, the majority of which represent genomic regions previously unmapped for these traits. To complement this strategy, we performed Bulk segregant RNA-seq (BSR-seq) analysis in segregants exhibiting extremely high and low ammonium consumption levels. This identified several QTL overlapping differentially expressed genes and refined the gene candidate search. Based on these approaches, we were able to validate ARO1, PDC1, CPS1, ASI2, LYP1, and ALP1 allelic variants underlying nitrogen consumption differences between strains, providing evidence of many genes with small phenotypic effects. Altogether, these variants significantly shape yeast nitrogen consumption with important implications for evolution, ecological, and quantitative genomics.

Adaptive divergence in wine yeasts and their wild relatives suggests a prominent role for introgressions and rapid evolution at noncoding sites

In Saccharomyces cerevisiae, the main yeast in wine fermentation, the opportunity to examine divergence at the molecular level between a domesticated lineage and its wild counterpart arose recently due to the identification of the closest relatives of wine strains, a wild population associated with Mediterranean oaks. As genomic data are available for a considerable number of representatives belonging to both groups, we used population genomics to estimate the degree and distribution of nucleotide variation between wine yeasts and their closest wild relatives. We found widespread genomewide divergence, particularly at noncoding sites, which, together with above average divergence in trans-acting DNA binding proteins, may suggest an important role for divergence at the level of transcriptional regulation. Nine outlier regions putatively under strong divergent selection were highlighted by a genomewide scan under stringent conditions. Several cases of introgressions, originating in the sibling species Saccharomyces paradoxus, were also identified in the Mediterranean oak population. FZF1 and SSU1, mostly known for conferring sulphite resistance in wine yeasts, were among the introgressed genes, although not fixed. Because the introgressions detected in our study are not found in wine strains, we hypothesize that ongoing divergent ecological selection segregates the two forms between the different niches. Together, our results provide a first insight into the extent and kind of divergence between wine yeasts and their closest wild relatives.

Natural gene expression variation studies in yeast

The rise of sequence information across different yeast species and strains is driving an increasing number of studies in the emerging field of genomics to associate polymorphic variants, mRNA abundance and phenotypic differences between individuals. Here, we gathered evidence from recent studies covering several layers that define the genotype-phenotype gap, such as mRNA abundance, allele-specific expression and translation efficiency to demonstrate how genetic variants co-evolve and define an individual's genome. Moreover, we exposed several antecedents where inter- and intra-specific studies led to opposite conclusions, probably owing to genetic divergence. Future studies in this area will benefit from the access to a massive array of well-annotated genomes and new sequencing technologies, which will allow the fine breakdown of the complex layers that delineate the genotype-phenotype map. Copyright © 2016 John Wiley & Sons, Ltd.

Sexual selection in hermaphrodites, sperm and broadcast spawners, plants and fungi

Darwin was the first to recognize that sexual selection is a strong evolutionary force. Exaggerated traits allow same-sex individuals to compete over access to mates and provide a mechanism by which mates are selected. It is relatively easy to appreciate how inter- and intrasexual selection work in organisms with the sensory capabilities to perceive physical or behavioural traits that signal mate quality or mate compatibility, and to assess the relative quality of competitors. It is therefore not surprising that most studies of sexual selection have focused on animals with separate sexes and obvious adaptations that function in the context of reproductive competition. Yet, many sexual organisms are both male and female at the same time, often lack sexual dimorphism and never come into direct contact at mating. How does sexual selection act in such species, and what can we learn from them? Here, we address these questions by exploring the potential for sexual selection in simultaneous hermaphrodites, sperm- and broadcast spawners, plants and fungi. Our review reveals a range of mechanisms of sexual selection, operating primarily after gametes have been released, which are common in many of these groups and also quite possibly in more familiar (internally fertilizing and sexually dimorphic) organisms.This article is part of the themed issue 'Weird sex: the underappreciated diversity of sexual reproduction'.

Exploiting budding yeast natural variation for industrial processes

For the last two decades, the natural variation of the yeast Saccharomyces cerevisiae has been massively exploited with the aim of understanding ecological and evolutionary processes. As a result, many new genetic variants have been uncovered, providing a large catalogue of alleles underlying complex traits. These alleles represent a rich genetic resource with the potential to provide new strains that can cope with the growing demands of industrial fermentation processes. When surveyed in detail, several of these variants have proven useful in wine and beer industries by improving nitrogen utilisation, fermentation kinetics, ethanol production, sulphite resistance and aroma production. Here, I illustrate how allele-specific expression and polymorphisms within the coding region of GDB1 underlie fermentation kinetic differences in synthetic wine must. Nevertheless, the genetic basis of how GDB1 variants and other natural alleles interact in foreign genetic backgrounds remains unclear. Further studies in large sets of strains, recombinant hybrids and multiple parental pairs will broaden our knowledge of the molecular and genetic basis of trait adaptation for utilisation in applied and industrial processes.