A CTCF-binding site in the Mdm1-Il22-Ifng locus shapes cytokine expression profiles and plays a critical role in early Th1 cell fate specification.

Cytokine expression during T cell differentiation is a highly regulated process that involves long-range promoter-enhancer and CTCF-CTCF contacts at cytokine loci. Here, we investigated the impact of dynamic chromatin loop formation within the topologically associating domain (TAD) in regulating the expression of interferon gamma (IFN-γ) and interleukin-22 (IL-22); these cytokine loci are closely located in the genome and are associated with complex enhancer landscapes, which are selectively active in type 1 and type 3 lymphocytes. In situ Hi-C analyses revealed inducible TADs that insulated Ifng and Il22 enhancers during Th1 cell differentiation. Targeted deletion of a 17 bp boundary motif of these TADs imbalanced Th1- and Th17-associated immunity, both in vitro and in vivo, upon Toxoplasma gondii infection. In contrast, this boundary element was dispensable for cytokine regulation in natural killer cells. Our findings suggest that precise cytokine regulation relies on lineage- and developmental stage-specific interactions of 3D chromatin architectures and enhancer landscapes.

Human lung cancer harbors spatially organized stem-immunity hubs associated with response to immunotherapy.

The organization of immune cells in human tumors is not well understood. Immunogenic tumors harbor spatially localized multicellular 'immunity hubs' defined by expression of the T cell-attracting chemokines CXCL10/CXCL11 and abundant T cells. Here, we examined immunity hubs in human pre-immunotherapy lung cancer specimens and found an association with beneficial response to PD-1 blockade. Critically, we discovered the stem-immunity hub, a subtype of immunity hub strongly associated with favorable PD-1-blockade outcome. This hub is distinct from mature tertiary lymphoid structures and is enriched for stem-like TCF7+PD-1+CD8+ T cells, activated CCR7+LAMP3+ dendritic cells and CCL19+ fibroblasts as well as chemokines that organize these cells. Within the stem-immunity hub, we find preferential interactions between CXCL10+ macrophages and TCF7-CD8+ T cells as well as between mature regulatory dendritic cells and TCF7+CD4+ and regulatory T cells. These results provide a picture of the spatial organization of the human intratumoral immune response and its relevance to patient immunotherapy outcomes.

Dysregulated Repeat Element Viral-like Immune Response in Hepatocellular Carcinoma.

Purpose: Dysregulation of viral-like repeat RNAs are a common feature across many malignancies that are linked with immunological response, but the characterization of these in hepatocellular carcinoma (HCC) is understudied. In this study, we performed RNA in situ hybridization (RNA-ISH) of different repeat RNAs, immunohistochemistry (IHC) for immune cell subpopulations, and spatial transcriptomics to understand the relationship of HCC repeat expression, immune response, and clinical outcomes. Experimental Design: RNA-ISH for LINE1, HERV-K, HERV-H, and HSATII repeats and IHC for T-cell, Treg, B-cell, macrophage, and immune checkpoint markers were performed on 43 resected HCC specimens. Spatial transcriptomics on tumor and vessel regions of interest was performed on 28 specimens from the same cohort. Results: High HERV-K and high LINE1 expression were both associated with worse overall survival. There was a positive correlation between LINE1 expression and FOXP3 T-regulatory cells (r = 0.51 p < 0.001) as well as expression of the TIM3 immune checkpoint (r = 0.34, p = 0.03). Spatial transcriptomic profiling of HERV-K high and LINE-1 high tumors identified elevated expression of multiple genes previously associated with epithelial mesenchymal transition, cellular proliferation, and worse overall prognosis in HCC including SSX1, MAGEC2, and SPINK1. Conclusion: Repeat RNAs may serve as useful prognostic biomarkers in HCC and may also serve as novel therapeutic targets. Additional study is needed to understand the mechanisms by which repeat RNAs impact HCC tumorigenesis.

High-resolution CTCF footprinting reveals impact of chromatin state on cohesin extrusion dynamics.

DNA looping is vital for establishing many enhancer-promoter interactions. While CTCF is known to anchor many cohesin-mediated loops, the looped chromatin fiber appears to predominantly exist in a poorly characterized actively extruding state. To better characterize extruding chromatin loop structures, we used CTCF MNase HiChIP data to determine both CTCF binding at high resolution and 3D contact information. Here we present FactorFinder, a tool that identifies CTCF binding sites at near base-pair resolution. We leverage this substantial advance in resolution to determine that the fully extruded (CTCF-CTCF) state is rare genome-wide with locus-specific variation from ~1-10%. We further investigate the impact of chromatin state on loop extrusion dynamics, and find that active enhancers and RNA Pol II impede cohesin extrusion, facilitating an enrichment of enhancer-promoter contacts in the partially extruded loop state. We propose a model of topological regulation whereby the transient, partially extruded states play active roles in transcription.

Integrative PTEN Enhancer Discovery Reveals a New Model of Enhancer Organization.

Enhancers possess both structural elements mediating promoter looping and functional elements mediating gene expression. Traditional models of enhancer-mediated gene regulation imply genomic overlap or immediate adjacency of these elements. We test this model by combining densely-tiled CRISPRa screening with nucleosome-resolution Region Capture Micro-C topology analysis. Using this integrated approach, we comprehensively define the cis-regulatory landscape for the tumor suppressor PTEN, identifying and validating 10 distinct enhancers and defining their 3D spatial organization. Unexpectedly, we identify several long-range functional enhancers whose promoter proximity is facilitated by chromatin loop anchors several kilobases away, and demonstrate that accounting for this spatial separation improves the computational prediction of validated enhancers. Thus, we propose a new model of enhancer organization incorporating spatial separation of essential functional and structural components.

Spatial transcriptomics reveals distinct tissue niches linked with steroid responsiveness in acute gastrointestinal GVHD.

Severe acute graft-versus-host disease (aGVHD) is associated with significant mortality and morbidity, especially in steroid-resistant (SR) cases. Spatial transcriptomic technology can elucidate tissue-based interactions in vivo and possibly identify predictors of treatment response. Tissue sections from 32 treatment-naïve patients with biopsy-confirmed lower gastrointestinal (GI) aGVHD were obtained. The GeoMx digital spatial profiler was used to capture transcriptome profiles of >18 000 genes from different foci of immune infiltrates, colonic epithelium, and vascular endothelium. Each tissue compartment sampled showed 2 distinct clusters that were analyzed for differential expression and spatially resolved correlation of gene signatures. Classic cell-mediated immunity signatures, normal differentiated epithelial cells, and inflamed vasculature dominated foci sampled from steroid-sensitive cases. In contrast, a neutrophil predominant noncanonical inflammation with regenerative epithelial cells and some indication of angiogenic endothelial response was overrepresented in areas from SR cases. Evaluation of potential prognostic biomarkers identified ubiquitin specific peptidase 17-like (USP17L) family of genes as being differentially expressed in immune cells from patients with worsened survival. In summary, we demonstrate distinct tissue niches with unique gene expression signatures within lower GI tissue from patients with aGVHD and provide evidence of a potential prognostic biomarker.

DNA hypomethylation silences anti-tumor immune genes in early prostate cancer and CTCs.

Cancer is characterized by hypomethylation-associated silencing of large chromatin domains, whose contribution to tumorigenesis is uncertain. Through high-resolution genome-wide single-cell DNA methylation sequencing, we identify 40 core domains that are uniformly hypomethylated from the earliest detectable stages of prostate malignancy through metastatic circulating tumor cells (CTCs). Nested among these repressive domains are smaller loci with preserved methylation that escape silencing and are enriched for cell proliferation genes. Transcriptionally silenced genes within the core hypomethylated domains are enriched for immune-related genes; prominent among these is a single gene cluster harboring all five CD1 genes that present lipid antigens to NKT cells and four IFI16-related interferon-inducible genes implicated in innate immunity. The re-expression of CD1 or IFI16 murine orthologs in immuno-competent mice abrogates tumorigenesis, accompanied by the activation of anti-tumor immunity. Thus, early epigenetic changes may shape tumorigenesis, targeting co-located genes within defined chromosomal loci. Hypomethylation domains are detectable in blood specimens enriched for CTCs.

Epigenetic clocks, aging, and cancer.

Global methylation changes in aging cells affect cancer risk and tissue homeostasis.

Polycomb-lamina antagonism partitions heterochromatin at the nuclear periphery.

The genome can be divided into two spatially segregated compartments, A and B, which partition active and inactive chromatin states. While constitutive heterochromatin is predominantly located within the B compartment near the nuclear lamina, facultative heterochromatin marked by H3K27me3 spans both compartments. How epigenetic modifications, compartmentalization, and lamina association collectively maintain heterochromatin architecture remains unclear. Here we develop Lamina-Inducible Methylation and Hi-C (LIMe-Hi-C) to jointly measure chromosome conformation, DNA methylation, and lamina positioning. Through LIMe-Hi-C, we identify topologically distinct sub-compartments with high levels of H3K27me3 and differing degrees of lamina association. Inhibition of Polycomb repressive complex 2 (PRC2) reveals that H3K27me3 is essential for sub-compartment segregation. Unexpectedly, PRC2 inhibition promotes lamina association and constitutive heterochromatin spreading into H3K27me3-marked B sub-compartment regions. Consistent with this repositioning, genes originally marked with H3K27me3 in the B compartment, but not the A compartment, remain largely repressed, suggesting that constitutive heterochromatin spreading can compensate for H3K27me3 loss at a transcriptional level. These findings demonstrate that Polycomb sub-compartments and their antagonism with lamina association are fundamental features of genome structure. More broadly, by jointly measuring nuclear position and Hi-C contacts, our study demonstrates how compartmentalization and lamina association represent distinct but interdependent modes of heterochromatin regulation.

CRISPR-Cas9 treatment partially restores amyloid-ß 42/40 in human fibroblasts with the Alzheimer's disease PSEN 1 M146L mutation.

Presenilin 1 (PS1) is a central component of γ-secretase, an enzymatic complex involved in the generation of the amyloid-ß (Aß) peptide that deposits as plaques in the Alzheimer's disease (AD) brain. The M146L mutation in the PS1 gene (PSEN1) leads to an autosomal dominant form of early-onset AD by promoting a relative increase in the generation of the more aggregation-prone Aß42. This change is evident not only in the brain but also in peripheral cells of mutation carriers. In this study we used the CRISPR-Cas9 system from Streptococcus pyogenes to selectively disrupt the PSEN1 M146L allele in human fibroblasts. A disruption of more than 50% of mutant alleles was observed in all CRISPR-Cas9-treated samples, resulting in reduced extracellular Aß42/40 ratios. Fluorescence resonance energy transfer-based conformation and western blot analyses indicated that CRISPR-Cas9 treatment also affects the overall PS1 conformation and reduces PS1 levels. Moreover, our guide RNA did not lead to any detectable editing at the highest-ranking candidate off-target sites identified by ONE-seq and CIRCLE-seq. Overall, our data support the effectiveness of CRISPR-Cas9 in selectively targeting the PSEN1 M146L allele and counteracting the AD-associated phenotype. We believe that this system could be developed into a therapeutic strategy for patients with this and other dominant mutations leading to early-onset AD.

Reverse Transcriptase Inhibition Disrupts Repeat Element Life Cycle in Colorectal Cancer.

Altered RNA expression of repetitive sequences and retrotransposition are frequently seen in colorectal cancer, implicating a functional importance of repeat activity in cancer progression. We show the nucleoside reverse transcriptase inhibitor 3TC targets activities of these repeat elements in colorectal cancer preclinical models with a preferential effect in p53-mutant cell lines linked with direct binding of p53 to repeat elements. We translate these findings to a human phase II trial of single-agent 3TC treatment in metastatic colorectal cancer with demonstration of clinical benefit in 9 of 32 patients. Analysis of 3TC effects on colorectal cancer tumorspheres demonstrates accumulation of immunogenic RNA:DNA hybrids linked with induction of interferon response genes and DNA damage response. Epigenetic and DNA-damaging agents induce repeat RNAs and have enhanced cytotoxicity with 3TC. These findings identify a vulnerability in colorectal cancer by targeting the viral mimicry of repeat elements. SIGNIFICANCE: Colorectal cancers express abundant repeat elements that have a viral-like life cycle that can be therapeutically targeted with nucleoside reverse transcriptase inhibitors (NRTI) commonly used for viral diseases. NRTIs induce DNA damage and interferon response that provide a new anticancer therapeutic strategy. This article is highlighted in the In This Issue feature, p. 1397.

Genome-wide DNA methylation patterns reveal clinically relevant predictive and prognostic subtypes in human osteosarcoma.

Aberrant methylation of genomic DNA has been reported in many cancers. Specific DNA methylation patterns have been shown to provide clinically useful prognostic information and define molecular disease subtypes with different response to therapy and long-term outcome. Osteosarcoma is an aggressive malignancy for which approximately half of tumors recur following standard combined surgical resection and chemotherapy. No accepted prognostic factor save tumor necrosis in response to adjuvant therapy currently exists, and traditional genomic studies have thus far failed to identify meaningful clinical associations. We studied the genome-wide methylation state of primary tumors and tested how they predict patient outcomes. We discovered relative genomic hypomethylation to be strongly predictive of response to standard chemotherapy. Recurrence and survival were also associated with genomic methylation, but through more site-specific patterns. Furthermore, the methylation patterns were reproducible in three small independent clinical datasets. Downstream transcriptional, in vitro, and pharmacogenomic analysis provides insight into the clinical translation of the methylation patterns. Our findings suggest the assessment of genomic methylation may represent a strategy for stratifying patients for the application of alternative therapies.

Congenital anemia reveals distinct targeting mechanisms for master transcription factor GATA1.

Master regulators, such as the hematopoietic transcription factor (TF) GATA1, play an essential role in orchestrating lineage commitment and differentiation. However, the precise mechanisms by which such TFs regulate transcription through interactions with specific cis-regulatory elements remain incompletely understood. Here, we describe a form of congenital hemolytic anemia caused by missense mutations in an intrinsically disordered region of GATA1, with a poorly understood role in transcriptional regulation. Through integrative functional approaches, we demonstrate that these mutations perturb GATA1 transcriptional activity by partially impairing nuclear localization and selectively altering precise chromatin occupancy by GATA1. These alterations in chromatin occupancy and concordant chromatin accessibility changes alter faithful gene expression, with failure to both effectively silence and activate select genes necessary for effective terminal red cell production. We demonstrate how disease-causing mutations can reveal regulatory mechanisms that enable the faithful genomic targeting of master TFs during cellular differentiation.

Defining genome-wide CRISPR-Cas genome-editing nuclease activity with GUIDE-seq.

Genome-wide unbiased identification of double-stranded breaks enabled by sequencing (GUIDE-seq) is a sensitive, unbiased, genome-wide method for defining the activity of genome-editing nucleases in living cells. GUIDE-seq is based on the principle of efficient integration of an end-protected double-stranded oligodeoxynucleotide tag into sites of nuclease-induced DNA double-stranded breaks, followed by amplification of tag-containing genomic DNA molecules and high-throughput sequencing. Here we describe a detailed GUIDE-seq protocol including cell transfection, library preparation, sequencing and bioinformatic analysis. The entire protocol including cell culture can be completed in 9 d. Once tag-integrated genomic DNA is isolated, library preparation, sequencing and analysis can be performed in 3 d. The result is a genome-wide catalog of off-target sites ranked by nuclease activity as measured by GUIDE-seq read counts. GUIDE-seq is one of the most sensitive cell-based methods for defining genome-wide off-target activity and has been broadly adopted for research and therapeutic use.

Augmenting and directing long-range CRISPR-mediated activation in human cells.

Epigenetic editing is an emerging technology that uses artificial transcription factors (aTFs) to regulate expression of a target gene. Although human genes can be robustly upregulated by targeting aTFs to promoters, the activation induced by directing aTFs to distal transcriptional enhancers is substantially less robust and consistent. Here we show that long-range activation using CRISPR-based aTFs in human cells can be made more efficient and reliable by concurrently targeting an aTF to the target gene promoter. We used this strategy to direct target gene choice for enhancers capable of regulating more than one promoter and to achieve allele-selective activation of human genes by targeting aTFs to single-nucleotide polymorphisms embedded in distally located sequences. Our results broaden the potential applications of the epigenetic editing toolbox for research and therapeutics.

Smart-RRBS for single-cell methylome and transcriptome analysis.

The integration of DNA methylation and transcriptional state within single cells is of broad interest. Several single-cell dual- and multi-omics approaches have been reported that enable further investigation into cellular heterogeneity, including the discovery and in-depth study of rare cell populations. Such analyses will continue to provide important mechanistic insights into the regulatory consequences of epigenetic modifications. We recently reported a new method for profiling the DNA methylome and transcriptome from the same single cells in a cancer research study. Here, we present details of the protocol and provide guidance on its utility. Our Smart-RRBS (reduced representation bisulfite sequencing) protocol combines Smart-seq2 and RRBS and entails physically separating mRNA from the genomic DNA. It generates paired epigenetic promoter and RNA-expression measurements for ~24% of protein-coding genes in a typical single cell. It also works for micro-dissected tissue samples comprising hundreds of cells. The protocol, excluding flow sorting of cells and sequencing, takes ~3 d to process up to 192 samples manually. It requires basic molecular biology expertise and laboratory equipment, including a PCR workstation with UV sterilization, a DNA fluorometer and a microfluidic electrophoresis system.

STAG2 loss rewires oncogenic and developmental programs to promote metastasis in Ewing sarcoma.

The core cohesin subunit STAG2 is recurrently mutated in Ewing sarcoma but its biological role is less clear. Here, we demonstrate that cohesin complexes containing STAG2 occupy enhancer and polycomb repressive complex (PRC2)-marked regulatory regions. Genetic suppression of STAG2 leads to a compensatory increase in cohesin-STAG1 complexes, but not in enhancer-rich regions, and results in reprogramming of cis-chromatin interactions. Strikingly, in STAG2 knockout cells the oncogenic genetic program driven by the fusion transcription factor EWS/FLI1 was highly perturbed, in part due to altered enhancer-promoter contacts. Moreover, loss of STAG2 also disrupted PRC2-mediated regulation of gene expression. Combined, these transcriptional changes converged to modulate EWS/FLI1, migratory, and neurodevelopmental programs. Finally, consistent with clinical observations, functional studies revealed that loss of STAG2 enhances the metastatic potential of Ewing sarcoma xenografts. Our findings demonstrate that STAG2 mutations can alter chromatin architecture and transcriptional programs to promote an aggressive cancer phenotype.

Extended-representation bisulfite sequencing of gene regulatory elements in multiplexed samples and single cells.

The biological roles of DNA methylation have been elucidated by profiling methods based on whole-genome or reduced-representation bisulfite sequencing, but these approaches do not efficiently survey the vast numbers of non-coding regulatory elements in mammalian genomes. Here we present an extended-representation bisulfite sequencing (XRBS) method for targeted profiling of DNA methylation. Our design strikes a balance between expanding coverage of regulatory elements and reproducibly enriching informative CpG dinucleotides in promoters, enhancers and CTCF binding sites. Barcoded DNA fragments are pooled before bisulfite conversion, allowing multiplex processing and technical consistency in low-input samples. Application of XRBS to single leukemia cells enabled us to evaluate genetic copy number variations and methylation variability across individual cells. Our analysis highlights heterochromatic H3K9me3 regions as having the highest cell-to-cell variability in their methylation, likely reflecting inherent epigenetic instability of these late-replicating regions, compounded by differences in cell cycle stages among sampled cells.