Skip to main content

Advertisement

Gene expression identifies heterogeneity of metastatic behavior among gastrointestinal stromal tumors

Journal of Translational Medicine volume 14, Article number: 51 (2016) Cite this article

Abstract

Background

Adjuvant imatinib is useful in patients with gastrointestinal stromal tumors (GIST) at high risk of recurrence. At present, the risk of recurrence is determined based on tumor size, mitotic rate, tumor site, and tumor rupture. Previous studies using various biochemical pathways identified gene expression patterns that distinguish two subsets of aggressive fibromatosis (AF), serous ovarian carcinoma (OVCA), and clear cell renal cell carcinoma (RCC). These gene sets separated soft tissue sarcomas into two groups with different probabilities of developing metastatic disease. The present study used these gene sets to identify GIST subgroups with different probabilities of developing metastatic disease.

Methods

We utilized these three gene sets, hierarchical clustering, and Kaplan–Meier analysis, to examine 60 primary resected GIST samples using Agilent chip expression profiling.

Results

Hierarchical clustering using both the combined and individual AF-, OVCA-, and RCC- gene sets identified differences in probabilities of developing metastatic disease between the clusters defined by the first branch point of the clustering dendrograms (p = 0.029 for the combined gene set, p = 0.003 for the AF-gene set, p < 0.001 for the OVCA-gene set, and p = 0.003 for the RCC-gene set).

Conclusions

Hierarchical clustering using these gene sets identified at least two subsets of GIST with distinct clinical behavior and risk of metastatic disease. The use of gene expression analysis along with other known prognostic factors may better predict the long-term outcome following surgery, and thus restrict the use of adjuvant therapy to high-risk GIST, and reduce heterogeneity among groups in clinical trials of new drugs.

Background

Gastrointestinal stromal tumors (GISTs) are the most common sarcoma of the gastrointestinal tract, occurring mostly in the muscular wall of the stomach or small bowel, where it is felt to arise from the interstitial cells of Cajal or similar cells [1, 2]. The primary treatment for GIST is surgical excision, but a significant number of cases recur [3, 4]. Adjuvant imatinib, a tyrosine kinase inhibitor, is useful in select cases of GIST based on risk of recurrence [5–8]. At present, the risk of recurrence is determined based on tumor size, mitotic rate, tumor site, and tumor rupture [1, 5, 8–13], as for example in the Miettinen risk score [11], but more accurate predictors would be useful to better direct therapy.

While most GISTs have mutations in the KIT gene, mutations in the platelet derived growth factor receptor alpha (PDGFRA) gene are also common [1, 2, 5, 8, 14]. In a small percentage of GISTs, mutations in other genes such as BRAF, succinate dehydrogenase (SDH), or neurofibromatosis (NF) may occur [1, 5, 15–20]. The type of KIT or PDGFRA mutation may affect the recurrence rate as well as response to imatinib [5, 8]. Despite the key role of activating mutations of KIT or PDGFRA, GIST biology is also dependent upon other genetic changes [1]. Cases of KIT-mutant GIST have been reported that present with coexisting downstream mutations [5, 8, 21, 22].

Gene expression patterns have been used to predict the development of metastases in soft tissue sarcoma [23–26]. Differences in the gene expression profiles of GISTs with different KIT- or PDGFRA-mutant tumors have been reported [27, 28], and several recent studies have explored the use of gene expression patterns to predict recurrence rate of GIST [29–35].

In previously published studies using various biochemical pathways, we derived gene expression profiles that identified two subgroups of aggressive fibromatosis (AF-gene set), ovarian carcinomas (OVCA-gene set), and clear cell renal cell carcinomas (RCC-gene set) [36–39]. We previously used a gene set derived from these three studies to separate 73 high grade soft tissue sarcoma into 2 or 4 groups with different propensities of metastasis [25]. In an independent study, these gene sets were used to separate 309 high-grade soft tissue sarcoma into 2 or 4 groups with different propensity of metastasis [26].

In the present study, we utilized our three gene sets to examine a group of 60 GISTs using Agilent chip based expression profiling [33]. These gene sets successfully separated the GIST samples into subsets with different probabilities of developing disease recurrence, and may be useful to better predict who would benefit from adjuvant imatinib.

Methods

Samples

Sixty primary tumor samples were obtained from patients who had surgical resection of a GIST, and patients were followed without treatment until tumors recurred as previously described [33]. Frozen samples from resected primary GISTs untreated until tumor recurrence were selected from the European GIST database CONTICAGIST (http://www.conticagist.org). According to French law at the time of the study, experiments were performed in agreement with the Bioethics Law 2004 800 and the Ethics Charter from the National Institute of Cancer; all subjects signed a non-opposition statement for research use of their sample. Total RNA was extracted from each frozen tumor sample, and analyzed on Agilent Whole human 44K Genome Oligo Array (Agilent Technologies) as previously described [33]. Patient characteristics were previously described [33]. These data were kindly provided by Dr. F. Chibon, Institute Bergonie, Bordeaux, France.

Gene sets

Three different previously described gene sets with limited overlap were used: the AF-gene set, OVCA-gene set, and RCC-gene set. These gene sets consist of 161, 173, and 138 known genes respectively [36–39]. The AF-gene set and RCC-gene set distinguished between two subgroups of AF samples and RCC samples, respectively. The OVCA-gene set distinguished borderline from invasive serous OVCA. These three gene sets were pooled resulting in a combined gene set.

Hierarchical clustering and fold-change analysis

The AF-, OVCA-, and RCC-gene sets were used individually or combined, to cluster the 60 primary GIST samples. For clustering, genes were median centered, normalized, and then clustered by complete hierarchical clustering using uncentered correlation with Eisen clustering software [40] and viewed using the TreeView software (http://www.rana.lbl.gov) [41].

Analysis of time to metastasis

For each data set, we used the Kaplan–Meier (K-M) method to calculate metastasis-free survival probabilities, and cumulative probabilities of metastasis (one minus survival probabilities) at critical time points (2, 4, 6, and 8 years). p values were calculated by using the log-rank test for comparing different groups. p values ≤0.05 were considered statistically significant. Analyses were performed in R version 3.0.1 [42].

Results

Analysis of GIST samples using the individual AF-, OVCA-, and RCC-gene sets

We analyzed 60 GIST samples with the individual AF-, RCC-, and OVCA-probe sets; patient characteristics have been previously reported [23]. Hierarchical clustering of the GIST samples using each individual gene set (Additional file 1: Figure S1) identified differences in time to metastasis when the GIST samples were analyzed as two groups defined by the first branch point of the clustering. For the AF-gene set, the probability of not developing metastases by 6 years in Group B was 0.54 while none of the patients in Group A recurred (Fig. 1a; Table 1A, p = 0.003). For the OVCA-gene set, the probability of not developing metastases at 6 years was 0.20 in Group B vs 0.97 in group A (Fig. 1b; Table 1B, p < 0.001). For the RCC-gene set, the probability of not developing metastases at 6 years was 0.46 for Group B and 0.90 for Group A (Fig. 1c; Table 1C, p = 0.003).

Fig. 1
figure1

Kaplan–Meier analysis of the time to development of metastases of two groups (groups A and B) defined by the first break point of the hierarchical clustering of 60 GIST samples. GIST samples were analyzed using the individual AF-gene set (a), OVCA-gene set (b), RCC-gene set (c), and the combined gene set (d). The time to development of metastasis differed between groups A and B (a, p = 0.003; b, p < 0.001; c, p = 0.003; and d, p = 0.029). GIST samples were also analyzed by separating the samples into 3 well-separated clusters using the combined gene set (e, p < 0.001), and by identifying samples as “good” or “poor” prognosis by any 2 of the 3 gene sets as described in the text (f, p < 0.001)

Full size image
Table 1 Probability (95 % CI) of no metastasis as a function of time
Full size table

Analysis of GIST samples using the combined gene set

Hierarchical clustering of the GIST samples was also performed using the combined gene set (AF-gene set, OVCA-gene set, and RCC-gene set) (Fig. 1d). Kaplan–Meier analysis was performed using the two sample sets defined by the first branch point. The probability of not developing a metastasis by 6 years was 0.59 for Group B, while none recurred in Group A (Fig. 1d; Table 1D, p = 0.029). In Group B, clustering was evident between two subgroups of sufficient sample size to analyze independently (Fig. 2). The probability of not developing a recurrence by 6 years was 0.39 for Group B2, while none of the patients in Group A or Group B1 recurred (Fig. 1e; Table 1E, p < 0.001 for comparisons between each of the 3 sets). We also grouped the samples into 2 groups defined as “good prognosis” or “poor prognosis”. Samples were defined as “good” prognosis if they were in Group A in the clustering by at least 2 of the 3 gene sets (AF-, OVCA-, or RCC; n = 30). The probability of not developing a metastasis by 6 years was 0.40 for the “poor” prognosis group, while none of the patients in the “good” prognosis group recurred (Fig. 1f; Table 1F, p < 0.001).

Fig. 2
figure2

Clustering of gene expression in the 60 GIST samples. The samples were clustered using the probes in the pooled gene set as described in the text. Groups A and B (= B1 + B2) are defined by the first branch point in the clustering. Groups B1 and B2 are defined by the second branch point in Group B

Full size image

Effect of Miettinen risk score on probability of developing recurrence

As some prognostic criteria correlate with recurrence of GIST, we questioned whether this scoring method might be improved by combining it with our clustering patterns. When the GIST samples were analyzed according to Miettinen risk status [11], none of the 29 patients in the low or very low risk groups recurred, yet 16 of 31 patients who scored in the high- or intermediate-risk score by Miettinen risk also did not recur (Fig. 3a; Table 2A).

Fig. 3
figure3

Kaplan–Meier analysis of the time to development of metastases. Kaplan–Meier analysis of the time to development of metastases for Miettinen risk groups of all GIST samples (a). Miettinen risk identified distinct risk groups (p < 0.001), panel a. When the 31 high- and intermediate-risk GIST samples from panel a were examined for their grouping using the individual AF-gene set (b), OVCA-gene set (c), RCC-gene set (d), and the combined gene set (e) in Fig. 1, the time to development of metastasis differed between groups A and B (b, p < 0.001; c, p = 0.003; d, p = 0.029; and e, p = 0.029)

Full size image
Table 2 Probability (95 % CI) of no metastasis as a function of time
Full size table

We went back to the hierarchical clustering performed with each of the gene sets in Fig. 1 to determine where these 31 patients with high- and intermediate-risk from the Miettinen score had been grouped, i.e. were they in Group A (good prognosis) or Group B (bad prognosis) (Fig. 1a–d). The probability of no metastasis for these 31 patients is shown in Fig. 3b–e and Table 2. Of interest is the finding that many of the 31 patients who were classified as high- or intermediate-risk by the Miettinen score were grouped as good prognosis using our gene sets and did not recur. The rate of recurrence in the “good risk” group was 0 % (0/8) for the AF-gene set, 14 % (2/14) for the OVCA-gene set, 18 % (3/17) for the RCC-gene set, and 0 % (0/5) for the pooled-gene set. Interestingly, among the “good” groups of high- and intermediate-risk samples defined by the AF-gene set, 0/3 high-risk and 0/5 intermediate-risk tumors recurred; 1/5 high- and 1/9 intermediate-risk tumors recurred in the “good” group identified by clustering with the OVCA-gene set, and 2/7 high- and 1/7 intermediate-risk tumors recurred when clustered with the RCC-gene set. Among the “bad” prognosis group defined by clustering with the AF-, OVCA-, and RCC- gene sets, 11/14 high- and 4/9 intermediate-risk, 10/12 high- and 3/5 intermediate-risk, and 9/10 high- and 3/7 intermediate-risk samples recurred, respectively.

Discussion

The biologic heterogeneity of GISTs, as with other soft tissue sarcomas, introduces complexities in deciding optimal treatment. This study used hierarchical clustering with gene sets derived from earlier studies of various biochemical pathways in aggressive fibromatosis, renal cell carcinoma, and ovarian carcinoma [36–39, 43] to examine 60 GIST samples using Agilent chip expression profiling. The analyses separated the GIST samples into at least two groups with different probabilities of developing metastatic disease. Although the gene sets were derived using biochemical pathways, we did not observe simple differences in biochemical pathways between the groups; possibly with a larger sample set, more detailed biochemical differences will become evident. Our data suggest that appreciation of these GIST subsets with distinct clinical behavior could be used to stratify GIST patients in clinical trials and in patient management. Miettinen risk group classification also identified distinct risk groups in our 60 GIST cases. In particular, our analysis also identified subsets of Miettinen high- and intermediate-risk samples that different in the risk of metastasis. When the high- and intermediate-risk GIST samples were examined without the low- and very low-risk samples, the individual AF-gene set, OVCA-gene set, RCC-gene set, and the combined gene set were associated with the time to development of metastasis. This finding suggests that further characterization of recurrence risk among samples classified as high- or intermediate-risk is possible. Furthermore, these results validate the potential role of the use of these gene sets in predicting the behavior of heterogeneous tumor sets.

These gene sets have also been shown to separate sets of soft tissue sarcoma samples into groups with different metastatic behavior [25, 26]. A gene set of 67 genes involved in mitosis and control of chromosome integrity, termed the complexity index in sarcomas (CINSARC), also predicts metastasis outcome in non-translocation dependent soft tissue sarcomas [23]. Both the gene sets used here and the CINSARC [23, 33] gene set identified subsets of the GIST samples that differed in time to recurrence. These data support the potential use of these gene sets to predict biological behavior in GIST as well as other soft tissue sarcomas. Only 11 of the 67 genes in the CINSARC gene set were also present in our pooled gene set.

Other methods of examining genetic heterogeneity may also be helpful. A recent study found that chromosomal changes detected by comparative genomic hybridization (CGH) were predictive of GIST outcome [33]. This study, as well as a second study, also found that a “genomic index” calculated from the number of chromosomal alterations (segmental gains and losses), and number of chromosomes involved was a strong predictor of recurrence as well [33, 44]. Another study using array-based analysis of gene copy number separated 42 GISTs into 4 groups with different survival rates [35].

Conclusions

Gene expression profiles may provide a useful technique to better predict long-term outcomes after surgery in patients with GIST and other sarcomas. Such information could be used to restrict the use of adjuvant therapy and reduce heterogeneity among groups in clinical trials. Due to the limited sample size of our study, we examined the identification of only two subsets of the GIST sample set with different metastatic propensity. The ability to detect multiple subgroups is highly dependent on the number of samples and the distribution of samples among the various groups. With larger sample sets, it may be possible to further refine classification and identify clinically useful heterogeneity. In addition, although gene expression analysis may provide a useful indicator of long-term outcomes, it should be used in combination with standard prognostic factors in order to have maximum predictive value [8]. For example, in this study, further characterization of recurrence risk among samples classified as high-or intermediate-risk was possible. These results also validate the potential role of the use of these gene sets in predicting the behavior of heterogeneous tumor sets. Several different gene sets appear to separate the samples into 2 groups with different behavior.

Abbreviations

GIST:

gastrointestinal stromal tumor

RCC:

clear cell renal cell carcinoma

OVCA:

ovarian cancer

AF:

aggressive fibromatosis

References

  1. 1.

    Corless CL, Barnett CM, Heinrich MC. Gastrointestinal stromal tumours: origin and molecular oncology. Nat Rev Cancer. 2011;11(12):865–78.

  2. 2.

    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–80.

  3. 3.

    Casali PG, Blay JY, ECECPO. Soft tissue sarcomas: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up. An Oncol. 2010;21(Suppl 5):v198–203.

  4. 4.

    Joensuu H, DeMatteo RP. The management of gastrointestinal stromal tumors: a model for targeted and multidisciplinary therapy of malignancy. Annu Rev Med. 2012;63:247–58.

  5. 5.

    Corless CL, Ballman KV, Antonescu CR, Kolesnikova V, Maki RG, Pisters PW, Blackstein ME, Blanke CD, Demetri GD, Heinrich MC, et al. Pathologic and molecular features correlate with long-term outcome after adjuvant therapy of resected primary GI stromal tumor: the ACOSOG Z9001 trial. J Clin Oncol. 2014;32(15):1563–70.

  6. 6.

    Dematteo RP, Ballman KV, Antonescu CR, Maki RG, Pisters PW, Demetri GD, Blackstein ME, Blanke CD, von Mehren M, Brennan MF, et al. Adjuvant imatinib mesylate after resection of localised, primary gastrointestinal stromal tumour: a randomised, double-blind, placebo-controlled trial. Lancet. 2009;373(9669):1097–104.

  7. 7.

    Joensuu H, Eriksson M, Sundby Hall K, Hartmann JT, Pink D, Schutte J, Ramadori G, Hohenberger P, Duyster J, Al-Batran SE, et al. One vs three years of adjuvant imatinib for operable gastrointestinal stromal tumor: a randomized trial. JAMA. 2012;307(12):1265–72.

  8. 8.

    Joensuu H, Rutkowski P, Nishida T, Steigen SE, Brabec P, Plank L, Nilsson B, Braconi C, Bordoni A, Magnusson MK, et al. KIT and PDGFRA Mutations and the Risk of GI Stromal Tumor Recurrence. J Clin Oncol. 2015;33(6):634–42.

  9. 9.

    Fletcher CD, Berman JJ, Corless C, Gorstein F, Lasota J, Longley BJ, Miettinen M, O’Leary TJ, Remotti H, Rubin BP, et al. Diagnosis of gastrointestinal stromal tumors: a consensus approach. Hum Pathol. 2002;33(5):459–65.

  10. 10.

    Joensuu H, Vehtari A, Riihimaki J, Nishida T, Steigen SE, Brabec P, Plank L, Nilsson B, Cirilli C, Braconi C, et al. Risk of recurrence of gastrointestinal stromal tumour after surgery: an analysis of pooled population-based cohorts. Lancet Oncol. 2012;13(3):265–74.

  11. 11.

    Miettinen M, Lasota J. Gastrointestinal stromal tumors: review on morphology, molecular pathology, prognosis, and differential diagnosis. Arch Pathol Lab Med. 2006;130(10):1466–78.

  12. 12.

    Patel S. Navigating risk stratification systems for the management of patients with GIST. Ann Surg Oncol. 2011;18(6):1698–704.

  13. 13.

    Rossi S, Miceli R, Messerini L, Bearzi I, Mazzoleni G, Capella C, Arrigoni G, Sonzogni A, Sidoni A, Toffolatti L, et al. Natural history of imatinib-naive GISTs: a retrospective analysis of 929 cases with long-term follow-up and development of a survival nomogram based on mitotic index and size as continuous variables. Am J Surg Pathol. 2011;35(11):1646–56.

  14. 14.

    Corless CL, Schroeder A, Griffith D, Town A, McGreevey L, Harrell P, Shiraga S, Bainbridge T, Morich J, Heinrich MC. PDGFRA mutations in gastrointestinal stromal tumors: frequency, spectrum and in vitro sensitivity to imatinib. J Clin Oncol. 2005;23(23):5357–64.

  15. 15.

    Agaram NP, Wong GC, Guo T, Maki RG, Singer S, Dematteo RP, Besmer P, Antonescu CR. Novel V600E BRAF mutations in imatinib-naive and imatinib-resistant gastrointestinal stromal tumors. Genes Chromosomes Cancer. 2008;47(10):853–9.

  16. 16.

    Andersson J, Sihto H, Meis-Kindblom JM, Joensuu H, Nupponen N, Kindblom LG. NF1-associated gastrointestinal stromal tumors have unique clinical, phenotypic, and genotypic characteristics. Am J Surg Pathol. 2005;29(9):1170–6.

  17. 17.

    Janeway KA, Kim SY, Lodish M, Nose V, Rustin P, Gaal J, Dahia PL, Liegl B, Ball ER, Raygada M, et al. Defects in succinate dehydrogenase in gastrointestinal stromal tumors lacking KIT and PDGFRA mutations. Proc Natl Acad Sci USA. 2011;108(1):314–8.

  18. 18.

    Killian JK, Kim SY, Miettinen M, Smith C, Merino M, Tsokos M, Quezado M, Smith WI Jr, Jahromi MS, Xekouki P, et al. Succinate dehydrogenase mutation underlies global epigenomic divergence in gastrointestinal stromal tumor. Cancer Discov. 2013;3(6):648–57.

  19. 19.

    Miettinen M, Wang ZF, Sarlomo-Rikala M, Osuch C, Rutkowski P, Lasota J. Succinate dehydrogenase-deficient GISTs: a clinicopathologic, immunohistochemical, and molecular genetic study of 66 gastric GISTs with predilection to young age. Am J Surg Pathol. 2011;35(11):1712–21.

  20. 20.

    Pantaleo MA, Lolli C, Nannini M, Astolfi A, Indio V, Saponara M, Urbini M, La Rovere S, Gill A, Goldstein D et al: Good survival outcome of metastatic SDH-deficient gastrointestinal stromal tumors harboring SDHA mutations. Genet Med. 2014;17(5):391–5.

  21. 21.

    Daniels M, Lurkin I, Pauli R, Erbstosser E, Hildebrandt U, Hellwig K, Zschille U, Luders P, Kruger G, Knolle J, et al. Spectrum of KIT/PDGFRA/BRAF mutations and Phosphatidylinositol-3-Kinase pathway gene alterations in gastrointestinal stromal tumors (GIST). Cancer Lett. 2011;312(1):43–54.

  22. 22.

    Miranda C, Nucifora M, Molinari F, Conca E, Anania MC, Bordoni A, Saletti P, Mazzucchelli L, Pilotti S, Pierotti MA, et al. KRAS and BRAF mutations predict primary resistance to imatinib in gastrointestinal stromal tumors. Clin Cancer Res. 2012;18(6):1769–76.

  23. 23.

    Chibon F, Lagarde P, Salas S, Perot G, Brouste V, Tirode F, Lucchesi C, de Reynies A, Kauffmann A, Bui B, et al. Validated prediction of clinical outcome in sarcomas and multiple types of cancer on the basis of a gene expression signature related to genome complexity. Nat Med. 2010;16(7):781–7.

  24. 24.

    Francis P, Namlos HM, Muller C, Eden P, Fernebro J, Berner JM, Bjerkehagen B, Akerman M, Bendahl PO, Isinger A, et al. Diagnostic and prognostic gene expression signatures in 177 soft tissue sarcomas: hypoxia-induced transcription profile signifies metastatic potential. BMC Genom. 2007;8:73.

  25. 25.

    Skubitz KM, Francis P, Skubitz AP, Luo X, Nilbert M. Gene expression identifies heterogeneity of metastatic propensity in high-grade soft tissue sarcomas. Cancer. 2012;118(17):4235–43.

  26. 26.

    Skubitz KM, Skubitz AP, Xu WW, Luo X, Lagarde P, Coindre JM, Chibon F. Gene expression identifies heterogeneity of metastatic behavior among high-grade non-translocation associated soft tissue sarcomas. J Transl Med. 2014;12:176.

  27. 27.

    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–74.

  28. 28.

    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–90.

  29. 29.

    Antonescu CR, Viale A, Sarran L, Tschernyavsky SJ, Gonen M, Segal NH, Maki RG, Socci ND, DeMatteo RP, Besmer P. Gene expression in gastrointestinal stromal tumors is distinguished by KIT genotype and anatomic site. Clin Cancer Res. 2004;10(10):3282–90.

  30. 30.

    Arne G, Kristiansson E, Nerman O, Kindblom LG, Ahlman H, Nilsson B, Nilsson O. Expression profiling of GIST: cD133 is associated with KIT exon 11 mutations, gastric location and poor prognosis. Int J Cancer. 2011;129(5):1149–61.

  31. 31.

    Bertucci F, Finetti P, Ostrowski J, Kim WK, Kim H, Pantaleo MA, Astolfi A, Polkowski M, Birnbaum D. Genomic Grade Index predicts postoperative clinical outcome of GIST. Br J Cancer. 2012;107(8):1433–41.

  32. 32.

    Koon N, Schneider-Stock R, Sarlomo-Rikala M, Lasota J, Smolkin M, Petroni G, Zaika A, Boltze C, Meyer F, Andersson L, et al. Molecular targets for tumour progression in gastrointestinal stromal tumours. Gut. 2004;53(2):235–40.

  33. 33.

    Lagarde P, Perot G, Kauffmann A, Brulard C, Dapremont V, Hostein I, Neuville A, Wozniak A, Sciot R, Schoffski P, et al. Mitotic checkpoints and chromosome instability are strong predictors of clinical outcome in gastrointestinal stromal tumors. Clin Cancer Res. 2012;18(3):826–38.

  34. 34.

    Yamaguchi U, Nakayama R, Honda K, Ichikawa H, Hasegawa T, Shitashige M, Ono M, Shoji A, Sakuma T, Kuwabara H, et al. Distinct gene expression-defined classes of gastrointestinal stromal tumor. J Clin Oncol. 2008;26(25):4100–8.

  35. 35.

    Ylipaa A, Hunt KK, Yang J, Lazar AJ, Torres KE, Lev DC, Nykter M, Pollock RE, Trent J, Zhang W. Integrative genomic characterization and a genomic staging system for gastrointestinal stromal tumors. Cancer. 2011;117(2):380–9.

  36. 36.

    Skubitz AP, Pambuccian SE, Argenta PA, Skubitz KM. Differential gene expression identifies subgroups of ovarian carcinoma. Transl Res. 2006;148(5):223–48.

  37. 37.

    Skubitz KM, Skubitz AP. Differential gene expression in renal-cell cancer. J Lab Clin Med. 2002;140(1):52–64.

  38. 38.

    Skubitz KM, Skubitz AP. Gene expression in aggressive fibromatosis. J Lab Clin Med. 2004;143(2):89–98.

  39. 39.

    Skubitz KM, Zimmermann W, Kammerer R, Pambuccian S, Skubitz AP. Differential gene expression identifies subgroups of renal cell carcinoma. J Lab Clin Med. 2006;147(5):250–67.

  40. 40.

    Eisen MB, Spellman PT, Brown PO, Botstein D. Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998;95(25):14863–8.

  41. 41.

    de Hoon MJ, Imoto S, Nolan J, Miyano S. Open source clustering software. Bioinformatics. 2004;20(9):1453–4.

  42. 42.

    Team RC. R: a language and environment for statistical computing. R Foundation for Statistical Computing. 2013.

  43. 43.

    Skubitz KM, Pambuccian S, Manivel JC, Skubitz AP. Identification of heterogeneity among soft tissue sarcomas by gene expression profiles from different tumors. J Transl Med. 2008;6:23.

  44. 44.

    Lartigue L, Neuville A, Lagarde P, Brulard C, Rutkowski P, Dei Tos P, Wardelmann E, Debiec-Rychter M, Italiano A, Coindre JM, et al. Genomic index predicts clinical outcome of intermediate-risk gastrointestinal stromal tumours, providing a new inclusion criterion for imatinib adjuvant therapy. Eur J Cancer. 2015;51(1):75–83.

Download references

Authors’ contributions

KMS, KG, APNS, WX participated in study design, data analysis, and helped draft the manuscript. JK participated in data analysis, and helped draft the manuscript. All authors read and approved the final manuscript.

Acknowledgements

Supported in part by the Karen Wyckoff Rein in Sarcoma Foundation, NIH P30 CA77598 for utilizing the support of the Biostatistics and Bioinformatics Core of the University of Minnesota Masonic Cancer Center shared resources, and the Minnesota Supercomputing Institute (MSI) at the University of Minnesota. We thank Dr. Frederick Chibon for generously providing the expression data and Michael Franklin for editorial assistance. Presented in part at the 20th annual meeting of the Connective Tissue Society, November 4–7, Salt Lake City, 2015.

Competing interests

The authors declare that they have no competing interests.

Author information

Affiliations

  1. Department of Medicine, The University of Minnesota Medical School, Minneapolis, MN, USA
    • Keith M. Skubitz
    •  & Kate Geschwind
  2. Masonic Cancer Center, The University of Minnesota Medical School, Minneapolis, MN, USA
    • Keith M. Skubitz
    • , Kate Geschwind
    • , Joseph S. Koopmeiners
    •  & Amy P. N. Skubitz
  3. Department of Biochemistry and Medical Genetics, Faculty of Medicine, University of Manitoba, The Research Institute of Oncology and Hematology, Cancer Care, Winnipeg, MA, Canada
    • Wayne W. Xu
  4. Division of Biostatistics, University of Minnesota School of Public Health, Minneapolis, USA
    • Joseph S. Koopmeiners
  5. Department of Laboratory Medicine and Pathology, University of Minnesota Medical School, Minneapolis, MN, USA
    • Amy P. N. Skubitz
Authors
  1. Search for Keith M. Skubitz in:
  2. Search for Kate Geschwind in:
  3. Search for Wayne W. Xu in:
  4. Search for Joseph S. Koopmeiners in:
  5. Search for Amy P. N. Skubitz in:

Corresponding author

Correspondence to Keith M. Skubitz.

Additional file

12967_2016_802_MOESM1_ESM.tif

Additional file 1: Figure S1. Clustering of gene expression in the 60 GIST samples. The samples were clustered using the probes in the indicated gene set as described in the text. Panel A, all probes; panel B, AF-gene set; panel C, OVCA-gene set; panel D, RCC-gene set.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Skubitz, K.M., Geschwind, K., Xu, W.W. et al. Gene expression identifies heterogeneity of metastatic behavior among gastrointestinal stromal tumors. J Transl Med 14, 51 (2016) doi:10.1186/s12967-016-0802-3

Download citation

Keywords

  • Microarray
  • Sarcoma
  • GIST
  • Gene expression
  • Heterogeneity
  • Subgroups
  • Metastasis
  • Prognosis

Journal of Translational Medicine

ISSN: 1479-5876

Contact us