Proteome-scale movements and compartment connectivity during the eukaryotic cell cycle.

Cell cycle progression relies on coordinated changes in the composition and subcellular localization of the proteome. By applying two distinct convolutional neural networks on images of millions of live yeast cells, we resolved proteome-level dynamics in both concentration and localization during the cell cycle, with resolution of ∼20 subcellular localization classes. We show that a quarter of the proteome displays cell cycle periodicity, with proteins tending to be controlled either at the level of localization or concentration, but not both. Distinct levels of protein regulation are preferentially utilized for different aspects of the cell cycle, with changes in protein concentration being mostly involved in cell cycle control and changes in protein localization in the biophysical implementation of the cell cycle program. We present a resource for exploring global proteome dynamics during the cell cycle, which will aid in understanding a fundamental biological process at a systems level.

The Parkinson's disease protein alpha-synuclein is a modulator of processing bodies and mRNA stability.

Alpha-synuclein (αS) is a conformationally plastic protein that reversibly binds to cellular membranes. It aggregates and is genetically linked to Parkinson's disease (PD). Here, we show that αS directly modulates processing bodies (P-bodies), membraneless organelles that function in mRNA turnover and storage. The N terminus of αS, but not other synucleins, dictates mutually exclusive binding either to cellular membranes or to P-bodies in the cytosol. αS associates with multiple decapping proteins in close proximity on the Edc4 scaffold. As αS pathologically accumulates, aberrant interaction with Edc4 occurs at the expense of physiologic decapping-module interactions. mRNA decay kinetics within PD-relevant pathways are correspondingly disrupted in PD patient neurons and brain. Genetic modulation of P-body components alters αS toxicity, and human genetic analysis lends support to the disease-relevance of these interactions. Beyond revealing an unexpected aspect of αS function and pathology, our data highlight the versatility of conformationally plastic proteins with high intrinsic disorder.

Translocon Declogger Ste24 Protects against IAPP Oligomer-Induced Proteotoxicity.

Aggregates of human islet amyloid polypeptide (IAPP) in the pancreas of patients with type 2 diabetes (T2D) are thought to contribute to ß cell dysfunction and death. To understand how IAPP harms cells and how this might be overcome, we created a yeast model of IAPP toxicity. Ste24, an evolutionarily conserved protease that was recently reported to degrade peptides stuck within the translocon between the cytoplasm and the endoplasmic reticulum, was the strongest suppressor of IAPP toxicity. By testing variants of the human homolog, ZMPSTE24, with varying activity levels, the rescue of IAPP toxicity proved to be directly proportional to the declogging efficiency. Clinically relevant ZMPSTE24 variants identified in the largest database of exomes sequences derived from T2D patients were characterized using the yeast model, revealing 14 partial loss-of-function variants, which were enriched among diabetes patients over 2-fold. Thus, clogging of the translocon by IAPP oligomers may contribute to ß cell failure.

Unraveling the biology of a fungal meningitis pathogen using chemical genetics.

The fungal meningitis pathogen Cryptococcus neoformans is a central driver of mortality in HIV/AIDS. We report a genome-scale chemical genetic data map for this pathogen that quantifies the impact of 439 small-molecule challenges on 1,448 gene knockouts. We identified chemical phenotypes for 83% of mutants screened and at least one genetic response for each compound. C. neoformans chemical-genetic responses are largely distinct from orthologous published profiles of Saccharomyces cerevisiae, demonstrating the importance of pathogen-centered studies. We used the chemical-genetic matrix to predict novel pathogenicity genes, infer compound mode of action, and to develop an algorithm, O2M, that predicts antifungal synergies. These predictions were experimentally validated, thereby identifying virulence genes, a molecule that triggers G2/M arrest and inhibits the Cdc25 phosphatase, and many compounds that synergize with the antifungal drug fluconazole. Our work establishes a chemical-genetic foundation for approaching an infection responsible for greater than one-third of AIDS-related deaths.

Repression of essential cell cycle genes increases cellular fitness.

A network of transcription factors (TFs) coordinates transcription with cell cycle events in eukaryotes. Most TFs in the network are phosphorylated by cyclin-dependent kinase (CDK), which limits their activities during the cell cycle. Here, we investigate the physiological consequences of disrupting CDK regulation of the paralogous repressors Yhp1 and Yox1 in yeast. Blocking Yhp1/Yox1 phosphorylation increases their levels and decreases expression of essential cell cycle regulatory genes which, unexpectedly, increases cellular fitness in optimal growth conditions. Using synthetic genetic interaction screens, we find that Yhp1/Yox1 mutations improve the fitness of mutants with mitotic defects, including condensin mutants. Blocking Yhp1/Yox1 phosphorylation simultaneously accelerates the G1/S transition and delays mitotic exit, without decreasing proliferation rate. This mitotic delay partially reverses the chromosome segregation defect of condensin mutants, potentially explaining their increased fitness when combined with Yhp1/Yox1 phosphomutants. These findings reveal how altering expression of cell cycle genes leads to a redistribution of cell cycle timing and confers a fitness advantage to cells.

Systems-level transcriptional regulation of Caenorhabditis elegans metabolism.

Metabolism is controlled to ensure organismal development and homeostasis. Several mechanisms regulate metabolism, including allosteric control and transcriptional regulation of metabolic enzymes and transporters. So far, metabolism regulation has mostly been described for individual genes and pathways, and the extent of transcriptional regulation of the entire metabolic network remains largely unknown. Here, we find that three-quarters of all metabolic genes are transcriptionally regulated in the nematode Caenorhabditis elegans. We find that many annotated metabolic pathways are coexpressed, and we use gene expression data and the iCEL1314 metabolic network model to define coregulated subpathways in an unbiased manner. Using a large gene expression compendium, we determine the conditions where subpathways exhibit strong coexpression. Finally, we develop "WormClust," a web application that enables a gene-by-gene query of genes to view their association with metabolic (sub)-pathways. Overall, this study sheds light on the ubiquity of transcriptional regulation of metabolism and provides a blueprint for similar studies in other organisms, including humans.

Dimensionality reduction methods for extracting functional networks from large-scale CRISPR screens.

CRISPR-Cas9 screens facilitate the discovery of gene functional relationships and phenotype-specific dependencies. The Cancer Dependency Map (DepMap) is the largest compendium of whole-genome CRISPR screens aimed at identifying cancer-specific genetic dependencies across human cell lines. A mitochondria-associated bias has been previously reported to mask signals for genes involved in other functions, and thus, methods for normalizing this dominant signal to improve co-essentiality networks are of interest. In this study, we explore three unsupervised dimensionality reduction methods-autoencoders, robust, and classical principal component analyses (PCA)-for normalizing the DepMap to improve functional networks extracted from these data. We propose a novel "onion" normalization technique to combine several normalized data layers into a single network. Benchmarking analyses reveal that robust PCA combined with onion normalization outperforms existing methods for normalizing the DepMap. Our work demonstrates the value of removing low-dimensional signals from the DepMap before constructing functional gene networks and provides generalizable dimensionality reduction-based normalization tools.

Skp, Cullin, F-box (SCF)-Met30 and SCF-Cdc4-Mediated Proteolysis of CENP-A Prevents Mislocalization of CENP-A for Chromosomal Stability in Budding Yeast.

Restricting the localization of the histone H3 variant CENP-A (Cse4 in yeast, CID in flies) to centromeres is essential for faithful chromosome segregation. Mislocalization of CENP-A leads to chromosomal instability (CIN) in yeast, fly and human cells. Overexpression and mislocalization of CENP-A has been observed in many cancers and this correlates with increased invasiveness and poor prognosis. Yet genes that regulate CENP-A levels and localization under physiological conditions have not been defined. In this study we used a genome-wide genetic screen to identify essential genes required for Cse4 homeostasis to prevent its mislocalization for chromosomal stability. We show that two Skp, Cullin, F-box (SCF) ubiquitin ligases with the evolutionarily conserved F-box proteins Met30 and Cdc4 interact and cooperatively regulate proteolysis of endogenous Cse4 and prevent its mislocalization for faithful chromosome segregation under physiological conditions. The interaction of Met30 with Cdc4 is independent of the D domain, which is essential for their homodimerization and ubiquitination of other substrates. The requirement for both Cdc4 and Met30 for ubiquitination is specifc for Cse4; and a common substrate for Cdc4 and Met30 has not previously been described. Met30 is necessary for the interaction between Cdc4 and Cse4, and defects in this interaction lead to stabilization and mislocalization of Cse4, which in turn contributes to CIN. We provide the first direct link between Cse4 mislocalization to defects in kinetochore structure and show that SCF-mediated proteolysis of Cse4 is a major mechanism that prevents stable maintenance of Cse4 at non-centromeric regions, thus ensuring faithful chromosome segregation. In summary, we have identified essential pathways that regulate cellular levels of endogenous Cse4 and shown that proteolysis of Cse4 by SCF-Met30/Cdc4 prevents mislocalization and CIN in unperturbed cells.

A method for benchmarking genetic screens reveals a predominant mitochondrial bias.

We present FLEX (Functional evaluation of experimental perturbations), a pipeline that leverages several functional annotation resources to establish reference standards for benchmarking human genome-wide CRISPR screen data and methods for analyzing them. FLEX provides a quantitative measurement of the functional information captured by a given gene-pair dataset and a means to explore the diversity of functions captured by the input dataset. We apply FLEX to analyze data from the diverse cell line screens generated by the DepMap project. We identify a predominant mitochondria-associated signal within co-essentiality networks derived from these data and explore the basis of this signal. Our analysis and time-resolved CRISPR screens in a single cell line suggest that the variable phenotypes associated with mitochondria genes across cells may reflect screen dynamics and protein stability effects rather than genetic dependencies. We characterize this functional bias and demonstrate its relevance for interpreting differential hits in any CRISPR screening context. More generally, we demonstrate the utility of the FLEX pipeline for performing robust comparative evaluations of CRISPR screens or methods for processing them.

EXO1 resection at G-quadruplex structures facilitates resolution and replication.

G-quadruplexes represent unique roadblocks to DNA replication, which tends to stall at these secondary structures. Although G-quadruplexes can be found throughout the genome, telomeres, due to their G-richness, are particularly predisposed to forming these structures and thus represent difficult-to-replicate regions. Here, we demonstrate that exonuclease 1 (EXO1) plays a key role in the resolution of, and replication through, telomeric G-quadruplexes. When replication forks encounter G-quadruplexes, EXO1 resects the nascent DNA proximal to these structures to facilitate fork progression and faithful replication. In the absence of EXO1, forks accumulate at stabilized G-quadruplexes and ultimately collapse. These collapsed forks are preferentially repaired via error-prone end joining as depletion of EXO1 diverts repair away from error-free homology-dependent repair. Such aberrant repair leads to increased genomic instability, which is exacerbated at chromosome termini in the form of dysfunction and telomere loss.

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.

Integrating Coexpression Networks with GWAS to Prioritize Causal Genes in Maize.

Genome-wide association studies (GWAS) have identified loci linked to hundreds of traits in many different species. Yet, because linkage equilibrium implicates a broad region surrounding each identified locus, the causal genes often remain unknown. This problem is especially pronounced in nonhuman, nonmodel species, where functional annotations are sparse and there is frequently little information available for prioritizing candidate genes. We developed a computational approach, Camoco, that integrates loci identified by GWAS with functional information derived from gene coexpression networks. Using Camoco, we prioritized candidate genes from a large-scale GWAS examining the accumulation of 17 different elements in maize (Zea mays) seeds. Strikingly, we observed a strong dependence in the performance of our approach based on the type of coexpression network used: expression variation across genetically diverse individuals in a relevant tissue context (in our case, roots that are the primary elemental uptake and delivery system) outperformed other alternative networks. Two candidate genes identified by our approach were validated using mutants. Our study demonstrates that coexpression networks provide a powerful basis for prioritizing candidate causal genes from GWAS loci but suggests that the success of such strategies can highly depend on the gene expression data context. Both the software and the lessons on integrating GWAS data with coexpression networks generalize to species beyond maize.

Improving Measures of Chemical Structural Similarity Using Machine Learning on Chemical-Genetic Interactions.

A common strategy for identifying molecules likely to possess a desired biological activity is to search large databases of compounds for high structural similarity to a query molecule that demonstrates this activity, under the assumption that structural similarity is predictive of similar biological activity. However, efforts to systematically benchmark the diverse array of available molecular fingerprints and similarity coefficients have been limited by a lack of large-scale datasets that reflect biological similarities of compounds. To elucidate the relative performance of these alternatives, we systematically benchmarked 11 different molecular fingerprint encodings, each combined with 13 different similarity coefficients, using a large set of chemical-genetic interaction data from the yeast Saccharomyces cerevisiae as a systematic proxy for biological activity. We found that the performance of different molecular fingerprints and similarity coefficients varied substantially and that the all-shortest path fingerprints paired with the Braun-Blanquet similarity coefficient provided superior performance that was robust across several compound collections. We further proposed a machine learning pipeline based on support vector machines that offered a fivefold improvement relative to the best unsupervised approach. Our results generally suggest that using high-dimensional chemical-genetic data as a basis for refining molecular fingerprints can be a powerful approach for improving prediction of biological functions from chemical structures.

Extensive rewiring and complex evolutionary dynamics in a C. elegans multiparameter transcription factor network.

Gene duplication results in two identical paralogs that diverge through mutation, leading to loss or gain of interactions with other biomolecules. Here, we comprehensively characterize such network rewiring for C. elegans transcription factors (TFs) within and across four newly delineated molecular networks. Remarkably, we find that even highly similar TFs often have different interaction degrees and partners. In addition, we find that most TF families have a member that is highly connected in multiple networks. Further, different TF families have opposing correlations between network connectivity and phylogenetic age, suggesting that they are subject to different evolutionary pressures. Finally, TFs that have similar partners in one network generally do not in another, indicating a lack of pressure to retain cross-network similarity. Our multiparameter analyses provide unique insights into the evolutionary dynamics that shaped TF networks.

Publisher Correction: Integrating yeast chemical genomics and mammalian cell pathway analysis.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Pathway-based discovery of genetic interactions in breast cancer.

Breast cancer is the second largest cause of cancer death among U.S. women and the leading cause of cancer death among women worldwide. Genome-wide association studies (GWAS) have identified several genetic variants associated with susceptibility to breast cancer, but these still explain less than half of the estimated genetic contribution to the disease. Combinations of variants (i.e. genetic interactions) may play an important role in breast cancer susceptibility. However, due to a lack of statistical power, the current tests for genetic interactions from GWAS data mainly leverage prior knowledge to focus on small sets of genes or SNPs that are known to have an association with breast cancer. Thus, many genetic interactions, particularly among novel variants, remain understudied. Reverse-genetic interaction screens in model organisms have shown that genetic interactions frequently cluster into highly structured motifs, where members of the same pathway share similar patterns of genetic interactions. Based on this key observation, we recently developed a method called BridGE to search for such structured motifs in genetic networks derived from GWAS studies and identify pathway-level genetic interactions in human populations. We applied BridGE to six independent breast cancer cohorts and identified significant pathway-level interactions in five cohorts. Joint analysis across all five cohorts revealed a high confidence consensus set of genetic interactions with support in multiple cohorts. The discovered interactions implicated the glutathione conjugation, vitamin D receptor, purine metabolism, mitotic prometaphase, and steroid hormone biosynthesis pathways as major modifiers of breast cancer risk. Notably, while many of the pathways identified by BridGE show clear relevance to breast cancer, variants in these pathways had not been previously discovered by traditional single variant association tests, or single pathway enrichment analysis that does not consider SNP-SNP interactions.

The highly buffered Arabidopsis immune signaling network conceals the functions of its components.

Plant immunity protects plants from numerous potentially pathogenic microbes. The biological network that controls plant inducible immunity must function effectively even when network components are targeted and disabled by pathogen effectors. Network buffering could confer this resilience by allowing different parts of the network to compensate for loss of one another's functions. Networks rich in buffering rely on interactions within the network, but these mechanisms are difficult to study by simple genetic means. Through a network reconstitution strategy, in which we disassemble and stepwise reassemble the plant immune network that mediates Pattern-Triggered-Immunity, we have resolved systems-level regulatory mechanisms underlying the Arabidopsis transcriptome response to the immune stimulant flagellin-22 (flg22). These mechanisms show widespread evidence of interactions among major sub-networks-we call these sectors-in the flg22-responsive transcriptome. Many of these interactions result in network buffering. Resolved regulatory mechanisms show unexpected patterns for how the jasmonate (JA), ethylene (ET), phytoalexin-deficient 4 (PAD4), and salicylate (SA) signaling sectors control the transcriptional response to flg22. We demonstrate that many of the regulatory mechanisms we resolved are not detectable by the traditional genetic approach of single-gene null-mutant analysis. Similar to potential pathogenic perturbations, null-mutant effects on immune signaling can be buffered by the network.

MOSAIC: a chemical-genetic interaction data repository and web resource for exploring chemical modes of action.

Summary: Chemical-genomic approaches that map interactions between small molecules and genetic perturbations offer a promising strategy for functional annotation of uncharacterized bioactive compounds. We recently developed a new high-throughput platform for mapping chemical-genetic (CG) interactions in yeast that can be scaled to screen large compound collections, and we applied this system to generate CG interaction profiles for more than 13 000 compounds. When integrated with the existing global yeast genetic interaction network, CG interaction profiles can enable mode-of-action prediction for previously uncharacterized compounds as well as discover unexpected secondary effects for known drugs. To facilitate future analysis of these valuable data, we developed a public database and web interface named MOSAIC. The website provides a convenient interface for querying compounds, bioprocesses (Gene Ontology terms) and genes for CG information including direct CG interactions, bioprocesses and gene-level target predictions. MOSAIC also provides access to chemical structure information of screened molecules, chemical-genomic profiles and the ability to search for compounds sharing structural and functional similarity. This resource will be of interest to chemical biologists for discovering new small molecule probes with specific modes-of-action as well as computational biologists interested in analysing CG interaction networks. Availability and implementation: MOSAIC is available at http://mosaic.cs.umn.edu. Contact: hisyo@riken.jp, yoshidam@riken.jp, charlie.boone@utoronto.ca or chadm@umn.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

Different Modes of Negative Regulation of Plant Immunity by Calmodulin-Related Genes.

Plant immune responses activated through the perception of microbe-associated molecular patterns, leading to pattern-triggered immunity, are tightly regulated. This results in low immune responses in the absence of pathogens and a rapid return to the resting state following an activation event. Here, we show that two CALMODULIN-LIKE genes, CML46 and CML47, negatively regulate salicylic acid accumulation and immunity in Arabidopsis (Arabidopsis thaliana). The double mutant cml46 cml47 is highly resistant to the pathogen Pseudomonas syringae pv maculicola (Pma). The effects of cml46 cml47 on Pma growth are genetically additive to that of cbp60a, a known negative regulator in the CALMODULIN-BINDING PROTEIN60 (CBP60) family. Transcriptome profiling revealed the effects of cbp60a and cml46 cml47 on both common and separate sets of genes, with the majorities of these differentially expressed genes being Pma responsive. CBP60g, a positive regulator of immunity in the CBP60 family, was found to be transcriptionally regulated by CBP60a, CML46, and CML47 Analysis of the flg22-induced mRNA levels of CBP60g in cbp60a and cml46 cml47 revealed that cml46 cml47 plants have higher induced expression while cbp60a plants retain elevated levels longer than wild-type plants. Assays for the effect of flg22 treatment on Pma growth showed that the effect is stronger in cml46 cml47 plants and lasts longer in cbp60a plants. Thus, the expression pattern of CBP60g is reflected in flg22-induced resistance to Pma.

Functional annotation of chemical libraries across diverse biological processes.

Chemical-genetic approaches offer the potential for unbiased functional annotation of chemical libraries. Mutations can alter the response of cells in the presence of a compound, revealing chemical-genetic interactions that can elucidate a compound's mode of action. We developed a highly parallel, unbiased yeast chemical-genetic screening system involving three key components. First, in a drug-sensitive genetic background, we constructed an optimized diagnostic mutant collection that is predictive for all major yeast biological processes. Second, we implemented a multiplexed (768-plex) barcode-sequencing protocol, enabling the assembly of thousands of chemical-genetic profiles. Finally, based on comparison of the chemical-genetic profiles with a compendium of genome-wide genetic interaction profiles, we predicted compound functionality. Applying this high-throughput approach, we screened seven different compound libraries and annotated their functional diversity. We further validated biological process predictions, prioritized a diverse set of compounds, and identified compounds that appear to have dual modes of action.