HNRNPA1-mediated 3' UTR length changes of HN1 contributes to cancer- and senescence-associated phenotypes

Transcriptomic analysis reveals regulation of adipogenesis via long non-coding RNA, alternative splicing, and alternative polyadenylation

Obesity is characterized by dysregulated adipogenesis that leads to increased number and/or size of adipocytes. Understanding the molecular mechanisms governing adipogenesis is therefore key to designing therapeutic interventions against obesity. In our study, we analyzed 3'-end sequencing data that we generated from human preadipocytes and adipocytes, as well as previously published RNA-seq datasets, to elucidate mechanisms of regulation via long non-coding RNA (lncRNA), alternative splicing (AS) and alternative polyadenylation (APA). We discovered lncRNAs that have not been previously characterized but may be key regulators of white adipogenesis. We also detected 100 AS events and, using motif enrichment analysis, identified RNA binding proteins (RBPs) that could mediate exon skipping-the most prevalent AS event. In addition, we show that usage of alternative poly(A) sites in introns or 3'-UTRs of key adipogenesis genes leads to isoform diversity, which can have significant biological consequences on differentiation efficiency. We also identified RBPs that may modulate APA and defined how 3'-UTR APA can regulate gene expression through gain or loss of specific microRNA binding sites. Taken together, our bioinformatics-based analysis reveals potential therapeutic avenues for obesity through manipulation of lncRNA levels and the profile of mRNA isoforms via alternative splicing and polyadenylation.

Single-cell RNA sequencing reveals blastomere heterogeneity of 2-cell embryos in pigs

In mammals, single blastomeres from as early as 2-cell embryos demonstrate heterogeneous developmental capacity and fate decision into different cell lineages. However, mechanisms underlying blastomere heterogeneity of 2-cell embryos remain largely unresolved. Here, we analysed the molecular heterogeneity of full-length mRNAs and their 3'UTR regions, based on the single-cell RNA-seq data of pig 2-cell embryos generated from in vivo fertilization (in vivo), in vitro fertilization (in vitro) and parthenogenetic activation (PA), respectively. First, unsupervised clustering helped discover two different groups of blastomeres for 2-cell pig embryos. Between these two groups of blastomeres in pig 2-cell embryos, 35, 301 and 428 full-length mRNAs respectively in in vivo, in vitro and PA embryo types were identified to be differentially expressed (padj  ≤ .05 and |log2 [fold change]| ≥1) (DE mRNAs), while 92, 89 and 42 mRNAs were shown to be with significantly different 3'UTR lengths (3'UTR DE) (padj  ≤ .05). Gene enrichment for both DE mRNAs and 3'UTR DE mRNAs found multiple signalling pathways, including cell cycle, RNA processing. Few numbers of common DE mRNAs and 3'UTR DE mRNAs existed between in vitro and in vivo blastomeres derived from 2-cell embryos, indicating the larger differences between in vitro and in vivo fertilized embryos. Integrative genomics viewer analysis further identified that 3'UTRs of HSDL2 and SGTA (in vivo), FAM204A and phosphoserine phosphatase (in vitro), PRPF40A and RPIA (PA) had >100 nt average length changes. Moreover, numbers and locations of regulatory elements (polyadenylation site, cytoplasmic polyadenylation element and microRNA binding sites) within 3'UTRs of these DE mRNAs were predicted. These results indicate that molecular heterogeneity existed among blastomeres from different types of pig 2-cell embryos, providing useful information and novel insights into future functional investigation on its relationship with the subsequent embryo development and differentiation.

How Driver Oncogenes Shape and Are Shaped by Alternative Splicing Mechanisms in Tumors

The development of RNA sequencing methods has allowed us to study and better understand the landscape of aberrant pre-mRNA splicing in tumors. Altered splicing patterns are observed in many different tumors and affect all hallmarks of cancer: growth signal independence, avoidance of apoptosis, unlimited proliferation, invasiveness, angiogenesis, and metabolism. In this review, we focus on the interplay between driver oncogenes and alternative splicing in cancer. On one hand, oncogenic proteins-mutant p53, CMYC, KRAS, or PI3K-modify the alternative splicing landscape by regulating expression, phosphorylation, and interaction of splicing factors with spliceosome components. Some splicing factors-SRSF1 and hnRNPA1-are also driver oncogenes. At the same time, aberrant splicing activates key oncogenes and oncogenic pathways: p53 oncogenic isoforms, the RAS-RAF-MAPK pathway, the PI3K-mTOR pathway, the EGF and FGF receptor families, and SRSF1 splicing factor. The ultimate goal of cancer research is a better diagnosis and treatment of cancer patients. In the final part of this review, we discuss present therapeutic opportunities and possible directions of further studies aiming to design therapies targeting alternative splicing mechanisms in the context of driver oncogenes.

High Expression of Heterogeneous Nuclear Ribonucleoprotein A1 Facilitates Hepatocellular Carcinoma Growth

Purpose

Hepatocellular carcinoma (HCC) represents one of the most common tumors in the world. Our study aims to explore new markers and therapeutic targets for HCC. Heterogeneous Nuclear ribonucleoprotein A1 (hnRNPA1) has recently been found to be involved in the progression of several types of cancer, but its role in HCC remains uncovered.

Methods

We performed bioinformatic analysis to preliminarily show the relationship between hnRNPA1 and liver cancer. Then the correlation of the hnRNPA1 gene expression with clinicopathological characteristics of HCC patients was verified by human liver cancer tissue microarrays. The functional role of this gene was evaluated by in vivo and vitro experiments.

Results

Results showed that the expression of hnRNPA1 was upregulated in HCC tissues and was associated with pathological stage of HCC patients. Knockdown of hnRNPA1 gene markedly inhibited tumor growth in vivo, and reversed the effects on proliferation, migration and invasion and promoted apoptosis in vitro. Furthermore, down-regulation of hnRNPA1 gene expression can inhibit the activity of the MEK/ERK pathway.

Conclusion

In our work, we combined bioinformatic analysis with in vivo and in vitro experiments to initially elucidate the function of hnRNPA1 in liver cancer, which may help to explore biomarkers and therapeutic targets for HCC patients.

Global 3'-UTRome of porcine immature Sertoli cells altered by acute heat stress

Alternative polyadenylation (APA) affects the composition of cis-elements in 3'-untranslated region (3'-UTR), to regulate gene expression and localization, and subsequently the downstream biological processes. Acute heat stress could change rapidly the cellular transcriptome, however the underlying molecular changes are less explored. Here, we systematically catalogued the global 3'-UTRome dynamics, by analyzing our previously reported transcriptome sequencing data of porcine immature Sertoli (iST) cells before (Control group), acute heat stress treatment at 43 °C for 0.5h (HS0.5 group), and 36h recovery culture (HS0.5-R36h group) after acute heat stress treatment. After three group comparisons (HS0.5 vs. Control, HS0.5-R36 vs. HS0.5, and HS0.5-R36 vs. Control), DaPars (dynamic analysis of alternative polyadenylation) identified 639, 464 and 290 mRNAs, and APAtrap (a tool to identify APA sites and detect changes of APA site usage) identified 713, 518 and 321 mRNAs, with significantly different 3'-UTRs (Padj.≤0.05), respectively. These genes with different 3'-UTR patterns were mainly enriched in P53, glycolysis/gluconeogenesis, HIF-1, apoptosis, PI3K-Akt and AMPK signaling pathways. Further analysis identified that average 3'-UTR lengths of Acss2, Inpp1 and Nr1h4 were more than 140 nt longer (HS0.5-R36 vs. HS0.5), and contained different cis-elements (PAS, CPE and microRNA binding sites). Moreover, Hsp70.2, Inhbb and Dhrs were identified to have extremely different 3'-UTR abundances. Further 3'RACE assays validated several 3'-UTRs of Nr1h4, and RT-qPCR confirmed the abundance changes of different 3'-UTR isoforms for Nr1h4 and Hsp70.2. Our findings provide useful information and resources to further uncover the molecular role of 3'-UTR, in regulating the response of porcine iST cells to acute heat stress.

A proteomic profile of the healthy human placenta

BACKGROUND

The placenta remains one of the least studied organs within the human body. Yet, placental dysfunction has been associated with various pregnancy complications leading to both maternal and fetal death and long-term health consequences. The aim of this study was to characterise the protein networks of healthy term placental sub-anatomical regions using label free quantification mass spectrometry.

METHODS

Three healthy placentae were sampled at five sample sites and each biopsy was dissected into maternal-, middle-, and fetal- sub-anatomical regions. Quadrupole-orbitrap mass spectrometer was used in data dependant analysis mode to identify 1859 unique proteins before detailed differential expression between regions.

RESULTS

Protein profiling identified 1081, 1086, and 1101 proteins in maternal, middle, and fetal sub-anatomical regions respectively. Differentially expressed proteins were identified considering the effect between sample site location and sub-anatomical region on protein expression. Of these, 374 differentially expressed proteins (Two-way ANOVA adjusted p-value < 0.05, HSD Tukey adjusted p-value 0.05) were identified between sample site locations and sub-anatomical regions. The placenta specific disease map NaviCenta ( https://www.sbi.uni-rostock.de/minerva/index.xhtml?id=NaviCenta ) was used to focus functional analysis results to the placenta specific context. Subsequently, functional analysis with a focus on senescence, and mitochondrial function were performed. Significant differences were observed between sub-anatomical regions in protein intensity and composition. A decrease in anti-senescent proteins within the maternal sub-anatomical region, and an increase in proteins associated with a switch from ATP to fatty acid consumption as a source of energy between middle and fetal sub-anatomical regions were observed.

CONCLUSION

These results suggest that normal proteomic variations exist within the anatomical structure of the placenta, thus recommending serial sectioning methodology for consistent placental research.

hnRNP A1 in RNA metabolism regulation and as a potential therapeutic target

Abnormal RNA metabolism, regulated by various RNA binding proteins, can have functional consequences for multiple diseases. Heterogeneous nuclear ribonucleoprotein A1 (hnRNP A1) is an important RNA binding protein, that regulates various RNA metabolic processes, including transcription, alternative splicing of pre-mRNA, translation, miRNA processing and mRNA stability. As a potent splicing factor, hnRNP A1 can regulate multiple splicing events, including itself, collaborating with other cooperative or antagonistical splicing factors by binding to splicing sites and regulatory elements in exons or introns. hnRNP A1 can modulate gene transcription by directly interacting with promoters or indirectly impacting Pol II activities. Moreover, by interacting with the internal ribosome entry site (IRES) or 3′-UTR of mRNAs, hnRNP A1 can affect mRNA translation. hnRNP A1 can alter the stability of mRNAs by binding to specific locations of 3′-UTR, miRNAs biogenesis and Nonsense-mediated mRNA decay (NMD) pathway. In this review, we conclude the selective sites where hnRNP A1 binds to RNA and DNA, and the co-regulatory factors that interact with hnRNP A1. Given the dysregulation of hnRNP A1 in diverse diseases, especially in cancers and neurodegeneration diseases, targeting hnRNP A1 for therapeutic treatment is extremely promising. Therefore, this review also provides the small-molecule drugs, biomedicines and novel strategies targeting hnRNP A1 for therapeutic purposes.

Tumor-Promoting Actions of HNRNP A1 in HCC Are Associated with Cell Cycle, Mitochondrial Dynamics, and Necroptosis

Hepatocellular carcinoma (HCC) is one of the most frequent malignancies in the world. Although increasing evidence supports the role of heterogeneous ribonucleoprotein particle A1 (HNRNP A1) in tumor progression, the function of HNRNP A1 in HCC remains unclear. Here, we focused on the role of HNRNP A1 in the development of HCC. In this study, we found HNRNP A1 participates in many aspects of HCC, such as progression and prognosis. Our results showed that HNRNP A1 is upregulated in human HCC tissues and cell lines. High expression of HNRNP A1 can promote the proliferation, migration, and invasion in HCC cells and accelerate tumor progression in mice. Moreover, we found that HNRNP A1 prevents the senescence process of HCC cells. Knocking down of HNRNP A1 promotes the expression of P16^INK4, which arrests the cell cycle and then induces the senescence phenotype in HCC cells. Furthermore, we found that HNRNP A1 regulated necroptosis and mitochondrial dynamics. In summary, our study indicates that HNRNP A1 promotes the development of HCC, which suggests a potential therapeutic target for HCC.

Distinct, opposing functions for CFIm59 and CFIm68 in mRNA alternative polyadenylation of Pten and in the PI3K/Akt signalling cascade

Precise maintenance of PTEN dosage is crucial for tumor suppression across a wide variety of cancers. Post-transcriptional regulation of Pten heavily relies on regulatory elements encoded by its 3'UTR. We previously reported the important diversity of 3'UTR isoforms of Pten mRNAs produced through alternative polyadenylation (APA). Here, we reveal the direct regulation of Pten APA by the mammalian cleavage factor I (CFIm) complex, which in turn contributes to PTEN protein dosage. CFIm consists of the UGUA-binding CFIm25 and APA regulatory subunits CFIm59 or CFIm68. Deep sequencing analyses of perturbed (KO and KD) cell lines uncovered the differential regulation of Pten APA by CFIm59 and CFIm68 and further revealed that their divergent functions have widespread impact for APA in transcriptomes. Differentially regulated genes include numerous factors within the phosphoinositide 3-kinase (PI3K)/protein kinase B (Akt) signalling pathway that PTEN counter-regulates. We further reveal a stratification of APA dysregulation among a subset of PTEN-driven cancers, with recurrent alterations among PI3K/Akt pathway genes regulated by CFIm. Our results refine the transcriptome selectivity of the CFIm complex in APA regulation, and the breadth of its impact in PTEN-driven cancers.

Role of noncoding RNAs and untranslated regions in cancer: A review

Cancer is one of the most prevalent diseases worldwide, and poses a threat to human health. Noncoding RNAs (ncRNAs) constitute most transcripts, but they cannot be translated into proteins. Studies have shown that ncRNAs can act as tumor suppressors or oncogenes. This review describes the role of several ncRNAs in various cancers, including microRNAs (miRNAs) such as the miR-34 family, let-7, miR-17-92 cluster, miR-210, and long noncoding RNAs (lncRNAs) such as HOX transcript antisense intergenic RNA (HOTAIR), Metastasis associated lung adenocarcinoma transcript 1 (MALAT1), H19, NF-κB-interacting lncRNA (NKILA), as well as circular RNAs (circRNAs) and untranslated regions (UTRs), highlighting their effects on cancer growth, invasion, metastasis, angiogenesis, and apoptosis. They function as tumor suppressors or oncogenes that interfere with different axes and pathways, including p53 and IL-6, which are involved in the progression of cancer. The characteristic expression of some ncRNAs in cancer also allows them to be used as biomarkers for early diagnosis and therapeutic candidates. There is a complex network of interactions between ncRNAs, with some lncRNAs and circRNAs acting as competitive endogenous RNAs (ceRNAs) to decoy miRNAs and repress their expression. The ceRNA network is a part of the ncRNA network and numerous ncRNAs work as nodes or hubs in the network, and disruption of their interactions can cause cancer development. Therefore, the balance and stabilization of this network are important for cancer diagnosis and treatment.

Global 3'-untranslated region landscape mediated by alternative polyadenylation during meiotic maturation of pig oocytes

Alternative polyadenylation affects the length and composition of 3'-untranslated region (3'-UTR) and regulates mRNA stability or translational activity to affect important biological processes. However, global 3'-UTR landscape and its relationship with gamete maturation remain less studied. Here, we analysed our previously reported single-cell RNA-seq data of germinal vesicle and metaphase II stage oocytes in pigs to systematically catalogue the 3'-UTR dynamics during oocyte maturation. Two softwares (DaPars and APAtrap) were employed and identified 110 and 228 mRNAs with significantly different 3'-UTRs (adjusted p ≤ .05), respectively. Gene enrichment analyses found signalling pathways related with biological processes of female gametophyte production, methyltransferase activity and mRNA surveillance (DaPars) and cell cycle process, regulation of ERK1 and ERK2 cascade, regulation of translation, spindle organization, kinetochore, condensed chromosome and progesterone-mediated oocyte maturation (APAtrap), respectively. Moreover, 18 of 110 mRNAs (|△PDUI| ≥ 0.25 and |log₂ PDUI ratio| ≥ 0.59) and 15 of 228 mRNAs (Perc. diff. ≥ 0.5) were with greater difference of 3'-UTR length or abundance, and integrative genomics viewer analysis further identified 4 (Alg10, Hadhb, Hsd17b4 and Sbds) of 18 mRNAs to be with 3'-UTR length differed ≥150 bp and 6 (Gcc1, Hnrnpa2b1, Lsm6, Prpf18, Sfr1 and Ust) of 15 mRNAs to be with 3'-UTR abundance extremely differed. Furthermore, the location, sequences and number of cis-elements were predicted, which were shown to derange cytoplasmic polyadenylation element, poly(A) site and microRNA binding sites within 3'-UTRs of Alg10, Hadhb, Hsd17b4 and Sbds mRNAs. Taken together, global 3'-UTR landscape changes dynamically with oocyte meiotic maturation, potentially involved in regulating oocyte meiotic process in pigs.

Alternative polyadenylation: An enigma of transcript length variation in health and disease

Alternative polyadenylation (APA) is a molecular mechanism during a pre-mRNA processing that involves usage of more than one polyadenylation site (PA-site) generating transcripts of varying length from a single gene. The location of a PA-site affects transcript length and coding potential of an mRNA contributing to both mRNA and protein diversification. This variation in the transcript length affects mRNA stability and translation, mRNA subcellular and tissue localization, and protein function. APA is now considered as an important regulatory mechanism in the pathophysiology of human diseases. An important consequence of the changes in the length of 3'-untranslated region (UTR) from disease-induced APA is altered protein expression. Yet, the relationship between 3'-UTR length and protein expression remains a paradox in a majority of diseases. Here, we review occurrence of APA, mechanism of PA-site selection, and consequences of transcript length variation in different diseases. Emerging evidence reveals coordinated involvement of core RNA processing factors including poly(A) polymerases in the PA-site selection in diseases-associated APAs. Targeting such APA regulators will be therapeutically significant in combating drug resistance in cancer and other complex diseases. This article is categorized under: RNA Processing > 3' End Processing RNA in Disease and Development > RNA in Disease Translation > Regulation.

Construction of Two Alternative Polyadenylation Signatures to Predict the Prognosis of Sarcoma Patients

Background

Increasing evidence indicates that alternative polyadenylation (APA) is associated with the prognosis of cancers.

Methods

We obtained gene expression and APA profiles of 259 sarcoma patients from the TCGA dataportal and TC3A database, respectively. The prognostic signatures, clinical nomograms, and regulatory networks were studied by integrated bioinformatics analyses. Then, the immune cell infiltration profile was obtained from the ImmuCellAI. The association between APA-based signature and immune cells was studied.

Results

A total of 61 and 38 APA events were identified as overall survival (OS)- and progress free-survival (PFS)-related biomarkers, respectively. Two signatures were generated. The area under the curves (AUC) values of OS signature were 0.900, 0.928, and 0.963 over 2-, 4-, and 6-years, respectively. And the AUC values of PFS signature at 2-, 4-, and 6-years were 0.826, 0.840, and 0.847, respectively. Overall and subgroup analyses indicated that high-risk patients had a worse prognosis than low-risk patients (all p-values < 0.05). In addition, immunomics analyses indicated that there are different patterns of immune cell infiltration between low- and high-risk patients. Furthermore, two clinical-APA nomograms were established and the C-indexes were 0.813 and 0.809 for OS nomogram and PFS nomogram, respectively. Finally, two APA regulatory networks were constructed. FIP1L1-VPS26B was identified as a key regulating relationship and validated in the pan-cancer analyses.

Conclusion

In this study, we identified prognostic predictors based on APA events with high accuracy for risk stratification in sarcoma patients and uncovered interesting regulatory networks in sarcoma that could be underlying mechanisms. This study not only provides novel potential prognostic biomarkers but promote precision medicine and provide potential novel research interests for immunotherapy.

Multiplexed In Situ Protein Profiling with High-Performance Cleavable Fluorescent Tyramide

Understanding the composition, function and regulation of complex cellular systems requires tools that quantify the expression of multiple proteins at their native cellular context. Here, we report a highly sensitive and accurate protein in situ profiling approach using off-the-shelf antibodies and cleavable fluorescent tyramide (CFT). In each cycle of this method, protein targets are stained with horseradish peroxidase (HRP) conjugated antibodies and CFT. Subsequently, the fluorophores are efficiently cleaved by mild chemical reagents, which simultaneously deactivate HRP. Through reiterative cycles of protein staining, fluorescence imaging, fluorophore cleavage, and HRP deactivation, multiplexed protein quantification in single cells in situ can be achieved. We designed and synthesized the high-performance CFT, and demonstrated that over 95% of the staining signals can be erased by mild chemical reagents while preserving the integrity of the epitopes on protein targets. Applying this method, we explored the protein expression heterogeneity and correlation in a group of genetically identical cells. With the high signal removal efficiency, this approach also enables us to accurately profile proteins in formalin-fixed paraffin-embedded (FFPE) tissues in the order of low to high and also high to low expression levels.

Alternative polyadenylation: methods, mechanism, function, and role in cancer

Occurring in over 60% of human genes, alternative polyadenylation (APA) results in numerous transcripts with differing 3'ends, thus greatly expanding the diversity of mRNAs and of proteins derived from a single gene. As a key molecular mechanism, APA is involved in various gene regulation steps including mRNA maturation, mRNA stability, cellular RNA decay, and protein diversification. APA is frequently dysregulated in cancers leading to changes in oncogenes and tumor suppressor gene expressions. Recent studies have revealed various APA regulatory mechanisms that promote the development and progression of a number of human diseases, including cancer. Here, we provide an overview of four types of APA and their impacts on gene regulation. We focus particularly on the interaction of APA with microRNAs, RNA binding proteins and other related factors, the core pre-mRNA 3'end processing complex, and 3'UTR length change. We also describe next-generation sequencing methods and computational tools for use in poly(A) signal detection and APA repositories and databases. Finally, we summarize the current understanding of APA in cancer and provide our vision for future APA related research.

RNA-biology ruling cancer progression? Focus on 3'UTRs and splicing

The protein-coding regions of mRNAs have the information to make proteins and hence have been at the center of attention for understanding altered protein functions in disease states, including cancer. Indeed, the discovery of genomic alterations and driver mutations that change protein levels and/or activity has been pivotal in our understanding of cancer biology. However, to better understand complex molecular mechanisms that are deregulated in cancers, we also need to look at non-coding parts of mRNAs, including 3'UTRs (untranslated regions), which control mRNA stability, localization, and translation efficiency. Recently, these rather overlooked regions of mRNAs are gaining attention as mounting evidence provides functional links between 3'UTRs, protein functions, and cancer-related molecular mechanisms. Here, roles of 3'UTRs in cancer biology and mechanisms that result in cancer-specific 3'-end isoform variants will be reviewed. An increased appreciation of 3'UTRs may help the discovery of new ways to explain as of yet unknown oncogene activation and tumor suppressor inactivation cases in cancers, and provide new avenues for diagnostic and therapeutic applications.

Highly Sensitive and Multiplexed Protein Imaging With Cleavable Fluorescent Tyramide Reveals Human Neuronal Heterogeneity

The ability to comprehensively profile proteins in intact tissues in situ is crucial for our understanding of health and disease. However, the existing methods suffer from low sensitivity and limited sample throughput. To address these issues, here we present a highly sensitive and multiplexed in situ protein analysis approach using cleavable fluorescent tyramide and off-the-shelf antibodies. Compared with the current methods, this approach enhances the detection sensitivity and reduces the imaging time by 1-2 orders of magnitude, and can potentially detect hundreds of proteins in intact tissues at the optical resolution. Applying this approach, we studied protein expression heterogeneity in a population of genetically identical cells, and performed protein expression correlation analysis to identify co-regulated proteins. We also profiled >6,000 neurons in a human formalin-fixed paraffin-embedded (FFPE) hippocampus tissue. By partitioning these neurons into varied cell clusters based on their multiplexed protein expression profiles, we observed different sub-regions of the hippocampus consist of neurons from distinct clusters.

Tempo-spatial alternative polyadenylation analysis reveals that 3' UTR lengthening of Mdm2 regulates p53 expression and cellular senescence in aged rat testis

Although tissue aging is accompanied with cellular senescence, it is much complicated than senescence given both types and number of cells change with age. Alternative polyadenylation (APA) had shown tissue specificity and APA-mediated 3' untranslated region (3' UTR) lengthening could regulate senescence-associated phenotypes. However, whether tissue aging shows similar trends remains unknown. Here, we performed a comprehensive analysis on RNA-seq datasets derived from multiple cells and rat tissues of young and old age. Although APA-mediated 3' UTR lengthening in various senescent cells reinforced the previous discovery, tissue aging showed much more complexity in APA. Interestingly, testis was the only tissue displaying dramatic 3' UTR lengthening and decreased expression trend of corresponding genes in aged rat. Genes with longer 3' UTR in aged testis were enriched in senescence-associated pathways, among which, Mdm2, encoding an E3 ligase of p53, favored distal poly(A) site resulting in lengthened 3' UTR and decreased expression. Longer 3' UTR of Mdm2 generated less protein, and decreased Mdm2 expression led to senescence-associated phenotypes along with increased p53 and p21 protein abundance, which could all be reversed by Mdm2 overexpression. Our work revealed complicated APA changes during tissue aging and discovered APA-mediated 3' UTR lengthening of Mdm2 is a hidden layer in regulating the well-known senescence-related p53-p21 signal axis during testis aging, and also has potential implications regarding declined male fertility along aging.