Disruption in circadian gene expression, whether due to genetic variation or environmental factors (e.g., light at night, shiftwork), is associated with increased incidence of breast, prostate, gastrointestinal and hematologic cancers and gliomas. Circadian genes are highly expressed in the ovaries where they regulate ovulation; circadian disruption is associated with several ovarian cancer risk factors (e.g., endometriosis). However, no studies have examined variation in germline circadian genes as predictors of ovarian cancer risk and invasiveness. The goal of the current study was to examine single nucleotide polymorphisms (SNPs) in circadian genes BMAL1, CRY2, CSNK1E, NPAS2, PER3, REV1 and TIMELESS and downstream transcription factors KLF10 and SENP3 as predictors of risk of epithelial ovarian cancer (EOC) and histopathologic subtypes. The study included a test set of 3,761 EOC cases and 2,722 controls and a validation set of 44,308 samples including 18,174 (10,316 serous) cases and 26,134 controls from 43 studies participating in the Ovarian Cancer Association Consortium (OCAC). Analysis of genotype data from 36 genotyped SNPs and 4600 imputed SNPs indicated that the most significant association was rs117104877 in BMAL1 (OR = 0.79, 95% CI = 0.68-0.90, p = 5.59 × 10-4]. Functional analysis revealed a significant down regulation of BMAL1 expression following cMYC overexpression and increasing transformation in ovarian surface epithelial (OSE) cells as well as alternative splicing of BMAL1 exons in ovarian and granulosa cells. These results suggest that variation in circadian genes, and specifically BMAL1, may be associated with risk of ovarian cancer, likely through disruption of hormonal pathways.

A locus at 19p13 is associated with breast cancer (BC) and ovarian cancer (OC) risk. Here we analyse 438 SNPs in this region in 46,451 BC and 15,438 OC cases, 15,252 BRCA1 mutation carriers and 73,444 controls and identify 13 candidate causal SNPs associated with serous OC (P=9.2 × 10(-20)), ER-negative BC (P=1.1 × 10(-13)), BRCA1-associated BC (P=7.7 × 10(-16)) and triple negative BC (P-diff=2 × 10(-5)). Genotype-gene expression associations are identified for candidate target genes ANKLE1 (P=2 × 10(-3)) and ABHD8 (P<2 × 10(-3)). Chromosome conformation capture identifies interactions between four candidate SNPs and ABHD8, and luciferase assays indicate six risk alleles increased transactivation of the ADHD8 promoter. Targeted deletion of a region containing risk SNP rs56069439 in a putative enhancer induces ANKLE1 downregulation; and mRNA stability assays indicate functional effects for an ANKLE1 3'-UTR SNP. Altogether, these data suggest that multiple SNPs at 19p13 regulate ABHD8 and perhaps ANKLE1 expression, and indicate common mechanisms underlying breast and ovarian cancer risk.

Common variants in 94 loci have been associated with breast cancer including 15 loci with genome-wide significant associations (P<5 × 10(-8)) with oestrogen receptor (ER)-negative breast cancer and BRCA1-associated breast cancer risk. In this study, to identify new ER-negative susceptibility loci, we performed a meta-analysis of 11 genome-wide association studies (GWAS) consisting of 4,939 ER-negative cases and 14,352 controls, combined with 7,333 ER-negative cases and 42,468 controls and 15,252 BRCA1 mutation carriers genotyped on the iCOGS array. We identify four previously unidentified loci including two loci at 13q22 near KLF5, a 2p23.2 locus near WDR43 and a 2q33 locus near PPIL3 that display genome-wide significant associations with ER-negative breast cancer. In addition, 19 known breast cancer risk loci have genome-wide significant associations and 40 had moderate associations (P<0.05) with ER-negative disease. Using functional and eQTL studies we implicate TRMT61B and WDR43 at 2p23.2 and PPIL3 at 2q33 in ER-negative breast cancer aetiology. All ER-negative loci combined account for ∼11% of familial relative risk for ER-negative disease and may contribute to improved ER-negative and BRCA1 breast cancer risk prediction.

Clinical genetic testing is commercially available for rs61764370, an inherited variant residing in a KRAS 3' UTR microRNA binding site, based on suggested associations with increased ovarian and breast cancer risk as well as with survival time. However, prior studies, emphasizing particular subgroups, were relatively small. Therefore, we comprehensively evaluated ovarian and breast cancer risks as well as clinical outcome associated with rs61764370.

Several common alleles have been shown to be associated with breast and/or ovarian cancer risk for BRCA1 and BRCA2 mutation carriers. Recent genome-wide association studies of breast cancer have identified eight additional breast cancer susceptibility loci: rs1011970 (9p21, CDKN2A/B), rs10995190 (ZNF365), rs704010 (ZMIZ1), rs2380205 (10p15), rs614367 (11q13), rs1292011 (12q24), rs10771399 (12p11 near PTHLH) and rs865686 (9q31.2).

Ovarian cancer represents the most lethal tumor type among malignancies of the female reproductive system. Overall survival rates remain low. In this study, we identify the serine protease inhibitor Kazal type 1 (SPINK1) as a potential therapeutic target for a subset of ovarian cancers. We show that SPINK1 drives ovarian cancer cell proliferation through activation of epidermal growth factor receptor (EGFR) signaling, and that SPINK1 promotes resistance to anoikis through a distinct mechanism involving protease inhibition. In analyses of ovarian tumor specimens from a Mayo Clinic cohort of 490 patients, we further find that SPINK1 immunostaining represents an independent prognostic factor for poor survival, with the strongest association in patients with nonserous histological tumor subtypes (endometrioid, clear cell, and mucinous). This study provides novel insight into the fundamental processes underlying ovarian cancer progression, and also suggests new avenues for development of molecularly targeted therapies.

Genome-wide association studies have reported 11 regions conferring risk of high-grade serous epithelial ovarian cancer (HGSOC). Expression quantitative trait locus (eQTL) analyses can identify candidate susceptibility genes at risk loci. Here we evaluate cis-eQTL associations at 47 regions associated with HGSOC risk (P≤10(-5)). For three cis-eQTL associations (P<1.4 × 10(-3), FDR<0.05) at 1p36 (CDC42), 1p34 (CDCA8) and 2q31 (HOXD9), we evaluate the functional role of each candidate by perturbing expression of each gene in HGSOC precursor cells. Overexpression of HOXD9 increases anchorage-independent growth, shortens population-doubling time and reduces contact inhibition. Chromosome conformation capture identifies an interaction between rs2857532 and the HOXD9 promoter, suggesting this SNP is a leading causal variant. Transcriptomic profiling after HOXD9 overexpression reveals enrichment of HGSOC risk variants within HOXD9 target genes (P=6 × 10(-10) for risk variants (P<10(-4)) within 10 kb of a HOXD9 target gene in ovarian cells), suggesting a broader role for this network in genetic susceptibility to HGSOC.

Epithelial ovarian cancer (EOC) is the fifth leading cause of cancer mortality in American women. Normal ovarian physiology is intricately connected to small GTP binding proteins of the Ras superfamily (Ras, Rho, Rab, Arf, and Ran) which govern processes such as signal transduction, cell proliferation, cell motility, and vesicle transport. We hypothesized that common germline variation in genes encoding small GTPases is associated with EOC risk. We investigated 322 variants in 88 small GTPase genes in germline DNA of 18,736 EOC patients and 26,138 controls of European ancestry using a custom genotype array and logistic regression fitting log-additive models. Functional annotation was used to identify biofeatures and expression quantitative trait loci that intersect with risk variants. One variant, ARHGEF10L (Rho guanine nucleotide exchange factor 10 like) rs2256787, was associated with increased endometrioid EOC risk (OR = 1.33, p = 4.46 x 10-6). Other variants of interest included another in ARHGEF10L, rs10788679, which was associated with invasive serous EOC risk (OR = 1.07, p = 0.00026) and two variants in AKAP6 (A-kinase anchoring protein 6) which were associated with risk of invasive EOC (rs1955513, OR = 0.90, p = 0.00033; rs927062, OR = 0.94, p = 0.00059). Functional annotation revealed that the two ARHGEF10L variants were located in super-enhancer regions and that AKAP6 rs927062 was associated with expression of GTPase gene ARHGAP5 (Rho GTPase activating protein 5). Inherited variants in ARHGEF10L and AKAP6, with potential transcriptional regulatory function and association with EOC risk, warrant investigation in independent EOC study populations.

Endometriosis is characterized by the presence of functional endometrial tissue outside of the uterine cavity. It affects 1 in 10 women of reproductive age. This chronic condition commonly leads to consequences such as pelvic pain, dysmenorrhea, infertility and an elevated risk of epithelial ovarian cancer. Despite the prevalence of endometriosis and its impact on women's lives, there are relatively few in vitro and in vivo models available for studying the complex disease biology, pathophysiology, and for use in the preclinical development of novel therapies. The goal of this study was to develop a novel three-dimensional (3D) cell culture model of ovarian endometriosis and to test whether it is more reflective of endometriosis biology than traditional two dimensional (2D) monolayer cultures.

Mammalian target of rapamycin (mTOR) represents a key downstream intermediate for a myriad of oncogenic receptor tyrosine kinases. In the case of the insulin-like growth factor (IGF) pathway, the mTOR complex (mTORC1) mediates IGF-1 receptor (IGF-1R)-induced estrogen receptor alpha (ERα) phosphorylation/activation and leads to increased proliferation and growth in breast cancer cells. As a result, the prevalence of mTOR inhibitors combined with hormonal therapy has increased in recent years. Conversely, activated mTORC1 provides negative feedback regulation of IGF signaling via insulin receptor substrate (IRS)-1/2 serine phosphorylation and subsequent proteasomal degradation. Thus, the IGF pathway may provide escape (e.g. de novo or acquired resistance) from mTORC1 inhibitors. It is therefore plausible that combined inhibition of mTORC1 and IGF-1R for select subsets of ER-positive breast cancer patients presents as a viable therapeutic option.

Sclerosing adenosis (SA), found in ¼ of benign breast disease (BBD) biopsies, is a histological feature characterized by lobulocentric proliferation of acini and stromal fibrosis and confers a two-fold increase in breast cancer risk compared to women in the general population. We evaluated a NanoString-based gene expression assay to model breast cancer risk using RNA derived from formalin-fixed, paraffin-embedded (FFPE) biopsies with SA.

Interest in preclinical drug development for ovarian cancer has stimulated development of patient-derived xenograft (PDX) or tumorgraft models. However, the unintended formation of human lymphoma in severe combined immunodeficiency (SCID) mice from Epstein-Barr virus (EBV)-infected human lymphocytes can be problematic. In this study, we have characterized ovarian cancer PDXs which developed human lymphomas and explore methods to suppress lymphoproliferative growth. Fresh human ovarian tumors from 568 patients were transplanted intraperitoneally in SCID mice. A subset of PDX models demonstrated atypical patterns of dissemination with mediastinal masses, hepatosplenomegaly, and CD45-positive lymphoblastic atypia without ovarian tumor engraftment. Expression of human CD20 but not CD3 supported a B-cell lineage, and EBV genomes were detected in all lymphoproliferative tumors. Immunophenotyping confirmed monoclonal gene rearrangements consistent with B-cell lymphoma, and global gene expression patterns correlated well with other human lymphomas. The ability of rituximab, an anti-CD20 antibody, to suppress human lymphoproliferation from a patient's ovarian tumor in SCID mice and prevent growth of an established lymphoma led to a practice change with a goal to reduce the incidence of lymphomas. A single dose of rituximab during the primary tumor heterotransplantation process reduced the incidence of CD45-positive cells in subsequent PDX lines from 86.3% (n = 117 without rituximab) to 5.6% (n = 160 with rituximab), and the lymphoma rate declined from 11.1% to 1.88%. Taken together, investigators utilizing PDX models for research should routinely monitor for lymphoproliferative tumors and consider implementing methods to suppress their growth.

KLK15 over-expression is reported to be a significant predictor of reduced progression-free survival and overall survival in ovarian cancer. Our aim was to analyse the KLK15 gene for putative functional single nucleotide polymorphisms (SNPs) and assess the association of these and KLK15 HapMap tag SNPs with ovarian cancer survival.

Breast cancer is a common disease partially caused by genetic risk factors. Germline pathogenic variants in DNA repair genes BRCA1, BRCA2, PALB2, ATM, and CHEK2 are associated with breast cancer risk. FANCM, which encodes for a DNA translocase, has been proposed as a breast cancer predisposition gene, with greater effects for the ER-negative and triple-negative breast cancer (TNBC) subtypes. We tested the three recurrent protein-truncating variants FANCM:p.Arg658*, p.Gln1701*, and p.Arg1931* for association with breast cancer risk in 67,112 cases, 53,766 controls, and 26,662 carriers of pathogenic variants of BRCA1 or BRCA2. These three variants were also studied functionally by measuring survival and chromosome fragility in FANCM -/- patient-derived immortalized fibroblasts treated with diepoxybutane or olaparib. We observed that FANCM:p.Arg658* was associated with increased risk of ER-negative disease and TNBC (OR = 2.44, P = 0.034 and OR = 3.79; P = 0.009, respectively). In a country-restricted analysis, we confirmed the associations detected for FANCM:p.Arg658* and found that also FANCM:p.Arg1931* was associated with ER-negative breast cancer risk (OR = 1.96; P = 0.006). The functional results indicated that all three variants were deleterious affecting cell survival and chromosome stability with FANCM:p.Arg658* causing more severe phenotypes. In conclusion, we confirmed that the two rare FANCM deleterious variants p.Arg658* and p.Arg1931* are risk factors for ER-negative and TNBC subtypes. Overall our data suggest that the effect of truncating variants on breast cancer risk may depend on their position in the gene. Cell sensitivity to olaparib exposure, identifies a possible therapeutic option to treat FANCM-associated tumors.

Critical developmental “master transcription factors” (MTFs) can be subverted during tumorigenesis to control oncogenic transcriptional programs. Current approaches to identifying MTFs rely on ChIP-seq data, which is unavailable for many cancers. We developed the CaCTS (Cancer Core Transcription factor Specificity) algorithm to prioritize candidate MTFs using pan-cancer RNA sequencing data. CaCTS identified candidate MTFs across 34 tumor types and 140 subtypes including predictions for cancer types/subtypes for which MTFs are unknown, including e.g. PAX8, SOX17, and MECOM as candidates in ovarian cancer (OvCa). In OvCa cells, consistent with known MTF properties, these factors are required for viability, lie proximal to superenhancers, co-occupy regulatory elements globally, co-bind loci encoding OvCa biomarkers, and are sensitive to pharmacologic inhibition of transcription. Our predictions of MTFs, especially for tumor types with limited understanding of transcriptional drivers, pave the way to therapeutic targeting of MTFs in a broad spectrum of cancers.

Although 70-80% of newly diagnosed ovarian cancer patients respond to first-line therapy, almost all relapse and five-year survival remains below 50%. One strategy to increase five-year survival is prolonging time to relapse by improving first-line therapy response. However, no biomarker today can accurately predict individual response to therapy. In this study, we present analytical and prospective clinical validation of a new test that utilizes primary patient tissue in 3D cell culture to make patient-specific response predictions prior to initiation of treatment in the clinic. Test results were generated within seven days of tissue receipt from newly diagnosed ovarian cancer patients obtained at standard surgical debulking or laparoscopic biopsy. Patients were followed for clinical response to chemotherapy. In a study population of 44, the 32 test-predicted Responders had a clinical response rate of 100% across both adjuvant and neoadjuvant treated populations with an overall prediction accuracy of 89% (39 of 44, p < 0.0001). The test also functioned as a prognostic readout with test-predicted Responders having a significantly increased progression-free survival compared to test-predicted Non-Responders, p = 0.01. This correlative accuracy establishes the test's potential to benefit ovarian cancer patients through accurate prediction of patient-specific response before treatment.

Predicting protein levels from genotypes for proteome-wide association studies (PWAS) may provide insight into the mechanisms underlying cancer susceptibility.

To determine early somatic changes in high-grade serous ovarian cancer (HGSOC), we performed whole genome sequencing on a rare collection of 16 low stage HGSOCs. The majority showed extensive structural alterations (one had an ultramutated profile), exhibited high levels of p53 immunoreactivity, and harboured a TP53 mutation, deletion or inactivation. BRCA1 and BRCA2 mutations were observed in two tumors, with nine showing evidence of a homologous recombination (HR) defect. Combined Analysis with The Cancer Genome Atlas (TCGA) indicated that low and late stage HGSOCs have similar mutation and copy number profiles. We also found evidence that deleterious TP53 mutations are the earliest events, followed by deletions or loss of heterozygosity (LOH) of chromosomes carrying TP53, BRCA1 or BRCA2. Inactivation of HR appears to be an early event, as 62.5% of tumours showed a LOH pattern suggestive of HR defects. Three tumours with the highest ploidy had little genome-wide LOH, yet one of these had a homozygous somatic frame-shift BRCA2 mutation, suggesting that some carcinomas begin as tetraploid then descend into diploidy accompanied by genome-wide LOH. Lastly, we found evidence that structural variants (SV) cluster in HGSOC, but are absent in one ultramutated tumor, providing insights into the pathogenesis of low stage HGSOC.

Inflammation is a common feature of epithelial ovarian cancer (EOC), and measurement of plasma markers of inflammation might identify candidate markers for use in screening or presurgical evaluation of patients with adnexal masses.

We previously identified a panel of genes associated with outcome of ovarian cancer. The purpose of the current study was to assess whether variants in these genes correlated with ovarian cancer risk.