Transcription Polyadenylation Research Articles

The MTR4/hnRNPK complex surveils aberrant polyadenylated RNAs with multiple exons

RNA surveillance systems degrade aberrant RNAs that result from defective transcriptional termination, splicing, and polyadenylation. Defective RNAs in the nucleus are recognized by RNA-binding proteins and MTR4, and are degraded by the RNA exosome complex. Here, we detect aberrant RNAs in MTR4-depleted cells using long-read direct RNA sequencing and 3′ sequencing. MTR4 destabilizes intronic polyadenylated transcripts generated by transcriptional read-through over one or more exons, termed 3′ eXtended Transcripts (3XTs). MTR4 also associates with hnRNPK, which recognizes 3XTs with multiple exons. Moreover, the aberrant protein translated from KCTD13 3XT is a target of the hnRNPK-MTR4-RNA exosome pathway and forms aberrant condensates, which we name KCTD13 3eXtended Transcript-derived protein (KeXT) bodies. Our results suggest that RNA surveillance in human cells inhibits the formation of condensates of a defective polyadenylated transcript-derived protein.

Uncovering the dynamics and consequences of RNA isoform changes during neuronal differentiation

Static gene expression programs have been extensively characterized in stem cells and mature human cells. However, the dynamics of RNA isoform changes upon cell-state-transitions during cell differentiation, the determinants and functional consequences have largely remained unclear. Here, we established an improved model for human neurogenesis in vitro that is amenable for systems-wide analyses of gene expression. Our multi-omics analysis reveals that the pronounced alterations in cell morphology correlate strongly with widespread changes in RNA isoform expression. Our approach identifies thousands of new RNA isoforms that are expressed at distinct differentiation stages. RNA isoforms mainly arise from exon skipping and the alternative usage of transcription start and polyadenylation sites during human neurogenesis. The transcript isoform changes can remodel the identity and functions of protein isoforms. Finally, our study identifies a set of RNA binding proteins as a potential determinant of differentiation stage-specific global isoform changes. This work supports the view of regulated isoform changes that underlie state-transitions during neurogenesis.

MAPP unravels frequent co-regulation of splicing and polyadenylation by RNA-binding proteins and their dysregulation in cancer

Maturation of eukaryotic pre-mRNAs via splicing and polyadenylation is modulated across cell types and conditions by a variety of RNA-binding proteins (RBPs). Although there exist over 1,500 RBPs in human cells, their binding motifs and functions still remain to be elucidated, especially in the complex environment of tissues and in the context of diseases. To overcome the lack of methods for the systematic and automated detection of sequence motif-guided pre-mRNA processing regulation from RNA sequencing (RNA-Seq) data we have developed MAPP (Motif Activity on Pre-mRNA Processing). Applying MAPP to RBP knock-down experiments reveals that many RBPs regulate both splicing and polyadenylation of nascent transcripts by acting on similar sequence motifs. MAPP not only infers these sequence motifs, but also unravels the position-dependent impact of the RBPs on pre-mRNA processing. Interestingly, all investigated RBPs that act on both splicing and 3’ end processing exhibit a consistently repressive or activating effect on both processes, providing a first glimpse on the underlying mechanism. Applying MAPP to normal and malignant brain tissue samples unveils that the motifs bound by the PTBP1 and RBFOX RBPs coordinately drive the oncogenic splicing program active in glioblastomas demonstrating that MAPP paves the way for characterizing pre-mRNA processing regulators under physiological and pathological conditions.

Therapeutic targeting of CPSF3-dependent transcriptional termination in ovarian cancer.

Transcriptional dysregulation is a recurring pathogenic hallmark and an emerging therapeutic vulnerability in ovarian cancer. Here, we demonstrated that ovarian cancer exhibited a unique dependency on the regulatory machinery of transcriptional termination, particularly, cleavage and polyadenylation specificity factor (CPSF) complex. Genetic abrogation of multiple CPSF subunits substantially hampered neoplastic cell viability, and we presented evidence that their indispensable roles converged on the endonuclease CPSF3. Mechanistically, CPSF perturbation resulted in lengthened 3'-untranslated regions, diminished intronic polyadenylation and widespread transcriptional readthrough, and consequently suppressed oncogenic pathways. Furthermore, we reported the development of specific CPSF3 inhibitors building upon the benzoxaborole scaffold, which exerted potent antitumor activity. Notably, CPSF3 blockade effectively exacerbated genomic instability by down-regulating DNA damage repair genes and thus acted in synergy with poly(adenosine 5'-diphosphate-ribose) polymerase inhibition. These findings establish CPSF3-dependent transcriptional termination as an exploitable driving mechanism of ovarian cancer and provide a promising class of boron-containing compounds for targeting transcription-addicted human malignancies.

Depletion of tRNA CCA-adding enzyme in Mycobacterium tuberculosis leads to polyadenylation of transcripts and precursor tRNAs

In reference to gene annotation, more than half of the tRNA species synthesized by Mycobacterium tuberculosis require the enzymatic addition of the cytosine-cytosine-adenine (CCA) tail, which is indispensable for amino acid charging and tRNA functionality. It makes the mycobacterial CCA-adding enzyme essential for survival of the bacterium and a potential target for novel pipelines in drug discovery avenues. Here, we described the rv3907c gene product, originally annotated as poly(A)polymerase (rv3907c, PcnA) as a functional CCA-adding enzyme (CCAMtb) essential for viability of M. tuberculosis. The depletion of the enzyme affected tRNAs maturation, inhibited bacilli growth, and resulted in abundant accumulation of polyadenylated RNAs. We determined the enzymatic activities displayed by the mycobacterial CCAMtb in vitro and studied the effects of inhibiting of its transcription in bacterial cells. We are the first to properly confirm the existence of RNA polyadenylation in mycobacteria, a previously controversial phenomenon, which we found promoted upon CCA-adding enzyme downexpression.

An RNA-Ligation-Based RACE-PAT Assay to Monitor Poly(A) Tail Length of mRNAs of Interest.

In eukaryotes, a non-templated poly-adenosine (poly(A)) tail is added co-transcriptionally to almost every messenger RNA (mRNA). The length of this poly(A) tail changes during the lifetime of mRNAs and has been shown in many circumstances to be an important factor controlling transcript fates. Yet, the measure of the length of this homogenous nucleotide sequence is technically challenging, making it difficult to assess its dynamic variation. In this chapter, we describe an RNA-ligation-based RACE-PAT (Rapid Amplification of cDNA End-Poly(A) Tail) assay to monitor the poly(A) tail length of mRNAs. In the first step, an RNA oligonucleotide is ligated to mRNA 3' ends providing an anchoring site to prime cDNA synthesis, avoiding the bias introduced by oligo(dT)-derived primers. Afterward, reverse transcription is performed with an anchor primer with a unique 5' extension. The choice of the oligonucleotide 3' end at this step allows further flexibility to amplify modified tails, for example, by uridylation. Next, short DNA fragments encompassing the poly(A) tails are amplified by Polymerase Chain Reaction (PCR) using as forward primer, a transcript-specific primer hybridizing close to the transcript polyadenylation signal, and as reverse primer, an oligonucleotide corresponding to the 5' extension of the primer used for cDNA synthesis, ensuring that only cDNAs are amplified. The resulting DNA fragments are then visualized after size fractionation by electrophoresis. This method does not provide exact nucleotide count and composition but has the advantage of allowing the processing of many samples in parallel at a low cost.

SINEs as Potential Expression Cassettes: Impact of Deletions and Insertions on Polyadenylation and Lifetime of B2 and Ves SINE Transcripts Generated by RNA Polymerase III.

Short Interspersed Elements (SINEs) are common in the genomes of most multicellular organisms. They are transcribed by RNA polymerase III from an internal promoter comprising boxes A and B. As transcripts of certain SINEs from mammalian genomes can be polyadenylated, such transcripts should contain the AATAAA sequence as well as those called β- and τ-signals. One of the goals of this work was to evaluate how autonomous and independent other SINE parts are β- and τ-signals. Extended regions outside of β- and τ-signals were deleted from SINEs B2 and Ves and the derived constructs were used to transfect HeLa cells in order to evaluate the relative levels of their transcripts as well as their polyadenylation efficiency. If the deleted regions affected boxes A and B, the 5'-flanking region of the U6 RNA gene with the external promoter was inserted upstream. Such substitution of the internal promoter in B2 completely restored its transcription. Almost all tested deletions/substitutions did not reduce the polyadenylation capacity of the transcripts, indicating a weak dependence of the function of β- and τ-signals on the neighboring sequences. A similar analysis of B2 and Ves constructs containing a 55-bp foreign sequence inserted between β- and τ-signals showed an equal polyadenylation efficiency of their transcripts compared to those of constructs without the insertion. The acquired poly(A)-tails significantly increased the lifetime and thus the cellular level of such transcripts. The data obtained highlight the potential of B2 and Ves SINEs as cassettes for the expression of relatively short sequences for various applications.

Elongation rate of RNA polymerase II affects pausing patterns across 3′ UTRs

Yeast mRNAs are polyadenylated at multiple sites in their 3' untranslated regions (3' UTRs), and poly(A) site usage is regulated by the rate of transcriptional elongation by RNA polymerase II (Pol II). Slow Pol II derivatives favor upstream poly(A) sites, and fast Pol II derivatives favor downstream poly(A) sites. Transcriptional elongation and polyadenylation are linked at the nucleotide level, presumably reflecting Pol II dwell time at each residue that influences the level of polyadenylation. Here, we investigate the effect of Pol II elongation rate on pausing patterns and the relationship between Pol II pause sites and poly(A) sites within 3' UTRs. Mutations that affect Pol II elongation rate alter sequence preferences at pause sites within 3' UTRs, and pausing preferences differ between 3' UTRs and coding regions. In addition, sequences immediately flanking the pause sites show preferences that are largely independent of Pol II speed. In wild-type cells, poly(A) sites are preferentially located< 50 nucleotides upstream from Pol II pause sites, but this spatial relationship is diminished in cells harboring Pol II speed mutants. Based on a random forest classifier, Pol II pause sites are modestly predicted by the distance to poly(A) sites but are better predicted by the chromatin landscape in Pol II speed derivatives. Transcriptional regulatory proteins can influence the relationship between Pol II pausing and polyadenylation but in a manner distinct from Pol II elongation rate derivatives. These results indicate a complex relationship between Pol II pausing and polyadenylation.

MRNA 3’UTR lengthening by alternative polyadenylation attenuates inflammatory responses and correlates with virulence of Influenza A virus

Changes of mRNA 3’UTRs by alternative polyadenylation (APA) have been associated to numerous pathologies, but the mechanisms and consequences often remain enigmatic. By combining transcriptomics, proteomics and recombinant viruses we show that all tested strains of IAV, including A/PR/8/34(H1N1) (PR8) and A/Cal/07/2009 (H1N1) (Cal09), cause APA. We mapped the effect to the highly conserved glycine residue at position 184 (G184) of the viral non-structural protein 1 (NS1). Unbiased mass spectrometry-based analyses indicate that NS1 causes APA by perturbing the function of CPSF4 and that this function is unrelated to virus-induced transcriptional shutoff. Accordingly, IAV strain PR8, expressing an NS1 variant with weak CPSF binding, does not induce host shutoff but only APA. However, recombinant IAV (PR8) expressing NS1(G184R) lacks binding to CPSF4 and thereby also the ability to cause APA. Functionally, the impaired ability to induce APA leads to an increased inflammatory cytokine production and an attenuated phenotype in a mouse infection model. Investigating diverse viral infection models showed that APA induction is a frequent ability of many pathogens. Collectively, we propose that targeting of the CPSF complex, leading to widespread alternative polyadenylation of host transcripts, constitutes a general immunevasion mechanism employed by a variety of pathogenic viruses.

The RNA-binding protein ZC3H11A interacts with the nuclear poly(A)-binding protein PABPN1 and alters polyadenylation of viral transcripts

Nuclear mRNA metabolism is regulated by multiple proteins, which either directly bind to RNA or form multi-protein complexes. The RNA-binding protein ZC3H11A is involved in nuclear mRNA export, NF-κB signaling and is essential during mouse embryo development. Furthermore, previous studies have shown that ZC3H11A is important for nuclear-replicating viruses. However, detailed biochemical characterization of the ZC3H11A protein has been lacking. In this study, we established the ZC3H11A protein interactome in human and mouse cells. We demonstrate that the nuclear poly(A)-binding protein PABPN1 interacts specifically with the ZC3H11A protein and controls ZC3H11A localization into nuclear speckles. We report that ZC3H11A specifically interacts with the human adenovirus type 5 (HAdV-5) capsid mRNA in a PABPN1-dependent manner. Notably, ZC3H11A uses the same zinc finger motifs to interact with PABPN1 and viral mRNA. Further, we demonstrate that the lack of ZC3H11A alters the polyadenylation of HAdV-5 capsid mRNA. Taken together, our results suggest that the ZC3H11A protein may act as a novel regulator of polyadenylation of nuclear mRNA.

Derepression of Y-linked multicopy protamine-like genes interferes with sperm nuclear compaction in D. melanogaster

Across species, sperm maturation involves the dramatic reconfiguration of chromatin into highly compact nuclei that enhance hydrodynamic ability and ensure paternal genomic integrity. This process is mediated by the replacement of histones by sperm nuclear basic proteins, also referred to as protamines. In humans, a carefully balanced dosage between two known protamine genes is required for optimal fertility. However, it remains unknown how their proper balance is regulated and how defects in balance may lead to compromised fertility. Here, we show that a nucleolar protein, modulo, a homolog of nucleolin, mediates the histone-to-protamine transition during Drosophila spermatogenesis. We find that modulo mutants display nuclear compaction defects during late spermatogenesis due to decreased expression of autosomal protamine genes (including Mst77F) and derepression of Y-linked multicopy Mst77F homologs (Mst77Y), leading to the mutant's known sterility. Overexpression of Mst77Y in a wild-type background is sufficient to cause nuclear compaction defects, similar to modulo mutant, indicating that Mst77Y is a dominant-negative variant interfering with the process of histone-to-protamine transition. Interestingly, ectopic overexpression of Mst77Y caused decompaction of X-bearing spermatids nuclei more frequently than Y-bearing spermatid nuclei, although this did not greatly affect the sex ratio of offspring. We further show that modulo regulates these protamine genes at the step of transcript polyadenylation. We conclude that the regulation of protamines mediated by modulo, ensuring the expression of functional ones while repressing dominant-negative ones, is critical for male fertility.

A post-transcriptional regulatory landscape of aging in the female mouse hippocampus.

Aging is associated with substantial physiological changes and constitutes a major risk factor for neurological disorders including dementia. Alterations in gene expression upon aging have been extensively studied; however, an in-depth characterization of post-transcriptional regulatory events remains elusive. Here, we profiled the age-related changes of the transcriptome and translatome in the female mouse hippocampus by RNA sequencing of total RNA and polysome preparations at four ages (3-, 6-, 12-, 20-month-old); and we implemented a variety of bioinformatics approaches to unravel alterations in transcript abundance, alternative splicing, and polyadenylation site selection. We observed mostly well-coordinated transcriptome and translatome expression signatures across age including upregulation of transcripts related to immune system processes and neuroinflammation, though transcripts encoding ribonucleoproteins or associated with mitochondrial functions, calcium signaling and the cell-cycle displayed substantial discordant profiles, suggesting translational control associated with age-related deficits in hippocampal-dependent behavior. By contrast, alternative splicing was less preserved, increased with age and was associated with distinct functionally-related transcripts encoding proteins acting at synapses/dendrites, RNA-binding proteins; thereby predicting regulatory roles for RBM3 and CIRBP. Only minor changes in polyadenylation site selection were identified, indicating pivotal 3'-end selection in young adults compared to older groups. Overall, our study provides a comprehensive resource of age-associated post-transcriptional regulatory events in the mouse hippocampus, enabling further examination of the molecular features underlying age-associated neurological diseases.

Improved Annotation of the Peach (Prunus persica) Genome and Identification of Tissue- or Development Stage-Specific Alternative Splicing through the Integration of Iso-Seq and RNA-Seq Data

Alternative splicing (AS) is an important way to generate notable regulatory and proteomic complexity in eukaryotes. However, accurate full-length splicing isoform discovery by second-generation sequencing (SGS) technologies is beset with the precise assembly of multiple isoforms from the same gene loci. In recent years, third-generation sequencing (TGS) technologies have been adopted to gain insight into different aspects of transcriptome complexity, such as complete sequences of mRNA, alternative splicing, fusion transcript, and alternative polyadenylation (APA). Here, we combined PacBio Iso-Seq and Illumina RNA-Seq technologies to decipher the full-length transcriptome of peach. In total, 40,477 nonredundant high-quality consensus transcript sequences were obtained from equally pooled libraries from 10 samples of 6 organs, including leaf, shoot, flower, fruit peel, fruit mesocarp, and fruit stone, of which 18,274 isoforms were novel isoforms of known genes and 546 isoforms were novel gene transcripts. We also discovered 148 fusion transcripts, 15,434 AS events, 508 potential lncRNAs, and 4368 genes with APA events. Of these AS events, the most abundant (62.48%) AS type was intron retention (IR). Moreover, the expression levels of different isoforms identified in this study were quantitatively evaluated, and highly tissue- or development stage-specific expression patterns were observed. The novel transcript isoforms and new characteristics of the peach transcriptome revealed by this study will facilitate the annotation of the peach genome and lay the foundations for functional research in the future.

Nucleotide-level linkage of transcriptional elongation and polyadenylation.

Alternative polyadenylation yields many mRNA isoforms whose 3' termini occur disproportionately in clusters within 3' untranslated regions. Previously, we showed that profiles of poly(A) site usage are regulated by the rate of transcriptional elongation by RNA polymerase (Pol) II (Geisberg et al., 2020). Pol II derivatives with slow elongation rates confer an upstream-shifted poly(A) profile, whereas fast Pol II strains confer a downstream-shifted poly(A) profile. Within yeast isoform clusters, these shifts occur steadily from one isoform to the next across nucleotide distances. In contrast, the shift between clusters - from the last isoform of one cluster to the first isoform of the next - is much less pronounced, even over large distances. GC content in a region 13-30 nt downstream from isoform clusters correlates with their sensitivity to Pol II elongation rate. In human cells, the upstream shift caused by a slow Pol II mutant also occurs continuously at single nucleotide resolution within clusters but not between them. Pol II occupancy increases just downstream of poly(A) sites, suggesting a linkage between reduced elongation rate and cluster formation. These observations suggest that (1) Pol II elongation speed affects the nucleotide-level dwell time allowing polyadenylation to occur, (2) poly(A) site clusters are linked to the local elongation rate, and hence do not arise simply by intrinsically imprecise cleavage and polyadenylation of the RNA substrate, (3) DNA sequence elements can affect Pol II elongation and poly(A) profiles, and (4) the cleavage/polyadenylation and Pol II elongation complexes are spatially, and perhaps physically, coupled so that polyadenylation occurs rapidly upon emergence of the nascent RNA from the Pol II elongation complex.

Genomic organization, intragenic tandem duplication, and expression analysis of chicken TGFBR2 gene

Transforming growth factor beta receptor Ⅱ (TGFBR2), a core member of the transforming growth factor-β (TGF-β) signaling pathway. To date, chicken TGFBR2 (cTGFBR2) genomic structure has not been fully explored. Here, the complete sequences of cTGFBR2 transcript isoforms were determined by 5′ and 3′ rapid amplification of cDNA ends (5′ & 3′ RACE) and reverse transcription polymerase chain reaction (RT-PCR); the tissue expression profiling of cTGFBR2 transcript isoforms was performed using quantitative real-time polymerase chain reaction (qRT-PCR). The results showed that cTGFBR2 gene produced 3 transcript isoforms though alternative transcription initiation, splicing, and polyadenylation, which were designated as cTGFBR2-1, cTGFBR2-2, and cTGFBR2-3, respectively. These 3 cTGFBR2 transcript isoforms encoded 3 protein isoforms: cTGFBR2-1, cTGFBR2-2, and cTGFBR2-3. Duplication analysis revealed that, unlike other animal species, cTGFBR2 gene harbored a 5.5-kb intragenic tandem duplication. Tissue expression profiling in the 4-wk-old Arbor Acres (AA) broiler chickens showed that cTGFBR2-1 was ubiquitously expressed, with high expression in abdominal fat, subcutaneous fat, lung, gizzard, and muscle; cTGFBR2-2 was highly expressed in heart, kidney, gizzard, and muscle; cTGFBR2-3 was weakly expressed in all the tested chicken tissues. Tissue expression profiling in the 7-wk-old broiler chickens of the fat and lean lines of Northeast Agricultural University broiler lines divergently selected for abdominal fat content (NEAUHLF) showed that cTGFBR2-1 was significantly differentially expressed in all the tested tissues except heart, cTGFBR2-2 was significantly differentially expressed in all the tested tissues except subcutaneous fat and liver, and cTGFBR2-3 was significantly differentially expressed in all the tested tissues between the lean and fat lines. Intriguingly, in the fat line, the 3 cTGFBR2 transcript isoforms were expressed to varying degrees in all the 3 tested fat tissues, while in the lean line, only cTGFBR2-1 was expressed in all the 3 tested fat tissues. This is the first report of intragenic tandem duplication within TGFBR2 gene. Our findings pave the way for further studies on the functions and regulation of cTGFBR2 gene.

CDK9 and PP2A regulate RNA polymerase II transcription termination and coupled RNA maturation

CDK9 is a kinase critical for the productive transcription of protein‐coding genes by RNA polymerase II (pol II). As part of P‐TEFb, CDK9 phosphorylates the carboxyl‐terminal domain (CTD) of pol II and elongation factors, which allows pol II to elongate past the early elongation checkpoint (EEC) encountered soon after initiation. We show that, in addition to halting pol II at the EEC, loss of CDK9 activity causes premature termination of transcription across the last exon, loss of polyadenylation factors from chromatin, and loss of polyadenylation of nascent transcripts. Inhibition of the phosphatase PP2A abrogates the premature termination and loss of polyadenylation caused by CDK9 inhibition, indicating that this kinase/phosphatase pair regulates transcription elongation and RNA processing at the end of protein‐coding genes. We also confirm the splicing factor SF3B1 as a target of CDK9 and show that SF3B1 in complex with polyadenylation factors is lost from chromatin after CDK9 inhibition. These results emphasize the important roles that CDK9 plays in coupling transcription elongation and termination to RNA maturation downstream of the EEC.

Zika Virus Infection Alters Gene Expression and Poly-Adenylation Patterns in Placental Cells.

Flaviviruses are small RNA viruses that are mainly transmitted via arthropod vectors and are found in tropic and sub-tropical regions. Most infections are asymptomatic (90–95%), but symptoms can be as severe as hemorrhagic fever and encephalitis. One recently emerged flavivirus is Zika virus (ZIKV), which was originally isolated from rhesus monkeys in Uganda roughly 70 years ago but has recently spread east, reaching S. America in 2015–2016. This outbreak was associated with the development of Guillain–Barré syndrome in adults and microcephaly in infants born to expectant mothers infected early in pregnancy. ZIKV must traverse the placenta to impact the development of the fetus, but the mechanisms responsible are unknown. While flaviviruses are known to disrupt splicing patterns in host cells, little is known about how flaviviruses such as ZIKV impact the alternative polyadenylation (APA) of host transcripts. This is important as APA is well-established as a mechanism in the regulation of mRNA metabolism and translation. Thus, we sought to characterize transcriptomic changes including APA in human placental (JEG3) cells in response to ZIKV infection using Poly(A)-ClickSeq (PAC-Seq). We used our differential Poly(A)-cluster (DPAC) analysis pipeline to characterize changes in differential gene expression, alternative poly-adenylation (APA) and the use of alternative terminal exons. We identified 98 upregulated genes and 28 downregulated genes. Pathway enrichment analysis indicated that many RNA processing and immune pathways were upregulated in ZIKV-infected JEG3 cells. We also updated DPAC to provide additional metrics of APA including the percentage-distal usage index (PDUI), which revealed that APA was extensive and the 3′ UTRs of 229 genes were lengthened while 269 were shortened. We further found that there were 214 upregulated and 59 downregulated poly(A)-clusters (PACs). We extracted the nucleotide sequences surrounding these PACs and found that the canonical signals for poly-adenylation (binding site for poly-A binding protein (PABP) upstream and a GU-rich region down-stream of the PAC) were only enriched in the downregulated PACs. These results indicate that ZIKV infection makes JEG3 cells more permissive to non-canonical poly-adenylation signals.

DNA- and RNA-Binding Proteins Linked Transcriptional Control and Alternative Splicing Together in a Two-Layer Regulatory Network System of Chronic Myeloid Leukemia.

DNA- and RNA-binding proteins (DRBPs) typically possess multiple functions to bind both DNA and RNA and regulate gene expression from more than one level. They are controllers for post-transcriptional processes, such as splicing, polyadenylation, transportation, translation, and degradation of RNA transcripts in eukaryotic organisms, as well as regulators on the transcriptional level. Although DRBPs are reported to play critical roles in various developmental processes and diseases, it is still unclear how they work with DNAs and RNAs simultaneously and regulate genes at the transcriptional and post-transcriptional levels. To investigate the functional mechanism of DRBPs, we collected data from a variety of databases and literature and identified 118 DRBPs, which function as both transcription factors (TFs) and splicing factors (SFs), thus called DRBP-SF. Extensive investigations were conducted on four DRBP-SFs that were highly expressed in chronic myeloid leukemia (CML), heterogeneous nuclear ribonucleoprotein K (HNRNPK), heterogeneous nuclear ribonucleoprotein L (HNRNPL), non-POU domain–containing octamer–binding protein (NONO), and TAR DNA-binding protein 43 (TARDBP). By integrating and analyzing ChIP-seq, CLIP-seq, RNA-seq, and shRNA-seq data in K562 using binding and expression target analysis and Statistical Utility for RBP Functions, we discovered a two-layer regulatory network system centered on these four DRBP-SFs and proposed three possible regulatory models where DRBP-SFs can connect transcriptional and alternative splicing regulatory networks cooperatively in CML. The exploration of the identified DRBP-SFs provides new ideas for studying DRBP and regulatory networks, holding promise for further mechanistic discoveries of the two-layer gene regulatory system that may play critical roles in the occurrence and development of CML.

DNA Damage Regulates the Functions of the RNA Binding Protein Sam68 through ATM-Dependent Phosphorylation.

Simple SummaryAlterations of the complex network of interactions between the DNA damage response pathway and RNA metabolism have been described in several tumors, and increasing efforts are devoted to the elucidation of the molecular mechanisms involved in this network. Previous large-scale proteomic studies identified the RNA binding protein Sam68 as a putative target of the ATM kinase. Herein, we demonstrate that ATM phosphorylates Sam68 upon DNA damage induction, and this post-translational modification regulates both the signaling function of Sam68 in the initial phase of the DNA damage response and its RNA processing activity. Thus, our study uncovers anew crosstalk between ATM and Sam68, which may represent a paradigm for the functional interaction between the DDR pathway and RNA binding proteins, and a possible actionabletarget in human cancers.Cancer cells frequently exhibit dysregulation of the DNA damage response (DDR), genomic instability, and altered RNA metabolism. Recent genome-wide studies have strongly suggested an interaction between the pathways involved in the cellular response to DDR and in the regulation of RNA metabolism, but the molecular mechanism(s) involved in this crosstalk are largely unknown. Herein, we found that activation of the DDR kinase ATM promotes its interaction with Sam68, leading to phosphorylation of this multifunctional RNA binding protein (RBP) on three residues: threonine 61, serine 388 and serine 390. Moreover, we demonstrate that ATM-dependent phosphorylation of threonine 61 promotes the function of Sam68 in the DDR pathway and enhances its RNA processing activity. Importantly, ATM-mediated phosphorylation of Sam68 in prostate cancer cells modulates alternative polyadenylation of transcripts that are targets of Sam68, supporting the notion that the ATM–Sam68 axis exerts a multifaceted role in the response to DNA damage. Thus, our work validates Sam68 as an ATM kinase substrate and uncovers an unexpected bidirectional interplay between ATM and Sam68, which couples the DDR pathway to modulation of RNA metabolism in response to genotoxic stress.

A high-resolution single-molecule sequencing-based Arabidopsis transcriptome using novel methods of Iso-seq analysis

BackgroundAccurate and comprehensive annotation of transcript sequences is essential for transcript quantification and differential gene and transcript expression analysis. Single-molecule long-read sequencing technologies provide improved integrity of transcript structures including alternative splicing, and transcription start and polyadenylation sites. However, accuracy is significantly affected by sequencing errors, mRNA degradation, or incomplete cDNA synthesis.ResultsWe present a new and comprehensive Arabidopsis thaliana Reference Transcript Dataset 3 (AtRTD3). AtRTD3 contains over 169,000 transcripts—twice that of the best current Arabidopsis transcriptome and including over 1500 novel genes. Seventy-eight percent of transcripts are from Iso-seq with accurately defined splice junctions and transcription start and end sites. We develop novel methods to determine splice junctions and transcription start and end sites accurately. Mismatch profiles around splice junctions provide a powerful feature to distinguish correct splice junctions and remove false splice junctions. Stratified approaches identify high-confidence transcription start and end sites and remove fragmentary transcripts due to degradation. AtRTD3 is a major improvement over existing transcriptomes as demonstrated by analysis of an Arabidopsis cold response RNA-seq time-series. AtRTD3 provides higher resolution of transcript expression profiling and identifies cold-induced differential transcription start and polyadenylation site usage.ConclusionsAtRTD3 is the most comprehensive Arabidopsis transcriptome currently. It improves the precision of differential gene and transcript expression, differential alternative splicing, and transcription start/end site usage analysis from RNA-seq data. The novel methods for identifying accurate splice junctions and transcription start/end sites are widely applicable and will improve single-molecule sequencing analysis from any species.