| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Advances in Brief |
Departments of Medical Oncology [D. T., P. G. F., J. S., P. W. K., M. M., W. R. S.], Radiation Oncology [M. E. L.], and Biostatistical Sciences [C. L.], Dana-Farber Cancer Institute; Departments of Medicine [P. W. K., W. R. S.] and Pathology [M. Lec., M. Lod.], Brigham and Womens Hospital; Departments of Medicine [P. G. F., P. W. K., W. R. S.] and Pathology [M. Lod., M. M.], Harvard Medical School; and Department of Biostatistics [M. Lin, C. L.], Harvard School of Public Health, Boston, Massachusetts 02115
| ABSTRACT |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
Recessive oncogenic alterations typically lead to the biallelic inactivation of tumor suppressor genes (1)
. Alterations such as homozygous and heterozygous deletions or gene conversions are thought to be among the most common genetic abnormalities in epithelial cancers and the detection using LOH5
analysis as being essential to the discovery of the genes targeted by such events (2)
. Oligonucleotide microarrays, capable of simultaneously determining the genotype of 1494 SNPs, have been used to map regions of LOH in small cell lung, breast, bladder, and prostate cancer (3, 4, 5, 6, 7)
. In small cell lung cancer, detection of regions of LOH using SNP arrays was shown to be comparable with LOH detection using microsatellite markers but required cancer cell purity of
90% (3)
. Similarly, LOH events detected by SNP arrays in bladder cancer were consistent with microsatellite-detected events (6)
. Here, we show that using this technology, the genotype of 1494 SNP alleles and genome-wide SNP-based LOH maps can be obtained from laser capture microdissected samples using as little as 5 ng of DNA obtained from paraffin-embedded prostate cancer samples. Furthermore, we demonstrate the utility of a new bioinformatic tool, dChipSNP, that automates the detection of shared regions of LOH, allows for the hierarchical clustering of cancers based on patterns of shared LOH, and allows the analysis of the relationship between clinical parameters and LOH genotypes. Using these methods, we show that subtypes of prostate cancer likely arise through independent genetic pathways.
| Materials and Methods |
|---|
|
|
|---|
DNA Isolation.
Half of each snap-frozen SV sample was fixed in 10% buffered Formalin acetate, mounted in Tissue-Tek OCT, H&E stained, and visualized by light microscopy to confirm the absence of cancer. Germ-line DNA was prepared from the remainder by proteinase K digestion in SDS buffer followed by two extractions in Tris-buffered phenol and chloroform extraction. DNA was precipitated with 1/10th volume of 3.5 M sodium acetate, 20 µg of glycogen, and 2.5 volume of ethanol; washed in 70% ethanol twice; and resuspended in 50 µl of TE.
Cancer DNA was obtained from the corresponding formaldehyde-fixed paraffin-embedded prostatectomy cancer tissue. In brief, cancer epithelial cells were identified by pathologists (M. Loda and M. Lechpammer) and retrieved by LCM using an Arcturus PixCell II system (Mountain View, CA). Before DNA extraction, LCM adhesive caps were visualized by light microscopy to ensure
90% purity of epithelial cells. Cancer cells adhered to the LCM cap were incubated in 50 µl of digestion buffer (1% Tween 20, 1 mg/ml Proteinase K, and 1 x TE) overnight at 37°C, heated to 95°C for 10 min, and clarified by centrifugation. DNA was quantified using the PicoGreen double-stranded DNA Quantitation Kit using 2 µl of DNA (Molecular Probes; P-11495). The sample volume was split in two and subjected to SNP PCR amplification according to the manufacturers protocol without further DNA purification.
HuSNP PCR Amplification.
Normal and cancer DNA was assayed using a modification of the GeneChip HuSNP protocol (Affymetrix). The total amount of input DNA was 120 ng for normal DNA and ranged from 4.4 to 47.1 ng (mean 26.4 ng) for cancer samples. For each genotype, 24 separate pools of primer pairs (
5060 loci/pool) were added to 5 ng of normal DNA (or 0.37 pg to 2 ng of cancer DNA), 5 mM MgCl2, 0.5 mM deoxynucleoside triphosphates, 1.25 units of AmpliTaq Gold (Applied Biosystems), and supplied buffer to a final reaction volume of 12.5 µl/pool. PCR amplification was carried out by denaturing at 95°C for 5 min, followed by 30 cycles of 95°C for 30 s, 58°C for 55 s, 72°C for 30 s, and a final extension of 7 min at 72°C. A 1-µl sample of each PCR reaction pool was diluted separately in 999 µl of H2O. From each of these dilutions, 2.5 µl were removed and added to 0.8 µM biotinylated-T3 and 0.8 µM biotinylated-T7 primers, 4 mM MgCl2, 0.4 mM deoxynucleoside triphosphates, 2.5 units of AmpliTaq Gold, and the supplied buffer for a final reaction volume of 25 µl and reamplified by denaturing at 95°C for 8 min, followed by 40 cycles of 95°C for 30 s, 55°C for 90 s, 72°C for 30 s, and a final extension of 7 min at 72°C. Successful PCR amplification was confirmed by resolving 1.8 µl of each pool in a 3% agarose gel. The 24 pools were combined and concentrated in a Microcon-YM10 centrifugal filter (Millipore Corp.), and the final sample volume was adjusted to 60 µl.
Hybridization, Washing, and Staining of the HU2K Oligonucleotide Array.
Thirty microliters of each sample were diluted in 3 M Tertramethylammouniom chloride, 2 nM oligonucleotide B1 (Affymetrix), 10 mM Tris-HCl (pH 7.8), 0.01% Tween 20, 5 mM EDTA (pH 8.0), 100 µg/ml herring sperm DNA, and 5 x Denhardts solution to a final volume of 135 µl; heated to 95°C for 10 min; and quenched on ice for 5 min. Samples were then hybridized to the GeneChip HuSNP arrays (Affymetrix) overnight at 44°C at 40 rpm. Arrays were washed twice with 6 x SSPE, 0.01% Triton X-100 at 25°C, six times with 4 x SSPE, 0.01% Triton X-100, and stained with 500 µl of 50 µg/ml R-Phycoerythrin Streptavidin (Molecular Probes), 5 µg/ml biotinylated-antistrepavidin antibody (Vector Labs), in 6 x SSPE, 1 x Denhardts solution, and 0.01% Triton X-100 for 30 min at 25°C. The arrays were washed six times with 6 x SSPE, 0.01% Triton X-100 at 25°C.
Scanning of the HuSNP Arrays and Assigning Genotypes.
After staining of the arrays, the chips were scanned using a HP GeneArray Scanner according to the GeneChip HuSNP Mapping Assay Manual (Affymetrix). Genotypes were assigned by the Affymetrix GeneChip 4.0 software. For each SNP, possible assigned genotypes included, homozygous for one allele (AA or BB), heterozygous (AB). If the software could not make a genotype determination, possible assignments included: (a) AB_A (implying the genotype is either AB or A); (b) AB_B (indicating the genotype is either AB or B); or (c) "no signal." The overall call rate was determined by the software to the number of SNPs assigned to AB, AA, or BB divided by the total number of SNPs on the microarrays (1494).
dChipSNP Analysis.
The detailed statistical methods used to define regions of LOH, assign LOH calls to each cancer, and cluster cancers based on LOH are discussed and presented in full in a forthcoming study.6
Briefly, within a 6 Mb window, a summary LOH score was derived across the entire set of cancer/normal comparisons. The summary LOH score was compared with the same score calculated from randomly permuted data. Regions of LOH differing statistically in the actual data from that derived from the permuted data were identified, and LOH assignments were made for each cancer based on the calls made for specific SNP alleles within the bounded region. Hierarchical clustering based on LOH calls for these specific regions was carried out as described previously (9)
.
| Results and Discussion |
|---|
|
|
|---|
The genotypes for 1494 SNPs were obtained using 120 ng of SV as the normal DNA. As described previously, each SNP allele was specifically amplified in 24 pooled highly multiplexed PCR reactions, and after reamplification and labeling, the genotype of each allele was determined by hybridization to the Hu2K SNP array (Affymetrix) using the standard HuSNP protocol (as described in Lindblad-Toh, 2000). The rate of alleles successfully called (call rate) on the array is an indicator of the quality of the DNA extraction, amplification, and hybridization procedures and determined using the MAS software (Affymetrix). The average call rate was 83.2% and ranged from 74.8 to 87.7% for assays run on the SV DNA. These results are consistent with data published previously (Fig. 1)
.
|
To enable the analysis of the SNP array data across multiple cancernormal comparisons, a new informatics platform known as dChipSNP was developed. This platform is based on dChip, a bioinformatics package for oligonucleotide array-based expression analysis (11 , 12) . The details of the statistical methods used in the dChipSNP program for LOH determination are described in a forthcoming study.5 dChipSNP integrates publicly available SNP, gene, and cytoband mapping information with automated methods for detecting statistically meaningful regions of LOH from a series of paired normal and cancer SNP genotypes. dChipSNP is available for download on the Internet.8
Regions of LOH, usually arising as a result of either hemizygous deletion or gene conversion events, are typically defined as stretches of chromosomal areas where all heterozygous and thereby informative alleles are rendered homozygous in the cancer. The boundaries of such regions of LOH are defined by either the presence of retained heterozygous alleles, the ends of chromosomal arms, or the centrosome. This classical definition assumes that all data points are completely accurate and that all polymorphic alleles are mapped correctly within the genome. Array-based methods of SNP detection may have a certain degree of inaccuracy ("noise"), and moreover, the precise genome mapping of each SNP is still not completely stable. Thus, "true" regions of LOH can be interrupted by apparently false positive "retained" SNP alleles. Conversely, true regions of retention of heterozygosity may be interrupted by false LOH calls. To take these concerns into consideration, a method for assessing the probability of LOH within a given set of normal cancer comparisons was used where within a sliding 6 Mb window, a summary LOH score was derived. To determine whether the measured summary LOH scores exceed those that might occur by chance alone, 1000 data sets were generated in which the sample designation of cancer versus normal was randomly permuted. Regions of deletion were identified where the summary LOH score exceeded that found in the permuted data sets. A P-value, corrected for multiple hypothesis testing, of 0.05 was used as the cutoff (as described in detail in Lin et al.). The predicted region of LOH was demarcated by dChipSNP. Fig. 2A
shows the application of this platform to this data set, specifically showing a demarcated region of deletion predicted on chromosome 8 (8p21).
|
Next, we sought to determine whether sample organization algorithms, such as hierarchical clustering, could be applied to this data set. Here, the direct use of the actual SNP allele data in clustering algorithms, much as one would use expression array data, would lead to sample organization based on the germ-line genetic similarities of the individuals rather than on the genetic similarities of the cancer. Moreover, LOH of the same chromosomal region can be marked by completely different sets of heterozygous and thus informative SNPs in distinct cancernormal pairs. Thus, it is necessary to render the SNP data for each cancernormal pair into an LOH-based data set. To this end, for each cancernormal pair, a designation of LOH, retained or uninformative, was made based on the actual genotype of the heterozygous SNPs within the seven regions of significant LOH identified above. Here, retention (RET) was defined as the presence of an informative allele within the region of interest in which heterozygosity was retained (RET) without any heterozygous alleles showing LOH. In dChipSNP, these data can be output as specific tab delimited data sets (see Supplemental Table 1). For purposes of display only, the regions of LOH or retention for each cancernormal pair are projected as a heat map of blue or yellow using a 50% extension to the next nearest heterozygous SNP. In certain instances, the heat map extending from an allele showing LOH extends into the region of LOH; however, if the actual SNP allele is not within the demarcated region, than the cancer is not designated as LOH for that region.
After such designation, one can cluster cancers based only on a single region of LOH in the chromosome view (Fig. 2)
or using all regions of loss defined as significant in the whole genome view. In the first instance, hierarchical clustering was driven by using the comparisons LOH-to-LOH, LOH-to-RET, and RET-to-RET. This resulted in either two or three clusters of tumors for each region of LOH, specifically a cluster with LOH, a cluster with retention, and a cluster of uninformative tumor (Fig. 2A
and Suppl. Fig. 1, AF). This analysis allows the automatic assignment of the genotype for this region to each individual tumor and can be then used in further downstream analysis pertinent to the specific region, e.g., we have looked for gene expression signatures associated with specific regions of LOH by using those tumors that are definitely scored as LOH and those that are definitively scored as RET, while ignoring those that are uninformative in a given region.
For many human cancers, it is thought that there is a sequential progression of accumulated genetic events, ultimately cumulating in a metastatic cancer. In prostate cancer, it remains to be seen whether the heterogeneous nature of the disease can be understood based on this model wherein early genetic events would be then accompanied or not accompanied by late genetic alterations, thus leading to differences in the disease, or whether an alternative, but not mutually exclusive, model in which distinct "parallel" sets of genetic alterations might occur, leading to unique genetic cancer subtypes that have distinct clinical behavior.
To determine whether genome-wide LOH mapping might provide an answer to this question, hierarchical clustering was used to look for separation or cosegregation of LOH events. In this analysis, clustering was enacted using the seven regions of LOH passing the P = 0.05 threshold. Again, each tumor was assigned a lost or retained designation based on SNP alleles within the boundaries of LOH. The comparisons of LOH-to-LOH and LOH-to-RET were used, whereas RET-to-RET was ignored. This allows clustering to be driven primarily by the similarity or difference in a deletion rather than the similarity in retention.
A number of cancers had either no informative alleles or did not show LOH in any of the seven regions and thus were not clustered in this analysis (single branches in Fig. 3
. The clustering of tumors with sufficient information revealed that samples were clustered into distinct branches and that these branches contained tumors enriched for specific nonoverlapping regions of heterozygosity. These data suggest the possibility that there are distinct genetic subsets of prostate cancer that can be defined based on LOH analysis. Validation of this observation on higher density and thus more sensitive SNP arrays is an important next step.
|
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
1 Supported by the Damon-Runyon Lilly Clinical Investigator Award, National Cancer Institute (P50CA09038 and U01CA84995), Association for the Cure of Cancer of the Prostate, and Linda and Arthur Gelb Center for Translational Research. M. M. was supported by the Flight Attendant Medical Research Institute, the Pew Scholars in the Biomedical Sciences, and the National Cancer Institute (U01CA84995). ![]()
2 Supplementary data for this article are available at Cancer Research Online (http://cancerres.aaccrjournals.org). ![]()
3 M. E. L. and M. L. contributed equally to this work. ![]()
4 To whom requests for reprints should be addressed, Dana-Farber Cancer Institute, Boston, MA 02115. Phone: (617) 632-4750; Fax: (617) 632-5417; E-mail: william_sellers{at}dfci.harvard.edu ![]()
5 The abbreviations used are: LOH, loss of heterozygosity; LCM, laser capture microdissection; SSPE, saline-sodium phosphate-EDTA; SNP, single nucleotide polymorphic allele; SV, seminal vesicle. ![]()
6 M. Lin et al., submitted for publication. ![]()
7 Internet address: http://research.dfci.harvard.edu/sellerslab/datasets/index.html. ![]()
8 Internet address: http://www.dchip.org/. ![]()
Received 6/11/03. Accepted 6/30/03.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
G. Neale, X. Su, C. L. Morton, D. Phelps, R. Gorlick, R. B. Lock, C. P. Reynolds, J. M. Maris, H. S. Friedman, J. Dome, et al. Molecular Characterization of the Pediatric Preclinical Testing Panel Clin. Cancer Res., July 15, 2008; 14(14): 4572 - 4583. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. W. Ross, P. D. Ouillette, C. M. Saddler, K. A. Shedden, and S. N. Malek Comprehensive Analysis of Copy Number and Allele Status Identifies Multiple Chromosome Defects Underlying Follicular Lymphoma Pathogenesis Clin. Cancer Res., August 15, 2007; 13(16): 4777 - 4785. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Fuchs, C. Nilsson, J. Kachergus, M. Munz, E.-M Larsson, B. Schule, J. W. Langston, F. A. Middleton, O. A. Ross, M. Hulihan, et al. Phenotypic variation in a large Swedish pedigree due to SNCA duplication and triplication Neurology, March 20, 2007; 68(12): 916 - 922. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Oosting, E. H. Lips, R. van Eijk, P. H.C. Eilers, K. Szuhai, C. Wijmenga, H. Morreau, and T. van Wezel High-resolution copy number analysis of paraffin-embedded archival tissue using SNP BeadArrays Genome Res., March 1, 2007; 17(3): 368 - 376. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Du, J. S. Rozowsky, J. O. Korbel, Z. D. Zhang, T. E. Royce, M. H. Schultz, M. Snyder, and M. Gerstein A supervised hidden markov model framework for efficiently segmenting tiling array data in transcriptional and chIP-chip experiments: systematically incorporating validated biological knowledge Bioinformatics, December 15, 2006; 22(24): 3016 - 3024. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. K. Thomas, B. Weir, and M. Meyerson Genomic approaches to lung cancer. Clin. Cancer Res., July 15, 2006; 12(14): 4384s - 4391s. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Ouyang, T. L. DeWeese, W. G. Nelson, and C. Abate-Shen Loss-of-Function of Nkx3.1 Promotes Increased Oxidative Damage in Prostate Carcinogenesis Cancer Res., August 1, 2005; 65(15): 6773 - 6779. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Nannya, M. Sanada, K. Nakazaki, N. Hosoya, L. Wang, A. Hangaishi, M. Kurokawa, S. Chiba, D. K. Bailey, G. C. Kennedy, et al. A Robust Algorithm for Copy Number Detection Using High-Density Oligonucleotide Single Nucleotide Polymorphism Genotyping Arrays Cancer Res., July 15, 2005; 65(14): 6071 - 6079. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Zhao, B. A. Weir, T. LaFramboise, M. Lin, R. Beroukhim, L. Garraway, J. Beheshti, J. C. Lee, K. Naoki, W. G. Richards, et al. Homozygous Deletions and Chromosome Amplifications in Human Lung Carcinomas Revealed by Single Nucleotide Polymorphism Array Analysis Cancer Res., July 1, 2005; 65(13): 5561 - 5570. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Rennert, C. M. Zeigler-Johnson, K. Addya, M. J. Finley, A. H. Walker, E. Spangler, D. G.B. Leonard, A. Wein, S. B. Malkowicz, and T. R. Rebbeck Association of Susceptibility Alleles in ELAC2/HPC2, RNASEL/HPC1, and MSR1 with Prostate Cancer Severity in European American and African American Men Cancer Epidemiol. Biomarkers Prev., April 1, 2005; 14(4): 949 - 957. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Halin, P. Wikstrom, S. H. Rudolfsson, P. Stattin, J. A. Doll, S. E. Crawford, and A. Bergh Decreased Pigment Epithelium-Derived Factor Is Associated with Metastatic Phenotype in Human and Rat Prostate Tumors Cancer Res., August 15, 2004; 64(16): 5664 - 5671. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. A. Rubin, S. Varambally, R. Beroukhim, S. A. Tomlins, D. R. Rhodes, P. L. Paris, M. D. Hofer, M. Storz-Schweizer, R. Kuefer, J. A. Fletcher, et al. Overexpression, Amplification, and Androgen Regulation of TPD52 in Prostate Cancer Cancer Res., June 1, 2004; 64(11): 3814 - 3822. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. G. Paez, M. Lin, R. Beroukhim, J. C. Lee, X. Zhao, D. J. Richter, S. Gabriel, P. Herman, H. Sasaki, D. Altshuler, et al. Genome coverage and sequence fidelity of {phi}29 polymerase-based multiple strand displacement whole genome amplification Nucleic Acids Res., May 18, 2004; 32(9): e71 - e71. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Zhao, C. Li, J. G. Paez, K. Chin, P. A. Janne, T.-H. Chen, L. Girard, J. Minna, D. Christiani, C. Leo, et al. An Integrated View of Copy Number and Allelic Alterations in the Cancer Genome Using Single Nucleotide Polymorphism Arrays Cancer Res., May 1, 2004; 64(9): 3060 - 3071. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |