Search

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.

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., (© 2024. The Author(s).)

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., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Agrawal et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Background: Current evidence suggests that cis-regulatory elements controlling gene expression may be the predominant target of natural selection in humans and other species. Detecting selection acting on these elements is critical to understanding evolution but remains challenging because we do not know which mutations will affect gene regulation., Results: To address this, we devise an approach to search for lineage-specific selection on three critical steps in transcriptional regulation: chromatin activity, transcription factor binding, and chromosomal looping. Applying this approach to lymphoblastoid cells from 831 individuals of either European or African descent, we find strong signals of differential chromatin activity linked to gene expression differences between ancestries in numerous contexts, but no evidence of functional differences in chromosomal looping. Moreover, we show that enhancers rather than promoters display the strongest signs of selection associated with sites of differential transcription factor binding., Conclusions: Overall, our study indicates that some cis-regulatory adaptation may be more easily detected at the level of chromatin than DNA sequence. This work provides a vast resource of genomic interaction data from diverse human populations and establishes a novel selection test that will benefit future study of regulatory evolution in humans and other species., (© 2024. The Author(s).)

Cell cycle (CC) facilitates cell division via robust, cyclical gene expression. Protective immunity requires the expansion of pathogen-responsive cell types, but whether CC confers unique gene expression programs that direct the subsequent immunological response remains unclear. Here, we demonstrate that single macrophages (MFs) adopt different plasticity states in CC, which leads to heterogeneous cytokine-induced polarization, priming, and repolarization programs. Specifically, MF plasticity to interferon gamma (IFNG) is substantially reduced during S-G2/M, whereas interleukin 4 (IL-4) induces S-G2/M-biased gene expression, mediated by CC-biased enhancers. Additionally, IL-4 polarization shifts the CC-phase distribution of MFs toward the G2/M phase, providing a subpopulation-specific mechanism for IL-4-induced, dampened IFNG responsiveness. Finally, we demonstrate CC-dependent MF responses in murine and human disease settings in vivo, including Th2-driven airway inflammation and pulmonary fibrosis, where MFs express an S-G2/M-biased tissue remodeling gene program. Therefore, MF inflammatory and regenerative responses are gated by CC in a cyclical, phase-dependent manner., Competing Interests: Declaration of interests J.A.B. is a consultant for Immunai. A.T.S. is a founder of Immunai and Cartography Biosciences and receives research funding from Allogene Therapeutics and Merck Research Laboratories. H.Y.C. is a member of the Advisory Board of Molecular Cell, is a co-founder of Accent Therapeutics, Boundless Bio, and Cartography Biosciences, and is an advisor to 10x Genomics, Arsenal Biosciences, and Spring Discovery., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Physiologic laminar shear stress (LSS) induces an endothelial gene expression profile that is vasculo-protective. In this report, we delineate how LSS mediates changes in the epigenetic landscape to promote this beneficial response. We show that under LSS, KLF4 interacts with the SWI/SNF nucleosome remodeling complex to increase accessibility at enhancer sites that promote the expression of homeostatic endothelial genes. By combining molecular and computational approaches we discover enhancers that loop to promoters of KLF4- and LSS-responsive genes that stabilize endothelial cells and suppress inflammation, such as BMPR2, SMAD5, and DUSP5. By linking enhancers to genes that they regulate under physiologic LSS, our work establishes a foundation for interpreting how non-coding DNA variants in these regions might disrupt protective gene expression to influence vascular disease., (© 2022. The Author(s).)

Although alterations in chromatin structure are known to exist in tumors, how these alterations relate to molecular phenotypes in cancer remains to be demonstrated. Multi-omics profiling of human tumors can provide insight into how alterations in chromatin structure are propagated through the pathway of gene expression to result in malignant protein expression. We applied multi-omics profiling of chromatin accessibility, RNA abundance, and protein abundance to 36 human thyroid cancer primary tumors, metastases, and patient-match normal tissue. Through quantification of chromatin accessibility associated with active transcription units and global protein expression, we identify a local chromatin structure that is highly correlated with coordinated RNA and protein expression. In particular, we identify enhancers located within gene-bodies as predictive of correlated RNA and protein expression, that is independent of overall transcriptional activity. To demonstrate the generalizability of these findings we also identify similar results in an independent cohort of human breast cancers. Taken together, these analyses suggest that local enhancers, rather than distal enhancers, are likely most predictive of cancer gene expression phenotypes. This allows for identification of potential targets for cancer therapeutic approaches and reinforces the utility of multi-omics profiling as a methodology to understand human disease., (© 2021. The Author(s).)

Genetic perturbations of cortical development can lead to neurodevelopmental disease, including autism spectrum disorder (ASD). To identify genomic regions crucial to corticogenesis, we mapped the activity of gene-regulatory elements generating a single-cell atlas of gene expression and chromatin accessibility both independently and jointly. This revealed waves of gene regulation by key transcription factors (TFs) across a nearly continuous differentiation trajectory, distinguished the expression programs of glial lineages, and identified lineage-determining TFs that exhibited strong correlation between linked gene-regulatory elements and expression levels. These highly connected genes adopted an active chromatin state in early differentiating cells, consistent with lineage commitment. Base-pair-resolution neural network models identified strong cell-type-specific enrichment of noncoding mutations predicted to be disruptive in a cohort of ASD individuals and identified frequently disrupted TF binding sites. This approach illustrates how cell-type-specific mapping can provide insights into the programs governing human development and disease., Competing Interests: Declaration of interests W.J.G. was a consultant for 10x Genomics and is named as an inventor on patents describing ATAC-seq methods. H.Y.C. is a co-founder of Accent Therapeutics and Boundless Bio and an advisor of 10x Genomics, Arsenal Biosciences, and Spring Discovery. A.S. is an employee of Insitro, Inc, and receives consulting fees from Myokardia, Inc. K.F. is an employee of Illumina, Inc., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Metastasis is the leading cause of cancer-related deaths and enables cancer cells to compromise organ function by expanding in secondary sites. Since primary tumours and metastases often share the same constellation of driver mutations, the mechanisms that drive their distinct phenotypes are unclear. Here we show that inactivation of the frequently mutated tumour suppressor gene LKB1 (encoding liver kinase B1) has evolving effects throughout the progression of lung cancer, which leads to the differential epigenetic re-programming of early-stage primary tumours compared with late-stage metastases. By integrating genome-scale CRISPR-Cas9 screening with bulk and single-cell multi-omic analyses, we unexpectedly identify LKB1 as a master regulator of chromatin accessibility in lung adenocarcinoma primary tumours. Using an in vivo model of metastatic progression, we further show that loss of LKB1 activates the early endoderm transcription factor SOX17 in metastases and a metastatic-like sub-population of cancer cells within primary tumours. The expression of SOX17 is necessary and sufficient to drive a second wave of epigenetic changes in LKB1-deficient cells that enhances metastatic ability. Overall, our study demonstrates how the downstream effects of an individual driver mutation can change throughout cancer development, with implications for stage-specific therapeutic resistance mechanisms and the gene regulatory underpinnings of metastatic evolution., (© 2021. The Author(s), under exclusive licence to Springer Nature Limited.)

Dysfunction in T cells limits the efficacy of cancer immunotherapy. We profiled the epigenome, transcriptome, and enhancer connectome of exhaustion-prone GD2-targeting HA-28z chimeric antigen receptor (CAR) T cells and control CD19-targeting CAR T cells, which present less exhaustion-inducing tonic signaling, at multiple points during their ex vivo expansion. We found widespread, dynamic changes in chromatin accessibility and three-dimensional (3D) chromosome conformation preceding changes in gene expression, notably at loci proximal to exhaustion-associated genes such as PDCD1 , CTLA4 , and HAVCR2 , and increased DNA motif access for AP-1 family transcription factors, which are known to promote exhaustion. Although T cell exhaustion has been studied in detail in mice, we find that the regulatory networks of T cell exhaustion differ between species and involve distinct loci of accessible chromatin and cis-regulated target genes in human CAR T cell exhaustion. Deletion of exhaustion-specific candidate enhancers of PDCD1 suppress the expression of PD-1 in an in vitro model of T cell dysfunction and in HA-28z CAR T cells, suggesting enhancer editing as a path forward in improving cancer immunotherapy., Competing Interests: Competing interest statement: H.Y.C. is a cofounder of Accent Therapeutics and Boundless Bio and is an advisor to 10× Genomics, Arsenal Bioscience, and Spring Discovery. C.L.M. is a cofounder of Lyell Immunopharma. R.C.L. is employed by and E.W.W. is a consultant for Lyell Immunopharma. A.T.S. is a cofounder of Immunai., (Copyright © 2021 the Author(s). Published by PNAS.)

The advent of single-cell chromatin accessibility profiling has accelerated the ability to map gene regulatory landscapes but has outpaced the development of scalable software to rapidly extract biological meaning from these data. Here we present a software suite for single-cell analysis of regulatory chromatin in R (ArchR; https://www.archrproject.com/ ) that enables fast and comprehensive analysis of single-cell chromatin accessibility data. ArchR provides an intuitive, user-focused interface for complex single-cell analyses, including doublet removal, single-cell clustering and cell type identification, unified peak set generation, cellular trajectory identification, DNA element-to-gene linkage, transcription factor footprinting, mRNA expression level prediction from chromatin accessibility and multi-omic integration with single-cell RNA sequencing (scRNA-seq). Enabling the analysis of over 1.2 million single cells within 8 h on a standard Unix laptop, ArchR is a comprehensive software suite for end-to-end analysis of single-cell chromatin accessibility that will accelerate the understanding of gene regulation at the resolution of individual cells.

Systemic sclerosis (SSc) is a disease at the intersection of autoimmunity and fibrosis. However, the epigenetic regulation and the contributions of diverse cell types to SSc remain unclear. Here we survey, using ATAC-seq, the active DNA regulatory elements of eight types of primary cells in normal skin from healthy controls, as well as clinically affected and unaffected skin from SSc patients. We find that accessible DNA elements in skin-resident dendritic cells (DCs) exhibit the highest enrichment of SSc-associated single-nucleotide polymorphisms (SNPs) and predict the degrees of skin fibrosis in patients. DCs also have the greatest disease-associated changes in chromatin accessibility and the strongest alteration of cell-cell interactions in SSc lesions. Lastly, data from an independent cohort of patients with SSc confirm a significant increase of DCs in lesioned skin. Thus, the DCs epigenome links inherited susceptibility and clinically apparent fibrosis in SSc skin, and can be an important driver of SSc pathogenesis.

Human dendritic cells (DCs) comprise subsets with distinct phenotypic and functional characteristics, but the transcriptional programs that dictate their identity remain elusive. Here, we analyze global chromatin accessibility profiles across resting and stimulated human DC subsets by means of the assay for transposase-accessible chromatin using sequencing (ATAC-seq). We uncover specific regions of chromatin accessibility for each subset and transcriptional regulators of DC function. By comparing plasmacytoid DC responses to IFN-I-producing and non-IFN-I-producing conditions, we identify genetic programs related to their function. Finally, by intersecting chromatin accessibility with genome-wide association studies, we recognize DC subset-specific enrichment of heritability in autoimmune diseases. Our results unravel the basis of human DC subset heterogeneity and provide a framework for their analysis in disease pathogenesis., Competing Interests: Declaration of Interests H.Y.C. is affiliated with Accent Therapeutics, Boundless Bio, 10x Genomics, Arsenal Biosciences, and Spring Discovery. A.T.S. is a scientific co-founder of Immunai and receives research funding from Arsenal Biosciences., (Copyright © 2020 The Author(s). Published by Elsevier Inc. All rights reserved.)

Forebrain development is characterized by highly synchronized cellular processes, which, if perturbed, can cause disease. To chart the regulatory activity underlying these events, we generated a map of accessible chromatin in human three-dimensional forebrain organoids. To capture corticogenesis, we sampled glial and neuronal lineages from dorsal or ventral forebrain organoids over 20 months in vitro. Active chromatin regions identified in human primary brain tissue were observed in organoids at different developmental stages. We used this resource to map genetic risk for disease and to explore evolutionary conservation. Moreover, we integrated chromatin accessibility with transcriptomics to identify putative enhancer-gene linkages and transcription factors that regulate human corticogenesis. Overall, this platform brings insights into gene-regulatory dynamics at previously inaccessible stages of human forebrain development, including signatures of neuropsychiatric disorders., (Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

We develop PIRCh-seq, a method which enables a comprehensive survey of chromatin-associated RNAs in a histone modification-specific manner. We identify hundreds of chromatin-associated RNAs in several cell types with substantially less contamination by nascent transcripts. Non-coding RNAs are found enriched on chromatin and are classified into functional groups based on the patterns of their association with specific histone modifications. We find single-stranded RNA bases are more chromatin-associated, and we discover hundreds of allele-specific RNA-chromatin interactions. These results provide a unique resource to globally study the functions of chromatin-associated lncRNAs and elucidate the basic mechanisms of chromatin-RNA interactions.

Identifying the causes of human diseases requires deconvolution of abnormal molecular phenotypes spanning DNA accessibility, gene expression and protein abundance 1-3 . We present a single-cell framework that integrates highly multiplexed protein quantification, transcriptome profiling and analysis of chromatin accessibility. Using this approach, we establish a normal epigenetic baseline for healthy blood development, which we then use to deconvolve aberrant molecular features within blood from patients with mixed-phenotype acute leukemia 4,5 . Despite widespread epigenetic heterogeneity within the patient cohort, we observe common malignant signatures across patients as well as patient-specific regulatory features that are shared across phenotypic compartments of individual patients. Integrative analysis of transcriptomic and chromatin-accessibility maps identified 91,601 putative peak-to-gene linkages and transcription factors that regulate leukemia-specific genes, such as RUNX1-linked regulatory elements proximal to the marker gene CD69. These results demonstrate how integrative, multiomic analysis of single cells within the framework of normal development can reveal both distinct and shared molecular mechanisms of disease from patient samples.

Oncogenes are commonly amplified on particles of extrachromosomal DNA (ecDNA) in cancer 1,2 , but our understanding of the structure of ecDNA and its effect on gene regulation is limited. Here, by integrating ultrastructural imaging, long-range optical mapping and computational analysis of whole-genome sequencing, we demonstrate the structure of circular ecDNA. Pan-cancer analyses reveal that oncogenes encoded on ecDNA are among the most highly expressed genes in the transcriptome of the tumours, linking increased copy number with high transcription levels. Quantitative assessment of the chromatin state reveals that although ecDNA is packaged into chromatin with intact domain structure, it lacks higher-order compaction that is typical of chromosomes and displays significantly enhanced chromatin accessibility. Furthermore, ecDNA is shown to have a significantly greater number of ultra-long-range interactions with active chromatin, which provides insight into how the structure of circular ecDNA affects oncogene function, and connects ecDNA biology with modern cancer genomics and epigenetics.

Understanding complex tissues requires single-cell deconstruction of gene regulation with precision and scale. Here, we assess the performance of a massively parallel droplet-based method for mapping transposase-accessible chromatin in single cells using sequencing (scATAC-seq). We apply scATAC-seq to obtain chromatin profiles of more than 200,000 single cells in human blood and basal cell carcinoma. In blood, application of scATAC-seq enables marker-free identification of cell type-specific cis- and trans-regulatory elements, mapping of disease-associated enhancer activity and reconstruction of trajectories of cellular differentiation. In basal cell carcinoma, application of scATAC-seq reveals regulatory networks in malignant, stromal and immune cells in the tumor microenvironment. Analysis of scATAC-seq profiles from serial tumor biopsies before and after programmed cell death protein 1 blockade identifies chromatin regulators of therapy-responsive T cell subsets and reveals a shared regulatory program that governs intratumoral CD8 + T cell exhaustion and CD4 + T follicular helper cell development. We anticipate that scATAC-seq will enable the unbiased discovery of gene regulatory factors across diverse biological systems.

Modular domains of long non-coding RNAs can serve as scaffolds to bring distant regions of the linear genome into spatial proximity. Here, we present HiChIRP, a method leveraging bio-orthogonal chemistry and optimized chromosome conformation capture conditions, which enables interrogation of chromatin architecture focused around a specific RNA of interest down to approximately ten copies per cell. HiChIRP of three nuclear RNAs reveals insights into promoter interactions (7SK), telomere biology (telomerase RNA component) and inflammatory gene regulation (lincRNA-EPS).

Tissue development results from lineage-specific transcription factors (TFs) programming a dynamic chromatin landscape through progressive cell fate transitions. Here, we define epigenomic landscape during epidermal differentiation of human pluripotent stem cells (PSCs) and create inference networks that integrate gene expression, chromatin accessibility, and TF binding to define regulatory mechanisms during keratinocyte specification. We found two critical chromatin networks during surface ectoderm initiation and keratinocyte maturation, which are driven by TFAP2C and p63, respectively. Consistently, TFAP2C, but not p63, is sufficient to initiate surface ectoderm differentiation, and TFAP2C-initiated progenitor cells are capable of maturing into functional keratinocytes. Mechanistically, TFAP2C primes the surface ectoderm chromatin landscape and induces p63 expression and binding sites, thus allowing maturation factor p63 to positively autoregulate its own expression and close a subset of the TFAP2C-initiated surface ectoderm program. Our work provides a general framework to infer TF networks controlling chromatin transitions that will facilitate future regenerative medicine advances., (Copyright © 2018. Published by Elsevier Inc.)

We present the genome-wide chromatin accessibility profiles of 410 tumor samples spanning 23 cancer types from The Cancer Genome Atlas (TCGA). We identify 562,709 transposase-accessible DNA elements that substantially extend the compendium of known cis-regulatory elements. Integration of ATAC-seq (the assay for transposase-accessible chromatin using sequencing) with TCGA multi-omic data identifies a large number of putative distal enhancers that distinguish molecular subtypes of cancers, uncovers specific driving transcription factors via protein-DNA footprints, and nominates long-range gene-regulatory interactions in cancer. These data reveal genetic risk loci of cancer predisposition as active DNA regulatory elements in cancer, identify gene-regulatory interactions underlying cancer immune evasion, and pinpoint noncoding mutations that drive enhancer activation and may affect patient survival. These results suggest a systematic approach to understanding the noncoding genome in cancer to advance diagnosis and therapy., (Copyright © 2018 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Understanding the genomic logic that underlies cellular diversity and developmental potential in the human pancreas will accelerate the growth of cell replacement therapies and reveal genetic risk mechanisms in diabetes. Here, we identified and characterized thousands of chromatin regions governing cell-specific gene regulation in human pancreatic endocrine and exocrine lineages, including islet β cells, α cells, duct, and acinar cells. Our findings have captured cellular ontogenies at the chromatin level, identified lineage-specific regulators potentially acting on these sites, and uncovered hallmarks of regulatory plasticity between cell types that suggest mechanisms to regenerate β cells from pancreatic endocrine or exocrine cells. Our work shows that disease risk variants related to pancreas are significantly enriched in these regulatory regions and reveals previously unrecognized links between endocrine and exocrine pancreas in diabetes risk., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Over 90% of genetic variants associated with complex human traits map to non-coding regions, but little is understood about how they modulate gene regulation in health and disease. One possible mechanism is that genetic variants affect the activity of one or more cis-regulatory elements leading to gene expression variation in specific cell types. To identify such cases, we analyzed ATAC-seq and RNA-seq profiles from stimulated primary CD4 + T cells in up to 105 healthy donors. We found that regions of accessible chromatin (ATAC-peaks) are co-accessible at kilobase and megabase resolution, consistent with the three-dimensional chromatin organization measured by in situ Hi-C in T cells. Fifteen percent of genetic variants located within ATAC-peaks affected the accessibility of the corresponding peak (local-ATAC-QTLs). Local-ATAC-QTLs have the largest effects on co-accessible peaks, are associated with gene expression and are enriched for autoimmune disease variants. Our results provide insights into how natural genetic variants modulate cis-regulatory elements, in isolation or in concert, to influence gene expression.

Human hematopoiesis involves cellular differentiation of multipotent cells into progressively more lineage-restricted states. While the chromatin accessibility landscape of this process has been explored in defined populations, single-cell regulatory variation has been hidden by ensemble averaging. We collected single-cell chromatin accessibility profiles across 10 populations of immunophenotypically defined human hematopoietic cell types and constructed a chromatin accessibility landscape of human hematopoiesis to characterize differentiation trajectories. We find variation consistent with lineage bias toward different developmental branches in multipotent cell types. We observe heterogeneity within common myeloid progenitors (CMPs) and granulocyte-macrophage progenitors (GMPs) and develop a strategy to partition GMPs along their differentiation trajectory. Furthermore, we integrated single-cell RNA sequencing (scRNA-seq) data to associate transcription factors to chromatin accessibility changes and regulatory elements to target genes through correlations of expression and regulatory element accessibility. Overall, this work provides a framework for integrative exploration of complex regulatory dynamics in a primary human tissue at single-cell resolution., (Copyright © 2018 Elsevier Inc. All rights reserved.)

T cells create vast amounts of diversity in the genes that encode their T cell receptors (TCRs), which enables individual clones to recognize specific peptide-major histocompatibility complex (MHC) ligands. Here we combined sequencing of the TCR-encoding genes with assay for transposase-accessible chromatin with sequencing (ATAC-seq) analysis at the single-cell level to provide information on the TCR specificity and epigenomic state of individual T cells. By using this approach, termed transcript-indexed ATAC-seq (T-ATAC-seq), we identified epigenomic signatures in immortalized leukemic T cells, primary human T cells from healthy volunteers and primary leukemic T cells from patient samples. In peripheral blood CD4 + T cells from healthy individuals, we identified cis and trans regulators of naive and memory T cell states and found substantial heterogeneity in surface-marker-defined T cell populations. In patients with a leukemic form of cutaneous T cell lymphoma, T-ATAC-seq enabled identification of leukemic and nonleukemic regulatory pathways in T cells from the same individual by allowing separation of the signals that arose from the malignant clone from the background T cell noise. Thus, T-ATAC-seq is a new tool that enables analysis of epigenomic landscapes in clonal T cells and should be valuable for studies of T cell malignancy, immunity and immunotherapy.

The advent of high-throughput epigenome mapping technologies has ushered in a new era of multiomics where powerful tools can now delineate and record different layers of genomic output. Integrating various components of the epigenome from these multiomics measurements allows the interrogation of cellular heterogeneity in addition to the discovery of molecular connectivity maps between the genome and its functional output. Mapping of chromatin accessibility dynamics and higher-order chromatin structure has enabled new levels of understanding of cell fate decisions, identity, and function in normal development, physiology, and disease. We provide a perspective on the progress of the epigenomics field and applications and anticipate an even greater revolution in our understanding of the human epigenome for years to come., (© 2018 American Heart Association, Inc.)

X-chromosome inactivation (XCI) is a critical epigenetic mechanism for balancing gene dosage between XY males and XX females in eutherian mammals. A long non-coding RNA (lncRNA), XIST, and its associated proteins orchestrate this multi-step process, resulting in the inheritable silencing of one of the two X-chromosomes in females. The XIST RNA is large and complex, exemplifying the unique challenges associated with the structural and functional analysis of lncRNAs. Recent technological advances in the analysis of macromolecular structure and interactions have enabled us to systematically dissect the XIST ribonucleoprotein complex, which is larger than the ribosome, and its place of action, the inactive X-chromosome. These studies shed light on key mechanisms of XCI, such as XIST coating of the X-chromosome, recruitment of DNA, RNA and histone modification enzymes, and compaction and compartmentalization of the inactive X. Here, we summarize recent studies on XCI, highlight the critical contributions of new technologies and propose a unifying model for XIST function in XCI where modular domains serve as the structural and functional units in both lncRNA-protein complexes and DNA-protein complexes in chromatin.This article is part of the themed issue 'X-chromosome inactivation: a tribute to Mary Lyon'., (© 2017 The Author(s).)

How transcription factors (TFs) reprogram one cell lineage to another remains unclear. Here, we define chromatin accessibility changes induced by the proneural TF Ascl1 throughout conversion of fibroblasts into induced neuronal (iN) cells. Thousands of genomic loci are affected as early as 12 hr after Ascl1 induction. Surprisingly, over 80% of the accessibility changes occur between days 2 and 5 of the 3-week reprogramming process. This chromatin switch coincides with robust activation of endogenous neuronal TFs and nucleosome phasing of neuronal promoters and enhancers. Subsequent morphological and functional maturation of iN cells is accomplished with relatively little chromatin reconfiguration. By integrating chromatin accessibility and transcriptome changes, we built a network model of dynamic TF regulation during iN cell reprogramming and identified Zfp238, Sox8, and Dlx3 as key TFs downstream of Ascl1. These results reveal a singular, coordinated epigenomic switch during direct reprogramming, in contrast to stepwise cell fate transitions in development., (Copyright © 2017 The Author(s). Published by Elsevier Inc. All rights reserved.)

Rationale: Recent advances have improved our ability to generate cardiomyocytes from human induced pluripotent stem cells (hiPSCs) and human embryonic stem cells (hESCs). However, our understanding of the transcriptional regulatory networks underlying early stages (ie, from mesoderm to cardiac mesoderm) of cardiomyocyte differentiation remains limited., Objective: To characterize transcriptome and chromatin accessibility during early cardiomyocyte differentiation from hiPSCs and hESCs., Methods and Results: We profiled the temporal changes in transcriptome and chromatin accessibility at genome-wide levels during cardiomyocyte differentiation derived from 2 hiPSC lines and 2 hESC lines at 4 stages: pluripotent stem cells, mesoderm, cardiac mesoderm, and differentiated cardiomyocytes. Overall, RNA sequencing analysis revealed that transcriptomes during early cardiomyocyte differentiation were highly concordant between hiPSCs and hESCs, and clustering of 4 cell lines within each time point demonstrated that changes in genome-wide chromatin accessibility were similar across hiPSC and hESC cell lines. Weighted gene co-expression network analysis (WGCNA) identified several modules that were strongly correlated with different stages of cardiomyocyte differentiation. Several novel genes were identified with high weighted connectivity within modules and exhibited coexpression patterns with other genes, including noncoding RNA LINC01124 and uncharacterized RNA AK127400 in the module related to the mesoderm stage; E-box-binding homeobox 1 ( ZEB1 ) in the module correlated with postcardiac mesoderm. We further demonstrated that ZEB1 is required for early cardiomyocyte differentiation. In addition, based on integrative analysis of both WGCNA and transcription factor motif enrichment analysis, we determined numerous transcription factors likely to play important roles at different stages during cardiomyocyte differentiation, such as T and eomesodermin ( EOMES ; mesoderm), lymphoid enhancer-binding factor 1 ( LEF1 ) and mesoderm posterior BHLH transcription factor 1 ( MESP1 ; from mesoderm to cardiac mesoderm), meis homeobox 1 ( MEIS1 ) and GATA-binding protein 4 ( GATA4 ) (postcardiac mesoderm), JUN and FOS families, and MEIS2 (cardiomyocyte)., Conclusions: Both hiPSCs and hESCs share similar transcriptional regulatory mechanisms underlying early cardiac differentiation, and our results have revealed transcriptional regulatory networks and new factors (eg, ZEB1) controlling early stages of cardiomyocyte differentiation., (© 2017 American Heart Association, Inc.)

Here, we define the landscape and dynamics of active regulatory DNA in cutaneous T cell lymphoma (CTCL) by ATAC-seq. Analysis of 111 human CTCL and control samples revealed extensive chromatin signatures that distinguished leukemic, host, and normal CD4 + T cells. We identify three dominant patterns of transcription factor (TF) activation that drive leukemia regulomes, as well as TF deactivations that alter host T cells in CTCL patients. Clinical response to histone deacetylase inhibitors (HDACi) is strongly associated with a concurrent gain in chromatin accessibility. HDACi causes distinct chromatin responses in leukemic and host CD4 + T cells, reprogramming host T cells toward normalcy. These results provide a foundational framework to study personal regulomes in human cancer and epigenetic therapy., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Although only 2% of the genome encodes protein, RNA is transcribed from the majority of the genetic sequence, suggesting a massive degree of cellular functionality is programmed in the noncoding genome. The mammalian genome contains tens of thousands of long noncoding RNAs (lncRNAs), many of which occur at disease-associated loci or are specifically expressed in cancer. Although the vast majority of lncRNAs have no known function, recurring molecular mechanisms for lncRNAs are now being observed in chromatin regulation and cancer pathways and emerging technologies are now providing tools to interrogate lncRNA molecular interactions and determine function of these abundant cellular macromolecules., (Copyright © 2017 Cold Spring Harbor Laboratory Press; all rights reserved.)

Spatial organization of the genome plays a central role in gene expression, DNA replication, and repair. But current epigenomic approaches largely map DNA regulatory elements outside of the native context of the nucleus. Here we report assay of transposase-accessible chromatin with visualization (ATAC-see), a transposase-mediated imaging technology that employs direct imaging of the accessible genome in situ, cell sorting, and deep sequencing to reveal the identity of the imaged elements. ATAC-see revealed the cell-type-specific spatial organization of the accessible genome and the coordinated process of neutrophil chromatin extrusion, termed NETosis. Integration of ATAC-see with flow cytometry enables automated quantitation and prospective cell isolation as a function of chromatin accessibility, and it reveals a cell-cycle dependence of chromatin accessibility that is especially dynamic in G1 phase. The integration of imaging and epigenomics provides a general and scalable approach for deciphering the spatiotemporal architecture of gene control.

Genome conformation is central to gene control but challenging to interrogate. Here we present HiChIP, a protein-centric chromatin conformation method. HiChIP improves the yield of conformation-informative reads by over 10-fold and lowers the input requirement over 100-fold relative to that of ChIA-PET. HiChIP of cohesin reveals multiscale genome architecture with greater signal-to-background ratios than those of in situ Hi-C.

We define the chromatin accessibility and transcriptional landscapes in 13 human primary blood cell types that span the hematopoietic hierarchy. Exploiting the finding that the enhancer landscape better reflects cell identity than mRNA levels, we enable 'enhancer cytometry' for enumeration of pure cell types from complex populations. We identify regulators governing hematopoietic differentiation and further show the lineage ontogeny of genetic elements linked to diverse human diseases. In acute myeloid leukemia (AML), chromatin accessibility uncovers unique regulatory evolution in cancer cells with a progressively increasing mutation burden. Single AML cells exhibit distinctive mixed regulome profiles corresponding to disparate developmental stages. A method to account for this regulatory heterogeneity identified cancer-specific deviations and implicated HOX factors as key regulators of preleukemic hematopoietic stem cell characteristics. Thus, regulome dynamics can provide diverse insights into hematopoietic development and disease., Competing Interests: The authors declare no competing financial interests.

ChIRP is a novel and easy-to-use technique for studying long noncoding RNA (lncRNA)-chromatin interactions. RNA and chromatin are cross-linked in vivo using formaldehyde or glutaraldehyde, and purified using biotinylated antisense oligonucleotides that hybridize to the target RNA. Co-precipitated DNA is then purified and analyzed by quantitative PCR (qPCR) or high-throughput sequencing.

Background: Open chromatin regions are correlated with active regulatory elements in development and are dysregulated in diseases. The BAF (SWI/SNF) complex is essential for development, and has been demonstrated to remodel reconstituted chromatin in vitro and to control the accessibility of a few individual regions in vivo. However, it remains unclear where and how BAF controls the open chromatin landscape to regulate developmental processes, such as human epidermal differentiation., Results: Using a novel "on-plate" ATAC-sequencing approach for profiling open chromatin landscapes with a low number of adherent cells, we demonstrate that the BAF complex is essential for maintaining 11.6 % of open chromatin regions in epidermal differentiation. These BAF-dependent open chromatin regions are highly cell-type-specific and are strongly enriched for binding sites for p63, a master epidermal transcription factor. The DNA sequences of p63 binding sites intrinsically favor nucleosome formation and are inaccessible in other cell types without p63 to prevent ectopic activation. In epidermal cells, BAF and p63 mutually recruit each other to maintain 14,853 open chromatin regions. We further demonstrate that BAF and p63 cooperatively position nucleosomes away from p63 binding sites and recruit transcriptional machinery to control tissue differentiation., Conclusions: BAF displays high specificity in controlling the open chromatin landscape during epidermal differentiation by cooperating with the master transcription factor p63 to maintain lineage-specific open chromatin regions.

Cell-to-cell variation is a universal feature of life that affects a wide range of biological phenomena, from developmental plasticity to tumour heterogeneity. Although recent advances have improved our ability to document cellular phenotypic variation, the fundamental mechanisms that generate variability from identical DNA sequences remain elusive. Here we reveal the landscape and principles of mammalian DNA regulatory variation by developing a robust method for mapping the accessible genome of individual cells by assay for transposase-accessible chromatin using sequencing (ATAC-seq) integrated into a programmable microfluidics platform. Single-cell ATAC-seq (scATAC-seq) maps from hundreds of single cells in aggregate closely resemble accessibility profiles from tens of millions of cells and provide insights into cell-to-cell variation. Accessibility variance is systematically associated with specific trans-factors and cis-elements, and we discover combinations of trans-factors associated with either induction or suppression of cell-to-cell variability. We further identify sets of trans-factors associated with cell-type-specific accessibility variance across eight cell types. Targeted perturbations of cell cycle or transcription factor signalling evoke stimulus-specific changes in this observed variability. The pattern of accessibility variation in cis across the genome recapitulates chromosome compartments de novo, linking single-cell accessibility variation to three-dimensional genome organization. Single-cell analysis of DNA accessibility provides new insight into cellular variation of the 'regulome'.

This unit describes Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq), a method for mapping chromatin accessibility genome-wide. This method probes DNA accessibility with hyperactive Tn5 transposase, which inserts sequencing adapters into accessible regions of chromatin. Sequencing reads can then be used to infer regions of increased accessibility, as well as to map regions of transcription-factor binding and nucleosome position. The method is a fast and sensitive alternative to DNase-seq for assaying chromatin accessibility genome-wide, or to MNase-seq for assaying nucleosome positions in accessible regions of the genome., (Copyright © 2015 John Wiley & Sons, Inc.)

Here we describe domain-specific chromatin isolation by RNA purification (dChIRP), a technique for dissecting the functional domains of a target RNA in situ. For an RNA of interest, dChIRP can identify domain-level intramolecular and intermolecular RNA-RNA, RNA-protein, and RNA-DNA interactions and maps the RNA's genomic binding sites with higher precision than domain-agnostic methods. We illustrate how this technique has been applied to the roX1 lncRNA to resolve its domain-level architecture, discover its protein- and chromatin-interacting domains, and map its occupancy on the X chromosome.

Little is known about the functional domain architecture of long noncoding RNAs (lncRNAs) because of a relative paucity of suitable methods to analyze RNA function at a domain level. Here we describe domain-specific chromatin isolation by RNA purification (dChIRP), a scalable technique to dissect pairwise RNA-RNA, RNA-protein and RNA-chromatin interactions at the level of individual RNA domains in living cells. dChIRP of roX1, a lncRNA essential for Drosophila melanogaster X-chromosome dosage compensation, reveals a 'three-fingered hand' ribonucleoprotein topology. Each RNA finger binds chromatin and the male-specific lethal (MSL) protein complex and can individually rescue male lethality in roX-null flies, thus defining a minimal RNA domain for chromosome-wide dosage compensation. dChIRP improves the RNA genomic localization signal by >20-fold relative to previous techniques, and these binding sites are correlated with chromosome conformation data, indicating that most roX-bound loci cluster in a nuclear territory. These results suggest dChIRP can reveal lncRNA architecture and function with high precision and sensitivity.

The cyclic AMP response element-binding protein (CREB) initiates transcriptional responses to a wide variety of stimuli. CREB activation involves its phosphorylation on Ser-133, which promotes interaction between the CREB kinase-inducible domain (KID) and the KID-interacting domain of the transcriptional coactivator, CREB-binding protein (CBP). The KID also contains a highly conserved phosphorylation cluster, termed the ATM/CK cluster, which is processively phosphorylated in response to DNA damage by the coordinated actions of ataxia-telangiectasia-mutated (ATM) and casein kinases (CKs) 1 and 2. The ATM/CK cluster phosphorylation attenuates CBP binding and CREB transcriptional activity. Paradoxically, it was recently reported that DNA damage activates CREB through homeodomain-interacting protein kinase 2-dependent phosphorylation of Ser-271 near the CREB bZIP DNA binding domain. In this study we sought to further clarify DNA damage-dependent CREB phosphorylation as well as to explore the possibility that the ATM/CK cluster and Ser-271 synergistically or antagonistically modulate CREB activity. We show that, rather than being induced by DNA damage, Ser-270 and Ser-271 of CREB cophosphorylated in a CDK1-dependent manner during G2/M phase. Functionally, we show that phosphorylation of CREB on Ser-270/Ser-271 during mitosis correlated with reduced CREB chromatin occupancy. Furthermore, CDK1-dependent phosphorylation of CREB in vitro inhibited its DNA binding activity. The combined results suggest that CDK1-dependent phosphorylation of CREB on Ser-270/Ser-271 facilitates its dissociation from chromatin during mitosis by reducing its intrinsic DNA binding potential.

Eukaryotic genomes are packaged into chromatin, where diverse histone modifications can demarcate chromatin domains that facilitate or block gene expression. While silent chromatin has been associated with long noncoding RNAs (lncRNAs) for some time, new studies suggest that noncoding RNAs also modulate the active chromatin state. Divergent, antisense, and enhancer-like intergenic noncoding RNAs can either activate or repress gene expression by altering histone H3 lysine 4 methylation. An emerging class of enhancer-like lncRNAs may link chromosome structure to chromatin state and establish active chromatin domains. The confluence of several new technologies promises to rapidly expand this fascinating topic of investigation., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

Long noncoding RNAs are key regulators of chromatin states for important biological processes such as dosage compensation, imprinting, and developmental gene expression (1,2,3,4,5,6,7). The recent discovery of thousands of lncRNAs in association with specific chromatin modification complexes, such as Polycomb Repressive Complex 2 (PRC2) that mediates histone H3 lysine 27 trimethylation (H3K27me3), suggests broad roles for numerous lncRNAs in managing chromatin states in a gene-specific fashion (8,9). While some lncRNAs are thought to work in cis on neighboring genes, other lncRNAs work in trans to regulate distantly located genes. For instance, Drosophila lncRNAs roX1 and roX2 bind numerous regions on the X chromosome of male cells, and are critical for dosage compensation (10,11). However, the exact locations of their binding sites are not known at high resolution. Similarly, human lncRNA HOTAIR can affect PRC2 occupancy on hundreds of genes genome-wide( 3,12,13), but how specificity is achieved is unclear. LncRNAs can also serve as modular scaffolds to recruit the assembly of multiple protein complexes. The classic trans-acting RNA scaffold is the TERC RNA that serves as the template and scaffold for the telomerase complex (14); HOTAIR can also serve as a scaffold for PRC2 and a H3K4 demethylase complex (13). Prior studies mapping RNA occupancy at chromatin have revealed substantial insights (15,16), but only at a single gene locus at a time. The occupancy sites of most lncRNAs are not known, and the roles of lncRNAs in chromatin regulation have been mostly inferred from the indirect effects of lncRNA perturbation. Just as chromatin immunoprecipitation followed by microarray or deep sequencing (ChIP-chip or ChIP-seq, respectively) has greatly improved our understanding of protein-DNA interactions on a genomic scale, here we illustrate a recently published strategy to map long RNA occupancy genome-wide at high resolution (17). This method, Chromatin Isolation by RNA Purification (ChIRP) (Figure 1), is based on affinity capture of target lncRNA:chromatin complex by tiling antisense-oligos, which then generates a map of genomic binding sites at a resolution of several hundred bases with high sensitivity and low background. ChIRP is applicable to many lncRNAs because the design of affinity-probes is straightforward given the RNA sequence and requires no knowledge of the RNA's structure or functional domains.

Long noncoding RNAs (lncRNAs) are key regulators of chromatin state, yet the nature and sites of RNA-chromatin interaction are mostly unknown. Here we introduce Chromatin Isolation by RNA Purification (ChIRP), where tiling oligonucleotides retrieve specific lncRNAs with bound protein and DNA sequences, which are enumerated by deep sequencing. ChIRP-seq of three lncRNAs reveal that RNA occupancy sites in the genome are focal, sequence-specific, and numerous. Drosophila roX2 RNA occupies male X-linked gene bodies with increasing tendency toward the 3' end, peaking at CES sites. Human telomerase RNA TERC occupies telomeres and Wnt pathway genes. HOTAIR lncRNA preferentially occupies a GA-rich DNA motif to nucleate broad domains of Polycomb occupancy and histone H3 lysine 27 trimethylation. HOTAIR occupancy occurs independently of EZH2, suggesting the order of RNA guidance of Polycomb occupancy. ChIRP-seq is generally applicable to illuminate the intersection of RNA and chromatin with newfound precision genome wide., (Copyright © 2011 Elsevier Inc. All rights reserved.)

The sirtuin Sirt6 is a NAD-dependent histone deacetylase that is implicated in gene regulation and lifespan control. Sirt6 can interact with the stress-responsive transcription factor NF-κB and regulate some NF-κB target genes, but the full scope of Sirt6 target genes as well as dynamics of Sirt6 occupancy on chromatin are not known. Here we map Sirt6 occupancy on mouse promoters genome-wide and show that Sirt6 occupancy is highly dynamic in response to TNF-α. More than half of Sirt6 target genes are only revealed upon stress-signaling. The majority of genes bound by NF-κB subunit RelA recruit Sirt6, and dynamic Sirt6 relocalization is largely driven in a RelA-dependent manner. Integrative analysis with global gene expression patterns in wild-type, Sirt6-/-, and double Sirt6-/- RelA-/- cells reveals the epistatic relationships between Sirt6 and RelA in shaping diverse temporal patterns of gene expression. Genes under the direct joint control of Sirt6 and RelA include several with prominent roles in cell senescence and organismal aging. These data suggest dynamic chromatin relocalization of Sirt6 as a key output of NF-κB signaling in stress response and aging., Competing Interests: The authors have declared that no competing interests exist.

Cellular pathways must be synergized, controlled and organized to manage homeostasis. To achieve high selectivity within the crowded cellular milieu the cell utilizes scaffolding complexes whose role is to bring molecules in proximity thereby controlling and enhancing intermolecular interactions and signaling events. To date, scaffolds have been shown to be composed of proteinaceous units; however, recent evidence has supported the idea that non-coding RNAs may also play a similar role. In this point of view article we discuss recent data on ncRNA scaffolds, with particular focus on ncRNA HOTAIR. Using our current knowledge of signaling networks we discuss the role that RNA may play in writing and regulating histone modifications and the information needed for correct gene expression. Further, we speculate on additional, yet undiscovered roles that ncRNAs may be playing as molecular scaffolds.

The genome is extensively transcribed into long intergenic noncoding RNAs (lincRNAs), many of which are implicated in gene silencing. Potential roles of lincRNAs in gene activation are much less understood. Development and homeostasis require coordinate regulation of neighbouring genes through a process termed locus control. Some locus control elements and enhancers transcribe lincRNAs, hinting at possible roles in long-range control. In vertebrates, 39 Hox genes, encoding homeodomain transcription factors critical for positional identity, are clustered in four chromosomal loci; the Hox genes are expressed in nested anterior-posterior and proximal-distal patterns colinear with their genomic position from 3' to 5'of the cluster. Here we identify HOTTIP, a lincRNA transcribed from the 5' tip of the HOXA locus that coordinates the activation of several 5' HOXA genes in vivo. Chromosomal looping brings HOTTIP into close proximity to its target genes. HOTTIP RNA binds the adaptor protein WDR5 directly and targets WDR5/MLL complexes across HOXA, driving histone H3 lysine 4 trimethylation and gene transcription. Induced proximity is necessary and sufficient for HOTTIP RNA activation of its target genes. Thus, by serving as key intermediates that transmit information from higher order chromosomal looping into chromatin modifications, lincRNAs may organize chromatin domains to coordinate long-range gene activation., (©2011 Macmillan Publishers Limited. All rights reserved)

Long intergenic noncoding RNAs (lincRNAs) regulate chromatin states and epigenetic inheritance. Here, we show that the lincRNA HOTAIR serves as a scaffold for at least two distinct histone modification complexes. A 5' domain of HOTAIR binds polycomb repressive complex 2 (PRC2), whereas a 3' domain of HOTAIR binds the LSD1/CoREST/REST complex. The ability to tether two distinct complexes enables RNA-mediated assembly of PRC2 and LSD1 and coordinates targeting of PRC2 and LSD1 to chromatin for coupled histone H3 lysine 27 methylation and lysine 4 demethylation. Our results suggest that lincRNAs may serve as scaffolds by providing binding surfaces to assemble select histone modification enzymes, thereby specifying the pattern of histone modifications on target genes.

Large intervening non-coding RNAs (lincRNAs) are pervasively transcribed in the genome yet their potential involvement in human disease is not well understood. Recent studies of dosage compensation, imprinting, and homeotic gene expression suggest that individual lincRNAs can function as the interface between DNA and specific chromatin remodelling activities. Here we show that lincRNAs in the HOX loci become systematically dysregulated during breast cancer progression. The lincRNA termed HOTAIR is increased in expression in primary breast tumours and metastases, and HOTAIR expression level in primary tumours is a powerful predictor of eventual metastasis and death. Enforced expression of HOTAIR in epithelial cancer cells induced genome-wide re-targeting of Polycomb repressive complex 2 (PRC2) to an occupancy pattern more resembling embryonic fibroblasts, leading to altered histone H3 lysine 27 methylation, gene expression, and increased cancer invasiveness and metastasis in a manner dependent on PRC2. Conversely, loss of HOTAIR can inhibit cancer invasiveness, particularly in cells that possess excessive PRC2 activity. These findings indicate that lincRNAs have active roles in modulating the cancer epigenome and may be important targets for cancer diagnosis and therapy.