Cybulski, Cezary

Abstract 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/. Supplementary information Supplementary data are available at Bioinformatics online.

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.

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.

Background The rarity of mutations in PALB2, CHEK2 and ATM make it difficult to estimate precisely associated cancer risks. Population-based family studies have provided evidence that at least some of these mutations are associated with breast cancer risk as high as those associated with rare BRCA2 mutations. We aimed to estimate the relative risks associated with specific rare variants in PALB2, CHEK2 and ATM via a multicentre case-control study. Methods We genotyped 10 rare mutations using the custom iCOGS array: PALB2 c.1592delT, c.2816T>G and c.3113G>A, CHEK2 c.349A>G, c.538C>T, c.715G>A, c.1036C>T, c.1312G>T, and c.1343T>G and ATM c.7271T>G. We assessed associations with breast cancer risk (42 671 cases and 42 164 controls), as well as prostate (22 301 cases and 22 320 controls) and ovarian (14 542 cases and 23 491 controls) cancer risk, for each variant. Results For European women, strong evidence of association with breast cancer risk was observed for PALB2 c.1592delT OR 3.44 (95% CI 1.39 to 8.52, p=7.1×10−5), PALB2 c.3113G>A OR 4.21 (95% CI 1.84 to 9.60, p=6.9×10−8) and ATM c.7271T>G OR 11.0 (95% CI 1.42 to 85.7, p=0.0012). We also found evidence of association with breast cancer risk for three variants in CHEK2, c.349A>G OR 2.26 (95% CI 1.29 to 3.95), c.1036C>T OR 5.06 (95% CI 1.09 to 23.5) and c.538C>T OR 1.33 (95% CI 1.05 to 1.67) (p≤0.017). Evidence for prostate cancer risk was observed for CHEK2 c.1343T>G OR 3.03 (95% CI 1.53 to 6.03, p=0.0006) for African men and CHEK2 c.1312G>T OR 2.21 (95% CI 1.06 to 4.63, p=0.030) for European men. No evidence of association with ovarian cancer was found for any of these variants. Conclusions This report adds to accumulating evidence that at least some variants in these genes are associated with an increased risk of breast cancer that is clinically important.

Abstract 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

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.

AbstractGenetic risk scores were used as unconfounded instruments for specific lipid traits (Mendelian randomization) to assess whether circulating lipids causally influence prostate cancer risk. Data from 22,249 prostate cancer cases and 22,133 controls from 22 studies within the international PRACTICAL consortium were analyzed. Allele scores based on single nucleotide polymorphisms (SNPs) previously reported to be uniquely associated with each of low‐density lipoprotein (LDL), high‐density lipoprotein (HDL), and triglyceride (TG) levels, were first validated in an independent dataset, and then entered into logistic regression models to estimate the presence (and direction) of any causal effect of each lipid trait on prostate cancer risk. There was weak evidence for an association between the LDL genetic score and cancer grade: the odds ratio (OR) per genetically instrumented standard deviation (SD) in LDL, comparing high‐ (≥7 Gleason score) versus low‐grade (<7 Gleason score) cancers was 1.50 (95% CI: 0.92, 2.46; P = 0.11). A genetically instrumented SD increase in TGs was weakly associated with stage: the OR for advanced versus localized cancer per unit increase in genetic risk score was 1.68 (95% CI: 0.95, 3.00; P = 0.08). The rs12916‐T variant in 3‐hydroxy‐3‐methylglutaryl‐CoA reductase (HMGCR) was inversely associated with prostate cancer (OR: 0.97; 95% CI: 0.94, 1.00; P = 0.03). In conclusion, circulating lipids, instrumented by our genetic risk scores, did not appear to alter prostate cancer risk. We found weak evidence that higher LDL and TG levels increase aggressive prostate cancer risk, and that a variant in HMGCR (that mimics the LDL lowering effect of statin drugs) reduces risk. However, inferences are limited by sample size and evidence of pleiotropy.

Abstract Background: While numerous susceptibility loci for epithelial ovarian cancer (EOC) have been identified, few associations have been reported with overall survival. In the absence of common prognostic genetic markers, we hypothesize that rare coding variants may be associated with overall EOC survival and assessed their contribution in two exome-based genotyping projects of the Ovarian Cancer Association Consortium (OCAC). Methods: The primary patient set (Set 1) included 14 independent EOC studies (4,293 patients) and 227,892 variants, and a secondary patient set (Set 2) included six additional EOC studies (1,744 patients) and 114,620 variants. Because power to detect rare variants individually is reduced, gene-level tests were conducted. Sets were analyzed separately at individual variants and by gene, and then combined with meta-analyses (73,203 variants and 13,163 genes overlapped). Results: No individual variant reached genome-wide statistical significance. A SNP previously implicated to be associated with EOC risk and, to a lesser extent, survival, rs8170, showed the strongest evidence of association with survival and similar effect size estimates across sets (Pmeta = 1.1E−6, HRSet1 = 1.17, HRSet2 = 1.14). Rare variants in ATG2B, an autophagy gene important for apoptosis, were significantly associated with survival after multiple testing correction (Pmeta = 1.1E−6; Pcorrected = 0.01). Conclusions: Common variant rs8170 and rare variants in ATG2B may be associated with EOC overall survival, although further study is needed. Impact: This study represents the first exome-wide association study of EOC survival to include rare variant analyses, and suggests that complementary single variant and gene-level analyses in large studies are needed to identify rare variants that warrant follow-up study. Cancer Epidemiol Biomarkers Prev; 25(3); 446–54. ©2016 AACR.

ABSTRACTEpithelial‐mesenchymal transition (EMT) is a process whereby epithelial cells assume mesenchymal characteristics to facilitate cancer metastasis. However, EMT also contributes to the initiation and development of primary tumors. Prior studies that explored the hypothesis that EMT gene variants contribute to epithelial ovarian carcinoma (EOC) risk have been based on small sample sizes and none have sought replication in an independent population. We screened 15,816 single‐nucleotide polymorphisms (SNPs) in 296 genes in a discovery phase using data from a genome‐wide association study of EOC among women of European ancestry (1,947 cases and 2,009 controls) and identified 793 variants in 278 EMT‐related genes that were nominally (P < 0.05) associated with invasive EOC. These SNPs were then genotyped in a larger study of 14,525 invasive‐cancer patients and 23,447 controls. A P‐value <0.05 and a false discovery rate (FDR) <0.2 were considered statistically significant. In the larger dataset, GPC6/GPC5 rs17702471 was associated with the endometrioid subtype among Caucasians (odds ratio (OR) = 1.16, 95% CI = 1.07–1.25, P = 0.0003, FDR = 0.19), whereas F8 rs7053448 (OR = 1.69, 95% CI = 1.27–2.24, P = 0.0003, FDR = 0.12), F8 rs7058826 (OR = 1.69, 95% CI = 1.27–2.24, P = 0.0003, FDR = 0.12), and CAPN13 rs1983383 (OR = 0.79, 95% CI = 0.69–0.90, P = 0.0005, FDR = 0.12) were associated with combined invasive EOC among Asians. In silico functional analyses revealed that GPC6/GPC5 rs17702471 coincided with DNA regulatory elements. These results suggest that EMT gene variants do not appear to play a significant role in the susceptibility to EOC.