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

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

Differential alternative polyadenylation of homoeologous genes of allohexaploid wheat ABD subgenomes during drought stress response

Because allohexaploid wheat genome contains ABD subgenomes, how the expression of homoeologous genes is coordinated remains largely unknown, particularly at the co-transcriptional level. Alternative polyadenylation (APA) is an important part of co-transcriptional regulation, which is crucial in developmental processes and stress responses. Drought stress is a major threat to the stable yield of wheat. Focusing on APA, we used poly(A) tag sequencing to track poly(A) site dynamics in wheat under drought stress. The results showed that drought stress led to extensive APA involving 37-47% of differentially expressed genes in wheat. Significant poly(A) site switching was found in stress-responsive genes. Interestingly, homoeologous genes exhibit unequal numbers of poly(A) sites, divergent APA patterns with tissue specificity and time-course dynamics, and distinct 3'-UTR length changes. Moreover, differentially expressed transcripts in leaves and roots used different poly(A) signals, the up- and downregulated isoforms had distinct preferences for non-canonical poly(A) sites. Genes that encode key polyadenylation factors showed differential expression patterns under drought stress. In summary, poly(A) signals and the changes in core poly(A) factors may widely affect the selection of poly(A) sites and gene expression levels during the response to drought stress, and divergent APA patterns among homoeologous genes add extensive plasticity to this responsive network. These results not only reveal the significant role of APA in drought stress response, but also provide a fresh perspective on how homoeologous genes contribute to adaptability through transcriptome diversity. In addition, this work provides information about the ends of transcripts for a better annotation of the wheat genome.

PABPN1 regulates mRNA alternative polyadenylation to inhibit bladder cancer progression

BACKGROUND

About 10-20% of patients with bladder cancer (BC) progress to muscle-invasive diseases, of which the underlying key molecular events have yet to be addressed.

RESULTS

Here, we identified poly(A) binding protein nuclear 1 (PABPN1), a general factor of alternative polyadenylation (APA), was downregulated in BC. Overexpression and knockdown of PABPN1 significantly decreased and increased BC aggressiveness, respectively. Mechanistically, we provide evidence that the preference of PABPN1-bound polyadenylation signals (PASs) depends on the relative location between canonical and non-canonical PASs. PABPN1 shapes inputs converging on Wnt signaling, cell cycle, and lipid biosynthesis.

CONCLUSIONS

Together, these findings provide insights into how PABPN1-mediated APA regulation contributes to BC progression, and suggest that pharmacological targeting PABPN1 might have therapeutic potential in patients with BC.

Transcript variants of long-chain acyl-CoA synthase 1 have distinct roles in sheep lipid metabolism

Mutton has recently been identified to be a consumer favorite, and intermuscular fat is the key factor in determining meat tenderness. Long-chain acyl-CoA synthetase 1 (ACSL1) is a vital subtype of the ACSL family that is involved in the synthesis of lipids from acyl-CoA and the oxidation of fatty acids. The amplification of the ACSL1 gene using rapid amplification of cDNA ends revealed that the alternative polyadenylation (APA) results in two transcripts of the ACSL1 gene. Exon 18 had premature termination, resulting in a shorter CDS region. In this study, the existence of two transcripts of varying lengths translated normally and designated ACSL1-a and ACSL1-b was confirmed. Overexpression of ACSL1-a can promote the synthesis of an intracellular diglyceride, while ACSL1-b can promote triglyceride synthesis. The transfection of ACSL1 shRNA knocks down both the transcripts, the triglyceride content was significantly reduced after differentiation and induction; and lipidome sequencing results exhibited a significant decrease in 14-22 carbon triglyceride metabolites. The results of the present study indicated that the ACSL1 gene played a crucial role in the synthesis of triglycerides. Furthermore, the two transcripts involved in various interactions in the triglyceride synthesis process may be the topic of interest for future research and provide a more theoretical basis for sheep breeding.

FNDC3B 3'-UTR shortening escapes from microRNA-mediated gene repression and promotes nasopharyngeal carcinoma progression

Alternative polyadenylation (APA), which induces shortening of the 3'-UTR, is emerging as an important feature in cancer development and progression. Nevertheless, the effects and mechanisms of APA-induced 3'-UTR shortening in nasopharyngeal carcinoma (NPC) remain largely unclear. Fibronectin type III domain containing 3B (FNDC3B) tended to use proximal polyadenylation site and produce shorter 3'-UTR according to our previous sequencing study. Herein, we found that FNDC3B with shorter 3'-UTR could escape from miRNA-mediated gene repression, and caused its increased expression in NPC. Knocking down of FNDC3B inhibited NPC cell proliferation, migration, invasion, and metastasis in vitro and in vivo. Overexpression of FNDC3B, especially those with shorter 3'-UTR, promoted NPC progression. Furthermore, the mechanism study revealed that FNDC3B could bind to and stabilize myosin heavy chain 9 (MYH9) to activate the Wnt/β-catenin signaling pathway. In addition, MYH9 could reverse the inhibitory effects of FNDC3B knockdown in NPC. Altogether, our results suggested that the 3'-UTR shortening of FNDC3B mRNA mediated its overexpression in NPC and promoted NPC progression by targeting MYH9. This newly identified FNDC3B-MYH9-Wnt/β-catenin axis could represent potential targets for individualized treatment in NPC.

A comprehensive analysis of core polyadenylation sequences and regulation by microRNAs in a set of cancer predisposition genes

Two core polyadenylation elements (CPE) located in the 3' untranslated region of eukaryotic pre-mRNAs play an essential role in their processing: the polyadenylation signal (PAS) AAUAAA and the cleavage site (CS), preferentially a CA dinucleotide. Herein, we characterized PAS and CS sequences in a set of cancer predisposition genes (CPGs) and performed an in silico investigation of microRNAs (miRNAs) regulation to identify potential tumor-suppressive and oncogenic miRNAs. NCBI and alternative polyadenylation databases were queried to characterize CPE sequences in 117 CPGs, including 81 and 17 known tumor suppressor genes and oncogenes, respectively. miRNA-mediated regulation analysis was performed using predicted and validated data sources. Based on NCBI analyses, we did not find an established PAS in 21 CPGs, and verified that the majority of PAS already described (74.4%) had the canonical sequence AAUAAA. Interestingly, "AA" dinucleotide was the most common CS (37.5%) associated with this set of genes. Approximately 90% of CPGs exhibited evidence of alternative polyadenylation (more than one functional PAS). Finally, the mir-192 family was significantly overrepresented as regulator of tumor suppressor genes (P < 0.01), which suggests a potential oncogenic function. Overall, this study provides a landscape of CPE in CPGs, which might be useful in development of future molecular analyses covering these frequently neglected regulatory sequences.

Differential alternative polyadenylation contributes to the developmental divergence between two rice subspecies, japonica and indica

Alternative polyadenylation (APA) is a widespread post-transcriptional mechanism that regulates gene expression through mRNA metabolism, playing a pivotal role in modulating phenotypic traits in rice (Oryza sativa L.). However, little is known about the APA-mediated regulation underlying the distinct characteristics between two major rice subspecies, indica and japonica. Using a poly(A)-tag sequencing approach, polyadenylation (poly(A)) site profiles were investigated and compared pairwise from germination to the mature stage between indica and japonica, and extensive differentiation in APA profiles was detected genome-wide. Genes with subspecies-specific poly(A) sites were found to contribute to subspecies characteristics, particularly in disease resistance of indica and cold-stress tolerance of japonica. In most tissues, differential usage of APA sites exhibited an apparent impact on the gene expression profiles between subspecies, and genes with those APA sites were significantly enriched in quantitative trait loci (QTL) related to yield traits, such as spikelet number and 1000-seed weight. In leaves of the booting stage, APA site-switching genes displayed global shortening of 3' untranslated regions with increased expression in indica compared with japonica, and they were overrepresented in the porphyrin and chlorophyll metabolism pathways. This phenomenon may lead to a higher chlorophyll content and photosynthesis in indica than in japonica, being associated with their differential growth rates and yield potentials. We further constructed an online resource for querying and visualizing the poly(A) atlas in these two rice subspecies. Our results suggest that APA may be largely involved in developmental differentiations between two rice subspecies, especially in leaf characteristics and the stress response, broadening our knowledge of the post-transcriptional genetic basis underlying the divergence of rice traits.

Cleavage factor 25 deregulation contributes to pulmonary fibrosis through alternative polyadenylation

Idiopathic pulmonary fibrosis (IPF) is a chronic and deadly disease with a poor prognosis and few treatment options. Pathological remodeling of the extracellular matrix (ECM) by myofibroblasts is a key factor that drives disease pathogenesis, although the underlying mechanisms remain unknown. Alternative polyadenylation (APA) has recently been shown to play a major role in cellular responses to stress by driving the expression of fibrotic factors and ECMs through altering microRNA sensitivity, but a connection to IPF has not been established. Here, we demonstrate that CFIm25, a global regulator of APA, is down-regulated in the lungs of patients with IPF and mice with pulmonary fibrosis, with its expression selectively reduced in alpha-smooth muscle actin (α-SMA) positive fibroblasts. Following the knockdown of CFIm25 in normal human lung fibroblasts, we identified 808 genes with shortened 3'UTRs, including those involved in the transforming growth factor-β signaling pathway, the Wnt signaling pathway, and cancer pathways. The expression of key pro-fibrotic factors can be suppressed by CFIm25 overexpression in IPF fibroblasts. Finally, we demonstrate that deletion of CFIm25 in fibroblasts or myofibroblast precursors using either the Col1a1 or the Foxd1 promoter enhances pulmonary fibrosis after bleomycin exposure in mice. Taken together, our results identified CFIm25 down-regulation as a novel mechanism to elevate pro-fibrotic gene expression in pulmonary fibrosis.

Differential genome-wide profiling of alternative polyadenylation sites in nasopharyngeal carcinoma by high-throughput sequencing

Background

Alternative polyadenylation (APA) is a widespread phenomenon in the posttranscriptional regulation of gene expression that generates mRNAs with alternative 3′-untranslated regions (3’UTRs). APA contributes to the pathogenesis of various diseases, including cancer. However, the potential role of APA in the development of nasopharyngeal carcinoma (NPC) remains largely unknown.

Methods

A strategy of sequencing APA sites (SAPAS) based on second-generation sequencing technology was carried out to explore the global patterns of APA sites and identify genes with tandem 3’UTRs in samples from 6 NPC and 6 normal nasopharyngeal epithelial tissue (NNET). Sequencing results were then validated using quantitative RT-PCR in a larger cohort of 16 NPC and 16 NNET samples.

Results

The sequencing data showed that the use of tandem APA sites was prevalent in NPC, and numerous genes with APA-switching events were discovered. In total, we identified 195 genes with significant differences in the tandem 3’UTR length between NPC and NNET; including 119 genes switching to distal poly (A) sites and 76 genes switching to proximal poly (A) sites. Several gene ontology (GO) terms were enriched in the list of genes with switched APA sites, including regulation of cell migration, macromolecule catabolic process, protein catabolic process, proteolysis, small conjugating protein ligase activity, and ubiquitin-protein ligase activity.

Conclusions

APA site-switching events are prevalent in NPC. APA-mediated regulation of gene expression may play an important role in the development of NPC, and more detailed studies targeting genes with APA-switching events may contribute to the development of novel future therapeutic strategies for NPC.

Electronic supplementary material

The online version of this article (10.1186/s12929-018-0477-6) contains supplementary material, which is available to authorized users.

Maximizing the Utility of Cancer Transcriptomic Data

Transcriptomic profiling has been applied to large numbers of cancer samples, by large-scale consortia, including The Cancer Genome Atlas, International Cancer Genome Consortium, and Cancer Cell Line Encyclopedia. Advances in mining cancer transcriptomic data enable us to understand the endless complexity of the cancer transcriptome and thereby to discover new biomarkers and therapeutic targets. In this paper, we review computational resources for deep mining of transcriptomic data to identify, quantify, and determine the functional effects and clinical utility of transcriptomic events, including noncoding RNAs, post-transcriptional regulation, exogenous RNAs, and transcribed genetic variants. These approaches can be applied to other complex diseases, thereby greatly leveraging the impact of this work.

Tandem 3' UTR Patterns and Gene Expression Profiles of Marc-145 Cells During PRRSV Infection

Porcine reproductive and respiratory syndrome virus (PRRSV) causes substantial economic losses to the global pig industry. Alternative polyadenylation (APA) is a mechanism that diversifies gene expression, which is important for tumorigenesis, development, and cell differentiation. However, it is unclear whether APA plays a role in the course of PRRSV infection. To address this issue, in this study we carried out a whole-genome transcriptome analysis of PRRSV-infected Marc-145 African green monkey kidney cells and identified 185 APA switching genes and 393 differentially expressed genes (DEGs). Most of these genes were involved in cellular process, metabolism, and biological regulation, and there was some overlap between the two gene sets. DEGs were found to be more directly involved in the antiviral response than APA genes. These findings provide insight into the dynamics of host gene regulation during PRRSV infection and a basis for elucidating the pathogenesis of PRRSV.

RNA Binding Proteins and Regulation of mRNA Translation in Erythropoiesis

Control of gene expression in erythropoiesis has to respond to signals that may emerge from intracellular processes or environmental factors. Control of mRNA translation allows for relatively rapid modulation of protein synthesis from the existing transcriptome. For instance, the protein synthesis rate needs to be reduced when reactive oxygen species or unfolded proteins accumulate in the cells, but also when iron supply is low or when growth factors are lacking in the environment. In addition, regulation of mRNA translation can be important as an additional layer of control on top of gene transcription, in which RNA binding proteins (RBPs) can modify translation of a set of transcripts to the cell’s actual protein requirement. The 5′ and 3′ untranslated regions of mRNA (5′UTR, 3′UTR) contain binding sites for general and sequence specific translation factors. They also contain secondary structures that may hamper scanning of the 5′UTR by translation complexes or may help to recruit translation factors. In addition, the term 5′UTR is not fully correct because many transcripts contain small open reading frames in their 5′UTR that are translated and contribute to regulation of mRNA translation. It is becoming increasingly clear that the transcriptome only partly predicts the proteome. The aim of this review is (i) to summarize how the availability of general translation initiation factors can selectively regulate transcripts because the 5′UTR contains secondary structures or short translated sequences, (ii) to discuss mechanisms that control the length of the mRNA poly(A) tail in relation to mRNA translation, and (iii) to give examples of sequence specific RBPs and their targets. We focused on transcripts and RBPs required for erythropoiesis. Whereas differentiation of erythroblasts to erythrocytes is orchestrated by erythroid transcription factors, the production of erythrocytes needs to respond to the availability of growth factors and nutrients, particularly the availability of iron.

MicroRNA-143 Targets ATG2B to Inhibit Autophagy and Increase Inflammatory Responses in Crohn's Disease

BACKGROUND

Dysfunctional autophagy is recognized as a contributing factor in many chronic inflammatory diseases, including Crohn's disease (CD). Genetic analyses have found that microRNA (miRNA) levels are altered in the intestinal tissues of CD patients.

METHODS

The Sequencing Alternative Poly-Adenylation Sites (SAPAS) method was used to compare the 3' end of the total mRNA sequence of 3 surgical specimens of CD patients (including inflamed tissues and corresponding noninflamed tissues in each case). The levels of autophagy-related 2B (ATG2B), LC3, and miR-143 were compared between inflamed tissues and noninflamed tissues using immunoblot and quantitative reverse transcription polymerase chain reaction. Luciferase assays were used to verify the interactions between miR-143 and ATG2B. Autophagy was measured by immunoblot analyses of LC3 and transmission electron microscopy. Inflammatory cytokines and IκBα were analyzed to evaluate the effect of miR-143 on inflammatory response.

RESULTS

The tandem repeat 3'-UTR of ATG2B was longer in inflamed tissues than in corresponding noninflamed tissues and contained an miR-143 target site. miR-143 expression was elevated, whereas ATG2B and LC3-II were downregulated in inflamed tissues. The direct interaction between miR-143 and ATG2B was verified by a 3'-UTR dual-luciferase reporter assay. Constitutive expression of miR-143 or depletion of ATG2B in cultured intestinal epithelial cells inhibited autophagy, reduced IκBα levels, and increased inflammatory responses.

CONCLUSIONS

miR-143 may induce bowel inflammation by regulating ATG2B and autophagy, suggesting that miR-143 might play a critical role in the development of CD. Therefore, miR-143 could be a promising novel target for gene therapy in CD patients.

Defining the 5΄ and 3΄ landscape of the Drosophila transcriptome with Exo-seq and RNaseH-seq

Cells regulate biological responses in part through changes in transcription start sites (TSS) or cleavage and polyadenylation sites (PAS). To fully understand gene regulatory networks, it is therefore critical to accurately annotate cell type-specific TSS and PAS. Here we present a simple and straightforward approach for genome-wide annotation of 5΄- and 3΄-RNA ends. Our approach reliably discerns bona fide PAS from false PAS that arise due to internal poly(A) tracts, a common problem with current PAS annotation methods. We applied our methodology to study the impact of temperature on the Drosophila melanogaster head transcriptome. We found hundreds of previously unidentified TSS and PAS which revealed two interesting phenomena: first, genes with multiple PASs tend to harbor a motif near the most proximal PAS, which likely represents a new cleavage and polyadenylation signal. Second, motif analysis of promoters of genes affected by temperature suggested that boundary element association factor of 32 kDa (BEAF-32) and DREF mediates a transcriptional program at warm temperatures, a result we validated in a fly line where beaf-32 is downregulated. These results demonstrate the utility of a high-throughput platform for complete experimental and computational analysis of mRNA-ends to improve gene annotation.

Dysregulation of Alternative Poly-adenylation as a Potential Player in Autism Spectrum Disorder

We present here the hypothesis that alternative poly-adenylation (APA) is dysregulated in the brains of individuals affected by Autism Spectrum Disorder (ASD), due to disruptions in the calcium signaling networks. APA, the process of selecting different poly-adenylation sites on the same gene, yielding transcripts with different-length 3′ untranslated regions (UTRs), has been documented in different tissues, stages of development and pathologic conditions. Differential use of poly-adenylation sites has been shown to regulate the function, stability, localization and translation efficiency of target RNAs. However, the role of APA remains rather unexplored in neurodevelopmental conditions. In the human brain, where transcripts have the longest 3′ UTRs and are thus likely to be under more complex post-transcriptional regulation, erratic APA could be particularly detrimental. In the context of ASD, a condition that affects individuals in markedly different ways and whose symptoms exhibit a spectrum of severity, APA dysregulation could be amplified or dampened depending on the individual and the extent of the effect on specific genes would likely vary with genetic and environmental factors. If this hypothesis is correct, dysregulated APA events might be responsible for certain aspects of the phenotypes associated with ASD. Evidence supporting our hypothesis is derived from standard RNA-seq transcriptomic data but we suggest that future experiments should focus on techniques that probe the actual poly-adenylation site (3′ sequencing). To address issues arising from the use of post-mortem tissue and low numbers of heterogeneous samples affected by confounding factors (such as the age, gender and health of the individuals), carefully controlled in vitro systems will be required to model the effect of calcium signaling dysregulation in the ASD brain.

A novel method for genome-wide profiling of dynamic host-pathogen interactions using 3' end enriched RNA-seq

Marek's disease is a contagious lymphoproliferative disease of chickens and typical model of viral oncogenesis. Mapping changes or different states over the course of infection for both host and pathogen would provide important insights into dynamic host-pathogen interactions. Here we introduced 3' end enriched RNA-seq as a novel method to study host-pathogen interactions in chicken embryo fibroblasts cells challenged with Marek's disease virus. The method allowed accurate profiling of gene expression and alternative polyadenylation sites for host and pathogen simultaneously. We totally identified 476 differentially expressed genes and 437 APA switching genes in host, including switching in tandem 3' UTRs and switching between coding region and 3' UTR. Most of these genes were related to innate immunity, apoptosis and metabolism, but two sets of genes overlapped a little, suggesting two complementary mechanisms in gene regulation during MDV infection. In summary, our results provided a relatively comprehensive insight into dynamic host-pathogen interactions in regulation of gene transcription during infection of Marek's disease virus and suggested that 3' end enriched RNA-seq was a promising method to investigate global host-pathogen interactions.

Cleavage factor Im (CFIm) as a regulator of alternative polyadenylation

Most mammalian protein coding genes are subject to alternative cleavage and polyadenylation (APA), which can generate distinct mRNA 3'UTRs with differing regulatory potential. Although this process has been intensely studied in recent years, it remains unclear how and to what extent cleavage site selection is regulated under different physiological conditions. The cleavage factor Im (CFIm) complex is a core component of the mammalian cleavage machinery, and the observation that its depletion causes transcriptome-wide changes in cleavage site use makes it a key candidate regulator of APA. This review aims to summarize current knowledge of the CFIm complex, and explores the evidence surrounding its potential contribution to regulation of APA.

Genome-Wide Analysis of Polyadenylation Events in Schmidtea mediterranea

In eukaryotes, 3' untranslated regions (UTRs) play important roles in regulating posttranscriptional gene expression. The 3'UTR is defined by regulated cleavage/polyadenylation of the pre-mRNA. The advent of next-generation sequencing technology has now enabled us to identify these events on a genome-wide scale. In this study, we used poly(A)-position profiling by sequencing (3P-Seq) to capture all poly(A) sites across the genome of the freshwater planarian, Schmidtea mediterranea, an ideal model system for exploring the process of regeneration and stem cell function. We identified the 3'UTRs for ∼14,000 transcripts and thus improved the existing gene annotations. We found 97 transcripts, which are polyadenylated within an internal exon, resulting in the shrinking of the ORF and loss of a predicted protein domain. Around 40% of the transcripts in planaria were alternatively polyadenylated (ApA), resulting either in an altered 3'UTR or a change in coding sequence. We identified specific ApA transcript isoforms that were subjected to miRNA mediated gene regulation using degradome sequencing. In this study, we also confirmed a tissue-specific expression pattern for alternate polyadenylated transcripts. The insights from this study highlight the potential role of ApA in regulating the gene expression essential for planarian regeneration.

Alternative Polyadenylation: Another Foe in Cancer

Advancements in sequencing and transcriptome analysis methods have led to seminal discoveries that have begun to unravel the complexity of cancer. These studies are paving the way toward the development of improved diagnostics, prognostic predictions, and targeted treatment options. However, it is clear that pieces of the cancer puzzle are still missing. In an effort to have a more comprehensive understanding of the development and progression of cancer, we have come to appreciate the value of the noncoding regions of our genomes, partly due to the discovery of miRNAs and their significance in gene regulation. Interestingly, the miRNA-mRNA interactions are not solely dependent on variations in miRNA levels. Instead, the majority of genes harbor multiple polyadenylation signals on their 3' UTRs (untranslated regions) that can be differentially selected on the basis of the physiologic state of cells, resulting in alternative 3' UTR isoforms. Deregulation of alternative polyadenylation (APA) has increasing interest in cancer research, because APA generates mRNA 3' UTR isoforms with potentially different stabilities, subcellular localizations, translation efficiencies, and functions. This review focuses on the link between APA and cancer and discusses the mechanisms as well as the tools available for investigating APA events in cancer. Overall, detection of deregulated APA-generated isoforms in cancer may implicate some proto-oncogene activation cases of unknown causes and may help the discovery of novel cases; thus, contributing to a better understanding of molecular mechanisms of cancer. Mol Cancer Res; 14(6); 507-17. ©2016 AACR.

Dynamic Regulation of Tandem 3' Untranslated Regions in Zebrafish Spleen Cells during Immune Response

Alternative polyadenylation (APA) has been found to be involved in tumorigenesis, development, and cell differentiation, as well as in the activation of several subsets of immune cells in vitro. Whether APA takes place in immune responses in vivo is largely unknown. We profiled the variation in tandem 3' untranslated regions (UTRs) in pathogen-challenged zebrafish and identified hundreds of APA genes with ∼ 10% being immune response genes. The detected immune response APA genes were enriched in TLR signaling, apoptosis, and JAK-STAT signaling pathways. A greater number of microRNA target sites and AU-rich elements were found in the extended 3' UTRs than in the common 3' UTRs of these APA genes. Further analysis suggested that microRNA and AU-rich element-mediated posttranscriptional regulation plays an important role in modulating the expression of APA genes. These results indicate that APA is extensively involved in immune responses in vivo, and it may be a potential new paradigm for immune regulation.

Expression profiling of the GBP1 gene as a candidate gene for porcine reproductive and respiratory syndrome resistance

A genomic region in pig chromosome 4 has been previously associated with higher viraemia levels and lower weight gain following porcine reproduction and respiratory syndrome virus (PRRSV) infection. The region includes the marker WUR1000125, a G>A polymorphism next to a putative polyadenylation site in the 3'-untranslated region (3'-UTR) of the guanylate-binding protein 1, interferon-induced (GBP1) gene. The protein encoded by GBP1 is a negative regulator of T-cell responses. We show here that GBP1 expression is lower in liver and tonsils of pigs carrying the WUR1000125-G allele due to differential allele expression (allele A expression is 1.9-fold higher than for allele G). We also show that the GBP1 gene has two active polyadenylation signals 421 bp apart and that polyadenylation usage is dependent on the WUR1000125 genotype. The distal site is the most prevalently used in all samples, but the presence of the A allele favours the generation of shorter transcripts from the proximal site. This is confirmed by a differential allele expression study in AG genotype liver and tonsil samples. The interaction between WUR1000125 and other mutations identified in the 5'- and 3'-UTR regions of this gene needs to be studied. In conclusion, our study indicates that the WUR1000125 mutation is associated with changes in the expression of the negative T-cell regulator GBP1 gene. However, the chromosome 4 locus for PRRSV viraemia levels and weight gain contains a cluster of four other GBP genes that remain to be studied as candidate genes for this QTL.

Evaluation of two statistical methods provides insights into the complex patterns of alternative polyadenylation site switching

Switching between different alternative polyadenylation (APA) sites plays an important role in the fine tuning of gene expression. New technologies for the execution of 3’-end enriched RNA-seq allow genome-wide detection of the genes that exhibit significant APA site switching between different samples. Here, we show that the independence test gives better results than the linear trend test in detecting APA site-switching events. Further examination suggests that the discrepancy between these two statistical methods arises from complex APA site-switching events that cannot be represented by a simple change of average 3’-UTR length. In theory, the linear trend test is only effective in detecting these simple changes. We classify the switching events into four switching patterns: two simple patterns (3’-UTR shortening and lengthening) and two complex patterns. By comparing the results of the two statistical methods, we show that complex patterns account for 1/4 of all observed switching events that happen between normal and cancerous human breast cell lines. Because simple and complex switching patterns may convey different biological meanings, they merit separate study. We therefore propose to combine both the independence test and the linear trend test in practice. First, the independence test should be used to detect APA site switching; second, the linear trend test should be invoked to identify simple switching events; and third, those complex switching events that pass independence testing but fail linear trend testing can be identified.

Advanced Datamining Using RNAseq Data

In recent years, RNA sequencing (RNAseq) technology has experienced a rapid rise in popularity. Often seen as a competitor of and the ultimate successor to microarray technology given its more accurate and quantitative gene expression measurement, RNAseq also offers a wealth of additional information that is often overlooked, and given the massive accumulation of RNAseq data available in public data repositories over the past few years, these data are ripe for discovery. Abundant opportunities exist for researchers to conduct in-depth, non-traditional analyses that take advantage of these secondary uses and for bioinformaticians to develop tools to make these data more accessible. This is discussed in this chapter.

Dynamic analyses of alternative polyadenylation from RNA-seq reveal a 3'-UTR landscape across seven tumour types

Alternative polyadenylation (APA) is a pervasive mechanism in the regulation of most human genes, and its implication in diseases including cancer is only beginning to be appreciated. Since conventional APA profiling has not been widely adopted, global cancer APA studies are very limited. Here we develop a novel bioinformatics algorithm (DaPars) for the de novo identification of dynamic APAs from standard RNA-seq. When applied to 358 TCGA Pan-Cancer tumour/normal pairs across seven tumour types, DaPars reveals 1,346 genes with recurrent and tumour-specific APAs. Most APA genes (91%) have shorter 3'-untranslated regions (3' UTRs) in tumours that can avoid microRNA-mediated repression, including glutaminase (GLS), a key metabolic enzyme for tumour proliferation. Interestingly, selected APA events add strong prognostic power beyond common clinical and molecular variables, suggesting their potential as novel prognostic biomarkers. Finally, our results implicate CstF64, an essential polyadenylation factor, as a master regulator of 3'-UTR shortening across multiple tumour types.

APASdb: a database describing alternative poly(A) sites and selection of heterogeneous cleavage sites downstream of poly(A) signals

Increasing amounts of genes have been shown to utilize alternative polyadenylation (APA) 3′-processing sites depending on the cell and tissue type and/or physiological and pathological conditions at the time of processing, and the construction of genome-wide database regarding APA is urgently needed for better understanding poly(A) site selection and APA-directed gene expression regulation for a given biology. Here we present a web-accessible database, named APASdb (http://mosas.sysu.edu.cn/utr), which can visualize the precise map and usage quantification of different APA isoforms for all genes. The datasets are deeply profiled by the sequencing alternative polyadenylation sites (SAPAS) method capable of high-throughput sequencing 3′-ends of polyadenylated transcripts. Thus, APASdb details all the heterogeneous cleavage sites downstream of poly(A) signals, and maintains near complete coverage for APA sites, much better than the previous databases using conventional methods. Furthermore, APASdb provides the quantification of a given APA variant among transcripts with different APA sites by computing their corresponding normalized-reads, making our database more useful. In addition, APASdb supports URL-based retrieval, browsing and display of exon-intron structure, poly(A) signals, poly(A) sites location and usage reads, and 3′-untranslated regions (3′-UTRs). Currently, APASdb involves APA in various biological processes and diseases in human, mouse and zebrafish.

Profiling of alternative polyadenylation sites in luminal B breast cancer using the SAPAS method

Breast cancer (BC) is a leading cause of cancer-related mortality in females and is recognized as a molecularly heterogeneous disease. Previous studies have suggested that alternative messenger RNA (mRNA) processing, particularly alternative polyadenylation [poly(A)] (APA), can be a powerful molecular biomarker with prognostic potential. Therefore, in the present study, we profiled APA sites in the luminal B subtype of BC by sequencing APA sites (SAPAS) method, in order to assess the relation of these APA site-switching events to the recognized molecular subtypes of BC, and to discover novel candidate genes and pathways in BC. Through comprehensive analysis, the trend of APA site-switching events in the 3' untranslated regions (3'UTRs) in the luminal B subtype of BC were found to be the same as that in MCF7 cell lines. Among the genes involved in the events, a significantly greater number of genes was found with shortened 3'UTRs in the samples, which were samples of primary cancer with relatively low proliferation. These findings may provide novel information for the clinical diagnosis and prognosis on a molecular level. Several potential biomarkers with significantly differential tandem 3'UTRs and expression were found and validated. The related biological progresses and pathways involved were partly confirmed by other studies. In conclusion, this study provides new insight into the diagnosis and prognosis of BC from the APA site profile aspect.

Pilot genome-wide study of tandem 3' UTRs in esophageal cancer using high-throughput sequencing

Regulatory regions within the 3' untranslated region (UTR) influence polyadenylation (polyA), translation efficiency, localization and stability of mRNA. Alternative polyA (APA) has been considered to have a key role in gene regulation since 2008. Esophageal carcinoma is the eighth most common type of cancer worldwide. The association between polyA and disease highlights the requirement for comprehensive characterization of genome-wide polyA profiles. In the present study, global polyA profiles were established using the sequencing APA sites (SAPAS) method in order to elucidate the interrelation between 3' UTR length and the development of esophageal cancer. PolyA profiles were analyzed in squamous cell carcinoma, with ~903 genes identified to have shortened 3' UTRs and 917 genes identified to use distal polyA sites. The genes with shortened 3' UTRs were primarily associated with adherens junctions and the cell cycle. Four differentially expressed genes were also found, among which three genes were observed to be upregulated in cancerous tissue and involved in the positive regulation of cell motion, migration and locomotion. One gene was found to be downregulated in cancerous tissue, and associated with oxidative phosphorylation. These findings suggest that esophagitis may have a key role in the development of esophageal carcinoma. Furthermore, the genes with tandem 3' UTRs and differential expression identified in the present study may have the potential to be used as biomarkers for the diagnosis and prognosis of esophageal cancer.

Alternative polyadenylation and differential expression of Shank mRNAs in the synaptic neuropil

The stability and dynamics of synapses rely on tight regulation of the synaptic proteome. Shank proteins, encoded by the three genes Shank1, Shank2 and Shank3 are scaffold molecules in the postsynaptic density of excitatory neurons that contribute to activity-dependent neuronal signalling. Mutations in the Shank genes are associated with neurological diseases. Using state-of-the-art technologies, we investigated the levels of expression of the Shank family messenger RNAs (mRNAs) within the synaptic neuropil of the rat hippocampus. We detected all three Shank transcripts in the neuropil of CA1 pyramidal neurons. We found Shank1 to be the most abundantly expressed among the three Shank mRNA homologues. We also examined the turnover of Shank mRNAs and predict the half-lives of Shank1, Shank2 and Shank3 mRNAs to be 18-28 h. Using 3'-end sequencing, we identified novel 3' ends for the Shank1 and Shank2 3' untranslated regions (3' UTRs) that may contribute to the diversity of alternative polyadenylation (APA) for the Shank transcripts. Our findings consolidate the view that the Shank molecules play a central role at the postsynaptic density. This study may shed light on synaptopathologies associated with disruption of local protein synthesis, perhaps linked to mutations in mRNA 3' UTRs or inappropriate 3' end processing.

Expression profiling of the sex-related gene Dmrt1 in adults of the Lusitanian toadfish Halobatrachus didactylus (Bloch and Schneider, 1801)

Doublesex and mab-3 related transcription factor 1 (Dmrt1) gene is a widely conserved gene involved in sex determination and differentiation across phyla. To gain insights on Dmrt1 implication for fish gonad cell differentiation and gametogenesis development, its mRNA was isolated from testis and ovary from the Lusitanian toadfish (Halobatrachus didactylus). The cDNA from Dmrt1 was synthesized and cloned, whereas its quantitative and qualitative gene expression, as well as its protein immunolocalization, were analyzed. A main product of 1.38 kb, which encodes a protein of 295 aa, was reported, but other minority Dmrt1 products were also identified by RACE-PCR. This gene is predominantly expressed in testis (about 20 times more than in other organs or tissues), specially in spermatogonia, spermatocytes and spermatids, as well as in somatic Sertoli cells, indicating that Dmrt1 plays an important role in spermatogenesis. Although Dmrt1 transcripts also seem to be involved in oogenesis development, and it cannot be excluded that toadfish Dmrt1 could be functionally involved in other processes not related to sex.

Means to an end: mechanisms of alternative polyadenylation of messenger RNA precursors

Expression of mature messenger RNAs (mRNAs) requires appropriate transcription initiation and termination, as well as pre-mRNA processing by capping, splicing, cleavage, and polyadenylation. A core 3'-end processing complex carries out the cleavage and polyadenylation reactions, but many proteins have been implicated in the selection of polyadenylation sites among the multiple alternatives that eukaryotic genes typically have. In recent years, high-throughput approaches to map both the 3'-end processing sites as well as the binding sites of proteins that are involved in the selection of cleavage sites and in the processing reactions have been developed. Here, we review these approaches as well as the insights into the mechanisms of polyadenylation that emerged from genome-wide studies of polyadenylation across a range of cell types and states.

Alternative cleavage and polyadenylation: extent, regulation and function

The 3' end of most protein-coding genes and long non-coding RNAs is cleaved and polyadenylated. Recent discoveries have revealed that a large proportion of these genes contains more than one polyadenylation site. Therefore, alternative polyadenylation (APA) is a widespread phenomenon, generating mRNAs with alternative 3' ends. APA contributes to the complexity of the transcriptome by generating isoforms that differ either in their coding sequence or in their 3' untranslated regions (UTRs), thereby potentially regulating the function, stability, localization and translation efficiency of target RNAs. Here, we review our current understanding of the polyadenylation process and the latest progress in the identification of APA events, mechanisms that regulate poly(A) site selection, and biological processes and diseases resulting from APA.

Polyadenylation site-specific differences in the activity of the neuronal βCstF-64 protein in PC-12 cells

Recent genome-wide analyses have implicated alternative polyadenylation - the process of regulated mRNA 3' end formation - as a critical mechanism that influences multiple steps of mRNA metabolism in addition to increasing the protein-coding capacity of the genome. Although the functional consequences of alternative polyadenylation are well known, protein factors that regulate this process are poorly characterized. Previously, we described an evolutionarily conserved family of neuronal splice variants of the CstF-64 mRNA, βCstF-64, that we hypothesized to function in alternative polyadenylation in the nervous system. In the present study, we show that βCstF-64 mRNA and protein expression increase in response to nerve growth factor (NGF), concomitant with differentiation of adrenal PC-12 cells into a neuronal phenotype, suggesting a role for βCstF-64 in neuronal gene expression. Using PC-12 cells as model, we show that βCstF-64 is a bona fide polyadenylation protein, as evidenced by its association with the CstF complex, and by its ability to stimulate polyadenylation of luciferase reporter mRNA. Using luciferase assays, we show that βCstF-64 stimulates polyadenylation equivalently at the two weak poly(A) sites of the β-adducin mRNA. Notably, we demonstrate that the activity of βCstF-64 is less than CstF-64 on a strong polyadenylation signal, suggesting polyadenylation site-specific differences in the activity of the βCstF-64 protein. Our data address the polyadenylation functions of βCstF-64 for the first time, and provide initial insights into the mechanism of alternative poly(A) site selection in the nervous system.

High throughput characterizations of poly(A) site choice in plants

The polyadenylation of mRNA in eukaryotes is an important biological process. In recent years, significant progress has been made in the field of mRNA polyadenylation owing to the advent of the next generation DNA sequencing technologies. The high-throughput sequencing capabilities have resulted in the direct experimental determinations of large numbers of polyadenylation sites, analysis of which has revealed a vast potential for the regulation of gene expression in eukaryotes. These collections have been generated using specialized sequencing methods that are targeted to the junction of 3'-UTR and the poly(A) tail. Here we present three variations of such a protocol that has been used for the analysis of alternative polyadenylation in plants. While all these methods use oligo-dT as an anchor to the 3'-end, they differ in the means of generating an anchor for the 5'-end in order to produce PCR products suitable for effective Illumina sequencing; the use of different methods to append 5' adapters expands the possible utility of these approaches. These methods are versatile, reproducible, and may be used for gene expression analysis as well as global determinations of poly(A) site choice.

The ticking tail: daily oscillations in mRNA poly(A) tail length drive circadian cycles in protein synthesis

In this issue of Genes & Development, Kojima and colleagues (pp. 2724-2736) examined the impact of mRNA poly(A) tail length on circadian gene expression. Their study demonstrates how dynamic changes in transcript poly(A) tail length can lead to rhythmic protein expression, irrespective of whether mRNA accumulation is circadian or constitutive.