Dynamics of alternative polyadenylation in single root cells of Arabidopsis thaliana.

Introduction: Single-cell RNA-seq (scRNA-seq) technologies have been widely used to reveal the diversity and complexity of cells, and pioneering studies on scRNA-seq in plants began to emerge since 2019. However, existing studies on plants utilized scRNA-seq focused only on the gene expression regulation. As an essential post-transcriptional mechanism for regulating gene expression, alternative polyadenylation (APA) generates diverse mRNA isoforms with distinct 3' ends through the selective use of different polyadenylation sites in a gene. APA plays important roles in regulating multiple developmental processes in plants, such as flowering time and stress response. Methods: In this study, we developed a pipeline to identify and integrate APA sites from different scRNA-seq data and analyze APA dynamics in single cells. First, high-confidence poly(A) sites in single root cells were identified and quantified. Second, three kinds of APA markers were identified for exploring APA dynamics in single cells, including differentially expressed poly(A) sites based on APA site expression, APA markers based on APA usages, and APA switching genes based on 3' UTR (untranslated region) length change. Moreover, cell type annotations of single root cells were refined by integrating both the APA information and the gene expression profile. Results: We comprehensively compiled a single-cell APA atlas from five scRNA-seq studies, covering over 150,000 cells spanning four major tissue branches, twelve cell types, and three developmental stages. Moreover, we quantified the dynamic APA usages in single cells and identified APA markers across tissues and cell types. Further, we integrated complementary information of gene expression and APA profiles to annotate cell types and reveal subtle differences between cell types. Discussion: This study reveals that APA provides an additional layer of information for determining cell identity and provides a landscape of APA dynamics during Arabidopsis root development.

Fateful Decisions of Where to Cut the Line: Pathology Associated with Aberrant 3' End Processing and Transcription Termination.

Aberrant gene expression lies at the heart of many pathologies. This review will point out how 3' end processing, the final mRNA-maturation step in the transcription cycle, is surprisingly prone to regulated as well as stochastic variations with a wide range of consequences. Whereas smaller variations contribute to the plasticity of gene expression, larger alternations to 3' end processing and coupled transcription termination can lead to pathological consequences. These can be caused by the local mutation of one gene or affect larger numbers of genes systematically, if aspects of the mechanisms of 3' end processing and transcription termination are altered.

WDR33 alternative polyadenylation is dependent on stochastic poly(a) site usage and splicing efficiencies.

Transcripts from the human WDR33 gene, which encodes a central component of the mRNA polyadenylation (PA) machinery, are subject to alternative polyadenylation (APA) within promoter-proximal introns/exons. This APA, which itself involves usage of multiple PA sites, results in the production of two non-canonical protein isoforms, V2 and V3, that are functionally completely unrelated to the full-length protein, with roles in innate immunity. The mechanism and regulation of WDR33 APA are unclear. Here, we report that levels of the PA factor CFIm25 modulate V2 and V3 expression, and that PA site usage of both V2 and V3 varies in distinct immune responses. Using newly developed assays to measure splicing and PA site strength, we show that splicing of V2-associated intron 6 is inefficient, allowing V2 to be produced using weak PA sites. Usage of V3's strong PA sites, on the other hand, is relatively low, reflecting the high efficiency of intron 7 splicing coupled with dependency on usage of an alternative 3' splice site within the intron. Overall, our findings demonstrate that usage of WDR33 alternative PA sites is stochastic, dependent on a complex interplay between splicing and PA, and thus provide new insights into mechanisms underlying APA.

Unraveling the role of MiR-181 in skin fibrosis pathogenesis by targeting NUDT21.

Systemic sclerosis (SSc) is a life-threatening autoimmune disease characterized by widespread fibrosis in the skin and several internal organs. Nudix Hydrolase 21 (NUDT2 or CFIm25) downregulation in fibroblasts is known to play detrimental roles in both skin and lung fibrosis. This study aims to investigate the upstream mechanisms that lead to NUDT21 repression in skin fibrosis. We identified transforming growth factor ß (TGFß1) as the primary cytokine that downregulated NUDT21 in normal skin fibroblasts. In the bleomycin-induced dermal fibrosis model, consistent with the peak activation of TGFß1 at the late fibrotic stage, NUDT21 was downregulated at this stage, and delayed NUDT21 knockdown during this fibrotic phase led to enhanced fibrotic response to bleomycin. Further investigation suggested TGFß downregulated NUDT21 through microRNA (miRNA) 181a and 181b induction. Both miR-181a and miR-181b were elevated in bleomycin-induced skin fibrosis in mice and primary fibroblasts isolated from SSc patients, and they directly targeted NUDT21 and led to its downregulation in skin fibroblasts. Functional studies demonstrated that miR-181a and miR-181b inhibitors attenuated bleomycin-induced skin fibrosis in mice in association with decreased NUDT21 expression, while miR-181a and miR-181b mimics promoted bleomycin-induced fibrosis. Overall, these findings suggest a novel role for miR-181a/b in SSc pathogenesis by repressing NUDT21 expression.

Crosstalk among Alternative Polyadenylation, Genetic Variants and Ubiquitin Modification Contribute to Lung Adenocarcinoma Risk.

Ubiquitin modification and alternative polyadenylation play crucial roles in the onset and progression of cancer. Hence, this study aims to comprehensively and deeply understand gene regulation and associated biological processes in lung adenocarcinoma (LUAD) by integrating both mechanisms. Alternative polyadenylation (APA)-related E3 ubiquitin ligases in LUAD were identified through multiple databases, and the association between selected genetic loci influencing gene expression (apaQTL-SNPs) and LUAD risk were evaluated through the GWAS database of the Female Lung Cancer Consortium in Asia (FLCCA). Subsequently, the interaction between RNF213 and ZBTB20, as well as their functional mechanisms in LUAD, were investigated using bioinformatics analysis, Western blot, co-immunoprecipitation, and colony formation experiments. A total of five apaQTL-SNPs (rs41301932, rs4494603, rs9890400, rs56066320, and rs41301932), located on RNF213, were significantly associated with LUAD risk (p < 0.05), and they inhibit tumor growth through ubiquitin-mediated degradation of ZBTB20.

Genetics of cell-type-specific post-transcriptional gene regulation during human neurogenesis.

The function of some genetic variants associated with brain-relevant traits has been explained through colocalization with expression quantitative trait loci (eQTL) conducted in bulk postmortem adult brain tissue. However, many brain-trait associated loci have unknown cellular or molecular function. These genetic variants may exert context-specific function on different molecular phenotypes including post-transcriptional changes. Here, we identified genetic regulation of RNA editing and alternative polyadenylation (APA) within a cell-type-specific population of human neural progenitors and neurons. More RNA editing and isoforms utilizing longer polyadenylation sequences were observed in neurons, likely due to higher expression of genes encoding the proteins mediating these post-transcriptional events. We also detected hundreds of cell-type-specific editing quantitative trait loci (edQTLs) and alternative polyadenylation QTLs (apaQTLs). We found colocalizations of a neuron edQTL in CCDC88A with educational attainment and a progenitor apaQTL in EP300 with schizophrenia, suggesting that genetically mediated post-transcriptional regulation during brain development leads to differences in brain function.

Involvement of oncomiRs miR-23, miR-24, and miR-27 in the regulation of alternative polyadenylation in glioblastoma via CFIm25 cleavage factor.

Glioblastoma multiforme (GBM) is a highly aggressive brain tumor with a poor prognosis. The cleavage factor Im 25 (CFIm25), a crucial component of the CFIm complex, plays a key role in regulating the length of the mRNA 3'-UTR and has been implicated in various cancers, including GBM. This study sought to investigate the regulatory influence of specific microRNAs (miRNAs) on CFIm25 expression in GBM, a highly aggressive brain tumor. Bioinformatics analysis identified miRNA candidates targeting CFIm25 mRNA, and gene expression profiles from the NCBI database (GSE90603) were used for further analysis. Expression levels of CFIm25 and selected miRNAs were assessed using qRT-PCR in GBM clinical samples (n = 20) and non-malignant brain tissues (n = 5). Additionally, the MTT assay was performed to examine the effect of miRNA overexpression on U251 cell viability. Lentivectors expressing the identified miRNAs were employed to experimentally validate their regulatory role on CFIm25 in U251 cell lines, and Western blot analysis was conducted to determine CFIm25 protein levels. We observed significantly increased levels of miR-23, miR-24, and miR-27 expression, associated with a marked reduction in CFIm25 expression in GBM samples compared to non-malignant brain tissues. In particular, overexpression of miR-23, miR-24, and miR-27 in U251 cells resulted in CFIm25 downregulation at both the mRNA and protein levels, while their inhibition increased CFIm25 and reduced cell proliferation. These observations strongly implicate miR-23, miR-24, and miR-27 in regulating CFIm25 expression in GBM, emphasizing their potential as promising therapeutic targets for enhancing treatment responses in glioblastoma.

Leveraging multi-omics data to infer regulators of mRNA 3' end processing in glioblastoma.

Alterations in mRNA 3' end processing and polyadenylation are widely implicated in the biology of many cancer types, including glioblastoma (GBM), one the most aggressive tumor types. Although several RNA-binding proteins (RBPs) responsible for alternative polyadenylation (APA) were identified from functional studies in cell lines, their contribution to the APA landscape in tumors in vivo was not thoroughly addressed. In this study we analyzed a large RNA-seq data set of glioblastoma (GBM) samples from The Cancer Genome Atlas (TCGA) to identify APA patterns differentiating the main molecular subtypes of GBM. We superimposed these to RBP footprinting data and to APA events occurring upon depletion of individual RBPs from a large panel tested by the ENCODE Consortium. Our analysis revealed 22 highly concordant and statistically significant RBP-APA associations, whereby changes in RBP expression were accompanied by APA in both TCGA and ENCODE datasets. Among these, we found a previously unknown PTBP1-regulated APA event in the PRRC2B gene and an HNRNPU-regulated event in the SC5D gene. Both of these were further supported by RNA-sequencing data of paired tumor center-periphery GBM samples obtained at the University Hospital of Basel. In addition, we validated the regulation of APA in PRRC2B by PTBP1 in siRNA-knockdown and overexpression experiments followed by RNA-sequencing in two glioblastoma cell lines. The transcriptome analysis workflow that we present here enables the identification of concordant RBP-APA associations in cancers.

Light regulates widespread plant alternative polyadenylation through the chloroplast.

Transcription of eukaryotic protein-coding genes generates immature mRNAs that are subjected to a series of processing events, including capping, splicing, cleavage, and polyadenylation (CPA), and chemical modifications of bases. Alternative polyadenylation (APA) greatly contributes to mRNA diversity in the cell. By determining the length of the 3' untranslated region, APA generates transcripts with different regulatory elements, such as miRNA and RBP binding sites, which can influence mRNA stability, turnover, and translation. In the model plant Arabidopsis thaliana, APA is involved in the control of seed dormancy and flowering. In view of the physiological importance of APA in plants, we decided to investigate the effects of light/dark conditions and compare the underlying mechanisms to those elucidated for alternative splicing (AS). We found that light controls APA in approximately 30% of Arabidopsis genes. Similar to AS, the effect of light on APA requires functional chloroplasts, is not affected in mutants of the phytochrome and cryptochrome photoreceptor pathways, and is observed in roots only when the communication with the photosynthetic tissues is not interrupted. Furthermore, mitochondrial and TOR kinase activities are necessary for the effect of light. However, unlike AS, coupling with transcriptional elongation does not seem to be involved since light-dependent APA regulation is neither abolished in mutants of the TFIIS transcript elongation factor nor universally affected by chromatin relaxation caused by histone deacetylase inhibition. Instead, regulation seems to correlate with changes in the abundance of constitutive CPA factors, also mediated by the chloroplast.

Neuronal RNA processing: cross-talk between transcriptional regulation and RNA-binding proteins.

In the nervous system, alternative RNA processing is particularly prevalent, which results in the expression of thousands of transcript variants found in no other tissue. Neuron-specific RNA-binding proteins co-transcriptionally regulate alternative splicing, alternative polyadenylation, and RNA editing, thereby shaping the RNA identity of nervous system cells. Recent evidence suggests that interactions between RNA-binding proteins and cis-regulatory elements such as promoters and enhancers play a role in the determination of neuron-specific expression profiles. Here, we discuss possible mechanisms through which transcription and RNA processing cross-talk to generate the uniquely complex neuronal transcriptome, with a focus on alternative 3'-end formation.

A developmental mechanism to regulate alternative polyadenylation in an adult stem cell lineage.

Alternative cleavage and polyadenylation (APA) often results in production of mRNA isoforms with either longer or shorter 3' UTRs from the same genetic locus, potentially impacting mRNA translation, localization, and stability. Developmentally regulated APA can thus make major contributions to cell type-specific gene expression programs as cells differentiate. During Drosophila spermatogenesis, ∼500 genes undergo APA when proliferating spermatogonia differentiate into spermatocytes, producing transcripts with shortened 3' UTRs, leading to profound stage-specific changes in the proteins expressed. The molecular mechanisms that specify usage of upstream polyadenylation sites in spermatocytes are thus key to understanding the changes in cell state. Here, we show that upregulation of PCF11 and Cbc, the two components of cleavage factor II (CFII), orchestrates APA during Drosophila spermatogenesis. Knockdown of PCF11 or cbc in spermatocytes caused dysregulation of APA, with many transcripts normally cleaved at a proximal site in spermatocytes now cleaved at their distal site, as in spermatogonia. Forced overexpression of CFII components in spermatogonia switched cleavage of some transcripts to the proximal site normally used in spermatocytes. Our findings reveal a developmental mechanism where changes in expression of specific cleavage factors can direct cell type-specific APA at selected genes.

Source leaves are regulated by sink strengths through non-coding RNAs and alternative polyadenylation in cucumber (Cucumis sativus L.).

BACKGROUND: The yield of major crops is generally limited by sink capacity and source strength. Cucumber is a typical raffinose family oligosaccharides (RFOs)-transporting crop. Non-coding RNAs and alternative polyadenylation (APA) play important roles in the regulation of growth process in plants. However, their roles on the sink‒source regulation have not been demonstrated in RFOs-translocating species. RESULTS: Here, whole-transcriptome sequencing was applied to compare the leaves of cucumber under different sink strength, that is, no fruit-carrying leaves (NFNLs) and fruit-carrying leaves (FNLs) at 12th node from the bottom. The results show that 1101 differentially expressed (DE) mRNAs, 79 DE long non-coding RNAs (lncRNAs) and 23 DE miRNAs were identified, which were enriched in photosynthesis, energy production and conversion, plant hormone signal transduction, starch and carbohydrate metabolism and protein synthesis pathways. Potential co-expression networks like, DE lncRNAs-DE mRNAs/ DE miRNAs-DE mRNAs, and competing endogenous RNA (ceRNA) regulation models (DE lncRNAs-DE miRNAs-DE mRNAs) associated with sink‒source allocation, were constructed. Furthermore, 37 and 48 DE genes, which enriched in MAPK signaling and plant hormone signal transduction pathway, exist differentially APA, and SPS (CsaV3_2G033300), GBSS1 (CsaV3_5G001560), ERS1 (CsaV3_7G029600), PNO1 (CsaV3_3G003950) and Myb (CsaV3_3G022290) may be regulated by both ncRNAs and APA between FNLs and NFNLs, speculating that ncRNAs and APA are involved in the regulation of gene expression of cucumber sink‒source carbon partitioning. CONCLUSIONS: These results reveal a comprehensive network among mRNAs, ncRNAs, and APA in cucumber sink-source relationships. Our findings also provide valuable information for further research on the molecular mechanism of ncRNA and APA to enhance cucumber yield.

Dynamic regulation of alternative polyadenylation by PQBP1 during neurogenesis.

Alternative polyadenylation (APA) is a critical post-transcriptional process that generates mRNA isoforms with distinct 3' untranslated regions (3' UTRs), thereby regulating mRNA localization, stability, and translational efficiency. Cell-type-specific APA extensively shapes the diversity of the cellular transcriptome, particularly during cell fate transition. Despite its recognized significance, the precise regulatory mechanisms governing cell-type-specific APA remain unclear. In this study, we uncover PQBP1 as an emerging APA regulator that actively maintains cell-specific APA profiles in neural progenitor cells (NPCs) and delicately manages the equilibrium between NPC proliferation and differentiation. Multi-omics analysis shows that PQBP1 directly interacts with the upstream UGUA elements, impeding the recruitment of the CFIm complex and influencing polyadenylation site selection within genes associated with the cell cycle. Our findings elucidate the molecular mechanism by which PQBP1 orchestrates dynamic APA changes during neurogenesis, providing valuable insights into the precise regulation of cell-type-specific APA and the underlying pathogenic mechanisms in neurodevelopmental disorders.

Single-cell landscape of alternative polyadenylation in human lymphoid hematopoiesis.

Alternative polyadenylation (APA) is an essential post-transcriptional process that produces mature mRNA isoforms by regulating the usage of polyadenylation sites (PASs). APA is involved in lymphocyte activation; however, its role throughout the entire differentiation trajectory remains elusive. Here, we analyzed single-cell 3'-end transcriptome data from healthy subjects to construct a dynamic-APA landscape from hematopoietic stem and progenitor cells (HSPCs) to terminally differentiated lymphocytes. This analysis covered 19973 cells of 12 clusters from five lineages (B cells, CD4+ T cells, CD8+ T cells, natural killer cells, and plasmacytoid dendritic cells). A total of 2364 genes exhibited differential 3'UTR PAS usage, and 3021 genes displayed differential intronic cleavage during lymphoid differentiation. We observed a global trend of 3'UTR shortening during lymphoid differentiation. Nevertheless, specific events of both 3'UTR shortening and lengthening were also identified within each cluster. The APA patterns delineated three differentiation stages: HSPCs, precursor cells, and mature cells. Moreover, we demonstrated that the conversion of naïve T cells to memory T cells was accompanied by dynamic APA in transcription factor-encoding genes (TCF7 and NFATC2IP), immune function-related genes (BCL2, CD5, CD28, GOLT1B, and TMEM59), and protein ubiquitination-related genes (UBE2G1, YPEL5, and SUMO3). These findings expand our understanding of the underlying molecular mechanisms of APA and facilitate studies on the regulatory role of APA in lymphoid hematopoiesis.

Dynamic Alternative Polyadenylation during Litopenaeus Vannamei Metamorphosis Development.

As an important mechanism in the post-transcriptional regulation of eukaryotic gene expression, alternative polyadenylation (APA) plays a key role in biological processes such as cell proliferation and differentiation. However, the role and dynamic pattern of APA during Litopenaeus vannamei metamorphosis are poorly understood. Here, RNA-seq data covering from the embryo to the maturation (16 time points) of L. vannamei were utilized. We identified 247 differentially expressed APA events between early and adult stages, and through fuzzy mean clustering analysis, we discovered five dynamic APA patterns. Among them, the gradual elongation of the 3'UTR is the major APA pattern that changes over time, and its genes are enriched in the pathways of protein and energy metabolism. Finally, we constructed mRNA-miRNA and PPI networks and detected several central miRNAs that may regulate L. vannamei development. Our results revealed the complex APA mechanisms in L. vannamei metamorphosis, shedding new light on post-transcriptional regulation of crustacean metamorphosis.

Life stage-specific poly(A) site selection regulated by Trypanosoma brucei DRBD18.

The kinetoplastid parasite, Trypanosoma brucei, undergoes a complex life cycle entailing slender and stumpy bloodstream forms in mammals and procyclic and metacyclic forms (MFs) in tsetse fly hosts. The numerous gene regulatory events that underlie T. brucei differentiation between hosts, as well as between active and quiescent stages within each host, take place in the near absence of transcriptional control. Rather, differentiation is controlled by RNA-binding proteins (RBPs) that associate with mRNA 3' untranslated regions (3'UTRs) to impact RNA stability and translational efficiency. DRBD18 is a multifunctional T. brucei RBP, shown to impact mRNA stability, translation, export, and processing. Here, we use single-cell RNAseq to characterize transcriptomic changes in cell populations that arise upon DRBD18 depletion, as well as to visualize transcriptome-wide alterations to 3'UTR length. We show that in procyclic insect stages, DRBD18 represses expression of stumpy bloodstream form and MF transcripts. Additionally, DRBD18 regulates the 3'UTR lengths of over 1,500 transcripts, typically promoting the use of distal polyadenylation sites, and thus the inclusion of 3'UTR regulatory elements. Remarkably, comparison of polyadenylation patterns in DRBD18 knockdowns with polyadenylation patterns in stumpy bloodstream forms shows numerous similarities, revealing a role for poly(A) site selection in developmental gene regulation, and indicating that DRBD18 controls this process for a set of transcripts. RNA immunoprecipitation supports a direct role for DRBD18 in poly(A) site selection. This report highlights the importance of alternative polyadenylation in T. brucei developmental control and identifies a critical RBP in this process.

SIZ1-mediated SUMOylation of CPSF100 promotes plant thermomorphogenesis by controlling alternative polyadenylation.

Under warm temperatures, plants adjust their morphologies for environmental adaption via precise gene expression regulation. However, the function and regulation of alternative polyadenylation (APA), an important fine-tuning of gene expression, remains unknown in plant thermomorphogenesis. In this study, we found that SUMOylation, a critical post-translational modification, is induced by a long-term treatment at warm temperatures via a SUMO ligase SIZ1 in Arabidopsis. Disruption of SIZ1 altered the global usage of polyadenylation signals and affected the APA dynamic of thermomorphogenesis-related genes. CPSF100, a key subunit of the CPSF complex for polyadenylation regulation, is SUMOylated by SIZ1. Importantly, we demonstrated that SUMOylation is essential for the function of CPSF100 in genome-wide polyadenylation site choice during thermomorphogenesis. Further analyses revealed that the SUMO conjugation on CPSF100 attenuates its interaction with two isoforms of its partner CPSF30, increasing the nuclear accumulation of CPSF100 for polyadenylation regulation. In summary, our study uncovers a regulatory mechanism of APA via SIZ1-mediated SUMOylation in plant thermomorphogenesis.

2'-O-methylation at internal sites on mRNA promotes mRNA stability.

2'-O-methylation (Nm) is a prominent RNA modification well known in noncoding RNAs and more recently also found at many mRNA internal sites. However, their function and base-resolution stoichiometry remain underexplored. Here, we investigate the transcriptome-wide effect of internal site Nm on mRNA stability. Combining nanopore sequencing with our developed machine learning method, NanoNm, we identify thousands of Nm sites on mRNAs with a single-base resolution. We observe a positive effect of FBL-mediated Nm modification on mRNA stability and expression level. Elevated FBL expression in cancer cells is associated with increased expression levels for 2'-O-methylated mRNAs of cancer pathways, implying the role of FBL in post-transcriptional regulation. Lastly, we find that FBL-mediated 2'-O-methylation connects to widespread 3' UTR shortening, a mechanism that globally increases RNA stability. Collectively, we demonstrate that FBL-mediated Nm modifications at mRNA internal sites regulate gene expression by enhancing mRNA stability.

Proximal termination generates a transcriptional state that determines the rate of establishment of Polycomb silencing.

The mechanisms and timescales controlling de novo establishment of chromatin-mediated transcriptional silencing by Polycomb repressive complex 2 (PRC2) are unclear. Here, we investigate PRC2 silencing at Arabidopsis FLOWERING LOCUS C (FLC), known to involve co-transcriptional RNA processing, histone demethylation activity, and PRC2 function, but so far not mechanistically connected. We develop and test a computational model describing proximal polyadenylation/termination mediated by the RNA-binding protein FCA that induces H3K4me1 removal by the histone demethylase FLD. H3K4me1 removal feeds back to reduce RNA polymerase II (RNA Pol II) processivity and thus enhance early termination, thereby repressing productive transcription. The model predicts that this transcription-coupled repression controls the level of transcriptional antagonism to PRC2 action. Thus, the effectiveness of this repression dictates the timescale for establishment of PRC2/H3K27me3 silencing. We experimentally validate these mechanistic model predictions, revealing that co-transcriptional processing sets the level of productive transcription at the locus, which then determines the rate of the ON-to-OFF switch to PRC2 silencing.

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.