Original Research ARTICLE
Serum Steroid Ratio Profiles in Prostate Cancer: A New Diagnostic Tool Toward a Personalized Medicine Approach
- 1IRCCS MultiMedica, Milan, Italy
- 2Department of Medicine and Surgery, University Milano-Bicocca, Milan, Italy
- 3Department of Oncology and Hemato-Oncology, University of Milan, Milan, Italy
- 4I.S.B.—Ion Source & Biotechnologies, Bresso, Italy
- 5Sant’Orsola Hospital, Bologna, Italy
- 6Department of Biotechnologies and Life Sciences, University of Insubria, Varese, Italy
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.
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
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 d8, Androstenedione 13C3, Corticosterone d4, Cortisol d4, 11-Deoxycortisol d5, DHEA d6, Progesterone 13C3, 17-OH-Progesterone 13C3, DHEAS d6, Testosterone 13C3). 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. Clinical characteristics of the (a) patients, (b) Benign prostate hypertrophy, and (c) controls subjects.
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.
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 103 counts/s were considered to improve the spectrum quality.
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:
We introduced a new coefficient in formula (I) that gives formula (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):
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.
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):
The Cn normalized data is calculated as a relative concentration with respect to Cv (Eq. 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).
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. 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. (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. 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.
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.
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.
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.
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.
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.
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.
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).
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
Keywords: steroids, steroids profile, prostate cancer, algorithm, EU law
Citation: Albini A, Bruno A, Bassani B, D’Ambrosio G, Pelosi G, Consonni P, Castellani L, Conti M, Cristoni S and Noonan DM (2018) Serum Steroid Ratio Profiles in Prostate Cancer: A New Diagnostic Tool Toward a Personalized Medicine Approach. Front. Endocrinol. 9:110. doi: 10.3389/fendo.2018.00110
Received: 15 January 2018; Accepted: 05 March 2018;
Published: 05 April 2018
Edited by:Gabriella Castoria, Università degli Studi della Campania “Luigi Vanvitelli” Naples, Italy
Reviewed by:Wen Zhou, Leonard M. Miller School of Medicine, United States
Erika Di Zazzo, Università degli Studi della Campania “Luigi Vanvitelli” Caserta, Italy
Luciana Marinelli, Università degli Studi di Napoli Federico II, Italy
Copyright: © 2018 Albini, Bruno, Bassani, D’Ambrosio, Pelosi, Consonni, Castellani, Conti, Cristoni and Noonan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Simone Cristoni, email@example.com