Mechanisms and consequences of alternative polyadenylation

A genome-wide perspective of the maternal mRNA translation program during oocyte development

Transcriptional and post-transcriptional regulations control gene expression in most cells. However, critical transitions during the development of the female gamete relies exclusively on regulation of mRNA translation in the absence of de novo mRNA synthesis. Specific temporal patterns of maternal mRNA translation are essential for the oocyte progression through meiosis, for generation of a haploid gamete ready for fertilization and for embryo development. In this review, we will discuss how mRNAs are translated during oocyte growth and maturation using mostly a genome-wide perspective. This broad view on how translation is regulated reveals multiple divergent translational control mechanisms required to coordinate protein synthesis with progression through the meiotic cell cycle and with development of a totipotent zygote.

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.

Nanopore Direct RNA Sequencing Reveals the Short-Term Salt Stress Response in Maize Roots

Transcriptome analysis, relying on the cutting-edge sequencing of cDNA libraries, has become increasingly prevalent within functional genome studies. However, the dependence on cDNA in most RNA sequencing technologies restricts their ability to detect RNA base modifications. To address this limitation, the latest Oxford Nanopore Direct RNA Sequencing (ONT DRS) technology was employed to investigate the transcriptome of maize seedling roots under salt stress. This approach aimed to unveil both the RNA transcriptional profiles and alterations in base modifications. The analysis of the differential expression revealed a total of 1398 genes and 2223 transcripts that exhibited significant variation within the maize root system following brief exposure to salt stress. Enrichment analyses, such as the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway assessments, highlighted the predominant involvement of these differentially expressed genes (DEGs) in regulating ion homeostasis, nitrogen metabolism, amino acid metabolism, and the phytohormone signaling pathways. The protein-protein interaction (PPI) analysis showed the participation of various proteins related to glycolytic metabolism, nitrogen metabolism, amino acid metabolism, abscisic acid signaling, and the jasmonate signaling pathways. It was through this intricate molecular network that these proteins collaborated to safeguard root cells against salt-induced damage. Moreover, under salt stress conditions, the occurrence of variable shear events (AS) in RNA modifications diminished, the average length of poly(A) tails underwent a slight decrease, and the number of genes at the majority of the variable polyadenylation (APA) sites decreased. Additionally, the levels of N5-methylcytosine (m5C) and N6-methyladenosine (m6A) showed a reduction. These results provide insights into the mechanisms of early salt tolerance in maize.

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 landscape of longissimus dorsi muscle with high and low intramuscular fat content in cattle

Intramuscular fat content is one of the most important factors affecting beef quality. However, the role of alternative polyadenylation (APA) in intramuscular fat deposition remains unclear. We compared APA events in muscle samples from high and low intramuscular fat (IMF) cattle, based on RNA-seq data. A total of 363 significant APAs were identified. Notably, the number of shortened 3'UTR events exceeded the number of lengthened 3'UTR events, and genes associated with shortened 3'UTR events were enriched in fatty acid metabolism-related pathways. Most APA events had alternative 3'UTR (aUTR) lengths of 200 to 300 bp. As the 3'UTR lengthened, the aUTR also lengthened (R2 = 0.79). These findings indicate that genes with longer 3'UTRs are more likely to be regulated by APA in the muscle of cattle with high IMF. To determine whether the identified APA events drove alterations in the expression of fat deposition-related genes, we analyzed the relationship between APA events and differentially expressed genes and identified several genes critical for fat deposition (e.g., PFKL and SLC1A5). Since miRNAs usually bind to the 3'UTR region of protein-coding genes and affect gene expression, we constructed an miRNA-APA network to detect several key miRNAs that may regulate fat deposition. We identified 10 important miRNAs that affect changes in IMF content, which may be gained (gained miRNA-binding sites) or lost (lost miRNA-binding sites) owing to 187 differential APA events. Our study characterized the APA profiles of cattle with high and low intramuscular fat content and provided further insights into the relationship between APA, miRNA, and fat deposition.

hnRNPA1 SUMOylation promotes cold hypersensitivity in chronic inflammatory pain by stabilizing TRPA1 mRNA

TRPA1 is pivotal in cold hypersensitivity, but its regulatory mechanisms in inflammatory cold hyperalgesia remain poorly understood. We show here that the upregulation of SUMO1-conjugated protein levels in a complete Freund's adjuvant (CFA)-induced inflammatory pain model enhances TRPA1 mRNA stability, ultimately leading to increased expression levels. We further demonstrate that hnRNPA1 binds to TRPA1 mRNA, and its SUMOylation, upregulated in CFA-induced inflammatory pain, contributes to stabilizing TRPA1 mRNA by accumulating hnRNPA1 in the cytoplasm. Moreover, we find that wild-type hnRNPA1 viral infection in dorsal root ganglia neurons, and not infection with the SUMOylation-deficient hnRNPA1 mutant, can rescue the reduced ability of hnRNPA1-knockdown mice to develop inflammatory cold pain hypersensitivity. These results suggest that hnRNPA1 is a regulator of TRPA1 mRNA stability, the capability of which is enhanced upon SUMO1 conjugation at lysine 3 in response to peripheral inflammation, and the increased expression of TRPA1 in turn underlies the development of chronic inflammatory cold pain hypersensitivity.

Long non-coding RNA generated from CDKN1A gene by alternative polyadenylation regulates p21 expression during DNA damage response

Alternative Polyadenylation (APA) is an emerging mechanism for dynamic changes in gene expression. Previously, we described widespread APA occurrence in introns during the DNA damage response (DDR). Here, we show that a DDR-activated APA event occurs in the first intron of CDKN1A, inducing an alternate last exon-containing lncRNA. We named this lncRNA SPUD (Selective Polyadenylation Upon DNA Damage). SPUD localizes to polysomes in the cytoplasm and is detectable as multiple isoforms in available high-throughput studies. SPUD has low abundance compared to the CDKN1A full-length isoform under non-stress conditions, and SPUD is induced in cancer and normal cells under a variety of DNA damaging conditions in part through p53. The RNA binding protein HuR binds to and promotes the stability of SPUD precursor RNA. SPUD induction increases p21 protein, but not mRNA levels, affecting p21 functions in cell-cycle, CDK2 expression and cell growth. Like CDKN1A full-length isoform, SPUD can bind two competitive p21 translational regulators, the inhibitor calreticulin and the activator CUGBP1; SPUD alters their association with CDKN1A full-length in a DDR-dependent manner, promoting CDKN1A translation. Together, these results show a new regulatory mechanism by which a lncRNA controls p21 expression post-transcriptionally, highlighting lncRNA relevance in DDR progression and cell-cycle.

Neuroepigenetics of ageing and neurodegeneration-associated dementia: An updated review

Gene expression is tremendously altered in the brain during memory acquisition, recall, and forgetfulness. However, non-genetic factors, including environmental elements, epigenetic changes, and lifestyle, have grabbed significant attention in recent years regarding the etiology of neurodegenerative diseases (NDD) and age-associated dementia. Epigenetic modifications are essential in regulating gene expression in all living organisms in a DNA sequence-independent manner. The genes implicated in ageing and NDD-related memory disorders are epigenetically regulated by processes such as DNA methylation, histone acetylation as well as messenger RNA editing machinery. The physiological and optimal state of the epigenome, especially within the CNS of humans, plays an intricate role in helping us adjust to the changing environment, and alterations in it cause many brain disorders, but the mechanisms behind it still need to be well understood. When fully understood, these epigenetic landscapes could act as vital targets for pharmacogenetic rescue strategies for treating several diseases, including neurodegeneration- and age-induced dementia. Keeping this objective in mind, this updated review summarises the epigenetic changes associated with age and neurodegeneration-associated dementia.

On a path toward a broad-spectrum anti-viral: inhibition of HIV-1 and coronavirus replication by SR kinase inhibitor harmine

IMPORTANCE

This study highlights the crucial role RNA processing plays in regulating viral gene expression and replication. By targeting SR kinases, we identified harmine as a potent inhibitor of HIV-1 as well as coronavirus (HCoV-229E and multiple SARS-CoV-2 variants) replication. Harmine inhibits HIV-1 protein expression and reduces accumulation of HIV-1 RNAs in both cell lines and primary CD4+ T cells. Harmine also suppresses coronavirus replication post-viral entry by preferentially reducing coronavirus sub-genomic RNA accumulation. By focusing on host factors rather than viral targets, our study offers a novel approach to combating viral infections that is effective against a range of unrelated viruses. Moreover, at doses required to inhibit virus replication, harmine had limited toxicity and minimal effect on the host transcriptome. These findings support the viability of targeting host cellular processes as a means of developing broad-spectrum anti-virals.

Alternative polyadenylation regulates acetyl-CoA carboxylase function in peanut

BACKGROUND

Polyadenylation is a crucial process that terminates mRNA molecules at their 3'-ends. It has been observed that alternative polyadenylation (APA) can generate multiple transcripts from a single gene locus, each with different polyadenylation sites (PASs). This leads to the formation of several 3' untranslated regions (UTRs) that vary in length and composition. APA has a significant impact on approximately 60-70% of eukaryotic genes and has far-reaching implications for cell proliferation, differentiation, and tumorigenesis.

RESULTS

In this study, we conducted long-read, single-molecule sequencing of mRNA from peanut seeds. Our findings revealed that over half of all peanut genes possess over two PASs, with older developing seeds containing more PASs. This suggesting that the PAS exhibits high tissue specificity and plays a crucial role in peanut seed maturation. For the peanut acetyl-CoA carboxylase A1 (AhACCA1) gene, we discovered four 3' UTRs referred to UTR1-4. RT-PCR analysis showed that UTR1-containing transcripts are predominantly expressed in roots, leaves, and early developing seeds. Transcripts containing UTR2/3 accumulated mainly in roots, flowers, and seeds, while those carrying UTR4 were constitutively expressed. In Nicotiana benthamiana leaves, we transiently expressed all four UTRs, revealing that each UTR impacted protein abundance but not subcellular location. For functional validation, we introduced each UTR into yeast cells and found UTR2 enhanced AhACCA1 expression compared to a yeast transcription terminator, whereas UTR3 did not. Furthermore, we determined ACC gene structures in seven plant species and identified 51 PASs for 15 ACC genes across four plant species, confirming that APA of the ACC gene family is universal phenomenon in plants.

CONCLUSION

Our data demonstrate that APA is widespread in peanut seeds and plays vital roles in peanut seed maturation. We have identified four 3' UTRs for AhACCA1 gene, each showing distinct tissue-specific expression patterns. Through subcellular location experiment and yeast transformation test, we have determined that UTR2 has a stronger impact on gene expression regulation compared to the other three UTRs.

Downregulation of apoptotic repressor AVEN exacerbates cardiac injury after myocardial infarction

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

Novel quinazolin-4(3H)-one based Cyclin K degraders regulate alternative polyadenylation activity

Phenotypic screening is gaining attention as a powerful method for identifying compounds that regulate cellular phenotypes of interest through novel mechanisms of action. Recently, a new modality of compounds, called molecular glues, which can induce the degradation of target proteins by forming ternary complexes of E3 ligases, has emerged from phenotypic screening. In this study, using global proteomic analysis, we identified a novel Cyclin K degrader, T4, which was previously discovered through phenotypic screening for alternative polyadenylation regulation. Our detailed mechanistic analysis revealed that T4 induced Cyclin K degradation, leading to the regulation of alternative polyadenylation. Additionally, we generated a more potent Cyclin K degrader, TR-213, through a structure-activity relationship study of T4. T4 and TR-213 are structurally distinct from other Cyclin K degraders and can be used as novel chemical tools to further analyze the degradation of Cyclin K and the regulation of alternative polyadenylation.

Recovery of missing single-cell RNA-sequencing data with optimized transcriptomic references

Single-cell RNA-sequencing (scRNA-seq) is an indispensable tool for characterizing cellular diversity and generating hypotheses throughout biology. Droplet-based scRNA-seq datasets often lack expression data for genes that can be detected with other methods. Here we show that the observed sensitivity deficits stem from three sources: (1) poor annotation of 3' gene ends; (2) issues with intronic read incorporation; and (3) gene overlap-derived read loss. We show that missing gene expression data can be recovered by optimizing the reference transcriptome for scRNA-seq through recovering false intergenic reads, implementing a hybrid pre-mRNA mapping strategy and resolving gene overlaps. We demonstrate, with a diverse collection of mouse and human tissue data, that reference optimization can substantially improve cellular profiling resolution and reveal missing cell types and marker genes. Our findings argue that transcriptomic references need to be optimized for scRNA-seq analysis and warrant a reanalysis of previously published datasets and cell atlases.

Bovine C-X-C Motif Chemokine Ligand 14 Expression Is Regulated by Alternative Polyadenylation and MicroRNAs

Alternative polyadenylation (APA), including APA that occurs only in the 3' UTR (3' UTR-APA), is an important post-transcriptional regulatory mechanism that leads to distinct 3' UTRs for some genes, increasing the complexity of the transcriptome. The post-transcriptional events regulating the expression of bovine, the C-X-C motif chemokine ligand 14 (CXCL14) gene, remains largely unknown. Here, we find that the bovine CXCL14 gene produces two different lengths of mRNA isoforms due to 3' UTR-APA, and the short and long 3' UTR is 126 bp and 1155 bp, respectively. We found that the expression level of the short isoform was significantly higher than that of the long isoform by luciferase assays and overexpression of different CXCL14 3' UTR-APA isoforms. Moreover, using luciferase assay and site-directed mutagenesis experiments, the results showed that the long CXCL14 3' UTR-APA isoform is downregulated by miR-17-5p, miR-150, and miR-217. However, because the short isoform lacks the true target of miR-17-5p, miR-150, and miR-217 in its 3' UTR and thus escapes the inhibitory effect of these microRNAs, its expression level is significantly higher than that of the long isoform. Finally, we demonstrate that the short CXCL14 3' UTR-APA isoform promotes preadipocyte proliferation by cell counting kit 8 (CCK8) assays. Collectively, our results show that the CXCL14 gene is post-transcriptionally regulated through APA and microRNAs.

Divergent molecular events underlying initial T-cell commitment in human prenatal and postnatal thymus

Introduction

Intrathymic T-cell development is a coordinated process accompanied by dynamic changes in gene expression. Although the transcriptome characteristics of developing T cells in both human fetal and postnatal thymus at single-cell resolution have been revealed recently, the differences between human prenatal and postnatal thymocytes regarding the ontogeny and early events of T-cell development still remain obscure. Moreover, the transcriptional heterogeneity and posttranscriptional gene expression regulation such as alternative polyadenylation at different stages are also unknown.

Method

In this study, we performed integrative single-cell analyses of thymocytes at distinct developmental stages.

Results

The subsets of prenatal CD4-CD8- double-negative (DN) cells, the most immature thymocytes responsible for T-cell lineage commitment, were characterized. By comprehensively comparing prenatal and postnatal DN cells, we revealed significant differences in some key gene expressions. Specifically, prenatal DN subpopulations exhibited distinct biological processes and markedly activated several metabolic programs that may be coordinated to meet the required bioenergetic demands. Although showing similar gene expression patterns along the developmental path, prenatal and postnatal thymocytes were remarkably varied regarding the expression dynamics of some pivotal genes for cell cycle, metabolism, signaling pathway, thymus homing, and T-cell commitment. Finally, we quantified the transcriptome-wide changes in alternative polyadenylation across T-cell development and found diverse preferences of polyadenylation site usage in divergent populations along the T-cell commitment trajectory.

Discussion

In summary, our results revealed transcriptional heterogeneity and a dynamic landscape of alternative polyadenylation during T-cell development in both human prenatal and postnatal thymus, providing a comprehensive resource for understanding T lymphopoiesis in human thymus.

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

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

Coupling of co-transcriptional splicing and 3' end Pol II pausing during termination in Arabidopsis

BACKGROUND

In Arabidopsis, RNA Polymerase II (Pol II) often pauses within a few hundred base pairs downstream of the polyadenylation site, reflecting efficient transcriptional termination, but how such pausing is regulated remains largely elusive.

RESULT

Here, we analyze Pol II dynamics at 3' ends by combining comprehensive experiments with mathematical modelling. We generate high-resolution serine 2 phosphorylated (Ser2P) Pol II positioning data specifically enriched at 3' ends and define a 3' end pause index (3'PI). The position but not the extent of the 3' end pause correlates with the termination window size. The 3'PI is not decreased but even mildly increased in the termination deficient mutant xrn3, indicating 3' end pause is a regulatory step early during the termination and before XRN3-mediated RNA decay that releases Pol II. Unexpectedly, 3'PI is closely associated with gene exon numbers and co-transcriptional splicing efficiency. Multiple exons genes often display stronger 3' end pauses and more efficient on-chromatin splicing than genes with fewer exons. Chemical inhibition of splicing strongly reduces the 3'PI and disrupts its correlation with exon numbers but does not globally impact 3' end readthrough levels. These results are further confirmed by fitting Pol II positioning data with a mathematical model, which enables the estimation of parameters that define Pol II dynamics.

CONCLUSION

Our work highlights that the number of exons via co-transcriptional splicing is a major determinant of Pol II pausing levels at the 3' end of genes in plants.

Development of Nanopore sequencing-based full-length transcriptome database toward functional genome annotation of the Pacific oyster, Crassostrea gigas

The Pacific oyster (Crassostrea gigas) is a widely cultivated shellfish in the world, while its transcriptome diversity remains less unexplored due to the limitation of short reads. In this study, we used Oxford Nanopore sequencing to develop the full-length transcriptome database of C. gigas. We identified 77,920 full-length transcripts from 21,523 genes, and uncovered 9668 alternative splicing events and 87,468 alternative polyadenylation sites. Notably, a total of 16,721 novel transcripts were annotated in this work. Furthermore, integrative analysis of 25 publicly available RNA-seq datasets revealed the transcriptome diversity involved in post-transcriptional regulation in C. gigas. We further developed a Drupal based webserver, Cgtdb, which can be used for transcriptome visualization, sequence alignment, and functional genome annotation analyses. This work provides valuable resources and a useful tool for integrative analysis of various transcriptome datasets in C. gigas, which will serve as an essential reference for functional annotation of the oyster genome.

Inhibitor AN3661 reveals biological functions of Arabidopsis CLEAVAGE and POLYADENYLATION SPECIFICITY FACTOR 73

Cleavage and polyadenylation specificity factor (CPSF) is a protein complex that plays an essential biochemical role in mRNA 3'-end formation, including poly(A) signal recognition and cleavage at the poly(A) site. However, its biological functions at the organismal level are mostly unknown in multicellular eukaryotes. The study of plant CPSF73 has been hampered by the lethality of Arabidopsis (Arabidopsis thaliana) homozygous mutants of AtCPSF73-I and AtCPSF73-II. Here, we used poly(A) tag sequencing to investigate the roles of AtCPSF73-I and AtCPSF73-II in Arabidopsis treated with AN3661, an antimalarial drug with specificity for parasite CPSF73 that is homologous to plant CPSF73. Direct seed germination on an AN3661-containing medium was lethal; however, 7-d-old seedlings treated with AN3661 survived. AN3661 targeted AtCPSF73-I and AtCPSF73-II, inhibiting growth through coordinating gene expression and poly(A) site choice. Functional enrichment analysis revealed that the accumulation of ethylene and auxin jointly inhibited primary root growth. AN3661 affected poly(A) signal recognition, resulted in lower U-rich signal usage, caused transcriptional readthrough, and increased the distal poly(A) site usage. Many microRNA targets were found in the 3' untranslated region lengthened transcripts; these miRNAs may indirectly regulate the expression of these targets. Overall, this work demonstrates that AtCPSF73 plays important part in co-transcriptional regulation, affecting growth, and development in Arabidopsis.

MicroRNA expression changes in the development of rotator cuff tendon injuries

Traumatic or degenerative rotator cuff (RC) tendon injuries are a leading cause of persistent shoulder pain and reduction of mobility with associated disability and dysfunction, which require each year more than 250,000 surgical repairs in the United States. MicroRNAs (miRNAs) are small noncoding RNAs, that in the posttranscriptional phase lead to the development and function of tissues. The aim of this review was to identify miRNA expression changes in patients with RC pathologies and to determine their relevance as a potential novel diagnostic and potentially therapeutic tool for RC disorders. Various miRNAs seemed to be key regulators in the muscle architecture, determining several modifications in muscle atrophy, skeletal muscle mechanical adaptation, lipid accumulation, and fibrosis in the presence of RC tears. The search was executed using PubMed, Medline, Scopus, and Cochrane Central. We included studies written in English that evaluated the role of miRNA in diagnosis, physiopathology, and potential therapeutic application of RC tendon injuries. We included 11 studies in this review. Many miRNAs emerged as key regulators in the pathogenesis of RC tears, inflammation, and muscle fatty degeneration. In fact, they are involved in the regulation of myogenesis, inflammatory cytokines, metalloproteases expression, muscle adaptation, adipogenesis, fibrogenic factors, and extracellular matrix synthesis. The gene expression may be altered in the pathological processes of tendon lesions. Therefore, the knowledge of all the gene mechanisms underlying RC tendinopathy should be achieved with future diagnostic and clinical studies.

Alternative polyadenylation regulation in cardiac development and cardiovascular disease

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

Transcriptome sequencing suggests that pre-mRNA splicing counteracts widespread intronic cleavage and polyadenylation

Alternative splicing (AS) and alternative polyadenylation (APA) are two crucial steps in the post-transcriptional regulation of eukaryotic gene expression. Protocols capturing and sequencing RNA 3'-ends have uncovered widespread intronic polyadenylation (IPA) in normal and disease conditions, where it is currently attributed to stochastic variations in the pre-mRNA processing. Here, we took advantage of the massive amount of RNA-seq data generated by the Genotype Tissue Expression project (GTEx) to simultaneously identify and match tissue-specific expression of intronic polyadenylation sites with tissue-specific splicing. A combination of computational methods including the analysis of short reads with non-templated adenines revealed that APA events are more abundant in introns than in exons. While the rate of IPA in composite terminal exons and skipped terminal exons expectedly correlates with splicing, we observed a considerable fraction of IPA events that lack AS support and attributed them to spliced polyadenylated introns (SPI). We hypothesize that SPIs represent transient byproducts of a dynamic coupling between APA and AS, in which the spliceosome removes the intron while it is being cleaved and polyadenylated. These findings indicate that cotranscriptional pre-mRNA splicing could serve as a rescue mechanism to suppress premature transcription termination at intronic polyadenylation sites.

stAPAminer: Mining Spatial Patterns of Alternative Polyadenylation for Spatially Resolved Transcriptomic Studies

Alternative polyadenylation (APA) contributes to transcriptome complexity and gene expression regulation and has been implicated in various cellular processes and diseases. Single-cell RNA sequencing (scRNA-seq) has enabled the profiling of APA at the single-cell level; however, the spatial information of cells is not preserved in scRNA-seq. Alternatively, spatial transcriptomics (ST) technologies provide opportunities to decipher the spatial context of the transcriptomic landscape. Pioneering studies have revealed potential spatially variable genes and/or splice isoforms; however, the pattern of APA usage in spatial contexts remains unappreciated. In this study, we developed a toolkit called stAPAminer for mining spatial patterns of APA from spatially barcoded ST data. APA sites were identified and quantified from the ST data. In particular, an imputation model based on the k-nearest neighbors algorithm was designed to recover APA signals, and then APA genes with spatial patterns of APA usage variation were identified. By analyzing well-established ST data of the mouse olfactory bulb (MOB), we presented a detailed view of spatial APA usage across morphological layers of the MOB. We compiled a comprehensive list of genes with spatial APA dynamics and obtained several major spatial expression patterns that represent spatial APA dynamics in different morphological layers. By extending this analysis to two additional replicates of the MOB ST data, we observed that the spatial APA patterns of several genes were reproducible among replicates. stAPAminer employs the power of ST to explore the transcriptional atlas of spatial APA patterns with spatial resolution. This toolkit is available at https://github.com/BMILAB/stAPAminer and https://ngdc.cncb.ac.cn/biocode/tools/BT007320.

A novel neurodevelopmental syndrome caused by loss-of-function of the Zinc Finger Homeobox 3 (ZFHX3) gene

Neurodevelopmental disorders (NDDs) result from impaired development and functioning of the brain. Here, we identify loss-of-function variation in ZFHX3 as a novel cause for syndromic intellectual disability (ID). ZFHX3, previously known as ATBF1, is a zinc-finger homeodomain transcription factor involved in multiple biological processes including cell differentiation and tumorigenesis. Through international collaboration, we collected clinical and morphometric data (Face2Gene) of 41 individuals with protein truncating variants (PTVs) or (partial) deletions of ZFHX3 . We used data mining, RNA and protein analysis to identify the subcellular localization and spatiotemporal expression of ZFHX3 in multiple in vitro models. We identified the DNA targets of ZFHX3 using ChIP seq. Immunoprecipitation followed by mass spectrometry indicated potential binding partners of endogenous ZFHX3 in neural stem cells that were subsequently confirmed by reversed co-immunoprecipitation and western blot. We evaluated a DNA methylation profile associated with ZFHX3 haploinsufficiency using DNA methylation analysis on whole blood extracted DNA of six individuals with ZFHX3 PTVs and four with a (partial) deletion of ZFHX3 . A reversed genetic approach characterized the ZFHX3 orthologue in Drosophila melanogaster . Loss-of-function variation of ZFHX3 consistently associates with (mild) ID and/or behavioural problems, postnatal growth retardation, feeding difficulties, and recognizable facial characteristics, including the rare occurrence of cleft palate. Nuclear abundance of ZFHX3 increases during human brain development and neuronal differentiation in neural stem cells and SH-SY5Y cells, ZFHX3 interacts with the chromatin remodelling BRG1/Brm-associated factor complex and the cleavage and polyadenylation complex. In line with a role for chromatin remodelling, ZFHX3 haploinsufficiency associates with a specific DNA methylation profile in leukocyte-derived DNA. The target genes of ZFHX3 are implicated in neuron and axon development. In Drosophila melanogaster , z fh2, considered to be the ZFHX3 orthologue, is expressed in the third instar larval brain. Ubiquitous and neuron-specific knockdown of zfh2 results in adult lethality underscoring a key role for zfh2 in development and neurodevelopment. Interestingly, ectopic expression of zfh2 as well as ZFHX3 in the developing wing disc results in a thoracic cleft phenotype. Collectively, our data shows that loss-of-function variants in ZFHX3 are a cause of syndromic ID, that associates with a specific DNA methylation profile. Furthermore, we show that ZFHX3 participates in chromatin remodelling and mRNA processing.

A novel scoring model based on RNA modification "writers" can predict the prognosis and guide immunotherapy in gastric cancer

Although extensive research has been carried out on the epigenetic regulation of single RNA modifications in gastric cancer, little is known regarding the crosstalk of four major RNA adenosine modifications, namely, m6A, m1A, alternative polyadenylation and adenosine-to-inosine RNA editing. By analyzing 26 RNA modification "writers" in 1750 gastric cancer samples, we creatively constructed a scoring model called the "Writers" of the RNA Modification Score (WRM_Score), which was able to quantify the RNA modification subtypes of individual patients. In addition, we explored the relationship between WRM_Score and transcriptional and posttranscriptional regulation, tumor microenvironment, clinical features and molecular subtypes. We constructed an RNA modification scoring model including two different subgroups: WRM_Score_low and WRM_Score_high. The former was associated with survival benefit and good efficacy of immune checkpoint inhibitors (ICIs) due to gene repair and immune activation, while the latter was related to poor prognosis and bad efficacy of ICIs because of stromal activation and immunosuppression. The WRM score based on immune and molecular characteristics of the RNA modification pattern is a reliable predictor of the prognosis of gastric cancer and the therapeutic efficacy of immune checkpoint inhibitors in gastric cancer.

Introns: the "dark matter" of the eukaryotic genome

The emergence of introns was a significant evolutionary leap that is a major distinguishing feature between prokaryotic and eukaryotic genomes. While historically introns were regarded merely as the sequences that are removed to produce spliced transcripts encoding functional products, increasingly data suggests that introns play important roles in the regulation of gene expression. Here, we use an intron-centric lens to review the role of introns in eukaryotic gene expression. First, we focus on intron architecture and how it may influence mechanisms of splicing. Second, we focus on the implications of spliceosomal snRNAs and their variants on intron splicing. Finally, we discuss how the presence of introns and the need to splice them influences transcription regulation. Despite the abundance of introns in the eukaryotic genome and their emerging role regulating gene expression, a lot remains unexplored. Therefore, here we refer to introns as the "dark matter" of the eukaryotic genome and discuss some of the outstanding questions in the field.

The ENCODE4 long-read RNA-seq collection reveals distinct classes of transcript structure diversity

The majority of mammalian genes encode multiple transcript isoforms that result from differential promoter use, changes in exonic splicing, and alternative 3' end choice. Detecting and quantifying transcript isoforms across tissues, cell types, and species has been extremely challenging because transcripts are much longer than the short reads normally used for RNA-seq. By contrast, long-read RNA-seq (LR-RNA-seq) gives the complete structure of most transcripts. We sequenced 264 LR-RNA-seq PacBio libraries totaling over 1 billion circular consensus reads (CCS) for 81 unique human and mouse samples. We detect at least one full-length transcript from 87.7% of annotated human protein coding genes and a total of 200,000 full-length transcripts, 40% of which have novel exon junction chains. To capture and compute on the three sources of transcript structure diversity, we introduce a gene and transcript annotation framework that uses triplets representing the transcript start site, exon junction chain, and transcript end site of each transcript. Using triplets in a simplex representation demonstrates how promoter selection, splice pattern, and 3' processing are deployed across human tissues, with nearly half of multi-transcript protein coding genes showing a clear bias toward one of the three diversity mechanisms. Evaluated across samples, the predominantly expressed transcript changes for 74% of protein coding genes. In evolution, the human and mouse transcriptomes are globally similar in types of transcript structure diversity, yet among individual orthologous gene pairs, more than half (57.8%) show substantial differences in mechanism of diversification in matching tissues. This initial large-scale survey of human and mouse long-read transcriptomes provides a foundation for further analyses of alternative transcript usage, and is complemented by short-read and microRNA data on the same samples and by epigenome data elsewhere in the ENCODE4 collection.

Genome-wide alternative polyadenylation dynamics underlying plant growth retardant-induced dwarfing of pomegranate

Dwarfed stature is a desired agronomic trait for pomegranate (Punica granatum L.), with its advantages such as lower cost and increased yield. A comprehensive understanding of regulatory mechanisms underlying the growth repression would provide a genetic foundation to molecular-assisted dwarfing cultivation of pomegranate. Our previous study induced dwarfed pomegranate seedlings via exogenous application of plant growth retardants (PGRs) and highlighted the important roles of differential expression of plant growth-related genes in eliciting the dwarfed phenotype of pomegranate. Alternative polyadenylation (APA) is an important post-transcriptional mechanism and has been demonstrated to act as a key regulator in plant growth and development. However, no attention has been paid to the role of APA in PGR-induced dwarfing in pomegranate. In this study, we characterized and compared APA-mediated regulation events underlying PGR-induced treatments and normal growth condition. Genome-wide alterations in the usage of poly(A) sites were elicited by PGR treatments, and these changes were involved in modulating the growth and development of pomegranate seedlings. Importantly, ample specificities were observed in APA dynamics among the different PGR treatments, which mirrors their distinct nature. Despite the asynchrony between APA events and differential gene expression, APA was found to regulate transcriptome via influencing microRNA (miRNA)-mediated mRNA cleavage or translation inhibition. A global preference for lengthening of 3' untranslated regions (3' UTRs) was observed under PGR treatments, which was likely to host more miRNA target sites in 3' UTRs and thus suppress the expression of the corresponding genes, especially those associated with developmental growth, lateral root branching, and maintenance of shoot apical meristem. Together, these results highlighted the key role of APA-mediated regulations in fine-tuning the PGR-induced dwarfed stature of pomegranate, which provides new insights into the genetic basis underlying the growth and development of pomegranate.

CPSF6-mediated XBP1 3'UTR shortening attenuates cisplatin-induced ER stress and elevates chemo-resistance in lung adenocarcinoma

Alternative polyadenylation (APA) is a widespread mechanism generating RNA molecules with alternative 3' ends. Herein, we discovered that TargetScan includes a novel XBP1 transcript with a longer 3' untranslated region (UTR) (XBP1-UL) than that included in NCBI. XBP1-UL exhibited a lowered level in blood samples from lung adenocarcinoma (LUAD) patients and in those after DDP treatment. Consistently, XBP1-UL was reduced in A549 cells compared to normal BEAS-2B cells, as well as in DDP-treated/resistant A549 cells relative to controls. Moreover, due to decreased usage of the distal polyadenylation site (PAS) in 3'UTR, XBP1-UL level was lowered in A549 cells and decreased further in DDP-resistant A549 (A549/DDP) cells. Importantly, use of the distal PAS (dPAS) and XBP1-UL level were gradually reduced in A549 cells under increasing concentrations of DDP, which was attributed to DDP-induced endoplasmic reticulum (ER) stress. Furthermore, XBP1 transcripts with shorter 3'UTR (XBP1-US) were more stable and presented stronger potentiation on DDP resistance. The choice of proximal PAS (pPAS) was attributed to CPSF6 elevation, which was caused by BRCA1-distrupted R-loop accumulation in CPSF6 5'end. DDP-induced nuclear LINC00221 also facilitated CPSF6-induced pPAS choice in the pre-XBP1 3'end. Finally, we found that unlike the unspliced XBP1 protein (XBP1-u), the spliced form XBP1-s retarded p53 degradation to facilitate DNA damage repair of LUAD cells. The current study provides new insights into tumor progression and DDP resistance in LUAD, which may contribute to improved LUAD treatment.

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

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

Iron restriction increases the expression of a cytotoxic cysteine proteinase TvCP2 by a novel mechanism of tvcp2 mRNA alternative polyadenylation in Trichomonas vaginalis

Trichomonas vaginalis TvCP2 (TVAG_057000) is a cytotoxic cysteine proteinase (CP) expressed under iron-limited conditions. This work aimed to identify one of the mechanisms of tvcp2 gene expression regulation by iron at the posttranscriptional level. We checked tvcp2 mRNA stability under both iron-restricted (IR) and high iron (HI) conditions in the presence of actinomycin D. Greater stability of the tvcp2 mRNA under the IR than in HI conditions was observed, as expected. In silico analysis of the 3' regulatory region showed the presence of two putative polyadenylation signals in the tvcp2 transcript. By 3'-RACE assays, we demonstrated the existence of two isoforms of the tvcp2 mRNA with different 3'-UTR that resulted in more TvCP2 protein under IR than in HI conditions detected by WB assays. Additionally, we searched for homologs of the trichomonad polyadenylation machinery by an in silico analysis in the genome database, TrichDB. 16 genes that encode proteins that could be part of the trichomonad polyadenylation machinery were found. qRT-PCR assays showed that most of these genes were positively regulated by iron. Thus, our results show the presence of alternative polyadenylation as a novel iron posttranscriptional regulatory mechanism in T. vaginalis for the tvcp2 gene expression.

Mechanism of STMN2 cryptic splice-polyadenylation and its correction for TDP-43 proteinopathies

Loss of nuclear TDP-43 is a hallmark of neurodegeneration in TDP-43 proteinopathies, including amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). TDP-43 mislocalization results in cryptic splicing and polyadenylation of pre-messenger RNAs (pre-mRNAs) encoding stathmin-2 (also known as SCG10), a protein that is required for axonal regeneration. We found that TDP-43 binding to a GU-rich region sterically blocked recognition of the cryptic 3' splice site in STMN2 pre-mRNA. Targeting dCasRx or antisense oligonucleotides (ASOs) suppressed cryptic splicing, which restored axonal regeneration and stathmin-2-dependent lysosome trafficking in TDP-43-deficient human motor neurons. In mice that were gene-edited to contain human STMN2 cryptic splice-polyadenylation sequences, ASO injection into cerebral spinal fluid successfully corrected Stmn2 pre-mRNA misprocessing and restored stathmin-2 expression levels independently of TDP-43 binding.

Random Forest model reveals the interaction between N6-methyladenosine modifications and RNA-binding proteins

RNA-binding proteins (RBPs) have critical roles in N6-methyladenosine (m6A) modification process. We designed a Random Forest (RF) model to systematically analyze the interaction among RBPs and m6A modifications by integrating the binding signals from hundreds of RBPs. Accurate prediction of m6A sites demonstrated significant connections between RBP bindings and m6A modifications. The relative importance of different RBPs from the model provided a quantitative metric to evaluate their interactions with m6A modifications. Redundancy analysis showed that several RBPs may have similar binding patterns with m6A sites. The RF model exhibited fairly high prediction accuracy across cell lines, suggesting a conservative RBP interaction network regulates m6A occupancy. Specific RBPs can engage to the corresponding regional m6A sites and deploy distinct regulatory processes, such as cleavage site selection of the alternative polyadenylation (APA). We also integrated histone modifications into our RF model, which demonstrated H3K36me3 and H3K27me3 as determining features for m6A distribution.

CAFuncAPA: a knowledgebase for systematic functional annotations of APA events in human cancers

Alternative polyadenylation (APA) is a widespread posttranscriptional regulation process. APA generates diverse mRNA isoforms with different 3' UTR lengths, affecting mRNA expression, miRNA binding regulation and alternative splicing events. Previous studies have demonstrated the important roles of APA in tumorigenesis and cancer progression through diverse aspects. Thus, a comprehensive functional landscape of diverse APA events would aid in a better understanding of the underlying mechanisms related to APA in human cancers. Here, we built CAFuncAPA (https://relab.xidian.edu.cn/CAFuncAPA/) to systematically annotate the functions of 15478 APA events in human pan-cancers. Specifically, we first identified APA events associated with cancer survival and tumor progression. We annotated the potential downstream effects of APA on genes/isoforms expression, regulation of miRNAs, RNA binding proteins (RBPs) and alternative splicing events. Moreover, we also identified up-regulators of APA events, including the effects of genetic variants on poly(A) sites and RBPs, as well as the effect of methylation phenotypes on APA events. These findings suggested that CAFuncAPA can be a helpful resource for a better understanding of APA regulators and potential functions in cancer biology.

Alternative polyadenylation alters protein dosage by switching between intronic and 3'UTR sites

Alternative polyadenylation (APA) creates distinct transcripts from the same gene by cleaving the pre-mRNA at poly(A) sites that can lie within the 3' untranslated region (3'UTR), introns, or exons. Most studies focus on APA within the 3'UTR; however, here, we show that CPSF6 insufficiency alters protein levels and causes a developmental syndrome by deregulating APA throughout the transcript. In neonatal humans and zebrafish larvae, CPSF6 insufficiency shifts poly(A) site usage between the 3'UTR and internal sites in a pathway-specific manner. Genes associated with neuronal function undergo mostly intronic APA, reducing their expression, while genes associated with heart and skeletal function mostly undergo 3'UTR APA and are up-regulated. This suggests that, under healthy conditions, cells toggle between internal and 3'UTR APA to modulate protein expression.

Alternative Splicing in the Regulatory Circuit of Plant Temperature Response

As sessile organisms, plants have evolved complex mechanisms to rapidly respond to ever-changing ambient temperatures. Temperature response in plants is modulated by a multilayer regulatory network, including transcriptional and post-transcriptional regulations. Alternative splicing (AS) is an essential post-transcriptional regulatory mechanism. Extensive studies have confirmed its key role in plant temperature response, from adjustment to diurnal and seasonal temperature changes to response to extreme temperatures, which has been well documented by previous reviews. As a key node in the temperature response regulatory network, AS can be modulated by various upstream regulations, such as chromatin modification, transcription rate, RNA binding proteins, RNA structure and RNA modifications. Meanwhile, a number of downstream mechanisms are affected by AS, such as nonsense-mediated mRNA decay (NMD) pathway, translation efficiency and production of different protein variants. In this review, we focus on the links between splicing regulation and other mechanisms in plant temperature response. Recent advances regarding how AS is regulated and the following consequences in gene functional modulation in plant temperature response will be discussed. Substantial evidence suggests that a multilayer regulatory network integrating AS in plant temperature response has been unveiled.

CPA-Perturb-seq: Multiplexed single-cell characterization of alternative polyadenylation regulators

Most mammalian genes have multiple polyA sites, representing a substantial source of transcript diversity that is governed by the cleavage and polyadenylation (CPA) regulatory machinery. To better understand how these proteins govern polyA site choice we introduce CPA-Perturb-seq, a multiplexed perturbation screen dataset of 42 known CPA regulators with a 3' scRNA-seq readout that enables transcriptome-wide inference of polyA site usage. We develop a statistical framework to specifically identify perturbation-dependent changes in intronic and tandem polyadenylation, and discover modules of co-regulated polyA sites exhibiting distinct functional properties. By training a multi-task deep neural network (APARENT-Perturb) on our dataset, we delineate a cis-regulatory code that predicts responsiveness to perturbation and reveals interactions between distinct regulatory complexes. Finally, we leverage our framework to re-analyze published scRNA-seq datasets, identifying new regulators that affect the relative abundance of alternatively polyadenylated transcripts, and characterizing extensive cellular heterogeneity in 3' UTR length amongst antibody-producing cells. Our work highlights the potential for multiplexed single-cell perturbation screens to further our understanding of post-transcriptional regulation in vitro and in vivo.

Alternative polyadenylation transcriptome-wide association study identifies APA-linked susceptibility genes in brain disorders

Alternative polyadenylation (APA) plays an essential role in brain development; however, current transcriptome-wide association studies (TWAS) largely overlook APA in nominating susceptibility genes. Here, we performed a 3' untranslated region (3'UTR) APA TWAS (3'aTWAS) for 11 brain disorders by combining their genome-wide association studies data with 17,300 RNA-seq samples across 2,937 individuals. We identified 354 3'aTWAS-significant genes, including known APA-linked risk genes, such as SNCA in Parkinson's disease. Among these 354 genes, ~57% are not significant in traditional expression- and splicing-TWAS studies, since APA may regulate the translation, localization and protein-protein interaction of the target genes independent of mRNA level expression or splicing. Furthermore, we discovered ATXN3 as a 3'aTWAS-significant gene for amyotrophic lateral sclerosis, and its modulation substantially impacted pathological hallmarks of amyotrophic lateral sclerosis in vitro. Together, 3'aTWAS is a powerful strategy to nominate important APA-linked brain disorder susceptibility genes, most of which are largely overlooked by conventional expression and splicing analyses.

Macrophage differentiation is marked by increased abundance of the mRNA 3' end processing machinery, altered poly(A) site usage, and sensitivity to the level of CstF64

Regulation of mRNA polyadenylation is important for response to external signals and differentiation in several cell types, and results in mRNA isoforms that vary in the amount of coding sequence or 3' UTR regulatory elements. However, its role in differentiation of monocytes to macrophages has not been investigated. Macrophages are key effectors of the innate immune system that help control infection and promote tissue-repair. However, overactivity of macrophages contributes to pathogenesis of many diseases. In this study, we show that macrophage differentiation is characterized by shortening and lengthening of mRNAs in relevant cellular pathways. The cleavage/polyadenylation (C/P) proteins increase during differentiation, suggesting a possible mechanism for the observed changes in poly(A) site usage. This was surprising since higher C/P protein levels correlate with higher proliferation rates in other systems, but monocytes stop dividing after induction of differentiation. Depletion of CstF64, a C/P protein and known regulator of polyadenylation efficiency, delayed macrophage marker expression, cell cycle exit, attachment, and acquisition of structural complexity, and impeded shortening of mRNAs with functions relevant to macrophage biology. Conversely, CstF64 overexpression increased use of promoter-proximal poly(A) sites and caused the appearance of differentiated phenotypes in the absence of induction. Our findings indicate that regulation of polyadenylation plays an important role in macrophage differentiation.

Long Non-Coding RNA Generated from CDKN1A Gene by Alternative Polyadenylation Regulates p21 Expression during DNA Damage Response

Alternative Polyadenylation (APA) is an emerging mechanism for dynamic changes in gene expression. Previously, we described widespread APA occurrence in introns during the DNA damage response (DDR). Here, we show that a DNA damage activated APA event occurs in the first intron of CDKN1A , inducing an alternate last exon (ALE)-containing lncRNA. We named this lncRNA SPUD (Selective Polyadenylation Upon Damage). SPUD localizes to polysomes in the cytoplasm and is detectable as multiple isoforms in available high throughput studies. SPUD has low abundance compared to the CDKN1A full-length isoform and is induced in cancer and normal cells under a variety of DNA damaging conditions in part through p53 transcriptional activation. RNA binding protein (RBP) HuR and the transcriptional repressor CTCF regulate SPUD levels. SPUD induction increases p21 protein, but not CDKN1A full-length levels, affecting p21 functions in cell-cycle, CDK2 expression, and cell viability. Like CDKN1A full-length isoform, SPUD can bind two competitive p21 translational regulators, the inhibitor calreticulin and the activator CUGBP1; SPUD can change their association with CDKN1A full-length in a DDR-dependent manner. Together, these results show a new regulatory mechanism by which a lncRNA controls p21 expression post-transcriptionally, highlighting lncRNA relevance in DDR progression and cellcycle.

Measuring the tail: Methods for poly(A) tail profiling

The 3'-end poly(A) tail is an important and potent feature of most mRNA molecules that affects mRNA fate and translation efficiency. Polyadenylation is a posttranscriptional process that occurs in the nucleus by canonical poly(A) polymerases (PAPs). In some specific instances, the poly(A) tail can also be extended in the cytoplasm by noncanonical poly(A) polymerases (ncPAPs). This epitranscriptomic regulation of mRNA recently became one of the most interesting aspects in the field. Advances in RNA sequencing technologies and software development have allowed the precise measurement of poly(A) tails, identification of new ncPAPs, expansion of the function of known enzymes, discovery and a better understanding of the physiological role of tail heterogeneity, and recognition of a correlation between tail length and RNA translatability. Here, we summarize the development of polyadenylation research methods, including classic low-throughput approaches, Illumina-based genome-wide analysis, and advanced state-of-art techniques that utilize long-read third-generation sequencing with Pacific Biosciences and Oxford Nanopore Technologies platforms. A boost in technical opportunities over recent decades has allowed a better understanding of the regulation of gene expression at the mRNA level. This article is categorized under: RNA Methods > RNA Analyses In Vitro and In Silico.

A novel RNA modification prognostic signature for predicting the characteristics of the tumor microenvironment in gastric cancer

Gastric cancer (GC) is one of the most common neoplastic malignancies, which permutes a fourth of cancer-related mortality globally. RNA modification plays a significant role in tumorigenesis, the underlying molecular mechanism of how different RNA modifications directly affect the tumor microenvironment (TME) in GC is unclear. Here, we profiled the genetic and transcriptional alterations of RNA modification genes (RMGs) in GC samples from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) cohorts. Through the unsupervised clustering algorithm, we identified three distinct RNA modification clusters and found that they participate in different biological pathways and starkly correlate with the clinicopathological characteristics, immune cell infiltration, and prognosis of GC patients. Subsequently, univariate Cox regression analysis unveiled 298 of 684 subtype-related differentially expressed genes (DEGs) are tightly interwoven to prognosis. In addition, we conducted the principal component analysis to develop the RM_Score system, which was used to quantify and predict the prognostic value of RNA modification in GC. Our analysis indicated that patients with high RM_Score were characterized by higher tumor mutational burden, mutation frequency, and microsatellite instability which were more susceptible to immunotherapy and had a favorable prognosis. Altogether, our study uncovered RNA modification signatures that may have a potential role in the TME and prediction of clinicopathological characteristics. Identification of these RNA modifications may provide a new understanding of immunotherapy strategies for gastric cancer.

RNA epitranscriptomics dysregulation: A major determinant for significantly increased risk of ASD pathogenesis

Autism spectrum disorders (ASDs) are perhaps the most severe, intractable and challenging child psychiatric disorders. They are complex, pervasive and highly heterogeneous and depend on multifactorial neurodevelopmental conditions. Although the pathogenesis of autism remains unclear, it revolves around altered neurodevelopmental patterns and their implications for brain function, although these cannot be specifically linked to symptoms. While these affect neuronal migration and connectivity, little is known about the processes that lead to the disruption of specific laminar excitatory and inhibitory cortical circuits, a key feature of ASD. It is evident that ASD has multiple underlying causes and this multigenic condition has been considered to also dependent on epigenetic effects, although the exact nature of the factors that could be involved remains unclear. However, besides the possibility for differential epigenetic markings directly affecting the relative expression levels of individual genes or groups of genes, there are at least three mRNA epitranscriptomic mechanisms, which function cooperatively and could, in association with both genotypes and environmental conditions, alter spatiotemporal proteins expression patterns during brain development, at both quantitative and qualitative levels, in a tissue-specific, and context-dependent manner. As we have already postulated, sudden changes in environmental conditions, such as those conferred by maternal inflammation/immune activation, influence RNA epitranscriptomic mechanisms, with the combination of these processes altering fetal brain development. Herein, we explore the postulate whereby, in ASD pathogenesis, RNA epitranscriptomics might take precedence over epigenetic modifications. RNA epitranscriptomics affects real-time differential expression of receptor and channel proteins isoforms, playing a prominent role in central nervous system (CNS) development and functions, but also RNAi which, in turn, impact the spatiotemporal expression of receptors, channels and regulatory proteins irrespective of isoforms. Slight dysregulations in few early components of brain development, could, depending upon their extent, snowball into a huge variety of pathological cerebral alterations a few years after birth. This may very well explain the enormous genetic, neuropathological and symptomatic heterogeneities that are systematically associated with ASD and psychiatric disorders at large.

Association of RNA-modification "writer" genes with prognosis and response to immunotherapy in patients with low-grade glioma

RNA modification is a key regulatory mechanism involved in tumorigenesis, tumor progression, and the immune response. However, the potential role of RNA modification "writer" genes in the immune microenvironment of gliomas and their effect on the response to immunotherapy remains unclear. The purpose of this study was to evaluate the role of RNA modification "writer" gene in the prognosis and immunotherapy response of low-grade glioma (LGG). The consensus non-negative matrix factorization (CNMF) method was used to identify different RNA modification subtypes. We used a novel eigengene screening method, the variable neighborhood learning Harris Hawks optimizer (VNLHHO), to screen for eigengenes among the RNA modification subtypes. We constructed a principal components analysis score(PCA_score)-based prognostic prediction model and validated it using an independent cohort. We also analyzed the association between PCA_score and the immune and molecular features of LGG. The results suggested that LGG can be divided into two different RNA modification-based subtypes with distinct prognostic and molecular features. High PCA_score was significantly associated with a poor prognosis in LGG and was an independent prognostic factor. A nomogram containing PCA_score and clinical features was constructed, and it showed a significant predictive value. PCA_score was negatively correlated with tumor purity and the abundance of CD4+ T cells in LGG patients. LGG patients with high PCA_score had lower Tumor Immune Dysfunction and Exclusion scores and showed an immunotherapy response. In conclusion, we report a novel RNA modification-based prognostic model for LGG that lays the foundation for evaluating LGG prognosis and developing more effective therapeutic strategies for these tumors.

RNA and Protein Interactomes of an RNA-Binding Protein Tagged with FLAG Epitopes Using Combinatory Approaches of Genome Engineering and Stable Transfection

To study the function of RNA-binding proteins (RBPs), an overexpression or knockout approach is generally used. However, as many RBPs are essential to cellular functions, the complete knockout of these proteins may be lethal to the cell. Overexpression of RBPs, on the other hand, may create an altered transcriptome and aberrant phenotypes that can mask their physiological function. Additionally, biochemical characterization of RBP often requires highly specific antibodies for efficient immunoprecipitation for downstream mass spectrometry or RNA footprinting profiling. To overcome these hurdles, we have developed a strategy to generate cellular systems either using a CRISPR-Cas9-mediated epitope tag knock-in approach or a two-step workflow to first stably express an exogenous Flag-tagged RBP and subsequently knockout the endogenous RBP using CRISPR-Cas9 gene editing. The generation of these cell lines will be beneficial for downstream RNA footprinting studies and mass spectrometry-mediated interactome studies.

Chicken GLUT4 undergoes complex alternative splicing events and its expression in striated muscle changes dramatically during development

Glucose transporter protein 4 (GLUT4) plays an important role in regulating insulin-mediated glucose homeostasis in mammals. Until now, studies on GLUT4 have focused on mammals mostly, while chicken GLUT4 has been rarely investigated. In this study, chicken GLUT4 mRNA sequences were obtained by combining conventional amplification, 5'- and 3'- rapid amplification of cDNA ends technique (RACE), then bioinformatics analysis on its genomic structure, splicing pattern, subcellular localization prediction and homologous comparisons were carried out. In addition, the distribution of GLUT4 was detected by RT-qPCR in bird's liver and striated muscles (cardiac muscle, pectoralis and leg muscle) at different ages, including embryonic day 14 (E14), E19, 7-day-old (D7), D21 and D49 (n = 3-4). Results showed that chicken GLUT4 gene produced at least 14 transcripts (GenBank accession No: OP491293-OP491306) through alternative splicing and polyadenylation, which predicted encoding 12 types of amino acid (AA) sequences (with length ranged from 65 AA to 519 AA). These proteins contain typical major facilitator superfamily domain of glucose transporters with length variations, sharing a common sequence of 59 AA, and were predicted to have distinct subcellular localization. The dominant transcript (named as T1) consists of 11 exons with an open reading frame being predicted encoding 519 AA. In addition, analyzing on the spatio-temporal expression of chicken GLUT4 showed it dominantly expressed in pectoralis, leg muscles and cardiac muscle, and the mRNA level of chicken GLUT4 dramatically fluctuated with birds' development in cardiac muscle, pectoralis and leg muscles, with the level at D21 significantly higher than that at E14, E19, and D49 (P < 0.05). These data indicated that chicken GLUT4 undergoes complex alternative splicing events, and GLUT4 expression level in striated muscle was subjected to dynamic regulation with birds' development. Results indicate these isoforms may play overlapping and distinct roles in chicken.

Slight Variations in the Sequence Downstream of the Polyadenylation Signal Significantly Increase Transgene Expression in HEK293T and CHO Cells

Compared to transcription initiation, much less is known about transcription termination. In particular, large-scale mutagenesis studies have, so far, primarily concentrated on promoter and enhancer, but not terminator sequences. Here, we used a massively parallel reporter assay (MPRA) to systematically analyze the influence of short (8 bp) sequence variants (mutations) located downstream of the polyadenylation signal (PAS) on the steady-state mRNA level of the upstream gene, employing an eGFP reporter and human HEK293T cells as a model system. In total, we evaluated 227,755 mutations located at different overlapping positions within +17..+56 bp downstream of the PAS for their ability to regulate the reporter gene expression. We found that the positions +17..+44 bp downstream of the PAS are more essential for gene upregulation than those located more distal to the PAS, and that the mutation sequences ensuring high levels of eGFP mRNA expression are extremely T-rich. Next, we validated the positive effect of a couple of mutations identified in the MPRA screening on the eGFP and luciferase protein expression. The most promising mutation increased the expression of the reporter proteins 13-fold and sevenfold on average in HEK293T and CHO cells, respectively. Overall, these findings might be useful for further improving the efficiency of production of therapeutic products, e.g., recombinant antibodies.

YMLA: A comparative platform to carry out functional enrichment analysis for multiple gene lists in yeast

Comparative analysis among multiple gene lists on their functional features is now a routine task due to the advancement of high-throughput experiments. Several enrichment analysis tools were developed in the past. However, these tools mainly focus on one gene list and contain only gene ontology or interaction features. What makes it worse, comparative investigation and customized feature set reanalysis are still unavailable. Therefore, we constructed the YMLA (Yeast Multiple List Analyzer) platform in this research. YMLA includes 39 yeast features and facilitates comparative analysis among multiple gene lists via tabular views, heatmaps, and network plots. Moreover, the customized feature set reanalysis function was implemented in YMLA to help form mechanism hypotheses based on a selected enriched feature subset. We demonstrated the biological applicability of YMLA via example lists consisting of genes with top/bottom translation efficiency values. The analysis results provided by YMLA reveal novel facts consistent with previous experiments. YMLA is available at https://cosbi7.ee.ncku.edu.tw/YMLA/.

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

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

Deciphering the impact of genetic variation on human polyadenylation using APARENT2

BACKGROUND

3'-end processing by cleavage and polyadenylation is an important and finely tuned regulatory process during mRNA maturation. Numerous genetic variants are known to cause or contribute to human disorders by disrupting the cis-regulatory code of polyadenylation signals. Yet, due to the complexity of this code, variant interpretation remains challenging.

RESULTS

We introduce a residual neural network model, APARENT2, that can infer 3'-cleavage and polyadenylation from DNA sequence more accurately than any previous model. This model generalizes to the case of alternative polyadenylation (APA) for a variable number of polyadenylation signals. We demonstrate APARENT2's performance on several variant datasets, including functional reporter data and human 3' aQTLs from GTEx. We apply neural network interpretation methods to gain insights into disrupted or protective higher-order features of polyadenylation. We fine-tune APARENT2 on human tissue-resolved transcriptomic data to elucidate tissue-specific variant effects. By combining APARENT2 with models of mRNA stability, we extend aQTL effect size predictions to the entire 3' untranslated region. Finally, we perform in silico saturation mutagenesis of all human polyadenylation signals and compare the predicted effects of [Formula: see text] million variants against gnomAD. While loss-of-function variants were generally selected against, we also find specific clinical conditions linked to gain-of-function mutations. For example, we detect an association between gain-of-function mutations in the 3'-end and autism spectrum disorder. To experimentally validate APARENT2's predictions, we assayed clinically relevant variants in multiple cell lines, including microglia-derived cells.

CONCLUSIONS

A sequence-to-function model based on deep residual learning enables accurate functional interpretation of genetic variants in polyadenylation signals and, when coupled with large human variation databases, elucidates the link between functional 3'-end mutations and human health.