- Research article
- Open Access
Fine-tuning of microRNA-mediated repression of mRNA by splicing-regulated and highly repressive microRNA recognition element
BMC Genomics volume 14, Article number: 438 (2013)
MicroRNAs are very small non-coding RNAs that interact with microRNA recognition elements (MREs) on their target messenger RNAs. Varying the concentration of a given microRNA may influence the expression of many target proteins. Yet, the expression of a specific target protein can be fine-tuned by alternative cleavage and polyadenylation to the corresponding mRNA.
This study showed that alternative splicing of mRNA is a fine-tuning mechanism in the cellular regulatory network. The splicing-regulated MREs are often highly repressive MREs. This phenomenon was observed not only in the hsa-miR-148a-regulated DNMT3B gene, but also in many target genes regulated by hsa-miR-124, hsa-miR-1, and hsa-miR-181a. When a gene contains multiple MREs in transcripts, such as the VEGF gene, the splicing-regulated MREs are again the highly repressive MREs. Approximately one-third of the analysable human MREs in MiRTarBase and TarBase can potentially perform the splicing-regulated fine-tuning. Interestingly, the high (+30%) repression ratios observed in most of these splicing-regulated MREs indicate associations with functions. For example, the MRE-free transcripts of many oncogenes, such as N-RAS and others may escape microRNA-mediated suppression in cancer tissues.
This fine-tuning mechanism revealed associations with highly repressive MRE. Since high-repression MREs are involved in many important biological phenomena, the described association implies that splicing-regulated MREs are functional. A possible application of this observed association is in distinguishing functionally relevant MREs from predicted MREs.
MicroRNAs, which are very short non-coding RNAs of 21–23 nucleotides in length, play an important role in the regulation of gene expression, disease progression , development [2, 3], differentiation  and cancer . This regulatory process is mainly mediated by targeting microRNA recognition elements (MREs) in the 3′ untranslated region (3′UTR)  or the coding region [7, 8] of one or more messenger RNAs (mRNAs). As a result, either the targeted messenger RNAs are cleaved or the translation of transcripts is repressed .
Legendre et al.  hypothesized that some of the shorter MRE-free transcripts may escape from microRNA-mediated inhibition and further observed that alternative cleavage and polyadenylation (APA) reduces the amount of MRE-containing isoforms in the corresponding microRNA-expressed tissue. This finding is consistent with the observation that approximately two-thirds of the MRE-containing genes have isoforms in the 3′-UTRs . This hypothesis is also supported by findings that activation of primary murine CD4+ T lymphocytes is associated with an increased number of isoforms that lack MREs . Restated, these MRE-free isoforms have escaped microRNA-mediated inhibition. Furthermore, in cancer cells, oncogene activation is often accompanied by APA of mRNAs, which can cause widespread loss of 3′UTR repressive elements . The above reports reveal not only that APA is a regulatory mechanism of microRNA-mediated inhibition of protein expression, but also that a functional MRE may reside in 3′-UTR. Yet not all functional MREs are present in the 3′-UTR . Whether this inhibition can also be regulated by another mechanism is yet to be explored.
In human multi-exon genes, 50% ~ 95% of primary transcripts generate isoforms by alternative splicing (AS) or by APA [15–19]. These events greatly increase the functional complexity of the human genome. Some isoforms apparently have active roles in various cell types or tissues [20, 21]. Therefore, although alternative splicing is an effective mechanism for generating isoforms, the isoforms may lack functional elements in either the coding or noncoding regions of a primary transcript. A recent systemic study of the Arabidopsis thaliana plant in Yang et al.  found a significantly higher frequency of alternative splicing in MREs compared to other regions. However, the relationship between microRNA-mediated inhibition and alternative splicing has not been systematically studied in humans. Moreover, no studies have investigated whether splicing-regulated MRE is highly repressive or functional. Therefore, this study analyzed not only the relationship between alternative splicing and microRNA-mediated protein inhibition, but also the role of splicing-regulated MREs in protein repression.
Discovering tissues enriched with specific isoforms derived by alternative splicing
Expressed sequence tag (EST) sequences have proven useful as a reference for identifying alternative splicing in mRNA [23, 24]. When sufficient EST sequences are available, the number of EST sequences is approximately proportional to the level of gene expression . Therefore, Legendre et al.  used the number of EST sequences to determine the differential expression of isoforms contained in EST libraries. They further used this method to show that isoforms may or may not contain conserved regulatory motifs. Since many microRNA recognition elements (MREs) reveal evolutionary conservation , MRE-containing or MRE-free transcripts should be detected by a similar approach.
Even all EST sequences used in this study were obtained from non-normalized cDNA libraries, the EST approach suffers from several limitations. For example, not observing an EST sequence does not mean that a gene is not expressed. Thus, all the genes analyzed herein have at least 10 EST sequences. Furthermore, the 3′end of a messenger RNA usually has more EST sequences, so the probabilities obtained by Fisher exact test  tended to be underestimated. This bias favours the argument made herein concerning the observation of a differential expression of the counts of an isoform between a given tissue and all the other tissues. Therefore, multiple libraries and additional EST sequences are needed for further discovery of significant MREs. From this point of view, the single cell RNA sequencing approach  will be better than pooling the information from different libraries. However, the depth of the public RNA sequencing data is not high enough to do a systematic survey.
Figure 1 shows the example of an MRE present in seven EST sequences but absent from three other EST sequences in brain tissue. However, MRE is present in one EST sequences but absent from six other EST sequences in “all other tissues” (colon, lung, and liver in this example). The two EST sequences from the liver library are excluded since they do not include the entire MRE region. All the aforementioned numbers can be summarized in a two-by-two contingency Table (Figure 1IV and subjected to the Fisher exact test, (Additional file 1: Table S1-1A).
The DNA methyltransferase 3b (DNMT3B) gene is known to contain two putative MREs that are targeted by the human microRNA hsa-miR-148a . The location of these two MRE sites was found in the coding regions by aligning the reported MRE sequences with the reference sequence of UniGene cluster Hs.713611 (Figure 2A). Transcript DNMT3B3 contains only MRE site#2 (nucleotide 1775–1797), but transcript DNMT3B1 contains both MRE site#1 (nucleotide 2739–2767) and MRE site#2 (Figure 2A and Additional file 1: Figure S1-1). The site#1 MRE-free transcript (i.e. DNMT3B3) encodes a protein, which has an in-frame deletion of about 83 amino acids compared to the protein encoded by the site#1 MRE-containing transcript (i.e. DNMT3B1). This deletion removes the catalytic site of the methylase activity , but it is not clear whether this shorter protein still have the DNA 5-hydroxymethylcytosine dehydroxymethylase activity .
The site#1 MRE-free DNMT3B3 transcript may escape translational suppression by hsa-miR-148a, which is expressed in many tissues, such as human embryonic stem cell , brain , cervix  and other tissues . Therefore, the different isoforms transcribed from the DNMT3B gene might be a good starting point to explore the effect of alternative splicing on microRNA-mediated suppression of protein expression.
Since DNMT3B is abundantly expressed in ES cells and early embryos , alternative splicing effects can be studied in embryonic tissue. To test whether the highly repressive (repressive ratio ~50%) site#1 MRE  is regulated by an alternative splicing event, gene expression was compared between site#1 MRE-containing isoforms and site#1 MRE-free isoforms. Figure 2B shows that, in embryonic tissue, site#1 MRE-containing (i.e. DNMT3B1) isoforms are represented by more EST sequences compared to site#1 MRE-free isoforms (i.e. DNMT3B3). Left-sided Fisher exact test (see “Methods” for details) showed that the expression in embryonic tissue significantly differed from those of all other tissues (P = 0.0019). In other words, embryonic tissue expressed a higher proportion of site#1 MRE-containing transcript than that in all other tissues and might be more responsive to the hsa-miR-148a-mediated protein repression than that in all other tissues.
In contrast, site#2 MRE-free and site#2 MRE-containing isoforms showed no significant differences between “embryonic tissue” and “all other tissues” based on Fisher exact test. Since the proportion of site#2 MRE-containing transcripts was not significantly changed by alternative splicing event, site#2 MRE was not expected to have a great difference in hsa-miR-148a-mediated protein repression.
Splicing-regulated site#1 MRE in DNMT3B transcripts is highly repressive
To examine whether alternative splicing may regulate miRNA-mediated protein repression in other tissues, Fisher exact tests were performed in several tissues with sufficient EST sequences. Figure 2C and Additional file 1: Table S1-1B show that the statistically significant P-value obtained by right-side Fisher exact test for brain tissue indicated preferential expression of site#1 MRE-free transcript (i.e. DNMT3B3) in the brain. This result in brain tissue was opposite to that in embryonic tissue. The observations based on EST approach were also supported by the single cell RNA sequencing data . For example, the tag counts of those site#1 MRE-containing and site#1 MRE-free isoforms in human embryonic and brain tissues were consistent with the observations by EST approach qualitatively (Additional file 2). Accordingly, alternative splicing is apparently able to regulate protein expression in different tissues.
According to TarBase  and miRTarBase , hsa-miR-148a has one and seven experimentally verified target messenger RNAs, respectively. In the DNMT3B3 gene, site#1 MRE-free transcript may escape repression of hsa-miR-148a without affecting expression of other targeted genes in the brain. Restated, by changing the properties of a transcript rather than by changing the concentration of a given microRNA, alternative splicing may be an additional regulatory mechanism in protein expression. This additional regulatory mechanism may be useful in fine-tuning the complex gene expression circuit.
The biological importance of microRNA-mediated regulation is typically associated with a highly repressive MRE. Comparisons of microRNA–mediated repression in different mRNA isoforms has been successfully used to associate MRE sites with functions in many studies, such as the prediction of target mRNAs [37–41], the influence of intron retention on human mRNA , cellular proliferation and differentiation [12, 43], and cancer . The effect of fine-tuning at the splicing level would be negligible if the associated MREs are not highly repressive. Therefore, we predicted that MREs that regulated by splicing are also highly repressive. The DNMT3B analysis in this study showed that only site#1 MRE is highly complementary to hsa-miR-148a and is highly repressive (repressive ratio ~50%) of protein expression . However, site#2 MRE has little no effect (repressive ratio ~10%), which is consistent not only with the above argument, but also with selection pressure in the evolutionary process. That is, if a highly repressive MRE for fine-tuning the microRNA-mediated protein suppression provided a selection advantage, the highly repressive MRE would be preserved in evolution. Therefore, the alternative splicing-modulated protein repression is likely to be associated with highly repressive MRE. In contrast, a weak repressive MRE would confer a smaller selection advantage because of the limited effect of splicing-mediated protein suppression by microRNA. Thus, the weak repressive MRE might no longer be associated with the splicing-level regulation in the evolutionary process.
Proteomic evidence that highly repressive MREs are generally associated with alternative splicing events
Until now, the association between alternative splicing and highly repressive MRE has been supported only by a single microRNA and a single target messenger RNA. Whether this evidence is applicable to other target messenger RNAs is unclear. Furthermore, gene expression is typically, but not always, proportional to protein expression . Therefore, further studies are needed to determine whether alternative splicing-regulated MREs are generally highly repressive of protein expression.
The global impact of a microRNA, such as miR-124, on the protein expression in the HeLa cell has been measured by stable isotopic labelling with amino acids in cell culture (SILAC) analysis , a mass spectrometry technique for using non-radioactive isotopic labelling to detect differences in protein abundance among samples. The SILAC results can also be used to calculate protein expression-based repression. The analysis herein identified 1,544 differentially expressed proteins whose expression resulted from expression of hsa-miR-124 in HeLa cells. However, only 1,486 have symbols in UniGene#230 (and are referred to as the “Proteomics” fraction, Additional file 1: Table S1-2A). Approximately half of these 1,486 proteins were repressed and, of these, approximately half were activated (Additional file 1: Table S1-2B and Text S1-1). These 1,486 proteins identified by proteomic analysis of HeLa cells were used as the starting point in a global analysis.
These proteins include the direct targets of hsa-miR-124 and the downstream proteins of the hsa-miR-124-regulated target genes . Target proteins were distinguished from downstream proteins by labelling the putative target messenger RNA as “seed”, which is a 7-nucleotide feature typically found at the 3′-end of an MRE and complementary to the 5′-end of a microRNA. The presence of a “seed” is thus considered evidence of a candidate MRE. Of these 1,486 proteins, the presence of the seed” region on the transcript indicated that only 408 proteins are directly regulated by hsa-miR-124 (Proteomics∩Seed(ALL), red line in Figure 3). Since the repression is too weak for depiction by histogram, the variation in repression in different states was compared by plotting a cumulative curve [12, 19, 46, 48] of the percentage of proteins in the given range of repression ratios. Figure 3 plots the cumulative curves based on the percentage of counted genes according to the full-length “seed” matches presented in Additional file 1: Table S1-2B and the results of the paired statistical analysis in Additional file 1: Table S1-2C.
The remaining 1,078 seed-free detected proteins were labelled “Proteomics∩Seed(X)” (blue line in Figure 3). The red line showed a significant left shift relative to the blue line (P < 0.001, in a Kolmogorov-Smirnov (KS) test with two independent samples, Additional file 1: Table S1-2C), indicating that repression was greater in the seed-annotated proteins compared to other differentially expressed proteins. However, the presence of a seed usually, but not always correlates with an MRE . Therefore, these “seed”-containing transcripts were further evaluated by the miRSVR algorithm  to predict the potential target of hsa-miR-124. These miRSVR-predicted MREs are considered putative MREs. In Figure 3, the cumulative curve of (Proteomics∩Seed(ALL)∩miRSVR) indicated by the purple line is significantly shifted to the left of the blue curve (P < 0.001, KS test, Additional file 1: Table S1-2C), indicating that the target prediction method is better than seed-free detected target messenger RNA. Alternative splicing was then analyzed in target genes (Proteomics∩Seed(ALL)∩miRSVR) encoding differentially expressed proteins. Additional file 1: Table S1-2 presents those proteins that were potentially controlled at the splicing level (Proteomics∩Seed(ALL)∩miRSVR∩AS). The cumulative curve of these splicing-regulated target proteins (green line in Figure 3) exhibited the highest repression, because the green curve shifted further to the left of the blue curve (seed-free target proteins) (P < 0.001, KS test, Additional file 1: Table S1-2C).
In addition to hsa-miR-124, the SILAC analysis  was also used to examine the on the global impact of hsa-miR-1 and hsa-miR-181a on protein expression in HeLa cell. This global trend has also been reported in studies of hsa-miR-1 and hsa-miR-181a (Additional file 3), which are consistent with the previous observation on hsa-miR-124. In summary, the ratio of repressive proteins was apparently higher in the alternative splicing-regulated target proteins compared to the remaining differentially expressed proteins. This global analysis of multiple target proteins of a single microRNA supports the association between alternative splicing and highly repressive MRE.
Splicing-regulated MREs are highly repressive in transcripts that contain multiple MREs
The DNMT3B gene is a good example of the highly repressive splicing-regulated MRE of hsa-miR-148a. The repressive effect is also observed in the target genes of hsa-miR-124 at the protein level. However, many genes are targeted by multiple miRNAs in the cells. For example, approximately 30 putative miRNAs regulate the vascular endothelial growth factor (VEGF) gene . To see the repressive effect and for a fair comparison, different MREs but with at least 90% overlap were grouped into a single “MRE region”. After alignment of the reported MRE regions with UniGene cluster Hs.73793, the splicing-regulated MREs were identified as described above in relation to Figure 1. Thirteen putative splicing-regulated MRE regions were identified (Figure 4) by comparing the locations of MRE regions with splicing sites. The repressive ratio data for these sites were collected in hypoxia-induced CNE cells (a nasopharyngeal carcinoma cell line) by introducing synthetic putative VEGF-regulative microRNA duplexes to this cell line  (Figure 4, bottom panel). Of these 13 putative splicing-regulated MRE regions, nine (approximately 70%) were considered highly repressive MRE regions (repressive ratios of +30%): regions 1, 2, 3, 4, 5, 7, 9, 11 and 12 (Table 1, Figure 4, bottom panel).
The cancer (top) and normal (middle) panels in Figure 4 plot the negative log10 transformed P values of all MREs. Four MRE regions showed significant splicing regulation (P < 0.05) in at least one cancer tissue. Two MRE regions appeared to be splicing-regulated in at least one normal tissue. These splicing-regulated MREs were marked with red stars in the two-tail test panels in Figure 4. The normal and cancerous tissues contained five splicing-regulated MRE regions. All five of the splicing-regulated MRE regions (according to Fisher exact test) are also highly repressive MRE regions even though not all highly repressive MREs are regulated by alternative splicing. This implies that analysis of splicing regulation may reveal highly repressive or “true” MREs.
Approximately one-third of experimentally confirmed MREs are associated with alternative splicing events
Table 1 shows that, in many miRNAs that target VEGF messenger RNA, splicing-regulated MREs were associated with highly repressive MREs. Despite the importance of this observation, this phenomenon is not biologically significant if splicing-regulated MREs are special cases. Therefore, alternative splicing analysis was repeated in experimentally supported human targets based on the reporter assay results . Random MREs were also analyzed to evaluate background effects. Notably, the true MRE was positively validated by the algorithm, the false MRE was negatively validated by the algorithm , and the random MRE was simply randomly selected from reference sequences. If our observation of DNMT3B is widely applicable, then the opportunity to observed splicing-regulated MRE would be greater for the experimentally verified MREs than for the false or random MREs.
The sequences of translationally repressed human MREs and false MREs were retrieved from both miRTarBase release 2.5  and TarBase v.4.0  and aligned to the corresponding reference sequences in UniGene database (see Methods). Only MREs with at least ten EST sequences were included in further analyses. The analysis included 256 unique MRE sites obtained from miRTarBase and 72 unique MRE sites and 14 unique false sites obtained from TarBase. From the 3′UTR and full-length sequences, 170,000 and 431,500 distinct random MREs were generated, respectively.
The number of unique MRE sites was determined by applying the following two rules. First, an MRE targeted by multiple microRNAs (e.g., hsa-miR-20a and hsa-miR-17 both targeting the same MRE in the E2F1 gene) was considered a single site. Second, significant tissue types (Additional file 1: Figure S1-2) identified in the alternative splicing analysis of a given MRE were considered OR significant tissue types (Additional file 1: Figure S1-2) identified in the alternative splicing analysis of a given MRE were defined as splicing-regulated. Alternative splicing in at least one tissue was identified in 95 out of 256 (37.1%) MRE sites in miRTarBase and in 22 out of 72 (30.6%) MRE sites in TarBase (Figure 5). These percentages were 2.2 and 2.6 times those of the false MRE sites in TarBase and miRTarBase, respectively. The percentages of both the true and false MRE sites significantly differed (average, 9.7%) from those of the random MRE sites selected from 3′ UTR or full-length sequence (P < 0.001, one-sample Wilcoxon signed-rank test). A possible explanation for the difference is that negative validation of the algorithm-predicted MRE has been reported for only one or a few tissues, which may significantly differ from other tissues in the EST libraries. The percentage of significant randomly selected MREs in this analysis (9.7%) also differed from those in Additional file 1: Figure S1-2B (2.8%) because the latter treated MREs as independent in different tissues.
The experimental results show that approximately one-third of the true MREs were associated with AS events. Additionally, reported repressive ratios exceeded 30% in almost all proteins in Figure 6. Restated, the association between splicing-regulated MREs and highly repressive MREs is again supported by other microRNAs and their target messenger RNAs.
This study showed that splicing-level regulation is a fine-tuning mechanism in the expressions of specific proteins. The fine-tuning process described above enables regulation of a single protein in the presence of all the other target proteins of a given microRNA. Briefly, specific examples of this process were extended to all analysable human targets in TarBase and miRTarBase. The microRNA-mediated regulation is quite complex. For example, a single microRNA may have multiple target genes, and a messenger RNA may be targeted by multiple microRNAs. Alternative splicing enables independent regulation of protein expression in a single gene. An MRE-containing isoform can be co-regulated with other target genes, but an MRE-free isoform may escape regulation of the matching microRNA.
Additionally, splicing-regulated MREs established here were highly repressive (repressive ratios +30%). This highly repressive character may enhance the effectiveness of alternative splicing. The fine-tuning process may exert a positive selection pressure in evolution. This effect may explain why almost all observed splicing-regulated MREs were highly repressive. However, not all genes require fine-tuning. Therefore, highly repressive MREs may not necessarily be regulated by alternative splicing. For example, the highly repressive (repressive ratio 75%) hsa-miR-137-specific MRE in the CDC42 gene  is not associated with alternative splicing (data not shown).
The discussed fine-tuning mechanism also suggests a means for cancer cells to escape regulation in a normal cell. As presented in Figure 6, enrichment of splicing-regulated MRE regions was significantly increased in cancer cells according to right-sided Fisher exact test. This right-side significance indicates an increased amount of an MRE-free transcript in cancer cells. Restated, genes such as CSDE1, ERBB2, ERBB3, GJA1, N-RAS, SLC16A, SMAD1, TMSB4X, TMSL3, TTK, and ZNF513 were up-regulated in cancer cells. This result is consistent with the notion that a loss of repressive elements located in 3′ untranslated regions (3′UTRs) by alternative cleavage and polyadenylation (APA) was associated with oncogene activation  and cellular activation . Thus, the above association is not only a common phenomenon, it also has important biological implications.
Regarding applications, identifying functional MREs is critical to determining the biological function of microRNAs. However, measuring the repressive ratio at an MRE site is laborious. Although the prediction algorithm considered many factors, the false positive rate was unacceptably high [3, 50, 51]. Therefore, predicting the functional MREs of a given microRNA remains challenging [52, 53]. Our experimental results suggest that alternative splicing can also be used to filter and identify functional MRE sites from a list of putative sites predicted by microRNA target identification algorithms. The false-negative rate of the suggested method of identifying novel targets of miRNAs is high since approximately two-thirds of the true MREs were unassociated with AS events. Although this method is not comprehensive, the false positive rate is low. More importantly, most of these splicing-regulated MREs are highly repressive MREs, which have functional implications.
In conclusion, the discovery of splicing-regulated MREs reveals an important fine-tuning mechanism in complex microRNA-mediated regulation. This phenomenon is widely observed in human cells and is accompanied by the strong repression of protein expression. This novel and biologically significant observation is useful for identifying MREs, which may have important biological functions.
EST data source and prediction of spliced sites
Publically available Expressed Sequence Tag (EST) sequences were used to discover isoforms of messenger RNA [24, 54, 55]. Alternatively spliced isoforms can be identified by aligning EST sequences with reference messenger RNA sequences. Tissues that express a particular EST sequence can also be annotated with NCBI information. The MRE-containing or MRE-free isoforms can then be identified by comparing the locations of functional MREs with alternative splicing (AS) sites. The UniGene (Human, Build #230) database was used to discover the alternative splicing sites by aligning each EST sequence with the reference sequence, which was the longest high-quality sequence in the NCBI (file “Hs.seq.uniq”). Alignment was performed in sim4 using the default parameters  as described by Huang, et al. . The relationship between a putative splicing site and an MRE site (see below for details) was determined by comparing the locations of the above two features.
Alignment of experimentally verified and predicted MREs with the reference sequence
The known microRNA genes were obtained from miRBase (version 17.0)  and the experimentally verified MREs were downloaded or retrieved from TarBase (version 4.0)  and miRTarBase (release 2.5) . For each predicted MRE, information obtained from microRNA.org (August, 2010 Release)  included potential binding sites identified by using the mirSVR  and miRanda  algorithms for a given microRNA. The MRE positions were then identified by using the BLASTN algorithm  to match each MRE with the reference sequences and then retaining only those with perfect matches over their full length.
Configuration of EST sequence and Fisher exact test
As in step (III) of Figure 1, information about alternative splicing, tissue, histology (normal vs. cancer), and MRE were integrated. Tissue and histology information were obtained from the EST library established by the Cancer Genome Anatomy Project [60, 61]. After the splicing and MRE information were integrated, each EST sequence had five possible configurations in relation to a specified MRE sequence: (i) partial or no overlap with MRE; (ii) absence of MRE in a putative exon; (iii) presence of MRE in a putative intron; (iv) partial MRE at the exon-intron junction, and (v) other. Each EST sequence was assigned a configuration number in this step (e.g., see the parenthesis after a given tissue name in step (III) of Figure 1) and was classified accordingly. Each configuration had a unique biological meaning. In configuration (i), the EST sequence was unrelated to alternative splicing. In configuration (ii), the EST sequence was in an MRE-containing isoform. In configuration (iii), the EST sequence was in an MRE-free isoform. In configuration (iv), the isoform contained only a partial, and potentially non-functional, MRE. Therefore, the isoform was classified as an MRE-free isoform.
The EST sequences were grouped according to whether they were MRE-containing or MRE-free isoforms in various tissues and histologies. In step (IV) (Figure 1), a two-way contingency table was established to determine whether the counts of MRE-free or MRE-containing isoforms were over-represented in a given tissue, in relation to “all other tissues”. The statistical analyses were limited to tissues with at least ten EST sequences. The P-value was computed by Fisher exact test (Figure 1 and Additional file 1: Table S1-1A) to identify non-random associations between two categorical variables. The null hypothesis was that, in the absence of natural selection, the proportion of (MRE-free isoforms/MRE-containing isoforms) in a given tissue is identical to that in all other tissues. Therefore, the hypergeometric distribution was used to calculate the probabilities of the observed data and all data sets with more extreme deviations. When a given tissue showed enriched MRE-containing isoforms, the left-side P-value was significant in one-tailed Fisher exact test. In contrast, when a given tissue showed increased MRE-free isoforms, the right side P-value was significant in a one-tailed Fisher exact test.
Preparing random sequences for Wilcoxon signed-rank test
To observe the background effect of the configuration in Figure 1, 170,000 and 431,500 different MREs were randomly generated from 3′UTR and full-length sequences, respectively. The numbers of random MREs sufficed for the purposes of this study since 2.8% of 100,000 random MREs were putative splicing-regulated MREs (P < =0.05, see red line in Additional file 1: Figure S1-2B indicating accumulated probability of all MREs randomly selected for statistical analysis). These reference sequences were randomly selected from the UniGene database for use as source sequences. The sequence source, the starting position of a random sequence, and the length (18 ~ 22 nucleotides) of a random sequence were all randomized. Once the starting site had been determined, an MRE sequence was randomly generated from the source sequence under the constraint of MRE length. Therefore, all selected MREs were in the source sequences and were included in the isoform enrichment analysis described above. Since the amounts of splicing-regulated random MREs were not normally distributed among these two sets (Figure 5, right panel), one-sample Wilcoxon signed-rank test was used to differentiate between true MREs and randomly generated MREs.
Kuchenbauer F, Morin RD, Argiropoulos B, Petriv OI, Griffith M, Heuser M, Yung E, Piper J, Delaney A, Prabhu AL, et al: In-depth characterization of the microRNA transcriptome in a leukemia progression model. Genome Res. 2008, 18 (11): 1787-1797. 10.1101/gr.077578.108.
Kloosterman WP, Plasterk RH: The diverse functions of microRNAs in animal development and disease. Dev Cell. 2006, 11 (4): 441-450. 10.1016/j.devcel.2006.09.009.
Zhu QH, Spriggs A, Matthew L, Fan L, Kennedy G, Gubler F, Helliwell C: A diverse set of microRNAs and microRNA-like small RNAs in developing rice grains. Genome Res. 2008, 18 (9): 1456-1465. 10.1101/gr.075572.107.
Bland CS, Cooper TA: Micromanaging alternative splicing during muscle differentiation. Dev Cell. 2007, 12 (2): 171-172. 10.1016/j.devcel.2007.01.014.
Blenkiron C, Miska EA: miRNAs in cancer: approaches, aetiology, diagnostics and therapy. Hum Mol Genet. 2007, 16: R106-R113. 10.1093/hmg/ddm056. Spec No 1
Lee RC, Feinbaum RL, Ambros V: The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993, 75 (5): 843-854. 10.1016/0092-8674(93)90529-Y.
Duursma AM, Kedde M, Schrier M, le Sage C, Agami R: miR-148 targets human DNMT3b protein coding region. RNA. 2008, 14 (5): 872-877. 10.1261/rna.972008.
Forman JJ, Legesse-Miller A, Coller HA: A search for conserved sequences in coding regions reveals that the let-7 microRNA targets Dicer within its coding sequence. Proc Natl Acad Sci USA. 2008, 105 (39): 14879-14884. 10.1073/pnas.0803230105.
Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5.
Legendre M, Ritchie W, Lopez F, Gautheret D: Differential repression of alternative transcripts: a screen for miRNA targets. PLoS Comput Biol. 2006, 2 (5): e43-10.1371/journal.pcbi.0020043.
Majoros WH, Ohler U: Spatial preferences of microRNA targets in 3′ untranslated regions. BMC Genomics. 2007, 8: 152-10.1186/1471-2164-8-152.
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-1647. 10.1126/science.1155390.
Mayr C, Bartel DP: Widespread shortening of 3′UTRs by alternative cleavage and polyadenylation activates oncogenes in cancer cells. Cell. 2009, 138 (4): 673-684. 10.1016/j.cell.2009.06.016.
Schnall-Levin M, Zhao Y, Perrimon N, Berger B: Conserved microRNA targeting in Drosophila is as widespread in coding regions as in 3′UTRs. Proc Natl Acad Sci USA. 2010, 107 (36): 15751-15756. 10.1073/pnas.1006172107.
Ben-Dov C, Hartmann B, Lundgren J, Valcarcel J: Genome-wide analysis of alternative pre-mRNA splicing. J Biol Chem. 2008, 283 (3): 1229-1233.
Xing Y, Lee C: Relating alternative splicing to proteome complexity and genome evolution. Adv Exp Med Biol. 2007, 623: 36-49. 10.1007/978-0-387-77374-2_3.
Lemischka IR, Pritsker M: Alternative splicing increases complexity of stem cell transcriptome. Cell Cycle. 2006, 5 (4): 347-351. 10.4161/cc.5.4.2424.
Singh P, Alley TL, Wright SM, Kamdar S, Schott W, Wilpan RY, Mills KD, Graber JH: Global changes in processing of mRNA 3′ untranslated regions characterize clinically distinct cancer subtypes. Cancer Res. 2009, 69 (24): 9422-9430. 10.1158/0008-5472.CAN-09-2236.
Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, Kingsmore SF, Schroth GP, Burge CB: Alternative isoform regulation in human tissue transcriptomes. Nature. 2008, 456 (7221): 470-476. 10.1038/nature07509.
Venables JP: Aberrant and alternative splicing in cancer. Cancer Res. 2004, 64 (21): 7647-7654. 10.1158/0008-5472.CAN-04-1910.
Wang Z, Lo HS, Yang H, Gere S, Hu Y, Buetow KH, Lee MP: Computational analysis and experimental validation of tumor-associated alternative RNA splicing in human cancer. Cancer Res. 2003, 63 (3): 655-657.
Yang X, Zhang H, Li L: Alternative mRNA processing increases the complexity of microRNA-based gene regulation in Arabidopsis. Plant J: for cell and molecular biology. 2012, 70 (3): 421-431. 10.1111/j.1365-313X.2011.04882.x.
Xu Q, Modrek B, Lee C: Genome-wide detection of tissue-specific alternative splicing in the human transcriptome. Nucleic Acids Res. 2002, 30 (17): 3754-3766. 10.1093/nar/gkf492.
Huang YH, Chen YT, Lai JJ, Yang ST, Yang UC: PALS db: Putative Alternative Splicing database. Nucleic Acids Res. 2002, 30 (1): 186-190. 10.1093/nar/30.1.186.
Ewing RM, Ben Kahla A, Poirot O, Lopez F, Audic S, Claverie JM: Large-scale statistical analyses of rice ESTs reveal correlated patterns of gene expression. Genome Res. 1999, 9 (10): 950-959. 10.1101/gr.9.10.950.
Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005, 120 (1): 15-20. 10.1016/j.cell.2004.12.035.
Ramskold D, Luo S, Wang YC, Li R, Deng Q, Faridani OR, Daniels GA, Khrebtukova I, Loring JF, Laurent LC, et al: Full-length mRNA-Seq from single-cell levels of RNA and individual circulating tumor cells. Nat Biotechnol. 2012, 30 (8): 777-782. 10.1038/nbt.2282.
Weisenberger DJ, Velicescu M, Cheng JC, Gonzales FA, Liang G, Jones PA: Role of the DNA methyltransferase variant DNMT3b3 in DNA methylation. Molecular Cancer Res: MCR. 2004, 2 (1): 62-72.
Chen CC, Wang KY, Shen CK: The mammalian de novo DNA methyltransferases DNMT3A and DNMT3B are also DNA 5-hydroxymethylcytosine dehydroxymethylases. J Biol Chem. 2012, 287 (40): 33116-33121. 10.1074/jbc.C112.406975.
Tzur G, Levy A, Meiri E, Barad O, Spector Y, Bentwich Z, Mizrahi L, Katzenellenbogen M, Ben-Shushan E, Reubinoff BE, et al: MicroRNA expression patterns and function in endodermal differentiation of human embryonic stem cells. PLoS One. 2008, 3 (11): e3726-10.1371/journal.pone.0003726.
Hua D, Mo F, Ding D, Li L, Han X, Zhao N, Foltz G, Lin B, Lan Q, Huang Q: A catalogue of glioblastoma and brain MicroRNAs identified by deep sequencing. Omics: a journal of integrative biology. 2012, 16 (12): 690-699. 10.1089/omi.2012.0069.
Witten D, Tibshirani R, Gu SG, Fire A, Lui WO: Ultra-high throughput sequencing-based small RNA discovery and discrete statistical biomarker analysis in a collection of cervical tumours and matched controls. BMC Biol. 2010, 8: 58-10.1186/1741-7007-8-58.
Hsu SD, Chu CH, Tsou AP, Chen SJ, Chen HC, Hsu PW, Wong YH, Chen YH, Chen GH, Huang HD: miRNAMap 2.0: genomic maps of microRNAs in metazoan genomes. Nucleic Acids Res. 2008, 36 (Database issue): D165-D169.
Okano M, Bell DW, Haber DA, Li E: DNA methyltransferases Dnmt3a and Dnmt3b are essential for de novo methylation and mammalian development. Cell. 1999, 99 (3): 247-257. 10.1016/S0092-8674(00)81656-6.
Sethupathy P, Corda B, Hatzigeorgiou AG: TarBase: A comprehensive database of experimentally supported animal microRNA targets. RNA. 2006, 12 (2): 192-197.
Hsu SD, Lin FM, Wu WY, Liang C, Huang WC, Chan WL, Tsai WT, Chen GZ, Lee CJ, Chiu CM, et al: miRTarBase: a database curates experimentally validated microRNA-target interactions. Nucleic Acids Res. 2011, 39 (Database issue): D163-D169.
Doench JG, Sharp PA: Specificity of microRNA target selection in translational repression. Genes Dev. 2004, 18 (5): 504-511. 10.1101/gad.1184404.
Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB: Prediction of mammalian microRNA targets. Cell. 2003, 115 (7): 787-798. 10.1016/S0092-8674(03)01018-3.
Robins H, Li Y, Padgett RW: Incorporating structure to predict microRNA targets. Proc Natl Acad Sci USA. 2005, 102 (11): 4006-4009. 10.1073/pnas.0500775102.
Long D, Lee R, Williams P, Chan CY, Ambros V, Ding Y: Potent effect of target structure on microRNA function. Nat Struct Mol Biol. 2007, 14 (4): 287-294. 10.1038/nsmb1226.
Kozomara A, Griffiths-Jones S: miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res. 2011, 39 (Database issue): D152-D157.
Tan S, Guo J, Huang Q, Chen X, Li-Ling J, Li Q, Ma F: Retained introns increase putative microRNA targets within 3′ UTRs of human mRNA. FEBS Lett. 2007, 581 (6): 1081-1086. 10.1016/j.febslet.2007.02.009.
Yeo GW, Coufal N, Aigner S, Winner B, Scolnick JA, Marchetto MC, Muotri AR, Carson C, Gage FH: Multiple layers of molecular controls modulate self-renewal and neuronal lineage specification of embryonic stem cells. Hum Mol Genet. 2008, 17 (R1): R67-R75. 10.1093/hmg/ddn065.
Johnson SM, Grosshans H, Shingara J, Byrom M, Jarvis R, Cheng A, Labourier E, Reinert KL, Brown D, Slack FJ: RAS is regulated by the let-7 microRNA family. Cell. 2005, 120 (5): 635-647. 10.1016/j.cell.2005.01.014.
de Sousa Abreu R, Penalva LO, Marcotte EM, Vogel C: Global signatures of protein and mRNA expression levels. Mol Biosyst. 2009, 5 (12): 1512-1526.
Baek D, Villen J, Shin C, Camargo FD, Gygi SP, Bartel DP: The impact of microRNAs on protein output. Nature. 2008, 455 (7209): 64-71. 10.1038/nature07242.
Li X, Jiang W, Li W, Lian B, Wang S, Liao M, Chen X, Wang Y, Lv Y, Yang L: Dissection of human MiRNA regulatory influence to subpathway. Brief Bioinform. 2012, 13 (2): 175-186. 10.1093/bib/bbr043.
Selbach M, Schwanhausser B, Thierfelder N, Fang Z, Khanin R, Rajewsky N: Widespread changes in protein synthesis induced by microRNAs. Nature. 2008, 455 (7209): 58-63. 10.1038/nature07228.
Betel D, Koppal A, Agius P, Sander C, Leslie C: Comprehensive modeling of microRNA targets predicts functional non-conserved and non-canonical sites. Genome Biol. 2010, 11 (8): R90-10.1186/gb-2010-11-8-r90.
Ye W, Lv Q, Wong CK, Hu S, Fu C, Hua Z, Cai G, Li G, Yang BB, Zhang Y: The effect of central loops in miRNA: MRE duplexes on the efficiency of miRNA-mediated gene regulation. PLoS One. 2008, 3 (3): e1719-10.1371/journal.pone.0001719.
Didiano D, Hobert O: Molecular architecture of a miRNA-regulated 3′ UTR. RNA. 2008, 14 (7): 1297-1317. 10.1261/rna.1082708.
Sethupathy P, Megraw M, Hatzigeorgiou AG: A guide through present computational approaches for the identification of mammalian microRNA targets. Nat Methods. 2006, 3 (11): 881-886. 10.1038/nmeth954.
Barbato C, Arisi I, Frizzo ME, Brandi R, Da Sacco L, Masotti A: Computational challenges in miRNA target predictions: to be or not to be a true target?. J Biomed Biotechnol. 2009, 2009: 803069-
Brett D, Pospisil H, Valcarcel J, Reich J, Bork P: Alternative splicing and genome complexity. Nat Genet. 2002, 30 (1): 29-30. 10.1038/ng803.
Chang YM, Juan HF, Lee TY, Chang YY, Yeh YM, Li WH, Shih AC: Prediction of human miRNAs using tissue-selective motifs in 3′ UTRs. Proc Natl Acad Sci USA. 2008, 105 (44): 17061-17066. 10.1073/pnas.0809151105.
Florea L, Hartzell G, Zhang Z, Rubin GM, Miller W: A computer program for aligning a cDNA sequence with a genomic DNA sequence. Genome Res. 1998, 8 (9): 967-974.
Betel D, Wilson M, Gabow A, Marks DS, Sander C: The microRNA.org resource: targets and expression. Nucleic Acids Res. 2008, 36 (Database issue): D149-D153.
John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS: Human MicroRNA targets. PLoS Biol. 2004, 2 (11): e363-10.1371/journal.pbio.0020363.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.
Wheeler DL, Barrett T, Benson DA, Bryant SH, Canese K, Chetvernin V, Church DM, Dicuccio M, Edgar R, Federhen S, et al: Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2008, 36 (Database issue): D13-D21.
Hess JL: The Cancer Genome Anatomy Project: power tools for cancer biologists. Cancer Invest. 2003, 21 (2): 325-326. 10.1081/CNV-120016428.
The authors would like to thank Dr. Chun-Houh Chen of Institute of Statistical Science, Academia Sinica; Dr. Chung-Cheng Liu, Dr. Ling-Mei Wang and other members of BDL, ITRI for stimulating discussions. Ted Knoy is appreciated for his editorial assistance.
CTW, CYC and HCC were partially supported by NSC99-3112-B-010-019 and NSC100-2319-B-010-002 (Taiwan Bioinformatics Consortium core facility of the National Core Facility Program for Biotechnology, National Science Council). The computing facility used by this work was supported by NSC100-2325-B-002-065 and a grant from Ministry of Education, Aim for the Top University Plan.
The authors declare that they have no competing interests.
CTW designed and performed the major analysis and wrote the first draft of this manuscript; CYC provided the alternative splicing data; HCC performed the single-cell level mRNA-seq analysis; UCY reorganized and revised the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 3: Global trend of target proteins that were repressed or activated after transfection of hsa-miR-1 or hsa-miR-181a to HeLa cells.(DOC 3 MB)
Authors’ original submitted files for images
About this article
Cite this article
Wu, CT., Chiou, CY., Chiu, HC. et al. Fine-tuning of microRNA-mediated repression of mRNA by splicing-regulated and highly repressive microRNA recognition element. BMC Genomics 14, 438 (2013). https://doi.org/10.1186/1471-2164-14-438
- Alternative splicing
- MicroRNA recognition element
- Repressive ratio