Citations

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.

AbstractBackground and AimTumor testing of colorectal cancers (CRC) for mismatch repair (MMR) deficiency is an effective approach to identify carriers of germline MMR gene mutation (Lynch syndrome). The aim of this study was to identify MMR gene mutation carriers in two cohorts of population‐based CRC utilizing a combination of tumor and germline testing approaches.MethodsColorectal cancers from 813 patients diagnosed with CRC < 60 years of age from the Australasian Colorectal Cancer Family Registry (ACCFR) and from 826 patients from the Melbourne Collaborative Cohort Study (MCCS) were tested for MMR protein expression using immunohistochemistry, microsatellite instability (MSI), BRAFV600E somatic mutation, and for MLH1 methylation. MMR gene mutation testing (Sanger sequencing and Multiplex Ligation Dependent Probe Amplification) was performed on germline DNA of patients with MMR‐deficient tumors and a subset of MMR‐proficient CRCs.ResultsOf the 813 ACCFR probands, 90 probands demonstrated tumor MMR deficiency (11.1%), and 42 had a MMR gene germline mutation (5.2%). For the MCCS, MMR deficiency was identified in the tumors of 103 probands (12.5%) and seven had a germline mutation (0.8%). All the mutation carriers were diagnosed prior to 70 years of age. Probands with a MMR‐deficient CRC without MLH1 methylation and a gene mutation were considered Lynch‐like and comprised 41.1% and 25.2% of the MMR‐deficient CRCs for the ACCFR and MCCS, respectively.ConclusionsIdentification of MMR gene mutation carriers in Australian CRC‐affected patients is optimized by immunohistochemistry screening of CRC diagnosed before 70 years of age. A significant proportion of MMR‐deficient CRCs will have unknown etiology (Lynch‐like) proving problematic for clinical management.

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.

Background BRCA1 interacting protein C-terminal helicase 1 (BRIP1) is one of the Fanconi Anaemia Complementation (FANC) group family of DNA repair proteins. Biallelic mutations in BRIP1 are responsible for FANC group J, and previous studies have also suggested that rare protein truncating variants in BRIP1 are associated with an increased risk of breast cancer. These studies have led to inclusion of BRIP1 on targeted sequencing panels for breast cancer risk prediction. Methods We evaluated a truncating variant, p.Arg798Ter (rs137852986), and 10 missense variants of BRIP1, in 48 144 cases and 43 607 controls of European origin, drawn from 41 studies participating in the Breast Cancer Association Consortium (BCAC). Additionally, we sequenced the coding regions of BRIP1 in 13 213 cases and 5242 controls from the UK, 1313 cases and 1123 controls from three population-based studies as part of the Breast Cancer Family Registry, and 1853 familial cases and 2001 controls from Australia. Results The rare truncating allele of rs137852986 was observed in 23 cases and 18 controls in Europeans in BCAC (OR 1.09, 95% CI 0.58 to 2.03, p=0.79). Truncating variants were found in the sequencing studies in 34 cases (0.21%) and 19 controls (0.23%) (combined OR 0.90, 95% CI 0.48 to 1.70, p=0.75). Conclusions These results suggest that truncating variants in BRIP1, and in particular p.Arg798Ter, are not associated with a substantial increase in breast cancer risk. Such observations have important implications for the reporting of results from breast cancer screening panels.

Individuals diagnosed with colorectal cancer (CRC) are at risk of developing a metachronous CRC. We examined the associations between personal, tumour‐related and lifestyle risk factors, and risk of metachronous CRC. A total of 7,863 participants with incident colon or rectal cancer who were recruited in the USA, Canada and Australia to the Colon Cancer Family Registry during 1997‐2012, except those identified as high‐risk, for example, Lynch syndrome, were followed up approximately every 5 years. We estimated the risk of metachronous CRC, defined as the first new primary CRC following an interval of at least one year after the initial CRC diagnosis. Observation time started at the age at diagnosis of the initial CRC and ended at the age at diagnosis of the metachronous CRC, last contact or death whichever occurred earliest, or were censored at the age at diagnosis of any metachronous colorectal adenoma. Cox regression was used to derive hazard ratios (HRs) and 95% confidence intervals (CIs). During a mean follow‐up of 6.6 years, 142 (1.81%) metachronous CRCs were diagnosed (mean age at diagnosis 59.8; incidence 2.7/1,000 person‐years). An increased risk of metachronous CRC was associated with the presence of a synchronous CRC (HR = 2.73; 95% CI: 1.30–5.72) and the location of cancer in the proximal colon at initial diagnosis (compared with distal colon or rectum, HR = 4.16; 95% CI: 2.80–6.18). The presence of a synchronous CRC and the location of the initial CRC might be useful for deciding the intensity of surveillance colonoscopy for individuals diagnosed with CRC.

Previous genome‐wide association studies among women of European ancestry identified two independent breast cancer susceptibility loci represented by single nucleotide polymorphisms (SNPs) rs13281615 and rs11780156 at 8q24. A fine‐mapping study across 2.06 Mb (chr8:127,561,724–129,624,067, hg19) in 55,540 breast cancer cases and 51,168 controls within the Breast Cancer Association Consortium was conducted. Three additional independent association signals in women of European ancestry, represented by rs35961416 (OR = 0.95, 95% CI = 0.93–0.97, conditional p = 5.8 × 10−6), rs7815245 (OR = 0.94, 95% CI = 0.91–0.96, conditional p = 1.1 × 10−6) and rs2033101 (OR = 1.05, 95% CI = 1.02–1.07, conditional p = 1.1 × 10−4) were found. Integrative analysis using functional genomic data from the Roadmap Epigenomics, the Encyclopedia of DNA Elements project, the Cancer Genome Atlas and other public resources implied that SNPs rs7815245 in Signal 3, and rs1121948 in Signal 5 (in linkage disequilibrium with rs11780156, r2 = 0.77), were putatively functional variants for two of the five independent association signals. The results highlighted multiple 8q24 variants associated with breast cancer susceptibility in women of European ancestry.

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.