Long-range DNA Interactions Research Articles

EPCO-46. CHARACTERIZING ATYPICAL TERATOID RHABDOID TUMOR SUBTYPE ENHANCER INTERACTION LANDSCAPE

Abstract Atypical Teratoid Rhabdoid Tumor (ATRT) is one of the most aggressive pediatric tumors of the central nervous system (CNS). The majority of ATRTs are characterized by the recurrent bi-allelic loss of the SWI/SNF nucleosome remodeling complex core component gene, SMARCB1, or SMARCA4 in rare cases. Despite its bland genome, three ATRT subgroups “SHH”, “TYR”, and “MYC” were discovered with unique molecular profiles. Also, distinct activating histone mark H3K27ac profiles between molecular subgroups were identified, suggesting subgroup-specific regulatory networks through enhancers. However, comprehensive cis-regulatory network descriptions in ATRT are lacking. Here, we aim to elucidate ATRT subgroup pathogenesis through its dysregulated enhancer-promoter interactions. We have created a regulatory mutational profile of the three ATRT subgroups using whole-genome sequencing method of 60 ATRT tumors and matching genome-wide H3K27ac profiles on 16 tumors. To investigate mutation burden on cis-regulations, we have employed high-throughput Hi-C chromatin immunoprecipitation (HiChIP) technique to directly map enhancer sequences involved in DNA-looping with high specificity and sensitivity. Recent findings from our lab suggest that activating interaction PDGFRB confer to therapeutics venerability in the MYC subgroup. Here, our preliminary analysis shows a hub of interactions near PDGFRB reveals an extensive interaction hub with long-range enhancer interactions extending as far as 1500Kb, suggesting a much more complicated regulatory picture than previously expected. Interestingly, we have identified long-range DNA interactions between PDGFRB and H3K27ac marks upstream of its gene body and overstretching into other gene elements. This result suggests a potential alternative activation mechanism for MYC cell lines. Our next steps are to perform performing additional HiChIP assays on other ATRT cell lines from other subgroups to capture the full range of ATRT’s heterogeneity. Comprehensive mapping of the non-coding genome and enhancer regulatory network would help identify critical pathogenic elements underlying ATRT biology and provide insights into therapeutic targets.

G-quadruplex DNA contributes to RNA polymerase II-mediated 3D chromatin architecture.

High-order chromatin organization plays an important role in biological processes and disease development. Previous studies revealed a widespread occurrence of guanine quadruplex (G4) structures in the human genome, with enrichment in gene regulatory regions, especially in promoters. However, it remains unclear whether G4 structures contribute to RNA polymerase II (RNAPII)-mediated long-range DNA interactions and transcription activity. In this study, we conducted an intuitive overlapping analysis of previously published RNAPII ChIA-PET (chromatin interaction analysis with paired-end tag) and BG4 ChIP-seq (chromatin immunoprecipitation followed by sequencing using a G4 structure-specific antibody) data. We observed a strong positive correlation between RNAPII-linked DNA loops and G4 structures in chromatin. Additionally, our RNAPII HiChIP-seq (in situ Hi-C followed by ChIP-seq) results showed that treatment of HepG2 cells with pyridostatin (PDS), a small-molecule G4-binding ligand, could diminish RNAPII-linked long-range DNA contacts, with more pronounced diminutions being observed for those contacts involving G4 structure loci. RNA sequencing data revealed that PDS treatment modulates the expression of not only genes with G4 structures in their promoters, but also those with promoters being connected with distal G4sthrough RNAPII-linked long-range DNA interactions. Together, our data substantiate the function of DNA G4s in RNAPII-associated DNA looping and transcription regulation.

Technologies, strategies, and cautions when deconvoluting genome-wide association signals: FTO in focus.

Genome-wide association studies have revealed a plethora of genetic variants that correlate with polygenic conditions. However, causal molecular mechanisms have proven challenging to fully define. Without such information, the associations are not physiologically useful or clinically actionable. By reviewing studies of the FTO locus in the genetic etiology of obesity, we wish to highlight advances in the field fueled by the evolution of technical and analytic strategies in assessing the molecular bases for genetic associations. Particular attention is drawn to extrapolating experimental findings from animal models and cell types to humans, as well as technical aspects used to identify long-range DNA interactions and their biological relevance with regard to the associated trait. A unifying model is proposed by which independent obesogenic pathways regulated by multiple FTO variants and genes are integrated at the primary cilium, a cellular antenna where signaling molecules that control energy balance convene.

Long-range DNA interactions: inter-molecular G-quadruplexes and their potential biological relevance.

Guanine-rich DNA sequences are known to fold into secondary structures called G-quadruplexes (G4s), which can form from either individual DNA strands (intra-molecular) or multiple DNA strands (inter-molecular, iG4s). Intra-molecular G4s have been the object of extensive biological investigation due to their enrichment in gene-promoters and telomers. On the other hand, iG4s have never been considered in biological contexts, as the interaction between distal sequences of DNA to form an iG4 in cells was always deemed as highly unlikely. In this feature article, we challenge this dogma by presenting our recent discovery of the first human protein (CSB) displaying astonishing picomolar affinity and binding selectivity for iG4s. These findings suggest potential for iG4 structures to form in cells and highlight the need of further studies to unravel the fundamental biological roles of these inter-molecular DNA structures. Furthermore, we discuss how the potential for formation of iG4s in neuronal cells, triggered by repeat expansions in the C9orf72 gene, can lead to the formation of nucleic-acids based pathological aggregates in neurodegenerative diseases like ALS and FTD. Finally, based on our recent work on short LNA-modified probes, we provide a prespective on how the rational design of G4-selective chemical tools can be leveraged to further elucidate the biological relevance of iG4 structures in the context of ageing-related diseases.

The connectome of neural crest enhancers reveals regulatory features of signaling systems.

Cell fate commitment is controlled by cis-regulatory elements often located in remote regions of the genome. To examine the role of long-range DNA interactions in early development, we generated a high-resolution contact map of active enhancers in avian neural crest cells. This analysis uncovered a diverse repertoire of enhancers that are part of the gene regulatory network underlying specification. We found that neural crest identity is largely regulated by cis-regulatory elements that propagate signaling inputs to network components. These genomic sensors display a combination of optimal and suboptimal TCF/LEF-binding sites, which allow cells to respond to Wnt signaling in a position-dependent manner. We propose that, rather than acting as upstream activators, signaling systems feed into regulatory circuits in a hub-and-spoke architecture. These results shed light on the tridimensional organization of the neural crest genome and define how signaling systems provide progenitors with spatial cues that transform their molecular identity.

Chromosomal Looping-Based Expression Activation System in Yeast.

There is growing evidence that 3D genome organization is a universal and significant mechanism of gene expression regulation. Tools to manipulate long-range DNA interactions can advance the accurate control of chromosomal architecture. However, simple eukaryotic systems available to engineer chromosomal looping and gene expression are very limited. This study has developed a tool designated as chromosomal looping-based expression activation system in yeast (CLEASY). Based on a modified yeast chromosome, it consists of conditionally interacting proteins, distal transcriptional regulatory elements, and a reporter gene. Exogenous chemical or light exposure induces the protein interaction, and results in the proximity of distal regulatory elements bound by these interacting proteins, and ultimately activates the reporter gene. In addition to this controllable induction, this system is compatible with the bivalent Cas9 complexes and their guide RNAs, to ensure target specificity and variability. Therefore, CLEASY can be utilized as a simplified eukaryotic model to engineer DNA looping machinery, and potentially serves as a fast platform to investigate looping mechanism and effective molecules.

Model-based analysis of chromatin interactions from dCas9-Based CAPTURE-3C-seq.

Deciphering long-range chromatin interactions is critical for understanding temporal and tissue-specific gene expression regulated by cis- and trans-acting factors. By combining the chromosome conformation capture (3C) and biotinylated dCas9 system, we previously established a method CAPTURE-3C-seq to unbiasedly identify high-resolution and locus-specific long-range DNA interactions. Here we present the statistical model and a flexible pipeline, C3S, for analysing CAPTURE-3C-seq or similar experimental data from raw sequencing reads to significantly interacting chromatin loci. C3S provides all steps for data processing, quality control and result illustration. It can automatically define the bin size based on the binding peak of the dCas9-targeted regions. Furthermore, it supports the analysis of intra- and inter-chromosomal interactions for different mammalian cell types. We successfully applied C3S across multiple datasets in human K562 cells and mouse embryonic stem cells (mESC) for detecting known and new chromatin interactions at multiple scales. Integrative and topological analysis of the interacted loci at the human β-globin gene cluster provides new insights into mechanisms in developmental gene regulation and network structure in local chromosomal architecture. Furthermore, computational results in mESCs reveal a role for chromatin interacting loops between enhancers and promoters in regulating alternative transcripts of the pluripotency gene OCT4.

PDTM-22. A C19MC-LIN28A-MYCN ONCOGENIC CIRCUIT DRIVEN BY HIJACKED SUPER-ENHANCERS IS A DISTINCT THERAPEUTIC VULNERABILITY IN ETMRS – A LETHAL BRAIN TUMOR

Abstract Embryonal tumors with multi-layered rosettes (ETMR) are aggressive brain cancers with characteristic C19MC oncomiR amplification and enrichment of pluripotency factor LIN28A. Here we investigated C19MC oncogenic mechanisms and discovered a potent C19MC-LIN28A-MYCN circuitry driven by multiple regulatory loops and super-enhancers resulting from long-range MYCN DNA interactions and C19MC gene fusions. C19MC and LIN28A targets respectively converge on critical cell cycle tumor suppressors and neo-embryonic DNMT3A/B isoforms. We identify a MYCN driven, core transcriptional network, conserved in early neural stem cells, that is potently abrogated by treatment with bromodomain inhibitor JQ1, leading to ETMR cell death. Our collective data suggest C19MC oncomiRs drive a malignant primitive cell state in ETMRs via entrapment of an early neural lineage network, which represents a critical therapeutic vulnerability for this recalcitrant disease.

A C19MC-LIN28A-MYCN Oncogenic Circuit Driven by Hijacked Super-enhancers Is a Distinct Therapeutic Vulnerability in ETMRs: A Lethal Brain Tumor

Embryonal tumors with multilayered rosettes (ETMRs) are highly lethal infant brain cancers with characteristic amplification of Chr19q13.41 miRNA cluster (C19MC) and enrichment of pluripotency factor LIN28A. Here we investigated C19MC oncogenic mechanisms and discovered a C19MC-LIN28A-MYCN circuit fueled by multiple complex regulatory loops including an MYCN core transcriptional network and super-enhancers resulting from long-range MYCN DNA interactions and C19MC gene fusions. Our data show that this powerful oncogenic circuit, which entraps an early neural lineage network, is potently abrogated by bromodomain inhibitor JQ1, leading to ETMR cell death.

CAPTURE: In Situ Analysis of Chromatin Composition of Endogenous Genomic Loci by Biotinylated dCas9

Cis-regulatory elements (CREs) play a pivotal role in spatiotemporal control of tissue-specific gene expression, yet the molecular composition of the vast majority of CREs in native chromatin remains unknown. In this article, we describe the clustered regularly interspaced short palindromic repeats (CRISPR) affinity purification in situ of regulatory elements (CAPTURE) approach to simultaneously identify locus-specific chromatin-regulating protein complexes and long-range DNA interactions. Using an in vivo biotinylated nuclease-deficient Cas9 (dCas9) protein and programmable single guide RNAs (sgRNAs), this approach allows for high-resolution and locus-specific isolation of protein complexes and long-range chromatin looping associated with single copy CREs in mammalian cells. Unbiased analysis of the compositional structure of developmentally regulated or disease-associated CREs identifies new features of transcriptional regulation. Hence, CAPTURE provides a versatile platform to study genomic locus-regulating chromatin composition in a mammalian genome. © 2018 by John Wiley & Sons, Inc.

In Situ Capture of Chromatin Interactions by Biotinylated dCas9

Cis-regulatory elements (CREs) are commonly recognized by correlative chromatin features, yet the molecular composition of the vast majority of CREs inchromatin remains unknown. Here, we describe aCRISPR affinity purification in situ of regulatory elements (CAPTURE) approach to unbiasedly identifylocus-specific chromatin-regulating protein complexes and long-range DNA interactions. Using an invivo biotinylated nuclease-deficient Cas9 protein and sequence-specific guide RNAs, we show high-resolution and selective isolation of chromatin interactions at a single-copy genomic locus. Purificationof human telomeres using CAPTURE identifies known and new telomeric factors. In situ capture of individual constituents of the enhancer cluster controlling human β-globin genes establishes evidence for composition-based hierarchical organization. Furthermore, unbiased analysis of chromatin interactions at disease-associated cis-elements and developmentally regulated super-enhancers reveals spatial features that causally control gene transcription. Thus, comprehensive and unbiased analysis oflocus-specific regulatory composition provides mechanistic insight into genome structure and function in development and disease.

Binding, sliding, and function of cohesin during transcriptional activation

The ring-shaped cohesin complex orchestrates long-range DNA interactions to mediate sister chromatid cohesion and other aspects of chromosome structure and function. In the yeast Saccharomyces cerevisiae, the complex binds discrete sites along chromosomes, including positions within and around genes. Transcriptional activity redistributes the complex to the 3' ends of convergently oriented gene pairs. Despite the wealth of information about where cohesin binds, little is known about cohesion at individual chromosomal binding sites and how transcription affects cohesion when cohesin complexes redistribute. In this study, we generated extrachromosomal DNA circles to study cohesion in response to transcriptional induction of a model gene, URA3. Functional cohesin complexes loaded onto the locus via a poly(dA:dT) tract in the gene promoter and mediated cohesion before induction. Upon transcription, the fate of these complexes depended on whether the DNA was circular or not. When gene activation occurred before DNA circularization, cohesion was lost. When activation occurred after DNA circularization, cohesion persisted. The presence of a convergently oriented gene also prevented transcription-driven loss of functional cohesin complexes, at least in M phase-arrested cells. The results are consistent with cohesin binding chromatin in a topological embrace and with transcription mobilizing functional complexes by sliding them along DNA.

Quantitative BrdU immunoprecipitation method demonstrates that Fkh1 and Fkh2 are rate-limiting activators of replication origins that reprogram replication timing in G1 phase.

The Saccharomyces cerevisiae Forkhead Box (FOX) proteins, Fkh1 and Fkh2, regulate diverse cellular processes including transcription, long-range DNA interactions during homologous recombination, and replication origin timing and long-range origin clustering. We hypothesized that, as stimulators of early origin activation, Fkh1 and Fkh2 abundance limits the rate of origin activation genome-wide. Existing methods, however, are not well-suited to quantitative, genome-wide measurements of origin firing between strains and conditions. To overcome this limitation, we developed qBrdU-seq, a quantitative method for BrdU incorporation analysis of replication dynamics, and applied it to show that overexpression of Fkh1 and Fkh2 advances the initiation timing of many origins throughout the genome resulting in a higher total level of origin initiations in early S phase. The higher initiation rate is accompanied by slower replication fork progression, thereby maintaining a normal length of S phase without causing detectable Rad53 checkpoint kinase activation. The advancement of origin firing time, including that of origins in heterochromatic domains, was established in late G1 phase, indicating that origin timing can be reset subsequently to origin licensing. These results provide novel insights into the mechanisms of origin timing regulation by identifying Fkh1 and Fkh2 as rate-limiting factors for origin firing that determine the ability of replication origins to accrue limiting factors and have the potential to reprogram replication timing late in G1 phase.

RNA Exosome-Regulated Long Non-Coding RNA Transcription Controls Super-Enhancer Activity

We have ablated the cellular RNA degradation machinery in differentiated B cells and pluripotent embryonic stem cells (ESCs) by conditional mutagenesis of core (Exosc3) and nuclear RNase (Exosc10) components of RNA exosome and identified a vast number of long non-coding RNAs (lncRNAs) and enhancer RNAs (eRNAs) with emergent functionality. Unexpectedly, eRNA-expressing regions accumulate R-loop structures upon RNA exosome ablation, thus demonstrating the role of RNA exosome in resolving deleterious DNA/RNA hybrids arising from active enhancers. We have uncovered a distal divergent eRNA-expressing element (lncRNA-CSR) engaged in long-range DNA interactions and regulating IgH 3' regulatory region super-enhancer function. CRISPR-Cas9-mediated ablation of lncRNA-CSR transcription decreases its chromosomal looping-mediated association with the IgH 3' regulatory region super-enhancer and leads to decreased class switch recombination efficiency. We propose that the RNA exosome protects divergently transcribed lncRNA expressing enhancers by resolving deleterious transcription-coupled secondary DNA structures, while also regulating long-range super-enhancer chromosomal interactions important for cellular function.

Chromatin insulator factors involved in long-range DNA interactions and their role in the folding of the Drosophila genome.

Chromatin insulators are genetic elements implicated in the organization of chromatin and the regulation of transcription. In Drosophila, different insulator types were characterized by their locus-specific composition of insulator proteins and co-factors. Insulators mediate specific long-range DNA contacts required for the three dimensional organization of the interphase nucleus and for transcription regulation, but the mechanisms underlying the formation of these contacts is currently unknown. Here, we investigate the molecular associations between different components of insulator complexes (BEAF32, CP190 and Chromator) by biochemical and biophysical means, and develop a novel single-molecule assay to determine what factors are necessary and essential for the formation of long-range DNA interactions. We show that BEAF32 is able to bind DNA specifically and with high affinity, but not to bridge long-range interactions (LRI). In contrast, we show that CP190 and Chromator are able to mediate LRI between specifically-bound BEAF32 nucleoprotein complexes in vitro. This ability of CP190 and Chromator to establish LRI requires specific contacts between BEAF32 and their C-terminal domains, and dimerization through their N-terminal domains. In particular, the BTB/POZ domains of CP190 form a strict homodimer, and its C-terminal domain interacts with several insulator binding proteins. We propose a general model for insulator function in which BEAF32/dCTCF/Su(HW) provide DNA specificity (first layer proteins) whereas CP190/Chromator are responsible for the physical interactions required for long-range contacts (second layer). This network of organized, multi-layer interactions could explain the different activities of insulators as chromatin barriers, enhancer blockers, and transcriptional regulators, and suggest a general mechanism for how insulators may shape the organization of higher-order chromatin during cell division.

A novel antisense long noncoding RNA within the IGF1R gene locus is imprinted in hematopoietic malignancies.

Dysregulation of the insulin-like growth factor type I receptor (IGF1R) has been implicated in the progression and therapeutic resistance of malignancies. In acute myeloid leukemia (AML) cells, IGF1R is one of the most abundantly phosphorylated receptor tyrosine kinases, promoting cell growth through the PI3K/Akt signaling pathway. However, little is known regarding the molecular mechanisms underlying IGF1R gene dysregulation in cancer. We discovered a novel intragenic long noncoding RNA (lncRNA) within the IGF1R locus, named IRAIN, which is transcribed in an antisense direction from an intronic promoter. The IRAIN lncRNA was expressed exclusively from the paternal allele, with the maternal counterpart being silenced. Using both reverse transcription-associated trap and chromatin conformation capture assays, we demonstrate that this lncRNA interacts with chromatin DNA and is involved in the formation of an intrachromosomal enhancer/promoter loop. Knockdown of IRAIN lncRNA with shRNA abolishes this intrachromosomal interaction. In addition, IRAIN was downregulated both in leukemia cell lines and in blood obtained from high-risk AML patients. These data identify IRAIN as a new imprinted lncRNA that is involved in long-range DNA interactions.

Long-range DNA interactions at the IL-1/IL-36/IL-37 gene cluster (2q13) are induced by activation of monocytes

The interleukin-1 gene cluster occupies a 360kb region of chromosome 2q13 and contains nine homologous genes. These include agonists and antagonists of the parallel IL-1 and IL-36 systems, and IL1F7, the gene encoding IL-37. As the genes of the cluster are structurally and functionally related and have similar mRNA kinetics, we have sought evidence for gene induction-specific looping of chromatin in the IL-1 cluster by chromatin conformation capture (3C).We show here that IL1A, IL1B and IL1F7 regulatory regions come in close proximity in LPS stimulated cells but not in resting human monocytes. This suggests that IL1A, IL1B and IL1F7 are likely transcribed by the same transcription factory.One cardinal function of transcriptional Locus Control Region (LCR) is bringing map-distant activated genes into close physical proximity within the transcription factory. Our data show distant intergenic DNA segments are also in close proximity to the regulatory regions of the three genes. This may indicate that they are co-regulated and raise the possibility of a LCR within the cluster.

CTCF and Rad21 Act as Host Cell Restriction Factors for Kaposi's Sarcoma-Associated Herpesvirus (KSHV) Lytic Replication by Modulating Viral Gene Transcription

Kaposi's sarcoma-associated herpesvirus (KSHV) is a human herpesvirus that causes Kaposi's sarcoma and is associated with the development of lymphoproliferative diseases. KSHV reactivation from latency and virion production is dependent on efficient transcription of over eighty lytic cycle genes and viral DNA replication. CTCF and cohesin, cellular proteins that cooperatively regulate gene expression and mediate long-range DNA interactions, have been shown to bind at specific sites in herpesvirus genomes. CTCF and cohesin regulate KSHV gene expression during latency and may also control lytic reactivation, although their role in lytic gene expression remains incompletely characterized. Here, we analyze the dynamic changes in CTCF and cohesin binding that occur during the process of KSHV viral reactivation and virion production by high resolution chromatin immunoprecipitation and deep sequencing (ChIP-Seq) and show that both proteins dissociate from viral genomes in kinetically and spatially distinct patterns. By utilizing siRNAs to specifically deplete CTCF and Rad21, a cohesin component, we demonstrate that both proteins are potent restriction factors for KSHV replication, with cohesin knockdown leading to hundred-fold increases in viral yield. High-throughput RNA sequencing was used to characterize the transcriptional effects of CTCF and cohesin depletion, and demonstrated that both proteins have complex and global effects on KSHV lytic transcription. Specifically, both proteins act as positive factors for viral transcription initially but subsequently inhibit KSHV lytic transcription, such that their net effect is to limit KSHV RNA accumulation. Cohesin is a more potent inhibitor of KSHV transcription than CTCF but both proteins are also required for efficient transcription of a subset of KSHV genes. These data reveal novel effects of CTCF and cohesin on transcription from a relatively small genome that resemble their effects on the cellular genome by acting as gene-specific activators of some promoters, but differ in acting as global negative regulators of transcription.

Long‐range transcriptional regulation of breast cancer genes

Breast cancer is a major health problem and understanding the genetic basis of this disease is crucial for predicting risk and developing effective targeted therapeutics. Several breast cancer predisposing genes have been identified, but mutations in the coding regions of these genes only accounts for a small proportion of risk. Research now suggests that combinations of multiple non-coding changes in breast cancer susceptibility genes, which cause moderate alterations in gene expression, will be responsible for the remaining inherited risk. These non-coding changes will include variants in proximal and distal transcriptional and post-transcriptional regulatory elements and may affect the levels and function of trans-acting factors, including proteins and RNAs, which act on these elements. Somatic changes in such elements and factors have also been associated with breast cancer progression. With the recent advent of techniques allowing the detection of long-range DNA interactions spanning the human genome, it has become increasingly clear that long-range regulatory elements constitute an important mechanism for gene regulation. Recent studies have identified several such elements that are important for regulating genes involved in breast cancer, raising the possibility that defects in these sequences may contribute to breast cancer predisposition and progression. In this review, we discuss the emerging functions of cis-regulatory elements and a subset of trans-acting factors in breast tumorigenesis. We also discuss some recent progress in our understanding of how dysregulation in these transcriptional components may contribute to breast cancer, and the potential implications for molecular diagnosis, prognosis prediction, and the treatment of this disease.

DNA replication timing and long-range DNA interactions predict mutational landscapes of cancer genomes

Somatic copy-number alterations (SCNA) are a hallmark of many cancer types, but the mechanistic basis underlying their genome-wide patterns remains incompletely understood. Here we integrate data on DNA replication timing, long-range interactions between genomic material, and 331,724 SCNAs from 2,792 cancer samples classified into 26 cancer types. We report that genomic regions of similar replication timing are clustered spatially in the nucleus, that the two boundaries of SCNAs tend to be found in such regions, and that regions replicated early and late display distinct patterns of frequencies of SCNA boundaries, SCNA size and a preference for deletions over insertions. We show that long-range interaction and replication timing data alone can identify a significant proportion of SCNAs in an independent test data set. We propose a model for the generation of SCNAs in cancer, suggesting that data on spatial proximity of regions replicating at the same time can be used to predict the mutational landscapes of cancer genomes.