http://orcid.org/0000-0002-5492-9355
Figures
Abstract
Novel biomarkers for prostate cancer (PCa) diagnosis and prognosis are necessary to improve the accuracy of current ones employed in clinic. We performed a retrospective study between the association of several polymorphisms in the main genes involved in the synthesis and metabolism of sex hormones and PCa risk and aggressiveness. A total of 311 Caucasian men (155 controls and 156 patients) were genotyped for 9 SNPs in AR, CYP17A1, LHCGR, ESR1 and ESR2 genes. Diagnostic PSA serum levels, Gleason score, tumor stage, D´Amico risk and data of clinical progression were obtained for patients at the moment of the diagnosis and after 54 months of follow-up. Chi-squared test were used for comparisons between clinical variables groups, logistic regression for clinical variables associations between SNPs; and Kaplan–Meier for the association between SNPs and time to biochemical progression. We found 5 variants (CYP17A1) rs743572, rs6162, rs6163; (LHCGR) rs2293275 and (ESR2) rs1256049 that were statistically significant according to clinical variables (PSA, D´Amico risk and T stage) on a case-case analysis. Moreover, the presence of A and G alleles in rs743572 and rs6162 respectively, increase the risk of higher PSA levels (>10 ng/μl). With respect to D´Amico risk rs743572 (AG-GG), rs6162 (AG-AA) and rs6163 (AC-AA) were associated with an increased risk; and last, AC and AA genotypes for rs6163 were associated with a shorter biochemical recurrence free survival (BRFS) in patients with radical prostatectomy. In multigene analysis, several variants in SNPs rs2293275, rs6152, rs1062577, rs6162, rs6163, rs1256049 and rs1004467 were described to be associated with a more aggressiveness in patients. However, none of the selected SNPs show significant values between patients and controls. In conclusion, this study identified inherited variants in genes CYP17A1, LHCGR and ESR2 related to more aggressiveness and/or a poor progression of the disease. According to this study, new promise PCa biomarkers for clinical management could be included in these previous SNPs.
Citation: Robles-Fernandez I, Martinez-Gonzalez LJ, Pascual-Geler M, Cozar JM, Puche-Sanz I, Serrano MJ, et al. (2017) Association between polymorphisms in sex hormones synthesis and metabolism and prostate cancer aggressiveness. PLoS ONE 12(10): e0185447. https://doi.org/10.1371/journal.pone.0185447
Editor: Chih-Pin Chuu, National Health Research Institutes, TAIWAN
Received: July 6, 2017; Accepted: September 12, 2017; Published: October 5, 2017
Copyright: © 2017 Robles-Fernandez et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: The author(s) received no specific funding for this work.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Prostate cancer (PCa) is one of the most prevalent cancers diagnosed in men with 1.1 million cases worldwide in 2012 [1]. However, few conclusive studies have been performed with regard to the genetics of this cancer. One of the main challenges is to find new specific biomarkers that allow clinicians to detect the disease at an early stage, refine risk stratification, and control the course of patients.
So far, only clinical risk classifications, such as that of D´Amico [2], have gained enough evidence to be implemented in routine practice. However, a limitation of these classifications is the lack of integration with other known risk factors and genomic data, which could provide a more personalized risk assessment. Although a number of genomic classifications, such as FGFR1, CDKN1A and PMP22 genes [3], have demonstrated their ability to differentiate between low and high-risk patients [3–5], none of them are currently adopted in routine clinical practice [6]. In the same way, there are not genetic biomarkers in the clinical to predict the outcome to different treatments used to attend PCa patients [7].
Androgens play a pivotal role in the development and function of prostate as well as in the pathogenesis and progression of PCa [8]. Moreover, experimental and epidemiological data have suggested that also estrogen signaling may contribute to PCa development and progression [9]. Although recent studies have evaluated polymorphisms in sex hormone metabolism genes, such as AR, SRD5A2, CYP17A1 or ESR as PCa risk factors [10–12] they did not provide conclusive results.
The aim of the present study was to analyze several SNPs of genes related to the synthesis and metabolism of steroid hormones, such as CYP17A1 gene (rs743572, rs6162, rs6163 and rs1004467); luteinizing hormone chorionic gonadotropin hormone receptor (LHCGR) gene (rs2293275); androgen receptor (AR) gene (rs6152 and rs9332696); and estrogen receptor (ESR1 and ESR2) genes (rs1062577 and rs1256049, respectively). We designed a genetic analysis in these 9 SNPs in order to find out genetic associations with PCa risk, biochemical recurrence and/or clinical stratification.
Materials and methods
Patients
From 2012 to 2014, a total of 311 subjects with PSA levels ≥4.0 ng/mL meeting the criteria for a prostate biopsy were included in this study (Table 1). Patients with positive biopsy were analyzed for T stage, serum PSA, Gleason score and were classified according to D´Amico risk classification (low, intermediate and high risk). Negative biopsy individuals were considered as controls. All individuals underwent a systematic 20-core ultrasound guided biopsy in order to limit the false negative rate. After primary therapy PSA was monitored every 3 or 6 months to evaluate the existence of biochemical recurrence. All subjects of the study provided a written informed consent to be enrolled, which was previously approved by the Research Ethics Committee of Granada Center (CEI-Granada) following Helsinki ethical declaration.
SNPs selection and genotyping
Immediately before the biopsy, a peripheral blood sample was extracted from patients, put into EDTA coated tubes and stored at -200 C until genomic DNA extraction. A standard organic extraction procedure by phenol/chloroform/isoamyl alcohol and proteinase K, followed by purification with Microcon H 100 filters (Millipore, Germany) was used. To determine extracted DNA purity and concentration a NanoDrop 2000c (Thermo Scientific, USA) was used. Thereafter, DNA was stored at -200 C until genotyping. Five genes involved both directly or indirectly, in androgen synthesis and/or its metabolism were selected: CYPI7A1, AR, LHCGR, ESR1 and ESR2 and the picked SNPs were rs6162, rs743572, rs6163 and rs1004467 for CYPI7A1; rs6152 and rs9332969 for AR; rs2293275 for LHCGR; rs1062577 for ESR1; and rs1256049 for ESR2 (S1 Table). SNPs in these genes were selected using The National Center for Biotechnology Information website [13].
DNA genotyping was performed using TaqMan® Genotyping Master Mix (Applied Biosystems, USA) which included all essential components (except probes, templates and water) for polymerase chain reaction (PCR). Allelic discrimination assays were carried out in a 7900HT Fast Real-Time PCR System (Applied Biosystems, USA). Results were read using SDS software v.2.4 (Applied Biosystems, USA). In order to warrant the results of genotyping we carried out a validation assay by Sanger sequencing (S2 Table). The validation cohort comprised a 3% of the genotyped samples (randomly selected) for each single SNP.
Statistical analysis
Software package SPSS v.22 was used for statistical analyses (IBM Corporation, USA). The analyses included chi-square and Fisher exact tests for small samples size. The association between clinical variables and SNPs were analyzed by a binary logistic regression using different genetic models. Binary logistic regression was adjusted for PSA levels at diagnosis time, Gleason score, T stage and/or age. The biochemical recurrence free survival (BRFS) interval was estimated with Kaplan–Meier analysis and significance was determined by log-rank test. Cox regression analysis was used to assess the association of genetic variants and BRFS adjusting for PSA level at diagnosis, Gleason score, T stage and age Statistical signification was considered with p values ≤ 0.05. Genotypes analyses as well as Hardy–Weinberg equilibrium and Linkage disequilibrium (LD) analyses were performed using the online SNPStats software [14]. SNPs were in LD when they had a value of r2 > 0.5.
Results
Nine SNPs across AR, CYP17A1, LHCGR, ESR1 and ESR2 genes were successfully genotyped in the total cohort (n = 311).
In the case of rs9332969 (AR gene is in X chromosome), all individuals were hemizygous for the G allele so it was discarded for later analyses.
Hardy-Weinberg analysis showed that all controls and patients were in equilibrium except for rs6152 (AR gene is in X chromosome). LD analysis showed a strong linkage between rs743572, rs6162 and rs6163, all of them in CYP17A1 gene (S3 Table).
Case-control study
Firstly, analysis was conducted as a case/control study, but none of the selected SNPs showed significant differences between PCa-confirmed patients and controls (data not shown). Genotyping and allelic distribution in this cohort is shown in Table 2.
Case-case study
Analysis of clinical variables.
Secondly, we focused the analysis in PCa-confirmed patient population. We found statistically significant data in relation to PSA values (4-10ng/ml; 10-20ng/ml and >20 ng/ml) in several SNPs of the CYP17A1 gene, such as rs743572 (p = 0.003); rs6162 (p = 0.015) and rs6163 (p = 0.010). Furthermore, we observed that the presence of A allele in rs743572 (p = 0.032; OR = 3.017, 95% CI (1.102–8.260); and the presence of G allele in rs6162 (p = 0.036; OR = 3.129; 95% CI (1.076–9.103) presented an increased risk of PSA values above 10 ng/ml versus GG and AA homozygous individuals, respectively. Concerning to Gleason score, we did not find any marker with statistically significant value, though rs1256049 in ESR2 gene was close to significance (p = 0.076) (S1 Fig). Regarding T stage, CT patients for rs2293275 in LHCGR gene were associated with more advanced stages, T3 and 4 (p = 0.037). See details in Tables 3 and 4.
We found statistical associations between D´Amico risk classification and the following SNPs in the CYP17A1 gene: rs743572 (p = 0.008), rs6162 (p = 0.019), rs6163 (p = 0.015), and rs1256049 (p = 0.029) (Table 3). Moreover G allele presence in rs743572 confers an increased risk of being in intermediate and high risk stratification versus AA patients (p = 0.016; OR = 3.856; 95% CI (1.281–11.603)). The presence of A allele in rs6162 confers an increased risk of being in intermediate and high risk stratification versus GG (p = 0.032; OR = 3.574; 95% CI (1.117–11.439)); and in rs6163 A allele confers an increased risk of being in intermediate and high risk stratification versus CC (p = 0.019; OR = 3.866; 95% CI (1.248–11.981)) (Table 4).
Analysis of biochemical recurrence.
A total of 148 patients were included in the biochemical recurrence analysis after initiation of primary treatment. Of them, 32 (21.6%) received ADT (androgen deprivation therapy); 73 (49.3%) and 38 (25.7%) patients underwent radiotherapy and radical prostatectomy, respectively, and only 5 patients (3.4%) remained in active surveillance. In general, within 54 months of observation after treatment, 123 (83.1%) patients did not show significant increases in PSA levels and 25 (16.9%) of them presented a biochemical recurrence.
After stratifying by the type of therapy received, we observed that the rate of biochemical recurrence was 31.3% for patients treated with ADT, 6.8% for those treated with radiotherapy and 26.3% for patients undergoing radical prostatectomy. None of the patients who remained in active surveillance manifested biochemical recurrence.
As SNPs rs743572, rs6162 and rs6163 showed an increased risk of D´Amico risk, we performed a Kaplan Meier analysis and log-rank test for ADT, radiotherapy and radical prostatectomy therapies. For ADT as well as for radiotherapy, none of the SNPs were significantly associated to BRFS. However, for radical prostatectomy AC and AA genotypes in rs6163 (CYPI7A1) were significantly associated with a shorter BRFS compared to the CC genotype (log-rank p = 0.039), 29.10 months vs 49.59 months, respectively (Fig 1). When a Cox regression multivariable analysis was performed, rs6163 was not an independent variable for risk to BRFS after radical prostatectomy (p = 0. 221; OR = 2.822; 95% CI (0.535–14.885)).
Kaplan-Meier curves of time to biochemical recurrence in patients treated with radical prostatectomy and stratified by rs6163 genotype (CC vs. AC+AA). P value obtained from log-rank t test.
Multigene analysis.
Multigene analysis of SNPs revealed that the genotype GTGTAACA for variants rs743572, rs2293275, rs6152, rs1062577, rs6162, rs6163, rs1256049 and rs1004467 (in this order) was associated with significant increased risk according to D´Amico classification (p = 0.0045; OR: 0.6, 95% CI(0.19–1.02)); higher PSA values (p = 0.0011; OR: 0.68, 95% CI (0.35–1.31)); and increased risk of Gleason score ≥7 (p = 0.026; OR: 5.21, 95% CI (1.23–22.02)). Regarding T stage, patients with genotype GCGTAACA (rs743572, rs2293275, rs6152, rs1062577, rs6162, rs6163, rs1256049 and rs1004467, in this order) showed an increased risk of having higher T stage values such as T3 and T4 (p = 0.045; OR: 6.34, 95% CI (1.06–37.88)) (S4 Table).
Conversely, TAAG genotype for SNPs rs1062577, rs6162, rs6163 and rs1004467 seem to protect from higher PSA values (p = 0.016; OR: - 0.43, 95% CI (-0.79–-0.08)) and ATGGCCA genotype for SNPs rs743572, rs2293275, rs6152, rs6162, rs6163, rs1256049 and rs1004467 showed a protective effect for higher Gleason scores (values ≥7) (p = <0.0001; OR: -0.38, 95% CI (-0.37–-0.40)). No genotypes were significantly associated with risk of PCa (S4 Table).
Discussion
PCa is a heterogeneous disease as evidenced by its variable clinical course [15]. PSA level, core biopsies, T stage and Gleason scores used for initial evaluation offer limited information to clinicians to determine diagnosis and harshness of disease [16]. Current clinical risk groups, used in clinical routine, seem to misclassify patients leading to over/undertreatment of these patients [16]. New and reliable tools are needed to improve the precision in diagnosis and stratification of PCa patients, and genetic markers could be the most suitable ones. In the last decade, development of high throughput technologies have favored the identification of genetic variations associated with PCa and their incorporation into clinical practice offers an opportunity to ease clinical decisions [17].
There is scarce data in relation to genetic germline biomarkers for PCa prognosis and stratification. Main researches are focused on RNAseL (locus HPC1), ELAC2 (locus HPC2), MSR1 (chromosome 8) [18] and BRCA1/2 genes [19], but no conclusive data have been reported. For instance, it is known that PCa patients harboring germline DNA repair alterations generally have a worse clinical evolution or earlier cancer events. However, little is known of a genetic stratification of the disease or even genetic responsiveness to radical prostatectomy, radiotherapy or hormonal therapies including ADT and second-generation hormonal agents (such as abiraterone and enzalutamide) [20]. PCa is a hormone-dependent cancer; the androgen receptor (AR) axis plays a pivotal role in both disease development and progression [21]; for this reason we focused on PCa sex hormones related genes, to perform a retrospective study for the relation among AR, CYP17A1, LHCGR and ESR polymorphisms with PCa predisposition and severity.
We could not prove any statistically significant difference between controls and patients but when the analysis was focused on PCa patients we proved these genetic markers had a predictive role on PCa aggressiveness characteristics (PSA, D´Amico risk and T stages). It is known that any change in androgen synthesis and metabolism genes can strongly affect the progression of PCa and the response to treatments [22]. Our aim is to try to discover optimal biomarkers associated to PCa aggressiveness such as recent patented genes like rs4054823 (17p12.) [23, 24].
Recent studies of NGS or Exome sequencing focused on finding new PCa biomarkers, have found that rs33999879 (SMC4) was a predictor for Gleason scores upgrade [25]. These NGS data also found that carrying any mutations at pathogenic germline variants (ATM, ATR, BRCA2, FANCL, MSR1, MUTYH, RB1, TSHR and WRN) were frequently observed in patients with metastatic CRPC (castration-resistant PCa) [26]. However, no data yet had analyzed the role of AR, CYP17A1, LHCGR and ESR genes in PCa.
For CYP17A1 gene we found that several SNPs (rs743572, rs6162 and rs6163) were statistically associated with a more aggressive PCa (PSA values > 10 ng/ml and a higher D´Amico risk). CYP17A1 gene encodes a key enzyme in the steroidogenic pathway that produces progestins, mineralocorticoids, glucocorticoids, androgens, and estrogens, and it has a critical function in PCa [27]. There is no published data that had previously established any relationship between clinical stages and CYP17A1 SNPs in African ancestry or Caucasian populations [28].
In silico analyses for CYP17A1 variants have showed that rs1004467 is on intronic region, rs743572 is on 5´ untranslated region; and rs6162 and rs6163, are both coding-synonymous (H46H and S65S, respectively) and with a tolerated phenotypic effect [29] (SITF score: 0.31 and 1, respectively). The presence of CYP17A1 SNPs is further associated with altered levels of circulating DHEA-S (dehydroepiandrosterone sulfate) in Caucasians, which likely modify steroid precursor levels available for intracrine conversion to more potent hormones in tissues and prostatic cells [7]. In addition, genetic polymorphisms in the CYP17A1 gene have been significantly associated with a risk of progression to CRPC [30], but there are no details on risk stratification. A recent meta-analysis carried out by Wang et al, did not find any significant association between rs743572 polymorphism and PCa risk but it was suggested that CYP17A1 rs743572 might modify the risk of PCa in the individuals of African origin [31].
In the present study only rs6163 showed a significant association with a shorter BRFS in patients with prostatectomy. There is not much data about genetic markers in androgen metabolism as indicators of outcomes after PCa therapies and scarce data in relation to CYP17A1 gene [7, 20]. Wright et al. investigated the genetic association between three SNPs in CYP17A1 (rs743572, rs10883783 and rs17115100) and their responses to treatment (principally prostatectomy), but they did not find any differences in the risk of recurrence/progression by this genotype analysis [32]. However SNPs in other genes such as SRD5A and HSD17B, both involved in androgen metabolism, were shown to be associated to biochemical recurrence in Caucasian and Asian PCa patients after prostatectomy [22, 33]. rs1004467 (CYP17A1) was previously found to be associated with PCa risk and disease progression after ADT among Japanese population [30]. It is proved that serum dihidrotestosterone (DHT) level and testosterone was significantly elevated in GG genotype for rs1004467 (CYP17A1) compared to A allele carriers [34]. Studies have shown that those individuals with higher levels of DHT were susceptible to CRPC development in the future. This could provide insights for an early PCa detection, diagnosis, management and potential therapeutic targets [34]. Last, our LD analysis showed a strong linkage between rs743572, rs6162 and rs6163, (CYP17A1). Similar results were supported by other studies proving that both, rs6162 and rs6163 were strongly linked with rs743572, so these SNPs exhibit their functions in coordination with rs743572 as a risk genotype [30].
In relation to ESR genes, we demonstrated a statistical association between rs1256049 (ESR2) and a higher D´Amico risk. Several studies have focused on the effect of ESR polymorphisms (like rs9340799 and rs1256049) on PCa development, but none of them have evaluated their relation to PCa aggressiveness. rs1256049 is a silent mutation in codon 328 and G>A change has a direct effect on modifying the secondary structure of the mRNA; leading to changes in mRNA stability and translation which makes it as a candidate polymorphism for PCa susceptibility in Caucasians [35, 36]. This SNP has been previously associated to an increased PCa risk, both in codominant and recessive genetic models [35].
Despite, ESR1 (rs1062577) was not associated with any clinical variable in the present study, it was previously related to changes in plasma steroid levels conducting cancer aggressiveness and efficacy of ADT [7]. We have proved with an in silico analysis that A allele (rs1062577) generates a new miRNA binding site for hsa-miR-3646, hsa-miR-3662 and hsa-miR-5585-3p [37]. In PCa there are similar expression patterns for miR-3646 and miR-3662, similarly happens in other hormonal dependant tumors like breast, ovarian or uterine [37]. For example, in breast cancer, miR-3646 produces cellular proliferation [38]. Moreover A allele in rs1062577 (ESR1) produces a loss of miR-186 and miR-6507-5p binding site [37]. miR-186 has a role as a PCa tumor suppressor, so any alteration on this SNP will be related to an oncological event [39].
We demonstrated in LHCGR gene, a relation between rs2293275 and high T stages. Other researchers had found an association between PCa incidence and this variant in different racial groups by no for stage of the disease [40]. rs2293275 produces a missense mutation with an amino acid change (N312S) [29]. Position 312 is located in exon 10 and this is important for receptor activation [40], but this amino acid change seems to not affect frameshift region so it has a tolerated phenotypic effect (SIFT score: 0.774) [29].
Finally, we did not find any significant association with PCa for AR SNPs (rs6152 and rs9332969). Likewise, a previous study in a Caucasian population did not find differences between case-control for rs6152 but it was described an association for this polymorphism with advances stages [41]. rs6152 is a coding-synonymous region (E213E) not frameshift alteration and with a tolerated phenotypic effect [29] (SITF score: 1). In relation to rs9332969, despite we not report any allelic variation (100% of patients and controls are G carriers), this SNP seems to be an attractive candidate biomarker for PCa aggressiveness because is a missense variant with an amino acid change (R841H) [29]. R841H variation causes androgen insensitivity syndrome [42]. This is due because AR protein with R841H variation alters the interaction with androgens which result in a partial AR functional disruption [42] and its frequency and effect had not been previously studied in the European population [43].
We fully acknowledge that our results must be interpreted with caution, as the sample size is limited. One core limitation of the study is the low specificity of prostate biopsy; although we tried to control this issue by performing systematic 20-core biopsies, the number of false negatives could still be high. Moreover, the observation period was relatively short (median follow-up of 34.18 month) and only one-third of the patients had a high risk of recurrence; and the number of patient in each treatment group is reduced. Nevertheless, this is the first time that the role of AR, CYP17A1, LHCGR and ESR polymorphisms have been studied in relation to PCa aggressiveness.
Conclusion
We describe the initial roles of CYP17A1, AR, LHCGR and ESR as risk disease biomarkers. We believe that rs743572, rs6162, rs6163, and rs1256049 (CYP17A1); and rs2293275 (LHCGR) are promise biomarkers for PCa aggressiveness. Future studies seem warranted in order to evaluate the real predictive and prognostic impact of CYP17A1, AR, LHCGR and ESR variations on each specific treatment response in a larger cohort of patients within longer observation periods.
Supporting information
S1 Fig. Association between rs1256049 (ESR2) and Gleason score.
https://doi.org/10.1371/journal.pone.0185447.s001
(TIF)
S1 Table. Description of SNPs included in the study.
Information was obtained from Ensembl and dbSNP data base, National Center for Biotechnology Information (NCBI).
https://doi.org/10.1371/journal.pone.0185447.s002
(XLSX)
S3 Table. Significant linkage disequilibrium for variants.
Test of linkage disequilibrium for all pairs of loci performed with SNPstats software. In blue, all polymorphism that are out of balance (significance level = 0.05). p = p-value; r2 = standardized disequilibrium values.
https://doi.org/10.1371/journal.pone.0185447.s004
(XLSX)
S4 Table. Multigene analysis.
All genotypes with a frequency less than 5% in the population have been omitted.
https://doi.org/10.1371/journal.pone.0185447.s005
(XLSX)
Acknowledgments
We want to thank to Gema Garcia-Rodriguez and Manuela Expósito-Ruiz for their collaboration and all the donors and nursery that make this study possible.
This paper will be part of Inmaculada Robles’s Doctoral Thesis performed in the Biomedicine Doctorate Program of the University of Granada, Spain
References
- 1. Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, et al. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer. 2015;136(5):E359–86. Epub 2014/09/16. pmid:25220842.
- 2. Narita S, Mitsuzuka K, Tsuchiya N, Koie T, Kawamura S, Ohyama C, et al. Reassessment of the risk factors for biochemical recurrence in D'Amico intermediate-risk prostate cancer treated using radical prostatectomy. Int J Urol. 2015;22(11):1029–35. Epub 2015/08/21. pmid:26290306.
- 3. Irshad S, Bansal M, Castillo-Martin M, Zheng T, Aytes A, Wenske S, et al. A molecular signature predictive of indolent prostate cancer. Sci Transl Med. 2013;5(202):202ra122. pmid:24027026; PubMed Central PMCID: PMCPMC3943244.
- 4. Lalonde E, Ishkanian AS, Sykes J, Fraser M, Ross-Adams H, Erho N, et al. Tumour genomic and microenvironmental heterogeneity for integrated prediction of 5-year biochemical recurrence of prostate cancer: a retrospective cohort study. Lancet Oncol. 2014;15(13):1521–32. pmid:25456371.
- 5. Ramos-Montoya A, Lamb AD, Russell R, Carroll T, Jurmeister S, Galeano-Dalmau N, et al. HES6 drives a critical AR transcriptional programme to induce castration-resistant prostate cancer through activation of an E2F1-mediated cell cycle network. EMBO Mol Med. 2014;6(5):651–61. pmid:24737870; PubMed Central PMCID: PMCPMC4023887.
- 6. Stelloo S, Nevedomskaya E, van der Poel HG, de Jong J, van Leenders GJ, Jenster G, et al. Androgen receptor profiling predicts prostate cancer outcome. EMBO Mol Med. 2015;7(11):1450–64. pmid:26412853; PubMed Central PMCID: PMCPMC4644377.
- 7. Lévesque É, Huang SP, Audet-Walsh É, Lacombe L, Bao BY, Fradet Y, et al. Molecular markers in key steroidogenic pathways, circulating steroid levels, and prostate cancer progression. Clin Cancer Res. 2013;19(3):699–709. Epub 2012/11/27. pmid:23186779.
- 8. Alukal JP, Lepor H. Testosterone Deficiency and the Prostate. Urol Clin North Am. 2016;43(2):203–8. pmid:27132577.
- 9. Lau KM, To KF. Importance of Estrogenic Signaling and Its Mediated Receptors in Prostate Cancer. Int J Mol Sci. 2016;17(9). Epub 2016/08/31. pmid:27589731; PubMed Central PMCID: PMCPMC5037713.
- 10. Weng H, Li S, Huang JY, He ZQ, Meng XY, Cao Y, et al. Androgen receptor gene polymorphisms and risk of prostate cancer: a meta-analysis. Sci Rep. 2017;7:40554. Epub 2017/01/16. pmid:28091563; PubMed Central PMCID: PMCPMC5238402.
- 11. El Ezzi AA, Baker MT, Zaidan WR, Hraiki KM, El Saidi MA, Kuddus RH. Association of Polymorphisms in the VDR, CYP17 and SRD5A2 Genes and Prostate Cancer Among Lebanese Men. Asian Pac J Cancer Prev. 2017;18(1):93–100. Epub 2017/01/01. pmid:28240015.
- 12. Beuten J, Gelfond JA, Franke JL, Weldon KS, Crandall AC, Johnson-Pais TL, et al. Single and multigenic analysis of the association between variants in 12 steroid hormone metabolism genes and risk of prostate cancer. Cancer Epidemiol Biomarkers Prev. 2009;18(6):1869–80. pmid:19505920.
- 13.
Information NCfB. NCBI. http://www.ncbi.nlm.nih.gov /.
- 14. Solé X, Guinó E, Valls J, Iniesta R, Moreno V. SNPStats: a web tool for the analysis of association studies. Bioinformatics. 2006;22(15):1928–9. Epub 2006/05/23. pmid:16720584.
- 15. Barbieri CE, Bangma CH, Bjartell A, Catto JW, Culig Z, Gronberg H, et al. The mutational landscape of prostate cancer. Eur Urol. 2013;64(4):567–76. Epub 2013/06/14. pmid:23759327; PubMed Central PMCID: PMC4342117.
- 16. Colicchia M, Morlacco A, Cheville JC, Karnes RJ. Genomic tests to guide prostate cancer management following diagnosis. Expert Rev Mol Diagn. 2017;17(4):367–77. Epub 2017/03/13. pmid:28277880.
- 17. Helfand BT, Catalona WJ, Xu J. A genetic-based approach to personalized prostate cancer screening and treatment. Curr Opin Urol. 2015;25(1):53–8. pmid:25405931; PubMed Central PMCID: PMCPMC4281884.
- 18. Alvarez-Cubero MJ, Pascual-Geler M, Martinez-Gonzalez LJ, Expósito Ruiz M, Saiz M, Cozar JM, et al. Association between RNASEL, MSR1, and ELAC2 single nucleotide polymorphisms and gene expression in prostate cancer risk. Urol Oncol. 2016;34(10):431.e1–8. Epub 2016/06/16. pmid:27318894.
- 19. Maia S, Cardoso M, Paulo P, Pinheiro M, Pinto P, Santos C, et al. The role of germline mutations in the BRCA1/2 and mismatch repair genes in men ascertained for early-onset and/or familial prostate cancer. Fam Cancer. 2016;15(1):111–21. Epub 2015/08/21. pmid:26289772.
- 20. Van den Broeck T, Joniau S, Clinckemalie L, Helsen C, Prekovic S, Spans L, et al. The role of single nucleotide polymorphisms in predicting prostate cancer risk and therapeutic decision making. Biomed Res Int. 2014;2014:627510. Epub 2014/02/19. pmid:24701578; PubMed Central PMCID: PMCPMC3950427.
- 21. Fujimoto N, Shiota M, Tomisaki I, Minato A. Gene Polymorphism-related Individual and Interracial Differences in the Outcomes of Androgen Deprivation Therapy for Prostate Cancer. Clin Genitourin Cancer. 2017;15(3):337–42. Epub 2017/02/12. pmid:28188049.
- 22. Audet-Walsh É, Bellemare J, Lacombe L, Fradet Y, Fradet V, Douville P, et al. The impact of germline genetic variations in hydroxysteroid (17-beta) dehydrogenases on prostate cancer outcomes after prostatectomy. Eur Urol. 2012;62(1):88–96. Epub 2011/12/21. pmid:22209174.
- 23. Xu J, Zheng SL, Isaacs SD, Wiley KE, Wiklund F, Sun J, et al. Inherited genetic variant predisposes to aggressive but not indolent prostate cancer. Proc Natl Acad Sci U S A. 2010;107(5):2136–40. Epub 2010/01/19. pmid:20080650; PubMed Central PMCID: PMC2836698.
- 24.
Xu J, Isaacs WB, Grönberg H. Methods and compositions for correlating genetic markers with risk of aggressive prostate cancer. Google Patents; 2017.
- 25. Oh JJ, Park S, Lee SE, Hong SK, Lee S, Choe G, et al. The use of exome genotyping to predict pathological Gleason score upgrade after radical prostatectomy in low-risk prostate cancer patients. PLoS One. 2014;9(8):e104146. Epub 2014/08/05. pmid:25093842; PubMed Central PMCID: PMCPMC4122442.
- 26. Hart SN, Ellingson MS, Schahl K, Vedell PT, Carlson RE, Sinnwell JP, et al. Determining the frequency of pathogenic germline variants from exome sequencing in patients with castrate-resistant prostate cancer. BMJ Open. 2016;6(4):e010332. Epub 2016/04/17. pmid:27084275; PubMed Central PMCID: PMC4838679.
- 27. Udhane SS, Dick B, Hu Q, Hartmann RW, Pandey AV. Specificity of anti-prostate cancer CYP17A1 inhibitors on androgen biosynthesis. Biochem Biophys Res Commun. 2016;477(4):1005–10. Epub 2016/07/06. pmid:27395338.
- 28. Brureau L, Moningo D, Emeville E, Ferdinand S, Punga A, Lufuma S, et al. Polymorphisms of Estrogen Metabolism-Related Genes and Prostate Cancer Risk in Two Populations of African Ancestry. PLoS One. 2016;11(4):e0153609. Epub 2016/04/13. pmid:27074016; PubMed Central PMCID: PMCPMC4830606.
- 29. Sim NL, Kumar P, Hu J, Henikoff S, Schneider G, Ng PC. SIFT web server: predicting effects of amino acid substitutions on proteins. Nucleic Acids Res. 2012;40(Web Server issue):W452–7. Epub 2012/06/11. pmid:22689647; PubMed Central PMCID: PMCPMC3394338.
- 30. Yamada T, Nakayama M, Shimizu T, Nonen S, Nakai Y, Nishimura K, et al. Genetic polymorphisms of CYP17A1 in steroidogenesis pathway are associated with risk of progression to castration-resistant prostate cancer in Japanese men receiving androgen deprivation therapy. Int J Clin Oncol. 2013;18(4):711–7. Epub 2012/06/21. pmid:22714708.
- 31. Wang Y, Zhang Y, Meng H, Hou X, Li Z, Liu Q, et al. Quantitative assessment of the association between CYP17 rs743572 polymorphism and prostate cancer risk. Cell Biochem Biophys. 2015;71(2):983–91. Epub 2014/10/18. pmid:25323563.
- 32. Wright JL, Kwon EM, Lin DW, Kolb S, Koopmeiners JS, Feng Z, et al. CYP17 polymorphisms and prostate cancer outcomes. Prostate. 2010;70(10):1094–101. pmid:20503394; PubMed Central PMCID: PMCPMC2878282.
- 33. Audet-Walsh E, Bellemare J, Nadeau G, Lacombe L, Fradet Y, Fradet V, et al. SRD5A polymorphisms and biochemical failure after radical prostatectomy. Eur Urol. 2011;60(6):1226–34. Epub 2011/06/29. pmid:21715084.
- 34. Poniah P, Mohamed Z, Apalasamy YD, Mohd Zain S, Kuppusamy S, Razack AH. Genetic polymorphisms in the androgen metabolism pathway and risk of prostate cancer in low incidence Malaysian ethnic groups. Int J Clin Exp Med. 2015;8(10):19232–40. Epub 2015/10/15. pmid:26770559; PubMed Central PMCID: PMCPMC4694459.
- 35. Dai ZJ, Wang BF, Ma YF, Kang HF, Diao Y, Zhao Y, et al. Current evidence on the relationship between rs1256049 polymorphism in estrogen receptor-β gene and cancer risk. Int J Clin Exp Med. 2014;7(12):5031–40. Epub 2014/12/15. pmid:25664002; PubMed Central PMCID: PMCPMC4307449.
- 36. Fu C, Dong WQ, Wang A, Qiu G. The influence of ESR1 rs9340799 and ESR2 rs1256049 polymorphisms on prostate cancer risk. Tumour Biol. 2014;35(8):8319–28. Epub 2014/05/24. pmid:24859835.
- 37. Gong J, Tong Y, Zhang HM, Wang K, Hu T, Shan G, et al. Genome-wide identification of SNPs in microRNA genes and the SNP effects on microRNA target binding and biogenesis. Hum Mutat. 2012;33(1):254–63. Epub 2011/11/23. pmid:22045659.
- 38. Tao S, Liu YB, Zhou ZW, Lian B, Li H, Li JP, et al. miR-3646 promotes cell proliferation, migration, and invasion via regulating G2/M transition in human breast cancer cells. Am J Transl Res. 2016;8(4):1659–77. Epub 2016/04/15. pmid:27186291; PubMed Central PMCID: PMCPMC4859896.
- 39. Hua X, Xiao Y, Pan W, Li M, Huang X, Liao Z, et al. miR-186 inhibits cell proliferation of prostate cancer by targeting GOLPH3. Am J Cancer Res. 2016;6(8):1650–60. Epub 2016/08/01. pmid:27648356; PubMed Central PMCID: PMCPMC5004070.
- 40. Ingles SA, Liu SV, Pinski J. LHRH and LHR genotypes and prostate cancer incidence and survival. Int J Mol Epidemiol Genet. 2013;4(4):228–34. Epub 2013/11/28. pmid:24319538; PubMed Central PMCID: PMCPMC3852642.
- 41. Hayes VM, Severi G, Eggleton SA, Padilla EJ, Southey MC, Sutherland RL, et al. The E211 G>A androgen receptor polymorphism is associated with a decreased risk of metastatic prostate cancer and androgenetic alopecia. Cancer Epidemiol Biomarkers Prev. 2005;14(4):993–6. Epub 2005/04/13. pmid:15824176.
- 42. Disciglio V, Devecchi A, Palumbo O, Carella M, Penso D, Milione M, et al. Whole exome sequencing and single nucleotide polymorphism array analyses to identify germline alterations in genes associated with testosterone metabolism in a patient with androgen insensitivity syndrome and early-onset colorectal cancer. Chin J Cancer. 2016;35(1):51. Epub 2016/06/07. pmid:27267075; PubMed Central PMCID: PMCPMC4897824.
- 43.
NCBI. db SNP. https://www.ncbi.nlm.nih.gov/projects/SNP/snp_ss.cgi?ss=ss12675702.