Dynamic enhancer-gene body contacts during transcription elongation

Resolving the three-dimensional interactome of Human Accelerated Regions during human and chimpanzee neurodevelopment

Human Accelerated Regions (HARs) are highly conserved across species but exhibit a significant excess of human-specific sequence changes, suggesting they may have gained novel functions in human evolution. HARs include transcriptional enhancers with human-specific activity and have been implicated in the evolution of the human brain. However, our understanding of how HARs contributed to uniquely human features of the brain is hindered by a lack of insight into the genes and pathways that HARs regulate. It is unclear whether HARs acted by altering the expression of gene targets conserved between HARs and their chimpanzee orthologs or by gaining new gene targets in human, a mechanism termed enhancer hijacking. We generated a high-resolution map of chromatin interactions for 1,590 HARs and their orthologs in human and chimpanzee neural stem cells (NSCs) to comprehensively identify gene targets in both species. HARs and their chimpanzee orthologs targeted a conserved set of 2,963 genes enriched for neurodevelopmental processes including neurogenesis and synaptic transmission. Changes in HAR enhancer activity were correlated with changes in conserved gene target expression. Conserved targets were enriched among genes differentially expressed between human and chimpanzee NSCs or between human and non-human primate developing and adult brain. Species-specific HAR gene targets did not converge on known biological functions and were not significantly enriched among differentially expressed genes, suggesting that HARs did not alter gene expression via enhancer hijacking. HAR gene targets, including differentially expressed targets, also showed cell type-specific expression patterns in the developing human brain, including outer radial glia, which are hypothesized to contribute to human cortical expansion. Our findings support that HARs influenced human brain evolution by altering the expression of conserved gene targets and provide the means to functionally link HARs with novel human brain features.

Physical and functional interaction among Irf8 enhancers during dendritic cell differentiation

The production of type 1 conventional dendritic cells (cDC1s) requires high expression of the transcription factor IRF8. Three enhancers at the Irf8 3' region function in a differentiation stage-specific manner. However, whether and how these enhancers interact physically and functionally remains unclear. Here, we show that the Irf8 3' enhancers directly interact with each other and contact the Irf8 gene body during cDC1 differentiation. The +56 kb enhancer, which functions from multipotent progenitor stages, activates the other 3' enhancers through an IRF8-dependent transcription factor program, that is, in trans. Then, the +32 kb enhancer, which operates in cDC1-committed cells, reversely acts in cis on the other 3' enhancers to maintain the high expression of Irf8. Indeed, mice with compound heterozygous deletion of the +56 and +32 kb enhancers are unable to generate cDC1s. These results illustrate how multiple enhancers cooperate to induce a lineage-determining transcription factor gene during cell differentiation.

Integrative analysis of transcriptomic and epigenomic data reveals distinct patterns for developmental and housekeeping gene regulation

BACKGROUND

Regulation of transcription is central to the emergence of new cell types during development, and it often involves activation of genes via proximal and distal regulatory regions. The activity of regulatory elements is determined by transcription factors (TFs) and epigenetic marks, but despite extensive mapping of such patterns, the extraction of regulatory principles remains challenging.

RESULTS

Here we study differentially and similarly expressed genes along with their associated epigenomic profiles, chromatin accessibility and DNA methylation, during lineage specification at gastrulation in mice. Comparison of the three lineages allows us to identify genomic and epigenomic features that distinguish the two classes of genes. We show that differentially expressed genes are primarily regulated by distal elements, while similarly expressed genes are controlled by proximal housekeeping regulatory programs. Differentially expressed genes are relatively isolated within topologically associated domains, while similarly expressed genes tend to be located in gene clusters. Transcription of differentially expressed genes is associated with differentially open chromatin at distal elements including enhancers, while that of similarly expressed genes is associated with ubiquitously accessible chromatin at promoters.

CONCLUSION

Based on these associations of (linearly) distal genes' transcription start sites (TSSs) and putative enhancers for developmental genes, our findings allow us to link putative enhancers to their target promoters and to infer lineage-specific repertoires of putative driver transcription factors, within which we define subgroups of pioneers and co-operators.

Multicenter integrated analysis of noncoding CRISPRi screens

The ENCODE Consortium's efforts to annotate noncoding cis-regulatory elements (CREs) have advanced our understanding of gene regulatory landscapes. Pooled, noncoding CRISPR screens offer a systematic approach to investigate cis-regulatory mechanisms. The ENCODE4 Functional Characterization Centers conducted 108 screens in human cell lines, comprising >540,000 perturbations across 24.85 megabases of the genome. Using 332 functionally confirmed CRE-gene links in K562 cells, we established guidelines for screening endogenous noncoding elements with CRISPR interference (CRISPRi), including accurate detection of CREs that exhibit variable, often low, transcriptional effects. Benchmarking five screen analysis tools, we find that CASA produces the most conservative CRE calls and is robust to artifacts of low-specificity single guide RNAs. We uncover a subtle DNA strand bias for CRISPRi in transcribed regions with implications for screen design and analysis. Together, we provide an accessible data resource, predesigned single guide RNAs for targeting 3,275,697 ENCODE SCREEN candidate CREs with CRISPRi and screening guidelines to accelerate functional characterization of the noncoding genome.

Transcription induces context-dependent remodeling of chromatin architecture during differentiation

Metazoan chromosomes are organized into discrete spatial domains (TADs), believed to contribute to the regulation of transcriptional programs. Despite extensive correlation between domain organization and gene activity, a direct mechanistic link is unclear, with perturbation studies often showing little effect. To follow chromatin architecture changes during development, we used Capture Hi-C to interrogate the domains around key differentially expressed genes during mouse thymocyte maturation, uncovering specific remodeling events. Notably, one TAD boundary was broadened to accommodate RNA polymerase elongation past the border, and subdomains were formed around some activated genes without changes in CTCF binding. The ectopic induction of some genes was sufficient to recapitulate domain formation in embryonic stem cells, providing strong evidence that transcription can directly remodel chromatin structure. These results suggest that transcriptional processes drive complex chromosome folding patterns that can be important in certain genomic contexts.

High-sensitive nascent transcript sequencing reveals BRD4-specific control of widespread enhancer and target gene transcription

Gene transcription by RNA polymerase II (Pol II) is under control of promoters and distal regulatory elements known as enhancers. Enhancers are themselves transcribed by Pol II correlating with their activity. How enhancer transcription is regulated and coordinated with transcription at target genes has remained unclear. Here, we developed a high-sensitive native elongating transcript sequencing approach, called HiS-NET-seq, to provide an extended high-resolution view on transcription, especially at lowly transcribed regions such as enhancers. HiS-NET-seq uncovers new transcribed enhancers in human cells. A multi-omics analysis shows that genome-wide enhancer transcription depends on the BET family protein BRD4. Specifically, BRD4 co-localizes to enhancer and promoter-proximal gene regions, and is required for elongation activation at enhancers and their genes. BRD4 keeps a set of enhancers and genes in proximity through long-range contacts. From these studies BRD4 emerges as a general regulator of enhancer transcription that may link transcription at enhancers and genes.

Glucocorticoid stimulation induces regionalized gene responses within topologically associating domains

Transcription-factor binding to cis-regulatory regions regulates the gene expression program of a cell, but occupancy is often a poor predictor of the gene response. Here, we show that glucocorticoid stimulation led to the reorganization of transcriptional coregulators MED1 and BRD4 within topologically associating domains (TADs), resulting in active or repressive gene environments. Indeed, we observed a bias toward the activation or repression of a TAD when their activities were defined by the number of regions gaining and losing MED1 and BRD4 following dexamethasone (Dex) stimulation. Variations in Dex-responsive genes at the RNA levels were consistent with the redistribution of MED1 and BRD4 at the associated cis-regulatory regions. Interestingly, Dex-responsive genes without the differential recruitment of MED1 and BRD4 or binding by the glucocorticoid receptor were found within TADs, which gained or lost MED1 and BRD4, suggesting a role of the surrounding environment in gene regulation. However, the amplitude of the response of Dex-regulated genes was higher when the differential recruitment of the glucocorticoid receptor and transcriptional coregulators was observed, reaffirming the role of transcription factor-driven gene regulation and attributing a lesser role to the TAD environment. These results support a model where a signal-induced transcription factor induces a regionalized effect throughout the TAD, redefining the notion of direct and indirect effects of transcription factors on target genes.

Recent advances in chromosome capture techniques unraveling 3D genome architecture in germ cells, health, and disease

In eukaryotes, the genome does not emerge in a specific shape but rather as a hierarchial bundle within the nucleus. This multifaceted genome organization consists of multiresolution cellular structures, such as chromosome territories, compartments, and topologically associating domains, which are frequently defined by architecture, design proteins including CTCF and cohesin, and chromatin loops. This review briefly discusses the advances in understanding the basic rules of control, chromatin folding, and functional areas in early embryogenesis. With the use of chromosome capture techniques, the latest advancements in technologies for visualizing chromatin interactions come close to revealing 3D genome formation frameworks with incredible detail throughout all genomic levels, including at single-cell resolution. The possibility of detecting variations in chromatin architecture might open up new opportunities for disease diagnosis and prevention, infertility treatments, therapeutic approaches, desired exploration, and many other application scenarios.

Transcription shapes 3D chromatin organization by interacting with loop extrusion

Cohesin folds mammalian interphase chromosomes by extruding the chromatin fiber into numerous loops. "Loop extrusion" can be impeded by chromatin-bound factors, such as CTCF, which generates characteristic and functional chromatin organization patterns. It has been proposed that transcription relocalizes or interferes with cohesin and that active promoters are cohesin loading sites. However, the effects of transcription on cohesin have not been reconciled with observations of active extrusion by cohesin. To determine how transcription modulates extrusion, we studied mouse cells in which we could alter cohesin abundance, dynamics, and localization by genetic "knockouts" of the cohesin regulators CTCF and Wapl. Through Hi-C experiments, we discovered intricate, cohesin-dependent contact patterns near active genes. Chromatin organization around active genes exhibited hallmarks of interactions between transcribing RNA polymerases (RNAPs) and extruding cohesins. These observations could be reproduced by polymer simulations in which RNAPs were moving barriers to extrusion that obstructed, slowed, and pushed cohesins. The simulations predicted that preferential loading of cohesin at promoters is inconsistent with our experimental data. Additional ChIP-seq experiments showed that the putative cohesin loader Nipbl is not predominantly enriched at promoters. Therefore, we propose that cohesin is not preferentially loaded at promoters and that the barrier function of RNAP accounts for cohesin accumulation at active promoters. Altogether, we find that RNAP is an extrusion barrier that is not stationary, but rather, translocates and relocalizes cohesin. Loop extrusion and transcription might interact to dynamically generate and maintain gene interactions with regulatory elements and shape functional genomic organization.

Gene silencing dynamics are modulated by transiently active regulatory elements

Transcriptional enhancers have been extensively characterized, but cis-regulatory elements involved in acute gene repression have received less attention. Transcription factor GATA1 promotes erythroid differentiation by activating and repressing distinct gene sets. Here, we study the mechanism by which GATA1 silences the proliferative gene Kit during murine erythroid cell maturation and define stages from initial loss of activation to heterochromatinization. We find that GATA1 inactivates a potent upstream enhancer but concomitantly creates a discrete intronic regulatory region marked by H3K27ac, short noncoding RNAs, and de novo chromatin looping. This enhancer-like element forms transiently and serves to delay Kit silencing. The element is ultimately erased via the FOG1/NuRD deacetylase complex, as revealed by the study of a disease-associated GATA1 variant. Hence, regulatory sites can be self-limiting by dynamic co-factor usage. Genome-wide analyses across cell types and species uncover transiently active elements at numerous genes during repression, suggesting that modulation of silencing kinetics is widespread.

Double-strand break repair and mis-repair in 3D

DNA double-strand breaks (DSBs) are lesions that arise frequently from exposure to damaging agents as well as from ongoing physiological DNA transactions. Mis-repair of DSBs leads to rearrangements and structural variations in chromosomes, including insertions, deletions, and translocations implicated in disease. The DNA damage response (DDR) limits pathologic mutations and large-scale chromosome rearrangements. DSB repair initiates in 2D at DNA lesions with the stepwise recruitment of repair proteins and local chromatin remodeling which facilitates break accessibility. More complex structures are then formed via protein assembly into nanodomains and via genome folding into chromatin loops. Subsequently, 3D reorganization of DSBs is guided by clustering forces which drive the assembly of repair domains harboring multiple lesions. These domains are further stabilized and insulated into condensates via liquid-liquid phase-separation. Here, we discuss the benefits and risks associated with this 3D reorganization of the broken genome.

Erythroid Cell Research: 3D Chromatin, Transcription Factors and Beyond

Studies of the regulatory networks and signals controlling erythropoiesis have brought important insights in several research fields of biology and have been a rich source of discoveries with far-reaching implications beyond erythroid cells biology. The aim of this review is to highlight key recent discoveries and show how studies of erythroid cells bring forward novel concepts and refine current models related to genome and 3D chromatin organization, signaling and disease, with broad interest in life sciences.

Spatial organization of transcribed eukaryotic genes

Despite the well-established role of nuclear organization in the regulation of gene expression, little is known about the reverse: how transcription shapes the spatial organization of the genome. Owing to the small sizes of most previously studied genes and the limited resolution of microscopy, the structure and spatial arrangement of a single transcribed gene are still poorly understood. Here we study several long highly expressed genes and demonstrate that they form open-ended transcription loops with polymerases moving along the loops and carrying nascent RNAs. Transcription loops can span across micrometres, resembling lampbrush loops and polytene puffs. The extension and shape of transcription loops suggest their intrinsic stiffness, which we attribute to decoration with multiple voluminous nascent ribonucleoproteins. Our data contradict the model of transcription factories and suggest that although microscopically resolvable transcription loops are specific for long highly expressed genes, the mechanisms underlying their formation could represent a general aspect of eukaryotic transcription.

MED12 and BRD4 cooperate to sustain cancer growth upon loss of mediator kinase

Mediator kinases (CDK8/19) are transcriptional regulators broadly implicated in cancer. Despite their central role in fine-tuning gene-expression programs, we find complete loss of CDK8/19 is tolerated in colorectal cancer (CRC) cells. Using orthogonal functional genomic and pharmacological screens, we identify BET protein inhibition as a distinct vulnerability in CDK8/19-depleted cells. Combined CDK8/19 and BET inhibition led to synergistic growth retardation in human and mouse models of CRC. Strikingly, depletion of CDK8/19 in these cells led to global repression of RNA polymerase II (Pol II) promoter occupancy and transcription. Concurrently, loss of Mediator kinase led to a profound increase in MED12 and BRD4 co-occupancy at enhancer elements and increased dependence on BET proteins for the transcriptional output of cell-essential genes. In total, this work demonstrates a synthetic lethal interaction between Mediator kinase and BET proteins and exposes a therapeutic vulnerability that can be targeted using combination therapies.

Exploring the role of eRNA in regulating gene expression

eRNAs as the products of enhancers can regulate gene expression via various possible ways, but which regulation way is more reasonable is debatable in biology, and in particular, how eRNAs impact gene expression remains unclear. Here we introduce a mechanistic model of gene expression to address these issues. This model considers three possible regulation ways of eRNA: Type-I by which eRNA regulates transcriptional activity by facilitating the formation of enhancer-promoter (E-P) loop, Type-II by which eRNA directly promotes the mRNA production rate, and mixed regulation (i.e., the combination of Type-I and Type-II). We show that with the increase of the E-P loop length, mRNA distribution can transition from unimodality to bimodality or vice versa in all the three regulation cases. However, in contrast to the other two regulations, Type-II regulation can lead to the highest mean mRNA level and the lowest mRNA noise, independent of the E-P loop length. These results would not only reveal the essential mechanism of how eRNA regulates gene expression, but also imply a new mechanism for phenotypic switching, namely the E-P loop can induce phenotypic switching.

Primordial super-enhancers: heat shock-induced chromatin organization in yeast

Specialized mechanisms ensure proper expression of critically important genes such as those specifying cell identity or conferring protection from environmental stress. Investigations of the heat shock response have been critical in elucidating basic concepts of transcriptional control. Recent studies demonstrate that in response to thermal stress, heat shock-responsive genes associate with high levels of transcriptional activators and coactivators and those in yeast intensely interact across and between chromosomes, coalescing into condensates. In mammalian cells, cell identity genes that are regulated by super-enhancers (SEs) are also densely occupied by transcriptional machinery that form phase-separated condensates. We suggest that the stress-remodeled yeast nucleome bears functional and structural resemblance to mammalian SEs, and will reveal fundamental mechanisms of gene control by transcriptional condensates.

Chromatin Conformation in Development and Disease

Chromatin domains and loops are important elements of chromatin structure and dynamics, but much remains to be learned about their exact biological role and nature. Topological associated domains and functional loops are key to gene expression and hold the answer to many questions regarding developmental decisions and diseases. Here, we discuss new findings, which have linked chromatin conformation with development, differentiation and diseases and hypothesized on various models while integrating all recent findings on how chromatin architecture affects gene expression during development, evolution and disease.

Collaboration through chromatin: motors of transcription and chromatin structure

Packaging of the eukaryotic genome into chromatin places fundamental physical constraints on transcription. Clarifying how transcription operates within these constraints is essential to understand how eukaryotic gene expression programs are established and maintained. Here we review what is known about the mechanisms of transcription on chromatin templates. Current models indicate that transcription through chromatin is accomplished by the combination of an inherent nucleosome disrupting activity of RNA polymerase and the action of ATP-dependent chromatin remodeling motors. Collaboration between these two types of molecular motors is proposed to occur at all stages of transcription through diverse mechanisms. Further investigation of how these two motors combine their basic activities is essential to clarify the interdependent relationship between genome structure and transcription.

Super-enhancer Acquisition Drives FOXC2 Expression in Middle Ear Cholesteatoma

Distinct histone modifications regulate gene expression in certain diseases, but little is known about histone epigenetics in middle ear cholesteatoma. It is known that histone acetylation destabilizes the nucleosome and chromatin structure and induces gene activation. The association of histone acetylation with chronic inflammatory diseases has been indicated in recent studies. In this study, we examined the localization of variously modified histone H3 acetylation at lysine 9, 14, 18, 23, and 27 in paraffin-embedded sections of human middle ear cholesteatoma (cholesteatoma) tissues and the temporal bones of an animal model of cholesteatoma immunohistochemically. As a result, we found that there was a significant increase of the expression levels of H3K27ac both in human cholesteatoma tissues and the animal model. In genetics, super-enhancers are clusters of enhancers that drive the transcription of genes involved in cell identity. Super-enhancers were originally defined using the H3K27ac signal, and then we used H3K27ac chromatin immunoprecipitation followed by sequencing to map the active cis-regulatory landscape in human cholesteatoma. Based on the results, we identified increased H3K27ac signals as super-enhancers of the FOXC2 loci, as well as increased protein of FOXC2 in cholesteatoma. Recent studies have indicated that menin-MLL inhibitor could suppress tumor growth through the control of histone H3 modification. In this study, we demonstrated that the expression of FOXC2 was inhibited by menin-MLL inhibitor in vivo. These findings indicate that FOXC2 expression under histone modifications promoted the pathogenesis of cholesteatoma and suggest that it may be a therapeutic target of cholesteatoma.

Mechanisms of enhancer action: the known and the unknown

Differential gene expression mechanisms ensure cellular differentiation and plasticity to shape ontogenetic and phylogenetic diversity of cell types. A key regulator of differential gene expression programs are the enhancers, the gene-distal cis-regulatory sequences that govern spatiotemporal and quantitative expression dynamics of target genes. Enhancers are widely believed to physically contact the target promoters to effect transcriptional activation. However, our understanding of the full complement of regulatory proteins and the definitive mechanics of enhancer action is incomplete. Here, we review recent findings to present some emerging concepts on enhancer action and also outline a set of outstanding questions.

EnhFFL: A database of enhancer mediated feed-forward loops for human and mouse

Feed-forward loops (FFLs) are thought to be one of the most common and important classes of transcriptional network motifs involved in various diseases. Enhancers are cis-regulatory elements that positively regulate protein-coding genes or microRNAs (miRNAs) by recruiting DNA-binding transcription factors (TFs). However, a comprehensive resource to identify, store, and analyze the FFLs of typical enhancer and super-enhancer FFLs is not currently available. Here, we present EnhFFL, an online database to provide a data resource for users to browse and search typical enhancer and super-enhancer FFLs. The current database covers 46 280/7000 TF-enhancer-miRNA FFLs, 9997/236 enhancer-miRNA-gene FFLs, 3 561 164/3 193 182 TF-enhancer-gene FFLs, and 1259/235 TF-enhancer feed-back loops (FBLs) across 91 tissues/cell lines of human and mouse, respectively. Users can browse loops by selecting species, types of tissue/cell line, and types of FFLs. EnhFFL supports searching elements including name/ID, genomic location, and the conservation of miRNA target genes. We also developed tools for users to screen customized FFLs using the threshold of q value as well as the confidence score of miRNA target genes. Disease and functional enrichment analysis showed that master miRNAs that are widely engaged in FFLs including TF-enhancer-miRNAs and enhancer-miRNA-genes are significantly involved in tumorigenesis. Database URL:http://lcbb.swjtu.edu.cn/EnhFFL/.

Enhancers predominantly regulate gene expression during differentiation via transcription initiation

Gene transcription occurs via a cycle of linked events, including initiation, promoter-proximal pausing, and elongation of RNA polymerase II (Pol II). A key question is how transcriptional enhancers influence these events to control gene expression. Here, we present an approach that evaluates the level and change in promoter-proximal transcription (initiation and pausing) in the context of differential gene expression, genome-wide. This combinatorial approach shows that in primary cells, control of gene expression during differentiation is achieved predominantly via changes in transcription initiation rather than via release of Pol II pausing. Using genetically engineered mouse models, deleted for functionally validated enhancers of the α- and β-globin loci, we confirm that these elements regulate Pol II recruitment and/or initiation to modulate gene expression. Together, our data show that gene expression during differentiation is regulated predominantly at the level of initiation and that enhancers are key effectors of this process.

Enhancer-Gene Interaction Analyses Identified the Epidermal Growth Factor Receptor as a Susceptibility Gene for Type 2 Diabetes Mellitus

Background

Genetic interactions are known to play an important role in the missing heritability problem for type 2 diabetes mellitus (T2DM). Interactions between enhancers and their target genes play important roles in gene regulation and disease pathogenesis. In the present study, we aimed to identify genetic interactions between enhancers and their target genes associated with T2DM.

Methods

We performed genetic interaction analyses of enhancers and protein-coding genes for T2DM in 2,696 T2DM patients and 3,548 controls of European ancestry. A linear regression model was used to identify single nucleotide polymorphism (SNP) pairs that could affect the expression of the protein-coding genes. Differential expression analyses were used to identify differentially expressed susceptibility genes in diabetic and nondiabetic subjects.

Results

We identified one SNP pair, rs4947941×rs7785013, significantly associated with T2DM (combined P=4.84×10-10). The SNP rs4947941 was annotated as an enhancer, and rs7785013 was located in the epidermal growth factor receptor (EGFR) gene. This SNP pair was significantly associated with EGFR expression in the pancreas (P=0.033), and the minor allele "A" of rs7785013 decreased EGFR gene expression and the risk of T2DM with an increase in the dosage of "T" of rs4947941. EGFR expression was significantly upregulated in T2DM patients, which was consistent with the effect of rs4947941×rs7785013 on T2DM and EGFR expression. A functional validation study using the Mouse Genome Informatics (MGI) database showed that EGFR was associated with diabetes-relevant phenotypes.

Conclusion

Genetic interaction analyses of enhancers and protein-coding genes suggested that EGFR may be a novel susceptibility gene for T2DM.

Molecular Basis of the Function of Transcriptional Enhancers

Transcriptional enhancers are major genomic elements that control gene activity in eukaryotes. Recent studies provided deeper insight into the temporal and spatial organization of transcription in the nucleus, the role of non-coding RNAs in the process, and the epigenetic control of gene expression. Thus, multiple molecular details of enhancer functioning were revealed. Here, we describe the recent data and models of molecular organization of enhancer-driven transcription.

Genome-wide analyses of chromatin interactions after the loss of Pol I, Pol II, and Pol III

BACKGROUND

The relationship between transcription and the 3D chromatin structure is debated. Multiple studies have shown that transcription affects global Cohesin binding and 3D genome structures. However, several other studies have indicated that inhibited transcription does not alter chromatin conformations.

RESULTS

We provide the most comprehensive evidence to date to demonstrate that transcription plays a relatively modest role in organizing the local, small-scale chromatin structures in mammalian cells. We show degraded Pol I, Pol II, and Pol III proteins in mESCs cause few or no changes in large-scale 3D chromatin structures, selected RNA polymerases with a high abundance of binding sites or active promoter-associated interactions appear to be relatively more affected after the degradation, transcription inhibition alters local, small loop domains, as indicated by high-resolution chromatin interaction maps, and loops with bound Pol II but without Cohesin or CTCF are identified and found to be largely unchanged after transcription inhibition. Interestingly, Pol II depletion for a longer time significantly affects the chromatin accessibility and Cohesin occupancy, suggesting that RNA polymerases are capable of affecting the 3D genome indirectly. These direct and indirect effects explain the previous inconsistent findings on the influence of transcription inhibition on the 3D genome.

CONCLUSIONS

We conclude that Pol I, Pol II, and Pol III loss alters local, small-scale chromatin interactions in mammalian cells, suggesting that the 3D chromatin structures are pre-established and relatively stable.

Transcriptional activation during cell reprogramming correlates with the formation of 3D open chromatin hubs

Chromosome structure is a crucial regulatory factor for a wide range of nuclear processes. Chromosome conformation capture (3C)-based experiments combined with computational modelling are pivotal for unveiling 3D chromosome structure. Here, we introduce TADdyn, a tool that integrates time-course 3C data, restraint-based modelling, and molecular dynamics to simulate the structural rearrangements of genomic loci in a completely data-driven way. We apply TADdyn on in situ Hi-C time-course experiments studying the reprogramming of murine B cells to pluripotent cells, and characterize the structural rearrangements that take place upon changes in the transcriptional state of 21 genomic loci of diverse expression dynamics. By measuring various structural and dynamical properties, we find that during gene activation, the transcription starting site contacts with open and active regions in 3D chromatin domains. We propose that these 3D hubs of open and active chromatin may constitute a general feature to trigger and maintain gene transcription.

Heat Shock Factor 1 Drives Intergenic Association of Its Target Gene Loci upon Heat Shock

Transcriptional induction of heat shock protein (HSP) genes is accompanied by dynamic changes in their 3D structure and spatial organization, yet the molecular basis for these phenomena remains unknown. Using chromosome conformation capture and single-cell imaging, we show that genes transcriptionally activated by Hsf1 specifically interact across chromosomes and coalesce into diffraction-limited intranuclear foci. Genes activated by the alternative stress regulators Msn2/Msn4, in contrast, do not interact among themselves nor with Hsf1 targets. Likewise, constitutively expressed genes, even those interposed between HSP genes, show no detectable interaction. Hsf1 forms discrete subnuclear puncta when stress activated, and these puncta dissolve in concert with transcriptional attenuation, paralleling the kinetics of HSP gene coalescence and dissolution. Nuclear Hsf1 and RNA Pol II are both necessary for intergenic HSP gene interactions, while DNA-bound Hsf1 is necessary and sufficient to drive heterologous gene coalescence. Our findings demonstrate that Hsf1 can dynamically restructure the yeast genome.

While gene repositioning is thought to be a general feature of transcription, Chowdhary et al. provide evidence that argues against this concept. The authors demonstrate that Hsf1-regulated genes in Saccharomyces cerevisiae distinctively coalesce into intranuclear foci upon their transcriptional activation, while those activated by alternative transcription factors do not.

Genome-wide analysis of aberrant methylation of enhancer DNA in human osteoarthritis

Background

Osteoarthritis is a chronic musculoskeletal disease characterized by age-related gradual thinning and a high risk in females. Recent studies have shown that DNA methylation plays important roles in osteoarthritis. However, the genome-wide pattern of methylation in enhancers in osteoarthritis remains unclear.

Methods

To explore the function of enhancers in osteoarthritis, we quantified CpG methylation in human enhancers based on a public dataset that included methylation profiles of 470,870 CpG probes in 108 samples from patients with hip and knee osteoarthritis and hip tissues from healthy individuals. Combining various bioinformatics analysis tools, we systematically analyzed aberrant DNA methylation of the enhancers throughout the genome in knee osteoarthritis and hip osteoarthritis.

Results

We identified 16,816 differentially methylated CpGs, and nearly half (8111) of them were from enhancers, suggesting major DNA methylation changes in both types of osteoarthritis in the enhancer regions. A detailed analysis of hip osteoarthritis identified 2426 differentially methylated CpGs in enhancers between male and female patients, and 84.5% of them were hypomethylated in female patients and enriched in phenotypes related to hip osteoarthritis in females. Next, we explored the enhancer methylation dynamics among patients with knee osteoarthritis and identified 280 differentially methylated enhancer CpGs that were enriched in the human phenotypes and disease ontologies related to osteoarthritis. Finally, a comparison of enhancer methylation between knee osteoarthritis and hip osteoarthritis revealed organ source-dependent differences in enhancer methylation.

Conclusion

Our findings indicate that aberrant methylation of enhancers is related to osteoarthritis phenotypes, and a comprehensive atlas of enhancer methylation is useful for further analysis of the epigenetic regulation of osteoarthritis and the development of clinical drugs for treatment of osteoarthritis.

A comprehensive resource for retrieving, visualizing, and integrating functional genomics data

This study presents qcGenomics, a user-friendly online resource for ultrafast retrieval, visualization, and comparative analysis of tens of thousands of genomics datasets to gain new functional insight from global or focused multi-dimensional data integration.

The enormous amount of freely accessible functional genomics data is an invaluable resource for interrogating the biological function of multiple DNA-interacting players and chromatin modifications by large-scale comparative analyses. However, in practice, interrogating large collections of public data requires major efforts for (i) reprocessing available raw reads, (ii) incorporating quality assessments to exclude artefactual and low-quality data, and (iii) processing data by using high-performance computation. Here, we present qcGenomics, a user-friendly online resource for ultrafast retrieval, visualization, and comparative analysis of tens of thousands of genomics datasets to gain new functional insight from global or focused multidimensional data integration.

Novel mRNAs 3' end-associated cis-regulatory elements with epigenomic signatures of mammalian enhancers in the Arabidopsis genome

The precise spatial and temporal control of gene expression requires the coordinated action of genomic cis-regulatory elements (CREs), including transcriptional enhancers. However, our knowledge of enhancers in plants remains rudimentary and only a few plant enhancers have been experimentally defined. Here, we screened the Arabidopsis thaliana genome and identified >1900 unique candidate CREs that carry the genomic signatures of mammalian enhancers. These were termed putative enhancer-like elements (PEs). Nearly all PEs are intragenic and, unexpectedly, most associate with the 3' ends of protein-coding genes. PEs are hotspots for transcription factor binding and harbor motifs resembling cleavage/polyadenylation signals, potentially coupling 3' end processing to the transcriptional regulation of other genes. Hi-C data showed that 24% of PEs are located at regions that can interact intrachromosomally with other protein-coding genes and, surprisingly, many of these target genes interact with PEs through their 3' UTRs. Examination of the genomes of 1135 sequenced Arabidopsis accessions showed that PEs are conserved. Our findings suggest that the identified PEs may serve as transcriptional enhancers and sites for mRNA 3' end processing, and constitute a novel group of CREs in Arabidopsis.

Comprehensive in vivo identification of the c-Myc mRNA protein interactome using HyPR-MS

Proteins bind mRNA through their entire life cycle from transcription to degradation. We analyzed c-Myc mRNA protein interactors in vivo using the HyPR-MS method to capture the crosslinked mRNA by hybridization and then analyzed the bound proteins using mass spectrometry proteomics. Using HyPR-MS, 229 c-Myc mRNA-binding proteins were identified, confirming previously proposed interactors, suggesting new interactors, and providing information related to the roles and pathways known to involve c-Myc. We performed structural and functional analysis of these proteins and validated our findings with a combination of RIP-qPCR experiments, in vitro results released in past studies, publicly available RIP- and eCLIP-seq data, and results from software tools for predicting RNA-protein interactions.

Transcription-driven genome organization: a model for chromosome structure and the regulation of gene expression tested through simulations

Current models for the folding of the human genome see a hierarchy stretching down from chromosome territories, through A/B compartments and topologically-associating domains (TADs), to contact domains stabilized by cohesin and CTCF. However, molecular mechanisms underlying this folding, and the way folding affects transcriptional activity, remain obscure. Here we review physical principles driving proteins bound to long polymers into clusters surrounded by loops, and present a parsimonious yet comprehensive model for the way the organization determines function. We argue that clusters of active RNA polymerases and their transcription factors are major architectural features; then, contact domains, TADs and compartments just reflect one or more loops and clusters. We suggest tethering a gene close to a cluster containing appropriate factors—a transcription factory—increases the firing frequency, and offer solutions to many current puzzles concerning the actions of enhancers, super-enhancers, boundaries and eQTLs (expression quantitative trait loci). As a result, the activity of any gene is directly influenced by the activity of other transcription units around it in 3D space, and this is supported by Brownian-dynamics simulations of transcription factors binding to cognate sites on long polymers.

Enhancer accessibility and CTCF occupancy underlie asymmetric TAD architecture and cell type specific genome topology

Cohesin and CTCF are master regulators of genome topology. How these ubiquitous proteins contribute to cell-type specific genome structure is poorly understood. Here, we explore quantitative aspects of topologically associated domains (TAD) between pluripotent embryonic stem cells (ESC) and lineage-committed cells. ESCs exhibit permissive topological configurations which manifest themselves as increased inter- TAD interactions, weaker intra-TAD interactions, and a unique intra-TAD connectivity whereby one border makes pervasive interactions throughout the domain. Such 'stripe' domains are associated with both poised and active chromatin landscapes and transcription is not a key determinant of their structure. By tracking the developmental dynamics of stripe domains, we show that stripe formation is linked to the functional state of the cell through cohesin loading at lineage-specific enhancers and developmental control of CTCF binding site occupancy. We propose that the unique topological configuration of stripe domains represents a permissive landscape facilitating both productive and opportunistic gene regulation and is important for cellular identity.

Chromosome conformation capture that detects novel cis- and trans-interactions in budding yeast

Chromosome Conformation Capture (3C) has emerged as a powerful approach for revealing the conformation and features of three-dimensional (3D) genomic organization. Yet attainment of higher resolution in organisms with compact genomes presents a challenge. Here, we describe modifications in the 3C technique that substantially enhance its resolution and sensitivity when applied to the 3D genome of budding yeast. Keys to our approach include use of a 4 bp cutter, Taq I, for cleaving the genome and quantitative PCR for measuring the frequency of ligation. Most importantly, we normalize the percent digestion at each restriction site to account for variation in accessibility of local chromatin structure under a given physiological condition. This strategy has led to the detection of physical interactions between regulatory elements and gene coding regions as well as intricate, stimulus-specific interchromosomal interactions between activated genes. We provide an algorithm that incorporates these and other modifications and allows quantitative determination of chromatin interaction frequencies in yeast under any physiological condition.

ChiCMaxima: a robust and simple pipeline for detection and visualization of chromatin looping in Capture Hi-C

Capture Hi-C (CHi-C) is a new technique for assessing genome organization based on chromosome conformation capture coupled to oligonucleotide capture of regions of interest, such as gene promoters. Chromatin loop detection is challenging because existing Hi-C/4C-like tools, which make different assumptions about the technical biases presented, are often unsuitable. We describe a new approach, ChiCMaxima, which uses local maxima combined with limited filtering to detect DNA looping interactions, integrating information from biological replicates. ChiCMaxima shows more stringency and robustness compared to previously developed tools. The tool includes a GUI browser for flexible visualization of CHi-C profiles alongside epigenomic tracks.

SRC-3 Coactivator Governs Dynamic Estrogen-Induced Chromatin Looping Interactions during Transcription

Enhancers are thought to activate transcription by physically contacting promoters via looping. However, direct assays demonstrating these contacts are required to mechanistically verify such cellular determinants of enhancer function. Here, we present versatile cell-free assays to further determine the role of enhancer-promoter contacts (EPCs). We demonstrate that EPC is linked to mutually stimulatory transcription at the enhancer and promoter in vitro. SRC-3 was identified as a critical looping determinant for the estradiol-(E2)-regulated GREB1 locus. Surprisingly, the GREB1 enhancer and promoter contact two internal gene body SRC-3 binding sites, GBS1 and GBS2, which stimulate their transcription. Utilizing time-course 3C assays, we uncovered SRC-3-dependent dynamic chromatin interactions involving the enhancer, promoter, GBS1, and GBS2. Collectively, these data suggest that the enhancer and promoter remain "poised" for transcription via their contacts with GBS1 and GBS2. Upon E2 induction, GBS1 and GBS2 disengage from the enhancer, allowing direct EPC for active transcription.

Regulatory functions of the Mediator kinases CDK8 and CDK19

The Mediator-associated kinases CDK8 and CDK19 function in the context of three additional proteins: CCNC and MED12, which activate CDK8/CDK19 kinase function, and MED13, which enables their association with the Mediator complex. The Mediator kinases affect RNA polymerase II (pol II) transcription indirectly, through phosphorylation of transcription factors and by controlling Mediator structure and function. In this review, we discuss cellular roles of the Mediator kinases and mechanisms that enable their biological functions. We focus on sequence-specific, DNA-binding transcription factors and other Mediator kinase substrates, and how CDK8 or CDK19 may enable metabolic and transcriptional reprogramming through enhancers and chromatin looping. We also summarize Mediator kinase inhibitors and their therapeutic potential. Throughout, we note conserved and divergent functions between yeast and mammalian CDK8, and highlight many aspects of kinase module function that remain enigmatic, ranging from potential roles in pol II promoter-proximal pausing to liquid-liquid phase separation.

Oestrogen promotes tumorigenesis of bladder cancer by inducing the enhancer RNA-eGREB1

In recent years, studies have shown that enhancer RNAs (eRNAs) can be transcribed from enhancers. Increasing evidence has revealed that eRNAs play critical roles in the development of various cancers. Oestrogen‐associated eRNAs are closely related to breast cancer. In view of the gender differences in bladder cancer (BCa), we suppose that oestrogen‐associated eRNAs are also involved in tumorigenesis of BCa. In our study, we first demonstrated that eGREB1 derived from the enhancer of an oestrogen‐responsive gene—GREB1 was up‐regulated in BCa tissues, and the expression level of eGREB1 is positively associated with the histological grade and TNM stage of BCa. Knockdown of eGREB1 by CRISPR‐Cas13a could inhibit cell proliferation, migration and invasion and induce apoptosis in BCa cells T24 and 5637. Besides, we exhibited the promoting effect of oestrogen on BCa cells. What's more, down‐regulation of eGREB1 could improve the malignant biological characteristics of BCa cells induced by oestrogen. In conclusion, our data indicated that eGREB1 plays oncogenic role and oestrogen may promote the occurrence and progression of BCa by inducing eGREB1 production. Our findings provide new insights into the prevention of BCa and develop a novel therapeutic target for the treatment of BCa.

Enhancer RNA - P2RY2e induced by estrogen promotes malignant behaviors of bladder cancer

Enhancers are transcriptional regulatory elements that increase target gene expression. It has reported that enhancers could universally transcribe into enhancer RNAs (eRNAs) with stimulation. Increasing evidence showed eRNAs participated in various disease processes including malignant tumors. P2RY2 enhancer RNA (P2RY2e) is an estrogen-responsive eRNA and involved in the development of breast cancer. However, the relationship between P2RY2e and bladder cancer (BCa) is unclear. In the study, we discovered that P2RY2e was upregulated in BCa tissues and estrogen-treated cells. Estrogen promoted the malignant abilities of BCa cells. P2RY2e knockdown by CRISPR-Cas13a inhibit the cell multiplication, invasion and migration. Additionally, the cell apoptosis was facilitated. What's more, downregulation of P2RY2e could weaken the cancer-promoting effects of estrogen on BCa. Our study revealed that P2RY2e played a carcinogenic role in BCa and estrogen might promote the initiation of BCa by inducing P2RY2e. We provide a potential therapeutic target for BCa and a new perspective for the tumorigenesis of bladder cancer.

Enhancer RNAs (eRNAs): New Insights into Gene Transcription and Disease Treatment

Enhancers are cis-acting elements that have the ability to increase the expression of target genes. Recent studies have shown that enhancers can act as transcriptional units for the production of enhancer RNAs (eRNAs), which are hallmarks of activity enhancers and are involved in the regulation of gene transcription. The in-depth study of eRNAs is of great significance for us to better understand enhancer function and transcriptional regulation in various diseases. Therefore, eRNAs may be a potential therapeutic target for diseases. Here, we review the current knowledge of the characteristics of eRNAs, the molecular mechanisms of eRNAs action, as well as diseases related to dysregulation of eRNAs.

Forces driving the three-dimensional folding of eukaryotic genomes

The last decade has radically renewed our understanding of higher order chromatin folding in the eukaryotic nucleus. As a result, most current models are in support of a mostly hierarchical and relatively stable folding of chromosomes dividing chromosomal territories into A‐ (active) and B‐ (inactive) compartments, which are then further partitioned into topologically associating domains (TADs), each of which is made up from multiple loops stabilized mainly by the CTCF and cohesin chromatin‐binding complexes. Nonetheless, the structure‐to‐function relationship of eukaryotic genomes is still not well understood. Here, we focus on recent work highlighting the biophysical and regulatory forces that contribute to the spatial organization of genomes, and we propose that the various conformations that chromatin assumes are not so much the result of a linear hierarchy, but rather of both converging and conflicting dynamic forces that act on it.

Imaging the dynamics of transcription loops in living chromosomes

When in the lampbrush configuration, chromosomes display thousands of visible DNA loops that are transcribed at exceptionally high rates by RNA polymerase II (pol II). These transcription loops provide unique opportunities to investigate not only the detailed architecture of pol II transcription sites but also the structural dynamics of chromosome looping, which is receiving fresh attention as the organizational principle underpinning the higher-order structure of all chromosome states. The approach described here allows for extended imaging of individual transcription loops and transcription units under conditions in which loop RNA synthesis continues. In intact nuclei from lampbrush-stage Xenopus oocytes isolated under mineral oil, highly specific targeting of fluorescent fusions of the RNA-binding protein CELF1 to nascent transcripts allowed functional transcription loops to be observed and their longevity assessed over time. Some individual loops remained extended and essentially static structures over time courses of up to an hour. However, others were less stable and shrank markedly over periods of 30-60 min in a manner that suggested that loop extension requires continued dense coverage with nascent transcripts. In stable loops and loop-derived structures, the molecular dynamics of the visible nascent RNP component were addressed using photokinetic approaches. The results suggested that CELF1 exchanges freely between the accumulated nascent RNP and the surrounding nucleoplasm, and that it exits RNP with similar kinetics to its entrance. Overall, it appears that on transcription loops, nascent transcripts contribute to a dynamic self-organizing structure that exemplifies a phase-separated nuclear compartment.

Regulation of the cohesin-loading factor NIPBL: Role of the lncRNA NIPBL-AS1 and identification of a distal enhancer element

Cohesin is crucial for genome stability, cell division, transcription and chromatin organization. Its functions critically depend on NIPBL, the cohesin-loader protein that is found to be mutated in >60% of the cases of Cornelia de Lange syndrome (CdLS). Other mutations are described in the cohesin subunits SMC1A, RAD21, SMC3 and the HDAC8 protein. In 25-30% of CdLS cases no mutation in the known CdLS genes is detected. Until now, functional elements in the noncoding genome were not characterized in the molecular etiology of CdLS and therefore are excluded from mutation screening, although the impact of such mutations has now been recognized for a wide range of diseases. We have identified different elements of the noncoding genome involved in regulation of the NIPBL gene. NIPBL-AS1 is a long non-coding RNA transcribed upstream and antisense to NIPBL. By knockdown and transcription blocking experiments, we could show that not the NIPBL-AS1 gene product, but its actual transcription is important to regulate NIPBL expression levels. This reveals a possibility to boost the transcriptional activity of the NIPBL gene by interfering with the NIPBL-AS1 lncRNA. Further, we have identified a novel distal enhancer regulating both NIPBL and NIPBL-AS1. Deletion of the enhancer using CRISPR genome editing in HEK293T cells reduces expression of NIPBL, NIPBL-AS1 as well as genes found to be dysregulated in CdLS.

Heat Shock Protein Genes Undergo Dynamic Alteration in Their Three-Dimensional Structure and Genome Organization in Response to Thermal Stress

Three-dimensional (3D) chromatin organization is important for proper gene regulation, yet how the genome is remodeled in response to stress is largely unknown. Here, we use a highly sensitive version of chromosome conformation capture in combination with fluorescence microscopy to investigate Heat Shock Protein (HSP) gene conformation and 3D nuclear organization in budding yeast. In response to acute thermal stress, HSP genes undergo intense intragenic folding interactions that go well beyond 5'-3' gene looping previously described for RNA polymerase II genes. These interactions include looping between upstream activation sequence (UAS) and promoter elements, promoter and terminator regions, and regulatory and coding regions (gene "crumpling"). They are also dynamic, being prominent within 60 s, peaking within 2.5 min, and attenuating within 30 min, and correlate with HSP gene transcriptional activity. With similarly striking kinetics, activated HSP genes, both chromosomally linked and unlinked, coalesce into discrete intranuclear foci. Constitutively transcribed genes also loop and crumple yet fail to coalesce. Notably, a missense mutation in transcription factor TFIIB suppresses gene looping, yet neither crumpling nor HSP gene coalescence is affected. An inactivating promoter mutation, in contrast, obviates all three. Our results provide evidence for widespread, transcription-associated gene crumpling and demonstrate the de novo assembly and disassembly of HSP gene foci.

Relationship between DNA Methylation in the 5' CpG Island of the SLC47A1 (Multidrug and Toxin Extrusion Protein MATE1) Gene and Interindividual Variability in MATE1 Expression in the Human Liver

Multidrug and toxin extrusion protein 1 (MATE1), which is encoded by solute carrier 47A1 (SLC47A1), mediates the excretion of organic cations into bile and urine. Some genetic variants in human MATE1 altered its transport function in in vitro experiments; however, differences in the pharmacokinetics of substrate drugs cannot be explained by genetic variations in humans. In this study, we investigated whether DNA methylation was involved in interindividual variability in MATE1 expression in the human liver. Approximately 20-fold interindividual variability in MATE1 mRNA expression levels was observed in liver samples and mRNA expression levels negatively correlated with methylation levels of the CpG island in the 27 kb upstream of SLC47A1 DNA demethylation by treatment with 5-aza-2'-deoxycytidine increased MATE1 mRNA expression in MATE1-negative cell lines. The luciferase reporter assay showed that the CpG island increased the transcriptional activity of the SLC47A1 promoter. MATE1 mRNA expression levels were significantly lower in CpG island knockout HepG2 cells than in control cells. These results suggest that the 5' CpG island of SLC47A1 acts as an enhancer for SLC47A1, and DNA methylation in the CpG island plays an important role in interindividual differences in hepatic MATE1 expression.

Histone Hypervariants H2A.Z.1 and H2A.Z.2 Play Independent and Context-Specific Roles in Neuronal Activity-Induced Transcription of Arc/Arg3.1 and Other Immediate Early Genes

The histone variant H2A.Z is an essential and conserved regulator of eukaryotic gene transcription. However, the exact role of this histone in the transcriptional process remains perplexing. In vertebrates, H2A.Z has two hypervariants, H2A.Z.1 and H2A.Z.2, that have almost identical sequences except for three amino acid residues. Due to such similarity, functional specificity of these hypervariants in neurobiological processes, if any, remain largely unknown. In this study with dissociated rat cortical neurons, we asked if H2A.Z hypervariants have distinct functions in regulating basal and activity-induced gene transcription. Hypervariant-specific RNAi and microarray analyses revealed that H2A.Z.1 and H2A.Z.2 regulate basal expression of largely nonoverlapping gene sets, including genes that code for several synaptic proteins. In response to neuronal activity, rapid transcription of our model gene Arc is impaired by depletion of H2A.Z.2, but not H2A.Z.1. This impairment is partially rescued by codepletion of the H2A.Z chaperone, ANP32E. In contrast, under a different context (after 48 h of tetrodotoxin, TTX), rapid transcription of Arc is impaired by depletion of either hypervariant. Such context-dependent roles of H2A.Z hypervariants, as revealed by our multiplexed gene expression assays, are also evident with several other immediate early genes, where regulatory roles of these hypervariants vary from gene to gene under different conditions. Together, our data suggest that H2A.Z hypervariants have context-specific roles that complement each other to mediate activity-induced neuronal gene transcription.