Proliferating cells express mRNAs with shortened 3' untranslated regions and fewer microRNA target sites

The transcriptional terminator XRN2 and the RNA-binding protein Sam68 link alternative polyadenylation to cell cycle progression in prostate cancer

Alternative polyadenylation (APA) yields transcripts differing in their 3'-end, and its regulation is altered in cancer, including prostate cancer. Here we have uncovered a mechanism of APA regulation impinging on the interaction between the exonuclease XRN2 and the RNA-binding protein Sam68, whose increased expression in prostate cancer is promoted by the transcription factor MYC. Genome-wide transcriptome profiling revealed a widespread impact of the Sam68/XRN2 complex on APA. XRN2 promotes recruitment of Sam68 to its target transcripts, where it competes with the cleavage and polyadenylation specificity factor for binding to strong polyadenylation signals at distal ends of genes, thus promoting usage of suboptimal proximal polyadenylation signals. This mechanism leads to 3' untranslated region shortening and translation of transcripts encoding proteins involved in G1/S progression and proliferation. Thus, our findings indicate that the APA program driven by Sam68/XRN2 promotes cell cycle progression and may represent an actionable target for therapeutic intervention.

The clinical utility of dysregulated microRNA expression in paediatric solid tumours

MicroRNAs (miRNAs) are short, non-protein-coding genes that regulate the expression of numerous protein-coding genes. Their expression is dysregulated in cancer, where they may function as oncogenes or tumour suppressor genes. As miRNAs are highly resistant to degradation, they are ideal biomarker candidates to improve the diagnosis and clinical management of cancer, including prognostication. Furthermore, miRNAs dysregulated in malignancy represent potential therapeutic targets. The use of miRNAs for these purposes is a particularly attractive option to explore for paediatric malignancies, where the mutational burden is typically low, in contrast to cancers affecting adult patients. As childhood cancers are rare, it has taken time to accumulate the necessary body of evidence showing the potential for miRNAs to improve clinical management across this group of tumours. Here, we review the current literature regarding the potential clinical utility of miRNAs in paediatric solid tumours, which is now both timely and justified. Exploring such avenues is warranted to improve the management and outcomes of children affected by cancer.

ReadZS detects cell type-specific and developmentally regulated RNA processing programs in single-cell RNA-seq

RNA processing, including splicing and alternative polyadenylation, is crucial to gene function and regulation, but methods to detect RNA processing from single-cell RNA sequencing data are limited by reliance on pre-existing annotations, peak calling heuristics, and collapsing measurements by cell type. We introduce ReadZS, an annotation-free statistical approach to identify regulated RNA processing in single cells. ReadZS discovers cell type-specific RNA processing in human lung and conserved, developmentally regulated RNA processing in mammalian spermatogenesis-including global 3' UTR shortening in human spermatogenesis. ReadZS also discovers global 3' UTR lengthening in Arabidopsis development, highlighting the usefulness of this method in under-annotated transcriptomes.

The relationship of mRNA with protein expression in CD8+ T cells associates with gene class and gene characteristics

T cells are key players in our defence against infections and malignancies. When T cells differentiate or become activated, they undergo substantial alterations in gene expression. Even though RNA expression levels are now well documented throughout different stages of T cells, it is not well understood how mRNA expression translates into the protein landscape. By combining paired RNA sequencing and mass spectrometry data of primary human CD8+ T cells, we report that mRNA expression is a poor proxy for the overall protein output, irrespective of the differentiation or activation status. Yet, gene class stratification revealed a function-specific correlation of mRNA with protein expression. This gene class-specific expression pattern associated with differences in gene characteristics such as sequence conservation and untranslated region (UTR) lengths. In addition, the presence of AU-rich elements in the 3'UTR associated with alterations in mRNA and protein abundance T cell activation dependent, gene class-specific manner. In conclusion, our study highlights the role of gene characteristics as a determinant for gene expression in T cells.

CPSF1 positively regulates NSDHL by alternative polyadenylation and promotes gastric cancer progression

Gastric cancer (GC) is a common malignancies with unfavourable prognosis. As one of the most common RNA modifications in nature, alternative polyadenylation (APA) plays a critical role in the progression of carcinomas. CPSF1 is a critical APA-related factor and is involved in many cancers. Nevertheless, the roles and underlying mechanisms of CPSF1 remain unclear in GC. In this work, we identified that CPSF1 is significantly upregulated in GC and that high CPSF1 expression indicates an unfavourable prognosis in GC patients. Moreover, CPSF1 expression levels were closely associated with tumour size, TNM stage and lymph node metastasis. CPSF1 depletion dramatically weakened GC cell proliferation and metastasis. We then performed RNA sequencing and found numerous downstream genes involved the regulation of CPSF1 with remarkable changes in 3'UTR length, among which NSDHL was positively regulated by CPSF1 and promoted GC progression. In addition, rescue assays demonstrated that NSDHL mediated the carcinogenic effect of CPSF1, and this process potentially involved APA. Therefore, this study showed that CPSF1 promotes GC progression, at least in part, by enhancing NSDHL and offered new insights into therapeutic targets for GC.

Characterization of Immune-Related Alternative Polyadenylation Events in Cancer Immunotherapy

Alternative polyadenylation (APA) is an important posttranscriptional modification commonly involved in tumor development. However, the functional roles of APA in tumor immunity remain largely unknown. Here, we performed an in-depth analysis of the 3'UTR usage of protein-coding genes and tumor immune response in 10,303 tumor samples across 31 cancer types to develop the immune-related APA event (ImmAPA) score pipeline, an integrated algorithm to characterize the regulatory landscape of APA events in cancer immunity-related pathways. Tumor-specific ImmAPAs that strongly correlate with immune cell infiltration and immune checkpoint blockade (ICB) treatment-related biomarkers were identified. Among these ImmAPAs, the top-ranking COL1A1 3'UTR usage was strongly associated with worse prognosis and tumor immune evasion. Furthermore, a machine learning approach to construct an ICB-related ImmAPA score model predicted immunotherapy efficacy. Overall, the characterization of immune-related APA that corresponds to tumor progression and tumor immunity highlights the clinical utility of APA events as potential biomarkers in cancer immunotherapy.

SIGNIFICANCE

Elucidation of the landscape of immune-related alternative polyadenylation in cancer identifies alternative polyadenylation events that may play a role in immune modulation and immunotherapy efficacy.

Introduction of the T315I gatekeeper mutation of BCR/ABL1 into a Philadelphia chromosome-positive lymphoid leukemia cell line using the CRISPR/Cas9 system

Imatinib and second-generation tyrosine kinase inhibitors (TKIs) have dramatically improved the prognosis of Philadelphia chromosome-positive (Ph+) acute lymphoblastic leukemia (ALL). However, overcoming TKI resistance due to the T315I gatekeeper mutation of BCR/ABL1 is crucial for further improving the prognosis. The clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 system is appropriate for establishing a human model of Ph+ ALL with the T315I mutation, because it can induce specific mutations via homologous recombination (HR) repair in cells with intact endogenous HR pathway. Here we used CRISPR/Cas9 to introduce the T315I mutation into the Ph+ lymphoid leukemia cell line KOPN55bi, which appeared to have an active HR pathway based on its resistance to a poly (ADP-Ribose) polymerase-1 inhibitor. Single-guide RNA targeting at codon 315 and single-strand oligodeoxynucleotide containing ACT to ATT nucleotide transition at codon 315 were electroporated with recombinant Cas9 protein. Dasatinib-resistant sublines were obtained after one-month selection with the therapeutic concentration of dasatinib, leading to T315I mutation acquisition through HR. T315I-acquired sublines were highly resistant to imatinib and second-generation TKIs but moderately sensitive to the therapeutic concentration of ponatinib. This authentic human model is helpful for developing new therapeutic strategies overcoming TKI resistance in Ph+ ALL due to T315I mutation.

APA-Scan: detection and visualization of 3'-UTR alternative polyadenylation with RNA-seq and 3'-end-seq data

Background

The eukaryotic genome is capable of producing multiple isoforms from a gene by alternative polyadenylation (APA) during pre-mRNA processing. APA in the 3′-untranslated region (3′-UTR) of mRNA produces transcripts with shorter or longer 3′-UTR. Often, 3′-UTR serves as a binding platform for microRNAs and RNA-binding proteins, which affect the fate of the mRNA transcript. Thus, 3′-UTR APA is known to modulate translation and provides a mean to regulate gene expression at the post-transcriptional level. Current bioinformatics pipelines have limited capability in profiling 3′-UTR APA events due to incomplete annotations and a low-resolution analyzing power: widely available bioinformatics pipelines do not reference actionable polyadenylation (cleavage) sites but simulate 3′-UTR APA only using RNA-seq read coverage, causing false positive identifications. To overcome these limitations, we developed APA-Scan, a robust program that identifies 3′-UTR APA events and visualizes the RNA-seq short-read coverage with gene annotations.

Methods

APA-Scan utilizes either predicted or experimentally validated actionable polyadenylation signals as a reference for polyadenylation sites and calculates the quantity of long and short 3′-UTR transcripts in the RNA-seq data. APA-Scan works in three major steps: (i) calculate the read coverage of the 3′-UTR regions of genes; (ii) identify the potential APA sites and evaluate the significance of the events among two biological conditions; (iii) graphical representation of user specific event with 3′-UTR annotation and read coverage on the 3′-UTR regions. APA-Scan is implemented in Python3. Source code and a comprehensive user’s manual are freely available at https://github.com/compbiolabucf/APA-Scan.

Result

APA-Scan was applied to both simulated and real RNA-seq datasets and compared with two widely used baselines DaPars and APAtrap. In simulation APA-Scan significantly improved the accuracy of 3′-UTR APA identification compared to the other baselines. The performance of APA-Scan was also validated by 3′-end-seq data and qPCR on mouse embryonic fibroblast cells. The experiments confirm that APA-Scan can detect unannotated 3′-UTR APA events and improve genome annotation.

Conclusion

APA-Scan is a comprehensive computational pipeline to detect transcriptome-wide 3′-UTR APA events. The pipeline integrates both RNA-seq and 3′-end-seq data information and can efficiently identify the significant events with a high-resolution short reads coverage plots.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04939-w.

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.

Shortening of the KHDRBS1 3'UTR by alternative cleavage and polyadenylation alters miRNA-mediated regulation and promotes gastric cancer progression

The shortening of the 3' untranslated regions (3'UTRs) due to alternative polyadenylation (APA) has become an important characteristic of cancer. However, the function of APA-induced 3'UTR shortening in gastric cancer (GC) remains unclear. KHDRBS1 (sam68), as an RNA-binding protein (RBP), is significantly upregulated in GC. In this study, we found that the 3'UTR of KHDRBS1 is generally shortened in GC tissues compared to paracancer tissues. Moreover, KHDRBS1 mRNA with a shortened 3'UTR can escape the inhibitory effect of miRNAs, resulting in its increased expression in GC. Overexpression of KHDRBS1, especially KHDRBS1 with a shortened 3'UTR, promotes the growth and metastasis of GC in vivo and in vitro. In conclusion, the experimental results show that shortening of the KHDRBS1 mRNA 3'UTR can mediate the overexpression of KHDRBS1 in GC cells and promote the progression of GC.

Gastrointestinal Biogeography of Luminal Microbiota and Short-Chain Fatty Acids in Sika Deer (Cervus nippon)

The gut microbiota of sika deer has been widely investigated, but the spatial distribution of symbiotic microbes among physical niches in the gastrointestinal tract remains to be established. While feces are the most commonly used biological samples in these studies, the accuracy of fecal matter as a proxy of the microbiome at other gastrointestinal sites is as yet unknown. In the present study, luminal contents obtained along the longitudinal axis of deer gastrointestinal tract (rumen, reticulum, omasum, abomasum, small intestine, cecum, colon, and rectum) were subjected to 16S rRNA gene sequencing for profiling of the microbial composition, and samples from the rumen, small intestine, and cecum were subjected to metabolomic analysis to evaluate short-chain fatty acid (SCFA) profiles. Prevotella bacteria were the dominant gastric core microbes, while Christensenellaceae_R-7_group was predominantly observed in the intestine. While the eight gastrointestinal sites displayed variations in microbial diversity, abundance, and function, they could be clustered into stomach, small intestine, and large intestine segments, and the results further highlighted a specific microbial niche of the small intestine. SCFA levels in the rumen, small intestine, and cecum were significantly different, with Bacteroidetes and Spirochaetes were shown to play a critical role in SCFA production. Finally, the rectal microbial composition was significantly correlated with colonic and cecum communities but not those of the small intestine and four gastric sites. Quantification of the compositions and biogeographic relationships between gut microbes and SCFAs in sika deer should provide valuable insights into the interactions contributing to microbial functions and metabolites. IMPORTANCE Feces or specific segments of the gastrointestinal tract (in particular, the rumen) were sampled to explore the gut microbiome. The gastrointestinal biogeography of the luminal microbiota in ruminants, which is critical to guide accurate sampling for different purposes, is poorly understood at present. The microbial community of the rectal sample (as a proxy of fecal sample) showed higher correlation with those of other large intestinal sites relative to the small intestine or stomach, suggesting that the microbial composition is specifically shaped by the unique physiological characteristics of different gastrointestinal niches. In addition, significant differences in microbiomes and SCFAs were observed among the different gastrointestinal sites.

Diverse cell-specific patterns of alternative polyadenylation in Drosophila

Most genes in higher eukaryotes express isoforms with distinct 3' untranslated regions (3' UTRs), generated by alternative polyadenylation (APA). Since 3' UTRs are predominant locations of post-transcriptional regulation, APA can render such programs conditional, and can also alter protein sequences via alternative last exon (ALE) isoforms. We previously used 3'-sequencing from diverse Drosophila samples to define multiple tissue-specific APA landscapes. Here, we exploit comprehensive single nucleus RNA-sequencing data (Fly Cell Atlas) to elucidate cell-type expression of 3' UTRs across >250 adult Drosophila cell types. We reveal the cellular bases of multiple tissue-specific APA/ALE programs, such as 3' UTR lengthening in differentiated neurons and 3' UTR shortening in spermatocytes and spermatids. We trace dynamic 3' UTR patterns across cell lineages, including in the male germline, and discover new APA patterns in the intestinal stem cell lineage. Finally, we correlate expression of RNA binding proteins (RBPs), miRNAs and global levels of cleavage and polyadenylation (CPA) factors in several cell types that exhibit characteristic APA landscapes, yielding candidate regulators of transcriptome complexity. These analyses provide a comprehensive foundation for future investigations of mechanisms and biological impacts of alternative 3' isoforms across the major cell types of this widely-studied model organism.

The Pleiotropy of PAX5 Gene Products and Function

PAX5, a member of the Paired Box (PAX) transcription factor family, is an essential factor for B-lineage identity during lymphoid differentiation. Mechanistically, PAX5 controls gene expression profiles, which are pivotal to cellular processes such as viability, proliferation, and differentiation. Given its crucial function in B-cell development, PAX5 aberrant expression also correlates with hallmark cancer processes leading to hematological and other types of cancer lesions. Despite the well-established association of PAX5 in the development, maintenance, and progression of cancer disease, the use of PAX5 as a cancer biomarker or therapeutic target has yet to be implemented. This may be partly due to the assortment of PAX5 expressed products, which layers the complexity of their function and role in various regulatory networks and biological processes. In this review, we provide an overview of the reported data describing PAX5 products, their regulation, and function in cellular processes, cellular biology, and neoplasm.

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

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

Potential of Establishing the Corresponding Human Microbial Community in Pseudo Germ-Free Mice through Fecal Microbe Transfer from Three Urolithin Metabotypes

Three urolithin metabotypes (UMs) have been defined in the population according to final urolithins converted by gut microbiota. Currently, it is difficult to establish the cause-and-effect relationship between urolithins and microbiota in human studies. Studies on the health effects of ellagic acid (EA) in animal models rarely consider the differences in the urolithin production. Therefore, the objective of this study is to establish human microbiota-associated (HMA) mice, imitating the microbiota composition of the three UMs. Antibiotic-induced pseudo germ-free mice were gavaged with fecal bacteria of the three UM donors for four weeks. The results showed that the ability to produce corresponding urolithins was successfully transferred from the donor of the three UMs to HMA mice. The three UM HMA mice adopted a humanized microbiota profile similar to their corresponding donor. The family Eggerthellaceae and genera Eggerthella and Gordonibacter were successfully transferred and colonized from UM-A/B donors to HMA mice. Overall, the three UM HMA mouse models were successfully established, which provide a basis for exploring the health effects of EA.

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.

Regulation of alternative polyadenylation by the C2H2-zinc-finger protein Sp1

Alternative polyadenylation (APA) enhances gene regulatory potential by increasing the diversity of mRNA transcripts. 3' UTR shortening through APA correlates with enhanced cellular proliferation and is a widespread phenomenon in tumor cells. Here, we show that the ubiquitously expressed transcription factor Sp1 binds RNA in vivo and is a common repressor of distal poly(A) site usage. RNA sequencing identified 2,344 genes (36% of the total mapped mRNA transcripts) with lengthened 3' UTRs upon Sp1 depletion. Sp1 preferentially binds the 3' UTRs of such lengthened transcripts and inhibits cleavage at distal sites by interacting with the subunits of the core cleavage and polyadenylation (CPA) machinery. The 3' UTR lengths of Sp1 target genes in breast cancer patient RNA-seq data correlate with Sp1 expression levels, implicating Sp1-mediated APA regulation in modulating tumorigenic properties. Taken together, our findings provide insights into the mechanism for dynamic APA regulation by unraveling a previously unknown function of the DNA-binding transcription factor Sp1.

Compartment-specific and ELAVL1-coordinated regulation of intronic polyadenylation isoforms by doxorubicin

Intronic polyadenylation (IPA) isoforms, which contain alternative last exons, are widely regulated in various biological processes and by many factors. However, little is known about their cytoplasmic regulation and translational status. In this study, we provide the first evidence that the genome-wide patterns of IPA isoform regulation during a biological process can be very distinct between the transcriptome and translatome, and between the nucleus and cytosol. Indeed, by 3'-seq analyses on breast cancer cells, we show that the genotoxic anticancer drug, doxorubicin, preferentially down-regulates the IPA to the last-exon (IPA:LE) isoform ratio in whole cells (as previously reported) but preferentially up-regulates it in polysomes. We further show that in nuclei, doxorubicin almost exclusively down-regulates the IPA:LE ratio, whereas in the cytosol, it preferentially up-regulates the isoform ratio, as in polysomes. Then, focusing on IPA isoforms that are up-regulated by doxorubicin in the cytosol and highly translated (up-regulated and/or abundant in polysomes), we identify several IPA isoforms that promote cell survival to doxorubicin. Mechanistically, by using an original approach of condition- and compartment-specific CLIP-seq (CCS-iCLIP) to analyze ELAVL1-RNA interactions in the nucleus and cytosol in the presence and absence of doxorubicin, as well as 3'-seq analyses upon ELAVL1 depletion, we show that the RNA-binding protein ELAVL1 mediates both nuclear down-regulation and cytosolic up-regulation of the IPA:LE isoform ratio in distinct sets of genes in response to doxorubicin. Altogether, these findings reveal differential regulation of the IPA:LE isoform ratio across subcellular compartments during drug response and its coordination by an RNA-binding protein.

CSTF2 Promotes Hepatocarcinogenesis and Hepatocellular Carcinoma Progression via Aerobic Glycolysis

Background

The shortening of 3’ untranslated regions (3’UTRs) of messenger RNAs(mRNAs) by alternative polyadenylation (APA) is an important mechanism for oncogene activation. Cleavage stimulation factor 2 (CSTF2), an important regulator of APA, has been reported to have a tumorigenic function in urothelial carcinoma of the bladder and lung cancers. However, the tumor-promoting role of CSTF2 in hepatocellular carcinoma (HCC) and its underlying molecular mechanism remains unclear.

Methods

Multiple databases were used to analyze the expression level and prognostic value of CSTF2 in HCC. Function enrichment analysis was used to investigate the molecular mechanism of CSTF2 for the occurrence and development of HCC. The biological function in HCC cell lines in vitro was determined by CCK8, colony formation, Transwell migration, and invasion assay. Moreover, the tumorigenic function of CSTF2 in vivo was measured by a subcutaneous tumor formation or injecting four plasmids into a mouse tail vein within 5–7 s in an immunocompetent HCC mouse model. In addition, aerobic glycolysis in HCC cells was determined by measuring the extracellular acid rate (ECAR) and extracellular glucose and lactate levels.

Results

Bioinformatics analysis revealed that CSTF2 was overexpressed in HCC tissues. The high expression of CSTF2 was correlated with a poor prognosis and high histological grades. CSTF2 knockout inhibited the proliferation, migration, and invasion of HCC cells. In addition, CSTF2 knockout HCC cells failed to form tumors by a subcutaneous graft experiment. Furthermore, endogenous CSTF2 knockout attenuated hepatocarcinogenesis in an immunocompetent HCC mouse model. Function enrichment analysis suggested that the high expression of CSTF2 was associated with enhanced glycolysis. Moreover, we found that CSTF2 knockout reduced the level of the short 3’ UTR isoform of hexokinase 2 and increased its level of long 3’UTR. Furthermore, CSTF2 knockout inhibited ECAR levels, glucose uptake, and lactate production.

Conclusion

Our results indicated that CSTF2 is highly expressed in HCC and is correlated with a poor prognosis and high histological grade. The knockout of CSTF2 inhibits the tumorigenesis and procession of HCC both in vitro and in vivo. Moreover, CSTF2 is associated with enhanced glycolysis. Therefore, this study suggests that CSTF2 might be a new prognostic biomarker and therapeutic target for HCC.

Intronic Polyadenylation in Acquired Cancer Drug Resistance Circumvented by Utilizing CRISPR/Cas9 with Homology-Directed Repair: The Tale of Human DNA Topoisomerase IIα

Simple Summary

DNA topoisomerase IIα (170 kDa, TOP2α/170) resolves nucleic acid topological entanglements by generating transient double-strand DNA breaks. TOP2α inhibitors/poisons stabilize TOP2α-DNA covalent complexes resulting in persistent DNA damage and are frequently utilized to treat a variety of cancers. Acquired resistance to these chemotherapeutic agents is often associated with decreased TOP2α/170 expression levels. Studies have demonstrated that a reduction in TOP2α/170 results from a type of alternative polyadenylation designated intronic polyadenylation (IPA). As a consequence of IPA, variant TOP2α mRNA transcripts have been characterized that have resulted in the translation of C-terminal truncated TOP2α isoforms with altered biological activities. In this paper, an example is discussed where circumvention of acquired TOP2α-mediated drug resistance was achieved by utilizing CRISPR/Cas9 specific gene editing of an exon/intron boundary through homology directed repair (HDR) to reduce TOP2α IPA. These results illustrate the therapeutic potential of CRISPR/Cas9/HDR to impact drug resistance associated with aberrant IPA.

Abstract

Intronic polyadenylation (IPA) plays a critical role in malignant transformation, development, progression, and cancer chemoresistance by contributing to transcriptome/proteome alterations. DNA topoisomerase IIα (170 kDa, TOP2α/170) is an established clinical target for anticancer agents whose efficacy is compromised by drug resistance often associated with a reduction of nuclear TOP2α/170 levels. In leukemia cell lines with acquired resistance to TOP2α-targeted drugs and reduced TOP2α/170 expression, variant TOP2α mRNA transcripts have been reported due to IPA that resulted in the translation of C-terminal truncated isoforms with altered nuclear-cytoplasmic distribution or heterodimerization with wild-type TOP2α/170. This review provides an overview of the various mechanisms regulating pre-mRNA processing and alternative polyadenylation, as well as the utilization of CRISPR/Cas9 specific gene editing through homology directed repair (HDR) to decrease IPA when splice sites are intrinsically weak or potentially mutated. The specific case of TOP2α exon 19/intron 19 splice site editing is discussed in etoposide-resistant human leukemia K562 cells as a tractable strategy to circumvent acquired TOP2α-mediated drug resistance. This example supports the importance of aberrant IPA in acquired drug resistance to TOP2α-targeted drugs. In addition, these results demonstrate the therapeutic potential of CRISPR/Cas9/HDR to impact drug resistance associated with aberrant splicing/polyadenylation.

Optimal CD8+ T cell effector function requires costimulation-induced RNA-binding proteins that reprogram the transcript isoform landscape

Boosting T cell activation through costimulation directs defense against cancer and viral infections. Despite multiple studies targeting costimulation in clinical trials, the increased potency and reprogramming of T cells endowed by costimulation is poorly understood. Canonical dogma states that transcription mediates T cell activation. Here, we show that the spliceosome, controlling post-transcriptional alternative splicing and alternative polyadenylation, is the most enriched pathway in T cells after CD134/CD137 costimulation. Costimulation of CD8+ T cells significantly increases expression of 29 RNA-binding proteins while RNA-seq uncovers over 1000 differential alternative splicing and polyadenylation events. Using in vivo mouse and in vitro human models, we demonstrate that RNA-binding protein Tardbp is required for effector cytokine production, CD8+ T cell clonal expansion, and isoform regulation after costimulation. The prospect of immune response optimization through reprogramming of mRNA isoform production offered herein opens new avenues for experimentally and therapeutically tuning the activities of T cells.

Neuronal Cells Display Distinct Stability Controls of Alternative Polyadenylation mRNA Isoforms, Long Non-Coding RNAs, and Mitochondrial RNAs

RNA stability plays an important role in gene expression. Here, using 3' end sequencing of newly made and pre-existing poly(A)+ RNAs, we compare transcript stability in multiple human cell lines, including HEK293T, HepG2, and SH-SY5Y. We show that while mRNA stability is generally conserved across the cell lines, specific transcripts having a high GC content and possibly more stable secondary RNA structures are relatively more stable in SH-SY5Y cells compared to the other 2 cell lines. These features also differentiate stability levels of alternative polyadenylation (APA) 3'UTR isoforms in a cell type-specific manner. Using differentiation of a neural stem cell line as a model, we show that mRNA stability difference could contribute to gene expression changes in neurogenesis and confirm the neuronal identity of SH-SY5Y cells at both gene expression and APA levels. In addition, compared to transcripts using 3'-most exon cleavage/polyadenylation sites (PASs), those using intronic PASs are generally less stable, especially when the PAS-containing intron is large and has a strong 5' splice site, suggesting that intronic polyadenylation mostly plays a negative role in gene expression. Interestingly, the differential mRNA stability among APA isoforms appears to buffer PAS choice in these cell lines. Moreover, we found that several other poly(A)+ RNA species, including promoter-associated long noncoding RNAs and transcripts encoded by the mitochondrial genome, are more stable in SH-SY5Y cells than the other 2 cell lines, further highlighting distinct RNA metabolism in neuronal cells. Together, our results indicate that distinct RNA stability control in neuronal cells may contribute to the gene expression and APA programs that define their cell identity.

Signatures of mRNA Alternative Polyadenylation in Arabidopsis Leaf Development

Alternative polyadenylation (APA) of pre-mRNA is an important co-transcriptional mechanism that modulates gene expression, leading to transcriptomic and functional diversities. The role of APA in Arabidopsis leaf development, however, remains elusive. We applied a poly(A)-tag sequencing (PAT-seq) technique to characterize APA-mediated regulation events in cotyledon and in five stages of true leaf development. Over 60% APA was identified in genes expressed in leaves, consistent with the results in previous publications. However, a reduced APA level was detected in younger leaves, reaching 44% in the 18th true leaf. Importantly, we also found that >70% of the poly(A) site usages were altered in the second true leaf relative to the cotyledon. Compared with the cotyledon, more genes in the second true leaf tended to use the distal site of 3′UTR, but this was not found in pairwise comparison among other true leaves. In addition, a significant APA gene was found to be decreased in a pairwise comparison among true leaves, including differentially expressed genes. The APA genes identified herein were associated with specific biological processes, including metabolic and cellular processes and response to stimuli and hormones. These results provide a new insight into the regulation of Arabidopsis leaf development through APA.

Precise gene models using long-read sequencing reveal a unique poly(A) signal in Giardia lamblia

During pre-mRNA processing, the poly(A) signal is recognized by a protein complex that ensures precise cleavage and polyadenylation of the nascent transcript. The location of this cleavage event establishes the length and sequence of the 3' UTR of an mRNA, thus determining much of its post-transcriptional fate. Using long-read sequencing, we characterize the polyadenylation signal and related sequences surrounding Giardia lamblia cleavage sites for over 2600 genes. We find that G. lamblia uses an AGURAA poly(A) signal, which differs from the mammalian AAUAAA. We also describe how G. lamblia lacks common auxiliary elements found in other eukaryotes, along with the proteins that recognize them. Further, we identify 133 genes with evidence of alternative polyadenylation. These results suggest that despite pared-down cleavage and polyadenylation machinery, 3' end formation still appears to be an important regulatory step for gene expression in G. lamblia.

Template-Independent Poly(A)-Tail Decay and RNASEL as Potential Cellular Biomarkers for Prostate Cancer Development

Simple Summary

The ultimate need in cancer tissue is to adapt translation machinery to accelerated protein synthesis in a rapidly proliferating environment. Our study was designed with the aim of integrating fundamental and clinical research to find new biomarkers for prostate cancer (PC) with clinical usefulness for the stratification prediction of healthy tissue transition into malignant phenotype. This study revealed: (i) an entirely novel mechanism of the regulatory influence of Poly(A) deadenylase in mRNAs translational activity and the 3′ mRNA untranslated region (3′UTR) length in cancer tissue and its regulation by the poly(A) decay; (ii) the RNASEL interrelationship with the inflammatory pattern of PC and corresponding tumor-adjacent and healthy tissue; and (iii) the sensitivity, specificity, and predictive value of these enzymes. The proposed manuscript is based on the use of specific biochemical and immunoassay methods with the principal research adapted for the use of tissue specimens.

Abstract

The post-transcriptional messenger RNA (mRNA) decay and turnover rate of the template-independent poly(A) tail, localized at the 3′-untranslated region (3′UTR) of mRNA, have been documented among subtle mechanisms of uncontrolled cancer tissue growth. The activity of Poly(A) deadenylase and the expression pattern of RNASEL have been examined. A total of 138 prostate tissue specimens from 46 PC patients (cancer specimens, corresponding adjacent surgically healthy tissues, and in their normal counterparts, at least 2 cm from carcinoma) were used. For the stratification prediction of healthy tissue transition into malignant phenotype, the enzyme activity of tumor-adjacent tissue was considered in relation to the presence of microfocal carcinoma. More than a four-times increase in specific enzyme activity (U/L g.prot) was registered in PC on account of both the dissociation of its inhibitor and genome reprogramming. The obtained ROC curve and Youden index showed that Poly(A) deadenylase identified PC with a sensitivity of 93.5% and a specificity of 94.6%. The RNASEL expression profile was raised significantly in PC, but the sensitivity was 40.5% and specificity was 86.9%. A significantly negative correlation between PC and control tissue counterparts with a higher expression pattern in lymphocyte-infiltrated samples were reported. In conclusion, significantly upregulated Poly(A) deadenylase activity may be a checkpoint for the transition of precancerous lesion to malignancy, while RNASEL may predict chronic inflammation.

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.

RNA-Seq Reveals miRNA and mRNA Co-regulate Muscle Differentiation in Fetal Leizhou Goats

Muscle differentiation is an essential link in animal growth and development, and microRNA and mRNA are indispensable in skeletal muscle differentiation. To improve the meat quality and production of the Leizhou goat, it is vital to understand the molecular mechanism by which its skeletal muscle differentiates. By RNA sequencing (RNA-SEQ), we established miRNA-mRNA profiles of Leizhou goats at three stages: fetal day 70, 90, and 120. There were 991 differently expressed mRNAs and 39 differentially expressed miRNAs found, with the differentially expressed mRNAs mainly enriched in calcium ion binding, ECM-receptor interaction, and Focal adhesion. CKM and MYH3, two muscle differentiation markers, were significantly differentially expressed during this period. In addition, we found that chi-miR-129-5p, chi-miR-433, and chi-miR-24-3p co-regulate muscle differentiation with their target genes. Finally, we can confirm that muscle differentiation occurred in Leizhou goat between 90 and 120 days of the fetus. This study is helpful to better explore the molecular mechanism of goat muscle differentiation.

Sorting and packaging of RNA into extracellular vesicles shape intracellular transcript levels

BACKGROUND

Extracellular vesicles (EVs) are released by nearly every cell type and have attracted much attention for their ability to transfer protein and diverse RNA species from donor to recipient cells. Much attention has been given so far to the features of EV short RNAs such as miRNAs. However, while the presence of mRNA and long noncoding RNA (lncRNA) transcripts in EVs has also been reported by multiple different groups, the properties and function of these longer transcripts have been less thoroughly explored than EV miRNA. Additionally, the impact of EV export on the transcriptome of exporting cells has remained almost completely unexamined. Here, we globally investigate mRNA and lncRNA transcripts in endothelial EVs in multiple different conditions.

RESULTS

In basal conditions, long RNA transcripts enriched in EVs have longer than average half-lives and distinctive stability-related sequence and structure characteristics including shorter transcript length, higher exon density, and fewer 3' UTR A/U-rich elements. EV-enriched long RNA transcripts are also enriched in HNRNPA2B1 binding motifs and are impacted by HNRNPA2B1 depletion, implicating this RNA-binding protein in the sorting of long RNA to EVs. After signaling-dependent modification of the cellular transcriptome, we observed that, unexpectedly, the rate of EV enrichment relative to cells was altered for many mRNA and lncRNA transcripts. This change in EV enrichment was negatively correlated with intracellular abundance, with transcripts whose export to EVs increased showing decreased abundance in cells and vice versa. Correspondingly, after treatment with inhibitors of EV secretion, levels of mRNA and lncRNA transcripts that are normally highly exported to EVs increased in cells, indicating a measurable impact of EV export on the long RNA transcriptome of the exporting cells. Compounds with different mechanisms of inhibition of EV secretion affected the cellular transcriptome differently, suggesting the existence of multiple EV subtypes with different long RNA profiles.

CONCLUSIONS

We present evidence for an impact of EV physiology on the characteristics of EV-producing cell transcriptomes. Our work suggests a new paradigm in which the sorting and packaging of transcripts into EVs participate, together with transcription and RNA decay, in controlling RNA homeostasis and shape the cellular long RNA abundance profile.

PABP prevents the untimely decay of select mRNA populations in human cells

Gene expression is tightly regulated at the levels of both mRNA translation and stability. The poly(A)-binding protein (PABP) is thought to play a role in regulating these processes by binding the mRNA 3' poly(A) tail and interacting with both the translation and mRNA deadenylation machineries. In this study, we directly investigate the impact of PABP on translation and stability of endogenous mRNAs in human cells. Remarkably, our transcriptome-wide analysis only detects marginal mRNA translation changes in PABP-depleted cells. In contrast, rapidly depleting PABP alters mRNA abundance and stability, albeit non-uniformly. Otherwise stable transcripts, including those encoding proteins with constitutive functions, are destabilized in PABP-depleted cells. In contrast, many unstable mRNAs, including those encoding proteins with regulatory functions, decay at similar rates in presence or absence of PABP. Moreover, PABP depletion-induced cell death can partially be suppressed by disrupting the mRNA decapping and 5'-3' decay machinery. Finally, we provide evidence that the LSM1-7 complex promotes decay of "stable" mRNAs in PABP-depleted cells. Taken together, these findings suggest that PABP plays an important role in preventing the untimely decay of select mRNA populations.

An Evolutionarily Conserved AU-Rich Element in the 3' Untranslated Region of a Transcript Misannotated as a Long Noncoding RNA Regulates RNA Stability

One of the primary mechanisms of post-transcriptional gene regulation is the modulation of RNA stability. We recently discovered that LINC00675, a transcript annotated as a long noncoding RNA (lncRNA), is transcriptionally regulated by FOXA1 and encodes a highly conserved small protein that localizes to the endoplasmic reticulum, hence renamed as FORCP (FOXA1-regulated conserved small protein). Here, we show that the endogenous FORCP transcript is rapidly degraded and rendered unstable as a result of 3'UTR-mediated degradation. Surprisingly, although the FORCP transcript is a canonical nonsense-mediated decay (NMD) and microRNA (miRNA) target, we found that it is not degraded by NMD or miRNAs. Targeted deletion of an evolutionarily conserved region in the FORCP 3'UTR using CRISPR/Cas9 significantly increased the stability of the FORCP transcript. Interestingly, this region requires the presence of an immediate downstream 55-nt-long sequence for transcript stability regulation. Functionally, colorectal cancer cells lacking this conserved region expressed from the endogenous FORCP locus displayed decreased proliferation and clonogenicity. These data demonstrate that the FORCP transcript is destabilized via conserved elements within its 3'UTR and emphasize the need to interrogate the function of a given 3'UTR in its native context.

CFIm-mediated alternative polyadenylation remodels cellular signaling and miRNA biogenesis

The mammalian cleavage factor I (CFIm) has been implicated in alternative polyadenylation (APA) in a broad range of contexts, from cancers to learning deficits and parasite infections. To determine how the CFIm expression levels are translated into these diverse phenotypes, we carried out a multi-omics analysis of cell lines in which the CFIm25 (NUDT21) or CFIm68 (CPSF6) subunits were either repressed by siRNA-mediated knockdown or over-expressed from stably integrated constructs. We established that >800 genes undergo coherent APA in response to changes in CFIm levels, and they cluster in distinct functional classes related to protein metabolism. The activity of the ERK pathway traces the CFIm concentration, and explains some of the fluctuations in cell growth and metabolism that are observed upon CFIm perturbations. Furthermore, multiple transcripts encoding proteins from the miRNA pathway are targets of CFIm-dependent APA. This leads to an increased biogenesis and repressive activity of miRNAs at the same time as some 3′ UTRs become shorter and presumably less sensitive to miRNA-mediated repression. Our study provides a first systematic assessment of a core set of APA targets that respond coherently to changes in CFIm protein subunit levels (CFIm25/CFIm68). We describe the elicited signaling pathways downstream of CFIm, which improve our understanding of the key role of CFIm in integrating RNA processing with other cellular activities.

MCL1 alternative polyadenylation is essential for cell survival and mitochondria morphology

Alternative polyadenylation in the 3' UTR (3' UTR-APA) is a mode of gene expression regulation, fundamental for mRNA stability, translation and localization. In the immune system, it was shown that upon T cell activation, there is an increase in the relative expression of mRNA isoforms with short 3' UTRs resulting from 3' UTR-APA. However, the functional significance of 3' UTR-APA remains largely unknown. Here, we studied the physiological function of 3' UTR-APA in the regulation of Myeloid Cell Leukemia 1 (MCL1), an anti-apoptotic member of the Bcl-2 family essential for T cell survival. We found that T cells produce two MCL1 mRNA isoforms (pA1 and pA2) by 3' UTR-APA. We show that upon T cell activation, there is an increase in both the shorter pA1 mRNA isoform and MCL1 protein levels. Moreover, the less efficiently translated pA2 isoform is downregulated by miR-17, which is also more expressed upon T cell activation. Therefore, by increasing the expression of the more efficiently translated pA1 mRNA isoform, which escapes regulation by miR-17, 3' UTR-APA fine tunes MCL1 protein levels, critical for activated T cells' survival. Furthermore, using CRISPR/Cas9-edited cells, we show that depletion of either pA1 or pA2 mRNA isoforms causes severe defects in mitochondria morphology, increases apoptosis and impacts cell proliferation. Collectively, our results show that MCL1 alternative polyadenylation has a key role in the regulation of MCL1 protein levels upon T cell activation and reveal an essential function for MCL1 3' UTR-APA in cell viability and mitochondria dynamics.

Differential Effects on the Translation of Immune-Related Alternatively Polyadenylated mRNAs in Melanoma and T Cells by eIF4A Inhibition

Targeting the translation initiation complex eIF4F, which binds the 5' cap of mRNAs, is a promising anti-cancer approach. Silvestrol, a small molecule inhibitor of eIF4A, the RNA helicase component of eIF4F, inhibits the translation of the mRNA encoding the signal transducer and activator of transcription 1 (STAT1) transcription factor, which, in turn, reduces the transcription of the gene encoding one of the major immune checkpoint proteins, i.e., programmed death ligand-1 (PD-L1) in melanoma cells. A large proportion of human genes produce multiple mRNAs differing in their 3'-ends through the use of alternative polyadenylation (APA) sites, which, when located in alternative last exons, can generate protein isoforms, as in the STAT1 gene. Here, we provide evidence that the STAT1α, but not STAT1β protein isoform generated by APA, is required for silvestrol-dependent inhibition of PD-L1 expression in interferon-γ-treated melanoma cells. Using polysome profiling in activated T cells we find that, beyond STAT1, eIF4A inhibition downregulates the translation of some important immune-related mRNAs, such as the ones encoding TIM-3, LAG-3, IDO1, CD27 or CD137, but with little effect on the ones for BTLA and ADAR-1 and no effect on the ones encoding CTLA-4, PD-1 and CD40-L. We next apply RT-qPCR and 3'-seq (RNA-seq focused on mRNA 3' ends) on polysomal RNAs to analyze in a high throughput manner the effect of eIF4A inhibition on the translation of APA isoforms. We identify about 150 genes, including TIM-3, LAG-3, AHNAK and SEMA4D, for which silvestrol differentially inhibits the translation of APA isoforms in T cells. It is therefore crucial to consider 3'-end mRNA heterogeneity in the understanding of the anti-tumor activities of eIF4A inhibitors.

Regulation of the Alternative Neural Transcriptome by ELAV/Hu RNA Binding Proteins

The process of alternative polyadenylation (APA) generates multiple 3' UTR isoforms for a given locus, which can alter regulatory capacity and on occasion change coding potential. APA was initially characterized for a few genes, but in the past decade, has been found to be the rule for metazoan genes. While numerous differences in APA profiles have been catalogued across genetic conditions, perturbations, and diseases, our knowledge of APA mechanisms and biology is far from complete. In this review, we highlight recent findings regarding the role of the conserved ELAV/Hu family of RNA binding proteins (RBPs) in generating the broad landscape of lengthened 3' UTRs that is characteristic of neurons. We relate this to their established roles in alternative splicing, and summarize ongoing directions that will further elucidate the molecular strategies for neural APA, the in vivo functions of ELAV/Hu RBPs, and the phenotypic consequences of these regulatory paradigms in neurons.

Cooperation of RNA-Binding Proteins – a Focus on Roquin Function in T Cells

Post-transcriptional gene regulation by RNA-binding proteins (RBPs) is important in the prevention of inflammatory and autoimmune diseases. With respect to T cell activation and differentiation, the RBPs Roquin-1/2 and Regnase-1 play pivotal roles by inducing degradation and/or translational silencing of target mRNAs. These targets encode important proinflammatory mediators and thus Roquin and Regnase-1 functions dampen cellular programs that can lead to inflammation and autoimmune disease. Recent findings demonstrate direct physical interaction of both RBPs. Here, we propose that cooperativity of trans-acting factors may be more generally used to reinforce the regulatory impact on selected targets and promote specific cell fate decisions. We develop this concept for Roquin and Regnase-1 function in resting and activated T cells and discuss the involvement in autoimmunity as well as how the therapeutic potential can be used in anti-tumor therapies.

MiR-138 is a potent regulator of the heterogenous MYC transcript population in cancers

3'UTR shortening in cancer has been shown to activate oncogenes, partly through the loss of microRNA-mediated repression. This suggests that many reported microRNA-oncogene target interactions may not be present in cancer cells. One of the most well-studied oncogenes is the transcription factor MYC, which is overexpressed in more than half of all cancers. MYC overexpression is not always accompanied by underlying genetic aberrations. In this study, we demonstrate that the MYC 3'UTR is shortened in colorectal cancer (CRC). Using unbiased computational and experimental approaches, we identify and validate microRNAs that target the MYC coding region. In particular, we show that miR-138 inhibits MYC expression and suppresses tumor growth of CRC and hepatocellular carcinoma (HCC) cell lines. Critically, the intravenous administration of miR-138 significantly impedes MYC-driven tumor growth in vivo. Taken together, our results highlight the previously uncharacterized shortening of the MYC 3'UTR in cancer, and identify miR-138 as a potent regulator of the heterogenous MYC transcript population.

Alternative polyadenylation dysregulation contributes to the differentiation block of acute myeloid leukemia

Posttranscriptional regulation has emerged as a driver for leukemia development and an avenue for therapeutic targeting. Among posttranscriptional processes, alternative polyadenylation (APA) is globally dysregulated across cancer types. However, limited studies have focused on the prevalence and role of APA in myeloid leukemia. Furthermore, it is poorly understood how altered poly(A) site usage of individual genes contributes to malignancy or whether targeting global APA patterns might alter oncogenic potential. In this study, we examined global APA dysregulation in patients with acute myeloid leukemia (AML) by performing 3' region extraction and deep sequencing (3'READS) on a subset of AML patient samples along with healthy hematopoietic stem and progenitor cells (HSPCs) and by analyzing publicly available data from a broad AML patient cohort. We show that patient cells exhibit global 3' untranslated region (UTR) shortening and coding sequence lengthening due to differences in poly(A) site (PAS) usage. Among APA regulators, expression of FIP1L1, one of the core cleavage and polyadenylation factors, correlated with the degree of APA dysregulation in our 3'READS data set. Targeting global APA by FIP1L1 knockdown reversed the global trends seen in patients. Importantly, FIP1L1 knockdown induced differentiation of t(8;21) cells by promoting 3'UTR lengthening and downregulation of the fusion oncoprotein AML1-ETO. In non-t(8;21) cells, FIP1L1 knockdown also promoted differentiation by attenuating mechanistic target of rapamycin complex 1 (mTORC1) signaling and reducing MYC protein levels. Our study provides mechanistic insights into the role of APA in AML pathogenesis and indicates that targeting global APA patterns can overcome the differentiation block in patients with AML.

The Role of Alternative Polyadenylation in the Regulation of Subcellular RNA Localization

Alternative polyadenylation (APA) is a widespread and conserved regulatory mechanism that generates diverse 3' ends on mRNA. APA patterns are often tissue specific and play an important role in cellular processes such as cell proliferation, differentiation, and response to stress. Many APA sites are found in 3' UTRs, generating mRNA isoforms with different 3' UTR contents. These alternate 3' UTR isoforms can change how the transcript is regulated, affecting its stability and translation. Since the subcellular localization of a transcript is often regulated by 3' UTR sequences, this implies that APA can also change transcript location. However, this connection between APA and RNA localization has only recently been explored. In this review, we discuss the role of APA in mRNA localization across distinct subcellular compartments. We also discuss current challenges and future advancements that will aid our understanding of how APA affects RNA localization and molecular mechanisms that drive these processes.

Alternative polyadenylation by sequential activation of distal and proximal PolyA sites

Analogous to alternative splicing, alternative polyadenylation (APA) has long been thought to occur independently at proximal and distal polyA sites. Using fractionation-seq, we unexpectedly identified several hundred APA genes in human cells whose distal polyA isoforms are retained in chromatin/nuclear matrix and whose proximal polyA isoforms are released into the cytoplasm. Global metabolic PAS-seq and Nanopore long-read RNA-sequencing provide further evidence that the strong distal polyA sites are processed first and the resulting transcripts are subsequently anchored in chromatin/nuclear matrix to serve as precursors for further processing at proximal polyA sites. Inserting an autocleavable ribozyme between the proximal and distal polyA sites, coupled with a Cleave-seq approach that we describe here, confirms that the distal polyA isoform is indeed the precursor to the proximal polyA isoform. Therefore, unlike alternative splicing, APA sites are recognized independently, and in many cases, in a sequential manner. This provides a versatile strategy to regulate gene expression in mammalian cells.

CRISPR-Mediated Knockout of Long 3' UTR mRNA Isoforms in mESC-Derived Neurons

Alternative cleavage and polyadenylation (APA) is pervasive, occurring for more than 70% of human and mouse genes. Distal poly(A) site selection to generate longer 3′ UTR mRNA isoforms is prevalent in the nervous system, affecting thousands of genes. Here, we establish mouse embryonic stem cell (mESC)-derived neurons (mES-neurons) as a suitable system to study long 3′ UTR isoforms. RNA-seq analysis revealed that mES-neurons show widespread 3′ UTR lengthening that closely resembles APA patterns found in mouse cortex. mESCs are highly amenable to genetic manipulation. We present a method to eliminate long 3′ UTR isoform expression using CRISPR/Cas9 editing. This approach can lead to clones with the desired deletion within several weeks. We demonstrate this strategy on the Mprip gene as a proof-of-principle. To confirm loss of long 3′ UTR expression and the absence of cryptic poly(A) site usage stemming from the CRISPR deletion, we present a simple and cost-efficient targeted long-read RNA-sequencing strategy using the Oxford Nanopore Technologies platform. Using this method, we confirmed specific loss of the Mprip long 3′ UTR isoform. CRISPR gene editing of mESCs thus serves as a highly relevant platform for studying the molecular and cellular functions of long 3′ UTR mRNA isoforms.

Predicting sequence and structural features of effective piRNA target binding sites

Piwi-interacting RNA (piRNA) targets are usually identified through base pairing between the piRNA seed region and complementary bases on the target mRNAs, which often results in false predictions. Crosslinking immunoprecipitation (CLIP) study emerges as a promising method that enables accurate identification of PIWI-clade-based targets containing RNA-binding sites. In the present study, we have analyzed the piRNA-target CLIP-seq datasets to uncover the additional characteristic features of piRNA targets. We studied important sequence and structural features using IP⁺ and IP^- set targets that might enhance the accuracy of target site predictions. Analysis has revealed substantial enrichment of AU in target sites as well as in and around the 30 nts upstream and downstream of target sites in IP⁺ set relative to IP^- set that might be contributing to lowering the minimal folding energy of target sites of IP⁺ + set that might be easing the base pairing between piRNA and their targets. We have also found a lower MFE threshold (en) and higher miRanda score for piRNA targets. Interestingly, we have found that majority of the target sites are residing within 3'UTR, suggesting 3'UTR as a preferential target site like that of miRNA targets. Thus, we hypothesize that our findings on additional key features of piRNA target sites might be valuable in identifying the potential targets of piRNA accurately, which will aid in decrypting their functional importance. This article is protected by copyright. All rights reserved.

The androgen receptor messenger RNA: what do we know?

The Androgen Receptor (AR), transcriptionally activated by its ligands, testosterone and dihydrotestosterone (DHT), is widely expressed in cells and tissues, influencing normal biology and disease states. The protein product of the AR gene is involved in the regulation of numerous biological functions, including the development and maintenance of the normal prostate gland and of the cardiovascular, musculoskeletal and immune systems. Androgen signalling, mediated by AR protein, plays a crucial role in the development of prostate cancer (PCa), and is presumed to be involved in other cancers including those of the breast, bladder, liver and kidney. Significant research and reviews have focused on AR protein function; however, inadequate research and literature exist to define the function of AR mRNA in normal and cancer cells. The AR mRNA transcript is nearly 11 Kb long and contains a long 3’ untranslated region (UTR), suggesting its biological role in post-transcriptional regulation, consequently affecting the overall functions of both normal and cancer cells. Research has demonstrated that many biological activities, including RNA stability, translation, cellular trafficking and localization, are associated with the 3’ UTRs of mRNAs. In this review, we describe the potential role of the AR 3’ UTR and summarize RNA-binding proteins (RBPs) that interact with the AR mRNA to regulate post-transcriptional metabolism. We highlight the importance of AR mRNA as a critical modulator of carcinogenesis and its important role in developing therapy-resistant prostate cancer.

Regulation of Alternative Splicing during Epithelial-Mesenchymal Transition

Alternative splicing is an essential mechanism of gene regulation, giving rise to remarkable protein diversity in higher eukaryotes. Epithelial-mesenchymal transition (EMT) is a developmental process that plays an essential role in metazoan embryogenesis. Recent studies have revealed that alternative splicing serves as a fundamental layer of regulation that governs cells to undergo EMT. In this review, we summarize recent findings on the functional impact of alternative splicing in EMT and EMT-associated activities. We then discuss the regulatory mechanisms that control alternative splicing changes during EMT.

Computational Detection of MicroRNA Targets

MicroRNAs (miRNAs) are small noncoding RNAs that are recognized as posttranscriptional regulators of gene expression. These molecules have been shown to play important roles in several cellular processes. MiRNAs act on their target by guiding the RISC complex and binding to the mRNA molecule. Thus, it is recognized that the function of a miRNA is determined by the function of its target (s). By using high-throughput methodologies, novel miRNAs are being identified, but their functions remain uncharted. Target validation is crucial to properly understand the specific role of a miRNA in a cellular pathway. However, molecular techniques for experimental validation of miRNA-target interaction are expensive, time-consuming, laborious, and can be not accurate in inferring true interactions. Thus, accurate miRNA target predictions are helpful to understand the functions of miRNAs. There are several algorithms proposed for target prediction and databases containing miRNA-target information. However, these available computational tools for prediction still generate a large number of false positives and fail to detect a considerable number of true targets, which indicates the necessity of highly confident approaches to identify bona fide miRNA-target interactions. This chapter focuses on tools and strategies used for miRNA target prediction, by providing practical insights and outlooks.

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.

MicroRNAs in Genetic Etiology of Human Diseases

Since their first discovery more than 20 years ago, miRNAs have been subject to deliberate research and analysis for revealing their physiological or pathological involvement. Regulatory roles of miRNAs in signal transduction, gene expression, and cellular processes in development, differentiation, proliferation, apoptosis, and homeostasis also imply their critical role in disease pathogenesis. Their roles in cancer, neurodegenerative diseases, and other systemic diseases have been studied broadly. In these regulatory pathways, their mutations and target sequence variations play critical roles to determine their functional repertoire. In this chapter, we summarize studies that investigated the role of mutations, polymorphisms, and other variations of miRNAs in respect to pathological processes.

Pax-5 Protein Expression Is Regulated by Transcriptional 3'UTR Editing

The Pax-5 gene encodes a transcription factor that is essential for B-cell commitment and maturation. However, Pax-5 deregulation is associated with various cancer lesions, notably hematopoietic cancers. Mechanistically, studies have characterized genetic alterations within the Pax-5 locus that result in either dominant oncogenic function or haploinsufficiency-inducing mutations leading to oncogenesis. Apart from these mutations, some examples of aberrant Pax-5 expression cannot be associated with genetic alterations. In the present study, we set out to elucidate potential alterations in post-transcriptional regulation of Pax-5 expression and establish that Pax-5 transcript editing represents an important means to aberrant expression. Upon the profiling of Pax-5 mRNA in leukemic cells, we found that the 3′end of the Pax-5 transcript is submitted to alternative polyadenylation (APA) and alternative splicing events. Using rapid amplification of cDNA ends (3′RACE) from polysomal fractions, we found that Pax-5 3′ untranslated region (UTR) shortening correlates with increased ribosomal occupancy for translation. These observations were also validated using reporter gene assays with truncated 3′UTR regions cloned downstream of a luciferase gene. We also showed that Pax-5 3′UTR editing has direct repercussions on regulatory elements such as miRNAs, which in turn impact Pax-5 protein expression. More importantly, we found that advanced staging of various hematopoietic cancer lesions relates to shorter Pax-5 3′UTRs. Altogether, our findings identify novel molecular mechanisms that account for aberrant expression and function of the Pax-5 oncogene in cancer cells. These findings also present new avenues for strategic intervention in Pax-5-mediated cancers.