The effect of sample size on polygenic hazard models for prostate cancer

Comparing Pruning and Thresholding with Continuous Shrinkage Polygenic Score Methods in a Large Sample of Ancestrally Diverse Adolescents from the ABCD Study®

Using individuals' genetic data researchers can generate Polygenic Scores (PS) that are able to predict risk for diseases, variability in different behaviors as well as anthropomorphic measures. This is achieved by leveraging models learned from previously published large Genome-Wide Association Studies (GWASs) associating locations in the genome with a phenotype of interest. Previous GWASs have predominantly been performed in European ancestry individuals. This is of concern as PS generated in samples with a different ancestry to the original training GWAS have been shown to have lower performance and limited portability, and many efforts are now underway to collect genetic databases on individuals of diverse ancestries. In this study, we compare multiple methods of generating PS, including pruning and thresholding and Bayesian continuous shrinkage models, to determine which of them is best able to overcome these limitations. To do this we use the ABCD Study, a longitudinal cohort with deep phenotyping on individuals of diverse ancestry. We generate PS for anthropometric and psychiatric phenotypes using previously published GWAS summary statistics and examine their performance in three subsamples of ABCD: African ancestry individuals (n = 811), European ancestry Individuals (n = 6703), and admixed ancestry individuals (n = 3664). We find that the single ancestry continuous shrinkage method, PRScs (CS), and the multi ancestry meta method, PRScsx Meta (CSx Meta), show the best performance across ancestries and phenotypes.

Polygenic risk of any, metastatic, and fatal prostate cancer in the Million Veteran Program

BACKGROUND

Genetic scores may provide an objective measure of prostate cancer risk and thus inform screening decisions. We evaluated whether a polygenic hazard score based on 290 genetic variants (PHS290) is associated with prostate cancer risk in a diverse population, including Black men, who have higher average risk of prostate cancer death but are often treated as a homogeneously high-risk group.

METHODS

This was a retrospective analysis of the Million Veteran Program, a national, population-based cohort study of US military veterans conducted 2011-2021. Cox proportional hazards analyses tested for association of genetic and other risk factors (including self-reported race and ethnicity and family history) with age at death from prostate cancer, age at diagnosis of metastatic (nodal or distant) prostate cancer, and age at diagnosis of any prostate cancer.

RESULTS

A total of 590 750 male participants were included. Median age at last follow-up was 69 years. PHS290 was associated with fatal prostate cancer in the full cohort and for each racial and ethnic group (P < .001). Comparing men in the highest 20% of PHS290 with those in the lowest 20% (based on percentiles from an independent training cohort), the hazard ratio for fatal prostate cancer was 4.42 (95% confidence interval = 3.91 to 5.02). When accounting for guideline-recommended risk factors (family history, race, and ethnicity), PHS290 remained a strong independent predictor of any, metastatic, and fatal prostate cancer.

CONCLUSIONS

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

Prostate cancer risk stratification improvement across multiple ancestries with new polygenic hazard score

BACKGROUND

Prostate cancer risk stratification using single-nucleotide polymorphisms (SNPs) demonstrates considerable promise in men of European, Asian, and African genetic ancestries, but there is still need for increased accuracy. We evaluated whether including additional SNPs in a prostate cancer polygenic hazard score (PHS) would improve associations with clinically significant prostate cancer in multi-ancestry datasets.

METHODS

In total, 299 SNPs previously associated with prostate cancer were evaluated for inclusion in a new PHS, using a LASSO-regularized Cox proportional hazards model in a training dataset of 72,181 men from the PRACTICAL Consortium. The PHS model was evaluated in four testing datasets: African ancestry, Asian ancestry, and two of European Ancestry-the Cohort of Swedish Men (COSM) and the ProtecT study. Hazard ratios (HRs) were estimated to compare men with high versus low PHS for association with clinically significant, with any, and with fatal prostate cancer. The impact of genetic risk stratification on the positive predictive value (PPV) of PSA testing for clinically significant prostate cancer was also measured.

RESULTS

The final model (PHS290) had 290 SNPs with non-zero coefficients. Comparing, for example, the highest and lowest quintiles of PHS290, the hazard ratios (HRs) for clinically significant prostate cancer were 13.73 [95% CI: 12.43-15.16] in ProtecT, 7.07 [6.58-7.60] in African ancestry, 10.31 [9.58-11.11] in Asian ancestry, and 11.18 [10.34-12.09] in COSM. Similar results were seen for association with any and fatal prostate cancer. Without PHS stratification, the PPV of PSA testing for clinically significant prostate cancer in ProtecT was 0.12 (0.11-0.14). For the top 20% and top 5% of PHS290, the PPV of PSA testing was 0.19 (0.15-0.22) and 0.26 (0.19-0.33), respectively.

CONCLUSIONS

We demonstrate better genetic risk stratification for clinically significant prostate cancer than prior versions of PHS in multi-ancestry datasets. This is promising for implementing precision-medicine approaches to prostate cancer screening decisions in diverse populations.

Performance of African-ancestry-specific polygenic hazard score varies according to local ancestry in 8q24

BACKGROUND

We previously developed an African-ancestry-specific polygenic hazard score (PHS46+African) that substantially improved prostate cancer risk stratification in men with African ancestry. The model consists of 46 SNPs identified in Europeans and 3 SNPs from 8q24 shown to improve model performance in Africans. Herein, we used principal component (PC) analysis to uncover subpopulations of men with African ancestry for whom the utility of PHS46+African may differ.

MATERIALS AND METHODS

Genotypic data were obtained from the PRACTICAL consortium for 6253 men with African genetic ancestry. Genetic variation in a window spanning 3 African-specific 8q24 SNPs was estimated using 93 PCs. A Cox proportional hazards framework was used to identify the pair of PCs most strongly associated with the performance of PHS46+African. A calibration factor (CF) was formulated using Cox coefficients to quantify the extent to which the performance of PHS46+African varies with PC.

RESULTS

CF of PHS46+African was strongly associated with the first and twentieth PCs. Predicted CF ranged from 0.41 to 2.94, suggesting that PHS46+African may be up to 7 times more beneficial to some African men than others. The explained relative risk for PHS46+African varied from 3.6% to 9.9% for individuals with low and high CF values, respectively. By cross-referencing our data set with 1000 Genomes, we identified significant associations between continental and calibration groupings.

CONCLUSION

We identified PCs within 8q24 that were strongly associated with the performance of PHS46+African. Further research to improve the clinical utility of polygenic risk scores (or models) is needed to improve health outcomes for men of African ancestry.

Genetic Stratification of Age-Dependent Parkinson's Disease Risk by Polygenic Hazard Score

BACKGROUND

Parkinson's disease (PD) is a highly age-related disorder, where common genetic risk variants affect both disease risk and age at onset. A statistical approach that integrates these effects across all common variants may be clinically useful for individual risk stratification. A polygenic hazard score methodology, leveraging a time-to-event framework, has recently been successfully applied in other age-related disorders.

OBJECTIVES

We aimed to develop and validate a polygenic hazard score model in sporadic PD.

METHODS

Using a Cox regression framework, we modeled the polygenic hazard score in a training data set of 11,693 PD patients and 9841 controls. The score was then validated in an independent test data set of 5112 PD patients and 5372 controls and a small single-study sample of 360 patients and 160 controls.

RESULTS

A polygenic hazard score predicts the onset of PD with a hazard ratio of 3.78 (95% confidence interval 3.49-4.10) when comparing the highest to the lowest risk decile. Combined with epidemiological data on incidence rate, we apply the score to estimate genetically stratified instantaneous PD risk across age groups.

CONCLUSIONS

We demonstrate the feasibility of a polygenic hazard approach in PD, integrating the genetic effects on disease risk and age at onset in a single model. In combination with other predictive biomarkers, the approach may hold promise for risk stratification in future clinical trials of disease-modifying therapies, which aim at postponing the onset of PD. © 2021 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Common genetic and clinical risk factors: association with fatal prostate cancer in the Cohort of Swedish Men

BACKGROUND

Clinical variables-age, family history, genetics-are used for prostate cancer risk stratification. Recently, polygenic hazard scores (PHS46, PHS166) were validated as associated with age at prostate cancer diagnosis. While polygenic scores are associated with all prostate cancer (not specific for fatal cancers), PHS46 was also associated with age at prostate cancer death. We evaluated if adding PHS to clinical variables improves associations with prostate cancer death.

METHODS

Genotype/phenotype data were obtained from a nested case-control Cohort of Swedish Men (n = 3279; 2163 with prostate cancer, 278 prostate cancer deaths). PHS and clinical variables (family history, alcohol intake, smoking, heart disease, hypertension, diabetes, body mass index) were tested via univariable Cox proportional hazards models for association with age at prostate cancer death. Multivariable Cox models with/without PHS were compared with log-likelihood tests.

RESULTS

Median age at last follow-up/prostate cancer death was 78.0 (IQR: 72.3-84.1) and 81.4 (75.4-86.3) years, respectively. On univariable analysis, PHS46 (HR 3.41 [95% CI 2.78-4.17]), family history (HR 1.72 [1.46-2.03]), alcohol (HR 1.74 [1.40-2.15]), diabetes (HR 0.53 [0.37-0.75]) were each associated with prostate cancer death. On multivariable analysis, PHS46 (HR 2.45 [1.99-2.97]), family history (HR 1.73 [1.48-2.03]), alcohol (HR 1.45 [1.19-1.76]), diabetes (HR 0.62 [0.42-0.90]) all remained associated with fatal disease. Including PHS46 or PHS166 improved multivariable models for fatal prostate cancer (p < 10^-15).

CONCLUSIONS

PHS had the most robust association with fatal prostate cancer in a multivariable model with common risk factors, including family history. Adding PHS to clinical variables may improve prostate cancer risk stratification strategies.

Polygenic Risk Scores in Prostate Cancer Risk Assessment and Screening

More than 40% of the risk of developing prostate cancer (PCa) is from genetic factors. Genome-wide association studies have led to the discovery of more than 140 variants associated with PCa risk. Polygenic risk scores (PRS) generated using these variants show promise in identifying individuals at much higher (and lower) lifetime risk than the average man. PCa PRS also improve the predictive value of prostate-specific antigen screening, may inform the age for starting PCa screening, and are informative for development of more aggressive tumors. Despite the promise, few clinical trials have evaluated the benefit of PCa PRS for clinical care.

Additional SNPs improve risk stratification of a polygenic hazard score for prostate cancer

BACKGROUND

Polygenic hazard scores (PHS) can identify individuals with increased risk of prostate cancer. We estimated the benefit of additional SNPs on performance of a previously validated PHS (PHS46).

MATERIALS AND METHOD

180 SNPs, shown to be previously associated with prostate cancer, were used to develop a PHS model in men with European ancestry. A machine-learning approach, LASSO-regularized Cox regression, was used to select SNPs and to estimate their coefficients in the training set (75,596 men). Performance of the resulting model was evaluated in the testing/validation set (6,411 men) with two metrics: (1) hazard ratios (HRs) and (2) positive predictive value (PPV) of prostate-specific antigen (PSA) testing. HRs were estimated between individuals with PHS in the top 5% to those in the middle 40% (HR95/50), top 20% to bottom 20% (HR80/20), and bottom 20% to middle 40% (HR20/50). PPV was calculated for the top 20% (PPV80) and top 5% (PPV95) of PHS as the fraction of individuals with elevated PSA that were diagnosed with clinically significant prostate cancer on biopsy.

RESULTS

166 SNPs had non-zero coefficients in the Cox model (PHS166). All HR metrics showed significant improvements for PHS166 compared to PHS46: HR95/50 increased from 3.72 to 5.09, HR80/20 increased from 6.12 to 9.45, and HR20/50 decreased from 0.41 to 0.34. By contrast, no significant differences were observed in PPV of PSA testing for clinically significant prostate cancer.

CONCLUSIONS

Incorporating 120 additional SNPs (PHS166 vs PHS46) significantly improved HRs for prostate cancer, while PPV of PSA testing remained the same.

Polygenic hazard score is associated with prostate cancer in multi-ethnic populations

Genetic models for cancer have been evaluated using almost exclusively European data, which could exacerbate health disparities. A polygenic hazard score (PHS1) is associated with age at prostate cancer diagnosis and improves screening accuracy in Europeans. Here, we evaluate performance of PHS2 (PHS1, adapted for OncoArray) in a multi-ethnic dataset of 80,491 men (49,916 cases, 30,575 controls). PHS2 is associated with age at diagnosis of any and aggressive (Gleason score ≥ 7, stage T3-T4, PSA ≥ 10 ng/mL, or nodal/distant metastasis) cancer and prostate-cancer-specific death. Associations with cancer are significant within European (n = 71,856), Asian (n = 2,382), and African (n = 6,253) genetic ancestries (p < 10-180). Comparing the 80th/20th PHS2 percentiles, hazard ratios for prostate cancer, aggressive cancer, and prostate-cancer-specific death are 5.32, 5.88, and 5.68, respectively. Within European, Asian, and African ancestries, hazard ratios for prostate cancer are: 5.54, 4.49, and 2.54, respectively. PHS2 risk-stratifies men for any, aggressive, and fatal prostate cancer in a multi-ethnic dataset.