Readthrough transcription: How are DoGs made and what do they do?

MYC Binding Near Transcriptional End Sites Regulates Basal Gene Expression, Read-Through Transcription and Intragenic Contacts

The MYC oncoprotein regulates numerous genes involved in cellular processes such as cell cycle and mitochondrial and ribosomal structure and function. This requires heterodimerization with its partner, MAX, and binding to specific promoter and enhancer elements. Here, we show that MYC and MAX also bind near transcriptional end sites (TESs) of over one-sixth of all annotated genes. These interactions are dose-dependent, evolutionarily conserved, stabilize the normally short-lived MYC protein and regulate expression both in concert with and independent of MYC's binding elsewhere. MYC's TES binding occurs in association with other transcription factors, alters the chromatin landscape, increases nuclease susceptibility and can alter transcriptional read-through, particularly in response to certain stresses. MYC-bound TESs can directly contact promoters and may fine-tune gene expression in response to both physiologic and pathologic stimuli. Collectively, these findings support a previously unrecognized role for MYC in regulating transcription and its read-through via direct intragenic contacts between TESs and promoters.

Downstream-of-gene (DoG) transcripts contribute to an imbalance in the cancer cell transcriptome

Downstream-of-gene (DoG) transcripts are an emerging class of noncoding RNAs. However, it remains largely unknown how DoG RNA production is regulated and whether alterations in DoG RNA signatures exist in major cancers. Here, through transcriptomic analyses of matched tumors and nonneoplastic tissues and cancer cell lines, we reveal a comprehensive catalog of DoG RNA signatures. Through separate lines of evidence, we support the biological importance of DoG RNAs in carcinogenesis. First, we show tissue-specific and stage-specific differential expression of DoG RNAs in tumors versus paired normal tissues with their respective host genes involved in tumor-promoting versus tumor-suppressor pathways. Second, we identify that differential DoG RNA expression is associated with poor patient survival. Third, we identify that DoG RNA induction is a consequence of treating colon cancer cells with the topoisomerase I (TOP1) poison camptothecin and following TOP1 depletion. Our results underlie the significance of DoG RNAs and TOP1-dependent regulation of DoG RNAs in diversifying and modulating the cancer transcriptome.

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.

Transcription termination and readthrough in African swine fever virus

Introduction

African swine fever virus (ASFV) is a nucleocytoplasmic large DNA virus (NCLDV) that encodes its own host-like RNA polymerase (RNAP) and factors required to produce mature mRNA. The formation of accurate mRNA 3' ends by ASFV RNAP depends on transcription termination, likely enabled by a combination of sequence motifs and transcription factors, although these are poorly understood. The termination of any RNAP is rarely 100% efficient, and the transcriptional "readthrough" at terminators can generate long mRNAs which may interfere with the expression of downstream genes. ASFV transcriptome analyses reveal a landscape of heterogeneous mRNA 3' termini, likely a combination of bona fide termination sites and the result of mRNA degradation and processing. While short-read sequencing (SRS) like 3' RNA-seq indicates an accumulation of mRNA 3' ends at specific sites, it cannot inform about which promoters and transcription start sites (TSSs) directed their synthesis, i.e., information about the complete and unprocessed mRNAs at nucleotide resolution.

Methods

Here, we report a rigorous analysis of full-length ASFV transcripts using long-read sequencing (LRS). We systematically compared transcription termination sites predicted from SRS 3' RNA-seq with 3' ends mapped by LRS during early and late infection.

Results

Using in-vitro transcription assays, we show that recombinant ASFV RNAP terminates transcription at polyT stretches in the non-template strand, similar to the archaeal RNAP or eukaryotic RNAPIII, unaided by secondary RNA structures or predicted viral termination factors. Our results cement this T-rich motif (U-rich in the RNA) as a universal transcription termination signal in ASFV. Many genes share the usage of the same terminators, while genes can also use a range of terminators to generate transcript isoforms varying enormously in length. A key factor in the latter phenomenon is the highly abundant terminator readthrough we observed, which is more prevalent during late compared with early infection.

Discussion

This indicates that ASFV mRNAs under the control of late gene promoters utilize different termination mechanisms and factors to early promoters and/or that cellular factors influence the viral transcriptome landscape differently during the late stages of infection.

Mechanisms and consequences of mRNA destabilization during viral infections

During viral infection there is dynamic interplay between the virus and the host to regulate gene expression. In many cases, the host induces the expression of antiviral genes to combat infection, while the virus uses "host shut-off" systems to better compete for cellular resources and to limit the induction of the host antiviral response. Viral mechanisms for host shut-off involve targeting translation, altering host RNA processing, and/or inducing the degradation of host mRNAs. In this review, we discuss the diverse mechanisms viruses use to degrade host mRNAs. In addition, the widespread degradation of host mRNAs can have common consequences including the accumulation of RNA binding proteins in the nucleus, which leads to altered RNA processing, mRNA export, and changes to transcription.

Transcriptional reprogramming at the intersection of the heat shock response and proteostasis

Cellular homeostasis is constantly challenged by a myriad of extrinsic and intrinsic stressors. To mitigate the stress-induced damage, cells activate transient survival programs. The heat shock response (HSR) is an evolutionarily well-conserved survival program that is activated in response to proteotoxic stress. The HSR encompasses a dual regulation of transcription, characterized by rapid activation of genes encoding molecular chaperones and concomitant global attenuation of non-chaperone genes. Recent genome-wide approaches have delineated the molecular depth of stress-induced transcriptional reprogramming. The dramatic rewiring of gene and enhancer networks is driven by key transcription factors, including heat shock factors (HSFs), that together with chromatin-modifying enzymes remodel the 3D chromatin architecture, determining the selection of either gene activation or repression. Here, we highlight the current advancements of molecular mechanisms driving transcriptional reprogramming during acute heat stress. We also discuss the emerging implications of HSF-mediated stress signaling in the context of physiological and pathological conditions.

The HSV-1 ICP22 protein selectively impairs histone repositioning upon Pol II transcription downstream of genes

Herpes simplex virus 1 (HSV-1) infection and stress responses disrupt transcription termination by RNA Polymerase II (Pol II). In HSV-1 infection, but not upon salt or heat stress, this is accompanied by a dramatic increase in chromatin accessibility downstream of genes. Here, we show that the HSV-1 immediate-early protein ICP22 is both necessary and sufficient to induce downstream open chromatin regions (dOCRs) when transcription termination is disrupted by the viral ICP27 protein. This is accompanied by a marked ICP22-dependent loss of histones downstream of affected genes consistent with impaired histone repositioning in the wake of Pol II. Efficient knock-down of the ICP22-interacting histone chaperone FACT is not sufficient to induce dOCRs in ΔICP22 infection but increases dOCR induction in wild-type HSV-1 infection. Interestingly, this is accompanied by a marked increase in chromatin accessibility within gene bodies. We propose a model in which allosteric changes in Pol II composition downstream of genes and ICP22-mediated interference with FACT activity explain the differential impairment of histone repositioning downstream of genes in the wake of Pol II in HSV-1 infection.

Intergenic transcription in in vivo developed bovine oocytes and pre-implantation embryos

Background

Intergenic transcription, either failure to terminate at the transcription end site (TES), or transcription initiation at other intergenic regions, is present in cultured cells and enhanced in the presence of stressors such as viral infection. Transcription termination failure has not been characterized in natural biological samples such as pre-implantation embryos which express more than 10,000 genes and undergo drastic changes in DNA methylation.

Results

Using Automatic Readthrough Transcription Detection (ARTDeco) and data of in vivo developed bovine oocytes and embryos, we found abundant intergenic transcripts that we termed as read-outs (transcribed from 5 to 15 kb after TES) and read-ins (transcribed 1 kb up-stream of reference genes, extending up to 15 kb up-stream). Read-throughs (continued transcription from TES of expressed reference genes, 4-15 kb in length), however, were much fewer. For example, the numbers of read-outs and read-ins ranged from 3,084 to 6,565 or 33.36-66.67% of expressed reference genes at different stages of embryo development. The less copious read-throughs were at an average of 10% and significantly correlated with reference gene expression (P < 0.05). Interestingly, intergenic transcription did not seem to be random because many intergenic transcripts (1,504 read-outs, 1,045 read-ins, and 1,021 read-throughs) were associated with common reference genes across all stages of pre-implantation development. Their expression also seemed to be regulated by developmental stages because many were differentially expressed (log2 fold change ≥ 2, P < 0.05). Additionally, while gradual but un-patterned decreases in DNA methylation densities 10 kb both up- and down-stream of the intergenic transcribed regions were observed, the correlation between intergenic transcription and DNA methylation was insignificant. Finally, transcription factor binding motifs and polyadenylation signals were found in 27.2% and 12.15% of intergenic transcripts, respectively, suggesting considerable novel transcription initiation and RNA processing.

Conclusion

In summary, in vivo developed oocytes and pre-implantation embryos express large numbers of intergenic transcripts, which are not related to the overall DNA methylation profiles either up- or down-stream.

Promoter-pervasive transcription causes RNA polymerase II pausing to boost DOG1 expression in response to salt

Eukaryotic genomes are pervasively transcribed by RNA polymerase II. Yet, the molecular and biological implications of such a phenomenon are still largely puzzling. Here, we describe noncoding RNA transcription upstream of the Arabidopsis thaliana DOG1 gene, which governs salt stress responses and is a key regulator of seed dormancy. We find that expression of the DOG1 gene is induced by salt stress, thereby causing a delay in seed germination. We uncover extensive transcriptional activity on the promoter of the DOG1 gene, which produces a variety of lncRNAs. These lncRNAs, named PUPPIES, are co-directionally transcribed and extend into the DOG1 coding region. We show that PUPPIES RNAs respond to salt stress and boost DOG1 expression, resulting in delayed germination. This positive role of pervasive PUPPIES transcription on DOG1 gene expression is associated with augmented pausing of RNA polymerase II, slower transcription and higher transcriptional burst size. These findings highlight the positive role of upstream co-directional transcription in controlling transcriptional dynamics of downstream genes.

Activation of HIV-1 proviruses increases downstream chromatin accessibility

It is unclear how the activation of HIV-1 transcription affects chromatin structure. We interrogated chromatin organization both genome-wide and nearby HIV-1 integration sites using Hi-C and ATAC-seq. In conjunction, we analyzed the transcription of the HIV-1 genome and neighboring genes. We found that long-range chromatin contacts did not differ significantly between uninfected cells and those harboring an integrated HIV-1 genome, whether the HIV-1 genome was actively transcribed or inactive. Instead, the activation of HIV-1 transcription changes chromatin accessibility immediately downstream of the provirus, demonstrating that HIV-1 can alter local cellular chromatin structure. Finally, we examined HIV-1 and neighboring host gene transcripts with long-read sequencing and found populations of chimeric RNAs both virus-to-host and host-to-virus. Thus, multiomics profiling revealed that the activation of HIV-1 transcription led to local changes in chromatin organization and altered the expression of neighboring host genes.

RNase L activation in the cytoplasm induces aberrant processing of mRNAs in the nucleus

The antiviral endoribonuclease, RNase L, is activated by the mammalian innate immune response to destroy host and viral RNA to ultimately reduce viral gene expression. Herein, we show that RNase L and RNase L-mediated mRNA decay are primarily localized to the cytoplasm. Consequently, RNA-binding proteins (RBPs) translocate from the cytoplasm to the nucleus upon RNase L activation due to the presence of intact nuclear RNA. The re-localization of RBPs to the nucleus coincides with global alterations to RNA processing in the nucleus. While affecting many host mRNAs, these alterations are pronounced in mRNAs encoding type I and type III interferons and correlate with their retention in the nucleus and reduction in interferon protein production. Similar RNA processing defects also occur during infection with either dengue virus or SARS-CoV-2 when RNase L is activated. These findings reveal that the distribution of RBPs between the nucleus and cytosol is dictated by the availability of RNA in each compartment. Thus, viral infections that trigger RNase L-mediated cytoplasmic RNA in the cytoplasm also alter RNA processing in the nucleus, resulting in an ingenious multi-step immune block to protein biogenesis.

It's a DoG-eat-DoG world-altered transcriptional mechanisms drive downstream-of-gene (DoG) transcript production

The past decade has revolutionized our understanding of regulatory noncoding RNAs (ncRNAs). Among the most recently identified ncRNAs are downstream-of-gene (DoG)-containing transcripts that are produced by widespread transcriptional readthrough. The discovery of DoGs has set the stage for future studies to address many unanswered questions regarding the mechanisms that promote readthrough transcription, RNA processing, and the cellular functions of the unique transcripts. In this review, we summarize current findings regarding the biogenesis, function, and mechanisms regulating this exciting new class of RNA molecules.

Who let the DoGs out? - biogenesis of stress-induced readthrough transcripts

Readthrough transcription caused by inefficient 3'-end cleavage of nascent mRNAs has emerged as a hallmark of the mammalian cellular stress response and results in the production of long noncoding RNAs known as downstream-of-gene (DoG)-containing transcripts. DoGs arise from around 10% of human protein-coding genes and are retained in the nucleus. They are produced minutes after cell exposure to stress and can be detected hours after stress removal. However, their biogenesis and the role(s) that DoGs or their production play in the cellular stress response are incompletely understood. We discuss findings that implicate host and viral proteins in the mechanisms underlying DoG production, as well as the transcriptional landscapes that accompany DoG induction under different stress conditions.

Integrator enforces the fidelity of transcriptional termination at protein-coding genes

Integrator regulates the 3′-end processing and termination of multiple classes of noncoding RNAs. Depletion of INTS11, the catalytic subunit of Integrator, or ectopic expression of its catalytic dead enzyme impairs the 3′-end processing and termination of a set of protein-coding transcripts termed Integrator-regulated termination (IRT) genes. This defect is manifested by increased RNA polymerase II (RNAPII) readthrough and occupancy of serine-2 phosphorylated RNAPII, de novo trimethylation of lysine-36 on histone H3, and a compensatory elevation of the cleavage and polyadenylation (CPA) complex beyond the canonical polyadenylation sites. 3′ RNA sequencing reveals that proximal polyadenylation site usage relies on the endonuclease activity of INTS11. The DNA sequence encompassing the transcription end sites of IRT genes features downstream polyadenylation motifs and an enrichment of GC content that permits the formation of secondary structures within the 3′UTR. Together, this study identifies a subset of protein-coding transcripts whose 3′ end processing requires the Integrator complex.

PRMT1-dependent regulation of RNA metabolism and DNA damage response sustains pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive cancer that has remained clinically challenging to manage. Here we employ an RNAi-based in vivo functional genomics platform to determine epigenetic vulnerabilities across a panel of patient-derived PDAC models. Through this, we identify protein arginine methyltransferase 1 (PRMT1) as a critical dependency required for PDAC maintenance. Genetic and pharmacological studies validate the role of PRMT1 in maintaining PDAC growth. Mechanistically, using proteomic and transcriptomic analyses, we demonstrate that global inhibition of asymmetric arginine methylation impairs RNA metabolism, which includes RNA splicing, alternative polyadenylation, and transcription termination. This triggers a robust downregulation of multiple pathways involved in the DNA damage response, thereby promoting genomic instability and inhibiting tumor growth. Taken together, our data support PRMT1 as a compelling target in PDAC and informs a mechanism-based translational strategy for future therapeutic development.Statement of significancePDAC is a highly lethal cancer with limited therapeutic options. This study identified and characterized PRMT1-dependent regulation of RNA metabolism and coordination of key cellular processes required for PDAC tumor growth, defining a mechanism-based translational hypothesis for PRMT1 inhibitors.

Tracking pre-mRNA maturation across subcellular compartments identifies developmental gene regulation through intron retention and nuclear anchoring

Steps of mRNA maturation are important gene regulatory events that occur in distinct cellular locations. However, transcriptomic analyses often lose information on the subcellular distribution of processed and unprocessed transcripts. We generated extensive RNA-seq data sets to track mRNA maturation across subcellular locations in mouse embryonic stem cells, neuronal progenitor cells, and postmitotic neurons. We find disparate patterns of RNA enrichment between the cytoplasmic, nucleoplasmic, and chromatin fractions, with some genes maintaining more polyadenylated RNA in chromatin than in the cytoplasm. We bioinformatically defined four regulatory groups for intron retention, including complete cotranscriptional splicing, complete intron retention in the cytoplasmic RNA, and two intron groups present in nuclear and chromatin transcripts but fully excised in cytoplasm. We found that introns switch their regulatory group between cell types, including neuronally excised introns repressed by polypyrimidine track binding protein 1 (PTBP1). Transcripts for the neuronal gamma-aminobutyric acid (GABA) B receptor, 1 (Gabbr1) are highly expressed in mESCs but are absent from the cytoplasm. Instead, incompletely spliced Gabbr1 RNA remains sequestered on chromatin, where it is bound by PTBP1, similar to certain long noncoding RNAs. Upon neuronal differentiation, Gabbr1 RNA becomes fully processed and exported for translation. Thus, splicing repression and chromatin anchoring of RNA combine to allow posttranscriptional regulation of Gabbr1 over development. For this and other genes, polyadenylated RNA abundance does not indicate functional gene expression. Our data sets provide a rich resource for analyzing many other aspects of mRNA maturation in subcellular locations and across development.

Hyperosmotic stress alters the RNA polymerase II interactome and induces readthrough transcription despite widespread transcriptional repression

Stress-induced readthrough transcription results in the synthesis of downstream-of-gene (DoG)-containing transcripts. The mechanisms underlying DoG formation during cellular stress remain unknown. Nascent transcription profiles during DoG induction in human cell lines using TT-TimeLapse sequencing revealed widespread transcriptional repression upon hyperosmotic stress. Yet, DoGs are produced regardless of the transcriptional level of their upstream genes. ChIP sequencing confirmed that stress-induced redistribution of RNA polymerase (Pol) II correlates with the transcriptional output of genes. Stress-induced alterations in the Pol II interactome are observed by mass spectrometry. While certain cleavage and polyadenylation factors remain Pol II associated, Integrator complex subunits dissociate from Pol II under stress leading to a genome-wide loss of Integrator on DNA. Depleting the catalytic subunit of Integrator using siRNAs induces hundreds of readthrough transcripts, whose parental genes partially overlap those of stress-induced DoGs. Our results provide insights into the mechanisms underlying DoG production and how Integrator activity influences DoG transcription.

Comprehensive Characterization of Transcriptional Activity during Influenza A Virus Infection Reveals Biases in Cap-Snatching of Host RNA Sequences

Infection with influenza A virus (IAV) infection is responsible for an estimated 500,000 deaths and up to 5 million cases of severe respiratory illness each year. In this study, we looked at human primary immune cells (macrophages) infected with IAV. Our method allows us to look at both the host and the virus in parallel. We used these data to explore a process known as “cap-snatching,” where IAV snatches a short nucleotide sequence from capped host RNA. This process was believed to be random. We demonstrate biased snatching of numerous host RNAs, including those associated with snRNA transcription, and avoidance of host transcripts encoding host ribosomal proteins, which are required by IAV for replication. We then describe the transcriptional landscape of the host response to IAV, observing new features, including a failure of IAV-treated MDMs to induce feedback inhibitors of inflammation, seen in response to other treatments.

Macrophages in the lung detect and respond to influenza A virus (IAV), determining the nature of the immune response. Using terminal-depth cap analysis of gene expression (CAGE), we quantified transcriptional activity of both host and pathogen over a 24-h time course of IAV infection in primary human monocyte-derived macrophages (MDMs). This method allowed us to observe heterogenous host sequences incorporated into IAV mRNA, “snatched” 5′ RNA caps, and corresponding RNA sequences from host RNAs. In order to determine whether cap-snatching is random or exhibits a bias, we systematically compared host sequences incorporated into viral mRNA (“snatched”) against a complete survey of all background host RNA in the same cells, at the same time. Using a computational strategy designed to eliminate sources of bias due to read length, sequencing depth, and multimapping, we were able to quantify overrepresentation of host RNA features among the sequences that were snatched by IAV. We demonstrate biased snatching of numerous host RNAs, particularly small nuclear RNAs (snRNAs), and avoidance of host transcripts encoding host ribosomal proteins, which are required by IAV for replication. We then used a systems approach to describe the transcriptional landscape of the host response to IAV, observing many new features, including a failure of IAV-treated MDMs to induce feedback inhibitors of inflammation, seen in response to other treatments.

IMPORTANCE Infection with influenza A virus (IAV) infection is responsible for an estimated 500,000 deaths and up to 5 million cases of severe respiratory illness each year. In this study, we looked at human primary immune cells (macrophages) infected with IAV. Our method allows us to look at both the host and the virus in parallel. We used these data to explore a process known as “cap-snatching,” where IAV snatches a short nucleotide sequence from capped host RNA. This process was believed to be random. We demonstrate biased snatching of numerous host RNAs, including those associated with snRNA transcription, and avoidance of host transcripts encoding host ribosomal proteins, which are required by IAV for replication. We then describe the transcriptional landscape of the host response to IAV, observing new features, including a failure of IAV-treated MDMs to induce feedback inhibitors of inflammation, seen in response to other treatments.

Gallop Racing Shifts Mature mRNA towards Introns: Does Exercise-Induced Stress Enhance Genome Plasticity?

Physical exercise is universally recognized as stressful. Among the "sport species", the horse is probably the most appropriate model for investigating the genomic response to stress due to the homogeneity of its genetic background. The aim of this work is to dissect the whole transcription modulation in Peripheral Blood Mononuclear Cells (PBMCs) after exercise with a time course framework focusing on unexplored regions related to introns and intergenic portions. PBMCs NGS from five 3 year old Sardinian Anglo-Arab racehorses collected at rest and after a 2000 m race was performed. Apart from differential gene expression ascertainment between the two time points the complexity of transcription for alternative transcripts was identified. Interestingly, we noted a transcription shift from the coding to the non-coding regions. We further investigated the possible causes of this phenomenon focusing on genomic repeats, using a differential expression approach and finding a strong general up-regulation of repetitive elements such as LINE. Since their modulation is also associated with the "exonization", the recruitment of repeats that act with regulatory functions, suggesting that there might be an active regulation of this transcriptional shift. Thanks to an innovative bioinformatic approach, our study could represent a model for the transcriptomic investigation of stress.

Using multiple reference genomes to identify and resolve annotation inconsistencies

BACKGROUND

Advances in sequencing technologies have led to the release of reference genomes and annotations for multiple individuals within more well-studied systems. While each of these new genome assemblies shares significant portions of synteny between each other, the annotated structure of gene models within these regions can differ. Of particular concern are split-gene misannotations, in which a single gene is incorrectly annotated as two distinct genes or two genes are incorrectly annotated as a single gene. These misannotations can have major impacts on functional prediction, estimates of expression, and many downstream analyses.

RESULTS

We developed a high-throughput method based on pairwise comparisons of annotations that detect potential split-gene misannotations and quantifies support for whether the genes should be merged into a single gene model. We demonstrated the utility of our method using gene annotations of three reference genomes from maize (B73, PH207, and W22), a difficult system from an annotation perspective due to the size and complexity of the genome. On average, we found several hundred of these potential split-gene misannotations in each pairwise comparison, corresponding to 3-5% of gene models across annotations. To determine which state (i.e. one gene or multiple genes) is biologically supported, we utilized RNAseq data from 10 tissues throughout development along with a novel metric and simulation framework. The methods we have developed require minimal human interaction and can be applied to future assemblies to aid in annotation efforts.

CONCLUSIONS

Split-gene misannotations occur at appreciable frequency in maize annotations. We have developed a method to easily identify and correct these misannotations. Importantly, this method is generic in that it can utilize any type of short-read expression data. Failure to account for split-gene misannotations has serious consequences for biological inference, particularly for expression-based analyses.

Using multiple reference genomes to identify and resolve annotation inconsistencies

BACKGROUND

Advances in sequencing technologies have led to the release of reference genomes and annotations for multiple individuals within more well-studied systems. While each of these new genome assemblies shares significant portions of synteny between each other, the annotated structure of gene models within these regions can differ. Of particular concern are split-gene misannotations, in which a single gene is incorrectly annotated as two distinct genes or two genes are incorrectly annotated as a single gene. These misannotations can have major impacts on functional prediction, estimates of expression, and many downstream analyses.

RESULTS

We developed a high-throughput method based on pairwise comparisons of annotations that detect potential split-gene misannotations and quantifies support for whether the genes should be merged into a single gene model. We demonstrated the utility of our method using gene annotations of three reference genomes from maize (B73, PH207, and W22), a difficult system from an annotation perspective due to the size and complexity of the genome. On average, we found several hundred of these potential split-gene misannotations in each pairwise comparison, corresponding to 3-5% of gene models across annotations. To determine which state (i.e. one gene or multiple genes) is biologically supported, we utilized RNAseq data from 10 tissues throughout development along with a novel metric and simulation framework. The methods we have developed require minimal human interaction and can be applied to future assemblies to aid in annotation efforts.

CONCLUSIONS

Split-gene misannotations occur at appreciable frequency in maize annotations. We have developed a method to easily identify and correct these misannotations. Importantly, this method is generic in that it can utilize any type of short-read expression data. Failure to account for split-gene misannotations has serious consequences for biological inference, particularly for expression-based analyses.

The landscape of eukaryotic mRNPs

The proper regulation of mRNA processing, localization, translation, and degradation occurs on mRNPs. However, the global principles of mRNP organization are poorly understood. We utilize the limited, but existing, information available to present a speculative synthesis of mRNP organization with the following key points. First, mRNPs form a compacted structure due to the inherent folding of RNA. Second, the ribosome is the principal mechanism by which mRNA regions are partially decompacted. Third, mRNPs are 50%-80% protein by weight, consistent with proteins modulating mRNP organization, but also suggesting the majority of mRNA sequences are not directly interacting with RNA-binding proteins. Finally, the ratio of mRNA-binding proteins to mRNAs is higher in the nucleus to allow effective RNA processing and limit the potential for nuclear RNA based aggregation. This synthesis of mRNP understanding provides a model for mRNP biogenesis, structure, and regulation with multiple implications.

Influenza Virus Mounts a Two-Pronged Attack on Host RNA Polymerase II Transcription

Influenza virus intimately associates with host RNA polymerase II (Pol II) and mRNA processing machinery. Here, we use mammalian native elongating transcript sequencing (mNET-seq) to examine Pol II behavior during viral infection. We show that influenza virus executes a two-pronged attack on host transcription. First, viral infection causes decreased Pol II gene occupancy downstream of transcription start sites. Second, virus-induced cellular stress leads to a catastrophic failure of Pol II termination at poly(A) sites, with transcription often continuing for tens of kilobases. Defective Pol II termination occurs independently of the ability of the viral NS1 protein to interfere with host mRNA processing. Instead, this termination defect is a common effect of diverse cellular stresses and underlies the production of previously reported downstream-of-gene transcripts (DoGs). Our work has implications for understanding not only host-virus interactions but also fundamental aspects of mammalian transcription.

•

Influenza virus infection dysregulates host transcription

•

Viral infection depletes Pol II from gene bodies downstream of the TSS

•

Virus-induced stress leads to a catastrophic failure of Pol II termination

•

Defective termination does not require viral NS1: host CPSF30 interaction

Increased grp78 transcription is correlated to reduced tlr4 transcription in patients surviving sepsis

Regulated transcriptional readthrough during stress maintains genome structure and ensures access to genes that are necessary for cellular recovery. A broad number of genes, including of the bacterial sensor Toll-like receptor 4 (TLR-4), are markedly transcribed on initiating the systemic inflammatory response. Here we study the transcriptional patterns of tlr4 and of its modulator grp78 during human sepsis, and establish their correlations with the outcome of patients. We measured the daily tlr4 and grp78 RNA expression levels in peripheral blood of septic patients, immediately after admission to intensive care, and modeled these RNA values with a sine damping function. We obtained negative correlations between the transcription of tlr4 and grp78 RNA in the survivor group. In contrast, such relation is lost in the deceased patients. Loss of transcriptional homeostasis predicted by our model within the initial 4 days of hospitalization was confirmed by death of those patients up to 28 days later.

A cell-based splicing reporter system to identify regulators of cis-splicing between adjacent genes

Chimeric RNAs generated by cis-splicing between adjacent genes (cis-SAGe) are increasingly recognized as a widespread phenomenon. These chimeric messenger RNAs are present in normal human cells, and are also detected in various cancers. The mechanisms for how this group of chimeras is formed are not yet clear, in part due to the lack of a tractable system for their experimental investigation. Here we developed a fast, easy and versatile cell-based reporter system to identify regulators of cis-SAGe. The reporter, consisting of four main cassettes, simultaneously measures the effects of a candidate regulator on cis-SAGe and canonical splicing. Using this cell-based assay, we screened 102 candidate factors involved in RNA pol II cleavage and termination, elongation, splicing, alternative splicing and R-loop formation. We discovered that two factors, SRRM1 and SF3B1, affect not only cis-SAGe chimeras, but also other types of chimeric RNAs in a genome-wide fashion. This system can be used for studying trans-acting factors and cis-acting sequence elements and factors, as well as for screening small molecule inhibitors.

DoGFinder: a software for the discovery and quantification of readthrough transcripts from RNA-seq

Background

Recent studies have described a widespread induction of transcriptional readthrough as a consequence of various stress conditions in mammalian cells. This novel phenomenon, initially identified from analysis of RNA-seq data, suggests intriguing new levels of gene expression regulation. However, the mechanism underlying naturally occurring transcriptional readthrough, as well as its regulatory consequences, still remain elusive. Furthermore, the readthrough response to stress has thus far not been investigated outside of mammalian species, and the occurrence of readthrough in many physiological and disease conditions remains to be explored.

Results

To facilitate a wider investigation into transcriptional readthrough, we created the DoGFinder software package, for the streamlined identification and quantification of readthrough transcripts, also known as DoGs (Downstream of Gene-containing transcripts), from any RNA-seq dataset. Using DoGFinder, we explore the dependence of DoG discovery potential on RNA-seq library depth, and show that stress-induced readthrough induction discovery is robust to sequencing depth, and input parameter settings. We further demonstrate the use of the DoGFinder software package on a new publically available RNA-seq dataset, and discover DoG induction in human PME cells following hypoxia – a previously unknown readthrough inducing stress type.

Conclusions

DoGFinder will enable users to explore, in a few simple steps, the readthrough phenomenon in any condition and organism.

DoGFinder is freely available at https://github.com/shalgilab/DoGFinder.

Electronic supplementary material

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

Transcriptional-Readthrough RNAs Reflect the Phenomenon of "A Gene Contains Gene(s)" or "Gene(s) within a Gene" in the Human Genome, and Thus Are Not Chimeric RNAs

Tens of thousands of chimeric RNAs, i.e., RNAs with sequences of two genes, have been identified in human cells. Most of them are formed by two neighboring genes on the same chromosome and are considered to be derived via transcriptional readthrough, but a true readthrough event still awaits more evidence and trans-splicing that joins two transcripts together remains as a possible mechanism. We regard those genomic loci that are transcriptionally read through as unannotated genes, because their transcriptional and posttranscriptional regulations are the same as those of already-annotated genes, including fusion genes formed due to genetic alterations. Therefore, readthrough RNAs and fusion-gene-derived RNAs are not chimeras. Only those two-gene RNAs formed at the RNA level, likely via trans-splicing, without corresponding genes as genomic parents, should be regarded as authentic chimeric RNAs. However, since in human cells, procedural and mechanistic details of trans-splicing have never been disclosed, we doubt the existence of trans-splicing. Therefore, there are probably no authentic chimeras in humans, after readthrough and fusion-gene derived RNAs are all put back into the group of ordinary RNAs. Therefore, it should be further determined whether in human cells all two-neighboring-gene RNAs are derived from transcriptional readthrough and whether trans-splicing truly exists.

Connections between Transcription Downstream of Genes and cis-SAGe Chimeric RNA

cis-Splicing between adjacent genes (cis-SAGe) is being recognized as one way to produce chimeric fusion RNAs. However, its detail mechanism is not clear. Recent study revealed induction of transcriptions downstream of genes (DoGs) under osmotic stress. Here, we investigated the influence of osmotic stress on cis-SAGe chimeric RNAs and their connection to DoGs. We found, the absence of induction of at least some cis-SAGe fusions and/or their corresponding DoGs at early time point(s). In fact, these DoGs and their cis-SAGe fusions are inversely correlated. This negative correlation was changed to positive at a later time point. These results suggest a direct competition between the two categories of transcripts when total pool of readthrough transcripts is limited at an early time point. At a later time point, DoGs and corresponding cis-SAGe fusions are both induced, indicating that total readthrough transcripts become more abundant. Finally, we observed overall enhancement of cis-SAGe chimeric RNAs in KCl-treated samples by RNA-Seq analysis.

Comparative analysis reveals genomic features of stress-induced transcriptional readthrough

Transcription is a highly regulated process, and stress-induced changes in gene transcription have been shown to play a major role in stress responses and adaptation. Genome-wide studies reveal prevalent transcription beyond known protein-coding gene loci, generating a variety of RNA classes, most of unknown function. One such class, termed downstream of gene-containing transcripts (DoGs), was reported to result from transcriptional readthrough upon osmotic stress in human cells. However, how widespread the readthrough phenomenon is, and what its causes and consequences are, remain elusive. Here we present a genome-wide mapping of transcriptional readthrough, using nuclear RNA-Seq, comparing heat shock, osmotic stress, and oxidative stress in NIH 3T3 mouse fibroblast cells. We observe massive induction of transcriptional readthrough, both in levels and length, under all stress conditions, with significant, yet not complete, overlap of readthrough-induced loci between different conditions. Importantly, our analyses suggest that stress-induced transcriptional readthrough is not a random failure process, but is rather differentially induced across different conditions. We explore potential regulators and find a role for HSF1 in the induction of a subset of heat shock-induced readthrough transcripts. Analysis of public datasets detected increases in polymerase II occupancy in DoG regions after heat shock, supporting our findings. Interestingly, DoGs tend to be produced in the vicinity of neighboring genes, leading to a marked increase in their antisense-generating potential. Finally, we examine genomic features of readthrough transcription and observe a unique chromatin signature typical of DoG-producing regions, suggesting that readthrough transcription is associated with the maintenance of an open chromatin state.

Mot1, Ino80C, and NC2 Function Coordinately to Regulate Pervasive Transcription in Yeast and Mammals

Pervasive transcription initiates from cryptic promoters and is observed in eukaryotes ranging from yeast to mammals. The Set2-Rpd3 regulatory system prevents cryptic promoter function within expressed genes. However, conserved systems that control pervasive transcription within intergenic regions have not been well established. Here we show that Mot1, Ino80 chromatin remodeling complex (Ino80C), and NC2 co-localize on chromatin and coordinately suppress pervasive transcription in S. cerevisiae and murine embryonic stem cells (mESCs). In yeast, all three proteins bind subtelomeric heterochromatin through a Sir3-stimulated mechanism and to euchromatin via a TBP-stimulated mechanism. In mESCs, the proteins bind to active and poised TBP-bound promoters along with promoters of polycomb-silenced genes apparently lacking TBP. Depletion of Mot1, Ino80C, or NC2 by anchor away in yeast or RNAi in mESCs leads to near-identical transcriptome phenotypes, with new subtelomeric transcription in yeast, and greatly increased pervasive transcription in both yeast and mESCs.

Alternative Start Sites Downstream of Non-Sense Mutations Drive Antigen Presentation and Tolerance Induction to C-Terminal Epitopes

CTL responses to the transgene product remain an active area of concern for the gene therapy field. A patient's underlying genetic mutation may influence the qualitative nature of these potentially destructive T cell responses. Individuals with a mutation that introduces a premature termination codon (PTC) that prevents synthesis of the full-length peptide are considered more likely to mount a transgene-specific T cell response because of a lack of immune tolerance to C-terminal epitopes as a consequence of absent endogenous Ag presentation. In this article, we demonstrate that a human ornithine transcarbamylase gene containing various PTC-inducing non-sense mutations is able to generate and present epitopes downstream of the termination codon. Generation of these epitopes occurs primarily from alternative translation start sites downstream of the stop codon. Furthermore, we show that expression of these genes from adeno-associated virus vectors in C57BL/6 mice is able to induce peripheral tolerance to epitopes downstream of the PTC. These results suggest that, despite the lack of full-length endogenous protein, patients with PTC-inducing non-sense mutations may still present T cell epitopes downstream of the premature termination site that may render the subject tolerant to wild-type transgene products.

Recent advances in understanding transcription termination by RNA polymerase II

Transcription termination is a fundamental process in which RNA polymerase ceases RNA chain extension and dissociates from the chromatin template, thereby defining the end of the transcription unit. Our understanding of the biological role and functional importance of termination by RNA polymerase II and the range of processes in which it is involved has grown significantly in recent years. A large set of nucleic acid-binding proteins and enzymes have been identified as part of the termination machinery. A greater appreciation for the coupling of termination to RNA processing and metabolism has been recognized. In addition to serving as an essential step at the end of the transcription cycle, termination is involved in the regulation of a broad range of cellular processes. More recently, a role for termination in pervasive transcription, non-coding RNA regulation, genetic stability, chromatin remodeling, the immune response, and disease has come to the fore. Interesting mechanistic questions remain, but the last several years have resulted in significant insights into termination and an increasing recognition of its biological importance.