Mechanisms and consequences of alternative polyadenylation

Analysis of C. elegans intestinal gene expression and polyadenylation by fluorescence-activated nuclei sorting and 3'-end-seq

Despite the many advantages of Caenorhabditis elegans, biochemical approaches to study tissue-specific gene expression in post-embryonic stages are challenging. Here, we report a novel experimental approach for efficient determination of tissue-specific transcriptomes involving the rapid release and purification of nuclei from major tissues of post-embryonic animals by fluorescence-activated nuclei sorting (FANS), followed by deep sequencing of linearly amplified 3'-end regions of transcripts (3'-end-seq). We employed these approaches to compile the transcriptome of the developed C. elegans intestine and used this to analyse tissue-specific cleavage and polyadenylation. In agreement with intestinal-specific gene expression, highly expressed genes have enriched GATA-elements in their promoter regions and their functional properties are associated with processes that are characteristic for the intestine. We systematically mapped pre-mRNA cleavage and polyadenylation sites, or polyA sites, including more than 3000 sites that have previously not been identified. The detailed analysis of the 3'-ends of the nuclear mRNA revealed widespread alternative polyA site use (APA) in intestinally expressed genes. Importantly, we found that intestinal polyA sites that undergo APA tend to have U-rich and/or A-rich upstream auxiliary elements that may contribute to the regulation of 3'-end formation in the intestine.

Extensive alternative polyadenylation during zebrafish development

The post-transcriptional fate of messenger RNAs (mRNAs) is largely dictated by their 3′ untranslated regions (3′ UTRs), which are defined by cleavage and polyadenylation (CPA) of pre-mRNAs. We used poly(A)-position profiling by sequencing (3P-seq) to map poly(A) sites at eight developmental stages and tissues in the zebrafish. Analysis of over 60 million 3P-seq reads substantially increased and improved existing 3′ UTR annotations, resulting in confidently identified 3′ UTRs for >79% of the annotated protein-coding genes in zebrafish. mRNAs from most zebrafish genes undergo alternative CPA, with those from more than a thousand genes using different dominant 3′ UTRs at different stages. These included one of the poly(A) polymerase genes, for which alternative CPA reinforces its repression in the ovary. 3′ UTRs tend to be shortest in the ovaries and longest in the brain. Isoforms with some of the shortest 3′ UTRs are highly expressed in the ovary, yet absent in the maternally contributed RNAs of the embryo, perhaps because their 3′ UTRs are too short to accommodate a uridine-rich motif required for stability of the maternal mRNA. At 2 h post-fertilization, thousands of unique poly(A) sites appear at locations lacking a typical polyadenylation signal, which suggests a wave of widespread cytoplasmic polyadenylation of mRNA degradation intermediates. Our insights into the identities, formation, and evolution of zebrafish 3′ UTRs provide a resource for studying gene regulation during vertebrate development.

Genome-wide analysis of pre-mRNA 3' end processing reveals a decisive role of human cleavage factor I in the regulation of 3' UTR length

Through alternative polyadenylation, human mRNAs acquire longer or shorter 3' untranslated regions, the latter typically associated with higher transcript stability and increased protein production. To understand the dynamics of polyadenylation site usage, we performed transcriptome-wide mapping of both binding sites of 3' end processing factors CPSF-160, CPSF-100, CPSF-73, CPSF-30, Fip1, CstF-64, CstF-64τ, CF I(m)25, CF I(m)59, and CF I(m)68 and 3' end processing sites in HEK293 cells. We found that although binding sites of these factors generally cluster around the poly(A) sites most frequently used in cleavage, CstF-64/CstF-64τ and CFI(m) proteins have much higher positional specificity compared to CPSF components. Knockdown of CF I(m)68 induced a systematic use of proximal polyadenylation sites, indicating that changes in relative abundance of a single 3' end processing factor can modulate the length of 3' untranslated regions across the transcriptome and suggesting a mechanism behind the previously observed increase in tumor cell invasiveness upon CF I(m)68 knockdown.

Translational control by changes in poly(A) tail length: recycling mRNAs

Beyond the well-known function of poly(A) tail length in mRNA stability, recent years have witnessed an explosion of information about how changes in tail length and the selection of alternative polyadenylation sites contribute to the translational regulation of a large portion of the genome. The mechanisms and factors mediating nuclear and cytoplasmic changes in poly(A) tail length have been studied in great detail, the targets of these mechanisms have been identified--in some cases by genome-wide screenings--and changes in poly(A) tail length are now implicated in a number of physiological and pathological processes. However, in very few cases have all three levels--mechanisms, targets and functions--been studied together.

MicroRNA-mediated posttranscriptional mechanisms of gene expression in proliferating and quiescent cancer cells

MicroRNAs are small non-coding RNA regulators of gene expression that play important roles in critical biological processes, including cell division, self-renewal and cell state maintenance. Their deregulation leads to extensive clinical consequences in tumorigenesis. Cancers demonstrate heterogeneity in their cell states implicated in their resistance and resurgence. Apart from proliferating cells, cancers harbor a small proportion of assorted quiescent cells that resist conventional therapeutics and contribute to cancer recurrence. MicroRNA expression, targets, microRNPs (microRNA-protein complexes) and their functions have been demonstrated to be regulated in distinct tumor cell states and as an adaptive response to stress signals in tumor-unfavorable environments. In turn, altered microRNPs and their modified post-transcriptional mechanisms of gene expression may contribute to tumor resistance and influence tumor progression. An understanding of distinct microRNA mechanisms in cancer cells would provide extensive insights into the versatile roles of microRNAs in the perpetuation of tumors and indicate potential therapeutic avenues.

PABPN1 shuts down alternative poly(A) sites

Although overlooked for many years, alternative cleavage and polyadenylation (APA) is now emerging as a major mechanism of gene regulation. A recent study identifies poly(A)-binding protein nuclear 1 (PABPN1), a general factor of polyadenylation, as a suppressor of alternative poly(A) sites.

Complex control of GABA(A) receptor subunit mRNA expression: variation, covariation, and genetic regulation

GABA type-A receptors are essential for fast inhibitory neurotransmission and are critical in brain function. Surprisingly, expression of receptor subunits is highly variable among individuals, but the cause and impact of this fluctuation remains unknown. We have studied sources of variation for all 19 receptor subunits using massive expression data sets collected across multiple brain regions and platforms in mice and humans. Expression of Gabra1, Gabra2, Gabrb2, Gabrb3, and Gabrg2 is highly variable and heritable among the large cohort of BXD strains derived from crosses of fully sequenced parents—C57BL/6J and DBA/2J. Genetic control of these subunits is complex and highly dependent on tissue and mRNA region. Remarkably, this high variation is generally not linked to phenotypic differences. The single exception is Gabrb3, a locus that is linked to anxiety. We identified upstream genetic loci that influence subunit expression, including three unlinked regions of chromosome 5 that modulate the expression of nine subunits in hippocampus, and that are also associated with multiple phenotypes. Candidate genes within these loci include, Naaa, Nos1, and Zkscan1. We confirmed a high level of coexpression for subunits comprising the major channel—Gabra1, Gabrb2, and Gabrg2—and identified conserved members of this expression network in mice and humans. Gucy1a3, Gucy1b3, and Lis1 are novel and conserved associates of multiple subunits that are involved in inhibitory signaling. Finally, proximal and distal regions of the 3′ UTRs of single subunits have remarkably independent expression patterns in both species. However, corresponding regions of different subunits often show congruent genetic control and coexpression (proximal-to-proximal or distal-to-distal), even in the absence of sequence homology. Our findings identify novel sources of variation that modulate subunit expression and highlight the extraordinary capacity of biological networks to buffer 4–100 fold differences in mRNA levels.

A molecular doorstop ensures a trickle through translational repression

Switching mRNA translation off and on is central to regulated gene expression, but what mechanisms moderate the extent of switch-off? Yao et al. describe how basal expression from interferon-gamma-induced transcripts is maintained during mRNA-specific translational repression. This antagonistic mechanism utilizes a truncated RNA-binding factor generated by a unique alternative polyadenylation event.

A quantitative atlas of polyadenylation in five mammals

We developed PolyA-seq, a strand-specific and quantitative method for high-throughput sequencing of 3′ ends of polyadenylated transcripts, and used it to globally map polyadenylation (polyA) sites in 24 matched tissues in human, rhesus, dog, mouse, and rat. We show that PolyA-seq is as accurate as existing RNA sequencing (RNA-seq) approaches for digital gene expression (DGE), enabling simultaneous mapping of polyA sites and quantitative measurement of their usage. In human, we confirmed 158,533 known sites and discovered 280,857 novel sites (FDR < 2.5%). On average 10% of novel human sites were also detected in matched tissues in other species. Most novel sites represent uncharacterized alternative polyA events and extensions of known transcripts in human and mouse, but primarily delineate novel transcripts in the other three species. A total of 69.1% of known human genes that we detected have multiple polyA sites in their 3′UTRs, with 49.3% having three or more. We also detected polyadenylation of noncoding and antisense transcripts, including constitutive and tissue-specific primary microRNAs. The canonical polyA signal was strongly enriched and positionally conserved in all species. In general, usage of polyA sites is more similar within the same tissues across different species than within a species. These quantitative maps of polyA usage in evolutionarily and functionally related samples constitute a resource for understanding the regulatory mechanisms underlying alternative polyadenylation.

Serotonin transporter polyadenylation polymorphism modulates the retention of fear extinction memory

Growing evidence suggests serotonin's role in anxiety and depression is mediated by its effects on learned fear associations. Pharmacological and genetic manipulations of serotonin signaling in mice alter the retention of fear extinction learning, which is inversely associated with anxious temperament in mice and humans. Here, we test whether genetic variation in serotonin signaling in the form of a common human serotonin transporter polyadenylation polymorphism (STPP/rs3813034) is associated with spontaneous fear recovery after extinction. We show that the risk allele of this polymorphism is associated with impaired retention of fear extinction memory and heightened anxiety and depressive symptoms. These STPP associations in humans mirror the phenotypic effects of serotonin transporter knockout in mice, highlighting the STPP as a potential genetic locus underlying interindividual differences in serotonin transporter function in humans. Furthermore, we show that the serotonin transporter polyadenylation profile associated with the STPP risk allele is altered through the chronic administration of fluoxetine, a treatment that also facilitates retention of extinction learning. The propensity to form persistent fear associations due to poor extinction recall may be an intermediate phenotype mediating the effects of genetic variation in serotonergic function on anxiety and depression. The consistency and specificity of these data across species provide robust support for this hypothesis and suggest that the little-studied STPP may be an important risk factor for mood and anxiety disorders in humans.

Coding region polyadenylation generates a truncated tRNA synthetase that counters translation repression

Posttranscriptional regulatory mechanisms superimpose "fine-tuning" control upon "on-off" switches characteristic of gene transcription. We have exploited computational modeling with experimental validation to resolve an anomalous relationship between mRNA expression and protein synthesis. The GAIT (gamma-interferon-activated inhibitor of translation) complex repressed VEGF-A synthesis to a low, constant rate independent of VEGF-A mRNA expression levels. Dynamic model simulations predicted an inhibitory GAIT-element-interacting factor to account for this relationship and led to the identification of a truncated form of glutamyl-prolyl tRNA synthetase (EPRS), a GAIT constituent that mediates binding to target transcripts. The truncated protein, EPRS(N1), shields GAIT-element-bearing transcripts from the inhibitory GAIT complex, thereby dictating a "translational trickle" of GAIT target proteins. EPRS(N1) mRNA is generated by polyadenylation-directed conversion of a Tyr codon in the EPRS-coding sequence to a stop codon (PAY(∗)). Genome-wide analysis revealed multiple candidate PAY(∗) targets, including the authenticated target RRM1, suggesting a general mechanism for production of C terminus-truncated regulatory proteins.

Genome-wide determination of a broad ESRP-regulated posttranscriptional network by high-throughput sequencing

Tissue-specific alternative splicing is achieved through the coordinated assembly of RNA binding proteins at specific sites to enhance or silence splicing at nearby splice sites. We used high-throughput sequencing (RNA-Seq) to investigate the complete spectrum of alternative splicing events that are regulated by the epithelium-specific splicing regulatory proteins ESRP1 and ESRP2. We also combined this analysis with direct RNA sequencing (DRS) to reveal ESRP-mediated regulation of alternative polyadenylation. To define binding motifs that mediate direct regulation of splicing and polyadenylation by ESRP, SELEX-Seq analysis was performed, coupling traditional SELEX with high-throughput sequencing. Identification and scoring of high-affinity ESRP1 binding motifs within ESRP target genes allowed the generation of RNA maps that define the position-dependent activity of the ESRPs in regulating cassette exons and alternative 3' ends. These extensive analyses provide a comprehensive picture of the functions of the ESRPs in an epithelial posttranscriptional gene expression program.

Genome-wide polyadenylation site mapping

Alternative polyadenylation site usage gives rise to variation in 3' ends of transcripts in diverse organisms ranging from yeast to human. Accurate mapping of polyadenylation sites of transcripts is of major biological importance, since the length of the 3'UTR can have a strong influence on transcript stability, localization, and translation. However, reads generated using total mRNA sequencing mostly lack the very 3' end of transcripts. Here, we present a method that allows simultaneous analysis of alternative 3' ends and transcriptome dynamics at high throughput. By using transcripts produced in vitro, the high precision of end mapping during the protocol can be controlled. This method is illustrated here for budding yeast. However, this method can be applied to any natural or artificially polyadenylated RNA.

Functional analysis of arylamine N-acetyltransferase 1 (NAT1) NAT1*10 haplotypes in a complete NATb mRNA construct

N-acetyltransferase 1 (NAT1) catalyzes N-acetylation of arylamines as well as the O-acetylation of N-hydroxylated arylamines. O-acetylation leads to the formation of electrophilic intermediates that result in DNA adducts and mutations. NAT1*10 is the most common variant haplotype and is associated with increased risk for numerous cancers. NAT1 is transcribed from a major promoter, NATb, and an alternative promoter, NATa, resulting in messenger RNAs (mRNAs) with distinct 5'-untranslated regions (UTRs). To best mimic in vivo metabolism and the effect of NAT1*10 polymorphisms on polyadenylation usage, pcDNA5/Flp recombination target plasmid constructs were prepared for transfection of full-length human mRNAs including the 5'-UTR derived from NATb, the open reading frame and 888 nucleotides of the 3'-UTR. Following stable transfection of NAT1*4, NAT1*10 and an additional NAT1*10 variant (termed NAT1*10B) into nucleotide excision repair-deficient Chinese hamster ovary cells, N- and O-acetyltransferase activity (in vitro and in situ), mRNA and protein expression were higher in cells transfected with NAT1*10 and NAT1*10B than in cells transfected with NAT1*4 (P < 0.05). Consistent with NAT1 expression and activity, cytotoxicity and hypoxanthine phosphoribosyl transferase mutants following 4-aminobiphenyl exposures were higher in NAT1*10 than in NAT1*4 transfected cells. Ribonuclease protection assays showed no difference between NAT1*4 and NAT1*10. However, protection of one probe by NAT1*10B was not observed with NAT1*4 or NAT1*10, suggesting additional mechanisms that regulate NAT1*10B. The higher mutants in cells transfected with NAT1*10 and NAT1*10B are consistent with an increased cancer risk for individuals possessing NAT1*10 haplotypes.

Transcriptional priming of cytoplasmic post-transcriptional regulation

Recent evidence suggests that post-transcriptional events that are executed in the cytoplasm can be predetermined during transcription. Here, I speculate that this is a widespread mode of regulation and discuss the potential mechanisms, advantages and implications of such a regulatory strategy.

RNA regulatory elements and polyadenylation in plants

Alternative poly(A) site choice (also known as alternative polyadenylation, or APA) has the potential to affect gene expression in qualitative and quantitative ways. APA may affect as many as 82% of all expressed genes in a plant. The consequences of APA include the generation of transcripts with differing 3'-UTRs (and thus differing regulatory potential) and of transcripts with differing protein-coding potential. Genome-wide studies of possible APA suggest a linkage with pre-mRNA splicing, and indicate a coincidence of and perhaps cooperation between RNA regulatory elements that affect splicing efficiency and the recognition of novel intronic poly(A) sites. These studies also raise the possibility of the existence of a novel class of polyadenylation-related cis elements that are distinct from the well-characterized plant polyadenylation signal. Many potential APA events, however, have not been associated with identifiable cis elements. The present state of the field reveals a broad scope of APA, and also numerous opportunities for research into mechanisms that govern both choice and regulation of poly(A) sites in plants.