Abstract
Alternative splicing is a post-transcriptional regulatory process that is attaining stronger recognition as a modulator of gene expression. Alternative splicing occurs when the primary RNA transcript is differentially processed into more than one mature RNAs. This is the result of a variable definition/inclusion of the exons, the sequences that are excised from the primary RNA to form the mature RNAs. Consequently, RNA expression can generate a collection of differentially spliced RNAs, which may distinctly influence subsequent biological events, such as protein synthesis or other biomolecular interactions. Still the mechanisms that control exon definition and exon inclusion are not fully clarified. This mini-review highlights advances in this field as well as the impact of single nucleotide polymorphisms in affecting splicing decisions. The Glioma-associated oncogene 1, GLI1, is taken as an example in addressing the role of nucleotide substitutions for splicing regulation.
REGULATION OF SPLICING – Cis ELEMENTS, TRANS FACTORS
Alternative splicing has, in recent years, been demonstrated to occur in almost all intron-containing human genes (; ). However, the regulatory mechanisms that control splice choices are not fully deciphered. Apart from the general problem of defining exonic sequences in a vast excess of intronic DNA (), a key issue concerns the mechanisms of the differential inclusion of exons in mature mRNAs. The prevailing views favor that this is controlled by the relative abundance of a limited number of splicing factors in a particular cell type (). Additionally, alternative promoter usage may affect splicing decisions, as this can influence the recruitment of splicing factors or the elongation rate of RNA polymerase II (; ). Furthermore, histone modifications are also thought to have a role in modulating splice choices (). Moreover, the RNA secondary structure may affect splicing decisions (), with this claim being suggested more than 25 years ago (). In fact, the expansion of splicing factors, including the SR and hnRNP proteins, in higher eukaryotes () is in agreement with the increased alternative splicing events in more complex organisms ().
In addition to this trans factor dependent regulation of splice choices, cis sequences are also important in influencing splicing decisions, as these are thought to act as binding sites for the splicing factors. To start with, exons are flanked by the invariant intronic dinucleotides AG and GT, with the introns also characterized by a polypyrimidine track near their 3′ end and the branch point adenosine residue used for lariat formation during intron excision. However, it is obvious that additional cis sequence information is required in defining exons and influencing their inclusion in mRNAs. Such auxiliary sequences, which stimulate splicing are found in both exons (exonic splicing enhancers, ESE), and introns (intronic splicing enhancers, ISE; ; ; ). Moreover, sequences that inhibit splicing have also been characterized, as exonic and intronic splicing silencers (ESS and ISS; ; ; ). Additionally, the sequence complexity of these splicing regulatory elements is highlighted by the in vivo analysis of all possible hexamers as exonic splicing regulators, which revealed that out of the 4,096 combinations more than half could act as either ESE or ESS ().
SINGLE NUCLEOTIDE POLYMORPHISMS IN HUMAN GENOMES AND ALTERNATIVE SPLICING
The advances in sequencing technologies have made possible the complete sequence of individual human genomes at a rapid pace (). Consequent is the expansion of the single nucleotide polymorphisms (SNP) entries in the dbSNP database, which currently reaches more than 40,000,000 validated RefSNPs1. Thus, on the average, a SNP is present in every 75 nucleotides within the three billion human genome sequence. Moreover, since the human protein-coding gene content is estimated to be 2% of the human genome, 800,000 SNPs can potentially influence gene expression and splicing, not to mention variations outside protein-coding genes that could alter gene regulation by mechanisms that are still poorly understood ().
What characterizes the SNPs that modulate splicing? Firstly, exons that are known to be alternatively spliced, in contrast to constitutive exons, appear to have a higher proportion of SNPs that affect exonic splicing regulators. Additionally, ESSs are more frequently modified by SNPs than ESEs (). A systematic analysis of the SNPs that are known to affect splicing revealed a number of characteristics, including a predominance of ESS gains relative to ESE losses (). These features have been introduced in a web-based tool allowing the prediction of the impact of SNPs on alternative splicing regulation, the Skippy software2. In this context it is worth noting that SNPs that affect the strength the 5′ splice site, but not the conserved GT dinucleotide, are compensated by other regulatory signals, resulting in minimal changes in splicing, arguing in favor of a certain robustness in maintaining splicing patterns ().
MUTATIONS AND GENE FUNCTION
It has generally been assumed that exonic mutations (the term mutation refers to a polymorphism occurring in less than 1% of the population) affect the protein product of the corresponding gene by introducing a stop codon (non-sense mutation), by altering an amino acid (missense mutation) or by changing the translation reading frame (frameshift mutation). However, this postulate had been challenged by the finding that nucleotide substitutions that do not change the encoded amino acid (synonymous mutation) may also alter protein structure, through an indirect way, namely by modulating splicing patterns (). This can be rationalized by either an activation of a cryptic splice site, as in the case of the LMNA gene (), or the presence of ESE and ESS motifs, which are affected by the introduction of a synonymous nucleotide substitution and consequently may alter the canonical splicing pattern. In fact, it has been convincingly argued that the genetic code can effectively carry additional information, such as splicing signals, besides its capacity to encode amino acids (). Moreover, an analysis of 27,681 non-sense and missense disease-causing mutations revealed that about 25% of these actually change exonic splicing regulators, providing evidence that alterations of splicing patterns is not a sole characteristic of synonymous substitutions ().
THE GLIOMA-ASSOCIATED ONCOGENE 1, GLI1 – ALTERNATIVE SPLICING
The GLI1 oncogene encompasses about 12 kb of genomic sequences in chromosome 12q13.2–q13.3 and is composed of two 5′ non-coding exons, exons 1 and 1A, and 11 coding exons, exons 2 to 12. Two splice variants have been identified in human tissues that affect coding exons, GLI1ΔN, which lacks exons 2 and 3 (), and tGLI1, which lacks exon 3 and part of exon 4 (). However, variants that do or do not include the non-coding exon 1A have also been described (; Figure 1A).
FIGURE 1
In the mouse, the Gli1 gene contains a third non-coding exon, exon 1B (Figure 1B). Interestingly, the 5′ splice site of exon 1B in the BALB/c mouse strain contains the canonical GT dinucleotide, while in the C57BL/6 mouse strain this is substituted with GC. The splice signal alteration confers a preferential inclusion of exon 1B in BALB/c mice resulting in a longer 5′ UTR region. Moreover, this increased 5′ UTR length was found to have a negative impact on the translation efficiency of the Gli1 mRNA, with a consequent reduced amount of synthesized protein (). These observations argue that SNPs that affect splicing may not necessarily alter the primary structure of the protein product but can also have an impact on the amount of produced protein, as a result of changes in the length/structure of the 5′ UTR region.
THE GLIOMA-ASSOCIATED ONCOGENE 1, GLI1 – SNPs
The dbSNP database revealed that there are 249 variants in the human GLI1 gene sequence (May 2012). From these 127 are intronic, 32 are synonymous, 59 are non-synonymous, and 14 are in the UTR regions.
By centering on the genomic region that engages in alternative splicing of protein-coding exons, namely exons 2 to 4 and the 100 bases flanking sequence, 27 SNPs, with four being synonymous and four non-synonymous, were identified. Moreover, the use of the Skippy software revealed that the highest score for exon skipping of these eight exonic SNPs is for the C to T missense rs143548857 in exon 4. Direct experimentation will be required to examine whether this predicted impact on exon 4 splicing is indeed taking place, as this will introduce a translational frameshift.
From the 36 synonymous and non-synonymous SNPs in the remaining protein-coding internals exons, namely exons 5 to 11, only the C to T synonymous rs7973381 in exon 8 has an equally high score for exon exclusion as rs143548857. Interestingly, the size of exon 8 is 150 nucleotides, highlighting that potential exon skipping will not initiate non-sense mediated decay ().
Since in the natural tGLI1 variant an alternative 3′ splice site is used in exon 4 (Figure 1A), SNPs that could potentially mediate activation of ectopic splice sites in exon 4 were examined using the Skippy software. Interestingly, the G to A synonymous rs139988338 in exon 4 creates a 3′ ectopic splice site with a high score. Moreover, none of the SNPs in the remaining internal exons of GLI1 create a 3′ ectopic splice site that has a higher score. Worth noting is that the strength of the rs139988338-mediated 3′ ectopic splice site in exon 4 is, in fact, higher than the natural tGLI1 alternative 3′ splice site, according to the MaxEnt splice site scoring program3.
Additionally, in the non-coding but alternatively spliced exon 1A three SNPs are identified, rs10783826, rs10783827, and rs118093490. These SNPs cannot effectively be evaluated by the Skippy software, as the program is not optimized for non-coding exons (Laura Elnitski, personal communication). Moreover, exon 1A is not fully annotated in the human Genome Build 37, which automatically excludes it from the software analysis, since only annotated exons can be used. On the other hand, manual inspection of the elicited base substitution by the three SNPs and taking advantage of the published strength of all hexamers as ESEs and ESSs () revealed that for rs10783826 and rs118093490 the nucleotide change increased the impact of the corresponding hexamer as a silencer, while for rs10783827 the change was in favor of an enhancer. These findings were observed irrespective of whether the hexamers were shifted by one base so the SNP was placed at either the third or the fourth nucleotide position within the hexamer. Moreover, the SROOGLE software4 (), which uses additional programs to identify ESEs and ESSs, revealed a similar trend. Thus, it is possible that the differential inclusion of exon 1A may be dependent of these SNPs, and the resulting variable length of the 5′ UTR region could have an impact on the amount of the translated protein product, as previously seen with the mouse Gli1.
CONCLUSION
Alternative splicing is a major post-transcriptional mechanism that diversifies gene expression. In addition to trans acting factors, cis regulatory sequences have also an impact in controlling splice choices, highlighting the potential significance of SNPs in modulating splicing decisions. SNPs that affect splicing may be present in both intronic and exonic regions, mapping at regulatory elements such as splicing enhancers and silencers. Moreover, for exonic SNPs, not only synonymous but also missense and non-sense SNPs can influence splice choices, providing evidence that the major biological impact of some of these exon variants may be on mRNA splicing rather than on protein synthesis.
Using the Glioma-associated oncogene 1, GLI1, as an example, SNPs are predicted to have a role in controlling splicing patterns. Interestingly, GLI1 exons that are known to engage in alternative splicing appear to preferentially encompass SNPs that may affect splice choices. Finally, although user-friendly prediction programs simplify the analysis of the role of SNPs in splicing regulation, the requirement of direct experimentation has not been fully substituted.
Statements
Acknowledgments
This work has been supported by the Swedish Childhood Cancer Foundation, the Swedish Cancer Fund, and the Swedish Research Council.
Conflict of interest
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Footnotes
1.^http://www.ncbi.nlm.nih.gov/projects/SNP/snp_summary.cgi
2.^http://research.nhgri.nih.gov/skippy
3.^http://genes.mit.edu/burgelab/maxent/Xmaxentscan_scoreseq_acc.html
REFERENCES
1
AznarezI.BarashY.ShaiO.HeD.ZielenskiJ.TsuiL. C.ParkinsonJ.FreyB. J.RommensJ. M.BlencoweB. J. (2008). A systematic analysis of intronic sequences downstream of 5′ splice sites reveals a widespread role for U-rich motifs and TIA1/TIAL1 proteins in alternative splicing regulation.Genome Res.181247–1258.
2
BrooksA. N.AspdenJ. L.PodgornaiaA. I.RioD. C.BrennerS. E. (2011). Identification and experimental validation of splicing regulatory elements in Drosophila melanogaster reveals functionally conserved splicing enhancers in metazoans.RNA171884–1894.
3
BuschA.HertelK. J. (2012). Evolution of SR protein and hnRNP splicing regulatory factors.Wiley Interdiscip. Rev. RNA31–12.
4
ChamaryJ. V.ParmleyJ. L.HurstL. D. (2006). Hearing silence: non-neutral evolution at synonymous sites in mammals.Nat. Rev. Genet.798–108.
5
ChenM.ManleyJ. L. (2009). Mechanisms of alternative splicing regulation: insights from molecular and genomics approaches.Nat. Rev. Mol. Cell Biol.10741–754.
6
de SouzaJ. E.RamalhoR. F.GalanteP. A.MeyerDde SouzaS. J. (2011). Alternative splicing and genetic diversity: silencers are more frequently modified by SNVs associated with alternative exon/intron borders.Nucleic Acids Res.394942–4948.
7
DurandS.Lykke-AndersenJ. (2011). SnapShot: nonsense-mediated mRNA decay.Cell145324–324.e2.
8
ErikssonM.BrownW. T.GordonL. B.GlynnM. W.SingerJ.ScottL.ErdosM. R.RobbinsC. M.MosesT. Y.BerglundP.DutraA.PakE.DurkinS.CsokaA. B.BoehnkeM.GloverT. W.CollinsF. S. (2003). Recurrent de novo point mutations in lamin A cause Hutchinson–Gilford progeria syndrome.Nature423293–298.
9
ItzkovitzS.AlonU. (2007). The genetic code is nearly optimal for allowing additional information within protein-coding sequences.Genome Res.17405–412.
10
KeS.ShangS.KalachikovS. M.MorozovaI.YuL.RussoJ. J.JuJ.ChasinL. A. (2011). Quantitative evaluation of all hexamers as exonic splicing elements.Genome Res.211360–1374.
11
KimE.MagenA.AstG. (2007). Different levels of alternative splicing among eukaryotes.Nucleic Acids Res.35125–131.
12
KornblihttA. R. (2005). Promoter usage and alternative splicing.Curr. Opin. Cell Biol.17262–268.
13
LoH. W.ZhuH.CaoX.AldrichA.Ali-OsmanF. (2009). A novel splice variant of GLI1 that promotes glioblastoma cell migration and invasion.Cancer Res.696790–6798.
14
LomelinD.JorgensonE.RischN. (2010). Human genetic variation recognizes functional elements in noncoding sequence.Genome Res.20311–319.
15
LuZ. X.JiangP.CaiJ. J.XingY. (2011). Context-dependent robustness to 5′ splice site polymorphisms in human populations.Hum. Mol. Genet.201084–1096.
16
LucoR. F.PanQ.TominagaK.BlencoweB. J.Pereira-SmithO. M.MisteliT. (2010). Regulation of alternative splicing by histone modifications.Science327996–1000.
17
McManusC. J.GraveleyB. R. (2011). RNA structure and the mechanisms of alternative splicing.Curr. Opin. Genet. Dev.21373–379.
18
MuñozM. J. Pérez SantangeloM. S.ParonettoM. P.de la MataM.PelischF.BoireauS.Glover-CutterK.Ben-DovC.BlausteinM.LozanoJ. J.BirdG.BentleyD.BertrandE.KornblihttA. R. (2009). DNA damage regulates alternative splicing through inhibition of RNA polymerase II elongation.Cell137708–720.
19
PalaniswamyR.TeglundS.LauthM.ZaphiropoulosP. G.ShimokawaT. (2010). Genetic variations regulate alternative splicing in the 5′ untranslated regions of the mouse glioma-associated oncogene 1, Gli1.BMC Mol. Biol.1132 10.1186/1471-2199-11-32
20
PanQ.ShaiO.LeeL. J.FreyB. J.BlencoweB. J. (2008). Deep surveying of alternative splicing complexity in the human transcriptome by high-throughput sequencing.Nat. Genet.401413–1415.
21
SchwartzS.HallE.AstG. (2009b). SROOGLE: webserver for integrative, user-friendly visualization of splicing signals.Nucleic Acids Res.37W189–W192.
22
SchwartzS.MeshorerE.AstG. (2009a). Chromatin organization marks exon–intron structure.Nat. Struct. Mol. Biol.16990–995.
23
ShimokawaT.TostarU.LauthM.PalaniswamyR.KasperM.ToftgårdR.ZaphiropoulosP. G. (2008). Novel human glioma-associated oncogene 1 (GLI1) splice variants reveal distinct mechanisms in the terminal transduction of the hedgehog signal.J. Biol. Chem.28314345–14354.
24
SolnickD. (1985). Alternative splicing caused by RNA secondary structure.Cell43667–676.
25
Sterne-WeilerT.HowardJ.MortM.CooperD. N.SanfordJ. R. (2011). Loss of exon identity is a common mechanism of human inherited disease.Genome Res.211563–1571.
26
ViselA.RubinE. M.PennacchioL. A. (2009). Genomic views of distant-acting enhancers.Nature461199–205.
27
WangE. T.SandbergR.LuoS.KhrebtukovaI.ZhangL.MayrC.KingsmoreS. F.SchrothG. P.BurgeC. B. (2008). Alternative isoform regulation in human tissue transcriptomes.Nature456470–476.
28
WangX. Q.RothnagelJ. A. (2001). Post-transcriptional regulation of the gli1 oncogene by the expression of alternative 5′ untranslated regions.J. Biol. Chem.2761311–1316.
29
WenJ.ChibaA.CaiX. (2010). Computational identification of tissue-specific alternative splicing elements in mouse genes from RNA-Seq.Nucleic Acids Res.387895–7907.
30
WoolfeA.MullikinJ. C.ElnitskiL. (2010). Genomic features defining exonic variants that modulate splicing.Genome Biol.11R20.
31
XiaJ.WangQ.JiaP.WangB.PaoW.ZhaoZ. (2012). NGS catalog: a database of next generation sequencing studies in humans.Hum. Mutat. 10.1002/humu.22096 [Epub ahead of print].
32
YuY.MaroneyP. A.DenkerJ. A.ZhangX. H.DybkovO.LührmannR.JankowskyE.ChasinL. A.NilsenT. W. (2008). Dynamic regulation of alternative splicing by silencers that modulate 5′ splice site competition.Cell1351224–1236.
33
ZhangC.LiW. H.KrainerA. R.ZhangM. Q. (2008). RNA landscape of evolution for optimal exon and intron discrimination.Proc. Natl. Acad. Sci. U.S.A.1055797–5802.
Summary
Keywords
exon skipping, splice signal, splicing enhancer, splicing silencer, cryptic splice site
Citation
Zaphiropoulos PG (2012) Genetic variations and alternative splicing: the Glioma associated oncogene 1, GLI1. Front. Gene. 3:119. doi: 10.3389/fgene.2012.00119
Received
06 May 2012
Accepted
13 June 2012
Published
06 July 2012
Volume
3 - 2012
Edited by
Peng Jiang, University of Iowa, USA
Reviewed by
Serdar Bozdag, National Institutes of Health, USA Paola Bonizzoni, Università di Milano-Bicocca, Italy
Copyright
© Zaphiropoulos.
This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and subject to any copyright notices concerning any third-party graphics etc.
*Correspondence: Peter G. Zaphiropoulos, Department of Biosciences and Nutrition, Karolinska Institutet, 141 83 Huddinge, Sweden. e-mail: peter.zaphiropoulos@ki.se
This article was submitted to Frontiers in Bioinformatics and Computational Biology, a specialty of Frontiers in Genetics.
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.