Serum metabolomic profiling of prostate cancer risk in the prostate, lung, colorectal, and ovarian cancer screening trial.

BACKGROUND: Two recent metabolomic analyses found serum lipid, energy, and other metabolites related to aggressive prostate cancer risk up to 20 years prior to diagnosis.
METHODS: We conducted a serum metabolomic investigation of prostate cancer risk in the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial that included annual serum total prostate-specific antigen measurement and digital rectal examination. This nested study included 380 cases diagnosed post-screening and 380 controls individually matched to cases on age, race, study centre, and blood-collection date (median time to diagnosis, 10 years (range 4.4-17 years)). Sera were analysed on a high-resolution accurate mass platform of ultrahigh-performance liquid and gas chromatography/mass spectroscopy that identified 695 known metabolites. Logistic regression conditioned on the matching factors estimated odds ratios (OR) and 95% confidence intervals of risk associated with an 80th percentile increase in the log-metabolite signal.
RESULTS: Twenty-seven metabolites were associated with prostate cancer at P<0.05. Pyroglutamine, gamma-glutamylphenylalanine, phenylpyruvate, N-acetylcitrulline, and stearoylcarnitine showed the strongest metabolite-risk signals (ORs=0.53, 0.51, 0.46, 0.58, and 1.74, respectively; 0.001⩽P⩽0.006). Findings were similar for aggressive disease (peptide chemical class, P=0.03). None of the P-values were below the threshold of Bonferroni correction, however.
CONCLUSIONS: A unique metabolomic profile associated with post-screening prostate cancer is identified that differs from that in a previously studied, unscreened population.

RCT Entities:   Population Interventions Outcomes

Prostate cancer is the most frequently diagnosed malignancy and the third leading cause of cancer mortality after lung and colorectal cancer among men in developed populations (Torre ). Its aetiology remains largely unknown, with few established risk factors, other than age, race, family history of prostate cancer, and a number of low-penetrance genetic variants (Bostwick ; Chan ; Wolk, 2005; Kote-Jarai ; Barbieri ; Allott ). Recent advances in liquid or gas chromatography, mass spectrometry, and nuclear magnetic resonance methods have facilitated measurement of hundreds or thousands of low molecular weight, 80–1000 Da metabolites in biological samples such as serum, urine, and tissue. Measurement of the metabolome provides an integrated assessment of exogenous and endogenous exposures, and the host response to them, and has the potential to elucidate novel disease associations and biological mechanisms.

Two recent studies (Mondul , 2015) in the Alpha-Tocopherol, Beta-Carotene Cancer Prevention (ATBC) cohort identified metabolites related to aggressive prostate cancer in sera collected two decades prior to diagnosis. Lipid and energy metabolites including alpha-ketoglutarate (AKG), citrate, inositol-1-phosphate, several glycerophospholipids, and fatty acids, were lower in aggressive cases than in controls (Mondul ). The present investigation examines the prospective metabolomic profile of prostate cancer risk in a case–control study nested within the Prostate, Lung, Colorectal, and Ovarian Cancer Screening (PLCO) Trial.

Materials and methods

Study population

The PLCO trial was a large randomised-controlled trial designed to evaluate the efficacy of screening methods for prostate, lung, colorectal, or ovarian cancer. Participants without a history of these malignancies, ages 55 to 74 years, were enroled from 10 centres in the USA between 1993 and 2001, and were randomly assigned to either the screening or non-screening arm (Gohagan ; Prorok ; Hayes ). Blood samples were collected from participants in the screening arm only, between 0700 hours and 1600 hours without regard to the fasting status, and were aliquoted within 2 h and stored at −70 °C. The cancer screenings included chest x-ray, flexible sigmoidoscopy, and annual serum total prostate-specific antigen (PSA) measurement (for 6 years) and digital rectal examination (DRE; for 4 years). Men with an elevated PSA of ⩾4 ng ml−1 or a DRE suspicious for prostate cancer were referred to their medical-care providers for further prostate cancer work-up. Information for prostate cancer diagnoses based on abnormal screening results was retrieved from medical and pathology records. Additional prostate cancer diagnoses, including post-screening trial cases were identified through annual self-reported questionnaire responses and through the National Death Index, followed by medical record review and confirmation. All men (n=28 243) were followed from their initial prostate cancer screen, to the date of prostate cancer diagnosis, date of death (National Death Index), censor date (26 April 2012), or loss to follow-up, whichever occurred first.

Participants completed self-administered questionnaires at enrolment that included data regarding sociodemographics, height, weight, smoking behaviour, family history of cancer, and other diseases, physical activity, use of selected medications, and recent history of screening exams. All participants provided written informed consent, and the PLCO study was approved by the institutional review boards of the U.S. National Cancer Institute and the 10 PLCO screening centres.

Case identification and control selection

In order to increase time from blood collection to clinical diagnoses, 380 prostate cancer cases were randomly chosen from participants diagnosed during the post-screening trial period (4.4–17.0 years after baseline). On the basis of previous study findings (Mondul ;, 2015), aggressive cases were oversampled (n=298) and defined as those diagnosed with stage III or IV based on the tumour-node-metastasis staging system (Fleming ), or a biopsy Gleason score ⩾8. Cases with a Gleason score sum of 7 and stage II were considered to be of an intermediate level of aggressiveness (n=82). Because baseline blood samples had been depleted by many previous prostate cancer studies, we analysed serum collected at the first or second screening year visit. Using the method of incidence-density sampling without replacement, we randomly selected 380 controls that were free from any cancer at the time of case diagnosis, and individually matched them to cases by age (within 5 years), race, study centre, study year, and date (within 30 days) of blood collection.

Metabolite assays

Serum metabolomic profiling was conducted on a high-resolution accurate mass platform of ultrahigh-performance liquid chromatography/mass spectroscopy and gas chromatograph/mass spectroscopy (GC-MS; Metabolon, Durham, NC, USA). Workflow including extraction of raw data, peak-identification, and quality control (QC) processed on the assay platform, has been described (Evans ; Dehaven ). Metabolon measured values for a total of 722 identified metabolites. We excluded 27 metabolites for which >500 participants (66%) had missing values (below limit of detection), leaving 695 identified compounds for analysis. Metabolites were categorised as belonging to one of eight mutually exclusive chemical classes: amino acids and amino acid derivatives (subsequently referred to as ‘amino acids'), carbohydrates, cofactors and vitamins, energy metabolites, lipids, nucleotides, peptides, or xenobiotics. Each batch included blinded quality control samples (9%) from two male individuals. To assess the technical reliability of data, coefficients of variation (median=0.2, interquartile range=0.12–0.32), and intraclass correlation coefficients (median=0.87, interquartile range=0.64–0.98) were calculated as measures of QC. The cases and their matched controls were assayed within the same batches in order to avoid any effect of batch differences on the risk estimates.

Statistical analysis

Baseline characteristics of cases and controls were compared by either Wilcoxon rank sum or χ2-tests, for continuous or categorical variables, respectively. The signal of each metabolite was normalised within a given batch to standardise the batch variability, and for missing (below limit of detection) values, we imputed the minimum non-missing value. Metabolite values were log-transformed for analysis. Conditional logistic regression was used to examine odds ratios (ORs) and their 95% confidence intervals (CI), for the association between prostate cancer and each log-metabolite signal with an 80th percentile increase. Only matching factors were included in the final models. In addition to the matching factors, we performed sensitivity analyses to include BMI (<25 kg m−2, 25–30 kg m−2, or ⩾30 kg m−2) in the model, in addition to several other potential confounding factors. These included smoking status (never, former, or current), diabetes (yes or no), height (<175 cm, 175–180 cm, or ⩾180 cm), physical activity (<1 h per week, 1–3 h per week, or ⩾4 h per week), alcohol consumption (<0.05 drinks per day, 0.05–1 drinks per day, or ⩾1 drinks per day), processed meat consumption (<6.6 g per day, 6.6–16.8 g per day, or ⩾16.8 g per day), red meat consumption (<20.4 g per day, 20.4–43.7 g per day, or ⩾43.7 g per day), and total fat intake (<60.3 g per day, 60.3–88 g per day, or ⩾88 g per day). According to a Bonferroni correction for 695 tests, the threshold for statistical significance in our analysis is P=0.000072. However, this threshold is highly stringent due to the inter-correlations between many metabolites. We therefore also used principle component analysis (Jolliffe, 2005), and explored whether the grouped metabolites can distinguish case–control status. The top 10 principle components of metabolite measurements were calculated, and the same approach of conditional logistic regression (log-level with an 80th percentile increase) was applied to examine whether these components were associated with overall or aggressive prostate cancer. We also used a false-discovery rate (FDR) of 20% to define the significance threshold. We used gene-set analysis (GSA), a standard pathway analysis, to examine whether pre-defined metabolic super- and sub-pathways were associated with prostate cancer (Subramanian ).

We performed additional analyses restricting to aggressive or non-aggressive cases, non-Hispanic white men, black men, and men of other races, and stratified by age at enrolment (<65 vs 65+ years) and follow-up time (<10 years, ⩾10 years). The analysis by race was conducted, in part, for comparison with prior studies of Caucasian populations (Mondul , 2015).

Analyses were performed with SAS software version 9.3 (SAS Institute, Cary, NC, USA), and the GSA analysis was performed with the R statistical language version 3.2.3 (Vienna, Austria). All reported P-values are two-sided.

Results

Median time from blood collection to prostate cancer diagnosis was 10 years (inter-decile range=7.2–13.0 years) and median age at diagnosis was 72 years. Compared with controls, cases were more likely to have higher serum PSA and a family history of prostate cancer, and were more physically active at age 40 years, but were similar with respect other characteristics at study baseline (Table 1).

Table 1

Selected baseline characteristics of the cases and controls, PLCO studya

	Controls	Cases	P-value
N	380	380
Age at blood collection, years	64	64	Matched
Race, %
White, non-Hispanic	90.8	90.8
Black	4.5	4.5
Other races	4.7	4.7	Matched
PSA at blood draw (ng ml−1)	1.4	2.1	<0.0001
Height, cm	178.4 (2.7)	178.0 (2.6)	0.371
Weight, kg	87.0 (14.3)	87.0 (13.7)	0.930
BMI, kg m−2	27.4 (4.0)	27.5 (4.1)	0.659
Tobacco smoking status, %
Never smoker	37.9	43.2
Former smoker	56.3	52.4
Current smoker	5.8	4.5	0.29
Tobacco smoking duration,b years	25.1 (12.9)	23.9 (13.3)	0.272
Family history of prostate cancer, %	3.6	5.4	0.073
Diabetes, %	6.1	7.1	0.56
Physical activity at age 40 (h per week), %
4+ h (vs <3 h)	29.7	38.4	0.017

Abbreviations: PLCO=Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial; PSA=prostate-specific antigen; BMI=body mass index.

Values are means and standard deviations unless otherwise indicated.

Tobacco smoking duration among current or former smokers.

Metabolites related to risk of overall, aggressive and non-aggressive prostate cancer with a nominal P-value of <0.05 are shown in Tables 2, 3, 4 sorted by chemical class, sub-pathway and P-values. All metabolite associations with overall prostate cancer risk are presented in Supplementary Table 1. None of the P-values were below the threshold of correction for multiple comparisons or with a FDR<20%. The amino acids pyroglutamine (pGLU), phenylpyruvate, and N-acetylcitrulline, as well as the peptide gamma-glutamylphenylalanine, yielded the strongest signals, being inversely associated with overall prostate cancer (0.001⩽P<0.004, Table 2), whereas the acylcarnitine metabolite stearoylcarnitine was positively associated (OR=1.74, Table 2). Findings were similar for aggressive disease, with the peptide class of compounds being inversely associated (Table 3, P=0.027), but not with non-aggressive cancer. For non-aggressive disease, alpha-tocopherol, primary bile acid, and steroid hormone metabolites were well-represented and inversely associated with risk, as was cyclic AMP (Table 4). Stratification based on median time to diagnosis of aggressive disease revealed long-chain fatty acids, monohydroxy fatty acids, acylcarnitines, monoacylglycerol lipids, and lysolipids being positively associated among cases diagnosed within 10 years, and citrate, sphingolipid, steroid hormone, and glutathione metabolites associated among cases diagnosed 10–20 years after blood collection (data not shown).

Table 2

OR and 95% CI from conditional logistic regression of overall prostate cancer comparision of the 90th and 10th percentiles for serum metabolites (P<0.05) sorted by chemical class, sub-pathway, and P-valuea , b

Metabolite	Sub-pathway	OR	95% CI	P-value	P-value for chemical class
Amino acids and amino acids derivatives					0.20
Pyroglutamine	Glutamate metabolism	0.53	0.36, 0.78	0.0013
Phenylpyruvate	Phenylalanine and tyrosine metabolism	0.46	0.28, 0.75	0.0022
N-Acetylcitrulline	Urea cycle; arginine and proline metabolism	0.58	0.41, 0.84	0.0034
Homocitrulline	Urea cycle; arginine and proline metabolism	0.69	0.47, 1.00	0.0480
N-Acetylarginine	Urea cycle; arginine and proline metabolism	0.69	0.48, 1.00	0.0490
N-Acetyltryptophan	Tryptophan metabolism	0.67	0.49, 0.92	0.0122
N-Acetylkynurenine	Tryptophan metabolism	0.62	0.40, 0.95	0.0267
Tryptophan	Tryptophan metabolism	0.72	0.51, 1.00	0.0490
5-Methylthioadenosine	Polyamine metabolism	0.64	0.45, 0.93	0.0193
Trans-urocanate	Histidine metabolism	1.50	1.01, 2.23	0.0467
Cofactors and vitamins					0.20
Gulonic acid	Ascorbate and aldarate metabolism	0.73	0.56, 0.96	0.0239
Energy					0.29
Citrate	TCA cycle	1.54	1.04, 2.28	0.0294
Lipids					0.43
Stearoylcarnitine	Fatty acid metabolism (acyl carnitine)	1.74	1.17, 2.59	0.0063
Palmitoylcarnitine	Fatty acid metabolism (acyl carnitine)	1.49	1.02, 2.16	0.0389
15-Methylpalmitate	Fatty acid, branched	1.63	1.12, 2.38	0.0112
17-Methylstearate	Fatty acid, branched	1.47	1.03, 2.11	0.0354
Hexadecanedioate	Fatty acid, dicarboxylate	1.45	1.07, 1.98	0.0167
Caproate (6:0)	Medium-chain fatty acid	1.44	1.05, 1.98	0.0248
1-Palmitoleoyl-2-linoleoyl-GPC (16:1/18:2)	Phospholipid metabolism	0.60	0.40, 0.88	0.0096
1-Palmitoyl-2-arachidonoyl-GPI (16:0/20:4)	Phospholipid metabolism	1.64	1.09, 2.47	0.0187
Trimethylamine-N-oxide	Phospholipid metabolism	0.68	0.47, 0.97	0.0327
1-Linoleoyl-2-arachidonoyl-GPC (18:2/20:4n6)	Phospholipid metabolism	1.48	1.01, 2.16	0.0425
Tauro-beta-muricholate	Primary bile acid metabolism	0.68	0.48, 0.97	0.0329
1-Arachidonoyl-GPC (20:4n6)	Lysolipid	1.46	1.01, 2.10	0.0457
Nucleotide					0.61
2′-Deoxyuridine	Pyrimidine metabolism	1.47	1.07, 2.03	0.0191
Peptides					0.092
Gamma-glutamylphenylalanine	Gamma-glutamyl amino acid	0.51	0.34, 0.78	0.0015
Leucylglutamine	Dipeptide	0.67	0.47, 0.97	0.0356

Abbreviations: CI=confidence interval; GPC=glycerophosphocholine; GPI=glycerophosphoinositol; OR=odds ratio; TCA=tricarboxylic acid cycle.

The odds ratio per 80th percentile increase in metabolite. On the basis of 380 case–control pairs.

All metabolites had detectable values in >90% of the study population with the exception of tauro-beta-muricholate and leucylglutamine (41% and 61%, respectively).

Table 3

OR and 95% CI from conditional logistic regression of aggressive prostate cancer comparision of the 90th and 10th percentiles for serum metabolites (P<0.05) sorted by chemical class, sub-pathway, and P-valuea , b

Metabolite	Sub-pathway	OR	95% CI	P-value	P-value for chemical class
Amino acids and amino acids derivatives					0.17
Pyroglutamine	Glutamate metabolism	0.48	0.31, 0.74	0.0011
Phenylpyruvate	Phenylalanine and tyrosine metabolism	0.41	0.23, 0.75	0.0035
Thyroxine	Phenylalanine and tyrosine metabolism	0.70	0.50, 0.99	0.0423
Tryptophan	Tryptophan metabolism	0.61	0.41, 0.90	0.0124
N-Acetylkynurenine	Tryptophan metabolism	0.55	0.34, 0.88	0.0137
N-Acetyl-cadaverine	Lysine metabolism	1.72	1.10, 2.69	0.0173
Glutarate (pentanedioate)	Lysine metabolism	1.68	1.07, 2.65	0.0240
N-Acetylcitrulline	Urea cycle; arginine and proline metabolism	0.62	0.41, 0.92	0.0182
Arginine	Urea cycle; arginine and proline metabolism	0.63	0.42, 0.96	0.0308
N-Acetylarginine	Urea cycle; arginine and proline metabolism	0.64	0.42, 0.98	0.0389
5-Methylthioadenosine	Polyamine metabolism	0.63	0.41, 0.96	0.0308
Trans-urocanate	Histidine metabolism	1.63	1.02, 2.60	0.0410
Lipids					0.17
Hexadecanedioate	Fatty acid, dicarboxylate	1.67	1.16, 2.42	0.0063
Tetradecanedioate	Fatty acid, dicarboxylate	1.53	1.05, 2.25	0.0278
Stearoylcarnitine	Fatty acid metabolism (acyl carnitine)	1.75	1.11, 2.76	0.0159
Palmitoylcarnitine	Fatty acid metabolism (acyl carnitine)	1.63	1.06, 2.52	0.0265
15-Methylpalmitate	Fatty acid, branched	1.69	1.10, 2.58	0.0163
17-Methylstearate	Fatty acid, branched	1.62	1.07, 2.44	0.0225
5Alpha-pregnan-3beta,20beta-diol monosulfate	Steroid	1.57	1.04, 2.37	0.0330
1-Oleoylglycerol (18:1)	Monoacylglycerol	1.58	1.03, 2.43	0.0375
1-Palmitoyl-2-arachidonoyl-GPI (16:0/20:4)	Phospholipid metabolism	1.65	1.03, 2.63	0.0376
1-Stearoyl-2-oleoyl-GPC (18:0/18:1)	Phospholipid metabolism	1.57	1.00, 2.45	0.0496
Eicosapentaenoate (20:5n3)	Polyunsaturated fatty acid (n3 and n6)	1.53	1.02, 2.30	0.0392
Nucleotide					0.81
2′-Deoxyuridine	Pyrimidine metabolism	1.50	1.05, 2.15	0.0267
Peptides					0.027
Leucylglutamine	Dipeptide	0.58	0.38, 0.88	0.0098
Gamma-glutamylphenylalanine	Gamma-glutamyl amino acid	0.57	0.36, 0.90	0.0152
Gamma-glutamylglycine	Gamma-glutamyl amino acid	0.64	0.42, 0.99	0.0430
Xenobiotics					0.58
7-Methylurate	Xanthine metabolism	0.60	0.37, 0.99	0.0462

Abbreviations: CI=confidence interval; GPC=glycerophosphocholine; GPI=glycerophosphoinositol; OR=odds ratio.

The odds ratio per 80th percentile increase in metabolite. On the basis of 298 case–control pairs.

All metabolites had detectable values in >90% of the study population with the exception of leucylglutamine and 7-methylurate (61% and 68%, respectively).

Table 4

OR and 95% CI from conditional logistic regression of non-aggressive prostate cancer comparision of the 90th and 10th percentiles for serum metabolites (P<0.05) sorted by chemical class, sub-pathway, and P-valuea , b

Metabolite	Sub-pathway	OR	95% CI	P-value	P-value for chemical class
Amino acids and amino acids derivatives					0.69
2-Aminoadipate	Lysine metabolism	0.63	0.41, 0.96	0.0311
Pipecolate	Lysine metabolism	0.68	0.48, 0.97	0.0332
3-Indoxyl sulfate	Tryptophan metabolism	0.70	0.49, 1.00	0.0483
Cofactors and vitamins					0.18
Alpha-CEHC glucuronide	Tocopherol metabolism	0.63	0.44, 0.90	0.0105
Gamma-tocopherol/beta-tocopherol	Tocopherol metabolism	1.51	1.07, 2.14	0.0196
Alpha-CEHC sulfate	Tocopherol metabolism	0.68	0.49, 0.95	0.0242
Alpha-tocopherol	Tocopherol metabolism	0.66	0.45, 0.96	0.0278
Lipids					0.32
1-Palmitoleoyl-2-linoleoyl-GPC (16:1/18:2)	Phospholipid metabolism	0.58	0.40, 0.86	0.0060
1-Palmitoyl-2-linolenoyl-GPC (16:0/18:3)	Phospholipid metabolism	1.62	1.05, 2.49	0.0293
1-Linoleoyl-2-arachidonoyl-GPC (18:2/20:4n6)	Phospholipid metabolism	1.50	1.04, 2.16	0.0303
Taurocholate	Primary bile acid metabolism	0.60	0.40, 0.88	0.0094
Glycocholate	Primary bile acid metabolism	0.69	0.49, 0.96	0.0272
4-Androsten-3beta,17beta-diol monosulfate	Steroid	0.62	0.43, 0.90	0.0117
5Alpha-androstan-3alpha,17beta-diol monosulfate	Steroid	0.61	0.40, 0.92	0.0187
Dehydroepiandrosterone sulfate	Steroid	0.64	0.44, 0.93	0.0204
17Alpha-hydroxypregnenolone sulfate	Steroid	0.70	0.51, 0.97	0.0303
4-Androsten-3beta,17beta-diol disulfate	Steroid	0.72	0.52, 0.99	0.0431
4-Androsten-3alpha,17alpha-diol monosulfate	Steroid	0.68	0.47, 0.99	0.0445
1-(1-Enyl-stearoyl)-2-arachidonoyl-GPE (P-18:0/20:4)	Plasmalogen	1.49	1.06, 2.10	0.0229
2-Hydroxyadipate	Fatty acid, dicarboxylate	0.66	0.45, 0.96	0.028
2-Stearoyl-GPE (18:0)	Lysolipid	0.67	0.46, 0.97	0.0347
Nucleotide					0.73
Adenosine 3′,5′-cyclic monophosphate	Purine metabolism, adenine containing	0.56	0.34, 0.92	0.0220
Peptides					0.29
Threonylphenylalanine	Dipeptide	1.47	1.05, 2.06	0.0256
Gamma-glutamylphenylalanine	Gamma-glutamyl amino acid	0.65	0.44, 0.96	0.0315
Gamma-glutamylglutamine	Gamma-glutamyl amino acid	1.41	1.03, 1.94	0.0339
Xenobiotics					0.51
1,3-Dimethylurate	Xanthine metabolism	0.69	0.48, 0.99	0.0426

Abbreviations: Alpha-CEHC=Alpha-carboxyethyl-hydroxychroman; CI=confidence interval; GPC=glycerophosphocholine; GPE=glycerophosphoethanolamine; OR=odds ratio.

The odds ratio per 80th percentile increase in metabolite. On the basis of 82 case–control pairs.

All metabolites had detectable values in >90% of the study population with the exception of alpha-CEHC glucuronide, alpha-CEHC sulfate, 1-palmitoyl-2-linolenoyl-GPC (16:0/18:3), 17alpha-hydroxypregnenolone sulfate, 2-hydroxyadipate, adenosine 3′,5′-cyclic monophosphate, and 1,3-dimethylurate (38%, 58%, 66%, 59%, 61%, 54% and 84%, respectively).

These findings were not materially changed when adjusted for BMI (Supplementary Table 2) or for BMI plus several additional factors including smoking status, diabetes, height, physical activity, alcohol consumption, processed and red meat consumption, and total fat intake (Supplementary Table 3). These findings were also not substantially altered when adjusted for serum total PSA (data not shown), and were essentially similar in non-Hispanic white, black, and other races combined groups, although with some differences in the latter two populations that are probably due to small sample sizes (Supplementary Tables 4–6). We found no material interactions between the top 27 metabolites and other risk factors, including age, BMI, and smoking (data not shown).

Although none of the top 10 principle components were significantly associated with overall or aggressive prostate cancer, with all tests having P-values>0.025 (P<0.005 was the significant level), the metabolite sub-pathway analysis revealed that branched and medium-chain fatty acid metabolites, as well as the tocopherol metabolites, were inversely associated with overall prostate cancer (P=0.016, 0.032, and 0.017, respectively, Table 5). Aggressive disease was positively associated with the branched-chain fatty acids P=0.015), and inversely related to the tryptophan and urea cycle/arginine/proline metabolites (P=0.037 and 0.046, respectively). The tocopherol and primary bile acid metabolite subclasses were inversely related to non-aggressive cancer (P=0.014 and 0.021, respectively, Table 5).

Table 5

GSA for sub-pathway of serum metabolites and prostate cancer (P<0.05)

Sub-pathway	No. of contributing metabolites	P-value from GSA
Overall prostate cancer
Fatty acid, branched	2	0.016
Tocopherol metabolism	6	0.017
Medium-chain fatty acid	7	0.032
Aggressive prostate cancer
Fatty acid, branched	2	0.015
Tryptophan metabolism	20	0.037
Urea cycle; arginine and proline metabolism	17	0.046
Non-aggressive prostate cancer
Tocopherol metabolism	6	0.014
Primary bile acid metabolism	10	0.021
Overall prostate cancer in Caucasians
Tocopherol metabolism	6	0.007
Fatty acid, branched	2	0.009
Medium-chain fatty acid	7	0.019

Abbreviation: GSA=gene-set analysis.

Among the top metabolites related to prostate cancer in our earlier investigation (Mondul ) that served as our a priori hypotheses for the present analysis, only three associations with overall prostate cancer replicated in the present data: the lipids 1-palmitoleoyl-2-linoleoyl-GPC (16:1/18:2) (OR=0.60, 95% CI=0.40–0.88, P=0.0096) and tauro-beta-muricholate (OR=0.68, 95% CI: 0.48–0.97, P=0.033), and the nucleotide 2′-deoxyuridine (OR=1.47, 95% CI: 1.07–2.03, P=0.019; Table 2). By contrast, the inverse association between aggressive prostate cancer, and the energy metabolites, (Mondul ) alpha-AKG or citrate, did not replicate, and appeared to have opposite associations (AKG OR=1.45, P=0.098, and citrate OR=1.48, P=0.077). It is noteworthy that these positive associations appeared stronger among cases diagnosed 10–20 years after blood collection (AKG OR=1.83, 95% CI 0.99–3.38, P=0.053; citrate OR=2.04, 95% CI 1.06–3.91, P=0.032) compared with those diagnosed within 10 years (OR=1.11, P=0.76 for AKG, and OR=1.13, P=0.68 for citrate). The present data also did not replicate the previous positive associations with aggressive prostate cancer for phenylpyruvate, thyroxine, and arginine (Table 3), and trimethylamine-N-oxide (data not shown).

Discussion

Top metabolite signals for overall prostate cancer risk observed here were pGLU, gamma-glutamylphenylalanine, phenylpyruvate, N-acetylcitrulline, and stearoylcarnitine. Peptide metabolites were associated with aggressive disease, whereas most of the previously reported lipid and energy metabolites were not. These profiles appeared similar across racial/ethnic groups, although our sample size was small for non-white men. By contrast, risk of non-aggressive disease was related to tocopherols, sex steroids, and primary bile acids.

Pyroglutamine is a key amino acid residue for thyrotrophin-releasing hormone (TRH) and TRH-like peptides (Huber ) that promote fertilisation by enhancing sperm capacitation (Cockle ; Huber ), act as paracrine factors influencing prostate cell growth and differentiation (Bilek ), and stimulate thyroid hormone secretion (Bilek ; Bilek, 2000; Maran ). Gamma-glutamylphenylalanine and gamma-glutamylglycine peptides were also inversely associated with aggressive prostate cancer, similar to the gamma-glutamylhistidine finding in the ATBC Study (Mondul ). These metabolites are glutathione (GSH) degradation products of gamma-glutamyl transpeptidase activity (Zhang ) that regenerates intracellular GSH (Wu ; Zhang ) for detoxification of reactive oxygen species (Wu ; Maher, 2005; Zhang ). Several N-acetyl amino acids including N-acetylcitrulline, N-acetylarginine, N-acetyltryptophan, and N-acetylkynurenine were also reduced in men later diagnosed with prostate cancer, possibly indicating increased post-translational protein N-terminal acetylation (Aksnes ) and aminoacetylase activity, the latter having been related to increased cell proliferation and carcinogenesis (Aksnes ), and histone-chromatin gene regulation (Eberharter and Becker, 2002). Also inversely associated were tryptophan metabolites that influence the role of inflammation in cancer progression (Prendergast ; Santhanam ). Tryptophan and kynurenine pathway metabolites have been previously associated with several, particularly gastrointestinal, cancers (Uyttenhove ; Witkiewicz ; Liu ; Koblish ; Balachandran ; Zhang ;, 2011b; Ferdinande ; Zhang ).

The proliferative anabolic state of prostate adenocarcinomas requires active biosynthesis of amino acids, nucleic acids, peptides, and lipids (Gang ), with dynamic metabolomic phenotype stages being likely. Aberrations in lipid metabolites and especially fatty acid metabolism have been extensively studied with respect to prostate carcinogenesis (Harvei ; Norrish ; Mannisto ; Chavarro ; Crowe ; Bassett ; Crowe ; Wu ). De novo fatty acid biosynthesis is upregulated in aggressive cancers (Nomura , 2011; Carracedo ), and the present data appear consistent with this only with respect to branched-chain fatty acids. The lower serum lipid metabolomic profile previously demonstrated for aggressive prostate cancer risk encompassed not only fatty acids, presumably to provide energy for cell proliferation and tumour growth, but also a large number of lysolipids, inositols, and sphingomyelins that implicated cell membrane biosynthesis and signalling (Kuemmerle ). Fatty acids rather than the latter compounds were more prominent in the present study and positively associated with risk, possibly indicating triglyceride mobilisation from adipose tissue and increased fatty acid synthase activity leading to increased acylcarnitine-mediated long-chain fatty acid transport into the mitochondria matrix for β-oxidation and ATP synthesis (Carter ; Crowe ). Despite laboratory evidence that eicosapentaenoic acid can inhibit prostate tumourigenesis (Cave, 1991; Rose, 1997; Astorg, 2004; Larsson ; Kobayashi ), we found this omega-3 fatty acid to be related to higher risk of aggressive prostate cancer, which is consistent with two studies (Crowe , 2014).

Our findings regarding non-aggressive prostate cancer are also of interest and differ substantially from those for aggressive disease. Notably, we detected an inverse risk association signal for alpha-tocopherol metabolites, consistent with both significantly lower prostate cancer incidence in one low-dose trial (Heinonen ) and a recent pooled analysis of 15 international cohorts (Key ). Several sex steroids also showed inverse risk associations with non-aggressive prostate cancer, including androgen precursors (e.g., DHEA-sulfate) and metabolites (e.g., 4-androstene-3beta,17 diol-sulfate). (Testosterone was also inversely associated, but more weakly (OR=0.76; P=0.11)). By contrast, only one sex steroid was related to aggressive disease and in a positive direction, more in line with the role androgens are thought to have to promote prostate carcinogenesis. Whether this stage difference represents true underlying biological variation in the influence of androgens on early vs aggressive disease will require further study.

Two other findings of inverse associations with non-aggressive disease risk deserve mention. 2-Aminoadipate (as well as pipecolate) is involved in lysine metabolism, and has been found elevated in malignant (T2 and T3) vs adjacent normal prostate tissue and related to earlier disease recurrence (Jung ). One of our strongest signals, cyclic AMP, is synthesised from ATP by plasma membrane-associated, hormonally responsive adenylate cyclase, and serves as an intracellular second messenger that mediates several hormones.

As compared with data from the ATBC Study that showed inverse prostate cancer associations with the energy and lipid compounds, including AKG, citrate, inositol-1-phosphate, glycerophospholipids, and fatty acids (Mondul ), the present findings portray a very different metabolite-risk profile that features primarily amino acids and peptides. Additional prospective investigations will be needed to interpret the two distinct patterns, but several differences in study designs and populations are likely responsible. Importantly, cases in PLCO were diagnosed within 15 years post-trial following up to four to six prostate screening examinations (i.e., DRE and serum PSA), whereas the cases in ATBC were clinically diagnosed within 20 years of serum collection without systematic PSA and DRE screening. Mean baseline total PSA concentrations of cases were 2 ng ml−1 and 8 ng ml−1 in PLCO and ATBC, respectively, and even the PLCO cases diagnosed with aggressive disease had a mean baseline PSA concentration of only 2 ng ml−1. These and other study differences likely resulted in examination of developmentally distinct metabolic phases of growing prostate adenocarcinomas related to primary tumour size and disease extent.

Time of serum collection and fasting status also differed markedly between the studies; that is, between 0700 hours and 1600 hours in PLCO without regard to time of last meal, vs 0800 hours to 1100 hours after an overnight fast in ATBC. Although the effects of non-fasting status on most of the circulating metabolites we examined have not been rigorously evaluated, some previous studies indicate minimal changes in PSA (Tuncel ), lipids (Langsted ; Mora ), and amino acids such as homocysteine (Fokkema ), and one report suggests that fasting status is not an important source of variability in measurements of over 200 metabolites (Townsend ). The different populations, US vs Finland, as well as the fact that ATBC included only smokers of at least five cigarettes per day at study entry may also have influenced the divergent findings.

The major strengths of this investigation include use of serum samples that were collected up to two decades prior to prostate cancer diagnosis and good laboratory reproducibility for several hundred metabolites. Our sample size was moderate, albeit the largest to date for a prospective study of the circulating metabolome and risk of prostate cancer. The impact of the PLCO prostate cancer screening on the stages of cases diagnosed and the non-fasting serologic status may have limited the robustness and sensitivity of replicating the previous ATBC findings. In addition, although we could not rule out the possibility of unmeasured confounding factors explaining some of the present findings, including some related to selection of trial post-screening cases, the materially unchanged results from sensitivity analyses that adjusted for several potential risk factors argues against such biases.

In conclusion, the present study demonstrates a unique metabolomic profile associated with post-screening prostate cancer in the PLCO cohort. The inverse associations for amino acid and peptide metabolites, and positive associations for lipids, contrast with the energy/lipid profile previously reported (Mondul , 2015). This is likely the result of studying prostate cancer cases with serum samples collected at an earlier point in the natural history of the disease due to the highly screened nature of the study population, rather than more clinically advanced, yet still undiagnosed, malignancies. The two distinct metabolite profiles may represent molecular species influencing risk of prostate cancers at different points in their development, including early initiating or tumour promoting metabolic states. Whether the inverse association with pGLU is related to a role for thyrotropin status of early prostate cancer will require further study. The metabolomic profiles identified to date should be re-examined in additional prospective analyses of prostate cancer in larger studies and consortia.