Alternative cleavage and polyadenylation: the long and short of it

Modulation of diverse biological processes by CPSF, the master regulator of mRNA 3' ends

The cleavage and polyadenylation specificity factor (CPSF) complex plays a central role in the formation of mRNA 3' ends, being responsible for the recognition of the poly(A) signal sequence, the endonucleolytic cleavage step, and recruitment of poly(A) polymerase. CPSF has been extensively studied for over three decades, and its functions and those of its individual subunits are becoming increasingly well-defined, with much current research focusing on the impact of these proteins on the normal functioning or disease/stress states of cells. In this review, we provide an overview of the general functions of CPSF and its subunits, followed by a discussion of how they exert their functions in a surprisingly diverse variety of biological processes and cellular conditions. These include transcription termination, small RNA processing, and R-loop prevention/resolution, as well as more generally cancer, differentiation/development, and infection/immunity.

Circumvention of Topoisomerase IIα Intron 19 Intronic Polyadenylation in Acquired Etoposide-Resistant Human Leukemia K562 Cells

DNA topoisomerase IIα (TOP2α; 170 kDa, TOP2α/170) is an essential enzyme for proper chromosome dysjunction by producing transient DNA double-stranded breaks and is an important target for DNA damage-stabilizing anticancer agents, such as etoposide. Therapeutic effects of TOP2α poisons can be limited due to acquired drug resistance. We previously demonstrated decreased TOP2α/170 levels in an etoposide-resistant human leukemia K562 subline, designated K/VP.5, accompanied by increased expression of a C-terminal truncated TOP2α isoform (90 kDa; TOP2α/90), which heterodimerized with TOP2α/170 and was a determinant of resistance by exhibiting dominant-negative effects against etoposide activity. Based on 3'-rapid amplification of cDNA ends, we confirmed TOP2α/90 as the translation product of a TOP2α mRNA in which a cryptic polyadenylation site (PAS) harbored in intron 19 (I19) was used. In this report, we investigated whether the resultant intronic polyadenylation (IPA) would be attenuated by blocking or mutating the I19 PAS, thereby circumventing acquired drug resistance. An antisense morpholino oligonucleotide was used to hybridize/block the PAS in TOP2α pre-mRNA in K/VP.5 cells, resulting in decreased TOP2α/90 mRNA/protein levels in K/VP.5 cells and partially circumventing drug resistance. Subsequently, CRISPR/CRISPR-associated protein 9 with homology-directed repair was used to mutate the cryptic I19 PAS (AATAAA→ACCCAA) to prevent IPA. Gene-edited clones exhibited increased TOP2α/170 and decreased TOP2α/90 mRNA/protein and demonstrated restored sensitivity to etoposide and other TOP2α-targeted drugs. Together, results indicated that blocking/mutating a cryptic I19 PAS in K/VP.5 cells reduced IPA and restored sensitivity to TOP2α-targeting drugs. SIGNIFICANCE STATEMENT: The results presented in this study indicate that CRISPR/CRISPR-associated protein 9 gene editing of a cryptic polyadenylation site (PAS) within I19 of the TOP2α gene results in the reversal of acquired resistance to etoposide and other TOP2-targeted drugs. An antisense morpholino oligonucleotide targeting the PAS also partially circumvented resistance.

Chromatin regulates alternative polyadenylation via the RNA polymerase II elongation rate

The RNA polymerase II (Pol II) elongation rate influences poly(A) site selection, with slow and fast Pol II derivatives causing upstream and downstream shifts, respectively, in poly(A) site utilization. In yeast, depletion of either of the histone chaperones FACT or Spt6 causes an upstream shift of poly(A) site use that strongly resembles the poly(A) profiles of slow Pol II mutant strains. Like slow Pol II mutant strains, FACT- and Spt6-depleted cells exhibit Pol II processivity defects, indicating that both Spt6 and FACT stimulate the Pol II elongation rate. Poly(A) profiles of some genes show atypical downstream shifts; this subset of genes overlaps well for FACT- or Spt6-depleted strains but is different from the atypical genes in Pol II speed mutant strains. In contrast, depletion of histone H3 or H4 causes a downstream shift of poly(A) sites for most genes, indicating that nucleosomes inhibit the Pol II elongation rate in vivo. Thus, chromatin-based control of the Pol II elongation rate is a potential mechanism, distinct from direct effects on the cleavage/polyadenylation machinery, to regulate alternative polyadenylation in response to genetic or environmental changes.

Roles for epithelial integrin α3β1 in regulation of the microenvironment during normal and pathological tissue remodeling

Integrin receptors for the extracellular matrix activate intracellular signaling pathways that are critical for tissue development, homeostasis, and regeneration/repair, and their loss or dysregulation contributes to many developmental defects and tissue pathologies. This review will focus on tissue remodeling roles for integrin α3β1, a receptor for laminins found in the basement membranes (BMs) that underlie epithelial cell layers. As a paradigm, we will discuss literature that supports a role for α3β1 in promoting ability of epidermal keratinocytes to modify their tissue microenvironment during skin development, wound healing, or tumorigenesis. Preclinical and clinical studies have shown that this role depends largely on ability of α3β1 to govern the keratinocyte's repertoire of secreted proteins, or the "secretome," including 1) matrix proteins and proteases involved in matrix remodeling and 2) paracrine-acting growth factors/cytokines that stimulate other cells with important tissue remodeling functions (e.g., endothelial cells, fibroblasts, inflammatory cells). Moreover, α3β1 signaling controls gene expression that helps epithelial cells carry out these functions, including genes that encode secreted matrix proteins, proteases, growth factors, or cytokines. We will review what is known about α3β1-dependent gene regulation through both transcription and posttranscriptional mRNA stability. Regarding the latter, we will discuss examples of α3β1-dependent alternative splicing (AS) or alternative polyadenylation (APA) that prevents inclusion of cis-acting mRNA sequences that would otherwise target the transcript for degradation via nonsense-mediated decay or destabilizing AU-rich elements (AREs) in the 3'-untranslated region (3'-UTR). Finally, we will discuss prospects and anticipated challenges of exploiting α3β1 as a clinical target for the treatment of cancer or wound healing.

Star-PAP controls oncogene expression through primary miRNA 3'-end formation to regulate cellular proliferation and tumour formation

Star-PAP is a non-canonical poly(A) polymerase that is down regulated in breast cancer. While Star-PAP down regulation impairs target mRNA polyadenylation, paradoxically, we see up regulation of a large number of oncogenes on Star-PAP knockdown. Using two breast cancer cells (MCF7 with high Star-PAP, and MDA-MB-231 with negligible Star-PAP level), we discover that Star-PAP negatively regulates oncogene expression and subsequently cellular proliferation. This regulation is compromised with Star-PAP mutant of 3'-end processing function (serine 6 to alanine, S6A phospho-mutation). Concomitantly, xenograft mice model using MDA-MB-231 cells reveals a reduction in the tumour formation on ectopic Star-PAP expression that is ameliorated by S6A mutation. We find that Star-PAP control of target oncogene expression is independent of Star-PAP-mediated alternative polyadenylation or target mRNA 3'-end formation. We demonstrate that Star-PAP regulates target oncogenes through cellular miRNAs (miR-421, miR-335, miR-424, miR-543, miR-205, miR-34a, and miR-26a) that are down regulated in breast cancer. Analysis of various steps in miRNA biogenesis pathway reveals that Star-PAP regulates 3'-end formation and synthesis of primary miRNA (host) transcripts that is dependent on S6 phosphorylation thus controlling mature miRNA generation. Using mimics and inhibitors of two target miRNAs (miR-421 and miR-424) after Star-PAP depletion in MCF7 or ectopic expression in MDA-MB-231 cells, we demonstrate that Star-PAP controls oncogene expression and cellular proliferation through targeting miRNAs that regulates tumour formation. Our study establishes a novel mechanism of oncogene expression independent of alternative polyadenylation through Star-PAP-mediated miRNA host transcript polyadenylation that regulates breast cancer progression.

How is polyadenylation restricted to 3'-untranslated regions?

Polyadenylation occurs at numerous sites within 3'-untranslated regions (3'-UTRs) but rarely within coding regions. How does Pol II travel through long coding regions without generating poly(A) sites, yet then permits promiscuous polyadenylation once it reaches the 3'-UTR? The cleavage/polyadenylation (CpA) machinery preferentially associates with 3'-UTRs, but it is unknown how its recruitment is restricted to 3'-UTRs during Pol II elongation. Unlike coding regions, 3'-UTRs have long AT-rich stretches of DNA that may be important for restricting polyadenylation to 3'-UTRs. Recognition of the 3'-UTR could occur at the DNA (AT-rich), RNA (AU-rich), or RNA:DNA hybrid (rU:dA- and/or rA:dT-rich) level. Based on the nucleic acid critical for 3'-UTR recognition, there are three classes of models, not mutually exclusive, for how the CpA machinery is selectively recruited to 3'-UTRs, thereby restricting where polyadenylation occurs: (1) RNA-based models suggest that the CpA complex directly (or indirectly through one or more intermediary proteins) binds long AU-rich stretches that are exposed after Pol II passes through these regions. (2) DNA-based models suggest that the AT-rich sequence affects nucleosome depletion or the elongating Pol II machinery, resulting in dissociation of some elongation factors and subsequent recruitment of the CpA machinery. (3) RNA:DNA hybrid models suggest that preferential destabilization of the Pol II elongation complex at rU:dA- and/or rA:dT-rich duplexes bridging the nucleotide addition and RNA exit sites permits preferential association of the CpA machinery with 3'-UTRs. Experiments to provide evidence for one or more of these models are suggested.

Emerging roles of RNA-binding proteins in fatty liver disease

A rampant and urgent global health issue of the 21st century is the emergence and progression of fatty liver disease (FLD), including alcoholic fatty liver disease and the more heterogenous metabolism-associated (or non-alcoholic) fatty liver disease (MAFLD/NAFLD) phenotypes. These conditions manifest as disease spectra, progressing from benign hepatic steatosis to symptomatic steatohepatitis, cirrhosis, and, ultimately, hepatocellular carcinoma. With numerous intricately regulated molecular pathways implicated in its pathophysiology, recent data have emphasized the critical roles of RNA-binding proteins (RBPs) in the onset and development of FLD. They regulate gene transcription and post-transcriptional processes, including pre-mRNA splicing, capping, and polyadenylation, as well as mature mRNA transport, stability, and translation. RBP dysfunction at every point along the mRNA life cycle has been associated with altered lipid metabolism and cellular stress response, resulting in hepatic inflammation and fibrosis. Here, we discuss the current understanding of the role of RBPs in the post-transcriptional processes associated with FLD and highlight the possible and emerging therapeutic strategies leveraging RBP function for FLD treatment. This article is categorized under: RNA in Disease and Development > RNA in Disease.

Gene expression investigation of four key regulators of polyadenylation and alternative adenylation in the periphery of late-onset Alzheimer's disease patients

BACKGROUND

Alzheimer's disease (AD) is a genetic and sporadic neurodegenerative disease considered by an archetypal cognitive impairment and a decrease in less common cognitive impairment. Notably, the discovery of goals in this paradigm is still a challenge, and understanding basic mechanisms is an important step toward improving disease management. Polyadenylation (PA) and alternative polyadenylation (APA) are two of the most critical RNA processing stages in 3'UTRs that influence various AD-related genes.

METHODS

In this study, we assessed Cleavage and polyadenylation specificity factors 1 and 6 (CPSF1 and CPSF6), cleavage stimulation factor 1 (CSTF1), and WD Repeat Domain 33 (WDR33) genes expression in the periphery of 50 AD patients and 50 healthy individuals with age and gender-matched by quantitative real-time PCR.

RESULTS

Comparing AD patients with healthy people using expression analysis revealed a substantial increase in CSTF1 (posterior beta = 0.773, adjusted P-value = 0.042). Significant positive correlations were found between CSTF1 and CPSF1 (r = 0.365, P < 0.001), WDR33 (r = 0.506, P < 0.001), and CPSF6 (r = 0.446, P < 0.001) expression levels.

CONCLUSION

Although further research is required to determine their potential contribution to AD, our findings offer a fresh perspective on molecular regulatory pathways associated with AD pathogenic mechanisms associated with PA and APA.

Single-Molecule Poly(A) Tail Sequencing (SM-PATseq) Using the PacBio Platform

The polyadenylation of the 3' ends of messenger RNAs is an important regulator of stability and translation. We developed the single-molecule poly(A) tail sequencing method, SM-PATseq, to assay tail lengths of the whole transcriptome at nucleotide resolution using long-read sequencing. This method generates cDNA using an oligo-dT 3' splint adaptor ligation to prime first-strand cDNA synthesis, followed by random hexamer priming for second-strand synthesis. By directly sequencing the cDNA on long-read platforms, we can resolve tail lengths at nucleotide resolution, identify non-A bases within the tail, and quantify transcript abundance analogous to traditional RNAseq methods. Here, we discuss the method for generating, sequencing, and primary analysis of poly(A) tail data from total RNA using the Pacific Biosciences Sequel platform.

Impact of MICA 3'UTR allelic variability on miRNA binding prediction, a bioinformatic approach

MicroRNAs (miRNAs) are small non-coding RNAs that participate as powerful genetic regulators. MiRNAs can interfere with cellular processes by interacting with a broad spectrum of target genes under physiological and pathological states, including cancer development and progression. Major histocompatibility complex major histocompatibility complex class I-related chain A (MICA) belongs to a family of proteins that bind the natural-killer group 2, member D (NKG2D) receptor on Natural Killer cells and other cytotoxic lymphocytes. MICA plays a crucial role in the host's innate immune response to several disease settings, including cancer. MICA harbors various single nucleotide polymorphisms (SNPs) located in its 3'-untranslated region (3'UTR), a characteristic that increases the complexity of MICA regulation, favoring its post-transcriptional modulation by miRNAs under physiological and pathological conditions. Here, we conducted an in-depth analysis of MICA 3'UTR sequences according to each MICA allele described to date using NCBI database. We also systematically evaluated interactions between miRNAs and their putative targets on MICA 3'UTR containing SNPs using in silico analysis. Our in silico results showed that MICA SNPs rs9266829, rs 1880, and rs9266825, located in the target sequence of miRNAs hsa-miR-106a-5p, hsa-miR-17-5p, hsa-miR-20a-5p, hsa-miR-20b-5p, hsa-miR-93, hsa-miR-1207.5p, and hsa-miR-711 could modify the binding free energy between -8.62 and -18.14 kcal/mol, which may affect the regulation of MICA expression. We believe that our results may provide a starting point for further exploration of miRNA regulatory effects depending on MICA allelic variability; they may also be a guide to conduct miRNA in silico analysis for other highly polymorphic genes.

Members of the CUGBP Elav-like family of RNA-binding proteins are expressed in distinct populations of primary sensory neurons

Primary sensory dorsal root ganglia (DRG) neurons are diverse, with distinct populations that respond to specific stimuli. Previously, we observed that functionally distinct populations of DRG neurons express mRNA transcript variants with different 3' untranslated regions (3'UTRs). 3'UTRs harbor binding sites for interaction with RNA-binding proteins (RBPs) for transporting mRNAs to subcellular domains, modulating transcript stability, and regulating the rate of translation. In the current study, analysis of publicly available single-cell RNA-sequencing data generated from adult mice revealed that 17 3'UTR-binding RBPs were enriched in specific populations of DRG neurons. This included four members of the CUG triplet repeat (CUGBP) Elav-like family (CELF): CELF2 and CELF4 were enriched in peptidergic, CELF6 in both peptidergic and nonpeptidergic, and CELF3 in tyrosine hydroxylase-expressing neurons. Immunofluorescence studies confirmed that 60% of CELF4+ neurons are small-diameter C fibers and 33% medium-diameter myelinated (likely Aδ) fibers and showed that CELF4 is distributed to peripheral termini. Coexpression analyses using transcriptomic data and immunofluorescence revealed that CELF4 is enriched in nociceptive neurons that express GFRA3, CGRP, and the capsaicin receptor TRPV1. Reanalysis of published transcriptomic data from macaque DRG revealed a highly similar distribution of CELF members, and reanalysis of single-nucleus RNA-sequencing data derived from mouse and rat DRG after sciatic injury revealed differential expression of CELFs in specific populations of sensory neurons. We propose that CELF RBPs may regulate the fate of mRNAs in populations of nociceptors, and may play a role in pain and/or neuronal regeneration following nerve injury.

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.

The predicted RNA-binding protein regulome of axonal mRNAs

Neurons are morphologically complex cells that rely on the compartmentalization of protein expression to develop and maintain their cytoarchitecture. The targeting of RNA transcripts to axons is one of the mechanisms that allows rapid local translation of proteins in response to extracellular signals. 3' Untranslated regions (UTRs) of mRNA are noncoding sequences that play a critical role in determining transcript localization and translation by interacting with specific RNA-binding proteins (RBPs). However, how 3' UTRs contribute to mRNA metabolism and the nature of RBP complexes responsible for these functions remains elusive. We performed 3' end sequencing of RNA isolated from cell bodies and axons of sympathetic neurons exposed to either nerve growth factor (NGF) or neurotrophin 3 (NTF3, also known as NT-3). NGF and NTF3 are growth factors essential for sympathetic neuron development through distinct signaling mechanisms. Whereas NTF3 acts mostly locally, NGF signal is retrogradely transported from axons to cell bodies. We discovered that both NGF and NTF3 affect transcription and alternative polyadenylation in the nucleus and induce the localization of specific 3' UTR isoforms to axons, including short 3' UTR isoforms found exclusively in axons. The integration of our data with CLIP sequencing data supports a model whereby long 3' UTR isoforms associate with RBP complexes in the nucleus and, upon reaching the axons, are remodeled locally into shorter isoforms. Our findings shed new light into the complex relationship between nuclear polyadenylation, mRNA localization, and local 3' UTR remodeling in developing neurons.

Pro-inflammatory polarization and colorectal cancer modulate alternative and intronic polyadenylation in primary human macrophages

Introduction

Macrophages are essential cells of the immune system that alter their inflammatory profile depending on their microenvironment. Alternative polyadenylation in the 3'UTR (3'UTR-APA) and intronic polyadenylation (IPA) are mechanisms that modulate gene expression, particularly in cancer and activated immune cells. Yet, how polarization and colorectal cancer (CRC) cells affect 3'UTR-APA and IPA in primary human macrophages was unclear.

Methods

In this study, we isolated primary human monocytes from healthy donors, differentiated and polarized them into a pro-inflammatory state and performed indirect co-cultures with CRC cells. ChrRNA-Seq and 3'RNA-Seq was performed to quantify gene expression and characterize new 3'UTR-APA and IPA mRNA isoforms.

Results

Our results show that polarization of human macrophages from naïve to a pro-inflammatory state causes a marked increase of proximal polyA site selection in the 3'UTR and IPA events in genes relevant to macrophage functions. Additionally, we found a negative correlation between differential gene expression and IPA during pro-inflammatory polarization of primary human macrophages. As macrophages are abundant immune cells in the CRC microenvironment that either promote or abrogate cancer progression, we investigated how indirect exposure to CRC cells affects macrophage gene expression and 3'UTR-APA and IPA events. Co-culture with CRC cells alters the inflammatory phenotype of macrophages, increases the expression of pro-tumoral genes and induces 3'UTR-APA alterations. Notably, some of these gene expression differences were also found in tumor-associated macrophages of CRC patients, indicating that they are physiologically relevant. Upon macrophage pro-inflammatory polarization, SRSF12 is the pre-mRNA processing gene that is most upregulated. After SRSF12 knockdown in M1 macrophages there is a global downregulation of gene expression, in particular in genes involved in gene expression regulation and in immune responses.

Discussion

Our results reveal new 3'UTR-APA and IPA mRNA isoforms produced during pro-inflammatory polarization of primary human macrophages and CRC co-culture that may be used in the future as diagnostic or therapeutic tools. Furthermore, our results highlight a function for SRSF12 in pro-inflammatory macrophages, key cells in the tumor response.

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.

The peptide woods are lovely, dark and deep: Hunting for novel cancer antigens

Harnessing the patient's immune system to control a tumor is a proven avenue for cancer therapy. T cell therapies as well as therapeutic vaccines, which target specific antigens of interest, are being explored as treatments in conjunction with immune checkpoint blockade. For these therapies, selecting the best suited antigens is crucial. Most of the focus has thus far been on neoantigens that arise from tumor-specific somatic mutations. Although there is clear evidence that T-cell responses against mutated neoantigens are protective, the large majority of these mutations are not immunogenic. In addition, most somatic mutations are unique to each individual patient and their targeting requires the development of individualized approaches. Therefore, novel antigen types are needed to broaden the scope of such treatments. We review high throughput approaches for discovering novel tumor antigens and some of the key challenges associated with their detection, and discuss considerations when selecting tumor antigens to target in the clinic.

Identification and characterization of the Doublesex gene and its mRNA isoforms in the brine shrimp Artemia franciscana

Doublesex (DSX) proteins are members of the Doublesex/mab-3-related (DMRT) protein family and play crucial roles in sex determination and differentiation among the animal kingdom. In the present study, we identified two Doublesex (Dsx)-like mRNA isoforms in the brine shrimp Artemia franciscana (Kellogg 1906), which are generated by the combination of alternative promoters, alternative splicing and alternative polyadenylation. The two transcripts exhibited sex-biased enrichment, which we termed AfrDsxM and AfrDsxF. They share a common region which encodes an identical N-terminal DNA-binding (DM) domain. RT-qPCR analyses showed that AfrDsxM is dominantly expressed in male Artemia while AfrDsxF is specifically expressed in females. Expression levels of both isoforms increased along with the developmental stages of their respective sexes. RNA interference with dsRNA showed that the knockdown of AfrDsxM in male larvae led to the appearance of female traits including an ovary-like structure in the original male reproductive system and an elevated expression of vitellogenin. However, silencing of AfrDsxF induced no clear phenotypic change in female Artemia. These results indicated that the male AfrDSXM may act as inhibiting regulator upon the default female developmental mode in Artemia. Furthermore, electrophoretic mobility shift assay analyses revealed that the unique DM domain of AfrDSXs can specifically bind to promoter segments of potential downstream target genes like AfrVtg. These data show that AfrDSXs play crucial roles in regulating sexual development in Artemia, and further provide insight into the evolution of sex determination/differentiation in sexual organisms.

Transcriptomics-proteomics Integration reveals alternative polyadenylation driving inflammation-related protein translation in patients with diabetic nephropathy

BACKGROUND

Diabetic nephropathy (DN) is a complex disease involving the upregulation of many inflammation-related proteins. Alternative polyadenylation (APA), a crucial post-transcriptional regulatory mechanism, has been proven to play vital roles in many inflammatory diseases. However, it is largely unknown whether and how APA exerts function in DN.

METHODS

We performed transcriptomics and proteomics analysis of glomeruli samples isolated from 50 biopsy-proven DN patients and 25 control subjects. DaPars and QAPA algorithms were adopted to identify APA events from RNA-seq data. The qRT-PCR analysis was conducted to verify 3'UTR length alteration. Short and long 3'UTRs isoforms were also overexpressed in podocytes under hyperglycemia condition for examining protein expression.

RESULTS

We detected transcriptome-wide 3'UTR APA events in DN, and found that APA-mediated 3'UTR lengthening of genes (APA genes) increased their expression at protein but not mRNA level. Increased protein level of 3'UTR lengthening gene was validated in podocytes under hyperglycemia condition. Pathway enrichment analysis showed that APA genes were enriched in inflammation-related biological processes including endoplasmic reticulum stress pathways, NF-κB signaling and autophagy. Further bioinformatics analysis demonstrated that 3'UTR APA of genes probably altered the binding sites for RNA-binding proteins, thus enhancing protein translation.

CONCLUSION

This study revealed for the first time that 3'UTR lengthening of APA genes contributed to the progression of DN by elevating the translation of corresponding proteins, providing new insight and a rich resource for investigating DN mechanisms.

3'UTR Diversity: Expanding Repertoire of RNA Alterations in Human mRNAs

Genomic information stored in the DNA is transcribed to the mRNA and translated to proteins. The 3' untranslated regions (3'UTRs) of the mRNA serve pivotal roles in posttranscriptional gene expression, regulating mRNA stability, translation, and localization. Similar to DNA mutations producing aberrant proteins, RNA alterations expand the transcriptome landscape and change the cellular proteome. Recent global analyses reveal that many genes express various forms of altered RNAs, including 3'UTR length variants. Alternative polyadenylation and alternative splicing are involved in diversifying 3'UTRs, which could act as a hidden layer of eukaryotic gene expression control. In this review, we summarize the functions and regulations of 3'UTRs and elaborate on the generation and functional consequences of 3'UTR diversity. Given that dynamic 3'UTR length control contributes to phenotypic complexity, dysregulated 3'UTR diversity might be relevant to disease development, including cancers. Thus, 3'UTR diversity in cancer could open exciting new research areas and provide avenues for novel cancer theragnostics.

PABPC1--mRNA stability, protein translation and tumorigenesis

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

Nucleotide-level linkage of transcriptional elongation and polyadenylation

Alternative polyadenylation yields many mRNA isoforms whose 3' termini occur disproportionately in clusters within 3' untranslated regions. Previously, we showed that profiles of poly(A) site usage are regulated by the rate of transcriptional elongation by RNA polymerase (Pol) II (Geisberg et al., 2020). Pol II derivatives with slow elongation rates confer an upstream-shifted poly(A) profile, whereas fast Pol II strains confer a downstream-shifted poly(A) profile. Within yeast isoform clusters, these shifts occur steadily from one isoform to the next across nucleotide distances. In contrast, the shift between clusters - from the last isoform of one cluster to the first isoform of the next - is much less pronounced, even over large distances. GC content in a region 13-30 nt downstream from isoform clusters correlates with their sensitivity to Pol II elongation rate. In human cells, the upstream shift caused by a slow Pol II mutant also occurs continuously at single nucleotide resolution within clusters but not between them. Pol II occupancy increases just downstream of poly(A) sites, suggesting a linkage between reduced elongation rate and cluster formation. These observations suggest that (1) Pol II elongation speed affects the nucleotide-level dwell time allowing polyadenylation to occur, (2) poly(A) site clusters are linked to the local elongation rate, and hence do not arise simply by intrinsically imprecise cleavage and polyadenylation of the RNA substrate, (3) DNA sequence elements can affect Pol II elongation and poly(A) profiles, and (4) the cleavage/polyadenylation and Pol II elongation complexes are spatially, and perhaps physically, coupled so that polyadenylation occurs rapidly upon emergence of the nascent RNA from the Pol II elongation complex.

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

Background

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

Methods

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

Result

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

Conclusion

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

Supplementary Information

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

A Survey on Methods for Predicting Polyadenylation Sites from DNA Sequences, Bulk RNA-seq, and Single-cell RNA-seq

Alternative polyadenylation (APA) plays important roles in modulating mRNA stability, translation, and subcellular localization, and contributes extensively to shaping eukaryotic transcriptome complexity and proteome diversity. Identification of poly(A) sites (pAs) on a genome-wide scale is a critical step toward understanding the underlying mechanism of APA-mediated gene regulation. A number of established computational tools have been proposed to predict pAs from diverse genomic data. Here we provided an exhaustive overview of computational approaches for predicting pAs from DNA sequences, bulk RNA sequencing (RNA-seq) data, and single-cell RNA sequencing (scRNA-seq) data. Particularly, we examined several representative tools using bulk RNA-seq and scRNA-seq data from peripheral blood mononuclear cells and put forward operable suggestions on how to assess the reliability of pAs predicted by different tools. We also proposed practical guidelines on choosing appropriate methods applicable to diverse scenarios. Moreover, we discussed in depth the challenges in improving the performance of pA prediction and benchmarking different methods. Additionally, we highlighted outstanding challenges and opportunities using new machine learning and integrative multi-omics techniques, and provided our perspective on how computational methodologies might evolve in the future for non-3' untranslated region, tissue-specific, cross-species, and single-cell pA prediction.

3' Untranslated Regions Are Modular Entities That Determine Polyadenylation Profiles

The 3' ends of eukaryotic mRNAs are generated by cleavage of nascent transcripts followed by polyadenylation, which occurs at numerous sites within 3' untranslated regions (3' UTRs) but rarely within coding regions. An individual gene can yield many 3'-mRNA isoforms with distinct half-lives. We dissect the relative contributions of protein-coding sequences (open reading frames [ORFs]) and 3' UTRs to polyadenylation profiles in yeast. ORF-deleted derivatives often display strongly decreased mRNA levels, indicating that ORFs contribute to overall mRNA stability. Poly(A) profiles, and hence relative isoform half-lives, of most (9 of 10) ORF-deleted derivatives are very similar to their wild-type counterparts. Similarly, in-frame insertion of a large protein-coding fragment between the ORF and 3' UTR has minimal effect on the poly(A) profile in all 15 cases tested. Last, reciprocal ORF/3'-UTR chimeric genes indicate that the poly(A) profile is determined by the 3' UTR. Thus, 3' UTRs are self-contained modular entities sufficient to determine poly(A) profiles and relative 3'-isoform half-lives. In the one atypical instance, ORF deletion causes an upstream shift of poly(A) sites, likely because juxtaposition of an unusually high AT-rich stretch directs polyadenylation closely downstream. This suggests that long AT-rich stretches, which are not encountered until after coding regions, are important for restricting polyadenylation to 3' UTRs.

Analysis of Polyadenylation Signal Usage with Full-Length Transcriptome in Spodoptera frugiperda (Lepidoptera: Noctuidae)

During the messenger RNA (mRNA) maturation process, RNA polyadenylation is a key step, and is coupled to the termination of transcription. Various cis-acting elements near the cleavage site and their binding factors would affect the process of polyadenylation, and AAUAAA, a highly conserved hexamer, was the most important polyadenylation signal (PAS). PAS usage is one of the critical modification determinants targeted at mRNA post-transcription. The full-length transcriptome has recently generated a massive amount of sequencing data, revealing poly(A) variation and alternative polyadenylation (APA) in Spodoptera frugiperda. We identified 50,616 polyadenylation signals in Spodoptera frugiperda via analysis of full-length transcriptome combined with expression Sequence Tags Technology (EST). The polyadenylation signal usage in Spodoptera frugiperda is conserved, and it is similar to that of flies and other animals. AAUAAA and AUUAAA are the most highly conserved polyadenylation signals of all polyadenylation signals we identified. Additionally, we found the U/GU-rich downstream sequence element (DSE) in the cleavage site. These results demonstrate that APA in Spodoptera frugiperda plays a significant role in root growth and development. This is the first polyadenylation signal usage analysis in agricultural pests, which can deepen our understanding of Spodoptera frugiperda and provide a theoretical basis for pest control.

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

Simple Summary

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

Abstract

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

It's a DoG-eat-DoG world-altered transcriptional mechanisms drive downstream-of-gene (DoG) transcript production

The past decade has revolutionized our understanding of regulatory noncoding RNAs (ncRNAs). Among the most recently identified ncRNAs are downstream-of-gene (DoG)-containing transcripts that are produced by widespread transcriptional readthrough. The discovery of DoGs has set the stage for future studies to address many unanswered questions regarding the mechanisms that promote readthrough transcription, RNA processing, and the cellular functions of the unique transcripts. In this review, we summarize current findings regarding the biogenesis, function, and mechanisms regulating this exciting new class of RNA molecules.

Full-length transcriptome reconstruction reveals genetic differences in hybrids of Oryza sativa and Oryza punctata with different ploidy and genome compositions

BACKGROUND

Allopolyploid breeding is an efficient technique for improving the low seed setting rate of autotetraploids in plant breeding and one of the most promising breeding methods. However, there have been few comprehensive studies of the posttranscriptional mechanism in allopolyploids.

RESULTS

By crossing cultivated rice (Oryza sativa, genome AA) with wild rice (Oryza punctata, genome BB), we created hybrid rice lines with different ploidy and genome compositions [diploid hybrid F01 (AB), allotetraploid hybrid F02 (AABB) and F03 (AAAB)]. The genetic differences of the hybrids and the mechanism of allopolyploid breeding dominance were revealed through morphological and cytological observations and single molecule real-time sequencing techniques. The tissues and organs of allotetraploid hybrid F02 exhibited "gigantism" and the highest levels of fertility. The numbers of non-redundant transcripts, gene loci and new isoforms in the polyploid rice lines were higher and the isoform lengths greater than those of the diploid line. Moreover, alternative splicing (AS) events occurred twice as often in the polyploid rice lines than the diploid line. During these events, intron retention dominated. Furthermore, a large number of new genes and isoforms specific to the lines of different ploidy were discovered.

CONCLUSIONS

The results indicated that alternative polyadenylation (APA) and AS events contributed to the complexity and superiority of polyploids in the activity of translation regulators, nucleic acid binding transcription factor activities and the regulation of molecular function. Therefore, these APA and AS events in allopolyploid rice were found to play a role in regulation. Our study provides new germplasm for polyploid rice breeding and reveals complex regulatory mechanisms that may be related to heterosis and fertility.

Translational control of E2f1 regulates the Drosophila cell cycle

E2F transcription factors are master regulators of the eukaryotic cell cycle. In Drosophila, the sole activating E2F, E2F1, is both required for and sufficient to promote G1→S progression. E2F1 activity is regulated both by binding to RB Family repressors and by posttranscriptional control of E2F1 protein levels by the EGFR and TOR signaling pathways. Here, we investigate cis-regulatory elements in the E2f1 messenger RNA (mRNA) that enable E2f1 translation to respond to these signals and promote mitotic proliferation of wing imaginal disc and intestinal stem cells. We show that small upstream open reading frames (uORFs) in the 5' untranslated region (UTR) of the E2f1 mRNA limit its translation, impacting rates of cell proliferation. E2f1 transgenes lacking these 5'UTR uORFs caused TOR-independent expression and excess cell proliferation, suggesting that TOR activity can bypass uORF-mediated translational repression. EGFR signaling also enhanced translation but through a mechanism less dependent on 5'UTR uORFs. Further, we mapped a region in the E2f1 mRNA that contains a translational enhancer, which may also be targeted by TOR signaling. This study reveals translational control mechanisms through which growth signaling regulates cell cycle progression.

Analysis Polyadenylation Signal Usage in Sus scrofa

RNA polyadenylation is an important step in the messenger RNA (mRNA) maturation process, and the first step is recognizing the polyadenylation signal (PAS). The PAS type and distribution is a key determinant of post-transcriptional mRNA modification and gene expression. However, little is known about PAS usage and alternative polyadenylation (APA) regulation in livestock species. Recently, sequencing technology has enabled the generation of a large amount of sequencing data revealing variation in poly(A) signals and APA regulation in Sus scrofa. We identified 62,491 polyadenylation signals in Sus scrofa using expressed sequence tag (EST) sequences combined with RNA-seq analysis. The composition and usage frequency of polyadenylation signal in Sus scrofa is similar with that of human and mouse. The most highly conserved polyadenylation signals are AAUAAA and AUUAAA, used for over 63.35% of genes. In addition, we also analyzed the U/GU-rich downstream sequence (DSE) element, located downstream of the cleavage site. Our results indicate that APA regulation was widely occurred in Sus scrofa, as in other organisms. Our result was useful for the accurate annotation of RNA 3' ends in Sus scrofa and the analysis of polyadenylation signal usage in Sus scrofa would give the new insights into the mechanisms of transcriptional regulation.

Tau Modulates mRNA Transcription, Alternative Polyadenylation Profiles of hnRNPs, Chromatin Remodeling and Spliceosome Complexes

Tau protein is a known contributor in several neurodegenerative diseases, including Alzheimer’s disease (AD) and frontotemporal dementia (FTD). It is well-established that tau forms pathological aggregates and fibrils in these diseases. Tau has been observed within the nuclei of neurons, but there is a gap in understanding regarding the mechanism by which tau modulates transcription. We are interested in the P301L mutation of tau, which has been associated with FTD and increased tau aggregation. Our study utilized tau-inducible HEK (iHEK) cells to reveal that WT and P301L tau distinctively alter the transcription and alternative polyadenylation (APA) profiles of numerous nuclear precursors mRNAs, which then translate to form proteins involved in chromatin remodeling and splicing. We isolated total mRNA before and after over-expressing tau and then performed Poly(A)-ClickSeq (PAC-Seq) to characterize mRNA expression and APA profiles. We characterized changes in Gene Ontology (GO) pathways using EnrichR and Gene Set Enrichment Analysis (GSEA). We observed that P301L tau up-regulates genes associated with reactive oxygen species responsiveness as well as genes involved in dendrite, microtubule, and nuclear body/speckle formation. The number of genes regulated by WT tau is greater than the mutant form, which indicates that the P301L mutation causes loss-of-function at the transcriptional level. WT tau up-regulates genes contributing to cytoskeleton-dependent intracellular transport, microglial activation, microtubule and nuclear chromatin organization, formation of nuclear bodies and speckles. Interestingly, both WT and P301L tau commonly down-regulate genes responsible for ubiquitin-proteosome system. In addition, WT tau significantly down-regulates several genes implicated in chromatin remodeling and nucleosome organization. Although there are limitations inherent to the model systems used, this study will improve understanding regarding the nuclear impact of tau at the transcriptional and post-transcriptional level. This study also illustrates the potential impact of P301L tau on the human brain genome during early phases of pathogenesis.

Integrator enforces the fidelity of transcriptional termination at protein-coding genes

Integrator regulates the 3′-end processing and termination of multiple classes of noncoding RNAs. Depletion of INTS11, the catalytic subunit of Integrator, or ectopic expression of its catalytic dead enzyme impairs the 3′-end processing and termination of a set of protein-coding transcripts termed Integrator-regulated termination (IRT) genes. This defect is manifested by increased RNA polymerase II (RNAPII) readthrough and occupancy of serine-2 phosphorylated RNAPII, de novo trimethylation of lysine-36 on histone H3, and a compensatory elevation of the cleavage and polyadenylation (CPA) complex beyond the canonical polyadenylation sites. 3′ RNA sequencing reveals that proximal polyadenylation site usage relies on the endonuclease activity of INTS11. The DNA sequence encompassing the transcription end sites of IRT genes features downstream polyadenylation motifs and an enrichment of GC content that permits the formation of secondary structures within the 3′UTR. Together, this study identifies a subset of protein-coding transcripts whose 3′ end processing requires the Integrator complex.

RNA Modifications and RNA Metabolism in Neurological Disease Pathogenesis

The intrinsic cellular heterogeneity and molecular complexity of the mammalian nervous system relies substantially on the dynamic nature and spatiotemporal patterning of gene expression. These features of gene expression are achieved in part through mechanisms involving various epigenetic processes such as DNA methylation, post-translational histone modifications, and non-coding RNA activity, amongst others. In concert, another regulatory layer by which RNA bases and sugar residues are chemically modified enhances neuronal transcriptome complexity. Similar RNA modifications in other systems collectively constitute the cellular epitranscriptome that integrates and impacts various physiological processes. The epitranscriptome is dynamic and is reshaped constantly to regulate vital processes such as development, differentiation and stress responses. Perturbations of the epitranscriptome can lead to various pathogenic conditions, including cancer, cardiovascular abnormalities and neurological diseases. Recent advances in next-generation sequencing technologies have enabled us to identify and locate modified bases/sugars on different RNA species. These RNA modifications modulate the stability, transport and, most importantly, translation of RNA. In this review, we discuss the formation and functions of some frequently observed RNA modifications—including methylations of adenine and cytosine bases, and isomerization of uridine to pseudouridine—at various layers of RNA metabolism, together with their contributions to abnormal physiological conditions that can lead to various neurodevelopmental and neurological disorders.

miRNA-Based Regulation of Alternative RNA Splicing in Metazoans

Alternative RNA splicing is an important regulatory process used by genes to increase their diversity. This process is mainly executed by specific classes of RNA binding proteins that act in a dosage-dependent manner to include or exclude selected exons in the final transcripts. While these processes are tightly regulated in cells and tissues, little is known on how the dosage of these factors is achieved and maintained. Several recent studies have suggested that alternative RNA splicing may be in part modulated by microRNAs (miRNAs), which are short, non-coding RNAs (~22 nt in length) that inhibit translation of specific mRNA transcripts. As evidenced in tissues and in diseases, such as cancer and neurological disorders, the dysregulation of miRNA pathways disrupts downstream alternative RNA splicing events by altering the dosage of splicing factors involved in RNA splicing. This attractive model suggests that miRNAs can not only influence the dosage of gene expression at the post-transcriptional level but also indirectly interfere in pre-mRNA splicing at the co-transcriptional level. The purpose of this review is to compile and analyze recent studies on miRNAs modulating alternative RNA splicing factors, and how these events contribute to transcript rearrangements in tissue development and disease.

Alternative polyadenylation: An enigma of transcript length variation in health and disease

Alternative polyadenylation (APA) is a molecular mechanism during a pre-mRNA processing that involves usage of more than one polyadenylation site (PA-site) generating transcripts of varying length from a single gene. The location of a PA-site affects transcript length and coding potential of an mRNA contributing to both mRNA and protein diversification. This variation in the transcript length affects mRNA stability and translation, mRNA subcellular and tissue localization, and protein function. APA is now considered as an important regulatory mechanism in the pathophysiology of human diseases. An important consequence of the changes in the length of 3'-untranslated region (UTR) from disease-induced APA is altered protein expression. Yet, the relationship between 3'-UTR length and protein expression remains a paradox in a majority of diseases. Here, we review occurrence of APA, mechanism of PA-site selection, and consequences of transcript length variation in different diseases. Emerging evidence reveals coordinated involvement of core RNA processing factors including poly(A) polymerases in the PA-site selection in diseases-associated APAs. Targeting such APA regulators will be therapeutically significant in combating drug resistance in cancer and other complex diseases. This article is categorized under: RNA Processing > 3' End Processing RNA in Disease and Development > RNA in Disease Translation > Regulation.

scAPAdb: a comprehensive database of alternative polyadenylation at single-cell resolution

Alternative polyadenylation (APA) is a widespread regulatory mechanism of transcript diversification in eukaryotes, which is increasingly recognized as an important layer for eukaryotic gene expression. Recent studies based on single-cell RNA-seq (scRNA-seq) have revealed cell-to-cell heterogeneity in APA usage and APA dynamics across different cell types in various tissues, biological processes and diseases. However, currently available APA databases were all collected from bulk 3′-seq and/or RNA-seq data, and no existing database has provided APA information at single-cell resolution. Here, we present a user-friendly database called scAPAdb (http://www.bmibig.cn/scAPAdb), which provides a comprehensive and manually curated atlas of poly(A) sites, APA events and poly(A) signals at the single-cell level. Currently, scAPAdb collects APA information from > 360 scRNA-seq experiments, covering six species including human, mouse and several other plant species. scAPAdb also provides batch download of data, and users can query the database through a variety of keywords such as gene identifier, gene function and accession number. scAPAdb would be a valuable and extendable resource for the study of cell-to-cell heterogeneity in APA isoform usages and APA-mediated gene regulation at the single-cell level under diverse cell types, tissues and species.

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

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

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

Single-molecule polyadenylated tail sequencing (SM-PAT-Seq) to measure polyA tail lengths transcriptome-wide

Polyadenylation of the 3' end of mRNAs is an important mechanism for regulating their stability and translation. We developed a nucleotide-resolution, transcriptome-wide, single-molecule SM-PAT-Seq method to accurately measure the polyA tail lengths of individual transcripts using long-read sequencing. The method generates cDNA using a double stranded splint adaptor targeting the far 3' end of the polyA tail for first strand synthesis along with random hexamers for second strand synthesis. This straight-forward method yields accurate polyA tail sequence lengths, can identify non-A residues in those tails, and quantitate transcript abundance.

Known sequence features can explain half of all human gene ends

Cleavage and polyadenylation (CPA) sites define eukaryotic gene ends. CPA sites are associated with five key sequence recognition elements: the upstream UGUA, the polyadenylation signal (PAS), and U-rich sequences; the CA/UA dinucleotide where cleavage occurs; and GU-rich downstream elements (DSEs). Currently, it is not clear whether these sequences are sufficient to delineate CPA sites. Additionally, numerous other sequences and factors have been described, often in the context of promoting alternative CPA sites and preventing cryptic CPA site usage. Here, we dissect the contributions of individual sequence features to CPA using standard discriminative models. We show that models comprised only of the five primary CPA sequence features give highest probability scores to constitutive CPA sites at the ends of coding genes, relative to the entire pre-mRNA sequence, for 41% of all human genes. U1-hybridizing sequences provide a small boost in performance. The addition of all known RBP RNA binding motifs to the model, however, increases this figure to 49%, and suggests an involvement of both known and suspected CPA regulators as well as potential new factors in delineating constitutive CPA sites. To our knowledge, this high effectiveness of established features to predict human gene ends has not previously been documented.

Multiple ways to a dead end: diverse mechanisms by which ALS mutant genes induce cell death

Amyotrophic Lateral Sclerosis (ALS) is a deadly neuromuscular disorder caused by progressive motor neuron loss in the brain and spinal cord. Over the past decades, a number of genetic mutations have been identified that cause or are associated with ALS disease progression. Numerous genes harbor ALS mutations, and they encode proteins displaying a wide range of physiological functions, with limited overlap. Despite the divergent functions, mutations in these genes typically trigger protein aggregation, which can confer gain- and/or loss-of-function to a number of essential cellular processes. Nuclear processes such as mRNA splicing and the response to DNA damage are significantly affected in ALS patients. Cytoplasmic organelles such as mitochondria are damaged by ALS mutant proteins. Processes that maintain cellular homeostasis such as autophagy, nonsense-mediated mRNA decay and nucleocytoplasmic transport, are also impaired by ALS mutations. Here, we review the multiple mechanisms by which mutations in major ALS-associated genes, such as TARDBP, C9ORF72 and FUS, lead to impairment of essential cellular processes.

Cytoplasmic cleavage of IMPA1 3' UTR is necessary for maintaining axon integrity

The 3' untranslated regions (3' UTRs) of messenger RNAs (mRNAs) are non-coding sequences involved in many aspects of mRNA metabolism, including intracellular localization and translation. Incorrect processing and delivery of mRNA cause severe developmental defects and have been implicated in many neurological disorders. Here, we use deep sequencing to show that in sympathetic neuron axons, the 3' UTRs of many transcripts undergo cleavage, generating isoforms that express the coding sequence with a short 3' UTR and stable 3' UTR-derived fragments of unknown function. Cleavage of the long 3' UTR of Inositol Monophosphatase 1 (IMPA1) mediated by a protein complex containing the endonuclease argonaute 2 (Ago2) generates a translatable isoform that is necessary for maintaining the integrity of sympathetic neuron axons. Thus, our study provides a mechanism of mRNA metabolism that simultaneously regulates local protein synthesis and generates an additional class of 3' UTR-derived RNAs.

Cross-species RNA-seq for deciphering host-microbe interactions

The human body is constantly exposed to microorganisms, which entails manifold interactions between human cells and diverse commensal or pathogenic bacteria. The cellular states of the interacting cells are decisive for the outcome of these encounters such as whether bacterial virulence programmes and host defence or tolerance mechanisms are induced. This Review summarizes how next-generation RNA sequencing (RNA-seq) has become a primary technology to study host-microbe interactions with high resolution, improving our understanding of the physiological consequences and the mechanisms at play. We illustrate how the discriminatory power and sensitivity of RNA-seq helps to dissect increasingly complex cellular interactions in time and space down to the single-cell level. We also outline how future transcriptomics may answer currently open questions in host-microbe interactions and inform treatment schemes for microbial disorders.

Epidermal progenitors suppress GRHL3-mediated differentiation through intronic polyadenylation promoted by CPSF-HNRNPA3 collaboration

In self-renewing somatic tissue such as skin epidermis, terminal differentiation genes must be suppressed in progenitors to sustain regenerative capacity. Here we show that hundreds of intronic polyadenylation (IpA) sites are differentially used during keratinocyte differentiation, which is accompanied by downregulation of the Cleavage and Polyadenylation Specificity Factor (CPSF) complex. Sustained CPSF expression in undifferentiated keratinocytes requires the contribution from the transcription factor MYC. In keratinocytes cultured in undifferentiation condition, CSPF knockdown induces premature differentiation and partially affects dynamically used IpA sites. These sites include an IpA site located in the first intron of the differentiation activator GRHL3. CRISPR knockout of GRHL3 IpA increased full-length GRHL3 mRNA expression. Using a targeted genetic screen, we identify that HNRNPA3 interacts with CPSF and enhances GRHL3 IpA. Our data suggest a model where the interaction between CPSF and RNA-binding proteins, such as HNRNPA3, promotes site-specific IpA and suppresses premature differentiation in progenitors.

Regulation of Kv11.1 Isoform Expression by Polyadenylate Binding Protein Nuclear 1

The Kv11.1 voltage-gated potassium channel, encoded by the KCNH2 gene, conducts the rapidly activating delayed rectifier current in the heart. KCNH2 pre-mRNA undergoes alternative polyadenylation to generate two C-terminal Kv11.1 isoforms in the heart. Utilization of a poly(A) signal in exon 15 produces the full-length, functional Kv11.1a isoform, while intron 9 polyadenylation generates the C-terminally truncated, nonfunctional Kv11.1a-USO isoform. The relative expression of Kv11.1a and Kv11.1a-USO isoforms plays an important role in the regulation of Kv11.1 channel function. In this study, we tested the hypothesis that the RNA polyadenylate binding protein nuclear 1 (PABPN1) interacts with a unique 22 nt adenosine stretch adjacent to the intron 9 poly(A) signal and regulates KCNH2 pre-mRNA alternative polyadenylation and the relative expression of Kv11.1a C-terminal isoforms. We showed that PABPN1 inhibited intron 9 poly(A) activity using luciferase reporter assays, tandem poly(A) reporter assays, and RNA pulldown assays. We also showed that PABPN1 increased the relative expression level of the functional Kv11.1a isoform using RNase protection assays, immunoblot analyses, and patch clamp recordings. Our present findings suggest a novel role for the RNA-binding protein PABPN1 in the regulation of functional and nonfunctional Kv11.1 isoform expression.

Advances in the Bioinformatics Knowledge of mRNA Polyadenylation in Baculovirus Genes

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

Alternative Polyadenylation: a new frontier in post transcriptional regulation

Polyadenylation of pre-messenger RNA (pre-mRNA) specific sites and termination of their downstream transcriptions are signaled by unique sequence motif structures such as AAUAAA and its auxiliary elements. Alternative polyadenylation (APA) is an important post-transcriptional regulatory mechanism that processes RNA products depending on its 3'-untranslated region (3'-UTR) specific sequence signal. APA processing can generate several mRNA isoforms from a single gene, which may have different biological functions on their target gene. As a result, cellular genomic stability, proliferation capability, and transformation feasibility could all be affected. Furthermore, APA modulation regulates disease initiation and progression. APA status could potentially act as a biomarker for disease diagnosis, severity stratification, and prognosis forecast. While the advance of modern throughout technologies, such as next generation-sequencing (NGS) and single-cell sequencing techniques, have enriched our knowledge about APA, much of APA biological process is unknown and pending for further investigation. Herein, we review the current knowledge on APA and how its regulatory complex factors (CFI/IIm, CPSF, CSTF, and RBPs) work together to determine RNA splicing location, cell cycle velocity, microRNA processing, and oncogenesis regulation. We also discuss various APA experiment strategies and the future direction of APA research.

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

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

The transcriptional elongation rate regulates alternative polyadenylation in yeast

Yeast cells undergoing the diauxic response show a striking upstream shift in poly(A) site utilization, with increased use of ORF-proximal poly(A) sites resulting in shorter 3' mRNA isoforms for most genes. This altered poly(A) pattern is extremely similar to that observed in cells containing Pol II derivatives with slow elongation rates. Conversely, cells containing derivatives with fast elongation rates show a subtle downstream shift in poly(A) sites. Polyadenylation patterns of many genes are sensitive to both fast and slow elongation rates, and a global shift of poly(A) utilization is strongly linked to increased purine content of sequences flanking poly(A) sites. Pol II processivity is impaired in diauxic cells, but strains with reduced processivity and normal Pol II elongation rates have normal polyadenylation profiles. Thus, Pol II elongation speed is important for poly(A) site selection and for regulating poly(A) patterns in response to environmental conditions.