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.

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.

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.

Grade and histotype of ovarian carcinomas are often used as surrogates of molecular subtypes. We examined factors affecting pathologists' reproducibility in two prospective studies.

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.

The Nurses' Health Study (NHS), Nurses' Health Study II (NHSII), Health Professionals Follow Up Study (HPFS) and the Physicians Health Study (PHS) have collected detailed longitudinal data on multiple exposures and traits for approximately 310,000 study participants over the last 35 years. Over 160,000 study participants across the cohorts have donated a DNA sample and to date, 20,691 subjects have been genotyped as part of genome-wide association studies (GWAS) of twelve primary outcomes. However, these studies utilized six different GWAS arrays making it difficult to conduct analyses of secondary phenotypes or share controls across studies. To allow for secondary analyses of these data, we have created three new datasets merged by platform family and performed imputation using a common reference panel, the 1,000 Genomes Phase I release. Here, we describe the methodology behind the data merging and imputation and present imputation quality statistics and association results from two GWAS of secondary phenotypes (body mass index (BMI) and venous thromboembolism (VTE)). We observed the strongest BMI association for the FTO SNP rs55872725 (β = 0.45, p = 3.48x10-22), and using a significance level of p = 0.05, we replicated 19 out of 32 known BMI SNPs. For VTE, we observed the strongest association for the rs2040445 SNP (OR = 2.17, 95% CI: 1.79-2.63, p = 2.70x10-15), located downstream of F5 and also observed significant associations for the known ABO and F11 regions. This pooled resource can be used to maximize power in GWAS of phenotypes collected across the cohorts and for studying gene-environment interactions as well as rare phenotypes and genotypes.

UDP-glucuronosyltransferase (UGT) and sulfotransferase (SULT) enzymes are involved in removing sex hormones from circulation. Polymorphic variation in five UGT and SULT genes - UGT1A1 ((TA)6/(TA)7), UGT2B4 (Asp458Glu), UGT2B7 (His268Tyr), UGT2B15 (Asp85Tyr), and SULT1A1 (Arg213His)--may be associated with circulating sex hormone concentrations, or the risk of an estrogen receptor-negative (ER-) or progesterone receptor-negative (PR-) tumor.

The Cohort Consortium Vitamin D Pooling Project of Rarer Cancers (VDPP), a consortium of 10 prospective cohort studies from the United States, Finland, and China, was formed to examine the associations between circulating 25-hydroxyvitamin D (25(OH)D) concentrations and the risk of rarer cancers. Cases (total n = 5,491) included incident primary endometrial (n = 830), kidney (n = 775), ovarian (n = 516), pancreatic (n = 952), and upper gastrointestinal tract (n = 1,065) cancers and non-Hodgkin lymphoma (n = 1,353) diagnosed in the participating cohorts. At least 1 control was matched to each case on age, date of blood collection (1974-2006), sex, and race/ethnicity (n = 6,714). Covariate data were obtained from each cohort in a standardized manner. The majority of the serum or plasma samples were assayed in a central laboratory using a direct, competitive chemiluminescence immunoassay on the DiaSorin LIAISON platform (DiaSorin, Inc., Stillwater, Minnesota). Masked quality control samples included serum standards from the US National Institute of Standards and Technology. Conditional logistic regression analyses were conducted using clinically defined cutpoints, with 50-<75 nmol/L as the reference category. Meta-analyses were also conducted using inverse-variance weights in random-effects models. This consortium approach permits estimation of the association between 25(OH)D and several rarer cancers with high accuracy and precision across a wide range of 25(OH)D concentrations.

Offspring of persons exposed to childhood abuse are at higher risk of neurodevelopmental and physical health disparities across the life course. Animal experiments have indicated that paternal environmental stressors can affect sperm DNA methylation and gene expression in an offspring. Childhood abuse has been associated with epigenetic marks in human blood, saliva, and brain tissue, with statistically significant methylation differences ranging widely. However, no studies have examined the association of childhood abuse with DNA methylation in gametes. We examined the association of childhood abuse with DNA methylation in human sperm. Combined physical, emotional, and sexual abuse in childhood was characterized as none, medium, or high. DNA methylation was assayed in 46 sperm samples from 34 men in a longitudinal non-clinical cohort using HumanMethylation450 BeadChips. We performed principal component analysis and examined the correlation of principal components with abuse exposure. Childhood abuse was associated with a component that captured 6.2% of total variance in DNA methylation (p < 0.05). Next, we investigated the regions differentially methylated by abuse exposure. We identified 12 DNA regions differentially methylated by childhood abuse, containing 64 probes and including sites on genes associated with neuronal function (MAPT, CLU), fat cell regulation (PRDM16), and immune function (SDK1). We examined adulthood health behaviors, mental health, and trauma exposure as potential mediators of an association between abuse and DNAm, and found that mental health and trauma exposure partly mediated the association. Finally, we constructed a parsimonious epigenetic marker for childhood abuse using a machine learning approach, which identified three probes that predicted high vs. no childhood abuse in 71% of participants. Our results suggested that childhood abuse is associated with sperm DNA methylation, which may have implications for offspring development. Larger samples are needed to identify with greater confidence specific genomic regions differentially methylated by childhood abuse.

Endometriosis is a painful gynecologic disease affecting one in ten reproductive aged women worldwide. Few studies have correlated this symptomatology with biomarker levels among women with and without endometriosis, and no studies correlating pain with biomarker levels have been performed in young patient populations. The purpose of this study was to examine whether CA125 correlates with different types and severity of pain among adolescents and young women with and without endometriosis and assess its performance as an endometriosis biomarker among those presenting with dysmenorrhea in this young population. Reproductive-aged women with laparoscopically-confirmed endometriosis (n = 282) and controls (n = 293) who participated in The Women's Health Study: From Adolescence to Adulthood (A2A), a cohort of adolescents and young women enrolled from 2012-2018, were included in this cross-sectional analysis. Plasma CA125 values were measured using WERF EPHect compliant blood samples collected at enrollment. Average CA125 were calculated by self-reported pain type (i.e. dysmenorrhea, non-cyclic/general pelvic pain, dyspareunia), severity, and frequency in endometriosis cases and controls. Median age at blood draw was 24 years in controls and 17 years in cases, with 68% and 89% non-Hispanic white, respectively. Most endometriosis cases (95%) were rASRM stage I/II. Average CA125 values were 12.5 U/mL in controls and 12.1 U/mL in cases adjusted for age. CA125 did not differ by pain type, its severity, or frequency in endometriosis cases or controls. Among participants who reported dysmenorrhea, CA125 did not discriminate endometriosis cases from controls using cutoff of 35 U/mL (AUC = 0.51, 95%CI = 0.50-0.53). Among adolescents and young adult women, CA125 did not correlate with pain type. CA125 did not efficiently discriminate endometriosis cases from controls even when accounting for pain symptomatology. Average CA125 values were low in adolescents and young women in both endometriosis cases and controls, suggesting cautious interpretation may be needed when measuring CA125 in this population.

Although ovarian cancer is a deadly disease, approximately a third of women survive ≥9 years after diagnosis. The factors associated with achieving long-term survival are not well understood. In this study, data from the Surveillance, Epidemiology, and End Results (SEER) program were used to determine predictors of survival trajectories among women with epithelial ovarian cancer and across histotype (high-grade serous carcinoma (HGSC) and non-HGSC).

Repeated exposure to the acute proinflammatory environment that follows ovulation at the ovarian surface and distal fallopian tube over a woman's reproductive years may increase ovarian cancer risk. To address this, analyses included individual-level data from 558,709 naturally menopausal women across 20 prospective cohorts, among whom 3,246 developed invasive epithelial ovarian cancer (2,045 serous, 319 endometrioid, 184 mucinous, 121 clear cell, 577 other/unknown). Cox models were used to estimate multivariable-adjusted HRs between lifetime ovulatory cycles (LOC) and its components and ovarian cancer risk overall and by histotype. Women in the 90th percentile of LOC (>514 cycles) were almost twice as likely to be diagnosed with ovarian cancer than women in the 10th percentile (<294) [HR (95% confidence interval): 1.92 (1.60-2.30)]. Risk increased 14% per 5-year increase in LOC (60 cycles) [(1.10-1.17)]; this association remained after adjustment for LOC components: number of pregnancies and oral contraceptive use [1.08 (1.04-1.12)]. The association varied by histotype, with increased risk of serous [1.13 (1.09-1.17)], endometrioid [1.20 (1.10-1.32)], and clear cell [1.37 (1.18-1.58)], but not mucinous [0.99 (0.88-1.10), P-heterogeneity = 0.01] tumors. Heterogeneity across histotypes was reduced [P-heterogeneity = 0.15] with adjustment for LOC components [1.08 serous, 1.11 endometrioid, 1.26 clear cell, 0.94 mucinous]. Although the 10-year absolute risk of ovarian cancer is small, it roughly doubles as the number of LOC rises from approximately 300 to 500. The consistency and linearity of effects strongly support the hypothesis that each ovulation leads to small increases in the risk of most ovarian cancers, a risk that cumulates through life, suggesting this as an important area for identifying intervention strategies. SIGNIFICANCE: Although ovarian cancer is rare, risk of most ovarian cancers doubles as the number of lifetime ovulatory cycles increases from approximately 300 to 500. Thus, identifying an important area for cancer prevention research.

By iCluster analysis, we found that the integrative molecular classification of the UM was primarily driven by DNA copy number variation on chromosomes 3, 6 and 8, differential methylation and expression of genes involved in the immune system, cell morphogenesis, movement and migration, and differential mutation of genes including GNA11, BAP1, EIF1AX, SF3B1 and GNAQ. Integrative analysis revealed that pathways including IL6/JAK/STAT3 signaling, angiogenesis, allograft rejection, inflammatory response and interferon gamma response were hypomethylated and up-regulated in the M3 iSubtype, which was associated with a worse overall survival, compared to the D3 iSubtype. Using two independent gene expression datasets, we demonstrated that the subtype-driving genes had an excellent prognostic power in classifying UM into high- or low-risk groups for metastasis. Integrative analysis of UM multi-omics data provided a comprehensive view of UM biology for understanding the underlying mechanism leading to UM metastasis. The concordant molecular alterations at multi-omics levels revealed by our integrative analysis could be used for patient stratification towards personalized management and surveillance.

Epithelial ovarian, fallopian tube, and primary peritoneal cancers have shared developmental pathways. Few studies have prospectively examined heterogeneity in risk factor associations across these three anatomic sites.

Aspirin use has been associated with reduced ovarian cancer risk, yet the underlying biological mechanisms are not fully understood. To gain mechanistic insights, we assessed the association between prediagnosis low and regular-dose aspirin use and gene expression profiles in ovarian tumors.

HNF1B is overexpressed in clear cell epithelial ovarian cancer, and we observed epigenetic silencing in serous epithelial ovarian cancer, leading us to hypothesize that variation in this gene differentially associates with epithelial ovarian cancer risk according to histological subtype. Here we comprehensively map variation in HNF1B with respect to epithelial ovarian cancer risk and analyse DNA methylation and expression profiles across histological subtypes. Different single-nucleotide polymorphisms associate with invasive serous (rs7405776 odds ratio (OR)=1.13, P=3.1 × 10(-10)) and clear cell (rs11651755 OR=0.77, P=1.6 × 10(-8)) epithelial ovarian cancer. Risk alleles for the serous subtype associate with higher HNF1B-promoter methylation in these tumours. Unmethylated, expressed HNF1B, primarily present in clear cell tumours, coincides with a CpG island methylator phenotype affecting numerous other promoters throughout the genome. Different variants in HNF1B associate with risk of serous and clear cell epithelial ovarian cancer; DNA methylation and expression patterns are also notably distinct between these subtypes. These findings underscore distinct mechanisms driving different epithelial ovarian cancer histological subtypes.

Genome-wide association studies (GWAS) have so far reported 12 loci associated with serous epithelial ovarian cancer (EOC) risk. We hypothesized that some of these loci function through nearby transcription factor (TF) genes and that putative target genes of these TFs as identified by coexpression may also be enriched for additional EOC risk associations.

We hypothesized that variants in genes expressed as a consequence of interactions between ovarian cancer cells and the host micro-environment could contribute to cancer susceptibility. We therefore used a two-stage approach to evaluate common single nucleotide polymorphisms (SNPs) in 173 genes involved in stromal epithelial interactions in the Ovarian Cancer Association Consortium (OCAC). In the discovery stage, cases with epithelial ovarian cancer (n=675) and controls (n=1,162) were genotyped at 1,536 SNPs using an Illumina GoldenGate assay. Based on Positive Predictive Value estimates, three SNPs-PODXL rs1013368, ITGA6 rs13027811, and MMP3 rs522616-were selected for replication using TaqMan genotyping in up to 3,059 serous invasive cases and 8,905 controls from 16 OCAC case-control studies. An additional 18 SNPs with Pper-allele<0.05 in the discovery stage were selected for replication in a subset of five OCAC studies (n=1,233 serous invasive cases; n=3,364 controls). The discovery stage associations in PODXL, ITGA6, and MMP3 were attenuated in the larger replication set (adj. Pper-allele>or=0.5). However genotypes at TERT rs7726159 were associated with ovarian cancer risk in the smaller, five-study replication study (Pper-allele=0.03). Combined analysis of the discovery and replication sets for this TERT SNP showed an increased risk of serous ovarian cancer among non-Hispanic whites [adj. ORper-allele 1.14 (1.04-1.24) p=0.003]. Our study adds to the growing evidence that, like the 8q24 locus, the telomerase reverse transcriptase locus at 5p15.33, is a general cancer susceptibility locus.

Small-fiber polyneuropathy (SFPN) affects unmyelinated and thinly myelinated peripheral axons. Several questionnaires have been developed to assess polyneuropathy from diabetes or chemotherapy, but none for SFPN from other or unknown causes. A comprehensive survey could help clinicians diagnose and assess treatment responses, define prevalence natural history and cures, and identify research subjects. Thus, we developed the 1-page Small-Fiber Symptom Survey, using input from patients and 21 medical/scientific experts. Participants comprised consenting consecutive patients evaluated for SFPN at the Massachusetts General Hospital plus normal control subjects. Participants SFPN status was stratified on the basis of the results of their objective diagnostic tests (distal leg skin biopsy and autonomic function testing). We measured internal consistency, test retest reliability, convergent validity, and performed a receiver operating curve analysis. The 179 participants averaged 46.6 ± 15.6 years old; they were 73.2% female and 92.2% Caucasian. Eighty-five had confirmed SFPN, mostly idiopathic. Principal component analysis revealed 5 symptom clusters. The questionnaire had good internal consistency (Cronbach α = .893), excellent test retest reliability (r = .927, P < .001) and good to fair convergent validity. Participants with confirmed SFPN had more severe symptoms than others (P = .009). The Small-Fiber Symptom Survey has satisfactory psychometric properties, indicating potential future utility for surveying patient-reported symptoms of SFPN regardless of its cause.