Mutations in PIK3CA are infrequent in neuroblastoma

Background Neuroblastoma is a frequently lethal pediatric cancer in which MYCN genomic amplification is highly correlated with aggressive disease. Deregulated MYC genes require co-operative lesions to foster tumourigenesis and both direct and indirect evidence support activated Ras signaling for this purpose in many cancers. Yet Ras genes and Braf, while often activated in cancer cells, are infrequent targets for activation in neuroblastoma. Recently, the Ras effector PIK3CA was shown to be activated in diverse human cancers. We therefore assessed PIK3CA for mutation in human neuroblastomas, as well as in neuroblastomas arising in transgenic mice with MYCN overexpressed in neural-crest tissues. In this murine model we additionally surveyed for Ras family and Braf mutations as these have not been previously reported. Methods Sixty-nine human neuroblastomas (42 primary tumors and 27 cell lines) were sequenced for PIK3CA activating mutations within the C2, helical and kinase domain "hot spots" where 80% of mutations cluster. Constitutional DNA was sequenced in cases with confirmed alterations to assess for germline or somatic acquisition. Additionally, Ras family members (Hras1, Kras2 and Nras) and the downstream effectors Pik3ca and Braf, were sequenced from twenty-five neuroblastomas arising in neuroblastoma-prone transgenic mice. Results We identified mutations in the PIK3CA gene in 2 of 69 human neuroblastomas (2.9%). Neither mutation (R524M and E982D) has been studied to date for effects on lipid kinase activity. Though both occurred in tumors with MYCN amplification the overall rate of PIK3CA mutations in MYCN amplified and single-copy tumors did not differ appreciably (2 of 31 versus 0 of 38, respectively). Further, no activating mutations were identified in a survey of Ras signal transduction genes (including Hras1, Kras2, Nras, Pik3ca, or Braf genes) in twenty-five neuroblastic tumors arising in the MYCN-initiated transgenic mouse model. Conclusion These data suggest that activating mutations in the Ras/Raf-MAPK/PI3K signaling cascades occur infrequently in neuroblastoma. Further, despite compelling evidence for MYC and RAS cooperation in vitro and in vivo to promote tumourigenesis, activation of RAS signal transduction does not constitute a preferred secondary pathway in neuroblastomas with MYCN deregulation in either human tumors or murine models.


Background
Neuroblastoma is a common childhood tumor arising from neural crest-derived cells of the sympathetic nervous system. It is one of the most extensively characterized solid tumors at the genomic level and provides the paradigm for the clinical utility of genomic alterations in defining tumor phenotype. Much of the clinical heterogeneity seen with this tumor (spanning spontaneous regression to relentless disease progression) can be attributed to distinct tumor-specific genetic alterations that correlate strongly with disease course. These include genomic loss within 1 p and 11 q, and gain of 17 q, among others. Still, little is known about the genes or cellular processes targeted by the majority of these genomic changes. Key regulatory genes frequently targeted for activation or inactivation in a plethora of adult human cancers (e.g., Ras, TP53, CDKN2A, RB1) are rarely, if ever, mutated in neuroblastoma. The sole exception to this is the MYCN proto-oncogene, which is markedly amplified and overexpressed in ~22% of neuroblastomas and independently correlates with advanced disease and adverse outcome [1,2].
Deregulated expression of MYC genes, including MYCN, inhibits cellular differentiation and promotes mitogenindependent proliferation [3]. However, sensitivity to apoptotic stressors is concurrently enhanced as a safeguard against neoplastic transformation, and therefore oncogenesis mediated by MYC genes requires co-operating lesions that either disable the MYC-death signaling axis [4,5] or repress it through constitutive survival signals. We are interested in defining these co-operating lesions, in addition to identifying oncogenic lesions that contribute to neuroblastoma tumourigenesis or progression independent of MYCN. Activating mutations in Ras family genes have been identified in 30% of human cancers [6]. Activated Ras (1) co-operates with MYC to transform primary cells in vitro and promote tumourigenesis in vivo [7]; (2) regulates myriad cellular processes including survival; and (3) acts to stabilize Myc proteins through post-translational modifications [8,9]. These data strongly support this pivotal pathway as a candidate for aberrant activation in human neuroblastoma, and specifically in co-operating with MYCN in tumors with MYCN amplification or deregulated expression.
Still, activating mutations have rarely been identified in HRAS1, KRAS2 or NRAS in human neuroblastoma [10][11][12]. Occasional Ras effector alterations have been identified but collectively they are infrequent. These include biallelic inactivation of NF1 and mutations in the Noonan-associated SHP2/PTPN11 gene [13][14][15][16]. Downstream of Ras are multiple signal transduction pathways responsible for the execution of Ras-mediated cellular effects [17] and cell-type and species specificity for transforma-tion has been noted [18]. The most widely validated as contributing to human tumorigenesis include the Raf-MAPK and phosphatidylinositol 3-kinase (PI3K) pathways [19,20]. Somatic activating mutations have been identified in multiple human cancers within these pathways, particularly in the BRAF and PIK3CA genes, respectively [21,22]. Activated PI3K signaling has been proposed to foster Myc/Ras co-operativity through Akt-mediated phosphorylation of FoxO transcription factors [23] to derepresses MYC target genes, and activated PI3K genes protect against Myc induced cell death sensitivity [24]. These Ras effector signaling components have not been studied in neuroblastoma. As PIK3CA has been shown to be activated through somatic mutation in >25% of colorectal, breast and ovarian carcinomas, as well as in human CNS tumors (including anaplastic oligodendrogliomas, glioblastoma multiforme and medulloblastoma), we hypothesized that PIK3CA mutations may substitute for activating RAS lesions in neuroblastoma tumourigenesis. Therefore, we performed a targeted sequencing analysis of the PIK3CA gene in primary neuroblastomas and neuroblastoma-derived cell lines.
To further explore Raf-MAPK and PI3K signaling in neuroblastic tumors we evaluated tumors arising in neuroblastoma-prone genetically engineered mice (GEM). In this model human MYCN gene overexpression is targeted to neural crest tissues using the rat tyrosine hydroxylase promoter [25]. Neuroblastic tumors arise stochastically within the sympathetic ganglia with near-complete penetrance. These tumors are genomically complex with chromosomal gains and losses at regions orthologous to regions similarly altered in human neuroblastoma [26,27]. Thus, the model faithfully recapitulates the human condition with deregulated MYCN as the initiating event biologically and genetically [25][26][27]. To determine whether alterations in Ras/Raf/PIK signaling were co-operating with MYCN in this model we sequenced Hras1, Kras2, Nras, Braf and Pik3ca to assess for activating mutations as these genes and pathways have not been studied previously using this model.

Human neuroblastoma specimens
Forty-two (42) primary neuroblastoma-derived DNA specimens were acquired from the Children's Oncology Group (COG; Arcadia, CA) for PIK3CA sequence analysis. Constitutional DNA samples were subsequently obtained from the COG Nucleic Acids bank for two patients with tumors demonstrating PIK3CA sequence alterations to determine whether these were tumor-specific (somatic) or germ line in origin. These studies were approved by the Institutional Review Board at The Children's Hospital of Philadelphia. Twenty-seven (27) additional neuroblastoma cell line DNA specimens, obtained from The Chil-dren's Hospital of Philadelphia cell line bank and the Wellcome Trust Sanger Institute and reported previously [28,29], were similarly investigated for PIK3CA mutations.

Murine neuroblastoma specimens
Twenty-five (25) tumors were harvested from genetically engineered mice overexpressing human MYCN in neuralcrest tissues under the rat tyrosine-hydroxylase promoter [TH-MYCN GEM model] [25]. In this model, tumors arise stochastically within the peripheral sympathetic nervous system (in the retroperitoneum or posterior thorax) with near-complete penetrance in mice homozygous for the MYCN transgene. Founders for establishment of the TH-MYCN GEM colony at The Children's Hospital of Philadelphia were generously provided by Dr. William Weiss (UCSF). Ten tumors in this study came from early backcrosses (received from UCSF) and fifteen additional tumors came from animals backcrossed to homozygosity in the 129X1/SvJ background. Tumors were flash frozen at necropsy and partitioned for subsequent DNA and RNA extraction using the DNeasy and RNeasy Midi Tissue kits, respectively, according to the manufacturer's instructions (Qiagen; Valencia, CA). All mouse experiments were approved by the Institutional Animal Care and Use Committee (IACUC) at The Children's Hospital of Philadelphia.

Murine Ras/Raf pathway genes
Primer sequences to amplify the entire coding region of murine Hras1 (GenBank:NM008284), Kras2 (Gen-Bank:NM021284), and Nras (GenBank:NM010937) from cDNA were kindly provided by Lewis Chodosh (University of Pennsylvania). These are: (1)  For sequencing from DNA, 80 ng were used as template input. For sequencing from cDNA, 2 μg RNA was reverse-transcribed in a final volume of 20 μl with the oligodeoxythymidylic acid primer using the Superscript II pre-amplification system according to the manufacturer's instructions (Life Technologies, Gaithersburg, MD). One μl of reverse transcribed cDNA was used as a template. The following PCR protocol was utilized for each sample and carried out in a Programmable Thermal Controller-100 (MJ Research, Inc., Waltham, MA): (1) 94°C × 12 minutes (polymerase activation), (2) thirty DNA polymerization cycles of 94°C × 45 seconds, 55°C × 45 seconds, 72°C × 45 seconds, followed by 72°C × 10 minutes and 4°C holding temperature. The individual PCR amplifications for each reaction were combined and electrophoresed through a 1% agarose gel in 1× TBE buffer for 40 minutes at 90 V. The PCR product was verified as a single specific band of predicted size then extracted and purified using the QIAQuick Gel Extraction Kit (Qiagen Inc., Valencia, CA).
All amplification products were sequenced bi-directionally in an automatic sequencer (Applied Biosystems, Fos-ter, CA) using the primary forward and reverse primers. This provided at least double coverage of all evaluated sequence in each tumor sample. Sequence analysis was carried out using Sequencher version 4.1.4 (Gene Codes Corporation, Ann Arbor, MI) in addition to a manual chromatogram review. Confirmatory resequencing from replicate PCR amplification reactions was performed for any sequence that was ambiguous or deviated from wildtype so that all abnormal sequences were confirmed in at least quadruplicate from replicate amplification reactions.

PIK3CA mutation screening in human neuroblastoma
We screened 69 human neuroblastomas (42 primary tumors and 27 neuroblastoma-derived cell lines) for PIK3CA gene mutations within exons 4, 5, 9 and 20 (Tables 1 and 2). PIK3CA is one of the most commonly mutated oncogenes in human cancers, occurring at high prevalence in tumors of the colon, breast, ovary, stomach and liver, and at lower prevalence in multiple other cancer types. Exons 9 and 20 encode portions of the helical and kinase domains, respectively, that ~80% of previously reported PIK3CA mutations cluster within [22,30]. These regions encode critical domains mediating lipid kinase activity [31]. Exons 4 and 5 encode the C2 region of the protein (calcium-dependent phospholipid binding domain) that is mutated in an additional ~5% of tumors, yet this region accounts for 50% of mutations reported in glioblastoma multiforme, a highly malignant neural tumor [22], so we included it in our neuroblastoma investigations ( Figure 1).
Overall, PIK3CA mutations were infrequently found, with non-synonymous coding sequence alterations being identified in 2 of 69 specimens (2.9%). An additional 16 tumors had sequence alterations identified in non-coding regions (including three specific alterations in introns and one in the 3'-untranslated region). Of the coding sequence mutations, one was identified within the helical domain encoded within exon 9, occurring in the MYCN amplified neuroblastoma cell line NLF. This heterozygous missense mutation (G to T) results in an arginine to methionine substitution (R524M). This mutation has not to our knowledge been reported previously and is of undetermined functional significance. No constitutional nor primary tumor DNA was available from the patient the cell line was derived from so it could not be directly confirmed that this was a somatic mutation, nor whether the mutation occurred in vivo or ex vivo during the course of tissue culture propagation. The second mutation was identified within the kinase domain (Figure 2A). This was a heterozygous G to T alteration leading to a glutamic acid to aspartic acid substitution (E982D) in primary tumor #1070 (INSS stage 3, MYCN amplified). Constitutional DNA from this individual was wild-type, confirming the kinase domain mutation to be tumor-specific (somatic). This specific mutation has not been previously described nor functionally characterized to date, although an E982K mutation has been found in anaplastic thryroid carcinoma [32].
In addition, four unique sequence alterations were identified in at least one specimen each in non-coding regions: two within intron 8 (C19339T and G19361A [nucleotide number refers to position in genomic sequence relative to the first nucleotide of exon 1]), one within intron 9 (T19623G), and one within the 3'-UTR (heterozygous T35577C). These were identified in 5 (7.2%), 1 (1.4%), 9 (13%) and 1 (1.4%) sample, respectively. The 3'-UTR alteration was identified in primary tumor #1218. Constitutional DNA from this patient was sequenced and demonstrated the same alteration, suggesting this represents a polymorphic variant ( Figure 2B). None of these alterations appears to disrupt splice donor or acceptor sites. However, none of these alterations represents a previously identified SNP based on a query of the dbSNP (build 125; http://www.ncbi.nlm.nih.gov). Finally, no mutations were identified within exons 4 or 5 within the C2 domain in any neuroblastoma tissue.
Thus, PIK3CA, like RAS genes, is an infrequent target for activating mutation in neuroblastoma. In our series of 69 human neuroblastomas, amino acid sequence alterations were identified in only two (one in a MYCN amplified cell line, and one in a MYCN amplified INSS stage 3 tumor), albeit in a survey restricted to domain "hotspots" (C2, helical, kinase) that harbor >80% of mutations identified to date [22]. Although these alterations occur within regions where mutations cluster, it is not known whether they augment lipid kinase activity and are functionally significant. Still, we did not identify any synonymous sequence alterations in the coding regions assessed so random nonselected alterations are infrequent. Of note, these findings are consistent with similar low frequencies of mutation (<5%) in other embryonal malignancies such as medulloblastoma [33] and leukemia/lymphomas [30], in contrast to the higher prevalence of mutations in many adult epithelial cancers. The prevalence of PIK3CA mutations in other pediatric cancers has not yet been reported.

Assessment of the Ras, Raf, PI3K pathway in murine neuroblastomas arising in the TH-MYCN genetically engineered mouse (GEM)
Neuroblastomas arise with near-complete penetrance in TH-MYCN GEM homozygous for the transgene. Tumors develop focally in the peripheral nervous system and demonstrate genomic complexity with high similarity to human neuroblastomas [26]. For this reason, we used this model to determine whether Ras, Raf and/or Pi3k signaling constituted a preferred secondary pathway with

MYCN deregulation. To assess for activating mutations in
Pi3kca we sequenced the helical and kinase encoding regions (exons 9 and 20) from 25 neuroblastomas harvested from heterozygous or homozygous transgene carrying mice. All 25 murine neuroblastomas harbored wildtype Pik3ca sequence (Table 3).
We next assessed genes in the Ras and Raf pathway in these same 25 murine neuroblastomas, including Hras1, Kras2, Nras and Braf, which are frequently mutated in human cancers and in MYC-initiated murine tumor models [34]. No non-synonymous sequence alterations were identified in these genes. One tumor (THMYCN-21) had three synonymous alterations identified: one each in the Hras1, Kras2 and Nras genes. One additional synonymous alterations was identified in THMYCN-17 in Kras2 (Table  3). These findings likely represent polymorphic variants and were found only in tumors from early mice back-crosses predicted to harbor strain specific polymorphisms. Tumors harvested following extensive backcrossing to homozygosity had no non-synonymous alterations identified. The low frequency of Pik3ca, Ras family, and Braf activating mutations in tumors arising in the TH-MYCN GEM supports that activated Ras-mediated signal transduction is not a pre-requisite for neuroblastoma development in concert with deregulated MYCN in this model. In this respect the GEM model is consistent with the human tumor.

Conclusion
Pik3ca, Ras family, and Braf activating mutations occur rarely in human neuroblastomas or in a murine model driven by neural-restricted MYCN overexpression. These data support that activated Ras-mediated signal transduction does not constitute a preferred secondary pathway for neuroblastic tumors in concert with deregulated MYCN.
However, it remains strictly possible that activating mutations in alternative Ras signaling constituents occur in neuroblastoma. Less studied family members such as MRAS [36] and RRAS [37] have oncogenic potential in mammalian cells and the former demonstrates restricted expression to neural tissues. These have not been investigated in neuroblastoma although both are expressed in primary neuroblastomas at the mRNA level (data not shown). In distinction to activating point mutations, pathway activation through PIK3CA gene dosage effects following low level genomic amplification have also been described [38]. However, these appear predominantly in those epithelial carcinomas (breast, ovarian, gastric carcinoma) in which activating mutations occur at higher frequency as well. In contrast, gain of 3q26 including the PIK3CA locus is infrequently reported in neuroblastoma cells, including in those cell lines studied herein [28]. Finally, additional effector arms to Ras-mediated signaling may play a role in tumourigenesis. These include Tiam-1, RalGDS, RASSF and others (reviewed in [17]). RASSF1A has been shown to be epigenetically silenced in a significant proportion of high-risk neuroblastomas [39] and elimination of this potentially pro-apoptotic pathway may be more critical than activation of proliferative or survival signaling in the establishment of embryonal cancers. Alternatively, Ras independent secondary pathways may be preferentially co-opted in neuroblastoma tumor initiation. These pathways remain elusive but certainly warrant further investigation, as defining the preferred pathways deregulated in neuroblastoma might provide suitable novel targets for directed therapies against this frequently lethal tumor.

Competing interests
The author(s) declare that they have no competing interests.