| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Imaging, Diagnosis, Prognosis |
Authors' Affiliations: 1 Department of Urology, Veterans Affairs Medical Center and University of California, San Francisco, California; 2 Department of Urology and Biochemistry, Shimane University, Izumo; and 3 Department of Urology, Kagoshima University, Kagoshima, Japan
Requests for reprints: Rajvir Dahiya, Urology Research Center, Department of Urology, Veterans Affairs Medical Center, University of California, San Francisco, CA 94121. Phone: 415-750-6964; Fax: 415-750-6639; E-mail: rdahiya{at}urol.ucsf.edu.
| Abstract |
|---|
|
|
|---|
Experimental Design: To test our hypothesis, prostate cancer samples (170) and benign prostatic hyperplasia samples (69) were examined by methylation-specific PCR for three genes: adenomatous polyposis coli (APC), glutathione S-transferase pi (GSTP1), and multidrug resistance 1 (MDR1). The methylation status of representative samples was confirmed by bisulfite DNA sequencing analysis. We further investigated whether methylation score (M score) can be used as a diagnostic and staging biomarker for prostate cancer. The M score of each sample was calculated as the sum of the corresponding log hazard ratio coefficients derived from multivariate logistic regression analysis of methylation status of various genes for benign prostatic hyperplasia and prostate cancer. The optimal sensitivity and specificity of the M score for diagnosis and for staging of prostate cancer was determined by receiver-operator characteristic (ROC) curve analysis. A pairwise comparison was employed to test for significance using the area under the ROC curve analysis. For each clinicopathologic finding, the association with prostate-specific antigen (PSA) failure-free probability was determined using Kaplan-Meier curves and a log-rank test was used to determine significance. The relationship between M score and clinicopathologic findings was analyzed by either the Mann-Whitney U test, Kruskal-Wallis test, or the Spearman rank correlation test.
Results: The frequency of positive methylation-specific PCR bands for APC, GSTP1, and MDR1 genes in prostate cancer samples was 64.1%, 54.0%, and 55.3%, respectively. In benign prostatic hyperplasia samples, it was 8.7%, 5.8%, and 11.6%, respectively. There was a significant correlation of M score with high pT category (P < 0.001), high Gleason sum (P < 0.001), high preoperative PSA (P = 0.027), and advanced pathologic features. For all patients, the M score had a sensitivity of 75.9% and a specificity of 84.1% as a diagnostic biomarker using a cutoff value of 1.0. In patients with low or borderline PSA levels (<10.0 ng/mL), the M score was significantly higher in prostate cancers than in benign prostatic hyperplasias (2.635 ± 0.200 and 0.357 ± 0.121, respectively). ROC curve analysis revealed that the M score had a sensitivity of 65.4% and a specificity of 94.2% when 1.0 was used as a cutoff value. For all patients, M score can distinguish organ-confined (
pT2) from locally advanced cancer (
pT3) with a sensitivity of 72.1% and a specificity of 67.8%. Moreover, considering patients with PSA levels of <10 ng/mL, the M score has a sensitivity of 67.1% and a specificity of 85.7%. The ROC curve analysis showed a significant difference between M score and PSA (P = 0.010).
Conclusions: This is the first report demonstrating that M score is a new method for multigene methylation analysis that can serve as a good diagnostic and staging biomarker for prostate cancer.
The ability to predict the outcome or pathologic stage of prostate cancer is a worthwhile goal which will enable physicians to make treatment recommendations for patients with prostate cancer. Nomograms, using a combination of three variables (PSA, biopsy Gleason score, and clinical stage) are currently distributed and used clinically (13, 14). However, the utility of multigene methylation analysis as a pretreatment staging biomarker has not been reported.
We have previously shown that the methylation status of GSTP1 or MDR1 is a good biomarker for detecting prostate cancer and correlates with clinicopathologic features (15, 16). We hypothesize that multigene methylation analysis can be a good diagnostic and staging biomarker prior to treatment. To this end, we did methylation analysis of the APC gene and combined the results with the GSTP1 and MDR1 data from the same 170 prostate cancer and 69 benign prostatic hyperplasia patients. We used the data to calculate a methylation score (M score) that is the sum of the log of the hazard ratios (HR) for each gene, analyzed by multivariate logistic regression analysis for pathology (benign prostatic hyperplasia versus prostate cancer). We also related the M score to clinical and pathologic outcome. Using receiver operator characteristic (ROC) curve analysis to determine the optimal cutoff value for the M score, we evaluated the sensitivity and specificity of M score as a staging biomarker compared with PSA or Gleason sum. In addition, we assessed PSA failure-free probability against the clinicopathologic features of prostate cancer.
| Materials and Methods |
|---|
|
|
|---|
0.2 ng/mL designated as PSA failure. Forty-six patients with prostate cancer were excluded from the PSA failure-free probability study because of adjuvant hormonal therapy immediately after radical prostatectomy. PSA failure-free probability was determined as the percentage of patients without PSA failure. Follow-up ranged from 0.7 to 91.4 months, with a median of 33.9 months. Written informed consent was obtained from all patients.
|
Nucleic acid extraction. Genomic DNA from all prostate samples was extracted using a commercial kit (Qiagen, Valencia, CA), and precipitated with ethanol. The concentration of DNA was determined with a spectrophotometer, and its integrity was checked by gel electrophoresis.
Methylation analysis. Genomic DNA from all prostate samples (100 ng) was subjected to sodium bisulfite modification using a CpGenome DNA Modification Kit (Intergen Co., Purchase, NY). Based on the functional promoter sequence of APC (21), methylation-specific PCR (MSP) and unmethylation-specific PCR (USP) primers were designed using MethPrimer software (http://www.urogene.org/methprimer) developed in our laboratory (22). Primers used for MSP and USP analysis are as follows: universal primers, 5'-TAATTTTTTTGTTTGTTGGGGATT-3' (sense), 5'-ACTACACCAATACAACCACATATC-3' (antisense); MSP primers, 5'-TATTGCGGAGTGCGGGTC-3' (sense), 5'-TCGACGAACTCCCGACGA-3' (antisense); USP primers, 5'-GTGTTTTATTGTGGAGTGTGGGTT-3' (sense), 5'-CCAATCAACAAACTCCCAACAA-3' (antisense). An initial PCR product was created with universal primers, which have no CpG sites in either forward or reverse primers, followed by a second nested PCR with primers specific for MSP or USP. For semiquantitative analysis, a preliminary suitable number of PCR cycles for each MSP and USP were carried out in order to determine the linear range of the reaction. To ensure this, at least one initial PCR was done using 32 cycles each for MSP and USP. Then, a suitable PCR cycle was chosen for each sample. The annealing temperature and PCR cycles used for MSP and USP primers were 64°C and 32 cycles, respectively. The sequences of primers for GSTP1 and MDR1, as well as their PCR conditions, were described previously (15, 16). In each assay, the absence of DNA template served as negative controls. The obtained MSP and USP products were analyzed by electrophoresis in 3% agarose gels and stained with ethidium bromide. With ImageJ software (http://rsb.info.nih.gov/ij), relative methylation levels (%) were calculated (15, 16, 23) by using the area under the curve corresponding to each band (MSP and USP). For methylation analysis of APC, we used 5.3% methylation as a cutoff value, which was the average percentage of methylation in 69 benign prostatic hyperplasia samples. Using these criteria, MSP positivity for APC was defined as those prostate cancers with a percentage of methylation of >5.3%, and negative methylation was <5.3%. The criteria for MSP positivity of GSTP1 and MDR1 were described previously (15, 16).
Bisulfite DNA sequencing analysis. Bisulfite-modified DNA (1 µL) was amplified using a pair of universal primers in a total volume of 20 µL. Direct bisulfite DNA sequencing of the PCR products using either forward universal primer or reverse universal primer was done according to the manufacturer's instructions (Applied BioSystems, Foster City, CA).
Statistical analysis. Using a previously reported analytic technique (24, 25), we calculated the M score for each sample, defined as the sum of the corresponding log HR coefficients, which were derived from multivariate logistic regression analysis of each methylated gene in the benign prostatic hyperplasia and prostate cancer samples (Table 2). The optimal sensitivity and specificity of the M score for diagnosis of prostate cancer and for staging was determined by ROC curve analysis using MedCalc Software (Mariakerke, Belgium). A pairwise comparison was employed to test for significance using the area under the ROC curve (AUC) analysis. For each clinicopathologic finding, the association with PSA failure-free probability was determined using Kaplan-Meier curves and a log-rank test was used to determine significance. The relationship between M score and clinicopathologic findings was blindly analyzed by either the Mann-Whitney U test, Kruskal-Wallis test, or the Spearman rank correlation test using StatView software (SAS Institute Inc., Cary, NC). A P value of <0.05 was regarded as statistically significant.
|
| Results |
|---|
|
|
|---|
|
|
10 ng/mL (2.64 ± 0.20; P = 0.027; Fig. 3C). Moreover, the M score was higher in advancing pathologic features as follows: in capsular invasion (Cap) [positive (3.80 ± 0.23) versus negative (2.47 ± 0.19; P < 0.001)], in seminal vesicle involvement [positive (4.88 ± 0.11) versus negative (2.80 ± 0.16; P = 0.002)], in pelvic lymph node metastasis (pN) [positive (4.60 ± 0.38) versus negative (2.57 ± 0.19; P = 0.001)], in venous involvement (v) [positive (3.92 ± 0.22) versus negative (2.44 ± 0.19; P < 0.001)], in lymphatic vessel involvement (Ly) [positive (3.77 ± 0.28) versus negative (2.61 ± 0.18; P = 0.001)], and in perineural invasion (PNI) [positive (3.41 ± 0.18) versus negative (1.86 ± 0.28; P < 0.001)] (Fig. 3C). We also observed that the M score was age-related in the total group of benign prostatic hyperplasias (P < 0.001) but not in the total group of prostate cancers (P = 0.108).
|
|
pT2) and locally advanced cancer (
pT3), the optimal cutoff value of M score, PSA, and Gleason sum was determined using ROC curve. As shown in Fig. 5A, for all patients with prostate cancer, the M score had a sensitivity of 72.1% and a specificity of 67.8% when 4.0 is employed as a cutoff value. PSA had 69.4% and 69.5%, respectively, with a cutoff value of 9.0 ng/mL; Gleason score had 60.0% and 65.5%, respectively, with a cutoff value of 7. The corresponding AUC for each was 0.721, 0.724, and 0.687, respectively. Pairwise analysis for AUC showed no statistical difference among these three markers. Looking at the patients with PSA levels of 10 ng/mL or less, the M score had a sensitivity of 67.1% and a specificity of 85.7% with 3.3 as a cutoff value. PSA had 75.6% and 38.1%, respectively, with a cutoff value of 7.7 ng/mL; Gleason score had 66.7% and 57.1%, respectively, with a cutoff value of 7. The corresponding AUC for each was 0.780, 0.550, and 0.663, respectively. There was a significant difference between M score and PSA (P = 0.010; Fig. 5B).
|
| Discussion |
|---|
|
|
|---|
PSA is the most sensitive diagnostic biomarker for prostate cancer detection thus far. However, its low specificity has forced unnecessary biopsy of patients in order to exclude prostate cancer. Using various kinds of PSA analysis such as free-PSA, complexed-PSA, or total-PSA, and its combinations% free-PSA or % complexed-PSA, many investigators have struggled to find better methods for prostate cancer detection (2, 3). However, the specificity of these tests is
60% at best with 80% sensitivity (2, 3). To make matters worse, in patients with low or borderline PSA levels, any PSA analysis is a poor diagnostic tool for prostate cancer detection because of much lower specificity (2, 3). In the current study, there was significant difference in M score between benign prostatic hyperplasia and prostate cancer (Fig. 2B, left). For prostate cancer detection, using a cutoff value of 1.0, the M score had a 75.9% sensitivity and 84.1% specificity, which is much higher compared with that reported with PSA (Fig. 2B, right). Moreover, in patients with low or borderline PSA levels (<10.0 ng/mL), the M score had high sensitivity of 65.4% and specificity of 94.2% for prostate cancer detection when 1.0 was used as a cutoff value (Fig. 2C, right). Thus, the M score can be a very useful and improved diagnostic biomarker for prostate cancer detection, even in patients with low or borderline PSA levels. Several investigators have already shown that GSTP1 hypermethylation can be readily detected in bodily fluids such as blood plasma, urine, or ejaculate from patients with prostate cancer (710). Taken together, these findings indicate that the M score is applicable for use as a diagnostic biomarker. Therefore, by using both the M score and PSA for prostate cancer screening, it may be possible to reduce the number of unnecessary biopsies, however, additional study is warranted to verify this hypothesis.
Predicting the probability of the final pathologic stage is a worthwhile goal so that physicians can make appropriate treatment recommendations for patients with prostate cancer. Nomograms using a combination of three variables (serum PSA, biopsy Gleason score and clinical stage) are already distributed and used clinically (13, 14) and have been verified in two studies (26, 27). Using ROC curve analysis, they analyzed staging probability based on the nomograms; however, the AUC reported in the two studies differ (0.787 versus 0.684). One reason for the discrepancy may be related to the fact that the clinical (tumor-node-metastasis) stage was subjectively determined by different individuals, thereby introducing a bias in the results. It would be best if pathologic staging were predicted without subjective variables. Interestingly, in this study, the M score showed significant correlation with worse clinicopathologic features such as higher pT and pN categories, higher Gleason sum, capsular extension, involvement of seminal vesicles, veins and lymphatic vessels, and higher preoperative PSA values (Fig. 3). Maruyama et al. also showed that there was significant correlation between their method for multigene methylation analysis and clinicopathologic features (12). However, in their study, the correlation with pathologic stage was less significant (P = 0.04) compared with our results (P < 0.001; Fig. 3) and no cutoff value was used in their analysis. Using ROC curve analysis, we were able to determine a cutoff value for the M score as a staging biomarker, which enabled us to distinguish organ-confined prostate cancer (
pT2) from locally advanced prostate cancer (
pT3). As shown in Fig. 5A, the M score has the sensitivity and specificity to serve as a good predictive staging biomarker like PSA or Gleason sum. Moreover, in patients with PSA levels of 10 ng/mL or less (Fig. 5B), the M score showed a sensitivity of 67.1% and a high specificity of 85.7% compared with corresponding PSA values of 38.1%. ROC curve analysis showed a significant difference between M score and PSA in this category. Our data also indicates that the M score can be a useful biomarker not only in distinguishing prostate cancer from benign prostatic hyperplasia, but also for predicting the final pathologic stage before medical treatment. Among patients with PSA levels of 10 ng/mL or less, the M score could predict the final pathologic stage more precisely than other biomarkers (Fig. 5B). However, a prospective study will be necessary to confirm this idea. We also analyzed PSA failure-free probability as disease-free survival. Of the clinicopathologic features considered, only Gleason sum was significantly associated with poor outcome (Fig. 4). Currently, there are no other reports demonstrating a significant correlation between multigene methylation analysis and PSA failure-free probability.
In conclusion, this is the first study to integrate the methylation status of multigenes using the M score which reflects the comprehensive methylation status of prostate tissues and is useful as a biomarker for detection and staging of prostate cancer. To elucidate the practical effect of M score in predicting prostate cancer outcomes, it will be necessary to include more genes from prostate tissue biopsy and body fluid samples in the future.
| Footnotes |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Received 3/23/05; revised 5/17/05; accepted 6/ 8/05.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
L. Richiardi, V. Fiano, L. Vizzini, L. De Marco, L. Delsedime, O. Akre, A. G. Tos, and F. Merletti Promoter Methylation in APC, RUNX3, and GSTP1 and Mortality in Prostate Cancer Patients J. Clin. Oncol., July 1, 2009; 27(19): 3161 - 3168. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. G. Bromer, J. Wu, Y. Zhou, and H. S. Taylor Hypermethylation of Homeobox A10 by in Utero Diethylstilbestrol Exposure: An Epigenetic Mechanism for Altered Developmental Programming Endocrinology, July 1, 2009; 150(7): 3376 - 3382. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Sunami, M. Shinozaki, C. S. Higano, R. Wollman, T. B. Dorff, S. J. Tucker, S. R. Martinez, F. R. Singer, and D. S.B. Hoon Multimarker Circulating DNA Assay for Assessing Blood of Prostate Cancer Patients Clin. Chem., March 1, 2009; 55(3): 559 - 567. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Roupret, V. Hupertan, D. R. Yates, J. W.F. Catto, I. Rehman, M. Meuth, S. Ricci, R. Lacave, G. Cancel-Tassin, A. de la Taille, et al. Molecular Detection of Localized Prostate Cancer Using Quantitative Methylation-Specific PCR on Urinary Cells Obtained Following Prostate Massage Clin. Cancer Res., March 15, 2007; 13(6): 1720 - 1725. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Urakami, H. Shiina, H. Enokida, H. Hirata, K. Kawamoto, T. Kawakami, N. Kikuno, Y. Tanaka, S. Majid, M. Nakagawa, et al. Wnt Antagonist Family Genes as Biomarkers for Diagnosis, Staging, and Prognosis of Renal Cell Carcinoma Using Tumor and Serum DNA Clin. Cancer Res., December 1, 2006; 12(23): 6989 - 6997. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Urakami, H. Shiina, H. Enokida, T. Kawakami, K. Kawamoto, H. Hirata, Y. Tanaka, N. Kikuno, M. Nakagawa, M. Igawa, et al. Combination analysis of hypermethylated Wnt-antagonist family genes as a novel epigenetic biomarker panel for bladder cancer detection. Clin. Cancer Res., April 1, 2006; 12(7 Pt 1): 2109 - 2116. [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 |