Genome-wide profiling reveals alternative polyadenylation of mRNA in human non-small cell lung cancer

ZFP36L1 and ZFP36L2 reduce cyclin D1 expression by decreasing expression of E2F1 and long 3'UTR isoform of CCND1 transcripts

The CCND1 mRNA possesses at least two distinct lengths of the 3'-untranslated region (3'UTR), with the long isoform containing multiple AU-rich elements (AREs). The tandem zinc finger (TZF) domains of human ZFP36 family members have the capacity to bind to AREs and promote mRNA degradation. Our previous study demonstrated that mutations in the TZF domain of ZFP36L1 or ZFP36L2 increased the CCND1 expression. In this study, we investigated whether ZFP36L1 and ZFP36L2 could downregulate the expression of the long 3'UTR isoform of CCND1 mRNA in human colorectal cancer (CRC) cells. Firstly, the Gene Expression Profiling Interactive Analysis 2 database indicated downregulation of ZFP36 and ZFP36L1, while E2F1 and CCND1 were upregulated in human CRC tissues compared to normal colorectal tissues. Overexpression of ZFP36L1 and/or ZFP36L2 in T-REx-293, DLD-1, and HCT116 cells led to a decrease in the total CCND1, long isoform ratio of CCND1 mRNA, and E2F1 expression. Conversely, knockdown of ZFP36L1 and ZFP36L2 in HCT116 cells resulted in an increase in total CCND1, long isoform ratio of CCND1 mRNA, and E2F1 expression. Knockdown of E2F1 decreased CCND1 expression, indicating a potential role for E2F1 in regulating CCND1 expression at the transcriptional level. These findings suggest that ZFP36L1 and ZFP36L2 play a negative role in CCND1 expression. The underlying mechanisms might involve the reduction of E2F1 transactivation at the transcriptional level and the promotion of AREs-mediated decay of the long 3'UTR isoform of CCND1 through posttranscriptional processes.

Expression of CSTF2 in oral squamous cell carcinoma and its relationship with immune infiltration and poor prognosis

Background

Oral squamous cell carcinoma (OSCC) is a prevalent and devastating malignancy of the oral cavity that profoundly affects patient survival and quality of life (QOL). Cleavage Stimulation Factor Subunit 2 (CSTF2) is known to influence tumor development across multiple cancer types. However, its specific association with patient prognosis and immune cell infiltration in OSCC remains insufficiently understood.

Methods

To assess the expression levels and prognostic implications of CSTF2 in OSCC, comprehensive data were acquired from The Cancer Genome Atlas (TCGA) and subsequently normalized. Immunohistochemical staining of tissue microarrays was performed to analyze CSTF2 expression in the OSCC samples. Differences in CSTF2 expression between OSCC and adjacent non-cancerous samples were evaluated using the Wilcoxon rank-sum test. Functional enrichment analyses have been performed to identify biological pathways and functions associated with CSTF2. The relationship between the infiltration of various immune cells and CSTF2 expression levels was assessed using single-sample gene set enrichment analysis (ssGSEA). Ultimately, the prognostic significance of CSTF2 was evaluated through Kaplan-Meier survival analysis, in conjunction with univariate and multivariate Cox regression analyses, as well as receiver operating characteristic (ROC) curves.

Results

High CSTF2 expression was observed in OSCC and associated with unfavorable clinicopathological variables, including histological grade and lymphnode neck dissection. Functional enrichment analysis indicated that CSTF2 plays a role in epidermal development and differentiation, immunoglobulin complexes, peptidases and endopeptidase inhibitor activity, and cytochrome P450 metabolic processes. Additionally, the overexpression of CSTF2 exhibited a negative correlation with the infiltration of immature dendritic cells (iDCs), cytotoxic cells, and plasmacytoid dendritic cells (pDCs). Notably, elevated CSTF2 expression is significantly associated with reduced patient outcomes.

Conclusion

Elevated CSTF2 expression in OSCC is associated with poor prognostic outcomes, highlighting its capacity to function as an innovative prognostic biomarker and a target for therapeutic interventions.

Alternative polyadenylation in cancer: Molecular mechanisms and clinical application

Alternative polyadenylation (APA) serves as a crucial mechanism for the posttranscriptional regulation of gene expression and influences gene expression by generating diverse mRNA isoforms. This process is regulated by a diverse array of RNA-binding proteins (RBPs), which selectively bind to specific sequences or structures within the pre-mRNA molecule. Dysregulation of APA and its associated RBPs has been implicated in numerous diseases, including cardiovascular diseases, nervous system disease, and cancer. For instance, aberrant APA events have been observed in several types of tumors, contributing to tumor heterogeneity and affecting key cellular pathways involved in cell proliferation, invasion, metastasis, and response to therapy. This review critically evaluates the current understanding of APA mechanisms and the multifaceted roles of RBPs in orchestrating this intricate process. We highlight recent advancements in high-throughput sequencing and bioinformatics tools that have enhanced our ability to study APA on a genome-wide scale. Moreover, we explored the pathological consequences of APA dysregulation, emphasizing its role in oncogenesis. By elucidating the intricate relationships between APA and RBPs, this review aims to underscore the potential of targeting the APA machinery and RBPs for therapeutic intervention. Understanding these molecular processes holds promise for developing novel diagnostic markers and treatment strategies for a range of human cancers.

Exploring vimentin's role in breast cancer via PICK1 alternative polyadenylation and the miR-615-3p-PICK1 interaction

Breast cancer continues to be a major health issue for women worldwide, with vimentin (VIM) identified as a crucial factor in its progression due to its role in cell migration and the epithelial-to-mesenchymal transition (EMT). This study focuses on elucidating VIM's regulatory mechanisms on the miR-615-3p/PICK1 axis. Utilizing the 4T1 breast cancer cell model, we first used RNA-seq and proteomics to investigate the changes in the APA of PICK1 following VIM knockout (KO). These high-throughput analyses aimed to uncover the underlying transcriptional and proteomic alterations associated with VIM's influence on breast cancer cells. RNA-seq and proteomic profiling revealed significant APA in PICK1 following VIM KO, suggesting a novel mechanism by which VIM regulates breast cancer progression. Validation experiments confirmed that VIM KO affects the miR-615-3p-PICK1 axis, with miR-615-3p's regulation of PICK1 being contingent upon the APA of PICK1. These findings highlight the complex interplay between VIM, miR-615-3p, and PICK1 in the regulation of breast cancer cell behavior. This study reveals that vimentin affects the miR-615-3p-PICK1 axis through APA, revealing the key role of VIM in cancer progression. Opened up new avenues for targeted cancer therapy, with a focus on regulating the interaction between APA and miR-615-3p-PICK1.

InPAS: An R/Bioconductor Package for Identifying Novel Polyadenylation Sites and Alternative Polyadenylation from Bulk RNA-seq Data

BACKGROUND

Alternative cleavage and polyadenylation (APA) is a crucial post-transcriptional gene regulation mechanism that regulates gene expression in eukaryotes by increasing the diversity and complexity of both the transcriptome and proteome. Despite the development of more than a dozen experimental methods over the last decade to identify and quantify APA events, widespread adoption of these methods has been limited by technical, financial, and time constraints. Consequently, APA remains poorly understood in most eukaryotes. However, RNA sequencing (RNA-seq) technology has revolutionized transcriptome profiling and recent studies have shown that RNA-seq data can be leveraged to identify and quantify APA events.

RESULTS

To fully capitalize on the exponentially growing RNA-seq data, we developed InPAS (Identification of Novel alternative PolyAdenylation Sites), an R/Bioconductor package for accurate identification of novel and known cleavage and polyadenylation sites (CPSs), as well as quantification of APA from RNA-seq data of various experimental designs. Compared to other APA analysis tools, InPAS offers several important advantages, including the ability to detect both novel proximal and distal CPSs, to fine tune positions of CPSs using a naïve Bayes classifier based on flanking sequence features, and to identify APA events from RNA-seq data of complex experimental designs using linear models. We benchmarked the performance of InPAS and other leading tools using simulated and experimental RNA-seq data with matched 3'-end RNA-seq data. Our results reveal that InPAS frequently outperforms existing tools in terms of precision, sensitivity, and specificity. Furthermore, we demonstrate its scalability and versatility by applying it to large, diverse RNA-seq datasets.

CONCLUSIONS

InPAS is an efficient and robust tool for identifying and quantifying APA events using readily accessible conventional RNA-seq data. Its versatility opens doors to explore APA regulation across diverse eukaryotic systems with various experimental designs. We believe that InPAS will drive APA research forward, deepening our understanding of its role in regulating gene expression, and potentially leading to the discovery of biomarkers or therapeutics for diseases.

Alternative polyadenylation determines the functional landscape of inverted Alu repeats

Inverted Alu repeats (IRAlus) are abundantly found in the transcriptome, especially in introns and 3' untranslated regions (UTRs). Yet, the biological significance of IRAlus embedded in 3' UTRs remains largely unknown. Here, we find that 3' UTR IRAlus silences genes involved in essential signaling pathways. We utilize J2 antibody to directly capture and map the double-stranded RNA structure of 3' UTR IRAlus in the transcriptome. Bioinformatic analysis reveals alternative polyadenylation as a major axis of IRAlus-mediated gene regulation. Notably, the expression of mouse double minute 2 (MDM2), an inhibitor of p53, is upregulated by the exclusion of IRAlus during UTR shortening, which is exploited to silence p53 during tumorigenesis. Moreover, the transcriptome-wide UTR lengthening in neural progenitor cells results in the global downregulation of genes associated with neurodegenerative diseases, including amyotrophic lateral sclerosis, via IRAlus inclusion. Our study establishes the functional landscape of 3' UTR IRAlus and its role in human pathophysiology.

Transitional Insight into the RNA-Based Oligonucleotides in Cancer Treatment

Conventional cancer therapies with chemodrugs suffer from various disadvantages, such as irreversible side effects on the skin, heart, liver, and nerves with even fatal consequences. RNA-based therapeutic is a novel technology which offers great potential as non-toxic, non-infectious, and well-tolerable platform. Herein, we introduce different RNA-based platforms with a special focus on siRNA, miRNA, and mRNA applications in cancer treatment in order to better understand the details of their therapeutic effects. Of note, the co-delivery of RNAs with other distinct RNA or drugs has provided safe, efficient, and novel treatment modalities for cancer treatment.

Transcriptome Responses to Different Salinity Conditions in Litoditis marina, Revealed by Long-Read Sequencing

The marine nematode Litoditis marina is widely distributed in intertidal zones around the globe, yet the mechanisms underlying its broad adaptation to salinity remain elusive. In this study, we applied ONT long-read sequencing technology to unravel the transcriptome responses to different salinity conditions in L. marina. Through ONT sequencing under 3‱, 30‱ and 60‱ salinity environments, we obtained 131.78 G clean data and 26,647 non-redundant long-read transcripts, including 6464 novel transcripts. The DEGs obtained from the current ONT lrRNA-seq were highly correlated with those identified in our previously reported Illumina short-read RNA sequencing data. When we compared the 30‱ to the 3‱ salinity condition, we found that GO terms such as oxidoreductase activity, cation transmembrane transport and ion transmembrane transport were shared between the ONT lrRNA-seq and Illumina data. Similarly, GO terms including extracellular space, structural constituents of cuticle, substrate-specific channel activity, ion transport and substrate-specific transmembrane transporter activity were shared between the ONT and Illumina data under 60‱ compared to 30‱ salinity. In addition, we found that 79 genes significantly increased, while 119 genes significantly decreased, as the salinity increased. Furthermore, through the GO enrichment analysis of 214 genes containing DAS, in 30‱ compared to 3‱ salinity, we found that GO terms such as cellular component assembly and coenzyme biosynthetic process were enriched. Additionally, we observed that GO terms such as cellular component assembly and coenzyme biosynthetic process were also enriched in 60‱ compared to 30‱ salinity. Moreover, we found that 86, 125, and 81 genes that contained DAS were also DEGs, in comparisons between 30‱ and 3‱, 60‱ and 30‱, and 60‱ and 3‱ salinity, respectively. In addition, we demonstrated the landscape of alternative polyadenylation in marine nematode under different salinity conditions This report provides several novel insights for the further study of the mechanisms by which euryhalinity formed and evolved, and it might also contribute to the investigation of salinity dynamics induced by global climate change.

The role of CSTF2 in cancer: from technology to clinical application

A protein called cleavage-stimulating factor subunit 2 (CSTF2, additionally called CSTF-64) binds RNA and is needed for the cleavage and polyadenylation of mRNA. CSTF2 is an important component subunit of the cleavage stimulating factor (CSTF), which is located on the X chromosome and encodes 557 amino acids. There is compelling evidence linking elevated CSTF2 expression to the pathological advancement of cancer and on its impact on the clinical aspects of the disease. The progression of cancers, including hepatocellular carcinoma, melanoma, prostate cancer, breast cancer, and pancreatic cancer, is correlated with the upregulation of CSTF2 expression. This review provides a fresh perspective on the investigation of the associations between CSTF2 and various malignancies and highlights current studies on the regulation of CSTF2. In particular, the mechanism of action and potential clinical applications of CSTF2 in cancer suggest that CSTF2 can serve as a new biomarker and individualized treatment target for a variety of cancer types.

A novel APA-based prognostic signature may predict the prognosis of lung adenocarcinoma in an East Asian population

The role of alternative polyadenylation (APA) in tumor development is becoming increasingly evident, but the impact of APA events on the prognosis of LUAD patients is unclear. Therefore, in the present study, we aimed to analyze specific APA events in LUAD to identify novel prognostic biomarkers for LUAD. We first identified prognostic candidate genes for LUAD associated with APA events and validated them in both the East Asian and the USA cohorts, finding that five genes (DCUN1D5, PSMC4, TFAM, THRA, and TMEM100) were of prognostic significance in both populations. Based on this, an APA-based prognostic signature was constructed for the East Asian population. The predictive accuracy of the prognostic signature was further evaluated by the time-dependent ROC, with 1-, 2-, and 3-year AUCs of 0.86, 0.81, and 0.71, respectively. This study may provide new markers for individualized diagnosis and prognostic assessment of LUAD and potential targets for precision treatment.

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.

The Transition from Cancer "omics" to "epi-omics" through Next- and Third-Generation Sequencing

Deciphering cancer etiopathogenesis has proven to be an especially challenging task since the mechanisms that drive tumor development and progression are far from simple. An astonishing amount of research has revealed a wide spectrum of defects, including genomic abnormalities, epigenomic alterations, disturbance of gene transcription, as well as post-translational protein modifications, which cooperatively promote carcinogenesis. These findings suggest that the adoption of a multidimensional approach can provide a much more precise and comprehensive picture of the tumor landscape, hence serving as a powerful tool in cancer research and precision oncology. The introduction of next- and third-generation sequencing technologies paved the way for the decoding of genetic information and the elucidation of cancer-related cellular compounds and mechanisms. In the present review, we discuss the current and emerging applications of both generations of sequencing technologies, also referred to as massive parallel sequencing (MPS), in the fields of cancer genomics, transcriptomics and proteomics, as well as in the progressing realms of epi-omics. Finally, we provide a brief insight into the expanding scope of sequencing applications in personalized cancer medicine and pharmacogenomics.

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.

CSTF2 Acts as a Prognostic Marker Correlated with Immune Infiltration in Hepatocellular Carcinoma

Background

Cleavage stimulation factor 2 (CSTF2) encodes a nuclear protein that is implicated in the development of various cancers. However, the role of CSTF2 in hepatocellular carcinoma (HCC) has not been understood. This study aims to explore the function of CSTF2 in HCC.

Methods

The expression, diagnostic capability, prognostic value, and immune cell effect of CSTF2 in HCC were explored using various databases. The expression level of CSTF2 were validated in our cell lines. The effect of CSTF2 on hepatocarcinogenesis was explored by CSTF2 silencing.

Results

CSTF2 expression was significantly elevated in HCC and correlated with multiple clinicopathological characteristics. CSTF2 exhibited good diagnostic capability in discriminating HCC samples from nontumorous samples. High CSTF2 expression was significantly related to poor overall survival. Univariate and multivariate Cox regression analyses suggested that CSTF2 expression was an independent risk factor for HCC. These results were validated in ICGC cohorts. In addition, the nomogram based on CSTF2 showed better predictive performance than the AJCC staging system in TCGA and ICGC cohorts. Functional enrichment analysis revealed that CSTF2-related genes were involved in DNA/RNA processing and the cell cycle. In addition, we found that CSTF2 expression was closely related to the levels of various infiltrating immune cells, especially neutrophils. Moreover, some immune checkpoints had positive relationships with CSTF2 expression. CSTF2 silencing inhibited proliferation, invasion and migration, and promoted apoptosis in HepG2 cells. Western blotting analysis revealed that CSTF2 silencing inactivated the Wnt/β-catenin signaling pathway.

Conclusion

High CSTF2 expression not only correlates with unfavorable outcomes but also affects immune cell infiltration and immune checkpoint expression in HCC. CSTF2 silencing can alleviate the malignant phenotypes of hepatic cancer cell by inactivating the Wnt/β-catenin signaling pathway. These results indicate that CSTF2 can serve as a promising prognostic marker and therapeutic target for HCC patients.

Alternative polyadenylation writer CSTF2 forms a positive loop with FGF2 to promote tubular epithelial-mesenchymal transition and renal fibrosis

Effective therapies for renal fibrosis, the common endpoint for most kidney diseases, are lacking. We previously reported that alternative polyadenylation (APA) drives transition from acute kidney injury to chronic kidney disease, suggesting a potential role for APA in renal fibrogenesis. Here, we found that among canonical APA writers, CSTF2 expression was upregulated in tubular epithelial cells (TEC) of fibrotic kidneys. CSTF2 was also identified as a TGF-β-inducible pro-fibrotic gene. Further analysis revealed that CSTF2 promoted epithelial-mesenchymal transition (EMT) and extracellular matrix (ECM) overproduction in TEC by inducing 3'UTR shortening and upregulation of the expression of basic fibroblast growth factor 2 (FGF2). Additionally, 3'UTR shortening stabilised FGF2 mRNA through miRNA evasion. Interestingly, FGF2 enhanced CSTF2 expression, leading to the forming of a CSTF2-FGF2 positive loop in TEC. Furthermore, CSTF2 knockdown alleviated unilateral ureteral obstruction-induced renal fibrosis in vivo. Finally, we developed a CSTF2-targeted antisense oligonucleotide (ASO) and validated its effectiveness in vitro. These results indicate that the expression of the APA writer, CSTF2, is upregulated by TGF-β and CSTF2 facilitates TGF-β-induced FGF2 overexpression, forming a TGF-β-CSTF2-FGF2 pro-fibrotic axis in TEC. CSTF2 is a potentially promising target for renal fibrosis that does not directly disrupt TGF-β.

Identification and Verification of an Alternative Polyadenylation-Related lncRNA Prognostic Signature for Glioma

Due to the high mortality and modality of glioma, it was urgently needed to develop a glioma prognostic assessment system. Previous studies demonstrated that alternative polyadenylation- (APA-) related genes are important in immune response and oncogenesis. mRNA and lncRNA expression information of glioma samples were acquired from CGGA and TCGA databases, and lncRNAs associated with APA were selected through correlation analysis. The prognosis model of APA-related lncRNAs was built by the univariate Cox, random forest, and univariate Cox regression analyses. Glioma samples were assigned into high- and low-risk groups. Independence and effectiveness of the prognostic model were evaluated by Kaplan-Meier analysis, ROC curve, and Cox regression analyses. GO, KEGG enrichment, and GSEA analyses showed that the mainly involved signaling pathways were enriched in cellular immunity and immune signal transduction. We further analyzed expression differences of negative immune regulatory genes and immune cell infiltration degree between two groups. Immune checkpoints CTLA4 and LAG3 and immune suppressors TGFB, IL10, NOS3, and IDO1 and immune cell infiltration were notably upregulated in the high-risk group. The PD1/PDL1 expression was significantly correlated with risk score, showing that the prognostic model of APA-related lncRNA could effectively assess the tumor immune suppression. In conclusion, we established a risk assessment model of APA-related lncRNA in glioma, which could effectively evaluate prognosis of patients with glioma and tumor immune suppression and could provide guidance for clinical treatment.

Investigation of the effect of meclofenamic acid on the proteome of LNCaP cells reveals changes in alternative polyadenylation and splicing machinery

Prostate cancer is the most common type of cancer among men, and there is still no definitively effective drug treatment. Thus, the search for novel drug agents that may be used for the effective treatment continues. Meclofenamic acid (MA), a non-steroidal anti-inflammatory drug, with anti-tumor effects in various types of cancers was used to investigate its effects on LNCaP cells, a prostate cancer cell line, at the proteome level. The cells were treated with 80 µM MA for 24 h and a comparative proteomic analysis was performed with their untreated control cells. Proteins were extracted from the cells and then were subjected to two-dimensional gel electrophoresis. Protein spots displaying changes in their regulation ratios for more than two-fold were excised from the gels and identified with MALDI-TOF/TOF mass spectrometry. Bioinformatics analysis of the differentially regulated proteins that we identified showed that they were all associated with and took part in related pathways. Glycolytic pathway, cytoskeletal formation, transport activity, protein metabolism, and most notably an mRNA processing pathway were affected by the MA treatment. In addition to presenting a detailed information for what is happening inside the cells upon MA treatment, the proteins affected by MA treatment hold the potential to be novel targets for prostate cancer treatment provided that further in vivo experiments are carried out.

Cleavage stimulation factor 2 promotes malignant progression of liver hepatocellular carcinoma by activating phosphatidylinositol 3'-kinase/protein kinase B/mammalian target of rapamycin pathway

Liver hepatocellular carcinoma (LIHC) is the most common type, comprising 75-85% of all liver malignancies. We investigated the roles of cleavage stimulation factor 2 (CSTF2) in LIHC and explored the underlying mechanisms. CSTF2 expression and its association with LIHC patient survival probability were analyzed with The Cancer Genome Atlas. CSTF2 expression in LIHC cells was assessed using western blot and quantitative real-time PCR. Alterations in CSTF2 expression were induced by cell transfection. Cell colony formation, apoptosis, proliferation, invasion, and migration were assessed using colony formation, flow cytometry, 5-ethynyl-2'-deoxyuridine, and transwell assays. Pathway enrichment analysis was performed using gene set enrichment analysis (GSEA). The expression of apoptosis-, metastasis-, and pathway-associated factors was determined via western blot. The pathway rescue assay was further performed using 740Y-P or Wortmannin. CSTF2 upregulation was observed in LIHC tissues and cells. Patients with high CSTF2 expression had a lower probability of overall survival. CSTF2 overexpression enhanced colony formation, proliferation, invasion and migration, while repressing apoptosis in LIHC cells. GSEA revealed that CSTF2 was mainly enriched in the phosphatidylinositol 3'-kinase/protein kinase B/mammalian target of rapamycin (PI3K/AKT/mTOR) pathway. Western blot analysis proved that CSTF2 overexpression activated this pathway. CSTF2 knockdown yielded the opposite effects. 740Y-P, a PI3K activator, reversed the CSTF2 knockdown-triggered effects on cell proliferation, apoptosis, invasion, and migration. Moreover, Wortmannin, a PI3K inhibitor, also reversed the CSTF2 overexpression-induced effects on cell proliferation, apoptosis, invasion, and migration. These results indicated that CSTF2 overexpression might exacerbate the malignant phenotypes of LIHC cells via activation of PI3K/AKT/mTOR pathway.

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.

Immunological Classification of Pancreatic Carcinomas to Identify Immune Index and Provide a Strategy for Patient Stratification

Background

Cancer immunotherapy has produced significant positive clinical effects in a variety of tumor types. However, pancreatic ductal adenocarcinoma (PDAC) is widely considered to be a "cold" cancer with poor immunogenicity. Our aim is to determine the detailed immune features of PDAC to seek new treatment strategies.

Methods

The immune cell abundance of PDAC patients was evaluated with the single-sample gene set enrichment analysis (ssGSEA) using 119 immune gene signatures. Based on these data, patients were classified into different immune subtypes (ISs) according to immune gene signatures. We analyzed their response patterns to immunotherapy in the datasets, then established an immune index to reflect the different degrees of immune infiltration through linear discriminant analysis (LDA). Finally, potential prognostic markers associated with the immune index were identified based on weighted correlation network analysis (WGCNA) that was functionally validated in vitro.

Results

Three ISs were identified in PDAC, of which IS3 had the best prognosis across all three cohorts. The different expressions of immune profiles among the three ISs indicated a distinct responsiveness to immunotherapies in PDAC subtypes. By calculating the immune index, we found that the IS3 represented higher immune infiltration, while IS1 represented lower immune infiltration. Among the investigated signatures, we identified ZNF185, FANCG, and CSTF2 as risk factors associated with immune index that could potentially facilitate diagnosis and could be therapeutic target markers in PDAC patients.

Conclusions

Our findings identified immunologic subtypes of PDAC with distinct prognostic implications, which allowed us to establish an immune index to represent the immune infiltration in each subtype. These results show the importance of continuing investigation of immunotherapy and will allow clinical workers to personalized treatment more effectively in PDAC patients.

Alternative poly-adenylation modulates α1-antitrypsin expression in chronic obstructive pulmonary disease

α1-anti-trypsin (A1AT), encoded by SERPINA1, is a neutrophil elastase inhibitor that controls the inflammatory response in the lung. Severe A1AT deficiency increases risk for Chronic Obstructive Pulmonary Disease (COPD), however, the role of A1AT in COPD in non-deficient individuals is not well known. We identify a 2.1-fold increase (p = 2.5x10-6) in the use of a distal poly-adenylation site in primary lung tissue RNA-seq in 82 COPD cases when compared to 64 controls and replicate this in an independent study of 376 COPD and 267 controls. This alternative polyadenylation event involves two sites, a proximal and distal site, 61 and 1683 nucleotides downstream of the A1AT stop codon. To characterize this event, we measured the distal ratio in human primary tissue short read RNA-seq data and corroborated our results with long read RNA-seq data. Integrating these results with 3' end RNA-seq and nanoluciferase reporter assay experiments we show that use of the distal site yields mRNA transcripts with over 50-fold decreased translation efficiency and A1AT expression. We identified seven RNA binding proteins using enhanced CrossLinking and ImmunoPrecipitation precipitation (eCLIP) with one or more binding sites in the SERPINA1 3' UTR. We combined these data with measurements of the distal ratio in shRNA knockdown experiments, nuclear and cytoplasmic fractionation, and chemical RNA structure probing. We identify Quaking Homolog (QKI) as a modulator of SERPINA1 mRNA translation and confirm the role of QKI in SERPINA1 translation with luciferase reporter assays. Analysis of single-cell RNA-seq showed differences in the distribution of the SERPINA1 distal ratio among hepatocytes, macrophages, αβ-Tcells and plasma cells in the liver. Alveolar Type 1,2, dendritic cells and macrophages also vary in their distal ratio in the lung. Our work reveals a complex post-transcriptional mechanism that regulates alternative polyadenylation and A1AT expression in COPD.

Dynamic Variations of 3'UTR Length Reprogram the mRNA Regulatory Landscape

This paper concerns 3′-untranslated regions (3′UTRs) of mRNAs, which are non-coding regulatory platforms that control stability, fate and the correct spatiotemporal translation of mRNAs. Many mRNAs have polymorphic 3′UTR regions. Controlling 3′UTR length and sequence facilitates the regulation of the accessibility of functional effectors (RNA binding proteins, miRNAs or other ncRNAs) to 3′UTR functional boxes and motifs and the establishment of different regulatory landscapes for mRNA function. In this context, shortening of 3′UTRs would loosen miRNA or protein-based mechanisms of mRNA degradation, while 3′UTR lengthening would strengthen accessibility to these effectors. Alterations in the mechanisms regulating 3′UTR length would result in widespread deregulation of gene expression that could eventually lead to diseases likely linked to the loss (or acquisition) of specific miRNA binding sites. Here, we will review the mechanisms that control 3′UTR length dynamics and their alterations in human disorders. We will discuss, from a mechanistic point of view centered on the molecular machineries involved, the generation of 3′UTR variability by the use of alternative polyadenylation and cleavage sites, of mutually exclusive terminal alternative exons (exon skipping) as well as by the process of exonization of Alu cassettes to generate new 3′UTRs with differential functional features.

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.

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.

Emerging Roles of RNA 3'-end Cleavage and Polyadenylation in Pathogenesis, Diagnosis and Therapy of Human Disorders

A crucial feature of gene expression involves RNA processing to produce 3′ ends through a process termed 3′ end cleavage and polyadenylation (CPA). This ensures the nascent RNA molecule can exit the nucleus and be translated to ultimately give rise to a protein which can execute a function. Further, alternative polyadenylation (APA) can produce distinct transcript isoforms, profoundly expanding the complexity of the transcriptome. CPA is carried out by multi-component protein complexes interacting with multiple RNA motifs and is tightly coupled to transcription, other steps of RNA processing, and even epigenetic modifications. CPA and APA contribute to the maintenance of a multitude of diverse physiological processes. It is therefore not surprising that disruptions of CPA and APA can lead to devastating disorders. Here, we review potential CPA and APA mechanisms involving both loss and gain of function that can have tremendous impacts on health and disease. Ultimately we highlight the emerging diagnostic and therapeutic potential CPA and APA offer.