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.

Germline-encoded amino acid-binding motifs drive immunodominant public antibody responses.

Despite the vast diversity of the antibody repertoire, infected individuals often mount antibody responses to precisely the same epitopes within antigens. The immunological mechanisms underpinning this phenomenon remain unknown. By mapping 376 immunodominant "public epitopes" at high resolution and characterizing several of their cognate antibodies, we concluded that germline-encoded sequences in antibodies drive recurrent recognition. Systematic analysis of antibody-antigen structures uncovered 18 human and 21 partially overlapping mouse germline-encoded amino acid-binding (GRAB) motifs within heavy and light V gene segments that in case studies proved critical for public epitope recognition. GRAB motifs represent a fundamental component of the immune system's architecture that promotes recognition of pathogens and leads to species-specific public antibody responses that can exert selective pressure on pathogens.

Gain-of-function genetic screening identifies the antiviral function of TMEM120A via STING activation.

Zika virus (ZIKV) infection can be associated with neurological pathologies, such as microcephaly in newborns and Guillain-Barre syndrome in adults. Effective therapeutics are currently not available. As such, a comprehensive understanding of virus-host interactions may guide the development of medications for ZIKV. Here we report a human genome-wide overexpression screen to identify host factors that regulate ZIKV infection and find TMEM120A as a ZIKV restriction factor. TMEM120A overexpression significantly inhibits ZIKV replication, while TMEM120A knockdown increases ZIKV infection in cell lines. Moreover, Tmem120a knockout in mice facilitates ZIKV infection in primary mouse embryonic fibroblasts (MEF) cells. Mechanistically, the antiviral activity of TMEM120A is dependent on STING, as TMEM120A interacts with STING, promotes the translocation of STING from the endoplasmic reticulum (ER) to ER-Golgi intermediate compartment (ERGIC) and enhances the phosphorylation of downstream TBK1 and IRF3, resulting in the expression of multiple antiviral cytokines and interferon-stimulated genes. In summary, our gain-of-function screening identifies TMEM120A as a key activator of the antiviral signaling of STING.

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.

Integrated proteogenetic analysis reveals the landscape of a mitochondrial-autophagosome synapse during PARK2-dependent mitophagy.

The PINK1 protein kinase activates the PARK2 ubiquitin ligase to promote mitochondrial ubiquitylation and recruitment of ubiquitin-binding mitophagy receptors typified by OPTN and TAX1BP1. Here, we combine proximity biotinylation of OPTN and TAX1BP1 with CRISPR-Cas9-based screens for mitophagic flux to develop a spatial proteogenetic map of PARK2-dependent mitophagy. Proximity labeling of OPTN allowed visualization of a "mitochondrial-autophagosome synapse" upon mitochondrial depolarization. Proximity proteomics of OPTN and TAX1BP1 revealed numerous proteins at the synapse, including both PARK2 substrates and autophagy components. Parallel mitophagic flux screens identified proteins with roles in autophagy, vesicle formation and fusion, as well as PARK2 targets, many of which were also identified via proximity proteomics. One protein identified in both approaches, HK2, promotes assembly of a high-molecular weight complex of PINK1 and phosphorylation of ubiquitin in response to mitochondrial damage. This work provides a resource for understanding the spatial and molecular landscape of PARK2-dependent mitophagy.

Measles virus infection diminishes preexisting antibodies that offer protection from other pathogens.

Measles virus is directly responsible for more than 100,000 deaths yearly. Epidemiological studies have associated measles with increased morbidity and mortality for years after infection, but the reasons why are poorly understood. Measles virus infects immune cells, causing acute immune suppression. To identify and quantify long-term effects of measles on the immune system, we used VirScan, an assay that tracks antibodies to thousands of pathogen epitopes in blood. We studied 77 unvaccinated children before and 2 months after natural measles virus infection. Measles caused elimination of 11 to 73% of the antibody repertoire across individuals. Recovery of antibodies was detected after natural reexposure to pathogens. Notably, these immune system effects were not observed in infants vaccinated against MMR (measles, mumps, and rubella), but were confirmed in measles-infected macaques. The reduction in humoral immune memory after measles infection generates potential vulnerability to future infections, underscoring the need for widespread vaccination.

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.