SGAtools: one-stop analysis and visualization of array-based genetic interaction screens

Ageing-associated long non-coding RNA extends lifespan and reduces translation in non-dividing cells

Genomes produce widespread long non-coding RNAs (lncRNAs) of largely unknown functions. We characterize aal1 (ageing-associated lncRNA), which is induced in quiescent fission yeast cells. Deletion of aal1 shortens the chronological lifespan of non-dividing cells, while ectopic overexpression prolongs their lifespan, indicating that aal1 acts in trans. Overexpression of aal1 represses ribosomal-protein gene expression and inhibits cell growth, and aal1 genetically interacts with coding genes functioning in protein translation. The aal1 lncRNA localizes to the cytoplasm and associates with ribosomes. Notably, aal1 overexpression decreases the cellular ribosome content and inhibits protein translation. The aal1 lncRNA binds to the rpl1901 mRNA, encoding a ribosomal protein. The rpl1901 levels are reduced ~2-fold by aal1, which is sufficient to extend lifespan. Remarkably, the expression of the aal1 lncRNA in Drosophila boosts fly lifespan. We propose that aal1 reduces the ribosome content by decreasing Rpl1901 levels, thus attenuating the translational capacity and promoting longevity. Although aal1 is not conserved, its effect in flies suggests that animals feature related mechanisms that modulate ageing, based on the conserved translational machinery.

Revealing host genome-microbiome networks underlying feed efficiency in dairy cows

Ruminants have the ability to digest human-inedible plant materials, due to the symbiotic relationship with the rumen microbiota. Rumen microbes supply short chain fatty acids, amino acids, and vitamins to dairy cows that are used for maintenance, growth, and lactation functions. The main goal of this study was to investigate gene-microbiome networks underlying feed efficiency traits by integrating genotypic, microbial, and phenotypic data from lactating dairy cows. Data consisted of dry matter intake (DMI), net energy secreted in milk, and residual feed intake (RFI) records, SNP genotype, and 16S rRNA rumen microbial abundances from 448 mid-lactation Holstein cows. We first assessed marginal associations between genotypes and phenotypic and microbial traits through genomic scans, and then, in regions with multiple significant hits, we assessed gene-microbiome-phenotype networks using causal structural learning algorithms. We found significant regions co-localizing the rumen microbiome and feed efficiency traits. Interestingly, we found three types of network relationships: (1) the cow genome directly affects both rumen microbial abundances and feed efficiency traits; (2) the cow genome (Chr3: 116.5 Mb) indirectly affects RFI, mediated by the abundance of Syntrophococcus, Prevotella, and an unknown genus of Class Bacilli; and (3) the cow genome (Chr7: 52.8 Mb and Chr11: 6.1-6.2 Mb) affects the abundance of Rikenellaceae RC9 gut group mediated by DMI. Our findings shed light on how the host genome acts directly and indirectly on the rumen microbiome and feed efficiency traits and the potential benefits of the inclusion of specific microbes in selection indexes or as correlated traits in breeding programs. Overall, the multistep approach described here, combining whole-genome scans and causal network reconstruction, allows us to reveal the relationship between genome and microbiome underlying dairy cow feed efficiency.

pSPObooster: A Plasmid System to Improve Sporulation Efficiency of Saccharomyces cerevisiae Lab Strains

Common Saccharomyces cerevisiae lab yeast strains derived from S288C have meiotic defects and therefore are poor sporulators. Here, we developed a plasmid system containing corrected alleles of the MKT1 and RME1 genes to rescue the meiotic defects and show that standard BY4741 and BY4742 strains containing the plasmid display faster and more efficient sporulation. The plasmid, pSPObooster, can be maintained as an episome and easily cured or stably integrated into the genome at a single locus. We demonstrate the use of pSPObooster in low- and high-throughput yeast genetic manipulations and show that it can expedite both procedures without impacting strain behavior.

Analysis of mitochondrial biogenesis and protein localization by genetic screens and automated imaging

Budding yeast is a laboratory model of a simple eukaryotic cell. Its compact genome is very easy to edit. This allowed to create systematic collections (libraries) of yeast strains where every gene is either perturbed or tagged. Here we review how such collections were used to study mitochondrial biology by doing genetic screens. First, we introduce the principles of yeast genome editing and the basics of its life cycle that are useful for genetic experiments. Then we overview what yeast strain collections were created over the past years. We also describe the creation and the usage of the new generation of SWAP-Tag (SWAT) collections that allow to create custom libraries. We outline the principles of changing the genetic background of whole collections in parallel, and the basics of synthetic genetic array (SGA) approach. Then we review the discoveries that were made using different types of genetic screens focusing on general mitochondrial functions, proteome, and protein targeting pathways. The development of new collections and screening techniques will continue to bring valuable insight into the function of mitochondria and other organelles.

Multi-omic analysis reveals genes and proteins integral to bioactivity of Echinochrome A isolated from the waste stream of the sea urchin industry in Aotearoa New Zealand

Evechinus chloroticus (commonly known as kina) is a sea urchin species endemic to New Zealand. Its roe is a culinary delicacy to the indigenous Māori and a globally exported food product. Echinochrome A (Ech A) is a bioactive compound isolated from the waste product of kina shells and spines; however, the molecular mechanisms of Ech A bioactivity are not well understood, partly due to Ech A never being studied using unbiased genome-wide analysis. To explore the high-value pharmaceutical potential of kina food waste, we obtained unbiased functional genomic and proteomic profiles of yeast cells treated with Echinochrome A. Abundance was measured for 4100 proteins every 30 min for four hours using fluorescent microscopy, resulting in the identification of 92 proteins with significant alterations in protein abundance caused by Ech A treatment that were over-represented with specific changes in DNA replication, repair and RNA binding after 30 min, followed by specific changes in the metabolism of metal ions (specifically iron and copper) from 60-240 min. Further analysis indicated that Ech A chelated iron, and that iron supplementation negated the growth inhibition caused by Ech A. Via a growth-based genome-wide analysis of 4800 gene deletion strains, 20 gene deletion strains were sensitive to Ech A in an iron-dependent manner. These genes were over-represented in the cellular response to oxidative stress, suggesting that Ech A suppressed growth inhibition caused by oxidative stress. Unexpectedly, genes integral to cardiolipin and inositol phosphate biosynthesis were required for Ech A bioactivity. Overall, these results identify genes, proteins, and cellular processes mediating the bioactivity of Ech A. Moreover, we demonstrate unbiased genomic and proteomic methodology that will be useful for characterizing bioactive compounds in food and food waste.

K29-linked free polyubiquitin chains affect ribosome biogenesis and direct ribosomal proteins to the intranuclear quality control compartment

Ribosome assembly requires precise coordination between the production and assembly of ribosomal components. Mutations in ribosomal proteins that inhibit the assembly process or ribosome function are often associated with ribosomopathies, some of which are linked to defects in proteostasis. In this study, we examine the interplay between several yeast proteostasis enzymes, including deubiquitylases (DUBs) Ubp2 and Ubp14, and E3 ligases Ufd4 and Hul5, and we explore their roles in the regulation of the cellular levels of K29-linked unanchored polyubiquitin (polyUb) chains. Accumulating K29-linked unanchored polyUb chains associate with maturing ribosomes to disrupt their assembly, activate the ribosome assembly stress response (RASTR), and lead to the sequestration of ribosomal proteins at the intranuclear quality control compartment (INQ). These findings reveal the physiological relevance of INQ and provide insights into mechanisms of cellular toxicity associated with ribosomopathies.

Genome-wide host-pathway interactions affecting cis-cis-muconic acid production in yeast

The success of forward metabolic engineering depends on a thorough understanding of the behaviour of a heterologous metabolic pathway within its host. We have recently described CRI-SPA, a high-throughput gene editing method enabling the delivery of a metabolic pathway to all strains of the Saccharomyces cerevisiae knock-out library. CRI-SPA systematically quantifies the effect of each modified gene present in the library on product synthesis, providing a complete map of host:pathway interactions. In its first version, CRI-SPA relied on the colour of the product betaxanthins to quantify strains synthesis ability. However, only a few compounds produce a visible or fluorescent phenotype limiting the scope of our approach. Here, we adapt CRI-SPA to onboard a biosensor reporting the interactions between host genes and the synthesis of the colourless product cis-cis-muconic acid (CCM). We phenotype >9,000 genotypes, including both gene knock-out and overexpression, by quantifying the fluorescence of yeast colonies growing in high-density agar arrays. We identify novel metabolic targets belonging to a broad range of cellular functions and confirm their positive impact on CCM biosynthesis. In particular, our data suggests a new interplay between CCM biosynthesis and cytosolic redox through their common interaction with the oxidative pentose phosphate pathway. Our genome-wide exploration of host:pathway interaction opens novel strategies for improved production of CCM in yeast cell factories.

Genetic screens in Saccharomyces cerevisiae identify a role for 40S ribosome recycling factors Tma20 and Tma22 in nonsense-mediated decay

The decay of messenger RNA with a premature termination codon by nonsense-mediated decay (NMD) is an important regulatory pathway for eukaryotes and an essential pathway in mammals. NMD is typically triggered by the ribosome terminating at a stop codon that is aberrantly distant from the poly-A tail. Here, we use a fluorescence screen to identify factors involved in NMD in Saccharomyces cerevisiae. In addition to the known NMD factors, including the entire UPF family (UPF1, UPF2, and UPF3), as well as NMD4 and EBS1, we identify factors known to function in posttermination recycling and characterize their contribution to NMD. These observations in S. cerevisiae expand on data in mammals indicating that the 60S recycling factor ABCE1 is important for NMD by showing that perturbations in factors implicated in 40S recycling also correlate with a loss of NMD.

A Matlab-based application for quantification of yeast cell growth on solid media

Quantitative assessment of growth and survival is a suitable technique in studying biochemical, genetic and physiological processes in the cells. The budding yeast Saccharomyces cerevisiae is one of the most widely used eukaryotic model organisms for studying cellular mechanisms and processes in evolutionarily distant species, including humans. Yeast growth can be evaluated on both liquid and solid media by measuring cell suspension turbidity and colony forming units, respectively. Several software tools utilizing different parameters have been proposed to quantify yeast growth on solid media. Here, we developed a Matlab-based application which provides a rapid and robust quantitative yeast growth analysis from spot plating assay. Spot plating assay is a typical procedure to evaluate yeast growth in low-throughput laboratory settings, including growth on different nutrient sources or treatment with specific stressors. The app has a one-step installation process, a self-explanatory interface and shorter analysis steps compared with previous established methods, providing a useful tool for both expert and non-expert yeast researchers.

The Involvement of YNR069C in Protein Synthesis in the Baker's Yeast, Saccharomyces cerevisiae

Maintaining translation fidelity is a critical step within the process of gene expression. It requires the involvement of numerous regulatory elements to ensure the synthesis of functional proteins. The efficient termination of protein synthesis can play a crucial role in preserving this fidelity. Here, we report on investigating a protein of unknown function, YNR069C (also known as BSC5), for its activity in the process of translation. We observed a significant increase in the bypass of premature stop codons upon the deletion of YNR069C. Interestingly, the genomic arrangement of this ORF suggests a compatible mode of expression reliant on translational readthrough, incorporating the neighboring open reading frame. We also showed that the deletion of YNR069C results in an increase in the rate of translation. Based on our results, we propose that YNR069C may play a role in translation fidelity, impacting the overall quantity and quality of translation. Our genetic interaction analysis supports our hypothesis, associating the role of YNR069C to the regulation of protein synthesis.

Mind your marker: the effect of common auxotrophic markers on complex traits in yeast

Yeast cells are extensively used as a key model organism owing to their highly conserved genome, metabolic pathways, and cell biology processes. To assist in genetic engineering and analysis, laboratory yeast strains typically harbor auxotrophic selection markers. When uncompensated, auxotrophic markers cause significant phenotypic bias compared to prototrophic strains and have different combinatorial influences on the metabolic network. Here, we used BY4741, a laboratory strain commonly used as a "wild type" strain in yeast studies, to generate a set of revertant strains, containing all possible combinations of four common auxotrophic markers (leu2∆, ura3∆, his3∆1, met15∆). We examined the effect of the auxotrophic combinations on complex phenotypes such as resistance to rapamycin, acetic acid, and ethanol. Among the markers, we found that leucine auxotrophy most significantly affected the phenotype. We analyzed the phenotypic bias caused by auxotrophy at the genomic level using a prototrophic version of a genome-wide deletion library and a decreased mRNA perturbation (DAmP) library. Prototrophy was found to suppress rapamycin sensitivity in many mutants previously annotated for the phenotype, raising a possible need for reevaluation of the findings in a native metabolic context. These results reveal a significant phenotypic bias caused by common auxotrophic markers and support the use of prototrophic wild-type strains in yeast research.

α-tubulin regulation by 5' introns in Saccharomyces cerevisiae

Across eukaryotic genomes, multiple α- and β-tubulin genes require regulation to ensure sufficient production of tubulin heterodimers. Features within these gene families that regulate expression remain underexplored. Here, we investigate the role of the 5' intron in regulating α-tubulin expression in Saccharomyces cerevisiae. We find that the intron in the α-tubulin, TUB1, promotes α-tubulin expression and cell fitness during microtubule stress. The role of the TUB1 intron depends on proximity to the TUB1 promoter and sequence features that are distinct from the intron in the alternative α-tubulin isotype, TUB3. These results lead us to perform a screen to identify genes that act with the TUB1 intron. We identified several genes involved in chromatin remodeling, α/β-tubulin heterodimer assembly, and the spindle assembly checkpoint. We propose a model where the TUB1 intron promotes expression from the chromosomal locus and that this may represent a conserved mechanism for tubulin regulation under conditions that require high levels of tubulin production.

Genome-wide identification of resistance genes and cellular analysis of key gene knockout strain under 5-hydroxymethylfurfural stress in Saccharomyces cerevisiae

In bioethanol production, the main by-product, 5-hydroxymethylfurfural (HMF), significantly hinders microbial fermentation. Therefore, it is crucial to explore genes related to HMF tolerance in Saccharomyces cerevisiae for enhancing the tolerance of ethanol fermentation strains. A comprehensive analysis was conducted using genome-wide deletion library scanning and SGAtools, resulting in the identification of 294 genes associated with HMF tolerance in S. cerevisiae. Further KEGG and GO enrichment analysis revealed the involvement of genes OCA1 and SIW14 in the protein phosphorylation pathway, underscoring their role in HMF tolerance. Spot test validation and subcellular structure observation demonstrated that, following a 3-h treatment with 60 mM HMF, the SIW14 gene knockout strain exhibited a 12.68% increase in cells with abnormal endoplasmic reticulum (ER) and a 22.41% increase in the accumulation of reactive oxygen species compared to the BY4741 strain. These findings indicate that the SIW14 gene contributes to the protection of the ER structure within the cell and facilitates the clearance of reactive oxygen species, thereby confirming its significance as a key gene for HMF tolerance in S. cerevisiae.

A Normalization Protocol Reduces Edge Effect in High-Throughput Analyses of Hydroxyurea Hypersensitivity in Fission Yeast

Edge effect denotes better growth of microbial organisms situated at the edge of the solid agar media. Although the precise reason underlying edge effect is unresolved, it is generally attributed to greater nutrient availability with less competing neighbors at the edge. Nonetheless, edge effect constitutes an unavoidable confounding factor that results in misinterpretation of cell fitness, especially in high-throughput screening experiments widely employed for genome-wide investigation using microbial gene knockout or mutant libraries. Here, we visualize edge effect in high-throughput high-density pinning arrays and report a normalization approach based on colony growth rate to quantify drug (hydroxyurea)-hypersensitivity in fission yeast strains. This normalization procedure improved the accuracy of fitness measurement by compensating cell growth rate discrepancy at different locations on the plate and reducing false-positive and -negative frequencies. Our work thus provides a simple and coding-free solution for a struggling problem in robotics-based high-throughput screening experiments.

Misregulation of cell cycle-dependent methylation of budding yeast CENP-A contributes to chromosomal instability

Centromere (CEN) identity is specified epigenetically by specialized nucleosomes containing evolutionarily conserved CEN-specific histone H3 variant CENP-A (Cse4 in Saccharomyces cerevisiae, CENP-A in humans), which is essential for faithful chromosome segregation. However, the epigenetic mechanisms that regulate Cse4 function have not been fully defined. In this study, we show that cell cycle-dependent methylation of Cse4-R37 regulates kinetochore function and high-fidelity chromosome segregation. We generated a custom antibody that specifically recognizes methylated Cse4-R37 and showed that methylation of Cse4 is cell cycle regulated with maximum levels of methylated Cse4-R37 and its enrichment at the CEN chromatin occur in the mitotic cells. Methyl-mimic cse4-R37F mutant exhibits synthetic lethality with kinetochore mutants, reduced levels of CEN-associated kinetochore proteins and chromosome instability (CIN), suggesting that mimicking the methylation of Cse4-R37 throughout the cell cycle is detrimental to faithful chromosome segregation. Our results showed that SPOUT methyltransferase Upa1 contributes to methylation of Cse4-R37 and overexpression of UPA1 leads to CIN phenotype. In summary, our studies have defined a role for cell cycle-regulated methylation of Cse4 in high-fidelity chromosome segregation and highlight an important role of epigenetic modifications such as methylation of kinetochore proteins in preventing CIN, an important hallmark of human cancers.

K29-linked unanchored polyubiquitin chains disrupt ribosome biogenesis and direct ribosomal proteins to the Intranuclear Quality control compartment (INQ)

Ribosome assembly requires precise coordination between the production and assembly of ribosomal components. Mutations in ribosomal proteins that inhibit the assembly process or ribosome function are often associated with Ribosomopathies, some of which are linked to defects in proteostasis. In this study, we examine the interplay between several yeast proteostasis enzymes, including deubiquitylases (DUBs), Ubp2 and Ubp14, and E3 ligases, Ufd4 and Hul5, and we explore their roles in the regulation of the cellular levels of K29-linked unanchored polyubiquitin (polyUb) chains. Accumulating K29-linked unanchored polyUb chains associate with maturing ribosomes to disrupt their assembly, activate the Ribosome assembly stress response (RASTR), and lead to the sequestration of ribosomal proteins at the Intranuclear Quality control compartment (INQ). These findings reveal the physiological relevance of INQ and provide insights into mechanisms of cellular toxicity associated with Ribosomopathies.

Network Analysis Reveals the Molecular Bases of Statin Pleiotropy That Vary with Genetic Background

Many approved drugs are pleiotropic: for example, statins, whose main cholesterol-lowering activity is complemented by anticancer and prodiabetogenic mechanisms involving poorly characterized genetic interaction networks. We investigated these using the Saccharomyces cerevisiae genetic model, where most genetic interactions known are limited to the statin-sensitive S288C genetic background. We therefore broadened our approach by investigating gene interactions to include two statin-resistant genetic backgrounds: UWOPS87-2421 and Y55. Networks were functionally focused by selection of HMG1 and BTS1 mevalonate pathway genes for detection of genetic interactions. Networks, multilayered by genetic background, were analyzed for key genes using network centrality (degree, betweenness, and closeness), pathway enrichment, functional community modules, and Gene Ontology. Specifically, we found modification genes related to dysregulated endocytosis and autophagic cell death. To translate results to human cells, human orthologues were searched for other drug targets, thus identifying candidates for synergistic anticancer bioactivity. IMPORTANCE Atorvastatin is a highly successful drug prescribed to lower cholesterol and prevent cardiovascular disease in millions of people. Though much of its effect comes from inhibiting a key enzyme in the cholesterol biosynthetic pathway, genes in this pathway interact with genes in other pathways, resulting in 15% of patients suffering painful muscular side effects and 50% having inadequate responses. Such multigenic complexity may be unraveled using gene networks assembled from overlapping pairs of genes that complement each other. We used the unique power of yeast genetics to construct genome-wide networks specific to atorvastatin bioactivity in three genetic backgrounds to represent the genetic variation and varying response to atorvastatin in human individuals. We then used algorithms to identify key genes and their associated FDA-approved drugs in the networks, which resulted in the distinction of drugs that may synergistically enhance the known anticancer activity of atorvastatin.

The SAGA histone acetyltransferase module targets SMC5/6 to specific genes

BACKGROUND

Structural Maintenance of Chromosomes (SMC) complexes are molecular machines driving chromatin organization at higher levels. In eukaryotes, three SMC complexes (cohesin, condensin and SMC5/6) play key roles in cohesion, condensation, replication, transcription and DNA repair. Their physical binding to DNA requires accessible chromatin.

RESULTS

We performed a genetic screen in fission yeast to identify novel factors required for SMC5/6 binding to DNA. We identified 79 genes of which histone acetyltransferases (HATs) were the most represented. Genetic and phenotypic analyses suggested a particularly strong functional relationship between the SMC5/6 and SAGA complexes. Furthermore, several SMC5/6 subunits physically interacted with SAGA HAT module components Gcn5 and Ada2. As Gcn5-dependent acetylation facilitates the accessibility of chromatin to DNA-repair proteins, we first analysed the formation of DNA-damage-induced SMC5/6 foci in the Δgcn5 mutant. The SMC5/6 foci formed normally in Δgcn5, suggesting SAGA-independent SMC5/6 localization to DNA-damaged sites. Next, we used Nse4-FLAG chromatin-immunoprecipitation (ChIP-seq) analysis in unchallenged cells to assess SMC5/6 distribution. A significant portion of SMC5/6 accumulated within gene regions in wild-type cells, which was reduced in Δgcn5 and Δada2 mutants. The drop in SMC5/6 levels was also observed in gcn5-E191Q acetyltransferase-dead mutant.

CONCLUSION

Our data show genetic and physical interactions between SMC5/6 and SAGA complexes. The ChIP-seq analysis suggests that SAGA HAT module targets SMC5/6 to specific gene regions and facilitates their accessibility for SMC5/6 loading.

DBP7 and YRF1-6 Are Involved in Cell Sensitivity to LiCl by Regulating the Translation of PGM2 mRNA

Lithium chloride (LiCl) has been widely researched and utilized as a therapeutic option for bipolar disorder (BD). Several pathways, including cell signaling and signal transduction pathways in mammalian cells, are shown to be regulated by LiCl. LiCl can negatively control the expression and activity of PGM2, a phosphoglucomutase that influences sugar metabolism in yeast. In the presence of galactose, when yeast cells are challenged by LiCl, the phosphoglucomutase activity of PGM2p is decreased, causing an increase in the concentration of toxic galactose metabolism intermediates that result in cell sensitivity. Here, we report that the null yeast mutant strains DBP7∆ and YRF1-6∆ exhibit increased LiCl sensitivity on galactose-containing media. Additionally, we demonstrate that DBP7 and YRF1-6 modulate the translational level of PGM2 mRNA, and the observed alteration in translation seems to be associated with the 5'-untranslated region (UTR) of PGM2 mRNA. Furthermore, we observe that DBP7 and YRF1-6 influence, to varying degrees, the translation of other mRNAs that carry different 5'-UTR secondary structures.

SWR1 chromatin remodeling complex prevents mitotic slippage during spindle position checkpoint arrest

Faithful chromosome segregation in budding yeast requires correct positioning of the mitotic spindle along the mother to daughter cell polarity axis. When the anaphase spindle is not correctly positioned, a surveillance mechanism, named as the spindle position checkpoint (SPOC), prevents the progression out of mitosis until correct spindle positioning is achieved. How SPOC works on a molecular level is not well understood. Here we performed a genome-wide genetic screen to search for components required for SPOC. We identified the SWR1 chromatin-remodeling complex (SWR1-C) among several novel factors that are essential for SPOC integrity. Cells lacking SWR1-C were able to activate SPOC upon spindle misorientation but underwent mitotic slippage upon prolonged SPOC arrest. This mitotic slippage required the Cdc14-early anaphase release pathway and other factors including the SAGA (Spt-Ada-Gcn5 acetyltransferase) histone acetyltransferase complex, proteasome components and the mitotic cyclin-dependent kinase inhibitor Sic1. Together, our data establish a novel link between SWR1-C chromatin remodeling and robust checkpoint arrest in late anaphase.

A synthetic genetic array screen for interactions with the RNA helicase DED1 during cell stress in budding yeast

During cellular stress it is essential for cells to alter their gene expression to adapt and survive. Gene expression is regulated at multiple levels, but translation regulation is both a method for rapid changes to the proteome and, as one of the most energy-intensive cellular processes, a way to efficiently redirect cellular resources during stress conditions. Despite this ideal positioning, many of the specifics of how translation is regulated, positively or negatively, during various types of cellular stress remain poorly understood. To further assess this regulation, we examined the essential translation factor Ded1, an RNA helicase that has been previously shown to play important roles in the translational response to cellular stress. In particular, ded1 mutants display an increased resistance to growth inhibition and translation repression induced by the TOR pathway inhibitor, rapamycin, suggesting that normal stress responses are partially defective in these mutants. To gain further insight into Ded1 translational regulation during stress, synthetic genetic array analysis was conducted in the presence of rapamycin with a ded1 mutant and a library of nonessential genes in Saccharomyces cerevisiae to identify positive and negative genetic interactions in an unbiased manner. Here, we report the results of this screen and subsequent network mapping and Gene Ontology-term analysis. Hundreds of candidate interactions were identified, which fell into expected categories, such as ribosomal proteins and amino acid biosynthesis, as well as unexpected ones, including membrane trafficking, sporulation, and protein glycosylation. Therefore, these results provide several specific directions for further comprehensive studies.

Using colony size to measure fitness in Saccharomyces cerevisiae

Competitive fitness assays in liquid culture have been a mainstay for characterizing experimental evolution of microbial populations. Growth of microbial strains has also been extensively characterized by colony size and could serve as a useful alternative if translated to per generation measurements of relative fitness. To examine fitness based on colony size, we established a relationship between cell number and colony size for strains of Saccharomyces cerevisiae robotically pinned onto solid agar plates in a high-density format. This was used to measure growth rates and estimate relative fitness differences between evolved strains and their ancestors. After controlling for edge effects through both normalization and agar-trimming, we found that colony size is a sensitive measure of fitness, capable of detecting 1% differences. While fitnesses determined from liquid and solid mediums were not equivalent, our results demonstrate that colony size provides a sensitive means of measuring fitness that is particularly well suited to measurements across many environments.

Genetic-interaction screens uncover novel biological roles and regulators of transcription factors in fission yeast

In Schizosaccharomyces pombe, systematic analyses of single transcription factor deletion or overexpression strains have made substantial advances in determining the biological roles and target genes of transcription factors, yet these characteristics are still relatively unknown for over a quarter of them. Moreover, the comprehensive list of proteins that regulate transcription factors remains incomplete. To further characterize Schizosaccharomyces pombe transcription factors, we performed synthetic sick/lethality and synthetic dosage lethality screens by synthetic genetic array. Examination of 2,672 transcription factor double deletion strains revealed a sick/lethality interaction frequency of 1.72%. Phenotypic analysis of these sick/lethality strains revealed potential cell cycle roles for several poorly characterized transcription factors, including SPBC56F2.05, SPCC320.03, and SPAC3C7.04. In addition, we examined synthetic dosage lethality interactions between 14 transcription factors and a miniarray of 279 deletion strains, observing a synthetic dosage lethality frequency of 4.99%, which consisted of known and novel transcription factor regulators. The miniarray contained deletions of genes that encode primarily posttranslational-modifying enzymes to identify putative upstream regulators of the transcription factor query strains. We discovered that ubiquitin ligase Ubr1 and its E2/E3-interacting protein, Mub1, degrade the glucose-responsive transcriptional repressor Scr1. Loss of ubr1⁺ or mub1⁺ increased Scr1 protein expression, which resulted in enhanced repression of flocculation through Scr1. The synthetic dosage lethality screen also captured interactions between Scr1 and 2 of its known repressors, Sds23 and Amk2, each affecting flocculation through Scr1 by influencing its nuclear localization. Our study demonstrates that sick/lethality and synthetic dosage lethality screens can be effective in uncovering novel functions and regulators of Schizosaccharomyces pombe transcription factors.

Distinct elongation stalls during translation are linked with distinct pathways for mRNA degradation

Key protein adapters couple translation to mRNA decay on specific classes of problematic mRNAs in eukaryotes. Slow decoding on non-optimal codons leads to codon-optimality-mediated decay (COMD) and prolonged arrest at stall sites leads to no-go decay (NGD). The identities of the decay factors underlying these processes and the mechanisms by which they respond to translational distress remain open areas of investigation. We use carefully designed reporter mRNAs to perform genetic screens and functional assays in Saccharomyces cerevisiae. We characterize the roles of Hel2, Syh1, and Smy2 in coordinating translational repression and mRNA decay on NGD reporter mRNAs, finding that Syh1 and, to a lesser extent its paralog Smy2, act in a distinct pathway from Hel2. This Syh1/Smy2-mediated pathway acts as a redundant, compensatory pathway to elicit NGD when Hel2-dependent NGD is impaired. Importantly, we observe that these NGD factors are not involved in the degradation of mRNAs enriched in non-optimal codons. Further, we establish that a key factor previously implicated in COMD, Not5, contributes modestly to the degradation of an NGD-targeted mRNA. Finally, we use ribosome profiling to reveal distinct ribosomal states associated with each reporter mRNA that readily rationalize the contributions of NGD and COMD factors to degradation of these reporters. Taken together, these results provide new insight into the role of Syh1 and Smy2 in NGD and into the ribosomal states that correlate with the activation of distinct pathways targeting mRNAs for degradation in yeast.

Genome-wide identification of resistance genes and transcriptome regulation in yeast to accommodate ammonium toxicity

Background

Ammonium is an important raw material for biomolecules and life activities, and the toxicity of ammonium is also an important ecological and agricultural issue. Ammonium toxicity in yeast has only recently been discovered, and information on its mechanism is limited. In recent years, environmental pollution caused by nitrogen-containing wastewater has been increasing. In addition, the use of yeast in bioreactors to produce nitrogen-containing compounds has been developed. Therefore, research on resistance mechanisms that allow yeast to grow under conditions of high concentrations of ammonium has become more and more important.

Results

To further understand the resistance mechanism of yeast to grow under high concentration of ammonium, we used NH₄Cl to screen a yeast non-essential gene-deletion library. We identified 61 NH₄Cl-sensitive deletion mutants from approximately 4200 mutants in the library, then 34 of them were confirmed by drop test analysis. Enrichment analysis of these 34 genes showed that biosynthesis metabolism, mitophagy, MAPK signaling, and other pathways may play important roles in NH₄Cl resistance. Transcriptome analysis under NH₄Cl stress revealed 451 significantly upregulated genes and 835 significantly downregulated genes. The genes are mainly enriched in: nitrogen compound metabolic process, cell wall, MAPK signaling pathway, mitophagy, and glycine, serine and threonine metabolism.

Conclusions

Our results present a broad view of biological pathways involved in the response to NH₄Cl stress, and thereby advance our understanding of the resistance genes and cellular transcriptional regulation under high concentration of ammonium.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08742-y.

Genetic background and mistranslation frequency determine the impact of mistranslating tRNASerUGG

Transfer RNA variants increase the frequency of mistranslation, the misincorporation of an amino acid not specified by the “standard” genetic code, to frequencies approaching 10% in yeast and bacteria. Cells cope with these variants by having multiple copies of each tRNA isodecoder and through pathways that deal with proteotoxic stress. In this study, we define the genetic interactions of the gene encoding tRNA^Ser_UGG,G26A, which mistranslates serine at proline codons. Using a collection of yeast temperature-sensitive alleles, we identify negative synthetic genetic interactions between the mistranslating tRNA and 109 alleles representing 91 genes, with nearly half of the genes having roles in RNA processing or protein folding and turnover. By regulating tRNA expression, we then compare the strength of the negative genetic interaction for a subset of identified alleles under differing amounts of mistranslation. The frequency of mistranslation correlated with the impact on cell growth for all strains analyzed; however, there were notable differences in the extent of the synthetic interaction at different frequencies of mistranslation depending on the genetic background. For many of the strains, the extent of the negative interaction with tRNA^Ser_UGG,G26A was proportional to the frequency of mistranslation or only observed at intermediate or high frequencies. For others, the synthetic interaction was approximately equivalent at all frequencies of mistranslation. As humans contain similar mistranslating tRNAs, these results are important when analyzing the impact of tRNA variants on disease, where both the individual’s genetic background and the expression of the mistranslating tRNA variant need to be considered.

Genome-Wide Identification of Cellular Pathways and Key Genes That Respond to Sodium Bicarbonate Stress in Saccharomyces cerevisiae

Sodium bicarbonate (NaHCO₃) is an important inorganic salt. It is not only widely used in industrial production and daily life, but is also the main stress in alkaline saline soil. NaHCO₃ has a strong ability to inhibit the growth of fungi in both natural environment and daily application. However, the mechanism by which fungi respond to NaHCO₃ stress is not fully understood. To further clarify the toxic mechanisms of NaHCO₃ stress and identify the specific cellular genes and pathways involved in NaHCO₃ resistance, we performed genome-wide screening with NaHCO₃ using a Saccharomyces cerevisiae deletion mutant library. A total of 33 deletion mutants with NaHCO₃ sensitivity were identified. Compared with wild-type strains, these mutants had significant growth defects in the medium containing NaHCO₃. Bioinformatics analysis found that the corresponding genes of these mutants are mainly enriched in the cell cycle, mitophagy, cell wall integrity, and signaling pathways. Further study using transcriptomic analysis showed that 309 upregulated and 233 downregulated genes were only responded to NaHCO₃ stress, when compared with yeast transcriptomic data under alkaline and saline stress. Upregulated genes were mainly concentrated in amino acid metabolism, steroid biosynthesis, and cell wall, while downregulated genes were enriched in various cellular metabolisms. In summary, we have identified the cellular pathways and key genes that respond to NaHCO₃ stress in the whole genome, providing resource and direction for understanding NaHCO₃ toxicity and cellular resistance mechanisms.

High-Throughput, High-Precision Colony Phenotyping with Pyphe

Colony fitness screens are powerful approaches for functional genomics and genetics. This protocol describes experimental and computational procedures for assaying the fitness of thousands of microbial strains in numerous conditions in parallel. Data analysis is based on pyphe, an all-in-one bioinformatics toolbox for scanning, image analysis, data normalization, and interpretation. We describe a standard protocol where endpoint colony areas are used as fitness proxy and two variations on this, one using colony growth curves and one using colony viability staining with phloxine B. Different strategies for experimental design, normalization and quality control are discussed. Using these approaches, it is possible to collect hundreds of thousands of data points, with low technical noise levels around 5%, in an experiment typically lasting 2 weeks or less.

Kel1 is a phosphorylation-regulated noise suppressor of the pheromone signaling pathway

Mechanisms have evolved that allow cells to detect signals and generate an appropriate response. The accuracy of these responses relies on the ability of cells to discriminate between signal and noise. How cells filter noise in signaling pathways is not well understood. Here, we analyze noise suppression in the yeast pheromone signaling pathway and show that the poorly characterized protein Kel1 serves as a major noise suppressor and prevents cell death. At the molecular level, Kel1 prevents spontaneous activation of the pheromone response by inhibiting membrane recruitment of Ste5 and Far1. Only a hypophosphorylated form of Kel1 suppresses signaling, reduces noise, and prevents pheromone-associated cell death, and our data indicate that the MAPK Fus3 contributes to Kel1 phosphorylation. Taken together, Kel1 serves as a phospho-regulated suppressor of the pheromone pathway to reduce noise, inhibit spontaneous activation of the pathway, regulate mating efficiency, and prevent pheromone-associated cell death.

Retromer Complex and PI3K Complex II-Related Genes Mediate the Yeast (Saccharomyces cerevisiae) Sodium Metabisulfite Resistance Response

Sodium metabisulfite (Na₂S₂O₅) is widely used as a preservative in the food and wine industry. However, it causes varying degrees of cellular damage to organisms. In order to improve our knowledge regarding its cyto-toxicity, a genome-wide screen using the yeast single deletion collection was performed. Additionally, a total of 162 Na₂S₂O₅-sensitive strains and 16 Na₂S₂O₅-tolerant strains were identified. Among the 162 Na₂S₂O₅ tolerance-related genes, the retromer complex was the top enriched cellular component. Further analysis demonstrated that retromer complex deletion leads to increased sensitivity to Na₂S₂O₅, and that Na₂S₂O₅ can induce mislocalization of retromer complex proteins. Notably, phosphatidylinositol 3-monophosphate kinase (PI3K) complex II, which is important for retromer recruitment to the endosome, might be a potential regulator mediating retromer localization and the yeast Na₂S₂O₅ tolerance response. Na₂S₂O₅ can decrease the protein expressions of Vps34, which is the component of PI3K complex. Therefore, Na₂S₂O₅-mediated retromer redistribution might be caused by the effects of decreased Vps34 expression levels. Moreover, both pharmaceutical inhibition of Vps34 functions and deletions of PI3K complex II-related genes affect cell tolerance to Na₂S₂O₅. The results of our study provide a global picture of cellular components required for Na₂S₂O₅ tolerance and advance our understanding concerning Na₂S₂O₅-induced cytotoxicity effects.

Using a chemical genetic screen to enhance our understanding of the antimicrobial properties of copper

The competitive toxic and stress inducing nature of copper necessitates systems that sequester and export this metal from the cytoplasm of bacterial cells. Several predicted mechanisms of toxicity include the production of reactive oxygen species, thiol depletion, DNA and iron-sulfur cluster disruption. Accompanying these mechanisms include pathways of homeostasis such as chelation, oxidation, and transport. Still, the mechanisms of copper resistance and sensitivity are not fully understood. Furthermore, studies fail to recognize that the response to copper is likely a result of numerous mechanisms, as in the case for homeostasis, in which proteins and enzymes work as a collective to maintain appropriate copper concentrations. In this study we used the Keio collection, an array of 3985 Escherichia coli mutants, each with a deleted non-essential gene, to gain a better understanding of prolonged copper exposure. In short, we recovered two copper homeostatic gene and genes involved in transporting and assembling to be involved in mediating prolonged copper stress under the conditions assessed. The gene coding for the protein TolC was uncovered as a sensitive hit and we demonstrated that tolC, an outer membrane efflux channel, is key in mitigating copper sensitivity. Additionally, the activity of tRNA processing was enriched and the deletion of several proteins involved in import generated copper tolerance. Lastly, key genes belonging to central carbon metabolism and nicotinamide adenine dinucleotide biosynthesis were uncovered as tolerant hits. Overall, this study shows that copper sensitivity and tolerance are a result of numerous mechanisms acting in combination within the cell.

Global analysis of Saccharomyces cerevisiae growth in mucin

Metagenomic profiling of the human gut microbiome has discovered DNA from dietary yeasts like Saccharomyces cerevisiae. However, it is unknown if the S. cerevisiae detected by common metagenomic methods are from dead dietary sources, or from live S. cerevisiae colonizing the gut similar to their close relative Candida albicans. While S. cerevisiae can adapt to minimal oxygen and acidic environments, it has not been explored whether this yeast can metabolize mucin, the large, gel-forming, highly glycosylated proteins representing a major source of carbon in the gut mucosa. We reveal that S. cerevisiae can utilize mucin as their main carbon source, as well as perform both a transcriptome analysis and a chemogenomic screen to identify biological pathways required for this yeast to grow optimally in mucin. In total, 739 genes demonstrate significant differential expression in mucin culture, and deletion of 21 genes impact growth in mucin. Both screens suggest that mitochondrial function is required for proper growth in mucin, and through secondary assays we determine that mucin exposure induces mitogenesis and cellular respiration. We further show that deletion of an uncharacterized ORF, YCR095W-A, led to dysfunction in mitochondrial morphology and oxygen consumption in mucin. Finally, we demonstrate that Yps7, an aspartyl protease and homolog to mucin-degrading proteins in C. albicans, is important for growth on mucin. Collectively, our work serves as the initial step toward establishing how this common dietary fungus can survive in the mucus environment of the human gut.

The amino acid substitution affects cellular response to mistranslation

Mistranslation, the misincorporation of an amino acid not specified by the "standard" genetic code, occurs in all organisms. tRNA variants that increase mistranslation arise spontaneously and engineered tRNAs can achieve mistranslation frequencies approaching 10% in yeast and bacteria. Interestingly, human genomes contain tRNA variants with the potential to mistranslate. Cells cope with increased mistranslation through multiple mechanisms, though high levels cause proteotoxic stress. The goal of this study was to compare the genetic interactions and the impact on transcriptome and cellular growth of two tRNA variants that mistranslate at a similar frequency but create different amino acid substitutions in Saccharomyces cerevisiae. One tRNA variant inserts alanine at proline codons whereas the other inserts serine for arginine. Both tRNAs decreased growth rate, with the effect being greater for arginine to serine than for proline to alanine. The tRNA that substituted serine for arginine resulted in a heat shock response. In contrast, heat shock response was minimal for proline to alanine substitution. Further demonstrating the significance of the amino acid substitution, transcriptome analysis identified unique up- and down-regulated genes in response to each mistranslating tRNA. Number and extent of negative synthetic genetic interactions also differed depending upon type of mistranslation. Based on the unique responses observed for these mistranslating tRNAs, we predict that the potential of mistranslation to exacerbate diseases caused by proteotoxic stress depends on the tRNA variant. Furthermore, based on their unique transcriptomes and genetic interactions, different naturally occurring mistranslating tRNAs have the potential to negatively influence specific diseases.

An Hsp90 co-chaperone links protein folding and degradation and is part of a conserved protein quality control

In this paper, we show that the essential Hsp90 co-chaperone Sgt1 is a member of a general protein quality control network that links folding and degradation through its participation in the degradation of misfolded proteins both in the cytosol and the endoplasmic reticulum (ER). Sgt1-dependent protein degradation acts in a parallel pathway to the ubiquitin ligase (E3) and ubiquitin chain elongase (E4), Hul5, and overproduction of Hul5 partly suppresses defects in cells with reduced Sgt1 activity. Upon proteostatic stress, Sgt1 accumulates transiently, in an Hsp90- and proteasome-dependent manner, with quality control sites (Q-bodies) of both yeast and human cells that co-localize with Vps13, a protein that creates organelle contact sites. Misfolding disease proteins, such as synphilin-1 involved in Parkinson's disease, are also sequestered to these compartments and require Sgt1 for their clearance.

Natural variants suppress mutations in hundreds of essential genes

The consequence of a mutation can be influenced by the context in which it operates. For example, loss of gene function may be tolerated in one genetic background, and lethal in another. The extent to which mutant phenotypes are malleable, the architecture of modifiers and the identities of causal genes remain largely unknown. Here, we measure the fitness effects of ~ 1,100 temperature-sensitive alleles of yeast essential genes in the context of variation from ten different natural genetic backgrounds and map the modifiers for 19 combinations. Altogether, fitness defects for 149 of the 580 tested genes (26%) could be suppressed by genetic variation in at least one yeast strain. Suppression was generally driven by gain-of-function of a single, strong modifier gene, and involved both genes encoding complex or pathway partners suppressing specific temperature-sensitive alleles, as well as general modifiers altering the effect of many alleles. The emerging frequency of suppression and range of possible mechanisms suggest that a substantial fraction of monogenic diseases could be managed by modulating other gene products.

A Genome-Wide Screen in Saccharomyces cerevisiae Reveals a Critical Role for Oxidative Phosphorylation in Cellular Tolerance to Lithium Hexafluorophosphate

Lithium hexafluorophosphate (LiPF₆) is one of the leading electrolytes in lithium-ion batteries, and its usage has increased tremendously in the past few years. Little is known, however, about its potential environmental and biological impacts. In order to improve our understanding of the cytotoxicity of LiPF₆ and the specific cellular response mechanisms to it, we performed a genome-wide screen using a yeast (Saccharomyces cerevisiae) deletion mutant collection and identified 75 gene deletion mutants that showed LiPF₆ sensitivity. Among these, genes associated with mitochondria showed the most enrichment. We also found that LiPF₆ is more toxic to yeast than lithium chloride (LiCl) or sodium hexafluorophosphate (NaPF₆). Physiological analysis showed that a high concentration of LiPF₆ caused mitochondrial damage, reactive oxygen species (ROS) accumulation, and ATP content changes. Compared with the results of previous genome-wide screening for LiCl-sensitive mutants, we found that oxidative phosphorylation-related mutants were specifically hypersensitive to LiPF₆. In these deletion mutants, LiPF₆ treatment resulted in higher ROS production and reduced ATP levels, suggesting that oxidative phosphorylation-related genes were important for counteracting LiPF₆-induced toxicity. Taken together, our results identified genes specifically involved in LiPF₆-modulated toxicity, and demonstrated that oxidative stress and ATP imbalance maybe the driving factors in governing LiPF₆-induced toxicity.

The chaperone-binding activity of the mitochondrial surface receptor Tom70 protects the cytosol against mitoprotein-induced stress

Most mitochondrial proteins are synthesized as precursors in the cytosol and post-translationally transported into mitochondria. The mitochondrial surface protein Tom70 acts at the interface of the cytosol and mitochondria. In vitro import experiments identified Tom70 as targeting receptor, particularly for hydrophobic carriers. Using in vivo methods and high-content screens, we revisit the question of Tom70 function and considerably expand the set of Tom70-dependent mitochondrial proteins. We demonstrate that the crucial activity of Tom70 is its ability to recruit cytosolic chaperones to the outer membrane. Indeed, tethering an unrelated chaperone-binding domain onto the mitochondrial surface complements most of the defects caused by Tom70 deletion. Tom70-mediated chaperone recruitment reduces the proteotoxicity of mitochondrial precursor proteins, particularly of hydrophobic inner membrane proteins. Thus, our work suggests that the predominant function of Tom70 is to tether cytosolic chaperones to the outer mitochondrial membrane, rather than to serve as a mitochondrion-specifying targeting receptor.

A yeast chemogenomic screen identifies pathways that modulate adipic acid toxicity

Adipic acid production by yeast fermentation is gaining attention as a renewable source of platform chemicals for making nylon products. However, adipic acid toxicity inhibits yeast growth and fermentation. Here, we performed a chemogenomic screen in Saccharomyces cerevisiae to understand the cellular basis of adipic acid toxicity. Our screen revealed that KGD1 (a key gene in the tricarboxylic acid cycle) deletion improved tolerance to adipic acid and its toxic precursor, catechol. Conversely, disrupting ergosterol biosynthesis as well as protein trafficking and vacuolar transport resulted in adipic acid hypersensitivity. Notably, we show that adipic acid disrupts the Membrane Compartment of Can1 (MCC) on the plasma membrane and impacts endocytosis. This was evidenced by the rapid internalization of Can1 for vacuolar degradation. As ergosterol is an essential component of the MCC and protein trafficking mechanisms are required for endocytosis, we highlight the importance of these cellular processes in modulating adipic acid toxicity.

•

Deletion of the TCA cycle gene KGD1 improves tolerance to adipic acid and catechol

•

Ergosterol and Pdr12 play non-overlapping roles protecting cell from adipic acid

•

Adipic acid-induced plasma membrane localization of Pdr12 is independent of ergosterol

•

Adipic acid disrupts the Membrane Compartment of Can1 (MCC) and induces endocytosis

Global network analysis in Schizosaccharomyces pombe reveals three distinct consequences of the common 1-kb deletion causing juvenile CLN3 disease

Juvenile CLN3 disease is a recessively inherited paediatric neurodegenerative disorder, with most patients homozygous for a 1-kb intragenic deletion in CLN3. The btn1 gene is the Schizosaccharomyces pombe orthologue of CLN3. Here, we have extended the use of synthetic genetic array (SGA) analyses to delineate functional signatures for two different disease-causing mutations in addition to complete deletion of btn1. We show that genetic-interaction signatures can differ for mutations in the same gene, which helps to dissect their distinct functional effects. The mutation equivalent to the minor transcript arising from the 1-kb deletion (btn1^102–208del) shows a distinct interaction pattern. Taken together, our results imply that the minor 1-kb deletion transcript has three consequences for CLN3: to both lose and retain some inherent functions and to acquire abnormal characteristics. This has particular implications for the therapeutic development of juvenile CLN3 disease. In addition, this proof of concept could be applied to conserved genes for other mendelian disorders or any gene of interest, aiding in the dissection of their functional domains, unpacking the global consequences of disease pathogenesis, and clarifying genotype–phenotype correlations. In doing so, this detail will enhance the goals of personalised medicine to improve treatment outcomes and reduce adverse events.

Methylation of rRNA as a host defense against rampant group II intron retrotransposition

BACKGROUND

Group II introns are mobile retroelements, capable of invading new sites in DNA. They are self-splicing ribozymes that complex with an intron-encoded protein to form a ribonucleoprotein that targets DNA after splicing. These molecules can invade DNA site-specifically, through a process known as retrohoming, or can invade ectopic sites through retrotransposition. Retrotransposition, in particular, can be strongly influenced by both environmental and cellular factors.

RESULTS

To investigate host factors that influence retrotransposition, we performed random insertional mutagenesis using the ISS1 transposon to generate a library of over 1000 mutants in Lactococcus lactis, the native host of the Ll.LtrB group II intron. By screening this library, we identified 92 mutants with increased retrotransposition frequencies (RTP-ups). We found that mutations in amino acid transport and metabolism tended to have increased retrotransposition frequencies. We further explored a subset of these RTP-up mutants, the most striking of which is a mutant in the ribosomal RNA methyltransferase rlmH, which exhibited a reproducible 20-fold increase in retrotransposition frequency. In vitro and in vivo experiments revealed that ribosomes in the rlmH mutant were defective in the m3Ψ modification and exhibited reduced binding to the intron RNA.

CONCLUSIONS

Taken together, our results reinforce the importance of the native host organism in regulating group II intron retrotransposition. In particular, the evidence from the rlmH mutant suggests a role for ribosome modification in limiting rampant retrotransposition.

Population genomics of the pathogenic yeast Candida tropicalis identifies hybrid isolates in environmental samples

Candida tropicalis is a human pathogen that primarily infects the immunocompromised. Whereas the genome of one isolate, C. tropicalis MYA-3404, was originally sequenced in 2009, there have been no large-scale, multi-isolate studies of the genetic and phenotypic diversity of this species. Here, we used whole genome sequencing and phenotyping to characterize 77 isolates of C. tropicalis from clinical and environmental sources from a variety of locations. We show that most C. tropicalis isolates are diploids with approximately 2-6 heterozygous variants per kilobase. The genomes are relatively stable, with few aneuploidies. However, we identified one highly homozygous isolate and six isolates of C. tropicalis with much higher heterozygosity levels ranging from 36-49 heterozygous variants per kilobase. Our analyses show that the heterozygous isolates represent two different hybrid lineages, where the hybrids share one parent (A) with most other C. tropicalis isolates, but the second parent (B or C) differs by at least 4% at the genome level. Four of the sequenced isolates descend from an AB hybridization, and two from an AC hybridization. The hybrids are MTLa/α heterozygotes. Hybridization, or mating, between different parents is therefore common in the evolutionary history of C. tropicalis. The new hybrids were predominantly found in environmental niches, including from soil. Hybridization is therefore unlikely to be associated with virulence. In addition, we used genotype-phenotype correlation and CRISPR-Cas9 editing to identify a genome variant that results in the inability of one isolate to utilize certain branched-chain amino acids as a sole nitrogen source.

LI Detector: a framework for sensitive colony-based screens regardless of the distribution of fitness effects

Microbial growth characteristics have long been used to investigate fundamental questions of biology. Colony-based high-throughput screens enable parallel fitness estimation of thousands of individual strains using colony growth as a proxy for fitness. However, fitness estimation is complicated by spatial biases affecting colony growth, including uneven nutrient distribution, agar surface irregularities, and batch effects. Analytical methods that have been developed to correct for these spatial biases rely on the following assumptions: (1) that fitness effects are normally distributed, and (2) that most genetic perturbations lead to minor changes in fitness. Although reasonable for many applications, these assumptions are not always warranted and can limit the ability to detect small fitness effects. Beneficial fitness effects, in particular, are notoriously difficult to detect under these assumptions. Here, we developed the linear interpolation-based detector (LI Detector) framework to enable sensitive colony-based screening without making prior assumptions about the underlying distribution of fitness effects. The LI Detector uses a grid of reference colonies to assign a relative fitness value to every colony on the plate. We show that the LI Detector is effective in correcting for spatial biases and equally sensitive toward increase and decrease in fitness. LI Detector offers a tunable system that allows the user to identify small fitness effects with unprecedented sensitivity and specificity. LI Detector can be utilized to develop and refine gene-gene and gene-environment interaction networks of colony-forming organisms, including yeast, by increasing the range of fitness effects that can be reliably detected.

τ-SGA: synthetic genetic array analysis for systematically screening and quantifying trigenic interactions in yeast

Systematic complex genetic interaction studies have provided insight into high-order functional redundancies and genetic network wiring of the cell. Here, we describe a method for screening and quantifying trigenic interactions from ordered arrays of yeast strains grown on agar plates as individual colonies. The protocol instructs users on the trigenic synthetic genetic array analysis technique, τ-SGA, for high-throughput screens. The steps describe construction of the double-mutant query strains and the corresponding single-mutant control query strains, which are screened in parallel in two replicates. The screening experimental set-up consists of sequential replica-pinning steps that enable automated mating, meiotic recombination and successive haploid selection steps for the generation of triple mutants, which are scored for colony size as a proxy for fitness, which enables the calculation of trigenic interactions. The procedure described here was used to conduct 422 trigenic interaction screens, which generated ~460,000 yeast triple mutants for trigenic interaction analysis. Users should be familiar with robotic equipment required for high-throughput genetic interaction screens and be proficient at the command line to execute the scoring pipeline. Large-scale screen computational analysis is achieved by using MATLAB pipelines that score raw colony size data to produce τ-SGA interaction scores. Additional recommendations are included for optimizing experimental design and analysis of smaller-scale trigenic interaction screens by using a web-based analysis system, SGAtools. This protocol provides a resource for those who would like to gain a deeper, more practical understanding of trigenic interaction screening and quantification methodology.

Functional and transcriptional profiling of non-coding RNAs in yeast reveal context-dependent phenotypes and in trans effects on the protein regulatory network

Non-coding RNAs (ncRNAs), including the more recently identified Stable Unannotated Transcripts (SUTs) and Cryptic Unstable Transcripts (CUTs), are increasingly being shown to play pivotal roles in the transcriptional and post-transcriptional regulation of genes in eukaryotes. Here, we carried out a large-scale screening of ncRNAs in Saccharomyces cerevisiae, and provide evidence for SUT and CUT function. Phenotypic data on 372 ncRNA deletion strains in 23 different growth conditions were collected, identifying ncRNAs responsible for significant cellular fitness changes. Transcriptome profiles were assembled for 18 haploid ncRNA deletion mutants and 2 essential ncRNA heterozygous deletants. Guided by the resulting RNA-seq data we analysed the genome-wide dysregulation of protein coding genes and non-coding transcripts. Novel functional ncRNAs, SUT125, SUT126, SUT035 and SUT532 that act in trans by modulating transcription factors were identified. Furthermore, we described the impact of SUTs and CUTs in modulating coding gene expression in response to different environmental conditions, regulating important biological process such as respiration (SUT125, SUT126, SUT035, SUT432), steroid biosynthesis (CUT494, SUT053, SUT468) or rRNA processing (SUT075 and snR30). Overall, these data capture and integrate the regulatory and phenotypic network of ncRNAs and protein-coding genes, providing genome-wide evidence of the impact of ncRNAs on cellular homeostasis.

A quarter of the yeast genome comprises non-coding RNA molecules (ncRNAs), which do not translate into proteins but are involved in the regulation of gene expression. ncRNAs can affect nearby genes by physically interfering with their transcription (cis mode of action), or they interact with DNA, proteins or other RNAs to regulate the expression of distant genes (trans mode of action). Examples of cis-acting ncRNAs have been broadly described, however, genome-wide studies to identify functional trans-acting ncRNAs involved in global gene regulation are still lacking. Here, we used a ncRNA yeast deletion collection to score ncRNA impact on cellular function in different environmental conditions. A group of 20 ncRNA deletion mutants with broad fitness diversity were selected to investigate the ncRNA effect on the protein and ncRNA expression network. We showed a high correlation between altered phenotypes and global transcriptional changes, in an environmental dependent manner. We confirmed the trans acting regulation of ncRNAs in the genome and their role in altering the expression of transcription factors. These findings support the notion of the involvement of ncRNAs in fine tuning cellular expression via regulation of transcription factors, as an advantageous RNA-mediated mechanism that can be fast and cost-effective for the cells.

Trigenic Synthetic Genetic Array (τ-SGA) Technique for Complex Interaction Analysis

Complex genetic interactions occur when mutant alleles of multiple genes combine to elicit an unexpected phenotype, which could not be predicted given the expectation based on the combination of phenotypes associated with individual mutant alleles. Trigenic Synthetic Genetic Array (τ-SGA) methodology was developed for the systematic analysis of complex interactions involving combinations of three gene perturbations. With a series of replica pinning steps of the τ-SGA procedure, haploid triple mutants are constructed through automated mating and meiotic recombination. For example, a double-mutant query strain carrying two mutant alleles of interest, such as a deletion allele of a nonessential gene and a conditional temperature-sensitive allele of an essential gene, is crossed to an input array of yeast mutants, such as the diagnostic array set of ~1200 mutants, to generate an output array of triple mutants. The colony-size measurements of the resulting triple mutants are used to estimate cellular fitness and quantify trigenic interactions by incorporating corresponding single- and double-mutant fitness estimates. Trigenic interaction networks can be further analyzed for functional modules using various clustering and enrichment analysis tools. Complex genetic interactions are rich in functional information and provide insight into the genotype-to-phenotype relationship, genome size, and speciation.

Chemical-Genetic Interactions with the Proline Analog L-Azetidine-2-Carboxylic Acid in Saccharomyces cerevisiae

Non-proteinogenic amino acids, such as the proline analog L-azetidine-2-carboxylic acid (AZC), are detrimental to cells because they are mis-incorporated into proteins and lead to proteotoxic stress. Our goal was to identify genes that show chemical-genetic interactions with AZC in Saccharomyces cerevisiae and thus also potentially define the pathways cells use to cope with amino acid mis-incorporation. Screening the yeast deletion and temperature sensitive collections, we found 72 alleles with negative chemical-genetic interactions with AZC treatment and 12 alleles that suppress AZC toxicity. Many of the genes with negative chemical-genetic interactions are involved in protein quality control pathways through the proteasome. Genes involved in actin cytoskeleton organization and endocytosis also had negative chemical-genetic interactions with AZC. Related to this, the number of actin patches per cell increases upon AZC treatment. Many of the same cellular processes were identified to have interactions with proteotoxic stress caused by two other amino acid analogs, canavanine and thialysine, or a mistranslating tRNA variant that mis-incorporates serine at proline codons. Alleles that suppressed AZC-induced toxicity functioned through the amino acid sensing TOR pathway or controlled amino acid permeases required for AZC uptake. Further suggesting the potential of genetic changes to influence the cellular response to proteotoxic stress, overexpressing many of the genes that had a negative chemical-genetic interaction with AZC suppressed AZC toxicity.

High-Throughput Identification of Nuclear Envelope Protein Interactions in Schizosaccharomyces pombe Using an Arrayed Membrane Yeast-Two Hybrid Library

The nuclear envelope (NE) contains a specialized set of integral membrane proteins that maintain nuclear shape and integrity and influence chromatin organization and gene expression. Advances in proteomics techniques and studies in model organisms have identified hundreds of proteins that localize to the NE. However, the function of many of these proteins at the NE remains unclear, in part due to a lack of understanding of the interactions that these proteins participate in at the NE membrane. To assist in the characterization of NE transmembrane protein interactions we developed an arrayed library of integral and peripheral membrane proteins from the fission yeast Schizosaccharomyces pombe for high-throughput screening using the split-ubiquitin based membrane yeast two -hybrid system. We used this approach to characterize protein interactions for three conserved proteins that localize to the inner nuclear membrane: Cut11/Ndc1, Lem2 and Ima1/Samp1/Net5. Additionally, we determined how the interaction network for Cut11 is altered in canonical temperature-sensitive cut11-ts mutants. This library and screening approach is readily applicable to characterizing the interactomes of integral membrane proteins localizing to various subcellular compartments.

Genetic profiling of protein burden and nuclear export overload

Overproduction (op) of proteins triggers cellular defects. One of the consequences of overproduction is the protein burden/cost, which is produced by an overloading of the protein synthesis process. However, the physiology of cells under a protein burden is not well characterized. We performed genetic profiling of protein burden by systematic analysis of genetic interactions between GFP-op, surveying both deletion and temperature-sensitive mutants in budding yeast. We also performed genetic profiling in cells with overproduction of triple-GFP (tGFP), and the nuclear export signal-containing tGFP (NES-tGFP). The mutants specifically interacted with GFP-op were suggestive of unexpected connections between actin-related processes like polarization and the protein burden, which was supported by morphological analysis. The tGFP-op interactions suggested that this protein probe overloads the proteasome, whereas those that interacted with NES-tGFP involved genes encoding components of the nuclear export process, providing a resource for further analysis of the protein burden and nuclear export overload.

Mistranslating tRNA identifies a deleterious S213P mutation in theSaccharomyces cerevisiaeeco1-1allele

Mistranslation occurs when an amino acid not specified by the standard genetic code is incorporated during translation. Since the ribosome does not read the amino acid, tRNA variants aminoacylated with a non-cognate amino acid or containing a non-cognate anticodon dramatically increase the frequency of mistranslation. In a systematic genetic analysis, we identified a suppression interaction between tRNASerUGG, G26A, which mistranslates proline codons by inserting serine, and eco1-1, a temperature sensitive allele of the gene encoding an acetyltransferase required for sister chromatid cohesion. The suppression was partial, with a tRNA that inserts alanine at proline codons and not apparent for a tRNA that inserts serine at arginine codons. Sequencing of the eco1-1 allele revealed a mutation that would convert the highly conserved serine 213 within β7 of the GCN5-related N-acetyltransferase core to proline. Mutation of P213 in eco1-1 back to the wild-type serine restored the function of the enzyme at elevated temperatures. Our results indicate the utility of mistranslating tRNA variants to identify functionally relevant mutations and identify eco1 as a reporter for mistranslation. We propose that mistranslation could be used as a tool to treat genetic disease.