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

We determined the effect of sample size on performance of polygenic hazard score (PHS) models in prostate cancer. Age and genotypes were obtained for 40,861 men from the PRACTICAL consortium. The dataset included 201,590 SNPs per subject, and was split into training and testing sets. Established-SNP models considered 65 SNPs that had been previously associated with prostate cancer. Discovery-SNP models used stepwise selection to identify new SNPs. The performance of each PHS model was calculated for random sizes of the training set. The performance of a representative Established-SNP model was estimated for random sizes of the testing set. Mean HR98/50 (hazard ratio of top 2% to average in test set) of the Established-SNP model increased from 1.73 [95% CI: 1.69-1.77] to 2.41 [2.40-2.43] when the number of training samples was increased from 1 thousand to 30 thousand. Corresponding HR98/50 of the Discovery-SNP model increased from 1.05 [0.93-1.18] to 2.19 [2.16-2.23]. HR98/50 of a representative Established-SNP model using testing set sample sizes of 0.6 thousand and 6 thousand observations were 1.78 [1.70-1.85] and 1.73 [1.71-1.76], respectively. We estimate that a study population of 20 thousand men is required to develop Discovery-SNP PHS models while 10 thousand men should be sufficient for Established-SNP models.

A Genetic Risk Score to Personalize Prostate Cancer Screening, Applied to Population Data

BACKGROUND

A polygenic hazard score (PHS), the weighted sum of 54 SNP genotypes, was previously validated for association with clinically significant prostate cancer and for improved prostate cancer screening accuracy. Here, we assess the potential impact of PHS-informed screening.

METHODS

United Kingdom population incidence data (Cancer Research United Kingdom) and data from the Cluster Randomized Trial of PSA Testing for Prostate Cancer were combined to estimate age-specific clinically significant prostate cancer incidence (Gleason score ≥7, stage T3-T4, PSA ≥10, or nodal/distant metastases). Using HRs estimated from the ProtecT prostate cancer trial, age-specific incidence rates were calculated for various PHS risk percentiles. Risk-equivalent age, when someone with a given PHS percentile has prostate cancer risk equivalent to an average 50-year-old man (50-year-standard risk), was derived from PHS and incidence data. Positive predictive value (PPV) of PSA testing for clinically significant prostate cancer was calculated using PHS-adjusted age groups.

RESULTS

The expected age at diagnosis of clinically significant prostate cancer differs by 19 years between the 1st and 99th PHS percentiles: men with PHS in the 1st and 99th percentiles reach the 50-year-standard risk level at ages 60 and 41, respectively. PPV of PSA was higher for men with higher PHS-adjusted age.

CONCLUSIONS

PHS provides individualized estimates of risk-equivalent age for clinically significant prostate cancer. Screening initiation could be adjusted by a man's PHS.

IMPACT

Personalized genetic risk assessments could inform prostate cancer screening decisions.

Polygenic hazard score to guide screening for aggressive prostate cancer: development and validation in large scale cohorts

OBJECTIVES

To develop and validate a genetic tool to predict age of onset of aggressive prostate cancer (PCa) and to guide decisions of who to screen and at what age.

DESIGN

Analysis of genotype, PCa status, and age to select single nucleotide polymorphisms (SNPs) associated with diagnosis. These polymorphisms were incorporated into a survival analysis to estimate their effects on age at diagnosis of aggressive PCa (that is, not eligible for surveillance according to National Comprehensive Cancer Network guidelines; any of Gleason score ≥7, stage T3-T4, PSA (prostate specific antigen) concentration ≥10 ng/L, nodal metastasis, distant metastasis). The resulting polygenic hazard score is an assessment of individual genetic risk. The final model was applied to an independent dataset containing genotype and PSA screening data. The hazard score was calculated for these men to test prediction of survival free from PCa.

SETTING

Multiple institutions that were members of international PRACTICAL consortium.

PARTICIPANTS

All consortium participants of European ancestry with known age, PCa status, and quality assured custom (iCOGS) array genotype data. The development dataset comprised 31 747 men; the validation dataset comprised 6411 men.

MAIN OUTCOME MEASURES

Prediction with hazard score of age of onset of aggressive cancer in validation set.

RESULTS

In the independent validation set, the hazard score calculated from 54 single nucleotide polymorphisms was a highly significant predictor of age at diagnosis of aggressive cancer (z=11.2, P<10-16). When men in the validation set with high scores (>98th centile) were compared with those with average scores (30th-70th centile), the hazard ratio for aggressive cancer was 2.9 (95% confidence interval 2.4 to 3.4). Inclusion of family history in a combined model did not improve prediction of onset of aggressive PCa (P=0.59), and polygenic hazard score performance remained high when family history was accounted for. Additionally, the positive predictive value of PSA screening for aggressive PCa was increased with increasing polygenic hazard score.

CONCLUSIONS

Polygenic hazard scores can be used for personalised genetic risk estimates that can predict for age at onset of aggressive PCa.

The 30 kb deletion in the APOBEC3 cluster decreases APOBEC3A and APOBEC3B expression and creates a transcriptionally active hybrid gene but does not associate with breast cancer in the European population

APOBEC3B, in addition to other members of the APOBEC3 gene family, has recently been intensively studied due to its identification as a gene whose activation in cancer is responsible for a specific pattern of massively occurring somatic mutations. It was recently shown that a common large deletion in the APOBEC3 cluster (the APOBEC3B deletion) may increase the risk of breast cancer. However, conflicting evidence regarding this association was also reported. In the first step of our study, using different approaches, including an in-house designed multiplex ligation-dependent probe amplification assay, we analyzed the structure of the deletion and showed that although the breakpoints are located in highly homologous regions, which may generate recurrent occurrence of similar but not identical deletions, there is no sign of deletion heterogeneity. This knowledge allowed us to distinguish transcripts of all affected genes, including the highly homologous canonical APOBEC3A and APOBEC3B, and the hybrid APOBEC3A/APOBEC3B gene. We unambiguously confirmed the presence of the hybrid transcript and showed that the APOBEC3B deletion negatively correlates with APOBEC3A and APOBEC3B expression and positively correlates with APOBEC3A/APOBEC3B expression, whose mRNA level is >10-fold and >1500-fold lower than the level of APOBEC3A and APOBEC3B, respectively. In the next step, we performed a large-scale association study in three different cohorts (2972 cases and 3682 controls) and showed no association of the deletion with breast cancer, familial breast cancer or ovarian cancer. Further, we conducted a meta-analysis that confirmed the lack of the association of the deletion with breast cancer in non-Asian populations.

Genome-Wide Meta-Analyses of Breast, Ovarian, and Prostate Cancer Association Studies Identify Multiple New Susceptibility Loci Shared by at Least Two Cancer Types

Breast, ovarian, and prostate cancers are hormone-related and may have a shared genetic basis, but this has not been investigated systematically by genome-wide association (GWA) studies. Meta-analyses combining the largest GWA meta-analysis data sets for these cancers totaling 112,349 cases and 116,421 controls of European ancestry, all together and in pairs, identified at P < 10(-8) seven new cross-cancer loci: three associated with susceptibility to all three cancers (rs17041869/2q13/BCL2L11; rs7937840/11q12/INCENP; rs1469713/19p13/GATAD2A), two breast and ovarian cancer risk loci (rs200182588/9q31/SMC2; rs8037137/15q26/RCCD1), and two breast and prostate cancer risk loci (rs5013329/1p34/NSUN4; rs9375701/6q23/L3MBTL3). Index variants in five additional regions previously associated with only one cancer also showed clear association with a second cancer type. Cell-type-specific expression quantitative trait locus and enhancer-gene interaction annotations suggested target genes with potential cross-cancer roles at the new loci. Pathway analysis revealed significant enrichment of death receptor signaling genes near loci with P < 10(-5) in the three-cancer meta-analysis.

SIGNIFICANCE

We demonstrate that combining large-scale GWA meta-analysis findings across cancer types can identify completely new risk loci common to breast, ovarian, and prostate cancers. We show that the identification of such cross-cancer risk loci has the potential to shed new light on the shared biology underlying these hormone-related cancers. Cancer Discov; 6(9); 1052-67. ©2016 AACR.This article is highlighted in the In This Issue feature, p. 932.

Atlas of prostate cancer heritability in European and African-American men pinpoints tissue-specific regulation

Although genome-wide association studies have identified over 100 risk loci that explain ∼33% of familial risk for prostate cancer (PrCa), their functional effects on risk remain largely unknown. Here we use genotype data from 59,089 men of European and African American ancestries combined with cell-type-specific epigenetic data to build a genomic atlas of single-nucleotide polymorphism (SNP) heritability in PrCa. We find significant differences in heritability between variants in prostate-relevant epigenetic marks defined in normal versus tumour tissue as well as between tissue and cell lines. The majority of SNP heritability lies in regions marked by H3k27 acetylation in prostate adenoc7arcinoma cell line (LNCaP) or by DNaseI hypersensitive sites in cancer cell lines. We find a high degree of similarity between European and African American ancestries suggesting a similar genetic architecture from common variation underlying PrCa risk. Our findings showcase the power of integrating functional annotation with genetic data to understand the genetic basis of PrCa.

Multiple novel prostate cancer susceptibility signals identified by fine-mapping of known risk loci among Europeans

Genome-wide association studies (GWAS) have identified numerous common prostate cancer (PrCa) susceptibility loci. We have fine-mapped 64 GWAS regions known at the conclusion of the iCOGS study using large-scale genotyping and imputation in 25 723 PrCa cases and 26 274 controls of European ancestry. We detected evidence for multiple independent signals at 16 regions, 12 of which contained additional newly identified significant associations. A single signal comprising a spectrum of correlated variation was observed at 39 regions; 35 of which are now described by a novel more significantly associated lead SNP, while the originally reported variant remained as the lead SNP only in 4 regions. We also confirmed two association signals in Europeans that had been previously reported only in East-Asian GWAS. Based on statistical evidence and linkage disequilibrium (LD) structure, we have curated and narrowed down the list of the most likely candidate causal variants for each region. Functional annotation using data from ENCODE filtered for PrCa cell lines and eQTL analysis demonstrated significant enrichment for overlap with bio-features within this set. By incorporating the novel risk variants identified here alongside the refined data for existing association signals, we estimate that these loci now explain ∼38.9% of the familial relative risk of PrCa, an 8.9% improvement over the previously reported GWAS tag SNPs. This suggests that a significant fraction of the heritability of PrCa may have been hidden during the discovery phase of GWAS, in particular due to the presence of multiple independent signals within the same region.