Intragenic Enhancers Attenuate Host Gene Expression

Super-enhancer omics in stem cell

The hallmarks of stem cells, such as proliferation, self-renewal, development, differentiation, and regeneration, are critical to maintain stem cell identity which is sustained by genetic and epigenetic factors. Super-enhancers (SEs), which consist of clusters of active enhancers, play a central role in maintaining stemness hallmarks by specifically transcriptional model. The SE-navigated transcriptional complex, including SEs, non-coding RNAs, master transcriptional factors, Mediators and other co-activators, forms phase-separated condensates, which offers a toggle for directing diverse stem cell fate. With the burgeoning technologies of multiple-omics applied to examine different aspects of SE, we firstly raise the concept of "super-enhancer omics", inextricably linking to Pan-omics. In the review, we discuss the spatiotemporal organization and concepts of SEs, and describe links between SE-navigated transcriptional complex and stem cell features, such as stem cell identity, self-renewal, pluripotency, differentiation and development. We also elucidate the mechanism of stemness and oncogenic SEs modulating cancer stem cells via genomic and epigenetic alterations hijack in cancer stem cell. Additionally, we discuss the potential of targeting components of the SE complex using small molecule compounds, genome editing, and antisense oligonucleotides to treat SE-associated organ dysfunction and diseases, including cancer. This review also provides insights into the future of stem cell research through the paradigm of SEs.

Loss of feedback regulation between FAM3B and androgen receptor driving prostate cancer progression

BACKGROUND

Although the fusion of the transmembrane serine protease 2 gene (TMPRSS2) with the erythroblast transformation-specific-related gene (ERG), or TMPRSS2-ERG, occurs frequently in prostate cancer, its impact on clinical outcomes remains controversial. Roughly half of TMPRSS2-ERG fusions occur through intrachromosomal deletion of interstitial genes and the remainder via insertional chromosomal rearrangements. Because prostate cancers with deletion-derived TMPRSS2-ERG fusions are more aggressive than those with insertional fusions, we investigated the impact of interstitial gene loss on prostate cancer progression.

METHODS

We conducted an unbiased analysis of transcriptome data from large collections of prostate cancer samples and employed diverse in vitro and in vivo models combined with genetic approaches to characterize the interstitial gene loss that imposes the most important impact on clinical outcome.

RESULTS

This analysis identified FAM3B as the top-ranked interstitial gene whose loss is associated with a poor prognosis. The association between FAM3B loss and poor clinical outcome extended to fusion-negative prostate cancers where FAM3B downregulation occurred through epigenetic imprinting. Importantly, FAM3B loss drives disease progression in prostate cancer. FAM3B acts as an intermediator of a self-governing androgen receptor feedback loop. Specifically, androgen receptor upregulates FAM3B expression by binding to an intronic enhancer to induce an enhancer RNA and facilitate enhancer-promoter looping. FAM3B, in turn, attenuates androgen receptor signaling.

CONCLUSION

Loss of FAM3B in prostate cancer, whether through the TMPRSS2-ERG translocation or epigenetic imprinting, causes an exit from this autoregulatory loop to unleash androgen receptor activity and prostate cancer progression. These findings establish FAM3B loss as a new driver of prostate cancer progression and support the utility of FAM3B loss as a biomarker to better define aggressive prostate cancer.

Mapping catalytically engaged TOP2B in neurons reveals the principles of topoisomerase action within the genome

We trapped catalytically engaged topoisomerase IIβ (TOP2B) in covalent DNA cleavage complexes (TOP2Bccs) and mapped their positions genome-wide in cultured mouse cortical neurons. We report that TOP2Bcc distribution varies with both nucleosome and compartmental chromosome organization. While TOP2Bccs in gene bodies correlate with their level of transcription, highly expressed genes that lack the usually associated chromatin marks, such as H3K36me3, show reduced TOP2Bccs, suggesting that histone posttranslational modifications regulate TOP2B activity. Promoters with high RNA polymerase II occupancy show elevated TOP2B chromatin immunoprecipitation sequencing signals but low TOP2Bccs, indicating that TOP2B catalytic engagement is curtailed at active promoters. Surprisingly, either poisoning or inhibiting TOP2B increases nascent transcription at most genes and enhancers but reduces transcription within long genes. These effects are independent of transcript length and instead correlate with the presence of intragenic enhancers. Together, these results clarify how cells modulate the catalytic engagement of topoisomerases to affect transcription.

TALE transcription factors: Cofactors no more

Exd/PBX, Hth/MEIS and PREP proteins belong to the TALE (three-amino-acid loop extension) superclass of transcription factors (TFs) with an atypical homedomain (HD). Originally discovered as "cofactors" to HOX proteins, revisiting their traditional role in light of genome-wide experiments reveals a strong and reproducible pattern of HOX and TALE co-occupancy across diverse embryonic tissues. While confirming that TALE increases HOX specificity and selectivity in vivo, this wider outlook also reveals novel aspects of HOX:TALE collaboration, namely that HOX TFs generally require pre-bound TALE factors to access their functional binding sites in vivo. In contrast to the restricted expression domains of HOX TFs, TALE factors are largely ubiquitous, and PBX and PREP are expressed at the earliest developmental stages. PBX and MEIS control development of many organs and tissues and their dysregulation is associated with congenital disease and cancer. Accordingly, many instances of TALE cooperation with non HOX TFs have been documented in various systems. The model that emerges from these studies is that TALE TFs create a permissive chromatin platform that is selected by tissue-restricted TFs for binding. In turn, HOX and other tissue-restricted TFs selectively convert a ubiquitous pool of low affinity TALE binding events into high confidence, tissue-restricted binding events associated with transcriptional activation. As a result, TALE:TF complexes are associated with active chromatin and domain/lineage-specific gene activity. TALE ubiquitous expression and broad genomic occupancy, as well as the increasing examples of TALE tissue-specific partners, reveal a universal and obligatory role for TALE in the control of tissue and lineage-specific transcriptional programs, beyond their initial discovery as HOX co-factors.

The chromatin signatures of enhancers and their dynamic regulation

Enhancers are cis-regulatory elements that can stimulate gene expression from distance, and drive precise spatiotemporal gene expression profiles during development. Functional enhancers display specific features including an open chromatin conformation, Histone H3 lysine 27 acetylation, Histone H3 lysine 4 mono-methylation enrichment, and enhancer RNAs production. These features are modified upon developmental cues which impacts their activity. In this review, we describe the current state of knowledge about enhancer functions and the diverse chromatin signatures found on enhancers. We also discuss the dynamic changes of enhancer chromatin signatures, and their impact on lineage specific gene expression profiles, during development or cellular differentiation.

eRNA co-expression network uncovers TF dependency and convergent cooperativity

Enhancer RNAs (eRNAs) are non-coding RNAs produced by transcriptional enhancers that are highly correlated with their activity. Using a capped nascent RNA sequencing (PRO-cap) dataset in human lymphoblastoid cell lines across 67 individuals, we identified inter-individual variation in the expression of over 80 thousand transcribed transcriptional regulatory elements (tTREs), in both enhancers and promoters. Co-expression analysis of eRNAs from tTREs across individuals revealed how enhancers are associated with each other and with promoters. Mid- to long-range co-expression showed a distance-dependent decay that was modified by TF occupancy. In particular, we found a class of "bivalent" TFs, including Cohesin, that both facilitate and isolate the interaction between enhancers and/or promoters, depending on their topology. At short distances, we observed strand-specific correlations between nearby eRNAs in both convergent and divergent orientations. Our results support a cooperative model of convergent eRNAs, consistent with eRNAs facilitating adjacent enhancers rather than interfering with each other. Therefore, our approach to infer functional interactions from co-expression analyses provided novel insights into the principles of enhancer interactions as a function of distance, orientation, and binding landscapes of TFs.

Hypoxic regulation of hypoxia inducible factor 1 alpha via antisense transcription

Impaired oxygen homeostasis is a frequently encountered pathophysiological factor in multiple complex diseases, including cardiovascular disease and cancer. While the canonical hypoxia response pathway is well characterized, less is known about the role of noncoding RNAs in this process. Here, we investigated the nascent and steady-state noncoding transcriptional responses in endothelial cells and their potential roles in regulating the hypoxic response. Notably, we identify a novel antisense long noncoding RNA that convergently overlaps the majority of the hypoxia inducible factor 1 alpha (HIF1A) locus, which is expressed across several cell types and elevated in atherosclerotic lesions. The antisense (HIF1A-AS) is produced as a stable, unspliced, and polyadenylated nuclear retained transcript. HIF1A-AS is highly induced in hypoxia by both HIF1A and HIF2A and exhibits anticorrelation with the coding HIF1A transcript and protein expression. We further characterized this functional relationship by CRISPR-mediated bimodal perturbation of the HIF1A-AS promoter. We provide evidence that HIF1A-AS represses the expression of HIF1a in cis by repressing transcriptional elongation and deposition of H3K4me3, and that this mechanism is dependent on the act of antisense transcription itself. Overall, our results indicate a critical regulatory role of antisense mediated transcription in regulation of HIF1A expression and cellular response to hypoxia.

Insight into the regulatory mechanism of dynamic chromatin 3D interactions during cardiomyocyte differentiation in human

Cardiogenesis is an extremely complicated process involved with DNA regulatory elements, and trans factors regulate gene expression pattern spatiotemporally. Enhancers, as the well-known DNA elements, activate target gene expression by transcription factors (TFs) occupied to organize dynamic three-dimensional (3D) interactions, which when affected or interrupted might cause heart defects or diseases. In this study, we integrated transcriptome, 3D genome, and regulatome to reorganize the global 3D genome in cardiomyogenesis, showing a gradually decreased trend of both chromatin interactions and topological associating domains (TADs) during cardiomyocyte differentiation. And almost all of the chromatin interactions occurred within the same or between adjacent TADs involved with enhancers, indicating that dynamical rewiring of enhancer-related chromatin interactions in the continuous expansive TADs is closely correlated to cardiogenesis. Moreover, we found stage-specific interactions activate stage-specific expression to be involved within corresponding biological functions, and the stage-specific combined regulations of enhancers and binding TFs form connected networks to control stage-specific expression and biological processes, which promote cardiomyocyte differentiation. Finally, we identified markers based on regulatory networks, which might drive cardiac development. This study demonstrates the power of enhancer interactome combined with active TFs to reveal insights into transcriptional regulatory networks during cardiomyogenesis.

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.

PAD2: interactive exploration of transcription factor genomic colocalization using ChIP-seq data

Characterizing transcription factor (TF) genomic colocalization is essential for identifying cooperative binding of TFs in controlling gene expression. Here, we introduce a protocol for using PAD2, an interactive web application that enables the investigation of colocalization of various TFs and chromatin-regulating proteins from mouse embryonic stem cells at various functional genomic regions. We describe steps for accessing and searching the PAD2 database and selecting and submitting genomic regions. We then detail protein colocalization analysis using heatmap and ranked correlation plot. For complete details on the use and execution of this protocol, please refer to Kim et al. (2022).¹.

Alternative promoters in CpG depleted regions are prevalently associated with epigenetic misregulation of liver cancer transcriptomes

Transcriptional regulation is commonly governed by alternative promoters. However, the regulatory architecture in alternative and reference promoters, and how they differ, remains elusive. In 100 CAGE-seq libraries from hepatocellular carcinoma patients, here we annotate 4083 alternative promoters in 2926 multi-promoter genes, which are largely undetected in normal livers. These genes are enriched in oncogenic processes and predominantly show association with overall survival. Alternative promoters are narrow nucleosome depleted regions, CpG island depleted, and enriched for tissue-specific transcription factors. Globally tumors lose DNA methylation. We show hierarchical retention of intragenic DNA methylation with CG-poor regions rapidly losing methylation, while CG-rich regions retain it, a process mediated by differential SETD2, H3K36me3, DNMT3B, and TET1 binding. This mechanism is validated in SETD2 knockdown cells and SETD2-mutated patients. Selective DNA methylation loss in CG-poor regions makes the chromatin accessible for alternative transcription. We show alternative promoters can control tumor transcriptomes and their regulatory architecture.

Cis-regulatory atlas of primary human CD4+ T cells

Cis-regulatory elements (CRE) are critical for coordinating gene expression programs that dictate cell-specific differentiation and homeostasis. Recently developed self-transcribing active regulatory region sequencing (STARR-Seq) has allowed for genome-wide annotation of functional CREs. Despite this, STARR-Seq assays are only employed in cell lines, in part, due to difficulties in delivering reporter constructs. Herein, we implemented and validated a STARR-Seq-based screen in human CD4+ T cells using a non-integrating lentiviral transduction system. Lenti-STARR-Seq is the first example of a genome-wide assay of CRE function in human primary cells, identifying thousands of functional enhancers and negative regulatory elements (NREs) in human CD4+ T cells. We find an unexpected difference in nucleosome organization between enhancers and NRE: enhancers are located between nucleosomes, whereas NRE are occupied by nucleosomes in their endogenous locations. We also describe chromatin modification, eRNA production, and transcription factor binding at both enhancers and NREs. Our findings support the idea of silencer repurposing as enhancers in alternate cell types. Collectively, these data suggest that Lenti-STARR-Seq is a successful approach for CRE screening in primary human cell types, and provides an atlas of functional CREs in human CD4+ T cells.

DDM1-mediated gene body DNA methylation is associated with inducible activation of defense-related genes in Arabidopsis

BACKGROUND

Plants memorize previous pathogen attacks and are "primed" to produce a faster and stronger defense response, which is critical for defense against pathogens. In plants, cytosines in transposons and gene bodies are reported to be frequently methylated. Demethylation of transposons can affect disease resistance by regulating the transcription of nearby genes during defense response, but the role of gene body methylation (GBM) in defense responses remains unclear.

RESULTS

Here, we find that loss of the chromatin remodeler decrease in DNA methylation 1 (ddm1) synergistically enhances resistance to a biotrophic pathogen under mild chemical priming. DDM1 mediates gene body methylation at a subset of stress-responsive genes with distinct chromatin properties from conventional gene body methylated genes. Decreased gene body methylation in loss of ddm1 mutant is associated with hyperactivation of these gene body methylated genes. Knockout of glyoxysomal protein kinase 1 (gpk1), a hypomethylated gene in ddm1 loss-of-function mutant, impairs priming of defense response to pathogen infection in Arabidopsis. We also find that DDM1-mediated gene body methylation is prone to epigenetic variation among natural Arabidopsis populations, and GPK1 expression is hyperactivated in natural variants with demethylated GPK1.

CONCLUSIONS

Based on our collective results, we propose that DDM1-mediated GBM provides a possible regulatory axis for plants to modulate the inducibility of the immune response.

Revealing eRNA interactions: TF dependency and convergent cooperativity

Enhancer RNAs (eRNAs) are non-coding RNAs produced from transcriptional enhancers that are highly correlated with their activities. Using capped nascent RNA sequencing (PRO-cap) dataset in human lymphoblastoid cell lines across individuals, we identified inter-individual variation of expression in over 80 thousand transcribed transcriptional regulatory elements (tTREs), in both enhancers and promoters. Co-expression analysis of eRNAs from tTREs across individuals revealed how enhancers interact with each other and with promoters. Mid-to-long range interactions showed distance-dependent decay, which was modified by TF occupancy. In particular, we found a class of 'bivalent' TFs, including Cohesin, which both facilitates and insulates the interaction between enhancers and/or promoters depending on the topology. In short ranges, we observed strand specific interactions between nearby eRNAs in both convergent or divergent orientations. Our finding supports a cooperative convergent eRNA model, which is compatible with eRNA remodeling neighboring enhancers rather than interfering with each other. Therefore, our approach to infer functional interactions from co-expression analyses provided novel insights into the principles of enhancer interactions depending on the distance, orientation, and the binding landscapes of TFs.

Dynamic interplay between non-coding enhancer transcription and gene activity in development

Non-coding transcription at the intergenic regulatory regions is a prevalent feature of metazoan genomes, but its biological function remains uncertain. Here, we devise a live-imaging system that permits simultaneous visualization of gene activity along with intergenic non-coding transcription at single-cell resolution in Drosophila. Quantitative image analysis reveals that elongation of RNA polymerase II across the internal core region of enhancers leads to suppression of transcriptional bursting from linked genes. Super-resolution imaging and genome-editing analysis further demonstrate that enhancer transcription antagonizes molecular crowding of transcription factors, thereby interrupting the formation of a transcription hub at the gene locus. We also show that a certain class of developmental enhancers are structurally optimized to co-activate gene transcription together with non-coding transcription effectively. We suggest that enhancer function is flexibly tunable through the modulation of hub formation via surrounding non-coding transcription during development.

Random allelic expression in the adult human body

Genes are typically assumed to express both parental alleles similarly, yet cell lines show random allelic expression (RAE) for many autosomal genes that could shape genetic effects. Thus, understanding RAE in human tissues could improve our understanding of phenotypic variation. Here, we develop a methodology to perform genome-wide profiling of RAE and biallelic expression in GTEx datasets for 832 people and 54 tissues. We report 2,762 autosomal genes with some RAE properties similar to randomly inactivated X-linked genes. We found that RAE is associated with rapidly evolving regions in the human genome, adaptive signaling processes, and genes linked to age-related diseases such as neurodegeneration and cancer. We define putative mechanistic subtypes of RAE distinguished by gene overlaps on sense and antisense DNA strands, aggregation in clusters near telomeres, and increased regulatory complexity and inputs compared with biallelic genes. We provide foundations to study RAE in human phenotypes, evolution, and disease.

TBX2 acts as a potent transcriptional silencer of tumour suppressor genes through interaction with the CoREST complex to sustain the proliferation of breast cancers

Chromosome 17q23 amplification occurs in 20% of primary breast tumours and is associated with poor outcome. The TBX2 gene is located on 17q23 and is often over-expressed in this breast tumour subset. TBX2 is an anti-senescence gene, promoting cell growth and survival through repression of Tumour Suppressor Genes (TSGs), such as NDRG1 and CST6. Previously we found that TBX2 cooperates with the PRC2 complex to repress several TSGs, and that PRC2 inhibition restored NDRG1 expression to impede cellular proliferation. Here, we now identify CoREST proteins, LSD1 and ZNF217, as novel interactors of TBX2. Genetic or pharmacological targeting of CoREST emulated TBX2 loss, inducing NDRG1 expression and abolishing breast cancer growth in vitro and in vivo. Furthermore, we uncover that TBX2/CoREST targeting of NDRG1 is achieved by recruitment of TBX2 to the NDRG1 promoter by Sp1, the abolishment of which resulted in NDRG1 upregulation and diminished cancer cell proliferation. Through ChIP-seq we reveal that 30% of TBX2-bound promoters are shared with ZNF217 and identify novel targets repressed by TBX2/CoREST; of these targets a lncRNA, LINC00111, behaves as a negative regulator of cell proliferation. Overall, these data indicate that inhibition of CoREST proteins represents a promising therapeutic intervention for TBX2-addicted breast tumours.

Large-scale multi-omics analysis suggests specific roles for intragenic cohesin in transcriptional regulation

Cohesin, an essential protein complex for chromosome segregation, regulates transcription through a variety of mechanisms. It is not a trivial task to assign diverse cohesin functions. Moreover, the context-specific roles of cohesin-mediated interactions, especially on intragenic regions, have not been thoroughly investigated. Here we perform a comprehensive characterization of cohesin binding sites in several human cell types. We integrate epigenomic, transcriptomic and chromatin interaction data to explore the context-specific functions of intragenic cohesin related to gene activation. We identify a specific subset of cohesin binding sites, decreased intragenic cohesin sites (DICs), which are negatively correlated with transcriptional regulation. A subgroup of DICs is enriched with enhancer markers and RNA polymerase II, while the others are more correlated to chromatin architecture. DICs are observed in various cell types, including cells from patients with cohesinopathy. We also implement machine learning to our data and identified genomic features for isolating DICs from all cohesin sites. These results suggest a previously unidentified function of cohesin on intragenic regions for transcriptional regulation.

Enhancer RNA: What we know and what we can achieve

Enhancers are important cis-acting elements that can regulate gene transcription and cell fate alongside promoters. In fact, many human cancers and diseases are associated with the malfunction of enhancers. Recent studies have shown that enhancers can produce enhancer RNAs (eRNAs) by RNA polymerase II. In this review, we discuss eRNA production, characteristics, functions and mechanics. eRNAs can determine chromatin accessibility, histone modification and gene expression by constructing a 'chromatin loop', thereby bringing enhancers to their target gene. eRNA can also be involved in the phase separation with enhancers and other proteins. eRNAs are abundant, and importantly, tissue-specific in tumours, various diseases and stem cells; thus, eRNAs can be a potential target for disease diagnosis and treatment. As eRNA is produced from the active transcription of enhancers and is involved in the regulation of cell fate, its manipulation will influence cell function, and therefore, it can be a new target for biological therapy.

Intragenic CpG Islands and Their Impact on Gene Regulation

The mammalian genome is depleted in CG dinucleotides, except at protected regions where they cluster as CpG islands (CGIs). CGIs are gene regulatory hubs and serve as transcription initiation sites and are as expected, associated with gene promoters. Advances in genomic annotations demonstrate that a quarter of CGIs are found within genes. Such intragenic regions are repressive environments, so it is surprising that CGIs reside here and even more surprising that some resist repression and are transcriptionally active within a gene. Hence, intragenic CGI positioning within genes is not arbitrary and is instead, selected for. As a wealth of recent studies demonstrate, intragenic CGIs are embedded within genes and consequently, influence ‘host’ gene mRNA isoform length and expand transcriptome diversity.

Placental DNA methylation marks are associated with maternal depressive symptoms during early pregnancy

Maternal depressive symptoms during pregnancy are a significant risk factor for adverse developmental and health outcomes of the offspring. The molecular mechanisms mediating the long-term effects of this exposure are not well understood. Previous studies have found association between prenatal exposure to maternal psychological distress and placental DNA methylation of candidate genes, which can influence placental barrier function and development of the fetus. Our objective in this study was to determine epigenome wide association of maternal depressive symptoms in early pregnancy with the placental DNA methylation. For this purpose we examined DNA methylomes of 92 placental samples by using reduced representation bisulfite sequencing. The placental samples were collected after deliveries of 39 girls and 59 boys, whose mothers had Edinburgh Postnatal Depression Score ranging from 0 to 19 at gestational week 14. According to our results maternal depressive symptoms are associated with DNA methylation of 2833 CpG sites, which are particularly over-represented in genic enhancers. The genes overlapping or nearest to these sites are functionally enriched for development of neurons and show expression enrichment in several regions of developing brain. The genomic regions harboring the DNA methylation marks are enriched for single nucleotide polymorphisms associated with mental disease trait class. Potential cellular signaling cascades mediating the effects include inflammatory and hormonal pathways. As a conclusion our results suggest that maternal depressive symptoms during early pregnancy are associated with DNA methylation marks in placenta in genes, which are important for the development and long-term health of the brain. Whether similar marks can be detected in exposed children remains to be elucidated in further studies.

Reactivation of tumour suppressor in breast cancer by enhancer switching through NamiRNA network

Dysfunction of Tumour Suppressor Genes (TSGs) is a common feature in carcinogenesis. Epigenetic abnormalities including DNA hypermethylation or aberrant histone modifications in promoter regions have been described for interpreting TSG inactivation. However, in many instances, how TSGs are silenced in tumours are largely unknown. Given that miRNA with low expression in tumours is another recognized signature, we hypothesize that low expression of miRNA may reduce the activity of TSG related enhancers and further lead to inactivation of TSG during cancer development. Here, we reported that low expression of miRNA in cancer as a recognized signature leads to loss of function of TSGs in breast cancer. In 157 paired breast cancer and adjacent normal samples, tumour suppressor gene GPER1 and miR-339 are both downregulated in Luminal A/B and Triple Negative Breast Cancer subtypes. Mechanistic investigations revealed that miR-339 upregulates GPER1 expression in breast cancer cells by switching on the GPER1 enhancer, which can be blocked by enhancer deletion through the CRISPR/Cas9 system. Collectively, our findings reveal novel mechanistic insights into TSG dysfunction in cancer development, and provide evidence that reactivation of TSG by enhancer switching may be a promising alternative strategy for clinical breast cancer treatment.

Antisense RNAs during early vertebrate development are divided in groups with distinct features

Long noncoding RNAs or lncRNAs are a class of non-protein-coding RNAs that are >200 nt in length. Almost 50% of lncRNAs during zebrafish development are transcribed in an antisense direction to a protein-coding gene. However, the role of these natural antisense transcripts (NATs) during development remains enigmatic. To understand NATs in early vertebrate development, we took a computational biology approach and analyzed existing as well as novel data sets. Our analysis indicates that zebrafish NATs can be divided into two major classes based on their coexpression patterns with respect to the overlapping protein-coding genes. Group 1 NATs have characteristics similar to maternally deposited RNAs in that their levels decrease as development progresses. Group 1 NAT levels are negatively correlated with that of overlapping sense-strand protein-coding genes. Conversely, Group 2 NATs are coexpressed with overlapping protein-coding genes. In contrast to Group 1, which is enriched in genes involved in developmental pathways, Group 2 protein-coding genes are enriched in housekeeping functions. Group 1 NATs also show larger overlap and higher complementarity with the sense-strand mRNAs compared to other NATs. In addition, our transcriptomics data, quantifying RNA levels from cytoplasmic and nuclear compartments, indicates that Group 1 NATs are more abundant in the cytosol. Based on their expression pattern, cytosolic nature, and their higher complementarity to the overlapping developmental mRNAs, we speculate that Group 1 NATs function post-transcriptionally to silence spurious expression of developmental genes.

Implications of Enhancer Transcription and eRNAs in Cancer

Despite extensive progress in developing anticancer therapies, therapy resistance remains a major challenge that promotes disease relapse. The changes that lead to therapy resistance can be intrinsically present or may be initiated during treatment. Genetic and epigenetic heterogeneity in tumors make it more challenging to deal with therapy resistance. Recent advances in genome-wide analyses have revealed that the deregulation of distal gene regulatory elements, such as enhancers, appears in several pathophysiological conditions, including cancer. Beyond the conventional function of enhancers in recruiting transcription factors to gene promoters, enhancer elements are also transcribed into noncoding RNAs known as enhancer RNAs (eRNA). Accumulating evidence suggests that uncontrolled enhancer activity with aberrant eRNA expression promotes oncogenesis. Interestingly, tissue-specific, transcribed eRNAs from active enhancers can serve as potential therapeutic targets or biomarkers in several cancer types. This review provides a comprehensive overview of the mechanisms of enhancer transcription and eRNAs as well as their potential roles in cancer and drug resistance.

Transcriptional Silencers: Driving Gene Expression with the Brakes On

Silencers are regulatory DNA elements that reduce transcription from their target promoters; they are the repressive counterparts of enhancers. Although discovered decades ago, and despite evidence of their importance in development and disease, silencers have been much less studied than enhancers. Recently, however, a series of papers have reported systematic studies of silencers in various model systems. Silencers are often bifunctional regulatory elements that can also act as enhancers, depending on cellular context, and are enriched for expression quantitative trait loci (eQTLs) and disease-associated variants. There is not yet evidence of a 'silencer chromatin signature', in the distribution of histone modifications or associated proteins, that is common to all silencers; instead, silencers may fall into various subclasses, acting by distinct (and possibly overlapping) mechanisms.

Temporal dissection of an enhancer cluster reveals distinct temporal and functional contributions of individual elements

Many genes are regulated by multiple enhancers that often simultaneously activate their target gene. However, how individual enhancers collaborate to activate transcription is not well understood. Here, we dissect the functions and interdependencies of five enhancer elements that together activate Fgf5 expression during exit from naive murine pluripotency. Four intergenic elements form a super-enhancer, and most of the elements contribute to Fgf5 induction at distinct time points. A fifth, poised enhancer located in the first intron contributes to Fgf5 expression at every time point by amplifying overall Fgf5 expression levels. Despite low individual enhancer activity, together these elements strongly induce Fgf5 expression in a super-additive fashion that involves strong accumulation of RNA polymerase II at the intronic enhancer. Finally, we observe a strong anti-correlation between RNA polymerase II levels at enhancers and their distance to the closest promoter, and we identify candidate elements with properties similar to the intronic enhancer.

Enhancer RNAs: transcriptional regulators and workmates of NamiRNAs in myogenesis

miRNAs are well known to be gene repressors. A newly identified class of miRNAs termed nuclear activating miRNAs (NamiRNAs), transcribed from miRNA loci that exhibit enhancer features, promote gene expression via binding to the promoter and enhancer marker regions of the target genes. Meanwhile, activated enhancers produce endogenous non-coding RNAs (named enhancer RNAs, eRNAs) to activate gene expression. During chromatin looping, transcribed eRNAs interact with NamiRNAs through enhancer-promoter interaction to perform similar functions. Here, we review the functional differences and similarities between eRNAs and NamiRNAs in myogenesis and disease. We also propose models demonstrating their mutual mechanism and function. We conclude that eRNAs are active molecules, transcriptional regulators, and partners of NamiRNAs, rather than mere RNAs produced during enhancer activation.

Regulatory Non-Coding RNAs: An Overview

The discovery of new classes of non-coding RNAs has always been preceded or accompanied by technological breakthroughs, and these outstanding progresses in transcriptomics approaches enabled to regularly add new members to the list. From the first detection of tRNAs, through the revolution of miRNAs discovery, to the recent identification of eRNAs or the identification of new functions for already known ncRNAs, this introductive review provides a very concise historical and functional overview of most prominent small regulatory non-coding RNA families.

Causes and consequences of RNA polymerase II stalling during transcript elongation

The journey of RNA polymerase II (Pol II) as it transcribes a gene is anything but a smooth ride. Transcript elongation is discontinuous and can be perturbed by intrinsic regulatory barriers, such as promoter-proximal pausing, nucleosomes, RNA secondary structures and the underlying DNA sequence. More substantial blocking of Pol II translocation can be caused by other physiological circumstances and extrinsic obstacles, including other transcribing polymerases, the replication machinery and several types of DNA damage, such as bulky lesions and DNA double-strand breaks. Although numerous different obstacles cause Pol II stalling or arrest, the cell somehow distinguishes between them and invokes different mechanisms to resolve each roadblock. Resolution of Pol II blocking can be as straightforward as temporary backtracking and transcription elongation factor S-II (TFIIS)-dependent RNA cleavage, or as drastic as premature transcription termination or degradation of polyubiquitylated Pol II and its associated nascent RNA. In this Review, we discuss the current knowledge of how these different Pol II stalling contexts are distinguished by the cell, how they overlap with each other, how they are resolved and how, when unresolved, they can cause genome instability.

Elongin A associates with actively transcribed genes and modulates enhancer RNA levels with limited impact on transcription elongation rate in vivo

Elongin A (EloA) is an essential transcription factor that stimulates the rate of RNA polymerase II (Pol II) transcription elongation in vitro. However, its role as a transcription factor in vivo has remained underexplored. Here we show that in mouse embryonic stem cells, EloA localizes to both thousands of Pol II transcribed genes with preference for transcription start site and promoter regions and a large number of active enhancers across the genome. EloA deletion results in accumulation of transcripts from a subset of enhancers and their adjacent genes. Notably, EloA does not substantially enhance the elongation rate of Pol II in vivo. We also show that EloA localizes to the nucleoli and associates with RNA polymerase I transcribed ribosomal RNA gene, Rn45s. EloA is a highly disordered protein, which we demonstrate forms phase-separated condensates in vitro, and truncation mutations in the intrinsically disordered regions (IDR) of EloA interfere with its targeting and localization to the nucleoli. We conclude that EloA broadly associates with transcribed regions, tunes RNA Pol II transcription levels via impacts on enhancer RNA synthesis, and interacts with the rRNA producing/processing machinery in the nucleolus. Our work opens new avenues for further investigation of the role of this functionally multifaceted transcription factor in enhancer and ribosomal RNA biology.

Population-scale study of eRNA transcription reveals bipartite functional enhancer architecture

Enhancer RNAs (eRNA) are unstable non-coding RNAs, transcribed bidirectionally from active regulatory sequences, whose expression levels correlate with enhancer activity. We use capped-nascent-RNA sequencing to efficiently capture bidirectional transcription initiation across several human lymphoblastoid cell lines (Yoruba population) and detect ~75,000 eRNA transcription sites with high sensitivity and specificity. The use of nascent-RNA sequencing sidesteps the confounding effect of eRNA instability. We identify quantitative trait loci (QTLs) associated with the level and directionality of eRNA expression. High-resolution analyses of these two types of QTLs reveal distinct positions of enrichment at the central transcription factor (TF) binding regions and at the flanking eRNA initiation regions, both of which are associated with mRNA expression QTLs. These two regions-the central TF-binding footprint and the eRNA initiation cores-define a bipartite architecture of enhancers, inform enhancer function, and can be used as an indicator of the significance of non-coding regulatory variants.

An Intronic Alu Element Attenuates the Transcription of a Long Non-coding RNA in Human Cell Lines

Alu elements are primate-specific repeats and represent the most abundant type of transposable elements (TE) in the human genome. Genome-wide analysis of the enrichment of histone post-translational modifications suggests that human Alu sequences could function as transcriptional enhancers; however, no functional experiments have evaluated the role of Alu sequences in the control of transcription in situ. The present study analyses the regulatory activity of a human Alu sequence from the AluSx family located in the second intron of the long intergenic non-coding RNA Linc00441, found in divergent orientation to the RB1 gene. We observed that the Alu sequence acts as an enhancer element based on reporter gene assays while CRISPR-Cas9 deletions of the Alu sequence in K562 cells resulted in a marked transcriptional upregulation of Linc00441 and a decrease in proliferation. Our results suggest that an intragenic Alu sequence with enhancer activity can act as a transcriptional attenuator of its host lincRNA.

Spirits in the Material World: Enhancer RNAs in Transcriptional Regulation

Responses to developmental and environmental cues depend on precise spatiotemporal control of gene transcription. Enhancers, which comprise DNA elements bound by regulatory proteins, can activate target genes in response to these external signals. Recent studies have shown that enhancers are transcribed to produce enhancer RNAs (eRNAs). Do eRNAs play a functional role in activating gene expression or are they non-functional byproducts of nearby transcription machinery? The unstable nature of eRNAs and over-reliance on knockdown approaches have made elucidating the possible functions of eRNAs challenging. We focus here on studies using cloned eRNAs to study their function as transcripts, revealing roles for eRNAs in enhancer-promoter looping, recruiting transcriptional machinery, and facilitating RNA polymerase pause-release to regulate gene expression.

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.

Transcription in Living Cells: Molecular Mechanisms of Bursting

Transcription in several organisms from certain bacteria to humans has been observed to be stochastic in nature: toggling between active and inactive states. Periods of active nascent RNA synthesis known as bursts represent individual gene activation events in which multiple polymerases are initiated. Therefore, bursting is the single locus illustration of both gene activation and repression. Although transcriptional bursting was originally observed decades ago, only recently have technological advances enabled the field to begin elucidating gene regulation at the single-locus level. In this review, we focus on how biochemical, genomic, and single-cell data describe the regulatory steps of transcriptional bursts.

Genomic Rewiring of SOX2 Chromatin Interaction Network during Differentiation of ESCs to Postmitotic Neurons

Cellular differentiation requires dramatic changes in chromatin organization, transcriptional regulation, and protein production. To understand the regulatory connections between these processes, we generated proteomic, transcriptomic, and chromatin accessibility data during differentiation of mouse embryonic stem cells (ESCs) into postmitotic neurons and found extensive associations between different molecular layers within and across differentiation time points. We observed that SOX2, as a regulator of pluripotency and neuronal genes, redistributes from pluripotency enhancers to neuronal promoters during differentiation, likely driven by changes in its protein interaction network. We identified ATRX as a major SOX2 partner in neurons, whose co-localization correlated with an increase in active enhancer marks and increased expression of nearby genes, which we experimentally confirmed for three loci. Collectively, our data provide key insights into the regulatory transformation of SOX2 during neuronal differentiation, and we highlight the significance of multi-omic approaches in understanding gene regulation in complex systems.

•

Complex interplay of RNA, protein, and chromatin during neuronal differentiation

•

Multi-omic profiling reveals divergent roles of SOX2 in stem cells and neurons

•

SOX2 on-chromatin interaction network changes from pluripotent to neuronal factors

•

ATRX interacts with SOX2 in neurons and co-binds highly expressed neuronal genes

Enhancer RNAs are an important regulatory layer of the epigenome

Noncoding RNAs (ncRNAs) direct a remarkable number of diverse functions in development and disease through their regulation of transcription, RNA processing and translation. Leading the charge in the RNA revolution is a class of ncRNAs that are synthesized at active enhancers, called enhancer RNAs (eRNAs). Here, we review recent insights into the biogenesis of eRNAs and the mechanisms underlying their multifaceted functions and consider how these findings could inform future investigations into enhancer transcription and eRNA function.

LINC01198 facilitates gliomagenesis through activating PI3K/AKT pathway

Glioma is the most malignant primary brain cancer which frequently occurred in adults. In recent years, long-non coding RNAs (lncRNAs) have been demonstrated to play pivotal roles in human cancers. However, the role of most lncRNAs in gliomagenesis has not been probed. Presently, through TCGA, a novel lncRNA LINC01198 was found to be up-regulated and associated with clinical outcomes in glioblastoma multiforme (GBM). In our study, LINC01198 was proved to be up-regulated in glioma cell lines, and silenced LINC01198 curbed glioma cell proliferation and accelerated cell apoptosis. Importantly, we corroborated that LINC01198 activated the PI3 K/AKT pathway to aggravate glioma progression by targeting PIK3 CA and PTEN. Subsequently, LINC01198 was validated to localize in both cytoplasm and nucleus of glioma cells. Through mechanistic exploration, we illustrated that LINC01198 increased PIK3CA expression by sponging miR-129-5p in the cytoplasm. Furthermore, PTEN was transcriptionally repressed by REST/RCOR1/HDAC2 complex. More importantly, LINC01198 accelerated the assembly of REST/RCOR1/HDAC2 complex and recruited such complex to PTEN promoter so as to impair PTEN expression in glioma. Finally, we further verified that LINC01198 hindered glioma tumour growth in vivo through AKT-dependent manner. Jointly, LINC01198 activates PI3 K/AKT signalling to exert oncogenic function in gliomagenesis by regulating PIK3CA and PTEN, which highlights a new approach for glioma treatment.

Towards a comprehensive catalogue of validated and target-linked human enhancers

The human gene catalogue is essentially complete, but we lack an equivalently vetted inventory of bona fide human enhancers. Hundreds of thousands of candidate enhancers have been nominated via biochemical annotations; however, only a handful of these have been validated and confidently linked to their target genes. Here we review emerging technologies for discovering, characterizing and validating human enhancers at scale. We furthermore propose a new framework for operationally defining enhancers that accommodates the heterogeneous and complementary results that are emerging from reporter assays, biochemical measurements and CRISPR screens.

Transcriptional network dynamics during the progression of pluripotency revealed by integrative statistical learning

The developmental potential of cells, termed pluripotency, is highly dynamic and progresses through a continuum of naive, formative and primed states. Pluripotency progression of mouse embryonic stem cells (ESCs) from naive to formative and primed state is governed by transcription factors (TFs) and their target genes. Genomic techniques have uncovered a multitude of TF binding sites in ESCs, yet a major challenge lies in identifying target genes from functional binding sites and reconstructing dynamic transcriptional networks underlying pluripotency progression. Here, we integrated time-resolved ‘trans-omic’ datasets together with TF binding profiles and chromatin conformation data to identify target genes of a panel of TFs. Our analyses revealed that naive TF target genes are more likely to be TFs themselves than those of formative TFs, suggesting denser hierarchies among naive TFs. We also discovered that formative TF target genes are marked by permissive epigenomic signatures in the naive state, indicating that they are poised for expression prior to the initiation of pluripotency transition to the formative state. Finally, our reconstructed transcriptional networks pinpointed the precise timing from naive to formative pluripotency progression and enabled the spatiotemporal mapping of differentiating ESCs to their in vivo counterparts in developing embryos.

Enhancers and MYC interplay in hematopoiesis

Transcription requires the fine interplay between enhancers and transcription factors. Enhancers are able to activate transcription of genes involved in normal cell biology, whereas aberrant enhancer activity leads to oncogenesis. MYC is a well-established proto-oncogene involved in half of human cancers amplifying the output of its targets. The crosstalk between MYC and enhancers is known for many years since the discovery of IgH enhancer juxtaposition with MYC in high-grade lymphomas. Here, we focus mainly in the enhancers surrounding MYC in the 8q24 locus. That region comprises several enhancers that associate with other transcription factors, transmembrane receptors, and fusion genes composing complex regulatory networks aberrantly expressed in almost all types of hematological malignancies. Understanding the nature of these interactions in normal blood cells and in leukemias/lymphomas will expand MYC targeting options in the armamentarium against hematological cancers.

Seq'ing identity and function in a repeat-derived noncoding RNA world

Innovations in high-throughout sequencing approaches are being marshaled to both reveal the composition of the abundant and heterogeneous noncoding RNAs that populate cell nuclei and lend insight to the mechanisms by which noncoding RNAs influence chromosome biology and gene expression. This review focuses on some of the recent technological developments that have enabled the isolation of nascent transcripts and chromatin-associated and DNA-interacting RNAs. Coupled with emerging genome assembly and analytical approaches, the field is poised to achieve a comprehensive catalog of nuclear noncoding RNAs, including those derived from repetitive regions within eukaryotic genomes. Herein, particular attention is paid to the challenges and advances in the sequence analyses of repeat and transposable element-derived noncoding RNAs and in ascribing specific function(s) to such RNAs.

A Novel Approach to Identify Enhancer lincRNAs by Integrating Genome, Epigenome, and Regulatome

LincRNAs enriched with high H3K4me1 and low H3K4me3 signals often have the enhancer-like features which are named as enhancer-associated lincRNAs (elincRNAs). ElincRNAs are considered to be indispensable for target gene transcription, which play important roles in development, signaling events, and even diseases. In this study, we developed a regularized regression model to identify elincRNAs by integrating the genomic, epigenomic, and regulatory data. Application of the proposed method to mouse ESCs reveals that besides the basic well-known epigenetic features H3K4me1 and H3K4me3, more specific epigenetic features, such as high DNA methylation, high H3K122ac, and H3K36me3 were contributed to mark elincRNAs with the best accuracy and precision. Finally, 3729 elincRNAs were identified in mouse ESCs. Furthermore, the elincRNAs and canonical lincRNAs exhibit distinct genomic features, and elincRNAs have the higher CGI enrichment and lower sequence conservation. Through the analysis of transcription regulation, we found that elincRNAs were significantly regulated by NANOG, POU5F1, SOX2 and ESRRB, and were involved in the core transcriptional regulatory circuitry controlling ES cell state Function enrichment analysis further discovered that elincRNAs tended to regulate specific embryonic development biological processes. These results indicated that these two types of lincRNAs had both specific epigenetic and transcriptional regulation mechanism and display distinct functional characters. In conclusion, we presented a credible computational model to prioritize novel elincRNAs, and depicted the atlas of elincRNAs in mouse ESCs, which would help dissect the function roles of lncRNAs during the mammalian development and diseases.

Diversity and Emerging Roles of Enhancer RNA in Regulation of Gene Expression and Cell Fate

Enhancers are cis-regulatory elements in the genome that cooperate with promoters to control target gene transcription. Unlike promoters, enhancers are not necessarily adjacent to target genes and can exert their functions regardless of enhancer orientations, positions and spatial segregations from target genes. Thus, for a long time, the question as to how enhancers act in a temporal and spatial manner attracted considerable attention. The recent discovery that enhancers are also abundantly transcribed raises interesting questions about the exact roles of enhancer RNA (eRNA) in gene regulation. In this review, we highlight the process of enhancer transcription and the diverse features of eRNA. We review eRNA functions, which include enhancer-promoter looping, chromatin modifying, and transcription regulating. As eRNA are transcribed from active enhancers, they exhibit tissue and lineage specificity, and serve as markers of cell state and function. Finally, we discuss the unique relationship between eRNA and super enhancers in phase separation wherein eRNA may contribute significantly to cell fate decisions.

The role of enhancer RNAs in epigenetic regulation of gene expression

Since the discovery that enhancers can support transcription, the roles of enhancer RNAs have remained largely elusive. We identified that enhancer RNAs interact with and augment bromodomain engagement with acetylated chromatin. Here, we discuss our recent findings and the potential mechanisms underlying the regulation and functions of enhancer RNA-bromodomain associations.

Coordinate regulation of ELF5 and EHF at the chr11p13 CF modifier region

E74-like factor 5 (ELF5) and ETS-homologous factor (EHF) are epithelial selective ETS family transcription factors (TFs) encoded by genes at chr11p13, a region associated with cystic fibrosis (CF) lung disease severity. EHF controls many key processes in lung epithelial function so its regulatory mechanisms are important. Using CRISPR/Cas9 technology, we removed three key cis-regulatory elements (CREs) from the chr11p13 region and also activated multiple open chromatin sites with CRISPRa in airway epithelial cells. Deletion of the CREs caused subtle changes in chromatin architecture and site-specific increases in EHF and ELF5. CRISPRa had most effect on ELF5 transcription. ELF5 levels are low in airway cells but higher in LNCaP (prostate) and T47D (breast) cancer cells. ATAC-seq in these lines revealed novel peaks of open chromatin at the 5' end of chr11p13 associated with an expressed ELF5 gene. Furthermore, 4C-seq assays identified direct interactions between the active ELF5 promoter and sites within the EHF locus, suggesting coordinate regulation between these TFs. ChIP-seq for ELF5 in T47D cells revealed ELF5 occupancy within EHF introns 1 and 6, and siRNA-mediated depletion of ELF5 enhanced EHF expression. These results define a new role for ELF5 in lung epithelial biology.

A Read/Write Mechanism Connects p300 Bromodomain Function to H2A.Z Acetylation

Acetylation of the histone variant H2A.Z (H2A.Zac) occurs at active regulatory regions associated with gene expression. Although the Tip60 complex is proposed to acetylate H2A.Z, functional studies suggest additional enzymes are involved. Here, we show that p300 acetylates H2A.Z at multiple lysines. In contrast, we found that although Tip60 does not efficiently acetylate H2A.Z in vitro, genetic inhibition of Tip60 reduces H2A.Zac in cells. Importantly, we found that interaction between the p300-bromodomain and H4 acetylation (H4ac) enhances p300-driven H2A.Zac. Indeed, H2A.Zac and H4ac show high genomic overlap, especially at active promoters. We also reveal unique chromatin features and transcriptional states at enhancers correlating with co-occurrence or exclusivity of H4ac and H2A.Zac. We propose that differential H4 and H2A.Z acetylation signatures can also define the enhancer state. In conclusion, we show both Tip60 and p300 contribute to H2A.Zac and reveal molecular mechanisms of writer/reader crosstalk between H2A.Z and H4 acetylation through p300.

•

p300 acetylates H2A.Z at multiple N-terminal lysine residues

•

Interaction of p300 bromodomain with H4ac enhances H2A.Zac

•

H2A.Zac and H4ac co-localize at active regulatory regions

•

H4ac and H2A.Zac differential signature can define the enhancer state