- Open Access
Gastrointestinal stromal tumors: a case-only analysis of single nucleotide polymorphisms and somatic mutations
Clinical Sarcoma Research volume 3, Article number: 12 (2013)
Gastrointestinal stromal tumors are rare soft tissue sarcomas that typically develop from mesenchymal cells with acquired gain-in-function mutations in KIT or PDGFRA oncogenes. These somatic mutations have been well-characterized, but little is known about inherited genetic risk factors. Given evidence that certain susceptibility loci and carcinogens are associated with characteristic mutations in other cancers, we hypothesized that these signature KIT or PDGFRA mutations may be similarly fundamental to understanding gastrointestinal stromal tumor etiology. Therefore, we examined associations between 522 single nucleotide polymorphisms and seven KIT or PDGFRA tumor mutations types. Candidate pathways included dioxin response, toxin metabolism, matrix metalloproteinase production, and immune and inflammatory response.
We estimated odds ratios and 95% confidence intervals for associations between each candidate SNP and tumor mutation type in 279 individuals from a clinical trial of adjuvant imatinib mesylate. We used sequence kernel association tests to look for pathway-level associations.
One variant, rs1716 on ITGAE, was significantly associated with KIT exon 11 non-codon 557–8 deletions (odds ratio = 2.86, 95% confidence interval: 1.71-4.78) after adjustment for multiple comparisons. Other noteworthy associations included rs3024498 (IL10) and rs1050783 (F13A1) with PDGFRA mutations, rs2071888 (TAPBP) with wild type tumors and several matrix metalloproteinase SNPs with KIT exon 11 codon 557–558 deletions. Several pathways were strongly associated with somatic mutations in PDGFRA, including defense response (p = 0.005) and negative regulation of immune response (p = 0.01).
This exploratory analysis offers novel insights into gastrointestinal stromal tumor etiology and provides a starting point for future studies of genetic and environmental risk factors for the disease.
Gastrointestinal stromal tumors, or GISTs, are soft tissue sarcomas that develop from mesenchymal connective tissue anywhere in the gastrointestinal tract, though they most frequently appear in the stomach or small intestine [1–3]. Recent advances in molecular biology have revealed a distinct subset of these tumors that express tyrosine kinase receptors or platelet-derived growth factor receptors [4–6]. The presence of these receptors, encoded by the KIT or PDGFRA oncogenes, respectively, is an indication that these tumors share a common origin of the interstitial cells of Cajal, the pacemaker cells of the gut. The Cajal cells normally express both CD117, the immunohistochemical (IHC) marker for tyrosine kinase receptors, and CD34, the IHC marker for platelet-derived growth factor receptors, and GISTs likely develop from Cajal cells with acquired gain-in-function KIT or PDGFRA mutations [4, 7, 8]. KIT exon 11 (50-60% of cases) and KIT exon 9 (5-10% of cases) are the most common mutation sites [2, 9, 10]. Such mutations can enable receptor activation in the absence of normal stem cell factor signaling mechanisms, thereby over-stimulating cell proliferation and leading to tumor development. Based on this discovery, a 2001 National Institutes of Health consensus panel  agreed to formally define GISTs as mesenchymal neoplasms of the gastrointestinal tract displaying positivity for CD117 or CD34, with some exceptions allowed for immunonegative tumors with otherwise consistent histology.
These relatively new and complex diagnostic criteria make disease surveillance and etiologic study difficult. In a recent evaluation of patients in the National Cancer Institute’s Surveillance, Epidemiology and End Results database, Rubin et al.  estimated an annual age-adjusted incidence rate of 0.32 cases per 100,000 individuals in the United States (US). The rarity of the disease makes it a difficult subject for population-based research or for prompt and unbiased assessment of non-genetic risk factors in any study population. An evaluation of the genetic determinants of GIST is much more feasible, as the germline DNA of individuals does not change over time or in response to disease processes. To date, no other research groups have published such evaluations.
To help fill this knowledge gap, we investigated the genetic determinants of GIST in a case-only analysis. Specifically, we examined the associations between select single nucleotide polymorphisms (SNPs), which are inherited variations in individuals’ DNA, and the several common types of acquired KIT and PDGFRA mutations present in GIST tissue. Certain susceptibility loci have been linked to characteristic mutations, or mutational “signatures”, in other cancers. These include associations between GSTM1-null genotype and TP53 transversion mutations among bladder cancer patients , and certain functional polymorphisms in XPD and G:C→T:A TP53 mutations among lung cancer patients . Similarly, we hypothesized that the characteristic somatic mutations in the KIT and PDGFRA genes in GIST tumors may be mutational signatures that are causally linked to specific mutagens or susceptibility loci. As such, identifying risk factors for the individual tumor subtypes may be fundamental to understanding the disease.
We conducted our evaluation in two phases. The first phase included genes previously linked to soft tissue sarcoma or to environmental risk factors for soft tissue sarcoma, such as dioxins, phenoxyherbicides, insecticides, and vinyl chloride [14–18], as well as genes previously linked to mutational signatures in other cancers [12, 13, 19–21]. We found that several SNPs were associated with GIST tumor subtypes, including SNPs on two xenobiotic metabolizing genes, CYP1B1 and GSTM1, and two DNA repair genes, RAD23B and ERCC2.
The present report includes results from the second phase of the study, in which we examined the relationship between these 7 somatic mutation categories and 522 additional candidate SNPs. These SNPs are located on genes that play secondary or less well-understood roles in dioxin response or toxin-metabolizing pathways, as well as SNPs on PDGFRA and 10 matrix metalloproteinase (MMP) genes, which are often over-expressed in GISTs and other soft tissue sarcomas and may be linked to tumor invasion and metastasis [23–25]. Based on previous evidence that immune and inflammation-related genes are associated with osteosarcomas and other gastrointestinal tumors [26–28], we also selected certain SNPs from the GeneChip® Human Immune and Inflammation SNP Kit genotyping panel designed by Affymetrix® [29, 30].
Our main objective for these analyses was to identify genes or gene pathways potentially related to GIST carcinogenesis. Therefore, in addition to assessing the effect of each individual SNP, we also assessed the joint effects of SNPs in the same functional categories. By conducting these exploratory analyses with a large and diverse gene panel, our goal is to identify specific variants, genes, or functional pathways meriting further investigation.
The study population consists of the first 279 individuals from the American College of Surgeons Oncology Group (ACOSOG) Z9001 clinical trial who provided blood and tumor tissue samples for ancillary research and had sufficient tumor tissue for mutation analysis. The Z9001 trial was a multicenter, Phase III, randomized, double-blind study of adjuvant imatinib (Gleevec™; Novartis Pharmaceuticals) versus placebo. To be eligible for the clinical trial, cases had to have a CD117+, resected, localized, primary GIST of at least 3 cm diagnosed between July 1, 2002 and April 18, 2007. IHC staining for CD117 was completed using the Dako antibody (DakoCytomation, Copenhagen, Denmark). Additional information on the Z9001 trial is published elsewhere . Institutional Review Boards at all participating institutions approved this ancillary study and all of the included participants consented to the use of their blood and tissue specimens.
Among the target genes described above, we identified SNPs within 2000 or 500 base pairs of the 5′ and 3′ ends of the coding regions, respectively, that could affect gene function and had at least a 10% minor allele frequency in the HapMap CEU population . This included nonsense, missense and splice site mutations, as well as SNPs overlapping with microRNA seed regions or transcription binding sites. If the selected SNPs did not meet design phase quality control standards (designability score <1 or final score <0.7), we selected a surrogate SNP in high linkage disequilibrium with the desired SNP. SNPs related to dioxin response or toxin metabolism (final n = 68), matrix metalloproteinase (final n = 24), or on PDGFRA (final n = 4) were selected in this manner.
We selected the remaining SNPs (final n = 426) from a pre-existing Affymetrix® panel designed to include potentially functional non-synonymous SNPs from 318 genes related to immune and inflammatory response. Genes were selected for inclusion based on their gene ontology (GO) categorization . GO categories were also used to separate the genes into functional subgroups (see Additional file 1: Table S1).
DNA for mutation analysis was extracted from snap-frozen tumor tissue and tested for KIT exon 11 mutations using polymerase chain reaction (PCR) analysis (Platinum TaqDNA Polymerase High Fidelity; Life Technologies, Inc., Gaithersburg, MD). The PCR conditions were as follows: 1) 94°C for 4 min; 2) 94°C for 30 sec, 3) the relevant annealing temperature for each primer set for 30 sec, 4) 72°C for 30 sec, (35 cycles); and 5) 72°C for 3 min. The PCR products were identified by agarose gel electrophoresis using a 2% MetaPhor™ agarose gel (BioWhittaker Applications, Rockland, ME). The PCR products were purified with the QIAquick™ PCR Purification Kit (Qiagen Inc., Valencia, CA) before sequencing. The sequencing reactions for each case were performed from both the forward and reverse directions. Tumors lacking exon 11 mutations were genotyped for mutations in KIT exons 9, 13, 14 and 17 and PDGFRA exons 12 and 18. A more detailed description of these assays can be found elsewhere [34, 35].
An initial 544 candidate SNPs were genotyped using the GoldenGate genotyping assay (Illumina Inc., San Diego, CA). Briefly, allele-specific oligos were hybridized directly to genomic DNA extracted from the blood samples. The hybridized DNA was extended and ligated to downstream locus-specific oligos and then amplified using universal PCR fluorescently labeled primers and allele-specific primers [30, 36]. After the resulting products were hybridized to their complementary bead types, the arrays were assessed using the BeadArray™ Reader.
Twenty-seven participants underwent duplicate genotype analysis for quality assurance purposes. Concordance for duplicate samples was 99.9%. After excluding SNPs that were mono-allelic (n = 3), had >5% missing data (n = 3), or showed poor clustering (n = 16) among our study subjects, we had 522 evaluable SNPs, as listed above. We retained three SNPs that showed evidence of possible copy number variation, but designated them accordingly.
We conducted descriptive analyses of selected demographic variables, tumor characteristics and genotypes. As this population includes some non-white participants, we calculated race-specific MAFs and compared genotype distributions across racial groups (white vs. other) using a Pearson χ2 test of association. Fisher’s exact test was used if one or more cells had less than 5 observations.
We categorized each individual’s tumor based on the presence or absence of each of the following outcomes: i) a deletion of KIT exon 11 codons 557–558, ii) any other (i.e. non-codon 557–8) deletion in KIT exon 11, iii) a KIT exon 11 insertion, iv) a KIT exon 11 point mutation, v) a KIT exon 9, exon 13, exon 14, or exon 17 mutation, vi) a PDGFRA exon 12 or 18 mutation, and vii) no KIT or PDGFRA mutation (wild type). KIT mutations in exons 9, 13, 14 and 17 were too rare for independent evaluation.
We obtained odds ratios (ORs), 95% confidence intervals (CIs) and p-values for each SNP-mutation combination using logistic regression. All models were adjusted for race, sex, and age at diagnosis. We assumed additive genetic models, denoting whether an individual had 0, 1 or 2 copies of the minor allele. P-values were calculated using trend tests and were corrected for multiple testing by controlling for a false discovery rate of 25%. This method is less conservative than a Bonferroni approach and is thus better suited for a hypothesis-generating study such as this .
We used a sequence kernel association test (SKAT) to assess the joint effect of a group of SNPs on an outcome [38, 39]. We grouped SNPs according to their functional category, as described above. This method is well-powered to detect associations when SNPs in a group are correlated with one another but only moderately associated with the outcome. Briefly, for individual i, the log odds of having the outcome given genotypes zi1 to zip and covariates xi1 to xim is modeled semi-parametrically using a logistic kernel-machine regression model:
where h(Z i ) is a function of a positive, definite kernel function, K(•,•) and some γi, …, γn:
We used an identity-by-state kernel: , as this does not require linearity assumptions and allows for epistasis . Assuming h follows an arbitrary distribution with a mean of 0 and variance τK, testing the null hypothesis H0: h(Z) = 0 is equivalent to testing H0: τ=0. This is accomplished using a modified variance-component score statistic:
Here, Q is comparable to a χ2 distribution with scale parameter κ and ν degrees of freedom, both of which are modified to account for correlation between SNPs in the same SNP-set (for calculations, see Appendix A in Wu et al. ).
Males and females were approximately equally represented in our study population, while 18% were non-white (Table 1). Median age at diagnosis was 58 years (range 18–85), though non-white participants tended to be younger (median age = 53 years). Non-white participants were also more likely to have a smaller tumor than white participants (median tumor sizes of 6.0 cm versus 6.5 cm), and more likely to have stomach tumors (74% versus 64%). Most patients had mutations in KIT exon 11 (70% overall), the largest proportion of which were codon 557–558 deletions (34% of exon 11 mutations). Demographic and tumor characteristics were very similar for males and females.
Compared with the larger ACOSOG Z9001 cohort, whites were somewhat over-represented in this ancillary study, which had otherwise similar characteristics (data not shown). Race-stratified MAF and association p-values for all 522 SNPs are displayed in Additional file 1: Table S2. As expected based on genotype distributions in ethnically diverse HapMap populations , genotype distributions in this study differed by race for many of the candidate polymorphisms.
The top 5 SNPs for each mutation subtype are displayed in Table 2. Only one SNP, rs1716 on ITGAE, was statistically significant after adjusting for multiple comparisons. This SNP was associated with KIT exon 11 non-codon 557–8 deletions (OR = 2.86, 95% CI: 1.71, 4.78; p = 6.4×10-5).
Though not statistically significant after adjusting for multiple comparisons, rs3024498 (IL10) and rs1050783 (F13A1) were strongly associated with PDGFRA mutations (OR = 0.31, 95% CI: 0.16,0.60 and OR = 0.31, 95% CI: 0.16,0.61, respectively) and rs2071888 (TAPBP) was strongly associated with wild type tumors (OR = 0.37, 95% CI: 0.20, 0.67). Additionally, several SNPs in matrix metalloproteinase genes were associated with tumor subtypes. MMP10 and MMP1 SNPs were associated with KIT exon 11 codon 557–558 deletions, 2 MMP7 SNPs were associated with KIT exon 11 point mutations and a MMP1 SNP was associated with KIT exon 9, 13, 14, or 17 mutations. Effect estimates and p-values for all 522 SNPs can be found in Additional file 1: Table S3. The relative magnitude of all SNP-subtype associations is depicted in Figure 1.
Despite strong SNP-level effects, the MMP pathway was not associated with any of the tumor subtypes in the SKAT analyses (minimum p-value =0.2 for KIT exon 11 point mutations; Additional file 1: Table S4). As seen in Table 3, the strongest pathway-level associations were in relation to somatic mutations in PDGFRA. This included defense response (p = 0.005), negative regulation of immune response (p = 0.01), protein phosphorylation (p = 0.02), positive regulation of immune response (p = 0.03), and AHR/dioxin response (0.03). Additionally, negative regulation of cell proliferation was associated with PDGFRA mutations (p = 0.04; Additional file 1: Table S4). In total, only 5 other pathways were associated with a tumor subtype at p < 0.05. These were AHR/dioxin response with non exon 11 KIT mutations (p = 0.01), humoral immune response with wild type mutations (p = 0.02), and response to stress, negative regulation of apoptosis, and protein tyrosine kinase activity with non-codon 557–8 KIT exon 11 deletions (p = 0.02, p = 0.03, and p = 0.04, respectively). No pathways were statistically significant after correcting for multiple testing. Log p-values for all pathway analyses can be seen in Figure 2.
In this exploratory analysis of genetic risk factors for GIST tumor subtypes, we identified one statistically significant association and a number of other potentially important associations for individual polymorphisms. We also identified several potentially relevant functional pathways. These novel findings offer clues about the etiology of these rare and poorly understood tumors.
The SNP with the strongest association with a tumor subtype was rs1716 on ITGAE. This SNP results in a missense mutation on ITGAE (also known as CD103), a gene involved in protein tyrosine phosphatase activity. This SNP was previously associated with increased risk of melanoma , as was another SNP in the gene. The CD103 protein is commonly expressed in intraepithelial lymphocytic T cells and hairy cell leukemia cells [41, 42].
The IL10 SNP associated with PDGFRA mutations, rs3024498, is located in a seed microRNA region. One previous study found an association between rs3024498 and colorectal cancer . The IL10 gene encodes a cytokine that plays a role in immunoregulation and inflammation, and has been previously linked to several cancers, including osteosarcoma , cervical cancer , and gastric cancer [45, 46].
rs1050783 in F13A1 is also in a seed microRNA region, but neither the SNP nor the gene has been previously linked to cancer. The same is true for rs2071888 in TAPBP, a missense mutation. As noted above, several studies have observed over-expression of matrix metalloproteinase genes in GISTs and other soft tissue sarcomas [23–25], though none of the evaluated SNPs have previously been associated with cancer risk.
Pathway analyses are common in cancer epidemiology, but there is little consistency in how the pathways are defined. We selected the well-documented and publicly available Gene Ontology  classification system to facilitate replication and follow-up studies. Although we did not identify any studies that examined the specific pathways included in the present analyses, numerous studies have observed associations of inflammatory or immune response genes with risk of sarcomas or gastrointestinal cancers [26, 28, 30, 47].
The present exploratory study was undertaken in light of evidence suggesting that some mutagens and susceptibility loci are associated with specific mutational “signatures,” i.e., characteristic mutation patterns [12, 13, 19–21]. At this point, GIST etiology has been insufficiently researched, and we do not know how well the mutation-based tumor classifications used in this study correspond to distinct carcinogenic processes. However, the existence of multiple types of mutations suggest that more than one mutagenic process could be involved, and we believe that identifying associations between germline genetic polymorphisms and unique tumor phenotypes could contribute valuable new information about disease etiology. This information could also help to elucidate environmental risk factors for this disease.
Because tumors with KIT exon 11 mutations were not assessed for other KIT or PDGFRA mutations, some tumors may be misclassified, though evidence from population-based studies suggests that few GISTs have more than one mutation type [2, 9]. Our study participants had similar mutation profiles to the individuals included in these population-based investigations, but the results from this predominately white clinical trials population may not be generalizable to all GIST patients. Lastly, this study had a small sample size. As such, we had limited power to detect true associations, particularly when the evaluated genotype and mutation type were rare.
In this novel study of genetic risk factors for GIST, we identified several SNPs and gene pathways associated with GIST mutation subtypes. This included SNPs involved in dioxin response, toxin metabolism, matrix metalloproteinase synthesis, and inflammatory or immune response. While only a single SNP was statistically significant after correcting for multiple comparisons, our overall findings provide an important starting point for future studies of genetic and environmental risk factors for this rare and poorly understood disease.
Gastrointestinal stromal tumors
Single nucleotide polymorphisms
Sequence kernel association test
American college of surgeons oncology group
Polymerase chain reaction
Minor allele frequency.
Gastrointestinal stromal tumors treatment (PDQ®) - national cancer institute.http://www.cancer.gov/cancertopics/pdq/treatment/gist/HealthProfessional,
Mazzola P, Spitale A, Banfi S, Mazzucchelli L, Frattini M, Bordoni A: Epidemiology and molecular biology of gastrointestinal stromal tumors (GISTs): a population-based study in the South of Switzerland, 1999–2005. Histol Histopathol. 2008, 23: 1379-1386.
Miranda C, Nucifora M, Molinari F, Conca E, Anania MC, Bordoni A, Saletti P, Mazzucchelli L, Pilotti S, Pierotti MA, Tamborini E, Greco A, Frattini M: KRAS and BRAF mutations predict primary resistance to imatinib in gastrointestinal stromal tumors. Clin Cancer Res. 2012, 18: 1769-1776. 10.1158/1078-0432.CCR-11-2230
Hirota S, Isozaki K, Moriyama Y, Hashimoto K, Nishida T, Ishiguro S, Kawano K, Hanada M, Kurata A, Takeda M, Tunio GM, Matsuzawa Y, Kanakura Y, Shinomura Y, Kitamura Y: Gain-of-function mutations of c-kit in human gastrointestinal stromal tumors. Science. 1998, 279: 577-580. 10.1126/science.279.5350.577
Nilsson B, Bümming P, Meis-Kindblom JM, Odén A, Dortok A, Gustavsson B, Sablinska K, Kindblom LG: Gastrointestinal stromal tumors: the incidence, prevalence, clinical course, and prognostication in the preimatinib mesylate era–a population-based study in western Sweden. Cancer. 2005, 103: 821-829. 10.1002/cncr.20862
Yan BM, Kaplan GG, Urbanski S, Nash CL, Beck PL: Epidemiology of gastrointestinal stromal tumors in a defined Canadian health region: a population-based study. Int J Surg Pathol. 2008, 16: 241-250. 10.1177/1066896907306967
Rubin BP, Heinrich MC, Corless CL: Gastrointestinal stromal tumour. Lancet. 2007, 369: 1731-1741. 10.1016/S0140-6736(07)60780-6
Berman J, O’Leary TJ: Gastrointestinal stromal tumor workshop. Hum Pathol. 2001, 32: 578-582. 10.1053/hupa.2001.25484
Cassier PA, Ducimetière F, Lurkin A, Ranchère-Vince D, Scoazec J, Bringuier P, Decouvelaere AV, Méeus P, Cellier D, Blay JY, Ray-Coquard I: A prospective epidemiological study of new incident GISTs during two consecutive years in Rhône alpes region: incidence and molecular distribution of GIST in a European region. Br J Cancer. 2010, 103: 165-170. 10.1038/sj.bjc.6605743
Steigen SE, Eide TJ: Trends in incidence and survival of mesenchymal neoplasm of the digestive tract within a defined population of northern Norway. APMIS. 2006, 114: 192-200. 10.1111/j.1600-0463.2006.apm_261.x
Rubin JL, Sanon M, Taylor DC, Coombs J, Bollu V, Sirulnik L: Epidemiology, survival, and costs of localized gastrointestinal stromal tumors. Int J Gen Med. 2011, 4: 121-130.
Ryk C, Berggren P, Kumar R, Hemminki K, Larsson P, Steineck G, Lambert B, Hou SM: Influence of GSTM1, GSTT1, GSTP1 and NAT2 genotypes on the p53 mutational spectrum in bladder tumours. Int J Cancer. 2005, 113: 761-768. 10.1002/ijc.20650
Mechanic LE, Marrogi AJ, Welsh JA, Bowman ED, Khan MA, Enewold L, Zheng YL, Chanock S, Shields PG, Harris CC: Polymorphisms in XPD and TP53 and mutation in human lung cancer. Carcinogenesis. 2005, 26: 597-604.
Kogevinas M, Becher H, Benn T, Bertazzi PA, Boffetta P, Bueno-de-Mesquita HB, Coggon D, Colin D, Flesch-Janys D, Fingerhut M, Green L, Kauppinen T, Littorin M, Lynge E, Matthews JD, Neuberger M, Pearce N, Saracci R: Cancer mortality in workers exposed to phenoxy herbicides, chlorophenols, and dioxins: an expanded and updated international cohort study. Am J Epidemiol. 1997, 145: 1061-1075. 10.1093/oxfordjournals.aje.a009069
Boffetta P, Matisane L, Mundt KA, Dell LD: Meta-analysis of studies of occupational exposure to vinyl chloride in relation to cancer mortality. Scand J Work Environ Health. 2003, 29: 220-229. 10.5271/sjweh.725
Hardell L: Pesticides, soft-tissue sarcoma and non-hodgkin lymphoma–historical aspects on the precautionary principle in cancer prevention. Acta Oncol. 2008, 47: 347-354. 10.1080/02841860701753697
Pahwa P, Karunanayake CP, Dosman JA, Spinelli JJ, McLaughlin JR, : Soft-tissue sarcoma and pesticides exposure in men: results of a Canadian case–control study. J Occup Environ Med. 2011, 53: 1279-1286. 10.1097/JOM.0b013e3182307845
Brady MS, Gaynor JJ, Brennan MF: Radiation-associated sarcoma of bone and soft tissue. Arch Surg. 1992, 127: 1379-1385. 10.1001/archsurg.1992.01420120013002
Pfeifer GP, Denissenko MF: Formation and repair of DNA lesions in the p53 gene: relation to cancer mutations?. Environ Mol Mutagen. 1998, 31: 197-205. 10.1002/(SICI)1098-2280(1998)31:3<197::AID-EM1>3.0.CO;2-I
Olivier M, Hussain SP, Caron De Fromentel C, Hainaut P, Harris CC: TP53 mutation spectra and load: a tool for generating hypotheses on the etiology of cancer. IARC Sci Publ. 2004, 157: 247-270.
Dixon K, Kopras E: Genetic alterations and DNA repair in human carcinogenesis. Semin Cancer Biol. 2004, 14: 441-448. 10.1016/j.semcancer.2004.06.007
O’Brien KM, Orlow I, Antonescu CR, Ballman KV, McCall L, DeMatteo RP, Engel LS: Gastrointestinal stromal tumors, somatic mutations and candidate genetic risk variants. PLoS One. 2013, 8: e62119- 10.1371/journal.pone.0062119
Alford SH, Vrana MS, Waite L, Heim-Hall J, Sylvia VL, Williams RP: Matrix metalloproteinase expression in high grade soft tissue sarcomas. Oncol Rep. 2007, 18: 1529-1536.
Nabeshima K, Iwasaki H, Nishio J, Koga K, Shishime M, Kikuchi M: Expression of emmprin and matrix metalloproteinases (MMPs) in peripheral nerve sheath tumors: emmprin and membrane-type (MT)1-MMP expressions are associated with malignant potential. Anticancer Res. 2006, 26: 1359-1367.
Nakagawa M, Nabeshima K, Asano S, Hamasaki M, Uesugi N, Tani H, Yamashita Y, Iwasaki H: Up-regulated expression of ADAM17 in gastrointestinal stromal tumors: coexpression with EGFR and EGFR ligands. Cancer Sci. 2009, 100: 654-662. 10.1111/j.1349-7006.2009.01089.x
Oliveira ID, Petrilli AS, Tavela MH, Zago MA, de Toledo SRC: TNF-alpha, TNF-beta, IL-6, IL-10, PECAM-1 and the MPO inflammatory gene polymorphisms in osteosarcoma. J Pediatr Hematol Oncol. 2007, 29: 293-297. 10.1097/MPH.0b013e3180587e69
Walczak A, Przybylowska K, Dziki L, Sygut A, Chojnacki C, Chojnacki J, Dziki A, Majsterek I: The lL-8 and IL-13 gene polymorphisms in inflammatory bowel disease and colorectal cancer. DNA Cell Biol. 2012, 31: 1431-1438. 10.1089/dna.2012.1692
Hussain SK, Mu L, Cai L, Chang S, Park SL, Oh SS, Wang Y, Goldstein BY, Ding BG, Jiang Q, Rao J, You NCY, Yu SZ, Papp JC, Zhao JK, Wang H, Zhang ZF: Genetic variation in immune regulation and DNA repair pathways and stomach cancer in China. Cancer Epidemiol Biomarkers Prev. 2009, 18: 2304-2309. 10.1158/1055-9965.EPI-09-0233
Cerhan JR, Ansell SM, Fredericksen ZS, Kay NE, Liebow M, Call TG, Dogan A, Cunningham JM, Wang AH, Liu-Mares W, Macon WR, Jelinek D, Witzig TE, Habermann TM, Slager SL: Genetic variation in 1253 immune and inflammation genes and risk of non-hodgkin lymphoma. Blood. 2007, 110: 4455-4463. 10.1182/blood-2007-05-088682
Dematteo RP, Ballman KV, Antonescu CR, Maki RG, Pisters PWT, Demetri GD, Blackstein ME, Blanke CD, von Mehren M, Brennan MF, Patel S, McCarter MD, Polikoff JA, Tan BR, Owzar K: Adjuvant imatinib mesylate after resection of localised, primary gastrointestinal stromal tumour: a randomised, double-blind, placebo-controlled trial. Lancet. 2009, 373: 1097-1104. 10.1016/S0140-6736(09)60500-6
The gene ontology.http://www.geneontology.org/,
Antonescu CR, Sommer G, Sarran L, Tschernyavsky SJ, Riedel E, Woodruff JM, Robson M, Maki R, Brennan MF, Ladanyi M, DeMatteo RP, Besmer P: Association of KIT exon 9 mutations with nongastric primary site and aggressive behavior: KIT mutation analysis and clinical correlates of 120 gastrointestinal stromal tumors. Clin Cancer Res. 2003, 9: 3329-3337.
Dematteo RP, Gold JS, Saran L, Gönen M, Liau KH, Maki RG, Singer S, Besmer P, Brennan MF, Antonescu CR: Tumor mitotic rate, size, and location independently predict recurrence after resection of primary gastrointestinal stromal tumor (GIST). Cancer. 2008, 112: 608-615. 10.1002/cncr.23199
Shen R, Fan J, Campbell D, Chang W, Chen J, Doucet D, Yeakley J, Bibikova M, Garcia EW, McBride C, Steemers F, Garcia F, Kermani BG, Gunderson K, Oliphant A: High-throughput SNP genotyping on universal bead arrays. Mutat Res. 2005, 573: 70-82. 10.1016/j.mrfmmm.2004.07.022
Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I: Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001, 125: 279-284. 10.1016/S0166-4328(01)00297-2
Wu MC, Kraft P, Epstein MP, Taylor DM, Chanock SJ, Hunter DJ, Lin X: Powerful SNP-set analysis for case–control genome-wide association studies. Am J Hum Genet. 2010, 86: 929-942. 10.1016/j.ajhg.2010.05.002
Wu MC, Lee S, Cai T, Li Y, Boehnke M, Lin X: Rare-variant association testing for sequencing data with the sequence kernel association test. Am J Hum Genet. 2011, 89: 82-93. 10.1016/j.ajhg.2011.05.029
Lenci RE, Rachakonda PS, Kubarenko AV, Weber ANR, Brandt A, Gast A, Sucker A, Hemminki K, Schadendorf D, Kumar R: Integrin genes and susceptibility to human melanoma. Mutagenesis. 2012, 27: 367-373. 10.1093/mutage/ger090
Anz D, Mueller W, Golic M, Kunz WG, Rapp M, Koelzer VH, Ellermeier J, Ellwart JW, Schnurr M, Bourquin C, Endres S: CD103 is a hallmark of tumor-infiltrating regulatory T cells. Int J Cancer. 2011, 129: 2417-2426. 10.1002/ijc.25902
Venkataraman G, Aguhar C, Kreitman RJ, Yuan CM, Stetler-Stevenson M: Characteristic CD103 and CD123 expression pattern defines hairy cell leukemia: usefulness of CD123 and CD103 in the diagnosis of mature B-cell lymphoproliferative disorders. Am J Clin Pathol. 2011, 136: 625-630. 10.1309/AJCPKUM9J4IXCWEU
Tsilidis KK, Helzlsouer KJ, Smith MW, Grinberg V, Hoffman-Bolton J, Clipp SL, Visvanathan K, Platz EA: Association of common polymorphisms in IL10, and in other genes related to inflammatory response and obesity with colorectal cancer. Cancer Causes Control. 2009, 20: 1739-1751. 10.1007/s10552-009-9427-7
Ni J, Ye Y, Teng F, Wu Q: Interleukin 10 polymorphisms and cervical cancer risk: a meta-analysis. Int J Gynecol Cancer. 2013, 23: 126-133. 10.1097/IGC.0b013e318274b1a2
Xue H, Wang Y, Lin B, An J, Chen L, Chen J, Fang JY: A meta-analysis of interleukin-10–592 promoter polymorphism associated with gastric cancer risk. PLoS One. 2012, 7: e39868- 10.1371/journal.pone.0039868
Pan F, Tian J, Pan Y, Zhang Y: Association of IL-10-1082 promoter polymorphism with susceptibility to gastric cancer: evidence from 22 case–control studies. Mol Biol Rep. 2012, 39: 7143-7154. 10.1007/s11033-012-1546-7
Castro FA, Koshiol J, Hsing AW, Gao Y, Rashid A, Chu LW: Inflammatory gene variants and the risk of biliary tract cancers and stones: a population-based study in China. BMC Cancer. 2012, 12: 468- 10.1186/1471-2407-12-468
We thank Aliaksandra Samoila and Kenneth Cheung for their technical assistance during the preparation of DNA specimens for genotyping and direct sequencing. The authors also wish to thank Michael Wu for his statistical advice. Novartis Pharmaceuticals Corporation provided funding for this research project with agreement #CSTI571BUS249.
This work was supported by Novartis Pharmaceuticals Corp. agreement #CSTI571BUS249. RD previously served as a paid consultant to Novartis Pharmaceuticals. All other authors declare that they have no conflict of interest.
LSE, IO, CRA and RD conceived the study design, selected the candidate genes, and supervised genetic analyses. KMO and LSE designed the statistical analysis and drafted the manuscript. KMO and LM preformed the statistical analyses, with contributions from LSE and KB. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Table S1: Genes and gene pathways. Table S2. Minor allele frequencies (MAF) and p-values for comparison of genotype frequencies: Z9001 genotyped whites (n=273) versus non-whites (n=58). Table S3. Odds Ratios (ORs), 95% Confidence Intervals (Cis) and p-values for the association between candidate SNPs and tumor mutation status in 279 GIST patients. Table S4. P-values for sequence kernel association tests (SKAT) of functional pathways and tumor mutation status. (DOCX 523 KB)
About this article
Cite this article
O’Brien, K.M., Orlow, I., Antonescu, C.R. et al. Gastrointestinal stromal tumors: a case-only analysis of single nucleotide polymorphisms and somatic mutations. Clin Sarcoma Res 3, 12 (2013). https://doi.org/10.1186/2045-3329-3-12
- Gastrointestinal stromal tumors
- Somatic mutations
- Single nucleotide polymorphisms