Email updates

Keep up to date with the latest news and content from BMC Cancer and BioMed Central.

Open Access Highly Accessed Research article

Functional features of gene expression profiles differentiating gastrointestinal stromal tumours according to KIT mutations and expression

Jerzy Ostrowski12*, Marcin Polkowski12, Agnieszka Paziewska1, Magdalena Skrzypczak1, Krzysztof Goryca1, Tymon Rubel27, Katarzyna Kokoszyñska78, Piotr Rutkowski3, Zbigniew I Nowecki3, Anna Jerzak Vel Dobosz4, Dorota Jarosz125, Wlodzimierz Ruka3 and Lucjan S Wyrwicz67

Author Affiliations

1 Department of Gastroenterology and Hepatology, Medical Center for Postgraduate Education, Warsaw, Poland

2 Department of Gastroenterology, M Sklodowska-Curie Memorial Cancer Center and Institute of Oncology, Warsaw, Poland

3 Department of Soft Tissue/Bone Sarcoma and Melanoma, M Sklodowska-Curie Memorial Cancer Center and Institute of Oncology, Warsaw, Poland

4 Department of Molecular Biology, M Sklodowska-Curie Memorial Cancer Center and Institute of Oncology, Warsaw, Poland

5 Department of Pathology, M Sklodowska-Curie Memorial Cancer Center and Institute of Oncology, Warsaw, Poland

6 Department of Colorectal Cancer, M Sklodowska-Curie Memorial Cancer Center and Institute of Oncology, Warsaw, Poland

7 Laboratory of Bioinformatics and Systems Biology M Sklodowska-Curie Memorial Cancer Center and Institute of Oncology, Warsaw, Poland

8 Bioinfobank Institute, Poznan, Poland

For all author emails, please log on.

BMC Cancer 2009, 9:413  doi:10.1186/1471-2407-9-413

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2407/9/413


Received:12 June 2009
Accepted:27 November 2009
Published:27 November 2009

© 2009 Ostrowski et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Gastrointestinal stromal tumours (GISTs) represent a heterogeneous group of tumours of mesenchymal origin characterized by gain-of-function mutations in KIT or PDGFRA of the type III receptor tyrosine kinase family. Although mutations in either receptor are thought to drive an early oncogenic event through similar pathways, two previous studies reported the mutation-specific gene expression profiles. However, their further conclusions were rather discordant. To clarify the molecular characteristics of differentially expressed genes according to GIST receptor mutations, we combined microarray-based analysis with detailed functional annotations.

Methods

Total RNA was isolated from 29 frozen gastric GISTs and processed for hybridization on GENECHIP® HG-U133 Plus 2.0 microarrays (Affymetrix). KIT and PDGFRA were analyzed by sequencing, while related mRNA levels were analyzed by quantitative RT-PCR.

Results

Fifteen and eleven tumours possessed mutations in KIT and PDGFRA, respectively; no mutation was found in three tumours. Gene expression analysis identified no discriminative profiles associated with clinical or pathological parameters, even though expression of hundreds of genes differentiated tumour receptor mutation and expression status. Functional features of genes differentially expressed between the two groups of GISTs suggested alterations in angiogenesis and G-protein-related and calcium signalling.

Conclusion

Our study has identified novel molecular elements likely to be involved in receptor-dependent GIST development and allowed confirmation of previously published results. These elements may be potential therapeutic targets and novel markers of KIT mutation status.

Background

Gastrointestinal stromal tumours (GISTs) arise from precursor cells shared with the interstitial cells of Cajal (ICC) [1,2] and encompass a group of heterogeneous neoplasms with different morphology, biologic behaviour, and genetic characteristics [3]. Histopathologically, GISTs are spindle-, epithelioid-, or mixed-cell tumours that usually develop in the wall of the gastrointestinal tract. GISTs can be classified as benign, borderline, or malignant tumours based on tumour size, mitotic index, and the invasion of surrounding tissues, and the majority of these tumours are clinically rather non-aggressive [3].

An early oncogenic event in the majority of GISTs is represented by gain-of-function mutations in either KIT or platelet-derived growth factor receptor α (PDGFRA). Both KIT and PDGFRA belong to the subclass III family of receptor tyrosine kinases [4-6]. The receptor-activating mutations lead to self-phosphorylation of a kinase domain, with the subsequent activation of the JAK/STAT, PI3K/AKT, Ras/ERK, and PLC-γ intracellular pathways in a ligand-independent manner, transmitting mitogenic signals [7-17].

Although mutations in KIT and PDGFRA contribute to tumour development through similar pathways, they correlate with certain clinicopathological features and different responses to imatinib treatment [3]. Moreover, GISTs with different mutation types exhibit differential gene expression at the mRNA [18,19] and protein [20] levels.

Two previous studies [18,19] reported differences between the gene expression profile and pattern of oncogenic mutations. Both studies and additional analyses have confirmed the unexpected observation that a mutation of KIT or PDGFRA is associated with its increased expression at the mRNA level, but in terms of further conclusions Subramanian et al. [19] and Kang et al. [18] are rather discordant. Subramanian and colleagues selected 1875 of almost 28 000 genes or ESTs (expressed sequence tags) clusters represented on cDNA microarray that passed filtering criteria and used it for further analysis. Of these selected genes, 338 were differentially expressed between GISTs assigned to a KIT exon 11 mutation and other types of mutations. A total of 270 genes were differentially expressed between GISTs with a PDGFRA mutation and other GISTs. Notably, a PDGFRA mutation was observed in only 8 of 26 analyzed samples.

In contrast, Kang et al. [18], using high-density spotted oligonucleotide microarrays, selected 4693 out of 18 664 oligonucleotides representing LEADS™ clusters. Among this set of pre-selected genes, only 70 were differentially expressed between GISTs exhibiting different mutation status. Of these 70, Subramanian et al. found only 13 (19%) to be differentially expressed. Both groups reported that on the basis of gene expression signatures, GISTs harbouring different types of mutations could not to be perfectly distinguished. Moreover, because of the far-from-complete coverage of the human genome using the methods in these studies, only limited functional annotations were reported. Thus, although these two important studies have been published, major questions about GIST biology remain open.

To clarify the molecular characteristics of differentially expressed genes according to receptor status, we combined microarray-based data with functional annotations. We selected a model of gastric GIST to obtain a balanced set of tumours with mutations in either KIT or PDGFRA [21]. Significant differences in the molecular makeup of the two groups of gastric GISTs allowed the development of novel functional hypotheses regarding the transduction of intracellular signalling contributing to GIST development.

Methods

Patients

Between April 2005 and March 2008, 31 patients with a diagnosis of gastric GIST were prospectively selected for the study. All patients underwent tumour surgical resection through laparotomy in the Department of Soft Tissue/Bone Sarcoma and Melanoma, and the final diagnosis was obtained from the analysis of clinicopathological findings (Table 1). The study protocol was approved by the Cancer Center Bioethical Committee, and all patients signed informed consent before inclusion. The morphological diagnosis was confirmed by standard H&E staining and immunoreactivity to KIT (CD117) [22,23] on routinely formalin-fixed paraffin-embedded specimens. One to two tumour fragments, depending on tumour size, were snap frozen and stored at -72°C until use. Then, collections of cryostat sections were prepared from different parts of each tumour fragment. Upper and lower sections from each cryosection collection were evaluated by the pathologist (DJ) to control the relative content of non-tumour cells, and the remaining internal portion of the specimen was used in the study if it contained ≥95% tumour cells. Genomic DNA from tissue samples was purified using the DNeasy Tissue Kit, and total RNA was isolated using the RNeasy Mini Kit (both kits from Qiagen GmbH, Hilden, Germany).

Table 1. Patient clinical, pathological and molecular characteristics of analyzed GIST samples

KIT/PDGFRA genotyping and real-time RT-PCR analysis

DNA samples were tested for hot-spot mutation sites of KIT (exons 9, 11, 13, 14, and 17) and PDGFRA (exons 12, 14, and 18) by PCR amplification using primers and annealing temperatures as previously described [24,25]. PCR products were sequenced in two directions by fluorescent dideoxysequencing on an ABI Prism 3100 Sequence Detection System (Applied Biosystems, Foster City, CA).

Specific RNA concentrations were determined by real-time reverse transcriptase (RT)-PCR. Total tissue RNA was isolated with the RNeasy Mini Kit and QIAshredder columns (Qiagen GmbH, Hilden, Germany). Reverse transcription was performed with the SuperScript II Reverse Transcriptase reagent set (Invitrogen Co., Carlsbad, CA) according to the manufacturer's instructions. Quantitative evaluation of mRNA was performed on an ABI Prism 7000 Sequence Detection System with a 25-μl reaction mixture containing 12.5 μl 2× SYBR Green PCR Master Mix (Applied Biosystems, Foster City, CA), 5 μl cDNA, and 50 nM primers. Oligonucleotide primers for the analyzed KIT/PDGFRA transcripts were designed using Primer Express Software (Applied Biosystems, Foster City, CA) and are listed in Supplementary Table S1 (see Additional file 1). For each run, standard curves were generated for a primer set by serial dilution of pooled cDNA to counterbalance variations in PCR reaction efficiency. Melting curves were generated after each reaction to verify the melting temperature of the amplicon. In addition, the purity of the RT-PCR product was verified by agarose gel electrophoresis. To normalize nonspecific variations in real-time PCR, the normalization factor was calculated as the geometric mean of RNA concentrations of three control genes, glyceraldehyde-3-phosphate dehydrogenase, ubiquitin C, and β-actin.

Additional file 1. Supplementary Tables S1-S6. Supplementary Table S1. Primers used in this study. Supplementary Table S2. Probesets discriminating samples according to KIT mutation status (adjusted p.val < 0,1; FC>2). Supplementary Table S3. Probesets discriminating samples according to KIT transcript levels (adjusted p.val < 0,1; FC>2). Supplementary Table S4. A. Gene ontology terms overrepresented in probesets differentiating samples according to mutation status. B. Gene ontology terms overrepresented in probesets differentiating samples according to expression of KIT. Supplementary Table S5. Probesets with expression changed at least 2 fold between samples with low and high KIT expression annotated to selected GO terms: A. G-protein interaction (GO:0007186), B. synaptic transmission (GO:0007268), C. blood vessel development (GO:0001568). Supplementary Table S6. Interactomes of KIT and PDGFRA receptors assembled from literature and protein interaction databases (data sources: Bg - BioGrid, B - BOND database, H - HPRD, P - PubMed). Common proteins are highlighted in yellow.

Format: DOC Size: 1.8MB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Gene expression analyses on microarrays

Gene expression profiling was carried out using Affymetrix oligonucleotide microarrays (GeneChip HG-U133plus2) as described previously [26].

To obtain gene expression measurements, the extraction of probe-level data was performed with a standard GC-RMA algorithm for background correction and summarization steps and least-variant set algorithm for normalization based on a least-variant set of probe sets. The calculations were performed using BioConductor (version 2.8.1) packages gcrma (version 2.14.1) and FLUSH.LVS.bundle (version 1.2.1, proportion = 0.6). To test the internal consistency of the data sets, we used principal component analysis, normalized unscaled standard error plots, and relative log expression plots.

The measured expression levels were log transformed (log2). For data filtration, we selected the probe sets exhibiting signal intensity above the threshold limit, which was established at the 95th percentile of the expression levels from Y-chromosome-linked probe set signals detectable in female samples. The low-expression (marginal) probe sets with levels below the threshold in at least 19 samples were rejected.

To establish gene expression profiles, differentially expressed probe sets in the pair-wise comparisons were identified using the Kruskal-Wallis test. The resulting P values were adjusted for testing of multiple hypotheses using the Benjamini-Hochberg procedure that controls a false discovery rate. The false discovery rate threshold was set to 0.1, and only probe sets exhibiting a minimum two-fold change in mean relative expression were included in the gene lists. Cluster analysis of probe sets exhibiting differential expression was also performed. The probe sets were divided into groups of distinct expression patterns by an evolutionary-driven k-means clustering algorithm with a distance metric derived from the Pearson correlation coefficient. Unsupervised average-linkage hierarchical clustering and PCA were used for a graphic summary and evaluation of relationships between samples.

Both statistical and clustering analyses were performed using a proprietary software working in the MATLAB (MathWorks) and Bioconductor (2.8.1) environments.

Functional analyses of gene expression by Gene Ontology

Differentially expressed probe sets were annotated with Gene Ontology (GO) terms (GO.db version 2.2.5) using the Bioconductor packages GOstats (version 2.8.0, Affymetrix HG-U133 Plus 2.0 Array Annotation Data) and package annotate (version 1.20.1). The significance of differential representation of GO terms between specified lists of probe sets was determined by the hypergeometric test implemented in GOstats (version 2.8.0). P values returned by GOstats were corrected for testing of multiple hypotheses with the Benjamini-Hochberg method implemented in an R environment (version 2.8.1, The R Foundation for Statistical Computing; http://www.r-project.org webcite). Adjusted P values of less than 0.1 were considered significant.

Models of KIT and PDGFRA (a and b subunits) signalling pathways were prepared on the basis of three databases: Biogrid (The Biological General Repository for Interaction Datasets) [27], HPRD (Human Protein Reference Database) [28], and BIND (Biomolecular Interaction Network Database) [29] and additional literature searches.

Results

KIT and PDGFRA mutation profiling

Gene expression profiles from 29 out of 31 primary gastric GISTs were selected for the molecular analysis. Two samples were rejected due to poor quality of extracted information after MAS5.0 testing of criteria suggested by the producer. Of the 29 cases, 24, 2, and 3 cases were classified as benign (very low or low risk), borderline (intermediate risk), and malignant (high risk), respectively (Table 1). Such a distribution is consistent with the general clinical picture of gastric GISTs [21]. Genotyping revealed that 15 tumours had KIT mutations (14 in exon 11 and one in exon 9), 11 tumours had PDGFRA mutations (9 in exon 18 and two in exon 12), and 3 were KIT/PDGFRA wild-type GISTs within the analyzed mutation hotspots (Table 1). A total of 23 tumours were CD117/KIT-positive on immunohistochemistry (IHC). Of these, 13 and 8 tumours had KIT or PDGFRA mutations, respectively. In addition, two wild-type tumours were also KIT-positive. Thus, expression of KIT in IHC did not correlate with receptor tyrosine kinase mutation status (Table 1).

Microarray results and mutation status

From 54 675 probe sets of the Affymetrix HGU133plus2 microarray, 16 880 probe sets passed the filtering procedure. Among the probe sets with detectable expression levels, the non-parametric Kruskal-Wallis test revealed 970 (311 with FC>2; FC = fold change) probe sets differentially expressed between tumours with KIT and PDGFRA mutations. Of these 311 (FC>2), 109 probe sets (81 genes) were upregulated, and 202 probe sets (143 genes) were downregulated in samples from tumours harbouring KIT mutations compared to those with PDGFRA mutations. Supplementary Table S2 (see Additional file 1) gives the complete list of differentially expressed probe sets. As expected, unsupervised hierarchical clustering (Figure 1) of differentially expressed genes distinguished GIST samples according to the mutation status.

thumbnailFigure 1. Unsupervised hierarchical clustering for the selection of differentially expressed genes in GIST tumours according to KIT mutation. Across the top, individual tumour samples are arrayed in a column (upper: blue - KIT mutation; red - PDGFRA mutation; gray- no mutation found; lower: green - low KIT expression/high PDGFRA expression; yellow- high KIT expression/low PDGFRA expression; dark gray- high KIT/PDGFRA expression); on the left side, 311 individual probe sets differentiating tumours in accordance with the mutation are shown in rows. The colour in each cell reflects the level of expression of the corresponding probe set in the corresponding array sample relative to its mean level of expression estimated for the entire set of samples. Red indicates expression levels greater than the mean, and green indicates lower than the mean.

As reported previously [18,19], increased expression of KIT and PDGFRA follows their mutation status. To evaluate the impact of mutation status on expression of genes encoding both receptors, quantitative RT-PCR analysis was performed simultaneously on the same RNA samples used in the microarray analysis. An overall good correlation was observed between quantitative changes in KIT expression levels obtained by microarrays and quantitative RT-PCR (not shown). However, while the microarray signal intensity of PDGFRA probe sets was in most arrays slightly above the threshold limit, RT-PCR allowed for more reliable quantification of PDGFRA transcript levels.

Overexpression of KIT and PDGFRA was closely related to receptor mutation status. As shown in Figure 2, in 14 out of 15 GISTs with a KIT mutation and in 2 out of 10 GISTs with a PDGFRA mutation, the relative expression of KIT was 1.47 arbitrary units (a.u.) or greater, while in the remaining tumours, it was 0.37 a.u. or lower. In contrast, in nine tumours with a PDGFRA mutation and one tumour with a KIT mutation, the expression of PDGFRA was ≤1.32 a.u., while in the remaining GISTs with KIT/PDGFRA mutations, it was ≤0.84 a.u. One tumour with a PDGFRA mutation exhibited overexpression of both KIT and PDGFRA. Among receptor wild-type GISTs, the expression profile for one tumour was typical for tumours with PDGFRA mutations and was typical in two others for tumours with a KIT mutation. These findings suggest the presence of a potential additional region undergoing oncogenic mutations in both analyzed genes.

thumbnailFigure 2. Relative mRNA expression of KIT (left panel) and PDGFRA (right panel) in GIST tumours according to sample mutation status. Black circles - KIT mutations; gray circles - PDGFRA mutations

When selection of differentially expressed genes was performed according to the discriminating KIT transcript levels, the microarray data from the two groups of 15 and 13 tumour samples with either high or low KIT expression were included in analyses. At this point, 3116 (680 with the FC>2) probe sets indicated genes that were differentially expressed. (see Additional file 1: Supplementary Table S3). Of these, 253 probe sets (171 genes) were upregulated, and 427 probe sets (299 genes) were downregulated in tumours with overexpressed KIT compared to those with low expression. Again, unsupervised hierarchical clustering showed that gene expression patterns of the individual tumour samples from these two groups clustered together (see Additional file 2: Supplementary Figure S1).

Additional file 2. Supplementary Figure S1. Unsupervised hierarchical clustering for the selection of differentially expressed genes in GIST tumours according to KIT mutation. Across the top, individual tumour samples are arrayed in a column (upper: green - low KIT expression/high PDGFRA expression; yellow- high KIT expression/low PDGFRA expression; dark gray- high KIT/PDGFRA expression; lower: blue - KIT mutation; red - PDGFRA mutation; gray- no mutation found); on the left side, 680 individual probe sets differentiating tumours in accordance with the KIT expression are shown in rows. The colour in each cell reflects the level of expression of the corresponding probe set in the corresponding array sample relative to its mean level of expression estimated for the entire set of samples. Red indicates expression levels greater than the mean, and green indicates lower than the mean.

Format: PDF Size: 104KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

No gene expression signature was found to be associated with tumour size, mitotic index, or risk category.

Gene enrichment analyzes

Differential expression was then analyzed on the level of GO categories and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways. Functional features of the 311 and 680 genes that differentiated tumours according to receptor mutations and expression allowed selection of respectively nine and seven subcategories, as analyzed by their annotation to GO terms (see Additional file 1: Supplementary Table S4). Overrepresentation (estimated by the highest level of significance) was represented by blood vessel development (GO:0001568), angiogenesis (GO:0001525), cell adhesion (GO:0007155), and G-protein coupled receptor protein signalling pathway (GO:0007186). Terms that exhibited lower levels of statistical significance were represented by elevation of cytosolic calcium ion concentration (GO:0007204), transmembrane receptor protein tyrosine kinase signalling pathway (GO:0007169), and cellular calcium ion homeostasis (GO:0006874). In similar pair-wise comparisons performed for KEGG pathways, differentially regulated pathways were represented by calcium signalling pathway and cytokine-cytokine receptor interaction.

Interestingly, the expression of 7 genes annotated to synaptic transmission, 15 genes annotated to blood vessel development, and 20 genes annotated to G-protein signalling were at least 2-fold higher in tumours with low KIT expression compared to those with high KIT expression (see Additional file 1: Supplementary Table S5).

Because previous studies indicated that KIT-mutant and PDGFRA-mutant GISTs may have features associated with activation of downstream pathways like ERK1/2, AKT, p70/85S6K, STAT1/STAT3, and PI3K/mTOR [5,7-17], the lists of differentially expressed genes were compared with both lists of interacting partners of KIT and PDGFR, summarized in Supplementary Table S5 (see Additional file 1). According to our model of PDGFRA signalling pathways, the final data set contained 44 and 52 proteins interacting with KIT and PDGFR, respectively, 13 of which were common for both pathways (see Additional file 1: Supplementary Table S6). Within the list of differentially expressed genes according to the mutation status only 3 genes corresponding to PDGFRA interactome were found. None of KIT genes encoding KIT-signalling pathway proteins was found within this list.

To further test other pathways described by Corless et al. [5], we also analyzed genes from selected KEGG pathways involved in mitogenic signal transduction (mTOR, Jak-STAT, MAPK, TGF-beta, calcium signalling, vascular endothelial growth factor (VEGF) signalling). Those genes were compared with lists of differentially expressed genes depending on KIT/PDGFRA mutation status. Apart from protein kinase C (PKC)-alpha, no other genes from such lists were found among the differentially expressed genes (data not shown).

Among PKC isoforms analyzed in this study by microarray (alpha, beta isoform 1, variant 1, eta, iota, zeta) and quantitative RT-PCR (alpha, beta isoform 1, variant 1 and 2, delta, epsilon, theta) (data not shown), expression of PKC-alpha was significantly lower (FC = 0.17; P < 0.00016), while expression of PKC-theta was significantly higher (FC = 2.34; P < 0.00016) in tumours with KIT mutations compared to those with PDGFRA mutations and wild-type tumours.

Discussion

GISTs express KIT, a 145-kD transmembrane glycoprotein that serves as the receptor for stem-cell factor [30,31]. KIT activates cellular signalling during embryogenesis [32,33] and is critical for the development of germ cells, hematopoietic progenitor cells, and mast cells [34]. Its function is closely related to that of other receptor tyrosine kinases for PDGF, macrophage colony stimulating factor, and FLT3 ligand [35]. Activating mutations in the KIT and PDGFRA genes result in ligand-independent activation of their receptor tyrosine kinase function, which may transmit early oncogenic signals in the majority of GISTs (Rutkowski et al., 2008).

Even though the biological consequences of KIT and PDGFRA mutations seem to be similar, our study has revealed hundreds of differentially expressed genes that group the tumours according to the receptor mutation status and receptor gene expression. However, although many of these genes may be involved in receptor-specific alterations of GIST intracellular signalling pathways, we identified no discriminative profiles of gene expression associated with clinical or pathological outcomes. Most of the discriminative genes were found to be upregulated in PDGFRA-mutated GISTS.

To further clarify if gene signatures that group GISTs according to KIT/PDGFRA mutation status, as described previously [18,19], may be defined also at the level of intracellular signalling pathways, we analyzed microarray data in the context of functional annotation.

Among terms and pathways in the current analysis with the highest overrepresentation in GISTs with mutated and/or overexpressed PDGFRA were "blood vessel development" and "angiogenesis". In fact, GISTs are highly vascularised tumours, and VEGF expression has been postulated to be a KIT-genotype-independent adverse prognostic indicator for early treatment failure and poor survival of GIST patients on imatinib therapy [36,37].

We also found that the functional features of genes differentially expressed between the two groups of GISTs were represented by the G-protein-coupled receptor protein signalling pathway (GO:0007186). The regulated secretion of transmitters and hormones, a characteristic event of neuroendocrine cells and tumours, is controlled by G-protein-coupled membrane receptors. Indirect evidence of neural or neuroendocrine phenotypes including high expression of this type of receptor have been described recently in GISTs [5].

As part of this study, we compiled lists of differentially expressed genes for comparison with lists of interacting partners of KIT and PDGFR. This analysis identified significantly lower expression of PKC-alpha and significantly higher expression of PKC-theta in tumours with KIT mutations compared to those with PDGFRA mutations or wild-type tumours. The PKC family consists of 10 related serine/threonine protein kinases, which are involved in regulation of cell proliferation, survival, and death. In addition, some are considered to be tumour promoters that may enhance multiple cellular oncogenic signalling pathways [38,39]. While the alpha, beta, epsilon, and atypical PKCs possess anti-apoptotic action, the delta and theta isoforms usually promote apoptosis [40]. Interestingly, PKC-theta has been selected previously as a sensitive marker of GISTs [41].

The enrichment analysis based on GO or similar annotations gives reliable results for sets of hundreds of genes, and such sets are likely to be identified with microarrays covering almost the whole transcriptome, as oligonucleotide microarrays which were used in this study. Experiments exploiting older technologies, like those used by Subramanian et al. [19] hold the middle ground between traditional single-gene techniques and high-throughput implementations and cannot provide sufficient data for systemic interpretation.

While the noise inherent in microarray technology often complicates the process of data interpretation, both the array quality and the choice of analytical processing methods have a major impact on differential expression analysis of microarray data [26]. Thus, stringent selection criteria are essential for identifying differentially expressed genes. In the case of a set of thousands of transcripts, P values and FC criteria are not sufficient if not coupled with an adequate P value correction method for simultaneous testing of multiple hypotheses (e.g., Bonferroni, Benjamini-Hochberg). An apparent lack of such a step in the study reported by Kang et al. (2005) is a likely reason for low concordance with other studies. With over 4000 genes tested, a P < 0.01 criterion leads to about 40 false positives (i.e., roughly 57% of a reported 70-gene data set). In the current work, we sought to address these issues through appropriate statistical adjustment for multiple hypotheses.

Conclusion

To summarized, our study has identified novel molecular mechanisms likely to be involved in receptor-dependent GIST development and allowed confirmation of previously published results. These observations may be useful for the development of molecular markers that might predict which GIST patients will experience an adequate response to proposed therapy. However, before health-care professionals can see benefits from molecular diagnostics, a full understanding of the biological processes underlying GIST development is required.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

JO, MP - planned the analysis; AP, MS - performed gene expression studies; KG, TR, KK, LSW, JO - analyzed the results; MP, PR, ZIN, WR - provided clinical material; AJD - performed mutation analysis; DJ - provided histopathological expertise; JO, MP, LSW, KG, KK - wrote the manuscript. All authors read and approved the final manuscript.

Acknowledgements

This work was supported by Ministry of Science and Higher Education grants 2P05B03329, PBZ-MNiI-2/1/2005, and N401 050 32/1181.

References

  1. Miettinen M, Lasota J: Gastrointestinal stromal tumors--definition, clinical, histological, immunohistochemical, and molecular genetic features and differential diagnosis.

    Virchows Arch 2001, 438(1):1-12. PubMed Abstract | Publisher Full Text OpenURL

  2. Kindblom LG, Remotti HE, Aldenborg F, Meis-Kindblom JM: Gastrointestinal pacemaker cell tumor (GIPACT): gastrointestinal stromal tumors show phenotypic characteristics of the interstitial cells of Cajal.

    Am J Pathol 1998, 152(5):1259-1269. PubMed Abstract | PubMed Central Full Text OpenURL

  3. Rutkowski P, Debiec-Rychter M, Ruka W: Gastrointestinal stromal tumors: key to diagnosis and choice of therapy.

    Mol Diagn Ther 2008, 12(3):131-143. PubMed Abstract OpenURL

  4. Hirota S, Isozaki K, Moriyama Y, Hashimoto K, Nishida T, Ishiguro S, Kawano K, Hanada M, Kurata A, Takeda M, et al.: Gain-of-function mutations of c-kit in human gastrointestinal stromal tumors.

    Science 1998, 279(5350):577-580. PubMed Abstract | Publisher Full Text OpenURL

  5. Corless CL, Fletcher JA, Heinrich MC: Biology of gastrointestinal stromal tumors.

    J Clin Oncol 2004, 22(18):3813-3825. PubMed Abstract | Publisher Full Text OpenURL

  6. Hirota S, Ohashi A, Nishida T, Isozaki K, Kinoshita K, Shinomura Y, Kitamura Y: Gain-of-function mutations of platelet-derived growth factor receptor alpha gene in gastrointestinal stromal tumors.

    Gastroenterology 2003, 125(3):660-667. PubMed Abstract | Publisher Full Text OpenURL

  7. Ihle JN, Kerr IM: Jaks and Stats in signaling by the cytokine receptor superfamily.

    Trends Genet 1995, 11(2):69-74. PubMed Abstract | Publisher Full Text OpenURL

  8. Schindler C, Darnell JE Jr: Transcriptional responses to polypeptide ligands: the JAK-STAT pathway.

    Annu Rev Biochem 1995, 64:621-651. PubMed Abstract | Publisher Full Text OpenURL

  9. Weiler SR, Mou S, DeBerry CS, Keller JR, Ruscetti FW, Ferris DK, Longo DL, Linnekin D: JAK2 is associated with the c-kit proto-oncogene product and is phosphorylated in response to stem cell factor.

    Blood 1996, 87(9):3688-3693. PubMed Abstract | Publisher Full Text OpenURL

  10. Brizzi MF, Dentelli P, Rosso A, Yarden Y, Pegoraro L: STAT protein recruitment and activation in c-Kit deletion mutants.

    J Biol Chem 1999, 274(24):16965-16972. PubMed Abstract | Publisher Full Text OpenURL

  11. Mui AL: The role of STATs in proliferation, differentiation, and apoptosis.

    Cell Mol Life Sci 1999, 55(12):1547-1558. PubMed Abstract | Publisher Full Text OpenURL

  12. Porfiri E, McCormick F: Regulation of epidermal growth factor receptor signaling by phosphorylation of the ras exchange factor hSOS1.

    J Biol Chem 1996, 271(10):5871-5877. PubMed Abstract | Publisher Full Text OpenURL

  13. Burack WR, Sturgill TW: The activating dual phosphorylation of MAPK by MEK is nonprocessive.

    Biochemistry 1997, 36(20):5929-5933. PubMed Abstract | Publisher Full Text OpenURL

  14. Zhang C, Schmidt M, von Eichei-Streiber C, Jakobs KH: Inhibition by toxin B of inositol phosphate formation induced by G protein-coupled and tyrosine kinase receptors in N1E-115 neuroblastoma cells: involvement of Rho proteins.

    Mol Pharmacol 1996, 50(4):864-869. PubMed Abstract | Publisher Full Text OpenURL

  15. Barker SA, Caldwell KK, Pfeiffer JR, Wilson BS: Wortmannin-sensitive phosphorylation, translocation, and activation of PLCgamma1, but not PLCgamma2, in antigen-stimulated RBL-2H3 mast cells.

    Mol Biol Cell 1998, 9(2):483-496. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Timokhina I, Kissel H, Stella G, Besmer P: Kit signaling through PI 3-kinase and Src kinase pathways: an essential role for Rac1 and JNK activation in mast cell proliferation.

    Embo J 1998, 17(21):6250-6262. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Schlessinger J: Cell signaling by receptor tyrosine kinases.

    Cell 2000, 103(2):211-225. PubMed Abstract | Publisher Full Text OpenURL

  18. Kang HJ, Nam SW, Kim H, Rhee H, Kim NG, Kim H, Hyung WJ, Noh SH, Kim JH, Yun CO, et al.: Correlation of KIT and platelet-derived growth factor receptor alpha mutations with gene activation and expression profiles in gastrointestinal stromal tumors.

    Oncogene 2005, 24(6):1066-1074. PubMed Abstract | Publisher Full Text OpenURL

  19. Subramanian S, West RB, Corless CL, Ou W, Rubin BP, Chu KM, Leung SY, Yuen ST, Zhu S, Hernandez-Boussard T, et al.: Gastrointestinal stromal tumors (GISTs) with KIT and PDGFRA mutations have distinct gene expression profiles.

    Oncogene 2004, 23(47):7780-7790. PubMed Abstract | Publisher Full Text OpenURL

  20. Kang HJ, Koh KH, Yang E, You KT, Kim HJ, Paik YK, Kim H: Differentially expressed proteins in gastrointestinal stromal tumors with KIT and PDGFRA mutations.

    Proteomics 2006, 6(4):1151-1157. PubMed Abstract | Publisher Full Text OpenURL

  21. Lasota J, Miettinen M: KIT and PDGFRA mutations in gastrointestinal stromal tumors (GISTs).

    Semin Diagn Pathol 2006, 23(2):91-102. PubMed Abstract | Publisher Full Text OpenURL

  22. Tran T, Davila JA, El-Serag HB: The epidemiology of malignant gastrointestinal stromal tumors: an analysis of 1,458 cases from 1992 to 2000.

    Am J Gastroenterol 2005, 100(1):162-168. PubMed Abstract | Publisher Full Text OpenURL

  23. Perez EA, Livingstone AS, Franceschi D, Rocha-Lima C, Lee DJ, Hodgson N, Jorda M, Koniaris LG: Current incidence and outcomes of gastrointestinal mesenchymal tumors including gastrointestinal stromal tumors.

    J Am Coll Surg 2006, 202(4):623-629. PubMed Abstract | Publisher Full Text OpenURL

  24. Lasota J, vel Dobosz AJ, Wasag B, Wozniak A, Kraszewska E, Michej W, Ptaszynski K, Rutkowski P, Sarlomo-Rikala M, Steigen SE, et al.: Presence of homozygous KIT exon 11 mutations is strongly associated with malignant clinical behavior in gastrointestinal stromal tumors.

    Lab Invest 2007, 87(10):1029-1041. PubMed Abstract | Publisher Full Text OpenURL

  25. Rutkowski P, Debiec-Rychter M, Nowecki ZI, Wozniak A, Michej W, Limon J, Siedlecki JA, Jerzak Vel Dobosz A, Grzesiakowska U, Nasierowska-Guttmejer A, et al.: Different factors are responsible for predicting relapses after primary tumors resection and for imatinib treatment outcomes in gastrointestinal stromal tumors.

    Med Sci Monit 2007, 13(11):CR515-522. PubMed Abstract | Publisher Full Text OpenURL

  26. Ostrowski J, Mikula M, Karczmarski J, Rubel T, Wyrwicz LS, Bragoszewski P, Gaj P, Dadlez M, Butruk E, Regula J: Molecular defense mechanisms of Barrett's metaplasia estimated by an integrative genomics.

    J Mol Med 2007, 85(7):733-743. PubMed Abstract | Publisher Full Text OpenURL

  27. Stark C, Breitkreutz BJ, Reguly T, Boucher L, Breitkreutz A, Tyers M: BioGRID: a general repository for interaction datasets.

    Nucleic Acids Res 2006, (34 Database):D535-539. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Peri S, Navarro JD, Kristiansen TZ, Amanchy R, Surendranath V, Muthusamy B, Gandhi TK, Chandrika KN, Deshpande N, Suresh S, et al.: Human protein reference database as a discovery resource for proteomics.

    Nucleic Acids Res 2004, (32 Database):D497-501. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Gilbert D: Biomolecular interaction network database.

    Brief Bioinform 2005, 6(2):194-198. PubMed Abstract | Publisher Full Text OpenURL

  30. Besmer P, Murphy JE, George PC, Qiu FH, Bergold PJ, Lederman L, Snyder HW Jr, Brodeur D, Zuckerman EE, Hardy WD: A new acute transforming feline retrovirus and relationship of its oncogene v-kit with the protein kinase gene family.

    Nature 1986, 320(6061):415-421. PubMed Abstract | Publisher Full Text OpenURL

  31. Yarden Y, Kuang WJ, Yang-Feng T, Coussens L, Munemitsu S, Dull TJ, Chen E, Schlessinger J, Francke U, Ullrich A: Human proto-oncogene c-kit: a new cell surface receptor tyrosine kinase for an unidentified ligand.

    Embo J 1987, 6(11):3341-3351. PubMed Abstract | PubMed Central Full Text OpenURL

  32. Besmer P: The kit ligand encoded at the murine Steel locus: a pleiotropic growth and differentiation factor.

    Curr Opin Cell Biol 1991, 3(6):939-946. PubMed Abstract | Publisher Full Text OpenURL

  33. Sanders KM, Ordog T, Koh SD, Torihashi S, Ward SM: Development and plasticity of interstitial cells of Cajal.

    Neurogastroenterol Motil 1999, 11(5):311-338. PubMed Abstract | Publisher Full Text OpenURL

  34. Nocka K, Majumder S, Chabot B, Ray P, Cervone M, Bernstein A, Besmer P: Expression of c-kit gene products in known cellular targets of W mutations in normal and W mutant mice--evidence for an impaired c-kit kinase in mutant mice.

    Genes Dev 1989, 3(6):816-826. PubMed Abstract | Publisher Full Text OpenURL

  35. Rousset D, Agnes F, Lachaume P, Andre C, Galibert F: Molecular evolution of the genes encoding receptor tyrosine kinase with immunoglobulinlike domains.

    J Mol Evol 1995, 41(4):421-429. PubMed Abstract | Publisher Full Text OpenURL

  36. Salto-Tellez M, Nga ME, Han HC, Wong AS, Lee CK, Anuar D, Ng SS, Ho M, Wee A, Chan YH, et al.: Tissue microarrays characterise the clinical significance of a VEGF-A protein expression signature in gastrointestinal stromal tumours.

    Br J Cancer 2007, 96(5):776-782. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. McAuliffe JC, Lazar AJ, Yang D, Steinert DM, Qiao W, Thall PF, Raymond AK, Benjamin RS, Trent JC: Association of intratumoral vascular endothelial growth factor expression and clinical outcome for patients with gastrointestinal stromal tumors treated with imatinib mesylate.

    Clin Cancer Res 2007, 13(22 Pt 1):6727-6734. PubMed Abstract | Publisher Full Text OpenURL

  38. Reyland ME: Protein kinase C isoforms: Multi-functional regulators of cell life and death.

    Front Biosci 2009, 14:2386-2399. PubMed Abstract | Publisher Full Text OpenURL

  39. Martiny-Baron G, Fabbro D: Classical PKC isoforms in cancer.

    Pharmacol Res 2007, 55(6):477-486. PubMed Abstract | Publisher Full Text OpenURL

  40. Gutcher I, Webb PR, Anderson NG: The isoform-specific regulation of apoptosis by protein kinase C.

    Cell Mol Life Sci 2003, 60(6):1061-1070. PubMed Abstract | Publisher Full Text OpenURL

  41. Lee HE, Kim MA, Lee HS, Lee BL, Kim WH: Characteristics of KIT-negative gastrointestinal stromal tumours and diagnostic utility of protein kinase C theta immunostaining.

    J Clin Pathol 2008, 61(6):722-729. PubMed Abstract | Publisher Full Text OpenURL

Pre-publication history

The pre-publication history for this paper can be accessed here:

http://www.biomedcentral.com/1471-2407/9/413/prepub