CRISPR/Cas9-Mediated Scanning for Regulatory Elements Required for HPRT1 Expression via Thousands of Large, Programmed Genomic Deletions

Engineering structural variants to interrogate genome function

Structural variation, such as deletions, duplications, inversions and complex rearrangements, can have profound effects on gene expression, genome stability, phenotypic diversity and disease susceptibility. Structural variants can encompass up to millions of bases and have the potential to rearrange substantial segments of the genome. They contribute considerably more to genetic diversity in human populations and have larger effects on phenotypic traits than point mutations. Until recently, our understanding of the effects of structural variants was driven mainly by studying naturally occurring variation. New genome-engineering tools capable of generating deletions, insertions, inversions and translocations, together with the discovery of new recombinases and advances in creating synthetic DNA constructs, now enable the design and generation of an extended range of structural variation. Here, we discuss these tools and examples of their application and highlight existing challenges that will need to be overcome to fully harness their potential.

Mapping cis- and trans-regulatory target genes of human-specific deletions

Deletion of functional sequence is predicted to represent a fundamental mechanism of molecular evolution1,2. Comparative genetic studies of primates2,3 have identified thousands of human-specific deletions (hDels), and the cis-regulatory potential of short (≤31 base pairs) hDels has been assessed using reporter assays4. However, how structural variant-sized (≥50 base pairs) hDels influence molecular and cellular processes in their native genomic contexts remains unexplored. Here, we design genome-scale libraries of single-guide RNAs targeting 7.2 megabases of sequence in 6,358 hDels and present a systematic CRISPR interference (CRISPRi) screening approach to identify hDels that modify cellular proliferation in chimpanzee pluripotent stem cells. By intersecting hDels with chromatin state features and performing single-cell CRISPRi (Perturb-seq) to identify their cis- and trans-regulatory target genes, we discovered 20 hDels controlling gene expression. We highlight two hDels, hDel_2247 and hDel_585, with tissue-specific activity in the brain. Our findings reveal a molecular and cellular role for sequences lost in the human lineage and establish a framework for functionally interrogating human-specific genetic variants.

DNA-binding factor footprints and enhancer RNAs identify functional non-coding genetic variants

BACKGROUND

Genome-wide association studies (GWAS) have revealed a multitude of candidate genetic variants affecting the risk of developing complex traits and diseases. However, the highlighted regions are typically in the non-coding genome, and uncovering the functional causative single nucleotide variants (SNVs) is challenging. Prioritization of variants is commonly based on genomic annotation with markers of active regulatory elements, but current approaches still poorly predict functional variants. To address this, we systematically analyze six markers of active regulatory elements for their ability to identify functional variants.

RESULTS

We benchmark against molecular quantitative trait loci (molQTL) from assays of regulatory element activity that identify allelic effects on DNA-binding factor occupancy, reporter assay expression, and chromatin accessibility. We identify the combination of DNase footprints and divergent enhancer RNA (eRNA) as markers for functional variants. This signature provides high precision, but with a trade-off of low recall, thus substantially reducing candidate variant sets to prioritize variants for functional validation. We present this as a framework called FINDER-Functional SNV IdeNtification using DNase footprints and eRNA.

CONCLUSIONS

We demonstrate the utility to prioritize variants using leukocyte count trait and analyze variants in linkage disequilibrium with a lead variant to predict a functional variant in asthma. Our findings have implications for prioritizing variants from GWAS, in development of predictive scoring algorithms, and for functionally informed fine mapping approaches.

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.

Modifying peptide/lipid-associated nucleic acids (PLANAs) for CRISPR/Cas9 ribonucleoprotein delivery

With the first reports on the possibility of genome editing by Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) and CRISPR-associated protein (Cas)9 surfacing in 2005, the enthusiasm for protein silencing via nucleic acid delivery experienced a resurgence following a period of diminished enthusiasm due to challenges in delivering small interfering RNAs (siRNA), especially in vivo. However, delivering the components necessary for this approach into the nucleus is challenging, maybe even more than the cytoplasmic delivery of siRNA. We previously reported the birth of peptide/lipid-associated nucleic acids (PLANAs) for siRNA delivery. This project was designed to investigate the efficiency of these nanoparticles for in vitro delivery of CRISPR/Cas9 ribonucleoproteins. Our initial experiments indicated higher toxicity for PLANAs with the more efficient reverse transfection method. Therefore, polyethylene glycol (PEG) was added to the composition for PEGylation of the nanoparticles by partially replacing two of the lipid components with the PEG-conjugated counterparts. The results indicated a more significant reduction in the toxicity of the nanoparticle, less compromise in encapsulation efficiency and more PEGylation of the surface of the nanoparticles using DOPE-PEG2000 at 50 % replacement of the naïve lipid. The cell internalization and transfection efficiency showed a comparable efficiency for the PEGylated and non-PEGylated PLANAs and the commercially available Lipofectamine™ CRISPRMAX™. Next Generation Sequencing of the cloned cells showed a variety of indels in the transfected cell population. Overall, our results indicate the efficiency and safety of PEGylated PLANAs for in vitro transfection with CRISPR/Cas9 ribonucleoproteins. PEGylation has been studied extensively for in vivo delivery, and PEGylated PLANAs will be candidates for future in vivo studies.

Multicenter integrated analysis of noncoding CRISPRi screens

The ENCODE Consortium's efforts to annotate noncoding cis-regulatory elements (CREs) have advanced our understanding of gene regulatory landscapes. Pooled, noncoding CRISPR screens offer a systematic approach to investigate cis-regulatory mechanisms. The ENCODE4 Functional Characterization Centers conducted 108 screens in human cell lines, comprising >540,000 perturbations across 24.85 megabases of the genome. Using 332 functionally confirmed CRE-gene links in K562 cells, we established guidelines for screening endogenous noncoding elements with CRISPR interference (CRISPRi), including accurate detection of CREs that exhibit variable, often low, transcriptional effects. Benchmarking five screen analysis tools, we find that CASA produces the most conservative CRE calls and is robust to artifacts of low-specificity single guide RNAs. We uncover a subtle DNA strand bias for CRISPRi in transcribed regions with implications for screen design and analysis. Together, we provide an accessible data resource, predesigned single guide RNAs for targeting 3,275,697 ENCODE SCREEN candidate CREs with CRISPRi and screening guidelines to accelerate functional characterization of the noncoding genome.

Mapping the functional impact of non-coding regulatory elements in primary T cells through single-cell CRISPR screens

BACKGROUND

Drug targets with genetic evidence are expected to increase clinical success by at least twofold. Yet, translating disease-associated genetic variants into functional knowledge remains a fundamental challenge of drug discovery. A key issue is that the vast majority of complex disease associations cannot be cleanly mapped to a gene. Immune disease-associated variants are enriched within regulatory elements found in T-cell-specific open chromatin regions.

RESULTS

To identify genes and molecular programs modulated by these regulatory elements, we develop a CRISPRi-based single-cell functional screening approach in primary human T cells. Our pipeline enables the interrogation of transcriptomic changes induced by the perturbation of regulatory elements at scale. We first optimize an efficient CRISPRi protocol in primary CD4+ T cells via CROPseq vectors. Subsequently, we perform a screen targeting 45 non-coding regulatory elements and 35 transcription start sites and profile approximately 250,000 T -cell single-cell transcriptomes. We develop a bespoke analytical pipeline for element-to-gene (E2G) mapping and demonstrate that our method can identify both previously annotated and novel E2G links. Lastly, we integrate genetic association data for immune-related traits and demonstrate how our platform can aid in the identification of effector genes for GWAS loci.

CONCLUSIONS

We describe "primary T cell crisprQTL" - a scalable, single-cell functional genomics approach for mapping regulatory elements to genes in primary human T cells. We show how this framework can facilitate the interrogation of immune disease GWAS hits and propose that the combination of experimental and QTL-based techniques is likely to address the variant-to-function problem.

Multiplexed screening reveals how cancer-specific alternative polyadenylation shapes tumor growth in vivo

Alternative polyadenylation (APA) is strikingly dysregulated in many cancers. Although global APA dysregulation is frequently associated with poor prognosis, the importance of most individual APA events is controversial simply because few have been functionally studied. Here, we address this gap by developing a CRISPR-Cas9-based screen to manipulate endogenous polyadenylation and systematically quantify how APA events contribute to tumor growth in vivo. Our screen reveals individual APA events that control mouse melanoma growth in an immunocompetent host, with concordant associations in clinical human cancer. For example, forced Atg7 3' UTR lengthening in mouse melanoma suppresses ATG7 protein levels, slows tumor growth, and improves host survival; similarly, in clinical human melanoma, a long ATG7 3' UTR is associated with significantly prolonged patient survival. Overall, our study provides an easily adaptable means to functionally dissect APA in physiological systems and directly quantifies the contributions of recurrent APA events to tumorigenic phenotypes.

Parallelized engineering of mutational models using piggyBac transposon delivery of CRISPR libraries

New technologies and large-cohort studies have enabled novel variant discovery and association at unprecedented scale, yet functional characterization of these variants remains paramount to deciphering disease mechanisms. Approaches that facilitate parallelized genome editing of cells of interest or induced pluripotent stem cells (iPSCs) have become critical tools toward this goal. Here, we developed an approach that incorporates libraries of CRISPR-Cas9 guide RNAs (gRNAs) together with inducible Cas9 into a piggyBac (PB) transposon system to engineer dozens to hundreds of genomic variants in parallel against isogenic cellular backgrounds. This method empowers loss-of-function (LoF) studies through the introduction of insertions or deletions (indels) and copy-number variants (CNVs), though generating specific nucleotide changes is possible with prime editing. The ability to rapidly establish high-quality mutational models at scale will facilitate the development of isogenic cellular collections and catalyze comparative functional genomic studies investigating the roles of hundreds of genes and mutations in development and disease.

A Q-marker screening strategy based on ADME studies and systems biology for Chinese herbal medicine, taking Qianghuo Shengshi decoction in treating rheumatoid arthritis as an example

Chinese herbal medicine (CHM) exhibits a broad spectrum of clinical applications and demonstrates favorable therapeutic efficacy. Nonetheless, elucidating the underlying mechanism of action (MOA) of CHM in disease treatment remains a formidable task due to its inherent characteristics of multi-level, multi-linked, and multi-dimensional non-linear synergistic actions. In recent years, the concept of a Quality marker (Q-marker) proposed by Liu et al. has significantly contributed to the monitoring and evaluation of CHM products, thereby fostering the advancement of CHM research. Within this study, a Q-marker screening strategy for CHM formulas has been introduced, particularly emphasising efficacy and biological activities, integrating absorption, distribution, metabolism, and excretion (ADME) studies, systems biology, and experimental verification. As an illustrative case, the Q-marker screening of Qianghuo Shengshi decoction (QHSSD) for treating rheumatoid arthritis (RA) has been conducted. Consequently, from a pool of 159 compounds within QHSSD, five Q-markers exhibiting significant in vitro anti-inflammatory effects have been identified. These Q-markers encompass notopterol, isoliquiritin, imperatorin, cimifugin, and glycyrrhizic acid. Furthermore, by employing an integrated analysis of network pharmacology and metabolomics, several instructive insights into pharmacological mechanisms have been gleaned. This includes the identification of key targets and pathways through which QHSSD exerts its crucial roles in the treatment of RA. Notably, the inhibitory effect of QHSSD on AKT1 and MAPK3 activation has been validated through western blot analysis, underscoring its potential to mitigate RA-related inflammatory responses. In summary, this research demonstrates the proposed strategy's feasibility and provides a practical reference model for the systematic investigation of CHM formulas.

CRISPR-based epigenome editing: mechanisms and applications

Epigenomic anomalies contribute significantly to the development of numerous human disorders. The development of epigenetic research tools is essential for understanding how epigenetic marks contribute to gene expression. A gene-editing technique known as CRISPR (clustered regularly interspaced short palindromic repeats) typically targets a particular DNA sequence using a guide RNA (gRNA). CRISPR/Cas9 technology has been remodeled for epigenome editing by generating a 'dead' Cas9 protein (dCas9) that lacks nuclease activity and juxtaposing it with an epigenetic effector domain. Based on fusion partners of dCas9, a specific epigenetic state can be achieved. CRISPR-based epigenome editing has widespread application in drug screening, cancer treatment and regenerative medicine. This paper discusses the tools developed for CRISPR-based epigenome editing and their applications.

Combined CRISPRi and proteomics screening reveal a cohesin-CTCF-bound allele contributing to increased expression of RUVBL1 and prostate cancer progression

Genome-wide association studies along with expression quantitative trait locus (eQTL) mapping have identified hundreds of single-nucleotide polymorphisms (SNPs) and their target genes in prostate cancer (PCa), yet functional characterization of these risk loci remains challenging. To screen for potential regulatory SNPs, we designed a CRISPRi library containing 9,133 guide RNAs (gRNAs) to cover 2,166 candidate SNP loci implicated in PCa and identified 117 SNPs that could regulate 90 genes for PCa cell growth advantage. Among these, rs60464856 was covered by multiple gRNAs significantly depleted in screening (FDR < 0.05). Pooled SNP association analysis in the PRACTICAL and FinnGen cohorts showed significantly higher PCa risk for the rs60464856 G allele (p value = 1.2 × 10-16 and 3.2 × 10-7, respectively). Subsequent eQTL analysis revealed that the G allele is associated with increased RUVBL1 expression in multiple datasets. Further CRISPRi and xCas9 base editing confirmed that the rs60464856 G allele leads to elevated RUVBL1 expression. Furthermore, SILAC-based proteomic analysis demonstrated allelic binding of cohesin subunits at the rs60464856 region, where the HiC dataset showed consistent chromatin interactions in prostate cell lines. RUVBL1 depletion inhibited PCa cell proliferation and tumor growth in a xenograft mouse model. Gene-set enrichment analysis suggested an association of RUVBL1 expression with cell-cycle-related pathways. Increased expression of RUVBL1 and activation of cell-cycle pathways were correlated with poor PCa survival in TCGA datasets. Our CRISPRi screening prioritized about one hundred regulatory SNPs essential for prostate cell proliferation. In combination with proteomics and functional studies, we characterized the mechanistic role of rs60464856 and RUVBL1 in PCa progression.

Dynamic network-guided CRISPRi screen identifies CTCF-loop-constrained nonlinear enhancer gene regulatory activity during cell state transitions

Comprehensive enhancer discovery is challenging because most enhancers, especially those contributing to complex diseases, have weak effects on gene expression. Our gene regulatory network modeling identified that nonlinear enhancer gene regulation during cell state transitions can be leveraged to improve the sensitivity of enhancer discovery. Using human embryonic stem cell definitive endoderm differentiation as a dynamic transition system, we conducted a mid-transition CRISPRi-based enhancer screen. We discovered a comprehensive set of enhancers for each of the core endoderm-specifying transcription factors. Many enhancers had strong effects mid-transition but weak effects post-transition, consistent with the nonlinear temporal responses to enhancer perturbation predicted by the modeling. Integrating three-dimensional genomic information, we were able to develop a CTCF-loop-constrained Interaction Activity model that can better predict functional enhancers compared to models that rely on Hi-C-based enhancer-promoter contact frequency. Our study provides generalizable strategies for sensitive and systematic enhancer discovery in both normal and pathological cell state transitions.

Evaluating the mouse neural precursor line, SN4741, as a suitable proxy for midbrain dopaminergic neurons

To overcome the ethical and technical limitations of in vivo human disease models, the broader scientific community frequently employs model organism-derived cell lines to investigate disease mechanisms, pathways, and therapeutic strategies. Despite the widespread use of certain in vitro models, many still lack contemporary genomic analysis supporting their use as a proxy for the affected human cells and tissues. Consequently, it is imperative to determine how accurately and effectively any proposed biological surrogate may reflect the biological processes it is assumed to model. One such cellular surrogate of human disease is the established mouse neural precursor cell line, SN4741, which has been used to elucidate mechanisms of neurotoxicity in Parkinson disease for over 25 years. Here, we are using a combination of classic and contemporary genomic techniques - karyotyping, RT-qPCR, single cell RNA-seq, bulk RNA-seq, and ATAC-seq - to characterize the transcriptional landscape, chromatin landscape, and genomic architecture of this cell line, and evaluate its suitability as a proxy for midbrain dopaminergic neurons in the study of Parkinson disease. We find that SN4741 cells possess an unstable triploidy and consistently exhibits low expression of dopaminergic neuron markers across assays, even when the cell line is shifted to the non-permissive temperature that drives differentiation. The transcriptional signatures of SN4741 cells suggest that they are maintained in an undifferentiated state at the permissive temperature and differentiate into immature neurons at the non-permissive temperature; however, they may not be dopaminergic neuron precursors, as previously suggested. Additionally, the chromatin landscapes of SN4741 cells, in both the differentiated and undifferentiated states, are not concordant with the open chromatin profiles of ex vivo, mouse E15.5 forebrain- or midbrain-derived dopaminergic neurons. Overall, our data suggest that SN4741 cells may reflect early aspects of neuronal differentiation but are likely not a suitable proxy for dopaminergic neurons as previously thought. The implications of this study extend broadly, illuminating the need for robust biological and genomic rationale underpinning the use of in vitro models of molecular processes.

Combined CRISPRi and proteomics screening reveal a cohesin-CTCF-bound allele contributing to increased expression of RUVBL1 and prostate cancer progression

Genome-wide association studies along with expression quantitative trait loci (eQTL) mapping have identified hundreds of single nucleotide polymorphisms (SNPs) and their target genes in prostate cancer (PCa), yet functional characterization of these risk loci remains challenging. To screen for potential regulatory SNPs, we designed a CRISPRi library containing 9133 guide RNAs (gRNAs) to target 2,166 candidate SNP sites implicated in PCa and identified 117 SNPs that could regulate 90 genes for PCa cell growth advantage. Among these, rs60464856 was covered by multiple gRNAs significantly depleted in the screening (FDR<0.05). Pooled SNP association analysis in the PRACTICAL and FinnGen cohorts showed significantly higher PCa risk for the rs60464856 G allele (pvalue=1.2E-16 and 3.2E-7). Subsequent eQTL analysis revealed that the G allele is associated with increased RUVBL1 expression in multiple datasets. Further CRISPRi and xCas9 base editing proved the rs60464856 G allele leading to an elevated RUVBL1 expression. Furthermore, SILAC-based proteomic analysis demonstrated allelic binding of cohesin subunits at the rs60464856 region, where HiC dataset showed consistent chromatin interactions in prostate cell lines. RUVBL1 depletion inhibited PCa cell proliferation and tumor growth in xenograft mouse model. Gene set enrichment analysis suggested an association of RUVBL1 expression with cell-cycle-related pathways. An increased expression of RUVBL1 and activations of cell-cycle pathways were correlated with poor PCa survival in TCGA datasets. Together, our CRISPRi screening prioritized about one hundred regulatory SNPs essential for prostate cell proliferation. In combination with proteomics and functional studies, we characterized the mechanistic role of rs60464856 and RUVBL1 in PCa progression.

CHD-associated enhancers shape human cardiomyocyte lineage commitment

Enhancers orchestrate gene expression programs that drive multicellular development and lineage commitment. Thus, genetic variants at enhancers are thought to contribute to developmental diseases by altering cell fate commitment. However, while many variant-containing enhancers have been identified, studies to endogenously test the impact of these enhancers on lineage commitment have been lacking. We perform a single-cell CRISPRi screen to assess the endogenous roles of 25 enhancers and putative cardiac target genes implicated in genetic studies of congenital heart defects (CHDs). We identify 16 enhancers whose repression leads to deficient differentiation of human cardiomyocytes (CMs). A focused CRISPRi validation screen shows that repression of TBX5 enhancers delays the transcriptional switch from mid- to late-stage CM states. Endogenous genetic deletions of two TBX5 enhancers phenocopy epigenetic perturbations. Together, these results identify critical enhancers of cardiac development and suggest that misregulation of these enhancers could contribute to cardiac defects in human patients.

Deletion mapping of regulatory elements for GATA3 in T cells reveals a distal enhancer involved in allergic diseases

GATA3 is essential for T cell differentiation and is surrounded by genome-wide association study (GWAS) hits for immune traits. Interpretation of these GWAS hits is challenging because gene expression quantitative trait locus (eQTL) studies lack power to detect variants with small effects on gene expression in specific cell types and the genome region containing GATA3 contains dozens of potential regulatory sequences. To map regulatory sequences for GATA3, we performed a high-throughput tiling deletion screen of a 2 Mb genome region in Jurkat T cells. This revealed 23 candidate regulatory sequences, all but one of which is within the same topological-associating domain (TAD) as GATA3. We then performed a lower-throughput deletion screen to precisely map regulatory sequences in primary T helper 2 (Th2) cells. We tested 25 sequences with ∼100 bp deletions and validated five of the strongest hits with independent deletion experiments. Additionally, we fine-mapped GWAS hits for allergic diseases in a distal regulatory element, 1 Mb downstream of GATA3, and identified 14 candidate causal variants. Small deletions spanning the candidate variant rs725861 decreased GATA3 levels in Th2 cells, and luciferase reporter assays showed regulatory differences between its two alleles, suggesting a causal mechanism for this variant in allergic diseases. Our study demonstrates the power of integrating GWAS signals with deletion mapping and identifies critical regulatory sequences for GATA3.

Genetic architecture of human brain evolution

Comparative studies of hominids have long sought to identify mutational events that shaped the evolution of the human nervous system. However, functional genetic differences are outnumbered by millions of nearly neutral mutations, and the developmental mechanisms underlying human nervous system specializations are difficult to model and incompletely understood. Candidate-gene studies have attempted to map select human-specific genetic differences to neurodevelopmental functions, but it remains unclear how to contextualize the relative effects of genes that are investigated independently. Considering these limitations, we discuss scalable approaches for probing the functional contributions of human-specific genetic differences. We propose that a systems-level view will enable a more quantitative and integrative understanding of the genetic, molecular and cellular underpinnings of human nervous system evolution.

Dynamic network-guided CRISPRi screen reveals CTCF loop-constrained nonlinear enhancer-gene regulatory activity in cell state transitions

Comprehensive enhancer discovery is challenging because most enhancers, especially those affected in complex diseases, have weak effects on gene expression. Our network modeling revealed that nonlinear enhancer-gene regulation during cell state transitions can be leveraged to improve the sensitivity of enhancer discovery. Utilizing hESC definitive endoderm differentiation as a dynamic transition system, we conducted a mid-transition CRISPRi-based enhancer screen. The screen discovered a comprehensive set of enhancers (4 to 9 per locus) for each of the core endoderm lineage-specifying transcription factors, and many enhancers had strong effects mid-transition but weak effects post-transition. Through integrating enhancer activity measurements and three-dimensional enhancer-promoter interaction information, we were able to develop a CTCF loop-constrained Interaction Activity (CIA) model that can better predict functional enhancers compared to models that rely on Hi-C-based enhancer-promoter contact frequency. Our study provides generalizable strategies for sensitive and more comprehensive enhancer discovery in both normal and pathological cell state transitions.

Translating non-coding genetic associations into a better understanding of immune-mediated disease

Genome-wide association studies have identified hundreds of genetic loci that are associated with immune-mediated diseases. Most disease-associated variants are non-coding, and a large proportion of these variants lie within enhancers. As a result, there is a pressing need to understand how common genetic variation might affect enhancer function and thereby contribute to immune-mediated (and other) diseases. In this Review, we first describe statistical and experimental methods to identify causal genetic variants that modulate gene expression, including statistical fine-mapping and massively parallel reporter assays. We then discuss approaches to characterise the mechanisms by which these variants modulate immune function, such as clustered regularly interspaced short palindromic repeats (CRISPR)-based screens. We highlight examples of studies that, by elucidating the effects of disease variants within enhancers, have provided important insights into immune function and uncovered key pathways of disease.

Evaluating the mouse neural precursor line, SN4741, as a suitable proxy for midbrain dopaminergic neurons

To overcome the ethical and technical limitations of in vivo human disease models, the broader scientific community frequently employs model organism-derived cell lines to investigate of disease mechanisms, pathways, and therapeutic strategies. Despite the widespread use of certain in vitro models, many still lack contemporary genomic analysis supporting their use as a proxy for the affected human cells and tissues. Consequently, it is imperative to determine how accurately and effectively any proposed biological surrogate may reflect the biological processes it is assumed to model. One such cellular surrogate of human disease is the established mouse neural precursor cell line, SN4741, which has been used to elucidate mechanisms of neurotoxicity in Parkinson disease for over 25 years. Here, we are using a combination of classic and contemporary genomic techniques - karyotyping, RT-qPCR, single cell RNA-seq, bulk RNA-seq, and ATAC-seq - to characterize the transcriptional landscape, chromatin landscape, and genomic architecture of this cell line, and evaluate its suitability as a proxy for midbrain dopaminergic neurons in the study of Parkinson disease. We find that SN4741 cells possess an unstable triploidy and consistently exhibits low expression of dopaminergic neuron markers across assays, even when the cell line is shifted to the non-permissive temperature that drives differentiation. The transcriptional signatures of SN4741 cells suggest that they are maintained in an undifferentiated state at the permissive temperature and differentiate into immature neurons at the non-permissive temperature; however, they may not be dopaminergic neuron precursors, as previously suggested. Additionally, the chromatin landscapes of SN4741 cells, in both the differentiated and undifferentiated states, are not concordant with the open chromatin profiles of ex vivo , mouse E15.5 forebrain- or midbrain-derived dopaminergic neurons. Overall, our data suggest that SN4741 cells may reflect early aspects of neuronal differentiation but are likely not a suitable a proxy for dopaminergic neurons as previously thought. The implications of this study extend broadly, illuminating the need for robust biological and genomic rationale underpinning the use of in vitro models of molecular processes.

Evaluating the mouse neural precursor line, SN4741, as a suitable proxy for midbrain dopaminergic neurons

To overcome the ethical and technical limitations of in vivo human disease models, the broader scientific community frequently employs model organism-derived cell lines to investigate of disease mechanisms, pathways, and therapeutic strategies. Despite the widespread use of certain in vitro models, many still lack contemporary genomic analysis supporting their use as a proxy for the affected human cells and tissues. Consequently, it is imperative to determine how accurately and effectively any proposed biological surrogate may reflect the biological processes it is assumed to model. One such cellular surrogate of human disease is the established mouse neural precursor cell line, SN4741, which has been used to elucidate mechanisms of neurotoxicity in Parkinson disease for over 25 years. Here, we are using a combination of classic and contemporary genomic techniques - karyotyping, RT-qPCR, single cell RNA-seq, bulk RNA-seq, and ATAC-seq - to characterize the transcriptional landscape, chromatin landscape, and genomic architecture of this cell line, and evaluate its suitability as a proxy for midbrain dopaminergic neurons in the study of Parkinson disease. We find that SN4741 cells possess an unstable triploidy and consistently exhibits low expression of dopaminergic neuron markers across assays, even when the cell line is shifted to the non-permissive temperature that drives differentiation. The transcriptional signatures of SN4741 cells suggest that they are maintained in an undifferentiated state at the permissive temperature and differentiate into immature neurons at the non-permissive temperature; however, they may not be dopaminergic neuron precursors, as previously suggested. Additionally, the chromatin landscapes of SN4741 cells, in both the differentiated and undifferentiated states, are not concordant with the open chromatin profiles of ex vivo , mouse E15.5 forebrain- or midbrain-derived dopaminergic neurons. Overall, our data suggest that SN4741 cells may reflect early aspects of neuronal differentiation but are likely not a suitable a proxy for dopaminergic neurons as previously thought. The implications of this study extend broadly, illuminating the need for robust biological and genomic rationale underpinning the use of in vitro models of molecular processes.

Norbornene-tetrazine ligation chemistry for controlling RNA-guided CRISPR systems

In the present study, norbornene-tetrazine ligation chemistry is harnessed for developing clickable RNA switches in biological contexts. This RNA control strategy is explored with a variety of applications. We further demonstrate the application of RNA-based norbornene-tetrazine ligation chemistry for controlling CRISPR systems. Moreover, the manipulation of gene editing in human cells is accomplished.

G-quadruplex-guided RNA engineering to modulate CRISPR-based genomic regulation

It is important to develop small moelcule-based methods to modulate gene editing and expression in human cells. The roles of the G-quadruplex (G4) in biological systems have been widely studied. Here, G4-guided RNA engineering is performed to generate guide RNA with G4-forming units (G4-gRNA). We further demonstrate that chemical targeting of G4-gRNAs holds promise as a general approach for modulating gene editing and expression in human cells. The rich structural diversity of RNAs offers a reservoir of targets for small molecules to bind, thus creating the potential to modulate RNA biology.

Functional genomic assays to annotate enhancer-promoter interactions genome wide

Enhancers are pivotal for regulating gene transcription that occurs at promoters. Identification of the interacting enhancer-promoter pairs and understanding the mechanisms behind how they interact and how enhancers modulate transcription can provide fundamental insight into gene regulatory networks. Recently, advances in high-throughput methods in three major areas-chromosome conformation capture assay, such as Hi-C to study basic chromatin architecture, ectopic reporter experiments such as self-transcribing active regulatory region sequencing (STARR-seq) to quantify promoter and enhancer activity, and endogenous perturbations such as clustered regularly interspaced short palindromic repeat interference (CRISPRi) to identify enhancer-promoter compatibility-have further our knowledge about transcription. In this review, we will discuss the major method developments and key findings from these assays.

Demystifying non-coding GWAS variants: an overview of computational tools and methods

Genome-wide association studies (GWAS) have found the majority of disease-associated variants to be non-coding. Major efforts into the charting of the non-coding regulatory landscapes have allowed for the development of tools and methods which aim to aid in the identification of causal variants and their mechanism of action. In this review, we give an overview of current tools and methods for the analysis of non-coding GWAS variants in disease. We provide a workflow that allows for the accumulation of in silico evidence to generate novel hypotheses on mechanisms underlying disease and prioritize targets for follow-up study using non-coding GWAS variants. Lastly, we discuss the need for comprehensive benchmarks and novel tools for the analysis of non-coding variants.

A Bcl6 Intronic Element Regulates T Follicular Helper Cell Differentiation

In response to an intracellular infectious agent, the immune system produces a specific cellular response as well as a T cell–dependent Ab response. Precursor T cells differentiate into effector T cells, including Th1 cells, and T follicular helper (TFH) cells. The latter cooperate with B cells to form germinal centers and induce the formation of Ab-forming plasmacytes. One major focal point for control of T cell differentiation is the transcription factor BCL6. In this study, we demonstrated that the Bcl6 gene is regulated by FOXO1-binding, cis-acting sequences located in a highly conserved region of the first Bcl6 intron. In both mouse and human T cells, deletion of the tandem FOXO1 binding sites increased the expression of BCL6 and enhanced the proportion of TFH cells. These results reveal a fundamental control point for cellular versus humoral immunity.

Functional inference of gene regulation using single-cell multi-omics

Cells require coordinated control over gene expression when responding to environmental stimuli. Here we apply scATAC-seq and single-cell RNA sequencing (scRNA-seq) in resting and stimulated human blood cells. Collectively, we generate ~91,000 single-cell profiles, allowing us to probe the cis-regulatory landscape of the immunological response across cell types, stimuli, and time. Advancing tools to integrate multi-omics data, we develop functional inference of gene regulation (FigR), a framework to computationally pair scA-TAC-seq with scRNA-seq cells, connect distal cis-regulatory elements to genes, and infer gene-regulatory networks (GRNs) to identify candidate transcription factor (TF) regulators. Utilizing these paired multi-omics data, we define domains of regulatory chromatin (DORCs) of immune stimulation and find that cells alter chromatin accessibility and gene expression at timescales of minutes. Construction of the stimulation GRN elucidates TF activity at disease-associated DORCs. Overall, FigR enables elucidation of regulatory interactions across single-cell data, providing new opportunities to understand the function of cells within tissues.

Recommendations for clinical interpretation of variants found in non-coding regions of the genome

BACKGROUND

The majority of clinical genetic testing focuses almost exclusively on regions of the genome that directly encode proteins. The important role of variants in non-coding regions in penetrant disease is, however, increasingly being demonstrated, and the use of whole genome sequencing in clinical diagnostic settings is rising across a large range of genetic disorders. Despite this, there is no existing guidance on how current guidelines designed primarily for variants in protein-coding regions should be adapted for variants identified in other genomic contexts.

METHODS

We convened a panel of nine clinical and research scientists with wide-ranging expertise in clinical variant interpretation, with specific experience in variants within non-coding regions. This panel discussed and refined an initial draft of the guidelines which were then extensively tested and reviewed by external groups.

RESULTS

We discuss considerations specifically for variants in non-coding regions of the genome. We outline how to define candidate regulatory elements, highlight examples of mechanisms through which non-coding region variants can lead to penetrant monogenic disease, and outline how existing guidelines can be adapted for the interpretation of these variants.

CONCLUSIONS

These recommendations aim to increase the number and range of non-coding region variants that can be clinically interpreted, which, together with a compatible phenotype, can lead to new diagnoses and catalyse the discovery of novel disease mechanisms.

Synthetic regulatory reconstitution reveals principles of mammalian Hox cluster regulation

Precise Hox gene expression is crucial for embryonic patterning. Intra-Hox transcription factor binding and distal enhancer elements have emerged as the major regulatory modules controlling Hox gene expression. However, quantifying their relative contributions has remained elusive. Here, we introduce "synthetic regulatory reconstitution," a conceptual framework for studying gene regulation, and apply it to the HoxA cluster. We synthesized and delivered variant rat HoxA clusters (130 to 170 kilobases) to an ectopic location in the mouse genome. We found that a minimal HoxA cluster recapitulated correct patterns of chromatin remodeling and transcription in response to patterning signals, whereas the addition of distal enhancers was needed for full transcriptional output. Synthetic regulatory reconstitution could provide a generalizable strategy for deciphering the regulatory logic of gene expression in complex genomes.

Designing libraries for pooled CRISPR functional screens of long noncoding RNAs

Human and other genomes encode tens of thousands of long noncoding RNAs (lncRNAs), the vast majority of which remain uncharacterised. High-throughput functional screening methods, notably those based on pooled CRISPR-Cas perturbations, promise to unlock the biological significance and biomedical potential of lncRNAs. Such screens are based on libraries of single guide RNAs (sgRNAs) whose design is critical for success. Few off-the-shelf libraries are presently available, and lncRNAs tend to have cell-type-specific expression profiles, meaning that library design remains in the hands of researchers. Here we introduce the topic of pooled CRISPR screens for lncRNAs and guide readers through the three key steps of library design: accurate annotation of transcript structures, curation of optimal candidate sets, and design of sgRNAs. This review is a starting point and reference for researchers seeking to design custom CRISPR screening libraries for lncRNAs.

Rational guide RNA engineering for small-molecule control of CRISPR/Cas9 and gene editing

It is important to control CRISPR/Cas9 when sufficient editing is obtained. In the current study, rational engineering of guide RNAs (gRNAs) is performed to develop small-molecule-responsive CRISPR/Cas9. For our purpose, the sequence of gRNAs are modified to introduce ligand binding sites based on the rational design of ligand-RNA pairs. Using short target sequences, we demonstrate that the engineered RNA provides an excellent scaffold for binding small molecule ligands. Although the 'stem-loop 1' variants of gRNA induced variable cleavage activity for different target sequences, all 'stem-loop 3' variants are well tolerated for CRISPR/Cas9. We further demonstrate that this specific ligand-RNA interaction can be utilized for functional control of CRISPR/Cas9 in vitro and in human cells. Moreover, chemogenetic control of gene editing in human cells transfected with all-in-one plasmids encoding Cas9 and designer gRNAs is demonstrated. The strategy may become a general approach for generating switchable RNA or DNA for controlling other biological processes.

Programmable deletion, replacement, integration and inversion of large DNA sequences with twin prime editing

The targeted deletion, replacement, integration or inversion of genomic sequences could be used to study or treat human genetic diseases, but existing methods typically require double-strand DNA breaks (DSBs) that lead to undesired consequences, including uncontrolled indel mixtures and chromosomal abnormalities. Here we describe twin prime editing (twinPE), a DSB-independent method that uses a prime editor protein and two prime editing guide RNAs (pegRNAs) for the programmable replacement or excision of DNA sequences at endogenous human genomic sites. The two pegRNAs template the synthesis of complementary DNA flaps on opposing strands of genomic DNA, which replace the endogenous DNA sequence between the prime-editor-induced nick sites. When combined with a site-specific serine recombinase, twinPE enabled targeted integration of gene-sized DNA plasmids (>5,000 bp) and targeted sequence inversions of 40 kb in human cells. TwinPE expands the capabilities of precision gene editing and might synergize with other tools for the correction or complementation of large or complex human pathogenic alleles.

Gene editing and its applications in biomedicine

The steady progress in genome editing, especially genome editing based on the use of clustered regularly interspaced short palindromic repeats (CRISPR) and programmable nucleases to make precise modifications to genetic material, has provided enormous opportunities to advance biomedical research and promote human health. The application of these technologies in basic biomedical research has yielded significant advances in identifying and studying key molecular targets relevant to human diseases and their treatment. The clinical translation of genome editing techniques offers unprecedented biomedical engineering capabilities in the diagnosis, prevention, and treatment of disease or disability. Here, we provide a general summary of emerging biomedical applications of genome editing, including open challenges. We also summarize the tools of genome editing and the insights derived from their applications, hoping to accelerate new discoveries and therapies in biomedicine.

Supramolecular CRISPR-OFF switches with host-guest chemistry

CRISPR (Clustered Regularly Interspaced Short Palindromic Repeat) technology is a powerful tool in biology and medicine. However, the safety and application of this technology is hampered by excessive activity of CRISPR machinery. It is particularly important to develop methods for switching off CRISPR activity in human cells. The current study demonstrates the concept of supramolecular CRISPR-OFF switches by employing host-guest chemistry. We demonstrate that the CRISPR systems show considerable tolerance to adamantoylation on guide RNAs (gRNAs), whereas supramolecular complexation tremendously affects the function of adamantoyl gRNAs. Host-guest chemistry is demonstrated to be novel and effective tools to reduce unwanted excessive activities of CRISPR complexes in human cells. This work indicates considerable potential of supramolecular strategy for controlling and enhancing CRISPR systems.

Precise genomic deletions using paired prime editing

Current methods to delete genomic sequences are based on clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 and pairs of single-guide RNAs (sgRNAs), but can be inefficient and imprecise, with errors including small indels as well as unintended large deletions and more complex rearrangements. In the present study, we describe a prime editing-based method, PRIME-Del, which induces a deletion using a pair of prime editing sgRNAs (pegRNAs) that target opposite DNA strands, programming not only the sites that are nicked but also the outcome of the repair. PRIME-Del achieves markedly higher precision than CRISPR-Cas9 and sgRNA pairs in programming deletions up to 10 kb, with 1-30% editing efficiency. PRIME-Del can also be used to couple genomic deletions with short insertions, enabling deletions with junctions that do not fall at protospacer-adjacent motif sites. Finally, extended expression of prime editing components can substantially enhance efficiency without compromising precision. We anticipate that PRIME-Del will be broadly useful for precise, flexible programming of genomic deletions, epitope tagging and, potentially, programming genomic rearrangements.

CRISPR-Based Approaches for the High-Throughput Characterization of Long Non-Coding RNAs

Over the last decade, tens of thousands of new long non-coding RNAs (lncRNAs) have been identified in the human genome. Nevertheless, except for a handful of genes, the genetic characteristics and functions of most of these lncRNAs remain elusive; this is partially due to their relatively low expression, high tissue specificity, and low conservation across species. A major limitation for determining the function of lncRNAs was the lack of methodologies suitable for studying these genes. The recent development of CRISPR/Cas9 technology has opened unprecedented opportunities to uncover the genetic and functional characteristics of the non-coding genome via targeted and high-throughput approaches. Specific CRISPR/Cas9-based approaches were developed to target lncRNA loci. Some of these approaches involve modifying the sequence, but others were developed to study lncRNAs by inducing transcriptional and epigenetic changes. The discovery of other programable Cas proteins broaden our possibilities to target RNA molecules with greater precision and accuracy. These approaches allow for the knock-down and characterization of lncRNAs. Here, we review how various CRISPR-based strategies have been used to characterize lncRNAs with important functions in different biological contexts and how these approaches can be further utilized to improve our understanding of the non-coding genome.

Assessing genome-wide dynamic changes in enhancer activity during early mESC differentiation by FAIRE-STARR-seq

Embryonic stem cells (ESCs) can differentiate into any given cell type and therefore represent a versatile model to study the link between gene regulation and differentiation. To quantitatively assess the dynamics of enhancer activity during the early stages of murine ESC differentiation, we analyzed accessible genomic regions using STARR-seq, a massively parallel reporter assay. This resulted in a genome-wide quantitative map of active mESC enhancers, in pluripotency and during the early stages of differentiation. We find that only a minority of accessible regions is active and that such regions are enriched near promoters, characterized by specific chromatin marks, enriched for distinct sequence motifs, and modeling shows that active regions can be predicted from sequence alone. Regions that change their activity upon retinoic acid-induced differentiation are more prevalent at distal intergenic regions when compared to constitutively active enhancers. Further, analysis of differentially active enhancers verified the contribution of individual TF motifs toward activity and inducibility as well as their role in regulating endogenous genes. Notably, the activity of retinoic acid receptor alpha (RARα) occupied regions can either increase or decrease upon the addition of its ligand, retinoic acid, with the direction of the change correlating with spacing and orientation of the RARα consensus motif and the co-occurrence of additional sequence motifs. Together, our genome-wide enhancer activity map elucidates features associated with enhancer activity levels, identifies regulatory regions disregarded by computational prediction tools, and provides a resource for future studies into regulatory elements in mESCs.

Parallel characterization of cis-regulatory elements for multiple genes using CRISPRpath

Current pooled CRISPR screens for cis-regulatory elements (CREs), based on transcriptional output changes, are typically limited to characterizing CREs of only one gene. Here, we describe CRISPRpath, a scalable screening strategy for parallelly characterizing CREs of genes linked to the same biological pathway and converging phenotypes. We demonstrate the ability of CRISPRpath for simultaneously identifying functional enhancers of six genes in the 6-thioguanine–induced DNA mismatch repair pathway using both CRISPR interference (CRISPRi) and CRISPR nuclease (CRISPRn) approaches. Sixty percent of the identified enhancers are known promoters with distinct epigenomic features compared to other active promoters, including increased chromatin accessibility and interactivity. Furthermore, by imposing different levels of selection pressure, CRISPRpath can distinguish enhancers exerting strong impact on gene expression from those exerting weak impact. Our results offer a nuanced view of cis-regulation and demonstrate that CRISPRpath can be leveraged for understanding the complex gene regulatory program beyond transcriptional output at scale.

Discovery of synthetic lethal and tumor suppressor paralog pairs in the human genome

CRISPR screens have accelerated the discovery of important cancer vulnerabilities. However, single-gene knockout phenotypes can be masked by redundancy among related genes. Paralogs constitute two-thirds of the human protein-coding genome, so existing methods are likely inadequate for assaying a large portion of gene function. Here, we develop paired guide RNAs for paralog genetic interaction mapping (pgPEN), a pooled CRISPR-Cas9 single- and double-knockout approach targeting more than 2,000 human paralogs. We apply pgPEN to two cell types and discover that 12% of human paralogs exhibit synthetic lethality in at least one context. We recover known synthetic lethal paralogs MEK1/MEK2, important drug targets CDK4/CDK6, and other synthetic lethal pairs including CCNL1/CCNL2. Additionally, we identify ten tumor suppressor paralog pairs whose compound loss promotes cell proliferation. These findings nominate drug targets and suggest that paralog genetic interactions could shape the landscape of positive and negative selection in cancer.

Androgen signaling connects short isoform production to breakpoint formation at Ewing sarcoma breakpoint region 1

Ewing sarcoma breakpoint region 1 (EWSR1) encodes a multifunctional protein that can cooperate with the transcription factor ERG to promote prostate cancer. The EWSR1 gene is also commonly involved in oncogenic gene rearrangements in Ewing sarcoma. Despite the cancer relevance of EWSR1, its regulation is poorly understood. Here we find that in prostate cancer, androgen signaling upregulates a 5′ EWSR1 isoform by promoting usage of an intronic polyadenylation site. This isoform encodes a cytoplasmic protein that can strongly promote cell migration and clonogenic growth. Deletion of an Androgen Receptor (AR) binding site near the 5′ EWSR1 polyadenylation site abolished androgen-dependent upregulation. This polyadenylation site is also near the Ewing sarcoma breakpoint hotspot, and androgen signaling promoted R-loop and breakpoint formation. RNase H overexpression reduced breakage and 5′ EWSR1 isoform expression suggesting an R-loop dependent mechanism. These data suggest that androgen signaling can promote R-loops internal to the EWSR1 gene leading to either early transcription termination, or breakpoint formation.

Peptide-tiling screens of cancer drivers reveal oncogenic protein domains and associated peptide inhibitors

Gene fragments derived from structural domains mediating physical interactions can modulate biological functions. Utilizing this, we developed lentiviral overexpression libraries of peptides comprehensively tiling high-confidence cancer driver genes. Toward inhibiting cancer growth, we assayed ~66,000 peptides, tiling 65 cancer drivers and 579 mutant alleles. Pooled fitness screens in two breast cancer cell lines revealed peptides, which selectively reduced cellular proliferation, implicating oncogenic protein domains important for cell fitness. Coupling of cell-penetrating motifs to these peptides enabled drug-like function, with peptides derived from EGFR and RAF1 inhibiting cell growth at IC50s of 27-63 μM. We anticipate that this peptide-tiling (PepTile) approach will enable rapid de novo mapping of bioactive protein domains and associated interfering peptides.

Systematic dissection of transcriptional regulatory networks by genome-scale and single-cell CRISPR screens

Millions of putative transcriptional regulatory elements (TREs) have been cataloged in the human genome, yet their functional relevance in specific pathophysiological settings remains to be determined. This is critical to understand how oncogenic transcription factors (TFs) engage specific TREs to impose transcriptional programs underlying malignant phenotypes. Here, we combine cutting edge CRISPR screens and epigenomic profiling to functionally survey ≈15,000 TREs engaged by estrogen receptor (ER). We show that ER exerts its oncogenic role in breast cancer by engaging TREs enriched in GATA3, TFAP2C, and H3K27Ac signal. These TREs control critical downstream TFs, among which TFAP2C plays an essential role in ER-driven cell proliferation. Together, our work reveals novel insights into a critical oncogenic transcription program and provides a framework to map regulatory networks, enabling to dissect the function of the noncoding genome of cancer cells.

Generation of 3'UTR knockout cell lines by CRISPR/Cas9-mediated genome editing

In addition to the protein code, messenger RNAs (mRNAs) also contain untranslated regions (UTRs). 3'UTRs span the region between the translational stop codon and the poly(A) tail. Sequence elements located in 3'UTRs are essential for pre-mRNA processing. 3'UTRs also contain elements that can regulate protein abundance, localization, and function. At least half of all human genes use alternative cleavage and polyadenylation (APA) to further diversify the regulatory potential of protein functions. Traditional gene editing approaches are designed to disrupt the production of functional proteins. Here, we describe a method that allows investigators to manipulate 3'UTR sequences of endogenous genes for both single- 3'UTR and multi-3'UTR genes. As 3'UTRs can regulate individual functions of proteins, techniques to manipulate 3'UTRs at endogenous gene loci will help to disentangle multi-functionality of proteins. Furthermore, the ability to directly examine the impact of gene regulatory elements in 3'UTRs will provide further insights into their functional significance.

Parallel genetics of regulatory sequences using scalable genome editing in vivo

How regulatory sequences control gene expression is fundamental for explaining phenotypes in health and disease. Regulatory elements must ultimately be understood within their genomic environment and development- or tissue-specific contexts. Because this is technically challenging, few regulatory elements have been characterized in vivo. Here, we use inducible Cas9 and multiplexed guide RNAs to create hundreds of mutations in enhancers/promoters and 3' UTRs of 16 genes in C. elegans. Our software crispr-DART analyzes indel mutations in targeted DNA sequencing. We quantify the impact of mutations on expression and fitness by targeted RNA sequencing and DNA sampling. When applying our approach to the lin-41 3' UTR, generating hundreds of mutants, we find that the two adjacent binding sites for the miRNA let-7 can regulate lin-41 expression independently of each other. Finally, we map regulatory genotypes to phenotypic traits for several genes. Our approach enables parallel analysis of regulatory sequences directly in animals.

CRISPR Screens in Synthetic Lethality and Combinatorial Therapies for Cancer

Cancer is a complex disease resulting from the accumulation of genetic dysfunctions. Tumor heterogeneity causes the molecular variety that divergently controls responses to chemotherapy, leading to the recurrent problem of cancer reappearance. For many decades, efforts have focused on identifying essential tumoral genes and cancer driver mutations. More recently, prompted by the clinical success of the synthetic lethality (SL)-based therapy of the PARP inhibitors in homologous recombinant deficient tumors, scientists have centered their novel research on SL interactions (SLI). The state of the art to find new genetic interactions are currently large-scale forward genetic CRISPR screens. CRISPR technology has rapidly evolved to be a common tool in the vast majority of laboratories, as tools to implement CRISPR screen protocols are available to all researchers. Taking advantage of SLI, combinatorial therapies have become the ultimate model to treat cancer with lower toxicity, and therefore better efficiency. This review explores the CRISPR screen methodology, integrates the up-to-date published findings on CRISPR screens in the cancer field and proposes future directions to uncover cancer regulation and individual responses to chemotherapy.

CRISPRi screens reveal a DNA methylation-mediated 3D genome dependent causal mechanism in prostate cancer

Prostate cancer (PCa) risk-associated SNPs are enriched in noncoding cis-regulatory elements (rCREs), yet their modi operandi and clinical impact remain elusive. Here, we perform CRISPRi screens of 260 rCREs in PCa cell lines. We find that rCREs harboring high risk SNPs are more essential for cell proliferation and H3K27ac occupancy is a strong indicator of essentiality. We also show that cell-line-specific essential rCREs are enriched in the 8q24.21 region, with the rs11986220-containing rCRE regulating MYC and PVT1 expression, cell proliferation and tumorigenesis in a cell-line-specific manner, depending on DNA methylation-orchestrated occupancy of a CTCF binding site in between this rCRE and the MYC promoter. We demonstrate that CTCF deposition at this site as measured by DNA methylation level is highly variable in prostate specimens, and observe the MYC eQTL in the 8q24.21 locus in individuals with low CTCF binding. Together our findings highlight a causal mechanism synergistically driven by a risk SNP and DNA methylation-mediated 3D genome architecture, advocating for the integration of genetics and epigenetics in assessing risks conferred by genetic predispositions.

Heart Enhancers: Development and Disease Control at a Distance

Bound by lineage-determining transcription factors and signaling effectors, enhancers play essential roles in controlling spatiotemporal gene expression profiles during development, homeostasis and disease. Recent synergistic advances in functional genomic technologies, combined with the developmental biology toolbox, have resulted in unprecedented genome-wide annotation of heart enhancers and their target genes. Starting with early studies of vertebrate heart enhancers and ending with state-of-the-art genome-wide enhancer discovery and testing, we will review how studying heart enhancers in metazoan species has helped inform our understanding of cardiac development and disease.

In search of lost time: Enhancers as modulators of timing in lymphocyte development and differentiation

Proper timing of gene expression is central to lymphocyte development and differentiation. Lymphocytes often delay gene activation for hours to days after the onset of signaling components, which act on the order of seconds to minutes. Such delays play a prominent role during the intricate choreography of developmental events and during the execution of an effector response. Though a number of mechanisms are sufficient to explain timing at short timescales, it is not known how timing delays are implemented over long timescales that may span several cell generations. Based on the literature, we propose that a class of cis-regulatory elements, termed "timing enhancers," may explain how timing delays are controlled over these long timescales. By considering chromatin as a kinetic barrier to state switching, the timing enhancer model explains experimentally observed dynamics of gene expression where other models fall short. In this review, we elaborate on features of the timing enhancer model and discuss the evidence for its generality throughout development and differentiation. We then discuss potential molecular mechanisms underlying timing enhancer function. Finally, we explore recent evidence drawing connections between timing enhancers and genetic risk for immunopathology. We argue that the timing enhancer model is a useful framework for understanding how cis-regulatory elements control the central dimension of timing in lymphocyte biology.

Enhancing CRISPR deletion via pharmacological delay of DNA-PKcs

CRISPR-Cas9 deletion (CRISPR-del) is the leading approach for eliminating DNA from mammalian cells and underpins a variety of genome-editing applications. Target DNA, defined by a pair of double-strand breaks (DSBs), is removed during nonhomologous end-joining (NHEJ). However, the low efficiency of CRISPR-del results in laborious experiments and false-negative results. By using an endogenous reporter system, we show that repression of the DNA-dependent protein kinase catalytic subunit (DNA-PKcs)—an early step in NHEJ—yields substantial increases in DNA deletion. This is observed across diverse cell lines, gene delivery methods, commercial inhibitors, and guide RNAs, including those that otherwise display negligible activity. We further show that DNA-PKcs inhibition can be used to boost the sensitivity of pooled functional screens and detect true-positive hits that would otherwise be overlooked. Thus, delaying the kinetics of NHEJ relative to DSB formation is a simple and effective means of enhancing CRISPR-deletion.