Lawrenson, Kate

Abstract Introduction: In females, X-chromosome inactivation (XCI) epigenetically silences transcription of one copy of the X chromosome. Which chromosome is silenced is randomly selected, and is tissue- and cell-specific. While some genes are known to escape XCI under normal conditions, aberrant XCI patterns are thought to occur in female-specific cancers, although the role of XCI in ovarian tumorigenesis and progression is largely unknown. The process of XCI is complex, and integration of gene expression, DNA methylation, and copy number data can inform the XCI status of individual genes and chromosome-wide XCI patterns for individual patients. Methods: We evaluated gene- and chromosome-level patterns of XCI by integrating RNA sequence, copy number alteration, genotype, and DNA methylation data to study XCI escape patterns in tumor samples from 99 ovarian cancer patients. We measured allele-specific expression (ASE) for 397 X-linked genes to identify the active alleles for each tumor. Combining ASE data with knowledge of copy number status, we used a Bayesian beta-binomial mixture model to estimate which genes escaped XCI for each patient, and validated our findings using DNA methylation data. To assess global XCI patterns, we performed cluster analyses on the ASE and methylation data, after adjusting for loss of heterozygosity. We examined the relationship between the clusters and clinical factors, including overall survival and time to recurrence. Results: DNA promoter methylation demonstrated inverse regional correlations with ASE. Cluster analyses using ASE and methylation data demonstrated evidence of two tumor clusters, representing normal XCI and global XCI dysregulation. The dysregulated XCI cluster (N=52) was associated with lower X-inactive specific transcript expression as expected (p<0.01). Patients with XCI dysregulated tumors were higher grade, stage, serous histology and were sub-optimally debulked (p<0.05). These patients also had shorter overall survival time (HR=1.87, p=0.02) and time to recurrence (HR=2.34, p<0.01), although associations were attenuated after covariate adjustment. In 45 tumor samples with sufficient data, we observed escape patterns largely consistent with previous reports of multiple tissue types. When comparing tumor to normal ovarian tissue, eight genes (CXorf23, CXorf36, BRWD3, ELF4, SLITRK4, GABRE, CLCN4, SH3BGRL) showed putative escape in the tumor and two genes (RBBP7, OFD1) showed discrepant tumor inactivation. Conclusions: We identified discrepant gene-level XCI tumor classifications compared to normal tissue and identified a group of patients with chromosome-wide XCI dysregulation associated with worse clinical prognosis. This provides evidence of the role of XCI in ovarian cancer and highlights the need to integrate multiple genomic data types to study XCI. Citation Format: Stacey J. Winham, Nicholas B. Larson, Sebastian M. Armasu, Zachary C. Fogarty, Melissa C. Larson, Kimberly R. Kalli, Kate Lawrenson, Simon Gayther, Brooke L. Fridley, Ellen L. Goode. Integrative analyses of gene expression, DNA methylation, genotype and copy number alterations characterize X-chromosome inactivation in ovarian cancer [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 2420. doi:10.1158/1538-7445.AM2017-2420

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

Abstract Epithelial ovarian cancer (EOC) is a rare but deadly disease for which there is notably poorer survival in women of African Ancestry (AA) compared to women of European Ancestry (EA). Previous EA-based genome-wide association studies (GWAS) have identified 30 common, low penetrant EOC susceptibility alleles. Using the custom-designed 533,631 SNP Illumina OncoArray and imputation to ~12 million genetic variants in the 1000 Genomes Phase 3, we conducted a GWAS in 755 AA EOC cases, including 537 high-grade serous ovarian cancer (HGSOC) cases, and 1,235 AA controls from the Ovarian Cancer Association Consortium (OCAC). We identified novel susceptibility loci with suggestive evidence of association with EOC, based on a threshold of P<1x10-6, at three distinct loci including 10p15.1 (lead SNP rs4525119, intronic to AKR1C3, P=4.9 x 10-7, minor allele frequency [MAF]=0.33), 3p25.3 (lead SNP rs7643459, intronic to LOC101927394, P=8.4 x 10-7, MAF=0.36), and 4q13.3 (lead SNP rs4286604, 15 kb 3’ of UGT2A1, P=8.5x10-7, MAF=0.27). In analysis of HGSOC, we identified suggestive evidence of association at six distinct loci including 5q11.2 (lead SNP rs37792, 132 kb 5’ of FST, P=6.0x10-8, MAF=0.34), Xq27.2 (lead SNP rs57403204, 81 kb 3' of MAGEC1, P=1.7x10-7, MAF=0.06), 10p15.1 (lead SNP rs79079890, LOC105376360 intronic, P=3.0x10-7, MAF=0.03), 17p25.1 (lead SNP rs66459581, 5 kb 5’ of PRPSAP1, P=5.1x10-7, MAF=0.23), 15p12 (lead SNP rs116046250, GABRG3 intronic, P=8.7x10-7, MAF=0.05), and 4q21.21 (lead SNP rs192876988, 30 kb 3’ of GK2, P=9.2x10-7, MAF=0.05). These SNPs showed no evidence of association with risk in the OCAC EA GWAS including up to 23,543 EOC cases and 29,444 controls. Of the SNPs identified with suggestive evidence of association in our GWAS, two are located near genes known to regulate hormones and previously reported in relation to diseases of the ovary (AKR1C3 and FST), two have been linked to cancer (AKR1C3 and MAGEC1) and one represents testis-specific expression (GK2). Of the 30 SNPs identified previously in a EA EOC GWAS, we observed nominally significant associations (P<0.05) with consistent direction of effect for eight SNPs in AA women. Our investigation presents evidence of (1) variants for EOC shared among EA and AA women and (2) novel EOC risk loci of importance for AA women. Citation Format: Ani Manichaikul, Lauren C. Peres, Xin-Qun Wang, Kate Lawrenson, Sarah Abbott, Ann G. Schwartz, Anna H. Wu, Edward Peters, Patricia G. Moorman, Michele L. Cote, Melissa Bondy, Linda Kelemen, Ellen L. Goode, Jill Barnholtz-Sloan, Simon A. Gather, Andrew Berchuck, Jennifer A. Doherty, Paul Pharoah, Joellen Schildkraut. Identification of novel epithelial ovarian cancer loci in women of African ancestry from the Ovarian Cancer Association Consortium [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 1307. doi:10.1158/1538-7445.AM2017-1307

Abstract Background: Genome-wide association studies (GWAS) have identified multiple loci associated with epithelial ovarian cancer (EOC) susceptibility, but further progress requires integration of epidemiology and biology to illuminate true risk loci below genome-wide significance levels (P < 5 × 10−8). Most risk SNPs lie within non–protein-encoding regions, and we hypothesize that long noncoding RNA (lncRNA) genes are enriched at EOC risk regions and represent biologically relevant functional targets. Methods: Using imputed GWAS data from about 18,000 invasive EOC cases and 34,000 controls of European ancestry, the GENCODE (v19) lncRNA database was used to annotate SNPs from 13,442 lncRNAs for permutation-based enrichment analysis. Tumor expression quantitative trait locus (eQTL) analysis was performed for sub-genome-wide regions (1 × 10−5 > P > 5 × 10−8) overlapping lncRNAs. Results: Of 5,294 EOC-associated SNPs (P < 1.0 × 10−5), 1,464 (28%) mapped within 53 unique lncRNAs and an additional 3,484 (66%) SNPs were correlated (r2 > 0.2) with SNPs within 115 lncRNAs. EOC-associated SNPs comprised 130 independent regions, of which 72 (55%) overlapped with lncRNAs, representing a significant enrichment (P = 5.0 × 10−4) that was more pronounced among a subset of 5,401 lncRNAs with active epigenetic regulation in normal ovarian tissue. EOC-associated lncRNAs and their putative promoters and transcription factors were enriched for biologically relevant pathways and eQTL analysis identified five novel putative risk regions with allele-specific effects on lncRNA gene expression. Conclusions: lncRNAs are significantly enriched at EOC risk regions, suggesting a mechanistic role for lncRNAs in driving predisposition to EOC. Impact: lncRNAs represent key candidates for integrative epidemiologic and functional studies. Further research on their biologic role in ovarian cancer is indicated. Cancer Epidemiol Biomarkers Prev; 26(1); 116–25. ©2016 AACR.

Abstract Background: Little transcriptomic research has compared epithelial ovarian cancer (EOC) histological subtypes. We set out to characterize the transcriptomes of high-grade serous carcinomas (HGSC) and endometrioid carcinomas (EC), which make up around 70% and 20% of EOC tumors, respectively, and have some histopathological similarities. Methods: Fresh frozen tumors from EOC patients seen at the Mayo Clinic (30 EC and 62 HGSC) were used. 1ug RNA riboZero was used for library preparation using the Illumina TruSeq kit and sequenced on a HiSeq 2000 machine. Reads were aligned using TopHat2 followed by quantification of abundances using RSEM and differential expression analysis with edgeR. We analyzed transcriptomes, conducted pathway analyses, and summarized key candidate gene sets. Expressed SNVs (eSNVs) from the RNA-seq data were determined using GATK and RVboost. Results: The analysis found 699 genes with FDR < 1×10-5 for differential expression between HGSC and EC, with most genes being up-regulated in EC. The top most associated genes were TPH1, MAP2K6, KLK2, ADAM23, TESC and TRAF3IP2 (p<10-22). Pathway analysis of the genes up-regulated in EC revealed enrichment of the “basal cell carcinoma signaling pathway” (p = 1.2×10-5). Within 1 Mb of the 25 known EOC risk loci, we observed higher expression in HGSC for RSPO1 and HPSE (p <5×10-7). For genes functionally related to EOC, we observed in HGSC up-regulation (*p < 10-5, ⁁p<0.003) for FOXM1⁁, CDKN2A⁁, CCNE1*, CCND2⁁, PIK3CA*, BRCA2⁁, BIRC5⁁, MMP9⁁, FANCD2⁁, and MAML2⁁. In contrast, we found up-regulation in EC for MDM2*, KLK4*, BCL2⁁, CCND1*, ANXA4*, CDH1⁁, MMP7⁁, and MAML3⁁. We also identified 204 eSNVs (44 non-synonymous) associated with EC v HGSC subtype (p<10-4); this included an exonic TRAF3IP2 eSNV (66% EC, 13% HGSC, p = 4×10-7, chr6:111877117). Discussion: Using one of the largest sets of identically processed fresh-frozen EOC tumors, some patterns emerged among the numerous EC v HGSC transcriptomic differences. TPH1, up-expressed in EC, is regulated by SOX4 which was also up-regulated in EC. Two sets of genes related to Kallikreins serine proteases were differentially expressed, including KLK2 which is known to regulate EGFR and pro-inflammatory cytokines and is regulated by MYC. Lastly, TRAF3IP2 encodes for a protein involved in regulating cytokines through members of the NFKB pathway. Conclusions: These findings suggest important biological insights into one of the rarer EOC histologies and may aid in the development of targeted treatment options. Research is on-going to incorporate additional features (e.g., DNA methylation, copy number) into a “systems biology” framework to better understand the molecular differences between EOC histologies. Citation Format: Brooke L. Fridley, Junqiang Dai, Rama Raghavan, Chen Wang, Pengcheng Lu, Stacey Winham, Madalene Earp, Kate Lawrenson, Simon A. Gayther, Kimberly R. Kalli, Ellen L. Goode. Transcriptome characterization of high grade serous and endometrioid epithelial ovarian cancer tumors. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 1833.

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.