- Research article
- Open Access
Evolutionarily emerged G tracts between the polypyrimidine tract and 3′ AG are splicing silencers enriched in genes involved in cancer
BMC Genomics volume 15, Article number: 1143 (2014)
The 3′ splice site (SS) at the end of pre-mRNA introns has a consensus sequence (Y)nNYAG for constitutive splicing of mammalian genes. Deviation from this consensus could change or interrupt the usage of the splice site leading to alternative or aberrant splicing, which could affect normal cell function or even the development of diseases. We have shown that the position “N” can be replaced by a CA-rich RNA element called CaRRE1 to regulate the alternative splicing of a group of genes.
Taking it a step further, we searched the human genome for purine-rich elements between the -3 and -10 positions of the 3′ splice sites of annotated introns. This identified several thousand such 3′SS; more than a thousand of them contain at least one copy of G tract. These sites deviate significantly from the consensus of constitutive splice sites and are highly associated with alterative splicing events, particularly alternative 3′ splice and intron retention. We show by mutagenesis analysis and RNA interference that the G tracts are splicing silencers and a group of the associated exons are controlled by the G tract binding proteins hnRNP H/F. Species comparison of a group of the 3′SS among vertebrates suggests that most (~87%) of the G tracts emerged in ancestors of mammals during evolution. Moreover, the host genes are most significantly associated with cancer.
We call these elements together with CaRRE1 regulatory RNA elements between the Py and 3′AG (REPA). The emergence of REPA in this highly constrained region indicates that this location has been remarkably permissive for the emergence of de novo regulatory RNA elements, even purine-rich motifs, in a large group of mammalian genes during evolution. This evolutionary change controls alternative splicing, likely to diversify proteomes for particular cellular functions.
Alternative splicing is the major source of proteomic diversity in vertebrates [1–3]. It occurs in more than 90% of human genes and plays important roles in cell function and the development of diseases [4–8]. The usage of alternative exons is tightly regulated in cells in a spatial and temporal manner by cis-acting RNA elements and trans-acting factors [9, 10]. Cis-acting elements have been identified around or overlapping with motifs of the splice sites. It is however very rare for a de novo regulatory element to exist between the polypyrimidine tract and 3′AG of the 3′ splice site (3′SS), apparently due to the constrained sequence and space of the motifs for proper splicing [11–13].
The 3′ splice sites at the end of introns have a highly conserved arrangement of consensus sequence (Y)nNYAG , where the polypyrimidine tract (Py) is close to the 3′ AG for the binding of the heterodimmeric U2AF65 and U2AF35 , respectively. The space between the (Y)n and 3′ AG is highly constrained: simply increasing the distance weakens splice site usage .
We have shown previously that a CA-rich element called CaRRE1 is within this 3′SS space and bound by hnRNP L to regulate depolarization-induced splicing [15–20], but it is not clear whether a purine-rich element could be tolerated at this location. One group of the well-characterized purine-rich splicing regulatory elements are guanidine (G) tracts containing a minimal functional G3, mostly as silencers in exons or enhancers in introns [21–23]. Here we searched the human genome for such potential elements and report the identification and characterization of a large group of G tracts between the (Y)n and 3′ AG, which we call with other elements at this special location REPA (regulatory RNA elements between the Py and 3′AG).
REPA G tracts within a group of human 3′SS
To determine if there are purine-rich REPA like the CaRRE1, we searched the 3′SS of annotated introns in the human genome for those with >60% A/G content within the -10 ~ -3 nucleotides of the intron end (Figure 1A-B), a length sufficient to harbor a splicing regulatory RNA element. This identified 5,041 unique 3′SS of 3′AG introns. On average, the A/G content between the -10 and -5 positions is about 70% compared to only about 20% for that of constitutive exons (Figure 1B-C). At the -4 and -3 positions, the A/G content is also higher, about 40% over the constitutive ones. Overall, the A/G content between the -10 and -3 positions is significantly higher (p = 0 for each position in hypergeometric test) than that of constitutive exons.
To identify consensus elements, we examined the nucleotides -15 to -3 nt to allow for a maximal 5 upstream nucleotides to form elements with the search target sequence. These sequences were subject to the MEME (multiple Em for motif elicitation) analysis as reported previously [24, 25]. This identified G-rich motifs, among others (Figure 2), as one of the top consensus elements.
The REPA G-tracts are highly associated with alternative splicing
To determine the strength of the splice sites containing the purine-rich elements, we calculated their MaxEntScan entropy scores, as described by Yeo and Burge . The scores are significantly lower than the 3847 3′SS randomly chosen from the human genome, with those containing the G pentamer and GGAAA among the lowest (Figure 2A, p = 3.8E-06 ~ 3.8E-39). Therefore, the purine-rich REPA likely weakens the 3′SS.To determine whether these 3′SS are associated with alternative splicing, we examined some of the REPA-containing 3′SS and their downstream exons of human mRNA or expressed sequence tags (EST) in the UCSC Genome Browser. The abundance of associated alternative splicing showed a general trend consistent with entropy scores, for example, higher abundance among exons with 3′SS containing the G pentamer (Figure 2B, p = 9E-07). However, further comparison among three groups of 3′SS (Figure 1A, gtggaa, ggggg and gaaaa groups) with similarly low entropy scores (around 0, p > 0.05 by t-test) indicated that not all of them were associated with higher abundance of alternative splicing. For instance, the gaaaa group had only about half of the G pentamer group in the abundance. Thus, specific purine sequences are likely also associated with the high abundance of alternative splicing.
We then chose to focus on the REPA G tracts. We examined in total 922 REPA G3–8-harboring 3′SS and found 627 of them (68%) associated with alternative splicing events (Additional file 1: Table S1). The first G mainly starts between the -11 and -6 positions (Figure 2C), and about 99% of them are G3–5(Figure 2D).
Major types of the associated alternative splicing events include alternative 3′SS or 5′SS usage, cassette exon and intron retention (Figure 3A). Of these, intron retention has the most enrichment, about 17 times of that surveyed in the human transcriptomes of different tissues . The other enriched type is alternative 3′ splice (1.6 times). Examples of these types of alternative splicing and their G tracts are shown in Figure 3B: the alternative 3′ splicing of the exon 5 of PAX8 (paired box 8), alternative 5′ splicing of the exon 22 of MYH13 (myosin, heavy chain 13, skeletal muscle), skipped exon (cassette) 2 of ABCC11 (ATP-binding cassette, sub-family C (CFTR/MRP), member 11) and retained intron 1 of TRERF1 (transcriptional regulating factor 1).
The REPA G-tracts are splicing silencers
To determine the role of the G-tracts in alternative splicing, we examined four REPA G tract-harboring 3′SS by transferring them to the upstream of a constitutive exon derived from the human beta-globin gene (Figure 4A, upper left). Three of them almost abolished the usage of the immediate downstream (proximal) 3′AG and activated substantial usage of a further downstream cryptic 3′SS within the exon or caused exon skipping (upper right, lanes 3, 5 and 7). The other one reduced the proximal 3′AG usage by about 22.4% (lane 9). Importantly, mutating the Gs to As of the REPA G tracts almost completely (lanes 4 and 10) or partially (lanes 6 and 8) restored the usage of the proximal 3′AG, accompanied by abolished (lane 4), reduced (lane 6) or increased (lane 8) cryptic splicing. This suggests that the G tracts are more than just a purine-rich spacer between the Py and 3′AG, which also weakens the 3′ splice site (see the MAXENT score and also Figure 2A). The G tracts inhibit the usage of proximal 3′AG but often activate distal 3′AG, consistent with their increased association with alternative 3′SS (Figure 3A) and the dependence of the effect on their location relative to the target splice site. Moreover, the ΔGs of the potential secondary structures caused by the mutations, as tested by Mfold , do not correlate with the splicing changes, suggesting that the mutation effect is not due to a secondary structural change. Therefore, we conclude that all four REPA G tracts tested are splicing silencers with the Gs playing critical roles.In order to examine the role of known G tract-binding factors hnRNP H/F in splicing through the REPA G tracts, we carried out RNA interference assays (Figure 4A, lower panel). Two of the reporters PAX8 and PAK4 almost completely lost the G tract silencer effect upon knockdown of hnRNP H/F (lanes 3 & 9) whereas the other two exhibited about 29.3% and 13.3% decrease in exon skipping (lanes 5 & 7). This is accompanied by abolished (lane 3), reduced (lane 5) or increased (lane 7) level of the cryptic splice product. Notably, the mutated mini-genes showed almost no change in exon skipping (lanes 4, 6 & 10) except one that showed about 12.7% decrease (lane 8) but still 10.4% less than its wild type. The mutants also showed reduced cryptic splicing in the cases of PRMT5 and GRK5 (lanes 6 and 8), suggesting that usage of the cryptic splice sites are also affected by other hnRNP H/F target element(s) in the reporter. Hence, the silencer effect of REPA G tract on exon skipping or its activation of cryptic 3′SS usage is mainly through hnRNP H/F.To further validate the G tract silencer effect on splicing of endogenous exons, we used RNA interference of hnRNP H/F in human HeLa cells and measured the effect of loss-of-function of these factors on exon usage (Figure 4B). We tested 8 exons of the above 4 and 4 additional genes. We were able to successfully detect expression and alternative splicing of 6 exons in this cell line (Figure 4B). All of them showed an increased exon inclusion (SLC44A2) or a reduced exon skipping (of the minor isoform) upon hnRNP H/F knockdown (p < 0.001). Thus, the G-tract-binding factors hnRNP H/F are essential repressors of all of the six exons of endogenous genes tested, consistent with the silencing effect of the G tracts in splicing reporters (Figure 4A).
The REPA G tracts emerged mostly in the ancestors of mammals during evolution
To determine the evolutionary origin of the G tracts, we analyzed 22 3′SS in multiple species ranging from fish to humans. We were able to determine the presence or absence of 15 of them in every class of animals analyzed (Table 1). G tracts of 10 of these 3′SS (~67%) were also found early as in marsupials, four of them (~27%) in other mammals and only two of them in birds, reptiles and fish. Therefore, it is likely that most of the G tracts (~87%) emerged in ancestors of mammals.As an example, the mammal-specific emergence of the G pentamer in CaMK1 is shown in Figure 5. Besides the full length CaMK1, there is a shorter transcript (GenBank accession # AB209395) with the retention of introns 8 and 9 within the 3′ untranslated region (UTR, Figure 5). The G pentamer GGGGG is at the upstream 3′SS of exon 9, where it might weaken the 3′SS to promote intron retention.
Alignment of the CaMK1 3′SS from 22 species indicates that the G pentamer (or tetramer in two species) is present in mammals from human to opossum but not in wallaby and platypus, nor is it present in chicken, lizard, frog or fish. Interestingly, there are three Gs in wallaby and two Gs in platypus, both interrupted by a single T. It thus seems that the G tracts emerged in mammalian ancestors but were lost in some marsupials.
The mostly mammalian emergence of the REPA G tracts is consistent with the conclusion that G triplets are more abundant in mammalian than in fish introns as previously reported in a genome-wide study .
Taken together, the examined REPA G tracts are mostly mammalian gene-specific, suggesting that they may have evolved to inhibit the splicing of a group of mammalian exons.
The REPA G tract-host genes are significantly associated with cancer
To determine whether these G tract-harboring 3′SS are associated with any specific biological functions or diseases, we analyzed 578 mapped genes of the 627 3′SS in the dataset in the Ingenuity Pathway Analyses. In this search, we also included 6 sets of randomly chosen human genes from the ENSEMBL genome database to control the specificity of the analysis. The result indicates that the 578 genes are markedly enriched for genes involved in cancer (444 genes in total, Figure 6, and see also Additional file 1: Tables S2-S3 for a complete list), by about 4 fold -log(p-value) over the 6 sets of randomly chosen genes. Particularly melanoma, carcinoma and solid tumors have the most significant enrichment of these genes, including ABCC11 (melanoma drug resistance), SNRNP200 (retinitis pigmentosa), TRERF1 (breast cancer) and WTX (Wilms tumor gene on X chromosome, Wilms tumor) [29–33]. SNRNP200 together with another seven genes encode splicing factors: HNRNPH1, HNRNPH3, KHSRP, MBNL1, PPAN, PRPF8, SNRPN. Particularly interesting is the hnRNP H1 and H3, which are G tract-binding proteins themselves, predicting an auto-regulatory circuit through alternative splicing. Moreover, most of the genes in Figures 3, 4 and 5 (except PCED1A) are involved in cancer as well (see also Additional file 1: Table S3). Furthermore, similar Ingenuity pathway analysis using the full set of the G3–8 3′SS-containing genes reached the same conclusion. Therefore, the REPA G tracts are significantly enriched in genes involved in cancer, likely to control their alternative splicing and cancer properties.
As an example of the genes involved in cancer, the WTX (or called AMER1, APC membrane recruitment protein 1) has a retained intron with a 3′SS G3 in its second exon  (Figure 6C). The full length transcript encodes an 1135aa protein that interacts with the transcription factor WT1 (Wilms tumor 1) and efficiently enhances its transcriptional activity [30, 35]. The intron splicing results in a shorter protein of 804aa lacking the WT1-interacting domain. The shorter isoform has reduced effect on WT1 transcriptional activity [30, 35]. Therefore, the 3′ REPA G tract-associated intron retention/splicing determine the high or low activity of the tumor suppressor WTX. Presence of the G tract appears to favor the full activity of the protein, likely by weakening the 3′ splice site and thus promoting intron retention.
Also enriched are biological functions associated with organismal survival (about 2.5 fold –log(p-value)), and to some extent cellular growth & proliferation (Figure 6A-B, and Additional file 1: Table S2).
Taken together, more than six hundred 3′SS that harbour REPA G tracts associated with alternative splicing are identified in the human genome. Instead of being enhancers as the other mammalian intronic G tracts, these specially localized G tracts appear to be intronic splicing silencers and mostly emerged in the ancestors of mammals. Moreover, their host genes are highly associated with cancer. These data suggest that mammals have evolved de novo RNA elements within this highly constrained space of 3′ splice site to alter the mRNA processing step for diverse proteins in particular cellular functions.
Regulatory RNA elements between Py and 3′AG (REPA)
The splicing regulatory elements can be outside or overlap with the splice site motifs. However, since the location and sequences of the polypyrimidine tract and 3′ AG have been more or less constrained for many splice sites, emergence of de novo regulatory elements between them has not been observed frequently. A G-rich sequence was identified in a group of drosophila genes but it was enriched further upstream (60nt) in the intron to inhibit splicing . While the previously reported CaRRE1 of the STREX exon support the existence of REPA elements , here the identification of a large group of 3′SS G tracts clearly demonstrate that the emergence of de novo regulatory elements do occur at this location, apparently for the regulation of alternative splicing mainly in mammals. Moreover, these elements are location-dependent silencers of the immediate downstream 3′AG (Figure 4A), instead of enhancers as found for the other mammalian intronic G tracts outside of the region between Py and 3′AG [21, 22, 24, 36]. Beyond G tracts, there are also other purine-rich elements identified in this search that await further detailed analysis to reveal their function in splicing regulation. Furthermore, this also implies the existence of other non-purine-rich elements. In fact, we have also observed UGCAUG elements at this location of alternative exons (Xie, unpublished observation). Therefore, the repertoires of REPA will apparently increase in the future. Overall, the REPA provides a distinct example for the variations of splice sites to cause alternative usage of exons .
Variations of splice sites for alternative splicing and proteomic diversity during evolution
Splicing regulatory RNA elements have evolved from fish to humans with dramatic emergence of some additional elements including G tracts . Here the REPA G tracts appear to have evolved mostly in mammalian ancestors, evident from the alignment of 3′SS of multiple species (Figure 5 and Table 1). The evolved elements likely increase the possibility for alternative splicing and cause functional changes of the proteins (Figure 6). Moreover, their recruitment of trans-acting factors, such as hnRNP H/F (Figure 4), would make it possible for upstream signaling to reach the target splice site , adding a further layer to the splicing regulation. With the help of more REPAs in the genome (Figures 1 and 2), more diverse protein isoforms and functions of mammalian genes could be generated.
Role of G tracts and alternative splicing in cancer
Accumulating evidence indicate that alternative splicing has been exploited by cancer cells to favor their survival [39, 40]. G tracts have been found to control the alternative splicing of genes involved in cancer, such as Bcl-x [22, 41]. Their trans-acting factors have also been found to regulate alternative splicing during tumorigenesis. For instance, upregulation of hnRNP H in gliomas drives the IG20 to anti-apoptotic MADD (MAP-kinase activating death domain protein) splicing in favor of invasiveness . In HeLa cells, hnRNP H/F upregulates the pro-apoptosis variant Bcl-xS . Multiple G-tracts can also form a stable four strand structure called a G quadruple that is less represented in exon sequences . They are formed in a cell-cycle-dependent way with the highest abundance in the S phase , and have been proposed as therapeutic targets for diseases including cancer, for example, by controlling oncogene transcription [45, 46]. They are also proposed to control alternative splicing of the tumor suppressor p53 , and the telomere RNA hTERT . Multiple copies of G tracts are also seen in some of the 3′SS examined here. For example the 3′SS of exon 4 of the matrix-remodeling associated 8 (MXRA8) gene has five copies of G2–5 within 50 nt of the intron end. Whether G quadruples may form and play a role in the control of alternative splicing of some of these cancer-related genes will need future studies.
In summary, the identification of hundreds of human 3′SS containing the REPA G tracts and the confirmed splicing silencer effect of representative ones suggest that the highly constrained 3′SS region has been permissive for the emergence of de novo regulatory elements for splicing control of a large group of mammalian genes. Their significant enrichment in genes involved in cancer implies that their regulation of splicing perhaps control certain properties of cancer cells.
A large group of human 3′SS have evolved purine-rich elements between the Py and 3′AG. Particularly the G tract elements are splicing silencers requiring hnRNP H/F for their activity. Such G tract-harboring human genes are significantly enriched for functions related to cancer. These data demonstrate that the often highly constrained space between the Py and 3′AG has been remarkably permissive for the evolution of splicing regulatory elements, likely to diversify protein functions in mammals. The presence of such elements at the 3′SS is also predictive of alternative splice sites.
Splicing reporters with 3′SS REPA or human G-tract mutant were created by replacing the corresponding upstream 3′SS of the middle constitutive exon of DUP175 as carried out previously .
Cell culture and transfection
HEK293T cells were grown in Dulbecco’s modified Eagle’s medium containing 10% new-born calf serum and 1% penicillin/streptomycin/glutamine solution (Invitrogen). Transfection of overnight cultures of HEK293T cells was performed with Lipofectamine 2000 (Invitrogen) according to manufacturer’s protocol, using 0.25 μg of splicing reporter minigene plasmids in 12-well plates. Transfection of siRNA into HeLa cells was done with Lipofectamine RNAiMAX (Invitrogen) according to manufacturer’s instructions.
To knock down hnRNP H/F in HeLa cells, we synthesized siRNA targeting the 20 nucleotide sequence GCGACCGAGAACGACAUUU . HeLa cells were transfected with 360 pmoles of siRNA twice in 6-well plates with 24-hour intervals. Cells were harvested after 72 hours of first transfection for protein and RNA analysis.
We performed semi-quantitative RT-PCR of endogenous genes based on previously described procedure . Specifically, in 10 ul reverse transcription reactions we used 250 ng of cytoplasmic RNA and used 0.75 ul of it in a 12.5 ul PCR reaction for 30 cycles (24 for mini-genes) with annealing temperature of 60°C. We quantified bands intensities using ImageJ (National Institutes of Health).
Human genome search
We wrote a bioperl script (Additional file 2: REPA_Perl) to search the 3′SS in the annotated EMSEMBL human genome for purine rich motifs within -10 ~ -3 position at intron ends. The consensus elements were identified using the -15 ~ -3 nucleotides by the MEME program as described previously . The MaxEntScan entropy scores of 3′SS were determined as described , using 23 nt (20 nt intron and 3 nt exon) of the 3′SS at the Burge lab website (http://genes.mit.edu/burgelab/maxent/Xmaxentscan_scoreseq_acc.html). Human genome sequence was downloaded from the ENSEMBL website (http://uswest.ensembl.org/index.html). HGNC symbols of genes containing identified 3′ spice sites were retrieved from Biomart (http://www.biomart.org/index.html) using ENSEMBL gene IDs. From these 3′SS, we determined the alternative exons by alignment with that in the UCSC Genome Browser (http://www.genome.ucsc.edu/index.html); those without evidence for the presence of splice sites were excluded from further analysis. We used ENSEMBL gene IDs to obtain the clusters of identified genes involved in same functions/diseases using Ingenuity Pathway Analysis application (http://www.ingenuity.com).
Entropy scores were compared with two-tailed Student’s t-test in Figure 2A. Ingenuity pathway analysis was by right-tailed Fisher’s exact test in Figure 6B. The abundance of purine nucleotides or alternative splicing was compared by hypergeometric test in Figures 1B and 2B, respectively.
No ethics approval was required for the research conducted for this manuscript.
Nilsen TW, Graveley BR: Expansion of the eukaryotic proteome by alternative splicing. Nature. 2010, 463 (7280): 457-463. 10.1038/nature08909.
Barbosa-Morais NL, Irimia M, Pan Q, Xiong HY, Gueroussov S, Lee LJ, Slobodeniuc V, Kutter C, Watt S, Colak R, Kim T, Misquitta-Ali CM, Wilson MD, Kim PM, Odom DT, Frey BJ, Blencowe BJ: The evolutionary landscape of alternative splicing in vertebrate species. Science. 2012, 338 (6114): 1587-1593. 10.1126/science.1230612.
Merkin J, Russell C, Chen P, Burge CB: Evolutionary dynamics of gene and isoform regulation in Mammalian tissues. Science. 2012, 338 (6114): 1593-1599. 10.1126/science.1228186.
Pan Q, Shai O, Lee LJ, Frey BJ, Blencowe BJ: Deep surveying of alternative splicing complexity in the human transcriptome by high-throughput sequencing. Nat Genet. 2008, 40 (12): 1413-1415. 10.1038/ng.259.
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.
Kelemen O, Convertini P, Zhang Z, Wen Y, Shen M, Falaleeva M, Stamm S: Function of alternative splicing. Gene. 2013, 514 (1): 1-30. 10.1016/j.gene.2012.07.083.
Feng D, Xie J: Aberrant Splicing in Neurological Diseases. Wiley Interdisciplinary Reviews RNA. 2013
Singh RK, Cooper TA: Pre-mRNA splicing in disease and therapeutics. Trends Mol Med. 2012, 18 (8): 472-482. 10.1016/j.molmed.2012.06.006.
Black DL: Mechanisms of alternative pre-messenger RNA splicing. Annu Rev Biochem. 2003, 72: 291-336. 10.1146/annurev.biochem.72.121801.161720.
Chen M, Manley JL: Mechanisms of alternative splicing regulation: insights from molecular and genomics approaches. Nat Rev Mol Cell Biol. 2009, 10 (11): 741-754.
Burge CB, Padgett RA, Sharp PA: Evolutionary fates and origins of U12-type introns. Mol Cell. 1998, 2 (6): 773-785. 10.1016/S1097-2765(00)80292-0.
Zhang MQ: Statistical features of human exons and their flanking regions. Hum Mol Genet. 1998, 7 (5): 919-932. 10.1093/hmg/7.5.919.
Moore MJ: Intron recognition comes of AGe. Nat Struct Biol. 2000, 7 (1): 14-16. 10.1038/71207.
Merendino L, Guth S, Bilbao D, Martinez C, Valcarcel J: Inhibition of msl-2 splicing by Sex-lethal reveals interaction between U2AF35 and the 3′ splice site AG. Nature. 1999, 402 (6763): 838-841. 10.1038/45602.
Xie J, Black DL: A CaMK IV responsive RNA element mediates depolarization-induced alternative splicing of ion channels. Nature. 2001, 410 (6831): 936-939. 10.1038/35073593.
Xie J, Jan C, Stoilov P, Park J, Black DL: A consensus CaMK IV-responsive RNA sequence mediates regulation of alternative exons in neurons. RNA. 2005, 11 (12): 1825-1834. 10.1261/rna.2171205.
Yu J, Hai Y, Liu G, Fang T, Kung SK, Xie J: The heterogeneous nuclear ribonucleoprotein L is an essential component in the Ca2+/calmodulin-dependent protein kinase IV-regulated alternative splicing through cytidine-adenosine repeats. J Biol Chem. 2009, 284 (3): 1505-1513.
Liu G, Razanau A, Hai Y, Yu J, Sohail M, Lobo VG, Chu J, Kung SK, Xie J: A conserved serine of heterogeneous nuclear ribonucleoprotein L (hnRNP L) mediates depolarization-regulated alternative splicing of potassium channels. J Biol Chem. 2012, 287 (27): 22709-22716. 10.1074/jbc.M112.357343.
Razanau A, Xie J: Emerging mechanisms and consequences of calcium regulation of alternative splicing in neurons and endocrine cells. Cell Mol Life Sci. 2013, 70 (23): 4527-4536. 10.1007/s00018-013-1390-5.
Liu G, Lei L, Yu J, Kung S, Xie J: Refinement of the spectra of exon usage by combined effects of extracellular stimulus and intracellular factors. Biochim Biophys Acta. 2014, 1839 (7): 537-545. 10.1016/j.bbagrm.2014.05.002.
McCullough AJ, Berget SM: G triplets located throughout a class of small vertebrate introns enforce intron borders and regulate splice site selection. Mol Cell Biol. 1997, 17 (8): 4562-4571.
Hai Y, Cao W, Liu G, Hong SP, Elela SA, Klinck R, Chu J, Xie J: A G-tract element in apoptotic agents-induced alternative splicing. Nucleic Acids Res. 2008, 36 (10): 3320-3331. 10.1093/nar/gkn207.
Wang Z, Rolish ME, Yeo G, Tung V, Mawson M, Burge CB: Systematic identification and analysis of exonic splicing silencers. Cell. 2004, 119 (6): 831-845. 10.1016/j.cell.2004.11.010.
Taliaferro JM, Alvarez N, Green RE, Blanchette M, Rio DC: Evolution of a tissue-specific splicing network. Genes Dev. 2011, 25 (6): 608-620. 10.1101/gad.2009011.
Li H, Liu G, Yu J, Cao W, Lobo VG, Xie J: In vivo selection of kinase-responsive RNA elements controlling alternative splicing. J Biol Chem. 2009, 284 (24): 16191-16201. 10.1074/jbc.M900393200.
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-394.
Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31 (13): 3406-3415. 10.1093/nar/gkg595.
Yeo G, Hoon S, Venkatesh B, Burge CB: Variation in sequence and organization of splicing regulatory elements in vertebrate genes. Proc Natl Acad Sci U S A. 2004, 101 (44): 15700-15705. 10.1073/pnas.0404901101.
Chen KG, Valencia JC, Gillet JP, Hearing VJ, Gottesman MM: Involvement of ABC transporters in melanogenesis and the development of multidrug resistance of melanoma. Pigment Cell Melanoma Res. 2009, 22 (6): 740-749. 10.1111/j.1755-148X.2009.00630.x.
Rivera MN, Kim WJ, Wells J, Driscoll DR, Brannigan BW, Han M, Kim JC, Feinberg AP, Gerald WL, Vargas SO, Chin L, Iafrate AJ, Bell DW, Haber DA: An X chromosome gene, WTX, is commonly inactivated in Wilms tumor. Science. 2007, 315 (5812): 642-645. 10.1126/science.1137509.
Gizard F, Robillard R, Barbier O, Quatannens B, Faucompre A, Revillion F, Peyrat JP, Staels B, Hum DW: TReP-132 controls cell proliferation by regulating the expression of the cyclin-dependent kinase inhibitors p21WAF1/Cip1 and p27Kip1. Mol Cell Biol. 2005, 25 (11): 4335-4348. 10.1128/MCB.25.11.4335-4348.2005.
Gizard F, Robillard R, Gross B, Barbier O, Revillion F, Peyrat JP, Torpier G, Hum DW, Staels B: TReP-132 is a novel progesterone receptor coactivator required for the inhibition of breast cancer cell growth and enhancement of differentiation by progesterone. Mol Cell Biol. 2006, 26 (20): 7632-7644. 10.1128/MCB.00326-06.
Zhang X, Lai TY, Chiang SW, Tam PO, Liu DT, Chan CK, Pang CP, Zhao C, Chen LJ: Contribution of SNRNP200 sequence variations to retinitis pigmentosa. Eye (Lond). 2013, 27 (10): 1204-1213. 10.1038/eye.2013.137.
Ota T, Suzuki Y, Nishikawa T, Otsuki T, Sugiyama T, Irie R, Wakamatsu A, Hayashi K, Sato H, Nagai K, Kimura K, Makita H, Sekine M, Obayashi M, Nishi T, Shibahara T, Tanaka T, Ishii S, Yamamoto J, Saito K, Kawai Y, Isono Y, Nakamura Y, Nagahari K, Murakami K, Yasuda T, Iwayanagi T, Wagatsuma M, Shiratori A, Sudo H: Complete sequencing and characterization of 21,243 full-length human cDNAs. Nat Genet. 2004, 36 (1): 40-45. 10.1038/ng1285.
Rivera MN, Kim WJ, Wells J, Stone A, Burger A, Coffman EJ, Zhang J, Haber DA: The tumor suppressor WTX shuttles to the nucleus and modulates WT1 activity. Proc Natl Acad Sci U S A. 2009, 106 (20): 8338-8343. 10.1073/pnas.0811349106.
Wang E, Dimova N, Cambi F: PLP/DM20 ratio is regulated by hnRNPH and F and a novel G-rich enhancer in oligodendrocytes. Nucleic Acids Res. 2007, 35 (12): 4164-4178. 10.1093/nar/gkm387.
Gilbert W: Why genes in pieces?. Nature. 1978, 271 (5645): 501-10.1038/271501a0.
Xie J: Differential evolution of signal-responsive RNA elements and upstream factors that control alternative splicing. Cell Mol Life Sci. 2014, 15: 1143-
David CJ, Manley JL: Alternative pre-mRNA splicing regulation in cancer: pathways and programs unhinged. Genes Dev. 2010, 24 (21): 2343-2364. 10.1101/gad.1973010.
Shkreta L, Bell B, Revil T, Venables JP, Prinos P, Elela SA, Chabot B: Cancer-associated perturbations in alternative Pre-messenger RNA splicing. Cancer Treat Res. 2013, 158: 41-94. 10.1007/978-3-642-31659-3_3.
Garneau D, Revil T, Fisette JF, Chabot B: Heterogeneous nuclear ribonucleoprotein F/H proteins modulate the alternative splicing of the apoptotic mediator Bcl-x. J Biol Chem. 2005, 280 (24): 22641-22650. 10.1074/jbc.M501070200.
Lefave CV, Squatrito M, Vorlova S, Rocco GL, Brennan CW, Holland EC, Pan YX, Cartegni L: Splicing factor hnRNPH drives an oncogenic splicing switch in gliomas. EMBO J. 2011, 30 (19): 4084-4097. 10.1038/emboj.2011.259.
Huppert JL, Balasubramanian S: Prevalence of quadruplexes in the human genome. Nucleic Acids Res. 2005, 33 (9): 2908-2916. 10.1093/nar/gki609.
Biffi G, Tannahill D, McCafferty J, Balasubramanian S: Quantitative visualization of DNA G-quadruplex structures in human cells. Nat Chem. 2013, 5 (3): 182-186. 10.1038/nchem.1548.
Balasubramanian S, Hurley LH, Neidle S: Targeting G-quadruplexes in gene promoters: a novel anticancer strategy?. Nat Rev Drug Discov. 2011, 10 (4): 261-275. 10.1038/nrd3428.
Balasubramanian S, Neidle S: G-quadruplex nucleic acids as therapeutic targets. Curr Opin Chem Biol. 2009, 13 (3): 345-353. 10.1016/j.cbpa.2009.04.637.
Marcel V, Tran PL, Sagne C, Martel-Planche G, Vaslin L, Teulade-Fichou MP, Hall J, Mergny JL, Hainaut P, Van Dyck E: G-quadruplex structures in TP53 intron 3: role in alternative splicing and in production of p53 mRNA isoforms. Carcinogenesis. 2011, 32 (3): 271-278. 10.1093/carcin/bgq253.
Gomez D, Lemarteleur T, Lacroix L, Mailliet P, Mergny JL, Riou JF: Telomerase downregulation induced by the G-quadruplex ligand 12459 in A549 cells is mediated by hTERT RNA alternative splicing. Nucleic Acids Res. 2004, 32 (1): 371-379. 10.1093/nar/gkh181.
Cao W, Razanau A, Feng D, Lobo VG, Xie J: Control of alternative splicing by forskolin through hnRNP K during neuronal differentiation. Nucleic Acids Res. 2012, 40 (16): 8059-8071. 10.1093/nar/gks504.
This work is supported by a discovery grant from the Natural Science and Engineering Research Council of Canada (NSERC) to JX. JX is also a recipient of Manitoba Research Chair fund from the Manitoba Health Research Council (MHRC).
The authors declare that they have no competing interests.
JX conceptualized the project. MS, WC and NM carried out the experiments. SPH wrote the bioperl script. All including MM contributed to data analysis. JX, MS, MM and WC prepared the manuscript. All authors read and approved the contents of the final manuscript.
Electronic supplementary material
Authors’ original submitted files for images
About this article
Cite this article
Sohail, M., Cao, W., Mahmood, N. et al. Evolutionarily emerged G tracts between the polypyrimidine tract and 3′ AG are splicing silencers enriched in genes involved in cancer. BMC Genomics 15, 1143 (2014) doi:10.1186/1471-2164-15-1143
- 3′ splice site
- Alternative splicing