Poly(A)-ClickSeq: click-chemistry for next-generation 3΄-end sequencing without RNA enrichment or fragmentation

InPAS: An R/Bioconductor Package for Identifying Novel Polyadenylation Sites and Alternative Polyadenylation from Bulk RNA-seq Data

BACKGROUND

Alternative cleavage and polyadenylation (APA) is a crucial post-transcriptional gene regulation mechanism that regulates gene expression in eukaryotes by increasing the diversity and complexity of both the transcriptome and proteome. Despite the development of more than a dozen experimental methods over the last decade to identify and quantify APA events, widespread adoption of these methods has been limited by technical, financial, and time constraints. Consequently, APA remains poorly understood in most eukaryotes. However, RNA sequencing (RNA-seq) technology has revolutionized transcriptome profiling and recent studies have shown that RNA-seq data can be leveraged to identify and quantify APA events.

RESULTS

To fully capitalize on the exponentially growing RNA-seq data, we developed InPAS (Identification of Novel alternative PolyAdenylation Sites), an R/Bioconductor package for accurate identification of novel and known cleavage and polyadenylation sites (CPSs), as well as quantification of APA from RNA-seq data of various experimental designs. Compared to other APA analysis tools, InPAS offers several important advantages, including the ability to detect both novel proximal and distal CPSs, to fine tune positions of CPSs using a naïve Bayes classifier based on flanking sequence features, and to identify APA events from RNA-seq data of complex experimental designs using linear models. We benchmarked the performance of InPAS and other leading tools using simulated and experimental RNA-seq data with matched 3'-end RNA-seq data. Our results reveal that InPAS frequently outperforms existing tools in terms of precision, sensitivity, and specificity. Furthermore, we demonstrate its scalability and versatility by applying it to large, diverse RNA-seq datasets.

CONCLUSIONS

InPAS is an efficient and robust tool for identifying and quantifying APA events using readily accessible conventional RNA-seq data. Its versatility opens doors to explore APA regulation across diverse eukaryotic systems with various experimental designs. We believe that InPAS will drive APA research forward, deepening our understanding of its role in regulating gene expression, and potentially leading to the discovery of biomarkers or therapeutics for diseases.

SARS-CoV-2 EndoU-ribonuclease regulates RNA recombination and impacts viral fitness

Coronaviruses (CoVs) maintain large RNA genomes that frequently undergoes mutations and recombination, contributing to their evolution and emergence. In this study, we find that SARS-CoV-2 has greater RNA recombination frequency than other human CoVs. In addition, coronavirus RNA recombination primarily occurs at uridine (U)-enriched RNA sequences. Therefore, we next evaluated the role of SARS-CoV-2 NSP15, a viral endonuclease that targets uridines (EndoU), in RNA recombination and virus infection. Using a catalytically inactivated EndoU mutant (NSP15H234A), we observed attenuated viral replication in vitro and in vivo. However, the loss of EndoU activity also dysregulated inflammation resulting in similar disease in vivo despite reduced viral loads. Next-generation sequencing (NGS) demonstrated that loss of EndoU activity disrupts SARS-CoV-2 RNA recombination by reducing viral sub-genomic message but increasing recombination events that contribute to defective viral genomes (DVGs). Overall, the study demonstrates that NSP15 plays a critical role in regulating RNA recombination and SARS-CoV-2 pathogenesis.

Regulation of the Drosophila transcriptome by Pumilio and the CCR4-NOT deadenylase complex

The sequence-specific RNA-binding protein Pumilio (Pum) controls Drosophila development; however, the network of mRNAs that it regulates remains incompletely characterized. In this study, we use knockdown and knockout approaches coupled with RNA-seq to measure the impact of Pum on the transcriptome of Drosophila cells in culture. We also use an improved RNA coimmunoprecipitation method to identify Pum-bound mRNAs in Drosophila embryos. Integration of these data sets with the locations of Pum-binding motifs across the transcriptome reveals novel direct Pum target genes involved in neural, muscle, wing, and germ cell development and in cellular proliferation. These genes include components of Wnt, TGF-β, MAPK/ERK, and Notch signaling pathways, DNA replication, and lipid metabolism. We identify the mRNAs regulated by the CCR4-NOT deadenylase complex, a key factor in Pum-mediated repression, and observe concordant regulation of Pum:CCR4-NOT target mRNAs. Computational modeling reveals that Pum binding, binding site number, clustering, and sequence context are important determinants of regulation. In contrast, we show that the responses of direct mRNA targets to Pum-mediated repression are not influenced by the content of optimal synonymous codons. Moreover, contrary to a prevailing model, we do not detect a role for CCR4-NOT in the degradation of mRNAs with low codon optimality. Together, the results of this work provide new insights into the Pum regulatory network and mechanisms and the parameters that influence the efficacy of Pum-mediated regulation.

PolyAMiner-Bulk is a deep learning-based algorithm that decodes alternative polyadenylation dynamics from bulk RNA-seq data

Alternative polyadenylation (APA) is a key post-transcriptional regulatory mechanism; yet, its regulation and impact on human diseases remain understudied. Existing bulk RNA sequencing (RNA-seq)-based APA methods predominantly rely on predefined annotations, severely impacting their ability to decode novel tissue- and disease-specific APA changes. Furthermore, they only account for the most proximal and distal cleavage and polyadenylation sites (C/PASs). Deconvoluting overlapping C/PASs and the inherent noisy 3' UTR coverage in bulk RNA-seq data pose additional challenges. To overcome these limitations, we introduce PolyAMiner-Bulk, an attention-based deep learning algorithm that accurately recapitulates C/PAS sequence grammar, resolves overlapping C/PASs, captures non-proximal-to-distal APA changes, and generates visualizations to illustrate APA dynamics. Evaluation on multiple datasets strongly evinces the performance merit of PolyAMiner-Bulk, accurately identifying more APA changes compared with other methods. With the growing importance of APA and the abundance of bulk RNA-seq data, PolyAMiner-Bulk establishes a robust paradigm of APA analysis.

Pan-cancer analysis of NUDT21 and its effect on the proliferation of human head and neck squamous cell carcinoma

BACKGROUND

Based on bioinformatics research of NUDT21 in pan-cancer, we aimed to clarify the mechanism of NUDT21 in HHNC by experiment.

METHODS

The correlation between differential expression of NUDT21 in pan-cancer and survival prognosis, genomic instability, tumor stemness, DNA repair, RNA methylation and with immune microenvironment were analyzed by the application of different pan-cancer analysis web databases. In addition, immunohistochemistry staining and genetic detection of NUDT21 in HHNCC tumor tissues by immunohistochemistry and qRT-PCR. Then, through in vitro cell experiments, NUDT21 was knocked down by lentivirus to detect the proliferation, cycle, apoptosis of FaDu and CNE-2Z cells, and finally by PathScan intracellular signaling array reagent to detect the apoptotic protein content.

RESULTS

Based on the pan-cancer analysis, we found that elevated expression of NUDT21 in most cancers was significantly correlated with TMB, MSI, neoantigens and chromosomal ploidy, and in epigenetics, elevated NUDT21 expression was strongly associated with genomic stability, mismatch repair genes, tumor stemness, and RNA methylation. Based on immunosuppressive score, we found that NUDT21 plays an essential role in the immunosuppressive environment by suppressing immune checkpointing effect in most cancers. In addition, using HHNSCC as a study target, PCR and pathological detection of NUDT21 in tumor tissues was significantly increased than that in paracancerous normal tissues. In vitro cellular assays, silencing NUDT21 inhibited proliferation and promoted apoptosis in FaDu and CNE-2Z cells, and blocked the cell cycle in the G2/M phase. Therefore, the experiments confirmed that NUDT21 promotes the proliferation of FaDu by suppressing the expression of apoptotic.

Regulation of the Drosophila transcriptome by Pumilio and CCR4-NOT deadenylase

The sequence-specific RNA-binding protein Pumilio controls development of Drosophila; however, the network of mRNAs that it regulates remains incompletely characterized. In this study, we utilize knockdown and knockout approaches coupled with RNA-Seq to measure the impact of Pumilio on the transcriptome of Drosophila cells. We also used an improved RNA co-immunoprecipitation method to identify Pumilio bound mRNAs in Drosophila embryos. Integration of these datasets with the content of Pumilio binding motifs across the transcriptome revealed novel direct Pumilio target genes involved in neural, muscle, wing, and germ cell development, and cellular proliferation. These genes include components of Wnt, TGF-beta, MAPK/ERK, and Notch signaling pathways, DNA replication, and lipid metabolism. Additionally, we identified the mRNAs regulated by the CCR4-NOT deadenylase complex, a key factor in Pumilio-mediated repression, and observed concordant regulation of Pumilio:CCR4-NOT target mRNAs. Computational modeling revealed that Pumilio binding, binding site number, density, and sequence context are important determinants of regulation. Moreover, the content of optimal synonymous codons in target mRNAs exhibits a striking functional relationship to Pumilio and CCR4-NOT regulation, indicating that the inherent translation efficiency and stability of the mRNA modulates their response to these trans-acting regulatory factors. Together, the results of this work provide new insights into the Pumilio regulatory network and mechanisms, and the parameters that influence the efficacy of Pumilio-mediated regulation.

NUDT21 alters glioma migration through differential alternative polyadenylation of LAMC1

PURPOSE

Gliomas and their surrounding microenvironment constantly interact to promote tumorigenicity, yet the underlying posttranscriptional regulatory mechanisms that govern this interplay are poorly understood.

METHODS

Utilizing our established PAC-seq approach and PolyAMiner bioinformatic analysis pipeline, we deciphered the NUDT21-mediated differential APA dynamics in glioma cells.

RESULTS

We identified LAMC1 as a critical NUDT21 alternative polyadenylation (APA) target, common in several core glioma-driving signaling pathways. qRT-PCR analysis confirmed that NUDT21-knockdown in glioma cells results in the preferred usage of the proximal polyA signal (PAS) of LAMC1. Functional studies revealed that NUDT21-knockdown-induced 3'UTR shortening of LAMC1 is sufficient to cause translational gain, as LAMC1 protein is upregulated in these cells compared to their respective controls. We demonstrate that 3'UTR shortening of LAMC1 after NUDT21 knockdown removes binding sites for miR-124/506, thereby relieving potent miRNA-based repression of LAMC1 expression. Remarkably, we report that the knockdown of NUDT21 significantly promoted glioma cell migration and that co-depletion of LAMC1 with NUDT21 abolished this effect. Lastly, we observed that LAMC1 3'UTR shortening predicts poor prognosis of low-grade glioma patients from The Cancer Genome Atlas.

CONCLUSION

This study identifies NUDT21 as a core alternative polyadenylation factor that regulates the tumor microenvironment through differential APA and loss of miR-124/506 inhibition of LAMC1. Knockdown of NUDT21 in GBM cells mediates 3'UTR shortening of LAMC1, contributing to an increase in LAMC1, increased glioma cell migration/invasion, and a poor prognosis.

Alternative polyadenylation regulation in cardiac development and cardiovascular disease

Cleavage and polyadenylation of pre-mRNAs is a necessary step for gene expression and function. Majority of human genes exhibit multiple polyadenylation sites, which can be alternatively used to generate different mRNA isoforms from a single gene. Alternative polyadenylation (APA) of pre-mRNAs is important for the proteome and transcriptome landscape. APA is tightly regulated during development and contributes to tissue-specific gene regulation. Mis-regulation of APA is linked to a wide range of pathological conditions. APA-mediated gene regulation in the heart is emerging as a new area of research. Here, we will discuss the impact of APA on gene regulation during heart development and in cardiovascular diseases. First, we will briefly review how APA impacts gene regulation and discuss molecular mechanisms that control APA. Then, we will address APA regulation during heart development and its dysregulation in cardiovascular diseases. Finally, we will discuss pre-mRNA targeting strategies to correct aberrant APA patterns of essential genes for the treatment or prevention of cardiovascular diseases. The RNA field is blooming due to advancements in RNA-based technologies. RNA-based vaccines and therapies are becoming the new line of effective and safe approaches for the treatment and prevention of human diseases. Overall, this review will be influential for understanding gene regulation at the RNA level via APA in the heart and will help design RNA-based tools for the treatment of cardiovascular diseases in the future.

The peptide woods are lovely, dark and deep: Hunting for novel cancer antigens

Harnessing the patient's immune system to control a tumor is a proven avenue for cancer therapy. T cell therapies as well as therapeutic vaccines, which target specific antigens of interest, are being explored as treatments in conjunction with immune checkpoint blockade. For these therapies, selecting the best suited antigens is crucial. Most of the focus has thus far been on neoantigens that arise from tumor-specific somatic mutations. Although there is clear evidence that T-cell responses against mutated neoantigens are protective, the large majority of these mutations are not immunogenic. In addition, most somatic mutations are unique to each individual patient and their targeting requires the development of individualized approaches. Therefore, novel antigen types are needed to broaden the scope of such treatments. We review high throughput approaches for discovering novel tumor antigens and some of the key challenges associated with their detection, and discuss considerations when selecting tumor antigens to target in the clinic.

U1 snRNP increases RNA Pol II elongation rate to enable synthesis of long genes

The expansion of introns within mammalian genomes poses a challenge for the production of full-length messenger RNAs (mRNAs), with increasing evidence that these long AT-rich sequences present obstacles to transcription. Here, we investigate RNA polymerase II (RNAPII) elongation at high resolution in mammalian cells and demonstrate that RNAPII transcribes faster across introns. Moreover, we find that this acceleration requires the association of U1 snRNP (U1) with the elongation complex at 5' splice sites. The role of U1 to stimulate elongation rate through introns reduces the frequency of both premature termination and transcriptional arrest, thereby dramatically increasing RNA production. We further show that changes in RNAPII elongation rate due to AT content and U1 binding explain previous reports of pausing or termination at splice junctions and the edge of CpG islands. We propose that U1-mediated acceleration of elongation has evolved to mitigate the risks that long AT-rich introns pose to transcript completion.

Relationship between NUDT21 mediated alternative polyadenylation process and tumor

Alternative polyadenylation (APA) is a molecular process that generates diversity at the 3' end of RNA polymerase II transcripts from over 60% of human genes. APA and microRNA regulation are both mechanisms of post-transcriptional regulation of gene expression. As a key molecular mechanism, Alternative polyadenylation often results in mRNA isoforms with the same coding sequence but different lengths of 3' UTRs, while microRNAs regulate gene expression by binding to specific mRNA 3' UTRs. Nudix Hydrolase 21 (NUDT21) is a crucial mediator involved in alternative polyadenylation (APA). Different studies have reported a dual role of NUDT21 in cancer (both oncogenic and tumor suppressor). The present review focuses on the functions of APA, miRNA and their interaction and roles in development of different types of tumors.NUDT21 mediated 3' UTR-APA changes can be used to generate specific signatures that can be used as potential biomarkers in development and disease. Due to the emerging role of NUDT21 as a regulator of the aforementioned RNA processing events, modulation of NUDT21 levels may be a novel viable therapeutic approach.

Quantitative trait locus mapping and identification of candidate genes for resistance to Verticillium wilt in four recombinant inbred line populations of Gossypium hirsutum

Improving resistance to Verticillium wilt is of great significance for achieving high and stable yields of Upland cotton (Gossypium hirsutum). To deeply understand the genetic basis of cotton resistance to Verticillium wilt, Verticillium wilt-resistant Upland Lumianyan 28 and four Verticillium wilt-susceptible Acala cotton cultivars were used to create four recombinant inbred line (RIL) populations of 469 families through nested hybridization. Phenotypic data collected in five stressful environments were used to select resistant and sensitive lines and create a mixed pool of extreme phenotypes for BSA-seq. A total of 8 QTLs associated with Verticillium wilt resistance were identified on 4 chromosomes, of which qVW-A12-5 was detected simultaneously in the RIL populations and in one of the RIL populations and was identified for the first time. According to the sequence comparison and transcriptome analysis of candidate genes in the QTL interval between parents and pools, 4 genes were identified in the qVW-A12-5 interval. qRT-PCR of parental and phenotypically extreme lines revealed that Gh_CPR30 was induced by and may be a candidate gene for resistance to Verticillium wilt in G. hirsutum. Furthermore, VIGS technology revealed that the disease severity index (DSI) of the Gh_CPR30-silenced plants was significantly higher than that of the control. These results indicate that the Gh_CPR30 gene plays an important role in the resistance of G. hirsutum to Verticillium wilt, and the study provides a molecular basis for analyzing the molecular mechanism underlying G. hirsutum resistance to Verticillium wilt.

Macrophage differentiation is marked by increased abundance of the mRNA 3' end processing machinery, altered poly(A) site usage, and sensitivity to the level of CstF64

Regulation of mRNA polyadenylation is important for response to external signals and differentiation in several cell types, and results in mRNA isoforms that vary in the amount of coding sequence or 3' UTR regulatory elements. However, its role in differentiation of monocytes to macrophages has not been investigated. Macrophages are key effectors of the innate immune system that help control infection and promote tissue-repair. However, overactivity of macrophages contributes to pathogenesis of many diseases. In this study, we show that macrophage differentiation is characterized by shortening and lengthening of mRNAs in relevant cellular pathways. The cleavage/polyadenylation (C/P) proteins increase during differentiation, suggesting a possible mechanism for the observed changes in poly(A) site usage. This was surprising since higher C/P protein levels correlate with higher proliferation rates in other systems, but monocytes stop dividing after induction of differentiation. Depletion of CstF64, a C/P protein and known regulator of polyadenylation efficiency, delayed macrophage marker expression, cell cycle exit, attachment, and acquisition of structural complexity, and impeded shortening of mRNAs with functions relevant to macrophage biology. Conversely, CstF64 overexpression increased use of promoter-proximal poly(A) sites and caused the appearance of differentiated phenotypes in the absence of induction. Our findings indicate that regulation of polyadenylation plays an important role in macrophage differentiation.

The Transition from Cancer "omics" to "epi-omics" through Next- and Third-Generation Sequencing

Deciphering cancer etiopathogenesis has proven to be an especially challenging task since the mechanisms that drive tumor development and progression are far from simple. An astonishing amount of research has revealed a wide spectrum of defects, including genomic abnormalities, epigenomic alterations, disturbance of gene transcription, as well as post-translational protein modifications, which cooperatively promote carcinogenesis. These findings suggest that the adoption of a multidimensional approach can provide a much more precise and comprehensive picture of the tumor landscape, hence serving as a powerful tool in cancer research and precision oncology. The introduction of next- and third-generation sequencing technologies paved the way for the decoding of genetic information and the elucidation of cancer-related cellular compounds and mechanisms. In the present review, we discuss the current and emerging applications of both generations of sequencing technologies, also referred to as massive parallel sequencing (MPS), in the fields of cancer genomics, transcriptomics and proteomics, as well as in the progressing realms of epi-omics. Finally, we provide a brief insight into the expanding scope of sequencing applications in personalized cancer medicine and pharmacogenomics.

Click Chemistry Enables Rapid Amplification of Full-Length Reverse Transcripts for Long-Read Third Generation Sequencing

Here we describe the development of a novel click chemistry-based method for the generation and amplification of full-length cDNA libraries from total RNA, while avoiding the need for problematic template-switching (TS) reactions. Compared with prior efforts, our method involves neither random priming nor stochastic cDNA termination, thus enabling amplification of transcripts that were previously inaccessible via related click chemistry-based RNA sequencing techniques. A key modification involving the use of PCR primers containing two overhanging 3'-nucleotides substantially improved the read-through compatibility of the 1,4-disubstituted 1,2,3-triazole-containing cDNA, where such modifications typically hinder amplification. This allowed us to more than double the possible insert size compared with the state-of-the art click chemistry-based technique, PAC-seq. Furthermore, our method performed on par with a commercially available PCR-cDNA RNA sequencing kit, as determined by Oxford Nanopore sequencing. Given the known advantages of PAC-seq, namely, suppression of PCR artifacts, we anticipate that our contribution could enable diverse applications including improved analyses of mRNA splicing variants and fusion transcripts.

Drought induces epitranscriptome and proteome changes in stem-differentiating xylem of Populus trichocarpa

Understanding gene expression and regulation requires insights into RNA transcription, processing, modification, and translation. However, the relationship between the epitranscriptome and the proteome under drought stress remains undetermined in poplar (Populus trichocarpa). In this study, we used Nanopore direct RNA sequencing and tandem mass tag-based proteomic analysis to examine epitranscriptomic and proteomic regulation induced by drought treatment in stem-differentiating xylem (SDX). Our results revealed a decreased full-length read ratio under drought treatment and, especially, a decreased association between transcriptome and proteome changes in response to drought. Epitranscriptome analysis of cellulose- and lignin-related genes revealed an increased N6-Methyladenosine (m6A) ratio, which was accompanied by decreased RNA abundance and translation, under drought stress. Interestingly, usage of the distal poly(A) site increased during drought stress. Finally, we found that transcripts of highly expressed genes tend to have shorter poly(A) tail length (PAL), and drought stress increased the percentage of transcripts with long PAL. These findings provide insights into the interplay among m6A, polyadenylation, PAL, and translation under drought stress in P. trichocarpa SDX.

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.

QTQTN motif upstream of the furin-cleavage site plays a key role in SARS-CoV-2 infection and pathogenesis

The furin cleavage site (FCS), an unusual feature in the SARS-CoV-2 spike protein, has been spotlighted as a factor key to facilitating infection and pathogenesis by increasing spike processing. Similarly, the QTQTN motif directly upstream of the FCS is also an unusual feature for group 2B coronaviruses (CoVs). The QTQTN deletion has consistently been observed in in vitro cultured virus stocks and some clinical isolates. To determine whether the QTQTN motif is critical to SARS-CoV-2 replication and pathogenesis, we generated a mutant deleting the QTQTN motif (ΔQTQTN). Here, we report that the QTQTN deletion attenuates viral replication in respiratory cells in vitro and attenuates disease in vivo. The deletion results in a shortened, more rigid peptide loop that contains the FCS and is less accessible to host proteases, such as TMPRSS2. Thus, the deletion reduced the efficiency of spike processing and attenuates SARS-CoV-2 infection. Importantly, the QTQTN motif also contains residues that are glycosylated, and disruption of its glycosylation also attenuates virus replication in a TMPRSS2-dependent manner. Together, our results reveal that three aspects of the S1/S2 cleavage site-the FCS, loop length, and glycosylation-are required for efficient SARS-CoV-2 replication and pathogenesis.

CFIm25 regulates human stem cell function independently of its role in mRNA alternative polyadenylation

It has recently been shown that CFIm25, a canonical mRNA 3’ processing factor, could play a variety of physiological roles through its molecular function in the regulation of mRNA alternative polyadenylation (APA). Here, we used CRISPR/Cas9-mediated gene editing approach in human embryonic stem cells (hESCs) for CFIm25, and obtained three gene knockdown/mutant cell lines. CFIm25 gene editing resulted in higher proliferation rate and impaired differentiation potential for hESCs, with these effects likely to be directly regulated by the target genes, including the pluripotency factor rex1. Mechanistically, we unexpected found that perturbation in CFIm25 gene expression did not significantly affect cellular mRNA 3’ processing efficiency and APA profile. Rather, we provided evidences that CFIm25 may impact RNA polymerase II (RNAPII) occupancy at the body of transcribed genes, and promote the expression level of a group of transcripts associated with cellular proliferation and/or differentiation. Taken together, these results reveal novel mechanisms underlying CFIm25ʹs modulation in determination of cell fate, and provide evidence that the process of mammalian gene transcription may be regulated by an mRNA 3’ processing factor.

Suppression of premature transcription termination leads to reduced mRNA isoform diversity and neurodegeneration

Tight regulation of mRNA isoform expression is essential for neuronal development, maintenance, and function; however, the repertoire of proteins that govern isoform composition and abundance remains incomplete. Here, we show that the RNA kinase CLP1 regulates mRNA isoform expression through suppression of proximal cleavage and polyadenylation. We found that human stem-cell-derived motor neurons without CLP1 or with the disease-associated CLP1 p.R140H variant had distinct patterns of RNA-polymerase-II-associated cleavage and polyadenylation complex proteins that correlated with polyadenylation site usage. These changes resulted in imbalanced mRNA isoform expression of long genes important for neuronal function that were recapitulated in vivo. Strikingly, we observed the same pattern of reduced mRNA isoform diversity in 3' end sequencing data from brain tissues of patients with neurodegenerative disease. Together, our results identify a previously uncharacterized role for CLP1 in mRNA 3' end formation and reveal an mRNA misprocessing signature in neurodegeneration that may suggest a common mechanism of disease.

Chikungunya Virus Infects the Heart and Induces Heart-Specific Transcriptional Changes in an Immunodeficient Mouse Model of Infection

Chikungunya virus (CHIKV) is a mosquito-transmitted pathogen in family Togaviridae, genus Alphavirus. Although CHIKV is well known for its ability to cause debilitating rheumatoid-like arthritis, it has been also been observed to cause cardiovascular symptoms such as arrhythmias. Here, using samples from a previous study, we sequenced RNA from serum, kidney, skeletal muscle, and cardiac muscle from CHIKV- and mock-infected IFN-αR-/- mice using two sequencing techniques to investigate heart-specific changes in virus mutational profiles and host gene expression. Mutation rates were similar across muscle tissues although heart tissue carried heart-specific CHIKV minority variants, one of which had a coding change in the nsP3 gene and another in the 3'UTR. Importantly, heart-specific transcriptional changes included differential expression of genes critical for ion transport and muscle contraction. These results demonstrate that CHIKV replicates in the hearts of immunodeficient mice and induce heart-specific mutations and host responses with implications for cardiac pathologies.

Sensitive and accurate analysis of gene expression signatures enabled by oligonucleotide-labelled cDNA

High-throughput RNA sequencing offers a comprehensive analysis of transcriptome complexity originated from regulatory events, such as differential gene expression, alternative polyadenylation and others, and allows the increase in diagnostic capacity and precision. For gene expression profiling applications that do not specifically require information on alternative splicing events, the mRNA 3' termini counting approach is a cost-effective alternative to whole transcriptome sequencing. Here, we report MTAS-seq (mRNA sequencing via terminator-assisted synthesis) - a novel RNA-seq library preparation method directed towards mRNA 3' termini. We demonstrate the specific enrichment for 3'-terminal regions by simple and quick single-tube protocol with built-in molecular barcoding to enable accurate estimation of transcript abundance. To achieve that, we synthesized oligonucleotide-modified dideoxynucleotides which enable the generation of cDNA libraries at the reverse transcription step. We validated the performance of MTAS-seq on well-characterized reference bulk RNA and further tested it with eukaryotic cell lysates.

RBFOX2 is critical for maintaining alternative polyadenylation patterns and mitochondrial health in rat myoblasts

RBFOX2, which has a well-established role in alternative splicing, is linked to heart diseases. However, it is unclear whether RBFOX2 has other roles in RNA processing that can influence gene expression in muscle cells, contributing to heart disease. Here, we employ both 3'-end and nanopore cDNA sequencing to reveal a previously unrecognized role for RBFOX2 in maintaining alternative polyadenylation (APA) signatures in myoblasts. RBFOX2-mediated APA modulates mRNA levels and/or isoform expression of a collection of genes, including contractile and mitochondrial genes. Depletion of RBFOX2 adversely affects mitochondrial health in myoblasts, correlating with disrupted APA of mitochondrial gene Slc25a4. Mechanistically, RBFOX2 regulation of Slc25a4 APA is mediated through consensus RBFOX2 binding motifs near the distal polyadenylation site, enforcing the use of the proximal polyadenylation site. In sum, our results unveil a role for RBFOX2 in fine-tuning expression of mitochondrial and contractile genes via APA in myoblasts relevant to heart diseases.

Tiled-ClickSeq for targeted sequencing of complete coronavirus genomes with simultaneous capture of RNA recombination and minority variants

High-throughput genomics of SARS-CoV-2 is essential to characterize virus evolution and to identify adaptations that affect pathogenicity or transmission. While single-nucleotide variations (SNVs) are commonly considered as driving virus adaption, RNA recombination events that delete or insert nucleic acid sequences are also critical. Whole genome targeting sequencing of SARS-CoV-2 is typically achieved using pairs of primers to generate cDNA amplicons suitable for next-generation sequencing (NGS). However, paired-primer approaches impose constraints on where primers can be designed, how many amplicons are synthesized and requires multiple PCR reactions with non-overlapping primer pools. This imparts sensitivity to underlying SNVs and fails to resolve RNA recombination junctions that are not flanked by primer pairs. To address these limitations, we have designed an approach called ‘Tiled-ClickSeq’, which uses hundreds of tiled-primers spaced evenly along the virus genome in a single reverse-transcription reaction. The other end of the cDNA amplicon is generated by azido-nucleotides that stochastically terminate cDNA synthesis, removing the need for a paired-primer. A sequencing adaptor containing a Unique Molecular Identifier (UMI) is appended to the cDNA fragment using click-chemistry and a PCR reaction generates a final NGS library. Tiled-ClickSeq provides complete genome coverage, including the 5’UTR, at high depth and specificity to the virus on both Illumina and Nanopore NGS platforms. Here, we analyze multiple SARS-CoV-2 isolates and clinical samples to simultaneously characterize minority variants, sub-genomic mRNAs (sgmRNAs), structural variants (SVs) and D-RNAs. Tiled-ClickSeq therefore provides a convenient and robust platform for SARS-CoV-2 genomics that captures the full range of RNA species in a single, simple assay.

Tiled-ClickSeq for targeted sequencing of complete coronavirus genomes with simultaneous capture of RNA recombination and minority variants

High-throughput genomics of SARS-CoV-2 is essential to characterize virus evolution and to identify adaptations that affect pathogenicity or transmission. While single-nucleotide variations (SNVs) are commonly considered as driving virus adaption, RNA recombination events that delete or insert nucleic acid sequences are also critical. Whole genome targeting sequencing of SARS-CoV-2 is typically achieved using pairs of primers to generate cDNA amplicons suitable for Next-Generation Sequencing (NGS). However, paired-primer approaches impose constraints on where primers can be designed, how many amplicons are synthesized and requires multiple PCR reactions with non-overlapping primer pools. This imparts sensitivity to underlying SNVs and fails to resolve RNA recombination junctions that are not flanked by primer pairs. To address these limitations, we have designed an approach called 'Tiled-ClickSeq', which uses hundreds of tiled-primers spaced evenly along the virus genome in a single reverse-transcription reaction. The other end of the cDNA amplicon is generated by azido-nucleotides that stochastically terminate cDNA synthesis, removing the need for a paired-primer. A sequencing adaptor containing a Unique Molecular Identifier (UMI) is appended to the cDNA fragment using click-chemistry and a PCR reaction generates a final NGS library. Tiled-ClickSeq provides complete genome coverage, including the 5'UTR, at high depth and specificity to the virus on both Illumina and Nanopore NGS platforms. Here, we analyze multiple SARS-CoV-2 isolates and clinical samples to simultaneously characterize minority variants, sub-genomic mRNAs (sgmRNAs), structural variants (SVs) and D-RNAs. Tiled-ClickSeq therefore provides a convenient and robust platform for SARS-CoV-2 genomics that captures the full range of RNA species in a single, simple assay.

QuantifyPoly(A): reshaping alternative polyadenylation landscapes of eukaryotes with weighted density peak clustering

The dynamic choice of different polyadenylation sites in a gene is referred to as alternative polyadenylation, which functions in many important biological processes. Large-scale messenger RNA 3' end sequencing has revealed that cleavage sites for polyadenylation are presented with microheterogeneity. To date, the conventional determination of polyadenylation site clusters is subjective and arbitrary, leading to inaccurate annotations. Here, we present a weighted density peak clustering method, QuantifyPoly(A), to accurately quantify genome-wide polyadenylation choices. Applying QuantifyPoly(A) on published 3' end sequencing datasets from both animals and plants, their polyadenylation profiles are reshaped into myriads of novel polyadenylation site clusters. Most of these novel polyadenylation site clusters show significantly dynamic usage across different biological samples or associate with binding sites of trans-acting factors. Upstream sequences of these clusters are enriched with polyadenylation signals UGUA, UAAA and/or AAUAAA in a species-dependent manner. Polyadenylation site clusters also exhibit species specificity, while plants ones generally show higher microheterogeneity than that of animals. QuantifyPoly(A) is broadly applicable to any types of 3' end sequencing data and species for accurate quantification and construction of the complex and dynamic polyadenylation landscape and enables us to decode alternative polyadenylation events invisible to conventional methods at a much higher resolution.

An atlas of alternative polyadenylation quantitative trait loci contributing to complex trait and disease heritability

Genome-wide association studies have identified thousands of noncoding variants associated with human traits and diseases. However, the functional interpretation of these variants is a major challenge. Here, we constructed a multi-tissue atlas of human 3'UTR alternative polyadenylation (APA) quantitative trait loci (3'aQTLs), containing approximately 0.4 million common genetic variants associated with the APA of target genes, identified in 46 tissues isolated from 467 individuals (Genotype-Tissue Expression Project). Mechanistically, 3'aQTLs can alter poly(A) motifs, RNA secondary structure and RNA-binding protein-binding sites, leading to thousands of APA changes. Our CRISPR-based experiments indicate that such 3'aQTLs can alter APA regulation. Furthermore, we demonstrate that mapping 3'aQTLs can identify APA regulators, such as La-related protein 4. Finally, 3'aQTLs are colocalized with approximately 16.1% of trait-associated variants and are largely distinct from other QTLs, such as expression QTLs. Together, our findings show that 3'aQTLs contribute substantially to the molecular mechanisms underlying human complex traits and diseases.

QPAT-seq, a rapid and deduplicatable method for quantification of poly(A) site usages

Alternative polyadenylation (APA) is an essential regulatory mechanism for gene expression. The next generation sequencing provides ample opportunity to precisely delineate APA sites genome-wide. Various methods for profiling transcriptome-wide poly(A) sites were developed. By comparing available methods, the ways for adding sequencing adaptors to fit with the Illumina sequencing platform are different. These methods have identified more than 50% genes that undergo APA in eukaryotes. However, due to the unbalanced PCR during library preparation, accurate quantification of poly(A) sites is still a challenge. Here, we describe an updated poly(A) tag sequencing method that incorporates unique molecular identifier (UMI) into the adaptor for removing quantification bias induced by PCR duplicates. Hence, quantification of poly(A) site usages can be achieved by counting UMIs. This protocol, quantifying poly(A) tag sequencing (QPAT-seq), can be finished in 1 day with reduced cost, and is particularly useful for application with a large number of samples.

The Detection and Bioinformatic Analysis of Alternative 3' UTR Isoforms as Potential Cancer Biomarkers

Alternative transcript cleavage and polyadenylation is linked to cancer cell transformation, proliferation and outcome. This has led researchers to develop methods to detect and bioinformatically analyse alternative polyadenylation as potential cancer biomarkers. If incorporated into standard prognostic measures such as gene expression and clinical parameters, these could advance cancer prognostic testing and possibly guide therapy. In this review, we focus on the existing methodologies, both experimental and computational, that have been applied to support the use of alternative polyadenylation as cancer biomarkers.

SAM homeostasis is regulated by CFIm-mediated splicing of MAT2A

S-adenosylmethionine (SAM) is the methyl donor for nearly all cellular methylation events. Cells regulate intracellular SAM levels through intron detention of MAT2A, the only SAM synthetase expressed in most cells. The N6-adenosine methyltransferase METTL16 promotes splicing of the MAT2A detained intron by an unknown mechanism. Using an unbiased CRISPR knock-out screen, we identified CFIm25 (NUDT21) as a regulator of MAT2A intron detention and intracellular SAM levels. CFIm25 is a component of the cleavage factor Im (CFIm) complex that regulates poly(A) site selection, but we show it promotes MAT2A splicing independent of poly(A) site selection. CFIm25-mediated MAT2A splicing induction requires the RS domains of its binding partners, CFIm68 and CFIm59 as well as binding sites in the detained intron and 3´ UTR. These studies uncover mechanisms that regulate MAT2A intron detention and reveal a previously undescribed role for CFIm in splicing and SAM metabolism.

Application and design considerations for 3'-end sequencing using click-chemistry

Over the past 15 years, investigations into alternative polyadenylation (APA) and its function in cellular physiology and pathology have greatly expanded due to the emergent appreciation of its key role in driving transcriptomic diversity. This growth has necessitated the development of new technologies capable of monitoring cleavage and polyadenylation events genome-wide. Advancements in approaches include both the creation of computational tools to re-analyze RNA-seq to identify APA events as well as targeted sequencing approaches customized to focus on the 3'-end of mRNA. Here we describe a streamlined protocol for polyA-Click-seq (PAC-seq), which utilizes click-chemistry to create mRNA 3'-ends sequencing libraries. Importantly, we offer additional considerations not present in our previous protocols including the use of spike-ins, unique molecular identifier primers, and guidance for appropriate depth of PAC-seq. In conjunction with the companion chapter on PolyA-miner (Yalamanchili et al., 2021) to computationally analyze PAC-seq data, we provide a complete experimental pipeline to analyze mRNA 3'-end usage in eukaryotic cells.

3' End sequencing of pA+ and pA- RNAs

The identity and metabolism of RNAs are often governed by their 5' and 3' ends. Single gene loci produce a variety of transcript isoforms, varying primarily in their RNA 3' end status and consequently facing radically different cellular fates. Knowledge about RNA termini is therefore key to understanding the diverse RNA output from individual transcription units. In addition, the 3' end of a nascent RNA at the catalytic center of RNA polymerase provides a precise and strand-specific measure of the transcription process. Here, we describe a modified RNA 3' end sequencing method, that utilizes the in vivo metabolic labeling of RNA followed by its purification and optional in vitro polyadenylation to provide a comprehensive view of all RNA 3' ends. The strategy offers the advantages of (i) nucleotide resolution mapping of RNA 3' ends, (ii) increased sequencing depth of lowly abundant RNA and (iii) inference of RNA 3' end polyadenylation status. We have used the method to study RNA decay and transcription termination mechanisms with the potential utility to a wider range of biological questions.

Human Pumilio proteins directly bind the CCR4-NOT deadenylase complex to regulate the transcriptome

Pumilio paralogs, PUM1 and PUM2, are sequence-specific RNA-binding proteins that are essential for vertebrate development and neurological functions. PUM1&2 negatively regulate gene expression by accelerating degradation of specific mRNAs. Here, we determined the repression mechanism and impact of human PUM1&2 on the transcriptome. We identified subunits of the CCR4-NOT (CNOT) deadenylase complex required for stable interaction with PUM1&2 and to elicit CNOT-dependent repression. Isoform-level RNA sequencing revealed broad coregulation of target mRNAs through the PUM-CNOT repression mechanism. Functional dissection of the domains of PUM1&2 identified a conserved amino-terminal region that confers the predominant repressive activity via direct interaction with CNOT. In addition, we show that the mRNA decapping enzyme, DCP2, has an important role in repression by PUM1&2 amino-terminal regions. Our results support a molecular model of repression by human PUM1&2 via direct recruitment of CNOT deadenylation machinery in a decapping-dependent mRNA decay pathway.

Alternative Polyadenylation in response to temperature stress contributes to gene regulation in Populus trichocarpa

BACKGROUND

Genome-wide change of polyadenylation (polyA) sites (also known as alternative polyadenylation, APA) is emerging as an important strategy of gene regulation in response to stress in plants. But little is known in woody perennials that are persistently dealing with multiple abiotic stresses.

RESULTS

Here, we performed a genome-wide profiling of polyadenylation sites under heat and cold treatments in Populus trichocarpa. Through a comprehensive analysis of polyA tail sequences, we identified 25,919 polyA-site clusters (PACs), and revealed 3429 and 3139 genes shifted polyA sites under heat and cold stresses respectively. We found that a small proportion of genes possessed APA that affected the open reading frames; and some shifts were commonly identified. Functional analysis of genes displaying shifted polyA tails suggested that pathways related to RNA metabolism were linked to regulate the APA events under both heat and cold stresses. Interestingly, we found that the heat stress induced a significantly more antisense PACs comparing to cold and control conditions. Furthermore, we showed that a unique cis-element (AAAAAA) was predominately enriched downstream of PACs in P. trichocarpa genes; and this sequence signal was only absent in shifted PACs under the heat condition, indicating a distinct APA mechanism responsive to heat tolerance.

CONCLUSIONS

This work provides a comprehensive picture of global polyadenylation patterns in response to temperatures stresses in trees. We show that the frequent change of polyA tail is a potential mechanism of gene regulation responsive to stress, which are associated with distinctive sequence signatures.

Alternative Polyadenylation: a new frontier in post transcriptional regulation

Polyadenylation of pre-messenger RNA (pre-mRNA) specific sites and termination of their downstream transcriptions are signaled by unique sequence motif structures such as AAUAAA and its auxiliary elements. Alternative polyadenylation (APA) is an important post-transcriptional regulatory mechanism that processes RNA products depending on its 3'-untranslated region (3'-UTR) specific sequence signal. APA processing can generate several mRNA isoforms from a single gene, which may have different biological functions on their target gene. As a result, cellular genomic stability, proliferation capability, and transformation feasibility could all be affected. Furthermore, APA modulation regulates disease initiation and progression. APA status could potentially act as a biomarker for disease diagnosis, severity stratification, and prognosis forecast. While the advance of modern throughout technologies, such as next generation-sequencing (NGS) and single-cell sequencing techniques, have enriched our knowledge about APA, much of APA biological process is unknown and pending for further investigation. Herein, we review the current knowledge on APA and how its regulatory complex factors (CFI/IIm, CPSF, CSTF, and RBPs) work together to determine RNA splicing location, cell cycle velocity, microRNA processing, and oncogenesis regulation. We also discuss various APA experiment strategies and the future direction of APA research.

PolyA-miner: accurate assessment of differential alternative poly-adenylation from 3'Seq data using vector projections and non-negative matrix factorization

Almost 70% of human genes undergo alternative polyadenylation (APA) and generate mRNA transcripts with varying lengths, typically of the 3′ untranslated regions (UTR). APA plays an important role in development and cellular differentiation, and its dysregulation can cause neuropsychiatric diseases and increase cancer severity. Increasing awareness of APA’s role in human health and disease has propelled the development of several 3′ sequencing (3′Seq) techniques that allow for precise identification of APA sites. However, despite the recent data explosion, there are no robust computational tools that are precisely designed to analyze 3′Seq data. Analytical approaches that have been used to analyze these data predominantly use proximal to distal usage. With about 50% of human genes having more than two APA isoforms, current methods fail to capture the entirety of APA changes and do not account for non-proximal to non-distal changes. Addressing these key challenges, this study demonstrates PolyA-miner, an algorithm to accurately detect and assess differential alternative polyadenylation specifically from 3′Seq data. Genes are abstracted as APA matrices, and differential APA usage is inferred using iterative consensus non-negative matrix factorization (NMF) based clustering. PolyA-miner accounts for all non-proximal to non-distal APA switches using vector projections and reflects precise gene-level 3′UTR changes. It can also effectively identify novel APA sites that are otherwise undetected when using reference-based approaches. Evaluation on multiple datasets—first-generation MicroArray Quality Control (MAQC) brain and Universal Human Reference (UHR) PolyA-seq data, recent glioblastoma cell line NUDT21 knockdown Poly(A)-ClickSeq (PAC-seq) data, and our own mouse hippocampal and human stem cell-derived neuron PAC-seq data—strongly supports the value and protocol-independent applicability of PolyA-miner. Strikingly, in the glioblastoma cell line data, PolyA-miner identified more than twice the number of genes with APA changes than initially reported. With the emerging importance of APA in human development and disease, PolyA-miner can significantly improve data analysis and help decode the underlying APA dynamics.

Emerging Roles of RNA 3'-end Cleavage and Polyadenylation in Pathogenesis, Diagnosis and Therapy of Human Disorders

A crucial feature of gene expression involves RNA processing to produce 3′ ends through a process termed 3′ end cleavage and polyadenylation (CPA). This ensures the nascent RNA molecule can exit the nucleus and be translated to ultimately give rise to a protein which can execute a function. Further, alternative polyadenylation (APA) can produce distinct transcript isoforms, profoundly expanding the complexity of the transcriptome. CPA is carried out by multi-component protein complexes interacting with multiple RNA motifs and is tightly coupled to transcription, other steps of RNA processing, and even epigenetic modifications. CPA and APA contribute to the maintenance of a multitude of diverse physiological processes. It is therefore not surprising that disruptions of CPA and APA can lead to devastating disorders. Here, we review potential CPA and APA mechanisms involving both loss and gain of function that can have tremendous impacts on health and disease. Ultimately we highlight the emerging diagnostic and therapeutic potential CPA and APA offer.

Partial loss of CFIm25 causes learning deficits and aberrant neuronal alternative polyadenylation

We previously showed that NUDT21-spanning copy-number variations (CNVs) are associated with intellectual disability (Gennarino et al., 2015). However, the patients' CNVs also included other genes. To determine if reduced NUDT21 function alone can cause disease, we generated Nudt21+/- mice to mimic NUDT21-deletion patients. We found that although these mice have 50% reduced Nudt21 mRNA, they only have 30% less of its cognate protein, CFIm25. Despite this partial protein-level compensation, the Nudt21+/- mice have learning deficits, cortical hyperexcitability, and misregulated alternative polyadenylation (APA) in their hippocampi. Further, to determine the mediators driving neural dysfunction in humans, we partially inhibited NUDT21 in human stem cell-derived neurons to reduce CFIm25 by 30%. This induced APA and protein level misregulation in hundreds of genes, a number of which cause intellectual disability when mutated. Altogether, these results show that disruption of NUDT21-regulated APA events in the brain can cause intellectual disability.

Alternative polyadenylation of mRNA and its role in cancer

Alternative polyadenylation (APA) is a molecular process that generates diversity at the 3′ end of RNA polymerase II transcripts from over 60% of human genes. APA is derived from the existence of multiple polyadenylation signals (PAS) within the same transcript, and results in the differential inclusion of sequence information at the 3′ end. While APA can occur between two PASs allowing for generation of transcripts with distinct coding potential from a single gene, most APA occurs within the untranslated region (3′UTR) and changes the length and content of these non-coding sequences. APA within the 3′UTR can have tremendous impact on its regulatory potential of the mRNA through a variety of mechanisms, and indeed this layer of gene expression regulation has profound impact on processes vital to cell growth and development. Recent studies have particularly highlighted the importance of APA dysregulation in cancer onset and progression. Here, we review the current knowledge of APA and its impacts on mRNA stability, translation, localization and protein localization. We also discuss the implications of APA dysregulation in cancer research and therapy.

NUDT21 inhibits bladder cancer progression through ANXA2 and LIMK2 by alternative polyadenylation

Purpose: Nudix Hydrolase 21 (NUDT21) is a crucial mediator involved in alternative polyadenylation (APA), and this molecule has been reported to be a tumor suppressor in human cancers. However, neither the role NUDT21 plays in bladder cancer (BC) nor the mechanisms which are involved have been investigated.

Methods: Expression levels of NUDT21 in BC were evaluated with real-time PCR, western blotting, and immunohistochemistry (IHC). In vitro and in vivo assays were performed to investigate the function of NUDT21 in tumorigenesis in bladder cancer cells. The TOP/FOP flash reporter assay, western blot, and global APA site profiling analysis were used to identify the pathway which mediates the biologic roles of NUDT21 in BC.

Results: NUDT21 expression is reduced in BC tissue and cells, and BC patients with lower NUDT21 expression have shorter overall and recurrent-free survival than patients with higher NUDT21 expression. NUDT21 ectopic expression or knockdown respectively profoundly inhibited or promoted the capacity of BC cells for proliferation, migration and invasion. We also identified a number of genes with shortened 3'UTRs through modulation of NUDT21 expression, and further characterized the NUDT21-regulated genes ANXA2 and LIMK2. We found NUDT21 modulates the expression of ANXA2 and LIMK2 in the Wnt/β-catenin and NF-κB signaling pathways.

Conclusions: These findings show NUDT21 plays a crucial role in BC progression, at least in part through ANXA2 and LIMK2 which act by alternative polyadenylation. NUDT21 may thus have potential as a diagnostic and therapeutic target in treatment of BC.

Alternative cleavage and polyadenylation in health and disease

Most human genes have multiple sites at which RNA 3' end cleavage and polyadenylation can occur, enabling the expression of distinct transcript isoforms under different conditions. Novel methods to sequence RNA 3' ends have generated comprehensive catalogues of polyadenylation (poly(A)) sites; their analysis using innovative computational methods has revealed how poly(A) site choice is regulated by core RNA 3' end processing factors, such as cleavage factor I and cleavage and polyadenylation specificity factor, as well as by other RNA-binding proteins, particularly splicing factors. Here, we review the experimental and computational methods that have enabled the global mapping of mRNA and of long non-coding RNA 3' ends, quantification of the resulting isoforms and the discovery of regulators of alternative cleavage and polyadenylation (APA). We highlight the different types of APA-derived isoforms and their functional differences, and illustrate how APA contributes to human diseases, including cancer and haematological, immunological and neurological diseases.

DPAC: A Tool for Differential Poly(A)-Cluster Usage from Poly(A)-Targeted RNAseq Data

Poly(A)-tail targeted RNAseq approaches, such as 3′READS, PAS-Seq and Poly(A)-ClickSeq, are becoming popular alternatives to random-primed RNAseq to focus sequencing reads just to the 3′ ends of polyadenylated RNAs to identify poly(A)-sites and characterize changes in their usage. Additionally, we and others have demonstrated that these approaches perform similarly to other RNAseq strategies for differential gene expression analysis, while saving on the volume of sequencing data required and providing a simpler library synthesis strategy. Here, we present DPAC (Differential Poly(A)-Clustering); a streamlined pipeline for the preprocessing of poly(A)-tail targeted RNAseq data, mapping of poly(A)-sites, poly(A)-site clustering and annotation, and determination of differential poly(A)-cluster usage using DESeq2. Changes in poly(A)-cluster usage is simultaneously used to report differential gene expression, differential terminal exon usage and alternative polyadenylation (APA).

Nudt21 regulates the alternative polyadenylation of Pak1 and is predictive in the prognosis of glioblastoma patients

Alternative polyadenylation (APA) has emerged as a prevalent feature associated with cancer development and progression. The advantage of APA to tumor progression is to induce oncogenes through 3'-UTR shortening, and to inactivate tumor suppressor genes via the re-routing of microRNA competition. We previously identified the Mammalian Cleavage Factor I-25 (CFIm25) (encoded by Nudt21 gene) as a master APA regulator whose expression levels directly impact the tumorigenicity of glioblastoma (GBM) in vitro and in vivo. Despite its importance, the role of Nudt21 in GBM development is not known and the genes subject to Nudt21 APA regulation that contribute to GBM progression have not been identified. Here, we find that Nudt21 is reduced in low grade glioma (LGG) and all four subtypes of high grade glioma (GBM). Reduced expression of Nudt21 associates with worse survival in TCGA LGG cohorts and two TCGA GBM cohorts. Moreover, although CFIm25 was initially identified as biochemically associated with both CFIm59 and CFIm68, we observed three CFIm distinct subcomplexes exist and CFIm59 protein level is dependent on Nudt21 expression in GBM cells, but CFIm68 is not, and that only CFIm59 predicts prognosis of GBM patients similar to Nudt21. Through the use of Poly(A)-Click-Seq to characterize APA, we define the mRNAs subject to 3'-UTR shortening upon Nudt21 depletion in GBM cells and observed enrichment in genes important in the Ras signaling pathway, including Pak1. Remarkably, we find that Pak1 expression is regulated by Nudt21 through its 3'-UTR APA, and the combination of Pak1 and Nudt21 expression generates an even stronger prognostic indicator of GBM survival versus either value used alone. Collectively, our data uncover Nudt21 and its downstream target Pak1 as a potential "combination biomarker" for predicting prognosis of GBM patients.

Development of Poly(A)-ClickSeq as a tool enabling simultaneous genome-wide poly(A)-site identification and differential expression analysis

The use of RNA-seq as a generalized tool to measure the differential expression of genes has essentially replaced the use of the microarray. Despite the acknowledged technical advantages to this approach, RNA-seq library preparation remains mostly conducted by core facilities rather than in the laboratory due to the infrastructure, expertise and time required per sample. We have recently described two 'click-chemistry' based library construction methods termed ClickSeq and Poly(A)-ClickSeq (PAC-seq) as alternatives to conventional RNA-seq that are both cost effective and rely on straightforward reagents readily available to most labs. ClickSeq is random-primed and can sequence any (unfragmented) RNA template, while PAC-seq is targeted to poly(A) tails of mRNAs. Here, we further develop PAC-seq as a platform that allows for simultaneous mapping of poly(A) sites and the measurement of differential expression of genes. We provide a detailed protocol, descriptions of appropriate computational pipelines, and a proof-of-principle dataset to illustrate the technique. PAC-seq offers a unique advantage over other 3' end mapping protocols in that it does not require additional purification, selection, or fragmentation steps allowing sample preparation directly from crude total cellular RNA. We have shown that PAC-seq is able to accurately and sensitively count transcripts for differential gene expression analysis, as well as identify alternative poly(A) sites and determine the precise nucleotides of the poly(A) tail boundaries.

Knockdown of NUDT21 inhibits proliferation and promotes apoptosis of human K562 leukemia cells through ERK pathway

Background

NUDT21 is a mammalian precursor mRNA(pre-mRNA) 3’ end processing factor and plays an important role in the selection of poly(A) sites in 3’-untranslated region (3’-UTR). NUDT21 links alternative polyadenylation with regulation of glioblastoma and osteosarcoma progression and is found to be related to drug resistance in childhood acute leukemia. However, the effect of NUDT21 on leukemia cells and the underlying mechanism are unknown.

Methods

We knocked down NUDT21 in K562 cells and applied qRT-PCR and western blotting to quantitate the mRNA and protein expression. Cell proliferating and apoptosis were investigated subsequently by flow cytometry, BrdU, Caspase3/7. RNA microarray and intracellular signaling array were used to determine the important cell signaling pathways.

Results

We clarified that the mRNA expression levels of NUDT21 are higher in primary chronic myelocytic leukemia patients and K562 leukemic cells compared with healthy controls and PBMCs. Downregulation of NUDT21 expression in K562 cells inhibits proliferation and promotes apoptosis. Screening by mRNA chip and intracellular signaling array, we found that MAPK/ERK pathway represented the main molecular mechanism underlying the effects of NUDT21 knockdown in K562 cells.

Conclusion

NUDT21 played an important role in promoting proliferation and inhibiting apoptosis in leukemia K562 cells. The underlying mechanisms involved the modulation of PTEN and a set of downstream molecules including ERK1/2.

Impact statement

The present work shows that the expression of NUDT21 was upregulated in chronic myelocytic leukemia and K562 cells. Silencing NUDT21 inhibited the proliferation and promoted the apoptosis of K562 cells. Subsequent experiments confirmed that NUDT21 promoted K562 proliferation through regulating the expression of p-ERK. Our findings may provide insights into the molecular mechanism underlying the effects of NUDT21 on leukemia cells and a novel strategy for the treatment of leukemia.

Alternative polyadenylation analysis in animals and plants: newly developed strategies for profiling, processing and validation

Alternative polyadenylation is an essential RNA processing event that contributes significantly to regulation of transcriptome diversity and functional dynamics in both animals and plants. Here we review newly developed next generation sequencing methods for genome-wide profiling of alternative polyadenylation (APA) sites, bioinformatics pipelines for data processing and both wet and dry laboratory approaches for APA validation. The library construction methods LITE-Seq (Low-Input 3'-Terminal sequencing) and PAC-seq (PolyA Click sequencing) tag polyA+ cDNA, while BAT-seq (BArcoded, three-prime specific sequencing) and PAPERCLIP (Poly(A) binding Protein-mediated mRNA 3'End Retrieval by CrossLinking ImmunoPrecipitation) enrich polyA+ RNA. Interestingly, only WTTS-seq (Whole Transcriptome Termini Site sequencing) targets both polyA+ RNA and polyA+ cDNA. Varieties of bioinformatics pipelines are well established to pursue read quality control, mapping, clustering, characterization and pathway analysis. The RHAPA (RNase H alternative polyadenylation assay) and 3'RACE-seq (3' rapid amplification of cDNA end sequencing) methods directly validate APA sites, while WTSS-seq (whole transcriptome start site sequencing), RNA-seq (RNA sequencing) and public APA databases can serve as indirect validation methods. We hope that these tools, pipelines and resources trigger huge waves of interest in the research community to investigate APA events underlying physiological, pathological and psychological changes and thus understand the information transfer events from genome to phenome relevant to economically important traits in both animals and plants.

ClickSeq: Replacing Fragmentation and Enzymatic Ligation with Click-Chemistry to Prevent Sequence Chimeras

We recently reported a fragmentation-free method for the synthesis of Next-Generation Sequencing libraries called "ClickSeq" that uses biorthogonal click-chemistry in place of enzymes for the ligation of sequencing adaptors. We found that this approach dramatically reduces artifactual chimera formation, allowing the study of rare recombination events that include viral replication intermediates and defective-interfering viral RNAs. ClickSeq illustrates how robust, bio-orthogonal chemistry can be harnessed in vitro to capture and dissect complex biological processes. Here, we describe an updated protocol for the synthesis of "ClickSeq" libraries.

An integrative model for alternative polyadenylation, IntMAP, delineates mTOR-modulated endoplasmic reticulum stress response

3'-untranslated regions (UTRs) can vary through the use of alternative polyadenylation sites during pre-mRNA processing. Multiple publically available pipelines combining high profiling technologies and bioinformatics tools have been developed to catalog changes in 3'-UTR lengths. In our recent RNA-seq experiments using cells with hyper-activated mammalian target of rapamycin (mTOR), we found that cellular mTOR activation leads to transcriptome-wide alternative polyadenylation (APA), resulting in the activation of multiple cellular pathways. Here, we developed a novel bioinformatics algorithm, IntMAP, which integrates RNA-Seq and PolyA Site (PAS)-Seq data for a comprehensive characterization of APA events. By applying IntMAP to the datasets from cells with hyper-activated mTOR, we identified novel APA events that could otherwise not be identified by either profiling method alone. Several transcription factors including Cebpg (CCAAT/enhancer binding protein gamma) were among the newly discovered APA transcripts, indicating that diverse transcriptional networks may be regulated by mTOR-coordinated APA. The prevention of APA in Cebpg using the CRISPR/cas9-mediated genome editing tool showed that mTOR-driven 3'-UTR shortening in Cebpg is critical in protecting cells from endoplasmic reticulum (ER) stress. Taken together, we present IntMAP as a new bioinformatics algorithm for APA analysis by which we expand our understanding of the physiological role of mTOR-coordinated APA events to ER stress response. IntMAP toolbox is available at http://compbio.cs.umn.edu/IntMAP/.

Untranslated Parts of Genes Interpreted: Making Heads or Tails of High-Throughput Transcriptomic Data via Computational Methods: Computational methods to discover and quantify isoforms with alternative untranslated regions

In this review we highlight the importance of defining the untranslated parts of transcripts, and present a number of computational approaches for the discovery and quantification of alternative transcription start and poly-adenylation events in high-throughput transcriptomic data. The fate of eukaryotic transcripts is closely linked to their untranslated regions, which are determined by the position at which transcription starts and ends at a genomic locus. Although the extent of alternative transcription starts and alternative poly-adenylation sites has been revealed by sequencing methods focused on the ends of transcripts, the application of these methods is not yet widely adopted by the community. We suggest that computational methods applied to standard high-throughput technologies are a useful, albeit less accurate, alternative to the expertise-demanding 5' and 3' sequencing and they are the only option for analysing legacy transcriptomic data. We review these methods here, focusing on technical challenges and arguing for the need to include better normalization of the data and more appropriate statistical models of the expected variation in the signal.