Abstract
Background: Early prediction of disease progression in men with very low-risk (VLR) prostate cancer who selected active surveillance (AS) rather than immediate treatment could reduce morbidity associated with overtreatment.
Methods: We evaluated the association of six biomarkers [Periostin, (−5, −7) proPSA, CACNA1D, HER2/neu, EZH2, and Ki-67] with different Gleason scores and biochemical recurrence (BCR) on prostate cancer TMAs of 80 radical prostatectomy (RP) cases. Multiplex tissue immunoblotting (MTI) was used to assess these biomarkers in cancer and adjacent benign areas of 5 μm sections. Multivariate logistic regression (MLR) was applied to model our results.
Results: In the RP cases, CACNA1D, HER2/neu, and Periostin expression were significantly correlated with aggressive phenotype in cancer areas. An MLR model in the cancer area yielded a ROC-AUC = 0.98, whereas in cancer-adjacent benign areas, yielded a ROC-AUC = 0.94. CACNA1D and HER2/neu expression combined with Gleason score in a MLR model yielded a ROC-AUC = 0.79 for BCR prediction. In the small biopsies from an AS cohort of 61 VLR cases, an MLR model for prediction of progressors at diagnosis retained (−5, −7) proPSA and CACNA1D, yielding a ROC-AUC of 0.78, which was improved to 0.82 after adding tPSA into the model.
Conclusions: The molecular profile of biomarkers is capable of accurately predicting aggressive prostate cancer on retrospective RP cases and identifying potential aggressive prostate cancer requiring immediate treatment on the AS diagnostic biopsy but limited in BCR prediction.
Impact: Comprehensive profiling of biomarkers using MTI predicts prostate cancer aggressive phenotype in RP and AS biopsies. Cancer Epidemiol Biomarkers Prev; 24(12); 1864–72. ©2015 AACR.
Introduction
Prostate cancer is the most common cancer in men in the United States and is the second most common cause of cancer death in men. Prostate cancer usually occurs after age 50 and the incidence increases with age. There are 233,000 estimated new cases and 29,480 estimated prostate cancer deaths in 2014 in the United States (1). Radical prostatectomy (RP) is an effective treatment for patients with organ-confined disease and has been demonstrated to reduce the risk of death from prostate cancer (2). Nearly 40% of prostate cancer patients who choose definitive therapy will undergo RP. In 38% to 52% of cases, advanced disease with potentially bad prognosis are found in surgical specimens (3). Each category of extraprostatic disease is associated with significantly increased risk of cancer recurrence and progression, measured at the earliest time with a detectable prostate-specific antigen (PSA; >0.20 ng/mL), or biochemical recurrence (BCR; ref. 4). The natural history of prostate cancer progression after BCR following surgery can be highly variable (i.e., 3–13 years); however, at least two thirds of BCR patients develop disease spread if left untreated and many will die because of distant progression (5).
In most cases, BCR is used as a surrogate measure for more clinically meaningful critical endpoints such as distant progression or cancer-specific mortality. Inaccurate risk classifications could result in inappropriate or unnecessary treatments. Unfortunately, existing tools (e.g., nomograms like Shariat, ref. 6; Swanson, ref. 7; and Capra, ref. 8) that rely solely on clinical variables such as PSA velocity, grade, and stage are unable to predict which men will go on to metastasis and ultimately prostate cancer-specific death. These tools may be unable to identify those men that have already demonstrated BCR. Therefore, a better method to predict whether a prostate cancer patient has a more aggressive phenotype at surgery will help to standardize treatments of the more aggressive cancer patients early and efficiently and spare patients of less aggressive cancer from the morbidity associated with adjuvant therapy.
PSA screening and prostate biopsies have resulted in overdiagnosis and overtreatment of prostate cancer (9–11). In fact, low-grade, low-stage prostate cancer of up to approximately 56% men remained undetected during their lifetime (9–12). In the present era, 30 men require treatment for prostate cancer in order to save one man (13). Such overtreatment can always create a chance of a decreased quality of life once sexual function and urinary function are compromised (14). The accurate identification of men with prostate cancer that are destined to progress to life threatening disease and who will benefit from curative intervention must be a high research priority, with the goal of reducing unnecessary treatments.
In 1995, Dr. Ballentine Carter started an active surveillance (AS) program with the delayed surgical intervention as a treatment option in Johns Hopkins Medical Institutions (JHMI, Baltimore, MD). Patients are enrolled if they meet Epstein inclusion very low-risk (VLR) criteria (15–17). These criteria include PSA density (PSAD) <0.15 ng/mL/g and favorable diagnostic needle biopsy characteristics (i.e., Gleason score ≤6, ≤2 cores involved with cancer, ≤50% of any core involved with cancer). Applications of these criteria for selection have been validated to identify men with indolent prostate cancer (15, 18–20). The current JHMI AS patients are followed-up semiannually with assessment of serum total PSA (tPSA), free PSA (fPSA), digital rectal examination (DRE), and annually with 12 core surveillance biopsy to decide requirement for immediate prostate cancer intervention under AS program (17). Furthermore, according to the published experience with AS from different groups, about one third of AS patients will require definitive treatment primarily due to upgrading overtime (21), and we cannot determine at diagnosis which patients will eventually fail from the program. Recently, there is a paradigm shift that stems from recent discoveries identifying important clinically and pathologic predictive pretreatment parameters that include an expanded AS criteria, including from very low- to low-risk prostate cancer (17, 22). In fact, the National Comprehensive Cancer Network (NCCN) now recognizes four separate risk groups that include VLR, low risk, intermediate risk, and high risk of newly diagnosed patients with prostate cancer (23).
This study evaluated six tissue molecular biomarkers (MB) to predict at entry into the AS program the likelihood of failure. The biomarkers evaluated include stromal protein (Periostin), membranous channels and transmembrane proteins (CACNA1D, HER2/neu), nuclear proliferating factor Ki-67, epigenetic regulators (EZH2), and prostate-specific protein (−5, −7) proPSA. They were chosen based on their location at cellular and tissue level, function, and their importance in previous biomarker research (see Discussion). We hypothesized that a combination of these molecular biomarkers could predict an aggressive prostate cancer phenotype in RP specimens. To test this hypothesis, we evaluated our biomarker panel using two JHMI TMA made up of 80 cases of prostate cancer stratified by Gleason score, 22 of which have BCR. A novel multiplex tissue immunoblotting (MTI) technology was applied to quantify the expression of these six biomarkers on single 5 μm TMA tissue sections of RP or biopsy specimens. Also, the biomarker research panel was applied to an AS cohort with VLR prostate cancer to see whether we could predict those men who failed quickly with unexpected aggressive outcomes (progressors). The overall objective of the study was to develop an integrated, quantitative molecular biomarker-based predictor for the early detection of a clinically aggressive prostate cancer based on biopsy specimens in men being considered for AS for VLR or LR prostate cancer.
Materials and Methods
Construction of TMA681 and TMA682
The two TMAs (TMA681 and TMA682) were prepared using a Beecher MT1 manual arrayer (Beecher Instruments) under the supervision of the Prostate Cancer Biorepository Network (PCBN). The formalin-fixed, paraffin embedded (FFPE) RP prostate cancer tissues and normal (benign) cancer-adjacent controls to be included in the two tissue TMAs were selected and reviewed by a JHMI pathologist (J.I. Epstein). Hematoxylin and eosin (H&E)-stained slides from all selected cases were reviewed and “mapped” by the pathologist: the adjacent normal-appearing along with staged and/or graded index tumor areas were identified and marked on the H&E slide for each case. Using these marked “template slides,” the tissue blocks were coordinately marked, and 0.60-mm diameter cores were punched from the normal-appearing and prostate tumor areas and then transferred to recipient blocks. Sections (5 μm thick) of these TMAs were cut and used for downstream MTI studies.
The two TMAs included 80 unique prostate cancer patients representing different Gleason scores (3+3, 3+4, 4+3, and ≥8) with quadruplicates of cancer and cancer-adjacent benign areas. Among these 80 cases, 22 had BCR based on PSA exceeding 0.2 ng/mL on follow-up. The average follow-up time after BCR is 6.6 years (range, 1–14 years). Among the 22 BCR cases, 12 cases had PSA increases only, one case had local recurrence, four cases had distant metastasis, three cases had both local recurrence and distant metastasis, and two cases received adjuvant radiation treatment and the PSA level decreased to < 0.2 ng/mL after radiation treatment. The detailed demographics of the total 80 cases stratified by Gleason scores of which 22 cases had BCR are shown in Table 1 and further subgrouped based on clinical features. In the 80 cases, one patient in TMA 682 died from noncancer cause; there is no follow-up data in five cases of TMA 681 and five cases of TMA 682.
. | . | Recurrence . | . | |
---|---|---|---|---|
Variable . | Subgroups . | No . | Yes . | Total . |
Age | 59.04 ± 6.57 | 58.82 ± 5.80 | ||
Gleason score | 6 | 7 | 1 | 8 |
7 (3+4) | 14 | 3 | 17 | |
7 (4+3) | 15 | 3 | 18 | |
≥ 8 | 11 | 15 | 26 | |
TNM stage | T2 | 26 | 2 | 28 |
T3A | 16 | 15 | 31 | |
T3B | 3 | 5 | 8 | |
Gland weight (g) | 51.94 ± 14.90 | 56.19 ± 21.38 | ||
PSA (ng/mL) | 8.29 ± 5.29 | 13.17 ± 10.15 | ||
PSAD (ng/mL/g) | 0.17 ± 0.11 | 0.25 ± 0.20 | ||
Exposure (y) | 5.06 ± 4.50 | 3.82 ± 2.36 | ||
Race | A | 6 | 3 | 9 |
H | 1 | 0 | 1 | |
O | 3 | 1 | 4 | |
C | 37 | 18 | 55 | |
Surgical margin | Negative | 41 | 16 | 57 |
Positive | 6 | 6 | 12 | |
Seminal vesicle status | Negative | 46 | 16 | 62 |
Positive | 6 | 6 | 12 | |
Lymph node status | Negative | 46 | 22 | 68 |
Positive | 1 | 0 | 1 | |
Capsular penetration (fcp) | Negative | 38 | 16 | 54 |
Positive | 8 | 6 | 14 | |
Capsular penetration (ecp) | Negative | 34 | 8 | 42 |
Positive | 13 | 14 | 27 | |
Organ confinement status | Nonconfined | 21 | 20 | 41 |
Confined | 26 | 2 | 28 |
. | . | Recurrence . | . | |
---|---|---|---|---|
Variable . | Subgroups . | No . | Yes . | Total . |
Age | 59.04 ± 6.57 | 58.82 ± 5.80 | ||
Gleason score | 6 | 7 | 1 | 8 |
7 (3+4) | 14 | 3 | 17 | |
7 (4+3) | 15 | 3 | 18 | |
≥ 8 | 11 | 15 | 26 | |
TNM stage | T2 | 26 | 2 | 28 |
T3A | 16 | 15 | 31 | |
T3B | 3 | 5 | 8 | |
Gland weight (g) | 51.94 ± 14.90 | 56.19 ± 21.38 | ||
PSA (ng/mL) | 8.29 ± 5.29 | 13.17 ± 10.15 | ||
PSAD (ng/mL/g) | 0.17 ± 0.11 | 0.25 ± 0.20 | ||
Exposure (y) | 5.06 ± 4.50 | 3.82 ± 2.36 | ||
Race | A | 6 | 3 | 9 |
H | 1 | 0 | 1 | |
O | 3 | 1 | 4 | |
C | 37 | 18 | 55 | |
Surgical margin | Negative | 41 | 16 | 57 |
Positive | 6 | 6 | 12 | |
Seminal vesicle status | Negative | 46 | 16 | 62 |
Positive | 6 | 6 | 12 | |
Lymph node status | Negative | 46 | 22 | 68 |
Positive | 1 | 0 | 1 | |
Capsular penetration (fcp) | Negative | 38 | 16 | 54 |
Positive | 8 | 6 | 14 | |
Capsular penetration (ecp) | Negative | 34 | 8 | 42 |
Positive | 13 | 14 | 27 | |
Organ confinement status | Nonconfined | 21 | 20 | 41 |
Confined | 26 | 2 | 28 |
NOTE: Age, gland weight, PSA, PSAD, exposure data were shown as mean ± SD. For race data, A, African American; H, Hispanic; O, Others; C, Caucasian. PSAD (ng/mL/g) = PSA (ng/mL)/gland weight (g).
Test TMAs for IHC optimization
A separate TMA (TMA 475), which included prostate cancer and normal prostate tissue, was kindly provided by PCBN and was used to optimize dilutions of antibodies of our biomarkers before applying the antibodies to TMA 681 and TMA 682.
Active surveillance biopsies
FFPE biopsies from a total of 61 patients (29 progressors and 32 nonprogressors) enrolled in the AS program at JHMI were used for MTI study. Progressors were defined as men who met all the AS criteria at entry but were found to have a tumor increased in volume and/or Gleason score >6 on follow-up biopsy/biopsies or retropubic radical prostatectomy criteria were met during AS. For each patient, continuous 5 μm tissue sections were cut and used either for H&E staining and pathologic assessment or for profiling biomarkers using MTI. The H&E slides were first reviewed by a pathologist and the cancer area were marked for reference of data analysis later on. Only slides with cancer and adjacent to the corresponding read H&E slides were used for MTI study.
All cases used for TMAs and AS biopsies were consented under an Institutional Review Board-approved protocol at Johns Hopkins University School of Medicine.
Multiplex tissue immunoblotting
MTI is a method that can be used to detect multiple molecular targets in FFPE tissues, while retaining both quantitative and histomorphologic diagnostic and prognostic features (24–26). Starting with a 5-μm FFPE tissue section on a standard glass slide, proteins were transferred from the tissue sections (TMA or biopsy) onto a series of overlapped thin membranes (P-Film, 20/20 GeneSystems, Inc.). Each membrane was probed with one of the six biomarkers [CACNA1D, Periostin, HER2/neu, EZH2, Ki67 and (−5, −7 proPSA)] followed by incubation with FITC-conjugated secondary antibodies. Total proteins collected on the blotted membranes were biotinylated, and followed by incubation with streptavidin-linked Cy5. The fluorescence signals of biomarkers and total proteins were acquired using a Typhoon 9410 imager (GE Healthcare) and quantified with ImageQuant5.2 software (GE Healthcare). The quantified signals of biomarkers were divided by that of total protein for normalization. The log-transformed ratio was used for downstream data analysis. Supplementary Table S1 provides the detailed information about the primary antibodies and dilutions used in the study. The quality of the antibodies was checked using Western blot analysis of prostate cancer cell lines (data not shown). This MTI method helps to preserve valuable tissue and while multiplexing prognostic biomarkers for prostate cancer.
Statistical analysis
GraphPad Prism 6 software was used to generate scatterplots of the biomarker expressions within the different Gleason score groups (3+3, 3+4, 4+3, and ≥8) as well as to generate the predictive probability plots. One-way ANOVA analysis followed by Dunnet multiple comparisons test was used to evaluate the biomarker expressions in the four Gleason score groups. STATA 13.0 (STATA, StataCorp LP) was used for all statistical modeling. Statistical significance was defined as a P < 0.05. Statistical modeling included standard and/or backwards stepwise multivariate logistic regression (MLR) to discriminate less aggressive from more aggressive prostate cancer. Decision curve analysis (DCA; refs. 27–29) was used to evaluate different MLR models for the prediction of more aggressive prostate cancer, patient that experienced BCR, and progressors. DCA is a method for evaluating and comparing different prediction models (27–29), which gives an expected net benefit per patient relative to the assumption that all patients are treated or not treated. For these evaluations, treatment was defined as the significant/interested changes in the patients (aggressiveness, BCR, progressors, etc.). The interpretation of net benefit is made by comparing a model curve to a baseline curve where all patients are considered as interested changes, and if the model curve is above the standard curve at a certain probability, then the model would be considered a better predictor of the outcome.
Results
Differential expression of CACNA1D, HER2/neu, and Periostin in prostate cancer with different Gleason scores
We optimized all MTI primary antibodies for six of our interested biomarkers on a test TMA (TMA 475) containing normal control and prostate cancer tissue samples. Our results showed that proteins on a single 5 μm FFPE tissue slide could be successfully transferred onto a series of six membranes and probed with six biomarkers simultaneously (Fig. 1).
After optimization of the six biomarkers, we performed MTI with these six biomarkers on our study, TMAs sections containing cancer and caner-adjacent benign tissue from 80 different prostate cancer patients and determined the relative expression of the biomarkers. These sections include four groups with different Gleason scores (3+3, 3+4, 4+3, and ≥ 8). Among the six biomarkers, CACNA1D, HER2/neu, and Periostin were differentially regulated among the four Gleason score groups in cancer area and cancer-adjacent benign areas. In the cancer areas, CACNA1D and HER2/neu relative expression were significantly lower in the Gleason score groups of 4+3 and ≥ 8 compared with the Gleason score group of 3+3 (P < 0.05). There was a trend of gradually decreased relative expression of both markers with increased Gleason scores (Fig. 2A and 2B). Pearson correlation analysis showed that CACNA1D had a strong negative correlation with Gleason score (correlation coefficient = −0.73, P < 0.001), whereas HER2/neu showed a moderate negative correlation (correlation coefficient = −0.39, P = 0.001). Interestingly, CACNA1D and HER2/neu expression behaved similarly in the cancer-adjacent benign areas (Fig. 2A and 2B vs. 2D and 2E), suggesting a possible field effect for these markers. Periostin expression was significantly higher in prostate cancer cases with Gleason score ≥ 8 compared with that of Gleason score 3+3 in cancer-adjacent benign areas (P < 0.05, Fig. 2F). Although Periostin expression was not statistically significantly different in cancer areas, there was a trend for the expression to increase with increasing Gleason score (Fig. 2C).
Biomarker expression in TMAs distinguishes less aggressive from more aggressive prostate cancer
Although the intermediate risk group of clinically localized prostate has a Gleason score of 7, it includes both Gleason 3+4 and Gleason 4+3, patients with the two pathologic statuses have significantly different treatment: most Gleason 3+4, if it is organ-confined, could be managed conservatively while all Gleason 4+3 necessitates immediate intervention with surgery and/or adjuvant therapy. Therefore, we grouped Gleason score 3+3 and 3+4 cases as a less aggressive phenotype and Gleason score 4+3 and ≥8 cases as a more aggressive phenotype. Using MLR, our results show that a model retaining CACNA1D, HER2/neu, and Ki67 expression distinguishes less aggressive and more aggressive prostate cancer in cancer areas (Fig. 3A and Supplementary Tables S2 and S3). Using this model, the predictive probability was close to 100% for most cases in the aggressive groups (Gleason 4+3 and ≥ 8; Fig. 3B). The ORs of Ki-67, CACNA1D, HER2/neu were 4.41e7, <0.01, and 0.02, respectively (Supplementary Table S2). Interestingly, an MLR model that retained CACNA1D and Periostin expression was able to distinguish less aggressive from more aggressive prostate cancer in cancer-adjacent benign areas (Fig. 3D and Supplementary Tables S2 and S3). The predictive probability was also close to 100% for most cases in the aggressive groups (Fig. 3E).The ORs for Periostin and CACNA1D were 1.14 and <0.01, respectively (Supplementary Table S2). The DCA curves showed that prediction of aggressiveness using a molecular biomarker (MB) model in both cancer areas (Fig. 3C) and cancer-adjacent benign areas (Fig. 3F) had superior advantage over that of treating all patients as more aggressive.
Biomarker expression in TMAs predicts BCR in prostate cancer
In the 80 TMA cases, standard logistic regression using Gleason score alone was able to predict BCR with a ROC-AUC of 0.71, while a stepwise MLR model that combines Gleason score with CACNA1D and HER2/neu yields a better ROC-AUC of 0.79 for the prediction of BCR (Fig. 4A and Supplementary Tables S2 and S3). However, there was no significant difference between these two models of BCR prediction (P = 0.15). The combination model could be used to predict recurrence with ORs of 6.34, 3.83, and 6.77 for Gleason score, CACNA1D, and HER2/neu, respectively. DCA curve showed only a limited power for prediction of BCR with Gleason score alone within the probability range of 0.10 to 0.50, while combination of Gleason score with molecular biomarkers showed improved prediction capability within a wider probability range (0.10–0.80; Fig. 4B). The predictive probability of BCR using Gleason score alone or combination of Gleason score with MBs is shown in Fig. 4C and 4D, respectively.
Biomarkers and tPSA separate progressors and nonprogressors in the active surveillance cohort
Using the biomarkers evaluated in the TMAs, we also evaluated the predictive ability of the biomarker molecular profiles in a total of 61 AS biopsy cases (32 nonprogressors and 29 progressors). In the AS program, progressors are defined as men that met all the above VLR criteria at entry (tumor increased in volume and/or Gleason score > 6); however, during monitoring, they were found to have a more aggressive prostate cancer than entry criteria on biopsy follow-up (i.e., AS event; or clinically aggressive prostate cancer). An MLR model that retained (−5, −7) proPSA and CACNA1D was able to distinguish progressors from nonprogressors in the AS cohort (Fig. 5A). This model had an accuracy of 72.13% in the prediction of progressors. The ORs for (−5, −7) proPSA and CACNA1D are 0.20 and 7.93, respectively. The predictive probability of this model is shown in Fig. 5D.
We further tried to separate progressors and nonprogressors in the AS cohort by using prebiopsy tPSA alone or adding the tPSA result to the biomarker model. We found that tPSA alone has limited prediction capability with an ROC-AUC of 0.73, whereas combining molecular biomarkers and tPSA improved the ROC-AUC of the model to 0.83, but without significant difference compared with the biomarker-only model (P = 0.26). The combination model (tPSA-MB) had a sensitivity of 79.31%, specificity of 78.13%, with the tPSA had an OR of 1.24 in this model. The combined model could correctly classify the two groups with an accuracy of 78.69% (Fig. 5A and Supplementary Tables S4 and S5). A plot of the predictive probabilities is shown in Fig. 5C and 5E. DCA curve showed that the tPSA only and the molecular biomarker-only model improved the net benefit within the probability range of 0.35 to 0.80 and 0 to 0.80, respectively, while the combined tPSA-MB model improved the net benefit within the probability range of 0 to 0.75. DCA curve comparing the three models to each other showed that the combined model was superior to the biomarker-only model within the probability range of 0.30 to 0.65 (Fig. 5B). This suggests that a subset of the biomarkers evaluated have the potential to differentiate progressors from nonprogressor patients that can remain annual monitoring in our preliminary data of 61 biopsy cases with high specificity. Adding tPSA improves the model, especially within the probability range of 0.30 to 0.65.
Discussion
The development and validation of sensitive and more accurate methods for early detection are pivotal in the management of prostate cancer, the overdiagnosis and overtreatment of which is a significant healthcare problem (13). AS has become the preferred method for the management of VLR and many LR patients. Currently, the new NCCN guidelines have included VLR and LR group of prostate cancer for AS (23). For VLR or LR cases, abnormal morphologic changes are minimal, and the future tumor progress could not be effectively predicted by pathologic grading. Meanwhile, molecular level changes do occur and may signal the potential disease progression. Characterization of prostate cancer necessitates a comprehensive panel of biomarkers considering the heterogeneity of prostate cancer biology (30). Therefore, we reason that a comprehensive panel of biomarkers involving the glandular acini and the stromal areas could be useful in the prediction of VLR and LR prostate cancer progression in an AS monitoring program.
To make full use of the valuable biopsy samples with limited size of cancer area, we adopted MTI as a method to profile the expression of six biomarkers on a single 5-μm section in TMAs of 80 RP cases. Using a backward stepwise MLR model, we found CACNA1D, HER2/neu, and Periostin expression could significantly distinguish less aggressive from more aggressive prostate cancer with high sensitivity and specificity in the cancer area. Interestingly, in the cancer-adjacent benign area, a different MLR model retains CACNA1D and Periostin could also distinguish less aggressive from more aggressive prostate cancer, implying the importance of these biomarkers in the prediction of aggressiveness in RP cases. Moreover, our data of the cancer-adjacent benign area support the concept of a possible field effect that involves changes of tumor progression at the molecular level before changes in glandular morphology. Our model also showed an improved prediction of BCR compared with the prediction using Gleason score alone. Furthermore, in a total of 61 AS cases, we successfully separated progressors and nonprogressors in the AS cohort with (−5, −7) proPSA and CACNA1D among these verified biomarkers, which are retained in a MLR model for molecular biomarkers to predict the AS cases that require definitive treatment.
It is worthwhile to note that our RP long-term follow-up MLR biomarker model differs from the cancer area and the cancer-adjacent benign area suggests that biology of molecular alterations in progression and changes in the field may differ in the cancer areas. This is supported by a systemic genomic analysis of prostate cancer and morphologically normal tissue (30). Notably, the changes of biomarkers in the cancer-adjacent benign area show a slightly different pattern compared with the cancer area, with changes of stromal protein Periostin expression being more significant. Further characterization of the molecular features of prostate cancer with VLR or LR biomarker expression patterns in these two histologic areas could provide an additional means to stratify the AS and RP cases to more accurately assess the progression of prostate cancer and effectively personalize the management of treatment.
Interestingly, CACNA1D stands out as a new useful marker in prostate cancer area, cancer-adjacent benign area, and also biopsies in the prediction of cancer progression (upgrading). CACNA1D is an ERG target gene upregulated in TMPRSS2-ERG–positive prostate cancer cells (31). Overexpression of CACNA1D in prostate cancer cells was reported to regulate the ligand independent activation of AR. ERG-induced expression of CACNA1D was reported to promote entry of calcium ions into cytosol (32). Here, we found that the expression of CACNA1D showed a trend of decrease with increasing Gleason score in RP specimens. However, the progressor AS biopsies have relatively higher expression of CACNA1D as an early event that consists with results in our RP model for 3+3 and 3+4 Gleason scores. It has been reported that blocking L-type calcium channel-mediated Ca2+ influx suppressed castration-induced apoptotic cell death in prostate epithelial cells (33, 34). We presume that the reduced expression of CACNA1D in more aggressive prostate cancer and progressors may contribute to decreased Ca2+ influx and results in the escape of apoptosis, thus suggesting that decreased CACNA1D expression could be a more aggressive phenotype in higher grade prostate cancer or progressors in AS.
The role of HER2/neu in prostate cancer remains controversial (35, 36). HER2/neu can often be overexpressed in prostate cancer (37, 38), but may depend on the antibody reagents used and two differentially expressed forms of HER2/neu transcripts were reported (39), which may need to be considered during interpretation of results. Antibodies detecting variable epitopes of HER2/neu may be one reason for the controversy. In the current study, we showed that in cancer area, along with CACNA1D and Gleason score, HER2/neu could be used to efficiently predict aggressiveness and recurrence of prostate cancer in RP or biopsy specimens.
Moreover, prostate cancer progression involves both epithelial changes and stromal changes. Periostin was identified as an important biomarker well correlating with the aggressive phenotype (40–42). Both stromal and epithelial Periostin expression was reported to be significantly increased in tumor tissues: stromal expression was significantly higher than epithelial expression as compared with normal tissue. Although high stromal expression was significantly associated with shorter survival, a low epithelial score significantly correlated with shorter PSA-free survival, suggesting that Periostin apparently plays an opposing biologic role depending on its tissue localization (43). Here, we included Periostin in our biomarker panel quantification in both epithelia and stroma since our MTI method could not separate the expression of these two areas. Our results suggest that along with CACNA1D, Periostin could predict the malignancy of prostate cancer based on the cancer-adjacent area (Fig. 3B).
Our MLR modeling system showed its power in the prediction of prostate cancer progression in biopsies, and RP cases with high sensitivity and specificity. Our early data were validated on 80 RP cases of our TMAs and 61 AS biopsies at diagnosis (i.e., entry to AS). In the latter case, the AS study was compromised by the fact that the total cancer surface area in the AS biopsies is often limited. Hence, due to this minimal cancer area in VLR prostate cancer, the expansion of the AS biopsy sample size is needed for obtaining additional results for the verification of the statistical stability of our models. Alternatively, the AS program can be altered to expand criteria for entry into the program, which is already being done at JHMI and other institutions (21, 44). In addition, MTI shows its advantage in multiplexing up to six biomarkers using the same 5-μm tissue section, it does have technical disadvantages: it relies on the relatively high expression of protein and the use of good quality antibodies. The detection of lowly expressed protein levels is limited, such as Ki-67 in AS biopsies. Finally, even if super-thin P-film was used in MTI, the protein transferred onto the membranes shows a gradual decrease from the one at the bottom closest to the tissue to the one on top. We normalized the expression of each biomarker to the total protein to balance these differences in total protein content on the same membrane.
Notably, we have developed a panel of six biomarkers profiled on the same 5-μm biopsy tissue slide using MTI, an integrated, high-throughput quantitative molecular biomarker-based method for predicting prostate cancer progression. Our model involving these biomarkers shows robust strength in the prediction of aggressiveness in the RP cases and more importantly, this method and biomarkers could be used to predict the progress of VLR and LR prostate cancer at diagnosis, which will reduce the overdiagnosis and overtreatment and greatly facilitate the effective management of prostate cancer patients. Future studies will significantly expand the number of biopsy cases studied and also include nuclear morphometry of Feulgen (DNA)-stained or H&E-stained biopsies of AS patients as a variable to more accurately predict the prostate cancer aggressive phenotype of progressors. Given the new problems of managing AS cohorts, new criteria for selection and new technology for multiplexing tissue biomarkers are needed to support the advancement of this concept to improve healthcare in the future for prostate cancer.
Disclosure of Potential Conflicts of Interest
C.S. Christudass has ownership interest in a patent filed in the United States (PCT/US2015/015074). No potential conflicts of interest were disclosed by the other authors.
Authors' Contributions
Conception and design: Z. Liu, C.S. Christudass, H. Zhang, S.M. Hewitt, R.W. Veltri
Development of methodology: G. Zhu, Z. Liu, C.S. Christudass, S.M. Hewitt, R.W. Veltri
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): G. Zhu, Z. Liu, J.I. Epstein, C. Davis, C.S. Christudass, H.B. Carter, P. Landis
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): G. Zhu, Z. Liu, C. Davis, M.C. Miller, R.W. Veltri
Writing, review, and/or revision of the manuscript: G. Zhu, Z. Liu, J.I. Epstein, C.S. Christudass, H. Zhang, J.-Y. Chung, M.C. Miller, R.W. Veltri
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): G. Zhu, C. Davis, H.B. Carter, P. Landis, J.-Y. Chung, S.M. Hewitt, M.C. Miller, R.W. Veltri
Study supervision: G. Zhu, H.B. Carter, H. Zhang, R.W. Veltri
Acknowledgments
The authors appreciate the help of Department of Pathology in Johns Hopkins Medical Institutions and the PCBN for the construction of TMA. The P-Films used in MTI were kindly provided by Vladimir Knezevic in 20/20 GeneSystems, Inc. The authors also appreciate the help of Nan Hu and Chaoyu Wang in Dr. Philip Taylor's lab at NCI for assisting in the conduction of experiments.
Grant Support
This work was supported by two grants: EDRN/NCI U01CA152813 grant (to H. Zhang P.I.) with an administrative supplement to Dr. Robert W. Veltri, Prostate Cancer Foundation and Patana Fund of The Brady Urological Institute at the Johns Hopkins University.
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.