Noncoding transcription at enhancers: general principles and functional models

Intragenic Enhancers Attenuate Host Gene Expression

Eukaryotic gene transcription is regulated at many steps, including RNA polymerase II (Pol II) recruitment, transcription initiation, promoter-proximal Pol II pause release, and transcription termination; however, mechanisms regulating transcription during productive elongation remain poorly understood. Enhancers, which activate gene transcription, themselves undergo Pol II-mediated transcription, but our understanding of enhancer transcription and enhancer RNAs (eRNAs) remains incomplete. Here we show that transcription at intragenic enhancers interferes with and attenuates host gene transcription during productive elongation. While the extent of attenuation correlates positively with nascent eRNA expression, the act of intragenic enhancer transcription alone, but not eRNAs, explains the attenuation. Through CRISPR/Cas9-mediated deletions, we demonstrate a physiological role for intragenic enhancer-mediated transcription attenuation in cell fate determination. We propose that intragenic enhancers not only enhance transcription of one or more genes from a distance but also fine-tune transcription of their host gene through transcription interference, facilitating differential utilization of the same regulatory element for disparate functions.

Sasquatch: predicting the impact of regulatory SNPs on transcription factor binding from cell- and tissue-specific DNase footprints

In the era of genome-wide association studies (GWAS) and personalized medicine, predicting the impact of single nucleotide polymorphisms (SNPs) in regulatory elements is an important goal. Current approaches to determine the potential of regulatory SNPs depend on inadequate knowledge of cell-specific DNA binding motifs. Here, we present Sasquatch, a new computational approach that uses DNase footprint data to estimate and visualize the effects of noncoding variants on transcription factor binding. Sasquatch performs a comprehensive k-mer-based analysis of DNase footprints to determine any k-mer's potential for protein binding in a specific cell type and how this may be changed by sequence variants. Therefore, Sasquatch uses an unbiased approach, independent of known transcription factor binding sites and motifs. Sasquatch only requires a single DNase-seq data set per cell type, from any genotype, and produces consistent predictions from data generated by different experimental procedures and at different sequence depths. Here we demonstrate the effectiveness of Sasquatch using previously validated functional SNPs and benchmark its performance against existing approaches. Sasquatch is available as a versatile webtool incorporating publicly available data, including the human ENCODE collection. Thus, Sasquatch provides a powerful tool and repository for prioritizing likely regulatory SNPs in the noncoding genome.

Mutations in ACTRT1 and its enhancer RNA elements lead to aberrant activation of Hedgehog signaling in inherited and sporadic basal cell carcinomas

Basal cell carcinoma (BCC), the most common human cancer, results from aberrant activation of the Hedgehog signaling pathway. Although most cases of BCC are sporadic, some forms are inherited, such as Bazex-Dupré-Christol syndrome (BDCS)-a cancer-prone genodermatosis with an X-linked, dominant inheritance pattern. We have identified mutations in the ACTRT1 gene, which encodes actin-related protein T1 (ARP-T1), in two of the six families with BDCS that were examined in this study. High-throughput sequencing in the four remaining families identified germline mutations in noncoding sequences surrounding ACTRT1. These mutations were located in transcribed sequences encoding enhancer RNAs (eRNAs) and were shown to impair enhancer activity and ACTRT1 expression. ARP-T1 was found to directly bind to the GLI1 promoter, thus inhibiting GLI1 expression, and loss of ARP-T1 led to activation of the Hedgehog pathway in individuals with BDCS. Moreover, exogenous expression of ACTRT1 reduced the in vitro and in vivo proliferation rates of cell lines with aberrant activation of the Hedgehog signaling pathway. In summary, our study identifies a disease mechanism in BCC involving mutations in regulatory noncoding elements and uncovers the tumor-suppressor properties of ACTRT1.

Gene regulation in the immune system by long noncoding RNAs

Long noncoding RNAs (lncRNAs) are emerging as critical regulators of gene expression in the immune system. Studies have shown that lncRNAs are expressed in a highly lineage-specific manner and control the differentiation and function of innate and adaptive cell types. In this Review, we focus on mechanisms used by lncRNAs to regulate genes encoding products involved in the immune response, including direct interactions with chromatin, RNA and proteins. In addition, we address new areas of lncRNA biology, such as the functions of enhancer RNAs, circular RNAs and chemical modifications to RNA in cellular processes. We emphasize critical gaps in knowledge and future prospects for the roles of lncRNAs in the immune system and autoimmune disease.

Single-cell profiling reveals that eRNA accumulation at enhancer-promoter loops is not required to sustain transcription

Enhancers are intergenic DNA elements that regulate the transcription of target genes in response to signaling pathways by interacting with promoters over large genomic distances. Recent studies have revealed that enhancers are bi-directionally transcribed into enhancer RNAs (eRNAs). Using single-molecule fluorescence in situ hybridization (smFISH), we investigated the eRNA-mediated regulation of transcription during estrogen induction in MCF-7 cells. We demonstrate that eRNAs are localized exclusively in the nucleus and are induced with similar kinetics as target mRNAs. However, eRNAs are mostly nascent at enhancers and their steady-state levels remain lower than those of their cognate mRNAs. Surprisingly, at the single-allele level, eRNAs are rarely co-expressed with their target loci, demonstrating that active gene transcription does not require the continuous transcription of eRNAs or their accumulation at enhancers. When co-expressed, sub-diffraction distance measurements between nascent mRNA and eRNA signals reveal that co-transcription of eRNAs and mRNAs rarely occurs within closed enhancer–promoter loops. Lastly, basal eRNA transcription at enhancers, but not E2-induced transcription, is maintained upon depletion of MLL1 and ERα, suggesting some degree of chromatin accessibility prior to signal-dependent activation of transcription. Together, our findings suggest that eRNA accumulation at enhancer–promoter loops is not required to sustain target gene transcription.

Profiles of Long Noncoding RNAs in Human Naive and Memory T Cells

We employed whole-genome RNA-sequencing to profile mRNAs and both annotated and novel long noncoding RNAs (lncRNAs) in human naive, central memory, and effector memory CD4+ T cells. Loci transcribing both lineage-specific annotated and novel lncRNA are adjacent to lineage-specific protein-coding genes in the genome. Lineage-specific novel lncRNA loci are transcribed from lineage-specific typical- and supertranscriptional enhancers and are not multiexonic, thus are more similar to enhancer RNAs. Novel enhancer-associated lncRNAs transcribed from the IFNG locus bind the transcription factor NF-κB and enhance binding of NF-κB to the IFNG genomic locus. Depletion of the annotated lncRNA, IFNG-AS1, or one IFNG enhancer-associated lncRNA abrogates IFNG expression by memory T cells, indicating these lncRNAs have biologic function.

Non-coding Transcripts from Enhancers: New Insights into Enhancer Activity and Gene Expression Regulation

Long non-coding RNAs (lncRNAs) have gained widespread interest in the past decade owing to their enormous amount and surprising functions implicated in a variety of biological processes. Some lncRNAs exert function as enhancers, i.e., activating gene transcription by serving as the cis-regulatory molecules. Furthermore, recent studies have demonstrated that many enhancer elements can be transcribed and produce RNA molecules, which are termed as enhancer RNAs (eRNAs). The eRNAs are not merely the by-product of the enhancer transcription. In fact, many of them directly exert or regulate enhancer activity in gene activation through diverse mechanisms. Here, we provide an overview of enhancer activity, transcription of enhancer itself, characteristics of eRNAs, as well as their roles in regulating enhancer activity and gene expression.

Enhancer-derived RNA: A Primer

Enhancer-derived RNAs (eRNAs) are a group of RNAs transcribed by RNA polymerase II from the domain of transcription enhancers, a major type of cis-regulatory elements in the genome. The correlation between eRNA production and enhancer activity has stimulated studies on the potential role of eRNAs in transcriptional regulation. Additionally, eRNA has also served as a marker for global identification of enhancers. Here I review the brief history and fascinating properties of eRNAs.

Expression of long non-coding RNAs in autoimmunity and linkage to enhancer function and autoimmune disease risk genetic variants

Genome-wide association studies have identified numerous genetic variants conferring autoimmune disease risk. Most of these genetic variants lie outside protein-coding genes hampering mechanistic explorations. Numerous mRNAs are also differentially expressed in autoimmune disease but their regulation is also unclear. The majority of the human genome is transcribed yet its biologic significance is incompletely understood. We performed whole genome RNA-sequencing [RNA-seq] to categorize expression of mRNAs, known and novel long non-coding RNAs [lncRNAs] in leukocytes from subjects with autoimmune disease and identified annotated and novel lncRNAs differentially expressed across multiple disorders. We found that loci transcribing novel lncRNAs were not randomly distributed across the genome but co-localized with leukocyte transcriptional enhancers, especially super-enhancers, and near genetic variants associated with autoimmune disease risk. We propose that alterations in enhancer function, including lncRNA expression, produced by genetics and environment, change cellular phenotypes contributing to disease risk and pathogenesis and represent attractive therapeutic targets.

Landscape of histone modifications in a sponge reveals the origin of animal cis-regulatory complexity

Combinatorial patterns of histone modifications regulate developmental and cell type-specific gene expression and underpin animal complexity, but it is unclear when this regulatory system evolved. By analysing histone modifications in a morphologically-simple, early branching animal, the sponge Amphimedonqueenslandica, we show that the regulatory landscape used by complex bilaterians was already in place at the dawn of animal multicellularity. This includes distal enhancers, repressive chromatin and transcriptional units marked by H3K4me3 that vary with levels of developmental regulation. Strikingly, Amphimedon enhancers are enriched in metazoan-specific microsyntenic units, suggesting that their genomic location is extremely ancient and likely to place constraints on the evolution of surrounding genes. These results suggest that the regulatory foundation for spatiotemporal gene expression evolved prior to the divergence of sponges and eumetazoans, and was necessary for the evolution of animal multicellularity.

Can We Predict Gene Expression by Understanding Proximal Promoter Architecture?

We review computational predictions of expression from the promoter architecture - the set of transcription factors that can bind the proximal promoter. We focus on spatial expression patterns in animals with complex body plans and many distinct tissue types. This field is ripe for change as functional genomics datasets accumulate for both expression and protein-DNA interactions. While there has been some success in predicting the breadth of expression (i.e., the fraction of tissue types a gene is expressed in), predicting tissue specificity remains challenging. We discuss how progress can be achieved through either machine learning or complementary combinatorial data mining. The likely impact of single-cell expression data is considered. Finally, we discuss the design of artificial promoters as a practical application.

Re-engineering the Pancreas Tumor Microenvironment: A "Regenerative Program" Hacked

The "hallmarks" of pancreatic ductal adenocarcinoma (PDAC) include proliferative, invasive, and metastatic tumor cells and an associated dense desmoplasia comprised of fibroblasts, pancreatic stellate cells, extracellular matrix, and immune cells. The oncogenically activated pancreatic epithelium and its associated stroma are obligatorily interdependent, with the resulting inflammatory and immunosuppressive microenvironment contributing greatly to the evolution and maintenance of PDAC. The peculiar pancreas-specific tumor phenotype is a consequence of oncogenes hacking the resident pancreas regenerative program, a tissue-specific repair mechanism regulated by discrete super enhancer networks. Defined as genomic regions containing clusters of multiple enhancers, super enhancers play pivotal roles in cell/tissue specification, identity, and maintenance. Hence, interfering with such super enhancer-driven repair networks should exert a disproportionately disruptive effect on tumor versus normal pancreatic tissue. Novel drugs that directly or indirectly inhibit processes regulating epigenetic status and integrity, including those driven by histone deacetylases, histone methyltransferase and hydroxylases, DNA methyltransferases, various metabolic enzymes, and bromodomain and extraterminal motif proteins, have shown the feasibility of disrupting super enhancer-dependent transcription in treating multiple tumor types, including PDAC. The idea that pancreatic adenocarcinomas rely on embedded super enhancer transcriptional mechanisms suggests a vulnerability that can be potentially targeted as novel therapies for this intractable disease. Clin Cancer Res; 23(7); 1647-55. ©2017 AACRSee all articles in this CCR Focus section, "Pancreatic Cancer: Challenge and Inspiration."

Comprehensive analysis of long non-coding RNAs highlights their spatio-temporal expression patterns and evolutional conservation in Sus scrofa

Despite modest sequence conservation and rapid evolution, long non-coding RNAs (lncRNAs) appear to be conserved in expression pattern and function. However, analysis of lncRNAs across tissues and developmental stages remains largely uncharacterized in mammals. Here, we systematically investigated the lncRNAs of the Guizhou miniature pig (Sus scrofa), which was widely used as biomedical model. We performed RNA sequencing across 9 organs and 3 developmental skeletal muscle, and developed a filtering pipeline to identify 10,813 lncRNAs (9,075 novel). Conservation patterns analysis revealed that 57% of pig lncRNAs showed homology to humans and mice based on genome alignment. 5,455 lncRNAs exhibited typical hallmarks of regulatory molecules, such as high spatio-temporal specificity. Notably, conserved lncRNAs exhibited higher tissue specificity than pig-specific lncRNAs and were significantly enriched in testis and ovary. Weighted co-expression network analysis revealed a set of conserved lncRNAs that are likely involved in postnatal muscle development. Based on the high degree of similarity in the structure, organization, and dynamic expression of pig lncRNAs compared with human and mouse lncRNAs, we propose that these lncRNAs play an important role in organ physiology and development in mammals. Our results provide a resource for studying animal evolution, morphological complexity, breeding, and biomedical research.

Retrieving Chromatin Patterns from Deep Sequencing Data Using Correlation Functions

Epigenetic modifications and other chromatin features partition the genome on multiple length scales. They define chromatin domains with distinct biological functions that come in sizes ranging from single modified DNA bases to several megabases in the case of heterochromatic histone modifications. Due to chromatin folding, domains that are well separated along the linear nucleosome chain can form long-range interactions in three-dimensional space. It has now become a routine task to map epigenetic marks and chromatin structure by deep sequencing methods. However, assessing and comparing the properties of chromatin domains and their positional relationships across data sets without a priori assumptions remains challenging. Here, we introduce multiscale correlation evaluation (MCORE), which uses the fluctuation spectrum of mapped sequencing reads to quantify and compare chromatin patterns over a broad range of length scales in a model-independent manner. We applied MCORE to map the chromatin landscape in mouse embryonic stem cells and differentiated neural cells. We integrated sequencing data from chromatin immunoprecipitation, RNA expression, DNA methylation, and chromosome conformation capture experiments into network models that reflect the positional relationships among these features on different genomic scales. Furthermore, we used MCORE to compare our experimental data to models for heterochromatin reorganization during differentiation. The application of correlation functions to deep sequencing data complements current evaluation schemes and will support the development of quantitative descriptions of chromatin networks.

The origins of developmental gene regulation

The leap from simple unicellularity to complex multicellularity remains one of life's major enigmas. The origins of metazoan developmental gene regulatory mechanisms are sought by analyzing gene regulation in extant eumetazoans, sponges, and unicellular organisms. The main hypothesis of this manuscript is that, developmental enhancers evolved from unicellular inducible promoters that diversified the expression of regulatory genes during metazoan evolution. Promoters and enhancers are functionally similar; both can regulate the transcription of distal promoters and both direct local transcription. Additionally, enhancers have experimentally characterized structural features that reveal their origin from inducible promoters. The distal co-operative regulation among promoters identified in unicellular opisthokonts possibly represents the precursor of distal regulation of promoters by enhancers. During metazoan evolution, constitutive-type promoters of regulatory genes would have acquired novel receptivity to distal regulatory inputs from promoters of inducible genes that eventually specialized as enhancers. The novel regulatory interactions would have caused constitutively expressed genes controlling differential gene expression in unicellular organisms to become themselves differentially expressed. The consequence of the novel regulatory interactions was that regulatory pathways of unicellular organisms became interlaced and ultimately evolved into the intricate developmental gene regulatory networks (GRNs) of extant metazoans.

Inferring a role for methylation of intergenic DNA in the regulation of genes aberrantly expressed in precursor B-cell acute lymphoblastic leukemia

A complete understanding of the mechanisms involved in the development of pre-B ALL is lacking. In this study, we integrated DNA methylation data and gene expression data to elucidate the impact of aberrant intergenic DNA methylation on gene expression in pre-B ALL. We found a subset of differentially methylated intergenic loci that were associated with altered gene expression in pre-B ALL patients. Notably, 84% of these regions were also bound by transcription factors (TF) known to play roles in differentiation and B-cell development in a lymphoblastoid cell line. Further, an overall downregulation of eRNA transcripts was observed in pre-B ALL patients and these transcripts were associated with the downregulation of putative target genes involved in B-cell migration, proliferation, and apoptosis. The identification of novel putative regulatory regions highlights the significance of intergenic DNA sequences and may contribute to the identification of new therapeutic targets for the treatment of pre-B ALL.

Long Noncoding RNA: Genome Organization and Mechanism of Action

For the last four decades, we have known that noncoding RNAs maintain critical housekeeping functions such as transcription, RNA processing, and translation. However, in the late 1990s and early 2000s, the advent of high-throughput sequencing technologies and computational tools to analyze these large sequencing datasets facilitated the discovery of thousands of small and long noncoding RNAs (lncRNAs) and their functional role in diverse biological functions. For example, lncRNAs have been shown to regulate dosage compensation, genomic imprinting, pluripotency, cell differentiation and development, immune response, etc. Here we review how lncRNAs bring about such copious functions by employing diverse mechanisms such as translational inhibition, mRNA degradation, RNA decoys, facilitating recruitment of chromatin modifiers, regulation of protein activity, regulating the availability of miRNAs by sponging mechanism, etc. In addition, we provide a detailed account of different mechanisms as well as general principles by which lncRNAs organize functionally different nuclear sub-compartments and their impact on nuclear architecture.

Biogenesis and Transcriptional Regulation of Long Noncoding RNAs in the Human Immune System

The central dogma of molecular biology states that DNA makes RNA makes protein. Discoveries over the last quarter of a century found that the process of DNA transcription into RNA gives rise to a diverse array of functional RNA species, including genes that code for protein and noncoding RNAs. For decades, the focus has been on understanding how protein-coding genes are regulated to influence protein expression. However, with the completion of the Human Genome Project and follow-up ENCODE data, it is now appreciated that only 2-3% of the genome codes for protein-coding gene exons and that the bulk of the transcribed genome, apart from ribosomal RNAs, is at the level of noncoding RNA genes. In this article, we focus on the biogenesis and regulation of a distinct class of noncoding RNA molecules termed long, noncoding RNAs in the context of the immune system.

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

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

MIST1 and PTF1 Collaborate in Feed-Forward Regulatory Loops That Maintain the Pancreatic Acinar Phenotype in Adult Mice

Much remains unknown regarding the regulatory networks formed by transcription factors in mature, differentiated mammalian cells in vivo, despite many studies of individual DNA-binding transcription factors. We report a constellation of feed-forward loops formed by the pancreatic transcription factors MIST1 and PTF1 that govern the differentiated phenotype of the adult pancreatic acinar cell. PTF1 is an atypical basic helix-loop-helix transcription factor complex of pancreatic acinar cells and is critical to acinar cell fate specification and differentiation. MIST1, also a basic helix-loop-helix transcription factor, enhances the formation and maintenance of the specialized phenotype of professional secretory cells. The MIST1 and PTF1 collaboration controls a wide range of specialized cellular processes, including secretory protein synthesis and processing, exocytosis, and homeostasis of the endoplasmic reticulum. PTF1 drives Mist1 transcription, and MIST1 and PTF1 bind and drive the transcription of over 100 downstream acinar genes. PTF1 binds two canonical bipartite sites within a 0.7-kb transcriptional enhancer upstream of Mist1 that are essential for the activity of the enhancer in vivo MIST1 and PTF1 coregulate target genes synergistically or additively, depending on the target transcriptional enhancer. The frequent close binding proximity of PTF1 and MIST1 in pancreatic acinar cell chromatin implies extensive collaboration although the collaboration is not dependent on a stable physical interaction.

Transcriptional enhancers: Transcription, function and flexibility

Active transcriptional enhancers are often transcribed to eRNAs, whose changing levels mirror those of the target gene mRNA. We discuss some of the reported functions of these eRNAs and their likely diversity to allow utilization of distinct cis regulatory regions to enhance transcription in diverse developmental and cellular contexts.

Epigenetic Regulation of Monocyte and Macrophage Function

SIGNIFICANCE

Monocytes and macrophages are key players in tissue homeostasis and immune responses. Epigenetic processes tightly regulate cellular functioning in health and disease. Recent Advances: Recent technical developments have allowed detailed characterizations of the transcriptional circuitry underlying monocyte and macrophage regulation. Upon differentiation and activation, enhancers are selected by lineage-determining and signal-dependent transcription factors. Enhancers are shown to be very dynamic and activation of these enhancers underlies the differences in gene transcription between monocytes and macrophages and their subtypes.

CRITICAL ISSUES

It has been shown that epigenetic enzymes regulate the functioning of these cells and targeting of epigenetic enzymes has been proven to be a valuable tool to dampen inflammatory responses. We give a comprehensive overview of recent developments and understanding of the epigenetic pathways that control monocyte and macrophage function and of the epigenetic enzymes involved in monocyte and macrophage differentiation and activation.

FUTURE DIRECTIONS

The key challenges in the upcoming years will be to study epigenetic changes in human disease and to better understand how epigenetic pathways control the inflammatory repertoire in disease. Antioxid. Redox Signal. 25, 758-774.

A generative model for the behavior of RNA polymerase

MOTIVATION

Transcription by RNA polymerase is a highly dynamic process involving multiple distinct points of regulation. Nascent transcription assays are a relatively new set of high throughput techniques that measure the location of actively engaged RNA polymerase genome wide. Hence, nascent transcription is a rich source of information on the regulation of RNA polymerase activity. To fully dissect this data requires the development of stochastic models that can both deconvolve the stages of polymerase activity and identify significant changes in activity between experiments.

RESULTS

We present a generative, probabilistic model of RNA polymerase that fully describes loading, initiation, elongation and termination. We fit this model genome wide and profile the enzymatic activity of RNA polymerase across various loci and following experimental perturbation. We observe striking correlation of predicted loading events and regulatory chromatin marks. We provide principled statistics that compute probabilities reminiscent of traveler's and divergent ratios. We finish with a systematic comparison of RNA Polymerase activity at promoter versus non-promoter associated loci.

AVAILABILITY AND IMPLEMENTATION

Transcription Fit (Tfit) is a freely available, open source software package written in C/C ++ that requires GNU compilers 4.7.3 or greater. Tfit is available from GitHub (https://github.com/azofeifa/Tfit).

CONTACT

robin.dowell@colorado.eduSupplementary information: Supplementary data are available at Bioinformatics online.

Chromatin remodeling effects on enhancer activity

During organism development, a diversity of cell types emerges with disparate, yet stable profiles of gene expression with distinctive cellular functions. In addition to gene promoters, the genome contains enhancer regulatory sequences, which are implicated in cellular specialization by facilitating cell-type and tissue-specific gene expression. Enhancers are DNA binding elements characterized by highly sophisticated and various mechanisms of action allowing for the specific interaction of general and tissue-specific transcription factors (TFs). However, eukaryotic organisms package their genetic material into chromatin, generating a physical barrier for TFs to interact with their cognate sequences. The ability of TFs to bind DNA regulatory elements is also modulated by changes in the chromatin structure, including histone modifications, histone variants, ATP-dependent chromatin remodeling, and the methylation status of DNA. Furthermore, it has recently been revealed that enhancer sequences are also transcribed into a set of enhancer RNAs with regulatory potential. These interdependent processes act in the context of a complex network of chromatin interactions, which together contributes to a renewed vision of how gene activation is coordinated in a cell-type-dependent manner. In this review, we describe the interplay between genetic and epigenetic aspects associated with enhancers and discuss their possible roles on enhancer function.

Enhancers reside in a unique epigenetic environment during early zebrafish development

BACKGROUND

Enhancers, not promoters, are the most dynamic in their DNA methylation status throughout development and differentiation. Generally speaking, enhancers that are primed to or actually drive gene expression are characterized by relatively low levels of DNA methylation (hypo-methylation), while inactive enhancers display hyper-methylation of the underlying DNA. The direct functional significance of the DNA methylation state of enhancers is, however, unclear for most loci.

RESULTS

In contrast to conventional epigenetic interactions at enhancers, we find that DNA methylation status and enhancer activity during early zebrafish development display very unusual correlation characteristics: hypo-methylation is a unique feature of primed enhancers whereas active enhancers are generally hyper-methylated. The hypo-methylated enhancers that we identify (hypo-enhancers) are enriched close to important transcription factors that act later in development. Interestingly, hypo-enhancers are de-methylated shortly before the midblastula transition and reside in a unique epigenetic environment. Finally, we demonstrate that hypo-enhancers do become active at later developmental stages and that they are physically associated with the transcriptional start site of target genes, irrespective of target gene activity.

CONCLUSIONS

We demonstrate that early development in zebrafish embodies a time window characterized by non-canonical DNA methylation-enhancer relationships, including global DNA hypo-methylation of inactive enhancers and DNA hyper-methylation of active enhancers.

Hierarchy within the mammary STAT5-driven Wap super-enhancer

Super-enhancers comprise of dense transcription factor platforms highly enriched for active chromatin marks. A paucity of functional data led us to investigate their role in the mammary gland, an organ characterized by exceptional gene regulatory dynamics during pregnancy. ChIP-Seq for the master regulator STAT5, the glucocorticoid receptor, H3K27ac and MED1, identified 440 mammary-specific super-enhancers, half of which were associated with genes activated during pregnancy. We interrogated the Wap super-enhancer, generating mice carrying mutations in STAT5 binding sites within its three constituent enhancers. Individually, only the most distal site displayed significant enhancer activity. However, combinatorial mutations showed that the 1,000-fold gene induction relied on all enhancers. Disabling the binding sites of STAT5, NFIB and ELF5 in the proximal enhancer incapacitated the entire super-enhancer, suggesting an enhancer hierarchy. The identification of mammary-specific super-enhancers and the mechanistic exploration of the Wap locus provide insight into the complexity of cell-specific and hormone-regulated genes.

Mechanisms of long noncoding RNA function in development and disease

Since decades it has been known that non-protein-coding RNAs have important cellular functions. Deep sequencing recently facilitated the discovery of thousands of novel transcripts, now classified as long noncoding RNAs (lncRNAs), in many vertebrate and invertebrate species. LncRNAs are involved in a wide range of cellular mechanisms, from almost all aspects of gene expression to protein translation and stability. Recent findings implicate lncRNAs as key players of cellular differentiation, cell lineage choice, organogenesis and tissue homeostasis. Moreover, lncRNAs are involved in pathological conditions such as cancer and cardiovascular disease, and therefore provide novel biomarkers and pharmaceutical targets. Here we discuss examples illustrating the versatility of lncRNAs in gene control, development and differentiation, as well as in human disease.

Long non-coding RNAs in cancer drug resistance development

The presence or emergence of chemoresistance in tumor cells is a major burden in cancer therapy. While drug resistance is a multifactorial phenomenon arising from altered membrane transport of drugs, altered drug metabolism, altered DNA repair, reduced apoptosis rate and alterations of drug metabolism, it can also be linked to genetic and epigenetic factors. Long non-coding RNAs (lncRNAs) have important regulatory roles in many aspects of genome function including gene transcription, splicing, and epigenetics as well as biological processes involved in cell cycle, cell differentiation, development, and pluripotency. As such, it may not be surprising that some lncRNAs have been recently linked to carcinogenesis and drug resistance/sensitivity. Research is accelerating to decipher the exact molecular mechanism of lncRNA-regulated drug resistance and its therapeutic implications. In this article, we will review the structure, biogenesis, and mode of action of lncRNAs. Then, the involvement of lncRNAs in drug resistance will be discussed in detail.

PEDLA: predicting enhancers with a deep learning-based algorithmic framework

Transcriptional enhancers are non-coding segments of DNA that play a central role in the spatiotemporal regulation of gene expression programs. However, systematically and precisely predicting enhancers remain a major challenge. Although existing methods have achieved some success in enhancer prediction, they still suffer from many issues. We developed a deep learning-based algorithmic framework named PEDLA (https://github.com/wenjiegroup/PEDLA), which can directly learn an enhancer predictor from massively heterogeneous data and generalize in ways that are mostly consistent across various cell types/tissues. We first trained PEDLA with 1,114-dimensional heterogeneous features in H1 cells, and demonstrated that PEDLA framework integrates diverse heterogeneous features and gives state-of-the-art performance relative to five existing methods for enhancer prediction. We further extended PEDLA to iteratively learn from 22 training cell types/tissues. Our results showed that PEDLA manifested superior performance consistency in both training and independent test sets. On average, PEDLA achieved 95.0% accuracy and a 96.8% geometric mean (GM) of sensitivity and specificity across 22 training cell types/tissues, as well as 95.7% accuracy and a 96.8% GM across 20 independent test cell types/tissues. Together, our work illustrates the power of harnessing state-of-the-art deep learning techniques to consistently identify regulatory elements at a genome-wide scale from massively heterogeneous data across diverse cell types/tissues.

Noncoding Transcription Is a Driving Force for Nucleosome Instability in spt16 Mutant Cells

FACT (facilitates chromatin transcription) consists of two essential subunits, Spt16 and Pob3, and functions as a histone chaperone. Mutation of spt16 results in a global loss of nucleosomes as well as aberrant transcription. Here, we show that the majority of nucleosome changes upon Spt16 depletion are alterations in nucleosome fuzziness and position shift. Most nucleosomal changes are suppressed by the inhibition of RNA polymerase II (Pol II) activity. Surprisingly, a small subgroup of nucleosome changes is resistant to transcriptional inhibition. Notably, Spt16 and distinct histone modifications are enriched at this subgroup of nucleosomes. We also report 1,037 Spt16-suppressed noncoding transcripts (SNTs) and found that the SNT start sites are enriched with the subgroup of nucleosomes resistant to Pol II inhibition. Finally, the nucleosomes at genes overlapping SNTs are more susceptible to changes upon Spt16 depletion than those without SNTs. Taken together, our results support a model in which Spt16 has a role in maintaining local nucleosome stability to inhibit initiation of SNT transcription, which once initiated drives additional nucleosome loss upon Spt16 depletion.

T-bet Activates Th1 Genes through Mediator and the Super Elongation Complex

The transcription factor T-bet directs Th1 cell differentiation, but the molecular mechanisms that underlie this lineage-specific gene regulation are not completely understood. Here, we show that T-bet acts through enhancers to allow the recruitment of Mediator and P-TEFb in the form of the super elongation complex (SEC). Th1 genes are occupied by H3K4me3 and RNA polymerase II in Th2 cells, while T-bet-mediated recruitment of P-TEFb in Th1 cells activates transcriptional elongation. P-TEFb is recruited to both genes and enhancers, where it activates enhancer RNA transcription. P-TEFb inhibition and Mediator and SEC knockdown selectively block activation of T-bet target genes, and P-TEFb inhibition abrogates Th1-associated experimental autoimmune uveitis. T-bet activity is independent of changes in NF-κB RelA and Brd4 binding, with T-bet- and NF-κB-mediated pathways instead converging to allow P-TEFb recruitment. These data provide insight into the mechanism through which lineage-specifying factors promote differentiation of alternative T cell fates.

12-HETER1/GPR31, a high-affinity 12(S)-hydroxyeicosatetraenoic acid receptor, is significantly up-regulated in prostate cancer and plays a critical role in prostate cancer progression

Previously we identified and deorphaned G-protein-coupled receptor 31 (GPR31) as the high-affinity 12(S)-hydroxyeicosatetraenoic acid [12(S)-HETE] receptor (12-HETER1). Here we have determined its distribution in prostate cancer tissue and its role in prostate tumorigenesis using in vitro and in vivo assays. Data-mining studies strongly suggest that 12-HETER1 expression positively correlates with the aggressiveness and progression of prostate tumors. This was corroborated with real-time PCR analysis of human prostate tumor tissue arrays that revealed the expression of 12-HETER1 positively correlates with the clinical stages of prostate cancers and Gleason scores. Immunohistochemistry analysis also proved that the expression of 12-HETER1 is positively correlated with the grades of prostate cancer. Knockdown of 12-HETER1 in prostate cancer cells markedly reduced colony formation and inhibited tumor growth in animals. To discover the regulatory factors, 5 candidate 12-HETER1 promoter cis elements were assayed as luciferase reporter fusions in Chinese hamster ovary (CHO) cells, where the putative cis element required for gene regulation was mapped 2 kb upstream of the 12-HETER1 transcriptional start site. The data implicate 12-HETER1 in a critical new role in the regulation of prostate cancer progression and offer a novel alternative target for therapeutic intervention.-Honn, K. V., Guo, Y., Cai, Y., Lee, M.-J., Dyson, G., Zhang, W., Tucker, S. C. 12-HETER1/GPR31, a high-affinity 12(S)-hydroxyeicosatetraenoic acid receptor, is significantly up-regulated in prostate cancer and plays a critical role in prostate cancer progression.

Molecular control of activation and priming in macrophages

In tissues, macrophages are exposed to metabolic, homeostatic and immunoregulatory signals of local or systemic origin that influence their basal functions and responses to danger signals. Signal-transduction pathways regulated by extracellular signals are coupled to distinct sets of broadly expressed stimulus-regulated transcription factors whose ability to elicit gene-expression changes is influenced by the accessibility of their binding sites in the macrophage genome. In turn, accessibility of macrophage-specific transcriptional regulatory elements (enhancers and promoters) is specified by transcription factors that determine the macrophage lineage or impose their tissue-specific properties. Here we review recent findings that advance the understanding of mechanisms underlying priming and signal-dependent activation of macrophages and discuss the effect of genetic variation on these processes.

Induction of Gnrh mRNA expression by the ω-3 polyunsaturated fatty acid docosahexaenoic acid and the saturated fatty acid palmitate in a GnRH-synthesizing neuronal cell model, mHypoA-GnRH/GFP

Gonadotropin-releasing hormone (GnRH) neurons coordinate reproduction. However, whether GnRH neurons directly sense free fatty acids (FFAs) is unknown. We investigated the individual effects of the FFAs docosahexaenoic acid (DHA), palmitate, palmitoleate, and oleate (100 μM each) on Gnrh mRNA expression in the mHypoA-GnRH/GFP neuronal cell model. We report that 2 h exposure to palmitate or DHA increases Gnrh transcription. Using the inhibitors AH7614, K252c, U0126, wortmannin, and LY294002, we demonstrate that the effect of DHA is mediated through GPR120 to downstream PKC/MAPK and PI3K signaling. Our results indicate that the effect of palmitate may depend on palmitoyl-coA synthesis and PI3K signaling. Finally, we demonstrate that both DHA and palmitate increase Gnrh enhancer-derived RNA levels. Overall, these studies provide evidence that GnRH neurons directly sense FFAs. This will advance our understanding of the mechanisms underlying FFA sensing in the brain and provides insight into the links between nutrition and reproductive function.

Long noncoding RNAs in kidney and cardiovascular diseases

Transcription of a large part of the human genome results in RNA transcripts that have limited or no protein-coding potential. These include long noncoding RNAs (lncRNAs), which are defined as being ≥200 nucleotides long. Unlike microRNAs, which have been extensively studied, little is known about the functional role of lncRNAs. However, studies over the past 5 years have shown that lncRNAs interfere with tissue homeostasis and have a role in pathological processes, including in the kidney and heart. The developmental expression of the microRNA sponge H19, for example, is altered in the kidneys of embryos carried by hyperglycaemic mothers, and the lncRNA Malat1 regulates hyperglycaemia-induced inflammation in endothelial cells. Putative roles for other lncRNAs have been identified in conditions such as heart failure, cardiac autophagy, hypertension, acute kidney injury, glomerular diseases, acute allograft rejection and renal cell carcinoma. This Review outlines our current understanding of the role and function of lncRNAs in kidney and cardiovascular disease as novel important regulators and potential therapeutic entry points of disease progression.

Unique features of long non-coding RNA biogenesis and function

Long non-coding RNAs (lncRNAs) are a diverse class of RNAs that engage in numerous biological processes across every branch of life. Although initially discovered as mRNA-like transcripts that do not encode proteins, recent studies have revealed features of lncRNAs that further distinguish them from mRNAs. In this Review, we describe special events in the lifetimes of lncRNAs - before, during and after transcription - and discuss how these events ultimately shape the unique characteristics and functional roles of lncRNAs.

Chromatin interactions and candidate genes at ten prostate cancer risk loci

Genome-wide association studies have identified more than 100 common single nucleotide polymorphisms (SNPs) that are associated with prostate cancer risk. However, the vast majority of these SNPs lie in noncoding regions of the genome. To test whether these risk SNPs regulate their target genes through long-range chromatin interactions, we applied capture-based 3C sequencing technology to investigate possible cis-interactions at ten prostate cancer risk loci in six cell lines. We identified significant physical interactions between risk regions and their potential target genes including CAPG at 2p11.2, C2orf43 at 2p24.1, RFX6 at 6q22.1, NFASC at 1q32.1, MYC at 8q24.1 and AGAP7P at 10q11.23. Most of the interaction peaks were co-localized to regions of active histone modification and transcription factor binding sites. Expression quantitative trait locus (eQTL) analysis showed suggestive eQTL signals at rs1446669, rs699664 and rs1078004 for CAPG (p < 0.004), rs13394027 for C2orf43 (p = 2.25E-27), rs10993994 and rs4631830 for AGAP7P (p < 8.02E-5). Further analysis revealed an enhancer activity at genomic region surrounding rs4631830 which was expected to disrupt HOXB-like DNA binding affinity. This study identifies a set of candidate genes and their potential regulatory variants, and provides additional evidence showing the role of long-range chromatin interactions in prostate cancer etiology.

Inhibition of Interleukin 1β (IL-1β) Expression by Anthrax Lethal Toxin (LeTx) Is Reversed by Histone Deacetylase 8 (HDAC8) Inhibition in Murine Macrophages

Many pathogenic microbes often release toxins that subvert the host's immune responses to render the environment suitable for their survival and proliferation. LeTx is one of the toxins causing immune paralysis by cleaving and inactivating the mitogen-activated protein kinase (MAPK) kinases (MEKs). Here, we show that inhibition of the histone deacetylase 8 (HDAC8) by either the HDAC8-specific inhibitor PCI-34051 or small interference (si)RNAs rendered LeTx-exposed murine macrophages responsive to LPS in pro-IL-1β production. HDAC8 selectively targeted acetylated histone H3 lysine 27 (H3K27Ac), which is known to associate with active enhancers. LeTx induced HDAC8 expression, in part through inhibiting p38 MAPK, which resulted in a decrease of H3K27Ac levels. Inhibition of HDAC8 increased H3K27Ac levels and enhanced NF-κB-mediated pro-IL-1β enhancer and messenger RNA production in LeTx-exposed macrophages. Collectively, this study demonstrates a novel role of HDAC8 in LeTx immunotoxicity and regulation of pro-IL-1β production likely through eRNAs. Targeting HDAC8 could be a strategy for enhancing immune responses in macrophages exposed to LeTx or other toxins that inhibit MAPKs.

Noncoding Regulatory RNAs in Hematopoiesis

Hematopoiesis is a dynamic process in which blood cells are continuously generated from hematopoietic stem cells (HSCs). The regulatory mechanisms controlling HSC fate have been studied extensively over the past several decades. Although many protein-coding genes have been shown to regulate hematopoietic differentiation, additional levels of HSC regulation are not well studied. Advances in deep sequencing have revealed many new classes of regulatory noncoding RNAs (ncRNAs), such as enhancer RNAs and antisense ncRNAs. Functional analysis of some of these ncRNAs has provided insights into the molecular mechanisms that regulate hematopoietic development and disease. In this review, we summarize recent advances in our understanding of functional regulatory ncRNAs associated with hematopoietic self-renewal and differentiation, as well as those dysregulated ncRNAs involved in hematologic malignancies.

Mitotic Transcriptional Activation: Clearance of Actively Engaged Pol II via Transcriptional Elongation Control in Mitosis

Although it is established that some general transcription factors are inactivated at mitosis, many details of mitotic transcription inhibition (MTI) and its underlying mechanisms are largely unknown. We have identified mitotic transcriptional activation (MTA) as a key regulatory step to control transcription in mitosis for genes with transcriptionally engaged RNA polymerase II (Pol II) to activate and transcribe until the end of the gene to clear Pol II from mitotic chromatin, followed by global impairment of transcription reinitiation through MTI. Global nascent RNA sequencing and RNA fluorescence in situ hybridization demonstrate the existence of transcriptionally engaged Pol II in early mitosis. Both genetic and chemical inhibition of P-TEFb in mitosis lead to delays in the progression of cell division. Together, our study reveals a mechanism for MTA and MTI whereby transcriptionally engaged Pol II can progress into productive elongation and finish transcription to allow proper cellular division.

Histone H3K4 trimethylation: dynamic interplay with pre-mRNA splicing

Histone H3 lysine 4 trimethylation (H3K4me3) is often stated as a mark of transcriptionally active promoters. However, closer study of the positioning of H3K4me3 shows the mark locating primarily after the first exon at the 5′ splice site and overlapping with a CpG island in mammalian cells. There are several enzyme complexes that are involved in the placement of the H3K4me3 mark, including multiple protein complexes containing SETD1A, SETD1B, and MLL1 enzymes (writers). CXXC1, which is associated with SETD1A and SETD1B, target these enzymes to unmethylated CpG islands. Lysine demethylases (KDM5 family members, erasers) demethylate H3K4me3. The H3K4me3 mark is recognized by several proteins (readers), including lysine acetyltransferase complexes, chromatin remodelers, and RNA bound proteins involved in pre-mRNA splicing. Interestingly, attenuation of H3K4me3 impacts pre-mRNA splicing, and inhibition of pre-mRNA splicing attenuates H3K4me3.

Lineage-specific enhancers activate self-renewal genes in macrophages and embryonic stem cells

Differentiated macrophages can self-renew in tissues and expand long term in culture, but the gene regulatory mechanisms that accomplish self-renewal in the differentiated state have remained unknown. Here we show that in mice, the transcription factors MafB and c-Maf repress a macrophage-specific enhancer repertoire associated with a gene network that controls self-renewal. Single-cell analysis revealed that, in vivo, proliferating resident macrophages can access this network by transient down-regulation of Maf transcription factors. The network also controls embryonic stem cell self-renewal but is associated with distinct embryonic stem cell-specific enhancers. This indicates that distinct lineage-specific enhancer platforms regulate a shared network of genes that control self-renewal potential in both stem and mature cells.

Noncoding RNAs and pancreatic cancer

Noncoding RNAs (ncRNAs) represent a class of RNA molecules that typically do not code for proteins. Emerging data suggest that ncRNAs play an important role in several physiological and pathological conditions such as cancer. The best-characterized ncRNAs are the microRNAs (miRNAs), which are short, approximately 22-nucleotide sequences of RNA of approximately 22-nucleotide in length that regulate gene expression at the posttranscriptional level, through transcript degradation or translational repression. MiRNAs can function as master gene regulators, impacting a variety of cellular pathways important to normal cellular functions as well as cancer development and progression. In addition to miRNAs, long ncRNAs, which are transcripts longer than 200 nucleotides, have recently emerged as novel drivers of tumorigenesis. However, the molecular mechanisms of their regulation and function, and the significance of other ncRNAs such as piwi-interacting RNAs in pancreas carcinogenesis are largely unknown. This review summarizes the growing body of evidence supporting the vital roles of ncRNAs in pancreatic cancer, focusing on their dysregulation through both genetic and epigenetic mechanisms, and highlighting the promise of ncRNAs in diagnostic and therapeutic applications of pancreatic cancer.

Long non-coding RNA expression in primary human monocytes

Long non-coding RNAs (lncRNAs) have been shown to possess a wide range of functions in both cellular and developmental processes including cancers. Although some of the lncRNAs have been implicated in the regulation of the immune response, the exact function of the large majority of lncRNAs still remains unknown. In this study, we characterized the lncRNAs in human primary monocytes, an essential component of the innate immune system. We performed RNA sequencing of monocytes from four individuals and combined our data with eleven other publicly available datasets. Our analysis led to identification of ~8000 lncRNAs of which >1000 have not been previously reported in monocytes. PCR-based validation of a subset of the identified novel long intergenic noncoding RNAs (lincRNAs) revealed distinct expression patterns. Our study provides a landscape of lncRNAs in monocytes, which could facilitate future experimental studies to characterize the functions of these molecules in the innate immune system.

Deposition of 5-Methylcytosine on Enhancer RNAs Enables the Coactivator Function of PGC-1α

The Peroxisome proliferator-activated receptor-gamma coactivator 1 alpha (PGC-1α) is a transcriptional co-activator that plays a central role in adapted metabolic responses. PGC-1α is dynamically methylated and unmethylated at the residue K779 by the methyltransferase SET7/9 and the Lysine Specific Demethylase 1A (LSD1), respectively. Interactions of methylated PGC-1α[K779me] with the Spt-Ada-Gcn5-acetyltransferase (SAGA) complex, the Mediator members MED1 and MED17, and the NOP2/Sun RNA methytransferase 7 (NSUN7) reinforce transcription, and are concomitant with the m(5)C mark on enhancer RNAs (eRNAs). Consistently, loss of Set7/9 and NSun7 in liver cell model systems resulted in depletion of the PGC-1α target genes Pfkl, Sirt5, Idh3b, and Hmox2, which was accompanied by a decrease in the eRNAs levels associated with these loci. Enrichment of m(5)C within eRNA species coincides with metabolic stress of fasting in vivo. Collectively, these findings illustrate the complex epigenetic circuitry imposed by PGC-1α at the eRNA level to fine-tune energy metabolism.