Mutated regions of nucleophosmin 1 elicit both CD4+ and CD8+ T-cell responses in patients with acute myeloid leukemia.

Mutations in the nucleophosmin gene (NPM1(mut)) are one of the most frequent molecular alterations in acute myeloid leukemia (AML), and immune responses may contribute to the favorable prognosis of AML patients with NPM1(mut). In the present study, we were able to demonstrate both CD4(+) and CD8(+) T-cell responses against NPM1(mut). Ten peptides derived from wild-type NPM1 and NPM1(mut) were subjected to ELISPOT analysis in 33 healthy volunteers and 27 AML patients. Tetramer assays against the most interesting epitopes were performed and Cr(51)-release assays were used to show the cytotoxicity of peptide-specific T cells. Moreover, HLA-DR-binding epitopes were used to test the role of CD4(+) T cells in NPM1 immunogenicity. Two epitopes (epitopes #1 and #3) derived from NPM1(mut) induced CD8(+) T-cell responses. A total of 33% of the NPM1(mut) AML patients showed immune responses against epitope #1 and 44% against epitope #3. Specific lysis of leukemic blasts was detected. To obtain robust immune responses against tumor cells, the activation of CD4(+) T cells is crucial. Therefore, overlapping (OL) peptides were analyzed in ELISPOT assays and OL8 was able to activate both CD8(+) and CD4(+) T cells. The results of the present study show that NPM1(mut) induces specific T-cell responses of CD4(+) and CD8(+) T cells and therefore is a promising target for specific immunotherapies in AML.

Zoledronate sensitizes endothelial cells to tumor necrosis factor-induced programmed cell death: evidence for the suppression of sustained activation of focal adhesion kinase and protein kinase B/Akt.

Bisphosphonates are potent inhibitors of osteoclast function widely used to treat conditions of excessive bone resorption, including tumor bone metastases. Recent evidence indicates that bisphosphonates have direct cytotoxic activity on tumor cells and suppress angiogenesis, but the associated molecular events have not been fully characterized. In this study we investigated the effects of zoledronate, a nitrogen-containing bisphosphonate, and clodronate, a non-nitrogen-containing bisphosphonate, on human umbilical vein endothelial cell (HUVEC) adhesion, migration, and survival, three events essential for angiogenesis. Zoledronate inhibited HUVEC adhesion mediated by integrin alphaVbeta3, but not alpha5beta1, blocked migration and disrupted established focal adhesions and actin stress fibers without modifying cell surface integrin expression level or affinity. Zoledronate treatment slightly decreased HUVEC viability and strongly enhanced tumor necrosis factor (TNF)-induced cell death. HUVEC treated with zoledronate and TNF died without evidence of enhanced annexin-V binding, chromatin condensation, or nuclear fragmentation and caspase dependence. Zoledronate inhibited sustained phosphorylation of focal adhesion kinase (FAK) and in combination with TNF, with and without interferon (IFN) gamma, of protein kinase B (PKB/Akt). Constitutive active PKB/Akt protected HUVEC from death induced by zoledronate and TNF/IFNgamma. Phosphorylation of c-Src and activation of NF-kappaB were not affected by zoledronate. Clodronate had no effect on HUVEC adhesion, migration, and survival nor did it enhanced TNF cytotoxicity. Taken together these data demonstrate that zoledronate sensitizes endothelial cells to TNF-induced, caspase-independent programmed cell death and point to the FAK-PKB/Akt pathway as a novel zoledronate target. These results have potential implications to the clinical use of zoledronate as an anti-angiogenic or anti-cancer agent.

Dramatic influence of V beta gene polymorphism on an antigen-specific CD8+ T cell response in vivo.

According to recent crystallographic studies, the TCR-alpha beta contacts MHC class I-bound antigenic peptides via the polymorphic V gene-encoded complementarity-determining region 1 beta (CDR1 beta) and the hypervariable (D)J-encoded CDR3 beta and CDR3 alpha domains. To evaluate directly the relative importance of CDR1 beta polymorphism on the fine specificity of T cell responses in vivo, we have taken advantage of congenic V beta a and V beta b mouse strains that differ by a CDR1 polymorphism in the V beta 10 gene segment. The V beta 10-restricted CD8+ T cell response to a defined immunodominant epitope was dramatically reduced in V beta a compared with V beta b mice, as measured either by the expansion of V beta 10+ cells or by the binding of MHC-peptide tetramers. These data indicate that V beta polymorphism has an important impact on TCR-ligand binding in vivo, presumably by modifying the affinity of CDR1 beta-peptide interactions.

Fission yeast Sec3 and Exo70 are transported on actin cables and localize the exocyst complex to cell poles.

The exocyst complex is essential for many exocytic events, by tethering vesicles at the plasma membrane for fusion. In fission yeast, polarized exocytosis for growth relies on the combined action of the exocyst at cell poles and myosin-driven transport along actin cables. We report here the identification of fission yeast Schizosaccharomyces pombe Sec3 protein, which we identified through sequence homology of its PH-like domain. Like other exocyst subunits, sec3 is required for secretion and cell division. Cells deleted for sec3 are only conditionally lethal and can proliferate when osmotically stabilized. Sec3 is redundant with Exo70 for viability and for the localization of other exocyst subunits, suggesting these components act as exocyst tethers at the plasma membrane. Consistently, Sec3 localizes to zones of growth independently of other exocyst subunits but depends on PIP(2) and functional Cdc42. FRAP analysis shows that Sec3, like all other exocyst subunits, localizes to cell poles largely independently of the actin cytoskeleton. However, we show that Sec3, Exo70 and Sec5 are transported by the myosin V Myo52 along actin cables. These data suggest that the exocyst holocomplex, including Sec3 and Exo70, is present on exocytic vesicles, which can reach cell poles by either myosin-driven transport or random walk.

Integrated Analysis of High-Resolution Gene Expression and Copy Number Profiling Identified Biallelic Deletion of CDKN2A/2B Tumor Suppressor Locus As the Most Frequent and Unique Genomic Abnormality in Diffuse Large B-Cell Lymphoma (DLBCL) with Strong Prognostic Value in Both GCB and ABC Subtypes and Not Overcome by a Dose-Intensive Immunochemotherapy Regimen Plus Rituximab. Results of a Prospective GELA Clinical Trial Program

Background and aim of the study: Genomic gains and losses play a crucial role in the development and progression of DLBCL and are closely related to gene expression profiles (GEP), including the germinal center B-cell like (GCB) and activated B-cell like (ABC) cell of origin (COO) molecular signatures. To identify new oncogenes or tumor suppressor genes (TSG) involved in DLBCL pathogenesis and to determine their prognostic values, an integrated analysis of high-resolution gene expression and copy number profiling was performed. Patients and methods: Two hundred and eight adult patients with de novo CD20+ DLBCL enrolled in the prospective multicentric randomized LNH-03 GELA trials (LNH03-1B, -2B, -3B, 39B, -5B, -6B, -7B) with available frozen tumour samples, centralized reviewing and adequate DNA/RNA quality were selected. 116 patients were treated by Rituximab(R)-CHOP/R-miniCHOP and 92 patients were treated by the high dose (R)-ACVBP regimen dedicated to patients younger than 60 years (y) in frontline. Tumour samples were simultaneously analysed by high resolution comparative genomic hybridization (CGH, Agilent, 144K) and gene expression arrays (Affymetrix, U133+2). Minimal common regions (MCR), as defined by segments that affect the same chromosomal region in different cases, were delineated. Gene expression and MCR data sets were merged using Gene expression and dosage integrator algorithm (GEDI, Lenz et al. PNAS 2008) to identify new potential driver genes. Results: A total of 1363 recurrent (defined by a penetrance > 5%) MCRs within the DLBCL data set, ranging in size from 386 bp, affecting a single gene, to more than 24 Mb were identified by CGH. Of these MCRs, 756 (55%) showed a significant association with gene expression: 396 (59%) gains, 354 (52%) single-copy deletions, and 6 (67%) homozygous deletions. By this integrated approach, in addition to previously reported genes (CDKN2A/2B, PTEN, DLEU2, TNFAIP3, B2M, CD58, TNFRSF14, FOXP1, REL...), several genes targeted by gene copy abnormalities with a dosage effect and potential physiopathological impact were identified, including genes with TSG activity involved in cell cycle (HACE1, CDKN2C) immune response (CD68, CD177, CD70, TNFSF9, IRAK2), DNA integrity (XRCC2, BRCA1, NCOR1, NF1, FHIT) or oncogenic functions (CD79b, PTPRT, MALT1, AUTS2, MCL1, PTTG1...) with distinct distribution according to COO signature. The CDKN2A/2B tumor suppressor locus (9p21) was deleted homozygously in 27% of cases and hemizygously in 9% of cases. Biallelic loss was observed in 49% of ABC DLBCL and in 10% of GCB DLBCL. This deletion was strongly correlated to age and associated to a limited number of additional genetic abnormalities including trisomy 3, 18 and short gains/losses of Chr. 1, 2, 19 regions (FDR < 0.01), allowing to identify genes that may have synergistic effects with CDKN2A/2B inactivation. With a median follow-up of 42.9 months, only CDKN2A/2B biallelic deletion strongly correlates (FDR p.value < 0.01) to a poor outcome in the entire cohort (4y PFS = 44% [32-61] respectively vs. 74% [66-82] for patients in germline configuration; 4y OS = 53% [39-72] vs 83% [76-90]). In a Cox proportional hazard prediction of the PFS, CDKN2A/2B deletion remains predictive (HR = 1.9 [1.1-3.2], p = 0.02) when combined with IPI (HR = 2.4 [1.4-4.1], p = 0.001) and GCB status (HR = 1.3 [0.8-2.3], p = 0.31). This difference remains predictive in the subgroup of patients treated by R-CHOP (4y PFS = 43% [29-63] vs. 66% [55-78], p=0.02), in patients treated by R-ACVBP (4y PFS = 49% [28-84] vs. 83% [74-92], p=0.003), and in GCB (4y PFS = 50% [27-93] vs. 81% [73-90], p=0.02), or ABC/unclassified (5y PFS = 42% [28-61] vs. 67% [55-82] p = 0.009) molecular subtypes (Figure 1). Conclusion: We report for the first time an integrated genetic analysis of a large cohort of DLBCL patients included in a prospective multicentric clinical trial program allowing identifying new potential driver genes with pathogenic impact. However CDKN2A/2B deletion constitutes the strongest and unique prognostic factor of chemoresistance to R-CHOP, regardless the COO signature, which is not overcome by a more intensified immunochemotherapy. Patients displaying this frequent genomic abnormality warrant new and dedicated therapeutic approaches.

Generation of short-term murine natural killer cell clones to analyze Ly49 gene expression.

Natural killer (NK) cellsexpress receptors specific for class I major histocompatibility complex (MHC) molecules. In the mouse, the class I specific receptors identified to date belong to the polymorphic Ly49 receptor family. Engagement of Ly49 receptors with their respective MHC ligands results in negative regulation of NK cell effector functions, consistent with a critical role of these receptors in "missing self" recognition. The Ly49 receptors analyzed so far are clonally distributed such that multiple distinct Ly49 receptors can be expressed by individual NK cells (for review see refs. 1-3). The finding that most NK cells that express the Ly49A receptor do so from a single Ly49A allele (whereby expression can occur from the maternal or the paternal chromosome) may thus reflect a putative receptor distribution process that restricts the number of Ly49 receptors expressed in a single NK cell (3-5).

BIRO1, a cell-permeable BH3 peptide, promotes mitochondrial fragmentation and death of retinoblastoma cells.

Retinoblastoma is the most common pediatric intraocular neoplasm. While retinoblastoma development requires the inactivation of both alleles of the retinoblastoma tumor suppressor gene (RB1) in the developing retina, additional genomic changes are involved in tumor progression, which progressively lead to resistance of tumor cells to death. Therapeutics acting at very downstream levels of death signaling pathways should therefore be interesting in killing retinoblastoma cells. The BH3-only proteins promote apoptosis by modulating the interaction between the pro- and antiapoptotic members of the BCL2 protein family, and this effect can be recapitulated by the BH3 domains. This report analyzes the effect of various BH3 peptides, corresponding to different BH3-only proteins, on two retinoblastoma cell lines, Y79 and WERI-Rb, as well as on the photoreceptor cell line 661W. The BH3 peptide BIRO1, derived from the BCL2L11 death domain, was very effective in promoting Y79 and WERI-Rb cell death without affecting the 661W photoreceptor cells. This cell death was efficient even in absence of BAX and was shown to be caspase independent. While ROS production or AIF release was not detected from mitochondria of treated cells, BIRO1 initiated mitochondria fragmentation in a short period of time following treatment. IMPLICATIONS: The BIRO1 peptide is highly effective at killing retinoblastoma cells and has potential as a peptidomimetic.

Expression of common acute lymphoblastic leukemia antigen (CALLA) on human malignant melanoma cell lines.

Fifteen human melanoma cells lines were tested by an antibody-binding radioimmunoassay using a monoclonal antibody (A12) directed against the common acute lymphoblastic leukemia antigen (CALLA). Cells from six melanoma lines were found to react with this antibody. The level of antigen and the percentage of positive cells in these six melanoma lines showed wide variation, as demonstrated by analysis in the fluorescence-activated cell sorter (FACS). Immunoprecipitation of solubilized 125I-labeled membrane proteins from CALLA positive melanoma cells with A12 monoclonal antibody revealed a major polypeptide chain with an apparent m.w. of 100,000 daltons, characteristic for CALLA as determined on SDS-polyacrylamide gel electrophoresis. The expression of CALLA on MP-6 melanoma cells was modulated when the cells were cultured in the presence of A12 antibody. Reexpression of CALLA on these cells occurred within 5 days after transfer of the modulated cells into medium devoid of monoclonal antibody.

Clinical, molecular, and prognostic significance of WHO type inv(3)(q21q26.2)/t(3;3)(q21;q26.2) and various other 3q abnormalities in acute myeloid leukemia.

PURPOSE: Acute myeloid leukemia (AML) with inv(3)(q21q26.2)/t(3;3)(q21;q26.2) [inv(3)/t(3;3)] is recognized as a distinctive entity in the WHO classification. Risk assignment and clinical and genetic characterization of AML with chromosome 3q abnormalities other than inv(3)/t(3;3) remain largely unresolved. PATIENTS AND METHODS: Cytogenetics, molecular genetics, therapy response, and outcome analysis were performed in 6,515 newly diagnosed adult AML patients. Patients were treated on Dutch-Belgian Hemato-Oncology Cooperative Group/Swiss Group for Clinical Cancer Research (HOVON/SAKK; n = 3,501) and German-Austrian Acute Myeloid Leukemia Study Group (AMLSG; n = 3,014) protocols. EVI1 and MDS1/EVI1 expression was determined by real-time quantitative polymerase chain reaction. RESULTS: 3q abnormalities were detected in 4.4% of AML patients (288 of 6,515). Four distinct groups were defined: A: inv(3)/t(3;3), 32%; B: balanced t(3q26), 18%; C: balanced t(3q21), 7%; and D: other 3q abnormalities, 43%. Monosomy 7 was the most common additional aberration in groups (A), 66%; (B), 31%; and (D), 37%. N-RAS mutations and dissociate EVI1 versus MDS1/EVI1 overexpression were associated with inv(3)/t(3;3). Patients with inv(3)/t(3;3) and balanced t(3q21) at diagnosis presented with higher WBC and platelet counts. In multivariable analysis, only inv(3)/t(3;3), but not t(3q26) and t(3q21), predicted reduced relapse-free survival (hazard ratio [HR], 1.99; P < .001) and overall survival (HR, 1.4; P = .006). This adverse prognostic impact of inv(3)/t(3;3) was enhanced by additional monosomy 7. Group D 3q aberrant AML also had a poor outcome related to the coexistence of complex and/or monosomal karyotypes and cryptic inv(3)/t(3;3). CONCLUSION: Various categories of 3q abnormalities in AML can be distinguished according to their clinical, hematologic, and genetic features. AML with inv(3)/t(3;3) represents a distinctive subgroup with unfavorable prognosis.

Geometric control of the cell cycle.

How do cells sense their own size and shape? And how does this information regulate progression of the cell cycle? Our group, in parallel to that of Paul Nurse, have recently demonstrated that fission yeast cells use a novel geometry-sensing mechanism to couple cell length perception with entry into mitosis. These rod-shaped cells measure their own length by using a medially-placed sensor, Cdr2, that reads a protein gradient emanating from cell tips, Pom1, to control entry into mitosis. Budding yeast cells use a similar molecular sensor to delay entry into mitosis in response to defects in bud morphogenesis. Metazoan cells also modulate cell proliferation in response to their own shape by sensing tension. Here I discuss the recent results obtained for the fission yeast system and compare them to the strategies used by these other organisms to perceive their own morphology.

Molecular modeling of peptide-MHC class I complexes : applications to peptide based immunotherapy

Summary The specific CD8+ T cell immune response against tumors relies on the recognition by the T cell receptor (TCR) on cytotoxic T lymphocytes (CTL) of antigenic peptides bound to the class I major histocompatibility complex (MHC) molecule. Such tumor associated antigenic peptides are the focus of tumor immunotherapy with peptide vaccines. The strategy for obtaining an improved immune response often involves the design of modified tumor associated antigenic peptides. Such modifications aim at creating higher affinity and/or degradation resistant peptides and require precise structures of the peptide-MHC class I complex. In addition, the modified peptide must be cross-recognized by CTLs specific for the parental peptide, i.e. preserve the structure of the epitope. Detailed structural information on the modified peptide in complex with MHC is necessary for such predictions. In this thesis, the main focus is the development of theoretical in silico methods for prediction of both structure and cross-reactivity of peptide-MHC class I complexes. Applications of these methods in the context of immunotherapy are also presented. First, a theoretical method for structure prediction of peptide-MHC class I complexes is developed and validated. The approach is based on a molecular dynamics protocol to sample the conformational space of the peptide in its MHC environment. The sampled conformers are evaluated using conformational free energy calculations. The method, which is evaluated for its ability to reproduce 41 X-ray crystallographic structures of different peptide-MHC class I complexes, shows an overall prediction success of 83%. Importantly, in the clinically highly relevant subset of peptide-HLAA*0201 complexes, the prediction success is 100%. Based on these structure predictions, a theoretical approach for prediction of cross-reactivity is developed and validated. This method involves the generation of quantitative structure-activity relationships using three-dimensional molecular descriptors and a genetic neural network. The generated relationships are highly predictive as proved by high cross-validated correlation coefficients (0.78-0.79). Together, the here developed theoretical methods open the door for efficient rational design of improved peptides to be used in immunotherapy. Résumé La réponse immunitaire spécifique contre des tumeurs dépend de la reconnaissance par les récepteurs des cellules T CD8+ de peptides antigéniques présentés par les complexes majeurs d'histocompatibilité (CMH) de classe I. Ces peptides sont utilisés comme cible dans l'immunothérapie par vaccins peptidiques. Afin d'augmenter la réponse immunitaire, les peptides sont modifiés de façon à améliorer l'affinité et/ou la résistance à la dégradation. Ceci nécessite de connaître la structure tridimensionnelle des complexes peptide-CMH. De plus, les peptides modifiés doivent être reconnus par des cellules T spécifiques du peptide natif. La structure de l'épitope doit donc être préservée et des structures détaillées des complexes peptide-CMH sont nécessaires. Dans cette thèse, le thème central est le développement des méthodes computationnelles de prédiction des structures des complexes peptide-CMH classe I et de la reconnaissance croisée. Des applications de ces méthodes de prédiction à l'immunothérapie sont également présentées. Premièrement, une méthode théorique de prédiction des structures des complexes peptide-CMH classe I est développée et validée. Cette méthode est basée sur un échantillonnage de l'espace conformationnel du peptide dans le contexte du récepteur CMH classe I par dynamique moléculaire. Les conformations sont évaluées par leurs énergies libres conformationnelles. La méthode est validée par sa capacité à reproduire 41 structures des complexes peptide-CMH classe I obtenues par cristallographie aux rayons X. Le succès prédictif général est de 83%. Pour le sous-groupe HLA-A*0201 de complexes de grande importance pour l'immunothérapie, ce succès est de 100%. Deuxièmement, à partir de ces structures prédites in silico, une méthode théorique de prédiction de la reconnaissance croisée est développée et validée. Celle-ci consiste à générer des relations structure-activité quantitatives en utilisant des descripteurs moléculaires tridimensionnels et un réseau de neurones couplé à un algorithme génétique. Les relations générées montrent une capacité de prédiction remarquable avec des valeurs de coefficients de corrélation de validation croisée élevées (0.78-0.79). Les méthodes théoriques développées dans le cadre de cette thèse ouvrent la voie du design de vaccins peptidiques améliorés.

Vascular continuity, cell axialization and auxin.

Phytohormones have been implicated in vascular development in various ways, but their precise function and the extent of their influence is still controversial. Recent results from experimental manipulation of developing organs and Arabidopsis developmental genetics support a role for polar auxin flow in cell axis formation within the vascular system and, interestingly, also in the embryonic establishment of the plant body axis. Vascular responses to auxin transport inhibition indicate patterns of auxin distribution during leaf development and new technologies may enable these predictions to be tested within the near future. Moreover, recently discovered Arabidopsis axialisation mutants seem to identify essential genes that relay auxin signals in vascular development. A first gene in this class, MONOPTEROS (MP) has been cloned and encodes a transcription factor capable of binding to auxin response elements in the control regions of auxin regulated genes. Molecular identification of further axialisation genes may provide access to a mechanistic understanding of plant cell axis formation.

Identification of novel and cell type enriched cofactors of the transcription activation domain of RelA (p65 NF-kappaB).

RelA (NF-kappaB) is a transcription factor inducible by distinct stimuli in many different cell types. To find new cell type specific cofactors of NF-kappaB dependent transcription, we isolated RelA transcription activation domain binding proteins from the nuclear extracts of three different cell types. Analysis by electrophoresis and liquid chromatography tandem mass spectrometry identified several novel putative molecular partners. Some were strongly enriched in the complex formed from the nuclear extracts of specific cell types.

Molecular imaging of matrix metalloproteinases in atherosclerotic plaques.

Ischaemic stroke and myocardial infarction often result from the sudden rupture of an atherosclerotic plaque. The subsequent arterial thrombosis occluding the vessel lumen has been widely indicated as the crucial acute event causing peripheral tissue ischaemia. A complex cross-talk between systemic and intraplaque inflammatory mediators has been shown to regulate maturation, remodeling and final rupture of an atherosclerotic plaque. Matrix metalloproteinases (MMPs) are proteolytic enzymes (released by several cell subsets within atherosclerotic plaques), which favour atherogenesis and increase plaque vulnerability. Thus, the assessment of intraplaque levels and activity of MMP might be of pivotal relevance in the evaluation of the risk of rupture. New imaging approaches, focused on the visualisation of inflammation in the vessel wall and plaque, may emerge as tools for individualised risk assessment and prevention of events. In this review, we summarize experimental findings of the currently available invasive and noninvasive imaging techniques, used to detect the presence and activity of MMPs in atherosclerotic plaques.

Molecular and Functional Characterization of Neuroblastoma-Initiating Cells

Neuroblastoma (NB) is the most common extracranial malignant tumor in young children and arises at any site of the sympathetic nervous system. The disease exhibits a remarkable phenotypic diversity ranging from spontaneous regression to fatal disease. Poor outcome results from a rapidly progressive, metastatic and drug-resistant disease. Recent studies have suggested that solid tumors may arise from a minor population of cancer stem cells (CSCs) with stem cell markers and typical properties such as self-renewal ability, asymmetric division and drug resistance. In this model, CSCs possess the exclusive ability to initiate and maintain the tumor, and to produce distant metastases. Tumor cell subpopulations with stem-like phenotypes have indeed been identified in several cancer including leukemia, breast, brain and colon cancers. CSC hypothesis still needs to be validated in the other cancers including NB.NB originates from neural crest-derived malignant sympatho-adrenal cells. We have identified rare cells that express markers in conformity with neural crest stem cells and their derived lineages within primary NB tissue and cell lines, leading us to postulate the existence of CSCs in NB tumors.In the absence of specific markers to isolate CSCs, we adapted to NB tumor cells the sphere functional assay, based on the ability of stem cells to grow as spheres in non-adherent conditions. By serial passages of spheres from bone marrow NB metastases, a subset of cells was gradually selected and its specific gene expression profile identified by micro-array time-course analysis. The differentially expressed genes in spheres are enriched in genes implicated in development including CD133, ABC-transporters, WNT and NOTCH genes, identified in others solid cancers as CSCs markers, and other new markers, all referred by us as the Neurosphere Expression Profile (NEP). We confirmed the presence of a cell subpopulation expressing a combination of the NEP markers within a few primary NB samples.The tumorigenic potential of NB spheres was assayed by in vivo tumor growth analyses using orthotopic (adrenal glands) implantations of tumor cells into immune-compromised mice. Tumors derived from the sphere cells were significantly more frequent and were detected earlier compared to whole tumor cells. However, NB cells expressing the neurosphere-associated genes and isolated from the bulk tumors did not recapitulate the CSC-like phenotype in the orthotopic model. In addition, the NB sphere cells lost their higher tumorigenic potential when implanted in a subcutaneous heterotopic in vivo model.These results highlighted the complex behavior of CSC functions and led us to consider the stem-like NB cells as a dynamic and heterogeneous cell population influenced by microenvironment signals.Our approach identified for the first time candidate genes that may be associated with NB self-renewal and tumorigenicity and therefore would establish specific functional targets for more effective therapies in aggressive NB.