The evolution, evolvability and engineering of gene regulatory DNA.

Mutations in non-coding regulatory DNA sequences can alter gene expression, organismal phenotype and fitness1-3. Constructing complete fitness landscapes, in which DNA sequences are mapped to fitness, is a long-standing goal in biology, but has remained elusive because it is challenging to generalize reliably to vast sequence spaces4-6. Here we build sequence-to-expression models that capture fitness landscapes and use them to decipher principles of regulatory evolution. Using millions of randomly sampled promoter DNA sequences and their measured expression levels in the yeast Saccharomyces cerevisiae, we learn deep neural network models that generalize with excellent prediction performance, and enable sequence design for expression engineering. Using our models, we study expression divergence under genetic drift and strong-selection weak-mutation regimes to find that regulatory evolution is rapid and subject to diminishing returns epistasis; that conflicting expression objectives in different environments constrain expression adaptation; and that stabilizing selection on gene expression leads to the moderation of regulatory complexity. We present an approach for using such models to detect signatures of selection on expression from natural variation in regulatory sequences and use it to discover an instance of convergent regulatory evolution. We assess mutational robustness, finding that regulatory mutation effect sizes follow a power law, characterize regulatory evolvability, visualize promoter fitness landscapes, discover evolvability archetypes and illustrate the mutational robustness of natural regulatory sequence populations. Our work provides a general framework for designing regulatory sequences and addressing fundamental questions in regulatory evolution.

Wild Patagonian yeast improve the evolutionary potential of novel interspecific hybrid strains for lager brewing.

Lager yeasts are limited to a few strains worldwide, imposing restrictions on flavour and aroma diversity and hindering our understanding of the complex evolutionary mechanisms during yeast domestication. The recent finding of diverse S. eubayanus lineages from Patagonia offers potential for generating new lager yeasts with different flavour profiles. Here, we leverage the natural genetic diversity of S. eubayanus and expand the lager yeast repertoire by including three distinct Patagonian S. eubayanus lineages. We used experimental evolution and selection on desirable traits to enhance the fermentation profiles of novel S. cerevisiae x S. eubayanus hybrids. Our analyses reveal an intricate interplay of pre-existing diversity, selection on species-specific mitochondria, de-novo mutations, and gene copy variations in sugar metabolism genes, resulting in high ethanol production and unique aroma profiles. Hybrids with S. eubayanus mitochondria exhibited greater evolutionary potential and superior fitness post-evolution, analogous to commercial lager hybrids. Using genome-wide screens of the parental subgenomes, we identified genetic changes in IRA2, IMA1, and MALX genes that influence maltose metabolism, and increase glycolytic flux and sugar consumption in the evolved hybrids. Functional validation and transcriptome analyses confirmed increased maltose-related gene expression, influencing greater maltotriose consumption in evolved hybrids. This study demonstrates the potential for generating industrially viable lager yeast hybrids from wild Patagonian strains. Our hybridization, evolution, and mitochondrial selection approach produced hybrids with high fermentation capacity and expands lager beer brewing options.

Climate change and population persistence in a hibernating marsupial.

Climate change has physiological consequences on organisms, ecosystems and human societies, surpassing the pace of organismal adaptation. Hibernating mammals are particularly vulnerable as winter survival is determined by short-term physiological changes triggered by temperature. In these animals, winter temperatures cannot surpass a certain threshold, above which hibernators arouse from torpor, increasing several fold their energy needs when food is unavailable. Here, we parameterized a numerical model predicting energy consumption in heterothermic species and modelled winter survival at different climate change scenarios. As a model species, we used the arboreal marsupial monito del monte (genus Dromiciops), which is recognized as one of the few South American hibernators. We modelled four climate change scenarios (from optimistic to pessimistic) based on IPCC projections, predicting that northern and coastal populations (Dromiciops bozinovici) will decline because the minimum number of cold days needed to survive the winter will not be attained. These populations are also the most affected by habitat fragmentation and changes in land use. Conversely, Andean and other highland populations, in cooler environments, are predicted to persist and thrive. Given the widespread presence of hibernating mammals around the world, models based on simple physiological parameters, such as this one, are becoming essential for predicting species responses to warming in the short term.

Nakazawaea atacamensis f.a., sp. nov. a novel nonconventional fermentative ascomycetous yeast species from the Atacama Desert.

In this study, we describe Nakazawaea atacamensis f. a., sp. nov., a novel species obtained from Neltuma chilensis plant samples in Chile's hyperarid Atacama Desert. In total, three strains of N. atacamensis were obtained from independent N. chilensis samples (synonym Prosopis chilensis, Algarrobo). Two strains were obtained from bark samples, while the third strain was obtained from bark-exuded gum from another tree. The novel species was defined using molecular characteristics and subsequently characterized with respect to morphological, physiological, and biochemical properties. A neighbor-joining analysis using the sequences of the D1/D2 domains of the large subunit ribosomal RNA gene revealed that N. atacamensis clustered with Nakazawaea pomicola. The sequence of N. atacamensis differed from closely related species by 1.3%-5.2% in the D1/D2 domains. A phylogenomic analysis based on single-nucleotide polymorphism's data confirms that the novel species belongs to the genus Nakazawaea, where N. atacamensis clustered with N. peltata. Phenotypic comparisons demonstrated that N. atacamensis exhibited distinct carbon assimilation patterns compared to its related species. Genome sequencing of the strain ATA-11A-BT revealed a genome size of approximately 12.4 Mbp, similar to other Nakazawaea species, with 5116 protein-coding genes annotated using InterProScan. In addition, N. atacamensis exhibited the capacity to ferment synthetic wine must, representing a potential new yeast for mono or co-culture wine fermentations. This comprehensive study expands our understanding of the genus Nakazawaea and highlights the ecological and industrial potential of N. atacamensis in fermentation processes. The holotype of N. atacamensis sp. nov. is CBS 18375T . The Mycobank number is MB 849680.

Taxogenomic analyses of Starmerella gilliamiae f.a, sp. nov. and Starmerella monicapupoae f.a., sp. nov., two novel species isolated from plant substrates and insects.

Four yeast isolates collected from flowers from different ecosystems in Brazil, one from fruit of Nothofagus alpina in Argentina, three from flowers of Neltuma chilensis in Chile and one obtained from the proventriculus of a female bumblebee in Canada were demonstred, by analysis of the sequences of the internal transcribed spacer (ITS) region and D1/D2 domains of the large subunit rRNA gene, to represent two novel species of the genus Starmerella. These species are described here as Starmerella gilliamiae f.a, sp. nov. (CBS 16166T; Mycobank MB 851206) and Starmerella monicapupoae f.a., sp. nov. (PYCC 8997T; Mycobank MB 851207). The results of a phylogenomic analysis using 1037 single-copy orthogroups indicated that S. gilliamiae is a member of a subclade that contains Starmerella opuntiae, Starmerella aceti and Starmerella apicola. The results also indicated that S. monicapupoae is phylogenetically related to Starmerella riodocensis. The two isolates of S. monicapupoae were obtained from flowers in Brazil and were probably vectored by insects that visit these substrates. Starmerella gilliamiae has a wide geographical distribution having been isolated in flowers from Brazil and Chile, fruit from Argentina and a bumblebee from Canada.

Modeling heterothermic fitness landscapes in a marsupial hibernator using changes in body composition.

Hibernation is an adaptive strategy that allows animals to enter a hypometabolic state, conserving energy and enhancing their fitness by surviving harsh environmental conditions. However, addressing the adaptive value of hibernation, at the individual level and in natural populations, has been challenging. Here, we applied a non-invasive technique, body composition analysis by quantitative magnetic resonance (qMR), to calculate energy savings by hibernation in a population of hibernating marsupials (Dromiciops gliroides). Using outdoor enclosures installed in a temperate rainforest, and measuring qMR periodically, we determined the amount of fat and lean mass consumed during a whole hibernation cycle. With this information, we estimated the daily energy expenditure of hibernation (DEEH) at the individual level and related to previous fat accumulation. Using model selection approaches and phenotypic selection analysis, we calculated linear (directional, ß), quadratic (stabilizing or disruptive, γ) and correlational (ρ) coefficients for DEEH and fat accumulation. We found significant, negative directional selection for DEEH (ßDEEH = - 0.58 ± 0.09), a positive value for fat accumulation (ßFAT = 0.34 ± 0.07), and positive correlational selection between both traits (ρDEEH × FAT = 0.24 ± 0.07). Then, individuals maximizing previous fat accumulation and minimizing DEEH were promoted by selection, which is visualized by a bi-variate selection surface estimated by generalized additive models. At the comparative level, results fall within the isometric allometry known for hibernation metabolic rate in mammals. Thus, by a combination of a non-invasive technique for body composition analysis and semi-natural enclosures, we were characterized the heterothermic fitness landscape in a semi-natural population of hibernators.

Addition of Saccharomyces eubayanus to SCOBY fermentations modulates the chemical and volatile compound profiles in kombucha.

Kombucha is a fermented beverage derived from a sweetened tea fermentation inoculated with a bacteria-yeast consortium referred to as Symbiotic Culture of Bacteria and Yeast (SCOBY). Different SCOBY cultures can impact the beverage's quality and make the whole process highly variable. Adding Saccharomyces yeast cultures to the fermentation process can avoid stalled fermentations, providing a reproducible beverage. Here, we explored using different Saccharomyces eubayanus strains together with SCOBY in the context of kombucha fermentation. Our results show that yeast x SCOBY co-cultures exhibited a robust fermentation profile, providing ethanol and acetic acid levels ranging from 0,18-1,81 %v/v and 0,35-1,15 g/L, respectively. The kombucha volatile compound profile of co-cultures was unique, where compounds such as Isopentyl acetate where only found in yeast x SCOBY fermentations. Metabarcoding revealed that the SCOBY composition was also dependent on the S. eubayanus genotype, where besides Saccharomyces, amplicon sequence variants belonging to Brettanomyces and Starmerella were detected. These differences concomitated global changes in transcript levels in S. eubayanus related to the metabolism of organic molecules used in kombucha fermentation. This study highlights the potential for exploring different S. eubayanus strains for kombucha fermentation, and the significant yeast genotype effect in the profile differentiation in this process.

An Out-of-Patagonia migration explains the worldwide diversity and distribution of Saccharomyces eubayanus lineages.

Population-level sampling and whole-genome sequences of different individuals allow one to identify signatures of hybridization, gene flow and potential molecular mechanisms of environmental responses. Here, we report the isolation of 160 Saccharomyces eubayanus strains, the cryotolerant ancestor of lager yeast, from ten sampling sites in Patagonia along 2,000 km of Nothofagus forests. Frequency of S. eubayanus isolates was higher towards southern and colder regions, demonstrating the cryotolerant nature of the species. We sequenced the genome of 82 strains and, together with 23 available genomes, performed a comprehensive phylogenetic analysis. Our results revealed the presence of five different lineages together with dozens of admixed strains. Various analytical methods reveal evidence of gene flow and historical admixture between lineages from Patagonia and Holarctic regions, suggesting the co-occurrence of these ancestral populations. Analysis of the genetic contribution to the admixed genomes revealed a Patagonian genetic origin of the admixed strains, even for those located in the North Hemisphere. Overall, the Patagonian lineages, particularly the southern populations, showed a greater global genetic diversity compared to Holarctic and Chinese lineages, in agreement with a higher abundance in Patagonia. Thus, our results are consistent with a likely colonization of the species from peripheral glacial refugia from South Patagonia. Furthermore, fermentative capacity and maltose consumption resulted negatively correlated with latitude, indicating better fermentative performance in northern populations. Our genome analysis, together with previous reports in the sister species S. uvarum suggests that a S. eubayanus ancestor was adapted to the harsh environmental conditions of Patagonia, a region that provides the ecological conditions for the diversification of these ancestral lineages.

Late Pleistocene-dated divergence between South Hemisphere populations of the non-conventional yeast L. cidri.

Most organisms belonging to the Saccharomycotina subphylum have high genetic diversity and a vast repertoire of metabolisms and lifestyles. Lachancea cidri is an ideal yeast model for exploring the interplay between genetics, ecological function and evolution. Lachancea cidri diverged from the Saccharomyces lineage before the whole-genome duplication and is distributed across the South Hemisphere, displaying an important ecological success. We applied phylogenomics to investigate the genetic variation of L. cidri isolates obtained from Australia and South America. Our approach revealed the presence of two main lineages according to their geographic distribution (Aus and SoAm). Estimation of the divergence time suggests that SoAm and Aus lineages diverged near the last glacial maximum event during the Pleistocene (64-8 KYA). Interestingly, we found that the French reference strain is closely related to the Australian strains, with a recent divergence (405-51 YA), likely associated to human movements. Additionally, we identified different lineages within the South American population, revealing that Patagonia contains a similar genetic diversity comparable to that of other lineages in S. cerevisiae. These findings support the idea of a Pleistocene-dated divergence between South Hemisphere lineages, where the Nothofagus and Araucaria ecological niches likely favoured the extensive distribution of L. cidri in Patagonia.

Identification of new ethanol-tolerant yeast strains with fermentation potential from central Patagonia.

The quest for new wild yeasts has increasingly gained attention because of their potential ability to provide unique organoleptic characters to fermented beverages. In this sense, Patagonia offers a wide diversity of ethanol-tolerant yeasts and stands out as a bioprospecting alternative. This study characterized the genetic and phenotypic diversity of yeast isolates obtained from Central Chilean Patagonia and analyzed their fermentation potential under different fermentative conditions. We recovered 125 colonies from Nothofagus spp. bark samples belonging to five yeast species: Saccharomyces eubayanus, Saccharomyces uvarum, Lachancea cidri, Kregervanrija delftensis, and Hanseniaspora valbyensis. High-throughput microcultivation assays demonstrated the extensive phenotypic diversity among Patagonian isolates, where Saccharomyces spp and L. cidri isolates exhibited the most outstanding fitness scores across the conditions tested. Fermentation performance assays under wine, mead, and beer conditions demonstrated the specific potential of the different species for each particular beverage. Saccharomyces spp. were the only isolates able to ferment beer wort. Interestingly, we found that L. cidri is a novel candidate species to ferment wine and mead, exceeding the fermentation capacity of a commercial strain. Unlike commercial strains, we found that L. cidri does not require nutritional supplements for efficient mead fermentation. In addition, L. cidri produces succinic and acetic acids, providing a distinct profile to the final fermented product. This work demonstrates the importance of bioprospecting efforts in Patagonia to isolate novel wild yeast strains with extraordinary biotechnological potential for the fermentation industry.

Yeasts from temperate forests.

Yeasts are ubiquitous in temperate forests. While this broad habitat is well-defined, the yeasts inhabiting it and their life cycles, niches, and contributions to ecosystem functioning are less understood. Yeasts are present on nearly all sampled substrates in temperate forests worldwide. They associate with soils, macroorganisms, and other habitats and no doubt contribute to broader ecosystem-wide processes. Researchers have gathered information leading to hypotheses about yeasts' niches and their life cycles based on physiological observations in the laboratory as well as genomic analyses, but the challenge remains to test these hypotheses in the forests themselves. Here, we summarize the habitat and global patterns of yeast diversity, give some information on a handful of well-studied temperate forest yeast genera, discuss the various strategies to isolate forest yeasts, and explain temperate forest yeasts' contributions to biotechnology. We close with a summary of the many future directions and outstanding questions facing researchers in temperate forest yeast ecology. Yeasts present an exciting opportunity to better understand the hidden world of microbial ecology in this threatened and global habitat.

Genetic Background Matters: Population-Based Studies in Model Organisms for Translational Research.

We are all similar but a bit different. These differences are partially due to variations in our genomes and are related to the heterogeneity of symptoms and responses to treatments that patients exhibit. Most animal studies are performed in one single strain with one manipulation. However, due to the lack of variability, therapies are not always reproducible when treatments are translated to humans. Panels of already sequenced organisms are valuable tools for mimicking human phenotypic heterogeneities and gene mapping. This review summarizes the current knowledge of mouse, fly, and yeast panels with insightful applications for translational research.

Anthocyanins from Aristotelia chilensis Prevent Olanzapine-Induced Hepatic-Lipid Accumulation but Not Insulin Resistance in Skeletal Muscle Cells.

Type 2 diabetes and obesity are major problems worldwide and dietary polyphenols have shown efficacy to ameliorate signs of these diseases. Anthocyanins from berries display potent antioxidants and protect against weight gain and insulin resistance in different models of diet-induced metabolic syndrome. Olanzapine is known to induce an accelerated form of metabolic syndrome. Due to the aforementioned, we evaluated whether delphinidin-3,5-O-diglucoside (DG) and delphinidin-3-O-sambubioside-5-O-glucoside (DS), two potent antidiabetic anthocyanins isolated from Aristotelia chilensis fruit, could prevent olanzapine-induced steatosis and insulin resistance in liver and skeletal muscle cells, respectively. HepG2 liver cells and L6 skeletal muscle cells were co-incubated with DG 50 µg/mL or DS 50 µg/mL plus olanzapine 50 µg/mL. Lipid accumulation was determined in HepG2 cells while the expression of p-Akt as a key regulator of the insulin-activated signaling pathways, mitochondrial function, and glucose uptake was assessed in L6 cells. DS and DG prevented olanzapine-induced lipid accumulation in liver cells. However, insulin signaling impairment induced by olanzapine in L6 cells was not rescued by DS and DG. Thus, anthocyanins modulate lipid metabolism, which is a relevant factor in hepatic tissue, but do not significantly influence skeletal muscle, where a potent antioxidant effect of olanzapine was found.

Bioprospecting for brewers: Exploiting natural diversity for naturally diverse beers.

The burgeoning interest in archaic, traditional, and novel beer styles has coincided with a growing appreciation of the role of yeasts in determining beer character as well as a better understanding of the ecology and biogeography of yeasts. Multiple studies in recent years have highlighted the potential of wild Saccharomyces and non-Saccharomyces yeasts for production of beers with novel flavour profiles and other desirable properties. Yeasts isolated from spontaneously fermented beers as well as from other food systems (wine, bread, and kombucha) have shown promise for brewing application, and there is evidence that such cross-system transfers have occurred naturally in the past. We review here the available literature pertaining to the use of nonconventional yeasts in brewing, with a focus on the origins of these yeasts, including methods of isolation. Practical aspects of utilizing nondomesticated yeasts are discussed, and modern methods to facilitate discovery of yeasts with brewing potential are highlighted.

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.

Extensive cis-regulatory variation robust to environmental perturbation in Arabidopsis.

cis- and trans-acting factors affect gene expression and responses to environmental conditions. However, for most plant systems, we lack a comprehensive map of these factors and their interaction with environmental variation. Here, we examined allele-specific expression (ASE) in an F1 hybrid to study how alleles from two Arabidopsis thaliana accessions affect gene expression. To investigate the effect of the environment, we used drought stress and developed a variance component model to estimate the combined genetic contributions of cis- and trans-regulatory polymorphisms, environmental factors, and their interactions. We quantified ASE for 11,003 genes, identifying 3318 genes with consistent ASE in control and stress conditions, demonstrating that cis-acting genetic effects are essentially robust to changes in the environment. Moreover, we found 1618 genes with genotype x environment (GxE) interactions, mostly cis x E interactions with magnitude changes in ASE. We found fewer trans x E interactions, but these effects were relatively less robust across conditions, showing more changes in the direction of the effect between environments; this confirms that trans-regulation plays an important role in the response to environmental conditions. Our data provide a detailed map of cis- and trans-regulation and GxE interactions in A. thaliana, laying the ground for mechanistic investigations and studies in other plants and environments.

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.

RIM15 antagonistic pleiotropy is responsible for differences in fermentation and stress response kinetics in budding yeast.

Different natural yeast populations have faced dissimilar selective pressures due to the heterogeneous fermentation substrates available around the world; this increases the genetic and phenotypic diversity in Saccharomyces cerevisiae In this context, we expect prominent differences between isolates when exposed to a particular condition, such as wine or sake musts. To better comprehend the mechanisms underlying niche adaptation between two S. cerevisiae isolates obtained from wine and sake fermentation processes, we evaluated fermentative and fungicide resistance phenotypes and identify the molecular origin of such adaptive variation. Multiple regions were associated with fermentation rate under different nitrogen conditions and fungicide resistance, with a single QTL co-localizing in all traits. Analysis around this region identified RIM15 as the causative locus driving fungicide sensitivity, together with efficient nitrogen utilization and glycerol production in the wine strain. A null RIM15 variant confers a greater fermentation rate through the utilization of available glucose instead of its storage. However, this variant has a detrimental effect on fungicide resistance since complex sugars are not synthesized and transported into the membrane. Together, our results reveal the antagonist pleiotropic nature of a RIM15 null variant, positively affecting a series of fermentation related phenotypes, but apparently detrimental in the wild.

Revealing the genetic structure of a trait by sequencing a population under selection.

One approach to understanding the genetic basis of traits is to study their pattern of inheritance among offspring of phenotypically different parents. Previously, such analysis has been limited by low mapping resolution, high labor costs, and large sample size requirements for detecting modest effects. Here, we present a novel approach to map trait loci using artificial selection. First, we generated populations of 10-100 million haploid and diploid segregants by crossing two budding yeast strains of different heat tolerance for up to 12 generations. We then subjected these large segregant pools to heat stress for up to 12 d, enriching for beneficial alleles. Finally, we sequenced total DNA from the pools before and during selection to measure the changes in parental allele frequency. We mapped 21 intervals with significant changes in genetic background in response to selection, which is several times more than found with traditional linkage methods. Nine of these regions contained two or fewer genes, yielding much higher resolution than previous genomic linkage studies. Multiple members of the RAS/cAMP signaling pathway were implicated, along with genes previously not annotated with heat stress response function. Surprisingly, at most selected loci, allele frequencies stopped changing before the end of the selection experiment, but alleles did not become fixed. Furthermore, we were able to detect the same set of trait loci in a population of diploid individuals with similar power and resolution, and observed primarily additive effects, similar to what is seen for complex trait genetics in other diploid organisms such as humans.