Optimized libraries for CRISPR-Cas9 genetic screens with multiple modalities

Allele-specific CRISPR-Cas9 editing inactivates a single nucleotide variant associated with collagen VI muscular dystrophy

The application of allele-specific gene editing tools can expand the therapeutic options for dominant genetic conditions, either via gene correction or via allelic gene inactivation in situations where haploinsufficiency is tolerated. Here, we used allele-targeted CRISPR-Cas9 guide RNAs (gRNAs) to introduce inactivating frameshifting indels at an SNV in the COL6A1 gene (c.868G>A; G290R), a variant that acts as dominant negative and that is associated with a severe form of congenital muscular dystrophy. We expressed SpCas9 along with allele-targeted gRNAs, without providing a repair template, in primary fibroblasts derived from four patients and one control subject. Amplicon deep sequencing for two gRNAs tested showed that single-nucleotide deletions accounted for the majority of indels introduced. While activity of the two gRNAs was greater at the G290R allele, both gRNAs were also active at the wild-type allele. To enhance allele selectivity, we introduced deliberate additional mismatches to one gRNA. One of these optimized gRNAs showed minimal activity at the WT allele, while generating productive edits and improving collagen VI matrix in cultured patient fibroblasts. This study strengthens the potential of gene editing to treat dominant-negative disorders, but also underscores the challenges in achieving allele selectivity with gRNAs.

Overexpression of human SAMD9 inhibits protein translation and alters MYC signaling resulting in cell cycle arrest

Inherited bone marrow failure syndromes often result from pathogenic mutations in genes that are important for ribosome function, namely, Diamond-Blackfan anemia, Shwachman-Diamond syndrome, and dyskeratosis congenita. Germline mutations in SAMD9 are a frequent genetic lesion resulting in an inherited bone marrow failure syndrome with monosomy 7; some patients have severe multisystem syndromes that include myelodysplasia. The association of germline SAMD9 mutations and bone marrow failure is clear; however, to date, there is no reliable method to predict whether a novel SAMD9 mutation is pathogenic unless it is accompanied by an obvious family history and/or clinical syndrome. The difficulty with pathogenicity prediction is, in part, due to the incomplete understanding of the biological functions of SAMD9. We used a SAMD9-targeted, inducible CRISPRa system and RNA sequencing to better understand the global transcriptional changes that result from transcriptional manipulation of SAMD9. Supporting recent discoveries that SAMD9 acts as a ACNase specific for phenylalanine tRNA (tRNA-Phe), we confirmed with crosslinking and solid-phase purification that SAMD9 is an RNA binding protein and analyzed how overexpression of tRNA-Phe may reverse transcriptomic changes caused by SAMD9 activation. Our data show that overexpression of SAMD9 from the endogenous locus results in decreased cell proliferation, cell cycle progression, and global protein translation. When SAMD9 contains a gain-of-function mutation (p.E1136Q), these functional phenotypes are exacerbated but only partially rescued with tRNA-Phe overexpression, suggesting additional molecular actions of SAMD9. Additionally, we demonstrate that gene expression pathways important for ribosome biogenesis and MYC signaling are the most significantly impacted by SAMD9 overexpression.

Altered RNA export by SF3B1 mutants confers sensitivity to nuclear export inhibition

SF3B1 mutations frequently occur in cancer yet lack targeted therapies. Clinical trials of XPO1 inhibitors, selinexor and eltanexor, in high-risk myelodysplastic neoplasms (MDS) revealed responders were enriched with SF3B1 mutations. Given that XPO1 (Exportin-1) is a nuclear exporter responsible for the export of proteins and multiple RNA species, this led to the hypothesis that SF3B1-mutant cells are sensitive to XPO1 inhibition, potentially due to altered splicing. Subsequent RNA sequencing after XPO1 inhibition in SF3B1 wildtype and mutant cells showed increased nuclear retention of RNA transcripts and increased alternative splicing in the SF3B1 mutant cells particularly of genes that impact apoptotic pathways. To identify novel drug combinations that synergize with XPO1 inhibition, a forward genetic screen was performed with eltanexor treatment implicating anti-apoptotic targets BCL2 and BCLXL, which were validated by functional testing in vitro and in vivo. These targets were tested in vivo using Sf3b1K700E conditional knock-in mice, which showed that the combination of eltanexor and venetoclax (BCL2 inhibitor) had a preferential sensitivity for SF3B1 mutant cells without excessive toxicity. In this study, we unveil the mechanisms underlying sensitization to XPO1 inhibition in SF3B1-mutant MDS and preclinically rationalize the combination of eltanexor and venetoclax for high-risk MDS.

Integrating genome-wide CRISPR screens and in silico drug profiling for targeted antidote development

Numerous toxins threaten humans, but specific antidotes are unavailable for most of them. Although CRISPR screening has aided the discovery of the mechanisms of some toxins, developing targeted antidotes remains a significant challenge. Recently, we established a systematic framework to develop antidotes by combining the identification of novel drug targets by using a genome-wide CRISPR screen with a virtual screen of drugs approved by the US Food and Drug Administration. This approach allows for a comprehensive understanding of toxin mechanisms at the whole-genome level and facilitates the identification of promising antidote drugs targeting specific molecules. Here, we present step-by-step instructions for executing genome-scale CRISPR-Cas9 knockout screens of toxins in HAP1 cells. We also provide detailed guidance for conducting an in silico drug screen and an in vivo drug validation. By using this protocol, it takes ~4 weeks to perform the genome-scale screen, 4 weeks for sequencing and data analysis, 4 weeks to validate candidate genes, 1 week for the virtual screen and 2 weeks for in vitro drug validation. This framework has the potential to accelerate the development of antidotes for a wide range of toxins and can rapidly identify promising drug candidates that are already known to be safe and effective. This could lead to the development of new antidotes much more quickly than traditional methods, protecting lives from diverse toxins and advancing human health.

Pharmacogenomic discovery of genetically targeted cancer therapies optimized against clinical outcomes

Despite the clinical success of dozens of genetically targeted cancer therapies, the vast majority of patients with tumors caused by loss-of-function (LoF) mutations do not have access to these treatments. This is primarily due to the challenge of developing a drug that treats a disease caused by the absence of a protein target. The success of PARP inhibitors has solidified synthetic lethality (SL) as a means to overcome this obstacle. Recent mapping of SL networks using pooled CRISPR-Cas9 screens is a promising approach for expanding this concept to treating cancers driven by additional LoF drivers. In practice, however, translating signals from cell lines, where these screens are typically conducted, to patient outcomes remains a challenge. We developed a pharmacogenomic (PGx) approach called "Clinically Optimized Driver Associated-PGx" (CODA-PGX) that accurately predicts genetically targeted therapies with clinical-stage efficacy in specific LoF driver contexts. Using approved targeted therapies and cancer drugs with available real-world evidence and molecular data from hundreds of patients, we discovered and optimized the key screening principles predictive of efficacy and overall patient survival. In addition to establishing basic technical conventions, such as drug concentration and screening kinetics, we found that replicating the driver perturbation in the right context, as well as selecting patients where those drivers are genuine founder mutations, were key to accurate translation. We used CODA-PGX to screen a diverse collection of clinical stage drugs and report dozens of novel LoF genetically targeted opportunities; many validated in xenografts and by real-world evidence. Notable examples include treating STAG2-mutant tumors with Carboplatin, SMARCB1-mutant tumors with Oxaliplatin, and TP53BP1-mutant tumors with Etoposide or Bleomycin.

CRISPR screening uncovers a long-range enhancer for ONECUT1 in pancreatic differentiation and links a diabetes risk variant

Functional enhancer annotation is critical for understanding tissue-specific transcriptional regulation and prioritizing disease-associated non-coding variants. However, unbiased enhancer discovery in disease-relevant contexts remains challenging. To identify enhancers pertinent to diabetes, we conducted a CRISPR interference (CRISPRi) screen in the human pluripotent stem cell (hPSC) pancreatic differentiation system. Among the enhancers identified, we focused on an enhancer we named ONECUT1e-664kb, ∼664 kb from the ONECUT1 promoter. Previous studies have linked ONECUT1 coding mutations to pancreatic hypoplasia and neonatal diabetes. We found that homozygous deletion of ONECUT1e-664kb in hPSCs leads to a near-complete loss of ONECUT1 expression and impaired pancreatic differentiation. ONECUT1e-664kb contains a type 2 diabetes-associated variant (rs528350911) disrupting a GATA motif. Introducing the risk variant into hPSCs reduced binding of key pancreatic transcription factors (GATA4, GATA6, and FOXA2), supporting its causal role in diabetes. This work highlights the utility of unbiased enhancer discovery in disease-relevant settings for understanding monogenic and complex disease.

Zika virus NS5 protein inhibits type I interferon signaling via CRL3 E3 ubiquitin ligase-mediated degradation of STAT2

The ZIKA virus (ZIKV) evades the host immune response by degrading STAT2 through its NS5 protein, thereby inhibiting type I interferon (IFN)-mediated antiviral immunity. However, the molecular mechanism underlying this process has remained elusive. In this study, we performed a genome-wide CRISPR/Cas9 screen, revealing that ZSWIM8 as the substrate receptor of Cullin3-RING E3 ligase is required for NS5-mediated STAT2 degradation. Genetic depletion of ZSWIM8 and CUL3 substantially impeded NS5-mediated STAT2 degradation. Biochemical analysis illuminated that NS5 enhances the interaction between STAT2 and the ZSWIM8-CUL3 E3 ligase complex, thereby facilitating STAT2 ubiquitination. Moreover, ZSWIM8 knockout endowed A549 and Huh7 cells with partial resistance to ZIKV infection and protected cells from the cytopathic effects induced by ZIKV, which was attributed to the restoration of STAT2 levels and the activation of IFN signaling. Subsequent studies in a physiologically relevant model, utilizing human neural progenitor cells, demonstrated that ZSWIM8 depletion reduced ZIKV infection, resulting from enhanced IFN signaling attributed to the sustained levels of STAT2. Our findings shed light on the role of ZIKV NS5, serving as the scaffold protein, reprograms the ZSWIM8-CUL3 E3 ligase complex to orchestrate STAT2 proteasome-dependent degradation, thereby facilitating evasion of IFN antiviral signaling. Our study provides unique insights into ZIKV-host interactions and holds promise for the development of antivirals and prophylactic vaccines.

Productive mRNA Chromatin Escape is Promoted by PRMT5 Methylation of SNRPB

Protein Arginine Methyltransferase 5 (PRMT5) regulates RNA splicing and transcription by symmetric dimethylation of arginine residues (Rme2s/SDMA) in many RNA binding proteins. However, the mechanism by which PRMT5 couples splicing to transcriptional output is unknown. Here, we demonstrate that a major function of PRMT5 activity is to promote chromatin escape of a novel, large class of mRNAs that we term Genomically Retained Incompletely Processed Polyadenylated Transcripts (GRIPPs). Using nascent and total transcriptomics, spike-in controlled fractionated cell transcriptomics, and total and fractionated cell proteomics, we show that PRMT5 inhibition and knockdown of the PRMT5 SNRP (Sm protein) adapter protein pICln (CLNS1A) -but not type I PRMT inhibition-leads to gross detention of mRNA, SNRPB, and SNRPD3 proteins on chromatin. Compared to most transcripts, these chromatin-trapped polyadenylated RNA transcripts have more introns, are spliced slower, and are enriched in detained introns. Using a combination of PRMT5 inhibition and inducible isogenic wildtype and arginine-mutant SNRPB, we show that arginine methylation of these snRNPs is critical for mediating their homeostatic chromatin and RNA interactions. Overall, we conclude that a major role for PRMT5 is in controlling transcript processing and splicing completion to promote chromatin escape and subsequent nuclear export.

Mammalian SWI/SNF complex activity regulates POU2F3 and constitutes a targetable dependency in small cell lung cancer

Small cell lung cancers (SCLCs) are composed of heterogeneous subtypes marked by lineage-specific transcription factors, including ASCL1, NEUROD1, and POU2F3. POU2F3-positive SCLCs, ∼12% of all cases, are uniquely dependent on POU2F3 itself; as such, approaches to attenuate POU2F3 expression may represent new therapeutic opportunities. Here using genome-scale screens for regulators of POU2F3 expression and SCLC proliferation, we define mSWI/SNF complexes as top dependencies specific to POU2F3-positive SCLC. Notably, chemical disruption of mSWI/SNF ATPase activity attenuates proliferation of all POU2F3-positive SCLCs, while disruption of non-canonical BAF (ncBAF) via BRD9 degradation is effective in pure non-neuroendocrine POU2F3-SCLCs. mSWI/SNF targets to and maintains accessibility over gene loci central to POU2F3-mediated gene regulatory networks. Finally, clinical-grade pharmacologic disruption of SMARCA4/2 ATPases and BRD9 decreases POU2F3-SCLC tumor growth and increases survival in vivo. These results demonstrate mSWI/SNF-mediated governance of the POU2F3 oncogenic program and suggest mSWI/SNF inhibition as a therapeutic strategy for POU2F3-positive SCLCs.

Translation efficiency covariation across cell types is a conserved organizing principle of mammalian transcriptomes

Characterization of shared patterns of RNA expression between genes across conditions has led to the discovery of regulatory networks and novel biological functions. However, it is unclear if such coordination extends to translation, a critical step in gene expression. Here, we uniformly analyzed 3,819 ribosome profiling datasets from 117 human and 94 mouse tissues and cell lines. We introduce the concept of Translation Efficiency Covariation (TEC), identifying coordinated translation patterns across cell types. We nominate potential mechanisms driving shared patterns of translation regulation. TEC is conserved across human and mouse cells and helps uncover gene functions. Moreover, our observations indicate that proteins that physically interact are highly enriched for positive covariation at both translational and transcriptional levels. Our findings establish translational covariation as a conserved organizing principle of mammalian transcriptomes.

Proximity analysis of native proteomes reveals phenotypic modifiers in a mouse model of autism and related neurodevelopmental conditions

One of the main drivers of autism spectrum disorder is risk alleles within hundreds of genes, which may interact within shared but unknown protein complexes. Here we develop a scalable genome-editing-mediated approach to target 14 high-confidence autism risk genes within the mouse brain for proximity-based endogenous proteomics, achieving the identification of high-specificity spatial proteomes. The resulting native proximity proteomes are enriched for human genes dysregulated in the brain of autistic individuals, and reveal proximity interactions between proteins from high-confidence risk genes with those of lower-confidence that may provide new avenues to prioritize genetic risk. Importantly, the datasets are enriched for shared cellular functions and genetic interactions that may underlie the condition. We test this notion by spatial proteomics and CRISPR-based regulation of expression in two autism models, demonstrating functional interactions that modulate mechanisms of their dysregulation. Together, these results reveal native proteome networks in vivo relevant to autism, providing new inroads for understanding and manipulating the cellular drivers underpinning its etiology.

RNA variant assessment using transactivation and transdifferentiation

Understanding the impact of splicing and nonsense variants on RNA is crucial for the resolution of variant classification as well as their suitability for precision medicine interventions. This is primarily enabled through RNA studies involving transcriptomics followed by targeted assays using RNA isolated from clinically accessible tissues (CATs) such as blood or skin of affected individuals. Insufficient disease gene expression in CATs does however pose a major barrier to RNA based investigations, which we show is relevant to 1,436 Mendelian disease genes. We term these "silent" Mendelian genes (SMGs), the largest portion (36%) of which are associated with neurological disorders. We developed two approaches to induce SMG expression in human dermal fibroblasts (HDFs) to overcome this limitation, including CRISPR-activation-based gene transactivation and fibroblast-to-neuron transdifferentiation. Initial transactivation screens involving 40 SMGs stimulated our development of a highly multiplexed transactivation system culminating in the 6- to 90,000-fold induction of expression of 20/20 (100%) SMGs tested in HDFs. Transdifferentiation of HDFs directly to neurons led to expression of 193/516 (37.4%) of SMGs implicated in neurological disease. The magnitude and isoform diversity of SMG expression following either transactivation or transdifferentiation was comparable to clinically relevant tissues. We apply transdifferentiation and/or gene transactivation combined with short- and long-read RNA sequencing to investigate the impact that variants in USH2A, SCN1A, DMD, and PAK3 have on RNA using HDFs derived from affected individuals. Transactivation and transdifferentiation represent rapid, scalable functional genomic solutions to investigate variants impacting SMGs in the patient cell and genomic context.

Genetic therapies and potential therapeutic applications of CRISPR activators in the eye

Conventional gene therapy involving supplementation only treats loss-of-function diseases and is limited by viral packaging sizes, precluding therapy of large genes. The discovery of CRISPR/Cas has led to a paradigm shift in the field of genetic therapy, with the promise of precise gene editing, thus broadening the range of diseases that can be treated. The initial uses of CRISPR/Cas have focused mainly on gene editing or silencing of abnormal variants via utilising Cas endonuclease to trigger the target cell endogenous non-homologous end joining. Subsequently, the technology has evolved to modify the Cas enzyme and even its guide RNA, leading to more efficient editing tools in the form of base and prime editing. Further advancements of this CRISPR/Cas technology itself have expanded its functional repertoire from targeted editing to programmable transactivation, shifting the therapeutic focus to precise endogenous gene activation or upregulation with the potential for epigenetic modifications. In vivo experiments using this platform have demonstrated the potential of CRISPR-activators (CRISPRa) to treat various loss-of-function diseases, as well as in regenerative medicine, highlighting their versatility to overcome limitations associated with conventional strategies. This review summarises the molecular mechanisms of CRISPRa platforms, the current applications of this technology in vivo, and discusses potential solutions to translational hurdles for this therapy, with a focus on ophthalmic diseases.

The extracellular matrix integrates mitochondrial homeostasis

Cellular homeostasis is intricately influenced by stimuli from the microenvironment, including signaling molecules, metabolites, and pathogens. Functioning as a signaling hub within the cell, mitochondria integrate information from various intracellular compartments to regulate cellular signaling and metabolism. Multiple studies have shown that mitochondria may respond to various extracellular signaling events. However, it is less clear how changes in the extracellular matrix (ECM) can impact mitochondrial homeostasis to regulate animal physiology. We find that ECM remodeling alters mitochondrial homeostasis in an evolutionarily conserved manner. Mechanistically, ECM remodeling triggers a TGF-β response to induce mitochondrial fission and the unfolded protein response of the mitochondria (UPRMT). At the organismal level, ECM remodeling promotes defense of animals against pathogens through enhanced mitochondrial stress responses. We postulate that this ECM-mitochondria crosstalk represents an ancient immune pathway, which detects infection- or mechanical-stress-induced ECM damage, thereby initiating adaptive mitochondria-based immune and metabolic responses.

Multiplexed single-cell characterization of alternative polyadenylation regulators

Most mammalian genes have multiple polyA sites, representing a substantial source of transcript diversity regulated by the cleavage and polyadenylation (CPA) machinery. To better understand how these proteins govern polyA site choice, we introduce CPA-Perturb-seq, a multiplexed perturbation screen dataset of 42 CPA regulators with a 3' scRNA-seq readout that enables transcriptome-wide inference of polyA site usage. We develop a framework to detect perturbation-dependent changes in polyadenylation and characterize modules of co-regulated polyA sites. We find groups of intronic polyA sites regulated by distinct components of the nuclear RNA life cycle, including elongation, splicing, termination, and surveillance. We train and validate a deep neural network (APARENT-Perturb) for tandem polyA site usage, delineating a cis-regulatory code that predicts perturbation response and reveals interactions between regulatory complexes. Our work highlights the potential for multiplexed single-cell perturbation screens to further our understanding of post-transcriptional regulation.

The neonatal Fc receptor (FcRn) is a pan-arterivirus receptor

Arteriviruses infect a variety of mammalian hosts, but the receptors used by these viruses to enter cells are poorly understood. We identified the neonatal Fc receptor (FcRn) as an important pro-viral host factor via comparative genome-wide CRISPR-knockout screens with multiple arteriviruses. Using a panel of cell lines and divergent arteriviruses, we demonstrate that FcRn is required for the entry step of arterivirus infection and serves as a molecular barrier to arterivirus cross-species infection. We also show that FcRn synergizes with another known arterivirus entry factor, CD163, to mediate arterivirus entry. Overexpression of FcRn and CD163 sensitizes non-permissive cells to infection and enables the culture of fastidious arteriviruses. Treatment of multiple cell lines with a pre-clinical anti-FcRn monoclonal antibody blocked infection and rescued cells from arterivirus-induced death. Altogether, this study identifies FcRn as a novel pan-arterivirus receptor, with implications for arterivirus emergence, cross-species infection, and host-directed pan-arterivirus countermeasure development.

Non-linear transcriptional responses to gradual modulation of transcription factor dosage

Genomic loci associated with common traits and diseases are typically non-coding and likely impact gene expression, sometimes coinciding with rare loss-of-function variants in the target gene. However, our understanding of how gradual changes in gene dosage affect molecular, cellular, and organismal traits is currently limited. To address this gap, we induced gradual changes in gene expression of four genes using CRISPR activation and inactivation. Downstream transcriptional consequences of dosage modulation of three master trans-regulators associated with blood cell traits (GFI1B, NFE2, and MYB) were examined using targeted single-cell multimodal sequencing. We showed that guide tiling around the TSS is the most effective way to modulate cis gene expression across a wide range of fold-changes, with further effects from chromatin accessibility and histone marks that differ between the inhibition and activation systems. Our single-cell data allowed us to precisely detect subtle to large gene expression changes in dozens of trans genes, revealing that many responses to dosage changes of these three TFs are non-linear, including non-monotonic behaviours, even when constraining the fold-changes of the master regulators to a copy number gain or loss. We found that the dosage properties are linked to gene constraint and that some of these non-linear responses are enriched for disease and GWAS genes. Overall, our study provides a straightforward and scalable method to precisely modulate gene expression and gain insights into its downstream consequences at high resolution.

Comprehensive transcription factor perturbations recapitulate fibroblast transcriptional states

Cell atlas projects have nominated recurrent transcriptional states as drivers of biological processes and disease, but their origins, regulation, and properties remain unclear. To enable complementary functional studies, we developed a scalable approach for recapitulating cell states in vitro using CRISPR activation (CRISPRa) Perturb-seq. Aided by a novel multiplexing method, we activated 1,836 transcription factors in two cell types. Measuring 21,958 perturbations showed that CRISPRa activated targets within physiological ranges, that epigenetic features predicted activatable genes, and that the protospacer seed region drove an off-target effect. Perturbations recapitulated in vivo fibroblast states, including universal and inflammatory states, and identified KLF4 and KLF5 as key regulators of the universal state. Inducing the universal state suppressed disease-associated states, highlighting its therapeutic potential. Our findings cement CRISPRa as a tool for perturbing differentiated cells and indicate that in vivo states can be elicited via perturbation, enabling studies of clinically relevant states ex vivo.

Epigenome editing technologies for discovery and medicine

Epigenome editing has rapidly evolved in recent years, with diverse applications that include elucidating gene regulation mechanisms, annotating coding and noncoding genome functions and programming cell state and lineage specification. Importantly, given the ubiquitous role of epigenetics in complex phenotypes, epigenome editing has unique potential to impact a broad spectrum of diseases. By leveraging powerful DNA-targeting technologies, such as CRISPR, epigenome editing exploits the heritable and reversible mechanisms of epigenetics to alter gene expression without introducing DNA breaks, inducing DNA damage or relying on DNA repair pathways.

Phosphorylation of GCN2 by mTOR confers adaptation to conditions of hyper-mTOR activation under stress

Adaptation to the shortage in free amino acids (AA) is mediated by 2 pathways, the integrated stress response (ISR) and the mechanistic target of rapamycin (mTOR). In response to reduced levels, primarily of leucine or arginine, mTOR in its complex 1 configuration (mTORC1) is suppressed leading to a decrease in translation initiation and elongation. The eIF2α kinase general control nonderepressible 2 (GCN2) is activated by uncharged tRNAs, leading to induction of the ISR in response to a broader range of AA shortage. ISR confers a reduced translation initiation, while promoting the selective synthesis of stress proteins, such as ATF4. To efficiently adapt to AA starvation, the 2 pathways are cross-regulated at multiple levels. Here we identified a new mechanism of ISR/mTORC1 crosstalk that optimizes survival under AA starvation, when mTORC1 is forced to remain active. mTORC1 activation during acute AA shortage, augmented ATF4 expression in a GCN2-dependent manner. Under these conditions, enhanced GCN2 activity was not dependent on tRNA sensing, inferring a different activation mechanism. We identified a labile physical interaction between GCN2 and mTOR that results in a phosphorylation of GCN2 on serine 230 by mTOR, which promotes GCN2 activity. When examined under prolonged AA starvation, GCN2 phosphorylation by mTOR promoted survival. Our data unveils an adaptive mechanism to AA starvation, when mTORC1 evades inhibition.

The proteogenomic landscape of multiple myeloma reveals insights into disease biology and therapeutic opportunities

Multiple myeloma (MM) is a plasma cell malignancy of the bone marrow. Despite therapeutic advances, MM remains incurable, and better risk stratification as well as new therapies are therefore highly needed. The proteome of MM has not been systematically assessed before and holds the potential to uncover insight into disease biology and improved prognostication in addition to genetic and transcriptomic studies. Here we provide a comprehensive multiomics analysis including deep tandem mass tag-based quantitative global (phospho)proteomics, RNA sequencing, and nanopore DNA sequencing of 138 primary patient-derived plasma cell malignancies encompassing treatment-naive MM, plasma cell leukemia and the premalignancy monoclonal gammopathy of undetermined significance, as well as healthy controls. We found that the (phospho)proteome of malignant plasma cells are highly deregulated as compared with healthy plasma cells and is both defined by chromosomal alterations as well as posttranscriptional regulation. A prognostic protein signature was identified that is associated with aggressive disease independent of established risk factors in MM. Integration with functional genetics and single-cell RNA sequencing revealed general and genetic subtype-specific deregulated proteins and pathways in plasma cell malignancies that include potential targets for (immuno)therapies. Our study demonstrates the potential of proteogenomics in cancer and provides an easily accessible resource for investigating protein regulation and new therapeutic approaches in MM.

Dissecting key regulators of transcriptome kinetics through scalable single-cell RNA profiling of pooled CRISPR screens

We present a combinatorial indexing method, PerturbSci-Kinetics, for capturing whole transcriptomes, nascent transcriptomes and single guide RNA (sgRNA) identities across hundreds of genetic perturbations at the single-cell level. Profiling a pooled CRISPR screen targeting various biological processes, we show the gene expression regulation during RNA synthesis, processing and degradation, miRNA biogenesis and mitochondrial mRNA processing, systematically decoding the genome-wide regulatory network that underlies RNA temporal dynamics at scale.

In vivo interaction screening reveals liver-derived constraints to metastasis

It is estimated that only 0.02% of disseminated tumour cells are able to seed overt metastases1. While this suggests the presence of environmental constraints to metastatic seeding, the landscape of host factors controlling this process remains largely unclear. Here, combining transposon technology2 and fluorescence niche labelling3, we developed an in vivo CRISPR activation screen to systematically investigate the interactions between hepatocytes and metastatic cells. We identify plexin B2 as a critical host-derived regulator of liver colonization in colorectal and pancreatic cancer and melanoma syngeneic mouse models. We dissect a mechanism through which plexin B2 interacts with class IV semaphorins on tumour cells, leading to KLF4 upregulation and thereby promoting the acquisition of epithelial traits. Our results highlight the essential role of signals from the liver parenchyma for the seeding of disseminated tumour cells before the establishment of a growth-promoting niche. Our findings further suggest that epithelialization is required for the adaptation of CRC metastases to their new tissue environment. Blocking the plexin-B2-semaphorin axis abolishes metastatic colonization of the liver and therefore represents a therapeutic strategy for the prevention of hepatic metastases. Finally, our screening approach, which evaluates host-derived extrinsic signals rather than tumour-intrinsic factors for their ability to promote metastatic seeding, is broadly applicable and lays a framework for the screening of environmental constraints to metastasis in other organs and cancer types.

CRISPR activation screens: navigating technologies and applications

Clustered regularly interspaced short palindromic repeats (CRISPR) activation (CRISPRa) has become an integral part of the molecular biology toolkit. CRISPRa genetic screens are an exciting high-throughput means of identifying genes the upregulation of which is sufficient to elicit a given phenotype. Activation machinery is continually under development to achieve greater, more robust, and more consistent activation. In this review, we offer a succinct technological overview of available CRISPRa architectures and a comprehensive summary of pooled CRISPRa screens. Furthermore, we discuss contemporary applications of CRISPRa across broad fields of research, with the aim of presenting a view of exciting emerging applications for CRISPRa screening.

The histone methyltransferase KMT2D maintains cellular glucocorticoid responsiveness by shielding the glucocorticoid receptor from degradation

Because of their ability to induce lymphocyte apoptosis, glucocorticoids (GC) are widely used to treat hematological malignancies such as lymphomas and multiple myeloma. Their effectiveness is often limited, however, due to the development of glucocorticoid resistance by a variety of molecular mechanisms. Here we performed an unbiased genome-wide CRISPR screen with the human T-cell leukemia cell line Jurkat to find previously unidentified genes required for GC-induced apoptosis. One such gene was KMT2D (also known as MLL2 or MLL4), which encodes a histone lysine methyltransferase whose mutations are associated with a variety of cancers, blood malignancies in particular, and are considered markers of poor prognosis. Knockout of KMT2D by CRISPR/Cas9 gene editing in Jurkat and several multiple myeloma cell lines downregulated GR protein expression. Surprisingly, this was not due to a reduction in GR transcripts, but rather to a decrease in the protein's half-life, primarily due to proteasomal degradation. Reconstitution of KMT2D expression restored GR levels. In contrast to the known ability of KMT2D to control gene transcription through covalent histone methylation, KMT2D-mediated upregulation of GR levels did not require its methyltransferase activity. Co-immunoprecipitation and proximity ligation assays found constitutive binding of KMT2D to the GR, which was enhanced in the presence of GC. These observations reveal KMT2D to be essential for the stabilization of cellular GR levels, and suggest a possible mechanism by which KMT2D mutations may lead to GC resistance in some malignancies.

The integrated stress response regulates 18S nonfunctional rRNA decay in mammals

18S nonfunctional rRNA decay (NRD) detects and eliminates translationally nonfunctional 18S rRNA. While this process is critical for ribosome quality control, the mechanisms underlying nonfunctional 18S rRNA turnover remain elusive. NRD was originally identified and has exclusively been studied in Saccharomyces cerevisiae. Here, we show that 18S NRD is conserved in mammals. Using genome-wide CRISPR genetic interaction screens, we find that mammalian NRD acts through the integrated stress response (ISR) via GCN2 and ribosomal protein ubiquitination by RNF10. Selective ribosome profiling reveals nonfunctional 18S rRNA induces translational arrest at start sites. Indeed, biochemical analyses demonstrate that ISR activation limits translation initiation and attenuates collisions between scanning 43S preinitiation complexes and nonfunctional 80S ribosomes arrested at start sites. Overall, the ISR promotes nonfunctional 18S rRNA and 40S ribosomal protein turnover by RNF10-mediated ubiquitination. These findings establish a dynamic feedback mechanism by which the GCN2-RNF10 axis surveils ribosome functionality at translation initiation.

Interface-guided phenotyping of coding variants in the transcription factor RUNX1

Single-gene missense mutations remain challenging to interpret. Here, we deploy scalable functional screening by sequencing (SEUSS), a Perturb-seq method, to generate mutations at protein interfaces of RUNX1 and quantify their effect on activities of downstream cellular programs. We evaluate single-cell RNA profiles of 115 mutations in myelogenous leukemia cells and categorize them into three functionally distinct groups, wild-type (WT)-like, loss-of-function (LoF)-like, and hypomorphic, that we validate in orthogonal assays. LoF-like variants dominate the DNA-binding site and are recurrent in cancer; however, recurrence alone does not predict functional impact. Hypomorphic variants share characteristics with LoF-like but favor protein interactions, promoting gene expression indicative of nerve growth factor (NGF) response and cytokine recruitment of neutrophils. Accessible DNA near differentially expressed genes frequently contains RUNX1-binding motifs. Finally, we reclassify 16 variants of uncertain significance and train a classifier to predict 103 more. Our work demonstrates the potential of targeting protein interactions to better define the landscape of phenotypes reachable by missense mutations.

NSD2 drives t(4;14) myeloma cell dependence on adenylate kinase 2 by diverting one-carbon metabolism to the epigenome

ABSTRACT

Chromosomal translocation (4;14), an adverse prognostic factor in multiple myeloma (MM), drives overexpression of the histone methyltransferase nuclear receptor binding SET domain protein 2 (NSD2). A genome-wide CRISPR screen in MM cells identified adenylate kinase 2 (AK2), an enzyme critical for high-energy phosphate transfer from the mitochondria, as an NSD2-driven vulnerability. AK2 suppression in t(4;14) MM cells decreased nicotinamide adenine dinucleotide phosphate (NADP[H]) critical for conversion of ribonucleotides to deoxyribonucleosides, leading to replication stress, DNA damage, and apoptosis. Driving a large genome-wide increase in chromatin methylation, NSD2 overexpression depletes S-adenosylmethionine, compromising the synthesis of creatine from its precursor, guanidinoacetate. Creatine supplementation restored NADP(H) levels, reduced DNA damage, and rescued AK2-deficient t(4;14) MM cells. As the creatine phosphate shuttle constitutes an alternative means for mitochondrial high-energy phosphate transport, these results indicate that NSD2-driven creatine depletion underlies the hypersensitivity of t(4;14) MM cells to AK2 loss. Furthermore, AK2 depletion in t(4;14) cells impaired protein folding in the endoplasmic reticulum, consistent with impaired use of mitochondrial adenosine triphosphate (ATP). Accordingly, AK2 suppression increased the sensitivity of MM cells to proteasome inhibition. These findings delineate a novel mechanism in which aberrant transfer of carbon to the epigenome creates a metabolic vulnerability, with direct therapeutic implications for t(4;14) MM.

MethNet: a robust approach to identify regulatory hubs and their distal targets from cancer data

Aberrations in the capacity of DNA/chromatin modifiers and transcription factors to bind non-coding regions can lead to changes in gene regulation and impact disease phenotypes. However, identifying distal regulatory elements and connecting them with their target genes remains challenging. Here, we present MethNet, a pipeline that integrates large-scale DNA methylation and gene expression data across multiple cancers, to uncover cis regulatory elements (CREs) in a 1 Mb region around every promoter in the genome. MethNet identifies clusters of highly ranked CREs, referred to as 'hubs', which contribute to the regulation of multiple genes and significantly affect patient survival. Promoter-capture Hi-C confirmed that highly ranked associations involve physical interactions between CREs and their gene targets, and CRISPR interference based single-cell RNA Perturb-seq validated the functional impact of CREs. Thus, MethNet-identified CREs represent a valuable resource for unraveling complex mechanisms underlying gene expression, and for prioritizing the verification of predicted non-coding disease hotspots.

RBM10 loss induces aberrant splicing of cytoskeletal and extracellular matrix mRNAs and promotes metastatic fitness

RBM10 modulates transcriptome-wide cassette exon splicing. Loss-of-function RBM10 mutations are enriched in thyroid cancers with distant metastases. Analysis of transcriptomes and genes mis-spliced by RBM10 loss showed pro-migratory and RHO/RAC signaling signatures. RBM10 loss increases cell velocity. Cytoskeletal and ECM transcripts subject to exon-inclusion events included vinculin (VCL), tenascin C (TNC) and CD44. Knockdown of the VCL exon inclusion transcript in RBM10-null cells reduced cell velocity, whereas knockdown of TNC and CD44 exon-inclusion isoforms reduced invasiveness. RAC1-GTP levels were increased in RBM10-null cells. Mouse Hras G12V /Rbm1O KO thyrocytes develop metastases that are reversed by RBM10 or by combined knockdown of VCL, CD44 and TNC inclusion isoforms. Thus, RBM10 loss generates exon inclusions in transcripts regulating ECM-cytoskeletal interactions, leading to RAC1 activation and metastatic competency. Moreover, a CRISPR-Cas9 screen for synthetic lethality with RBM10 loss identified NFkB effectors as central to viability, providing a therapeutic target for these lethal thyroid cancers.

A whole-genome CRISPR screen identifies the spindle accessory checkpoint as a locus of nab-paclitaxel resistance in a pancreatic cancer cell line

Pancreatic adenocarcinoma is one of the most aggressive and lethal forms of cancer. Chemotherapy is the primary treatment for pancreatic cancer, but resistance to the drugs used remains a major challenge. A genome-wide CRISPR interference and knockout screen in the PANC-1 cell line with the drug nab-paclitaxel has identified a group of spindle assembly checkpoint (SAC) genes that enhance survival in nab-paclitaxel. Knockdown of these SAC genes (BUB1B, BUB3, and TTK) attenuates paclitaxel-induced cell death. Cells treated with the small molecule inhibitors BAY 1217389 or MPI 0479605, targeting the threonine tyrosine kinase (TTK), also enhance survival in paclitaxel. Overexpression of these SAC genes does not affect sensitivity to paclitaxel. These discoveries have helped to elucidate the mechanisms behind paclitaxel cytotoxicity. The outcomes of this investigation may pave the way for a deeper comprehension of the diverse responses of pancreatic cancer to therapies including paclitaxel. Additionally, they could facilitate the formulation of novel treatment approaches for pancreatic cancer.

Epistatic interactions between NMD and TRP53 control progenitor cell maintenance and brain size

Mutations in human nonsense-mediated mRNA decay (NMD) factors are enriched in neurodevelopmental disorders. We show that deletion of key NMD factor Upf2 in mouse embryonic neural progenitor cells causes perinatal microcephaly but deletion in immature neurons does not, indicating NMD's critical roles in progenitors. Upf2 knockout (KO) prolongs the cell cycle of radial glia progenitor cells, promotes their transition into intermediate progenitors, and leads to reduced upper-layer neurons. CRISPRi screening identified Trp53 knockdown rescuing Upf2KO progenitors without globally reversing NMD inhibition, implying marginal contributions of most NMD targets to the cell cycle defect. Integrated functional genomics shows that NMD degrades selective TRP53 downstream targets, including Cdkn1a, which, without NMD suppression, slow the cell cycle. Trp53KO restores the progenitor cell pool and rescues the microcephaly of Upf2KO mice. Therefore, one physiological role of NMD in the developing brain is to degrade selective TRP53 targets to control progenitor cell cycle and brain size.

CRISPR-Cas9 Knockout Screens Identify DNA Damage Response Pathways and BTK as Essential for Cisplatin Response in Diffuse Large B-Cell Lymphoma

The recurrence of diffuse large B-cell lymphoma (DLBCL) has been observed in 40% of cases. The standard of care for refractory/relapsed DLBCL (RR-DLBCL) is platinum-based treatment prior to autologous stem cell transplantation; however, the prognosis for RR-DLBCL patients remains poor. Thus, to identify genes affecting the cisplatin response in DLBCL, cisplatin-based whole-genome CRISPR-Cas9 knockout screens were performed in this study. We discovered DNA damage response (DDR) pathways as enriched among identified sensitizing CRISPR-mediated gene knockouts. In line, the knockout of the nucleotide excision repair genes XPA and ERCC6 sensitized DLBCL cells to platinum drugs irrespective of proliferation rate, thus documenting DDR as essential for cisplatin sensitivity in DLBCL. Functional analysis revealed that the loss of XPA and ERCC6 increased DNA damage levels and altered cell cycle distribution. Interestingly, we also identified BTK, which is involved in B-cell receptor signaling, to affect cisplatin response. The knockout of BTK increased cisplatin sensitivity in DLBCL cells, and combinatory drug screens revealed a synergistic effect of the BTK inhibitor, ibrutinib, with platinum drugs at low concentrations. Applying local and external DLBCL cohorts, we addressed the clinical relevance of the genes identified in the CRISPR screens. BTK was among the most frequently mutated genes with a frequency of 3-5%, and XPA and ERCC6 were also mutated, albeit at lower frequencies. Furthermore, 27-54% of diagnostic DLBCL samples had mutations in pathways that can sensitize cells to cisplatin. In conclusion, this study shows that XPA and ERCC6, in addition to BTK, are essential for the response to platinum-based drugs in DLBCL.

Splice variants of CK1α and CK1α-like: Comparative analysis of subcellular localization, kinase activity, and function in the Wnt signaling pathway

Members of the casein kinase 1 (CK1) family are important regulators of multiple signaling pathways. CK1α is a well-known negative regulator of the Wnt/β-catenin pathway, which promotes the degradation of β-catenin via its phosphorylation of Ser45. In contrast, the closest paralog of CK1α, CK1α-like, is a poorly characterized kinase of unknown function. In this study, we show that the deletion of CK1α, but not CK1α-like, resulted in a strong activation of the Wnt/β-catenin pathway. Wnt-3a treatment further enhanced the activation, which suggests there are at least two modes, a CK1α-dependent and Wnt-dependent, of β-catenin regulation. Rescue experiments showed that only two out of ten naturally occurring splice CK1α/α-like variants were able to rescue the augmented Wnt/β-catenin signaling caused by CK1α deficiency in cells. Importantly, the ability to phosphorylate β-catenin on Ser45 in the in vitro kinase assay was required but not sufficient for such rescue. Our compound CK1α and GSK3α/β KO models suggest that the additional nonredundant function of CK1α in the Wnt pathway beyond Ser45-β-catenin phosphorylation includes Axin phosphorylation. Finally, we established NanoBRET assays for the three most common CK1α splice variants as well as CK1α-like. Target engagement data revealed comparable potency of known CK1α inhibitors for all CK1α variants but not for CK1α-like. In summary, our work brings important novel insights into the biology of CK1α, including evidence for the lack of redundancy with other CK1 kinases in the negative regulation of the Wnt/β-catenin pathway at the level of β-catenin and Axin.

CRISPR Screening Uncovers a Long-Range Enhancer for ONECUT1 in Pancreatic Differentiation and Links a Diabetes Risk Variant

Functional enhancer annotation is a valuable first step for understanding tissue-specific transcriptional regulation and prioritizing disease-associated non-coding variants for investigation. However, unbiased enhancer discovery in physiologically relevant contexts remains a major challenge. To discover regulatory elements pertinent to diabetes, we conducted a CRISPR interference screen in the human pluripotent stem cell (hPSC) pancreatic differentiation system. Among the enhancers uncovered, we focused on a long-range enhancer ∼664 kb from the ONECUT1 promoter, since coding mutations in ONECUT1 cause pancreatic hypoplasia and neonatal diabetes. Homozygous enhancer deletion in hPSCs was associated with a near-complete loss of ONECUT1 gene expression and compromised pancreatic differentiation. This enhancer contains a confidently fine-mapped type 2 diabetes associated variant (rs528350911) which disrupts a GATA motif. Introduction of the risk variant into hPSCs revealed substantially reduced binding of key pancreatic transcription factors (GATA4, GATA6 and FOXA2) on the edited allele, accompanied by a slight reduction of ONECUT1 transcription, supporting a causal role for this risk variant in metabolic disease. This work expands our knowledge about transcriptional regulation in pancreatic development through the characterization of a long-range enhancer and highlights the utility of enhancer discovery in disease-relevant settings for understanding monogenic and complex disease.

CRISPR activation screens identify the SWI/SNF ATPases as suppressors of ferroptosis

Ferroptosis is an iron-dependent cell death mechanism characterized by the accumulation of toxic lipid peroxides and cell membrane rupture. GPX4 (glutathione peroxidase 4) prevents ferroptosis by reducing these lipid peroxides into lipid alcohols. Ferroptosis induction by GPX4 inhibition has emerged as a vulnerability of cancer cells, highlighting the need to identify ferroptosis regulators that may be exploited therapeutically. Through genome-wide CRISPR activation screens, we identify the SWI/SNF (switch/sucrose non-fermentable) ATPases BRM (SMARCA2) and BRG1 (SMARCA4) as ferroptosis suppressors. Mechanistically, they bind to and increase chromatin accessibility at NRF2 target loci, thus boosting NRF2 transcriptional output to counter lipid peroxidation and confer resistance to GPX4 inhibition. We further demonstrate that the BRM/BRG1 ferroptosis connection can be leveraged to enhance the paralog dependency of BRG1 mutant cancer cells on BRM. Our data reveal ferroptosis induction as a potential avenue for broadening the efficacy of BRM degraders/inhibitors and define a specific genetic context for exploiting GPX4 dependency.

Epigenomic signatures of sarcomatoid differentiation to guide the treatment of renal cell carcinoma

Renal cell carcinoma with sarcomatoid differentiation (sRCC) is associated with poor survival and a heightened response to immune checkpoint inhibitors (ICIs). Two major barriers to improving outcomes for sRCC are the limited understanding of its gene regulatory programs and the low diagnostic yield of tumor biopsies due to spatial heterogeneity. Herein, we characterized the epigenomic landscape of sRCC by profiling 107 epigenomic libraries from tissue and plasma samples from 50 patients with RCC and healthy volunteers. By profiling histone modifications and DNA methylation, we identified highly recurrent epigenomic reprogramming enriched in sRCC. Furthermore, CRISPRa experiments implicated the transcription factor FOSL1 in activating sRCC-associated gene regulatory programs, and FOSL1 expression was associated with the response to ICIs in RCC in two randomized clinical trials. Finally, we established a blood-based diagnostic approach using detectable sRCC epigenomic signatures in patient plasma, providing a framework for discovering epigenomic correlates of tumor histology via liquid biopsy.

Variants in tubule epithelial regulatory elements mediate most heritable differences in human kidney function

Kidney disease is highly heritable; however, the causal genetic variants, the cell types in which these variants function, and the molecular mechanisms underlying kidney disease remain largely unknown. To identify genetic loci affecting kidney function, we performed a GWAS using multiple kidney function biomarkers and identified 462 loci. To begin to investigate how these loci affect kidney function, we generated single-cell chromatin accessibility (scATAC-seq) maps of the human kidney and identified candidate cis-regulatory elements (cCREs) for kidney podocytes, tubule epithelial cells, and kidney endothelial, stromal, and immune cells. Kidney tubule epithelial cCREs explained 58% of kidney function SNP-heritability and kidney podocyte cCREs explained an additional 6.5% of SNP-heritability. In contrast, little kidney function heritability was explained by kidney endothelial, stromal, or immune cell-specific cCREs. Through functionally informed fine-mapping, we identified putative causal kidney function variants and their corresponding cCREs. Using kidney scATAC-seq data, we created a deep learning model (which we named ChromKid) to predict kidney cell type-specific chromatin accessibility from sequence. ChromKid and allele specific kidney scATAC-seq revealed that many fine-mapped kidney function variants locally change chromatin accessibility in tubule epithelial cells. Enhancer assays confirmed that fine-mapped kidney function variants alter tubule epithelial regulatory element function. To map the genes which these regulatory elements control, we used CRISPR interference (CRISPRi) to target these regulatory elements in tubule epithelial cells and assessed changes in gene expression. CRISPRi of enhancers harboring kidney function variants regulated NDRG1 and RBPMS expression. Thus, inherited differences in tubule epithelial NDRG1 and RBPMS expression may predispose to kidney disease in humans. We conclude that genetic variants affecting tubule epithelial regulatory element function account for most SNP-heritability of human kidney function. This work provides an experimental approach to identify the variants, regulatory elements, and genes involved in polygenic disease.

Thiophenyl Derivatives of Nicotinamide Are Metabolized by the NAD Salvage Pathway into Unnatural NAD Derivatives That Inhibit IMPDH and Are Toxic to Peripheral Nerve Cancers

N-Pyridinylthiophene carboxamide (compound 21) displays activity against peripheral nerve sheath cancer cells and mouse xenografts by an unknown mechanism. Through medicinal chemistry, we identified a more active derivative, compound 9, and found that only analogues with structures similar to nicotinamide retained activity. Genetic screens using compound 9 found that both NAMPT and NMNAT1, enzymes in the NAD salvage pathway, are necessary for activity. Compound 9 is metabolized by NAMPT and NMNAT1 into an adenine dinucleotide (AD) derivative in a cell-free system, cultured cells, and mice, and inhibition of this metabolism blocked compound activity. AD analogues derived from compound 9 inhibit IMPDH in vitro and cause cell death by inhibiting IMPDH in cells. These findings nominate these compounds as preclinical candidates for the development of tumor-activated IMPDH inhibitors to treat neuronal cancers.

Defining cis-regulatory elements and transcription factors that control human cortical interneuron development

Although human cortical interneurons (cINs) are a minority population in the cerebral cortex, disruption of interneuron development is a frequent contributor to neurodevelopmental disorders. Here, we utilized a model for deriving cINs from human embryonic stem cells to profile chromatin state changes and generate an atlas of cis-regulatory elements (CREs) controlling human cIN development. We used these data to define candidate transcription factors (TFs) that may bind these CREs to regulate interneuron progenitor specification. Among these were RFX3 and RFX4, risk genes for autism spectrum disorder (ASD) with uncharacterized roles in human neuronal development. Using RFX3 and RFX4 knockdown models, we demonstrated new requirements for both genes in interneuron progenitor specification, with RFX3 deficiency causing precocious neuronal differentiation while RFX4 deficiency instead resulted in cessation of progenitor cell proliferation. Together, this work both defined central features of cis-regulatory control and identified new TF requirements for human interneuron development.

Genome-Wide CRISPR Screen Identifies an NF2-Adherens Junction Mechanistic Dependency for Cardiac Lineage

BACKGROUND

Cardiomyocyte differentiation involves a stepwise clearance of repressors and fate-restricting regulators through the modulation of BMP (bone morphogenic protein)/Wnt-signaling pathways. However, the mechanisms and how regulatory roadblocks are removed with specific developmental signaling pathways remain unclear.

METHODS

We conducted a genome-wide CRISPR screen to uncover essential regulators of cardiomyocyte specification in human embryonic stem cells using a myosin heavy chain 6 (MYH6)-GFP (green fluorescence protein) reporter system. After an independent secondary single guide ribonucleic acid validation of 25 candidates, we identified NF2 (neurofibromin 2), a moesin-ezrin-radixin like (MERLIN) tumor suppressor, as an upstream driver of early cardiomyocyte lineage specification. Independent monoclonal NF2 knockouts were generated using CRISPR-Cas9, and cell states were inferred through bulk RNA sequencing and protein expression analysis across differentiation time points. Terminal lineage differentiation was assessed by using an in vitro 2-dimensional-micropatterned gastruloid model, trilineage differentiation, and cardiomyocyte differentiation. Protein interaction and post-translation modification of NF2 with its interacting partners were assessed using site-directed mutagenesis, coimmunoprecipitation, and proximity ligation assays.

RESULTS

Transcriptional regulation and trajectory inference from NF2-null cells reveal the loss of cardiomyocyte identity and the acquisition of nonmesodermal identity. Sustained elevation of early mesoderm lineage repressor SOX2 and upregulation of late anticardiac regulators CDX2 and MSX1 in NF2 knockout cells reflect a necessary role for NF2 in removing regulatory roadblocks. Furthermore, we found that NF2 and AMOT (angiomotin) cooperatively bind to YAP (yes-associated protein) during mesendoderm formation, thereby preventing YAP activation, independent of canonical MST (mammalian sterile 20-like serine-threonine protein kinase)-LATS (large tumor suppressor serine-threonine protein kinase) signaling. Mechanistically, cardiomyocyte lineage identity was rescued by wild-type and NF2 serine-518 phosphomutants, but not NF2 FERM (ezrin-radixin-meosin homology protein) domain blue-box mutants, demonstrating that the critical FERM domain-dependent formation of the AMOT-NF2-YAP scaffold complex at the adherens junction is required for early cardiomyocyte lineage differentiation.

CONCLUSIONS

These results provide mechanistic insight into the essential role of NF2 during early epithelial-mesenchymal transition by sequestering the repressive effect of YAP and relieving regulatory roadblocks en route to cardiomyocytes.

CRISPR-Cas9 screens reveal common essential miRNAs in human cancer cell lines

BACKGROUND

Genome-wide functional screening using the CRISPR-Cas9 system is a powerful tool to uncover tumor-specific and common genetic dependencies across cancer cell lines. Current CRISPR-Cas9 knockout libraries, however, primarily target protein-coding genes. This limits functional genomics-based investigations of miRNA function.

METHODS

We designed a novel CRISPR-Cas9 knockout library (lentiG-miR) of 8107 distinct sgRNAs targeting a total of 1769 human miRNAs and benchmarked its single guide RNA (sgRNA) composition, predicted on- and off-target activity, and screening performance against previous libraries. Using a total of 45 human cancer cell lines, representing 16 different tumor entities, we performed negative selection screens to identify miRNA fitness genes. Fitness miRNAs in each cell line were scored using a combination of supervised and unsupervised essentiality classifiers. Common essential miRNAs across distinct cancer cell lines were determined using the 90th percentile method. For subsequent validation, we performed knockout experiments for selected common essential miRNAs in distinct cancer cell lines and gene expression profiling.

RESULTS

We found significantly lower off-target activity for protein-coding genes and a higher miRNA gene coverage for lentiG-miR as compared to previously described miRNA-targeting libraries, while preserving high on-target activity. A minor fraction of miRNAs displayed robust depletion of targeting sgRNAs, and we observed a high level of consistency between redundant sgRNAs targeting the same miRNA gene. Across 45 human cancer cell lines, only 217 (12%) of all targeted human miRNAs scored as a fitness gene in at least one model, and fitness effects for most miRNAs were confined to small subsets of cell lines. In contrast, we identified 49 common essential miRNAs with a homogenous fitness profile across the vast majority of all cell lines. Transcriptional profiling verified highly consistent gene expression changes in response to knockout of individual common essential miRNAs across a diverse set of cancer cell lines.

CONCLUSIONS

Our study presents a miRNA-targeting CRISPR-Cas9 knockout library with high gene coverage and optimized on- and off-target activities. Taking advantage of the lentiG-miR library, we define a catalogue of miRNA fitness genes in human cancer cell lines, providing the foundation for further investigation of miRNAs in human cancer.

NRP1 is a receptor for mammalian orthoreovirus engaged by distinct capsid subunits

Mammalian orthoreovirus (reovirus) is a nonenveloped virus that establishes primary infection in the intestine and disseminates to sites of secondary infection, including the CNS. Reovirus entry involves multiple engagement factors, but how the virus disseminates systemically and targets neurons remains unclear. In this study, we identified murine neuropilin 1 (mNRP1) as a receptor for reovirus. mNRP1 binds reovirus with nanomolar affinity using a unique mechanism of virus-receptor interaction, which is coordinated by multiple interactions between distinct reovirus capsid subunits and multiple NRP1 extracellular domains. By exchanging essential capsid protein-encoding gene segments, we determined that the multivalent interaction is mediated by outer-capsid protein σ3 and capsid turret protein λ2. Using capsid mutants incapable of binding NRP1, we found that NRP1 contributes to reovirus dissemination and neurovirulence in mice. Collectively, our results demonstrate that NRP1 is an entry receptor for reovirus and uncover mechanisms by which NRPs promote viral entry and pathogenesis.

Identification of genetic modifiers of Huntington's disease somatic CAG repeat instability by in vivo CRISPR-Cas9 genome editing

Huntington's disease (HD), one of >50 inherited repeat expansion disorders (Depienne and Mandel, 2021), is a dominantly-inherited neurodegenerative disease caused by a CAG expansion in HTT (The Huntington's Disease Collaborative Research Group, 1993). Inherited CAG repeat length is the primary determinant of age of onset, with human genetic studies underscoring that the property driving disease is the CAG length-dependent propensity of the repeat to further expand in brain (Swami et al ., 2009; GeM-HD, 2015; Hensman Moss et al ., 2017; Ciosi et al ., 2019; GeM-HD, 2019; Hong et al ., 2021). Routes to slowing somatic CAG expansion therefore hold great promise for disease-modifying therapies. Several DNA repair genes, notably in the mismatch repair (MMR) pathway, modify somatic expansion in HD mouse models (Wheeler and Dion, 2021). To identify novel modifiers of somatic expansion, we have used CRISPR-Cas9 editing in HD knock-in mice to enable in vivo screening of expansion-modifier candidates at scale. This has included testing of HD onset modifier genes emerging from human genome-wide association studies (GWAS), as well as interactions between modifier genes, thereby providing new insight into pathways underlying CAG expansion and potential therapeutic targets.

Germline variation contributes to false negatives in CRISPR-based experiments with varying burden across ancestries

Reducing disparities is vital for equitable access to precision treatments in cancer. Socioenvironmental factors are a major driver of disparities, but differences in genetic variation likely also contribute. The impact of genetic ancestry on prioritization of cancer targets in drug discovery pipelines has not been systematically explored due to the absence of pre-clinical data at the appropriate scale. Here, we analyze data from 611 genome-scale CRISPR/Cas9 viability experiments in human cell line models to identify ancestry-associated genetic dependencies essential for cell survival. Surprisingly, we find that most putative associations between ancestry and dependency arise from artifacts related to germline variants. Our analysis suggests that for 1.2-2.5% of guides, germline variants in sgRNA targeting sequences reduce cutting by the CRISPR/Cas9 nuclease, disproportionately affecting cell models derived from individuals of recent African descent. We propose three approaches to mitigate this experimental bias, enabling the scientific community to address these disparities.

Therapeutic validation of MMR-associated genetic modifiers in a human ex vivo model of Huntington disease

The pathological huntingtin (HTT) trinucleotide repeat underlying Huntington disease (HD) continues to expand throughout life. Repeat length correlates both with earlier age at onset (AaO) and faster progression, making slowing its expansion an attractive therapeutic approach. Genome-wide association studies have identified candidate variants associated with altered AaO and progression, with many found in DNA mismatch repair (MMR)-associated genes. We examine whether lowering expression of these genes affects the rate of repeat expansion in human ex vivo models using HD iPSCs and HD iPSC-derived striatal medium spiny neuron-enriched cultures. We have generated a stable CRISPR interference HD iPSC line in which we can specifically and efficiently lower gene expression from a donor carrying over 125 CAG repeats. Lowering expression of each member of the MMR complexes MutS (MSH2, MSH3, and MSH6), MutL (MLH1, PMS1, PMS2, and MLH3), and LIG1 resulted in characteristic MMR deficiencies. Reduced MSH2, MSH3, and MLH1 slowed repeat expansion to the largest degree, while lowering either PMS1, PMS2, or MLH3 slowed it to a lesser degree. These effects were recapitulated in iPSC-derived striatal cultures where MutL factor expression was lowered. CRISPRi-mediated lowering of key MMR factor expression to levels feasibly achievable by current therapeutic approaches was able to effectively slow the expansion of the HTT CAG tract. We highlight members of the MutL family as potential targets to slow pathogenic repeat expansion with the aim to delay onset and progression of HD and potentially other repeat expansion disorders exhibiting somatic instability.

TMPRSS13 promotes the cell entry of swine acute diarrhea syndrome coronavirus

Swine acute diarrhea syndrome coronavirus (SADS-CoV) has caused severe intestinal diseases in pigs. It originates from bat coronaviruses HKU2 and has a potential risk of cross-species transmission, raising concerns about its zoonotic potential. Viral entry-related host factors are critical determinants of susceptibility to cells, tissues, or species, and remain to be elucidated for SADS-CoV. Type II transmembrane serine proteases (TTSPs) family is involved in many coronavirus infections and has trypsin-like catalytic activity. Here we examine all 18 members of the TTSPs family through CRISPR-based activation of endogenous protein expression in cells, and find that, in addition to TMPRSS2 and TMPRSS4, TMPRSS13 significantly facilitates SADS-CoV infection. This is confirmed by ectopic expression of TMPRSS13, and specific to trypsin-dependent SADS-CoV. Infection with pseudovirus bearing SADS-CoV spike protein indicates that TMPRSS13 acts at the entry step and is sensitive to serine protease inhibitor Camostat. Moreover, both human and pig TMPRSS13 are able to enhance the cell-cell membrane fusion and cleavage of spike protein. Overall, we demonstrate that TMPRSS13 is another host serine protease promoting the membrane-fusion entry of SADS-CoV, which may expand its host tropism by using diverse TTSPs.

Integrated enhancer regulatory network by enhancer-promoter looping in gastric cancer

Enhancer cis-regulatory elements play critical roles in gene regulation at many stages of cell growth. Enhancers in cancer cells also regulate the transcription of oncogenes. In this study, we performed a comprehensive analysis of long-range chromatin interactions, histone modifications, chromatin accessibility and expression in two gastric cancer (GC) cell lines compared to normal gastric epithelial cells. We found that GC-specific enhancers marked by histone modifications can activate a population of genes, including some oncogenes, by interacting with their proximal promoters. In addition, motif analysis of enhancer-promoter interacting enhancers showed that GC-specific transcription factors are enriched. Among them, we found that MYB is crucial for GC cell growth and activated by the enhancer with an enhancer-promoter loop and TCF7 upregulation. Clinical GC samples showed epigenetic activation of enhancers at the MYB locus and significant upregulation of TCF7 and MYB, regardless of molecular GC subtype and clinicopathological factors. Single-cell RNA sequencing of gastric mucosa with intestinal metaplasia showed high expression of TCF7 and MYB in intestinal stem cells. When we inactivated the loop-forming enhancer at the MYB locus using CRISPR interference (dCas9-KRAB), GC cell growth was significantly inhibited. In conclusion, we identified MYB as an oncogene activated by a loop-forming enhancer and contributing to GC cell growth.

Targetable leukaemia dependency on noncanonical PI3Kγ signalling

Phosphoinositide-3-kinase-γ (PI3Kγ) is implicated as a target to repolarize tumour-associated macrophages and promote antitumour immune responses in solid cancers1-4. However, cancer cell-intrinsic roles of PI3Kγ are unclear. Here, by integrating unbiased genome-wide CRISPR interference screening with functional analyses across acute leukaemias, we define a selective dependency on the PI3Kγ complex in a high-risk subset that includes myeloid, lymphoid and dendritic lineages. This dependency is characterized by innate inflammatory signalling and activation of phosphoinositide 3-kinase regulatory subunit 5 (PIK3R5), which encodes a regulatory subunit of PI3Kγ5 and stabilizes the active enzymatic complex. We identify p21 (RAC1)-activated kinase 1 (PAK1) as a noncanonical substrate of PI3Kγ that mediates this cell-intrinsic dependency and find that dephosphorylation of PAK1 by PI3Kγ inhibition impairs mitochondrial oxidative phosphorylation. Treatment with the selective PI3Kγ inhibitor eganelisib is effective in leukaemias with activated PIK3R5. In addition, the combination of eganelisib and cytarabine prolongs survival over either agent alone, even in patient-derived leukaemia xenografts with low baseline PIK3R5 expression, as residual leukaemia cells after cytarabine treatment have elevated G protein-coupled purinergic receptor activity and PAK1 phosphorylation. Together, our study reveals a targetable dependency on PI3Kγ-PAK1 signalling that is amenable to near-term evaluation in patients with acute leukaemia.

A transcriptional enhancer regulates cardiac maturation

Cardiomyocyte maturation is crucial for generating adult cardiomyocytes and the application of human pluripotent stem cell-derived cardiomyocytes (hPSC-CMs). However, regulation at the cis-regulatory element level and its role in heart disease remain unclear. Alpha-actinin 2 (ACTN2) levels increase during CM maturation. In this study, we investigated a clinically relevant, conserved ACTN2 enhancer's effects on CM maturation using hPSC and mouse models. Heterozygous ACTN2 enhancer deletion led to abnormal CM morphology, reduced function and mitochondrial respiration. Transcriptomic analyses in vitro and in vivo showed disrupted CM maturation and upregulated anabolic mammalian target for rapamycin (mTOR) signaling, promoting senescence and hindering maturation. As confirmation, ACTN2 enhancer deletion induced heat shock protein 90A expression, a chaperone mediating mTOR activation. Conversely, targeting the ACTN2 enhancer via enhancer CRISPR activation (enCRISPRa) promoted hPSC-CM maturation. Our studies reveal the transcriptional enhancer's role in cardiac maturation and disease, offering insights into potentially fine-tuning gene expression to modulate cardiomyocyte physiology.