Systematic identification of gene annotation errors in the widely used yeast mutation collections

Control of telomere length in yeast by SUMOylated PCNA and the Elg1 PCNA unloader

Telomeres cap and protect the linear eukaryotic chromosomes. Telomere length is determined by an equilibrium between positive and negative regulators of telomerase activity. A systematic screen for yeast mutants that affect telomere length maintenance in the yeast Saccharomyces cerevisiae revealed that mutations in any of ~500 genes affects telomere length. One of the genes that, when mutated, causes telomere elongation is ELG1, which encodes an unloader of PCNA, the processivity factor for replicative DNA polymerases. PCNA can undergo SUMOylation on two conserved residues, K164 and K127, or ubiquitination at lysine 164. These modifications have already been implicated in genome stability processes. We report that SUMOylated PCNA acts as a signal that positively regulates telomerase activity. We also uncovered physical interactions between Elg1 and the CST (Cdc13-Stn1-Ten) complex and addressed the mechanism by which Elg1 and Stn1 negatively regulates telomere elongation, coordinated by SUMO. We discuss these results with respect to how chromosomal replication and telomere elongation are coordinated.

Copy number variation alters local and global mutational tolerance

Copy number variants (CNVs), duplications and deletions of genomic sequences, contribute to evolutionary adaptation but can also confer deleterious effects and cause disease. Whereas the effects of amplifying individual genes or whole chromosomes (i.e., aneuploidy) have been studied extensively, much less is known about the genetic and functional effects of CNVs of differing sizes and structures. Here, we investigated Saccharomyces cerevisiae (yeast) strains that acquired adaptive CNVs of variable structures and copy numbers following experimental evolution in glutamine-limited chemostats. Although beneficial in the selective environment, CNVs result in decreased fitness compared with the euploid ancestor in rich media. We used transposon mutagenesis to investigate mutational tolerance and genome-wide genetic interactions in CNV strains. We find that CNVs increase mutational target size, confer increased mutational tolerance in amplified essential genes, and result in novel genetic interactions with unlinked genes. We validated a novel genetic interaction between different CNVs and BMH1 that was common to multiple strains. We also analyzed global gene expression and found that transcriptional dosage compensation does not affect most genes amplified by CNVs, although gene-specific transcriptional dosage compensation does occur for ∼12% of amplified genes. Furthermore, we find that CNV strains do not show previously described transcriptional signatures of aneuploidy. Our study reveals the extent to which local and global mutational tolerance is modified by CNVs with implications for genome evolution and CNV-associated diseases, such as cancer.

Global analysis of the yeast knockout phenome

Genome-wide phenotypic screens in the budding yeast Saccharomyces cerevisiae, enabled by its knockout collection, have produced the largest, richest, and most systematic phenotypic description of any organism. However, integrative analyses of this rich data source have been virtually impossible because of the lack of a central data repository and consistent metadata annotations. Here, we describe the aggregation, harmonization, and analysis of ~14,500 yeast knockout screens, which we call Yeast Phenome. Using this unique dataset, we characterized two unknown genes (YHR045W and YGL117W) and showed that tryptophan starvation is a by-product of many chemical treatments. Furthermore, we uncovered an exponential relationship between phenotypic similarity and intergenic distance, which suggests that gene positions in both yeast and human genomes are optimized for function.

Functional Clustering of Metabolically Related Genes Is Conserved across Dikarya

Transcriptional regulation is vital for organismal survival, with many layers and mechanisms collaborating to balance gene expression. One layer of this regulation is genome organization, specifically the clustering of functionally related, co-expressed genes along the chromosomes. Spatial organization allows for position effects to stabilize RNA expression and balance transcription, which can be advantageous for a number of reasons, including reductions in stochastic influences between the gene products. The organization of co-regulated gene families into functional clusters occurs extensively in Ascomycota fungi. However, this is less characterized within the related Basidiomycota fungi despite the many uses and applications for the species within this clade. This review will provide insight into the prevalence, purpose, and significance of the clustering of functionally related genes across Dikarya, including foundational studies from Ascomycetes and the current state of our understanding throughout representative Basidiomycete species.

Transgene insertion into the plastid genome alters expression of adjacent native chloroplast genes at the transcriptional and translational levels

In plant biotechnology and basic research, chloroplasts have been used as chassis for the expression of various transgenes. However, potential unintended side effects of transgene insertion and high-level transgene expression on the expression of native chloroplast genes are often ignored and have not been studied comprehensively. Here, we examined expression of the chloroplast genome at both the transcriptional and translational levels in five transplastomic tobacco (Nicotiana tabacum) lines carrying the identical aadA resistance marker cassette in diverse genomic positions. Although none of the lines exhibits a pronounced visible phenotype, the analysis of three lines that contain the aadA insertion in different locations within the petL-petG-psaJ-rpl33-rps18 transcription unit demonstrates that transcriptional read-through from the aadA resistance marker is unavoidable, and regularly causes overexpression of downstream sense-oriented chloroplast genes at the transcriptional and translational levels. Investigation of additional lines that harbour the aadA intergenically and outside of chloroplast transcription units revealed that expression of the resistance marker can also cause antisense effects by interference with transcription/transcript accumulation and/or translation of downstream antisense-oriented genes. In addition, we provide evidence for a previously suggested role of genomically encoded tRNAs in chloroplast transcription termination and/or transcript processing. Together, our data uncover principles of neighbouring effects of chloroplast transgenes and suggest general strategies for the choice of transgene insertion sites and expression elements to minimize unintended consequences of transgene expression on the transcription and translation of native chloroplast genes.

Bidirectional promoter activity from expression cassettes can drive off-target repression of neighboring gene translation

Targeted selection-based genome-editing approaches have enabled many fundamental discoveries and are used routinely with high precision. We found, however, that replacement of DBP1 with a common selection cassette in budding yeast led to reduced expression and function for the adjacent gene, MRP51, despite all MRP51 coding and regulatory sequences remaining intact. Cassette-induced repression of MRP51 drove all mutant phenotypes detected in cells deleted for DBP1. This behavior resembled the 'neighboring gene effect' (NGE), a phenomenon of unknown mechanism whereby cassette insertion at one locus reduces the expression of a neighboring gene. Here, we leveraged strong off-target mutant phenotypes resulting from cassette replacement of DBP1 to provide mechanistic insight into the NGE. We found that the inherent bidirectionality of promoters, including those in expression cassettes, drives a divergent transcript that represses MRP51 through combined transcriptional interference and translational repression mediated by production of a long undecoded transcript isoform (LUTI). Divergent transcript production driving this off-target effect is general to yeast expression cassettes and occurs ubiquitously with insertion. Despite this, off-target effects are often naturally prevented by local sequence features, such as those that terminate divergent transcripts between the site of cassette insertion and the neighboring gene. Thus, cassette-induced off-target effects can be eliminated by the insertion of transcription terminator sequences into the cassette, flanking the promoter. Because the driving features of this off-target effect are broadly conserved, our study suggests it should be considered in the design and interpretation of experiments using integrated expression cassettes in other eukaryotic systems, including human cells.

A standard knockout procedure alters expression of adjacent loci at the translational level

The Saccharomyces cerevisiae gene deletion collection is widely used for functional gene annotation and genetic interaction analyses. However, the standard G418-resistance cassette used to produce knockout mutants delivers strong regulatory elements into the target genetic loci. To date, its side effects on the expression of neighboring genes have never been systematically assessed. Here, using ribosome profiling data, RT-qPCR, and reporter expression, we investigated perturbations induced by the KanMX module. Our analysis revealed significant alterations in the transcription efficiency of neighboring genes and, more importantly, severe impairment of their mRNA translation, leading to changes in protein abundance. In the ‘head-to-head’ orientation of the deleted and neighboring genes, knockout often led to a shift of the transcription start site of the latter, introducing new uAUG codon(s) into the expanded 5′ untranslated region (5′ UTR). In the ‘tail-to-tail’ arrangement, knockout led to activation of alternative polyadenylation signals in the neighboring gene, thus altering its 3′ UTR. These events may explain the so-called neighboring gene effect (NGE), i.e. false genetic interactions of the deleted genes. We estimate that in as much as ∼1/5 of knockout strains the expression of neighboring genes may be substantially (>2-fold) deregulated at the level of translation.

Overexpression approaches to advance understanding of Candida albicans

Candida albicans is an opportunistic fungal pathogen that is responsible for infections linked to high mortality. Loss‐of‐function approaches, taking advantage of gene knockouts or inducible down‐regulation, have been successfully used in this species in order to understand gene function. However, overexpression of a gene provides an alternative, powerful tool to elucidate gene function and identify novel phenotypes. Notably, overexpression can identify pathway components that might remain undetected using loss‐of‐function approaches. Several repressible or inducible promoters have been developed which allow to shut off or turn on the expression of a gene in C. albicans upon growth in the presence of a repressor or inducer. In this review, we summarize recent overexpression approaches used to study different aspects of C. albicans biology, including morphogenesis, biofilm formation, drug tolerance, and commensalism.

In vivo imaging of fluorescent single-walled carbon nanotubes within C. elegans nematodes in the near-infrared window

Caenorhabditis elegans (C. elegans) nematodes serve as a model organism for eukaryotes, especially due to their genetic similarity. Although they have many advantages like their small size and transparency, their autofluorescence in the entire visible wavelength range poses a challenge for imaging and tracking fluorescent proteins or dyes using standard fluorescence microscopy. Herein, near-infrared (NIR) fluorescent single-walled carbon nanotubes (SWCNTs) are utilized for in vivo imaging within the gastrointestinal track of C. elegans. The SWCNTs are biocompatible, and do not affect the worms' viability nor their reproduction ability. The worms do not show any autofluorescence in the NIR range, thus enabling the spectral separation between the SWCNT NIR fluorescence and the strong autofluorescence of the worm gut granules. The worms are fed with ssDNA-SWCNT which are visualized mainly in the intestine lumen. The NIR fluorescence is used in vivo to track the contraction and relaxation in the area of the pharyngeal valve at the anterior of the terminal bulb. These biocompatible, non-photobleaching, NIR fluorescent nanoparticles can advance in vivo imaging and tracking within C. elegans and other small model organisms by overcoming the signal-to-noise challenge stemming from the wide-range visible autofluorescence.

Transcription at a Distance in the Budding Yeast Saccharomyces cerevisiae

Proper transcriptional regulation depends on the collaboration of multiple layers of control simultaneously. Cells tightly balance cellular resources and integrate various signaling inputs to maintain homeostasis during growth, development and stressors, among other signals. Many eukaryotes, including the budding yeast Saccharomyces cerevisiae, exhibit a non-random distribution of functionally related genes throughout their genomes. This arrangement coordinates the transcription of genes that are found in clusters, and can occur over long distances. In this work, we review the current literature pertaining to gene regulation at a distance in budding yeast.

A case of convergent-gene interference in the budding yeast knockout library causing chromosome instability

To maintain genome stability, organisms depend on faithful chromosome segregation, a process affected by diverse genetic pathways, some of which are not directly linked to mitosis. In this study, we set out to explore one such pathway represented by an undercharacterized gene, SNO1, identified previously in screens of the yeast knockout (YKO) library for mitotic fidelity genes. We found that the causative factor increasing mitotic error rate in the sno1Δ mutant is not loss of the Sno1 protein, but rather perturbation to the mRNA of the neighboring convergent gene, CTF13, encoding an essential component for forming the yeast kinetochore. This is caused by a combination of the Kanamycin resistance gene and the transcriptional terminator used in the YKO library affecting the CTF13 mRNA level and quality . We further provide a list of gene pairs potentially subjected to this artifact, which may be useful for accurate phenotypic interpretation of YKO mutants.

Genomic clustering within functionally related gene families in Ascomycota fungi

Multiple mechanisms collaborate for proper regulation of gene expression. One layer of this regulation is through the clustering of functionally related genes at discrete loci throughout the genome. This phenomenon occurs extensively throughout Ascomycota fungi and is an organizing principle for many gene families whose proteins participate in diverse molecular functions throughout the cell. Members of this phylum include organisms that serve as model systems and those of interest medically, pharmaceutically, and for industrial and biotechnological applications. In this review, we discuss the prevalence of functional clustering through a broad range of organisms within the phylum. Position effects on transcription, genomic locations of clusters, transcriptional regulation of clusters, and selective pressures contributing to the formation and maintenance of clusters are addressed, as are common methods to identify and characterize clusters.

Comparing the utility of in vivo transposon mutagenesis approaches in yeast species to infer gene essentiality

In vivo transposon mutagenesis, coupled with deep sequencing, enables large-scale genome-wide mutant screens for genes essential in different growth conditions. We analyzed six large-scale studies performed on haploid strains of three yeast species (Saccharomyces cerevisiae, Schizosaccaromyces pombe, and Candida albicans), each mutagenized with two of three different heterologous transposons (AcDs, Hermes, and PiggyBac). Using a machine-learning approach, we evaluated the ability of the data to predict gene essentiality. Important data features included sufficient numbers and distribution of independent insertion events. All transposons showed some bias in insertion site preference because of jackpot events, and preferences for specific insertion sequences and short-distance vs long-distance insertions. For PiggyBac, a stringent target sequence limited the ability to predict essentiality in genes with few or no target sequences. The machine learning approach also robustly predicted gene function in less well-studied species by leveraging cross-species orthologs. Finally, comparisons of isogenic diploid versus haploid S. cerevisiae isolates identified several genes that are haplo-insufficient, while most essential genes, as expected, were recessive. We provide recommendations for the choice of transposons and the inference of gene essentiality in genome-wide studies of eukaryotic haploid microbes such as yeasts, including species that have been less amenable to classical genetic studies.

Systematic analysis of nuclear gene function in respiratory growth and expression of the mitochondrial genome in S. cerevisiae

The production of metabolic energy in form of ATP by oxidative phosphorylation depends on the coordinated action of hundreds of nuclear-encoded mitochondrial proteins and a handful of proteins encoded by the mitochondrial genome (mtDNA). We used the yeast Saccharomyces cerevisiae as a model system to systematically identify the genes contributing to this process. Integration of genome-wide high-throughput growth assays with previously published large data sets allowed us to define with high confidence a set of 254 nuclear genes that are indispensable for respiratory growth. Next, we induced loss of mtDNA in the yeast deletion collection by growth on ethidium bromide-containing medium and identified twelve genes that are essential for viability in the absence of mtDNA (i.e. petite-negative). Replenishment of mtDNA by cytoduction showed that respiratory-deficient phenotypes are highly variable in many yeast mutants. Using a mitochondrial genome carrying a selectable marker, ARG8^m, we screened for mutants that are specifically defective in maintenance of mtDNA and mitochondrial protein synthesis. We found that up to 176 nuclear genes are required for expression of mitochondria-encoded proteins during fermentative growth. Taken together, our data provide a comprehensive picture of the molecular processes that are required for respiratory metabolism in a simple eukaryotic cell.

The kinase Isr1 negatively regulates hexosamine biosynthesis in S. cerevisiae

The S. cerevisiae ISR1 gene encodes a putative kinase with no ascribed function. Here, we show that Isr1 acts as a negative regulator of the highly-conserved hexosamine biosynthesis pathway (HBP), which converts glucose into uridine diphosphate N-acetylglucosamine (UDP-GlcNAc), the carbohydrate precursor to protein glycosylation, GPI-anchor formation, and chitin biosynthesis. Overexpression of ISR1 is lethal and, at lower levels, causes sensitivity to tunicamycin and resistance to calcofluor white, implying impaired protein glycosylation and reduced chitin deposition. Gfa1 is the first enzyme in the HBP and is conserved from bacteria and yeast to humans. The lethality caused by ISR1 overexpression is rescued by co-overexpression of GFA1 or exogenous glucosamine, which bypasses GFA1's essential function. Gfa1 is phosphorylated in an Isr1-dependent fashion and mutation of Isr1-dependent sites ameliorates the lethality associated with ISR1 overexpression. Isr1 contains a phosphodegron that is phosphorylated by Pho85 and subsequently ubiquitinated by the SCF-Cdc4 complex, largely confining Isr1 protein levels to the time of bud emergence. Mutation of this phosphodegron stabilizes Isr1 and recapitulates the overexpression phenotypes. As Pho85 is a cell cycle and nutrient responsive kinase, this tight regulation of Isr1 may serve to dynamically regulate flux through the HBP and modulate how the cell's energy resources are converted into structural carbohydrates in response to changing cellular needs.

A Genome-Wide Screen for Genes Affecting Spontaneous Direct-Repeat Recombination in Saccharomyces cerevisiae

Homologous recombination is an important mechanism for genome integrity maintenance, and several homologous recombination genes are mutated in various cancers and cancer-prone syndromes. However, since in some cases homologous recombination can lead to mutagenic outcomes, this pathway must be tightly regulated, and mitotic hyper-recombination is a hallmark of genomic instability. We performed two screens in Saccharomyces cerevisiae for genes that, when deleted, cause hyper-recombination between direct repeats. One was performed with the classical patch and replica-plating method. The other was performed with a high-throughput replica-pinning technique that was designed to detect low-frequency events. This approach allowed us to validate the high-throughput replica-pinning methodology independently of the replicative aging context in which it was developed. Furthermore, by combining the two approaches, we were able to identify and validate 35 genes whose deletion causes elevated spontaneous direct-repeat recombination. Among these are mismatch repair genes, the Sgs1-Top3-Rmi1 complex, the RNase H2 complex, genes involved in the oxidative stress response, and a number of other DNA replication, repair and recombination genes. Since several of our hits are evolutionarily conserved, and repeated elements constitute a significant fraction of mammalian genomes, our work might be relevant for understanding genome integrity maintenance in humans.

Genomic Considerations for the Modification of Saccharomyces cerevisiae for Biofuel and Metabolite Biosynthesis

The growing global population and developing world has put a strain on non-renewable natural resources, such as fuels. The shift to renewable sources will, thus, help meet demands, often through the modification of existing biosynthetic pathways or the introduction of novel pathways into non-native species. There are several useful biosynthetic pathways endogenous to organisms that are not conducive for the scale-up necessary for industrial use. The use of genetic and synthetic biological approaches to engineer these pathways in non-native organisms can help ameliorate these challenges. The budding yeast Saccharomyces cerevisiae offers several advantages for genetic engineering for this purpose due to its widespread use as a model system studied by many researchers. The focus of this review is to present a primer on understanding genomic considerations prior to genetic modification and manipulation of S. cerevisiae. The choice of a site for genetic manipulation can have broad implications on transcription throughout a region and this review will present the current understanding of position effects on transcription.

Single-Gene Deletions Contributing to Loss of Heterozygosity in Saccharomyces cerevisiae: Genome-Wide Screens and Reproducibility

Loss of heterozygosity (LOH) is a phenomenon commonly observed in cancers; the loss of chromosomal regions can be both causal and indicative of underlying genome instability. Yeast has long been used as a model organism to study genetic mechanisms difficult to study in mammalian cells. Studying gene deletions leading to increased LOH in yeast aids our understanding of the processes involved, and guides exploration into the etiology of LOH in cancers. Yet, before in-depth mechanistic studies can occur, candidate genes of interest must be identified. Utilizing the heterozygous Saccharomyces cerevisiae deletion collection (≈ 6500 strains), 217 genes whose disruption leads to increased LOH events at the endogenously heterozygous mating type locus were identified. Our investigation to refine this list of genes to candidates with the most definite impact on LOH includes: secondary testing for LOH impact at an additional locus, gene ontology analysis to determine common gene characteristics, and positional gene enrichment studies to identify chromosomal regions important in LOH events. Further, we conducted extensive comparisons of our data to screens with similar, but distinct methodologies, to further distinguish genes that are more likely to be true contributors to instability due to their reproducibility, and not just identified due to the stochastic nature of LOH. Finally, we selected nine candidate genes and quantitatively measured their impact on LOH as a benchmark for the impact of genes identified in our study. Our data add to the existing body of work and strengthen the evidence of single-gene knockdowns contributing to genome instability.

Genome architecture and stability in the Saccharomyces cerevisiae knockout collection

Despite major progress in defining the functional roles of genes, a complete understanding of their influences is far from being realized, even in relatively simple organisms. A major milestone in this direction arose via the completion of the yeast Saccharomyces cerevisiae gene-knockout collection (YKOC), which has enabled high-throughput reverse genetics, phenotypic screenings and analyses of synthetic-genetic interactions1-3. Ensuing experimental work has also highlighted some inconsistencies and mistakes in the YKOC, or genome instability events that rebalance the effects of specific knockouts4-6, but a complete overview of these is lacking. The identification and analysis of genes that are required for maintaining genomic stability have traditionally relied on reporter assays and on the study of deletions of individual genes, but whole-genome-sequencing technologies now enable-in principle-the direct observation of genome instability globally and at scale. To exploit this opportunity, we sequenced the whole genomes of nearly all of the 4,732 strains comprising the homozygous diploid YKOC. Here, by extracting information on copy-number variation of tandem and interspersed repetitive DNA elements, we describe-for almost every single non-essential gene-the genomic alterations that are induced by its loss. Analysis of this dataset reveals genes that affect the maintenance of various genomic elements, highlights cross-talks between nuclear and mitochondrial genome stability, and shows how strains have genetically adapted to life in the absence of individual non-essential genes.

Connectivity problems on heterogeneous graphs

BACKGROUND

Network connectivity problems are abundant in computational biology research, where graphs are used to represent a range of phenomena: from physical interactions between molecules to more abstract relationships such as gene co-expression. One common challenge in studying biological networks is the need to extract meaningful, small subgraphs out of large databases of potential interactions. A useful abstraction for this task turned out to be the Steiner Network problems: given a reference "database" graph, find a parsimonious subgraph that satisfies a given set of connectivity demands. While this formulation proved useful in a number of instances, the next challenge is to account for the fact that the reference graph may not be static. This can happen for instance, when studying protein measurements in single cells or at different time points, whereby different subsets of conditions can have different protein milieu.

RESULTS AND DISCUSSION

We introduce the condition Steiner Network problem in which we concomitantly consider a set of distinct biological conditions. Each condition is associated with a set of connectivity demands, as well as a set of edges that are assumed to be present in that condition. The goal of this problem is to find a minimal subgraph that satisfies all the demands through paths that are present in the respective condition. We show that introducing multiple conditions as an additional factor makes this problem much harder to approximate. Specifically, we prove that for C conditions, this new problem is NP-hard to approximate to a factor of C - ϵ , for every C ≥ 2 and ϵ > 0 , and that this bound is tight. Moving beyond the worst case, we explore a special set of instances where the reference graph grows monotonically between conditions, and show that this problem admits substantially improved approximation algorithms. We also developed an integer linear programming solver for the general problem and demonstrate its ability to reach optimality with instances from the human protein interaction network.

CONCLUSION

Our results demonstrate that in contrast to most connectivity problems studied in computational biology, accounting for multiplicity of biological conditions adds considerable complexity, which we propose to address with a new solver. Importantly, our results extend to several network connectivity problems that are commonly used in computational biology, such as Prize-Collecting Steiner Tree, and provide insight into the theoretical guarantees for their applications in a multiple condition setting.

Functionally Related Genes Cluster into Genomic Regions That Coordinate Transcription at a Distance in Saccharomyces cerevisiae

The two-dimensional, physical positioning of genes along a chromosome can impact proper transcriptional regulation throughout a genomic region. The transcription of neighboring genes is correlated in a genome-wide manner, which is a characteristic of eukaryotes. Many coregulated gene families can be found clustered with another member of the same set—which can result in adjacent gene coregulation of the pair. Due to the myriad gene families that exhibit a nonrandom genomic distribution, there are likely multiple mechanisms working in concert to properly regulate transcriptional coordination of functionally clustered genes. In this study, we utilized budding yeast in an attempt to elucidate mechanisms that underlie this coregulation: testing and empirically validating the enhancer-promoter hypothesis in this species and reporting that functionally related genes cluster to genomic regions that are more conducive to transcriptional regulation at a distance. These clusters rely, in part, on chromatin maintenance and remodelers to maintain proper transcriptional coordination. Our work provides insight into the mechanisms underlying adjacent gene coregulation.

Balancing gene expression is a fundamental challenge of all cell types. To properly regulate transcription on a genome-wide level, there are myriad mechanisms employed by the cell. One layer to this regulation is through spatial positioning, with particular chromosomal loci exerting an influence on transcription throughout a region. Many coregulated gene families utilize spatial positioning to coordinate transcription, with functionally related genes clustering together which can allow coordinated expression via adjacent gene coregulation. The mechanisms underlying this process have not been elucidated, though there are many coregulated gene families that exhibit this genomic distribution. In the present study, we tested for a role for the enhancer-promoter (EP) hypothesis, which demonstrates that regulatory elements can exert transcriptional effects over a broad distance, in coordinating transcriptional coregulation using budding yeast, Saccharomyces cerevisiae. We empirically validated the EP model, finding that the genomic distance a promoter can affect varies by locus, which can profoundly affect levels of transcription, phenotype, and the extent of transcriptional disruption throughout a genomic region. Using the nitrogen metabolism, ribosomal protein, toxin response, and heat shock gene families as our test case, we report functionally clustered genes localize to genomic loci that are more conducive to transcriptional regulation at a distance compared to the unpaired members of the same families. Furthermore, we report that the coregulation of functional clusters is dependent, in part, on chromatin maintenance and remodeling, providing one mechanism underlying adjacent gene coregulation.

IMPORTANCE The two-dimensional, physical positioning of genes along a chromosome can impact proper transcriptional regulation throughout a genomic region. The transcription of neighboring genes is correlated in a genome-wide manner, which is a characteristic of eukaryotes. Many coregulated gene families can be found clustered with another member of the same set—which can result in adjacent gene coregulation of the pair. Due to the myriad gene families that exhibit a nonrandom genomic distribution, there are likely multiple mechanisms working in concert to properly regulate transcriptional coordination of functionally clustered genes. In this study, we utilized budding yeast in an attempt to elucidate mechanisms that underlie this coregulation: testing and empirically validating the enhancer-promoter hypothesis in this species and reporting that functionally related genes cluster to genomic regions that are more conducive to transcriptional regulation at a distance. These clusters rely, in part, on chromatin maintenance and remodelers to maintain proper transcriptional coordination. Our work provides insight into the mechanisms underlying adjacent gene coregulation.

Cis and trans interactions between genes encoding PAF1 complex and ESCRT machinery components in yeast

Saccharomyces cerevisiae is a commonly used model organism for understanding eukaryotic gene function. However, the close proximity between yeast genes can complicate the interpretation of yeast genetic data, particularly high-throughput data. In this study, we examined the interplay between genes encoding components of the PAF1 complex and VPS36, the gene located next to CDC73 on chromosome XII. The PAF1 complex (Cdc73, Paf1, Ctr9, Leo1, and Rtf1, in yeast) affects RNA levels by affecting transcription, histone modifications, and post-transcriptional RNA processing. The human PAF1 complex is linked to cancer, and in yeast, it has been reported to play a role in telomere biology. Vps36, part of the ESCRT-II complex, is involved in sorting proteins for vacuolar/lysosomal degradation. We document a complex set of genetic interactions, which include an adjacent gene effect between CDC73 and VPS36 and synthetic sickness between vps36Δ and cdc73Δ, paf1Δ, or ctr9Δ. Importantly, paf1Δ and ctr9Δ are synthetically lethal with deletions of other components of the ESCRT-II (SNF8 and VPS25), ESCRT-I (STP22), or ESCRT-III (SNF7) complexes. We found that RNA levels of VPS36, but not other ESCRT components, are positively regulated by all components of the PAF1 complex. Finally, we show that deletion of ESCRT components decreases the telomere length in the S288C yeast genetic background, but not in the W303 background. Together, our results outline complex interactions, in cis and in trans, between genes encoding PAF1 and ESCRT-II complex components that affect telomere function and cell viability in yeast.

AnnotationBustR: an R package to extract subsequences from GenBank annotations

BACKGROUND

DNA sequences are pivotal for a wide array of research in biology. Large sequence databases, like GenBank, provide an amazing resource to utilize DNA sequences for large scale analyses. However, many sequence records on GenBank contain more than one gene or are portions of genomes. Inconsistencies in the way genes are annotated and the numerous synonyms a single gene may be listed under provide major challenges for extracting large numbers of subsequences for comparative analysis across taxa. At present, there is no easy way to extract portions from many GenBank accessions based on annotations where gene names may vary extensively.

RESULTS

The R package AnnotationBustR allows users to extract sequences based on GenBank annotations through the ACNUC retrieval system given search terms of gene synonyms and accession numbers. AnnotationBustR extracts subsequences of interest and then writes them to a FASTA file for users to employ in their research endeavors.

CONCLUSION

FASTA files of extracted subsequences and accession tables generated by AnnotationBustR allow users to quickly find and extract subsequences from GenBank accessions. These sequences can then be incorporated in various analyses, like the construction of phylogenies to test a wide range of ecological and evolutionary hypotheses.

Systematic Identification, Characterization, and Conservation of Adjacent-Gene Coregulation in the Budding Yeast Saccharomyces cerevisiae

The spatial positioning of genes throughout the genome arrangement can alter their expression in many eukaryotic organisms. Often this results in a genomic context-specific effect on transcription. One example of this is through the clustering of functionally related genes, which results in adjacent-gene coregulation in the budding yeast Saccharomyces cerevisiae. In the present study, we set out to systematically characterize the prevalence of this phenomenon, finding the genomic organization of functionally related genes into clusters is a characteristic of myriad gene families. These arrangements are found in many evolutionarily divergent fungi and thus represent a widespread, yet underappreciated, layer of transcriptional regulation.

It is essential that cells orchestrate gene expression for the specific niche that they occupy, and this often requires coordination of the expression of large sets of genes. There are multiple regulatory systems that exist for modulation of gene expression, including the adjacent-gene coregulation of the rRNA and ribosome biogenesis and ribosomal protein families. Both gene families exhibit a nonrandom genomic distribution, often clustered directly adjacent to another member of the same family, which results in a tighter transcriptional coordination among adjacent paired genes than that of the unpaired genes within each regulon and can result in a shared promoter that coordinates expression of the pairs. This nonrandom genomic distribution has been seen in a few functionally related gene families, and many of these functional pairings are conserved across divergent fungal lineages. To date, the significance of these observations has not been extended in a systematic way to characterize how prevalent the role of adjacent-gene coregulation is in transcriptional regulation. In the present study, we systematically analyzed the transcriptional coherence of the functional pairs compared to the singletons within all gene families defined by the Gene Ontology Slim designation, using Saccharomyces cerevisiae as a model system, finding that clusters exhibit a tighter transcriptional correlation under specific contexts. We found that the longer a functional pairing is conserved the tighter its response to broad stress and nutritional responses, that roughly 25% of gene families exhibit a nonrandom genomic distribution, and that many of these clusters are conserved. This suggests that adjacent-gene coregulation is a widespread, yet underappreciated, transcriptional mechanism.

IMPORTANCE The spatial positioning of genes throughout the genome arrangement can alter their expression in many eukaryotic organisms. Often this results in a genomic context-specific effect on transcription. One example of this is through the clustering of functionally related genes, which results in adjacent-gene coregulation in the budding yeast Saccharomyces cerevisiae. In the present study, we set out to systematically characterize the prevalence of this phenomenon, finding the genomic organization of functionally related genes into clusters is a characteristic of myriad gene families. These arrangements are found in many evolutionarily divergent fungi and thus represent a widespread, yet underappreciated, layer of transcriptional regulation.

High-throughput creation and functional profiling of DNA sequence variant libraries using CRISPR-Cas9 in yeast

Construction and characterization of large genetic variant libraries is essential for understanding genome function, but remains challenging. Here, we introduce a Cas9-based approach for generating pools of mutants with defined genetic alterations (deletions, substitutions, and insertions) with an efficiency of 80–100% in yeast, along with methods for tracking their fitness en masse. We demonstrate the utility of our approach by characterizing the DNA helicase SGS1 with small tiling deletion mutants that span the length of the protein and a series of point mutations against highly conserved residues in the protein. In addition, we created a genome-wide library targeting 315 poorly characterized small open reading frames (smORFs, <100 amino acids in length) scattered throughout the yeast genome, and assessed which are vital for growth under various environmental conditions. Our strategy allows fundamental biological questions to be investigated in a high-throughput manner with precision.

Marker-free genetic manipulations in yeast using CRISPR/CAS9 system

The budding yeast is currently one of the major model organisms for the study of a wide variety of biological processes. Genetic manipulation of yeast involves the extensive usage of selectable markers that can lead to undesired effects. Thus, marker-free genetic manipulation in yeast is highly desirable for gene/promoter replacement and various other applications. Here we combine the power of selectable markers followed by CRISPR/CAS9 genome editing for common genetic manipulations in yeast in a marker-free manner. We demonstrate our approach for whole gene and promoter replacements and for high-efficiency operator array integration. Our approach allows the utilization of many thousands of existing strains including library strains for the generation of significant genetic changes in yeast in a marker-free and cloning-free fashion.

Using genetic buffering relationships identified in fission yeast to reveal susceptibilities in cells lacking hamartin or tuberin function

Tuberous sclerosis complex is an autosomal dominant disorder characterized by benign tumors arising from the abnormal activation of mTOR signaling in cells lacking TSC1 (hamartin) or TSC2 (tuberin) activity. To expand the genetic framework surrounding this group of growth regulators, we utilized the model eukaryote Schizosaccharomyces pombe to uncover and characterize genes that buffer the phenotypic effects of mutations in the orthologous tsc1 or tsc2 loci. Our study identified two genes: fft3 (encoding a DNA helicase) and ypa1 (encoding a peptidyle-prolyl cis/trans isomerase). While the deletion of fft3 or ypa1 has little effect in wild-type fission yeast cells, their loss in tsc1Δ or tsc2Δ backgrounds results in severe growth inhibition. These data suggest that the inhibition of Ypa1p or Fft3p might represent an 'Achilles' heel' of cells defective in hamartin/tuberin function. Furthermore, we demonstrate that the interaction between tsc1/tsc2 and ypa1 can be rescued through treatment with the mTOR inhibitor, torin-1, and that ypa1Δ cells are resistant to the glycolytic inhibitor, 2-deoxyglucose. This identifies ypa1 as a novel upstream regulator of mTOR and suggests that the effects of ypa1 loss, together with mTOR activation, combine to result in a cellular maladaptation in energy metabolism that is profoundly inhibitory to growth.

Proteomic identification of histone post-translational modifications and proteins enriched at a DNA double-strand break

Here, we use ChAP-MS (chromatin affinity purification with mass spectrometry), for the affinity purification of a sequence-specific single-copy endogenous chromosomal locus containing a DNA double-strand break (DSB). We found multiple new histone post-translational modifications enriched on chromatin bearing a DSB from budding yeast. One of these, methylation of histone H3 on lysine 125, has not previously been reported. Among over 100 novel proteins enriched at a DSB were the phosphatase Sit4, the RNA pol II degradation factor Def1, the mRNA export protein Yra1 and the HECT E3 ligase Tom1. Each of these proteins was required for resistance to radiomimetics, and many were required for resistance to heat, which we show here to cause a defect in DSB repair in yeast. Yra1 and Def1 were required for DSB repair per se, while Sit4 was required for rapid inactivation of the DNA damage checkpoint after DSB repair. Thus, our unbiased proteomics approach has led to the unexpected discovery of novel roles for these and other proteins in the DNA damage response.

Bug mapping and fitness testing of chemically synthesized chromosome X

Debugging a genome sequence is imperative for successfully building a synthetic genome. As part of the effort to build a designer eukaryotic genome, yeast synthetic chromosome X (synX), designed as 707,459 base pairs, was synthesized chemically. SynX exhibited good fitness under a wide variety of conditions. A highly efficient mapping strategy called pooled PCRTag mapping (PoPM), which can be generalized to any watermarked synthetic chromosome, was developed to identify genetic alterations that affect cell fitness ("bugs"). A series of bugs were corrected that included a large region bearing complex amplifications, a growth defect mapping to a recoded sequence in FIP1, and a loxPsym site affecting promoter function of ATP2 PoPM is a powerful tool for synthetic yeast genome debugging and an efficient strategy for phenotype-genotype mapping.

Network-assisted target identification for haploinsufficiency and homozygous profiling screens

Chemical genomic screens have recently emerged as a systematic approach to drug discovery on a genome-wide scale. Drug target identification and elucidation of the mechanism of action (MoA) of hits from these noisy high-throughput screens remain difficult. Here, we present GIT (Genetic Interaction Network-Assisted Target Identification), a network analysis method for drug target identification in haploinsufficiency profiling (HIP) and homozygous profiling (HOP) screens. With the drug-induced phenotypic fitness defect of the deletion of a gene, GIT also incorporates the fitness defects of the gene's neighbors in the genetic interaction network. On three genome-scale yeast chemical genomic screens, GIT substantially outperforms previous scoring methods on target identification on HIP and HOP assays, respectively. Finally, we showed that by combining HIP and HOP assays, GIT further boosts target identification and reveals potential drug's mechanism of action.

A mutated dph3 gene causes sensitivity of Schizosaccharomyces pombe cells to cytotoxic agents

Dph3 is involved in diphthamide modification of the eukaryotic translation elongation factor eEF2 and in Elongator-mediated modifications of tRNAs, where a 5-methoxycarbonyl-methyl moiety is added to wobble uridines. Lack of such modifications affects protein synthesis due to inaccurate translation of mRNAs at ribosomes. We have discovered that integration of markers at the msh3 locus of Schizosaccharomyces pombe impaired the function of the nearby located dph3 gene. Such integrations rendered cells sensitive to the cytotoxic drugs hydroxyurea and methyl methanesulfonate. We constructed dph3 and msh3 strains with mutated ATG start codons (ATGmut), which allowed investigating drug sensitivity without potential interference by marker insertions. The dph3-ATGmut and a dph3::loxP-ura4-loxM gene disruption strain, but not msh3-ATGmut, turned out to be sensitive to hydroxyurea and methyl methanesulfonate, likewise the strains with cassettes integrated at the msh3 locus. The fungicide sordarin, which inhibits diphthamide modified eEF2 of Saccharomyces cerevisiae, barely affected survival of wild type and msh3Δ S. pombe cells, while the dph3Δ mutant was sensitive. The msh3-ATG mutation, but not dph3Δ or the dph3-ATG mutation caused a defect in mating-type switching, indicating that the ura4 marker at the dph3 locus did not interfere with Msh3 function. We conclude that Dph3 is required for cellular resistance to the fungicide sordarin and to the cytotoxic drugs hydroxyurea and methyl methanesulfonate. This is likely mediated by efficient translation of proteins in response to DNA damage and replication stress.

TheCellMap.org: A Web-Accessible Database for Visualizing and Mining the Global Yeast Genetic Interaction Network

Providing access to quantitative genomic data is key to ensure large-scale data validation and promote new discoveries. TheCellMap.org serves as a central repository for storing and analyzing quantitative genetic interaction data produced by genome-scale Synthetic Genetic Array (SGA) experiments with the budding yeast Saccharomyces cerevisiae. In particular, TheCellMap.org allows users to easily access, visualize, explore, and functionally annotate genetic interactions, or to extract and reorganize subnetworks, using data-driven network layouts in an intuitive and interactive manner.

Characterization of contiguous gene deletions in COL4A6 and COL4A5 in Alport syndrome-diffuse leiomyomatosis

Alport syndrome-diffuse leiomyomatosis (AS-DL, OMIM: 308940) is a rare variant of the X-linked Alport syndrome that shows overgrowth of visceral smooth muscles in the gastrointestinal, respiratory and female reproductive tracts in addition to renal symptoms. AS-DL results from deletions that encompass the 5' ends of the COL4A5 and COL4A6 genes, but deletion breakpoints between COL4A5 and COL4A6 have been determined in only four cases. Here, we characterize deletion breakpoints in five AS-DL patients and show a contiguous COL4A6/COL4A5 deletion in each case. We also demonstrate that eight out of nine deletion alleles involved sequences homologous between COL4A5 and COL4A6. Most breakpoints took place in recognizable transposed elements, including long and short interspersed repeats, DNA transposons and long-terminal repeat retrotransposons. Because deletions involved the bidirectional promoter region in each case, we suggest that the occurrence of leiomyomatosis in AS-DL requires inactivation of both genes. Altogether, our study highlights the importance of homologous recombination involving multiple transposed elements for the development of this continuous gene syndrome and other atypical loss-of-function phenotypes.

Targeted proteome analysis of single-gene deletion strains of Saccharomyces cerevisiae lacking enzymes in the central carbon metabolism

Central carbon metabolism is controlled by modulating the protein abundance profiles of enzymes that maintain the essential systems in living organisms. In this study, metabolic adaptation mechanisms in the model organism Saccharomyces cerevisiae were investigated by direct determination of enzyme abundance levels in 30 wild type and mutant strains. We performed a targeted proteome analysis using S. cerevisiae strains that lack genes encoding the enzymes responsible for central carbon metabolism. Our analysis revealed that at least 30% of the observed variations in enzyme abundance levels could be explained by global regulatory mechanisms. A enzyme-enzyme co-abundance analysis revealed that the abundances of enzyme proteins involved in the trehalose metabolism and glycolysis changed in a coordinated manner under the control of the transcription factors for global regulation. The remaining variations were derived from local mechanisms such as a mutant-specific increase in the abundances of remote enzymes. The proteome data also suggested that, although the functional compensation of the deficient enzyme was attained by using more resources for protein biosynthesis, available resources for the biosynthesis of the enzymes responsible for central metabolism were not abundant in S. cerevisiae cells. These results showed that global and local regulation of enzyme abundance levels shape central carbon metabolism in S. cerevisiae by using a limited resource for protein biosynthesis.

Gene-nutrient interaction markedly influences yeast chronological lifespan

PURPOSE

Research into the genetic mechanisms of aging has expanded rapidly over the past two decades. This has in part been the result of the use of model organisms (particularly yeast, worms and flies) and high-throughput technologies, combined with a growing interest in aging research. Despite this progress, widespread consensus regarding the pathways that are fundamental to the modulation of cellular aging and lifespan for all organisms has been limited due to discrepancies between different studies. We have compared results from published genome-wide, chronological lifespan (CLS) screens of individual gene deletion strains in Saccharomyces cerevisiae in order to identify gene deletion strains with consistent influences on longevity as possible indicators of fundamental aging processes from this single-celled, eukaryotic model organism.

METHODS

Three previous reports have described genetic modifiers of chronological aging in the budding yeast (S. cerevisiae) using the yeast gene deletion strain collection. We performed a comparison among the data sets using correlation and decile distribution analysis to describe concordance between screens and identify strains that consistently increased or decreased CLS. We used gene enrichment analysis in an effort to understand the biology underlying genes identified in multiple studies. We attempted to replicate the different experimental conditions employed by the screens to identify potential sources of variability in CLS worth further investigating.

RESULTS

Among 3209 strains present in all three screens, nine deletions strains were in common in the longest-lived decile (2.80%) and thirteen were in common in the shortest-lived decile (4.05%) of all three screens. Similarly, pairwise overlap between screens was low. When the same comparison was extended to three deciles to include more mutants studied in common between the three screens, enrichment of cellular processes based on gene ontology analysis in the long-lived strains remained very limited. To test the hypothesis that different parental strain auxotrophic requirements or media formulations employed by the respective genome-wide screens might contribute to the lack of concordance, different CLS assay conditions were assessed in combination with strains having different ploidy and auxotrophic requirements (all relevant to differences in the way the three genome-wide CLS screens were performed). This limited but systematic analysis of CLS with respect to auxotrophy, ploidy, and media revealed several instances of gene-nutrient interaction.

CONCLUSIONS

There is surprisingly little overlap between the results of three independently performed genome-wide screens of CLS in S. cerevisiae. However, differences in strain genetic background (ploidy and specific auxotrophic requirements) were present, as well as different media and experimental conditions (e.g., aeration and pooled vs. individual culturing), which, along with stochastic effects such as genetic drift or selection of secondary mutations that suppress the loss of function from gene deletion, could in theory account for some of the lack of consensus between results. Considering the lack of overlap in CLS phenotypes among the set of genes reported by all three screens, and the results of a CLS experiment that systematically tested (incorporating extensive controls) for interactions between variables existing between the screens, we propose that discrepancies can be reconciled through deeper understanding of the influence of cell intrinsic factors such as auxotrophic requirements ploidy status, extrinsic factors such as media composition and aeration, as well as interactions that may occur between them, for example as a result of different pooling vs. individually aging cultures. Such factors may have a more significant impact on CLS outcomes than previously realized. Future studies that systematically account for these contextual factors, and can thus clarify the interactions between genetic and nutrient factors that alter CLS phenotypes, should aid more complete understanding of the underlying biology so that genetic principles of CLS in yeast can be extrapolated to differential cellular aging observed in animal models.

Functional characterisation of long intergenic non-coding RNAs through genetic interaction profiling in Saccharomyces cerevisiae

BACKGROUND

Transcriptome studies have revealed that many eukaryotic genomes are pervasively transcribed producing numerous long non-coding RNAs (lncRNAs). However, only a few lncRNAs have been ascribed a cellular role thus far, with most regulating the expression of adjacent genes. Even less lncRNAs have been annotated as essential hence implying that the majority may be functionally redundant. Therefore, the function of lncRNAs could be illuminated through systematic analysis of their synthetic genetic interactions (GIs).

RESULTS

Here, we employ synthetic genetic array (SGA) in Saccharomyces cerevisiae to identify GIs between long intergenic non-coding RNAs (lincRNAs) and protein-coding genes. We first validate this approach by demonstrating that the telomerase RNA TLC1 displays a GI network that corresponds to its well-described function in telomere length maintenance. We subsequently performed SGA screens on a set of uncharacterised lincRNAs and uncover their connection to diverse cellular processes. One of these lincRNAs, SUT457, exhibits a GI profile associating it to telomere organisation and we consistently demonstrate that SUT457 is required for telomeric overhang homeostasis through an Exo1-dependent pathway. Furthermore, the GI profile of SUT457 is distinct from that of its neighbouring genes suggesting a function independent to its genomic location. Accordingly, we show that ectopic expression of this lincRNA suppresses telomeric overhang accumulation in sut457Δ cells assigning a trans-acting role for SUT457 in telomere biology.

CONCLUSIONS

Overall, our work proposes that systematic application of this genetic approach could determine the functional significance of individual lncRNAs in yeast and other complex organisms.

Direct screening for chromatin status on DNA barcodes in yeast delineates the regulome of H3K79 methylation by Dot1

Given the frequent misregulation of chromatin in cancer, it is important to understand the cellular mechanisms that regulate chromatin structure. However, systematic screening for epigenetic regulators is challenging and often relies on laborious assays or indirect reporter read-outs. Here we describe a strategy, Epi-ID, to directly assess chromatin status in thousands of mutants. In Epi-ID, chromatin status on DNA barcodes is interrogated by chromatin immunoprecipitation followed by deep sequencing, allowing for quantitative comparison of many mutants in parallel. Screening of a barcoded yeast knock-out collection for regulators of histone H3K79 methylation by Dot1 identified all known regulators as well as novel players and processes. These include histone deposition, homologous recombination, and adenosine kinase, which influences the methionine cycle. Gcn5, the acetyltransferase within the SAGA complex, was found to regulate histone methylation and H2B ubiquitination. The concept of Epi-ID is widely applicable and can be readily applied to other chromatin features.

DOI:http://dx.doi.org/10.7554/eLife.18919.001

To fit into the nucleus of eukaryotic cells (which include plant, animal and yeast cells), DNA wraps around histone proteins to form a structure called chromatin. Histones can be modified by a variety of chemical tags, which affect how easily nearby DNA can be accessed by other molecules in the cell. These modifications therefore help to control the activity of the genes encoded in the DNA and other key processes such as DNA repair.

If histone modifications are not regulated correctly, diseases such as cancer may result. Enzymes generally perform the actual modification, but there is another layer of regulation that controls the activity of these enzymes that not much is known about.

The activity of an enzyme that performs a histone modification known as H3K79 methylation (which involves a methyl chemical group being added to a particular region of a particular histone protein) has been linked to some forms of leukemia. Collections of mutant yeast cells can be used to identify the factors that regulate histone modifications in both yeast and human cells. However, current methods that screen for these regulators are time consuming.

To make the search for histone modification regulators more efficient, Vlaming et al. developed a new screening procedure called Epi-ID that can measure the amount of a specific histone modification in thousands of budding yeast mutants at the same time. In Epi-ID, each mutant yeast cell has a unique DNA sequence, or “barcode”. The mutant cells are mixed together and the barcodes that are modified by a particular histone modification – such as H3K79 methylation – are isolated and then counted using a DNA sequencing technique. A high barcode count of a certain mutant indicates that more of the histone modification occurs in that mutant.

Using Epi-ID to survey H3K79 methylation enabled Vlaming et al. to successfully identify all previously known H3K79 methylation regulators, as well several new ones. These new regulators included enzymes that deposit histones on DNA, that carry out DNA repair, and that modify or de-modify histone proteins. To move forward with the newly identified regulators, it will be important to analyze how they control H3K79 methylation in yeast cells and to determine whether the regulators also control H3K79 methylation in human cells.

Finally, Epi-ID can be used to identify regulators of other types of histone modifications. A better understanding of chromatin regulation – and H3K79 methylation regulation in particular – can increase our understanding of diseases in which chromatin is deregulated, and may yield new strategies for the treatment of such diseases.

DOI:http://dx.doi.org/10.7554/eLife.18919.002

High-Throughput Identification of Adaptive Mutations in Experimentally Evolved Yeast Populations

High-throughput sequencing has enabled genetic screens that can rapidly identify mutations that occur during experimental evolution. The presence of a mutation in an evolved lineage does not, however, constitute proof that the mutation is adaptive, given the well-known and widespread phenomenon of genetic hitchhiking, in which a non-adaptive or even detrimental mutation can co-occur in a genome with a beneficial mutation and the combined genotype is carried to high frequency by selection. We approximated the spectrum of possible beneficial mutations in Saccharomyces cerevisiae using sets of single-gene deletions and amplifications of almost all the genes in the S. cerevisiae genome. We determined the fitness effects of each mutation in three different nutrient-limited conditions using pooled competitions followed by barcode sequencing. Although most of the mutations were neutral or deleterious, ~500 of them increased fitness. We then compared those results to the mutations that actually occurred during experimental evolution in the same three nutrient-limited conditions. On average, ~35% of the mutations that occurred during experimental evolution were predicted by the systematic screen to be beneficial. We found that the distribution of fitness effects depended on the selective conditions. In the phosphate-limited and glucose-limited conditions, a large number of beneficial mutations of nearly equivalent, small effects drove the fitness increases. In the sulfate-limited condition, one type of mutation, the amplification of the high-affinity sulfate transporter, dominated. In the absence of that mutation, evolution in the sulfate-limited condition involved mutations in other genes that were not observed previously—but were predicted by the systematic screen. Thus, gross functional screens have the potential to predict and identify adaptive mutations that occur during experimental evolution.

Experimental evolution allows us to observe evolution in real time. New advances in genome sequencing make it trivial to discover the mutations that have arisen in evolved cultures; however, linking those mutations to particular adaptive traits remains difficult. We evaluated the fitness impacts of thousands of single-gene losses and amplifications in yeast. We discovered that only a fraction of the hundreds of possible beneficial mutations were actually detected in evolution experiments performed previously. Our results provide evidence that 35% of the mutations identified in experimentally evolved populations are advantageous and that the distribution of beneficial fitness effects depends on the genetic background and the selective conditions. Furthermore, we show that it is possible to select for alternative mutations that improve fitness by blocking particularly high-fitness routes to adaptation.

Functional Profiling Using the Saccharomyces Genome Deletion Project Collections

The ability to measure and quantify the fitness of an entire organism requires considerably more complex approaches than simply using traditional "omic" methods that examine, for example, the abundance of RNA transcripts, proteins, or metabolites. The yeast deletion collections represent the only systematic, comprehensive set of null alleles for any organism in which such fitness measurements can be assayed. Generated by the Saccharomyces Genome Deletion Project, these collections allow the systematic and parallel analysis of gene functions using any measurable phenotype. The unique 20-bp molecular barcodes engineered into the genome of each deletion strain facilitate the massively parallel analysis of individual fitness. Here, we present functional genomic protocols for use with the yeast deletion collections. We describe how to maintain, propagate, and store the deletion collections and how to perform growth fitness assays on single and parallel screening platforms. Phenotypic fitness analyses of the yeast mutants, described in brief here, provide important insights into biological functions, mechanisms of drug action, and response to environmental stresses. It is important to bear in mind that the specific assays described in this protocol represent some of the many ways in which these collections can be assayed, and in this description particular attention is paid to maximizing throughput using growth as the phenotypic measure.

Genome-wide screen uncovers novel pathways for tRNA processing and nuclear-cytoplasmic dynamics

In this resource, Wu et al. present the first comprehensive unbiased analysis of the role of nearly an entire proteome in tRNA biology and describe 162 novel and 12 previously known Saccharomyces cerevisiae gene products that function in tRNA processing, turnover, and subcellular movement. The findings from this genome-wide screen describe putative novel pathways for tRNA nuclear export and extensive links between tRNA biology and other aspects of cell physiology.

Transfer ribonucleic acids (tRNAs) are essential for protein synthesis. However, key gene products involved in tRNA biogenesis and subcellular movement remain to be discovered. We conducted the first comprehensive unbiased analysis of the role of nearly an entire proteome in tRNA biology and describe 162 novel and 12 previously known Saccharomyces cerevisiae gene products that function in tRNA processing, turnover, and subcellular movement. tRNA nuclear export is of particular interest because it is essential, but the known tRNA exporters (Los1 [exportin-t] and Msn5 [exportin-5]) are unessential. We report that mutations of CRM1 (Exportin-1), MEX67/MTR2 (TAP/p15), and five nucleoporins cause accumulation of unspliced tRNA, a hallmark of defective tRNA nuclear export. CRM1 mutation genetically interacts with los1Δ and causes altered tRNA nuclear–cytoplasmic distribution. The data implicate roles for the protein and mRNA nuclear export machineries in tRNA nuclear export. Mutations of genes encoding actin cytoskeleton components and mitochondrial outer membrane proteins also cause accumulation of unspliced tRNA, likely due to defective splicing on mitochondria. Additional gene products, such as chromatin modification enzymes, have unanticipated effects on pre-tRNA end processing. Thus, this genome-wide screen uncovered putative novel pathways for tRNA nuclear export and extensive links between tRNA biology and other aspects of cell physiology.

Mixed Integer Linear Programming based machine learning approach identifies regulators of telomerase in yeast

Understanding telomere length maintenance mechanisms is central in cancer biology as their dysregulation is one of the hallmarks for immortalization of cancer cells. Important for this well-balanced control is the transcriptional regulation of the telomerase genes. We integrated Mixed Integer Linear Programming models into a comparative machine learning based approach to identify regulatory interactions that best explain the discrepancy of telomerase transcript levels in yeast mutants with deleted regulators showing aberrant telomere length, when compared to mutants with normal telomere length. We uncover novel regulators of telomerase expression, several of which affect histone levels or modifications. In particular, our results point to the transcription factors Sum1, Hst1 and Srb2 as being important for the regulation of EST1 transcription, and we validated the effect of Sum1 experimentally. We compiled our machine learning method leading to a user friendly package for R which can straightforwardly be applied to similar problems integrating gene regulator binding information and expression profiles of samples of e.g. different phenotypes, diseases or treatments.

Chromosome-Specific and Global Effects of Aneuploidy in Saccharomyces cerevisiae

Aneuploidy, an unbalanced karyotype in which one or more chromosomes are present in excess or reduced copy number, causes an array of known phenotypes including proteotoxicity, genomic instability, and slowed proliferation. However, the molecular consequences of aneuploidy are poorly understood and an unbiased investigation into aneuploid cell biology is lacking. We performed high-throughput screens for genes the deletion of which has a synthetic fitness cost in aneuploidy Saccharomyces cerevisiae cells containing single extra chromosomes. This analysis identified genes that, when deleted, decrease the fitness of specific disomic strains as well as those that impair the proliferation of a broad range of aneuploidies. In one case, a chromosome-specific synthetic growth defect could be explained fully by the specific duplication of a single gene on the aneuploid chromosome, highlighting the ability of individual dosage imbalances to cause chromosome-specific phenotypes in aneuploid cells. Deletion of other genes, particularly those involved in protein transport, however, confers synthetic sickness on a broad array of aneuploid strains. Indeed, aneuploid cells, regardless of karyotype, exhibit protein secretion and cell-wall integrity defects. Thus, we were able to use this screen to identify novel cellular consequences of aneuploidy, dependent on both specific chromosome imbalances and caused by many different aneuploid karyotypes. Interestingly, the vast majority of cancer cells are highly aneuploid, so this approach could be of further use in identifying both karyotype-specific and nonspecific stresses exhibited by cancer cells as potential targets for the development of novel cancer therapeutics.

Understanding the regulation of coding and noncoding transcription in cell populations

Whole transcriptome analyses have unveiled the uncomfortable truth that we know less about how transcription is regulated then we thought. In addition to its role in classic promoter-driven transcription of coding RNA, it is now clear that RNA Pol II also drives abundant expression of noncoding RNA. For the majority of this the functional significance remains unclear. Moreover, its regulation and impact are hard to predict because it often proceeds in unexpected ways from cryptic promoters, including by driving convergent antisense transcription from within 3' UTRs. This review suggests that its time to rethink how we envisage gene expression by inclusion of the regulatory architecture of the full genetic locus, and expanding our thinking to encompass the fact that we generally study cells within heterogeneous populations.

Systematic identification and correction of annotation errors in the genetic interaction map of Saccharomyces cerevisiae

The yeast mutant collections are a fundamental tool in deciphering genomic organization and function. Over the last decade, they have been used for the systematic exploration of ∼6 000 000 double gene mutants, identifying and cataloging genetic interactions among them. Here we studied the extent to which these data are prone to neighboring gene effects (NGEs), a phenomenon by which the deletion of a gene affects the expression of adjacent genes along the genome. Analyzing ∼90,000 negative genetic interactions observed to date, we found that more than 10% of them are incorrectly annotated due to NGEs. We developed a novel algorithm, GINGER, to identify and correct erroneous interaction annotations. We validated the algorithm using a comparative analysis of interactions from Schizosaccharomyces pombe. We further showed that our predictions are significantly more concordant with diverse biological data compared to their mis-annotated counterparts. Our work uncovered about 9500 new genetic interactions in yeast.

Annotation inconsistencies beyond sequence similarity-based function prediction - phylogeny and genome structure

The function annotation process in computational biology has increasingly shifted from the traditional characterization of individual biochemical roles of protein molecules to the system-wide detection of entire metabolic pathways and genomic structures. The so-called genome-aware methods broaden misannotation inconsistencies in genome sequences beyond protein function assignments, encompassing phylogenetic anomalies and artifactual genomic regions. We outline three categories of error propagation in databases by providing striking examples – at various levels of appreciation by the community from traditional to emerging, thus raising awareness for future solutions.

A Genetic Screen for Saccharomyces cerevisiae Mutants That Fail to Enter Quiescence

Budding yeast begin the transition to quiescence by prolonging G1 and accumulating limited nutrients. They undergo asymmetric cell divisions, slow cellular expansion, acquire significant stress tolerance and construct elaborate cell walls. These morphologic changes give rise to quiescent (Q) cells, which can be distinguished from three other cell types in a stationary phase culture by flow cytometry. We have used flow cytometry to screen for genes that are required to obtain the quiescent cell fraction. We find that cell wall integrity is critical and these genes may help define quiescence-specific features of the cell wall. Genes required to evade the host innate immune response are common. These may be new targets for antifungal drugs. Acquired thermotolerance is also a common property, and we show that the stress-response transcription factors Msn2 and Msn4 promote quiescence. Many other pathways also contribute, including a subset of genes involved in autophagy, ubiquitin-mediated proteolysis, DNA replication, bud site selection, and cytokinesis.

Screening of high-level 4-hydroxy-2 (or 5)-ethyl-5 (or 2)-methyl-3(2H)-furanone-producing strains from a collection of gene deletion mutants of Saccharomyces cerevisiae

4-Hydroxy-2 (or 5)-ethyl-5 (or 2)-methyl-3(2H)-furanone (HEMF) is an important flavor compound that contributes to the sensory properties of many natural products, particularly soy sauce and soybean paste. The compound exhibits a caramel-like aroma and several important physiological activities, such as strong antioxidant activity. HEMF is produced by yeast species in soy sauce manufacturing; however, the enzymes involved in HEMF production remain unknown, hindering efforts to breed yeasts with high-level HEMF production. In this study, we identified high-level HEMF-producing mutants among a Saccharomyces cerevisiae gene deletion mutant collection. Fourteen deletion mutants were screened as high-level HEMF-producing mutants, and the ADH1 gene deletion mutant (adh1Δ) exhibited the maximum HEMF production capacity. Further investigations of the adh1Δ mutant implied that acetaldehyde accumulation contributes to HEMF production, agreeing with previous findings. Therefore, acetaldehyde might be a precursor for HEMF. The ADH1 gene deletion mutant of Zygosaccharomyces rouxii, which is the dominant strain of yeast found during soy sauce fermentation, also produces HEMF effectively, suggesting that acetaldehyde accumulation might be a benchmark for breeding industrial yeasts with excellent HEMF production abilities.

The yeast deletion collection: a decade of functional genomics

The yeast deletion collections comprise >21,000 mutant strains that carry precise start-to-stop deletions of ∼6000 open reading frames. This collection includes heterozygous and homozygous diploids, and haploids of both MATa and MATα mating types. The yeast deletion collection, or yeast knockout (YKO) set, represents the first and only complete, systematically constructed deletion collection available for any organism. Conceived during the Saccharomyces cerevisiae sequencing project, work on the project began in 1998 and was completed in 2002. The YKO strains have been used in numerous laboratories in >1000 genome-wide screens. This landmark genome project has inspired development of numerous genome-wide technologies in organisms from yeast to man. Notable spinoff technologies include synthetic genetic array and HIPHOP chemogenomics. In this retrospective, we briefly describe the yeast deletion project and some of its most noteworthy biological contributions and the impact that these collections have had on the yeast research community and on genomics in general.

Interplay between nonsense-mediated mRNA decay and DNA damage response pathways reveals that Stn1 and Ten1 are the key CST telomere-cap components

A large and diverse set of proteins, including CST complex, nonsense mediated decay (NMD), and DNA damage response (DDR) proteins, play important roles at the telomere in mammals and yeast. Here, we report that NMD, like the DDR, affects single-stranded DNA (ssDNA) production at uncapped telomeres. Remarkably, we find that the requirement for Cdc13, one of the components of CST, can be efficiently bypassed when aspects of DDR and NMD pathways are inactivated. However, identical genetic interventions do not bypass the need for Stn1 and Ten1, the partners of Cdc13. We show that disabling NMD alters the stoichiometry of CST components at telomeres and permits Stn1 to bind telomeres in the absence of Cdc13. Our data support a model that Stn1 and Ten1 can function in a Cdc13-independent manner and have implications for the function of CST components across eukaryotes.