Exploring genetic suppression interactions on a global scale

Genome-wide conditional degron libraries for functional genomics

Functional genomics with libraries of knockout alleles is limited to non-essential genes and convoluted by the potential accumulation of suppressor mutations in knockout backgrounds, which can lead to erroneous functional annotations. To address these limitations, we constructed genome-wide libraries of conditional alleles based on the auxin-inducible degron (AID) system for inducible degradation of AID-tagged proteins in the budding yeast Saccharomyces cerevisiae. First, we determined that N-terminal tagging is at least twice as likely to inadvertently impair protein function across the proteome. We thus constructed two libraries with over 5,600 essential and non-essential proteins fused at the C-terminus with an AID tag and an optional fluorescent protein. Approximately 90% of AID-tagged proteins were degraded in the presence of the auxin analog 5-Ph-IAA, with initial protein abundance and tag accessibility as limiting factors. Genome-wide screens for DNA damage response factors revealed a role for the glucose signaling factor GSF2 in resistance to hydroxyurea, highlighting how the AID libraries extend the yeast genetics toolbox.

SNARE chaperone Sly1 directly mediates close-range vesicle tethering

The essential Golgi protein Sly1 is a member of the Sec1/mammalian Unc-18 (SM) family of SNARE chaperones. Sly1 was originally identified through remarkable gain-of-function alleles that bypass requirements for diverse vesicle tethering factors. Employing genetic analyses and chemically defined reconstitutions of ER-Golgi fusion, we discovered that a loop conserved among Sly1 family members is not only autoinhibitory but also acts as a positive effector. An amphipathic lipid packing sensor (ALPS)-like helix within the loop directly binds high-curvature membranes. Membrane binding is required for relief of Sly1 autoinhibition and also allows Sly1 to directly tether incoming vesicles to the Qa-SNARE on the target organelle. The SLY1-20 mutation bypasses requirements for diverse tethering factors but loses this ability if the tethering activity is impaired. We propose that long-range tethers, including Golgins and multisubunit tethering complexes, hand off vesicles to Sly1, which then tethers at close range to initiate trans-SNARE complex assembly and fusion in the early secretory pathway.

Yeast eIF2A has a minimal role in translation initiation and uORF-mediated translational control in vivo

Initiating translation of most eukaryotic mRNAs depends on recruitment of methionyl initiator tRNA (Met-tRNAi) in a ternary complex (TC) with GTP-bound eukaryotic initiation factor 2 (eIF2) to the small (40S) ribosomal subunit, forming a 43S preinitiation complex (PIC) that attaches to the mRNA and scans the 5'-untranslated region (5' UTR) for an AUG start codon. Previous studies have implicated mammalian eIF2A in GTP-independent binding of Met-tRNAi to the 40S subunit and its recruitment to specialized mRNAs that do not require scanning, and in initiation at non-AUG start codons, when eIF2 function is attenuated by phosphorylation of its α-subunit during stress. The role of eIF2A in translation in vivo is poorly understood however, and it was unknown whether the conserved ortholog in budding yeast can functionally substitute for eIF2. We performed ribosome profiling of a yeast deletion mutant lacking eIF2A and isogenic wild-type (WT) cells in the presence or absence of eIF2α phosphorylation induced by starvation for amino acids isoleucine and valine. Whereas starvation of WT confers changes in translational efficiencies (TEs) of hundreds of mRNAs, the eIF2AΔ mutation conferred no significant TE reductions for any mRNAs in non-starved cells, and it reduced the TEs of only a small number of transcripts in starved cells containing phosphorylated eIF2α. We found no evidence that eliminating eIF2A altered the translation of mRNAs containing putative internal ribosome entry site (IRES) elements, or harboring uORFs initiated by AUG or near-cognate start codons, in non-starved or starved cells. Thus, very few mRNAs (possibly only one) appear to employ eIF2A for Met-tRNAi recruitment in yeast cells, even when eIF2 function is attenuated by stress.

Meta-analysis of dispensable essential genes and their interactions with bypass suppressors

Genes have been historically classified as essential or non-essential based on their requirement for viability. However, genomic mutations can sometimes bypass the requirement for an essential gene, challenging the binary classification of gene essentiality. Such dispensable essential genes represent a valuable model for understanding the incomplete penetrance of loss-of-function mutations often observed in natural populations. Here, we compiled data from multiple studies on essential gene dispensability in Saccharomyces cerevisiae to comprehensively characterize these genes. In analyses spanning different evolutionary timescales, dispensable essential genes exhibited distinct phylogenetic properties compared with other essential and non-essential genes. Integration of interactions with suppressor genes that can bypass the gene essentiality revealed the high functional modularity of the bypass suppression network. Furthermore, dispensable essential and bypass suppressor gene pairs reflected simultaneous changes in the mutational landscape of S. cerevisiae strains. Importantly, species in which dispensable essential genes were non-essential tended to carry bypass suppressor mutations in their genomes. Overall, our study offers a comprehensive view of dispensable essential genes and illustrates how their interactions with bypass suppressors reflect evolutionary outcomes.

Pharmacological interventions for lipid transport disorders

The recent discovery that defects in inter-organelle lipid transport are at the heart of several neurological and neurodegenerative disorders raises the challenge of identifying therapeutic strategies to correct lipid transport defects. This perspective highlights two potential strategies suggested by the study of lipid transport in budding yeast. In the first approach, small molecules are proposed that enhance the lipid transfer activity of VPS13 proteins and thereby compensate for reduced transport. In the second approach, molecules that act as inter-organelle tethers could be used to create artificial contact sites and bypass the loss of endogenous contacts.

Yeast eIF2A has a minimal role in translation initiation and uORF-mediated translational control in vivo

Initiating translation of most eukaryotic mRNAs depends on recruitment of methionyl initiator tRNA (Met-tRNAi) in a ternary complex (TC) with GTP-bound eukaryotic initiation factor 2 (eIF2) to the small (40S) ribosomal subunit, forming a 43S preinitiation complex (PIC) that attaches to the mRNA and scans the 5'-untranslated region (5' UTR) for an AUG start codon. Previous studies have implicated mammalian eIF2A in GTP-independent binding of Met-tRNAi to the 40S subunit and its recruitment to specialized mRNAs that do not require scanning, and in initiation at non-AUG start codons, when eIF2 function is attenuated by phosphorylation of its α-subunit during stress. The role of eIF2A in translation in vivo is poorly understood however, and it was unknown whether the conserved ortholog in budding yeast can functionally substitute for eIF2. We performed ribosome profiling of a yeast deletion mutant lacking eIF2A and isogenic wild-type (WT) cells in the presence or absence of eIF2α phosphorylation induced by starvation for amino acids isoleucine and valine. Whereas starvation of WT confers changes in translational efficiencies (TEs) of hundreds of mRNAs, the eIF2AΔ mutation conferred no significant TE reductions for any mRNAs in non-starved cells, and it reduced the TEs of only a small number of transcripts in starved cells containing phosphorylated eIF2α. We found no evidence that eliminating eIF2A altered the translation of mRNAs containing putative IRES elements, or harboring uORFs initiated by AUG or near-cognate start codons, in non-starved or starved cells. Thus, very few mRNAs (possibly only one) appear to employ eIF2A for Met-tRNAi recruitment in yeast cells, even when eIF2 function is attenuated by stress.

The role of lipid scramblases in regulating lipid distributions at cellular membranes

Glycerophospholipids, sphingolipids and cholesterol assemble into lipid bilayers that form the scaffold of cellular membranes, in which proteins are embedded. Membrane composition and membrane protein profiles differ between plasma and intracellular membranes and between the two leaflets of a membrane. Lipid distributions between two leaflets are mediated by lipid translocases, including flippases and scramblases. Flippases use ATP to catalyze the inward movement of specific lipids between leaflets. In contrast, bidirectional flip-flop movements of lipids across the membrane are mediated by scramblases in an ATP-independent manner. Scramblases have been implicated in disrupting the lipid asymmetry of the plasma membrane, protein glycosylation, autophagosome biogenesis, lipoprotein secretion, lipid droplet formation and communications between organelles. Although scramblases in plasma membranes were identified over 10 years ago, most progress about scramblases localized in intracellular membranes has been made in the last few years. Herein, we review the role of scramblases in regulating lipid distributions in cellular membranes, focusing primarily on intracellular membrane-localized scramblases.

Global analysis of suppressor mutations that rescue human genetic defects

BACKGROUND

Genetic suppression occurs when the deleterious effects of a primary "query" mutation, such as a disease-causing mutation, are rescued by a suppressor mutation elsewhere in the genome.

METHODS

To capture existing knowledge on suppression relationships between human genes, we examined 2,400 published papers for potential interactions identified through either genetic modification of cultured human cells or through association studies in patients.

RESULTS

The resulting network encompassed 476 unique suppression interactions covering a wide spectrum of diseases and biological functions. The interactions frequently linked genes that operate in the same biological process. Suppressors were strongly enriched for genes with a role in stress response or signaling, suggesting that deleterious mutations can often be buffered by modulating signaling cascades or immune responses. Suppressor mutations tended to be deleterious when they occurred in absence of the query mutation, in apparent contrast with their protective role in the presence of the query. We formulated and quantified mechanisms of genetic suppression that could explain 71% of interactions and provided mechanistic insight into disease pathology. Finally, we used these observations to predict suppressor genes in the human genome.

CONCLUSIONS

The global suppression network allowed us to define principles of genetic suppression that were conserved across diseases, model systems, and species. The emerging frequency of suppression interactions among human genes and range of underlying mechanisms, together with the prevalence of suppression in model organisms, suggest that compensatory mutations may exist for most genetic diseases.

An integrated computational strategy to predict personalized cancer drug combinations by reversing drug resistance signatures

Drug resistance currently poses the greatest barrier to cancer treatments. To overcome drug resistance, drug combination therapy has been proposed as a promising treatment strategy. Herein, we present Re-Sensitizing Drug Prediction (RSDP), a novel computational strategy, for predicting the personalized cancer drug combination A + B by reversing the resistance signature of drug A. The process integrates multiple biological features using a robust rank aggregation algorithm, including Connectivity Map, synthetic lethality, synthetic rescue, pathway, and drug target. Bioinformatics assessments revealed that RSDP achieved a relatively accurate prediction performance for identifying personalized combinational re-sensitizing drug B against cell line-specific intrinsic resistance, cell line-specific acquired resistance, and patient-specific intrinsic resistance to drug A. In addition, we developed the largest resource of cell line-specific cancer drug resistance signatures, including intrinsic and acquired resistance, as a byproduct of the proposed strategy. The findings indicate that personalized drug resistance signature reversal is a promising strategy for identifying personalized drug combinations, which may guide future clinical decisions regarding personalized medicine.

Rapid evolutionary repair by secondary perturbation of a primary disrupted transcriptional network

The discrete steps of transcriptional rewiring have been proposed to occur neutrally to ensure steady gene expression under stabilizing selection. A conflict-free switch of a regulon between regulators may require an immediate compensatory evolution to minimize deleterious effects. Here, we perform an evolutionary repair experiment on the Lachancea kluyveri yeast sef1Δ mutant using a suppressor development strategy. Complete loss of SEF1 forces cells to initiate a compensatory process for the pleiotropic defects arising from misexpression of TCA cycle genes. Using different selective conditions, we identify two adaptive loss-of-function mutations of IRA1 and AZF1. Subsequent analyses show that Azf1 is a weak transcriptional activator regulated by the Ras1-PKA pathway. Azf1 loss-of-function triggers extensive gene expression changes responsible for compensatory, beneficial, and trade-off phenotypes. The trade-offs can be alleviated by higher cell density. Our results not only indicate that secondary transcriptional perturbation provides rapid and adaptive mechanisms potentially stabilizing the initial stage of transcriptional rewiring but also suggest how genetic polymorphisms of pleiotropic mutations could be maintained in the population.

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.

A tractable physical model for the yeast polarity predicts epistasis and fitness

Accurate phenotype prediction based on genetic information has numerous societal applications, such as crop design or cellular factories. Epistasis, when biological components interact, complicates modelling phenotypes from genotypes. Here we show an approach to mitigate this complication for polarity establishment in budding yeast, where mechanistic information is abundant. We coarse-grain molecular interactions into a so-called mesotype, which we combine with gene expression noise into a physical cell cycle model. First, we show with computer simulations that the mesotype allows validation of the most current biochemical polarity models by quantitatively matching doubling times. Second, the mesotype elucidates epistasis emergence as exemplified by evaluating the predicted mutational effect of key polarity protein Bem1p when combined with known interactors or under different growth conditions. This example also illustrates how unlikely evolutionary trajectories can become more accessible. The tractability of our biophysically justifiable approach inspires a road-map towards bottom-up modelling complementary to statistical inferences.

This article is part of the theme issue ‘Interdisciplinary approaches to predicting evolutionary biology’.

DNA polymerase δ: A single Pol31 polymorphism suppresses the strain background-specific lethality of Pol32 inactivation in Saccharomyces cerevisiae

The evolutionarily conserved DNA polymerase delta (Polδ) plays several essential roles in eukaryotic DNA replication and repair, responsible for the synthesis of the lagging-strand, lower replicative mutagenesis via its proof-reading exonuclease activity and synthetizes both strands during break-induced replication. In Saccharomyces cerevisiae, the Polδ protein complex consists of three subunits encoded by the POL3, POL31 and POL32 genes. Surprisingly, in contrast to POL3 and POL31, the POL32 gene deletion was found to be viable but lethal in all other eukaryotes, raising the question to which extent the viability of the POL32 deletion in S. cerevisiae was species specific. To address this issue, we inactivated the POL32 gene in 10 evolutionary close or distant S. cerevisiae strains and found that POL32 was either essential (3 strains including SK1), non-essential (5 strains including the reference S288C strain) or confers a slow-growth phenotype (2 strains). Whole-genome sequencing of S288C/SK1 pol32∆ meiotic segregants identified the lethal/suppressor effect of the single Pol31-C43Y polymorphism. Consistently, the introduction of the Pol31-43C allele in the SK1 and West African (WA) pol32∆ mutants was sufficient to restore cell viability and wild-type growth upon introduction of two copies of POL31-43C in the SK1 haploid strain. Reciprocally, introduction of the SK1 POL31-43Y allele in the S288C pol32∆ mutant was lethal. Sequence analyses of the POL31 polymorphisms in the 1,011 yeasts genome dataset correlates with the strict occurrence of the POL31-43Y allele in the yeast African palm wine clade. Differently, the single Pol31-E400G polymorphism confers pol32∆ lethality in the Malaysian strain. In the yeast two-hybrid assay, we observed a weakened interaction between Pol3 and Pol31-43Y versus Pol31-43C suggesting an insufficient level of the Polδ holoenzyme stability/activity. Thus, the enigmatic non-essentiality of Pol32 in S. cerevisiae results from single Pol31 amino acid polymorphism and is clade rather than species specific.

Experimental approaches to study evolutionary cell biology using yeasts

The past century has witnessed tremendous advances in understanding how cells function. Nevertheless, how cellular processes have evolved is still poorly understood. Many studies have highlighted surprising molecular diversity in how cells from diverse species execute the same processes, and advances in comparative genomics are likely to reveal much more molecular diversity than was believed possible until recently. Extant cells remain therefore the product of an evolutionary history that we vastly ignore. Evolutionary cell biology has emerged as a discipline aiming to address this knowledge gap by combining evolutionary, molecular, and cellular biology thinking. Recent studies have shown how even essential molecular processes, such as DNA replication, can undergo fast adaptive evolution under certain laboratory conditions. These developments open new lines of research where the evolution of cellular processes can be investigated experimentally. Yeasts naturally find themselves at the forefront of this research line. Not only do they allow the observation of fast evolutionary adaptation, but they also provide numerous genomic, synthetic, and cellular biology tools already developed by a large community. Here we propose that yeasts can serve as an "evolutionary cell lab" to test hypotheses, principles, and ideas in evolutionary cell biology. We discuss various experimental approaches available for this purpose, and how biology at large can benefit from them.

Epistasis and evolutionary dependencies in human cancers

Cancer evolution is driven by the concerted action of multiple molecular alterations, which emerge and are selected during tumor progression. An alteration is selected when it provides an advantage to the tumor cell. However, the advantage provided by a specific alteration depends on the tumor lineage, cell epigenetic state, and presence of additional alterations. In this case, we say that an evolutionary dependency exists between an alteration and what influences its selection. Epistatic interactions between altered genes lead to evolutionary dependencies (EDs), by favoring or vetoing specific combinations of events. Large-scale cancer genomics studies have discovered examples of such dependencies, and showed that they influence tumor progression, disease phenotypes, and therapeutic response. In the past decade, several algorithmic approaches have been proposed to infer EDs from large-scale genomics datasets. These methods adopt diverse strategies to address common challenges and shed new light on cancer evolutionary trajectories. Here, we review these efforts starting from a simple conceptualization of the problem, presenting the tackled and still unmet needs in the field, and discussing the implications of EDs in cancer biology and precision oncology.

A comparative analysis of telomere length maintenance circuits in fission and budding yeast

The natural ends of the linear eukaryotic chromosomes are protected by telomeres, which also play an important role in aging and cancer development. Telomere length varies between species, but it is strictly controlled in all organisms. The process of Telomere Length Maintenance (TLM) involves many pathways, protein complexes and interactions that were first discovered in budding and fission yeast model organisms (Saccharomyces cerevisiae, Schizosaccharomyces pombe). In particular, large-scale systematic genetic screens in budding yeast uncovered a network of ≈500 genes that, when mutated, cause telomeres to lengthen or to shorten. In contrast, the TLM network in fission yeast remains largely unknown and systematic data is still lacking. In this work we try to close this gap and develop a unified interpretable machine learning framework for TLM gene discovery and phenotype prediction in both species. We demonstrate the utility of our framework in pinpointing the pathways by which TLM homeostasis is maintained and predicting novel TLM genes in fission yeast. The results of this study could be used for better understanding of telomere biology and serve as a step towards the adaptation of computational methods based on telomeric data for human prognosis.

yEvo: experimental evolution in high school classrooms selects for novel mutations that impact clotrimazole resistance in Saccharomyces cerevisiae

Antifungal resistance in pathogenic fungi is a growing global health concern. Nonpathogenic laboratory strains of Saccharomyces cerevisiae are an important model for studying mechanisms of antifungal resistance that are relevant to understanding the same processes in pathogenic fungi. We have developed a series of laboratory modules in which high school students used experimental evolution to study antifungal resistance by isolating azole-resistant S. cerevisiae mutants and examining the genetic basis of resistance. We have sequenced 99 clones from these experiments and found that all possessed mutations previously shown to impact azole resistance, validating our approach. We additionally found recurrent mutations in an mRNA degradation pathway and an uncharacterized mitochondrial protein (Csf1) that have possible mechanistic connections to azole resistance. The scale of replication in this initiative allowed us to identify candidate epistatic interactions, as evidenced by pairs of mutations that occur in the same clone more frequently than expected by chance (positive epistasis) or less frequently (negative epistasis). We validated one of these pairs, a negative epistatic interaction between gain-of-function mutations in the multidrug resistance transcription factors Pdr1 and Pdr3. This high school-university collaboration can serve as a model for involving members of the broader public in the scientific process to make meaningful discoveries in biomedical research.

Viral protein engagement of GBF1 induces host cell vulnerability through synthetic lethality

Viruses co-opt host proteins to carry out their lifecycle. Repurposed host proteins may thus become functionally compromised; a situation analogous to a loss-of-function mutation. We term such host proteins as viral-induced hypomorphs. Cells bearing cancer driver loss-of-function mutations have successfully been targeted with drugs perturbing proteins encoded by the synthetic lethal (SL) partners of cancer-specific mutations. Similarly, SL interactions of viral-induced hypomorphs can potentially be targeted as host-based antiviral therapeutics. Here, we use GBF1, which supports the infection of many RNA viruses, as a proof-of-concept. GBF1 becomes a hypomorph upon interaction with the poliovirus protein 3A. Screening for SL partners of GBF1 revealed ARF1 as the top hit, disruption of which selectively killed cells that synthesize 3A alone or in the context of a poliovirus replicon. Thus, viral protein interactions can induce hypomorphs that render host cells selectively vulnerable to perturbations that leave uninfected cells otherwise unscathed. Exploiting viral-induced vulnerabilities could lead to broad-spectrum antivirals for many viruses, including SARS-CoV-2.

A bacterial pan-genome makes gene essentiality strain-dependent and evolvable

Many bacterial species are represented by a pan-genome, whose genetic repertoire far outstrips that of any single bacterial genome. Here we investigate how a bacterial pan-genome might influence gene essentiality and whether essential genes that are initially critical for the survival of an organism can evolve to become non-essential. By using Transposon insertion sequencing (Tn-seq), whole-genome sequencing and RNA-seq on a set of 36 clinical Streptococcus pneumoniae strains representative of >68% of the species' pan-genome, we identify a species-wide 'essentialome' that can be subdivided into universal, core strain-specific and accessory essential genes. By employing 'forced-evolution experiments', we show that specific genetic changes allow bacteria to bypass essentiality. Moreover, by untangling several genetic mechanisms, we show that gene essentiality can be highly influenced by and/or be dependent on: (1) the composition of the accessory genome, (2) the accumulation of toxic intermediates, (3) functional redundancy, (4) efficient recycling of critical metabolites and (5) pathway rewiring. While this functional characterization underscores the evolvability potential of many essential genes, we also show that genes with differential essentiality remain important antimicrobial drug target candidates, as their inactivation almost always has a severe fitness cost in vivo.

Exploring conditional gene essentiality through systems genetics approaches in yeast

An essential gene encodes for a cellular function that is required for viability. Although viability is a straightforward phenotype to analyze in yeast, defining a gene as essential is not always trivial. Gene essentiality has generally been studied in specific laboratory strains and under standard growth conditions, however, essentiality can vary across species, strains, and environments. Recent systematic studies of gene essentiality revealed that two sets of essential genes exist: core essential genes that are always required for viability and conditional essential genes that vary in essentiality in different genetic and environmental contexts. Here, we review recent advances made in the systematic analysis of gene essentiality in yeast and discuss the properties that distinguish core from context-dependent essential genes.

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.

Functional analysis of the Candida albicans kinome reveals Hrr25 as a regulator of antifungal susceptibility

Candida albicans is a leading cause of death due to systemic fungal infections. Poor patient outcomes are attributable to the limited number of antifungal classes and the increasing prevalence of drug resistance. Protein kinases have emerged as rewarding targets in the development of drugs for diverse diseases, yet kinases remain untapped in the quest for new antifungals. Here, we performed a comprehensive analysis of the C. albicans kinome to identify genes for which loss-of-function confers hypersensitivity to the two most widely deployed antifungals, echinocandins and azoles. Through this analysis, we found a role for the casein kinase 1 (CK1) homologue Hrr25 in regulating tolerance to both antifungals as well as target-mediated echinocandin resistance. Follow-up investigations established that Hrr25 regulates these responses through its interaction with the SBF transcription factor. Thus, we provide insights into the circuitry governing cellular responses to antifungals and implicate Hrr25 as a key mediator of drug resistance.

•

Screening Candida albicans kinase mutants reveals 47 regulators of antifungal tolerance

•

Hrr25 is important for growth and cell wall/membrane stress tolerance

•

Hrr25 enables target-mediated echinocandin resistance

•

Hrr25 interacts with the SBF transcription factor complex

Genome-wide CRISPR screen reveals CLPTM1L as a lipid scramblase required for efficient glycosylphosphatidylinositol biosynthesis

Scramblases translocate lipids across the lipid bilayer without consumption of ATP, thereby regulating lipid distributions in cellular membranes. Cytosol-to-lumen translocation across the endoplasmic reticulum (ER) membrane is a common process among lipid glycoconjugates involved in posttranslational protein modifications in eukaryotes. These translocations are thought to be mediated by specific ER-resident scramblases, but the identity of these proteins and the underlying molecular mechanisms have been elusive. Here, we show that CLPTM1L, an integral membrane protein with eight putative transmembrane domains, is the major lipid scramblase involved in efficient glycosylphosphatidylinositol biosynthesis in the ER membrane. Our results validate the long-standing hypothesis that lipid scramblases ensure the efficient translocations of lipid glycoconjugates across the ER membrane for protein glycosylation pathways.

Glycosylphosphatidylinositols (GPIs) are complex glycolipids that act as membrane anchors of many eukaryotic cell surface proteins. Biosynthesis of GPIs is initiated at the cytosolic face of the endoplasmic reticulum (ER) by generation of N-acetylglucosaminyl-phosphatidylinositol (GlcNAc-PI). The second intermediate, glucosaminyl-phosphatidylinositol (GlcN-PI), is translocated across the membrane to the luminal face for later biosynthetic steps and attachment to proteins. The mechanism of the luminal translocation of GlcN-PI is unclear. Here, we report a genome-wide CRISPR knockout screen of genes required for rescuing GPI-anchored protein expression after addition of chemically synthesized GlcNAc-PI to PIGA-knockout cells that cannot synthesize GlcNAc-PI. We identified CLPTM1L (cleft lip and palate transmembrane protein 1-like), an ER-resident multipass membrane protein, as a GlcN-PI scramblase required for efficient biosynthesis of GPIs. Knockout of CLPTM1L in PIGA-knockout cells impaired the efficient utilization of chemically synthesized GlcNAc-PI and GlcN-PI for GPI biosynthesis. Purified CLPTM1L scrambled GlcN-PI, GlcNAc-PI, PI, and several other phospholipids in vitro. CLPTM1L, a member of the PQ-loop family of proteins, represents a type of lipid scramblase having no structural similarity to known lipid scramblases. Knockout of CLPTM1L in various wild-type mammalian cultured cells partially decreased the level of GPI-anchored proteins. These results suggest that CLPTM1L is the major lipid scramblase involved in cytosol-to-lumen translocation of GlcN-PI across the ER membrane for efficient GPI biosynthesis.

Functional network motifs defined through integration of protein-protein and genetic interactions

Cells are enticingly complex systems. The identification of feedback regulation is critically important for understanding this complexity. Network motifs defined as small graphlets that occur more frequently than expected by chance have revolutionized our understanding of feedback circuits in cellular networks. However, with their definition solely based on statistical over-representation, network motifs often lack biological context, which limits their usefulness. Here, we define functional network motifs (FNMs) through the systematic integration of genetic interaction data that directly inform on functional relationships between genes and encoded proteins. Occurring two orders of magnitude less frequently than conventional network motifs, we found FNMs significantly enriched in genes known to be functionally related. Moreover, our comprehensive analyses of FNMs in yeast showed that they are powerful at capturing both known and putative novel regulatory interactions, thus suggesting a promising strategy towards the systematic identification of feedback regulation in biological networks. Many FNMs appeared as excellent candidates for the prioritization of follow-up biochemical characterization, which is a recurring bottleneck in the targeting of complex diseases. More generally, our work highlights a fruitful avenue for integrating and harnessing genomic network data.

Epistatic genetic interactions govern morphogenesis during sexual reproduction and infection in a global human fungal pathogen

Cellular development is orchestrated by evolutionarily conserved signaling pathways, which are often pleiotropic and involve intra- and interpathway epistatic interactions that form intricate, complex regulatory networks. Cryptococcus species are a group of closely related human fungal pathogens that grow as yeasts yet transition to hyphae during sexual reproduction. Additionally, during infection they can form large, polyploid titan cells that evade immunity and develop drug resistance. Multiple known signaling pathways regulate cellular development, yet how these are coordinated and interact with genetic variation is less well understood. Here, we conducted quantitative trait locus (QTL) analyses of a mapping population generated by sexual reproduction of two parents, only one of which is unisexually fertile. We observed transgressive segregation of the unisexual phenotype among progeny, as well as a large-cell phenotype under mating-inducing conditions. These large-cell progeny were found to produce titan cells both in vitro and in infected animals. Two major QTLs and corresponding quantitative trait genes (QTGs) were identified: RIC8 (encoding a guanine-exchange factor) and CNC06490 (encoding a putative Rho-GTPase activator), both involved in G protein signaling. The two QTGs interact epistatically with each other and with the mating-type locus in phenotypic determination. These findings provide insights into the complex genetics of morphogenesis during unisexual reproduction and pathogenic titan cell formation and illustrate how QTL analysis can be applied to identify epistasis between genes. This study shows that phenotypic outcomes are influenced by the genetic background upon which mutations arise, implicating dynamic, complex genotype-to-phenotype landscapes in fungal pathogens and beyond.

YPK9 and WHI2 Negatively Interact during Oxidative Stress

Yeast PARK9 (YPK9) shares homology with human ATP13A2, which encodes a polyamine transporter implicated in juvenile forms of Parkinson's disease. We used YPK9 to gain insight into how ATP13A2 affects cell growth and sensitivity to oxidative stress. Surprisingly, the YPK9 deletion strain from the Saccharomyces cerevisiae deletion collection (YKO) in wildtype BY4741 (mating type a) grew faster and was more resistant to hydrogen peroxide than a commercial, putative parental BY4741 wildtype strain (BY4741^COM). In contrast, deleting YPK9 from BY4741^COM rendered it very sensitive to hydrogen peroxide, suggesting its background is different from that of the deletion collection. Whole-genome sequencing revealed that BY4741^COM and BY4741^COMypk9∆ contain a novel premature stop codon near the 3' end of WHI2 (WHI2^G1324T), whereas the collection's YPK9 deletion strain contains WHI2, which encodes a 486 amino acid protein, Whi2p. Replacing full-length WHI2 with the sequence coding for the predicted truncation (Whi2p^E442*) rendered strains more sensitive to hydrogen peroxide, whereas the converse replacement rendered them more resistant. The sequences of WHI2 in 20 randomly chosen strains from the collection encode the full-length protein, indicating that the putative parental BY4741 WHI2^G1324T strain's genetic background differs from that of the deletion collection. Examination of WHI2 sequences in several commonly used wildtype S. cerevisiae strains and isolates revealed other Whi2p truncations that might yield altered phenotypes. Together, these results demonstrate a novel premature stop codon in WHI2 that renders yeast sensitive to hydrogen peroxide; they also reveal a negative genetic interaction between WHI2 and YPK9 in the presence of hydrogen peroxide in the BY4741 background.

The metabolic growth limitations of petite cells lacking the mitochondrial genome

Eukaryotic cells can survive the loss of their mitochondrial genome, but consequently suffer from severe growth defects. 'Petite yeasts', characterized by mitochondrial genome loss, are instrumental for studying mitochondrial function and physiology. However, the molecular cause of their reduced growth rate remains an open question. Here we show that petite cells suffer from an insufficient capacity to synthesize glutamate, glutamine, leucine and arginine, negatively impacting their growth. Using a combination of molecular genetics and omics approaches, we demonstrate the evolution of fast growth overcomes these amino acid deficiencies, by alleviating a perturbation in mitochondrial iron metabolism and by restoring a defect in the mitochondrial tricarboxylic acid cycle, caused by aconitase inhibition. Our results hence explain the slow growth of mitochondrial genome-deficient cells with a partial auxotrophy in four amino acids that results from distorted iron metabolism and an inhibited tricarboxylic acid cycle.

Structural basis of the P4B ATPase lipid flippase activity

P4 ATPases are lipid flippases that are phylogenetically grouped into P4A, P4B and P4C clades. The P4A ATPases are heterodimers composed of a catalytic α-subunit and accessory β-subunit, and the structures of several heterodimeric flippases have been reported. The S. cerevisiae Neo1 and its orthologs represent the P4B ATPases, which function as monomeric flippases without a β-subunit. It has been unclear whether monomeric flippases retain the architecture and transport mechanism of the dimeric flippases. Here we report the structure of a P4B ATPase, Neo1, in its E1-ATP, E2P-transition, and E2P states. The structure reveals a conserved architecture as well as highly similar functional intermediate states relative to dimeric flippases. Consistently, structure-guided mutagenesis of residues in the proposed substrate translocation path disrupted Neo1’s ability to establish membrane asymmetry. These observations indicate that evolutionarily distant P4 ATPases use a structurally conserved mechanism for substrate transport.

The P4 ATPase lipid flippases play a crucial role in membrane biogenesis. Here the authors report the structure of the monomeric P4B ATPase Neo1 in several states, clarifying the mechanism of substrate transport.

Decoding 'Unnecessary Complexity': A Law of Complexity and a Concept of Hidden Variation Behind "Missing Heritability" in Precision Medicine

The high hopes for the Human Genome Project and personalized medicine were not met because the relationship between genotypes and phenotypes turned out to be more complex than expected. In a previous study we laid the foundation of a theory of complexity and showed that because of the blind nature of evolution, and molecular and historical contingency, cells have accumulated unnecessary complexity, complexity beyond what is necessary and sufficient to describe an organism. Here we provide empirical evidence and show that unnecessary complexity has become integrated into the genome in the form of redundancy and is relevant to molecular evolution of phenotypic complexity. Unnecessary complexity creates uncertainty between molecular and phenotypic complexity, such that phenotypic complexity (CP) is higher than molecular complexity (CM), which is higher than DNA complexity (CD). The qualitative inequality in complexity is based on the following hierarchy: CP > CM > CD. This law-like relationship holds true for all complex traits, including complex diseases. We present a hypothesis of two types of variation, namely open and closed (hidden) systems, show that hidden variation provides a hitherto undiscovered "third source" of phenotypic variation, beside genotype and environment, and argue that "missing heritability" for some complex diseases is likely to be a case of "diluted heritability". There is a need for radically new ways of thinking about the principles of genotype-phenotype relationship. Understanding how cells use hidden, pathway variation to respond to stress can shed light on why two individuals who share the same risk factors may not develop the same disease, or how cancer cells escape death.

Landscape of functional interactions of human processive ribonucleases revealed by high-throughput siRNA screenings

Processive exoribonucleases are executors of RNA decay. In humans, their physical but not functional interactions were thoughtfully investigated. Here we have screened cells deficient in DIS3, XRN2, EXOSC10, DIS3L, and DIS3L2 with a custom siRNA library and determined their genetic interactions (GIs) with diverse pathways of RNA metabolism. We uncovered a complex network of positive interactions that buffer alterations in RNA degradation and reveal reciprocal cooperation with genes involved in transcription, RNA export, and splicing. Further, we evaluated the functional distinctness of nuclear DIS3 and cytoplasmic DIS3L using a library of all known genes associated with RNA metabolism. Our analysis revealed that DIS3 mutation suppresses RNA splicing deficiency, while DIS3L GIs disclose the interplay of cytoplasmic RNA degradation with nuclear RNA processing. Finally, genome-wide DIS3 GI map uncovered relations with genes not directly involved in RNA metabolism, like microtubule organization or regulation of telomerase activity.

The cost of evolved constitutive lac gene expression is usually, but not always, maintained during evolution of generalist populations

Beneficial mutations can become costly following an environmental change. Compensatory mutations can relieve these costs, while not affecting the selected function, so that the benefits are retained if the environment shifts back to be similar to the one in which the beneficial mutation was originally selected. Compensatory mutations have been extensively studied in the context of antibiotic resistance, responses to specific genetic perturbations, and in the determination of interacting gene network components. Few studies have focused on the role of compensatory mutations during more general adaptation, especially as the result of selection in fluctuating environments where adaptations to different environment components may often involve trade-offs. We examine whether costs of a mutation in lacI, which deregulated the expression of the lac operon in evolving populations of Escherichia coli bacteria, were compensated. This mutation occurred in multiple replicate populations selected in environments that fluctuated between growth on lactose, where the mutation was beneficial, and on glucose, where it was deleterious. We found that compensation for the cost of the lacI mutation was rare, but, when it did occur, it did not negatively affect the selected benefit. Compensation was not more likely to occur in a particular evolution environment. Compensation has the potential to remove pleiotropic costs of adaptation, but its rarity indicates that the circumstances to bring about the phenomenon may be peculiar to each individual or impeded by other selected mutations.

Mutational sources of trans-regulatory variation affecting gene expression in Saccharomyces cerevisiae

Heritable variation in a gene’s expression arises from mutations impacting cis- and trans-acting components of its regulatory network. Here, we investigate how trans-regulatory mutations are distributed within the genome and within a gene regulatory network by identifying and characterizing 69 mutations with trans-regulatory effects on expression of the same focal gene in Saccharomyces cerevisiae. Relative to 1766 mutations without effects on expression of this focal gene, we found that these trans-regulatory mutations were enriched in coding sequences of transcription factors previously predicted to regulate expression of the focal gene. However, over 90% of the trans-regulatory mutations identified mapped to other types of genes involved in diverse biological processes including chromatin state, metabolism, and signal transduction. These data show how genetic changes in diverse types of genes can impact a gene’s expression in trans, revealing properties of trans-regulatory mutations that provide the raw material for trans-regulatory variation segregating within natural populations.

Investigation of RNA metabolism through large-scale genetic interaction profiling in yeast

Gene deletion and gene expression alteration can lead to growth defects that are amplified or reduced when a second mutation is present in the same cells. We performed 154 genetic interaction mapping (GIM) screens with query mutants related with RNA metabolism and estimated the growth rates of about 700 000 double mutant Saccharomyces cerevisiae strains. The tested targets included the gene deletion collection and 900 strains in which essential genes were affected by mRNA destabilization (DAmP). To analyze the results, we developed RECAP, a strategy that validates genetic interaction profiles by comparison with gene co-citation frequency, and identified links between 1471 genes and 117 biological processes. In addition to these large-scale results, we validated both enhancement and suppression of slow growth measured for specific RNA-related pathways. Thus, negative genetic interactions identified a role for the OCA inositol polyphosphate hydrolase complex in mRNA translation initiation. By analysis of suppressors, we found that Puf4, a Pumilio family RNA binding protein, inhibits ribosomal protein Rpl9 function, by acting on a conserved UGUAcauUA motif located downstream the stop codon of the RPL9B mRNA. Altogether, the results and their analysis should represent a useful resource for discovery of gene function in yeast.

Using Budding Yeast to Identify Molecules That Block Cancer Cell 'Mitotic Slippage' Only in the Presence of Mitotic Poisons

Research on the budding yeast Saccharomyces cerevisiae has yielded fundamental discoveries on highly conserved biological pathways and yeast remains the best-studied eukaryotic cell in the world. Studies on the mitotic cell cycle and the discovery of cell cycle checkpoints in budding yeast has led to a detailed, although incomplete, understanding of eukaryotic cell cycle progression. In multicellular eukaryotic organisms, uncontrolled aberrant cell division is the defining feature of cancer. Some of the most successful classes of anti-cancer chemotherapeutic agents are mitotic poisons. Mitotic poisons are thought to function by inducing a mitotic spindle checkpoint-dependent cell cycle arrest, via the assembly of the highly conserved mitotic checkpoint complex (MCC), leading to apoptosis. Even in the presence of mitotic poisons, some cancer cells continue cell division via 'mitotic slippage', which may correlate with a cancer becoming refractory to mitotic poison chemotherapeutic treatments. In this review, knowledge about budding yeast cell cycle control is explored to suggest novel potential drug targets, namely, specific regions in the highly conserved anaphase-promoting complex/cyclosome (APC/C) subunits Apc1 and/or Apc5, and in a specific N-terminal region in the APC/C co-factor cell division cycle 20 (Cdc20), which may yield molecules which block 'mitotic slippage' only in the presence of mitotic poisons.

An atlas of protein-protein interactions across mouse tissues

Cellular processes arise from the dynamic organization of proteins in networks of physical interactions. Mapping the interactome has therefore been a central objective of high-throughput biology. However, the dynamics of protein interactions across physiological contexts remain poorly understood. Here, we develop a quantitative proteomic approach combining protein correlation profiling with stable isotope labeling of mammals (PCP-SILAM) to map the interactomes of seven mouse tissues. The resulting maps provide a proteome-scale survey of interactome rewiring across mammalian tissues, revealing more than 125,000 unique interactions at a quality comparable to the highest-quality human screens. We identify systematic suppression of cross-talk between the evolutionarily ancient housekeeping interactome and younger, tissue-specific modules. Rewired proteins are tightly regulated by multiple cellular mechanisms and are implicated in disease. Our study opens up new avenues to uncover regulatory mechanisms that shape in vivo interactome responses to physiological and pathophysiological stimuli in mammalian systems.

Compensatory epistasis explored by molecular dynamics simulations

A non-negligible proportion of human pathogenic variants are known to be present as wild type in at least some non-human mammalian species. The standard explanation for this finding is that molecular mechanisms of compensatory epistasis can alleviate the mutations' otherwise pathogenic effects. Examples of compensated variants have been described in the literature but the interacting residue(s) postulated to play a compensatory role have rarely been ascertained. In this study, the examination of five human X-chromosomally encoded proteins (FIX, GLA, HPRT1, NDP and OTC) allowed us to identify several candidate compensated variants. Strong evidence for a compensated/compensatory pair of amino acids in the coagulation FIXa protein (involving residues 270 and 271) was found in a variety of mammalian species. Both amino acid residues are located within the 60-loop, spatially close to the 39-loop that performs a key role in coagulation serine proteases. To understand the nature of the underlying interactions, molecular dynamics simulations were performed. The predicted conformational change in the 39-loop consequent to the Glu270Lys substitution (associated with hemophilia B) appears to impair the protein's interaction with its substrate but, importantly, such steric hindrance is largely mitigated in those proteins that carry the compensatory residue (Pro271) at the neighboring amino acid position.

Mapping the multiscale structure of biological systems

Biological systems are by nature multiscale, consisting of subsystems that factor into progressively smaller units in a deeply hierarchical structure. At any level of the hierarchy, an ever-increasing diversity of technologies can be applied to characterize the corresponding biological units and their relations, resulting in large networks of physical or functional proximities-e.g., proximities of amino acids within a protein, of proteins within a complex, or of cell types within a tissue. Here, we review general concepts and progress in using network proximity measures as a basis for creation of multiscale hierarchical maps of biological systems. We discuss the functionalization of these maps to create predictive models, including those useful in translation of genotype to phenotype, along with strategies for model visualization and challenges faced by multiscale modeling in the near future. Collectively, these approaches enable a unified hierarchical approach to biological data, with application from the molecular to the macroscopic.

Genetic Dissection of Vps13 Regulation in Yeast Using Disease Mutations from Human Orthologs

The VPS13 family of proteins have emerged as key players in intracellular lipid transport and human health. Humans have four different VPS13 orthologs, the dysfunction of which leads to different diseases. Yeast has a single VPS13 gene, which encodes a protein that localizes to multiple different membrane contact sites. The yeast vps13Δ mutant is pleiotropic, exhibiting defects in sporulation, protein trafficking, endoplasmic reticulum (ER)-phagy and mitochondrial function. Non-null alleles resulting from missense mutations can be useful reagents for understanding the multiple functions of a gene. The exceptionally large size of Vps13 makes the identification of key residues challenging. As a means to identify critical residues in yeast Vps13, amino acid substitution mutations from VPS13A, B, C and D, associated with human disease, were introduced at the cognate positions of yeast VPS13, some of which created separation-of-function alleles. Phenotypic analyses of these mutants have revealed that the promotion of ER-phagy is a fourth, genetically separable role of VPS13 and provide evidence that co-adaptors at the endosome mediate the activity of VPS13 in vacuolar sorting.

Common polymorphic OTC variants can act as genetic modifiers of enzymatic activity

Understanding the role of common polymorphisms in modulating the clinical phenotype when they co‐occur with a disease‐causing lesion is of critical importance in medical genetics. We explored the impact of apparently neutral common polymorphisms, using the gene encoding the urea cycle enzyme, ornithine transcarbamylase (OTC), as a model system. Distinct combinations of genetic backgrounds embracing two missense polymorphisms were created in cis with the pathogenic p.Arg40His replacement. In vitro enzymatic assays revealed that the polymorphic variants were able to modulate OTC activity both in the presence or absence of the pathogenic lesion. First, we found that the combination of the minor alleles of polymorphisms p.Lys46Arg and p.Gln270Arg significantly enhanced enzymatic activity in the wild‐type protein. Second, enzymatic assays revealed that the minor allele of the p.Gln270Arg polymorphism was capable of ameliorating OTC activity when combined in cis with the pathogenic p.Arg40His replacement. Structural analysis predicted that the minor allele of the p.Gln270Arg polymorphism would serve to stabilize the OTC wild‐type protein, thereby corroborating the results of the experimental assays. Our findings demonstrate the potential importance of cis‐interactions between common polymorphic variants and pathogenic missense mutations and illustrate how standing genetic variation can modulate protein function.

Characterization of micron-scale protein-depleted plasma membrane domains in phosphatidylserine-deficient yeast cells

Membrane phase separation to form micron-scale domains of lipids and proteins occurs in artificial membranes; however, a similar large-scale phase separation has not been reported in the plasma membrane of the living cells. We show here that a stable micron-scale protein-depleted region is generated in the plasma membrane of yeast mutants lacking phosphatidylserine at high temperatures. We named this region the 'void zone'. Transmembrane proteins and certain peripheral membrane proteins and phospholipids are excluded from the void zone. The void zone is rich in ergosterol, and requires ergosterol and sphingolipids for its formation. Such properties are also found in the cholesterol-enriched domains of phase-separated artificial membranes, but the void zone is a novel membrane domain that requires energy and various cellular functions for its formation. The formation of the void zone indicates that the plasma membrane in living cells has the potential to undergo phase separation with certain lipid compositions. We also found that void zones were frequently in contact with vacuoles, in which a membrane domain was also formed at the contact site.

Environmental robustness of the global yeast genetic interaction network

Phenotypes associated with genetic variants can be altered by interactions with other genetic variants (GxG), with the environment (GxE), or both (GxGxE). Yeast genetic interactions have been mapped on a global scale, but the environmental influence on the plasticity of genetic networks has not been examined systematically. To assess environmental rewiring of genetic networks, we examined 14 diverse conditions and scored 30,000 functionally representative yeast gene pairs for dynamic, differential interactions. Different conditions revealed novel differential interactions, which often uncovered functional connections between distantly related gene pairs. However, the majority of observed genetic interactions remained unchanged in different conditions, suggesting that the global yeast genetic interaction network is robust to environmental perturbation and captures the fundamental functional architecture of a eukaryotic cell.

Analysing the yeast complexome-the Complex Portal rising to the challenge

The EMBL-EBI Complex Portal is a knowledgebase of macromolecular complexes providing persistent stable identifiers. Entries are linked to literature evidence and provide details of complex membership, function, structure and complex-specific Gene Ontology annotations. Data are freely available and downloadable in HUPO-PSI community standards and missing entries can be requested for curation. In collaboration with Saccharomyces Genome Database and UniProt, the yeast complexome, a compendium of all known heteromeric assemblies from the model organism Saccharomyces cerevisiae, was curated. This expansion of knowledge and scope has led to a 50% increase in curated complexes compared to the previously published dataset, CYC2008. The yeast complexome is used as a reference resource for the analysis of complexes from large-scale experiments. Our analysis showed that genes coding for proteins in complexes tend to have more genetic interactions, are co-expressed with more genes, are more multifunctional, localize more often in the nucleus, and are more often involved in nucleic acid-related metabolic processes and processes where large machineries are the predominant functional drivers. A comparison to genetic interactions showed that about 40% of expanded co-complex pairs also have genetic interactions, suggesting strong functional links between complex members.

A deep dive into the RecQ interactome: something old and something new

RecQ family helicases are found in all domains of life and play roles in multiple processes that underpin genomic integrity. As such, they are often referred to as guardians or caretakers of the genome. Despite their importance, however, there is still much we do not know about their basic functions in vivo, nor do we fully understand how they interact in organisms that encode more than one RecQ family member. We recently took a multi-omics approach to better understand the Saccharomyces cerevisiae Hrq1 helicase and its interaction with Sgs1, with these enzymes being the functional homologs of the disease-linked RECQL4 and BLM helicases, respectively. Using synthetic genetic array analyses, immuno-precipitation coupled to mass spectrometry, and RNA-seq, we found that Hrq1 and Sgs1 likely participate in many pathways outside of the canonical DNA recombination and repair functions for which they are already known. For instance, connections to transcription, ribosome biogenesis, and chromatin/chromosome organization were uncovered. These recent results are briefly detailed with respect to current knowledge in the field, and possible follow-up experiments are suggested. In this way, we hope to gain a wholistic understanding of these RecQ helicases and how their mutation leads to genomic instability.

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.

Alternative glycosylation controls endoplasmic reticulum dynamics and tubular extension in mammalian cells

The endoplasmic reticulum (ER) is a central eukaryotic organelle with a tubular network made of hairpin proteins linked by hydrolysis of guanosine triphosphate nucleotides. Among posttranslational modifications initiated at the ER level, glycosylation is the most common reaction. However, our understanding of the impact of glycosylation on the ER structure remains unclear. Here, we show that exostosin-1 (EXT1) glycosyltransferase, an enzyme involved in N-glycosylation, is a key regulator of ER morphology and dynamics. We have integrated multiomics and superresolution imaging to characterize the broad effect of EXT1 inactivation, including the ER shape-dynamics-function relationships in mammalian cells. We have observed that inactivating EXT1 induces cell enlargement and enhances metabolic switches such as protein secretion. In particular, suppressing EXT1 in mouse thymocytes causes developmental dysfunctions associated with the ER network extension. Last, our data illuminate the physical and functional aspects of the ER proteome-glycome-lipidome structure axis, with implications in biotechnology and medicine.

Identification and Characterization of Rapidly Accumulating sch9Δ Suppressor Mutations in Saccharomyces cerevisiae

Nutrient sensing is important for cell growth, aging, and longevity. In Saccharomyces cerevisiae, Sch9, an AGC-family protein kinase, is a major nutrient sensing kinase homologous to mammalian Akt and S6 kinase. Sch9 integrates environmental cues with cell growth by functioning downstream of TORC1 and in parallel with the Ras/PKA pathway. Mutations in SCH9 lead to reduced cell growth in dextrose medium; however, reports on the ability of sch9Δ mutants to utilize non-fermentable carbon sources are inconsistent. Here we show that sch9Δ mutant strains cannot grow on non-fermentable carbon sources and rapidly accumulate suppressor mutations, which reverse growth defects of sch9Δ mutants. sch9Δ induces gene expression of three transcription factors required for utilization of non-fermentable carbon sources, Cat8, Adr1, and Hap4, while sch9Δ suppressor mutations, termed sns1 and sns2, strongly decrease the gene expression of those transcription factors. Despite the genetic suppression interactions, both sch9Δ and sns1 (or sns2) homozygous mutants have severe defects in meiosis. By screening mutants defective in sporulation, we identified additional sch9Δ suppressor mutants with mutations in GPB1, GPB2, and MCK1. Using library complementation and genetic analysis, we identified SNS1 and SNS2 to be IRA2 and IRA1, respectively. Furthermore, we discovered that lifespan extension in sch9Δ mutants is dependent on IRA2 and that PKA inactivation greatly increases basal expression of CAT8, ADR1, and HAP4. Our results demonstrate that sch9Δ leads to complete loss of growth on non-fermentable carbon sources and mutations in MCK1 or genes encoding negative regulators of the Ras/PKA pathway reverse sch9Δ mutant phenotypes.

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.

GREB1L as a candidate gene of Mayer-Rokitansky-Küster-Hauser Syndrome

Mayer-Rokitansky-Küster-Hauser (MRKH) Syndrome is a sex development disorder that affects 1 in every 4500 46, XX live births. At least a subset of MRKH syndrome is genetically related to which various candidate genes have been identified. The growth regulation by estrogen in breast cancer 1-like gene (GREB1L) is an androgen-regulated gene reported to be a co-activator of the retinoic acid receptor gene (RAR). Thus expression levels of GREB1L have implications on renal system cellular differentiation, morphogenesis, and homeostasis in vertebrates. Variants of GREB1L have been reported in familial and sporadic MRKH Syndrome and more importantly, in a three-generation family ofMRKH syndrome propositae. Much the same way, Mutants of GREB1L have also been identified in isolated bilateral renal agenesis and deafness both of which are extra-genital tract anomalies in MRKH type 2. Again, renal agenesis transgenic mice have been produced from an E13.5 CRISPR/cas9 GREB1L mutagenesis. Though no GREB1L mutation has been reported in cardiac malformation, there is evidence that GREB1L is involved in ventricular development. Here, we intorigate evidence that projects GREB1L as a candidate gene of Mayer-Rokitansky-Küster-Hauser Syndrome and propose that functional validation analysis to that effect is imparative.

Emergence and propagation of epistasis in metabolic networks

Epistasis is often used to probe functional relationships between genes, and it plays an important role in evolution. However, we lack theory to understand how functional relationships at the molecular level translate into epistasis at the level of whole-organism phenotypes, such as fitness. Here, I derive two rules for how epistasis between mutations with small effects propagates from lower- to higher-level phenotypes in a hierarchical metabolic network with first-order kinetics and how such epistasis depends on topology. Most importantly, weak epistasis at a lower level may be distorted as it propagates to higher levels. Computational analyses show that epistasis in more realistic models likely follows similar, albeit more complex, patterns. These results suggest that pairwise inter-gene epistasis should be common, and it should generically depend on the genetic background and environment. Furthermore, the epistasis coefficients measured for high-level phenotypes may not be sufficient to fully infer the underlying functional relationships.

τ-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.