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.

Epidemiological, biological, and molecular data suggest links between endometriosis and endometrial cancer, with recent epidemiological studies providing evidence for an association between a previous diagnosis of endometriosis and risk of endometrial cancer. We used genetic data as an alternative approach to investigate shared biological etiology of these two diseases. Genetic correlation analysis of summary level statistics from genomewide association studies (GWAS) using LD Score regression revealed moderate but significant genetic correlation (rg  = 0.23, P = 9.3 × 10-3 ), and SNP effect concordance analysis provided evidence for significant SNP pleiotropy (P = 6.0 × 10-3 ) and concordance in effect direction (P = 2.0 × 10-3 ) between the two diseases. Cross-disease GWAS meta-analysis highlighted 13 distinct loci associated at P ≤ 10-5 with both endometriosis and endometrial cancer, with one locus (SNP rs2475335) located within PTPRD associated at a genomewide significant level (P = 4.9 × 10-8 , OR = 1.11, 95% CI = 1.07-1.15). PTPRD acts in the STAT3 pathway, which has been implicated in both endometriosis and endometrial cancer. This study demonstrates the value of cross-disease genetic analysis to support epidemiological observations and to identify biological pathways of relevance to multiple diseases.

The CCNE1 locus, which encodes cyclin E1, is amplified in many types of cancer cells and is activated in hepatocellular carcinomas (HCCs) from patients infected with hepatitis B virus or adeno-associated virus type 2, due to integration of the virus nearby. We investigated cell-cycle and oncogenic effects of cyclin E1 overexpression in tissues of mice.

Reduced BRCA1 expression causes homologous recombination (HR) repair defects in high-grade serous ovarian cancers (HGSOCs). Here, we demonstrate that BRCA1 is transcriptionally activated by a previously unknown function of ZC3H18. We show that ZC3H18 is a DNA-binding protein that interacts with an E2F site in the BRCA1 promoter where it facilitates recruitment of E2F4 to an adjacent E2F site to promote BRCA1 transcription. Consistent with ZC3H18 role in activating BRCA1 expression, ZC3H18 depletion induces BRCA1 promoter methylation, reduces BRCA1 expression, disrupts HR, and sensitizes cells to DNA crosslinkers and poly(ADP-ribose) polymerase inhibitors. Moreover, in patient-derived xenografts and primary HGSOC tumors, ZC3H18 and E2F4 mRNA levels are positively correlated with BRCA1 mRNA levels, further supporting ZC3H18 role in regulating BRCA1. Given that ZC3H18 lies within 16q24.2, a region with frequent copy number loss in HGSOC, these findings suggest that ZC3H18 copy number losses could contribute to HR defects in HGSOC.

We analyzed whole exome sequencing data in germline DNA from 412 high grade serous ovarian cancer (HGSOC) cases from The Cancer Genome Atlas Project and identified 5,517 genes harboring a predicted deleterious germline coding mutation in at least one HGSOC case. Gene-set enrichment analysis showed enrichment for genes involved in DNA repair (p = 1.8×10-3). Twelve DNA repair genes - APEX1, APLF, ATX, EME1, FANCL, FANCM, MAD2L2, PARP2, PARP3, POLN, RAD54L and SMUG1 - were prioritized for targeted sequencing in up to 3,107 HGSOC cases, 1,491 cases of other epithelial ovarian cancer (EOC) subtypes and 3,368 unaffected controls of European origin. We estimated mutation prevalence for each gene and tested for associations with disease risk. Mutations were identified in both cases and controls in all genes except MAD2L2, where we found no evidence of mutations in controls. In FANCM we observed a higher mutation frequency in HGSOC cases compared to controls (29/3,107 cases, 0.96 percent; 13/3,368 controls, 0.38 percent; P=0.008) with little evidence for association with other subtypes (6/1,491, 0.40 percent; P=0.82). The relative risk of HGSOC associated with deleterious FANCM mutations was estimated to be 2.5 (95% CI 1.3 - 5.0; P=0.006). In summary, whole exome sequencing of EOC cases with large-scale replication in case-control studies has identified FANCM as a likely novel susceptibility gene for HGSOC, with mutations associated with a moderate increase in risk. These data may have clinical implications for risk prediction and prevention approaches for high-grade serous ovarian cancer in the future and a significant impact on reducing disease mortality.

We previously identified associations with ovarian cancer outcome at five genetic loci. To identify putatively causal genetic variants and target genes, we prioritized two ovarian outcome loci (1q22 and 19p12) for further study. Bioinformatic and functional genetic analyses indicated that MEF2D and ZNF100 are targets of candidate outcome variants at 1q22 and 19p12, respectively. At 19p12, the chromatin interaction of a putative regulatory element with the ZNF100 promoter region correlated with candidate outcome variants. At 1q22, putative regulatory elements enhanced MEF2D promoter activity and haplotypes containing candidate outcome variants modulated these effects. In a public dataset, MEF2D and ZNF100 expression were both associated with ovarian cancer progression-free or overall survival time. In an extended set of 6,162 epithelial ovarian cancer patients, we found that functional candidates at the 1q22 and 19p12 loci, as well as other regional variants, were nominally associated with patient outcome; however, no associations reached our threshold for statistical significance (p<1×10-5). Larger patient numbers will be needed to convincingly identify any true associations at these loci.

During the past decade, the age-adjusted mortality rate for endometrial cancer (EC) increased 1.9% annually with TP53 mutant (TP53-mu) EC disproportionally represented in advanced disease and deaths. Therefore, we aimed to assess pivotal molecular parameters that differentiate clinical outcomes in high- and low-risk EC. Using the Cancer Genome Atlas, we analyzed EC specimens with available DNA sequences and quantitative gene-specific RNA expression data. After polymerase ɛ (POLE)-mutant specimens were excluded, differential gene-specific mutations and mRNA expressions were annotated and integrated. Consequent to TP53-mu failure to induce p21, derepression of multiple oncogenes harboring promoter p21 repressive sites was observed, including CCNA2 and FOXM1 (P < .001 compared with TP53 wild type [TP53-wt]). TP53-wt EC with high CCNA2 expression (CCNA2-H) had a targeted transcriptomic profile similar to that of TP53-mu EC, suggesting CCNA2 is a seminal determinant for both TP53-wt and TP53-mu EC. CCNA2 enhances E2F1 function, upregulating FOXM1 and CIP2A, as observed in TP53-mu and CCNA2-H TP53-wt EC (P < .001). CIP2A inhibits protein phosphatase 2A, leading to AKT inactivation of GSK3β and restricted oncoprotein degradation; PPP2R1A and FBXW7 mutations yield similar results. Upregulation of FOXM1 and failed degradation of FOXM1 is evidenced by marked upregulation of multiple homologous recombination genes (P < .001). Integrating these molecular aberrations generated a molecular biomarker panel with significant prognostic discrimination (P = 5.8×10-7); adjusting for age, histology, grade, myometrial invasion, TP53 status, and stage, only CCNA2-H/E2F1-H (P = .0003), FBXW7-mu/PPP2R1A-mu (P = .0002), and stage (P = .017) were significant. The generated prognostic molecular classification system identifies dissimilar signaling aberrations potentially amenable to targetable therapeutic options.

High stage and recurrent ovarian clear cell carcinoma (OCC) are associated with poor prognosis and resistance to chemotherapy. A distinguishing histological feature of OCC is abundant cytoplasmic stores of glucose, in the form of glycogen, that can be mobilized for cellular metabolism. Here, we report the effect on preclinical models of OCC of disrupting glycogen utilization using the glucose analogue 2-deoxy-D-glucose (2DG). At concentrations significantly lower than previously reported for other cancers, 2DG markedly improves the efficacy in vitro of carboplatin chemotherapy against chemo-sensitive TOV21G and chemo-resistant OVTOKO OCC cell lines, and this is accompanied by the depletion of glycogen. Of note, 2DG doses-of more than 10-fold lower than previously reported for other cancers-significantly improve the efficacy of carboplatin against cell line and patient-derived xenograft models in mice that mimic the chemo-responsiveness of OCC. These findings are encouraging, in that 2DG doses, which are substantially lower than previously reported to cause adverse events in cancer patients, can safely and significantly improve the efficacy of carboplatin against OCC. Our results thus justify clinical trials to evaluate whether low dose 2DG improves the efficacy of carboplatin in OCC patients.

Somatic loss of the tumour suppressor RB1 is a common event in tubo-ovarian high-grade serous carcinoma (HGSC), which frequently co-occurs with alterations in homologous recombination DNA repair genes including BRCA1 and BRCA2 (BRCA). We examined whether tumour expression of RB1 was associated with survival across ovarian cancer histotypes (HGSC, endometrioid (ENOC), clear cell (CCOC), mucinous (MOC), low-grade serous carcinoma (LGSC)), and how co-occurrence of germline BRCA pathogenic variants and RB1 loss influences long-term survival in a large series of HGSC.

Poly(ADP-ribose) polymerase (PARP) inhibitors have yielded encouraging responses in high-grade serous ovarian carcinomas (HGSOCs), but the optimal treatment setting remains unknown. We assessed the effect of niraparib on HGSOC patient-derived xenograft (PDX) models as well as the relationship between certain markers of homologous recombination (HR) status, including BRCA1/2 mutations and formation of RAD51 foci after DNA damage, and response of these PDXs to niraparib in vivo.

The mechanisms of recurrence have been under-studied in rare histologies of invasive epithelial ovarian cancer (EOC) (endometrioid, clear cell, mucinous, and low-grade serous). We hypothesised the existence of an expression signature predictive of outcome in the rarer histologies.

RNA editing in mammals is a form of post-transcriptional modification in which adenosine is converted to inosine by the adenosine deaminases acting on RNA (ADAR) family of enzymes. Based on evidence of altered ADAR expression in epithelial ovarian cancers (EOC), we hypothesized that single nucleotide polymorphisms (SNPs) in ADAR genes modify EOC susceptibility, potentially by altering ovarian tissue gene expression. Using directly genotyped and imputed data from 10,891 invasive EOC cases and 21,693 controls, we evaluated the associations of 5,303 SNPs in ADAD1, ADAR, ADAR2, ADAR3, and SND1. Unconditional logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI), with adjustment for European ancestry. We conducted gene-level analyses using the Admixture Maximum Likelihood (AML) test and the Sequence-Kernel Association test for common and rare variants (SKAT-CR). Association analysis revealed top risk-associated SNP rs77027562 (OR (95% CI)= 1.39 (1.17-1.64), P=1.0x10-4) in ADAR3 and rs185455523 in SND1 (OR (95% CI)= 0.68 (0.56-0.83), P=2.0x10-4). When restricting to serous histology (n=6,500), the magnitude of association strengthened for rs185455523 (OR=0.60, P=1.0x10-4). Gene-level analyses revealed that variation in ADAR was associated (P<0.05) with EOC susceptibility, with PAML=0.022 and PSKAT-CR=0.020. Expression quantitative trait locus analysis in EOC tissue revealed significant associations (P<0.05) with ADAR expression for several SNPs in ADAR, including rs1127313 (G/A), a SNP in the 3' untranslated region. In summary, germline variation involving RNA editing genes may influence EOC susceptibility, warranting further investigation of inherited and acquired alterations affecting RNA editing.

Epithelial ovarian cancer (EOC) is the fifth leading cause of cancer death among women in the United States (5 % of cancer deaths). The standard treatment for patients with advanced EOC is initial debulking surgery followed by carboplatin-paclitaxel combination chemotherapy. Unfortunately, with chemotherapy most patients relapse and die resulting in a five-year overall survival around 45 %. Thus, finding novel therapeutics for treating EOC is essential. Connectivity Mapping (CMAP) has been used widely in cancer drug discovery and generally has relied on cancer cell line gene expression and drug phenotype data. Therefore, we took a CMAP approach based on tumor information and clinical endpoints from high grade serous EOC patients.

Ovarian carcinosarcoma is a rare subtype of ovarian cancer with poor clinical outcomes. The low incidence of this disease makes accrual to large clinical trials challenging. However, studies have shown that treatment responses in patient-derived xenograft (PDX) models correlate with matched-patient responses in the clinic, supporting their use for preclinical testing of standard and novel therapies. An ovarian carcinosarcoma PDX is presented herein and showed resistance to carboplatin and paclitaxel (similar to the patient) but exhibited significant sensitivity to ifosfamide and paclitaxel. The PDX demonstrated overexpression of EGFR mRNA and gene amplification by array comparative genomic hybridization (log2 ratio 0.399). EGFR phosphorylation was also detected. Angiogensis and insulin-like growth factor pathways were also implicated by overexpression of VEGFC and IRS1. In order to improve response to chemotherapy, the PDX was treated with carboplatin/paclitaxel with or without a pan-HER and VEGF inhibitor (BMS-690514) but there was no tumor growth inhibition or improved animal survival, which may be explained by a KRAS mutation. Resistance was also observed when the IGF-1R inhibitor BMS-754807 was combined with carboplatin/paclitaxel. Because poly (ADP-ribose) polymerase inhibitors have activity in ovarian cancer patients, with and without BRCA mutations, ABT-888 was also tested but found to have no activity. Pathogenic mutations were also detected in TP53 and PIK3CA. In conclusion, ifosfamide/paclitaxel was superior to carboplatin/paclitaxel in this ovarian carcinosarcoma PDX and gene overexpression or amplification alone was not sufficient to predict response to targeted therapy. Better predictive markers of response are needed.

Genome-wide association studies (GWAS) have identified 12 epithelial ovarian cancer (EOC) susceptibility alleles. The pattern of association at these loci is consistent in BRCA1 and BRCA2 mutation carriers who are at high risk of EOC. After imputation to 1000 Genomes Project data, we assessed associations of 11 million genetic variants with EOC risk from 15,437 cases unselected for family history and 30,845 controls and from 15,252 BRCA1 mutation carriers and 8,211 BRCA2 mutation carriers (3,096 with ovarian cancer), and we combined the results in a meta-analysis. This new study design yielded increased statistical power, leading to the discovery of six new EOC susceptibility loci. Variants at 1p36 (nearest gene, WNT4), 4q26 (SYNPO2), 9q34.2 (ABO) and 17q11.2 (ATAD5) were associated with EOC risk, and at 1p34.3 (RSPO1) and 6p22.1 (GPX6) variants were specifically associated with the serous EOC subtype, all with P < 5 × 10(-8). Incorporating these variants into risk assessment tools will improve clinical risk predictions for BRCA1 and BRCA2 mutation carriers.

Mammographic density reflects the amount of stromal and epithelial tissues in relation to adipose tissue in the breast and is a strong risk factor for breast cancer. Here we report the results from meta-analysis of genome-wide association studies (GWAS) of three mammographic density phenotypes: dense area, non-dense area and percent density in up to 7,916 women in stage 1 and an additional 10,379 women in stage 2. We identify genome-wide significant (P<5 × 10(-8)) loci for dense area (AREG, ESR1, ZNF365, LSP1/TNNT3, IGF1, TMEM184B and SGSM3/MKL1), non-dense area (8p11.23) and percent density (PRDM6, 8p11.23 and TMEM184B). Four of these regions are known breast cancer susceptibility loci, and four additional regions were found to be associated with breast cancer (P<0.05) in a large meta-analysis. These results provide further evidence of a shared genetic basis between mammographic density and breast cancer and illustrate the power of studying intermediate quantitative phenotypes to identify putative disease-susceptibility loci.

ABCB1 encodes the multi-drug efflux pump P-glycoprotein (P-gp) and has been implicated in multi-drug resistance. We comprehensively evaluated this gene and flanking regions for an association with clinical outcome in epithelial ovarian cancer (EOC).

ATIC, SHMT2, and SLC46A1 have essential roles in one-carbon (1-C) transfer. The authors examined whether associations between ovarian carcinoma and 15 variants in these genes are modified by regular multivitamin use, a source of 1-C donors, among Caucasian participants from two US case-control studies. Using a phased study design, variant-by-multivitamin interactions were tested, and associations between variants and ovarian carcinoma were reported stratified by multivitamin supplement use. Per-allele risk associations were modified by multivitamin use at six variants among 655 cases and 920 controls (Phase 1). In a larger sample of 968 cases and 1,265 controls (Phases 1 and 2), interactions were significant (P ≤ 0.03) for two variants, particularly among regular multivitamin users: ATIC rs7586969 [odds ratio (OR) = 0.7, 95% confidence interval (CI) = 0.6-0.9] and ATIC rs16853834 (OR = 1.5, 95% CI = 1.1-2.0). The two ATIC single nucleotide polymorphisms (SNPs) did not share the same haplotype; however, the haplotypes they comprised mirrored their SNP risk associations among regular multivitamin supplement users. A multi-variant analysis was also performed by comparing the observed likelihood ratio test statistic from adjusted models with and without the two ATIC variant-by-multivitamin interaction terms with a null distribution of test statistics generated by permuting case status 10,000 times. The corresponding observed P value of 0.001 was more extreme than the permutation-derived P value of 0.009, suggesting rejection of the null hypothesis of no association. In summary, there is little statistical evidence that the 15 variants are independently associated with risk of ovarian carcinoma. However, the statistical interaction of ATIC variants with regular multivitamin intake, when evaluated at both the SNP and gene level, may support these findings as relevant to ovarian health and disease processes.

Host antiviral genes are important regulators of antiviral immunity and plausible genetic determinants of immune response heterogeneity after vaccination. We genotyped and analyzed 307 common candidate tagSNPs from 12 antiviral genes in a cohort of 745 schoolchildren immunized with two doses of measles-mumps-rubella (MMR) vaccine. Associations between SNPs/haplotypes and measles virus-specific immune outcomes were assessed using linear regression methodologies in Caucasians and African-Americans. Genetic variants within the DDX58/RIG-I gene, including a coding polymorphism (rs3205166/Val800Val), were associated as single-SNPs (p≤0.017; although these SNPs did not remain significant after correction for false discovery rate/FDR) and in haplotype-level analysis, with measles-specific antibody variations in Caucasians (haplotype allele p-value=0.021; haplotype global p-value=0.076). Four DDX58 polymorphisms, in high LD, demonstrated also associations (after correction for FDR) with variations in both measles-specific IFN-γ and IL-2 secretion in Caucasians (p≤0.001, q=0.193). Two intronic OAS1 polymorphisms, including the functional OAS1 SNP rs10774671 (p=0.003), demonstrated evidence of association with a significant allele-dose-related increase in neutralizing antibody levels in African-Americans. Genotype and haplotype-level associations demonstrated the role of ADAR genetic variants, including a non-synonymous SNP (rs2229857/Arg384Lys; p=0.01), in regulating measles virus-specific IFN-γ Elispot responses in Caucasians (haplotype global p-value=0.017). After correction for FDR, 15 single-SNP associations (11 SNPs in Caucasians and 4 SNPs in African-Americans) still remained significant at the q-value<0.20. In conclusion, our findings strongly point to genetic variants/genes, involved in antiviral sensing and antiviral control, as critical determinants, differentially modulating the adaptive immune responses to live attenuated measles vaccine in Caucasians and African-Americans.

The role of the innate immune response in colorectal cancer is understudied. We examined the survival of colorectal cancer patients in relation to eosinophils, innate immune cells, infiltrating the tumor. Tissue microarrays were constructed from paraffin-embedded tumor tissues collected between 1986 and 2002 from 441 post-menopausal women diagnosed with colorectal cancer in the Iowa Women's Health Study. Tissue microarrays were stained with an eosinophil peroxidase antibody. Eosinophils in epithelial and stromal tissues within the tumor (called epithelial and stromal eosinophils, hereafter) were counted and scored into three and four categories, respectively. In addition, the degree of eosinophil degranulation (across epithelial and stromal tissues combined) was quantified and similarly categorized. We used Cox regression to estimate the hazard ratios and 95% confidence interval for all-cause and colorectal cancer death during 5-year follow-up after diagnosis and during follow-up through 2011 ('total follow-up'). The hazard ratios associated with eosinophil scores were adjusted for age of diagnosis, SEER (Surveillance, Epidemiology, and End Results) stage, tumor grade, body mass, and smoking history. High tumor stromal eosinophil score was inversely correlated with age and stage, and was associated with a decreased risk for all-cause and colorectal cancer death: hazard ratios (95% confidence intervals) were 0.61 (0.36-1.02; P-trend=0.02) and 0.48 (0.24-0.93; P-trend=0.01), respectively, during the 5-year follow-up for the highest vs lowest category. The inverse associations also existed for total follow-up for all-cause and colorectal cancer death for the highest vs lowest stromal eosinophil score: hazard ratios (95% confidence intervals) were 0.72 (0.48-1.08; P-trend=0.04) and 0.61 (0.34-1.12; P-trend=0.04), respectively. Further adjustment for treatment, comorbidities, additional lifestyle factors, tumor location, or molecular markers did not markedly change the associations, while adjustment for cytotoxic T cells slightly attenuated all associations. The infiltration of tumors with eosinophils, especially in stromal tissue, may be an important prognostic factor in colorectal cancer.