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

RNA processing in innate immunity: regulation by RNA-binding proteins

RNA processing is an important but often overlooked process in determining protein expression. Alternative polyadenylation and regulated mRNA stability control the amount and duration of protein expression, while alternative splicing also controls protein identity and function. Much work in innate immunity has focused on the activation of transcription factors and the downstream consequences in gene expression. However, there is increasing evidence indicating that regulation of RNA processing by RNA-binding proteins (RBPs) contributes significantly to tuning the innate immune response. Herein we review work highlighting the impact of RNA processing in innate immunity and describe the RBPs and mechanisms driving this regulation. We conclude with a discussion of unanswered questions to motivate continued research in this important and understudied field.

Post-transcriptional regulation of diabetic wound healing by junctional adhesion molecule A/miR-106b axis

BACKGROUND

Despite advancements in molecular science and biomaterial technology, the mechanisms underlying the impaired healing of diabetic wounds remain unclear. In this study, we investigated the post-transcriptional regulation of diabetic wound healing using JAM-A.

METHODS

Mouse wound models, hematoxylin and eosin staining analysis, and scratch wound assays were used to investigate the effects of JAM-A 3'-UTR on the re-epithelialization of diabetic wounds, whereas RNA pulldown, microRNA-seq, and bioinformatics analyses were performed to identify key miRNA players and predict their target genes. In situ hybridization, immunohistochemistry, western blotting, and polymerase chain reaction (PCR) were used to confirm the alternative splicing of JAM-A 3'-UTR in diabetic conditions. CCK-8 proliferation assays, scratch wound assays, PCR, western blotting, and dual-luciferase assays were performed to study the changes in cell proliferation and migration induced by miR-106b-5p modification and confirm the target gene.

RESULTS

JAM-A 3'-UTR accelerated re-epithelialization in diabetic mouse wounds. Shortened splicing was found in the 3'-UTR of JAM-A under diabetic conditions, leading to the excessive release of miR-106b-5p while the promoter of miR-106b was activated. Furthermore, upregulated miR-106b-5p over-activated cell proliferation and inhibited cell migration in diabetic wound keratinocytes by suppressing the target gene PTEN/TIAM1 and regulating the AKT and RAC1 pathways, thereby impairing wound re-epithelialization.

CONCLUSIONS

We identified alternative splicing of JAM-A 3'-UTR in diabetic conditions, which caused the excessive release of miR-106b. Upregulation of miR-106b reduced the expression of its target genes, PTEN and TIAM1, which led to hyperactive proliferation and impaired migration of keratinocytes, thereby dysregulating wound re-epithelialization.

The Dynamic Landscape of 3'-UTR Alternative Polyadenylation Across Mouse Fetal Development and Anatomy

Alternative cleavage and polyadenylation (APA) in the 3'-untranslated region (3'-UTR) of mRNA produces transcripts with varied 3'-UTR and plays a key role in development and organogenesis. This work characterizes 3'-UTR APA using 85 high-quality RNA-seq datasets encompassing 12 tissue types and eight developmental stages of mouse fetuses. Results show that 46.4% of expressed genes undergo APA in a tissue-specific manner. Changes in polyadenylation site (pAS) usage often operate beyond transcriptional control, revealing APA as an additional layer of gene regulation. Sequence analysis demonstrates that pAS selection, governed by polyadenylation signal strength and adenine preferences, is evolutionarily conserved between mice and humans. Intriguingly, brain tissues display complex 3'-UTR APA dynamics during development, potentially regulated by RNA-binding proteins such as Rbm38, potentially impacting 3' UTR extension by restricting distal pAS usage. These APA events are associated with a depletion of conserved miRNA binding sites and an enrichment of transposable elements within alternative 3' UTRs. To facilitate further research, this work develops APApedia (http://xozhanglab.com/apapedia/), a comprehensive database cataloging identified 3'-UTR APA events, which serves as a valuable resource for the community to study APA in development and tissue-specific regulation. This comprehensive resource aids in deciphering the functional implications of APA in mouse fetal development.

Srsf3-Dependent APA Drives Macrophage Maturation and Limits Atherosclerosis

BACKGROUND

Circulating monocytes largely contribute to macrophage buildup in atheromata, which is crucial for clearing subendothelial LDLs (low-density lipoproteins) and dead cells; however, the transitional trajectory from monocytes to macrophages in atherosclerotic plaques and the underlying regulatory mechanism remain unclear. Moreover, the role of alternative polyadenylation, a posttranscriptional regulator of cell fate, in monocyte/macrophage fate decisions during atherogenesis is not entirely understood.

METHODS

To identify monocyte/macrophage subtypes in atherosclerotic lesions and the effect of alternative polyadenylation on these subtypes and atherogenesis, single-cell RNA sequencing, 3'-end sequencing, flow cytometric, and histopathologic analyses were performed on plaques obtained from Apoe-/- mouse arteries with or without myeloid deletion of Srsf3 (serine/arginine-rich splicing factor 3). Cell fractionation, polysome profiling, L-azidohomoalanine metabolic labeling assay, and metabolomic profiling were conducted to disclose the underlying mechanisms. Reprogramming of widespread alternative polyadenylation patterns was estimated in human plaques via bulk RNA sequencing.

RESULTS

We identified a subset of lesional cells in a monocyte-to-macrophage transitional state, which exhibited high expression of chemokines in mice. Srsf3 deletion caused a maturation delay of these transitional cells and phagocytic impairment of lesional macrophages, aggravating atherosclerosis. Mechanistically, Srsf3 deficiency shortened 3' untranslated regions of mitochondria-associated Aars2 (alanyl-tRNA synthetase 2), disrupting its translation. The resultant impairment of protein synthesis in mitochondria led to mitochondrial dysfunction with declined NAD+ (nicotinamide adenine dinucleotide, oxidized form) levels, activation of the integrated stress response, and metabolic reprogramming in macrophages. Administering an NAD+ precursor nicotinamide mononucleotide or the integrated stress response inhibitor partially restored Srsf3-deficient macrophage maturation, and nicotinamide mononucleotide treatment mitigated the proatherosclerotic effects of Srsf3 deficiency. Consistently, Srsf3 downregulation, global 3' untranslated region shortening, and accumulation of these transitional macrophages were associated with atherosclerosis progression in humans.

CONCLUSIONS

Our study reveals that Srsf3-dependent generation of long 3' untranslated region is required for efficient mitochondrial translation, which promotes mature phagocytic macrophage formation, thereby playing a protective role in atherosclerosis.

Morphogenesis, starvation, and light responses in a mushroom-forming fungus revealed by long-read sequencing and extensive expression profiling

Mushroom-forming fungi (Agaricomycetes) are emerging as pivotal players in several fields of science and industry. Genomic data for Agaricomycetes are accumulating rapidly; however, this is not paralleled by improvements of gene annotations, which leave gene function notoriously poorly understood. We set out to improve our functional understanding of the model mushroom Coprinopsis cinerea by integrating a new, chromosome-level assembly, high-quality gene predictions, and functional information derived from broad gene-expression profiling data. The new annotation includes 5' and 3' untranslated regions (UTRs), polyadenylation sites (PASs), upstream open reading frames (uORFs), splicing isoforms, and microexons, as well as core gene sets corresponding to carbon starvation, light response, and hyphal differentiation. As a result, the genome of C. cinerea has now become the most comprehensively annotated genome among mushroom-forming fungi, which will contribute to multiple rapidly expanding fields, including research on their life history, light and stress responses, as well as multicellular development.

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.

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

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

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.

Testis- and ovary-expressed polo-like kinase transcripts and gene duplications affect male fertility when expressed in the Drosophila melanogaster germline

Polo-like kinases (Plks) are essential for spindle attachment to the kinetochore during prophase and the subsequent dissociation after anaphase in both mitosis and meiosis. There are structural differences in the spindle apparatus among mitosis, male meiosis, and female meiosis. It is therefore possible that alleles of Plk genes could improve kinetochore attachment or dissociation in spermatogenesis or oogenesis, but not both. These opposing effects could result in sexually antagonistic selection at Plk loci. In addition, Plk genes have been independently duplicated in many different evolutionary lineages within animals. This raises the possibility that Plk gene duplication may resolve sexual conflicts over mitotic and meiotic functions. We investigated this hypothesis by comparing the evolution, gene expression, and functional effects of the single Plk gene in Drosophila melanogaster (polo) and the duplicated Plks in D. pseudoobscura (Dpse-polo and Dpse-polo-dup1). Dpse-polo-dup1 is expressed primarily in testis, while other Drosophila Plk genes have broader expression profiles. We found that the protein-coding sequence of Dpse-polo-dup1 is evolving significantly faster than a canonical polo gene across all functional domains, yet the essential structure of the encoded protein has been retained. We present additional evidence that the faster evolution of Dpse-polo-dup1 is driven by the adaptive fixation of amino acid substitutions. We also found that over or ectopic expression of polo or Dpse-polo in the D. melanogaster male germline resulted in greater male infertility than expression of Dpse-polo-dup1. Last, expression of Dpse-polo or an ovary-derived transcript of polo in the male germline caused males to sire female-biased broods, suggesting that some Plk transcripts can affect the meiotic transmission of the sex chromosomes in the male germline. However, there was no sex bias in the progeny when Dpse-polo-dup1 was ectopically expressed, or a testis-derived transcript of polo was overexpressed in the D. melanogaster male germline. Our results therefore suggest that Dpse-polo-dup1 may have experienced positive selection to improve its regulation of the male meiotic spindle, resolving sexual conflict over meiotic Plk functions. Alternatively, Dpse-polo-dup1 may encode a hypomorphic Plk that has reduced deleterious effects when overexpressed in the male germline. Similarly, testis transcripts of D. melanogaster polo may be optimized for regulating the male meiotic spindle, and we provide evidence that the untranslated regions of the polo transcript may be involved in sex-specific germline functions.

SINE retrotransposons import polyadenylation signals to 3'UTRs in dog (Canis familiaris)

BACKGROUND

Messenger RNA 3' untranslated regions (3'UTRs) control many aspects of gene expression and determine where the transcript will terminate. The polyadenylation signal (PAS) AAUAAA (AATAAA in DNA) is a key regulator of transcript termination and this hexamer, or a similar sequence, is very frequently found within 30 bp of 3'UTR ends. Short interspersed element (SINE) retrotransposons are found throughout genomes in high copy numbers. When inserted into genes they can disrupt expression, alter splicing, or cause nuclear retention of mRNAs. The genomes of the domestic dog and other carnivores carry hundreds of thousands of Can-SINEs, a tRNA-related SINE with transcription termination potential. Because of this we asked whether Can-SINEs may terminate transcript in some dog genes.

RESULTS

Each of the dog's nine Can-SINE consensus sequences carry an average of three AATAAA PASs on their sense strands but zero on their antisense strands. Consistent with the idea that Can-SINEs can terminate transcripts, we find that sense-oriented Can-SINEs are approximately ten times more frequent at 3' ends of 3'UTRs compared to further upstream within 3'UTRs. Furthermore, the count of AATAAA PASs on head-to-tail SINE sequences differs significantly between sense and antisense-oriented retrotransposons in transcripts. Can-SINEs near 3'UTR ends are likely to carry an AATAAA motif on the mRNA sense strand while those further upstream are not. We identified loci where Can-SINE insertion has truncated or altered a 3'UTR of the dog genome (dog 3'UTR) compared to the human ortholog. Dog 3'UTRs have peaks of AATAAA PAS frequency at 28, 32, and 36 bp from the end. The periodicity is partly explained by TAAA(n) repeats within Can-SINE AT-rich tails. We annotated all repeat-masked Can-SINE copies in the Boxer reference genome and found that the young SINEC_Cf type has a mode of 15 bp length for target site duplications (TSDs). All dog Can-SINE types favor integration at TSDs beginning with A(4).

CONCLUSION

Dog Can-SINE retrotransposition has imported AATAAA PASs into gene transcripts and led to alteration of 3'UTRs. AATAAA sequences are selectively removed from Can-SINEs in introns and upstream 3'UTR regions but are retained at the far downstream end of 3'UTRs, which we infer reflects their role as termination sequences for these transcripts.

Post-transcriptional dynamics and RNA homeostasis in autophagy and cancer

Autophagy is an essential recycling and quality control pathway which preserves cellular and organismal homeostasis. As a catabolic process, autophagy degrades damaged and aged intracellular components in response to conditions of stress, including nutrient deprivation, oxidative and genotoxic stress. Autophagy is a highly adaptive and dynamic process which requires an intricately coordinated molecular control. Here we provide an overview of how autophagy is regulated post-transcriptionally, through RNA processing events, epitranscriptomic modifications and non-coding RNAs. We further discuss newly revealed RNA-binding properties of core autophagy machinery proteins and review recent indications of autophagy's ability to impact cellular RNA homeostasis. From a physiological perspective, we examine the biological implications of these emerging regulatory layers of autophagy, particularly in the context of nutrient deprivation and tumorigenesis.

Fateful Decisions of Where to Cut the Line: Pathology Associated with Aberrant 3' End Processing and Transcription Termination

Aberrant gene expression lies at the heart of many pathologies. This review will point out how 3' end processing, the final mRNA-maturation step in the transcription cycle, is surprisingly prone to regulated as well as stochastic variations with a wide range of consequences. Whereas smaller variations contribute to the plasticity of gene expression, larger alternations to 3' end processing and coupled transcription termination can lead to pathological consequences. These can be caused by the local mutation of one gene or affect larger numbers of genes systematically, if aspects of the mechanisms of 3' end processing and transcription termination are altered.

Single-cell landscape of alternative polyadenylation in human lymphoid hematopoiesis

Alternative polyadenylation (APA) is an essential post-transcriptional process that produces mature mRNA isoforms by regulating the usage of polyadenylation sites (PASs). APA is involved in lymphocyte activation; however, its role throughout the entire differentiation trajectory remains elusive. Here, we analyzed single-cell 3'-end transcriptome data from healthy subjects to construct a dynamic-APA landscape from hematopoietic stem and progenitor cells (HSPCs) to terminally differentiated lymphocytes. This analysis covered 19973 cells of 12 clusters from five lineages (B cells, CD4+ T cells, CD8+ T cells, natural killer cells, and plasmacytoid dendritic cells). A total of 2364 genes exhibited differential 3'-untranslated region (3'UTR) PAS usage, and 3021 genes displayed differential intronic cleavage during lymphoid differentiation. We observed a global trend of 3'UTR shortening during lymphoid differentiation. Nevertheless, specific events of both 3'UTR shortening and lengthening were also identified within each cluster. The APA patterns delineated three differentiation stages: HSPCs, precursor cells, and mature cells. Moreover, we demonstrated that the conversion of naïve T cells to memory T cells was accompanied by dynamic APA in transcription factor-encoding genes (TCF7 and NFATC2IP), immune function-related genes (BCL2, CD5, CD28, GOLT1B, and TMEM59), and protein ubiquitination-related genes (UBE2G1, YPEL5, and SUMO3). These findings expand our understanding of the underlying molecular mechanisms of APA and facilitate studies on the regulatory role of APA in lymphoid hematopoiesis.

Alternative mRNA polyadenylation regulates macrophage hyperactivation via the autophagy pathway

Macrophage hyperactivation is a hallmark of inflammatory diseases, yet the role of alternative polyadenylation (APA) of mRNAs in regulating innate immunity remains unclear. In this study, we focused on 3'UTR-APA and demonstrated that Nudt21, a crucial RNA-binding component of the 3'UTR-APA machinery, is significantly upregulated in various inflammatory conditions. By utilizing myeloid-specific Nudt21-deficient mice, we revealed a protective effect of Nudt21 depletion against colitis and severe hyperinflammation, primarily through diminished production of proinflammatory cytokines. Notably, Nudt21 regulates the mRNA stability of key autophagy-related genes, Map1lc3b and Ulk2, by mediating selective 3'UTR polyadenylation in activated macrophages. As a result, Nudt21-deficient macrophages display increased autophagic activity, which leads to reduced cytokine secretion. Our findings highlight an unexplored role of Nudt21-mediated 3'UTR-APA in modulating macrophage autophagy and offer new insights into the modulation of inflammation and disease progression.

Understanding isoform expression by pairing long-read sequencing with single-cell and spatial transcriptomics

RNA isoform diversity, produced via alternative splicing, and alternative usage of transcription start and poly(A) sites, results in varied transcripts being derived from the same gene. Distinct isoforms can play important biological roles, including by changing the sequences or expression levels of protein products. The first single-cell approaches to RNA sequencing-and later, spatial approaches-which are now widely used for the identification of differentially expressed genes, rely on short reads and offer the ability to transcriptomically compare different cell types but are limited in their ability to measure differential isoform expression. More recently, long-read sequencing methods have been combined with single-cell and spatial technologies in order to characterize isoform expression. In this review, we provide an overview of the emergence of single-cell and spatial long-read sequencing and discuss the challenges associated with the implementation of these technologies and interpretation of these data. We discuss the opportunities they offer for understanding the relationships between the distinct variable elements of transcript molecules and highlight some of the ways in which they have been used to characterize isoforms' roles in development and pathology. Single-nucleus long-read sequencing, a special case of the single-cell approach, is also discussed. We attempt to cover both the limitations of these technologies and their significant potential for expanding our still-limited understanding of the biological roles of RNA isoforms.

YTHDF2 upregulation and subcellular localization dictate CD8 T cell polyfunctionality in anti-tumor immunity

RNA methylation is an important regulatory process to determine immune cell function but how it affects the anti-tumor activity of CD8 T cells is not fully understood. Here we show that the N6-methyladenosine (m6A) RNA reader YTHDF2 is highly expressed in early effector or effector-like CD8 T cells. We find that YTHDF2 facilitates nascent RNA synthesis, and m6A recognition is fundamental for this distinctively nuclear function of the protein, which also reinforces its autoregulation at the RNA level. Loss of YTHDF2 in T cells exacerbates tumor progression and confers unresponsiveness to PD-1 blockade in mice and in humans. In addition to initiating RNA decay that is necessary for mitochondrial fitness, YTHDF2 orchestrates chromatin changes that promote T cell polyfunctionality. YTHDF2 interacts with IKZF1/3, which is important for sustained transcription of their target genes. Accordingly, immunotherapy-induced efficacy could be largely restored in YTHDF2-deficient T cells through combinational use of IKZF1/3 inhibitor lenalidomide in a mouse model. Thus, YTHDF2 coordinates epi-transcriptional and transcriptional networks to potentiate T cell immunity, which could inform therapeutic intervention.

Aberrant pre-mRNA processing in cancer

Dysregulation of the flow of information from genomic DNA to RNA to protein occurs within all cancer types. In this review, we described the current state of understanding of how RNA processing is dysregulated in cancer with a focus on mutations in the RNA splicing factor machinery that are highly prevalent in hematologic malignancies. We discuss the downstream effects of these mutations highlighting both individual genes as well as common pathways that they perturb. We highlight examples of how alterations in RNA processing have been harnessed for therapeutic intent as well as to promote the selective toxicity of cancer cells.

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

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

APC mutations dysregulate alternative polyadenylation in cancer

BACKGROUND

Alternative polyadenylation (APA) affects most human genes and is recurrently dysregulated in all studied cancers. However, the mechanistic origins of this dysregulation are incompletely understood.

RESULTS

We describe an unbiased analysis of molecular regulators of poly(A) site selection across The Cancer Genome Atlas and identify that colorectal adenocarcinoma is an outlier relative to all other cancer subtypes. This distinction arises from the frequent presence of loss-of-function APC mutations in colorectal adenocarcinoma, which are strongly associated with long 3' UTR expression relative to tumors lacking APC mutations. APC knockout similarly dysregulates APA in human colon organoids. By mining previously published APC eCLIP data, we show that APC preferentially binds G- and C-rich motifs just upstream of proximal poly(A) sites. Lastly, we find that reduced APC expression is associated with APA dysregulation in tumor types lacking recurrent APC mutations.

CONCLUSIONS

As APC has been previously identified as an RNA-binding protein that preferentially binds 3' UTRs during mouse neurogenesis, our results suggest that APC promotes proximal poly(A) site use and that APC loss and altered expression contribute to pervasive APA dysregulation in cancers.

Involvement of oncomiRs miR-23, miR-24, and miR-27 in the regulation of alternative polyadenylation in glioblastoma via CFIm25 cleavage factor

Glioblastoma multiforme (GBM) is a highly aggressive brain tumor with a poor prognosis. The cleavage factor Im 25 (CFIm25), a crucial component of the CFIm complex, plays a key role in regulating the length of the mRNA 3'-UTR and has been implicated in various cancers, including GBM. This study sought to investigate the regulatory influence of specific microRNAs (miRNAs) on CFIm25 expression in GBM, a highly aggressive brain tumor. Bioinformatics analysis identified miRNA candidates targeting CFIm25 mRNA, and gene expression profiles from the NCBI database (GSE90603) were used for further analysis. Expression levels of CFIm25 and selected miRNAs were assessed using qRT-PCR in GBM clinical samples (n = 20) and non-malignant brain tissues (n = 5). Additionally, the MTT assay was performed to examine the effect of miRNA overexpression on U251 cell viability. Lentivectors expressing the identified miRNAs were employed to experimentally validate their regulatory role on CFIm25 in U251 cell lines, and Western blot analysis was conducted to determine CFIm25 protein levels. We observed significantly increased levels of miR-23, miR-24, and miR-27 expression, associated with a marked reduction in CFIm25 expression in GBM samples compared to non-malignant brain tissues. In particular, overexpression of miR-23, miR-24, and miR-27 in U251 cells resulted in CFIm25 downregulation at both the mRNA and protein levels, while their inhibition increased CFIm25 and reduced cell proliferation. These observations strongly implicate miR-23, miR-24, and miR-27 in regulating CFIm25 expression in GBM, emphasizing their potential as promising therapeutic targets for enhancing treatment responses in glioblastoma.

The landscape and clinical relevance of intronic polyadenylation in human cancers

Intronic polyadenylation (IPA) is an RNA 3' end processing event which has been reported to play important roles in cancer development. However, the comprehensive landscape of IPA events across various cancer types is lacking. Here, we apply IPAFinder to identify and quantify IPA events in 10,383 samples covering all 33 cancer types from The Cancer Genome Atlas (TCGA) project. We identify a total of 21,835 IPA events, almost half of which are ubiquitously expressed. We identify 2761 unique dynamically changed IPA events across cancer types. Furthermore, we observe 8855 non-redundant clinically relevant IPA events, which could potentially be used as prognostic indicators. Our analysis also reveals that dynamic IPA usage within cancer signaling pathways may affect drug response. Finally, we develop a user-friendly data portal, IPACancer Atlas (http://www.tingni-lab.com/Pancan_IPA/), to search and explore IPAs in cancer.

Leveraging multi-omics data to infer regulators of mRNA 3' end processing in glioblastoma

Alterations in mRNA 3' end processing and polyadenylation are widely implicated in the biology of many cancer types, including glioblastoma (GBM), one the most aggressive tumor types. Although several RNA-binding proteins (RBPs) responsible for alternative polyadenylation (APA) were identified from functional studies in cell lines, their contribution to the APA landscape in tumors in vivo was not thoroughly addressed. In this study we analyzed a large RNA-seq data set of glioblastoma (GBM) samples from The Cancer Genome Atlas (TCGA) to identify APA patterns differentiating the main molecular subtypes of GBM. We superimposed these to RBP footprinting data and to APA events occurring upon depletion of individual RBPs from a large panel tested by the ENCODE Consortium. Our analysis revealed 22 highly concordant and statistically significant RBP-APA associations, whereby changes in RBP expression were accompanied by APA in both TCGA and ENCODE datasets. Among these, we found a previously unknown PTBP1-regulated APA event in the PRRC2B gene and an HNRNPU-regulated event in the SC5D gene. Both of these were further supported by RNA-sequencing data of paired tumor center-periphery GBM samples obtained at the University Hospital of Basel. In addition, we validated the regulation of APA in PRRC2B by PTBP1 in siRNA-knockdown and overexpression experiments followed by RNA-sequencing in two glioblastoma cell lines. The transcriptome analysis workflow that we present here enables the identification of concordant RBP-APA associations in cancers.

Upstream alternative polyadenylation in SCN5A produces a short transcript isoform encoding a mitochondria-localized NaV1.5 N-terminal fragment that influences cardiomyocyte respiration

SCN5A encodes the cardiac voltage-gated Na+ channel, NaV1.5, that initiates action potentials. SCN5A gene variants cause arrhythmias and increased heart failure risk. Mechanisms controlling NaV1.5 expression and activity are not fully understood. We recently found a well-conserved alternative polyadenylation (APA) signal downstream of the first SCN5A coding exon. This yields a SCN5A-short transcript isoform expressed in several species (e.g. human, pig, and cat), though rodents lack this upstream APA. Reanalysis of transcriptome-wide cardiac APA-seq and mRNA-seq data shows reductions in both upstream APA usage and short/full-length SCN5A mRNA ratios in failing hearts. Knock-in of the human SCN5A APA sequence into mice is sufficient to enable expression of SCN5A -short transcript, while significantly decreasing expression of full-length SCN5A mRNA. Notably, SCN5A -short transcript encodes a novel protein (NaV1.5-NT), composed of an N-terminus identical to NaV1.5 and a unique C-terminus derived from intronic sequence. AAV9 constructs were able to achieve stable NaV1.5-NT expression in mouse hearts, and western blot of human heart tissues showed bands co-migrating with NaV1.5-NT transgene-derived bands. NaV1.5-NT is predicted to contain a mitochondrial targeting sequence and localizes to mitochondria in cultured cardiomyocytes and in mouse hearts. NaV1.5-NT expression in cardiomyocytes led to elevations in basal oxygen consumption rate, ATP production, and mitochondrial ROS, while depleting NADH supply. Native PAGE analyses of mitochondria lysates revealed that NaV1.5-NT expression resulted in increased levels of disassembled complex V subunits and accumulation of complex I-containing supercomplexes. Overall, we discovered that APA-mediated regulation of SCN5A produces a short transcript encoding NaV1.5-NT. Our data support that NaV1.5-NT plays a multifaceted role in influencing mitochondrial physiology: 1) by increasing basal respiration likely through promoting complex V conformations that enhance proton leak, and 2) by increasing overall respiratory efficiency and NADH consumption by enhancing formation and/or stability of complex I-containing respiratory supercomplexes, though the specific molecular mechanisms underlying each of these remain unresolved.

Multiplexed single-cell characterization of alternative polyadenylation regulators

Most mammalian genes have multiple polyA sites, representing a substantial source of transcript diversity regulated by the cleavage and polyadenylation (CPA) machinery. To better understand how these proteins govern polyA site choice, we introduce CPA-Perturb-seq, a multiplexed perturbation screen dataset of 42 CPA regulators with a 3' scRNA-seq readout that enables transcriptome-wide inference of polyA site usage. We develop a framework to detect perturbation-dependent changes in polyadenylation and characterize modules of co-regulated polyA sites. We find groups of intronic polyA sites regulated by distinct components of the nuclear RNA life cycle, including elongation, splicing, termination, and surveillance. We train and validate a deep neural network (APARENT-Perturb) for tandem polyA site usage, delineating a cis-regulatory code that predicts perturbation response and reveals interactions between regulatory complexes. Our work highlights the potential for multiplexed single-cell perturbation screens to further our understanding of post-transcriptional regulation.

Dynamic Alternative Polyadenylation during Litopenaeus Vannamei Metamorphosis Development

As an important mechanism in the post-transcriptional regulation of eukaryotic gene expression, alternative polyadenylation (APA) plays a key role in biological processes such as cell proliferation and differentiation. However, the role and dynamic pattern of APA during Litopenaeus vannamei metamorphosis are poorly understood. Here, RNA-seq data covering from the embryo to the maturation (16 time points) of L. vannamei were utilized. We identified 247 differentially expressed APA events between early and adult stages, and through fuzzy mean clustering analysis, we discovered five dynamic APA patterns. Among them, the gradual elongation of the 3'UTR is the major APA pattern that changes over time, and its genes are enriched in the pathways of protein and energy metabolism. Finally, we constructed mRNA-miRNA and PPI networks and detected several central miRNAs that may regulate L. vannamei development. Our results revealed the complex APA mechanisms in L. vannamei metamorphosis, shedding new light on post-transcriptional regulation of crustacean metamorphosis.

Integration mapping of cardiac fibroblast single-cell transcriptomes elucidates cellular principles of fibrosis in diverse pathologies

Single-cell technology has allowed researchers to probe tissue complexity and dynamics at unprecedented depth in health and disease. However, the generation of high-dimensionality single-cell atlases and virtual three-dimensional tissues requires integrated reference maps that harmonize disparate experimental designs, analytical pipelines, and taxonomies. Here, we present a comprehensive single-cell transcriptome integration map of cardiac fibrosis, which underpins pathophysiology in most cardiovascular diseases. Our findings reveal similarity between cardiac fibroblast (CF) identities and dynamics in ischemic versus pressure overload models of cardiomyopathy. We also describe timelines for commitment of activated CFs to proliferation and myofibrogenesis, profibrotic and antifibrotic polarization of myofibroblasts and matrifibrocytes, and CF conservation across mouse and human healthy and diseased hearts. These insights have the potential to inform knowledge-based therapies.

MicroRNA-mediated regulation of nonsense-mediated mRNA decay factors: Insights into microRNA prediction tools and profiling techniques

Nonsense-mediated mRNA decay (NMD) stands out as a prominent RNA surveillance mechanism within eukaryotes, meticulously overseeing both RNA abundance and integrity by eliminating aberrant transcripts. These defective transcripts are discerned through the concerted efforts of translating ribosomes, eukaryotic release factors (eRFs), and trans-acting NMD factors, with Up-Frameshift 3 (UPF3) serving as a noteworthy component. Remarkably, in humans, UPF3 exists in two paralogous forms, UPF3A (UPF3) and UPF3B (UPF3X). Beyond its role in quality control, UPF3 wields significant influence over critical cellular processes, including neural development, synaptic plasticity, and axon guidance. However, the precise regulatory mechanisms governing UPF3 remain elusive. MicroRNAs (miRNAs) emerge as pivotal post-transcriptional gene regulators, exerting substantial impact on diverse pathological and physiological pathways. This comprehensive review encapsulates our current understanding of the intricate regulatory nexus between NMD and miRNAs, with particular emphasis on the essential role played by UPF3B in neurodevelopment. Additionally, we bring out the significance of the 3'-untranslated region (3'-UTR) as the molecular bridge connecting NMD and miRNA-mediated gene regulation. Furthermore, we provide an in-depth exploration of diverse computational tools tailored for the prediction of potential miRNA targets. To complement these computational approaches, we delineate experimental techniques designed to validate predicted miRNA-mRNA interactions, empowering readers with the knowledge necessary to select the most appropriate methodology for their specific research objectives.

Quantifying 3'UTR length from scRNA-seq data reveals changes independent of gene expression

Although more than half of all genes generate transcripts that differ in 3'UTR length, current analysis pipelines only quantify the amount but not the length of mRNA transcripts. 3'UTR length is determined by 3' end cleavage sites (CS). We map CS in more than 200 primary human and mouse cell types and increase CS annotations relative to the GENCODE database by 40%. Approximately half of all CS are used in few cell types, revealing that most genes only have one or two major 3' ends. We incorporate the CS annotations into a computational pipeline, called scUTRquant, for rapid, accurate, and simultaneous quantification of gene and 3'UTR isoform expression from single-cell RNA sequencing (scRNA-seq) data. When applying scUTRquant to data from 474 cell types and 2134 perturbations, we discover extensive 3'UTR length changes across cell types that are as widespread and coordinately regulated as gene expression changes but affect mostly different genes. Our data indicate that mRNA abundance and mRNA length are two largely independent axes of gene regulation that together determine the amount and spatial organization of protein synthesis.

Alternative 3'UTR expression induced by T cell activation is regulated in a temporal and signal dependent manner

The length of 3' untranslated regions (3'UTR) is highly regulated during many transitions in cell state, including T cell activation, through the process of alternative polyadenylation (APA). However, the regulatory mechanisms and functional consequences of APA remain largely unexplored. Here we present a detailed analysis of the temporal and condition-specific regulation of APA following activation of primary human CD4+ T cells. We find that global APA changes are regulated temporally and CD28 costimulatory signals enhance a subset of these changes. Most APA changes upon T cell activation involve 3'UTR shortening, although a set of genes enriched for function in the mTOR pathway exhibit 3'UTR lengthening. While upregulation of the core polyadenylation machinery likely induces 3'UTR shortening following prolonged T cell stimulation; a significant program of APA changes occur prior to cellular proliferation or upregulation of the APA machinery. Motif analysis suggests that at least a subset of these early changes in APA are driven by upregulation of RBM3, an RNA-binding protein which competes with the APA machinery for binding. Together this work expands our understanding of the impact and mechanisms of APA in response to T cell activation and suggests new mechanisms by which APA may be regulated.

Nuclear export is a limiting factor in eukaryotic mRNA metabolism

The eukaryotic mRNA life cycle includes transcription, nuclear mRNA export and degradation. To quantify all these processes simultaneously, we perform thiol-linked alkylation after metabolic labeling of RNA with 4-thiouridine (4sU), followed by sequencing of RNA (SLAM-seq) in the nuclear and cytosolic compartments of human cancer cells. We develop a model that reliably quantifies mRNA-specific synthesis, nuclear export, and nuclear and cytosolic degradation rates on a genome-wide scale. We find that nuclear degradation of polyadenylated mRNA is negligible and nuclear mRNA export is slow, while cytosolic mRNA degradation is comparatively fast. Consequently, an mRNA molecule generally spends most of its life in the nucleus. We also observe large differences in the nuclear export rates of different 3'UTR transcript isoforms. Furthermore, we identify genes whose expression is abruptly induced upon metabolic labeling. These transcripts are exported substantially faster than average mRNAs, suggesting the existence of alternative export pathways. Our results highlight nuclear mRNA export as a limiting factor in mRNA metabolism and gene regulation.

Translation-dependent and -independent mRNA decay occur through mutually exclusive pathways defined by ribosome density during T cell activation

mRNA translation and decay are tightly interconnected processes both in the context of mRNA quality-control pathways and for the degradation of functional mRNAs. Cotranslational mRNA degradation through codon usage, ribosome collisions, and the recruitment of specific proteins to ribosomes is an important determinant of mRNA turnover. However, the extent to which translation-dependent mRNA decay (TDD) and translation-independent mRNA decay (TID) pathways participate in the degradation of mRNAs has not been studied yet. Here we describe a comprehensive analysis of basal and signal-induced TDD and TID in mouse primary CD4+ T cells. Our results indicate that most cellular transcripts are decayed to some extent in a translation-dependent manner. Our analysis further identifies the length of untranslated regions, the density of ribosomes, and GC3 content as important determinants of TDD magnitude. Consistently, all transcripts that undergo changes in ribosome density within their coding sequence upon T cell activation display a corresponding change in their TDD level. Moreover, we reveal a dynamic modulation in the relationship between GC3 content and TDD upon T cell activation, with a reversal in the impact of GC3- and AU3-rich codons. Altogether, our data show a strong and dynamic interconnection between mRNA translation and decay in mammalian primary cells.

InPACT: a computational method for accurate characterization of intronic polyadenylation from RNA sequencing data

Alternative polyadenylation can occur in introns, termed intronic polyadenylation (IPA), has been implicated in diverse biological processes and diseases, as it can produce noncoding transcripts or transcripts with truncated coding regions. However, a reliable method is required to accurately characterize IPA. Here, we propose a computational method called InPACT, which allows for the precise characterization of IPA from conventional RNA-seq data. InPACT successfully identifies numerous previously unannotated IPA transcripts in human cells, many of which are translated, as evidenced by ribosome profiling data. We have demonstrated that InPACT outperforms other methods in terms of IPA identification and quantification. Moreover, InPACT applied to monocyte activation reveals temporally coordinated IPA events. Further application on single-cell RNA-seq data of human fetal bone marrow reveals the expression of several IPA isoforms in a context-specific manner. Therefore, InPACT represents a powerful tool for the accurate characterization of IPA from RNA-seq data.

T helper cells exhibit a dynamic and reversible 3'-UTR landscape

3' untranslated regions (3' UTRs) are critical elements of messenger RNAs, as they contain binding sites for RNA-binding proteins (RBPs) and microRNAs that affect various aspects of the RNA life cycle including transcript stability and cellular localization. In response to T cell receptor activation, T cells undergo massive expansion during the effector phase of the immune response and dynamically modify their 3' UTRs. Whether this serves to directly regulate the abundance of specific mRNAs or is a secondary effect of proliferation remains unclear. To study 3'-UTR dynamics in T helper cells, we investigated division-dependent alternative polyadenylation (APA). In addition, we generated 3' end UTR sequencing data from naive, activated, memory, and regulatory CD4+ T cells. 3'-UTR length changes were estimated using a nonnegative matrix factorization approach and were compared with those inferred from long-read PacBio sequencing. We found that APA events were transient and reverted after effector phase expansion. Using an orthogonal bulk RNA-seq data set, we did not find evidence of APA association with differential gene expression or transcript usage, indicating that APA has only a marginal effect on transcript abundance. 3'-UTR sequence analysis revealed conserved binding sites for T cell-relevant microRNAs and RBPs in the alternative 3' UTRs. These results indicate that poly(A) site usage could play an important role in the control of cell fate decisions and homeostasis.

Utilizing non-coding RNA-mediated regulation of ATP binding cassette (ABC) transporters to overcome multidrug resistance to cancer chemotherapy

Multidrug resistance (MDR) is one of the primary factors that produces treatment failure in patients receiving cancer chemotherapy. MDR is a complex multifactorial phenomenon, characterized by a decrease or abrogation of the efficacy of a wide spectrum of anticancer drugs that are structurally and mechanistically distinct. The overexpression of the ATP-binding cassette (ABC) transporters, notably ABCG2 and ABCB1, are one of the primary mediators of MDR in cancer cells, which promotes the efflux of certain chemotherapeutic drugs from cancer cells, thereby decreasing or abolishing their therapeutic efficacy. A number of studies have suggested that non-coding RNAs (ncRNAs), particularly microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs), play a pivotal role in mediating the upregulation of ABC transporters in certain MDR cancer cells. This review will provide updated information about the induction of ABC transporters due to the aberrant regulation of ncRNAs in cancer cells. We will also discuss the measurement and biological profile of circulating ncRNAs in various body fluids as potential biomarkers for predicting the response of cancer patients to chemotherapy. Sequence variations, such as alternative polyadenylation of mRNA and single nucleotide polymorphism (SNPs) at miRNA target sites, which may indicate the interaction of miRNA-mediated gene regulation with genetic variations to modulate the MDR phenotype, will be reviewed. Finally, we will highlight novel strategies that could be used to modulate ncRNAs and circumvent ABC transporter-mediated MDR.

Nuclear localization of Argonaute 2 is affected by cell density and may relieve repression by microRNAs

Argonaute protein is associated with post-transcriptional control of cytoplasmic gene expression through miRNA-induced silencing complexes (miRISC). Specific cellular and environmental conditions can trigger AGO protein to accumulate in the nucleus. Localization of AGO is central to understanding miRNA action, yet the consequences of AGO being in the nucleus are undefined. We show nuclear enrichment of AGO2 in HCT116 cells grown in two-dimensional culture to high density, HCT116 cells grown in three-dimensional tumor spheroid culture, and human colon tumors. The shift in localization of AGO2 from cytoplasm to nucleus de-represses cytoplasmic AGO2-eCLIP targets that were candidates for canonical regulation by miRISC. Constitutive nuclear localization of AGO2 using an engineered nuclear localization signal increases cell migration. Critical RNAi factors also affect the localization of AGO2. Knocking out an enzyme essential for miRNA biogenesis, DROSHA, depletes mature miRNAs and restricts AGO2 localization to the cytoplasm, while knocking out the miRISC scaffolding protein, TNRC6, results in nuclear localization of AGO2. These data suggest that AGO2 localization and miRNA activity can be regulated depending on environmental conditions, expression of mature miRNAs, and expression of miRISC cofactors. Localization and expression of core miRISC protein machinery should be considered when investigating the roles of miRNAs.

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.

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

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

Alternative cleavage and polyadenylation of the Ccnb1 mRNA defines accumulation of cyclin protein during the meiotic cell cycle

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

Alternative polyadenylation and dynamic 3' UTR length is associated with polysome recruitment throughout the cardiomyogenic differentiation of hESCs

Alternative polyadenylation (APA) increases transcript diversity through the generation of isoforms with varying 3' untranslated region (3' UTR) lengths. As the 3' UTR harbors regulatory element target sites, such as miRNAs or RNA-binding proteins, changes in this region can impact post-transcriptional regulation and translation. Moreover, the APA landscape can change based on the cell type, cell state, or condition. Given that APA events can impact protein expression, investigating translational control is crucial for comprehending the overall cellular regulation process. Revisiting data from polysome profiling followed by RNA sequencing, we investigated the cardiomyogenic differentiation of pluripotent stem cells by identifying the transcripts that show dynamic 3' UTR lengthening or shortening, which are being actively recruited to ribosome complexes. Our findings indicate that dynamic 3' UTR lengthening is not exclusively associated with differential expression during cardiomyogenesis but rather with recruitment to polysomes. We confirm that the differentiated state of cardiomyocytes shows a preference for shorter 3' UTR in comparison to the pluripotent stage although preferences vary during the days of the differentiation process. The most distinct regulatory changes are seen in day 4 of differentiation, which is the mesoderm commitment time point of cardiomyogenesis. After identifying the miRNAs that would target specifically the alternative 3' UTR region of the isoforms, we constructed a gene regulatory network for the cardiomyogenesis process, in which genes related to the cell cycle were identified. Altogether, our work sheds light on the regulation and dynamic 3' UTR changes of polysome-recruited transcripts that take place during the cardiomyogenic differentiation of pluripotent stem cells.

Biocompatibility and proteomic profiling of DMSA-coated iron nanocubes in a human glioblastoma cell line

Background: Superparamagnetic iron core iron oxide shell nanocubes have previously shown superior performance in magnetic resonance imaging T2 contrast enhancement compared with spherical nanoparticles. Methods: Iron core iron oxide shell nanocubes were synthesized, stabilized with dimercaptosuccinic acid (DMSA-NC) and physicochemically characterized. MRI contrast enhancement and biocompatibility were assessed in vitro. Results: DMSA-NC showed a transverse relaxivity of 122.59 mM-1·s-1 Fe. Treatment with DMSA-NC did not induce cytotoxicity or oxidative stress in U-251 cells, and electron microscopy demonstrated DMSA-NC localization within endosomes and lysosomes in cells following internalization. Global proteomics revealed dysregulation of iron storage, transport, transcription and mRNA processing proteins. Conclusion: DMSA-NC is a promising T2 MRI contrast agent which, in this preliminary investigation, demonstrates favorable biocompatibility with an astrocyte cell model.

Multiplexed screening reveals how cancer-specific alternative polyadenylation shapes tumor growth in vivo

Alternative polyadenylation (APA) is strikingly dysregulated in many cancers. Although global APA dysregulation is frequently associated with poor prognosis, the importance of most individual APA events is controversial simply because few have been functionally studied. Here, we address this gap by developing a CRISPR-Cas9-based screen to manipulate endogenous polyadenylation and systematically quantify how APA events contribute to tumor growth in vivo. Our screen reveals individual APA events that control mouse melanoma growth in an immunocompetent host, with concordant associations in clinical human cancer. For example, forced Atg7 3' UTR lengthening in mouse melanoma suppresses ATG7 protein levels, slows tumor growth, and improves host survival; similarly, in clinical human melanoma, a long ATG7 3' UTR is associated with significantly prolonged patient survival. Overall, our study provides an easily adaptable means to functionally dissect APA in physiological systems and directly quantifies the contributions of recurrent APA events to tumorigenic phenotypes.

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

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

mRNA initiation and termination are spatially coordinated

The expression of a precise mRNA transcriptome is crucial for establishing cell identity and function, with dozens of alternative isoforms produced for a single gene sequence. The regulation of mRNA isoform usage occurs by the coordination of co-transcriptional mRNA processing mechanisms across a gene. Decisions involved in mRNA initiation and termination underlie the largest extent of mRNA isoform diversity, but little is known about any relationships between decisions at both ends of mRNA molecules. Here, we systematically profile the joint usage of mRNA transcription start sites (TSSs) and polyadenylation sites (PASs) across tissues and species. Using both short and long read RNA-seq data, we observe that mRNAs preferentially using upstream TSSs also tend to use upstream PASs, and congruently, the usage of downstream sites is similarly paired. This observation suggests that mRNA 5' end choice may directly influence mRNA 3' ends. Our results suggest a novel "Positional Initiation-Termination Axis" (PITA), in which the usage of alternative terminal sites are coupled based on the order in which they appear in the genome. PITA isoforms are more likely to encode alternative protein domains and use conserved sites. PITA is strongly associated with the length of genomic features, such that PITA is enriched in longer genes with more area devoted to regions that regulate alternative 5' or 3' ends. Strikingly, we found that PITA genes are more likely than non-PITA genes to have multiple, overlapping chromatin structural domains related to pairing of ordinally coupled start and end sites. In turn, PITA coupling is also associated with fast RNA Polymerase II (RNAPII) trafficking across these long gene regions. Our findings indicate that a combination of spatial and kinetic mechanisms couple transcription initiation and mRNA 3' end decisions based on ordinal position to define the expression mRNA isoforms.

Alternative Polyadenylation of Malic Enzyme 1 Is Essential for Accelerated Adipogenesis

Understanding the mechanism of adipogenesis is an important basis for improving meat quality traits of livestock. Alternative polyadenylation (APA) is a vital mechanism to regulate the expression of eukaryotic genes. However, how the individual APA functions in adipogenesis remains elusive. This study was intended to investigate the effect of malic enzyme 1 (ME1) APA on adipogenesis. Here, intracellular lipid droplets were stained using Oil red O. 3' RACE was used to verify APA events of the ME1 gene. Interactions between ME1 3' untranslated region (3' UTR)-APA isoforms and miRNAs, as well as differential expression of isoforms, were examined using dual-luciferase reporter and molecular experiments. The mechanism of ME1 APA on adipogenesis was explored by gain and loss of function assays. In this study, two ME1 isoforms with different 3' UTR lengths were detected during adipogenesis. Moreover, the ME1 isoform with a short 3' UTR was significantly upregulated compared with the one with a long 3' UTR. Mechanistically, only the long ME1 isoform was targeted by miR-153-3p to attenuate adipogenesis, while the short one escaped the regulation of miR-153-3p to accelerate adipogenesis. Our results reveal a novel mechanism of ME1 APA in regulating adipogenesis.

Immune-response 3'UTR alternative polyadenylation quantitative trait loci contribute to variation in human complex traits and diseases

Genome-wide association studies (GWASs) have identified thousands of non-coding variants that are associated with human complex traits and diseases. The analysis of such GWAS variants in different contexts and physiological states is essential for deciphering the regulatory mechanisms underlying human disease. Alternative polyadenylation (APA) is a key post-transcriptional modification for most human genes that substantially impacts upon cell behavior. Here, we mapped 9,493 3'-untranslated region APA quantitative trait loci in 18 human immune baseline cell types and 8 stimulation conditions (immune 3'aQTLs). Through the comparison between baseline and stimulation data, we observed the high responsiveness of 3'aQTLs to immune stimulation (response 3'aQTLs). Co-localization and mendelian randomization analyses of immune 3'aQTLs identified 678 genes where 3'aQTL are associated with variation in complex traits, 27.3% of which were derived from response 3'aQTLs. Overall, these analyses reveal the role of immune 3'aQTLs in the determination of complex traits, providing new insights into the regulatory mechanisms underlying disease etiologies.

GCLiPP: global crosslinking and protein purification method for constructing high-resolution occupancy maps for RNA binding proteins

GCLiPP is a global RNA interactome capture method that detects RNA-binding protein (RBP) occupancy transcriptome-wide. GCLiPP maps RBP-occupied sites at a higher resolution than phase separation-based techniques. GCLiPP sequence tags correspond with known RBP binding sites and are enriched for sites detected by RBP-specific crosslinking immunoprecipitation (CLIP) for abundant cytosolic RBPs. Comparison of human Jurkat T cells and mouse primary T cells uncovers shared peaks of GCLiPP signal across homologous regions of human and mouse 3' UTRs, including a conserved mRNA-destabilizing cis-regulatory element. GCLiPP signal overlapping with immune-related SNPs uncovers stabilizing cis-regulatory regions in CD5, STAT6, and IKZF1.

Regulation of alternative splicing and polyadenylation in neurons

Cell-type-specific gene expression is a fundamental feature of multicellular organisms and is achieved by combinations of regulatory strategies. Although cell-restricted transcription is perhaps the most widely studied mechanism, co-transcriptional and post-transcriptional processes are also central to the spatiotemporal control of gene functions. One general category of expression control involves the generation of multiple transcript isoforms from an individual gene, whose balance and cell specificity are frequently tightly regulated via diverse strategies. The nervous system makes particularly extensive use of cell-specific isoforms, specializing the neural function of genes that are expressed more broadly. Here, we review regulatory strategies and RNA-binding proteins that direct neural-specific isoform processing. These include various classes of alternative splicing and alternative polyadenylation events, both of which broadly diversify the neural transcriptome. Importantly, global alterations of splicing and alternative polyadenylation are characteristic of many neural pathologies, and recent genetic studies demonstrate how misregulation of individual neural isoforms can directly cause mutant phenotypes.

Therapeutic targeting of CPSF3-dependent transcriptional termination in ovarian cancer

Transcriptional dysregulation is a recurring pathogenic hallmark and an emerging therapeutic vulnerability in ovarian cancer. Here, we demonstrated that ovarian cancer exhibited a unique dependency on the regulatory machinery of transcriptional termination, particularly, cleavage and polyadenylation specificity factor (CPSF) complex. Genetic abrogation of multiple CPSF subunits substantially hampered neoplastic cell viability, and we presented evidence that their indispensable roles converged on the endonuclease CPSF3. Mechanistically, CPSF perturbation resulted in lengthened 3'-untranslated regions, diminished intronic polyadenylation and widespread transcriptional readthrough, and consequently suppressed oncogenic pathways. Furthermore, we reported the development of specific CPSF3 inhibitors building upon the benzoxaborole scaffold, which exerted potent antitumor activity. Notably, CPSF3 blockade effectively exacerbated genomic instability by down-regulating DNA damage repair genes and thus acted in synergy with poly(adenosine 5'-diphosphate-ribose) polymerase inhibition. These findings establish CPSF3-dependent transcriptional termination as an exploitable driving mechanism of ovarian cancer and provide a promising class of boron-containing compounds for targeting transcription-addicted human malignancies.

Effects of Combinations of Untranslated-Region Sequences on Translation of mRNA

mRNA-based therapeutics have been found to be a promising treatment strategy in immunotherapy, gene therapy, and cancer treatments. Effectiveness of mRNA therapeutics depends on the level and duration of a desired protein's expression, which is determined by various cis- and trans-regulatory elements of the mRNA. Sequences of 5' and 3' untranslated regions (UTRs) are responsible for translational efficiency and stability of mRNA. An optimal combination of the regulatory sequences allows researchers to significantly increase the target protein's expression. Using both literature data and previously obtained experimental data, we chose six sequences of 5'UTRs (adenoviral tripartite leader [TPL], HBB, rabbit β-globin [Rabb], H4C2, Moderna, and Neo2) and five sequences of 3'UTRs (mtRNR-EMCV, mtRNR-AES, mtRNR-mtRNR, BioNTech, and Moderna). By combining them, we constructed 30 in vitro transcribed RNAs encoding firefly luciferase with various combinations of 5'- and 3'UTRs, and the resultant bioluminescence was assessed in the DC2.4 cell line at 4, 8, 24, and 72 h after transfection. The cellular data enabled us to identify the best seven combinations of 5'- and 3'UTRs, whose translational efficiency was then assessed in BALB/c mice. Two combinations of 5'- and 3'UTRs (5'Rabb-3'mtRNR-EMCV and 5'TPL-3'Biontech) led to the most pronounced increase in the luciferase amount in the in vivo experiment in mice. Subsequent analysis of the stability of the mRNA indicated that the increase in luciferase expression is explained primarily by the efficiency of translation, not by the number of RNA molecules. Altogether, these findings suggest that 5'UTR-and-3'UTR combinations 5'Rabb-3'mtRNR- EMCV and 5'TPL-3'Biontech lead to high expression of target proteins and may be considered for use in preventive and therapeutic modalities based on mRNA.

Deep learning of human polyadenylation sites at nucleotide resolution reveals molecular determinants of site usage and relevance in disease

The genomic distribution of cleavage and polyadenylation (polyA) sites should be co-evolutionally optimized with the local gene structure. Otherwise, spurious polyadenylation can cause premature transcription termination and generate aberrant proteins. To obtain mechanistic insights into polyA site optimization across the human genome, we develop deep/machine learning models to identify genome-wide putative polyA sites at unprecedented nucleotide-level resolution and calculate their strength and usage in the genomic context. Our models quantitatively measure position-specific motif importance and their crosstalk in polyA site formation and cleavage heterogeneity. The intronic site expression is governed by the surrounding splicing landscape. The usage of alternative polyA sites in terminal exons is modulated by their relative locations and distance to downstream genes. Finally, we apply our models to reveal thousands of disease- and trait-associated genetic variants altering polyadenylation activity. Altogether, our models represent a valuable resource to dissect molecular mechanisms mediating genome-wide polyA site expression and characterize their functional roles in human diseases.