Comparing Platforms for Messenger RNA Expression Profiling of Archival Formalin-Fixed, Paraffin-Embedded Tissues

Intratumoral vitamin D signaling and lethal prostate cancer

High circulating vitamin D levels and supplementation may lower prostate cancer mortality. To probe for direct effects of vitamin D signaling in the primary tumor, we assessed how activation of intratumoral vitamin D signaling in prostate cancer is associated with lethal prostate cancer during long-term follow-up. Among 404 participants with primary prostate cancer in the Health Professionals Follow-up Study and the Physicians' Health Study, we defined a gene score of expected activated intratumoral vitamin D signaling consisting of transcriptionally upregulated (CYP27A1, CYP2R1, RXRA, RXRB, and VDR) and downregulated genes (CYP24A1 and DHCR7). We contrasted vitamin D signaling in tumors that progressed to lethal disease (metastases/prostate cancer-specific death, n = 119) over up to three decades of follow-up with indolent tumors that remained nonmetastatic for >8 years post-diagnosis (n = 285). The gene score was downregulated in tumor tissue compared with tumor-adjacent histologically normal tissue of the same men. Higher vitamin D gene scores were inversely associated with lethal prostate cancer (odds ratio for highest versus lowest quartile: 0.46, 95% confidence interval: 0.21-0.99) in a dose-response fashion and after adjusting for clinical and pathologic factors. This association appeared strongest among men with high predicted plasma 25-hydroxyvitamin D3 and men with body mass index ≥25 kg/m2. Findings were replicated with broader gene sets. These data support the hypothesis that active intratumoral vitamin D signaling is associated with better prostate cancer outcomes and provide further rationale for testing how vitamin D-related interventions after diagnosis could improve prostate cancer survival through effects on the tumor.

LINCATRA: Two-cycle method to amplify RNA for transcriptome analysis from formalin-fixed paraffin-embedded tissue

Whole transcriptome analysis (WTA) using RNA extracted from Formalin Fixed Paraffin Embedded (FFPE) tissue is an invaluable tool to understand the molecular pathology of disease. RNA extracted from FFPE tissue is either degraded and/or in very low quantities hampering gene expression analysis. Earlier studies described protocols applied for cellular RNA using poly-A primer-based linear amplification. The current study describes a method, LINCATRA (LINear amplifiCAtion of RNA for whole TRAnscriptome analysis). It employs random nonamer primer based method which can amplify short, fragmented RNA with high fidelity from as low as 5 ng to obtain enough material for WTA. The two-cycle method significantly amplified RNA at ∼1000 folds (p < 0.0001) improving the mean read lengths (p < 0.05) in WTA. Overall, increased mean read length positively correlated with on-target reads (Pearson's r = 0.71, p < 0.0001) in both amplified and unamplified RNA-seq analysis. Gene expression analysis compared between unamplified and amplified group displayed substantial overlap of the differentially expressed genes (DEGs) (log2 fold change cut-off < -2 and >2, p < 0.05) identified between lung cancer and asthma cohorts validating the method developed. This method is applicable in clinical molecular pathology field for both diagnostics and elucidation of disease mechanisms.

Prognostic significance of EGFR, AREG and EREG amplification and gene expression in muscle invasive bladder cancer

Introduction

Muscle invasive bladder cancer (MIBC) remains a prevalent cancer with limited therapeutic options, obviating the need for innovative therapies. The epidermal growth factor receptor (EGFR) is a linchpin in tumor progression and presents a potential therapeutic target in MIBC. Additionally, the EGFR ligands AREG and EREG have shown associations with response to anti-EGFR therapy and improved progression-free survival in colorectal carcinoma.

Materials and methods

We investigated the prognostic significance of EGFR, AREG, and EREG in MIBC. Gene expression and copy number analyses were performed via qRT-PCR on tissue samples from 100 patients with MIBC who underwent radical cystectomy at the University Hospital Mannheim (MA; median age 72, interquartile range [IQR] 64-78 years, 25% female). Results were validated in 361 patients from the 2017 TCGA MIBC cohort (median age 69, IQR 60-77 years, 27% female), in the Chungbuk and MDACC cohort. Gene expressions were correlated with clinicopathologic parameters using the Mann-Whitney test, Kruskal-Wallis- test and Spearman correlation. For overall survival (OS), cancer-specific survival (CSS) and disease-free survival (DFS) gene expression was analyzed with Kaplan-Meier and Cox-proportional hazard models.

Results

Significant gene expression differences in EGFR, AREG, and EREG could be detected in all cohorts. In the TCGA cohort, EGFR expression was significantly elevated in patients with EGFR amplification and KRAS wildtype. High AREG expression independently predicted longer OS (HR = 0.35, CI 0.19 - 0.63, p = 0.0004) and CSS (HR = 0.42, CI 0.18 - 0.95, p = 0.0378) in the MA cohort. In the TCGA cohort, high EGFR, AREG, and EREG expression correlated with shorter OS (AREG: HR = 1.57, CI 1.12 - 2.20, p = 0.0090) and DFS (EGFR: HR = 1.91, CI 1.31 - 2.8, p = 0.0008). EGFR amplification was also associated with reduced DFS.

Discussion

High EGFR and EREG indicate worse survival in patients with MIBC. The prognostic role of AREG should further be investigated in large, prospective series. Divergent survival outcomes between the four cohorts should be interpreted cautiously, considering differences in analysis methods and demographics. Further in vitro investigations are necessary to elucidate the functional mechanisms underlying the associations observed in this study.

Transcriptomes of Prostate Cancer with TMPRSS2:ERG and Other ETS Fusions

The most common somatic event in primary prostate cancer is a fusion between the androgen-related TMPRSS2 gene and the ERG oncogene. Tumors with these fusions, which occur early in carcinogenesis, have a distinctive etiology. A smaller subset of other tumors harbor fusions between TMPRSS2 and members of the ETS transcription factor family other than ERG. To assess the genomic similarity of tumors with non-ERG ETS fusions and those with fusions involving ERG, this study derived a transcriptomic signature of non-ERG ETS fusions and assessed this signature and ERG-related gene expression in 1,050 men with primary prostate cancer from three independent population-based and hospital-based studies. Although non-ERG ETS fusions involving ETV1, ETV4, ETV5, or FLI1 were individually rare, they jointly accounted for one in seven prostate tumors. Genes differentially regulated between non-ERG ETS tumors and tumors without ETS fusions showed similar differential expression when ERG tumors and tumors without ETS fusions were compared (differences explained: R2 = 69-77%), including ETS-related androgen receptor (AR) target genes. Differences appeared to result from similarities among ETS tumors rather than similarities among non-ETS tumors. Gene sets associated with ERG fusions were consistent with gene sets associated with non-ERG ETS fusions, including fatty acid and amino acid metabolism, an observation that was robust across cohorts.

IMPLICATIONS

Considering ETS fusions jointly may be useful for etiologic studies on prostate cancer, given that the transcriptome is profoundly impacted by ERG and non-ERG ETS fusions in a largely similar fashion, most notably genes regulating metabolic pathways.

Identification of Novel Diagnostic Markers for Malignant Pleural Mesothelioma Using a Reverse Translational Approach Based on a Rare Synchronous Tumor

Although the routine use of immunohistochemistry has improved the accuracy of histopathologic diagnosis in clinical practice, new methods for discovering novel diagnostic markers are still needed. We sought new diagnostic markers for malignant pleural mesothelioma (MPM) using a reverse translational approach with limited archival tissues from a very rare case. Total RNA extracted from formalin-fixed paraffin-embedded (FFPE) tissues of a synchronous collision tumor consisting of MPM and pulmonary adenocarcinoma (PAC) was employed for gene expression profiling (GEP) analysis. Among the 54 genes selected by GEP analysis, we finally identified the following two candidate MPM marker genes: PHGDH and TRIM29. Immunohistochemical analysis of 48 MM and 20 PAC cases showed that both PHGDH and TRIM29 had sensitivity and specificity almost equivalent to those of calretinin (sensitivity 50% and 46% vs. 63%, and specificity 95% and 100% vs. 100%, respectively). Importantly, of the 23 epithelioid MMs, all 3 calretinin-negative cases were positive for TRIM29. These two markers may be diagnostically useful for immunohistochemical distinction between MPMs and PACs. This successful reverse translational approach based on FFPE samples from one very rare case encourages the further use of such samples for the development of novel diagnostic markers.

High IL-22RA1 gene expression is associated with poor outcome in muscle invasive bladder cancer

BACKGROUND

The cell surface interleukin 22 (IL-22) receptor complex is mainly expressed in epithelial and tissue cells like pancreatitis cells. Recent studies described that IL-22R was overexpressed in malignant diseases and was associated with a poor overall survival (OS). The role of IL-22RA1 gene expression in muscle invasive bladder cancer (MIBC) has not been investigated, yet.

OBJECTIVES

The aim of this study was to analyze the role of IL-22RA1 gene expression in patients with MIBC.

METHODS

In a cohort of 114 patients with MIBC who underwent radical cystectomy, IL-22RA1 gene expression was analyzed with qRT-PCR and correlated with clinical parameters. Furthermore, Kaplan-Meier and Cox regression analysis were performed. For validation, an in silico dataset (TCGA 2017, n=407) was reanalyzed.

RESULTS

IL-22RA1 gene expression was independent of clinicopathological parameters like age (P=0.2681), T stage (P=0.2130), nodal status (P=0.3238) and lymph vascular invasion (LVI, P=0.5860) in patients with MIBC. A high expression of IL-22RA1 was associated with a shorter OS (P=0.0040) and disease-specific survival (P=0.0385). Furthermore, a shorter disease-free survival (DFS) was also associated with a high expression of IL-22RA1 (P=0.0102). In the multivariable analysis, IL-22RA1 expression was an independent prognostic predictors regarding OS (P=0.0096, HR=0.48). In the TCGA cohort, IL-22RA1 expression was independent regarding to OS and DFS.

CONCLUSION

A high IL-22RA1 gene expression was associated with worse outcome. Furthermore, IL-22RA1 represented an independent predictor regarding OS in our cohort and therefore might be used for risk stratification in patients with MIBC.

Gene Expression Pathways in Prostate Tissue Associated with Vigorous Physical Activity in Prostate Cancer

BACKGROUND

Men engaged in high physical activity have lower risks of advanced and fatal prostate cancer. Mechanisms underlying this association are not well understood but may include systemic and tumor-specific effects. We investigated potential mechanisms linking physical activity and gene expression in prostate tissue from men with prostate cancer.

METHODS

We included a subset of 118 men in the Health Professionals Follow-up Study diagnosed with prostate cancer between 1986 and 2005 with whole-transcriptome gene expression profiling on tumor and adjacent normal prostate tissue and physical activity data. Long-term vigorous physical activity was self-reported as the average time spent engaged in various forms of recreational physical activity at baseline and biennially until prostate cancer diagnosis. Gene set enrichment analysis was performed among KEGG and Hallmark gene sets to identify pathways with differential expression based on vigorous physical activity.

RESULTS

In adjacent normal tissue, we identified 25 KEGG gene sets enriched (downregulated) in the highest compared with lowest quintile of vigorous physical activity at an FDR <0.10, including a number of cancer- and immune-related pathways. Although no gene sets reached statistical significance in tumor tissue, top gene sets differentially expressed included TGF beta, apoptosis, and p53 signaling pathways.

CONCLUSIONS

These findings suggest that physical activity may influence the tumor microenvironment. Future studies are needed to confirm these findings and further investigate potential mechanisms linking physical activity to lethal prostate cancer.

IMPACT

Identification of gene expression alterations in the prostate associated with physical activity can improve our understanding of prostate cancer etiology.

Suppression of ARID1A associated with decreased CD8 T cells improves cell survival of ovarian clear cell carcinoma

Objective

AT-rich interactive domain 1A (ARID1A) plays an important role as a tumor suppressor gene in ovarian clear cell carcinoma (OCCC), but the clinical application of ARID1A remains unclear. The aim of this study was to analyze clinicopathological parameters, molecular interactions and immune-infiltration in patients with low ARID1A expression and to provide candidate target drugs.

Methods

We investigated the clinicopathologic parameters, specific gene sets/genes, and immunological relevance according to ARID1A expression in 998 OCCC patients from 12 eligible studies (using meta-analyses); 30 OCCC patients from the Hanyang University Guri Hospital (HYGH) cohort; and 52 OCCC patients from gene set enrichment (GSE) 65986 (25 patients), 63885 (9 patients), and 54809 (6 patients and 12 healthy people) of the Gene Expression Omnibus (GEO). We analyzed network-based pathways based on gene set enrichment analysis (GSEA) and performed in vitro drug screening.

Results

Low ARID1A expression was associated with poor survival in OCCC from the meta-analysis, HYGH cohort and GEO data. In GSEA, low ARID1A expression was related to the tumor invasion process as well as a low immune-infiltration. In silico cytometry showed that CD8 T cells were decreased with low ARID1A expression. In pathway analysis, ARID1A was associated with angiogenic endothelial cell signaling. In vitro drug screening revealed that cabozantinib and bicalutamide effectively inhibited specific hub genes, such as vascular endothelial growth factor-A and androgen receptor, in OCCC cells with low ARID1A expression.

Conclusions

Therapeutic strategies making use of low ARID1A could contribute to better clinical management/research for patients with OCCC.

Extensive Changes in Transcriptomic "Fingerprints" and Immunological Cells in the Large Organs of Patients Dying of Acute Septic Shock and Multiple Organ Failure Caused by Neisseria meningitidis

Background: Patients developing meningococcal septic shock reveal levels of Neisseria meningitidis (10⁶-10⁸/mL) and endotoxin (10¹-10³ EU/mL) in the circulation and organs, leading to acute cardiovascular, pulmonary and renal failure, coagulopathy and a high case fatality rate within 24 h. Objective: To investigate transcriptional profiles in heart, lungs, kidneys, liver, and spleen and immunostain key inflammatory cells and proteins in post mortem formalin-fixed, paraffin-embedded (FFPE) tissue samples from meningococcal septic shock patients. Patients and Methods: Total RNA was isolated from FFPE and fresh frozen (FF) tissue samples from five patients and two controls (acute non-infectious death). Differential expression of genes was detected using Affymetrix microarray analysis. Lung and heart tissue samples were immunostained for T-and B cells, macrophages, neutrophils and the inflammatory markers PAI-1 and MCP-1. Inflammatory mediators were quantified in lysates from FF tissues. Results: The transcriptional profiles showed a complex pattern of protein-coding and non-coding RNAs with significant regulation of pathways associated with organismal death, cell death and survival, leukocyte migration, cellular movement, proliferation of cells, cell-to-cell signaling, immune cell trafficking, and inflammatory responses in an organ-specific clustering manner. The canonical pathways including acute phase response-, EIF2-, TREM1-, IL-6-, HMBG1-, PPAR signaling, and LXR/RXR activation were associated with acute heart, pulmonary, and renal failure. Fewer genes were regulated in the liver and particularly in the spleen. The main upstream regulators were TNF, IL-1β, IL-6, RICTOR, miR-6739-3p, and CD3. Increased numbers of inflammatory cells (CD68+, MPO+, CD3+, and CD20+) were found in lungs and heart. PAI-1 inhibiting fibrinolysis and MCP-1 attracting leukocyte were found significantly present in the septic tissue samples compared to the controls. Conclusions: FFPE tissue samples can be suitable for gene expression studies as well as immunostaining of specific cells or molecules. The most pronounced gene expression patterns were found in the organs with highest levels of Neisseria meningitidis DNA. Thousands of protein-coding and non-coding RNA transcripts were altered in lungs, heart and kidneys. We identified specific biomarker panels both protein-coding and non-coding RNA transcripts, which differed from organ to organ. Involvement of many genes and pathways add up and the combined effect induce organ failure.

Expression Profiling Reveals Involvement of WNT Pathway in the Malignant Progression of Sessile Serrated Adenomas

Approximately 15% to 20% of colorectal cancers are developed through the serrated pathway of tumorigenesis, which is associated with BRAF mutation, CpG island methylation phenotype, and MLH1 methylation. However, the detailed process of progression from sessile serrated adenoma (SSA) to dysplasia and carcinoma has not been elucidated. To further characterize mechanisms involved in the dysplastic progression of SSA, we investigated differential expressions of mRNAs between areas with and without dysplasia within the same SSA polyps. Significantly dysregulated genes in paired samples were applied for functional annotation and biological significance. The same lysates from a subset of matched samples were subjected for miRNA expression profiling. Differentially expressed miRNAs were determined, and their targeted mRNAs were compared in parallel to the list of differentially expressed mRNAs from an RNA sequencing study. Fourteen common mRNA targets were identified, which include AXIN2, a known indicator of WNT/β-catenin pathway activation. Together, in this study, different genes, pathways, and biological processes involved in the initiation and progression of dysplasia in the serrated pathway are documented. One of the most significant findings is the involvement of the WNT/β-catenin pathway in the dysplastic progression of SSAs with different genes being targeted in early versus advanced dysplasia.

Expression Concordance of 325 Novel RNA Biomarkers between Data Generated by NanoString nCounter and Affymetrix GeneChip

Background

With the development of new drug combinations and targeted treatments for multiple types of cancer, the ability to stratify categories of patient populations and to develop companion diagnostics has become increasingly important. A panel of 325 RNA biomarkers was selected based on cancer-related biological processes of healthy cells and gene expression changes over time during nonmalignant epithelial cell organization. This "cancer in reverse" approach resulted in a panel of biomarkers relevant for at least 7 cancer types, providing gene expression profiles representing key cellular signaling pathways beyond mutations in "driver genes." Objective. To further investigate this biomarker panel, the objective of the current study is to (1) validate the assay reproducibility for the 325 RNA biomarkers and (2) compare gene expression profiles side by side using two technology platforms.

Methods and Results

We have mapped the 325 RNA transcripts and in a custom NanoString nCounter expression panel to be compared to all potential probe sets in the Affymetrix Human Genome U133 Plus 2.0. The experiments were conducted with 10 unique biological formalin-fixed paraffin-embedded (FFPE) breast tumor samples. Each site extracted RNA from four sections of 10-micron thick FFPE tissue over three different days by two different operators using an optimized standard operating procedure and quality control criteria. Samples were analyzed using mas5 in BioConductor and NanoStringNorm in R. Pearson correlation showed reproducibility between sites for all 60 samples with r = 0.995 for Affymetrix and r = 0.999 for NanoString. Correlation in multiple days and multiple users was for Affymetrix r = (0.962 - 0.999) and for NanoString r = (0.982 - 0.991).

Conclusion

The 325 RNA biomarkers showed reproducibility in two technology platforms with moderate to high concordance. Future directions include performing clinical validation studies and generating rationale for patient selection in clinical trials using the technically validated assay.

Single-nucleotide polymorphisms in DNMT3B gene and DNMT3B mRNA expression in association with prostate cancer mortality

BACKGROUND

Germline variants in DNA methyltransferase 3B (DNMT3B) may influence DNMT3B enzymatic activity, which, in turn, may affect cancer aggressiveness by altering DNA methylation.

METHODS

The study involves two Italian cohorts (NTAT cohort, n = 157, and 1980s biopsy cohort, n = 182) and two U.S. cohorts (Health Professionals Follow-Up Study, n = 214, and Physicians' Health Study, n = 298) of prostate cancer (PCa) patients, and a case-control study of lethal (n = 113) vs indolent (n = 290) PCa with DNMT3B mRNA expression data nested in the U.S. cohorts. We evaluated the association between: three selected DNMT3B variants and global DNA methylation using linear regression in the NTAT cohort, the three DNMT3B variants and PCa mortality using Cox proportional hazards regression in all cohorts, and DNMT3B expression and lethal PCa using logistic regression, with replication in publicly available databases (TCGA, n = 492 and MSKCC, n = 140).

RESULTS

The TT genotype of rs1569686 was associated with LINE-1 hypomethylation in tumor tissue (β = -2.71, 95% CI: -5.41, -0.05). There was no evidence of association between DNMT3B variants and PCa mortality. DNMT3B expression was consistently associated with lethal PCa in the two U.S. cohorts (3rd vs 1st tertile, combined cohorts: OR = 2.04, 95% CI: 1.13, 3.76); the association was replicated in TCGA and MSKCC data (3rd vs 1st tertile, TCGA: HR = 3.00, 95% CI: 1.78, 5.06; MSKCC: HR = 2.22, 95% CI: 1.01, 4.86).

CONCLUSIONS

Although there was no consistent evidence of an association between DNMT3B variants and PCa mortality, the TT genotype of rs1569686 was associated with LINE-1 hypomethylation in tumor tissue and DNMT3B mRNA expression was associated with an increased risk of lethal PCa.

ARHGEF38 as a novel biomarker to predict aggressive prostate cancer

Prostate cancer (PCa) metastasis is considered the leading cause of cancer death in males. Therapeutic strategies and diagnosis for stage-specific PCa have not been well understood. Rho guanine nucleotide exchange factor 38 (ARHGEF38) is related to tumor cell polarization and is frequently expressed in PCa. Microarray data of PCa were downloaded from GEO and TCGA databases. A total of 243 DEGs were screened, of which, 32 genes were upregulated. The results of enrichment analysis showed the participation of these DEGs in the tumor cell metastasis pathway. ARHGEF38 was significantly up-regulated in the four most prevalent cancers worldwide (p < 0.05), and its expression was higher in the tumor samples with higher Gleason score (GS). IHC, qRT-PCR, and western-blot analyses showed the higher expression of ARHGEF38 in PCa than benign prostatic hyperplasia (BPH). In addition, IHC results demonstrated a higher expression of ARHGEF38 in high-grade PCa than the low-grade PCa.

Differential Gene Expression in Prostate Tissue According to Ejaculation Frequency

In a prospective study of 31 925 men with 18 yr of follow-up, higher ejaculation frequency (EF) throughout adulthood was associated with lower rates of prostate cancer. To further explore this association, we evaluated whole transcriptome gene expression in the prostate tissue from study participants who developed prostate cancer between 1992 and 2004 (n=157 tumor tissue, n=85 adjacent normal). We tested for trends in gene expression according to the level of EF as self-reported in 1992 for ages 20-29 yr, 40-49 yr, and the year prior to the questionnaire, 1991. There were no associations between EF and gene expression in areas of tumor after accounting for multiple testing. In contrast, in the adjacent normal tissue, 409 genes and six pathways were differentially expressed at a false discovery rate ≤0.2 across categories of EF in 1991. These results suggest that ejaculation affects the expression of genes in the normal prostate tissue. The identified genes and pathways provide potential biological links between EF and prostate tumorigenesis.

PATIENT SUMMARY

To explore previous findings that men who ejaculate more frequently have lower risk of prostate cancer, we evaluated molecular alterations in the prostate tissue according to each man's frequency of ejaculation prior to diagnosis. We identified biological processes that could link ejaculation frequency and prostate cancer.

Gene expression profiling of prostate tissue identifies chromatin regulation as a potential link between obesity and lethal prostate cancer

BACKGROUND

Obese men are at higher risk of advanced prostate cancer and cancer-specific mortality; however, the biology underlying this association remains unclear. This study examined gene expression profiles of prostate tissue to identify biological processes differentially expressed by obesity status and lethal prostate cancer.

METHODS

Gene expression profiling was performed on tumor (n = 402) and adjacent normal (n = 200) prostate tissue from participants in 2 prospective cohorts who had been diagnosed with prostate cancer from 1982 to 2005. Body mass index (BMI) was calculated from the questionnaire immediately preceding cancer diagnosis. Men were followed for metastases or prostate cancer-specific death (lethal disease) through 2011. Gene Ontology biological processes differentially expressed by BMI were identified using gene set enrichment analysis. Pathway scores were computed by averaging the signal intensities of member genes. Odds ratios (ORs) for lethal prostate cancer were estimated with logistic regression.

RESULTS

Among 402 men, 48% were healthy weight, 31% were overweight, and 21% were very overweight/obese. Fifteen gene sets were enriched in tumor tissue, but not normal tissue, of very overweight/obese men versus healthy-weight men; 5 of these were related to chromatin modification and remodeling (false-discovery rate < 0.25). Patients with high tumor expression of chromatin-related genes had worse clinical characteristics (Gleason grade > 7, 41% vs 17%; P = 2 × 10^-4 ) and an increased risk of lethal disease that was independent of grade and stage (OR, 5.26; 95% confidence interval, 2.37-12.25).

CONCLUSIONS

This study improves our understanding of the biology of aggressive prostate cancer and identifies a potential mechanistic link between obesity and prostate cancer death that warrants further study. Cancer 2017;123:4130-4138. © 2017 American Cancer Society.

The role of tumor metabolism as a driver of prostate cancer progression and lethal disease: results from a nested case-control study

Background

Understanding the biologic mechanisms underlying the development of lethal prostate cancer is critical for improved therapeutic and prevention strategies. In this study we explored the role of tumor metabolism in prostate cancer progression using mRNA expression profiling of seven metabolic pathways; fatty acid metabolism, glycolysis/gluconeogenesis, oxidative phosphorylation, pentose phosphate, purine metabolism, pyrimidine metabolism and the tricarboxylic acid cycle.

Methods

The study included 404 men with archival formalin-fixed, paraffin-embedded prostate tumor tissue from the prospective Health Professionals Follow-up Study and Physicians’ Health Study. Lethal cases (n = 113) were men who experienced a distant metastatic event or died of prostate cancer during follow-up. Non-lethal controls (n = 291) survived at least 8 years post-diagnosis without metastases. Of 404 men, 202 additionally had matched normal tissue (140 non-lethal, 62 lethal). Analyses compared expression levels between tumor and normal tissue, by Gleason grade and by lethal status. Secondary analyses considered the association with biomarkers of cell proliferation, apoptosis and angiogenesis.

Results

Oxidative phosphorylation and pyrimidine metabolism were identified as the most dysregulated pathways in lethal tumors (p < 0.007), and within these pathways, a number of novel differentially expressed genes were identified including POLR2K and APT6V1A. The associations were tumor specific as there was no evidence any pathways were altered in the normal tissue of lethal compared to non-lethal cases.

Conclusions

The results suggest prostate cancer progression and lethal disease are associated with alterations in key metabolic signaling pathways. Pathways supporting proliferation appeared to be of particular importance in prostate tumor aggressiveness.

Electronic supplementary material

The online version of this article (doi:10.1186/s40170-016-0161-9) contains supplementary material, which is available to authorized users.

Early Molecular Stratification of High-risk Primary Biliary Cholangitis

High-risk primary biliary cholangitis (PBC), defined by inadequate response at one year to Ursodeoxycholic acid (UDCA), is associated with disease progression and liver transplantation. Stratifying high-risk patients early would facilitate improved approaches to care. Using long-term follow-up data to define risk at presentation, 6 high-risk PBC patients and 8 low-risk patients were identified from biopsy, transplant and biochemical archival records. Formalin-fixed paraffin-embedded (FFPE) liver biopsies taken at presentation were graded (Scheuer and Nakanuma scoring) and gene expression analysed using the NanoString® nCounter PanCancer Immunity 770-gene panel. Principle component analysis (PCA) demonstrated discrete gene expression clustering between controls and high- and low-risk PBC. High-risk PBC was characterised by up-regulation of genes linked to T-cell activation and apoptosis, INF-γ signalling and leukocyte migration and down-regulation of those linked to the complement pathway. CDKN1a, up-regulated in high-risk PBC, correlated with significantly increased expression of its gene product, the senescence marker p21^WAF1/Cip, by biliary epithelial cells. Our findings suggest high- and low-risk PBC are biologically different from disease outset and senescence an early feature in high-risk disease. Identification of a high-risk 'signal' early from standard FFPE tissue sections has clear clinical utility allowing for patient stratification and second-line therapeutic intervention.

Prognostic Utility of a New mRNA Expression Signature of Gleason Score

PURPOSE

Gleason score strongly predicts prostate cancer mortality; however, scoring varies among pathologists, and many men are diagnosed with intermediate-risk Gleason score 7. We previously developed a 157-gene signature for Gleason score using a limited gene panel. Using a new whole-transcriptome expression dataset, we verified the previous signature's performance and developed a new Gleason signature to improve lethal outcome prediction among men with Gleason score 7.

EXPERIMENTAL DESIGN

We generated mRNA expression data from prostate tumor tissue from men in the Physicians' Health Study and Health Professionals Follow-Up Study (N = 404) using the Affymetrix Human Gene 1.0 ST microarray. The Prediction Analysis for Microarrays method was used to develop a signature to distinguish high (≥8) versus low (≤6) Gleason score. We evaluated the signature's ability to improve prediction of lethality among men with Gleason score 7, adjusting for 3 + 4/4 + 3 status, by quantifying the area under the receiver operating characteristic (ROC) curve (AUC).

RESULTS

We identified a 30-gene signature that best distinguished Gleason score ≤6 from ≥8. The AUC to predict lethal disease among Gleason score 7 men was 0.76 [95% confidence interval (CI), 0.67-0.84] compared with 0.68 (95% CI, 0.59-0.76) using 3 + 4/4 + 3 status alone (P = 0.0001). This signature was a nonsignificant (P = 0.09) improvement over our previous signature (AUC = 0.72).

CONCLUSIONS

Our new 30-gene signature improved prediction of lethality among men with Gleason score 7. This signature can potentially become a useful prognostic tool for physicians to improve treatment decision making. Clin Cancer Res; 23(1); 81-87. ©2016 AACRSee related commentary by Yin et al., p. 6.

Stress-Related Signaling Pathways in Lethal and Nonlethal Prostate Cancer

PURPOSE

Recent data suggest that neuroendocrine signaling may influence progression in some cancers. We aimed to determine whether genes within the five major stress-related signaling pathways are differentially expressed in tumor tissue when comparing prostate cancer patients with lethal and nonlethal disease.

EXPERIMENTAL DESIGN

We measured mRNA expression of 51 selected genes involved in predetermined stress-related signaling pathways (adrenergic, glucocorticoid, dopaminergic, serotoninergic, and muscarinic systems) in tumor tissue and normal prostate tissue collected from prostate cancer patients in the Physicians' Health Study (n = 150; n = 82 with normal) and the Health Professionals Follow-Up Study (n = 254; n = 120 with normal). We assessed differences in pathway expression in relation to prostate cancer lethality as the primary outcome and to biomarkers as secondary outcomes.

RESULTS

Differential mRNA expression of genes within the adrenergic (P = 0.001), glucocorticoid (P < 0.0001), serotoninergic (P = 0.0019), and muscarinic (P = 0.0045) pathways in tumor tissue was associated with the risk of lethality. The adrenergic pathway was also statistically significant (P = 0.001) when comparing against differential expression of genes not involved in the pathways. In adjacent normal prostate tissue, none of the pathways was clearly differentially expressed between lethal and nonlethal prostate cancer. The glucocorticoid and adrenergic pathways were associated with cell proliferation, while the glucocorticoid pathway was additionally associated with angiogenesis and perineural invasion.

CONCLUSIONS

Our study suggests that stress-related signaling pathways, particularly the adrenergic and glucocorticoid, may be dysregulated in the tumors of men whose prostate cancer proves to be lethal, and motivates further investigation of these pathways in functional studies.