Widespread Expansion of Protein Interaction Capabilities by Alternative Splicing.

While alternative splicing is known to diversify the functional characteristics of some genes, the extent to which protein isoforms globally contribute to functional complexity on a proteomic scale remains unknown. To address this systematically, we cloned full-length open reading frames of alternatively spliced transcripts for a large number of human genes and used protein-protein interaction profiling to functionally compare hundreds of protein isoform pairs. The majority of isoform pairs share less than 50% of their interactions. In the global context of interactome network maps, alternative isoforms tend to behave like distinct proteins rather than minor variants of each other. Interaction partners specific to alternative isoforms tend to be expressed in a highly tissue-specific manner and belong to distinct functional modules. Our strategy, applicable to other functional characteristics, reveals a widespread expansion of protein interaction capabilities through alternative splicing and suggests that many alternative "isoforms" are functionally divergent (i.e., "functional alloforms").

Active growth signaling promotes senescence and cancer cell sensitivity to CDK7 inhibition.

Tumor growth is driven by continued cellular growth and proliferation. Cyclin-dependent kinase 7's (CDK7) role in activating mitotic CDKs and global gene expression makes it therefore an attractive target for cancer therapies. However, what makes cancer cells particularly sensitive to CDK7 inhibition (CDK7i) remains unclear. Here, we address this question. We show that CDK7i, by samuraciclib, induces a permanent cell-cycle exit, known as senescence, without promoting DNA damage signaling or cell death. A chemogenetic genome-wide CRISPR knockout screen identified that active mTOR (mammalian target of rapamycin) signaling promotes samuraciclib-induced senescence. mTOR inhibition decreases samuraciclib sensitivity, and increased mTOR-dependent growth signaling correlates with sensitivity in cancer cell lines. Reverting a growth-promoting mutation in PIK3CA to wild type decreases sensitivity to CDK7i. Our work establishes that enhanced growth alone promotes CDK7i sensitivity, providing an explanation for why some cancers are more sensitive to CDK inhibition than normally growing cells.

Widespread macromolecular interaction perturbations in human genetic disorders.

How disease-associated mutations impair protein activities in the context of biological networks remains mostly undetermined. Although a few renowned alleles are well characterized, functional information is missing for over 100,000 disease-associated variants. Here we functionally profile several thousand missense mutations across a spectrum of Mendelian disorders using various interaction assays. The majority of disease-associated alleles exhibit wild-type chaperone binding profiles, suggesting they preserve protein folding or stability. While common variants from healthy individuals rarely affect interactions, two-thirds of disease-associated alleles perturb protein-protein interactions, with half corresponding to "edgetic" alleles affecting only a subset of interactions while leaving most other interactions unperturbed. With transcription factors, many alleles that leave protein-protein interactions intact affect DNA binding. Different mutations in the same gene leading to different interaction profiles often result in distinct disease phenotypes. Thus disease-associated alleles that perturb distinct protein activities rather than grossly affecting folding and stability are relatively widespread.

A proteome-scale map of the human interactome network.

Just as reference genome sequences revolutionized human genetics, reference maps of interactome networks will be critical to fully understand genotype-phenotype relationships. Here, we describe a systematic map of ?14,000 high-quality human binary protein-protein interactions. At equal quality, this map is ?30% larger than what is available from small-scale studies published in the literature in the last few decades. While currently available information is highly biased and only covers a relatively small portion of the proteome, our systematic map appears strikingly more homogeneous, revealing a "broader" human interactome network than currently appreciated. The map also uncovers significant interconnectivity between known and candidate cancer gene products, providing unbiased evidence for an expanded functional cancer landscape, while demonstrating how high-quality interactome models will help "connect the dots" of the genomic revolution.

Genome-Wide Screens Reveal that Resveratrol Induces Replicative Stress in Human Cells.

Resveratrol is a natural product associated with wide-ranging effects in animal and cellular models, including lifespan extension. To identify the genetic target of resveratrol in human cells, we conducted genome-wide CRISPR-Cas9 screens to pinpoint genes that confer sensitivity or resistance to resveratrol. An extensive network of DNA damage response and replicative stress genes exhibited genetic interactions with resveratrol and its analog pterostilbene. These genetic profiles showed similarity to the response to hydroxyurea, an inhibitor of ribonucleotide reductase that causes replicative stress. Resveratrol, pterostilbene, and hydroxyurea caused similar depletion of nucleotide pools, inhibition of replication fork progression, and induction of replicative stress. The ability of resveratrol to inhibit cell proliferation and S phase transit was independent of the histone deacetylase sirtuin 1, which has been implicated in lifespan extension by resveratrol. These results establish that a primary impact of resveratrol on human cell proliferation is the induction of low-level replicative stress.

Genetic enhancers of partial PLK1 inhibition reveal hypersensitivity to kinetochore perturbations.

Polo-like kinase 1 (PLK1) is a serine/threonine kinase required for mitosis and cytokinesis. As cancer cells are often hypersensitive to partial PLK1 inactivation, chemical inhibitors of PLK1 have been developed and tested in clinical trials. However, these small molecule inhibitors alone are not completely effective. PLK1 promotes numerous molecular and cellular events in the cell division cycle and it is unclear which of these events most crucially depend on PLK1 activity. We used a CRISPR-based genome-wide screening strategy to identify genes whose inactivation enhances cell proliferation defects upon partial chemical inhibition of PLK1. Genes identified encode proteins that are functionally linked to PLK1 in multiple ways, most notably factors that promote centromere and kinetochore function. Loss of the kinesin KIF18A or the outer kinetochore protein SKA1 in PLK1-compromised cells resulted in mitotic defects, activation of the spindle assembly checkpoint and nuclear reassembly defects. We also show that PLK1-dependent CENP-A loading at centromeres is extremely sensitive to partial PLK1 inhibition. Our results suggest that partial inhibition of PLK1 compromises the integrity and function of the centromere/kinetochore complex, rendering cells hypersensitive to different kinetochore perturbations. We propose that KIF18A is a promising target for combinatorial therapies with PLK1 inhibitors.

The microprotein Nrs1 rewires the G1/S transcriptional machinery during nitrogen limitation in budding yeast.

Commitment to cell division at the end of G1 phase, termed Start in the budding yeast Saccharomyces cerevisiae, is strongly influenced by nutrient availability. To identify new dominant activators of Start that might operate under different nutrient conditions, we screened a genome-wide ORF overexpression library for genes that bypass a Start arrest caused by absence of the G1 cyclin Cln3 and the transcriptional activator Bck2. We recovered a hypothetical gene YLR053c, renamed NRS1 for Nitrogen-Responsive Start regulator 1, which encodes a poorly characterized 108 amino acid microprotein. Endogenous Nrs1 was nuclear-localized, restricted to poor nitrogen conditions, induced upon TORC1 inhibition, and cell cycle-regulated with a peak at Start. NRS1 interacted genetically with SWI4 and SWI6, which encode subunits of the main G1/S transcription factor complex SBF. Correspondingly, Nrs1 physically interacted with Swi4 and Swi6 and was localized to G1/S promoter DNA. Nrs1 exhibited inherent transactivation activity, and fusion of Nrs1 to the SBF inhibitor Whi5 was sufficient to suppress other Start defects. Nrs1 appears to be a recently evolved microprotein that rewires the G1/S transcriptional machinery under poor nitrogen conditions.

The BioGRID interaction database: 2019 update.

The Biological General Repository for Interaction Datasets (BioGRID: https://thebiogrid.org) is an open access database dedicated to the curation and archival storage of protein, genetic and chemical interactions for all major model organism species and humans. As of September 2018 (build 3.4.164), BioGRID contains records for 1 598 688 biological interactions manually annotated from 55 809 publications for 71 species, as classified by an updated set of controlled vocabularies for experimental detection methods. BioGRID also houses records for >700 000 post-translational modification sites. BioGRID now captures chemical interaction data, including chemical-protein interactions for human drug targets drawn from the DrugBank database and manually curated bioactive compounds reported in the literature. A new dedicated aspect of BioGRID annotates genome-wide CRISPR/Cas9-based screens that report gene-phenotype and gene-gene relationships. An extension of the BioGRID resource called the Open Repository for CRISPR Screens (ORCS) database (https://orcs.thebiogrid.org) currently contains over 500 genome-wide screens carried out in human or mouse cell lines. All data in BioGRID is made freely available without restriction, is directly downloadable in standard formats and can be readily incorporated into existing applications via our web service platforms. BioGRID data are also freely distributed through partner model organism databases and meta-databases.

The Sam68 STAR RNA-binding protein regulates mTOR alternative splicing during adipogenesis.

We report that mice ablated for the Sam68 RNA-binding protein exhibit a lean phenotype as a result of increased energy expenditure, decreased commitment to early adipocyte progenitors, and defects in adipogenic differentiation. The Sam68(-/-) mice were protected from obesity, insulin resistance, and glucose intolerance induced with a high-fat diet. To identify the alternative splice events regulated by Sam68, genome-wide exon usage profiling in white adipose tissue was performed. Adipocytes from Sam68(-/-) mice retained intron 5 within the mTOR transcript introducing a premature termination codon, leading to an unstable mRNA. Consequently, Sam68-depleted cells had reduced mTOR levels resulting in lower levels of insulin-stimulated S6 and Akt phosphorylation leading to defects in adipogenesis, and this defect was rescued by the exogenous expression of full-length mTOR. Sam68 bound intronic splice elements within mTOR intron 5 required for the usage of the 5' splice site. We propose that Sam68 regulates alternative splicing during adipogenesis.

RNA-Seq identifies SPGs as a ventral skeletal patterning cue in sea urchins.

The sea urchin larval skeleton offers a simple model for formation of developmental patterns. The calcium carbonate skeleton is secreted by primary mesenchyme cells (PMCs) in response to largely unknown patterning cues expressed by the ectoderm. To discover novel ectodermal cues, we performed an unbiased RNA-Seq-based screen and functionally tested candidates; we thereby identified several novel skeletal patterning cues. Among these, we show that SLC26a2/7 is a ventrally expressed sulfate transporter that promotes a ventral accumulation of sulfated proteoglycans, which is required for ventral PMC positioning and skeletal patterning. We show that the effects of SLC perturbation are mimicked by manipulation of either external sulfate levels or proteoglycan sulfation. These results identify novel skeletal patterning genes and demonstrate that ventral proteoglycan sulfation serves as a positional cue for sea urchin skeletal patterning.

CLMSVault: A Software Suite for Protein Cross-Linking Mass-Spectrometry Data Analysis and Visualization.

Protein cross-linking mass spectrometry (CL-MS) enables the sensitive detection of protein interactions and the inference of protein complex topology. The detection of chemical cross-links between protein residues can identify intra- and interprotein contact sites or provide physical constraints for molecular modeling of protein structure. Recent innovations in cross-linker design, sample preparation, mass spectrometry, and software tools have significantly improved CL-MS approaches. Although a number of algorithms now exist for the identification of cross-linked peptides from mass spectral data, a dearth of user-friendly analysis tools represent a practical bottleneck to the broad adoption of the approach. To facilitate the analysis of CL-MS data, we developed CLMSVault, a software suite designed to leverage existing CL-MS algorithms and provide intuitive and flexible tools for cross-platform data interpretation. CLMSVault stores and combines complementary information obtained from different cross-linkers and search algorithms. CLMSVault provides filtering, comparison, and visualization tools to support CL-MS analyses and includes a workflow for label-free quantification of cross-linked peptides. An embedded 3D viewer enables the visualization of quantitative data and the mapping of cross-linked sites onto PDB structural models. We demonstrate the application of CLMSVault for the analysis of a noncovalent Cdc34-ubiquitin protein complex cross-linked under different conditions. CLMSVault is open-source software (available at https://gitlab.com/courcelm/clmsvault.git ), and a live demo is available at http://democlmsvault.tyerslab.com/ .

The anticancer potential of the CLK kinases inhibitors 1C8 and GPS167 revealed by their impact on the epithelial-mesenchymal transition and the antiviral immune response.

The diheteroarylamide-based compound 1C8 and the aminothiazole carboxamide-related compound GPS167 inhibit the CLK kinases, and affect the proliferation of a broad range of cancer cell lines. A chemogenomic screen previously performed with GPS167 revealed that the depletion of components associated with mitotic spindle assembly altered sensitivity to GPS167. Here, a similar screen performed with 1C8 also established the impact of components involved in mitotic spindle assembly. Accordingly, transcriptome analyses of cells treated with 1C8 and GPS167 indicated that the expression and RNA splicing of transcripts encoding mitotic spindle assembly components were affected. The functional relevance of the microtubule connection was confirmed by showing that subtoxic concentrations of drugs affecting mitotic spindle assembly increased sensitivity to GPS167. 1C8 and GPS167 impacted the expression and splicing of transcripts in pathways relevant to tumor progression, including MYC targets and the epithelial mesenchymal transition (EMT). Finally, 1C8 and GPS167 altered the expression and alternative splicing of transcripts involved in the antiviral immune response. Consistent with this observation, depleting the double-stranded RNA sensor DHX33 suppressed GPS167-mediated cytotoxicity on HCT116 cells. Our study uncovered molecular mechanisms through which 1C8 and GPS167 affect cancer cell proliferation as well as processes critical for metastasis.

Regulatory network structure as a dominant determinant of transcription factor evolutionary rate.

The evolution of transcriptional regulatory networks has thus far mostly been studied at the level of cis-regulatory elements. To gain a complete understanding of regulatory network evolution we must also study the evolutionary role of trans-factors, such as transcription factors (TFs). Here, we systematically assess genomic and network-level determinants of TF evolutionary rate in yeast, and how they compare to those of generic proteins, while carefully controlling for differences of the TF protein set, such as expression level. We found significantly distinct trends relating TF evolutionary rate to mRNA expression level, codon adaptation index, the evolutionary rate of physical interaction partners, and, confirming previous reports, to protein-protein interaction degree and regulatory in-degree. We discovered that for TFs, the dominant determinants of evolutionary rate lie in the structure of the regulatory network, such as the median evolutionary rate of target genes and the fraction of species-specific target genes. Decomposing the regulatory network by edge sign, we found that this modular evolution of TFs and their targets is limited to activating regulatory relationships. We show that fast evolving TFs tend to regulate other TFs and niche-specific processes and that their targets show larger evolutionary expression changes than targets of other TFs. We also show that the positive trend relating TF regulatory in-degree and evolutionary rate is likely related to the species-specificity of the transcriptional regulation modules. Finally, we discuss likely causes for TFs' different evolutionary relationship to the physical interaction network, such as the prevalence of transient interactions in the TF subnetwork. This work suggests that positive and negative regulatory networks follow very different evolutionary rules, and that transcription factor evolution is best understood at a network- or systems-level.

Fine-scale variation and genetic determinants of alternative splicing across individuals.

Recently, thanks to the increasing throughput of new technologies, we have begun to explore the full extent of alternative pre-mRNA splicing (AS) in the human transcriptome. This is unveiling a vast layer of complexity in isoform-level expression differences between individuals. We used previously published splicing sensitive microarray data from lymphoblastoid cell lines to conduct an in-depth analysis on splicing efficiency of known and predicted exons. By combining publicly available AS annotation with a novel algorithm designed to search for AS, we show that many real AS events can be detected within the usually unexploited, speculative majority of the array and at significance levels much below standard multiple-testing thresholds, demonstrating that the extent of cis-regulated differential splicing between individuals is potentially far greater than previously reported. Specifically, many genes show subtle but significant genetically controlled differences in splice-site usage. PCR validation shows that 42 out of 58 (72%) candidate gene regions undergo detectable AS, amounting to the largest scale validation of isoform eQTLs to date. Targeted sequencing revealed a likely causative SNP in most validated cases. In all 17 incidences where a SNP affected a splice-site region, in silico splice-site strength modeling correctly predicted the direction of the micro-array and PCR results. In 13 other cases, we identified likely causative SNPs disrupting predicted splicing enhancers. Using Fst and REHH analysis, we uncovered significant evidence that 2 putative causative SNPs have undergone recent positive selection. We verified the effect of five SNPs using in vivo minigene assays. This study shows that splicing differences between individuals, including quantitative differences in isoform ratios, are frequent in human populations and that causative SNPs can be identified using in silico predictions. Several cases affected disease-relevant genes and it is likely some of these differences are involved in phenotypic diversity and susceptibility to complex diseases.

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.

Vesicular trafficking is a key determinant of the statin response in acute myeloid leukemia.

Cholesterol homeostasis has been proposed as one mechanism contributing to chemoresistance in AML and hence, inclusion of statins in therapeutic regimens as part of clinical trials in AML has shown encouraging results. Chemical screening of primary human AML specimens by our group led to the identification of lipophilic statins as potent inhibitors of AMLs from a wide range of cytogenetic groups. Genetic screening to identify modulators of the statin response uncovered the role of protein geranylgeranylation and of RAB proteins, coordinating various aspect of vesicular trafficking, in mediating the effects of statins on AML cell viability. We further show that statins can inhibit vesicle-mediated transport in primary human specimens, and that statins sensitive samples show expression signatures reminiscent of enhanced vesicular trafficking. Overall, this study sheds light into the mechanism of action of statins in AML and identifies a novel vulnerability for cytogenetically diverse AML.

A novel p53 regulator, C16ORF72/TAPR1, buffers against telomerase inhibition.

Telomere erosion in cells with insufficient levels of the telomerase reverse transcriptase (TERT), contributes to age-associated tissue dysfunction and senescence, and p53 plays a crucial role in this response. We undertook a genome-wide CRISPR screen to identify gene deletions that sensitized p53-positive human cells to telomerase inhibition. We uncovered a previously unannotated gene, C16ORF72, which we term Telomere Attrition and p53 Response 1 (TAPR1), that exhibited a synthetic-sick relationship with TERT loss. A subsequent genome-wide CRISPR screen in TAPR1-disrupted cells reciprocally identified TERT as a sensitizing gene deletion. Cells lacking TAPR1 or TERT possessed elevated p53 levels and transcriptional signatures consistent with p53 upregulation. The elevated p53 response in TERT- or TAPR1-deficient cells was exacerbated by treatment with the MDM2 inhibitor and p53 stabilizer nutlin-3a and coincided with a further reduction in cell fitness. Importantly, the sensitivity to treatment with nutlin-3a in TERT- or TAPR1-deficient cells was rescued by loss of p53. These data suggest that TAPR1 buffers against the deleterious consequences of telomere erosion or DNA damage by constraining p53. These findings identify C16ORF72/TAPR1 as new regulator at the nexus of telomere integrity and p53 regulation.

UM171 Preserves Epigenetic Marks that Are Reduced in Ex Vivo Culture of Human HSCs via Potentiation of the CLR3-KBTBD4 Complex.

Human hematopoietic stem cells (HSCs) exhibit attrition of their self-renewal capacity when cultured ex vivo, a process that is partially reversed upon treatment with epigenetic modifiers, most notably inhibitors of histone deacetylases (HDACs) or lysine-specific demethylase LSD1. A recent study showed that the human HSC self-renewal agonist UM171 modulates the CoREST complex, leading to LSD1 degradation, whose inhibition mimics the activity of UM171. The mechanism underlying the UM171-mediated loss of CoREST function remains undetermined. We now report that UM171 potentiates the activity of a CULLIN3-E3 ubiquitin ligase (CRL3) complex whose target specificity is dictated by the poorly characterized Kelch/BTB domain protein KBTBD4. CRL3KBTBD4 targets components of the LSD1/RCOR1 corepressor complex for proteasomal degradation, hence re-establishing H3K4me2 and H3K27ac epigenetic marks, which are rapidly decreased upon ex vivo culture of human HSCs.

Identification and optimization of molecular glue compounds that inhibit a noncovalent E2 enzyme-ubiquitin complex.

Pharmacological control of the ubiquitin-proteasome system (UPS) is of intense interest in drug discovery. Here, we report the development of chemical inhibitors of the ubiquitin-conjugating (E2) enzyme CDC34A (also known as UBE2R1), which donates activated ubiquitin to the cullin-RING ligase (CRL) family of ubiquitin ligase (E3) enzymes. A FRET-based interaction assay was used to screen for novel compounds that stabilize the noncovalent complex between CDC34A and ubiquitin, and thereby inhibit the CDC34A catalytic cycle. An isonipecotamide hit compound was elaborated into analogs with ~1000-fold increased potency in stabilizing the CDC34A-ubiquitin complex. These analogs specifically inhibited CDC34A-dependent ubiquitination in vitro and stabilized an E2~ubiquitin thioester reaction intermediate in cells. The x-ray crystal structure of a CDC34A-ubiquitin-inhibitor complex uncovered the basis for analog structure-activity relationships. The development of chemical stabilizers of the CDC34A-ubiquitin complex illustrates a general strategy for de novo discovery of molecular glue compounds that stabilize weak protein interactions.

A novel class of inhibitors that target SRSF10 and promote p53-mediated cytotoxicity on human colorectal cancer cells.

The elevated expression of the splicing regulator SRSF10 in metastatic colorectal cancer (CRC) stimulates the production of the pro-tumorigenic BCLAF1-L splice variant. We discovered a group of small molecules with an aminothiazole carboxamide core (GPS167, GPS192 and others) that decrease production of BCLAF1-L. While additional alternative splicing events regulated by SRSF10 are affected by GPS167/192 in HCT116 cells (e.g. in MDM4, WTAP, SLK1 and CLK1), other events are shifted in a SRSF10-independent manner (e.g. in MDM2, NAB2 and TRA2A). GPS167/192 increased the interaction of SRSF10 with the CLK1 and CLK4 kinases, leading us to show that GPS167/192 can inhibit CLK kinases preferentially impacting the activity of SRSF10. Notably, GPS167 impairs the growth of CRC cell lines and organoids, inhibits anchorage-independent colony formation, cell migration, and promotes cytoxicity in a manner that requires SRSF10 and p53. In contrast, GPS167 only minimally affects normal colonocytes and normal colorectal organoids. Thus, GPS167 reprograms the tumorigenic activity of SRSF10 in CRC cells to elicit p53-dependent apoptosis.