Skip to main content
  • Research article
  • Open access
  • Published:

Somatic loss of function mutations in neurofibromin 1 and MYC associated factor X genes identified by exome-wide sequencing in a wild-type GIST case

Abstract

Background

Approximately 10–15 % of gastrointestinal stromal tumors (GISTs) lack gain of function mutations in the KIT and platelet-derived growth factor receptor alpha (PDGFRA) genes. An alternate mechanism of oncogenesis through loss of function of the succinate-dehydrogenase (SDH) enzyme complex has been identified for a subset of these “wild type” GISTs.

Methods

Paired tumor and normal DNA from an SDH-intact wild-type GIST case was subjected to whole exome sequencing to identify the pathogenic mechanism(s) in this tumor. Selected findings were further investigated in panels of GIST tumors through Sanger DNA sequencing, quantitative real-time PCR, and immunohistochemical approaches.

Results

A hemizygous frameshift mutation (p.His2261Leufs*4), in the neurofibromin 1 (NF1) gene was identified in the patient’s GIST; however, no germline NF1 mutation was found. A somatic frameshift mutation (p.Lys54Argfs*31) in the MYC associated factor X (MAX) gene was also identified. Immunohistochemical analysis for MAX on a large panel of GISTs identified loss of MAX expression in the MAX-mutated GIST and in a subset of mainly KIT-mutated tumors.

Conclusion

This study suggests that inactivating NF1 mutations outside the context of neurofibromatosis may be the oncogenic mechanism for a subset of sporadic GIST. In addition, loss of function mutation of the MAX gene was identified for the first time in GIST, and a broader role for MAX in GIST progression was suggested.

Peer Review reports

Background

Gastrointestinal stromal tumor (GIST) is a mesenchymal neoplasm that originates throughout the GI tract, primarily in the stomach (>50 %) and small intestine (~30 %) [1]. GIST generally presents in older adults, while ~2 % of cases are children [2, 3]. Originally thought to be of smooth muscle origin, immunohistochemical and ultrastructural studies suggest that GIST is related to spindle-shaped cells of the GI tract known as the interstitial cells of Cajal (ICC) [4, 5]. ICC and the majority (95 %) of GIST express the type III receptor tyrosine kinase KIT (CD117), and variably exhibit myoid or neural features. The majority of GISTs exhibit gain of function mutations in KIT or in the related receptor PDGFRA [6, 7]. A subset (~10–15 %) of GISTs in adults lack mutations in the KIT and PDGFRA genes, as do almost all pediatric cases [8, 9]. The commonly used label of “wild type” (WT) GIST belies the epidemiological, clinico-pathological and molecular heterogeneity that define these tumors. WT GIST occurs in the context of several multitumor syndromes, including the inherited Carney-Stratakis Syndrome (CSS) and the non-familial Carney triad (CT). Manifestations of CSS and CT include gastric GIST and paraganglioma (PGL), a neural crest-derived tumor, while the spectrum of CT neoplasms includes pulmonary chondroma and several other neoplasms [10, 11]. CSS results from loss of function mutations in subunit genes of the succinate-dehydrogenase (SDH) enzyme complex, SDHB, SDHC, and SDHD [12]. Inactivation of the SDHA gene subunit has recently been implicated as a causative factor in a subset of apparently sporadic adult WT GIST (reviewed in [13]). GISTs from CT patients do not manifest SDHX mutations; however, these tumors are also SDH-deficient, and the molecular underpinning of CT GIST has been attributed to epigenetic silencing of the SDHC gene [14]. Pediatric GIST patients share hallmarks of CT-associated GIST, namely early-onset, multi-focal, gastric disease with a predilection towards females [8], and pediatric GIST cases have also been associated with SDHC epimutation [15]. Thus the identification of SDH-deficient GIST, also referred to as “type 2”, helps distinguish between KIT/PDGFRA mutant, or type I GIST, and a majority of wild type GIST.

SDH-intact WT GISTs with alternate oncogenic events have been described. Mutations in the serine-threonine kinase gene BRAF have recently been identified in approximately 5–15 % of sporadic wild type GIST [16]. These tumors are generally KIT-positive with spindle cell or mixed morphology, and are found primarily in the small intestine in adult cases. Approximately 1–2 % of GISTs occur in the context of neurofibromatosis type I (NF1) [1], an autosomal dominant disorder with skin and ophthalmologic manifestations that predisposes to a variety of benign and malignant tumors. GIST in NF1 individuals also present typically in the small bowel with spindle-cell morphology, are found in men and women at a younger median age than KIT/PDGFRA mutant GIST, and are often multifocal [17, 18]. Neurofibromatosis is due to germline mutations in neurofibromin 1, a RAS-GAP protein and negative regulator for RAS signaling, and germline NF1 mutations accompanied by somatic events have been identified in NF1 GIST cases [19].

In this report we describe whole exome sequencing (WES) of a particularly complex, SDH-intact wild type GIST case. The WES analysis identified for the first time the somatic inactivation of NF1 in a GIST outside the context of NF1 syndrome. A novel somatic loss of function mutation in the MYC-associated factor X (MAX) gene was also identified. Immunohistochemical studies of a panel of GISTs identified deficiencies in MAX expression in a number of tumors. Implications for these and other identified mutations are discussed.

Methods

Preparation of genomic DNA and total RNA

De-identified tumor samples and normal blood were obtained following written informed consent from the Fox Chase Cancer Center Biosample Repository. The protocol was approved by the Fox Chase Cancer Center Institutional Review Board (#03-848). The isolation and characterization of genomic DNA and total RNA from frozen tumor specimens embedded in optimum-temperature cutting medium has been described [20].

Whole exome sequencing data analysis

Exome-enriched genomic libraries (TruSeq, Illumina, San Diaego CA) from normal and tumor DNA were subjected to paired-end 100 bp sequencing on the Illumina HiSeq 2000 instrument. Reads were mapped to the 1000Genome Project reference human genome (Hg19 corresponding v37) using the BWA aligner [21] and mapped reads were sorted, merged, and de-duplicated (Picard), yielding an average of 51.6 million unique mapped reads per sample. GATK realignment was used to realign reads locally in areas surrounding insertions and deletions (indels) [22, 23]. Variant calling and filtering was performed using GATK UnifiedGenotyper [22, 23] and single nucleotide variants (SNVs) annotated with modified ANNOVAR [24]. This pipeline yielded an average SNV rate of ~ 0.34 % per sample. The downstream analysis of SNVs and indels was done by custom Perl scripts. Non-synonymous, potentially deleterious coding region variants, splice-site mutations, and indels that were predicted to be present in the tumor only, were visually confirmed on the Integrative Genomics Viewer (IGV) [25], and confirmed by exon-based Sanger sequencing. Confirmed somatic indels, and deleterious missense mutations predicted by the SIFT algorithm [26] and confirmed by a consensus approach [27] are listed in Table 1. Mutation nomenclature conforms to the recommendations of the Human Genome Variation Society [28].

Table 1 Confirmed somatic mutations

Sanger sequencing

Primers for amplification and sequencing of KIT (exons 9, 11, 13, 17), PDGFRA (exons 12, 14, 18), and BRAF (exons 11,15) have been described [29], as have primers for SDHA [30] and SDHB-D [31]. Primer sequences for confirmation of mutations listed in Table 1 and MAX genomic sequencing are shown in Additional file 1: Table S1. Relevant exons were PCR-amplified from genomic DNA and subjected to Sanger sequencing (Beckman Coulter Genomics).

Immunohistochemical analysis

GIST tissue microarrays (TMAs) were constructed in conjunction with the FCCC Biosample Repository. H&E-stained sections from paraffin-embedded tissue blocks were evaluated by a pathologist for tumor content and cellularity, and two cores from each block were selected for the TMA. Each TMA consists of ~ 30 GIST specimens along with normal tissue sections. IHC for MAX was performed with the SC-197 antibody (Santa Cruz Biotechnology, Dallas TX) at a 1:400 dilution with antigen retrieval. Aperio Digital Pathology (Leica Biosystems, Buffalo Grove, IL) was used to capture and quantify MAX-stained TMAs using the nuclear algorithm. MAX-deficient cases were confirmed on whole-tissue sections, as were a subset of MAX-positive cases. IHC for KIT was performed as described [32].

Gene expression analysis

Random-primed cDNA was prepared from 2 μg total RNA using the High Capacity cDNA Reverse Transcription KIT (Life Technologies). RNA expression was measured by real-time PCR (qRT-PCR) on an ABI PRISM 7900 HT Sequence Detection System using fluorescein phosphoramidite (FAM) primer/probe sets (Applied Biosystems). RNA expression data for MAX were normalized using hypoxanthine guanine phosphoribosyl transferase 1 (HPRT1) and glucuronidase beta (GUSB). Taqmen sets used were Hs99999909_m1 (HPRT1), Hs99999908_m1 (GUSB), and Hs00811069_g1 (MAX).

Results

The patient first presented at the age of 54 with a high-risk GIST of the small bowel. Resections of several local and distant recurrences were documented over several years, and the patient was treated with imatinib, sunitinib, and several additional targeted agents. A locally recurrent 1.5 cm small bowel tumor was resected, and a small portion of flash-frozen tumor and a whole blood sample were provided following informed consent. Formalin-fixed paraffin embedded (FFPE) tissue from an earlier small bowel resection was also available. Histologic evaluation of these specimens indicated a high-cellularity spindle-cell tumor with a high mitotic index (5-10/50 HPF) and no necrosis.

Sanger sequencing of DNA from the flash-frozen tumor indicated that no mutations were present in the “hotspot” exons of KIT, PDGFRA, or BRAF. Exon-based sequencing of the SDH complex subunit genes SDHA-D identified no SDHX mutations, and the tumor was immunohistochemically positive for SDHB expression [20]. These analyses suggested that the case belonged to the small subset of SDH-intact and KIT, PDGFA, BRAF wild-type GIST, for which no clear molecular pathogenic mechanism has been established. DNA from this GIST and from the patient’s blood was therefore subjected to whole-exome sequencing (WES) (see Materials and Methods).

WES analysis identified a two-base insertion (c.6781_6782insTT; p.His2240Leufs*4) in exon 44 in the Neurofibromatosis type I gene (NF1) that was confirmed by Sanger sequencing (Fig. 1a and Table 1). The mutation is not seen in the patient’s germline DNA, and the wild type allele is not represented in the tumor in either the WES or Sanger analysis. A previously reported single-nucleotide polymorphism (SNP) array analysis of this GIST (case 26 [29]) identified copy-number loss of the region encompassing the NF1 gene locus, suggesting somatic NF1 gene inactivation through the frameshift mutation combined with loss of the wild type gene. This particular NF1 mutation was not found in the COSMIC [33], ClinVar [34] or Leiden Open Variation databases [35]. To our knowledge this is the first reported example of GIST with an inactivating NF1 mutation outside the context of the NF1 syndrome. This sporadic GIST does share certain characteristics with GISTs from NF1 patients, such as small bowel origin, spindled-cell morphology, and immunopositivity for KIT and SDHB [36]. It is reasonable to suggest that somatic NF1 gene inactivation may be a causative factor in formation of the patient’s disease.

Fig. 1
figure 1

a WES (top) and Sanger (bottom) sequencing showing the two-base (TT) insertion in NF1. A subset of reads visualized on the Integrative Genomics Viewer (IGV) shows the insertion represented by the purple bar in 100 % of the tumor reads, as confirmed by the chromatogram below. b Top and bottom panels show the single-base (c) deletion in the MAX gene in a majority (~90 %) of reads, again confirmed by the chromatogram below. Red arrow indicates direction of transcription for MAX

WES analysis also identified a truncating frameshift mutation (c.160delC; p.Gln54Lysfs*10) in exon 3 of the MYC-associated factor (MAX) gene in the tumor DNA (Fig. 1b and Table 1). This mutation was found in 91 % of the WES reads. MAX is a basic helix-loop-helix (H-L-H) leucine zipper (LZ) transcription factor and a key member of the MYC/MAX/MXD network [37]. This truncating mutation is predicted to disrupt domains responsible for MAX homo-dimerization and hetero-dimerization [38]. Inactivating MAX mutations have recently been reported in inherited and sporadic PGL and pheochromocytomas (PCC) [39, 40], and in small cell lung cancer (SCLC) specimens [41]. No mutations in MAX have been reported in GIST, and we did not find additional MAX mutations in a sample set of 16 wild-type tumors from 11 patients. However, these WT GIST are all SDH-negative tumors, and exhibit other characteristics (e.g. gastric location, epithelioid cell morphology, lack of genome complexity) [20] that are not found in our index case. An earlier report documented a significant reduction in MAX expression in association with copy-number loss surrounding the MAX gene locus in a set of kinase-mutant GISTs [42]. We hypothesized that reduction or loss of MAX expression may be associated with mutant GIST. To test this hypothesis, immunohistochemistry from MAX was carried out on a series of ~80 GIST specimens contained on 3 GIST tissue microarrays (TMAs). The antibody was first tested against human seminal vesicle tissue, which exhibited strong nuclear staining (Fig. 2, panel a). Staining of the GIST from the patient’s earlier resection confirmed complete absence of MAX and strong plasma membrane staining for KIT in the tumor (Fig. 2, panels b and c respectively). IHC analysis of the 3 TMAs identified a wide range of MAX staining in GIST sections, ranging from strong, widely distributed nuclear staining to complete or near-complete absence of staining. Images were captured and quantified (Materials and Methods). Nuclear staining intensity (0–3) and distribution (0–100) were combined to generate nuclear H-Scores, with a potential range of 0–300. The mean nuclear H-score across the GIST samples on the TMA was 99.7 (0–252). Visual re-examination of the stained spots suggested H-scores <20 as a reasonable cutoff below which very little positive nuclear staining was seen. Whole tissue sections from a number of these tumors were re-stained with MAX, and the TMA results confirmed in 10 of these cases (Table 2 and Fig. 2). Several GISTs that exhibited strong or intermediate nuclear MAX staining are shown in Fig. 2 (panels d-f, with H-scores of 210, 192, and 84.8, respectively). In Fig. 2, panels g-p correspond to Cases 1–10 in Table 2. These images generally corroborate the low nuclear MAX H-scores from the TMA, although some tissues show isolated nuclear staining (e.g. panel h).

Fig. 2
figure 2

Immunohistochemistry for MAX and KIT. a Control staining of nuclear MAX in seminal vesicles. b Negative staining for MAX in index case. c Strong plasma membrane straining for KIT/CD117 in index case. d-f Strong nuclear staining for MAX in MAX-positive GISTs. g-p Mainly negative nuclear staining for MAX in GIST cases 1–10. Red bar: 10 μM

Table 2 Description of MAX-negative cases

Table 2 lists the clinico-pathological characteristics of the 10 MAX-deficient cases. Nine cases harbor KIT exon 11 mutations, while one patient’s tumor exhibited a mutation in exon 18 of PDGFRA. The cases are mainly males (7 of 9 available), with an average age of presentation of 56.9 years. Most were of gastric origin, and nine were high- or intermediate-risk GIST. None of these parameters varied significantly however when compared to the rest of the sample set (Additional file 2: Table S2). Exon-based Sanger sequencing of eight of these cases did not identify any mutations in MAX.

Table 1 lists several other confirmed somatic mutations in our index case: two indels and five missense mutations that were predicted to affect protein structure or function. The variant calls represent ~30–60 % of the total reads at these positions, and exon-based Sanger sequencing confirmed that these mutations were at most heterozygous in the tumors. A frameshift deletion in reticulon 4 (RTN4) was identified, along with an in-frame insertion in the coiled-coiled domain-containing protein CCDC66. Missense mutations were identified in 5 genes: the gene encoding the enzyme mevalonate decarboxylase (MVD) that catalyzes an early step in cholesterol synthesis; MAFA (musculoaponeurotic fibrosarcoma oncogene family, protein A), a transcription factor that controls insulin gene expression in the pancreas; the ring-finger protein gene RNF123, which acts as a ubiquitin ligase towards the cyclin-dependent kinase inhibitor KIP1; a member of the spindlin family of chromatin readers (SPIN4); and the SELP (selectin P) gene, encoding a calcium-dependent receptor that mediates the interaction of activated endothelial cells or platelets with leukocytes.

Discussion

This study is the first report of somatic inactivation of the NF1 gene in a sporadic SDH-intact GIST lacking gain of function mutations in KIT, PDGFR, and BRAF. Although tissue from the primary tumor from this patient was not available for analysis, it is reasonable to suggest that somatic inactivation of NF1 may have been an early causative event in this case. NF1 patients are ~45X more likely to develop GIST than the general population [17], and it has been estimated that 1–2 % of GIST arise in patients with NF1 [1]. GISTs in NF1 patients commonly lack activating mutations in the KIT and PDGFRA receptors [17, 18, 43], and may owe their incidence to germline NF1 mutation coupled with somatic second hits, as has been demonstrated in some cases [19, 44, 45]. It is perhaps not surprising to find NF1 gene inactivation in a sporadic wild-type GIST, as NF1 somatic mutations have been identified in a number of non-NF1-associated tumor types (reviewed in [44]). Moreover, while PCCs are known to occur in the context of NF1 [46], NF1 somatic mutations were also identified in a high percentage (21/61) of sporadic PCC selected for specific gene expression patterns or low levels of NF1 gene expression [47]. The finding of NF1 gene inactivation in sporadic GIST has diagnostic implications, as the molecular identification of mutations in the large and complex NF1 gene is a challenging task. A comprehensive approach combining NF1 transcript and genomic sequencing with multiplex ligation dependent probe amplification and other techniques for detection of gene duplications and deletions has been used to detect mutations in up to 95% of NF1 cases [48]. Immunohistochemical approaches using available anti-NF1 antibodies have been largely unsuccessful in identifying NF1-mutated PCC with a high degree of sensitivity or specificity [47, 49]. Recently, WES approaches have been used to identify the germline and somatic NF1 events in various tumors from an NF1 patient [50]. Whether NF1 inactivation is a common event in sporadic, SDH-intact wild type GIST is an open question. In a recent transcriptome-sequencing study [51], no NF1 mutations were identified in two SDH intact wild-type GIST, although it has been shown that only a portion of exonic variants are typically captured by RNA-seq approaches [52]. Interestingly, next-generation sequencing of eight SDH-negative GIST cases using a targeted cancer-associated gene capture library identified a low-frequency (8 %) frameshift NF1 mutation in a GIST that also harbored an activating KRAS gene mutation (G12V) [53].

The identification of a loss of function mutation in the MAX gene is a novel finding in GIST. As a hetero-dimeric partner for MYC, MAX was originally thought to be required for oncogenic pathways initiated by MYC over-expression [54]. However, in MAX-deficient rat PC12 PCC cell lines it has been shown that MAX is dispensable for MYC transcriptional regulation [55]. Intriguingly, MAX inactivation has recently been implicated in both inherited and sporadic PCC and PGL cases [39, 40]. Mutation of MAX is a relatively rare event in these tumors, accounting for 1.12 % of hereditary PCC/PGL lacking mutations in other susceptibility genes and 1.65 % of sporadic cases. In the familial PCC cases preferential transmission of the disease from the paternal allele was observed, along with a tendency towards aggressive behavior. A recent report also describes MAX inactivation in 6 % of primary SCLC specimens [41]. The authors suggest that SCLC, like PCC, may arise from neuroendocrine cells or differentiate towards neural features, which may explain the shared mechanism of MAX-associated oncogenesis. Similarly, sub-populations of GISTs may also variably exhibit neural properties or markers, and in a recent report a set of SDH-intact WT GIST was shown to exhibit high relative expression of neural markers, along with expression of members of the insulin-like growth family network [56]. In SCLC, inactivating MAX mutations were found to be mutually exclusive with amplifications of hetero-dimeric partners MYC, MYCL1, and MYCN, and mutations in BRG1, which encodes an ATPase of the SWI/SNF chromatin -remodeling complex that regulates expression of MYC, MYC target genes, as well as MAX. In GIST, the contributions of the MYC/MAX/MXD network to pathogenesis have not been extensively described. There have been descriptions of amplification of the MYC gene locus on chromosome 8q [57, 58], and reduced mRNA expression of the MAX gene associated with copy number loss of chromosome14q [42]. These secondary chromosome aberrations are common in KIT/PDGFRA-mutated GIST [59]. In this report we used immunohistochemical approaches to identify reduced/absent MAX nuclear staining in 10/78 (~13 %) of GIST cases analyzed, in addition to the index case. We found no additional MAX mutations in these tumors, and MAX RNA expression was only marginally and not significantly reduced (1.3-fold, P = 0.47) compared to the remaining MAX-positive cases. Further investigations into the mechanism(s) of MAX dysregulation and its contribution to pathogenesis in GIST are warranted.

In addition to inactivating NF1 and MAX mutations found in our index case, we also identified a heterozygous 4-base deletion in RTN4 and an in-frame insertion in CCD66, as well as potentially pathogenic substitutions in the MVD, MAFA, RNF123, SPIN4, and SELP genes. Although to our knowledge these genes have not been studied in GIST, a recent systems biological approach to identifying key transcriptional regulators in GIST and leiomyosarcoma identified nine differentially expressed genes, including the MYC gene, the MAF gene which encodes another basic leucine zipper (bZIP)-containing transcription factor closely related to MAFA, and another coiled-coil domain containing transcription factor gene, CCDC6 [60]. The MAF proteins are members of the AP1 family: the large MAF proteins contain transactivation domains and are considered onco-proteins by virtue of their ability to transform primary cells and induce tumors in various animal models (reviewed in [61]). Interestingly, the Q257L substitution we identified is located in the MAFA bZIP domain, very close to the predicted DNA-binding domain that has been shown to be required for MAFA transformation activity in avian fibroblasts [62]. The substitution we identified may affect the specificity or avidity of MAFA binding to its target sequences.

Conclusions

In conclusion, next-generation sequencing of an SDH-intact, KIT, PDGFR, BRAF wild type GIST identified for the first time somatic loss of function mutations in two tumor-suppressor genes, NF1 and MAX. Somatic inactivation of neurofibromin should be explored as a potential oncogenic mechanism in this subset of GIST. The identification of MAX inactivation provides another etiological link between GIST and PGL/PCC, in addition to mutations of the NF1 gene and mutations in the subunit genes of the SDH complex that have been identified in these tumors.

Availability of supporting data

The data supporting the results of this article are included within the article and its additional files.

Abbreviations

CSS:

Carney-Stratakis syndrome

CT:

Carney triad

GIST:

gastrointestinal stromal tumors

H&E:

hematoxylin and eosin

ICC:

interstitial cells of Cajal

IHC:

immunohistochemistry

Indels:

insertions and deletions

MAX:

Myc-associated factor X

NF1:

neurofibromatosis type 1

PDGFRA:

platelet-derived growth factor receptor alpha

PGL:

paraganglioma

PCC:

pheochromocytoma

qRT-PCR:

quantitative real-time polymerase chain reaction

SCLC:

small-cell lung cancer

SDH:

succinate dehydrogenase

SNV:

single nucleotide variant

SNP:

single nucleotide polymorphism

TMA:

tissue microarray

WES:

whole exome sequencing

References

  1. Miettinen M, Lasota J. Gastrointestinal stromal tumors. Gastroenterol Clin North Am. 2013;42(2):399–415.

    Article  PubMed Central  PubMed  Google Scholar 

  2. Miettinen M, Lasota J, Sobin LH. Gastrointestinal stromal tumors of the stomach in children and young adults: a clinicopathologic, immunohistochemical, and molecular genetic study of 44 cases with long-term follow-up and review of the literature. Am J Surg Pathol. 2005;29(10):1373–81.

    Article  PubMed  Google Scholar 

  3. Prakash S, Sarran L, Socci N, DeMatteo RP, Eisenstat J, Greco AM, et al. Gastrointestinal stromal tumors in children and young adults: a clinicopathologic, molecular, and genomic study of 15 cases and review of the literature. J Pediatr Hematol Oncol. 2005;27(4):179–87.

    Article  PubMed  Google Scholar 

  4. 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–69.

    PubMed Central  CAS  PubMed  Google Scholar 

  5. Sircar K, Hewlett BR, Huizinga JD, Chorneyko K, Berezin I, Riddell RH. Interstitial cells of Cajal as precursors of gastrointestinal stromal tumors. Am J Surg Pathol. 1999;23(4):377–89.

    Article  CAS  PubMed  Google Scholar 

  6. Hirota S, Isozaki K, Moriyama Y, Hashimoto K, Nishida T, Ishiguro S, et al. Gain-of-function mutations of c-kit in human gastrointestinal stromal tumors. Science. 1998;279(5350):577–80.

    Article  CAS  PubMed  Google Scholar 

  7. Heinrich MC, Corless CL, Duensing A, McGreevey L, Chen CJ, Joseph N, et al. PDGFRA activating mutations in gastrointestinal stromal tumors. Science. 2003;299(5607):708–10.

    Article  CAS  PubMed  Google Scholar 

  8. Rink L, Godwin AK. Clinical and molecular characteristics of gastrointestinal stromal tumors in the pediatric and young adult population. Curr Oncol Rep. 2009;11(4):314–21.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Pappo AS, Janeway KA. Pediatric gastrointestinal stromal tumors. Hematol Oncol Clin North Am. 2009;23(1):15–34. vii.

    Article  PubMed  Google Scholar 

  10. Stratakis CA, Carney JA. The triad of paragangliomas, gastric stromal tumours and pulmonary chondromas (Carney triad), and the dyad of paragangliomas and gastric stromal sarcomas (Carney-Stratakis syndrome): molecular genetics and clinical implications. J Intern Med. 2009;266(1):43–52.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Zhang L, Smyrk TC, Young Jr WF, Stratakis CA, Carney JA. Gastric stromal tumors in Carney triad are different clinically, pathologically, and behaviorally from sporadic gastric gastrointestinal stromal tumors: findings in 104 cases. Am J Surg Pathol. 2010;34(1):53–64.

    Article  PubMed Central  PubMed  Google Scholar 

  12. Pasini B, McWhinney SR, Bei T, Matyakhina L, Stergiopoulos S, Muchow M, et al. Clinical and molecular genetics of patients with the Carney-Stratakis syndrome and germline mutations of the genes coding for the succinate dehydrogenase subunits SDHB, SDHC, and SDHD. Eur J Hum Genet. 2008;16(1):79–88.

    Article  CAS  PubMed  Google Scholar 

  13. Belinsky MG, Rink L, von Mehren M. Succinate dehydrogenase deficiency in pediatric and adult gastrointestinal stromal tumors. Front Oncol. 2013;3:117.

    Article  PubMed Central  PubMed  Google Scholar 

  14. Haller F, Moskalev EA, Faucz FR, Barthelmess S, Wiemann S, Bieg M, et al. Aberrant DNA hypermethylation of SDHC: a novel mechanism of tumor development in Carney triad. Endocr Relat Cancer. 2014;21(4):567–77.

    Article  CAS  PubMed  Google Scholar 

  15. Killian JK, Miettinen M, Walker RL, Wang Y, Zhu YJ, Waterfall JJ, et al. Recurrent epimutation of SDHC in gastrointestinal stromal tumors. Sci Transl Med. 2014;6(268):268ra177.

    Article  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Miettinen M, Fetsch JF, Sobin LH, Lasota J. Gastrointestinal stromal tumors in patients with neurofibromatosis 1: a clinicopathologic and molecular genetic study of 45 cases. Am J Surg Pathol. 2006;30(1):90–6.

    Article  PubMed  Google Scholar 

  18. 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.

    Article  PubMed  Google Scholar 

  19. Maertens O, Prenen H, Debiec-Rychter M, Wozniak A, Sciot R, Pauwels P, et al. Molecular pathogenesis of multiple gastrointestinal stromal tumors in NF1 patients. Hum Mol Genet. 2006;15(6):1015–23.

    Article  CAS  PubMed  Google Scholar 

  20. Belinsky MG, Rink L, Flieder DB, Jahromi MS, Schiffman JD, Godwin AK, et al. Overexpression of insulin-like growth factor 1 receptor and frequent mutational inactivation of SDHA in wild-type SDHB-negative gastrointestinal stromal tumors. Genes Chromosomes Cancer. 2012.

  21. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25(14):1754–60.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011;43(5):491–8.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38(16), e164.

    Article  PubMed Central  PubMed  Google Scholar 

  25. Robinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, Getz G, et al. Integrative genomics viewer. Nat Biotechnol. 2011;29(1):24–6.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Kumar P, Henikoff S, Ng PC. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat Protoc. 2009;4(7):1073–81.

    Article  CAS  PubMed  Google Scholar 

  27. Kumar S, Sanderford M, Gray VE, Ye J, Liu L. Evolutionary diagnosis method for variants in personal exomes. Nat Methods. 2012;9(9):855–6.

    Article  CAS  PubMed  Google Scholar 

  28. den Dunnen JT, Antonarakis SE. Mutation nomenclature extensions and suggestions to describe complex mutations: a discussion. Hum Mutat. 2000;15(1):7–12.

    Article  Google Scholar 

  29. Belinsky MG, Skorobogatko YV, Rink L, Pei J, Cai KQ, Vanderveer LA, et al. High density DNA array analysis reveals distinct genomic profiles in a subset of gastrointestinal stromal tumors. Genes Chromosomes Cancer. 2009;48(10):886–96.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Burnichon N, Briere JJ, Libe R, Vescovo L, Riviere J, Tissier F, et al. SDHA is a tumor suppressor gene causing paraganglioma. Hum Mol Genet. 2010;19(15):3011–20.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Wood LD, Parsons DW, Jones S, Lin J, Sjoblom T, Leary RJ, et al. The genomic landscapes of human breast and colorectal cancers. Science. 2007;318(5853):1108–13.

    Article  CAS  PubMed  Google Scholar 

  32. Tarn C, Rink L, Merkel E, Flieder D, Pathak H, Koumbi D, et al. Insulin-like growth factor 1 receptor is a potential therapeutic target for gastrointestinal stromal tumors. Proc Natl Acad Sci U S A. 2008;105(24):8387–92.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Forbes SA, Tang G, Bindal N, Bamford S, Dawson E, Cole C, et al. COSMIC (the Catalogue of Somatic Mutations in Cancer): a resource to investigate acquired mutations in human cancer. Nucleic Acids Res. 2010;38(Database issue):D652–7.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Landrum MJ, Lee JM, Riley GR, Jang W, Rubinstein WS, Church DM, et al. ClinVar: public archive of relationships among sequence variation and human phenotype. Nucleic Acids Res. 2014;42(Database issue):D980–5.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Fokkema IF, den Dunnen JT, Taschner PE. LOVD: easy creation of a locus-specific sequence variation database using an “LSDB-in-a-box” approach. Hum Mutat. 2005;26(2):63–8.

    Article  CAS  PubMed  Google Scholar 

  36. Wang JH, Lasota J, Miettinen M. Succinate Dehydrogenase Subunit B (SDHB) Is Expressed in Neurofibromatosis 1-Associated Gastrointestinal Stromal Tumors (Gists): Implications for the SDHB Expression Based Classification of Gists. J Cancer. 2011;2:90–3.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Diolaiti D, McFerrin L, Carroll PA, Eisenman RN. Functional interactions among members of the MAX and MLX transcriptional network during oncogenesis. Biochim Biophys Acta. 2014;1849(5):484–500.

    Article  PubMed  Google Scholar 

  38. Reddy CD, Dasgupta P, Saikumar P, Dudek H, Rauscher 3rd FJ, Reddy EP. Mutational analysis of Max: role of basic, helix-loop-helix/leucine zipper domains in DNA binding, dimerization and regulation of Myc-mediated transcriptional activation. Oncogene. 1992;7(10):2085–92.

    CAS  PubMed  Google Scholar 

  39. Burnichon N, Cascon A, Schiavi F, Morales NP, Comino-Mendez I, Abermil N, et al. MAX mutations cause hereditary and sporadic pheochromocytoma and paraganglioma. Clin Cancer Res. 2012;18(10):2828–37.

    Article  CAS  PubMed  Google Scholar 

  40. Comino-Mendez I, Gracia-Aznarez FJ, Schiavi F, Landa I, Leandro-Garcia LJ, Leton R, et al. Exome sequencing identifies MAX mutations as a cause of hereditary pheochromocytoma. Nat Genet. 2011;43(7):663–7.

    Article  CAS  PubMed  Google Scholar 

  41. Romero OA, Torres-Diz M, Pros E, Savola S, Gomez A, Moran S, et al. MAX inactivation in small cell lung cancer disrupts MYC-SWI/SNF programs and is synthetic lethal with BRG1. Cancer Discov. 2014;4(3):292–303.

    Article  CAS  PubMed  Google Scholar 

  42. Wozniak A, Sciot R, Guillou L, Pauwels P, Wasag B, Stul M, et al. Array CGH analysis in primary gastrointestinal stromal tumors: cytogenetic profile correlates with anatomic site and tumor aggressiveness, irrespective of mutational status. Genes, Chromosomes & Cancer. 2007;46(3):261–76.

    Article  CAS  Google Scholar 

  43. Kinoshita K, Hirota S, Isozaki K, Ohashi A, Nishida T, Kitamura Y, et al. Absence of c-kit gene mutations in gastrointestinal stromal tumours from neurofibromatosis type 1 patients. J Pathol. 2004;202(1):80–5.

    Article  CAS  PubMed  Google Scholar 

  44. Laycock-van Spyk S, Thomas N, Cooper DN, Upadhyaya M. Neurofibromatosis type 1-associated tumours: their somatic mutational spectrum and pathogenesis. Hum Genomics. 2011;5(6):623–90.

    Article  PubMed  Google Scholar 

  45. Stewart DR, Corless CL, Rubin BP, Heinrich MC, Messiaen LM, Kessler LJ, et al. Mitotic recombination as evidence of alternative pathogenesis of gastrointestinal stromal tumours in neurofibromatosis type 1. J Med Genet. 2007;44(1), e61.

    Article  PubMed Central  PubMed  Google Scholar 

  46. Walther MM, Herring J, Enquist E, Keiser HR, Linehan WM. von Recklinghausen’s disease and pheochromocytomas. J Urol. 1999;162(5):1582–6.

    Article  CAS  PubMed  Google Scholar 

  47. Burnichon N, Buffet A, Parfait B, Letouze E, Laurendeau I, Loriot C, et al. Somatic NF1 inactivation is a frequent event in sporadic pheochromocytoma. Hum Mol Genet. 2012;21(26):5397–405.

    Article  CAS  PubMed  Google Scholar 

  48. Messiaen LM, Wimmer K: NF1 mutational spectrum. In: Neurofibromatoses. Edited by D K, vol. 16; 2008: 63–77.

  49. Stenman A, Svahn F, Welander J, Gustavson B, Soderkvist P, Gimm O, Juhlin CC: Immunohistochemical NF1 analysis does not predict NF1 gene mutation status in pheochromocytoma. Endocr Pathol 2015, 26(1):9–14.

    Article  CAS  PubMed  Google Scholar 

  50. Hirbe A, Dahiya SM, Miller C, Li T, Fulton R, Zhang X, et al. Whole exome sequencing reveals the order of genetic changes during malignant transformation and metastasis in a single patient with NF1-plexiform neurofibroma. Clin Cancer Res. 2015.

  51. Nannini M, Astolfi A, Urbini M, Indio V, Santini D, Heinrich MC, et al. Integrated genomic study of quadruple-WT GIST (KIT/PDGFRA/SDH/RAS pathway wild-type GIST). BMC Cancer. 2014;14:685.

    Article  PubMed Central  PubMed  Google Scholar 

  52. Cirulli ET, Singh A, Shianna KV, Ge D, Smith JP, Maia JM, et al. Screening the human exome: a comparison of whole genome and whole transcriptome sequencing. Genome Biol. 2010;11(5):R57.

    Article  PubMed Central  PubMed  Google Scholar 

  53. Hechtman JF, Zehir A, Mitchell T, Borsu L, Singer S, Tap W, et al. Novel oncogene and tumor suppressor mutations in KIT and PDGFRA wild type gastrointestinal stromal tumors revealed by next generation sequencing. Genes Chromosomes Cancer. 2014;54(3):177–84.

    Article  PubMed  Google Scholar 

  54. Amati B, Brooks MW, Levy N, Littlewood TD, Evan GI, Land H. Oncogenic activity of the c-Myc protein requires dimerization with Max. Cell. 1993;72(2):233–45.

    Article  CAS  PubMed  Google Scholar 

  55. Ribon V, Leff T, Saltiel AR. c-Myc does not require max for transcriptional activity in PC-12 cells. Mol Cell Neurosci. 1994;5(3):277–82.

    Article  CAS  PubMed  Google Scholar 

  56. Beadling C, Patterson J, Justusson E, Nelson D, Pantaleo MA, Hornick JL, et al. Gene expression of the IGF pathway family distinguishes subsets of gastrointestinal stromal tumors wild type for KIT and PDGFRA. Cancer Med. 2013;2(1):21–31.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  57. Yang J, Du X, Lazar AJ, Pollock R, Hunt K, Chen K, et al. Genetic aberrations of gastrointestinal stromal tumors. Cancer. 2008;113(7):1532–43.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  58. Assamaki R, Sarlomo-Rikala M, Lopez-Guerrero JA, Lasota J, Andersson LC, Llombart-Bosch A, et al. Array comparative genomic hybridization analysis of chromosomal imbalances and their target genes in gastrointestinal stromal tumors. Genes, Chromosomes & Cancer. 2007;46(6):564–76.

    Article  CAS  Google Scholar 

  59. Gunawan B, von Heydebreck A, Sander B, Schulten HJ, Haller F, Langer C, et al. An oncogenetic tree model in gastrointestinal stromal tumours (GISTs) identifies different pathways of cytogenetic evolution with prognostic implications. J Pathol. 2007;211(4):463–70.

    Article  CAS  PubMed  Google Scholar 

  60. Ylipaa A, Yli-Harja O, Zhang W, Nykter M. A systems biological approach to identify key transcription factors and their genomic neighborhoods in human sarcomas. Chin J Cancer. 2011;30(1):27–40.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  61. Eychene A, Rocques N, Pouponnot C. A new MAFia in cancer. Nat Rev Cancer. 2008;8(9):683–93.

    Article  CAS  PubMed  Google Scholar 

  62. Kataoka K, Nishizawa M, Kawai S. Structure-function analysis of the maf oncogene product, a member of the b-Zip protein family. J Virol. 1993;67(4):2133–41.

    PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to acknowledge Dr. Mark Andrake of the FCCC Molecular Modeling facility, along with the FCCC Biosample Repository, Histopathology Facility, Genotyping and Real-Time PCR facility, and the DNA Sequencing Facility for work contributing to this manuscript. The authors acknowledge support from the Kansas Bioscience Authority Eminent Scholar Program (awarded to A.K. Godwin). Dr. Godwin is the Chancellors Distinguished Chair in Biomedical Sciences endowment at University of Kansas. The authors also acknowledge support from the GIST Cancer Research Fund and the D’Amico family fund.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Martin G. Belinsky.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

MB: contributed to design and execution of the study and drafted the manuscript; LR: contributed to execution of the study and involved in drafting the manuscript; KQ: carried out IHC analysis and involved in drafting the manuscript; SC: performed validation work for WES analysis; YH and JC: carried out the WES bioinformatics analysis and data interpretation; AG: contributed to design and execution of the study and involved in revising the manuscript; MvM: contributed to design and execution of the study; involved in drafting the manuscript, gave approval of the final manuscript. All authors read and approved the final manuscript.

Additional files

Additional file 1: Table S1.

Primers for Sanger sequencing. Nucleotide sequences are listed for primers used for exon-based validation of somatic mutations listed in Table 1, and for all exons of the MAX gene. (DOC 42 kb)

Additional file 2: Table S2.

Characteristics of MAX-negative and MAX-positive GIST cases. Selected molecular, demographic, and clinical characteristics of GIST sample sets stratified by MAX immunohistochemical expression. (DOC 60 kb)

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Belinsky, M.G., Rink, L., Cai, K.Q. et al. Somatic loss of function mutations in neurofibromin 1 and MYC associated factor X genes identified by exome-wide sequencing in a wild-type GIST case. BMC Cancer 15, 887 (2015). https://doi.org/10.1186/s12885-015-1872-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-015-1872-y

Keywords