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

Convergence of the dysregulated regulome in schizophrenia with polygenic risk and evolutionarily constrained enhancers

Enhancers play an essential role in the etiology of schizophrenia; however, the dysregulation of enhancer activity and its impact on the regulome in schizophrenia remains understudied. To address this gap in our knowledge, we assessed enhancer and gene expression in 1,382 brain samples comprising cases with schizophrenia and unaffected controls. Dysregulation of enhancer expression was concordant with changes in gene expression, and was more closely associated with schizophrenia polygenic risk, suggesting that enhancer dysregulation is proximal to the genetic etiology of the disease. Modeling the shared variance of cis-coordinated genes and enhancers revealed a gene regulatory program that was highly associated with genetic vulnerability to schizophrenia. By integrating coordinated factors with evolutionary constraints, we found that enhancers acquired during human evolution are more likely to regulate genes that are implicated in neuropsychiatric disorders and, thus, hold potential as therapeutic targets. Our analysis provides a systematic view of regulome dysregulation in schizophrenia and highlights its convergence with schizophrenia polygenic risk and human-gained enhancers.

Histone deacetylases maintain expression of the pluripotent gene network via recruitment of RNA polymerase II to coding and noncoding loci

Histone acetylation is a dynamic modification regulated by the opposing actions of histone acetyltransferases (HATs) and histone deacetylases (HDACs). Deacetylation of histone tails results in chromatin tightening, and therefore, HDACs are generally regarded as transcriptional repressors. Counterintuitively, simultaneous deletion of Hdac1 and Hdac2 in embryonic stem cells (ESCs) reduces expression of the pluripotency-associated transcription factors Pou5f1, Sox2, and Nanog (PSN). By shaping global histone acetylation patterns, HDACs indirectly regulate the activity of acetyl-lysine readers, such as the transcriptional activator BRD4. Here, we use inhibitors of HDACs and BRD4 (LBH589 and JQ1, respectively) in combination with precision nuclear run-on and sequencing (PRO-seq) to examine their roles in defining the ESC transcriptome. Both LBH589 and JQ1 cause a marked reduction in the pluripotent gene network. However, although JQ1 treatment induces widespread transcriptional pausing, HDAC inhibition causes a reduction in both paused and elongating polymerase, suggesting an overall reduction in polymerase recruitment. Using enhancer RNA (eRNA) expression to measure enhancer activity, we find that LBH589-sensitive eRNAs are preferentially associated with superenhancers and PSN binding sites. These findings suggest that HDAC activity is required to maintain pluripotency by regulating the PSN enhancer network via the recruitment of RNA polymerase II.

An Epigenomic fingerprint of human cancers by landscape interrogation of super enhancers at the constituent level

Super enhancers (SE), large genomic elements that activate transcription and drive cell identity, have been found with cancer-specific gene regulation in human cancers. Recent studies reported the importance of understanding the cooperation and function of SE internal components, i.e., the constituent enhancers (CE). However, there are no pan-cancer studies to identify cancer-specific SE signatures at the constituent level. Here, by revisiting pan-cancer SE activities with H3K27Ac ChIP-seq datasets, we report fingerprint SE signatures for 28 cancer types in the NCI-60 cell panel. We implement a mixture model to discriminate active CEs from inactive CEs by taking into consideration ChIP-seq variabilities between cancer samples and across CEs. We demonstrate that the model-based estimation of CE states provides improved functional interpretation of SE-associated regulation. We identify cancer-specific CEs by balancing their active prevalence with their capability of encoding cancer type identities. We further demonstrate that cancer-specific CEs have the strongest per-base enhancer activities in independent enhancer sequencing assays, suggesting their importance in understanding critical SE signatures. We summarize fingerprint SEs based on the cancer-specific statuses of their component CEs and build an easy-to-use R package to facilitate the query, exploration, and visualization of fingerprint SEs across cancers.

Participation of the TBP-associated factors (TAFs) in cell differentiation

The understanding of the mechanisms that regulate gene expression to establish differentiation programs and determine cell lineages, is one of the major challenges in Developmental Biology. Besides the participation of tissue-specific transcription factors and epigenetic processes, the role of general transcription factors has been ignored. Only in recent years, there have been scarce studies that address this issue. Here, we review the studies on the biological activity of some TATA-box binding protein (TBP)-associated factors (TAFs) during the proliferation of stem/progenitor cells and their involvement in cell differentiation. Particularly, the accumulated evidence suggests that TAF4, TAF4b, TAF7L, TAF8, TAF9, and TAF10, among others, participate in nervous system development, adipogenesis, myogenesis, and epidermal differentiation; while TAF1, TAF7, TAF15 may be involved in the regulation of stem cell proliferative abilities and cell cycle progression. On the other hand, evidence suggests that TBP variants such as TBPL1 and TBPL2 might be regulating some developmental processes such as germ cell maturation and differentiation, myogenesis, or ventral specification during development. Our analysis shows that it is necessary to study in greater depth the biological function of these factors and its participation in the assembly of specific transcription complexes that contribute to the differential gene expression that gives rise to the great diversity of cell types existing in an organism. The understanding of TAFs' regulation might lead to the development of new therapies for patients which suffer from mutations, alterations, and dysregulation of these essential elements of the transcriptional machinery.

TE-TSS: an integrated data resource of human and mouse transposable element (TE)-derived transcription start site (TSS)

Transposable elements (TEs) are abundant in the genome and serve as crucial regulatory elements. Some TEs function as epigenetically regulated promoters, and these TE-derived transcription start sites (TSSs) play a crucial role in regulating genes associated with specific functions, such as cancer and embryogenesis. However, the lack of an accessible database that systematically gathers TE-derived TSS data is a current research gap. To address this, we established TE-TSS, an integrated data resource of human and mouse TE-derived TSSs (http://xozhanglab.com/TETSS). TE-TSS has compiled 2681 RNA sequencing datasets, spanning various tissues, cell lines and developmental stages. From these, we identified 5768 human TE-derived TSSs and 2797 mouse TE-derived TSSs, with 47% and 38% being experimentally validated, respectively. TE-TSS enables comprehensive exploration of TSS usage in diverse samples, providing insights into tissue-specific gene expression patterns and transcriptional regulatory elements. Furthermore, TE-TSS compares TE-derived TSS regions across 15 mammalian species, enhancing our understanding of their evolutionary and functional aspects. The establishment of TE-TSS facilitates further investigations into the roles of TEs in shaping the transcriptomic landscape and offers valuable resources for comprehending their involvement in diverse biological processes.

DNA-sequence and epigenomic determinants of local rates of transcription elongation

Across all branches of life, transcription elongation is a crucial, regulated phase in gene expression. Many recent studies in eukaryotes have focused on the regulation of promoter-proximal pausing of RNA Polymerase II (Pol II), but rates of productive elongation also vary substantially throughout the gene body, both within and across genes. Here, we introduce a probabilistic model for systematically evaluating potential determinants of the local elongation rate based on nascent RNA sequencing (NRS) data. Our model is derived from a unified model for both the kinetics of Pol II movement along the DNA template and the generation of NRS read counts at steady state. It allows for a continuously variable elongation rate along the gene body, with the rate at each nucleotide defined by a generalized linear relationship with nearby genomic and epigenomic features. High-dimensional feature vectors are accommodated through a sparse-regression extension. We show with simulations that the model allows accurate detection of associated features and accurate prediction of local elongation rates. In an analysis of public PRO-seq and epigenomic data, we identify several features that are strongly associated with reductions in the local elongation rate, including DNA methylation, splice sites, RNA stem-loops, CTCF binding sites, and several histone marks, including H3K36me3 and H4K20me1. By contrast, low-complexity sequences and H3K79me2 marks are associated with increases in elongation rate. In an analysis of DNA k -mers, we find that cytosine nucleotides are strongly associated with reductions in local elongation rate, particularly when preceded by guanines and followed by adenines or thymines. Increases in elongation rate are associated with thymines and A+T-rich k -mers. These associations are generally shared across cell types, and by considering them our model is effective at predicting features of held-out PRO-seq data. Overall, our analysis is the first to permit genome-wide predictions of relative nucleotide-specific elongation rates based on complex sets of genomic and epigenomic covariates. We have made predictions available for the K562, CD14+, MCF-7, and HeLa-S3 cell types in a UCSC Genome Browser track.

The human genome contains over a million autonomous exons

Mammalian mRNA and lncRNA exons are often small compared to introns. The exon definition model predicts that exons splice autonomously, dependent on proximal exon sequence features, explaining their delineation within large introns. This model has not been examined on a genome-wide scale, however, leaving open the question of how often mRNA and lncRNA exons are autonomous. It is also unknown how frequently such exons can arise by chance. Here, we directly assayed large fragments (500-1000 bp) of the human genome by exon trapping, which detects exons spliced into a heterologous transgene, here designed with a large intron context. We define the trapped exons as "autonomous." We obtained ∼1.25 million trapped exons, including most known mRNA and well-annotated lncRNA internal exons, demonstrating that human exons are predominantly autonomous. mRNA exons are trapped with the highest efficiency. Nearly a million of the trapped exons are unannotated, most located in intergenic regions and antisense to mRNA, with depletion from the forward strand of introns. These exons are not conserved, suggesting they are nonfunctional and arose from random mutations. They are nonetheless highly enriched with known splicing promoting sequence features that delineate known exons. Novel autonomous exons are more numerous than annotated lncRNA exons, and computational models also indicate they will occur with similar frequency in any randomly generated sequence. These results show that most human coding exons splice autonomously, and provide an explanation for the existence of many unconserved lncRNAs, as well as a new annotation and inclusion levels of spliceable loci in the human genome.

Regulatory effect of NK homeobox 1 (NKX2.1) on melanocortin 4 receptor (Mc4r) promoter in Mandarin fish

The melanocortin 4 receptor (MC4R) is a G protein-coupled transporter that mediates the regulation of thyroid hormones and leptin on energy balance and food intake. However, the mechanisms of transcriptional regulation of Mc4r by thyroid hormone and leptin in fish have been rarely reported. The messenger RNA expression of Mc4r gene was significantly higher in brain than those in other tissues of mandarin fish. We analyzed the structure and function of a 2029 bp sequence of Mc4r promoter. Meanwhile, overexpression of NKX2.1 and incubation with leptin significantly increased Mc4r promoter activity, but triiodothyronine showed the opposite effect. In addition, mutations in the NKX2.1 binding site abolished not only the activation of Mc4r promoter activity by leptin but also the inhibitory effect of thyroid hormones on Mc4r promoter activity. In summary, these results suggested that thyroid hormones and leptin might regulate the transcriptional expression of Mc4r through NKX2.1.

Finding Needles in the Haystack: Strategies for Uncovering Noncoding Regulatory Variants

Despite accumulating evidence implicating noncoding variants in human diseases, unraveling their functionality remains a significant challenge. Systematic annotations of the regulatory landscape and the growth of sequence variant data sets have fueled the development of tools and methods to identify causal noncoding variants and evaluate their regulatory effects. Here, we review the latest advances in the field and discuss potential future research avenues to gain a more in-depth understanding of noncoding regulatory variants.

Global identification of SWI/SNF targets reveals compensation by EP400

Mammalian SWI/SNF chromatin remodeling complexes move and evict nucleosomes at gene promoters and enhancers to modulate DNA access. Although SWI/SNF subunits are commonly mutated in disease, therapeutic options are limited by our inability to predict SWI/SNF gene targets and conflicting studies on functional significance. Here, we leverage a fast-acting inhibitor of SWI/SNF remodeling to elucidate direct targets and effects of SWI/SNF. Blocking SWI/SNF activity causes a rapid and global loss of chromatin accessibility and transcription. Whereas repression persists at most enhancers, we uncover a compensatory role for the EP400/TIP60 remodeler, which reestablishes accessibility at most promoters during prolonged loss of SWI/SNF. Indeed, we observe synthetic lethality between EP400 and SWI/SNF in cancer cell lines and human cancer patient data. Our data define a set of molecular genomic features that accurately predict gene sensitivity to SWI/SNF inhibition in diverse cancer cell lines, thereby improving the therapeutic potential of SWI/SNF inhibitors.

CHEX-seq detects single-cell genomic single-stranded DNA with catalytical potential

Genomic DNA (gDNA) undergoes structural interconversion between single- and double-stranded states during transcription, DNA repair and replication, which is critical for cellular homeostasis. We describe "CHEX-seq" which identifies the single-stranded DNA (ssDNA) in situ in individual cells. CHEX-seq uses 3'-terminal blocked, light-activatable probes to prime the copying of ssDNA into complementary DNA that is sequenced, thereby reporting the genome-wide single-stranded chromatin landscape. CHEX-seq is benchmarked in human K562 cells, and its utilities are demonstrated in cultures of mouse and human brain cells as well as immunostained spatially localized neurons in brain sections. The amount of ssDNA is dynamically regulated in response to perturbation. CHEX-seq also identifies single-stranded regions of mitochondrial DNA in single cells. Surprisingly, CHEX-seq identifies single-stranded loci in mouse and human gDNA that catalyze porphyrin metalation in vitro, suggesting a catalytic activity for genomic ssDNA. We posit that endogenous DNA enzymatic activity is a function of genomic ssDNA.

eRNA co-expression network uncovers TF dependency and convergent cooperativity

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

Hypoxic regulation of hypoxia inducible factor 1 alpha via antisense transcription

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

Characterizing intergenic transcription at RNA polymerase II binding sites in normal and cancer tissues

Intergenic transcription in normal and cancerous tissues is pervasive but incompletely understood. To investigate this, we constructed an atlas of over 180,000 consensus RNA polymerase II (RNAPII)-bound intergenic regions from 900 RNAPII chromatin immunoprecipitation sequencing (ChIP-seq) experiments in normal and cancer samples. Through unsupervised analysis, we identified 51 RNAPII consensus clusters, many of which mapped to specific biotypes and revealed tissue-specific regulatory signatures. We developed a meta-clustering methodology to integrate our RNAPII atlas with active transcription across 28,797 RNA sequencing (RNA-seq) samples from The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx), and Encyclopedia of DNA Elements (ENCODE). This analysis revealed strong tissue- and disease-specific interconnections between RNAPII occupancy and transcriptional activity. We demonstrate that intergenic transcription at RNAPII-bound regions is a novel per-cancer and pan-cancer biomarker. This biomarker displays genomic and clinically relevant characteristics, distinguishing cancer subtypes and linking to overall survival. Our results demonstrate the effectiveness of coherent data integration to uncover intergenic transcriptional activity in normal and cancer tissues.

Epromoters are new players in the regulatory landscape with potential pleiotropic roles

Precise spatiotemporal control of gene expression during normal development and cell differentiation is achieved by the combined action of proximal (promoters) and distal (enhancers) cis-regulatory elements. Recent studies have reported that a subset of promoters, termed Epromoters, works also as enhancers to regulate distal genes. This new paradigm opened novel questions regarding the complexity of our genome and raises the possibility that genetic variation within Epromoters has pleiotropic effects on various physiological and pathological traits by differentially impacting multiple proximal and distal genes. Here, we discuss the different observations pointing to an important role of Epromoters in the regulatory landscape and summarize the evidence supporting a pleiotropic impact of these elements in disease. We further hypothesize that Epromoter might represent a major contributor to phenotypic variation and disease.

Enhancers associated with unstable RNAs are rare in plants

Unstable transcripts have emerged as markers of active enhancers in vertebrates and shown to be involved in many cellular processes and medical disorders. However, their prevalence and role in plants is largely unexplored. Here, we comprehensively captured all actively initiating ("nascent") transcripts across diverse crops and other plants using capped small (cs)RNA-seq. We discovered that unstable transcripts are rare, unlike in vertebrates, and often originate from promoters. Additionally, many "distal" elements in plants initiate tissue-specific stable transcripts and are likely bone fide promoters of yet-unannotated genes or non-coding RNAs, cautioning against using genome annotations to infer "enhancers" or transcript stability. To investigate enhancer function, we integrated STARR-seq data. We found that annotated promoters, and other regions that initiate stable transcripts rather than unstable transcripts, function as stronger enhancers in plants. Our findings underscore the blurred line between promoters and enhancers and suggest that cis-regulatory elements encompass diverse structures and mechanisms in eukaryotes.

Single-cell nascent RNA sequencing using click-chemistry unveils coordinated transcription

Transcription is the primary regulatory step in gene expression. Divergent transcription initiation from promoters and enhancers produces stable RNAs from genes and unstable RNAs from enhancers1-5. Nascent RNA capture and sequencing assays simultaneously measure gene and enhancer activity in cell populations6-9. However, fundamental questions in the temporal regulation of transcription and enhancer-gene synchrony remain unanswered primarily due to the absence of a single-cell perspective on active transcription. In this study, we present scGRO-seq - a novel single-cell nascent RNA sequencing assay using click-chemistry - and unveil the coordinated transcription throughout the genome. scGRO-seq demonstrates the episodic nature of transcription, and estimates burst size and frequency by directly quantifying transcribing RNA polymerases in individual cells. It reveals the co-transcription of functionally related genes and leverages the replication-dependent non-polyadenylated histone genes transcription to elucidate cell-cycle dynamics. The single-nucleotide spatial and temporal resolution of scGRO-seq identifies networks of enhancers and genes and indicates that the bursting of transcription at super-enhancers precedes the burst from associated genes. By imparting insights into the dynamic nature of transcription and the origin and propagation of transcription signals, scGRO-seq demonstrates its unique ability to investigate the mechanisms of transcription regulation and the role of enhancers in gene expression.

KMT2A associates with PHF5A-PHF14-HMG20A-RAI1 subcomplex in pancreatic cancer stem cells and epigenetically regulates their characteristics

Pancreatic cancer (PC), one of the most aggressive and life-threatening human malignancies, is known for its resistance to cytotoxic therapies. This is increasingly ascribed to the subpopulation of undifferentiated cells, known as pancreatic cancer stem cells (PCSCs), which display greater evolutionary fitness than other tumor cells to evade the cytotoxic effects of chemotherapy. PCSCs are crucial for tumor relapse as they possess 'stem cell-like' features that are characterized by self-renewal and differentiation. However, the molecular mechanisms that maintain the unique characteristics of PCSCs are poorly understood. Here, we identify the histone methyltransferase KMT2A as a physical binding partner of an RNA polymerase-associated PHF5A-PHF14-HMG20A-RAI1 protein subcomplex and an epigenetic regulator of PCSC properties and functions. Targeting the protein subcomplex in PCSCs with a KMT2A-WDR5 inhibitor attenuates their self-renewal capacity, cell viability, and in vivo tumorigenicity.

Delineating highly transcribed noncoding elements landscape in breast cancer

Highly transcribed noncoding elements (HTNEs) are critical noncoding elements with high levels of transcriptional capacity in particular cohorts involved in multiple cellular biological processes. Investigation of HTNEs with persistent aberrant expression in abnormal tissues could be of benefit in exploring their roles in disease occurrence and progression. Breast cancer is a highly heterogeneous disease for which early screening and prognosis are exceedingly crucial. In this study, we developed a HTNE identification framework to systematically investigate HTNE landscapes in breast cancer patients and identified over ten thousand HTNEs. The robustness and rationality of our framework were demonstrated via public datasets. We revealed that HTNEs had significant chromatin characteristics of enhancers and long noncoding RNAs (lncRNAs) and were significantly enriched with RNA-binding proteins as well as targeted by miRNAs. Further, HTNE-associated genes were significantly overexpressed and exhibited strong correlations with breast cancer. Ultimately, we explored the subtype-specific transcriptional processes associated with HTNEs and uncovered the HTNE signatures that could classify breast cancer subtypes based on the properties of hormone receptors. Our results highlight that the identified HTNEs as well as their associated genes play crucial roles in breast cancer progression and correlate with subtype-specific transcriptional processes of breast cancer.

Combining TSS-MPRA and sensitive TSS profile dissimilarity scoring to study the sequence determinants of transcription initiation

Cis-regulatory elements (CREs) can be classified by the shapes of their transcription start site (TSS) profiles, which are indicative of distinct regulatory mechanisms. Massively parallel reporter assays (MPRAs) are increasingly being used to study CRE regulatory mechanisms, yet the degree to which MPRAs replicate individual endogenous TSS profiles has not been determined. Here, we present a new low-input MPRA protocol (TSS-MPRA) that enables measuring TSS profiles of episomal reporters as well as after lentiviral reporter chromatinization. To sensitively compare MPRA and endogenous TSS profiles, we developed a novel dissimilarity scoring algorithm (WIP score) that outperforms the frequently used earth mover's distance on experimental data. Using TSS-MPRA and WIP scoring on 500 unique reporter inserts, we found that short (153 bp) MPRA promoter inserts replicate the endogenous TSS patterns of ∼60% of promoters. Lentiviral reporter chromatinization did not improve fidelity of TSS-MPRA initiation patterns, and increasing insert size frequently led to activation of extraneous TSS in the MPRA that are not active in vivo. We discuss the implications of our findings, which highlight important caveats when using MPRAs to study transcription mechanisms. Finally, we illustrate how TSS-MPRA and WIP scoring can provide novel insights into the impact of transcription factor motif mutations and genetic variants on TSS patterns and transcription levels.

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

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

Enhancer in cancer pathogenesis and treatment

Enhancers are essential cis-acting regulatory elements that determine cell identity and tumor progression. Enhancer function is dependent on the physical interaction between the enhancer and its target promoter inside its local chromatin environment. Enhancer reprogramming is an important mechanism in cancer pathogenesis and can be driven by both cis and trans factors. Super enhancers are acquired at oncogenes in numerous cancer types and represent potential targets for cancer treatment. BET and CDK inhibitors act through mechanisms of enhancer function and have shown promising results in therapy for various types of cancer. Genome editing is another way to reprogram enhancers in cancer treatment. The relationship between enhancers and cancer has been revised by several authors in the past few years, which mainly focuses on the mechanisms by which enhancers can impact cancer. Here, we emphasize SE's role in cancer pathogenesis and the new therapies involving epigenetic regulators (BETi and CDKi). We suggest that understanding mechanisms of activity would aid clinical success for these anti-cancer agents.

Enhancer RNAs: mechanisms in transcriptional regulation and functions in diseases

In recent years, increasingly more non-coding RNAs have been detected with the development of high-throughput sequencing technology, including microRNAs (miRNAs), long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), small nucleolar RNAs (snoRNAs), and piwi-interacting RNA (piRNAs). The discovery of enhancer RNAs (eRNAs) in 2010 has further broadened the range of non-coding RNAs revealed. eRNAs are non-coding RNA molecules produced by the transcription of DNA cis-acting elements, enhancer fragments. Recent studies revealed that the transcription of eRNAs may be a biological marker responding to enhancer activity that can participate in the regulation of coding gene transcription. In this review, we discussed the biological characteristics of eRNAs, their functions in transcriptional regulation, the regulation factors of eRNAs production, and the research progress of eRNAs in different diseases. Video Abstract.

Evaluating the role of the nuclear microenvironment in gene function by population-based modeling

The nuclear folding of chromosomes relative to nuclear bodies is an integral part of gene function. Here, we demonstrate that population-based modeling-from ensemble Hi-C data-provides a detailed description of the nuclear microenvironment of genes and its role in gene function. We define the microenvironment by the subnuclear positions of genomic regions with respect to nuclear bodies, local chromatin compaction, and preferences in chromatin compartmentalization. These structural descriptors are determined in single-cell models, thereby revealing the structural variability between cells. We demonstrate that the microenvironment of a genomic region is linked to its functional potential in gene transcription, replication, and chromatin compartmentalization. Some chromatin regions feature a strong preference for a single microenvironment, due to association with specific nuclear bodies in most cells. Other chromatin shows high structural variability, which is a strong indicator of functional heterogeneity. Moreover, we identify specialized nuclear microenvironments, which distinguish chromatin in different functional states and reveal a key role of nuclear speckles in chromosome organization. We demonstrate that our method produces highly predictive three-dimensional genome structures, which accurately reproduce data from a variety of orthogonal experiments, thus considerably expanding the range of Hi-C data analysis.

RNA polymerase II dynamics shape enhancer-promoter interactions

How enhancers control target gene expression over long genomic distances remains an important unsolved problem. Here we investigated enhancer-promoter communication by integrating data from nucleosome-resolution genomic contact maps, nascent transcription and perturbations affecting either RNA polymerase II (Pol II) dynamics or the activity of thousands of candidate enhancers. Integration of new Micro-C experiments with published CRISPRi data demonstrated that enhancers spend more time in close proximity to their target promoters in functional enhancer-promoter pairs compared to nonfunctional pairs, which can be attributed in part to factors unrelated to genomic position. Manipulation of the transcription cycle demonstrated a key role for Pol II in enhancer-promoter interactions. Notably, promoter-proximal paused Pol II itself partially stabilized interactions. We propose an updated model in which elements of transcriptional dynamics shape the duration or frequency of interactions to facilitate enhancer-promoter communication.

Activation, decommissioning, and dememorization: enhancers in a life cycle

Spatiotemporal regulation of cell type-specific gene expression is essential to convert a zygote into a complex organism that contains hundreds of distinct cell types. A class of cis-regulatory elements called enhancers, which have the potential to enhance target gene transcription, are crucial for precise gene expression programs during development. Following decades of research, many enhancers have been discovered and how enhancers become activated has been extensively studied. However, the mechanisms underlying enhancer silencing are less well understood. We review current understanding of enhancer decommissioning and dememorization, both of which enable enhancer silencing. We highlight recent progress from genome-wide perspectives that have revealed the life cycle of enhancers and how its dynamic regulation underlies cell fate transition, development, cell regeneration, and epigenetic reprogramming.

Toward a comprehensive catalog of regulatory elements

Regulatory elements are the genomic regions that interact with transcription factors to control cell-type-specific gene expression in different cellular environments. A precise and complete catalog of functional elements encoded by the human genome is key to understanding mammalian gene regulation. Here, we review the current state of regulatory element annotation. We first provide an overview of assays for characterizing functional elements, including genome, epigenome, transcriptome, three-dimensional chromatin interaction, and functional validation assays. We then discuss computational methods for defining regulatory elements, including peak-calling and other statistical modeling methods. Finally, we introduce several high-quality lists of regulatory element annotations and suggest potential future directions.

Prioritizing genes associated with brain disorders by leveraging enhancer-promoter interactions in diverse neural cells and tissues

BACKGROUND

Prioritizing genes that underlie complex brain disorders poses a considerable challenge. Despite previous studies have found that they shared symptoms and heterogeneity, it remained difficult to systematically identify the risk genes associated with them.

METHODS

By using the CAGE (Cap Analysis of Gene Expression) read alignment files for 439 human cell and tissue types (including primary cells, tissues and cell lines) from FANTOM5 project, we predicted enhancer-promoter interactions (EPIs) of 439 cell and tissue types in human, and examined their reliability. Then we evaluated the genetic heritability of 17 diverse brain disorders and behavioral-cognitive phenotypes in each neural cell type, brain region, and developmental stage. Furthermore, we prioritized genes associated with brain disorders and phenotypes by leveraging the EPIs in each neural cell and tissue type, and analyzed their pleiotropy and functionality for different categories of disorders and phenotypes. Finally, we characterized the spatiotemporal expression dynamics of these associated genes in cells and tissues.

RESULTS

We found that identified EPIs showed activity specificity and network aggregation in cell and tissue types, and enriched TF binding in neural cells played key roles in synaptic plasticity and nerve cell development, i.e., EGR1 and SOX family. We also discovered that most neurological disorders exhibit heritability enrichment in neural stem cells and astrocytes, while psychiatric disorders and behavioral-cognitive phenotypes exhibit enrichment in neurons. Furthermore, our identified genes recapitulated well-known risk genes, which exhibited widespread pleiotropy between psychiatric disorders and behavioral-cognitive phenotypes (i.e., FOXP2), and indicated expression specificity in neural cell types, brain regions, and developmental stages associated with disorders and phenotypes. Importantly, we showed the potential associations of brain disorders with brain regions and developmental stages that have not been well studied.

CONCLUSIONS

Overall, our study characterized the gene-enhancer regulatory networks and genetic mechanisms in the human neural cells and tissues, and illustrated the value of reanalysis of publicly available genomic datasets.

JUN-induced super-enhancer RNA forms R-loop to promote nasopharyngeal carcinoma metastasis

Oncogenic super-enhancers (SEs) generate noncoding enhancer/SE RNAs (eRNAs/seRNAs) that exert a critical function in malignancy through powerful regulation of target gene expression. Herein, we show that a JUN-mediated seRNA can form R-loop to regulate target genes to promote metastasis of nasopharyngeal carcinoma (NPC). A combination of global run-on sequencing, chromatin-immunoprecipitation sequencing, and RNA sequencing was used to screen seRNAs. A specific seRNA associated with NPC metastasis (seRNA-NPCM) was identified as a transcriptional regulator for N-myc downstream-regulated gene 1 (NDRG1). JUN was found to regulate seRNA-NPCM through motif binding. seRNA-NPCM was elevated in NPC cancer tissues and highly metastatic cell lines, and promoted the metastasis of NPC cells in vitro and in vivo. Mechanistically, the 3' end of seRNA-NPCM hybridizes with the SE region to form an R-loop, and the middle segment of seRNA-NPCM binds to heterogeneous nuclear ribonucleoprotein R (hnRNPR) at the promoter of distal gene NDRG1 and neighboring gene tribbles pseudokinase 1 (TRIB1). These structures promote chromatin looping and long-distance chromatin interactions between SEs and promoters, thus facilitating NDRG1 and TRIB1 transcription. Furthermore, the clinical analyses showed that seRNA-NPCM and NDRG1 were independent prognostic factors for NPC patients. seRNA-NPCM plays a critical role in orchestrating target gene transcription to promote NPC metastasis.

Sequence basis of transcription initiation in human genome

Transcription initiation is an essential process for ensuring proper function of any gene, however, a unified understanding of sequence patterns and rules that determine transcription initiation sites in human genome remains elusive. By explaining transcription initiation at basepair resolution from sequence with a deep learning-inspired explainable modeling approach, here we show that simple rules can explain the vast majority of human promoters. We identified key sequence patterns that contribute to human promoter function, each activating transcription with a distinct position-specific effect curve that likely reflects its mechanism of promoting transcription initiation. Most of these position-specific effects have not been previously characterized, and we verified them using experimental perturbations of transcription factors and sequences. We revealed the sequence basis of bidirectional transcription at promoters and links between promoter selectivity and gene expression variation across cell types. Additionally, by analyzing 241 mammalian genomes and mouse transcription initiation site data, we showed that the sequence determinants are conserved across mammalian species. Taken together, we provide a unified model of the sequence basis of transcription initiation at the basepair level that is broadly applicable across mammalian species, and shed new light on basic questions related to promoter sequence and function.

Defining the fine structure of promoter activity on a genome-wide scale with CISSECTOR

Classic promoter mutagenesis strategies can be used to study how proximal promoter regions regulate the expression of particular genes of interest. This is a laborious process, in which the smallest sub-region of the promoter still capable of recapitulating expression in an ectopic setting is first identified, followed by targeted mutation of putative transcription factor binding sites. Massively parallel reporter assays such as survey of regulatory elements (SuRE) provide an alternative way to study millions of promoter fragments in parallel. Here we show how a generalized linear model (GLM) can be used to transform genome-scale SuRE data into a high-resolution genomic track that quantifies the contribution of local sequence to promoter activity. This coefficient track helps identify regulatory elements and can be used to predict promoter activity of any sub-region in the genome. It thus allows in silico dissection of any promoter in the human genome to be performed. We developed a web application, available at cissector.nki.nl, that lets researchers easily perform this analysis as a starting point for their research into any promoter of interest.

Chromatin alternates between A and B compartments at kilobase scale for subgenic organization

Nuclear compartments are prominent features of 3D chromatin organization, but sequencing depth limitations have impeded investigation at ultra fine-scale. CTCF loops are generally studied at a finer scale, but the impact of looping on proximal interactions remains enigmatic. Here, we critically examine nuclear compartments and CTCF loop-proximal interactions using a combination of in situ Hi-C at unparalleled depth, algorithm development, and biophysical modeling. Producing a large Hi-C map with 33 billion contacts in conjunction with an algorithm for performing principal component analysis on sparse, super massive matrices (POSSUMM), we resolve compartments to 500 bp. Our results demonstrate that essentially all active promoters and distal enhancers localize in the A compartment, even when flanking sequences do not. Furthermore, we find that the TSS and TTS of paused genes are often segregated into separate compartments. We then identify diffuse interactions that radiate from CTCF loop anchors, which correlate with strong enhancer-promoter interactions and proximal transcription. We also find that these diffuse interactions depend on CTCF's RNA binding domains. In this work, we demonstrate features of fine-scale chromatin organization consistent with a revised model in which compartments are more precise than commonly thought while CTCF loops are more protracted.

Functional analysis of a random-sequence chromosome reveals a high level and the molecular nature of transcriptional noise in yeast cells

We measure transcriptional noise in yeast by analyzing chromatin structure and transcription of an 18-kb region of DNA whose sequence was randomly generated. Nucleosomes fully occupy random-sequence DNA, but nucleosome-depleted regions (NDRs) are much less frequent, and there are fewer well-positioned nucleosomes and shorter nucleosome arrays. Steady-state levels of random-sequence RNAs are comparable to yeast mRNAs, although transcription and decay rates are higher. Transcriptional initiation from random-sequence DNA occurs at numerous sites, indicating very low intrinsic specificity of the RNA Pol II machinery. In contrast, poly(A) profiles of random-sequence RNAs are roughly comparable to those of yeast mRNAs, suggesting limited evolutionary restraints on poly(A) site choice. Random-sequence RNAs show higher cell-to-cell variability than yeast mRNAs, suggesting that functional elements limit variability. These observations indicate that transcriptional noise occurs at high levels in yeast, and they provide insight into how chromatin and transcription patterns arise from the evolved yeast genome.

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

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

Decoding enhancer complexity with machine learning and high-throughput discovery

Enhancers are genomic DNA elements controlling spatiotemporal gene expression. Their flexible organization and functional redundancies make deciphering their sequence-function relationships challenging. This article provides an overview of the current understanding of enhancer organization and evolution, with an emphasis on factors that influence these relationships. Technological advancements, particularly in machine learning and synthetic biology, are discussed in light of how they provide new ways to understand this complexity. Exciting opportunities lie ahead as we continue to unravel the intricacies of enhancer function.

Mistimed origin licensing and activation stabilize common fragile sites under tight DNA-replication checkpoint activation

Genome integrity requires replication to be completed before chromosome segregation. The DNA-replication checkpoint (DRC) contributes to this coordination by inhibiting CDK1, which delays mitotic onset. Under-replication of common fragile sites (CFSs), however, escapes surveillance, resulting in mitotic chromosome breaks. Here we asked whether loose DRC activation induced by modest stresses commonly used to destabilize CFSs could explain this leakage. We found that tightening DRC activation or CDK1 inhibition stabilizes CFSs in human cells. Repli-Seq and molecular combing analyses showed a burst of replication initiations implemented in mid S-phase across a subset of late-replicating sequences, including CFSs, while the bulk genome was unaffected. CFS rescue and extra-initiations required CDC6 and CDT1 availability in S-phase, implying that CDK1 inhibition permits mistimed origin licensing and firing. In addition to delaying mitotic onset, tight DRC activation therefore supports replication completion of late origin-poor domains at risk of under-replication, two complementary roles preserving genome stability.

Histone Deacetylases (HDACs) maintain expression of the pluripotent gene network via recruitment of RNA polymerase II to coding and non-coding loci

Histone acetylation is a dynamic modification regulated by the opposing actions of histone acetyltransferases (HATs) and histone deacetylases (HDACs). Deacetylation of histone tails results in chromatin tightening and therefore HDACs are generally regarded as transcriptional repressors. Counterintuitively, simultaneous deletion of Hdac1 and Hdac2 in embryonic stem cells (ESC) reduced expression of pluripotent transcription factors, Oct4, Sox2 and Nanog (OSN). By shaping global histone acetylation patterns, HDACs indirectly regulate the activity of acetyl-lysine readers, such as the transcriptional activator, BRD4. We used inhibitors of HDACs and BRD4 (LBH589 and JQ1 respectively) in combination with precision nuclear run-on and sequencing (PRO-seq) to examine their roles in defining the ESC transcriptome. Both LBH589 and JQ1 caused a marked reduction in the pluripotent network. However, while JQ1 treatment induced widespread transcriptional pausing, HDAC inhibition caused a reduction in both paused and elongating polymerase, suggesting an overall reduction in polymerase recruitment. Using enhancer RNA (eRNA) expression to measure enhancer activity we found that LBH589-sensitive eRNAs were preferentially associated with super-enhancers and OSN binding sites. These findings suggest that HDAC activity is required to maintain pluripotency by regulating the OSN enhancer network via the recruitment of RNA polymerase II.

Considerations and caveats for analyzing chromatin compartments

Genomes are organized into nuclear compartments, separating active from inactive chromatin. Chromatin compartments are readily visible in a large number of species by experiments that map chromatin conformation genome-wide. When analyzing these maps, a common step is the identification of genomic intervals that interact within A (active) and B (inactive) compartments. It has also become increasingly common to identify and analyze subcompartments. We review different strategies to identify A/B and subcompartment intervals, including a discussion of various machine-learning approaches to predict these features. We then discuss the strengths and limitations of current strategies and examine how these aspects of analysis may have impacted our understanding of chromatin compartments.

Insights into Enhancer RNAs: Biogenesis and Emerging Role in Brain Diseases

Enhancers are cis-acting elements that control the transcription of target genes and are transcribed into a class of noncoding RNAs (ncRNAs) termed enhancer RNAs (eRNAs). eRNAs have shorter half-lives than mRNAs and long noncoding RNAs; however, the frequency of transcription of eRNAs is close to that of mRNAs. eRNA expression is associated with a high level of histone mark H3K27ac and a low level of H3K27me3. Although eRNAs only account for a small proportion of ncRNAs, their functions are important. eRNAs can not only increase enhancer activity by promoting the formation of enhancer-promoter loops but also regulate transcriptional activation. Increasing numbers of studies have found that eRNAs play an important role in the occurrence and development of brain diseases; however, further research into eRNAs is required. This review discusses the concept, characteristics, classification, function, and potential roles of eRNAs in brain diseases.

Revealing eRNA interactions: TF dependency and convergent cooperativity

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

Profiling Chromatin Accessibility on Replicated DNA with repli-ATAC-Seq

Open or accessible chromatin typifies euchromatic regions and helps define cell type-specific transcription programs. DNA replication massively disorders chromatin composition and structure, and how accessible regions are affected by and recover from this disruption has been unclear. Here, we present repli-ATAC-seq, a protocol to profile accessible chromatin genome-wide on replicated DNA starting from 100,000 cells. In this method, replicated DNA is labeled with a short 5-ethynyl-2'-deoxyuridine (EdU) pulse in cultured cells and isolated from a population of tagmented fragments for amplification and next-generation sequencing. Repli-ATAC-seq provides high-resolution information on chromatin dynamics after DNA replication and reveals new insights into the interplay between DNA replication, transcription, and the chromatin landscape.

Dynamic activity in cis-regulatory elements of leukocytes identifies transcription factor activation and stratifies COVID-19 severity in ICU patients

Transcription factor programs mediating the immune response to coronavirus disease 2019 (COVID-19) are not fully understood. Capturing active transcription initiation from cis-regulatory elements such as enhancers and promoters by capped small RNA sequencing (csRNA-seq), in contrast to capturing steady-state transcripts by conventional RNA-seq, allows unbiased identification of the underlying transcription factor activity and regulatory pathways. Here, we profile transcription initiation in critically ill COVID-19 patients, identifying transcription factor motifs that correlate with clinical lung injury and disease severity. Unbiased clustering reveals distinct subsets of cis-regulatory elements that delineate the cell type, pathway-specific, and combinatorial transcription factor activity. We find evidence of critical roles of regulatory networks, showing that STAT/BCL6 and E2F/MYB regulatory programs from myeloid cell populations are activated in patients with poor disease outcomes and associated with COVID-19 susceptibility genetic variants. More broadly, we demonstrate how capturing acute, disease-mediated changes in transcription initiation can provide insight into the underlying molecular mechanisms and stratify patient disease severity.

PRO-IP-seq Tracks Molecular Modifications of Engaged Pol II Complexes at Nucleotide Resolution

RNA Polymerase II (Pol II) is a multi-subunit complex that undergoes covalent modifications as transcription proceeds through genes and enhancers. Rate-limiting steps of transcription control Pol II recruitment, site and degree of initiation, pausing duration, productive elongation, nascent transcript processing, transcription termination, and Pol II recycling. Here, we developed Precision Run-On coupled to Immuno-Precipitation sequencing (PRO-IP-seq) and tracked phosphorylation of Pol II C-terminal domain (CTD) at nucleotide-resolution. We uncovered precise positional control of Pol II CTD phosphorylation as transcription proceeds from the initiating nucleotide, through early and late promoter-proximal pause, and into productive elongation. Pol II CTD was predominantly unphosphorylated in the early pause-region, whereas serine-2- and serine-5-phosphorylations occurred preferentially in the later pause-region. Serine-7-phosphorylation dominated after the pause-release in a region where Pol II accelerates to its full elongational speed. Interestingly, tracking transcription upon heat-induced reprogramming demonstrated that Pol II with phosphorylated CTD remains paused on heat-repressed genes.

ATM suppresses c-Myc overexpression in the mammary epithelium in response to estrogen

ATM gene mutation carriers are predisposed to estrogen-receptor-positive breast cancer (BC). ATM prevents BC oncogenesis by activating p53 in every cell; however, much remains unknown about tissue-specific oncogenesis after ATM loss. Here, we report that ATM controls the early transcriptional response to estrogens. This response depends on topoisomerase II (TOP2), which generates TOP2-DNA double-strand break (DSB) complexes and rejoins the breaks. When TOP2-mediated ligation fails, ATM facilitates DSB repair. After estrogen exposure, TOP2-dependent DSBs arise at the c-MYC enhancer in human BC cells, and their defective repair changes the activation profile of enhancers and induces the overexpression of many genes, including the c-MYC oncogene. CRISPR/Cas9 cleavage at the enhancer also causes c-MYC overexpression, indicating that this DSB causes c-MYC overexpression. Estrogen treatment induced c-Myc protein overexpression in mammary epithelial cells of ATM-deficient mice. In conclusion, ATM suppresses the c-Myc-driven proliferative effects of estrogens, possibly explaining such tissue-specific oncogenesis.

Specificity Guides Interpretation: On H3K4 Methylation at Enhancers and Broad Promoters

In 2018, we used internally calibrated chromatin immunoprecipitation (ICeChIP) to find that many of the most commonly used antibodies against H3K4 methylforms had significant off-target binding, which compromised the findings of at least eight literature paradigms that used these antibodies for ChIP-seq (Shah et al., 2018). In many cases, we were able to recapitulate the prior findings in K562 cells with the original, low-quality antibody, only to find that the models did not hold up to scrutiny with highly specific reagents and quantitative calibration. In a recent preprint originally prepared as a Letter to the Editor of Molecular Cell, though they agree with our overarching conclusions, Pekowska and colleagues take issue with analyses presented for two relatively minor points of the paper (Pekowska et al., 2023). We are puzzled by the assertion that these two points constitute the "bulk" of our findings, nor is it clear which components of our "analytical design" they find problematic. We feel their critique, however mild, is misguided.

The regulation of transcription elongation in embryonic stem cells

Transcription elongation is a fundamental molecular process which is accurately regulated to ensure proper gene expression in cellular activities whereas its malfunction is associated with impaired cellular functions. Embryonic stem cells (ESCs) have significant value in regenerative medicine due to their self-renewal ability and their potential to differentiate to almost all types of cells. Therefore, dissection of the exact regulatory mechanism of transcription elongation in ESCs is crucial for both basic research and their clinical applications. In this review, we discuss the current understanding on the regulatory mechanisms of transcription elongation mediated by transcription factors and epigenetic modifications in ESCs.

Transcription factor antagonism regulates heterogeneity in embryonic stem cell states

Gene expression heterogeneity underlies cell states and contributes to developmental robustness. While heterogeneity can arise from stochastic transcriptional processes, the extent to which it is regulated is unclear. Here, we characterize the regulatory program underlying heterogeneity in murine embryonic stem cell (mESC) states. We identify differentially active and transcribed enhancers (DATEs) across states. DATEs regulate differentially expressed genes and are distinguished by co-binding of transcription factors Klf4 and Zfp281. In contrast to other factors that interact in a positive feedback network stabilizing mESC cell-type identity, Klf4 and Zfp281 drive opposing transcriptional and chromatin programs. Abrogation of factor binding to DATEs dampens variation in gene expression, and factor loss alters kinetics of switching between states. These results show antagonism between factors at enhancers results in gene expression heterogeneity and formation of cell states, with implications for the generation of diverse cell types during development.

MEPP: more transparent motif enrichment by profiling positional correlations

Score-based motif enrichment analysis (MEA) is typically applied to regulatory DNA to infer transcription factors (TFs) that may modulate transcription and chromatin state in different conditions. Most MEA methods determine motif enrichment independent of motif position within a sequence, even when those sequences harbor anchor points that motifs and their bound TFs may functionally interact with in a distance-dependent fashion, such as other TF binding motifs, transcription start sites (TSS), sequencing assay cleavage sites, or other biologically meaningful features. We developed motif enrichment positional profiling (MEPP), a novel MEA method that outputs a positional enrichment profile of a given TF's binding motif relative to key anchor points (e.g. transcription start sites, or other motifs) within the analyzed sequences while accounting for lower-order nucleotide bias. Using transcription initiation and TF binding as test cases, we demonstrate MEPP's utility in determining the sequence positions where motif presence correlates with measures of biological activity, inferring positional dependencies of binding site function. We demonstrate how MEPP can be applied to interpretation and hypothesis generation from experiments that quantify transcription initiation, chromatin structure, or TF binding measurements. MEPP is available for download from https://github.com/npdeloss/mepp.

A common transcriptional mechanism involving R-loop and RNA abasic site regulates an enhancer RNA of APOE

RNA is modified by hundreds of chemical reactions and folds into innumerable shapes. However, the regulatory role of RNA sequence and structure and how dysregulation leads to diseases remain largely unknown. Here, we uncovered a mechanism where RNA abasic sites in R-loops regulate transcription by pausing RNA polymerase II. We found an enhancer RNA, AANCR, that regulates the transcription and expression of apolipoprotein E (APOE). In some human cells such as fibroblasts, AANCR is folded into an R-loop and modified by N-glycosidic cleavage; in this form, AANCR is a partially transcribed nonfunctional enhancer and APOE is not expressed. In contrast, in other cell types including hepatocytes and under stress, AANCR does not form a stable R-loop as its sequence is not modified, so it is transcribed into a full-length enhancer that promotes APOE expression. DNA sequence variants in AANCR are associated significantly with APOE expression and Alzheimer's Disease, thus AANCR is a modifier of Alzheimer's Disease. Besides AANCR, thousands of noncoding RNAs are regulated by abasic sites in R-loops. Together our data reveal the essentiality of the folding and modification of RNA in cellular regulation and demonstrate that dysregulation underlies common complex diseases such as Alzheimer's disease.

Identification of transcription factors dictating blood cell development using a bidirectional transcription network-based computational framework

Advanced computational methods exploit gene expression and epigenetic datasets to predict gene regulatory networks controlled by transcription factors (TFs). These methods have identified cell fate determining TFs but require large amounts of reference data and experimental expertise. Here, we present an easy to use network-based computational framework that exploits enhancers defined by bidirectional transcription, using as sole input CAGE sequencing data to correctly predict TFs key to various human cell types. Next, we applied this Analysis Algorithm for Networks Specified by Enhancers based on CAGE (ANANSE-CAGE) to predict TFs driving red and white blood cell development, and THP-1 leukemia cell immortalization. Further, we predicted TFs that are differentially important to either cell line- or primary- associated MLL-AF9-driven gene programs, and in primary MLL-AF9 acute leukemia. Our approach identified experimentally validated as well as thus far unexplored TFs in these processes. ANANSE-CAGE will be useful to identify transcription factors that are key to any cell fate change using only CAGE-seq data as input.