Brainwide silencing of prion protein by AAV-mediated delivery of an engineered compact epigenetic editor.

Prion disease is caused by misfolding of the prion protein (PrP) into pathogenic self-propagating conformations, leading to rapid-onset dementia and death. However, elimination of endogenous PrP halts prion disease progression. In this study, we describe Coupled Histone tail for Autoinhibition Release of Methyltransferase (CHARM), a compact, enzyme-free epigenetic editor capable of silencing transcription through programmable DNA methylation. Using a histone H3 tail-Dnmt3l fusion, CHARM recruits and activates endogenous DNA methyltransferases, thereby reducing transgene size and cytotoxicity. When delivered to the mouse brain by systemic injection of adeno-associated virus (AAV), Prnp-targeted CHARM ablates PrP expression across the brain. Furthermore, we have temporally limited editor expression by implementing a kinetically tuned self-silencing approach. CHARM potentially represents a broadly applicable strategy to suppress pathogenic proteins, including those implicated in other neurodegenerative diseases.

Nucleotide depletion promotes cell fate transitions by inducing DNA replication stress.

Control of cellular identity requires coordination of developmental programs with environmental factors such as nutrient availability, suggesting that perturbing metabolism can alter cell state. Here, we find that nucleotide depletion and DNA replication stress drive differentiation in human and murine normal and transformed hematopoietic systems, including patient-derived acute myeloid leukemia (AML) xenografts. These cell state transitions begin during S phase and are independent of ATR/ATM checkpoint signaling, double-stranded DNA break formation, and changes in cell cycle length. In systems where differentiation is blocked by oncogenic transcription factor expression, replication stress activates primed regulatory loci and induces lineage-appropriate maturation genes despite the persistence of progenitor programs. Altering the baseline cell state by manipulating transcription factor expression causes replication stress to induce genes specific for alternative lineages. The ability of replication stress to selectively activate primed maturation programs across different contexts suggests a general mechanism by which changes in metabolism can promote lineage-appropriate cell state transitions.

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.

Deciphering cell states and genealogies of human haematopoiesis.

The human blood system is maintained through the differentiation and massive amplification of a limited number of long-lived haematopoietic stem cells (HSCs)1. Perturbations to this process underlie diverse diseases, but the clonal contributions to human haematopoiesis and how this changes with age remain incompletely understood. Although recent insights have emerged from barcoding studies in model systems2-5, simultaneous detection of cell states and phylogenies from natural barcodes in humans remains challenging. Here we introduce an improved, single-cell lineage-tracing system based on deep detection of naturally occurring mitochondrial DNA mutations with simultaneous readout of transcriptional states and chromatin accessibility. We use this system to define the clonal architecture of HSCs and map the physiological state and output of clones. We uncover functional heterogeneity in HSC clones, which is stable over months and manifests as both differences in total HSC output and biases towards the production of different mature cell types. We also find that the diversity of HSC clones decreases markedly with age, leading to an oligoclonal structure with multiple distinct clonal expansions. Our study thus provides a clonally resolved and cell-state-aware atlas of human haematopoiesis at single-cell resolution, showing an unappreciated functional diversity of human HSC clones and, more broadly, paving the way for refined studies of clonal dynamics across a range of tissues in human health and disease.

Systematic functional interrogation of SARS-CoV-2 host factors using Perturb-seq.

Genomic and proteomic screens have identified numerous host factors of SARS-CoV-2, but efficient delineation of their molecular roles during infection remains a challenge. Here we use Perturb-seq, combining genetic perturbations with a single-cell readout, to investigate how inactivation of host factors changes the course of SARS-CoV-2 infection and the host response in human lung epithelial cells. Our high-dimensional data resolve complex phenotypes such as shifts in the stages of infection and modulations of the interferon response. However, only a small percentage of host factors showed such phenotypes upon perturbation. We further identified the NF-κB inhibitor IκBα (NFKBIA), as well as the translation factors EIF4E2 and EIF4H as strong host dependency factors acting early in infection. Overall, our study provides massively parallel functional characterization of host factors of SARS-CoV-2 and quantitatively defines their roles both in virus-infected and bystander cells.

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.

Microfluidics-free single-cell genomics with templated emulsification.

Current single-cell RNA-sequencing approaches have limitations that stem from the microfluidic devices or fluid handling steps required for sample processing. We develop a method that does not require specialized microfluidic devices, expertise or hardware. Our approach is based on particle-templated emulsification, which allows single-cell encapsulation and barcoding of cDNA in uniform droplet emulsions with only a vortexer. Particle-templated instant partition sequencing (PIP-seq) accommodates a wide range of emulsification formats, including microwell plates and large-volume conical tubes, enabling thousands of samples or millions of cells to be processed in minutes. We demonstrate that PIP-seq produces high-purity transcriptomes in mouse-human mixing studies, is compatible with multiomics measurements and can accurately characterize cell types in human breast tissue compared to a commercial microfluidic platform. Single-cell transcriptional profiling of mixed phenotype acute leukemia using PIP-seq reveals the emergence of heterogeneity within chemotherapy-resistant cell subsets that were hidden by standard immunophenotyping. PIP-seq is a simple, flexible and scalable next-generation workflow that extends single-cell sequencing to new applications.

MTCH2 is a mitochondrial outer membrane protein insertase.

In the mitochondrial outer membrane, α-helical transmembrane proteins play critical roles in cytoplasmic-mitochondrial communication. Using genome-wide CRISPR screens, we identified mitochondrial carrier homolog 2 (MTCH2), and its paralog MTCH1, and showed that it is required for insertion of biophysically diverse tail-anchored (TA), signal-anchored, and multipass proteins, but not outer membrane ß-barrel proteins. Purified MTCH2 was sufficient to mediate insertion into reconstituted proteoliposomes. Functional and mutational studies suggested that MTCH2 has evolved from a solute carrier transporter. MTCH2 uses membrane-embedded hydrophilic residues to function as a gatekeeper for the outer membrane, controlling mislocalization of TAs into the endoplasmic reticulum and modulating the sensitivity of leukemia cells to apoptosis. Our identification of MTCH2 as an insertase provides a mechanistic explanation for the diverse phenotypes and disease states associated with MTCH2 dysfunction.

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.

Cellular barcoding to decipher clonal dynamics in disease.

Cellular barcodes are distinct DNA sequences that enable one to track specific cells across time or space. Recent advances in our ability to detect natural or synthetic cellular barcodes, paired with single-cell readouts of cell state, have markedly increased our knowledge of clonal dynamics and genealogies of the cells that compose a variety of tissues and organs. These advances hold promise to redefine our view of human disease. Here, we provide an overview of cellular barcoding approaches, discuss applications to gain new insights into disease mechanisms, and provide an outlook on future applications. We discuss unanticipated insights gained through barcoding in studies of cancer and blood cell production and describe how barcoding can be applied to a growing array of medical fields, particularly with the increasing recognition of clonal contributions in human diseases.

Adaptive exchange sustains cullin-RING ubiquitin ligase networks and proper licensing of DNA replication.

Cop9 signalosome (CSN) regulates the function of cullin-RING E3 ubiquitin ligases (CRLs) by deconjugating the ubiquitin-like protein NEDD8 from the cullin subunit. To understand the physiological impact of CSN function on the CRL network and cell proliferation, we combined quantitative mass spectrometry and genome-wide CRISPR interference (CRISPRi) and CRISPR activation (CRISPRa) screens to identify factors that modulate cell viability upon inhibition of CSN by the small molecule CSN5i-3. CRL components and regulators strongly modulated the antiproliferative effects of CSN5i-3, and in addition we found two pathways involved in genome integrity, SCFFBXO5-APC/C-GMNN and CUL4DTL-SETD8, that contribute substantially to the toxicity of CSN inhibition. Our data highlight the importance of CSN-mediated NEDD8 deconjugation and adaptive exchange of CRL substrate receptors in sustaining CRL function and suggest approaches for leveraging CSN inhibition for the treatment of cancer.

Ras-mutant cancers are sensitive to small molecule inhibition of V-type ATPases in mice.

Mutations in Ras family proteins are implicated in 33% of human cancers, but direct pharmacological inhibition of Ras mutants remains challenging. As an alternative to direct inhibition, we screened for sensitivities in Ras-mutant cells and discovered 249C as a Ras-mutant selective cytotoxic agent with nanomolar potency against a spectrum of Ras-mutant cancers. 249C binds to vacuolar (V)-ATPase with nanomolar affinity and inhibits its activity, preventing lysosomal acidification and inhibiting autophagy and macropinocytosis pathways that several Ras-driven cancers rely on for survival. Unexpectedly, potency of 249C varies with the identity of the Ras driver mutation, with the highest potency for KRASG13D and G12V both in vitro and in vivo, highlighting a mutant-specific dependence on macropinocytosis and lysosomal pH. Indeed, 249C potently inhibits tumor growth without adverse side effects in mouse xenografts of KRAS-driven lung and colon cancers. A comparison of isogenic SW48 xenografts with different KRAS mutations confirmed that KRASG13D/+ (followed by G12V/+) mutations are especially sensitive to 249C treatment. These data establish proof-of-concept for targeting V-ATPase in cancers driven by specific KRAS mutations such as KRASG13D and G12V.

A campaign targeting a conserved Hsp70 binding site uncovers how subcellular localization is linked to distinct biological activities.

The potential of small molecules to localize within subcellular compartments is rarely explored. To probe this question, we measured the localization of Hsp70 inhibitors using fluorescence microscopy. We found that even closely related analogs had dramatically different distributions, with some residing predominantly in the mitochondria and others in the ER. CRISPRi screens supported this idea, showing that different compounds had distinct chemogenetic interactions with Hsp70s of the ER (HSPA5/BiP) and mitochondria (HSPA9/mortalin) and their co-chaperones. Moreover, localization seemed to determine function, even for molecules with conserved binding sites. Compounds with distinct partitioning have distinct anti-proliferative activity in breast cancer cells compared with anti-viral activity in cellular models of Dengue virus replication, likely because different sets of Hsp70s are required in these processes. These findings highlight the contributions of subcellular partitioning and chemogenetic interactions to small molecule activity, features that are rarely explored during medicinal chemistry campaigns.

Variant to function mapping at single-cell resolution through network propagation.

Genome-wide association studies in combination with single-cell genomic atlases can provide insights into the mechanisms of disease-causal genetic variation. However, identification of disease-relevant or trait-relevant cell types, states and trajectories is often hampered by sparsity and noise, particularly in the analysis of single-cell epigenomic data. To overcome these challenges, we present SCAVENGE, a computational algorithm that uses network propagation to map causal variants to their relevant cellular context at single-cell resolution. We demonstrate how SCAVENGE can help identify key biological mechanisms underlying human genetic variation, applying the method to blood traits at distinct stages of human hematopoiesis, to monocyte subsets that increase the risk for severe Coronavirus Disease 2019 (COVID-19) and to intermediate lymphocyte developmental states that predispose to acute leukemia. Our approach not only provides a framework for enabling variant-to-function insights at single-cell resolution but also suggests a more general strategy for maximizing the inferences that can be made using single-cell genomic data.

Lineage tracing reveals the phylodynamics, plasticity, and paths of tumor evolution.

Tumor evolution is driven by the progressive acquisition of genetic and epigenetic alterations that enable uncontrolled growth and expansion to neighboring and distal tissues. The study of phylogenetic relationships between cancer cells provides key insights into these processes. Here, we introduced an evolving lineage-tracing system with a single-cell RNA-seq readout into a mouse model of Kras;Trp53(KP)-driven lung adenocarcinoma and tracked tumor evolution from single-transformed cells to metastatic tumors at unprecedented resolution. We found that the loss of the initial, stable alveolar-type2-like state was accompanied by a transient increase in plasticity. This was followed by the adoption of distinct transcriptional programs that enable rapid expansion and, ultimately, clonal sweep of stable subclones capable of metastasizing. Finally, tumors develop through stereotypical evolutionary trajectories, and perturbing additional tumor suppressors accelerates progression by creating novel trajectories. Our study elucidates the hierarchical nature of tumor evolution and, more broadly, enables in-depth studies of tumor progression.

Mapping transcriptomic vector fields of single cells.

Single-cell (sc)RNA-seq, together with RNA velocity and metabolic labeling, reveals cellular states and transitions at unprecedented resolution. Fully exploiting these data, however, requires kinetic models capable of unveiling governing regulatory functions. Here, we introduce an analytical framework dynamo (https://github.com/aristoteleo/dynamo-release), which infers absolute RNA velocity, reconstructs continuous vector fields that predict cell fates, employs differential geometry to extract underlying regulations, and ultimately predicts optimal reprogramming paths and perturbation outcomes. We highlight dynamo's power to overcome fundamental limitations of conventional splicing-based RNA velocity analyses to enable accurate velocity estimations on a metabolically labeled human hematopoiesis scRNA-seq dataset. Furthermore, differential geometry analyses reveal mechanisms driving early megakaryocyte appearance and elucidate asymmetrical regulation within the PU.1-GATA1 circuit. Leveraging the least-action-path method, dynamo accurately predicts drivers of numerous hematopoietic transitions. Finally, in silico perturbations predict cell-fate diversions induced by gene perturbations. Dynamo, thus, represents an important step in advancing quantitative and predictive theories of cell-state transitions.

Variant to function mapping at single-cell resolution through network propagation.

With burgeoning human disease genetic associations and single-cell genomic atlases covering a range of tissues, there are unprecedented opportunities to systematically gain insights into the mechanisms of disease-causal variation. However, sparsity and noise, particularly in the context of single-cell epigenomic data, hamper the identification of disease- or trait-relevant cell types, states, and trajectories. To overcome these challenges, we have developed the SCAVENGE method, which maps causal variants to their relevant cellular context at single-cell resolution by employing the strategy of network propagation. We demonstrate how SCAVENGE can help identify key biological mechanisms underlying human genetic variation including enrichment of blood traits at distinct stages of human hematopoiesis, defining monocyte subsets that increase the risk for severe coronavirus disease 2019 (COVID-19), and identifying intermediate lymphocyte developmental states that are critical for predisposition to acute leukemia. Our approach not only provides a framework for enabling variant-to-function insights at single-cell resolution, but also suggests a more general strategy for maximizing the inferences that can be made using single-cell genomic data.

BRD2 inhibition blocks SARS-CoV-2 infection by reducing transcription of the host cell receptor ACE2.

SARS-CoV-2 infection of human cells is initiated by the binding of the viral Spike protein to its cell-surface receptor ACE2. We conducted a targeted CRISPRi screen to uncover druggable pathways controlling Spike protein binding to human cells. Here we show that the protein BRD2 is required for ACE2 transcription in human lung epithelial cells and cardiomyocytes, and BRD2 inhibitors currently evaluated in clinical trials potently block endogenous ACE2 expression and SARS-CoV-2 infection of human cells, including those of human nasal epithelia. Moreover, pharmacological BRD2 inhibition with the drug ABBV-744 inhibited SARS-CoV-2 replication in Syrian hamsters. We also found that BRD2 controls transcription of several other genes induced upon SARS-CoV-2 infection, including the interferon response, which in turn regulates the antiviral response. Together, our results pinpoint BRD2 as a potent and essential regulator of the host response to SARS-CoV-2 infection and highlight the potential of BRD2 as a therapeutic target for COVID-19.

High-content CRISPR screening.

CRISPR screens are a powerful source of biological discovery, enabling the unbiased interrogation of gene function in a wide range of applications and species. In pooled CRISPR screens, various genetically encoded perturbations are introduced into pools of cells. The targeted cells proliferate under a biological challenge such as cell competition, drug treatment or viral infection. Subsequently, the perturbation-induced effects are evaluated by sequencing-based counting of the guide RNAs that specify each perturbation. The typical results of such screens are ranked lists of genes that confer sensitivity or resistance to the biological challenge of interest. Contributing to the broad utility of CRISPR screens, adaptations of the core CRISPR technology make it possible to activate, silence or otherwise manipulate the target genes. Moreover, high-content read-outs such as single-cell RNA sequencing and spatial imaging help characterize screened cells with unprecedented detail. Dedicated software tools facilitate bioinformatic analysis and enhance reproducibility. CRISPR screening has unravelled various molecular mechanisms in basic biology, medical genetics, cancer research, immunology, infectious diseases, microbiology and other fields. This Primer describes the basic and advanced concepts of CRISPR screening and its application as a flexible and reliable method for biological discovery, biomedical research and drug development - with a special emphasis on high-content methods that make it possible to obtain detailed biological insights directly as part of the screen.