Alternative Polyadenylation: Methods, Findings, and Impacts

Genetic analysis of PALB2 gene WD40 domain in canine mammary tumour patients

BACKGROUND

DNA repair mechanisms are essential for tumorigenesis and disruption of HR mechanism is an important predisposing factor of human breast cancers (BC). PALB2 is an important part of the HR. There are similarities between canine mammary tumours (CMT) and BCs. As its human counterpart, PALB2 mutations could be a predisposing factor of CMT.

OBJECTIVES

In this study, we aimed to investigate the impacts of PALB2 variants on tumorigenesis and canine mammary tumor (CMT) malignancy.

METHODS

We performed Sanger sequencing to detect germline mutations in the WD40 domain of the canine PALB2 gene in CMT patients. We conducted in silico analysis to investigate the variants, and compared the germline PALB2 mutations in humans that cause breast cancer (BC) with the variants detected in dogs with CMT.

RESULTS

We identified an intronic (c.3096+8C>G) variant, two exonic (p.A1050V and p.R1354R) variants, and a 3' UTR variant (c.4071T>C). Of these, p.R1354R and c.4071T>C novel variants were identified for the first time in this study. We found that the p.A1050V mutation had a significant effect. However, we could not determine sufficient similarity due to the differences in nucleotide/amino acid sequences between two species. Nonetheless, possible variants of human sequences in the exact location as their dog counterparts are associated with several cancer types, implying that the variants could be crucial for tumorigenesis in dogs. Our results did not show any effect of the variants on tumor malignancy.

CONCLUSIONS

The current project is the first study investigating the relationship between the PALB2 gene WD40 domain and CMTs. Our findings will contribute to a better understanding of the pathogenic mechanism of the PALB2 gene in CMTs. In humans, variant positions in canines have been linked to cancer-related phenotypes such as familial BC, endometrial tumor, and hereditary cancer predisposition syndrome. The results of bioinformatics analyses should be investigated through functional tests or case-control studies.

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.

Identification of sleep and circadian alternative polyadenylation sites associated with APA-linked human brain disorders

Sleep and circadian rhythm disruptions are comorbid features of many pathologies and can negatively influence numerous health conditions, including degenerative diseases, metabolic illnesses, cancer, and various neurological disorders. Genetic association studies linking sleep and circadian disturbances with disease susceptibility have mainly focused on changes in gene expression due to mutations, such as single-nucleotide polymorphisms. Thus, associations between sleep and/or circadian rhythm and alternative polyadenylation (APA), particularly in the context of other health challenges, are largely undescribed. APA is a process that generates various transcript isoforms from the same gene, resulting in effects on mRNA translation, stability, localization, and subsequent function. Here, we have identified unique APAs in rat brain that exhibit time-of-day-dependent oscillations in expression as well as APAs that are altered by sleep deprivation and the subsequent recovery period. Genes affected by APA usage include Mapt/Tau, Ntrk2, Homer1A, Sin3band Sorl. Sorl1 has two APAs which cycle with a 24 h period, one additional APA cycles with a 12 h period and one more that is reduced during recovery sleep. Finally, we compared sleep- or circadian-associated APAs with recently described APA-linked brain disorder susceptibility genes and found 46 genes in common.

Alternative polyadenylation reprogramming of MORC2 induced by NUDT21 loss promotes KIRC carcinogenesis

Alternative polyadenylation (APA), a posttranscriptional mechanism of gene expression via determination of 3'UTR length, has an emerging role in carcinogenesis. Although abundant APA reprogramming is found in kidney renal clear cell carcinoma (KIRC), which is one of the major malignancies, whether APA functions in KIRC remains unknown. Herein, we found that chromatin modifier MORC2 gained oncogenic potential in KIRC among the genes with APA reprogramming, and moreover, its oncogenic potential was enhanced by 3'UTR shortening through stabilization of MORC2 mRNA. MORC2 was found to function in KIRC by downregulating tumor suppressor DAPK1 via DNA methylation. Mechanistically, MORC2 recruited DNMT3A to facilitate hypermethylation of the DAPK1 promoter, which was strengthened by 3'UTR shortening of MORC2. Furthermore, loss of APA regulator NUDT21, which was induced by DNMT3B-mediated promoter methylation, was identified as responsible for 3'UTR shortening of MORC2 in KIRC. Additionally, NUDT21 was confirmed to act as a tumor suppressor mainly depending on downregulation of MORC2. Finally, we designed an antisense oligonucleotide (ASO) to enhance NUDT21 expression and validated its antitumor effect in vivo and in vitro. This study uncovers the DNMT3B/NUDT21/APA/MORC2/DAPK1 regulatory axis in KIRC, disclosing the role of APA in KIRC and the crosstalk between DNA methylation and APA.

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.

Oncogenic SRSF3 in health and diseases

Serine/arginine rich splicing factor 3 (SRSF3) is an important multi-functional splicing factor, and has attracted increasing attentions in the past thirty years. The importance of SRSF3 is evidenced by its impressively conserved protein sequences in all animals and alternative exon 4 which represents an autoregulatory mechanism to maintain its proper cellular expression level. New functions of SRSF3 have been continuously discovered recently, especially its oncogenic function. SRSF3 plays essential roles in many cellular processes by regulating almost all aspects of RNA biogenesis and processing of many target genes, and thus, contributes to tumorigenesis when overexpressed or disregulated. This review updates and highlights the gene, mRNA, and protein structure of SRSF3, the regulatory mechanisms of SRSF3 expression, and the characteristics of SRSF3 targets and binding sequences that contribute to SRSF3's diverse molecular and cellular functions in tumorigenesis and human diseases.

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.

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.

Coupling epigenetics and RNA polyadenylation: missing links

Precise regulation of gene expression is crucial for plant survival. As a cotranscriptional regulatory mechanism, pre-mRNA polyadenylation is essential for fine-tuning gene expression. Polyadenylation can be alternatively projected at various sites of a transcript, which contributes to transcriptome diversity. Epigenetic modification is another mechanism of transcriptional control. Recent studies have uncovered crosstalk between cotranscriptional polyadenylation processes and both epigenomic and epitranscriptomic markers. Genetic analyses have demonstrated that DNA methylation, histone modifications, and epitranscriptomic modification are involved in regulating polyadenylation in plants. Here we summarize current understanding of the links between epigenetics and polyadenylation and their novel biological efficacy for plant development and environmental responses. Unresolved issues and future directions are discussed to shed light on the field.

Integrative analysis of Iso-Seq and RNA-seq reveals dynamic changes of alternative promoter, alternative splicing and alternative polyadenylation during Angiotensin II-induced senescence in rat primary aortic endothelial cells

In eukaryotes, alternative promoter (AP), alternative splicing (AS), and alternative polyadenylation (APA) are three crucial regulatory mechanisms that modulate message RNA (mRNA) diversity. Although AP, AS and APA are involved in diverse biological processess, whether they have dynamic changes in Angiotensin II (Ang II) induced senescence in rat primary aortic endothelial cells (RAECs), an important cellular model for studying cardiovascular disease, remains unclear. Here we integrated both PacBio single-molecule long-read isoform sequencing (Iso-Seq) and Illumina short-read RNA sequencing (RNA-seq) to analyze the changes of AP, AS and APA in Ang II-induced senescent RAECs. Iso-Seq generated 36,278 isoforms from 10,145 gene loci and 65.81% of these isoforms are novel, which were further cross-validated by public data obtained by other techonologies such as CAGE, PolyA-Seq and 3'READS. APA contributed most to novel isoforms, followed by AS and AP. Further investigation showed that AP, AS and APA could all contribute to the regulation of isoform, but AS has more dynamic changes compared to AP and APA upon Ang II stimulation. Genes undergoing AP, AS and APA in Ang II-treated cells are enriched in various pathways related to aging or senescence, suggesting that these molecular changes are involved in functional alterations during Ang II-induced senescence. Together, the present study largely improved the annotation of rat genome and revealed gene expression changes at isoform level, extending the understanding of the complexity of gene regulation in Ang II-treated RAECs, and also provided novel clues for discovering the regulatory mechanism undelying Ang II caused vascular senescence and diseases.

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

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

Genomic organization, intragenic tandem duplication, and expression analysis of chicken TGFBR2 gene

Transforming growth factor beta receptor Ⅱ (TGFBR2), a core member of the transforming growth factor-β (TGF-β) signaling pathway. To date, chicken TGFBR2 (cTGFBR2) genomic structure has not been fully explored. Here, the complete sequences of cTGFBR2 transcript isoforms were determined by 5′ and 3′ rapid amplification of cDNA ends (5′ & 3′ RACE) and reverse transcription polymerase chain reaction (RT-PCR); the tissue expression profiling of cTGFBR2 transcript isoforms was performed using quantitative real-time polymerase chain reaction (qRT-PCR). The results showed that cTGFBR2 gene produced 3 transcript isoforms though alternative transcription initiation, splicing, and polyadenylation, which were designated as cTGFBR2-1, cTGFBR2-2, and cTGFBR2-3, respectively. These 3 cTGFBR2 transcript isoforms encoded 3 protein isoforms: cTGFBR2-1, cTGFBR2-2, and cTGFBR2-3. Duplication analysis revealed that, unlike other animal species, cTGFBR2 gene harbored a 5.5-kb intragenic tandem duplication. Tissue expression profiling in the 4-wk-old Arbor Acres (AA) broiler chickens showed that cTGFBR2-1 was ubiquitously expressed, with high expression in abdominal fat, subcutaneous fat, lung, gizzard, and muscle; cTGFBR2-2 was highly expressed in heart, kidney, gizzard, and muscle; cTGFBR2-3 was weakly expressed in all the tested chicken tissues. Tissue expression profiling in the 7-wk-old broiler chickens of the fat and lean lines of Northeast Agricultural University broiler lines divergently selected for abdominal fat content (NEAUHLF) showed that cTGFBR2-1 was significantly differentially expressed in all the tested tissues except heart, cTGFBR2-2 was significantly differentially expressed in all the tested tissues except subcutaneous fat and liver, and cTGFBR2-3 was significantly differentially expressed in all the tested tissues between the lean and fat lines. Intriguingly, in the fat line, the 3 cTGFBR2 transcript isoforms were expressed to varying degrees in all the 3 tested fat tissues, while in the lean line, only cTGFBR2-1 was expressed in all the 3 tested fat tissues. This is the first report of intragenic tandem duplication within TGFBR2 gene. Our findings pave the way for further studies on the functions and regulation of cTGFBR2 gene.

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.

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

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

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.

Identification and analysis of HD-Zip genes involved in the leaf development of Liriodendron chinense using multidimensional analysis

Homeodomain-leucine zipper (HD-Zip) proteins are plant-specific transcription factors that play important roles in different biological processes, especially leaf development. However, no studies to date have identified the HD-Zip genes in Liriodendron chinense nor characterized their functions. We identified the HD-Zip genes in L. chinense by analysing the phylogeny, chromosome location, structure, conserved motif, cis-regulatory elements, synteny, post-transcriptional regulation and expression patterns of these genes during leaf development. A total of 36 LcHD-Zip genes were identified and divided into four subfamilies (HD-Zip I to IV). Synteny analysis revealed that segmental duplication was the main force driving the expansion of LcHD-Zip genes. These 36 LcHD-Zip genes exhibited 11 different expression patterns. Pattern 1, 2, 3, 4, 6, 7, 8 and 9 genes may play important roles in leaf development, such as leaf initiation, leaf polarity establishment, leaf shape development, phytohormone-mediated leaf growth and leaf epidermal structure formation. Four HD-Zip III genes were targeted by microRNAs (miRNAs), and the miR165/166a-HD-Zip regulatory module formed regulated leaf initiation and leaf polarity establishment. Overall, LcHD-Zip genes play key roles in leaf development of L. chinense. This work provides a foundation for the functional verification of HD-Zip genes identified in this study.

scAPAmod: Profiling Alternative Polyadenylation Modalities in Single Cells from Single-Cell RNA-Seq Data

Alternative polyadenylation (APA) is a key layer of gene expression regulation, and APA choice is finely modulated in cells. Advances in single-cell RNA-seq (scRNA-seq) have provided unprecedented opportunities to study APA in cell populations. However, existing studies that investigated APA in single cells were either confined to a few cells or focused on profiling APA dynamics between cell types or identifying APA sites. The diversity and pattern of APA usages on a genomic scale in single cells remains unappreciated. Here, we proposed an analysis framework based on a Gaussian mixture model, scAPAmod, to identify patterns of APA usage from homogeneous or heterogeneous cell populations at the single-cell level. We systematically evaluated the performance of scAPAmod using simulated data and scRNA-seq data. The results show that scAPAmod can accurately identify different patterns of APA usages at the single-cell level. We analyzed the dynamic changes in the pattern of APA usage using scAPAmod in different cell differentiation and developmental stages during mouse spermatogenesis and found that even the same gene has different patterns of APA usages in different differentiation stages. The preference of patterns of usages of APA sites in different genomic regions was also analyzed. We found that patterns of APA usages of the same gene in 3′ UTRs (3′ untranslated region) and non-3′ UTRs are different. Moreover, we analyzed cell-type-specific APA usage patterns and changes in patterns of APA usages across cell types. Different from the conventional analysis of single-cell heterogeneity based on gene expression profiling, this study profiled the heterogeneous pattern of APA isoforms, which contributes to revealing the heterogeneity of single-cell gene expression with higher resolution.

3’UTR-Seq analysis of chicken abdominal adipose tissue reveals widespread intron retention in 3’UTR and provides insight into molecular basis of feed efficiency

Feed efficiency (FE) is an important trait in the broiler industry due to its direct correlation to efficient muscle growth instead of fat deposition. The present study characterized and compared gene expression profiles in abdominal fat from broiler chickens of different FE levels to enhance the understanding of FE biology. Specifically, traditional whole-transcript RNA-sequencing (RNA-seq) and 3’ UTR-sequencing (3’ UTR-seq) were applied to 22 and 61 samples, respectively. Overall, these two sequencing techniques shared a high correlation (0.76) between normalized counts, although 3’ UTR-seq showed a higher variance in sequencing and mapping performance statistics across samples and a lower rate of uniquely mapped reads. A higher percentage of 3’ UTR-seq reads mapped to introns suggested the frequent presence of cleavage sites in introns, thus warranting future research to study its regulatory function. Differential expression analysis identified 1198 differentially expressed genes (DEGs) between high FE (HFE) and intermediate FE (IFE) chickens with False Discovery Rate < 0.05 and fold change > 1.2. The processes that were significantly enriched by the DEGs included extracellular matrix remodeling and mechanisms impacting gene expression at the transcriptional and translational levels. Gene ontology enrichment analysis suggested that the divergence in fat deposition and FE in broiler chickens could be associated with peroxisome and lipid metabolism possibly regulated by G0/G1 switch gene 2 (G0S2).

RNA polyadenylation patterns in the human transcriptome

The eukaryotic transcriptome undergoes various post-transcriptional modifications which assists gene expression. Polyadenylation is a molecular process occurring at the 3'-end of the RNA molecule which involves the poly(A) polymerase attaching adenine monophosphate molecules in a chain-like fashion to assemble a poly(A) tail. Multiple RNA isoforms are produced with differing 3'-UTR and exonic compositions through alternative polyadenylation (APA) which enhances the diversification of alternatively spliced mRNA transcripts. To study polyadenylation patterns, novel methods have been developed using short-read and long-read sequencing technologies to analyse the 3'-ends of the transcript. Recent studies have identified unique polyadenylation patterns in different cellular functions, including oncogenic activity, which could prove valuable in the understanding of medical genetics, particularly in the discovery of biomarkers in diseased states. We present a review of current literature reporting on polyadenylation and the biological relevance in the mammalian transcriptome, with a focus on the human transcriptome. Additionally, we have explored the various methods available to detect polyadenylation patterns using second and third generation sequencing technologies.

MDM4 alternative splicing and implication in MDM4 targeted cancer therapies

The oncogenic MDM4, initially named MDMX, has been identified as a p53-interacting protein and a key upstream negative regulator of the tumor suppressor p53. Accumulating evidence indicates that MDM4 plays critical roles in the initiation and progression of multiple human cancers. MDM4 is frequently amplified and upregulated in human cancers, contributing to overgrowth and apoptosis inhibition by blocking the expression of downstream target genes of p53 pathway. Disruptors for MDM4-p53 interaction have been shown to restore the anti-tumor activity of p53 in cancer cells. MDM4 possesses multiple splicing isoforms whose expressions are driven by the presence of oncogenes in cancer cells. Some of the MDM4 splicing isoforms lack p53 binding domain and may exhibit p53-independent oncogenic functions. These features render MDM4 to be an attractive therapeutic target for cancer therapy. In the present review, we primarily focus on the detailed molecular structure of MDM4 splicing isoforms, candidate regulators for initiating MDM4 splicing, deregulation of MDM4 isoforms in cancer and potential therapy strategies by targeting splicing isoforms of MDM4.

Genome-Wide Analysis of Alternative Splicing during Host-Virus Interactions in Chicken

The chicken is a model animal for the study of evolution, immunity and development. In addition to their use as a model organism, chickens also represent an important agricultural product. Pathogen invasion has already been shown to modulate the expression of hundreds of genes, but the role of alternative splicing in avian virus infection remains unclear. We used RNA-seq data to analyze virus-induced changes in the alternative splicing of Gallus gallus, and found that a large number of alternative splicing events were induced by virus infection both in vivo and in vitro. Virus-responsive alternative splicing events preferentially occurred in genes involved in metabolism and transport. Many of the alternatively spliced transcripts were also expressed from genes with a function relating to splicing or immune response, suggesting a potential impact of virus infection on pre-mRNA splicing and immune gene regulation. Moreover, exon skipping was the most frequent AS event in chickens during virus infection. This is the first report describing a genome-wide analysis of alternative splicing in chicken and contributes to the genomic resources available for studying host-virus interaction in this species. Our analysis fills an important knowledge gap in understanding the extent of genome-wide alternative splicing dynamics occurring during avian virus infection and provides the impetus for the further exploration of AS in chicken defense signaling and homeostasis.

Alternative polyadenylated mRNAs behave as asynchronous rhythmic transcription in Arabidopsis

Alternative polyadenylation (APA) is a widespread post-transcriptional modification method that changes the 3' ends of transcripts by altering poly(A) site usage. However, the longitudinal transcriptomic 3' end profile and its mechanism of action are poorly understood. We applied diurnal time-course poly(A) tag sequencing (PAT-seq) for Arabidopsis and identified 3284 genes that generated both rhythmic and arrhythmic transcripts. These two classes of transcripts appear to exhibit dramatic differences in expression and translation activisty. The asynchronized transcripts derived by APA are embedded with different poly(A) signals, especially for rhythmic transcripts, which contain higher AAUAAA and UGUA signal proportions. The Pol II occupancy maximum is reached upstream of rhythmic poly(A) sites, while it is present directly at arrhythmic poly(A) sites. Integrating H3K9ac and H3K4me3 time-course data analyses revealed that transcriptional activation of histone markers may be involved in the differentiation of rhythmic and arrhythmic APA transcripts. These results implicate an interplay between histone modification and RNA 3'-end processing, shedding light on the mechanism of transcription rhythm and alternative polyadenylation.

Alternative splicing landscape of small brown planthopper and different response of JNK2 isoforms to rice stripe virus infection

Alternative splicing (AS) is a frequent posttranscriptional regulatory event occurring in response to various endogenous and exogenous stimuli in most eukaryotic organisms. However, little is known about the effects of insect-transmitted viruses on AS events in insect vectors. The present study used third-generation sequencing technology and RNA sequencing (RNA-Seq) to evaluate the AS response in the small brown planthopper Laodelphax striatellus to rice stripe virus (RSV). The full-length transcriptome of L. striatellus was obtained using single-molecule real-time sequencing technology (SMRT). Posttranscriptional regulatory events, including AS, alternative polyadenylation, and fusion transcripts, were analyzed. A total of 28,175 nonredundant transcript isoforms included 24,950 transcripts assigned to 8,500 annotated genes of L. striatellus, and 5,000 of these genes (58.8%) had AS events. RNA-Seq of the gut samples of insects infected by RSV for 8 d identified 3,458 differentially expressed transcripts (DETs); 2,185 of these DETs were transcribed from 1,568 genes that had AS events, indicating that 31.4% of alternatively spliced genes responded to RSV infection of the gut. One of the c-Jun N-terminal kinase (JNK) genes, JNK2, experienced exon skipping, resulting in three transcript isoforms. These three isoforms differentially responded to RSV infection during development and in various organs. Injection of double-stranded RNAs targeting all or two isoforms indicated that three or at least two JNK2 isoforms facilitated RSV accumulation in planthoppers. These results implied that AS events could participate in the regulation of complex relationships between viruses and insect vectors. Importance Alternative splicing (AS) is a regulatory mechanism that occurs after gene transcription. AS events can enrich protein diversity to promote the reactions of the organisms to various endogenous and exogenous stimulations. It is not known how insect vectors exploit AS events to cope with transmitted viruses. The present study used third-generation sequencing technology to obtain the profile of AS events in the small brown planthopper Laodelphax striatellus, which is an efficient vector for rice stripe virus (RSV). The results indicated that 31.4% of alternatively spliced genes responded to RSV infection in the gut of planthoppers. One of the c-Jun N-terminal kinase (JNK) genes, JNK2, produced three transcript isoforms by AS. These three isoforms showed different responses to RSV infection, and at least two isoforms facilitated viral accumulation in planthoppers. These results implied that AS events could participate in the regulation of complex relationships between viruses and insect vectors.

Lung cancer cells expressing a shortened CDK16 3'UTR escape senescence through impaired miR-485-5p targeting

Inducing senescence in cancer cells is an emerging strategy for cancer therapy. The dysregulation and mutation of genes encoding cyclin‐dependent kinases (CDKs) have been implicated in various human cancers. However, whether CDK can induce cancer cell senescence remains poorly understood. We observed that CDK16 expression was high in multiple cancer types, including lung cancer, whereas various replicative senescence models displayed low CDK16 expression. CDK16 knockdown caused senescence‐associated phenotypes in lung cancer cell lines. Interestingly, the CDK16 3′ UTR was shortened in cancer and lengthened in senescence models, which was regulated by alternative polyadenylation (APA). The longer 3′UTR [using the distal polyA (pA) site] generated less protein than the shorter one (using the proximal pA site). Since microRNAs (miRNAs) usually bind to the 3′UTR of target genes to suppress their expression, we investigated whether miRNAs targeting the region between the shortened and longer 3′UTR are responsible for the reduced expression. We found that miR‐485‐5p targeted the 3′UTR between the distal and proximal pA site and caused senescence‐associated phenotypes by reducing protein production from the longer CDK16 transcript. Of note, CDK16 knockdown led to a reduced expression of MYC proto‐oncogene, bHLH transcription factor (MYC) and CD274 molecule (PD‐L1), which in turn enhanced the tumor‐suppressive effects of senescent cancer cells. The present study discovered that CDK16, whose expression is under the regulation of APA and miR‐485‐5p, is a potential target for prosenescence therapy for lung cancer.

Combining protein and RNA quantification to evaluate promoter activity by using dual-color fluorescent reporting systems

Herein, Broccoli/mCherry and an EGFP/mCherry dual-color fluorescent reporting systems have been established to quantify the promoter activity at transcription and translation levels in eukaryotic cells. Based on those systems, four commonly used promoters (CMV and SV40 of Pol II and U6, H1 of Pol III) were accurately evaluated at both the transcriptional and translational levels by combining accurate protein and RNA quantification. Furthermore, we verified that Pol III promoters can induce proteins expression, and Pol II promoter can be applied to express RNA molecules with defined length by combining a self-cleaving ribozyme and an artificial poly(A) tail. The dual-color fluorescence reporting systems described here could play a significant role in evaluating other gene expression regulators for gene therapy.

Transcriptomic analysis to elucidate the effects of high stocking density on grass carp (Ctenopharyngodon idella)

Background

Grass carp (Ctenopharyngodon idella) is one of the most widely cultivated fishes in China. High stocking density can reportedly affect fish growth and immunity. Herein we performed PacBio long-read single-molecule real-time (SMRT) sequencing and Illumina RNA sequencing to evaluate the effects of high stocking density on grass carp transcriptome.

Results

SMRT sequencing led to the identification of 33,773 genes (14,946 known and 18,827 new genes). From the structure analysis, 8,009 genes were detected with alternative splicing events, 10,219 genes showed alternative polyadenylation sites and 15,521 long noncoding RNAs. Further, 1,235, 962, and 213 differentially expressed genes (DEGs) were identified in the intestine, muscle, and brain tissues, respectively. We performed functional enrichment analyses of DEGs, and they were identified to be significantly enriched in nutrient metabolism and immune function. The expression levels of several genes encoding apolipoproteins and activities of enzymes involved in carbohydrate enzymolysis were found to be upregulated in the high stocking density group, indicating that lipid metabolism and carbohydrate decomposition were accelerated. Besides, four isoforms of grass carp major histocompatibility complex class II antigen alpha and beta chains in the aforementioned three tissue was showed at least a 4-fold decrease.

Conclusions

The results suggesting that fish farmed at high stocking densities face issues associated with the metabolism and immune system. To conclude, our results emphasize the importance of maintaining reasonable density in grass carp aquaculture.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07924-4.

QuantifyPoly(A): reshaping alternative polyadenylation landscapes of eukaryotes with weighted density peak clustering

The dynamic choice of different polyadenylation sites in a gene is referred to as alternative polyadenylation, which functions in many important biological processes. Large-scale messenger RNA 3' end sequencing has revealed that cleavage sites for polyadenylation are presented with microheterogeneity. To date, the conventional determination of polyadenylation site clusters is subjective and arbitrary, leading to inaccurate annotations. Here, we present a weighted density peak clustering method, QuantifyPoly(A), to accurately quantify genome-wide polyadenylation choices. Applying QuantifyPoly(A) on published 3' end sequencing datasets from both animals and plants, their polyadenylation profiles are reshaped into myriads of novel polyadenylation site clusters. Most of these novel polyadenylation site clusters show significantly dynamic usage across different biological samples or associate with binding sites of trans-acting factors. Upstream sequences of these clusters are enriched with polyadenylation signals UGUA, UAAA and/or AAUAAA in a species-dependent manner. Polyadenylation site clusters also exhibit species specificity, while plants ones generally show higher microheterogeneity than that of animals. QuantifyPoly(A) is broadly applicable to any types of 3' end sequencing data and species for accurate quantification and construction of the complex and dynamic polyadenylation landscape and enables us to decode alternative polyadenylation events invisible to conventional methods at a much higher resolution.

Construction of Two Alternative Polyadenylation Signatures to Predict the Prognosis of Sarcoma Patients

Background

Increasing evidence indicates that alternative polyadenylation (APA) is associated with the prognosis of cancers.

Methods

We obtained gene expression and APA profiles of 259 sarcoma patients from the TCGA dataportal and TC3A database, respectively. The prognostic signatures, clinical nomograms, and regulatory networks were studied by integrated bioinformatics analyses. Then, the immune cell infiltration profile was obtained from the ImmuCellAI. The association between APA-based signature and immune cells was studied.

Results

A total of 61 and 38 APA events were identified as overall survival (OS)- and progress free-survival (PFS)-related biomarkers, respectively. Two signatures were generated. The area under the curves (AUC) values of OS signature were 0.900, 0.928, and 0.963 over 2-, 4-, and 6-years, respectively. And the AUC values of PFS signature at 2-, 4-, and 6-years were 0.826, 0.840, and 0.847, respectively. Overall and subgroup analyses indicated that high-risk patients had a worse prognosis than low-risk patients (all p-values < 0.05). In addition, immunomics analyses indicated that there are different patterns of immune cell infiltration between low- and high-risk patients. Furthermore, two clinical-APA nomograms were established and the C-indexes were 0.813 and 0.809 for OS nomogram and PFS nomogram, respectively. Finally, two APA regulatory networks were constructed. FIP1L1-VPS26B was identified as a key regulating relationship and validated in the pan-cancer analyses.

Conclusion

In this study, we identified prognostic predictors based on APA events with high accuracy for risk stratification in sarcoma patients and uncovered interesting regulatory networks in sarcoma that could be underlying mechanisms. This study not only provides novel potential prognostic biomarkers but promote precision medicine and provide potential novel research interests for immunotherapy.

QPAT-seq, a rapid and deduplicatable method for quantification of poly(A) site usages

Alternative polyadenylation (APA) is an essential regulatory mechanism for gene expression. The next generation sequencing provides ample opportunity to precisely delineate APA sites genome-wide. Various methods for profiling transcriptome-wide poly(A) sites were developed. By comparing available methods, the ways for adding sequencing adaptors to fit with the Illumina sequencing platform are different. These methods have identified more than 50% genes that undergo APA in eukaryotes. However, due to the unbalanced PCR during library preparation, accurate quantification of poly(A) sites is still a challenge. Here, we describe an updated poly(A) tag sequencing method that incorporates unique molecular identifier (UMI) into the adaptor for removing quantification bias induced by PCR duplicates. Hence, quantification of poly(A) site usages can be achieved by counting UMIs. This protocol, quantifying poly(A) tag sequencing (QPAT-seq), can be finished in 1 day with reduced cost, and is particularly useful for application with a large number of samples.

A computational pipeline to infer alternative poly-adenylation from 3' sequencing data

An increasing number of investigations have established alternative polyadenylation (APA) as a key mechanism of gene regulation through altering the length of 3' untranslated region (UTR) and generating distinct mRNA termini. Further, appreciation for the significance of APA in disease contexts propelled the development of several 3' sequencing techniques. While these RNA sequencing technologies have advanced APA analysis, the intrinsic limitation of 3' read coverage and lack of appropriate computational tools constrain precise mapping and quantification of polyadenylation sites. Notably, Poly(A)-ClickSeq (PAC-seq) overcomes limiting factors such as poly(A) enrichment and 3' linker ligation steps using click-chemistry. Here we provide an updated PolyA-miner protocol, a computational approach to analyze PAC-seq or other 3'-Seq datasets. As a key practical constraint, we also provide a detailed account on the impact of sequencing depth on the number of detected polyadenylation sites and APA changes. This protocol is also updated to handle unique molecular identifiers used to address PCR duplication potentially observed in PAC-seq.

Analysis of alternative polyadenylation from single-cell RNA-seq using scDaPars reveals cell subpopulations invisible to gene expression

Alternative polyadenylation (APA) is a major mechanism of post-transcriptional regulation in various cellular processes including cell proliferation and differentiation, but the APA heterogeneity among single cells remains largely unknown. Single-cell RNA sequencing (scRNA-seq) has been extensively used to define cell subpopulations at the transcription level. Yet, most scRNA-seq data have not been analyzed in an "APA-aware" manner. Here, we introduce scDaPars (Dynamic Analysis of Alternative PolyAdenylation from Single-cell RNA-seq), a bioinformatics algorithm to accurately quantify APA events at both single-cell and single-gene resolution using either 3' end (10x Chromium) or full-length (Smart-seq2) scRNA-seq data. Validations in both real and simulated data indicate that scDaPars can robustly recover missing APA events caused by the low amounts of mRNA sequenced in single cells. When applied to cancer and human endoderm differentiation data, scDaPars not only revealed cell type-specific APA regulation but also identified cell subpopulations that are otherwise invisible to conventional gene expression analysis. Thus, scDaPars will enable us to understand cellular heterogeneity at the post-transcriptional APA level.

The Detection and Bioinformatic Analysis of Alternative 3' UTR Isoforms as Potential Cancer Biomarkers

Alternative transcript cleavage and polyadenylation is linked to cancer cell transformation, proliferation and outcome. This has led researchers to develop methods to detect and bioinformatically analyse alternative polyadenylation as potential cancer biomarkers. If incorporated into standard prognostic measures such as gene expression and clinical parameters, these could advance cancer prognostic testing and possibly guide therapy. In this review, we focus on the existing methodologies, both experimental and computational, that have been applied to support the use of alternative polyadenylation as cancer biomarkers.

Alternative Polyadenylation in Stem Cell Self-Renewal and Differentiation

Cellular function is shaped by transcriptional and post-transcriptional mechanisms, including alternative polyadenylation (APA). By directly controlling 3'- untranslated region (UTR) length and the selection of the last exon, APA regulates up to 70% of all cellular transcripts influencing RNA stability, output, and protein isoform expression. Cell-state-dependent 3'-UTR shortening has been identified as a hallmark of cellular proliferation. Hence, quiescent/dormant stem cells are characterized by long 3'-UTRs, whereas proliferative stem/progenitor cells exhibit 3'-UTR shortening. Here, the latest studies analyzing the role of APA in regulating stem cell state, self-renewal, differentiation, and metabolism are reviewed. The new role of APA in controlling stem cell fate opens novel potential therapeutic avenues in the field of regenerative medicine.

Systematic refinement of gene annotations by parsing mRNA 3' end sequencing datasets

Alternative cleavage and polyadenylation generates mRNA 3' isoforms in a cell type-specific manner. Due to finite available RNA sequencing data of organisms with vast cell type complexity, currently available gene annotation resources are incomplete, which poses significant challenges to the comprehensive interpretation and quantification of transcriptomes. In this chapter, we introduce 3'GAmES, a stand-alone computational pipeline for the identification and quantification of novel mRNA 3'end isoforms from 3'mRNA sequencing data. 3'GAmES expands available repositories and improves comprehensive gene-tag counting by cost-effective 3' mRNA sequencing, faithfully mirroring whole-transcriptome RNAseq measurements. By employing R and bash shell scripts (assembled in a Singularity container) 3'GAmES systematically augments cell type-specific 3' ends of RNA polymerase II transcripts and increases the sensitivity of quantitative gene expression profiling by 3' mRNA sequencing. Public access: https://github.com/AmeresLab/3-GAmES.git.

A Tissue-Specific Landscape of Alternative Polyadenylation, lncRNAs, TFs, and Gene Co-expression Networks in Liriodendron chinense

Liriodendron chinense is an economically and ecologically important deciduous tree species. Although the reference genome has been revealed, alternative polyadenylation (APA), transcription factors (TFs), long non-coding RNAs (lncRNAs), and co-expression networks of tissue-specific genes remain incompletely annotated. In this study, we used the bracts, petals, sepals, stamens, pistils, leaves, and shoot apex of L. chinense as materials for hybrid sequencing. On the one hand, we improved the annotation of the genome. We detected 13,139 novel genes, 7,527 lncRNAs, 1,791 TFs, and 6,721 genes with APA sites. On the other hand, we found that tissue-specific genes play a significant role in maintaining tissue characteristics. In total, 2,040 tissue-specific genes were identified, among which 9.2% of tissue-specific genes were affected by APA, and 1,809 tissue-specific genes were represented in seven specific co-expression modules. We also found that bract-specific hub genes were associated plant defense, leaf-specific hub genes were involved in energy metabolism. Moreover, we also found that a stamen-specific hub TF Lchi25777 may be involved in the determination of stamen identity, and a shoot-apex-specific hub TF Lchi05072 may participate in maintaining meristem characteristic. Our study provides a landscape of APA, lncRNAs, TFs, and tissue-specific gene co-expression networks in L. chinense that will improve genome annotation, strengthen our understanding of transcriptome complexity, and drive further research into the regulatory mechanisms of tissue-specific genes.

Aptamer-Mediated Control of Polyadenylation for Gene Expression Regulation in Mammalian Cells

Small aptamer-based regulatory devices can be designed to control a range of RNA-dependent cellular processes and emerged as promising tools for fine-tuning gene expression in synthetic biology. Here, we design a conceptually new riboswitch device that allows for the conditional regulation of polyadenylation. By making use of ligand-induced sequence occlusion, the system efficiently controls the accessibility of the eukaryotic polyadenylation signal. Undesirable 3'-extended read-through products are counteracted by the downstream insertion of a microRNA target site. We demonstrate the modularity of the system with regard to sensor aptamers and polyadenylation signals used and combine the newly designed riboswitch with well-known aptazymes to yield superior composite systems. In addition, we show that the switches can be used to control alternative polyadenylation. The presented genetic switches require very little coding space and can be easily optimized by rational adjustments of the thermodynamic stability. The polyadenylation riboswitch extends the repertoire of RNA-based regulators and opens new possibilities for the generation of complex synthetic circuits.

scAPAtrap: identification and quantification of alternative polyadenylation sites from single-cell RNA-seq data

Alternative polyadenylation (APA) generates diverse mRNA isoforms, which contributes to transcriptome diversity and gene expression regulation by affecting mRNA stability, translation and localization in cells. The rapid development of 3' tag-based single-cell RNA-sequencing (scRNA-seq) technologies, such as CEL-seq and 10x Genomics, has led to the emergence of computational methods for identifying APA sites and profiling APA dynamics at single-cell resolution. However, existing methods fail to detect the precise location of poly(A) sites or sites with low read coverage. Moreover, they rely on priori genome annotation and can only detect poly(A) sites located within or near annotated genes. Here we proposed a tool called scAPAtrap for detecting poly(A) sites at the whole genome level in individual cells from 3' tag-based scRNA-seq data. scAPAtrap incorporates peak identification and poly(A) read anchoring, enabling the identification of the precise location of poly(A) sites, even for sites with low read coverage. Moreover, scAPAtrap can identify poly(A) sites without using priori genome annotation, which helps locate novel poly(A) sites in previously overlooked regions and improve genome annotation. We compared scAPAtrap with two latest methods, scAPA and Sierra, using scRNA-seq data from different experimental technologies and species. Results show that scAPAtrap identified poly(A) sites with higher accuracy and sensitivity than competing methods and could be used to explore APA dynamics among cell types or the heterogeneous APA isoform expression in individual cells. scAPAtrap is available at https://github.com/BMILAB/scAPAtrap.

Single-cell RNA cap and tail sequencing (scRCAT-seq) reveals subtype-specific isoforms differing in transcript demarcation

The differences in transcription start sites (TSS) and transcription end sites (TES) among gene isoforms can affect the stability, localization, and translation efficiency of mRNA. Gene isoforms allow a single gene diverse functions across different cell types, and isoform dynamics allow different functions over time. However, methods to efficiently identify and quantify RNA isoforms genome-wide in single cells are still lacking. Here, we introduce single cell RNA Cap And Tail sequencing (scRCAT-seq), a method to demarcate the boundaries of isoforms based on short-read sequencing, with higher efficiency and lower cost than existing long-read sequencing methods. In conjunction with machine learning algorithms, scRCAT-seq demarcates RNA transcripts with unprecedented accuracy. We identified hundreds of previously uncharacterized transcripts and thousands of alternative transcripts for known genes, revealed cell-type specific isoforms for various cell types across different species, and generated a cell atlas of isoform dynamics during the development of retinal cones.

scDAPA: detection and visualization of dynamic alternative polyadenylation from single cell RNA-seq data

MOTIVATION

Alternative polyadenylation (APA) plays a key post-transcriptional regulatory role in mRNA stability and functions in eukaryotes. Single cell RNA-seq (scRNA-seq) is a powerful tool to discover cellular heterogeneity at gene expression level. Given 3' enriched strategy in library construction, the most commonly used scRNA-seq protocol-10× Genomics enables us to improve the study resolution of APA to the single cell level. However, currently there is no computational tool available for investigating APA profiles from scRNA-seq data.

RESULTS

Here, we present a package scDAPA for detecting and visualizing dynamic APA from scRNA-seq data. Taking bam/sam files and cell cluster labels as inputs, scDAPA detects APA dynamics using a histogram-based method and the Wilcoxon rank-sum test, and visualizes candidate genes with dynamic APA. Benchmarking results demonstrated that scDAPA can effectively identify genes with dynamic APA among different cell groups from scRNA-seq data.

AVAILABILITY AND IMPLEMENTATION

The scDAPA package is implemented in Shell and R, and is freely available at https://scdapa.sourceforge.io.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

PolyA-miner: accurate assessment of differential alternative poly-adenylation from 3'Seq data using vector projections and non-negative matrix factorization

Almost 70% of human genes undergo alternative polyadenylation (APA) and generate mRNA transcripts with varying lengths, typically of the 3′ untranslated regions (UTR). APA plays an important role in development and cellular differentiation, and its dysregulation can cause neuropsychiatric diseases and increase cancer severity. Increasing awareness of APA’s role in human health and disease has propelled the development of several 3′ sequencing (3′Seq) techniques that allow for precise identification of APA sites. However, despite the recent data explosion, there are no robust computational tools that are precisely designed to analyze 3′Seq data. Analytical approaches that have been used to analyze these data predominantly use proximal to distal usage. With about 50% of human genes having more than two APA isoforms, current methods fail to capture the entirety of APA changes and do not account for non-proximal to non-distal changes. Addressing these key challenges, this study demonstrates PolyA-miner, an algorithm to accurately detect and assess differential alternative polyadenylation specifically from 3′Seq data. Genes are abstracted as APA matrices, and differential APA usage is inferred using iterative consensus non-negative matrix factorization (NMF) based clustering. PolyA-miner accounts for all non-proximal to non-distal APA switches using vector projections and reflects precise gene-level 3′UTR changes. It can also effectively identify novel APA sites that are otherwise undetected when using reference-based approaches. Evaluation on multiple datasets—first-generation MicroArray Quality Control (MAQC) brain and Universal Human Reference (UHR) PolyA-seq data, recent glioblastoma cell line NUDT21 knockdown Poly(A)-ClickSeq (PAC-seq) data, and our own mouse hippocampal and human stem cell-derived neuron PAC-seq data—strongly supports the value and protocol-independent applicability of PolyA-miner. Strikingly, in the glioblastoma cell line data, PolyA-miner identified more than twice the number of genes with APA changes than initially reported. With the emerging importance of APA in human development and disease, PolyA-miner can significantly improve data analysis and help decode the underlying APA dynamics.

Gene expression dynamics are a proxy for selective pressures on alternatively polyadenylated isoforms

Alternative polyadenylation (APA) produces isoforms with distinct 3′-ends, yet their functional differences remain largely unknown. Here, we introduce the APA-seq method to detect the expression levels of APA isoforms from 3′-end RNA-Seq data by exploiting both paired-end reads for gene isoform identification and quantification. We detected the expression levels of APA isoforms in individual Caenorhabditis elegans embryos at different stages throughout embryogenesis. Examining the correlation between the temporal profiles of isoforms led us to distinguish two classes of genes: those with highly correlated isoforms (HCI) and those with lowly correlated isoforms (LCI) across time. We hypothesized that variants with similar expression profiles may be the product of biological noise, while the LCI variants may be under tighter selection and consequently their distinct 3′ UTR isoforms are more likely to have functional consequences. Supporting this notion, we found that LCI genes have significantly more miRNA binding sites, more correlated expression profiles with those of their targeting miRNAs and a relative lack of correspondence between their transcription and protein abundances. Collectively, our results suggest that a lack of coherence among the regulation of 3′ UTR isoforms is a proxy for selective pressures acting upon APA usage and consequently for their functional relevance.

Self-Reporting Transposons Enable Simultaneous Readout of Gene Expression and Transcription Factor Binding in Single Cells

Cellular heterogeneity confounds in situ assays of transcription factor (TF) binding. Single-cell RNA sequencing (scRNA-seq) deconvolves cell types from gene expression, but no technology links cell identity to TF binding sites (TFBS) in those cell types. We present self-reporting transposons (SRTs) and use them in single-cell calling cards (scCC), a novel assay for simultaneously measuring gene expression and mapping TFBS in single cells. The genomic locations of SRTs are recovered from mRNA, and SRTs deposited by exogenous, TF-transposase fusions can be used to map TFBS. We then present scCC, which map SRTs from scRNA-seq libraries, simultaneously identifying cell types and TFBS in those same cells. We benchmark multiple TFs with this technique. Next, we use scCC to discover BRD4-mediated cell-state transitions in K562 cells. Finally, we map BRD4 binding sites in the mouse cortex at single-cell resolution, establishing a new method for studying TF biology in situ.

HNRNPA1-mediated 3' UTR length changes of HN1 contributes to cancer- and senescence-associated phenotypes

Cellular senescence has been regarded as a mechanism of tumor suppression. Studying the regulation of gene expression at various levels in cell senescence will shed light on cancer therapy. Alternative polyadenylation (APA) regulates gene expression by altering 3′ untranslated regions (3′ UTR) and plays important roles in diverse biological processes. However, whether APA of a specific gene functions in both cancer and senescence remains unclear. Here, we discovered that 3′ UTR of HN1 (or JPT1) showed shortening in cancers and lengthening in senescence, correlated well with its high expression in cancer cells and low expression in senescent cells, respectively. HN1 transcripts with longer 3′ UTR were less stable and produced less protein. Down-regulation of HN1 induced senescence-associated phenotypes in both normal and cancer cells. Patients with higher HN1 expression had lower survival rates in various carcinomas. Interestingly, down-regulating the splicing factor HNRNPA1 induced 3′ UTR lengthening of HN1 and senescence-associated phenotypes, which could be partially reversed by overexpressing HN1. Together, we revealed for the first time that HNRNPA1-mediated APA of HN1 contributed to cancer- and senescence-related phenotypes. Given senescence is a cancer prevention mechanism, our discovery indicates the HNRNPA1-HN1 axis as a potential target for cancer treatment.

Involved microRNAs in alternative polyadenylation intervene in breast cancer via regulation of cleavage factor "CFIm25"

Cleavage factor “CFIm25”, as a key repressor at proximal poly (A) site, negatively correlates to cell proliferation and tumorigenicity in various cancers. Hence, understanding CFIm25 mechanism of action in breast cancer would be a great benefit. To this aim four steps were designed. First, potential miRNAs that target 3′-UTR of CFIm25 mRNA, retrieved from Targetscan web server. Second, screened miRNAs were profiled in 100 breast cancer and 100 normal adjacent samples. Third, miRNAs that their expression was inversely correlated to the CFIm25, overexpressed in MDA-MB-231 cell line, and their effect on proliferation and migration monitored via MTT and wound healing assays, respectively. Fourth, interaction of miRNAs of interest with 3′-UTR of CFIm25 confirmed via luciferase assay and western blot. Our results indicate that CFIm25 considerably down-regulates in human breast cancer tissue. qRT-PCR assay, luciferase test, and western blotting confirm that CFIm25 itself could be directly regulated by oncomiRs such as miR-23, -24, -27, -135, -182 and -374. Besides, according to MTT and wound healing assays of cell lines, CFIm25 knockdown intensifies cell growth, proliferation and migration. Our results also confirm indirect impact of CFIm25 on regulation of mRNA’s 3′–UTR length, which then control corresponding miRNAs’ action. miRNAs directly control CFIm25 expression level, which then tunes expression of the oncogenes and tumor proliferation. Therefore, regulation of CFIm25 expression level via miRNAs is expected to improve treatment responses in breast cancer.

Nuclear eIF4E Stimulates 3'-End Cleavage of Target RNAs

The eukaryotic translation initiation factor eIF4E is nuclear and cytoplasmic where it plays roles in export and translation of specific transcripts, respectively. When we were studying its mRNA export activity, we unexpectedly discovered that eIF4E drives the protein expression of elements of the 3′-end core cleavage complex involved in cleavage and polyadenylation (CPA), including CPSF3, the enzyme responsible for cleavage, as well as its co-factors CPSF1, CPSF2, CPSF4, Symplekin, WDR33, and FIP1L1. Using multiple strategies, we demonstrate that eIF4E stimulates 3′-end cleavage of selected RNAs. eIF4E physically interacts with CPSF3, CPSF1, and uncleaved target RNA, suggesting it acts directly and indirectly on the pathway. Through these effects, eIF4E can generate better substrates for its mRNA export and translation activities. Thus, we identified an unanticipated function for eIF4E in 3′-end processing of specific target RNAs, and this function could potentially affect the expression of a broad range of oncoproteins.

Davis et al. demonstrate that the eukaryotic translation initiation factor eIF4E, which is usually associated with nuclear export and translation of specific transcripts, also acts in 3′-end processing of selected RNAs. Through these effects, eIF4E can generate better substrates for its export and translation activities and, thus, modulate the proteome.

Alternative polyadenylation dependent function of splicing factor SRSF3 contributes to cellular senescence

Down-regulated splicing factor SRSF3 is known to promote cellular senescence, an important biological process in preventing cancer and contributing to individual aging, via its alternative splicing dependent function in human cells. Here we discovered alternative polyadenylation (APA) dependent function of SRSF3 as a novel mechanism explaining SRSF3 downregulation induced cellular senescence. Knockdown of SRSF3 resulted in preference usage of proximal poly(A) sites and thus global shortening of 3′ untranslated regions (3′ UTRs) of mRNAs. SRSF3-depletion also induced senescence-related phenotypes in both human and mouse cells. These 3′ UTR shortened genes were enriched in senescence-associated pathways. Shortened 3′ UTRs tended to produce more proteins than the longer ones. Simulating the effects of 3′ UTR shortening by overexpression of three candidate genes (PTEN, PIAS1 and DNMT3A) all led to senescence-associated phenotypes. Mechanistically, SRSF3 has higher binding density near proximal poly(A) site than distal one in 3′ UTR shortened genes. Further, upregulation of PTEN by either ectopic overexpression or SRSF3-knockdown induction both led to reduced phosphorylation of AKT and ultimately senescence-associated phenotypes. We revealed for the first time that reduced SRSF3 expression could promote cellular senescence through its APA-dependent function, largely extending our mechanistic understanding in splicing factor regulated cellular senescence.

Aire-dependent genes undergo Clp1-mediated 3'UTR shortening associated with higher transcript stability in the thymus

The ability of the immune system to avoid autoimmune disease relies on tolerization of thymocytes to self-antigens whose expression and presentation by thymic medullary epithelial cells (mTECs) is controlled predominantly by Aire at the transcriptional level and possibly regulated at other unrecognized levels. Aire-sensitive gene expression is influenced by several molecular factors, some of which belong to the 3'end processing complex, suggesting they might impact transcript stability and levels through an effect on 3'UTR shortening. We discovered that Aire-sensitive genes display a pronounced preference for short-3'UTR transcript isoforms in mTECs, a feature preceding Aire's expression and correlated with the preferential selection of proximal polyA sites by the 3'end processing complex. Through an RNAi screen and generation of a lentigenic mouse, we found that one factor, Clp1, promotes 3'UTR shortening associated with higher transcript stability and expression of Aire-sensitive genes, revealing a post-transcriptional level of control of Aire-activated expression in mTECs.

Discovery of alternative polyadenylation dynamics from single cell types

Alternative polyadenylation (APA) occurs in the process of mRNA maturation by adding a poly(A) tail at different locations, resulting increased diversity of mRNA isoforms and contributing to the complexity of gene regulatory network. Benefit from the development of high-throughput sequencing technologies, we could now delineate APA profiles of transcriptomes at an unprecedented pace. Especially the single cell RNA sequencing (scRNA-seq) technologies provide us opportunities to interrogate biological details of diverse and rare cell types. Despite increasing evidence showing that APA is involved in the cell type-specific regulation and function, efficient and specific laboratory methods for capturing poly(A) sites at single cell resolution are underdeveloped to date. In this review, we summarize existing experimental and computational methods for the identification of APA dynamics from diverse single cell types. A future perspective is also provided.

The transcription and export complex THO/TREX contributes to transcription termination in plants

Transcription termination has important regulatory functions, impacting mRNA stability, localization and translation potential. Failure to appropriately terminate transcription can also lead to read-through transcription and the synthesis of antisense RNAs which can have profound impact on gene expression. The Transcription-Export (THO/TREX) protein complex plays an important role in coupling transcription with splicing and export of mRNA. However, little is known about the role of the THO/TREX complex in the control of transcription termination. In this work, we show that two proteins of the THO/TREX complex, namely TREX COMPONENT 1 (TEX1 or THO3) and HYPER RECOMBINATION1 (HPR1 or THO1) contribute to the correct transcription termination at several loci in Arabidopsis thaliana. We first demonstrate this by showing defective termination in tex1 and hpr1 mutants at the nopaline synthase (NOS) terminator present in a T-DNA inserted between exon 1 and 3 of the PHO1 locus in the pho1-7 mutant. Read-through transcription beyond the NOS terminator and splicing-out of the T-DNA resulted in the generation of a near full-length PHO1 mRNA (minus exon 2) in the tex1 pho1-7 and hpr1 pho1-7 double mutants, with enhanced production of a truncated PHO1 protein that retained phosphate export activity. Consequently, the strong reduction of shoot growth associated with the severe phosphate deficiency of the pho1-7 mutant was alleviated in the tex1 pho1-7 and hpr1 pho1-7 double mutants. Additionally, we show that RNA termination defects in tex1 and hpr1 mutants leads to 3’UTR extensions in several endogenous genes. These results demonstrate that THO/TREX complex contributes to the regulation of transcription termination.

Production of messenger RNAs (mRNAs) involves numerous steps including initiation of transcription, elongation, splicing, termination, as well as export out of the nucleus. All these steps are highly coordinated and failure in any steps has a profound impact on the level and identity of mRNAs produced. The THO/TREX protein complex is associated with nascent RNAs and contributes to several mRNA biogenesis steps, including splicing and export. However, the contribution of the THO/TREX complex to mRNA termination was poorly defined. We have identified a role for two components of the THO/TREX complex, namely the proteins TEX1 and HPR1, in the control of transcription termination in the plant Arabidopsis thaliana. We show that the tex1 and hpr1 mutants have defects in terminating mRNA at the nopaline synthase (NOS) terminator found in a T-DNA insertion mutant leading to the transcriptional read-through pass the NOS terminator. We also show that tex1 and hpr1 mutants have defects in mRNA termination at several endogenous genes, leading to the production of 3’UTR extensions. Together, these results highlight a role for the THO/TREX complex in mRNA termination.

Characterization and Functional Analysis of Polyadenylation Sites in Fast and Slow Muscles

Many increasing documents have proved that alternative polyadenylation (APA) events with different polyadenylation sites (PAS) contribute to posttranscriptional regulation. However, little is known about the detailed molecular features of PASs and its role in porcine fast and slow skeletal muscles through microRNAs (miRNAs) and RNA binding proteins (RBPs). In this study, we combined single-molecule real-time sequencing and Illumina RNA-seq datasets to comprehensively analyze polyadenylation in pigs. We identified a total of 10,334 PASs, of which 8734 were characterized by reference genome annotation. 32.86% of PAS-associated genes were determined to have more than one PAS. Further analysis demonstrated that tissue-specific PASs between fast and slow muscles were enriched in skeletal muscle development pathways. In addition, we obtained 1407 target genes regulated by APA events through potential binding 69 miRNAs and 28 RBPs in variable 3′ UTR regions and some are involved in myofiber transformation. Furthermore, the de novo motif search confirmed that the most common usage of canonical motif AAUAAA and three types of PASs may be related to the strength of motifs. In summary, our results provide a useful annotation of PASs for pig transcriptome and suggest that APA may serve as a role in fast and slow muscle development under the regulation of miRNAs and RBPs.