Dynamic reorganization of the AC16 cardiomyocyte transcriptome in response to TNFα signaling revealed by integrated genomic analyses

Relaxin Modulates the Genomic Actions and Biological Effects of Estrogen in the Myometrium by Reducing Estrogen Receptor Alpha Phosphorylation

Estradiol (E2) and relaxin (Rln) are steroid and polypeptide hormones, respectively, with important roles in the female reproductive tract, including myometrium. Some actions of Rln, which are mediated by its membrane receptor RXFP1, require or are augmented by E2 signaling through its cognate nuclear steroid receptor, estrogen receptor alpha (ERα). In contrast, other actions of Rln act in opposition to the effects of E2. Here we explore the molecular and genomic mechanisms that underlie the functional interplay between E2 and Rln in the myometrium. We used both ovariectomized female mice and immortalized human myometrial cells expressing wild type or mutant ERα (hTERT-HM-ERα cells). Our results indicate that Rln attenuates the genomic actions and biological effects of estrogen in the myometrium and myometrial cells by reducing phosphorylation ERα on serine 118 (S118). Interestingly, we observed a potent inhibitory effect of Rln on the E2-dependent binding of ERα across the genome. The reduction in ERα binding was associated with changes in the hormone-regulated transcriptome, including a decrease in the E2-dependent expression of neighboring genes. The inhibitory effects of Rln cotreatment on the E2-dependent phosphorylation of ERα required the nuclear dual-specificity phosphatases DUSP1 and DUSP5. Moreover, the inhibitory effects of Rln were reflected in a concomitant inhibition of the E2-dependent contraction of myometrial cells. Collectively, our results identify a pathway that integrates Rln/RXFP1 and E2/ERα signaling, resulting in a convergence of membrane and nuclear signaling pathways to control genomic and biological outcomes.

A transcription factor (TF) inference method that broadly measures TF activity and identifies mechanistically distinct TF networks

TF profiler is a method of inferring transcription factor regulatory activity, i.e. when a TF is present and actively regulating transcription, directly directly from nascent sequencing assays such as PRO-seq and GRO-seq. Transcription factors orchestrate transcription and play a critical role in cellular maintenance, identity and response to external stimuli. While ChIP assays have measured DNA localization, they fall short of identifying when and where transcription factors are actively regulating transcription. Our method, on the other hand, uses RNA polymerase activity to infer TF activity across hundreds of data sets and transcription factors. Based on these classifications we identify three distinct classes of transcription factors: ubiquitous factors that play roles in cellular homeostasis, driving basal gene programs across tissues and cell types, tissue specific factors that act almost exclusively at enhancers and are themselves regulated at transcription, and stimulus responsive TFs which are regulated post-transcriptionally but act predominantly at enhancers. TF profiler is broadly applicable, providing regulatory insights on any PRO-seq sample for any transcription factor with a known binding motif.

Application of nucleoside or nucleotide analogues in RNA dynamics and RNA-binding protein analysis

Cellular RNAs undergo dynamic changes during RNA biological processes, which are tightly orchestrated by RNA-binding proteins (RBPs). Yet, the investigation of RNA dynamics is hurdled by highly abundant steady-state RNAs, which make the signals of dynamic RNAs less detectable. Notably, the exert of nucleoside or nucleotide analogue-based RNA technologies has provided a remarkable platform for RNA dynamics research, revealing diverse unnoticed features in RNA metabolism. In this review, we focus on the application of two types of analogue-based RNA sequencing, antigen-/antibody- and click chemistry-based methodologies, and summarize the RNA dynamics features revealed. Moreover, we discuss emerging single-cell newly transcribed RNA sequencing methodologies based on nucleoside analogue labeling, which provides novel insights into RNA dynamics regulation at single-cell resolution. On the other hand, we also emphasize the identification of RBPs that interact with polyA, non-polyA RNAs, or newly transcribed RNAs and also their associated RNA-binding domains at genomewide level through ultraviolet crosslinking and mass spectrometry in different contexts. We anticipated that further modification and development of these analogue-based RNA and RBP capture technologies will aid in obtaining an unprecedented understanding of RNA biology. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry RNA Methods > RNA Analyses in Cells.

The Effect of Daily Methylsulfonylmethane (MSM) Consumption on High-Density Lipoprotein Cholesterol in Healthy Overweight and Obese Adults: A Randomized Controlled Trial

Interventions to decrease inflammation and improve metabolic function hold promise for the prevention of obesity-related diseases. Methylsulfonylmethane (MSM) is a naturally occurring compound that demonstrates antioxidant and anti-inflammatory effects. Improvements in measures of metabolic health have been observed in mouse models of obesity and diabetes following MSM treatment. However, the effects of MSM on obesity-related diseases in humans have not been investigated. Therefore, the purpose of this investigation was to determine whether MSM supplementation improves cardiometabolic health, and markers of inflammation and oxidative status. A randomized, double-blind, placebo-controlled design was utilized with a total of 22 overweight or obese adults completing the study. Participants received either a placebo (white rice flour) or 3 g MSM daily for 16 weeks. Measurements occurred at baseline and after 4, 8, and 16 weeks. Outcome measures included fasting glucose, insulin, blood lipids, blood pressure, body composition, metabolic rate, and markers of inflammation and oxidative status. The primary finding of this work shows that high-density lipoprotein cholesterol was elevated at 8 and 16 weeks of daily MSM consumption compared to baseline, (p = 0.008, p = 0.013). Our findings indicate that MSM supplementation may improve the cholesterol profile by resulting in higher levels of high-density lipoprotein cholesterol.

Genomic Landscapes of Noncoding RNAs Regulating VEGFA and VEGFC Expression in Endothelial Cells

Vascular endothelial growth factors (VEGFs) are best known as key regulators of angiogenesis and lymphangiogenesis. Although VEGFs have been promising therapeutic targets for various cardiovascular diseases, their regulatory landscape in endothelial cells remains elusive. Several studies have highlighted the involvement of noncoding RNAs (ncRNAs) in the modulation of VEGF expression. In this study, we investigated the role of two classes of ncRNAs, long ncRNAs (lncRNAs) and enhancer RNAs (eRNAs), in the transcriptional regulation of VEGFA and VEGFC. By integrating genome-wide global run-on sequencing (GRO-Seq) and chromosome conformation capture (Hi-C) data, we identified putative lncRNAs and eRNAs associated with VEGFA and VEGFC genes in endothelial cells. A subset of the identified putative enhancers demonstrated regulatory activity in a reporter assay. Importantly, we demonstrate that deletion of enhancers and lncRNAs by CRISPR/Cas9 promoted significant changes in VEGFA and VEGFC expression. Transcriptome sequencing (RNA-Seq) data from lncRNA deletions showed downstream factors implicated in VEGFA- and VEGFC-linked pathways, such as angiogenesis and lymphangiogenesis, suggesting functional roles for these lncRNAs. Our study uncovers novel lncRNAs and eRNAs regulating VEGFA and VEGFC that can be targeted to modulate the expression of these important molecules in endothelial cells.

Transcription factor enrichment analysis (TFEA) quantifies the activity of multiple transcription factors from a single experiment

Detecting changes in the activity of a transcription factor (TF) in response to a perturbation provides insights into the underlying cellular process. Transcription Factor Enrichment Analysis (TFEA) is a robust and reliable computational method that detects positional motif enrichment associated with changes in transcription observed in response to a perturbation. TFEA detects positional motif enrichment within a list of ranked regions of interest (ROIs), typically sites of RNA polymerase initiation inferred from regulatory data such as nascent transcription. Therefore, we also introduce muMerge, a statistically principled method of generating a consensus list of ROIs from multiple replicates and conditions. TFEA is broadly applicable to data that informs on transcriptional regulation including nascent transcription (eg. PRO-Seq), CAGE, histone ChIP-Seq, and accessibility data (e.g., ATAC-Seq). TFEA not only identifies the key regulators responding to a perturbation, but also temporally unravels regulatory networks with time series data. Consequently, TFEA serves as a hypothesis-generating tool that provides an easy, rigorous, and cost-effective means to broadly assess TF activity yielding new biological insights.

Quantitative Proteomic Profiling of Mitochondrial Toxicants in a Human Cardiomyocyte Cell Line

Mitochondria are essential cellular organelles that participate in important cellular processes, including bioenergetics, metabolism, and signaling. Recent functional and proteomic studies have revealed the remarkable complexity of mitochondrial protein organization. Mitochondrial protein machineries with diverse functions such as protein translocation, respiration, metabolite transport, protein quality control and the control of membrane architecture interact with each other in dynamic networks. The goal of this study was to identify protein expression changes in a human cardiomyocyte cell line treated with several mitochondrial toxicants which inhibit mitochondrial membrane potential (MMP) and mitochondrial respiration. AC16 human cardiomyocyte cells were treated with carbonyl cyanide p-(trifluoromethoxy)phenylhydrazone (FCCP), dinoterb, picoxystrobin, pinacyanol, and triclocarban for 18 h around the IC₅₀ values generated from MMP assay. The samples were harvested and labeled with tandem mass tags with different mass isotopes. Peptide assignment was performed in Proteome Discoverer. Each dataset was analyzed in Ingenuity Pathway Analysis (IPA). In the proteomic profile, these compounds showed dysregulation of a group of mitochondrial proteins (e.g., NDUA, NDUB, BCS1, CYB5B, and SDHF2), as well as proteins involved in lipid metabolism (e.g., CPT, MECR, and LPGAT1), cytoskeleton protein changes (e.g., CROCC, LAMC3, FBLN1, and FMN2) and stress response (e.g., IKBKG, IKBB, SYVN1, SOD2, and CPIN1). Proteomic data from the current study provides key insights into chemical induced cellular pathway dysregulation, supporting the use of proteomic profiling as a sensitive method to further explore molecular functions and disease pathogenesis upon exposure to environmental chemicals.

Characterization of basal and estrogen-regulated antisense transcription in breast cancer cells: Role in regulating sense transcription

Estrogen-responsive breast cancer cells exhibit both basal and estrogen-regulated transcriptional programs, which lead to the transcription of many different transcription units (i.e., genes), including those that produce coding and non-coding sense (e.g., mRNA, lncRNA) and antisense (i.e., asRNA) transcripts. We have previously characterized the global basal and estrogen-regulated transcriptomes in estrogen receptor alpha (ERα)-positive MCF-7 breast cancer cells. Herein, we have mined genomic data to define three classes of antisense transcription in MCF-7 cells based on where their antisense transcription termination sites reside relative to their cognate sense mRNA and lncRNA genes. These three classes differ in their response to estrogen treatment, the enrichment of a number of genomic features associated with active promoters (H3K4me3, RNA polymerase II, open chromatin architecture), and the biological functions of their cognate sense genes as analyzed by DAVID gene ontology. We further characterized two estrogen-regulated antisense transcripts arising from the MYC gene in MCF-7 cells, showing that these antisense transcripts are 5'-capped, 3'-polyadenylated, and localized to different compartments of the cell. Together, our analyses have revealed distinct classes of antisense transcription correlated to different biological processes and response to estrogen stimulation, uncovering another layer of hormone-regulated gene regulation.

Nucleoporins in cardiovascular disease

Cardiovascular disease is a pressing health problem with significant global health, societal, and financial burdens. Understanding the molecular basis of polygenic cardiac pathology is thus essential to devising novel approaches for management and treatment. Recent identification of uncharacterized regulatory functions for a class of nuclear envelope proteins called nucleoporins offers the opportunity to understand novel putative mechanisms of cardiac disease development and progression. Consistent reports of nucleoporin deregulation associated with ischemic and dilated cardiomyopathies, arrhythmias and valvular disorders suggests that nucleoporin impairment may be a significant but understudied variable in cardiopathologic disorders. This review discusses and converges existing literature regarding nuclear pore complex proteins and their association with cardiac pathologies, and proposes a role for nucleoporins as facilitators of cardiac disease.

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.

The Toxoplasma effector TEEGR promotes parasite persistence by modulating NF-κB signalling via EZH2

The protozoan parasite Toxoplasma gondii has co-evolved with its homeothermic hosts (humans included) strategies that drive its quasi-asymptomatic persistence in hosts, hence optimizing the chance of transmission to new hosts. Persistence, which starts with a small subset of parasites that escape host immune killing and colonize the so-called immune privileged tissues where they differentiate into a low replicating stage, is driven by the interleukin 12 (IL-12)-interferon-γ (IFN-γ) axis. Recent characterization of a family of Toxoplasma effectors that are delivered into the host cell, in which they rewire the host cell gene expression, has allowed the identification of regulators of the IL-12-IFN-γ axis, including repressors. We now report on the dense granule-resident effector, called TEEGR (Toxoplasma E2F4-associated EZH2-inducing gene regulator) that counteracts the nuclear factor-κB (NF-κB) signalling pathway. Once exported into the host cell, TEEGR ends up in the nucleus where it not only complexes with the E2F3 and E2F4 host transcription factors to induce gene expression, but also promotes shaping of a non-permissive chromatin through its capacity to switch on EZH2. Remarkably, EZH2 fosters the epigenetic silencing of a subset of NF-κB-regulated cytokines, thereby strongly contributing to the host immune equilibrium that influences the host immune response and promotes parasite persistence in mice.

Chromatin run-on and sequencing maps the transcriptional regulatory landscape of glioblastoma multiforme

The human genome encodes a variety of poorly understood RNA species that remain challenging to identify using existing genomic tools. We developed chromatin run-on and sequencing (ChRO-seq) to map the location of RNA polymerase for almost any input sample, including samples with degraded RNA that are intractable to RNA sequencing. We used ChRO-seq to map nascent transcription in primary human glioblastoma (GBM) brain tumors. Enhancers identified in primary GBMs resemble open chromatin in the normal human brain. Rare enhancers that are activated in malignant tissue drive regulatory programs similar to the developing nervous system. We identified enhancers that regulate groups of genes that are characteristic of each known GBM subtype and transcription factors that drive them. Finally we discovered a core group of transcription factors that control the expression of genes associated with clinical outcomes. This study characterizes the transcriptional landscape of GBM and introduces ChRO-seq as a method to map regulatory programs that contribute to complex diseases.

Metabolic regulation of transcription through compartmentalized NAD+ biosynthesis

NAD⁺ (nicotinamide adenine dinucleotide in its oxidized state) is an essential molecule for a variety of physiological processes. It is synthesized in distinct subcellular compartments by three different synthases (NMNAT-1, -2, and -3). We found that compartmentalized NAD⁺ synthesis by NMNATs integrates glucose metabolism and adipogenic transcription during adipocyte differentiation. Adipogenic signaling rapidly induces cytoplasmic NMNAT-2, which competes with nuclear NMNAT-1 for the common substrate, nicotinamide mononucleotide, leading to a precipitous reduction in nuclear NAD⁺ levels. This inhibits the catalytic activity of poly[adenosine diphosphate (ADP)-ribose] polymerase-1 (PARP-1), a NAD⁺-dependent enzyme that represses adipogenic transcription by ADP-ribosylating the adipogenic transcription factor C/EBPβ. Reversal of PARP-1-mediated repression by NMNAT-2-mediated nuclear NAD⁺ depletion in response to adipogenic signals drives adipogenesis. Thus, compartmentalized NAD⁺ synthesis functions as an integrator of cellular metabolism and signal-dependent transcriptional programs.

Detecting Differential Transcription Factor Activity from ATAC-Seq Data

Transcription factors are managers of the cellular factory, and key components to many diseases. Many non-coding single nucleotide polymorphisms affect transcription factors, either by directly altering the protein or its functional activity at individual binding sites. Here we first briefly summarize high-throughput approaches to studying transcription factor activity. We then demonstrate, using published chromatin accessibility data (specifically ATAC-seq), that the genome-wide profile of TF recognition motifs relative to regions of open chromatin can determine the key transcription factor altered by a perturbation. Our method of determining which TFs are altered by a perturbation is simple, is quick to implement, and can be used when biological samples are limited. In the future, we envision that this method could be applied to determine which TFs show altered activity in response to a wide variety of drugs and diseases.

Enhancer RNA profiling predicts transcription factor activity

Transcription factors (TFs) exert their regulatory influence through the binding of enhancers, resulting in coordination of gene expression programs. Active enhancers are often characterized by the presence of short, unstable transcripts termed enhancer RNAs (eRNAs). While their function remains unclear, we demonstrate that eRNAs are a powerful readout of TF activity. We infer sites of eRNA origination across hundreds of publicly available nascent transcription data sets and show that eRNAs initiate from sites of TF binding. By quantifying the colocalization of TF binding motif instances and eRNA origins, we derive a simple statistic capable of inferring TF activity. In doing so, we uncover dozens of previously unexplored links between diverse stimuli and the TFs they affect.

Dynamic evolution of regulatory element ensembles in primate CD4+ T cells

How evolutionary changes at enhancers affect the transcription of target genes remains an important open question. Previous comparative studies of gene expression have largely measured the abundance of messenger RNA, which is affected by post-transcriptional regulatory processes, hence limiting inferences about the mechanisms underlying expression differences. Here, we directly measured nascent transcription in primate species, allowing us to separate transcription from post-transcriptional regulation. We used precision run-on and sequencing to map RNA polymerases in resting and activated CD4⁺ T cells in multiple human, chimpanzee and rhesus macaque individuals, with rodents as outgroups. We observed general conservation in coding and non-coding transcription, punctuated by numerous differences between species, particularly at distal enhancers and non-coding RNAs. Genes regulated by larger numbers of enhancers are more frequently transcribed at evolutionarily stable levels, despite reduced conservation at individual enhancers. Adaptive nucleotide substitutions are associated with lineage-specific transcription and at one locus, SGPP2, we predict and experimentally validate that multiple substitutions contribute to human-specific transcription. Collectively, our findings suggest a pervasive role for evolutionary compensation across ensembles of enhancers that jointly regulate target genes.

Enhancer transcription reveals subtype-specific gene expression programs controlling breast cancer pathogenesis

Noncoding transcription is a defining feature of active enhancers, linking transcription factor (TF) binding to the molecular mechanisms controlling gene expression. To determine the relationship between enhancer activity and biological outcomes in breast cancers, we profiled the transcriptomes (using GRO-seq and RNA-seq) and epigenomes (using ChIP-seq) of 11 different human breast cancer cell lines representing five major molecular subtypes of breast cancer, as well as two immortalized (“normal”) human breast cell lines. In addition, we developed a robust and unbiased computational pipeline that simultaneously identifies putative subtype-specific enhancers and their cognate TFs by integrating the magnitude of enhancer transcription, TF mRNA expression levels, TF motif P-values, and enrichment of H3K4me1 and H3K27ac. When applied across the 13 different cell lines noted above, the Total Functional Score of Enhancer Elements (TFSEE) identified key breast cancer subtype-specific TFs that act at transcribed enhancers to dictate gene expression patterns determining growth outcomes, including Forkhead TFs, FOSL1, and PLAG1. FOSL1, a Fos family TF, (1) is highly enriched at the enhancers of triple negative breast cancer (TNBC) cells, (2) acts as a key regulator of the proliferation and viability of TNBC cells, but not Luminal A cells, and (3) is associated with a poor prognosis in TNBC breast cancer patients. Taken together, our results validate our enhancer identification pipeline and reveal that enhancers transcribed in breast cancer cells direct critical gene regulatory networks that promote pathogenesis.

PARP-1 Controls the Adipogenic Transcriptional Program by PARylating C/EBPβ and Modulating Its Transcriptional Activity

Poly(ADP-ribosyl)ation (PARylation) is a post-translational modification of proteins mediated by PARP family members, such as PARP-1. Although PARylation has been studied extensively, few examples of definitive biological roles for site-specific PARylation have been reported. Here we show that C/EBPβ, a key pro-adipogenic transcription factor, is PARylated by PARP-1 on three amino acids in a conserved regulatory domain. PARylation at these sites inhibits C/EBPβ's DNA binding and transcriptional activities and attenuates adipogenesis in various genetic and cell-based models. Interestingly, PARP-1 catalytic activity drops precipitously during the first 48 hr of differentiation, corresponding to a release of C/EBPβ from PARylation-mediated inhibition. This promotes the binding of C/EBPβ at enhancers controlling the expression of adipogenic target genes and continued differentiation. Depletion or chemical inhibition of PARP-1, or mutation of the PARylation sites on C/EBPβ, enhances these early adipogenic events. Collectively, our results provide a clear example of how site-specific PARylation drives biological outcomes.

Dynamic assembly and activation of estrogen receptor α enhancers through coregulator switching

Although many features of active transcriptional enhancers have been defined by genomic assays, we lack a clear understanding of the order of events leading to enhancer formation and activation as well as the dynamics of coregulator interactions within the enhancer complex. Here, we used selective loss- or gain-of-function mutants of estrogen receptor α (ERα) to define two distinct phases of ligand-dependent enhancer formation. In the first phase (0-20 min), p300 is recruited to ERα by Mediator as well as p300's acetylhistone-binding bromodomain to promote initial enhancer formation, which is not competent for sustained activation. In the second phase (20-45 min), p300 is recruited to ERα by steroid receptor coregulators (SRCs) for enhancer maturation and maintenance. Successful transition between these two phases ("coregulator switching") is required for proper enhancer function. Failure to recruit p300 during either phase leads to abortive enhancer formation and a lack of target gene expression. Our results reveal an ordered and cooperative assembly of ERα enhancers requiring functional interplay among p300, Mediator, and SRCs, which has implications for hormone-dependent gene regulation in breast cancers. More broadly, our results demonstrate the unexpectedly dynamic nature of coregulator interactions within enhancer complexes, which are likely to be a defining feature of all enhancers.

Genome-Wide Dynamics of Nascent Noncoding RNA Transcription in Porcine Heart After Myocardial Infarction

Background—

Microarrays and RNA sequencing are widely used to profile transcriptome remodeling during myocardial ischemia. However, the steady-state RNA analysis lacks in sensitivity to detect all noncoding RNA species and does not provide separation between transcriptional and post-transcriptional regulations. Here, we provide the first comprehensive analysis of nascent RNA profiles of mRNAs, primary micro-RNAs, long noncoding RNAs, and enhancer RNAs in a large animal model of acute infarction.

Methods and Results—

Acute infarction was induced by cardiac catheterization of domestic swine. Nuclei isolated from healthy, border zone, and ischemic regions of the affected heart were subjected to global run-on sequencing. Global run-on sequencing analysis indicated that half of affected genes are regulated at the level of transcriptional pausing. A gradient of induction of inflammatory mediators and repression of peroxisome proliferator-activated receptor signaling and oxidative phosphorylation was detected when moving from healthy toward infarcted area. In addition, we interrogated the transcriptional regulation of primary micro-RNAs and provide evidence that several arrhythmia-related target genes exhibit repression at post-transcriptional level. We identified 450 long noncoding RNAs differently regulated by ischemia, including novel conserved long noncoding RNAs expressed in antisense orientation to myocardial transcription factors GATA-binding protein 4, GATA-binding protein 6, and Krüppel-like factor 6. Finally, characterization of enhancers exhibiting differential expression of enhancer RNAs pointed a central role for Krüppel-like factor, MEF2C, ETS, NFY, ATF, E2F2, and NRF1 transcription factors in determining transcriptional responses to ischemia.

Conclusions—

Global run-on sequencing allowed us to follow the gradient of gene expression occurring in the ischemic heart and identify novel noncoding RNAs regulated by oxygen deprivation. These findings highlight potential new targets for diagnosis and treatment of myocardial ischemia.

Computational Approaches for Mining GRO-Seq Data to Identify and Characterize Active Enhancers

Transcriptional enhancers are DNA regulatory elements that are bound by transcription factors and act to positively regulate the expression of nearby or distally located target genes. Enhancers have many features that have been discovered using genomic analyses. Recent studies have shown that active enhancers recruit RNA polymerase II (Pol II) and are transcribed, producing enhancer RNAs (eRNAs). GRO-seq, a method for identifying the location and orientation of all actively transcribing RNA polymerases across the genome, is a powerful approach for monitoring nascent enhancer transcription. Furthermore, the unique pattern of enhancer transcription can be used to identify enhancers in the absence of any information about the underlying transcription factors. Here, we describe the computational approaches required to identify and analyze active enhancers using GRO-seq data, including data pre-processing, alignment, and transcript calling. In addition, we describe protocols and computational pipelines for mining GRO-seq data to identify active enhancers, as well as known transcription factor binding sites that are transcribed. Furthermore, we discuss approaches for integrating GRO-seq-based enhancer data with other genomic data, including target gene expression and function. Finally, we describe molecular biology assays that can be used to confirm and explore further the function of enhancers that have been identified using genomic assays. Together, these approaches should allow the user to identify and explore the features and biological functions of new cell type-specific enhancers.

Chemical genetic discovery of PARP targets reveals a role for PARP-1 in transcription elongation

Poly[adenosine diphosphate (ADP)-ribose] polymerases (PARPs) are a family of enzymes that modulate diverse biological processes through covalent transfer of ADP-ribose from the oxidized form of nicotinamide adenine dinucleotide (NAD(+)) onto substrate proteins. Here we report a robust NAD(+) analog-sensitive approach for PARPs, which allows PARP-specific ADP-ribosylation of substrates that is suitable for subsequent copper-catalyzed azide-alkyne cycloaddition reactions. Using this approach, we mapped hundreds of sites of ADP-ribosylation for PARPs 1, 2, and 3 across the proteome, as well as thousands of PARP-1-mediated ADP-ribosylation sites across the genome. We found that PARP-1 ADP-ribosylates and inhibits negative elongation factor (NELF), a protein complex that regulates promoter-proximal pausing by RNA polymerase II (Pol II). Depletion or inhibition of PARP-1 or mutation of the ADP-ribosylation sites on NELF-E promotes Pol II pausing, providing a clear functional link between PARP-1, ADP-ribosylation, and NELF. This analog-sensitive approach should be broadly applicable across the PARP family and has the potential to illuminate the ADP-ribosylated proteome and the molecular mechanisms used by individual PARPs to mediate their responses to cellular signals.

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

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

groHMM: a computational tool for identifying unannotated and cell type-specific transcription units from global run-on sequencing data

Background

Global run-on coupled with deep sequencing (GRO-seq) provides extensive information on the location and function of coding and non-coding transcripts, including primary microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and enhancer RNAs (eRNAs), as well as yet undiscovered classes of transcripts. However, few computational tools tailored toward this new type of sequencing data are available, limiting the applicability of GRO-seq data for identifying novel transcription units.

Results

Here, we present groHMM, a computational tool in R, which defines the boundaries of transcription units de novo using a two state hidden-Markov model (HMM). A systematic comparison of the performance between groHMM and two existing peak-calling methods tuned to identify broad regions (SICER and HOMER) favorably supports our approach on existing GRO-seq data from MCF-7 breast cancer cells. To demonstrate the broader utility of our approach, we have used groHMM to annotate a diverse array of transcription units (i.e., primary transcripts) from four GRO-seq data sets derived from cells representing a variety of different human tissue types, including non-transformed cells (cardiomyocytes and lung fibroblasts) and transformed cells (LNCaP and MCF-7 cancer cells), as well as non-mammalian cells (from flies and worms). As an example of the utility of groHMM and its application to questions about the transcriptome, we show how groHMM can be used to analyze cell type-specific enhancers as defined by newly annotated enhancer transcripts.

Conclusions

Our results show that groHMM can reveal new insights into cell type-specific transcription by identifying novel transcription units, and serve as a complete and useful tool for evaluating functional genomic elements in cells.

Electronic supplementary material

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

A novel, dynamic pattern-based analysis of NF-κB binding during the priming phase of liver regeneration reveals switch-like functional regulation of target genes

Following partial hepatectomy, a coordinated series of molecular events occurs to regulate hepatocyte entry into the cell cycle to recover lost mass. In rats during the first 6 h following resection, hepatocytes are primed by a tightly controlled cytokine response to prepare hepatocytes to begin replication. Although it appears to be a critical element driving regeneration, the cytokine response to resection has not yet been fully characterized. Specifically, the role of one of the key response elements to cytokine signaling (NF-κB) remains incompletely characterized. In this study, we present a novel, genome-wide, pattern-based analysis characterizing NF-κB binding during the priming phase of liver regeneration. We interrogated the dynamic regulation of priming by NF-κB through categorizing NF-κB binding in different temporal profiles: immediate sustained response, early transient response, and delayed response to partial hepatectomy. We then identified functional regulation of NF-κB binding by relating the temporal response profile to differential gene expression. We found that NF-κB bound genes govern negative regulation of cell growth and inflammatory response immediately following hepatectomy. NF-κB also transiently regulates genes responsible for lipid biosynthesis and transport as well as induction of apoptosis following hepatectomy. By the end of the priming phase, NF-κB regulation of genes involved in inflammatory response, negative regulation of cell death, and extracellular structure organization became prominent. These results suggest that NF-κB regulates target genes through binding and unbinding in immediate, transient, and delayed patterns. Such dynamic switch-like patterns of NF-κB binding may govern different functional transitions that drive the onset of regeneration.

A unified architecture of transcriptional regulatory elements

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

Identification of active transcriptional regulatory elements from GRO-seq data

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

Long noncoding RNAs in cardiovascular diseases

In recent year, increasing evidence suggests that noncoding RNAs play important roles in the regulation of tissue homeostasis and pathophysiological conditions. Besides small noncoding RNAs (eg, microRNAs), >200-nucleotide long transcripts, namely long noncoding RNAs (lncRNAs), can interfere with gene expressions and signaling pathways at various stages. In the cardiovascular system, studies have detected and characterized the expression of lncRNAs under normal physiological condition and in disease states. Several lncRNAs are regulated during acute myocardial infarction (eg, Novlnc6) and heart failure (eg, Mhrt), whereas others control hypertrophy, mitochondrial function and apoptosis of cardiomyocytes. In the vascular system, the endothelial-expressed lncRNAs (eg, MALAT1 and Tie-1-AS) can regulate vessel growth and function, whereas the smooth-muscle-expressed lncRNA smooth muscle and endothelial cell-enriched migration/differentiation-associated long noncoding RNA was recently shown to control the contractile phenotype of smooth muscle cells. This review article summarizes the data on lncRNA expressions in mouse and human and highlights identified cardiovascular lncRNAs that might play a role in cardiovascular diseases. Although our understanding of lncRNAs is still in its infancy, these examples may provide helpful insights how lncRNAs interfere with cardiovascular diseases.

TNFα signaling exposes latent estrogen receptor binding sites to alter the breast cancer cell transcriptome

The interplay between mitogenic and proinflammatory signaling pathways plays key roles in determining the phenotypes and clinical outcomes of breast cancers. Using GRO-seq in MCF-7 cells, we defined the immediate transcriptional effects of crosstalk between estradiol (E2) and TNFα, identifying a large set of target genes whose expression is rapidly altered with combined E2 + TNFα treatment, but not with either agent alone. The pleiotropic effects on gene transcription in response to E2 + TNFα are orchestrated by extensive remodeling of the ERα enhancer landscape in an NF-κB- and FoxA1-dependent manner. In addition, expression of the de novo and synergistically regulated genes is strongly associated with clinical outcomes in breast cancers. Together, our genomic and molecular analyses indicate that TNFα signaling, acting in pathways culminating in the redistribution of NF-κB and FoxA1 binding sites across the genome, creates latent ERα binding sites that underlie altered patterns of gene expression and clinically relevant cellular responses.

From discovery to function: the expanding roles of long noncoding RNAs in physiology and disease

Long noncoding RNAs (lncRNAs) are a relatively poorly understood class of RNAs with little or no coding capacity transcribed from a set of incompletely annotated genes. They have received considerable attention in the past few years and are emerging as potentially important players in biological regulation. Here we discuss the evolving understanding of this new class of molecular regulators that has emerged from ongoing research, which continues to expand our databases of annotated lncRNAs and provide new insights into their physical properties, molecular mechanisms of action, and biological functions. We outline the current strategies and approaches that have been employed to identify and characterize lncRNAs, which have been instrumental in revealing their multifaceted roles ranging from cis- to trans-regulation of gene expression and from epigenetic modulation in the nucleus to posttranscriptional control in the cytoplasm. In addition, we highlight the molecular and biological functions of some of the best characterized lncRNAs in physiology and disease, especially those relevant to endocrinology, reproduction, metabolism, immunology, neurobiology, muscle biology, and cancer. Finally, we discuss the tremendous diagnostic and therapeutic potential of lncRNAs in cancer and other diseases.

4sUDRB-seq: measuring genomewide transcriptional elongation rates and initiation frequencies within cells

Although transcriptional elongation by RNA polymerase II is coupled with many RNA-related processes, genomewide elongation rates remain unknown. We describe a method, called 4sUDRB-seq, based on reversible inhibition of transcription elongation coupled with tagging newly transcribed RNA with 4-thiouridine and high throughput sequencing to measure simultaneously with high confidence genome-wide transcription elongation rates in cells. We find that most genes are transcribed at about 3.5 Kb/min, with elongation rates varying between 2 Kb/min and 6 Kb/min. 4sUDRB-seq can facilitate genomewide exploration of the involvement of specific elongation factors in transcription and the contribution of deregulated transcription elongation to various pathologies.