Zheng, Wei

Abstract Background Although 20 pancreatic cancer susceptibility loci have been identified through genome-wide association studies in individuals of European ancestry, much of its heritability remains unexplained and the genes responsible largely unknown. Methods To discover novel pancreatic cancer risk loci and possible causal genes, we performed a pancreatic cancer transcriptome-wide association study in Europeans using three approaches: FUSION, MetaXcan, and Summary-MulTiXcan. We integrated genome-wide association studies summary statistics from 9040 pancreatic cancer cases and 12 496 controls, with gene expression prediction models built using transcriptome data from histologically normal pancreatic tissue samples (NCI Laboratory of Translational Genomics [n = 95] and Genotype-Tissue Expression v7 [n = 174] datasets) and data from 48 different tissues (Genotype-Tissue Expression v7, n = 74–421 samples). Results We identified 25 genes whose genetically predicted expression was statistically significantly associated with pancreatic cancer risk (false discovery rate < .05), including 14 candidate genes at 11 novel loci (1p36.12: CELA3B; 9q31.1: SMC2, SMC2-AS1; 10q23.31: RP11-80H5.9; 12q13.13: SMUG1; 14q32.33: BTBD6; 15q23: HEXA; 15q26.1: RCCD1; 17q12: PNMT, CDK12, PGAP3; 17q22: SUPT4H1; 18q11.22: RP11-888D10.3; and 19p13.11: PGPEP1) and 11 at six known risk loci (5p15.33: TERT, CLPTM1L, ZDHHC11B; 7p14.1: INHBA; 9q34.2: ABO; 13q12.2: PDX1; 13q22.1: KLF5; and 16q23.1: WDR59, CFDP1, BCAR1, TMEM170A). The association for 12 of these genes (CELA3B, SMC2, and PNMT at novel risk loci and TERT, CLPTM1L, INHBA, ABO, PDX1, KLF5, WDR59, CFDP1, and BCAR1 at known loci) remained statistically significant after Bonferroni correction. Conclusions By integrating gene expression and genotype data, we identified novel pancreatic cancer risk loci and candidate functional genes that warrant further investigation.

Multiple myeloma (MM) is a neoplasm thought to arise from a damaged germinal center B-cell that progresses to a plasma cell clone arising in bone marrow. MM comprises 20% of all hematologic cancer deaths. Persons of African ancestry (AA) have a 1.5 to 2-fold higher risk compared to individuals of European ancestry (EA). Genetically driven differences in hematopoiesis may lead to variation in the levels white blood cell (WBC) subsets which could, in turn, be associated with MM etiology. There are differences in genetic determinants of WBC traits between EA and AA populations, with possible implications for the racial disparity in risk. We tested the above hypothesis using Mendelian randomization (MR), an approach that leverages genetic determinants of specific traits (i.e. WBC counts) to estimate their effects on the risk of an outcome; and a transcriptome-wide association study (TWAS), which utilizes genetic predictors of gene expression to identify susceptibility genes. These analytic approaches were applied to data from the African American Multiple Myeloma Study (AAMMS) consisting of 1813 cases and 8871 AA cancer-free controls to examine how differences in heritable WBC gene expression profiles influence MM risk. Genetic determinants of variation in WBC subsets in AA were obtained from the literature and supplemented with new genome-wide association findings in AA subjects from the UK Biobank cohort (n=6108). Odds ratios (OR) for MM per 1 standard deviation (SD) increase in each WBC phenotype were estimated using independent (linkage disequilibrium (LD) r2<0.10) variants with P<10-6 as genetic instruments. Analyses based on variants associated with WBC traits in AA populations did not identify any statistically-significant associations between MM risk and WBC overall (p=0.81, using 15 SNPs) or subsets (lymphocytes, monocytes, eosinophils, neutrophils and basophils, p>0.05 for each). However, when we applied genetic determinants of WBC identified in 330,000 cancer-free EA UK Biobank participants (P<10-8, replication P<0.05, LD r2<0.05), a statistically significant inverse relationship emerged between increasing lymphocyte counts and MM risk (OR=0.80, 95% CI: 0.66-0.97, p=0.02, using 385 SNPs), as well as increasing basophil counts (OR=0.63, 95% CI: 0.41-0.96, p=0.03, using 140 SNPs). Next, we examined the association between WBC gene expression profiles and MM risk in AAMMS data. We applied published and validated ancestry-specific models developed using the PrediXcan approach, which leverage germline genetic and transcriptomic data from the Multi-Ethnic Study of Atherosclerosis (MESA) (Mogil et al. PMID: 30096133). The primary TWAS used gene expression models trained in AA subjects (n=233), with sensitivity analyses using models developed in AA and Hispanic subjects (n=585). The TWAS significance threshold was based on the number of genes with significant germline prediction models (p<0.05 and R2 ≥0.05) in AA, corresponding to P < 0.05/2700 = 1.85×10-5. The expression of two genes was significantly associated with MM risk: KANK1at 9p24.3 (P = 1.01×10-5) and DNAJC27at 2p23.3 (P = 1.56×10-5). KANK1is a candidate tumor suppressor gene for renal cell carcinoma and has recently been associated with MM risk in AA [Du, Blood 2017 130:3058]. Here we provide additional evidence for its role in MM etiology via gene expression-mediated mechanisms. DNAJC27 (previously known as RBJ) is a novel MM risk gene linked to constitutive activation of ERK in solid tumors. We also identified two suggestively associated genes: PRR14 (P = 1.34×10-4; combined AA-Hispanic sample: P = 1.56×10-6), which has been linked to MM risk in EA populations, and PARP16 (P = 9.46×10-5). To our knowledge this is the first study to comprehensively examine variation in WBC traits and gene expression profiles with respect to MM risk in AA. Our TWAS analysis leveraged data from the largest collection of genetic and gene expression data in AA, enabling ancestry-matched inference and identification of two novel risk genes. Although the limited availability of genetic instruments for WBC limited the power of MR analysis, findings using variants identified in European populations may offer some insight into trans-ethnic etiologic pathways and contribute to risk stratification strategies using genetic and blood cell count biomarkers. Future studies, particularly with MGUS-free controls, are needed to validate these results. Disclosures Song: Millennium Pharmaceuticals Inc: Employment. Rand:Ancestry.com: Employment. Ailawadhi:Cellectar: Research Funding; Janssen: Consultancy, Research Funding; Celgene: Consultancy; Pharmacyclics: Research Funding; Amgen: Consultancy, Research Funding; Takeda: Consultancy. Nooka:Takeda: Honoraria, Other: advisory board participation; Janssen: Honoraria, Other: advisory board participation; GSK: Honoraria, Other: advisory board participation; Spectrum pharmaceuticals: Honoraria, Other: advisory board participation; Adaptive technologies: Honoraria, Other: advisory board participation; Amgen: Honoraria, Other: advisory board participation; Celgene: Honoraria, Other: advisory board participation; BMS: Honoraria, Other: advisory board participation. Singhal:Bureau of Millennium/Takeda, Celgene, Janssen, Celgene, Bristol-Myers Squibb and Bluebird: Speakers Bureau. van Rhee:Takeda: Consultancy; Sanofi Genzyme: Consultancy; Castleman Disease Collaborative Network: Consultancy; EUSA: Consultancy; Adicet Bio: Consultancy; Kite Pharma: Consultancy; Karyopharm Therapeutics: Consultancy. Mehta:Millennium/Takeda, Celgene; stock in Celgene, Bristol-Myers Squibb and Bluebird: Speakers Bureau. Wolf:Takeda: Consultancy; Celgene: Consultancy; Novartis: Consultancy; Janssen: Consultancy; Amgen: Consultancy. Martin:Roche and Juno: Consultancy; Amgen, Sanofi, Seattle Genetics: Research Funding. Fiala:Incyte: Research Funding. Terebelo:Jannsen: Speakers Bureau; Celgene: Honoraria; Newland Medical Asociates: Employment. Anderson:Bristol-Myers Squibb: Other: Scientific Founder; Oncopep: Other: Scientific Founder; Amgen: Consultancy, Speakers Bureau; Janssen: Consultancy, Speakers Bureau; Takeda: Consultancy, Speakers Bureau; Celgene: Consultancy, Speakers Bureau; Sanofi-Aventis: Other: Advisory Board. Vij:Bristol-Myers Squibb: Honoraria, Research Funding; Celgene: Honoraria, Research Funding; Genentech: Honoraria; Janssen: Honoraria; Karyopharm: Honoraria; Sanofi: Honoraria; Takeda: Honoraria, Research Funding. Bernal-Mizrachi:TAKEDA: Research Funding; Kodikas Therapeutic Solutions, Inc: Equity Ownership; Winship Cancer Institute: Employment, Patents & Royalties. Morgan:Amgen, Roche, Abbvie, Takeda, Celgene, Janssen: Honoraria, Membership on an entity's Board of Directors or advisory committees; Celgene: Other: research grant, Research Funding. Zonder:Celgene Corporation: Consultancy, Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; Takeda: Membership on an entity's Board of Directors or advisory committees; Amgen: Consultancy, Membership on an entity's Board of Directors or advisory committees; Intellia: Consultancy, Membership on an entity's Board of Directors or advisory committees; Caelum: Consultancy, Membership on an entity's Board of Directors or advisory committees; Alnylam: Consultancy, Membership on an entity's Board of Directors or advisory committees; BMS: Consultancy, Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; Janssen: Consultancy, Membership on an entity's Board of Directors or advisory committees; Oncopeptides: Consultancy, Membership on an entity's Board of Directors or advisory committees. Huff:Member of Safety Monitoring Board for Johnson and Johnson: Membership on an entity's Board of Directors or advisory committees; Karyopharm, Sanofi, MiDiagnostics: Consultancy. Lonial:Karyopharm: Consultancy; Takeda: Consultancy, Research Funding; Amgen: Consultancy; BMS: Consultancy; Janssen: Consultancy, Research Funding; GSK: Consultancy; Celgene Corporation: Consultancy, Research Funding; Genentech: Consultancy. Orlowski:Poseida Therapeutics, Inc.: Research Funding.

Abstract Multiple myeloma (MM) is 2-3 times more common among African-Americans compared to non-Hispanic whites. The 2-3-fold increased risk among family members of cases suggests a genetic contribution to risk. Genome-wide association studies (GWAS) in populations of European ancestry have identified seven novel risk loci at 2p23.3 (rs6746082), 3q26.2 (rs10936599), 3p22.1 (rs1052501), 6p21.32 (rs2285803), 7p15.3 (rs4487645), 17p11.2 (rs4273077) and 22q13.1 (rs877529) (Broderick, et al. Nat Genet, 2011, Chubb, et al. Nat Genet, 2013), three of which were replicated in another European series (Martino et al., Br J Haematol, 2012). Here we examined the index signals and conducted fine-mapping for each locus in a case-control study of 1,049 multiple myeloma cases and 7,084 controls of African ancestry to identify better markers of risk and novel independent loci in seven previously reported regions in this high risk population. Incident cases were recruited from 10 clinical centers and SEER cancer registries from 2011 to 2013 and genotyped using the Illumina HumanCore GWAS array. Control data were obtained from previous genome-wide studies of breast and prostate cancer, genotyped using the Illumina 1M-Duo in 4425 male controls from the African Ancestry Prostate Cancer Consortium (consisting of 14 independent studies) and 2632 female controls from a breast cancer GWAS of African-American women (consisting of 9 independent studies). Imputation to 1000 Genomes (March 2012 release) was conducted for regions around six of the previously identified single nucleotide polymorphisms [SNPs] (the HLA region harboring rs2285803 is still being imputed, results will be presented). A case-control analysis of SNPs/indels >1% frequency within 250 kb of each index variant was conducted using unconditional multivariable logistic regression adjusting for age, sex and five leading principal components. Region-specific alpha levels were determined through permutation tests. The minimum alpha level across the six regions was α=0.002. All previously reported risk variants were common in African-Americans (minor allele frequency [MAF] >0.05). For five of the six SNPs, we had ≥94% power to detect the same effect observed in non-Hispanic whites, and 64% power for the less common variant rs10936599 (MAF=0.07). We observed directionally consistent effects (odds ratio [OR]>1) for the six risk variants tested, with three replicating at p≤0.05 (7p15.3, p=1.4x10-7; 17p11.2, p=0.05; 22q13.1, p=0.02). For three of the six regions, we observed better markers of risk in African-Americans that were correlated with the index SNP in Europeans (7p15.3, rs56333627, p=1.5x10-5, r2=0.89; 17p11.2, rs34562254, p=2.9x10-3, r2=0.90; 22q13.1, rs2092410, p=1.1x10-4 r2=.71). The missense variant identified in the 17p11.2 region (rs34562254, Pro251Leu) is located in TNFRSF13B, which encodes the protein TACI, a B cell surface receptor which plays a role in B cell maturation, apoptosis and antibody production by inducing activation of transcription factors including NFAT and NFκβ. In addition, there is evidence suggesting that TACI is involved in MM pathogenesis. Our results demonstrate that many of the risk loci for MM found in European ancestry populations are also risk loci in men and women of African ancestry and that by fine-mapping, we are able to identify variants that better capture risk in populations of African ancestry. Disclosures Terebelo: Celgene Corp: Membership on an entity's Board of Directors or advisory committees. Lonial:Millennium: The Takeda Oncology Company: Consultancy, Research Funding; Celgene: Consultancy, Research Funding; Novartis: Consultancy, Research Funding; Bristol-Myers Squibb: Consultancy, Research Funding; Onyx Pharmaceuticals: Consultancy, Research Funding.

Abstract Epithelial ovarian cancer (EOC) is the most common gynecological cancer, and ranks fifth in cancer deaths among women in the United States. Previous studies have identified several circulating proteins biomarkers for EOC risk, such as CA125 antigen and insulin-like growth factor-1 (IGF-1). However, previous studies on protein biomarkers for EOC were limited with small sample sizes, and only a few proteins investigated. In this study, we integrated genomics and proteomics data to search for novel circulating protein biomarkers for EOC risk. We constructed study instruments using genetic variants for protein quantitative trait loci (pQTL) for over 1,400 circulating proteins. Beta coefficients and standard errors of the selected pQTL variants in relation to EOC risk were obtained from 25,509 EOC cases and 40,941 controls from the Ovarian Cancer Association Consortium (OCAC). An inverse-variance weighted method was used to evaluate the association of these genetically predicted protein levels and risk of overall and serous EOC. We identified 35 protein biomarkers, for which their genetically predicted levels were significantly associated with the risk of either overall or serous EOC (false discovery rate &lt;0.05). There are 25 proteins showing significant associations with both overall and serous EOC; most of the associations were stronger for serous EOC than overall EOC. The association for 31 proteins could be attributed, in part, to variants at ABO 9q34.2, a known ovarian cancer susceptibility locus. Five ABO locus-related proteins are in the PI3K-Akt signaling pathway, including insulin-like growth factor 1 receptor, hepatocyte growth factor receptor, interleukin-3 receptor subunit alpha, insulin receptor, and vascular endothelial growth factor receptor 2. Four of these proteins are also in the Rap1 and Ras signaling pathway. These five proteins are all inversely associated with the risk of overall EOC, with odds ratios ranging from 0.64 to 0.93 per unit increase in genetic risk scores (p values: 7.55×10-7 to 5.66×10-2). Four EOC associated proteins which were not influenced by variants in 9q34.2 are microfibrillar-associated protein 2, protein delta homolog 1, semenogelin-2, and netrin-G1, all of which showed significant associations with overall EOC, but not the serous subtype. Our study identifies potential novel biomarkers for ovarian cancer and suggests that genetic variants at locus 9q34.2 may play a role in the relation between levels of these proteins and risk of ovarian cancer. Citation Format: Guochong Jia, Xiang Shu, Xiao-Ou Shu, Jirong Long, Qiuyin Cai, Kirsten B. Moysich, Ellen L. Goode, Thomas A. Sellers, Wei Zheng. Genetically predicted concentration of circulating proteins and epithelial ovarian cancer risk [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 598.

Abstract Pancreatic Ductal Adenocarcinoma (PDAC) is the third leading cause of cancer-related deaths in United States with a 5-year survival rate of only 8%. Inherited predisposition plays an important role in PDAC risk. Rare, moderately to highly penetrant mutations in hereditary cancer and pancreatitis genes, identified in families with a high incidence of disease, account for a small fraction of cases. At the other end of the spectrum, common risk variants of low penetrance have been discovered through genome-wide association studies (GWAS). Currently, over 20 common PDAC risk signals have been identified through GWAS in populations of European ancestry. However, these loci explain only a small fraction of genetic heritability for PDAC, and the causal genes for the large majority of these are unknown. Most susceptibility alleles discovered through GWAS are located in noncoding regions of the genome and are likely to function through allele specific regulation of gene expression. Transcriptome-wide association studies integrate expression quantitative trait loci (eQTL) with GWAS data to discover genes whose cis-regulated expression is associated with complex traits. To identify novel susceptibility loci and likely causal genes, we applied this approach to PDAC. Firstly, we used two reference datasets generated from histologically normal pancreatic tissue samples from individuals of European ancestry, the DCEG/LTG (n=95) and Genotype-Tissue Expression (GTEx v7, n=174), to build robust predictors for gene expression (n=2,872-5,906; r2&gt;0.01), using 5 different models (Top1, Lasso, Enet, Blup and Bslmm) in FUSION and PrediXcan. Gene expression predictors (n=2,043-21,422) were also generated by Enet in 48 different tissues (n=74-421 samples) from 608 individuals of European ancestry (GTEx). We then applied these predictor models to PDAC GWAS summary statistics from 9,040 PC cases and 12,496 controls. We discovered 24 genes whose genetically predicted expression is significantly associated with PDAC risk (FDR &lt; 0.05), including 13 candidate genes at 10 novel loci (CELA3B, SMC2, SMC2-AS1, RP11-80H5.9, SMUG1, BTBD6, RCCD1, PNMT, CDK12, PGAP3, SUPT4H1, RP11-888D10.3 and PGPEP1) and 11 at 5 known risk loci (TERT, CLPTM1L, ZDHHC11B, INHBA, ABO, PDX1, KLF5, WDR59, CFDP1, BCAR1 and TMEM170A). Pathway enrichment analysis for genes differentially expressed in the top vs. bottom quartile of CELA3B expression highlighted higher expression of inflammatory and immune response genes in both the DCEG/LTG and GTEx pancreatic transcriptome datasets, indicating that low expression of this gene may be associated with an inflammatory state in the pancreas. Our findings showcase the utility of integrating multiple datasets to generate hypotheses and provide novel clues for the genetic mechanisms underlying pancreatic cancer risk. Citation Format: Jun Zhong, Ashley Jermusyk, Lang Wu, Jason W. Hoskins, Irene Collins, Mingfeng Zhang, Song Lei, Charles C. Chung, Tongwu Zhang, Wenming Xiao, Rachael Stolzenberg-Solomon, Alison P. Klein, Brian M. Wolpin, Xiao-Ou Shu, Stephen J. Chanock, Sara Olson, Nilanjan Chatterjee, Jill Smith, Jianxin Shi, Peter Kraft, Gloria M. Petersen, Wei Zheng, Laufey T. Amundadottir. Large-scale transcriptome-wide association study (TWAS) identifies novel candidate susceptibility genes for pancreatic cancer [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 1591.

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 &lt; 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 &lt; 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 &lt; 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.

AbstractBackground: Pathogenic germline mutations in the CDKN2A tumor suppressor gene are rare and associated with highly penetrant familial melanoma and pancreatic cancer in non-Hispanic whites (NHW). To date, the prevalence and impact of CDKN2A rare coding variants (RCV) in racial minority groups remain poorly characterized. We examined the role of CDKN2A RCVs on the risk of pancreatic cancer among minority subjects.Methods: We sequenced CDKN2A in 220 African American (AA) pancreatic cancer cases, 900 noncancer AA controls, and 183 Nigerian controls. RCV frequencies were determined for each group and compared with that of 1,537 NHW patients with pancreatic cancer. Odds ratios (OR) and 95% confidence intervals (CI) were calculated for both a case–case comparison of RCV frequencies in AAs versus NHWs, and case–control comparison between AA cases versus noncancer AA controls plus Nigerian controls. Smaller sets of Hispanic and Native American cases and controls also were sequenced.Results: One novel missense RCV and one novel frameshift RCV were found among AA patients: 400G&gt;A and 258_278del. RCV carrier status was associated with increased risk of pancreatic cancer among AA cases (11/220; OR, 3.3; 95% CI, 1.5–7.1; P = 0.004) compared with AA and Nigerian controls (17/1,083). Further, AA cases had higher frequency of RCVs: 5.0% (OR, 13.4; 95% CI, 4.9–36.7; P &lt; 0.001) compared with NHW cases (0.4%).Conclusions: CDKN2A RCVs are more common in AA than in NHW patients with pancreatic cancer and associated with moderately increased pancreatic cancer risk among AAs.Impact: RCVs in CDKN2A are frequent in AAs and are associated with risk for pancreatic cancer. Cancer Epidemiol Biomarkers Prev; 27(11); 1364–70. ©2018 AACR.