Hildebrandt, Michelle A.T.

Background Genome-wide association studies (GWAS) conducted in populations of European ancestry (EA) have identified and confirmed 23 germline susceptibility loci for multiple myeloma (MM). The effect sizes of single nucleotide polymorphisms (SNPs) at these loci are small, therefore combining them into a single summary measure, known as a polygenic risk score (PRS), may provide a more meaningful risk factor. We have previously shown a PRS comprised of the 23 SNPs for MM contributes to increased risk of MM, with a 2.7-fold increase for highest vs. lowest PRS quintiles. Whether the MM-PRS is also associated with overall survival (OS) in MM cases has not been evaluated. We examined the association between MM-PRS and OS in two EA studies. Methods The first study consisted of 2,179 EA MM cases from ten studies included in the Multiple Myeloma Working Group within the International Lymphoma Consortium (InterLymph). Cases were diagnosed between 1970 and 2015 and genotyped using multiple platforms (Oncoarray, Affymetrix, Human660W-quad Beadchip, and Illumina arrays); 885 cases also had stage [based on International Staging System (ISS)] available. Each of the GWAS was subjected to rigorous standard quality control independently (prior to imputation via the Michigan imputation server based on the Haplotype Reference Consortium (HRC). The second study consisted of 515 newly diagnosed EA MM cases from CoMMpass (Relating Clinical Outcomes in Multiple Myeloma to Personal Assessment of Genetic Profile), diagnosed from 2011-2013, who had whole genome sequencing (WGS) performed on germline DNA. The WGS data was used to call common germline genetic variants through the Mayo Clinic bioinformatics pipeline. Briefly, genetic variants were detected with GenomeGPS, aligned to the hg19 reference genome, called using the GATK (V3.6) Haplotype Caller, and merged for multiple-sample joint calling. To reduce the false positive variants, variant quality score recalibration (VQSR) was applied for both SNPs and INDELs. After quality control, 458 EA samples remained. Follow-up was available for both studies and consisted of time from MM diagnosis date until death or date of last known follow-up. The PRS was constructed from the 23 MM SNPs using the published per allele odds ratio associated with MM risk. The published log odds ratios for each SNP were multiplied by the number of risk alleles (0, 1, 2) for the corresponding SNP, and summed, resulting in a unique score per person. Kaplan-Meier curves and Cox proportional hazard models were used to assess the association between PRS with MM OS considering two models: 1) adjusted for age, sex, study and 2) additional adjustment by stage (ISS). Hazard ratios (OR) and 95% confidence intervals (CI) were estimated. The PRS was evaluated both as a continuous variable, per standard deviation (SD), and as a categorical variable (quintiles). Results MM cases (N=2,179) in the InterLymph study were 59% male and 41% female and the median age was 61.0 years (26-90 years). Median follow-up time was 57.2 months (1.0-509.0 months) with 868 reported deaths. MM cases with stage information available consisted of 20% stage I (n=178), 53% stage II (n=466), and 27% stage III (n=241). No association was observed between PRS and OS in MM patients regardless of adjustment for stage (continuous PRS (HR: 1.03, 95% CI: 0.83-1.28, P=0.80) or by quintile PRS (p>0.05)) (Table). The CoMMpass EA MM cases (n=458) had similar distributions for sex (61% male and 39% females) but were slightly older 65 years (27-93 years) and had shorter follow-up time (median=39.75 months (0.13-77.2)) with 117 deaths. Stage was available for 96% of CoMMpass cases including 36% stage I (n=159), 33% stage II (n=146), and 31% stage III (n=134). We also observed no association of PRS and OS in the CoMMpass study (HR=1.02, 95% CI: 0.72 -1.46, P= 0.89), adjusted for age, sex, and stage (Table). Discussion A PRS score for MM risk is not associated with OS for MM cases in two EA populations. Given that prior studies have shown association of genetic variation with MM survival, efforts to identify additional loci associated with OS or MM specific survival are warranted. Future studies should also consider germline variants impact on molecular subtypes, specific therapies, and outcomes. Disclosures Kumar: Celgene: Consultancy, Research Funding; Janssen: Consultancy, Research Funding; Takeda: Research Funding.

Abstract We previously identified five single nucleotide polymorphisms (SNPs) at four susceptibility loci for diffuse large B-cell lymphoma (DLBCL) in individuals of European ancestry through a large genome-wide association study (GWAS). To further elucidate genetic susceptibility to DLBCL, we sought to validate two loci at 3q13.33 and 3p24.1 that were suggestive in the original GWAS with additional genotyping. In the meta-analysis (5662 cases and 9237 controls) of the four original GWAS discovery scans and three replication studies, the 3q13.33 locus (rs9831894; minor allele frequency [MAF] = 0.40) was associated with DLBCL risk [odds ratio (OR) = 0.83, P = 3.62 × 10−13]. rs9831894 is in linkage disequilibrium (LD) with additional variants that are part of a super-enhancer that physically interacts with promoters of CD86 and ILDR1. In the meta-analysis (5510 cases and 12 817 controls) of the four GWAS discovery scans and four replication studies, the 3p24.1 locus (rs6773363; MAF = 0.45) was also associated with DLBCL risk (OR = 1.20, P = 2.31 × 10−12). This SNP is 29 426-bp upstream of the nearest gene EOMES and in LD with additional SNPs that are part of a highly lineage-specific and tumor-acquired super-enhancer that shows long-range interaction with AZI2 promoter. These loci provide additional evidence for the role of immune function in the etiology of DLBCL, the most common lymphoma subtype.

Abstract Genome-wide association studies (GWAS) conducted among populations of European ancestry (EA) have identified 23 common single nucleotide polymorphisms (SNPs) associated with multiple myeloma (MM) risk. We hypothesize that the combination of these SNPs in a polygenic risk score (PRS) is likely to be a strong risk factor for MM. However, it is unclear whether the genetic variation associated with MM susceptibility also predisposes to monoclonal gammopathy of undetermined significance (MGUS). Thus, we calculated a PRS and evaluated the association with risk of MM and its precursor, MGUS. We pooled genotype data for 2434 MM and 3446 controls from ten MM GWAS of individuals of EA within the Interlymph Consortium, for 23 MM risk SNPs identified by prior GWAS. An additional 754 MGUS cases were ascertained from Mayo Clinic and MD Anderson clinical practices. To calculate the PRS, we used the risk estimates corresponding to the 23 SNP associations from the largest published MM GWAS. The log of the odds ratio (OR) for each SNP was multiplied by the respective number of risk alleles and summed to generate a PRS for each individual. The PRS was examined continuously, per one standard deviation (SD), and as quintiles, based on the PRS distribution in the controls. Associations of PRS with MM and MGUS risk were examined separately, using multivariable logistic regression assuming an additive model to assess ORs and 95% confidence intervals adjusted for age, sex, and site. We also evaluated age and sex stratified models. The distribution of sex within MM cases, MGUS cases and controls were each ~60% male and ~40% female. The median age was 61, 66, and 66 years for MM cases, MGUS cases and controls, respectively. PRS ranged from 1.52-4.91, with a median PRS of 3.21 for MM cases, 3.19 for MGUS cases, and 3.05 for controls. PRS was significantly associated with MM risk when assessed continuously (OR=1.19 per SD, p=2.2x10-16) and categorically; compared with the middle quintile (Q3), individuals in the highest quintile (Q5) had a 66% increased MM risk (OR=1.66, p=2.3x10-9) and those in the lowest quintile (Q1) had a 38% decreased MM risk (OR=0.62, p=1.3x10-6). PRS was also significantly associated with MGUS risk (OR=1.19 per SD, p=1.7x10-11); individuals with the highest PRS (Q5) had a 77% increased risk (OR=1.77, p=4.0x10-4) and those with lowest PRS (Q1) had 30% decreased risk (OR=0.70, p=0.04), compared with Q3. When stratified by age and sex, similar associations and trends were found. Using an independent sample of MM / MGUS cases and controls, we showed that a PRS constructed from 23 common genetic variants for MM risk is associated with risk of both MM and MGUS, regardless of age or sex. A future direction of this work is testing associations with PRS and clinical characteristics of the MM cases, as well as differences between MGUS cases that progress and those that do not. Our results suggest that common genetic variation may predispose to MGUS as the precursor to MM. Citation Format: Alyssa I. Clay-Gilmour, Michelle A. Hildebrandt, Nicola J. Camp, Elad Ziv, Elizabeth E. Brown, Jonathan N. Hofmann, John J. Spinelli, Graham G. Giles, Parveen Bhatti, Wendy Cozen, Xifeng Wu, Dennis P. Robinson, Aaron D. Norman, Jason P. Sinnwell, Shaji K. Kumar, S Vincent Rajkumar, Susan L. Slager, Celine M. Vachon. Associations between a polygenic risk score and risk of multiple myeloma and its precursor [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 2686.

Abstract DNA methylation is instrumental for gene regulation. Global changes in the epigenetic landscape have been recognized as a hallmark of cancer. However, the role of DNA methylation in epithelial ovarian cancer (EOC) remains unclear. In this study, high-density genetic and DNA methylation data in white blood cells from the Framingham Heart Study (N = 1,595) were used to build genetic models to predict DNA methylation levels. These prediction models were then applied to the summary statistics of a genome-wide association study (GWAS) of ovarian cancer including 22,406 EOC cases and 40,941 controls to investigate genetically predicted DNA methylation levels in association with EOC risk. Among 62,938 CpG sites investigated, genetically predicted methylation levels at 89 CpG were significantly associated with EOC risk at a Bonferroni-corrected threshold of P < 7.94 × 10−7. Of them, 87 were located at GWAS-identified EOC susceptibility regions and two resided in a genomic region not previously reported to be associated with EOC risk. Integrative analyses of genetic, methylation, and gene expression data identified consistent directions of associations across 12 CpG, five genes, and EOC risk, suggesting that methylation at these 12 CpG may influence EOC risk by regulating expression of these five genes, namely MAPT, HOXB3, ABHD8, ARHGAP27, and SKAP1. We identified novel DNA methylation markers associated with EOC risk and propose that methylation at multiple CpG may affect EOC risk via regulation of gene expression. Significance: Identification of novel DNA methylation markers associated with EOC risk suggests that methylation at multiple CpG may affect EOC risk through regulation of gene expression.

Abstract Large-scale genome-wide association studies (GWAS) have identified approximately 35 loci associated with epithelial ovarian cancer (EOC) risk. The majority of GWAS-identified disease susceptibility variants are located in noncoding regions, and causal genes underlying these associations remain largely unknown. Here, we performed a transcriptome-wide association study to search for novel genetic loci and plausible causal genes at known GWAS loci. We used RNA sequencing data (68 normal ovarian tissue samples from 68 individuals and 6,124 cross-tissue samples from 369 individuals) and high-density genotyping data from European descendants of the Genotype-Tissue Expression (GTEx V6) project to build ovarian and cross-tissue models of genetically regulated expression using elastic net methods. We evaluated 17,121 genes for their cis-predicted gene expression in relation to EOC risk using summary statistics data from GWAS of 97,898 women, including 29,396 EOC cases. With a Bonferroni-corrected significance level of P < 2.2 × 10−6, we identified 35 genes, including FZD4 at 11q14.2 (Z = 5.08, P = 3.83 × 10−7, the cross-tissue model; 1 Mb away from any GWAS-identified EOC risk variant), a potential novel locus for EOC risk. All other 34 significantly associated genes were located within 1 Mb of known GWAS-identified loci, including 23 genes at 6 loci not previously linked to EOC risk. Upon conditioning on nearby known EOC GWAS-identified variants, the associations for 31 genes disappeared and three genes remained (P < 1.47 × 10−3). These data identify one novel locus (FZD4) and 34 genes at 13 known EOC risk loci associated with EOC risk, providing new insights into EOC carcinogenesis. Significance: Transcriptomic analysis of a large cohort confirms earlier GWAS loci and reveals FZD4 as a novel locus associated with EOC risk. Cancer Res; 78(18); 5419–30. ©2018 AACR.

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

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.

Abstract Background: Genome-wide association studies (GWAS) have so far reported 12 loci associated with serous epithelial ovarian cancer (EOC) risk. We hypothesized that some of these loci function through nearby transcription factor (TF) genes and that putative target genes of these TFs as identified by coexpression may also be enriched for additional EOC risk associations. Methods: We selected TF genes within 1 Mb of the top signal at the 12 genome-wide significant risk loci. Mutual information, a form of correlation, was used to build networks of genes strongly coexpressed with each selected TF gene in the unified microarray dataset of 489 serous EOC tumors from The Cancer Genome Atlas. Genes represented in this dataset were subsequently ranked using a gene-level test based on results for germline SNPs from a serous EOC GWAS meta-analysis (2,196 cases/4,396 controls). Results: Gene set enrichment analysis identified six networks centered on TF genes (HOXB2, HOXB5, HOXB6, HOXB7 at 17q21.32 and HOXD1, HOXD3 at 2q31) that were significantly enriched for genes from the risk-associated end of the ranked list (P &lt; 0.05 and FDR &lt; 0.05). These results were replicated (P &lt; 0.05) using an independent association study (7,035 cases/21,693 controls). Genes underlying enrichment in the six networks were pooled into a combined network. Conclusion: We identified a HOX-centric network associated with serous EOC risk containing several genes with known or emerging roles in serous EOC development. Impact: Network analysis integrating large, context-specific datasets has the potential to offer mechanistic insights into cancer susceptibility and prioritize genes for experimental characterization. Cancer Epidemiol Biomarkers Prev; 24(10); 1574–84. ©2015 AACR.