Plasma Metabolome Profiling for the Diagnosis of Catecholamine Producing Tumors

Context: Pheochromocytomas and paragangliomas (PPGL) cause catecholamine excess leading to a characteristic clinical phenotype. Intra-individual changes at metabolome level have been described after surgical PPGL removal. The value of metabolomics for the diagnosis of PPGL has not been studied yet.

Objective: Evaluation of quantitative metabolomics as a diagnostic tool for PPGL.

Design: Targeted metabolomics by liquid chromatography-tandem mass spectrometry of plasma specimens and statistical modeling using ML-based feature selection approaches in a clinically well characterized cohort study.

Patients: Prospectively enrolled patients (n=36, 17 female) from the Prospective Monoamine-producing Tumor Study (PMT) with hormonally active PPGL and 36 matched controls in whom PPGL was rigorously excluded.

Results: Among 188 measured metabolites, only without considering false discovery rate, 4 exhibited statistically significant differences between patients with PPGL and controls (histidine p=0.004, threonine p=0.008, lyso PC a C28:0 p=0.044, sum of hexoses p=0.018). Weak, but significant correlations for histidine, threonine and lyso PC a C28:0 with total urine catecholamine levels were identified. Only the sum of hexoses (reflecting glucose) showed significant correlations with plasma metanephrines.

By using ML-based feature selection approaches, we identified diagnostic signatures which all exhibited low accuracy and sensitivity. The best predictive value (sensitivity 87.5%, accuracy 67.3%) was obtained by using Gradient Boosting Machine Modelling.

Conclusions: The diabetogenic effect of catecholamine excess dominates the plasma metabolome in PPGL patients. While curative surgery for PPGL led to normalization of catecholamine-induced alterations of metabolomics in individual patients, plasma metabolomics are not useful for diagnostic purposes, most likely due to inter-individual variability.

Introduction

Pheochromocytomas and paragangliomas (PPGL) are defined as catecholamine-producing tumors that arise from chromaffin cells (1). Pheochromocytomas represent more than 80% of all PPGL and are located in the adrenal medulla whereas paraganglioma arise from paravertebral sympathetic ganglia and are most frequently located in the abdomen, chest, and pelvis (2). Paragangliomas deriving from parasympathetic tissue in the head and neck rarely produce hormones (1–3). Predisposing germline mutations, extra-adrenal location, and dopaminergic phenotype are the most relevant risk factors for malignancy (4, 5). Current data suggest that germline mutations are present in up to 40% of all patients with PPGL, with 18 susceptibility genes identified so far. Mutations are most frequently found in genes encoding subunits of succinate dehydrogenase (SDH), von Hippel-Lindau gene (VHL) and rearranged during transfection (RET) gene. As the presence of a germ line mutation was found to be an important factor of prognosis of affected patients, testing is recommended (5–8).

Catecholamine excess leads to a variety of well-known but unspecific symptoms such as hypertension, palpitation, headache and pallor (1, 9) and causes cardio- and cerebrovascular complications. The measurement of plasma free metanephrine (MN), normetanephrine (NMN) and methoxytyramine (MTY) is now a cornerstone of diagnosis and follow-up in clinical practice, providing high diagnostic accuracy when adequate pre-analytics, analytics and reference ranges are applied (2, 10, 11). In recent years mass spectrometry has become the gold standard due to its high analytic sensitivity and specificity (12) not limited on quantification of established markers but showing additionally its usefulness to identify tissue metabolomic profiles via MALDI-MSI (13). Nevertheless, the diagnosis of PPGL remains challenging and is often delayed due to lack of consideration of PPGL (4). In addition there is a high risk of false positive test results when strict pre-analytical conditions are not followed.

Metabolomics is the screening for characteristic substances in body fluids and tissue, which serve as direct marker of biochemical activity because they are not exposed to epigenetic regulation and post-translational modifications like proteins or genes and therefore reflect the individual phenotype (14). Untargeted metabolomics allows the identification of numerous molecules without prior knowledge of their presence in predefined groups but has the disadvantage of generating mostly qualitative information on target molecules. On the other hand, quantitation of previously specified molecules is possible by targeted metabolomics. Still, the number of metabolites is typically limited to substances that are precisely characterized by their chemical structure and molecular mass.

In a targeted metabolomics approach, we recently identified significant intra-individual metabolic alterations in patients with PPGL before vs. after tumor removal and demonstrated that several of those are related to cardiovascular risk (15). Characterization of the metabolic profile in patients with PPGL might help to understand the metabolic effects of excessive catecholamine levels and harbor additional diagnostic potential.

The aim of our study was to characterize differences in plasma metabolic profile between patients with PPGL and controls with consideration of the secretory phenotype. We applied tandem mass spectrometry using a targeted metabolomics approach and logistic regression modeling to identify discriminative pattern potentially useful for diagnostic workup of PPGL.

Subjects and Methods

Subjects

Patients with suspected PPGL were recruited from a single center participating in the Prospective Monoamine-Producing tumor (PMT) study, which has been described in detail previously (2). The diagnosis of PPGL was based on biochemical assessment, imaging, and histology. Follow-up ruled out the presence of PPGLs in patients who served as controls. The latter were matched for sex and age at the date of sampling according to patient data. The study protocol was approved by the Ethics Committee of the University Hospital Würzburg (104/11). All patients provided written informed consent.

Sample Collection

Plasma samples were collected as described elsewhere (2). Briefly, blood was drawn in the morning after an overnight fast for at least 8h and in a supine position for at least 30 minutes. Patients were instructed to refrain from alcohol, nicotine, decaffeinated and caffeinated beverages for 12 hours as well as avoid acetaminophen five days before sample collection (16). Blood was collected into EDTA or heparinized tubes and placed on ice before centrifugation at 20°C for five minutes at 4000 rpm. Plasma was aliquoted and the samples were stored at -80°C until assayed. Urine collection was performed according to the PMT protocol.

Mass Spectrometry

Plasma free metanephrines and urine catecholamines were measured as previously described (17–19).

Targeted metabolomics was performed by using the AbsoluteIDQ™-p180 Kit (Biocrates Life Sciences AG, Innsbruck, Austria). The method has been described in detail previously (15, 20, 21) and complies with EMA “Guideline on bioanalytical method validation” (July 21st 2011).The measurement consists of a ultra-high performance liquid chromatography (UHPLC) separation step and a flow injection analysis (FIA) step, both followed by mass spectrometry analyses (LC-MS/MS and FIA-MS/MS). This method enables for measurement of a total of 188 metabolites, of which 42 are included in the LC-MS/MS part (21 amino acids, 21 biogenic amines) and 146 metabolites in the FIA-MS/MS protocol (40 acylcarnitines including free carnitine, 38 phosphatidylcholines with acyl/acyl side chains [PCaa], 38 phosphatidylcholines with acyl/alkyl side chains [PCae], 14 lysophosphatidylcholines [lysoPC], 15 sphingolipids [SM] and the sum of hexoses [H1]).

A volume of 10 µl plasma was used and prepared according to the manufacturer’s manual. Internal standards served as reference for quantification, human reference plasma was included into each batch to ensure quality control, comparability between batch measurements, and normalization of the data (20). Metabolite concentrations are given in µmol/l. LC-MS/MS and FIA-MS/MS were performed by using SCIEX QTRAP^® 4500MD MS-system (SCIEX, Darmstadt, Germany) coupled to an Agilent 1290 Infinity UHPLC-system (Agilent, Santa Clara, USA). Analyst^® software version 1.6.3 MD (SCIEX, Darmstadt, Germany) was used for data procession. Data was validated and processed with MetIDQ™ software version 5.5.4-DB100 Boron-2623 (Biocrates Life Sciences AG, Innsbruck, Austria).

Genetics

Genetic data were retrieved from patient records or provided by the CNIO institute in Madrid as a part of the PMT study. Targeted next generation sequencing assay, Sanger sequencing and multiplex ligation-dependent probe amplification or custom array comparative genomic hybridization for deletion detection (22, 23) were applied as appropriate.

Statistical Analysis

Baseline data are shown as frequencies for categorical variables and as medians with interquartile range (IQR) for numerical variables. Malignancy was defined as the presence of metastases in non-chromaffin organs. The secretory phenotype of PPGL was characterized as noradrenergic and adrenergic according to an established algorithm which has been descripted in detail elsewhere (24). Metabolites with more than 40% of concentrations below the lower limit of quantification (LLOQ) and samples with more than 40% of analytes lower than the lower limit of detection (LOD) were excluded from further analysis. In the remaining metabolites, the non-valid values measured below LLOQ and LOD were left unchanged and included in further analyses. Values with no detectable signal were replaced by (LOD/ √2) x (random number between 0.75-1.25) (25). To detect metabolic alterations associated with catecholamine excess, metabolite values in PPGL and controls were compared. Subgroups were analyzed after stratification for sex (males vs. females), BMI (≤ 25 kg/m² vs > 25 kg/m²) and secretory phenotype (adrenergic vs. noradrenergic). Comparisons between groups were analyzed using the Mann-Whitney-U test, significance was defined as p-value <0.05. The calculation of false discovery rate (FDR) corrected p-values was performed according to the method of Benjamini and Hochberg (26). Spearman test was used for correlations between metabolite and catecholamine concentrations. Statistical analyses were performed by SPSS version 25 (IMP, New York, USA) and Prism 7.05 (GraphPad, San Diego, CA, USA), for principal component analysis, MetaboAnalyst (4.0) was used.

Logistic Regression Modelling

Feature selection models were developed by applying the machine learning methods Elastic net (ELA), Support Vector Machine (SVM) and Gradient Boosting Machine (GBM) using the caret package version 6.0.84. The analysis was done in R (version 3.5.3; script is given online). The data set was normalized using the PreProcess function of the caret package (version 6.0.84). We split the dataset in a training (80%) and validation/test (20%) dataset. We tested each model using repeated 10-fold cross-validation. The variables were selected using the impact that they had on the predictive power of the different models. The models were compared using the predictive values accuracy (correct classification) and kappa (inter-rater reliability; classification including random chance normalization). Identified variables were further analysed using a Wilcoxon test to determine if there is a systematic difference between the conditions (class: PPGL vs. control).

Results

Patients Characteristics

The study workflow is depicted in Figure 1. 36 patients with confirmed PPGL prior any specific treatment and controls matched for sex and age at date of sample were selected (Table 1). In controls, PPGL was suspected based on the incidental finding of an adrenal mass upon imaging for an unrelated condition (n=21), signs and symptoms suggestive for PPGL (n=9) or therapy resistant hypertension (n=6) but excluded by normal follow-up biochemistry, negative imaging, resolved signs and symptoms or an alternative diagnosis (17). Other endocrinological causes for resistant hypertension or adrenal incidentalomas such as hyperaldosteronism, acromegaly, hyperthyroidism, hyperparathyroidism, or Cushing’s syndrome were excluded in all subjects.

FIGURE 1

Figure 1 Summary of the workflow leading to identification of significant metabolites via liquid chromatography tandem mass spectrometry (LC-MS/MS), in pheochromocytoma and paraganglioma (PPGL),The bioinformatic approach included Elastic net (ELA), Gradient Boosting Machine (GBM), Support Vector Machine (SVM) and Principal component analysis (PCA) approaches.

TABLE 1

Table 1 Patient characteristics stratified by patients and controls.

Complete urinary catecholamine data (31/36 controls, 27/36 PPGL) and genetic data (33/36 PPGL) were available in a subset of individuals. There were no statistically significant differences between groups concerning age and time between sampling and measurement of metabolomics while a statistically significant difference in Body mass index (in kg/m²) was present in PPGL (25.2 [23.6–26.7]) vs controls (28.6 [26.0–31.2], p=0.043). Plasma markers of catecholamine excess were significantly increased in PPGL.

Targeted Metabolomics PPGL vs. Controls

Overall, 130 of 188 measured metabolites were included in the statistical analysis (Supplemental Data). However, only when p-values were not corrected for FDR, four of them showed significantly different concentrations between the two groups. In PPGL vs. controls (Figure 2) the amino acids histidine (75.40 [61.03-87.05] vs. 86.40 [75.63-96.35] µmol/l, p=0.004) and threonine (105.00 [88.57-125.00] vs. 128.00 [93.32-147.50] µmol/l, p=0.008) were significantly lower, while lyso PC a C28:0 (0.11 [0.10-0.12] vs. 0.12 [0.11-0.14] µmol/l, p=0.044) was only slightly decreased. On the opposite, the sum of hexoses was significantly higher in PPGL patients compared to controls (4844.00 [4325.50-5364.50] vs. 4215.50 [3791.00-5086.00] µmol/l, p=0.018). The plasma concentrations of biogenic amines, acylcarnitines, and sphingolipids were comparable between PPGL and controls.

FIGURE 2

Figure 2 Scatter plot of median plasma levels from metabolites measured by LC-MS/MS with significant differences between PPGL patients and controls. (Mann-Whitney-U-test, p < 0.05). H1, sum of hexoses; lysoPC, lysophosphatidylcholine; PPGL, pheochromocytoma and paraganglioma. *<0.05, **<0.01.

Correlation With Plasma Metanephrine and Urine Catecholamine Values

The association of metabolic changes with catecholamine excretion has already been demonstrated (15). Therefore, we correlated altered metabolites with urinary catecholamines, which represent the biologically active form, and MN, NMN and MTY in plasma as non-functional disease markers (Table 2 and Figure 3). Histidine showed a significant negative correlation with plasma NMN and plasma MTY, as well as with urine free epinephrine (EPI) and urine free dopamine (DA). Threonine was negatively correlated with plasma MTY, urine free norepinephrine (NE), and urine free EPI. LysoPC a C28:0 was negatively associated with urine free DA, whereas the sum of hexoses showed a positive correlation with all plasma metanephrines and with urine free NE. Furthermore, Histidine, threonine and lysoPC C:28 revealed a negative correlation with the total urinary catecholamines.

TABLE 2

Table 2 Plasma levels of significant altered metabolites (p ≤ 0.05) in patients with PPGL in comparison to controls including subgroup analysis and the correlation with free plasma metanephrines and 24h urinary free catecholamine excretion values.

FIGURE 3

Figure 3 Correlation between plasma concentrations of metanephrine, normetanephrine, and methoxytyramine, 24h urine concentrations of catecholamines, and plasma concentrations of significantly altered metabolites in PPGL patients. Spearmen-coefficient r_s is presented by color coding (positive correlation: red; negative correlation: blue). An Asterisk indicates a statistically significant correlation at the level of (p < 0.05). DA, dopamine; EPI, epinephrine; MN, metanephrine; MTY, methoxytyramine; NE, norepinephrine; TotalUrineCat, Total urine catecholamines.

Subgroup Analyses

Despite the small group sizes, we also explored differences in subgroups. However, these were significant only without FDR correction. If only males were taken into account (n=19), patients with PPGL had lower levels of threonine (102.00 [85.50-127.00] vs. 132.00 [104.00-156.00] µmol/l, p=0.008) and higher levels of H1 (4845.00 [4489.00-5285.00] vs. 4079.00 [3789.00-5131.00] µmol/l, p=0.050) than controls. In addition, alterations in three additional metabolites were present: Lyso PC a C16:1 showed lower levels in PPGL (1.68 [1.09-1.79] vs. 1.89 [1.54-2.16] µmol/l, p=0.025), whilst PC ae C30:2 and SM (OH) C14:1 had significantly higher levels compared to controls (0.07 [0.06-0.08] vs. 0.06 [0.06-0.07] µmol/l, p=0.040 and 4.77 [3.20-5.43] vs. 3.46 [2.59-4.38] µmol/l, p=0.040).

In the subgroup with BMI below or equal 25 kg/m² the PPGL patients (n=22) exhibited lower level of histidine (75.20 [64.35-85.90] vs. 87.50 [82.40-103.00] μmol/l, p=0.006) and higher level of H1 (4773.00 [4357.00-5311.50] vs. 4146.00 [3789.00-4428.00] µmol/l, p=0.003) than controls (n=15). In the subgroup with BMI above 25 kg/m² PPGL patients (n=12) had higher level of SM OH C22:1 (7.09 [5.49-8.29] vs. 5.88 [4.81-6.64] µmol/l, p=0.048) and SM OH C22:2 (5.66 [5.16-6.81] vs. 4.90 [3.94-5.51] µmol/l, p=0.036). Octadecenoylcarnitine (0.099 [0.088-0.131] vs. 0.135 [0.102-0.203] µmol/l, p=0.048), octadecadienylcarnitine 0.029 [0.023-0.034] vs. 0.035 [0.029-0.064] µmol/l, p=0.044) and histidine 71.10 [58.80-85.17] vs. 84.20 [71.50-93.80] µmol/l, p=0.036), were lower in PPGL patients than in controls (n=20), The same applies for ornithine (60.75 [54.10-71.30] vs. 92.60 [58.10-124.50] µmol/l, p=0.036) and threonine (98.85 [85.50-122.25] vs. 127.00 [93.33-153.50] µmol/l, p=0.036).

Female PPGL patients exhibited higher levels of lysoPC a C20:4 (3.95 [3.50-4.95] vs. 2.62 [2.44-3.54] µmol/l, p=0.006), PC aa C36:4 (203.00 [165.00-249.50] vs. 160.00 [124.00-211.00] µmol/l, p=0.041) and PC aa C38:4 (94.10 [77.35-107.00] vs. 82.60 [64.40-91.55] µmol/l, p=0.049) as well as lower values of PC ae C38:1 (0.37 [0.28-1.24] vs. 0.71 [0.47-1.69] µmol/l, p=0.022) and histidine (64.10 [56.70-78.30] vs. 85.90 [71.70-104.50] µmol/l, p=0.001).

Analyzing only the subgroup of adrenergic phenotypes, we found increased levels of glycin (221.00 [182.00-373.00] vs. 167.00 [145.00-192.00] µmol/l, p=0.007) and lysoPC a C20:4 (4.00 [3.37-5.18] vs. 3.10 [2.46-4.04] µmol/l, p=0.019] combined with decreased levels of lysoPC a C28:0 (0.11 [0.08-0.11] vs. 0.13 [0.11-0.14] µmol/l, p=0.026) in patients with PPGL compared to controls. Histidine, which had lower levels in the entire group as well as in females, showed also lower concentrations in adrenergic phenotypes (69.00 [58.40-76.10] vs. 85.10 [71.30-93.60] µmol/l, p=0.004). In PPGL with noradrenergic phenotype we found lower concentrations of C0 (35.60 [28.70-41.15] vs. 43.00 [34.05-48.70] µmol/l, p=0.042), asparagine (37.20 [32.00-41.95] vs. 42.00 [38.65-48.05] µmol/l, p=0.013), threonine (112.00 [96.20-127.50] vs. 133.00 [120.00-151.00] µmol/l, p=0.002) and ADMA (0.51 [0.39-0.80] vs. 0.62 [0.55-0.91] µmol/l, p=0.048).

Feature Selection Using Machine Learning Techniques and Principal Component Analysis

The machine learning models GBM, ELA, and SVM were run to determine features which are important to class prediction. Each of the models was run with 10-fold cross-validation on the training dataset, the features which contributed most to class prediction were obtained and the models were compared based on their estimated performance. The GBM, ELA and SVM had an estimated accuracy of 0.67, 0.53 and 0.53 respectively (Kappa: GBM 0.33, ELA 0.06 and SVM 0.08; Figure 4A). ELA and SVM selected 20 variables, whereas the GBM only selected 9 variables (Figure 4B). The GBM had the best estimated predictive value of which hexose showed the largest contribution to diagnosis, calculated using an out-of-bag estimate of the improvement in predictive performance. Comparing the selected variables, only H1 was shared between the 3 modelling algorithms. Furthermore, for GBM each of the selected variables was evaluated for difference (Supplemental Data), which found that only the H1 predictor has a significant difference between the two classes. Testing for the class prediction in the validation dataset showed limited predictive value for all models.

FIGURE 4

Figure 4 Variable selection in the training dataset. (A) Accuracy and Kappa comparisons of the three models used in feature selection. (B) Importance of the variable selection and out-of-bag predictive performance estimation.

Similar analyses were also preformed to select predictors of catecholamine producing tumor types (pheo vs. PGL vs. control) phenotypes (adrenergic vs. noradrenergic vs. control) and malignancy (benign vs. malignant vs. control). The predictors that were selected in all these cases were the same as from the initial analysis, but further investigation of the distribution of the predictors showed a much weaker difference between the groups (diagnosis, phenotype and malignancy) when compared with the distribution between the difference classes (PPGL vs. control). GBM selected the same variables in each section (Supplemental Data). However, accuracy was always lower than 60%. The most commonly selected variable H1 showed difference between different analyses in phenotype p=0.10 between control and adrenergic/noradrenergic, but p=0.68 between the two states. In malignancy we calculated p=0.11 between benign and control and p=0.08 between benign and malignant, but p=0.76 between control and malignant.

The score plots obtained from principal component analysis models after logarithmic normalization of the entire dataset and subgroup (adrenergic, noradrenergic) demonstrated that the groups were not well discriminated (Supplemental Data).

Discussion

In our study 36 PPGL patients and matched controls were analyzed with targeted metabolomics. Despite the highly standardized sampling conditions and quantitative liquid chromatography-tandem mass spectrometry, we only found metabolites with significant differences between PPGL patients and matched controls when no correction for FDR was performed. After correction for multiple comparisons, the statistical significance was not retained for all metabolites. Classifying substances as significant without FDR-correction was reasonable in our preliminary setting as the focus was to identify potentially relevant metabolites, which have to be validated in further studies. In this approach it may be preferable to explore leads that may turn out to be wrong than loosing promising markers in early state by stringent statistical criteria, as has been argued by others (27). By using machine learning, we failed to establish metabolic signatures associated with PPGL diagnosis, indicating little value of such targeted metabolomics approach for diagnostic purposes in PPGL at variance to plasma metanephrin, normetanephrin and 3-methoxytyramin which have proven excellent sensitivity and specificity when performed with appropriate preanalytics, analytics and reference intervals.

It is well known that catecholamine excess leads to a diabetogenic state (15, 28–30) and we accordingly found higher levels of hexoses in PPGL patients. This reflects increased glycogen catabolism, glucagon release and gluconeogenesis finally leading to the high prevalence of diabetes mellitus (21-37%) in PPGL patients.

Erlic et al. compared the same metabolite spectrum in PPGL patients before and post-surgery and found lower histidine levels in preoperative samples (15). A low histidine level is linked to type 2 diabetes, increased inflammation and cardiovascular disease (31–33), potentially explaining (at least in part) such effects in PPGL patients (31–34).

The negative correlation of histidine with plasma NMN and MTY, urine free EPI and DA in our cohort suggests a catabolic phenotype which has been linked to proinflammatory mediators as well (35). In this regard it is noteworthy that BMI was significantly lower in patients compared to controls. While one may argue that this reflects an imbalance in the matching of base line characteristics, it may rather reflect the catabolic effect of PPGL.

Sex specific differences in metabolic pattern found by Erlic et al. by comparing intra-individual metabolite profiles prior and after surgical tumor removal were confirmed here even if the results should be treated with caution due to the small number of patients in the subgroups (15).

Significant changes between PPGL and controls observed in threonine and histidine are in accordance with other studies that focused on cancer (36, 37). For example, Miyagi et al. also showed decreased levels of histidine in patients with gastric, colorectal, lung and breast cancer, while threonine was lower in gastric as well as colorectal cancer and higher in bronchial carcinoma. The mechanistic background and clinical significance are still subject of discussion (37).

ML approaches imply the need for a substantial amount of data for the development of clinical useful diagnostic models, thus limits the application potential in rare diseases such as PPGL. We overcome this by combining statistical correlation analysis with 3 different feature selection algorithm and applied several validation steps such as data splitting and 10-fold cross-validation on a non-linear and large dimensional variable (130 metabolites) dataset (38). However, even the best performing model (GBM with 9 selected metabolites) showed low predictive values within the test dataset. The GBM model shows a higher sensitivity (87.5%) compared to the other models but was outperformed by ELA and SVM in terms of specificity. Given the low accuracy and no significant p-value, the models show low ability to distinguish between PPGL and control. Similarly, model selection based on diagnosis, phenotype and malignancy predictors showed a much weaker difference between the groups (accuracy always <60%). Given the small size of the disease subsets and the poor significance found in the distributions between each group subset, the confidence with which the features can be used to differentiate between the different group subsets was limited.

It might be argued that a larger study has enabled us to identify a distinct phenotype using ML-based selection approaches which then could possibly include a larger number of features under study. We do not share this point of view because any test applied in PPGL for diagnostic purposes requires a much higher sensitivity and specificity compared to that found in our pilot study to be clinically meaningful. This is particularly true when the generally low pre-test likelihood of a PPGL and the prevalence of diabetes and catabolism in an unselected population is considered. Consequently, a further analysis of the total cohort of the PMT study regarding this pilot study does not seem reasonable.

PPGL-associated mutations in genes involved in the Krebs cycle and electron transport chain have shown to translate into characteristic tumoral metabolic changes (39). Tumors caused by mutations in the SDH genes that belong to this cluster 1 tumors exhibit an increased succinate fumarate ratio in tumor tissue (23, 40, 41) which leads to a pseudo-hypoxic phenotype that downstream activates hypoxia induced factor signaling and angiogenesis.

Recently, Wallace et al. showed that assessment of Krebs cycle related metabolites by using LC-MS/MS in addition to immunohistochemistry improved the diagnosis of SDH impairment at a functional level (42). These tumoral metabolic pathway alterations translate into characteristic secretory pattern of catecholamine metabolites which contribute to the diagnosis of malignancy (39, 43–46). Of note, they also appear to have therapeutic relevance (39, 43–45).

Taken together, we confirmed previous findings of metabolic alterations caused by PPGL related catecholamine excess by comparing PPGL patients and controls. We applied machine learning algorithms, but these failed to provide feature-selection signatures that may be useful for PPGL diagnosis in clinical routine. Still, our study broadens and complements the understanding of changes in the metabolic profile of patients with PPGL.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by Ethics committee of the Medical Faculty, University of Würzburg. Written informed consent to participate in this study was obtained from all study participants and, in the case of minors, from their legal representatives.

Author Contributions

JM: Data acquisition, data analysis, data interpretation, and writing of first paper draft. MK: Conception of the work, data acquisition, data analysis, writing of first paper draft, and revising the work critically for important intellectual content. OR-L: Data analysis, data interpretation, writing of first paper draft, and revising the work. TD: Data interpretation and revising the work critically for important intellectual content. MP: Data acquisition, data analysis, data interpretation, and revising the work critically for important intellectual content. CP: Data acquisition and revising the work critically for important intellectual content. DW: Revising the work critically for important intellectual content. MR: Data acquisition and revising the work critically for important intellectual content. JA: Data acquisition and revising the work critically for important intellectual content. MF: Planning and design of the work and revising the work critically for important intellectual content. MKu: Data analysis, data interpretation, and drafting and revising it critically for important intellectual content. MKr: Conception, planning and design of the work, data analysis, data interpretation, drafting and revising it critically for important intellectual content, important intellectual content. MKu and MKr contributed equally to this study. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the German Research Council (DFG) project 314061271 (SFB/CRC Transregio 205: The adrenal – central relay in health and disease).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We thank Sabine Kendl for expert technical assistance. Funding by the German Federal Ministry of Education and Research (BMBF) is acknowledged (CompLS program grant 031L0262C and Era-Net grant 01KT1801 to MKu). OR-L acknowledges the Schickedanz Kinderkrebsstiftung.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2021.722656/full#supplementary-material

Abbreviations

a, acyl; aa, diacyl; ae, acyl-alkyl; C x:y, indicates the lipid chain composition where “x” is the number of carbons and “y” the number of bonds; C, carbon; DA, dopamine; ELA, Elastic net; EPI, epinephrine; FDR, false discovery rate; FIA, flow injection analysis; GBM, Gradient Boosting Machine; H1, sum of hexoses; IQR, interquartile range; LC-MS/MS, liquid chromatography tandem mass spectrometry; LLOQ, lower limit of quantification; LOD, lower limit of detection; lysoPC, lysophosphatidylcholine; MN, metanephrine; MTY, methoxytyramine; NE, norepinephrine; NMN, normetanephrine; PC, phosphatidylcholine; PMT, Prospective Monoamine-producing Tumor Study; PPGL, pheochromocytoma/paraganglioma; SDH, succinate dehydrogenase; SM, sphingolipids; SVM, Support Vector Machine; UHLPC, ultra-high performance liquid chromatography.

References

1. Lenders JW, Eisenhofer G, Mannelli M, Pacak K. Phaeochromocytoma. Lancet (2005) 366(9486):665–75. doi: 10.1016/S0140-6736(05)67139-5

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Eisenhofer G, Prejbisz A, Peitzsch M, Pamporaki C, Masjkur J, Rogowski-Lehmann N, et al. Biochemical Diagnosis of Chromaffin Cell Tumors in Patients at High and Low Risk of Disease: Plasma Versus Urinary Free or Deconjugated O-Methylated Catecholamine Metabolites. Clin Chem (2018) 64(11):1646–56. doi: 10.1373/clinchem.2018.291369

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Whalen RK, Althausen AF, Daniels GH. Extra-Adrenal Pheochromocytoma. J Urol (1992) 147(1):1–10. doi: 10.1016/S0022-5347(17)37119-7

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Lenders JWM, Eisenhofer G. Update on Modern Management of Pheochromocytoma and Paraganglioma. Endocrinol Metab (Seoul) (2017) 32(2):152–61. doi: 10.3803/EnM.2017.32.2.152

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Lenders JWM, Kerstens MN, Amar L, Prejbisz A, Robledo M, Taieb D, et al. Genetics, Diagnosis, Management and Future Directions of Research of Phaeochromocytoma and Paraganglioma: A Position Statement and Consensus of the Working Group on Endocrine Hypertension of the European Society of Hypertension. J Hypertens (2020) 38(8):1443–56. doi: 10.1097/HJH.0000000000002438

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Favier J, Amar L, Gimenez-Roqueplo AP. Paraganglioma and Phaeochromocytoma: From Genetics to Personalized Medicine. Nat Rev Endocrinol (2015) 11(2):101–11. doi: 10.1038/nrendo.2014.188

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Dahia PL. Pheochromocytoma and Paraganglioma Pathogenesis: Learning From Genetic Heterogeneity. Nat Rev Cancer (2014) 14(2):108–19. doi: 10.1038/nrc3648

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Plouin PF, Amar L, Dekkers OM, Fassnacht M, Gimenez-Roqueplo AP, Lenders JW, et al. European Society of Endocrinology Clinical Practice Guideline for Long-Term Follow-Up of Patients Operated on for a Phaeochromocytoma or a Paraganglioma. Eur J Endocrinol (2016) 174(5):G1–G10. doi: 10.1530/EJE-16-0033

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Geroula A, Deutschbein T, Langton K, Masjkur J, Pamporaki C, Peitzsch M, et al. Pheochromocytoma and Paraganglioma: Clinical Feature-Based Disease Probability in Relation to Catecholamine Biochemistry and Reason for Disease Suspicion. Eur J Endocrinol (2019) 181(4):409–20. doi: 10.1530/EJE-19-0159

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Rao D, Peitzsch M, Prejbisz A, Hanus K, Fassnacht M, Beuschlein F, et al. Plasma Methoxytyramine: Clinical Utility With Metanephrines for Diagnosis of Pheochromocytoma and Paraganglioma. Eur J Endocrinol (2017) 177(2):103–13. doi: 10.1530/EJE-17-0077

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Eisenhofer G. Screening for Pheochromocytomas and Paragangliomas. Curr Hypertens Rep (2012) 14(2):130–7. doi: 10.1007/s11906-012-0246-y

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Weismann D, Peitzsch M, Raida A, Prejbisz A, Gosk M, Riester A, et al. Measurements of Plasma Metanephrines by Immunoassay vs Liquid Chromatography With Tandem Mass Spectrometry for Diagnosis of Pheochromocytoma. Eur J Endocrinol (2015) 172(3):251–60. doi: 10.1530/EJE-14-0730

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Murakami M, Sun N, Greunke C, Feuchtinger A, Kircher S, Deutschbein T, et al. Mass Spectrometry Imaging Identifies Metabolic Patterns Associated With Malignant Potential in Pheochromocytoma and Paraganglioma. Eur J Endocrinol (2021) 185(1):179–91. doi: 10.1530/EJE-20-1407

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Patti GJ, Yanes O, Siuzdak G. Innovation: Metabolomics: The Apogee of the Omics Trilogy. Nat Rev Mol Cell Biol (2012) 13(4):263–9. doi: 10.1038/nrm3314

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Erlic Z, Kurlbaum M, Deutschbein T, Nolting S, Prejbisz A, Timmers H, et al. Metabolic Impact of Pheochromocytoma/Paraganglioma: Targeted Metabolomics in Patients Before and After Tumor Removal. Eur J Endocrinol (2019) 181(6):647–57. doi: 10.1530/EJE-19-0589

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Galetta F, Franzoni F, Bernini G, Poupak F, Carpi A, Cini G, et al. Cardiovascular Complications in Patients With Pheochromocytoma: A Mini-Review. BioMed Pharmacother (2010) 64(7):505–9. doi: 10.1016/j.biopha.2009.09.014

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Peitzsch M, Prejbisz A, Kroiss M, Beuschlein F, Arlt W, Januszewicz A, et al. Analysis of Plasma 3-Methoxytyramine, Normetanephrine and Metanephrine by Ultraperformance Liquid Chromatography-Tandem Mass Spectrometry: Utility for Diagnosis of Dopamine-Producing Metastatic Phaeochromocytoma. Ann Clin Biochem (2013) 50(Pt 2):147–55. doi: 10.1258/acb.2012.012112

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Peitzsch M, Pelzel D, Glockner S, Prejbisz A, Fassnacht M, Beuschlein F, et al. Simultaneous Liquid Chromatography Tandem Mass Spectrometric Determination of Urinary Free Metanephrines and Catecholamines, With Comparisons of Free and Deconjugated Metabolites. Clin Chim Acta (2013) 418:50–8. doi: 10.1016/j.cca.2012.12.031

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Peitzsch M, Adaway JE, Eisenhofer G. Interference From 3-O-Methyldopa With Ultra-High Performance LC-MS/MS Measurements of Plasma Metanephrines: Chromatographic Separation Remains Important. Clin Chem (2015) 61(7):993–6. doi: 10.1373/clinchem.2015.239962

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Zukunft S, Sorgenfrei M, Prehn C, Moller G, Adamski J. Targeted Metabolomics of Dried Blood Spot Extracts. Chromatographia (2013) 76(19-20):1295–305. doi: 10.1007/s10337-013-2429-3

CrossRef Full Text | Google Scholar

21. Siskos AP, Jain P, Romisch-Margl W, Bennett M, Achaintre D, Asad Y, et al. Interlaboratory Reproducibility of a Targeted Metabolomics Platform for Analysis of Human Serum and Plasma. Anal Chem (2017) 89(1):656–65. doi: 10.1021/acs.analchem.6b02930

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Curras-Freixes M, Pineiro-Yanez E, Montero-Conde C, Apellaniz-Ruiz M, Calsina B, Mancikova V, et al. Pheoseq: A Targeted Next-Generation Sequencing Assay for Pheochromocytoma and Paraganglioma Diagnostics. J Mol Diagn (2017) 19(4):575–88. doi: 10.1016/j.jmoldx.2017.04.009

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Richter S, Gieldon L, Pang Y, Peitzsch M, Huynh T, Leton R, et al. Metabolome-Guided Genomics to Identify Pathogenic Variants in Isocitrate Dehydrogenase, Fumarate Hydratase, and Succinate Dehydrogenase Genes in Pheochromocytoma and Paraganglioma. Genet Med (2019) 21(3):705–17. doi: 10.1038/s41436-018-0106-5

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Pamporaki C, Hamplova B, Peitzsch M, Prejbisz A, Beuschlein F, Timmers H, et al. Characteristics of Pediatric vs Adult Pheochromocytomas and Paragangliomas. J Clin Endocrinol Metab (2017) 102(4):1122–32. doi: 10.1210/jc.2016-3829

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Xia J, Psychogios N, Young N, Wishart DS. Metaboanalyst: A Web Server for Metabolomic Data Analysis and Interpretation. Nucleic Acids Res (2009) 37(Web Server issue):W652–60. doi: 10.1093/nar/gkp356

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I. Controlling the False Discovery Rate in Behavior Genetics Research. Behav Brain Res (2001) 125(1-2):279–84. doi: 10.1016/S0166-4328(01)00297-2

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Rothman KJ. No Adjustments Are Needed for Multiple Comparisons. Epidemiology (1990) 1(1):43–6. doi: 10.1097/00001648-199001000-00010

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Leroy ED, James JR, Semans H, Howard JE. Adrenal Medullary Tumor (Pheochromocytoma) and Diabetes Mellitus; Disappearance of Diabetes After Removal of the Tumor. Ann Internal Med (1944) 20(5):815–21. doi: 10.003-4819-20-5-815

Google Scholar

29. McCullagh EP, Engel WJ. Pheochromocytoma With Hypermetabolism: Report of Two Cases. Ann Surg (1942) 116(1):61–75. doi: 10.1097/00000658-194207000-00008

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Mesmar B, Poola-Kella S, Malek R. The Physiology Behind Diabetes Mellitus in Patients With Pheochromocytoma: A Review of the Literature. Endocrine Pract (2017) 23(8):999–1005. doi: 10.4158/EP171914.RA

CrossRef Full Text | Google Scholar

31. Du SS, Sun SH, Liu L, Zhang Q, Guo FC, Li CL, et al. Effects of Histidine Supplementation on Global Serum and Urine H-1 NMR-Based Metabolomics and Serum Amino Acid Profiles in Obese Women From a Randomized Controlled Study. J Proteome Res (2017) 16(6):2221–30. doi: 10.1021/acs.jproteome.7b00030

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Mong MC, Chao CY, Yin MC. Histidine and Carnosine Alleviated Hepatic Steatosis in Mice Consumed High Saturated Fat Diet. Eur J Pharmacol (2011) 653(1-3):82–8. doi: 10.1016/j.ejphar.2010.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Hasegawa S, Ichiyama T, Sonaka I, Ohsaki A, Okada S, Wakiguchi H, et al. Cysteine, Histidine and Glycine Exhibit Anti-Inflammatory Effects in Human Coronary Arterial Endothelial Cells. Clin Exp Immunol (2012) 167(2):269–74. doi: 10.1111/j.1365-2249.2011.04519.x

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Bosanska L, Petrak O, Zelinka T, Mraz M, Widimsky J, Haluzik M. The Effect of Pheochromocytoma Treatment on Subclinical Inflammation and Endocrine Function of Adipose Tissue. Physiol Res (2009) 58(3):319–25. doi: 10.33549/physiolres.931483

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Torres KCL, Antonelli LRV, Souza ALS, Teixeira MM, Dutra WO, Gollob KJ. Norepinephrine, Dopamine and Dexamethasone Modulate Discrete Leukocyte Subpopulations and Cytokine Profiles From Human PBMC. J Neuroimmunol (2005) 166(1-2):144–57. doi: 10.1016/j.jneuroim.2005.06.006

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Kobayashi T, Nishiumi S, Ikeda A, Yoshie T, Sakai A, Matsubara A, et al. A Novel Serum Metabolomics-Based Diagnostic Approach to Pancreatic Cancer. Cancer Epidemiol Biomarkers Prev (2013) 22(4):571–9. doi: 10.1158/1055-9965.EPI-12-1033

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Miyagi Y, Higashiyama M, Gochi A, Akaike M, Ishikawa T, Miura T, et al. Plasma Free Amino Acid Profiling of Five Types of Cancer Patients and Its Application for Early Detection. PloS One (2011) 6(9). doi: 10.1371/journal.pone.0024143

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Vey J, Kapsner LA, Fuchs M, Unberath P, Veronesi G, Kunz M. A Toolbox for Functional Analysis and the Systematic Identification of Diagnostic and Prognostic Gene Expression Signatures Combining Meta-Analysis and Machine Learning. Cancers (Basel) (2019) 11(10). doi: 10.3390/cancers11101606

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Buffet A, Burnichon N, Favier J, Gimenez-Roqueplo AP. An Overview of 20 Years of Genetic Studies in Pheochromocytoma and Paraganglioma. Best Pract Res Clin Endocrinol Metab (2020) 34(2). doi: 10.1016/j.beem.2020.101416

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Lendvai N, Pawlosky R, Bullova P, Eisenhofer G, Patocs A, Veech RL, et al. Succinate-to-Fumarate Ratio as a New Metabolic Marker to Detect the Presence of SDHB/D-Related Paraganglioma: Initial Experimental and Ex Vivo Findings. Endocrinology (2014) 155(1):27–32. doi: 10.1210/en.2013-1549

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Richter S, Peitzsch M, Rapizzi E, Lenders JW, Qin N, de Cubas AA, et al. Krebs Cycle Metabolite Profiling for Identification and Stratification of Pheochromocytomas/Paragangliomas Due to Succinate Dehydrogenase Deficiency. J Clin Endocrinol Metab (2014) 99(10):3903–11. doi: 10.1210/jc.2014-2151

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Wallace PW, Conrad C, Bruckmann S, Pang Y, Caleiras E, Murakami M, et al. Metabolomics, Machine Learning and Immunohistochemistry to Predict Succinate Dehydrogenase Mutational Status in Phaeochromocytomas and Paragangliomas. J Pathol (2020) 251(4):378–87. doi: 10.1002/path.5472

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Tevosian SG, Ghayee HK. Pheochromocytomas and Paragangliomas. Endocrinol Metab Clinics North America (2019) 48(4):727–50. doi: 10.1016/j.ecl.2019.08.006

CrossRef Full Text | Google Scholar

44. Matlac DM, Vanova KH, Bechmann N, Richter S, Folberth J, Ghayee HK, et al. Succinate Mediates Tumorigenic Effects via Succinate Receptor 1: Potential for New Targeted Treatment Strategies in Succinate Dehydrogenase Deficient Paragangliomas. Front Endocrinol (2021) 12:129. doi: 10.3389/fendo.2021.589451

CrossRef Full Text | Google Scholar

45. Fankhauser M, Bechmann MN, Lauseker M, Goncalves J, Favier J, Klink B, et al. Synergistic Highly Potent Targeted Drug Combinations in Different Pheochromocytoma Models Including Human Tumor Cultures. Endocrinology (2019) 160(11):2600–17. doi: 10.1210/en.2019-00410

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Eisenhofer G, Klink B, Richter S, Lenders JW, Robledo M. Metabologenomics of Phaeochromocytoma and Paraganglioma: An Integrated Approach for Personalised Biochemical and Genetic Testing. Clin Biochem Rev (2017) 38(2):69–100.

PubMed Abstract | Google Scholar