Li, Bingshan

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

Abstract Background: Large-scale genome-wide association studies (GWAS) have identified ~35 loci associated with epithelial ovarian cancer risk. However, these genetic loci together explain only a small portion of the heritability of this malignancy. The large majority of the GWAS-identified variants are located in non-coding regions, thus possible causal genes in these loci remain largely unknown. We performed a transcriptome-wide association study (TWAS) to search for novel genetic loci for ovarian cancer risk and plausible causal genes at GWAS-identified loci. Method: Transcriptome data from normal ovarian tissue samples (n=68) and all tissue samples (n = 369), along with their high-density genotyping data, obtained from only European descendants included in the Genotype-Tissue Expression Project (GTEx), were used to build ovarian linear prediction models and cross-tissue models (to increase statistical power) using the elastic net method. Based on model performance, we evaluated 17,121 genes for their cis-predicted gene expressions in relation to ovarian cancer risk using summary statistics data generated in GWAS of ovarian cancer from 25,509 cases and 40,941 controls. MetaXcan was used to integrate gene expression prediction model with summary statistics. Results: We identified 35 genes with predicted expression levels associated with ovarian cancer risk at P value < 2.2 × 10-6, the Bonferroni corrected significance level for multiple comparisons. Of these, 12 genes at 4 genetic loci are located ≥500kb away from risk SNPs previously reported in GWAS, representing potential novel genetic loci for ovarian cancer risk. The remaining 23 genes at 12 loci are located within known ovarian cancer risk loci. Fifteen of these genes at the 12 loci have not been reported in previous studies. Analyses by ovarian cancer histological subtypes showed that the majority of these 35 genes are associated with serous invasive carcinoma. Several new associations were identified in histological subtype analyses. Conclusion: In this TWAS we identified multiple genes with predicted expressions related to ovarian cancer risk and provide substantial new information to enhance the understanding of ovarian cancer biology and genetics. Citation Format: Yingchang Lu, Joellen M. Schildkraut, Thomas A. Sellers, Lang Wu, Xingyi Guo, Bingshan Li, Y. Ann Chen, Jennifer B. Doherty, Simon Gayther, Ellen L. Goode, Hae Kyung Im, Siddhartha Kar, Kate Lawrenson, Ani W. Manichaikul, Jennifer B. Permuth, Brett M. Reid, Jamie K. Teer, Paul Pharoah, Wei Zheng, Jirong Long. Transcriptome-wide association study among 66,450 women to identify candidate susceptible genes for ovarian cancer risk [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1308. doi:10.1158/1538-7445.AM2017-1308