Metabolic imaging across scales reveals distinct prostate cancer phenotypes =========================================================================== * Nikita Sushentsev * Gregory Hamm * Lucy Flint * Daniel Birtles * Aleksandr Zakirov * Jack Richings * Stephanie Ling * Jennifer Y. Tan * Mary A. McLean * Vinay Ayyappan * Ines Horvat Menih * Cara Brodie * Jodi L. Miller * Ian G. Mills * Vincent J. Gnanapragasam * Anne Y. Warren * Simon T. Barry * Richard J.A. Goodwin * Tristan Barrett * Ferdia A. Gallagher ## Abstract Hyperpolarised magnetic resonance imaging (HP-13C-MRI) has shown promise as a clinical tool for detecting and characterising prostate cancer. Here we have used a range of spatially resolved histological techniques to identify the biological mechanisms underpinning differential [1-13C]lactate labelling between benign and malignant prostate, as well as tumours containing cribriform and non-cribriform Gleason pattern 4 disease. The elevated hyperpolarised [1-13C]lactate signal in prostate cancer compared to the benign prostate is primarily driven by increased tumour epithelial cell density and vascularity, rather than differences in epithelial lactate concentration between tumour and normal. We also demonstrate that tumours of the cribriform subtype may lack [1-13C]lactate labelling, which is explained by their lower epithelial lactate dehydrogenase expression, higher mitochondrial pyruvate carrier density, and increased lipid abundance compared to lactate-rich non-cribriform lesions. These findings highlight the potential of combining spatial metabolic imaging tools across scales to identify novel metabolic phenotypes in prostate cancer. ## Introduction Hyperpolarised [1-13C]pyruvate MRI (HP-13C-MRI) is an emerging non-ionising metabolic imaging technique to probe tumour metabolism1. The method offers several advantages over routine [18F]-2-fluoro-2-deoxy-D-glucose ([18F]FDG) positron emission tomography (PET) including discrimination of the injected substrate from its metabolic products to quantify cellular metabolism within specific tumour compartments such as the epithelium2, as well as sensitivity to lesions with low glucose uptake including prostate cancer (PCa)3. Several clinical reports have shown the potential of HP-13C-MRI to detect PCa, assess disease aggressiveness, and detect tumour response to therapy2–7. Importantly, some of these studies have attempted to explain the mechanism of the observed tumour [1-13C]lactate labelling patterns by correlating imaging data with tissue-based metabolic biomarkers2,8–10. However, a more comprehensive biological validation of HP-13C-MRI is required to assess its translational potential and navigate its use in specific clinical scenarios11. In addition to PCa detection by differentiating it from the healthy prostate, an area of particular unmet need is the non-invasive phenotyping of intermediate-risk disease12, which includes tumours comprised of Gleason pattern 4 (GP4) glands of cribriform and non-cribriform morphology. In 2019, both the International Society of Urological Pathology (ISUP) and Genitourinary Pathology Society (GUPS) recommended the reporting of invasive cribriform carcinoma (ICC) in biopsy specimens due to its higher recurrence rate and increased PCa-specific mortality compared to non-cribriform lesions13,14. While recent studies in this area have mostly focused on transcriptomic features of ICC15–17, little is known about its metabolic properties in contrast to other PCa phenotypes. In addition, the appearances of this tumour variant on HP-13C-MRI are unknown, which is of high translational relevance since up to 80% of pure cribriform tumours can be missed on conventional MRI18. The current understanding of the biological determinants of [1-13C]lactate labelling detected with HP-13C-MRI is that it reflects a complex interplay between biological factors which are summarised in **Fig. 1**. The primary determinant is hyperpolarised [1-13C]pyruvate delivery to the tissue of interest, which is a function of perfusion, microvascular density, and endothelial permeability19–21. The second major factor is the cellular capacity for [1-13C]pyruvate uptake, which is primarily mediated by the monocarboxylate transporter 1 (MCT1)4,22–24. The third factor relates to the intracellular metabolic fate of [1-13C]pyruvate, which is a function of its enzymatic exchange into [1-13C]lactate, and mitochondrial flux as measured through the formation of 13C-bicarbonate25. The fourth factor represents the natural abundance of the endogenous intracellular lactate pool which accepts the 13C label from the imported [1-13C]pyruvate26: previous clinical studies have inferred this indirectly from lactate dehydrogenase (LDH) expression2,5 but it can also be measured directly using tissue-based spatial metabolomics within specific cellular regions of the tumour27. Finally, the fifth factor is the capacity of the tissue-of-interest to generate a sufficient [1-13C]lactate signal-to-noise ratio (SNR) for detection with HP-13C-MRI, which in part is influenced by the density of MCT1-expressing and lactate-abundant cells. ![Fig. 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/10/02/2023.09.29.23296228/F1.medium.gif) [Fig. 1:](http://medrxiv.org/content/early/2023/10/02/2023.09.29.23296228/F1) Fig. 1: Study design. This clinical study included two prospective surgical cohorts of PCa patients whose imaging data and surgical specimens were analysed to measure five biological factors that can influence clinical [1-13C]lactate labelling. In the Hyperpolarised MRI cohort (n = 8), quantitative 1H-MRI, IHC, and RNAscope data were used to infer tissue delivery, epithelial uptake, and the intracellular metabolic fate of [1-13C]pyruvate in 15 tumour and 15 benign areas. This was complemented by spatial metabolomic analysis of a set of histologically matched fresh-frozen benign (n = 61) and tumour (n = 56) samples from a DESI-MSI cohort (n = 13), which enabled us to assess the endogenous lactate pool as a measure of the cellular capacity for [1-13C]labelling. Digital pathology data from both cohorts were used to quantify the density of epithelial cells as a measure of tissue capacity for generating detectable [1-13C]lactate signal. Tumour characteristics and methods are detailed in **Table 1** and **Methods**. View this table: [Table 1.](http://medrxiv.org/content/early/2023/10/02/2023.09.29.23296228/T1) Table 1. Key histopathological characteristics of tumours included in this study. a These lesions were excluded from the tissue-based analysis due to interval ADT prior to RARP2. b These lesions were HP-13C-MRI occult, with [1-13C]lactate SNR <5.0 in both cases. Here, we applied this five-factor biological framework to prospectively collect imaging and spatially resolved multi-modal data from two matched intermediate-risk PCa patient cohorts. We show that tumour [1-13C]lactate labelling is primarily driven by the increased epithelial cell density and vascularity of PCa areas compared to the benign prostate, while the endogenous epithelial lactate concentration is similar between the two tissue types. However, we also demonstrate that some cribriform tumours lack [1-13C]lactate signal on HP-13C-MRI despite sufficient cellularity and vascularity, a finding potentially explained by their lower LDH expression and higher mitochondrial pyruvate flux compared to equally cellular and well-perfused non-cribriform lesions. Using spatial metabolomics acquired from mass spectrometry imaging (MSI), we then show that ICC epithelium has a low endogenous lactate pool in the presence of abundant unsaturated fatty acids, presenting a distinct metabolic phenotype compared to lactate-rich non-cribriform GP4 glands. These findings indicate that although [1-13C]lactate labelling may be a powerful tool for non-invasively imaging tumour metabolism, it may provide negative contrast for ICC detection, representing a limitation for characterisation of intermediate-risk PCa. This study also highlights the potential of tissue-based spatial metabolomics to guide future developments in clinical metabolic imaging techniques that are specific to the metabolic features of selected tumours. ## Results As described in **Fig. 1**, this study included two prospective cohorts of PCa patients who underwent robot-assisted radical prostatectomy (RARP) in our centre. The HP-13C-MRI cohort included 8 patients who underwent successful HP-13C-MRI prior to RARP and represented a subset of a previously described patient population2. Histopathological examination of formalin-fixed paraffin-embedded (FFPE) whole-mount slides revealed the presence of 15 lesions, and the detailed morphological characteristics of these are presented in **Table 1**. To complement the use of tissue-based biomarkers for inferring cellular metabolic phenotype, we sought to directly assess endogenous epithelial metabolite abundance by means of spatially resolved desorption electrospray ionisation mass spectrometry imaging (DESI-MSI). To avoid the potentially detrimental effect of using FFPE with DESI-MSI28,29, fresh-frozen RARP samples from a cohort of 13 patients were used for this analysis (spatial metabolomics cohort in **Fig. 1**). A total of 117 tissue cores were obtained, of which 61 and 56 were derived from 15 tumours and 15 benign tissue areas, respectively. For comparability, tumours from the two cohorts were matched for the key histopathological parameters presented in **Table 1**. ### Increased tumour epithelial cell density, not lactate abundance, explains the ability of HP-13C-MRI to distinguish PCa from the benign prostate Our first step in the biological validation of HP-13C-MRI was aimed at dissecting the specific mechanisms behind its ability to detect PCa by distinguishing it from the healthy prostate. In the HP-13C-MRI cohort, both hyperpolarised [1-13C]pyruvate and [1-13C]lactate were measured in histopathologically-proven tumour areas (n = 13) (*P* < 0.001 for both; **Fig. 2b**), with results in line with similar findings described in prior clinical reports3,4,30. A similar trend was noted in the DESI-MSI cohort, where the whole-core measurements of absolute lactate abundance were significantly higher in tumour samples compared to benign specimens (*P* < 0.0001; **Fig. 2c**), in keeping with prior studies using bulk MSI31,32. However, these results can be partially explained by tissue density, which is known to increase in PCa compared to the benign tissue due to the obstruction of empty glandular spaces by neoplastic epithelium33. Hence, when the whole-core lactate abundance was corrected for tissue density, which was significantly higher in tumour samples (*P* = 0.007; **Fig. 2c**), the difference in endogenous lactate between the two tissue types reduced but remained significant (*P* = 0.043; **Fig. 2c**). ![Fig. 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/10/02/2023.09.29.23296228/F2.medium.gif) [Fig. 2:](http://medrxiv.org/content/early/2023/10/02/2023.09.29.23296228/F2) Fig. 2: Biological validation of the ability of HP-13C-MRI to detect PCa by differentiating it from the healthy prostate. **a,** Representative whole-mount H&E section of a surgical FFPE sample obtained from a patient who harboured a single focus of ISUP GG2 PCa (red outline) that showed increased [1-13C]pyruvate and [1-13C]lactate signal on corresponding HP-13C-MRI maps compared to contralateral benign prostate (green outline). **b,** Plots comparing [1-13C]pyruvate and [1-13C]lactate SNR derived from the ROIs encompassing areas of HP-13C-MRI-visible PCa (n = 13) and contralateral benign prostate (n = 13) from the hyperpolarised MRI cohort. **c,** Representative fused H&E and DESI-MSI lactate maps of benign and tumour FF cores from the spatial metabolomics cohort. Plots on the left compare the absolute and tissue-density-corrected whole-core DESI-MSI derived lactate abundance between the benign (n = 61) and tumour (n = 38) cores; plots on the right compare the whole-core tissue density between the two specimen types. **d,** DenseNet tissue classifier outputs overlaid on the H&E images of representative benign and tumour cores, with the below plots comparing epithelial and stromal cell fractions both within and across the benign and tumour cores. **e,** H&E maps of representative benign and tumour cores with manually segmented areas of benign epithelium (green), tumour epithelium (red), and stroma (blue) used to derive DESI-MSI measured of endogenous lactate abundance compared in the below plots between and across benign and tumour cores. **f-l,** Plots and representative images comparing epithelial cell density (**f**), 1H-MRI-derived ADC (**g**), CD31 density (**i**), 1H-MRI-derived *K*trans (**j**), and epithelial MCT1 density (**l**) between the benign (n = 11) and tumour (n = 11) areas from the hyperpolarised MRI cohorts. **h** and **k** present Spearman’s correlation plots comparing ADC values with epithelial cell density (**h**), as well as CD31 with *K*trans (**k**) derived from both benign and tumour areas. **m,** Plot comparing TCGA-PRAD-derived *SLC16A1* mRNA expression between the benign and tumour prostatectomy samples. All plots are scatterplots with bars (lines are median values, bars are interquartile ranges) with *P* derived using the Mann-Whitney *U* test. Scale bars denote 50 μm. Given the established metabolic compartmentalisation within the human prostate2,34–36, it was then important to identify the cellular source for the increased lactate abundance in PCa samples. To do so, we first trained a DenseNet model to segment each tissue section into stromal and epithelial compartments, allowing us to quantify their proportions in our specimens. As expected in this age group37, benign cores were mostly comprised of stromal tissue, whereas tumour samples were predominantly epithelial (*P* < 0.005 for both; **Fig. 2d**). While the epithelial lactate abundance was similar between benign and tumour cores (*P* = 0.54; **Fig. 2e**), it was the significant increase in the stromal lactate (*P* < 0.0001; **Fig. 2e**) that explained the difference in whole-core lactate measurements observed in **Fig. 2c**. However, despite the high concentration of lactate in the stroma, this is unlikely to make a significant contribution to the imaged HP-13C-MRI signal due to negligible MCT1 expression2 and therefore slow pyruvate uptake. Consequently, the imaged tumour [1-13C]lactate labelling noted in the HP-13C-MRI cohort is likely to be heavily determined by the increased epithelial cell density within the malignant regions. In the HP-13C-MRI cohort, tumour epithelial cell density was significantly increased compared to the benign areas (*P* = 0.1; **Fig. 2f**) and was positively correlated with [1-13C]lactate SNR (πs = 0.76; *P* = 0.1; **Extended Data Fig. 1**). On standard-of-care 1H-MRI, this high tumour epithelial cell density was reflected by a significantly decreased tumour apparent diffusion coefficient (ADC) on diffusion-weighted imaging (*P* < 0.0001; **Fig. 2g**), with a strong negative correlation observed between the two parameters (πs = -0.67; *P* = 0.0006; **Fig. 2h**). Given the previously reported negative correlation between tumour ADC and [1-13C]lactate SNR2, these data further support the hypothesis that increased epithelial cell density is an important biological driver of increased [1-13C]lactate labelling in PCa compared to the benign prostate. Another potential driver of tissue [1-13C]lactate labelling is the delivery of hyperpolarised [1-13C]pyruvate to the organ of interest. Tumours in the HP-13C-MRI cohort had a significantly higher CD31-derived microvessel density (MVD) on immunohistochemistry (IHC; *P* = 0.034; **Fig. 2i**), which was associated with a significant increase in tumour-derived *K*trans (*P* = 0.003; **Fig. 2j**), a 1H-MRI-derived measure of vascular perfusion and permeability acquired after the intravenous administration of contrast media38. *K*trans and CD31 showed a moderate positive correlation between each other (ρs = 0.55, *P* = 0.02; **Fig. 2k**), which is in agreement with the established link between vascularity and contrast agent delivery to the tissue-of-interest39. While neither of these markers showed a significant relationship with tumour [1-13C]lactate SNR2 (**Extended Data Fig. 2**), increased tumour vascularity may further contribute to the [1-13C]lactate labelling identified in the tumour compared to the benign prostate. Finally, given the importance of MCT1 for cellular [1-13C]pyruvate uptake and subsequent [1-13C]lactate labelling4,22–24, we compared its expression between the benign and tumour samples in the HP-13C-MRI cohort. While MCT1 staining was predominantly epithelial2, it was also heterogeneously distributed, with the resulting epithelial MCT1 density being similar between the benign and tumour areas (*P* = 0.3; **Fig. 2l**). While this agrees with some prior studies investigating prostate MCT1 expression using IHC40,41, the TCGA-PRAD-derived mRNA data42 showed a significant overexpression of *SLC16A1* (MCT1) in tumour samples compared to benign surgical specimens (*P* = 0.009; **Fig. 2m**). The differences between the IHC and TCGA-PRAD results could be explained by limitations in the sample size of the former, and discrepancies between mRNA and protein expression which are well established43. Nonetheless, as no correlation between tumour MCT1 expression and [1-13C]lactate SNR was noted in the HP-13C-MRI cohort2, epithelial MCT1 expression is unlikely to explain the increased tumour [1-13C]lactate labelling observed in this study. ### Increased pyruvate mitochondrial import supports PCa metabolic reprogramming and correlates negatively with clinical [1-13C]lactate labelling Having dissected the potential influence of the four aforementioned biological factors on differential [1-13C]lactate labelling between the benign and malignant prostate, we then explored the dominant intracellular metabolic fate of [1-13C]pyruvate in the two tissue types. As shown in **Fig. 3a**, glycolytic pyruvate formation is followed by either its enzymatic conversion to lactate or mitochondrial flux resulting in the synthesis of acetyl-CoA, which undergoes a condensation reaction with oxaloacetate to form mitochondrial citrate. In the benign prostatic epithelium, further citrate oxidation within the tricarboxylic acid (TCA) cycle is truncated, with the accumulating citrate transported outside the cell through its mitochondrial and plasma membrane transporters44,45. Conversely, early-stage PCa shows fully restored TCA cycle, which has the capacity to generate sufficient NADH for mitochondrial ATP production44,46. Instead of being exported outside the cell, cytosolic tumour citrate is cleaved by ATP citrate lyase to form acetyl-CoA and fuel lipogenesis, to sustain membrane biosynthesis in the proliferating cells47–50. The molecules involved in the rate-limiting reactions of this metabolic reprogramming include mitochondrial pyruvate carriers (MPCs), mitochondrial pyruvate dehydrogenase (PDH), and cytosolic fatty acid synthase (FASN)47–50, and their expression is directly induced by the androgen receptor (AR). Given the similar lactate concentrations in benign and malignant prostate tissue, we hypothesised that the key change in the tumour intracellular pyruvate fate is its increased mitochondrial import, rather than cytosolic reduction to lactate. ![Fig. 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/10/02/2023.09.29.23296228/F3.medium.gif) [Fig. 3:](http://medrxiv.org/content/early/2023/10/02/2023.09.29.23296228/F3) Fig. 3: The role of mitochondrial pyruvate import in PCa metabolic reprogramming and its impact on clinical [1-13C]lactate labelling. **a,** Schematic representation of the proposed differences in the metabolic fate of [1-13C]pyruvate between the benign and malignant prostatic epithelium, with the latter showing increased mitochondrial pyruvate uptake via AR-regulated MPCs to fuel both the restored TCA cycle and FASN-catalysed fatty acid biosynthesis. **b,** Representative fluorescent RNAscope images of epithelial mRNA *LDHA* (gold) and *LDHB* (white) expression, along with IHC images of epithelial MPC2, FASN, and AR expression in the benign and malignant glands. **c,** Scatterplots with bars comparing the expression of epithelial LDH, MPC2, FASN, and AR density between the benign (n = 11) and tumour (n = 11) areas from the HP-13C-MRI cohort. **d,** Top: Plots comparing the log-transformed epithelial and stromal MPC1 and MPC2 densities in the benign (n = 11) and tumour (n = 11) areas from the HP-13C-MRI cohort. Bottom: Spearman’s correlation plots comparing tumour epithelial MPC2 density against tumour epithelial FASN and AR densities, as well as HP-13C-MRI-derived [1-13C]lactate SNR. **e,** Mixed box-and-whisker and scatterplots comparing TCGA-PRAD-derived bulk mRNA expression of *LDHA*, *LDHB*, total LDH, *MPC2*, *FASN*, and *AR* between benign and tumour prostatectomy samples. **f**, Average expression dot plots comparing single-cell RNA-seq epithelial expression of the same genes from publicly available EGAS00001005787 (left) and GSE176031 (right) datasets. **g**, H&E map of a surgical specimen including areas of benign and malignant prostate with corresponding spatial transcriptomics maps demonstrating the expression of total LDH, *MPC2*, and *FASN* obtained from a publicly available EGAS00001006124 dataset. In panels **c-**e, lines are median values and bars are interquartile ranges, with *P* derived using the Mann-Whitney *U* test. Scale bars denote 5-10 μm To test this hypothesis, a random forest tissue classifier2,51 was applied to the same segmented benign and tumour ROIs on histopathology used for HP-13C-MRI data extraction (**Fig. 2a**): total epithelial LDH mRNA density was derived, which we have previously shown to correlate strongly with tumour [1-13C]lactate labelling in the primary cohort2. We found that this total epithelial LDH density was significantly reduced in tumour ROIs compared to the benign areas (*P* = 0.024; **Fig. 3b,c**), which is in keeping with findings derived from open-source TCGA-PRAD bulk mRNA sequencing data42 (**Fig. 3e**), the epithelial single-cell mRNA sequencing data obtained from two recent publications by Tuong *et al.*52 and Song *et al.*53 (**Fig. 3f**), and spatial transcriptomics data from a study by Erickson *et al*.54 (**Fig. 3g**). While this did not have an apparent effect on the intracellular lactate pool, a reduction in the total LDH expression in tumour epithelium may be functionally consistent with its increased capacity for mitochondrial [1-13C]pyruvate import. To evaluate the role of mitochondrial import, we used IHC to measure the expression of MPC1 and MPC2, with the latter previously shown to be significantly overexpressed in the human prostate under the direct transcriptional control of nuclear AR49. Here, we show almost exclusive epithelial MPC1 and MPC2 localisation in both benign and malignant prostate (*P* < 0.05 for both; **Fig. 3d**), with tumour epithelial MPC2 density being significantly higher compared to the healthy tissue (*P* = 0.005; **Fig. 3b,c**), an observation also supported by bulk42, single-cell52,53, and spatial54 mRNA sequencing data presented in **Fig. 3e-g**. Importantly, tumour epithelial MPC2 density showed a strong negative correlation with clinical [1-13C]lactate labelling (ρs = -0.80, *P* = 0.005; **Fig. 3d**), which aligns with a prior preclinical report showing an increase in glycolytic intermediates following MPC inhibition55. Taken together with our previous work which showed strong positive associations between tumour [1-13C]lactate SNR, total epithelial LDH density, and epithelial *LDHA*/*PDHA1* ratio2, these results indicate that tumour [1-13C]lactate labelling is likely to be a function of both the epithelial capacity for cytosolic LDH-catalysed pyruvate reduction, and MPC-driven mitochondrial pyruvate import. We then explored the potential role of MPC2 in supporting AR-driven tumour reprogramming towards a lipogenic phenotype, by measuring the expression of both AR and FASN in our samples. As expected, tumour areas exhibited a significant increase in epithelial FASN density (*P* = 0.003; **Fig. 3b** and **3c**; supported by **Fig. 3e-g**), which correlated strongly with the epithelial MPC2 density (ρs = 0.93, *P* = 0.001; **Fig. 3d**). Both MPC2 and FASN also showed strong positive correlations with epithelial nuclear AR density (ρs = 0.89 and 0.68, respectively; *P* < 0.05 for both; **Fig. 3d**), corroborating prior preclinical reports49,50 identifying MPC2 and FASN as AR transcriptional targets. Overall, these findings suggest the potential role of mitochondrial pyruvate flux for enabling PCa metabolic rewiring and highlight its importance for clinical HP-13C-MRI interpretation. ### Spatial metabolomics combined with a deep learning metabolic classifier suggest differences in pyruvate metabolism between glycolytic benign prostate and lipogenic cancer To further assess the role of mitochondrial pyruvate import in supporting lipogenic PCa metabolic reprogramming, we interrogated the DESI-MSI data derived from benign (n = 695) and tumour (n = 468) epithelial ROIs drawn manually on haematoxylin-and-eosin (H&E) slides obtained from the DESI-MSI cohort (**Fig. 2e** and **Fig. 4e**). We used the ROIs to extract DESI-MSI-derived metabolites from three relevant pathways within the KEGG database: glycolysis/gluconeogenesis, TCA cycle, and fatty acid biosynthesis pathways (**Supplementary Table 1**). Subsequent metabolic pathway enrichment analysis (MPEA) showed significant enrichment of the TCA cycle pathway in the benign ROIs (*P* = 0.01; **Fig. 4a**), which at the level of individual metabolites was primarily driven by citrate and phosphoenolpyruvate (PEP) (*P* < 0.001 for both; **Fig. 4b,c**). Significant citrate enrichment in the benign prostate is well established44,46 and therefore serves as a validation of the DESI-MSI approach. In turn, PEP enrichment in the benign tissue may be explained by both its glycolytic production and anaplerotic synthesis from cytosolic oxaloacetate, a process aimed at regenerating the TCA cycle carbons56. In addition, benign prostatic epithelium showed a non-significant enrichment of the glycolysis/gluconeogenesis pathway relative to PCa (*P* = 0.27; **Fig. 4a**), with PEP and glucose (*P* < 0.0001; **Fig. 4b,c**) being the two metabolite hits. As in the case of citrate, the observed significant enrichment of glucose in the benign epithelium compared to intermediate-risk organ-confined PCa aligns with the published literature and explains the well-known limited diagnostic utility of [18F]-FDG-PET in PCa57. Importantly, of the three pathways analysed, only metabolites involved in fatty acid biosynthesis were significantly enriched in tumour epithelial ROIs (*P* < 0.0001; **Fig. 4a**), including tetradecanoic, oleic, and palmitoleic acids being the three metabolite hits (*P* < 0.0001 for all; **Fig. 4b,c**). Notably, another significantly enriched tumour epithelial metabolite was S-Acetyldihydrolipoamide-E, which is a pyruvate-derived precursor of mitochondrial acetyl-CoA formed as part of the PDH-catalysed reaction58. Increased tumour epithelial S-Acetyldihydrolipoamide-E enrichment supports the notion of increased mitochondrial oxidation as the primary metabolic fate of epithelial pyruvate in human PCa, which has not been reported previously. ![Fig. 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/10/02/2023.09.29.23296228/F4.medium.gif) [Fig. 4:](http://medrxiv.org/content/early/2023/10/02/2023.09.29.23296228/F4) Fig. 4: Spatially resolved metabolic profiling of the benign and malignant prostatic epithelium. **a**, Outputs of the DESI-MSI-derived MPEA demonstrating comparative enrichment of KEGG glycolysis, fatty acid biosynthesis, and TCA cycle pathways in benign (n = 695) and tumour (n = 468) ROIs from the spatial metabolomics cohort. **b**, Volcano plot showing individual differentially enriched metabolites from the three KEGG pathways between the benign and tumour ROIs; *P* were derived using FDR-corrected Wilcoxon rank sum test. **c,** Mixed box-and-whisker and scatterplots comparing epithelial abundance of the key differentially enriched metabolites between benign and tumour epithelial ROIs; the data are presented as median and interquartile range; *P* were derived using the Mann-Whitney *U* test. For illustrative purposes, outliers were removed using the ROUT method60 with Q = 5%. **d**, tSNE plot of DESI-MSI data acquired from the benign and tumour ROIs focusing on ions corresponding to metabolites related to the three KEGG pathways. **e**, Summary diagram describing the key steps in developing a deep-learning-based metabolic tissue classifier using DESI-MSI-derived metabolites from the three KEGG pathways to discriminate between benign and tumour epithelial ROIs; the SHAP plot lists the ten most important metabolites used by the final model to achieve a median performance of 0.91 presented in an AUC plot in **f**. An alternative approach to testing the same hypothesis in a more clinically applicable way was to assess a list of differentially enriched metabolites by building a DESI-MSI-based metabolic classifier of benign and malignant epithelial regions. Using metabolite data from the same three KEGG pathways, we first noted a clear unsupervised separation between the benign and tumour epithelial cell clusters shown in **Fig. 4d**. We then trained and tested a deep-learning-based metabolic classifier (**Fig. 4e**; **link to code availability**) that achieved a median performance of 0.91 for differentiating between benign and tumour epithelial ROIs following five-fold cross-validation (**Fig. 4f**). Using Shapley additive explanations59, we identified the list of the ten most important metabolites for the classification, all presented in **Fig. 4e**. In keeping with the MPEA results, six of these were fatty acids, and these were more abundant in correctly classified tumour epithelial ROIs. The remaining four metabolites were glucose, malate, citrate, and malonate, and in comparison, were in high abundance within benign ROIs. Unsurprisingly, lactate had little impact on the developed classifier, ranking 28/31 in the list of the most important metabolites, with a mean absolute SHAP value of 0.009 (for comparison, the same value for oleic acid was 0.056). Overall, these results complement both the MPEA findings and tissue biomarker expression from the HP-13C-MRI cohort, suggesting that increased tumour mitochondrial pyruvate flux is an important metabolic feature of human PCa that can be probed by both *in vivo* and *ex vivo* using metabolic imaging techniques. However, the impact of epithelial mitochondrial pyruvate import on hyperpolarised [1-13C]lactate labelling can only be assessed in tumour but not benign tissue, due to the lack of discernible HP-13C-MRI signal in the less cellular and vascular benign prostate. ### Cribriform PCa may show no [1-13C]lactate signal due to high MPC density and low LDH expression Having explored the mechanisms behind the ability of clinical HP-13C-MRI to detect PCa by distinguishing it from the healthy tissue, we then attempted to assess the differential [1-13C]lactate labelling between tumours harbouring cribriform and non-cribriform GP4 disease. As shown in **Table 1** and **Fig. 5a-c**, 2/3 ISUP grade group (GG) 3 ICC lesions from the HP-13C-MRI cohort had no detectable [1-13C]lactate signal despite sufficient [1-13C]pyruvate delivery and restricted diffusion on 1H-MRI ADC indicating high cellularity. Notably, these lesions were found in patients who also harboured two contralateral small-volume ISUP GG2 non-cribriform tumours, which were both detectable on HP-13C-MRI (**Fig. 5a,b**), therefore suggesting that the absence of [1-13C]lactate signal in cribriform disease was not due to technical issues. ![Fig. 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/10/02/2023.09.29.23296228/F5.medium.gif) [Fig. 5:](http://medrxiv.org/content/early/2023/10/02/2023.09.29.23296228/F5) Fig. 5: Comparative assessment of biological factors underpinning differential HP-13C-MRI-visibility of cribriform and non-cribriform PCa. **a-b**, Whole-mount H&E, ADC, along with [1-13C]pyruvate and [1-13C]lactate SNR maps demonstrating the presence of large, cellular ISUP GG3 lesions with dominant ICC component that were HP-13C-MRI occult compared to contralateral small-volume foci of ISUP GG2 disease with <5% non-cribriform GP4 glands. **c**, Comparator case of a HP-13C-MRI-visible large cribriform ISUP GG3 tumour that was also visible on 1H-MRI ADC. **d**, Representative H&E slides, along with corresponding IHC-derived CD31, MCT1, MPC2, FASN, AR, and RNAscope-derived total LDH images obtained from HP-13C-MRI-visible and HP-13C-MRI-occult lesions shown in panels **a-c**. **e**, Mixed box-and-whisker and scatterplots comparing [1-13C]pyruvate and [1-13C]lactate SNR, epithelial cell density, CD31 density, as well as epithelial MCT1, LDH, MPC2, FASN, and AR density between HP-13C-MRI-visible (n = 11) and HP-13C-MRI-occult (n = 2) lesions from the hyperpolarised MRI cohort. In **e**, the data are presented as median with interquartile range; *P* were derived using the Mann-Whitney *U* test. Scale bars in **d** denote 5-50 μm. To explain these findings within the same biological validation framework applied to the prior research questions above, we compared the metabolic and non-metabolic drivers of [1-13C]lactate labelling between the HP-13C-MRI-visible and HP-13C-MRI-occult tumours. Interestingly, while epithelial cell density, CD31-derived MVD, and epithelial MCT1 expression were similar between the two tumour types (*P* > 0.05 for all; **Fig. 5d,e**), ICC lesions with no discernible [1-13C]lactate labelling showed almost negligible epithelial LDH density together with a significant increase in epithelial MPC2, FASN, and AR expression (*P* = 0.026 for all; **Fig. 5d,e**). These findings provide additional support for the hypothesis that tumour [1-13C]lactate labelling reflects the differential intracellular metabolic fate of pyruvate between LDH-catalysed reduction to lactate, and MPC-driven mitochondrial import to facilitate oxidative metabolism. Moreover, this observation highlights the multifactorial nature of signal generation with HP-13C-MRI, including both metabolic and physical or physiological tissue characteristics. Specifically, if non-metabolic factors such as epithelial cell density and vascularity can be corrected for, then tumour [1-13C]lactate labelling patterns may reflect the true metabolic differences between lesions of different histological phenotypes ### Gleason pattern 4 PCa consists of lipogenic cribriform and glycolytic non-cribriform metabolic phenotypes Given the absence of HP-13C-MRI signal in most cribriform PCa, we conducted a granular comparative metabolic analysis of different Gleason pattern glands included in this study. We previously demonstrated a significant positive correlation between tumour [1-13C]lactate labelling and the percentage of GP4 (%GP4) disease in the HP-13C-MRI cohort2. However, since all but one HP-13C-MRI-visible lesion consisted of non-cribriform GP4 glands, this reported relationship actually existed between [1-13C]lactate SNR and the percentage of non-cribriform GP4 disease. When all 15 lesions were divided using the %GP4 cut-off of 10%, no difference in [1-13C]lactate labelling was noted between high %GP4 and low %GP4 tumours (*P* = 0.46; **Fig. 6a**). Conversely, when the high %GP4 group was further subdivided by the dominant GP4 subtype, ICC or non-cribriform, the three resulting tumour types all showed significantly different [1-13C]lactate labelling (*P* < 0.05 for all; **Fig. 6a**) in the absence of notable changes in tumour ADC, *K*trans, and epithelial MCT1 density (*P* > 0.05 for all; **Fig. 6a**). ![Fig. 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/10/02/2023.09.29.23296228/F6.medium.gif) [Fig. 6:](http://medrxiv.org/content/early/2023/10/02/2023.09.29.23296228/F6) Fig. 6: Comparative metabolic characterisation of intermediate-risk PCa with varying percentage and phenotype of GP4 disease. **a**, Scatterplots with bars comparing clinical [1-13C]lactate labelling between tumours divided based on the percentage (far left; n = 7 for low %GP4 and n = 8 for high %GP4 lesions, respectively) and histological subtype (second from left; n = 7 for low %GP4, n = 5 for non-cribriform high %GP4, and n = 3 for cribriform high %GP4 lesions, respectively) of GP4 disease. Intergroup comparisons of 1H-MRI-derived tumour ADC and *K*trans, as well as tissue-based epithelial MCT1 density are also presented. **b**, Mixed box-and-whisker and scatterplots comparing tissue-based total epithelial LDH, nuclear HIF-1α, MPC2, FASN, and nuclear AR between ROIs harbouring individual GP3 (n = 50), non-cribriform GP4 (n = 22), and GP4 ICC (n = 17) glands. **c**, Representative H&E, RNAscope, and IHC images illustrating differential expression of tissue-based metabolic biomarkers between the three Gleason pattern glands. In **a-b**, the data are presented as median (denoted by the bars or boxes) with interquartile range (denoted by the bars or whiskers). *P* were derived using the Mann-Whitney *U* test. We subsequently measured the expression of key metabolic biomarkers at the level of individual GP3 (n = 50), GP4 non-cribriform (n = 22), and GP4 ICC (n = 17) glands. Total epithelial LDH expression was highest in non-cribriform GP4 glands compared to both GP3 and GP4 ICC (*P* < 0.0001 for both; **Fig. 6b,c**), with no difference noted between the latter two morphological subtypes (*P* = 0.4; **Fig. 6b,c**), mirroring the trends noted in **Fig. 6a** for [1-13C]lactate labelling across histological types. At the transcriptional level, this observation was also in line with a similar pattern of nuclear HIF-1α expression, which was also significantly increased in non-cribriform GP4 glands compared to the other glands (*P* < 0.05 for both; **Fig. 6b,c**). Conversely, both MPC2 and FASN increased significantly in a stepwise fashion across the three Gleason phenotypes, with a similar trend noted for the epithelial nuclear AR expression (*P* < 0.05 for all; **Fig. 6b** and **c**). Overall, these results confirm that clinical [1-13C]lactate labelling reflects the tumour-specific pyruvate metabolic fate, which appears to favour a reduction to lactate in non-cribriform lesions and towards mitochondrial flux in cribriform GP4 glands. These tissue-based metabolic metrics were then compared with endogenous metabolite measurements by divided lesions from the DESI-MSI cohort into similar groups depending on the percentage and morphology of GP4 disease. In line with the [1-13C]lactate labelling pattern shown in **Fig. 6a**, DESI-MSI-derived epithelial lactate was similar between cores derived from low (n = 19) and high (n = 19) %GP4 lesions (*P* = 0.93; **Fig. 7a**). Moreover, when the cribriform and non-cribriform lesions were separated, the epithelial lactate pool was also found to be significantly increased in non-cribriform high %GP4 lesions compared to both high %GP4 ICC tumours and low %GP4 disease (*P* < 0.05 for all; **Fig. 7a** and **b**). In contrast to the HP-13C-MRI cohort, high %GP4 ICC lesions had the lowest epithelial lactate abundance and showed decreased epithelial cellularity, which was significantly reduced compared to high %GP4 non-cribriform tumours and was similar to low %GP4 PCa (*P* = 0.04 and 0.26, respectively; **Fig. 7a**). Compared to cribriform lesions from the HP-13C-MRI cohort, which required a sufficient epithelial cell density to be 1H-MRI-visible and therefore eligible for study inclusion2, this is more representative of routine clinical practice where only 17% of pure cribriform lesions are detected on conventional imaging due to their loose cellular composition18,61. Hence, considering the importance of epithelial cell density in generating clinical HP-13C-MRI signal, as well as the low epithelial LDH density in 2/3 cribriform lesions in the primary cohort, the low tissue lactate pool in ICC lesion highlights the potential limitation of using HP-13C-MRI for detecting cribriform PCa and distinguishing it from less aggressive tumour subtypes. ![Fig. 7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/10/02/2023.09.29.23296228/F7.medium.gif) [Fig. 7:](http://medrxiv.org/content/early/2023/10/02/2023.09.29.23296228/F7) Fig. 7: Spatially resolved metabolic phenotyping of intermediate-risk human PCa. **a**, Mixed box-and-whisker plots comparing DESI-MSI derived epithelial lactate abundance between cores derived from lesions with low %GP4 (n = 19) and high %GP4 (n = 19) (far left), as well as between cores sub-stratified by the division of high %GP4 lesions into those with dominant non-cribriform (n = 10) and cribriform (n = 9) GP4 component (second left). Plots demonstrating intergroup comparison of epithelial palmitoleic acid pool and epithelial cell fraction are also presented; the data are presented as median (boxes) with interquartile ranges (whiskers); *P* were derived using the Mann-Whitney *U* test. **b**, Representative fused H&E and DESI-MSI lactate and palmitoleic acid images demonstrating differential metabolite abundance between GP3 (yellow), GP4 non-cribriform (red), and GP4 cribriform (black) glands. **c**, MPEA outputs demonstrating differential enrichment of glycolysis, TCA cycle, and fatty acid biosynthesis KEGG pathways and individual metabolites between ROIs classified as benign and GP3 (n = 360), non-cribriform GP4 (n = 70), and cribriform GP4 (n = 36) glands; *P* were derived using FDR-corrected Wilcoxon rank sum test. **d** Diagram summarising the key results of this study and highlighting the presence of lactate-rich non-cribriform and lactate-poor, fatty acid-rich cribriform GP4 phenotypes within intermediate-risk PCa. Conversely, in line with the FASN expression trend, epithelial palmitoleic acid showed a similar stepwise increase across the three tumour phenotypes that was independent of epithelial cellularity (*P* < 0.05 for all; **Fig. 7a** and **b**). This observation aligned with the MPEA results (**Fig. 7c**), which compared the enrichment of KEGG glycolysis, TCA cycle, and fatty acid biosynthesis pathways across individual benign/GP3 glands (n = 360), non-cribriform GP4 glands (n = 70), and GP4 ICC glands (n = 36). Specifically, fatty acid biosynthesis was consistently enriched in more aggressive disease, and highest in GP4 ICC glands (**Fig. 7c**), in line with a recent study involving transcriptomic profiling of cribriform lesions16. Conversely, epithelial lactate as part of the glycolytic pathway was significantly enriched in non-cribriform GP4 tumours compared to both benign and GP3 glands, as well as GP4 ICC lesions (*P* < 0.05 for both; **Fig. 7c**). As summarised in **Fig. 7d**, these results indicate that GP4 PCa is comprised of two distinct metabolic phenotypes: lactate-rich, less aggressive non-cribriform tumours and lipid-rich, more aggressive ICC lesions. ## Discussion HP-13C-MRI is an emerging clinical tool to non-invasively image tissue metabolism. This exploratory clinical study deployed an array of *in vivo* and *ex vivo* metabolic imaging techniques to explore the biological mechanisms underpinning the ability of HP-13C-MRI to study human PCa. This has demonstrated how metabolic signatures within PCa can vary between histological subtypes and has revealed two metabolically distinct phenotypes within intermediate-risk GP4 disease. Our results highlight several important considerations for future clinical translation of HP-13C-MRI both in PCa and in other tumour types and could have important implications for characterising tumours in the future. First, we show that clinical [1-13C]lactate labelling offers a measure of absolute tissue lactate abundance, which can be influenced by both metabolic and non-metabolic (e.g. cell density and blood flow) tissue properties. For example, the increase in the absolute tissue lactate pool in HP-13C-MRI-visible PCa compared to either occult PCa or benign prostate can be partly explained by the increased number and density of tumour epithelial cells. However, further factors play a role as two HP-13C-MRI-cribriform lesions were occult on imaging suggesting that epithelial cell density and vascularity are insufficient for generating clinical [1-13C]lactate labelling in tumours in the context of low epithelial LDH expression, reduced endogenous lactate pool, and increased capacity for mitochondrial pyruvate uptake. These results highlight the importance of using a similar biological validation framework in future clinical HP-13C-MRI studies to understand the exact mechanisms underlying the observed changes in [1-13C]lactate labelling in different cancers and various clinical scenarios. For example, since MPC2 is directly regulated by AR, early metabolic response to androgen deprivation therapy (ADT) could involve an increased tumour [1-13C]lactate labelling due to reduced [1-13C]pyruvate mitochondrial uptake; we have previously shown a similar early increase following treatment has been shown in human breast cancer9,62. Conversely, a later reduction in tumour [1-13C]lactate signal demonstrated in a case report by Aggarwal *et al*.6 may reflect a combination of impaired tumour perfusion due to ADT-induced endothelial damage63 and decreased tumour epithelial cell density64 due to the loss of glandular architecture and apoptotic cell death65. Similarly, if HP-13C-MRI is used for imaging response to poly (ADP-ribose) polymerases (PARP) inhibitors, early post-treatment scans may show increased [1-13C]lactate labelling due to replenished NAD+ levels, 90% of which are consumed by DNA-damage activated PARPs9,66. The lack of [1-13C]lactate signal in aggressive cribriform lesions, which are known to harbour increased genomic instability67,may also be explained by a decreased NAD+ pool in addition to low LDH expression, and an increase in MPC-driven mitochondrial [1-13C]pyruvate uptake. While non-invasive assessment of the percentage of lactate-rich non-cribriform GP4 glands may be helpful in some clinical settings such as active surveillance (AS), accurate detection and quantification of cribriform disease is of high diagnostic importance for several reasons. First, it has been shown that patients with non-cribriform ISUP GG2 disease have the same prognosis as men with ISUP GG1 lesions, whereas those with an ICC component had more frequent surgical and radiotherapy failure68,69. The European Association of Urology PCa guidelines recommend against offering AS to patients with even small amounts of ICC detected on diagnostic biopsies70,71. The absence of [1-13C]lactate signal in high-probability lesions as shown in **Fig. 5a** and **5b** might raise the suspicion for the underlying cribriform morphology. However, as most cribriform lesions are difficult to visualise on standard 1H-MRI18,72, the potential reduction in their [1-13C]lactate labelling may limit the utility of HP-13C-MRI for increasing their detection rate. Furthermore, in cases of mixed lesions with other Gleason patterns present, as well as intraductal carcinoma (IDC) that is often composed of a cribriform component, lactate-poor ICC glands may artificially decrease tumour [1-13C]lactate labelling with potential implications for reducing the accuracy of lesion characterisation. Conversely, the results here suggest that future development of novel 13C-labelled probes targeting fatty acid metabolism may yield more consistent results given the continuous upregulation of lipid metabolism noted throughout PCa development73,74. Future studies could further assess the role of PET agents such as 11C-acetate and 11C-choline for detecting ICC and IDC, as well the applicability of prostate-specific membrane antigen (PSMA) PET for detecting cribriform PCa where there are currently conflicting results75–77. This study adds a new dimension to the rapidly growing body of literature investigating the molecular heterogeneity of different GP4 phenotypes by demonstrating the clinical potential for metabolic characterisation of intermediate-risk disease. Recent studies including bulk16,78 and single-cell15 mRNA sequencing have not focused on metabolic gene sets, and clinical studies involving spatial metabolomics techniques79–83 have not distinguished between GP4 subtypes. Our future work will build on the current results by identifying transcriptomic drivers behind the described metabolic reprogramming in different GP4 phenotypes, as well as GP3 and GP5 glands. If correlated with clinical outcomes and imaging findings, this work may yield important insights into metabolic characteristics of clinically significant PCa, thereby navigating the development of more specific metabolic imaging probes beyond [1-13C]pyruvate. The limitations of this study include the lack of a direct correlation between HP-13C-MRI and spatial metabolomics data in the same cohort, which was due to the challenges in using formalin-fixed and paraffin-embedded tissue for DESI-MSI. The current study has compensated for this by careful histopathological matching between the two prospective patient cohorts, with the key results being consistent both with the two patient groups and with previously published data. In addition, this study did not include patients with GP5 disease or IDC, which will be comprehensively profiled in the future. Importantly, future studies will include larger patient cohorts, helping to assess the generalisability of our exploratory findings, as well as to investigate the role of other metabolic pathways in differentiating different PCa subtypes. In conclusion, this study demonstrates the impact of metabolic and non-metabolic tissue characteristics on the characterisation of PCa using HP-13C-MRI. Using a combination of macroscopic and microscopic imaging techniques, we show that clinical [1-13C]lactate labelling is a function of epithelial cell density, vascularity, cellular capacity for pyruvate reduction to lactate, and MPC-driven mitochondrial pyruvate import. We have demonstrated that these characteristics differ significantly between both benign and malignant prostate tissue, and between intermediate-risk tumours containing cribriform and non-cribriform GP4 glands. If validated in larger cohorts and expanded to other tumour subtypes, these results have important implications for the development and clinical translation of PCa metabolic imaging techniques. ## Methods ### Patient selection and ethics This clinical study analysed prospectively collected clinical, imaging, and tissue-based data obtained from patients from two separate surgical cohorts, as shown in **Fig.1**. The HP-13C-MRI cohort was a subset of a previously reported2 population of consecutive male patients with at least one MR-visible (> 1cm, Likert/PI-RADS 4-5) histologically proven PCa focus (ISUP GG ≥2) scheduled for RARP in our centre. All patients provided written consent to participate in the original prospective study (Molecular Imaging and Spectroscopy with Stable Isotopes in Oncology and Neurology – Imaging metabolism in prostate [MISSION-Prostate] protocol), retrospective analysis of which was approved by the institutional review board (National Research Ethics Service Committee East of England, Cambridge South, Research Ethics Committee number 16/EE/0205). Of the ten patients recruited originally, eight were included in the current analysis due to technical failure of HP-13C-MRI in the remaining two patients2. Following a review of post-surgical whole-mount histological specimens, 15 lesions were reported, with their key morphological characteristics reported in **Table 1**. To validate the key findings obtained in the HP-13C-MRI cohort by directly measuring tissue metabolite concentrations in fresh-frozen prostatectomy samples, we then retrospectively identified a secondary cohort of treatment-naïve surgical patients who had been originally recruited into an ethically-approved prospective national study (DIAMOND, National Research Ethics Service Committee East of England, Cambridge South, Research Ethics Committee number 03/018) following the provision of written consent. To ensure comparability with the HP-13C-MRI cohort, a total of 13 patients with 15 lesions matching the key histopathological characteristics of tumours from the primary cohort (**Table 1**) were selected for analysis as part of this study. Further description of tumour histopathological assessment and matching is provided in the corresponding section below. ### HP-13C-MRI acquisition and analysis In the primary cohort, pre-surgical HP-13C-MRI was acquired on a clinical 3.0 T MR system (MR750, GE Healthcare, Waukesha WI, USA) using a bespoke 1H/13C endorectal receive coil84 following the protocol detailed in the original MISSION-Prostate report2. In brief, samples containing 1.47 g of [1-13C]pyruvic acid (Sigma Aldrich, St Louis MO, USA) and 15 mM electron paramagnetic agent (EPA) were hyperpolarised using a clinical hyperpolariser (SPINlab; 5T Research Circle Technology, GE Healthcare, Waukesha WI, USA) by microwave irradiation at 139 GHz at ∼0.8 K for ∼3 h followed by rapid dissolution in 38 mL of superheated sterile water and filtration to remove EPA to a concentration below ≤3 μM2,8. Radiofrequency pulses with a nominal flip angle of 15° were applied using a clamshell RF coil (GE Healthcare, Waukesha WI, USA) to acquire a 20×20 cm2 FOV with a matrix size of 32×32 and a temporal resolution of 4 s for 20 time points. Images had a true in-plane resolution of 12.5×12.5 mm2 and were reconstructed with a resolution of 128×128. The imaging data were reconstructed in MATLAB (MathWorks, Natick MA, USA)2,8, with [1-13C]pyruvate and [1-13C]lactate SNR derived from ROIs containing areas of histopathologically-proven PCa (n = 15) and contralateral benign (n = 15) tissues. The ROIs were originally drawn on H&E whole-mount pathology (WMP) maps by an experienced genitourinary pathologist (A.Y.W.) (**Extended Data Fig. 3**). Corresponding lesions were then outlined on HP-13C-MRI metabolite maps (**Extended Data Fig. 3**) in consensus by two fellowship-trained consultant radiologists (T.B., F.A.G.) and a research fellow (N.S.) with 13-, 14-, and 5-years’ experience in prostate MRI, respectively85. ### 1H-MRI acquisition and analysis Either before or immediately after the HP-13C-MRI acquisition, all patients in the primary cohort underwent standard-of-care multiparametric MRI of the prostate. The protocol, detailed in the prior cohort description2, included axial T1-weighted fast spin echo (FSE), high-resolution T2-weighted 2D fast recovery fast spin echo (FRFSE, echo time [TE] 98-107 ms, field-of-view [FOV] 22 × 22 cm2, acquisition matrix 320– 384 × 256, slice thickness 3 mm with 0 mm gap, 3 signal averages, repetition time [TR] 3000–5000 ms, echo train length 16, receiver bandwidth ± 31.25 or ± 41.67 kHz), diffusion-weighted imaging (DWI, spin-echo echo-planar imaging pulse sequence with b values of 150, 550, 750, 1000, and 1400, with a separate high b value acquisition of 2000 s/mm2), and dynamic contrast-enhanced MRI (DCE-MRI, axial 3D fast spoiled gradient echo [FSPGR], TR/TE 4.1/1.8 ms, FOV 24 × 24 cm2, following bolus injection of gadobutrol [Gadovist, Bayer Healthcare, Berlin, Germany] via a power injector, rate 3 mL/s [dose 0.1 mmol/kg], temporal resolution 7 s). ADC maps were calculated automatically, representing the quantitative assessment of tissue density. The mean ADC values corresponding to the areas of benign (n = 15) and malignant (n = 15) prostate were extracted from the same ROIs used for the previously described HP-13C-MRI analysis (**Extended Data Fig. 3**). The same ROIs were also used for extracting *K*trans values as measures of vascular perfusion and permeability (**Extended Data Fig. 3**), with DCE-MRI analysis performed in an in-house-developed MATLAB software that was used to generate B1 maps2,8. These maps were then transferred to MIStar (Apollo Medical Imaging, Melbourne, Australia) to generate B1-corrected T1 maps, to perform motion correction of the DCE-MRI data using a 3D affine model, and for pharmacokinetic modelling using the standard Tofts model2,8. In patient 6 (**Table 1**), DCE-MRI was not performed for technical reasons. ### Histopathological assessment and digital pathology analysis In the HP-13C-MRI cohort, prospectively obtained whole-mount surgical sections were formalin fixed and paraffin-embedded following the routine clinical procedure. The resulting FFPE slides were stained with H&E and digitised for subsequent histopathological assessment at 40x using Aperio CS2 Digital Pathology Scanner (Aperio Tech., Oxford, UK)86. The latter was performed by a fellowship-trained genitourinary pathologist who was a member of the 2019 ISUP Consensus Working Group on PCa grading14 and has more than 20 years’ experience of clinical PCa histopathological assessment. First, whole-mount slides were assessed for the presence of PCa foci and assigning their ISUP GG, quantifying %GP4, and identifying the dominant GP4 subtype in ISUP GG 2-3 lesions. As described in the previous sections, tumour (n = 15) and contralateral benign tissue ROIs (n = 15) were used for HP-13C-MRI and 1H-MRI segmentation (**Extended Data Fig. 3**). The same ROIs were subsequently transposed onto IHC and RNAscope whole-mount slides for analyses described in the corresponding sections below. In addition, as described in the prior publication2, whole-mount H&E maps were used to draw standardised random ROIs encompassing tumour foci represented by clear GP3 (n = 50), non-cribriform GP4 (n = 22; poorly formed glands, fused glands, glomeruloid glands), and cribriform GP4 (n = 17; invasive cribriform carcinoma) glands (**Extended Data Fig. 4**); the identification of the latter was performed according to the 2021 ISUP consensus definition of cribriform pattern PCa87. These ROIs were also transposed onto IHC and RNAscope slides for the corresponding analyses. The key histopathological characteristics of tumours from the HP-13C-MRI cohort (ISUP GG, %GP4, dominant GP4 subtype) were used to guide the search for a matching set of 15 lesions to form the secondary DESI-MSI cohort. The original fresh-frozen surgical sample collection and preparation was performed by the same expert genitourinary pathologist (A.Y.W.) according to the previously described protocol88. In brief, RARP specimens were transferred on ice to the laboratory within 30 min of surgical resection, with multiple punch biopsies removed using a standard 4-6 mm skin punch. The sites of the punched cores were marked on a “map” diagram (**Fig. 1**), from which cores including benign (n = 61) and tumour tissues (n = 56) were selected for this study. H&E staining was performed following the FF sample preparation described in the DESI-MSI section below, with the resulting slides scanned at 40x using the same Aperio Digital Pathology Scanner. As in the HP-13C-MRI cohort, tissue assessment and annotation were performed in QuPath 0.2389 by the same genitourinary pathologist (A.Y.W.). First, depending on the presence and amount of tumour tissue, each core was labelled as either benign (no tumour tissue, n = 61), mixed (<50% of tumour tissue, n = 23), and tumour (≥50% of tumour tissue, n = 38) (**Extended Data Fig. 5**) to ensure accurate comparison between benign and malignant prostatic tissues during the whole-core lactate abundance assessment detailed below. Second, the epithelial component of each individual core was annotated manually as shown in **Fig. 2e** and labelled according to its morphological characteristics (benign, GP3, GP4 non-cribriform, and GP4 ICC). This enabled us to assign the overall ISUP GG and %GP4, as well as identify the dominant GP4 subtype both at the level of an individual core and for the combination of cores obtained from a single large tumour focus. Both whole-core (n = 69) and individual epithelial ROIs (n = 695 benign and n = 468 tumour) were used to extract the DESI-MSI data as described in the relevant section below. In both cohorts, spatially resolved tissue analysis at the level of epithelial and stromal compartments was critical due to the known metabolic compartmentalisation between the two tissue types in the human prostate. In the HP-13C-MRI cohort, we used a previously described random forest tissue classifier2 embedded into the HALO v3.2.1851.266 (Indica Labs, Albuquerque, NM, USA) software. In short, the classifier was trained using manual annotations of epithelial and stromal areas (**Extended Data Fig. 6**) according to the previously described methodological pipeline51. In the DESI-MSI cohort, we used the HALO AI platform to train a similar deep neural network (DenseNet V2, path size: 100 × 100 pixels) tissue classifier to again distinguish between epithelial and stromal cellular compartments. Since the DESI-MSI data was extracted from manually drawn epithelial ROIs, here the classification was only required for quantifying the relative distribution of epithelial and stromal compartments within a given core (quantified as the area of epithelial/stromal cells per the total tissue area, both in mm2). To train the classifier, we used a subset of pathologist-defined annotations (**Extended Data Fig. 6**) of epithelial, stromal, and background regions from n = 133 separate tissue regions with a total area of 171.4 mm2. Classification outputs were visualised as mark-up images (**Extended Data Fig. 6**) and sent back to the expert pathologist (A.Y.W.) for corrections as part of an active learning process. The classifier performance was cross-validated against ground truth annotations using the unseen dataset of the validation cohort. ### Immunohistochemistry To assess biological factors such as epithelial capacity for [1-13C]pyruvate uptake, tissue vascularity, and inferred [1-13C]pyruvate metabolic fate, we used IHC to stain the HP-13C-MRI cohort samples for the following targets: MCT1 (membrane pyruvate importer; Cat. No. HPA003324, Atlas Antibodies, Bromma, Sweden), MPC1 and MPC2 (mitochondrial pyruvate carriers; Cat. No. PAB28306, Abnova, Taipei, Taiwan for MPC1; Cat. No. D417G, Cell Signaling Technology, Danvers MA, USA for MPC2), FASN (cytosolic fatty acid synthase; Cat. No. 3180, Cell Signaling Technology, Danvers MA, USA), AR (nuclear receptor regulating the expression of MPCs and FASN; Cat. No. NCL-AR-318, Novocastra, Newcastle, UK), HIF-1α (nuclear transcription factor regulating the expression of *LDHA*; Cat. No. ab51608, Abcam, Cambridge, UK), and CD31 (endothelial biomarker; Cat. No. M0823, Dako, Santa Clara CA, USA). All antibodies have been previously validated in our centre using positive and negative tissue controls under the supervision of specialist pathologists. The staining was performed on FFPE prostatectomy tumour blocks using Leica’s Polymer Refine Detection System (DS9800) in combination with their Bond automated system (Leica Biosystems Newcastle Ltd, Newcastle, UK) following a previously published protocol2. Briefly, sections were cut to 4 μm thickness and baked for 1 hour at 60°C ahead of deparaffinisation and rehydration, as standard, on the ST5020 Multistainer (Leica Biosystems). Subsequent immunohistochemical staining was carried out on Leica’s automated Bond III platform (Leica Biosystems) in conjunction with their Polymer Refine Detection System (Cat. No. DS9800, Leica Biosystems). Sections stained for HIF-1α were pre-treated with Epitope Retrieval Solution 1 (Cat. No. AR9961, Leica Biosystems) and those stained for other antibodies were pre-treated with Epitope Retrieval Solution 2 (Cat. No. AR9640, Leica Biosystems). Incubation was for 20 minutes at 99°C. Antibodies were diluted to 23.36 µg/mL (MCT1), 0.3 µg/mL (MPC1), 4.5 µg/mL (MPC2), 1:100 (FASN), 1:50 (AR), 0.6 µg/ml (HIF-1α), and 4.1 µg/ml (CD31), respectively. Endogenous peroxidase activity was quenched using 3-4% (v/v) hydrogen peroxide and primary antibody was detected using Anti-rabbit Poly-HRP-IgG (<25 μg/mL; part of Leica Biosystems Polymer Refine Detection System) containing 10% (v/v) animal serum in tris-buffered saline/0.09% ProClin 950.The complex was visualized using 66 mM 3,3’-Diaminobenzidine tetrahydrochloride hydrate in a stabilizer solution and ≤0.1% (v/v) Hydrogen Peroxide. DAB Enhancer (Cat. No. AR9432, Leica Biosystems) was used to intensify the signal. Cell nuclei were counterstained with <0.1% haematoxylin. As described previously2, HALO v3.3.2541.405 (Indica Labs, Albuquerque, NM, USA) Membrane v1.7 (MCT1), multiplex IHC v2.3.4 (MPC1, MPC2, HIF-1α), area quantification v2.2.1 (CD31), multiplex IHC v3.1.4 (AR), and multiplex IHC v2.3.4 (FASN) modules were used for automated analysis of scanned sections. Optical densities for weak, moderate, and strong stains were: MCT1, 0.1602, 0.2302, 0.4037; HIF-1α nuclear 0.1958, 0.7522, 0.885; CD31, 0.2164, 0.2721, 0.3832; AR nuclear, 0.0987,0.3761,0.544; FASN, 0.1881, 0.5949, 0.8407; MPC1, positive optical density threshold 0.0759; MPC2, positive optical density threshold 0.0758. Given the prior identification of the epithelial compartment as the source of clinical [1-13C]lactate labelling in the HP-13C-MRI cohort2, we first quantified epithelial cell density of each of the above proteins by dividing the number of positive epithelial cells by the ROI areas measured in mm2. This approach helped us achieve a closer comparison between the limited spatial resolution of absolute [1-13C]lactate signal quantification on HP-13C-MRI and the density of cells that could contribute to the observed signal. In addition, we also assessed the percentage of positive epithelial cells (MPC2, FASN) and nuclei (HIF-1α, AR) in individual epithelial glands harbouring GP3, non-cribriform GP4, and GP4 ICC disease to obtain a more refined metric of the biomarker expression at a cellular level. ### Spatial transcriptomics To infer epithelial capacity for LDH-catalysed [1-13C]pyruvate-to-[1-13C]lactate conversion, we used RNAscope to derive the total epithelial LDH density, measured as the total number of epithelial *LDHA* and *LDHB* mRNA copies per mm2. The analysis was performed using the RNAscope spatial transcriptomics technique according to the previously described protocol2. Briefly, FFPE sections were cut to 4 μm thickness and baked for 1 hour at 60°C before loading onto a Bond RX instrument (Leica Biosystems Newcastle Ltd, Newcastle, UK). Slides were deparaffinized and rehydrated on board prior to pre-treatments using Epitope Retrieval Solution 2 (Cat No. AR9640, Leica Biosystems) at 95°C for 15 minutes, and ACD Enzyme from the Multiplex Reagent kit at 40oC for 15 minutes. Probes (*LDHA*, *LDHB*) were visualised using Opal fluorophores diluted to 1:1000 using RNAscope LS Multiplex TSA Buffer. Probe hybridisation, signal amplification, and detection were all performed on the Bond Rx according to the ACD protocol. Slides were removed from the Bond Rx and mounted using Prolong Diamond (Cat. No. P36965, ThermoFisher Scientific, Watham, MA, USA). In all patients, simultaneous detection of human *LDHA* and *LDHB* was performed using Advanced Cell Diagnostics (ACD, Bio-Techne, Abingdon, UK) RNAscope 2.5 LS Multiplex Reagent Kit (Cat No. 322800), and RNAscope 2.5 LS probes (ACD, Hayward, CA, USA). The slides were imaged on the AxioScan (Carl-Zeiss-Stiftung, Stuttgart, Germany) to create whole-slide images. Images were captured at 40x magnification, with a resolution of 0.25 microns per pixel. HALO v3.2.1851.266 and the FISH v2.2.0 modules were used for the automated analysis of scanned RNAscope sections. ### Spatial metabolomics acquisition and analysis Fresh-frozen RARP specimens obtained from the DESI-MSI cohort were transferred on ice to the laboratory within 30min of surgical resection, with multiple punch biopsies collected as described above. For the DESI-MSI experiment, the punched cores (n = 117) were embedded and prepared according to a previously reported sample preparation workflow90. In brief, the punched cores were co-embedded in a (Hydroxypropyl)-methylcellulose (HPMC) and Polyvinylpyrrolidone (PVP) hydrogel to enable time-efficient sectioning under comparable conditions for all specimens analysed in one experiment. A total of 18 punched cores were placed upright in peel-away moulds (Thermo Scientific, Waltham, MA, USA) pre-filled with ice-cold embedding medium. Snap freezing of the filled mould was performed in dry ice-chilled isopropanol followed by a wash in dry ice chilled iso-pentane to wash off the excess of isopropanol. The frozen moulds were kept on dry ice to allow evaporation of the adherent iso-pentane before sectioning. The resulting blocks (n = 10) were sectioned to 10µm thickness using a CM3050 cryo-microtome (Leica Biosystems, Nussloch, Germany) and thaw-mounted onto Superfrost slides (Fisher Scientific, Loughborough, UK) for DESI-MSI and H&E histological examination. DESI-MSI was carried out using an automated 2D DESI source (Prosolia Inc, Indianapolis, IN, USA) with home-built sprayer assembly mounted to a Q-Exactive FTMS instrument (Thermo Scientific, Bremen, Germany). Analyses were performed at spatial resolutions of 65μm in negative ion mode and mass spectra were collected in the mass range of 80−600 Da with mass resolving power set to 70000 at *m/z* 200 and an S-Lens setting of 100. Methanol/water (95:5 v/v) was used as the electrospray solvent at a flow rate of 1.0 μL/min and a spray voltage of -4.5kV. Distance between DESI sprayer to MS inlet was 7mm, while distance between sprayer tip to sample surface was 1.5mm at an angle of 75°. Nitrogen N4.8 was used as nebulizing gas at a pressure of 6.5 bar. Omnispray 2D (Prosolia, Indianapolis, USA) and Xcalibur (Thermo Fisher Scientific) software were used for MS data acquisition. Individual line scans were converted into centroided. mzML format using MSConvert (ProteoWizard toolbox version 3.0.4043) and subsequently into .imzML using imzML converter v1.3. Raw metabolite spectra were extracted using the pyimzML python package and aligned to a common mass axis using a fixed bin size of 0.01 Da. Spectra were normalised by the root mean square of the intensity across the entire mass range, to compensate for signal instabilities and to facilitate comparison between experiments. The analysis of metabolites within pathologist-annotated H&E ROIs required the co-registration of whole-slide histological images with DESI-MSI data. To do so, we first extracted a low-resolution version of each whole-slide H&E-stained image, downsampling it by a factor of 32. A binary tissue-background mask was computed for each H&E slide by thresholding the saturation of the image using Otsu thresholding. We then removed all objects and holes in the binary mask smaller than 64 pixels. For the DESI-MSI data, we computed a total ion current (TIC) image, normalised to lie in the range 0-1. The H&E binary mask was then resized to have the same dimensions as the MSI TIC image. Using the resized H&E mask and the MSI TIC image, we computed an affine transform using the Elastix package, treating H&E mask as the moving image and the MSI image as the fixed image. Pathologist annotations were converted to binary masks, resized to the same dimensions as the MSI total ion current image, and then transformed using the previously calculated affine transformation. This yielded a mask identifying the pixels in the MSI image corresponding to the region annotated by the pathologist on the whole slide H&E image. We then computed a single spectrum per ROI by taking the mean of all MSI pixels inside the region. To resolve the metabolic phenotype of benign and tumour prostatic epithelial cells within the framework of this study, we used individual benign (n = 695) and tumour (n = 468) epithelial ROIs to derive relative abundance values of endogenous metabolites contributing to the three key metabolic pathways of interest (glycolysis, TCA cycle, fatty acid biosynthesis; individual KEGG metabolites are listed in **Supplementary Table 1**). The relative abundance values of all metabolites that have been mapped to the reference m/z values (within 5 ppm accuracy) in the KEGG database were compared between the groups of interest. The corresponding P values have been calculated for each metabolite using non-parametric Wilcoxon rank sum test. To perform the metabolic pathway enrichment analysis (MPEA), the metabolites enriched (p > 0.05; Log2FC > 1) in either of the comparison groups were subjected to overrepresentation analysis. Only the pathways with an overall size of more than two metabolites and having more than two enriched metabolites contributing to the named pathway (“hits”) were retained. The produced output of MPEA contains a P value based on a hypergeometric test, as well as the pathway enrichment score (PES) for each pathway. PES represents a log2 transformed ratio between the number of observed enriched metabolites in the dataset and the number of metabolites expected to be enriched by random chance. In addition to MPEA, we also sought to ascertain the individual contribution of constituent KEGG metabolites on the ability of DESI-MSI to distinguish between benign and tumour epithelial cells. To do so, we trained a small neural network to distinguish the spectra of epithelial ROIs annotated as either tumour or benign (**Fig. 4e**). The full dataset consistent of a total of 1163 ROIs (n = 695 benign and n = 468 tumour) that had also been used for the MPEA. Each sample, corresponding to an ROI, consisted of the mean value of 31 selected metabolites corresponding to KEGG glycolysis and glyconeogenesis, TCA cycle, and fatty acid synthesis metabolic pathways (**Supplementary Table 1**). As a first step, we normalised the data, such that the mean intensity of each metabolite across all ROIs was equal to 0, with a standard deviation of 1. A suite of models using 5-fold cross validation were trained. The splits were determined at the patient level, such that no samples from patients in the test set were present in the training set. This was done to ensure that the model learned features which distinguished benign regions from tumour regions and did not memorise patient-specific features common to both tumour and benign regions. The neural network had a multi-layer perceptron architecture, with 64 neurons in each of the two hidden layers. We used ReLu activations in the hidden layers, with a single neuron in the final layer which used a sigmoid activation. We used a binary cross entropy loss function, the Adam optimiser with a learning rate of 10-4 and we used a batch size of 32. To prevent overfitting, we applied data augmentation during training by adding a small amount of Gaussian noise to each sample. The amplitude of the noise was set to be 10% of the absolute value of each metabolite. In addition, we employed an early stopping protocol to determine the number of epochs for training, based on the validation loss. To evaluate the importance of its constituent features on the final model performance, we combined Shapley additive explanations (SHAP) for each test set in our 5-fold cross validation training scheme59. The median performance of the tissue classifier for discriminating between benign and tumour ROI classes across the five different networks was assessed using the area under the ROC curve. ### Open-source spatial transcriptomics data analysis To evaluate the generalisability of the key findings from the HP-13C-MRI cohort, we analysed publicly available mRNA data related to the expression of *SLC16A1, LDHA, LDHB, MPC2, FASN,* and *AR* in benign and tumour surgical specimens. First, bulk mRNA sequencing outputs (log-transformed FPKM values) of the TCGA-PRAD study43 were accessed through the NCI GDC Data Portal92 and compared between the two sample classes using the Mann-Whitney *U* test. To allow comparability with the HP-13C-MRI cohort, only data from patients with primary Gleason patterns of 3 and 4 were analysed. In addition, single-cell and spatial RNA-sequencing analyses were performed on two datasets (EGAS00001005787 and GSE176031). For each, data were all processed using the Seurat toolkit ([https://satijalab.org/seurat/](https://satijalab.org/seurat/))91–94. Data, obtained as matrices comprising feature (gene) counts as rows and barcoded cells as columns (with annotations regarding cell type and other metadata as available), were first preprocessed by filtering out cells with over 2,500 or under 200 unique feature/gene counts as well as those with >5% of counts including mitochondrial DNA. Data were then log-normalised and scaled by a factor of 10,000. Highly variable features were selected using Seurat’s ‘vst ’selection method, and a linear scaling transformation was performed before dimensionality reduction was performed via principal component analysis (PCA). Jackstraw plots were generated95 to select an appropriate number of dimensions, which was also validated by ranking principal components by the percentage of variance explained by adding each principal component. The “appropriate” number of dimensions was that at which no appreciable increase in signal/percentage of variance explained was obtained by adding an additional dimension. Cells were subsequently clustered using the Louvain clustering algorithm96, set to a clustering resolution of 0.5. Nonlinear dimensionality reduction was performed with the UMAP method97. Differential expression among clusters was analysed using Seurat’s default differential expression analysis function, ‘FindAllMarkers’. Finally, for analysis of gene signatures comprising multiple genes, a new “meta-feature” was created using Seurat’s ‘MetaFeature ’function, which calculates the relative contribution of a gene set to each cell in the dataset. For integration of spatial transcriptomics data, 10x Visium data were used and loaded to R ([https://www.r-project.org/](https://www.r-project.org/)) as an object comprising an image of the tissue slice along with spot-by-spot gene expression data. As above, data were preprocessed by normalisation to account for spot-by-spot differences in sequencing depth. SCTransform98 was used to normalise data and detect high-variance features. ### Statistical analysis Statistical analyses were conducted using GraphPad Prism (version 9.0.2, GraphPad Software, San Diego, CA, USA). Normal distribution of the data was assessed using the D’Agostino-Pearson test (threshold *P* ≥ 0.05). All intergroup comparisons were performed using the two-sided Mann-Whitney *U* test. Correlation analysis was conducted using the Spearman’s rank correlation test since at least one variable was always non-normally distributed. In the HP-13C-MRI cohort, intergroup comparison of HP-13C-MRI-derived [1-13C]pyruvate and [1-13C]lactate labelling was conducted between n = 13 HP-13C-MRI-visible tumours (total carbon SNR threshold > 5.0 according to the Rose criterion85,99) and n = 13 contralateral areas of benign tissue. Simultaneously, due to interval androgen deprivation therapy (ADT) administered to patient 5 in **Table 1**, his two lesions were excluded from intergroup comparison of tissue-based parameters, with the remaining sample size including n = 11 tumour and n = 11 benign areas. In addition, values of MPC1 and MPC2 epithelial density were log-transformed due to their high range. In the DESI-MSI cohort, the comparison of whole-core epithelial lactate abundance, tissue density, epithelial and stromal cell fractions, as well as compartmentalised lactate pool was made between cores including purely benign (n = 61) and >50% tumour (n = 38) tissues (**Extended Data Fig. 5**). MPEA and deep-learning-based tissue classifier were built using the ROIs derived from all tissue cores included in the study, including n = 23 mixed cores. For intergroup comparisons and correlation analyses, no multiplicity correction was applied, and therefore all significant tests should be interpreted as exploratory rather than confirmatory. ## Supporting information Supplementary Table 1 [[supplements/296228_file02.pdf]](pending:yes) ## Data Availability All data produced in the present study are available upon reasonable request to the authors and subject to appropriate data transfer agreements. ## Extended Data Figures ![Extended Data Fig. 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/10/02/2023.09.29.23296228/F8.medium.gif) [Extended Data Fig. 1.](http://medrxiv.org/content/early/2023/10/02/2023.09.29.23296228/F8) Extended Data Fig. 1. Relationship between tumour epithelial cell density and clinical [1-13C]lactate labelling. Spearman’s correlation plot comparing tumour-derived HP-13C-MRI-derived [1-13C]lactate SNR and digital-pathology-derived epithelial cell density quantified using the random forest tissue classifier described in the main manuscript and prior publications1,2. The plot includes the rank correlation coefficient and a corresponding *P* value. The analysis included data from n = 11 HP-13C-MRI-visible lesions. ![Extended Data Fig. 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/10/02/2023.09.29.23296228/F9.medium.gif) [Extended Data Fig. 2.](http://medrxiv.org/content/early/2023/10/02/2023.09.29.23296228/F9) Extended Data Fig. 2. Relationship between tumour vascular density, perfusion and permeability, and clinical [1-13C]lactate labelling. **a,** Spearman’s correlation plot comparing tumour-derived HP-13C-MRI-derived [1-13C]lactate SNR and 1H-MRI-derived *K*trans. **b,** Spearman’s correlation plot comparing tumour-derived HP-13C-MRI-derived [1-13C]lactate SNR and tissue-based CD31 density. Both panels include individual rank correlation coefficients and corresponding *P* values. Both analyses included data from n = 11 HP-13C-MRI-visible lesions. ![Extended Data Fig. 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/10/02/2023.09.29.23296228/F10.medium.gif) [Extended Data Fig. 3.](http://medrxiv.org/content/early/2023/10/02/2023.09.29.23296228/F10) Extended Data Fig. 3. Segmentation of benign and malignant prostate areas in the HP-13C-MRI cohort. Representative whole-mount pathology, 1H-MRI-derived ADC and *K*trans maps, along with HP-13C-MRI-derived [1-13C]pyruvate and [1-13C]lactate SNR maps overlaid on conventional T2-weighted images with overlaid tumour (red) and contralateral benign (green) ROIs. Drawn by an expert genitourinary pathologist, these ROIs were transposed onto IHC and RNAscope images to derive the expression of proteins and mRNA molecules of interest. ![Extended Data Fig. 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/10/02/2023.09.29.23296228/F11.medium.gif) [Extended Data Fig. 4.](http://medrxiv.org/content/early/2023/10/02/2023.09.29.23296228/F11) Extended Data Fig. 4. Single-gland segmentation of benign and malignant prostate in the HP-13C-MRI cohort. Representative H&E images including standardised ROIs encompassing areas of pure benign prostatic epithelium, as well as GP3 PCa, non-cribriform GP4 disease (glomeruloid glands), and cribriform GP4 PCa (invasive cribriform carcinoma). These ROIs were also transposed onto IHC and RNAscope slides to derive biomarker expression specific to a particular histological entity. ![Extended Data Fig. 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/10/02/2023.09.29.23296228/F12.medium.gif) [Extended Data Fig. 5.](http://medrxiv.org/content/early/2023/10/02/2023.09.29.23296228/F12) Extended Data Fig. 5. Fresh-frozen specimen classification in the DESI-MSI cohort. Representative H&E images of benign (0% tumour tissue), mixed (<50% tumour tissue), and tumour (>80% tumour tissue) fresh-frozen specimen from the spatial metabolomics cohort. Scale bars denote 1mm. ![Extended Data Fig. 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/10/02/2023.09.29.23296228/F13.medium.gif) [Extended Data Fig. 6.](http://medrxiv.org/content/early/2023/10/02/2023.09.29.23296228/F13) Extended Data Fig. 6. Tissue classifiers used in this study. **a,** H&E image with epithelial (red) and stromal (green) ROIs used to train the random forest classifier used in the HP-13C-MRI cohort. **b,** H&E images with stromal (green), epithelial (red), and luminal (yellow) ROIs used to train the DenseNet tissue classifier used in the DESI-MSI cohort. ## Footnotes * †† These authors jointly supervised the work * Received September 29, 2023. * Revision received September 29, 2023. * Accepted October 2, 2023. * © 2023, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. 1.Kurhanewicz, J. et al. Hyperpolarized 13C MRI: Path to Clinical Translation in Oncology. Neoplasia (United States*)* vol. 21 1–16 (2019). 2. 2.Sushentsev, N. et al. Hyperpolarised 13C-MRI identifies the emergence of a glycolytic cell population within intermediate-risk human prostate cancer. Nat. Commun. *2022 131* 13, 1–12 (2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-021-27838-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 3. 3.Nelson, S. J. et al. Metabolic imaging of patients with prostate cancer using hyperpolarized [1-13C]pyruvate. Science Translational Medicine vol. 5 198ra108 (2013). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6InNjaXRyYW5zbWVkIjtzOjU6InJlc2lkIjtzOjE0OiI1LzE5OC8xOThyYTEwOCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzEwLzAyLzIwMjMuMDkuMjkuMjMyOTYyMjguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 4. 4.Granlund, K. L. et al. Hyperpolarized MRI of Human Prostate Cancer Reveals Increased Lactate with Tumor Grade Driven by Monocarboxylate Transporter 1. Cell Metab. 31, 105–114.e3 (2020). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31564440&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 5. 5.Chen, H. Y. et al. Hyperpolarized 13C-pyruvate MRI detects real-time metabolic flux in prostate cancer metastases to bone and liver: a clinical feasibility study. Prostate Cancer Prostatic Dis. 23, 269–276 (2020). 6. 6.Aggarwal, R., Vigneron, D. B. & Kurhanewicz, J. Hyperpolarized 1-[13C]-Pyruvate Magnetic Resonance Imaging Detects an Early Metabolic Response to Androgen Ablation Therapy in Prostate Cancer. European Urology vol. 72 1028–1029 (2017). 7. 7.Kouchkovsky, I. de et al. Hyperpolarized 1-[13C]-Pyruvate Magnetic Resonance Imaging Detects an Early Metabolic Response to Immune Checkpoint Inhibitor Therapy in Prostate Cancer. Eur. Urol. , (2021). 8. 8.Gallagher, F. A. et al. Imaging breast cancer using hyperpolarized carbon-13 MRI. Proc. Natl. Acad. Sci. U. S. A. 117, 2092–2098 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMDoiMTE3LzQvMjA5MiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzEwLzAyLzIwMjMuMDkuMjkuMjMyOTYyMjguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 9. 9.Woitek, R. et al. Hyperpolarized Carbon-13 MRI for Early Response Assessment of Neoadjuvant Chemotherapy in Breast Cancer Patients. Cancer Res. 81, 6004 (2021). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiY2FucmVzIjtzOjU6InJlc2lkIjtzOjEwOiI4MS8yMy82MDA0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMTAvMDIvMjAyMy4wOS4yOS4yMzI5NjIyOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 10. 10.Ursprung, S. et al. Hyperpolarized 13C-Pyruvate Metabolism as a Surrogate for Tumor Grade and Poor Outcome in Renal Cell Carcinoma—A Proof of Principle Study. Cancers *2022*, Vol. 14, *Page* 335 **14**, 335 (2022). 11. 11.O’Connor, J. P. B. et al. Imaging biomarker roadmap for cancer studies. Nat. Rev. Clin. Oncol. *2016 143* 14, 169–186 (2016). 12. 12.Preisser, F. et al. Intermediate-risk Prostate Cancer: Stratification and Management. *Eur*. Urol. Oncol. 3, 270–280 (2020). 13. 13.Kweldam, C. F. et al. Cribriform growth is highly predictive for postoperative metastasis and disease-specific death in Gleason score 7 prostate cancer. Mod. Pathol. 28, 457–464 (2015). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 14. 14.van Leenders, G. J. L. H., et al. The 2019 International Society of Urological Pathology (ISUP) Consensus Conference on Grading of Prostatic Carcinoma. Am. J. Surg. Pathol. 44, e87–e99 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/PAS.0000000000001497&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 15. 15.Wong, H. Y. et al. Single cell analysis of cribriform prostate cancer reveals cell intrinsic and tumor microenvironmental pathways of aggressive disease. Nat. Commun. 2022 131 13, 1–21 (2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-021-27838-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 16. 16.Seasor, T. et al. Tumor transcriptomic profile in patients (pts) with advanced prostate cancer (PCa) with cribriform and signet ring-like histologies. doi:10.1200/JCO.2023.41.6_suppl.259 41, 259–259 (2023). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1200/JCO.2023.41.6_suppl.259&link_type=DOI) 17. 17.Chappidi, M. R. et al. Transcriptomic Heterogeneity of Expansile Cribriform and Other Gleason Pattern 4 Prostate Cancer Subtypes. *Eur*. Urol. Oncol. (2023) doi:10.1016/J.EUO.2023.06.007. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/J.EUO.2023.06.007&link_type=DOI) 18. 18.Truong, M. et al. A Comprehensive Analysis of Cribriform Morphology on Magnetic Resonance Imaging/Ultrasound Fusion Biopsy Correlated with Radical Prostatectomy Specimens. J. Urol. 199, 106–113 (2018). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 19. 19.Bankson, J. A. et al. Kinetic Modeling and Constrained Reconstruction of Hyperpolarized [1-13C]-Pyruvate Offers Improved Metabolic Imaging of Tumors. Cancer Res. 75, 4708 (2015). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiY2FucmVzIjtzOjU6InJlc2lkIjtzOjEwOiI3NS8yMi80NzA4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMTAvMDIvMjAyMy4wOS4yOS4yMzI5NjIyOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 20. 20.Gordon, J. W., Vigneron, D. B. & Larson, P. E. Z. Development of a symmetric echo planar imaging framework for clinical translation of rapid dynamic hyperpolarized 13 C imaging. Magn. Reson. Med. 77, 826–832 (2017). 21. 21.Park, I., Lupo, J. M. & Nelson, S. J. Correlation of Tumor Perfusion Between Carbon-13 Imaging with Hyperpolarized Pyruvate and Dynamic Susceptibility Contrast MRI in Pre-Clinical Model of Glioblastoma. Mol. Imaging Biol. 21, 626–632 (2019). 22. 22.Harris, T., Eliyahu, G., Frydman, L. & Degani, H. Kinetics of hyperpolarized 13C1-pyruvate transport and metabolism in living human breast cancer cells. Proc. Natl. Acad. Sci. U. S. A. 106, 18131–18136 (2009). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTA2LzQzLzE4MTMxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMTAvMDIvMjAyMy4wOS4yOS4yMzI5NjIyOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 23. 23.Rao, Y. et al. Hyperpolarized [1-13C]pyruvate-to-[1-13C]lactate conversion is rate-limited by monocarboxylate transporter-1 in the plasma membrane. Proc. Natl. Acad. Sci. U. S. A. 117, 22378–22389 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzM2LzIyMzc4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMTAvMDIvMjAyMy4wOS4yOS4yMzI5NjIyOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 24. 24.Sriram, R. et al. Elevated Tumor Lactate and Efflux in High-grade Prostate Cancer demonstrated by Hyperpolarized 13C Magnetic Resonance Spectroscopy of Prostate Tissue Slice Cultures. Cancers (Basel*).* 12, 537 (2020). 25. 25.Wang, Z. J. et al. Hyperpolarized 13C MRI: State of the Art and Future Directions. Radiology 291, 273–284 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1148/radiol.2019182391&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30835184&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 26. 26.Brindle, K. M., Bohndiek, S. E., Gallagher, F. A. & Kettunen, M. I. Tumor imaging using hyperpolarized 13C magnetic resonance spectroscopy. Magn. Reson. Med. 66, 505–519 (2011). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mrm.22999&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21661043&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 27. 27.Sun, C. et al. Spatially resolved multi-omics highlights cell-specific metabolic remodeling and interactions in gastric cancer. Nat. Commun. 14, 2692 (2023). 28. 28.Bruinen, A. L. et al. Mass Spectrometry Imaging of Drug Related Crystal-Like Structures in Formalin-Fixed Frozen and Paraffin-Embedded Rabbit Kidney Tissue Sections. J. Am. Soc. Mass Spectrom. 27, 117–123 (2016). 29. 29.Goodwin, R. J. A. Sample preparation for mass spectrometry imaging: small mistakes can lead to big consequences. J. Proteomics. 75, 4893–4911 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jprot.2012.04.012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22554910&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 30. 30.Chowdhury, R. et al. Quantification of Prostate Cancer Metabolism Using 3D Multiecho bSSFP and Hyperpolarized [1-13 C] Pyruvate: Metabolism Differs Between Tumors of the Same Gleason Grade. J. Magn. Reson. Imaging (2022) doi:10.1002/JMRI.28467. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/JMRI.28467&link_type=DOI) 31. 31.Tessem, M. B. et al. Evaluation of Lactate and Alanine as Metabolic Biomarkers of Prostate Cancer Using 1H HR-MAS Spectroscopy of Biopsy Tissues. Magn. Reson. Med. 60, 510 (2008). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mrm.21694&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18727052&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000259053900002&link_type=ISI) 32. 32.Albers, M. J. et al. Hyperpolarized 13C Lactate, Pyruvate, and Alanine: Noninvasive Biomarkers for Prostate Cancer Detection and Grading. Cancer Res. 68, 8607–8615 (2008). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiY2FucmVzIjtzOjU6InJlc2lkIjtzOjEwOiI2OC8yMC84NjA3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMTAvMDIvMjAyMy4wOS4yOS4yMzI5NjIyOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 33. 33.Bailey, C. et al. VERDICT MRI validation in fresh and fixed prostate specimens using patient-specific moulds for histological and MR alignment. NMR Biomed. 32, (2019). 34. 34.Ippolito, L. et al. Lactate Rewires Lipid Metabolism and Sustains a Metabolic-Epigenetic Axis in Prostate Cancer. Cancer Res. 82, 1267–1282 (2022). 35. 35.Andersen, M. K. et al. Spatial differentiation of metabolism in prostate cancer tissue by MALDI-TOF MSI. Cancer Metab. 2021 91 9, 1–13 (2021). 36. 36.Andersen, S. et al. Organized metabolic crime in prostate cancer: The coexpression of MCT1 in tumor and MCT4 in stroma is an independent prognosticator for biochemical failure. Urol. Oncol. Semin. Orig. Investig. (2015) doi:10.1016/j.urolonc.2015.05.013. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.urolonc.2015.05.013&link_type=DOI) 37. 37.Zaichick, V. & Zaichick, S. Age-related histological and zinc content changes in adult nonhyperplastic prostate glands. Age (Omaha*).* 36, 167 (2014). 38. 38.Cuenod, C. A. & Balvay, D. Perfusion and vascular permeability: Basic concepts and measurement in DCE-CT and DCE-MRI. Diagn. Interv. Imaging 94, 1187–1204 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.diii.2013.10.010&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24211260&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 39. 39.Oto, A. et al. Diffusion-weighted and dynamic contrast-enhanced MRI of prostate cancer: correlation of quantitative MR parameters with Gleason score and tumor angiogenesis. AJR. Am. J. Roentgenol. 197, 1382–1390 (2011). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2214/AJR.11.6861&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22109293&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 40. 40.Pertega-Gomes, N. et al. A glycolytic phenotype is associated with prostate cancer progression and aggressiveness: a role for monocarboxylate transporters as metabolic targets for therapy. J. Pathol. 236, 517–530 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/path.4547&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25875424&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 41. 41.Pértega-Gomes, N. et al. Monocarboxylate transporter 4 (MCT4) and CD147 overexpression is associated with poor prognosis in prostate cancer. BMC Cancer 11, (2011). 42. 42.Abeshouse, A. et al. The molecular taxonomy of primary prostate cancer. Cell 63, 1011 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2015.10.025&link_type=DOI) 43. 43.Latonen, L. et al. Integrative proteomics in prostate cancer uncovers robustness against genomic and transcriptomic aberrations during disease progression. Nat. Commun. *2018 91* 9, 1–13 (2018). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-018-02974-x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29317637&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 44. 44.Costello, L. C. & Franklin, R. B. A comprehensive review of the role of zinc in normal prostate function and metabolism$\mathsemicolon$ and its implications in prostate cancer. Arch. Biochem. Biophys. 611, 100–112 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.abb.2016.04.014&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27132038&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 45. 45.Costello, L. C. & Franklin, R. B. Concepts of citrate production and secretion by prostate 1. Metabolic relationships. Prostate (1991) doi:10.1002/pros.2990180104. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/pros.2990180104&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1987578&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 46. 46.Bader, D. A. & McGuire, S. E. Tumour metabolism and its unique properties in prostate adenocarcinoma. Nat. Rev. Urol. *2020 174* 17, 214–231 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41585-020-0288-x&link_type=DOI) 47. 47.Swinnen, J. V., Van Veldhoven, P. P., Esquenet, M., Heyns, W. & Verhoeven, G. Androgens markedly stimulate the accumulation of neutral lipids in the human prostatic adenocarcinoma cell line LNCaP. Endocrinology 137, 4468– 4474 (1996). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/en.137.10.4468&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8828509&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996VN28800053&link_type=ISI) 48. 48.Swinnen, J. V et al. Selective activation of the fatty acid synthesis pathway in human prostate cancer. Int. J. Cancer 88, 176–179 (2000). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/1097-0215(20001015)88:2<176::AID-IJC5>3.0.CO;2-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11004665&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000089629700005&link_type=ISI) 49. 49.Bader, D. A. et al. Mitochondrial pyruvate import is a metabolic vulnerability in androgen receptor-driven prostate cancer. Nat. Metab. 1, 70–85 (2019). 50. 50.Massie, C. E. et al. The androgen receptor fuels prostate cancer by regulating central metabolism and biosynthesis. EMBO J. 30, 2719–2733 (2011). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZW1ib2pubCI7czo1OiJyZXNpZCI7czoxMDoiMzAvMTMvMjcxOSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzEwLzAyLzIwMjMuMDkuMjkuMjMyOTYyMjguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 51. 51.Brodie, C. Overcoming Autofluorescence (AF) and Tissue Variation in Image Analysis of In Situ Hybridization. Methods Mol. Biol. 2148, 19–32 (2020). 52. 52.Tuong, Z. K. et al. Resolving the immune landscape of human prostate at a single-cell level in health and cancer. Cell Rep. 37, 110132 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.celrep.2021.110132&link_type=DOI) 53. 53.Song, H. et al. Single-cell analysis of human primary prostate cancer reveals the heterogeneity of tumor-associated epithelial cell states. Nat. Commun. *2022 131* 13, 1–20 (2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-021-27838-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 54. 54.Erickson, A. et al. Spatially resolved clonal copy number alterations in benign and malignant tissue. Nat. *2022 6087922* 608, 360–367 (2022). 55. 55.Fernandez-Caggiano, M. et al. Mitochondrial pyruvate carrier abundance mediates pathological cardiac hypertrophy. Nat. Metab. *2020 211* 2, 1223–1231 (2020). 56. 56.Yang, J., Kalhan, S. C. & Hanson, R. W. What Is the Metabolic Role of Phosphoenolpyruvate Carboxykinase? J. Biol. Chem. 284, 27025 (2009). [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEyOiIyODQvNDAvMjcwMjUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8xMC8wMi8yMDIzLjA5LjI5LjIzMjk2MjI4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 57. 57.Chetan, M. R., Barrett, T. & Gallagher, F. A. Clinical significance of prostate (18)F-labelled fluorodeoxyglucose uptake on positron emission tomography/computed tomography: A five-year review. World J. Radiol. 9, 350–358 (2017). 58. 58.O’Connor, T. P., Roche, T. E. & Paukstelis, J. V. 13C nuclear magnetic resonance study of the pyruvate dehydrogenase-catalyzed acetylation of dihydrolipoamide. J. Biol. Chem. 257, 3110–3112 (1982). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEwOiIyNTcvNi8zMTEwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMTAvMDIvMjAyMy4wOS4yOS4yMzI5NjIyOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 59. 59.Lundberg, S. & Lee, S.-I. A Unified Approach to Interpreting Model Predictions. (2017). 60. 60.Motulsky, H. J. & Brown, R. E. Detecting outliers when fitting data with nonlinear regression - A new method based on robust nonlinear regression and the false discovery rate. BMC Bioinformatics 7, 1–20 (2006). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2105-7-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16393334&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000236027800001&link_type=ISI) 61. 61.Miyai, K. et al. Histological differences in cancer cells, stroma, and luminal spaces strongly correlate with in vivo MRI-detectability of prostate cancer. Mod. Pathol. 32, 1536–1543 (2019). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 62. 62.Woitek, R. et al. Hyperpolarized 13C MRI of Tumor Metabolism Demonstrates Early Metabolic Response to Neoadjuvant Chemotherapy in Breast Cancer. Radiol. Imaging cancer 2, (2020). 63. 63.Barrett, T. et al. DCE and DW MRI in monitoring response to androgen deprivation therapy in patients with prostate cancer: a feasibility study. Magn. Reson. Med. 67, 778–785 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mrm.23062&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22135228&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 64. 64.Hötker, A. M. et al. Prostate Cancer: Assessing the effects of androgen-deprivation therapy using quantitative multi-parametric MRI. Eur. Radiol. 25, 2665 (2015). 65. 65.Reuter, V. E. Pathological changes in benign and malignant prostatic tissue following androgen deprivation therapy. Urology 49, 16–22 (1997). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0090-4295(97)00164-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9123731&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997WQ71200003&link_type=ISI) 66. 66.Xie, N. et al. NAD+ metabolism: pathophysiologic mechanisms and therapeutic potential. Signal Transduct. Target. Ther. 5, (2020). 67. 67.Böttcher, R. et al. Cribriform and intraductal prostate cancer are associated with increased genomic instability and distinct genomic alterations. BMC Cancer 18, (2018). 68. 68.van Leenders, G. J. L. H. et al. Improved Prostate Cancer Biopsy Grading by Incorporation of Invasive Cribriform and Intraductal Carcinoma in the 2014 Grade Groups. Eur. Urol. 77, 191–198 (2020). 69. 69.Tom, M. C. et al. Impact of Cribriform Pattern and Intraductal Carcinoma on Gleason 7 Prostate Cancer Treated with External Beam Radiotherapy. J. Urol. 202, 710–715 (2019). 70. 70.Ericson, K. J. et al. Diagnostic Accuracy of Prostate Biopsy for Detecting Cribriform Gleason Pattern 4 Carcinoma and Intraductal Carcinoma in Paired Radical Prostatectomy Specimens: Implications for Active Surveillance. J. Urol. 203, 311–317 (2020). 71. 71.Mottet, N. et al. EAU-EANM-ESTRO-ESUR-SIOG Guidelines on Prostate Cancer—2020 Update. Part 1: Screening, Diagnosis, and Local Treatment with Curative Intent. Eur. Urol. 79, 243–262 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.eururo.2020.09.042&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33172724&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 72. 72.Norris, J. M. et al. Conspicuity of cribriform prostate cancer on multiparametric magnetic resonance imaging: the jury is still out. BJU Int. 127, 169–170 (2021). 73. 73.Bader, D. A. & McGuire, S. E. Tumour metabolism and its unique properties in prostate adenocarcinoma. Nature Reviews Urology vol. 17 214–231 (2020). 74. 74.Ahmad, F., Cherukuri, M. K. & Choyke, P. L. Metabolic reprogramming in prostate cancer. Br. J. Cancer 125, 1185–1196 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41416-021-01435-5&link_type=DOI) 75. 75.Gao, J. et al. Diagnostic performance of 68Ga-PSMA PET/CT for identification of aggressive cribriform morphology in prostate cancer with whole-mount sections. Eur. J. Nucl. Med. Mol. Imaging 46, 1531–1541 (2019). 76. 76.Heetman, J. G. et al. Is cribriform pattern in prostate biopsy a risk factor for metastatic disease on 68Ga-PSMA-11 PET/CT? World J. Urol. (2023) doi:10.1007/S00345-023-04467-Z. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/S00345-023-04467-Z&link_type=DOI) 77. 77.Rüschoff, J. H. et al. What’s behind 68Ga-PSMA-11 uptake in primary prostate cancer PET? Investigation of histopathological parameters and immunohistochemical PSMA expression patterns. Eur. J. Nucl. Med. Mol. Imaging 48, 4042–4053 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00259-021-05501-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34386839&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 78. 78.Chappidi, M. R. et al. Transcriptomic Heterogeneity of Expansile Cribriform and Other Gleason Pattern 4 Prostate Cancer Subtypes. *Eur*. Urol. Oncol. (2023) doi:10.1016/J.EUO.2023.06.007. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/J.EUO.2023.06.007&link_type=DOI) 79. 79.Mutuku, S. M. et al. Unravelling Prostate Cancer Heterogeneity Using Spatial Approaches to Lipidomics and Transcriptomics. Cancers (Basel*).* 14, 1702 (2022). 80. 80.Banerjee, S. et al. Diagnosis of prostate cancer by desorption electrospray ionization mass spectrometric imaging of small metabolites and lipids. Proc. Natl. Acad. Sci. U. S. A. 114, 3334–3339 (2017). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTE0LzEzLzMzMzQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8xMC8wMi8yMDIzLjA5LjI5LjIzMjk2MjI4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 81. 81.Andersen, M. K. et al. Spatial differentiation of metabolism in prostate cancer tissue by MALDI-TOF MSI. Cancer Metab. *2021 91* 9, 1–13 (2021). 82. 82.Butler, L. M. et al. Lipidomic profiling of clinical prostate cancer reveals targetable alterations in membrane lipid composition. Cancer Res. 81, 4981– 4993 (2021). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiY2FucmVzIjtzOjU6InJlc2lkIjtzOjEwOiI4MS8xOS80OTgxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMTAvMDIvMjAyMy4wOS4yOS4yMzI5NjIyOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 83. 83.Randall, E. C. et al. Molecular Characterization of Prostate Cancer with Associated Gleason Score Using Mass Spectrometry Imaging. Mol. Cancer Res. 17, 1155–1165 (2019). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToibW9sY2FucmVzIjtzOjU6InJlc2lkIjtzOjk6IjE3LzUvMTE1NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzEwLzAyLzIwMjMuMDkuMjkuMjMyOTYyMjguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 84. 84.Barrett, T. et al. Molecular imaging of the prostate: Comparing total sodium concentration quantification in prostate cancer and normal tissue using dedicated 13 C and 23 Na endorectal coils. J. Magn. Reson. Imaging 51, 90–97 (2020). 85. 85.Sushentsev, N. et al. The potential of hyperpolarised 13C-MRI to target glycolytic tumour core in prostate cancer. Eur. Radiol. 1–8 (2022) doi:10.1007/s00330-022-08929-7. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00330-022-08929-7&link_type=DOI) 86. 86.Lawrence, E. M. et al. Evaluating prostate cancer using fractional tissue composition of radical prostatectomy specimens and pre-operative diffusional kurtosis magnetic resonance imaging. PLoS One 11, (2016). 87. 87.Van Der Kwast, T. H. et al. ISUP Consensus Definition of Cribriform Pattern Prostate Cancer. Am. J. Surg. Pathol. 45, 1118–1126 (2021). 88. 88.Warren, A. Y. et al. Method for sampling tissue for research which preserves pathological data in radical prostatectomy. Prostate 73, 194–202 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/pros.22556&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22806573&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000312810900010&link_type=ISI) 89. 89.Bankhead, P. et al. QuPath: Open source software for digital pathology image analysis. Sci. Reports *2017 71* 7, 1–7 (2017). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/srep39967&link_type=DOI) 90. 90.Dannhorn, A. et al. Universal Sample Preparation Unlocking Multimodal Molecular Tissue Imaging. Anal. Chem. 92, 11080–11088 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1021/acs.analchem.0c00826&link_type=DOI) 91. 91.Hao, Y. et al. Integrated analysis of multimodal single-cell data. Cell 184, 3573–3587.e29 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/J.CELL.2021.04.048&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 92. 92.Stuart, T. et al. Comprehensive Integration of Single-Cell Data. Cell 177, 1888–1902.e21 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2019.05.031&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 93. 93.Butler, A., Hoffman, P., Smibert, P., Papalexi, E. & Satija, R. Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat. Biotechnol. 36, 411–420 (2018). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nbt.4096&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29608179&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 94. 94.Satija, R., Farrell, J. A., Gennert, D., Schier, A. F. & Regev, A. Spatial reconstruction of single-cell gene expression data. Nat. Biotechnol. 33, 495– 502 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nbt.3192&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25867923&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 95. 95.Macosko, E. Z. et al. Highly Parallel Genome-wide Expression Profiling of Individual Cells Using Nanoliter Droplets. Cell 161, 1202–1214 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2015.05.002&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26000488&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 96. 96.Nguyen, L. Van et al. Fast unfolding of communities in large networks. J. Stat. Mech. Theory Exp. 2008, P10008 (2008). 97. 97.Becht, E. et al. Dimensionality reduction for visualizing single-cell data using UMAP. Nat. Biotechnol. 37, 38–47 (2018). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nbt.4314&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F02%2F2023.09.29.23296228.atom) 98. 98.Hafemeister, C. & Satija, R. Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression. Genome Biol. 20, (2019). 99. 99.Burgess, A. E. The Rose model, revisited. JOSA A, *Vol.* 16, *Issue* 3, *pp.* 633-646 16, 633–646 (1999).