No retraction, correction or concern found when last checked on Oct 24, 2026.
Cleavage and polyadenylation: Ending the message expands gene regulation
Jonathan Neve, Radhika Patel, Zhiqiao Wang and 2 more
RNA Biology · 2017 · Review · CC BY
doi:10.1080/15476286.2017.1306171
Abstract
Cleavage and polyadenylation (pA) is a fundamental step that is required for the maturation of primary protein encoding transcripts into functional mRNAs that can be exported from the nucleus and translated in the cytoplasm. 3'end processing is dependent on the assembly of a multiprotein processing complex on the pA signals that reside in the pre-mRNAs. Most eukaryotic genes have multiple pA signals, resulting in alternative cleavage and polyadenylation (APA), a widespread phenomenon that is important to establish cell state and cell type specific transcriptomes. Here, we review how pA sites are recognized and comprehensively summarize how APA is regulated and creates mRNA isoform profiles that are characteristic for cell types, tissues, cellular states and disease.
Full text
From PubMed Central, U.S. National Library of Medicine, as current on Oct 24, 2026. Paragraph numbers are ours, so an answer can point at one.
Introduction
To express genes, the information present in the DNA blueprint is transcribed by RNA polymerases into RNA molecules. These initial or primary transcripts are often subjected to modifications and processing reactions that mature them into functional molecules. While the processing of some noncoding RNAs such as rRNA and tRNAs are also essential in prokaryotes, extensive processing of protein encoding pre-mRNA transcripts is a unique feature of eukaryotic gene expression.
In eukaryotes, 3 major processing reactions are responsible for the maturation of pre-mRNA molecules into mRNAs that are competent for nuclear cytoplasmic export and subsequent translation. Firstly, the 5′ends of RNAs are modified by the addition of a cap structure in the form of a guanosine molecule that is enzymatically attached to the first nucleotide of the transcribed mRNA by an unusual 5′–5′ triphosphate linkage and subsequently modified by methylation. The genes that encode proteins are mostly not continuous, but separated into alternating exons (coding regions) and introns, which are noncoding regions that need to be removed. This occurs during splicing, where a large multi-protein-RNA complex, known as the spliceosome, directs precise excision of introns from the primary transcript and fuses the exons together. This generates an mRNA with a continuous open reading frame (ORF) that serves as a template for protein synthesis in the cytoplasm. In order for the transcript to be exported to the cytoplasm, the transcript needs to be detached from RNA Polymerase, which occurs by a 2-step cleavage and polyadenylation reaction. This forms the typical uniform polyadenylated 3′end of almost all of the protein encoding nuclear mRNAs.
While splicing has long been recognized to play an important role in the regulation of gene expression in eukaryotes, the scope for gene regulation at the level of cleavage and polyadenylation has only recently been fully recognized.
Core cleavage and polyadenylation sequence elements
With the exception of replication dependent histone genes, metazoan protein encoding mRNAs contain a uniform 3′end consisting of a stretch of adenosines. This poly(A) tail is synthesized in a template independent fashion and is the result of an RNA processing reaction, where a multiprotein complex assembles on specific sequences on the pre-mRNAs, called the cleavage and polyadenylation signals (pA signals). pA signals consist of sequences that flank either side of where the pre-mRNA is endonucleolytically cleaved and subsequently polyadenylated (Fig. 1). The core pA signal is a bipartite sequence element that constitutes an A-rich hexameric sequence, found 21 nucleotides upstream of the cleavage site and a U and/or G/U rich sequence, located 10–30 nucleotides downstream of the cleavage site. Whereas the upstream hexamer motif for most (70–75%) pA signals conforms to the canonical A[A/U]UAAA hexamer, the nucleotide composition of the downstream sequence element (DSE) can vary significantly between pA signals. Although no clear consensus can be defined, often a GU rich sequence followed by a stretch of 3 or more uridines is present within a 40 nucleotide wide window downstream of the cleavage site.
Core sequences upstream of the cleavage site
The importance of the nucleotide composition in the hexamer sequence element is underpinned by numerous in vivo and in vitro experimental examples. Single nucleotide changes in the A[A/U]UAAA sequence, or its entire removal, have proven to severely impair cleavage efficiency. These laboratory based observations are further strengthened by the naturally occurring polymorphisms such as found in the human globin genes, which are associated with the thalassemia phenotype. However, these single nucleotide polymorphisms (SNPs) rarely occur in the pA hexamer, further emphasizing the importance of the hexamer nucleotide composition. However, naturally occurring non-canonical variants, with a deviation of one or more nucleotides from the A[A/U]UAAA hexamer sequence consensus, represent about 10-20% of all hexamer signals.
Core sequences downstream of the cleavage site
In contrast to the hexamer, the DSE has long been considered to be much more tolerant to single nucleotide changes. However, at least for a subset of genes, in particular those with degenerated hexamer sequences, the opposite is true: while mutations in the hexamer had little effect, base changes in the DSE reduced 3′end processing significantly. That the DSE can be a decisive element in regulating efficiency of cleavage is further supported by clinically relevant examples where changes in the DSE force a gain of function and increase the processing capacity of a particular pA site. In addition, knockdown of CstF64, the component of the 3′end processing machinery known to interact with the DSE (see below), affects usage of pA sites that are characterized by a non-canonical hexamer and a GUKKU “type” DSE. Furthermore, for a subset of genes that have pA sites with degenerated DSEs, CstF64 is not required for 3′end processing. The mRNAs of these genes are polyadenylated by Star-PAP, a non-canonical poly(A) polymerase. These observations suggest that for some pA sites, the DSE may represent the critical core element that drives cleavage efficiency and directs how the 3′end machinery is assembled at the pA site.
The cleavage site
The point of cleavage occurs between the hexamer and DSE. In mammals, CA and UA appear to be the most frequent dinucleotides that precede the actual site of cleavage, however, the actual site where cleavage occurs is known to be heterogeneous. Nucleotide composition at the cleavage site has a complex impact on 3′end processing. Initial in vitro based studies revealed that mutations of the dinucleotides at the cleavage site generally resulted in relatively moderate effects on cleavage efficiency. Conversely, a prominent SNP occurring at a 1% frequency in the Caucasian population, altering the cleavage site of the prothrombin gene from CG to CA, significantly increases cleavage efficiency of the pre-mRNA, and has been linked to a higher tendency toward thrombophilia.
Auxiliary elements
In addition to the core sequence elements, many pA sites possess additional auxiliary sequences that influence their overall strength. These are diverse in location and nucleotide composition, but in most cases feature U/G rich sequences. The presence of these auxiliary elements correlates with deviation from the consensus A[A/U]UAAA hexamer, perhaps compensating for the loss of core sequence integrity. These “processing enhancers” can be found both upstream and downstream of the cleavage sites.
G-rich auxiliary sequences are some of the best characterized downstream enhancers that may fold into G-quadruplex structures. Several viral and mammalian pA sites have been described where these G-rich sequences significantly influence cleavage efficiency at the pA site by interacting with factors, such as hnRNPH, that stimulate processing. Because the G-rich sequences can exert their influence over distances that span several hundred bases, the true extent to which pA sites depend on these distant enhancers may be undervalued.
U-rich auxiliary regions are frequently found upstream and downstream of the cleavage site. The upstream U-rich sequences can be divided into 2 categories. The first category represents U-rich elements that are located between the hexamer and the actual cleavage site. They can function as binding sites for the CPSF subunit Fip1, which can stimulate the polyadenylation step and control pA site selection. The second category includes U-rich sequences that are positioned in the UTR upstream of the hexamer. These upstream elements (USEs) were first identified in several viral transcripts including SV40, Cauliflower Mosaic Virus (CaMV) and HIV-1 In the HIV-1 transcript the USE, identified between 56 and 93 nucleotides upstream of the hexamer, was found to interact with CPSF160 and aid its interaction with the poly(A) signal, proving to be critical for the cleavage reaction. USEs generally have no consensus sequences, although tend to be U-rich. A number of USEs have now been identified in transcripts of the human genome. In the lamin B2 and C2 complement genes, auxiliary USEs are required for interaction with PTB and CstF64 respectively, which is critical for the cleavage and polyadenylation reaction. A USE controls the usage of the poly(A) sites in the prothrombin gene in a stress response, and efficiency of polyadenylation in collagen genes has also been shown to be modulated by USEs. In addition, the usage of the proximal poly(A) site of the cyclo-oxygenase 2 gene is controlled by 3 U-rich USEs that bind a protein complex to recruit the core polyadenylation factors.
A UGUA motif has been identified as a USE, which is bound by the cleavage factor I (CFI) complex. UGUA mediated CFIm recruitment to the pre-mRNA can direct cleavage and polyadenylation at non-canonical sites, by anchoring Fip1 and the poly(A) polymerase (PAP) to the pre-mRNA. More recent whole transcriptome based analyses identified a more global role for this interaction and recognized CFI as a key regulator of pA site choice during alternative cleavage and polyadenylation (APA) for a significant number of genes. Finally, some of the U-rich elements have also been shown to mediate an interaction with hnRNPC similarly influencing pA site selection during APA.
AUA auxiliary element. A subset of mRNAs are polyadenylated by the non-canonical poly(A) polymerase Star-PAP (see below), which binds an AUA motif in the 3′UTR of its target pre-mRNAs.
The core poly(A) machinery
In mammals there are more than 80 proteins associated with the 3′end processing machinery, but this can be condensed to around 20 proteins that constitute the core factors (Fig. 2). Part of this core machinery includes 4 multi-subunit protein complexes which are highly conserved: the cleavage and polyadenylation specificity factor (CPSF), the cleavage stimulation factor (CstF), and mammalian cleavage factors I and II (CFIm and CFIIm). Additional core factors include symplekin, the poly(A) polymerase (PAP), nuclear poly(A) binding protein (PABPN1) and the C-terminal domain (CTD) of the largest subunit of RNA polymerase II (Pol II). The process of 3′end formation is initiated by the recognition of the respective cis-elements on the nascent transcripts by CPSF and CstF complexes in a cooperative manner. Depending on the pA site, this initial step is supported and often dependent on the interaction of other core factors (for example CFI) and also auxiliary factors (for example hnRNP H). After a functional complex is assembled cleavage and subsequent polyadenylation of the pre-mRNA can commence.
CPSF
The CPSF protein complex features 6 protein subunits; WDR33, CPSF30, CPSF160, hFip1, CPSF100 and CPSF73, which are required for cleavage of the pre-mRNA. In contrast, only 4 CPSF subunits, CPSF160, CPSF30, WDR33 and hFip1 together with PAP are sufficient to direct AAUAAA dependent poly(A) tail addition in vitro. The different dependencies on the CPSF subunits raises the intriguing possibility that different sub-complexes may be specific for the cleavage or the polyadenylation steps. This idea is supported by the finding that different CPSF sub-complexes have been identified that cater either for pA or histone 3′end processing.
WDR33 (WDR33) is a large 146 kDa protein, and despite the essential role its homolog Pfs2p plays in 3′end formation in yeast, its role during mammalian cleavage and polyadenylation has only recently been uncovered. The mammalian WDR33 was first shown to be essential in cleavage and polyadenylation in vitro, and has subsequently been identified as the CPSF subunit that directly binds the AAUAAA hexamer. Its importance for pA site recognition is further underpinned by PAR-CLIP data, which shows a strong interaction between WDR33 and the pre-mRNA directly upstream (−16 to −18 nt) of cleavage sites, close to the region where hexamers peak in transcripts. It is currently unclear how its specific interaction with AAUAAA is achieved, but the WD40 domain, which is principally associated with protein-protein interactions, also has RNA-binding properties.
CPSF30 (CPSF4) is the smallest subunit of the CPSF complex. The involvement of a 30 kDa protein in the recognition of AAUAAA has long been recognized, but it took more than 25 years to show that this protein, together with WDR33, is responsible for the specific interaction of CPSF with the AAUAAA hexamer in mammals. Interestingly, CPSF30 has also been shown to play a central role in pA site selection in plants. The RNA-binding properties of CPSF30 reside in the 5 zinc finger (ZF) domains, primarily mediated by ZF2 and ZF3. The binding of CPSF30 to RNA is required for the entire CPSF complex to interact with the RNA substrate. Notably, the CPSF30-AAUAAA interaction is specifically inhibited by the influenza protein NS1A, that blocks ZF2 and ZF3, and so globally abrogates host 3′end processing. These observations clearly highlight the central role of CPSF30 in pA site recognition.
CPSF160 (CPSF1) was considered to be the major component that mediates the binding of the CPSF complex to AAUAAA. This was primarily based on the observation that a 160 kDa protein can be specifically cross linked to AAUAAA-containing RNA substrates and the Yhh1/Cft1p, the yeast homolog of CPSF160, also bound RNA. Furthermore, pull-down assays confirmed that recombinant CPSF160 interacts with AAUAAA-containing RNAs. However, as mentioned above, it is now clear that AAUAAA-specific binding of the CPSF complex falls to WDR33 and CPSF30 and thus the precise role of CPSF160 for the interaction with the pre-mRNA needs to be re-evaluated. Nevertheless, CPSF160 is a critical component of the cleavage and polyadenylation machinery and establishes important protein-protein interactions. These include interactions between CPSF and CstF by directly contacting the CstF77 subunit. This interaction may be integral to the cooperative nature of pA site recognition. In addition, CPSF160 also interacts functionally with the PAP.
hFip1 (FIP1L1), the human factor interacting with PAP, was discovered based on the sequence similarity to the yeast polyadenylation factor Fip1p. True to its name, Fip1, and hFip1, interact with and recruit PAP to the pA site where it can simultaneously bind to U-rich regions that are often located either upstream or downstream of the pA hexamer. The RNA-binding activity of hFip1 lies within the arginine-rich C-terminus. A role in pA site recognition is further supported by the fact that hFip1 levels change in stem cells compared with differentiated cells. This directs APA during stem cell self-renewal. Higher levels of hFip1 in ESCs and iPSCs promote the recognition of weaker, hFip1 dependent pA sites.
CPSF100 (CPSF2) and CPSF73 (CPSF3), which both contain metallo-β-lactamase and β-CASP domains, form a heterodimer in mammals. The precise role of CPSF100 for cleavage and polyadenylation in mammals is still unclear. There is however evidence that CPSF100 is involved in THOC5 (a member of the TREX export complex) mediated control of 3′end processing of immediate early genes. THOC5 recruits CPSF100 to the 3′end of immediate early genes that are regulated by THOC5. The CASP domain in CPSF73 provides the zinc-dependent endonuclease activity during 3′end formation for both the cleavage and polyadenylation and the histone 3′end processing machineries. While they play a central role during cleavage, both proteins are superfluous for the polyadenylation step.
CstF
CstF is the component of the cleavage and polyadenylation machinery that specifically recognizes the DSE region in pA sites. It is a heterotrimeric protein complex containing the subunits CstF77, CstF64 and CstF50. The CstF complex is essential for the cleavage reaction but is not required for the polyadenylation step.
CstF64 (CSTF2) mediates the binding to the DSE in the pre-mRNA. The interaction interface resides in an RNA recognition motif (RRM) at its N-terminus, which specifically recognizes U/GU-rich regions downstream of the cleavage sites. How CstF64 is able to recognize diverse DSEs in pre-mRNAs remains uncertain, but consecutive Us are required for strong CstF-GU interactions. Interestingly, RNA binding also induces structural changes in CstF64 that may be critical to facilitate assembly of larger complexes. A global analysis that mapped transcriptome wide CstF64-RNA interactions in vivo suggests a role for CstF64 in modulating the selection of pA sites during APA. This promotes the usage of non-canonical pA sites that feature a GUKKU motif in their DSE.
CstF77 (CSTF3) bridges the CstF64 and CstF50 subunits and thus plays a key role for the assembly of the CstF complex. CstF77 appears to be an elongated dimer, suggesting that CstF may form a hexomeric complex in cells. CstF77 also appears to be critical for the establishment of protein interactions between CstF and other pA complexes. Cross linking experiments demonstrate a strong interaction between CPSF160 and CstF77, which may provide the basis for the observed CPSF-CstF cooperative RNA binding during the assembly process.
CstF50 (CSTF1), like CstF77, exists as a homodimer providing further evidence for the hexameric nature of functional CstF complexes. CstF50 contains 7 repeats of the WD40 motif, which mediates protein-protein interactions with CstF77 and BRCA1 associated RING domain 1 (BARD1). The latter interaction is inhibitory and believed to suppress erroneous 3′end formation possibly during transcription coupled DNA repair. CstF50 also interacts with the CTD of Pol II and this interaction may be established via the N-terminal part of CstF50.
CFIm is a complex is made up of 4 polypeptides: CFIm25 (CPSF5/NUDT21), CFIm68 (CPSF6), CFIm59 (CPSF7) and CFIm72 (alternative spliced isoform of CFIm68). It has been suggested that different forms of CFIm complex exist in vivo including CFIm25-68, CFIm25-59 and possibly CFIm25-72. As shown by SELEX experiments and site directed mutagenesis, CFIm recognizes a UGUAN motif in RNA through the RRM in CFIm25. The RRM in CFIm68 is critical for the interaction with CFIm25 and enhances RNA-binding of the complex. CFIm68 also facilitates RNA looping, which may play a role in the regulation of APA. CFIm has originally been shown to stabilize the binding of CPSF complex to the pre-mRNA and then emerged as a critical factor for the recognition of pA sites that lack the A[A/U]UAAA hexamer. To date, CFI is arguably the best characterized factor that regulates APA at the point of cleavage. This is demonstrated by numerous examples where a depletion of CFIm results in the preferred usage of the proximal pA sites. CFIm has also been implicated in regulating alternative usage of non-canonical pA sites during spermatogenesis. There are also contacts between CFIm68 and THOC5, a component of the transcription export complex (TREX). THOC5 appears to be important for the recruitment of CFIm68 to distal pA sites.
CFIIm, mammalian CFII, is only required for the cleavage step and is the least well characterized complex of the core component. To date, this complex has only been partly purified and can be separated into 2 fractions: CFIIAm and CFIIBm. The CFIIAm fraction contains the essential components for the cleavage reaction while CFIIBm is not required but has a stimulatory role for cleavage. One of the 15 polypeptides associated with CFIIAm, hClp1 has been shown to interact with CFIm and CPSF by immunoprecipitation experiments, suggesting that it bridges these 2 during the cleavage process. An additional component of CFII is protein 1 of cleavage factor 1, hPcf11, which has been shown to be critical in the control of Pol II transcription termination, by mediating the degradation of the polymerase associated 3′ RNA product after cleavage at the pA site. More recently, hPcf11 has also been implicated in regulating APA.
PAP (PAPOLA) is a monomer and responsible for the addition of the 3′ polyadenosine tail to a newly synthesized pre-mRNA molecule by catalyzing the chemical reaction ATP + RNA-3′OH → pyrophosphate + RNApA-3′OH. PAP adds 200–250 adenosines in a template-independent manner to the 3′end of the cleaved pre-mRNA. PAP is recruited to the 3′end processing complex by CPSF through interactions mediated by CPSF160 and hFip1. There are several PAP isoforms present in human cells which can chiefly be separated into canonical PAPs and non-canonical PAPs. Most RNAs are polyadenylated by canonical PAPs. Non-canonical PAPs, among other functions, are associated with mitochondrial mRNA polyadenylation (hmtPAP), cytoplasmic polyadenylation (hGld2) and miRNA biogenesis. A prominent non-canonical PAP called Star-PAP, controls the polyadenylation of a subset of mRNAs that encode for proteins that are associated with DNA damage induced apoptosis and stress responses.
PABPN1 (PABPN1) has been shown to contain a single RRM and an arginine rich CTD. The rate at which PAP adds adenosine nucleotides is dependent on the presence of PABPN1. The first few nucleotides added by PAP are added with low efficiency until the short polyadenosine tail is bound by PABPN1, which accelerates the rate of adenosine addition by PAP. While PAP adds 200–250 adenosine nucleotides to the 3′end of the mRNAs, the median length of the poly(A) tail in cellular mRNAs is much shorter, between 50 to 100 nucleotides long due to the action of cytoplasmic deadenylases. In addition to its primary function, PABPN1 has also been implicated in controlling APA.
Symplekin (SYMPK) is a protein that forms a high-molecular weight complex with CPSF and CstF by interacting with CPSF73, CPSF100 and CstF64. It is suggested to serve as a scaffold for recruiting other factors to the cleavage and polyadenylation complex. It also participates in the assembly of a processing complex that matures histone mRNA 3′ends, which do not undergo polyadenylation. Symplekin has also been found to form a complex with heat shock transcription factor 1 (HSF1) after stress treatment. This complex may be critical for the recruitment of the 3′end processing machinery to the heat shock protein HSB pre-mRNA during stress and so safeguard its expression. Symplekin appears to provide a scaffold around which different CPSF sub-complexes are assembled. Perhaps these different sub-complexes provide the necessary scope and flexibility to the 3′end machinery that is required to associate with pA sequences that differ widely in their architecture.
Pol II CTD (RBP1) is the key platform that couples pre-mRNA processing to transcription and it is required for efficient cleavage at the pA site. The human Pol II CTD features 21 consensus YSPTSPS and 31 non-consensus heptad repeats which form and present a dynamic interaction platform for a variety of pre-mRNA processing components including cleavage and polyadenylation factors such as Pcf11, CstF77 and CstF50.
The assembly of a functional cleavage and polyadenylation complex
The initiating step of the 3′end processing complex assembly (Fig. 2) is the coordinated recognition of the hexamer signal AAUAAA and the DSE on the nascent RNA by CPSF and CstF. The recruitment of CPSF is central, as it constitutes the core processing complex required for both the cleavage and the subsequent polyadenylation reaction. It binds the RNA directly by the association of WDR33 and CPSF30 with the hexamer. Whether the 2 subunits bind to pre-mRNAs concurrently or in a time and/or pA site type dependent manner remains unknown. It may well be that different CPSF complexes provide a flexible machinery that can recognize different types of pA signals. Unlike CPSF, the CstF complex is essential only for the cleavage reaction. CstF associates with di-uridine pockets in the U/GU rich DSE via its RRM in CstF64. CstF77 is a critical component that bridges the CPSF and CstF components, facilitating the cooperative binding. With CPSF and CstF anchored to the pA site, additional factors join and complete the assembly of a functional 3′end processing complex. The complex can then catalyze a 2-step reaction comprising of an endonucleolyic cleavage at the cleavage site and the polyadenylation of the 3′end of the 5′cleavage product. RNA cleavage occurs between the hexamer and the DSE, generally within a window of 20 nucleotides, either side of these core sequences.
The speed of assembly of a functional cleavage and polyadenylation complex is dependent on the strength of the pA site, but is generally achieved within about 10 seconds. The absolute strength of a given pA site is highly complex and depends on the sequence architecture of the pA site (see above) and the relative availability of core and auxiliary promoting or inhibitory factors.
Interconnection between pA cleavage, capping and splicing
The assembly process of the functional machinery is further complicated by the interconnection between 3′end processing with capping, pre-mRNA splicing and transcription. The relationship between capping and 3′end formation was first suggested after capped RNAs proved to be better substrates for 3′end processing in nuclear extracts. These initial observations were confirmed and shown to be mediated by the cap binding complex, CBC.
The relationship between splicing and 3′end formation is extensive and principally relies on direct interactions between splicing and cleavage and polyadenylation components. The positive influence of splicing factors on 3′end formation was first reported over 25 years ago and has since been extensively documented. Intron insertion into a replication dependent histone gene reporter is one example of this interconnection, where a splicing event in the histone pre-mRNA is concomitant with the activation of cryptic pA sites. As a result, cleavage and polyadenylation in the spliced transcripts is favored over U7 snRNP dependent 3′end formation, which is typical for replication dependent histone genes (reviewed in Romeo and Schümperli, 2016). Splicing has not only been shown to activate cryptic pA sites but generally stimulates cleavage and polyadenylation of the downstream pA sites. This stimulation is achieved by interactions between components of the splicing machinery that associate with the 3' splice sites (3SS) and several poly(A) factors notably between the 65 kDa subunit of U2AF, PAP and CFIm59. In addition, physical interactions of CPSF components with several subunits of the U2 snRNP, including SF3b155, SF3b130 and SF3b49 have been found. Correspondingly, mutations of U2 snRNA-binding sites significantly reduced the cleavage efficiency in reporter genes. The stimulatory effect through these interactions is mutual as mutations of the hexamer, can inhibit terminal intron removal and the depletion of CPSF100 impaired upstream splicing in an in vitro coupled splicing and 3′ processing system. CPSF and symplekin have been implicated in the promotion of alternative splicing on a global scale which further emphasizes the close connections of splicing and 3′end formation.
While splicing components associated with the 3SS appear to stimulate 3′end formation, the opposite is generally true for splicing components that bind to the 5SS, notably the U1 snRNP. The suppression of pA sites by U1 snRNP has proven to be instrumental for the inhibition of late gene expression in the bovine papilloma type 1 virus and in HIV-1, where the U1 snRNP bound to 5SS acts as a suppressor of upstream pA site usage. In the latter example, suppression of the upstream pA site is crucial to enable transcription of viral protein-encoding mRNA and the genomic RNA. This U1 snRNP-mediated type of suppression of 3′end formation is not restricted to viruses but was suggested to prevent premature transcription termination in polycistronic transcription units at a global scale in Caenorhabditis elegans by suppressing the recognition of the pA sites in upstream positioned genes in polycistronic mRNAs by Pol II associated poly(A) factors. Importantly, the suppressive role of U1 snRNP for 3′end formation emerged as a general mechanism by which cryptic pA sites in eukaryotic introns are suppressed globally, thereby protecting the cell from aberrant pA usage and premature transcription termination. Physiologically, U1 snRNP levels are very high which may ensure that levels never fall below a critical threshold that would result in the activation of these cryptic pA sites. However, U1 snRNP levels in cells can fluctuate, for example as a result of UV-induced DNA damage which causes increased usage of intronic alternative pA sites.
The polyadenylation step
The second part of the reaction comprises of the non-templated addition of about 200 adenosines, which in reconstituted systems only requires the cleaved pre-mRNA template, CPSF, PAP and the poly(A) binding protein (PABPN1, PABPII). PAP is initially tethered to the cleaved pre-mRNA via interactions with the CPSF subunits 160 kDa and Fip1. This loose interaction leads to regular dissociation of PAP and adenylation is restricted to intermittent addition of a few adenosines. Subsequently, the nuclear poly(A) binding protein, PABPN1, binds these short oligoadenosine tails added by PAP in its distributive mode. PABPN1-binding provides an additional anchor point for PAP and increases its affinity to the RNA 80-fold. Once this quaternary complex is established on the cleaved transcript, CPSF and PABPN1 cooperatively shift PAP into a processive mode, enabling it to catalyze the full length poly(A) tail without dissociation. The actual number of adenosines added to transcripts is uniform within species, ranging from 70 to 90 in Saccharomyces cerevisiae to around 200 in mammals. Extending the poly(A) tail to an appropriate length appears to be critical for appropriate gene expression, as both hypo- and hyperadenylation negatively influence the fate of the mRNA. Hypoadenylated mammalian transcripts have been shown to be retained in the nucleus, and inefficient polyadenylation in yeast is linked to nuclear mRNA degradation. In addition, hyperadenylated host transcripts observed in herpesvirus-infected cells are retained in the nucleus and funnelled into the nuclear mRNA degradation pathway. Furthermore, inhibition of RNA export forces hyperadenylation of nuclear mRNAs.
How the precise length of poly(A) tail synthesis is achieved by the polyadenylation complex has been the subject of intensive investigations and much has been learned from reconstituted in vitro systems. Central to the model derived from these observations is that the decoration of the poly(A) tail with PABPN1 provides some means of measuring and controlling the number of adenosines added to the 3′end of the cleaved mRNA. The current model envisages that the growing poly(A) tail coated with PABPN1 folds into a spherical complex that promotes contacts between CPSF, PAP and PABPN1. When the tail extends to 250 adenosines, new PABPN1 proteins can no longer be integrated into the sphere and this leads to the disruption of the interactions between CPSF and PAP. Subsequently PAP remains tethered to the RNA only by its interaction with PABPN1 causing it to readopt the distributive mode and ultimately terminate polyadenylation. A direct involvement of PABPN1 in poly(A) length control is further supported by in vivo evidence. siRNA mediated depletion of PABPN1 in mouse myoblasts triggered shortening of poly(A) tails and subsequent accumulation of these transcripts in the nucleus. Similarly, nuclear retained transcripts with shortened adenosine tails are characteristic for influenza A virus infected cells. The viral NS1 protein appears to restrict the length of poly(A) tails on host transcripts in infected cells to 12 nucleotides by sequestering PABPN1.
While CPSF, PAP and PABPN1 are sufficient to direct poly(A) tail length control in reconstituted systems, it appears that in vivo, more factors are required. The protein nucleophosmin (NPM1) has been found to associate with many mRNAs only after proper termination of polyadenylation. Intriguingly, RNAi knockdown of NPM1 leads to hyperadenylation and nuclear retention of mRNAs suggesting a role of NPM1 in poly(A) tail measurement. The precise molecular mechanisms are not clear but NPM1 somehow contributes to poly(A) tail sizing by regulating the dissociation of the quaternary complex. Perhaps NPM1 functions as a gatekeeper and associates with correctly polyadenylated and export competent transcripts.
Alternative cleavage and polyadenylation
The complexity of the pA site architecture combined with the multitude of core and auxiliary interacting factors and interconnection with splicing and transcription provides ample opportunities to regulate pA site recognition. It is therefore perhaps not surprising that alternative pA site usage is prevalent in the mammalian transcriptomes with ∼70% of all human genes subjected to what is commonly known as alternative cleavage and polyadenylation (APA).
APA occurs when a single gene has multiple pA sites, compared with constitutive polyadenylation, where one gene has one sole pA site. Like alternative splicing, APA results in multiple distinct RNA transcripts being produced from a single gene. The relative positions of these pA sites in a gene and their usage will ultimately determine the coding and regulatory sequence elements that are present in the different transcript isoforms. Depending on where the different pA sites are located in genes, APA events are separated into 2 major categories: untranslated region APA and coding region APA (Fig. 3).
Untranslated region APA (UTR-APA) occurs when multiple pA sites are available downstream of the terminal exon of a particular gene. As they are positioned outside the coding regions, differential usage of such pA sites does not alter the coding capacity of the resulting transcript isoforms but instead, alters the lengths of their 3′UTRs. The sequence that is absent in the short isoforms is commonly referred to as the alternative UTR (aUTR), relative to the common UTR (cUTR), which is present in both long and short APA isoforms. The presence or absence of an aUTR may equip a particular mRNA isoform with regulatory elements that influence the fate of the transcript regarding its stability, availability to the translational machinery and subsequent protein output, subcellular localization and even final protein destination. While aUTRs provide scope for potential regulation, it is not always clear what proportion of UTR-APA events have a physiologic impact.
Coding region APA (CR-APA) occurs when multiple pA sites are present at positions upstream of the 3′UTR. Thus CR-APA events can produce APA mRNA isoforms that differ in their coding capacity of the final protein product. Cryptic pA sites are also abundant in introns but, as mentioned above, they are actively suppressed by the U1snRNP.
In comparison to UTR-APA, it is easier to delineate a physiologic consequence for CR-APA events as they directly impact on the cellular proteome. Despite this, in the literature global analysis of APA is heavily bias toward UTR-APA.
Transcriptome wide APA profiles
The sum of all UTR-APA and CR-APA isoforms constitute the cellular APA profile. Advances in next generation sequencing have allowed cellular APA profiles to be analyzed on a transcriptome wide scale and several different pipelines have been developed and are available to investigate APA events (Table 1). UTR-APA in particular has attracted significant attention in the last decade and a large number of APA profiles have since been determined and are accessible in various repositories including: [http://www.polyasite.unibas.ch] [http://genome.bucm.edu.cn/utr/]. The analysis of these large data sets revealed that APA is highly dynamic and that the profiles can change depending on the context particular cells are studied in.
APA profiles in the context of proliferating and activated cells
A general pattern of shifts toward the proximal pA sites (UTR shortening) has been seen upon increased cellular proliferation. In the case of murine T-lymphocyte activation, the UTR-APA events show a clear pattern of a switch to proximal pA sites, compared with CR-APA, where there appeared to be equal movements toward both proximal and distal pA sites upon cell activation. Comparison with an extensive range of both murine and human cell lines revealed a similar pattern whereby the proliferation rate was inversely correlated with the 3′UTR length.
The activation of neurons by depolarising agents also induces shifts from the distal to proximal pA sites in a subset of genes that are regulated by the MEF2 family of transcription factors. Unlike the examples described above, many of these APA shifts are CR-APA events and thus result in the expression of functionally distinct proteins that may be critical for synapse development. The generality of APA profile shifts in response to cell activation is further shown by shifts that are apparent when astrocytes are exposed to eukaryotic growth factor EGF, T-cells are activated by interleukin-2 and B cells are cultured in the presence of anti-IgM and CpG. Interestingly there appears to be little overlap in genes that are undergoing stimulus dependent APA. This indicates that these shifts are stimulus and also cell type specific, but how they are achieved is unclear in most cases.
Similar to activation, stressed cells tend to enhance the use of intergenic pA sites and produce 3′extended transcripts. While no clear trend toward shortening or lengthening in UTR-APA has been observed, intronic pA sites are suppressed upon exposure to stress.
APA profiles in the context of development, tissue and cell types and different species
APA profiles have also been found to vary as a function of the differentiation status of cells. This was first described during mouse embryonic development, where the differentiation of C2C12 myoblast cells to myotubes showed a general shift to distal pA sites (UTR lengthening). Conversely, if differentiation is reversed, and induced pluripotent stem cells are generated from fully differentiated somatic cells, a general 3′UTR shortening is observed. Widespread CR-APA and UTR-APA have been reported to be critical during spermatogenesis. Genes associated with sperm maturation and testis specific genes compared with ubiquitously expressed genes tend to undergo 3′UTR shortening early in spermatogenesis avoiding destabilising and transposable elements that are located in the aUTRs of these genes.
When specific cell types or tissues are compared, they appear to have characteristic APA profiles. The most striking differences in APA profiles are those between mammalian and invertebrate testes and brain. While in testis the proximal pA sites are generally favored, in the brain the opposite is true and a preferred shift toward distal pA sites is characteristic. The prevalence of long 3′UTRs in neuronal tissue was reported using northern blot analysis over 20 y ago.
In addition to simple switches between common pA sites, there are several examples of unique pA sites that are only activated in specific tissues such as testis and the brain. In testis, these sites appear to have a lower incidence of the AAUAAA canonical hexamer, as well as having unique upstream and downstream elements. Since these sites are not used efficiently in somatic cells, it has been suggested that testis-specific pA factors may be responsible for this pA site switching, such as the CstF-64 variant, τCstF-64.
Similarly, tissue specific pA sites are also used in the brain where they generate specific APA mRNA isoforms with unusually long 3′UTRs ( >10 kb). These tissue specific APA events potentially provide a huge repertoire of diverse regulatory elements that can be plugged into tissue specific post-transcriptional regulatory networks.
The conservation of APA sites between species is limited and is restricted to fewer than 500 genes when tandem APA sites are compared between mouse and humans. If human pA sites are compared with rhesus, mouse, rat and dog, fewer than 1% are found in similar regions between them, indicating that APA is chiefly tissue specific, rather than species specific.
Experiments in yeast have compared the conservation of APA between different yeast species including Saccharomyces cerevisiae, Kluyveromyces lactis, and Debaryomyces hansenii. Here, despite the similarities in the general sequence preferences surrounding the pA sites, the APA profiles show species specificity. This led to the conclusion that much of the heterogeneity seen in the pA sites is likely to be biologic noise as defined by “biochemical events that occur in vivo as a result of a low-specificity process that has not been subjected to evolutionary optimisation.”
Recently, single cell studies have started to expose the heterogeneity of pA site usage within a population of homogenous stem cells. This revealed that individual cells differ significantly in their choice of pA site usage. Interestingly, it appears that the heterogeneity in pA site choice observed in these cells surpasses that of what would be expected if pA site choice was random, thereby concluding that this variability has the potential to contribute to the functional cell-to-cell heterogeneity.
APA profiles in the context of cancer
Several studies have observed shifts to usage of proximal pA sites upon cellular transformation and carcinogenesis. It is thought that shifts to these proximal pA sites in cancer function to exclude instability elements present in the aUTR, avoiding particular regulatory networks, such as those involving miRNA-mediated destabilisation of specific transcripts. Accordingly, in many cases shortening can result in increased mRNA and protein levels of oncogenes, for example in CCND1. However, it is noteworthy that in the cancer cells that show general 3′UTR shortening, there are also a significant number of genes that undergo 3′UTR lengthening. This clearly highlights that APA patterns in transformed cells are complex making the interpretation of the profiles in the context of the underlying mechanisms and the physiologic impact of the shifts more difficult.
How cancer specific APA profiles are established is largely unclear, but the dysregulation of CFIm in some cancer cells can be linked to 3′UTR shortening and provides a mechanistic explanation for this cancer-specific APA feature.
In addition to further the understanding of how cancer specific gene expression profiles are established, logging changes in APA profiles from tumor cells may also provide a valuable resource that can be exploited to develop tools to type cancer subtypes, provide information for diagnosis, and prognosis, and perhaps even treatment stratification. Indeed, APA profiles have previously been proposed to serve as biomarkers for cancer progression.
APA profiles in the context of other disease phenotypes
Many individual APA events have been linked with specific disease phenotypes such as thalassemia and IPEX syndrome. In these examples, pathogenesis is generally the result of mutations of cis-elements in one specific gene that compromises its 3′end processing efficacy. These mutations can result either in loss of function, such as in thalassemia and IPEX syndrome, or gain of function, as exemplified in thrombophilia and mantle cell lymphoma. In the latter, a point mutation in the CCND1 gene causes activation of a novel pA site leading to 3′UTR shortening and transcript and protein overexpression.
If proteins that regulate cleavage and polyadenylation are mutated, 3′end processing at many different pA sites can be altered and this can lead to changes in APA profiles. For example, the APA profiles between normal and failing hearts differ. While no clear global shifts in 3′UTR lengths are obvious, smaller cohorts of genes in the diseased compared with the normal heart show distinct changes in their APA profiles. Mechanistically the shifts in these APA profiles may be caused by the reduced expression levels of several cleavage and polyadenylation factors, including PABPN1. Interestingly, the levels of functional PABPN1 have before been linked to global shifts from distal to proximal pA sites in cells from oculopharyngeal muscular dystrophy (OPMD) patients. In OPMD cells the availability of functional PABPN1 is compromised by a triplet repeat expansion mutation of an alanine repeat in the PABPN1 gene. It has been suggested that in healthy cells, high levels of PABPN1 are available and can bind and suppress the proximal pA site in a significant number of genes. Conversely, in OPMD cells, the depletion of functional protein results in the activation of proximal pA sites causing 3′UTR shortening which may contribute to the pathology.
Myotonic dystrophy, a neuromuscular disease, represents an additional example where disease specific APA profiles are dependent on the availability of functional proteins. The Muscleblind-like (MBNL) protein family are key regulators of alternative splicing and mutations that compromise their function are associated with myotonic dystrophy. Mutations that cause loss of function in MBNL proteins also result in large scale bidirectional APA shifts.
Extensive APA defects that occur in the cerebellum of amyotrophic lateral sclerosis (ALS) patients have been reported. ALS patients that feature a repeat expansion in the C9orf72, have a higher frequency of short 3′UTR APA isoforms. The cause and the physiologic impact of these APA changes however is unknown.
For all these examples it is important to keep in mind that it is not clear whether the shifts in APA profiles are driving disease development or progression, or whether they are a reflection of the disease dependent changes in the cells.
Issues surrounding the global analysis and interpretation of APA profiles
As outlined above, APA profiles of cells appear to have characteristic signatures that are dynamic and can change in response to specific cues and disease. To elucidate the physiologic role of global APA changes, it is imperative to unravel the mechanisms and networks that are responsible for shaping these profiles. To do this, a distinction has to be drawn between alterations of APA profiles that are due to specific selection of one pA site over another at the point of cleavage, and changes in APA profiles that are due to post-transcriptional mechanisms, such as degradation of specific isoforms.
APA profile changes that occur as a result of cleavage site selection can be regarded as “active” APA (Fig. 4, left panel). Active APA can be achieved by changing the availability of factors that either enhance or suppress the usage of a particular pA site over another in response to a specific cue or change in the state of the cell. Conversely, profile changes that are the result of post-transcriptional events can be described as “passive” APA and may occur in response to a specific cue that alters the availability of factors, for example miRNA levels, that act post-transcriptionally (Fig. 4, right panel). This includes processes that alter the APA profile by interfering with the nuclear export rates or the stability of specific APA mRNA isoforms.
As most APA analysis is performed using whole-cell RNA, active and passive APA events are difficult to distinguish (Fig. 4). This may complicate the evaluation of particular profiles and may lead to misinterpretations of the underlying mechanisms and thus the physiologic relevance of APA. For example, shortening events in the context of cancer are often attributed to the selection of proximal sites at the point of cleavage to avoid post-transcriptional control executed by miRNAs. This interpretation however overlooks the fact that miRNA levels are known to be subjected to dramatic changes during oncogenic transformation and that different cancer subtypes have distinctive miRNA signatures. Thus a significant number of the shortening or lengthening events may be due to passive APA acting at the post-transcriptional level. Liaw et. al. highlight this possibility by cross-referencing miRNA profiles with previous APA data showing apparent 3′UTR shortening in the breast cancer cell line, MCF-7. This analysis showed that genes with target sites for upregulated miRNAs in their aUTR exhibited a larger degree of 3′UTR shortening relative to those without targets, and the extent of this was influenced by the number of target sites. In addition, the positioning of target sites for miRNAs upregulated in MCF-7 cells showed that they were enriched in the aUTRs of genes where 3′UTR shortening was seen, relative to downregulated miRNAs. Combining this data, it was concluded that selective degradation of solely the longer transcript isoforms by upregulated miRNAs has a significant impact that contributed to the observed changes in APA profiles in MCF-7 cells.
Perhaps many of these passive APA events in cancer cells are “collateral damage” caused by the upregulation or downregulation of particular miRNAs that aim to modulate a handful of specific genes. However, due to the nature of miRNA-mediated regulation, this will also affect many other genes that have little, if any, consequence for the progression or establishment of tumorigenesis.
The most common method used to untangle active and passive APA events is to look at the expression levels of all transcript isoforms. This approach argues that if observed shifts in APA patterns are due to a change in pA site choice, the increased proportion of the shorter transcripts leads to a lower overall decay rate, as 3′UTRs overall are of a repressive nature.
The expression level observed would, therefore, increase in this situation. Conversely, if the shift in APA pattern is due to an increase in the decay rate of solely the longer transcript, a decrease in expression would be seen. This method has its limitations, as any number of mechanisms, aside from APA, can influence transcript levels of a gene. To address this issue more directly, APA profiling of different subcellular fractionations have been used. In particular, the comparison between nuclear and cytoplasmic APA events can be used to get a better resolution on specific APA events and determine whether they are due to a change in cleavage site choice, or due to regulation at the post-transcriptional level including nuclear retention and transcript stability (Fig. 4). If subcellular APA is performed between different states of cells it may reveal how APA events are integrated with other regulatory networks that dynamically respond to changes in the state of cells.
Factors that regulate active APA by influencing pA sites choice at the point of cleavage
Despite the large number of reported shifts between pA sites in different cells and cell states, the precise mechanism of defining pA site choice in active APA remains ill defined. This is partly due to the complexity by which the strength of a pA site is defined and the many ways this can be modulated. In addition, as mentioned above, the difficulty to discriminate between active and passive APA further complicates the identification of factors. As the cis-elements influencing polyadenylation efficiency do not change, the switching between pA sites can only be achieved through the fluctuation of trans-factors or the presentation of pA sites to the processing machinery during transcription. The former is exemplified by the transplant of large chromosomal segments from Debaryomyces hansenii into Saccharomyces cerevisiae and the concomitant adoption of the host APA profile by the recipient. Several trans-factors have been shown to influence pA site choice (Fig. 5), including core pA factors (Table 2) and several RNA-binding proteins of diverse function (Table 3).
As can be seen from Tables 2, 3 and 4, there is a plethora of trans-factors that, when depleted or compromised, cause shifts in APA profiles.
A number of RNA binding proteins have been associated with influencing pA site usage, that primarily have other nuclear functions, as outlined in Table 3. Many of these mediate splicing, again highlighting the interconnection between these 2 pre-mRNA processing mechanisms.
As summarised in Table 4, several features that modulate gene transcription have also been linked to active APA. It has been long known that pre-mRNA processing and transcription are tightly coupled. Cleavage and polyadenylation are essential for transcription termination, as well as the recycling of transcriptional factors and the re-initiation at upstream promoters. The phosphorylated CTD of Pol II also provides an essential platform for the recruitment of several pA factors to the pre-mRNA. Specifically, CPSF and CstF, together with Pcf11, are brought to the pre-mRNA through the interaction with Pol II. The pausing of Pol II during transcription, and the subsequent recruitment of CstF to the transcription elongation complex, are decisive in pA site choice. Even after the Pol II has passed the pA site, the RNA tether between the 2 is necessary for the assembly of the mature protein complex, to perform efficient cleavage and polyadenylation.
It is therefore not surprising that APA is also governed by parameters that compromise the kinetics of transcription, including Pol II transcriptional rates, chromatin structure, and histone modifications. This is beautifully exemplified in mutant Drosophila strains that reduce transcription elongation rates of Pol II which results in the preferential use of proximal pA sites in several alternatively polyadenylated genes. This proposed mechanism is analogous to the one proposed for exon skipping, whereby a faster Pol II reduces the latency period between when the proximal and distal pA sites are transcribed, thus reducing the probability of using the generally weaker proximal pA site. Furthermore, the conformation of chromatin has been shown to influence APA. Studies looking at nucleosome positioning surrounding pA sites show a general depletion of nucleosomes immediately downstream of high-usage pA sites, but then an enrichment past ∼100 bp downstream of the pA site, when compared with low-usage pA sites. The influence of transcription on cleavage and polyadenylation is not limited to the 3′end of the gene. The presence or absence of specific transcriptional activators can stimulate transcription-coupled cleavage and polyadenylation in yeast through Paf1c, which is thought to facilitate the formation of the pA factors in the elongation complex. In mammalian cells, the Paf1 complex is also implicated in pA site regulation. Here, the depletion of some subunits (Paf1, Cdc73, Ski8) caused pervasive transcript shortening of UTRs, and reduction in Paf1c or Cdc73 also increased usage of pA site located in upstream introns/exons. It appears that Paf1c depletion affects Pol II progression through the gene affecting pA site selection. Absence of Paf1c function causes slowing of the polymerase which would favor the more proximal positioned sites. Furthermore, the promoters of specific genes in Drosophila, have been shown to be essential in the recruitment of ELAV to Pol II, which subsequently acts at the 3′ends of those genes to inhibit proximal pA site, resulting in 3′UTR extension. The upregulation of specific transcription factors in response to cell activation, for example E2F, can influence APA profiles indirectly by increasing the gene expression of key 3′end processing factors.
The impact of active APA on the regulation of gene expression
Modulating the length of the UTRs by APA has the potential to produce mRNA isoforms that can be subjected to distinct post-transcriptional regulation (Fig. 6). As 3′UTRs are thought to have an overall repressive role, cleavage and polyadenylation of a transcript at a pA site more proximal to the stop codon is generally considered to result in the removal of potentially repressive cis-elements in the 3′UTR, summarised in Table 5. This can, therefore, stabilize the transcript, or promote its translation, thus increasing the overall expression of that gene (Fig. 6). However, several RNA binding proteins have been identified to bind to the 3′UTR and stabilize the transcript, for example the Hu family of proteins (including HuR), which bind to AREs. The neuronal specific HuD protein stabilizes transcripts required in neurite outgrowth enabling neuronal differentiation in several cell lines.
The impact of APA on transcript isoform stability and translation
Global comparisons of cis-elements located in the aUTRs compared with cUTRs revealed a bias toward a higher representation of conserved miRNA seed regions in aUTRs. The miRNA target sites in aUTRs are also located in regions with higher AU content. This reduces the possibility that these target sites are involved in secondary structures and thus would represent better targets for miRNAs.
Individual gene examples such as PAX3, CCND1, CCND2, CDC6, MAPT and MGMT support that the biased distribution of miRNA binding sites in APA mRNA isoforms can act as a potent layer to govern miRNA mediated regulation of gene expression. Similarly, APA in ELAVL1 and PDCL regulates their expression by controlling the presence or absence of destabilising or stabilizing RNA binding protein recognition motifs (Table 6).
However, despite these individual gene examples and the general biased distribution of cis-elements, the actual global impact of APA on isoform stability and transcript levels has been found to be fairly modest. This was further supported by more detailed analyses addressing the impact of UTR shortening on mRNA isoform half-life times in proliferating T-cells and mouse 3T3 cells, where aUTRs showed a limited influence on the stability of APA mRNA isoforms. In addition, contrary to expectations, a significant proportion of APA mRNA isoforms with extended UTRs have been shown to have increased stability in 3T3 cells, a phenomenon that was also observed in adipocyte stem cells. Following this trend, a recent analysis of nuclear and cytoplasmic fractions confirmed that the impact of UTR-APA on the stability is modest. In HEK293 cells, around 10 percent of all APA events are found to be subject to post-transcriptional regulation and around 3% of all those cytoplasmic events are destabilised by miRNAs. Interestingly, this study further showed that miRNAs target both short and long UTR-APA isoforms for degradation.
The lack of a global impact of APA on mRNA stability in mammalian cells however is not universal. Contrary to mammals, in yeast, changes in UTR length have comprehensive implications on transcript stability and even changes at a single nucleotide level can have a significant large scale impact on transcript stability.
The impact of different APA mRNA isoforms on translation efficiency follows a similar pattern. While the APA mediated shortening of UTRs in several mRNAs including CDC6, HSPA2, ECE-1, CCND2 and DICER1 was concomitant with a higher translational output, the opposite was true for many genes in 3T3 cells. Most notably in the SERT gene (serotonin receptor), the distal isoform is translated at a higher rate. As for the impact on stability, APA controlled UTR shortening appears to have a limited global impact on translation efficiency.
While there is not a clear global trend, it is nevertheless clear from these examples that the differential usage of pA sites can alter the stability of selected transcripts and affect their overall protein expression, consequently having a significant impact on major biochemical processes.
UTR-APA and the control of subcellular localization of transcripts
UTRs can play an important role in directing the subcellular localization of transcripts and it is therefore not surprising that UTR-APA influences subcellular distribution of mRNA isoforms. 3′UTR lengthening in the brain is one of the most striking examples of tissue-specific APA and this also appears to provide an important platform to control the localization of APA mRNA isoforms. BDNF, Ranbp1, Impa1, MKK7 and KPNB1 are examples where APA isoforms with long 3′UTRs include localization signals resulting in targeted transport and translation of the respective transcripts. Similarly, 3′UTR extension in α-synuclein transcripts, which are associated with Parkinson disease pathology, not only increases the translational output but also affects their localization away from the synaptic-terminals to mitochondria. On a more global scale, a recent comparison between RNA isolated from the dendrites and the stroma in neurons found UTR-APA isoforms that display differential subcellular localization. Distal APA isoforms that are induced during neuronal differentiation are more likely to be localized in neurite projections than their proximal APA counterparts.
APA mediated subcellular localization is not just restricted to highly polarized cells. At the individual gene level, it has been shown that the presence of Alu repeats in the aUTRs of the Nicolin1 gene is incompatible with nuclear cytoplasmic export in HEK293 cells. In a more recent global study, a significant number of nuclear retained APA mRNA isoforms have been identified when nuclear and cytoplasmic APA profiles in HEK293 cells were compared. At least for some of these APA isoforms incomplete splicing resulting in intron retention is instrumental to prevent export into the cytoplasm. In addition, APA isoforms with short 3′UTRs tend to be overrepresented in the cytoplasm compared with the nucleus and while this phenomenon is conserved between different cell types, the genes affected are cell type specific. Interestingly, 2 prominent splicing factors, SRSF3 and SRSF7 have recently been linked to both pA site selection and selective nuclear cytoplasmic transport. SRSF3 promotes the usage of distal pA sites and by recruiting the export factor NXF1, controls the export of long 3′UTR APA isoforms. In contrast, SRSF7 causes a shift toward proximal pA site usage and subsequently supports the export of the resulting short UTR-APA isoforms.
UTR-APA controlling localized translation has also been seen with some membrane bound proteins including CD47. The distal APA isoform of CD47 allows binding of HuR and SET resulting in the translated protein to be localized to the plasma membrane. Conversely, the proximal APA isoform translated CD47 locates primarily to the endoplasmic reticulum.
The impact of CR-APA
In contrast to UTR-APA, CR-APA by definition, will affect the protein identity, as introducing a pA site upstream of a translational stop codon will result in a protein with a different C-terminus. This may be a simple truncation in the case of exonic CR-APA, or more commonly a different C-terminus in composite-intronic CR-APA or skipped-intronic CR-APA (Fig. 1). In the case of skipped-intronic APA, polyadenylation is coupled to an alternative splicing event of the alternative terminal intron and is therefore in dynamic competition with splicing. Thousands of these intronic pA sites have been identified, but largely remain dormant through inhibition by factors such as U1 snRNP (see above). However, these intronic pA sites can indeed be activated, as is seen during increased proliferation. In around 7.5% of breast cancer patients, a short version of MAG1I3 as a result of premature polyadenylation at a cryptic pA site located in intron 10 of the gene is present. The shortened MAHI3 protein product promotes mammary cell transformation. However, the activation of this cryptic pA site is unclear as it is not dependent on mutations of local cis-elements and thus is more likely to involve an imbalance of trans-factors.
The classic example of physiologic relevant CR-APA is that of the switch between the membrane-bound form of IgM in B-cells to the secreted form in plasma cells. This switch is thought to be governed by changes in the CstF-64 levels, whereby an elevation of CstF-64 in plasma cells activates the weaker proximal pA site, triggering the switch to the upstream pA site, ultimately removing the membrane anchoring C-terminus.
One important constraint of CR-APA isoforms is that they do not result in the introduction of premature translational termination codons. Therefore, they are not subject to nonsense-mediated decay and have the capability of producing alternative protein isoforms with potentially distinct physiologic roles. CR-APA events also include examples where the polyadenylation event within the coding region of an exon can convert a tyrosine codon (TAT) into a stop codon (TAA). This mechanism regulates expression of a truncated form of glutamyl-prolyl tRNA synthetase that enables the escape of inflammatory genes from the GAIT complex mediated translational silencing.
Outlook
The advent of high-throughput sequencing technologies has no doubt made a huge contribution to our understanding of 3′end formation and APA. It has allowed us to map pA sites on a transcriptome wide scale and compare their differential usage in tissues and diverse cellular states. While in the past the focus was on mapping “global” trends such as 3′UTR lengthening and shortening in different cellular contexts, in the future the focus has to shift more toward identifying APA events that have a proven physiologic impact. This requires a clear distinction to be made between active and passive APA events. Furthermore, rather than just bioinformatics correlations, the physiologic consequence of such APA events also has to be confirmed by direct experimentation at the individual gene level.
The appreciation of APA as a global regulator of gene expression in recent years has also revived a renewed focus on the mechanistic side of pA site recognition. Despite the simplicity of the reaction, it is still unclear how the pA machinery can assemble into functional complexes at pA signals that deviate considerably from the consensus sequences and differ widely in their architecture. Developing methods that can assess the likelihood of a pA site being subjected to regulation and further dissecting the mechanistic details of pA site recognition will also be critical to fully understand APA and its physiologic impact.
References 277
- 1.
Evguenieva-Hackenberg E, Klug G New aspects of RNA processing in prokaryotes. Curr Opin Microbiol 2011; 14:587-92; PMID:21945217; https://doi.org/
Not in our library10.1016/j.mib.2011.07.025
- 2.
Proudfoot NJ, Furger A, Dye MJ Integrating mRNA processing with transcription. Cell 2002; 108:501-12; https://doi.org/
Not in our library10.1016/s0092-8674(02)00617-7
- 3.
Gruber AJ, Schmidt R, Gruber AR, Martin G, Ghosh S, Belmadani M, Keller W, Zavolan M A comprehensive analysis of 3′ end sequencing data sets reveals novel polyadenylation signals and the repressive role of heterogeneous ribonucleoprotein C on cleavage and polyadenylation. Genome Res 2016; 26:1145-59; PMID:27382025; https://doi.org/
Not in our library10.1101/gr.202432.115
- 4.
Beaudoing E, Freier S, Wyatt JR, Claverie JM, Gautheret D Patterns of variant polyadenylation signal usage in human genes. Genome Res 2000; 10:1001-10; PMID:10899149; https://doi.org/
Not in our library10.1101/gr.10.7.1001
- 5.
Tian B, Hu J, Zhang H, Lutz CS A large-scale analysis of mRNA polyadenylation of human and mouse genes. Nucleic Acids Res 2005; 33:201-12; PMID:15647503; https://doi.org/
Not in our library10.1093/nar/gki158
- 6.
Zarudnaya MI, Kolomiets IM, Potyahaylo AL, Hovorun DM Downstream elements of mammalian pre-mRNA polyadenylation signals: Primary, secondary and higher-order structures. Nucleic Acids Res 2003; 31:1375-86; PMID:12595544; https://doi.org/
Not in our library10.1093/nar/gkg241
- 7.
Hu J, Lutz CS, Wilusz J, Tian B Bioinformatic identification of candidate cis-regulatory elements involved in human mRNA polyadenylation. RNA 2005; 11:1485-93; PMID:16131587; https://doi.org/
Not in our library10.1261/rna.2107305
- 8.
Fitzgerald M, Shenk T The sequence 5′-AAUAAA-3′forms parts of the recognition site for polyadenylation of late SV40 mRNAs. Cell 1981; 24:251-60; PMID:6113054; https://doi.org/
Not in our library10.1016/0092-8674(81)90521-3
- 9.
Montell C, Fisher EF, Caruthers MH, Berk AJ Inhibition of RNA cleavage but not polyadenylation by a point mutation in mRNA 3′ consensus sequence AAUAAA. Nature 1983; 305:600-5; PMID:6194440; https://doi.org/
Not in our library10.1038/305600a0
- 10.
Wickens M, Stephenson P Role of the conserved AAUAAA sequence: Four AAUAAA point mutants prevent messenger RNA 3′ end formation. Science (80-) 1984; 226:1045-51; PMID:6208611; https://doi.org/
Not in our library10.1126/science.6208611
- 11.
Orkin SH, Cheng TC, Antonarakis SE, Kazazian HH Thalassemia due to a mutation in the cleavage-polyadenylation signal of the human beta-globin gene. EMBO J 1985; 4:453-6; PMID:4018033
Not in our library10.1002/j.1460-2075.1985.tb03650.x
- 12.
Chen J-M, Férec C, Cooper D A systematic analysis of disease-associated variants in the 3′ regulatory regions of human protein-coding genes I: General principles and overview. Hum Genet 2006; 120:1-21; PMID:16645853; https://doi.org/
Not in our library10.1007/s00439-006-0180-7
- 13.
Danckwardt S, Hentze MW, Kulozik AE 3′ end mRNA processing: Molecular mechanisms and implications for health and disease. Embo J 2008; 27:482-98; PMID:18256699; https://doi.org/ Available from http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&dopt=Citation&list_uids=18256699
Not in our library10.1038/sj.emboj.7601932
- 14.
Higgs DR, Goodbourn SE, Lamb J, Clegg JB, Weatherall DJ, Proudfoot NJ Alpha-thalassaemia caused by a polyadenylation signal mutation. Nature 1983; 306:398-400; PMID:6646217 Available from http://www.ncbi.nlm.nih.gov/pubmed/6646217
Not in our library10.1038/306398a0
- 15.
Ferreira PG, Oti M, Barann M, Wieland T, Ezquina S, Friedländer MR, Rivas MA, Esteve-Codina A, Rosenstiel P, Strom TM, et al. Sequence variation between 462 human individuals fine-tunes functional sites of RNA processing. Sci Rep 2016; 6:32406; PMID:27617755; https://doi.org/
Not in our library10.1038/srep32406
- 16.
Neve J, Burger K, Li W, Hoque M, Patel R, Tian B, Gullerova M, Furger A Subcellular RNA profiling links splicing and nuclear DICER1 to alternative cleavage and polyadenylation. Genome Res 2016; 26:24-35; PMID:26546131; https://doi.org/ Available from http://genome.cshlp.org/lookup/doi/10.1101/gr.193995.115
Not in our library10.1101/gr.193995.115
- 17.
Zhao J, Hyman L, Moore C Formation of mRNA 3′ ends in eukaryotes: Mechanism, regulation, and interrelationships with other steps in mRNA synthesis. Microbiol Mol Biol Rev 1999; 63:405-45; PMID:10357856 Available from http://www.ncbi.nlm.nih.gov/pubmed/10357856
Not in our library10.1128/mmbr.63.2.405-445.1999
- 18.
Nunes NM, Li W, Tian B, Furger A A functional human poly(A) site requires only a potent DSE and an A-rich upstream sequence. EMBO J 2010; 29:1523-36; PMID:20339349; https://doi.org/
Not in our library10.1038/emboj.2010.42
- 19.
Danckwardt S, Gehring NH, Neu-Yilik G, Hundsdoerfer P, Pforsich M, Frede U, Hentze MW, Kulozik AE The prothrombin 3′end formation signal reveals a unique architecture that is sensitive to thrombophilic gain-of-function mutations. Blood 2004; 104:428-35; PMID:15059842; https://doi.org/ Available from http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&dopt=Citation&list_uids=15059842
Not in our library10.1182/blood-2003-08-2894
- 20.
Hwang H-W, Park CY, Goodarzi H, Fak JJ, Mele A, Moore MJ, Saito Y, Darnell RB PAPERCLIP identifies MicroRNA targets and a role of CstF64/64tau in promoting non-canonical poly(A) site usage. Cell Rep 2016; 15:423-35; PMID:27050522; https://doi.org/ Available from http://www.sciencedirect.com/science/article/pii/S2211124716302649
Not in our library10.1016/j.celrep.2016.03.023
- 21.
Kandala DT, Mohan N, A V, AP S, G R, Laishram RS CstF-64 and 3′-UTR cis-element determine Star-PAP specificity for target mRNA selection by excluding PAPα. Nucleic Acids Res 2016; 44:811-23; PMID:26496945; https://doi.org/ Available from http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4737136/
Not in our library10.1093/nar/gkv1074
- 22.
Li X-Q, Du D RNA polyadenylation sites on the genomes of microorganisms, animals, and plants. PLoS One 2013; 8:e79511; PMID:24260238; https://doi.org/ Available from https://doi.org/10.1371%252Fjournal.pone.0079511
Not in our library10.1371/journal.pone.0079511
- 23.
Gehring NH, Frede U, Neu-Yilik G, Hundsdoerfer P, Vetter B, Hentze MW, Kulozik AE Increased efficiency of mRNA 3[prime] end formation: A new genetic mechanism contributing to hereditary thrombophilia. Nat Genet 2001; 28:389-92; PMID:11443298; https://doi.org/
Not in our library10.1038/ng578
- 24.
Dalziel M, Nunes NM, Furger A Two G-rich regulatory elements located adjacent to and 440 nucleotides downstream of the core poly(A) site of the intronless melanocortin receptor 1 gene are critical for efficient 3′ end processing. Mol Cell Biol 2007; 27:1568-80; PMID:17189425; https://doi.org/
Not in our library10.1128/mcb.01821-06
- 25.
Arhin GK, Boots M, Bagga PS, Milcarek C, Wilusz J Downstream sequence elements with different affinities for the hnRNP H/H’ protein influence the processing efficiency of mammalian polyadenylation signals. Nucleic Acids Res 2002; 30:1842-50; PMID:11937639 Available from http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&dopt=Citation&list_uids=11937639
Not in our library10.1093/nar/30.8.1842
- 26.
Dalziel M, Kolesnichenko M, das Neves RP, Iborra F, Goding C, Furger A {alpha}-MSH regulates intergenic splicing of MC1R and TUBB3 in human melanocytes. Nucleic Acids Res 2011; 39:2378-92; PMID:21071418; https://doi.org/ Available from http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&dopt=Citation&list_uids=21071418
Not in our library10.1093/nar/gkq1125
- 27.
Oberg D, Fay J, Lambkin H, Schwartz S A downstream polyadenylation element in human papillomavirus type 16 L2 encodes multiple GGG motifs and interacts with hnRNP H. J Virol 2005; 79:9254-69; PMID:15994820; https://doi.org/ Available from http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&dopt=Citation&list_uids=15994820
Not in our library10.1128/jvi.79.14.9254-9269.2005
- 28.
Ozsolak F, Kapranov P, Foissac S, Kim SW, Fishilevich E, Monaghan AP, John B, Milos PM Comprehensive polyadenylation site maps in yeast and human reveal pervasive alternative polyadenylation. Cell 2011; 143:1018-29. Available from http://linkinghub.elsevier.com/retrieve/pii/S0092867410013000
Not in our library10.1016/j.cell.2010.11.020
- 29.
Kaufmann I, Martin G, Friedlein A, Langen H, Keller W Human Fip1 is a subunit of CPSF that binds to U-rich RNA elements and stimulates poly(A) polymerase. Embo J 2004; 23:616-26; PMID:14749727; https://doi.org/ Available from http://www.ncbi.nlm.nih.gov/pubmed/14749727
Not in our library10.1038/sj.emboj.7600070
- 30.
Lackford B, Yao C, Charles GM, Weng L, Zheng X, Choi EA, Xie X, Wan J, Xing Y, Freudenberg JM, et al. Fip1 regulates mRNA alternative polyadenylation to promote stem cell self-renewal. Embo J 2014; 33:878-89; PMID:24596251; https://doi.org/ Available from http://www.ncbi.nlm.nih.gov/pubmed/24596251
Not in our library10.1002/embj.201386537
- 31.
- 32.
Schek N, Cooke C, Alwine JC Definition of the upstream efficiency element of the simian virus 40 late polyadenylation signal by using in vitro analyses. Mol Cell Biol 1992; 12:5386-93; PMID:1333042 Available from http://mcb.asm.org/content/12/12/5386.abstract
Not in our library10.1128/mcb.12.12.5386-5393.1992
- 33.
Rothnie HM, Reid J, Hohn T The contribution of AAUAAA and the upstream element UUUGUA to the efficiency of mRNA 3′-end formation in plants. EMBO J 1994; 13:2200-10; PMID:8187773 Available from http://www.ncbi.nlm.nih.gov/pmc/articles/PMC395075/
Not in our library10.1002/j.1460-2075.1994.tb06497.x
- 34.
Valsamakis A, Schek N, Alwine JC Elements upstream of the AAUAAA within the human immunodeficiency virus polyadenylation signal are required for efficient polyadenylation in vitro. Mol Cell Biol 1992; 12:3699-705; PMID:1508176 Available from http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&dopt=Citation&list_uids=1508176
Not in our library10.1128/mcb.12.9.3699-3705.1992
- 35.
Gilmartin GM, Fleming ES, Oetjen J, Graveley BR CPSF recognition of an HIV-1 mRNA 3′-processing enhancer: Multiple sequence contacts involved in poly(A) site definition. Genes Dev 1995; 9:72-83; PMID:7828853 Available from http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&dopt=Citation&list_uids=7828853
Not in our library10.1101/gad.9.1.72
- 36.
Brackenridge S, Proudfoot NJ Recruitment of a basal polyadenylation factor by the upstream sequence element of the human lamin B2 polyadenylation signal. Mol Cell Biol 2000; 20:2660-9; PMID:10733568 Available from http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&dopt=Citation&list_uids=10733568
Not in our library10.1128/mcb.20.8.2660-2669.2000
- 37.
Moreira A, Takagaki Y, Brackenridge S, Wollerton M, Manley JL, Proudfoot NJ The upstream sequence element of the C2 complement poly(A) signal activates mRNA 3′ end formation by two distinct mechanisms. Genes Dev 1998; 12:2522-34; PMID:9716405 Available from http://www.ncbi.nlm.nih.gov/pubmed/9716405
Not in our library10.1101/gad.12.16.2522
- 38.
Danckwardt S, Kaufmann I, Gentzel M, Foerstner KU, Gantzert AS, Gehring NH, Neu-Yilik G, Bork P, Keller W, Wilm M, et al. Splicing factors stimulate polyadenylation via USEs at non-canonical 3′ end formation signals. Embo J 2007; 26:2658-69; PMID:17464285; https://doi.org/ Available from http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&dopt=Citation&list_uids=17464285
Not in our library10.1038/sj.emboj.7601699
- 39.
Natalizio BJ, Muniz LC, Arhin GK, Wilusz J, Lutz CS Upstream elements present in the 3′-untranslated region of collagen genes influence the processing efficiency of overlapping polyadenylation signals. J Biol Chem 2002; 277:42733-40; PMID:12200454; https://doi.org/ Available from http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&dopt=Citation&list_uids=12200454
Not in our library10.1074/jbc.m208070200
- 40.
Hall-Pogar T, Zhang H, Tian B, Lutz CS Alternative polyadenylation of cyclooxygenase-2. Nucleic Acids Res 2005; 33:2565-79. Available from http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&dopt=Citation&list_uids=15872218
Not in our library10.1093/nar/gki544
Reads
Read by people on this page
None yet.
Checked or read by AI through naulon
None yet.
A person counts once a day per paper. We keep no cookie and no address: a reader is a hash that changes every day. Any month under five reads shows no number.
Are you an author?
Is this your paper?
Claim it to add your note and be paid when agents cite it. Free.