Differential genome-wide profiling of tandem 3' UTRs among human breast cancer and normal cells by high-throughput sequencing

The polyadenylation code: a unified model for the regulation of mRNA alternative polyadenylation

The majority of eukaryotic genes produce multiple mRNA isoforms with distinct 3' ends through a process called mRNA alternative polyadenylation (APA). Recent studies have demonstrated that APA is dynamically regulated during development and in response to environmental stimuli. A number of mechanisms have been described for APA regulation. In this review, we attempt to integrate all the known mechanisms into a unified model. This model not only explains most of previous results, but also provides testable predictions that will improve our understanding of the mechanistic details of APA regulation. Finally, we briefly discuss the known and putative functions of APA regulation.

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.

Kleat: cleavage site analysis of transcriptomes

In eukaryotic cells, alternative cleavage of 3' untranslated regions (UTRs) can affect transcript stability, transport and translation. For polyadenylated (poly(A)) transcripts, cleavage sites can be characterized with short-read sequencing using specialized library construction methods. However, for large-scale cohort studies as well as for clinical sequencing applications, it is desirable to characterize such events using RNA-seq data, as the latter are already widely applied to identify other relevant information, such as mutations, alternative splicing and chimeric transcripts. Here we describe KLEAT, an analysis tool that uses de novo assembly of RNA-seq data to characterize cleavage sites on 3' UTRs. We demonstrate the performance of KLEAT on three cell line RNA-seq libraries constructed and sequenced by the ENCODE project, and assembled using Trans-ABySS. Validating the KLEAT predictions with matched ENCODE RNA-seq and RNA-PET libraries, we show that the tool has over 90% positive predictive value when there are at least three RNA-seq reads supporting a poly(A) tail and requiring at least three RNA-PET reads mapping within 100 nucleotides as validation. We also compare the performance of KLEAT with other popular RNA-seq analysis pipelines that reconstruct 3' UTR ends, and show that it performs favourably, based on an ROC-like curve.

VAAPA: a web platform for visualization and analysis of alternative polyadenylation

Polyadenylation [poly(A)] is an essential process during the maturation of most mRNAs in eukaryotes. Alternative polyadenylation (APA) as an important layer of gene expression regulation has been increasingly recognized in various species. Here, a web platform for visualization and analysis of alternative polyadenylation (VAAPA) was developed. This platform can visualize the distribution of poly(A) sites and poly(A) clusters of a gene or a section of a chromosome. It can also highlight genes with switched APA sites among different conditions. VAAPA is an easy-to-use web-based tool that provides functions of poly(A) site query, data uploading, downloading, and APA sites visualization. It was designed in a multi-tier architecture and developed based on Smart GWT (Google Web Toolkit) using Java as the development language. VAAPA will be a valuable addition to the community for the comprehensive study of APA, not only by making the high quality poly(A) site data more accessible, but also by providing users with numerous valuable functions for poly(A) site analysis and visualization.

Global and quantitative profiling of polyadenylated RNAs using PAS-seq

mRNA alternative polyadenylation (APA) has been increasingly recognized as a widespread and evolutionarily conserved mechanism for eukaryotic gene regulation. Here we describe a method called poly(A) site sequencing that can not only map RNA polyadenylation sites on a transcriptome level but also provide quantitative information on the relative abundance of polyadenylated RNAs. This method has been successfully used for both global APA analysis and digital gene expression profiling.

Implications of polyadenylation in health and disease

Polyadenylation is the RNA processing step that completes the maturation of nearly all eukaryotic mRNAs. It is a two-step nuclear process that involves an endonucleolytic cleavage of the pre-mRNA at the 3'-end and the polymerization of a polyadenosine (polyA) tail, which is fundamental for mRNA stability, nuclear export and efficient translation during development. The core molecular machinery responsible for the definition of a polyA site includes several recognition, cleavage and polyadenylation factors that identify and act on a given polyA signal present in a pre-mRNA, usually an AAUAAA hexamer or similar sequence. This mechanism is tightly regulated by other cis-acting elements and trans-acting factors, and its misregulation can cause inefficient gene expression and may ultimately lead to disease. The majority of genes generate multiple mRNAs as a result of alternative polyadenylation in the 3'-untranslated region. The variable lengths of the 3' untranslated regions created by alternative polyadenylation are a recognizable target for differential regulation and clearly affect the fate of the transcript, ultimately modulating the expression of the gene. Over the past few years, several studies have highlighted the importance of polyadenylation and alternative polyadenylation in gene expression and their impact in a variety of physiological conditions, as well as in several illnesses. Abnormalities in the 3'-end processing mechanisms thus represent a common feature among many oncological, immunological, neurological and hematological disorders, but slight imbalances can lead to the natural establishment of a specific cellular state. This review addresses the key steps of polyadenylation and alternative polyadenylation in different cellular conditions and diseases focusing on the molecular effectors that ensure a faultless pre-mRNA 3' end formation.

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.

Alternative polyadenylation is associated with lower expression of PABPN1 and poor prognosis in non-small cell lung cancer

Alternative polyadenylation (APA), which induces shortening of the 3′UTR, is emerging as an important phenomenon in gene regulation. APA is involved in development, cancer and cell proliferation. APA may lead to disruption of microRNA-mediated gene silencing in cancer cells via detachment of microRNA binding sites. We studied the correlation between the APA profile and the tumor aggressiveness in cases of lung cancer. We selected the top 10 genes showing significant 3′UTR shortening in lung cancer, using the package of the Bioconductor for probe-level analyses of expression microarrays. We established and evaluated the APA score by quantitative RT-PCR in 147 clinical specimens of non-small cell lung cancer and compared the results with the clinical outcomes and expression levels of APA-related genes, including PABPN1, CPEB1, E2F1 and proliferation markers (MKI67, TOP2A and MCM2). High APA scores were correlated with an advanced tumor stage and a poor prognosis (P < 0.001). Multivariate analysis identified the APA score as an independent prognostic factor (hazard ratio, 3.0; P = 0.03). Both lower expression of PABPN1 and higher expression of the proliferation markers were correlated with high APA scores and a poor prognosis, with suppression of PABPN1 exerting its influence independent of gain of the proliferation markers. Moreover, the APA score was correlated with the maximum standardized uptake value of the tumors on positron emission tomography (r = 0.53; P < 0.001). Our results indicate that the loss of PABPN1, a suppressor of APA, might promote tumor aggressiveness by releasing the cancer cells from microRNA-mediated gene regulation.

mRNA 3'end processing: A tale of the tail reaches the clinic

Recent advances reveal mRNA 3′end processing as a highly regulated process that fine-tunes posttranscriptional gene expression. This process can affect the site and/or the efficiency of 3′end processing, controlling the quality and the quantity of substrate mRNAs. The regulation of 3′end processing plays a central role in fundamental physiology such as blood coagulation and innate immunity. In addition, errors in mRNA 3′end processing have been associated with a broad spectrum of human diseases, including cancer. We summarize and discuss the paradigmatic shift in the understanding of 3′end processing as a mechanism of posttranscriptional gene regulation that has reached clinical medicine.

Introduction to microRNAs in biological systems

MicroRNAs are 20-24-nucleotide-long noncoding RNAs that bind to the 3' UTR (untranslated region) of target mRNAs. Since their discovery, microRNAs have been gaining attention for their ability to contribute to gene expression regulation under various physiological conditions. Consequently, deregulated expression of microRNAs has been linked to different disease states. Here, a brief overview of the canonical and alternative microRNA biogenesis pathways and microRNA functions in biological systems is given based on recent developments. In addition, newly emerging regulatory mechanisms, such as alternative polyadenylation, in connection with microRNA-dependent gene expression regulation are discussed.

3'UTR shortening identifies high-risk cancers with targeted dysregulation of the ceRNA network

Competing endogenous RNA (ceRNA) interactions form a multilayered network that regulates gene expression in various biological pathways. Recent studies have demonstrated novel roles of ceRNA interactions in tumorigenesis, but the dynamics of the ceRNA network in cancer remain unexplored. Here, we examine ceRNA network dynamics in prostate cancer from the perspective of alternative cleavage and polyadenylation (APA) and reveal the principles of such changes. Analysis of exon array data revealed that both shortened and lengthened 3′UTRs are abundant. Consensus clustering with APA data stratified cancers into groups with differing risks of biochemical relapse and revealed that a ceRNA subnetwork enriched with cancer genes was specifically dysregulated in high-risk cancers. The novel connection between 3′UTR shortening and ceRNA network dysregulation was supported by the unusually high number of microRNA response elements (MREs) shared by the dysregulated ceRNA interactions and the significantly altered 3′UTRs. The dysregulation followed a fundamental principle in that ceRNA interactions connecting genes that show opposite trends in expression change are preferentially dysregulated. This targeted dysregulation is responsible for the majority of the observed expression changes in genes with significant ceRNA dysregulation and represents a novel mechanism underlying aberrant oncogenic expression.

RNA sequencing: from tag-based profiling to resolving complete transcript structure

Technological advances in the sequencing field support in-depth characterization of the transcriptome. Here, we review genome-wide RNA sequencing methods used to investigate specific aspects of gene expression and its regulation, from transcription to RNA processing and translation. We discuss tag-based methods for studying transcription, alternative initiation and polyadenylation events, shotgun methods for detection of alternative splicing, full-length RNA sequencing for the determination of complete transcript structures, and targeted methods for studying the process of transcription and translation. With the ensemble of technologies available, it is now possible to obtain a comprehensive view on transcriptome complexity and the regulation of transcript diversity.

A change-point model for identifying 3'UTR switching by next-generation RNA sequencing

MOTIVATION

Next-generation RNA sequencing offers an opportunity to investigate transcriptome in an unprecedented scale. Recent studies have revealed widespread alternative polyadenylation (polyA) in eukaryotes, leading to various mRNA isoforms differing in their 3' untranslated regions (3'UTR), through which, the stability, localization and translation of mRNA can be regulated. However, very few, if any, methods and tools are available for directly analyzing this special alternative RNA processing event. Conventional methods rely on annotation of polyA sites; yet, such knowledge remains incomplete, and identification of polyA sites is still challenging. The goal of this article is to develop methods for detecting 3'UTR switching without any prior knowledge of polyA annotations.

RESULTS

We propose a change-point model based on a likelihood ratio test for detecting 3'UTR switching. We develop a directional testing procedure for identifying dramatic shortening or lengthening events in 3'UTR, while controlling mixed directional false discovery rate at a nominal level. To our knowledge, this is the first approach to analyze 3'UTR switching directly without relying on any polyA annotations. Simulation studies and applications to two real datasets reveal that our proposed method is powerful, accurate and feasible for the analysis of next-generation RNA sequencing data.

CONCLUSIONS

The proposed method will fill a void among alternative RNA processing analysis tools for transcriptome studies. It can help to obtain additional insights from RNA sequencing data by understanding gene regulation mechanisms through the analysis of 3'UTR switching.

AVAILABILITY AND IMPLEMENTATION

The software is implemented in Java and can be freely downloaded from http://utr.sourceforge.net/.

CONTACT

zhiwei@njit.edu or hongzhe@mail.med.upenn.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Tandem alternative polyadenylation events of genes in non-eosinophilic nasal polyp tissue identified by high-throughput sequencing analysis

Nasal polyps (NP) is highly associated with the disorder of immune cells. Alternative polyadenylation (APA) produces mRNA isoforms with different length of 3′-untranslated region (UTR) and regulates gene expression. It has been proven that this APA-mediated regulation of 3′UTR length is an immune-associated phenomenon. The aim of this study was to investigate the genome-wide alternative tandem 3′UTR length switching events in non-eosinophilic nasal polyp tissue. Thirteen patients diagnosed as having non-eosinophilic nasal polyps were included in this study. Nasal polyp tissue and control mucosa were collected during surgery. The 3′ end library of cDNA was constructed. The recovered libraries were sequenced with second sequencing technology, and the sequencing data were analyzed by an in-house bioinformatics pipeline. Tandem 3′UTR length switching between samples was detected by a test of linear trend alternative to independence. We found a significant alteration in the tandem 3′UTR length in 1,920 genes in nasal polyp samples. Functional annotation results showed that several gene ontology (GO) terms were enriched in the list of genes with switched APA sites, including regulation of transcription, macromolecule catabolic localization and mRNA processing. The results suggested that APA-mediated alternative 3′UTR regulation plays an important role in the post-transcriptional regulation of gene expression in non-eosinophilic nasal polyps.

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.

TAIL-seq: genome-wide determination of poly(A) tail length and 3' end modifications

Global investigation of the 3' extremity of mRNA (3'-terminome), despite its importance in gene regulation, has not been feasible due to technical challenges associated with homopolymeric sequences and relative paucity of mRNA. We here develop a method, TAIL-seq, to sequence the very end of mRNA molecules. TAIL-seq allows us to measure poly(A) tail length at the genomic scale. Median poly(A) length is 50-100 nt in HeLa and NIH 3T3 cells. Poly(A) length correlates with mRNA half-life, but not with translational efficiency. Surprisingly, we discover widespread uridylation and guanylation at the downstream of poly(A) tail. The U tails are generally attached to short poly(A) tails (<25 nt), while the G tails are found mainly on longer poly(A) tails (>40 nt), implicating their generic roles in mRNA stability control. TAIL-seq is a potent tool to dissect dynamic control of mRNA turnover and translational control, and to discover unforeseen features of RNA cleavage and tailing.

Simultaneous mapping of transcript ends at single-nucleotide resolution and identification of widespread promoter-associated non-coding RNA governed by TATA elements

Understanding the relationships between regulatory factor binding, chromatin structure, cis-regulatory elements and RNA-regulation mechanisms relies on accurate information about transcription start sites (TSS) and polyadenylation sites (PAS). Although several approaches have identified transcript ends in yeast, limitations of resolution and coverage have remained, and definitive identification of TSS and PAS with single-nucleotide resolution has not yet been achieved. We developed SMORE-seq (simultaneous mapping of RNA ends by sequencing) and used it to simultaneously identify the strongest TSS for 5207 (90%) genes and PAS for 5277 (91%) genes. The new transcript annotations identified by SMORE-seq showed improved distance relationships with TATA-like regulatory elements, nucleosome positions and active RNA polymerase. We found 150 genes whose TSS were downstream of the annotated start codon, and additional analysis of evolutionary conservation and ribosome footprinting suggests that these protein-coding sequences are likely to be mis-annotated. SMORE-seq detected short non-coding RNAs transcribed divergently from more than a thousand promoters in wild-type cells under normal conditions. These divergent non-coding RNAs were less evident at promoters containing canonical TATA boxes, suggesting a model where transcription initiation at promoters by RNAPII is bidirectional, with TATA elements serving to constrain the directionality of initiation.

Multiplex analysis of polyA-linked sequences (MAPS): an RNA-seq strategy to profile poly(A+) RNA

We summarize 12 experimental methods that have been developed for profiling gene expression by focusing on the 3'-end of poly(A+) mRNA, distilling both common and unique features. Of this family of methods, we provide detailed protocol for MAPS, a method we believe is the simplest and most cost-effective for profiling gene expression and quantifying alternative polyadenylation events by oligo-dT priming followed by random priming and deep sequencing. This method also enables library multiplexing by using a set of bar coding primers during PCR amplification. We also provide a general guideline for analysis of the data generated by MAPS by using the software package maps3end.

RNA-binding proteins in regulation of alternative cleavage and polyadenylation

Almost all eukaryotic pre-mRNAs are processed at the 3' end by the cleavage and polyadenylation (C/P) reaction, which preludes termination of transcription and gives rise to the poly(A) tail of mature mRNA. Genomic studies in recent years have indicated that most eukaryotic mRNA genes have multiple cleavage and polyadenylation sites (pAs), leading to alternative cleavage and polyadenylation (APA) products. APA isoforms generally differ in their 3' untranslated regions (3' UTRs), but can also have different coding sequences (CDSs). APA expands the repertoire of transcripts expressed from the genome, and is highly regulated under various physiological and pathological conditions. Growing lines of evidence have shown that RNA-binding proteins (RBPs) play important roles in regulation of APA. Some RBPs are part of the machinery for C/P; others influence pA choice through binding to adjacent regions. In this chapter, we review cis elements and trans factors involved in C/P, the significance of APA, and increasingly elucidated roles of RBPs in APA regulation. We also discuss analysis of APA using transcriptome-wide techniques as well as molecular biology approaches.

Alternative polyadenylation of tumor suppressor genes in small intestinal neuroendocrine tumors

The tumorigenesis of small intestinal neuroendocrine tumors (SI-NETs) is poorly understood. Recent studies have associated alternative polyadenylation (APA) with proliferation, cell transformation, and cancer. Polyadenylation is the process in which the pre-messenger RNA is cleaved at a polyA site and a polyA tail is added. Genes with two or more polyA sites can undergo APA. This produces two or more distinct mRNA isoforms with different 3' untranslated regions. Additionally, APA can also produce mRNAs containing different 3'-terminal coding regions. Therefore, APA alters both the repertoire and the expression level of proteins. Here, we used high-throughput sequencing data to map polyA sites and characterize polyadenylation genome-wide in three SI-NETs and a reference sample. In the tumors, 16 genes showed significant changes of APA pattern, which lead to either the 3' truncation of mRNA coding regions or 3' untranslated regions. Among these, 11 genes had been previously associated with cancer, with 4 genes being known tumor suppressors: DCC, PDZD2, MAGI1, and DACT2. We validated the APA in three out of three cases with quantitative real-time-PCR. Our findings suggest that changes of APA pattern in these 16 genes could be involved in the tumorigenesis of SI-NETs. Furthermore, they also point to APA as a new target for both diagnostic and treatment of SI-NETs. The identified genes with APA specific to the SI-NETs could be further tested as diagnostic markers and drug targets for disease prevention and treatment.

Quantitative polyadenylation site mapping with single-molecule direct RNA sequencing

The known regulatory role of 3' untranslated regions (3'UTRs) and poly(A) tails in RNA localization, stability, and translation, and polyadenylation regulation defects leading to human diseases such as oculopharyngeal muscular dystrophy, thalassemias, thrombophilia, and IPEX syndrome underline the need to fully characterize genome-wide polyadenylation states and mechanisms across normal physiological and disease states. This chapter outlines the quantitative polyadenylation site mapping and analysis strategies developed with the single-molecule direct RNA sequencing technology.

Genome-wide identification and predictive modeling of tissue-specific alternative polyadenylation

Motivation: Pre-mRNA cleavage and polyadenylation are essential steps for 3′-end maturation and subsequent stability and degradation of mRNAs. This process is highly controlled by cis-regulatory elements surrounding the cleavage/polyadenylation sites (polyA sites), which are frequently constrained by sequence content and position. More than 50% of human transcripts have multiple functional polyA sites, and the specific use of alternative polyA sites (APA) results in isoforms with variable 3′-untranslated regions, thus potentially affecting gene regulation. Elucidating the regulatory mechanisms underlying differential polyA preferences in multiple cell types has been hindered both by the lack of suitable data on the precise location of cleavage sites, as well as of appropriate tests for determining APAs with significant differences across multiple libraries.

Results: We applied a tailored paired-end RNA-seq protocol to specifically probe the position of polyA sites in three human adult tissue types. We specified a linear-effects regression model to identify tissue-specific biases indicating regulated APA; the significance of differences between tissue types was assessed by an appropriately designed permutation test. This combination allowed to identify highly specific subsets of APA events in the individual tissue types. Predictive models successfully classified constitutive polyA sites from a biologically relevant background (auROC = 99.6%), as well as tissue-specific regulated sets from each other. We found that the main cis-regulatory elements described for polyadenylation are a strong, and highly informative, hallmark for constitutive sites only. Tissue-specific regulated sites were found to contain other regulatory motifs, with the canonical polyadenylation signal being nearly absent at brain-specific polyA sites. Together, our results contribute to the understanding of the diversity of post-transcriptional gene regulation.

Availability: Raw data are deposited on SRA, accession numbers: brain SRX208132, kidney SRX208087 and liver SRX208134. Processed datasets as well as model code are published on our website: http://www.genome.duke.edu/labs/ohler/research/UTR/

Contact:uwe.ohler@duke.edu

Transcriptome dynamics and diversity in the early zebrafish embryo

Recent years advances in high-throughput sequencing have improved our understanding of how transcripts regulate early vertebrate development. Here, we review the transcriptome dynamics and diversity during early stages of zebrafish embryogenesis. Transcriptome dynamics is characterized by different patterns of mRNA degradation, activation of dormant transcripts and onset of transcription. Several studies have shown a striking diversity of both coding and non-coding transcripts. However, in the aftermath of this immense increase in data, functional studies of both protein-coding and non-coding transcripts are lagging behind. We anticipate that the forthcoming years will see studies relying on different high-throughput sequencing technologies and genomic tools developed for zebrafish embryos to further pin down yet un-annotated transcript-function relationships.

U1 snRNP-Dependent Suppression of Polyadenylation: Physiological Role and Therapeutic Opportunities in Cancer

Pre-mRNA splicing and polyadenylation are critical steps in the maturation of eukaryotic mRNA. U1 snRNP is an essential component of the splicing machinery and participates in splice-site selection and spliceosome assembly by base-pairing to the 5' splice site. U1 snRNP also plays an additional, nonsplicing global function in 3' end mRNA processing; it actively suppresses the polyadenylation machinery from using early, mostly intronic polyadenylation signals which would lead to aberrant, truncated mRNAs. Thus, U1 snRNP safeguards pre-mRNA transcripts against premature polyadenylation and contributes to the regulation of alternative polyadenylation. Here, we review the role of U1 snRNP in 3' end mRNA processing, outline the evidence that led to the recognition of its physiological, general role in inhibiting polyadenylation, and finally highlight the possibility of manipulating this U1 snRNP function for therapeutic purposes in cancer.

Ubiquitously transcribed genes use alternative polyadenylation to achieve tissue-specific expression

A majority of human genes use alternative cleavage and polyadenylation to generate mRNA transcripts that differ in the lengths of their 3′ untranslated regions (UTRs). Here, Lianoglou et al. develop a sequencing method, 3′-seq, to measure 3′ UTR isoform expression across diverse human tissues and isogenic transformation systems. The analyses reveal that during transformation and differentiation, single-UTR genes typically change their mRNA abundance levels, while multi-UTR genes change 3′ UTR isoform ratios to achieve tissue specificity. This study offers surprising new insights into how cell type-specific gene expression is achieved.

More than half of human genes use alternative cleavage and polyadenylation (ApA) to generate mRNA transcripts that differ in the lengths of their 3′ untranslated regions (UTRs), thus altering the post-transcriptional fate of the message and likely the protein output. The extent of 3′ UTR variation across tissues and the functional role of ApA remain poorly understood. We developed a sequencing method called 3′-seq to quantitatively map the 3′ ends of the transcriptome of diverse human tissues and isogenic transformation systems. We found that cell type-specific gene expression is accomplished by two complementary programs. Tissue-restricted genes tend to have single 3′ UTRs, whereas a majority of ubiquitously transcribed genes generate multiple 3′ UTRs. During transformation and differentiation, single-UTR genes change their mRNA abundance levels, while multi-UTR genes mostly change 3′ UTR isoform ratios to achieve tissue specificity. However, both regulation programs target genes that function in the same pathways and processes that characterize the new cell type. Instead of finding global shifts in 3′ UTR length during transformation and differentiation, we identify tissue-specific groups of multi-UTR genes that change their 3′ UTR ratios; these changes in 3′ UTR length are largely independent from changes in mRNA abundance. Finally, tissue-specific usage of ApA sites appears to be a mechanism for changing the landscape targetable by ubiquitously expressed microRNAs.

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.

Role of alternative polyadenylation during adipogenic differentiation: an in silico approach

Post-transcriptional regulation of stem cell differentiation is far from being completely understood. Changes in protein levels are not fully correlated with corresponding changes in mRNAs; the observed differences might be partially explained by post-transcriptional regulation mechanisms, such as alternative polyadenylation. This would involve changes in protein binding, transcript usage, miRNAs and other non-coding RNAs. In the present work we analyzed the distribution of alternative transcripts during adipogenic differentiation and the potential role of miRNAs in post-transcriptional regulation. Our in silico analysis suggests a modest, consistent, bias in 3'UTR lengths during differentiation enabling a fine-tuned transcript regulation via small non-coding RNAs. Including these effects in the analyses partially accounts for the observed discrepancies in relative abundance of protein and mRNA.

3' UTR-isoform choice has limited influence on the stability and translational efficiency of most mRNAs in mouse fibroblasts

Variation in protein output across the genome is controlled at several levels, but the relative contributions of different regulatory mechanisms remain poorly understood. Here, we obtained global measurements of decay and translation rates for mRNAs with alternative 3′ untranslated regions (3′ UTRs) in murine 3T3 cells. Distal tandem isoforms had slightly but significantly lower mRNA stability and greater translational efficiency than proximal isoforms on average. The diversity of alternative 3′ UTRs also enabled inference and evaluation of both positively and negatively acting cis-regulatory elements. The 3′ UTR elements with the greatest implied influence were microRNA complementary sites, which were associated with repression of 32% and 4% at the stability and translational levels, respectively. Nonetheless, both the decay and translation rates were highly correlated for proximal and distal 3′ UTR isoforms from the same genes, implying that in 3T3 cells, alternative 3′ UTR sequences play a surprisingly small regulatory role compared to other mRNA regions.

Distinct polyadenylation landscapes of diverse human tissues revealed by a modified PA-seq strategy

Background

Polyadenylation is a key regulatory step in eukaryotic gene expression and one of the major contributors of transcriptome diversity. Aberrant polyadenylation often associates with expression defects and leads to human diseases.

Results

To better understand global polyadenylation regulation, we have developed a polyadenylation sequencing (PA-seq) approach. By profiling polyadenylation events in 13 human tissues, we found that alternative cleavage and polyadenylation (APA) is prevalent in both protein-coding and noncoding genes. In addition, APA usage, similar to gene expression profiling, exhibits tissue-specific signatures and is sufficient for determining tissue origin. A 3′ untranslated region shortening index (USI) was further developed for genes with tandem APA sites. Strikingly, the results showed that different tissues exhibit distinct patterns of shortening and/or lengthening of 3′ untranslated regions, suggesting the intimate involvement of APA in establishing tissue or cell identity.

Conclusions

This study provides a comprehensive resource to uncover regulated polyadenylation events in human tissues and to characterize the underlying regulatory mechanism.

Accurate identification of polyadenylation sites from 3' end deep sequencing using a naive Bayes classifier

MOTIVATION

3' end processing is important for transcription termination, mRNA stability and regulation of gene expression. To identify 3' ends, most techniques use an oligo-dT primer to construct deep sequencing libraries. However, this approach can lead to identification of artifactual polyadenylation sites due to internal priming in homopolymeric stretches of adenines. Although heuristic filters have been applied in these cases, they typically result in a high proportion of both false-positive and -negative classifications. Therefore, there is a need to develop improved algorithms to better identify mis-priming events in oligo-dT primed sequences.

RESULTS

By analyzing sequence features flanking 3' ends derived from oligo-dT-based sequencing, we developed a naïve Bayes classifier to classify them as true or false/internally primed. The resulting algorithm is highly accurate, outperforms previous heuristic filters and facilitates identification of novel polyadenylation sites.

Comprehensive analysis of microRNA (miRNA) targets in breast cancer cells

MicroRNAs (miRNAs) regulate mRNA stability and translation through the action of the RNAi-induced silencing complex. In this study, we systematically identified endogenous miRNA target genes by using AGO2 immunoprecipitation (AGO2-IP) and microarray analyses in two breast cancer cell lines, MCF7 and MDA-MB-231, representing luminal and basal-like breast cancer, respectively. The expression levels of ∼70% of the AGO2-IP mRNAs were increased by DROSHA or DICER1 knockdown. In addition, integrated analysis of miRNA expression profiles, mRNA-AGO2 interaction, and the 3'-UTR of mRNAs revealed that >60% of the AGO2-IP mRNAs were putative targets of the 50 most abundantly expressed miRNAs. Together, these results suggested that the majority of the AGO2-associated mRNAs were bona fide miRNA targets. Functional enrichment analysis uncovered that the AGO2-IP mRNAs were involved in regulation of cell cycle, apoptosis, adhesion/migration/invasion, stress responses (e.g. DNA damage and endoplasmic reticulum stress and hypoxia), and cell-cell communication (e.g. Notch and Ephrin signaling pathways). A role of miRNAs in regulating cell migration/invasion and stress response was further defined by examining the impact of DROSHA knockdown on cell behaviors. We demonstrated that DROSHA knockdown enhanced cell migration and invasion, whereas it sensitized cells to cell death induced by suspension culture, glucose depletion, and unfolding protein stress. Data from an orthotopic xenograft model showed that DROSHA knockdown resulted in reduced growth of primary tumors but enhanced lung metastasis. Taken together, these results suggest that miRNAs collectively function to promote survival of tumor cells under stress but suppress cell migration/invasion in breast cancer cells.

Impaired pre-mRNA processing and altered architecture of 3' untranslated regions contribute to the development of human disorders

The biological fate of each mRNA and consequently, the protein to be synthesised, is highly dependent on the nature of the 3' untranslated region. Despite its non-coding character, the 3' UTR may affect the final mRNA stability, the localisation, the export from the nucleus and the translation efficiency. The conserved regulatory sequences within 3' UTRs and the specific elements binding to them enable gene expression control at the posttranscriptional level and all these processes reflect the actual state of the cell including proliferation, differentiation, cellular stress or tumourigenesis. Through this article, we briefly outline how the alterations in the establishment and final architecture of 3' UTRs may contribute to the development of various disorders in humans.

Dynamic expression of 3' UTRs revealed by Poisson hidden Markov modeling of RNA-Seq: implications in gene expression profiling

RNA sequencing (RNA-Seq) allows for the identification of novel exon-exon junctions and quantification of gene expression levels. We show that from RNA-Seq data one may also detect utilization of alternative polyadenylation (APA) in 3' untranslated regions (3' UTRs) known to play a critical role in the regulation of mRNA stability, cellular localization and translation efficiency. Given the dynamic nature of APA, it is desirable to examine the APA on a sample by sample basis. We used a Poisson hidden Markov model (PHMM) of RNA-Seq data to identify potential APA in human liver and brain cortex tissues leading to shortened 3' UTRs. Over three hundred transcripts with shortened 3' UTRs were detected with sensitivity >75% and specificity >60%. Tissue-specific 3' UTR shortening was observed for 32 genes with a q-value ≤ 0.1. When compared to alternative isoforms detected by Cufflinks or MISO, our PHMM method agreed on over 100 transcripts with shortened 3' UTRs. Given the increasing usage of RNA-Seq for gene expression profiling, using PHMM to investigate sample-specific 3' UTR shortening could be an added benefit from this emerging technology.

Global patterns of tissue-specific alternative polyadenylation in Drosophila

We analyzed the usage and consequences of alternative cleavage and polyadenylation (APA) in Drosophila melanogaster by using >1 billion reads of stranded mRNA-seq across a variety of dissected tissues. Beyond demonstrating that a majority of fly transcripts are subject to APA, we observed broad trends for 3' untranslated region (UTR) shortening in the testis and lengthening in the central nervous system (CNS); the latter included hundreds of unannotated extensions ranging up to 18 kb. Extensive northern analyses validated the accumulation of full-length neural extended transcripts, and in situ hybridization indicated their spatial restriction to the CNS. Genes encoding RNA binding proteins (RBPs) and transcription factors were preferentially subject to 3' UTR extensions. Motif analysis indicated enrichment of miRNA and RBP sites in the neural extensions, and their termini were enriched in canonical cis elements that promote cleavage and polyadenylation. Altogether, we reveal broad tissue-specific patterns of APA in Drosophila and transcripts with unprecedented 3' UTR length in the nervous system.

Alternative CHRNB4 3'-UTRs mediate the allelic effects of SNP rs1948 on gene expression

Common genetic factors strongly contribute to both nicotine, the main addictive component of tobacco, and alcohol use. Several lines of evidence suggest nicotinic acetylcholine receptors as common sites of action for nicotine and alcohol. Specifically, rs1948, a single-nucleotide polymorphism (SNP) located in the CHRNB4 3′-untranslated region (UTR), has been associated to early age of initiation for both alcohol and tobacco use. To determine the allelic effects of rs1948 on gene expression, two rs1948-containing sequences of different lengths corresponding to the CHRNB4 3′-UTR were cloned into pGL3-promoter luciferase reporter vectors. Data obtained showed that the allelic effects of SNP rs1948 on luciferase expression are mediated by the length and species of transcripts generated. In addition, it was found that miR-3157 increased the overall luciferase expression while miR-138, a microRNA known to play a role in neuroadaptation to drug abuse, decreased luciferase expression when compared to basal conditions. These findings demonstrate the importance of SNP rs1948 on the regulation of CHRNB4 expression and provide the first evidence of CHRNB4 down-regulation by miR-138.

αCP Poly(C) binding proteins act as global regulators of alternative polyadenylation

We have previously demonstrated that the KH-domain protein αCP binds to a 3' untranslated region (3'UTR) C-rich motif of the nascent human alpha-globin (hα-globin) transcript and enhances the efficiency of 3' processing. Here we assess the genome-wide impact of αCP RNA-protein (RNP) complexes on 3' processing with a specific focus on its role in alternative polyadenylation (APA) site utilization. The major isoforms of αCP were acutely depleted from a human hematopoietic cell line, and the impact on mRNA representation and poly(A) site utilization was determined by direct RNA sequencing (DRS). Bioinformatic analysis revealed 357 significant alterations in poly(A) site utilization that could be specifically linked to the αCP depletion. These APA events correlated strongly with the presence of C-rich sequences in close proximity to the impacted poly(A) addition sites. The most significant linkage was the presence of a C-rich motif within a window 30 to 40 bases 5' to poly(A) signals (AAUAAA) that were repressed upon αCP depletion. This linkage is consistent with a general role for αCPs as enhancers of 3' processing. These findings predict a role for αCPs in posttranscriptional control pathways that can alter the coding potential and/or levels of expression of subsets of mRNAs in the mammalian transcriptome.

Alternative cleavage and polyadenylation: the long and short of it

Cleavage and polyadenylation (C/P) of nascent transcripts is essential for maturation of the 3' ends of most eukaryotic mRNAs. Over the past three decades, biochemical studies have elucidated the machinery responsible for the seemingly simple C/P reaction. Recent genomic analyses have indicated that most eukaryotic genes have multiple cleavage and polyadenylation sites (pAs), leading to transcript isoforms with different coding potentials and/or variable 3' untranslated regions (UTRs). As such, alternative cleavage and polyadenylation (APA) is an important layer of gene regulation impacting mRNA metabolism. Here, we review our current understanding of APA and recent progress in this field.

Genome-wide maps of polyadenylation reveal dynamic mRNA 3'-end formation in mammalian cell lineages

Post-transcriptional regulation, often mediated by miRNAs and RNA-binding proteins at the 3' untranslated regions (UTRs) of mRNAs, is implicated in important roles in the output of transcriptome. To decipher this layer of gene regulation, it is essential to measure global mRNA expression quantitatively in a 3'-UTR-specific manner. Here we establish an experimental and bioinformatics pipeline that simultaneously determines 3'-end formation by leveraging local nucleotide composition and quantitatively measures mRNA expression by sequencing polyadenylated transcripts. When applied to purified mouse embryonic skin stem cells and their daughter lineages, we identify 18,060 3' UTRs representing 12,739 distinct mRNAs that are abundantly expressed in the skin. We determine that ∼78% of UTRs are formed by using canonical A[A/U]UAAA polyadenylation signals, whereas ∼22% of UTRs use alternative signals. By comparing to relative and absolute mRNA abundance determined by qPCR, our RNA-seq approach can precisely measure mRNA fold-change and accurately determine the expression of mRNAs over four orders of magnitude. Surprisingly, only 829 out of 12,739 genes show differential 3'-end usage between embryonic skin stem cells and their immediate daughter cells, whereas the numbers increase to 933 genes when comparing embryonic skin stem cells with the more remotely related hair follicle cells. This suggests an evolving diversity instead of switch-like dynamics in 3'-end formation during development. Finally, core components of the miRNA pathway including Dicer, Dgcr8, Xpo5, and Argonautes show dynamic 3'-UTR formation patterns, indicating a self-regulatory mechanism. Together, our quantitative analysis reveals a dynamic picture of mRNA 3'-end formation in tissue stem cell lineages in vivo.

Differential microRNA regulation correlates with alternative polyadenylation pattern between breast cancer and normal cells

Alternative polyadenylation (APA) could result in mRNA isoforms with variable lengths of 3′ UTRs. Gain of microRNA target sites in the 3′ UTR of a long mRNA isoform may cause different regulation from the corresponding short isoform. It has been known that cancer cells globally exhibit a lower ratio of long and short isoforms (LSR); that is, they tend to express larger amounts of short isoforms. The objective of this study is to illustrate the relationship between microRNA differential regulation and LSR. We retrieved public APA annotations and isoform expression profiles of breast cancer and normal cells from a high-throughput sequencing method study specific for the mRNA 3′ end. Combining microRNA expression profiles, we performed statistical analysis to reveal and estimate microRNA regulation on APA patterns in a global scale. First, we found that the amount of microRNA target sites in the alternative UTR (aUTR), the region only present in long isoforms, could affect the LSR of the target genes. Second, we observed that the genes whose aUTRs were targeted by up-regulated microRNAs in cancer cells had an overall lower LSR. Furthermore, the target sites of up-regulated microRNAs tended to appear in aUTRs. Finally, we demonstrated that the amount of target sites for up-regulated microRNAs in aUTRs correlated with the LSR change between cancer and normal cells. The results indicate that up-regulation of microRNAs might cause lower LSRs of target genes in cancer cells through degradation of their long isoforms. Our findings provide evidence of how microRNAs might play a crucial role in APA pattern shifts from normal to cancerous or proliferative states.

Alternative polyadenylation in glioblastoma multiforme and changes in predicted RNA binding protein profiles

Alternative polyadenylation (APA) is widely present in the human genome and plays a key role in carcinogenesis. We conducted a comprehensive analysis of the APA products in glioblastoma multiforme (GBM, one of the most lethal brain tumors) and normal brain tissues and further developed a computational pipeline, RNAelements (http://sysbio.zju.edu.cn/RNAelements/), using covariance model from known RNA binding protein (RBP) targets acquired by RNA Immunoprecipitation (RIP) analysis. We identified 4530 APA isoforms for 2733 genes in GBM, and found that 182 APA isoforms from 148 genes showed significant differential expression between normal and GBM brain tissues. We then focused on three genes with long and short APA isoforms that show inconsistent expression changes between normal and GBM brain tissues. These were myocyte enhancer factor 2D, heat shock factor binding protein 1, and polyhomeotic homolog 1 (Drosophila). Using the RNAelements program, we found that RBP binding sites were enriched in the alternative regions between the first and the last polyadenylation sites, which would result in the short APA forms escaping regulation from those RNA binding proteins. To the best of our knowledge, this report is the first comprehensive APA isoform dataset for GBM and normal brain tissues. Additionally, we demonstrated a putative novel APA-mediated mechanism for controlling RNA stability and translation for APA isoforms. These observations collectively lay a foundation for novel diagnostics and molecular mechanisms that can inform future therapeutic interventions for GBM.

All's well that ends well: alternative polyadenylation and its implications for stem cell biology

Stem cell quiescence, activation, and differentiation are governed by a complex network of molecular pathways. There has been a growing recognition that posttranscriptional modifications, such as alternative polyadenylation (APA) of transcripts, play an important role in regulating gene expression and function. Recent analyses of stem cell populations have suggested that APA controls stem cell fate and behavior. Here, we review recent developments that have shaped our understanding of the control of stem cell behavior by APA and we highlight promising areas for future investigation.

Alterations in polyadenylation and its implications for endocrine disease

INTRODUCTION

Polyadenylation is the process in which the pre-mRNA is cleaved at the poly(A) site and a poly(A) tail is added - a process necessary for normal mRNA formation. Genes with multiple poly(A) sites can undergo alternative polyadenylation (APA), producing distinct mRNA isoforms with different 3' untranslated regions (3' UTRs) and in some cases different coding regions. Two thirds of all human genes undergo APA. The efficiency of the polyadenylation process regulates gene expression and APA plays an important part in post-transcriptional regulation, as the 3' UTR contains various cis-elements associated with post-transcriptional regulation, such as target sites for micro-RNAs and RNA-binding proteins. Implications of alterations in polyadenylation for endocrine disease: Alterations in polyadenylation have been found to be causative of neonatal diabetes and IPEX (immune dysfunction, polyendocrinopathy, enteropathy, X-linked) and to be associated with type I and II diabetes, pre-eclampsia, fragile X-associated premature ovarian insufficiency, ectopic Cushing syndrome, and many cancer diseases, including several types of endocrine tumor diseases.

PERSPECTIVES

Recent developments in high-throughput sequencing have made it possible to characterize polyadenylation genome-wide. Antisense elements inhibiting or enhancing specific poly(A) site usage can induce desired alterations in polyadenylation, and thus hold the promise of new therapeutic approaches.

SUMMARY

This review gives a detailed description of alterations in polyadenylation in endocrine disease, an overview of the current literature on polyadenylation and summarizes the clinical implications of the current state of research in this field.

Analysis of alternative cleavage and polyadenylation by 3' region extraction and deep sequencing

Alternative cleavage and polyadenylation (APA) leads to mRNA isoforms with different coding sequences (CDS) and/or 3′ untranslated regions (3′UTRs). Using 3′ Region Extraction And Deep Sequencing (3′READS), a method which addresses the internal priming and oligo(A) tail issues that commonly plague polyA site (pA) identification, we comprehensively mapped pAs in the mouse genome, thoroughly annotating 3′ ends of genes and revealing over five thousand pAs (~8% of total) flanked by A-rich sequences, which have hitherto been overlooked. About 79% of mRNA genes and 66% of long non-coding RNA (lncRNA) genes have APA; but these two gene types have distinct usage patterns for pAs in introns and upstream exons. Promoter-distal pAs become relatively more abundant during embryonic development and cell differentiation, a trend affecting pAs in both 3′-most exons and upstream regions. Upregulated isoforms generally have stronger pAs, suggesting global modulation of the 3′ end processing activity in development and differentiation.

Alternative polyadenylation: new insights from global analyses

Recent studies have revealed widespread mRNA alternative polyadenylation (APA) in eukaryotes and its dynamic spatial and temporal regulation. APA not only generates proteomic and functional diversity, but also plays important roles in regulating gene expression. Global deregulation of APA has been demonstrated in a variety of human diseases. Recent exciting advances in the field have been made possible in a large part by high throughput analyses using newly developed experimental tools. Here I review the recent progress in global studies of APA and the insights that have emerged from these and other studies that use more conventional methods.

Cleavage factor Im is a key regulator of 3' UTR length

In eukaryotes, the 3' ends of RNA polymerase II-transcribed RNAs are generated in the majority of cases by site-specific endonucleolytic cleavage, followed by the addition of a poly(A) tail. Through alternative polyadenylation, a gene can give rise to multiple mRNA isoforms that differ in the length of their 3' UTRs and hence in their susceptibility to post-transcriptional regulatory factors such as microRNAs. A series of recently conducted high-throughput studies of poly(A) site usage revealed an extensive tissue-specific control and drastic changes in the length of mRNA 3' UTRs upon induction of proliferation in resting cells. To understand the dynamics of poly(A) site choice, we recently identified binding sites of the major pre-mRNA 3' end processing factors - cleavage and polyadenylation specificity factor (CPSF), cleavage stimulation factor (CstF), and cleavage factor Im (CF Im) - and mapped polyadenylation sites in HEK293 cells. Our present study extends previous findings on the role of CF Im in alternative polyadenylation and reveals that subunits of the CF Im complex generally control 3' UTR length. More specifically, we demonstrate that the loss-of-function of CF Im 68 and CF Im 25 but not of CF Im 59 leads to a transcriptome-wide increase in the use of proximal polyadenylation sites in HEK293 cells.

A global analysis of tandem 3'UTRs in eosinophilic chronic rhinosinusitis with nasal polyps

Background

Alternative polyadenylation (APA) is emerging as a widespread mechanism of gene regulation. The usage of APA sites allows a single gene to encode multiple mRNA transcripts with different 3′-untranslated region (3′UTR) lengths. Many disease processes reflect the importance of the regulation of APA site switching. The objective of this study was to explore the profiling of tandem APA sites in nasal polyps compared with nasal uncinate process mucosa.

Methods

Sequencing of APA sites (SAPAS) based on second-generation sequencing technology was undertaken to investigate the use of tandem APA sites and identify gene expression patterns in samples from the nasal polyps and nasal uncinate process mucosa of two patients with chronic rhinosinusitis with nasal polyps. The findings of the SAPAS analysis were validated via quantitative reverse-transcription polymerase chain reaction (qRT-PCR).

Results

First, the results showed a switching of 3′UTR lengths in nasal polyps compared with nasal uncinate process mucosa. From the two patients, 105 genes that were detected in both patients in the nasal polyps were switched to distal poly(A) sites, and 90 such genes were switched to proximal poly(A) sites. Several Gene Ontology terms were enriched in the list of genes with switched APA sites, including transcription regulation, cell cycle, apoptosis, and metabolism. Second, we detected genes that showed differential expression with at least a 3-fold difference between nasal polyp tissue and nasal uncinate process mucosa. Between the two sample types, 627 genes exhibited differential expression. The qRT-PCR results confirmed our SAPAS results.

Conclusion

APA site-switching events of 3′UTRs are prevalent in nasal polyp tissue, and the regulation of gene expression mediated by APA may play an important role in the formation and persistence of nasal polyps. Our results may provide new insights into the possible pathophysiologic processes involved in nasal polyps.

Regulation of eukaryotic gene expression by the untranslated gene regions and other non-coding elements

There is now compelling evidence that the complexity of higher organisms correlates with the relative amount of non-coding RNA rather than the number of protein-coding genes. Previously dismissed as "junk DNA", it is the non-coding regions of the genome that are responsible for regulation, facilitating complex temporal and spatial gene expression through the combinatorial effect of numerous mechanisms and interactions working together to fine-tune gene expression. The major regions involved in regulation of a particular gene are the 5' and 3' untranslated regions and introns. In addition, pervasive transcription of complex genomes produces a variety of non-coding transcripts that interact with these regions and contribute to regulation. This review discusses recent insights into the regulatory roles of the untranslated gene regions and non-coding RNAs in the control of complex gene expression, as well as the implications of this in terms of organism complexity and evolution.