Alternative splicing of apoptosis-related genes in imatinib-treated K562 cells identified by exon array analysis

  • Authors:
    • Jing Liu
    • Yun Xiao
    • Huo-Mei Xiong
    • Jing Li
    • Bo Huang
    • Hai-Bin Zhang
    • Dan-Qin Feng
    • Xi-Min Chen
    • Xiao-Zhong Wang
  • View Affiliations

  • Published online on: December 29, 2011     https://doi.org/10.3892/ijmm.2011.872
  • Pages: 690-698
Metrics: Total Views: 0 (Spandidos Publications: | PMC Statistics: )
Total PDF Downloads: 0 (Spandidos Publications: | PMC Statistics: )


Abstract

Imatinib is the therapeutic standard for newly diagnosed patients with chronic myeloid leukemia (CML). In these patients, imatinib has been shown to induce an apoptotic response specifically in cells expressing the oncogenic fusion protein BCR-ABL. Previous studies in our lab revealed that imatinib-induced apoptosis in K562 cells involves a shift in production of Bcl-x splice isoforms towards the pro-apoptotic Bcl-xs splice variant. Here, we report the findings from our subsequent study to identify other apoptosis-related genes that are differentially spliced in response to imatinib treatment. Gene expression profiling of imatinib-treated K562 cells was performed by the Affymetrix GeneChip® Human Exon 1.0 ST array, and differences in exon-level expression and alternative splicing were analyzed using the easyExon software. Detailed analysis by reverse transcription-PCR (RT-PCR) and sequencing of key genes confirmed the experimental results of the exon array. Our results suggest that imatinib treatment of K562 cells causes a transcriptional shift towards alternative splicing in a large number of apoptotic genes. The present study provides insight into the molecular character of apoptotic leukemia cells and may help to improve the mechanism of imatinib therapy in patients with CML.

Introduction

Chronic myeloid leukemia (CML) is a hematopoietic stem cell disorder characterized by marked increase in granulocytes, bone marrow hyperplasia, and splenomegaly. The current annual incidence of CML has been reported to be 1 to 2 cases per 100,000 people, and accounts for 15% of adult leukemias diagnosed worldwide (1). The molecular hallmark of CML is the Philadelphia (Ph) chromosome, a reciprocal translocation t (9:22) (q34; q11) between the proto-oncogene ABL on chromosome 9 and the BCR gene on chromosome 22. This abnormality has been clinically detected in 95% of CML patients, and is considered the primary risk factor for CML development (2). The Ph translocation creates a fusion BCR-ABL oncogene, which encodes a functional protein with constitutive tyrosine kinase activity that aberrantly modulates tumorigenic pathways.

Imatinib is the first chemotherapeutic agent described as a specific BCR-ABL tyrosine kinase inhibitor. Its proven efficacy for the treatment of patients with CML has led to its preferred use as first-line treatment, supplanting the traditional single agents busulphan, hydroxyurea and interferon-α (3). The molecular mechanism of imatinib has been described: imatinib potently binds to the inactive form of BCR-ABL, thereby blocking the ATP-binding site and its enzymatic activity (4). Imatinib treatment of CML patients has produced estimated cumulative best rates of complete hematologic response (CHR) and complete cytogenetic response (CCyR) of 98 and 87%, respectively (5). Despite the significant improvements in the outcome for patients with CML following treatment with imatinib, the exact mechanism of its anti-leukemic effect remains to be completely elucidated (6).

In a previous study, we discovered that imatinib treatment of K562 leukemia cells produced a significant change in the ratio of Bcl-x splice variants. The Bcl-x gene undergoes alternative splicing at exon 2 to produce two isoforms with distinctive anti- and pro-apoptotic properties. The Bcl-xS isoform promotes apoptosis, while the Bcl-xL inhibits apoptosis. Imatinib treatment of K562 cells led to a transcriptional shift towards the Bcl-xS isoform, which in turn promoted apoptosis in the treated cells. Accordingly, this finding suggested that imatinib could regulate the alternative splicing of other apoptotic genes in K562 cells.

Alternative splicing is a common mechanism used by the human system to generate diversity in the transcriptome and proteome. Through alternative splicing, multiple different transcripts can be generated from a single mRNA precursor. Recent studies have suggested that up to 94% of all human genes undergo alternative splicing, and aberrant pre-mRNA splicing has been implicated in the pathogenesis of neoplasia (7,8). Likewise, studies have demonstrated that most of the currently used chemotherapeutic agents not only induce an apoptotic response, but also affect the alternative splicing of some genes. Shkreta et al (9) determined that 20 of the mainstream anticancer drugs could influence the production of Bcl-x splice isoforms. In addition, these drugs could alter a subset of alternative splicing events in some cell lines, suggesting that the influence of anticancer drugs on alternative splicing is a very common phenomenon.

In recent years, powerful techniques for genome-wide identification and analysis of alternative splicing isoforms have been developed. These large-scale high-throughput analytical methods have been successfully applied for the identification of differential splicing events in various cancer tissues (10). Exon arrays, which contain both known and predicted exons, are rapidly gaining popularity and becoming a standard for both gene- and exon-level expression analyses (11). The Affymetrix GeneChip® Human Exon 1.0 ST array (Human Exon 1.0 ST array) is the latest product in the family of exon arrays. It contains approximately 5.4 million probes grouped into 1.4 million probe sets, allowing the interrogation of over 1 million exon clusters, which are exon annotations from various sources that overlap by genomic location. For each gene, the average number of probes on the exon array is 35, and these probes are usually distributed along the entire transcript sequence. The commercial availability of specific software for exon-level expression analysis has provided a powerful tool for the study of alternative splicing for every known and predicted exon. At the same time, exon arrays provide robust gene-level expression analysis (Affymetrix Technical notes, Identifying and Validating Alternative Splicing Events) facilitating identification of relatively subtle differences in the isoform expression ratios.

In our current study, we employed the Human Exon 1.0 ST Array to investigate the comprehensive profile of apoptotic genes in leukemia K562 cells that are differentially spliced in response to imatinib treatment.

Materials and methods

Cell culture

Human myeloid leukemia cells K562 were grown in Roswell Park Memorial Institute 1640 (RPMI-1640; Hyclone, USA) media supplemented with 10% fetal bovine serum (FBS; Hyclone) in a 5% CO2 humidified atmosphere at 37°C. Experimental cultures were initiated at a density of 2×105 cells/ml. Imatinib (Novartis Co., China) was added from a 1 mM stock solution in DMSO to achieve various final concentrations in culture. The control non-treated cells were generated by adding an equal volume of DMSO alone. Cells were incubated for up to 24 h before cell proliferation was assessed by WST-1 reagent based assay kit (Millipore, USA) and apoptosis was assessed by flow cytometry. For details of the two methods, refer to our previous study (12).

RNA extraction and array hybridization

Total-RNA was extracted from 107 K562 cells by using the TRIzol reagent (Tiangen, China) following the manufacturer’s instructions. The extracted RNA was further purified by DNase I treatment using an RNeasy mini kit and an RNase-Free DNase Set (Qiagen, China), according to the manufacturer’s protocols. The concentration and purity of total-RNA were determined by spectrophotometric measurements at an optical density (OD) of 260/280 nm and electrophoresis on 1% agarose gel. An OD260/OD280 ratio of 1.8–2.0 was considered to indicate pure RNA. RNA samples were immediately frozen and stored at −80°C until required for use in further microarray studies.

Two 1-μg RNA samples from untreated control K562 cells and imatinib-treated K562 cells (1 μM of imatinib for 24 h) were labeled according to the GeneChip whole transcript (WT) Sense Target Labeling Assay as provided by the manufacturer (http://www.affymetrix.com). Labeled RNA was hybridized to the Human Exon 1.0 ST Arrays for 16 h, following the Affymetrix protocol. Afterwards, arrays were scanned using the Affymetrix GCS 3000 7G Scanner and digital images were processed by the GeneChip® Operating Software v1.4 to produce raw CEL intensity files for analysis of differential expression.

Data analysis

Data from the exon arrays were processed using the easyExon software (http://microarray.ym.edu.tw/easyexon). This tool provides a user-friendly, platform-independent and efficient processing method to explore Affymetrix exon array data. Two parallel analyses (gene-level and exon-level) were carried out by this software. To improve the processing efficiency, we downloaded the related human gene database and meta probeset files from the software homepage (http://microarray.ym.edu.tw:8080/easyexon/index.jsp?mode=support) and installed each as a local copy.

After data loading of the CEL files, gene/exon signal summarizations were calculated automatically by the Affymetrix Power Tools (APT) package, which can execute the binary file ‘apt-probeset-summarize’ to generate both exon-level and gene-level summary files. Two of the most commonly used signal estimation algorithms, RMA and PLIER, were implemented to combine information from probes belonging to the same transcript, or exon, to generate a summarized expression signal value of the single represented gene or exon (13). Background noise was detected using the ‘Detection above Background (DABG)’ algorithm. Normalization was performed by the ‘quantile normalization’ algorithm for both the exon- and gene-level.

Reverse transcription-PCR (RT-PCR) and sequencing validation

RT-PCR and direct sequencing were used to confirm the differentially expressed alternatively spliced genes detected by the Human Exon 1.0 ST Arrays. Reverse transcription was performed with 1 μg of total-RNA using a reverse transcription kit (Takara Bio, Inc., Japan), following the manufacturer’s instructions. For subsequent PCR analysis, 1 μl of cDNA was used. The primers and reaction conditions for each transcript are shown in Table I. The amplified products were electrophoresed through 1.5–2.5% agarose gels, according to amplified fragment size requirements for resolution. Images of the gels were analyzed by the Quantity One 4.6.2 software (Bio-Rad, USA). Ratios between splice variants expressed in untreated and imatinib-treated cells were calculated. Distinct gel bands were purified with a gel extraction kit (Tiangen), sequenced using the same primers as for PCR and the BigDye Terminator v3.1 Cycle Sequencing kit (Applied Biosystems, USA), and analyzed on an ABI PRISM 3100 Genetic Analyzer (Applied Biosystems).

Table I

Primers and amplification conditions used in this study.

Table I

Primers and amplification conditions used in this study.

Gene nameSequence (5′→3′)Reaction conditionSplicing isoformProduct size (bp)
Bcl-xF: ATGGCAGCAGTAAAGCAAGCG32 cycles at 94°C for 30 sec, 55°C for 30 sec, and 72°C for 1 min Bcl-xL456
R: TCATTTCCGACTGAAGAGTGA Bcl-xS267
RNF34F: TGAGTTTCCTGGTAGAGCC35 cycles at 94°C for 30 sec, 55°C for 30 sec, and 72°C for 1 minVariant L173
R: TTCCCATGACTTCATTCAGCVariant S107
SNCAF: GTTGGAGGAGCAGTGGTG35 cycles at 94°C for 30 sec, 64.5°C for 45 sec, and 72°C for 1 minVariant L313
R: ATGACTGGGCACATTGGAVariant S229
β-actinF: CGGGAAATCGTGCGTGAC35 cycles at 94°C for 30 sec, 55°C for 30 sec, and 72°C for 1 min-443
R: TGGAAGGTGGACAGCGAGG

Results

Exon-level expression differences induced by imatinib treatment

The pro-apoptotic activity of imatinib on K562 cells was first confirmed by a WST-1 assay and a double fluorescence assay to analyze cell proliferation and apoptosis, respectively. Our results revealed that imatinib exposure rapidly induced proliferation inhibition and apoptosis of K562 cells in a time-dependent manner [data not shown; identical to that of our previous study (12)].

The expression of alternatively spliced genes in the K562 cells in response to treatment with imatinib were then assessed via exon array and easyExon analysis. Before performing statistical filtration, the limitation number of probesets in a transcript cluster was set from 4 to 40. For any transcript cluster with probesets greater than the setting number, only the setting number of probesets was included in the latter analysis. Three statistical methods for exon-level filtration in the easyExon software were available: Affymetrix MIDAS (Microarray Detection of Alternative Splicing), Partek AS ANOVA (Alternative Splice Analysis of Variance), and PAC (Pattern-Based Correlation) for exon analysis. We used the default Affymetrix MIDAS. For selecting alternatively spliced genes, the criteria for probeset filtration were set as follows, fold-change (FC) ≥2 or ≤0.5, with a statistical significance of P<0.05.

A total of 331 genes represented by one or several probesets met the above criteria. Most of these probesets presented change in expression of less than 3-fold. The genes affiliated with these probesets were investigated by the easyExon transcript cluster filtration tool to determine which were related to cell apoptosis. Briefly, the easyExon ‘GO Biological Process’ parameter was set for ‘apoptosis’ related to transcript-affiliated gene name or accession number. The program identified 175 apoptosis-related genes from the original set of 331. Each of these genes contained at least one probeset meeting the above criteria, for a total of 430 probesets. Tables II and III list the 20 most upregulated and most downregulated probesets.

Table II

Alternative splice probesets that were significantly upregulated in K562 cells in response to imatinib treatment (20 most differentially expressed)a.

Table II

Alternative splice probesets that were significantly upregulated in K562 cells in response to imatinib treatment (20 most differentially expressed)a.

Gene probe IDExon probe IDSplice indexGene titleGene symbolAccession no.
4007617400763415.12Pim-2 oncogene Pimatinib2NM_006875
263867626386929.74ELL associated factor 2EAF2NM_018456
393954539395499.47Macrophage migration inhibitory factor (glycosylation-inhibiting factor)MIFNM_002415
257475225747858.85Excision repair cross-complementing rodent repair deficiency, complementation group 3 (xeroderma pigmentosum group B complementing)ERCC3NM_000122
381846838184807.39Tumor necrosis factor (ligand) superfamily, member 9TNFSF9NM_003811
293903429390566.77Serpin peptidase inhibitor, clade B (ovalbumin), member 6 SERPINB6NM_004568
312770331277316.05Tumor necrosis factor receptor superfamily, member 10b TNFRSF10BNM_003842
373629037362926.00Baculoviral IAP repeat-containing 5 (survivin)BIRC5NM_001168
252474325247695.58FAST kinase domains 2FASTKD2NM_014929
381650938165155.30Growth arrest and DNA-damage-inducible, βGADD45BNM_015675
345466234546755.05Family with sequence similarity 130, member A1 FAM130A1NM_030809
386818338682424.80Interleukin 4 induced 1IL4I1NM_152899
379070437907074.62 Phorbol-12-myristate-13-acetate-induced protein 1PMAIP1NM_021127
395643339564834.58CHK2 checkpoint homolog (S. pombe)CHEK2NM_001005735
294092029409594.50Eukaryotic translation elongation factor 1 epsilon 1EEF1E1NM_004280
388670438867114.50Serine/threonine kinase 4STK4NM_006282
342625734262614.47Suppressor of cytokine signaling 2SOCS2NM_003877
263867626386904.36ELL associated factor 2EAF2NM_018456
294092029409624.29Eukaryotic translation elongation factor 1 epsilon 1EEF1E1NM_004280
258921425892154.27Protein kinase, interferon-inducible double stranded RNA dependent activatorPRKRANM_003690

a A more comprehensive list of alternatively spliced apoptotic genes in K562 cells treated with imatinib is available upon request.

Table III

Alternative splice probesets that were significantly down-regulated in K562 cells in response to imatinib treatment (20 most differentially expressed).

Table III

Alternative splice probesets that were significantly down-regulated in K562 cells in response to imatinib treatment (20 most differentially expressed).

Gene probe IDExon probe IDSplice indexGene titleGene symbolAccession no.
374278337428000.21NLR family, pyrin domain containing 1NLRP1NM_033004
348141034814440.21Tumor necrosis factor receptor superfamily, member 19 TNFRSF19NM_018647
342625734262790.24Suppressor of cytokine signaling 2SOCS2NM_003877
387316038731770.26Tribbles homolog 3 (Drosophila)TRIB3NM_021158
386530138653280.27Excision repair cross-complementing rodent repair deficiency, complementation group 2 (xeroderma pigmentosum D)ERCC2NM_000400
374278337428080.27NLR family, pyrin domain containing 1NLRP1NM_033004
387316038731790.29Tribbles homolog 3 (Drosophila)TRIB3NM_021158
342625734262800.3Suppressor of cytokine signaling 2SOCS2NM_003877
348141034814590.3Tumor necrosis factor receptor superfamily, member 19 TNFRSF19NM_018647
387316038731780.3Tribbles homolog 3 (Drosophila)TRIB3NM_021158
383879538388060.31Bcl-2-like 12 (proline rich)BCL2L12NM_138639
369950836995570.31Craniofacial development protein 1CFDP1NM_006324
392414439241800.31Collagen, type XVIII, α 1COL18A1NM_030582
387316038731760.31Tribbles homolog 3 (Drosophila)TRIB3NM_021158
374278337427920.32NLR family, pyrin domain containing 1NLRP1NM_033004
346540934654130.33B-cell translocation gene 1, anti-proliferativeBTG1NM_001731
344184934418610.33Tumor necrosis factor receptor superfamily, member 1A TNFRSF1ANM_001065
344184934418630.33Tumor necrosis factor receptor superfamily, member 1A TNFRSF1ANM_001065
319489631949520.33TNF receptor-associated factor 2TRAF2NM_021138
248741224874470.34Annexin A4ANXA4NM_001153
Detection of alternative splicing events in response to imatinib treatment

EasyExon provides a module to visualize and interpret selected subsets of exons and transcripts that meet specified criteria. The normalized, log-transformed and variance stabilized probeset intensities of a transcript are plotted as signal mean ± standard error in an intensity plot. Probesets marked with asterisks by the program represent statistically significant alternative splicing events which pass the defined MIDAS/Partek AS ANOVA threshold. Probesets colored in bold by the program indicate cases in which the intensity FC between groups was greater than the threshold (by at least 2-fold). The degree of alternative splicing can also be visualized by the splicing index (SI) plot in the same panel. Different probesets for the same exon group together and are visually indicated by the program using a triangle or trapezoid. Adjacent exons are distinguished by different colors. From the graphic presentation of analysis results, we were able to identify whether a specific isoform was up or downregulated in the imatinib-treated cells (Fig. 1B) (13).

After identifying probesets that were differentially spliced in response to imatinib treatment, we randomly chose three internal exons (RNF34, DAPK1 and SNCA) for experimental validation. In addition, we included Bcl-x as the positive control, which was shown to be differentially spliced in K562 cells treated with imatinib in our previous study (12).

We introduce here a workflow for detecting alternative splicing events between untreated K562 cells and imatinib-treated K562 cells by using the easyExon program. In Bcl-x (ID 3902489), the probeset (ID 3902522, in bold) targeting the 3′ region of exon 2 was found to be downregulated in imatinib-treated K562 cells. The expression fold-change of this probeset was larger than 2 (labeled in black) and the difference was statistically significant (marked by an asterisk by the program). This indicated that there were two Bcl-x isoforms in K562 cells, one with both the 5′ and 3′ region of exon 2 and another with only the 5′ region of exon 2, and the latter one was upregulated in imatinib-treated K562 cells (Fig. 1B). In RNF34 (ID 3434823), the second probeset (ID 3434832, in bold) targeting exon 2 was found to be downregulated in imatinib-treated K562 cells. In SNCA (ID 2777714), the probeset (ID 2777722) targeting exon 5 was found to be upregulated in imatinib-treated K562 cells. Both RNF34 and SNCA had two spliceosomes in K562 cells and the alternative splicing was affected by imatinib (Figs. 2B and 3B). Similarly, we were able to identify other transcripts that undergo alternative splicing by using this process.

RT-PCR and sequencing validation

The expression of Bcl-x, RNF34 and SNCA mRNAs in K562 cells before and after treatment with imatinib were evaluated by RT-PCR using primers targeting the exons that flank the exon of interest. Amplification results were visualized by gel electrophoresis. When two alternatively spliced isoforms were present, two different bands would be produced by the same primer pair between the two groups (Figs. 1A and C, 2A and C, 3A and C). The relative quantity of the transcripts was calculated by the Quantity One software. The results indicated a difference in the expression ratio of the isoforms in response to imatinib treatment for probesets 3902522 (Bcl-x), 3434832 (RNF34) and 2777722 (SNCA). The relative quantity analyses for Bcl-x, RNF34 and SNCA were consistent with the trend of SI values calculated from the exon array data (Figs. 1C, 2C and 3C). To further validate the differences between the spliceosomes, we performed sequencing validation of the distinct bands resolved by gel electrophoresis. The sequencing results for RNF34 and SNCA confirmed the alternative splicing events that were predicted by the exon array data (Figs. 2D and 3D).

For DAPK1, the sequencing validation did not support the presence of alternatively spliced isoforms even though two gel bands that appeared to be in the correct position were detected (data not shown).

Discussion

Alternative splicing of mRNA precursors is a nearly ubiquitous and extremely flexible point of gene control in humans (14). Abnormal alternative splicing has been shown to be associated with many human cancers (15). In addition, previous studies have demonstrated that some anticancer drugs can affect alternative splicing of genes in cancer cells (9,16). These results may provide a novel rationale for the effectiveness of anticancer drugs.

Apoptosis is a major pathway involved in the complex homeostatic balance between cellular proliferation and cell death. Previous studies have also demonstrated that transcripts from a significant number of genes involved in apoptosis are alternatively spliced, often resulting in isoforms with opposing roles in promoting or preventing cell death (1719). Examples of this type are found in every category of apoptotic regulators, from transmembrane receptors (Fas, Fas ligand and LARD) and adaptor molecules (Bcl-x, Bak, Apaf1, survivin, Mcl-1 and TRAF2) to caspases (caspase-1, -2, -6, -7, -8 and-9) and executors (FLIP and ICAD) (20). Although alternative splicing is an important mechanism regulating apoptosis, the effect of chemotherapeutic agents on the alternative splicing of apoptotic genes has rarely been examined (9).

In previous studies, we have demonstrated that one of the most commonly used chemotherapeutic agents, imatinib, which is considered the most effective and a relatively safe drug for the treatment of the chronic phase of CML, could regulate alternative splicing of the apoptotic gene Bcl-x in K562 cells through the activation of PP1 (12). In the present study, we employed the Human Exon 1.0 ST Array, together with the easyExon software, to assess the comprehensive profile of alternatively spliced apoptotic genes in imatinib-treated K562 cells. A total of 175 (175/661, 26.48%) differentially spliced apoptosis-related transcripts were identified, several of which were successfully validated by RT-PCR and sequencing. Our results revealed that imatinib treatment led to aberrant alternative splicing events in K562 cells.

Table IV summarizes the previous studies in the literature that have demonstrated the association of some of these alternatively spliced isoforms with apoptotic events. These genes belong to different cell death protein families, including adaptor proteins, Bcl-2 family, and caspases, in which the splicing mechanism has been shown to modulate both intrinsic (mitochondrial) and extrinsic (death receptor) apoptosis pathways. These two signaling routes to cell death are linked via the Bcl-2 family member Bid, whose alternative splicing was also detected by the exon array in our study. Thus, it is likely that imatinib has wide-ranging effects on apoptosis of K562 cells. In addition, our results suggest the existence of a novel mechanism by which imatinib can induce cell apoptosis, which is the regulation of alternative splicing of a series of apoptotic genes.

Table IV

Apoptotic genes previously shown to be regulated by alternative splicing.

Table IV

Apoptotic genes previously shown to be regulated by alternative splicing.

Gene symbol (ref.)Gene titleAccession no.Gene probe IDExon probe IDSplice index
Adaptor proteins and regulators
TRAF2 (22)TNF receptor-associated factor 2NM_021138319489631949190.50
31949500.50
31949440.48
31949520.33
Apaf-1 (23)Apoptotic peptidase activating factor 1NM_181868342787634278773.68
34279363.00
34279372.46
34278812.26
34278892.08
Bcl-2 family
Bcl-x (24)Bcl-2-like 1NM_001191390248939025230.4
Bak (25) Bcl-2-antagonist/killer 1NM_001188295075329507682.14
Bax (26)Bcl2-associated X proteinNM_004324383806738380840.48
Bid (18)BH3 interacting domain death agonistNM_197966395192739519503.91
39519293.64
39519282.63
39519302.15
Bmf (17)Bcl-2 modifying factorNM_001003940361922936192423.59
Caspases and caspase-like proteins
Caspase-3 (19)Caspase 3, apoptosis-related cysteine peptidaseNM_004346279648427964992.00
27965052.00
Caspase-10 (27)Caspase 10, apoptosis-related cysteine peptidaseNM_032977252269325227192.33
25227072.07
25227152.00

Human Exon 1.0 ST Arrays, which are comprised of probesets that cover the entire length of all known transcripts, facilitate the study of gene expression at the exon level (probeset expression) and predict the likelihood of encountering alternative splicing for a given gene. This platform has also enabled detection of specific changes at the gene expression level (transcript cluster expression). Exon array technology has several key advantages over conventional analytical methods, including high-throughput, updated content, probesets designed to span the entire length of a given transcript, and expression analysis at the exon level. In the present study, using RT-PCR and sequencing, two of the three randomly-chosen exons (67%) from the 430 probesets were able to be validated for cell differences in abundance of the respective transcript isoforms. Though the exon array technology is reliable and sensitive enough to detect differential expression at the exon level, there are some important limitations that exist and should be considered when analyzing the exon array results. The sensitivity and specificity of the Affymetrix exon arrays for detecting exon splicing across the whole genome has not yet been defined, at least based on the current published data. Thus, confirmation of the detected exon splicing events requires PCR with primers that flank the exons of interest or RNA deep sequencing to demonstrate consistent differential splicing (21).

In conclusion, using exon arrays, we have discovered that imatinib can produce a wide-ranging effect on the alternative splicing of apoptotic factors in K562 leukemic cells. This information may help to improve the mechanism of imatinib therapy in patients with CML. Even though several of the imatinib-induced alternative splicing events were successfully validated by RT-PCR and sequencing in our study, a more comprehensive validation will be necessary to provide a more accurate estimation of the current findings.

Acknowledgements

This study was supported by a grant from the National Science Foundation of China (no. 30700338).

References

1 

S FaderlM TalpazZ EstrovS O’BrienR KurzrockHM KantarjianThe biology of chronic myeloid leukemiaN Engl J Med341164172199910.1056/NEJM199907153410306

2 

CL SawyersChronic myeloid leukemiaN Engl J Med34013301340199910.1056/NEJM19990429340170610219069

3 

WY AuPB CaguioaC ChuahChronic myeloid leukemia in AsiaInt J Hematol891423200910.1007/s12185-008-0230-019101781

4 

T SchindlerW BornmannP PellicenaWT MillerB ClarksonJ KuriyanStructural mechanism for STI-571 inhibition of abelson tyrosine kinaseScience28919381942200010.1126/science.289.5486.193810988075

5 

JM GoldmanHow I treat chronic myeloid leukemia in the imatinib eraBlood11028282837200710.1182/blood-2007-04-03894317626839

6 

A PuissantP ColosettiG RobertJP CassutoS RaynaudP AubergerCathepsin B release after imatinib-mediated lysosomal membrane permeabilization triggers BCR-ABL cleavage and elimination of chronic myelogenous leukemia cellsLeukemia24115124201010.1038/leu.2009.233

7 

ET WangR SandbergS LuoAlternative isoform regulation in human tissue transcriptomesNature456470476200810.1038/nature0750918978772

8 

JP VenablesAberrant and alternative splicing in cancerCancer Res6476477654200410.1158/0008-5472.CAN-04-191015520162

9 

L ShkretaU FroehlichER PaquetJ ToutantSA ElelaB ChabotAnticancer drugs affect the alternative splicing of Bcl-x and other human apoptotic genesMol Cancer Ther713981409200810.1158/1535-7163.MCT-08-019218566212

10 

JA CalarcoAL SaltzmanJY IpBJ BlencoweTechnologies for the global discovery and analysis of alternative splicingAdv Exp Med Biol6236484200710.1007/978-0-387-77374-2_518380341

11 

H ChenY GuoM HuW DuanG ChangC LiDifferential expression and alternative splicing of genes in lumbar spinal cord of an amyotrophic lateral sclerosis mouse modelBrain Res13405269201010.1016/j.brainres.2010.03.07520362558

12 

Y XiaoH XiongJ LiImatinib regulates the alternative pre-mRNA splicing of Bcl-x in K562 cellsAsian Biomed(In press)

13 

TY ChangYY LiCH JeneasyExon - a Java-based GUI tool for processing and visualization of Affymetrix exon array dataBMC Bioinformatics9432200810.1186/1471-2105-9-43218851762

14 

CJ DavidJL ManleyAlternative pre-mRNA splicing regulation in cancer: pathways and programs unhingedGenes Dev2423432364201010.1101/gad.197301021041405

15 

Z KalninaP ZayakinK SilinaA LineAlterations of pre-mRNA splicing in cancerGenes Chromosomes Cancer42342357200510.1002/gcc.2015615648050

16 

JG ChangDM YangWH ChangSmall molecule amiloride modulates oncogenic RNA alternative splicing to devitalize human cancer cellsPLoS One6e18643201110.1371/journal.pone.001864321694768

17 

AA MoralesA OlssonF CelsingA OsterborgM JondalLM OsorioExpression and transcriptional regulation of functionally distinct Bmf isoforms in B-chronic lymphocytic leukemia cellsLeukemia184147200410.1038/sj.leu.240318314574334

18 

SA RenshawCE DempseyFA BarnesThree novel Bid proteins generated by alternative splicing of the human Bid geneJ Biol Chem27928462855200410.1074/jbc.M30976920014583606

19 

Y HuangNH ShinY SunKK WangMolecular cloning and characterization of a novel caspase-3 variant that attenuates apoptosis induced by proteasome inhibitionBiochem Biophys Res Commun283762769200110.1006/bbrc.2001.487111350049

20 

C SchwerkK Schulze-OsthoffRegulation of apoptosis by alternative pre-mRNA splicingMol Cell19113200510.1016/j.molcel.2005.05.02615989960

21 

Y TianIH LiaoX ZhanExon expression and alternatively spliced genes in Tourette SyndromeAm J Med Genet B Neuropsychiatr Genet156B7278201110.1002/ajmg.b.3114021184586

22 

R BrinkHF LodishTumor necrosis factor receptor (TNFR)-associated factor 2A (TRAF2A), a TRAF2 splice variant with an extended RING finger domain that inhibits TNFR2-mediated NF-kappaB activationJ Biol Chem27341294134199810.1074/jbc.273.7.41299461607

23 

MA BenedictY HuN InoharaG NunezExpression and functional analysis of Apaf-1 isoforms. Extra Wd-40 repeat is required for cytochrome c binding and regulated activation of procaspase-9J Biol Chem27584618468200010.1074/jbc.275.12.846110722681

24 

LH BoiseM Gonzalez-GarciaCE Postemabcl-x, a bcl-2-related gene that functions as a dominant regulator of apoptotic cell deathCell74597608199310.1016/0092-8674(93)90508-N8358789

25 

YF SunLY YuM SaarmaT TimmuskU ArumaeNeuron-specific Bcl-2 homology 3 domain-only splice variant of Bak is anti-apoptotic in neurons, but pro-apoptotic in non-neuronal cellsJ Biol Chem2761624016247200110.1074/jbc.M01041920011278671

26 

ZN OltvaiCL MillimanSJ KorsmeyerBcl-2 heterodimerizes in vivo with a conserved homolog, Bax, that accelerates programmed cell deathCell74609619199310.1016/0092-8674(93)90509-O8358790

27 

PW NgAG PorterRU JanickeMolecular cloning and characterization of two novel pro-apoptotic isoforms of caspase-10J Biol Chem2741030110308199910.1074/jbc.274.15.1030110187817

Related Articles

Journal Cover

April 2012
Volume 29 Issue 4

Print ISSN: 1107-3756
Online ISSN:1791-244X

Sign up for eToc alerts

Recommend to Library

Copy and paste a formatted citation
x
Spandidos Publications style
Liu J, Xiao Y, Xiong H, Li J, Huang B, Zhang H, Feng D, Chen X and Wang X: Alternative splicing of apoptosis-related genes in imatinib-treated K562 cells identified by exon array analysis. Int J Mol Med 29: 690-698, 2012
APA
Liu, J., Xiao, Y., Xiong, H., Li, J., Huang, B., Zhang, H. ... Wang, X. (2012). Alternative splicing of apoptosis-related genes in imatinib-treated K562 cells identified by exon array analysis. International Journal of Molecular Medicine, 29, 690-698. https://doi.org/10.3892/ijmm.2011.872
MLA
Liu, J., Xiao, Y., Xiong, H., Li, J., Huang, B., Zhang, H., Feng, D., Chen, X., Wang, X."Alternative splicing of apoptosis-related genes in imatinib-treated K562 cells identified by exon array analysis". International Journal of Molecular Medicine 29.4 (2012): 690-698.
Chicago
Liu, J., Xiao, Y., Xiong, H., Li, J., Huang, B., Zhang, H., Feng, D., Chen, X., Wang, X."Alternative splicing of apoptosis-related genes in imatinib-treated K562 cells identified by exon array analysis". International Journal of Molecular Medicine 29, no. 4 (2012): 690-698. https://doi.org/10.3892/ijmm.2011.872