Tissue-specific RNA expression marks distant-acting developmental enhancers

DNA-binding factor footprints and enhancer RNAs identify functional non-coding genetic variants

BACKGROUND

Genome-wide association studies (GWAS) have revealed a multitude of candidate genetic variants affecting the risk of developing complex traits and diseases. However, the highlighted regions are typically in the non-coding genome, and uncovering the functional causative single nucleotide variants (SNVs) is challenging. Prioritization of variants is commonly based on genomic annotation with markers of active regulatory elements, but current approaches still poorly predict functional variants. To address this, we systematically analyze six markers of active regulatory elements for their ability to identify functional variants.

RESULTS

We benchmark against molecular quantitative trait loci (molQTL) from assays of regulatory element activity that identify allelic effects on DNA-binding factor occupancy, reporter assay expression, and chromatin accessibility. We identify the combination of DNase footprints and divergent enhancer RNA (eRNA) as markers for functional variants. This signature provides high precision, but with a trade-off of low recall, thus substantially reducing candidate variant sets to prioritize variants for functional validation. We present this as a framework called FINDER-Functional SNV IdeNtification using DNase footprints and eRNA.

CONCLUSIONS

We demonstrate the utility to prioritize variants using leukocyte count trait and analyze variants in linkage disequilibrium with a lead variant to predict a functional variant in asthma. Our findings have implications for prioritizing variants from GWAS, in development of predictive scoring algorithms, and for functionally informed fine mapping approaches.

CWAS-Plus: estimating category-wide association of rare noncoding variation from whole-genome sequencing data with cell-type-specific functional data

Variants in cis-regulatory elements link the noncoding genome to human pathology; however, detailed analytic tools for understanding the association between cell-level brain pathology and noncoding variants are lacking. CWAS-Plus, adapted from a Python package for category-wide association testing (CWAS), enhances noncoding variant analysis by integrating both whole-genome sequencing (WGS) and user-provided functional data. With simplified parameter settings and an efficient multiple testing correction method, CWAS-Plus conducts the CWAS workflow 50 times faster than CWAS, making it more accessible and user-friendly for researchers. Here, we used a single-nuclei assay for transposase-accessible chromatin with sequencing to facilitate CWAS-guided noncoding variant analysis at cell-type-specific enhancers and promoters. Examining autism spectrum disorder WGS data (n = 7280), CWAS-Plus identified noncoding de novo variant associations in transcription factor binding sites within conserved loci. Independently, in Alzheimer's disease WGS data (n = 1087), CWAS-Plus detected rare noncoding variant associations in microglia-specific regulatory elements. These findings highlight CWAS-Plus's utility in genomic disorders and scalability for processing large-scale WGS data and in multiple-testing corrections. CWAS-Plus and its user manual are available at https://github.com/joonan-lab/cwas/ and https://cwas-plus.readthedocs.io/en/latest/, respectively.

Genome organization and botanical diversity

The rich diversity of angiosperms, both the planet's dominant flora and the cornerstone of agriculture, is integrally intertwined with a distinctive evolutionary history. Here, we explore the interplay between angiosperm genome organization and botanical diversity, empowered by genomic approaches ranging from genetic linkage mapping to analysis of gene regulation. Commonality in the genetic hardware of plants has enabled robust comparative genomics that has provided a broad picture of angiosperm evolution and implicated both general processes and specific elements in contributing to botanical diversity. We argue that the hardware of plant genomes-both in content and in dynamics-has been shaped by selection for rather substantial differences in gene regulation between plants and animals such as maize and human, organisms of comparable genome size and gene number. Their distinctive genome content and dynamics may reflect in part the indeterminate development of plants that puts strikingly different demands on gene regulation than in animals. Repeated polyploidization of plant genomes and multiplication of individual genes together with extensive rearrangement and differential retention provide rich raw material for selection of morphological and/or physiological variations conferring fitness in specific niches, whether natural or artificial. These findings exemplify the burgeoning information available to employ in increasing knowledge of plant biology and in modifying selected plants to better meet human needs.

CWAS-Plus: Estimating category-wide association of rare noncoding variation from whole-genome sequencing data with cell-type-specific functional data

Variants in cis-regulatory elements link the noncoding genome to human brain pathology; however, detailed analytic tools for understanding the association between cell-level brain pathology and noncoding variants are lacking. CWAS-Plus, adapted from a Python package for category-wide association testing (CWAS) employs both whole-genome sequencing and user-provided functional data to enhance noncoding variant analysis, with a faster and more efficient execution of the CWAS workflow. Here, we used single-nuclei assay for transposase-accessible chromatin with sequencing to facilitate CWAS-guided noncoding variant analysis at cell-type specific enhancers and promoters. Examining autism spectrum disorder whole-genome sequencing data (n = 7,280), CWAS-Plus identified noncoding de novo variant associations in transcription factor binding sites within conserved loci. Independently, in Alzheimer's disease whole-genome sequencing data (n = 1,087), CWAS-Plus detected rare noncoding variant associations in microglia-specific regulatory elements. These findings highlight CWAS-Plus's utility in genomic disorders and scalability for processing large-scale whole-genome sequencing data and in multiple-testing corrections. CWAS-Plus and its user manual are available at https://github.com/joonan-lab/cwas/ and https://cwas-plus.readthedocs.io/en/latest/, respectively.

Constructing eRNA-mediated gene regulatory networks to explore the genetic basis of muscle and fat-relevant traits in pigs

BACKGROUND

Enhancer RNAs (eRNAs) play a crucial role in transcriptional regulation. While significant progress has been made in understanding epigenetic regulation mediated by eRNAs, research on the construction of eRNA-mediated gene regulatory networks (eGRN) and the identification of critical network components that influence complex traits is lacking.

RESULTS

Here, employing the pig as a model, we conducted a comprehensive study using H3K27ac histone ChIP-seq and RNA-seq data to construct eRNA expression profiles from multiple tissues of two distinct pig breeds, namely Enshi Black (ES) and Duroc. In addition to revealing the regulatory landscape of eRNAs at the tissue level, we developed an innovative network construction and refinement method by integrating RNA-seq, ChIP-seq, genome-wide association study (GWAS) signals and enhancer-modulating effects of single nucleotide polymorphisms (SNPs) measured by self-transcribing active regulatory region sequencing (STARR-seq) experiments. Using this approach, we unraveled eGRN that significantly influence the growth and development of muscle and fat tissues, and identified several novel genes that affect adipocyte differentiation in a cell line model.

CONCLUSIONS

Our work not only provides novel insights into the genetic basis of economic pig traits, but also offers a generalizable approach to elucidate the eRNA-mediated transcriptional regulation underlying a wide spectrum of complex traits for diverse organisms.

Identification and characterization of CHD4-associated eRNA as a novel modulator of fetal hemoglobin levels in β-thalassemia

Fetal-to-adult hemoglobin switching is controlled by programmed silencing of γ-globin while the re-activation of fetal hemoglobin (HbF) is an effective strategy for ameliorating the clinical severity of β-thalassemia and sickle cell disease. The identification of enhancer RNAs (eRNAs) related to the fetal (α2γ2) to adult hemoglobin (α2β2) switching remains incomplete. In this study, the transcriptomes of GYPA+ cells from six β-thalassemia patients with extreme HbF levels were sequenced to identify differences in patterns of noncoding RNA expression. It is interesting that an enhancer upstream of CHD4, an HbF-related core subunit of the NuRD complex, was differentially transcribed. We found a significantly positive correlation of eRNA-CHD4 enhancer-gene interaction using the public database of FANTOM5. Specifically, the eRNA-CHD4 expression was found to be significantly higher in both CD34+ HSPCs and HUDEP-2 than those in K562 cells which commonly expressed high level of HbF, suggesting a correlation between eRNA and HbF expression. Furthermore, prediction of transcription binding sites of cis-eQTLs and the CHD4 genomic region revealed a putative interaction site between rs73264846 and ZNF410, a known transcription factor regulating HbF expression. Moreover, in-vitro validation showed that the inhibition of eRNA could reduce the expression of HBG expression in HUDEP-2 cells. Taken together, the findings of this study demonstrate that a distal enhancer contributes to stage-specific silencing of γ-globin genes through direct modulation of CHD4 expression and provide insights into the epigenetic mechanisms of NuRD-mediated hemoglobin switching.

Pig-eRNAdb: a comprehensive enhancer and eRNA dataset of pigs

Enhancers and the enhancer RNAs (eRNAs) have been strongly implicated in regulations of transcriptions. Based the multi-omics data (ATAC-seq, ChIP-seq and RNA-seq) from public databases, Pig-eRNAdb is a dataset that comprehensively integrates enhancers and eRNAs for pigs using the machine learning strategy, which incorporates 82,399 enhancers and 37,803 eRNAs from 607 samples across 15 tissues of pigs. This user-friendly dataset covers a comprehensive depth of enhancers and eRNAs annotation for pigs. The coordinates of enhancers and the expression patterns of eRNAs are downloadable. Besides, thousands of regulators on eRNAs, the target genes of eRNAs, the tissue-specific eRNAs, and the housekeeping eRNAs are also accessible as well as the sequence similarity of eRNAs with humans. Moreover, the tissue-specific eRNA-trait associations encompass 652 traits are also provided. It will crucially facilitate investigations on enhancers and eRNAs with Pig-eRNAdb as a reference dataset in pigs.

Global mapping of RNA-chromatin contacts reveals a proximity-dominated connectivity model for ncRNA-gene interactions

Non-coding RNAs (ncRNAs) are transcribed throughout the genome and provide regulatory inputs to gene expression through their interaction with chromatin. Yet, the genomic targets and functions of most ncRNAs are unknown. Here we use chromatin-associated RNA sequencing (ChAR-seq) to map the global network of ncRNA interactions with chromatin in human embryonic stem cells and the dynamic changes in interactions during differentiation into definitive endoderm. We uncover general principles governing the organization of the RNA-chromatin interactome, demonstrating that nearly all ncRNAs exclusively interact with genes in close three-dimensional proximity to their locus and provide a model predicting the interactome. We uncover RNAs that interact with many loci across the genome and unveil thousands of unannotated RNAs that dynamically interact with chromatin. By relating the dynamics of the interactome to changes in gene expression, we demonstrate that activation or repression of individual genes is unlikely to be controlled by a single ncRNA.

Identification and characterization of transcribed enhancers during cerebellar development through enhancer RNA analysis

BACKGROUND

The development of the brain requires precise coordination of molecular processes across many cell-types. Underpinning these events are gene expression programs which require intricate regulation by non-coding regulatory sequences known as enhancers. In the context of the developing brain, transcribed enhancers (TEs) regulate temporally-specific expression of genes critical for cell identity and differentiation. Transcription of non-coding RNAs at active enhancer sequences, known as enhancer RNAs (eRNAs), is tightly associated with enhancer activity and has been correlated with target gene expression. TEs have been characterized in a multitude of developing tissues, however their regulatory role has yet to be described in the context of embryonic and early postnatal brain development. In this study, eRNA transcription was analyzed to identify TEs active during cerebellar development, as a proxy for the developing brain. Cap Analysis of Gene Expression followed by sequencing (CAGE-seq) was conducted at 12 stages throughout embryonic and early postnatal cerebellar development.

RESULTS

Temporal analysis of eRNA transcription identified clusters of TEs that peak in activity during either embryonic or postnatal times, highlighting their importance for temporally specific developmental events. Functional analysis of putative target genes identified molecular mechanisms under TE regulation revealing that TEs regulate genes involved in biological processes specific to neurons. We validate enhancer activity using in situ hybridization of eRNA expression from TEs predicted to regulate Nfib, a gene critical for cerebellar granule cell differentiation.

CONCLUSIONS

The results of this analysis provide a valuable dataset for the identification of cerebellar enhancers and provide insight into the molecular mechanisms critical for brain development under TE regulation. This dataset is shared with the community through an online resource ( https://goldowitzlab.shinyapps.io/trans-enh-app/ ).

Long non-coding RNAs: definitions, functions, challenges and recommendations

Genes specifying long non-coding RNAs (lncRNAs) occupy a large fraction of the genomes of complex organisms. The term 'lncRNAs' encompasses RNA polymerase I (Pol I), Pol II and Pol III transcribed RNAs, and RNAs from processed introns. The various functions of lncRNAs and their many isoforms and interleaved relationships with other genes make lncRNA classification and annotation difficult. Most lncRNAs evolve more rapidly than protein-coding sequences, are cell type specific and regulate many aspects of cell differentiation and development and other physiological processes. Many lncRNAs associate with chromatin-modifying complexes, are transcribed from enhancers and nucleate phase separation of nuclear condensates and domains, indicating an intimate link between lncRNA expression and the spatial control of gene expression during development. lncRNAs also have important roles in the cytoplasm and beyond, including in the regulation of translation, metabolism and signalling. lncRNAs often have a modular structure and are rich in repeats, which are increasingly being shown to be relevant to their function. In this Consensus Statement, we address the definition and nomenclature of lncRNAs and their conservation, expression, phenotypic visibility, structure and functions. We also discuss research challenges and provide recommendations to advance the understanding of the roles of lncRNAs in development, cell biology and disease.

Dynamics of histone acetylation during human early embryogenesis

It remains poorly understood about the regulation of gene and transposon transcription during human early embryogenesis. Here, we report that broad H3K27ac domains are genome-widely distributed in human 2-cell and 4-cell embryos and transit into typical peaks in the 8-cell embryos. The broad H3K27ac domains in early embryos before zygotic genome activation (ZGA) are also observed in mouse. It suggests that broad H3K27ac domains play conserved functions before ZGA in mammals. Intriguingly, a large portion of broad H3K27ac domains overlap with broad H3K4me3 domains. Further investigation reveals that histone deacetylases are required for the removal or transition of broad H3K27ac domains and ZGA. After ZGA, the number of typical H3K27ac peaks is dynamic, which is associated with the stage-specific gene expression. Furthermore, P300 is important for the establishment of H3K27ac peaks and the expression of associated genes in early embryos after ZGA. Our data also indicate that H3K27ac marks active transposons in early embryos. Interestingly, H3K27ac and H3K18ac signals rather than H3K9ac signals are enriched at ERVK elements in mouse embryos after ZGA. It suggests that different types of histone acetylations exert distinct roles in the activation of transposons. In summary, H3K27ac modification undergoes extensive reprogramming during early embryo development in mammals, which is associated with the expression of genes and transposons.

Cell-specific expression of the FAP gene is regulated by enhancer elements

Fibroblast activation protein (FAP) is an integral membrane serine protease that acts as both dipeptidyl peptidase and collagenase. In recent years, FAP has attracted considerable attention due to its specific upregulation in multiple types of tumor cell populations, including cancer cells in various cancer types, making FAP a potential target for therapy. However, relatively few papers pay attention to the mechanisms driving the cell-specific expression of the FAP gene. We found no correlation between the activities of the two FAP promoter variants (short and long) and the endogenous FAP mRNA expression level in several cell lines with different FAP expression levels. This suggested that other mechanisms may be responsible for specific transcriptional regulation of the FAP gene. We analyzed the distribution of known epigenetic and structural chromatin marks in FAP-positive and FAP-negative cell lines and identified two potential enhancer-like elements (E1 and E2) in the FAP gene locus. We confirmed the specific enrichment of H3K27ac in the putative enhancer regions in FAP-expressing cells. Both the elements exhibited enhancer activity independently of each other in the functional test by increasing the activity of the FAP promoter variants to a greater extent in FAP-expressing cell lines than in FAP-negative cell lines. The transcription factors AP-1, CEBPB, and STAT3 may be involved in FAP activation in the tumors. We hypothesized the existence of a positive feedback loop between FAP and STAT3, which may have implications for developing new approaches in cancer therapy.

The integrated landscape of eRNA in gastric cancer reveals distinct immune subtypes with prognostic and therapeutic relevance

The comprehensive regulation effect of eRNA on tumor immune cell infiltration and the outcome remains obscure. We comprehensively identify the eRNA-mediated immune infiltration patterns of gastric cancer (GC) samples. We creatively proposed a random forest machine-learning (ML) algorithm to map eRNA to mRNA expression patterns. The eRNA score was constructed using principal component analysis algorithms and validated in an independent cohort. Three subtypes with distinct eRNA expression patterns were determined in GC. There were significant differences between the three subtypes in the overall survival rate, immune cell infiltration characteristics, and immunotherapy response indicators. The patients in the high eRNA score group have a higher overall survival rate and might benefit from immunotherapy. This work revealed that eRNA regulation might be a new prognostic index and might offer a potential biomarker in the response of immunotherapy. Evaluating the eRNA regulation manner of GC will contribute to guiding more effective immunotherapy strategies.

Oct4 differentially regulates chromatin opening and enhancer transcription in pluripotent stem cells

The transcription factor Oct4 is essential for the maintenance and induction of stem cell pluripotency, but its functional roles are not fully understood. Here, we investigate the functions of Oct4 by depleting and subsequently recovering it in mouse embryonic stem cells (ESCs) and conducting a time-resolved multiomics analysis. Oct4 depletion leads to an immediate loss of its binding to enhancers, accompanied by a decrease in mRNA synthesis from its target genes that are part of the transcriptional network that maintains pluripotency. Gradual decrease of Oct4 binding to enhancers does not immediately change the chromatin accessibility but reduces transcription of enhancers. Conversely, partial recovery of Oct4 expression results in a rapid increase in chromatin accessibility, whereas enhancer transcription does not fully recover. These results indicate different concentration-dependent activities of Oct4. Whereas normal ESC levels of Oct4 are required for transcription of pluripotency enhancers, low levels of Oct4 are sufficient to retain chromatin accessibility, likely together with other factors such as Sox2.

Pregnancy, preeclampsia and maternal aging: From epidemiology to functional genomics

Women live longer than men but experience greater disability and a longer period of illness as they age. Despite clear sex differences in aging, the impact of pregnancy and its complications, such as preeclampsia, on aging is an underexplored area of geroscience. This review summarizes our current knowledge about the complex links between pregnancy and age-related diseases, including evidence from epidemiology, clinical research, and genetics. We discuss the relationship between normal and pathological pregnancy and maternal aging, using preeclampsia as a primary example. We review the results of human genetics studies of preeclampsia, including genome wide association studies (GWAS), and attempted to catalog genes involved in preeclampsia as a gateway to mechanisms underlying an increased risk of later life cardio- and neuro- vascular events. Lastly, we discuss challenges in interpreting the GWAS of preeclampsia and provide a functional genomics framework for future research needed to fully realize the promise of GWAS in identifying targets for geroprotective prevention and therapeutics against preeclampsia.

Integrative Analysis of Epigenome and Transcriptome Data Reveals Aberrantly Methylated Promoters and Enhancers in Hepatocellular Carcinoma

DNA methylation is a key transcription regulator, whose aberration was ubiquitous and important in most cancers including hepatocellular carcinoma (HCC). Whole-genome bisulfite sequencing (WGBS) was conducted for comparison of DNA methylation in tumor and adjacent tissues from 33 HCC patients, accompanying RNA-seq to determine differentially methylated region-associated, differentially expressed genes (DMR-DEGs), which were independently replicated in the TCGA-LIHC cohort and experimentally validated via 5-aza-2-deoxycytidine (5-azadC) demethylation. A total of 9,867,700 CpG sites showed significantly differential methylation in HCC. Integrations of mRNA-seq, histone ChIP-seq, and WGBS data identified 611 high-confidence DMR-DEGs. Enrichment analysis demonstrated activation of multiple molecular pathways related to cell cycle and DNA repair, accompanying repression of several critical metabolism pathways such as tyrosine and monocarboxylic acid metabolism. In TCGA-LIHC, we replicated about 53% of identified DMR-DEGs and highlighted the prognostic significance of combinations of methylation and expression of nine DMR-DEGs, which were more efficient prognostic biomarkers than considering either type of data alone. Finally, we validated 22/23 (95.7%) DMR-DEGs in 5-azadC-treated LO2 and/or HepG2 cells. In conclusion, integration of epigenome and transcriptome data depicted activation of multiple pivotal cell cycle-related pathways and repression of several metabolic pathways triggered by aberrant DNA methylation of promoters and enhancers in HCC.

Identification and Functional Characterization of Two Noncoding RNAs Transcribed from Putative Active Enhancers in Hepatocellular Carcinoma

Enhancers have been conventionally perceived as cis-acting elements that provide binding sites for trans-acting factors. However, recent studies have shown that enhancers are transcribed and that these transcripts, called enhancer RNAs (eRNAs), have a regulatory function. Here, we identified putative eRNAs by profiling and determining the overlap between noncoding RNA expression loci and eRNA-associated histone marks such as H3K27ac and H3K4me1 in hepatocellular carcinoma (HCC) cell lines. Of the 132 HCC-derived noncoding RNAs, 74 overlapped with the eRNA loci defined by the FANTOM consortium, and 65 were located in the proximal regions of genes differentially expressed between normal and tumor tissues in TCGA dataset. Interestingly, knockdown of two selected putative eRNAs, THUMPD3-AS1 and LINC01572, led to downregulation of their target mRNAs and to a reduction in the proliferation and migration of HCC cells. Additionally, the expression of these two noncoding RNAs and target mRNAs was elevated in tumor samples in the TCGA dataset, and high expression was associated with poor survival of patients. Collectively, our study suggests that noncoding RNAs such as THUMPD3-AS1 and LINC01572 (i.e., putative eRNAs) can promote the transcription of genes involved in cell proliferation and differentiation and that the dysregulation of these noncoding RNAs can cause cancers such as HCC.

Current Advances on the Important Roles of Enhancer RNAs in Molecular Pathways of Cancer

Enhancers are critical genomic elements that can cooperate with promoters to regulate gene transcription in both normal and cancer cells. Recent studies reveal that enhancer regions are transcribed to produce a class of noncoding RNAs referred to as enhancer RNAs (eRNAs). Emerging evidence shows that eRNAs play important roles in enhancer activation and enhancer-driven gene regulation, and the expression of eRNAs may be a critical factor in tumorigenesis. The important roles of eRNAs in cancer signaling pathways are also gradually unveiled, providing a new insight into cancer therapy. Here, we review the roles of eRNAs in regulating cancer signaling pathways and discuss the potential of eRNA-targeted therapy for human cancers.

Nuclear Dishevelled targets gene regulatory regions and promotes tumor growth

Dishevelled (DVL) critically regulates Wnt signaling and contributes to a wide spectrum of diseases and is important in normal and pathophysiological settings. However, how it mediates diverse cellular functions remains poorly understood. Recent discoveries have revealed that constitutive Wnt pathway activation contributes to breast cancer malignancy, but the mechanisms by which this occurs are unknown and very few studies have examined the nuclear role of DVL. Here, we have performed DVL3 ChIP-seq analyses and identify novel target genes bound by DVL3. We show that DVL3 depletion alters KMT2D binding to novel targets and changes their epigenetic marks and mRNA levels. We further demonstrate that DVL3 inhibition leads to decreased tumor growth in two different breast cancer models in vivo. Our data uncover new DVL3 functions through its regulation of multiple genes involved in developmental biology, antigen presentation, metabolism, chromatin remodeling, and tumorigenesis. Overall, our study provides unique insight into the function of nuclear DVL, which helps to define its role in mediating aberrant Wnt signaling.

An interdependent network of functional enhancers regulates transcription and EZH2 loading at the INK4a/ARF locus

The INK4a/ARF locus encodes important cell-cycle regulators p14ARF, p15INK4b, and p16INK4a. The neighboring gene desert to this locus is the most reproducible GWAS hotspot that harbors one of the densest enhancer clusters in the genome. However, how multiple enhancers that overlap with GWAS variants regulate the INK4a/ARF locus is unknown, which is an important step in linking genetic variation with associated diseases. Here, we show that INK4a/ARF promoters interact with a subset of enhancers in the cluster, independent of their H3K27ac and eRNA levels. Interacting enhancers transcriptionally control each other and INK4a/ARF promoters over long distances as an interdependent single unit. The deletion of even a single interacting enhancer results in an unexpected collapse of the entire enhancer cluster and leads to EZH2 enrichment on promoters in an ANRIL-independent manner. Dysregulated genes genome-wide mimic 9p21-associated diseases under these scenarios. Our results highlight intricate dependencies of promoter-interacting enhancers on each other.

Characterization of transcripts emanating from enhancer Eβ of the murine TCRβ locus

Enhancers are well established as critical regulators of gene expression, but the mechanisms underlying the molecular basis of their specificity and activity are only partly understood. One of the most exciting recent observations is the discovery of enhancer RNA (eRNA), a class of noncoding RNAs derived from enhancer regions. Transcription of developmentally regulated enhancers has been observed to be associated with their active state. The nature of transcripts (eRNA) and their functional attributes are diverse and context dependent. The majority of eRNA are nonpolyadenylated and present in low abundance owing to their low stability, and may represent transcriptional noise. However, some eRNAs have been reported to be reasonably long and stable, are enriched in nuclei, exhibit tissue-specific expression and may contribute to enhancer function. Transcription of enhancers has been postulated to mediate enhancer function through either the act of transcription or via the transcribed RNA per se and is a useful feature to be analysed to understand mechanisms underlying enhancer activity. Enhancer Eβ at the murine TCRβ locus has been reported to exhibit enhanced occupancy of RNA polymerase II in developing thymocytes. Here, we investigated the transcriptional potential of Eβ in developing thymocytes and detected overlapping bidirectional transcripts at Eβ ranging between 0.7 and 1.7 kb. These noncoding transcripts are capped, polyadenylated, nuclear and expressed specifically in thymocytes. Delineation of these characteristics is important to further investigate their functional roles in mediating enhancer activity.

Intronic enhancer region governs transcript-specific Bdnf expression in rodent neurons

Brain-derived neurotrophic factor (BDNF) controls the survival, growth, and function of neurons both during the development and in the adult nervous system. Bdnf is transcribed from several distinct promoters generating transcripts with alternative 5' exons. Bdnf transcripts initiated at the first cluster of exons have been associated with the regulation of body weight and various aspects of social behavior, but the mechanisms driving the expression of these transcripts have remained poorly understood. Here, we identify an evolutionarily conserved intronic enhancer region inside the Bdnf gene that regulates both basal and stimulus-dependent expression of the Bdnf transcripts starting from the first cluster of 5' exons in mouse and rat neurons. We further uncover a functional E-box element in the enhancer region, linking the expression of Bdnf and various pro-neural basic helix–loop–helix transcription factors. Collectively, our results shed new light on the cell-type- and stimulus-specific regulation of the important neurotrophic factor BDNF.

Transcription activation depends on the length of the RNA polymerase II C-terminal domain

Eukaryotic RNA polymerase II (Pol II) contains a tail‐like, intrinsically disordered carboxy‐terminal domain (CTD) comprised of heptad‐repeats, that functions in coordination of the transcription cycle and in coupling transcription to co‐transcriptional processes. The CTD repeat number varies between species and generally increases with genome size, but the reasons for this are unclear. Here, we show that shortening the CTD in human cells to half of its length does not generally change pre‐mRNA synthesis or processing in cells. However, CTD shortening decreases the duration of promoter‐proximal Pol II pausing, alters transcription of putative enhancer elements, and delays transcription activation after stimulation of the MAP kinase pathway. We suggest that a long CTD is required for efficient enhancer‐dependent recruitment of Pol II to target genes for their rapid activation.

Transcriptionally active enhancers in human cancer cells

The growth of human cancer cells is driven by aberrant enhancer and gene transcription activity. Here, we use transient transcriptome sequencing (TT-seq) to map thousands of transcriptionally active putative enhancers in fourteen human cancer cell lines covering seven types of cancer. These enhancers were associated with cell type-specific gene expression, enriched for genetic variants that predispose to cancer, and included functionally verified enhancers. Enhancer-promoter (E-P) pairing by correlation of transcription activity revealed ~ 40,000 putative E-P pairs, which were depleted for housekeeping genes and enriched for transcription factors, cancer-associated genes, and 3D conformational proximity. The cell type specificity and transcription activity of target genes increased with the number of paired putative enhancers. Our results represent a rich resource for future studies of gene regulation by enhancers and their role in driving cancerous cell growth.

eRNAs and Superenhancer lncRNAs Are Functional in Human Prostate Cancer

Prostate cancer (PCa) is one of the most commonly diagnosed cancers in males worldwide. lncRNAs (long noncoding RNAs) play a significant role in the occurrence and development of PCa. eRNAs (enhancer RNAs) and SE-lncRNAs (superenhancer lncRNAs) are important elements of lncRNAs, but the role of eRNAs and SE-lncRNAs in PCa remains largely unclear. In this work, we identified 681 eRNAs and 292 SE-lncRNAs that were expressed differentially in PCa using a microarray. We also found that eRNAs transcribed from active open chromatin had significantly higher expression than those from active closed chromatin, and SE-lncRNAs had a little higher expression than eRNAs. Next, we constructed a transcriptional regulation network that eRNA-related enhancer and the target genes shared the same TF-binding motifs. Further, we investigated whether CTCF played a role in mediating the transcriptional regulation network. eRNAs, especially those that regulate androgen response genes, may be candidates for prognostic biomarkers and therapy targets. Our work provides a new perspective for developing medical treatments and therapies for prostate cancer.

Versatile interactions and bioinformatics analysis of noncoding RNAs

Advances in RNA sequencing technologies and computational methodologies have provided a huge impetus to noncoding RNA (ncRNA) study. Once regarded as inconsequential results of transcriptional promiscuity, ncRNAs were later found to exert great roles in various aspects of biological functions. They are emerging as key players in gene regulatory networks by interacting with other biomolecules (DNA, RNA or protein). Here, we provide an overview of ncRNA repertoire and highlight recent discoveries of their versatile interactions. To better investigate the ncRNA-mediated regulation, it is necessary to make full use of innovative sequencing techniques and computational tools. We further describe a comprehensive workflow for in silico ncRNA analysis, providing up-to-date platforms, databases and tools dedicated to ncRNA identification and functional annotation.

Total Functional Score of Enhancer Elements Identifies Lineage-Specific Enhancers That Drive Differentiation of Pancreatic Cells

The differentiation of embryonic stem cells into various lineages is highly dependent on the chromatin state of the genome and patterns of gene expression. To identify lineage-specific enhancers driving the differentiation of progenitors into pancreatic cells, we used a previously described computational framework called Total Functional Score of Enhancer Elements (TFSEE), which integrates multiple genomic assays that probe both transcriptional and epigenomic states. First, we evaluated and compared TFSEE as an enhancer-calling algorithm with enhancers called using GRO-seq-defined enhancer transcripts (method 1) versus enhancers called using histone modification ChIP-seq data (method 2). Second, we used TFSEE to define the enhancer landscape and identify transcription factors (TFs) that maintain the multipotency of a subpopulation of endodermal stem cells during differentiation into pancreatic lineages. Collectively, our results demonstrate that TFSEE is a robust enhancer-calling algorithm that can be used to perform multilayer genomic data integration to uncover cell type-specific TFs that control lineage-specific enhancers.

Blood disease-causing and -suppressing transcriptional enhancers: general principles and GATA2 mechanisms

Intensive scrutiny of human genomes has unveiled considerable genetic variation in coding and noncoding regions. In cancers, including those of the hematopoietic system, genomic instability amplifies the complexity and functional consequences of variation. Although elucidating how variation impacts the protein-coding sequence is highly tractable, deciphering the functional consequences of variation in noncoding regions (genome reading), including potential transcriptional-regulatory sequences, remains challenging. A crux of this problem is the sheer abundance of gene-regulatory sequence motifs (cis elements) mediating protein-DNA interactions that are intermixed in the genome with thousands of look-alike sequences lacking the capacity to mediate functional interactions with proteins in vivo. Furthermore, transcriptional enhancers harbor clustered cis elements, and how altering a single cis element within a cluster impacts enhancer function is unpredictable. Strategies to discover functional enhancers have been innovated, and human genetics can provide vital clues to achieve this goal. Germline or acquired mutations in functionally critical (essential) enhancers, for example at the GATA2 locus encoding a master regulator of hematopoiesis, have been linked to human pathologies. Given the human interindividual genetic variation and complex genetic landscapes of hematologic malignancies, enhancer corruption, creation, and expropriation by new genes may not be exceedingly rare mechanisms underlying disease predisposition and etiology. Paradigms arising from dissecting essential enhancer mechanisms can guide genome-reading strategies to advance fundamental knowledge and precision medicine applications. In this review, we provide our perspective of general principles governing the function of blood disease-linked enhancers and GATA2-centric mechanisms.

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.

Epigenetic and Transcriptional Networks Underlying Atrial Fibrillation

Genome-wide association studies have uncovered over a 100 genetic loci associated with atrial fibrillation (AF), the most common arrhythmia. Many of the top AF-associated loci harbor key cardiac transcription factors, including PITX2, TBX5, PRRX1, and ZFHX3. Moreover, the vast majority of the AF-associated variants lie within noncoding regions of the genome where causal variants affect gene expression by altering the activity of transcription factors and the epigenetic state of chromatin. In this review, we discuss a transcriptional regulatory network model for AF defined by effector genes in Genome-wide association studies loci. We describe the current state of the field regarding the identification and function of AF-relevant gene regulatory networks, including variant regulatory elements, dose-sensitive transcription factor functionality, target genes, and epigenetic states. We illustrate how altered transcriptional networks may impact cardiomyocyte function and ionic currents that impact AF risk. Last, we identify the need for improved tools to identify and functionally test transcriptional components to define the links between genetic variation, epigenetic gene regulation, and atrial function.

Characterization of Arabidopsis thaliana Promoter Bidirectionality and Antisense RNAs by Inactivation of Nuclear RNA Decay Pathways

In animals, RNA polymerase II initiates transcription bidirectionally from gene promoters to produce pre-mRNAs on the forward strand and promoter upstream transcripts (PROMPTs) on the reverse strand. PROMPTs are degraded by the nuclear exosome. Previous studies based on nascent RNA approaches concluded that Arabidopsis (Arabidopsis thaliana) does not produce PROMPTs. Here, we used steady-state RNA sequencing in mutants defective in nuclear RNA decay including the exosome to reassess the existence of Arabidopsis PROMPTs. While they are rare, we identified ∼100 cases of exosome-sensitive PROMPTs in Arabidopsis. Such PROMPTs are sources of small interfering RNAs in exosome-deficient mutants, perhaps explaining why plants have evolved mechanisms to suppress PROMPTs. In addition, we found ∼200 long, unspliced and exosome-sensitive antisense RNAs that arise from transcription start sites within parts of the genome encoding 3'-untranslated regions on the sense strand. The previously characterized noncoding RNA that regulates expression of the key seed dormancy regulator, DELAY OF GERMINATION1, is a typical representative of this class of RNAs. Transcription factor genes are overrepresented among loci with exosome-sensitive antisense RNAs, suggesting a potential for widespread control of gene expression via this class of noncoding RNAs. Lastly, we assess the use of alternative promoters in Arabidopsis and compare the accuracy of existing TSS annotations.

Resolving mechanisms of immune-mediated disease in primary CD4 T cells

Deriving mechanisms of immune-mediated disease from GWAS data remains a formidable challenge, with attempts to identify causal variants being frequently hampered by strong linkage disequilibrium. To determine whether causal variants could be identified from their functional effects, we adapted a massively parallel reporter assay for use in primary CD4 T cells, the cell type whose regulatory DNA is most enriched for immune-mediated disease SNPs. This enabled the effects of candidate SNPs to be examined in a relevant cellular context and generated testable hypotheses into disease mechanisms. To illustrate the power of this approach, we investigated a locus that has been linked to six immune-mediated diseases but cannot be fine-mapped. By studying the lead expression-modulating SNP, we uncovered an NF-κB-driven regulatory circuit which constrains T-cell activation through the dynamic formation of a super-enhancer that upregulates TNFAIP3 (A20), a key NF-κB inhibitor. In activated T cells, this feedback circuit is disrupted-and super-enhancer formation prevented-by the risk variant at the lead SNP, leading to unrestrained T-cell activation via a molecular mechanism that appears to broadly predispose to human autoimmunity.

Evidence for craniofacial enhancer variation underlying nonsyndromic cleft lip and palate

Nonsyndromic cleft lip with or without cleft palate (NSCLP) is a common birth defect for which only ~ 20% of the underlying genetic variation has been identified. Variants in noncoding regions have been increasingly suggested to contribute to the missing heritability. In this study, we investigated whether variation in craniofacial enhancers contributes to NSCLP. Candidate enhancers were identified using VISTA Enhancer Browser and previous publications. Prioritization was based on patterning defects in knockout mice, deletion/duplication of craniofacial genes in animal models and results of whole exome/whole genome sequencing studies. This resulted in 20 craniofacial enhancers to be investigated. Custom amplicon-based sequencing probes were designed and used for sequencing 380 NSCLP probands (from multiplex and simplex families of non-Hispanic white (NHW) and Hispanic ethnicities) using Illumina MiSeq. The frequencies of identified variants were compared to ethnically matched European (CEU) and Los Angeles Mexican (MXL) control genomes and used for association analyses. Variants in mm427/MSX1 and hs1582/SPRY1 showed genome-wide significant association with NSCLP (p ≤ 6.4 × 10^-11). In silico analysis showed that these enhancer variants may disrupt important transcription factor binding sites. Haplotypes involving these enhancers and also mm435/ABCA4 were significantly associated with NSCLP, especially in NHW (p ≤ 6.3 × 10^-7). Importantly, groupwise burden analysis showed several enhancer combinations significantly over-represented in NSCLP individuals, revealing novel NSCLP pathways and supporting a polygenic inheritance model. Our findings support the role of craniofacial enhancer sequence variation in the etiology of NSCLP.

Millefy: visualizing cell-to-cell heterogeneity in read coverage of single-cell RNA sequencing datasets

BACKGROUND

Read coverage of RNA sequencing data reflects gene expression and RNA processing events. Single-cell RNA sequencing (scRNA-seq) methods, particularly "full-length" ones, provide read coverage of many individual cells and have the potential to reveal cellular heterogeneity in RNA transcription and processing. However, visualization tools suited to highlighting cell-to-cell heterogeneity in read coverage are still lacking.

RESULTS

Here, we have developed Millefy, a tool for visualizing read coverage of scRNA-seq data in genomic contexts. Millefy is designed to show read coverage of all individual cells at once in genomic contexts and to highlight cell-to-cell heterogeneity in read coverage. By visualizing read coverage of all cells as a heat map and dynamically reordering cells based on diffusion maps, Millefy facilitates discovery of "local" region-specific, cell-to-cell heterogeneity in read coverage. We applied Millefy to scRNA-seq data sets of mouse embryonic stem cells and triple-negative breast cancers and showed variability of transcribed regions including antisense RNAs, 3 ' UTR lengths, and enhancer RNA transcription.

CONCLUSIONS

Millefy simplifies the examination of cellular heterogeneity in RNA transcription and processing events using scRNA-seq data. Millefy is available as an R package (https://github.com/yuifu/millefy) and as a Docker image for use with Jupyter Notebook (https://hub.docker.com/r/yuifu/datascience-notebook-millefy).

Splicing of enhancer-associated lincRNAs contributes to enhancer activity

Transcription is common at active mammalian enhancers sometimes giving rise to stable enhancer-associated long intergenic noncoding RNAs (elincRNAs). Expression of elincRNA is associated with changes in neighboring gene product abundance and local chromosomal topology, suggesting that transcription at these loci contributes to gene expression regulation in cis Despite the lack of evidence supporting sequence-dependent functions for most elincRNAs, splicing of these transcripts is unexpectedly common. Whether elincRNA splicing is a mere consequence of cognate enhancer activity or if it directly impacts enhancer function remains unresolved. Here, we investigate the association between elincRNA splicing and enhancer activity in mouse embryonic stem cells. We show that multi-exonic elincRNAs are enriched at conserved enhancers, and the efficient processing of elincRNAs is strongly associated with their cognate enhancer activity. This association is supported by their enrichment in enhancer-specific chromatin signatures; elevated binding of co-transcriptional regulators; increased local intra-chromosomal DNA contacts; and strengthened cis-regulation on target gene expression. Our results support the role of efficient RNA processing of enhancer-associated transcripts to cognate enhancer activity.

Enhancer transcription identifies cis-regulatory elements for photoreceptor cell types

Identification of cell type-specific cis-regulatory elements (CREs) is crucial for understanding development and disease, although identification of functional regulatory elements remains challenging. We hypothesized that context-specific CREs could be identified by context-specific non-coding RNA (ncRNA) profiling, based on the observation that active CREs produce ncRNAs. We applied ncRNA profiling to identify rod and cone photoreceptor CREs from wild-type and mutant mouse retinas, defined by presence or absence, respectively, of the rod-specific transcription factor (TF) NrlNrl-dependent ncRNA expression strongly correlated with epigenetic profiles of rod and cone photoreceptors, identified thousands of candidate rod- and cone-specific CREs, and identified motifs for rod- and cone-specific TFs. Colocalization of NRL and the retinal TF CRX correlated with rod-specific ncRNA expression, whereas CRX alone favored cone-specific ncRNA expression, providing quantitative evidence that heterotypic TF interactions distinguish cell type-specific CRE activity. We validated the activity of novel Nrl-dependent ncRNA-defined CREs in developing cones. This work supports differential ncRNA profiling as a platform for the identification of cell type-specific CREs and the discovery of molecular mechanisms underlying TF-dependent CRE activity.

Invited Review: Epigenetics in neurodevelopment

Neural development requires the orchestration of dynamic changes in gene expression to regulate cell fate decisions. This regulation is heavily influenced by epigenetics, heritable changes in gene expression not directly explained by genomic information alone. An understanding of the complexity of epigenetic regulation is rapidly emerging through the development of novel technologies that can assay various features of epigenetics and gene regulation. Here, we provide a broad overview of several commonly investigated modes of epigenetic regulation, including DNA methylation, histone modifications, noncoding RNAs, as well as epitranscriptomics that describe modifications of RNA, in neurodevelopment and diseases. Rather than functioning in isolation, it is being increasingly appreciated that these various modes of gene regulation are dynamically interactive and coordinate the complex nature of neurodevelopment along multiple axes. Future work investigating these interactions will likely utilize 'multi-omic' strategies that assay cell fate dynamics in a high-dimensional and high-throughput fashion. Novel human neurodevelopmental models including iPSC and cerebral organoid systems may provide further insight into human-specific features of neurodevelopment 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 therapeutic and diagnostic potential of regulatory noncoding RNAs in medulloblastoma

Medulloblastoma, a central nervous system tumor that predominantly affects children, always requires aggressive therapy. Nevertheless, it frequently recurs as resistant disease and is associated with high morbidity and mortality. While recent efforts to subclassify medulloblastoma based on molecular features have advanced our basic understanding of medulloblastoma pathogenesis, optimal targets to increase therapeutic efficacy and reduce side effects remain largely undefined. Noncoding RNAs (ncRNAs) with known regulatory roles, particularly long noncoding RNAs (lncRNAs) and microRNAs (miRNAs), are now known to participate in medulloblastoma biology, although their functional significance remains obscure in many cases. Here we review the literature on regulatory ncRNAs in medulloblastoma. In providing a comprehensive overview of ncRNA studies, we highlight how different lncRNAs and miRNAs have oncogenic or tumor suppressive roles in medulloblastoma. These ncRNAs possess subgroup specificity that can be exploited to personalize therapy by acting as theranostic targets. Several of the already identified ncRNAs appear specific to medulloblastoma stem cells, the most difficult-to-treat component of the tumor that drives metastasis and acquired resistance, thereby providing opportunities for therapy in relapsing, disseminating, and therapy-resistant disease. Delivering ncRNAs to tumors remains challenging, but this limitation is gradually being overcome through the use of advanced technologies such as nanotechnology and rational biomaterial design.

Landscape of Enhancer-Enhancer Cooperative Regulation during Human Cardiac Commitment

Although accumulating evidence has demonstrated the key roles of enhancers in gene expression regulation, the contribution of genome-wide enhancer-enhancer interactions to developmental decisions remains unclear. Here we explored the cooperative regulation patterns among enhancers to understand their regulatory mechanism. We first filtered robust enhancers in embryonic stem cells (ESCs) through integrating bidirectional transcription, genomic location, and epigenetic modification. Genome-wide enhancer-enhancer interactions were then identified based on enhancer-promoter relationships that were derived from Hi-C data. We further explored the interacting principles of the identified enhancer-enhancer interactions. The results revealed that the observed cooperativity occurred mainly between enhancers distributed within a 1-kb to 10-Mb distance across the genome. In addition, enhancer-enhancer pairs had higher expression correlations than non-interacting pairs. Finally, we identified robust enhancers during human cardiac commitment, and we found that enhancers exhibited strong stage-specific expression patterns. We further inferred the enhancer-enhancer interactions based on RNA sequencing (RNA-seq) data in heart development, according to the regulatory principles characterized from Hi-C data. The identified enhancer-enhancer interaction networks (EEINs) presented highly dynamic linkages. Moreover, enhancers cooperatively targeted many marker genes in each developmental stage to regulate stage-specific functions, which contribute to the organization of cell identity in heart development. Our work will increase the understanding of enhancer regulation in human heart development.

The evolution of Great Apes has shaped the functional enhancers' landscape in human embryonic stem cells

High-throughput functional assays of enhancer activity have recently enabled the genome-scale definition of molecular, structural, and biochemical features of these genomic regulatory regions. To infer the evolutionary origin of DNA sequences operating as functional enhancers in human embryonic stem cells (hESC), we examined the patterns of evolutionary conservation and divergence in the genome-wide functional enhancers' landscape of hESC. We show that a prominent majority (up to 94%) of DNA sequences identified in hESC as functional enhancers are conserved in humans and our closest evolutionary relatives, Chimpanzee and Bonobo. More than 91% of functional enhancers that are highly conserved in both Chimpanzee and Bonobo, are conserved among other Great Apes and >75% are conserved in the Rhesus genome. In striking contrast, <5% of DNA sequences operating in hESC as functional enhancers are conserved in rodents. Conserved in primates enhancers' sequences are complemented by 1619 sequences of enhancers that are specific to humans. Enhancers that harbor human-specific sequences appear enriched among the invariant enhancer module maintaining activity in different pluripotent states and these regions are associated with pluripotency- and embryonic-lineage-related genes. However, functional enhancers make up only a minority of all conserved in primates or human-specific transcription factor binding sites. Our analyses revealed that sequences that are conserved during ~8 million years of primate evolution dominate the genomic landscape of functional enhancers in both primed and naïve hESC. Collectively, these observations revealed thousands of evolutionarily conserved sequences that function as a core regulatory network in human embryonic stem cells which has recently undergone further extension after divergence of modern humans from our closest relatives, Chimpanzee and Bonobo.

RNA element discovery from germ cell to blastocyst

Recent studies have shown that tissue-specific transcriptomes contain multiple types of RNAs that are transcribed from intronic and intergenic sequences. The current study presents a tool for the discovery of transcribed, unannotated sequence elements from RNA-seq libraries. This RNA Element (RE) discovery algorithm (REDa) was applied to a spectrum of tissues and cells representing germline, embryonic, and somatic tissues and examined as a function of differentiation through the first set of cell divisions of human development. This highlighted extensive transcription throughout the genome, yielding previously unidentified human spermatogenic RNAs. Both exonic and novel X-chromosome REs were subject to robust meiotic sex chromosome inactivation, although an extensive de-repression occurred in the post-meiotic stages of spermatogenesis. Surprisingly, 2.4% of the 10,395 X chromosome exonic REs were present in mature sperm. Transcribed genomic repetitive sequences, including simple centromeric repeats, HERVE and HSAT1, were also shown to be associated with RE expression during spermatogenesis. These results suggest that pervasive intergenic repetitive sequence expression during human spermatogenesis may play a role in regulating chromatin dynamics. Repetitive REs switching repeat classes during differentiation upon fertilization and embryonic genome activation was evident.

Transcriptionally induced enhancers in the macrophage immune response to Mycobacterium tuberculosis infection

Background

Tuberculosis is a life-threatening infectious disease caused by Mycobacterium tuberculosis (M.tb). M.tb subverts host immune responses to build a favourable niche and survive inside of host macrophages. Macrophages can control or eliminate the infection, if acquire appropriate functional phenotypes. Transcriptional regulation is a key process that governs the activation and maintenance of these phenotypes. Among the factors orchestrating transcriptional regulation during M.tb infection, transcriptional enhancers still remain unexplored.

Results

We analysed transcribed enhancers in M.tb-infected mouse bone marrow-derived macrophages. We established a link between known M.tb-responsive transcription factors and transcriptional activation of enhancers and their target genes. Our data suggest that enhancers might drive macrophage response via transcriptional activation of key immune genes, such as Tnf, Tnfrsf1b, Irg1, Hilpda, Ccl3, and Ccl4. We report enhancers acquiring transcription de novo upon infection. Finally, we link highly transcriptionally induced enhancers to activation of genes with previously unappreciated roles in M.tb infection, such as Fbxl3, Tapt1, Edn1, and Hivep1.

Conclusions

Our findings suggest the importance of macrophage host transcriptional enhancers during M.tb infection. Our study extends current knowledge of the regulation of macrophage responses to M.tb infection and provides a basis for future functional studies on enhancer-gene interactions in this process.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5450-6) contains supplementary material, which is available to authorized users.

Enhancer RNAs: a missing regulatory layer in gene transcription

Enhancers and super-enhancers exert indispensable roles in maintaining cell identity through spatiotemporally regulating gene transcription. Meanwhile, active enhancers and super-enhancers also produce transcripts termed enhancer RNAs (eRNAs) from their DNA elements. Although enhancers have been identified for more than 30 years, widespread transcription from enhancers are just discovered by genome-wide sequencing and considered as the key to understand longstanding questions in gene transcription. RNA-transcribed enhancers are marked by histone modifications such as H3K4m1/2 and H3K27Ac, and enriched with transcription regulatory factors such as LDTFs, P300, CBP, BRD4 and MED1. Those regulatory factors might constitute a Mega-Trans-like complex to potently activate enhancers. Compared to mRNAs, eRNAs are quite unstable and play roles at local. Functionally, it has been shown that eRNAs promote formation of enhancer-promoter loops. Several studies also demonstrated that eRNAs help the binding of RNA polymerase II (RNAPII) or transition of paused RNAPII by de-association of the negative elongation factor (NELF) complex. Nevertheless, these proposed mechanisms are not universally accepted and still under controversy. Here, we comprehensively summarize the reported findings and make perspectives for future exploration. We also believe that super-enhancer derived RNAs (seRNAs) might be informative to understand the nature of super-enhancers.

Genome-Wide Maps of Transcription Regulatory Elements and Transcription Enhancers in Development and Disease

Gene expression is regulated by numerous elements including enhancers, insulators, transcription factors, and architectural proteins. Regions of DNA distal to the transcriptional start site, called enhancers, play a central role in the temporal and tissue-specific regulation of gene expression through RNA polymerase II. The identification of enhancers and other cis regulatory elements has largely been possible due to advances in next generation sequencing technologies. Enhancers regulate gene expression through chromatin loops mediated by architectural proteins such as YY1, CTCF, the cohesin complex, and LDB1. Additionally, enhancers can be transcribed to produce noncoding RNAs termed enhancer RNAs that likely participate in transcriptional regulation. The central role of enhancers in regulating gene expression implicates them in both normal physiology but also many disease states. The importance of enhancers is evident by the suggested role of SNPs, duplications, and other alterations of enhancer function in many diseases, ranging from cancer to atherosclerosis to chronic kidney disease. Although much progress has been made in recent years, the field of enhancer biology and our knowledge of the cis regulome remains a work in progress. This review will highlight recent seminal studies which demonstrate the role of enhancers in normal physiology and disease pathogenesis. © 2019 American Physiological Society. Compr Physiol 9:439-455, 2019.

Enhancers active in dopamine neurons are a primary link between genetic variation and neuropsychiatric disease

Enhancers function as DNA logic gates and may control specialized functions of billions of neurons. Here we show a tailored program of noncoding genome elements active in situ in physiologically distinct dopamine neurons of the human brain. We found 71,022 transcribed noncoding elements, many of which were consistent with active enhancers and with regulatory mechanisms in zebrafish and mouse brains. Genetic variants associated with schizophrenia, addiction, and Parkinson's disease were enriched in these elements. Expression quantitative trait locus analysis revealed that Parkinson's disease-associated variants on chromosome 17q21 cis-regulate the expression of an enhancer RNA in dopamine neurons. This study shows that enhancers in dopamine neurons link genetic variation to neuropsychiatric traits.

Integrated omics approaches to characterize a nuclear receptor corepressor-associated histone deacetylase in mouse skeletal muscle

Nuclear receptors regulate gene expression by differentially binding to coactivators or corepressors in a ligand-dependent manner, which further recruits a set of epigenome-modifying enzymes that remodel chromatin conformation. Histone acetylation is a major epigenomic change controlled by histone acetyltransferases (HATs) and histone deacetylases (HDACs). HDAC3 is the only HDAC that confers the enzymatic activity to the complexes nucleated by nuclear receptor corepressors NCoR and SMRT. To address the metabolic function of HDAC3, we have deleted it specifically in mouse skeletal muscles. We have performed the following omics profiling in skeletal muscles of these mice: (1) RNA-seq profiling of total RNA; (2) Global nuclear run-on (GRO-seq) analysis of nascent RNAs; (3) Chromatin immuno-precipitation (ChIP-seq) of HDAC3 at both early evening and early morning; (4) proteomics profiling with mass spectrometry; (5) snap-shot metabolomics profiling of water-soluble metabolites at the basal condition; (6) snap-shot metabolomics profiling of lipid species at the basal condition; (7) kinetic fluxomics analysis of glucose utilization using ¹³C₆-glucose In vivo during treadmill running exercise. These approaches have provided several novel insights into how nuclear receptors regulate circadian rhythm of skeletal muscle fuel metabolism, which has been published elsewhere. Here we present the original datasets and technical details during the execution, analysis, and interpretation of these omics studies.

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.