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

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.

Genome-wide kinetic profiling of pre-mRNA 3' end cleavage

Cleavage and polyadenylation is necessary for the formation of mature mRNA molecules. The rate at which this process occurs can determine the temporal availability of mRNA for subsequent function throughout the cell and is likely tightly regulated. Despite advances in high-throughput approaches for global kinetic profiling of RNA maturation, genome-wide 3' end cleavage rates have never been measured. Here, we describe a novel approach to estimate the rates of cleavage, using metabolic labeling of nascent RNA, high-throughput sequencing, and mathematical modeling. Using in silico simulations of nascent RNA-seq data, we show that our approach can accurately and precisely estimate cleavage half-lives for both constitutive and alternative sites. We find that 3' end cleavage is fast on average, with half-lives under a minute, but highly variable across individual sites. Rapid cleavage is promoted by the presence of canonical sequence elements and an increased density of polyadenylation signals near a cleavage site. Finally, we find that cleavage rates are associated with the localization of RNA polymerase II at the end of a gene, and faster cleavage leads to quicker degradation of downstream readthrough RNA. Our findings shed light on the features important for efficient 3' end cleavage and the regulation of transcription termination.

Transcription readthrough is prevalent in healthy human tissues and associated with inherent genomic features

Transcription termination is a crucial step in the production of conforming mRNAs and functional proteins. Under cellular stress conditions, the transcription machinery fails to identify the termination site and continues transcribing beyond gene boundaries, a phenomenon designated as transcription readthrough. However, the prevalence and impact of this phenomenon in healthy human tissues remain unexplored. Here, we assessed transcription readthrough in almost 3000 transcriptome profiles representing 23 human tissues and found that 34% of the expressed protein-coding genes produced readthrough transcripts. The production of readthrough transcripts was restricted in genomic regions with high transcriptional activity and was associated with inefficient splicing and increased chromatin accessibility in terminal regions. In addition, we showed that these transcripts contained several binding sites for the same miRNA, unravelling a potential role as miRNA sponges. Overall, this work provides evidence that transcription readthrough is pervasive and non-stochastic, not only in abnormal conditions but also in healthy tissues. This suggests a potential role for such transcripts in modulating normal cellular functions.

Transcriptional Stress Induces the Generation of DoGs in Cancer Cells

A characteristic of the cellular response to stress is the production of RNAs generated from a readthrough transcription of genes, called downstream-of-gene-(DoG)-containing transcripts. Additionally, transcription inhibitor drugs are candidates for fighting cancer. In this work, we report the results of a bioinformatic analysis showing that one of the responses to transcription inhibition is the generation of DoGs in cancer cells. Although some genes that form DoGs were shared between the two cancer lines, there did not appear to be a functional correlation between them. However, our findings show that DoGs are generated as part of the cellular response to transcription inhibition like other types of cellular stress, suggesting that they may be part of the defense against transcriptional stress.

Protocol for quantitative analysis of RNA 3'-end processing induced by disassociated subunits using chromatin-associated RNA-seq data

Sequencing chromatin-associated RNA using libraries from the chromatin fraction makes it possible to characterize RNA processing driven by disassociated subunits. Here, we present an experimental strategy and computational pipeline for processing chromatin-associated RNA-seq data to detect and quantify readthrough transcripts. We describe steps for constructing degron mouse embryonic stem cells, detecting readthrough genes, data processing, and data analysis. This protocol can be adapted to various biological scenarios and other types of nascent RNA-seq, such as TT-seq. For complete details on the use and execution of this protocol, please refer to Li et al. (2023).1.

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.

RNA Pol II preferentially regulates ribosomal protein expression by trapping disassociated subunits

RNA polymerase II (RNA Pol II) has been recognized as a passively regulated multi-subunit holoenzyme. However, the extent to which RNA Pol II subunits might be important beyond the RNA Pol II complex remains unclear. Here, fractions containing disassociated RPB3 (dRPB3) were identified by size exclusion chromatography in various cells. Through a unique strategy, i.e., "specific degradation of disassociated subunits (SDDS)," we demonstrated that dRPB3 functions as a regulatory component of RNA Pol II to enable the preferential control of 3' end processing of ribosomal protein genes directly through its N-terminal domain. Machine learning analysis of large-scale genomic features revealed that the little elongation complex (LEC) helps to specialize the functions of dRPB3. Mechanistically, dRPB3 facilitates CBC-PCF11 axis activity to increase the efficiency of 3' end processing. Furthermore, RPB3 is dynamically regulated during development and diseases. These findings suggest that RNA Pol II gains specific regulatory functions by trapping disassociated subunits in mammalian cells.

UBE3D Regulates mRNA 3'-End Processing and Maintains Adipogenic Potential in 3T3-L1 Cells

We have previously described the role of an essential Saccharomyces cerevisiae gene, important for cleavage and polyadenylation 1 (IPA1), in the regulation of gene expression through its interaction with Ysh1, the endonuclease subunit of the mRNA 3'-end processing complex. Through a similar mechanism, the mammalian homolog ubiquitin protein ligase E3D (UBE3D) promotes the migratory and invasive potential of breast cancer cells, but its role in the regulation of gene expression during normal cellular differentiation has not previously been described. In this study, we show that CRISPR/Cas9-mediated knockout of Ube3d in 3T3-L1 cells blocks their ability to differentiate into mature adipocytes. Consistent with previous studies in other cell types, Ube3d knockout leads to decreased levels of CPSF73 and global changes in cellular mRNAs indicative of a loss of 3'-end processing capacity. Ube3d knockout cells also display decreased expression of known preadipogenic markers. Overexpression of either UBE3D or CPSF73 rescues the differentiation defect and partially restores protein levels of these markers. These results support a model in which UBE3D is necessary for the maintenance of the adipocyte-committed state via its regulation of the mRNA 3'-end processing machinery.

Stress-induced transcriptional readthrough into neighboring genes is linked to intron retention

Exposure to certain stresses leads to readthrough transcription. Using polyA-selected RNA-seq in mouse fibroblasts subjected to heat shock, oxidative, or osmotic stress, we found that readthrough transcription can proceed into proximal downstream genes, in a phenomenon previously termed "read-in." We found that read-in genes share distinctive genomic characteristics; they are GC-rich and extremely short , with genomic features conserved in human. Using ribosome profiling, we found that read-in genes show significantly reduced translation. Strikingly, read-in genes demonstrate marked intron retention, mostly in their first introns, which could not be explained solely by their short introns and GC-richness, features often associated with intron retention. Finally, we revealed H3K36me3 enrichment upstream to read-in genes. Moreover, demarcation of exon-intron junctions by H3K36me3 was absent in read-in first introns. Our data portray a relationship between read-in and intron retention, suggesting they may have co-evolved to facilitate reduced translation of read-in genes during stress.

Targeted protein degradation reveals RNA Pol II heterogeneity and functional diversity

RNA polymerase II (RNA Pol II) subunits are thought to be involved in various transcription-associated processes, but it is unclear whether they play different regulatory roles in modulating gene expression. Here, we performed nascent and mature transcript sequencing after the acute degradation of 12 mammalian RNA Pol II subunits and profiled their genomic binding sites and protein interactomes to dissect their molecular functions. We found that RNA Pol II subunits contribute differently to RNA Pol II cellular localization and transcription processes and preferentially regulate RNA processing (such as RNA splicing and 3' end maturation). Genes sensitive to the depletion of different RNA Pol II subunits tend to be involved in diverse biological functions and show different RNA half-lives. Sequences, associated protein factors, and RNA structures are correlated with RNA Pol II subunit-mediated differential gene expression. These findings collectively suggest that the heterogeneity of RNA Pol II and different genes appear to depend on some of the subunits.

Targeting the mRNA endonuclease CPSF73 inhibits breast cancer cell migration, invasion, and self-renewal

Cleavage by the endonuclease CPSF73 and polyadenylation of nascent RNA is an essential step in co-transcriptional mRNA maturation. Recent work has surprisingly identified CPSF73 as a promising drug target for inhibiting the growth of specific cancers, triggering further studies on understanding CPSF73 regulation and functions in cells. Here, we report that a HECT-like E3 ligase, UBE3D, participates in stabilizing CPFS73 protein by preventing its ubiquitin-mediated degradation by the proteasome. Depletion of UBE3D leads to CPSF73 downregulation, a pre-mRNA cleavage defect, and dysregulated gene expression in cells. UBE3D dysfunction or chemical inactivation of CPSF73 inhibited migration and invasion as well as stem cell renewal phenotypes in vitro in triple-negative breast cancer cells. In addition, genetic overexpression of CPSF73 promoted breast cancer stemness and knocking down CPSF73 inhibited stem cell renewal properties. Together, our findings indicate that targeting the pre-mRNA processing nuclease CPSF73 has potential for breast cancer therapy.

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.

Comprehensive Detection of Pseudogenes Transcribed by Readthrough

Transcription termination is a critical stage for the production of legitimate mRNAs, and consequently functional proteins. However, the transcription machinery can ignore the stop signs and continue elongating beyond gene boundaries, invading downstream neighboring genes. Such phenomenon, designated transcription readthrough, can trigger the expression of pseudogenes usually silenced or lacking the proper regulatory signals. Due to the sequence similarity to parental genes, readthrough transcribed pseudogenes can regulate relevant protein-coding genes and impact biological functions. Here, we describe a computational pipeline that employs already existent bioinformatic tools to detect readthrough transcribed pseudogenes from expression profiles. We also unveil that combining strand-specific transcriptome data and epigenetic profiles can enhance and corroborate the results. By applying such approach to renal cancer biopsies, we show that pseudogenes can be readthrough transcribed as part of unspliced transcripts or processed RNA chimeras. Overall, our pipeline allows us to scrutinize transcriptome profiles to detect a diversity of readthrough events leading to expression of pseudogenes.

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.

ARTDeco: automatic readthrough transcription detection

BACKGROUND

Mounting evidence suggests several diseases and biological processes target transcription termination to misregulate gene expression. Disruption of transcription termination leads to readthrough transcription past the 3' end of genes, which can result in novel transcripts, changes in epigenetic states and altered 3D genome structure.

RESULTS

We developed Automatic Readthrough Transcription Detection (ARTDeco), a tool to detect and analyze multiple features of readthrough transcription from RNA-seq and other next-generation sequencing (NGS) assays that profile transcriptional activity. ARTDeco robustly quantifies the global severity of readthrough phenotypes, and reliably identifies individual genes that fail to terminate (readthrough genes), are aberrantly transcribed due to upstream termination failure (read-in genes), and novel transcripts created as a result of readthrough (downstream of gene or DoG transcripts). We used ARTDeco to characterize readthrough transcription observed during influenza A virus (IAV) infection, validating its specificity and sensitivity by comparing its performance in samples infected with a mutant virus that fails to block transcription termination. We verify ARTDeco's ability to detect readthrough as well as identify read-in genes from different experimental assays across multiple experimental systems with known defects in transcriptional termination, and show how these results can be leveraged to improve the interpretation of gene expression and downstream analysis. Applying ARTDeco to a gene expression data set from IAV-infected monocytes from different donors, we find strong evidence that read-in gene-associated expression quantitative trait loci (eQTLs) likely regulate genes upstream of read-in genes. This indicates that taking readthrough transcription into account is important for the interpretation of eQTLs in systems where transcription termination is blocked.

CONCLUSIONS

ARTDeco aids researchers investigating readthrough transcription in a variety of systems and contexts.