- Open Access
DNA copy number changes in high-grade malignant peripheral nerve sheath tumors by array CGH
Molecular Cancervolume 7, Article number: 48 (2008)
Malignant peripheral nerve sheath tumors (MPNSTs) are rare and highly aggressive soft tissue tumors showing complex chromosomal aberrations. In order to identify recurrent chromosomal regions of gain and loss, and thereby novel gene targets of potential importance for MPNST development and/or progression, we have analyzed DNA copy number changes in seven high-grade MPNSTs using microarray-based comparative genomic hybridization (array CGH).
Considerable more gains than losses were observed, and the most frequent minimal recurrent regions of gain included 1q24.1-q24.2, 1q24.3-q25.1, 8p23.1-p12, 9q34.11-q34.13 and 17q23.2-q25.3, all gained in five of seven samples. The 17q23.2-q25.3 region was gained in all five patients with poor outcome and not in the two patients with disease-free survival. cDNA microarray analysis and quantitative real-time reverse transcription PCR were used to investigate expression of genes located within these regions. The gene lysyl oxidase-like 2 (LOXL2) was identified as a candidate target for the 8p23.1-p12 gain. Within 17q, the genes topoisomerase II-α (TOP2A), ets variant gene 4 (E1A enhancer binding protein, E1AF) (ETV4) and baculoviral IAP repeat-containing 5 (survivin) (BIRC5) showed increased expression in all samples compared to two benign tumors. Increased expression of these genes has previously been associated with poor survival in other malignancies, and for TOP2A, in MPNSTs as well. In addition, we have analyzed the expression of five micro RNAs located within the 17q23.2-q25.3 region, but none of them showed high expression levels compared to the benign tumors.
Our study shows the potential of using DNA copy number changes obtained by array CGH to predict the prognosis of MPNST patients. Although no clear correlations between the expression level and patient outcome were observed, the genes TOP2A, ETV4 and BIRC5 are interesting candidate targets for the 17q gain associated with poor survival.
Malignant peripheral nerve sheath tumors (MPNSTs) are rare tumors that arise sporadically or as part of the neurofibromatosis type 1 (NF1) or -2 (NF2) autosomal inherited disorder. The NF1/von Recklinghausen neurofibromatosis, caused by germ line mutations of the NF1 tumor suppressor gene, is one of the most common autosomal dominant inherited disorders, occurring at a frequency of one in every 4,000 individuals . Patients with this disease have an increased risk of benign and malignant tumors . In contrast to other soft tissue malignancies, the majority of MPNSTs derives from previously existing neurofibromas .
Cytogenetically, MPNSTs have complex karyotypes with multiple losses frequently observed in chromosome regions 1p, 9p, 17q and 22 [4, 5]. Several comparative genomic hybridization (CGH) studies have revealed a higher frequency of gains compared to losses, involving chromosome regions 5p, 7p, 7q, 8q and 17q [6–10]. Recurrent gain of 7p15-p21 and 17q22-qter has been associated with poor overall survival , and increased copy number and expression of topoisomerase II-α (TOP2A) in 17q21.2 have been associated with poor cancer-specific survival and presence of metastasis .
In order to identify specific genomic events and candidate targets that may play a role in MPNST development and/or progression, we have used microarray-based CGH (array CGH) to map the distribution and frequency of DNA copy number changes in seven high-grade MPNSTs. cDNA microarray analysis and quantitative real-time reverse transcription PCR (RT-PCR) were used to investigate expression of genes located within the most frequently altered chromosomal regions. In addition, micro RNA (miRNA) expression in a recurrent region of gain was determined using quantitative real-time RT-PCR.
Recurrently altered chromosomal regions in MPNSTs
DNA copy number changes in seven high-grade MPNSTs (Table 1) were analyzed using a 1 Mb resolution bacterial- and P1 artificial chromosome (BAC and PAC) genomic microarray supplemented with the tiling-path between 1q12 and the beginning of 1q25. A heat map of DNA copy number ratios of the tumor samples is shown in Figure 1.
Regions with significant DNA copy number changes in each sample were identified using the "Analysis of Copy Errors" (ACE) algorithm in CGH-Explorer. The resulting frequency plot of gains and losses is shown in Figure 2a, and a representative ratio plot for this type of tumors in Figure 2b. Genome-wide ratio plots for all samples are shown in Additional file 1. Minimal recurrent regions of alteration identified by ACE in at least three of seven (≥ 43%) samples are presented in Table 2. The complete list of data of all defined regions of gain and loss from the ACE analysis is presented in Additional file 2.
The tumors showed considerably more recurrent gains than losses. Thirty-five of the 40 identified recurrent regions of alteration were gains, compared with only five regions of loss. The most frequent regions of increased copy number were in 1q, 8p, 9q and 17q, all detected in five of seven tumors. In 1q, three minimal recurrent regions of gain were identified; 1q24.1-q24.2 (3.7 Mb) and 1q24.3-q25.1 (3.5 Mb) in five of seven tumors and 1q21.1 (1.2 Mb) in three tumors. The other regions gained in five tumors were 8p23.1-p12 (22.9 Mb), 9q34.11-q34.13 (1.0 Mb) and 17q23.2-q25.3 (26.2 Mb). Ten minimal recurrent regions of gain were observed in four tumors; 2q11.2-q13 (11.0 Mb), 3p26.2-p25.1 (12.6 Mb), 5q34-q35.3 (15.2 Mb), 7q11.23-q21.11 (7.4 Mb), 9q21.32-q22.33 (16 Mb), 12q13.3-q15 (14.2 Mb), 13q22.1-q22.2 (2.1 Mb), 16p13.3-p13.2 (7.9 Mb), 16p13.12-p13.11 (0.9 Mb) and 19p13.3-p13.2 (4.9 Mb). In addition, 20 regions of gain were identified in three tumors (see Table 2). High-level amplification (log2 ratio > 1) was observed in some of the tumors, mainly of regions in 1q and 12q [see Additional file 2].
Three of five identified minimal recurrent regions of loss were located in chromosome 11. Four tumors showed loss of 11p13 (1.5 Mb), whereas 11q22.3-q23.1 (7.7 Mb) and 11q23.2-q23.3 (7.0 Mb) were lost in three tumors. In addition, loss of 9p22.3-p21.2 (12.8 Mb) and 14q21.3-q23.3 (15.3 Mb) were observed in three tumors. Homozygous deletion (log2 ratio < -1) was observed in some of the tumors, mainly of regions in 9p [see Additional file 2].
Gene expression in frequently altered chromosomal regions
Gene expression has previously been analyzed using cDNA microarrays in a panel of soft tissue sarcomas, including six of the MPNSTs studied here (using the xenograft of MPNST2) . In order to identify candidate target genes for the DNA copy number changes, the expression level of genes located within the most frequently altered regions (gained in five of seven samples) was investigated. Genes with increased expression relative to the median for soft tissue sarcomas (log2 ratio > 1) in two or more of the six MPNSTs analyzed were identified.
Within the minimal recurrent region of gain in 8p23.1-p12, six genes showed increased expression . Two genes were over-expressed in three samples; lysyl oxidase-like 2 (LOXL2) and zinc finger protein 395 (ZNF395). In two of the samples, the genes mitochondrial tumor suppressor 1 (MTUS1), leucine zipper, putative tumor suppressor 1 (LZTS1), scavenger receptor class A, member 3 (SCARA3) and UBX domain-containing protein 6 (Reproduction 8 protein) (Protein Rep-8) (UBXD6) were over-expressed (data not shown).
The expression level of LOXL2 and ZNF395 was in addition determined using quantitative real-time RT-PCR in six of the samples (using the xenograft of MPNST2 and -13 and the patient sample of MPNST8x). The expression level was compared to the average expression in two benign tumors, one benign schwannoma (BS) and one neurofibroma (NF) (Table 1). Figure 3 shows the relative expression level of LOXL2 and ZNF395. LOXL2 showed increased expression in MPNST1, -7 and -9x, whereas ZNF395 showed approximately similar expression levels in the MPNSTs and the benign tumors.
In the minimal recurrent regions of gain in 1q24.1-q24.2, 1q24.3-q25.1 and 9q34.11-q34.13, no genes present in the cDNA microarray showed increased expression in more than one of the samples.
Homozygous deletion (log2 ratio < -1) of the region containing the tumor suppressor gene cyclin-dependent kinase inhibitor 2A (melanoma, p16, inhibits CDK4) (CDKN2A) was observed in two samples, and heterozygous deletion in one sample [see Additional file 2]. The expression level of CDKN2A was determined using quantitative real-time RT-PCR in six of the samples (Figure 3). The three samples with deletion of the region, MPNST2x, -7 and -9x, showed no or low expression compared to the benign tumors, whereas the other three MPNSTs showed increased expression.
Characterization of the 17q gain
The minimal recurrent region 17q23.2-q25.3 was gained in all patients who died of the disease (five of seven patients), but not in the two patients with disease-free survival (MPNST1 and -13, with follow-up of 100 and 216 months, respectively). Figure 4a shows the copy number of chromosome 17 for all the tumor samples.
The expression level of genes located within the minimal recurrent region, as well as the proximal region where only four of seven samples showed gain (17q21.2-q23.2), was examined. By cDNA microarray analysis, ten genes showed increased expression relative to the median for soft tissue sarcomas (log2 ratio > 1) in two or more of the six MPNSTs analyzed . Three genes were over-expressed in three MPNSTs; baculoviral IAP repeat-containing 5 (survivin) (BIRC5), ets variant gene 4 (E1A enhancer binding protein, E1AF) (ETV4) and homeobox B7 (HOXB7). In two of the samples, the genes RAB5C, member RAS oncogene family (RAB5C), vesicle amine transport protein 1 homolog (T. californica) (VAT1), distal-less homeobox 4 (DLX4), DEAD (Asp-Glu-Ala-Asp) box polypeptide 5 (DDX5), thymidine kinase 1, soluble (TK1), stimulated by retinoic acid 13 homolog (mouse) (STRA13) and solute carrier family 16, member 3 (monocarboxylic acid transporter 4) (SLC16A3) were over-expressed (data not shown).
The expression level of ETV4, HOXB7 and BIRC5 was in addition determined using quantitative real-time RT-PCR in six of the samples. The expression level of TOP2A was also determined, since increased expression of TOP2A has been associated with poor cancer-specific survival in MPNSTs . Figure 4b shows the relative expression level of TOP2A, ETV4, HOXB7 and BIRC5. Of these four genes, only BIRC5 is located within the minimal recurrent region of gain. TOP2A, ETV4 and BIRC5 showed increased expression in all MPNST samples compared to the benign tumors, in most cases more than 20-fold, whereas HOXB7 showed approximately similar expression levels in the MPNSTs and the benign tumors. ETV4 showed the highest level of expression; its expression was more than 140-fold higher in MPNST9x and more than 90-fold higher in MPNST1 compared to the benign tumors. No clear differences in the expression levels were seen between MPNSTs from patients with poor outcome and MPNSTs from patients who showed disease-free survival (MPNST1 and -13).
The expression levels of miRNAs located within 17q were in addition investigated, in order to see if their expression was increased due to the genomic gain. At the time the analyses were performed, five miRNAs were identified within the minimal recurrent region of gain according to Ensembl; hsa-miR-142p-3p, -142-5p, -301, -21 and -338. The expression level of these miRNAs was determined using quantitative real-time RT-PCR in six of the samples. Figure 4c shows the relative expression level of hsa-miR-142p-3p, -142-5p, -301, -21 and -338. The two miRNAs hsa-miR-142-3p and -142-5p, originating from the different arms of the stem-loop precursor, showed a similar expression pattern. In general, none of the miRNAs showed high expression levels in the MPNSTs compared to the benign tumors, and no clear differences in the miRNA expression levels were seen between MPNSTs from patients with poor outcome and MPNSTs from patients who showed disease-free survival (MPNST1 and -13).
We have used array CGH to analyze DNA copy number changes in a small panel of high-grade MPNSTs, in order to identify recurrent copy number alterations at high-resolution and thereby novel candidate oncogenes and/or tumor suppressor genes. After a recent pathological revision of our tumor collection, seven samples were classified as MPNSTs based on the current classification standard and included in the study. Although this is a small number of samples, the analyses revealed results of potential importance for this malignancy.
Considerably more gains than losses were observed in the seven MPNSTs. Only five of the 40 identified minimal recurrent regions of alteration were losses (see Figure 2a and Table 2), similar to previous observations by others [7, 10, 13]. However, other studies have reported far more frequent losses than gains in MPNSTs as well [4, 14]. Differences in the chromosomal regions altered in sporadic versus neurofibromatosis-associated MPNSTs have previously been demonstrated [6, 9, 10, 13], but the only difference observed here was loss of 11q23.2-q23.3. This region was deleted in the three sporadic MPNSTs, but not in the four neurofibromatosis-associated MPNSTs. However, two of the neurofibromatosis-associated MPNSTs showed loss of other parts of chromosome 11, and the two other minimal recurrent regions of loss in chromosome 11, 11p13 and 11q22.3-q23.1, were deleted in both sporadic and neurofibromatosis-associated MPNSTs [see Additional file 2].
It has previously been reported that patients with neurofibromatosis-associated MPNSTs have a worse survival than patients with sporadic MPNSTs [15, 16], but this has not been a consistent finding [17, 18]. In this panel of tumors, the two longest surviving patients who showed no evidence of the MPNST tumor after 100 and 216 months, respectively, have neurofibromatosis (Table 1).
The most frequent alterations observed in the seven MPNSTs were gains of regions in 1q, 8p, 9q and 17q, all present in five of the seven tumors analyzed. Within 1q, three minimal recurrent regions were identified. 1q24.1-q24.2 (3.7 Mb) and 1q24.3-q25.1 (3.5 Mb) were gained in five of seven samples, whereas 1q21.1 (1.2 Mb) was gained in three samples. Gain of regions in 1q has not frequently been reported in MPNSTs previously, but it is a recurrent finding in soft tissue sarcomas [19, 20] and other malignancies. Since gene expression of six of the MPNSTs studied here has previously been analyzed using cDNA microarrays , the expression level of genes located within these regions was investigated, but none of the genes present on the cDNA microarray were over-expressed compared to the median for soft tissue sarcomas (log2 ratio > 1) in more than one sample. This was also the case for genes located within the minimal recurrent region in 9q34.11-q34.13 (1.0 Mb).
A minimal recurrent region in 8p23.1-p12 (22.9 Mb) was gained in five of the samples. Within this region, six genes showed increased expression in two or more of the samples . LOXL2 and ZNF395 were over-expressed in three of the samples. The expression level of these two genes was in addition determined using quantitative real-time RT-PCR, and the expression level was compared to two benign tumors (Figure 3). LOXL2 showed more than 3-fold increased expression in three MPNSTs, whereas ZNF395 showed approximately similar expression levels in the MPNSTs and the benign tumors. Increased expression of LOXL2, a member of the lysyl oxidase family, has previously been shown in colon- and esophageal cancer , and it has also been associated with breast cancer tumor grade . Thus, LOXL2 may be a candidate target for the 8p23.1-p12 gain in MPNSTs.
Scattered high-level amplification (log2 ratio > 1) and homozygous deletion (log2 ratio < -1) were observed in some of the tumors [see Additional file 2]. Within 9p22.3-p21.2, homozygous deletion of the region harboring the gene CDKN2A was observed in two samples, and heterozygous deletion in one sample [see Additional file 2]. Inactivation of this region is a frequent finding in MPNSTs [23, 24], as well as other cancer types . The expression level of CDKN2A was determined using quantitative real-time RT-PCR (Figure 3), and the three samples with deletion of the region showed no or very low expression compared to the benign tumors, whereas the other three samples showed increased expression. Notably, the expression level of CDKN2A was lower in the four patients with poor outcome than the two patients with disease-free survival (Figure 3).
We observed that all five patients with poor outcome showed gain of the distal part of 17q, whereas the two patients with disease-free survival did not (Figure 4a). Several other studies have also reported that gain of 17q is frequent in MPNSTs [6, 7, 9, 10, 13, 26], and this alteration has been associated with poor outcome  and development of metastasis . These CGH studies have reported the region of gain to be from either 17q22 or 17q24 to 17q25/17qtel. The minimal recurrent region identified in this study, using microarrays, covered 17q23.2-q25.3 (qtel). This region was gained in five tumors, and four of these showed also gain of the proximal region (17q21.2-q23.2).
The expression level of genes located within the minimal recurrent region, as well as the proximal region where four of seven samples showed gain, was investigated. Ten genes showed increased expression relative to the median for soft tissue sarcomas in two or more of the six tumors analyzed , although it should be noted that the cDNA microarray only contained probes for about half of the genes in the region (including novel and predicted genes). Three genes were over-expressed in three of the tumors; ETV4, HOXB7 and BIRC5, and the expression levels of these genes, as well as TOP2A, were in addition determined using quantitative real-time RT-PCR (Figure 4b). Of these four genes, only BIRC5 is located within the minimal recurrent region of gain. TOP2A, ETV4 and BIRC5 showed increased expression in all MPNST samples, in most cases more than 20-fold, whereas HOXB7 showed approximately similar expression levels in the MPNSTs and the benign tumors. Even though only tumors from patients who died of the disease showed gain of this region, increased expression of these genes was also seen in tumors from the patients who showed no evidence of the disease. Thus, there were no correlations between the expression level of the genes and patient survival.
BIRC5, also known as survivin, is an inhibitor of apoptosis that has been shown to be highly expressed in the majority of cancers, including soft tissue sarcomas and osteosarcomas [27, 28]. Increased expression of BIRC5 has been associated with chemotherapy resistance, enhanced proliferation, increased tumor recurrence and shorter patient survival . BIRC5 is located within a 2 Mb region in 17q25 previously shown to be commonly amplified in MPNSTs , and increased expression of BIRC5 in MPNSTs compared to neurofibromas and benign schwannomas has been reported [29, 30]. Hence, there is considerable evidence suggesting that BIRC5 may be involved in MPNST tumorigenesis.
Increased expression of TOP2A has been previously associated with poor cancer-specific survival in MPNSTs , whereas increased expression of ETV4, a member of the Ets family of transcription factors, has been associated with shorter patient survival in colorectal cancer  and gastric cancer . Although these two genes were not located within the minimal recurrent region of gain, their associations with poor outcome in MPNSTs (for TOP2A) and other malignancies suggest that increased expression of these genes may play a role in MPNST tumorigenesis as well.
Recently, it has been shown that expression of miRNAs can be deregulated in cancer, and that miRNAs may act as oncogenes and tumor suppressor genes . miRNAs are frequently found in regions of DNA copy number aberrations, and a general correlation between miRNA copy number and expression level has been reported [34, 35]. In order to investigate if miRNAs may be candidate targets for the 17q gain, we analyzed the expression level of five miRNAs present in the minimal recurrent region (Figure 4c). Compared to the benign tumors, none of the miRNAs showed high expression levels, and no clear differences in the miRNA expression levels were seen between MPNSTs from patients with poor outcome and MPNSTs from patients who showed disease-free survival. Thus, no clear correlation could be found between DNA copy number and miRNA expression in 17q in our samples.
Our study has identified recurrent copy number alterations in MPNSTs at high resolution, and shows the potential of using DNA copy number changes obtained by array CGH to predict the prognosis of MPNST patients. LOXL2 was identified as a candidate target gene for the 8p23.1-p12 gain. Although no clear correlations between the expression level and patient outcome were observed here, the genes TOP2A, ETV4 and BIRC5 are interesting candidate targets for the 17q gain associated with poor outcome, but further validation is required on a larger tumor set.
Six human sarcomas classified as MPNSTs were selected from a tumor collection at the Department of Tumor Biology at the Norwegian Radium Hospital. One additional sample initially diagnosed as malignant fibrous histiocytoma was included in the study after reclassification to MPNST. Two benign tumors, one BS and one NF, were used as a reference for the quantitative real-time RT-PCR analyses. All tumors were revised at the time of the study by the pathologist (B.B.) and diagnosed according to the current World Health Organization classification. The informed consent used and the collection of samples were approved by the ethical committee of Southern Norway.
Clinical samples were collected immediately after surgery, cut into small pieces, frozen in liquid nitrogen and stored at -70°C until use. Some of the samples were grown subcutaneously in immunodeficient mice as xenografts (suffix x). Animal care was in accordance with the institution's guidelines. Clinical data for all samples are given in Table 1.
The genomic microarray used contained 4,549 BAC and PAC clones representing the human genome at approximately 1 Mb resolution, as well as the minimal tiling-path between 1q12 and the beginning of 1q25. Detailed information on the construction and preparation of the microarray has been previously described . The microarrays were provided by the Norwegian Microarray Consortium.
Array CGH was performed essentially as described previously . In brief, approximately 500 ng of Dpn II-digested total genomic DNA was labeled by random priming using BioPrime DNA Labeling System (Invitrogen, California, USA) and Cy3-dCTP (tumor) or Cy5-dCTP (reference) (PerkinElmer, Massachusetts, USA). Labeled tumor and reference DNA were combined together with 135 μg human Cot-1 DNA (Invitrogen). Hybridization was performed using an automated hybridization station, GeneTAC (Genomic Solutions/PerkinElmer), agitating the hybridization solution for 42–46 hours at 37°C. The arrays were scanned using an Agilent G2565BA scanner (Agilent Technologies, California, USA), and the images were segmented using GenePix Pro 6.0 (Axon Laboratories, California, USA). Further data processing, including filtering and normalization, was performed using M-CGH as previously described [36, 37].
Array CGH data analysis
The complete array CGH dataset for the seven MPNSTs can be viewed in the ArrayExpress microarray database (accession number E-MEXP-869). Clones belonging to chromosomes 1–22 with known unique chromosomal location in Ensembl (v33, Sep 2005) were considered for analysis (3,351 clones). Due to experimental variation in normal control experiments, 22 clones (0.7%) were discarded as described previously . In addition, clones with missing values in three or more of the seven samples were discarded, leaving 3,167 clones for analysis. The remaining missing values were imputed via a K-Nearest Neighbor algorithm normalization using "Significance Analysis of Microarrays" .
In order to determine copy number changes, CGH-Explorer v. 2.55 was used . ACE was performed using a false discovery rate of 0.0000. Chromosomal segments showing gains or losses in at least three of seven MPNSTs (≥ 43%) were used to identify minimal recurrent regions of alteration.
Quantitative real-time RT-PCR
Quantitative real-time RT-PCR was performed using TaqMan Gene Expression and MicroRNA Assays (Applied Biosystems, California, USA). The expression level was determined for the genes LOXL2 (assay ID Hs00158757_m1), ZNF395 (assay ID Hs00608626_m1), CDKN2A (assay ID Hs00924091_m1), TOP2A (assay ID Hs01032127_g1), ETV4 (assay ID Hs00385910_m1), HOXB7 (assay ID Hs00270131_m1) and BIRC5 (assay ID Hs00153353_m1). The genes beta-2-microglobulin (B2M, assay ID Hs99999907_m1), glyceraldehyde-3-phosphate dehydrogenase (GAPDH, assay ID Hs99999905_m1) and TATA-box binding protein (TBP, assay ID Hs99999910_m1) were used as endogenous controls for normalization of gene expression. The expression level was determined for the miRNAs hsa-miR-142-3p (assay ID 464), hsa-miR-142-5p (assay ID 465), hsa-miR-301 (assay ID 528), hsa-miR-21 (assay ID 397) and hsa-miR-338 (assay ID 548). The small nuclear RNAs RNU6B (assay ID 1093) and RNU24 (assay ID 1001) were used as endogenous controls for normalization of miRNA expression.
Frozen tumor tissue was pulverized in liquid nitrogen, and total RNA was extracted using Trizol (Invitrogen) according to the manufacturer's instructions. The total RNA was further purified using the RNeasy Mini Kit (QIAGEN, California, USA) as described by the manufacturer, with a few modifications in order to preserve the miRNAs. In brief, after addition of buffer RLT and vortexing, 3.5 volumes of 100% ethanol was added and mixed by vortexing. The sample was subsequently applied to the RNeasy Mini column. Washing with buffer RW1 was not performed. Universal Human Reference RNA (Stratagene, California, USA) was used as a reference for the gene expression assays. cDNA synthesis and real-time PCR were performed essentially as described in the protocols supplied by the manufacturer (Applied Biosystems). The PCR amplification was performed using the ABI 7500 Real Time PCR System (Applied Biosystems).
Each gene expression assay included (in duplicate) a standard curve of four serial dilutions of the Universal Human Reference RNA cDNA (ranging from 50 ng to 50 pg), 10 ng of each tumor cDNA and a no-template control. The expression levels were determined from the standard curves as described by the manufacturer. The expression level of each gene was normalized with the average expression of the three endogenous controls. The expression level of each gene in the MPNSTs was determined relative to the average expression of the benign tumors.
Each miRNA expression assay included (in duplicate) 0.3 ng of each tumor cDNA and a no-template control. The experiments were done twice, and the average values were used. The expression levels were determined using the comparative CT method as described by the manufacturer, and the expression level of the miRNAs was normalized with the average expression of the two endogenous controls. The expression level of each miRNA in the MPNSTs was determined relative to the average expression of the benign tumors.
Friedman JM: Epidemiology of neurofibromatosis type 1. Am J Med Genet. 1999, 89 (1): 1-6. 10.1002/(SICI)1096-8628(19990326)89:1<1::AID-AJMG3>3.0.CO;2-8
Sorensen SA, Mulvihill JJ, Nielsen A: Long-term follow-up of von Recklinghausen neurofibromatosis. Survival and malignant neoplasms. N Engl J Med. 1986, 314 (16): 1010-1015.
Woodruff JM: Pathology of the major peripheral nerve sheath neoplasms. Monogr Pathol. 1996, 38: 129-161.
Mertens F, Dal Cin P, De Wever I, Fletcher CD, Mandahl N, Mitelman F, Rosai J, Rydholm A, Sciot R, Tallini G, van Den Berghe H, Vanni R, Willen H: Cytogenetic characterization of peripheral nerve sheath tumours: a report of the CHAMP study group. J Pathol. 2000, 190 (1): 31-38. 10.1002/(SICI)1096-9896(200001)190:1<31::AID-PATH505>3.0.CO;2-#
Mertens F, Rydholm A, Bauer HF, Limon J, Nedoszytko B, Szadowska A, Willen H, Heim S, Mitelman F, Mandahl N: Cytogenetic findings in malignant peripheral nerve sheath tumors. Int J Cancer. 1995, 61 (6): 793-798. 10.1002/ijc.2910610609
Lothe RA, Karhu R, Mandahl N, Mertens F, Saeter G, Heim S, Borresen-Dale AL, Kallioniemi OP: Gain of 17q24-qter detected by comparative genomic hybridization in malignant tumors from patients with von Recklinghausen's neurofibromatosis. Cancer Res. 1996, 56 (20): 4778-4781.
Mechtersheimer G, Otano-Joos M, Ohl S, Benner A, Lehnert T, Willeke F, Moller P, Otto HF, Lichter P, Joos S: Analysis of chromosomal imbalances in sporadic and NF1-associated peripheral nerve sheath tumors by comparative genomic hybridization. Genes Chromosomes Cancer. 1999, 25 (4): 362-369. 10.1002/(SICI)1098-2264(199908)25:4<362::AID-GCC8>3.0.CO;2-Q
Nakagawa Y, Yoshida A, Numoto K, Kunisada T, Wai D, Ohata N, Takeda K, Kawai A, Ozaki T: Chromosomal imbalances in malignant peripheral nerve sheath tumor detected by metaphase and microarray comparative genomic hybridization. Oncol Rep. 2006, 15 (2): 297-303.
Schmidt H, Taubert H, Meye A, Wurl P, Bache M, Bartel F, Holzhausen HJ, Hinze R: Gains in chromosomes 7, 8q, 15q and 17q are characteristic changes in malignant but not in benign peripheral nerve sheath tumors from patients with Recklinghausen's disease. Cancer Lett. 2000, 155 (2): 181-190. 10.1016/S0304-3835(00)00426-2
Schmidt H, Wurl P, Taubert H, Meye A, Bache M, Holzhausen HJ, Hinze R: Genomic imbalances of 7p and 17q in malignant peripheral nerve sheath tumors are clinically relevant. Genes Chromosomes Cancer. 1999, 25 (3): 205-211. 10.1002/(SICI)1098-2264(199907)25:3<205::AID-GCC2>3.0.CO;2-B
Skotheim RI, Kallioniemi A, Bjerkhagen B, Mertens F, Brekke HR, Monni O, Mousses S, Mandahl N, Saeter G, Nesland JM, Smeland S, Kallioniemi OP, Lothe RA: Topoisomerase-II alpha is upregulated in malignant peripheral nerve sheath tumors and associated with clinical outcome. J Clin Oncol. 2003, 21 (24): 4586-4591. 10.1200/JCO.2003.07.067
Francis P, Namlos HM, Muller C, Eden P, Fernebro J, Berner JM, Bjerkehagen B, Akerman M, Bendahl PO, Isinger A, Rydholm A, Myklebost O, Nilbert M: Diagnostic and prognostic gene expression signatures in 177 soft tissue sarcomas: hypoxia-induced transcription profile signifies metastatic potential. BMC genomics. 2007, 8: 73- 10.1186/1471-2164-8-73
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 (1): 98-107. 10.1002/path.1101
Bridge RS, Bridge JA, Neff JR, Naumann S, Althof P, Bruch LA: Recurrent chromosomal imbalances and structurally abnormal breakpoints within complex karyotypes of malignant peripheral nerve sheath tumour and malignant triton tumour: a cytogenetic and molecular cytogenetic study. J Clin Pathol. 2004, 57 (11): 1172-1178. 10.1136/jcp.2004.019026
Evans DG, Baser ME, McGaughran J, Sharif S, Howard E, Moran A: Malignant peripheral nerve sheath tumours in neurofibromatosis 1. J Med Genet. 2002, 39 (5): 311-314. 10.1136/jmg.39.5.311
Hagel C, Zils U, Peiper M, Kluwe L, Gotthard S, Friedrich RE, Zurakowski D, von Deimling A, Mautner VF: Histopathology and clinical outcome of NF1-associated vs. sporadic malignant peripheral nerve sheath tumors. J Neurooncol. 2007, 82 (2): 187-192. 10.1007/s11060-006-9266-2
Anghileri M, Miceli R, Fiore M, Mariani L, Ferrari A, Mussi C, Lozza L, Collini P, Olmi P, Casali PG, Pilotti S, Gronchi A: Malignant peripheral nerve sheath tumors: prognostic factors and survival in a series of patients treated at a single institution. Cancer. 2006, 107 (5): 1065-1074. 10.1002/cncr.22098
Cashen DV, Parisien RC, Raskin K, Hornicek FJ, Gebhardt MC, Mankin HJ: Survival data for patients with malignant schwannoma. Clin Orthop Relat Res. 2004, 69-73.
Forus A, Weghuis DO, Smeets D, Fodstad O, Myklebost O, van Kessel AG: Comparative genomic hybridization analysis of human sarcomas: I. Occurrence of genomic imbalances and identification of a novel major amplicon at 1q21-q22 in soft tissue sarcomas. Genes Chromosomes Cancer. 1995, 14 (1): 8-14. 10.1002/gcc.2870140103
Tarkkanen M, Huuhtanen R, Virolainen M, Wiklund T, Asko-Seljavaara S, Tukiainen E, Lepantalo M, Elomaa I, Knuutila S: Comparison of genetic changes in primary sarcomas and their pulmonary metastases. Genes Chromosomes Cancer. 1999, 25 (4): 323-331. 10.1002/(SICI)1098-2264(199908)25:4<323::AID-GCC3>3.0.CO;2-5
Fong SF, Dietzsch E, Fong KS, Hollosi P, Asuncion L, He Q, Parker MI, Csiszar K: Lysyl oxidase-like 2 expression is increased in colon and esophageal tumors and associated with less differentiated colon tumors. Genes Chromosomes Cancer. 2007, 46 (7): 644-655. 10.1002/gcc.20444
Akiri G, Sabo E, Dafni H, Vadasz Z, Kartvelishvily Y, Gan N, Kessler O, Cohen T, Resnick M, Neeman M, Neufeld G: Lysyl oxidase-related protein-1 promotes tumor fibrosis and tumor progression in vivo. Cancer Res. 2003, 63 (7): 1657-1666.
Berner JM, Sorlie T, Mertens F, Henriksen J, Saeter G, Mandahl N, Brogger A, Myklebost O, Lothe RA: Chromosome band 9p21 is frequently altered in malignant peripheral nerve sheath tumors: studies of CDKN2A and other genes of the pRB pathway. Genes Chromosomes Cancer. 1999, 26 (2): 151-160. 10.1002/(SICI)1098-2264(199910)26:2<151::AID-GCC7>3.0.CO;2-A
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 (11): 4132-4138.
Sharpless NE: INK4a/ARF: a multifunctional tumor suppressor locus. Mutat Res. 2005, 576 (1-2): 22-38.
Storlazzi CT, Brekke HR, Mandahl N, Brosjo O, Smeland S, Lothe RA, Mertens F: Identification of a novel amplicon at distal 17q containing the BIRC5/SURVIVIN gene in malignant peripheral nerve sheath tumours. J Pathol. 2006, 209 (4): 492-500. 10.1002/path.1998
Fukuda S, Pelus LM: Survivin, a cancer target with an emerging role in normal adult tissues. Mol Cancer Ther. 2006, 5 (5): 1087-1098. 10.1158/1535-7163.MCT-05-0375
Osaka E, Suzuki T, Osaka S, Yoshida Y, Sugita H, Asami S, Tabata K, Sugitani M, Nemoto N, Ryu J: Survivin expression levels as independent predictors of survival for osteosarcoma patients. J Orthop Res. 2007, 25 (1): 116-121. 10.1002/jor.20291
Karube K, Nabeshima K, Ishiguro M, Harada M, Iwasaki H: cDNA microarray analysis of cancer associated gene expression profiles in malignant peripheral nerve sheath tumours. J Clin Pathol. 2006, 59 (2): 160-165. 10.1136/jcp.2004.023598
Levy P, Vidaud D, Leroy K, Laurendeau I, Wechsler J, Bolasco G, Parfait B, Wolkenstein P, Vidaud M, Bieche I: Molecular profiling of malignant peripheral nerve sheath tumors associated with neurofibromatosis type 1, based on large-scale real-time RT-PCR. Mol Cancer. 2004, 3: 20- 10.1186/1476-4598-3-20
Horiuchi S, Yamamoto H, Min Y, Adachi Y, Itoh F, Imai K: Association of ets-related transcriptional factor E1AF expression with tumour progression and overexpression of MMP-1 and matrilysin in human colorectal cancer. J Pathol. 2003, 200 (5): 568-576. 10.1002/path.1387
Yamamoto H, Horiuchi S, Adachi Y, Taniguchi H, Nosho K, Min Y, Imai K: Expression of ets-related transcriptional factor E1AF is associated with tumor progression and over-expression of matrilysin in human gastric cancer. Carcinogenesis. 2004, 25 (3): 325-332. 10.1093/carcin/bgh011
Garzon R, Fabbri M, Cimmino A, Calin GA, Croce CM: MicroRNA expression and function in cancer. Trends Mol Med. 2006, 12 (12): 580-587. 10.1016/j.molmed.2006.10.006
Calin GA, Croce CM: MicroRNAs and chromosomal abnormalities in cancer cells. Oncogene. 2006, 25 (46): 6202-6210. 10.1038/sj.onc.1209910
Zhang L, Huang J, Yang N, Greshock J, Megraw MS, Giannakakis A, Liang S, Naylor TL, Barchetti A, Ward MR, Yao G, Medina A, O'Brien-Jenkins A, Katsaros D, Hatzigeorgiou A, Gimotty PA, Weber BL, Coukos G: microRNAs exhibit high frequency genomic alterations in human cancer. Proc Natl Acad Sci USA. 2006, 103 (24): 9136-9141. 10.1073/pnas.0508889103
Meza-Zepeda LA, Kresse SH, Barragan-Polania AH, Bjerkehagen B, Ohnstad HO, Namlos HM, Wang J, Kristiansen BE, Myklebost O: Array comparative genomic hybridization reveals distinct DNA copy number differences between gastrointestinal stromal tumors and leiomyosarcomas. Cancer Res. 2006, 66 (18): 8984-8993. 10.1158/0008-5472.CAN-06-1972
Wang J, Meza-Zepeda LA, Kresse SH, Myklebost O: M-CGH: analysing microarray-based CGH experiments. BMC bioinformatics. 2004, 5: 74- 10.1186/1471-2105-5-74
Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98 (9): 5116-5121. 10.1073/pnas.091062498
Lingjaerde OC, Baumbusch LO, Liestol K, Glad IK, Borresen-Dale AL: CGH-Explorer: a program for analysis of array-CGH data. Bioinformatics . 2005, 21 (6): 821-822. 10.1093/bioinformatics/bti113
We thank Erik B. Paulsen for excellent technical assistance and Bjørn E. Kristiansen for assistance on microarray production.
This work was supported by the University of Oslo (EMBIO), the Norwegian Cancer Society and the functional genomics programme (FUGE) of the Research Council of Norway. The genomic microarrays were provided by the Norwegian Microarray Consortium (NMC) at the national technology platform, supported by the functional genomics programme (FUGE) of the Research Council of Norway.
The authors declare that they have no competing interests.
SHK performed the array CGH experiments, all analyses and drafted the manuscript, MS performed the quantitative real-time RT-PCR experiments, HOO collected the clinical data, HMN performed the microarray expression profiling, BB revised pathologically all the clinical samples, OM and LAMZ conceived of the study, participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
About this article
- Soft Tissue Sarcoma
- Malignant Peripheral Nerve Sheath Tumor
- Array Comparative Genomic Hybridization
- Gene TOP2A
- Recurrent Copy Number Alteration