CRISPRi screens identify the lncRNA, LOUP, as a multifunctional locus regulating macrophage differentiation and inflammatory signaling.

Long noncoding RNAs (lncRNAs) account for the largest portion of RNA from the transcriptome, yet most of their functions remain unknown. Here, we performed two independent high-throughput CRISPRi screens to understand the role of lncRNAs in monocyte function and differentiation. The first was a reporter-based screen to identify lncRNAs that regulate TLR4-NFkB signaling in human monocytes and the second screen identified lncRNAs involved in monocyte to macrophage differentiation. We successfully identified numerous noncoding and protein-coding genes that can positively or negatively regulate inflammation and differentiation. To understand the functional roles of lncRNAs in both processes, we chose to further study the lncRNA LOUP [lncRNA originating from upstream regulatory element of SPI1 (also known as PU.1)], as it emerged as a top hit in both screens. Not only does LOUP regulate its neighboring gene, the myeloid fate-determining factor SPI1, thereby affecting monocyte to macrophage differentiation, but knockdown of LOUP leads to a broad upregulation of NFkB-targeted genes at baseline and upon TLR4-NFkB activation. LOUP also harbors three small open reading frames capable of being translated and are responsible for LOUP's ability to negatively regulate TLR4/NFkB signaling. This work emphasizes the value of high-throughput screening to rapidly identify functional lncRNAs in the innate immune system.

FET fusion oncoproteins disrupt physiologic DNA repair networks and induce ATR synthetic lethality in cancer.

The genetic principle of synthetic lethality is clinically validated in cancers with loss of specific DNA damage response (DDR) pathway genes (i.e. BRCA1/2 tumor suppressor mutations). The broader question of whether and how oncogenes create tumor-specific vulnerabilities within DDR networks remains unanswered. Native FET protein family members are among the earliest proteins recruited to DNA double-strand breaks (DSBs) during the DDR, though the function of both native FET proteins and FET fusion oncoproteins in DSB repair remains poorly defined. Here we focus on Ewing sarcoma (ES), an EWS-FLI1 fusion oncoprotein-driven pediatric bone tumor, as a model for FET rearranged cancers. We discover that the EWS-FLI1 fusion oncoprotein is recruited to DNA DSBs and interferes with native EWS function in activating the DNA damage sensor ATM. Using preclinical mechanistic approaches and clinical datasets, we establish functional ATM deficiency as a principal DNA repair defect in ES and the compensatory ATR signaling axis as a collateral dependency and therapeutic target in FET rearranged cancers. Thus, aberrant recruitment of a fusion oncoprotein to sites of DNA damage can disrupt normal DSB repair, revealing a mechanism for how oncogenes can create cancer-specific synthetic lethality within DDR networks.

FET fusion oncoproteins disrupt physiologic DNA repair networks in cancer.

While oncogenes promote cancer cell growth, unrestrained proliferation represents a significant stressor to cellular homeostasis networks such as the DNA damage response (DDR). To enable oncogene tolerance, many cancers disable tumor suppressive DDR signaling through genetic loss of DDR pathways and downstream effectors (e.g., ATM or p53 tumor suppressor mutations). Whether and how oncogenes can help "self-tolerize" by creating analogous functional deficiencies in physiologic DDR networks is not known. Here we focus on Ewing sarcoma, a FET fusion oncoprotein (EWS-FLI1) driven pediatric bone tumor, as a model for the class of FET rearranged cancers. Native FET protein family members are among the earliest factors recruited to DNA double-strand breaks (DSBs) during the DDR, though the function of both native FET proteins and FET fusion oncoproteins in DNA repair remains to be defined. Using preclinical mechanistic studies of the DDR and clinical genomic datasets from patient tumors, we discover that the EWS-FLI1 fusion oncoprotein is recruited to DNA DSBs and interferes with native FET (EWS) protein function in activating the DNA damage sensor ATM. As a consequence of FET fusion-mediated interference with the DDR, we establish functional ATM deficiency as the principal DNA repair defect in Ewing sarcoma and the compensatory ATR signaling axis as a collateral dependency and therapeutic target in multiple FET rearranged cancers. More generally, we find that aberrant recruitment of a fusion oncoprotein to sites of DNA damage can disrupt physiologic DSB repair, revealing a mechanism for how growth-promoting oncogenes can also create a functional deficiency within tumor suppressive DDR networks.

Single-cell multi-scale footprinting reveals the modular organization of DNA regulatory elements.

Cis-regulatory elements control gene expression and are dynamic in their structure, reflecting changes to the composition of diverse effector proteins over time1-3. Here we sought to connect the structural changes at cis-regulatory elements to alterations in cellular fate and function. To do this we developed PRINT, a computational method that uses deep learning to correct sequence bias in chromatin accessibility data and identifies multi-scale footprints of DNA-protein interactions. We find that multi-scale footprints enable more accurate inference of TF and nucleosome binding. Using PRINT with single-cell multi-omics, we discover wide-spread changes to the structure and function of candidate cis-regulatory elements (cCREs) across hematopoiesis, wherein nucleosomes slide, expose DNA for TF binding, and promote gene expression. Activity segmentation using the co-variance across cell states identifies "sub-cCREs" as modular cCRE subunits of regulatory DNA. We apply this single-cell and PRINT approach to characterize the age-associated alterations to cCREs within hematopoietic stem cells (HSCs). Remarkably, we find a spectrum of aging alterations among HSCs corresponding to a global gain of sub-cCRE activity while preserving cCRE accessibility. Collectively, we reveal the functional importance of cCRE structure across cell states, highlighting changes to gene regulation at single-cell and single-base-pair resolution.

Dual genome-wide coding and lncRNA screens in neural induction of induced pluripotent stem cells.

Human chromosomes are pervasively transcribed, but systematic understanding of coding and lncRNA genome function in cell differentiation is lacking. Using CRISPR interference (CRISPRi) in human induced pluripotent stem cells, we performed dual genome-wide screens - assessing 18,905 protein-coding and 10,678 lncRNA loci - and identified 419 coding and 201 lncRNA genes that regulate neural induction. Integrative analyses revealed distinct properties of coding and lncRNA genome function, including a 10-fold enrichment of lncRNA genes for roles in differentiation compared to proliferation. Further, we applied Perturb-seq to obtain granular insights into neural induction phenotypes. While most coding hits stalled or aborted differentiation, lncRNA hits were enriched for the genesis of diverse cellular states, including those outside the neural lineage. In addition to providing a rich resource (danlimlab.shinyapps.io/dualgenomewide) for understanding coding and lncRNA gene function in development, these results indicate that the lncRNA genome regulates lineage commitment in a manner fundamentally distinct from coding genes.

GATOR2-dependent mTORC1 activity is a therapeutic vulnerability in FOXO1 fusion-positive rhabdomyosarcoma.

Oncogenic FOXO1 gene fusions drive a subset of rhabdomyosarcoma (RMS) with poor survival; to date, these cancer drivers are therapeutically intractable. To identify new therapies for this disease, we undertook an isogenic CRISPR-interference screen to define PAX3-FOXO1-specific genetic dependencies and identified genes in the GATOR2 complex. GATOR2 loss in RMS abrogated aa-induced lysosomal localization of mTORC1 and consequent downstream signaling, slowing G1-S cell cycle transition. In vivo suppression of GATOR2 impaired the growth of tumor xenografts and favored the outgrowth of cells lacking PAX3-FOXO1. Loss of a subset of GATOR2 members can be compensated by direct genetic activation of mTORC1. RAS mutations are also sufficient to decouple mTORC1 activation from GATOR2, and indeed, fusion-negative RMS harboring such mutations exhibit aa-independent mTORC1 activity. A bisteric, mTORC1-selective small molecule induced tumor regressions in fusion-positive patient-derived tumor xenografts. These findings highlight a vulnerability in FOXO1 fusion-positive RMS and provide rationale for the clinical evaluation of bisteric mTORC1 inhibitors, currently in phase I testing, to treat this disease. Isogenic genetic screens can, thus, identify potentially exploitable vulnerabilities in fusion-driven pediatric cancers that otherwise remain mostly undruggable.

Genome-wide CRISPRi screening identifies OCIAD1 as a prohibitin client and regulatory determinant of mitochondrial Complex III assembly in human cells.

Dysfunction of the mitochondrial electron transport chain (mETC) is a major cause of human mitochondrial diseases. To identify determinants of mETC function, we screened a genome-wide human CRISPRi library under oxidative metabolic conditions with selective inhibition of mitochondrial Complex III and identified ovarian carcinoma immunoreactive antigen (OCIA) domain-containing protein 1 (OCIAD1) as a Complex III assembly factor. We find that OCIAD1 is an inner mitochondrial membrane protein that forms a complex with supramolecular prohibitin assemblies. Our data indicate that OCIAD1 is required for maintenance of normal steady-state levels of Complex III and the proteolytic processing of the catalytic subunit cytochrome c1 (CYC1). In OCIAD1 depleted mitochondria, unprocessed CYC1 is hemylated and incorporated into Complex III. We propose that OCIAD1 acts as an adaptor within prohibitin assemblies to stabilize and/or chaperone CYC1 and to facilitate its proteolytic processing by the IMMP2L protease.

High-content imaging-based pooled CRISPR screens in mammalian cells.

CRISPR (clustered regularly interspaced short palindromic repeats)-based gene inactivation provides a powerful means for linking genes to particular cellular phenotypes. CRISPR-based screening typically uses large genomic pools of single guide RNAs (sgRNAs). However, this approach is limited to phenotypes that can be enriched by chemical selection or FACS sorting. Here, we developed a microscopy-based approach, which we name optical enrichment, to select cells displaying a particular CRISPR-induced phenotype by automated imaging-based computation, mark them by photoactivation of an expressed photoactivatable fluorescent protein, and then isolate the fluorescent cells using fluorescence-activated cell sorting (FACS). A plugin was developed for the open source software µManager to automate the phenotypic identification and photoactivation of cells, allowing ∼1.5 million individual cells to be screened in 8 h. We used this approach to screen 6,092 sgRNAs targeting 544 genes for their effects on nuclear size regulation and identified 14 bona fide hits. These results present a scalable approach to facilitate imaging-based pooled CRISPR screens.

Genome-Scale Perturbation of Long Noncoding RNA Expression Using CRISPR Interference.

CRISPR-mediated interference (CRISPRi), a robust and specific system for programmably repressing transcription, provides a versatile tool for systematically characterizing the function of long noncoding RNAs (lncRNAs). When used with highly parallel, lentiviral pooled screening approaches, CRISPRi enables the targeted knockdown of tens of thousands of lncRNA-expressing loci in a single screen. Here we describe the use of CRISPRi to target lncRNA loci in a pooled screen, using cell growth and proliferation as an example of a phenotypic readout. Considerations for custom lncRNA-targeting libraries, alternative phenotypic readouts, and orthogonal validation approaches are also discussed.

Pharmaceutical-Grade Rigosertib Is a Microtubule-Destabilizing Agent.

We recently used CRISPRi/a-based chemical-genetic screens and cell biological, biochemical, and structural assays to determine that rigosertib, an anti-cancer agent in phase III clinical trials, kills cancer cells by destabilizing microtubules. Reddy and co-workers (Baker et al., 2020, this issue of Molecular Cell) suggest that a contaminating degradation product in commercial formulations of rigosertib is responsible for the microtubule-destabilizing activity. Here, we demonstrate that cells treated with pharmaceutical-grade rigosertib (>99.9% purity) or commercially obtained rigosertib have qualitatively indistinguishable phenotypes across multiple assays. The two formulations have indistinguishable chemical-genetic interactions with genes that modulate microtubule stability, both destabilize microtubules in cells and in vitro, and expression of a rationally designed tubulin mutant with a mutation in the rigosertib binding site (L240F TUBB) allows cells to proliferate in the presence of either formulation. Importantly, the specificity of the L240F TUBB mutant for microtubule-destabilizing agents has been confirmed independently. Thus, rigosertib kills cancer cells by destabilizing microtubules, in agreement with our original findings.

CRISPRi-based radiation modifier screen identifies long non-coding RNA therapeutic targets in glioma.

BACKGROUND: Long non-coding RNAs (lncRNAs) exhibit highly cell type-specific expression and function, making this class of transcript attractive for targeted cancer therapy. However, the vast majority of lncRNAs have not been tested as potential therapeutic targets, particularly in the context of currently used cancer treatments. Malignant glioma is rapidly fatal, and ionizing radiation is part of the current standard-of-care used to slow tumor growth in both adult and pediatric patients. RESULTS: We use CRISPR interference (CRISPRi) to screen 5689 lncRNA loci in human glioblastoma (GBM) cells, identifying 467 hits that modify cell growth in the presence of clinically relevant doses of fractionated radiation. Thirty-three of these lncRNA hits sensitize cells to radiation, and based on their expression in adult and pediatric gliomas, nine of these hits are prioritized as lncRNA Glioma Radiation Sensitizers (lncGRS). Knockdown of lncGRS-1, a primate-conserved, nuclear-enriched lncRNA, inhibits the growth and proliferation of primary adult and pediatric glioma cells, but not the viability of normal brain cells. Using human brain organoids comprised of mature neural cell types as a three-dimensional tissue substrate to model the invasive growth of glioma, we find that antisense oligonucleotides targeting lncGRS-1 selectively decrease tumor growth and sensitize glioma cells to radiation therapy. CONCLUSIONS: These studies identify lncGRS-1 as a glioma-specific therapeutic target and establish a generalizable approach to rapidly identify novel therapeutic targets in the vast non-coding genome to enhance radiation therapy.

Titrating gene expression using libraries of systematically attenuated CRISPR guide RNAs.

A lack of tools to precisely control gene expression has limited our ability to evaluate relationships between expression levels and phenotypes. Here, we describe an approach to titrate expression of human genes using CRISPR interference and series of single-guide RNAs (sgRNAs) with systematically modulated activities. We used large-scale measurements across multiple cell models to characterize activities of sgRNAs containing mismatches to their target sites and derived rules governing mismatched sgRNA activity using deep learning. These rules enabled us to synthesize a compact sgRNA library to titrate expression of ~2,400 genes essential for robust cell growth and to construct an in silico sgRNA library spanning the human genome. Staging cells along a continuum of gene expression levels combined with single-cell RNA-seq readout revealed sharp transitions in cellular behaviors at gene-specific expression thresholds. Our work provides a general tool to control gene expression, with applications ranging from tuning biochemical pathways to identifying suppressors for diseases of dysregulated gene expression.

New factors for protein transport identified by a genome-wide CRISPRi screen in mammalian cells.

Protein and membrane trafficking pathways are critical for cell and tissue homeostasis. Traditional genetic and biochemical approaches have shed light on basic principles underlying these processes. However, the list of factors required for secretory pathway function remains incomplete, and mechanisms involved in their adaptation poorly understood. Here, we present a powerful strategy based on a pooled genome-wide CRISPRi screen that allowed the identification of new factors involved in protein transport. Two newly identified factors, TTC17 and CCDC157, localized along the secretory pathway and were found to interact with resident proteins of ER-Golgi membranes. In addition, we uncovered that upon TTC17 knockdown, the polarized organization of Golgi cisternae was altered, creating glycosylation defects, and that CCDC157 is an important factor for the fusion of transport carriers to Golgi membranes. In conclusion, our work identified and characterized new actors in the mechanisms of protein transport and secretion and opens stimulating perspectives for the use of our platform in physiological and pathological contexts.

Exploring genetic interaction manifolds constructed from rich single-cell phenotypes.

How cellular and organismal complexity emerges from combinatorial expression of genes is a central question in biology. High-content phenotyping approaches such as Perturb-seq (single-cell RNA-sequencing pooled CRISPR screens) present an opportunity for exploring such genetic interactions (GIs) at scale. Here, we present an analytical framework for interpreting high-dimensional landscapes of cell states (manifolds) constructed from transcriptional phenotypes. We applied this approach to Perturb-seq profiling of strong GIs mined from a growth-based, gain-of-function GI map. Exploration of this manifold enabled ordering of regulatory pathways, principled classification of GIs (e.g., identifying suppressors), and mechanistic elucidation of synergistic interactions, including an unexpected synergy between CBL and CNN1 driving erythroid differentiation. Finally, we applied recommender system machine learning to predict interactions, facilitating exploration of vastly larger GI manifolds.

Cellular response to small molecules that selectively stall protein synthesis by the ribosome.

Identifying small molecules that inhibit protein synthesis by selectively stalling the ribosome constitutes a new strategy for therapeutic development. Compounds that inhibit the translation of PCSK9, a major regulator of low-density lipoprotein cholesterol, have been identified that reduce LDL cholesterol in preclinical models and that affect the translation of only a few off-target proteins. Although some of these compounds hold potential for future therapeutic development, it is not known how they impact the physiology of cells or ribosome quality control pathways. Here we used a genome-wide CRISPRi screen to identify proteins and pathways that modulate cell growth in the presence of high doses of a selective PCSK9 translational inhibitor, PF-06378503 (PF8503). The two most potent genetic modifiers of cell fitness in the presence of PF8503, the ubiquitin binding protein ASCC2 and helicase ASCC3, bind to the ribosome and protect cells from toxic effects of high concentrations of the compound. Surprisingly, translation quality control proteins Pelota (PELO) and HBS1L sensitize cells to PF8503 treatment. In genetic interaction experiments, ASCC3 acts together with ASCC2, and functions downstream of HBS1L. Taken together, these results identify new connections between ribosome quality control pathways, and provide new insights into the selectivity of compounds that stall human translation that will aid the development of next-generation selective translation stalling compounds to treat disease.

A high-throughput screen of real-time ATP levels in individual cells reveals mechanisms of energy failure.

Insufficient or dysregulated energy metabolism may underlie diverse inherited and degenerative diseases, cancer, and even aging itself. ATP is the central energy carrier in cells, but critical pathways for regulating ATP levels are not systematically understood. We combined a pooled clustered regularly interspaced short palindromic repeats interference (CRISPRi) library enriched for mitochondrial genes, a fluorescent biosensor, and fluorescence-activated cell sorting (FACS) in a high-throughput genetic screen to assay ATP concentrations in live human cells. We identified genes not known to be involved in energy metabolism. Most mitochondrial ribosomal proteins are essential in maintaining ATP levels under respiratory conditions, and impaired respiration predicts poor growth. We also identified genes for which coenzyme Q10 (CoQ10) supplementation rescued ATP deficits caused by knockdown. These included CoQ10 biosynthetic genes associated with human disease and a subset of genes not linked to CoQ10 biosynthesis, indicating that increasing CoQ10 can preserve ATP in specific genetic contexts. This screening paradigm reveals mechanisms of metabolic control and genetic defects responsive to energy-based therapies.

Mapping the Genetic Landscape of Human Cells.

Seminal yeast studies have established the value of comprehensively mapping genetic interactions (GIs) for inferring gene function. Efforts in human cells using focused gene sets underscore the utility of this approach, but the feasibility of generating large-scale, diverse human GI maps remains unresolved. We developed a CRISPR interference platform for large-scale quantitative mapping of human GIs. We systematically perturbed 222,784 gene pairs in two cancer cell lines. The resultant maps cluster functionally related genes, assigning function to poorly characterized genes, including TMEM261, a new electron transport chain component. Individual GIs pinpoint unexpected relationships between pathways, exemplified by a specific cholesterol biosynthesis intermediate whose accumulation induces deoxynucleotide depletion, causing replicative DNA damage and a synthetic-lethal interaction with the ATR/9-1-1 DNA repair pathway. Our map provides a broad resource, establishes GI maps as a high-resolution tool for dissecting gene function, and serves as a blueprint for mapping the genetic landscape of human cells.

Promoter of lncRNA Gene PVT1 Is a Tumor-Suppressor DNA Boundary Element.

Noncoding mutations in cancer genomes are frequent but challenging to interpret. PVT1 encodes an oncogenic lncRNA, but recurrent translocations and deletions in human cancers suggest alternative mechanisms. Here, we show that the PVT1 promoter has a tumor-suppressor function that is independent of PVT1 lncRNA. CRISPR interference of PVT1 promoter enhances breast cancer cell competition and growth in vivo. The promoters of the PVT1 and the MYC oncogenes, located 55 kb apart on chromosome 8q24, compete for engagement with four intragenic enhancers in the PVT1 locus, thereby allowing the PVT1 promoter to regulate pause release of MYC transcription. PVT1 undergoes developmentally regulated monoallelic expression, and the PVT1 promoter inhibits MYC expression only from the same chromosome via promoter competition. Cancer genome sequencing identifies recurrent mutations encompassing the human PVT1 promoter, and genome editing verified that PVT1 promoter mutation promotes cancer cell growth. These results highlight regulatory sequences of lncRNA genes as potential disease-associated DNA elements.

Identification of a transporter complex responsible for the cytosolic entry of nitrogen-containing bisphosphonates.

Nitrogen-containing-bisphosphonates (N-BPs) are a class of drugs widely prescribed to treat osteoporosis and other bone-related diseases. Although previous studies have established that N-BPs function by inhibiting the mevalonate pathway in osteoclasts, the mechanism by which N-BPs enter the cytosol from the extracellular space to reach their molecular target is not understood. Here, we implemented a CRISPRi-mediated genome-wide screen and identified SLC37A3 (solute carrier family 37 member A3) as a gene required for the action of N-BPs in mammalian cells. We observed that SLC37A3 forms a complex with ATRAID (all-trans retinoic acid-induced differentiation factor), a previously identified genetic target of N-BPs. SLC37A3 and ATRAID localize to lysosomes and are required for releasing N-BP molecules that have trafficked to lysosomes through fluid-phase endocytosis into the cytosol. Our results elucidate the route by which N-BPs are delivered to their molecular target, addressing a key aspect of the mechanism of action of N-BPs that may have significant clinical relevance.