Impact Factor 4.848 | CiteScore 3.5
More on impact ›

Original Research ARTICLE

Front. Oncol., 22 October 2020 | https://doi.org/10.3389/fonc.2020.583625

Evaluation of Polygenic Risk Scores for Prediction of Prostate Cancer in Korean Men

Jong Jin Oh1,2†, Eunae Kim3†, Eunjin Woo3, Sang Hun Song1, Jung Kwon Kim1, Hakmin Lee1, Sangchul Lee1, Sung Kyu Hong1,2 and Seok-Soo Byun1,2*
  • 1Department of Urology, Seoul National University Bundang Hospital, Seongnam, South Korea
  • 2Department of Urology, Seoul National University College of Medicine, Seoul, South Korea
  • 3Procagen, Seongnam, South Korea

Aims: The purpose of this study is to evaluate an aggregate influence of prostate cancer (PCa) susceptibility variants on the development of PCa in Korean men by using the polygenic risk score (PRS) approach.

Methods: An analysis of 1,001 cases of PCa and 2,641 controls was performed to: (i) identify potential PCa-related risk loci in Koreans and (ii) validate the cumulative association between these loci and PCa using the PRS. Subgroup analyses based on risk stratification were conducted to better characterize the potential correlation to key PCa-related clinical outcomes (e.g., Gleason score, prostate-specific antigen levels). The results were replicated using 514 cases of PCa and 548 controls from an independent cohort.

Results: Genome-wide association analysis from our discovery cohort revealed 11 candidate single-nucleotide polymorphisms (SNPs) associated with PCa showing statistical significance of p < 5.0 × 10–5. Seven variants were located at 8q24.21 (rs1016343, rs16901979, and rs13252298 in PRNCR1; rs4242384, rs7837688, and rs1447295 in CASC8; and rs1512268 in NKX3). Two variants located within HNF1B (rs7501939 and rs4430796) had a significant negative association with PCa risk [odds ratio (OR) = 0.717 and 0.747, p = 6.42 × 10–7 and 3.67 × 10–6, respectively]. Of the six independent SNPs that remained after linkage disequilibrium (LD) pruning, the top four SNPs best predicted PCa risk with an area under the receiver operating characteristic curve (AUC) of 0.637 (95% CI: 0.582–0.692). Those with top 25% polygenic risk had a 4.2-fold increased risk of developing PCa compared with those with low risk.

Conclusion: Eleven PCa risk variants in Korean men were identified; PRSs of a subset of these variants could help predict PCa susceptibility.

Introduction

Prostate cancer (PCa) is the second most common cancer in men worldwide (1). In South Korea, the number of men diagnosed with PCa increased rapidly between 1999 and 2009, with an annual increase of 12.9% (2). Between 2009 and 2017, the incidence of PCa grew by 0.8% per year, resulting in the fourth highest incidence and the third highest prevalence of cancer types in South Korean men. Interestingly, while the incidence of PCa continues to increase in South Korea, it has gradually declined for other common cancers (e.g., stomach, lung, colon, liver, and thyroid) (2). Importantly, despite having a lower incidence compared with Western populations, the proportion of advanced-stage PCa in Koreans and other Asian populations is higher (3). Given the associated public health burden caused by PCa, there is a growing need to identify high-risk groups to generate effective screening and prevention strategies for PCa.

Genetic profiling can be a useful clinical instrument to help determine an individual’s risk for PCa. Multiple large-scale genome-wide association studies (GWASs) have led to the identification of more than 170 single-nucleotide polymorphisms (SNPs) underlying susceptibility to PCa (4). Multiple PCa-associated SNPs have been estimated to explain 33% of the risk of developing PCa (5). SNPs associated with PCa can vary greatly by population; several attempts have focused on identifying SNPs that may be associated with PCa in Asian populations. One large-scale meta-analysis reported Asian-specific PCa-associated SNPs from two Asian (Japanese and Chinese) populations (6). Additionally, a Korean population-based exome-wide study identified five significant SNPs across four distinct loci (7). Each of the common SNPs identified by the GWAS confers small-to-modest effects on the development of PCa (8).

An aggregate influence of SNPs can be assessed by generating a polygenic risk score (PRS)—a measure of the cumulative contribution of individual SNPs carried by a particular person (9). Even if individual variants have only small effect sizes, their cumulative impact on risk of PCa can be significant, thus making the PRS a potentially powerful tool for the prediction of PCa (9, 10). Eeles et al. (11) generated a PRS using 68 established PCa risk variants and reported that men in the top 1% of the risk distribution had more than a four-fold increased risk for PCa compared with those in the average risk distribution range. Another PRS study demonstrated that men in the top 10% of the risk distribution had a 3.19-fold higher risk of PCa compared with those with average risk (12). In this study, we assessed the cumulative impact of PCa-related genetic variants in predicting the risk of PCa using weighted PRS in the Korean male population.

Materials and Methods

Ethics Statement

After approval by our institutional review board (B-1312/232-302), all analyses were performed following the Declaration of Helsinki. All study participants provided written informed consent.

Study Population and Genome-Wide Association Study Genotyping

For the discovery of PCa-associated candidate SNPs, we initially obtained genotype data of 1,001 PCa samples from a single tertiary hospital and 2,210 controls from the Korean Association Resource (KARE) study as part of the Korean Genome and Epidemiology Study (KoGES). Gleason score (GS) was identified with ≥12 core transrectal ultrasound prostate biopsy and/or radical prostatectomy (RP) specimens; specimens were reviewed by a single experienced uro-pathologist. Controls were: (i) men who had never been diagnosed with PCa, (ii) residing in the cities of Ansung or Ansan, (iii) recruited between 2001 and 2002, and (iv) 60 years of age or older. More detailed information about the cohort is available in a previously published article (13).

Genotyping of PCa samples was performed with the HumanExome BeadChip 12v1-1 system (Illumina, Inc., San Diego, CA, United States) (14, 15). Details about SNP content and selection strategies can be found at http://genome.sph.umich.edu/wiki/Exome_Chip_Design. Genotype calling was performed using Illumina’s GenTrain version 2.0 clustering algorithm with the GenomeStudio software (V2011.1). Cluster boundaries were determined using Illumina’s standard cluster file. To improve the accuracy of variant calling, manual reclustering and visual inspection were conducted for genotypes based on the CHARGE clustering method (14). Sample quality control was carried out to exclude samples with genotyping rates < 95%, heterozygosity, and cryptic relatedness. Markers were excluded based on the following criteria: (i) monomorphic in our samples, (ii) with missing call rate > 5%, (iii) with minor allele frequency (MAF) < 5%, or (iv) significantly deviated from the Hardy–Weinberg equilibrium (p < 1.0 × 10–6) using PLINK 1.9. After quality control, 24,023 variants from 984 PCa cases and 2,194 control subjects remained for subsequent analysis. To evaluate allelic associations with PCa development, logistic regression analysis was performed given case/control status after adjusting for age as a covariate using the PLINK software.

For the evaluation of polygenic risk in an independent cohort, 516 cases of PCa and 546 controls from Chungbuk National University Hospital were initially obtained. Candidate PCa-associated SNPs showing suggestive significance were included for analyses. The genotyping of these SNPs was performed using the Fluidigm 192.24 Dynamic Array TM IFC and Biomark HD systems. Duplicates and negative controls were included in each 96-well plate for quality control. The average concordance rate between the duplicate samples was >99%. After sample quality control and exclusion of control samples <60 years of age, we included 514 cases and 124 controls for further analysis.

For PRS construction, the most significant SNP in each linkage disequilibrium (LD) block was further selected via LD clumping (R2 < 0.001). The PRS was calculated by summing effect size-weighted counts of risk alleles for PCa-associated SNPs. An individual i’s PRS is defined as follows:

S i = ( j = 1 ) M X i j β ^ j n ,

where Xj is the number of risk alleles (0, 1, or 2) for the variant j, βj, is weighting [log(OR)] of the variant j, which is obtained from the discovery set, and n is the total number of the variants included. For those SNPs whose minor alleles showed protective effects on PCa, we converted their minor alleles to major alleles as risk alleles, which results in positive weight values for all variants. As an increasing number of top SNPs were included (n = 3–6), predictive abilities of their PRS were compared using the area under the receiver operating characteristic (ROC) curve (AUC) (16). Improvement in AUC between ROC curves were tested using Delong’s method (17). Youden’s Index (J), which is defined as (sensitivity + specificity-1), was acquired to capture predictive performance and also to determine the cutoff PRS at the maximum J (18). The analyses were performed using the R package “pROC.”

Results

For cases in the discovery set, the median age was 68 years and the median prostate-specific antigen (PSA) level was 9.19 ng/ml. Most of the cases (90%) were diagnosed with PCa with a GS of 7 or higher. The controls were 3 years younger (median age of 65 years) compared to the case group (p < 0.0001; Table 1). The average body mass index (BMI) was significantly higher in cases compared to controls by 1.1 kg/m2 (p < 0.0001). For the validation set, the median ages of cases and controls are 69 and 67 years (p = 0.04). Eighty percent of the cases in the validation set are diagnosed with PCa with a GS ≥ 7.

TABLE 1
www.frontiersin.org

Table 1. Baseline characteristics of the study population.

Genome-wide association analysis from the discovery cohort revealed 11 SNPs associated with PCa, showing genome-wide suggestive significance of p < 5.0 × 10–5 (Table 2). Seven variants were located at 8q24.21 (rs1016343, rs16901979, and rs13252298 in PRNCR1; rs4242384, rs7837688, and rs1447295 in CASC8; rs1512268 in NKX3). The most significantly associated SNP was rs1016343 [odds ratio (OR) = 1.598 (1.424–1.793); p = 1.46 × 10–15] located within the PRNCR1 locus. Two variants located within HNF1B (rs7501939 and rs4430796) had a significant negative association with PCa risk (OR = 0.717 and 0.746, p = 6.42 × 10–7, and 3.67 × 10–6, respectively). Of the top 11 SNPs, six variants were positively associated with PCa risk (OR > 1), while five showed negative associations (OR < 1). Of the three variants in PRNCR1, one variant (rs13252298) showed protective effects compared to two other susceptibility variants (rs1016343 and rs16901979).

TABLE 2
www.frontiersin.org

Table 2. PCa-associated SNPs (genome-wide significance <5.0 × 10–5) identified by GWAS (n = 10).

After LD clumping of the top 11 SNPs, the remaining six variants (Table 3) were included for PRS calculation. The mean PRS was significantly higher in PCa cases compared to controls across all SNP sets (Table 4 and Figure 1A). Polygenic risk of PCa was best predicted with an AUC of 0.637 (95% CI: 0.582–0.692) when calculated upon the top four independent SNPs (Figure 2A). The PRS model built with four SNPs showed significant improvement in AUC compared to the model with three SNPs according to DeLong’s tests (p = 0.005; Table 4 and Figure 2B). However, it was not significantly different from the AUCs of PRS models computed upon five or six SNPs. Meanwhile, the set of the top four SNPs yielded the second highest maximum Youden’s index (J = 0.221), with sensitivity and specificity of 0.543 and 0.677, respectively, which was slightly lower than that of the top five SNPs (J = 0.227; Table 4). Overall, the specificity of the PRS at the maximum Youden’s index was higher compared to its sensitivity, with the exception of the model built with three SNPs (sensitivity: 0.607, specificity: 0.565).

TABLE 3
www.frontiersin.org

Table 3. Candidate PCa-associated SNPs retained after LD clumping for construction of PRS (n = 6).

TABLE 4
www.frontiersin.org

Table 4. Comparison of the predictive performance of PRS according to different SNP sets.

FIGURE 1
www.frontiersin.org

Figure 1. Comparison of prostate cancer (PCa) polygenic risk between PCa cases and controls. (A) Distribution of the polygenic risk score (PRS) in PCa cases and controls (black solid line represents cutoff value at the maximum J). (B) Distribution of the PRS in PCa cases and controls (each dashed line represents high-risk PRS group—top 2.5, 5, 10, and 20%). (C) Distribution of cases and controls according to PRS quantiles. (D) Odds ratio for developing PCa according to PRS quantiles.

FIGURE 2
www.frontiersin.org

Figure 2. Receiver operating characteristic (ROC) curve of the polygenic risk score (PRS) for the prediction of prostate cancer (PCa). (A) Comparison of ROC curves according to increasing number of single-nucleotide polymorphisms (SNPs) included for PRS calculation (n = 3–6). (B) ROC curves of the PRS using four SNPs.

When the PRS was computed upon the top four independent SNPs, the upper quantile (top 25%) of the PRS had a 4.2-fold greater risk of developing PCa compared with the reference group (Q1) [OR: 4.2 (2.32–7.98)] (Table 5 and Figure 1C). With a larger number of PCa cases observed in the upper quantiles of the PRS, the Q2–Q4 groups carried significantly higher risk compared to the reference (Q1) group (Table 5 and Figure 1D). High-risk group defined by top 5% to 25% in the distribution had also significantly elevated risk of PCa compared to the remaining population: the top 10% of the PRS had a 3.08-fold risk, and the top 5% had a 3.71-fold risk of developing PCa compared to the remaining population (Table 6 and Figure 1B).

TABLE 5
www.frontiersin.org

Table 5. Distribution of cases and controls according to PRS* quintiles.

TABLE 6
www.frontiersin.org

Table 6. Risk of high PRS groups for development of PCa.

Discussion

In the development of PCa, genetic susceptibility plays an important role. Stratification of individuals based on their inherited genetic risk can be important in screening and prevention strategies of PCa. The present study used multiple PCa-associated SNPs detected in Korean men to evaluate their predictive ability using the weighted PRS approach.

We identified 11 SNPs associated with PCa showing statistical significance of p < 5.0 × 10–5. Since the ultimate goal of this study was to evaluate the utility of the PRS as a predictor of PCa risk rather than identifying causal PCa SNPs, we applied a lenient statistical threshold to select candidate SNPs for PRS calculation and thus compared the performance by different PRS models. Those PCa-associated variants included rs1016343, rs16901979, and rs13252298 located at 8q24.21 within PRNCR1. The association of rs1016343, which was most significantly associated in this study, has been previously reported in other populations (19, 20); the effect of the variant on PCa in the Korean population (OR = 1.598) were greater compared to that of other populations of European ancestry. Four other variants (rs7837688, rs4242384, rs4242382, and rs1447295) in CASC8 were also replicated; the effect sizes were comparable to those of other studies. The significant association of the rs1447295 variant has been reported in Japanese and Chinese populations (21, 22). The results presented here are of importance in that they focus specifically on the Korean population.

We identified several variants protective for PCa risk. One variant (rs13252298) located within PRNCR1 showed protective effects, while the other two PRNCR1 loci conferred susceptibility to PCa (rs1016343 and rs16901979). It was contrasted with other regions harboring multiple significant SNPs such as CASC9, HNF1B, and RFX6, within which consistent directions of effects were observed. Two variants (rs7501939 and rs4430796) at 17q12 in HNF1B showed protective effects on PCa risk, with ORs of 0.717 (0.629–0.817) and 0.746 (0.659–0.844), respectively. Previous studies also reported protective effects for these variants in the European (23) and Korean population (24). However, the associations reported by Gudmundsson et al. (25) were not consistent with this and other studies of European populations, showing susceptibility to PCa with ORs of 1.19 (1.12–1.26) and 1.22 (1.15–1.30) for rs7501939 and rs4430796, respectively. The HNF1B gene (formerly known as transcription factor TCF2) is known to, at least in part, regulate the levels of metabolic and hormonal factors in PCa.

Since the previous report on strong cumulative effects of five SNPs on PCa in Korean men (10), many studies have reported significant associations between the PRS and the risk of PCa (19, 26, 27). In the present study, we applied weighted PRS models after excluding controls <60 years of age from the previously available cohort to reduce possible confounding effects due to ages and adjusted for ages as a covariate in analyses. We also compared polygenic risk by different weighted PRS models, which produced an AUC of 0.637 using the top four independent SNPs compared to the AUC of 0.605 using non-weighted PRS using the top five SNPs in the previous study (10). Despite the improvement in predictive performance, our model shows modest performance compared to other large-scale studies such as an international PCa genetics consortium [Prostate Cancer Association Group to Investigate Cancer Associated Alterations in the Genome (PRACTICAL)] that yielded an AUC of 0.67 using 68 genetic variants (5). Still, this study holds significant implications as an efficient tool for screening high-risk group using only a small number of SNPs for PRS calculation: those with the top 25% PRS showed a 4.2-fold increase of developing PCa compared to the low-risk group. It was comparable impact of the PRS built with 68 variants reported by Eeles et al. (11), in which men in the top 1% of the risk distribution had a 4.4-fold increased polygenic risk for PCa compared with those with average risk. Latino men in the top 10% PRS stratum had a 3.19-fold elevated risk compared with those in the average range (12). Despite multiple SNPs commonly found across populations, such discrepancy in genetic risk between populations arises from variations in risk allele frequencies, LD structures, and effect sizes on PCa.

Individually calculated PRSs can facilitate stratification of disease risk, which can be utilized for screening and prevention in clinical practice (8, 9). The conventional PRS approach has been recently extended to polygenic hazard score (PHS), which informs the onset age and age-specific genetic risk of certain diseases (28). Applications of PHS may add valuable information for personalized life planning and disease screening. Furthermore, many previous studies showed that, despite genetics being a non-modifiable risk factor, those at higher polygenic risk have been shown to have reduced risk by lifestyle modifications or clinical interventions in other diseases (2932). To introduce the PRS into clinical practice, there is a need to improve its predictive value by (i) obtaining more valid summary statistics to be applied for computing the PRS using larger-scale data sets, (ii) applying advanced modeling technologies (such as machine learning and other artificial intelligence methods), (iii) incorporating family history, lifestyle factors, and other clinical factors into modeling, and (iv) including super healthy controls strictly screened for diseases.

This study provides valuable scientific evidence to develop screening and prevention strategies that can identify genetic high-risk groups in Korean men. Nonetheless, we acknowledge that our study had several limitations. First, our sample size is modest compared with other large cohort-based studies, especially those on European populations. Another large-scale study may be necessary for deriving more valid summary statistics used for PRS calculation and further testing its performance. Second, although only a small number of variants were able to screen high-risk PCa groups, the inclusion of even more SNPs may achieve a superior predictive ability. Finally, the validation set is not representative of the Korean male population, as the PCa prevalence is 84% in the validation set, which is considerably higher than that of the Korean male population. Although we defined the lowest PRS quantile group as the reference group to represent the general average-risk population, it could have led to underestimation of polygenic risk in high-risk groups due to the possibly elevated PRS in the reference group. For evaluation of predictive utility and generalizability of our findings, an external validation set representative of the general Korean male population may be warranted.

In conclusion, we identified 11 PCa risk variants in Korean men and report that PRSs using a subset of these variants may be useful for determining an individual’s risk of developing PCa. The addition of individually calculated PRSs effectively increased the accuracy of predicting PCa. Future studies on modifications of polygenic risk by lifestyle factors may add valuable scientific evidence in preventing the development of PCa in which genetics plays a critical role.

Data Availability Statement

The datasets generated for this study are available upon request to the corresponding author.

Ethics Statement

The studies involving human participants were reviewed and approved by Seoul National University Bundang Hospital institutional review board. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

JO, EK, and S-SB conceived, designed, and supervised the experiments. SS, JK, HL, SL, and SH performed the experiments. JO, EK, and EW analyzed the data. JO and EK wrote the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by Grant No. 13-2016-007, 13-2017-006, 14-2020-030, and 14-2020-032 from the Seoul National University Bundang Hospital (SNUBH) Research Fund and Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (2020R1A2C1100011).

Conflict of Interest

EK and EW was employed by the company Procagen, South Korea.

The remaining 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.

References

1. Rawla P. Epidemiology of prostate cancer. World J Oncol. (2019) 10:63–89.

Google Scholar

2. National Cancer Center Annual Report of Korean Central Cancer Registry 2017. Goyang-si: National Cancer Center (2020).

Google Scholar

3. Ito K. Prostate cancer in Asian men. Nat Rev Urol. (2014) 11:197–212. doi: 10.1038/nrurol.2014.42

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Schumacher FR, Al Olama AA, Berndt SI, Benlloch S, Ahmed M, Saunders EJ, et al. Association analyses of more than 140,000 men identify 63 new prostate cancer susceptibility loci. Nat Genet. (2018) 50:928–36.

Google Scholar

5. Szulkin R, Whitington T, Eklund M, Aly M, Eeles RA, Easton D, et al. Prediction of individual genetic risk to prostate cancer using a polygenic score. Prostate. (2015) 75:1467–74.

Google Scholar

6. Wang M, Takahashi A, Liu F, Ye D, Ding Q, Qin C, et al. Large-scale association analysis in Asians identifies new susceptibility loci for prostate cancer. Nat Commun. (2015) 6:8469.

Google Scholar

7. Oh JJ, Lee SJ, Hwang JY, Kim D, Lee SE, Hong SK, et al. Exome-based genome-wide association study and risk assessment using genetic risk score to prostate cancer in the Korean population. Oncotarget. (2017) 8:43934–43. doi: 10.18632/oncotarget.16540

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Bush WS, Moore JH. Chapter 11: genome-wide association studies. PLoS Comput Biol. (2012) 8:e1002822. doi: 10.1371/journal.pcbi.1002822

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Torkamani A, Wineinger NE, Topol EJ. The personal and clinical utility of polygenic risk scores. Nat Rev Genet. (2018) 19:581–90. doi: 10.1038/s41576-018-0018-x

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Zheng SL, Sun J, Wiklund F, Smith S, Stattin P, Li G, et al. Cumulative association of five genetic variants with prostate cancer. N Engl J Med. (2008) 358:910–9.

Google Scholar

11. Eeles RA, Olama AA, Benlloch S, Saunders EJ, Leongamornlert DA, Tymrakiewicz M, et al. Identification of 23 new prostate cancer susceptibility loci using the iCOGS custom genotyping array. Nat Genet. (2013) 45:385–91.

Google Scholar

12. Du Z, Hopp H, Ingles SA, Huff C, Sheng X, Weaver B, et al. A genome-wide association study of prostate cancer in Latinos. Int J Cancer. (2020) 146:1819–26.

Google Scholar

13. Cho YS, Go MJ, Kim YJ, Heo JY, Oh JH, Ban HJ, et al. A large-scale genome-wide association study of Asian populations uncovers genetic factors influencing eight quantitative traits. Nat Genet. (2009) 41:527–34. doi: 10.1038/ng.357

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Grove ML, Yu B, Cochran BJ, Haritunians T, Bis JC, Taylor KD, et al. Best practices and joint calling of the HumanExome BeadChip: the CHARGE consortium. PLoS One. (2013) 8:e68095. doi: 10.1371/journal.pone.0068095

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Huyghe JR, Jackson AU, Fogarty MP, Buchkovich ML, Stancakova A, Stringham HM, et al. Exome array analysis identifies new loci and low-frequency variants influencing insulin processing and secretion. Nat Genet. (2013) 45:197–201. doi: 10.1038/ng.2507

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Hanley JA, McNeil BJ. The meaning and use of the area under a receiver operating characteristic (ROC) curve. Radiology. (1982) 143:29–36. doi: 10.1148/radiology.143.1.7063747

PubMed Abstract | CrossRef Full Text | Google Scholar

17. DeLong ER, DeLong DM, Clarke-Pearson DL. Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics. (1988) 44:837–45. doi: 10.2307/2531595

CrossRef Full Text | Google Scholar

18. Fluss R, Faraggi D, Reiser B. Estimation of the Youden Index and its associated cutoff point. Biom J. (2005) 47:458–72. doi: 10.1002/bimj.200410135

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Chen H, Liu X, Brendler CB, Ankerst DP, Leach RJ, Goodman PJ, et al. Adding genetic risk score to family history identifies twice as many high-risk men for prostate cancer: results from the prostate cancer prevention trial. Prostate. (2016) 76:1120–9. doi: 10.1002/pros.23200

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Schumacher FR, Berndt SI, Siddiq A, Jacobs KB, Wang Z, Lindstrom S, et al. Genome-wide association study identifies new prostate cancer susceptibility loci. Hum Mol Genet. (2011) 20:3867–75.

Google Scholar

21. Yamada H, Penney KL, Takahashi H, Katoh T, Yamano Y, Yamakado M, et al. Replication of prostate cancer risk loci in a Japanese case-control association study. J Natl Cancer Inst. (2009) 101:1330–6. doi: 10.1093/jnci/djp287

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Zheng SL, Hsing AW, Sun J, Chu LW, Yu K, Li G, et al. Association of 17 prostate cancer susceptibility loci with prostate cancer risk in Chinese men. Prostate. (2010) 70:425–32. doi: 10.1002/pros.21076

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Thomas G, Jacobs KB, Yeager M, Kraft P, Wacholder S, Orr N, et al. Multiple loci identified in a genome-wide association study of prostate cancer. Nat Genet. (2008) 40:310–5.

Google Scholar

24. Kim HJ, Bae JS, Lee J, Chang IH, Kim KD, Shin HD, et al. HNF1B polymorphism associated with development of prostate cancer in Korean patients. Urology. (2011) 78:969.e1–6. doi: 10.1016/j.urology.2011.06.045

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Gudmundsson J, Sulem P, Steinthorsdottir V, Bergthorsson JT, Thorleifsson G, Manolescu A, et al. Two variants on chromosome 17 confer prostate cancer risk, and the one in TCF2 protects against type 2 diabetes. Nat Genet. (2007) 39:977–83.

Google Scholar

26. Hoffmann TJ, Van Den Eeden SK, Sakoda LC, Jorgenson E, Habel LA, Graff RE, et al. A large multiethnic genome-wide association study of prostate cancer identifies novel risk variants and substantial ethnic differences. Cancer Discov. (2015) 5:878–91.

Google Scholar

27. Kader AK, Sun J, Reck BH, Newcombe PJ, Kim ST, Hsu FC, et al. Potential impact of adding genetic markers to clinical parameters in predicting prostate biopsy outcomes in men following an initial negative biopsy: findings from the REDUCE trial. Eur Urol. (2012) 62:953–61. doi: 10.1016/j.eururo.2012.05.006

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Seibert TM, Fan CC, Wang Y, Zuber V, Karunamuni R, Parsons JK, et al. Polygenic hazard score to guide screening for aggressive prostate cancer: development and validation in large scale cohorts. BMJ. (2018) 360: j5757.

Google Scholar

29. Khera AV, Emdin CA, Drake I, Natarajan P, Bick AG, Cook NR, et al. Genetic risk, adherence to a healthy lifestyle, and coronary disease. N Engl J Med. (2016) 375:2349–58.

Google Scholar

30. Natarajan P, Young R, Stitziel NO, Padmanabhan S, Baber U, Mehran R, et al. Polygenic risk score identifies subgroup with higher burden of atherosclerosis and greater relative benefit from statin therapy in the primary prevention setting. Circulation. (2017) 135:2091–101. doi: 10.1161/circulationaha.116.024436

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Rutten-Jacobs LC, Larsson SC, Malik R, Rannikmäe K, Sudlow CL, Dichgans M, et al. Genetic risk, incident stroke, and the benefits of adhering to a healthy lifestyle: cohort study of 306?473 UK Biobank participants. BMJ. (2018) 363:k4168.

Google Scholar

32. Fang J, Gong C, Wan Y, Xu Y, Tao F, Sun Y. Polygenic risk, adherence to a healthy lifestyle, and childhood obesity. Pediatr Obes. (2019) 14:e12489.

Google Scholar

Keywords: prediction, Genomics, polygenic risk score, Koreans, prostate cancer

Citation: Oh JJ, Kim E, Woo E, Song SH, Kim JK, Lee H, Lee S, Hong SK and Byun S-S (2020) Evaluation of Polygenic Risk Scores for Prediction of Prostate Cancer in Korean Men. Front. Oncol. 10:583625. doi: 10.3389/fonc.2020.583625

Received: 15 July 2020; Accepted: 16 September 2020;
Published: 22 October 2020.

Edited by:

Michael T. Schweizer, University of Washington, United States

Reviewed by:

Benjamin A. Teply, University of Nebraska Medical Center, United States
Laura Graham, University of Washington, United States
Mark C. Markowski, Johns Hopkins Medicine, United States

Copyright © 2020 Oh, Kim, Woo, Song, Kim, Lee, Lee, Hong and Byun. 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(s) 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: Seok-Soo Byun, ssbyun@snubh.org

These authors have contributed equally to this work