Schleutker, Johanna

Abstract 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.

Abstract 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.

AbstractMotivationThe 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.ResultsWe 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 implementationThe ‘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 informationSupplementary data are available at Bioinformatics online.

ABSTRACTNext‐generation sequencing technologies have afforded unprecedented characterization of low‐frequency and rare genetic variation. Due to low power for single‐variant testing, aggregative methods are commonly used to combine observed rare variation within a single gene. Causal variation may also aggregate across multiple genes within relevant biomolecular pathways. Kernel‐machine regression and adaptive testing methods for aggregative rare‐variant association testing have been demonstrated to be powerful approaches for pathway‐level analysis, although these methods tend to be computationally intensive at high‐variant dimensionality and require access to complete data. An additional analytical issue in scans of large pathway definition sets is multiple testing correction. Gene set definitions may exhibit substantial genic overlap, and the impact of the resultant correlation in test statistics on Type I error rate control for large agnostic gene set scans has not been fully explored. Herein, we first outline a statistical strategy for aggregative rare‐variant analysis using component gene‐level linear kernel score test summary statistics as well as derive simple estimators of the effective number of tests for family‐wise error rate control. We then conduct extensive simulation studies to characterize the behavior of our approach relative to direct application of kernel and adaptive methods under a variety of conditions. We also apply our method to two case‐control studies, respectively, evaluating rare variation in hereditary prostate cancer and schizophrenia. Finally, we provide open‐source R code for public use to facilitate easy application of our methods to existing rare‐variant analysis results.

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

ABSTRACTRare variants (RVs) have been shown to be significant contributors to complex disease risk. By definition, these variants have very low minor allele frequencies and traditional single‐marker methods for statistical analysis are underpowered for typical sequencing study sample sizes. Multimarker burden‐type approaches attempt to identify aggregation of RVs across case‐control status by analyzing relatively small partitions of the genome, such as genes. However, it is generally the case that the aggregative measure would be a mixture of causal and neutral variants, and these omnibus tests do not directly provide any indication of which RVs may be driving a given association. Recently, Bayesian variable selection approaches have been proposed to identify RV associations from a large set of RVs under consideration. Although these approaches have been shown to be powerful at detecting associations at the RV level, there are often computational limitations on the total quantity of RVs under consideration and compromises are necessary for large‐scale application. Here, we propose a computationally efficient alternative formulation of this method using a probit regression approach specifically capable of simultaneously analyzing hundreds to thousands of RVs. We evaluate our approach to detect causal variation on simulated data and examine sensitivity and specificity in instances of high RV dimensionality as well as apply it to pathway‐level RV analysis results from a prostate cancer (PC) risk case‐control sequencing study. Finally, we discuss potential extensions and future directions of this work.