Skip to main content


Molecular profiling of malignant peripheral nerve sheath tumors associated with neurofibromatosis type 1, based on large-scale real-time RT-PCR



Neurofibromatosis type 1 (NF1) is an autosomal dominant disorder with a complex range of clinical symptoms. The hallmark of NF1 is the onset of heterogeneous (dermal or plexiform) benign neurofibromas. Plexiform neurofibromas can give rise to malignant peripheral nerve sheath tumors (MPNSTs), and the underlying molecular mechanisms are largely unknown.


To obtain further insight into the molecular pathogenesis of MPNSTs, we used real-time quantitative RT-PCR to quantify the mRNA expression of 489 selected genes in MPNSTs, in comparison with plexiform neurofibromas.

The expression of 28 (5.7%) of the 489 genes was significantly different between MPNSTs and plexiform neurofibromas; 16 genes were upregulated and 12 were downregulated in MPNSTs.

The altered genes were mainly involved in cell proliferation (MKI67, TOP2A, CCNE2), senescence (TERT, TERC), apoptosis (BIRC5/Survivin, TP73) and extracellular matrix remodeling (MMP13, MMP9, TIMP4, ITGB4). More interestingly, other genes were involved in the Ras signaling pathway (RASSF2, HMMR/RHAMM) and the Hedgehog-Gli signaling pathway (DHH, PTCH2). Several of the down-regulated genes were Schwann cell-specific (L1CAM, MPZ, S100B, SOX10, ERBB3) or mast cell-specific (CMA1, TPSB), pointing to a depletion and/or dedifferentiation of Schwann cells and mast cells during malignant transformation of plexiform neurofibromas.


These data suggest that a limited number of signaling pathways, and particularly the Hedgehog-Gli signaling pathway, may be involved in malignant transformation of plexiform neurofibromas. Some of the relevant genes or their products warrant further investigation as potential therapeutic targets in NF1.


Neurofibromatosis type 1 (NF1) is an autosomal dominant neurocutaneous disorder affecting 1 in 3000 individuals worldwide [1]. The NF1 gene, located on chromosome 17q11.2, was identified by positional cloning, and its protein product, neurofibromin, functions as a tumor suppressor [2, 3]. Neurofibromin contains a central domain homologous to the Ras-GTPase-activating protein family (Ras-GAPs), which function as negative regulators of Ras proteins [4].

The main clinical features of NF1 are café au lait macules, skinfold freckling and iris Lisch nodules. Patients are at an increased risk of both benign and malignant tumors, and NF1 is thus classified as a tumor predisposition syndrome. The most common tumors are benign peripheral nerve sheath tumors (neurofibromas), which vary greatly in both number and size, and may be dermal or plexiform [5]. In contrast to dermal neurofibromas, which are typically small and grow as discrete lesions in the dermis, plexiform neurofibromas can develop internally along the plexus of major peripheral nerves and become quite large [6]. Both dermal and plexiform neurofibromas are heterogeneous tumors mainly composed of Schwann cells (60–80%), together with neurons, fibroblasts, mast cells and other cells.

About 5% of patients with NF1, neurofibromas (mainly plexiform neurofibromas) progress to malignant peripheral nerve sheath tumors (MPNSTs). More than 80% of MPNSTs are high-grade malignant tumors, corresponding to WHO grade III-IV [7]. MPNSTs are resistant to conventional therapies, and their deep-seated position and locally invasive growth hinder complete surgical resection. The 5-year survival rate among patients with MPNSTs ranges from 30–50%. Schwann cells are considered to be the progenitors of both neurofibromas and MPNSTs, but recent data suggest that other cell types may contribute to the development of these tumors [8].

The molecular mechanisms responsible for malignant progression of neurofibromas are largely unknown. Only a few relevant genetic alterations have so far been identified [9, 10]. In keeping with its role as a classical tumor suppressor gene, NF1 loss of heterozygosity (LOH) has been found in NF1-associated MPNSTs (but also in benign neurofibromas) [11]. TP53 mutations have been identified in MPNSTs but not in benign neurofibromas, indicating that the p53-mediated pathway is involved in tumor progression [1214]. Consistent with the role of p53 in the progression of MPNSTs, mice that harbor both NF1 and TP53 mutations develop MPNSTs [15]. Alterations of other genes (p16/CDKN2A, p14/ARF, p27/KIP1, EGFR) are frequent in MPNSTs but not in neurofibromas [1620]. Taken together, these studies suggest that loss of NF1 initiates tumor formation, and that malignant progression requires additional genetic lesions.

The recent development of efficient tools for large-scale analysis of gene expression has provided new insights into the involvement of gene networks and regulatory pathways in various tumoral processes [21]. These methods include cDNA microarrays, which can be used to analyze the expression of thousands of genes at a time, and real-time RT-PCR assays for more accurate and quantitative expression analysis of smaller numbers of candidate genes [22].

To obtain further insight into the molecular pathogenesis of MPNSTs, we used real-time quantitative RT-PCR to quantify the mRNA expression of a large number of selected genes in pooled MPNST samples, in comparison with pooled plexiform neurofibroma samples. We assessed the expression level of 489 genes involved in various cellular and molecular phenomena associated with tumorigenesis. We focused particularly on the expression of genes related to MPNSTs, and genes expressed during Schwann cell differentiation. Forty genes of interest were further investigated in nine individual MPNSTs (all arising from plexiform neurofibromas), in comparison with 14 plexiform neurofibromas.


We first quantified the mRNA expression level of the 489 genes (see list in annex; Additional file 1) in an MPNST pool, a plexiform neurofibroma pool and a dermal neurofibroma pool. We then selected for further study those genes whose expression in the MPNST pool differed markedly (> 10-fold) from that in both the plexiform neurofibroma pool and the dermal neurofibroma pool. These robust selection criteria (a cutoff of 10-fold expression difference in the MPNST pool) ensure identification of gene with marked interest. The mRNA expression of the genes thus identified was then determined in 9 individual MPNSTs (Table 1) in comparison with 14 plexiform neurofibromas.

Table 1 Clinical and histological characteristics of the 9 patients with MPNST

mRNA expression of 489 genes in the MPNST pool relative to the plexiform and dermal neurofibroma pools

The MPNST, plexiform neurofibroma and dermal neurofibroma pools were each prepared by mixing identical amounts of tumor RNA from four patients. The mean TBP gene Ct (threshold cycle) values for the four tumor samples were 26.18 ± 0.39 (MPNST pool), 26.23 ± 1.45 (plexiform neurofibroma pool) and 26.30 ± 0.47 (dermal neurofibroma pool).

Very low levels of target gene mRNA, that were only detectable but not reliably quantifiable by means of real-time quantitative RT-PCR assays, mainly based on fluorescence SYBR Green methodology (Ct >32), were observed for 50 (10.2%) of the 489 genes in the MPNST, plexiform and dermal neurofibroma pools.

Forty (9.1%) of the 439 remaining genes were expressed at a different level (> 10-fold) in the MPNST pool compared to both the dermal neurofibroma and plexiform neurofibroma pools; 27 (6.1%) genes were upregulated and 13 (3.0%) were down-regulated. The NF1 gene expression level did not differ between the dermal neurofibroma, plexiform neurofibroma and MPNST pools: the NNF1 values (see Patients and Methods) were 1.00, 1.39 and 1.45, respectively. This lack of difference was probably due to the fact that the NF1 gene is ubiquitously expressed and therefore expressed in the different cell components of the neurofibromas and MPNSTs, and that NF1-/- Schwann cells represent only a fraction of the total Schwann cell population in tumor samples.

mRNA expression of the 27 upregulated genes in 9 MPNSTs and 14 plexiform neurofibromas

The expression level of the 27 upregulated genes identified by pooled sample analysis was then determined individually in 9 MPNSTs and 14 plexiform neurofibromas. Sixteen (59.3%) of the 27 genes were significantly upregulated in the 9 MPNSTs (P < 0.05; Table 2).

Table 2 List of the significantly up-regulated genes in the MPNSTs relative to the plexiform neurofibromas.

The 16 upregulated genes were mainly involved in cell proliferation (MKI67, TOP2A, CCNE2), senescence (TERT, TERC/hTR), apoptosis (BIRC5/Survivin, TP73) and extracellular matrix remodeling (MMP13, MMP9).

The capacity of each of these 16 genes to discriminate between MPNSTs and plexiform neurofibromas was then tested by ROC curve analysis. The overall diagnostic value of the 16 molecular markers was assessed in terms of the AUCs (Table 2). Figure 1 shows the mRNA levels of the three most discriminatory genes, namely MKI67 (AUC-ROC, 1.000), BIRC5/Survivin (AUC-ROC, 0.984), and SPP1/Osteopontin (AUC-ROC, 0.984), in each MPNST and plexiform neurofibroma sample. For information, Figure 1 also shows the mRNA levels of these three genes in 10 dermal neurofibromas.

Figure 1

mRNA levels of MKI67, BIRC5/Survivin and SPP1 in 10 individual dermal neurofibromas (white bars), 14 plexiform neurofibromas (gray bars) and 9 MPNSTs (black bars). Median values (and ranges) are indicated for each tumor subgroup.

mRNA expression of the 13 down-regulated genes in 9 MPNSTs and 14 plexiform neurofibromas

Twelve (92.3%) of the 13 genes were significantly down-regulated in the 9 MPNSTs (P < 0.05; Table 3).

Table 3 List of the significantly down-regulated genes in the MPNSTs relative to the plexiform neurofibromas

The 12 down-regulated genes mainly were cell type-specific, and included Schwann cell-specific genes (L1CAM, MPZ, S100B, SOX10) and mast cell-specific genes (CMA1, TPSB). The others down-regulated genes were involved in extracellular matrix remodeling (ITGB4, TIMP4) and in the Hedgehog-Gli signaling pathway (DHH, PTCH2).

The capacity of each of these 12 genes to discriminate between MPNSTs and plexiform neurofibromas was then tested by ROC curve analysis. The overall diagnostic value of the 12 molecular markers was assessed in terms of the AUCs (Table 3). Figure 2 shows the mRNA levels of the three most discriminatory genes, namely ITGB4 (AUC-ROC, 1.000), CMA1/Chymase 1 (AUC-ROC, 1.000), and L1CAM (AUC-ROC, 1.000), in each MPNST and plexiform neurofibroma sample, and also in each dermal neurofibroma sample.

Figure 2

mRNA levels of ITGB4, CMA1/Chymase 1 and L1CAM in 10 individual dermal neurofibromas (white bars), 14 plexiform neurofibromas (gray bars) and 9 MPNSTs (black bars). Median values (and ranges) are indicated for each tumor subgroup.

The mRNA levels indicated in Tables 2 and 3 (calculated as described in Materials and Methods) are expressed relative to the endogenous control TBP mRNA level, to control for the starting amount and quality of total RNA. Similar results were obtained with a second endogenous control, RPLP0 (also known as 36B4). Indeed, the 16 upregulated genes and the 12 down-regulated genes were also significantly up-regulated or down-regulated in the MPNSTs relative to the plexiform neurofibromas.


We used real-time quantitative RT-PCR to quantify the mRNA expression of 489 selected genes in pooled MPNST samples, in comparison with pooled plexiform neurofibroma and pooled dermal neurofibroma samples. Forty genes of interest were then investigated in 9 individual MPNSTs and 14 plexiform neurofibromas. Comparison of pool values with the mean of the corresponding individual values showed that RNA pooling was an appropriate initial screening approach, significantly limiting the required number of PCR experiments. Using the same approach, we have also shown the involvement of several molecular pathways in the genesis of plexiform neurofibroma [23], and other human diseases [24].

Real-time quantitative RT-PCR is a promising alternative to cDNA microarrays for molecular tumor profiling, being far more precise, reproducible and quantitative. Real-time RT-PCR is more useful for analyzing weakly expressed genes (such as TERT, CXCL5, FOXM1 and FOXA2/HNF3B in the present study). Finally, real-time RT-PCR requires smaller amounts of total RNA (about 2 ng per target gene), and is therefore suitable for analyzing small tumors (as in the present study) or microdissected samples.

We included a number of genes known to be involved in various cellular and molecular mechanisms associated with tumorigenesis, and known to be altered (mainly at the transcriptional level) in various cancers. These genes encode proteins involved in cell cycle control, cell-cell interactions, signal transduction pathways, apoptosis, angiogenesis, etc. (about 10–20 genes were selected per pathway). After scrutinizing the literature, we also included most genes reported to be involved in neurofibromas and MPNSTs, and also genes expressed during Schwann cell differentiation.

Among the 489 genes analyzed, 28 (5.7%) showed a significantly different level of expression in MPNSTs relative to plexiform neurofibromas, suggesting that several signaling pathways are specifically involved in MPNST (Tables 2 and 3).

Some of our results support those reported in the literature on NF1-associated MPNSTs. First, several genes associated with cell cycle control (MKI67, TOP2A, CCNE2) were over-expressed in MPNSTs, suggesting higher cell proliferation rates than in plexiform neurofibromas, as previously reported by Kindblom et al. [25]. MKI67 encodes Ki-67, a large protein of unknown function used as a classical histopathological marker of cell proliferation. Other proliferation-associated genes (CCND1, CCNE1, etc) were also upregulated in MPNSTs, but less markedly than MKI67, TOP2A and CCNE2. Second, most of mast cell-specific genes (CMA1/Chymase 1 and TPSB/Tryptase beta) and Schwann cell-specific genes (L1CAM, MPZ, S100B, SOX10) tested here were markedly under-expressed in MPNSTs, probably owing either to a lower abundance of a particular cell type (likely mast cells) in MPNSTs relative to plexiform neurofibromas, or to dedifferentiation of a particular cell type (likely Schwann cells). Schwann cell dedifferentiation could also explain the observed under-expression of ITGB4 and ERBB3, key genes in the Schwann cell lineage [26, 27]. DeClue et al. [19] suggested that, following acquisition of ERBB1/EGFR overexpression, Schwann cells may begin to dedifferentiate and lose markers such as ERBB3. In agreement with these data, we observed slight over-expression of ERBB1/EGFR in MPNSTs, whereas the two other ERBB family members (ERBB2 and ERBB4) were normally expressed (Figure 3A). It is noteworthy that, in Schwann cells, ERBB3 could be regulated by SOX10 (found here to be under-expressed in MPNSTs), as in peripheral glial cells [28].

Figure 3

Expression level of various genes in the dermal neurofibroma pool (white bars), the plexiform neurofibroma pool (gray bars) and the MPNST pool (black bars). A, ERBB family member genes; B, matrix metalloproteinase genes; C, inhibitor of metalloproteinase genes; D, chemokine genes; E, Ras signaling pathway genes; F, Hedgehog-Gli signaling pathway genes; G, Gli-regulated genes.

Some of our results for MPNSTs are new, but are in keeping with general concepts of tumorigenesis, in which altered genes are involved in senescence, apoptosis and extracellular matrix remodeling. First, TERT and TERC/hTR, the two main components of human telomerase, were upregulated in MPNSTs. The expression patterns of TERT and TERC/hTR were very similar (r = +0.793, P = 0.00001; Spearman rank correlation test). TERT (human telomerase revere transcriptase) is the rate-limiting factor for telomerase activity, whereas TERC/hTR (human telomerase RNA) also seems to correlate to a certain extent with telomerase reactivation [24]. TERT and TERC/hTR are upregulated in almost all malignancies but not in benign tumors [29].

We also found that two genes involved in apoptosis – BIRC5/Survivin and TP73 – were upregulated in MPNSTs. BIRC5/Survivin encodes an antiapoptotic protein overexpressed in common human cancers [30]. It is noteworthy that coexpression of BIRC5/Survivin and TERT transcripts is associated with a high risk of tumor-related death among patients with soft-tissue sarcomas [31]. TP73, the second apoptosis-related candidate gene identified here, encodes two different proteins that are expressed under the control of two independent promoters and have opposite activities, namely the transcriptionally active full-length protein TAp73, and the amino-terminally truncated dominant-negative protein Δ Np73 [32]. Unlike TP53, TP73 is mainly regulated at the transcriptional level. TAp73 induces cell-cycle arrest and apoptosis, whereas ΔNp73 inhibits both TAp73-induced and p53-induced apoptosis. Furthermore, ΔNp73 is induced by TAp73 and p53, in a dominant-negative feedback loop that regulates p53 and p73 function [32]. Interesting, we found that TAp73 transcription was significantly upregulated in MPNSTs (P < 0.01; AUC-ROC, 0.952), while ΔNp73 transcription was significantly down-regulated (P < 0.01; AUC-ROC, 0.917) (data not shown).

MPNSTs showed upregulation of two matrix metalloproteinase genes (MMP9 and MMP13) and down-regulation of a tissue inhibitor of metalloproteinases (TIMP4). MMPs are associated with extracellular matrix turnover, a process that plays a very active role in tumor invasion and metastasis [33]. For their part, TIMPs play a key role in extracellular matrix homeostasis by regulating MMP activity [34]. Many reports suggest that an MMP-TIMP imbalance may contribute to the malignant phenotype [33, 34]. The other extracellular matrix-related genes tested in the present studies showed no statistically differentially expression levels between MPNSTs and plexiform neurofibromas (Figures 3B and 3C). We have previously identified MMP9 as a major upregulated gene in plexiform neurofibromas compared to dermal neurofibromas [23] (Figure 3B), suggesting that MMP9 upregulation may be an early event in MPNST tumorigenesis.

Finally, our results suggest that inappropriate activation of molecular signaling pathways occurs specifically in NF1-associated MPNST or in limited human tumors, particularly brain and skin tumors. MPNSTs showed upregulation of CXCL5, the gene encoding the chemokine ENA-78 (epithelial cell-derived neutrophil-activating peptide) [35]. Chemokines play fundamental roles in the development, homeostasis and functioning of the immune system, but also have effects on endothelial cells and cells of the nervous system [36]. CXCL5 is physically linked, in chromosomal region 4q12-q13, to other chemokine genes with markedly similar nucleotide sequences, including CXCL6, IL8 and GRO1/CXCL1 [37]. We found that CXCL6 was slightly upregulated (trend toward statistical significance) in MPNSTs. We have previously detected IL8 and GRO1/CXCL1 upregulation in the early stages of plexiform neurofibroma tumorigenesis [23] (Figure 3D). Our results thus point to a role of paracrine and autocrine signaling defects involving these chemokines, located in 4q12-q13, in the tumorigenesis of plexiform neurofibromas and/or MPNSTs. Comparative genomic hybridization analysis has shown gains of chromosome arm 4q12-ter in both sporadic and NF1-associated MPNSTs [38], suggesting co-overexpression of these chemokines by a DNA amplification mechanism.

We identified two genes (RASSF2, HMMR/RHAMM) putatively involved in the Ras signaling pathway; RASSF2 was down-regulated, whereas HMMR/RHAMM was upregulated. RASSF2 is a member of the RASSF family genes that encode for proteins which bind directly to Ras protein in a GTP-dependent manner via a Ras effector domain. This family currently has six members, RASSF1, RASSF2, RASSF3, RIN2, AD037 and NORE1. RASSF proteins promote both cell cycle arrest and apoptosis; they behave as tumor suppressor genes, and their down-regulation might play a key role in tumorigenesis. RASSF genes are a new class of tumor suppressor genes for which epigenetic silencing is an overwhelming mechanism of inactivation, and somatic mutations are rare [39, 40]. Vos et al. [41]recently obtained evidence that RASSF2 is specific to K-Ras, as only a weak interaction with H-ras was detected. The other five RASSF genes were also tested here, and showed no difference in their expression levels between MPNSTs and plexiform neurofibromas (Figure 3E). The second Ras signaling-associated gene identified here, HMMR/RHAMM, is an oncogene that regulates signaling through Ras and controls mitogen-activated protein kinase (ERK protein) expression [42]. However, HMMR/RHAMM has also recently been characterized as a centrosomal protein that interacts with dynein and maintains spindle pole stability [43]. Finally, HMMR/RHAMM has been identified as a major cell-cycle-regulated gene by cDNA microarray analysis [44].

Interesting, seven of the 28 genes showing significant differential expression between MPNSTs and plexiform neurofibromas are involved in the Hedgehog-Gli signaling pathway, namely DHH, PATCH2 and five downstream target genes of Gli transcription factors (SPP1/Osteopontin, FOXA2/HNF3B, FOXM1, OSF-2/Periostin and EPHA7) [4548]. The Hedgehog-Gli signaling pathway is important in regulating patterning, proliferation, survival and growth in both embryos and adults [49]. Inappropriate activation of the Hedgehog-Gli signaling pathway occurs in several tumor types, including brain and skin tumors [50]. PTCH2 functions as a tumor suppressor gene, as its product normally inhibits SMOH (which functions as an oncogene), resulting in Gli inhibition. Mutations of human PTCH1 and PTCH2 have been detected in basal cell carcinomas and medulloblastomas, and result in Gli signaling activation [50]. Other genes associated with the Hedgehog-Gli signaling pathway were slightly (but not significantly) upregulated in MPNSTs, including SMOH, HHIP (a Hedgehog-interacting protein), three human GLI genes (GLI2, GLI3 and GLI4; Figures 3F) and the major Gli-regulated genes IGF2, PTHR1 and H19 (Figures 3G) [47, 48, 51]. The eight Gli-regulated genes that we found altered in MPNSTs encode growth factors (IGF2, SPP1/Osteopontin, OSF-2/Periostin), growth factor receptors (PTHR1, EPHA7) and transcription factors (FOXA2/HNF3B, FOXM1), whereas H19 yields an untranslated mRNA. Several of these genes are altered in various tumors [45, 5254]. The positive correlation between H19 and IGF2 expression in our tumor series argues against abnormal IGF2 imprinting as a source of IGF2 over-expression in MPNSTs. Indeed, a reduction or loss of H19 expression would be expected to result in increased IGF2 transcription, as expression of these two genes is functionally linked by their competition for a common enhancer [55].

We have previously found that SHH, GLI1, IGF2 and H19 are up-regulated in plexiform neurofibromas [23] (Figures 3F and 3G). Taken together, our results suggest that activation of the Hedgehog-Gli signaling pathway is not only an early event in the genesis of benign plexiform neurofibromas, but that it is also required for malignant transformation.


In conclusion, this study points to the involvement of several altered molecular pathways, and especially the Hedgehog-Gli signaling pathway, in the tumorigenesis of NF1-associated MPNST. Further studies are necessary to elucidate the genetic (or epigenetic) mechanisms responsible for the altered gene expression. It will be of interest to study the genes identified here in sporadic MPNSTs and in other neurological tumors (medulloblastomas, oligodendrogliomas, astrocytomas, neuroblastomas and schwannomas).


Patients and Samples

Samples of 14 plexiform neurofibromas and 9 MPNSTs were obtained by surgical excision from patients with NF1 at Henri Mondor hospital (Creteil, France).

The plexiform neurofibromas (deep lesions involving a plexus of nerves) were large, had a nodular aspect, and severely deformed the affected tissues. They were all S100-positive by immunostaining.

The main clinical and histological characteristics of the 9 patients with MPNSTs are shown in Table 1. The MPNSTs all arose from plexiform neurofibromas and showed very weak S100 immunostaining.

Ten dermal neurofibromas were used as "normal" controls, as they are not at risk of developing into malignant MPNSTs. Neurofibromas are heterogeneous benign tumors composed of Schwann cells, neurons, fibroblasts, mast cells and other cells, and have no "normal" tissue equivalent. The 489 gene expression levels in plexiform neurofibromas and MPNSTs were thus expressed relative to the expression levels in dermal neurofibromas. The dermal neurofibromas were obtained by laser excision from patients free of plexiform neurofibromas. They affected the dermis and subcutis, and were soft, slightly elevated, painless and smaller than 20 mm.

Immediately after surgery the tumor samples were flash-frozen in liquid nitrogen and stored at -80°C until RNA extraction.

Real-time RT-PCR

Theoretical basis

Reactions are characterized by the point during cycling when amplification of the PCR product is first detected, rather than the amount of PCR product accumulated after a fixed number of cycles. The larger the starting quantity of the target molecule, the earlier a significant increase in fluorescence is observed. The parameter Ct (threshold cycle) is defined as the fractional cycle number at which the fluorescence generated by cleavage of a TaqMan probe (or by SYBR green dye-amplicon complex formation) passes a fixed threshold above baseline. The increase in fluorescent signal associated with exponential growth of PCR products is detected by the laser detector of the ABI Prism 7700 Sequence Detection System (Perkin-Elmer Applied Biosystems, Foster City, CA), using PE Biosystems analysis software according to the manufacturer's manuals.

The precise amount of total RNA added to each reaction mix (based on optical density) and its quality (i.e. lack of extensive degradation) are both difficult to assess. We therefore also quantified transcripts of two endogenous RNA control genes involved in two cellular metabolic pathways, namely TBP (Genbank accession NM_003194), which encodes the TATA box-binding protein (a component of the DNA-binding protein complex TFIID), and RPLP0 (also known as 36B4; NM_001002), which encodes human acidic ribosomal phosphoprotein P0. Each sample was normalized on the basis of its TBP (or RPLPO) content.

Results, expressed as N-fold differences in target gene expression relative to the TBP (or RPLPO) gene, and termed "Ntarget", were determined as , where the ΔCt value of the sample was determined by subtracting the average Ct value of the target gene from the average Ct value of the TBP (or RPLP0) gene [56, 57].

The Ntarget values of the samples were subsequently normalized such that the mean of the dermal neurofibroma Ntarget values was 1.

Primers and controls

Primers for TBP, RPLP0 and the 489 target genes (List in annex; supplementary information) were chosen with the assistance of the Oligo 5.0 computer program (National Biosciences, Plymouth, MN).

We conducted searches in dbEST, htgs and nr databases to confirm the total gene specificity of the nucleotide sequences chosen as primers, and the absence of single nucleotide polymorphisms. In particular, the primer pairs were selected to be unique relative to the sequences of closely related family member genes or of the corresponding retropseudogenes. To avoid amplification of contaminating genomic DNA, one of the two primers was placed at the junction between two exons, if possible. In general, amplicons were between 70 and 120 nucleotides long. Gel electrophoresis was used to verify the specificity of PCR amplicons.

For each primer pair, we performed no-template control (NTC) and no-reverse-transcriptase control (RT negative) assays, which produced negligible signals (usually > 40 in Ct value), suggesting that primer-dimer formation and genomic DNA contamination effects were negligible.

RNA extraction

Total RNA was extracted from frozen tumor samples by using the acid-phenol guanidinium method. The quality of the RNA samples was determined by electrophoresis through agarose gels and staining with ethidium bromide, the 18S and 28S RNA bands being visualized under ultraviolet light.

cDNA Synthesis

Total RNA was reverse transcribed in a final volume of 20 μl containing 1X RT buffer (500 μM each dNTP, 3 mM MgCl2, 75 mM KCl, 50 mM Tris-HCl pH 8.3), 20 units of RNasin RNase inhibitor (Promega, Madison, WI), 10 mM DDT, 100 units of Superscript II RNase H-reverse transcriptase (Invitrogen, Cergy Pontoise, France), 3 μM random hexamers (Pharmacia, Uppsala, Sweden) and 100 ng of total RNA. The samples were incubated at 20°C for 10 min and 42°C for 30 min, and reverse transcriptase was inactivated by heating at 99°C for 5 min and cooling at 5°C for 5 min.

PCR amplification

All PCR reactions were performed using an ABI Prism 7700 Sequence Detection System (Perkin-Elmer Applied Biosystems) and either the TaqMan® PCR Core REAGENTS Kit or the SYBR® Green PCR Core Reagents kit (Perkin-Elmer Applied Biosystems). Ten microliters of diluted sample cDNA (produced from 2 ng of total RNA) was added to 15 microliters of the PCR master-mix.

The thermal cycling conditions comprised an initial denaturation step at 95°C for 10 min, and 50 cycles at 95°C for 15 s and 65°C for 1 min.

Statistical Analysis

As the mRNA levels did not fit a Gaussian distribution, (a) the mRNA levels in each subgroup of samples were characterized by their median values and ranges, rather than their mean values and coefficients of variation, and (b) relationships between the molecular markers and clinical and histological parameters were tested using the non parametric Mann-Whitney U test [58]. Differences between two populations were judged significant at confidence levels greater than 95% (p < 0.05).

To visualize the capacity of a given molecular marker to discriminate between two populations (in the absence of an arbitrary cutoff value), we summarized the data in a ROC (receiver operating characteristic) curve [59]. This curve plots the sensitivity (true positives) on the Y axis against 1-specificity (false positives) on the X axis, considering each value as a possible cutoff. The AUC (area under curve) was calculated as a single measure for the discriminatory capacity of each molecular marker. When a molecular marker had no discriminatory value, the ROC curve lies close to the diagonal and the AUC is close to 0.5. In contrast, when a molecular marker has strong discriminatory value, the ROC curve moves to the upper left-hand corner and the AUC is close to 1.0.


  1. 1.

    Friedman JM: Epidemiology of neurofibromatosis type 1. Am J Med Genet. 1999, 89: 1-6. 10.1002/(SICI)1096-8628(19990326)89:1<1::AID-AJMG3>3.3.CO;2-#

  2. 2.

    Cawthon RM, Weiss R, Xu GF, Viskochil D, Culver M, Stevens J, Robertson M, Dunn D, Gesteland R, O'Connell P, : A major segment of the neurofibromatosis type 1 gene: cDNA sequence, genomic structure, and point mutations. Cell. 1990, 62: 193-201. 10.1016/0092-8674(90)90253-B

  3. 3.

    Wallace MR, Marchuk DA, Andersen LB, Letcher R, Odeh HM, Saulino AM, Fountain JW, Brereton A, Nicholson J, Mitchell AL, : Type 1 neurofibromatosis gene: identification of a large transcript disrupted in three NF1 patients. Science. 1990, 249: 181-186.

  4. 4.

    DeClue JE, Papageorge AG, Fletcher JA, Diehl SR, Ratner N, Vass WC, Lowy DR: Abnormal regulation of mammalian p21ras contributes to malignant tumor growth in von Recklinghausen (type 1) neurofibromatosis. Cell. 1992, 69: 265-273. 10.1016/0092-8674(92)90407-4

  5. 5.

    Ferner RE, Gutmann DH: International consensus statement on malignant peripheral nerve sheath tumors in neurofibromatosis. Cancer Res. 2002, 62: 1573-1577.

  6. 6.

    Scheithauer B, Woodruff J, Erlandson R: Tumors of the Peripheral Nervous System. 1999, Washington, DC, Armed Forces Institute of Pathology: Atlas of Tumor Pathology. Third Series, Fascicle 24, 3

  7. 7.

    Woodruff JM: Pathology of tumors of the peripheral nerve sheath in type 1 neurofibromatosis. Am J Med Genet. 1999, 89: 23-30. 10.1002/(SICI)1096-8628(19990326)89:1<23::AID-AJMG6>3.0.CO;2-#

  8. 8.

    Zhu Y, Ghosh P, Charnay P, Burns DK, Parada LF: Neurofibromas in NF1: Schwann cell origin and role of tumor environment. Science. 2002, 296: 920-922. 10.1126/science.1068452

  9. 9.

    Cichowski K, Jacks T: NF1 tumor suppressor gene function: narrowing the GAP. Cell. 2001, 104: 593-604. 10.1016/S0092-8674(01)00245-8

  10. 10.

    Dasgupta B, Gutmann DH: Neurofibromatosis 1: closing the GAP between mice and men. Curr Opin Genet Dev. 2003, 13: 20-27. 10.1016/S0959-437X(02)00015-1

  11. 11.

    Legius E, Marchuk DA, Collins FS, Glover TW: Somatic deletion of the neurofibromatosis type 1 gene in a neurofibrosarcoma supports a tumour suppressor gene hypothesis. Nat Genet. 1993, 3: 122-126.

  12. 12.

    Legius E, Dierick H, Wu R, Hall BK, Marynen P, Cassiman JJ, Glover TW: TP53 mutations are frequent in malignant NF1 tumors. Genes Chromosomes Cancer. 1994, 10: 250-255.

  13. 13.

    Leroy K, Dumas V, Martin-Garcia N, Falzone MC, Voisin MC, Wechsler J, Revuz J, Creange A, Levy E, Lantieri L, Zeller J, Wolkenstein P: Malignant peripheral nerve sheath tumors associated with neurofibromatosis type 1: a clinicopathologic and molecular study of 17 patients. Arch Dermatol. 2001, 137: 908-913.

  14. 14.

    Menon AG, Anderson KM, Riccardi VM, Chung RY, Whaley JM, Yandell DW, Farmer GE, Freiman RN, Lee JK, Li FP, : Chromosome 17p deletions and p53 gene mutations associated with the formation of malignant neurofibrosarcomas in von Recklinghausen neurofibromatosis. Proc Natl Acad Sci U S A. 1990, 87: 5435-5439.

  15. 15.

    Cichowski K, Shih TS, Schmitt E, Santiago S, Reilly K, McLaughlin ME, Bronson RT, Jacks T: Mouse models of tumor development in neurofibromatosis type 1. Science. 1999, 286: 2172-2176. 10.1126/science.286.5447.2172

  16. 16.

    Kourea HP, Orlow I, Scheithauer BW, Cordon-Cardo C, Woodruff JM: Deletions of the INK4A gene occur in malignant peripheral nerve sheath tumors but not in neurofibromas. Am J Pathol. 1999, 155: 1855-1860.

  17. 17.

    Kourea HP, Cordon-Cardo C, Dudas M, Leung D, Woodruff JM: Expression of p27(kip) and other cell cycle regulators in malignant peripheral nerve sheath tumors and neurofibromas: the emerging role of p27(kip) in malignant transformation of neurofibromas. Am J Pathol. 1999, 155: 1885-1891.

  18. 18.

    Nielsen GP, Stemmer-Rachamimov AO, Ino Y, Moller MB, Rosenberg AE, Louis DN: Malignant transformation of neurofibromas in neurofibromatosis 1 is associated with CDKN2A/p16 inactivation. Am J Pathol. 1999, 155: 1879-1884.

  19. 19.

    DeClue JE, Heffelfinger S, Benvenuto G, Ling B, Li S, Rui W, Vass WC, Viskochil D, Ratner N: Epidermal growth factor receptor expression in neurofibromatosis type 1-related tumors and NF1 animal models. J Clin Invest. 2000, 105: 1233-1241.

  20. 20.

    Perrone F, Tabano S, Colombo F, Dagrada G, Birindelli S, Gronchi A, Colecchia M, Pierotti MA, Pilotti S: p15INK4b, p14ARF, and p16INK4a inactivation in sporadic and neurofibromatosis type 1-related malignant peripheral nerve sheath tumors. Clin Cancer Res. 2003, 9: 4132-4138.

  21. 21.

    DeRisi J, Penland L, Brown PO, Bittner ML, Meltzer PS, Ray M, Chen Y, Su YA, Trent JM: Use of a cDNA microarray to analyse gene expression patterns in human cancer. Nat Genet. 1996, 14: 457-460.

  22. 22.

    Paradis V, Bieche I, Dargere D, Laurendeau I, Nectoux J, Degott C, Belghiti J, Vidaud M, Bedossa P: A quantitative gene expression study suggests a role for angiopoietins in focal nodular hyperplasia. Gastroenterology. 2003, 124: 651-659. 10.1053/gast.2003.50104

  23. 23.

    Levy P, Bieche I, Leroy K, Parfait B, Wechsler J, Laurendeau I, Wolkenstein P, Vidaud M, Vidaud D: Molecular profiles of neurofibromatosis type 1-associated plexiform neurofibromas: identification of a gene expression signature of poor prognosis. Clin Cancer Res. 2004, 10: 3763-3771.

  24. 24.

    Latil A, Bieche I, Chene L, Laurendeau I, Berthon P, Cussenot O, Vidaud M: Gene expression profiling in clinically localized prostate cancer: a four-gene expression model predicts clinical behavior. Clin Cancer Res. 2003, 9: 5477-5485.

  25. 25.

    Kindblom LG, Ahlden M, Meis-Kindblom JM, Stenman G: Immunohistochemical and molecular analysis of p53, MDM2, proliferating cell nuclear antigen and Ki67 in benign and malignant peripheral nerve sheath tumours. Virchows Arch. 1995, 427: 19-26.

  26. 26.

    Garratt AN, Britsch S, Birchmeier C: Neuregulin, a factor with many functions in the life of a schwann cell. Bioessays. 2000, 22: 987-996. 10.1002/1521-1878(200011)22:11<987::AID-BIES5>3.3.CO;2-X

  27. 27.

    Feltri ML, Scherer SS, Nemni R, Kamholz J, Vogelbacker H, Scott MO, Canal N, Quaranta V, Wrabetz L: Beta 4 integrin expression in myelinating Schwann cells is polarized, developmentally regulated and axonally dependent. Development. 1994, 120: 1287-1301.

  28. 28.

    Britsch S, Goerich DE, Riethmacher D, Peirano RI, Rossner M, Nave KA, Birchmeier C, Wegner M: The transcription factor Sox10 is a key regulator of peripheral glial development. Genes Dev. 2001, 15: 66-78. 10.1101/gad.186601

  29. 29.

    Neumann AA, Reddel RR: Telomere maintenance and cancer -- look, no telomerase. Nat Rev Cancer. 2002, 2: 879-884. 10.1038/nrc929

  30. 30.

    Altieri DC: Validating survivin as a cancer therapeutic target. Nat Rev Cancer. 2003, 3: 46-54. 10.1038/nrc968

  31. 31.

    Wurl P, Kappler M, Meye A, Bartel F, Kohler T, Lautenschlager C, Bache M, Schmidt H, Taubert H: Co-expression of survivin and TERT and risk of tumour-related death in patients with soft-tissue sarcoma. Lancet. 2002, 359: 943-945. 10.1016/S0140-6736(02)07990-4

  32. 32.

    Melino G, De Laurenzi V, Vousden KH: p73: Friend or foe in tumorigenesis. Nat Rev Cancer. 2002, 2: 605-615. 10.1038/nrc861

  33. 33.

    Stamenkovic I: Matrix metalloproteinases in tumor invasion and metastasis. Semin Cancer Biol. 2000, 10: 415-433. 10.1006/scbi.2000.0379

  34. 34.

    Jiang Y, Goldberg ID, Shi YE: Complex roles of tissue inhibitors of metalloproteinases in cancer. Oncogene. 2002, 21: 2245-2252. 10.1038/sj.onc.1205291

  35. 35.

    Chang MS, McNinch J, Basu R, Simonet S: Cloning and characterization of the human neutrophil-activating peptide (ENA-78) gene. J Biol Chem. 1994, 269: 25277-25282.

  36. 36.

    Fujioka T, Kolson DL, Rostami AM: Chemokines and peripheral nerve demyelination. J Neurovirol. 1999, 5: 27-31.

  37. 37.

    O'Donovan N, Galvin M, Morgan JG: Physical mapping of the CXC chemokine locus on human chromosome 4. Cytogenet Cell Genet. 1999, 84: 39-42.

  38. 38.

    Koga T, Iwasaki H, Ishiguro M, Matsuzaki A, Kikuchi M: Frequent genomic imbalances in chromosomes 17, 19, and 22q in peripheral nerve sheath tumours detected by comparative genomic hybridization analysis. J Pathol. 2002, 197: 98-107. 10.1002/path.1101

  39. 39.

    Hesson L, Dallol A, Minna JD, Maher ER, Latif F: NORE1A, a homologue of RASSF1A tumour suppressor gene is inactivated in human cancers. Oncogene. 2003, 22: 947-954. 10.1038/sj.onc.1206191

  40. 40.

    Dammann R, Li C, Yoon JH, Chin PL, Bates S, Pfeifer GP: Epigenetic inactivation of a RAS association domain family protein from the lung tumour suppressor locus 3p21.3. Nat Genet. 2000, 25: 315-319. 10.1038/77083

  41. 41.

    Vos MD, Ellis CA, Elam C, Ulku AS, Taylor BJ, Clark GJ: RASSF2 is a novel K-Ras-specific effector and potential tumor suppressor. J Biol Chem. 2003, 278: 28045-28051. 10.1074/jbc.M300554200

  42. 42.

    Zhang S, Chang MC, Zylka D, Turley S, Harrison R, Turley EA: The hyaluronan receptor RHAMM regulates extracellular-regulated kinase. J Biol Chem. 1998, 273: 11342-11348. 10.1074/jbc.273.18.11342

  43. 43.

    Maxwell CA, Keats JJ, Crainie M, Sun X, Yen T, Shibuya E, Hendzel M, Chan G, Pilarski LM: RHAMM is a centrosomal protein that interacts with dynein and maintains spindle pole stability. Mol Biol Cell. 2003, 14: 2262-2276. 10.1091/mbc.E02-07-0377

  44. 44.

    Cho RJ, Huang M, Campbell MJ, Dong H, Steinmetz L, Sapinoso L, Hampton G, Elledge SJ, Davis RW, Lockhart DJ: Transcriptional regulation and function during the human cell cycle. Nat Genet. 2001, 27: 48-54.

  45. 45.

    Teh MT, Wong ST, Neill GW, Ghali LR, Philpott MP, Quinn AG: FOXM1 is a downstream target of Gli1 in basal cell carcinomas. Cancer Res. 2002, 62: 4773-4780.

  46. 46.

    Sasaki H, Hui C, Nakafuku M, Kondoh H: A binding site for Gli proteins is essential for HNF-3beta floor plate enhancer activity in transgenics and can respond to Shh in vitro. Development. 1997, 124: 1313-1322.

  47. 47.

    Yoon JW, Kita Y, Frank DJ, Majewski RR, Konicek BA, Nobrega MA, Jacob H, Walterhouse D, Iannaccone P: Gene expression profiling leads to identification of GLI1-binding elements in target genes and a role for multiple downstream pathways in GLI1-induced cell transformation. J Biol Chem. 2002, 277: 5548-5555. 10.1074/jbc.M105708200

  48. 48.

    Louro ID, Bailey EC, Li X, South LS, McKie-Bell PR, Yoder BK, Huang CC, Johnson MR, Hill AE, Johnson RL, Ruppert JM: Comparative gene expression profile analysis of GLI and c-MYC in an epithelial model of malignant transformation. Cancer Res. 2002, 62: 5867-5873.

  49. 49.

    Ruiz i Altaba A, Sanchez P, Dahmane N: Gli and hedgehog in cancer: tumours, embryos and stem cells. Nat Rev Cancer. 2002, 2: 361-372. 10.1038/nrc796

  50. 50.

    Villavicencio EH, Walterhouse DO, Iannaccone PM: The sonic hedgehog-patched-gli pathway in human development and disease. Am J Hum Genet. 2000, 67: 1047-1054.

  51. 51.

    Hahn H, Wojnowski L, Specht K, Kappler R, Calzada-Wack J, Potter D, Zimmer A, Muller U, Samson E, Quintanilla-Martinez L: Patched target Igf2 is indispensable for the formation of medulloblastoma and rhabdomyosarcoma. J Biol Chem. 2000, 275: 28341-28344. 10.1074/jbc.C000352200

  52. 52.

    Evliyaoglu C, Carroll R, Folkerth R, Bello L, Bruns DE, Black PM: Parathyroid hormone-related protein and its receptor in human glial tumors. Acta Neurochir (Wien). 2000, 142: 871-878. 10.1007/s007010070072

  53. 53.

    Saitoh Y, Kuratsu J, Takeshima H, Yamamoto S, Ushio Y: Expression of osteopontin in human glioma. Its correlation with the malignancy. Lab Invest. 1995, 72: 55-63.

  54. 54.

    Watson MA, Gutmann DH, Peterson K, Chicoine MR, Kleinschmidt-DeMasters BK, Brown HG, Perry A: Molecular characterization of human meningiomas by gene expression profiling using high-density oligonucleotide microarrays. Am J Pathol. 2002, 161: 665-672.

  55. 55.

    Joyce JA, Lam WK, Catchpoole DJ, Jenks P, Reik W, Maher ER, Schofield PN: Imprinting of IGF2 and H19: lack of reciprocity in sporadic Beckwith-Wiedemann syndrome. Hum Mol Genet. 1997, 6: 1543-1548. 10.1093/hmg/6.9.1543

  56. 56.

    Bieche I, Parfait B, Le Doussal V, Olivi M, Rio MC, Lidereau R, Vidaud M: Identification of CGA as a novel estrogen receptor-responsive gene in breast cancer: an outstanding candidate marker to predict the response to endocrine therapy. Cancer Res. 2001, 61: 1652-1658.

  57. 57.

    Bieche I, Onody P, Laurendeau I, Olivi M, Vidaud D, Lidereau R, Vidaud M: Real-time reverse transcription-PCR assay for future management of ERBB2-based clinical applications. Clin Chem. 1999, 45: 1148-1156.

  58. 58.

    Mann H, Whitney D: On a test of whether one of two random variables is stochastically larger than the other. Annals of Mathematical Statistics. 1947, 18: 50-60.

  59. 59.

    Hanley JA, McNeil BJ: The meaning and use of the area under a receiver operating characteristic (ROC) curve. Radiology. 1982, 143: 29-36.

Download references


We thank Juliette Nectoux and Eric Pasmant for the design of oligonucleotide primers. This work was supported by Association pour la Recherche sur le Cancer, Association Neurofibromatoses et Recklinghausen and Ministère de l'Enseignement Supérieur et de la Recherche.

Author information

Correspondence to Ivan Bièche.

Additional information

Authors' contributions

PW was responsible for patients screening. KL and JW determined the histological diagnosis. PL, IL, BP and GB designed the 489 primers. Real-time RT-PCR have been carried out by IL. IB, IL and PL interpreted the result, performed bioinformatics and statistical analyses. IB managed the limiting factors. This study has been supervised by DV, MV and IB.

Electronic supplementary material

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Lévy, P., Vidaud, D., Leroy, K. et al. Molecular profiling of malignant peripheral nerve sheath tumors associated with neurofibromatosis type 1, based on large-scale real-time RT-PCR. Mol Cancer 3, 20 (2004).

Download citation


  • Schwann Cell
  • Neurofibroma
  • Malignant Peripheral Nerve Sheath Tumor
  • Extracellular Matrix Remodel
  • Plexiform Neurofibroma