#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Systematic Profiling of Poly(A)+ Transcripts Modulated by Core 3’ End Processing and Splicing Factors Reveals Regulatory Rules of Alternative Cleavage and Polyadenylation


A gene can express multiple isoforms varying in the 3’ end, a phenomenon called alternative cleavage and polyadenylation, or APA. Previous studies have indicated that most eukaryotic genes display APA and the APA profile changes under different physiological and pathological conditions. However, how APA is regulated in the cell is unclear. Here using gene knockdown and high throughput sequencing we examine how APA is regulated by factors in the machinery responsible for cleavage and polyadenylation as well as factors that play essential roles in splicing. We identify several factors that play significant roles in APA in the last exon, including CFI-25/68, PABPN1, PABPC1, Fip1 and Pcf11. We also elucidate how cleavage and polyadenylation events are regulated in introns and near the transcription start site. We uncover a group of APA events that are highly regulated by core factors as well as in cell differentiation and development. We present an APA code where an APA event in a given cellular context is regulated by a number of parameters, including relative location to the transcription start site, splicing context, distance between competing pAs, surrounding cis elements and concentrations of core cleavage and polyadenylation factors.


Published in the journal: . PLoS Genet 11(4): e32767. doi:10.1371/journal.pgen.1005166
Category: Research Article
doi: https://doi.org/10.1371/journal.pgen.1005166

Summary

A gene can express multiple isoforms varying in the 3’ end, a phenomenon called alternative cleavage and polyadenylation, or APA. Previous studies have indicated that most eukaryotic genes display APA and the APA profile changes under different physiological and pathological conditions. However, how APA is regulated in the cell is unclear. Here using gene knockdown and high throughput sequencing we examine how APA is regulated by factors in the machinery responsible for cleavage and polyadenylation as well as factors that play essential roles in splicing. We identify several factors that play significant roles in APA in the last exon, including CFI-25/68, PABPN1, PABPC1, Fip1 and Pcf11. We also elucidate how cleavage and polyadenylation events are regulated in introns and near the transcription start site. We uncover a group of APA events that are highly regulated by core factors as well as in cell differentiation and development. We present an APA code where an APA event in a given cellular context is regulated by a number of parameters, including relative location to the transcription start site, splicing context, distance between competing pAs, surrounding cis elements and concentrations of core cleavage and polyadenylation factors.

Introduction

Pre-mRNA cleavage and polyadenylation (C/P) is a 3’ end processing mechanism employed in eukaryotes for expression of almost all protein-coding transcripts and long non-coding RNAs by RNA polymerase II (RNAPII)[1, 2]. The site for C/P, commonly known as the polyA site or pA, is defined by both upstream and downstream cis elements [3, 4]. As with core RNAPII promoters [5], core C/P signals are proving to be complex. In mammals, upstream elements include the polyadenylation signal (PAS), such as AAUAAA, AUUAAA, or close variants, located within ~40 nucleotides (nt) from the pA; UGUA elements, typically located upstream of the PAS; and U-rich elements located around the PAS. Downstream elements include U- and GU-rich elements, located within ~100 nt downstream of the pA.

Most mammalian genes express alternative cleavage and polyadenylation (APA) isoforms [69]. APA in the 3’ untranslated region (3’UTR), called 3’UTR-APA, leads to isoforms with different 3’UTR lengths. Because the 3’UTR plays an important role in aspects of mRNA metabolism, such as subcellular localization, stability, and translation, 3’UTR-APA can impact the post-transcriptional control of gene expression. In addition, about 40% of mammalian genes display APA in upstream introns and internal exons [10], leading to changes of both coding sequences (CDSs) and 3’UTRs. This type of APA is called CDS-APA herein. Studies have shown that the APA pattern of genes is tissue-specific [7, 9, 11, 12], and is regulated under various conditions, such as cell proliferation, differentiation and development [8, 1316] and response to extracellular signals [17]. However, despite recent advances, the molecular mechanisms that regulate APA are still poorly understood (see [4, 1820] for reviews).

The C/P machinery in mammalian cells is composed of over 20 core factors [21]. Some form subcomplexes, including the Cleavage and Polyadenylation Specificity Factor (CPSF), containing CPSF-160, CPSF-100, CPSF-73, CPSF-30, Fip1, and WDR33; the Cleavage stimulation Factor (CstF), containing CstF-77, CstF-64, and CstF-50; the Cleavage Factor I (CFI), containing CFI-68 or CFI-59 and CFI-25; and the Cleavage Factor II (CFII), containing Pcf11 and Clp1. Single proteins involved in C/P include Symplekin, poly(A) polymerase (PAP), nuclear poly(A) binding protein (PABPN), RBBP6, and RNAPII (specifically the C-terminal domain of its largest subunit). In addition, protein phosphatase 1α (PP1α) and protein phosphatase 1β (PP1β) are present in the C/P complex and homologous to yeast C/P factors [22], but their functions in 3’ end processing are yet to be established. Since the initial study indicating that CstF-64 can regulate APA of the IgM heavy chain gene pre-mRNA during B cell maturation [23], a growing number of core C/P factors have been shown to impact pA choice, such as CFI factors [2426], PABPN1 [27, 28], CstF-64τ [29, 30], Fip1 [31] and RBBP6 [32]. Whether other C/P factors are involved in APA is not known, and to what extent different C/P factors differentially modulate APA has not been systematically examined.

C/P can also be regulated by splicing (reviewed in [33]), which has long been thought to help define the 3’ terminal exon [34]. U1 snRNP (or U1) has been shown to suppress cryptic pA usage near the transcription start site (TSS)[35], which may be attributable to its inhibitory activity on poly(A) polymerase α (PAPα)[36]. The interplay between U1 and C/P has recently been implicated in controlling expression of sense vs. antisense transcripts from bidirectional promoters [37, 38]. Interestingly, mild attenuation of U1 that does not inhibit splicing was found to regulate 3’UTR length via APA [39]. This mechanism, dubbed telescripting, has been implicated in pre-mRNA shortening during rapid and transient transcriptional upregulation in neuronal activation when there is shortage of U1. By contrast, no global regulation of APA has been associated with U2 snRNP (or U2), despite various interactions between U2 factors and the C/P machinery [40, 41].

Here using C2C12 cells and siRNA-based knockdown (KD) of individual core C/P factors, we examine the role of C/P factors in global APA regulation. We also reduce the expression and/or activity of several core splicing factors in order to understand the role of splicing in APA. In addition, we investigate expression of sense and antisense transcripts using pAs near the TSS when C/P or splicing factors are inhibited. Lastly, we examine how APA events regulated by C/P factor KDs are related to those taking place in cell differentiation. Altogether, our results reveal multiple regulatory rules of APA and show that modulation of core RNA processing factor levels provides a powerful mechanism to control global APA and that different factors have distinct impacts on pA usage.

Results

Global analysis of APA upon inhibition of core C/P and splicing factors

To examine the role of core C/P factors in APA regulation, we set out to knock down by siRNA the expression of each of the known factors (Fig 1 and S1 Table, see Introduction), including CPSF-160, CPSF-100, CPSF-73, CPSF-30, Fip1, WDR33, RBBP6, CstF-77, CstF-64, CstF-64τ, CstF-50, CFI-25, CFI-59, CFI-68, Pcf11, Clp1, PAPα, PAPγ, PABPN1, and Symplekin. In addition, we designed siRNAs for several C/P complex-associated factors [22], including PP1α, PP1β, and poly(A) binding protein C1 (PABPC1). Since a large fraction of pAs are located in introns and are regulated in development and differentiation [10], we also included siRNAs for several core splicing factors that have been previously implicated in regulation of C/P through various mechanisms [33], including U1 factor U1-70K, U2 factor SF3b155 and U2-associated factor U2AF65, Moreover, because functional inhibition of U1 was previously shown to cause global regulation of APA in both introns [35] and 3’UTRs [39], we included in our experiments an oligonucleotide (oligo) that can base pair with the 5’ end region of U1 snRNA involved in 5’ splice site (5’SS) recognition (Fig 1B). This oligo, named U1 domain oligo (U1D), is composed of locked nucleic acids and 2’-O-methylated nucleotides, and has been shown to base pair efficiently with U1 snRNA, thereby functionally knocking down U1 [42]. An oligo that had the same chemical composition but contained two mutations at critical positions for base pairing with U1 snRNA, named mutant U1D oligo (mU1D, Fig 1B), was used as a control.

Fig. 1. Systematic analysis of APA events modulated by C/P and splicing factors.
Systematic analysis of APA events modulated by C/P and splicing factors.
(A) Experimental design. Proliferating C2C12 myoblast cells were transfected with siRNAs against 23 C/P factors and 3 splicing factors for 48 hr, followed by APA analysis using 3’ region extraction and deep sequencing (3’READS). Sequences of the siRNAs used for knockdown or control are listed in S1 Table. For direct inhibition of U1 snRNP (U1) activity, C2C12 cells were transfected with U1D or mutant U1D (mU1D) oligonucleotides (oligos, shown in (B)) for 8 hr or 24 hr, followed by 3’READS analysis. (B) Oligos used to study U1 functions. Top, the 5’ end region of U1 snRNA and a consensus sequence surrounding 5’ splice site (5’SS, -3 nt to +6 nt, based on all annotated 5’SS in the mouse genome), are shown to illustrate the base-pairing between U1 snRNA and 5’SS sequence; Bottom, U1D or mutant U1D (mU1D) sequences. Locked nucleic acids are in upper case and 2'-O-methyl nucleotide in lower case. Nucleotides shown in red correspond to the -3 to +6 region surrounding the 5’SS. mU1D differs from U1D in two nucleotides as indicated. (C) Gene expression of different factors after 48 hr of siRNA knockdown (KD), as measured by reverse transcription-quantitative PCR (RT-qPCR) or Western Blot (WB). Error bars for RT-qPCR are standard deviation based on two technical replicates. Some KD samples were not analyzed by WB, as indicated. See S1 Fig for WB images.

We chose C2C12 myoblast cells for this study because we previously found that during differentiation of C2C12 cells into myotubes, APA isoforms using proximal pAs (relative to the 5’ end of gene), including those in 3’UTRs and introns, are relatively downregulated, while distal pAs are relatively upregulated [10, 14]. Therefore, by using C2C12 cells, we could compare APA events regulated by C/P or splicing factors with those taking place naturally during cell differentiation. For knockdown (KD) experiments, we selected siRNAs that could reduce their target mRNA expression by >30% in C2C12 cells, as determined by reverse transcription-quantitative PCR (RT-qPCR)(Fig 1C; primers listed in S3 Table). For the factors for which we had antibodies, Western Blot analysis was also carried out to confirm success of KD (Fig 1C and S1 Fig; antibody information listed in S2 Table). We extracted total RNA 48 hr after siRNA transfection, or 8 or 24 hr after transfection with U1D or mU1D. To obtain a global view of APA, we applied 3’ Region Extraction And Deep Sequencing (3’READS), a method recently developed by our lab to sequence the 3’ end region of poly(A)+ RNAs genome-wide [10]. The statistics of sequencing reads aligned to pAs are shown in S4 Table.

Regulation of 3’UTR-APA

We first examined APA changes in 3’UTRs, which results in differential expression of isoforms with different 3’UTR lengths (illustrated in Fig 2A). To simplify analysis, we focused on the relative expression changes of the two most abundant 3’UTR isoforms of each gene. Using relative expression difference (RED) between proximal and distal pA isoforms in KD vs. control samples (RED = difference in log2(ratio) of read numbers of two pA isoforms between two samples; illustrated in Fig 2A), we found that several KD samples showed significant upregulation of proximal pAs (negative RED median of all genes), including those of siCFI-68, siCFI-25, siPABPN1, and siPABPC1, and several others had the opposite trend (positive RED median), including the siFip1 and siPcf11 samples (Fig 2B). Differentiation of C2C12 also had a positive RED median, consistent with our previous finding of 3’UTR lengthening in this process [10, 14]. To statistically examine APA regulation, we developed a method named Significance Analysis of Alternative Polyadenylation (SAAP), which evaluated the significance of a RED score using the read distribution of APA isoforms in two comparing samples. A q-value was calculated based on data randomization (illustrated in S2 Fig and see Materials and Methods for detail), which indicated the false discovery rate (FDR). An example gene Timp2 (tissue inhibitor of metalloproteinase 2) is shown in Fig 2C, whose APA was strongly affected by KDs of several factors: siCFI-25, siCFI-68, siPABPN1, and siPABPC1 resulted in upregulation of the proximal pA isoform (3’UTR size = 110 nt), whereas siPcf11 and siFip1 led to upregulation of the distal pA isoform (3’UTR size = 2.6 kb). SAAP analysis of the two isoforms showed that APA regulations in these samples, as well as in C2C12 differentiation, were all significant (q-value = 0, Fig 2C, middle). This result was confirmed by reverse transcription-quantitative PCR (RT-qPCR) using amplicons before and after the proximal pA (Fig 2C, right). Notably, APA of human Timp2 has also been reported in several studies [25, 26], suggesting the importance of its regulation across species.

Fig. 2. 3’UTR-APA.
3’UTR-APA.
(A) Schematic of 3’UTR-APA that results in alternative 3’UTR isoforms. The region between two pAs, pA1 and pA2, is called alternative UTR (aUTR). AAAn, poly(A) tail. CDS and 3’UTR are indicated. Regulation of APA is represented by RED (relative expression difference), whose formula is shown in the figure. #pA1 and #pA2 are numbers of poly(A) site-supporting (PASS) reads for pA1 and pA2, respectively. Test and Ctrl are test and control samples, respectively. (B) RED values for different samples. The top two most abundant APA isoforms based on the number of PASS reads of each gene were analyzed. Only genes with > = 20 PASS reads for proximal and distal pAs combined were used. The median RED value of each sample is shown as a thick black line and interquartile range (between the 25th and 75th percentiles) is indicated by a gray box. (C) An example of 3’UTR-APA regulation. Left, APA isoforms of Timp2. The gene structure and the zoomed-in 3’-most exon are shown on the top. Two pAs are indicated, and their 3’UTR lengths and RPM values in proliferating and differentiated C2C12 cells are shown. Sequence conservation of the shown region (based on mammals) is indicated, with the height of line reflecting the degree of conservation. RT-qPCR amplicons to study APA regulation are indicated. Middle, regulation of the two 3’UTR-APA isoforms of Timp2 in several samples as indicated. RED values (distal pA vs. proximal pA, knockdown (KD) vs. siCtrl or differentiated cells vs. proliferating cells), q-values (Significance Analysis of Alternative Polyadenylation, SAAP, see Materials and Methods for detail), and direction of 3’UTR change (shortened, sh; or lengthened, Le) are indicated. Right, RT-qPCR analysis of 3’UTR-APA isoforms of Timp2 in several samples. RT-qPCR ΔΔCycle threshold (Ct) values were calculated by comparing Ct difference between proximal and distal amplicons in test vs. ctrl samples. Test is a sample from knockdown or differentiated cells, and ctrl is a sample from siCtrl or proliferating cells. Note the samples used for RT-qPCR analysis were not the same as those used for 3’READS. Error bars are standard deviation based on two replicates. (D) Normalized number of genes with regulated 3’UTR-APA in each sample as examined by GAAP. Normalized number is based on (observed value—expected value), and zero is used if a value is negative. Red and blue bars represent genes with lengthened 3’UTRs (Le, distal pA isoform upregulated) and shortened 3’UTRs (Sh, proximal pA isoform upregulated), respectively. Q-value < 0.05 (SAAP) was used to select genes with significant 3’UTR-APA regulation. Only the top two most abundant APA isoforms (based on the number of PASS reads) of each gene were used for this analysis. Error bars are standard deviation based on 20 times of bootstrapping. Samples are sorted by the total number of genes with 3’UTR-APA changes. Data for C2C12 differentiation are shown at the bottom for comparison. Log2(Le/Sh) is log2(ratio) of the number of Le genes to the number of Sh genes. (E) Correlation between normalized no. of genes with 3’UTR-APA (Le+Sh in (D)) and percent of genes with >2 pAs regulated (see S3 Fig). C2C12 differentiation data and siPABPN1 data are shown in red and others in gray. R2 is based on linear regression of all dots. (F) Correlation between log2(Le/Sh) and median aUTR size change. Data from all samples shown in (D) except the siCFI-68 and siCF-25 samples are plotted, with C2C12 differentiation data and siPABPN1 data shown in red diamond and square, respectively, and others in gray dots. R2 is based on linear regression of all gray dots. Median aUTR size change was calculated by the median aUTR size of lengthened 3’UTRs minus that of shortened 3’UTRs. Each gene was calculated once.

We next wanted to compare the extent of APA regulation across samples. We carried out KD experiments in several batches, each with a negative control. However, because samples from different batches were processed and sequenced at different times, systematic biases (“batch effect”) could be introduced into the data. In addition, samples with different levels of sequencing depth could lead to variable sensitivity for APA detection, making them not directly comparable. To address these issues, we developed a computational pipeline, named Global Analysis of Alternative Polyadenylation (GAAP), in which we randomly sampled reads to control for sequencing depth and permutated treatment and control data to obtain expected data (S2C Fig). Both expected and observed data were subject to SAAP for identification of significantly regulated APA events. The number of regulated APA events from the observed data was subtracted by that from the expected data to obtain normalized number of APA events. As such, each sample was internally controlled and cross-batch comparison was more feasible.

Using GAAP, we found that almost all KD samples, including inhibition of U1, resulted in considerable changes of 3’UTR-APA (Fig 2D). The top ten treatments with respect to normalized number of genes with altered 3’UTR-APA were, in the order of number, siCFI-68, siCFI-25, siCPSF-160, siFip1, siSF3b155, siPABPC1, siPABPN1, siPcf11, U1D (24 hr) and siRBBP6. To gauge the overall trend of APA direction, i.e., 3’UTR lengthening or shortening, we calculated the ratio of the number of genes with lengthened 3’UTRs (Le) to the number of genes with shortened 3’UTRs (Sh), or log2(Le/Sh)(Fig 2D). We found that siCFI-25 and siCFI-68 led to the most substantial 3’UTR shortening, with log2(Le/Sh) values being negative infinite after normalization to expected values (Fig 2D). This result is consistent with previous reports showing significant 3’UTR shortening after CFI-25 or CFI-68 KD [24, 43]. KDs of the two PABPs also led to considerable 3’UTR shortening, with log2(Le/Sh) = -2.8 and -2.3, respectively, for siPABPN1 and siPABPC1. While the siPABPN1 result is in line with previous reports [27, 28], this is the first time PABPC1 is found to play a global role in APA and the extent of its modulation of 3’UTR-APA appeared to be similar to that of PABPN1. We found that siFip1 and siPcf11 led to the most substantial 3’UTR lengthening, with log2(Le/Sh) = 1.5 and 1.4, respectively. In comparison, the log2(Le/Sh) value was 1.8 for C2C12 differentiation (Fig 2D). Consistent with previous reports [39], we also found that reduction of U1 activity by U1D (24 hr) or siU1-70K led to 3’UTR shortening. However, the 3’UTR-APA changes in these samples appeared modest compared to the more significant C/P factor KD samples described above. It is worth noting that a similar analysis using the two most significantly regulated pA isoforms of each gene gave essentially identical results (S3A Fig), confirming the robustness of our findings.

As expected, the GAAP result (Fig 2D), which indicates the number of genes with significant APA regulation, correlated with the RED result (Fig 2B), which indicates the extent of APA regulation. We further found that the GAAP result also correlated with the number of significantly regulated pA isoforms per gene (Fig 2E and S3 Fig), i.e., samples with more genes having 3’UTR-APA regulation tended to have a greater number of regulated isoforms per gene. For example, in siCFI-25 and siCFI-68 samples, ~30% of genes had more than two pA isoforms significantly regulated (S3 Fig). This result indicates that multiple pAs are interrelated in regulation of their usage.

The region between two regulated pAs is called alternative UTR (aUTR; illustrated in Fig 2A). We found that when there was global 3’UTR lengthening, the lengthened aUTRs were generally longer than the shortened ones, and vice versa (S3 Fig). Consistently, there was a good correlation between log2(Le/Sh) and median aUTR size difference between lengthened and shortened 3’UTRs (R2 = 0.57, excluding siPABPN1, Fig 2F; note that siCFI-68 and siCFI-25 were not used for this analysis because of their negative infinite log2(Le/Sh) values). The 3’UTR regulation in C2C12 differentiation also followed this trend (Fig 2F). However, the siPABPN1 sample was a notable exception, in which shortened aUTRs were generally shorter than lengthened ones despite a global trend of 3’UTR shortening (Fig 2F and S3 Fig). This result indicates that change of 3’UTR size is generally coupled with the extent of 3’UTR-APA regulation, and the APA regulatory mechanism by PABPN1 is distinct from that of other factors.

Regulation of CDS-APA

We next examined CDS-APA events, as defined by their pA location in introns or exons upstream of the 3’-most exon (illustrated in Fig 3A). Using GAAP, we found that KDs of several splicing factors, such as by siU1-70K and siSF3b155, as well as inhibition of U1 by U1D (8hr and 24 hr), led to substantial changes of the CDS-APA pattern (Fig 3B). By contrast, siU2AF65 resulted in much fewer regulated events. All splicing factor KD samples showed overall upregulation of CDS-APA isoforms, as indicated by their log2(UP/DN) values (UP, number of genes with upregulated CDS-APA isoforms; DN, number of genes with downregulated CDS-APA isoforms)(Fig 3B). Notably, the upregulated CDS-APA isoforms were generally expressed at low levels in control cells (relative abundance <10%, S4 Fig). Taken together, these results indicate that C/P in regions upstream of the 3’-most exon is generally inhibited by splicing and this mechanism is in effect under normal conditions.

Fig. 3. CDS-APA.
CDS-APA.
(A) Schematic of CDS-APA. (B) Normalized number of genes with regulated CDS-APA as examined by GAAP. Red and blue bars represent genes with upregulated CDS-APA isoforms (UP) and downregulated isoforms (DN), respectively. All CDS-APA isoforms were combined and compared to all 3’-most exon isoforms combined by SAAP. Q-value < 0.05 (SAAP) was used to select genes with significant CDS-APA regulation. Error bars are standard deviation based on 20 times of bootstrapping. Samples are sorted by the total number of genes with CDS-APA changes. Data for C2C12 differentiation are shown at the bottom for comparison. Log2(UP/DN) is log2(ratio) of the number of UP genes to the number of DN genes. (C) Normalized expression changes of intronic pA isoforms in several samples. Introns were divided into first (+1), second (+2), last (-1), and second to last (-2), and middle (between +2 and -2 introns) groups. Only genes with ≥4 introns and only pA isoforms with ≥10 PASS reads in two comparing samples combined were analyzed. Expression changes are log2(ratio) of PASS reads in test sample vs. control sample. Values for five intron groups were normalized by mean-centering to reveal bias of intron location. Error bars are standard error of mean. (D) Features of introns containing pAs of upregulated isoforms, including intron size, and 5’ and 3’ splice site (SS) strengths. Numbers are significance score (SS), which was calculated by –log10(P)*S, where P was based on the Wilcoxon rank sum test comparing an intron set of interest with a background set, and S = 1 when the intron set of interest had a larger median value (intron size, 5’SS strength or 3’SS strength) than the background set or -1 otherwise. The background set was derived from introns that contained detectable pA isoform expression in control samples. Introns were divided into five groups based on location, as in (C). The SS data are colored according to the color scheme shown in the graph.

KDs of several C/P factors also led to regulation of many CDS-APA events, including siCFI-25, siCFI-68, siPABPN1, siPABPC1, siFip1 and siPcf11 (Fig 3B). While siCFI-25, siCFI-68, siPABPN1, and siPABPC1 resulted in overall upregulation of CDS-APA isoforms, siFip1 and siPcf11 samples showed the opposite trend. This difference between these two sets of KD samples mirrors that for 3’UTR-APA regulation (Fig 2), implying that regulations of CDS-APA and 3’UTR-APA events by these factors are mechanistically related.

Most pAs of CDS-APA isoforms were located in introns (>90%, S4 Fig). We next asked whether introns containing regulated APA events had special features. Focusing on KD samples with substantial regulations, we analyzed expression changes of the isoforms using pAs in the first (+1), second (+2), second to last (-2), last (-1), or middle (M, not the first two or last two) introns. We found that, in U1D (8 hr), U1D (24 hr), siU1-70K and siPABPN1 samples, APA isoforms using pAs in 5’ end introns tended to be relatively upregulated compared to those in 3’ end introns (Fig 3C). By contrast, isoforms using pAs in the last intron tended to be upregulated in siCFI-68 and siCFI-25 samples (Fig 3C). On the other hand, no obvious location bias could be discerned for upregulated intronic APA isoforms in the siSF3b155, siPABPC1, or siU2AF65 samples (Fig 3C), nor for downregulated intronic APA isoforms in siFip1 or siPcf11 samples (S5 Fig).

We found that the middle introns (not the first two or last two introns) that contained pAs of upregulated isoforms by U1 inhibition or siSF3b155 tended to be much smaller than other middle introns (P = 1x10-8-10-18, Wilcoxon test, Fig 3D). In addition, the middle and last introns containing pAs upregulated by siSF3b155 tended to have stronger 5’SS (P = 1x10-4 and 1x10-5, respectively) than their respective control introns (Fig 3D). Because intron size and 5’SS strength are relevant to splicing kinetics, these results indicate that splicing activity has a substantial influence on C/P in introns. By contrast, 3’SS in general did not appear to play a major role in intronic C/P in the samples analyzed in this study, except for some modest significance for the last introns containing upregulated pAs by siSF3b155 and siU2AF65 (Fig 3D). No significant intron features could be identified for intronic APA isoforms downregulated by siFip1 or siPcf11 (S5 Fig), suggesting that the primary mechanism for their regulation is through modulation of C/P activity.

Regulation of C/P events near the TSS

The enrichment of pAs in the 5’ end introns for upregulated isoforms in siPABPN1 and U1D samples prompted us to examine C/P events near the promoter. Since a large fraction of RNAPII promoters in mammalian cells are bidirectional, leading to both sense and antisense transcripts (reviewed in [44, 45] and illustrated in Fig 4A), we examined C/P events in both sense and antisense orientations around the TSS. Consistent with previous reports [37, 38], we found that in our control samples, the pAs of upstream antisense transcripts, termed uaRNAs or PROMPTs, were distributed within 2 kb from the TSS, peaking around -700 nt (Fig 4B). A similar peak was found in the sense direction with a similar mode of distribution (Fig 4B). For simplicity, we called transcripts using pAs within 2 kb from the TSS in the sense direction (excluding those in 3’-most exons or in single exon genes) sense proximal RNAs (spRNAs).

Fig. 4. C/P events around the transcriptional start site (TSS).
C/P events around the transcriptional start site (TSS).
(A) Schematic of C/P events around the TSS. uaRNA and spRNA are upstream antisense and sense proximal RNAs, respectively. (B) Distribution of uaRNA and spRNA pAs utilized in control C2C12 cells. In this study, we required that the pA of a uaRNA/spRNA was within 2 kb from the TSS, a uaRNA did not overlap with any known protein-coding genes, and the pA of an spRNA was not in the 3’-most exon or in a single-exon gene. (C) Nucleotide frequency profiles around pAs of uaRNAs (left) and spRNAs (middle) and around 3’-most pAs (right). Dotted lines indicate the 25% value. (D) Normalized number of genes with regulated uaRNA expression. Red and blue bars represent genes with upregulated (UP) and downregulated (DN) uaRNA expression, respectively. Log2(UP/DN) is log2(ratio) of the number of UP genes to the number of DN genes. All uaRNAs were combined and compared by SAAP with all sense strand transcripts whose pAs were beyond 2 kb from the TSS. Q-value < 0.05 (SAAP) was used to select genes with a significant uaRNA expression difference. Error bars are standard deviation based on 20 times of bootstrapping. Samples are sorted by the total number of genes with uaRNA expression change. (E) Normalized number of genes with regulated spRNA expression. Red and blue bars represent genes with upregulated (UP) and downregulated (DN) spRNA expression, respectively. All spRNAs were combined and compared to all sense strand transcripts whose pAs were beyond 2 kb from the TSS by SAAP. Q-value < 0.05 (SAAP) was used to select genes with significant spRNA expression difference. Error bars are standard error of mean based on 20 times of bootstrap sampling. Samples are sorted by the total number of genes with spRNA expression changes. Log2(UP/DN) is log2(ratio) of the number of UP genes to the number of DN genes. (F-H) Metagene plots of uaRNA and spRNA expression in siPABPN1 (F), U1D (8 hr and 24 hr, G), and siRrp44 + siRrp6 (H) samples. Expression is represented by reads per million (RPM, poly(A) site-supporting reads only) at pA positions.

Nucleotide frequency analysis of the regions around pAs of uaRNAs and spRNAs indicated that they each had distinct nucleotide profiles compared to those of 3’-most pAs (Fig 4C). While similar U-rich and A-rich upstream peaks and a U-rich downstream peak were found for all types of pAs, the upstream regions of both uaRNA and spRNA pAs had a higher GC content, in line with their location close to promoters with CpG islands [46]. In addition, compared to uaRNA and 3’-most exon pAs, spRNA pAs had a lower A-content in surrounding regions (Fig 4C).

We next examined our KD samples with the goal of finding factors involved in expression of uaRNAs or spRNAs. Using isoforms whose pAs were beyond 2 kb from the TSS as a reference, we examined uaRNA and spRNA regulations in the KD samples by GAAP (Fig 4D and 4E). Significant upregulations of both uaRNA and spRNA were observed with the siPABPN1 sample, with upregulation of uaRNAs being greater than that of spRNAs (Fig 4D and 4E). By contrast, consistent with the findings of Almada et al. [37], regulation by U1 inhibition, such as U1D (8 hr), U1D (24 hr) or siU1-70K, showed the opposite trend, with spRNAs being more significantly upregulated than uaRNAs (Fig 4D and 4E). The siSF3b155 sample also showed significant regulation of uaRNAs, but the number of upregulated events was similar to that of downregulated ones (Fig 4D). In contrast, spRNAs were generally upregulated by siSF3b155 (Fig 4E). The difference between siPABPN1 and U1 inhibition in regulating spRNAs and uaRNAs can also be seen with the isoform expression profiles around the TSS (Fig 4F and 4G).

While upregulation of spRNAs by U1 inhibition may be related to the role of U1 in suppressing C/P near the TSS [37], the regulation of spRNAs and uaRNAs by PABPN1 is completely not clear. Studies have shown that the nuclear exosome is involved in degrading uaRNAs [38] and PABPN1 is involved in stability of nuclear RNAs [47, 48]. We thus asked whether the regulation of uaRNA and spRNA expression by PABPN1 was related to exosome-mediated RNA decay. To this end, we knocked down Rrp44 and Rrp6 together by siRNAs, two nucleases associated with the exosome, followed by 3’READS. As expected, their KD led to higher abundance of uaRNAs and spRNAs, presumably through stabilization of these transcripts (Fig 4H). Importantly, the uaRNA and spRNA profiles in the siRrp44 + siRrp6 sample resembled those of siPABPN1, suggesting a functional connection between PABPN1 and the exosome in regulating uaRNAs and spRNAs.

Detailed study of five C/P factors

Since 3’UTR-APA had much greater regulations than CDS-APA or expression of uaRNAs or spRNAs, as indicated by the number of genes with significant APA changes in KD samples (GAAP analysis results in Figs 2, 3 and 4), we wanted to further examine how different factors regulate 3’UTR-APA. While our GAAP analysis based on random sampling of genes enabled comparison of global APA patterns across samples, this approach is not suitable for comparing individual pA usage across samples due to the batch effect. We therefore repeated KDs of several key C/P factors in one batch, which showed substantial APA regulations in our initial study, including CFI-68, PABPN1, PABPC1, Fip1 and Pcf11. We did not include CFI-25 because the APA profile of its KD was highly similar to that of siCFI-68 (S6 Fig). To mitigate indirect effects and the influence of mRNA decay in cytoplasm, we harvested cells 32 hr after siRNA transfection and extracted both total and nuclear RNAs for 3’READS analysis (Fig 5A). Western Blot analysis indicated that the KD efficiencies in these samples were >60% at the time of cell harvest (Fig 5B).

Fig. 5. Detailed analysis of five C/P factors.
Detailed analysis of five C/P factors.
(A) Schematic of the experimental design. Proliferating C2C12 cells were harvested 32 hr after knockdown (KD), and both total and nuclear RNAs were extracted for 3’READS analysis. (B) Western Blot analysis of protein expression after 32 hr of KD. The percent of expression in KD cells compared to siCtrl cells for each KD is indicated. (C) Normalized number of genes with regulated 3’UTR-APA in each sample as examined by GAAP. See Fig 2D for details of the plot. Both total and nuclear RNA data are shown. (D) Cluster analysis of 3’UTR-APA regulation by the five factors. RED scores using the top two most abundant APA isoforms (based on the number of poly(A) site-supporting reads in all samples) of each gene were used for this analysis. Only pA isoforms with read number ≥5 in all samples were used. A RED score is difference in relative expression of distal pA isoform vs. proximal pA isoform between KD and siCtrl cells, as illustrated in Fig 2A. RED scores are represented in a heatmap using the color scheme shown in the graph. Positive and negative RED scores indicate lengthened and shortened 3’UTRs, respectively. RED scores for APA events were set to 0 when q-value > 0.05 (SAAP). Pearson correlation was used as metric for hierarchical clustering. (E) Venn diagrams comparing genes with significant 3’UTR-APA regulations by the five factors using total RNA (left) or nuclear RNA (right) data. (F) Relationship between extent of 3’UTR-APA regulation and aUTR size. Genes were divided into five bins based on the aUTR size (distance between the pAs of top two most abundant APA isoforms). The aUTR size range for each bin is shown in the graph. The extent of 3’UTR-APA regulation is represented by average RED scores, based on the data shown in (D). Only genes with ≥20 PASS reads (proximal and distal pAs combined) in both KD and siCtrl samples were used for RED calculation. RED scores of genes in bin #1 were compared with those in bin #5 by the Wilcoxon rank sum test for each sample, and p-values are shown.

We found that using total or nuclear RNAs from 32 hr KD samples gave rise to similar results to using total RNAs from 48 hr KD samples (Fig 5C), with siCFI-68, siPABPN1, and siPABPC1 leading to 3’UTR shortening and siPcf11 and siFip1 causing 3’UTR lengthening. This result indicates that our initial observations were not obstructed, at least not substantially, by indirect effects (which could be introduced over time), or by mRNA decay in cytoplasm (which could lead to different APA profiles between nuclear and total RNAs). This notion was further confirmed by detailed comparisons of total RNAs from 32 hr vs. 48 hr KD samples (S7 Fig) and total vs. nuclear RNAs from 32 hr KD samples (S8 Fig), which showed significant correlations of APA changes of genes. Moreover, gene expression changes were largely uncoupled from 3’UTR changes using both nuclear or total RNAs (S9 Fig), indicating that APA isoforms are not likely to differ in mRNA stability in general. However, a mild difference in gene expression could be discerned between genes with shortened and lengthened 3’UTRs in the siPcf11 sample (S9 Fig), suggesting a potential role of Pcf11 in post-transcriptional regulation of gene expression. This will be explored in the future.

Because the five factor KD samples were processed in the same batch, we directly compared their APA profiles. Using RED scores based on the top two most abundant 3’UTR isoforms of each gene, we performed cluster analysis (Fig 5D). We found that total and nuclear RNA data were clustered together for all KD samples, consistent with the notion that cytoplasmic mRNA decay did not substantially alter APA profiles in our samples. In addition, samples involving global 3’UTR shortening, i.e., siCFI-68, siPABPN1 and siPABPC1 samples, were separated from those involving global lengthening, i.e., siPcf11 and siFip1 samples, suggesting that some common sets of pAs were regulated by different factors. On the other hand, Venn diagram analysis of APA events regulated by the five factors also indicated that a fraction of pAs were distinctly regulated by these five factors (Fig 5E): Each factor regulated a set of unique APA events and the number of APA events regulated by all five factors was quite small (27 using total RNA data, and 52 using nuclear RNA data, Fig 5E), indicating distinct regulatory mechanisms among the factors.

Because of the relevance of aUTR length to APA regulation (Fig 2F), we next examined the extent of APA regulation for genes with different aUTR lengths (Fig 5F). Genes were first divided into five equally sized bins based on their aUTR length (Fig 5F). The extent of APA regulation was represented by the RED score. We found that genes that had longer aUTRs tended to have greater APA regulation than genes with shorter ones in all samples except siPABPN1. This is supported by 1) the trend of RED score across gene bins and 2) the p-values indicating the difference in RED scores between the first and fifth gene bins (with the shortest and longest aUTRs, respectively)(Fig 5F). Intriguingly, this analysis also revealed that in the siFip1 sample, genes with the shortest aUTRs (bin 1) showed 3’UTR shortening in general whereas genes in other bins showed the opposite trend (Fig 5F). This dichotomous APA pattern with respect to aUTR size was not seen with other samples, indicating a unique aUTR size-dependent mechanism of Fip1 in 3’UTR-APA regulation (see below for more analyses).

Cis elements in APA regulation by the five C/P factors

We next asked whether cis elements around the pA contributed to differential APA regulation in different KD samples. Since proximal and distal pAs tend to be surrounded by different cis elements [49], we analyzed these two pA groups separately (illustrated in Fig 6A). For each group, we compared pAs that were regulated in one of the five KD samples with pAs regulated in other samples. As such, the identified cis elements should be specific for the factor under investigation, and should not be caused by pA relative locations. We examined three sub-regions around the pA, namely, -100 to -41 nt, -40 to -1 nt and +1 to +100 nt (the pA was set to position 0), for significantly enriched and depleted K-mers (4-mer or 6-mer, P < 0.001, Fisher’s exact test). Fig 6B shows the number of significant 4-mers for each region, reflecting the extent to which cis elements in the region were involved in pA regulation.

Fig. 6. Cis elements associated with regulated pAs.
Cis elements associated with regulated pAs.
(A) Schematic showing the analysis method. As indicated, pAs were divided into proximal and distal pA groups, and pAs of regulated isoforms were compared only with other pAs in the same group. As such, proximal pAs were only compared with proximal pAs and so were distal pAs. pAs of upregulated and downregulated isoforms were analyzed separately. Only data of nuclear RNA samples were used for analysis, because they were expected to have less post-transcriptional effects than total RNA samples. (B) Number of 4-mers with significantly biased frequency of occurrence (P < 0.001, Fisher’s exact test) near regulated pAs. Regulated pAs were those with q-value < 0.05 (SAAP). Three regions around the pA were analyzed, including -100 to -41 nt, -40 to -1 nt and +1 to +100 nt. Data for all 4-mers and top 6-mers are shown in S5 Table and S6 Table, respectively. (C) Significant 4-mers enriched for or depleted from pAs regulated by siCFI-68. Only top five 4-mers for the regions with ≥5 significant 4-mers are shown. Numbers are significance score (SS), which was calculated by –log10(P)*S, where P was based on the Fisher’s exact test and S = 1 for enrichment and -1 for depletion. (D) As in (C), significant 4-mers enriched for or depleted from pAs regulated by siFip1. (E) Regulation of different types of pAs by siFip1, as shown by Cumulative Distribution Function (CDF) curves of RED scores. Genes were divided into four groups based on i) distance between proximal and distal pAs (<120 nt or ≥120 nt), and 2) whether or not there was AAUAAA within 100 nt downstream of the proximal pA. These groups are also illustrated in the graph. The number of genes and the median RED score of each group are shown in a table. The differences between groups are indicated by p-values (Kolmogorov–Smirnov test).

We found much greater 4-mer biases around pAs regulated by siCFI-68 or siFip1 than those by siPcf11, siPABPC1 or siPABPN1 (Fig 6B and S5 Table). In the siCFI-68 sample, major biases were found in regions surrounding both proximal and distal pAs of genes with shortened 3’UTRs (Fig 6B and 6C), including depletion of TGTA in the -100 to -41 nt region of proximal pAs (significance score (SS) = -11.1; SS = -log10(p-value)*S, where p-value was based on the Fisher’s exact test and S was 1 for enrichment or -1 for depletion) and enrichment of this motif in the same region of distal pAs (SS = 21.0), suggesting that the presence or absence of TGTA in the -100 to -41 nt region is a major reason for APA regulation by CFI-68. This result is consistent with the binding site for CFI complex [24, 43, 50]. Other motifs were also found to be significantly biased in these regions, as well as in the +1 to +100 nt region of proximal pAs, and the -40 to -1 nt and +1 to +100 nt regions of distal pAs, albeit with less significance (Fig 6C).

In the siFip1 sample, both pAs of shortened 3’UTRs and of lengthened 3’UTRs displayed motif biases (Fig 6D). For lengthened 3’UTRs, there was an enrichment of TTTT in the -100 to -41 nt region of proximal pAs, and an depletion of the motif in the same region of distal pAs, suggesting that U-rich elements play a role in APA regulation mediated by Fip1. This result is consistent with the reported U-rich sequence binding for Fip1 [51] and is in good agreement with the binding locations reported by Lackford et al. [31]. Several 4-mers displayed a strong bias in the +1 to +100 nt region of proximal pAs of shortened 3’UTRs, including TAAA, AATA and ATAA. Hexamer analysis indicated that these motifs were derived from AATAAA (S6 Table). This result suggests that Fip1 inhibits usage of pAs with downstream AATAAA. To test this hypothesis further, we specifically examined pAs with or without downstream AATAAA and analyzed the influence of distance between pAs (Fig 6E). We found that, in the siFip1 sample, when the distance between proximal and distal pAs was < 120 nt, proximal pAs with downstream AATAAA within 100 nt (group 1, Fig 6E) tended to be much more used than those without the motif (group 2), as indicated by their RED scores (-0.44 vs. -0.03, P < 1x10-6, Kolmogorov–Smirnov, or K-S, test,). This trend could also been seen when the distance was > = 120 nt (0.11 vs. 0.30, group 3 vs. group 4, P = 0.14). Thus, Fip1 plays a role in selection of adjacent pAs, favoring AATAAA-associated downstream pAs. This result offers an explanation as to why genes with short aUTRs tended to have upregulated proximal pAs in the siFip1 sample (Fig 5F).

Some cis element biases were also found in siPcf11 and siPABPN1 samples, with lower statistical significance (Fig 6B, S5 Table, and S6 Table). For example, TATT and TTAT were enriched for the +1 to +100 nt region of proximal pAs downregulated by siPcf11, and several TA-rich motifs were depleted from the -100 to -41 nt region of downregulated distal pAs but enriched for the same region of upregulated distal pAs in the siPABPN1 sample. Of all samples, siPABPC1 displayed the least cis element bias around regulated pAs, suggesting that its regulation of C/P does not involve specific cis elements.

Relevance of the APA regulation by the five C/P factors to cell differentiation and development

We next asked whether the regulated APA events by five C/P factors were related to those taking place during C2C12 differentiation [14]. We examined 3’UTR-APA in C2C12 differentiation for genes that showed shortened or lengthened 3’UTRs in different KD samples (32 hr KD, total RNA). Using RED scores for differentiated vs. proliferating C2C12 samples, we found that genes regulated by siCFI-68 and siPcf11 showed significant RED differences in C2C12 differentiation (Fig 7A): genes with shortened 3’UTRs in the siCFI-68 sample and those with lengthened 3’UTRs in the siPcf11 sample were more likely to have lengthened 3’UTRs in C2C12 differentiation than genes with the opposite 3’UTR regulation in their respective samples (P = 4x10-3 and 7x10-5, respectively, K-S test).

Fig. 7. Comparison of APA events regulated by the five C/P factors with those regulated in C2C12 differentiation.
Comparison of APA events regulated by the five C/P factors with those regulated in C2C12 differentiation.
(A) APA regulation in C2C12 differentiation for genes that showed shortened 3’UTRs (Sh, blue line) or lengthened 3’UTRs (Le, red line) in different knockdown (KD) samples (total RNA only). APA regulation in C2C12 differentiation is represented by Cumulative Distribution Function (CDF) curves of RED scores. Only genes with q-value < 0.1 (SAAP) were used. The numbers of Sh genes and Le genes are shown in blue or red, respectively. The Δ value is difference between the median RED scores of Le and Sh genes. Median RED scores are also indicated by vertical dotted lines. P-values (Kolmogorov–Smirnov test) are shown to indicate RED score difference between Le and Sh genes in C2C12 differentiation. (B) Ten groups of genes based on 3’UTR-APA regulation by siCFI-68 and siPcf11 (q-value < 0.1). Group 9 contained genes whose 3’UTRs were regulated by siPABPN1, siPABPC1, or siFip1. Group 10 contains other genes whose 3’UTR-APA isoforms were detectable in C2C12 cells. Le, 3’UTR lengthened; Sh, 3’UTR shortened. Number of genes and median aUTR size are also shown. (C) Median RED scores for gene groups shown in (B) in C2C12 differentiation. Group 6 is not included because of a very small number of genes (5) in the group. (D) RED scores of group 3 genes were compared with those of genes in other groups in C2C12 differentiation (left), 3T3-L1 differentiation (middle), and embryonic development (15 day vs. 11 day, right). As in (A), difference between the median RED scores of two sets (Δ) and p-value (Kolmogorov–Smirnov test) comparing the two sets are also shown.

Based on APA regulations by siCFI-68 and siPcf11, we next divided genes into 10 groups (Fig 7B and S7 Table) and asked how different sets of genes were regulated in C2C12 differentiation. Genes in group 3, whose 3’UTRs were shortened by siCFI-68 and lengthened by siPcf11 showed greatest 3’UTR lengthening in C2C12 differentiation as indicated by their RED median (Fig 7C). They were followed by group 2 genes, whose 3’UTRs were shortened by siCFI-68 but not regulated by siPcf11, and group 5 genes, whose 3’UTRs were lengthened by siPcf11 but not by regulated by siCFI-68 (Fig 7C). Interestingly, genes in groups 3, 2 and 5 had longer aUTRs than genes in other groups, with median size of 1,618, 1,058 and 866 nt, respectively (Fig 7B), highlighting the importance of distance between pAs for APA regulation in C2C12 differentiation. Further analysis using our previous 3’READS data from 3T3-L1 pre-adipocyte differentiation and embryonic development (15 day vs. 11 day)[10] indicated that group 3 genes also tended to have significantly greater 3’UTR lengthening than other genes in these processes (P = 6x10-8 and 1x10-2, K-S test comparing group 3 genes with all other genes, Fig 7D), similar to their regulation in C2C12 differentiation (P = 4x10-7, Fig 7D). Interestingly, group 3 genes were significantly associated with several Gene Ontology (GO) terms (S10 Fig), such as “single-organism membrane organization”, “spermatogenesis”, “actomyosin structure organization”, “intracellular protein transport”, “cell body” and “cytoplasmic vesicle.” Notably, several different GOs were also found to be associated with group 2 genes, such as “endosome to lysosome transport”, “transforming growth factor beta receptor signaling pathway”, “cellular response to endogenous stimulus”, “cell leading edge”, “coated vesicle membrane”, “cell junction”, and “endomembrane system.” Intriguingly, group 1 genes, whose 3’UTRs were shortened by both siCFI-68 and siPcf11, also showed strong association with several GO terms, such as “regulation of phosphorylation”, “cell death”, “intracellular signal transduction”, etc. Taken together, these data indicate that genes in different functional groups may be differentially regulated by APA in cell differentiation due to distinct pA and aUTR features.

Discussion

Here we present a systematic study examining modulation of several types of APA by all known core C/P factors and several key splicing factors. Using proliferating C2C12 cells, we were able to compare the impacts of different factors on APA in a single system, and correlate the regulations with those taking place during C2C12 differentiation. We also developed a statistical method SAAP to assess the significance of APA regulation and employed a computational approach GAAP to address batch effects and normalize sequencing depth. Comparisons between KD samples at different time points and between total and nuclear RNA samples indicated that our results were not substantially affected by indirect effects. However, off-target effects by siRNAs are possible, despite the use of multiple siRNAs in most samples. In addition, different KD efficiencies by different siRNAs can lead to variable effects on APA, making it difficult to compare KD samples directly. However, we did not find a correlation between the level of KD and the extent of 3’UTR-APA regulation (S11 Fig), suggesting our observations were not severely affected by variable KD efficiencies. Nevertheless, future work using more specific KD or knockout methods may be desirable for more precise comparisons.

Notwithstanding the potential technical limitations, we were able to make a number of key findings concerning the contributions of specific C/P and splicing factors to APA, and to address how they pertain to APA regulation in cell proliferation/differentiation. Below we discuss these findings; some of the conclusions are illustrated in Fig 8.

Fig. 8. APA models.
APA models.
(A) Regulation of C/P events near the transcription start site (TSS) and introns. Both U1 and PABPN1 inhibit polyA(+) transcript expression near the TSS. For U1, inhibition of sense strand RNA is more prominent than that of upstream antisense RNA (uaRNA) due to higher frequency of its binding sites, i.e., 5’ splice site or related sequences, in the sense strand. The mechanism is likely to be inhibition of polyadenylation. PABPN1 has the opposite trend, with suppression of uaRNA expression being more significant. This function of PABPN1 is likely to involve exosome-mediated RNA decay. U1 and U2 inhibit intronic C/P because C/P is in a kinetic competition with splicing. (B) Regulation of pA usage in the last intron and the 3’-most exon. CFI-25/68 promotes usage of distal pAs through binding to the UGUA element, leading to longer 3’UTRs and selection of downstream terminal exons. PABPN1 and PABPC1 also promote distal pA usage, whereas Pcf11 and Fip1 promote proximal pA usage. Note that if two pAs are close to each other and there is an AAUAAA motif downstream of the proximal pA, Fip1 helps select the distal pA (not indicated in the graph). aUTR size is indicated to highlight its importance in APA regulation.

CPSF

Among the CPSF subunits, Fip1 levels were found to have the greatest impact on APA, suggesting that Fip1 recruitment is a key step for CPSF’s function in C/P, at least in the context of C2C12 cells. Given the essential role of CPSF in the C/P reaction, it is not surprising that inhibition of Fip1 leads to 3’UTR lengthening: downregulated C/P activity would favor distal pAs, which in general have stronger C/P cis elements than proximal ones [6]. Interestingly, a group of genes had 3’UTR shortening after Fip1 KD; their proximal pAs tend to have the AAUAAA motif in the downstream region and the distance between proximal and distal pAs is short. This result is largely in line with what was observed by the Shi group with Fip1 KD in embryonic stem cells [31]. As proposed by Lackford et al., there exist two modes of APA regulation by Fip1 depending upon pA strength and distance between pAs [31]. We additionally found that while selection of distal pAs in the siFip1 sample involves U-rich elements, activation of proximal pAs does not, suggesting that U-rich element binding by Fip1 plays a role in pA selection only when competing pAs are far apart. Recent studies have revealed that CPSF30 and WDR300 directly bind PAS, [52, 53]. Intriguingly, KD of any one of these factors had a modest effect on APA. Whether these factors can compensate each other in PAS interaction is to be examined in the future.

CFI

Consistent with the findings by the Zavolan and Keller groups [24, 43], KD of CFI-25 or CFI-68, but not CFI-59, led to a substantial shift to proximal pA usage in 3’UTRs, the extent of which is more significant than other samples examined in this study. The CFI complex, composed of either CFI-25/CFI-68 or CFI-25/CFI-59, has been shown to interact with UGUA elements [50], which are typically located upstream of the PAS [49]. Consistently, distal pAs downregulated by siCFI-25/siCFI-68 were enriched with UGUA element(s) in the -100 to -41 nt region whereas upregulated proximal pAs were depleted of them in the same region. Some other cis elements were also identified around regulated pAs (Fig 6C), such as UAUU, CCUC. Whether they simply piggyback with UGUA or are actively engaged in binding with some proteins that interact with CFI need to be further studied in the future.

We also found that CFI-25 or CFI-68 KD led to general upregulation of isoforms using pAs in the last intron, a feature that has not been reported in previous studies [24, 43], suggesting that CFI-25/68 also play a role in 3’ terminal exon selection. One possibility is that removal of the last intron is slow, creating a time window in which pAs in the last intron compete with 3’-most exon pAs for C/P. It is also possible that CFI-25/CFI-68 may facilitate the removal of the last intron, e.g., through interaction with splicing factors [33], thereby inhibiting pA usage in the last intron.

CFII

Our data indicate that Pcf11 has a substantial impact on APA, which has not been detected previously in metazoan cells. In S. cerevisiae, deletion of the mRNA export adaptor Yra1, which inhibits C/P through competition with Clp1 for Pcf11 binding, leads to widespread APA events [54]. The same study also suggested that Clp1 and Pcf11 are not necessarily recruited as a complete unit in yeast. Our data supports this view, because siPcf11 and siClp1 had different effects on APA, with respect to extent and direction of regulation. This notion is also in line with the report by Shi et al. showing that Clp1 was not present in the C/P complex purified from mammalian cells [22]. Given the significant role of Pcf11 in APA, it is possible that its recruitment to the C/P complex is a rate-limiting step for the C/P reaction.

PABPs

Inhibition of PABPN1 has been shown to elicit global shortening of 3’UTRs via APA [27, 28]. It was suggested that this regulation may be relevant to the etiology of human disease oculopharyngeal muscular dystrophy (OPMD), which is caused by an expansion mutation in the polyalanine repeat in the N-terminus of PABPN1. Our results are largely in line with these findings but also revealed that modulation of 3’UTR-APA by PABPN1 is quite different than other C/P factors. For example, there is no relationship between the extent of 3’UTR-APA regulation (based on number of genes with APA changes) and 3’UTR size changes. We also found that PABPN1 plays a very significant role in inhibition of uaRNA and spRNA expression (with regulation of uaRNAs being more significant). This property may be related to PABPN1’s role in RNA stability [47, 48]. Consistently, the uaRNA and spRNA profiles in the siPABPN1 sample are similar to those in the siRrp44 + siRrp6 sample. A key question thus is whether the 3’UTR-APA regulation by PABPN1 is related to its functions around the TSS. We did not find significant 3’UTR-APA in the siRrp44/siRrp6 sample (S12 Fig) and genes with 3’UTR-APA regulation by siPABPN1 did not appear to have significant expression changes overall (S9 Fig), suggesting that 3’UTR-APA regulation by PABPN1 is not related to its function in uaRNA and spRNA regulation. How PABPN1 exerts distinct functions at the 5’ and 3’ ends of genes respectively awaits further experimentation.

Interestingly, PABPC1, which binds the poly(A) tail in the nucleus [55] and shuttles between nucleus and cytoplasm [56], appears to be as potent as PABPN1 in regulation of 3’UTR-APA, raising the possibility that PABPs in general can regulate the C/P reaction. However, it is also notable that the APA events regulated by PABPC1 are quite different than those by PABPN1, with respect to uaRNA and spRNA expression changes and the role of aUTR size in regulation. There are no obvious cis elements enriched for pAs regulated by PABPC1 but aUTR size appears to be related to its APA regulation, suggesting that PABPC1 may have a general impact on APA without sequence specificity. Future work will be needed to delineate the respective roles of these PABPs as well as other poly(A)-binding proteins [57] in APA.

CstF and other C/P factors

CstF subunits appear to have modest impacts on APA in this study. While siCstF-50 and siCstF-77 led to global 3’UTR lengthening, consistent with their roles in C/P, siCstF-64 did not result in a directional change of APA and siCstF-64τ in fact led to mild general 3’UTR shortening. These puzzling results may be due to the fact that CstF-64 and CstF-64τ have overlapping functions [29] and can regulate each other’s expression [58]. As such, KD of one factor to different levels may give rise to different overall activities, leading to complex results. Thus, the single factor KD-based approach used in this study would be unwieldy in elucidating the exact roles of CstF-64 and CstF-64τ in APA. This issue may also apply to other factors which showed modest APA regulation in this study. For example, there are several PAPs in the cell and KD of one PAP, such as PAPα or PAPγ, may be compensated functionally by other PAPs.

Splicing in APA

Splicing is believed to be intimately involved in C/P (reviewed in [33]). Our data define two general modes of APA regulation involving splicing. First, U1 plays a significant role in suppressing C/P in 5’ end introns, which is consistent with the findings made by the Dreyfuss group [35] and is in line with the inhibitory activity of U1 on C/P [36]. In agreement with the telescripting model proposed by the Dreyfuss group [39], inhibition of intronic pA usage by U1 displays 5’ to 3’ polarity (Fig 3C). But we only observed mild 3’UTR lengthening in U1 inhibition samples, suggesting a modest role of telescripting in 3’UTR-APA. However, we cannot rule out the possibility that difference in cell type and experimental conditions may lead to this discrepancy. Second, as supported by the siSF3b155 data, splicing activity in general plays a role in suppressing intronic C/P (Fig 3B). It remains to be seen how other splicing factors, particularly those involved in alternative splicing, globally regulate intronic APA and impact selection of 3’ terminal exons.

APA code in proliferation/differentiation

Our data suggest that an APA event in a given cellular condition is regulated by a number of parameters, including relative location to the TSS, splicing context, distance between competing pAs, surrounding cis elements, and concentrations of C/P factors. In the context of C2C12 differentiation, which involves global 3’UTR lengthening, almost all C/P factors showed downregulation, at least at the mRNA level (S13 Fig). However, given the diverse consequences of different C/P factor KDs, it would be too simplistic to attribute the 3’UTR lengthening to downregulation of C/P factors as a whole. On the other hand, several lines of evidence support similarities in 3’UTR-APA regulation between C/P factor KD and C2C12 differentiation. First, aUTR size appears to be an important factor in APA regulation in both KD cells (with the exception of siPABPN1) and cell differentiation. In general, a longer aUTR confers more regulability. Since longer aUTRs could contain more cis regulatory elements for mRNA metabolism, this result suggests that genes with highly regulatable APA tend to be more controlled post-transcriptionally as well. Second, groups of genes with significant regulation by siCFI-68 and siPcf11 are also substantially regulated in differentiation, implying similar mechanisms in these KD conditions and cell differentiation. Importantly, these genes also displayed significant 3’UTR-APA in differentiation of 3T3-L1 cells and embryonic development, suggesting a general APA code in cell proliferation/differentiation. Future studies need to test the APA code with more perturbations, such as overexpression of different factors, and to explore the input of different RNA-binding proteins in condition- and tissue-specific APA regulations [59].

Materials and Methods

Cell culture, transfection, RT-qPCR and western blot

C2C12 cells were maintained in Dulbecco's Modified Eagles Medium (DMEM) supplemented with 10% fetal bovine serum (FBS). Differentiation of C2C12 cells was induced by switching cell media to DMEM+ 2% horse serum (Sigma) when cells were ~90–100% confluent. All media were also supplemented with 100 units/ml penicillin and 100 μg/ml streptomycin. Differentiated C2C12 cells in this study were harvested four days after differentiation. Transfection with siRNAs was carried out by Lipofectamine 2000 (Invitrogen) according to manufacturer’s recommendations. Transfection was carried out for 48 hr or 32 hr. siRNA sequences are shown in S1 Table. The U1D oligo (5'-gCcAgGuAaGuau) and mutant U1D (mU1D) oligo (5'-gCcAgGcAcGuau), where locked nucleic acid (LNA) residues are in uppercase and 2’-OMe RNA bases in lowercase, were previously described in [42]. These oligos were transfected into C2C12 cells at 35 nM using Lipofectamine 2000 when the confluency of cells was about 50%. Cells were harvested 8 hr or 24 hr after transfection. For nuclear RNA extraction, cells collected by a scrapper were suspended in cell lysis buffer (10 mM Tris pH 7.4, 10 mM NaCl, 0.5% NP40, 1 mM DTT), followed by vortexing for 10 sec and incubation on ice for 10 min. After centrifugation of the lysate at 500 x g for 5 min at 4°C, the pellet was re-suspended in the cell lysis buffer for nuclear RNA extraction. Both total and nuclear RNAs were extracted using Trizol (Invitrogen) according to manufacturer's protocol. RNA quality was examined in an Agilent Bioanalyzer using the RNA pico600 kit before processing for deep sequencing. For RT-qPCR, mRNA was reverse-transcribed using an oligo(dT) primer, and qPCR was carried out with Syber-Green I as dye. Primer sequences are listed in S3 Table. Antibodies used for Western Blot analysis are listed in S2 Table.

3’READS

The 3’READS method used in this study was previously described [60]. Briefly, 25 μg of input RNA was used for each sample, and poly(A)+ RNA was selected using oligo d(T)25 magnetic beads (NEB), followed by on-bead fragmentation using RNase III (NEB). Poly(A)+ RNA fragments were then selected using the chimeric U5 and T45 (CU5T45) oligo conjugated on streptavidin beads, followed by RNase H (NEB) digestion. Eluted RNA fragments were ligated with 5’ and 3’ adapters, followed by RT and PCR (15x) to obtain cDNA libraries for sequencing on the Illumina platform. All data can be obtained from the NCBI GEO database (GSE62001). Processing of 3’READS data was carried out as previously described (Hoque et al. 2013). Briefly, reads were mapped to the mouse genome using bowtie 2 (Langmead and Salzberg 2012). Reads with ≥2 unaligned Ts at the 5’ end are called poly(A) site-supporting (PASS) reads, which were used to identify pAs. pAs located within 24 nt from each other were clustered together. The number of PASS reads generated in each sample is listed in S4 Table.

Significance analysis of alternative polyadenylation

We developed a randomization-based method to statistically assess the significance of difference between two samples for each APA event, called Significance Analysis of Alternative Polyadenylation, or SAAP. The method is illustrated in S2 Fig. Briefly, for two pAs (or two pA sets) from two comparing samples, a Relative Expression Difference (RED) score is first calculated (S2 Fig). The PASS reads are then sampled based on the assumption that the relative abundance of each pA isoform is the same in two samples. Sampling is preformed m times (20 in this study) to obtain mean and standard deviation of RED. The observed and expected RED values are standardized to Z scores (minus mean and divided by standard deviation). False Discovery Rate (FDR) and q-value are calculated by comparing observed Z (Zo) and expected Z (Ze) for a given Z cutoff value (Zc).

For 3’UTR-APA, we either selected the two most abundant pA isoforms for analysis, or used all pA isoforms for analysis (one vs. others) and then selected the top two most regulated isoforms. For CDS-APA, we combined all isoforms using pAs in upstream regions of the 3’-most exon and compared their expression change with that of isoforms using pAs in the 3’-most exon. Individual intronic pAs were also analyzed by comparing to all other pA isoforms of the gene. For uaRNA analysis, we combined all antisense transcripts using pAs within 2 kb upstream of the transcription start site (TSS), excluding those mapped to other mRNA genes, and compared them to all sense strand transcripts, excluding pAs located within 2 kb downstream of the TSS. For spRNA analysis, we grouped pAs located within 2 kb downstream of the TSS, excluding those located in 3’-most exons or in single exon genes, and compared them with other sense strand isoforms. We used q-value < 0.05 (SAAP) to select significantly regulated APA events.

Global analysis of alternative polyadenylation

We developed an approach named Global Analysis of Alternative Polyadenylation (GAAP) to normalize sequencing depth and to obtain expected values for a given data set. The method is illustrated in S2 Fig. Let A and B be two 3’READS data sets for two samples that are processed at the same time. For example, A is a KD sample and B is a control sample. Let a and b be two data sets sampled by bootstrapping from A and B, respectively. Bootstrapping is carried out n times with p number of PASS reads sampled each time. In this study, n = 20, and p = 1.5 M. Let A’ and B’ be randomly permutated data sets based on A and B. In permutation, PASS reads from A and B are first combined and then sampled without replacement to obtain A’ and B’, with the total number of PASS reads of the permutated sets A’ and B’ being the same as those of the original sets A and B, respectively. Permutation is carried out m times. Two sets after each permutation (a’ and b’) are sampled with q number of reads from each set by bootstrapping. In this study, m = n, and q = p. The permutated sets provide expected values (Exp), whereas the original set provides observed values (Obs). APA analyses of a vs. b and a’ vs. b’ were carried out by SAAP.

K-mer analysis

To identify cis elements biased to or against pAs regulated in a KD sample, we compared k-mer frequencies around the pAs regulated in the sample with other pAs regulated in other samples. To mitigate the possibility that identified cis elements are related to location, we first grouped all proximal and distal pAs together into two separate sets (illustrated in Fig 6A). For each KD sample, regulated proximal pAs were compared to other proximal pAs to identify overrepresented and underrepresented k-mers. The same approach was used for distal pAs. We examined three regions around the pA, i.e., -100 to -41 nt, -40 to -1 nt and +1 to +100nt. For each region, the Fisher’s exact test was used to examine whether a k-mer (4-mer or 6-mer) was enriched for or depleted from a set of pAs vs. other pAs.

Analysis of introns

The intron location was based on the RefSeq database, with all RefSeq splicing isoforms combined. Distribution of regulated intronic pA isoforms was compared to that of background set, which was derived from all detected intronic pAs in all control samples in this study (isoform relative abundance ≥5% in at least two samples and read count ≥2 in at least two samples). To calculate 5’SS or 3’SS strength, we used all 5’SS or 3’SS supported by mouse RefSeq sequences. The maximum entropy scores were calculated by MaxEntScan [61]. The 5’SS or 3’SS strength of introns containing regulated pAs was compared to that of background introns with the same relative location in the gene by the Wilcoxon rank sum test.

Other analyses

Gene expression changes were based on PASS reads mapped to the 3’-most exon of a gene, represented by the reads per million total PASS reads (RPM) value. Venn diagrams were generated by VennDIS [62]. Gene Ontology (GO) analysis was carried out using the Fisher’s exact test. GO annotation of genes was obtained from the NCBI Gene database.

Supporting Information

Attachment 1

Attachment 2

Attachment 3

Attachment 4

Attachment 5

Attachment 6

Attachment 7

Attachment 8

Attachment 9

Attachment 10

Attachment 11

Attachment 12

Attachment 13

Attachment 14

Attachment 15

Attachment 16

Attachment 17

Attachment 18

Attachment 19

Attachment 20


Zdroje

1. 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(2):405–45. 10357856

2. Colgan DF, Manley JL. Mechanism and regulation of mRNA polyadenylation. Genes Dev. 1997;11(21):2755–66. 9353246

3. Tian B, Graber JH. Signals for pre-mRNA cleavage and polyadenylation. WIREs RNA. 2012, 3: 385–396. doi: 10.1002/wrna.116 22012871

4. Proudfoot NJ. Ending the message: poly(A) signals then and now. Genes Dev. 2011;25(17):1770–82. doi: 10.1101/gad.17268411 21896654

5. Kadonaga JT. Perspectives on the RNA polymerase II core promoter. Wiley interdisciplinary reviews Developmental biology. 2012;1(1):40–51. doi: 10.1002/wdev.21 23801666

6. 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(1):201–12. 15647503

7. Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, et al. Alternative isoform regulation in human tissue transcriptomes. Nature. 2008;456(7221):470–6. doi: 10.1038/nature07509 18978772

8. Shepard PJ, Choi EA, Lu J, Flanagan LA, Hertel KJ, Shi Y. Complex and dynamic landscape of RNA polyadenylation revealed by PAS-Seq. Rna. 2011;17(4):761–72. doi: 10.1261/rna.2581711 21343387

9. Derti A, Garrett-Engele P, Macisaac KD, Stevens RC, Sriram S, Chen R, et al. A quantitative atlas of polyadenylation in five mammals. Genome Res. 2012;22(6):1173–83. doi: 10.1101/gr.132563.111 22454233

10. Hoque M, Ji Z, Zheng D, Luo W, Li W, You B, et al. Analysis of alternative cleavage and polyadenylation by 3' region extraction and deep sequencing. Nat Methods. 2013;10(2):133–9. doi: 10.1038/nmeth.2288 23241633

11. Zhang H, Lee JY, Tian B. Biased alternative polyadenylation in human tissues. Genome Biol. 2005;6:R100. 16356263

12. Lianoglou S, Garg V, Yang JL, Leslie CS, Mayr C. Ubiquitously transcribed genes use alternative polyadenylation to achieve tissue-specific expression. Genes Dev. 2013;27(21):2380–96. doi: 10.1101/gad.229328.113 24145798

13. Sandberg R, Neilson JR, Sarma A, Sharp PA, Burge CB. Proliferating cells express mRNAs with shortened 3' untranslated regions and fewer microRNA target sites. Science. 2008;320(5883):1643–7. doi: 10.1126/science.1155390 18566288

14. Ji Z, Lee JY, Pan Z, Jiang B, Tian B. Progressive lengthening of 3' untranslated regions of mRNAs by alternative polyadenylation during mouse embryonic development. Proc Natl Acad Sci U S A. 2009;106:7028–33. doi: 10.1073/pnas.0900028106 19372383

15. Ji Z, Tian B. Reprogramming of 3' untranslated regions of mRNAs by alternative polyadenylation in generation of pluripotent stem cells from different cell types. PLoS One. 2009;4(12):e8419. doi: 10.1371/journal.pone.0008419 20037631

16. Mayr C, Bartel DP. Widespread shortening of 3'UTRs by alternative cleavage and polyadenylation activates oncogenes in cancer cells. Cell. 2009;138(4):673–84. doi: 10.1016/j.cell.2009.06.016 19703394

17. Flavell SW, Kim TK, Gray JM, Harmin DA, Hemberg M, Hong EJ, et al. Genome-wide analysis of MEF2 transcriptional program reveals synaptic target genes and neuronal activity-dependent polyadenylation site selection. Neuron. 2008;60(6):1022–38. doi: 10.1016/j.neuron.2008.11.029 19109909

18. Tian B, Manley JL. Alternative cleavage and polyadenylation: the long and short of it. Trends Biochem Sci. 2013;38(6):312–20. doi: 10.1016/j.tibs.2013.03.005 23632313

19. Elkon R, Ugalde AP, Agami R. Alternative cleavage and polyadenylation: extent, regulation and function. Nat Rev Genet. 2013;14(7):496–506. doi: 10.1038/nrg3482 23774734

20. Shi Y. Alternative polyadenylation: new insights from global analyses. Rna. 2012;18(12):2105–17. doi: 10.1261/rna.035899.112 23097429

21. Mandel CR, Bai Y, Tong L. Protein factors in pre-mRNA 3'-end processing. Cell Mol Life Sci. 2008;65(7–8):1099–122. doi: 10.1007/s00018-008-8530-3 18989620

22. Shi Y, Di Giammartino DC, Taylor D, Sarkeshik A, Rice WJ, Yates JR 3rd, et al. Molecular architecture of the human pre-mRNA 3' processing complex. Mol Cell. 2009;33(3):365–76. doi: 10.1016/j.molcel.2008.12.028 19217410

23. Takagaki Y, Seipelt RL, Peterson ML, Manley JL. The polyadenylation factor CstF-64 regulates alternative processing of IgM heavy chain pre-mRNA during B cell differentiation. Cell. 1996;87(5):941–52. 8945520

24. Martin G, Gruber AR, Keller W, Zavolan M. 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. Cell Rep. 2012;1(6):753–63. doi: 10.1016/j.celrep.2012.05.003 22813749

25. Kubo T, Wada T, Yamaguchi Y, Shimizu A, Handa H. Knock-down of 25 kDa subunit of cleavage factor Im in Hela cells alters alternative polyadenylation within 3'-UTRs. Nucleic Acids Res. 2006;34(21):6264–71. 17098938

26. Masamha CP, Xia Z, Yang J, Albrecht TR, Li M, Shyu AB, et al. CFIm25 links alternative polyadenylation to glioblastoma tumour suppression. Nature. 2014;510(7505):412–6. doi: 10.1038/nature13261 24814343

27. Jenal M, Elkon R, Loayza-Puch F, van Haaften G, Kuhn U, Menzies FM, et al. The poly(A)-binding protein nuclear 1 suppresses alternative cleavage and polyadenylation sites. Cell. 2012;149(3):538–53. doi: 10.1016/j.cell.2012.03.022 22502866

28. de Klerk E, Venema A, Anvar SY, Goeman JJ, Hu O, Trollet C, et al. Poly(A) binding protein nuclear 1 levels affect alternative polyadenylation. Nucleic Acids Res. 2012;40(18):9089–101. doi: 10.1093/nar/gks655 22772983

29. Yao C, Choi EA, Weng L, Xie X, Wan J, Xing Y, et al. Overlapping and distinct functions of CstF64 and CstF64tau in mammalian mRNA 3' processing. Rna. 2013. 19(12):1781–90. doi: 10.1261/rna.042317.113 24149845

30. Li W, Yeh HJ, Shankarling GS, Ji Z, Tian B, Macdonald CC. The tauCstF-64 Polyadenylation Protein Controls Genome Expression in Testis. PLoS One. 2012;7(10):e48373. doi: 10.1371/journal.pone.0048373 23110235

31. Lackford B, Yao C, Charles GM, Weng L, Zheng X, Choi EA, et al. Fip1 regulates mRNA alternative polyadenylation to promote stem cell self-renewal. Embo j. 2014;33(8):878–89. doi: 10.1002/embj.201386537 24596251

32. Di Giammartino DC, Li W, Ogami K, Yashinskie JJ, Hoque M, Tian B, et al. RBBP6 isoforms regulate the human polyadenylation machinery and modulate expression of mRNAs with AU-rich 3' UTRs. Genes Dev. 2014;28(20):2248–60. doi: 10.1101/gad.245787.114 25319826

33. Martinson HG. An active role for splicing in 3'-end formation. Wiley Interdiscip Rev RNA. 2011;2(4):459–70. doi: 10.1002/wrna.68 21957037

34. Niwa M, Rose SD, Berget SM. In vitro polyadenylation is stimulated by the presence of an upstream intron. Genes Dev. 1990;4(9):1552–9. 1701407

35. Kaida D, Berg MG, Younis I, Kasim M, Singh LN, Wan L, et al. U1 snRNP protects pre-mRNAs from premature cleavage and polyadenylation. Nature. 2010;468(7324):664–8. doi: 10.1038/nature09479 20881964

36. Gunderson SI, Polycarpou-Schwarz M, Mattaj IW. U1 snRNP inhibits pre-mRNA polyadenylation through a direct interaction between U1 70K and poly(A) polymerase. Mol Cell. 1998;1(2):255–64. 9659922

37. Almada AE, Wu X, Kriz AJ, Burge CB, Sharp PA. Promoter directionality is controlled by U1 snRNP and polyadenylation signals. Nature. 2013;499(7458):360–3. doi: 10.1038/nature12349 23792564

38. Ntini E, Jarvelin AI, Bornholdt J, Chen Y, Boyd M, Jorgensen M, et al. Polyadenylation site-induced decay of upstream transcripts enforces promoter directionality. Nat Struct Mol Biol. 2013;20(8):923–8. doi: 10.1038/nsmb.2640 23851456

39. Berg MG, Singh LN, Younis I, Liu Q, Pinto AM, Kaida D, et al. U1 snRNP determines mRNA length and regulates isoform expression. Cell. 2012;150(1):53–64. doi: 10.1016/j.cell.2012.05.029 22770214

40. Kyburz A, Friedlein A, Langen H, Keller W. Direct interactions between subunits of CPSF and the U2 snRNP contribute to the coupling of pre-mRNA 3' end processing and splicing. Mol Cell. 2006;23(2):195–205. 16857586

41. Millevoi S, Loulergue C, Dettwiler S, Karaa SZ, Keller W, Antoniou M, et al. An interaction between U2AF 65 and CF I(m) links the splicing and 3' end processing machineries. Embo J. 2006;25(20):4854–64. 17024186

42. Goraczniak R, Behlke MA, Gunderson SI. Gene silencing by synthetic U1 adaptors. Nat Biotechnol. 2009;27(3):257–63. doi: 10.1038/nbt.1525 19219028

43. Gruber AR, Martin G, Keller W, Zavolan M. Cleavage factor Im is a key regulator of 3' UTR length. RNA Biol. 2012;9(12):1405–12. doi: 10.4161/rna.22570 23187700

44. Richard P, Manley JL. How bidirectional becomes unidirectional. Nat Struct Mol Biol. 2013;20(9):1022–4. doi: 10.1038/nsmb.2657 24008561

45. Seila AC, Core LJ, Lis JT, Sharp PA. Divergent transcription: a new feature of active promoters. Cell Cycle. 2009;8(16):2557–64. 19597342

46. Seila AC, Calabrese JM, Levine SS, Yeo GW, Rahl PB, Flynn RA, et al. Divergent transcription from active promoters. Science. 2008;322(5909):1849–51. doi: 10.1126/science.1162253 19056940

47. Bresson SM, Conrad NK. The human nuclear poly(a)-binding protein promotes RNA hyperadenylation and decay. PLoS Genet. 2013;9(10):e1003893. doi: 10.1371/journal.pgen.1003893 24146636

48. Beaulieu YB, Kleinman CL, Landry-Voyer AM, Majewski J, Bachand F. Polyadenylation-dependent control of long noncoding RNA expression by the poly(A)-binding protein nuclear 1. PLoS Genet. 2012;8(11):e1003078. doi: 10.1371/journal.pgen.1003078 23166521

49. Hu J, Lutz CS, Wilusz J, Tian B. Bioinformatic identification of candidate cis-regulatory elements involved in human mRNA polyadenylation. RNA. 2005;11(10):1485–93. 16131587

50. Venkataraman K, Brown KM, Gilmartin GM. Analysis of a noncanonical poly(A) site reveals a tripartite mechanism for vertebrate poly(A) site recognition. Genes Dev. 2005;19(11):1315–27. 15937220

51. 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 J2004. p. 616–26. 14749727

52. Chan SL, Huppertz I, Yao C, Weng L, Moresco JJ, Yates JR 3rd, et al. CPSF30 and Wdr33 directly bind to AAUAAA in mammalian mRNA 3' processing. Genes Dev. 2014;28(21):2370–80. doi: 10.1101/gad.250993.114 25301780

53. Schonemann L, Kuhn U, Martin G, Schafer P, Gruber AR, Keller W, et al. Reconstitution of CPSF active in polyadenylation: recognition of the polyadenylation signal by WDR33. Genes Dev. 2014;28(21):2381–93. doi: 10.1101/gad.250985.114 25301781

54. Johnson SA, Kim H, Erickson B, Bentley DL. The export factor Yra1 modulates mRNA 3' end processing. Nat Struct Mol Biol. 2011;18(10):1164–71. doi: 10.1038/nsmb.2126 21947206

55. Hosoda N, Lejeune F, Maquat LE. Evidence that poly(A) binding protein C1 binds nuclear pre-mRNA poly(A) tails. Mol Cell Biol. 2006;26(8):3085–97. 16581783

56. Afonina E, Stauber R, Pavlakis GN. The human poly(A)-binding protein 1 shuttles between the nucleus and the cytoplasm. J Biol Chem. 1998;273(21):13015–21. 9582337

57. Wigington CP, Williams KR, Meers MP, Bassell GJ, Corbett AH. Poly(A) RNA-binding proteins and polyadenosine RNA: new members and novel functions. Wiley Interdiscip Rev RNA. 2014;5(5):601–22. doi: 10.1002/wrna.1233 24789627

58. Yao C, Biesinger J, Wan J, Weng L, Xing Y, Xie X, et al. Transcriptome-wide analyses of CstF64-RNA interactions in global regulation of mRNA alternative polyadenylation. Proc Natl Acad Sci U S A. 2012;109(46):18773–8. doi: 10.1073/pnas.1211101109 23112178

59. Zheng D, Tian B. RNA-binding proteins in regulation of alternative cleavage and polyadenylation. Adv Exp Med Biol. 2014;825:97–127. doi: 10.1007/978-1-4939-1221-6_3 25201104

60. Hoque M, Li W, Tian B. Accurate mapping of cleavage and polyadenylation sites by 3' region extraction and deep sequencing. Methods Mol Biol. 2014;1125:119–29. doi: 10.1007/978-1-62703-971-0_10 24590784

61. Yeo G, Burge CB. Maximum entropy modeling of short sequence motifs with applications to RNA splicing signals. J Comput Biol. 2004;11(2–3):377–94. 15700406

62. Ignatchenko V, Ignatchenko A, Sinha A, Boutros PC, Kislinger T. VennDIS: A JavaFX-based Venn and Euler Diagram Software to Generate Publication Quality Figures. Proteomics. 2014. E-pub ahead of print. doi: 10.1002/pmic.201400320.

Štítky
Genetika Reprodukční medicína

Článek vyšel v časopise

PLOS Genetics


2015 Číslo 4
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

plice
INSIGHTS from European Respiratory Congress
nový kurz

Současné pohledy na riziko v parodontologii
Autoři: MUDr. Ladislav Korábek, CSc., MBA

Svět praktické medicíny 3/2024 (znalostní test z časopisu)

Kardiologické projevy hypereozinofilií
Autoři: prof. MUDr. Petr Němec, Ph.D.

Střevní příprava před kolonoskopií
Autoři: MUDr. Klára Kmochová, Ph.D.

Všechny kurzy
Kurzy Podcasty Doporučená témata Časopisy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#