Single-cell chromatin accessibility reveals malignant regulatory programs in primary human cancers.

To identify cancer-associated gene regulatory changes, we generated single-cell chromatin accessibility landscapes across eight tumor types as part of The Cancer Genome Atlas. Tumor chromatin accessibility is strongly influenced by copy number alterations that can be used to identify subclones, yet underlying cis-regulatory landscapes retain cancer type-specific features. Using organ-matched healthy tissues, we identified the "nearest healthy" cell types in diverse cancers, demonstrating that the chromatin signature of basal-like-subtype breast cancer is most similar to secretory-type luminal epithelial cells. Neural network models trained to learn regulatory programs in cancer revealed enrichment of model-prioritized somatic noncoding mutations near cancer-associated genes, suggesting that dispersed, nonrecurrent, noncoding mutations in cancer are functional. Overall, these data and interpretable gene regulatory models for cancer and healthy tissue provide a framework for understanding cancer-specific gene regulation.

Enhancer activation from transposable elements in extrachromosomal DNA.

Extrachromosomal DNA (ecDNA) is a hallmark of aggressive cancer, contributing to both oncogene amplification and tumor heterogeneity. Here, we used Hi-C, super-resolution imaging, and long-read sequencing to explore the nuclear architecture of MYC-amplified ecDNA in colorectal cancer cells. Intriguingly, we observed frequent spatial proximity between ecDNA and 68 repetitive elements which we called ecDNA-interacting elements or EIEs. To characterize a potential regulatory role of EIEs, we focused on a fragment of the L1M4a1#LINE/L1 which we found to be co-amplified with MYC on ecDNA, gaining enhancer-associated chromatin marks in contrast to its normally silenced state. This EIE, in particular, existed as a naturally occurring structural variant upstream of MYC, gaining oncogenic potential in the transcriptionally permissive ecDNA environment. This EIE sequence is sufficient to enhance MYC expression and is required for cancer cell fitness. These findings suggest that silent repetitive genomic elements can be reactivated on ecDNA, leading to functional cooption and amplification. Repeat element activation on ecDNA represents a mechanism of accelerated evolution and tumor heterogeneity and may have diagnostic and therapeutic potential.

Deciphering regulation of FOXP3 expression in human conventional T cells.

FOXP3 is a lineage-defining transcription factor that controls differentiation and maintenance of suppressive function of regulatory T cells (Tregs). Foxp3 is exclusively expressed in Tregs in mice. However, in humans, FOXP3 is not only constitutively expressed in Tregs; it is also transiently expressed in stimulated CD4+CD25- conventional T cells (Tconvs)1-3. Mechanisms governing the expression of FOXP3 in human Tconvs are not understood. Here, we performed CRISPR interference (CRISPRi) screens using a 15K-member gRNA library tiling 39 kb downstream of the FOXP3 transcriptional start site (TSS) to 85 kb upstream of the TSS in Treg and Tconvs. The FOXP3 promoter and conserved non-coding sequences (CNS0, CNS1, CNS2 and CNS3), characterized as enhancer elements in murine Tregs, were required for maintenance of FOXP3 in human Tregs. In contrast, FOXP3 in human Tconvs depended on regulation at CNS0 and a novel Tconv-specific noncoding sequence (TcNS+) located upstream of CNS0. Arrayed validations of these sites identified an additional repressive cis-element overlapping with the PPP1R3F promoter (TcNS-). Pooled CRISPR knockouts revealed multiple transcription factors required for proper expression of FOXP3 in Tconvs, including GATA3, STAT5, IRF4, ETS1 and DNA methylation-associated regulators DNMT1 and MBD2. Analysis of ChIP-seq and ATAC-seq paired with knock-out (KO) of GATA3, STAT5, IRF4, and ETS1 revealed regulation of CNS0 and TcNS+ accessibility. Collectively, this work identified Treg-shared and Tconv-specific cis-elements and the trans-factors that interact with them, building a network of regulators controlling FOXP3 expression in human Tconvs.

Disparate pathways for extrachromosomal DNA biogenesis and genomic DNA repair.

Oncogene amplification on extrachromosomal DNA (ecDNA) is a pervasive driver event in cancer, yet our understanding of how ecDNA forms is limited. Here, we couple a CRISPR-based method for ecDNA induction with extensive characterization of newly formed ecDNA to examine their biogenesis. We find that DNA circularization is efficient, irrespective of 3D genome context, with formation of 800kb, 1 Mb, and 1.8 Mb ecDNAs reaching or exceeding 15%. We show non-homologous end joining and microhomology-mediated end joining both contribute to ecDNA formation, while inhibition of DNA-PKcs and ATM have opposing impacts on ecDNA formation. EcDNA and the corresponding chromosomal excision scar can form at significantly different rates and respond differently to DNA-PKcs and ATM inhibition. Taken together, our results support a model of ecDNA formation in which double strand break ends dissociate from their legitimate ligation partners prior to joining of illegitimate ends to form the ecDNA and excision scar.

Cohesin prevents cross-domain gene coactivation.

The contrast between the disruption of genome topology after cohesin loss and the lack of downstream gene expression changes instigates intense debates regarding the structure-function relationship between genome and gene regulation. Here, by analyzing transcriptome and chromatin accessibility at the single-cell level, we discover that, instead of dictating population-wide gene expression levels, cohesin supplies a general function to neutralize stochastic coexpression tendencies of cis-linked genes in single cells. Notably, cohesin loss induces widespread gene coactivation and chromatin co-opening tens of million bases apart in cis. Spatial genome and protein imaging reveals that cohesin prevents gene co-bursting along the chromosome and blocks spatial mixing of transcriptional hubs. Single-molecule imaging shows that cohesin confines the exploration of diverse enhancer and core promoter binding transcriptional regulators. Together, these results support that cohesin arranges nuclear topology to control gene coexpression in single cells.

A human autoimmune organoid model reveals IL-7 function in coeliac disease.

In vitro models of autoimmunity are constrained by an inability to culture affected epithelium alongside the complex tissue-resident immune microenvironment. Coeliac disease (CeD) is an autoimmune disease in which dietary gluten-derived peptides bind to the major histocompatibility complex (MHC) class II human leukocyte antigen molecules (HLA)-DQ2 or HLA-DQ8 to initiate immune-mediated duodenal mucosal injury1-4. Here, we generated air-liquid interface (ALI) duodenal organoids from intact fragments of endoscopic biopsies that preserve epithelium alongside native mesenchyme and tissue-resident immune cells as a unit without requiring reconstitution. The immune diversity of ALI organoids spanned T cells, B and plasma cells, natural killer (NK) cells and myeloid cells, with extensive T-cell and B-cell receptor repertoires. HLA-DQ2.5-restricted gluten peptides selectively instigated epithelial destruction in HLA-DQ2.5-expressing organoids derived from CeD patients, and this was antagonized by blocking MHC-II or NKG2C/D. Gluten epitopes stimulated a CeD organoid immune network response in lymphoid and myeloid subsets alongside anti-transglutaminase 2 (TG2) autoantibody production. Functional studies in CeD organoids revealed that interleukin-7 (IL-7) is a gluten-inducible pathogenic modulator that regulates CD8+ T-cell NKG2C/D expression and is necessary and sufficient for epithelial destruction. Furthermore, endogenous IL-7 was markedly upregulated in patient biopsies from active CeD compared with remission disease from gluten-free diets, predominantly in lamina propria mesenchyme. By preserving the epithelium alongside diverse immune populations, this human in vitro CeD model recapitulates gluten-dependent pathology, enables mechanistic investigation and establishes a proof of principle for the organoid modelling of autoimmunity.

Clonal inactivation of TERT impairs stem cell competition.

Telomerase is intimately associated with stem cells and cancer, because it catalytically elongates telomeres-nucleoprotein caps that protect chromosome ends1. Overexpression of telomerase reverse transcriptase (TERT) enhances the proliferation of cells in a telomere-independent manner2-8, but so far, loss-of-function studies have provided no evidence that TERT has a direct role in stem cell function. In many tissues, homeostasis is shaped by stem cell competition, a process in which stem cells compete on the basis of inherent fitness. Here we show that conditional deletion of Tert in the spermatogonial stem cell (SSC)-containing population in mice markedly impairs competitive clone formation. Using lineage tracing from the Tert locus, we find that TERT-expressing SSCs yield long-lived clones, but that clonal inactivation of TERT promotes stem cell differentiation and a genome-wide reduction in open chromatin. This role for TERT in competitive clone formation occurs independently of both its reverse transcriptase activity and the canonical telomerase complex. Inactivation of TERT causes reduced activity of the MYC oncogene, and transgenic expression of MYC in the TERT-deleted pool of SSCs efficiently rescues clone formation. Together, these data reveal a catalytic-activity-independent requirement for TERT in enhancing stem cell competition, uncover a genetic connection between TERT and MYC and suggest that a selective advantage for stem cells with high levels of TERT contributes to telomere elongation in the male germline during homeostasis and ageing.

Transcriptional immune suppression and up-regulation of double-stranded DNA damage and repair repertoires in ecDNA-containing tumors.

Extrachromosomal DNA is a common cause of oncogene amplification in cancer. The non-chromosomal inheritance of ecDNA enables tumors to rapidly evolve, contributing to treatment resistance and poor outcome for patients. The transcriptional context in which ecDNAs arise and progress, including chromosomally-driven transcription, is incompletely understood. We examined gene expression patterns of 870 tumors of varied histological types, to identify transcriptional correlates of ecDNA. Here, we show that ecDNA-containing tumors impact four major biological processes. Specifically, ecDNA-containing tumors up-regulate DNA damage and repair, cell cycle control, and mitotic processes, but down-regulate global immune regulation pathways. Taken together, these results suggest profound alterations in gene regulation in ecDNA-containing tumors, shedding light on molecular processes that give rise to their development and progression.

hnRNPM protects against the dsRNA-mediated interferon response by repressing LINE-associated cryptic splicing.

RNA splicing is pivotal in post-transcriptional gene regulation, yet the exponential expansion of intron length in humans poses a challenge for accurate splicing. Here, we identify hnRNPM as an essential RNA-binding protein that suppresses cryptic splicing through binding to deep introns, maintaining human transcriptome integrity. Long interspersed nuclear elements (LINEs) in introns harbor numerous pseudo splice sites. hnRNPM preferentially binds at intronic LINEs to repress pseudo splice site usage for cryptic splicing. Remarkably, cryptic exons can generate long dsRNAs through base-pairing of inverted ALU transposable elements interspersed among LINEs and consequently trigger an interferon response, a well-known antiviral defense mechanism. Significantly, hnRNPM-deficient tumors show upregulated interferon-associated pathways and elevated immune cell infiltration. These findings unveil hnRNPM as a guardian of transcriptome integrity by repressing cryptic splicing and suggest that targeting hnRNPM in tumors may be used to trigger an inflammatory immune response, thereby boosting cancer surveillance.

Pathways for macrophage uptake of cell-free circular RNAs.

Circular RNAs (circRNAs) are stable RNAs present in cell-free RNA, which may comprise cellular debris and pathogen genomes. Here, we investigate the phenomenon and mechanism of cellular uptake and intracellular fate of exogenous circRNAs. Human myeloid cells and B cells selectively internalize extracellular circRNAs. Macrophage uptake of circRNA is rapid, energy dependent, and saturable. CircRNA uptake can lead to translation of encoded sequences and antigen presentation. The route of internalization influences immune activation after circRNA uptake, with distinct gene expression programs depending on the route of RNA delivery. Genome-scale CRISPR screens and chemical inhibitor studies nominate macrophage scavenger receptor MSR1, Toll-like receptors, and mTOR signaling as key regulators of receptor-mediated phagocytosis of circRNAs, a dominant pathway to internalize circRNAs in parallel to macropinocytosis. These results suggest that cell-free circRNA serves as an "eat me" signal and danger-associated molecular pattern, indicating orderly pathways of recognition and disposal.

Annotation of nuclear lncRNAs based on chromatin interactions.

The human genome is pervasively transcribed and produces a wide variety of long non-coding RNAs (lncRNAs), constituting the majority of transcripts across human cell types. Some specific nuclear lncRNAs have been shown to be important regulatory components acting locally. As RNA-chromatin interaction and Hi-C chromatin conformation data showed that chromatin interactions of nuclear lncRNAs are determined by the local chromatin 3D conformation, we used Hi-C data to identify potential target genes of lncRNAs. RNA-protein interaction data suggested that nuclear lncRNAs act as scaffolds to recruit regulatory proteins to target promoters and enhancers. Nuclear lncRNAs may therefore play a role in directing regulatory factors to locations spatially close to the lncRNA gene. We provide the analysis results through an interactive visualization web portal at https://fantom.gsc.riken.jp/zenbu/reports/#F6_3D_lncRNA.

Organ- and Cell-Selective Delivery of mRNA In Vivo Using Guanidinylated Serinol Charge-Altering Releasable Transporters.

Selective RNA delivery is required for the broad implementation of RNA clinical applications, including prophylactic and therapeutic vaccinations, immunotherapies for cancer, and genome editing. Current polyanion delivery relies heavily on cationic amines, while cationic guanidinium systems have received limited attention due in part to their strong polyanion association, which impedes intracellular polyanion release. Here, we disclose a general solution to this problem in which cationic guanidinium groups are used to form stable RNA complexes upon formulation but at physiological pH undergo a novel charge-neutralization process, resulting in RNA release. This new delivery system consists of guanidinylated serinol moieties incorporated into a charge-altering releasable transporter (GSer-CARTs). Significantly, systematic variations in structure and formulation resulted in GSer-CARTs that exhibit highly selective mRNA delivery to the lung (∼97%) and spleen (∼98%) without targeting ligands. Illustrative of their breadth and translational potential, GSer-CARTs deliver circRNA, providing the basis for a cancer vaccination strategy, which in a murine model resulted in antigen-specific immune responses and effective suppression of established tumors.

Bidirectional epigenetic editing reveals hierarchies in gene regulation.

CRISPR perturbation methods are limited in their ability to study non-coding elements and genetic interactions. In this study, we developed a system for bidirectional epigenetic editing, called CRISPRai, in which we apply activating (CRISPRa) and repressive (CRISPRi) perturbations to two loci simultaneously in the same cell. We developed CRISPRai Perturb-seq by coupling dual perturbation gRNA detection with single-cell RNA sequencing, enabling study of pooled perturbations in a mixed single-cell population. We applied this platform to study the genetic interaction between two hematopoietic lineage transcription factors, SPI1 and GATA1, and discovered novel characteristics of their co-regulation on downstream target genes, including differences in SPI1 and GATA1 occupancy at genes that are regulated through different modes. We also studied the regulatory landscape of IL2 (interleukin-2) in Jurkat T cells, primary T cells and chimeric antigen receptor (CAR) T cells and elucidated mechanisms of enhancer-mediated IL2 gene regulation. CRISPRai facilitates investigation of context-specific genetic interactions, provides new insights into gene regulation and will enable exploration of non-coding disease-associated variants.

The CD8+ T cell tolerance checkpoint triggers a distinct differentiation state defined by protein translation defects.

Peripheral CD8+ T cell tolerance is a checkpoint in both autoimmune disease and anti-cancer immunity. Despite its importance, the relationship between tolerance-induced states and other CD8+ T cell differentiation states remains unclear. Using flow cytometric phenotyping, single-cell RNA sequencing (scRNA-seq), and chromatin accessibility profiling, we demonstrated that in vivo peripheral tolerance to a self-antigen triggered a fundamentally distinct differentiation state separate from exhaustion, memory, and functional effector cells but analogous to cells defectively primed against tumors. Tolerant cells diverged early and progressively from effector cells, adopting a transcriptionally and epigenetically distinct state within 60 h of antigen encounter. Breaching tolerance required the synergistic actions of strong T cell receptor (TCR) signaling and inflammation, which cooperatively induced gene modules that enhanced protein translation. Weak TCR signaling during bystander infection failed to breach tolerance due to the uncoupling of effector gene expression from protein translation. Thus, tolerance engages a distinct differentiation trajectory enforced by protein translation defects.

Escape from X inactivation is directly modulated by levels of Xist non-coding RNA.

In placental females, one copy of the two X chromosomes is largely silenced during a narrow developmental time window, in a process mediated by the non-coding RNA Xist1. Here, we demonstrate that Xist can initiate X-chromosome inactivation (XCI) well beyond early embryogenesis. By modifying its endogenous level, we show that Xist has the capacity to actively silence genes that escape XCI both in neuronal progenitor cells (NPCs) and in vivo, in mouse embryos. We also show that Xist plays a direct role in eliminating TAD-like structures associated with clusters of escapee genes on the inactive X chromosome, and that this is dependent on Xist's XCI initiation partner, SPEN2. We further demonstrate that Xist's function in suppressing gene expression of escapees and topological domain formation is reversible for up to seven days post-induction, but that sustained Xist up-regulation leads to progressively irreversible silencing and CpG island DNA methylation of facultative escapees. Thus, the distinctive transcriptional and regulatory topologies of the silenced X chromosome is actively, directly - and reversibly - controlled by Xist RNA throughout life.

Epigenetic signature and key transcriptional regulators of human antigen-specific type 1 regulatory T cells.

Human adaptive immunity is orchestrated by effector and regulatory T (Treg) cells. Natural Tregs arise in the thymus where they are shaped to recognize self-antigens, while type 1 Tregs or Tr1 cells are induced from conventional peripheral CD4 + T cells in response to peripheral antigens, such as alloantigens and allergens. Tr1 cells have been developed as a potential therapy for inducing antigen-specific tolerance, because they can be rapidly differentiated in vitro in response to a target antigen. However, the epigenetic landscape and the identity of transcription factors (TFs) that regulate differentiation, phenotype, and functions of human antigen-specific Tr1 cells is largely unknown, hindering Tr1 research and broader clinical development. Here, we reveal the unique epigenetic signature of antigen-specific Tr1 cells, and TFs that regulate their differentiation, phenotype and function. We showed that in vitro induced antigen-specific Tr1 cells are distinct both clonally and transcriptionally from natural Tregs and other conventional CD4 + T cells on a single-cell level. An integrative analysis of Tr1 cell epigenome and transcriptome identified a TF signature unique to antigen-specific Tr1 cells, and predicted that IRF4, BATF, and MAF act as their transcriptional regulators. Using functional genomics, we showed that each of these TFs play a non-redundant role in regulating Tr1 cell differentiation, suppressive function, and expression of co-inhibitory and cytotoxic proteins. By using the Tr1-specific TF signature as a molecular fingerprint, we tracked Tr1 cells in peripheral blood of recipients of allogeneic hematopoietic stem cell transplantation treated with adoptive Tr1 cell therapy. Furthermore, the same signature identified Tr1 cells in resident CD4 + T cells in solid tumors. Altogether, these results reveal the epigenetic signature and the key transcriptional regulators of human Tr1 cells. These data will guide mechanistic studies of human Tr1 cell biology and the development and optimization of adoptive Tr1 cell therapies.

Regulation of immune signal integration and memory by inflammation-induced chromosome conformation.

3-dimensional (3D) genome conformation is central to gene expression regulation, yet our understanding of its contribution to rapid transcriptional responses, signal integration, and memory in immune cells is limited. Here, we study the molecular regulation of the inflammatory response in primary macrophages using integrated transcriptomic, epigenomic, and chromosome conformation data, including base pair-resolution Micro-Capture C. We demonstrate that interleukin-4 (IL-4) primes the inflammatory response in macrophages by stably rewiring 3D genome conformation, juxtaposing endotoxin-, interferon-gamma-, and dexamethasone-responsive enhancers in close proximity to their cognate gene promoters. CRISPR-based perturbations of enhancer-promoter contacts or CCCTC-binding factor (CTCF) boundary elements demonstrated that IL-4-driven conformation changes are indispensable for enhanced and synergistic endotoxin-induced transcriptional responses, as well as transcriptional memory following stimulus removal. Moreover, transcriptional memory mediated by changes in chromosome conformation often occurred in the absence of changes in chromatin accessibility or histone modifications. Collectively, these findings demonstrate that rapid and memory transcriptional responses to immunological stimuli are encoded in the 3D genome.

CoRAL accurately resolves extrachromosomal DNA genome structures with long-read sequencing.

Extrachromosomal DNA (ecDNA) is a central mechanism for focal oncogene amplification in cancer, occurring in approximately 15% of early stage cancers and 30% of late-stage cancers. EcDNAs drive tumor formation, evolution, and drug resistance by dynamically modulating oncogene copy-number and rewiring gene-regulatory networks. Elucidating the genomic architecture of ecDNA amplifications is critical for understanding tumor pathology and developing more effective therapies. Paired-end short-read (Illumina) sequencing and mapping have been utilized to represent ecDNA amplifications using a breakpoint graph, where the inferred architecture of ecDNA is encoded as a cycle in the graph. Traversals of breakpoint graph have been used to successfully predict ecDNA presence in cancer samples. However, short-read technologies are intrinsically limited in the identification of breakpoints, phasing together of complex rearrangements and internal duplications, and deconvolution of cell-to-cell heterogeneity of ecDNA structures. Long-read technologies, such as from Oxford Nanopore Technologies, have the potential to improve inference as the longer reads are better at mapping structural variants and are more likely to span rearranged or duplicated regions. Here, we propose CoRAL (Complete Reconstruction of Amplifications with Long reads), for reconstructing ecDNA architectures using long-read data. CoRAL reconstructs likely cyclic architectures using quadratic programming that simultaneously optimizes parsimony of reconstruction, explained copy number, and consistency of long-read mapping. CoRAL substantially improves reconstructions in extensive simulations and 9 datasets from previously-characterized cell-lines as compared to previous short-read-based tools. As long-read usage becomes wide-spread, we anticipate that CoRAL will be a valuable tool for profiling the landscape and evolution of focal amplifications in tumors.

Approaches to probe and perturb long noncoding RNA functions in diseases.

Long noncoding RNAs (lncRNAs) are a class of RNA molecules exceeding 200 nucleotides in length that lack long open-reading frames. Transcribed predominantly by RNA polymerase II (>500nt), lncRNAs can undergo splicing and are produced from various regions of the genome, including intergenic regions, introns, and in antisense orientation to protein-coding genes. Aberrations in lncRNA expression or function have been associated with a wide variety of diseases, including cancer, cardiovascular diseases, diabetes, and neurodegeneration. Despite the growing recognition of select lncRNAs as key players in cellular processes and diseases, several challenges obscure a comprehensive understanding of their functional landscape. Recent technological innovations, such as in sequencing, affinity-based techniques, imaging, and RNA perturbation, have advanced functional characterization and mechanistic understanding of disease-associated lncRNAs.

Xist ribonucleoproteins promote female sex-biased autoimmunity.

Autoimmune diseases disproportionately affect females more than males. The XX sex chromosome complement is strongly associated with susceptibility to autoimmunity. Xist long non-coding RNA (lncRNA) is expressed only in females to randomly inactivate one of the two X chromosomes to achieve gene dosage compensation. Here, we show that the Xist ribonucleoprotein (RNP) complex comprising numerous autoantigenic components is an important driver of sex-biased autoimmunity. Inducible transgenic expression of a non-silencing form of Xist in male mice introduced Xist RNP complexes and sufficed to produce autoantibodies. Male SJL/J mice expressing transgenic Xist developed more severe multi-organ pathology in a pristane-induced lupus model than wild-type males. Xist expression in males reprogrammed T and B cell populations and chromatin states to more resemble wild-type females. Human patients with autoimmune diseases displayed significant autoantibodies to multiple components of XIST RNP. Thus, a sex-specific lncRNA scaffolds ubiquitous RNP components to drive sex-biased immunity.