Analysis of nascent RNA identifies a unified architecture of initiation regions at mammalian promoters and enhancers

Understanding Celiac Disease by Genomics

Celiac disease (CeD) is a complex immune-mediated disease. Genetic studies have implicated 43 predisposing loci that collectively explain some 50% of the genetic variance in CeD. More than ∼90% of CeD-associated single nucleotide polymorphisms (SNPs) localize to the non-coding genome, which we need to better understand to translate genetic knowledge into clinical practice. New genomic technologies and resources are permitting a systematic analysis of the functional elements in the non-coding part of the genome. Here we explain how investigating the regulatory and epigenomic landscape will help to pinpoint the cell types involved in CeD, and the driver genes and gene regulatory networks that are affected by CeD-associated SNPs.

Enhancers as non-coding RNA transcription units: recent insights and future perspectives

Networks of regulatory enhancers dictate distinct cell identities and cellular responses to diverse signals by instructing precise spatiotemporal patterns of gene expression. However, 35 years after their discovery, enhancer functions and mechanisms remain incompletely understood. Intriguingly, recent evidence suggests that many, if not all, functional enhancers are themselves transcription units, generating non-coding enhancer RNAs. This observation provides a fundamental insight into the inter-regulation between enhancers and promoters, which can both act as transcription units; it also raises crucial questions regarding the potential biological roles of the enhancer transcription process and non-coding enhancer RNAs. Here, we review research progress in this field and discuss several important, unresolved questions regarding the roles and mechanisms of enhancers in gene regulation.

Mammalian NET-seq analysis defines nascent RNA profiles and associated RNA processing genome-wide

The transcription cycle of RNA polymerase II (Pol II) correlates with changes to the phosphorylation state of its large subunit C-terminal domain (CTD). We recently developed Native Elongation Transcript sequencing using mammalian cells (mNET-seq), which generates single-nucleotide-resolution genome-wide profiles of nascent RNA and co-transcriptional RNA processing that are associated with different CTD phosphorylation states. Here we provide a detailed protocol for mNET-seq. First, Pol II elongation complexes are isolated with specific phospho-CTD antibodies from chromatin solubilized by micrococcal nuclease digestion. Next, RNA derived from within the Pol II complex is size fractionated and Illumina sequenced. Using mNET-seq, we have previously shown that Pol II pauses at both ends of protein-coding genes but with different CTD phosphorylation patterns, and we have also detected phosphorylation at serine 5 (Ser5-P) CTD-specific splicing intermediates and Pol II accumulation over co-transcriptionally spliced exons. With moderate biochemical and bioinformatic skills, mNET-seq can be completed in ∼6 d, not including sequencing and data analysis.

7SK-BAF axis controls pervasive transcription at enhancers

RNA functions at enhancers remain mysterious. Here we show that the 7SK small nuclear RNA (snRNA) inhibits enhancer transcription by modulating nucleosome position. 7SK occupies enhancers and super enhancers genome-wide in mouse and human cells, and 7SK is required to limit eRNA initiation and synthesis in a manner distinct from promoter pausing. Clustered elements at super enhancers uniquely require 7SK to prevent convergent transcription and DNA damage signaling. 7SK physically interacts with the BAF chromatin remodeling complex, recruit BAF to enhancers, and inhibits enhancer transcription by modulating chromatin structure. In turn, 7SK occupancy at enhancers coincides with Brd4 and is exquisitely sensitive to the bromodomain inhibitor JQ1. Thus, 7SK employs distinct mechanisms to counteract diverse consequences of pervasive transcription that distinguish super enhancers, enhancers, and promoters.

Genomic Organization of Human Transcription Initiation Complexes

A repertoire of transcription initiation factors engage the core promoter of mRNA genes to recruit RNA polymerase (Pol) II to initiate transcription, yet their precise spatial organization remains unclear. Using ChIP-exo, here we detail the interactions and genomic organization of initiation factors TBP, TFIIB, and Pol II at mRNA genes and within CpG islands. We find that when Pol II moves into a transcriptionally paused state, TBP/TFIIB remain at the promoter. We show that TBP and TFIIB bound to the core promoter at two separate, resolvable locations that coincided with sites of divergent transcription initiation. We also examine the precise binding of TBP at Pol III transcribed tRNA genes. We find that TBP crosslinked to tRNA genes in a similar manner as at Pol II transcribed genes. This comprehensive and high resolution genome-wide detection of the initiation machinery produces a consolidated view of transcription initiation events humans at Pol II coding and Pol III transcribed tRNA genes.

Regulation of mammalian transcription and splicing by Nuclear RNAi

RNA interference (RNAi) is well known as a mechanism for controlling mammalian mRNA translation in the cytoplasm, but what would be the consequences if it also functions in cell nuclei? Although RNAi has also been found in nuclei of plants, yeast, and other organisms, there has been relatively little progress towards understanding the potential involvement of mammalian RNAi factors in nuclear processes including transcription and splicing. This review summarizes evidence for mammalian RNAi factors in cell nuclei and mechanisms that might contribute to the control of gene expression. When RNAi factors bind small RNAs, they form ribonucleoprotein complexes that can be selective for target sequences within different classes of nuclear RNA substrates. The versatility of nuclear RNAi may supply a previously underappreciated layer of regulation to transcription, splicing, and other nuclear processes.

Genomic Views of Transcriptional Enhancers: Essential Determinants of Cellular Identity and Activity-Dependent Responses in the CNS

Sprinkled throughout the genome are a million regulatory sequences called transcriptional enhancers that activate gene promoters in the right cells, at the right time. Enhancers endow the brain with its incredible diversity of cell types and also translate neural activity into gene induction. Thanks to rapid advances in genomic technologies, it is now possible to identify thousands of enhancers rapidly, test their transcriptional function en masse, and address their neurobiological functions via genome editing. Enhancers also promise to be a great technological opportunity for neuroscience, offering the potential for cell-type-specific genetic labeling and manipulation without the need for transgenesis. The objective of this review and the accompanying 2015 SfN mini-symposium is to highlight the use of new and emerging genomic technologies to probe enhancer function in the nervous system.

SIGNIFICANCE STATEMENT

Transcriptional enhancers turn on genes in the right cells, at the right time. Enhancers are also the genomic sequences that encode the incredible diversity of cell types in the brain and enable the brain to turn genes on in response to new experiences. New technology enables enhancers to be found and manipulated. The study of enhancers promises to inform our understanding of brain development and function. The application of enhancer technology holds promise in accelerating basic neuroscience research and enabling gene therapies to be targeted to specific cell types in the brain.

Promoter architectures and developmental gene regulation

Core promoters are minimal regions sufficient to direct accurate initiation of transcription and are crucial for regulation of gene expression. They are highly diverse in terms of associated core promoter motifs, underlying sequence composition and patterns of transcription initiation. Distinctive features of promoters are also seen at the chromatin level, including nucleosome positioning patterns and presence of specific histone modifications. Recent advances in identifying and characterizing promoters using next-generation sequencing-based technologies have provided the basis for their classification into functional groups and have shed light on their modes of regulation, with important implications for transcriptional regulation in development. This review discusses the methodology and the results of genome-wide studies that provided insight into the diversity of RNA polymerase II promoter architectures in vertebrates and other Metazoa, and the association of these architectures with distinct modes of regulation in embryonic development and differentiation.

Chromatin structure and gene regulation: a dynamic view of enhancer function

Localized chromatin organization is now recognized as an important determinant of cell identity and developmental pathways. Recent studies have demonstrated that these epigenetic states are unexpectedly dynamic and malleable. In this Extra view we will highlight the transient nature of stimulus-induced enhancer accessibility and its importance for transcription regulation. Using glucocorticoid receptor (GR) as a model system we will discuss spatiotemporal relationships between receptor/chromatin interactions, lifetimes of the DNase I hypersensitivity sites (DHSs), long-range interactions, and gene regulation. We propose that differential temporal activation and utilization of distal regulatory elements plays a role in directing divergent stimulus-induced transcriptional programs.

The life history of retrocopies illuminates the evolution of new mammalian genes

New genes contribute substantially to adaptive evolutionary innovation, but the functional evolution of new mammalian genes has been little explored at a broad scale. Previous work established mRNA-derived gene duplicates, known as retrocopies, as models for the study of new gene origination. Here we combine mammalian transcriptomic and epigenomic data to unveil the processes underlying the evolution of stripped-down retrocopies into complex new genes. We show that although some robustly expressed retrocopies are transcribed from preexisting promoters, most evolved new promoters from scratch or recruited proto-promoters in their genomic vicinity. In particular, many retrocopy promoters emerged from ancestral enhancers (or bivalent regulatory elements) or are located in CpG islands not associated with other genes. We detected 88–280 selectively preserved retrocopies per mammalian species, illustrating that these mechanisms facilitated the birth of many functional retrogenes during mammalian evolution. The regulatory evolution of originally monoexonic retrocopies was frequently accompanied by exon gain, which facilitated co-option of distant promoters and allowed expression of alternative isoforms. While young retrogenes are often initially expressed in the testis, increased regulatory and structural complexities allowed retrogenes to functionally diversify and evolve somatic organ functions, sometimes as complex as those of their parents. Thus, some retrogenes evolved the capacity to temporarily substitute for their parents during the process of male meiotic X inactivation, while others rendered parental functions superfluous, allowing for parental gene loss. Overall, our reconstruction of the “life history” of mammalian retrogenes highlights retroposition as a general model for understanding new gene birth and functional evolution.

Identifying transcription start sites and active enhancer elements using BruUV-seq

BruUV-seq utilizes UV light to introduce transcription-blocking DNA lesions randomly in the genome prior to bromouridine-labeling and deep sequencing of nascent RNA. By inhibiting transcription elongation, but not initiation, pre-treatment with UV light leads to a redistribution of transcription reads resulting in the enhancement of nascent RNA signal towards the 5′-end of genes promoting the identification of transcription start sites (TSSs). Furthermore, transcripts associated with arrested RNA polymerases are protected from 3′–5′ degradation and thus, unstable transcripts such as putative enhancer RNA (eRNA) are dramatically increased. Validation of BruUV-seq against GRO-cap that identifies capped run-on transcripts showed that most BruUV-seq peaks overlapped with GRO-cap signal over both TSSs and enhancer elements. Finally, BruUV-seq identified putative enhancer elements induced by tumor necrosis factor (TNF) treatment concomitant with expression of nearby TNF-induced genes. Taken together, BruUV-seq is a powerful new approach for identifying TSSs and active enhancer elements genome-wide in intact cells.

Progress and challenges in bioinformatics approaches for enhancer identification

Enhancers are cis-acting DNA elements that play critical roles in distal regulation of gene expression. Identifying enhancers is an important step for understanding distinct gene expression programs that may reflect normal and pathogenic cellular conditions. Experimental identification of enhancers is constrained by the set of conditions used in the experiment. This requires multiple experiments to identify enhancers, as they can be active under specific cellular conditions but not in different cell types/tissues or cellular states. This has opened prospects for computational prediction methods that can be used for high-throughput identification of putative enhancers to complement experimental approaches. Potential functions and properties of predicted enhancers have been catalogued and summarized in several enhancer-oriented databases. Because the current methods for the computational prediction of enhancers produce significantly different enhancer predictions, it will be beneficial for the research community to have an overview of the strategies and solutions developed in this field. In this review, we focus on the identification and analysis of enhancers by bioinformatics approaches. First, we describe a general framework for computational identification of enhancers, present relevant data types and discuss possible computational solutions. Next, we cover over 30 existing computational enhancer identification methods that were developed since 2000. Our review highlights advantages, limitations and potentials, while suggesting pragmatic guidelines for development of more efficient computational enhancer prediction methods. Finally, we discuss challenges and open problems of this topic, which require further consideration.

Multiple Hepatic Regulatory Variants at the GALNT2 GWAS Locus Associated with High-Density Lipoprotein Cholesterol

Genome-wide association studies (GWASs) have identified more than 150 loci associated with blood lipid and cholesterol levels; however, the functional and molecular mechanisms for many associations are unknown. We examined the functional regulatory effects of candidate variants at the GALNT2 locus associated with high-density lipoprotein cholesterol (HDL-C). Fine-mapping and conditional analyses in the METSIM study identified a single locus harboring 25 noncoding variants (r(2) > 0.7 with the lead GWAS variants) strongly associated with total cholesterol in medium-sized HDL (e.g., rs17315646, p = 3.5 × 10(-12)). We used luciferase reporter assays in HepG2 cells to test all 25 variants for allelic differences in regulatory enhancer activity. rs2281721 showed allelic differences in transcriptional activity (75-fold [T] versus 27-fold [C] more than the empty-vector control), as did a separate 780-bp segment containing rs4846913, rs2144300, and rs6143660 (49-fold [AT(-) haplotype] versus 16-fold [CC(+) haplotype] more). Using electrophoretic mobility shift assays, we observed differential CEBPB binding to rs4846913, and we confirmed this binding in a native chromatin context by performing chromatin-immunoprecipitation (ChIP) assays in HepG2 and Huh-7 cell lines of differing genotypes. Additionally, sequence reads in HepG2 DNase-I-hypersensitivity and CEBPB ChIP-seq signals spanning rs4846913 showed significant allelic imbalance. Allelic-expression-imbalance assays performed with RNA from primary human hepatocyte samples and expression-quantitative-trait-locus (eQTL) data in human subcutaneous adipose tissue samples confirmed that alleles associated with increased HDL-C are associated with a modest increase in GALNT2 expression. Together, these data suggest that at least rs4846913 and rs2281721 play key roles in influencing GALNT2 expression at this HDL-C locus.

Architectural and Functional Commonalities between Enhancers and Promoters

With the explosion of genome-wide studies of regulated transcription, it has become clear that traditional definitions of enhancers and promoters need to be revisited. These control elements can now be characterized in terms of their local and regional architecture, their regulatory components, including histone modifications and associated binding factors, and their functional contribution to transcription. This Review discusses unifying themes between promoters and enhancers in transcriptional regulatory mechanisms.

Stimulus-specific combinatorial functionality of neuronal c-fos enhancers

The c-fos gene is induced by a broad range of stimuli, and has been commonly used as a reliable marker for neural activity. Its induction mechanism and available reporter mouse lines are exclusively based on the c-fos promoter activity. Here, we demonstrate that multiple enhancers surrounding the c-fos gene are critical for ensuring robust c-fos response to various stimuli. Membrane depolarization, brain-derived neurotrophic factor (BDNF), and Forskolin activate distinct subsets of the enhancers to induce c-fos transcription in neurons, suggesting that stimulus-specific combinatorial activation of multiple enhancers underlies the broad inducibility of the c-fos gene. Accordingly, the functional requirement of key transcription factors varies depending on the type of stimulation. Combinatorial enhancer activation also occurs in the brain. Providing a comprehensive picture of the c-fos induction mechanism beyond the minimal promoter, our study should help in understanding the physiological nature of c-fos induction in relation to neural activity and plasticity.

In the loop: promoter-enhancer interactions and bioinformatics

Enhancer-promoter regulation is a fundamental mechanism underlying differential transcriptional regulation. Spatial chromatin organization brings remote enhancers in contact with target promoters in cis to regulate gene expression. There is considerable evidence for promoter-enhancer interactions (PEIs). In the recent years, genome-wide analyses have identified signatures and mapped novel enhancers; however, being able to precisely identify their target gene(s) requires massive biological and bioinformatics efforts. In this review, we give a short overview of the chromatin landscape and transcriptional regulation. We discuss some key concepts and problems related to chromatin interaction detection technologies, and emerging knowledge from genome-wide chromatin interaction data sets. Then, we critically review different types of bioinformatics analysis methods and tools related to representation and visualization of PEI data, raw data processing and PEI prediction. Lastly, we provide specific examples of how PEIs have been used to elucidate a functional role of non-coding single-nucleotide polymorphisms. The topic is at the forefront of epigenetic research, and by highlighting some future bioinformatics challenges in the field, this review provides a comprehensive background for future PEI studies.

RNA:DNA hybrids in the human genome have distinctive nucleotide characteristics, chromatin composition, and transcriptional relationships

Background

RNA:DNA hybrids represent a non-canonical nucleic acid structure that has been associated with a range of human diseases and potential transcriptional regulatory functions. Mapping of RNA:DNA hybrids in human cells reveals them to have a number of characteristics that give insights into their functions.

Results

We find RNA:DNA hybrids to occupy millions of base pairs in the human genome. A directional sequencing approach shows the RNA component of the RNA:DNA hybrid to be purine-rich, indicating a thermodynamic contribution to their in vivo stability. The RNA:DNA hybrids are enriched at loci with decreased DNA methylation and increased DNase hypersensitivity, and within larger domains with characteristics of heterochromatin formation, indicating potential transcriptional regulatory properties. Mass spectrometry studies of chromatin at RNA:DNA hybrids shows the presence of the ILF2 and ILF3 transcription factors, supporting a model of certain transcription factors binding preferentially to the RNA:DNA conformation.

Conclusions

Overall, there is little to indicate a dependence for RNA:DNA hybrids forming co-transcriptionally, with results from the ribosomal DNA repeat unit instead supporting the intriguing model of RNA generating these structures intrans. The results of the study indicate heterogeneous functions of these genomic elements and new insights into their formation and stability in vivo.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-015-0040-6) contains supplementary material, which is available to authorized users.

The identification of cis-regulatory elements: A review from a machine learning perspective

The majority of the human genome consists of non-coding regions that have been called junk DNA. However, recent studies have unveiled that these regions contain cis-regulatory elements, such as promoters, enhancers, silencers, insulators, etc. These regulatory elements can play crucial roles in controlling gene expressions in specific cell types, conditions, and developmental stages. Disruption to these regions could contribute to phenotype changes. Precisely identifying regulatory elements is key to deciphering the mechanisms underlying transcriptional regulation. Cis-regulatory events are complex processes that involve chromatin accessibility, transcription factor binding, DNA methylation, histone modifications, and the interactions between them. The development of next-generation sequencing techniques has allowed us to capture these genomic features in depth. Applied analysis of genome sequences for clinical genetics has increased the urgency for detecting these regions. However, the complexity of cis-regulatory events and the deluge of sequencing data require accurate and efficient computational approaches, in particular, machine learning techniques. In this review, we describe machine learning approaches for predicting transcription factor binding sites, enhancers, and promoters, primarily driven by next-generation sequencing data. Data sources are provided in order to facilitate testing of novel methods. The purpose of this review is to attract computational experts and data scientists to advance this field.

Promoter-distal RNA polymerase II binding discriminates active from inactive CCAAT/ enhancer-binding protein beta binding sites

Transcription factors (TFs) bind to thousands of DNA sequences in mammalian genomes, but most of these binding events appear to have no direct effect on gene expression. It is unclear why only a subset of TF bound sites are actively involved in transcriptional regulation. Moreover, the key genomic features that accurately discriminate between active and inactive TF binding events remain ambiguous. Recent studies have identified promoter-distal RNA polymerase II (RNAP2) binding at enhancer elements, suggesting that these interactions may serve as a marker for active regulatory sequences. Despite these correlative analyses, a thorough functional validation of these genomic co-occupancies is still lacking. To characterize the gene regulatory activity of DNA sequences underlying promoter-distal TF binding events that co-occur with RNAP2 and TF sites devoid of RNAP2 occupancy using a functional reporter assay, we performed cis-regulatory element sequencing (CRE-seq). We tested more than 1000 promoter-distal CCAAT/enhancer-binding protein beta (CEBPB)-bound sites in HepG2 and K562 cells, and found that CEBPB-bound sites co-occurring with RNAP2 were more likely to exhibit enhancer activity. CEBPB-bound sites further maintained substantial cell-type specificity, indicating that local DNA sequence can accurately convey cell-type-specific regulatory information. By comparing our CRE-seq results to a comprehensive set of genome annotations, we identified a variety of genomic features that are strong predictors of regulatory element activity and cell-type-specific activity. Collectively, our functional assay results indicate that RNAP2 occupancy can be used as a key genomic marker that can distinguish active from inactive TF bound sites.

Chromatin Architecture of the Pitx2 Locus Requires CTCF- and Pitx2-Dependent Asymmetry that Mirrors Embryonic Gut Laterality

Expression of Pitx2 on the left side of the embryo patterns left-right (LR) organs including the dorsal mesentery (DM), whose asymmetric cell behavior directs gut looping. Despite the importance of organ laterality, chromatin-level regulation of Pitx2 remains undefined. Here, we show that genes immediately neighboring Pitx2 in chicken and mouse, including a long noncoding RNA (Pitx2 locus-asymmetric regulated RNA or Playrr), are expressed on the right side and repressed by Pitx2. CRISPR/Cas9 genome editing of Playrr, 3D fluorescent in situ hybridization (FISH), and variations of chromatin conformation capture (3C) demonstrate that mutual antagonism between Pitx2 and Playrr is coordinated by asymmetric chromatin interactions dependent on Pitx2 and CTCF. We demonstrate that transcriptional and morphological asymmetries driving gut looping are mirrored by chromatin architectural asymmetries at the Pitx2 locus. We propose a model whereby Pitx2 auto-regulation directs chromatin topology to coordinate LR transcription of this locus essential for LR organogenesis.

RNA polymerase II pausing as a context-dependent reader of the genome

The RNA polymerase II (Pol II) transcribes all mRNA genes in eukaryotes and is among the most highly regulated enzymes in the cell. The classic model of mRNA gene regulation involves recruitment of the RNA polymerase to gene promoters in response to environmental signals. Higher eukaryotes have an additional ability to generate multiple cell types. This extra level of regulation enables each cell to interpret the same genome by committing to one of the many possible transcription programs and executing it in a precise and robust manner. Whereas multiple mechanisms are implicated in cell type-specific transcriptional regulation, how one genome can give rise to distinct transcriptional programs and what mechanisms activate and maintain the appropriate program in each cell remains unclear. This review focuses on the process of promoter-proximal Pol II pausing during early transcription elongation as a key step in context-dependent interpretation of the metazoan genome. We highlight aspects of promoter-proximal Pol II pausing, including its interplay with epigenetic mechanisms, that may enable cell type-specific regulation, and emphasize some of the pertinent questions that remain unanswered and open for investigation.

Analysis of paired end Pol II ChIP-seq and short capped RNA-seq in MCF-7 cells

While a role of promoter-proximal RNA Polymerase II (Pol II) pausing in regulation of eukaryotic gene expression is implied, the mechanisms and dynamics of this process are poorly understood. We performed genome-wide analysis of short capped RNAs (scRNAs) and Pol II chromatin immunoprecipitation sequencing (ChIP-seq) in human breast cancer MCF-7 cells to better understand Pol II pausing (Samarakkody, A., Abbas, A., Scheidegger, A., Warns, J., Nnoli, O., Jokinen, B., Zarns, K., Kubat, B., Dhasarathy, A. and Nechaev, S. (2015) RNA polymerase II pausing can be retained or acquired during activation of genes involved in the epithelial to mesenchymal transition. Nucleic Acids Res43, 3938–3949). The data are available at the NCBI Gene Expression Omnibus under accession number GSE67041. For both ChIP and scRNA samples, we used paired end sequencing on the Illumina MiSeq instrument. For ChIP-seq, the use of paired end sequencing allowed us to avoid ambiguities in center-read definition. For scRNA seq, this allowed us to identify both the 5′-end and the 3′-end in the same run that represent, respectively, the transcription start sites and the locations of Pol II pausing. The sharpening of Pol II ChIP-seq metagene profiles when aligned against 5′-ends of scRNAs indicates that these RNAs can be used to define the start sites for the majority of mRNA transcription events.

Exploiting genomics and natural genetic variation to decode macrophage enhancers

The mammalian genome contains on the order of a million enhancer-like regions that are required to establish the identities and functions of specific cell types. Here, we review recent studies in immune cells that have provided insight into the mechanisms that selectively activate certain enhancers in response to cell lineage and environmental signals. We describe a working model wherein distinct classes of transcription factors define the repertoire of active enhancers in macrophages through collaborative and hierarchical interactions, and discuss important challenges to this model, specifically providing examples from T cells. We conclude by discussing the use of natural genetic variation as a powerful approach for decoding transcription factor combinations that play dominant roles in establishing the enhancer landscapes, and the potential that these insights have for advancing our understanding of the molecular causes of human disease.

An Overview of Genome Organization and How We Got There: from FISH to Hi-C

In humans, nearly two meters of genomic material must be folded to fit inside each micrometer-scale cell nucleus while remaining accessible for gene transcription, DNA replication, and DNA repair. This fact highlights the need for mechanisms governing genome organization during any activity and to maintain the physical organization of chromosomes at all times. Insight into the functions and three-dimensional structures of genomes comes mostly from the application of visual techniques such as fluorescence in situ hybridization (FISH) and molecular approaches including chromosome conformation capture (3C) technologies. Recent developments in both types of approaches now offer the possibility of exploring the folded state of an entire genome and maybe even the identification of how complex molecular machines govern its shape. In this review, we present key methodologies used to study genome organization and discuss what they reveal about chromosome conformation as it relates to transcription regulation across genomic scales in mammals.

Why the activity of a gene depends on its neighbors

Sixty years ago, the position of a gene on a chromosome was seen to be a major determinant of gene activity; however, position effects are rarely central to current discussions of gene expression. We describe a comprehensive and simplifying view of how position in 1D sequence and 3D nuclear space underlies expression. We suggest that apparently-different regulatory motifs including enhancers, silencers, insulators, barriers, and boundaries act similarly - they are active promoters that tether target genes close to, or distant from, appropriate transcription sites or 'factories'. We also suggest that any active transcription unit regulates the firing of its neighbors - and thus can be categorized as one or other type of motif; this is consistent with expression quantitative trait loci (eQTLs) being widely dispersed.

THZ1 Reveals Roles for Cdk7 in Co-transcriptional Capping and Pausing

The Cdk7 subunit of TFIIH phosphorylates RNA polymerase II (Pol II) during initiation, and, while recent studies show that inhibition of human Cdk7 negatively influences transcription, the mechanisms involved are unclear. Using in vitro transcription with nuclear extract, we demonstrate that THZ1, a covalent Cdk7 inhibitor, causes defects in Pol II phosphorylation, co-transcriptional capping, promoter proximal pausing, and productive elongation. THZ1 does not affect initiation but blocks essentially all Pol II large subunit C-terminal domain (CTD) phosphorylation. We found that guanylylation of nascent RNAs is length dependent and modulated by a THZ1-sensitive factor present in nuclear extract. THZ1 impacts pausing through a capping-independent block of DSIF and NELF loading. The P-TEFb-dependent transition into productive elongation was also inhibited by THZ1, likely due to loss of DSIF. Capping and pausing were also reduced in THZ1-treated cells. Our results provide mechanistic insights into THZ1 action and how Cdk7 broadly influences transcription and capping.

Ready, pause, go: regulation of RNA polymerase II pausing and release by cellular signaling pathways

Promoter-proximal pausing by RNA polymerase II (Pol II) is a well-established mechanism to control the timing, rate, and possibly the magnitude of transcriptional responses. Recent studies have shown that cellular signaling pathways can regulate gene transcription and signaling outcomes by controlling Pol II pausing in a wide array of biological systems. Identification of the proteins and small molecules that affect the establishment and release of paused Pol II is shedding new light on the mechanisms and biology of Pol II pausing. This review focuses on the interplay between cellular signaling pathways and Pol II pausing during normal development and under disease conditions.

Coexpression networks identify brain region-specific enhancer RNAs in the human brain

Despite major progress in identifying enhancer regions on a genome-wide scale, the majority of available data are limited to model organisms and human transformed cell lines. We have identified a robust set of enhancer RNAs (eRNAs) expressed in the human brain and constructed networks assessing eRNA-gene coexpression interactions across human fetal brain and multiple adult brain regions. Our data identify brain region-specific eRNAs and show that enhancer regions expressing eRNAs are enriched for genetic variants associated with autism spectrum disorders.

Differential and coherent processing patterns from small RNAs

Post-transcriptional processing events related to short RNAs are often reflected in their read profile patterns emerging from high-throughput sequencing data. MicroRNA arm switching across different tissues is a well-known example of what we define as differential processing. Here, short RNAs from the nine cell lines of the ENCODE project, irrespective of their annotation status, were analyzed for genomic loci representing differential or coherent processing. We observed differential processing predominantly in RNAs annotated as miRNA, snoRNA or tRNA. Four out of five known cases of differentially processed miRNAs that were in the input dataset were recovered and several novel cases were discovered. In contrast to differential processing, coherent processing is observed widespread in both annotated and unannotated regions. While the annotated loci predominantly consist of ~24 nt short RNAs, the unannotated loci comparatively consist of ~17 nt short RNAs. Furthermore, these ~17 nt short RNAs are significantly enriched for overlap to transcription start sites and DNase I hypersensitive sites (p-value < 0.01) that are characteristic features of transcription initiation RNAs. We discuss how the computational pipeline developed in this study has the potential to be applied to other forms of RNA-seq data for further transcriptome-wide studies of differential and coherent processing.

Genome-wide analysis of human global and transcription-coupled excision repair of UV damage at single-nucleotide resolution

Hu et al. developed a method for genome-wide mapping of DNA excision repair named XR-seq (excision repair sequencing) and used it to produce stranded, nucleotide-resolution maps of repair of two UV-induced DNA damages in human cells. XR-seq and the resulting repair maps will facilitate studies of the effects of genomic location, chromatin context, transcription, and replication on DNA repair in human cells.

We developed a method for genome-wide mapping of DNA excision repair named XR-seq (excision repair sequencing). Human nucleotide excision repair generates two incisions surrounding the site of damage, creating an ∼30-mer. In XR-seq, this fragment is isolated and subjected to high-throughput sequencing. We used XR-seq to produce stranded, nucleotide-resolution maps of repair of two UV-induced DNA damages in human cells: cyclobutane pyrimidine dimers (CPDs) and (6-4) pyrimidine–pyrimidone photoproducts [(6-4)PPs]. In wild-type cells, CPD repair was highly associated with transcription, specifically with the template strand. Experiments in cells defective in either transcription-coupled excision repair or general excision repair isolated the contribution of each pathway to the overall repair pattern and showed that transcription-coupled repair of both photoproducts occurs exclusively on the template strand. XR-seq maps capture transcription-coupled repair at sites of divergent gene promoters and bidirectional enhancer RNA (eRNA) production at enhancers. XR-seq data also uncovered the repair characteristics and novel sequence preferences of CPDs and (6-4)PPs. XR-seq and the resulting repair maps will facilitate studies of the effects of genomic location, chromatin context, transcription, and replication on DNA repair in human cells.

A unified architecture of transcriptional regulatory elements

Gene expression is precisely controlled in time and space through the integration of signals that act at gene promoters and gene-distal enhancers. Classically, promoters and enhancers are considered separate classes of regulatory elements, often distinguished by histone modifications. However, recent studies have revealed broad similarities between enhancers and promoters, blurring the distinction: active enhancers often initiate transcription, and some gene promoters have the potential to enhance transcriptional output of other promoters. Here, we propose a model in which promoters and enhancers are considered a single class of functional element, with a unified architecture for transcription initiation. The context of interacting regulatory elements and the surrounding sequences determine local transcriptional output as well as the enhancer and promoter activities of individual elements.

Protein-DNA binding in high-resolution

Recent advances in experimental and computational methodologies are enabling ultra-high resolution genome-wide profiles of protein-DNA binding events. For example, the ChIP-exo protocol precisely characterizes protein-DNA cross-linking patterns by combining chromatin immunoprecipitation (ChIP) with 5' → 3' exonuclease digestion. Similarly, deeply sequenced chromatin accessibility assays (e.g. DNase-seq and ATAC-seq) enable the detection of protected footprints at protein-DNA binding sites. With these techniques and others, we have the potential to characterize the individual nucleotides that interact with transcription factors, nucleosomes, RNA polymerases and other regulatory proteins in a particular cellular context. In this review, we explain the experimental assays and computational analysis methods that enable high-resolution profiling of protein-DNA binding events. We discuss the challenges and opportunities associated with such approaches.

Eukaryotic enhancers: common features, regulation, and participation in diseases

Enhancers are positive DNA regulatory sequences controlling temporal and tissue-specific gene expression. These elements act independently of their orientation and distance relative to the promoters of target genes. Enhancers act through a variety of transcription factors that ensure their correct match with target promoters and consequent gene activation. There is a growing body of evidence on association of enhancers with transcription factors, co-activators, histone chromatin marks, and lncRNAs. Alterations in enhancers lead to misregulation of gene expression, causing a number of human diseases. In this review, we focus on the common characteristics of enhancers required for transcription stimulation.

Bidirectional Transcription Arises from Two Distinct Hubs of Transcription Factor Binding and Active Chromatin

Anti-sense transcription originating upstream of mammalian protein-coding genes is a well-documented phenomenon, but remarkably little is known about the regulation or function of anti-sense promoters and the non-coding RNAs they generate. Here we define at nucleotide resolution the divergent transcription start sites (TSSs) near mouse mRNA genes. We find that coupled sense and anti-sense TSSs precisely define the boundaries of a nucleosome-depleted region (NDR) that is highly enriched in transcription factor (TF) motifs. Notably, as the distance between sense and anti-sense TSSs increases, so does the size of the NDR, the level of signal-dependent TF binding, and gene activation. We further discover a group of anti-sense TSSs in macrophages with an enhancer-like chromatin signature. Interestingly, this signature identifies divergent promoters that are activated during immune challenge. We propose that anti-sense promoters serve as platforms for TF binding and establishment of active chromatin to further regulate or enhance sense-strand mRNA expression.

Identification of active transcriptional regulatory elements from GRO-seq data

Modifications to the global run-on and sequencing (GRO-seq) protocol that enrich for 5'-capped RNAs can be used to reveal active transcriptional regulatory elements (TREs) with high accuracy. Here, we introduce discriminative regulatory-element detection from GRO-seq (dREG), a sensitive machine learning method that uses support vector regression to identify active TREs from GRO-seq data without requiring cap-based enrichment (https://github.com/Danko-Lab/dREG/). This approach allows TREs to be assayed together with gene expression levels and other transcriptional features in a single experiment. Predicted TREs are more enriched for several marks of transcriptional activation—including expression quantitative trait loci, disease-associated polymorphisms, acetylated histone 3 lysine 27 (H3K27ac) and transcription factor binding—than those identified by alternative functional assays. Using dREG, we surveyed TREs in eight human cell types and provide new insights into global patterns of TRE function.

The core promoter: At the heart of gene expression

The identities of different cells and tissues in multicellular organisms are determined by tightly controlled transcriptional programs that enable accurate gene expression. The mechanisms that regulate gene expression comprise diverse multiplayer molecular circuits of multiple dedicated components. The RNA polymerase II (Pol II) core promoter establishes the center of this spatiotemporally orchestrated molecular machine. Here, we discuss transcription initiation, diversity in core promoter composition, interactions of the basal transcription machinery with the core promoter, enhancer-promoter specificity, core promoter-preferential activation, enhancer RNAs, Pol II pausing, transcription termination, Pol II recycling and translation. We further discuss recent findings indicating that promoters and enhancers share similar features and may not substantially differ from each other, as previously assumed. Taken together, we review a broad spectrum of studies that highlight the importance of the core promoter and its pivotal role in the regulation of metazoan gene expression and suggest future research directions and challenges.

Enhancer-associated RNAs as therapeutic targets

INTRODUCTION

Regulation of gene expression involves a variety of mechanisms driven by a complex regulatory network of factors. Control of transcription is an important step in gene expression regulation, which integrates the function of cis-acting and trans-acting elements. Among cis-regulatory elements, enhancer RNA (eRNA)-producing domains recently emerged as widespread and potent regulators of transcription and cell fate decision. Thus, manipulation of eRNA levels becomes a novel and appealing avenue for the design of new therapeutic treatments.

AREAS COVERED

In this review, we focus on eRNA-producing domains. We describe mechanisms involved in their cell-type specific selection and activation as well as their epigenetic features. In addition, we present their function and the growing evidences of their deregulation in human diseases. Finally, we discuss eRNAs as potential therapeutic targets.

EXPERT OPINION

As key factors in the control of transcription, eRNAs appear to possess a great potential for the establishment of new therapy options. However, thorough testing as well as providing the genetic toolbox to target eRNAs will be needed to fully assess the practical and clinical possibilities.

RNA polymerase II pausing can be retained or acquired during activation of genes involved in the epithelial to mesenchymal transition

Promoter-proximal RNA polymerase II (Pol II) pausing is implicated in the regulation of gene transcription. However, the mechanisms of pausing including its dynamics during transcriptional responses remain to be fully understood. We performed global analysis of short capped RNAs and Pol II Chromatin Immunoprecipitation sequencing in MCF-7 breast cancer cells to map Pol II pausing across the genome, and used permanganate footprinting to specifically follow pausing during transcriptional activation of several genes involved in the epithelial to mesenchymal transition (EMT). We find that the gene for EMT master regulator Snail (SNAI1), but not Slug (SNAI2), shows evidence of Pol II pausing before activation. Transcriptional activation of the paused SNAI1 gene is accompanied by a further increase in Pol II pausing signal, whereas activation of non-paused SNAI2 gene results in the acquisition of a typical pausing signature. The increase in pausing signal reflects increased transcription initiation without changes in Pol II pausing. Activation of the heat shock HSP70 gene involves pausing release that speeds up Pol II turnover, but does not change pausing location. We suggest that Pol II pausing is retained during transcriptional activation and can further undergo regulated release in a signal-specific manner.

Getting up to speed with transcription elongation by RNA polymerase II

Recent advances in sequencing techniques that measure nascent transcripts and that reveal the positioning of RNA polymerase II (Pol II) have shown that the pausing of Pol II in promoter-proximal regions and its release to initiate a phase of productive elongation are key steps in transcription regulation. Moreover, after the release of Pol II from the promoter-proximal region, elongation rates are highly dynamic throughout the transcription of a gene, and vary on a gene-by-gene basis. Interestingly, Pol II elongation rates affect co-transcriptional processes such as splicing, termination and genome stability. Increasing numbers of factors and regulatory mechanisms have been associated with the steps of transcription elongation by Pol II, revealing that elongation is a highly complex process. Elongation is thus now recognized as a key phase in the regulation of transcription by Pol II.

The selection and function of cell type-specific enhancers

The human body contains several hundred cell types, all of which share the same genome. In metazoans, much of the regulatory code that drives cell type-specific gene expression is located in distal elements called enhancers. Although mammalian genomes contain millions of potential enhancers, only a small subset of them is active in a given cell type. Cell type-specific enhancer selection involves the binding of lineage-determining transcription factors that prime enhancers. Signal-dependent transcription factors bind to primed enhancers, which enables these broadly expressed factors to regulate gene expression in a cell type-specific manner. The expression of genes that specify cell type identity and function is associated with densely spaced clusters of active enhancers known as super-enhancers. The functions of enhancers and super-enhancers are influenced by, and affect, higher-order genomic organization.

Core promoters in transcription: old problem, new insights

Early studies established that transcription initiates within an approximately 50 bp DNA segment capable of nucleating the assembly of RNA polymerase II (Pol II) and associated general transcription factors (GTFs) necessary for transcriptional initiation; this region is called a core promoter. Subsequent analyses identified a series of conserved DNA sequence elements, present in various combinations or not at all, in core promoters. Recent genome-wide analyses have provided further insights into the complexity of core promoter architecture and function. Here we review recent studies that delineate the active role of core promoters in the transcriptional regulation of diverse physiological systems.

A shared architecture for promoters and enhancers

A new study detects unstable nascent RNAs and uncovers thousands of transcription initiation sites in promoters and enhancers. Detailed analysis shows that these initiation sites have a similar architecture and that they are differentiated by post-transcriptional regulation rather than transcription initiation.

Sequencing of first-strand cDNA library reveals full-length transcriptomes

Massively parallel strand-specific sequencing of RNA (ssRNA-seq) has emerged as a powerful tool for profiling complex transcriptomes. However, many current methods for ssRNA-seq suffer from the underrepresentation of both the 5' and 3' ends of RNAs, which can be attributed to second-strand cDNA synthesis. The 5' and 3' ends of RNA harbour crucial information for gene regulation; namely, transcription start sites (TSSs) and polyadenylation sites. Here we report a novel ssRNA-seq method that does not involve second-strand cDNA synthesis, as we Directly Ligate sequencing Adaptors to the First-strand cDNA (DLAF). This novel method with fewer enzymatic reactions results in a higher quality of the libraries than the conventional method. Sequencing of DLAF libraries followed by a novel analysis pipeline enables the profiling of both 5' ends and polyadenylation sites at near-base resolution. Therefore, DLAF offers the first genomics tool to obtain the 'full-length' transcriptome with a single library.

A method for calculating probabilities of fitness consequences for point mutations across the human genome

We describe a novel computational method for estimating the probability that a point mutation at each position in a genome will influence fitness. These fitness consequence (fitCons) scores serve as evolution-based measures of potential genomic function. Our approach is to cluster genomic positions into groups exhibiting distinct “fingerprints” based on high-throughput functional genomic data, then to estimate a probability of fitness consequences for each group from associated patterns of genetic polymorphism and divergence. We have generated fitCons scores for three human cell types based on public data from ENCODE. Compared with conventional conservation scores, fitCons scores show considerably improved prediction power for cis-regulatory elements. In addition, fitCons scores indicate that 4.2–7.5% of nucleotides in the human genome have influenced fitness since the human-chimpanzee divergence, and they suggest that recent evolutionary turnover has had limited impact on the functional content of the genome.