Identification of orphan ligand-receptor relationships using a cell-based CRISPRa enrichment screening platform.

Secreted proteins, which include cytokines, hormones, and growth factors, are extracellular ligands that control key signaling pathways mediating cell-cell communication within and between tissues and organs. Many drugs target secreted ligands and their cell surface receptors. Still, there are hundreds of secreted human proteins that either have no identified receptors ('orphans') or are likely to act through cell surface receptors that have not yet been characterized. Discovery of secreted ligand-receptor interactions by high-throughput screening has been problematic, because the most commonly used high-throughput methods for protein-protein interaction (PPI) screening are not optimized for extracellular interactions. Cell-based screening is a promising technology for the deorphanization of ligand-receptor interactions, because multimerized ligands can enrich for cells expressing low affinity cell surface receptors, and such methods do not require purification of receptor extracellular domains. Here, we present a proteo-genomic cell-based CRISPR activation (CRISPRa) enrichment screening platform employing customized pooled cell surface receptor sgRNA libraries in combination with a magnetic bead selection-based enrichment workflow for rapid, parallel ligand-receptor deorphanization. We curated 80 potentially high-value orphan secreted proteins and ultimately screened 20 secreted ligands against two cell sgRNA libraries with targeted expression of all single-pass (TM1) or multi-pass transmembrane (TM2+) receptors by CRISPRa. We identified previously unknown interactions in 12 of these screens, and validated several of them using surface plasmon resonance and/or cell binding assays. The newly deorphanized ligands include three receptor protein tyrosine phosphatase (RPTP) ligands and a chemokine-like protein that binds to killer immunoglobulin-like receptors (KIRs). These new interactions provide a resource for future investigations of interactions between the human-secreted and membrane proteomes.

Phage Particles of Controlled Length and Genome for In Vivo Targeted Glioblastoma Imaging and Therapeutic Delivery.

M13 bacteriophage (phage) are versatile, genetically tunable nanocarriers that have been recently adapted for use as diagnostic and therapeutic platforms. Applying p3 capsid chlorotoxin fusion with the "inho" circular single-stranded DNA (cssDNA) gene packaging system, we produced miniature chlorotoxin inho (CTX-inho) phage particles with a minimum length of 50 nm that can target intracranial orthotopic patient-derived GBM22 glioblastoma tumors in the brains of mice. Systemically administered indocyanine green conjugated CTX-inho phage accumulated in brain tumors, facilitating shortwave infrared detection. Furthermore, we show that our inho phage can carry cssDNA that are transcriptionally active when delivered to GBM22 glioma cells in vitro. The ability to modulate the capsid display, surface loading, phage length, and cssDNA gene content makes the recombinant M13 phage particle an ideal delivery platform.

Roadmap for the use of base editors to decipher drug mechanism of action.

CRISPR base editors are powerful tools for large-scale mutagenesis studies. This kind of approach can elucidate the mechanism of action of compounds, a key process in drug discovery. Here, we explore the utility of base editors in an early drug discovery context focusing on G-protein coupled receptors. A pooled mutagenesis screening framework was set up based on a modified version of the CRISPR-X base editor system. We determine optimized experimental conditions for mutagenesis where sgRNAs are delivered by cell transfection or viral infection over extended time periods (>14 days), resulting in high mutagenesis produced in a short region located at -4/+8 nucleotides with respect to the sgRNA match. The ß2 Adrenergic Receptor (B2AR) was targeted in this way employing a 6xCRE-mCherry reporter system to monitor its response to isoproterenol. The results of our screening indicate that residue 184 of B2AR is crucial for its activation. Based on our experience, we outline the crucial points to consider when designing and performing CRISPR-based pooled mutagenesis screening, including the typical technical hurdles encountered when studying compound pharmacology.

Genome-wide CRISPR screens reveal a specific ligand for the glycan-binding immune checkpoint receptor Siglec-7.

Glyco-immune checkpoint receptors, molecules that inhibit immune cell activity following binding to glycosylated cell-surface antigens, are emerging as attractive targets for cancer immunotherapy. Defining biologically relevant ligands that bind and activate such receptors, however, has historically been a significant challenge. Here, we present a CRISPRi genomic screening strategy that allowed unbiased identification of the key genes required for cell-surface presentation of glycan ligands on leukemia cells that bind the glyco-immune checkpoint receptors Siglec-7 and Siglec-9. This approach revealed a selective interaction between Siglec-7 and the mucin-type glycoprotein CD43. Further work identified a specific N-terminal glycopeptide region of CD43 containing clusters of disialylated O-glycan tetrasaccharides that form specific Siglec-7 binding motifs. Knockout or blockade of CD43 in leukemia cells relieves Siglec-7-mediated inhibition of immune killing activity. This work identifies a potential target for immune checkpoint blockade therapy and represents a generalizable approach to dissection of glycan-receptor interactions in living cells.

High-Throughput Discovery and Characterization of Human Transcriptional Effectors.

Thousands of proteins localize to the nucleus; however, it remains unclear which contain transcriptional effectors. Here, we develop HT-recruit, a pooled assay where protein libraries are recruited to a reporter, and their transcriptional effects are measured by sequencing. Using this approach, we measure gene silencing and activation for thousands of domains. We find a relationship between repressor function and evolutionary age for the KRAB domains, discover that Homeodomain repressor strength is collinear with Hox genetic organization, and identify activities for several domains of unknown function. Deep mutational scanning of the CRISPRi KRAB maps the co-repressor binding surface and identifies substitutions that improve stability/silencing. By tiling 238 proteins, we find repressors as short as ten amino acids. Finally, we report new activator domains, including a divergent KRAB. These results provide a resource of 600 human proteins containing effectors and demonstrate a scalable strategy for assigning functions to protein domains.

LRRC8A:C/E Heteromeric Channels Are Ubiquitous Transporters of cGAMP.

Extracellular 2'3'-cyclic-GMP-AMP (cGAMP) is an immunotransmitter exported by diseased cells and imported into host cells to activate the innate immune STING pathway. We previously identified SLC19A1 as a cGAMP importer, but its use across human cell lines is limited. Here, we identify LRRC8A heteromeric channels, better known as volume-regulated anion channels (VRAC), as widely expressed cGAMP transporters. LRRC8A forms complexes with LRRC8C and/or LRRC8E, depending on their expression levels, to transport cGAMP and other 2'3'-cyclic dinucleotides. In contrast, LRRC8D inhibits cGAMP transport. We demonstrate that cGAMP is effluxed or influxed via LRRC8 channels, as dictated by the cGAMP electrochemical gradient. Activation of LRRC8A channels, which can occur under diverse stresses, strongly potentiates cGAMP transport. We identify activator sphingosine 1-phosphate and inhibitor DCPIB as chemical tools to manipulate channel-mediated cGAMP transport. Finally, LRRC8A channels are key cGAMP transporters in resting primary human vasculature cells and universal human cGAMP transporters when activated.

Transcriptomic signatures across human tissues identify functional rare genetic variation.

Rare genetic variants are abundant across the human genome, and identifying their function and phenotypic impact is a major challenge. Measuring aberrant gene expression has aided in identifying functional, large-effect rare variants (RVs). Here, we expanded detection of genetically driven transcriptome abnormalities by analyzing gene expression, allele-specific expression, and alternative splicing from multitissue RNA-sequencing data, and demonstrate that each signal informs unique classes of RVs. We developed Watershed, a probabilistic model that integrates multiple genomic and transcriptomic signals to predict variant function, validated these predictions in additional cohorts and through experimental assays, and used them to assess RVs in the UK Biobank, the Million Veterans Program, and the Jackson Heart Study. Our results link thousands of RVs to diverse molecular effects and provide evidence to associate RVs affecting the transcriptome with human traits.

Metabolic precision labeling enables selective probing of O-linked N-acetylgalactosamine glycosylation.

Protein glycosylation events that happen early in the secretory pathway are often dysregulated during tumorigenesis. These events can be probed, in principle, by monosaccharides with bioorthogonal tags that would ideally be specific for distinct glycan subtypes. However, metabolic interconversion into other monosaccharides drastically reduces such specificity in the living cell. Here, we use a structure-based design process to develop the monosaccharide probe N-(S)-azidopropionylgalactosamine (GalNAzMe) that is specific for cancer-relevant Ser/Thr(O)-linked N-acetylgalactosamine (GalNAc) glycosylation. By virtue of a branched N-acylamide side chain, GalNAzMe is not interconverted by epimerization to the corresponding N-acetylglucosamine analog by the epimerase N-acetylgalactosamine-4-epimerase (GALE) like conventional GalNAc-based probes. GalNAzMe enters O-GalNAc glycosylation but does not enter other major cell surface glycan types including Asn(N)-linked glycans. We transfect cells with the engineered pyrophosphorylase mut-AGX1 to biosynthesize the nucleotide-sugar donor uridine diphosphate (UDP)-GalNAzMe from a sugar-1-phosphate precursor. Tagged with a bioorthogonal azide group, GalNAzMe serves as an O-glycan-specific reporter in superresolution microscopy, chemical glycoproteomics, a genome-wide CRISPR-knockout (CRISPR-KO) screen, and imaging of intestinal organoids. Additional ectopic expression of an engineered glycosyltransferase, "bump-and-hole" (BH)-GalNAc-T2, boosts labeling in a programmable fashion by increasing incorporation of GalNAzMe into the cell surface glycoproteome. Alleviating the need for GALE-KO cells in metabolic labeling experiments, GalNAzMe is a precision tool that allows a detailed view into the biology of a major type of cancer-relevant protein glycosylation.

Retro-2 protects cells from ricin toxicity by inhibiting ASNA1-mediated ER targeting and insertion of tail-anchored proteins.

The small molecule Retro-2 prevents ricin toxicity through a poorly-defined mechanism of action (MOA), which involves halting retrograde vesicle transport to the endoplasmic reticulum (ER). CRISPRi genetic interaction analysis revealed Retro-2 activity resembles disruption of the transmembrane domain recognition complex (TRC) pathway, which mediates post-translational ER-targeting and insertion of tail-anchored (TA) proteins, including SNAREs required for retrograde transport. Cell-based and in vitro assays show that Retro-2 blocks delivery of newly-synthesized TA-proteins to the ER-targeting factor ASNA1 (TRC40). An ASNA1 point mutant identified using CRISPR-mediated mutagenesis abolishes both the cytoprotective effect of Retro-2 against ricin and its inhibitory effect on ASNA1-mediated ER-targeting. Together, our work explains how Retro-2 prevents retrograde trafficking of toxins by inhibiting TA-protein targeting, describes a general CRISPR strategy for predicting the MOA of small molecules, and paves the way for drugging the TRC pathway to treat broad classes of viruses known to be inhibited by Retro-2.

Mitigation of off-target toxicity in CRISPR-Cas9 screens for essential non-coding elements.

Pooled CRISPR-Cas9 screens are a powerful method for functionally characterizing regulatory elements in the non-coding genome, but off-target effects in these experiments have not been systematically evaluated. Here, we investigate Cas9, dCas9, and CRISPRi/a off-target activity in screens for essential regulatory elements. The sgRNAs with the largest effects in genome-scale screens for essential CTCF loop anchors in K562 cells were not single guide RNAs (sgRNAs) that disrupted gene expression near the on-target CTCF anchor. Rather, these sgRNAs had high off-target activity that, while only weakly correlated with absolute off-target site number, could be predicted by the recently developed GuideScan specificity score. Screens conducted in parallel with CRISPRi/a, which do not induce double-stranded DNA breaks, revealed that a distinct set of off-targets also cause strong confounding fitness effects with these epigenome-editing tools. Promisingly, filtering of CRISPRi libraries using GuideScan specificity scores removed these confounded sgRNAs and enabled identification of essential regulatory elements.

SLC19A1 Is an Importer of the Immunotransmitter cGAMP.

2'3'-cyclic-GMP-AMP (cGAMP) is a second messenger that activates the antiviral stimulator of interferon genes (STING) pathway. We recently identified a novel role for cGAMP as a soluble, extracellular immunotransmitter that is produced and secreted by cancer cells. Secreted cGAMP is then sensed by host cells, eliciting an antitumoral immune response. Due to the antitumoral effects of cGAMP, other CDN-based STING agonists are currently under investigation in clinical trials for metastatic solid tumors. However, it is unknown how cGAMP and other CDNs cross the cell membrane to activate intracellular STING. Using a genome-wide CRISPR screen, we identified SLC19A1 as the first known importer of cGAMP and other CDNs, including the investigational new drug 2'3'-bisphosphosphothioate-cyclic-di-AMP (2'3'-CDAS). These discoveries will provide insight into cGAMP's role as an immunotransmitter and aid in the development of more targeted CDN-based cancer therapeutics.

A CRISPR-based screen for Hedgehog signaling provides insights into ciliary function and ciliopathies.

Primary cilia organize Hedgehog signaling and shape embryonic development, and their dysregulation is the unifying cause of ciliopathies. We conducted a functional genomic screen for Hedgehog signaling by engineering antibiotic-based selection of Hedgehog-responsive cells and applying genome-wide CRISPR-mediated gene disruption. The screen can robustly identify factors required for ciliary signaling with few false positives or false negatives. Characterization of hit genes uncovered novel components of several ciliary structures, including a protein complex that contains δ-tubulin and ε-tubulin and is required for centriole maintenance. The screen also provides an unbiased tool for classifying ciliopathies and showed that many congenital heart disorders are caused by loss of ciliary signaling. Collectively, our study enables a systematic analysis of ciliary function and of ciliopathies, and also defines a versatile platform for dissecting signaling pathways through CRISPR-based screening.

The impact of rare variation on gene expression across tissues.

Rare genetic variants are abundant in humans and are expected to contribute to individual disease risk. While genetic association studies have successfully identified common genetic variants associated with susceptibility, these studies are not practical for identifying rare variants. Efforts to distinguish pathogenic variants from benign rare variants have leveraged the genetic code to identify deleterious protein-coding alleles, but no analogous code exists for non-coding variants. Therefore, ascertaining which rare variants have phenotypic effects remains a major challenge. Rare non-coding variants have been associated with extreme gene expression in studies using single tissues, but their effects across tissues are unknown. Here we identify gene expression outliers, or individuals showing extreme expression levels for a particular gene, across 44 human tissues by using combined analyses of whole genomes and multi-tissue RNA-sequencing data from the Genotype-Tissue Expression (GTEx) project v6p release. We find that 58% of underexpression and 28% of overexpression outliers have nearby conserved rare variants compared to 8% of non-outliers. Additionally, we developed RIVER (RNA-informed variant effect on regulation), a Bayesian statistical model that incorporates expression data to predict a regulatory effect for rare variants with higher accuracy than models using genomic annotations alone. Overall, we demonstrate that rare variants contribute to large gene expression changes across tissues and provide an integrative method for interpretation of rare variants in individual genomes.

Methods and Applications of CRISPR-Mediated Base Editing in Eukaryotic Genomes.

The past several years have seen an explosion in development of applications for the CRISPR-Cas9 system, from efficient genome editing, to high-throughput screening, to recruitment of a range of DNA and chromatin-modifying enzymes. While homology-directed repair (HDR) coupled with Cas9 nuclease cleavage has been used with great success to repair and re-write genomes, recently developed base-editing systems present a useful orthogonal strategy to engineer nucleotide substitutions. Base editing relies on recruitment of cytidine deaminases to introduce changes (rather than double-stranded breaks and donor templates) and offers potential improvements in efficiency while limiting damage and simplifying the delivery of editing machinery. At the same time, these systems enable novel mutagenesis strategies to introduce sequence diversity for engineering and discovery. Here, we review the different base-editing platforms, including their deaminase recruitment strategies and editing outcomes, and compare them to other CRISPR genome-editing technologies. Additionally, we discuss how these systems have been applied in therapeutic, engineering, and research settings. Lastly, we explore future directions of this emerging technology.

Static and Dynamic DNA Loops form AP-1-Bound Activation Hubs during Macrophage Development.

The three-dimensional arrangement of the human genome comprises a complex network of structural and regulatory chromatin loops important for coordinating changes in transcription during human development. To better understand the mechanisms underlying context-specific 3D chromatin structure and transcription during cellular differentiation, we generated comprehensive in situ Hi-C maps of DNA loops in human monocytes and differentiated macrophages. We demonstrate that dynamic looping events are regulatory rather than structural in nature and uncover widespread coordination of dynamic enhancer activity at preformed and acquired DNA loops. Enhancer-bound loop formation and enhancer activation of preformed loops together form multi-loop activation hubs at key macrophage genes. Activation hubs connect 3.4 enhancers per promoter and exhibit a strong enrichment for activator protein 1 (AP-1)-binding events, suggesting that multi-loop activation hubs involving cell-type-specific transcription factors represent an important class of regulatory chromatin structures for the spatiotemporal control of transcription.

Genome-scale measurement of off-target activity using Cas9 toxicity in high-throughput screens.

CRISPR-Cas9 screens are powerful tools for high-throughput interrogation of genome function, but can be confounded by nuclease-induced toxicity at both on- and off-target sites, likely due to DNA damage. Here, to test potential solutions to this issue, we design and analyse a CRISPR-Cas9 library with 10 variable-length guides per gene and thousands of negative controls targeting non-functional, non-genic regions (termed safe-targeting guides), in addition to non-targeting controls. We find this library has excellent performance in identifying genes affecting growth and sensitivity to the ricin toxin. The safe-targeting guides allow for proper control of toxicity from on-target DNA damage. Using this toxicity as a proxy to measure off-target cutting, we demonstrate with tens of thousands of guides both the nucleotide position-dependent sensitivity to single mismatches and the reduction of off-target cutting using truncated guides. Our results demonstrate a simple strategy for high-throughput evaluation of target specificity and nuclease toxicity in Cas9 screens.

Population- and individual-specific regulatory variation in Sardinia.

Genetic studies of complex traits have mainly identified associations with noncoding variants. To further determine the contribution of regulatory variation, we combined whole-genome and transcriptome data for 624 individuals from Sardinia to identify common and rare variants that influence gene expression and splicing. We identified 21,183 expression quantitative trait loci (eQTLs) and 6,768 splicing quantitative trait loci (sQTLs), including 619 new QTLs. We identified high-frequency QTLs and found evidence of selection near genes involved in malarial resistance and increased multiple sclerosis risk, reflecting the epidemiological history of Sardinia. Using family relationships, we identified 809 segregating expression outliers (median z score of 2.97), averaging 13.3 genes per individual. Outlier genes were enriched for proximal rare variants, providing a new approach to study large-effect regulatory variants and their relevance to traits. Our results provide insight into the effects of regulatory variants and their relationship to population history and individual genetic risk.

Synergistic drug combinations for cancer identified in a CRISPR screen for pairwise genetic interactions.

Identification of effective combination therapies is critical to address the emergence of drug-resistant cancers, but direct screening of all possible drug combinations is infeasible. Here we introduce a CRISPR-based double knockout (CDKO) system that improves the efficiency of combinatorial genetic screening using an effective strategy for cloning and sequencing paired single guide RNA (sgRNA) libraries and a robust statistical scoring method for calculating genetic interactions (GIs) from CRISPR-deleted gene pairs. We applied CDKO to generate a large-scale human GI map, comprising 490,000 double-sgRNAs directed against 21,321 pairs of drug targets in K562 leukemia cells and identified synthetic lethal drug target pairs for which corresponding drugs exhibit synergistic killing. These included the BCL2L1 and MCL1 combination, which was also effective in imatinib-resistant cells. We further validated this system by identifying known and previously unidentified GIs between modifiers of ricin toxicity. This work provides an effective strategy to screen synergistic drug combinations in high-throughput and a CRISPR-based tool to dissect functional GI networks.

Directed evolution using dCas9-targeted somatic hypermutation in mammalian cells.

Engineering and study of protein function by directed evolution has been limited by the technical requirement to use global mutagenesis or introduce DNA libraries. Here, we develop CRISPR-X, a strategy to repurpose the somatic hypermutation machinery for protein engineering in situ. Using catalytically inactive dCas9 to recruit variants of cytidine deaminase (AID) with MS2-modified sgRNAs, we can specifically mutagenize endogenous targets with limited off-target damage. This generates diverse libraries of localized point mutations and can target multiple genomic locations simultaneously. We mutagenize GFP and select for spectrum-shifted variants, including EGFP. Additionally, we mutate the target of the cancer therapeutic bortezomib, PSMB5, and identify known and novel mutations that confer bortezomib resistance. Finally, using a hyperactive AID variant, we mutagenize loci both upstream and downstream of transcriptional start sites. These experiments illustrate a powerful approach to create complex libraries of genetic variants in native context, which is broadly applicable to investigate and improve protein function.