Long, Jirong

Genome-wide association studies (GWAS) have identified more than 200 susceptibility loci for breast cancer, but these variants explain less than a fifth of the disease risk. Although gene–environment interactions have been proposed to account for some of the remaining heritability, few studies have empirically assessed this. We obtained genotype and risk factor data from 46,060 cases and 47,929 controls of European ancestry from population-based studies within the Breast Cancer Association Consortium (BCAC). We built gene expression prediction models for 4,864 genes with a significant (P < 0.01) heritable component using the transcriptome and genotype data from the Genotype-Tissue Expression (GTEx) project. We leveraged predicted gene expression information to investigate the interactions between gene-centric genetic variation and 14 established risk factors in association with breast cancer risk, using a mixed-effects score test. After adjusting for number of tests using Bonferroni correction, no interaction remained statistically significant. The strongest interaction observed was between the predicted expression of the C13orf45 gene and age at first full-term pregnancy (PGXE = 4.44 × 10−6). In this transcriptome-informed genome-wide gene–environment interaction study of breast cancer, we found no strong support for the role of gene expression in modifying the associations between established risk factors and breast cancer risk. Our study suggests a limited role of gene–environment interactions in breast cancer risk.

AbstractBackground:The role of methylation in pancreatic cancer risk remains unclear. We integrated genome and methylome data to identify CpG sites (CpG) with the genetically predicted methylation to be associated with pancreatic cancer risk. We also studied gene expression to understand the identified associations.Methods:Using genetic data and white blood cell methylation data from 1,595 subjects of European descent, we built genetic models to predict DNA methylation levels. After internal and external validation, we applied prediction models with satisfactory performance to the genetic data of 8,280 pancreatic cancer cases and 6,728 controls of European ancestry to investigate the associations of predicted methylation with pancreatic cancer risk. For associated CpGs, we compared their measured levels in pancreatic tumor versus benign tissue.Results:We identified 45 CpGs at nine loci showing an association with pancreatic cancer risk, including 15 CpGs showing an association independent from identified risk variants. We observed significant correlations between predicted methylation of 16 of the 45 CpGs and predicted expression of eight adjacent genes, of which six genes showed associations with pancreatic cancer risk. Of the 45 CpGs, we were able to compare measured methylation of 16 in pancreatic tumor versus benign pancreatic tissue. Of them, six showed differentiated methylation.Conclusions:We identified methylation biomarker candidates associated with pancreatic cancer using genetic instruments and added additional insights into the role of methylation in regulating gene expression in pancreatic cancer development.Impact:A comprehensive study using genetic instruments identifies 45 CpG sites at nine genomic loci for pancreatic cancer risk.

Abstract Background: Whether circulating polyunsaturated fatty acid (PUFA) levels are associated with pancreatic cancer risk is uncertain. Mendelian randomization (MR) represents a study design using genetic instruments to better characterize the relationship between exposure and outcome. Methods: We utilized data from genome-wide association studies within the Pancreatic Cancer Cohort Consortium and Pancreatic Cancer Case–Control Consortium, involving approximately 9,269 cases and 12,530 controls of European descent, to evaluate associations between pancreatic cancer risk and genetically predicted plasma n-6 PUFA levels. Conventional MR analyses were performed using individual-level and summary-level data. Results: Using genetic instruments, we did not find evidence of associations between genetically predicted plasma n-6 PUFA levels and pancreatic cancer risk [estimates per one SD increase in each PUFA-specific weighted genetic score using summary statistics: linoleic acid odds ratio (OR) = 1.00, 95% confidence interval (CI) = 0.98–1.02; arachidonic acid OR = 1.00, 95% CI = 0.99–1.01; and dihomo-gamma-linolenic acid OR = 0.95, 95% CI = 0.87–1.02]. The OR estimates remained virtually unchanged after adjustment for covariates, using individual-level data or summary statistics, or stratification by age and sex. Conclusions: Our results suggest that variations of genetically determined plasma n-6 PUFA levels are not associated with pancreatic cancer risk. Impact: These results suggest that modifying n-6 PUFA levels through food sources or supplementation may not influence risk of pancreatic cancer.

Abstract Pancreatic ductal adenocarcinoma (PDAC) is one of most lethal malignancies with few known risk factors and biomarkers. Identification of disease biomarkers is critical for understanding the pathogenesis of this cancer and identifying high risk individuals for close surveillance. Several blood protein biomarkers have been linked to PDAC in previous studies, but these studies have assessed only a limited number of biomarkers usually in small samples. To identify novel circulating protein biomarkers of PDAC, we studied 8,280 cases and 6,728 controls of European descent from the Pancreatic Cancer Cohort Consortium and the Pancreatic Cancer Case-Control Consortium, by using genetic instruments. Protein quantitative trait loci (pQTLs) for 1,226 plasma proteins identified in a large INTERVAL study of 3,301 healthy European descendants were used as instruments to evaluate associations between genetically predicted protein levels and PDAC. For proteins showing a significant association, we further conducted conditional analysis with adjustments for previously identified risk variants to assess whether the observed associations between genetically predicted protein concentrations and PDAC risk were independent of the risk variants identified in genome-wide association studies (GWAS). Furthermore, for the proteins that were associated with PDAC risk, we performed an enrichment analysis of the genes encoding these proteins to examine whether they are enriched in specific pathways, functions or networks. We observed associations between predicted concentrations of 38 proteins and PDAC risk at a false discovery rate of < 0.05, including those of 23 proteins that showed a significant association even after Bonferroni correction (4.08 × 10−5). These include Histo-blood group ABO system transferase encoded by ABO, which has been previously implicated as a potential target gene of PDAC risk variant identified in GWAS. Eight of the identified proteins (Beta-crystallin B2, Dedicator of cytokinesis protein 9, VIP36-like protein, Erythrocyte band 7 integral membrane protein, Tensin-2, Transmembrane protease serine 11D, Alcohol dehydrogenase 1B, and C-X-C motif chemokine 10) were associated with PDAC risk after conditioning on previously reported pancreatic cancer risk variants (odds ratios ranged from 0.79 to 1.52, P-values from 1.28 × 10−3 to 6.47 × 10−4). Pathway enrichment analysis showed that the encoding genes for the implicated proteins were significantly enriched in cancer-related pathways, such as STAT3 and IL-15 production. In conclusion, we identified 38 protein biomarker candidates for PDAC risk, which if validated by additional studies, may contribute to the etiological understanding of PDAC tumor development. Citation Format: Jingjing Zhu, Xiang Shu, Xingyi Guo, Duo Liu, Jiandong bao, Roger Milne, Graham G Giles, Chong Wu, Mengmeng Du, Emily White, Harvey A Risch, Nuria Malats, Eric J. Duell, Phyllis J. Goodman, Donghui Li, Paige Bracci, Verena Katzke, Rachel E Neale, Steven Gallinger, Stephen Van Den Eeden, Alan Arslan, Federico Canzian, Charles Kooperberg, Brian Wolpin, Laura Beane-Freeman, Ghislaine Scelo, Kala Visvanatha, Christopher A. Haiman, Loïc Le Marchand, Herbert Yu, Gloria M Petersen, Rachael Stolzenberg-Solomon, Alison P Klein, Laufey T Amundadottir, Qiuyin Cai, Jirong Long, Xiao-Ou Shu, Wei Zheng, Lang Wu. Associations between genetically predicted blood protein biomarkers and pancreatic ductal adenocarcinoma risk [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 1200.

Abstract Blood lipids have been associated with the development of a range of cancers, including breast, lung and colorectal cancer. For endometrial cancer, observational studies have reported inconsistent associations between blood lipids and cancer risk. To reduce biases from unmeasured confounding, we performed a bidirectional, two‐sample Mendelian randomization analysis to investigate the relationship between levels of three blood lipids (low‐density lipoprotein [LDL] and high‐density lipoprotein [HDL] cholesterol, and triglycerides) and endometrial cancer risk. Genetic variants associated with each of these blood lipid levels ( P < 5 × 10 −8 ) were identified as instrumental variables, and assessed using genome‐wide association study data from the Endometrial Cancer Association Consortium (12 906 cases and 108 979 controls) and the Global Lipids Genetic Consortium (n = 188 578). Mendelian randomization analyses found genetically raised LDL cholesterol levels to be associated with lower risks of endometrial cancer of all histologies combined, and of endometrioid and non‐endometrioid subtypes. Conversely, higher genetically predicted HDL cholesterol levels were associated with increased risk of non‐endometrioid endometrial cancer. After accounting for the potential confounding role of obesity (as measured by genetic variants associated with body mass index), the association between genetically predicted increased LDL cholesterol levels and lower endometrial cancer risk remained significant, especially for non‐endometrioid endometrial cancer. There was no evidence to support a role for triglycerides in endometrial cancer development. Our study supports a role for LDL and HDL cholesterol in the development of non‐endometrioid endometrial cancer. Further studies are required to understand the mechanisms underlying these findings.

Abstract Background: Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal malignancies, with few known risk factors and biomarkers. Several blood protein biomarkers have been linked to PDAC in previous studies, but these studies have assessed only a limited number of biomarkers, usually in small samples. In this study, we evaluated associations of circulating protein levels and PDAC risk using genetic instruments. Methods: To identify novel circulating protein biomarkers of PDAC, we studied 8,280 cases and 6,728 controls of European descent from the Pancreatic Cancer Cohort Consortium and the Pancreatic Cancer Case-Control Consortium, using genetic instruments of protein quantitative trait loci. Results: We observed associations between predicted concentrations of 38 proteins and PDAC risk at an FDR of &lt; 0.05, including 23 of those proteins that showed an association even after Bonferroni correction. These include the protein encoded by ABO, which has been implicated as a potential target gene of PDAC risk variant. Eight of the identified proteins (LMA2L, TM11D, IP-10, ADH1B, STOM, TENC1, DOCK9, and CRBB2) were associated with PDAC risk after adjusting for previously reported PDAC risk variants (OR ranged from 0.79 to 1.52). Pathway enrichment analysis showed that the encoding genes for implicated proteins were significantly enriched in cancer-related pathways, such as STAT3 and IL15 production. Conclusions: We identified 38 candidates of protein biomarkers for PDAC risk. Impact: This study identifies novel protein biomarker candidates for PDAC, which if validated by additional studies, may contribute to the etiologic understanding of PDAC development.

AbstractPrevious transcriptome‐wide association studies (TWAS) have identified breast cancer risk genes by integrating data from expression quantitative loci and genome‐wide association studies (GWAS), but analyses of breast cancer subtype‐specific associations have been limited. In this study, we conducted a TWAS using gene expression data from GTEx and summary statistics from the hitherto largest GWAS meta‐analysis conducted for breast cancer overall, and by estrogen receptor subtypes (ER+ and ER−). We further compared associations with ER+ and ER− subtypes, using a case‐only TWAS approach. We also conducted multigene conditional analyses in regions with multiple TWAS associations. Two genes, STXBP4 and HIST2H2BA, were specifically associated with ER+ but not with ER– breast cancer. We further identified 30 TWAS‐significant genes associated with overall breast cancer risk, including four that were not identified in previous studies. Conditional analyses identified single independent breast‐cancer gene in three of six regions harboring multiple TWAS‐significant genes. Our study provides new information on breast cancer genetics and biology, particularly about genomic differences between ER+ and ER− breast cancer.

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