GRO-seq, A Tool for Identification of Transcripts Regulating Gene Expression

Herpes simplex virus 1 immediate early transcription initiation, pause-release, elongation, and termination in the presence and absence of ICP4

IMPORTANCE

Infection with herpes simplex virus 1 (HSV-1) leads to lifelong infection due to the virus's remarkable ability to control transcription of its own genome, resulting in two transcriptional programs: lytic (highly active) and latent (restricted). The lytic program requires immediate early (IE) proteins to first repress transcription of late viral genes, which then undergo sequential de-repression, leading to a specific sequence of gene expression. Here, we show that the IE ICP4 functions to regulate the cascade by limiting RNA polymerase initiation at immediate early times. However, late viral genes that initiate too early in the absence of ICP4 do not yield mRNA as transcription stalls within gene bodies. It follows that other regulatory steps intercede to prevent elongation of genes at the incorrect time, demonstrating the precise control HSV-1 exerts over its own transcription.

Experimental Validation and Prediction of Super-Enhancers: Advances and Challenges

Super-enhancers (SEs) are cis-regulatory elements of the human genome that have been widely discussed since the discovery and origin of the term. Super-enhancers have been shown to be strongly associated with the expression of genes crucial for cell differentiation, cell stability maintenance, and tumorigenesis. Our goal was to systematize research studies dedicated to the investigation of structure and functions of super-enhancers as well as to define further perspectives of the field in various applications, such as drug development and clinical use. We overviewed the fundamental studies which provided experimental data on various pathologies and their associations with particular super-enhancers. The analysis of mainstream approaches for SE search and prediction allowed us to accumulate existing data and propose directions for further algorithmic improvements of SEs' reliability levels and efficiency. Thus, here we provide the description of the most robust algorithms such as ROSE, imPROSE, and DEEPSEN and suggest their further use for various research and development tasks. The most promising research direction, which is based on topic and number of published studies, are cancer-associated super-enhancers and prospective SE-targeted therapy strategies, most of which are discussed in this review.

Tn5 tagments and transposes oligos to single-stranded DNA for strand-specific RNA sequencing

Tn5 transposon tagments double-stranded DNA and RNA/DNA hybrids to generate nucleic acids that are ready to be amplified for high-throughput sequencing. The nucleic acid substrates for the Tn5 transposon must be explored to increase the applications of Tn5. Here, we found that the Tn5 transposon can transpose oligos into the 5' end of single-stranded DNA longer than 140 nucleotides. Based on this property of Tn5, we developed a tagmentation-based and ligation-enabled single-stranded DNA sequencing method called TABLE-seq. Through a series of reaction temperature, time, and enzyme concentration tests, we applied TABLE-seq to strand-specific RNA sequencing, starting with as little as 30 pg of total RNA. Moreover, compared with traditional dUTP-based strand-specific RNA sequencing, this method detects more genes, has a higher strand specificity, and shows more evenly distributed reads across genes. Together, our results provide insights into the properties of Tn5 transposons and expand the applications of Tn5 in cutting-edge sequencing techniques.

Predictive model of transcriptional elongation control identifies trans regulatory factors from chromatin signatures

Promoter-proximal Polymerase II (Pol II) pausing is a key rate-limiting step for gene expression. DNA and RNA-binding trans-acting factors regulating the extent of pausing have been identified. However, we lack a quantitative model of how interactions of these factors determine pausing, therefore the relative importance of implicated factors is unknown. Moreover, previously unknown regulators might exist. Here we address this gap with a machine learning model that accurately predicts the extent of promoter-proximal Pol II pausing from large-scale genome and transcriptome binding maps and gene annotation and sequence composition features. We demonstrate high accuracy and generalizability of the model by validation on an independent cell line which reveals the model's cell line agnostic character. Model interpretation in light of prior knowledge about molecular functions of regulatory factors confirms the interconnection of pausing with other RNA processing steps. Harnessing underlying feature contributions, we assess the relative importance of each factor, quantify their predictive effects and systematically identify previously unknown regulators of pausing. We additionally identify 16 previously unknown 7SK ncRNA interacting RNA-binding proteins predictive of pausing. Our work provides a framework to further our understanding of the regulation of the critical early steps in transcriptional elongation.

Typical Enhancers, Super-Enhancers, and Cancers

Simple Summary

The cancer genome has been exhaustively studied upon the advent of Next-Generation Sequencing technologies. Coding and non-coding sequences have been defined as hotspots of genomic variations that affect the naïve gene expression programs established in normal cells, thus working as endogenous drivers of carcinogenesis. In this review, we comprehensively summarize fundamental aspects of gene expression regulation, with emphasis on the impact of sequence and structural variations mapped across non-coding cis-acting elements of genes encoding for tumor-related transcription factors. Chromatin architecture, epigenome reprogramming, transcriptional enhancers and Super-enhancers, oncogene regulation, cutting-edge technologies, and pharmacological treatment are substantially highlighted.

Abstract

Non-coding segments of the human genome are enriched in cis-regulatory modules that constitute functional elements, such as transcriptional enhancers and Super-enhancers. A hallmark of cancer pathogenesis is the dramatic dysregulation of the “archetype” gene expression profiles of normal human cells. Genomic variations can promote such deficiencies when occurring across enhancers and Super-enhancers, since they affect their mechanistic principles, their functional capacity and specificity, and the epigenomic features of the chromatin microenvironment across which these regulatory elements reside. Here, we comprehensively describe: fundamental mechanisms of gene expression dysregulation in cancers that involve genomic abnormalities within enhancers’ and Super-enhancers’ (SEs) sequences, which alter the expression of oncogenic transcription factors (TFs); cutting-edge technologies applied for the analysis of variation-enriched hotspots of the cancer genome; and pharmacological approaches for the treatment of Super-enhancers’ aberrant function. Finally, we provide an intratumor meta-analysis, which highlights that genomic variations in transcription-factor-driven tumors are accompanied overexpression of genes, a portion of which encodes for additional cancer-related transcription factors.

Which field of research would Gregor Mendel choose in the 21st century?

Gregor Mendel's work on segregation of traits in plants established the basic methodology and rules of genetics. The interruption of Mendel's research activities in 1870 impeded the immediate recognition of the value of his work until the dawn of the 20th century. Only then were his founding laws of genetics validated, propelling the development of biological research toward the birth of molecular biology in the second half of the 20th century. While molecular plant genetics can be viewed as the spiritual heir of Mendel's research, one might wonder whether in the 21st century Gregor Mendel would prefer to practice scientific approaches other than molecular genetics such as population genetics, comparative genomics, or the emerging field of evo-chromo. In this perspective, I review aspects of these fields that might have attracted or perplexed a 21st century Mendel.

Global-run on sequencing identifies Gm11967 as an Akt-dependent long noncoding RNA involved in insulin sensitivity

The insulin responsive Akt and FoxO1 signaling axis is a key regulator of the hepatic transcriptional response to nutrient intake. Here, we used global run-on sequencing (GRO-seq) to measure the nascent transcriptional response to fasting and refeeding as well as define the specific role of hepatic Akt and FoxO1 signaling in mediating this response. We identified 599 feeding-regulated transcripts, as well as over 6,000 eRNAs, and mapped their dependency on Akt and FoxO1 signaling. Further, we identified several feeding-regulated lncRNAs, including the lncRNA Gm11967, whose expression was dependent upon the liver Akt-FoxO1 axis. Restoring Gm11967 expression in mice lacking liver Akt improved insulin sensitivity and induced glucokinase protein expression, indicating that Akt-dependent control of Gm11967 contributes to the translational control of glucokinase. More broadly, we have generated a unique genome-wide dataset that defines the feeding and Akt/FoxO1-dependent transcriptional changes in response to nutrient availability.

SMYD5 catalyzes histone H3 lysine 36 trimethylation at promoters

Histone marks, carriers of epigenetic information, regulate gene expression. In mammalian cells, H3K36me3 is mainly catalyzed by SETD2 at gene body regions. Here, we find that in addition to gene body regions, H3K36me3 is enriched at promoters in primary cells. Through screening, we identify SMYD5, which is recruited to chromatin by RNA polymerase II, as a methyltransferase catalyzing H3K36me3 at promoters. The enzymatic activity of SMYD5 is dependent on its C-terminal glutamic acid-rich domain. Overexpression of full-length Smyd5, but not the C-terminal domain-truncated Smyd5, restores H3K36me3 at promoters in Smyd5 knockout cells. Furthermore, elevated Smyd5 expression contributes to tumorigenesis in liver hepatocellular carcinoma. Together, our findings identify SMYD5 as the H3K36me3 methyltransferase at promoters that regulates gene expression, providing insights into the localization and function of H3K36me3.

ICP22 of Herpes Simplex Virus 1 Decreases RNA Polymerase Processivity

To determine the role of ICP22 in transcription, we performed precise nuclear run-on followed by deep sequencing (PRO-seq) and global nuclear run-on with sequencing (GRO-seq) in cells infected with a viral mutant lacking the entire ICP22-encoding α22 (US1/US1.5) gene and a virus derived from this mutant bearing a restored α22 gene. At 3 h postinfection (hpi), the lack of ICP22 reduced RNA polymerase (Pol) promoter proximal pausing (PPP) on the immediate early α4, α0, and α27 genes. Diminished PPP at these sites accompanied increased Pol processivity across the entire herpes simplex virus 1 (HSV-1) genome in GRO-seq assays, resulting in substantial increases in antisense and intergenic transcription. The diminished PPP on α gene promoters at 3 hpi was distinguishable from effects caused by treatment with a viral DNA polymerase inhibitor at this time. The ICP22 mutant had multiple defects at 6 hpi, including lower viral DNA replication, reduced Pol activity on viral genes, and increased Pol activity on cellular genes. The lack of ICP22 also increased PPP release from most cellular genes, while a minority of cellular genes exhibited decreased PPP release. Taken together, these data indicate that ICP22 acts to negatively regulate transcriptional elongation on viral genes in part to limit antisense and intergenic transcription on the highly compact viral genome. This regulatory function directly or indirectly helps to retain Pol activity on the viral genome later in infection. IMPORTANCE The longstanding observation that ICP22 reduces RNA polymerase II (Pol II) serine 2 phosphorylation, which initiates transcriptional elongation, is puzzling because this phosphorylation is essential for viral replication. The current study helps explain this apparent paradox because it demonstrates significant advantages in negatively regulating transcriptional elongation, including the reduction of antisense and intergenic transcription. Delays in elongation would be expected to facilitate the ordered assembly and functions of transcriptional initiation, elongation, and termination complexes. Such limiting functions are likely to be important in herpesvirus genomes that are otherwise highly transcriptionally active and compact, comprising mostly short, intronless genes near neighboring genes of opposite sense and containing numerous 3'-nested sets of genes that share transcriptional termination signals but differ at transcriptional start sites on the same template strand.

Decoding non-canonical mRNA decay by the endoplasmic-reticulum stress sensor IRE1α

Inositol requiring enzyme 1 (IRE1) mitigates endoplasmic-reticulum (ER) stress by orchestrating the unfolded-protein response (UPR). IRE1 spans the ER membrane, and signals through a cytosolic kinase-endoribonuclease module. The endoribonuclease generates the transcription factor XBP1s by intron excision between similar RNA stem-loop endomotifs, and depletes select cellular mRNAs through regulated IRE1-dependent decay (RIDD). Paradoxically, in mammals RIDD seems to target only mRNAs with XBP1-like endomotifs, while in flies RIDD exhibits little sequence restriction. By comparing nascent and total IRE1α-controlled mRNAs in human cells, we identify not only canonical endomotif-containing RIDD substrates, but also targets without such motifs-degraded by a process we coin RIDDLE, for RIDD lacking endomotif. IRE1α displays two basic endoribonuclease modalities: highly specific, endomotif-directed cleavage, minimally requiring dimers; and more promiscuous, endomotif-independent processing, requiring phospho-oligomers. An oligomer-deficient IRE1α mutant fails to support RIDDLE in vitro and in cells. Our results advance current mechanistic understanding of the UPR.

Enhancers as potential targets for engineering salinity stress tolerance in crop plants

Enhancers represent noncoding regulatory regions of the genome located distantly from their target genes. They regulate gene expression programs in a context-specific manner via interacting with promoters of one or more target genes and are generally associated with transcription factor binding sites and epi(genomic)/chromatin features, such as regions of chromatin accessibility and histone modifications. The enhancers are difficult to identify due to the modularity of their associated features. Although enhancers have been studied extensively in human and animals, only a handful of them has been identified in few plant species till date due to nonavailability of plant-specific experimental and computational approaches for their discovery. Being an important regulatory component of the genome, enhancers represent potential targets for engineering agronomic traits, including salinity stress tolerance in plants. Here, we provide a review of the available experimental and computational approaches along with the associated sequence and chromatin/epigenetic features for the discovery of enhancers in plants. In addition, we provide insights into the challenges and future prospects of enhancer research in plant biology with emphasis on potential applications in engineering salinity stress tolerance in crop plants.

A Transcription Start Site Map in Human Pancreatic Islets Reveals Functional Regulatory Signatures

Identifying the tissue-specific molecular signatures of active regulatory elements is critical to understand gene regulatory mechanisms. Here, we identify transcription start sites (TSS) using cap analysis of gene expression (CAGE) across 57 human pancreatic islet samples. We identify 9,954 reproducible CAGE tag clusters (TCs), ∼20% of which are islet specific and occur mostly distal to known gene TSS. We integrated islet CAGE data with histone modification and chromatin accessibility profiles to identify epigenomic signatures of transcription initiation. Using a massively parallel reporter assay, we validated the transcriptional enhancer activity for 2,279 of 3,378 (∼68%) tested islet CAGE elements (5% false discovery rate). TCs within accessible enhancers show higher enrichment to overlap type 2 diabetes genome-wide association study (GWAS) signals than existing islet annotations, which emphasizes the utility of mapping CAGE profiles in disease-relevant tissue. This work provides a high-resolution map of transcriptional initiation in human pancreatic islets with utility for dissecting active enhancers at GWAS loci.

An integrated comparative physiology and molecular approach pinpoints mediators of breath-hold capacity in dolphins

Background and objectives

Ischemic events, such as ischemic heart disease and stroke, are the number one cause of death globally. Ischemia prevents blood, carrying essential nutrients and oxygen, from reaching tissues, leading to cell and tissue death, and eventual organ failure. While humans are relatively intolerant to ischemic events, other species, such as marine mammals, have evolved a unique tolerance to chronic ischemia/reperfusion during apneic diving. To identify possible molecular features of an increased tolerance for apnea, we examined changes in gene expression in breath-holding dolphins.

Methodology

Here, we capitalized on the adaptations possesed by bottlenose dolphins (Tursiops truncatus) for diving as a comparative model of ischemic stress and hypoxia tolerance to identify molecular features associated with breath holding. Given that signals in the blood may influence physiological changes during diving, we used RNA-Seq and enzyme assays to examine time-dependent changes in gene expression in the blood of breath-holding dolphins.

Results

We observed time-dependent upregulation of the arachidonate 5-lipoxygenase (ALOX5) gene and increased lipoxygenase activity during breath holding. ALOX5 has been shown to be activated during hypoxia in rodent models, and its metabolites, leukotrienes, induce vasoconstriction.

Conclusions and implications

The upregulation of ALOX5 mRNA occurred within the calculated aerobic dive limit of the species, suggesting that ALOX5 may play a role in the dolphin’s physiological response to diving, particularly in a pro-inflammatory response to ischemia and in promoting vasoconstriction. These observations pinpoint a potential molecular mechanism by which dolphins, and perhaps other marine mammals, respond to the prolonged breath holds associated with diving.

Selective Translation of Low Abundance and Upregulated Transcripts in Halobacterium salinarum

Our findings demonstrate conclusively that low abundance and upregulated transcripts are preferentially translated, potentially by environment-specific translation systems with distinct ribosomal protein composition. We show that a complex interplay of transcriptional and posttranscriptional regulation underlies the conditional and modular regulatory programs that generate ribosomes of distinct protein composition. The modular regulation of ribosomal proteins with other transcription, translation, and metabolic genes is generalizable to bacterial and eukaryotic microbes. These findings are relevant to how microorganisms adapt to unfavorable environments when they transition from active growth to quiescence by generating proteins from upregulated transcripts that are in considerably lower abundance relative to transcripts associated with the previous physiological state. Selective translation of transcripts by distinct ribosomes could form the basis for adaptive evolution to new environments through a modular regulation of the translational systems.

When organisms encounter an unfavorable environment, they transition to a physiologically distinct, quiescent state wherein abundant transcripts from the previous active growth state continue to persist, albeit their active transcription is downregulated. In order to generate proteins for the new quiescent physiological state, we hypothesized that the translation machinery must selectively translate upregulated transcripts in an intracellular milieu crowded with considerably higher abundance transcripts from the previous active growth state. Here, we have analyzed genome-wide changes in the transcriptome (RNA sequencing [RNA-seq]), changes in translational regulation and efficiency by ribosome profiling across all transcripts (ribosome profiling [Ribo-seq]), and protein level changes in assembled ribosomal proteins (sequential window acquisition of all theoretical mass spectra [SWATH-MS]) to investigate the interplay of transcriptional and translational regulation in Halobacterium salinarum as it transitions from active growth to quiescence. We have discovered that interplay of regulatory processes at different levels of information processing generates condition-specific ribosomal complexes to translate preferentially pools of low abundance and upregulated transcripts. Through analysis of the gene regulatory network architecture of H. salinarum, Escherichia coli, and Saccharomyces cerevisiae, we demonstrate that this conditional, modular organization of regulatory programs governing translational systems is a generalized feature across all domains of life.

IMPORTANCE Our findings demonstrate conclusively that low abundance and upregulated transcripts are preferentially translated, potentially by environment-specific translation systems with distinct ribosomal protein composition. We show that a complex interplay of transcriptional and posttranscriptional regulation underlies the conditional and modular regulatory programs that generate ribosomes of distinct protein composition. The modular regulation of ribosomal proteins with other transcription, translation, and metabolic genes is generalizable to bacterial and eukaryotic microbes. These findings are relevant to how microorganisms adapt to unfavorable environments when they transition from active growth to quiescence by generating proteins from upregulated transcripts that are in considerably lower abundance relative to transcripts associated with the previous physiological state. Selective translation of transcripts by distinct ribosomes could form the basis for adaptive evolution to new environments through a modular regulation of the translational systems.

Enhancer occlusion transcripts regulate the activity of human enhancer domains via transcriptional interference: a computational perspective

Analysis of ENCODE long RNA-Seq and ChIP-seq (Chromatin Immunoprecipitation Sequencing) datasets for HepG2 and HeLa cell lines uncovered 1647 and 1958 transcripts that interfere with transcription factor binding to human enhancer domains. TFBSs (Transcription Factor Binding Sites) intersected by these 'Enhancer Occlusion Transcripts' (EOTrs) displayed significantly lower relative transcription factor (TF) binding affinities compared to TFBSs for the same TF devoid of EOTrs. Expression of most EOTrs was regulated in a cell line specific manner; analysis for the same TFBSs across cell lines, i.e. in the absence or presence of EOTrs, yielded consistently higher relative TF/DNA-binding affinities for TFBSs devoid of EOTrs. Lower activities of EOTr-associated enhancer domains coincided with reduced occupancy levels for histone tail modifications H3K27ac and H3K9ac. Similarly, the analysis of EOTrs with allele-specific expression identified lower activities for alleles associated with EOTrs. ChIA-PET (Chromatin Interaction Analysis by Paired-End Tag Sequencing) and 5C (Carbon Copy Chromosome Conformation Capture) uncovered that enhancer domains associated with EOTrs preferentially interacted with poised gene promoters. Analysis of EOTr regions with GRO-seq (Global run-on) data established the correlation of RNA polymerase pausing and occlusion of TF-binding. Our results implied that EOTr expression regulates human enhancer domains via transcriptional interference.

Bridging the gap between theory and practice in elucidating modular gene regulatory sequence organisation within genomes

Changes to promoter regions probably have been responsible for many morphological evolutionary transitions, especially in animals. This idea is becoming testable, as data from genome projects amass and enable bioinformaticians to conduct comparative sequence analyses and test for correlations between genotypic similarities or differences and phenotypic likeness or disparity. Although such practical pursuits have initiated some theoretical considerations, a conceptual framework for understanding promoter region evolution, potentially effecting morphological evolution, is only starting to emerge, predominantly resulting from computational research. We contribute to this framework by specifying three big problems for promoter region research; reviewing computational research on promoter region evolution; and exemplifying a topic for future promoter region research - module evolution.

Interpreting and integrating big data in non-coding RNA research

In the last two decades, we have witnessed an impressive crescendo of non-coding RNA studies, due to both the development of high-throughput RNA-sequencing strategies and an ever-increasing awareness of the involvement of newly discovered ncRNA classes in complex regulatory networks. Together with excitement for the possibility to explore previously unknown layers of gene regulation, these advancements led to the realization of the need for shared criteria of data collection and analysis and for novel integrative perspectives and tools aimed at making biological sense of very large bodies of molecular information. In the last few years, efforts to respond to this need have been devoted mainly to the regulatory interactions involving ncRNAs as direct or indirect regulators of protein-coding mRNAs. Such efforts resulted in the development of new computational tools, allowing the exploitation of the information spread in numerous different ncRNA data sets to interpret transcriptome changes under physiological and pathological cell responses. While experimental validation remains essential to identify key RNA regulatory interactions, the integration of ncRNA big data, in combination with systematic literature mining, is proving to be invaluable in identifying potential new players, biomarkers and therapeutic targets in cancer and other diseases.

An Overview of Methodologies in Studying lncRNAs in the High-Throughput Era: When Acronyms ATTACK!

The discovery of pervasive transcription in eukaryotic genomes provided one of many surprising (and perhaps most surprising) findings of the genomic era and led to the uncovering of a large number of previously unstudied transcriptional events. This pervasive transcription leads to the production of large numbers of noncoding RNAs (ncRNAs) and thus opened the window to study these diverse, abundant transcripts of unclear relevance and unknown function. Since that discovery, recent advances in high-throughput sequencing technologies have identified a large collection of ncRNAs, from microRNAs to long noncoding RNAs (lncRNAs). Subsequent discoveries have shown that many lncRNAs play important roles in various eukaryotic processes; these discoveries have profoundly altered our understanding of the regulation of eukaryotic gene expression. Although the identification of ncRNAs has become a standard experimental approach, the functional characterization of these diverse ncRNAs remains a major challenge. In this chapter, we highlight recent progress in the methods to identify lncRNAs and the techniques to study the molecular function of these lncRNAs and the application of these techniques to the study of plant lncRNAs.

Thinking BIG rheumatology: how to make functional genomics data work for you

High-throughput sequencing assays have become an increasingly common part of biological research across multiple fields. Even as the resulting sequences pile up in public databases, it is not always obvious how to make use of these data sets. Functional genomics offers approaches to integrate these "big" data into our understanding of rheumatic diseases. This review aims to provide a primer on thinking about big data from functional genomics in the context of rheumatology, using examples from the field's literature as well as the author's own work to illustrate the execution of functional genomics research. Study design is crucial to ensure the right samples are used to address the question of interest. In addition, sequencing assays produce a variety of data types, from gene expression to 3D chromatin structure and single-cell technologies, that can be integrated into a model of the underlying gene regulatory networks. The best approach for this analysis uses the scientific process: bioinformatic methods should be used in an iterative, hypothesis-driven manner to uncover the disease mechanism. Finally, the future of functional genomics will see big data fully integrated into rheumatology, leading to computationally trained researchers and interactive databases. The goal of this review is not to provide a manual, but to enhance the familiarity of readers with functional genomic approaches and provide a better sense of the challenges and possibilities.