Human promoters are intrinsically directional

Mediator kinase inhibition drives myometrial stem cell differentiation and the uterine fibroid phenotype through super-enhancer reprogramming

Uterine fibroids (UFs) are the most common non-cutaneous tumors in women worldwide. UFs arise from genetic alterations in myometrial stem cells (MM SCs) that trigger their transformation into tumor-initiating cells (UF SCs). Mutations in the RNA polymerase II Mediator subunit MED12 are dominant drivers of UFs, accounting for 70% of these clinically significant lesions. Biochemically, UF driver mutations in MED12 disrupt CDK8/19 kinase activity in Mediator, but how Mediator kinase disruption triggers MM SC transformation remains unknown. Here, we show that pharmacologic inhibition of CDK8/19 in MM SCs removes a barrier to myogenic differentiation down an altered pathway characterized by molecular phenotypes characteristic of UFs, including oncogenic growth and extracellular matrix (ECM) production. These perturbations appear to be induced by transcriptomic changes, arising in part through epigenomic alteration and super-enhancer reprogramming, that broadly recapitulate those found in MED12-mutant UFs. Altogether, these findings provide new insights concerning the biological role of CDK8/19 in MM SC biology and UF formation. KEY MESSAGES: Mediator kinase inhibition in myometrial stem cells (MM SCs) induces spontaneous differentiation. Transcriptional changes upon Mediator kinase inhibition recapitulate those of MED12 mutant uterine fibroids (UFs). Such transcriptional changes are partially mediated by super-enhancer reprogramming. Mediator kinase functions to enforce cell states and its loss induces cellular plasticity.

Novel insights into RNA polymerase II transcription regulation: transcription factors, phase separation, and their roles in cardiovascular diseases

Transcription factors (TFs) are specialized proteins that bind DNA in a sequence-specific manner and modulate RNA polymerase II (Pol II) in multiple steps of the transcription process. Phase separation is a spontaneous or driven process that can form membrane-less organelles called condensates. By creating different liquid phases at active transcription sites, the formation of transcription condensates can reduce the water content of the condensate and lower the dielectric constant in biological systems, which in turn alters the structure and function of proteins and nucleic acids in the condensate. In RNA Pol II transcription, phase separation formation shortens the time at which TFs bind to target DNA sites and promotes transcriptional bursting. RNA Pol II transcription is engaged in developing several diseases, such as cardiovascular disease, by regulating different TFs and mediating the occurrence of phase separation. This review aims to summarize the advances in the molecular mechanisms of RNA Pol II transcriptional regulation, in particular the effect of TFs and phase separation. The role of RNA Pol II transcriptional regulation in cardiovascular disease will be elucidated, providing potential therapeutic targets for the management and treatment of cardiovascular disease.

Mediator kinase inhibition drives myometrial stem cell differentiation and the uterine fibroid phenotype through super-enhancer reprogramming

Uterine fibroids (UFs) are the most common non-cutaneous tumors in women worldwide. UFs arise from genetic alterations in myometrial stem cells (MM SCs) that trigger their transformation into tumor initiating cells (UF SCs). Mutations in the RNA polymerase II Mediator subunit MED12 are dominant drivers of UFs, accounting for 70% of these clinically significant lesions. Biochemically, UF driver mutations in MED12 disrupt CDK8/19 kinase activity in Mediator, but how Mediator kinase disruption triggers MM SC transformation remains unknown. Here, we show that pharmacologic inhibition of CDK8/19 in MM SCs removes a barrier to myogenic differentiation down an altered pathway characterized by molecular phenotypes characteristic of UFs, including oncogenic growth and extracellular matrix (ECM) production. These perturbations appear to be induced by transcriptomic changes, arising in part through epigenomic alteration and super-enhancer reprogramming, that broadly recapitulate those found in MED12-mutant UFs. Altogether these findings provide new insights concerning the biological role of CDK8/19 in MM SC biology and UF formation.

SRSF2 safeguards efficient transcription of DNA damage and repair genes

The serine-/arginine-rich splicing factor 2 (SRSF2) plays pivotal roles in pre-mRNA processing and gene transcription. Recurrent mutations, particularly a proline-to-histidine substitution at position 95 (P95H), are common in neoplastic diseases. Here, we assess SRSF2's diverse functions in squamous cell carcinoma. We show that SRSF2 deletion or homozygous P95H mutation both cause extensive DNA damage leading to cell-cycle arrest. Mechanistically, SRSF2 regulates efficient bi-directional transcription of DNA replication and repair genes, independent from its function in splicing. Further, SRSF2 haploinsufficiency induces DNA damage without halting the cell cycle. Exposing mouse skin to tumor-promoting carcinogens enhances the clonal expansion of heterozygous Srsf2 P95H epidermal cells but unexpectedly inhibits tumor formation. To survive carcinogen treatment, Srsf2 P95H+/- cells undergo substantial transcriptional rewiring and restore bi-directional gene expression. Thus, our study underscores SRSF2's importance in regulating transcription to orchestrate the cell cycle and the DNA damage response.

Dissecting the cis-regulatory syntax of transcription initiation with deep learning

Despite extensive characterization of mammalian Pol II transcription, the DNA sequence determinants of transcription initiation at a third of human promoters and most enhancers remain poorly understood. We trained and interpreted a neural network called ProCapNet that accurately models base-resolution initiation profiles from PRO-cap experiments using local DNA sequence. ProCapNet learns sequence motifs with distinct effects on initiation rates and TSS positioning and uncovers context-specific cryptic initiator elements intertwined within other TF motifs. ProCapNet annotates predictive motifs in nearly all actively transcribed regulatory elements across multiple cell-lines, revealing a shared cis-regulatory logic across promoters and enhancers and a highly epistatic sequence syntax of cooperative and competitive motif interactions. ProCapNet models of steady-state RAMPAGE profiles distill initiation signals on par with models trained directly on PRO-cap profiles. ProCapNet learns a largely cell-type-agnostic cis-regulatory code of initiation complementing sequence drivers of cell-type-specific chromatin state critical for accurate prediction of cell-type-specific transcription initiation.

mosGraphGen: a novel tool to generate multi-omics signaling graphs to facilitate integrative and interpretable graph AI model development

Motivation

Multi-omics data, i.e. genomics, epigenomics, transcriptomics, proteomics, characterize cellular complex signaling systems from multi-level and multi-view and provide a holistic view of complex cellular signaling pathways. However, it remains challenging to integrate and interpret multi-omics data for mining critical biomarkers. Graph AI models have been widely used to analyze graph-structure datasets, and are ideal for integrative multi-omics data analysis because they can naturally integrate and represent multi-omics data as a biologically meaningful multi-level signaling graph and interpret multi-omics data via graph node and edge ranking analysis. Nevertheless, it is nontrivial for graph-AI model developers to pre-analyze multi-omics data and convert the data into biologically meaningful graphs, which can be directly fed into graph-AI models.

Results

To resolve this challenge, we developed mosGraphGen (multi-omics signaling graph generator), generating Multi-omics Signaling graphs (mos-graph) of individual samples by mapping multi-omics data onto a biologically meaningful multi-level background signaling network with data normalization by aggregating measurements and aligning to the reference genome. With mosGraphGen, AI model developers can directly apply and evaluate their models using these mos-graphs. In the results, mosGraphGen was used and illustrated using two widely used multi-omics datasets of The Cancer Genome Atlas (TCGA) and Alzheimer's disease (AD) samples.

Availability and implementation

The code of mosGraphGen is open-source and publicly available via GitHub: https://github.com/FuhaiLiAiLab/mosGraphGen.

mosGraphGen: a novel tool to generate multi-omics signaling graphs to facilitate integrative and interpretable graph AI model development

Multi-omics data, i.e., genomics, epigenomics, transcriptomics, proteomics, characterize cellular complex signaling systems from multi-level and multi-view and provide a holistic view of complex cellular signaling pathways. However, it remains challenging to integrate and interpret multi-omics data for mining key disease targets and signaling pathways. Graph AI models have been widely used to analyze graph-structure datasets, and are ideal for integrative multi-omics data analysis because they can naturally integrate and represent multi-omics data as a biologically meaningful multi-level signaling graph and interpret multi-omics data via graph node and edge ranking analysis. However, it is non-trivial for graph-AI model developers to pre-analyze multi-omics data and convert the data into biologically meaningful graphs, which can be directly fed into graph-AI models. To resolve this challenge, we developed mosGraphGen (multi-omics signaling graph generator), generating Multi-omics Signaling graphs (mos-graph) of individual samples by mapping multi-omics data onto a biologically meaningful multi-level background signaling network with data normalization by aggregating measurements and aligning to the reference genome. With mosGraphGen, AI model developers can directly apply and evaluate their models using these mos-graphs. In the results, mosGraphGen was used and illustrated using two widely used multi-omics datasets of TCGA and Alzheimer's disease (AD) samples. The code of mosGraphGen is open-source and publicly available via GitHub: https://github.com/FuhaiLiAiLab/mosGraphGen.

Transcription directionality is licensed by Integrator at active human promoters

A universal characteristic of eukaryotic transcription is that the promoter recruits RNA polymerase II (RNAPII) to produce both precursor mRNAs (pre-mRNAs) and short unstable promoter upstream transcripts (PROMPTs) toward the opposite direction. However, how the transcription machinery selects the correct direction to produce pre-mRNAs is largely unknown. Here, through multiple acute auxin-inducible degradation systems, we show that rapid depletion of an RNAPII-binding protein complex, Integrator, results in robust PROMPT accumulation throughout the genome. Interestingly, the accumulation of PROMPTs is compensated by the reduction of pre-mRNA transcripts in actively transcribed genes. Consistently, Integrator depletion alters the distribution of polymerase between the sense and antisense directions, which is marked by increased RNAPII-carboxy-terminal domain Tyr1 phosphorylation at PROMPT regions and a reduced Ser2 phosphorylation level at transcription start sites. Mechanistically, the endonuclease activity of Integrator is critical to suppress PROMPT production. Furthermore, our data indicate that the presence of U1 binding sites on nascent transcripts could counteract the cleavage activity of Integrator. In this process, the absence of robust U1 signal at most PROMPTs allows Integrator to suppress the antisense transcription and shift the transcriptional balance in favor of the sense direction.

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 sequencing. We discovered that unstable transcripts are rare in plants, unlike in vertebrates, and when present, often originate from promoters. In addition, many 'distal' elements in plants initiate tissue-specific stable transcripts and are likely bona fide promoters of as-yet-unannotated genes or non-coding RNAs, cautioning against using reference genome annotations to infer putative enhancer sites. To investigate enhancer function, we integrated data from self-transcribing active regulatory region (STARR) sequencing. We found that annotated promoters and other regions that initiate stable transcripts, but not those marked by unstable or bidirectional unstable transcripts, showed stronger enhancer activity in this assay. Our findings underscore the blurred line between promoters and enhancers and suggest that cis-regulatory elements can encompass diverse structures and mechanisms in eukaryotes, including humans.

Molecular models of bidirectional promoter regulation

Approximately 11% of human genes are transcribed by a bidirectional promoter (BDP), defined as two genes with <1 kb between their transcription start sites. Despite their evolutionary conservation and enrichment for housekeeping genes and oncogenes, the regulatory role of BDPs remains unclear. BDPs have been suggested to facilitate gene coregulation and/or decrease expression noise. This review discusses these potential regulatory functions through the context of six prospective underlying mechanistic models: a single nucleosome free region, shared transcription factor/regulator binding, cooperative negative supercoiling, bimodal histone marks, joint activation by enhancer(s), and RNA-mediated recruitment of regulators. These molecular mechanisms may act independently and/or cooperatively to facilitate the coregulation and/or decreased expression noise predicted of BDPs.

Characterizing nascent transcription patterns of PROMPTs, eRNAs, and readthrough transcripts in the ENCODE4 deeply profiled cell lines

Arising as co-products of canonical gene expression, transcription-associated lincRNAs, such as promoter upstream transcripts (PROMPTs), enhancer RNAs (eRNAs), and readthrough (RT) transcripts, are often regarded as byproducts of transcription, although they may be important for the expression of nearby genes. We identified regions of nascent expression of these lincRNA in 16 human cell lines using Bru-seq techniques, and found distinctly regulated patterns of PROMPT, eRNA, and RT transcription using the diverse biochemical approaches in the ENCODE4 deeply profiled cell lines collection. Transcription of these lincRNAs was influenced by sequence-specific features and the local or 3D chromatin landscape. However, these sequence and chromatin features do not describe the full spectrum of lincRNA expression variability we identify, highlighting the complexity of their regulation. This may suggest that transcription-associated lincRNAs are not merely byproducts, but rather that the transcript itself, or the act of its transcription, is important for genomic function.

Cis-regulatory effect of HPV integration is constrained by host chromatin architecture in cervical cancers

Human papillomavirus (HPV) infections are the primary drivers of cervical cancers, and often HPV DNA gets integrated into the host genome. Although the oncogenic impact of HPV encoded genes is relatively well known, the cis-regulatory effect of integrated HPV DNA on host chromatin structure and gene regulation remains less understood. We investigated genome-wide patterns of HPV integrations and associated host gene expression changes in the context of host chromatin states and topologically associating domains (TADs). HPV integrations were significantly enriched in active chromatin regions and depleted in inactive ones. Interestingly, regardless of chromatin state, genomic regions flanking HPV integrations showed transcriptional upregulation. Nevertheless, upregulation (both local and long-range) was mostly confined to TADs with integration, but not affecting adjacent TADs. Few TADs showed recurrent integrations associated with overexpression of oncogenes within them (e.g. MYC, PVT1, TP63 and ERBB2) regardless of proximity. Hi-C and 4C-seq analyses in cervical cancer cell line (HeLa) demonstrated chromatin looping interactions between integrated HPV and MYC/PVT1 regions (~ 500 kb apart), leading to allele-specific overexpression. Based on these, we propose HPV integrations can trigger multimodal oncogenic activation to promote cancer progression.

Sequence basis of transcription initiation in the human genome

Transcription initiation is a process that is essential to ensuring the proper function of any gene, yet we still lack a unified understanding of sequence patterns and rules that explain most transcription start sites in the human genome. By predicting transcription initiation at base-pair resolution from sequences with a deep learning-inspired explainable model called Puffin, we show that a small set of simple rules can explain transcription initiation at most human promoters. We identify key sequence patterns that contribute to human promoter activity, each activating transcription with distinct position-specific effects. Furthermore, we explain the sequence basis of bidirectional transcription at promoters, identify the links between promoter sequence and gene expression variation across cell types, and explore the conservation of sequence determinants of transcription initiation across mammalian species.

The role of histone acetylation in transcriptional regulation and seed development

Histone acetylation is highly conserved across eukaryotes and has been linked to gene activation since its discovery nearly 60 years ago. Over the past decades, histone acetylation has been evidenced to play crucial roles in plant development and response to various environmental cues. Emerging data indicate that histone acetylation is one of the defining features of "open chromatin," while the role of histone acetylation in transcription remains controversial. In this review, we briefly describe the discovery of histone acetylation, the mechanism of histone acetylation regulating transcription in yeast and mammals, and summarize the research progress of plant histone acetylation. Furthermore, we also emphasize the effect of histone acetylation on seed development and its potential use in plant breeding. A comprehensive knowledge of histone acetylation might provide new and more flexible research perspectives to enhance crop yield and stress resistance.

The melon Fom-1-Prv resistance gene pair: Correlated spatial expression and interaction with a viral protein

The head-to-head oriented pair of melon resistance genes, Fom-1 and Prv, control resistance to Fusarium oxysporum races 0 and 2 and papaya ringspot virus (PRSV), respectively. They encode, via several RNA splice variants, TIR-NBS-LRR proteins, and Prv has a C-terminal extra domain with a second NBS homologous sequence. In other systems, paired R-proteins were shown to operate by "labor division," with one protein having an extra integrated domain that directly binds the pathogen's Avr factor, and the second protein executing the defense response. We report that the expression of the two genes in two pairs of near-isogenic lines was higher in the resistant isoline and inducible by F. oxysporum race 2 but not by PRSV. The intergenic DNA region separating the coding sequences of the two genes acted as a bi-directional promoter and drove GUS expression in transgenic melon roots and transgenic tobacco plants. Expression of both genes was strong in melon root tips, around the root vascular cylinder, and the phloem and xylem parenchyma of tobacco stems and petioles. The pattern of GUS expression suggests coordinated expression of the two genes. In agreement with the above model, Prv's extra domain was shown to interact with the cylindrical inclusion protein of PRSV both in yeast cells and in planta.

Pig H3K4me3, H3K27ac, and gene expression profiles reveal reproductive tissue-specific activity of transposable elements

Regulatory sequences and transposable elements (TEs) account for a large proportion of the genomic sequences of species; however, their roles in gene transcription, especially tissue-specific expression, remain largely unknown. Pigs serve as an excellent animal model for studying genomic sequence biology due to the extensive diversity among their wild and domesticated populations. Here, we conducted an integrated analysis using H3K27ac ChIP-seq, H3K4me3 ChIP-seq, and RNA-seq data from 10 different tissues of seven fetuses and eight closely related adult pigs. We aimed to annotate the regulatory elements and TEs to elucidate their associations with histone modifications and mRNA expression across different tissues and developmental stages. Based on correlation analysis between mRNA expression and H3K27ac and H3K4me3 peak activity, results indicated that H3K27ac exhibited stronger associations with gene expression than H3K4me3. Furthermore, 1.45% of TEs overlapped with either the H3K27ac or H3K4me3 peaks, with the majority displaying tissue-specific activity. Notably, a TE subfamily (LTR4C_SS), containing binding motifs for SIX1 and SIX4, showed specific enrichment in the H3K27ac peaks of the adult and fetal ovaries. RNA-seq analysis also revealed widespread expression of TEs in the exons or promoters of genes, including 4 688 TE-containing transcripts with distinct development stage-specific and tissue-specific expression. Of note, 1 967 TE-containing transcripts were enriched in the testes. We identified a long terminal repeat (LTR), MLT1F1, acting as a testis-specific alternative promoter in SRPK2 (a cell cycle-related protein kinase) in our pig dataset. This element was also conserved in humans and mice, suggesting either an ancient integration of TEs in genes specifically expressed in the testes or parallel evolutionary patterns. Collectively, our findings demonstrate that TEs are deeply embedded in the genome and exhibit important tissue-specific biological functions, particularly in the reproductive organs.

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.

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.

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.

Engineering and Implementation of Synthetic Molecular Tools in the Basidiomycete Fungus Ustilago maydis

The basidiomycete Ustilago maydis is a well-characterized model organism for studying pathogen-host interactions and of great interest for a broad spectrum of biotechnological applications. To facilitate research and enable applications, in this study, three luminescence-based and one enzymatic quantitative reporter were implemented and characterized. Several dual-reporter constructs were generated for ratiometric normalization that can be used as a fast-screening platform for reporter gene expression, applicable to in vitro and in vivo detection. Furthermore, synthetic bidirectional promoters that enable bicisitronic expression for gene expression studies and engineering strategies were constructed and implemented. These noninvasive, quantitative reporters and expression tools will significantly widen the application range of biotechnology in U. maydis and enable the in planta detection of fungal infection.

CTCF blocks antisense transcription initiation at divergent promoters

Transcription at most promoters is divergent, initiating at closely spaced oppositely oriented core promoters to produce sense transcripts along with often unstable upstream antisense transcripts (uasTrx). How antisense transcription is regulated and to what extent it is coordinated with sense transcription is not well understood. Here, by combining acute degradation of the multi-functional transcription factor CTCF and nascent transcription measurements, we find that CTCF specifically suppresses antisense but not sense transcription at hundreds of divergent promoters. Primary transcript RNA-FISH shows that CTCF lowers burst fraction but not burst intensity of uasTrx and that co-bursting of sense and antisense transcripts is disfavored. Genome editing, chromatin conformation studies and high-resolution transcript mapping revealed that precisely positioned CTCF directly suppresses the initiation of uasTrx, in a manner independent of its architectural function. In sum, CTCF shapes the transcriptional landscape in part by suppressing upstream antisense transcription.

CTCF blocks antisense transcription initiation at divergent promoters

Transcription at most promoters is divergent, initiating at closely spaced oppositely oriented core promoters to produce sense transcripts along with often unstable upstream antisense transcripts (uasTrx). How antisense transcription is regulated and to what extent it is coordinated with sense transcription is not well understood. Here, by combining acute degradation of the multi-functional transcription factor CTCF and nascent transcription measurements, we find that CTCF specifically suppresses antisense but not sense transcription at hundreds of divergent promoters. Primary transcript RNA-FISH shows that CTCF lowers burst fraction but not burst intensity of uasTrx and that co-bursting of sense and antisense transcripts is disfavored. Genome editing, chromatin conformation studies and high-resolution transcript mapping revealed that precisely positioned CTCF directly suppresses the initiation of uasTrx, in a manner independent of its architectural function. In sum, CTCF shapes the transcriptional landscape in part by suppressing upstream antisense transcription.

Nuclear microRNAs release paused Pol II via the DDX21-CDK9 complex

RNA activation (RNAa) is an uncharacterized mechanism of transcriptional activation mediated by small RNAs, such as microRNAs (miRNAs). A critical issue in RNAa research is that it is difficult to distinguish between changes in gene expression caused indirectly by post-transcriptional regulation and direct induction of gene expression by RNAa. Therefore, in this study, we seek to identify a key factor involved in RNAa, using the induction of ZMYND10 by miR-34a as a system to evaluate RNAa. We identify the positive transcription elongation factors CDK9 and DDX21, which form a complex with nuclear AGO and TNRC6A, as important transcriptional activators of RNAa. In addition, we find that inhibition of DDX21 suppresses RNAa by miR-34a and other miRNAs without inhibiting post-transcriptional regulation. Our findings reveal a strong connection between RNAa and release of paused Pol II, facilitating RNAa research by making it possible to separately analyze post-transcriptional regulation and RNAa.

Decoding Transcription Regulatory Mechanisms Associated with Coccidioides immitis Phase Transition Using Total RNA

New or emerging infectious diseases are commonly caused by pathogens that cannot be readily manipulated or studied under common laboratory conditions. These limitations hinder standard experimental approaches and our abilities to define the fundamental molecular mechanisms underlying pathogenesis. The advance of capped small RNA sequencing (csRNA-seq) now enables genome-wide mapping of actively initiated transcripts from genes and other regulatory transcribed start regions (TSRs) such as enhancers at a precise moment from total RNA. As RNA is nonpathogenic and can be readily isolated from inactivated infectious samples, csRNA-seq can detect acute changes in gene regulation within or in response to a pathogen with remarkable sensitivity under common laboratory conditions. Studying valley fever (coccidioidomycosis), an emerging endemic fungal infection that increasingly impacts livestock, pet, and human health, we show how csRNA-seq can unravel transcriptional programs driving pathogenesis. Performing csRNA-seq on RNA isolated from different stages of the valley fever pathogen Coccidioides immitis revealed alternative promoter usage, connected cis-regulatory domains, and a WOPR family transcription factor, which are known regulators of virulence in other fungi, as being critical for pathogenic growth. We further demonstrate that a C. immitis WOPR homologue, CIMG_02671, activates transcription in a WOPR motif-dependent manner. Collectively, these findings provide novel insights into valley fever pathogenesis and provide a proof of principle for csRNA-seq as a powerful means to determine the genes, regulatory mechanisms, and transcription factors that control the pathogenesis of highly infectious agents. IMPORTANCE Infectious pathogens like airborne viruses or fungal spores are difficult to study; they require high-containment facilities, special equipment, and expertise. As such, establishing approaches such as genome editing or other means to identify the factors and mechanisms underlying caused diseases, and, thus, promising drug targets, is costly and time-intensive. These obstacles particularly hinder the analysis of new, emerging, or rare infectious diseases. We recently developed a method termed capped small RNA sequencing (csRNA-seq) that enables capturing acute changes in active gene expression from total RNA. Prior to csRNA-seq, such an analysis was possible only by using living cells or nuclei, in which pathogens are highly infectious. The process of RNA purification, however, inactivates pathogens and thus enables the analysis of gene expression during disease progression under standard laboratory conditions. As a proof of principle, here, we use csRNA-seq to unravel the gene regulatory programs and factors likely critical for the pathogenesis of valley fever, an emerging endemic fungal infection that increasingly impacts livestock, pet, and human health.

Single-cell-resolved dynamics of chromatin architecture delineate cell and regulatory states in zebrafish embryos

DNA accessibility of cis-regulatory elements (CREs) dictates transcriptional activity and drives cell differentiation during development. While many genes regulating embryonic development have been identified, the underlying CRE dynamics controlling their expression remain largely uncharacterized. To address this, we produced a multimodal resource and genomic regulatory map for the zebrafish community, which integrates single-cell combinatorial indexing assay for transposase-accessible chromatin with high-throughput sequencing (sci-ATAC-seq) with bulk histone PTMs and Hi-C data to achieve a genome-wide classification of the regulatory architecture determining transcriptional activity in the 24-h post-fertilization (hpf) embryo. We characterized the genome-wide chromatin architecture at bulk and single-cell resolution, applying sci-ATAC-seq on whole 24-hpf stage zebrafish embryos, generating accessibility profiles for ∼23,000 single nuclei. We developed a genome segmentation method, ScregSeg (single-cell regulatory landscape segmentation), for defining regulatory programs, and candidate CREs, specific to one or more cell types. We integrated the ScregSeg output with bulk measurements for histone post-translational modifications and 3D genome organization and identified new regulatory principles between chromatin modalities prevalent during zebrafish development. Sci-ATAC-seq profiling of npas4l/cloche mutant embryos identified novel cellular roles for this hematovascular transcriptional master regulator and suggests an intricate mechanism regulating its expression. Our work defines regulatory architecture and principles in the zebrafish embryo and establishes a resource of cell-type-specific genome-wide regulatory annotations and candidate CREs, providing a valuable open resource for genomics, developmental, molecular, and computational biology.

RN7SK small nuclear RNA controls bidirectional transcription of highly expressed gene pairs in skin

Pausing of RNA polymerase II (Pol II) close to promoters is a common regulatory step in RNA synthesis, and is coordinated by a ribonucleoprotein complex scaffolded by the noncoding RNA RN7SK. The function of RN7SK-regulated gene transcription in adult tissue homoeostasis is currently unknown. Here, we deplete RN7SK during mouse and human epidermal stem cell differentiation. Unexpectedly, loss of this small nuclear RNA specifically reduces transcription of numerous cell cycle regulators leading to cell cycle exit and differentiation. Mechanistically, we show that RN7SK is required for efficient transcription of highly expressed gene pairs with bidirectional promoters, which in the epidermis co-regulated cell cycle and chromosome organization. The reduction in transcription involves impaired splicing and RNA decay, but occurs in the absence of chromatin remodelling at promoters and putative enhancers. Thus, RN7SK is directly required for efficient Pol II transcription of highly transcribed bidirectional gene pairs, and thereby exerts tissue-specific functions, such as maintaining a cycling cell population in the epidermis.

Investigating crosstalk between H3K27 acetylation and H3K4 trimethylation in CRISPR/dCas-based epigenome editing and gene activation

Epigenome editing methods enable the precise manipulation of epigenetic modifications, such as histone posttranscriptional modifications (PTMs), for uncovering their biological functions. While histone PTMs have been correlated with certain gene expression status, the causalities remain elusive. Histone H3 Lysine 27 acetylation (H3K27ac) and histone H3 Lysine 4 trimethylation (H3K4me3) are both associated with active genes, and located at active promoters and enhancers or around transcriptional start sites (TSSs). Although crosstalk between histone lysine acetylation and H3K4me3 has been reported, relationships between specific epigenetic marks during transcriptional activation remain largely unclear. Here, using clustered regularly interspaced short palindromic repeats (CRISPR)/dCas-based epigenome editing methods, we discovered that the ectopic introduction of H3K27ac in the promoter region lead to H3K4me3 enrichment around TSS and transcriptional activation, while H3K4me3 installation at the promoter cannot induce H3K27ac increase and failed to activate gene expression. Blocking the reading of H3K27ac by BRD proteins using inhibitor JQ1 abolished H3K27ac-induced H3K4me3 installation and downstream gene activation. Furthermore, we uncovered that BRD2, not BRD4, mediated H3K4me3 installation and gene activation upon H3K27ac writing. Our studies revealed the relationships between H3K27ac and H3K4me3 in gene activation process and demonstrated the application of CRISPR/dCas-based epigenome editing methods in elucidating the crosstalk between epigenetic mechanisms.

The evolutionary acquisition and mode of functions of promoter-associated non-coding RNAs (pancRNAs) for mammalian development

Increasing evidence has shown that many long non-coding RNAs (lncRNAs) are involved in gene regulation in a variety of ways such as transcriptional, post-transcriptional and epigenetic regulation. Promoter-associated non-coding RNAs (pancRNAs), which are categorized into the most abundant single-copy lncRNA biotype, play vital regulatory roles in finely tuning cellular specification at the epigenomic level. In short, pancRNAs can directly or indirectly regulate downstream genes to participate in the development of organisms in a cell-specific manner. In this review, we will introduce the evolutionarily acquired characteristics of pancRNAs as determined by comparative epigenomics and elaborate on the research progress on pancRNA-involving processes in mammalian embryonic development, including neural differentiation.

A Chinese hamster transcription start site atlas that enables targeted editing of CHO cells

Chinese hamster ovary (CHO) cells are widely used for producing biopharmaceuticals, and engineering gene expression in CHO is key to improving drug quality and affordability. However, engineering gene expression or activating silent genes requires accurate annotation of the underlying regulatory elements and transcription start sites (TSSs). Unfortunately, most TSSs in the published Chinese hamster genome sequence were computationally predicted and are frequently inaccurate. Here, we use nascent transcription start site sequencing methods to revise TSS annotations for 15 308 Chinese hamster genes and 3034 non-coding RNAs based on experimental data from CHO-K1 cells and 10 hamster tissues. We further capture tens of thousands of putative transcribed enhancer regions with this method. Our revised TSSs improves upon the RefSeq annotation by revealing core sequence features of gene regulation such as the TATA box and the Initiator and, as exemplified by targeting the glycosyltransferase gene Mgat3, facilitate activating silent genes by CRISPRa. Together, we envision our revised annotation and data will provide a rich resource for the CHO community, improve genome engineering efforts and aid comparative and evolutionary studies.

PEPPRO: quality control and processing of nascent RNA profiling data

Nascent RNA profiling is growing in popularity; however, there is no standard analysis pipeline to uniformly process the data and assess quality. Here, we introduce PEPPRO, a comprehensive, scalable workflow for GRO-seq, PRO-seq, and ChRO-seq data. PEPPRO produces uniformly processed output files for downstream analysis and assesses adapter abundance, RNA integrity, library complexity, nascent RNA purity, and run-on efficiency. PEPPRO is restartable and fault-tolerant, records copious logs, and provides a web-based project report. PEPPRO can be run locally or using a cluster, providing a portable first step for genomic nascent RNA analysis.

Genome-wide strand asymmetry in massively parallel reporter activity favors genic strands

Massively parallel reporter assays (MPRAs) are useful tools to characterize regulatory elements in human genomes. An aspect of MPRAs that is not typically the focus of analysis is their intrinsic ability to differentiate activity levels for a given sequence element when placed in both of its possible orientations relative to the reporter construct. Here, we describe pervasive strand asymmetry of MPRA signals in data sets from multiple reporter configurations in both published and newly reported data. These effects are reproducible across different cell types and in different treatments within a cell type and are observed both within and outside of annotated regulatory elements. From elements in gene bodies, MPRA strand asymmetry favors the sense strand, suggesting that function related to endogenous transcription is driving the phenomenon. Similarly, we find that within Alu mobile element insertions, strand asymmetry favors the transcribed strand of the ancestral retrotransposon. The effect is consistent across the multiplicity of Alu elements in human genomes and is more pronounced in less diverged Alu elements. We find sequence features driving MPRA strand asymmetry and show its prediction from sequence alone. We see some evidence for RNA stabilization and transcriptional activation mechanisms and hypothesize that the effect is driven by natural selection favoring efficient transcription. Our results indicate that strand asymmetry is a pervasive and reproducible feature in MPRA data. More importantly, the fact that MPRA asymmetry favors naturally transcribed strands suggests that it stems from preserved biological functions that have a substantial, global impact on gene and genome evolution.

Structural insights into preinitiation complex assembly on core promoters

Transcription factor IID (TFIID) recognizes core promoters and supports preinitiation complex (PIC) assembly for RNA polymerase II (Pol II)-mediated eukaryotic transcription. We determined the structures of human TFIID-based PIC in three stepwise assembly states and revealed two-track PIC assembly: stepwise promoter deposition to Pol II and extensive modular reorganization on track I (on TATA-TFIID-binding element promoters) versus direct promoter deposition on track II (on TATA-only and TATA-less promoters). The two tracks converge at an ~50-subunit holo PIC in identical conformation, whereby TFIID stabilizes PIC organization and supports loading of cyclin-dependent kinase (CDK)-activating kinase (CAK) onto Pol II and CAK-mediated phosphorylation of the Pol II carboxyl-terminal domain. Unexpectedly, TBP of TFIID similarly bends TATA box and TATA-less promoters in PIC. Our study provides structural visualization of stepwise PIC assembly on highly diversified promoters.

Inferring time series chromatin states for promoter-enhancer pairs based on Hi-C data

BACKGROUND

Co-localized combinations of histone modifications ("chromatin states") have been shown to correlate with promoter and enhancer activity. Changes in chromatin states over multiple time points ("chromatin state trajectories") have previously been analyzed at promoter and enhancers separately. With the advent of time series Hi-C data it is now possible to connect promoters and enhancers and to analyze chromatin state trajectories at promoter-enhancer pairs.

RESULTS

We present TimelessFlex, a framework for investigating chromatin state trajectories at promoters and enhancers and at promoter-enhancer pairs based on Hi-C information. TimelessFlex extends our previous approach Timeless, a Bayesian network for clustering multiple histone modification data sets at promoter and enhancer feature regions. We utilize time series ATAC-seq data measuring open chromatin to define promoters and enhancer candidates. We developed an expectation-maximization algorithm to assign promoters and enhancers to each other based on Hi-C interactions and jointly cluster their feature regions into paired chromatin state trajectories. We find jointly clustered promoter-enhancer pairs showing the same activation patterns on both sides but with a stronger trend at the enhancer side. While the promoter side remains accessible across the time series, the enhancer side becomes dynamically more open towards the gene activation time point. Promoter cluster patterns show strong correlations with gene expression signals, whereas Hi-C signals get only slightly stronger towards activation. The code of the framework is available at https://github.com/henriettemiko/TimelessFlex .

CONCLUSIONS

TimelessFlex clusters time series histone modifications at promoter-enhancer pairs based on Hi-C and it can identify distinct chromatin states at promoter and enhancer feature regions and their changes over time.

Transcription shapes genome-wide histone acetylation patterns

Histone acetylation is a ubiquitous hallmark of transcription, but whether the link between histone acetylation and transcription is causal or consequential has not been addressed. Using immunoblot and chromatin immunoprecipitation-sequencing in S. cerevisiae, here we show that the majority of histone acetylation is dependent on transcription. This dependency is partially explained by the requirement of RNA polymerase II (RNAPII) for the interaction of H4 histone acetyltransferases (HATs) with gene bodies. Our data also confirms the targeting of HATs by transcription activators, but interestingly, promoter-bound HATs are unable to acetylate histones in the absence of transcription. Indeed, HAT occupancy alone poorly predicts histone acetylation genome-wide, suggesting that HAT activity is regulated post-recruitment. Consistent with this, we show that histone acetylation increases at nucleosomes predicted to stall RNAPII, supporting the hypothesis that this modification is dependent on nucleosome disruption during transcription. Collectively, these data show that histone acetylation is a consequence of RNAPII promoting both the recruitment and activity of histone acetyltransferases.

Transcriptionally active enhancers in human cancer cells

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

Crosstalk of promoter and terminator during RNA polymerase II transcription cycle

The transcription cycle of RNAPII is comprised of three consecutive steps; initiation, elongation and termination. It has been assumed that the initiation and termination steps occur in spatial isolation, essentially as independent events. A growing body of evidence, however, has challenged this dogma. First, factors involved in initiation and termination exhibit both a genetic and a physical interaction during transcription. Second, the initiation and termination factors have been found to occupy both ends of a transcribing gene. Third, physical interaction of initiation and termination factors occupying distal ends of a gene sometime results in the entire terminator region of a genes looping back and contact its cognate promoter, thereby forming a looped gene architecture during transcription. A logical interpretation of these findings is that the initiation and termination steps of transcription do not occur in isolation. There is extensive communication of factors occupying promoter and terminator ends of a gene during transcription cycle. This review entails a discussion of the promoter-terminator crosstalk and its implication in the context of transcription.

Identification of the human DPR core promoter element using machine learning

The RNA polymerase II (Pol II) core promoter is the strategic site of convergence of the signals that lead to transcription initiation^1-5, but the downstream core promoter in humans has been difficult to decipher^1-3. Here, we analyze the human Pol II core promoter and use machine learning to generate predictive models for the downstream core promoter region (DPR) and the TATA box. We developed a method termed HARPE (high-throughput analysis of randomized promoter elements) to create hundreds of thousands of DPR (or TATA box) variants that are each of known transcriptional strength. We then analyzed the HARPE data by support vector regression (SVR) to provide comprehensive models for the sequence motifs, and found that the SVR-based approach is more effective than a consensus-based method for predicting transcriptional activity. These studies revealed that the DPR is a functionally important core promoter element that is widely used in human promoters. Importantly, there appears to be a duality between the DPR and TATA box, as many promoters contain one or the other element. More broadly, these findings show that functional DNA motifs can be identified by machine learning analysis of a comprehensive set of sequence variants.

Widespread Transcriptional Scanning in the Testis Modulates Gene Evolution Rates

The testis expresses the largest number of genes of any mammalian organ, a finding that has long puzzled molecular biologists. Our single-cell transcriptomic data of human and mouse spermatogenesis provide evidence that this widespread transcription maintains DNA sequence integrity in the male germline by correcting DNA damage through a mechanism we term transcriptional scanning. We find that genes expressed during spermatogenesis display lower mutation rates on the transcribed strand and have low diversity in the population. Moreover, this effect is fine-tuned by the level of gene expression during spermatogenesis. The unexpressed genes, which in our model do not benefit from transcriptional scanning, diverge faster over evolutionary timescales and are enriched for sensory and immune-defense functions. Collectively, we propose that transcriptional scanning shapes germline mutation signatures and modulates mutation rates in a gene-specific manner, maintaining DNA sequence integrity for the bulk of genes but allowing for faster evolution in a specific subset.

Universal promoter scanning by Pol II during transcription initiation in Saccharomyces cerevisiae

BACKGROUND

The majority of eukaryotic promoters utilize multiple transcription start sites (TSSs). How multiple TSSs are specified at individual promoters across eukaryotes is not understood for most species. In Saccharomyces cerevisiae, a pre-initiation complex (PIC) comprised of Pol II and conserved general transcription factors (GTFs) assembles and opens DNA upstream of TSSs. Evidence from model promoters indicates that the PIC scans from upstream to downstream to identify TSSs. Prior results suggest that TSS distributions at promoters where scanning occurs shift in a polar fashion upon alteration in Pol II catalytic activity or GTF function.

RESULTS

To determine the extent of promoter scanning across promoter classes in S. cerevisiae, we perturb Pol II catalytic activity and GTF function and analyze their effects on TSS usage genome-wide. We find that alterations to Pol II, TFIIB, or TFIIF function widely alter the initiation landscape consistent with promoter scanning operating at all yeast promoters, regardless of promoter class. Promoter architecture, however, can determine the extent of promoter sensitivity to altered Pol II activity in ways that are predicted by a scanning model.

CONCLUSIONS

Our observations coupled with previous data validate key predictions of the scanning model for Pol II initiation in yeast, which we term the shooting gallery. In this model, Pol II catalytic activity and the rate and processivity of Pol II scanning together with promoter sequence determine the distribution of TSSs and their usage.

Structure and mechanism of the RNA polymerase II transcription machinery

RNA polymerase II (Pol II) transcribes all protein-coding genes and many noncoding RNAs in eukaryotic genomes. Although Pol II is a complex, 12-subunit enzyme, it lacks the ability to initiate transcription and cannot consistently transcribe through long DNA sequences. To execute these essential functions, an array of proteins and protein complexes interact with Pol II to regulate its activity. In this review, we detail the structure and mechanism of over a dozen factors that govern Pol II initiation (e.g., TFIID, TFIIH, and Mediator), pausing, and elongation (e.g., DSIF, NELF, PAF, and P-TEFb). The structural basis for Pol II transcription regulation has advanced rapidly in the past decade, largely due to technological innovations in cryoelectron microscopy. Here, we summarize a wealth of structural and functional data that have enabled a deeper understanding of Pol II transcription mechanisms; we also highlight mechanistic questions that remain unanswered or controversial.

Lessons from eRNAs: understanding transcriptional regulation through the lens of nascent RNAs

Nascent transcription assays, such as global run-on sequencing (GRO-seq) and precision run-on sequencing (PRO-seq), have uncovered a myriad of unstable RNAs being actively produced from numerous sites genome-wide. These transcripts provide a more complete and immediate picture of the impact of regulatory events. Transcription factors recruit RNA polymerase II, effectively initiating the process of transcription; repressors inhibit polymerase recruitment. Efficiency of recruitment is dictated by sequence elements in and around the RNA polymerase loading zone. A combination of sequence elements and RNA binding proteins subsequently influence the ultimate stability of the resulting transcript. Some of these transcripts are capable of providing feedback on the process, influencing subsequent transcription. By monitoring RNA polymerase activity, nascent assays provide insights into every step of the regulated process of transcription.

Nuclear sorting of RNA

The majority of the mammalian genome is transcribed by RNA polymerase II, yielding a vast amount of noncoding RNA (ncRNA) in addition to the standard production of mRNA. The typical nuclear biogenesis of mRNA relies on the tightly controlled coupling of co- and post-transcriptional processing events, which ultimately results in the export of transcripts into the cytoplasm. These processes are subject to surveillance by nuclear RNA decay pathways to prevent the export of aberrant, or otherwise "non-optimal," transcripts. However, unlike mRNA, many long ncRNAs are nuclear retained and those that maintain enduring functions must employ precautions to evade decay. Proper sorting and localization of RNA is therefore an essential activity in eukaryotic cells and the formation of ribonucleoprotein complexes during early stages of RNA synthesis is central to deciding such transcript fate. This review details our current understanding of the pathways and factors that direct RNAs towards a particular destiny and how transcripts combat the adverse conditions of the nucleus. This article is categorized under: RNA Export and Localization > Nuclear Export/Import RNA Turnover and Surveillance > Turnover/Surveillance Mechanisms RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications.

R-Loops Promote Antisense Transcription across the Mammalian Genome

Widespread antisense long noncoding RNA (lncRNA) overlap with many protein-coding genes in mammals and emanate from gene promoter, enhancer, and termination regions. However, their origin and biological purpose remain unclear. We show that these antisense lncRNA can be generated by R-loops that form when nascent transcript invades the DNA duplex behind elongating RNA polymerase II (Pol II). Biochemically, R-loops act as intrinsic Pol II promoters to induce de novo RNA synthesis. Furthermore, their removal across the human genome by RNase H1 overexpression causes the selective reduction of antisense transcription. Consequently, we predict that R-loops act to facilitate the synthesis of many gene proximal antisense lncRNA. Not only are R-loops widely associated with DNA damage and repair, but we now show that they have the capacity to promote de novo transcript synthesis that may have aided the evolution of gene regulation.

•

R-loops formed within plasmids promote antisense transcription in nuclear extracts

•

TSS of lncRNA and eRNA are often near R-loop structures and sensitive to RNase H1

•

Preinitiation complexes associated with lncRNA synthesis are R-loop dependent

•

Many mammalian lncRNA derive from R-loop promoter activity

Identification and dynamic quantification of regulatory elements using total RNA

The spatial and temporal regulation of transcription initiation is pivotal for controlling gene expression. Here, we introduce capped-small RNA-seq (csRNA-seq), which uses total RNA as starting material to detect transcription start sites (TSSs) of both stable and unstable RNAs at single-nucleotide resolution. csRNA-seq is highly sensitive to acute changes in transcription and identifies an order of magnitude more regulated transcripts than does RNA-seq. Interrogating tissues from species across the eukaryotic kingdoms identified unstable transcripts resembling enhancer RNAs, pri-miRNAs, antisense transcripts, and promoter upstream transcripts in multicellular animals, plants, and fungi spanning 1.6 billion years of evolution. Integration of epigenomic data from these organisms revealed that histone H3 trimethylation (H3K4me3) was largely confined to TSSs of stable transcripts, whereas H3K27ac marked nucleosomes downstream from all active TSSs, suggesting an ancient role for posttranslational histone modifications in transcription. Our findings show that total RNA is sufficient to identify transcribed regulatory elements and capture the dynamics of initiated stable and unstable transcripts at single-nucleotide resolution in eukaryotes.

Transcription start site analysis reveals widespread divergent transcription in D. melanogaster and core promoter-encoded enhancer activities

Mammalian gene promoters and enhancers share many properties. They are composed of a unified promoter architecture of divergent transcripton initiation and gene promoters may exhibit enhancer function. However, it is currently unclear how expression strength of a regulatory element relates to its enhancer strength and if the unifying architecture is conserved across Metazoa. Here we investigate the transcription initiation landscape and its associated RNA decay in Drosophila melanogaster. We find that the majority of active gene-distal enhancers and a considerable fraction of gene promoters are divergently transcribed. We observe quantitative relationships between enhancer potential, expression level and core promoter strength, providing an explanation for indirectly related histone modifications that are reflecting expression levels. Lowly abundant unstable RNAs initiated from weak core promoters are key characteristics of gene-distal developmental enhancers, while the housekeeping enhancer strengths of gene promoters reflect their expression strengths. The seemingly separable layer of regulation by gene promoters with housekeeping enhancer potential is also indicated by chromatin interaction data. Our results suggest a unified promoter architecture of many D. melanogaster regulatory elements, that is universal across Metazoa, whose regulatory functions seem to be related to their core promoter elements.

A step-by-step guide to analyzing CAGE data using R/Bioconductor

Cap Analysis of Gene Expression (CAGE) is one of the most popular 5'-end sequencing methods. In a single experiment, CAGE can be used to locate and quantify the expression of both Transcription Start Sites (TSSs) and enhancers. This is workflow is a case study on how to use the CAGEfightR package to orchestrate analysis of CAGE data within the Bioconductor project. This workflow starts from BigWig-files and covers both basic CAGE analyses such as identifying, quantifying and annotating TSSs and enhancers, advanced analysis such as finding interacting TSS-enhancer pairs and enhancer clusters, to differential expression analysis and alternative TSS usage. R-code, discussion and references are intertwined to help provide guidelines for future CAGE studies of the same kind.

Identification of De Novo Enhancers Activated by TGFβ to Drive Expression of CDKN2A and B in HeLa Cells

Disruption of the CDKN2A (INK4A/ARF) and B (INK4B) genes, which encode three function-independent tumor suppressors, is one of the most common events in human cancer. Because their relative importance in tumor prevention appears to be species- and context-specific, studying their regulation can shed light on mechanisms by which they are bypassed in malignant transformation. We previously unveiled a new pathway in which TGFβ selectively induces Arf at mouse Cdkn2a in eye development and cultured fibroblasts. As TGFβ signaling is often derailed in cancer development or progression, we investigated its control of CDKN2A/B in human cancer. Computational analyses of sequencing and array data from nearly 11,000 patients with cancer in TCGA showed discordant expression of ARF and INK4A in most cancer subtypes, with gene copy-number loss and promoter methylation involved in only a subset. Using HeLa cells as a model, we found that exogenous TGFβ induced ARF mRNA and protein, and ARF knockdown limited TGFβ-mediated growth suppression. TGFβ-mediated ARF mRNA induction required SMAD2/3, p38MAPK, and SP1, and ARF mRNA was induced without added RNAPII recruitment. Chromatin immunoprecipitation unveiled a remote enhancer element engaged by TGFβ by a mechanism that partially depended on p38MAPK. CRISPR-based editing of this enhancer limited induction of ARF and INK4B by TGFβ, but not by oncogenic RAS. IMPLICATIONS: Our findings reveal new molecular mechanisms by which CDKN2A/B regulation is coupled to external cues, and those findings represent entry points to further explore pharmacologic strategies to restore their expression in cancer.

Differential regulation of CpG island methylation within divergent and unidirectional promoters in colorectal cancer

The silencing of tumor suppressor genes by promoter CpG island (CGI) methylation is an important cause of oncogenesis. Silencing of MLH1 and BRCA1, two examples of oncogenic events, results from promoter CGI methylation. Interestingly, both MLH1 and BRCA1 have a divergent promoter, from which another gene on the opposite strand is also transcribed. Although studies have shown that divergent transcription is an important factor in transcriptional regulation, little is known about its implication in aberrant promoter methylation in cancer. In this study, we analyzed the methylation status of CGI in divergent promoters using a recently enriched transcriptome database. We measured the extent of CGI methylation in 119 colorectal cancer (CRC) clinical samples (65 microsatellite instability high [MSI‐H] CRC with CGI methylator phenotype, 28 MSI‐H CRC without CGI methylator phenotype and 26 microsatellite stable CRC) and 21 normal colorectal tissues using Infinium MethylationEPIC BeadChip. We found that CGI within divergent promoters are less frequently methylated than CGI within unidirectional promoters in normal cells. In the genome of CRC cells, CGI within unidirectional promoters are more vulnerable to aberrant methylation than CGI within divergent promoters. In addition, we identified three DNA sequence motifs that correlate with methylated CGI. We also showed that methylated CGI are associated with genes whose expression is low in normal cells. Thus, we here provide fundamental observations regarding the methylation of divergent promoters that are essential for the understanding of carcinogenesis and development of cancer prevention strategies.

Identification of regulatory elements from nascent transcription using dREG

Our genomes encode a wealth of transcription initiation regions (TIRs) that can be identified by their distinctive patterns of actively elongating RNA polymerase. We previously introduced dREG to identify TIRs using PRO-seq data. Here, we introduce an efficient new implementation of dREG that uses PRO-seq data to identify both uni- and bidirectionally transcribed TIRs with 70% improvement in accuracy, three- to fourfold higher resolution, and >100-fold increases in computational efficiency. Using a novel strategy to identify TIRs based on their statistical confidence reveals extensive overlap with orthogonal assays, yet also reveals thousands of additional weakly transcribed TIRs that were not identified by H3K27ac ChIP-seq or DNase-seq. Novel TIRs discovered by dREG were often associated with RNA polymerase III initiation, bound by pioneer transcription factors, or located in broad domains marked by repressive chromatin modifications. Our results suggest that transcription initiation can be a powerful tool for expanding the catalog of functional elements.

Systematic interrogation of human promoters

Despite much research, our understanding of the architecture and cis-regulatory elements of human promoters is still lacking. Here, we devised a high-throughput assay to quantify the activity of approximately 15,000 fully designed sequences that we integrated and expressed from a fixed location within the human genome. We used this method to investigate thousands of native promoters and preinitiation complex (PIC) binding regions followed by in-depth characterization of the sequence motifs underlying promoter activity, including core promoter elements and TF binding sites. We find that core promoters drive transcription mostly unidirectionally and that sequences originating from promoters exhibit stronger activity than those originating from enhancers. By testing multiple synthetic configurations of core promoter elements, we dissect the motifs that positively and negatively regulate transcription as well as the effect of their combinations and distances, including a 10-bp periodicity in the optimal distance between the TATA and the initiator. By comprehensively screening 133 TF binding sites, we find that in contrast to core promoters, TF binding sites maintain similar activity levels in both orientations, supporting a model by which divergent transcription is driven by two distinct unidirectional core promoters sharing bidirectional TF binding sites. Finally, we find a striking agreement between the effect of binding site multiplicity of individual TFs in our assay and their tendency to appear in homotypic clusters throughout the genome. Overall, our study systematically assays the elements that drive expression in core and proximal promoter regions and sheds light on organization principles of regulatory regions in the human genome.