Lineage-Specific Genome Architecture Links Enhancers and Non-coding Disease Variants to Target Gene Promoters.

Long-range interactions between regulatory elements and gene promoters play key roles in transcriptional regulation. The vast majority of interactions are uncharted, constituting a major missing link in understanding genome control. Here, we use promoter capture Hi-C to identify interacting regions of 31,253 promoters in 17 human primary hematopoietic cell types. We show that promoter interactions are highly cell type specific and enriched for links between active promoters and epigenetically marked enhancers. Promoter interactomes reflect lineage relationships of the hematopoietic tree, consistent with dynamic remodeling of nuclear architecture during differentiation. Interacting regions are enriched in genetic variants linked with altered expression of genes they contact, highlighting their functional role. We exploit this rich resource to connect non-coding disease variants to putative target promoters, prioritizing thousands of disease-candidate genes and implicating disease pathways. Our results demonstrate the power of primary cell promoter interactomes to reveal insights into genomic regulatory mechanisms underlying common diseases.

High-resolution simulations of chromatin folding at genomic rearrangements in malignant B cells provide mechanistic insights into proto-oncogene deregulation.

Genomic rearrangements are known to result in proto-oncogene deregulation in many cancers, but the link to 3D genome structure remains poorly understood. Here, we used the highly predictive heteromorphic polymer (HiP-HoP) model to predict chromatin conformations at the proto-oncogene CCND1 in healthy and malignant B cells. After confirming that the model gives good predictions of Hi-C data for the nonmalignant human B cell-derived cell line GM12878, we generated predictions for two cancer cell lines, U266 and Z-138. These possess genome rearrangements involving CCND1 and the immunoglobulin heavy locus (IGH), which we mapped using targeted genome sequencing. Our simulations showed that a rearrangement in U266 cells where a single IGH super-enhancer is inserted next to CCND1 leaves the local topologically associated domain (TAD) structure intact. We also observed extensive changes in enhancer-promoter interactions within the TAD, suggesting that it is the downstream chromatin remodeling which gives rise to the oncogene activation, rather than the presence of the inserted super-enhancer DNA sequence per se. Simulations of the IGH-CCND1 reciprocal translocation in Z-138 cells revealed that an oncogenic fusion TAD is created, encompassing CCND1 and the IGH super-enhancers. We predicted how the structure and expression of CCND1 changes in these different cell lines, validating this using qPCR and fluorescence in situ hybridization microscopy. Our work demonstrates the power of polymer simulations to predict differences in chromatin interactions and gene expression for different translocation breakpoints.

The HDAC7-TET2 epigenetic axis is essential during early B lymphocyte development.

Correct B cell identity at each stage of cellular differentiation during B lymphocyte development is critically dependent on a tightly controlled epigenomic landscape. We previously identified HDAC7 as an essential regulator of early B cell development and its absence leads to a drastic block at the pro-B to pre-B cell transition. More recently, we demonstrated that HDAC7 loss in pro-B-ALL in infants associates with a worse prognosis. Here we delineate the molecular mechanisms by which HDAC7 modulates early B cell development. We find that HDAC7 deficiency drives global chromatin de-condensation, histone marks deposition and deregulates other epigenetic regulators and mobile elements. Specifically, the absence of HDAC7 induces TET2 expression, which promotes DNA 5-hydroxymethylation and chromatin de-condensation. HDAC7 deficiency also results in the aberrant expression of microRNAs and LINE-1 transposable elements. These findings shed light on the mechanisms by which HDAC7 loss or misregulation may lead to B cell-based hematological malignancies.

MAX mutant small-cell lung cancers exhibit impaired activities of MGA-dependent noncanonical polycomb repressive complex.

The MYC axis is disrupted in cancer, predominantly through activation of the MYC family oncogenes but also through inactivation of the MYC partner MAX or of the MAX partner MGA. MGA and MAX are also members of the polycomb repressive complex, ncPRC1.6. Here, we use genetically modified MAX-deficient small-cell lung cancer (SCLC) cells and carry out genome-wide and proteomics analyses to study the tumor suppressor function of MAX. We find that MAX mutant SCLCs have ASCL1 or NEUROD1 or combined ASCL1/NEUROD1 characteristics and lack MYC transcriptional activity. MAX restitution triggers prodifferentiation expression profiles that shift when MAX and oncogenic MYC are coexpressed. Although ncPRC1.6 can be formed, the lack of MAX restricts global MGA occupancy, selectively driving its recruitment toward E2F6-binding motifs. Conversely, MAX restitution enhances MGA occupancy to repress genes involved in different functions, including stem cell and DNA repair/replication. Collectively, these findings reveal that MAX mutant SCLCs have either ASCL1 or NEUROD1 or combined characteristics and are MYC independent and exhibit deficient ncPRC1.6-mediated gene repression.

From Loops to Looks: Transcription Factors and Chromatin Organization Shaping Terminal B Cell Differentiation.

B lymphopoiesis is tightly regulated at the level of gene transcription. In recent years, investigators have shed light on the transcription factor networks and the epigenetic machinery involved at all differentiation steps of mammalian B cell development. During terminal differentiation, B cells undergo dramatic changes in gene transcriptional programs to generate germinal center B cells, plasma cells and memory B cells. Recent evidence indicates that mature B cell formation involves an essential contribution from 3D chromatin conformations through its interplay with transcription factors and epigenetic machinery. Here, we provide an up-to-date overview of the coordination between transcription factors, epigenetic changes, and chromatin architecture during terminal B cell differentiation, focusing on recent discoveries and technical advances for studying 3D chromatin structures.

p53 rapidly restructures 3D chromatin organization to trigger a transcriptional response.

Activation of the p53 tumor suppressor triggers a transcriptional program to control cellular response to stress. However, the molecular mechanisms by which p53 controls gene transcription are not completely understood. Here, we uncover the critical role of spatio-temporal genome architecture in this process. We demonstrate that p53 drives direct and indirect changes in genome compartments, topologically associating domains, and DNA loops prior to one hour of its activation, which escort the p53 transcriptional program. Focusing on p53-bound enhancers, we report 340 genes directly regulated by p53 over a median distance of 116 kb, with 74% of these genes not previously identified. Finally, we showcase that p53 controls transcription of distal genes through newly formed and pre-existing enhancer-promoter loops in a cohesin dependent manner. Collectively, our findings demonstrate a previously unappreciated architectural role of p53 as regulator at distinct topological layers and provide a reliable set of new p53 direct target genes that may help designs of cancer therapies.

Changes in the pattern of DNA methylation associate with twin discordance in systemic lupus erythematosus.

Monozygotic (MZ) twins are partially concordant for most complex diseases, including autoimmune disorders. Whereas phenotypic concordance can be used to study heritability, discordance suggests the role of non-genetic factors. In autoimmune diseases, environmentally driven epigenetic changes are thought to contribute to their etiology. Here we report the first high-throughput and candidate sequence analyses of DNA methylation to investigate discordance for autoimmune disease in twins. We used a cohort of MZ twins discordant for three diseases whose clinical signs often overlap: systemic lupus erythematosus (SLE), rheumatoid arthritis, and dermatomyositis. Only MZ twins discordant for SLE featured widespread changes in the DNA methylation status of a significant number of genes. Gene ontology analysis revealed enrichment in categories associated with immune function. Individual analysis confirmed the existence of DNA methylation and expression changes in genes relevant to SLE pathogenesis. These changes occurred in parallel with a global decrease in the 5-methylcytosine content that was concomitantly accompanied with changes in DNA methylation and expression levels of ribosomal RNA genes, although no changes in repetitive sequences were found. Our findings not only identify potentially relevant DNA methylation markers for the clinical characterization of SLE patients but also support the notion that epigenetic changes may be critical in the clinical manifestations of autoimmune disease.

An Integrated Workflow to Study the Promoter-Centric Spatio-Temporal Genome Architecture in Scarce Cell Populations.

Spatiotemporal gene transcription is tightly regulated by distal regulatory elements, such as enhancers and silencers, which rely on physical proximity with their target gene promoters to control transcription. Although these regulatory elements are easy to identify, their target genes are difficult to predict, since most of them are cell-type specific and may be separated by hundreds of kilobases in the linear genome sequence, skipping over other non-target genes. For several years, Promoter Capture Hi-C (PCHi-C) has been the gold standard for the association of distal regulatory elements to their target genes. However, PCHi-C relies on the availability of millions of cells, prohibiting the study of rare cell populations such as those commonly obtained from primary tissues. To overcome this limitation, low input Capture Hi-C (liCHi-C), a cost-effective and customizable method to identify the repertoire of distal regulatory elements controlling each gene of the genome, has been developed. liCHi-C relies on a similar experimental and computational framework as PCHi-C, but by employing minimal tube changes, modifying the reagent concentration and volumes, and swapping or eliminating steps, it accounts for minimal material loss during library construction. Collectively, liCHi-C enables the study of gene regulation and spatiotemporal genome organization in the context of developmental biology and cellular function.

Low input capture Hi-C (liCHi-C) identifies promoter-enhancer interactions at high-resolution.

Long-range interactions between regulatory elements and promoters are key in gene transcriptional control; however, their study requires large amounts of starting material, which is not compatible with clinical scenarios nor the study of rare cell populations. Here we introduce low input capture Hi-C (liCHi-C) as a cost-effective, flexible method to map and robustly compare promoter interactomes at high resolution. As proof of its broad applicability, we implement liCHi-C to study normal and malignant human hematopoietic hierarchy in clinical samples. We demonstrate that the dynamic promoter architecture identifies developmental trajectories and orchestrates transcriptional transitions during cell-state commitment. Moreover, liCHi-C enables the identification of disease-relevant cell types, genes and pathways potentially deregulated by non-coding alterations at distal regulatory elements. Finally, we show that liCHi-C can be harnessed to uncover genome-wide structural variants, resolve their breakpoints and infer their pathogenic effects. Collectively, our optimized liCHi-C method expands the study of 3D chromatin organization to unique, low-abundance cell populations, and offers an opportunity to uncover factors and regulatory networks involved in disease pathogenesis.

Intra- and interchromosomal contact mapping reveals the Igh locus has extensive conformational heterogeneity and interacts with B-lineage genes.

To produce a diverse antibody repertoire, immunoglobulin heavy-chain (Igh) loci undergo large-scale alterations in structure to facilitate juxtaposition and recombination of spatially separated variable (VH), diversity (DH), and joining (JH) genes. These chromosomal alterations are poorly understood. Uncovering their patterns shows how chromosome dynamics underpins antibody diversity. Using tiled Capture Hi-C, we produce a comprehensive map of chromatin interactions throughout the 2.8-Mb Igh locus in progenitor B cells. We find that the Igh locus folds into semi-rigid subdomains and undergoes flexible looping of the VH genes to its 3' end, reconciling two views of locus organization. Deconvolution of single Igh locus conformations using polymer simulations identifies thousands of different structures. This heterogeneity may underpin the diversity of V(D)J recombination events. All three immunoglobulin loci also participate in a highly specific, developmentally regulated network of interchromosomal interactions with genes encoding B cell-lineage factors. This suggests a model of interchromosomal coordination of B cell development.

Neural crest-related NXPH1/α-NRXN signaling opposes neuroblastoma malignancy by inhibiting organotropic metastasis.

Neuroblastoma is a pediatric cancer that can present as low- or high-risk tumors (LR-NBs and HR-NBs), the latter group showing poor prognosis due to metastasis and strong resistance to current therapy. Whether LR-NBs and HR-NBs differ in the way they exploit the transcriptional program underlying their neural crest, sympatho-adrenal origin remains unclear. Here, we identified the transcriptional signature distinguishing LR-NBs from HR-NBs, which consists mainly of genes that belong to the core sympatho-adrenal developmental program and are associated with favorable patient prognosis and with diminished disease progression. Gain- and loss-of-function experiments revealed that the top candidate gene of this signature, Neurexophilin-1 (NXPH1), has a dual impact on NB cell behavior in vivo: whereas NXPH1 and its receptor α-NRXN1 promote NB tumor growth by stimulating cell proliferation, they conversely inhibit organotropic colonization and metastasis. As suggested by RNA-seq analyses, these effects might result from the ability of NXPH1/α-NRXN signalling to restrain the conversion of NB cells from an adrenergic state to a mesenchymal one. Our findings thus uncover a transcriptional module of the sympatho-adrenal program that opposes neuroblastoma malignancy by impeding metastasis, and pinpoint NXPH1/α-NRXN signaling as a promising target to treat HR-NBs.

The onset of PI3K-related vascular malformations occurs during angiogenesis and is prevented by the AKT inhibitor miransertib.

Low-flow vascular malformations are congenital overgrowths composed of abnormal blood vessels potentially causing pain, bleeding and obstruction of different organs. These diseases are caused by oncogenic mutations in the endothelium, which result in overactivation of the PI3K/AKT pathway. Lack of robust in vivo preclinical data has prevented the development and translation into clinical trials of specific molecular therapies for these diseases. Here, we demonstrate that the Pik3caH1047R activating mutation in endothelial cells triggers a transcriptome rewiring that leads to enhanced cell proliferation. We describe a new reproducible preclinical in vivo model of PI3K-driven vascular malformations using the postnatal mouse retina. We show that active angiogenesis is required for the pathogenesis of vascular malformations caused by activating Pik3ca mutations. Using this model, we demonstrate that the AKT inhibitor miransertib both prevents and induces the regression of PI3K-driven vascular malformations. We confirmed the efficacy of miransertib in isolated human endothelial cells with genotypes spanning most of human low-flow vascular malformations.

A new epigenetic challenge: systemic lupus erythematosus.

In recent years, compelling evidence has been gathered that supports a role for epigenetic alterations in the pathogenesis of systemic lupus erythematosus (SLE). Different blood cell populations of SLE patients are characterized by a global loss of DNA methylation. This process is associated with defects in ERK pathway signalling and consequent DNMT 1 downregulation. Hypomethylation of gene promoters has been described, which permits transcriptional activation and therefore functional changes in the cells and also hypomethylation of the ribosomal RNA gene cluster. Among the identified targets undergoing demethylation are genes involved in autoreactivity (ITGAL), osmotic lysis and apoptosis (PRF1, MMP14 and LCN2), antigen presentation (CSF3R), inflammation(MMP 14), B- T-cell interaction (CD70 and CD40LG) and cytokine pathways (CSF3R, IL-4, IL-6 and IFNGR2). DNA methylation inhibitors are also known to induce autoreactivity in vitro and cause a lupus-like disease in vivo. Further, altered patterns of histone modifications have been described in SLE. CD4+ lymphocytes undergo global histone H3 and H4 deacetylation and consequent skewed gene expression. Although multiple lines of evidence highlight the contribution of epigenetic alterations to the pathogenesis of lupus in genetically predisposed individuals, many questions remain to be answered. Attaining a deeper understanding of these matters will create opportunities in the promising area of epigenetic treatments.

The Genome in a Three-Dimensional Context: Deciphering the Contribution of Noncoding Mutations at Enhancers to Blood Cancer.

Associations between blood cancer and genetic predisposition, including both inherited variants and acquired mutations and epimutations, have been well characterized. However, the majority of these variants affect noncoding regions, making their mechanisms difficult to hypothesize and hindering the translation of these insights into patient benefits. Fueled by unprecedented progress in next-generation sequencing and computational integrative analysis, studies have started applying combinations of epigenetic, genome architecture, and functional assays to bridge the gap between noncoding variants and blood cancer. These complementary tools have not only allowed us to understand the potential malignant role of these variants but also to differentiate key variants, cell-types, and conditions from misleading ones. Here, we briefly review recent studies that have provided fundamental insights into our understanding of how noncoding mutations at enhancers predispose and promote blood malignancies in the context of spatial genome architecture.

Epigenomics and transcriptomics of systemic sclerosis CD4+ T cells reveal long-range dysregulation of key inflammatory pathways mediated by disease-associated susceptibility loci.

BACKGROUND: Systemic sclerosis (SSc) is a genetically complex autoimmune disease mediated by the interplay between genetic and epigenetic factors in a multitude of immune cells, with CD4+ T lymphocytes as one of the principle drivers of pathogenesis. METHODS: DNA samples exacted from CD4+ T cells of 48 SSc patients and 16 healthy controls were hybridized on MethylationEPIC BeadChip array. In parallel, gene expression was interrogated by hybridizing total RNA on Clariom™ S array. Downstream bioinformatics analyses were performed to identify correlating differentially methylated CpG positions (DMPs) and differentially expressed genes (DEGs), which were then confirmed utilizing previously published promoter capture Hi-C (PCHi-C) data. RESULTS: We identified 9112 and 3929 DMPs and DEGs, respectively. These DMPs and DEGs are enriched in functional categories related to inflammation and T cell biology. Furthermore, correlation analysis identified 17,500 possible DMP-DEG interaction pairs within a window of 5 Mb, and utilizing PCHi-C data, we observed that 212 CD4+ T cell-specific pairs of DMP-DEG also formed part of three-dimensional promoter-enhancer networks, potentially involving CTCF. Finally, combining PCHi-C data with SSc GWAS data, we identified four important SSc-associated susceptibility loci, TNIP1 (rs3792783), GSDMB (rs9303277), IL12RB1 (rs2305743), and CSK (rs1378942), that could potentially interact with DMP-DEG pairs cg17239269-ANXA6, cg19458020-CCR7, cg10808810-JUND, and cg11062629-ULK3, respectively. CONCLUSION: Our study unveils a potential link between genetic, epigenetic, and transcriptional deregulation in CD4+ T cells of SSc patients, providing a novel integrated view of molecular components driving SSc pathogenesis.

Highly interconnected enhancer communities control lineage-determining genes in human mesenchymal stem cells.

Adipocyte differentiation is driven by waves of transcriptional regulators that reprogram the enhancer landscape and change the wiring of the promoter interactome. Here, we use high-throughput chromosome conformation enhancer capture to interrogate the role of enhancer-to-enhancer interactions during differentiation of human mesenchymal stem cells. We find that enhancers form an elaborate network that is dynamic during differentiation and coupled with changes in enhancer activity. Transcription factors (TFs) at baited enhancers amplify TF binding at target enhancers, a phenomenon we term cross-interaction stabilization of TFs. Moreover, highly interconnected enhancers (HICE) act as integration hubs orchestrating differentiation by the formation of three-dimensional enhancer communities, inside which, HICE, and other enhancers, converge on phenotypically important gene promoters. Collectively, these results indicate that enhancer interactions play a key role in the regulation of enhancer function, and that HICE are important for both signal integration and compartmentalization of the genome.

Human pancreatic islet three-dimensional chromatin architecture provides insights into the genetics of type 2 diabetes.

Genetic studies promise to provide insight into the molecular mechanisms underlying type 2 diabetes (T2D). Variants associated with T2D are often located in tissue-specific enhancer clusters or super-enhancers. So far, such domains have been defined through clustering of enhancers in linear genome maps rather than in three-dimensional (3D) space. Furthermore, their target genes are often unknown. We have created promoter capture Hi-C maps in human pancreatic islets. This linked diabetes-associated enhancers to their target genes, often located hundreds of kilobases away. It also revealed >1,300 groups of islet enhancers, super-enhancers and active promoters that form 3D hubs, some of which show coordinated glucose-dependent activity. We demonstrate that genetic variation in hubs impacts insulin secretion heritability, and show that hub annotations can be used for polygenic scores that predict T2D risk driven by islet regulatory variants. Human islet 3D chromatin architecture, therefore, provides a framework for interpretation of T2D genome-wide association study (GWAS) signals.

Author Correction: Promoter interactome of human embryonic stem cell-derived cardiomyocytes connects GWAS regions to cardiac gene networks.

In the original version of the Article, the gene symbol for tissue factor pathway inhibitor was inadvertently given as 'TFP1' instead of 'TFPI'. This has now been corrected in both the PDF and HTML versions of the Article.

Promoter interactome of human embryonic stem cell-derived cardiomyocytes connects GWAS regions to cardiac gene networks.

Long-range chromosomal interactions bring distal regulatory elements and promoters together to regulate gene expression in biological processes. By performing promoter capture Hi-C (PCHi-C) on human embryonic stem cell-derived cardiomyocytes (hESC-CMs), we show that such promoter interactions are a key mechanism by which enhancers contact their target genes after hESC-CM differentiation from hESCs. We also show that the promoter interactome of hESC-CMs is associated with expression quantitative trait loci (eQTLs) in cardiac left ventricular tissue; captures the dynamic process of genome reorganisation after hESC-CM differentiation; overlaps genome-wide association study (GWAS) regions associated with heart rate; and identifies new candidate genes in such regions. These findings indicate that regulatory elements in hESC-CMs identified by our approach control gene expression involved in ventricular conduction and rhythm of the heart. The study of promoter interactions in other hESC-derived cell types may be of utility in functional investigation of GWAS-associated regions.

The reference epigenome and regulatory chromatin landscape of chronic lymphocytic leukemia.

Chronic lymphocytic leukemia (CLL) is a frequent hematological neoplasm in which underlying epigenetic alterations are only partially understood. Here, we analyze the reference epigenome of seven primary CLLs and the regulatory chromatin landscape of 107 primary cases in the context of normal B cell differentiation. We identify that the CLL chromatin landscape is largely influenced by distinct dynamics during normal B cell maturation. Beyond this, we define extensive catalogues of regulatory elements de novo reprogrammed in CLL as a whole and in its major clinico-biological subtypes classified by IGHV somatic hypermutation levels. We uncover that IGHV-unmutated CLLs harbor more active and open chromatin than IGHV-mutated cases. Furthermore, we show that de novo active regions in CLL are enriched for NFAT, FOX and TCF/LEF transcription factor family binding sites. Although most genetic alterations are not associated with consistent epigenetic profiles, CLLs with MYD88 mutations and trisomy 12 show distinct chromatin configurations. Furthermore, we observe that non-coding mutations in IGHV-mutated CLLs are enriched in H3K27ac-associated regulatory elements outside accessible chromatin. Overall, this study provides an integrative portrait of the CLL epigenome, identifies extensive networks of altered regulatory elements and sheds light on the relationship between the genetic and epigenetic architecture of the disease.