| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Molecular and Cellular Biology, January 2008, p. 883-895, Vol. 28, No. 2
0270-7306/08/$08.00+0 doi:10.1128/MCB.01345-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Cancer Institute of New Jersey, Department of Pharmacology, Robert Wood Johnson Medical School, University of Medicine and Dentistry of New Jersey, New Brunswick, New Jersey 08901
Received 26 July 2007/ Returned for modification 24 August 2007/ Accepted 2 November 2007
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Decades of elegant studies have defined the basic splicing mechanism, which requires the interaction of spliceosomal core components with splicing signals resident within the pre-mRNA sequences. However, the regulation of spliceosome recruitment and splice site recognition is highly complex and only partly understood. Current models suggest that splice site recognition is under the control of both specific RNA regulatory motifs and at least two classes of RNA binding proteins. The first class includes the heterogenous nuclear ribonucleoproteins (hnRNPs) (33), which generally function as splicing repressors through their recognition of exonic and intronic silencer elements. Interactions of hnRNPs with silencer elements repress proximal exon recognition by inhibiting the use of adjacent splice sites, either by antagonizing the function of positive regulatory elements or by recruiting factors that directly interfere with the splicing machinery (9, 58). Indeed, the human genome is rich in silencing elements (10), suggesting a prevalent role of negative regulation in exon definition (54, 55).
The second class of RNA binding proteins is the serine-arginine-rich (SR) protein family (31). SR proteins are essential for multiple steps of spliceosome assembly at both the 3' and 5' splice sites and can function in both constitutive and alternative splicing (14), although how they accomplish this remains largely unclear. What has been shown is that SR proteins bind to either exonic or intronic splicing enhancers, promoting "exon recognition" by directly recruiting the splicing machinery and/or by antagonizing the action of nearby silencer elements (58). Many alternatively spliced exons, small exons, or exons with weak splice sites rely upon the activity of enhancer-bound SR proteins for their inclusion (56). Therefore, SR proteins such as ASF/SF2 and SC35 are generally considered splicing activators (14).
Both SR and hnRNP proteins recognize a large repertoire of target sequences. They have also been shown to modulate alternative splice site selection in a concentration-dependent manner (2, 15); competition between SR proteins and hnRNPs for interaction with regulatory sequences and/or constitutive spliceosomal components, achieved through either qualitative or quantitative changes in these splicing factors, can result in a subtle shift in the balance between positive and negative regulation (3). Adding to this subtle rheostat for splicing control are a number of auxiliary proteins that help to regulate the regulators. Finally, under some circumstances, SR proteins can act as repressors and hnRNPs can activate splicing, adding an additional layer of complexity to this already elaborate process (6, 18).
Despite these complexities, considerable progress has been made in the identification of trans-acting regulators and cis-acting elements involved in the choice of splice sites for a few specific genes, and recent genome-wide computational annotations on expressed sequence tags (35, 43, 54) as well as SELEX analyses (16, 25, 49) have begun to provide a broader view of splice site regulation. Nevertheless, the detailed molecular mechanism of alternative splicing regulation remains largely unclear, hampered by a loose definition of splice site sequences in higher eukaryotes, a high degeneracy of RNA regulatory elements, a limited understanding of RNA secondary structure, the extreme complexity of spliceosomal components and functions and, particularly, a paucity of means to manipulate and control alternative splicing in an experimental system.
In this report, we identify a subset of alternative splicing events in cancer-related genes that are coordinately regulated by the highly consumed psychoactive drug caffeine. The particular splice variants of these genes appear to be tumor specific, in that the alternate transcripts have not been identified in normal cells. Caffeine-mediated induction is rapid, time and concentration dependent, and reversible. Focusing on the tumor suppressor gene KLF6 as a prototype for this caffeine-regulated subset, we show that the SR protein SC35 is required for alternative splicing of this gene. Importantly, caffeine induces SC35 levels, implicating this protein in caffeine-mediated alternative splice site selection. That a small highly diffusible molecule widely present in the human diet can have such profound effects on alternative splicing is an entirely novel and highly significant observation. We propose that caffeine is mimicking an endogenous pathway that is activated in certain cancer cells to induce alternative splicing of a specific group of genes; thus, these studies identify both a unique tool and a unique subset of genes to investigate the mechanisms underlying alternative splicing associated with the cancer phenotype.
| MATERIALS AND METHODS |
|---|
|
|
|---|
-minimal essential medium. Human cervical carcinoma HeLa cells, human fetal lung fibroblast WI38 cells, the rat glioma cell line RG-2, and the mouse melanoma cell line B16-F0 were maintained in Dulbecco's modified Eagle's medium. In all cases, medium was supplemented with 10% (vol/vol) fetal bovine serum, 2.0 mM glutamine, 100 units/ml penicillin, and 100 µg/ml streptomycin. One day prior to treatment, 3A cells were transferred from 37°C to 40°C to induce a trophoblast-like phenotype. Semiquantitative RT-PCR assays. Cells were seeded in six-well plates at a density of 2.0 x 105 to 2.5 x 105 cells/well. After 18 to 24 h, cells were either untreated or treated with caffeine (8 to 14 mM) for 1 to 24 h as indicated in the figure legends. For inhibition of protein synthesis, 40 µg/ml of cycloheximide was added to cells in the presence or absence of 14 mM caffeine. Either 75 ng or 300 ng of total RNA, prepared using TRIzol reagent (Invitrogen), was analyzed using SS one-step reverse transcription-PCR reverse transcription-PCR (RT-PCR) reagents (Invitrogen) according to the manufacturer's recommendations. The number of amplification cycles used for each reaction was predetermined to ensure that transcript amplification was within a linear range. Primers used for KLF6 RT-PCR are shown below. The primers for β2-microglobulin have been described previously (42). Primers for microarray data validation are available upon request.
Real-time RT-PCR. Total RNA was extracted from either caffeine-treated or untreated cells and reverse transcribed using MultiScribe reverse transcriptase (Applied Biosystems, CA). TaqMan real-time PCR was performed using primers and probes either specific to total KLF6 transcripts (primer F, 5'-CCTCCACGCCTCCATCTT-3'; R, 5'-GCACGCAACCCCACAGT-3'; probe, 5'-FAM [6-carboxyfluorescein]-TCCGGAACTGAGCAGGGAACCTTCT-BHQ1 [Black Hole Quencher 1])-3'), or specific to SpKLF6 (F, 5'-GCTAGAAAGTTTGAAATGGTGGG-3'; R, 5'-CGTTCCAGCTCTAGGCAGG-3'; probe, 5'-FAM-TGCTGATCTCGGCAGTTGTAGCTTCTTG-BHQ1-3') using the ABI Prism 7700 sequence detection system (Applied Biosystems, Foster City, CA). Each sample was thermocycled in duplicate at three different concentrations of cDNA (10, 1, and 0.1 ng) to verify the linearity of the assay. Each experiment was repeated twice. Data were analyzed using ABI Prism SDS software.
RNA interference assay. A total of 3 x 104 HeLa cells in 0.5 ml of Dulbecco's modified Eagle's medium without antibiotics were seeded into each well of a 24-well plate. After 18 h, 25 to 100 nM small interfering RNA (siRNA; siGENOME SMARTpool reagent; Dharmacon) was mixed with Oligofectamine reagent (Invitrogen), incubated at room temperature for 20 min, and added drop-wise to cells. At 48, 72, and 96 h post-siRNA transfection, cells were harvested for either whole-cell lysates or RNA isolation. The concentration and time required for optimal downregulation were determined empirically by both RT-PCR and Western blot analysis.
Western blot analysis. Cells were washed twice with cold phosphate-buffered saline and lysed in radioimmunoprecipitation buffer (50 mM Tris-HCl, pH 7.4, 150 mM NaCl, 1.0% Triton X-100, 0.1% sodium dodecyl sulfate, 1% deoxycholate, Na salt) plus protease inhibitor cocktail (Roche Diagnostics), 1 mM sodium fluoride, 1 mM sodium orthovanadate, and 100 mg/ml phenylmethylsulfonyl fluoride. The protein concentration of cell lysates was determined using the Pierce bicinchoninic acid protein assay. Equal amounts of total protein (15 to 25 µg) were analyzed by 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis followed by immunoblotting using goat polyclonal antibody against hUpf1 (1:1,000; Abcam), mouse monoclonal antibody against SC35 (1:20; provided by Cyril Bourgois and James Stevenin), mouse anti-c-Myc monoclonal antibody (1:1,000; Abcam), or goat polyclonal anti-β-actin antibody (1:1,000; Santa Cruz Biotechnology). The secondary antibody was either horseradish peroxidase-conjugated donkey anti-goat immunoglobulin G (IgG; Santa Cruz Biotechnology) or horseradish peroxidase-conjugated goat-anti-mouse IgG (1:2,500; Amersham). Immunoreactive bands were visualized using an enhanced chemiluminescent system (Amersham) according to the manufacturer's recommendations.
Microarray analysis. The Cancer-Specific Splice Variant array was purchased from ArrayIT (Sunnyvale, CA; this specific product has been discontinued but similar products can be obtained from ExonHit Therapeutics [Gaithersburg, MD]). This array contains 524 pairs of oligonucleotides homologous to the intron-exon boundary sequence of the 524 pairs of isoforms; one member of each pair has been identified in normal cells, while the second member of the pair was found in tumor cells. The sequences of these oligos can be found on the website of the UCLA Alternative Splicing Annotation Project (http://bioinfo.mbi.ucla.edu/ASAP/). HeLa cells were seeded in 100-mm dishes, incubated overnight, and then treated with 14 mM caffeine. Total RNA was harvested 24 h later in TRIzol reagent (Invitrogen) following the manufacturer's protocol until the chloroform extraction step. The postextraction aqueous phase was mixed with 1 volume of 70% ethanol and immediately subjected to the RNeasy mini column (Qiagen), purified, and labeled as follows: 15 µg of total RNA from either untreated or caffeine-treated HeLa cells was labeled with either Cy5 or Cy3 and subjected to array hybridization. The probe labeling, array hybridization, washing, scanning, and data retrieval steps were performed by the Gene Expression Core Facility (RWJMS—Cancer Institute of New Jersey). Positive controls were included on the array to establish criteria for data mining, and signal differences of <2-fold were considered background.
Minigene and minigene deletion/mutation constructs.
The wild-type KLF6 minigene was created by inserting a PCR-generated fragment containing exon 1 through exon 2 of the KLF6 gene between the KpnI and BamHI sites of pcDNA 3.1(+) (Invitrogen). A PCR was performed using genomic DNA purified from HeLa cells; the primers for this PCR are available upon request. To prepare deletion constructs, the KLF6 minigene plasmid DNA was digested with ClaI and PmlI to remove a large portion of introns 1 and 1a and exon 1a, leaving 146 nucleotides (nt) downstream of exon 1 and 24 nt upstream of exon 2 within the pcDNA 3.1 backbone (pKLF6-
intron 1/1a/exon 1a) (see Fig. 5, below). Various sizes of fragments containing exon 1a and different lengths of flanking intronic sequences (indicated in Fig. 5B, below) were generated by PCR and cloned into pKLF6-
intron 1/1a/exon 1a at the ClaI/PmlI sites. Point mutations were introduced into the deletion construct (De-4) using the QuikChange site-directed mutagenesis kit (Stratagene) according to the manufacturer's recommendations. All the constructs were verified by sequencing.
|
Minigene transfection and splicing assays. A total of 2.0 x 105 to 2.5 x 105 Hela cells/well were seeded into six-well plates. After 18 h, 0.1 to 0.2 µg of minigene construct supplemented with salmon sperm DNA for a total amount of 2 µg was transfected into cells using Lipofectamine 2000 reagent (Invitrogen). Following a 24-h incubation, cells were shifted to caffeine-containing medium and incubated for an additional 24 h prior to RNA isolation using TRIzol reagent (Invitrogen). KLF6 splice variants were assayed by RT-PCR using either primers containing both vector and KLF6 sequences (specific to KLF6 minigenes) or β-globin primers located in exons 1 and 2 and specific to the minigene-generated transcript.
SC35 mammalian expression vector. SC35 cDNA was amplified by RT-PCR from total RNA prepared from HeLa cells using the forward primer 5'-CGGGAATTCATGAGCTACGGCCGCCCCCC-3' and the reverse primer 5'-CGCCGGATCCAAAGAGGACACCGCTCCTTCCTC-3'. This PCR fragment was purified and digested with BamHI and EcoRI and then cloned into pcDNA3.1/myc-His(–)-B at the BamHI/EcoRI sites. The c-Myc and His tags are in frame with SC35 at the C terminus of SC35 cDNA. The construct (pcDNA3.1-SC35-c-Myc/His) was verified by sequencing.
To generate the SC35 tet-on construct, the SC35-c-Myc-His DNA fragment was amplified by PCR from the pcDNA3.1-SC35-c-Myc/His plasmid while adaptive restrictive enzyme sites EcoRI/ClaI were introduced at the ends. This PCR fragment was purified, inserted into pTRE-tight vector (Clontech), and verified by sequencing. The HeLa SC35 tet-on inducible clone 6-4 and vector (pTRE-tight)-only clone 4-1 were established following instructions from the manufacturer (Clontech).
SC35 overexpression and minigene splicing assay. HeLa cells were seeded at a density of 3.0 x 105 cells/well in six-well plates and maintained for 6 h prior to transient transfection with increasing concentrations (0.5 to 1.5 µg) of pcDNA3.1-SC35-c-Myc/His and 0.2 µg of KLF6 minigene. Transfection was performed using Lipofectamine 2000 reagent (Invitrogen) according to the manufacturer's recommendations. Either whole-cell lysates or total RNA was prepared as described above at 26 h posttransfection. Whole-cell lysates were subjected to Western blot analysis using anti-c-Myc antibody (1:1,000; Abcam). Total RNA was assayed by RT-PCR using either primers containing both vector and KLF6 sequences (specific to KLF6 minigenes) or β-globin primers located in exons 1 and 2 and specific to the minigene-generated transcript.
Nuclear extract, [32P]RNA probes, and RNA electrophoretic mobility shift assay (EMSA).
The HeLa SC35 tet-on clone 6-4 and vector-only tet-on clone 4-1 were treated with 1 µg/ml doxycycline for 48 h, and nuclear extracts were purified from these cells as described previously (21). The protein concentration was determined using a Bio-Rad protein assay kit. Expression of SC35-c-Myc was examined by Western blot analysis (see Fig. 8E). Wild-type and mutant RNA probes were transcribed from single-stranded oligonucleotide templates containing the 22-nt T7 RNA polymerase promoter annealed to a complementary primer. The T7 MEGAshortscript kit (Ambion) was used to transcribe RNAs with the following changes: the reaction mixture contained 750 µM each of CTP, UTP, and GTP, 75 µM of ATP, and 50 µCi of [
-32P]ATP (GE-Amersham Biosciences) for labeling reactions. Transcripts were purified using G-25 Sephadex columns (Roche). Each 30-µl RNA-protein binding reaction mixture contained 20 µg of nuclear extract incubated with 5 x 105 cpm [32P]RNA probe at 25°C for 30 min in HEPES binding buffer (10 mM HEPES at pH 7.5, 50 mM KCl, 5 mM MgCl2, 1 mM dithiothreitol, 1 mM phenylmethylsulfonyl fluoride, 16 units SUPERase·In RNase inhibitor [Ambion], 120 µg/ml yeast tRNA, 5% glycerol) in the presence or absence of an
10-fold excess of unlabeled transcript. For supershift assays, 4 to 8 µg of either control IgG or anti-c-Myc antibody was added to the reaction mixture and incubated on ice for 10 min prior to the binding reaction. Reaction mixtures were electrophoresed on 8% nondenaturing polyacrylamide gels at 200 V for 2 h. Radiolabeled complexes were detected by autoradiography on screen-K using the Molecular Imager FX (Bio-Rad).
|
| RESULTS |
|---|
|
|
|---|
300 bp; hereafter designated SpKLF6) created at the expense of the KLF6 transcript (223 bp). This result was confirmed by independent real-time RT-PCR analysis (Fig. 1A, right panel). Sequence analysis of the upper band indicated that it was derived from a transcript identical to KLF6 except for the inclusion of an additional 66 nucleotides that mapped to the intron between exons 1 and 2 (Fig. 1B). Examination of this sequence revealed the presence of consensus splice sites at both the 5' and 3' borders (Fig. 1B); while the 5' splice site fit the consensus, the 3' sequence delimited a "weak" splice site, consistent with a cryptic exon (46) (hereafter designated exon 1a).
|
The induction of SpKLF6 by caffeine was concentration dependent, with modest induction observed at 8 mM and significant expression induced at 14 mM (Fig. 2A). A time course analysis indicated that caffeine-mediated SpKLF6 induction occurred rapidly (as early as 1 hour following caffeine exposure) and was reversible, since SpKLF6 levels were reduced following the removal of caffeine from the medium (Fig. 2B). To determine the generality of this observation, experiments were repeated using several other human tumor cell lines as well as two nontransformed "normal" cell lines, with similar results (Fig. 2C).
|
|
Using RNA isolated from vehicle- and caffeine-treated HeLa cells, we performed an initial probe of this microarray (see Materials and Methods). Interestingly, caffeine induced alternative splicing in
40 of the 524 alternative spliced genes (unpublished data), indicating that caffeine does not impact the general splicing mechanism but is specific for a subset of genes. To validate this observation, semiquantitative RT-PCR was performed on total RNA samples harvested from both caffeine-treated and vehicle-treated HeLa cells, using primers designed to distinguish between the variants. The results for several genes are shown in Fig. 4. Chaperonin-containing TCP1 subunit 3 (CCT3) expresses two transcripts containing either exon 2 or exon 2a, exclusively. While the exon 2a-containing transcript is predominant in untreated HeLa cells, exposure to caffeine resulted in an increase in levels of the exon 2-containing transcript accompanied by a decrease in the level of the exon 2a-containing transcript (this increase in the splice variant at the expense of the normal transcript was similar to what was observed for KLF6/SpKLF6 [Fig. 1A]). Similar changes were observed for the asparagine synthetase (ASNS), COMM-domain containing 5 (COMMD5), and ATP binding cassette subfamily F member 2 (ABCF2) genes. Changes were also observed in the levels of transcripts containing exclusively spliced exons for the yippee-like 5 gene (YPEL5, exons 2 and 3). Finally, similar to what was observed for SpKLF6, exon inclusion events (inclusion of exon 2 in the YPEL5 transcript and exon 2b in the SLC39A1/ZIRTL transcript) were also induced by caffeine.
|
cis-acting elements located within and adjacent to exon 1a are essential for SpKLF6 induction. Although KLF6 splice variants have been identified, there have been no studies that investigated the alternative splicing regulation of this gene. To determine the cis elements regulating KLF6 and SpKLF6 splicing, particularly the element(s) responsive to caffeine, a KLF6 minigene and minigene deletion mutants were generated. These minigenes were transfected into HeLa cells either untreated or treated with caffeine. RT-PCR analyses were performed on total RNA harvested 24 h after caffeine treatment to monitor the splicing pattern of the KLF6 minigenes. Two sets of primers were used. The first set, the same as the one used in Fig. 1A, was complementary to the KLF6 sequence and detected both KLF6 and SpKLF6 transcripts but could not distinguish between endogenous and minigene-generated RNA. The second set of primers, complementary to the junction regions of the vector and KLF6 sequences, was designed to detect only transcripts generated from the minigene. As shown in Fig. 5A (right panel, lane 3), the minigene accurately recapitulated the endogenous KLF6 gene, i.e., KLF6 was the predominant transcript generated from the minigene in untreated cells, and caffeine induced the expression of minigene-derived SpKLF6 transcripts. To delimit intronic sequences required for induction of SpKLF6, minigene deletion constructs that sequentially eliminated exon 1a flanking sequence were analyzed (De-1, -2, -3, -4, -5, and -6) (Fig. 5B). Interestingly, even a deletion construct with only nt –65 to +63 adjacent to exon 1a (De-4) was permissive for SpKLF6 induction. Indeed, when exon 1a and the flanking –65/+63 sequence were inserted into a heterologous context (Dup4-1, human β-globin exons 1 and 2 [34]), this sequence was sufficient to support caffeine-mediated induction of exon1a inclusion [Fig. 5C, DUP4-1 (KLF6)], suggesting that the cis-acting regulatory elements involved in exon 1a inclusion lie within this region. When the exon 1a flanking region was further delimited to –37/+63 (De-5), a higher basal level of SpKLF6 transcript was expressed (Fig. 5B, compare lane 11 to lane 1), which was further induced by caffeine treatment (Fig. 5B, compare lane 12 to lane 11). This result suggests the presence of exon1a silencing elements that are disrupted upon deletion of the sequences from –65 to –37. Notably, De-5 still supports alternative splicing in the presence of caffeine. In contrast, deletion of sequences from +38 to +63 (De-6) inhibited SpKLF6 expression, even in the presence of caffeine (Fig. 5B, lane 14), suggesting the presence of an enhancing element(s) within this region that is likely to be responsive to caffeine.
To test whether the +38 to +63 region contains a caffeine response element (CafRE) that is capable of promoting exon inclusion, we took an approach that was successfully employed in studies of the c-src N1 exon enhancer, using an engineered DUP4-1 minigene (34). The DUP4-1 minigene contains exons 1 and 2 of the human β-globin gene as well as an embedded 33-nt chimeric exon comprised of a segment of the first β-globin exon fused to the second exon (Fig. 5C). DUP4-1 was initially constructed to test a relationship between exon length and constitutive exon skipping; the study showed that the middle 33-nt exon was constitutively skipped until it reached a length of >51 nt (7). Taking advantage of this system, Black and colleagues inserted a c-src N1 enhancing element in triplicate downstream of the 33-nt exon in DUP4-1 and demonstrated the activity of the N1 enhancer in promoting inclusion of the 33-nt middle exon in a complete heterologous context (34). Using the same strategy, the putative CafRE was inserted into the DUP4-1 minigene to generate a heterologous minigene, DUP(en) (Fig. 5C). DUP(en) was transfected into HeLa cells, and minigene splicing assays were performed using total RNA from untreated or caffeine-treated cells. As shown in Fig. 5C, panel 2, caffeine induced inclusion of the heterologous 33-nt hybrid exon in a subset of β-globin transcripts derived from the DUP(en) minigene (7.6%). Although this rate of exon inclusion was modest, it was highly reproducible and is similar to what was observed in analyses of the role of the c-src exon N1 enhancer; in that study, DUP4-1 containing three copies of the c-src exon N1 enhancer supported a low level of exon inclusion when analyzed in HEK cells (34). Overall, our data strongly suggest that activation by caffeine occurs via the CafRE.
To further delineate the CafRE, scanning mutations were performed in the region between +38 and +63 in the context of De-4. As shown in Fig. 6A, mutations in Mu(2) led to dominant expression of SpKLF6, while those in Mu(4) resulted in a significant reduction of exon 1a inclusion in the presence of caffeine, suggesting one silencing motif and one enhancing motif in close proximity within this region. Essential nucleotides for caffeine-mediated exon 1a inclusion were localized within the TTGAG motif. Interestingly, sequence elements with high homology to this motif were found in other caffeine-regulated genes identified in the microarray; notably, all of these elements were found within the flanking introns of the exons regulated by caffeine (Fig. 6B). Extensive experimental analyses are currently under way to test the functional relevance of these putative caffeine response elements.
|
SC35-specific siRNA was transfected into HeLa cells, resulting in a greater-than-80% decrease in SC35 levels by 72 h posttransfection (Fig. 7A, lower panel). Under these conditions of SC35 depletion, the induction of SpKLF6 by caffeine was markedly reduced (Fig. 7A, upper panel), implicating this SR protein in the activation of exon 1a inclusion. To determine whether caffeine had a direct effect on SC35, Western analysis was performed on caffeine-treated cells. As shown in Fig. 7B, the baseline (untreated) level of SC35 remained constant for 24 h (lanes 1, 3, and 5); in contrast, treatment with caffeine resulted in a marked (>6-fold) increase in SC35 protein (lanes 2, 4, and 6). Notably, this increase occurred at a time and concentration commensurate with the induction of SpKLF6. Western blot analysis using a pan-anti-SR antibody (that recognizes most SR domain-containing proteins except ASF/SF2 and SC35), monoclonal antibody MAb104 (that recognizes most phosphorylated SR proteins), or an anti-ASF/SF2 antibody revealed no changes in these proteins following caffeine treatment (data not shown).
|
SC35 binds to the CafRE in vitro. An in vitro EMSA analysis was performed to demonstrate interactions between the CafRE and SC35. As shown in Fig. 8D, when a radio labeled wild-type RNA probe containing the KLF6 CafRE was incubated with nuclear extract from either vector-only or SC35-overexpressing HeLa cells, a unique band was observed compared to what was observed with probe alone (Fig. 8D, compare lane 3 to lane 2). Notably, this band was markedly increased using extracts from SC35-overexpressing cells compared to mock-transfected cell extract (compare lane 3 to lane 4; see Fig. 8E for verification of overexpression). Excess cold wild-type RNA transcript competed effectively for the novel band (lane 5), while a nonspecific cold RNA transcript did not (lane 6), verifying the specificity of the binding. Moreover, when a probe containing the point mutations shown to inhibit induction of SpKLF6 [as in minigene Mu(4)] was incubated with SC35-enriched nuclear extract, no specific band was observed (lane 7). To determine whether SC35 was part of the novel complex, supershift assays were preformed using anti-c-Myc antibody. As shown in lanes 9 and 11, incubation with the specific antibody resulted in a decrease in the unique band and the appearance of a supershifted band, while the control IgG did not affect the specific complex. Taken together, these data indicate that overexpression of SC35 results in increased complex formation at the CafRE. Whether SC35 interacts directly or indirectly with the CafRE is currently under investigation.
| DISCUSSION |
|---|
|
|
|---|
8% of these genes is regulated by caffeine. To begin to evaluate the mechanism underlying caffeine-mediated splicing of this gene subset, we have focused on the induction of SpKLF6, an alternatively spliced form of the KLF6 tumor suppressor gene. Minigene deletion and mutation analyses identified both positive- and negative-acting elements regulating SpKLF6 induction. Notably, the activating element (CafRE) was able to support caffeine-mediated exon inclusion in a heterologous context, and sequences with high homology to the KLF6 CafRE were found in the regions flanking the caffeine-regulated exons of other genes in the subset. Using an RNA interference approach, we have further shown that the SR splicing factor SC35 is required for caffeine-mediated splicing. Strikingly, caffeine treatment led to an increase in SC35 levels commensurate with the induction of SpKLF6, and overexpression of SC35 mimicked the induction of SpKLF6 by caffeine. Furthermore, an RNA EMSA demonstrated increased complex formation at the CafRE with extracts from cells overexpressing SC35, validating a role for this factor in this splicing event. Taken together, our data inform a model of caffeine-regulated alternative splicing of KLF6 and, by inference, a subset of alternatively spliced genes. Caffeine, alternative splicing of KLF6, and SC35. KLF6, a member of the Kruppel-like family of transcription factors, has recently been identified as a key tumor suppressor gene in a variety of cancers, where it is often inactivated by allelic loss and/or somatic mutation (19, 20, 26, 40, 45). Notably, it has been shown that mice with targeted deletion of one KLF6 allele (KLF6+/–) display increased liver mass, suggesting that haploinsufficiency at the KLF6 locus may be sufficient to deregulate cell growth control (41). In the present study, we have observed a caffeine-mediated induction of the SpKLF6 variant at the expense of the wild-type KLF6 transcript, effectively reducing the KLF6 transcript by about half (Fig. 1A). Moreover, while the function of the SpKLF6 splice variant is still under investigation, it is interesting that its novel exon1a includes an early termination codon, predicting an SpKLF6 protein product lacking the DNA binding domain. In all likelihood, this product would either lack transcription activator function or act as a dominant negative regulator. In either case, the resultant haploinsufficiency or loss of KLF6 function, respectively, could very well impact tumor cell growth.
Given the clear role of KLF6 in cancer etiology, and the probability that SpKLF6 will impact that role, KLF6 was chosen as the prototype to investigate the regulation of alternative splicing of cancer-related variants. We have identified two classes of cis elements that regulate expression of the SpKLF6 variant. One class functions as silencers, since mutation of these elements led to constitutive overexpression of SpKLF6. The other element, which we have termed the CafRE, acts as an enhancer mediating caffeine-regulated induction of SpKLF6 (Fig. 6A). In light of these results, we propose a model in which the inclusion of exon 1a is normally repressed through protein interactions at the silencer element(s). Upon caffeine treatment, this repression is relieved by a quantitative or qualitative change in protein interactions at the CafRE. This appears to be a relatively common mechanism for the regulation of alternative splicing; for example, there have been several reports in which splicing activator proteins antagonize the silencing activity of repressor proteins, resulting in the inclusion of relatively weak or cryptic exons (56, 58). In the majority of these cases, the choice between activation and silencing of exon inclusion is most often determined by the relative abundance of the activator and silencer proteins, although posttranslational modifications may also play a role (2, 15).
Our observation that the SR protein SC35 is also induced by caffeine (Fig. 7B) is unique and is consistent with our model. SC35 is a prototypical SR protein initially identified as an essential splicing factor involved in multiple steps of spliceosome assembly (12) and thereafter shown to affect alternative splicing as well (13, 28). Regulation of alternative splicing by SC35 is both complex and context dependent. Like other SR proteins, the binding of SC35 to enhancer elements generally promotes exon splicing, often by counteracting the activity of negatively acting hnRNPs. Our observations that caffeine induces SC35 (Fig. 7B), SC35 overexpression leads to CafRE-dependent cryptic exon inclusion in the context of the KLF6 gene as well as in a heterologous globin gene construct (Fig. 8), and overexpression of SC35 results in increased complex formation at the CafRE in vitro (Fig. 8D) allow us to extend our model to state that caffeine regulates the CafRE-dependent induction of SpKLF6, at least in part, by inducing the expression of SC35, thereby altering the activator/repressor ratio. While computational analysis did not identify an SC35 consensus binding site (29, 49) within the CafRE, alternative splicing regulatory elements are known to be highly degenerative. Therefore, we are currently investigating whether SC35 acts directly by binding to a nonconsensus site within the CafRE or indirectly by modulating or interacting with a CafRE binding protein complex. In either case, we propose that the CafRE complex antagonizes the silencing elements. In this regard, preliminary studies have identified two hnRNPs that are involved in exon 1a repression (data not shown). Future studies will be directed at determining whether these repressor proteins interact with the KLF6 silencing elements and whether their activity is impacted by caffeine.
Caffeine regulates alternative splicing in a subset of genes. Our data from individual gene analysis (KLF6) and splice variant array analysis/validation have clearly demonstrated that caffeine can regulate alternative splicing in a subset of genes. Whether caffeine accomplishes this through activation of SC35 alone or whether caffeine has additional effects on other factors remains to be determined. However, as a first step toward identifying a common mechanism for regulation of this gene subset, we have examined a number of caffeine-regulated genes for the presence of sequences with high homology to the enhancing motif identified within the KLF6 CafRE. Given the high degeneracy of splicing factor binding sites in general, we were encouraged to find matches to the sequence motif within close proximity to the caffeine-regulated exons, although more extensive analyses will be needed to determine the functional relevance of these putative CafREs.
Alternative splicing, SC35 expression, and cancer. Within the past several years, considerable interest has been paid to the role of aberrant alternative splicing in cancer (52). Isoform imbalance is a common occurrence in tumors, and several splice variants are serving as cancer biomarkers due to their correlation with malignancy (47); cancer-specific splice variants may also provide targets for tumor-targeted toxin-conjugated antibodies (24). It has recently been shown that alternative splicing produces oncogenic variants of KLF6 (SV1, SV2, and SV3); these variants are associated with a germ line single-nucleotide polymorphism downstream of exon 1a and have been found to be overexpressed in prostate cancers (38, 39). While these KLF6 splice variants do not appear to be regulated by caffeine (data not shown), these studies do support the notion that KLF6 splice variants can impact the cancer phenotype and mark the introns surrounding exon1a as critical sites for alternative splicing regulation.
Factors involved in alternative splicing pathways are also under investigation as potential therapeutic targets (51). Indeed, SF2/ASF has been recently identified as a proto-oncogene that activates several key targets contributing to particular transformed phenotypes (23). Aberrant increases in SC35 have also been associated with the cancer phenotype (27). For instance, transformation of mouse fibroblasts with either large T antigen or c-erbB2 caused an increase in SC35 levels when compared to untransformed cells (30), and a marked induction of SC35 has been observed in malignant ovarian tissues relative to normal ovarian tissue in certain malignancies (11). It is intriguing to speculate that SC35 is also a proto-oncogene whose transforming potential lies in its ability to regulate some of the genes identified in this study.
Conclusion. In this report we show, for the first time, that the splicing of a subset of genes is regulated by caffeine and identify SC35 as both a mediator of caffeine-regulated splicing and a target of caffeine action. We hypothesize that caffeine is mimicking an endogenous pathway that is aberrantly activated in certain cancer cells to induce alternatively expressed splicing variants of a subset of genes. Thus, our studies have identified both a unique tool and a unique cadre of genes to investigate the mechanisms underlying alternative splicing associated with the cancer phenotype. Given the dearth of information regarding the regulation of alternative splicing, particularly for the many genes whose alternative splice variants have been implicated in the cancer phenotype, we believe that caffeine, and the subset of genes that it targets, provides us with a novel and powerful system by which to gain essential insights into this complicated and potentially clinically relevant process.
| ACKNOWLEDGMENTS |
|---|
This work was supported by NIH-CA06927 (FCCC), NIH-P30-CA072720 (Cancer Institute of New Jersey), NCI-R01 CA 122573 (K.W.S.), and the UMDNJ Foundation grant program (J.S.).
| FOOTNOTES |
|---|
Published ahead of print on 19 November 2007. ![]()
| REFERENCES |
|---|
|
|
|---|
2. Caceres, J. F., T. Misteli, G. R. Screaton, D. L. Spector, and A. R. Krainer. 1997. Role of the modular domains of SR proteins in subnuclear localization and alternative splicing specificity. J. Cell Biol. 138:225-238.
3. Caputi, M., and A. M. Zahler. 2002. SR proteins and hnRNP H regulate the splicing of the HIV-1 tev-specific exon 6D. EMBO J. 21:845-855.[CrossRef][Medline]
4. Carter, M. S., J. Doskow, P. Morris, S. Li, R. P. Nhim, S. Sandstedt, and M. F. Wilkinson. 1995. A regulatory mechanism that detects premature nonsense codons in T-cell receptor transcripts in vivo is reversed by protein synthesis inhibitors in vitro. J. Biol. Chem. 270:28995-29003.
5. Cho, Y. G., B. J. Choi, J. W. Song, S. Y. Kim, S. W. Nam, S. H. Lee, N. J. Yoo, J. Y. Lee, and W. S. Park. 2006. Aberrant expression of krÜppel-like factor 6 protein in colorectal cancers. World J. Gastroenterol. 12:2250-2253.[Medline]
6. Chou, M. Y., N. Rooke, C. W. Turck, and D. L. Black. 1999. hnRNP H is a component of a splicing enhancer complex that activates a c-src alternative exon in neuronal cells. Mol. Cell. Biol. 19:69-77.
7. Dominski, Z., and R. Kole. 1991. Selection of splice sites in pre-mRNAs with short internal exons. Mol. Cell. Biol. 11:6075-6083.
8. Ebbinghaus, C., J. Scheuermann, D. Neri, and G. Elia. 2004. Diagnostic and therapeutic applications of recombinant antibodies: targeting the extra-domain B of fibronectin, a marker of tumor angiogenesis. Curr. Pharm. Des. 10:1537-1549.[CrossRef][Medline]
9. Eperon, I. C., O. V. Makarova, A. Mayeda, S. H. Munroe, J. F. Caceres, D. G. Hayward, and A. R. Krainer. 2000. Selection of alternative 5' splice sites: role of U1 snRNP and models for the antagonistic effects of SF2/ASF and hnRNP A1. Mol. Cell. Biol. 20:8303-8318.
10. Fairbrother, W. G., and L. A. Chasin. 2000. Human genomic sequences that inhibit splicing. Mol. Cell. Biol. 20:6816-6825.
11. Fischer, D. C., K. Noack, I. B. Runnebaum, D. O. Watermann, D. G. Kieback, S. Stamm, and E. Stickeler. 2004. Expression of splicing factors in human ovarian cancer. Oncol. Rep. 11:1085-1090.[Medline]
12. Fu, X. D., and T. Maniatis. 1990. Factor required for mammalian spliceosome assembly is localized to discrete regions in the nucleus. Nature 343:437-441.[CrossRef][Medline]
13. Fu, X. D., A. Mayeda, T. Maniatis, and A. R. Krainer. 1992. General splicing factors SF2 and SC35 have equivalent activities in vitro, and both affect alternative 5' and 3' splice site selection. Proc. Natl. Acad. Sci. USA 89:11224-11228.
14. Graveley, B. R. 2000. Sorting out the complexity of SR protein functions. RNA 6:1197-1211.[CrossRef][Medline]
15. Hanamura, A., J. F. Caceres, A. Mayeda, B. R. Franza, Jr., and A. R. Krainer. 1998. Regulated tissue-specific expression of antagonistic pre-mRNA splicing factors. RNA 4:430-444.[Abstract]
16. Hui, J., L. H. Hung, M. Heiner, S. Schreiner, N. Neumuller, G. Reither, S. A. Haas, and A. Bindereif. 2005. Intronic CA-repeat and CA-rich elements: a new class of regulators of mammalian alternative splicing. EMBO J. 24:1988-1998.[CrossRef][Medline]
17. Hui, L., X. Zhang, X. Wu, Z. Lin, Q. Wang, Y. Li, and G. Hu. 2004. Identification of alternatively spliced mRNA variants related to cancers by genome-wide ESTs alignment. Oncogene 23:3013-3023.[CrossRef][Medline]
18. Ibrahim E. C., T. D. Schaal, K. J. Hertel, R. Reed, and T. Maniatis. 2005. Serine/arginine-rich protein-dependent suppression of exon skipping by exonic splicing enhancers. Proc. Natl. Acad. Sci. USA 102:5002-5007.
19. Ito, G., M. Uchiyama, M. Kondo, S. Mori, N. Usami, O. Maeda, T. Kawabe, Y. Hasegawa, K. Shimokata, and Y. Sekido. 2004. Kruppel-like factor 6 is frequently down-regulated and induces apoptosis in non-small cell lung cancer cells. Cancer Res. 64:3838-3843.
20. Jeng, Y. M., and H. C. Hsu. 2003. KLF6, a putative tumor suppressor gene, is mutated in astrocytic gliomas. Int. J. Cancer 105:625-629.[CrossRef][Medline]
21. Jin, S., and K. W. Scotto. 1998. Transcriptional regulation of the MDR1 gene by histone acetyltransferase and deacetylase is mediated by NF-Y. Mol. Cell. Biol. 18:4377-4384.
22. Kalnina, Z., P. Zayakin, K. Silina, and A. Line. 2005. Alterations of pre-mRNA splicing in cancer. Genes Chromosomes Cancer 42:342-357.[CrossRef][Medline]
23. Karni, R., E. de Stanchina, S. W. Lowe, R. Sinha, D. Mu, and A. R. Krainer. 2007. The gene encoding the splicing factor SF2/ASF is a proto-oncogene. Nat. Struct. Mol. Biol. 14:185-193.[CrossRef][Medline]
24. Kessler, K. R., and R. Benecke. 1997. The EBD test: a clinical test for the detection of antibodies to botulinum toxin type A. Mov. Disord. 12:95-99.[CrossRef][Medline]
25. Kim, S., H. Shi, D. K. Lee, and J. T. Lis. 2003. Specific SR protein-dependent splicing substrates identified through genomic SELEX. Nucleic Acids Res. 31:1955-1961.
26. Kimmelman, A. C., R. F. Qiao, G. Narla, A. Banno, N. Lau, P. D. Bos, N. Nunez Rodriguez, B. C. Liang, A. Guha, J. A. Martignetti, S. L. Friedman, and A. M. Chan. 2004. Suppression of glioblastoma tumorigenicity by the Kruppel-like transcription factor KLF6. Oncogene 23:5077-5083.[CrossRef][Medline]
27. Kirschbaum-Slager, N., G. M. Lopes, P. A. Galante, G. J. Riggins, and S. J. de Souza. 2004. Splicing factors are differentially expressed in tumors. Genet. Mol. Res. 3:512-520.[Medline]
28. Lin, S., R. Xiao, P. Sun, X. Xu, and X. D. Fu. 2005. Dephosphorylation-dependent sorting of SR splicing factors during mRNP maturation. Mol. Cell 20:413-425.[CrossRef][Medline]
29. Liu, H. X., S. L. Chew, L. Cartegni, M. Q. Zhang, and A. R. Krainer. 2000. Exonic splicing enhancer motif recognized by human SC35 under splicing conditions. Mol. Cell. Biol. 20:1063-1071.
30. Maeda, T., and S. Furukawa. 2001. Transformation-associated changes in gene expression of alternative splicing regulatory factors in mouse fibroblast cells. Oncol. Rep. 8:563-566.[Medline]
31. Manley, J. L., and R. Tacke. 1996. SR proteins and splicing control. Genes Dev. 10:1569-1579.
32. Matlin, A. J., F. Clark, and C. W. Smith. 2005. Understanding alternative splicing: towards a cellular code. Nat. Rev. Mol. Cell Biol. 6:386-398.[CrossRef][Medline]
33. Mayeda, A., S. H. Munroe, J. F. Caceres, and A. R. Krainer. 1994. Function of conserved domains of hnRNP A1 and other hnRNP A/B proteins. EMBO J. 13:5483-5495.[Medline]
34. Modafferi, E. F., and D. L. Black. 1997. A complex intronic splicing enhancer from the c-src pre-mRNA activates inclusion of a heterologous exon. Mol. Cell. Biol. 17:6537-6545.[Abstract]
35. Modrek, B., and C. J. Lee. 2003. Alternative splicing in the human, mouse and rat genomes is associated with an increased frequency of exon creation and/or loss. Nat. Genet. 34:177-180.[CrossRef][Medline]
36. Modrek, B., A. Resch, C. Grasso, and C. Lee. 2001. Genome-wide detection of alternative splicing in expressed sequences of human genes. Nucleic Acids Res. 29:2850-2859.
37. Nagy, E., and L. E. Maquat. 1998. A rule for termination-codon position within intron-containing genes: when nonsense affects RNA abundance. Trends Biochem. Sci. 23:198-199.[CrossRef][Medline]
38. Narla, G., A. Difeo, H. L. Reeves, D. J. Schaid, J. Hirshfeld, E. Hod, A. Katz, W. B. Isaacs, S. Hebbring, A. Komiya, S. K. McDonnell, K. E. Wiley, S. J. Jacobsen, S. D. Isaacs, P. C. Walsh, S. L. Zheng, B. L. Chang, D. M. Friedrichsen, J. L. Stanford, E. A. Ostrander, A. M. Chinnaiyan, M. A. Rubin, J. Xu, S. N. Thibodeau, S. L. Friedman, and J. A. Martignetti. 2005. A germline DNA polymorphism enhances alternative splicing of the KLF6 tumor suppressor gene and is associated with increased prostate cancer risk. Cancer Res. 65:1213-1222.
39. Narla, G., A. DiFeo, S. Yao, A. Banno, E. Hod, H. L. Reeves, R. F. Qiao, O. Camacho-Vanegas, A. Levine, A. Kirschenbaum, A. M. Chan, S. L. Friedman, and J. A. Martignetti. 2005. Targeted inhibition of the KLF6 splice variant, KLF6 SV1, suppresses prostate cancer cell growth and spread. Cancer Res. 65:5761-5768.
40. Narla, G., K. E. Heath, H. L. Reeves, D. Li, L. E. Giono, A. C. Kimmelman, M. J. Glucksman, J. Narla, F. J. Eng, A. M. Chan, A. C. Ferrari, J. A. Martignetti, and S. L. Friedman. 2001. KLF6, a candidate tumor suppressor gene mutated in prostate cancer. Science 294:2563-2566.
41. Narla, G., S. Kremer-Tal, N. Matsumoto, X. Zhao, S. Yao, K. Kelley, M. Tarocchi, and S. L. Friedman. 2007. In vivo regulation of p21 by the Kruppel-like factor 6 tumor-suppressor gene in mouse liver and human hepatocellular carcinoma. Oncogene 26:4428-4434.[CrossRef][Medline]
42. Noonan, K. E., C. Beck, T. A. Holzmayer, J. E. Chin, J. S. Wunder, I. L. Andrulis, A. F. Gazdar, C. L. Willman, B. Griffith, D. D. Von Hoff, et al. 1990. Quantitative analysis of MDR1 (multidrug resistance) gene expression in human tumors by polymerase chain reaction. Proc. Natl. Acad. Sci. USA 87:7160-7164.
43. Ohler, U., N. Shomron, and C. B. Burge. 2005. Recognition of unknown conserved alternatively spliced exons. PLoS Comput. Biol. 1:113-122.[Medline]
44. Pan, X. C., Z. Chen, F. Chen, X. H. Chen, C. Zhou, and Z. G. Yang. 2006. Mutations of the tumor suppressor Kruppel-like factor 6 (KLF6) gene in hepatocellular carcinoma and its effect of growth suppression on human hepatocellular carcinoma cell line HepG2. Zhonghua Gan Zang Bing Za Zhi 14:109-113. (In Chinese.)[Medline]
45. Reeves, H. L., G. Narla, O. Ogunbiyi, A. I. Haq, A. Katz, S. Benzeno, E. Hod, N. Harpaz, S. Goldberg, S. Tal-Kremer, F. J. Eng, M. J. Arthur, J. A. Martignetti, and S. L. Friedman. 2004. Kruppel-like factor 6 (KLF6) is a tumor-suppressor gene frequently inactivated in colorectal cancer. Gastroenterology 126:1090-1103.[CrossRef][Medline]
46. Sadusky, T., A. J. Newman, and N. J. Dibb. 2004. Exon junction sequences as cryptic splice sites: implications for intron origin. Curr. Biol. 14:505-509.[Medline]
47. Smid-Koopman, E., L. J. Blok, T. J. Helmerhorst, S. Chadha-Ajwani, C. W. Burger, A. O. Brinkmann, and F. J. Huikeshoven. 2004. Gene expression profiling in human endometrial cancer tissue samples: utility and diagnostic value. Gynecol. Oncol. 93:292-300.[CrossRef][Medline]