The genetic architecture of low-temperature adaptation in the wine yeast Saccharomyces cerevisiae

[Background] Low-temperature growth and fermentation of wine yeast can enhance wine aroma and make them highly desirable traits for the industry. Elucidating response to cold in Saccharomyces cerevisiae is, therefore, of paramount importance to select or genetically improve new wine strains. As most enological traits of industrial importance in yeasts, adaptation to low temperature is a polygenic trait regulated by many interacting loci.[Results] In order to unravel the genetic determinants of low-temperature fermentation, we mapped quantitative trait loci (QTLs) by bulk segregant analyses in the F13 offspring of two Saccharomyces cerevisiae industrial strains with divergent performance at low temperature. We detected four genomic regions involved in the adaptation at low temperature, three of them located in the subtelomeric regions (chromosomes XIII, XV and XVI) and one in the chromosome XIV. The QTL analysis revealed that subtelomeric regions play a key role in defining individual variation, which emphasizes the importance of these regions’ adaptive nature.[Conclusions] The reciprocal hemizygosity analysis (RHA), run to validate the genes involved in low-temperature fermentation, showed that genetic variation in mitochondrial proteins, maintenance of correct asymmetry and distribution of phospholipid in the plasma membrane are key determinants of low-temperature adaptation.This work has been financially supported from the Spanish Government through MINECO and FEDER funds (AGL2013-47300-C3-3-R and PCIN-2015-143 grants) and from Generalitat Valenciana through PROMETEOII/2014/042 grant, awarded to JMG. This study has been carried out in the context of the European Project ERA-IB “YeastTempTation” EGR thanks the Spanish government for an FPI grant BES-2011-044498 and MM also thanks the Generalitat Valenciana for a VALi+d ACIF/2015/194 grant. We acknowledge support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).Peer reviewe

Population genomics of domestic and wild yeasts

The natural genetics of an organism is determined by the distribution of sequences of its genome. Here we present one- to four-fold, with some deeper, coverage of the genome sequences of over seventy isolates of the domesticated baker's yeast, _Saccharomyces cerevisiae_, and its closest relative, the wild _S. paradoxus_, which has never been associated with human activity. These were collected from numerous geographic locations and sources (including wild, clinical, baking, wine, laboratory and food spoilage). These sequences provide an unprecedented view of the population structure, natural (and artificial) selection and genome evolution in these species. Variation in gene content, SNPs, indels, copy numbers and transposable elements provide insights into the evolution of different lineages. Phenotypic variation broadly correlates with global genome-wide phylogenetic relationships however there is no correlation with source. _S. paradoxus_ populations are well delineated along geographic boundaries while the variation among worldwide _S. cerevisiae_ isolates show less differentiation and is comparable to a single _S. paradoxus_ population. Rather than one or two domestication events leading to the extant baker's yeasts, the population structure of _S. cerevisiae_ shows a few well defined geographically isolated lineages and many different mosaics of these lineages, supporting the notion that human influence provided the opportunity for outbreeding and production of new combinations of pre-existing variation

The ribosomal P-stalk couples amino acid starvation to GCN2 activation in mammalian cells

The eukaryotic translation initiation factor 2α (eIF2α) kinase GCN2 is activated by amino acid starvation to elicit a rectifying physiological program known as the Integrated Stress Response (ISR). A role for uncharged tRNAs as activating ligands of yeast GCN2 is supported experimentally. However, mouse GCN2 activation has recently been observed in circumstances associated with ribosome stalling with no global increase in uncharged tRNAs. We report on a mammalian CHO cell-based CRISPR-Cas9 mutagenesis screen for genes that contribute to ISR activation by amino acid starvation. Disruption of genes encoding components of the ribosome P-stalk, uL10 and P1, selectively attenuated GCN2-mediated ISR activation by amino acid starvation or interference with tRNA charging without affecting the endoplasmic reticulum unfolded protein stress-induced ISR, mediated by the related eIF2α kinase PERK. Wildtype ribosomes isolated from CHO cells, but not those with P-stalk lesions, stimulated GCN2-dependent eIF2α phosphorylation in vitro. These observations support a model whereby lack of a cognate charged tRNA exposes a latent capacity of the ribosome P-stalk to activate GCN2 in cells and help explain the emerging link between ribosome stalling and ISR activation

Type II and type IV toxin–antitoxin systems show different evolutionary patterns in the global Klebsiella pneumoniae population

The Klebsiella pneumoniae species complex includes important opportunistic pathogens which have become public health priorities linked to major hospital outbreaks and the recent emergence of multidrug-resistant hypervirulent strains. Bacterial virulence and the spread of multidrug resistance have previously been linked to toxin-antitoxin (TA) systems. TA systems encode a toxin that disrupts essential cellular processes, and a cognate antitoxin which counteracts this activity. Whilst associated with the maintenance of plasmids, they also act in bacterial immunity and antibiotic tolerance. However, the evolutionary dynamics and distribution of TA systems in clinical pathogens are not well understood. Here we present a comprehensive survey and description of the diversity of TA systems in 259 clinically relevant genomes of K. pneumoniae. We show that TA systems are highly prevalent with a median of 20 loci per strain. Importantly, these toxins differ substantially in their distribution patterns and in their range of cognate antitoxins. Classification along these properties suggests different roles of TA systems and highlights the association and co-evolution of toxins and antitoxins

SLING: a tool to search for linked genes in bacterial datasets.

Gene arrays and operons that encode functionally linked proteins form the most basic unit of transcriptional regulation in bacteria. Rules that govern the order and orientation of genes in these systems have been defined; however, these were based on a small set of genomes that may not be representative. The growing availability of large genomic datasets presents an opportunity to test these rules, to define the full range and diversity of these systems, and to understand their evolution. Here we present SLING, a tool to Search for LINked Genes by searching for a single functionally essential gene, along with its neighbours in a rule-defined proximity (https://github.com/ghoresh11/sling/wiki). Examining this subset of genes enables us to understand the basic diversity of these genetic systems in large datasets. We demonstrate the utility of SLING on a clinical collection of enteropathogenic Escherichia coli for two relevant operons: toxin antitoxin (TA) systems and RND efflux pumps. By examining the diversity of these systems, we gain insight on distinct classes of operons which present variable levels of prevalence and ability to be lost or gained. The importance of this analysis is not limited to TA systems and RND pumps, and can be expanded to understand the diversity of many other relevant gene arrays

The miRNA Profile of Human Pancreatic Islets and Beta-Cells and Relationship to Type 2 Diabetes Pathogenesis

Recent advances in the understanding of the genetics of type 2 diabetes (T2D) susceptibility have focused attention on the regulation of transcriptional activity within the pancreatic beta-cell. MicroRNAs (miRNAs) represent an important component of regulatory control, and have proven roles in the development of human disease and control of glucose homeostasis. We set out to establish the miRNA profile of human pancreatic islets and of enriched beta-cell populations, and to explore their potential involvement in T2D susceptibility. We used Illumina small RNA sequencing to profile the miRNA fraction in three preparations each of primary human islets and of enriched beta-cells generated by fluorescence-activated cell sorting. In total, 366 miRNAs were found to be expressed (i.e. >100 cumulative reads) in islets and 346 in beta-cells; of the total of 384 unique miRNAs, 328 were shared. A comparison of the islet-cell miRNA profile with those of 15 other human tissues identified 40 miRNAs predominantly expressed (i.e. >50% of all reads seen across the tissues) in islets. Several highly-expressed islet miRNAs, such as miR-375, have established roles in the regulation of islet function, but others (e.g. miR-27b-3p, miR-192-5p) have not previously been described in the context of islet biology. As a first step towards exploring the role of islet-expressed miRNAs and their predicted mRNA targets in T2D pathogenesis, we looked at published T2D association signals across these sites. We found evidence that predicted mRNA targets of islet-expressed miRNAs were globally enriched for signals of T2D association (p-values <0.01, q-values <0.1). At six loci with genome-wide evidence for T2D association (AP3S2, KCNK16, NOTCH2, SCL30A8, VPS26A, and WFS1) predicted mRNA target sites for islet-expressed miRNAs overlapped potentially causal variants. In conclusion, we have described the miRNA profile of human islets and beta-cells and provide evidence linking islet miRNAs to T2D pathogenesis

Predicting the mutations generated by repair of Cas9-induced double-strand breaks.

The DNA mutation produced by cellular repair of a CRISPR-Cas9-generated double-strand break determines its phenotypic effect. It is known that the mutational outcomes are not random, but depend on DNA sequence at the targeted location. Here we systematically study the influence of flanking DNA sequence on repair outcome by measuring the edits generated by >40,000 guide RNAs (gRNAs) in synthetic constructs. We performed the experiments in a range of genetic backgrounds and using alternative CRISPR-Cas9 reagents. In total, we gathered data for >109 mutational outcomes. The majority of reproducible mutations are insertions of a single base, short deletions or longer microhomology-mediated deletions. Each gRNA has an individual cell-line-dependent bias toward particular outcomes. We uncover sequence determinants of the mutations produced and use these to derive a predictor of Cas9 editing outcomes. Improved understanding of sequence repair will allow better design of gene editing experiments

New genetic loci link adipose and insulin biology to body fat distribution.

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, here we conduct genome-wide association meta-analyses of traits related to waist and hip circumferences in up to 224,459 individuals. We identify 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (BMI), and an additional 19 loci newly associated with related waist and hip circumference measures (P < 5 × 10(-8)). In total, 20 of the 49 waist-to-hip ratio adjusted for BMI loci show significant sexual dimorphism, 19 of which display a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms

Quantitative CRISPR interference screens in yeast identify chemical-genetic interactions and new rules for guide RNA design

BACKGROUND: Genome-scale CRISPR interference (CRISPRi) has been used in human cell lines; however, the features of effective guide RNAs (gRNAs) in different organisms have not been well characterized. Here, we define rules that determine gRNA effectiveness for transcriptional repression in Saccharomyces cerevisiae. RESULTS: We create an inducible single plasmid CRISPRi system for gene repression in yeast, and use it to analyze fitness effects of gRNAs under 18 small molecule treatments. Our approach correctly identifies previously described chemical-genetic interactions, as well as a new mechanism of suppressing fluconazole toxicity by repression of the ERG25 gene. Assessment of multiple target loci across treatments using gRNA libraries allows us to determine generalizable features associated with gRNA efficacy. Guides that target regions with low nucleosome occupancy and high chromatin accessibility are clearly more effective. We also find that the best region to target gRNAs is between the transcription start site (TSS) and 200 bp upstream of the TSS. Finally, unlike nuclease-proficient Cas9 in human cells, the specificity of truncated gRNAs (18 nt of complementarity to the target) is not clearly superior to full-length gRNAs (20 nt of complementarity), as truncated gRNAs are generally less potent against both mismatched and perfectly matched targets. CONCLUSIONS: Our results establish a powerful functional and chemical genomics screening method and provide guidelines for designing effective gRNAs, which consider chromatin state and position relative to the target gene TSS. These findings will enable effective library design and genome-wide programmable gene repression in many genetic backgrounds. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s13059-016-0900-9) contains supplementary material, which is available to authorized users