Kalli, Kimberly R.

Abstract Background: Endometrioid carcinoma (EC) and clear cell carcinoma (CC) histotypes of epithelial ovarian cancer are understudied compared with the more common high-grade serous carcinomas (HGSC). We therefore sought to characterize EC and CC transcriptomes in relation to HGSC. Methods: Following bioinformatics processing and gene abundance normalization, differential expression analysis of RNA sequence data collected on fresh-frozen tumors was completed with nonparametric statistical analysis methods (55 ECs, 19 CCs, 112 HGSCs). Association of gene expression with progression-free survival (PFS) was completed with Cox proportional hazards models. Eight additional multi-histotype expression array datasets (N = 852 patients) were used for replication. Results: In the discovery set, tumors generally clustered together by histotype. Thirty-two protein-coding genes were differentially expressed across histotype (P < 1 × 10−10) and showed similar associations in replication datasets, including MAP2K6, KIAA1324, CDH1, ENTPD5, LAMB1, and DRAM1. Nine genes associated with PFS (P < 0.0001) showed similar associations in replication datasets. In particular, we observed shorter PFS time for CC and EC patients with high gene expression for CCNB2, CORO2A, CSNK1G1, FRMD8, LIN54, LINC00664, PDK1, and PEX6, whereas, the converse was observed for HGSC patients. Conclusions: The results suggest important histotype differences that may aid in the development of treatment options, particularly those for patients with EC or CC. Impact: We present replicated findings on transcriptomic differences and how they relate to clinical outcome for two of the rarer ovarian cancer histotypes of EC and CC, along with comparison with the common histotype of HGSC. Cancer Epidemiol Biomarkers Prev; 27(9); 1101–9. ©2018 AACR.

Abstract Purpose: Folate receptor alpha (FR) is overexpressed in several cancers. Endogenous immunity to the FR has been demonstrated in patients and suggests the feasibility of targeting FR with vaccine or other immune therapies. CD4 helper T cells are central to the development of coordinated immunity, and prior work shows their importance in protecting against relapse. Our previous identification of degenerate HLA-class II epitopes from human FR led to the development of a broad coverage epitope pool potentially useful in augmenting antigen-specific immune responses in most patients. Patients and Methods: We conducted a phase I clinical trial testing safety and immunogenicity of this vaccine, enrolling patients with ovarian cancer or breast cancer who completed conventional treatment and who showed no evidence of disease. Patients were initially treated with low-dose cyclophosphamide and then vaccinated 6 times, monthly. Immunity and safety were examined during the vaccine period and up to 1 year later. Results: Vaccination was well tolerated in all patients. Vaccine elicited or augmented immunity in more than 90% of patients examined. Unlike recall immunity to tetanus toxoid (TT), FR T-cell responses developed slowly over the course of vaccination with a median time to maximal immunity in 5 months. Despite slow development of immunity, responsiveness appeared to persist for at least 12 months. Conclusions: The results demonstrate that it is safe to augment immunity to the FR tumor antigen, and the developed vaccine is testable for therapeutic activity in most patients whose tumors express FR, regardless of HLA genotype. Clin Cancer Res; 24(13); 3014–25. ©2018 AACR.

High-grade serous ovarian cancer (HGSOC) is a complex disease in which initiation and progression have been associated with copy number alterations, epigenetic processes, and, to a lesser extent, germline variation. We hypothesized that, when summarized at the gene level, tumor methylation and germline genetic variation, alone or in combination, influence tumor gene expression in HGSOC. We used Elastic Net (ENET) penalized regression method to evaluate these associations and adjust for somatic copy number in 3 independent data sets comprising tumors from more than 470 patients. Penalized regression models of germline variation, with or without methylation, did not reveal a role in HGSOC gene expression. However, we observed significant association between regional methylation and expression of 5 genes ( WDPCP, KRT6C, BRCA2, EFCAB13, and ZNF283). CpGs retained in ENET model for BRCA2 and ZNF283 appeared enriched in several regulatory elements, suggesting that regularized regression may provide a novel utility for integrative genomic analysis.

ABSTRACTX‐chromosome inactivation (XCI) epigenetically silences transcription of an X chromosome in females; patterns of XCI are thought to be aberrant in women's cancers, but are understudied due to statistical challenges. We develop a two‐stage statistical framework to assess skewed XCI and evaluate gene‐level patterns of XCI for an individual sample by integration of RNA sequence, copy number alteration, and genotype data. Our method relies on allele‐specific expression (ASE) to directly measure XCI and does not rely on male samples or paired normal tissue for comparison. We model ASE using a two‐component mixture of beta distributions, allowing estimation for a given sample of the degree of skewness (based on a composite likelihood ratio test) and the posterior probability that a given gene escapes XCI (using a Bayesian beta‐binomial mixture model). To illustrate the utility of our approach, we applied these methods to data from tumors of ovarian cancer patients. Among 99 patients, 45 tumors were informative for analysis and showed evidence of XCI skewed toward a particular parental chromosome. For 397 X‐linked genes, we observed tumor XCI patterns largely consistent with previously identified consensus states based on multiple normal tissue types. However, 37 genes differed in XCI state between ovarian tumors and the consensus state; 17 genes aberrantly escaped XCI in ovarian tumors (including many oncogenes), whereas 20 genes were unexpectedly inactivated in ovarian tumors (including many tumor suppressor genes). These results provide evidence of the importance of XCI in ovarian cancer and demonstrate the utility of our two‐stage analysis.

Abstract Ligation of programmed cell death-1 (PD-1) in the tumor microenvironment is known to inhibit effective adaptive antitumor immunity. Blockade of PD-1 in humans has resulted in impressive, durable regression responses in select tumor types. However, durable responses have been elusive in ovarian cancer patients. PD-1 was recently shown to be expressed on and thereby impair the functions of tumor-infiltrating murine and human myeloid dendritic cells (TIDC) in ovarian cancer. In the present work, we characterize the regulation of PD-1 expression and the effects of PD-1 blockade on TIDC. Treatment of TIDC and bone marrow–derived dendritic cells (DC) with IL10 led to increased PD-1 expression. Both groups of DCs also responded to PD-1 blockade by increasing production of IL10. Similarly, treatment of ovarian tumor–bearing mice with PD-1 blocking antibody resulted in an increase in IL10 levels in both serum and ascites. While PD-1 blockade or IL10 neutralization as monotherapies were inefficient, combination of these two led to improved survival and delayed tumor growth; this was accompanied by augmented antitumor T- and B-cell responses and decreased infiltration of immunosuppressive MDSC. Taken together, our findings implicate compensatory release of IL10 as one of the adaptive resistance mechanisms that undermine the efficacy of anti–PD-1 (or anti–PD-L1) monotherapies and prompt further studies aimed at identifying such resistance mechanisms. Cancer Res; 77(23); 6667–78. ©2017 AACR.

Abstract High grade serous ovarian cancer (HGSOC) is a complex disease in which initiation and progression have been associated with gene mutation, DNA methylation changes, genetic variation, and epigenetic processes. Variation in several susceptibility regions and cancer-typical global methylation patterns have been observed in HGSOC; however, this knowledge has not been compelling in understanding HGSOC intiation or progression. As ingetration of genomic, epigenomic, and transcriptomic data has increased mechanistic understanding in other cancers, we hypothesized that tumor methylation alone or in combination with germline genetic variation influences tumor gene expression in HGSOC. We examined three nested models using an Elastnic Net (ENET) penalized regression method while adjusting for somatic copy number (CNV): a) germline genotype and tumor DNA methylation (full model), b) genotype only, and c) DNA methylation only. We included 339 cases from The Cancer Genome Atlas (TCGA), 54 cases from Mayo Clinic, and 78 cases from the Australian Ovarian Cancer Study (AOCS). Genotyping and copy number calls on germline DNA, expression, methylation and copy number on somatic samples were collected and analyzed on different platforms separately at each study site. We excluded genes with low overall expression and thus analyzed a total of 11,922 genes available in three datasets ( Ensembl IDs, 500kb window up- and downstream). In general, combining genomic data in HGSOC did not reveal a role for germline genetic variation in altering gene expression. However, in methylation only models 79 genes were associated with differential expression in the TCGA cases (permutation multiple testing adjusted p-val <0.05), in the Mayo cases (unadjusted p-val <0.05) and AOCS cases (unadjusted p-val <0.05). A known tummor suppressor (FBXW7) was associated with differential expression in the three datasets at p-val <0.01. This work demonstrates the feasibility, utility, and statistical power of ENET gene-level analyses incoporating maximal genomic information. Citation Format: Yanina Natanzon, Madalene Earp, Julie M. Cunningham, Kimberly R. Kalli, Stacey J. Winham, Sebastian M. Armasu, Melissa C. Larson, Chen Wang, David D. Bowtell, Dale W. Garsed, Ellen Goode. Omics data integration analysis in high grade serous ovarian cancer: results from three studies [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 3368. doi:10.1158/1538-7445.AM2017-3368

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