Nudt21 Controls Cell Fate by Connecting Alternative Polyadenylation to Chromatin Signaling

Multiplexed single-cell characterization of alternative polyadenylation regulators

Most mammalian genes have multiple polyA sites, representing a substantial source of transcript diversity regulated by the cleavage and polyadenylation (CPA) machinery. To better understand how these proteins govern polyA site choice, we introduce CPA-Perturb-seq, a multiplexed perturbation screen dataset of 42 CPA regulators with a 3' scRNA-seq readout that enables transcriptome-wide inference of polyA site usage. We develop a framework to detect perturbation-dependent changes in polyadenylation and characterize modules of co-regulated polyA sites. We find groups of intronic polyA sites regulated by distinct components of the nuclear RNA life cycle, including elongation, splicing, termination, and surveillance. We train and validate a deep neural network (APARENT-Perturb) for tandem polyA site usage, delineating a cis-regulatory code that predicts perturbation response and reveals interactions between regulatory complexes. Our work highlights the potential for multiplexed single-cell perturbation screens to further our understanding of post-transcriptional regulation.

The role of alternative polyadenylation in breast cancer

Breast cancer (BC), as a highly prevalent malignant tumor worldwide, is still unclear in its pathogenesis and has poor therapeutic outcomes. Alternative polyadenylation (APA) is a post-transcriptional regulatory mechanism widely found in eukaryotes. Precursor mRNA (pre-mRNA) undergoes the APA process to generate multiple mRNA isoforms with different coding regions or 3'UTRs, thereby greatly increasing the diversity and complexity of the eukaryotic transcriptome and proteome. Studies have shown that APA is involved in the progression of various diseases, including cancer, and plays a crucial role. Therefore, clarifying the biological mechanisms of APA and its regulators in breast cancer will help to comprehensively understand the pathogenesis of breast cancer and provide new ideas for its prevention and treatment.

Quantifying 3'UTR length from scRNA-seq data reveals changes independent of gene expression

Although more than half of all genes generate transcripts that differ in 3'UTR length, current analysis pipelines only quantify the amount but not the length of mRNA transcripts. 3'UTR length is determined by 3' end cleavage sites (CS). We map CS in more than 200 primary human and mouse cell types and increase CS annotations relative to the GENCODE database by 40%. Approximately half of all CS are used in few cell types, revealing that most genes only have one or two major 3' ends. We incorporate the CS annotations into a computational pipeline, called scUTRquant, for rapid, accurate, and simultaneous quantification of gene and 3'UTR isoform expression from single-cell RNA sequencing (scRNA-seq) data. When applying scUTRquant to data from 474 cell types and 2134 perturbations, we discover extensive 3'UTR length changes across cell types that are as widespread and coordinately regulated as gene expression changes but affect mostly different genes. Our data indicate that mRNA abundance and mRNA length are two largely independent axes of gene regulation that together determine the amount and spatial organization of protein synthesis.

InPACT: a computational method for accurate characterization of intronic polyadenylation from RNA sequencing data

Alternative polyadenylation can occur in introns, termed intronic polyadenylation (IPA), has been implicated in diverse biological processes and diseases, as it can produce noncoding transcripts or transcripts with truncated coding regions. However, a reliable method is required to accurately characterize IPA. Here, we propose a computational method called InPACT, which allows for the precise characterization of IPA from conventional RNA-seq data. InPACT successfully identifies numerous previously unannotated IPA transcripts in human cells, many of which are translated, as evidenced by ribosome profiling data. We have demonstrated that InPACT outperforms other methods in terms of IPA identification and quantification. Moreover, InPACT applied to monocyte activation reveals temporally coordinated IPA events. Further application on single-cell RNA-seq data of human fetal bone marrow reveals the expression of several IPA isoforms in a context-specific manner. Therefore, InPACT represents a powerful tool for the accurate characterization of IPA from RNA-seq data.

Downregulation of CPSF6 leads to global mRNA 3' UTR shortening and enhanced antiviral immune responses

Alternative polyadenylation (APA) is a widespread mechanism of gene regulation that generates mRNA isoforms with alternative 3' untranslated regions (3' UTRs). Our previous study has revealed the global 3' UTR shortening of host mRNAs through APA upon viral infection. However, how the dynamic changes in the APA landscape occur upon viral infection remains largely unknown. Here we further found that, the reduced protein abundance of CPSF6, one of the core 3' processing factors, promotes the usage of proximal poly(A) sites (pPASs) of many immune related genes in macrophages and fibroblasts upon viral infection. Shortening of the 3' UTR of these transcripts may improve their mRNA stability and translation efficiency, leading to the promotion of type I IFN (IFN-I) signalling-based antiviral immune responses. In addition, dysregulated expression of CPSF6 is also observed in many immune related physiological and pathological conditions, especially in various infections and cancers. Thus, the global APA dynamics of immune genes regulated by CPSF6, can fine-tune the antiviral response as well as the responses to other cellular stresses to maintain the tissue homeostasis, which may represent a novel regulatory mechanism for antiviral immunity.

Lysine 68 Methylation-Dependent SOX9 Stability Control Modulates Chondrogenic Differentiation in Dental Pulp Stem Cells

Dental pulp stem cells (DPSCs), characterized by easy availability, multi-lineage differentiation ability, and high proliferation ability, are ideal seed cells for cartilage tissue engineering. However, the epigenetic mechanism underlying chondrogenesis in DPSCs remains elusive. Herein, it is demonstrated that KDM3A and G9A, an antagonistic pair of histone-modifying enzymes, bidirectionally regulate the chondrogenic differentiation of DPSCs by controlling SOX9 (sex-determining region Y-type high-mobility group box protein 9) degradation through lysine methylation. Transcriptomics analysis reveals that KDM3A is significantly upregulated during the chondrogenic differentiation of DPSCs. In vitro and in vivo functional analyses further indicate that KDM3A promotes chondrogenesis in DPSCs by boosting the SOX9 protein level, while G9A hinders the chondrogenic differentiation of DPSCs by reducing the SOX9 protein level. Furthermore, mechanistic studies indicate that KDM3A attenuates the ubiquitination of SOX9 by demethylating lysine (K) 68 residue, which in turn enhances SOX9 stability. Reciprocally, G9A facilitates SOX9 degradation by methylating K68 residue to increase the ubiquitination of SOX9. Meanwhile, BIX-01294 as a highly specific G9A inhibitor significantly induces the chondrogenic differentiation of DPSCs. These findings provide a theoretical basis to ameliorate the clinical use of DPSCs in cartilage tissue-engineering therapies.

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.

Transcriptional buffering and 3'UTR lengthening are shaped during human neurodevelopment by shifts in mRNA stability and microRNA load

The contribution of mRNA half-life is commonly overlooked when examining changes in mRNA abundance during development. mRNA levels of some genes are regulated by transcription rate only, but others may be regulated by mRNA half-life only shifts. Furthermore, transcriptional buffering is predicted when changes in transcription rates have compensating shifts in mRNA half-life resulting in no change to steady-state levels. Likewise, transcriptional boosting should result when changes in transcription rate are accompanied by amplifying half-life shifts. During neurodevelopment there is widespread 3'UTR lengthening that could be shaped by differential shifts in the stability of existing short or long 3'UTR transcript isoforms. We measured transcription rate and mRNA half-life changes during induced human Pluripotent Stem Cell (iPSC)-derived neuronal development using RATE-seq. During transitions to progenitor and neuron stages, transcriptional buffering occurred in up to 50%, and transcriptional boosting in up to 15%, of genes with changed transcription rates. The remaining changes occurred by transcription rate only or mRNA half-life only shifts. Average mRNA half-life decreased two-fold in neurons relative to iPSCs. Short gene isoforms were more destabilized in neurons and thereby increased the average 3'UTR length. Small RNA sequencing captured an increase in microRNA copy number per cell during neurodevelopment. We propose that mRNA destabilization and 3'UTR lengthening are driven in part by an increase in microRNA load in neurons. Our findings identify mRNA stability mechanisms in human neurodevelopment that regulate gene and isoform level abundance and provide a precedent for similar post-transcriptional regulatory events as other tissues develop.

CPA-Perturb-seq: Multiplexed single-cell characterization of alternative polyadenylation regulators

Most mammalian genes have multiple polyA sites, representing a substantial source of transcript diversity that is governed by the cleavage and polyadenylation (CPA) regulatory machinery. To better understand how these proteins govern polyA site choice we introduce CPA-Perturb-seq, a multiplexed perturbation screen dataset of 42 known CPA regulators with a 3' scRNA-seq readout that enables transcriptome-wide inference of polyA site usage. We develop a statistical framework to specifically identify perturbation-dependent changes in intronic and tandem polyadenylation, and discover modules of co-regulated polyA sites exhibiting distinct functional properties. By training a multi-task deep neural network (APARENT-Perturb) on our dataset, we delineate a cis-regulatory code that predicts responsiveness to perturbation and reveals interactions between distinct regulatory complexes. Finally, we leverage our framework to re-analyze published scRNA-seq datasets, identifying new regulators that affect the relative abundance of alternatively polyadenylated transcripts, and characterizing extensive cellular heterogeneity in 3' UTR length amongst antibody-producing cells. Our work highlights the potential for multiplexed single-cell perturbation screens to further our understanding of post-transcriptional regulation in vitro and in vivo.

PolyAMiner-Bulk: A Machine Learning Based Bioinformatics Algorithm to Infer and Decode Alternative Polyadenylation Dynamics from bulk RNA-seq data

More than half of human genes exercise alternative polyadenylation (APA) and generate mRNA transcripts with varying 3' untranslated regions (UTR). However, current computational approaches for identifying cleavage and polyadenylation sites (C/PASs) and quantifying 3'UTR length changes from bulk RNA-seq data fail to unravel tissue- and disease-specific APA dynamics. Here, we developed a next-generation bioinformatics algorithm and application, PolyAMiner-Bulk, that utilizes an attention-based machine learning architecture and an improved vector projection-based engine to infer differential APA dynamics accurately. When applied to earlier studies, PolyAMiner-Bulk accurately identified more than twice the number of APA changes in an RBM17 knockdown bulk RNA-seq dataset compared to current generation tools. Moreover, on a separate dataset, PolyAMiner-Bulk revealed novel APA dynamics and pathways in scleroderma pathology and identified differential APA in a gene that was identified as being involved in scleroderma pathogenesis in an independent study. Lastly, we used PolyAMiner-Bulk to analyze the RNA-seq data of post-mortem prefrontal cortexes from the ROSMAP data consortium and unraveled novel APA dynamics in Alzheimer's Disease. Our method, PolyAMiner-Bulk, creates a paradigm for future alternative polyadenylation analysis from bulk RNA-seq data.

CFIm-mediated alternative polyadenylation safeguards the development of mammalian pre-implantation embryos

Alternative polyadenylation (APA) gives rise to transcripts with distinct 3' untranslated regions (3' UTRs), thereby affecting the fate of mRNAs. APA is strongly associated with cell proliferation and differentiation status, and thus likely plays a critical role in the embryo development. However, the pattern of APA in mammalian early embryos is still unknown. Here, we analyzed the 3' UTR lengths in human and mouse pre-implantation embryos using available single cell RNA-seq datasets and explored the underlying mechanism driving the changes. Although human and mouse early embryos displayed distinct patterns of 3' UTR changing, RNA metabolism pathways were involved in both species. The 3' UTR lengths are likely determined by the abundance of the cleavage factor I complex (CFIm) components NUDT21 and CPSF6 in the nucleus. Importantly, depletion of either component resulted in early embryo development arrest and 3' UTR shortening. Collectively, these data highlight an essential role for APA in the development of mammalian early embryos.

Context-specific regulation and function of mRNA alternative polyadenylation

Alternative cleavage and polyadenylation (APA) is a widespread mechanism to generate mRNA isoforms with alternative 3' untranslated regions (UTRs). The expression of alternative 3' UTR isoforms is highly cell type specific and is further controlled in a gene-specific manner by environmental cues. In this Review, we discuss how the dynamic, fine-grained regulation of APA is accomplished by several mechanisms, including cis-regulatory elements in RNA and DNA and factors that control transcription, pre-mRNA cleavage and post-transcriptional processes. Furthermore, signalling pathways modulate the activity of these factors and integrate APA into gene regulatory programmes. Dysregulation of APA can reprogramme the outcome of signalling pathways and thus can control cellular responses to environmental changes. In addition to the regulation of protein abundance, APA has emerged as a major regulator of mRNA localization and the spatial organization of protein synthesis. This role enables the regulation of protein function through the addition of post-translational modifications or the formation of protein-protein interactions. We further discuss recent transformative advances in single-cell RNA sequencing and CRISPR-Cas technologies, which enable the mapping and functional characterization of alternative 3' UTRs in any biological context. Finally, we discuss new APA-based RNA therapeutics, including compounds that target APA in cancer and therapeutic genome editing of degenerative diseases.

NUDT21 limits CD19 levels through alternative mRNA polyadenylation in B cell acute lymphoblastic leukemia

B cell progenitor acute lymphoblastic leukemia (B-ALL) treatment has been revolutionized by T cell-based immunotherapies-including chimeric antigen receptor T cell therapy (CAR-T) and the bispecific T cell engager therapeutic, blinatumomab-targeting surface glycoprotein CD19. Unfortunately, many patients with B-ALL will fail immunotherapy due to 'antigen escape'-the loss or absence of leukemic CD19 targeted by anti-leukemic T cells. In the present study, we utilized a genome-wide CRISPR-Cas9 screening approach to identify modulators of CD19 abundance on human B-ALL blasts. These studies identified a critical role for the transcriptional activator ZNF143 in CD19 promoter activation. Conversely, the RNA-binding protein, NUDT21, limited expression of CD19 by regulating CD19 messenger RNA polyadenylation and stability. NUDT21 deletion in B-ALL cells increased the expression of CD19 and the sensitivity to CD19-specific CAR-T and blinatumomab. In human B-ALL patients treated with CAR-T and blinatumomab, upregulation of NUDT21 mRNA coincided with CD19 loss at disease relapse. Together, these studies identify new CD19 modulators in human B-ALL.

Diverse cell-specific patterns of alternative polyadenylation in Drosophila

Most genes in higher eukaryotes express isoforms with distinct 3' untranslated regions (3' UTRs), generated by alternative polyadenylation (APA). Since 3' UTRs are predominant locations of post-transcriptional regulation, APA can render such programs conditional, and can also alter protein sequences via alternative last exon (ALE) isoforms. We previously used 3'-sequencing from diverse Drosophila samples to define multiple tissue-specific APA landscapes. Here, we exploit comprehensive single nucleus RNA-sequencing data (Fly Cell Atlas) to elucidate cell-type expression of 3' UTRs across >250 adult Drosophila cell types. We reveal the cellular bases of multiple tissue-specific APA/ALE programs, such as 3' UTR lengthening in differentiated neurons and 3' UTR shortening in spermatocytes and spermatids. We trace dynamic 3' UTR patterns across cell lineages, including in the male germline, and discover new APA patterns in the intestinal stem cell lineage. Finally, we correlate expression of RNA binding proteins (RBPs), miRNAs and global levels of cleavage and polyadenylation (CPA) factors in several cell types that exhibit characteristic APA landscapes, yielding candidate regulators of transcriptome complexity. These analyses provide a comprehensive foundation for future investigations of mechanisms and biological impacts of alternative 3' isoforms across the major cell types of this widely-studied model organism.

Enhancers regulate 3' end processing activity to control expression of alternative 3'UTR isoforms

Multi-UTR genes are widely transcribed and express their alternative 3'UTR isoforms in a cell type-specific manner. As transcriptional enhancers regulate mRNA expression, we investigated if they also regulate 3'UTR isoform expression. Endogenous enhancer deletion of the multi-UTR gene PTEN did not impair transcript production but prevented 3'UTR isoform switching which was recapitulated by silencing of an enhancer-bound transcription factor. In reporter assays, enhancers increase transcript production when paired with single-UTR gene promoters. However, when combined with multi-UTR gene promoters, they change 3'UTR isoform expression by increasing 3' end processing activity of polyadenylation sites. Processing activity of polyadenylation sites is affected by transcription factors, including NF-κB and MYC, transcription elongation factors, chromatin remodelers, and histone acetyltransferases. As endogenous cell type-specific enhancers are associated with genes that increase their short 3'UTRs in a cell type-specific manner, our data suggest that transcriptional enhancers integrate cellular signals to regulate cell type-and condition-specific 3'UTR isoform expression.

CFIm-mediated alternative polyadenylation remodels cellular signaling and miRNA biogenesis

The mammalian cleavage factor I (CFIm) has been implicated in alternative polyadenylation (APA) in a broad range of contexts, from cancers to learning deficits and parasite infections. To determine how the CFIm expression levels are translated into these diverse phenotypes, we carried out a multi-omics analysis of cell lines in which the CFIm25 (NUDT21) or CFIm68 (CPSF6) subunits were either repressed by siRNA-mediated knockdown or over-expressed from stably integrated constructs. We established that >800 genes undergo coherent APA in response to changes in CFIm levels, and they cluster in distinct functional classes related to protein metabolism. The activity of the ERK pathway traces the CFIm concentration, and explains some of the fluctuations in cell growth and metabolism that are observed upon CFIm perturbations. Furthermore, multiple transcripts encoding proteins from the miRNA pathway are targets of CFIm-dependent APA. This leads to an increased biogenesis and repressive activity of miRNAs at the same time as some 3′ UTRs become shorter and presumably less sensitive to miRNA-mediated repression. Our study provides a first systematic assessment of a core set of APA targets that respond coherently to changes in CFIm protein subunit levels (CFIm25/CFIm68). We describe the elicited signaling pathways downstream of CFIm, which improve our understanding of the key role of CFIm in integrating RNA processing with other cellular activities.

Low BACH2 Expression Predicts Adverse Outcome in Chronic Lymphocytic Leukaemia

Chronic lymphocytic leukaemia (CLL) is a heterogeneous disease with a highly variable clinical outcome. There are well-established CLL prognostic biomarkers that have transformed treatment and improved the understanding of CLL biology. Here, we have studied the clinical significance of two crucial B cell regulators, BACH2 (BTB and CNC homology 1, basic leucine zipper transcription factor 2) and BCL6 (B-cell CLL/lymphoma 6), in a cohort of 102 CLL patients and determined the protein interaction networks that they participate in using MEC-1 CLL cells. We observed that CLL patients expressing low levels of BCL6 and BACH2 RNA had significantly shorter overall survival (OS) than high BCL6- and BACH2-expressing cases. Notably, their low expression specifically decreased the OS of immunoglobulin heavy chain variable region-mutated (IGHV-M) CLL patients, as well as those with 11q and 13q deletions. Similar to the RNA data, a low BACH2 protein expression was associated with a significantly shorter OS than a high expression. There was no direct interaction observed between BACH2 and BCL6 in MEC-1 CLL cells, but they shared protein networks that included fifty different proteins. Interestingly, a prognostic index (PI) model that we generated, using integrative risk score values of BACH2 RNA expression, age, and 17p deletion status, predicted patient outcomes in our cohort. Taken together, these data have shown for the first time a possible prognostic role for BACH2 in CLL and have revealed protein interaction networks shared by BCL6 and BACH2, indicating a significant role for BACH2 and BCL6 in key cellular processes, including ubiquitination mediated B-cell receptor functions, nucleic acid metabolism, protein degradation, and homeostasis in CLL biology.

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.

PAS-seq 2: A fast and sensitive method for global profiling of polyadenylated RNAs

Alternative polyadenylation (APA) is a widespread phenomenon in eukaryotes that contributes to regulating gene expression and generating proteomic diversity. APA plays critical roles in development and its mis-regulation has been implicated in a wide variety of human diseases, including cancer. To study APA on the transcriptome-wide level, numerous deep sequencing methods that capture 3' end of mRNAs have been developed in the past decade, but they generally require a large amount of hands-on time and/or high RNA input. Here, we introduce PAS-seq 2, a fast and sensitive method for global and quantitative profiling of polyadenylated RNAs. Compared to our original PAS-seq, this method takes less time and requires much lower total RNA input due to improvement in the reverse transcription process. PAS-seq 2 can be applied to both APA and differential gene expression analyses.

NUDT21 Links Mitochondrial IPS-1 to RLR-Containing Stress Granules and Activates Host Antiviral Defense

Viral RNA in the cytoplasm of mammalian host cells is recognized by retinoic acid-inducible protein-I-like receptors (RLRs), which localize to cytoplasmic stress granules (SGs). Activated RLRs associate with the mitochondrial adaptor protein IPS-1, which activates antiviral host defense mechanisms, including type I IFN induction. It has remained unclear, however, how RLRs in SGs and IPS-1 in the mitochondrial outer membrane associate physically and engage in information transfer. In this study, we show that NUDT21, an RNA-binding protein that regulates alternative transcript polyadenylation, physically associates with IPS-1 and mediates its localization to SGs in response to transfection with polyinosinic-polycytidylic acid [poly(I:C)], a mimic of viral dsRNA. We found that despite its well-established function in the nucleus, a fraction of NUDT21 localizes to mitochondria in resting cells and becomes localized to SGs in response to poly(I:C) transfection. NUDT21 was also found to be required for efficient type I IFN induction in response to viral infection in both human HeLa cells and mouse macrophage cell line RAW264.7 cells. Our results together indicate that NUDT21 links RLRs in SGs to mitochondrial IPS-1 and thereby activates host defense responses to viral infection.

NUDT21 knockdown inhibits proliferation and promotes apoptosis of pancreatic ductal adenocarcinoma through EIF2 signaling

The NUDT family is thought to play an important role in cancer growth and progression. However, the clinicopathologic significance and potential role of nucleotide diphosphate-linked X-component motif 21, NM_007006 (NUDT21) in pancreatic ductal adenocarcinoma (PDAC) remains largely unknown. In this study, we observed that NUDT21 was frequently up-expressed in PDAC. Clinical data revealed that its level positively correlated with poor survival of patients with PDAC. We found that knockdown of NUDT21 significantly inhibited cell proliferation and promoted apoptosis both in vitro and in vivo. Screening by microarray analysis and verifying by Western blot, we found that the EIF2 signaling pathway represented the main molecular mechanism underlying the effects of NUDT21 knockdown in PANC-1 cells, and PKR, HSPA5, EIF4E and DDIT3 may be its target genes. Thus, our results revealed for the first time that NUDT21, a valuable marker of PDAC prognosis, promotes tumor proliferation, inhibits cells apoptosis and might represent a potential target for gene-based therapy.

ELAV and FNE Determine Neuronal Transcript Signatures through EXon-Activated Rescue

The production of alternative RNA variants contributes to the tissue-specific regulation of gene expression. In the animal nervous system, a systematic shift toward distal sites of transcription termination produces transcript signatures that are crucial for neuron development and function. Here, we report that, in Drosophila, the highly conserved protein ELAV globally regulates all sites of neuronal 3' end processing and directly binds to proximal polyadenylation sites of target mRNAs in vivo. We uncover an endogenous strategy of functional gene rescue that safeguards neuronal RNA signatures in an ELAV loss-of-function context. When not directly repressed by ELAV, the transcript encoding the ELAV paralog FNE acquires a mini-exon, generating a new protein able to translocate to the nucleus and rescue ELAV-mediated alternative polyadenylation and alternative splicing. We propose that exon-activated functional rescue is a more widespread mechanism that ensures robustness of processes regulated by a hierarchy, rather than redundancy, of effectors.

Downregulation of CFIm25 amplifies dermal fibrosis through alternative polyadenylation

This study implicates the key regulator of alternative polyadenylation, CFIm25 in dermal fibrosis and in systemic sclerosis (scleroderma) pathogenesis. CFIm25 downregulation promotes the expression of profibrotic factors, exaggerates bleomycin-induced skin fibrosis, while CFIm25 restoration attenuates skin fibrosis.

Systemic sclerosis (SSc; scleroderma) is a multisystem fibrotic disease. The mammalian cleavage factor I 25-kD subunit (CFIm25; encoded by NUDT21) is a key regulator of alternative polyadenylation, and its depletion causes predominantly 3′UTR shortening through loss of stimulation of distal polyadenylation sites. A shortened 3′UTR will often lack microRNA target sites, resulting in increased mRNA translation due to evasion of microRNA-mediated repression. Herein, we report that CFlm25 is downregulated in SSc skin, primary dermal fibroblasts, and two murine models of dermal fibrosis. Knockdown of CFIm25 in normal skin fibroblasts is sufficient to promote the 3′UTR shortening of key TGFβ-regulated fibrotic genes and enhance their protein expression. Moreover, several of these fibrotic transcripts show 3′UTR shortening in SSc skin. Finally, mice with CFIm25 deletion in fibroblasts show exaggerated skin fibrosis upon bleomycin treatment, and CFIm25 restoration attenuates bleomycin-induced skin fibrosis. Overall, our data link this novel RNA-processing mechanism to dermal fibrosis and SSc pathogenesis.

Regulation of alternative polyadenylation in the yeast Saccharomyces cerevisiae by histone H3K4 and H3K36 methyltransferases

Adjusting DNA structure via epigenetic modifications, and altering polyadenylation (pA) sites at which precursor mRNA is cleaved and polyadenylated, allows cells to quickly respond to environmental stress. Since polyadenylation occurs co-transcriptionally, and specific patterns of nucleosome positioning and chromatin modifications correlate with pA site usage, epigenetic factors potentially affect alternative polyadenylation (APA). We report that the histone H3K4 methyltransferase Set1, and the histone H3K36 methyltransferase Set2, control choice of pA site in Saccharomyces cerevisiae, a powerful model for studying evolutionarily conserved eukaryotic processes. Deletion of SET1 or SET2 causes an increase in serine-2 phosphorylation within the C-terminal domain of RNA polymerase II (RNAP II) and in the recruitment of the cleavage/polyadenylation complex, both of which could cause the observed switch in pA site usage. Chemical inhibition of TOR signaling, which causes nutritional stress, results in Set1- and Set2-dependent APA. In addition, Set1 and Set2 decrease efficiency of using single pA sites, and control nucleosome occupancy around pA sites. Overall, our study suggests that the methyltransferases Set1 and Set2 regulate APA induced by nutritional stress, affect the RNAP II C-terminal domain phosphorylation at Ser2, and control recruitment of the 3′ end processing machinery to the vicinity of pA sites.

Stem cell programs in cancer initiation, progression, and therapy resistance

Over the past few decades, substantial evidence has convincingly revealed the existence of cancer stem cells (CSCs) as a minor subpopulation in cancers, contributing to an aberrantly high degree of cellular heterogeneity within the tumor. CSCs are functionally defined by their abilities of self-renewal and differentiation, often in response to cues from their microenvironment. Biological phenotypes of CSCs are regulated by the integrated transcriptional, post-transcriptional, metabolic, and epigenetic regulatory networks. CSCs contribute to tumor progression, therapeutic resistance, and disease recurrence through their sustained proliferation, invasion into normal tissue, promotion of angiogenesis, evasion of the immune system, and resistance to conventional anticancer therapies. Therefore, elucidation of the molecular mechanisms that drive cancer stem cell maintenance, plasticity, and therapeutic resistance will enhance our ability to improve the effectiveness of targeted therapies for CSCs. In this review, we highlight the key features and mechanisms that regulate CSC function in tumor initiation, progression, and therapy resistance. We discuss factors for CSC therapeutic resistance, such as quiescence, induction of epithelial-to-mesenchymal transition (EMT), and resistance to DNA damage-induced cell death. We evaluate therapeutic approaches for eliminating therapy-resistant CSC subpopulations, including anticancer drugs that target key CSC signaling pathways and cell surface markers, viral therapies, the awakening of quiescent CSCs, and immunotherapy. We also assess the impact of new technologies, such as single-cell sequencing and CRISPR-Cas9 screening, on the investigation of the biological properties of CSCs. Moreover, challenges remain to be addressed in the coming years, including experimental approaches for investigating CSCs and obstacles in therapeutic targeting of CSCs.

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.

NUDT21 suppresses the growth of small cell lung cancer by modulating GLS1 splicing

The mRNA precursor 3'-end modification factor NUDT21 is a major regulator of 3'UTR shortening and an important component of pre-mRNA cleavage and polyadenylation. However, its role in pathologic progress of small cell lung cancer (SCLC) remains unclear. In this study, we observed that NUDT21 expression is downregulated in SCLC tissues. Hypoxia-induced down-regulation of NUDT21 through HIF-1α. NUDT21 shRNA transduction promotes proliferation and inhibits apoptosis of A549 cells. NUDT21 inhibition also promotes tumor growth in a mouse xenograft model. Furthermore, we clarified that HIF-1α mediated NUDT21 downregulation which altered the expression patterns of two isoforms of GLS1, GAC and KGA. These results link the hypoxic tumor environments to aberrant glutamine metabolism which is important for cellular energy in SCLC cells. Therefore, NUDT21 could be considered as a potential target for the treatment of SCLC.

3' UTRs Regulate Protein Functions by Providing a Nurturing Niche during Protein Synthesis

Messenger RNAs (mRNAs) are the templates for protein synthesis as the coding region is translated into the amino acid sequence. mRNAs also contain 3' untranslated regions (3' UTRs) that harbor additional elements for the regulation of protein function. If the amino acid sequence of a protein is necessary and sufficient for its function, we call it 3' UTR-independent. In contrast, functions that are accomplished by protein complexes whose formation requires the presence of a specific 3' UTR are 3' UTR-dependent protein functions. We showed that 3' UTRs can regulate protein activity without affecting protein abundance, and alternative 3' UTRs can diversify protein functions. We currently think that the regulation of protein function by 3' UTRs is facilitated by the local environment at the site of protein synthesis, which we call the nurturing niche for nascent proteins. This niche is composed of the mRNA and the bound proteins that consist of RNA-binding proteins and recruited proteins. It enables the formation of specific protein complexes, as was shown for TIS granules, a recently discovered cytoplasmic membraneless organelle. This finding suggests that changing the niche for nascent proteins will alter protein activity and function, implying that cytoplasmic membraneless organelles can regulate protein function in a manner that is independent of protein abundance.

What Are 3' UTRs Doing?

3' untranslated regions (3' UTRs) of messenger RNAs (mRNAs) are best known to regulate mRNA-based processes, such as mRNA localization, mRNA stability, and translation. In addition, 3' UTRs can establish 3' UTR-mediated protein-protein interactions (PPIs), and thus can transmit genetic information encoded in 3' UTRs to proteins. This function has been shown to regulate diverse protein features, including protein complex formation or posttranslational modifications, but is also expected to alter protein conformations. Therefore, 3' UTR-mediated information transfer can regulate protein features that are not encoded in the amino acid sequence. This review summarizes both 3' UTR functions-the regulation of mRNA and protein-based processes-and highlights how each 3' UTR function was discovered with a focus on experimental approaches used and the concepts that were learned. This review also discusses novel approaches to study 3' UTR functions in the future by taking advantage of recent advances in technology.

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.

Gain of Additional BIRC3 Protein Functions through 3'-UTR-Mediated Protein Complex Formation

Alternative 3' untranslated regions (3' UTRs) are widespread, but their functional roles are largely unknown. We investigated the function of the long BIRC3 3' UTR, which is upregulated in leukemia. The 3' UTR does not regulate BIRC3 protein localization or abundance but is required for CXCR4-mediated B cell migration. We established an experimental pipeline to study the mechanism of regulation and used mass spectrometry to identify BIRC3 protein interactors. In addition to 3'-UTR-independent interactors involved in known BIRC3 functions, we detected interactors that bind only to BIRC3 protein encoded from the mRNA with the long 3' UTR. They regulate several functions, including CXCR4 trafficking. We further identified RNA-binding proteins differentially bound to the alternative 3' UTRs and found that cooperative binding of Staufen and HuR mediates 3'-UTR-dependent complex formation. We show that the long 3' UTR is required for the formation of specific protein complexes that enable additional functions of BIRC3 protein beyond its 3'-UTR-independent functions.

Single cell RNA-seq identifies the origins of heterogeneity in efficient cell transdifferentiation and reprogramming

Forced transcription factor expression can transdifferentiate somatic cells into other specialised cell types or reprogram them into induced pluripotent stem cells (iPSCs) with variable efficiency. To better understand the heterogeneity of these processes, we used single-cell RNA sequencing to follow the transdifferentation of murine pre-B cells into macrophages as well as their reprogramming into iPSCs. Even in these highly efficient systems, there was substantial variation in the speed and path of fate conversion. We predicted and validated that these differences are inversely coupled and arise in the starting cell population, with Myc^high large pre-BII cells transdifferentiating slowly but reprogramming efficiently and Myc^low small pre-BII cells transdifferentiating rapidly but failing to reprogram. Strikingly, differences in Myc activity predict the efficiency of reprogramming across a wide range of somatic cell types. These results illustrate how single cell expression and computational analyses can identify the origins of heterogeneity in cell fate conversion processes.

Maternal Proteins That Are Phosphoregulated upon Egg Activation Include Crucial Factors for Oogenesis, Egg Activation and Embryogenesis in Drosophila melanogaster

Egg activation is essential for the successful transition from a mature oocyte to a developmentally competent egg. It consists of a series of events including the resumption and completion of meiosis, initiation of translation of some maternal mRNAs and destruction of others, and changes to the vitelline envelope. This major change of cell state is accompanied by large scale alteration in the oocyte’s phosphoproteome. We hypothesize that the cohort of proteins that are subject to phosphoregulation during egg activation are functionally important for processes before, during, or soon after this transition, potentially uniquely or as proteins carrying out essential cellular functions like those they do in other (somatic) cells. In this study, we used germline-specific RNAi to examine the function of 189 maternal proteins that are phosphoregulated during egg activation in Drosophila melanogaster. We identified 53 genes whose knockdown reduced or abolished egg production and caused a range of defects in ovarian morphology, as well as 51 genes whose knockdown led to significant impairment or abolishment of the egg hatchability. We observed different stages of developmental arrest in the embryos and various defects in spindle morphology and aberrant centrosome activities in the early arrested embryos. Our results, validated by the detection of multiple genes with previously-documented maternal effect phenotypes among the proteins we tested, revealed 15 genes with newly discovered roles in egg activation and early embryogenesis in Drosophila. Given that protein phosphoregulation is a conserved characteristic of this developmental transition, we suggest that the phosphoregulated proteins may provide a rich pool of candidates for the identification of important players in the egg-to-embryo transition.