|Home | About | Journals | Submit | Contact Us | Français|
The current focus on biomarker discovery is a result of an improved understanding of the biological basis for carcinogenesis and advances in technology. Biomarkers can aid in diagnosis, prognosis, treatment selection, and drug development. There is an urgent need for high-resolution tools that perform well using archived tissue for biomarker discovery and tools that can translate into the clinic.
Oligonucleotide array comparative genomic hybridization (oCGH) was compared to BAC-based aCGH using unamplified total genomic DNA from formalin fixed paraffin-embedded (FFPE) prostate tissue.
The copy number aberrations detected with the BAC and oligonucleotide arrays were highly correlated in cases where the arrays contained probes in similar genomic locations. The oligonucleotide array platform provided more precise mapping due to the higher density of oligonucleotide probes.
These results demonstrate the utility of high-resolution oligonucleotide arrays designed to use genomic DNA for CGH measurements using archived tissue samples for discovery and clinic based assays.
Array comparative genomic hybridization (aCGH) is a valuable tool for identifying DNA copy number changes in tumor genomes. Genomic copy number alterations can lead to altered expression of oncogenes and tumor suppressor genes. Moreover, copy number abnormalities can be associated with the clinical course of the disease and used for selection of therapy. Detailed mapping of amplicons and deletions localizes potential therapeutic targets. Formalin-fixed paraffin-embedded (FFPE) tumor specimens provide a rich source of patient samples for these studies, and often with corresponding long follow-up data. FFPE samples yield DNA that is often degraded and therefore more challenging to use for aCGH. Our published DNA extraction methods, utilized in this manuscript, allows for routine processing of FFPE tissue from a variety of sources . The use of FFPE prostate tissue on BAC-based arrays has been demonstrated previously by our laboratory . Recently, oligonucleotide array based CGH platforms have emerged with the potential for more flexible array design and higher-resolution copy number mapping [2–4] thus, significantly increasing the power of studies aimed at discovery of new therapeutic and diagnostic targets. Combined with biomarkers of aggressive disease, such technologies could translate directly into the clinic. This is of particular relevance to prostate cancer where over detection and over treatment are significant despite a protracted and non-life threatening natural history in many men . Therefore, it is important to couple biomarkers and technologies that will work well with FFPE biopsy specimens. We report on the use of genomic DNA from FFPE prostate tumors using the Agilent oligonucleotide CGH (oCGH) platform and the comparison of these results to those obtained using UCSF scanning BAC arrays. Application of Agilent arrays for analysis of FFPE biopsy specimens is also demonstrated.
DNA from four different radical prostatectomy cases (designated 13P, 33P, 34P, and 41P) was isolated from FFPE tissue. A single pathologist outlined areas of greater than 75% tumor for macrodissection with a scalpel. DNA was extracted from the tissue scrapings using the Puregene DNA isolation kit (Gentra) as per the manufacturer’s instructions. Two phenol:-chloroform extractions followed by an ethanol precipitation were performed after the Gentra kit’s final elution step. The biopsy specimen from a routine fine needle, biopsy was embedded in an FFPE block. An H&E guide slide was used to macrodissect tumor tissue (>90%) from 10 slides of 10 μm thickness. The tissue was digested using proteinase K for 48 hr at 56μC. Genomic DNA was isolated using the Qiagen (Valencia, CA) QIAamp DNA Micropurification kit using the manufacturer’s protocol for biopsy specimens. This was followed by an ethanol precipitation. DNA was quantitated using the Nanodrop spectrophotometer. DNA quality was assessed by the 260:280 ratio and its integrity by agarose gel ethidium bromide visualization. The DNAs were found to be of typical quality in terms of purity and the amount of degradation observable by gel electrophoresis.
aCGH was performed using BAC arrays containing 2,460 BAC clones printed at UCSF as well as Agilent Human Genome 44 or 244 K 60 mer oligonucleotide arrays containing approximately 40,000 probes with an average spatial resolution of ~35 kb or 244,000 probes with an average resolution of ~9 kb. The BAC aCGH was performed as described in Paris et al.  with a male reference DNA (Promega). The standard oligonucleotide oCGH experiments were performed as per the manufacturer’s instructions (http://chem.agilent.com/scripts/LiteraturePDF.asp?iWHID=39980) with the following exceptions: the input DNA was lessened from 2 to ~1 μg, the hybridization time was increased from 40 to 65 hr and wash 3 was not performed. Prostate samples 13P and 33P were hybridized with female reference DNA (Promega) in dye flip pairs to provide an additional source of confidence in copy number calls. Additional oligonucleotide array experiments were performed for sample 13P, 34P, and 41P and the biopsy using only 500 ng DNA to ascertain the utility of working with low genomic yield samples. Components of labeling and hybridization were identical to standard reactions, except the DNA was digested with restriction enzymes Alu1 and Rsa1 for only 2 hr followed immediately (without cleanup) by fluorescent labeling for 1 hr, and hybridization was carried out for 40 hr. The 500 ng prostate sample was hybridized with male reference DNA. Agilent Feature Extraction software version 188.8.131.52 was used to extract feature level data from the Agilent Microarray Scanner files.
Regions of copy number gain and loss for the BAC aCGH data were identified by creating sample specific thresholds [6,7]. The clones with log2 ratios above or below +/− a tumor sample’s threshold value were considered as gains or losses, respectively. The aberration detection module-1 (ADM-1) aberration detection algorithm from Agilent’s CGH Analytics software was used to identify regions of copy number gain or loss from both the oligonucleotide and BAC CGH data . Briefly, ADM-1 identifies aberrant genomic intervals based on a statistical score. This score is calculated as the average log2 ratio of probes in the interval, multiplied by the square root of the number of such probes and divided by the derivative log2 ratio spread of the array. Aberration calls were made for each experiment after initial preprocessing of the data that included combining log2 ratios for replicate probes and centralization. In the centralization step, all log2 ratios are shifted by an array-specific constant, such that ADM-1 applied to shifted log2 ratios calls the minimum number of probes aberrant. The centralization step was applied only to the oligonucleotide data.
aCGH and oCGH were performed with DNA isolated from FFPE preserved samples from four different radical prostatectomy cases (designated 13P, 33P, 34P, 41P) using both BAC arrays containing 2,460 BAC clones printed at UCSF as well as Agilent’s Human Genome 44 or 244 K 60 mer oligonucleotide arrays containing approximately 40,000 and 244,000 probes, respectively. Samples were selected to represent both good (≥12 kb) and poor quality (≥500 bp) DNA that is obtained from FFPE tissue. All samples were archived for 8 years prior to processing, except sample 41P which was 6 years old.
Detailed comparison of aberration calls made from the 44 K oligonucleotide and BAC array data was performed for FFPE samples 13P and 33P, with the focus being regions where both platforms had overlapping probes. Dye flip experiments for the oligonucleotide data were combined, and aberrations were called using CGH Analytics aberration detection algorithm ADM-1 with a threshold of 6. Aberrations in BAC data were identified using the standard sample specific threshold method [6,7] and also using ADM-1. Aberrant intervals identified in the BAC data showed good agreement with corresponding intervals from the oligonucleotide data (Fig. 1 and Table I). We also computed the fraction of BAC probes above or below the sample specific threshold showing gain or loss that were inside oligonucleotide amplified or deleted intervals. For samples 13P and 33P, 82 and 90%, respectively, of BAC probes showing gain or loss were identified within oligonucleotide aberrant intervals.
There was good concordance observed between the overall copy number profiles across the genome obtained from the BAC and oligonucleotide array platforms for both samples 13P and 33P both in terms of the genomic position of the gains and losses and in the magnitude of the copy number differences (Fig. 1). The Pearson correlation of the average log2 ratios in matching aberrant regions was 0.87 and 0.96 for samples 13P and 33P, respectively. For example, deletion of the entire p arm on chromosome 8 for sample 13P was identified by an average log2 ratio in both the BAC and oligonucleotide data as 0.56. This deletion call was based on 59 BAC probes and 556 oligonucleotide probes, extending from 0 to 38 Mb in the BAC data and 0 to 43 Mb in the oligonucleotide data. A detailed visualization of the aberrations observed on chromosome 8 is highlighted in Figure 2A. Chromosome 8 was chosen because 8p is known to be commonly deleted in prostate cancer [9–11].
In some cases, the oligonucleotide array platform provided more precise mapping of aberration boundaries due to the higher density of oligonucleotide probes. The higher density of oligonucleotide probes can also add statistical confidence to copy number calls, especially where only one BAC probe maps to an aberration. An example highlighting this was found in sample 13P on chromosome 12 (p 12.1) where the same copy number change was observed in both the BAC and oligonucleotide array data. However, the aberration breakpoints were more precisely mapped in the oligonucleotide data due to the higher density of probes (Fig. 2B), thereby narrowing the number of candidate genes.
Oligonucleotide arrays developed specifically for CGH [2–4] and whole genome BAC tiling arrays  have the potential to provide very high-resolution copy number measurements. The goal of these experiments was to assess the quality of CGH data obtained with whole genome oligonucleotide arrays using clinically relevant, but challenging DNA from archived prostate tissue, that has been shown previously to work well using BAC arrays. The results from the oligonucleotide array platform correlated very well with the BAC array results although there were some instances of aberrations detected more robustly and with more precise mapping with the oligonucleotide arrays due to the higher probe density on these arrays. The only considerable differences were attributable to regions where clones were absent in the BAC data. It is note worthy that Agilent has selected probes biased toward genes, in particular cancer related genes (represented by a minimum of two probes), ensuring adequate coverage in the most commonly studied genomic regions.
In our hands, DNA extracted from frozen or FFPE tissue using our published protocols are equivalent for BAC-based aCGH and this is true regardless of the laboratory of source’s fixation protocol [1,7,13]. This is consistent with a study published by Little et al. comparing frozen and FFPE DNA for CGH on BAC arrays. In the present study matched fresh frozen and fixed specimens could not be directly compared. To overcome this limitation, we embedded DUI45 prostate cancer cells to mimic routine frozen and FFPE archiving and compared extracted DNA on 244 K oCGH arrays to DNA extracted from unfixed DU145 cells on BAC arrays. The frozen and FFPE DNA produced concordant copy number profiles on the Agilent arrays and produced copy number profiles essentially identical to each other (Fig. 3). In addition, these profiles match our unpublished and others published BAC aCGH data for DU145 [14,15]. We chose to focus on FFPE material in this manuscript because it is commonly believed to be more difficult to work with than frozen material and because of its importance for translational research.
DNA yields from FFPE specimens may be small because many of the most informative experiments and clinical applications will need to begin with needle biopsies where yields may be in the range of 500 ng. Thus, it is significant that we obtained comparable oCGH results using only 500 ng of FFPE (Fig. 4). Figure 4A shows an overlay of BAC-based aCGH and oCGH using 1 μg DNA and oCGH using 500 ng DNA from sample 13P (Fig. 4A). Qualitative assessment of two additional oCGH 500 ng samples in Figure 4 shows great similarity with the corresponding aCGH data using 1 μg of DNA. The average log2 ratio standard deviation of the replicate probes randomly dispersed on the array, which serve as a measure of the quality of the array result, was 0.028 (34P) and 0.071 (41P). Next we extracted DNA from an FFPE prostate tumor biopsy and for analysis on the 244K oCGH platform. The average standard deviation of the log2 ratios for the replicate probes on the biopsy array was 0.039. A penetrance plot is shown for the copy number changes detected by oCGH for the biopsy and its matched primary tumor (Fig. 5). It may of interest to note the similarity between the two copy number plots despite the fact that distinct foci of the same tumor were analyzed.
To the best of our knowledge, this is the first report on the use of commercial oligonucleotide arrays with unamplified FFPE prostate biopsy and radical prostatectomy tissue. Van den Ijssel et al.  recently published a paper reporting the use of their in-house spotted oligonucleotide arrays with DNA from a single FFPE, gastric tumor. These authors report qualitatively similar results with their oligonucleotide platform and BAC arrays. In our report, we have conducted qualitative and quantitative comparisons and have used a commercially available oligonucleotide platform. This latter point is particularly relevant for future clinical applications since the FDA just recently approved an assay system that included the Agilent microarray (Agendia MammaPrint).
The use of archival tissue is very important for prostate cancer research due to its long protracted natural history. The value of PSA screening has long been a matter of debate. A recent study found that men who have been screened for prostate cancer by the most commonly used tests (PSA and DRE) have no greater chance of surviving the disease than those who have not been screened at all . This highlights the need for the identification of new prognostic and predictive biomarkers that will compliment, and therefore improve upon existing clinical standards to help physicians and their patients make decisions regarding treatment. To the best of our knowledge this is the first time oCGH copy number data has been obtained using unamplified DNA from a FFPE tumor biopsy and its matched primary tumor. This should encourage utilization of biopsy specimens for basic and translational studies. oCGH’s higher resolution will allow for better characterization of tumor genomes and therefore expedite the identification of genes driving progression and candidate therapeutic targets. The DNA extraction protocol followed in this paper is an essential step to obtaining DNA from FFPE tissue that is useable for aCGH and oCGH. Although we cannot generalize to all cancer tissue types, this DNA protocol has allowed us to obtain copy number profiles from DNA extracted from FFPE tissue obtained from multiple institutions around the world, regardless of the age of the sample (up to 16 years old). (1, 7, and unpublished data). Based on our experience, the UV–Vis values (260/280 ~1.8, 260/230 ~2) and DNA integrity (>500 bp) visualized by gel chromatography can determine whether a sample will perform well on a CGH array. A sample seems to fail due to the case, rather than the archival method (e.g., frozen vs. FFPE). The Agilent oligonucleotide platform and methods used in these studies provide high quality, reproducible oCGH copy number profiles from small amounts of FFPE extracted DNA and therefore, represents a valuable tool for biomarker discovery and possibly development of clinical assays.
Grant sponsor: UCSF Prostate Cancer SPORE, NIH; Grant number: P50CA89520.
This work was supported by the UCSF Prostate Cancer SPORE, NIH Grant P50CA89520. We would like to thank Dr. Peter Carroll, Dr. Jeffry Simko, and Dr. Giancarlo Albo for assistance with sample acquisition, pathology and DNA extraction, respectively. We are very appreciative of the programming assistance from Emmanuel Yera. In addition, the authors also thank Amir Ben-Dor and Zohar Yakhini from Agilent Laboratories for helpful discussions on data analysis.