Widespread variation in molecular interactions and regulatory properties among transcription factor isoforms.

Most human Transcription factors (TFs) genes encode multiple protein isoforms differing in DNA binding domains, effector domains, or other protein regions. The global extent to which this results in functional differences between isoforms remains unknown. Here, we systematically compared 693 isoforms of 246 TF genes, assessing DNA binding, protein binding, transcriptional activation, subcellular localization, and condensate formation. Relative to reference isoforms, two-thirds of alternative TF isoforms exhibit differences in one or more molecular activities, which often could not be predicted from sequence. We observed two primary categories of alternative TF isoforms: "rewirers" and "negative regulators", both of which were associated with differentiation and cancer. Our results support a model wherein the relative expression levels of, and interactions involving, TF isoforms add an understudied layer of complexity to gene regulatory networks, demonstrating the importance of isoform-aware characterization of TF functions and providing a rich resource for further studies.

Systematically testing human HMBS missense variants to reveal mechanism and pathogenic variation.

Defects in hydroxymethylbilane synthase (HMBS) can cause acute intermittent porphyria (AIP), an acute neurological disease. Although sequencing-based diagnosis can be definitive, ∼⅓ of clinical HMBS variants are missense variants, and most clinically reported HMBS missense variants are designated as "variants of uncertain significance" (VUSs). Using saturation mutagenesis, en masse selection, and sequencing, we applied a multiplexed validated assay to both the erythroid-specific and ubiquitous isoforms of HMBS, obtaining confident functional impact scores for >84% of all possible amino acid substitutions. The resulting variant effect maps generally agreed with biochemical expectations and provide further evidence that HMBS can function as a monomer. Additionally, the maps implicated specific residues as having roles in active site dynamics, which was further supported by molecular dynamics simulations. Most importantly, these maps can help discriminate pathogenic from benign HMBS variants, proactively providing evidence even for yet-to-be-observed clinical missense variants.

A comprehensive map of human glucokinase variant activity.

BACKGROUND: Glucokinase (GCK) regulates insulin secretion to maintain appropriate blood glucose levels. Sequence variants can alter GCK activity to cause hyperinsulinemic hypoglycemia or hyperglycemia associated with GCK-maturity-onset diabetes of the young (GCK-MODY), collectively affecting up to 10 million people worldwide. Patients with GCK-MODY are frequently misdiagnosed and treated unnecessarily. Genetic testing can prevent this but is hampered by the challenge of interpreting novel missense variants. RESULT: Here, we exploit a multiplexed yeast complementation assay to measure both hyper- and hypoactive GCK variation, capturing 97% of all possible missense and nonsense variants. Activity scores correlate with in vitro catalytic efficiency, fasting glucose levels in carriers of GCK variants and with evolutionary conservation. Hypoactive variants are concentrated at buried positions, near the active site, and at a region of known importance for GCK conformational dynamics. Some hyperactive variants shift the conformational equilibrium towards the active state through a relative destabilization of the inactive conformation. CONCLUSION: Our comprehensive assessment of GCK variant activity promises to facilitate variant interpretation and diagnosis, expand our mechanistic understanding of hyperactive variants, and inform development of therapeutics targeting GCK.

Systematically testing human HMBS missense variants to reveal mechanism and pathogenic variation.

Defects in hydroxymethylbilane synthase (HMBS) can cause Acute Intermittent Porphyria (AIP), an acute neurological disease. Although sequencing-based diagnosis can be definitive, ~⅓ of clinical HMBS variants are missense variants, and most clinically-reported HMBS missense variants are designated as "variants of uncertain significance" (VUS). Using saturation mutagenesis, en masse selection, and sequencing, we applied a multiplexed validated assay to both the erythroid-specific and ubiquitous isoforms of HMBS, obtaining confident functional impact scores for >84% of all possible amino-acid substitutions. The resulting variant effect maps generally agreed with biochemical expectation. However, the maps showed variants at the dimerization interface to be unexpectedly well tolerated, and suggested residue roles in active site dynamics that were supported by molecular dynamics simulations. Most importantly, these HMBS variant effect maps can help discriminate pathogenic from benign variants, proactively providing evidence even for yet-to-be-observed clinical missense variants.

Shifting landscapes of human MTHFR missense-variant effects.

Most rare clinical missense variants cannot currently be classified as pathogenic or benign. Deficiency in human 5,10-methylenetetrahydrofolate reductase (MTHFR), the most common inherited disorder of folate metabolism, is caused primarily by rare missense variants. Further complicating variant interpretation, variant impacts often depend on environment. An important example of this phenomenon is the MTHFR variant p.Ala222Val (c.665C>T), which is carried by half of all humans and has a phenotypic impact that depends on dietary folate. Here we describe the results of 98,336 variant functional-impact assays, covering nearly all possible MTHFR amino acid substitutions in four folinate environments, each in the presence and absence of p.Ala222Val. The resulting atlas of MTHFR variant effects reveals many complex dependencies on both folinate and p.Ala222Val. MTHFR atlas scores can distinguish pathogenic from benign variants and, among individuals with severe MTHFR deficiency, correlate with age of disease onset. Providing a powerful tool for understanding structure-function relationships, the atlas suggests a role for a disordered loop in retaining cofactor at the active site and identifies variants that enable escape of inhibition by S-adenosylmethionine. Thus, a model based on eight MTHFR variant effect maps illustrates how shifting landscapes of environment- and genetic-background-dependent missense variation can inform our clinical, structural, and functional understanding of MTHFR deficiency.

Systematic analysis of bypass suppression of essential genes.

Essential genes tend to be highly conserved across eukaryotes, but, in some cases, their critical roles can be bypassed through genetic rewiring. From a systematic analysis of 728 different essential yeast genes, we discovered that 124 (17%) were dispensable essential genes. Through whole-genome sequencing and detailed genetic analysis, we investigated the genetic interactions and genome alterations underlying bypass suppression. Dispensable essential genes often had paralogs, were enriched for genes encoding membrane-associated proteins, and were depleted for members of protein complexes. Functionally related genes frequently drove the bypass suppression interactions. These gene properties were predictive of essential gene dispensability and of specific suppressors among hundreds of genes on aneuploid chromosomes. Our findings identify yeast's core essential gene set and reveal that the properties of dispensable essential genes are conserved from yeast to human cells, correlating with human genes that display cell line-specific essentiality in the Cancer Dependency Map (DepMap) project.

A reference map of the human binary protein interactome.

Global insights into cellular organization and genome function require comprehensive understanding of the interactome networks that mediate genotype-phenotype relationships1,2. Here we present a human 'all-by-all' reference interactome map of human binary protein interactions, or 'HuRI'. With approximately 53,000 protein-protein interactions, HuRI has approximately four times as many such interactions as there are high-quality curated interactions from small-scale studies. The integration of HuRI with genome3, transcriptome4 and proteome5 data enables cellular function to be studied within most physiological or pathological cellular contexts. We demonstrate the utility of HuRI in identifying the specific subcellular roles of protein-protein interactions. Inferred tissue-specific networks reveal general principles for the formation of cellular context-specific functions and elucidate potential molecular mechanisms that might underlie tissue-specific phenotypes of Mendelian diseases. HuRI is a systematic proteome-wide reference that links genomic variation to phenotypic outcomes.

Highly Combinatorial Genetic Interaction Analysis Reveals a Multi-Drug Transporter Influence Network.

Many traits are complex, depending non-additively on variant combinations. Even in model systems, such as the yeast S. cerevisiae, carrying out the high-order variant-combination testing needed to dissect complex traits remains a daunting challenge. Here, we describe "X-gene" genetic analysis (XGA), a strategy for engineering and profiling highly combinatorial gene perturbations. We demonstrate XGA on yeast ABC transporters by engineering 5,353 strains, each deleted for a random subset of 16 transporters, and profiling each strain's resistance to 16 compounds. XGA yielded 85,648 genotype-to-resistance observations, revealing high-order genetic interactions for 13 of the 16 transporters studied. Neural networks yielded intuitive functional models and guided exploration of fluconazole resistance, which was influenced non-additively by five genes. Together, our results showed that highly combinatorial genetic perturbation can functionally dissect complex traits, supporting pursuit of analogous strategies in human cells and other model systems.

Mapping DNA damage-dependent genetic interactions in yeast via party mating and barcode fusion genetics.

Condition-dependent genetic interactions can reveal functional relationships between genes that are not evident under standard culture conditions. State-of-the-art yeast genetic interaction mapping, which relies on robotic manipulation of arrays of double-mutant strains, does not scale readily to multi-condition studies. Here, we describe barcode fusion genetics to map genetic interactions (BFG-GI), by which double-mutant strains generated via en masse "party" mating can also be monitored en masse for growth to detect genetic interactions. By using site-specific recombination to fuse two DNA barcodes, each representing a specific gene deletion, BFG-GI enables multiplexed quantitative tracking of double mutants via next-generation sequencing. We applied BFG-GI to a matrix of DNA repair genes under nine different conditions, including methyl methanesulfonate (MMS), 4-nitroquinoline 1-oxide (4NQO), bleomycin, zeocin, and three other DNA-damaging environments. BFG-GI recapitulated known genetic interactions and yielded new condition-dependent genetic interactions. We validated and further explored a subnetwork of condition-dependent genetic interactions involving MAG1, SLX4, and genes encoding the Shu complex, and inferred that loss of the Shu complex leads to an increase in the activation of the checkpoint protein kinase Rad53.

CRISPR/Cas9 System as a Valuable Genome Editing Tool for Wine Yeasts with Application to Decrease Urea Production.

An extensive repertoire of molecular tools is available for genetic analysis in laboratory strains of S. cerevisiae. Although this has widely contributed to the interpretation of gene functionality within haploid laboratory isolates, the genetics of metabolism in commercially-relevant polyploid yeast strains is still poorly understood. Genetic engineering in industrial yeasts is undergoing major changes due to Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) and CRISPR-associated protein (Cas) engineering approaches. Here we apply the CRISPR/Cas9 system to two commercial "starter" strains of S. cerevisiae (EC1118, AWRI796), eliminating the CAN1 arginine permease pathway to generate strains with reduced urea production (18.5 and 35.5% for EC1118 and AWRI796, respectively). In a wine-model environment based on two grape musts obtained from Chardonnay and Cabernet Sauvignon cultivars, both S. cerevisiae starter strains and CAN1 mutants completed the must fermentation in 8-12 days. However, recombinant strains carrying the can1 mutation failed to produce urea, suggesting that the genetic modification successfully impaired the arginine metabolism. In conclusion, the reduction of urea production in a wine-model environment confirms that the CRISPR/Cas9 system has been successfully established in S. cerevisiae wine yeasts.

Open Source Drug Discovery: Highly Potent Antimalarial Compounds Derived from the Tres Cantos Arylpyrroles.

The development of new antimalarial compounds remains a pivotal part of the strategy for malaria elimination. Recent large-scale phenotypic screens have provided a wealth of potential starting points for hit-to-lead campaigns. One such public set is explored, employing an open source research mechanism in which all data and ideas were shared in real time, anyone was able to participate, and patents were not sought. One chemical subseries was found to exhibit oral activity but contained a labile ester that could not be replaced without loss of activity, and the original hit exhibited remarkable sensitivity to minor structural change. A second subseries displayed high potency, including activity within gametocyte and liver stage assays, but at the cost of low solubility. As an open source research project, unexplored avenues are clearly identified and may be explored further by the community; new findings may be cumulatively added to the present work.

Exploring genetic suppression interactions on a global scale.

Genetic suppression occurs when the phenotypic defects caused by a mutation in a particular gene are rescued by a mutation in a second gene. To explore the principles of genetic suppression, we examined both literature-curated and unbiased experimental data, involving systematic genetic mapping and whole-genome sequencing, to generate a large-scale suppression network among yeast genes. Most suppression pairs identified novel relationships among functionally related genes, providing new insights into the functional wiring diagram of the cell. In addition to suppressor mutations, we identified frequent secondary mutations,in a subset of genes, that likely cause a delay in the onset of stationary phase, which appears to promote their enrichment within a propagating population. These findings allow us to formulate and quantify general mechanisms of genetic suppression.

Pooled-matrix protein interaction screens using Barcode Fusion Genetics.

High-throughput binary protein interaction mapping is continuing to extend our understanding of cellular function and disease mechanisms. However, we remain one or two orders of magnitude away from a complete interaction map for humans and other major model organisms. Completion will require screening at substantially larger scales with many complementary assays, requiring further efficiency gains in proteome-scale interaction mapping. Here, we report Barcode Fusion Genetics-Yeast Two-Hybrid (BFG-Y2H), by which a full matrix of protein pairs can be screened in a single multiplexed strain pool. BFG-Y2H uses Cre recombination to fuse DNA barcodes from distinct plasmids, generating chimeric protein-pair barcodes that can be quantified via next-generation sequencing. We applied BFG-Y2H to four different matrices ranging in scale from ~25 K to 2.5 M protein pairs. The results show that BFG-Y2H increases the efficiency of protein matrix screening, with quality that is on par with state-of-the-art Y2H methods.

Genome-wide RNAi analysis reveals that simultaneous inhibition of specific mevalonate pathway genes potentiates tumor cell death.

The mevalonate (MVA) pathway is often dysregulated or overexpressed in many cancers suggesting tumor dependency on this classic metabolic pathway. Statins, which target the rate-limiting enzyme of this pathway, 3-hydroxy-3-methylglutaryl-CoA reductase (HMGCR), are promising agents currently being evaluated in clinical trials for anti-cancer efficacy. To uncover novel targets that potentiate statin-induced apoptosis when knocked down, we carried out a pooled genome-wide short hairpin RNA (shRNA) screen. Genes of the MVA pathway were amongst the top-scoring targets, including sterol regulatory element binding transcription factor 2 (SREBP2), 3-hydroxy-3-methylglutaryl-coenzyme A synthase 1 (HMGCS1) and geranylgeranyl diphosphate synthase 1 (GGPS1). Each gene was independently validated and shown to significantly sensitize A549 cells to statin-induced apoptosis when knocked down. SREBP2 knockdown in lung and breast cancer cells completely abrogated the fluvastatin-induced upregulation of sterol-responsive genes HMGCR and HMGCS1. Knockdown of SREBP2 alone did not affect three-dimensional growth of lung and breast cancer cells, yet in combination with fluvastatin cell growth was disrupted. Taken together, these results show that directly targeting multiple levels of the MVA pathway, including blocking the sterol-feedback loop initiated by statin treatment, is an effective and targetable anti-tumor strategy.

Select microtubule inhibitors increase lysosome acidity and promote lysosomal disruption in acute myeloid leukemia (AML) cells.

To identify new biological vulnerabilities in acute myeloid leukemia, we screened a library of natural products for compounds cytotoxic to TEX leukemia cells. This screen identified the novel small molecule Deoxysappanone B 7,4' dimethyl ether (Deox B 7,4), which possessed nanomolar anti-leukemic activity. To determine the anti-leukemic mechanism of action of Deox B 7,4, we conducted a genome-wide screen in Saccharomyces cerevisiae and identified enrichment of genes related to mitotic cell cycle as well as vacuolar acidification, therefore pointing to microtubules and vacuolar (V)-ATPase as potential drug targets. Further investigations into the mechanisms of action of Deox B 7,4 and a related analogue revealed that these compounds were reversible microtubule inhibitors that bound near the colchicine site. In addition, Deox B 7,4 and its analogue increased lysosomal V-ATPase activity and lysosome acidity. The effects on microtubules and lysosomes were functionally important for the anti-leukemic effects of these drugs. The lysosomal effects were characteristic of select microtubule inhibitors as only the Deox compounds and nocodazole, but not colchicine, vinca alkaloids or paclitaxel, altered lysosome acidity and induced lysosomal disruption. Thus, our data highlight a new mechanism of action of select microtubule inhibitors on lysosomal function.

Genes required for survival in microgravity revealed by genome-wide yeast deletion collections cultured during spaceflight.

Spaceflight is a unique environment with profound effects on biological systems including tissue redistribution and musculoskeletal stresses. However, the more subtle biological effects of spaceflight on cells and organisms are difficult to measure in a systematic, unbiased manner. Here we test the utility of the molecularly barcoded yeast deletion collection to provide a quantitative assessment of the effects of microgravity on a model organism. We developed robust hardware to screen, in parallel, the complete collection of ~4800 homozygous and ~5900 heterozygous (including ~1100 single-copy deletions of essential genes) yeast deletion strains, each carrying unique DNA that acts as strain identifiers. We compared strain fitness for the homozygous and heterozygous yeast deletion collections grown in spaceflight and ground, as well as plus and minus hyperosmolar sodium chloride, providing a second additive stressor. The genome-wide sensitivity profiles obtained from these treatments were then queried for their similarity to a compendium of drugs whose effects on the yeast collection have been previously reported. We found that the effects of spaceflight have high concordance with the effects of DNA-damaging agents and changes in redox state, suggesting mechanisms by which spaceflight may negatively affect cell fitness.

Mapping the cellular response to small molecules using chemogenomic fitness signatures.

Genome-wide characterization of the in vivo cellular response to perturbation is fundamental to understanding how cells survive stress. Identifying the proteins and pathways perturbed by small molecules affects biology and medicine by revealing the mechanisms of drug action. We used a yeast chemogenomics platform that quantifies the requirement for each gene for resistance to a compound in vivo to profile 3250 small molecules in a systematic and unbiased manner. We identified 317 compounds that specifically perturb the function of 121 genes and characterized the mechanism of specific compounds. Global analysis revealed that the cellular response to small molecules is limited and described by a network of 45 major chemogenomic signatures. Our results provide a resource for the discovery of functional interactions among genes, chemicals, and biological processes.

A compendium of nucleosome and transcript profiles reveals determinants of chromatin architecture and transcription.

Nucleosomes in all eukaryotes examined to date adopt a characteristic architecture within genes and play fundamental roles in regulating transcription, yet the identity and precise roles of many of the trans-acting factors responsible for the establishment and maintenance of this organization remain to be identified. We profiled a compendium of 50 yeast strains carrying conditional alleles or complete deletions of genes involved in transcriptional regulation, histone biology, and chromatin remodeling, as well as compounds that target transcription and histone deacetylases, to assess their respective roles in nucleosome positioning and transcription. We find that nucleosome patterning in genes is affected by many factors, including the CAF-1 complex, Spt10, and Spt21, in addition to previously reported remodeler ATPases and histone chaperones. Disruption of these factors or reductions in histone levels led genic nucleosomes to assume positions more consistent with their intrinsic sequence preferences, with pronounced and specific shifts of the +1 nucleosome relative to the transcription start site. These shifts of +1 nucleosomes appear to have functional consequences, as several affected genes in Ino80 mutants exhibited altered expression responses. Our parallel expression profiling compendium revealed extensive transcription changes in intergenic and antisense regions, most of which occur in regions with altered nucleosome occupancy and positioning. We show that the nucleosome-excluding transcription factors Reb1, Abf1, Tbf1, and Rsc3 suppress cryptic transcripts at their target promoters, while a combined analysis of nucleosome and expression profiles identified 36 novel transcripts that are normally repressed by Tup1/Cyc8. Our data confirm and extend the roles of chromatin remodelers and chaperones as major determinants of genic nucleosome positioning, and these data provide a valuable resource for future studies.

Lysosomal disruption preferentially targets acute myeloid leukemia cells and progenitors.

Despite efforts to understand and treat acute myeloid leukemia (AML), there remains a need for more comprehensive therapies to prevent AML-associated relapses. To identify new therapeutic strategies for AML, we screened a library of on- and off-patent drugs and identified the antimalarial agent mefloquine as a compound that selectively kills AML cells and AML stem cells in a panel of leukemia cell lines and in mice. Using a yeast genome-wide functional screen for mefloquine sensitizers, we identified genes associated with the yeast vacuole, the homolog of the mammalian lysosome. Consistent with this, we determined that mefloquine disrupts lysosomes, directly permeabilizes the lysosome membrane, and releases cathepsins into the cytosol. Knockdown of the lysosomal membrane proteins LAMP1 and LAMP2 resulted in decreased cell viability, as did treatment of AML cells with known lysosome disrupters. Highlighting a potential therapeutic rationale for this strategy, leukemic cells had significantly larger lysosomes compared with normal cells, and leukemia-initiating cells overexpressed lysosomal biogenesis genes. These results demonstrate that lysosomal disruption preferentially targets AML cells and AML progenitor cells, providing a rationale for testing lysosomal disruption as a novel therapeutic strategy for AML.

Effects of the Paf1 complex and histone modifications on snoRNA 3'-end formation reveal broad and locus-specific regulation.

Across diverse eukaryotes, the Paf1 complex (Paf1C) plays critical roles in RNA polymerase II transcription elongation and regulation of histone modifications. Beyond these roles, the human and Saccharomyces cerevisiae Paf1 complexes also interact with RNA 3'-end processing components to affect transcript 3'-end formation. Specifically, the Saccharomyces cerevisiae Paf1C functions with the RNA binding proteins Nrd1 and Nab3 to regulate the termination of at least two small nucleolar RNAs (snoRNAs). To determine how Paf1C-dependent functions regulate snoRNA formation, we used high-density tiling arrays to analyze transcripts in paf1Δ cells and uncover new snoRNA targets of Paf1. Detailed examination of Paf1-regulated snoRNA genes revealed locus-specific requirements for Paf1-dependent posttranslational histone modifications. We also discovered roles for the transcriptional regulators Bur1-Bur2, Rad6, and Set2 in snoRNA 3'-end formation. Surprisingly, at some snoRNAs, this function of Rad6 appears to be primarily independent of its role in histone H2B monoubiquitylation. Cumulatively, our work reveals a broad requirement for the Paf1C in snoRNA 3'-end formation in S. cerevisiae, implicates the participation of transcriptional proteins and histone modifications in this process, and suggests that the Paf1C contributes to the fine tuning of nuanced levels of regulation that exist at individual loci.