Thibodeau, Stephen N.

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.

A cancer-promoting histone protein Mutations in the chromatin protein histone H3 are found in a number of pediatric cancers. The lysine-36–to–methionine (K36M) “oncohistone” mutation is seen in almost all chondroblastomas. Fang et al. show that the K36M mutant histones inhibit the normal methylation of this same residue in wild-type H3 histones. They do so by interfering with the enzymes that normally methylate this residue. The altered chromatin methylation patterns alter the expression of known cancer-related genes and impart cancer-related characteristics to the chondrocyte cells. Science , this issue p. 1344

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.

Abstract African American Hereditary Prostate Cancer Study (AAHPC) was developed as a national collaboration to explore the role of genetics in the causation of hereditary prostate cancer (HPC) in African American (AA) men. AAHPC is in partnership with the International Consortium for Prostate Cancer Genetics (ICPCG), which conducts collaborative studies of HPC genetics in multiplex families. As part of an ICPCG sequencing study of 539 affected individuals from 366 HPC pedigrees, we performed whole exome sequencing in 21 ICPCG AA families, of which there were 14 AAHPC affected men from 11 pedigrees. The combined ICPCG AA cohort consisted of N=26 affected members. Post-variant calling quality control (QC) was implemented using Golden Helix SVS 8 software with filters set for removal of variants with Read Depth &lt; 10, Quality Score &lt; 20, Quality Score: Read Depth Ratio &lt; 0.5, Call Rate &lt; 0.75. Variants were additionally filtered by MAF based on the NHLBI ESP650051-V2 exomes variant frequencies for the AA population using a MAF threshold of 5%. Following QC, 176/299 SNVs and 20/39 INDELs remained for further analysis. In these analyses, we focused on 13 known cancer genes (MSR1, AR, BRCA1, BRCA2, BTNL2, EPHB2, CDH1, RNASEL, ELAC2, HOXB13, CHEK2, TP53 and NBN). Three sequenced families had &gt; 1 affected members sequenced (2 or 3 per family) and the remaining 18 families had one member sequenced. Under the dominant model, our preliminary results show that no rare variants in the 13 candidate genes were found in 3/3 affecteds in two families. Rare SNVs in seven candidate genes (AR, CDH1, ELAC2, HOXB13, RNASEL, BRCA2 and EPHB2) were found in 2/3 affecteds for two families and 2/2 affecteds in one family. Several of the remaining 18 affected men (1 sequenced per pedigree) shared the same rare SNV in these candidate genes. For INDELs, rare variants in three candidate genes were found in pedigrees with ≥ 2 affecteds. Several of the remaining 18 affected men (one sequenced per pedigree) shared the same rare INDEL. Additional QC is underway to validate these variants and bioinformatic analyses are being used to predict effects of the variants in an effort to unravel the complex genetic heterogeneity of HPC in AA. Citation Format: Cheryl D. Cropp, Shannon K. McDonnell, Sumit Middha, Melissa DeRycke, Danielle M. Karyadi, Daniel Schaid, Stephen N. Thibodeau, William B. Isaacs, Elaine A. Ostrander, Janet Stanford, Kathleen A. Cooney, Joan E. Bailey-Wilson, John D. Carpten. Rare variant discovery in known cancer genes from whole-exome sequencingof African American hereditary prostate cancer families. [abstract]. In: Proceedings of the Eighth AACR Conference on The Science of Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; Nov 13-16, 2015; Atlanta, GA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2016;25(3 Suppl):Abstract nr B40.

Germline mutations in the DNA base excision repair gene MUTYH are known to increase a carrier's risk of colorectal cancer. However, the risks of other (extracolonic) cancers for MUTYH mutation carriers are not well defined. We identified 266 probands (91% Caucasians) with a MUTYH mutation (41 biallelic and 225 monoallelic) from the Colon Cancer Family Registry. Mutation status, sex, age and histories of cancer from their 1,903 first‐ and 3,255 second‐degree relatives were analyzed using modified segregation analysis conditioned on the ascertainment criteria. Compared with incidences for the general population, hazard ratios (HRs) (95% confidence intervals [CIs]) for biallelic MUTYH mutation carriers were: urinary bladder cancer 19 (3.7–97) and ovarian cancer 17 (2.4–115). The HRs (95% CI) for monoallelic MUTYH mutation carriers were: gastric cancer 9.3 (6.7–13); hepatobiliary cancer 4.5 (2.7–7.5); endometrial cancer 2.1 (1.1–3.9) and breast cancer 1.4 (1.0–2.0). There was no evidence for an increased risk of cancers at the other sites examined (brain, pancreas, kidney or prostate). Based on the USA population incidences, the estimated cumulative risks (95% CI) to age 70 years for biallelic mutation carriers were: bladder cancer 25% (5–77%) for males and 8% (2–33%) for females and ovarian cancer 14% (2–65%). The cumulative risks (95% CI) for monoallelic mutation carriers were: gastric cancer 5% (4–7%) for males and 2.3% (1.7–3.3%) for females; hepatobiliary cancer 3% (2–5%) for males and 1.4% (0.8–2.3%) for females; endometrial cancer 3% (2%–6%) and breast cancer 11% (8–16%). These unbiased estimates of both relative and absolute risks of extracolonic cancers for people, mostly Caucasians, with MUTYH mutations will be important for their clinical management.

Abstract Identifying genetic variants with pleiotropic associations can uncover common pathways influencing multiple cancers. We took a two-stage approach to conduct genome-wide association studies for lung, ovary, breast, prostate, and colorectal cancer from the GAME-ON/GECCO Network (61,851 cases, 61,820 controls) to identify pleiotropic loci. Findings were replicated in independent association studies (55,789 cases, 330,490 controls). We identified a novel pleiotropic association at 1q22 involving breast and lung squamous cell carcinoma, with eQTL analysis showing an association with ADAM15/THBS3 gene expression in lung. We also identified a known breast cancer locus CASP8/ALS2CR12 associated with prostate cancer, a known cancer locus at CDKN2B-AS1 with different variants associated with lung adenocarcinoma and prostate cancer, and confirmed the associations of a breast BRCA2 locus with lung and serous ovarian cancer. This is the largest study to date examining pleiotropy across multiple cancer-associated loci, identifying common mechanisms of cancer development and progression. Cancer Res; 76(17); 5103–14. ©2016 AACR.

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: Unnecessary intervention and overtreatment of indolent disease are common challenges in clinical management of prostate cancer. Improved tools to distinguish lethal from indolent disease are critical. Methods: We performed a genome-wide survival analysis of cause-specific death in 24,023 prostate cancer patients (3,513 disease-specific deaths) from the PRACTICAL and BPC3 consortia. Top findings were assessed for replication in a Norwegian cohort (CONOR). Results: We observed no significant association between genetic variants and prostate cancer survival. Conclusions: Common genetic variants with large impact on prostate cancer survival were not observed in this study. Impact: Future studies should be designed for identification of rare variants with large effect sizes or common variants with small effect sizes. Cancer Epidemiol Biomarkers Prev; 24(11); 1796–800. ©2015 AACR.