Evaluating a polygenic hazard score to predict risk of developing metastatic or fatal prostate cancer in the multi-ancestry Million Veteran Program cohort

Meghana S. Pagadala, Julie Lynch, Roshan Karunamuni, Patrick R. Alba, Kyung Min Lee, Fatai Y. Agiri, Tori Anglin, Hannah Carter, J. Michael Gaziano, Guneet Kaur Jasuja, Rishi Deka, Brent S. Rose, Matthew S. Panizzon, Richard L. Hauger, Tyler M. Seibert
doi: https://doi.org/10.1101/2021.09.24.21264093

Abstract

Background Genetic scores may provide an objective measure of a man’s risk of dying from prostate cancer and thus inform screening decisions, especially in men of African ancestry, who have a higher average risk of prostate cancer death but are often treated as a homogeneous group.

Objective Determine whether a polygenic hazard score based on 290 genetic variants (PHS290) is associated with risk of fatal or metastatic prostate cancer in a diverse population.

Design, Setting, and Participants Retrospective analysis of Million Veteran Program (MVP), a national, population-based cohort study of United States military veterans conducted 2011-2021.

Exposure(s) Genotype data to calculate the genetic score, PHS290. Family history of prostate cancer and ancestry group (European, African, Hispanic, or Asian) were also studied.

Outcome Measurements and Statistical Analysis Primary outcome: age at death from prostate cancer. Key secondary outcome: age at diagnosis of prostate cancer metastases. Statistical analysis: Cox proportional hazards.

Results 582,515 MVP participants. Median age at last follow-up: 69 years. PHS290 was associated with fatal prostate cancer in the full cohort and for each ancestry group (p<10-16). Comparing men in the highest 20% of PHS290 to those in the lowest 20%, the hazard ratio for death from prostate cancer was 4.41 [95% CI: 3.9-5.02]. Corresponding hazard ratios for European, African, Hispanic, and Asian subsets were 4.26 [3.66-4.9], 2.4 [1.77-3.23], 4.72 [2.68-8.87], and 10.46 [2.01-101.0]. When accounting for family history and ancestry group, PHS290 remained a strong independent predictor of fatal prostate cancer. PHS290 was also associated with metastasis.

Conclusions PHS290 stratified US veterans of diverse ancestry for lifetime risk of metastatic or fatal prostate cancer. Predicting genetic risk of lethal prostate cancer with PHS290 might inform individualized decisions about prostate cancer screening.

Patient Summary In a large, diverse population (including over 100,000 African Americans), we evaluated whether a genetic score (PHS290) was associated with a man’s risk of fatal or metastatic prostate cancer. We conclude that men with a high PHS290 genetic score, independent of ancestry or family history, are more likely to die of prostate cancer. This information might be useful to guide which men might benefit from screening for prostate cancer.

Introduction

Prostate cancer is the most diagnosed and second deadliest cancer in men1. Despite the enormous mortality from this disease, early detection of prostate cancer remains controversial. Screening all men via prostate-specific antigen (PSA) testing, regardless of underlying risk, has been shown to reduce prostate cancer deaths by 27% but also results in frequent overdiagnosis of indolent prostate cancer that may never have become symptomatic2–4. These overdiagnoses often lead to unnecessary treatment, with attendant side effects and societal costs. A better strategy is to target PSA screening to those men at higher risk of developing metastatic or fatal prostate cancer.

As one of the most heritable cancers5, genetic risk stratification is a promising approach for identifying individuals at higher risk of developing metastatic or fatal prostate cancer1,3,6.Measures of genetic risk have proven highly effective for predicting lifetime risk of being diagnosed with prostate cancer, outperforming family history or other clinical risk factors7–10. Rather than only predicting lifetime risk, however, an ideal genetic test would focus on clinically significant prostate cancer and estimate age-specific risk. Prostate cancer is highly age dependent, with very low incidence before 50 years of age and increasing exponentially as men get older11,12. Absolute incidence of aggressive prostate cancer also increases with age11,12. Meanwhile, some men with high genetic risk develop aggressive prostate cancer at a younger age and are at particular risk of dying from this disease. Age-specific genetic risk could inform individualized decisions about PSA testing, in the context of a given man’s overall health and competing causes of mortality.

A major limitation of early studies of polygenic risk was an exclusive focus on men of European ancestry13,14. Such systematic bias may exacerbate existing health disparities in prostate cancer incidence and health outcomes15,16. This is particularly worrisome for men of African ancestry, who have a higher overall incidence of metastatic and fatal prostate cancer than men of European or Asian ancestry17,18.

Our group has developed a novel risk prediction tool called a polygenic hazard score (PHS) that identifies men who are likely to develop clinically significant prostate cancers at younger ages. This score, which can be calculated from a single saliva sample at any point in a man’s life, was strongly associated with age at diagnosis of clinically significant prostate cancer in large datasets10,19,20. The score also improved the accuracy of conventional screening with PSA7,10,12. We subsequently expanded the model to optimize performance in men of all ancestries, particularly men with African ancestry19–21. Here, we seek to validate the ability of the PHS to identify men at risk of metastatic or fatal prostate cancer within the Million Veteran Program (MVP), one of the largest and most racially and ethnically diverse populations studied to date22.

Methods

Participants

We retrospectively obtained data from the MVP, composed of individuals between ages 19 to over 100 years who were recruited from 63 Veterans Affairs Medical Centers across the United States (US). Recruitment for the MVP started in 2011, and all veterans were eligible for participation. Consent to participate and permission to re-contact was provided after counseling by research staff and mailing of informational materials. Study participation included consenting to access the participant’s electronic health records for research purposes. The MVP received ethical and study protocol approval from the VA Central Institutional Review Board in accordance with the principles outlined in the Declaration of Helsinki.

Only men were included in this prostate cancer study, comprising 582,515 individuals of European (73.3%), African (17.2%), Hispanic (8.2%), and Asian ancestry (1.3%) (Table 1). There were no inclusion or exclusion criteria for age. Median age at last follow-up was 69 years (interquartile range 59-74 years). Men not meeting the endpoint for each analysis were censored at age at last follow-up. Clinical information used for analyses was retrieved as described below in the Clinical Data Extraction section.

View this table:
Table 1. Participant characteristics, n=582,515

Genotype Data

All study participants provided blood samples for DNA extraction and genotyping. Researchers are provided data that is de-identified except for dates. Blood samples were collected by phlebotomists and banked at the VA Central Biorepository in Boston, MA, where DNA was extracted and shipped to two external centers for genotyping. DNA extracted from buffy coat was genotyped using a custom Affymetrix Axiom biobank array.

The MVP 1.0 genotyping array contains a total of 723,305 single nucleotide polymorphisms (SNPs), enriched for low frequency variants in African and Hispanic populations, and variants associated with diseases common to the VA population22.

Harmonized Ancestry and Race/Ethnicity (HARE)

The MVP has previously developed an approach called Harmonized Ancestry and Race/Ethnicity (HARE) to identify ancestry groups using a machine learning algorithm23. HARE utilizes a support vector machine to estimate probabilities of an individual belonging to one of four ancestry groups using both self-identified race/ethnicity and genetic ancestry23. Information on race and ethnicity was obtained based on self-report through centralized VA data collection methods using standardized survey forms or using information from the VA Corporate Data Warehouse or Observational Medical Outcomes Partnership data. All but 9,989 (1.52%) MVP participants were assigned to a HARE ancestry group. The support vector machine was trained with the top 30 principal components of population stratification analysis and self-identified ancestry. Regularization constant C and inverse variance of kernel were optimized through 2-dimensional grid searching and 5-fold cross-fold validation. Individuals were categorized as predominantly European, Hispanic, African, or Asian based on output probabilities.

Clinical Data Extraction

Each participant’s electronic health record is integrated into the MVP biorepository. These records include International Classification of Diseases (ICD) diagnosis codes (ICD-9-CM and ICD-10-CM), procedure codes (ICD, Current Procedural Terminology, and Healthcare Common Procedure Coding (HCPCS)), laboratory values, medications, and clinical notes documenting VA care (inpatient and outpatient) and non-VA care paid for by the VA.

Prostate cancer diagnosis, age at diagnosis, and date of last follow-up were retrieved from the VA Corporate Data Warehouse based on ICD codes and VA Central Cancer Registry data. Age at diagnosis of metastasis was determined via a validated natural language processing tool and a search of individual participant’s medical records in the Veterans Affairs system, as described previously24. This tool was developed using data from over 1 million VA patients with prostate cancer; compared to manual chart review, the natural language processing tool had 92% sensitivity and 98% specificity for diagnosis of metastatic prostate cancer. Cause and date of death was collected from National Death Index. Participants with ICD10 code “C61” as underlying cause of death were considered to have died from prostate cancer. Age of death was determined from difference between year of death and year of birth.

Polygenic Hazard Score (PHS290)

The most recent version of the PHS, called PHS290, was calculated as the vector product of participants’ genotype dosage (Xi) for 290 SNPs and the corresponding parameter estimates (βi) from Cox proportional hazards regression: Embedded Image

The development of this score has been described elsewhere21. Briefly, 299 previously identified SNPs associated with prostate cancer risk (in single-ancestry or all-ancestry analyses) were simultaneously evaluated using a machine-learning least absolute shrinkage and selection operator (LASSO) approach to generate an optimal combined model for association with age at prostate cancer diagnosis.

We calculated PHS290 for each MVP participant. Distributions were visualized using histograms for each ancestry group. Differences in mean PHS290 between ancestry groups were assessed via ANOVA. In all statistical analyses, significance was set at a two-tailed alpha of 0.01. As in prior studies, p-values less than 10-16 were truncated at this value, as comparison of miniscule values is not likely to be meaningful7,10,12,20.

Cox Proportional Hazards Analysis

We evaluated association of PHS290 with two important clinical endpoints extracted from clinical data: age at death from prostate cancer (i.e., lifetime prostate-cancer-specific mortality) and age at diagnosis of metastases from prostate cancer (i.e., lifetime distant-metastasis-free survival). Secondarily, we also tested for association with age at diagnosis of any prostate cancer. To visualize the association in the full dataset, we generated cause-specific cumulative incidence curves for each endpoint and each of several PHS290 risk groups. Cox proportional hazards models were used to assess these associations in the full dataset and in each ancestral group (European, African, Hispanic, Asian). Individuals not meeting the endpoint of interest were censored at age at last follow-up. Sample-weight corrections were applied to all Cox models to correct for potential bias due to a higher number of prostate cancer cases compared to a general population and permit direct comparisons to other studies7,9,19,25. Although underlying population incidence may vary (e.g., across ancestry groups), a constant correction factor was used for all analyses (in this case, based on previously reported population data from Sweden, for consistency with prior work), as effect size estimates have been shown to be robust to wide variation in population incidence19.

Effect sizes were estimated using hazard ratios (HRs) between risk strata, as described previously7,9,10,12,19,20,26. Percentiles of genetic risk were calculated using percentile thresholds defined in a prior study of men of European ancestry less than 70 years old and with no diagnosis of prostate cancer21. HRs for each ancestry group were calculated to make the following comparisons: HR80/20, men in the highest 20% vs. lowest 20%; HR95/50, men in the highest 5% of genetic risk vs. those with average risk (30–70th percentile); and HR20/50, men in the lowest 20% vs. those with average risk. These risk groups were chosen to mirror prior work and permit direct comparison of HRs7,19,20; the same risk groups were chosen as the strata for incidence curves for each endpoint. Reference thresholds for hazards ratio estimation were: 9.004659 (20th quantile), 9.123500 (30th quantile), 9.519703 (70th quantile), 9.639068 (80th quantile), 9.946332 (95th quantile)21.

Ancestry, Family History, and PHS290

To assess the added value of PHS290 beyond commonly used clinical risk factors, we tested a multivariable Cox proportional hazards model with ancestry group, family history, and PHS2907,9,19. This combined model was limited to the 374,455 participants who provided family history information in baseline survey data. Family history was recorded as either the presence or absence of (one or more) first-degree relatives with prostate cancer. Cox proportional hazards models tested associations with fatal, metastatic, or any prostate cancer. For PHS290, the effect size was illustrated via the hazard ratio for the highest 20% vs. lowest 20% of genetic risk. Hazard ratios for ancestry groups were estimated using European as the reference.

A univariable Cox proportional hazards model was applied to test for association of ancestry group with prostate cancer risk. Similarly, a univariable model tested for association of family history alone. The @anova function from the R ‘survival’ package (version 3.2-13; Therneau 2021) was used to compare the nested Cox models (multivariable vs. univariable), based on the log partial likelihood of the model fits. Significance was set at a two-tailed alpha of 0.01 for the test of whether the multivariable model performed better than either univariable model alone.

Results

PHS290 Score

The distribution of PHS290 in the European ancestry group was similar to that reported previously for men of European ancestry (mean=9.37, SD=0.37)21. Mean PHS290 did vary by ancestry group with statistically significant differences between all groups (ANOVA p < 10-16; all pair-wise t-tests also p < 10-16). The distribution for the Hispanic ancestry group overlapped closely with that of the European group (mean=9.35, SD=0.37), while PHS290 tended to be lower among men of Asian ancestry (mean=9.17, SD=0.36) and higher among men of African ancestry (mean=9.56, SD=0.34, t-test) (Figure 1, Supplementary Figure 1).

Figure 1.
Figure 1. PHS290 Score Density Plots in Million Veteran Program.

PHS290 score density plot in European, African, Hispanic, and Asian ancestry groups.

Association of PHS290 with Fatal Prostate Cancer

PHS290 was associated with age at death from prostate cancer in the full dataset and in each of the four ancestry groups (Table 2). Comparing 80th and 20th percentiles of genetic risk in the full dataset, men with higher PHS290 had an HR80/20 of 4.41 [95% CI: 3.9-5.02]. Cause-specific cumulative incidence curves for various PHS290 percentile groups demonstrated risk stratification (Figure 2). Hazard ratios quantified significant risk stratification using PHS290 in the full MVP dataset and in each ancestry group (though confidence intervals were large in the Asian ancestry group). For European, African, Hispanic, and Asian men, HRs80/20 were 4.26 [95% CI: 3.66-4.9], 2.4 [1.77-3.23], 4.72 [2.68-8.87] and 10.46 [2.01-101.0], respectively.

View this table:
Table 2: Association of PHS290 with fatal prostate cancer.

Cox Proportional-Hazards results from association with age at death from prostate cancer. P-values reported are from univariable models using PHS290 as the sole predictor variable. Hazard ratios (HRs) are shown comparing men in various percentiles of genetic risk. HR80/20: highest 20% (≥80th percentile of PHS290, using previously published thresholds for men <70 years old and no diagnosis of cancer) vs. average risk (30-70th percentile). HR20/50: lowest 20% (≤20th percentile) vs. average risk. HR80/50: highest 20% vs. average risk. HR95/50: highest 5% (≥95th percentile) vs. average risk. Numbers in brackets are 95% confidence intervals.

Figure 2.
Figure 2. Million Veteran Program (MVP) Cause-specific Cumulative Incidence.

Cause-specific cumulative incidence within MVP, stratified by PHS290, for (A) fatal prostate cancer, (B) metastatic prostate cancer, and (C) prostate cancer. PHS290 percentile groups shown for each endpoint: 0-20th, 30-60th, 80-100th, and 95-100th. The y-axis was adjusted for secondary analysis (C) to show the higher incidence values for all prostate cancer.

Association of PHS290 with Metastatic Prostate Cancer

PHS290 was associated with age at diagnosis of metastases from prostate cancer in the full dataset and in each of the four ancestry groups (Supplementary Table 1). Comparing 80th and 20th percentiles of genetic risk, men with higher PHS290 had an HR80/20 of 4.94 [95% CI: 4.58-5.28]. For European, African, Hispanic, and Asian men, HRs80/20 were 4.64 [4.28-5.03], 3.02 [2.6-3.57], 3.95 [2.87-5.37] and 6.85 [2.84-17.12], respectively.

Association of PHS290 with Prostate Cancer

PHS290 was associated with age at prostate cancer diagnosis in all 4 ancestry groups (Supplementary Table 2). Comparing 80th and 20th percentiles of genetic risk, men with higher PHS290 had an HR80/20 of 6.29 [95% CI: 6.14-6.46]. For European, African, Hispanic, and Asian men, HRs80/20 were 6.19 [6.01-6.38], 3.83 [3.61-4.08], 4.75 [4.22-5.32] and 5.52 [3.98-7.56], respectively.

Ancestry, Family History, and PHS290

Ancestry group, alone, was associated with differential risk of fatal prostate cancer (Supplementary Table 3). These associations were largely driven by an increased risk from African ancestry. Similar patterns were seen for age at diagnosis of prostate cancer and for age at diagnosis of prostate cancer metastasis (Supplementary Table

3). Compared to the European group, men in the African ancestry group had a HR of 2.65 [2.37-2.96] for dying of prostate cancer.

Family history, alone, was also associated with fatal prostate cancer, as well as with diagnosis of prostate cancer and with age at diagnosis of prostate cancer metastasis (Supplementary Table 4). Compared to men with no family history of prostate cancer, men with one or more first-degree relatives who had prostate cancer had a HR of 1.84 [1.54-2.17] for dying of prostate cancer.

PHS290 added significant value beyond ancestry group or family history in a multivariable model that included all three variables and tested for association with age at prostate cancer death (Table 3). The multivariable model improved performance over the common risk factors alone (ANOVA p < 10-16). Similarly, the combination proved optimal when evaluating age at diagnosis of prostate cancer diagnosis or age at diagnosis of prostate cancer metastasis (ANOVA p<10-16). Independent of ancestry and family history, a high PHS290 (top 20%) approximately quadrupled a man’s risk of death from prostate cancer, compared to a low PHS290 (bottom 20%) (Table 3).

View this table:
Table 3: Multivariable models combining Ancestry, Family History, and PHS290 for three prostate cancer clinical endpoints.

Cox proportional hazards results for association with age at death from prostate cancer, at diagnosis of metastatic prostate cancer, and age at diagnosis with prostate cancer. P-values reported are from multivariable models using HARE ancestry group, family history, and PHS290. For PHS290, effect size was illustrated via the hazard ratio (HR80/20) for the highest 20% vs. lowest 20% of genetic risk. Hazard ratios for ancestry groups were estimated using European as the reference. Hazard ratios for family history were for one or more first-degree relatives diagnosed with prostate cancer. This analysis was limited to the 374,455 participants who provided family history information in baseline survey data. Numbers in brackets are 95% confidence intervals. Significant predictors in the multivariable model are indicated by * (p<0.01) and *** (p<10-16).

To ensure that genetic ancestry was accounted for, we repeated the multivariable analysis with the top 10 principal components from population stratification analysis (Supplementary Table 5). PHS290 was still significantly associated with all three prostate cancer clinical endpoints.

Discussion

PHS290 was associated with lifetime prostate-cancer-specific mortality in this large and diverse dataset. Even when accounting for family history and ancestry group, PHS290 remained a strong independent predictor of dying from prostate cancer. The genomic score was also associated with age at diagnosis of metastasis from prostate cancer and at age at diagnosis of any prostate cancer, consistent with previous reports that common genetic markers for overall prostate cancer risk overlap with those for aggressive prostate cancer risk27,28. Metastatic prostate cancer has poor prognostic outcomes and is major driver of pain, disability and aggressive medical therapy29. To our knowledge, this study is the first to show the association of a genomic score with lifetime risk of metastatic prostate cancer. This study also represents the largest and most ancestry-diverse independent validation of polygenic association with lifetime risk of fatal prostate cancer.

Men of African ancestry in the US are substantially more likely to develop metastatic disease and to die from prostate cancer18. The causes of this disparity are likely a combination of genetic, environmental, and social factors, including systemic racism30–34. National guidelines recommend consideration of prostate screening in men of African ancestry at a younger age and that screening occur at more frequent intervals35. The results of the present study demonstrate that men of African ancestry have highly variable levels of lifetime risk and should not be treated as a homogeneous group. This finding is consistent with prior results and with the known admixture and genetic diversity of the African American population9,26. Moreover, PHS290 can identify those more likely to develop lethal prostate cancer and may facilitate personalized screening recommendations.

Intriguingly, typical PHS290 scores differed between ancestry groups, with the mean PHS290 slightly higher in the African ancestry group and slightly lower in Asian ancestry group than in the European or Hispanic groups. These shifts in PHS290 distribution are consistent with reported differences in prostate cancer incidence across racial groups36–40. Higher overall PHS290 scores in African ancestry group may point to true differences in prostate cancer risk but could also be inflated by minor allele frequency (MAF) differences between ancestry groups. Incorporating approaches for local ancestry and admixture can also boost genetic model performance and should be explored further to improve the predictive accuracy of polygenic scores41.

Along with race/ancestry, family history is another important clinical consideration in prostate cancer screening decisions35,42–45. Prior studies have found polygenic scores to be the most important risk factor for prostate cancer, with family history sometimes offering modest improvement in multivariable models, possibly by capturing yet unknown genetic factors and/or shared familial environmental factors7,9,19,46. Among the subset of MVP participants who provided family history information, family history of prostate cancer was independently associated with prostate cancer risk in a multivariable model that included PHS290 and ancestry group. The relationship of environmental exposures, family history, and prostate cancer risk are worth further investigation9, particularly in groups like veterans who may have been exposed to rare carcinogens47.

The present study builds on prior work that reported the performance of polygenic scores in non-Europeans19,21,26,27,48 and is consistent with those prior studies in showing a strong association of polygenic scores with prostate cancer risk, including death from prostate cancer9,19,46. Polygenic hazard scores designed to incorporate the strong age-dependence of prostate cancer have also been shown to increase the accuracy of conventional prostate cancer screening7,10,12. Population-level analyses of benefit, harm, and cost-effectiveness support incorporation of genomic risk into screening3,6. The present study adds to the literature an independent validation in a dataset of over half a million men with diverse ancestry. Current clinical guidelines try to achieve targeted, or risk-stratified, screening by recommending each man discuss his individual risk factors, emphasizing racial/ethnic background35,43–45. It is particularly important, therefore, that this study was able to combine ancestry and genetic risk to estimate the relative impact of each and to demonstrate that a polygenic score adds considerable information beyond ancestry alone, for a man’s individual risk of metastasis or death from prostate cancer.

While PHS290 performed well in the present study to stratify men by genetic prostate cancer risk, the effect sizes estimated here are lower than those reported in previous studies. For example, HR80/20 for fatal prostate cancer using PHS290 was 7.73 [95% CI 6.45, 9.27] in a population-based Swedish cohort21, compared to 4.41 [3.9-5.02] in all participants and 4.26 [3.66-4.9] in the European ancestry group in the present study. A similar pattern of smaller effect size was seen when comparing the strength of association with age at diagnosis of prostate cancer within MVP ancestry groups to the hazard ratios previously reported for European, Asian, and African genetic ancestry groups21. Most likely, the discrepancy arises in differences in the populations studied; for example, the MVP dataset comes exclusively from a population of US veterans, with many receiving care in a single US-based healthcare system, whereas the prior study used data from multiple countries and widely varying recruitment strategies. Patterns of screening, detection, and treatment of prostate cancer in the present dataset could be different from clinical trial and case-control datasets used in previous work. Some of the difference in performance could also be explained by the fact that the testing datasets in the prior report for PHS290 had been included in the discovery of a majority of the 290 SNPs in the model; on the other hand, the testing datasets represented a very small proportion of the discovery datasets, and the model weights were estimated in an independent training dataset.

The present work shows that adding PHS290 to guideline-recommended risk factors improves risk stratification for meaningful clinical endpoints of death or metastases from prostate cancer. Men at highest risk of metastatic or fatal prostate cancer are potentially those most likely to benefit from screening. Prior studies have further suggested genetic scores could also add value even after results from screening or diagnostic tests are already available, but this needs further investigation49,50. For example, one early detection strategy with strong evidence is early baseline PSA (e.g., at age 45-49)51,52. This strategy has not yet been widely adopted in the U.S.53, but future studies should evaluate whether PHS290 adds value in men where early baseline PSA is known and whether PHS290, if known prior to PSA testing, should inform the decision of whether to obtain an early baseline PSA test. Another compelling avenue for future studies is whether high genetic risk can be mitigated by lifestyle or other preventive intervention54,55.

Limitations of this study include heterogeneity of phenotyping and smaller sample sizes for Asian ancestry group. Heterogeneity of prostate cancer screening and diagnostic pathways by clinicians across VA and other hospitals in the US could potentially introduce noise, although this heterogeneity likely leads to underestimation of associations with prostate cancer. Large confidence intervals in the Asian ancestry group may be due to relatively smaller sample sizes, but they also suggest increased heterogeneity compared to prior datasets (including Asian ancestry) of similar sample size19,21,56. Finally, we acknowledge that while HARE ancestry groups may be a reasonable attempt to harmonize genetic ancestry, race, and ethnicity, these groups cannot account for—much less, disentangle—the complex web of biological and social factors associated with these categories. Further work will attempt to incorporate agnostic genetic ancestry groups and address impacts of admixture and local/regional genetic ancestry on risk stratification with PHS19.

We show that PHS290 stratified US men for lifetime risk of metastatic or fatal prostate cancer. Critically, this genetic risk stratification was successful within each of four ancestry groups in this diverse dataset. PHS290 was higher, on average, among men with African ancestry, who were also at higher risk from prostate cancer.

The combination of ancestry, family history, and PHS290 performed better than any variable, alone, in identifying men at highest risk of prostate cancer metastasis and death. Predicting genetic risk of lethal prostate cancer with PHS290 might inform individualized decisions about screening and early cancer detection.

Data Availability

Requests regarding data access may be directed to MVPLOI{at}va.gov

Competing Interests

TMS reports honoraria from Varian Medical Systems and WebMD; he has an equity interest in CorTechs Labs, Inc. and serves on its Scientific Advisory Board; he has received in-kind research support from GE Healthcare via a research agreement with the University of California San Diego. These companies might potentially benefit from the research results. The terms of this arrangement have been reviewed and approved by the University of California San Diego in accordance with its conflict-of-interest policies.

Data Availability Statement

Requests regarding data access may be directed to MVPLOI{at}va.gov.

Take-Home Message

PHS290 stratified US veterans of diverse ancestry for lifetime risk of fatal or metastatic prostate cancer, even when accounting for family history and ancestry. Predicting genetic risk of lethal prostate cancer might inform individualized decisions about prostate cancer screening.

Acknowledgements

This research used data from the Million Veteran Program, Office of Research and Development, Veterans Health Administration. This research was supported by the Million Veteran Program MVP022 award # I01 CX001727 (PI: Richard L. Hauger MD). This publication does not represent the views of the Department of Veterans Affairs or the United States Government. Dr. Hauger was additionally funded by the VISN-22 VA Center of Excellence for Stress and Mental Health (CESAMH) and National Institute of Aging RO1 grant AG050595 (The VETSA Longitudinal Twin Study of Cognition and Aging VETSA 4). This research was supported by VA MVP022. Meghana S. Pagadala was supported by the National Institutes of Health (#1F30CA247168, #T32CA067754). Tyler Seibert and Roshan Karunamuni were supported by the National Institutes of Health (NIH/NIBIB #K08EB026503), the Prostate Cancer Foundation, and the University of California (#C21CR2060).

Footnotes

  • Correction of a typo in number of participants. Addition of supplementary analysis re: top 10 principal components. Expansion of Discussion section.

References

  1. 1.↵
    Torkamani A, Wineinger NE, Topol EJ. The personal and clinical utility of polygenic risk scores. Nat Rev Genet. 2018;19(9):581–590.
    OpenUrlCrossRefPubMed
  2. 2.↵
    Schröder FH, Hugosson J, Roobol MJ, et al. Screening and prostate cancer mortality: results of the European Randomised Study of Screening for Prostate Cancer (ERSPC) at 13 years of follow-up. Lancet. 2014;384(9959):2027–2035.
    OpenUrlCrossRefPubMedWeb of Science
  3. 3.↵
    Callender T, Emberton M, Morris S, et al. Polygenic risk-tailored screening for prostate cancer: A benefit-harm and cost-effectiveness modelling study. PLoS Med. 2019;16(12):e1002998.
    OpenUrlPubMed
  4. 4.↵
    Loeb S, Bjurlin MA, Nicholson J, et al. Overdiagnosis and overtreatment of prostate cancer. Eur Urol. 2014;65(6):1046–1055.
    OpenUrlCrossRefPubMed
  5. 5.↵
    Mucci LA, Hjelmborg JB, Harris JR, et al. Familial risk and heritability of cancer among twins in Nordic countries. JAMA. 2016;315(1):68–76.
    OpenUrlCrossRefPubMed
  6. 6.↵
    Callender T, Emberton M, Morris S, Pharoah PDP, Pashayan N. Benefit, harm, and cost-effectiveness associated with magnetic resonance imaging before biopsy in age-based and risk-stratified screening for prostate cancer. JAMA Netw Open. 2021;4(3):e2037657.
    OpenUrl
  7. 7.↵
    Seibert TM, Fan CC, Wang Y, et al. Polygenic hazard score to guide screening for aggressive prostate cancer: development and validation in large scale cohorts. BMJ. 2018;2018:j5757.
    OpenUrl
  8. 8.
    Schumacher FR, Olama AAA, Berndt SI, et al. Author Correction: Association analyses of more than 140,000 men identify 63 new prostate cancer susceptibility loci. Nat Genet. 2019;51(2):363.
    OpenUrl
  9. 9.↵
    Huynh-Le MP, Karunamuni R, Fan CC, et al. Common genetic and clinical risk factors: association with fatal prostate cancer in the Cohort of Swedish Men. Prostate Cancer Prostatic Dis. Published online March 15, 2021. doi:10.1038/s41391-021-00341-4
    OpenUrlCrossRef
  10. 10.↵
    Karunamuni RA, Huynh-Le MP, Fan CC, et al. Additional SNPs improve risk stratification of a polygenic hazard score for prostate cancer. Prostate Cancer Prostatic Dis. 2021;24(2):532–541.
    OpenUrl
  11. 11.↵
    Huynh-Le MP, Myklebust TÅ, Feng CH, et al. Age dependence of modern clinical risk groups for localized prostate cancer-A population-based study. Cancer. 2020;126(8):1691–1699.
    OpenUrl
  12. 12.↵
    Huynh-Le MP, Fan CC, Karunamuni R, et al. A genetic risk score to personalize Prostate Cancer screening, applied to population data. Cancer Epidemiol Biomarkers Prev. 2020;29(9):1731–1738.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    Duncan L, Shen H, Gelaye B, et al. Analysis of polygenic risk score usage and performance in diverse human populations. Nat Commun. 2019;10(1):3328.
    OpenUrlPubMed
  14. 14.↵
    Petrovski S, Goldstein DB. Unequal representation of genetic variation across ancestry groups creates healthcare inequality in the application of precision medicine. Genome Biol. 2016;17(1). doi:10.1186/s13059-016-1016-y
    OpenUrlCrossRef
  15. 15.↵
    Martin AR, Kanai M, Kamatani Y, Okada Y, Neale BM, Daly MJ. Clinical use of current polygenic risk scores may exacerbate health disparities. Nat Genet. 2019;51(4):584–591.
    OpenUrlCrossRefPubMed
  16. 16.↵
    Popejoy AB, Fullerton SM. Genomics is failing on diversity. Nature. 2016;538(7624):161–164.
    OpenUrlCrossRefPubMed
  17. 17.↵
    DeSantis CE, Miller KD, Goding Sauer A, Jemal A, Siegel RL. Cancer statistics for African Americans, 2019. CA Cancer J Clin. 2019;69(3):211–233.
    OpenUrlPubMed
  18. 18.↵
    Tsodikov A, Gulati R, de Carvalho TM, et al. Is prostate cancer different in black men? Answers from 3 natural history models. Cancer. 2017;123(12):2312–2319.
    OpenUrlCrossRef
  19. 19.↵
    Huynh-Le MP, Fan CC, Karunamuni R, et al. Polygenic hazard score is associated with prostate cancer in multi-ethnic populations. Nat Commun. 2021;12(1):1236.
    OpenUrl
  20. 20.↵
    Karunamuni RA, Huynh-Le MP, Fan CC, et al. African-specific improvement of a polygenic hazard score for age at diagnosis of prostate cancer. Int J Cancer. 2021;148(1):99–105.
    OpenUrl
  21. 21.↵
    Huynh-Le MP, Karunamuni R, Fan CC, et al. Prostate cancer risk stratification improved across multiple ancestries with new polygenic hazard score. Published online August 18, 2021. doi:10.1101/2021.08.14.21261931
    OpenUrlAbstract/FREE Full Text
  22. 22.↵
    Gaziano JM, Concato J, Brophy M, et al. Million Veteran Program: A mega-biobank to study genetic influences on health and disease. J Clin Epidemiol. 2016;2016:214–223.
    OpenUrl
  23. 23.↵
    Fang H, Hui Q, Lynch J, et al. Harmonizing genetic ancestry and self-identified race/ethnicity in genome-wide association studies. Am J Hum Genet. 2019;105(4):763–772.
    OpenUrlCrossRef
  24. 24.↵
    Alba PR, Gao A, Lee KM, et al. Ascertainment of Veterans with Metastatic Prostate Cancer in Electronic Health Records: Demonstrating the Case for Natural Language Processing. JCO Clin Cancer Inform. 2021;2021:1005–1014.
    OpenUrl
  25. 25.↵
    Therneau TM, Li H. Computing the Cox model for case cohort designs. Lifetime Data Anal. 1999;5(2):99–112.
    OpenUrlCrossRefPubMedWeb of Science
  26. 26.↵
    Karunamuni RA, Huynh-Le MP, Fan CC, et al. Performance of African-ancestry-specific polygenic hazard score varies according to local ancestry in 8q24. Prostate Cancer Prostatic Dis. Published online June 14, 2021. doi:10.1038/s41391-021-00403-7
    OpenUrlCrossRef
  27. 27.↵
    Conti DV, Darst BF, Moss LC, et al. Trans-ancestry genome-wide association meta-analysis of prostate cancer identifies new susceptibility loci and informs genetic risk prediction. Nat Genet. 2021;53(1):65–75.
    OpenUrlCrossRef
  28. 28.↵
    Schumacher FR, Al Olama AA, Berndt SI, et al. Association analyses of more than 140,000 men identify 63 new prostate cancer susceptibility loci. Nat Genet. 2018;50(7):928–936.
    OpenUrlCrossRefPubMed
  29. 29.↵
    Sartor O, de Bono JS. Metastatic prostate cancer. N Engl J Med. 2018;378(7):645–657.
    OpenUrlCrossRefPubMed
  30. 30.↵
    Barocas DA, Grubb R 3rd., Black A, et al. Association between race and follow-up diagnostic care after a positive prostate cancer screening test in the prostate, lung, colorectal, and ovarian cancer screening trial. Cancer. 2013;119(12):2223–2229.
    OpenUrlCrossRefPubMed
  31. 31.
    Han Y, Rand KA, Hazelett DJ, et al. Prostate cancer susceptibility in men of African ancestry at 8q24. J Natl Cancer Inst. 2016;108(7):djv431.
    OpenUrlCrossRefPubMed
  32. 32.
    Mahal BA, Aizer AA, Ziehr DR, et al. Trends in disparate treatment of African American men with localized prostate cancer across National Comprehensive Cancer Network risk groups. Urology. 2014;84(2):386–392.
    OpenUrlPubMed
  33. 33.
    Yamoah K, Johnson MH, Choeurng V, et al. Novel biomarker signature that may predict aggressive disease in African American men with prostate cancer. J Clin Oncol. 2015;33(25):2789–2796.
    OpenUrlAbstract/FREE Full Text
  34. 34.↵
    Zhang H, Messing EM, Travis LB, et al. Age and racial differences among PSA-detected (AJCC stage T1cN0M0) prostate cancer in the U.s.: A population-based study of 70,345 men. Front Oncol. 2013;2013:312.
    OpenUrl
  35. 35.↵
    Carroll PR, Parsons JK, Box G, Carlsson S, Catalona WJ, Drake BF. NCCN Guidelines Prostate Cancer Early Detection Version 1. 2021. Published online 2021.
  36. 36.↵
    Siegel DA, O’Neil ME, Richards TB, Dowling NF, Weir HK. Prostate cancer incidence and survival, by stage and race/ethnicity - United States, 2001-2017. MMWR Morb Mortal Wkly Rep. 2020;69(41):1473–1480.
    OpenUrl
  37. 37.
    Cook MB, Rosenberg PS, McCarty FA, et al. Racial disparities in prostate cancer incidence rates by census division in the United States, 1999-2008. Prostate. 2015;75(7):758–763.
    OpenUrlPubMed
  38. 38.
    Atlas of cancer mortality in the United States, 1950-94. Int J Epidemiol. 2000;29(3):602-a-602.
    OpenUrlCrossRef
  39. 39.
    Jemal A, Ward E, Wu X, Martin HJ, McLaughlin CC, Thun MJ. Geographic patterns of prostate cancer mortality and variations in access to medical care in the United States. Cancer Epidemiol Biomarkers Prev. 2005;14(3):590–595.
    OpenUrlAbstract/FREE Full Text
  40. 40.↵
    Curry SJ, Krist AH, Owens DK. Annual report to the nation on the status of cancer, part II: Recent changes in prostate cancer trends and disease characteristics. Cancer. 2019;125(2):317–318.
    OpenUrl
  41. 41.↵
    Marnetto D, Pärna K, Läll K, et al. Ancestry deconvolution and partial polygenic score can improve susceptibility predictions in recently admixed individuals. Nat Commun. 2020;11(1):1628.
    OpenUrlCrossRef
  42. 42.↵
    Powell IJ. The precise role of ethnicity and family history on aggressive prostate cancer: a review analysis. Arch Esp Urol. 2011;64(8):711–719.
    OpenUrlPubMed
  43. 43.↵
    Wolf AMD, Wender RC, Etzioni RB, et al. American Cancer Society guideline for the early detection of prostate cancer: update 2010. CA Cancer J Clin. 2010;60(2):70–98.
    OpenUrlCrossRefPubMedWeb of Science
  44. 44.
    Horwich A, Hugosson J, de Reijke T, et al. Prostate cancer: ESMO Consensus Conference Guidelines 2012. Ann Oncol. 2013;24(5):1141–1162.
    OpenUrlCrossRefPubMedWeb of Science
  45. 45.↵
    Qaseem A, Barry MJ, Denberg TD, Owens DK, Shekelle P, Clinical Guidelines Committee of the American College of Physicians. Screening for prostate cancer: a guidance statement from the Clinical Guidelines Committee of the American College of Physicians. Ann Intern Med. 2013;158(10):761–769.
    OpenUrlCrossRefPubMedWeb of Science
  46. 46.↵
    Shi Z, Platz EA, Wei J, et al. Performance of three inherited risk measures for predicting prostate cancer incidence and mortality: A population-based prospective analysis. Eur Urol. 2021;79(3):419–426.
    OpenUrl
  47. 47.↵
    Schecter A, Needham L, Pavuk M, et al. Agent Orange exposure, Vietnam war veterans, and the risk of prostate cancer. Cancer. 2009;115(14):3369–3371.
    OpenUrlCrossRefPubMed
  48. 48.
    Plym A, Penney KL, Kalia S, et al. Evaluation of a multiethnic polygenic risk score model for prostate cancer. J Natl Cancer Inst. Published online April 1, 2021. doi:10.1093/jnci/djab058
    OpenUrlCrossRef
  49. 49.↵
    Eklund M, Jäderling F, Discacciati A, et al. MRI-targeted or standard biopsy in prostate cancer screening. N Engl J Med. 2021;385(10):908–920.
    OpenUrl
  50. 50.↵
    Jiang Y, Meyers TJ, Emeka AA, et al. Genetic factors associated with prostate cancer conversion from active surveillance to treatment. HGG Adv. 2022;3(1):100070.
    OpenUrl
  51. 51.↵
    Vickers AJ, Ulmert D, Sjoberg DD, et al. Strategy for detection of prostate cancer based on relation between prostate specific antigen at age 40-55 and long term risk of metastasis: case-control study. BMJ. 2013;346(apr15 5):f2023.
    OpenUrlAbstract/FREE Full Text
  52. 52.↵
    Lilja H, Ulmert D, Björk T, et al. Long-term prediction of prostate cancer up to 25 years before diagnosis of prostate cancer using prostate kallikreins measured at age 44 to 50 years. J Clin Oncol. 2007;25(4):431–436.
    OpenUrlAbstract/FREE Full Text
  53. 53.↵
    US Preventive Services Task Force, Grossman DC, Curry SJ, et al. Screening for prostate cancer. JAMA. 2018;319(18):1901.
    OpenUrlCrossRefPubMed
  54. 54.↵
    Plym A, Zhang Y, Stopsack K, et al. Abstract 822: Can the genetic risk of prostate cancer be attenuated by a healthy lifestyle. In: Epidemiology. American Association for Cancer Research; 2021. doi:10.1158/1538-7445.am2021-822
    OpenUrlCrossRef
  55. 55.↵
    Ahmed M, Goh C, Saunders E, et al. Germline genetic variation in prostate susceptibility does not predict outcomes in the chemoprevention trials PCPT and SELECT. Prostate Cancer Prostatic Dis. 2020;23(2):333–342.
    OpenUrl
  56. 56.↵
    Karunamuni RA, Huynh-Le MP, Fan CC, et al. The effect of sample size on polygenic hazard models for prostate cancer. Eur J Hum Genet. 2020;28(10):1467–1475.
    OpenUrl
Back to top
Posted February 16, 2022.
Download PDF

Supplementary Material

Data/Code
Email
Evaluating a polygenic hazard score to predict risk of developing metastatic or fatal prostate cancer in the multi-ancestry Million Veteran Program cohort
Meghana S. Pagadala, Julie Lynch, Roshan Karunamuni, Patrick R. Alba, Kyung Min Lee, Fatai Y. Agiri, Tori Anglin, Hannah Carter, J. Michael Gaziano, Guneet Kaur Jasuja, Rishi Deka, Brent S. Rose, Matthew S. Panizzon, Richard L. Hauger, Tyler M. Seibert
medRxiv 2021.09.24.21264093; doi: https://doi.org/10.1101/2021.09.24.21264093
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area