Using gene and gene-set association tests to identify lethal prostate cancer genes

Background: Recent advances in the detection and treatment of prostate cancer (PCa) have reduced morbidity and mortality from this common cancer. Despite these improvements, PCa remains the second leading cause of cancer death in men in the United States. Further understanding of the genetic underpinnings of lethal PCa is required to drive risk detection and prevention and ultimately reduce mortality. We therefore set out to identify germline variants associated with cases of lethal prostate cancer (LPCa).

Methods: Using a two-stage study design, we compared whole-exome sequencing data of 550 LPCa patients to 488 healthy male controls. Men were classified as having LPCa based on medical record review. Candidate genes were identified using gene- and gene-set-based rare truncating variant association tests. Case-control burden testing through Firth's penalized logistic regression and case-gnomAD allelic burden testing through a one-sided mid-p Fisher's exact test were conducted. Each gene's p-values from these tests were combined into an omnibus p-value for candidate gene selection. In the subsequent validation stage, genes were assessed using the UK Biobank and Firth's penalized logistic regression for each ancestry, combined through meta-analysis.

Results: Gene-based rare variant association tests identified 12 genes nominally associated with LPCa. Rare-variant association tests identified a gene set with a significantly higher burden of truncating germline mutations in LPCa patients than controls. Combining gene- and gene-set test results, four nominally significant genes (PPP1R3A, TG, PPFIBP2, and BTN3A3) were selected as candidates. Subsequent validation using the UK Biobank found that PPP1R3A was significantly associated with LPCa risk (odds ratio 2.34, CI 1.20-4.59). Specifically, pGln662ArgfsTer7 was identified as the predominant variant in PPP1R3A among LPCa patients in our dataset.

Conclusions: Both individual gene and gene-set analyses identified candidates associated with LPCa. The novel association of PPP1R3A and LPCa risk merits further investigation.