Genome-Wide Microarray Expression and Genomic Alterations by Array-CGH Analysis in Neuroblastoma Stem-Like Cells

Abstract

Neuroblastoma has a very diverse clinical behaviour: from spontaneous regression to a very aggressive malignant progression and resistance to chemotherapy. This heterogeneous clinical behaviour might be due to the existence of Cancer Stem Cells (CSC), a subpopulation within the tumor with stem-like cell properties: a significant proliferation capacity, a unique self-renewal capacity, and therefore, a higher ability to form new tumors. We enriched the CSC-like cell population content of two commercial neuroblastoma cell lines by the use of conditioned cell culture media for neurospheres, and compared genomic gains and losses and genome expression by array-CGH and microarray analysis, respectively (in CSC-like versus standard tumor cells culture). Despite the array-CGH did not show significant differences between standard and CSC-like in both analyzed cell lines, the microarray expression analysis highlighted some of the most relevant biological processes and molecular functions that might be responsible for the CSC-like phenotype. Some signalling pathways detected seem to be involved in self-renewal of normal tissues (Wnt, Notch, Hh and TGF-β) and contribute to CSC phenotype. We focused on the aberrant activation of TGF-β and Hh signalling pathways, confirming the inhibition of repressors of TGF-β pathway, as SMAD6 and SMAD7 by RT-qPCR. The analysis of the Sonic Hedgehog pathway showed overexpression of PTCH1, GLI1 and SMO. We found overexpression of CD133 and CD15 in SIMA neurospheres, confirming that this cell line was particularly enriched in stem-like cells. This work shows a cross-talk among different pathways in neuroblastoma and its importance in CSC-like cells.

About the Authors:

Raquel Ordóñez

Contributed equally to this work with: Raquel Ordóñez, Gabriel Gallo-Oller

Affiliation: Department of Biochemistry and Genetics, University of Navarra School of Sciences, Pamplona, Spain

Gabriel Gallo-Oller

Contributed equally to this work with: Raquel Ordóñez, Gabriel Gallo-Oller

Affiliation: Department of Biochemistry and Genetics, University of Navarra School of Sciences, Pamplona, Spain

Soledad Martínez-Soto

Affiliation: Department of Biochemistry and Genetics, University of Navarra School of Sciences, Pamplona, Spain

Sheila Legarra

Affiliation: Department of Biochemistry and Genetics, University of Navarra School of Sciences, Pamplona, Spain

Noémie Pata-Merci

Affiliation: Institut Gustave Roussy, Villejuif, France

Justine Guegan

Affiliation: Institut Gustave Roussy, Villejuif, France

Giselle Danglot

Affiliation: Institut Gustave Roussy, Villejuif, France

Alain Bernheim

Affiliation: Institut Gustave Roussy, Villejuif, France

Bárbara Meléndez

Affiliation: Molecular Pathology Research Unit, Department of Pathology, Virgen de la Salud Hospital, Toledo, Spain

Juan A. Rey

Affiliation: IdiPaz Research Unit, La Paz University Hospital, Madrid, Spain

Javier S. Castresana

* E-mail: [email protected]

Affiliation: Department of Biochemistry and Genetics, University of Navarra School of Sciences, Pamplona, Spain

Introduction

Neuroblastoma is a neuroendocrine tumor of unknown etiology, derived from primordial neural crest cells which afterwards develop into adrenal medulla and sympathetic ganglia [1]. It is one the most prevalent cancers in childhood and nearly 50% of the cases take place in children younger than two years. The estimated 5-years survival rate is 90% and 50% for patients with non-high-risk and high-risk neuroblastoma, respectively. Most neuroblastomas occur sporadically and develop with very diverse clinical behavior, from spontaneous regression to aggressive malignant progression and resistance to chemotherapy [2], [3].

The actual treatment of neuroblastoma depends on the clinic stage of the tumor, but it is commonly based on radiopharmaceutical therapy in combination with surgery [4]. Current studies and clinical trials are combining conventional chemotherapy with monoclonal antibodies, stem cell transplants and retinoids, but because of the complexity of this pathology, progress remains extremely slow [5].

Some authors propose that this diverse clinical behavior of neuroblastoma might be due to molecular differences in cell subpopulations [6]. The cancer stem cells (CSC) model might be an explanation for this heterogeneous behavior [7]. This model proposes that only a small subpopulation with characteristics of stem cells within the tumor has the ability to proliferate and maintain its growth. Even if a tumoral mass shows a substantial decrease in size in response to therapy, if the CSC are spared, it will regrow leading to a relapse [8], [9]. Some studies propose that CSC operate with the machinery and developmental programs expressed in normal stem cells [10], [11]. There is growing evidence of some signalling pathways involved in self-renewal of both normal and tumor tissues as Wnt [12], [13], Sonic Hedgehog (Hh) [14], [15], Notch [16], [17] and Transforming Growth Factor Beta (TGF-β) [18]–[20] signalling pathways, that might contribute to tumorigenesis when deregulated [14], [21].

In this line, the development of new therapies based in molecular targets may be of great value for the treatment of neuroblastoma. The observation of some pathways acting on multiple levels to promote the development of neuroblastoma and CSC subpopulation, has prompted new therapeutic strategies to treat not only this neoplasm but other brain and nervous system tumors [2], [18], [22].

Therefore, this project performs a genomic analysis of CSC-like by array-CGH and expression array, with special focus on altered signalling pathways that might explain the stem cell phenotype of the CSC subpopulation. The aim of this study is to provide a powerful tool to open up new targets for therapy or redirect current cancer treatments towards CSC in order to achieve total elimination of tumor cell population and improve treatment effectiveness.

Materials and Methods

Cell lines culture

Two commercial neuroblastoma-derived cell lines were used: SK-N-DZ cell line (ATCC N°CRL-2149) provided by the American Type Culture Collection (ATCC, Manassas, VA, USA) and SIMA cell line (DSMZ N° ACC 164) provided by the German Collection of Microorganisms and Cell Cultures (DSMZ, Braunschweig, Germany). Both were grown as monolayer with DMEM/GlutaMAX medium supplemented with 10% fetal bovine serum, 5% non-essential aminoacids, 1% penicillin/streptomycin and 0.1% amphotericin B. To obtain the neurospheres cultures in order to enrich them in CSC-like cells, after chemical dissociation of SK-N-DZ and SIMA cell lines, 5×10⁵ cells were transferred to 25 cm² flasks (positioned vertically) and grown in suspension with 5 ml of selective medium: DMEM-F12/GlutaMAX plus 1% penicillin/streptomycin and 0.1% amphotericin B, supplemented with Epidermal Growth Factor (EGF) (20 ng/ml), Fibroblast Growth Factor Basic (FGFb) (20 ng/ml) and B27 Supplement (1X). All cultures were maintained at 37°C in a humidified atmosphere of 5% CO₂/95% air.

DNA/RNA isolation and reverse transcription

DNA from 1×10⁶ cells pellet was extracted using Wizard Genomic DNA Purification Kit (Promega, Madison, WI, USA) and total RNA from 2×10⁶ cells pellet, employing Illustra RNAspin Mini Kit (GE Healthcare, Buckinghamshire, UK) following manufacturer's instructions in both cases. DNA and RNA samples were stored at −20°C and −80°C until utilization, respectively. DNA and RNA concentrations were measured with a SmartSpec Spectophotometer (Bio-Rad, Hercules, CA, USA). Reverse transcription was carried out from 2 µg of total RNA (previously denatured at 65°C for 5 min) with a mixture of 0.5 µg of random primers, 0.5 mM dNTPs, 10 mM DTT, First-Strand buffer 1X and 200U SuperScript II Reverse Transcriptase (Invitrogen, Carlsbad, CA, USA) in a total volume of 40 µl. The mixture was incubated for 10 min at 25°C, 50 min at 42°C and 15 min at 70°C, followed by chilling on ice. The cDNA was diluted 1/5 and stored at −80°C until utilization.

Array-CGH

DNA samples were analyzed using a 244K microarray (Agilent Technologies, Santa Clara, CA, USA). Oligonucleotide aCGH processing was performed as detailed in the manufacturer's protocol (version 4.0; http://www.agilent.com). Data were extracted from scanned images using feature extraction software (version A.8.5.3, Agilent). Raw data text files from the latter were then imported for analysis into CGH Analytics 3.4.40. Aberrations were detected with the ADM2 algorithm and filtering options of a minimum of 5 probes and abs(log2Ratio)>0.3. Aberration segments were individually reviewed using build 35, hg17 of UCSC. Anomalies that were localized to regions with high-copy repetitive or GC-rich DNA sequences including telomeric regions were excluded. Gains and losses for the oligonucleotide dataset were defined as a linear ratio ≥1.2 or ≤0.8, respectively. High and low-level amplification events were defined as a linear ratio ≥4 or 2<ratio<4, respectively. The data are described in accordance with MIAME guidelines and have been deposited in ArrayExpress under E-MTAB-2866 accession number. All the analysis and statistics related with the array-CGH was performed at Institut Gustave Roussy by the bioinformatics team.

Microarray experiment, data normalization and analysis

SurePrint G3 human Gene Expression Microarray Kit (Agilent, ID:028004) was used to analyze the transcriptional profiles of all samples. This array contains 42,404 60-mer oligonucleotides targeting 27,988 Entrez Gene RNAs and 7,419 lincRNAs (long intergenic non coding RNAs) in an 8×60K format slide (each slide contains 8 arrays of 60,000 features with 1440 spots reserved for internal quality control). The transcripts used can be checked in the EMBL-EBI Database (http://www.ebi.ac.uk/arrayexpress/arrays/A-GEOD-14550).

Array performance and analysis were carried out at the Institut Gustave Roussy, using the Agilent Feature Extraction software. Data sets were normalized using quantile normalization, and fold change between cell line and neurospheres was calculated for each gene. The whole dataset was filtered by intensity of the probes. Spots with too high intensities (saturating intensities) and irregular spots (spots inside whose neighbouring pixel showed very different intensities) were not further considered. Likewise, genes with a signal-to-noise ratio lower than 6 fold the average standard deviation observed for all negative control probes were excluded. Finally, >1 fold-change between normalized gene expression measured in neurospheres (or CSC-like cells) versus standard cell lines was used as the criteria to select differentially expressed genes.

The list of differentially expressed genes in both cell lines was imported and classified by DAVID (Database for Annotation, Visualization and Integrated Discovery, http://david.abcc.ncifcrf.gov/), version 6.7 [23]–[25]. Enriched categories in biological processes and molecular function were defined by a p value <0.05 and more than 2 genes differentially expressed in each category. After the enrichment analysis, a classification by PANTHER (Protein ANalysis THrough Evolutionary Relationships, www.pantherdb.org/) Data Base [26]–[28] was performed. Finally, a deeper pathway analysis was carried out by comparing the same set of genes with the NCBI BioSystems Database (http://www.ncbi.nlm.nih.gov/biosystems). Microarray expression data are available at ArrayExpress with accession number E-MTAB-2867.

Real-Time polymerase chain reaction (RT-qPCR)

A first screening of our primer library was performed according to the microarray data and pathway analysis. 44 pairs of primers were selected and aligned with the mRNA and DNA sequences obtained from RefSeq (http://www.ncbi.nlm.nih.gov/RefSeq/) and UCSC Genome Browser (http://genome.ucsc.edu/). A deeper analysis of each pair of primers was performed with Primer3 (http://primer3.sourceforge.net/) and Primer Blast (http://www.ncbi.nlm.nih.gov/tools/primer-blast/) in order to obtain product length, %GC, melting temperature (Tm), self-complementarity or possible products on unintended templates. After the computational analysis, 17 out of 44 primers (Table 1) were selected for further studies. RT-qPCRs were performed in order to define the optimal Tm for each pair of primers. The melting curves obtained were analyzed to check amplification efficiency in each temperature, and detect additional peaks displaced from the desired amplicon peak (determined by agarose gels checking size and identifying possible intron amplifications). For the quantitative PCR reaction, 2.5 µl of the generated cDNA was added to a PCR mix containing 12.5 µl SYBR Green Supermix (Bio-Rad, Hercules, CA, USA) and 12.5 pM forward and reverse primers, in a total volume of 25 µl. The reaction was performed according to the following protocol: 10 min at 95°C, followed by 40 cycles of 15 s at 95°C, 30 s at the specific Tm for each primer (Table 1) and 20 s at 72°C. A melting curve was added at the end of the protocol as a quality control. All samples were run in triplicate and three independent experiments were carried out. Gene expression levels between samples were normalized using the expression levels of the HPRT1 gene. Relative gene expression was analyzed according to the 2^−ΔΔCt method [29].

[Figure omitted. See PDF.]

Table 1. Primer sequences.

https://doi.org/10.1371/journal.pone.0113105.t001

Statistical analysis

The statistical analysis of arrays are detailed in the description for each technique. The data of gene expression measured by RT-qPCR was graph as the mean ± standard deviation. The 2^−ΔΔCt values were analysed for normality distribution and homogeneity of variance. Since the data did not show normality distribution, the statistical significance between groups was compared by pairs (the expression level of each gene in standard cultures vs neurospheres cultures) with the non-parametric Mann-Whitney U test. The significance level was set as p<0.05. The statistical program GraphPad Prism for Windows version 5.04 was used.

Results

Genomic profile by array-CGH analysis

To identify the potential common chromosomal alterations indicating possible zones implicated in the CSC-like phenotype, an array-CGH analysis was performed. In the case of SK-N-DZ, we only found six chromosomal regions amplified and no deleted region. SIMA cell line showed more regions of gain and loss, specifically 10 and 30 regions, respectively (data summarized in Table 2). Notwithstanding we did not find any region that was lost or gained in both cell lines and the significance was found only in some points of the genome, we identified numerous chromosomal sequences with copy number variants (CNV) both in SK-N-DZ and SIMA cell lines (Table 2). Whole genome profiles are presented in Figure 1. Among the genes that were included in these regions identified by the array-CGH, the analysis by PANTHER did not show any pattern or group of genes directly related with signalling pathways involved in the CSC-like phenotype. Only a few genes were associated with pathways already described in tumorigenesis and stemness (as p53) but were not principal effectors and its position was downstream, including genes that were related with different regulation points and pathways (data not shown) making the analysis difficult. Despite these results, we noticed that several amplified areas included CpG islands (Table 2). Using the miRBase of the University of Manchester (www.mirbase.org), we were able to identify the microRNAs included in the gained areas of the chromosomes 1, 2 and 12 (Table 3).

[Figure omitted. See PDF.]

Figure 1. Whole chromosome plots.

Array-CGH from SK-N-DZ (A) and SIMA (B) cell lines. The X-axis represents the chromosomes while the Y-axis represent the normalize log2 ratio fluorescence intensity thresholds −1 (loss) and 1 (gain). The results show gains and losses of small chromosomal regions.

https://doi.org/10.1371/journal.pone.0113105.g001

[Figure omitted. See PDF.]

Table 2. Overview on gained (G) and lost (L) chromosomal regions detected by array-CGH.

https://doi.org/10.1371/journal.pone.0113105.t002

[Figure omitted. See PDF.]

Table 3. MicroRNA identified in array-CGH analysis.

https://doi.org/10.1371/journal.pone.0113105.t003

Gene expression profiles by microarray analysis

In order to characterize the differentially expressed genes in CSC-like compared to the standard tumor cell line, mRNA expression was analyzed by microarray analysis. After the first intensity filter, 25,368 genes out of the whole dataset of 42,405 (∼60%), reached the minimum quality threshold and were selected for the differential expression analysis. When setting a threshold value of >1 fold-change, 4,831 genes (∼11.4%) were identified as differentially expressed in SK-N-DZ neurospheres and 6,613 genes (∼15.6%) in SIMA neurospheres. Figure 2A shows the number of genes differentially up and downregulated in each cell line. Among these gene sets, only 757 genes have significantly changed expression in both cell lines (Figure 2B). The sets of genes up and dowregulated for each cell line are summarized in Table S1 and the genes shared by both cell lines are represented in Tables S2, S3 and S4.

[Figure omitted. See PDF.]

Figure 2. Upregulated (UR, red) and downregulated (DR, green) genes (>1 fold-change) in neurospheres compared to standard cells lines SK-N-DZ and SIMA (A), represented with MAplots, were the Y-axis represents the red/green intensity ratio “M” and the X-axis the average intensity “A”.

Venn diagrams illustrating the number of differentially expressed genes in SK-N-DZ and SIMA neurospheres compared to standard cell lines (B). The overlapping area represents the set of genes altered in both cell lines.

https://doi.org/10.1371/journal.pone.0113105.g002

Differentially expressed genes with >1 fold-change in microarray analysis were categorized according to the biological processes they are involved in, and to the molecular functions they code for, by using the DAVID database. Among all categories, those that appeared significantly altered or enriched (p<0.05 and more than 2 genes differentially expressed in both cell lines by category) were selected for further studies (Table 4).

[Figure omitted. See PDF.]

Table 4. DAVID generated categories.

https://doi.org/10.1371/journal.pone.0113105.t004

Further analysis of the significantly enriched categories and signalling pathways involved was performed with PANTHER database and information from the NCBI BioSystems database. As a result, we observed significant alteration in representative signalling pathways known to be also involved in stem cell self-renewal, as Wnt, Notch, Hh and TGF-β (Figure 3).

[Figure omitted. See PDF.]

Figure 3. PANTHER classification by signalling pathway.

The differentially expressed genes in both SK-N-DZ and SIMA CSC-like cells were classified by PANTHER and graphed. The percentage represents the number of genes altered against total number of genes involved in each pathway.

https://doi.org/10.1371/journal.pone.0113105.g003

Biological processes and signalling pathways expression in CSC-like subpopulation

The expression microarray experiment showed up a high percentage of genes altered in neurospheres (Figures 2 and 3). A complete analysis of the data highlighted some of the most relevant biological processes and molecular functions affected (Table 4) that might be responsible for the CSC-like phenotype of this subpopulation. First, we observed an apparent alteration of ion transporters. In our analysis, several genes of Ca²⁺ and K⁺ channels showed up as differentially expressed, as KCNMA1 and CACN1AG. Secondly, we detected alterations in some of the most commonly studied signalling pathways in brain tumors and stem cell development: Wnt, Notch, Hh and TGF-β.

Analyzing Wnt signalling pathway, we observed downregulation of some of the upstream regulators of the pathway: FZD8 (Frizzled receptor), PRKCH (protein kinase C) and FRAT1, inhibitor of GSK3-mediated phosphorylation of β-catenin, and therefore a positive regulator of the pathway that could act as a proto oncogene [30], [31]. The upregulation of POU5F1 as a downstream target of the pathway was also noticed.

In Notch signalling pathway, overexpression of the activator ligand of the pathway, JAG1 was identified. We also detected upregulation of HES7, one of the final targets of the pathway [32], [33]. This evidence showed an apparent activation of the pathway in CSC-like comparing with the standard cell line.

Expression microarray data showed upregulation of PTHC1, the main inhibitor of Hh signalling pathway, but also one of the gene targets of the pathway which transcription would be increased if the pathway is active, as we confirmed with the RT-qPCR data (Figure 4A).

[Figure omitted. See PDF.]

Figure 4. RT-qPCR for Sonic Hedgehog (A) and TGF-β pathways analysis (B) in SIMA CSC-like cells.

The graphs represent the 2^−ΔΔCt values obtained by RT-qPCR for neurospheres. The dotted line indicates the 2^−ΔΔCt control (cell line) value equal to 1. Significance against control: p<0.05 (*); p<0.01 (**); p<0.001 (***) and p<0.0001 (****). Selected genes cover the most relevant components of each pathway. Results confirmed aberrant activation of Hh (A) and TGF-β (B) signalling pathways in CSC-like cells.

https://doi.org/10.1371/journal.pone.0113105.g004

Likewise, when analyzing TGF-β signalling pathway, the expression microarray data suggested an activation of the pathway (Figure 4B). Ligands like TGFB3 and BMP8B were upregulated, as well as the signal transducer SHC. Quite the contrary, we found downregulation of inhibitory SMADs (like SMAD6 and SMAD7) and the pseudoreceptor BAMBI, which could sequester ligands and inhibit signal transduction by different mechanisms [34].

Parallel behaviour in SK-N-DZ and SIMA stem-like cells

RT-qPCR was performed in order to corroborate the expression microarray results. In both cell lines, we confirmed the decrease in the expression of SMAD6, SMAD7 and BAMBI, as the increase observed in JAG1. The results suggested that the CSC-like subpopulation of both cell lines could have a similar expression profile after enrichment with the neurosphere assay (Figure 5A).

[Figure omitted. See PDF.]

Figure 5. Gene expression by RT-qPCR.

The graphs represent the 2^−ΔΔCt values obtained by RT-qPCR for neurospheres. The dotted line indicates the 2^−ΔΔCt control (cell line) value equal to 1. Significance against control: p<0.05 (*); p<0.01 (**); p<0.001 (***) and p<0.0001 (****). Corroboration of expression array (A). Both SK-N-DZ (▪) and SIMA (▴) CSC-like cells showed a similar expression profile. Stem Cell Markers expression (B). Results confirmed a high overexpression of CSC markers in SIMA cell line after the neurosphere formation assay.

https://doi.org/10.1371/journal.pone.0113105.g005

Hh and TGF-β analysis in SIMA stem-like cells

Based on our previous experience in both Hh and TGF-β signalling pathways, further expression analysis by RT-qPCR was performed on some of the most relevant components of both signalling networks in SIMA cell line (Figure 4). Results confirmed Hh aberrant activation by overexpression of its key components as PTCH1, GLI1 and SMO. Likewise, TGF-β seemed to be activated by overexpression of its receptors (TGFBR1 and TGFBR2) and ligands (TGFB2), (despite TGFB3 is downregulated even though it is also a ligand of the pathway), and downregulation of its repressors (SMAD6, SMAD7 and BAMBI). We confirmed by RT-qPCR the increase of the levels of mRNA of JAG1, a Notch activating ligand that is related with TGF-β pathway [35]. Its upregulation supports Notch activation as aforementioned.

Expression of stem-cell markers

To confirm the enrichment of the isolated population in CSC-like cells, analysis of CSC-specific cell surface markers was performed. RT-qPCR on CD133 and CD15 was carried out in both SK-N-DZ and SIMA neurospheres (Figure 5B). On the basis of these results, SIMA culture was highly enriched in CSC-like cells as shown by the significant increase in CSC surface markers expression. An interesting point is that SK-N-DZ subpopulation did not show any overexpression of the stem-cell markers analyzed.

Discussion

The heterogeneous clinical behavior of neuroblastoma tumors has been widely studied [36]–[39]. It is already described that genomic rearrangements cannot explain per se this heterogeneity neither the CSC-like subpopulation [40]. In this work, we hypothesized that this heterogeneity might be due to genomic rearrangements together with changes in the expression profile of tumor cell subpopulations, particularly in CSC-like cells. In connection with this model, the genomic profile of CSC-like cells and the regulatory mechanisms or signalling pathways involved in CSC promotion could explain the heterogeneity and biology of this subpopulation and help to its identification in neuroblastoma.

Chromosomal aberrations are common in neuroblastoma, including numerical whole chromosomal gains, segmental chromosomal gains and losses and somatic mutations [2]. However, array-CGH did not show large modified areas in the genome of CSC-like subpopulations neither individually nor shared by both cell lines (Table 2). Despite this result, we identified other genetic elements involved in the biology of neuroblastoma and in the regulation of several tumors as CNV [41]–[43] and CpG islands [44]–[46] (Table 2). Both elements are being investigated as possible therapeutic targets in neuroblastoma [5] and were already identified in other tumors confirming their therapeutic potential [47]–[51]. These results suggest the fact that these two elements could be involved in tumorigenesis and CSC-like cells generation in neuroblastoma, and, as in other tumors, this might be a starting point to develop or improve current therapies. This is the case of decitabine (5-aza-2′-deoxycytidine). This agent interferes with DNA methylation but the preclinical assays show that it must be administrated at doses that are not tolerable in order to produce a biological effect [52]. Our findings propose a genetic basis for the effectivity of this kind of compounds in neuroblastoma highlighting the importance for further research in demethylating agents.

In addition, some microRNAs were identified in lost and gained areas (Table 3). In the last years, these elements have become promising therapeutic targets in neuroblastoma [53]. For example, in amplified loci of SIMA neurospheres, hsa-mir-4254 and hsa-mir-4262 were indentified (Table 3). These two microRNAs were already related with the stem cell phenotype [54], supporting the hypothesis that the CSC might be using the same regulatory machinery as normal stem cells [10], [11]. Other microRNA identified were hsa-mir-4420, hsa-mir-4429 and hsa-mir-4497, all described as possible new microRNAs involved in malignant human B cells regulation [55], and hsa-mir-3605, a novel microRNA linked to human cervical cancer [56].

In contrast, the expression array results indicated a modified expression pattern in CSC-like cells including changes in different processes and cell functions (Figures 2, 3 and Table 4). Different studies indicated some of the most important pathways involved in the development and maintenance of CSC subpopulation. Among them, the most remarkable are Wnt, Notch, Hh and TGF-β, which have been already described for different tumors and especially in brain tumors [14]. In neuroblastoma, Wnt and Notch genes have been identified as possible CSC-like markers. The difficulty in the identification of specific genes and accurate markers of CSC-like lies in the heterogeneity of this tumor and the poor knowledge related with this subpopulation, since no specific and definitive gene expression profile has been proposed for neuroblastoma CSC-like cells [57].

Interestingly, the expression pattern in CSC-like cells involves not only Wnt and Notch but also Hh and TGF-β signalling pathways. For example, deregulation of several genes of Wnt pathway was detected. The results indicated a decrease in the expression of upstream regulators as FZD8, PRKCH and FRAT1 and an increase in the downstream target POU5F1, suggesting the inactivation of Wnt signalling pathway. POU5F1 was described in several tumors (urothelial, prostatic, cervical, breast and lung) and is implicated in different functions as chemoresistance, prognosis marker, proliferation, apoptosis, migration and invasion [58]–[61]. In addition to all these functions, this gene plays an important role in stemness and self-renewal of stem cells in normal and tumor tissues [60], [62], [63] and was found to be overexpressed and correlated with the progression of neuroblastoma [64]. This gene could interact with members from other pathways as TGF-β, Notch and Wnt [65], [66] working together in the regulation of stem cell pluripotency. Our results show alterations not only in POU5F1 but also in members of TGF-β, Notch and Wnt signalling pathways indicating the potential of this gene to be used as a new therapeutical target. Our findings are supported by other works that described its potential in other tumors [62], [64], [67] and as a possible new CSC marker in neuroblastoma [67]–[69].

In this context, an apparent activation of Notch pathway was detected in neurospheres. It has been demonstrated that inhibition of Notch in cancer cells has the potential to slow down cell proliferation and induce apoptosis, despite the extensive crosstalk of this pathway with other major cancer pathways as Ras, Akt, NF-κB, Wnt, Hh and TGF-β [70]–[72]. This might happen because in neuroblastoma as in other tumors, the CSC subpopulation seems to be especially sensitive to inhibition of stem cell pathways as Notch [73].

Additionally, Hh was activated in CSC-like cells by upregulation of transcription targets and activators (PTCH1, SMO and GLI1) regardless downregulation of GLI2. It is known that GLI2 is a bi-functional transcription factor, with activator and repressor regions. Therefore the pathway could still be activated if GLI2 is silenced depending on the context [74], [75]. In neuroblastoma CSC-like the Gli factors act through cooperative functional interactions in target gene regulation [76] although in our case it seems that GLI2 does not play an important role. Alike, the TGF-β signalling pathway is activated in neurospheres cultures (Figure 4B) despite the downregulation of TGFB3, a ligand of the pathway. A possible explanation could be that the downregulation of BAMBI (Figure 4B) would let a greater proportion of ligands free, increasing the concentration of extracellular ligands and consequently facilitating the activation of the pathway or acting as a negative regulator of the TGF-β pathway [34]. Interestingly, JAG1, one of the targets of TGF-β [77] and also the main activator of Notch signalling pathway is upregulated, supporting Notch signalling pathway activation. Kurpinski et al described JAG1 as a crosstalk point between these two pathways in the regulation of muscular stem cells [78] which probably suggests an important role of this gene in the generation of CSC-like cells in neuroblastoma. These findings also give us the possibility of studying targets as JAG1 that participate in several genetic pathways involved in CSC-like cells phenotype.

It is clear that the expression profile that we observed in the CSC-like cells is similar to other profiles described for this subpopulation in different tumors including neuroblastoma. The identification of CSC is based on the expression of different markers, but in the case of neuroblastoma it remains unclear which markers are accurate to identify this subpopulation, due to its heterogeneity. Our findings show differences between the two cell lines used (Figure 5B) in two CSC markers widely described: CD133 and CD15. This result could highlight the limitations of neurospheres culture [79] since it was proved that differences in the culture conditions could affect the enrichment and isolation of CSC-like cells [80]. However, since the genes related with CSC followed a similar pattern in both cell lines, the more plausible explanation is that the phenotype of CSC-like cells is not yet fully characterized in neuroblastoma [57]. Coulon et al [57] questioned the expression of CD133 to identify CSC in neuroblastoma, highlighting the need of characterization of other CSC markers or expression profiles that allow the recognition of this subpopulation. The choice of CSC markers is still controversial and the search of specific genes that help to identify it is the real challenge for the characterization and research of new therapeutical targets in neuroblastoma [81].

Regardless the growing evidence in new approaches and in the identification of new therapeutic targets, the treatment for neuroblastoma is still evasive. Nowadays, it includes surgery, radiation and/or chemotherapy depending on the patient's stage and risk stratification of the disease [5]. Only some new compounds are included in the schema as 13-cis-retinoic acid together with anti-GD2 antibodies and interleukin-2, a combination that has increased the progression-free survival [82]. The literature includes an increasing number of new compounds with great potential to be introduced in the treatment. This is the case of drugs that targeted MYCN and ALK inhibitors that currently are in preclinical evaluation [83]. The evidence of possible targets involved in the development and regulation of other tumors could be an approach to develop new strategies in neuroblastoma. A good example are the different pathways involved in tumorigenesis and CSC regulation, as Wnt [62], [64], [67]–[69], Notch [84], Hh [85]–[87] and TGF-β [18], [88], [89]. In this connection, an interesting case are the ion transporters, described as essential for cell proliferation and which appear to have a role in cancer development [90]–[92]. Evidence is particularly extensive for Ca²⁺ and K⁺ channels [93] and Romania et al already described their importance in neuroblastoma [94]. In our analysis, several genes coding for Ca²⁺ and K⁺ channels were altered in CSC-like subpopulation like KCNMA1, previously described in prostate cancer [95] and in breast cancer [96], [97], and CACN1AG also altered in gastric cancers, colorectal cancers and acute myeloid leukemia (AML) [98]. Several reports propose that ion transporters have clinical potential not only as therapeutical targets but also as prognosis markers or to improve actual treatments [99]–[103], as blocking channel activity seems to impair the growth of some tumors, including neuroblastoma [99], [104]. The presence of altered genes of Ca²⁺ and K⁺ channels in neuroblastoma CSC-like subpopulations confirms and opens the possibility to develop new compounds that could act at different levels and subpopulations within neuroblastoma. This is only one example of the importance of redirecting treatments towards more effective CSC molecular therapies to achieve total elimination of the tumor cell population and avoid relapses [40], [105].

In conclusion, we found that the expression profile of neuroblastoma CSC-like cells differs significantly from the expression profile in standard neuroblastoma cell lines, which is seen as altered signalling pathways involved in stem cell proliferation pathways as Wnt, Notch, Hh and TGF-β, suggesting a cross talk among them and with other pathways. The CSC markers analysis reveals no correlation between the cell lines, displaying the heterogeneity of neuroblastoma tumors and CSC-like population but opening the possibility to identify other potential markers. Taking together, these results confirm the importance of different pathways in CSC regulation and the identification of possible candidate targets for molecular CSC therapies in neuroblastoma. The results presented highlight new information for the CSC-like phenotype in neuroblastoma and indicate the importance of redirecting current cancer treatments towards CSC molecular therapies to achieve total elimination of the tumor cell population and improve treatment effectiveness.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

Supplementary data. Up and downregulated genes list from the expression array in each cell line.

https://doi.org/10.1371/journal.pone.0113105.s001

(XLS)

Table S2.

Supplementary data. Upregulated genes list shared in both cell lines.

https://doi.org/10.1371/journal.pone.0113105.s002

(XLS)

Table S3.

Supplementary data. Downregulated genes list shared in both cell lines.

https://doi.org/10.1371/journal.pone.0113105.s003

(XLS)

Table S4.

Supplementary data. Genes list not found in PANTHER classification.

https://doi.org/10.1371/journal.pone.0113105.s004

(XLS)

Author Contributions

Conceived and designed the experiments: GGO AB JSC. Performed the experiments: RO GGO SM SL NP JG GD. Analyzed the data: GGO BM JAR JSC. Contributed reagents/materials/analysis tools: AB BM JAR JSC. Wrote the paper: RO GGO JSC.

Citation: Ordóñez R, Gallo-Oller G, Martínez-Soto S, Legarra S, Pata-Merci N, Guegan J, et al. (2014) Genome-Wide Microarray Expression and Genomic Alterations by Array-CGH Analysis in Neuroblastoma Stem-Like Cells. PLoS ONE 9(11): e113105. https://doi.org/10.1371/journal.pone.0113105

References

1. Dhir S, Wheeler K (2010) Neonatal neuroblastoma. Early Hum Dev 86: 601–605.

2. Cheung NK, Dyer MA (2013) Neuroblastoma: developmental biology, cancer genomics and immunotherapy. Nat Rev Cancer 13: 397–411.

3. Park JR, Bagatell R, London WB, Maris JM, Cohn SL, et al. (2013) Children's Oncology Group's 2013 blueprint for research: neuroblastoma. Pediatr Blood Cancer 60: 985–993.

4. Maris JM (2010) Recent advances in neuroblastoma. N Engl J Med 362: 2202–2211.

5. Davidoff AM (2012) Neuroblastoma. Semin Pediatr Surg 21: 2–14.

6. Matthay KK, George RE, Yu AL (2012) Promising therapeutic targets in neuroblastoma. Clin Cancer Res 18: 2740–2753.

7. Ross RA, Spengler BA (2007) Human neuroblastoma stem cells. Semin Cancer Biol 17: 241–247.

8. Singh SK, Clarke ID, Hide T, Dirks PB (2004) Cancer stem cells in nervous system tumors. Oncogene 23: 7267–7273.

9. Reya T, Morrison SJ, Clarke MF, Weissman IL (2001) Stem cells, cancer, and cancer stem cells. Nature 414: 105–111.

10. Alison MR, Lin WR, Lim SM, Nicholson LJ (2012) Cancer stem cells: in the line of fire. Cancer Treat Rev 38: 589–598.

11. Sumer-Turanligil NC, Cetin EO, Uyanikgil Y (2013) A contemporary review of molecular candidates for the development and treatment of childhood medulloblastoma. Childs Nerv Syst 29: 381–388.

12. Gong A, Huang S (2012) FoxM1 and Wnt/beta-catenin signaling in glioma stem cells. Cancer Res 72: 5658–5662.

13. Holland JD, Klaus A, Garratt AN, Birchmeier W (2013) Wnt signaling in stem and cancer stem cells. Curr Opin Cell Biol 25: 254–264.

14. Faigle R, Song H (2013) Signaling mechanisms regulating adult neural stem cells and neurogenesis. Biochim Biophys Acta 1830: 2435–2448.

15. Huang FT, Zhuan-Sun YX, Zhuang YY, Wei SL, Tang J, et al. (2012) Inhibition of hedgehog signaling depresses self-renewal of pancreatic cancer stem cells and reverses chemoresistance. Int J Oncol 41: 1707–1714.

16. Sethi N, Kang Y (2011) Notch signalling in cancer progression and bone metastasis. Br J Cancer 105: 1805–1810.

17. Shao H, Huang Q, Liu ZJ, Keiran SMS (2012) Chapter Seven - Targeting Notch Signaling for Cancer Therapeutic Intervention. Advances in Pharmacology: Academic Press. pp. 191–234.

18. Joseph JV, Balasubramaniyan V, Walenkamp A, Kruyt FA (2013) TGF-beta as a therapeutic target in high grade gliomas - promises and challenges. Biochem Pharmacol 85: 478–485.

19. Park KS (2011) TGF-beta Family Signaling in Embryonic Stem Cells. Int J Stem Cells 4: 18–23.

20. Sakaki-Yumoto M, Katsuno Y, Derynck R (2013) TGF-beta family signaling in stem cells. Biochim Biophys Acta 1830: 2280–2296.

21. Zhou BB, Zhang H, Damelin M, Geles KG, Grindley JC, et al. (2009) Tumour-initiating cells: challenges and opportunities for anticancer drug discovery. Nat Rev Drug Discov 8: 806–823.

22. Wainwright DA, Nigam P, Thaci B, Dey M, Lesniak MS (2012) Recent developments on immunotherapy for brain cancer. Expert Opin Emerg Drugs 17: 181–202.

23. Huang da W, Sherman BT, Stephens R, Baseler MW, Lane HC, et al. (2008) DAVID gene ID conversion tool. Bioinformation 2: 428–430.

24. Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, et al. (2003) DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol 4: P3.

25. Huang da W, Sherman BT, Tan Q, Collins JR, Alvord WG, et al. (2007) The DAVID Gene Functional Classification Tool: a novel biological module-centric algorithm to functionally analyze large gene lists. Genome Biol 8: R183.

26. Mi H, Dong Q, Muruganujan A, Gaudet P, Lewis S, et al. (2010) PANTHER version 7: improved phylogenetic trees, orthologs and collaboration with the Gene Ontology Consortium. Nucleic Acids Res 38: D204–210.

27. Mi H, Muruganujan A, Thomas PD (2013) PANTHER in 2013: modeling the evolution of gene function, and other gene attributes, in the context of phylogenetic trees. Nucleic Acids Res 41: D377–386.

28. Thomas PD, Campbell MJ, Kejariwal A, Mi H, Karlak B, et al. (2003) PANTHER: a library of protein families and subfamilies indexed by function. Genome Res 13: 2129–2141.

29. Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods 25: 402–408.

30. Thomas GM, Frame S, Goedert M, Nathke I, Polakis P, et al. (1999) A GSK3-binding peptide from FRAT1 selectively inhibits the GSK3-catalysed phosphorylation of axin and beta-catenin. FEBS Lett 458: 247–251.

31. van Amerongen R, Nawijn MC, Lambooij JP, Proost N, Jonkers J, et al. (2010) Frat oncoproteins act at the crossroad of canonical and noncanonical Wnt-signaling pathways. Oncogene 29: 93–104.

32. Bessho Y, Miyoshi G, Sakata R, Kageyama R (2001) Hes7: a bHLH-type repressor gene regulated by Notch and expressed in the presomitic mesoderm. Genes Cells 6: 175–185.

33. Kochert K, Ullrich K, Kreher S, Aster JC, Kitagawa M, et al. (2011) High-level expression of Mastermind-like 2 contributes to aberrant activation of the NOTCH signaling pathway in human lymphomas. Oncogene 30: 1831–1840.

34. Yan X, Lin Z, Chen F, Zhao X, Chen H, et al. (2009) Human BAMBI cooperates with Smad7 to inhibit transforming growth factor-beta signaling. J Biol Chem 284: 30097–30104.

35. Zavadil J, Cermak L, Soto-Nieves N, Bottinger EP (2004) Integration of TGF-beta/Smad and Jagged1/Notch signalling in epithelial-to-mesenchymal transition. EMBO J 23: 1155–1165.

36. Speleman F, De Preter K, Vandesompele J (2011) Neuroblastoma genetics and phenotype: a tale of heterogeneity. Semin Cancer Biol 21: 238–244.

37. Capasso M, Diskin SJ (2010) Genetics and genomics of neuroblastoma. Cancer Treat Res 155: 65–84.

38. Ohira M, Nakagawara A (2010) Global genomic and RNA profiles for novel risk stratification of neuroblastoma. Cancer Sci 101: 2295–2301.

39. Villamon E, Berbegall AP, Piqueras M, Tadeo I, Castel V, et al. (2013) Genetic instability and intratumoral heterogeneity in neuroblastoma with MYCN amplification plus 11q deletion. PLoS One 8: e53740.

40. Schatton T, Frank NY, Frank MH (2009) Identification and targeting of cancer stem cells. Bioessays 31: 1038–1049.

41. Kuiper RP, Ligtenberg MJ, Hoogerbrugge N, Geurts van Kessel A (2010) Germline copy number variation and cancer risk. Curr Opin Genet Dev 20: 282–289.

42. Dear PH (2009) Copy-number variation: the end of the human genome? Trends Biotechnol 27: 448–454.

43. Fanciulli M, Petretto E, Aitman TJ (2010) Gene copy number variation and common human disease. Clin Genet 77: 201–213.

44. Ghavifekr Fakhr M, Farshdousti Hagh M, Shanehbandi D, Baradaran B (2013) DNA Methylation Pattern as Important Epigenetic Criterion in Cancer. Genet Res Int 2013: 317569.

45. Nejman D, Straussman R, Steinfeld I, Ruvolo M, Roberts D, et al. (2014) Molecular rules governing de novo methylation in cancer. Cancer Res 74: 1475–1483.

46. Jones PA, Baylin SB (2007) The epigenomics of cancer. Cell 128: 683–692.

47. Liu Y, Cope L, Sun W, Wang Y, Prasad N, et al. (2013) DNA copy number variations characterize benign and malignant thyroid tumors. J Clin Endocrinol Metab 98: E558–566.

48. Kim KI, Kim TK, Kim IW, Ahn KS, Yoon SS, et al. (2012) Copy number variations in normal karyotype acute myeloid leukaemia and their association with treatment response. Basic Clin Pharmacol Toxicol 111: 317–324.

49. Jones PA (2014) At the tipping point for epigenetic therapies in cancer. J Clin Invest 124: 14–16.

50. Sigalotti L, Fratta E, Coral S, Maio M (2013) Epigenetic drugs as immunomodulators for combination therapies in solid tumors. Pharmacol Ther 142: 339–350.

51. Montenegro MF, Sanchez-Del-Campo L, Fernandez-Perez MP, Saez-Ayala M, Cabezas-Herrera J, et al. (2014) Targeting the epigenetic machinery of cancer cells. Oncogene

52. George RE, Lahti JM, Adamson PC, Zhu K, Finkelstein D, et al. (2010) Phase I study of decitabine with doxorubicin and cyclophosphamide in children with neuroblastoma and other solid tumors: a Children's Oncology Group study. Pediatr Blood Cancer 55: 629–638.

53. Stallings RL, Foley NH, Bryan K, Buckley PG, Bray I (2010) Therapeutic targeting of miRNAs in neuroblastoma. Expert Opin Ther Targets 14: 951–962.

54. Goff LA, Davila J, Swerdel MR, Moore JC, Cohen RI, et al. (2009) Ago2 immunoprecipitation identifies predicted microRNAs in human embryonic stem cells and neural precursors. PLoS One 4: e7192.

55. Jima DD, Zhang J, Jacobs C, Richards KL, Dunphy CH, et al. (2010) Deep sequencing of the small RNA transcriptome of normal and malignant human B cells identifies hundreds of novel microRNAs. Blood 116: e118–127.

56. Witten D, Tibshirani R, Gu SG, Fire A, Lui WO (2010) Ultra-high throughput sequencing-based small RNA discovery and discrete statistical biomarker analysis in a collection of cervical tumours and matched controls. BMC Biol 8: 58.

57. Coulon A, Flahaut M, Muhlethaler-Mottet A, Meier R, Liberman J, et al. (2011) Functional sphere profiling reveals the complexity of neuroblastoma tumor-initiating cell model. Neoplasia 13: 991–1004.

58. Wezel F, Pearson J, Kirkwood LA, Southgate J (2013) Differential expression of Oct4 variants and pseudogenes in normal urothelium and urothelial cancer. Am J Pathol 183: 1128–1136.

59. de Resende MF, Chinen LT, Vieira S, Jampietro J, da Fonseca FP, et al. (2013) Prognostication of OCT4 isoform expression in prostate cancer. Tumour Biol 34: 2665–2673.

60. Wang YD, Cai N, Wu XL, Cao HZ, Xie LL, et al. (2013) OCT4 promotes tumorigenesis and inhibits apoptosis of cervical cancer cells by miR-125b/BAK1 pathway. Cell Death Dis 4: e760.

61. Ezeh UI, Turek PJ, Reijo RA, Clark AT (2005) Human embryonic stem cell genes OCT4, NANOG, STELLAR, and GDF3 are expressed in both seminoma and breast carcinoma. Cancer 104: 2255–2265.

62. Lu Y, Zhu H, Shan H, Lu J, Chang X, et al. (2013) Knockdown of Oct4 and Nanog expression inhibits the stemness of pancreatic cancer cells. Cancer Lett 340: 113–123.

63. Tantin D (2013) Oct transcription factors in development and stem cells: insights and mechanisms. Development 140: 2857–2866.

64. Yang S, Zheng J, Ma Y, Zhu H, Xu T, et al. (2012) Oct4 and Sox2 are overexpressed in human neuroblastoma and inhibited by chemotherapy. Oncol Rep 28: 186–192.

65. Jerabek S, Merino F, Scholer HR, Cojocaru V (2013) OCT4: Dynamic DNA binding pioneers stem cell pluripotency. Biochim Biophys Acta 1839: 138–154.

66. Faunes F, Hayward P, Descalzo SM, Chatterjee SS, Balayo T, et al. (2013) A membrane-associated beta-catenin/Oct4 complex correlates with ground-state pluripotency in mouse embryonic stem cells. Development 140: 1171–1183.

67. Xin YH, Bian BS, Yang XJ, Cui W, Cui HJ, et al. (2013) POU5F1 enhances the invasiveness of cancer stem-like cells in lung adenocarcinoma by upregulation of MMP-2 expression. PLoS One 8: e83373.

68. Yang L, Zheng J, Xu T, Xiao X (2013) Downregulation of OCT4 promotes differentiation and inhibits growth of BE (2)-C human neuroblastoma I-type cells. Oncol Rep 29: 2191–2196.

69. Shi J, Shi W, Ni L, Xu X, Su X, et al. (2013) OCT4 is epigenetically regulated by DNA hypomethylation of promoter and exon in primary gliomas. Oncol Rep 30: 201–206.

70. Bailey JM, Leach SD (2012) Signaling pathways mediating epithelial- mesenchymal crosstalk in pancreatic cancer: Hedgehog, Notch and TGFbeta. In: Grippo PJ, Munshi HG, editors. Pancreatic Cancer and Tumor Microenvironment: Transworld Research Network, Trivandrum, India.

71. Bertrand FE, Angus CW, Partis WJ, Sigounas G (2012) Developmental pathways in colon cancer: crosstalk between WNT, BMP, Hedgehog and Notch. Cell Cycle 11: 4344–4351.

72. Wu WK, Wang XJ, Cheng AS, Luo MX, Ng SS, et al. (2013) Dysregulation and crosstalk of cellular signaling pathways in colon carcinogenesis. Crit Rev Oncol Hematol 86: 251–277.

73. Purow B (2012) Notch inhibition as a promising new approach to cancer therapy. Adv Exp Med Biol 727: 305–319.

74. Buttitta L, Mo R, Hui CC, Fan CM (2003) Interplays of Gli2 and Gli3 and their requirement in mediating Shh-dependent sclerotome induction. Development 130: 6233–6243.

75. McDermott A, Gustafsson M, Elsam T, Hui CC, Emerson CP Jr, et al. (2005) Gli2 and Gli3 have redundant and context-dependent function in skeletal muscle formation. Development 132: 345–357.

76. Nguyen V, Chokas AL, Stecca B, Ruiz i Altaba A (2005) Cooperative requirement of the Gli proteins in neurogenesis. Development 132: 3267–3279.

77. Guo X, Wang XF (2009) Signaling cross-talk between TGF-beta/BMP and other pathways. Cell Res 19: 71–88.

78. Kurpinski K, Lam H, Chu J, Wang A, Kim A, et al. (2010) Transforming growth factor-beta and notch signaling mediate stem cell differentiation into smooth muscle cells. Stem Cells 28: 734–742.

79. Wan F, Zhang S, Xie R, Gao B, Campos B, et al. (2010) The utility and limitations of neurosphere assay, CD133 immunophenotyping and side population assay in glioma stem cell research. Brain Pathol 20: 877–889.

80. Jiang Y, Boije M, Westermark B, Uhrbom L (2011) PDGF-B Can sustain self-renewal and tumorigenicity of experimental glioma-derived cancer-initiating cells by preventing oligodendrocyte differentiation. Neoplasia 13: 492–503.

81. Cournoyer S, Nyalendo C, Addioui A, Belounis A, Beaunoyer M, et al. (2012) Genotype analysis of tumor-initiating cells expressing CD133 in neuroblastoma. Genes Chromosomes Cancer 51: 792–804.

82. Yu AL, Gilman AL, Ozkaynak MF, London WB, Kreissman SG, et al. (2010) Anti-GD2 antibody with GM-CSF, interleukin-2, and isotretinoin for neuroblastoma. N Engl J Med 363: 1324–1334.

83. Barone G, Anderson J, Pearson AD, Petrie K, Chesler L (2013) New strategies in neuroblastoma: Therapeutic targeting of MYCN and ALK. Clin Cancer Res 19: 5814–5821.

84. Takebe N, Nguyen D, Yang SX (2013) Targeting Notch signaling pathway in cancer: Clinical development advances and challenges. Pharmacol Ther 141: 140–149.

85. Amakye D, Jagani Z, Dorsch M (2013) Unraveling the therapeutic potential of the Hedgehog pathway in cancer. Nat Med 19: 1410–1422.

86. Michaud NR, Wang Y, McEachern KA, Jordan JJ, Mazzola AM, et al. (2014) Novel Neutralizing Hedgehog Antibody MEDI-5304 Exhibits Antitumor Activity by Inhibiting Paracrine Hedgehog Signaling. Mol Cancer Ther 13: 386–398.

87. Xie J, Bartels CM, Barton SW, Gu D (2013) Targeting hedgehog signaling in cancer: research and clinical developments. Onco Targets Ther 6: 1425–1435.

88. Liu Z, Bandyopadhyay A, Nichols RW, Wang L, Hinck AP, et al. (2012) Blockade of Autocrine TGF-beta Signaling Inhibits Stem Cell Phenotype, Survival, and Metastasis of Murine Breast Cancer Cells. J Stem Cell Res Ther 2: 1–8.

89. Wakefield LM, Hill CS (2013) Beyond TGFbeta: roles of other TGFbeta superfamily members in cancer. Nat Rev Cancer 13: 328–341.

90. Schwab A, Stock C (2014) Ion channels and transporters in tumour cell migration and invasion. Philos Trans R Soc Lond B Biol Sci 369: 20130102.

91. Hoffmann EK, Lambert IH (2014) Ion channels and transporters in the development of drug resistance in cancer cells. Philos Trans R Soc Lond B Biol Sci 369: 20130109.

92. Becchetti A, Munaron L, Arcangeli A (2013) The role of ion channels and transporters in cell proliferation and cancer. Front Physiol 4: 312.

93. Arcangeli A, Becchetti A (2010) New Trends in Cancer Therapy: Targeting Ion Channels and Transporters. Pharmaceuticals 3: 1202–1224.

94. Romania P, Castellano A, Surace C, Citti A, De Ioris MA, et al. (2013) High-resolution array CGH profiling identifies Na/K transporting ATPase interacting 2 (NKAIN2) as a predisposing candidate gene in neuroblastoma. PLoS One 8: e78481.

95. Kunzelmann K (2005) Ion channels and cancer. J Membr Biol 205: 159–173.

96. Khaitan D, Sankpal UT, Weksler B, Meister EA, Romero IA, et al. (2009) Role of KCNMA1 gene in breast cancer invasion and metastasis to brain. BMC Cancer 9: 258.

97. Oeggerli M, Tian Y, Ruiz C, Wijker B, Sauter G, et al. (2012) Role of KCNMA1 in breast cancer. PLoS One 7: e41664.

98. Arcangeli A, Crociani O, Lastraioli E, Masi A, Pillozzi S, et al. (2009) Targeting ion channels in cancer: a novel frontier in antineoplastic therapy. Curr Med Chem 16: 66–93.

99. Pardo LA, Stuhmer W (2014) The roles of K(+) channels in cancer. Nat Rev Cancer 14: 39–48.

100. Ningaraj NS, Sankpal UT, Khaitan D, Meister EA, Vats TS (2009) Modulation of KCa channels increases anticancer drug delivery to brain tumors and prolongs survival in xenograft model. Cancer Biol Ther 8: 1924–1933.

101. Lang F, Stournaras C (2014) Ion channels in cancer: future perspectives and clinical potential. Philos Trans R Soc Lond B Biol Sci 369: 20130108.

102. Azimi I, Roberts-Thomson SJ, Monteith GR (2014) Calcium influx pathways in breast cancer: opportunities for pharmacological intervention. Br J Pharmacol 171: 945–960.

103. Kondratskyi A, Yassine M, Kondratska K, Skryma R, Slomianny C, et al. (2013) Calcium-permeable ion channels in control of autophagy and cancer. Front Physiol 4: 272.

104. Weaver AK, Liu X, Sontheimer H (2004) Role for calcium-activated potassium channels (BK) in growth control of human malignant glioma cells. J Neurosci Res 78: 224–234.

105. Rahman R, Heath R, Grundy R (2009) Cellular immortality in brain tumours: an integration of the cancer stem cell paradigm. Biochim Biophys Acta 1792: 280–288.