Different distribution of histone modifications in genes with unidirectional and bidirectional transcription and a role of CTCF and cohesin in directing transcription

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.

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.

Pan-Cancer Analysis Reveals the Prognostic Potential of the THAP9/THAP9-AS1 Sense-Antisense Gene Pair in Human Cancers

Human THAP9, which encodes a domesticated transposase of unknown function, and lncRNA THAP9-AS1 (THAP9-antisense1) are arranged head-to-head on opposite DNA strands, forming a sense and antisense gene pair. We predict that there is a bidirectional promoter that potentially regulates the expression of THAP9 and THAP9-AS1. Although both THAP9 and THAP9-AS1 are reported to be involved in various cancers, their correlative roles on each other's expression has not been explored. We analyzed the expression levels, prognosis, and predicted biological functions of the two genes across different cancer datasets (TCGA, GTEx). We observed that although the expression levels of the two genes, THAP9 and THAP9-AS1, varied in different tumors, the expression of the gene pair was strongly correlated with patient prognosis; higher expression of the gene pair was usually linked to poor overall and disease-free survival. Thus, THAP9 and THAP9-AS1 may serve as potential clinical biomarkers of tumor prognosis. Further, we performed a gene co-expression analysis (using WGCNA) followed by a differential gene correlation analysis (DGCA) across 22 cancers to identify genes that share the expression pattern of THAP9 and THAP9-AS1. Interestingly, in both normal and cancer samples, THAP9 and THAP9-AS1 often co-express; moreover, their expression is positively correlated in each cancer type, suggesting the coordinated regulation of this H2H gene pair.

Co-regulation and function of FOXM1/RHNO1 bidirectional genes in cancer

The FOXM1 transcription factor is an oncoprotein and a top biomarker of poor prognosis in human cancer. Overexpression and activation of FOXM1 is frequent in high-grade serous carcinoma (HGSC), the most common and lethal form of human ovarian cancer, and is linked to copy number gains at chromosome 12p13.33. We show that FOXM1 is co-amplified and co-expressed with RHNO1, a gene involved in the ATR-Chk1 signaling pathway that functions in the DNA replication stress response. We demonstrate that FOXM1 and RHNO1 are head-to-head (i.e., bidirectional) genes (BDG) regulated by a bidirectional promoter (BDP) (named F/R-BDP). FOXM1 and RHNO1 each promote oncogenic phenotypes in HGSC cells, including clonogenic growth, DNA homologous recombination repair, and poly-ADP ribosylase inhibitor resistance. FOXM1 and RHNO1 are one of the first examples of oncogenic BDG, and therapeutic targeting of FOXM1/RHNO1 BDG is a potential therapeutic approach for ovarian and other cancers.

A Role for the Mre11-Rad50-Xrs2 Complex in Gene Expression and Chromosome Organization

Mre11-Rad50-Xrs2 (MRX) is a highly conserved complex with key roles in various aspects of DNA repair. Here, we report a new function for MRX in limiting transcription in budding yeast. We show that MRX interacts physically and colocalizes on chromatin with the transcriptional co-regulator Mediator. MRX restricts transcription of coding and noncoding DNA by a mechanism that does not require the nuclease activity of Mre11. MRX is required to tether transcriptionally active loci to the nuclear pore complex (NPC), and it also promotes large-scale gene-NPC interactions. Moreover, MRX-mediated chromatin anchoring to the NPC contributes to chromosome folding and helps to control gene expression. Together, these findings indicate that MRX has a role in transcription and chromosome organization that is distinct from its known function in DNA repair.

LncRNA-MAP3K4 regulates vascular inflammation through the p38 MAPK signaling pathway and cis-modulation of MAP3K4

Chronic vascular inflammation plays a key role in the pathogenesis of atherosclerosis. Long non-coding RNAs (lncRNAs) have emerged as essential inflammation regulators. We identify a novel lncRNA termed lncRNA-MAP3K4 that is enriched in the vessel wall and regulates vascular inflammation. In the aortic intima, lncRNA-MAP3K4 expression was reduced by 50% during the progression of atherosclerosis (chronic inflammation) and 70% during endotoxemia (acute inflammation). lncRNA-MAP3K4 knockdown reduced the expression of key inflammatory factors (eg, ICAM-1, E-selectin, MCP-1) in endothelial cells or vascular smooth muscle cells and decreased monocytes adhesion to endothelium, as well as reducing TNF-α, IL-1β, COX2 expression in macrophages. Mechanistically, lncRNA-MAP3K4 regulates inflammation through the p38 MAPK signaling pathway. lncRNA-MAP3K4 shares a bidirectional promoter with MAP3K4, an upstream regulator of the MAPK signaling pathway, and regulates its transcription in cis. lncRNA-MAP3K4 and MAP3K4 show coordinated expression in response to inflammation in vivo and in vitro. Similar to lncRNA-MAP3K4, MAP3K4 knockdown reduced the expression of inflammatory factors in several different vascular cells. Furthermore, lncRNA-MAP3K4 and MAP3K4 knockdown showed cooperativity in reducing inflammation in endothelial cells. Collectively, these findings unveil the role of a novel lncRNA in vascular inflammation by cis-regulating MAP3K4 via a p38 MAPK pathway.

Genome-wide identification of histone methylation (H3K9me2) and acetylation (H4K12ac) marks in two ecotypes of switchgrass (Panicum virgatum L.)

BACKGROUND

Histone modifications play a significant role in the regulation of transcription and various biological processes, such as development and regeneration. Though a few genomic (including DNA methylation patterns) and transcriptomic studies are currently available in switchgrass, the genome-wide distribution of histone modifications has not yet been studied to help elucidate gene regulation and its application to switchgrass improvement.

RESULTS

This study provides a comprehensive epigenomic analyses of two contrasting switchgrass ecotypes, lowland (AP13) and upland (VS16), by employing chromatin immunoprecipitation sequencing (ChIP-Seq) with two histone marks (suppressive- H3K9me2 and active- H4K12ac). In this study, most of the histone binding was in non-genic regions, and the highest enrichment was seen between 0 and 2 kb regions from the transcriptional start site (TSS). Considering the economic importance and potential of switchgrass as a bioenergy crop, we focused on genes, transcription factors (TFs), and pathways that were associated with C4-photosynthesis, biomass, biofuel production, biotic stresses, and abiotic stresses. Using quantitative real-time PCR (qPCR) the relative expression of five genes selected from the phenylpropanoid-monolignol pathway showed preferential binding of acetylation marks in AP13 rather than in VS16.

CONCLUSIONS

The genome-wide histone modifications reported here can be utilized in understanding the regulation of genes important in the phenylpropanoid-monolignol biosynthesis pathway, which in turn, may help understand the recalcitrance associated with conversion of biomass to biofuel, a major roadblock in utilizing lignocellulosic feedstocks.

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.

Integrative analysis of single-cell expression data reveals distinct regulatory states in bidirectional promoters

Background

Bidirectional promoters (BPs) are prevalent in eukaryotic genomes. However, it is poorly understood how the cell integrates different epigenomic information, such as transcription factor (TF) binding and chromatin marks, to drive gene expression at BPs. Single-cell sequencing technologies are revolutionizing the field of genome biology. Therefore, this study focuses on the integration of single-cell RNA-seq data with bulk ChIP-seq and other epigenetics data, for which single-cell technologies are not yet established, in the context of BPs.

Results

We performed integrative analyses of novel human single-cell RNA-seq (scRNA-seq) data with bulk ChIP-seq and other epigenetics data. scRNA-seq data revealed distinct transcription states of BPs that were previously not recognized. We find associations between these transcription states to distinct patterns in structural gene features, DNA accessibility, histone modification, DNA methylation and TF binding profiles.

Conclusions

Our results suggest that a complex interplay of all of these elements is required to achieve BP-specific transcriptional output in this specialized promoter configuration. Further, our study implies that novel statistical methods can be developed to deconvolute masked subpopulations of cells measured with different bulk epigenomic assays using scRNA-seq data.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0236-7) contains supplementary material, which is available to authorized users.

Bidirectional promoters exhibit characteristic chromatin modification signature associated with transcription elongation in both sense and antisense directions

Background

In contrast to unidirectional promoters wherein antisense transcription results in short transcripts which are rapidly degraded, bidirectional promoters produce mature transcripts in both sense and antisense orientation. To understand the molecular mechanism of how productive bidirectional transcription is regulated, we focused on delineating the chromatin signature of bidirectional promoters.

Results

We report generation and utility of a reporter system that enables simultaneous scoring of transcriptional activity in opposite directions. Testing of putative bidirectional promoters in this system demonstrates no measurable bias towards any one direction of transcription. We analyzed the NUP26L-PIH1D3 bidirectional gene pair during Retinoic acid mediated differentiation of embryonic carcinoma cells. In their native context, we observed that the chromatin landscape at and around the transcription regulatory region between the pair of bidirectional genes is modulated in concordance with transcriptional activity of each gene in the pair. We then extended this analysis to 974 bidirectional gene pairs in two different cell lines, H1 human embryonic stem cells and CD4 positive T cells using publicly available ChIP-Seq and RNA-Seq data. Bidirectional gene pairs were classified based on the intergenic distance separating the two TSS of the transcripts analyzed as well as the relative expression of each transcript in a bidirectional gene pair. We report that for the entire range of intergenic distance separating bidirectional genes, the expression profile of such genes (symmetric or asymmetric) matches the histone modification profile of marks associated with active transcription initiation and elongation.

Conclusions

We demonstrate unique distribution of histone modification marks that correlate robustly with the transcription status of genes regulated by bidirectional promoters. These findings strongly imply that occurrence of these marks might signal the transcription machinery to drive maturation of antisense transcription from the bidirectional promoters.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4697-7) contains supplementary material, which is available to authorized users.

Widespread activation of antisense transcription of the host genome during herpes simplex virus 1 infection

Background

Herpesviruses can infect a wide range of animal species. Herpes simplex virus 1 (HSV-1) is one of the eight herpesviruses that can infect humans and is prevalent worldwide. Herpesviruses have evolved multiple ways to adapt the infected cells to their needs, but knowledge about these transcriptional and post-transcriptional modifications is sparse.

Results

Here, we show that HSV-1 induces the expression of about 1000 antisense transcripts from the human host cell genome. A subset of these is also activated by the closely related varicella zoster virus. Antisense transcripts originate either at gene promoters or within the gene body, and they show different susceptibility to the inhibition of early and immediate early viral gene expression. Overexpression of the major viral transcription factor ICP4 is sufficient to turn on a subset of antisense transcripts. Histone marks around transcription start sites of HSV-1-induced and constitutively transcribed antisense transcripts are highly similar, indicating that the genetic loci are already poised to transcribe these novel RNAs. Furthermore, an antisense transcript overlapping with the BBC3 gene (also known as PUMA) transcriptionally silences this potent inducer of apoptosis in cis.

Conclusions

We show for the first time that a virus induces widespread antisense transcription of the host cell genome. We provide evidence that HSV-1 uses this to downregulate a strong inducer of apoptosis. Our findings open new perspectives on global and specific alterations of host cell transcription by viruses.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1329-5) contains supplementary material, which is available to authorized users.

MeCP2 and CTCF: enhancing the cross-talk of silencers

This paper provides a brief introductory review of the most recent advances in our knowledge about the structural and functional aspects of two transcriptional regulators: MeCP2, a protein whose mutated forms are involved in Rett syndrome; and CTCF, a constitutive transcriptional insulator. This is followed by a description of the PTMs affecting these two proteins and an analysis of their known interacting partners. A special emphasis is placed on the recent studies connecting these two proteins, focusing on the still poorly understood potential structural and functional interactions between the two of them on the chromatin substrate. An overview is provided for some of the currently known genes that are dually regulated by these two proteins. Finally, a model is put forward to account for their possible involvement in their regulation of gene expression.

Genome-wide mapping of transcriptional enhancer candidates using DNA and chromatin features in maize

Background

While most cells in multicellular organisms carry the same genetic information, in each cell type only a subset of genes is being transcribed. Such differentiation in gene expression depends, for a large part, on the activation and repression of regulatory sequences, including transcriptional enhancers. Transcriptional enhancers can be located tens of kilobases from their target genes, but display characteristic chromatin and DNA features, allowing their identification by genome-wide profiling. Here we show that integration of chromatin characteristics can be applied to predict distal enhancer candidates in Zea mays, thereby providing a basis for a better understanding of gene regulation in this important crop plant.

Result

To predict transcriptional enhancers in the crop plant maize (Zea mays L. ssp. mays), we integrated available genome-wide DNA methylation data with newly generated maps for chromatin accessibility and histone 3 lysine 9 acetylation (H3K9ac) enrichment in young seedling and husk tissue. Approximately 1500 intergenic regions, displaying low DNA methylation, high chromatin accessibility and H3K9ac enrichment, were classified as enhancer candidates. Based on their chromatin profiles, candidate sequences can be classified into four subcategories. Tissue-specificity of enhancer candidates is defined based on the tissues in which they are identified and putative target genes are assigned based on tissue-specific expression patterns of flanking genes.

Conclusions

Our method identifies three previously identified distal enhancers in maize, validating the new set of enhancer candidates and enlarging the toolbox for the functional characterization of gene regulation in the highly repetitive maize genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1273-4) contains supplementary material, which is available to authorized users.

miR-34a Targets HDAC1-Regulated H3K9 Acetylation on Lipid Accumulation Induced by Homocysteine in Foam Cells

Hyperhomocysteinemia (HHcy) promotes atherogenesis by modification of histone acetylation patterns and regulation of miRNA expression while the underlying molecular mechanisms are not well known. In this study, we investigated the effects of homocysteine (Hcy) on the expression of histone deacetylase 1 (HDAC1) and tested our hypothesis that Hcy-induced atherosclerosis is mediated by increased HDAC1 expression, which is regulated by miR-34a. The expression of HDAC1 increased and acetylation of histone H3 at lysine 9 (H3K9ac) decreased in the aorta of ApoE^-/- mice fed with high methionine diet, whereas miR-34a expression was inhibited. Over-expression of HDAC1 inhibited H3K9ac level and promoted the accumulation of total cholesterol, free cholesterol, and triglycerides in the foam cells. Furthermore, up-regulation of miR-34a reduced HDAC1 expression and inhibited the accumulation of total cholesterol (TC), free cholesterol (FC), and triglycerides (TG) in the foam cells. These data suggest that HDAC1-related H3K9ac plays a key role in Hcy-mediated lipid metabolism disorders, and that miR-34a may be a novel therapeutic target in Hcy-related atherosclerosis. J. Cell. Biochem. 118: 4617-4627, 2017. © 2017 Wiley Periodicals, Inc.

Histone modifications facilitate the coexpression of bidirectional promoters in rice

Background

Bidirectional gene pairs are highly abundant and mostly co-regulated in eukaryotic genomes. The structural features of bidirectional promoters (BDPs) have been well studied in yeast, humans and plants. However, the underlying mechanisms responsible for the coexpression of BDPs remain understudied, especially in plants.

Results

Here, we characterized chromatin features associated with rice BDPs. Several unique chromatin features were present in rice BDPs but were missing from unidirectional promoters (UDPs), including overrepresented active histone marks, canonical nucleosomes and underrepresented H3K27me3. In particular, overrepresented active marks (H3K4ac, H4K12ac, H4K16ac, H3K4me2 and H3K36me3) were truly overrepresented in type I BDPs but not in the other two BDPs, based on a Kolmogorov-Smirnov test.

Conclusions

Our analyses indicate that active marks (H3K4ac, H4K12ac, H4K16ac, H3K4me3, H3K9ac and H3K27ac) may coordinate with repressive marks (H3K27me3 and H3K9me1/3) to build a unique chromatin structure that favors the coregulation of bidirectional gene pairs. Thus, our findings help to enhance the understanding of unique epigenetic mechanisms that regulate bidirectional gene pairs and may improve the manipulation of gene pairs for crop bioengineering.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3125-0) contains supplementary material, which is available to authorized users.

CTCF: making the right connections

Ghirlando and Felsenfeld review recent major advances in understanding the multiple roles of CTCF in gene regulation and genome organization and especially in how such domains are generated and organized.

The role of the zinc finger protein CTCF in organizing the genome within the nucleus is now well established. Widely separated sites on DNA, occupied by both CTCF and the cohesin complex, make physical contacts that create large loop domains. Additional contacts between loci within those domains, often also mediated by CTCF, tend to be favored over contacts between loci in different domains. A large number of studies during the past 2 years have addressed the questions: How are these loops generated? What are the effects of disrupting them? Are there rules governing large-scale genome organization? It now appears that the strongest and evolutionarily most conserved of these CTCF interactions have specific rules for the orientation of the paired CTCF sites, implying the existence of a nonequilibrium mechanism of generation. Recent experiments that invert, delete, or inactivate one of a mating CTCF pair result in major changes in patterns of organization and gene expression in the surrounding regions. What remain to be determined are the detailed molecular mechanisms for re-establishing loop domains and maintaining them after replication and mitosis. As recently published data show, some mechanisms may involve interactions with noncoding RNAs as well as protein cofactors. Many CTCF sites are also involved in other functions such as modulation of RNA splicing and specific regulation of gene expression, and the relationship between these activities and loop formation is another unanswered question that should keep investigators occupied for some time.

Functional characterization of open chromatin in bidirectional promoters of rice

Bidirectional gene pairs tend to be highly coregulated and function in similar biological processes in eukaryotic genomes. Structural features and functional consequences of bidirectional promoters (BDPs) have received considerable attention among diverse species. However, the underlying mechanisms responsible for the bidirectional transcription and coexpression of BDPs remain poorly understood in plants. In this study, we integrated DNase-seq, RNA-seq, ChIP-seq and MNase-seq data and investigated the effect of physical DNase I hypersensitive site (DHS) positions on the transcription of rice BDPs. We found that the physical position of a DHS relative to the TSS of bidirectional gene pairs can affect the expression of the corresponding genes: the closer a DHS is to the TSS, the higher is the expression level of the genes. Most importantly, we observed that the distribution of DHSs plays a significant role in the regulation of transcription and the coexpression of gene pairs, which are possibly mediated by orchestrating the positioning of histone marks and canonical nucleosomes around BDPs. Our results demonstrate that the combined actions of chromatin structures with DHSs, which contain functional cis-elements for interaction with transcriptional machinery, may play an important role in the regulation of the bidirectional transcription or coexpression in rice BDPs. Our findings may help to enhance the understanding of DHSs in the regulation of bidirectional gene pairs.

The Determinants of Directionality in Transcriptional Initiation

A new paradigm has emerged in recent years characterizing transcription initiation as a bidirectional process encompassing a larger proportion of the genome than previously thought. Past concepts of coding genes thinly scattered among a vast background of transcriptionally inert noncoding DNA have been abandoned. A richer picture has taken shape, integrating transcription of coding genes, enhancer RNAs (eRNAs), and various other noncoding transcriptional events. In this review we give an overview of recent studies detailing the mechanisms of RNA polymerase II (RNA Pol II)-based transcriptional initiation and discuss the ways in which transcriptional direction is established as well as its functional implications.