Massively parallel in vivo Perturb-seq reveals cell-type-specific transcriptional networks in cortical development.

Leveraging AAVs' versatile tropism and labeling capacity, we expanded the scale of in vivo CRISPR screening with single-cell transcriptomic phenotyping across embryonic to adult brains and peripheral nervous systems. Through extensive tests of 86 vectors across AAV serotypes combined with a transposon system, we substantially amplified labeling efficacy and accelerated in vivo gene delivery from weeks to days. Our proof-of-principle in utero screen identified the pleiotropic effects of Foxg1, highlighting its tight regulation of distinct networks essential for cell fate specification of Layer 6 corticothalamic neurons. Notably, our platform can label >6% of cerebral cells, surpassing the current state-of-the-art efficacy at <0.1% by lentivirus, to achieve analysis of over 30,000 cells in one experiment and enable massively parallel in vivo Perturb-seq. Compatible with various phenotypic measurements (single-cell or spatial multi-omics), it presents a flexible approach to interrogate gene function across cell types in vivo, translating gene variants to their causal function.

An evolved AAV variant enables efficient genetic engineering of murine T cells.

Precise targeting of large transgenes to T cells using homology-directed repair has been transformative for adoptive cell therapies and T cell biology. Delivery of DNA templates via adeno-associated virus (AAV) has greatly improved knockin efficiencies, but the tropism of current AAV serotypes restricts their use to human T cells employed in immunodeficient mouse models. To enable targeted knockins in murine T cells, we evolved Ark313, a synthetic AAV that exhibits high transduction efficiency in murine T cells. We performed a genome-wide knockout screen and identified QA2 as an essential factor for Ark313 infection. We demonstrate that Ark313 can be used for nucleofection-free DNA delivery, CRISPR-Cas9-mediated knockouts, and targeted integration of large transgenes. Ark313 enables preclinical modeling of Trac-targeted CAR-T and transgenic TCR-T cells in immunocompetent models. Efficient gene targeting in murine T cells holds great potential for improved cell therapies and opens avenues in experimental T cell immunology.

Mitigation of chromosome loss in clinical CRISPR-Cas9-engineered T cells.

CRISPR-Cas9 genome editing has enabled advanced T cell therapies, but occasional loss of the targeted chromosome remains a safety concern. To investigate whether Cas9-induced chromosome loss is a universal phenomenon and evaluate its clinical significance, we conducted a systematic analysis in primary human T cells. Arrayed and pooled CRISPR screens revealed that chromosome loss was generalizable across the genome and resulted in partial and entire loss of the targeted chromosome, including in preclinical chimeric antigen receptor T cells. T cells with chromosome loss persisted for weeks in culture, implying the potential to interfere with clinical use. A modified cell manufacturing process, employed in our first-in-human clinical trial of Cas9-engineered T cells (NCT03399448), reduced chromosome loss while largely preserving genome editing efficacy. Expression of p53 correlated with protection from chromosome loss observed in this protocol, suggesting both a mechanism and strategy for T cell engineering that mitigates this genotoxicity in the clinic.

Large-scale multiplexed mosaic CRISPR perturbation in the whole organism.

Here, we report inducible mosaic animal for perturbation (iMAP), a transgenic platform enabling in situ CRISPR targeting of at least 100 genes in parallel throughout the mouse body. iMAP combines Cre-loxP and CRISPR-Cas9 technologies and utilizes a germline-transmitted transgene carrying a large array of individually floxed, tandemly linked gRNA-coding units. Cre-mediated recombination triggers expression of all the gRNAs in the array but only one of them per cell, converting the mice to mosaic organisms suitable for phenotypic characterization and also for high-throughput derivation of conventional single-gene perturbation lines via breeding. Using gRNA representation as a readout, we mapped a miniature Perturb-Atlas cataloging the perturbations of 90 genes across 39 tissues, which yields rich insights into context-dependent gene functions and provides a glimpse of the potential of iMAP in genome decoding.

Emerging enterococcus pore-forming toxins with MHC/HLA-I as receptors.

Enterococci are a part of human microbiota and a leading cause of multidrug resistant infections. Here, we identify a family of Enterococcus pore-forming toxins (Epxs) in E. faecalis, E. faecium, and E. hirae strains isolated across the globe. Structural studies reveal that Epxs form a branch of ß-barrel pore-forming toxins with a ß-barrel protrusion (designated the top domain) sitting atop the cap domain. Through a genome-wide CRISPR-Cas9 screen, we identify human leukocyte antigen class I (HLA-I) complex as a receptor for two members (Epx2 and Epx3), which preferentially recognize human HLA-I and homologous MHC-I of equine, bovine, and porcine, but not murine, origin. Interferon exposure, which stimulates MHC-I expression, sensitizes human cells and intestinal organoids to Epx2 and Epx3 toxicity. Co-culture with Epx2-harboring E. faecium damages human peripheral blood mononuclear cells and intestinal organoids, and this toxicity is neutralized by an Epx2 antibody, demonstrating the toxin-mediated virulence of Epx-carrying Enterococcus.

Lrp1 is a host entry factor for Rift Valley fever virus.

Rift Valley fever virus (RVFV) is a zoonotic pathogen with pandemic potential. RVFV entry is mediated by the viral glycoprotein (Gn), but host entry factors remain poorly defined. Our genome-wide CRISPR screen identified low-density lipoprotein receptor-related protein 1 (mouse Lrp1/human LRP1), heat shock protein (Grp94), and receptor-associated protein (RAP) as critical host factors for RVFV infection. RVFV Gn directly binds to specific Lrp1 clusters and is glycosylation independent. Exogenous addition of murine RAP domain 3 (mRAPD3) and anti-Lrp1 antibodies neutralizes RVFV infection in taxonomically diverse cell lines. Mice treated with mRAPD3 and infected with pathogenic RVFV are protected from disease and death. A mutant mRAPD3 that binds Lrp1 weakly failed to protect from RVFV infection. Together, these data support Lrp1 as a host entry factor for RVFV infection and define a new target to limit RVFV infections.

Genome-wide CRISPR Screens Reveal Host Factors Critical for SARS-CoV-2 Infection.

Identification of host genes essential for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection may reveal novel therapeutic targets and inform our understanding of coronavirus disease 2019 (COVID-19) pathogenesis. Here we performed genome-wide CRISPR screens in Vero-E6 cells with SARS-CoV-2, Middle East respiratory syndrome CoV (MERS-CoV), bat CoV HKU5 expressing the SARS-CoV-1 spike, and vesicular stomatitis virus (VSV) expressing the SARS-CoV-2 spike. We identified known SARS-CoV-2 host factors, including the receptor ACE2 and protease Cathepsin L. We additionally discovered pro-viral genes and pathways, including HMGB1 and the SWI/SNF chromatin remodeling complex, that are SARS lineage and pan-coronavirus specific, respectively. We show that HMGB1 regulates ACE2 expression and is critical for entry of SARS-CoV-2, SARS-CoV-1, and NL63. We also show that small-molecule antagonists of identified gene products inhibited SARS-CoV-2 infection in monkey and human cells, demonstrating the conserved role of these genetic hits across species. This identifies potential therapeutic targets for SARS-CoV-2 and reveals SARS lineage-specific and pan-CoV host factors that regulate susceptibility to highly pathogenic CoVs.

Multimodal Analysis of Composition and Spatial Architecture in Human Squamous Cell Carcinoma.

To define the cellular composition and architecture of cutaneous squamous cell carcinoma (cSCC), we combined single-cell RNA sequencing with spatial transcriptomics and multiplexed ion beam imaging from a series of human cSCCs and matched normal skin. cSCC exhibited four tumor subpopulations, three recapitulating normal epidermal states, and a tumor-specific keratinocyte (TSK) population unique to cancer, which localized to a fibrovascular niche. Integration of single-cell and spatial data mapped ligand-receptor networks to specific cell types, revealing TSK cells as a hub for intercellular communication. Multiple features of potential immunosuppression were observed, including T regulatory cell (Treg) co-localization with CD8 T cells in compartmentalized tumor stroma. Finally, single-cell characterization of human tumor xenografts and in vivo CRISPR screens identified essential roles for specific tumor subpopulation-enriched gene networks in tumorigenesis. These data define cSCC tumor and stromal cell subpopulations, the spatial niches where they interact, and the communicating gene networks that they engage in cancer.

Systematic Immunotherapy Target Discovery Using Genome-Scale In Vivo CRISPR Screens in CD8 T Cells.

CD8 T cells play essential roles in anti-tumor immune responses. Here, we performed genome-scale CRISPR screens in CD8 T cells directly under cancer immunotherapy settings and identified regulators of tumor infiltration and degranulation. The in vivo screen robustly re-identified canonical immunotherapy targets such as PD-1 and Tim-3, along with genes that have not been characterized in T cells. The infiltration and degranulation screens converged on an RNA helicase Dhx37. Dhx37 knockout enhanced the efficacy of antigen-specific CD8 T cells against triple-negative breast cancer in vivo. Immunological characterization in mouse and human CD8 T cells revealed that DHX37 suppresses effector functions, cytokine production, and T cell activation. Transcriptomic profiling and biochemical interrogation revealed a role for DHX37 in modulating NF-κB. These data demonstrate high-throughput in vivo genetic screens for immunotherapy target discovery and establishes DHX37 as a functional regulator of CD8 T cells.

CCAR1 promotes DNA repair via alternative splicing.

DNA repair is directly performed by hundreds of core factors and indirectly regulated by thousands of others. We massively expanded a CRISPR inhibition and Cas9-editing screening system to discover factors indirectly modulating homology-directed repair (HDR) in the context of ∼18,000 individual gene knockdowns. We focused on CCAR1, a poorly understood gene that we found the depletion of reduced both HDR and interstrand crosslink repair, phenocopying the loss of the Fanconi anemia pathway. CCAR1 loss abrogated FANCA protein without substantial reduction in the level of its mRNA or that of other FA genes. We instead found that CCAR1 prevents inclusion of a poison exon in FANCA. Transcriptomic analysis revealed that the CCAR1 splicing modulatory activity is not limited to FANCA, and it instead regulates widespread changes in alternative splicing that would damage coding sequences in mouse and human cells. CCAR1 therefore has an unanticipated function as a splicing fidelity factor.

Genome-scale exon perturbation screens uncover exons critical for cell fitness.

CRISPR-Cas technology has transformed functional genomics, yet understanding of how individual exons differentially shape cellular phenotypes remains limited. Here, we optimized and conducted massively parallel exon deletion and splice-site mutation screens in human cell lines to identify exons that regulate cellular fitness. Fitness-promoting exons are prevalent in essential and highly expressed genes and commonly overlap with protein domains and interaction interfaces. Conversely, fitness-suppressing exons are enriched in nonessential genes, exhibiting lower inclusion levels, and overlap with intrinsically disordered regions and disease-associated mutations. In-depth mechanistic investigation of the screen-hit TAF5 alternative exon-8 revealed that its inclusion is required for assembly of the TFIID general transcription initiation complex, thereby regulating global gene expression output. Collectively, our orthogonal exon perturbation screens established a comprehensive repository of phenotypically important exons and uncovered regulatory mechanisms governing cellular fitness and gene expression.

MTHFD2 is a metabolic checkpoint controlling effector and regulatory T cell fate and function.

Antigenic stimulation promotes T cell metabolic reprogramming to meet increased biosynthetic, bioenergetic, and signaling demands. We show that the one-carbon (1C) metabolism enzyme methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) regulates de novo purine synthesis and signaling in activated T cells to promote proliferation and inflammatory cytokine production. In pathogenic T helper-17 (Th17) cells, MTHFD2 prevented aberrant upregulation of the transcription factor FoxP3 along with inappropriate gain of suppressive capacity. MTHFD2 deficiency also promoted regulatory T (Treg) cell differentiation. Mechanistically, MTHFD2 inhibition led to depletion of purine pools, accumulation of purine biosynthetic intermediates, and decreased nutrient sensor mTORC1 signaling. MTHFD2 was also critical to regulate DNA and histone methylation in Th17 cells. Importantly, MTHFD2 deficiency reduced disease severity in multiple in vivo inflammatory disease models. MTHFD2 is thus a metabolic checkpoint to integrate purine metabolism with pathogenic effector cell signaling and is a potential therapeutic target within 1C metabolism pathways.

A Genome-wide CRISPR Screen in Toxoplasma Identifies Essential Apicomplexan Genes.

Apicomplexan parasites are leading causes of human and livestock diseases such as malaria and toxoplasmosis, yet most of their genes remain uncharacterized. Here, we present the first genome-wide genetic screen of an apicomplexan. We adapted CRISPR/Cas9 to assess the contribution of each gene from the parasite Toxoplasma gondii during infection of human fibroblasts. Our analysis defines ∼200 previously uncharacterized, fitness-conferring genes unique to the phylum, from which 16 were investigated, revealing essential functions during infection of human cells. Secondary screens identify as an invasion factor the claudin-like apicomplexan microneme protein (CLAMP), which resembles mammalian tight-junction proteins and localizes to secretory organelles, making it critical to the initiation of infection. CLAMP is present throughout sequenced apicomplexan genomes and is essential during the asexual stages of the malaria parasite Plasmodium falciparum. These results provide broad-based functional information on T. gondii genes and will facilitate future approaches to expand the horizon of antiparasitic interventions.

FACS-based genome-wide CRISPR screens define key regulators of DNA damage signaling pathways.

DNA damage-activated signaling pathways are critical for coordinating multiple cellular processes, which must be tightly regulated to maintain genome stability. To provide a comprehensive and unbiased perspective of DNA damage response (DDR) signaling pathways, we performed 30 fluorescence-activated cell sorting (FACS)-based genome-wide CRISPR screens in human cell lines with antibodies recognizing distinct endogenous DNA damage signaling proteins to identify critical regulators involved in DDR. We discovered that proteasome-mediated processing is an early and prerequisite event for cells to trigger camptothecin- and etoposide-induced DDR signaling. Furthermore, we identified PRMT1 and PRMT5 as modulators that regulate ATM protein level. Moreover, we discovered that GNB1L is a key regulator of DDR signaling via its role as a co-chaperone specifically regulating PIKK proteins. Collectively, these screens offer a rich resource for further investigation of DDR, which may provide insight into strategies of targeting these DDR pathways to improve therapeutic outcomes.

Unbiased interrogation of functional lysine residues in human proteome.

CRISPR screens have empowered the high-throughput dissection of gene functions; however, more explicit genetic elements, such as codons of amino acids, require thorough interrogation. Here, we establish a CRISPR strategy for unbiasedly probing functional amino acid residues at the genome scale. By coupling adenine base editors and barcoded sgRNAs, we target 215,689 out of 611,267 (35%) lysine codons, involving 85% of the total protein-coding genes. We identify 1,572 lysine codons whose mutations perturb human cell fitness, with many of them implicated in cancer. These codons are then mirrored to gene knockout screen data to provide functional insights into the role of lysine residues in cellular fitness. Mining these data, we uncover a CUL3-centric regulatory network in which lysine residues of CUL3 CRL complex proteins control cell fitness by specifying protein-protein interactions. Our study offers a general strategy for interrogating genetic elements and provides functional insights into the human proteome.

ZNF410 Uniquely Activates the NuRD Component CHD4 to Silence Fetal Hemoglobin Expression.

Metazoan transcription factors typically regulate large numbers of genes. Here we identify via a CRISPR-Cas9 genetic screen ZNF410, a pentadactyl DNA-binding protein that in human erythroid cells directly activates only a single gene, the NuRD component CHD4. Specificity is conveyed by two highly evolutionarily conserved clusters of ZNF410 binding sites near the CHD4 gene with no counterparts elsewhere in the genome. Loss of ZNF410 in adult-type human erythroid cell culture systems and xenotransplantation settings diminishes CHD4 levels and derepresses the fetal hemoglobin genes. While previously known to be silenced by CHD4, the fetal globin genes are exposed here as among the most sensitive to reduced CHD4 levels.. In vitro DNA binding assays and crystallographic studies reveal the ZNF410-DNA binding mode. ZNF410 is a remarkably selective transcriptional activator in erythroid cells, and its perturbation might offer new opportunities for treatment of hemoglobinopathies.

Zmat3 Is a Key Splicing Regulator in the p53 Tumor Suppression Program.

Although TP53 is the most commonly mutated gene in human cancers, the p53-dependent transcriptional programs mediating tumor suppression remain incompletely understood. Here, to uncover critical components downstream of p53 in tumor suppression, we perform unbiased RNAi and CRISPR-Cas9-based genetic screens in vivo. These screens converge upon the p53-inducible gene Zmat3, encoding an RNA-binding protein, and we demonstrate that ZMAT3 is an important tumor suppressor downstream of p53 in mouse KrasG12D-driven lung and liver cancers and human carcinomas. Integrative analysis of the ZMAT3 RNA-binding landscape and transcriptomic profiling reveals that ZMAT3 directly modulates exon inclusion in transcripts encoding proteins of diverse functions, including the p53 inhibitors MDM4 and MDM2, splicing regulators, and components of varied cellular processes. Interestingly, these exons are enriched in NMD signals, and, accordingly, ZMAT3 broadly affects target transcript stability. Collectively, these studies reveal ZMAT3 as a novel RNA-splicing and homeostasis regulator and a key component of p53-mediated tumor suppression.

Synthetic Lethal and Resistance Interactions with BET Bromodomain Inhibitors in Triple-Negative Breast Cancer.

BET bromodomain inhibitors (BBDIs) are candidate therapeutic agents for triple-negative breast cancer (TNBC) and other cancer types, but inherent and acquired resistance to BBDIs limits their potential clinical use. Using CRISPR and small-molecule inhibitor screens combined with comprehensive molecular profiling of BBDI response and resistance, we identified synthetic lethal interactions with BBDIs and genes that, when deleted, confer resistance. We observed synergy with regulators of cell cycle progression, YAP, AXL, and SRC signaling, and chemotherapeutic agents. We also uncovered functional similarities and differences among BRD2, BRD4, and BRD7. Although deletion of BRD2 enhances sensitivity to BBDIs, BRD7 loss leads to gain of TEAD-YAP chromatin binding and luminal features associated with BBDI resistance. Single-cell RNA-seq, ATAC-seq, and cellular barcoding analysis of BBDI responses in sensitive and resistant cell lines highlight significant heterogeneity among samples and demonstrate that BBDI resistance can be pre-existing or acquired.

HNRNPH1 regulates the neuroprotective cold-shock protein RBM3 expression through poison exon exclusion.

Enhanced expression of the cold-shock protein RNA binding motif 3 (RBM3) is highly neuroprotective both in vitro and in vivo. Whilst upstream signalling pathways leading to RBM3 expression have been described, the precise molecular mechanism of RBM3 cold induction remains elusive. To identify temperature-dependent modulators of RBM3, we performed a genome-wide CRISPR-Cas9 knockout screen using RBM3-reporter human iPSC-derived neurons. We found that RBM3 mRNA and protein levels are robustly regulated by several splicing factors, with heterogeneous nuclear ribonucleoprotein H1 (HNRNPH1) being the strongest positive regulator. Splicing analysis revealed that moderate hypothermia significantly represses the inclusion of a poison exon, which, when retained, targets the mRNA for nonsense-mediated decay. Importantly, we show that HNRNPH1 mediates this cold-dependent exon skipping via its thermosensitive interaction with a G-rich motif within the poison exon. Our study provides novel mechanistic insights into the regulation of RBM3 and provides further targets for neuroprotective therapeutic strategies.

A large-scale CRISPR screen reveals context-specific genetic regulation of retinal ganglion cell regeneration.

Many genes are known to regulate retinal regeneration after widespread tissue damage. Conversely, genes controlling regeneration after limited cell loss, as per degenerative diseases, are undefined. As stem/progenitor cell responses scale to injury levels, understanding how the extent and specificity of cell loss impact regenerative processes is important. Here, transgenic zebrafish enabling selective retinal ganglion cell (RGC) ablation were used to identify genes that regulate RGC regeneration. A single cell multiomics-informed screen of 100 genes identified seven knockouts that inhibited and 11 that promoted RGC regeneration. Surprisingly, 35 out of 36 genes known and/or implicated as being required for regeneration after widespread retinal damage were not required for RGC regeneration. The loss of seven even enhanced regeneration kinetics, including the proneural factors neurog1, olig2 and ascl1a. Mechanistic analyses revealed that ascl1a disruption increased the propensity of progenitor cells to produce RGCs, i.e. increased 'fate bias'. These data demonstrate plasticity in the mechanism through which Müller glia convert to a stem-like state and context specificity in how genes function during regeneration. Increased understanding of how the regeneration of disease-relevant cell types is specifically controlled will support the development of disease-tailored regenerative therapeutics.