- Research article
- Open Access
Genetic variations regulate alternative splicing in the 5' untranslated regions of the mouse glioma-associated oncogene 1, Gli1
BMC Molecular Biology volume 11, Article number: 32 (2010)
Alternative splicing is one of the key mechanisms that generate biological diversity. Even though alternative splicing also occurs in the 5' and 3' untranslated regions (UTRs) of mRNAs, the understanding of the significance and the regulation of these variations is rather limited.
We investigated 5' UTR mRNA variants of the mouse Gli1 oncogene, which is the terminal transcriptional effector of the Hedgehog (HH) signaling pathway. In addition to identifying novel transcription start sites, we demonstrated that the expression ratio of the Gli1 splice variants in the 5' UTR is regulated by the genotype of the mouse strain analyzed. The GT allele, which contains the consensus intronic dinucleotides at the 5' splice site of intron 1B, favors exon 1B inclusion, while the GC allele, having a weaker 5' splice site sequence, promotes exon 1B skipping. Moreover, the alternative Gli1 5' UTRs had an impact on translational capacity, with the shorter and the exon 1B-skipped mRNA variants being most effective.
Our findings implicate novel, genome-based mechanisms as regulators of the terminal events in the mouse HH signaling cascade.
Alternative splicing and transcriptional initiation are key mechanisms, which generate diversity both at the mRNA and protein levels. Recently, several independent research efforts revealed that more than 90% of human genes are alternatively spliced [1, 2], and about 50% of both human and mouse genes have multiple alternative promoters . Additionally, a genome-wide screening of alternative splicing and transcriptional initiation estimated that a significant number of genes are differentially spliced within 5' and 3' untranslated regions (UTRs) . Moreover, another genome-wide analysis identified 324 out of 17897 genes that display associations between flanking single nucleotide polymorphisms (SNPs) and gene expression/alternative transcription, demonstrating the regulatory effects of genetic variation in human populations . These non-bias/genome-wide analyses highlight the importance of alternative splicing/promoter usage as general mechanisms of regulation control in mammalian cells.
UTRs are considered to influence gene expression by modulating mRNA stability and/or translational efficiency. Consequently, UTR heterogeneity for a specific gene is likely to have a differential impact on protein expression . Analysis of variable 5' UTRs in the TGF-beta, BRCA1 and MDM2 genes, have indicated that the shorter UTR variants are translated more efficiently .
The Hedgehog (HH) signaling pathway plays a central role in embryonic development and adult tissue homeostasis . Abnormal activation of the pathway has been associated with various cancers in skin, brain, lung, digestive tract, prostate and pancreas [9–12]. The mechanistic details of the HH signaling pathway, which is generally thought to be well conserved in evolution, have mostly emerged from studies in Drosophila. In the absence of HH ligands, the PTCH receptor inhibits the activity of the 7-pass transmembrane protein Smoothened (SMO), which acts as a positive regulator of the pathway. Interaction of HH ligands with PTCH relieves the inhibitory control of PTCH on SMO, allowing the GLI transcription factor to dissociate from the negative regulator Suppressor of Fused (SUFU) and translocate into the nucleus, activating target genes. In mammals there are three paralogues of HH proteins, Sonic HH, Indian HH and Desert HH, two PTCH receptors, PTCH1 and PTCH2, and three GLI transcription factors, GLI1, GLI2 and GLI3. Additionally, splice variants of several HH signaling components have been identified , in line with an earlier report, which, by the use of genome-wide RNAi, highlighted the importance of alternative splicing in this pathway .
GLI1 was originally identified as a highly expressed gene in human glioma  and acts as a downstream effector of the HH signaling cascade, mediating the transcriptional response . GLI1 is also a target gene of the HH pathway, resulting in a positive feedback loop. Moreover, overexpression of GLI1 in transgenic mice leads to the induction of basal cell carcinomas and trichoepitheliomas . Alternative splicing in the GLI1 5' UTR regions of human and mouse has earlier been reported by Wang and Rothnagel . Recently, we and others have demonstrated the occurrence of functional, differential splicing events in the coding regions of human GLI1 [19, 20].
In this study, novel mouse Gli1 mRNA variants, whose transcriptional initiation is further upstream of the reported exon 1 sequence, were identified. Additionally, we obtained evidence that genetic variation is a key determining factor for the alternative splicing events in the Gli1 5' UTRs, which have functional implications on translational efficiency.
Identification of novel transcription start sites in the mouse Gli1 gene
In a previous report, we have shown the presence of novel human GLI1 splice variants, which skip exons 2 and 3 . Moreover, Wang and Rothnagel had earlier reported the skipping of exon 1A and/or exon 1B in mouse Gli1 . To clarify the splicing pattern of the mouse Gli1 gene at its 5' end, R apid A mplificaton of c DNA E nds (RACE) analysis using reverse primers within exon 4 (Table 1) was performed on RNA extracted from NIH3T3 cells. Surprisingly, we did not detect any RACE product below the expected size of a transcript that includes exon 1, 1A, 1B, 2, 3 and 4 (Figure 1A, white triangle). On the contrary, the major products observed were significantly longer. This finding implies the presence of other exon(s) and/or the extension of known exons. Sequence-verification of the RACE products demonstrated the existence of novel transcriptional start sites (TSSs) for exon 1, which are located approximately 350 (TSS-L) or 290 (TSS-M) nucleotides upstream of the reported Gli1 TSS, TSS-S (Figure 1B, DDBJ: AB234616 and AB234617). Additionally, we analyzed Gli1 exon-skipping events in the HH signaling responsive NIH3T3 cells, with or without treatment by the S mo ag onist SAG , the constitutively activated Ptch1-/- mouse embryonic fibroblast (MEF) cell line, and C57BL/6 mouse embryos at 8.5 and 9.5 days post coitum (d.p.c.) by nested PCR with primer sets within exon 1 and 4 (Table 1). Inclusion of all exons (1-1A-1B-2-3-4) and exon 1B skipping were the two splicing events detected in the analyzed samples (Figure 1C). Interestingly, the ratio of these two variants was different depending on context. Moreover, we could neither detect exon 2 and 3 skipping, as observed in human , nor skipping of both exons 1A and 1B  in these samples or in a mouse cDNA panel (data not shown). Additionally, the recently reported exon 3-partial exon 4 skipping in human , was not observed, even though the reverse primers used hybridized to the retained exon 4 segment. Note that all mouse Gli1 mRNA variants observed are characterized by unique 5' untranslated regions (5' UTRs), but a common translational start site in exon 2.
Expression profiles of Gli1 variants in embryos, cell lines, and medulloblastoma tumors
To compare the expression pattern of the 5' Gli1 variants, we performed real-time RT-PCR with PCR primers designed to specifically detect the individual Gli1 transcripts (Figure 2A and Table 2). As a control, we also used another primer set, from exon 11 to exon 12, to detect the 3' end of the Gli1 mRNAs.
Initially, we confirmed Gli1 variant expression in the C57BL/6 mouse embryos at 8.5 and 9.5 d.p.c. by real-time RT-PCR (Figure 2B). Similarly to the nested PCR analysis shown in figure 1B, Gli1 variants skipping exon 1B (ΔE1B) were more abundant to the ones including exon 1B (E1B) in both embryonic stages. Furthermore, the C ycle t hreshold (Ct) values of Gli1 variants that were transcribed from TSS-L and TSS-M (L+M) were nearly equal to the Ct values of either the ΔE1B or the exon 11 to 12 (E11-12) transcripts. This finding highlighted that the majority of the mouse Gli1 mRNAs initiate at upstream TSSs in the embryo, in agreement with the previous observations of the NIH3T3 cell line (Figure 1B). Additionally, to elucidate the in vivo distribution of these Gli1 transcripts in developing embryonic compartments, we analyzed 9.5 d.p.c. mouse embryos by whole mount in situ hybridization (See additional file 1). All probes used showed comparable distribution patterns of the variants, as also did a Gli1 3' end control riboprobe. These data suggested that the expression of all the alternative Gli1 mRNAs was likely to be controlled by similar developmental mechanisms.
To test the impact of HH signaling activation on the expression of the Gli1 variants, real-time RT-PCR analysis was performed in the HH signaling responsive NIH3T3 cells and wild type MEFs (wtMEFs), and the constitutively active Sufu-/- and Ptch1-/- MEFs (Figure 3A and 3B). In NIH3T3 cells, the E1B variant was more abundant than the ΔE1B, and HH signaling activation, elicited by SAG, could not alter the predominance of E1B. On the other hand, wtMEFs expressed the ΔE1B variant as the major Gli1 transcript and, in line with NHI3T3 cells, exhibited a parallel upregulation of the alternative Gli1 variants by HH signaling activation, resulting in an expression pattern similar to that of Sufu-/- MEFs. Moreover, Ptch1-/- MEFs equally expressed both the E1B and the ΔE1B variants, even though the Ct value of the L+M transcripts was similar to that of SAG-treated wtMEFs or Sufu-/- MEFs. Importantly, the wtMEFs and Sufu-/- MEFs, which were ΔE1B variant-dominant, were both established from C57BL/6 mice, NIH3T3 cells, which were E1B variant-dominant were from BALB/c mice, while Ptch1-/- MEFs, which had comparable expression of both variants, were from mice with mixed genetic background .
Overactivated HH signaling and Gli1 overexpression play a central role in medulloblastoma tumorigenesis . This prompted us to investigate the expression profile of the Gli1 variants in medulloblastomas from the tumor-prone Sufu+/-Trp53-/- mice maintained on a C57BL/6 genetic background (Heby-Henricson, K, Bergström, Å, Rozell, B, Toftgård, R, Teglund, S, unpublished) by using real-time RT-PCR (Figure 3C). As anticipated, the expression of all Gli1 variants was remarkably upregulated in the tumor samples compared to normal cerebellum, but the relative expression pattern of the variants was not apparently influenced by tumorigenesis. Additionally, the ΔE1B variant was found to be expressed at higher levels than E1B, as in wtMEFs and Sufu-/- MEFs. Similar upregulation patterns of the Gli1 transcripts were also observed in medulloblastomas from Ptch1+/- mice maintained on a C57BL/6 genetic background (See additional file 2).
These results are therefore suggesting that neither HH signaling activation nor tumorigenesis preferentially affect the regulation of expression of the individual Gli1 variants. On the other hand there is a strong influence of the mouse strain/genotype. C57BL/6 mice predominantly express the ΔE1B variant, BALB/c mice the E1B variant, while mice with mixed genetic background have comparable expression of both. In line with these observations, a mouse cDNA panel from BALB/c mice consistently showed increased expression of the E1B relative to the ΔE1B variant (data not shown).
Identification of SNP and SINE polymorphisms involved in mouse Gli1 exon 1B inclusion/skipping
Since our results clearly indicated that the relative expression pattern of the Gli1 variants was not affected by either HH signaling or tumorigenesis, but instead the ΔE1B to E1B ratio was well correlated with the mouse strains analyzed, variations in the genomic DNA sequence, such as single nucleotide polymorphisms (SNPs), could be involved. To identify differences among the mouse strains that influence the expression of these variants, we PCR amplified and sequenced the Gli1 exon 1B genomic region from NHI3T3, wtMEFs and Ptch1-/- MEFs. Three SNPs in the vicinity of the donor site of the intron 1B were identified, one of which affected the conserved 5' intronic dinucleotides (Figure 4, DDBJ: AB520987 - AB520989). NIH3T3 cells, which predominantly express the exon 1B-included variants, have the canonical GT dinucleotide sequence at the intronic boundary. On the other hand, wtMEFs, as well as Sufu-/- MEFs (data not shown), which predominantly express the exon 1B-skipped variants, have a GC intronic dinucleotide, implying a less efficient exon definition and consequently an increased exon skipping . Interestingly, Ptch1-/- MEFs are apparently heterozygous, containing both the GT and the GC alleles, in line with the comparable expression of the exon 1B-included and -skipped variants. Moreover, to further assess the effects of these SNPs on splicing efficiency, the GT and GC alleles were analyzed with two different splice site predictors, Human Splicing Finder  and NNSPLICE 0.9 . Both programs demonstrated that the SNPs affect the prediction of splice sites, with the GT allele predicted to contain the expected 5' splice site, while this was not the case for the GC allele.
Additionally, C3H/10T1/2 cells, which are generally used for HH signaling-dependent differentiation assays , were analyzed. Interestingly, no exon 1B inclusion could be detected by either real-time RT-PCR (Figure 5A) or nested PCR with primers within exon 1 and exon 4 (data not shown). Moreover, sequencing the exon 1B genomic region identified two insertions, a SINE B1 repeat in intron 1A and importantly, a SINE B2 repeat in exon 1B, which expanded the exonic sequence from 119 to 328 nucleotides (Figure 5B and 5C, DDBJ: AB520990). Consequently, the expanded exon 1B sequence, in combination with the presence of the GC intronic dinucleotide, further reduced exon definition, resulting in complete exon skipping . These findings are therefore providing support to the notion of the importance of sequence polymorphisms in defining splicing patterns.
Functional analysis of the Gli1 5' untranslated regions
Since the expression pattern of Gli1 variants was depended on the genetic background and was maintained as a fingerprint in different mouse strains, the possibility of regulatory roles for the alternative 5' UTRs was considered . This was further supported by the fact that the 5' UTRs of mouse Gli1 have high GC content and are predicted to form multi- and long stem-loop structures (See additional file 3). To address the functional impact of these untranslated regions on the Gli1 mRNA transcripts, six 5' UTRs, namely L, M, S, LΔ1B, MΔ1B and SΔ1B (Figure 6A), were cloned in front of the Renilla luciferase gene of the psiCHECK2 vector. These constructs were transiently transfected into wtMEFs or NIH3T3 cells, with or without signaling activation by SAG treatment, and into Ptch1-/- or Sufu-/- MEFs. In all cases, elimination of 5' segments, from L to M and then to S resulted in higher luciferase activity (Figure 6B, C, D and See additional file 4). Additionally, skipping of exon 1B had a significant positive effect on the luciferase activity in the case of the L and S constructs, but was less pronounced for the M construct. This implies that not only 5' UTR length but also alterations in secondary structure, elicited by exon 1B skipping, may influence heterologous protein production. Moreover, HH signaling activation had no apparent influence on the observed luciferase activities. Collectively, these data suggest that the choice of transcription start site, in combination with exon 1B inclusion/skipping, influences the capacity of the Gli1 variants to produce a protein product.
Finally, we investigated whether this differential capacity of the 5' UTRs is due to modulations of the mRNA stability or the translational efficiency. To achieve this, the endogenous Gli1 levels in wtMEFs were enriched by SAG, followed by inhibition of the process of transcription with actinomycin D treatment and subsequent mRNA measurements. The results revealed no significant differences in mRNA stability among the Gli1 variants (See additional file 4). Thus, it may be suggested that the alternative 5' UTRs are likely to regulate Gli1 protein levels via translational mechanisms.
In this report the 5' UTR regions of the mouse Gli1 gene were analyzed, resulting in the identification of two novel transcription start sites and sequence polymorphisms, which control exon 1B inclusion/skipping. Our findings demonstrate that SNPs in the 5' splice site of intron 1B and a SINE B2 insertion in exon 1B have major effects in determining the expression ratio of the 5' UTR variant mRNAs, which distinctly influence translational efficiency (Figure 7).
Initially, 5' RACE and RT-PCR analysis clearly showed 5' variations in the Gli1 transcripts (Figure 1). However, the previously reported exon 2 - exon 3 skipping in human normal tissues/tumor cell lines , exon 3 - partial exon 4 skipping in human cancer cells and especially glioblastomas  and exon 1A - exon 1B skipping in mouse  could not be detected. Moreover, the 5' RACE analysis revealed the presence of some minor transcripts, apparently transcribed from alternative TSSs (data not shown, DDBJ: AB232673 - AB232676). As their expression levels were quite low, we chose to focus on the relatively abundant Gli1 variants in order to analyze the impact of alternative transcripts on HH signaling.
We demonstrated that exon 1B inclusion/skipping is intimately related with SNPs in the donor site of intron 1B (Figures 1, 2, 3, 4 and 5). The GT allele contains a canonical 5' splice site sequence, allowing efficient exon 1B definition and consequently enhances the inclusion of exon 1B, as observed in NIH3T3 cells. On the other hand, the GC allele has a weaker 5' splice site sequence, resulting in a less efficient exon 1B definition and therefore promotes the skipping of exon 1B, as observed in wtMEFs and Sufu-/- MEFs. Moreover, Ptch1-/- MEFs, which are heterozygous and contain both the GT and GC alleles, expressed comparable levels of the exon 1B-included and -skipped variants.
Since the relative expression pattern of the Gli1 variants was not affected by either HH signaling activity (Figure 3A and 3B) or tumorigenesis (Figure 3C)/embryogenesis (Figure 2B), the identified polymorphisms are primary determinants in controlling the splicing regulation of the Gli1 5' UTRs. SNPs/mutations affecting splicing have been reported for consensus splice sites, but were also identified at significant distances from splice junctions [28, 29]. Our own analysis of the 284 reported mutations of PTCH1 revealed the presence of 20 intronic splice changes . These findings highlight the importance of SNPs/mutations in altering splicing patterns, as these may occur not only at canonical splice sites but also at exonic/intronic splicing enhancers/silencers.
Additionally, evolutionary comparisons indicated that GC is likely to represent the ancestral intronic dinucleotide at the 5' splice site of intron 1B, since it is conserved in rat and other species, with the GT substitution and the SINE B2 insertion occurring at later stages (See additional file 5 and Figure 5). Interestingly, in the primate lineage the exon 1B region is characterized by an Alu insertion. These facts suggest that this genomic segment might be genetically unstable and a hotspot for transposon insertion.
5' UTRs are known to regulate protein expression via modulation of mRNA stability and/or translational efficiency . We analyzed the alternative Gli1 5' UTRs in various cell lines and found that 5' end shortening as well as skipping of exon 1B increased their capacity for heterologous protein expression (Figure 6). These observations are in line with previous claims purporting that shorter 5' UTRs of Gli1 are more capable of efficient translation , and support the notion that alternative events in 5' UTRs of mammalian genes are likely to contribute to the regulation of translation . Moreover, mRNA stability assays of the alternative Gli1 transcripts revealed that variations in the 5' UTRs did not affect the pattern of RNA degradation, and consequently, these untranslated sequences regulate Gli1 protein levels by modulating translational efficiency. Interestingly and in line with these observations, Hedgehog signaling-dependent mouse models for medulloblastoma development are apparently influenced by the Gli1 genotype. Deletion of one Ptch1 allele in C57BL/6 mice, which are homozygous for the GC allele, results in a higher incidence of medulloblastomas compared to mice with a mixed genetic background .
Secondary structure prediction of the Gli1 5' UTRs by Mfold highlighted the presence of a long stem-loop structure in L and M that is retained in MΔ1B but not in LΔ1B and may have a role in the differential translatabilities of MΔ1B versus LΔ1B (See additional file 3). However, free energy, ΔG, calculations were not fully in line with the capacity for translation. In addition, we examined whether two other parameters, G-quadruplex (G4) DNA/RNA structures and upstream ORF (upORFs), might influence the translational efficiency (See additional file 6). Nucleotide region 251-294 has a G4 motif sequence G3-N1-7-G3-N1-7-G3-N1-7-G3, and nucleotide region 22-55 has a similar sequence with a potential G4 structure. Moreover, the genome wide G4 DNA database QuadBase  predicted an antisense G4 motif at nucleotide region 322-360. G4-structures are formed not only on DNA but also on RNA , and a G4 motif on the NRAS mRNA was reported to suppress translation . Additionally, the exon 1B-included variants, L, M and S, have an upORF in the same frame as the Gli1 ORF (fully upstream, 46 encoded amino acids), while the exon 1B-skipped variants, LΔ1B, MΔ1B and SΔ1B, have an upORF (38 encoded amino acids) that overlaps with the Gli1 ORF. Although no significant differences between fully upstream and overlapping upORF could be identified in a recent report, long cap-to-upORF distances were found to increase translational inhibition . Wang and Rothnagel have used a 5' UTR construct (alfa-UTR), which is almost equivalent to the S construct in this report, mutated at four ATG codons and apparently eliminating 46 amino acids upORF that we have identified, and observed increased reporter activity, in line with the above predictions . Thus, the combinations of upORFs with G4 structures in the Gli1 5' UTRs are likely to have a role as mediators of the observed patterns of translation.
Our findings highlight the complex posttranscriptional regulation of the mouse Gli1 oncogene. mRNA variants with alternative 5' UTRs were identified, mechanisms that control their expression levels were dissected, and the differential impact of the 5' UTRs on protein synthesis was determined. Moreover, the demonstrated strain differences in regulatory controls of this oncogene suggest that these may have a role in modulating tumor susceptibility in mouse models.
RACE and PCR
5' RACE was performed by using the GeneRacer kit (Life Technologies, CA, USA), with mouse Gli1 exon 4 reverse primers (MWG-Biotech, Ebergsberg, Germany) (Table 1). The RACE products were analyzed on a 4% NuSieve 3:1 agarose gel (FMC BioProducts, ME, USA) and verified by PCR direct sequencing or sequencing of TA-clones in the pGEM-T vector (Promega, WI, USA). Pairs of initial and nested primers were also designed within mouse Gli1 exon 1, as shown in Table 1, and used in combination with the RACE exon 4 primers. The nested PCR analysis and sequence-verification were carried out as described in previous reports [37, 38].
The murine fibroblast cell lines NIH3T3, Ptch1-/- MEFs , wtMEFs, Sufu-/- MEFs and C3H/10T1/2 were cultured as described before [37, 40, 41]. Cells were treated with the S mo ag onist SAG at a concentration of 100 nM, with the medium changed to low serum (0.5% FBS or 1% FBS for wtMEFs), and allowed to grow for an additional 2 days.
The use of animals was approved by the Stockholm South Animal Ethics Committee. The mice were kept at the animal facility of the Karolinska University Hospital, according to local and national regulations. The Sufu+/-Trp53-/- mice were generated by intercrossing Sufu+/-  and Trp53+/- mice . The Ptch1+/- mouse strain has been described previously . Both the Sufu+/-Trp53-/- and the Ptch1+/- strains were maintained on a C57BL/6 genetic background.
Isolation of cerebellum cells
Normal cerebella and medulloblastoma tumors from Sufu+/-Trp53-/- mice were digested with papain, triturated to obtain single-cell suspensions and then centrifuged through a 35%-65% Percoll gradient. Cells from the 35%-65% interface were suspended in Neurobasal medium (Life Technologies, CA, USA). Isolated granule cells were counted and checked with a microscope.
RNA isolation, and real-time RT-PCR
Total RNA was isolated from cells, tissues and mouse embryos, using the RNeasy kit (Qiagen GmbH, Hilden, Germany) according to the manufacturer's protocol. Real-time RT-PCR was performed as described before . Dissociation curves were generated after each PCR run to ensure that a single, specific product was amplified. The results were analyzed with the comparative C ycle t hreshold (Ct) method. For normalization, we used the expression level of Glyceraldehyde-3-phosphate dehydrogenase (Gapdh), and/or Acidic ribosomal protein (Arp). The PCR primers are shown in Table 2.
Analysis of polymorphic variants
Genomic DNA was amplified with pairs of initial and nested primers, which are flanking the exon 1B region and are listed in Table 3. The PCR products were purified (Qiagen), and directly sequenced (MWG-Biotech). The obtained sequences were compared with the Mus musculus chromosome 10 genomic contig, strain C57BL/6J (GenBank: NT_039500).
Functional analysis of 5' UTRs
For generating the 5' UTR constructs, we PCR amplified the selected regions by using specific RACE products as templates, with the primers (MWG-Biotech) listed in Table 4. The PCR products were then digested with the NheI restriction enzyme and cloned upstream of the Renilla luciferase open reading frame of the psiCHECK2 vector (Promega, GenBank: AY535007), which is under the control of the SV40 early enhancer/promoter. All constructs were verified by sequencing using BigDye Terminator v1.1 Cycle Sequencing Kits and an ABI prism DNA sequencer (Life Technologies).
Two hundred ng of each of the 5' UTR constructs were transfected into NIH3T3 or wtMEF cells, with or without SAG treatment initiated 24 hours after transfection, and Ptch1-/- or Sufu-/- MEFs using the FuGENE6 (Roche Diagnostics, Basel, Switzerland) transfection reagent. The activities of Renilla and Firefly luciferases were determined by using the dual-luciferase reporter assay system (Promega) with a FB12 Luminometer (Berthold Detection System, Pforzheim, Germany) or an Infinite M200 (Tecan, Männedorf, Switzerland) according to the manufacturer's recommendations. The experiments and individual measurements were performed at least twice.
Glioma associated oncogene 1
Suppressor of Fused
acid ribosomal protein
days post coitum
long 5' UTR transcript
medium 5' UTR transcript
short 5' UTR transcript
L transcript with skipped exon 1B
M transcript with skipped exon 1B
S transcript with skipped exon 1B
single nucleotide polymorphism
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
Davuluri RV, Suzuki Y, Sugano S, Plass C, Huang TH: The functional consequences of alternative promoter use in mammalian genomes. Trends Genet. 2008, 24 (4): 167-177. 10.1016/j.tig.2008.01.008
Nagasaki H, Arita M, Nishizawa T, Suwa M, Gotoh O: Species-specific variation of alternative splicing and transcriptional initiation in six eukaryotes. Gene. 2005, 364: 53-62. 10.1016/j.gene.2005.07.027
Kwan T, Benovoy D, Dias C, Gurd S, Provencher C, Beaulieu P, Hudson TJ, Sladek R, Majewski J: Genome-wide analysis of transcript isoform variation in humans. Nat Genet. 2008, 40 (2): 225-231. 10.1038/ng.2007.57
Hughes TA: Regulation of gene expression by alternative untranslated regions. Trends Genet. 2006, 22 (3): 119-122. 10.1016/j.tig.2006.01.001
Pickering BM, Willis AE: The implications of structured 5' untranslated regions on translation and disease. Semin Cell Dev Biol. 2005, 16 (1): 39-47. 10.1016/j.semcdb.2004.11.006
Jiang J, Hui CC: Hedgehog signaling in development and cancer. Dev Cell. 2008, 15 (6): 801-812. 10.1016/j.devcel.2008.11.010
Hahn H, Wicking C, Zaphiropoulos PG, Gailani MR, Shanley S, Chidambaram A, Vorechovsky I, Holmberg E, Unden AB, Gillies S, et al: Mutations of the human homolog of Drosophila patched in the nevoid basal cell carcinoma syndrome. Cell. 1996, 85 (6): 841-851. 10.1016/S0092-8674(00)81268-4
Berman DM, Karhadkar SS, Maitra A, Montes De Oca R, Gerstenblith MR, Briggs K, Parker AR, Shimada Y, Eshleman JR, Watkins DN, et al: Widespread requirement for Hedgehog ligand stimulation in growth of digestive tract tumours. Nature. 2003, 425 (6960): 846-851. 10.1038/nature01972
Sanchez P, Hernandez AM, Stecca B, Kahler AJ, DeGueme AM, Barrett A, Beyna M, Datta MW, Datta S, Ruiz i Altaba A: Inhibition of prostate cancer proliferation by interference with SONIC HEDGEHOG-GLI1 signaling. Proc Natl Acad Sci USA. 2004, 101 (34): 12561-12566. 10.1073/pnas.0404956101
Yuan Z, Goetz JA, Singh S, Ogden SK, Petty WJ, Black CC, Memoli VA, Dmitrovsky E, Robbins DJ: Frequent requirement of hedgehog signaling in non-small cell lung carcinoma. Oncogene. 2007, 26 (7): 1046-1055. 10.1038/sj.onc.1209860
Zaphiropoulos PG, Shimokawa T: Splicing variations in components of Patched/Hedgehog signaling. Trends Cell Mol Biol. 2008, 3: 31-35.
Nybakken K, Vokes SA, Lin TY, McMahon AP, Perrimon N: A genome-wide RNA interference screen in Drosophila melanogaster cells for new components of the Hh signaling pathway. Nat Genet. 2005, 37 (12): 1323-1332. 10.1038/ng1682
Kinzler KW, Ruppert JM, Bigner SH, Vogelstein B: The GLI gene is a member of the Kruppel family of zinc finger proteins. Nature. 1988, 332 (6162): 371-374. 10.1038/332371a0
Ruiz i Altaba A: The works of GLI and the power of hedgehog. Nat Cell Biol. 1999, 1 (6): E147-148. 10.1038/14099
Nilsson M, Unden AB, Krause D, Malmqwist U, Raza K, Zaphiropoulos PG, Toftgård R: Induction of basal cell carcinomas and trichoepitheliomas in mice overexpressing GLI-1. Proc Natl Acad Sci USA. 2000, 97 (7): 3438-3443. 10.1073/pnas.050467397
Wang XQ, Rothnagel JA: Post-transcriptional regulation of the gli1 oncogene by the expression of alternative 5' untranslated regions. J Biol Chem. 2001, 276 (2): 1311-1316. 10.1074/jbc.M005191200
Shimokawa T, Tostar U, Lauth M, Palaniswamy R, Kasper M, Toftgård R, Zaphiropoulos PG: Novel human glioma-associated oncogene 1 (GLI1) splice variants reveal distinct mechanisms in the terminal transduction of the hedgehog signal. J Biol Chem. 2008, 283 (21): 14345-14354. 10.1074/jbc.M800299200
Lo HW, Zhu H, Cao X, Aldrich A, Ali-Osman F: A Novel splice variant of GLI1 that promotes glioblastoma cell migration and invasion. Cancer Res. 2009, 69 (17): 6790-6798. 10.1158/0008-5472.CAN-09-0886
Chen JK, Taipale J, Young KE, Maiti T, Beachy PA: Small molecule modulation of Smoothened activity. Proc Natl Acad Sci USA. 2002, 99 (22): 14071-14076. 10.1073/pnas.182542899
Goodrich LV, Milenkovic L, Higgins KM, Scott MP: Altered neural cell fates and medulloblastoma in mouse patched mutants. Science. 1997, 277 (5329): 1109-1113. 10.1126/science.277.5329.1109
Kimura H, Stephen D, Joyner A, Curran T: Gli1 is important for medulloblastoma formation in Ptc1+/- mice. Oncogene. 2005, 24 (25): 4026-4036.
Sterner DA, Carlo T, Berget SM: Architectural limits on split genes. Proc Natl Acad Sci USA. 1996, 93 (26): 15081-15085. 10.1073/pnas.93.26.15081
Desmet FO, Hamroun D, Lalande M, Collod-Beroud G, Claustres M, Beroud C: Human Splicing Finder: an online bioinformatics tool to predict splicing signals. Nucleic Acids Res. 2009, 37 (9): e67- 10.1093/nar/gkp215
Reese MG, Eeckman FH, Kulp D, Haussler D: Improved splice site detection in Genie. J Comput Biol. 1997, 4 (3): 311-323. 10.1089/cmb.1997.4.311
Nakamura T, Aikawa T, Iwamoto-Enomoto M, Iwamoto M, Higuchi Y, Pacifici M, Kinto N, Yamaguchi A, Noji S, Kurisu K, et al: Induction of osteogenic differentiation by hedgehog proteins. Biochem Biophys Res Commun. 1997, 237 (2): 465-469. 10.1006/bbrc.1997.7156
Cooper TA, Wan L, Dreyfuss G: RNA and disease. Cell. 2009, 136 (4): 777-793. 10.1016/j.cell.2009.02.011
Hull J, Campino S, Rowlands K, Chan MS, Copley RR, Taylor MS, Rockett K, Elvidge G, Keating B, Knight J, et al: Identification of common genetic variation that modulates alternative splicing. PLoS Genet. 2007, 3 (6): e99- 10.1371/journal.pgen.0030099
Lindström E, Shimokawa T, Toftgård R, Zaphiropoulos PG: PTCH mutations: distribution and analyses. Hum Mutat. 2006, 27 (3): 215-219. 10.1002/humu.20296
Resch AM, Ogurtsov AY, Rogozin IB, Shabalina SA, Koonin EV: Evolution of alternative and constitutive regions of mammalian 5'UTRs. BMC Genomics. 2009, 10: 162- 10.1186/1471-2164-10-162
Svärd J, Rozell B, Toftgård R, Teglund S: Tumor suppressor gene cooperativity in compound Patched1 and suppressor of fused heterozygous mutant mice. Mol Carcinog. 2009, 48 (5): 408-419. 10.1002/mc.20479
Yadav VK, Abraham JK, Mani P, Kulshrestha R, Chowdhury S: QuadBase: genome-wide database of G4 DNA--occurrence and conservation in human, chimpanzee, mouse and rat promoters and 146 microbes. Nucleic Acids Res. 2008, D381-385. 36 Database
Patel DJ, Phan AT, Kuryavyi V: Human telomere, oncogenic promoter and 5'-UTR G-quadruplexes: diverse higher order DNA and RNA targets for cancer therapeutics. Nucleic Acids Res. 2007, 35 (22): 7429-7455. 10.1093/nar/gkm711
Kumari S, Bugaut A, Huppert JL, Balasubramanian S: An RNA G-quadruplex in the 5' UTR of the NRAS proto-oncogene modulates translation. Nat Chem Biol. 2007, 3 (4): 218-221. 10.1038/nchembio864
Calvo SE, Pagliarini DJ, Mootha VK: Upstream open reading frames cause widespread reduction of protein expression and are polymorphic among humans. Proc Natl Acad Sci USA. 2009, 106 (18): 7507-7512. 10.1073/pnas.0810916106
Shimokawa T, Rahnama F, Zaphiropoulos PG: A novel first exon of the Patched1 gene is upregulated by Hedgehog signaling resulting in a protein with pathway inhibitory functions. FEBS Lett. 2004, 578 (1-2): 157-162. 10.1016/j.febslet.2004.11.006
Shimokawa T, Svärd J, Heby-Henricson K, Teglund S, Toftgård R, Zaphiropoulos PG: Distinct roles of first exon variants of the tumor-suppressor Patched1 in Hedgehog signaling. Oncogene. 2007, 26 (34): 4889-4896. 10.1038/sj.onc.1210301
Taipale J, Chen JK, Cooper MK, Wang B, Mann RK, Milenkovic L, Scott MP, Beachy PA: Effects of oncogenic mutations in Smoothened and Patched can be reversed by cyclopamine. Nature. 2000, 406 (6799): 1005-1009. 10.1038/35023008
Rahnama F, Shimokawa T, Lauth M, Finta C, Kogerman P, Teglund S, Toftgård R, Zaphiropoulos PG: Inhibition of GLI1 gene activation by Patched1. Biochem J. 2006, 394 (Pt 1): 19-26.
Svärd J, Heby-Henricson K, Persson-Lek M, Rozell B, Lauth M, Bergström A, Ericson J, Toftgård R, Teglund S: Genetic elimination of Suppressor of fused reveals an essential repressor function in the mammalian Hedgehog signaling pathway. Dev Cell. 2006, 10 (2): 187-197. 10.1016/j.devcel.2005.12.013
Donehower LA, Harvey M, Slagle BL, McArthur MJ, Montgomery CA, Butel JS, Bradley A: Mice deficient for p53 are developmentally normal but susceptible to spontaneous tumours. Nature. 1992, 356 (6366): 215-221. 10.1038/356215a0
We thank J. Svärd for help with the Ptch+/- medulloblastoma samples. Professor Rune Toftgård is acknowledged for generous support. We are grateful to V. Jaks, M. Kasper, I. Sur, C. Finta, R. Saraswathi and Å. Bergström for helpful discussions and technical support. This study was supported by the Swedish Cancer Fund, the Swedish Research Council, the Swedish Childhood Cancer Foundation and the Magnus Bergvalls Foundation. RP was partially supported by Karolinska Institutet' s funding for postgraduate students. TS was supported by a Marie Curie International Incoming Fellowship.
PZ and TS planned the work and supervised RP. RP was responsible for all experiments. ST and ML contributed to the experiments using mice. All authors read and approved the final version of the manuscript.
Electronic supplementary material
Additional file 1: Expression pattern of Gli1 variants in mouse embryos. Additional figure 1, additional method, additional reference and additional table. (PDF 330 KB)
Additional file 3: Secondary structure characteristics of Gli1 5' UTRs. Additional figure 3 and additional reference. (PDF 781 KB)
Additional file 5: Comparison of the exon 1B/intron 1B junctions among different species. Additional figure 5. (PDF 67 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Palaniswamy, R., Teglund, S., Lauth, M. et al. Genetic variations regulate alternative splicing in the 5' untranslated regions of the mouse glioma-associated oncogene 1, Gli1. BMC Molecular Biol 11, 32 (2010). https://doi.org/10.1186/1471-2199-11-32
- NIH3T3 Cell
- Mixed Genetic Background
- Splice Site Sequence
- Relative Expression Pattern