Differential genome-wide profiling of tandem 3' UTRs among human breast cancer and normal cells by high-throughput sequencing

SLC27A5 inhibits cancer stem cells by inducing alternative polyadenylation of METTL14 in hepatocellular carcinoma

Solute carrier family 27 member 5 (SLC27A5/FATP5) is a liver-specific metabolic enzyme that plays a crucial role in fatty acid transport and bile acid metabolism. Deficiency of SLC27A5 promotes the progression of hepatocellular carcinoma (HCC) and is strongly associated with a poor prognosis. SLC27A5 exhibits noncanonical functions beyond its metabolic role; however, its specific mechanisms in hepatocarcinogenesis remain elusive and are therefore investigated in this study. Immunoprecipitation-mass spectrometry analysis showed that SLC27A5-interacting proteins were significantly enriched in alternative polyadenylation (APA). RNA-sequencing data provided evidence that SLC27A5 plays a global role in regulating APA events in HCC. Mechanistically, SLC27A5 facilitates the usage of the proximal polyadenylation site of METTL14 by downregulating the expression of the APA-associated factor PABPC1, resulting in the shortening of the METTL14-3'UTR and the conversion of METTL14-UL to METTL14-US. In contrast to METTL14-UL, METTL14-US escapes the inhibitory effect of miRNA targeting, leading to increased METTL14 expression. METTL14-US upregulation by SLC27A5 suppressed the stemness of HCC. Therefore, low levels of SLC27A5 and METTL14 may serve as reliable biomarkers for identifying poor prognosis in patients with HCC. In conclusion, SLC27A5/PABPC1 inhibits HCC stemness via APA-regulated expression of METTL14, providing potential avenues for the development of novel therapeutic strategies.

Transcriptomic profiling of severe and critical COVID-19 patients reveals alterations in expression, splicing and polyadenylation

Coronavirus disease 2019 (COVID-19) is a multi-systemic illness that became a pandemic in March 2020. Although environmental factors and comorbidities can influence disease progression, there is a lack of prognostic markers to predict the severity of COVID-19 illness. Identifying these markers is crucial for improving patient outcomes and appropriately allocating scarce resources. Here, an RNA-sequencing study was conducted on blood samples from unvaccinated, hospitalized patients divided by disease severity; 367 moderate, 173 severe, and 199 critical. Using a bioinformatics approach, we identified differentially expressed genes (DEGs), alternative splicing (AS) and alternative polyadenylation (APA) events that were severity-dependent. In the severe group, we observed a higher expression of kappa immunoglobulins compared to the moderate group. In the critical cohort, a majority of AS events were mutually exclusive exons and APA genes mostly had longer 3'UTRs. Interestingly, multiple genes associated with cytoskeleton, TUBA4A, NRGN, BSG, and CD300A, were differentially expressed, alternatively spliced and polyadenylated in the critical group. Furthermore, several inflammation-related pathways were observed predominantly in critical vs. moderate. We demonstrate that integrating multiple downstream analyses of transcriptomics, from moderate, severe, and critical patients confers a significant advantage in identifying relevant dysregulated genes and pathways.

E3 ligase SYVN1-mediated polyubiquitination of CPSF6 promotes alternative polyadenylation and antivirus effects of macrophages

Transcriptome-wide alternative polyadenylation (APA) is involved in both innate and adaptive immune responses of immune cells. Downregulation of the CPSF6 protein, one of the 3' end-processing factors, mediates APA in macrophages with responses to virus infection and plays an important role in its anti-virus effect. However, the signaling pathway and molecular mechanism underlying the downregulation of the CPSF6 protein remain elusive. Here, we found that MAVS triggers the nuclear import of the E3 ligase SYVN1 mediated by NUP153 in response to vesicular stomatitis virus infection. Then, SYVN1 catalyzes K48-linked polyubiquitination of CPSF6, resulting in degradation of CPSF6 via the proteasome and then transcriptome-wide APA and anti-virus effects. Our results identify an antiviral mechanism via APA regulation based on ubiquitination modification of the CPSF6 protein, which may serve as a target for developing immune interventions.

CPSF1 inhibition promotes widespread use of intergenic polyadenylation sites and impairs glycolysis in prostate cancer cells

Localized prostate cancer can be cured by radiation or surgery, but advanced prostate cancer continues to be a clinical challenge. Altered alternative polyadenylation occurs in numerous cancers and can downregulate tumor-suppressor genes and upregulate oncogenes. We found that the cleavage and polyadenylation specificity factor (CPSF) complex factor CPSF1 is upregulated in patients with advanced prostate cancer, with high CPSF1 expression correlating with worse progression-free survival. Knockdown of CPSF1 selectively inhibited the growth of prostate cancer cells and reduced glycolytic output. Evaluating the changes in global poly(A) site usage in prostate cancer cells following CPSF1 knockdown revealed widespread usage of intergenic poly(A) sites distal to annotated 3' UTRs, which lengthened 3' UTRs and decreased levels of thousands of mRNAs, including key glycolysis genes. These findings uncover a role for CPSF1 in the suppression of intergenic poly(A) sites in prostate cancer and nominate CPSF1 as a therapeutic target in advanced prostate cancer.

CircSugp1 interacts with CPSF6 to modulate intestinal mucosa repair by regulating alternative polyadenylation-mediated shortening of the Wdr89 3'UTR

Circular RNAs are a single-stranded non-coding RNAs and play an important role in the development of many diseases. Alternative polyadenylation (APA) regulates the gene 3'UTR length for controlling gene expressions. Although the APA mechanism has been widely studied in the development of diseases, there is no data on its role in the burned intestinal mucosa. We thus herein assessed the role of the circSugp1-initiating APA mechanism in the burned intestinal mucosa. CircSugp1 was downregulated in the intestinal mucosa of burned mice. CircSugp1 promoted proliferation and migration in vitro and in vivo. CircSugp1 promotes the expression of CPSF6; the overexpression of CPSF6 can shorten the gene 3'UTR within the transcript APA range. The promoting effect of circSugp1 on value-added migration was mediated by the APA regulation of the Wdr89 short 3'UTR isoform. CircSugp1 targeted the upregulation of the expression of CPSF6, followed by upregulation of the expression of Wdr89 through APA, promoting the repair of intestinal mucosal damage in burned mice.

Coupling of alternative splicing and alternative polyadenylation

RNA splicing and 3'-cleavage and polyadenylation (CPA) are essential processes for the maturation of RNA. There have been extensive independent studies of these regulated processing events, including alternative splicing (AS) and alternative polyadenylation (APA). However, growing evidence suggests potential crosstalk between splicing and 3'-end processing in regulating AS or APA. Here, we first provide a brief overview of the molecular machines involved in splicing and 3'-end processing events, and then review recent studies on the functions and mechanisms of the crosstalk between the two processes. On the one hand, 3'-end processing can affect splicing, as 3'-end processing factors and CPA-generated polyA tail promote the splicing of the last intron. Beyond that, 3'-end processing factors can also influence the splicing of internal and terminal exons. Those 3'-end processing factors can also interact with different RNA-binding proteins (RBPs) to exert their effects on AS. The length of 3' untranslated region (3' UTR) can affect the splicing of upstream exons. On the other hand, splicing and CPA may compete within introns in generating different products. Furthermore, splicing within the 3' UTR is a significant factor contributing to 3' UTR diversity. Splicing also influences 3'-end processing through the actions of certain splicing factors. Interestingly, some classical RBPs play dual roles in both splicing and 3'-end processing. Finally, we discuss how long-read sequencing technologies aid in understanding the coordination of AS-APA events and envision that these findings may potentially promote the development of new strategies for disease diagnosis and treatment.

Regulation of alternative polyadenylation isoforms of Timp2 is an effector event of RAS signaling in cell transformation

Alternative polyadenylation (APA) generates mRNA isoforms with different lengths of the 3' untranslated region (3' UTR). The tissue inhibitor of metalloproteinase 2 (TIMP2) plays a key role in extracellular matrix remodeling under various developmental and disease conditions. Both human and mouse genes encoding TIMP2 contain two highly conserved 3'UTR APA sites, leading to mRNA isoforms that differ substantially in 3'UTR size. APA of Timp2 is one of the most significantly regulated events in multiple cell differentiation lineages. Here we show that Timp2 APA is highly regulated in transformation of NIH3T3 cells by the oncogene HRAS G12V . Perturbations of isoform expression with long 3'UTR isoform-specific knockdown or genomic removal of the alternative UTR (aUTR) region indicate that the long 3'UTR isoform contributes to the secreted Timp2 protein much more than the short 3'UTR isoform. The short and long 3'UTR isoforms differ in subcellular localization to endoplasmic reticulum (ER). Strikingly, Timp2 aUTR enhances secreted protein expression but no effect on intracellular proteins in reporter assays. Furthermore, downregulation of Timp2 long isoform mitigates gene expression changes elicited by HRAS G12V . Together, our data indicate that regulation of Timp2 protein expression through APA isoform changes is an integral part of RAS-mediated cell transformation and 3'UTR isoforms of Timp2 can have distinct impacts on expression of secreted vs. intracellular proteins.

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.

Disruption of PABPN1 phase separation by SNRPD2 drives colorectal cancer cell proliferation and migration through promoting alternative polyadenylation of CTNNBIP1

Generally shortened 3' UTR due to alternative polyadenylation (APA) is widely observed in cancer, but its regulation mechanisms for cancer are not well characterized. Here, with profiling of APA in colorectal cancer tissues and poly(A) signal editing, we firstly identified that the shortened 3' UTR of CTNNIBP1 in colorectal cancer promotes cell proliferation and migration. We found that liquid-liquid phase separation (LLPS) of PABPN1 is reduced albeit with higher expression in cancer, and the reduction of LLPS leads to the shortened 3' UTR of CTNNBIP1 and promotes cell proliferation and migration. Notably, the splicing factor SNRPD2 upregulated in colorectal cancer, can interact with glutamic-proline (EP) domain of PABPN1, and then disrupt LLPS of PABPN1, which attenuates the repression effect of PABPN1 on the proximal poly(A) sites. Our results firstly reveal a new regulation mechanism of APA by disruption of LLPS of PABPN1, suggesting that regulation of APA by interfering LLPS of 3' end processing factor may have the potential as a new way for the treatment of cancer.

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 cleavage and polyadenylation of the Ccnb1 mRNA defines accumulation of cyclin protein during the meiotic cell cycle

Progression through the mitotic and meiotic cell cycle is driven by fluctuations in the levels of cyclins, the regulatory subunits controlling the localization and activity of CDK1 kinases. Cyclin levels are regulated through a precise balance of synthesis and degradation. Here we demonstrate that the synthesis of Cyclin B1 during the oocyte meiotic cell cycle is defined by the selective translation of mRNA variants generated through alternative cleavage and polyadenylation (APA). Using gene editing in mice, we introduced mutations into the proximal and distal polyadenylation elements of the 3' untranslated region (UTR) of the Ccnb1 mRNA. Through in vivo loss-of-function experiments, we demonstrate that the translation of mRNA with a short 3' UTR specifies Cyclin B1 protein levels that set the timing of meiotic re-entry. In contrast, translation directed by a long 3' UTR is necessary to direct Cyclin B1 protein accumulation during the MI/MII transition. These findings establish that the progression through the cell cycle is dependent on the selective translation of multiple mRNA variants generated by APA.

Comparative alternative polyadenylation profiles in differentiated adipocytes of subcutaneous and intramuscular fat tissue in cattle

Alternative polyadenylation (APA) is a key molecular mechanism involved in the post-transcriptional regulation of gene expression, which has been proven to play a critical role in cell differentiation. In the present study, we performed IVT-SAPAS sequencing to profile the dynamic changes of APA sites in bovine subcutaneous preadipocytes and intramuscular preadipocytes during adipogenesis. A total of 52621 high quality APA sites were identified in preadipocytes and adipocytes. Compared with preadipocytes, the increased usage of canonical AATAAA was observed in the cell-biased APA sites of adipocytes. Furthermore, 1933 and 2140 differentially expressed APA (DE-APA) sites, as well as 341 and 337 untranslated region-APA (UTR-APA) switching genes were identified in subcutaneous preadipocytes and intramuscular preadipocytes during adipogenesis, respectively. The UTR-APA switching genes showed divergent trends in preadipocytes, among which UTR-APA switching genes in intramuscular preadipocytes tended to use shorter 3'UTR for differentiation into mature adipocytes. APA events mediated by UTR-APA switching in intramuscular adipocytes were enriched in lipid synthesis and adipocyte differentiation. TRIB3, WWTR1, and INSIG1 played important roles in the differentiation of intramuscular preadipocytes. Briefly, our results provided new insights into understanding the mechanisms of bovine adipocyte differentiation.

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

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

CPSF6 regulates alternative polyadenylation and proliferation of cancer cells through phase separation

Cancer cells usually exhibit shortened 3' untranslated regions (UTRs) due to alternative polyadenylation (APA) to promote cell proliferation and migration. Upregulated CPSF6 leads to a systematic prolongation of 3' UTRs, but CPSF6 expression in tumors is typically higher than that in healthy tissues. This contradictory observation suggests that it is necessary to investigate the underlying mechanism by which CPSF6 regulates APA switching in cancer. Here, we find that CPSF6 can undergo liquid-liquid phase separation (LLPS), and elevated LLPS is associated with the preferential usage of the distal poly(A) sites. CLK2, a kinase upregulated in cancer cells, destructs CPSF6 LLPS by phosphorylating its arginine/serine-like domain. The reduction of CPSF6 LLPS can lead to a shortened 3' UTR of cell-cycle-related genes and accelerate cell proliferation. These results suggest that CPSF6 LLPS, rather than its expression level, may be responsible for APA regulation in cancer cells.

Downregulation of apoptotic repressor AVEN exacerbates cardiac injury after myocardial infarction

Myocardial infarction (MI) is a leading cause of heart failure (HF), associated with morbidity and mortality worldwide. As an essential part of gene expression regulation, the role of alternative polyadenylation (APA) in post-MI HF remains elusive. Here, we revealed a global, APA-mediated, 3' untranslated region (3' UTR)-lengthening pattern in both human and murine post-MI HF samples. Furthermore, the 3' UTR of apoptotic repressor gene, AVEN, is lengthened after MI, contributing to its downregulation. AVEN knockdown increased cardiomyocyte apoptosis, whereas restoration of AVEN expression substantially improved cardiac function. Mechanistically, AVEN 3' UTR lengthening provides additional binding sites for miR-30b-5p and miR-30c-5p, thus reducing AVEN expression. Additionally, PABPN1 (poly(A)-binding protein 1) was identified as a potential regulator of AVEN 3' UTR lengthening after MI. Altogether, our findings revealed APA as a unique mechanism regulating cardiac injury in response to MI and also indicated that the APA-regulated gene, AVEN, holds great potential as a critical therapeutic target for treating post-MI HF.

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.

Alternative 3'-untranslated regions regulate high-salt tolerance of Spartina alterniflora

High-salt stress continues to challenge the growth and survival of many plants. Alternative polyadenylation (APA) produces mRNAs with different 3'-untranslated regions (3' UTRs) to regulate gene expression at the post-transcriptional level. However, the roles of alternative 3' UTRs in response to salt stress remain elusive. Here, we report the function of alternative 3' UTRs in response to high-salt stress in S. alterniflora (Spartina alterniflora), a monocotyledonous halophyte tolerant of high-salt environments. We found that high-salt stress induced global APA dynamics, and ∼42% of APA genes responded to salt stress. High-salt stress led to 3' UTR lengthening of 207 transcripts through increasing the usage of distal poly(A) sites. Transcripts with alternative 3' UTRs were mainly enriched in salt stress-related ion transporters. Alternative 3' UTRs of HIGH-AFFINITY K+ TRANSPORTER 1 (SaHKT1) increased RNA stability and protein synthesis in vivo. Regulatory AU-rich elements were identified in alternative 3' UTRs, boosting the protein level of SaHKT1. RNAi-knock-down experiments revealed that the biogenesis of 3' UTR lengthening in SaHKT1 was controlled by the poly(A) factor CLEAVAGE AND POLYADENYLATION SPECIFICITY FACTOR 30 (SaCPSF30). Over-expression of SaHKT1 with an alternative 3' UTR in rice (Oryza sativa) protoplasts increased mRNA accumulation of salt-tolerance genes in an AU-rich element-dependent manner. These results suggest that mRNA 3' UTR lengthening is a potential mechanism in response to high-salt stress. These results also reveal complex regulatory roles of alternative 3' UTRs coupling APA and regulatory elements at the post-transcriptional level in plants.

Relationship between NUDT21 mediated alternative polyadenylation process and tumor

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

U1 snRNP proteins promote proximal alternative polyadenylation sites by directly interacting with 3' end processing core factors

In eukaryotic cells, both alternative splicing and alternative polyadenylation (APA) play essential roles in the gene regulation network. U1 small ribonucleoprotein particle (U1 snRNP) is a major component of spliceosome, and U1 snRNP complex can suppress proximal APA sites through crosstalking with 3' end processing factors. However, here we show that both knockdown and overexpression of SNRPA, SNRPC, SNRNP70, and SNRPD2, the U1 snRNP proteins, promote the usage of proximal APA sites at the transcriptome level. SNRNP70 can drive the phase transition of PABPN1 from droplet to aggregate, which may reduce the repressive effects of PABPN1 on the proximal APA sites. Additionally, SNRNP70 can also promote the proximal APA sites by recruiting CPSF6, suggesting that the function of CPSF6 on APA is related with other RNA-binding proteins and cell context-dependent. Consequently, these results reveal that, on the contrary to U1 snRNP complex, the free proteins of U1 snRNP complex can promote proximal APA sites through the interaction with 3' end processing machinery.

PABPC1--mRNA stability, protein translation and tumorigenesis

Mammalian poly A-binding proteins (PABPs) are highly conserved multifunctional RNA-binding proteins primarily involved in the regulation of mRNA translation and stability, of which PABPC1 is considered a central regulator of cytoplasmic mRNA homing and is involved in a wide range of physiological and pathological processes by regulating almost every aspect of RNA metabolism. Alterations in its expression and function disrupt intra-tissue homeostasis and contribute to the development of various tumors. There is increasing evidence that PABPC1 is aberrantly expressed in a variety of tumor tissues and cancers such as lung, gastric, breast, liver, and esophageal cancers, and PABPC1 might be used as a potential biomarker for tumor diagnosis, treatment, and clinical application in the future. In this paper, we review the abnormal expression, functional role, and molecular mechanism of PABPC1 in tumorigenesis and provide directions for further understanding the regulatory role of PABPC1 in tumor cells.

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.

CYCLIN K down-regulation induces androgen receptor gene intronic polyadenylation, variant expression and PARP inhibitor vulnerability in castration-resistant prostate cancer

Expression of androgen receptor variants (AR-Vs) is implicated in the development of castration-resistant prostate cancer (PCa). Others have shown that androgen depletion or antiandrogen treatment induces AR-V expression in PCa cell lines, xenografts, and patient samples, although the underlying mechanism remains unclear. Our findings reveal that hormonal therapy–induced CYCLIN K down-regulation represents a key mechanism that drives intronic polyadenylation (IPA) usage in the AR gene and AR-V expression and castration resistance in PCa, and that this mechanism of action can be therapeutically targeted by the PARP inhibitor.

Androgen receptor (AR) messenger RNA (mRNA) alternative splicing variants (AR-Vs) are implicated in castration-resistant progression of prostate cancer (PCa), although the molecular mechanism underlying the genesis of AR-Vs remains poorly understood. The CDK12 gene is often deleted or mutated in PCa and CDK12 deficiency is known to cause homologous recombination repair gene alteration or BRCAness via alternative polyadenylation (APA). Here, we demonstrate that pharmacological inhibition or genetic inactivation of CDK12 induces AR gene intronic (intron 3) polyadenylation (IPA) usage, AR-V expression, and PCa cell resistance to the antiandrogen enzalutamide (ENZ). We further show that AR binds to the CCNK gene promoter and up-regulates CYCLIN K expression. In contrast, ENZ decreases AR occupancy at the CCNK gene promoter and suppresses CYCLIN K expression. Similar to the effect of the CDK12 inhibitor, CYCLIN K degrader or ENZ treatment promotes AR gene IPA usage, AR-V expression, and ENZ-resistant growth of PCa cells. Importantly, we show that targeting BRCAness induced by CYCLIN K down-regulation with the PARP inhibitor overcomes ENZ resistance. Our findings identify CYCLIN K down-regulation as a key driver of IPA usage, hormonal therapy–induced AR-V expression, and castration resistance in PCa. These results suggest that hormonal therapy–induced AR-V expression and therapy resistance are vulnerable to PARP inhibitor treatment.

Developmentally regulated alternate 3' end cleavage of nascent transcripts controls dynamic changes in protein expression in an adult stem cell lineage

Alternative polyadenylation (APA) generates transcript isoforms that differ in the position of the 3' cleavage site, resulting in the production of mRNA isoforms with different length 3' UTRs. Although widespread, the role of APA in the biology of cells, tissues, and organisms has been controversial. We identified >500 Drosophila genes that express mRNA isoforms with a long 3' UTR in proliferating spermatogonia but a short 3' UTR in differentiating spermatocytes due to APA. We show that the stage-specific choice of the 3' end cleavage site can be regulated by the arrangement of a canonical polyadenylation signal (PAS) near the distal cleavage site but a variant or no recognizable PAS near the proximal cleavage site. The emergence of transcripts with shorter 3' UTRs in differentiating cells correlated with changes in expression of the encoded proteins, either from off in spermatogonia to on in spermatocytes or vice versa. Polysome gradient fractionation revealed >250 genes where the long 3' UTR versus short 3' UTR mRNA isoforms migrated differently, consistent with dramatic stage-specific changes in translation state. Thus, the developmentally regulated choice of an alternative site at which to make the 3' end cut that terminates nascent transcripts can profoundly affect the suite of proteins expressed as cells advance through sequential steps in a differentiation lineage.

Dynamic and widespread control of poly(A) tail length during macrophage activation

The poly(A) tail enhances translation and transcript stability, and tail length is under dynamic control during cell state transitions. Tail regulation plays essential roles in translational timing and fertilization in early development, but poly(A) tail dynamics have not been fully explored in post-embryonic systems. Here, we examined the landscape and impact of tail length control during macrophage activation. Upon activation, more than 1500 mRNAs, including proinflammatory genes, underwent distinctive changes in tail lengths. Increases in tail length correlated with mRNA levels regardless of transcriptional activity, and many mRNAs that underwent tail extension encode proteins necessary for immune function and post-transcriptional regulation. Strikingly, we found that ZFP36, whose protein product destabilizes target transcripts, undergoes tail extension. Our analyses indicate that many mRNAs undergoing tail lengthening are, in turn, degraded by elevated levels of ZFP36, constituting a post-transcriptional feedback loop that ensures transient regulation of transcripts integral to macrophage activation. Taken together, this study establishes the complexity, relevance, and widespread nature of poly(A) tail dynamics, and the resulting post-transcriptional regulation during macrophage activation.

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.

TMAO-Activated Hepatocyte-Derived Exosomes Impair Angiogenesis via Repressing CXCR4

Objective: Trimethylamine-N-oxide (TMAO) was found to play crucial roles in vascular endothelial function. However, the exact molecular mechanisms are not yet entirely clear. Recently, we found that exosomes (Exos) isolated from TMAO-treated hepatocytes (TMAO-Exos) contained a distinctive profile of miRNAs compared to those from the TMAO-free group (Control-Exos). Furthermore, TMAO-Exos could notably promote inflammation, damage vascular endothelial cells (VECs), and impair endothelium-dependent vasodilation. This study aimed to further evaluate the effects of TMAO-Exos on VECs and explore the underlying mechanisms.

Methods: Exos were isolated from the hepatocyte culture supernatant with or without TMAO, using differential centrifugation. Then, VECs were treated with these Exos for 48 h and subjected to RNA-sequencing for detecting the changes of alternative polyadenylation (APA) and mRNA. After validation by qPCR and western blotting, the recombinant viruses were used to mediate the overexpression of C-X-C motif chemokine receptor 4 (CXCR4). The in vitro VEC function was evaluated by cell migration and tube formation, and in vivo angiogenesis was investigated in hindlimb ischemia models.

Results: Exos released from hepatocytes were differentially regulated by TMAO; both could be taken up by VECs; and furthermore, TMAO-Exos significantly reduced cell migration and tube formation in vitro and impaired perfusion recovery and angiogenesis after hindlimb ischemia, by down-regulating the CXCR4 expression. However, TMAO-Exos failed to regulate the splicing events, at least in this experimental setting, which suggested that TMAO-Exos may affect CXCR4 expression via an APA-independent manner.

Conclusions: Our findings revealed a novel indirect mechanism by which TMAO impaired endothelial function through stimulating hepatocytes to produce Exos that possessed distinctive activity. The crosstalk between the liver and vascular endothelial mediated by these Exos may offer a new target for restraining the harmful effects induced by TMAO.

Codon optimality in cancer

A key characteristic of cancer cells is their increased proliferative capacity, which requires elevated levels of protein synthesis. The process of protein synthesis involves the translation of codons within the mRNA coding sequence into a string of amino acids to form a polypeptide chain. As most amino acids are encoded by multiple codons, the nucleotide sequence of a coding region can vary dramatically without altering the polypeptide sequence of the encoded protein. Although mutations that do not alter the final amino acid sequence are often thought of as silent/synonymous, these can still have dramatic effects on protein output. Because each codon has a distinct translation elongation rate and can differentially impact mRNA stability, each codon has a different degree of 'optimality' for protein synthesis. Recent data demonstrates that the codon preference of a transcriptome matches the abundance of tRNAs within the cell and that this supply and demand between tRNAs and mRNAs varies between different cell types. The largest observed distinction is between mRNAs encoding proteins associated with proliferation or differentiation. Nevertheless, precisely how codon optimality and tRNA expression levels regulate cell fate decisions and their role in malignancy is not fully understood. This review describes the current mechanistic understanding on codon optimality, its role in malignancy and discusses the potential to target codon optimality therapeutically in the context of cancer.

Suppression of CPSF6 Enhances Apoptosis Through Alternative Polyadenylation-Mediated Shortening of the VHL 3'UTR in Gastric Cancer Cells

Alternative polyadenylation (APA) is an important RNA post-transcriptional process, which can generate diverse mRNA isoforms. Increasing evidence shows that APA is involved in cell self-renewal, development, immunity, and cancer. CPSF6 is one of the core proteins of CFIm complex and can modulate the APA process. Although it has been reported to play oncogenic roles in cancer, the underlying mechanisms remain unclear. The aim of the present study was to characterize CPSF6 in human gastric cancer (GC). We observed that CPSF6 was upregulated in GC. Knockdown of CPSF6 inhibited proliferation and enhanced apoptosis of GC cells both in vitro and in vivo. Global APA site profiling analysis revealed that knockdown of CPSF6 induced widespread 3′UTR shortening of genes in GC cells, including VHL. We also found CPSF6 negatively regulated the expression of VHL through APA and VHL short-3′UTR isoform enhanced apoptosis and inhibited cell growth in GC cells. Our data suggested that CPSF6-induced cell proliferation and inhibition of apoptosis were mediated by the preferential usage of poly(A) in VHL. Our data provide insights into the function of CPSF6 and may imply potential therapeutic targets against GC.

mRNA alternative polyadenylation (APA) in regulation of gene expression and diseases

The mRNA polyadenylation plays essential function in regulation of mRNA metabolism. Mis-regulations of mRNA polyadenylation are frequently linked with aberrant gene expression and disease progression. Under the action of polyadenylate polymerase, poly(A) tail is synthesized after the polyadenylation signal (PAS) sites on the mRNAs. Alternative polyadenylation (APA) often occurs in mRNAs with multiple poly(A) sites, producing different 3' ends for transcript variants, and therefore plays important functions in gene expression regulation. In this review, we first summarize the classical process of mRNA 3'-terminal formation and discuss the length control mechanisms of poly(A) in nucleus and cytoplasm. Then we review the research progress on alternative polyadenylation regulation and the APA site selection mechanism. Finally, we summarize the functional roles of APA in the regulation of gene expression and diseases including cancers.

Alternative polyadenylation: An untapped source for prostate cancer biomarkers and therapeutic targets?

OBJECTIVE

To review alternative polyadenylation (APA) as a mechanism of gene regulation and consider potential roles for APA in prostate cancer (PCa) biology and treatment.

METHODS

An extensive review of mRNA polyadenylation, APA, and PCa literature was performed. This review article introduces APA and its association with human disease, outlines the mechanisms and components of APA, reviews APA in cancer biology, and considers whether APA may contribute to PCa progression and/or produce novel biomarkers and therapeutic targets for PCa.

RESULTS

Eukaryotic mRNA 3'-end cleavage and polyadenylation play a critical role in gene expression. Most human genes encode more than one polyadenylation signal, and produce more than one transcript isoform, through APA. Polyadenylation can occur throughout the gene body to generate transcripts with differing 3'-termini and coding sequence. Differences in 3'-untranslated regions length can modify post-transcriptional gene regulation by microRNAs and RNA binding proteins, and alter mRNA stability, translation efficiency, and subcellular localization. Distinctive APA patterns are associated with human diseases, tissue origins, and changes in cellular proliferation rate and differentiation state. APA events may therefore generate unique mRNA biomarkers or therapeutic targets in certain cancer types or phenotypic states.

CONCLUSIONS

The full extent of cancer-associated and tissue-specific APA events have yet to be defined, and the mechanisms and functional consequences of APA in cancer remain incompletely understood. There is evidence that APA is active in PCa, and that it may be an untapped resource for PCa biomarkers or therapeutic targets.

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.

The landscape of alternative polyadenylation in single cells of the developing mouse embryo

3′ untranslated regions (3′ UTRs) post-transcriptionally regulate mRNA stability, localization, and translation rate. While 3′-UTR isoforms have been globally quantified in limited cell types using bulk measurements, their differential usage among cell types during mammalian development remains poorly characterized. In this study, we examine a dataset comprising ~2 million nuclei spanning E9.5–E13.5 of mouse embryonic development to quantify transcriptome-wide changes in alternative polyadenylation (APA). We observe a global lengthening of 3′ UTRs across embryonic stages in all cell types, although we detect shorter 3′ UTRs in hematopoietic lineages and longer 3′ UTRs in neuronal cell types within each stage. An analysis of RNA-binding protein (RBP) dynamics identifies ELAV-like family members, which are concomitantly induced in neuronal lineages and developmental stages experiencing 3′-UTR lengthening, as putative regulators of APA. By measuring 3′-UTR isoforms in an expansive single cell dataset, our work provides a transcriptome-wide and organism-wide map of the dynamic landscape of alternative polyadenylation during mammalian organogenesis.

Alternative polyadenylation regulates localization, half-life and translation of mRNA isoforms. Here the authors investigate alternative polyadenylation using single cell RNA sequencing data from mouse embryos and identify 3’-UTR isoforms that are regulated across cell types and developmental time.

Alternative polyadenylation events differ dramatically between Tongcheng and Large White pigs in response to PRRSV infection

Alternative polyadenylation (APA) is a widespread post-transcriptional regulation mechanism that increases the biological complexity of transcriptome and proteome. However, it is unclear whether APA regulation plays a role in genetic resistance to porcine reproductive and respiratory syndrome virus (PRRSV). Here, we reported genome-wide APA regulation of porcine alveolar macrophages in PRRSV-resistant Tongcheng (TC) pigs and PRRSV-susceptible Large White (LW) pigs upon PRRSV infection. Using 3' mRNA sequencing strategy, we detected 75 981 high-quality APA sites in porcine alveolar macrophages of TC and LW pigs. Furthermore, 1202 and 1089 differentially expressed APA sites, as well as 79 and 117 untranslated region-APA switching genes were identified in TC pigs and LW pigs upon PRRSV infection respectively. The APA events in TC pigs and LW pigs were involved in different biological pathways, while APA events in TC pigs are directly associated with the immune response to PRRSV infection. In addition, we identified genetic variations affecting polyadenylation signal between TC pigs and LW pigs. These findings would provide helpful information on APA regulation for further understanding of genetic resistance to PRRSV.

Alternative 3' UTRs play a widespread role in translation-independent mRNA association with the endoplasmic reticulum

Transcripts encoding membrane and secreted proteins are known to associate with the endoplasmic reticulum (ER) through translation. Here, using cell fractionation, polysome profiling, and 3' end sequencing, we show that transcripts differ substantially in translation-independent ER association (TiERA). Genes in certain functional groups, such as cell signaling, tend to have significantly higher TiERA potentials than others, suggesting the importance of ER association for their mRNA metabolism, such as localized translation. The TiERA potential of a transcript is determined largely by size, sequence content, and RNA structures. Alternative polyadenylation (APA) isoforms can have distinct TiERA potentials because of changes in transcript features. The widespread 3' UTR lengthening in cell differentiation leads to greater transcript association with the ER, especially for genes that are capable of expressing very long 3' UTRs. Our data also indicate that TiERA is in dynamic competition with translation-dependent ER association, suggesting limited space on the ER for mRNA association.

TRENDseq-A highly multiplexed high throughput RNA 3' end sequencing for mapping alternative polyadenylation

Alternative polyadenylation (APA) is a widespread and highly dynamic mechanism of gene regulation. It affects more than 70% of all genes, resulting in transcript isoforms with distinct 3' end termini. APA thereby considerably expands the diversity of the transcriptome 3' end (TREND). This leads to mRNA isoforms with profoundly different physiological effects, by affecting protein output, production of distinct protein isoforms, or modulating protein localization. APA is globally regulated in various conditions, including developmental and adaptive programs. Since perturbations of APA can disrupt biological processes, ultimately resulting in most devastating disorders, querying the APA landscape is crucial to decipher underlying mechanisms, resulting consequences and potential diagnostic and therapeutic implications. Here we provide a detailed step-by-step protocol for TRENDseq, a method for transcriptome-wide high-throughput sequencing of polyadenylated RNA 3' ends in a highly multiplexed fashion. TRENDseq exploits linear amplification of the starting material to improve sensitivity while significantly reducing the amount of input material. It thereby represents a powerful tool to study APA in numerous experimental set-ups and/or limited human samples in a highly multiplexed and reproducible manner.

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

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

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

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

PAPOLA contributes to cyclin D1 mRNA alternative polyadenylation and promotes breast cancer cell proliferation

Poly(A) polymerases add the poly(A) tail at the 3' end of nearly all eukaryotic mRNA, and are associated with proliferation and cancer. To elucidate the role of the most-studied mammalian poly(A) polymerase, poly(A) polymerase α (PAPOLA), in cancer, we assessed its expression in 221 breast cancer samples and found it to correlate strongly with the aggressive triple-negative subtype. Silencing PAPOLA in MCF-7 and MDA-MB-231 breast cancer cells reduced proliferation and anchorage-independent growth by decreasing steady-state cyclin D1 (CCND1) mRNA and protein levels. Whereas the length of the CCND1 mRNA poly(A) tail was not affected, its 3' untranslated region (3'UTR) lengthened. Overexpressing PAPOLA caused CCND1 mRNA 3'UTR shortening with a concomitant increase in the amount of corresponding transcript and protein, resulting in growth arrest in MCF-7 cells and DNA damage in HEK-293 cells. Such overexpression of PAPOLA promoted proliferation in the p53 mutant MDA-MB-231 cells. Our data suggest that PAPOLA is a possible candidate target for the control of tumor growth that is mostly relevant to triple-negative tumors, a group characterized by PAPOLA overexpression and lack of alternative targeted therapies.

Streamlining differential exon and 3' UTR usage with diffUTR

Background

Despite the importance of alternative poly-adenylation and 3′ UTR length for a variety of biological phenomena, there are limited means of detecting UTR changes from standard transcriptomic data.

Results

We present the diffUTR Bioconductor package which streamlines and improves upon differential exon usage (DEU) analyses, and leverages existing DEU tools and alternative poly-adenylation site databases to enable differential 3′ UTR usage analysis. We demonstrate the diffUTR features and show that it is more flexible and more accurate than state-of-the-art alternatives, both in simulations and in real data.

Conclusions

diffUTR enables differential 3′ UTR analysis and more generally facilitates DEU and the exploration of their results.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04114-7.

On the function and relevance of alternative 3'-UTRs in gene expression regulation

Messanger RNA (mRNA) isoforms with alternative 3'-untranslated regions (3'-UTRs) are produced by alternative polyadenylation (APA), which occurs during transcription in most eukaryotic genes. APA fine-tunes gene expression in a cell-type- and cellular state-dependent manner. Selection of an APA site entails the binding of core cleavage and polyadenylation factors to a particular polyadenylation site localized in the pre-mRNA and is controlled by multiple regulatory determinants, including transcription, pre-mRNA cis-regulatory sequences, and protein factors. Alternative 3'-UTRs serve as platforms for specific RNA binding proteins and microRNAs, which regulate gene expression in a coordinated manner by controlling mRNA fate and function in the cell. Genome-wide studies illustrated the full extent of APA prevalence and revealed that specific 3'-UTR profiles are associated with particular cellular states and diseases. Generally, short 3'-UTRs are associated with proliferative and cancer cells, and long 3'-UTRs are mostly found in polarized and differentiated cells. Fundamental new insights on the physiological consequences of this widespread event and the molecular mechanisms involved have been revealed through single-cell studies. Publicly available comprehensive databases that cover all APA mRNA isoforms identified in many cellular states and diseases reveal specific APA signatures. Therapies tackling APA mRNA isoforms or APA regulators may be regarded as innovative and attractive tools for diagnostics or treatment of several pathologies. We highlight the function of APA and alternative 3'-UTRs in gene expression regulation, the control of these mechanisms, their physiological consequences, and their potential use as new biomarkers and therapeutic tools. This article is categorized under: RNA Processing > 3' End Processing RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA in Disease and Development > RNA in Disease.

Cleavage and Polyadenylation Specific Factor 1 Promotes Tumor Progression via Alternative Polyadenylation and Splicing in Hepatocellular Carcinoma

Alternative polyadenylation (APA) is an important post-transcriptional regulatory mechanism required for cleavage and polyadenylation (CPA) of the 3′ untranslated region (3′ UTR) of mRNAs. Several aberrant APA events have been reported in hepatocellular carcinoma (HCC). However, the regulatory mechanisms underlying APA remain unclear. In this study, we found that the expression of cleavage and polyadenylation specific factor 1 (CPSF1), a major component of the CPA complex, was significantly increased in HCC tissues and correlated with unfavorable survival outcomes. Knockdown of CPSF1 inhibited HCC cell proliferation and migration, whereas overexpression of CPSF1 caused the opposite effect. Based on integrative analysis of Iso-Seq and RNA-seq data from HepG2.2.15 cells, we identified a series of transcripts with differential 3′ UTR lengths following the knockdown of CPSF1. These transcripts were related to the biological functions of gene transcription, cytoskeleton maintenance, and endomembrane system transportation. Moreover, knockdown of CPSF1 induced an increase in alternative splicing (AS) events in addition to APA. Taken together, this study provides new insights into our understanding of the post-transcriptional regulatory mechanisms in HCC and implies that CPSF1 may be a potential prognostic biomarker and therapeutic target for HCC.

CPSF6 links alternative polyadenylation to metabolism adaption in hepatocellular carcinoma progression

BACKGROUND

Alternative polyadenylation (APA) is an important mechanism of gene expression regulation through generation of RNA isoforms with distinct 3' termini. Increasing evidence has revealed that APA is actively involved in development and disease, including hepatocellular carcinoma (HCC). However, how APA functions in tumor formation and progression remains elusive. In this study, we investigated the role of cleavage factor I (CFIm) subunit CPSF6 in human hepatocellular carcinoma (HCC).

METHODS

Expression levels of CPSF6 in clinical tissues and cell lines were determined by qRT-PCR and western blot. Functional assays, including the cell number, MTT, colony formation and transwell, were used to determine the oncogenic role of CPSF6 in HCC. Animal experiments were used to determine the role of CPSF6 in HCC tumorigenicity in vivo. Deep sequencing-based 3 T-seq was used to profile the transcriptome-wide APA sites in both HCC cells and CPSF6 knockdown HCC cells. The function of CPSF6-affected target NQO1 with distinct 3'UTRs was characterized by metabolism assays.

RESULTS

We observed CPSF6 was upregulated in HCC and the high expression of CPSF6 was associated with poor prognosis in patients. Overexpression of CPSF6 promoted proliferation, migration and invasion of HCC cells in vitro and in vivo. Transcriptome-wide APA profiling analysis indicated that high expression of CPSF6 promoted the favorable usage of the proximal poly(A) site in the 3'UTR of NQO1. We demonstrated CPSF6-induced tumorigenic activities were mediated by the NQO1 isoform with short 3'UTR. Furthermore, we found that CPSF6 induced metabolic alterations in liver cells through NQO1.

CONCLUSION

CPSF6 plays a critical role in HCC progression by upregulating NQO1 expression through APA. These findings provide evidence to demonstrate that APA of NQO1 contributes to HCC progression and may have implications for developing new therapeutic strategy against this disease.

Profiling of alternative polyadenylation and gene expression in PEDV-infected IPEC-J2 cells

Since 2010, porcine epidemic diarrhea virus (PEDV) has received global attention with the emergence of variant strains characterized with high pathogenicity. The pathogen-host interaction after PEDV infection is still unclear. To investigate this issue, high-throughput-based sequencing technology is one of the optimal choices. In this study, we used in vitro transcription sequencing alternative polyadenylation sites (IVT-SAPAS) method, which allowed accurate profiling of gene expression and alternative polyadenylation (APA) sites to profile APA switching genes and differentially expressed genes (DEGs) in IPEC-J2 cells during PEDV variant strain infection. We found 804 APA switching genes, including switching in tandem 3' UTRs and switching between coding region and 3' UTR, and 1,677 DEGs in host after PEDV challenge. These genes participated in variety of biological processes such as cellular process, metabolism and immunity reactions. Moreover, 413 genes, most of which are the "focus" genes in interaction networks, were found to be involved in both APA switching genes and DEGs, suggesting these genes were synchronously regulated by different mechanisms. In summary, our results gave a relatively comprehensive insight into dynamic host-pathogen interactions in the regulation of host gene transcripts during PEDV infection.

Advances in the Bioinformatics Knowledge of mRNA Polyadenylation in Baculovirus Genes

Baculoviruses are a group of insect viruses with large circular dsDNA genomes exploited in numerous biotechnological applications, such as the biological control of agricultural pests, the expression of recombinant proteins or the gene delivery of therapeutic sequences in mammals, among others. Their genomes encode between 80 and 200 proteins, of which 38 are shared by all reported species. Thanks to multi-omic studies, there is remarkable information about the baculoviral proteome and the temporality in the virus gene expression. This allows some functional elements of the genome to be very well described, such as promoters and open reading frames. However, less information is available about the transcription termination signals and, consequently, there are still imprecisions about what are the limits of the transcriptional units present in the baculovirus genomes and how is the processing of the 3′ end of viral mRNA. Regarding to this, in this review we provide an update about the characteristics of DNA signals involved in this process and we contribute to their correct prediction through an exhaustive analysis that involves bibliography information, data mining, RNA structure and a comprehensive study of the core gene 3′ ends from 180 baculovirus genomes.

Emerging roles of the MAGE protein family in stress response pathways

The melanoma antigen (MAGE) proteins all contain a MAGE homology domain. MAGE genes are conserved in all eukaryotes and have expanded from a single gene in lower eukaryotes to ∼40 genes in humans and mice. Whereas some MAGEs are ubiquitously expressed in tissues, others are expressed in only germ cells with aberrant reactivation in multiple cancers. Much of the initial research on MAGEs focused on exploiting their antigenicity and restricted expression pattern to target them with cancer immunotherapy. Beyond their potential clinical application and role in tumorigenesis, recent studies have shown that MAGE proteins regulate diverse cellular and developmental pathways, implicating them in many diseases besides cancer, including lung, renal, and neurodevelopmental disorders. At the molecular level, many MAGEs bind to E3 RING ubiquitin ligases and, thus, regulate their substrate specificity, ligase activity, and subcellular localization. On a broader scale, the MAGE genes likely expanded in eutherian mammals to protect the germline from environmental stress and aid in stress adaptation, and this stress tolerance may explain why many cancers aberrantly express MAGEs Here, we present an updated, comprehensive review on the MAGE family that highlights general characteristics, emphasizes recent comparative studies in mice, and describes the diverse functions exerted by individual MAGEs.

scAPAtrap: identification and quantification of alternative polyadenylation sites from single-cell RNA-seq data

Alternative polyadenylation (APA) generates diverse mRNA isoforms, which contributes to transcriptome diversity and gene expression regulation by affecting mRNA stability, translation and localization in cells. The rapid development of 3' tag-based single-cell RNA-sequencing (scRNA-seq) technologies, such as CEL-seq and 10x Genomics, has led to the emergence of computational methods for identifying APA sites and profiling APA dynamics at single-cell resolution. However, existing methods fail to detect the precise location of poly(A) sites or sites with low read coverage. Moreover, they rely on priori genome annotation and can only detect poly(A) sites located within or near annotated genes. Here we proposed a tool called scAPAtrap for detecting poly(A) sites at the whole genome level in individual cells from 3' tag-based scRNA-seq data. scAPAtrap incorporates peak identification and poly(A) read anchoring, enabling the identification of the precise location of poly(A) sites, even for sites with low read coverage. Moreover, scAPAtrap can identify poly(A) sites without using priori genome annotation, which helps locate novel poly(A) sites in previously overlooked regions and improve genome annotation. We compared scAPAtrap with two latest methods, scAPA and Sierra, using scRNA-seq data from different experimental technologies and species. Results show that scAPAtrap identified poly(A) sites with higher accuracy and sensitivity than competing methods and could be used to explore APA dynamics among cell types or the heterogeneous APA isoform expression in individual cells. scAPAtrap is available at https://github.com/BMILAB/scAPAtrap.

Poly(A)-DG: A deep-learning-based domain generalization method to identify cross-species Poly(A) signal without prior knowledge from target species

In eukaryotes, polyadenylation (poly(A)) is an essential process during mRNA maturation. Identifying the cis-determinants of poly(A) signal (PAS) on the DNA sequence is the key to understand the mechanism of translation regulation and mRNA metabolism. Although machine learning methods were widely used in computationally identifying PAS, the need for tremendous amounts of annotation data hinder applications of existing methods in species without experimental data on PAS. Therefore, cross-species PAS identification, which enables the possibility to predict PAS from untrained species, naturally becomes a promising direction. In our works, we propose a novel deep learning method named Poly(A)-DG for cross-species PAS identification. Poly(A)-DG consists of a Convolution Neural Network-Multilayer Perceptron (CNN-MLP) network and a domain generalization technique. It learns PAS patterns from the training species and identifies PAS in target species without re-training. To test our method, we use four species and build cross-species training sets with two of them and evaluate the performance of the remaining ones. Moreover, we test our method against insufficient data and imbalanced data issues and demonstrate that Poly(A)-DG not only outperforms state-of-the-art methods but also maintains relatively high accuracy when it comes to a smaller or imbalanced training set.

Profiling Alternative 3' Untranslated Regions in Sorghum using RNA-seq Data

Sorghum is an important crop widely used for food, feed, and fuel. Transcriptome-wide studies of 3′ untranslated regions (3′UTR) using regular RNA-seq remain scarce in sorghum, while transcriptomes have been characterized extensively using Illumina short-read sequencing platforms for many sorghum varieties under various conditions or developmental contexts. 3′UTR is a critical regulatory component of genes, controlling the translation, transport, and stability of messenger RNAs. In the present study, we profiled the alternative 3′UTRs at the transcriptome level in three genetically related but phenotypically contrasting lines of sorghum: Rio, BTx406, and R9188. A total of 1,197 transcripts with alternative 3′UTRs were detected using RNA-seq data. Their categorization identified 612 high-confidence alternative 3′UTRs. Importantly, the high-confidence alternative 3′UTR genes significantly overlapped with the genesets that are associated with RNA N⁶-methyladenosine (m⁶A) modification, suggesting a clear indication between alternative 3′UTR and m⁶A methylation in sorghum. Moreover, taking advantage of sorghum genetics, we provided evidence of genotype specificity of alternative 3′UTR usage. In summary, our work exemplifies a transcriptome-wide profiling of alternative 3′UTRs using regular RNA-seq data in non-model crops and gains insights into alternative 3′UTRs and their genotype specificity.

A comprehensive rat transcriptome built from large scale RNA-seq-based annotation

The rat is an important model organism in biomedical research for studying human disease mechanisms and treatments, but its annotated transcriptome is far from complete. We constructed a Rat Transcriptome Re-annotation named RTR using RNA-seq data from 320 samples in 11 different organs generated by the SEQC consortium. Totally, there are 52 807 genes and 114 152 transcripts in RTR. Transcribed regions and exons in RTR account for ∼42% and ∼6.5% of the genome, respectively. Of all 73 074 newly annotated transcripts in RTR, 34 213 were annotated as high confident coding transcripts and 24 728 as high confident long noncoding transcripts. Different tissues rather than different stages have a significant influence on the expression patterns of transcripts. We also found that 11 715 genes and 15 852 transcripts were expressed in all 11 tissues and that 849 house-keeping genes expressed different isoforms among tissues. This comprehensive transcriptome is freely available at http://www.unimd.org/rtr/. Our new rat transcriptome provides essential reference for genetics and gene expression studies in rat disease and toxicity models.