Defining E3 ligase-substrate relationships through multiplex CRISPR screening.

Specificity within the ubiquitin-proteasome system is primarily achieved through E3 ubiquitin ligases, but for many E3s their substrates-and in particular the molecular features (degrons) that they recognize-remain largely unknown. Current approaches for assigning E3s to their cognate substrates are tedious and low throughput. Here we developed a multiplex CRISPR screening platform to assign E3 ligases to their cognate substrates at scale. A proof-of-principle multiplex screen successfully performed ~100 CRISPR screens in a single experiment, refining known C-degron pathways and identifying an additional pathway through which Cul2FEM1B targets C-terminal proline. Further, by identifying substrates for Cul1FBXO38, Cul2APPBP2, Cul3GAN, Cul3KLHL8, Cul3KLHL9/13 and Cul3KLHL15, we demonstrate that the approach is compatible with pools of full-length protein substrates of varying stabilities and, when combined with site-saturation mutagenesis, can assign E3 ligases to their cognate degron motifs. Thus, multiplex CRISPR screening will accelerate our understanding of how specificity is achieved within the ubiquitin-proteasome system.

High-resolution epitope mapping by AllerScan reveals relationships between IgE and IgG repertoires during peanut oral immunotherapy.

Peanut allergy can result in life-threatening reactions and is a major public health concern. Oral immunotherapy (OIT) induces desensitization to food allergens through administration of increasing amounts of allergen. To dissect peanut-specific immunoglobulin E (IgE) and IgG responses in subjects undergoing OIT, we have developed AllerScan, a method that leverages phage-display and next-generation sequencing to identify the epitope targets of peanut-specific antibodies. We observe a striking diversification and boosting of the peanut-specific IgG repertoire after OIT and a reduction in pre-existing IgE levels against individual epitopes. High-resolution epitope mapping reveals shared recognition of public epitopes in Ara h 1, 2, 3, and 7. In individual subjects, OIT-induced IgG specificities overlap extensively with IgE and exhibit strikingly similar antibody footprints, suggesting related clonal lineages or convergent evolution of peanut-specific IgE and IgG B cells. Individual differences in epitope recognition identified via AllerScan could inform safer and more effective personalized immunotherapy.

The adaptive immune system is a major driver of selection for tumor suppressor gene inactivation.

During tumorigenesis, tumors must evolve to evade the immune system and do so by disrupting the genes involved in antigen processing and presentation or up-regulating inhibitory immune checkpoint genes. We performed in vivo CRISPR screens in syngeneic mouse tumor models to examine requirements for tumorigenesis both with and without adaptive immune selective pressure. In each tumor type tested, we found a marked enrichment for the loss of tumor suppressor genes (TSGs) in the presence of an adaptive immune system relative to immunocompromised mice. Nearly one-third of TSGs showed preferential enrichment, often in a cancer- and tissue-specific manner. These results suggest that clonal selection of recurrent mutations found in cancer is driven largely by the tumor's requirement to avoid the adaptive immune system.

CARM1 Inhibition Enables Immunotherapy of Resistant Tumors by Dual Action on Tumor Cells and T Cells.

A number of cancer drugs activate innate immune pathways in tumor cells but unfortunately also compromise antitumor immune function. We discovered that inhibition of CARM1, an epigenetic enzyme and cotranscriptional activator, elicited beneficial antitumor activity in both cytotoxic T cells and tumor cells. In T cells, Carm1 inactivation substantially enhanced their antitumor function and preserved memory-like populations required for sustained antitumor immunity. In tumor cells, Carm1 inactivation induced a potent type 1 interferon response that sensitized resistant tumors to cytotoxic T cells. Substantially increased numbers of dendritic cells, CD8 T cells, and natural killer cells were present in Carm1-deficient tumors, and infiltrating CD8 T cells expressed low levels of exhaustion markers. Targeting of CARM1 with a small molecule elicited potent antitumor immunity and sensitized resistant tumors to checkpoint blockade. Targeting of this cotranscriptional regulator thus offers an opportunity to enhance immune function while simultaneously sensitizing resistant tumor cells to immune attack. SIGNIFICANCE: Resistance to cancer immunotherapy remains a major challenge. Targeting of CARM1 enables immunotherapy of resistant tumors by enhancing T-cell functionality and preserving memory-like T-cell populations within tumors. CARM1 inhibition also sensitizes resistant tumor cells to immune attack by inducing a tumor cell-intrinsic type 1 interferon response.This article is highlighted in the In This Issue feature, p. 1861.

Antibody responses to endemic coronaviruses modulate COVID-19 convalescent plasma functionality.

SARS-CoV-2 (CoV2) antibody therapies, including COVID-19 convalescent plasma (CCP), monoclonal antibodies, and hyperimmune globulin, are among the leading treatments for individuals with early COVID-19 infection. The functionality of convalescent plasma varies greatly, but the association of antibody epitope specificities with plasma functionality remains uncharacterized. We assessed antibody functionality and reactivities to peptides across the CoV2 and the 4 endemic human coronavirus (HCoV) genomes in 126 CCP donations. We found strong correlation between plasma functionality and polyclonal antibody targeting of CoV2 spike protein peptides. Antibody reactivity to many HCoV spike peptides also displayed strong correlation with plasma functionality, including pan-coronavirus cross-reactive epitopes located in a conserved region of the fusion peptide. After accounting for antibody cross-reactivity, we identified an association between greater alphacoronavirus NL63 antibody responses and development of highly neutralizing antibodies against CoV2. We also found that plasma preferentially reactive to the CoV2 spike receptor binding domain (RBD), versus the betacoronavirus HKU1 RBD, had higher neutralizing titer. Finally, we developed a 2-peptide serosignature that identifies plasma donations with high anti-spike titer, but that suffer from low neutralizing activity. These results suggest that analysis of coronavirus antibody fine specificities may be useful for selecting desired therapeutics and understanding the complex immune responses elicited by CoV2 infection.

Antibody responses to endemic coronaviruses modulate COVID-19 convalescent plasma functionality.

COVID-19 convalescent plasma, particularly plasma with high-titer SARS-CoV-2 (CoV2) antibodies, has been successfully used for treatment of COVID-19. The functionality of convalescent plasma varies greatly, but the association of antibody epitope specificities with plasma functionality remains uncharacterized. We assessed antibody functionality and reactivities to peptides across the CoV2 and the four endemic human coronavirus (HCoV) genomes in 126 COVID-19 convalescent plasma donations. We found strong correlation between plasma functionality and polyclonal antibody targeting of CoV2 spike protein peptides. Antibody reactivity to many HCoV spike peptides also displayed strong correlation with plasma functionality, including pan-coronavirus cross-reactive epitopes located in a conserved region of the fusion peptide. After accounting for antibody cross-reactivity, we identified an association between greater alphacoronavirus NL63 antibody responses and development of highly neutralizing antibodies to SARS-CoV-2. We also found that plasma preferentially reactive to the CoV2 receptor binding domain (RBD), versus the betacoronavirus HKU1 RBD, had higher neutralizing titer. Finally, we developed a two-peptide serosignature that identifies plasma donations with high anti-S titer but that suffer from low neutralizing activity. These results suggest that analysis of coronavirus antibody fine specificities may be useful for selecting therapeutic plasma with desired functionalities.

Viral epitope profiling of COVID-19 patients reveals cross-reactivity and correlates of severity.

Understanding humoral responses to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is critical for improving diagnostics, therapeutics, and vaccines. Deep serological profiling of 232 coronavirus disease 2019 (COVID-19) patients and 190 pre-COVID-19 era controls using VirScan revealed more than 800 epitopes in the SARS-CoV-2 proteome, including 10 epitopes likely recognized by neutralizing antibodies. Preexisting antibodies in controls recognized SARS-CoV-2 ORF1, whereas only COVID-19 patient antibodies primarily recognized spike protein and nucleoprotein. A machine learning model trained on VirScan data predicted SARS-CoV-2 exposure history with 99% sensitivity and 98% specificity; a rapid Luminex-based diagnostic was developed from the most discriminatory SARS-CoV-2 peptides. Individuals with more severe COVID-19 exhibited stronger and broader SARS-CoV-2 responses, weaker antibody responses to prior infections, and higher incidence of cytomegalovirus and herpes simplex virus 1, possibly influenced by demographic covariates. Among hospitalized patients, males produce stronger SARS-CoV-2 antibody responses than females.

Diversified Application of Barcoded PLATO (PLATO-BC) Platform for Identification of Protein Interactions.

Proteins usually associate with other molecules physically to execute their functions. Identifying these interactions is important for the functional analysis of proteins. Previously, we reported the parallel analysis of translated ORFs (PLATO) to couple ribosome display of full-length ORFs with affinity enrichment of mRNA/protein/ribosome complexes for the "bait" molecules, followed by the deep sequencing analysis of mRNA. However, the sample processing, from extraction of precipitated mRNA to generation of DNA libraries, includes numerous steps, which is tedious and may cause the loss of materials. Barcoded PLATO (PLATO-BC), an improved platform was further developed to test its application for protein interaction discovery. In this report, we tested the antisera-antigen interaction using serum samples from patients with inclusion body myositis (IBM). Tripartite motif containing 21 (TRIM21) was identified as a potentially new IBM autoantigen. We also expanded the application of PLATO-BC to identify protein interactions for JQ1, single ubiquitin peptide, and NS5 protein of Zika virus. From PLATO-BC analyses, we identified new protein interactions for these "bait" molecules. We demonstrate that Ewing sarcoma breakpoint region 1 (EWSR1) binds to JQ1 and their interactions may interrupt the EWSR1 association with acetylated histone H4. RIO kinase 3 (RIOK3), a newly identified ubiquitin-binding protein, is preferentially associated with K63-ubiquitin chain. We also find that Zika NS5 protein interacts with two previously unreported host proteins, par-3 family cell polarity regulator (PARD3) and chromosome 19 open reading frame 53 (C19orf53), whose attenuated expression benefits the replication of Zika virus. These results further demonstrate that PLATO-BC is capable of identifying novel protein interactions for various types of "bait" molecules.

A Druggable Genome Screen Identifies Modifiers of α-Synuclein Levels via a Tiered Cross-Species Validation Approach.

Accumulation of α-Synuclein (α-Syn) causes Parkinson's disease (PD) as well as other synucleopathies. α-Syn is the major component of Lewy bodies and Lewy neurites, the proteinaceous aggregates that are a hallmark of sporadic PD. In familial forms of PD, mutations or copy number variations in SNCA (the α-Syn gene) result in a net increase of its protein levels. Furthermore, common risk variants tied to PD are associated with small increases of wild-type α-Syn levels. These findings are further bolstered by animal studies which show that overexpression of α-Syn is sufficient to cause PD-like features. Thus, increased α-Syn levels are intrinsically tied to PD pathogenesis and underscore the importance of identifying the factors that regulate its levels. In this study, we establish a pooled RNAi screening approach and validation pipeline to probe the druggable genome for modifiers of α-Syn levels and identify 60 promising targets. Using a cross-species, tiered validation approach, we validate six strong candidates that modulate α-Syn levels and toxicity in cell lines, Drosophila, human neurons, and mouse brain of both sexes. More broadly, this genetic strategy and validation pipeline can be applied for the identification of therapeutic targets for disorders driven by dosage-sensitive proteins.SIGNIFICANCE STATEMENT We present a research strategy for the systematic identification and validation of genes modulating the levels of α-Synuclein, a protein involved in Parkinson's disease. A cell-based screen of the druggable genome (>7,500 genes that are potential therapeutic targets) yielded many modulators of α-Synuclein that were subsequently confirmed and validated in Drosophila, human neurons, and mouse brain. This approach has broad applicability to the multitude of neurological diseases that are caused by mutations in genes whose dosage is critical for brain function.

The Eukaryotic Proteome Is Shaped by E3 Ubiquitin Ligases Targeting C-Terminal Degrons.

Degrons are minimal elements that mediate the interaction of proteins with degradation machineries to promote proteolysis. Despite their central role in proteostasis, the number of known degrons remains small, and a facile technology to characterize them is lacking. Using a strategy combining global protein stability (GPS) profiling with a synthetic human peptidome, we identify thousands of peptides containing degron activity. Employing CRISPR screening, we establish that the stability of many proteins is regulated through degrons located at their C terminus. We characterize eight Cullin-RING E3 ubiquitin ligase (CRL) complex adaptors that regulate C-terminal degrons, including six CRL2 and two CRL4 complexes, and computationally implicate multiple non-CRLs in end recognition. Proteome analysis revealed that the C termini of eukaryotic proteins are depleted for C-terminal degrons, suggesting an E3-ligase-dependent modulation of proteome composition. Thus, we propose that a series of "C-end rules" operate to govern protein stability and shape the eukaryotic proteome.

Profound Tissue Specificity in Proliferation Control Underlies Cancer Drivers and Aneuploidy Patterns.

Genomics has provided a detailed structural description of the cancer genome. Identifying oncogenic drivers that work primarily through dosage changes is a current challenge. Unrestrained proliferation is a critical hallmark of cancer. We constructed modular, barcoded libraries of human open reading frames (ORFs) and performed screens for proliferation regulators in multiple cell types. Approximately 10% of genes regulate proliferation, with most performing in an unexpectedly highly tissue-specific manner. Proliferation drivers in a given cell type showed specific enrichment in somatic copy number changes (SCNAs) from cognate tumors and helped predict aneuploidy patterns in those tumors, implying that tissue-type-specific genetic network architectures underlie SCNA and driver selection in different cancers. In vivo screening confirmed these results. We report a substantial contribution to the catalog of SCNA-associated cancer drivers, identifying 147 amplified and 107 deleted genes as potential drivers, and derive insights about the genetic network architecture of aneuploidy in tumors.

A Role for Mitochondrial Translation in Promotion of Viability in K-Ras Mutant Cells.

Activating mutations in the KRAS oncogene are highly prevalent in tumors, especially those of the colon, lung, and pancreas. To better understand the genetic dependencies that K-Ras mutant cells rely upon for their growth, we employed whole-genome CRISPR loss-of-function screens in two isogenic pairs of cell lines. Since loss of essential genes is uniformly toxic in CRISPR-based screens, we also developed a small hairpin RNA (shRNA) library targeting essential genes. These approaches uncovered a large set of proteins whose loss results in the selective reduction of K-Ras mutant cell growth. Pathway analysis revealed that many of these genes function in the mitochondria. For validation, we generated isogenic pairs of cell lines using CRISPR-based genome engineering, which confirmed the dependency of K-Ras mutant cells on these mitochondrial pathways. Finally, we found that mitochondrial inhibitors reduce the growth of K-Ras mutant tumors in vivo, aiding in the advancement of strategies to target K-Ras-driven malignancy.

A genetic interaction analysis identifies cancer drivers that modify EGFR dependency.

A large number of cancer drivers have been identified through tumor sequencing efforts, but how they interact and the degree to which they can substitute for each other have not been systematically explored. To comprehensively investigate how cancer drivers genetically interact, we searched for modifiers of epidermal growth factor receptor (EGFR) dependency by performing CRISPR, shRNA, and expression screens in a non-small cell lung cancer (NSCLC) model. We elucidated a broad spectrum of tumor suppressor genes (TSGs) and oncogenes (OGs) that can genetically modify proliferation and survival of cancer cells when EGFR signaling is altered. These include genes already known to mediate EGFR inhibitor resistance as well as many TSGs not previously connected to EGFR and whose biological functions in tumorigenesis are not well understood. We show that mutation of PBRM1, a subunit of the SWI/SNF complex, attenuates the effects of EGFR inhibition in part by sustaining AKT signaling. We also show that mutation of Capicua (CIC), a transcriptional repressor, suppresses the effects of EGFR inhibition by partially restoring the EGFR-promoted gene expression program, including the sustained expression of Ets transcription factors such as ETV1 Together, our data provide strong support for the hypothesis that many cancer drivers can substitute for each other in certain contexts and broaden our understanding of EGFR regulation.

Systematic autoantigen analysis identifies a distinct subtype of scleroderma with coincident cancer.

Scleroderma is a chronic autoimmune rheumatic disease associated with widespread tissue fibrosis and vasculopathy. Approximately two-thirds of all patients with scleroderma present with three dominant autoantibody subsets. Here, we used a pair of complementary high-throughput methods for antibody epitope discovery to examine patients with scleroderma with or without known autoantibody specificities. We identified a specificity for the minor spliceosome complex containing RNA Binding Region (RNP1, RNA recognition motif) Containing 3 (RNPC3) that is found in patients with scleroderma without known specificities and is absent in unrelated autoimmune diseases. We found strong evidence for both intra- and intermolecular epitope spreading in patients with RNA polymerase III (POLR3) and the minor spliceosome specificities. Our results demonstrate the utility of these technologies in rapidly identifying antibodies that can serve as biomarkers of disease subsets in the evolving precision medicine era.

Sources of Error in Mammalian Genetic Screens.

Genetic screens are invaluable tools for dissection of biological phenomena. Optimization of such screens to enhance discovery of candidate genes and minimize false positives is thus a critical aim. Here, we report several sources of error common to pooled genetic screening techniques used in mammalian cell culture systems, and demonstrate methods to eliminate these errors. We find that reverse transcriptase-mediated recombination during retroviral replication can lead to uncoupling of molecular tags, such as DNA barcodes (BCs), from their associated library elements, leading to chimeric proviral genomes in which BCs are paired to incorrect ORFs, shRNAs, etc This effect depends on the length of homologous sequence between unique elements, and can be minimized with careful vector design. Furthermore, we report that residual plasmid DNA from viral packaging procedures can contaminate transduced cells. These plasmids serve as additional copies of the PCR template during library amplification, resulting in substantial inaccuracies in measurement of initial reference populations for screen normalization. The overabundance of template in some samples causes an imbalance between PCR cycles of contaminated and uncontaminated samples, which results in a systematic artifactual depletion of GC-rich library elements. Elimination of contaminating plasmid DNA using the bacterial endonuclease Benzonase can restore faithful measurements of template abundance and minimize GC bias.

Functional genomics reveals that tumors with activating phosphoinositide 3-kinase mutations are dependent on accelerated protein turnover.

Activating mutations in the phosphoinositide 3-kinase (PI3K) signaling pathway are frequently identified in cancer. To identify pathways that support PI3K oncogenesis, we performed a genome-wide RNAi screen in isogenic cell lines harboring wild-type or mutant PIK3CA to search for PI3K synthetic-lethal (SL) genes. A combined analysis of these results with a meta-analysis of two other large-scale RNAi screening data sets in PI3K mutant cancer cell lines converged on ribosomal protein translation and proteasomal protein degradation as critical nononcogene dependencies for PI3K-driven tumors. Genetic or pharmacologic inhibition of either pathway alone, but not together, selectively killed PI3K mutant tumor cells in an mTOR-dependent manner. The expression of ribosomal and proteasomal components was significantly up-regulated in primary human colorectal tumors harboring PI3K pathway activation. Importantly, a PI3K SL gene signature containing the top hits of the SL genes identified in our meta-analysis robustly predicted overall patient survival in colorectal cancer, especially among patients with tumors with an activated PI3K pathway. These results suggest that disruption of protein turnover homeostasis via ribosome or proteasome inhibition may be a novel treatment strategy for PI3K mutant human tumors.

The DNA damage response induces inflammation and senescence by inhibiting autophagy of GATA4.

Cellular senescence is a terminal stress-activated program controlled by the p53 and p16(INK4a) tumor suppressor proteins. A striking feature of senescence is the senescence-associated secretory phenotype (SASP), a pro-inflammatory response linked to tumor promotion and aging. We have identified the transcription factor GATA4 as a senescence and SASP regulator. GATA4 is stabilized in cells undergoing senescence and is required for the SASP. Normally, GATA4 is degraded by p62-mediated selective autophagy, but this regulation is suppressed during senescence, thereby stabilizing GATA4. GATA4 in turn activates the transcription factor NF-κB to initiate the SASP and facilitate senescence. GATA4 activation depends on the DNA damage response regulators ATM and ATR, but not on p53 or p16(INK4a). GATA4 accumulates in multiple tissues, including the aging brain, and could contribute to aging and its associated inflammation.

Viral immunology. Comprehensive serological profiling of human populations using a synthetic human virome.

The human virome plays important roles in health and immunity. However, current methods for detecting viral infections and antiviral responses have limited throughput and coverage. Here, we present VirScan, a high-throughput method to comprehensively analyze antiviral antibodies using immunoprecipitation and massively parallel DNA sequencing of a bacteriophage library displaying proteome-wide peptides from all human viruses. We assayed over 10(8) antibody-peptide interactions in 569 humans across four continents, nearly doubling the number of previously established viral epitopes. We detected antibodies to an average of 10 viral species per person and 84 species in at least two individuals. Although rates of specific virus exposure were heterogeneous across populations, antibody responses targeted strongly conserved "public epitopes" for each virus, suggesting that they may elicit highly similar antibodies. VirScan is a powerful approach for studying interactions between the virome and the immune system.

STOP gene Phactr4 is a tumor suppressor.

Cancer develops through genetic and epigenetic alterations that allow unrestrained proliferation and increased survival. Using a genetic RNAi screen, we previously identified hundreds of suppressors of tumorigenesis and/or proliferation (STOP) genes that restrain normal cell proliferation. Our STOP gene set was significantly enriched for known and putative tumor suppressor genes. Here, we report a tumor-suppressive role for one STOP gene, phosphatase and actin regulator 4 (PHACTR4). Phactr4 is one of four members of the largely uncharacterized Phactr family of protein phosphatase 1 (PP1)-and actin-binding proteins. Our work suggests that Phactr4 restrains normal cell proliferation and transformation. Depletion of Phactr4 with multiple shRNAs leads to increased proliferation and soft agar colony formation. Phactr4 acts, in part, through an Rb-dependent pathway, because Rb phosphorylation is maintained upon growth factor withdrawal in Phactr4-depleted cells. Examination of tumor copy number analysis and sequencing revealed that PHACTR4 is significantly deleted and mutant in many tumor subtypes. Furthermore,cancer cell lines with reduced Phactr4 expression exhibit tumor suppressor hypersensitivity upon Phactr4 complementation,leading to reduced proliferation, transformation, and tumor formation. Thus, Phactr4 acts as a tumor suppressor that is deleted and mutant in several cancers.