Characterizing prostate cancer risk through multi-ancestry genome-wide discovery of 187 novel risk variants.

The transferability and clinical value of genetic risk scores (GRSs) across populations remain limited due to an imbalance in genetic studies across ancestrally diverse populations. Here we conducted a multi-ancestry genome-wide association study of 156,319 prostate cancer cases and 788,443 controls of European, African, Asian and Hispanic men, reflecting a 57% increase in the number of non-European cases over previous prostate cancer genome-wide association studies. We identified 187 novel risk variants for prostate cancer, increasing the total number of risk variants to 451. An externally replicated multi-ancestry GRS was associated with risk that ranged from 1.8 (per standard deviation) in African ancestry men to 2.2 in European ancestry men. The GRS was associated with a greater risk of aggressive versus non-aggressive disease in men of African ancestry (P = 0.03). Our study presents novel prostate cancer susceptibility loci and a GRS with effective risk stratification across ancestry groups.

Trans-ancestry genome-wide association meta-analysis of prostate cancer identifies new susceptibility loci and informs genetic risk prediction.

Prostate cancer is a highly heritable disease with large disparities in incidence rates across ancestry populations. We conducted a multiancestry meta-analysis of prostate cancer genome-wide association studies (107,247 cases and 127,006 controls) and identified 86 new genetic risk variants independently associated with prostate cancer risk, bringing the total to 269 known risk variants. The top genetic risk score (GRS) decile was associated with odds ratios that ranged from 5.06 (95% confidence interval (CI), 4.84-5.29) for men of European ancestry to 3.74 (95% CI, 3.36-4.17) for men of African ancestry. Men of African ancestry were estimated to have a mean GRS that was 2.18-times higher (95% CI, 2.14-2.22), and men of East Asian ancestry 0.73-times lower (95% CI, 0.71-0.76), than men of European ancestry. These findings support the role of germline variation contributing to population differences in prostate cancer risk, with the GRS offering an approach for personalized risk prediction.

Publisher Correction: Shared heritability and functional enrichment across six solid cancers.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Author Correction: Germline variation at 8q24 and prostate cancer risk in men of European ancestry.

The original version of this Article contained an error in the spelling of the author Manuela Gago-Dominguez, which was incorrectly given as Manuela G. Dominguez. This has now been corrected in both the PDF and HTML versions of the Article.

Author Correction: Association analyses of more than 140,000 men identify 63 new prostate cancer susceptibility loci.

In the version of this article initially published, the name of author Manuela Gago-Dominguez was misspelled as Manuela Gago Dominguez. The error has been corrected in the HTML and PDF version of the article.

Shared heritability and functional enrichment across six solid cancers.

Quantifying the genetic correlation between cancers can provide important insights into the mechanisms driving cancer etiology. Using genome-wide association study summary statistics across six cancer types based on a total of 296,215 cases and 301,319 controls of European ancestry, here we estimate the pair-wise genetic correlations between breast, colorectal, head/neck, lung, ovary and prostate cancer, and between cancers and 38 other diseases. We observed statistically significant genetic correlations between lung and head/neck cancer (rg = 0.57, p = 4.6 × 10-8), breast and ovarian cancer (rg = 0.24, p = 7 × 10-5), breast and lung cancer (rg = 0.18, p =1.5 × 10-6) and breast and colorectal cancer (rg = 0.15, p = 1.1 × 10-4). We also found that multiple cancers are genetically correlated with non-cancer traits including smoking, psychiatric diseases and metabolic characteristics. Functional enrichment analysis revealed a significant excess contribution of conserved and regulatory regions to cancer heritability. Our comprehensive analysis of cross-cancer heritability suggests that solid tumors arising across tissues share in part a common germline genetic basis.

Circulating Metabolic Biomarkers of Screen-Detected Prostate Cancer in the ProtecT Study.

BACKGROUND: Whether associations between circulating metabolites and prostate cancer are causal is unknown. We report on the largest study of metabolites and prostate cancer (2,291 cases and 2,661 controls) and appraise causality for a subset of the prostate cancer-metabolite associations using two-sample Mendelian randomization (MR). METHODS: The case-control portion of the study was conducted in nine UK centers with men ages 50-69 years who underwent prostate-specific antigen screening for prostate cancer within the Prostate Testing for Cancer and Treatment (ProtecT) trial. Two data sources were used to appraise causality: a genome-wide association study (GWAS) of metabolites in 24,925 participants and a GWAS of prostate cancer in 44,825 cases and 27,904 controls within the Association Group to Investigate Cancer Associated Alterations in the Genome (PRACTICAL) consortium. RESULTS: Thirty-five metabolites were strongly associated with prostate cancer (P < 0.0014, multiple-testing threshold). These fell into four classes: (i) lipids and lipoprotein subclass characteristics (total cholesterol and ratios, cholesterol esters and ratios, free cholesterol and ratios, phospholipids and ratios, and triglyceride ratios); (ii) fatty acids and ratios; (iii) amino acids; (iv) and fluid balance. Fourteen top metabolites were proxied by genetic variables, but MR indicated these were not causal. CONCLUSIONS: We identified 35 circulating metabolites associated with prostate cancer presence, but found no evidence of causality for those 14 testable with MR. Thus, the 14 MR-tested metabolites are unlikely to be mechanistically important in prostate cancer risk. IMPACT: The metabolome provides a promising set of biomarkers that may aid prostate cancer classification.

Germline variation at 8q24 and prostate cancer risk in men of European ancestry.

Chromosome 8q24 is a susceptibility locus for multiple cancers, including prostate cancer. Here we combine genetic data across the 8q24 susceptibility region from 71,535 prostate cancer cases and 52,935 controls of European ancestry to define the overall contribution of germline variation at 8q24 to prostate cancer risk. We identify 12 independent risk signals for prostate cancer (p < 4.28 × 10-15), including three risk variants that have yet to be reported. From a polygenic risk score (PRS) model, derived to assess the cumulative effect of risk variants at 8q24, men in the top 1% of the PRS have a 4-fold (95%CI = 3.62-4.40) greater risk compared to the population average. These 12 variants account for ~25% of what can be currently explained of the familial risk of prostate cancer by known genetic risk factors. These findings highlight the overwhelming contribution of germline variation at 8q24 on prostate cancer risk which has implications for population risk stratification.

AA9int: SNP interaction pattern search using non-hierarchical additive model set.

Motivation: The use of single nucleotide polymorphism (SNP) interactions to predict complex diseases is getting more attention during the past decade, but related statistical methods are still immature. We previously proposed the SNP Interaction Pattern Identifier (SIPI) approach to evaluate 45 SNP interaction patterns/patterns. SIPI is statistically powerful but suffers from a large computation burden. For large-scale studies, it is necessary to use a powerful and computation-efficient method. The objective of this study is to develop an evidence-based mini-version of SIPI as the screening tool or solitary use and to evaluate the impact of inheritance mode and model structure on detecting SNP-SNP interactions. Results: We tested two candidate approaches: the 'Five-Full' and 'AA9int' method. The Five-Full approach is composed of the five full interaction models considering three inheritance modes (additive, dominant and recessive). The AA9int approach is composed of nine interaction models by considering non-hierarchical model structure and the additive mode. Our simulation results show that AA9int has similar statistical power compared to SIPI and is superior to the Five-Full approach, and the impact of the non-hierarchical model structure is greater than that of the inheritance mode in detecting SNP-SNP interactions. In summary, it is recommended that AA9int is a powerful tool to be used either alone or as the screening stage of a two-stage approach (AA9int+SIPI) for detecting SNP-SNP interactions in large-scale studies. Availability and implementation: The 'AA9int' and 'parAA9int' functions (standard and parallel computing version) are added in the SIPI R package, which is freely available at https://linhuiyi.github.io/LinHY_Software/. Supplementary information: Supplementary data are available at Bioinformatics online.

Association analyses of more than 140,000 men identify 63 new prostate cancer susceptibility loci.

Genome-wide association studies (GWAS) and fine-mapping efforts to date have identified more than 100 prostate cancer (PrCa)-susceptibility loci. We meta-analyzed genotype data from a custom high-density array of 46,939 PrCa cases and 27,910 controls of European ancestry with previously genotyped data of 32,255 PrCa cases and 33,202 controls of European ancestry. Our analysis identified 62 novel loci associated (P < 5.0 × 10-8) with PrCa and one locus significantly associated with early-onset PrCa (≤55 years). Our findings include missense variants rs1800057 (odds ratio (OR) = 1.16; P = 8.2 × 10-9; G>C, p.Pro1054Arg) in ATM and rs2066827 (OR = 1.06; P = 2.3 × 10-9; T>G, p.Val109Gly) in CDKN1B. The combination of all loci captured 28.4% of the PrCa familial relative risk, and a polygenic risk score conferred an elevated PrCa risk for men in the ninetieth to ninety-ninth percentiles (relative risk = 2.69; 95% confidence interval (CI): 2.55-2.82) and first percentile (relative risk = 5.71; 95% CI: 5.04-6.48) risk stratum compared with the population average. These findings improve risk prediction, enhance fine-mapping, and provide insight into the underlying biology of PrCa1.

Fine-mapping of prostate cancer susceptibility loci in a large meta-analysis identifies candidate causal variants.

Prostate cancer is a polygenic disease with a large heritable component. A number of common, low-penetrance prostate cancer risk loci have been identified through GWAS. Here we apply the Bayesian multivariate variable selection algorithm JAM to fine-map 84 prostate cancer susceptibility loci, using summary data from a large European ancestry meta-analysis. We observe evidence for multiple independent signals at 12 regions and 99 risk signals overall. Only 15 original GWAS tag SNPs remain among the catalogue of candidate variants identified; the remainder are replaced by more likely candidates. Biological annotation of our credible set of variants indicates significant enrichment within promoter and enhancer elements, and transcription factor-binding sites, including AR, ERG and FOXA1. In 40 regions at least one variant is colocalised with an eQTL in prostate cancer tissue. The refined set of candidate variants substantially increase the proportion of familial relative risk explained by these known susceptibility regions, which highlights the importance of fine-mapping studies and has implications for clinical risk profiling.

Gene and pathway level analyses of germline DNA-repair gene variants and prostate cancer susceptibility using the iCOGS-genotyping array.

This corrects the article DOI: 10.1038/bjc.2016.50.

SNP interaction pattern identifier (SIPI): an intensive search for SNP-SNP interaction patterns.

Motivation: Testing SNP-SNP interactions is considered as a key for overcoming bottlenecks of genetic association studies. However, related statistical methods for testing SNP-SNP interactions are underdeveloped. Results: We propose the SNP Interaction Pattern Identifier (SIPI), which tests 45 biologically meaningful interaction patterns for a binary outcome. SIPI takes non-hierarchical models, inheritance modes and mode coding direction into consideration. The simulation results show that SIPI has higher power than MDR (Multifactor Dimensionality Reduction), AA_Full, Geno_Full (full interaction model with additive or genotypic mode) and SNPassoc in detecting interactions. Applying SIPI to the prostate cancer PRACTICAL consortium data with approximately 21 000 patients, the four SNP pairs in EGFR-EGFR , EGFR-MMP16 and EGFR-CSF1 were found to be associated with prostate cancer aggressiveness with the exact or similar pattern in the discovery and validation sets. A similar match for external validation of SNP-SNP interaction studies is suggested. We demonstrated that SIPI not only searches for more meaningful interaction patterns but can also overcome the unstable nature of interaction patterns. Availability and Implementation: The SIPI software is freely available at http://publichealth.lsuhsc.edu/LinSoftware/ . Contact: hlin1@lsuhsc.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

Investigating the possible causal role of coffee consumption with prostate cancer risk and progression using Mendelian randomization analysis.

Coffee consumption has been shown in some studies to be associated with lower risk of prostate cancer. However, it is unclear if this association is causal or due to confounding or reverse causality. We conducted a Mendelian randomisation analysis to investigate the causal effects of coffee consumption on prostate cancer risk and progression. We used two genetic variants robustly associated with caffeine intake (rs4410790 and rs2472297) as proxies for coffee consumption in a sample of 46,687 men of European ancestry from 25 studies in the PRACTICAL consortium. Associations between genetic variants and prostate cancer case status, stage and grade were assessed by logistic regression and with all-cause and prostate cancer-specific mortality using Cox proportional hazards regression. There was no clear evidence that a genetic risk score combining rs4410790 and rs2472297 was associated with prostate cancer risk (OR per additional coffee increasing allele: 1.01, 95% CI: 0.98,1.03) or having high-grade compared to low-grade disease (OR: 1.01, 95% CI: 0.97,1.04). There was some evidence that the genetic risk score was associated with higher odds of having nonlocalised compared to localised stage disease (OR: 1.03, 95% CI: 1.01, 1.06). Amongst men with prostate cancer, there was no clear association between the genetic risk score and all-cause mortality (HR: 1.00, 95% CI: 0.97,1.04) or prostate cancer-specific mortality (HR: 1.03, 95% CI: 0.98,1.08). These results, which should have less bias from confounding than observational estimates, are not consistent with a substantial effect of coffee consumption on reducing prostate cancer incidence or progression.

Alcohol consumption and prostate cancer incidence and progression: A Mendelian randomisation study.

Prostate cancer is the most common cancer in men in developed countries, and is a target for risk reduction strategies. The effects of alcohol consumption on prostate cancer incidence and survival remain unclear, potentially due to methodological limitations of observational studies. In this study, we investigated the associations of genetic variants in alcohol-metabolising genes with prostate cancer incidence and survival. We analysed data from 23,868 men with prostate cancer and 23,051 controls from 25 studies within the international PRACTICAL Consortium. Study-specific associations of 68 single nucleotide polymorphisms (SNPs) in 8 alcohol-metabolising genes (Alcohol Dehydrogenases (ADHs) and Aldehyde Dehydrogenases (ALDHs)) with prostate cancer diagnosis and prostate cancer-specific mortality, by grade, were assessed using logistic and Cox regression models, respectively. The data across the 25 studies were meta-analysed using fixed-effect and random-effects models. We found little evidence that variants in alcohol metabolising genes were associated with prostate cancer diagnosis. Four variants in two genes exceeded the multiple testing threshold for associations with prostate cancer mortality in fixed-effect meta-analyses. SNPs within ALDH1A2 associated with prostate cancer mortality were rs1441817 (fixed effects hazard ratio, HRfixed = 0.78; 95% confidence interval (95%CI):0.66,0.91; p values = 0.002); rs12910509, HRfixed = 0.76; 95%CI:0.64,0.91; p values = 0.003); and rs8041922 (HRfixed = 0.76; 95%CI:0.64,0.91; p values = 0.002). These SNPs were in linkage disequilibrium with each other. In ALDH1B1, rs10973794 (HRfixed = 1.43; 95%CI:1.14,1.79; p values = 0.002) was associated with prostate cancer mortality in men with low-grade prostate cancer. These results suggest that alcohol consumption is unlikely to affect prostate cancer incidence, but it may influence disease progression.

Prostate cancer risk regions at 8q24 and 17q24 are differentially associated with somatic TMPRSS2:ERG fusion status.

Molecular and epidemiological differences have been described between TMPRSS2:ERG fusion-positive and fusion-negative prostate cancer (PrCa). Assuming two molecularly distinct subtypes, we have examined 27 common PrCa risk variants, previously identified in genome-wide association studies, for subtype specific associations in a total of 1221 TMPRSS2:ERG phenotyped PrCa cases. In meta-analyses of a discovery set of 552 cases with TMPRSS2:ERG data and 7650 unaffected men from five centers we have found support for the hypothesis that several common risk variants are associated with one particular subtype rather than with PrCa in general. Risk variants were analyzed in case-case comparisons (296 TMPRSS2:ERG fusion-positive versus 256 fusion-negative cases) and an independent set of 669 cases with TMPRSS2:ERG data was established to replicate the top five candidates. Significant differences (P < 0.00185) between the two subtypes were observed for rs16901979 (8q24) and rs1859962 (17q24), which were enriched in TMPRSS2:ERG fusion-negative (OR = 0.53, P = 0.0007) and TMPRSS2:ERG fusion-positive PrCa (OR = 1.30, P = 0.0016), respectively. Expression quantitative trait locus analysis was performed to investigate mechanistic links between risk variants, fusion status and target gene mRNA levels. For rs1859962 at 17q24, genotype dependent expression was observed for the candidate target gene SOX9 in TMPRSS2:ERG fusion-positive PrCa, which was not evident in TMPRSS2:ERG negative tumors. The present study established evidence for the first two common PrCa risk variants differentially associated with TMPRSS2:ERG fusion status. TMPRSS2:ERG phenotyping of larger studies is required to determine comprehensive sets of variants with subtype-specific roles in PrCa.

Polyunsaturated fatty acids and prostate cancer risk: a Mendelian randomisation analysis from the PRACTICAL consortium.

BACKGROUND: Prostate cancer is a common cancer worldwide with no established modifiable lifestyle factors to guide prevention. The associations between polyunsaturated fatty acids (PUFAs) and prostate cancer risk have been inconsistent. Using Mendelian randomisation, we evaluated associations between PUFAs and prostate cancer risk. METHODS: We used individual-level data from a consortium of 22 721 cases and 23 034 controls of European ancestry. Externally-weighted PUFA-specific polygenic risk scores (wPRSs), with explanatory variation ranging from 0.65 to 33.07%, were constructed and used to evaluate associations with prostate cancer risk per one standard deviation (s.d.) increase in genetically-predicted plasma PUFA levels using multivariable-adjusted unconditional logistic regression. RESULTS: No overall association was observed between the genetically-predicted PUFAs evaluated in this study and prostate cancer risk. However, risk reductions were observed for short-chain PUFAs, linoleic (ORLA=0.95, 95%CI=0.92, 0.98) and α-linolenic acids (ORALA=0.96, 95%CI=0.93, 0.98), among men <62 years; whereas increased risk was found among men ⩾62 years for LA (ORLA=1.04, 95%CI=1.01, 1.07). For long-chain PUFAs (i.e., arachidonic, eicosapentaenoic, and docosapentaenoic acids), increased risks were observed among men <62 years (ORAA=1.05, 95%CI=1.02, 1.08; OREPA=1.04, 95%CI=1.01, 1.06; ORDPA=1.05, 95%CI=1.02, 1.08). CONCLUSION: Results from this study suggest that circulating ω-3 and ω-6 PUFAs may have a different role in the aetiology of early- and late-onset prostate cancer.

Assessing the role of insulin-like growth factors and binding proteins in prostate cancer using Mendelian randomization: Genetic variants as instruments for circulating levels.

Circulating insulin-like growth factors (IGFs) and their binding proteins (IGFBPs) are associated with prostate cancer. Using genetic variants as instruments for IGF peptides, we investigated whether these associations are likely to be causal. We identified from the literature 56 single nucleotide polymorphisms (SNPs) in the IGF axis previously associated with biomarker levels (8 from a genome-wide association study [GWAS] and 48 in reported candidate genes). In ∼700 men without prostate cancer and two replication cohorts (N ∼ 900 and ∼9,000), we examined the properties of these SNPS as instrumental variables (IVs) for IGF-I, IGF-II, IGFBP-2 and IGFBP-3. Those confirmed as strong IVs were tested for association with prostate cancer risk, low (< 7) vs. high (≥ 7) Gleason grade, localised vs. advanced stage, and mortality, in 22,936 controls and 22,992 cases. IV analysis was used in an attempt to estimate the causal effect of circulating IGF peptides on prostate cancer. Published SNPs in the IGFBP1/IGFBP3 gene region, particularly rs11977526, were strong instruments for IGF-II and IGFBP-3, less so for IGF-I. Rs11977526 was associated with high (vs. low) Gleason grade (OR per IGF-II/IGFBP-3 level-raising allele 1.05; 95% CI: 1.00, 1.10). Using rs11977526 as an IV we estimated the causal effect of a one SD increase in IGF-II (∼265 ng/mL) on risk of high vs. low grade disease as 1.14 (95% CI: 1.00, 1.31). Because of the potential for pleiotropy of the genetic instruments, these findings can only causally implicate the IGF pathway in general, not any one specific biomarker.

Pubertal development and prostate cancer risk: Mendelian randomization study in a population-based cohort.

BACKGROUND: Epidemiological studies have observed a positive association between an earlier age at sexual development and prostate cancer, but markers of sexual maturation in boys are imprecise and observational estimates are likely to suffer from a degree of uncontrolled confounding. To obtain causal estimates, we examined the role of pubertal development in prostate cancer using genetic polymorphisms associated with Tanner stage in adolescent boys in a Mendelian randomization (MR) approach. METHODS: We derived a weighted genetic risk score for pubertal development, combining 13 SNPs associated with male Tanner stage. A higher score indicated a later puberty onset. We examined the association of this score with prostate cancer risk, stage and grade in the UK-based ProtecT case-control study (n = 2,927), and used the PRACTICAL consortium (n = 43,737) as a replication sample. RESULTS: In ProtecT, the puberty genetic score was inversely associated with prostate cancer grade (odds ratio (OR) of high- vs. low-grade cancer, per tertile of the score: 0.76; 95 % CI, 0.64-0.89). In an instrumental variable estimation of the causal OR, later physical development in adolescence (equivalent to a difference of one Tanner stage between pubertal boys of the same age) was associated with a 77 % (95 % CI, 43-91 %) reduced odds of high Gleason prostate cancer. In PRACTICAL, the puberty genetic score was associated with prostate cancer stage (OR of advanced vs. localized cancer, per tertile: 0.95; 95 % CI, 0.91-1.00) and prostate cancer-specific mortality (hazard ratio amongst cases, per tertile: 0.94; 95 % CI, 0.90-0.98), but not with disease grade. CONCLUSIONS: Older age at sexual maturation is causally linked to a reduced risk of later prostate cancer, especially aggressive disease.