Alternative polyadenylation of mRNA and its role in cancer

RNA kinetics influence the response to transcriptional perturbation in leukaemia cell lines

Therapeutic targeting of dysregulated transcription has emerged as a promising strategy for the treatment of cancers, such as leukaemias. The therapeutic response to small molecule inhibitors of Bromodomain-Containing Proteins (BRD), such as BRD2 and BRD4, P300/cAMP-response element binding protein (CBP) and Cyclin Dependent Kinases (CDKs), is generally attributed to the selective disruption of oncogenic gene expression driven by enhancers, super-enhancers (SEs) and lineage-specific transcription factors (TFs), including the c-MYC oncogene. The selectivity of compounds targeting the transcriptional machinery may be further shaped by post-transcriptional processes. To quantitatively assess the contribution of post-transcriptional regulation in responses to transcription inhibition, we performed multi-omics analyses to accurately measure mRNA production and decay kinetics. We demonstrate that it is not only the selective disruption of mRNA production, but rather mRNA decay rates that largely influence the selectivity associated with transcriptional inhibition. Accordingly, genes down-regulated with transcriptional inhibitors are largely characterized by extremely rapid mRNA production and turnover. In line with this notion, stabilization of the c-MYC transcript through swapping of its 3' untranslated region (UTR) rendered c-MYC insensitive to transcriptional targeting. This failed to negate the impact on c-MYC downstream targets and did not abrogate therapeutic responses. Finally, we provide evidence that modulating post-transcriptional pathways, such as through ELAVL1 targeting, can sensitize long-lived mRNAs to transcriptional inhibition and be considered as a combination therapy approach in leukaemia. Taken together, these data demonstrate that mRNA kinetics influence the therapeutic response to transcriptional perturbation and can be modulated for novel therapeutic outcomes using transcriptional agents in leukaemia.

TRIM65 deficiency alleviates renal fibrosis through NUDT21-mediated alternative polyadenylation

Chronic kidney disease (CKD) is a major global health concern and the third leading cause of premature death. Renal fibrosis is the primary process driving the progression of CKD, but the mechanisms behind it are not fully understood, making treatment options limited. Here, we find that the E3 ligase TRIM65 is a positive regulator of renal fibrosis. Deletion of TRIM65 results in a reduction of pathological lesions and renal fibrosis in mouse models of kidney fibrosis induced by unilateral ureteral obstruction (UUO)- and folic acid. Through screening with a yeast-hybrid system, we identify a new interactor of TRIM65, the mammalian cleavage factor I subunit CFIm25 (NUDT21), which plays a crucial role in fibrosis through alternative polyadenylation (APA). TRIM65 interacts with NUDT21 to induce K48-linked polyubiquitination of lysine 56 and proteasomal degradation, leading to the inhibition of TGF-β1-mediated SMAD and ERK1/2 signaling pathways. The degradation of NUDT21 subsequently altered the length and sequence content of the 3'UTR (3'UTR-APA) of several pro-fibrotic genes including Col1a1, Fn-1, Tgfbr1, Wnt5a, and Fzd2. Furthermore, reducing NUDT21 expression via hydrodynamic renal pelvis injection of adeno-associated virus 9 (AAV9) exacerbated UUO-induced renal fibrosis in the normal mouse kidneys and blocked the protective effect of TRIM65 deletion. These findings suggest that TRIM65 promotes renal fibrosis by regulating NUDT21-mediated APA and highlight TRIM65 as a potential target for reducing renal fibrosis in CKD patients.

A transposable element prevents severe hemophilia B and provides insights into the evolution of new- and old world primates

Alu-elements comprise a large part of the human genome and some insertions have been shown to cause diseases. Here, we illuminate the protective role of an Alu-element in the 3'UTR of the human Factor 9 gene and its ability to ameliorate a poly(A) site mutation in a hemophilia B patient, preventing him from developing a severe disease. Using a minigene, we examined the disease-causing mutation and the modifying effect of the transposon in cellulo. Further, we simulated evolutionary scenarios regarding alternative polyadenylation before and after Alu insertion. A sequence analysis revealed that Old World monkeys displayed a highly conserved polyadenylation sites in this Alu-element, whereas New World monkeys lacked this motif, indicating a selective pressure. We conclude that this transposon has inserted shortly before the separation of Old and New World monkeys and thus also serves as a molecular landmark in primate evolution.

The roles and mechanisms of coding and noncoding RNA variations in cancer

Functional variations in coding and noncoding RNAs are crucial in tumorigenesis, with cancer-specific alterations often resulting from chemical modifications and posttranscriptional processes mediated by enzymes. These RNA variations have been linked to tumor cell proliferation, growth, metastasis, and drug resistance and are valuable for identifying diagnostic or prognostic cancer biomarkers. The diversity of posttranscriptional RNA modifications, such as splicing, polyadenylation, methylation, and editing, is particularly significant due to their prevalence and impact on cancer progression. Additionally, other modifications, including RNA acetylation, circularization, miRNA isomerization, and pseudouridination, are recognized as key contributors to cancer development. Understanding the mechanisms underlying these RNA modifications in cancer can enhance our knowledge of cancer biology and facilitate the development of innovative therapeutic strategies. Targeting these RNA modifications and their regulatory enzymes may pave the way for novel RNA-based therapies, enabling tailored interventions for specific cancer subtypes. This review provides a comprehensive overview of the roles and mechanisms of various coding and noncoding RNA modifications in cancer progression and highlights recent advancements in RNA-based therapeutic applications.

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.

A pan-cancer analysis of the core pre-mRNA 3' end processing factors, and their association with prognosis, tumor microenvironment, and potential targets

Alternative polyadenylation (APA) is a crucial mechanism for regulating gene expression during pre-mRNA 3' processing. Pre-mRNA 3' end processing factors is the main factor involved in this process. However, pre-mRNA 3' end processing factors in different cancer expression profiles and the relationship between pre-mRNA 3' end processing factors and tumor microenvironment and the prognosis of the same patient is still unclear. In this study, we conducted a comprehensive exploration of the core pre-mRNA 3' end processing factors across various cancer types by utilizing common cancer database, and revealing a robust correlation between the expression of these core factors and tumor characteristics. Leveraging advanced bioinformatics databases, we evaluated the expression levels and prognostic relevance of pre-mRNA 3' end processing factors across pan-cancer tissues. Our extensive pan-cancer analysis revealed unique expression patterns of pre-mRNA 3' end processing factors in both tumor and adjacent non-tumorous tissues. Notably, we found a significant correlation between the expression levels of pre-mRNA 3' end processing factors and patient prognosis. Furthermore, we identified strong associations between pre-mRNA 3' end processing factors expression and various factors, such as stromal, immune, RNA stemness, and DNA stemness scores across pan-cancer tissues. Our data also highlighted a link between the expression of pre-mRNA 3' end processing factors and sensitivity to specific drugs, including pyrazoloacndine, amonaflide, and chelerythrinede, among others. We found four key pre-mRNA 3' end processing factors that play a crucial role in mRNA preprocessing. Our study illuminates the potential promotion and inhibition role of pre-mRNA 3' end processing regulators in the progression of cancer, CPSF2, CPSF3, CSTF2, SYMPK offering valuable insights for future research investigations on these regulators as diagnostic markers and therapeutic targets across pan-cancer.

A pan-cancer interrogation of intronic polyadenylation and its association with cancer characteristics

3'UTR-APAs have been extensively studied, but intronic polyadenylations (IPAs) remain largely unexplored. We characterized the profiles of 22 260 IPAs in 9679 patient samples across 32 cancer types from the Cancer Genome Atlas cohort. By comparing tumor and paired normal tissues, we identified 180 ~ 4645 dysregulated IPAs in 132 ~ 2249 genes in each of 690 patient tumors from 22 cancer types that showed consistent patterns within individual cancer types. We selected 2741 genes that showed consistently patterns across cancer types, including 1834 pan-cancer tumor-enriched and 907 tumor-depleted IPA genes; the former were amply represented in the functional pathways such as deoxyribonucleic acid damage repair. Expression of IPA isoforms was associated with tumor mutation burden and patient characteristics (e.g. sex, race, cancer stages, and subtypes) in cancer-specific and feature-specific manners, and could be a more accurate prognostic marker than gene expression (summary of all isoforms). In summary, our study reveals the roles and the clinical relevance of tumor-associated IPAs.

Cleavage and polyadenylation machinery as a novel targetable vulnerability for human cancer

The role of alternative polyadenylation of mRNA in sustaining aggressive features of tumors is quite well established, as it is responsible for the 3'UTR shortening of oncogenes and subsequent relief from miRNA-mediated repression observed in cancer cells. However, the information regarding the vulnerability of cancer cells to the inhibition of cleavage and polyadenylation (CPA) machinery is very scattered. Only few recent reports show the antitumor activity of pharmacological inhibitors of CPSF3, one among CPA factors. More in general, the fact that deregulated CPA can be seen as a new hallmark of cancer and as a potential reservoir of novel therapeutic targets has never been formalized. Here, to extend our view on the potential of CPA inhibition (CPAi) approaches as anticancer therapies, we systematically tested the fitness of about one thousand cell lines of different cancer types upon depletion of all known CPA factors by interrogating genome-scale CRISPR and RNAi dependency maps of the DepMap project. Our analysis confirmed core and accessory CPA factors as novel vulnerabilities for human cancer, thus highlighting the potential of CPAi as anticancer therapy. Among all, CPSF1 appeared as a promising actionable candidate for drug development, as it showed low dependency scores pancancer and particularly in highly proliferating cells. In a personalized medicine perspective, the observed differential vulnerability of cancer cell lines to selected CPA factors may be used to build up signatures to predict response of individual human tumors to CPAi approaches.

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.

Genetic analysis of PALB2 gene WD40 domain in canine mammary tumour patients

BACKGROUND

DNA repair mechanisms are essential for tumorigenesis and disruption of HR mechanism is an important predisposing factor of human breast cancers (BC). PALB2 is an important part of the HR. There are similarities between canine mammary tumours (CMT) and BCs. As its human counterpart, PALB2 mutations could be a predisposing factor of CMT.

OBJECTIVES

In this study, we aimed to investigate the impacts of PALB2 variants on tumorigenesis and canine mammary tumor (CMT) malignancy.

METHODS

We performed Sanger sequencing to detect germline mutations in the WD40 domain of the canine PALB2 gene in CMT patients. We conducted in silico analysis to investigate the variants, and compared the germline PALB2 mutations in humans that cause breast cancer (BC) with the variants detected in dogs with CMT.

RESULTS

We identified an intronic (c.3096+8C>G) variant, two exonic (p.A1050V and p.R1354R) variants, and a 3' UTR variant (c.4071T>C). Of these, p.R1354R and c.4071T>C novel variants were identified for the first time in this study. We found that the p.A1050V mutation had a significant effect. However, we could not determine sufficient similarity due to the differences in nucleotide/amino acid sequences between two species. Nonetheless, possible variants of human sequences in the exact location as their dog counterparts are associated with several cancer types, implying that the variants could be crucial for tumorigenesis in dogs. Our results did not show any effect of the variants on tumor malignancy.

CONCLUSIONS

The current project is the first study investigating the relationship between the PALB2 gene WD40 domain and CMTs. Our findings will contribute to a better understanding of the pathogenic mechanism of the PALB2 gene in CMTs. In humans, variant positions in canines have been linked to cancer-related phenotypes such as familial BC, endometrial tumor, and hereditary cancer predisposition syndrome. The results of bioinformatics analyses should be investigated through functional tests or case-control studies.

Roles for epithelial integrin α3β1 in regulation of the microenvironment during normal and pathological tissue remodeling

Integrin receptors for the extracellular matrix activate intracellular signaling pathways that are critical for tissue development, homeostasis, and regeneration/repair, and their loss or dysregulation contributes to many developmental defects and tissue pathologies. This review will focus on tissue remodeling roles for integrin α3β1, a receptor for laminins found in the basement membranes (BMs) that underlie epithelial cell layers. As a paradigm, we will discuss literature that supports a role for α3β1 in promoting ability of epidermal keratinocytes to modify their tissue microenvironment during skin development, wound healing, or tumorigenesis. Preclinical and clinical studies have shown that this role depends largely on ability of α3β1 to govern the keratinocyte's repertoire of secreted proteins, or the "secretome," including 1) matrix proteins and proteases involved in matrix remodeling and 2) paracrine-acting growth factors/cytokines that stimulate other cells with important tissue remodeling functions (e.g., endothelial cells, fibroblasts, inflammatory cells). Moreover, α3β1 signaling controls gene expression that helps epithelial cells carry out these functions, including genes that encode secreted matrix proteins, proteases, growth factors, or cytokines. We will review what is known about α3β1-dependent gene regulation through both transcription and posttranscriptional mRNA stability. Regarding the latter, we will discuss examples of α3β1-dependent alternative splicing (AS) or alternative polyadenylation (APA) that prevents inclusion of cis-acting mRNA sequences that would otherwise target the transcript for degradation via nonsense-mediated decay or destabilizing AU-rich elements (AREs) in the 3'-untranslated region (3'-UTR). Finally, we will discuss prospects and anticipated challenges of exploiting α3β1 as a clinical target for the treatment of cancer or wound healing.

PABPN1 loss-of-function causes APA-shift in oculopharyngeal muscular dystrophy

Alternative polyadenylation (APA) at the 3' UTR of transcripts contributes to the cell transcriptome. APA is suppressed by the nuclear RNA-binding protein PABPN1. Aging-associated reduced PABPN1 levels in skeletal muscles lead to muscle wasting. Muscle weakness in oculopharyngeal muscular dystrophy (OPMD) is caused by short alanine expansion in PABPN1 exon1. The expanded PABPN1 forms nuclear aggregates, an OPMD hallmark. Whether the expanded PABPN1 affects APA and how it contributes to muscle pathology is unresolved. To investigate these questions, we developed a procedure including RNA library preparation and a simple pipeline calculating the APA-shift ratio as a readout for PABPN1 activity. Comparing APA-shift results to previously published PAS utilization and APA-shift results, we validated this procedure. The procedure was then applied on the OPMD cell model and on RNA from OPMD muscles. APA-shift was genome-wide in the mouse OPMD model, primarily affecting muscle transcripts. In OPMD individuals, APA-shift was enriched with muscle transcripts. In an OPMD cell model APA-shift was not significant. APA-shift correlated with reduced expression levels of a subset of PABPN1 isoforms, whereas the expression of the expanded PABPN1 did not correlate with APA-shift. PABPN1 activity is not affected by the expression of expanded PABPN1, but rather by reduced PABPN1 expression levels. In muscles, PABPN1 activity initially affects muscle transcripts. We suggest that muscle weakness in OPMD is caused by PABPN1 loss-of-function leading to APA-shift that primarily affects in muscle transcripts.

An AluYa5 Insertion in the 3'UTR of COL4A1 and Cerebral Small Vessel Disease

Importance

Cerebral small vessel diseases (CSVDs) account for one-fifth of stroke cases. Numerous familial cases remain unresolved after routine screening of known CSVD genes.

Objective

To identify novel genes and mechanisms associated with familial CSVD.

Design, Setting, and Participants

This 2-stage study involved linkage analysis and a case-control study; linkage analysis and whole exome and genome sequencing were used to identify candidate gene variants in 2 large families with CSVD (9 patients with CSVD). Then, a case-control analysis was conducted on 246 unrelated probands, including probands from these 2 families and 244 additional probands. All probands (clinical onset

Main Outcomes and Measures

A pathogenic AluYa5 insertion was identified within the COL4A1 3'UTR in the 2 large families with CSVD. Reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR), Western blot, and long-read RNA sequencing were used to investigate outcomes associated with the insertion using patient fibroblasts. Clinical and magnetic resonance imaging features of probands with variants and available relatives were assessed.

Results

Among 246 probands (141 females [57.3%]; median [IQR] age at referral, 56 [49-64] years), 7 patients of French ancestry carried the insertion. This insertion was absent in 467 healthy French individuals in a control group (odds ratio, ∞; 95% CI, 2.78 to ∞; P = 5 × 10-4) and 10 847 individuals from the gnomAD structural variant database (odds ratio, ∞; 95% CI, 64.77 to ∞; P = 2.42 × 10-12). In these 7 patients' families, 19 family members with CSVD carried the insertion. RT-qPCR and Western blot showed an upregulation of COL4A1 mRNA (10.6-fold increase; 95% CI, 1.4-fold to 17.1-fold increase) and protein levels (2.8-fold increase; 95% CI, 2.1-fold to 3.5-fold increase) in patient vs control group fibroblasts. Long-read RNA sequencing data showed that the insertion was associated with perturbation in the use of canonical COL4A1 polyadenylation signals (approximately 87% of isoforms transcribed from the wild type allele vs 5% of isoforms transcribed from the allele with the insertion used the 2 distal canonical polyadenylation signals). The main clinical feature of individuals with CSVD was the recurrence of pontine ischemic lesions starting at an early age (17 of 19 patients [89.5%]).

Conclusions and relevance

This study found a novel mechanism associated with COL4A1 upregulation and a highly penetrant adult-onset CSVD. These findings suggest that quantitative alterations of the cerebrovascular matrisome are associated with CSVD pathogenesis, with diagnostic and therapeutic implications.

Collapse

Key Words Collapse

MESH Headings

• Adult
• Female
• Humans
• Middle Aged
• 3' Untranslated Regions/genetics
• Alleles
• Case-Control Studies
• Cerebral Small Vessel Diseases/genetics
• Collagen Type IV/metabolism
• Protein Isoforms
• Mutagenesis, Insertional

Collapse

Grants Collapse

Affiliation(s)

Chaker Aloui NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France
Lisa Neumann NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France
Françoise Bergametti NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France
Eric Sartori Service de Neurologie, Centre Hospitalier Bretagne Sud, Lorient, France
Marc Herbreteau Service de Neurologie, Centre Hospitalier Bretagne Sud, Lorient, France
Arnaud Maillard Assistance Publique-Hôpitaux de Paris, Service de Génétique Moléculaire Neurovasculaire, Hôpital Saint-Louis, Paris, France
Thibault Coste NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France Assistance Publique-Hôpitaux de Paris, Service de Génétique Moléculaire Neurovasculaire, Hôpital Saint-Louis, Paris, France
Hélène Morel NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France Assistance Publique-Hôpitaux de Paris, Service de Génétique Moléculaire Neurovasculaire, Hôpital Saint-Louis, Paris, France
Dominique Hervé NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France Assistance Publique-Hôpitaux de Paris, Service de Neurologie, Hôpital Lariboisière, Paris, France
Hugues Chabriat NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France Assistance Publique-Hôpitaux de Paris, Service de Neurologie, Hôpital Lariboisière, Paris, France
Serge Timsit Service de Neurologie Vasculaire, Centre Hospitalier Régional Universitaire de Brest, Brest, France
Irina Viakhireva Service de Neurologie Vasculaire, Centre Hospitalier Régional Universitaire de Brest, Brest, France
Yves Denoyer Service de Neurologie, Centre Hospitalier Bretagne Sud, Lorient, France Université de Rennes, Laboratoire Traitement du Signal et de l'Image, Institut National de la Santé Et de la Recherche Médicale Unité Mixte de Recherche 1099, Rennes, France
Rémi Allibert Service de Neurologie, Unité Neurovasculaire, Centre Hospitalier Universitaire de Saint Etienne, Saint Etienne, France
Florence Demurger Service de Neurologie, Unité Neurovasculaire, Centre Hospitalier Bretagne Atlantique, Vannes, France
Cedric Gollion Service de Neurologie, Centre Hospitalier Universitaire de Toulouse, Toulouse, France
Patrick Vermersch Univ. Lille, Institut National de la Santé Et de la Recherche Médicale Unité Mixte de Recherche 1172 LilNCog, Centre Hospitalier Universitaire Lille, Fédérations Hospitalo-Universitaire Precise, Lille, France
Florence Marchelli Assistance Publique-Hôpitaux de Paris, Service de Génétique Moléculaire Neurovasculaire, Hôpital Saint-Louis, Paris, France
Corinne Blugeon GenomiqueENS, Institut de Biologie de l’Ecole Normale Supérieur, Département de biologie, École Normale Supérieure, Centre National de la Recherche Scientifique, Institut National de la Santé Et de la Recherche Médicale, Université Paris Sciences et Lettres, Paris, France
Sophie Lemoine GenomiqueENS, Institut de Biologie de l’Ecole Normale Supérieur, Département de biologie, École Normale Supérieure, Centre National de la Recherche Scientifique, Institut National de la Santé Et de la Recherche Médicale, Université Paris Sciences et Lettres, Paris, France
Claire Tourtier-Bellosta Service de Neurologie, Centre Hospitalier Bretagne Sud, Lorient, France
Alexis Brouazin Service de neurologie, Centre Hospitalier de Cornouaille, Quimper, France
Anne-Louise Leutenegger NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France
Eva Pipiras NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France Assistance Publique-Hôpitaux de Paris, Hôpitaux Jean Verdier et Armand Trousseau, Université Sorbonne Paris Nord, Bobigny, France
Elisabeth Tournier-Lasserve NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France Assistance Publique-Hôpitaux de Paris, Service de Génétique Moléculaire Neurovasculaire, Hôpital Saint-Louis, Paris, France

Collapse

Star-PAP controls oncogene expression through primary miRNA 3'-end formation to regulate cellular proliferation and tumour formation

Star-PAP is a non-canonical poly(A) polymerase that is down regulated in breast cancer. While Star-PAP down regulation impairs target mRNA polyadenylation, paradoxically, we see up regulation of a large number of oncogenes on Star-PAP knockdown. Using two breast cancer cells (MCF7 with high Star-PAP, and MDA-MB-231 with negligible Star-PAP level), we discover that Star-PAP negatively regulates oncogene expression and subsequently cellular proliferation. This regulation is compromised with Star-PAP mutant of 3'-end processing function (serine 6 to alanine, S6A phospho-mutation). Concomitantly, xenograft mice model using MDA-MB-231 cells reveals a reduction in the tumour formation on ectopic Star-PAP expression that is ameliorated by S6A mutation. We find that Star-PAP control of target oncogene expression is independent of Star-PAP-mediated alternative polyadenylation or target mRNA 3'-end formation. We demonstrate that Star-PAP regulates target oncogenes through cellular miRNAs (miR-421, miR-335, miR-424, miR-543, miR-205, miR-34a, and miR-26a) that are down regulated in breast cancer. Analysis of various steps in miRNA biogenesis pathway reveals that Star-PAP regulates 3'-end formation and synthesis of primary miRNA (host) transcripts that is dependent on S6 phosphorylation thus controlling mature miRNA generation. Using mimics and inhibitors of two target miRNAs (miR-421 and miR-424) after Star-PAP depletion in MCF7 or ectopic expression in MDA-MB-231 cells, we demonstrate that Star-PAP controls oncogene expression and cellular proliferation through targeting miRNAs that regulates tumour formation. Our study establishes a novel mechanism of oncogene expression independent of alternative polyadenylation through Star-PAP-mediated miRNA host transcript polyadenylation that regulates breast cancer progression.

Alternative polyadenylation quantitative trait methylation mapping in human cancers provides clues into the molecular mechanisms of APA

Genetic variants are involved in the orchestration of alternative polyadenylation (APA) events, while the role of DNA methylation in regulating APA remains unclear. We generated a comprehensive atlas of APA quantitative trait methylation sites (apaQTMs) across 21 different types of cancer (1,612 to 60,219 acting in cis and 4,448 to 142,349 in trans). Potential causal apaQTMs in non-cancer samples were also identified. Mechanistically, we observed a strong enrichment of cis-apaQTMs near polyadenylation sites (PASs) and both cis- and trans-apaQTMs in proximity to transcription factor (TF) binding regions. Through the integration of ChIP-signals and RNA-seq data from cell lines, we have identified several regulators of APA events, acting either directly or indirectly, implicating novel functions of some important genes, such as TCF7L2, which is known for its involvement in type 2 diabetes and cancers. Furthermore, we have identified a vast number of QTMs that share the same putative causal CpG sites with five different cancer types, underscoring the roles of QTMs, including apaQTMs, in the process of tumorigenesis. DNA methylation is extensively involved in the regulation of APA events in human cancers. In an attempt to elucidate the potential underlying molecular mechanisms of APA by DNA methylation, our study paves the way for subsequent experimental validations into the intricate biological functions of DNA methylation in APA regulation and the pathogenesis of human cancers. To present a comprehensive catalog of apaQTM patterns, we introduce the Pancan-apaQTM database, available at https://pancan-apaqtm-zju.shinyapps.io/pancanaQTM/.

The Transcriptional Landscape of Immune-Response 3'-UTR Alternative Polyadenylation in Melanoma

The prognosis of patients with malignant melanoma has been improved in recent decades due to advancements in immunotherapy. However, a considerable proportion of patients are refractory to treatment, particularly at advanced stages. This underscores the necessity of developing a new strategy to improve it. Alternative polyadenylation (APA), as a marker of crucial posttranscriptional regulation, has emerged as a major new type of epigenetic marker involved in tumorigenesis. However, the potential roles of APA in shaping the tumor microenvironment (TME) are largely unexplored. Herein, we collected two cohorts comprising melanoma patients who received immune checkpoint inhibitor (ICI) immunotherapy to quantify transcriptome-wide discrepancies in APA. We observed a global change in 3'-UTRs between responders and non-responders, which might involve DNA damage response, angiogenesis, PI3K-AKT signaling pathways, etc. Ten putative master APA regulatory factors for those APA events were detected via a network analysis. Notably, we established an immune response-related APA scoring system (IRAPAss), which exhibited a great performance of predicting immunotherapy response in multiple cohorts. Furthermore, we examined the correlation of APA with TME at the single-cell level using four single-cell immune profiles of tumor-infiltrating lymphocytes (TILs), which revealed an overall discrepancy in 3'-UTR length across diverse T cell populations, probably contributing to immunoregulation in melanoma. In conclusion, our study provides a transcriptional landscape of APA implicated in immunoregulation, which might lay the foundation for developing a new strategy for improving immunotherapy response for melanoma patients by targeting APA.

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

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

CPSF3 inhibition blocks pancreatic cancer cell proliferation through disruption of core histone mRNA processing

Pancreatic ductal adenocarcinoma (PDAC) is a lethal disease with limited effective treatment options, potentiating the importance of uncovering novel drug targets. Here, we target cleavage and polyadenylation specificity factor 3 (CPSF3), the 3' endonuclease that catalyzes mRNA cleavage during polyadenylation and histone mRNA processing. We find that CPSF3 is highly expressed in PDAC and is associated with poor prognosis. CPSF3 knockdown blocks PDAC cell proliferation and colony formation in vitro and tumor growth in vivo. Chemical inhibition of CPSF3 by the small molecule JTE-607 also attenuates PDAC cell proliferation and colony formation, while it has no effect on cell proliferation of nontransformed immortalized control pancreatic cells. Mechanistically, JTE-607 induces transcriptional readthrough in replication-dependent histones, reduces core histone expression, destabilizes chromatin structure, and arrests cells in the S-phase of the cell cycle. Therefore, CPSF3 represents a potential therapeutic target for the treatment of PDAC.

Genome-wide kinetic profiling of pre-mRNA 3' end cleavage

Cleavage and polyadenylation is necessary for the formation of mature mRNA molecules. The rate at which this process occurs can determine the temporal availability of mRNA for subsequent function throughout the cell and is likely tightly regulated. Despite advances in high-throughput approaches for global kinetic profiling of RNA maturation, genome-wide 3' end cleavage rates have never been measured. Here, we describe a novel approach to estimate the rates of cleavage, using metabolic labeling of nascent RNA, high-throughput sequencing, and mathematical modeling. Using in silico simulations of nascent RNA-seq data, we show that our approach can accurately and precisely estimate cleavage half-lives for both constitutive and alternative sites. We find that 3' end cleavage is fast on average, with half-lives under a minute, but highly variable across individual sites. Rapid cleavage is promoted by the presence of canonical sequence elements and an increased density of polyadenylation signals near a cleavage site. Finally, we find that cleavage rates are associated with the localization of RNA polymerase II at the end of a gene, and faster cleavage leads to quicker degradation of downstream readthrough RNA. Our findings shed light on the features important for efficient 3' end cleavage and the regulation of transcription termination.

The implication of integrative multiple RNA modification-based subtypes in gastric cancer immunotherapy and prognosis

Previous studies have focused on the impact of individual RNA modifications on tumor development. This study comprehensively investigated the effects of multiple RNA modifications, including m6A, alternative polyadenylation, pseudouridine, adenosine-to-inosine editing, and uridylation, on gastric cancer (GC). By analyzing 1,946 GC samples from eleven independent cohorts, we identified distinct clusters of RNA modification genes with varying survival rates and immunological characteristics. We assessed the chromatin activity of these RNA modification clusters through regulon enrichment analysis. A prognostic model was developed using Stepwise Regression and Random Survival Forest algorithms and validated in ten independent datasets. Notably, the low-risk group showed a more favorable prognosis and positive response to immune checkpoint blockade therapy. Single-cell RNA sequencing confirmed the abundant expression of signature genes in B cells and plasma cells. Overall, our findings shed light on the potential significance of multiple RNA modifications in GC prognosis, stemness development, and chemotherapy resistance.

The role of m6A epigenetic modifications in tumor coding and non-coding RNA processing

BACKGROUND

Epigenetic modifications of RNA significantly contribute to the regulatory processes in tumors and have, thus, received considerable attention. The m6A modification, known as N6-methyladenosine, is the predominant epigenetic alteration found in both eukaryotic mRNAs and ncRNAs.

MAIN BODY

m6A methylation modifications are dynamically reversible and are catalyzed, removed, and recognized by the complex of m6A methyltransferase (MTases), m6A demethylase, and m6A methyl recognition proteins (MRPs). Published evidence suggests that dysregulated m6A modification results in abnormal biological behavior of mature mRNA, leading to a variety of abnormal physiological processes, with profound implications for tumor development in particular.

CONCLUSION

Abnormal RNA processing due to dysregulation of m6A modification plays an important role in tumor pathogenesis and potential mechanisms of action. In this review, we comprehensively explored the mechanisms by which m6A modification regulates mRNA and ncRNA processing, focusing on their roles in tumors, and aiming to understand the important regulatory function of m6A modification, a key RNA epigenetic modification, in tumor cells, with a view to providing theoretical support for tumor diagnosis and treatment. Video Abstract.

The rosetta stone of successful ageing: does oral health have a role?

Ageing is an inevitable aspect of life and thus successful ageing is an important focus of recent scientific efforts. The biological process of ageing is mediated through the interaction of genes with environmental factors, increasing the body's susceptibility to insults. Elucidating this process will increase our ability to prevent and treat age-related disease and consequently extend life expectancy. Notably, centenarians offer a unique perspective on the phenomenon of ageing. Current research highlights several age-associated alterations on the genetic, epigenetic and proteomic level. Consequently, nutrient sensing and mitochondrial function are altered, resulting in inflammation and exhaustion of regenerative ability.Oral health, an important contributor to overall health, remains underexplored in the context of extreme longevity. Good masticatory function ensures sufficient nutrient uptake, reducing morbidity and mortality in old age. The relationship between periodontal disease and systemic inflammatory pathologies is well established. Diabetes, rheumatoid arthritis and cardiovascular disease are among the most significant disease burdens influenced by inflammatory oral health conditions. Evidence suggests that the interaction is bi-directional, impacting progression, severity and mortality. Current models of ageing and longevity neglect an important factor in overall health and well-being, a gap that this review intends to illustrate and inspire avenues for future research.

PABPN1 promotes clear cell renal cell carcinoma progression by suppressing the alternative polyadenylation of SGPL1 and CREG1

Alternative polyadenylation (APA) is an important post-transcriptional regulatory mechanism in cancer development and progression. Poly(A) binding protein nuclear 1 (PABPN1) is a gene that encodes abundant nuclear protein, binds with high affinity to nascent poly(A) tails, and is crucial for 3'-UTR (3'-untranslated region) APA. Although PABPN1 has been recently reported as a dominant master APA regulator in clear cell renal cell carcinoma (ccRCC), the underlying functional mechanism remain unclear and the genes subject to PABPN1 regulation that contribute to ccRCC progression have not been identified. Here, we found that PABPN1 is upregulated in ccRCC, and its expression is highly associated with the clinical prognosis of ccRCC patients. PABPN1 promotes ccRCC cell proliferation, migration, invasion, and exerts an influence on sphingolipid metabolism and cell cycle. Moreover, PABPN1 depletion significantly suppressed cancer cell growth via induction of cell cycle arrest and apoptosis. In particular, we characterized PABPN1-regulated 3'-UTR APA of sphingosine-1-phosphate lyase 1 (SGPL1) and cellular repressor of E1A stimulated genes 1 (CREG1), which contribute to ccRCC progression. Collectively, our data revealed that PABPN1 promotes ccRCC progression at least in part, by suppressing SGPL1 and CREG1. Thus, PABPN1 may be a potential therapeutic target in ccRCC.

Pharmacological and therapeutic potentials of cordycepin in hematological malignancies

Hematological malignancies(HMs) are highly heterogeneous diseases with globally rising incidence. Despite major improvements in the management of HMs, conventional therapies have limited efficacy, and relapses with high mortality rates are still frequent. Cordycepin, a nucleoside analog extracted from Cordyceps species, represents a wide range of therapeutic effects, including anti-inflammatory, anti-tumor, and anti-metastatic activities. Cordycepin induces apoptosis in different subtypes of HMs by triggering adenosine receptors, death receptors, and several vital signaling pathways such as MAPK, ERK, PI3K, AKT, and GSK-3β/β-catenin. This review article summarizes the impact of utilizing cordycepin on HMs, and highlights its potential as a promising avenue for future cancer research based on evidence from in vitro and in vivo studies, as well as clinical trials.

Alternative polyadenylation reprogramming of MORC2 induced by NUDT21 loss promotes KIRC carcinogenesis

Alternative polyadenylation (APA), a posttranscriptional mechanism of gene expression via determination of 3'UTR length, has an emerging role in carcinogenesis. Although abundant APA reprogramming is found in kidney renal clear cell carcinoma (KIRC), which is one of the major malignancies, whether APA functions in KIRC remains unknown. Herein, we found that chromatin modifier MORC2 gained oncogenic potential in KIRC among the genes with APA reprogramming, and moreover, its oncogenic potential was enhanced by 3'UTR shortening through stabilization of MORC2 mRNA. MORC2 was found to function in KIRC by downregulating tumor suppressor DAPK1 via DNA methylation. Mechanistically, MORC2 recruited DNMT3A to facilitate hypermethylation of the DAPK1 promoter, which was strengthened by 3'UTR shortening of MORC2. Furthermore, loss of APA regulator NUDT21, which was induced by DNMT3B-mediated promoter methylation, was identified as responsible for 3'UTR shortening of MORC2 in KIRC. Additionally, NUDT21 was confirmed to act as a tumor suppressor mainly depending on downregulation of MORC2. Finally, we designed an antisense oligonucleotide (ASO) to enhance NUDT21 expression and validated its antitumor effect in vivo and in vitro. This study uncovers the DNMT3B/NUDT21/APA/MORC2/DAPK1 regulatory axis in KIRC, disclosing the role of APA in KIRC and the crosstalk between DNA methylation and APA.

Cell-type-specific alternative polyadenylation promotes oncogenic gene expression in non-small cell lung cancer progression

Disrupted alternative polyadenylation (APA) is frequently involved in tumorigenesis and cancer progression by regulating the gene expression of oncogenes and tumor suppressors. However, limited knowledge of tumor-type- and cell-type-specific APA events may lead to novel APA events and their functions being overlooked. Here, we compared APA events across different cell types in non-small cell lung cancer (NSCLC) and normal tissues and identified functionally related APA events in NSCLC. We found several cell-specific 3'-UTR alterations that regulate gene expression changes showed prognostic value in NSCLC. We further investigated the function of APA-mediated 3'-UTR shortening through loss of microRNA (miRNA)-binding sites, and we identified and experimentally validated several oncogene-miRNA-tumor suppressor axes. According to our analyses, we found SPARC as an APA-regulated oncogene in cancer-associated fibroblasts in NSCLC. Knockdown of SPARC attenuates lung cancer cell invasion and metastasis. Moreover, we found high SPARC expression associated with resistance to several drugs except cisplatin. NSCLC patients with high SPARC expression could benefit more compared to low-SPARC-expression patients with cisplatin treatment. Overall, our comprehensive analysis of cell-specific APA events shed light on the regulatory mechanism of cell-specific oncogenes and provided opportunities for combination of APA-regulated therapeutic target and cell-specific therapy development.

Alternative polyadenylation regulation in cardiac development and cardiovascular disease

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

Lengthening of 3' Untranslated Regions of mRNAs by Alternative Polyadenylation Is Associated With Tumor Progression and Poor Prognosis of Clear Cell Renal Cell Carcinoma

Alternative polyadenylation (APA) is emerging as a major posttranscriptional mechanism for gene regulation in cancer. A prevailing hypothesis is that shortening of the 3' untranslated region (3'UTR) increases oncoprotein expression because of the loss of miRNA-binding sites (MBSs). We showed that the longer 3'UTR is associated with a more advanced tumor stage in patients with clear cell renal cell carcinoma (ccRCC). More surprisingly, 3'UTR shortening is correlated with better overall survival in patients with ccRCC. Furthermore, we identified a mechanism by which longer transcripts lead to increased oncogenic protein and decreased tumor-suppressive protein expression compared to the shorter transcripts. In our model, shortening of 3'UTRs by APA may increase the mRNA stability of the majority of the potential tumor-suppressor genes due to the loss of MBSs and AU-rich elements (AREs). Unlike potential tumor-suppressor genes, the potential oncogenes display much lower MBS and ARE density and globally much higher m6A density in distal 3'UTRs. As a result, 3'UTRs shortening decreases the mRNA stability of potential oncogenes and enhances the mRNA stability of potential tumor-suppressor genes. Our findings highlight the cancer-specific pattern of APA regulation and extend our understanding of the mechanism of APA-mediated 3'UTR length changes in cancer biology.

Splicing Characterization and Isoform Switch Events in Human Keratinocytes Carrying Oncogenes from High-Risk HPV-16 and Low-Risk HPV-84

Infection of epithelial cells with high-risk HPV (HR-HPV) types, followed by expression of virus oncogenic proteins (E5, E6, and E7), leads to genomic imbalance, suppression of tumor inhibitors, and induction of oncogenes. Low-risk HPV (LR-HPV) may slow the rate at which cervical cancer spreads to an invasive stage since co-infection with LR-HPV is linked to a decreased risk of future invasive cancer than infection with HR-HPV alone. We then propose that cancer-progressing changes may be distinguished through identifying the functional differences between LR-HPV and HR-HPV. Lentiviral strategies were followed to establish HaCaT cells with constitutive expression of HPV oncogenes. RNAseq experiments were designed to analyze the transcriptome modulations caused by each of the E5, E6, and E7 oncogenes of HPV-16 and HPV-84 in HaCaT cells. We identified enhanced RNA degradation, spliceosome, and RNA polymerase pathways related to mRNA processing. ATTS (alternative transcription termination site) was discovered to be more prevalent in cells with HPV-16E5 than HPV-84E5. In HPV-16E6-infected cells, ATTS gain was significantly higher than ATTS loss. Cells with HPV-16E7 had more isoforms with intron retention (IR) than those with HPV-84E7. We identified switches in ADAM10, CLSPN, and RNPS1 that led to greater expression of the coding isoforms in HR-HPV. The results of this work highlight differences between LR-HPV and HR-HPV in mRNA processing. Moreover, crucial cervical cancer-related switch events were detected.

Regulation of the hypertonic stress response by the 3' mRNA cleavage and polyadenylation complex

Maintenance of osmotic homeostasis is one of the most aggressively defended homeostatic set points in physiology. One major mechanism of osmotic homeostasis involves the upregulation of proteins that catalyze the accumulation of solutes called organic osmolytes. To better understand how osmolyte accumulation proteins are regulated, we conducted a forward genetic screen in Caenorhabditis elegans for mutants with no induction of osmolyte biosynthesis gene expression (Nio mutants). The nio-3 mutant encoded a missense mutation in cpf-2/CstF64, while the nio-7 mutant encoded a missense mutation in symk-1/Symplekin. Both cpf-2 and symk-1 are nuclear components of the highly conserved 3' mRNA cleavage and polyadenylation complex. cpf-2 and symk-1 block the hypertonic induction of gpdh-1 and other osmotically induced mRNAs, suggesting they act at the transcriptional level. We generated a functional auxin-inducible degron (AID) allele for symk-1 and found that acute, post-developmental degradation in the intestine and hypodermis was sufficient to cause the Nio phenotype. symk-1 and cpf-2 exhibit genetic interactions that strongly suggest they function through alterations in 3' mRNA cleavage and/or alternative polyadenylation. Consistent with this hypothesis, we find that inhibition of several other components of the mRNA cleavage complex also cause a Nio phenotype. cpf-2 and symk-1 specifically affect the osmotic stress response since heat shock-induced upregulation of a hsp-16.2::GFP reporter is normal in these mutants. Our data suggest a model in which alternative polyadenylation of 1 or more mRNAs is essential to regulate the hypertonic stress response.

AR-A014418 regulates intronic polyadenylation and transcription of PD-L1 through inhibiting CDK12 and CDK13 in tumor cells

BACKGROUND

Immune checkpoint molecules, especially programmed death 1 (PD-1) and its ligand, programmed death ligand 1 (PD-L1), protect tumor cells from T cell-mediated killing. Immune checkpoint inhibitors, designed to restore the antitumor immunosurveillance, have exhibited significant clinical benefits for patients with certain cancer types. Nevertheless, the relatively low response rate and acquisition of resistance greatly limit their clinical applications. A deeper understanding of the regulatory mechanisms of PD-L1 protein expression and activity will help to develop more effective therapeutic strategies.

METHODS

The effects of AR-A014418 and THZ531 on PD-L1 expression were detected by western blot, reverse transcription-quantitative polymerase chain reaction (RT-qPCR) and flow cytometry. In vitro kinase assays with recombinant proteins were performed to confirm that AR-A014418 functioned as a CDK12 and CDK13 dual inhibitor. The roles of CDK12 and CDK13 in intronic polyadenylation (IPA) and transcription of PD-L1 were determined via RNA interference or protein overexpression. T-cell cytotoxicity assays were used to validate the activation of antitumor immunity by AR-A014418 and THZ531.

RESULTS

AR-A014418 inhibits CDK12 to enhance the IPA, and inhibits CDK13 to repress the transcription of PD-L1. IPA generates a secreted PD-L1 isoform (PD-L1-v4). The extent of IPA was not enough to reduce full-length PD-L1 expression obviously. Only the superposition of enhancing IPA and repressing transcription (dual inhibition of CDK12 and CDK13) dramatically suppresses full-length PD-L1 induction by interferon-γ. AR-A014418 and THZ531 could potentiate T-cell cytotoxicity against tumor cells.

CONCLUSIONS

Our work identifies a new regulatory pathway for PD-L1 expression and discovers CDK12 and CDK13 as promising drug targets for immune modulation and combined therapeutic strategies.

Genome-wide screening for genetic variants in polyadenylation signal (PAS) sites in mouse selection lines for fatness and leanness

Alternative polyadenylation (APA) determines mRNA stability, localisation, translation and protein function. Several diseases, including obesity, have been linked to APA. Studies have shown that single nucleotide polymorphisms in polyadenylation signals (PAS-SNPs) can influence APA and affect phenotype and disease susceptibility. However, these studies focussed on associations between single PAS-SNP alleles with very large effects and phenotype. Therefore, we performed a genome-wide screening for PAS-SNPs in the polygenic mouse selection lines for fatness and leanness by whole-genome sequencing. The genetic variants identified in the two lines were overlapped with locations of PAS sites obtained from the PolyASite 2.0 database. Expression data for selected genes were extracted from the microarray expression experiment performed on multiple tissue samples. In total, 682 PAS-SNPs were identified within 583 genes involved in various biological processes, including transport, protein modifications and degradation, cell adhesion and immune response. Moreover, 63 of the 583 orthologous genes in human have been previously associated with human diseases, such as nervous system and physical disorders, and immune, endocrine, and metabolic diseases. In both lines, PAS-SNPs have also been identified in genes broadly involved in APA, such as Polr2c, Eif3e and Ints11. Five PAS-SNPs within 5 genes (Car, Col4a1, Itga7, Lat, Nmnat1) were prioritised as potential functional variants and could contribute to the phenotypic disparity between the two selection lines. The developed PAS-SNPs catalogue presents a key resource for planning functional studies to uncover the role of PAS-SNPs in APA, disease susceptibility and fat deposition.

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.

Circular RNAs and Untranslated Regions in Acute Myeloid Leukemia

Before the advent of next-generation sequencing, research on acute myeloid leukemia (AML) mostly centered on protein-coding genes. In recent years, breakthroughs in RNA sequencing technologies and whole transcriptome analysis have led to the discovery that approximately 97.5% of the human genome is transcribed into non-coding RNAs (ncRNAs). This paradigm shift has led to an explosion of research interest in different classes of non-coding RNAs, such as circular RNAs (circRNAs) as well as non-coding untranslated regions (UTRs) of protein-coding messenger RNAs. The critical roles of circRNAs and UTRs in AML pathogenesis have become increasingly apparent. In this review, we discuss the cellular mechanisms of circRNAs and summarize recent studies that reveal their biological roles in AML. Furthermore, we also review the contribution of 3'UTRs to disease progression. Finally, we discuss the potential of circRNAs and 3'UTRs as new biomarkers for disease stratification and/or the prediction of treatment response and targets for the development of RNA-directed therapeutic applications.

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.

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.

Advances of CRISPR-Cas13 system in COVID-19 diagnosis and treatment

The ongoing global pandemic of coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has resulted in over 570 million infections and 6 million deaths worldwide. Early detection and quarantine are essential to arrest the spread of the highly contagious COVID-19. High-risk groups, such as older adults and individuals with comorbidities, can present severe symptoms, including pyrexia, pertussis, and acute respiratory distress syndrome, on SARS-CoV-2 infection that can prove fatal, demonstrating a clear need for high-throughput and sensitive platforms to detect and eliminate SARS-CoV-2. CRISPR-Cas13, an emerging CRISPR system targeting RNA with high specificity and efficiency, has recently drawn much attention for COVID-19 diagnosis and treatment. Here, we summarized the current research progress on CRISPR-Cas13 in COVID-19 diagnosis and treatment and highlight the challenges and future research directions of CRISPR-Cas13 for effectively counteracting COVID-19.

UBE3D Regulates mRNA 3'-End Processing and Maintains Adipogenic Potential in 3T3-L1 Cells

We have previously described the role of an essential Saccharomyces cerevisiae gene, important for cleavage and polyadenylation 1 (IPA1), in the regulation of gene expression through its interaction with Ysh1, the endonuclease subunit of the mRNA 3'-end processing complex. Through a similar mechanism, the mammalian homolog ubiquitin protein ligase E3D (UBE3D) promotes the migratory and invasive potential of breast cancer cells, but its role in the regulation of gene expression during normal cellular differentiation has not previously been described. In this study, we show that CRISPR/Cas9-mediated knockout of Ube3d in 3T3-L1 cells blocks their ability to differentiate into mature adipocytes. Consistent with previous studies in other cell types, Ube3d knockout leads to decreased levels of CPSF73 and global changes in cellular mRNAs indicative of a loss of 3'-end processing capacity. Ube3d knockout cells also display decreased expression of known preadipogenic markers. Overexpression of either UBE3D or CPSF73 rescues the differentiation defect and partially restores protein levels of these markers. These results support a model in which UBE3D is necessary for the maintenance of the adipocyte-committed state via its regulation of the mRNA 3'-end processing machinery.

Nuclear m6 A reader YTHDC1 suppresses proximal alternative polyadenylation sites by interfering with the 3' processing machinery

N6-methyladenosine (m6 A) and alternative polyadenylation (APA) are important regulators of gene expression in eukaryotes. Recently, it was found that m6 A is closely related to APA. However, the molecular mechanism of this new APA regulation remains elusive. Here, we show that YTHDC1, a nuclear m6 A reader, can suppress proximal APA sites and produce longer 3' UTR transcripts by binding to their upstream m6 A sites. YTHDC1 can directly interact with the 3' end processing factor FIP1L1 and interfere with its ability to recruit CPSF4. Binding to the m6 A sites can promote liquid-liquid phase separation of YTHDC1 and FIP1L1, which may play an important role in their interaction and APA regulation. Collectively, YTHDC1 as an m6 A "reader" links m6 A modification with pre-mRNA 3' end processing, providing a new mechanism for APA regulation.

A Novel apaQTL-SNP for the Modification of Non-Small-Cell Lung Cancer Susceptibility across Histological Subtypes

Background: Alternative polyadenylation (APA) events may be modulated by single nucleotide polymorphisms (SNPs). Therefore, this study aims to evaluate the association between APA quantitative trait loci (apaQTLs)-related SNPs (apaQTL-SNPs) and non-small-cell lung cancer (NSCLC) risk. Methods: APA-related genes associated with NSCLC (LUAD and LUSC) were first identified, and the respective apaQTL-SNPs of those genes were selected. Then, a two-phase case-control study was performed to evaluate the association between candidate apaQTL-SNPs and NSCLC risk. Results: A total of 7 LUAD- and 21 LUSC-associated apaQTL-SNPs were selected. In the first phase, the apaQTL-SNP rs10138506 was significantly associated with LUAD risk (p < 0.05), whereas the other two apaQTL-SNPs (rs1130698 and rs1130719) were significantly associated with LUSC risk (p < 0.05). In the second phase, the variant G allele of rs10138506 was still significantly associated with an increased risk of LUAD (OR = 1.42, 95%CI = 1.02−1.98, p = 0.038). Functional annotation indicated that the variant G allele of rs10138506 was significantly associated with a higher PDUI value of CHURC1. Meanwhile, 3′RACE experiments verified the presence of two poly(A) sites (proximal and distal) in CHURC1, while qRT-PCR results indicated that different genotypes of rs1127968 which, in perfect LD with rs10138506, can mediate changes in the lengths of the 3′UTR of CHURC1 isoforms. Conclusion: The variant G allele of rs10138506 in CHURC1 was correlated with a longer 3′UTR of CHURC1 mRNA and an increased LUAD risk. Further studies should evaluate the interaction between rs10138506 and different 3′UTR lengths of CHURC1 that regulate LUAD development.

HSV-1 and influenza infection induce linear and circular splicing of the long NEAT1 isoform

The herpes simplex virus 1 (HSV-1) virion host shut-off (vhs) protein cleaves both cellular and viral mRNAs by a translation-initiation-dependent mechanism, which should spare circular RNAs (circRNAs). Here, we show that vhs-mediated degradation of linear mRNAs leads to an enrichment of circRNAs relative to linear mRNAs during HSV-1 infection. This was also observed in influenza A virus (IAV) infection, likely due to degradation of linear host mRNAs mediated by the IAV PA-X protein and cap-snatching RNA-dependent RNA polymerase. For most circRNAs, enrichment was not due to increased circRNA synthesis but due to a general loss of linear RNAs. In contrast, biogenesis of a circRNA originating from the long isoform (NEAT1_2) of the nuclear paraspeckle assembly transcript 1 (NEAT1) was induced both in HSV-1 infection-in a vhs-independent manner-and in IAV infection. This was associated with induction of novel linear splicing of NEAT1_2 both within and downstream of the circRNA. NEAT1_2 forms a scaffold for paraspeckles, nuclear bodies located in the interchromatin space, must likely remain unspliced for paraspeckle assembly and is up-regulated in HSV-1 and IAV infection. We show that NEAT1_2 splicing and up-regulation can be induced by ectopic co-expression of the HSV-1 immediate-early proteins ICP22 and ICP27, potentially linking increased expression and splicing of NEAT1_2. To identify other conditions with NEAT1_2 splicing, we performed a large-scale screen of published RNA-seq data. This uncovered both induction of NEAT1_2 splicing and poly(A) read-through similar to HSV-1 and IAV infection in cancer cells upon inhibition or knockdown of CDK7 or the MED1 subunit of the Mediator complex phosphorylated by CDK7. In summary, our study reveals induction of novel circular and linear NEAT1_2 splicing isoforms as a common characteristic of HSV-1 and IAV infection and highlights a potential role of CDK7 in HSV-1 or IAV infection.

New hallmarks of ageing: a 2022 Copenhagen ageing meeting summary

Genomic instability, telomere attrition, epigenetic alterations, mitochondrial dysfunction, loss of proteostasis, deregulated nutrient-sensing, cellular senescence, stem cell exhaustion, and altered intercellular communication were the original nine hallmarks of ageing proposed by López-Otín and colleagues in 2013. The proposal of these hallmarks of ageing has been instrumental in guiding and pushing forward research on the biology of ageing. In the nearly past 10 years, our in-depth exploration on ageing research has enabled us to formulate new hallmarks of ageing which are compromised autophagy, microbiome disturbance, altered mechanical properties, splicing dysregulation, and inflammation, among other emerging ones. Amalgamation of the 'old' and 'new' hallmarks of ageing may provide a more comprehensive explanation of ageing and age-related diseases, shedding light on interventional and therapeutic studies to achieve healthy, happy, and productive lives in the elderly.

PCF11, a Novel CD44-Downstream Transcriptional Target, Linking Its 3'-End Polyadenylation Function to Tumor Cell Metastasis

Breast Cancer (BC) is the most common and the major health issue in women worldwide. Metastasis, a multistep process, is the worst aspect of cancer and tumor cell invasion is the defining step. Tumor cell invasion requires cell adhesion molecules (CAMs), and alterations in CAMs is considered as an initiating event in metastasis. Among CAMs, CD44 is a large family of more than 100 isoform, and its precise function was initially controversial in BC. Therefore, we have previously established a (Tet)-off inducible expression system of CD44 in MCF-7 primary BC cell line, and showed that CD44 promoted BC invasion/metastasis both in vitro and in vivo. A microarray gene expression profiling revealed more than 200 CD44-downstream potential transcriptional target genes, mediating its role in BC cell invasion and metastasis. Among these CD44-target genes, the Pre-mRNA cleavage complex 2 protein (PCF11) was upregulated upon the activation of CD44 by its major ligand hyaluronan (HA); This prompted us to hypothesize PCF11 as a potential novel transcriptional target of CD44-promoted BC cell invasion and metastasis. A large body of evidence from the literature supports our hypothesis that CD44 might regulate PCF11 via MAPK/ERK pathway. This review aims to discuss these findings from the literature that support our hypothesis, and further provide possible mechanisms linking CD44-promoted cell invasion through regulation of its potential target PCF11.

Chicken CDS2 isoforms presented distinct spatio-temporal expression pattern and regulated by insulin in a breed-specific manner

The cytidine diphosphate diacylglycerol synthases (CDSs) gene encodes the cytidine diphosphate-diacylglycerol (CDP-DAG) synthase enzyme that catalyzes the formation of CDP-diacylglycerol from phosphatidic acid. At present, there are no reports of CDS2 in birds. Here, we identified chicken CDS2 transcripts by combining conventional RT-PCR amplification, 5′ rapid amplification of cDNA ends (RACE), and 3′ RACE, explored the spatio-temporal expression profiles of total CDS2 and the longest transcript variant CDS2-4, and investigated the effect of exogenous insulin on the mRNA level of total CDS2 via quantitative RT-PCR. Four transcripts of chicken CDS2 (CDS2-1, -2, -3, and -4) were identified, which were alternatively spliced at the 3′-untranslated region (UTR). Both total CDS2 and CDS2-4 were prominently expressed in adipose tissue, and exhibited low expression in liver and pectoralis of 49-day-old chickens. Regarding the spatio-temporal expression patterns of CDS2 in chicken, total CDS2 exhibited a similar temporal expression tendency with a high level in the later period of incubation (embryonic day 19 [E19] or 1-day-old) in the brain, liver, and pectoralis. While CDS2-4 presented a distinct temporal expression pattern in these tissues, CDS2-4 levels peaked at 21 d in the brain and pectoralis, while liver CDS2-4 mRNA levels were highest at the early stage of hatching (E10). Total CDS2 (P < 0.001) and CDS2-4 (P = 0.0090) mRNA levels in the liver were differentially regulated throughout the development of the chicken. Total CDS2 levels in the liver of Silky chickens were higher than that of the broiler in the basal state and after insulin stimulation. Exogenous insulin significantly down-regulated the level of total CDS2 at 240 min in the pectoralis of Silky chickens (P < 0.01). In conclusion, chicken CDS2 isoforms with variation at the 3′-UTR were identified, which was prominently expressed in adipose tissue. Total CDS2 and CDS2-4 presented distinct spatio-temporal expression patterns, that is they were differentially regulated with age in brain, liver, and pectoralis. Insulin could regulate chicken CDS2 levels in a breed- and tissue-specific manner.

RNA polyadenylation patterns in the human transcriptome

The eukaryotic transcriptome undergoes various post-transcriptional modifications which assists gene expression. Polyadenylation is a molecular process occurring at the 3'-end of the RNA molecule which involves the poly(A) polymerase attaching adenine monophosphate molecules in a chain-like fashion to assemble a poly(A) tail. Multiple RNA isoforms are produced with differing 3'-UTR and exonic compositions through alternative polyadenylation (APA) which enhances the diversification of alternatively spliced mRNA transcripts. To study polyadenylation patterns, novel methods have been developed using short-read and long-read sequencing technologies to analyse the 3'-ends of the transcript. Recent studies have identified unique polyadenylation patterns in different cellular functions, including oncogenic activity, which could prove valuable in the understanding of medical genetics, particularly in the discovery of biomarkers in diseased states. We present a review of current literature reporting on polyadenylation and the biological relevance in the mammalian transcriptome, with a focus on the human transcriptome. Additionally, we have explored the various methods available to detect polyadenylation patterns using second and third generation sequencing technologies.

The role of RNA binding proteins in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the leading causes of overall cancer deaths worldwide with limited therapeutic options. Due to the heterogeneity of HCC pathogenesis, the molecular mechanisms underlying HCC development are not fully understood. Emerging evidence indicates that RNA-binding proteins (RBPs) play a vital role throughout hepatocarcinogenesis. Thus, a deeper understanding of how RBPs contribute to HCC progression will provide new tools for early diagnosis and prognosis of this devastating disease. In this review, we summarize the tumor suppressive and oncogenic roles of RBPs and their roles in hepatocarcinogenesis. The diagnostic and therapeutic potential of RBPs in HCC, including their limitations, are also discussed.

The fellowship of the RING: BRCA1, its partner BARD1 and their liaison in DNA repair and cancer

The breast cancer type 1 susceptibility protein (BRCA1) and its partner - the BRCA1-associated RING domain protein 1 (BARD1) - are key players in a plethora of fundamental biological functions including, among others, DNA repair, replication fork protection, cell cycle progression, telomere maintenance, chromatin remodeling, apoptosis and tumor suppression. However, mutations in their encoding genes transform them into dangerous threats, and substantially increase the risk of developing cancer and other malignancies during the lifetime of the affected individuals. Understanding how BRCA1 and BARD1 perform their biological activities therefore not only provides a powerful mean to prevent such fatal occurrences but can also pave the way to the development of new targeted therapeutics. Thus, through this review work we aim at presenting the major efforts focused on the functional characterization and structural insights of BRCA1 and BARD1, per se and in combination with all their principal mediators and regulators, and on the multifaceted roles these proteins play in the maintenance of human genome integrity.

The Identification of Alternative Polyadenylation in Stomach Adenocarcinomas Using the Genotype-Tissue Expression Project and the Cancer Genome Atlas- Stomach Adenocarcinoma Profiles

Objective

Alternative polyadenylation (APA) is a common mechanism that is present in most human genes and determines the length of the messenger ribonucleic acid (mRNA) three prime untranslated region (3ʹ-UTR), which can give rise to changes in mRNA stability and localization. However, little is known about the specific changes related to APA in stomach adenocarcinomas (STADs).

Methods

We integrated RNA sequencing data from The Cancer Genome Atlas and Genotype-Tissue Expression project to comprehensively analyze APA events in 289 cases of STAD.

Results

Our results showed that APA events were widespread in patients with STAD and were rich in genes related to known STAD pathways. The APA events result in the loss of tumor-suppressing micro-ribonucleic acid (miRNA) binding sites and increased heterogeneity in the length of the 3ʹ-UTR altered genes. Survival analysis revealed that specific subsets of 3ʹ-UTR-altered genes independently characterized a poor prognostic cohort among patients with STAD, thereby indicating the potential of APA as a new prognostic biomarker.

Conclusion

Our single-cancer analysis showed that by losing miRNA regulation, APA can become a driving factor for regulating the expression of oncogenic genes in STAD and promote its development. Our research revealed that APA events regulated STAD genes that were functionally related, thereby providing a new approach for gaining a better understanding of the progress of STADs and a means for identifying new drug targets as avenues of treatment.

Emerging roles of alternative cleavage and polyadenylation (APA) in human disease

In the messenger RNA (mRNA) maturation process, the 3'-end of pre-mRNA is cleaved and a poly(A) sequence is added, this is an important determinant of mRNA stability and its cellular functions. More than 60%-70% of human genes have three or more polyadenylation (APA) sites and can be cleaved at different sites, generating mRNA transcripts of varying lengths. This phenomenon is termed as alternative cleavage and polyadenylation (APA) and it plays role in key biological processes like gene regulation, cell proliferation, senescence, and also in various human diseases. Loss of regulatory microRNA binding sites and interactions with RNA-binding proteins leading to APA are largely investigated in human diseases. However, the functions of the core APA machinery and related factors during disease conditions remain largely unknown. In this review, we discuss the roles of polyadenylation machinery in relation to brain disease, cardiac failure, pulmonary fibrosis, cancer, infectious conditions, and other human diseases. Collectively, we believe this review will be a useful avenue for understanding the emerging role of APA in the pathobiology of various human diseases.

LABRAT reveals association of alternative polyadenylation with transcript localization, RNA binding protein expression, transcription speed, and cancer survival

BACKGROUND

The sequence content of the 3' UTRs of many mRNA transcripts is regulated through alternative polyadenylation (APA). The study of this process using RNAseq data, though, has been historically challenging.

RESULTS

To combat this problem, we developed LABRAT, an APA isoform quantification method. LABRAT takes advantage of newly developed transcriptome quantification techniques to accurately determine relative APA site usage and how it varies across conditions. Using LABRAT, we found consistent relationships between gene-distal APA and subcellular RNA localization in multiple cell types. We also observed connections between transcription speed and APA site choice as well as tumor-specific transcriptome-wide shifts in APA isoform abundance in hundreds of patient-derived tumor samples that were associated with patient prognosis. We investigated the effects of APA on transcript expression and found a weak overall relationship, although many individual genes showed strong correlations between relative APA isoform abundance and overall gene expression. We interrogated the roles of 191 RNA-binding proteins in the regulation of APA isoforms, finding that dozens promote broad, directional shifts in relative APA isoform abundance both in vitro and in patient-derived samples. Finally, we find that APA site shifts in the two classes of APA, tandem UTRs and alternative last exons, are strongly correlated across many contexts, suggesting that they are coregulated.

CONCLUSIONS

We conclude that LABRAT has the ability to accurately quantify APA isoform ratios from RNAseq data across a variety of sample types. Further, LABRAT is able to derive biologically meaningful insights that connect APA isoform regulation to cellular and molecular phenotypes.