Serum Steroid Ratio Profiles in Prostate Cancer: A New Diagnostic Tool Toward a Personalized Medicine Approach

Background: Serum steroids are crucial molecules altered in prostate cancer (PCa). Mass spectrometry (MS) is currently the elected technology for the analysis of steroids in diverse biological samples. Steroids have complex biological pathways and stoichiometry and it is important to evaluate their quantitative ratio. MS applications to patient hormone profiling could lead to a diagnostic approach.

Methods: Here, we employed the Surface Activated Chemical Ionization-Electrospray-NIST (SANIST) developed in our laboratories, to obtain quantitative serum steroid ratio relationship profiles with a machine learning Bayesian model to discriminate patients with PCa. The approach is focused on steroid relationship profiles and disease association.

Results: A pilot study on patients affected by PCa, benign prostate hypertrophy (BPH), and control subjects [prostate-specific antigen (PSA) lower than 2.5 ng/mL] was done in order to investigate the classification performance of the SANIST platform. The steroid profiles of 71 serum samples (31 controls, 20 patients with PCa and 20 subjects with benign prostate hyperplasia) were evaluated. The levels of 10 steroids were quantitated on the SANIST platform: Aldosterone, Corticosterone, Cortisol, 11-deoxycortisol, Androstenedione, Testosterone, dehydroepiandrosterone, dehydroepiandrosterone sulfate (DHEAS), 17-OH-Progesterone and Progesterone. We performed both traditional and a machine learning analysis.

Conclusion: We show that the machine learning approach based on the steroid relationships developed here was much more accurate than the PSA, DHEAS, and direct absolute value match method in separating the PCa, BPH and control subjects, increasing the sensitivity to 90% and specificity to 84%. This technology, if applied in the future to a larger number of samples will be able to detect the individual enzymatic disequilibrium associated with the steroid ratio and correlate it with the disease. This learning machine approach could be valid in a personalized medicine setting.

Introduction

Hormonal physiology plays an important role in the function of the prostate. Development, growth, and functions of this gland depend on the actions of testosterone and its metabolite dihydrotestosterone (DHT), which enable the growth and proliferation of the glandular component of the prostate through binding and activation of androgen receptor (AR) within the cytoplasm of prostatic epithelial cells (1). In contrast to the chemical and virus-induced tumors, the hormone-related cancers shared a quite different mechanism of carcinogenesis: hormones, both endogenous and exogenous, drive cell proliferation, increasing the number of cell divisions and the opportunity for random genetic errors (2–4). In some cases, hormones can be metabolized to genotoxic agents by enzymes, which can form adducts with DNA and lead to mutations (5). The actions of androgen hormones are very controversial in prostatic hyperplasia [benign prostate hypertrophy (BPH)] and prostate cancer (PCa). Several studies showed that the activation of AR results in proliferative growth of prostatic epithelium, favoring the transition from prostatic hyperplasia (BPH) to PCa (4, 6, 7). However, lower levels of testosterone and DHT are associated with BPH and PCa (8–11). As men age, the free serum testosterone declines by 2–3% approximately annually (12); yet, there is increase in the incidence of BPH, and longevity increases the number at risk. Also, men who have high body mass index have larger prostates but lower testosterone levels (13). There is a potential explanation of lower levels of testosterone and DHT with BPH, as explained through a murine model where the authors blocked AR activation in the luminal cells, leading to inflammation of the prostate gland in these animals (14). Inflammation and subsequent angiogenesis lead to PCa progression (15), and we identified a fragment of complement C3f in low prostate-specific antigen (PSA) samples (16). PSA is a specific prostate marker, but only marginally influences mortality from PCa, and a recent study of the available PCa biomarkers suggested more candidate biomarkers (17).

Most of the studies measuring testosterone and DHT (cited above) used immunoenzymatic (ELISA) and radioimmuno-assay (RIA) to measure total testosterone, then applied Vermeulen’s formula (based on sex hormone-binding globulin and albumin) to calculate free testosterone. Analysis of testosterone in serum is quite complicated, especially in ELISA and RIA conditions, where the low amount and enzymatic cross-reactions seriously affect their quantitation (18). Recently, mass spectrometry (MS) has been elected as the reference technology for steroid dosage (19–23) since the sensitivity and specificity of this technology avoids the quantification problems in ELISA and RIA conditions. Steroids dosage also represents a crucial task for the diagnosis of different diseases (22, 24–27), such congenital adrenal hyperplasia (28), cardiovascular (22), and other degenerative diseases (24, 26, 27).

Liquid chromatography tandem MS operating in multiple reaction monitoring conditions (LC-MS/MS-MRM) is actually the most employed technology for steroids analyses (23, 29). Different methods have been developed in both electrospray (ESI) (30) and Atmospheric Pressure Ionization conditions (APCI) (31) combined to triple quadrupole (QQQ) as mass analyzer (32). QQQ has gained the consensus of clinical institutes mainly due to its higher quantitation accuracy in target analysis with respect to other analyzers (29, 32).

In addition to quantification of steroids by LC-MS/MS, other investigations show that the clinical information correlated with the disease is not only contained in the absolute steroid quantitation, but rather in their ratios (33, 34). This reflects their synthesis starting from the precursor cholesterol and the enzymes involved (Figure S1 in Supplementary Material). For example, a deficit in 21-hydroxylase leads to a decrease in mineralocorticoids and glucocorticoids and to an increase of androgens and testosterone (4). Thus, the ratio can be an informative and complementary diagnostic tool.

Surface Activated Chemical Ionization-Electrospray-NIST Bayesian model database search (SANIST) platform has been recently proposed by us as a powerful tool for personalized medicine (28, 35). Basically, it is an open platform based on three core points:

(a) Standardized kits for sample preparation to maximize the inter-laboratory data reproducibility.

(b) Use of SACI/ESI ionization source to reduce the instrument chemical noise and increase the number of detectable analytes.

(c) The SANIST data elaboration technology that is based on an open model: different applications are progressively added depending to the analysis target.

Concerning the point (c) SANIST has been currently implemented with two applications:

(a) An application to be used for biomarker discovery, able to classify samples and diagnose diseases on the basis of Bayesian machine learning approach (35). This technology makes possible to consider the individual variability and correct the diagnose bias obtained using fixed biomarkers (36, 37).

(b) An application to identify the analytes in conformity with the guidelines suggested by the European Community (EU directive 2002/657/EC) (35).

We describe here the development and application of a new SANIST platform to quantify steroids for diagnostic purposes going toward a personalized medicine approach. In detail, 10 steroids [aldosterone, corticosterone, cortisol, 11-deoxycortisol, androstenedione, testosterone, dehydroepiandrosterone (DHEA), dehydroepiandrosterone sulfate (DHEAS), 17-OH-progesterone and progesterone] have been quantitatively detected by means of SANIST technology and their profile has been analyzed for PCa diagnosis. The profiles of 71 serum samples (31 controls, 20 patients with PCa, and 20 subjects with benign prostate hyperplasia) were evaluated. The ability of the steroids profiling to discriminate among the different groups was tested both considering the individual percent concentration and the relative ratio relationships, showing that the predictive performances are better when using the relative relationship conditions.

Materials and Methods

Chemicals

A CE-IVD certified and standardized kit for steroids analysis was purchased from ISBN (Varese, Italy). The kit contains LC mobile phases (phases A and B), the sample precipitation reagent E and the sample resuspension reagent D. Aldosterone, Corticosterone, Cortisol, 11-deoxycortisol, Androstenedione, Testosterone, DHEA, DHEAS, 17-OH-Progesterone and Progesterone standards were purchased from (ISBN laboratory, Varese, Italy) together with their internal standard (Aldosterone d₈, Androstenedione ¹³C₃, Corticosterone d₄, Cortisol d₄, 11-Deoxycortisol d₅, DHEA d₆, Progesterone ¹³C₃, 17-OH-Progesterone ¹³C₃, DHEAS d₆, Testosterone ¹³C₃). Analyte controls at defined concentrations were purchased from (ISBN laboratory, Varese, Italy). Table S1 in Supplementary Material reports the internal standard concentrations while Table S2 in Supplementary Material has the known calibrate analyte concentrations and the calculated concentrations of both known and unknown samples.

Sample Selection, Collection, and Preparation

Samples were selected among 50-year-old male patients being subjected to a prostate biopsy for PCa diagnosis at the Urology Unit of the MultiMedica Castellanza, Varese, Italy, without considering those affected by autoimmune diseases, hypersensitivities, and other immune-mediated physical states. Multiple serum PSA analyses were performed on patients along with digital rectal exam and ultrasound. Many different clinical parameters were used to classify patients among which the histological analysis of the biopsy, total PSA, free-PSA, PSA-velocity (changes in PSA levels over time), ultrasound and rectal analysis. Controls included healthy subjects undergoing PSA dosage at our Urology unit that had a PSA lower than 2.5 ng/mL (28 patients) or less than 4 ng/mL and free PSA greater than 15% (3 patients). The study was approved by the institutional review board ethics committee (Protocol N 10 2 10/2011) and, according to the Helsinki Declaration of 1975 as revised in 2013, informed consent was provided by all patients. 71 serum samples (31 controls, 20 patients with PCa, and 20 subjects with benign prostate hyperplasia) were collected (sample number, age, PSA, AR positivity, and Gleason score are shown in Table 1). Blood samples were treated to obtain sera: they were collected in siliconized tubes, left to coagulate for 30 min and then subjected to centrifuge for 30 min at 3,000 rpm at 4°C to reduce variability through the collection of the samples (38). 1.8-mL cryotube vials (Thermo Fisher Scientific, Rodano, Milan, Italy) were used to collect sera and store them at −80°C until the time of analysis. 10 µL of steroids isotopically labeled internal standard at known concentration (Table S1 in Supplementary Material) were added to 290 µL of serum to normalize them. 600 µL of precipitating reagent E were added to each sample. They were then centrifuged for 1 min at 13,000 × g/min using a (Biospin, Eppendorf, Germany). 700 µL of the supernatant were collected in another Eppendorf tube and dried with nitrogen gas. The dried pellet was resuspended using 70 µL of reagent D. The analyte calibrants (Table S2 in Supplementary Material) were treated with the same procedure employed for the sample preparation. Quantitative analysis stability was checked by inserting standard controls at different known concentrations enclosed in the kit employed for the dosage (Table S1 in Supplementary Material). Limit of detection (LOD) and quantitation (LOQ) were calculated using the standard calibration curves (7 points).

TABLE 1

Table 1. Clinical characteristics of the (a) patients, (b) Benign prostate hypertrophy, and (c) controls subjects.

Chromatography

To study the selected analytes, a NeXera UPLC liquid chromatography apparatus (Shimatzu, San Jose, CA, USA) was employed. A Hypersylgold C18 2.1 mm × 50 mm 1.8 µm column was used. The ISBN kit has two mobile phases that were employed (Phases A and B). The LC apparatus worked in a binary gradient: 5% of B was kept for 0.5 min, after 3.5 min, B was increased to 45%, after other 8 min, B was raised to 65% and in 0.1 min to 100%. This percent was maintained for 1.9 min, then, 5% of B was achieved in 0.1 min and the column equilibrium was re-established reaching the starting conditions for 2.9 min. The total time of the analysis was 17 min. A 0.55 mL/min chromatographic flow was employed with a 20 µL injection volume.

Mass Spectrometry

Molecule assessments were achieved using 120 series triple quadrupole (PerkinElmer, USA) operating in SANIST mode. Heated SACI-ESI capillary had a 1,500 V voltage, SACI-ESI surface had a 47 V voltage, dry gas was kept to 2 L/min, 80 psi Nebulizer, and 40°C temperature were set. Collision-induced dissociation conditions were established for tandem MS experiments using Ar as the collision gas. The first quadrupole had an isolation window of ±0.3 m/z and the second one of 0.1 m/z.

Spectra Stability Evaluation

Chromatographic peak spectral average was performed to obtain the higher spectrum quality. Only signal intensities higher than 10³ counts/s were considered to improve the spectrum quality.

General Normalization

As a normalization factor, since the concentrations of cortisol and DHEAS were far above any of the other steroid values, they were reduced by dividing 10 and 100 times, respectively.

SANIST Data Elaboration Platform

The disease classification among controls, PCa, and subjects with benign prostate hyperplasia were obtained using two approaches: (a) singular absolute concentration match and (b) steroid ratio relationship match (Figure S2 in Supplementary Material). The SANIST tool acquires the selected potential steroid quantitation value and calculates the different combinatorial concentration ratio values. The steroids concentration ratio relationships are processed with the novel Bayesian data analysis-based SANIST mathematical model. SANIST data elaboration system calculates the probability that the detected steroid profile was associated to known disease (e.g., PCa). Biomarker profiles of serum from biopsy positive (PCa) and negative (BPH) subjects were analyzed and inserted to previously instruct the system. The probability in which an event occurs is determined by the Bayes’ theorem. Based on this interpretation, Bayes’ theorem is defined as the relationship between P(A) (probability that an acquired profile matches a PCa patient), P(B) (probability that a “fingerprint” biomarker profile corresponds to a biopsy negative subject), P(A | B) and P(B | A) for any events considering A and B in the same event space. In our case A and B indicate the different steroid concentration ratios. Following the Bayesian interpretation, probability, or uncertainty measures the confidence that something is true. Under this theory, Bayes’ theorem relates uncertainty before and after observing proof. The initial uncertainty in A is the prior, P(A). The uncertainty having accounted for evidence B is the posterior, P(A | B). The degree of support B provides for A is represented by the P(B | A)/P(B) ratio. Formula I shows the relationship among the probabilistic function:

$\text{P}(\text{A|B})\text{=}\frac{\text{P}(\text{B|A})\text{P}(\text{A})}{\text{P}(\text{B})}$ (I)

We introduced a new coefficient in formula (I) that gives formula (II):

$\text{P}(\text{A|B})\text{=}\frac{\text{P}(\text{B|A})\text{P}(\text{A})}{\text{P}(\text{B})}\ast C(mz,rt,i)$ (II)

The correction factor C (mz, rt, i) enables to progressively rise the disease classification accuracy, on the basis of the number of uncorrected sample classification. In our case, the correct classification of the samples as PCa biopsy positive or negative with the SANIST tool. At least 10,000 sequential classification steps are required to reach this result. It is possible to evaluate the probability that the data observed in the serum are correctly classifiable with the theorem considering the steroid levels and different combinatorial concentration ratios.

We used two models for classification of the individual patient profiles:

(a) Singular absolute concentration match. In this model, a vector database containing the analyte m/z ratio vs the singular analytes concentration and a percent value is obtained. The control, PCa, or benign prostate hyperplasia samples were classified, using the SANIST algorithm (35), by matching their concentration vector again the machine learning database excluding the identity match.

(b) Steroids ratio relationships match. In this model, all the ratios of the concentration of the analyzed steroids are obtained and a vector of ratio index [R(i)] is obtained for each sample. The vector machine learning database is obtained by classifying each vector as a control, patient with PCa, or with benign prostate hyperplasia. Each subject was classified, using the SANIST algorithm (39), by similarity of its ratio vector versus the database by excluding the identity match.

SANIST Application for Steroids Relationship Calculation Spectral Matching Algorithm

SANIST proprietary data elaboration platform configuration has shown in a previous publication (35). SANIST application converts the sample amount in ratio index using the formula (1):

$R(i,k)(k=\text{1 to }(n-\text{1}))[(i=(k+\text{1})\text{ to }(n))]=C(k)/C(i+\text{1})$ (1)

where n is the total number of analyte, i is an analyte counter, and R is the analyte concentration ratio given by the C concentration factor. The formula (1) is repeated for (n − 1) cycles so to calculate a number or R(i) ratio following the k index. The obtained ratios are employed in order to obtain a comparison vector with the (k;i) vs R(i;k). The number of R ratios depends on the number of analytes (n) and is given by the Eq. 2.

$N(R)((k=\text{1 to }(\text{}n-\text{1}))=\text{Sum }(n-i)$ (2)

SANIST Application for Steroid Singular Absolute Concentration Match

Classical data comparison was obtained by calculating the comparison vector considering the m/z vs C normalized (Cn) data. A m/z virtual parameter having a value of 1,000 was inserted in the vector as normalizing factor and its C value (Cv) was given by the sum of all the concentration parameters (see Eq. 3):

$\text{Cv}=\text{Sum }(i=\text{1 to }n)(Ci)$ (3)

The Cn normalized data is calculated as a relative concentration with respect to Cv (Eq. 4):

$\text{Cn}=(\text{Ci}/\text{Cv})\ast \text{1}00$ (4)

The final comparison vector is so given by m/z vs Cn pairs.

Results and Discussion

The LC-MS/MS-MRM data acquisition was optimized by monitoring the transitions reported in Table 2. Figure 1 reports the chromatographic profiles obtained by analyzing the 10 selected steroids injecting 20 µL of the calibration 7 standard solution (see Table S2 in Supplementary Material). Aldosterone and Cortisol are the only co-eluting steroids. Their retention time is 5.56 min. However, they are discriminated based on of their m/z parent ion difference (m/z 359.2 for aldosterone and m/z 363.1 for cortisol) and selective MRM transition (m/z 359.2—>189.1 and 359.2—>331.3 for aldosterone; m/z 363.1—>121, and 363.1—>115.1 for cortisol). DHEA and DHEAS exhibit the same parent ion, due to the fact that DHEAS loses its sulfate groups in the source. The discrimination between the two compounds occurs on the basis of LC retention time (9.30 min for DHEA and 6.89 min for DHEAS).

TABLE 2

Table 2. Table reporting the mass spectrometry (MS)/MS MRM analyte transitions.

FIGURE 1

Figure 1. Liquid chromatographic (LC) profiles obtained using the standards for the 10 selected steroids (20 µL) that were subjected to mass spectrometry (MS)/MS-MRM. Only two of the analytes eluted at the same time, but these were easily distinguished by the MS/MS-MRM.

Stability and performance of the method mainly in terms of matrix effect, LOD, LOQ, and linearity range were tested. Matrix effect was evaluated both based on internal standards stability and on precision and accuracy parameters (Table 3). Internal standard inter-analysis variation was lower than 5% while intra-day precision and accuracy error were lower than 6 and 8%, respectively. Even interday data show good reproducibility with precision and accuracy error lower than 8 and 9%, respectively. Table 3 reports also the LOD, LOQ, and linearity range for each analyte. The LOD are among 0.01 and 1 ng/mL depending by the analytes while LOQ among 0.05 and 5 ng/mL. Linearity range covers four orders of magnitude for each analyte.

TABLE 3

Table 3. Precision and accuracy error %, limit of detection (LOD), limit of quantitation (LOQ), and linearity range for each steroid.

After verifying the method stability, patient samples were analyzed. The absolute quantitation values for each steroid are reported in Table S2 in Supplementary Material together with the quantitation curves. The steroids concentrations are in different ranges based on the analyzed molecules, which are given for each patient in Table S2 in Supplementary Material. In particular, 11-deoxycortisol, 17-Oh-progesterone, aldosterone, androstenedione, corticosterone, DHEA, progesterone, testosterone concentration is under 50 ng/mL, while DHEAS and cortisol are much higher, in some cases can be over 1,000 ng/mL.

Since the PSA levels were used to determine the control group, this would introduce a bias into the system for determining the performance of PSA. Using PSA for the BPH and the PCa patients, the PSA receiver-operating characteristic (ROC) curve was calculated and area under the curve was 0.7118. We determined a threshold based on the lower limit of the 95 confidence intervals of the geometric mean of the PSA of the PCa patients, lead to a sensitivity (true positives) of 63% and specificity (true negatives) of 59%. The PSA, Gleason score, and AR receptor data of the cancer patients, benign prostate hypertrophy and control subjects are reported in Table 1A–C, respectively.

Most of the steroids were lower in the PCa patients than in the controls (both healthy and BPH) and are shown graphically in Figure 2. This is in keeping with the lower levels of free testosterone that are associated with “progression” or “reclassification” of PCa (9, 11). However, these data do not coincide with the lower levels of free testosterone in patients with BPH (8, 10) and the inflammatory state with the lack of AR in luminal cells in murine models (14). These studies used RIA conditions, which could affect the level of accuracy. Since more and more clinical studies will be using LC/MS-MS for analyzing androgen levels, our data will be possibly confirmed. Both the enzymes that synthesize androgens (40, 41) and the breakdown of androgens (42, 43) are present in PCa cells and may affect circulating androgen and steroid levels.

FIGURE 2

Figure 2. (A) Classic analysis of absolute concentrations of the steroids. (B) Receiver-operating characteristic (ROC) curve of dehydroepiandrosterone sulfate based on the classification in Table 1. (C) Analysis of the output of the machine learning (based on the classification in Table 1) and ROC of the values obtained (*p < 0.05; **p < 0.01; p < 0.001 one way ANOVA, Graphpad Prism 5).

There were several statistically significant differences between the hormones that separated the patients with prostate carcinoma versus those of benign prostate hyperplasia and controls. The only hormonal value separating statistically all three cases was the one of DHEAS. We determined a threshold based on the upper limit of the 95 confidence intervals of the geometric mean of the DHEAS of the PCa patients (Table 1). This identified 11/20 (55%) of the PCa patients, 16/20 (80%) of the BPH patients, and 19/31 (62%) of the healthy controls. The ROC curve was calculated and area under the curve (separating the PCa from BPH and HC) was 0.702 (Figure 2). The sensitivity was 55% and the specificity was 59%. We also performed a Principal Component Analysis (PCA) (44) on the steroid absolute values (Figure 3). PCA is a method where a multivariate data set is linearly transformed into a set of uncorrelated variables that are ordered in descending manner, the first few components that often explain large amount of the variation (44). BPH patients and controls overlapped extensively (Figure 3). The vast majority of the PCa patients (Figure 3) were grouped into a restricted subset of the area of the BPH patients and controls.

FIGURE 3

Figure 3. Principal component analysis (PCA) of the absolute (normalized) concentrations of all the steroids. Points labeled in red are PCa patients, in blue are the benign prostate hypertrophy patients, green are the healthy controls. Principal component (PC) 1 and PC 2 explain 34.4 and 18.5% of the total variance, respectively. The prediction ellipses are such that with a probability 0.95, a new observation from the same group will fall inside the ellipse.

We then proceeded in calculating the comparison vector employed to build the progressively updated supervised database that will be used for sample classification on the basis of individual steroid enzymatic profiles. When similarity vectors are calculated following the SANIST algorithm (39), the cortisol and DHEAS exhibits a high weight on the calculation function and the other concentration parameters area weakly considered by the discrimination function. To avoid this problem, we “normalized” the weight of the cortisol and DHEAS by dividing their value by 10 and 100, respectively. All the samples were subjected to this normalization model.

The database vectors are matched with the vectors of unclassified samples in order to obtain a diagnosis. The algorithm then provides for each subject the steroid ratio relationships with respect to the others, and a complex relational map is obtained by using this approach (Figure 4). The discriminating power is not given by the position of the nodes but from the density of their relationships that is measured in a normalized scale between 1 and 100.

FIGURE 4

Figure 4. Vector diagram of the steroid ratio relationships.

Two database search profile match models were employed to predict the sample classification as control, PCa, or benign prostate hyperplasia. A cartoon of the classical approach and the ratio approach used in the matches are shown in Figure S2 in Supplementary Material. In the case of steroid ratio concentration match, the angular difference among the vectors depends on the steroid concentration ratios, in the case of singular absolute concentration, it is correlated with the absolute concentration only.

We excluded each one of the samples and the learning machine used the other samples as a database to obtain the vectors. Then the sample that was excluded was categorized based on the vectors produced (Table 1). The average output of this learning machine is given in Table 1. There is a clear and significant separation of the PCa (18/20, 90%, all the false negatives were classified as HC), BPH (16/20, 80%, false positives were all PCa), and HC (25/31, 81%, false positives were all PCa) groups (Table 1; Figure 2). The ROC curves were calculated from these data and the area under the curve (again separating the PCa patients from BPH and HC) was 0.8353. The better performance was obtained using the steroid ratio relationship match model exhibiting a sensitivity of 90% and a specificity of 84%. In the case of singular absolute concentration match model the sensitivity (51%) and specificity (52%) are lower. There was no correlation with Gleason score or AR receptor status of any of the classification methods.

Another important information available through the steroid ratio concentration match model is that related to the couple of steroids that could lead to the enzyme classification. It is possible to extrapolate that by minimizing the Geometric Chebyshev distance by means of the SANIST technology (35, 39) and considering a ratio relational coefficient higher than 40 in relative value. Therefore, it is possible to correlate the altered steroids concentration ratios with the biological pathways and consequently with the enzyme potentially responsible for the disease. For example, Table S3 in Supplementary Material shows two steroid relationship ratios that can selectively discriminate patients with prostate carcinoma from benign prostate hyperplasia and healthy controls, together with the associated enzyme (P450C17). This fact opens a way to a future prospective in personalized medicine applications that will provide detailed information to the clinical specialist together with the enzymes that can be targets of certain drugs. This information makes it possible to formulate personalized therapies that re-equilibrate the biochemical disequilibrium. To enrich database with the most relevant individual enzymatic pathways to be correlated with the cancer development and progression, an extensive multi-center study on more subjects is necessary and will be the target of future works.

Conclusion

Here, we provide a method for quantifying accurately 10 steroids in serum and calculate their ratios. We noted that the levels in the serum of PCa patients were lower than that of BPH and healthy controls, in particular testosterone. Based on clinical and experimental data (8–11, 14), this would imply that supplementation of testosterone in aged men might not be harmful. We also provide a novel SANIST data elaboration application to elaborate steroid profile data in PCa patients and to correlate them with known disease in a learning machine approach. It is based on the evaluation of steroids ratio relationships instead of the simple absolute concentrations. Steroids ratio relationships performed better (sensitivity 90%, specificity 84%) than PSA (sensitivity 63%, specificity 59%) or individual steroids concentrations (sensitivity 51%, specificity 52%). Our study has one limitations, a small sample size. However, it opens a new way to provide to the clinic a diagnostic hypothesis based on metabolic profile together with the steroid concentrations and their ratio relationships that eventually indirectly indicate the enzymes involved in endocrine and metabolic disorders.

Ethics Statement

The institutional review board ethics committee of IRCCS MultiMedica approved the study (Protocol N 10 2 10/2011). All patients enrolled in the study signed informed consent and accordingly to the Helsinki Declaration of 1975 as revised in 2013.

Author Contributions

AA: designed the study, conceived the experiments, performed the analysis, wrote the manuscript. AB: conceived the experiments, processed the clinical samples, performed the statistical analysis, and wrote the manuscript. BB: processed the clinical samples, managed the clinical database, and performed the statistical analysis. GA: provided clinical support-patients histological classification GP: provided clinical support, histology evaluation. PC: provided clinical samples, clinical data and the clinical support, and wrote the manuscript. LC: contributed to the MS data analysis. MC: provided important chemical-physical contribution for steroid mass spectra data elaboration. SC: conceived the experiments, run the MS, performed the statistical and bioinformatic analysis, wrote the manuscript. DN: conceived the experiments, performed the statistical analysis, and wrote the manuscript. All authors confirmed they have contributed to the intellectual content of this paper and have met the following three requirements: (a) significant contributions to the conception and design, acquisition of data, or analysis and interpretation of data; (b) drafting or revising the article for intellectual content; and (c) final approval of the published article.

Conflict of Interest Statement

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.

The reviewer EZ and handling Editor declared their shared affiliation.

Acknowledgments

The authors thank Sig. Maria Florio for her support and Dr. Rosanna Falbo for English editing. We thank Simone Libertini and Paola Corradino for bibliography and assistance with the manuscript. We are grateful to Professor Giancarlo Comeri and Dott.ssa Stefania Comini for support and collaboration to the clinical study.

Funding

These studies were supported by a grant from the AIRC (Associazione Italiana per la Ricerca sul Cancro, IG15895 to DN) and a grant of the Italian Ministry of Health (MS2010 to Paola Erba and AA). AB was a Fondazione Umberto Veronesi (FUV) fellow. AB is part of a Ph.D. iTalent program. The funding sources had no influence on study design, data collection, analysis and interpretation, writing of the report and on the decision to submit the article for publication.

Supplementary Material

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

Figure S1. Schematic representation for steroid synthesis starting from the precursor cholesterol and the enzymes involved.

Figure S2. Scheme of singular absolute concentration match and the steroids ratio concentration match models.

Table S1. The internal standard concentrations.

Table S2. The calibrate analyte concentrations and the calculated concentrations of both known and unknown samples.

Table S3. Two steroid relationship ratios that can selectively discriminate patients with prostate carcinoma from benign prostate hyperplasia and healthy controls, together with the associated enzyme (P450C17).

References

1. Hayward SW, Cunha GR. The prostate: development and physiology. Radiol Clin North Am (2000) 38:1–14. doi:10.1016/S0033-8389(05)70146-9

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Black A, Pinsky PF, Grubb RL III, Falk RT, Hsing AW, Chu L, et al. Sex steroid hormone metabolism in relation to risk of aggressive prostate cancer. Cancer Epidemiol Biomarkers Prev (2014) 23:2374–82. doi:10.1158/1055-9965.EPI-14-0700

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Schenk JM, Till C, Hsing AW, Stanczyk FZ, Gong Z, Neuhouser ML, et al. Serum androgens and prostate cancer risk: results from the placebo arm of the prostate cancer prevention trial. Cancer Causes Control (2016) 27:175–82. doi:10.1007/s10552-015-0695-0

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Srigley JR. Benign mimickers of prostatic adenocarcinoma. Mod Pathol (2004) 17:328–48. doi:10.1038/modpathol.3800055

PubMed Abstract | CrossRef Full Text | Google Scholar

5. D’Uva G, Baci D, Albini A, Noonan DM. Cancer chemoprevention revisited: cytochrome P450 family 1B1 as a target in the tumor and the microenvironment. Cancer Treat Rev (2017) 63:1–18. doi:10.1016/j.ctrv.2017.10.013

CrossRef Full Text | Google Scholar

6. Oto A, Kayhan A, Jiang Y, Tretiakova M, Yang C, Antic T, et al. Prostate cancer: differentiation of central gland cancer from benign prostatic hyperplasia by using diffusion-weighted and dynamic contrast-enhanced MR imaging. Radiology (2010) 257:715–23. doi:10.1148/radiol.10100021

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Schenk JM, Kristal AR, Arnold KB, Tangen CM, Neuhouser ML, Lin DW, et al. Association of symptomatic benign prostatic hyperplasia and prostate cancer: results from the prostate cancer prevention trial. Am J Epidemiol (2011) 173:1419–28. doi:10.1093/aje/kwq493

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Kristal AR, Schenk JM, Song Y, Arnold KB, Neuhouser ML, Goodman PJ, et al. Serum steroid and sex hormone-binding globulin concentrations and the risk of incident benign prostatic hyperplasia: results from the prostate cancer prevention trial. Am J Epidemiol (2008) 168:1416–24. doi:10.1093/aje/kwn272

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Biolchi V, Silva Neto B, Pianta DB, Koff WJ, Berger M, Brum IS. Androgen receptor GGC polymorphism and testosterone levels associated with high risk of prostate cancer and benign prostatic hyperplasia. Mol Biol Rep (2013) 40:2749–56. doi:10.1007/s11033-012-2293-5

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Jarvis TR, Chughtai B, Kaplan SA. Testosterone and benign prostatic hyperplasia. Asian J Androl (2015) 17:212–6. doi:10.4103/1008-682X.140966

PubMed Abstract | CrossRef Full Text | Google Scholar

11. San Francisco IF, Rojas PA, DeWolf WC, Morgentaler A. Low free testosterone levels predict disease reclassification in men with prostate cancer undergoing active surveillance. BJU Int (2014) 114:229–35. doi:10.1111/bju.12682

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Harman SM, Metter EJ, Tobin JD, Pearson J, Blackman MR; Baltimore Longitudinal Study of Aging. Longitudinal effects of aging on serum total and free testosterone levels in healthy men. Baltimore longitudinal study of aging. J Clin Endocrinol Metab (2001) 86:724–31. doi:10.1210/jcem.86.2.7219

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Klaassen Z, Howard LE, Moreira DM, Andriole GL Jr, Terris MK, Freedland SJ. Association of obesity-related hemodilution of prostate-specific antigen, dihydrotestosterone, and testosterone. Prostate (2017) 77:466–70. doi:10.1002/pros.23285

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Zhang B, Kwon OJ, Henry G, Malewska A, Wei X, Zhang L, et al. Non-cell-autonomous regulation of prostate epithelial homeostasis by androgen receptor. Mol Cell (2016) 63:976–89. doi:10.1016/j.molcel.2016.07.025

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Araldi EM, Dell’aica I, Sogno I, Lorusso G, Garbisa S, Albini A. Natural and synthetic agents targeting inflammation and angiogenesis for chemoprevention of prostate cancer. Curr Cancer Drug Targets (2008) 8:146–55. doi:10.2174/156800908783769382

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Fania C, Sogno I, Vasso M, Torretta E, Leone R, Bruno A, et al. A PSA-guided approach for a better diagnosis of prostatic adenocarcinoma based on MALDI profiling and peptide identification. Clin Chim Acta (2015) 439:42–9. doi:10.1016/j.cca.2014.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Amaro A, Esposito AI, Gallina A, Nees M, Angelini G, Albini A, et al. Validation of proposed prostate cancer biomarkers with gene expression data: a long road to travel. Cancer Metastasis Rev (2014) 33:657–71. doi:10.1007/s10555-013-9470-4

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Tate J, Ward G. Interferences in immunoassay. Clin Biochem Rev (2004) 25:105–20.

PubMed Abstract | Google Scholar

19. Adaway JE, Keevil BG, Owen LJ. Liquid chromatography tandem mass spectrometry in the clinical laboratory. Ann Clin Biochem (2015) 52:18–38. doi:10.1177/0004563214557678

CrossRef Full Text | Google Scholar

20. Griffiths WJ, Abdel-Khalik J, Yutuc E, Morgan AH, Gilmore I, Hearn T, et al. Cholesterolomics: an update. Anal Biochem (2017) 524:56–67. doi:10.1016/j.ab.2017.01.009

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Sabbagh B, Mindt S, Neumaier M, Findeisen P. Clinical applications of MS-based protein quantification. Proteomics Clin Appl (2016) 10:323–45. doi:10.1002/prca.201500116

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Yeap BB. Sex steroids and cardiovascular disease. Asian J Androl (2014) 16:239–47. doi:10.4103/1008-682X.122357

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Soldin SJ, Soldin OP. Steroid hormone analysis by tandem mass spectrometry. Clin Chem (2009) 55:1061–6. doi:10.1373/clinchem.2007.100008

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Pannu PK, Calton EK, Soares MJ. Calcium and vitamin D in obesity and related chronic disease. Adv Food Nutr Res (2016) 77:57–100. doi:10.1016/bs.afnr.2015.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Peskoe SB, Joshu CE, Rohrmann S, McGlynn KA, Nyante SJ, Bradwin G, et al. Circulating total testosterone and PSA concentrations in a nationally representative sample of men without a diagnosis of prostate cancer. Prostate (2015) 75:1167–76. doi:10.1002/pros.22998

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Pongkan W, Chattipakorn SC, Chattipakorn N. Roles of testosterone replacement in cardiac ischemia-reperfusion injury. J Cardiovasc Pharmacol Ther (2016) 21:27–43. doi:10.1177/1074248415587977

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Wallace TC, Slavin M, Frankenfeld CL. Systematic review of anthocyanins and markers of cardiovascular disease. Nutrients (2016) 8:E32. doi:10.3390/nu8010032

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Cristoni S, Cuccato D, Sciannamblo M, Bernardi LR, Biunno I, Gerthoux P, et al. Analysis of 21-deoxycortisol, a marker of congenital adrenal hyperplasia, in blood by atmospheric pressure chemical ionization and electrospray ionization using multiple reaction monitoring. Rapid Commun Mass Spectrom (2004) 18:77–82. doi:10.1002/rcm.1284

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Wozniak B, Matraszek-Zuchowska I, Zmudzki J. LC-MS/MS fast analysis of androgenic steroids in urine. Anal Bioanal Chem (2012) 403:2965–72. doi:10.1007/s00216-012-5859-0

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Higashi T, Ogawa S. Chemical derivatization for enhancing sensitivity during LC/ESI-MS/MS quantification of steroids in biological samples: a review. J Steroid Biochem Mol Biol (2016) 162:57–69. doi:10.1016/j.jsbmb.2015.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Yu T, Xu H, Wang W, Li S, Chen Z, Deng H. Determination of endogenous corticosterone in rodent’s blood, brain and hair with LC-APCI-MS/MS. J Chromatogr B Analyt Technol Biomed Life Sci (2015) 1002:267–76. doi:10.1016/j.jchromb.2015.08.035

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Raro M, Portoles T, Pitarch E, Sancho JV, Hernandez F, Garrostas L, et al. Potential of atmospheric pressure chemical ionization source in gas chromatography tandem mass spectrometry for the screening of urinary exogenous androgenic anabolic steroids. Anal Chim Acta (2016) 906:128–38. doi:10.1016/j.aca.2015.11.041

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Lee SH, Lee N, Hong Y, Chung BC, Choi MH. Simultaneous analysis of free and sulfated steroids by liquid chromatography/mass spectrometry with selective mass spectrometric scan modes and polarity switching. Anal Chem (2016) 88:11624–30. doi:10.1021/acs.analchem.6b03183

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Zhang Q, Chen Z, Chen S, Xu Y, Deng H. Intraindividual stability of cortisol and cortisone and the ratio of cortisol to cortisone in saliva, urine and hair. Steroids (2017) 118:61–7. doi:10.1016/j.steroids.2016.12.008

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Cristoni S, Dusi G, Brambilla P, Albini A, Conti M, Brambilla M, et al. SANIST: optimization of a technology for compound identification based on the European union directive with applications in forensic, pharmaceutical and food analyses. J Mass Spectrom (2017) 52:16–21. doi:10.1002/jms.3895

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Cruz JA, Wishart DS. Applications of machine learning in cancer prediction and prognosis. Cancer Inform (2007) 2:59–77.

Google Scholar

37. Menden MP, Iorio F, Garnett M, McDermott U, Benes CH, Ballester PJ, et al. Machine learning prediction of cancer cell sensitivity to drugs based on genomic and chemical properties. PLoS One (2013) 8:e61318. doi:10.1371/journal.pone.0061318

PubMed Abstract | CrossRef Full Text | Google Scholar

38. D’Angelo G, Villa C. Comparison between siliconized evacuated glass and plastic blood collection tubes for prothrombin time and activated partial thromboplastin time assay in normal patients, patients on oral anticoagulant therapy and patients with unfractioned heparin therapy. Int J Lab Hematol (2011) 33:219–25. doi:10.1111/j.1751-553X.2010.01271.x

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Albini A, Briga D, Conti M, Bruno A, Farioli D, Canali S, et al. SANIST: a rapid mass spectrometric SACI/ESI data acquisition and elaboration platform for verifying potential candidate biomarkers. Rapid Commun Mass Spectrom (2015) 29:1703–10. doi:10.1002/rcm.7270

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Hofland J, van Weerden WM, Dits NF, Steenbergen J, van Leenders GJ, Jenster G, et al. Evidence of limited contributions for intratumoral steroidogenesis in prostate cancer. Cancer Res (2010) 70:1256–64. doi:10.1158/0008-5472.CAN-09-2092

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Nakamura Y, Suzuki T, Nakabayashi M, Endoh M, Sakamoto K, Mikami Y, et al. In situ androgen producing enzymes in human prostate cancer. Endocr Relat Cancer (2005) 12:101–7. doi:10.1677/erc.1.00914

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Belledant A, Hovington H, Garcia L, Caron P, Brisson H, Villeneuve L, et al. The UGT2B28 sex-steroid inactivation pathway is a regulator of steroidogenesis and modifies the risk of prostate cancer progression. Eur Urol (2016) 69:601–9. doi:10.1016/j.eururo.2015.06.054

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Khvostova EP, Otpuschennikov AA, Pustylnyak VO, Gulyaeva LF. Gene expression of androgen metabolising enzymes in benign and malignant prostatic tissues. Horm Metab Res (2015) 47:119–24. doi:10.1055/s-0034-1374631

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Metsalu T, Vilo J. ClustVis: a web tool for visualizing clustering of multivariate data using principal component analysis and heatmap. Nucleic Acids Res (2015) 43:W566–70. doi:10.1093/nar/gkv468

PubMed Abstract | CrossRef Full Text | Google Scholar