Email updates

Keep up to date with the latest news and content from BMC Cancer and BioMed Central.

Open Access Research article

Genome-wide array comparative genomic hybridization analysis reveals distinct amplifications in osteosarcoma

Tsz-Kwong Man1, Xin-Yan Lu1, Kim Jaeweon2, Laszlo Perlaky1, Charles P Harris1, Shishir Shah2, Marc Ladanyi3, Richard Gorlick4, Ching C Lau1 and Pulivarthi H Rao1*

  • * Corresponding author: Pulivarthi H Rao phrao@txccc.org

  • † Equal contributors

Author Affiliations

1 Texas Children's Cancer Center, Baylor College of Medicine, Houston, TX, USA

2 Spectral Genomics, Houston, TX, USA

3 Department of Pathology, Memorial Sloan Kettering Cancer Center, New York, New York, USA

4 Department of Pediatrics, Memorial Sloan Kettering Cancer Center, New York, New York, USA

For all author emails, please log on.

BMC Cancer 2004, 4:45  doi:10.1186/1471-2407-4-45


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2407/4/45


Received:4 May 2004
Accepted:7 August 2004
Published:7 August 2004

© 2004 Man et al; licensee BioMed Central Ltd.

This is an open-access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Osteosarcoma is a highly malignant bone neoplasm of children and young adults. It is characterized by extremely complex karyotypes and high frequency of chromosomal amplifications. Currently, only the histological response (degree of necrosis) to therapy represent gold standard for predicting the outcome in a patient with non-metastatic osteosarcoma at the time of definitive surgery. Patients with lower degree of necrosis have a higher risk of relapse and poor outcome even after chemotherapy and complete resection of the primary tumor. Therefore, a better understanding of the underlying molecular genetic events leading to tumor initiation and progression could result in the identification of potential diagnostic and therapeutic targets.

Methods

We used a genome-wide screening method – array based comparative genomic hybridization (array-CGH) to identify DNA copy number changes in 48 patients with osteosarcoma. We applied fluorescence in situ hybridization (FISH) to validate some of amplified clones in this study.

Results

Clones showing gains (79%) were more frequent than losses (66%). High-level amplifications and homozygous deletions constitute 28.6% and 3.8% of tumor genome respectively. High-level amplifications were present in 238 clones, of which about 37% of them showed recurrent amplification. Most frequently amplified clones were mapped to 1p36.32 (PRDM16), 6p21.1 (CDC5L, HSPCB, NFKBIE), 8q24, 12q14.3 (IFNG), 16p13 (MGRN1), and 17p11.2 (PMP22 MYCD, SOX1,ELAC27). We validated some of the amplified clones by FISH from 6p12-p21, 8q23-q24, and 17p11.2 amplicons. Homozygous deletions were noted for 32 clones and only 7 clones showed in more than one case. These 7 clones were mapped to 1q25.1 (4 cases), 3p14.1 (4 cases), 13q12.2 (2 cases), 4p15.1 (2 cases), 6q12 (2 cases), 6q12 (2 cases) and 6q16.3 (2 cases).

Conclusions

This study clearly demonstrates the utility of array CGH in defining high-resolution DNA copy number changes and refining amplifications. The resolution of array CGH technology combined with human genome database suggested the possible target genes present in the gained or lost clones.

Background

Osteosarcoma (OS) is a primary malignant tumor of bone arising from primitive bone-forming mesenchymal cells and it accounts for approximately 60% of malignant bone tumors in the first two decades of life [1]. These tumors typically arise in the metaphyseal regions of long bones, with the distal femur, proximal tibia and proximal humerus. A significant number of osteosarcomas are of conventional type which can be subdivided into three major categories based on their predominant differentiation of tumor cells: osteoblastic, chondroblastic, and fibroblastic. Currently, only the histological response (degree of necrosis) to therapy represent gold standard for predicting the outcome in a patient with non-metastatic osteosarcoma at the time of definitive surgery [2]. Patients with lower degree of necrosis have a higher risk of relapse and poor outcome even after chemotherapy and complete resection of the primary tumor. Therefore, a better understanding of the underlying molecular genetic events leading to tumor initiation and progression could result in the identification of potential diagnostic and therapeutic targets.

Chromosomal aberrations in osteosarcoma are highly complex and characterized by high frequency of amplifications. These amplifications may result in the overexpression of genes and contribute to the genomic instability in osteosarcoma. The identification of genes within the amplified sites is crucial for understanding the biology and clinical behavior of osteosarcoma. Until, recently gene amplification has been detected by PCR, southern blot analysis or FISH-based approach using gene specific probes. These techniques are inherently restricted to the previously known amplified genes in the genome. In contrast, genome-wide screening of amplified chromosomal regions with CGH has become an important tool for the detection of amplified regions in the tumor genome. So far published chromosomal CGH studies in osteosarcoma have identified several high-level chromosomal amplifications at 1p22, 1p31, 1q21, 1q23, 2q24, 3p25, 3q26, 6q24.3, 4q12, 5p14-p15, 5q33, 6p12-p21, 6q24.3, 7p21-p22, 8q12-q23, 10p21, 10q11.1, 10q22, 11q13, 11q23, 12p13, 12q12-q15, 17p11.2, 17q21, 18q22, 19p13.1 and 20p11.2 [3-7]. However, conventional CGH has limited sensitivity and resolution (~10–15 megabases) because of its dependence on the morphology of metaphase chromosomes. In addition, extensive follow-up work is required to identify candidate genes after regions of gain or loss have been identified. Recently, novel method termed as array-based comparative genomic hybridization (array CGH) has been described, which enables high throughput quantitative measurement of high-resolution DNA copy number changes throughout the genome [8]. This method is based on hybridization of differentially labeled test and reference DNAs to an array of mapped human genomic DNA fragments (~100–200 kb) and has been recently applied to human and mouse tumors [9-14]. To identify high-resolution copy number, we used array CGH to the panel of 48 tumors. The resolution of array CGH technology combined with human genome database not only allowed a precise identification of amplicons but also suggested the possible target genes within the amplicons.

Methods

Patient samples

A total of 48 tumors from 42 patients (20 males and 22 females) were collected from the Texas Children's Cancer Center, Houston, TX (tumors 193, 196, 204, 209, 221, 226, 248, 274, 295, 311, 326, 341, 345, 360, 400, 464, 481, 501, 527, 591 and 606) and Memorial Sloan Kettering Cancer Center, New York (tumors 06, 15, 24, 25, 27, 29, 32, 34, 40, 48, 68, 76, 78, 79, 80, 82, 83, 85, 88, 95, 98, 99, 102, 123, 423, 425, and 474). All tissues in this study were obtained after IRB approved informed consents were signed. The age at diagnosis ranged from 5 years to 71 years at diagnosis. The histological information of 42 patients is presented in Table 1.

Table 1. Histological information on 48 osteosarcoma samples

Array CGH

The array used in this study consists of 967 human BACs, which were spaced approximately 3 megabase across the whole genome. These arrays were obtained from Spectral Genomics, Houston, TX. The experiments were performed according to the manufacturer's protocol. Arrays were pre-hybridized with human Cot-I DNA (GIBCO Invitrogen, Carlsbad, CA) and salmon testes DNA to block the repetitive sequences on BACs. One microgram of normal DNA (reference) and tumor DNA (test) was labeled with cy5-dUTP and cy3-dUTP respectively, by random priming. To avoid dye bias, we performed dye swap experiments for each sample. The probe mixture is dissolved in hybridization mixture, denatured, cooled, and mounted with 22 × 60 mm coverslip. Hybridizations were performed in sealed chambers for 16–20 hours at 60°C. After post hybridization washes, arrays were rinsed, dried with compressed air, and scanned into two 16-bit TIFF image files using Gene Pix 4000A two-color fluorescent scanner (Axon Instruments, Inc., Union City, CA) and quantitated using GenePix software (Axon Instruments, Union City, CA).

Data processing and analysis

After scanning of the slide, the fluorescent intensities of the green and red channels were background subtracted. The resulting values were normalized by intensity based local weighted regression method (Lowess) to correct for systematic bias in dye labeling and fluorescent intensity [15]. Then the ratio of the red/green channel of each clone was calculated and log base 2 transformed (log ratios). Each experiment was repeated once with dye reversal to addressing the confounding effect of the dye and experiment. The average of the dye-reversal experiment pair was calculated by reversing the sign of one experiment so that the log ratio reflects the tumor versus normal ratio.

We developed a new analytical method called invariant analysis to define the significant copy number changes. This method is designed to: i) increase the power of the analysis by combining all the cases in our dataset to define an invariant population (unchanged population); and, ii) to address the signal to noise differences among individual cases due to sample and hybridization variability. Our goal is to define a set of unchanged clones that can be used to calculate the upper and lower bound thresholds of the log ratios for the unchanged population in each experiment. First, we calculated the variance of each clone from all the experiments. We computed the p-values of the each clone by comparing to the clone with median variance using chi-square distribution http://linus.nci.nih.gov/~brb/TechReport.htm webcite. The clones that have p-value greater than preset cutoff 0.9 were considered as invariant clone set, i.e. clones that do not vary significantly in all experiments. Then the mean and standard deviation of the log ratios of these invariant clones in each experiment were calculated. The clones with log ratios that exceed mean +/- 2 × SD of the invariant set were considered gains and losses, respectively. For amplification and homozygous deletions, clones were defined to have at least 2 fold of the upper bound threshold and 4-fold of lower bound threshold, respectively. The gene(s) present in the clones were identified using UCSC browser http://genome.ucsc.edu/cgi-bin/hgGateway ?org=human webcite by downloading gene table (refFlat) from human gene assembly, July 2003. We search the candidate genes based on linear mapping position, which include 100 kb up and downstream from the clone center position. The supplemental data for this article is available at: https://bigfile.bcm.tmc.edu/outbound/prao/supplemental%20data%20Man%20et%20al.xls webcite

Statistical analysis

Significant clones in 6p, 8q, 12q and 17p amplicons were calculated using 2-sample t-test with randomized variance model http://linus.nci.nih.gov/~brb/TechReport.htm webcite. The experiments in each of the two groups, amplification and normal, used for comparison were defined based on the invariant analysis (see above). The clones that have p <0.001 were considered as significant. We chose a stringent cutoff to minimize the multiple testing problem.

FISH

FISH was performed to validate and quantify chromosomal amplicons using clones from 6p12-p21 (RP11-91E11, AL391415, RP11-81F7, RP11-79I2, RP11-90H17 and RP11-79F13), 8q24.3 (RP11-89K10), and 17p11.2-p12 (RP11-64B12, RP11-89K6 and RP11-189D22 on tumors metaphase/interphase cells from cases 274, 364, 425, 426, 527 and 628. We confirmed the map positions of all clones on normal human metaphase cells by FISH. The bacterial artificial chromosome (BAC) clones, and centromeric clone from 6 (pEDZ6) were labeled with Spectrum Red or Spectrum Green (Vysis, Downers Grove, IL) by nick translation. Hybridization and FISH analysis was performed as described previously [16].

Results

To define the gains and losses in our experiments, we used invariant analysis for the first time to describe genomic changes by array CGH. In this method, we defined an invariant clone set that has low variance of log ratios among all the array experiments. After the mean and standard deviation of the log ratios in the invariant set of each experiment were calculated, clones that have higher or lower log ratios than the mean +/- 2SD of the invariant set (upper bound and lower bound) were used to define gains and losses. We chose to use this method because it addresses some of the shortcomings of the modeling method, such as using all information provided in a set of experiment to determine the unchanged population instead of using one experiment at a time. However, the variation of each experiment is accounted for because the thresholds are calculated using the invariant set from each experiment. It also does not require a separate reference set for comparison. Finally, it provides an adjustable cutoff to optimize the thresholds to the training data, if provided.

The amplified and homozygously deleted clones were defined to have at least 2 fold of the upper bound and 4-fold of lower bound, respectively. Figure 1 summarizes the high-resolution DNA copy number changes identified by array CGH in 48 osteosarcomas derived from 42 patients. Copy number changes were detected involving small genomic regions, whole chromosomes, and chromosomal arms showing homozygous deletions and high-level amplifications.

thumbnailFigure 1. Frequency of DNA copy number changes identified by array CGH in 48 osteosarcomas. The gains and losses are depicted as green and red color bars respectively. Clones are arranged from chromosome 1 to 22 and within each chromosome on the basis of UCSF mapping position.

Overview of genomic profiles

Copy number changes excluding clones from sex chromosomes were involved in a significant fraction of most tumor genome. The estimated average genomic distance between clones was ~3–4 Mb. The frequency of clones showing gains (79%) was greater than losses (66%). High-level amplifications and homozygous deletions constitute 28.6% and 3.8% of tumor genome respectively. The most frequently deleted clones were identified from the chromosomal regions 2q31.1, 3p14.1, 4p16.2, 6q12, 6q21, 6q27, 7q35, 10p15.1, 10q22-q23, 10q25-q26, 11q25, 13q12.2, 13q14.3, 13q22.1, 17p13.3 and 17q12 (Table 2). Most frequently gained clones were mapped to chromosome 1p36, 4p16, 6p12-p21, 8q21, 8q23-q24, 12q14.3, 16p13.3, 16q24.3, 17p11-p12, 19p13.3 and 21q22.3 (Table 3). We explored the possible statistical relationship between copy number alterations and histological and clinical parameters. We found no significant relationship between copy number changes and primary/metastatic disease, or histological type or histological response. This may be due to the involvement of large number of genomic loci and insufficient sample size.

Table 2. Most frequently lost clones.

Table 3. Most frequently gained clones.

Homozygous deletions were noted for 32 clones (3.8%). Recurrent homozygous deletions were noted for 7 clones that are were mapped to 1q25.1 (4 cases), 3p14.1 (4 cases), 13q12.2 (2 cases), 4p15.1 (2 cases), 6q12 (2 cases), 6q12 (2 cases) and 6q16.3 (2 cases). Figure 2A is showing a homozygous deletion at 3p14.1 in tumor 06. Loss of 6q12 region was noted in 35% of the osteosarcomas. This region was covered with four clones spanning ~4.2 Mb. Two tumors (tumor 27 and 345) showed low intensity ratios indicting homozygous deletions in this region, one tumor (tumor 345) showed all 4 deleted clones spanning ~4.2 Mb with RP1-129L7 having the lowest ratio intensity decrease. In another case (tumor 27), two clones (RP1-46B1 and RP1-129L7) showed decreased intensity ratios indicating homozygous deletions. Both these clones spanning approximately 2.6 Mb of 6q12 region.

thumbnailFigure 2. A representative chromosome profiles showing homozygous deletions (A, E) and high-level amplifications (B- E). The clones showing homozygous deletions from 3p14.1 (RP11-89A12) and 17q12 (RP11-79O9) in tumor 06 (A, E) and high-level amplifications from 6p12-p21 in tumor 248(B), 8q12 in tumors 06 and 341 and 8q23-q24 in tumor 06, (C), 12q13-q15 in tumor 48, (D) and 17p11.2 in tumor 06(E), are shown as log ratios (Y-axis). The clones on each chromosome are arranged (pter to qter) on the basis of UCSC mapping positions.

Amplification is a frequent phenomenon in osteosarcoma

Previous studies using CGH have identified several chromosomal amplification sites in osteosarcoma. Because of the limitation of the method, it fails to pinpoint the precise site of amplicon. However, the present study by array CGH has identified 238 clones (28.6%) with high-level amplifications. Recurrent amplifications were noted in ~37% of the total amplified clones (Figure 3). These amplified clones were mapped to 1p22, 1p31.1 (ROR1), 1p36.1 (PRDM16), 1q21, 1q23 (TNFF6), 2q24, 3p25, 3q26.1, 4p16.3, 5p14, 5q33, 6p11.2-p21, 7p21, 8q12.1, 8q24.13, 10p21, 10q11.1, 10q22 (KCNMA1), 11q13, 11q23 (GRIK4), 12q12, 12q13-q15, 12q21-q21.33, 17p11.2-p12, 17q21 (NGFR), 18q22, and 19p13.1 (NFAT). Of these amplified sites, 6p11.2-p21, 8q12.1, 8q24.13, 12q12, 12q13-q15, 12q21-q21.33, 16p13 and 17p11.2-p12 were frequent.

thumbnailFigure 3. Ideogram showing recurrent homozygous deletions (left) and high-level amplifications (right) identified by array CGH in 48 cases. The gene(s) contained within the BAC clone are shown in parentheses of the respective clones. Gene(s) present in the BAC clones were identified using UCSC browser http://genome.ucsc.edu/cgi-bin/hgGateway?org=human webcite by downloading gene table (refFlat) from human gene assembly, July 2003. The ideograms for chromosomes 6, 8, 12 and 17 are shown separately.

Gain of clones from 6p12-p21 regions was noted in 33/48 (~65%) cases analyzed. High-level amplification of the clones from same region was noted in 25% of the cases by array CGH. We found that most of the cases with amplification of 6p12-p21 displayed either increased or slightly varying degree of copy number increase across the 6p12-p21 region. The combined log ratios from all the cases defined the boundaries of amplification between RP3-329A5 and RP11-79F13. The amplicon spans approximately 9.4 Mb with amplification peak for clone RP11-81F7. Further, we used FISH to validate 6p amplicon on tumor metaphase and interphase cells from cases 274, 364, 426 and 527. Increased copy numbers for clones RP11-91E11, AL391415, RP11-81F7, RP11-79I2, RP11-90H17 and RP11-79F13 were noted in interphase cells with maximum copy number increase for clone RP11-81F7 (Figure 4A). This was consistent with amplification peak for clone RP11-81F7 in the tumors profiled by array CGH (Figure 2B). In addition, we used 2-sample t-test with randomized variance model to define significant clones from 6p12-p21 amplicon. By this method, we identified RP11-79F13 (p = 0.00000007), RP11-79I2 (p = 0.00000007) and RP11-81F7 (p = 0.00000007) as statistically significant clones.

thumbnailFigure 4. FISH validation of some of the high-level amplifications (6p12.1, 8q24.3 and 17p11.2) identified by array CGH. Interphase cells hybridized with centromere 6 (red)/RP11-81F7 (green) in case 274 (A), RP11-89K10 (red) in case 527 (B) and RP11-189D22 (red) in case 364 (C). The ploidy of these cases was determined based on the modal chromosome number of the respective cases, e.g. diploid (case 426) triploid (cases 274 and 364), and tetraploid (case 527).

Most cases with 8q gain, displayed varying degree of copy number increase predominantly from 8q12.1 (16.9%), 8q21.13 (29%), and 8q24.3 (35%). High-level amplifications were also noted from 8q12.1 (RP11-550I15 – 6.3%; Figure 2C), 8q21.13 (RP11-89H1 – 6.3%), 8q24.3 (RP11-89K10 – 6.3%) and RP11-637F16 (12.5%). FISH using clone RP11-89K10 (p = 0.00049) on interphase cells from case 527 confirmed the amplification (10–12 copies) (Figure 4B).

Amplification of 12q was noted in 14/51 (~27%) tumors analyzed by array CGH. Three distinct amplicons – AMP1 (12q12), AMP2 (12q14.1) and AMP3 (12q21.33) were noted across the entire long arm of chromosome 12 (Figure 2D). Of these 14 cases, four of them (80, 123, 248, 341) displayed all three amplicons. The AMP1 was noted in 10 cases covering 1.8 Mb region between RP11-91K15 and RP11-90I21 with peak amplification for clone RP11-91K15 (p = 0.00000004). Another amplicon (AMP2) was noted 24.48 Mb distal to AMP1 between RP11-91K23 and RP11-89P15. The AMP3, which was 23.3 Mb distal to AMP2 containing RP11-89F6.

Amplification of 17p11.2 was noted in 27% of the cases analyzed by array CGH. The amplicon was composed of three clones RP11-64B12 (p = 0.0000014), RP11-89K6 (p = 0.00000005) and RP11-189D22 (p = 0.0000001) and covering 3.7 Mb region on the short arm of chromosome 17 (Figure 2E). We used these three clones as FISH probes to validate 17p amplicon in tumors 274, 364, 425 and 628 on interphase/ metaphase cells. The distribution of copy number for this amplicon in all the cases ranged from 4–14 copies with peak amplification for clone RP11-189D22 (10–14 copies), followed by and RP11-89K6 (8–10 copies) RP11-64B12 (6–8 copies) (Figure 4C).

Discussion

This study represents the first application of genome-wide copy number changes by array CGH in osteosarcoma. Recent studies in breast, renal and bladder cancer showed the potential assessment of this technology in detecting high-resolution copy number changes [9,11,14]. This approach will augment the identification of cancer causing genes by relating the clone information directly with sequence information from human genome database. In this study, we used array CGH to screen for high-resolution DNA copy number changes and precise identification of amplifications in a panel of 48 osteosarcomas.

Gene amplification is an important genetic mechanism in human cancers, as it clearly associated with tumor progression and has a prognostic significance and has even provided a target for therapeutics [17,18]. These amplifications are often seen at the cytogenetic level as homozygously staining regions (hsrs) or double minute chromosomes (dms). However, cytogenetic recognition of amplifications doesn't contribute to the mapping and identification of amplified DNA sequences. The advent of CGH points an ever-increasing number of chromosomal amplifications in various tumors. These amplifications contribute to the genomic instability in tumors. We have recently shown that the mutation of p53 significantly correlates with genome-wide DNA instability and seems to represent a major genetic factor contributing to the extremely high levels of genomic instability found in high-grade osteosarcomas [19].

Our analysis have identified frequently amplified clones from 6p11.2-p21, 8q12.1, 8q24.13, 12q12, 12q13-q15, 12q21-q21.33, 16p13 and 17p11.2-p12. Amplification of clones from 6p12-p21 region was noted in 25% of the cases analyzed. This was consistent with the previously published results by CGH. By array CGH, we refined the 6p amplicon to 9.4 Mb with amplification peak for clone RP11-81F7. We recently demonstrated the origin of 6p amplicon as consequence of tandem duplication of clones RP11-81F7 and RP11-79F13 [7]. Based on combined array CGH and FISH analysis suggest CDC5L, HSPCB, and NFKBIE, and HGNC and MRPL14 are the target genes from 6p12-p21 amplicon. Of these genes, CDC5L may be an important gene in cancer because of its role as a positive cell cycle regulator for G2/M transition[20]. Consistent with our analysis, overexpression of HSPCB was shown recently by cDNA microarray studies on osteosarcoma [21]. This protein was shown to play an important role in assemble/disassembly of tubulin by inhibiting tubulin polymerization.

High-level amplifications were also noted from 8q12.1 (RP11-550I15 – 6.3%), 8q21.13 (RP11-89H1 – 6.3%), 8q24.3 (RP11-89K10 – 6.3%) and RP11-637F16 (12.5%). There were no candidate genes present in clones RP11-550I15, RP11-89H1 and RP11-637F16, but clone RP11-89K10 contained NSE2 (breast cancer membrane protein 101 kDa) gene.

High-level amplification of clones on 12q revealed three distinct sites of amplifications – AMP1 (12q12), AMP2 (12q14.1) and AMP3 (12q21.33). Pervious studies have shown the amplification GLI, CHOP, SAS, HMGI-C, CDK4, HDM2, and PRIM1 from 12q13-q15 region in osteosarcoma [22,23]. The present array CGH analysis identified a possible target gene IFNG from AMP2 (RP11-298M11; p = 0.0000001), which is physically mapped close to the HDM2 oncogene locus[24]. Previous studies demonstrated that T-cell production of IFNG strongly suppresses osteoclastogenesis by interfering with the RANKL-RANK signaling pathway. IFNG induces rapid degradation of the RANK adaptor protein, TRAF6, resulting in strong inhibition of the RANKL-induced activation of the transcription factor NFKB and JNK [25]. The AMP3, which was 23.3 Mb distal to AMP2 containing RP11-89F6. Our analysis from AMP3 revealed two interesting candidate genes: transcription factor ELK3 and PCTAIRE protein kinase 2 (PCTK2). ELK3 is a member of the ETS-domain transcription factor family and the protein is activated by signal-induced phosphorylation [26]. The protein encoded by PCTK2 belongs to the cdc2/cdkx subfamily of the ser/thr family of protein kinases and play an important role in the regulation of the mammalian cell cycle [27]. High-level amplification of three clones from 12p13 was noted in case 27 and the amplicon span 4.6 Mb with peak amplification for clone RP11-89D16. No candidate genes contained with in this BAC. Amplification 12p has been reported previously in 9/19 high-grade osteosarcomas by CGH. Recent FISH analysis has identified the amplification of CCND2, ETV6, and KRAS2 from 12p region [28].

Amplification of 17p11.2 was noted in 27% of the cases analyzed by array CGH. Our array CGH analysis has identified three clones with high-level amplifications that spans ~3.7 Mb region on 17p11.2. Several candidate genes were identified within these clones (TPP3A, SMCR5, DRG2, FL11, MYCD, SOX 17, ELAC2, and PMP22). Recent studies have shown the amplification of some of the genes identified in the present study (PMP22, and TOP3A) from 17p11.2-p12 in high-grade OS by semi-quantitative PCR and cDNA microarrays [29,30].

The present array CGH analysis has identified seven recurrent clones exhibiting homozygous deletions from 1q25.1, 3p14.1, 13q12.2, 4p15.1, 6q12, 6q12 and 6q16.3. These chromosomal regions were consistent with previously reported studies by loss of heterozygosity (LOH) and CGH [3-7,31]. The clone, RP11-90M15 (13q12.2) contain possible candidate gene MTMR6, a protein-tyrosine phosphatase gene and shown to be present within a cloned region that encompasses a translocation breakpoint t(8;13) in an atypical myoproliferative disorder [32]. Homozygous deletions of two clones spanning approximately 2.6 Mb of 6q12 region containing candidate genes – nuclear fragile X mental retardation protein interacting protein 1 pseudogene (NUFIP1P) and BAI3 gene (brain-specific angiogenesis inhibitor gene), which is to homologous to BAI1 and shown to suppress glioblastoma [33].

Conclusions

In summary, high resolution array-based CGH revealed large number of chromosomal aberrations previously identified in osteosarcoma by chromosomal CGH and conventional cytogenetic methods. The present study allowed precise identification of smaller DNA copy number alterations, which suggest the presence of specific target genes in osteosarcoma. Although this study suggested several possible target genes from amplified regions from 6p, 8q, 12q and 17p, but these genes should be validated by other molecular and immunohistochemical approaches on well-defined large patient samples. Further, interaction or association studies between small genomic losses and gains will facilitate the identification of new genetic pathways in the pathogenesis of osteosarcoma.

Competing interests

None declared.

Authors contributions

TKM and KJ have contributed towards the data analysis. LP, ML, RG, and CL were assisted in sample collection and clinical information of the patients. X-YL has involved in array CGH experiments and data collection. CPH has involved in extracting the gene information from BAC clones. SS has provided the arrays used in this study. PHR was involved in the planning, and organization of the project.

Acknowledgements

We thank Prof Rocchi, University of Bari, Italy for providing chromosome 6 centromeric clone and Dr. Richard Simon and Amy Peng for BRB ArrayTools software. This work was supported in part by NIH grant CA88126 and by grants from the Dunn Foundation and the Kleberg Foundation to C.C.L.

References

  1. Huvos A: Bone Tumors: Diagnosis, Treatment and Prognosis,. Edited by Huvos A. Philadelphia, WB Saunders; 1991.

  2. Provisor AJ, Ettinger LJ, Nachman JB, Krailo MD, Makley JT, Yunis EJ, Huvos AG, Betcher DL, Baum ES, Kisker CT, Miser JS: Treatment of nonmetastatic osteosarcoma of the extremity with preoperative and postoperative chemotherapy: a report from the children's cancer group.

    J Clin Oncol 1997, 15:76-84. PubMed Abstract OpenURL

  3. Tarkkanen M, Karhu R, Kallioniemi A, Elomaa I, Kivioja AH, Nevalainen J, Bohling T, Karaharju E, Hyytinen E, Knuutila S: Gains and losses of DNA sequences in osteosarcomas by comparative genome hybridization.

    Cancer Res 1995, 55:1334-1338. PubMed Abstract OpenURL

  4. Tarkkanen M, Elomaa I, Blomqvist C, Kivioja AH, Kellokumpu-Lehtinen P, Bohling T, Valle J, Knuutila S: DNA sequence copy number increase at 8q: a potential new prognostic marker in high-grade osteosarcoma.

    Int J Cancer 1999, 84:114-121. PubMed Abstract | Publisher Full Text OpenURL

  5. Stock C, Kager L, Fink F-M, Gadner H, Ambros PF: Chromosomal regions involved in the pathogenesis of osteosarcoma.

    Genes Chromosomes Cancer 2000, 28:329-336. PubMed Abstract | Publisher Full Text OpenURL

  6. Zielenska M, Bayani J, Pandira A, Toledo S, Marrano P, Andrade J, Petrilli A, Thorner P, Sorensen P, Squire JA: Comparative genomic hybridization analysis identifies gains of 1p35-36 and chromosome 19 in osteosarcoma.

    Cancer Genet Cytogenet 2001, 130:14-21. PubMed Abstract | Publisher Full Text OpenURL

  7. Lau CC, Harris CP, Lu XY, Perlaky L, Gogineni S, Chintagumpala M, Hicks J, Johnson ME, Davino NA, Huvos AG, Meyers PA, Healy JH, Gorlick R, Rao PH: Frequent amplification and rearrangement of chromosomal bands 6p12-p21 and 17p11.2 in osteosarcoma.

    Genes Chromosomes Cancer 2004, 39:11-21. PubMed Abstract | Publisher Full Text OpenURL

  8. Pinkel D, Segraves R, Sudar D, Clark S, Poole I, Koebel D, Collins C, Kuo WL, Chen Zhai Y, Dairkee SH, Ljung BM, Gray JW, Albertson D: High resolution analysis of DNA copy number variation using comparative genomic hybridization to microarrays.

    Nature Genet 1998, 20:207-211. PubMed Abstract | Publisher Full Text OpenURL

  9. Albertson DG, Ylstra B, Segraves R, Collins C, Dairkee SH, Kowbel D, Kuo WL, Gray JW, Pinkel D: Quantitative mapping of amplicon structure by array CGH identifies CYP24 as a candidate oncogene.

    Nat Genet 2000, 25:144-146. PubMed Abstract | Publisher Full Text OpenURL

  10. Cai WW, Mao JH, Chow CW, Damani S, Balmain A, Bradley A: Genome-wide detection of chromosomal imbalances in tumors using BAC microarrays.

    Nat Biotechnol 2002, 20:393-6. PubMed Abstract | Publisher Full Text OpenURL

  11. Wilhelm M, Veltman JA, Olshen AB, Jain AN, Moore DH, Presti JC Jr, Kovacs G, Waldman FM: Array-based comparative genomic hybridization for the differential diagnosis of renal cell cancer.

    Cancer Res 2002, 62:957-60. PubMed Abstract | Publisher Full Text OpenURL

  12. Hackett CS, Hodgson JG, Law ME, Fridlyand J, Osoegawa K, de Jong PJ, Nowak NJ, Pinkel D, Albertson DG, Jain A, Jenkins R, Gray JW, Weiss WA: Genome-wide array CGH analysis of murine neuroblastoma reveals distinct genomic aberrations, which parallel those in human tumors.

    Cancer Res 2003, 63:5266-73. PubMed Abstract | Publisher Full Text OpenURL

  13. O'Hagan RC, Brennan CW, Strahs A, Zhang X, Kannan K, Donovan M, Cauwels C, Sharpless NE, Wong WH, Chin L: Array comparative genome hybridization for tumor classification and gene discovery in mouse models of malignant melanoma.

    Cancer Res 2003, 63:5352-5356. PubMed Abstract | Publisher Full Text OpenURL

  14. Veltman JA, Fridlyand J, Pejavar S, Olshen AB, Korkola JE, DeVries S, Carroll P, Kuo WL, Pinkel D, Albertson D, Cordon-Cardo C, Jain AN, Waldman FM: Array-based comparative genomic hybridization for genome-wide screening of DNA copy number in bladder tumors.

    Cancer Res 2003, 63:2872-80. PubMed Abstract | Publisher Full Text OpenURL

  15. Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP: Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation.

    Nucleic Acids Res 2002, 30:e15. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Singh B, Gogineni SK, Sacks PG, Shaha AR, Shah JP, Stoffel A, Rao PH: Molecular cytogenetic characterization of head and neck squamous cell carcinoma and refinement of 3q amplification.

    Cancer Res 2001, 61:4506-4513. PubMed Abstract | Publisher Full Text OpenURL

  17. Lengauer C, Kinzler KW, Vogelstein B: Genetic instabilities in human cancers.

    Nature 1998, 396:643-649. PubMed Abstract | Publisher Full Text OpenURL

  18. Schwab M: Oncogene amplification in solid tumors.

    Review Semin Cancer Biol 1999, 9:319-325. Publisher Full Text OpenURL

  19. Overholtzer M, Rao PH, Favis R, Lu XY, Elowitz MB, Barany F, Ladanyi M, Gorlick R, Levine AJ: The presence of p53 mutations in human osteosarcomas correlates with high levels of genomic instability.

    Proc Natl Acad Sci U S A 2003, 100:11547-52. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Bernstein HS, Coughlin SR: A mammalian homolog of fission yeast Cdc5 regulates G2 progression and mitotic entry.

    J Biol Chem 1998, 273:4666-4671. PubMed Abstract | Publisher Full Text OpenURL

  21. Wolf M, El-Rifai W, Tarkkanen M, Kononen J, Serra M, Eriksen EF, Elomaa I, Kallioniemi A, Kallioniemi OP, Knuutila S: Novel findings in gene expression detected in human osteosarcoma by cDNA microarray.

    Cancer Genet Cytogenet 2000, 123:128-32. PubMed Abstract | Publisher Full Text OpenURL

  22. Wei G, Lonardo F, Ueda T, Kim T, Huvos AG, Healey JH, Ladanyi M: CDK4 gene amplification in osteosarcoma: reciprocal relationship with INK4A gene alterations and mapping of 12q13 amplicons.

    Int J Cancer 1999, 80:199-204. PubMed Abstract | Publisher Full Text OpenURL

  23. Yotov WV, Hamel H, Rivard GE, Champange MA, Russo PA, Leclecrc JM, Bernstein ML, Levy E: Amplification of DNA Primase 1 (PRIM1) in human osteosarcoma.

    Genes Chromosomes Cancer 1999, 26:62-69. PubMed Abstract | Publisher Full Text OpenURL

  24. Bureau JF, Bihl F, Brahic M, Le Paslier D: The gene coding for interferon-gamma is linked to the D12S335 and D12S313 microsatellites and to the MDM2 gene.

    Genomics 1995, 28:109-112. PubMed Abstract | Publisher Full Text OpenURL

  25. Takayanagi H, Ogasawara K, Hida S, Chiba T, Murata S, Sato K, Takaoka A, Yokochi T, Oda H, Tanaka K, Nakamura K, Taniguchi T: T-cell-mediated regulation of osteoclastogenesis by signaling cross-talk between RANKL and IFN-gamma.

    Nature 2000, 408:600-605. PubMed Abstract | Publisher Full Text OpenURL

  26. Lopez M, Oettgen P, Akbarali Y, Dendorfer U, Libermann TA: ERP, a new member of the ets transcription factor/oncoprotein family. cloning, characterization, and differential expression during B-lymphocyte development.

    Molec Cell Biol 1994, 14:3292-3309. PubMed Abstract | PubMed Central Full Text OpenURL

  27. Meyerson M, Enders GH, Wu C-L, Su L-K, Gorka C, Nelson C, Harlow E, Tsai L-H: A family of human cdc2-related protein kinases.

    EMBO J 1992, 11:2909-2917. PubMed Abstract OpenURL

  28. Gisselsson D, Palsson E, Hoglund M, Domanski H, Mertens F, Pandis N, Sciot R, Dal Cin P, Bridge JA, Mandahl N: Differentially amplified chromosome 12 sequences in low-and high-grade osteosarcoma.

    Genes Chromosomes Cancer 2002, 33:133-40. PubMed Abstract | Publisher Full Text OpenURL

  29. van Dartel M, Cornelissen PW, Redeker S, Tarkkanen M, Knuutila S, Hogendoorn PC, Westerveld A, Gomes I, Bras J, Hulsebos TJ: Amplification of 17p11.2 approximately p12, including PMP22, TOP3A, and MAPK7, in high-grade osteosarcoma.

    Cancer Genet Cytogenet 2002, 139:91-96. PubMed Abstract | Publisher Full Text OpenURL

  30. Squire JA, Pei J, Marrano P, Beheshti B, Bayani J, Lim G, Moldovan L, Zielenska M: High-resolution mapping of amplifications and deletions in pediatric osteosarcoma by use of CGH analysis of cDNA microarrays.

    Genes Chromosomes Cancer 2003, 38:215-25. PubMed Abstract | Publisher Full Text OpenURL

  31. Yamaguchi T, Toguchida J, Yamamuro T, Kotoura Y, Takada N, Kawaguchi N, Kaneko Y, Nakamura Y, Sasaki MS, Ishizaki K: Allelotype analysis in osteosarcoma: frequent allele loss on 3q, 13q, 17q, and 18q.

    Cancer Res 1992, 52:2419-2423. PubMed Abstract OpenURL

  32. Laporte J, Blondeau F, Buj-Bello A, Tentler D, Kretz C, Dahl N, Mandel J-L: Characterization of the myotubularin dual specificity phosphatase gene family from yeast to human.

    Hum Molec Genet 1998, 7:1703-1712. PubMed Abstract | Publisher Full Text OpenURL

  33. Shiratsuchi T, Nishimori H, Ichise H, Nakamura Y, Tokino T: Cloning and characterization of BAI2 and BAI3, novel genes homologous to brain-specific angiogenesis inhibitor 1 (BAI1).

    Cytogenet Cell Gene 1997, 79:103-108. OpenURL

Pre-publication history

The pre-publication history for this paper can be accessed here:

http://www.biomedcentral.com/1471-2407/4/45/prepub