Metabolomic profile of prostate cancer-specific survival among 1812 Finnish men
BMC Medicine volume 20, Article number: 362 (2022)
Abnormal metabolism and perturbations in metabolic pathways play significant roles in the development and progression of prostate cancer; however, comprehensive metabolomic analyses of human data are lacking and needed to elucidate the interrelationships.
We examined the serum metabolome in relation to prostate cancer survival in a cohort of 1812 cases in the Alpha-Tocopherol, Beta-Carotene Cancer Prevention (ATBC) Study. Using an ultrahigh-performance LC-MS/MS platform, we identified 961 known metabolites in prospectively collected serum. Median survival time from diagnosis to prostate cancer-specific death (N=472) was 6.6 years (interquartile range=2.9–11.1 years). Cox proportional hazards regression models estimated hazard ratios and 95% confidence intervals of the associations between the serum metabolites (in quartiles) and prostate cancer death, adjusted for age at baseline and diagnosis, disease stage, and Gleason sum. In order to calculate risk scores, we first randomly divided the metabolomic data into a discovery set (70%) and validated in a replication set (30%).
Overall, 49 metabolites were associated with prostate cancer survival after Bonferroni correction. Notably, higher levels of the phospholipid choline, amino acid glutamate, long-chain polyunsaturated fatty acid (n6) arachidonate (20:4n6), and glutamyl amino acids gamma-glutamylglutamate, gamma-glutamylglycine, and gamma-glutamylleucine were associated with increased risk of prostate cancer-specific mortality (fourth versus first quartile HRs=2.07–2.14; P-values <5.2×10−5). By contrast, the ascorbate/aldarate metabolite oxalate, xenobiotics S-carboxymethyl-L-cysteine, fibrinogen cleavage peptides ADpSGEGDFXAEGGGVR and fibrinopeptide B (1-12) were related to reduced disease-specific mortality (fourth versus first quartile HRs=0.82–0.84; P-value <5.2×10−5). Further adjustment for years from blood collection to cancer diagnosis, body mass index, smoking intensity and duration, and serum total and high-density lipoprotein cholesterol did not alter the results. Participants with a higher metabolic score based on the discovery set had an elevated risk of prostate cancer-specific mortality in the replication set (fourth versus first quartile, HR=3.9, P-value for trend<0.0001).
The metabolic traits identified in this study, including for choline, glutamate, arachidonate, gamma-glutamyl amino acids, fibrinopeptides, and endocannabinoid and redox pathways and their composite risk score, corroborate our previous analysis of fatal prostate cancer and provide novel insights and potential leads regarding the molecular basis of prostate cancer progression and mortality.
With a global estimate of 1.4 million new cases and 375,000 deaths, prostate cancer is the second most commonly diagnosed malignancy and the fifth most frequent cause of cancer death among men in 2020, contributing to a considerable health burden among men both in the USA and worldwide . Despite this reality, prostate cancer etiology and effective means for its prevention remain incompletely understood. The well-recognized risk factors for prostate cancer include older age, African ancestry, family history of prostate cancer, and low-penetrant genetic variants, none of which are modifiable . It is well-accepted that many prostate cancer cases, particularly those diagnosed through prostate-specific antigen (PSA) screening, have microscopic foci of limited (i.e., non-life-threatening) clinical relevance, and research suggests that more aggressive and fatal prostate malignancies have etiologies distinct from indolent disease [2,3,4]. Thus, there is a need for further study of prostate cancer survival.
A poor understanding of the etiology of prostate cancer survival limits our ability to make meaningful recommendations regarding prevention and effective risk stratification for men diagnosed with prostate cancer. Alterations in metabolism and perturbations in metabolic pathways can play significant roles in the development and progression of prostate cancer; however, substantial human data are needed to elucidate these associations [5, 6]. Recent advances of high-throughput metabolomics approaches enable the systematic measurement of a broad spectrum of small compounds, and the unique metabolic profiles can reflect an integrated final molecular phenotype of endogenous and exogenous exposure, including influence from gene-environment interactions, that will facilitate the identification of novel risk factors and metabolic pathways associated with disease progression [7,8,9].
Metabolomic analyses of nested case-control sets within three major cohorts have explored associations between pre-diagnostic levels of circulating metabolites and risk of prostate cancer [10,11,12,13,14]. Associations between urine and tissue metabolites and risk of prostate cancer have also been reported in some clinical studies [15, 16]. To our knowledge, only two untargeted metabolomic analyses have been conducted examining fatal prostate cancer specifically. One investigation included 523 prostate cancer cases in the Alpha-Tocopherol, Beta-Carotene Cancer Prevention (ATBC) Study  and showed risk associations with glutathione-reactive oxygen species (ROS) metabolites including thioproline, cysteine, and cystine, while another smaller study found endocannabinoid N-oleoyl taurine and sterol/steroids metabolites related to prostate cancer survival . Thus, there remains a substantial research gap for the investigation of circulating metabolite profiles in relation to prostate cancer survival, that urgently warrants investigation.
We conducted a prospective untargeted metabolomic analysis with the aim to identify pre-diagnostic serum metabolic profiles and metabolic pathways that are associated with survival among 1812 prostate cancer patients within the ATBC Study. We hypothesize that unique metabolomic profiles would be associated with prostate cancer-specific survival.
The ATBC Study was a double-blind, 2×2 factorial, randomized, placebo-controlled trial that enrolled participants from 1985 to 1988, and the primary aim of this intervention study was to test whether supplementation of alpha-tocopherol and beta-carotene could reduce cancer incidence, especially for lung cancer. The ATBC Study enrolled 29,133 Finnish male smokers, who were aged 50 to 69 years old from southwest Finland, and the eligible participants were assigned into 4 groups to receive the supplement daily (alpha-tocopherol, beta-carotene, both, or placebo) for 5–8 years until the end of the intervention (April 30, 1993). At enrollment, participants completed questionnaires that included detailed information on lifestyle risk factors and a separate validated food frequency questionnaire. Height, weight, heart rate, and blood pressure were measured, and overnight fasting blood samples were collected by research nurses. The blood samples were then aliquoted and stored at −70 °C until assayed.
Via linkage to the Finnish Cancer Registry by using the International Classification of Disease (ICD-9) code of 185, we identified 1812 incident prostate cancer cases diagnosed through December 31, 2012, with available non-thawed baseline serum and non-missing cancer stage data. We used the underlying cause of death from the Statistics Finland Death Registry to identify prostate cancer-specific mortality, based on the code 185 of ICD-9 and C61 of ICD-10. In the present study, we found 472 prostate cancer deaths that occurred during the observational follow-up period through December 31, 2016. The survival time was computed from the date of prostate cancer diagnosis to the date of prostate cancer-specific mortality, or to the censor date (date of death due to other causes or December 31, 2016), whichever occurred first.
Metabolomic profiling assays
Serum metabolite profiling was conducted on a high-resolution accurate mass (HRAM) platform, ultrahigh-performance liquid chromatography/tandem mass spectrometry (LC-MS/MS) at Metabolon Inc. We measured 1195 metabolites, of which 234 were excluded due to being unidentified metabolites or metabolites for which more than 90% of subjects had values below the limit of detection, leaving 961 metabolites in the final analysis. For those metabolites remaining after these exclusions were applied, undetectable values were assigned to one-half minimum value of the metabolite. Based on biochemical information from the Kyoto Encyclopedia of Genes and Genomics (KEGG) database, and the Human Metabolome Database, identified known metabolites were classified into different chemical classes: amino acids or amino acid derivatives (further referred to as “amino acids”), carbohydrates, cofactors and vitamins, energy biochemicals, lipids, nucleotides, peptides or xenobiotics. In order to assess the reliability and reproducibility of the biochemical data, each metabolite’s intraclass correlation coefficient (ICC) and coefficient of variation (CV) were computed based on blinded quality control samples embedded in each batch (6%). The median ICC and CV for measured metabolites were 0.95 (interquartile range=0.87 to 0.98) and 0.12 (interquartile range=0.07 to 0.23), respectively.
Serum levels of alpha-tocopherol and retinol quantified by metabolomic profiling were highly correlated with concentrations originally assayed for the entire cohort at baseline using the reversed-phase isocratic high-performance liquid chromatography (HPLC) method, suggesting excellent validity and reproducibility for the present metabolomic assay platform (alpha-tocopherol: r=0.75, P-value <10−8; retinol, r=0.77, P-value <10−8).
In order to control for batch variation, metabolites were batch-normalized and divided by the batch median value. Then each metabolite was log-transformed and normalized with a mean of 0 and a variance of 1. We used Cox proportional hazards regression models to calculate the hazard ratios (HRs) and 95% confidence intervals (CIs) of the associations between the serum metabolites and prostate cancer death, adjusted for age at baseline, age at prostate cancer diagnosis, disease stage, and Gleason sum, and we modeled standardized metabolites as both categorical (quartiles) and continuous variables (per 1-standard deviation [SD] log-metabolite increment). In a separate Cox proportional regression model, we tested for linear trend by assigning the ordinal value of the quartile and treating this as a continuous variable. To control for multiple hypothesis testing , metabolite-risk associations with a Bonferroni corrected threshold (a P-value of 5.2×10−5 or less, 0.05/961) were considered statistically significant.
In additional analyses, we further adjusted for body mass index (BMI, weight (kg)/height (m)2), number of cigarettes smoked daily, years of smoking, total serum cholesterol (mmol/L), serum high-density lipoprotein cholesterol (mmol/L), and calendar year of cancer diagnosis (all as continuous variables). We also performed analyses stratified by time from blood collection to prostate cancer diagnosis (median split), disease stage at diagnosis (stage I–II and stage III–IV) and BMI (<25, 25–<30, and ≥30 kg/m2). In the entire dataset, we conducted a forward stepwise Cox proportional hazards regression to determine the set of metabolites that were conditionally independently associated with prostate cancer death, with an alpha-error, the Bonferroni corrected threshold of .05/961 to enter the model (i.e., top 49 signals) and P=.05 to remain in the model. For comparison of the ability for different clinical factors and selected metabolites to distinguish between prostate cancer-specific mortality and non-prostate cancer-specific mortality, we utilized adjusted receiver operating characteristics (ROC) and adjusted area under the curve (AUC) analyses. Differences between AUC were tested using the method proposed by DeLong et al. ,
Pathway analyses evaluated the associations between chemical sub-classes of metabolites and risk of prostate cancer-specific mortality. Due to the correlations between metabolites in the sub-class pathways, we used a parametric bootstrap method to calculate the P-value from each pathway level (based on 100,000 permutations), with the assumption that the Z-statistics of the component metabolites followed a multivariate normal distribution with a mean of 0 and estimated covariance matrix [20, 21]. For pathways associated with prostate cancer-specific mortality, we conducted principal component analysis (PCA) using the varimax rotation method and defined the “pathway score” by the first generated principal component. Cox proportional hazards regression models were applied to estimate the HRs of prostate cancer-specific mortality for the resultant PCA score.
In order to calculate risk scores, we first randomly divided the metabolomic data into a discovery set (70%) and a replication set (30%). Within the discovery set, we used Cox proportional hazards regression models to select metabolites (Bonferroni corrected threshold, a P-value of 5.2×10−5 or less, 0.05/961) associated with prostate cancer-specific mortality (number of metabolites=5). We then created a risk score (sum) of these metabolites in the replication set weighted by their coefficients as obtained from the Cox regression models in the discovery set. The metabolite risk score was then normalized and modeled as both quartiles and as a continuous variable (with a mean=0 and SD=1) to estimate the association with prostate cancer-specific mortality, according to the Cox regression model and adjusted for the above-mentioned covariates.
All analyses were conducted using SAS version 9.4 (SAS Institute, Cary, NC), and R version 3.6.1 (R Development Core Team, Vienna, Austria). The reported statistical tests were two-sided.
Baseline and clinical characteristics of the 1812 prostate cancer cases being studied are presented in Table 1. Through the end of follow-up, we observed 1522 deaths (84.0%) of which 472 died from prostate cancer (26.1% of all cases). Comparing cases who died from prostate cancer to those who did not, the median age at diagnosis was 69.0 years (interquartile range=65.0 to 73.0) and 71.0 years (interquartile range= 68.0 to 75.0), respectively, and the median time from baseline blood collection to cancer diagnosis was 11.0 years (interquartile range=8.0 to 15.0) and 14.0 (interquartile range=11.0 to 18.0), respectively. The cancer stage distribution among cases who died from prostate cancer was 39.4% (cancer stage I or II), 15.3% (locally advanced disease, stage III), and 45.3% (metastatic disease, stage IV), whereas for those who did not die from prostate cancer it was 81.9%, 12.6%, and 5.5%, respectively.
Cox proportional regression models found that 49 of the 961 identified metabolites were associated with prostate cancer mortality at the Bonferroni corrected P-value threshold of 5.2×10−5 or less, including 18 peptides, 14 amino acids, 12 lipids, two cofactors/vitamins, and one each of nucleotides, carbohydrates, and xenobiotics (Table 2). The two strongest metabolite associations with prostate cancer death were for the lipid choline (fourth versus first quartile, HR=2.07, 95% CI:1.60–2.69, Ptrend across quartiles=7.2×10−11), and the amino acid glutamate (HR=2.31, 95% CI: 1.77–3.02, Ptrend=6.7×10−10). Other statistically significant positive associations included the following: 11 gamma-glutamyl amino acids (fourth versus first quartile, HRs=1.58–2.14; 1.3×10−9≤P≤3.1×10−5), three endocannabinoids (HRs=1.62–1.69; 3.1×10−5≤P≤4.6×10−5), three dipeptides (HRs=1.83–1.97; 3.0×10−7≤P≤4.9×10−6), two metabolites in the glycine-serine-threonine metabolism pathway (HRs=1.81–1.95; 3.7×10−8≤P≤3.3×10−6), two metabolites in the methionine-cysteine-SAM-taurine pathway (HRs=1.96–2.06; 2.8×10−8≤P≤2.5×10−7), two lysophospholipids (HRs=1.61–1.76; 2.7×10−6≤P≤2.7×10−5), one fibrinogen cleavage peptide, fibrinopeptide B (1-9) (HR=1.80, 95% CI:1.38–2.36; Ptrend=4.0×10−6), and the glutathione metabolite 5-oxoproline (HR=1.67, 95% CI:1.30–2.16; Ptrend=2.3×10−6). We also observed significant inverse associations for three fibrinogen cleavage peptides (fourth versus first quartile, HRs=0.54–0.60; 5.1×10−6≤P≤5.2×10−5), the three glutathione metabolites cysteine-glutathione disulfide, Cys-Gly, and oxidized cysteinylglycine (HRs=0.53–0.59; 1.4 ×10−5≤P≤4.1×10−5), and two ascorbate-aldarate metabolites, oxalate and threonate (HRs=0.75 and 0.80; Ptrend=3.0×10−6 and 3.2×10−5, respectively) (Table 2).
Further model adjustment for BMI, number of cigarettes smoked daily, years of smoking, serum total cholesterol, serum high-density lipoprotein cholesterol and calendar year of prostate cancer diagnosis did not alter the observed associations (Additional File 1: Tables S1).
Figure 1 shows the heat-map of Pearson correlation coefficients for metabolites that were significantly associated with prostate cancer-specific mortality after Bonferroni correction. Stronger positive correlations were found within the gamma-glutamyl amino acid chemical sub-class pathway.
Stratifying cases by cancer stage at diagnosis, the observed metabolite associations with prostate cancer mortality were generally stronger for men with earlier stage cancers (i.e., I/II), with statistically significant stronger associations for 5-oxoproline (1-SD HR 1.54 and 1.17 for earlier versus late cancer stage), fibrinopeptide B (1–12) (HR 0.70 and 0.90), fibrinopeptide B (1–11) (HR 0.70 and 0.91), gamma-glutamylisoleucine (HR 1.57 and 1.17), and gamma-glutamylleucine (HR 1.55 and 1.18) (all P for interaction <3.4×10−4, Bonferroni correction for 49×3 tests, Additional File 1: Tables S2). There was no statistically significant effect modification of metabolite associations by BMI or time from blood collection to cancer diagnosis (Additional File 1: Tables S1 ).
Using the forward stepwise selection analysis, we identified five metabolites that were independently statistically significantly associated with prostate cancer death up to model step five, including gamma-glutamylglutamate, heptenedioate (C7:1-DC), oxalate, arachidonate (20:4n6), and linoleoyl ethanolamide (Additional File 1: Tables S3). The corresponding adjusted ROC had an adjusted AUC of 0.82 (95% CI = 0.80 to 0.85) for clinical factors including age at blood collection, age at diagnosis, cancer stage at diagnosis, Gleason score at diagnosis, cigarettes smoked per day, and BMI, and 0.86 (95% CI = 0.84 to 0.88) for those clinical factors and five metabolites selected from the forward stepwise regression models, including gamma-glutamylglutamate, heptenedioate (C7:1-DC), oxalate, arachidonate (20:4n6) and linoleoyl ethanolamide. The adjusted AUC was significantly influenced and improved from 0.82 to 0.86 by the five selected metabolites (AUC 0.86 versus 0.82, P value = 0.04; Fig. 2).
Pathway and principal component analyses highlighted the potentially relevant importance of eight metabolic chemical pathways that were significantly associated with prostate cancer mortality at the Bonferroni-corrected threshold of 5.3×10−4, including endocannabinoid, glycine-serine-and-threonine, gamma-glutamyl amino acid, glutamate, eicosanoid, sphingolipid synthesis, glutathione, and dipeptide (Table 3). Seven of these pathways were positively associated with prostate cancer-specific mortality (fourth versus first quartile: HRs=1.61–1.94, 4.5×10−7≤P-value≤2.0×10−4), and only the glutathione metabolism pathway was inversely associated (fourth versus first quartile: HR=0.62, P-value=1.7×10−4).
We combined the five metabolites selected from Cox regression model in the 70% discovery set with those achieving Bonferroni corrected threshold (<5.2×10−5) to generate a metabolite risk score and observed statistically significantly increased mortality in the replication set for cases having a higher score (per 1-SD HR=1.48, 95% CI: 1.30 to 1.68, Table 4). In the categorized analysis of the risk score quartiles, we found cases in the second through fourth quartiles experienced an 83%, 85%, and 290% higher risk of prostate cancer-specific mortality, when compared with those in the lowest quartile (P for trend<0.0001, Table 4).
In this large prospective metabolomic analysis, we found 49 metabolites that were statistically significantly associated with prostate cancer mortality, with top signals including choline, glutamate, arachidonate, and gamma-glutamyl amino acids. Additionally, we identified metabolites in the endocannabinoid, glycine-serine-and-threonine-metabolism, eicosanoid, sphingolipid synthesis, dipeptide, and redox pathways being most strongly associated with risk of prostate cancer mortality. Of note, the risk in men with prostate cancer was elevated fourfold in those in the highest metabolite risk score category.
As the top signal in this present metabolomic analysis, serum choline was positively associated with risk of prostate cancer mortality. Abnormal choline metabolism has been characterized as a cholinic phenotype and a metabolic hallmark, related to increased oncogenic signaling and tumor progression . In a case-only survival analysis within the Health Professionals Follow-Up Study, the authors showed that post-diagnostic intake of dietary choline was related to an increased risk of lethal prostate cancer, with a 70% increased risk for men in the highest quintile of intake . Furthermore, a prospective nested case-control study within the ATBC cohort showed that men with greater serum choline levels had a 3-fold increased risk of developing colorectal cancer .
The second top signal identified was glutamate, with higher levels related to an elevated risk of prostate cancer mortality. Glutamate is catabolized from glutamine and in turn incorporated into the tricarboxylic acid (TCA) cycle and lipogenesis . A large portion of glutamine-derived glutamate is secreted from cells via cystine/glutamate transporters in exchange for cystine, which plays vital roles in redox homeostasis . Moreover, the secreted glutamate can stimulate cancer cell proliferation via glutamate receptors , consistent with serum glutamate being correlated with Gleason sum in prostate cancer patients  and in vivo studies demonstrating that reduced glutamate excretion inhibited tumor xenograft growth, including prostate cancer [27, 28]. In addition, expression of the cystine/glutamate transporter SLC7A11 is upregulated in different cancers including prostate [26, 29].
Consistent with our previous analysis of fatal prostate cancer risk in a smaller set of 523 cases-control sets , serum dipeptides including leucylglycine, valylleucine, glycylvaline, and gamma-glutamyl amino acids were associated with an increased risk of prostate cancer-specific mortality here. As a key enzyme in the gamma-glutamyl cycle, gamma-glutamyl peptidase (GGT) is known for degradation of extracellular glutathione and liberation of free gamma-glutamyl peptides and is characterized as a clinical indicator of prostate cancer . Experimental data also demonstrate that GGT can be involved in cellular detoxification and apoptotic balance pathways . For example, in vitro and in vivo studies report that inhibition of the GGT pathway can impact cancer cell growth through cell-cycle regulation, and enhance the sensitivity to chemotherapy . Consistent with our findings, earlier human data showed that circulating gamma-glutamyl peptides were associated with an increased risk of hepatocellular carcinoma, and elevated serum GGT was related to an increased risk of overall and site-specific cancers, including prostate cancer .
Perturbations in endocannabinoid and glycine-serine-threonine metabolism have been related to the risk of prostate cancer survival. Endocannabinoids are endogenous ligands of cannabinoid receptors (CB1 and CB2) , and dysregulation of the endocannabinoid system has been associated with multiple health outcomes including cardiovascular disease, neurodegenerative disorders, and cancer . Several components of the endocannabinoid system were overexpressed in prostate cancer cell lines and correlated with prostate cancer grade and progression in epithelial cells of prostate cancer tissue samples . Our current and previous findings show positive associations between serum endocannabinoids (i.e., N-oleoyl taurine, N-stearoyl serine, N-oleoyl serine, and linoleoyl ethanolamide) and risk of prostate cancer survival . Furthermore, biosynthesis and metabolism of serine, aspartate, and glycine are required for cancer cell growth and progression [35, 36], for example, malignant cells upregulate de novo serine and glycine synthesis and secretion [37, 38] which may help maintain a supportive tumor microenvironment .
The positive association we observed between eicosanoids and prostate cancer mortality is also biologically plausible in that eicosanoids have been found to promote cancer cell growth in an autocrine fashion by activating their receptors implicated in the regulation of cell proliferation, apoptosis, migration, and invasion . Also, eicosanoid signaling inhibitors (e.g., CAY10404 and celecoxib) decrease prostate cancer cell proliferation , and eicosanoids are reported to have immunosuppressive effects which alter the tumor microenvironment and T cell profiles, potentially serving as a therapeutic target for prostate cancer .
Alterations in sphingolipid biosynthesis and metabolism have been associated with elevated prostate cancer mortality, and serum sphingolipids concentrations are related to colorectal cancer and hepatocellular carcinoma [41, 42]. Sphingolipids with their fatty acid tail and sphingoid ring base have significant cell membrane roles in regulating various cellular biological processes including cell proliferation, apoptosis, and inflammatory response [43, 44], and emerging evidence supports their potential use as cancer therapeutic targets . The essential amino acid methionine also plays a critical role in tumor cell proliferation and metabolism , including glutathione formation, polyamine synthesis, and methyl group donation . Reduced methionine status through dietary restriction or enzymatic degradation can inhibit prostate tumor growth . Another of our metabolite signals, arginine, also plays multiple roles in cellular physiology, including nitric oxide and polyamines synthesis, cell proliferation, cellular signaling, and immune system regulation [47, 48], and can activate rapamycin complex 1 (mTORC1), a nutrient-sensing kinase directly involved in tumorigenesis .
We also observed associations between elevated fatty acid metabolites arachidonate, 9-hydroxyoctadecadienoic acid (HODE), and 13-HODE and increased prostate cancer mortality. 9-HODE and 13-HODE can be synthesized by cyclooxygenase (COX) and 15-lipoxygenase 1(15-LOX 1) from linoleic acid [49, 50], and activation of the 15-lipoxygenase 1/13(S)-HODE axis can promote prostate cancer cell growth through the epidermal growth factor receptor (EGFR) signaling pathway . Consistent with our findings, a nested case-control study in the Prostate, Lung, Colorectal, and Ovarian (PLCO) Cancer Screening Trial cohort reported positive associations between serum concentrations of 9-HODE and 13-HODE and increased risk of developing ovarian cancer, highlighting the possible roles of linoleic acid metabolites in cancer etiology . In addition, N6-methyladenosine (m6A) is one of the most prevalent internal posttranscriptional modifications of RNA and participates in RNA metabolism, including mRNA translation, splicing, folding, and degradation . Growing evidence suggests a key role for m6A in tumorigenesis, including for prostate cancer  and its progression , and tumor data from the Cancer Genome Atlas database indicate overexpression of m6A methylation regulators in the aggressive prostate tumor tissues . These findings are in line with the current data suggesting a positive association between serum m6A and prostate cancer mortality.
We have validated findings from our previous metabolomics studies including gamma-glutamyl amino acids, fibrinopeptides, endocannabinoids, and redox metabolites. Additionally, based on the larger sample size of the present study, we identified new metabolites associated with poor prostate cancer survival, including choline, glutamate, and arachidonate. When we compared the present findings with the 34 metabolites identified in our previous study of lethal prostate cancer in 523 case-control sets , 25 metabolites were validated in the present study, three metabolites were not measured in the current panel (histidylalanine, valyglycine, and leucylglutamine), four metabolites were marginally associated in the same direction but with lower statistical significance (P value range of 0.08 to 0.11; C-glycosyltryptophan, 1-linoleoyl-GPC (18:2), 5,6-dihydrouridine and 5-methyluridine), and only two metabolites were not validated by the present study (cystine and pseudouridine, HRs=1). Differences in our findings between the current analysis and our previous studies may be due to its larger sample size and different study design; e.g., our previous study compared lethal prostate cancer cases to cancer-free controls, whereas the present study examined prostate cancer survival among men diagnosed with the disease. Further studies will be needed to understand these nuances.
Our study has several important and novel strengths, including its quite large sample size, a prospective study design, and complete, long-term follow-up identification of prostate cancer diagnoses and mortality outcomes. We used a stringent statistical threshold, namely the Bonferroni correction, to reduce the potential for false-positive findings. Importantly, metabolomic profiles were measured in serum samples that were collected up to a median of 14 years (interquartile range = 10 to 17 years) prior to prostate cancer diagnosis, thus minimizing the possibility of reverse causality. In addition, the untargeted metabolic platform with excellent laboratory reproducibility was highly reliable and enabled us to measure more than 900 known metabolites simultaneously reflecting a wide range of biological pathways relevant to prostate carcinogenesis. Our study is one of the few prospective cohort studies that collected blood samples from participants after an overnight fast, minimizing post-prandial fluctuations in the metabolome which affect most other studies. It should also be noted that the ATBC Study was conducted in Finland, where population-based PSA screening was not implemented, resulting in most of the cases being clinically diagnosed which allowed us to examine which metabolites were related to fatal disease compared to other clinically relevant cases. In most US cohorts, by contrast, a great deal of investigational “noise” is introduced through the inclusion of indolent, PSA-detected cases. As most of the studies have been focused on the metabolomics and prostate cancer incidence, our findings highlight multiple potentially important novel metabolites and biological pathways that may play key roles in the progression and mortality of prostate cancer. Some study limitations should also be mentioned. We did not have an external validation set and relied on the internal replication set; to our knowledge, however, other studies have not measured this large a number of metabolites using an HRAM LC-MS/MS platform. Prostate cancer molecular sub-type data were not available for examination of unique metabolite relationships (i.e., ETS-related gene [ERG]-positive or negative tumors). We cannot preclude residual confounding from other factors that may have influenced the observed associations, although we have controlled for several potential confounding factors in our models, and found no strong evidence of confounding. The homogenous nature of the study population of male Finnish smokers aged 50–69 years may limit the generalizability of our findings to non-smokers, younger men, and other racial/ethnic groups.
Our prospective serum metabolomic investigation of 1812 men with prostate cancer elucidated a panel of metabolites and a composite metabolite risk score from blood collected prior to diagnosis that were associated with later risk of prostate cancer-specific mortality. These findings provide novel, biologically based insights and potential leads regarding the molecular basis of prostate cancer progression and mortality, and may identify a high-risk group of men in whom prevention strategies may be implemented and for whom more specific screening recommendations may be developed.
Availability of data and materials
Because of previously enacted EU General Data Protection Regulation privacy rules and an existing data use agreement between Finland and the US National Cancer Institute, the ATBC Study data and materials described in the article may not be made publicly available for the purpose of reproducing the findings. The principal investigators of the ATBC Study can be contacted with specific data requests (https://atbcstudy.cancer.gov/). To minimize the possibility of unintentionally sharing information that can be used to identify private information, a subset of the data generated for this study will be made available first on reasonable request.
Alpha-Tocopherol, Beta-Carotene Cancer Prevention Study
Area under the curve
Body mass index
Coefficient of variation
Epidermal growth factor receptor
High-performance liquid chromatography
High-resolution accurate mass
Intraclass correlation coefficient
International Classification of Disease
Kyoto Encyclopedia of Genes and Genomics
Liquid chromatography/tandem mass spectrometry
Principal component analysis
Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial cohort
Receiver operating characteristics
Reactive oxygen species
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin. 2021;71(3):209–49.
Hurwitz LM, Agalliu I, Albanes D, Barry KH, Berndt SI, Cai Q, et al. Recommended Definitions of Aggressive Prostate Cancer for Etiologic Epidemiologic Research. J Natl Cancer Inst. 2021;113(6):727–34.
Prensner JR, Rubin MA, Wei JT, Chinnaiyan AM. Beyond PSA: the next generation of prostate cancer biomarkers. Sci Transl Med. 2012;4(127):127rv123.
Platz EA, De Marzo AM, Giovannucci E. Prostate cancer association studies: pitfalls and solutions to cancer misclassification in the PSA era. J Cell Biochem. 2004;91(3):553–71.
Salciccia S, Capriotti AL, Lagana A, Fais S, Logozzi M, De Berardinis E, et al. Biomarkers in Prostate Cancer Diagnosis: From Current Knowledge to the Role of Metabolomics and Exosomes. Int J Mol Sci. 2021;22(9):4367.
Kelly RS, Vander Heiden MG, Giovannucci E, Mucci LA. Metabolomic Biomarkers of Prostate Cancer: Prediction, Diagnosis, Progression, Prognosis, and Recurrence. Cancer Epidemiol Biomarkers Prev. 2016;25(6):887–906.
Nicholson G, Rantalainen M, Maher AD, Li JV, Malmodin D, Ahmadi KR, et al. Human metabolic profiles are stably controlled by genetic and environmental variation. Mol Syst Biol. 2011;7:525.
Playdon MC, Joshi AD, Tabung FK, Cheng S, Henglin M, Kim A, et al. Metabolomics Analytics Workflow for Epidemiological Research: Perspectives from the Consortium of Metabolomics Studies (COMETS). Metabolites. 2019;9(7):145.
Dunn WB, Broadhurst DI, Atherton HJ, Goodacre R, Griffin JL. Systems level studies of mammalian metabolomes: the roles of mass spectrometry and nuclear magnetic resonance spectroscopy. Chem Soc Rev. 2011;40(1):387–426.
Mondul AM, Moore SC, Weinstein SJ, Karoly ED, Sampson JN, Albanes D. Metabolomic analysis of prostate cancer risk in a prospective cohort: The alpha-tocolpherol, beta-carotene cancer prevention (ATBC) study. Int J Cancer. 2015;137(9):2124–32.
Rohnisch HE, Kyro C, Olsen A, Thysell E, Hallmans G, Moazzami AA. Identification of metabolites associated with prostate cancer risk: a nested case-control study with long follow-up in the Northern Sweden Health and Disease Study. BMC Med. 2020;18(1):187.
Huang J, Mondul AM, Weinstein SJ, Koutros S, Derkach A, Karoly E, et al. Serum metabolomic profiling of prostate cancer risk in the prostate, lung, colorectal, and ovarian cancer screening trial. Br J Cancer. 2016;115(9):1087–95.
Huang J, Mondul AM, Weinstein SJ, Derkach A, Moore SC, Sampson JN, et al. Prospective serum metabolomic profiling of lethal prostate cancer. Int J Cancer. 2019;145(12):3231–43.
Schmidt JA, Fensom GK, Rinaldi S, Scalbert A, Appleby PN, Achaintre D, et al. Pre-diagnostic metabolite concentrations and prostate cancer risk in 1077 cases and 1077 matched controls in the European Prospective Investigation into Cancer and Nutrition. BMC Med. 2017;15(1):122.
Cerrato A, Bedia C, Capriotti AL, Cavaliere C, Gentile V, Maggi M, et al. Untargeted metabolomics of prostate cancer zwitterionic and positively charged compounds in urine. Anal Chim Acta. 2021;1158:338381.
Lima AR, Pinto J, Amaro F, Bastos ML, Carvalho M, Guedes de Pinho P. Advances and Perspectives in Prostate Cancer Biomarker Discovery in the Last 5 Years through Tissue and Urine Metabolomics. Metabolites. 2021;11(3):181.
Huang J, Weinstein SJ, Moore SC, Derkach A, Hua X, Mondul AM, et al. Pre-diagnostic Serum Metabolomic Profiling of Prostate Cancer Survival. J Gerontol A Biol Sci Med Sci. 2019;74(6):853–9.
Bland JM, Altman DG. Multiple significance tests: the Bonferroni method. BMJ. 1995;310(6973):170.
DeLong ER, DeLong DM, Clarke-Pearson DL. Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics. 1988;44(3):837–45.
Davison AC, Hinkley DV. Bootstrap Methods and Their Application. Cambridge University Press; 1997.
Chen H, Lumley T, Brody J, Heard-Costa NL, Fox CS, Cupples LA, et al. Sequence kernel association test for survival traits. Genet Epidemiol. 2014;38(3):191–7.
Glunde K, Bhujwalla ZM, Ronen SM. Choline metabolism in malignant transformation. Nat Rev Cancer. 2011;11(12):835–48.
Richman EL, Kenfield SA, Stampfer MJ, Giovannucci EL, Zeisel SH, Willett WC, et al. Choline intake and risk of lethal prostate cancer: incidence and survival. Am J Clin Nutr. 2012;96(4):855–63.
Guertin KA, Li XS, Graubard BI, Albanes D, Weinstein SJ, Goedert JJ, et al. Serum Trimethylamine N-oxide, Carnitine, Choline, and Betaine in Relation to Colorectal Cancer Risk in the Alpha Tocopherol, Beta Carotene Cancer Prevention Study. Cancer Epidemiol Biomarkers Prev. 2017;26(6):945–52.
Koochekpour S, Majumdar S, Azabdaftari G, Attwood K, Scioneaux R, Subramani D, et al. Serum glutamate levels correlate with Gleason score and glutamate blockade decreases proliferation, migration, and invasion and induces apoptosis in prostate cancer cells. Clin Cancer Res. 2012;18(21):5888–901.
Fazzari J, Lin H, Murphy C, Ungard R, Singh G. Inhibitors of glutamate release from breast cancer cells; new targets for cancer-induced bone-pain. Sci Rep. 2015;5:8380.
Lo M, Wang YZ, Gout PW. The x(c)- cystine/glutamate antiporter: a potential target for therapy of cancer and other diseases. J Cell Physiol. 2008;215(3):593–602.
Nilsson A, Haanstra JR, Engqvist M, Gerding A, Bakker BM, Klingmuller U, et al. Quantitative analysis of amino acid metabolism in liver cancer links glutamate excretion to nucleotide synthesis. Proc Natl Acad Sci U S A. 2020;117(19):10294–304.
Martin L, Gardner LB. Stress-induced inhibition of nonsense-mediated RNA decay regulates intracellular cystine transport and intracellular glutathione through regulation of the cystine/glutamate exchanger SLC7A11. Oncogene. 2015;34(32):4211–8.
Grimm C, Hofstetter G, Aust S, Mutz-Dehbalaie I, Bruch M, Heinze G, et al. Association of gamma-glutamyltransferase with severity of disease at diagnosis and prognosis of ovarian cancer. Br J Cancer. 2013;109(3):610–4.
Bansal A, Sanchez DJ, Nimgaonkar V, Sanchez D, Riscal R, Skuli N, et al. Gamma-Glutamyltransferase 1 Promotes Clear Cell Renal Cell Carcinoma Initiation and Progression. Mol Cancer Res. 2019;17(9):1881–92.
Kunutsor SK, Laukkanen JA. Gamma-glutamyltransferase and risk of prostate cancer: Findings from the KIHD prospective cohort study. Int J Cancer. 2017;140(4):818–24.
Laezza C, Pagano C, Navarra G, Pastorino O, Proto MC, Fiore D, et al. The Endocannabinoid System: a target for cancer treatment. Int J Mol Sci. 2020;21(3):747. https://doi.org/10.3390/ijms21030747.
Diaz-Laviada I. The endocannabinoid system in prostate cancer. Nat Rev Urol. 2011;8(10):553–61.
Wei Z, Liu X, Cheng C, Yu W, Yi P. Metabolism of Amino Acids in Cancer. Front Cell Dev Biol. 2020;8:603837.
Bader DA, McGuire SE. Tumour metabolism and its unique properties in prostate adenocarcinoma. Nat Rev Urol. 2020;17(4):214–31.
Geeraerts SL, Heylen E, De Keersmaecker K, Kampen KR. The ins and outs of serine and glycine metabolism in cancer. Nat Metab. 2021;3(2):131–41.
Amelio I, Cutruzzola F, Antonov A, Agostini M, Melino G. Serine and glycine metabolism in cancer. Trends Biochem Sci. 2014;39(4):191–8.
Wang D, Dubois RN. Eicosanoids and cancer. Nat Rev Cancer. 2010;10(3):181–93.
Panagiotopoulos AA, Kalyvianaki K, Castanas E, Kampa M. Eicosanoids in prostate cancer. Cancer Metastasis Rev. 2018;37(2-3):237–43.
Separovic D, Shields AF, Philip PA, Bielawski J, Bielawska A, Pierce JS, et al. Altered Levels of Serum Ceramide, Sphingosine and Sphingomyelin Are Associated with Colorectal Cancer: A Retrospective Pilot Study. Anticancer Res. 2017;37(3):1213–8.
Jiang Y, Tie C, Wang Y, Bian D, Liu M, Wang T, et al. Upregulation of Serum Sphingosine (d18:1)-1-P Potentially Contributes to Distinguish HCC Including AFP-Negative HCC From Cirrhosis. Front Oncol. 2020;10:1759.
Ogretmen B. Sphingolipid metabolism in cancer signalling and therapy. Nat Rev Cancer. 2018;18(1):33–50.
Kroll A, Cho HE, Kang MH. Antineoplastic Agents Targeting Sphingolipid Pathways. Front Oncol. 2020;10:833.
Wanders D, Hobson K, Ji X. Methionine Restriction and Cancer Biology. Nutrients. 2020;12:3.
Lu WC, Saha A, Yan W, Garrison K, Lamb C, Pandey R, et al. Enzyme-mediated depletion of serum l-Met abrogates prostate cancer growth via multiple mechanisms without evidence of systemic toxicity. Proc Natl Acad Sci U S A. 2020;117(23):13000–11.
Patil MD, Bhaumik J, Babykutty S, Banerjee UC, Fukumura D. Arginine dependence of tumor cells: targeting a chink in cancer's armor. Oncogene. 2016;35(38):4957–72.
Chen CL, Hsu SC, Ann DK, Yen Y, Kung HJ. Arginine Signaling and Cancer Metabolism. Cancers (Basel). 2021;13:14.
Kelavkar UP, Hutzley J, McHugh K, Allen KG, Parwani A. Prostate tumor growth can be modulated by dietarily targeting the 15-lipoxygenase-1 and cyclooxygenase-2 enzymes. Neoplasia. 2009;11(7):692–9.
Kelavkar UP, Nixon JB, Cohen C, Dillehay D, Eling TE, Badr KF. Overexpression of 15-lipoxygenase-1 in PC-3 human prostate cancer cells increases tumorigenesis. Carcinogenesis. 2001;22(11):1765–73.
Hada M, Edin ML, Hartge P, Lih FB, Wentzensen N, Zeldin DC, et al. Prediagnostic Serum Levels of Fatty Acid Metabolites and Risk of Ovarian Cancer in the Prostate, Lung, Colorectal, and Ovarian (PLCO) Cancer Screening Trial. Cancer Epidemiol Biomarkers Prev. 2019;28(1):189–97.
He L, Li H, Wu A, Peng Y, Shu G, Yin G. Functions of N6-methyladenosine and its role in cancer. Mol Cancer. 2019;18(1):176.
Wu Q, Xie X, Huang Y, Meng S, Li Y, Wang H, et al. N6-methyladenosine RNA methylation regulators contribute to the progression of prostate cancer. J Cancer. 2021;12(3):682–92.
Ji G, Huang C, He S, Gong Y, Song G, Li X, et al. Comprehensive analysis of m6A regulators prognostic value in prostate cancer. Aging (Albany NY). 2020;12(14):14863–84.
We appreciate all participants in the ATBC Study cohort for their contribution to the study.
The ATBC Study is supported by the Intramural Research Program of the U.S. National Cancer Institute, National Institutes of Health, Department of Health and Human Services. This work was funded by the Department of Defense award number W81XWH-18-1-0275.
Ethics approval and consent to participate
The study was approved by institutional review boards from both the Finnish National Institute for Health and Welfare and the U.S. National Cancer Institute (reference ID: OH95CN012). Written informed consent was provided by all participants.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
[Hazard Ratios and 95% Confidence Intervals for the Association Between Prostate Cancer Mortality and Prediagnostic Serum Metabolites Reaching the Bonferroni Corrected Threshold Based on 1,812 Prostate Cancer Cases in the ATBC Study, with further adjustment in the models]. Table S2. [Hazard Ratios and 95% Confidence Intervals for the Association Between Prostate Cancer Mortality and Prediagnostic Serum Metabolites Reaching the Bonferroni Corrected Threshold Stratified by Selected Factors in the ATBC Study]. Table S3. [Hazard Ratios and 95% Confidence Intervals for Serum Metabolites Achieving the Bonferroni Corrected Threshold Using Stepwise Cox Proportional Hazards Regression Analysis in the ATBC Study].
About this article
Cite this article
Huang, J., Zhao, B., Weinstein, S.J. et al. Metabolomic profile of prostate cancer-specific survival among 1812 Finnish men. BMC Med 20, 362 (2022). https://doi.org/10.1186/s12916-022-02561-4
- Metabolomic profile
- Prostate cancer survival
- Gamma-glutamyl amino acids
- Redox metabolites