The identification of trans-associations between prostate cancer GWAS SNPs and RNA expression differences in tumor-adjacent stroma.

Here we tested the hypothesis that SNPs associated with prostate cancer risk, might differentially affect RNA expression in prostate cancer stroma. The most significant 35 SNP loci were selected from Genome Wide Association (GWA) studies of ~40,000 patients. We also selected 4030 transcripts previously associated with prostate cancer diagnosis and prognosis. eQTL analysis was carried out by a modified BAYES method to analyze the associations between the risk variants and expressed transcripts jointly in a single model. We observed 47 significant associations between eight risk variants and the expression patterns of 46 genes. This is the first study to identify associations between multiple SNPs and multiple in trans gene expression differences in cancer stroma. Potentially, a combination of SNPs and associated expression differences in prostate stroma may increase the power of risk assessment for individuals, and for cancer progression.

Generation of "virtual" control groups for single arm prostate cancer adjuvant trials.

It is difficult to construct a control group for trials of adjuvant therapy (Rx) of prostate cancer after radical prostatectomy (RP) due to ethical issues and patient acceptance. We utilized 8 curve-fitting models to estimate the time to 60%, 65%, … 95% chance of progression free survival (PFS) based on the data derived from Kattan post-RP nomogram. The 8 models were systematically applied to a training set of 153 post-RP cases without adjuvant Rx to develop 8 subsets of cases (reference case sets) whose observed PFS times were most accurately predicted by each model. To prepare a virtual control group for a single-arm adjuvant Rx trial, we first select the optimal model for the trial cases based on the minimum weighted Euclidean distance between the trial case set and the reference case set in terms of clinical features, and then compare the virtual PFS times calculated by the optimum model with the observed PFSs of the trial cases by the logrank test. The method was validated using an independent dataset of 155 post-RP patients without adjuvant Rx. We then applied the method to patients on a Phase II trial of adjuvant chemo-hormonal Rx post RP, which indicated that the adjuvant Rx is highly effective in prolonging PFS after RP in patients at high risk for prostate cancer recurrence. The method can accurately generate control groups for single-arm, post-RP adjuvant Rx trials for prostate cancer, facilitating development of new therapeutic strategies.

Expression differences between African American and Caucasian prostate cancer tissue reveals that stroma is the site of aggressive changes.

In prostate cancer, race/ethnicity is the highest risk factor after adjusting for age. African Americans have more aggressive tumors at every clinical stage of the disease, resulting in poorer prognosis and increased mortality. A major barrier to identifying crucial gene activity differences is heterogeneity, including tissue composition variation intrinsic to the histology of prostate cancer. We hypothesized that differences in gene expression in specific tissue types would reveal mechanisms involved in the racial disparities of prostate cancer. We examined 17 pairs of arrays for AAs and Caucasians that were formed by closely matching the samples based on the known tissue type composition of the tumors. Using pair-wise t-test we found significantly altered gene expression between AAs and CAs. Independently, we performed multiple linear regression analyses to associate gene expression with race considering variation in percent tumor and stroma tissue. The majority of differentially expressed genes were associated with tumor-adjacent stroma rather than tumor tissue. Extracellular matrix, integrin family and signaling mediators of the epithelial-to-mesenchymal transition (EMT) pathways were all downregulated in stroma of AAs. Using MetaCore (GeneGo) analysis, we observed that 35% of significant (p < 10(-3)) pathways identified EMT and 25% identified immune response pathways especially for interleukins-2, -4, -5, -6, -7, -10, -13, -15 and -22 as the major changes. Our studies reveal that altered immune and EMT processes in tumor-adjacent stroma may be responsible for the aggressive nature of prostate cancer in AAs.

An accurate prostate cancer prognosticator using a seven-gene signature plus Gleason score and taking cell type heterogeneity into account.

One of the major challenges in the development of prostate cancer prognostic biomarkers is the cellular heterogeneity in tissue samples. We developed an objective Cluster-Correlation (CC) analysis to identify gene expression changes in various cell types that are associated with progression. In the Cluster step, samples were clustered (unsupervised) based on the expression values of each gene through a mixture model combined with a multiple linear regression model in which cell-type percent data were used for decomposition. In the Correlation step, a Chi-square test was used to select potential prognostic genes. With CC analysis, we identified 324 significantly expressed genes (68 tumor and 256 stroma cell expressed genes) which were strongly associated with the observed biochemical relapse status. Significance Analysis of Microarray (SAM) was then utilized to develop a seven-gene classifier. The Classifier has been validated using two independent Data Sets. The overall prediction accuracy and sensitivity is 71% and 76%, respectively. The inclusion of the Gleason sum to the seven-gene classifier raised the prediction accuracy and sensitivity to 83% and 76% respectively based on independent testing. These results indicated that our prognostic model that includes cell type adjustments and using Gleason score and the seven-gene signature has some utility for predicting outcomes for prostate cancer for individual patients at the time of prognosis. The strategy could have applications for improving marker performance in other cancers and other diseases.

Expression changes in the stroma of prostate cancer predict subsequent relapse.

Biomarkers are needed to address overtreatment that occurs for the majority of prostate cancer patients that would not die of the disease but receive radical treatment. A possible barrier to biomarker discovery may be the polyclonal/multifocal nature of prostate tumors as well as cell-type heterogeneity between patient samples. Tumor-adjacent stroma (tumor microenvironment) is less affected by genetic alteration and might therefore yield more consistent biomarkers in response to tumor aggressiveness. To this end we compared Affymetrix gene expression profiles in stroma near tumor and identified a set of 115 probe sets for which the expression levels were significantly correlated with time-to-relapse. We also compared patients that chemically relapsed shortly after prostatectomy (<1 year), and patients that did not relapse in the first four years after prostatectomy. We identified 131 differentially expressed microarray probe sets between these two categories. 19 probe sets (15 genes overlapped between the two gene lists with p<0.0001). We developed a PAM-based classifier by training on samples containing stroma near tumor: 9 rapid relapse patient samples and 9 indolent patient samples. We then tested the classifier on 47 different samples, containing 90% or more stroma. The classifier predicted the risk status of patients with an average accuracy of 87%. This is the first general tumor microenvironment-based prognostic classifier. These results indicate that the prostate cancer microenvironment exhibits reproducible changes useful for predicting outcomes for patients.

Diagnosis of prostate cancer using differentially expressed genes in stroma.

More than one million prostate biopsies are performed in the United States every year. A failure to find cancer is not definitive in a significant percentage of patients due to the presence of equivocal structures or continuing clinical suspicion. We have identified gene expression changes in stroma that can detect tumor nearby. We compared gene expression profiles of 13 biopsies containing stroma near tumor and 15 biopsies from volunteers without prostate cancer. About 3,800 significant expression changes were found and thereafter filtered using independent expression profiles to eliminate possible age-related genes and genes expressed at detectable levels in tumor cells. A stroma-specific classifier for nearby tumor was constructed on the basis of 114 candidate genes and tested on 364 independent samples including 243 tumor-bearing samples and 121 nontumor samples (normal biopsies, normal autopsies, remote stroma, as well as stroma within a few millimeters of tumor). The classifier predicted the tumor status of patients using tumor-free samples with an average accuracy of 97% (sensitivity = 98% and specificity = 88%) whereas classifiers trained with sets of 100 randomly generated genes had no diagnostic value. These results indicate that the prostate cancer microenvironment exhibits reproducible changes useful for categorizing the presence of tumor in patients when a prostate sample is derived from near the tumor but does not contain any recognizable tumor.

In silico estimates of tissue components in surgical samples based on expression profiling data.

Tissue samples from many diseases have been used for gene expression profiling studies, but these samples often vary widely in the cell types they contain. Such variation could confound efforts to correlate expression with clinical parameters. In principle, the proportion of each major tissue component can be estimated from the profiling data and used to triage samples before studying correlations with disease parameters. Four large gene expression microarray data sets from prostate cancer, whose tissue components were estimated by pathologists, were used to test the performance of multivariate linear regression models for in silico prediction of major tissue components. Ten-fold cross-validation within each data set yielded average differences between the pathologists' predictions and the in silico predictions of 8% to 14% for the tumor component and 13% to 17% for the stroma component. Across independent data sets that used similar platforms and fresh frozen samples, the average differences were 11% to 12% for tumor and 12% to 17% for stroma. When the models were applied to 219 arrays of "tumor-enriched" samples in the literature, almost one quarter were predicted to have 30% or less tumor cells. Furthermore, there was a 10.5% difference in the average predicted tumor content between 37 recurrent and 42 nonrecurrent cancer patients. As a result, genes that correlated with tissue percentage generally also correlated with recurrence. If such a correlation is not desired, then some samples might be removed to rebalance the data set or tissue percentages might be incorporated into the prediction algorithm. A web service, "CellPred," has been designed for the in silico prediction of sample tissue components based on expression data.

Bcl-B expression in human epithelial and nonepithelial malignancies.

PURPOSE: Apoptosis plays an important role in neoplastic processes. Bcl-B is an antiapoptotic Bcl-2 family member, which is known to change its phenotype upon binding to Nur77/TR3. The expression pattern of this protein in human malignancies has not been reported. EXPERIMENTAL DESIGN: We investigated Bcl-B expression in normal human tissues and several types of human epithelial and nonepithelial malignancy by immunohistochemistry, correlating results with tumor stage, histologic grade, and patient survival. RESULTS: Bcl-B protein was strongly expressed in all normal plasma cells but found in only 18% of multiple myelomas (n = 133). Bcl-B immunostaining was also present in normal germinal center centroblasts and centrocytes and in approximately half of diffuse large B-cell lymphoma (n = 48) specimens, whereas follicular lymphomas (n = 57) did not contain Bcl-B. In breast (n = 119), prostate (n = 66), gastric (n = 180), and colorectal (n = 106) adenocarcinomas, as well as in non-small cell lung cancers (n = 82), tumor-specific overexpression of Bcl-B was observed. Bcl-B expression was associated with variables of poor prognosis, such as high tumor grade in breast cancer (P = 0.009), microsatellite stability (P = 0.0002), and left-sided anatomic location (P = 0.02) of colorectal cancers, as well as with greater incidence of death from prostate cancer (P = 0.005) and shorter survival of patients with small cell lung cancer (P = 0.009). Conversely, although overexpressed in many gastric cancers, Bcl-B tended to correlate with better outcome (P = 0.01) and more differentiated tumor histology (P < 0.0001). CONCLUSIONS: Tumor-specific alterations in Bcl-B expression may define subsets of nonepithelial and epithelial neoplasms with distinct clinical behaviors.

In silico dissection of cell-type-associated patterns of gene expression in prostate cancer.

Prostate tumors are complex entities composed of malignant cells mixed and interacting with nonmalignant cells. However, molecular analyses by standard gene expression profiling are limited because spatial information and nontumor cell types are lost in sample preparation. We scored 88 prostate specimens for relative content of tumor, benign hyperplastic epithelium, stroma, and dilated cystic glands. The proportions of these cell types were then linked in silico to gene expression levels determined by microarray analysis, revealing unique cell-specific profiles. Gene expression differences for malignant and nonmalignant epithelial cells (tumor versus benign hyperplastic epithelium) could be identified without being confounded by contributions from stroma that dominate many samples or sacrificing possible paracrine influences. Cell-specific expression of selected genes was validated by immunohistochemistry and quantitative PCR. The results provide patterns of gene expression for these three lineages with relevance to pathogenetic, diagnostic, and therapeutic considerations.