Targeted next-generation sequencing of 21 candidate genes in hereditary ovarian cancer patients from the Republic of Bashkortostan

About 5-10% of all ovarian cancer cases show familial clustering, and some 15-25% of familial ovarian cancer cases are mediated by high-penetrance mutations in the BRCA1 and BRCA2 genes. Only few other genes have been identified for familial ovarian cancer.We conducted targeted next-generation sequencing of the protein coding region of 21 candidate genes, including UTR regions, in genomic DNA samples of 48 patients with familial ovarian cancer from the Republic of Bashkortostan. We identified deleterious variants in BRCA1, BRCA2, CHEK2, MSH6 and NBN in a total of 16 patients (33%). The NBN truncating variant, p.W143X, had not previously been reported. Seven patients (15%) were carriers of the c.5266dupC variant in BRCA1, supporting a Russian origin of this founder allele. An additional 15 variants of uncertain clinical significance were observed. We conclude that our gene panel explains about one-third of familial ovarian cancer risk in the Republic of Bashkortostan.

Individual radiosensitivity reflected by γ-H2AX and 53BP1 foci predicts outcome in PSMA-targeted radioligand therapy

PURPOSE

γ-H2AX and 53BP1 are fundamental for cellular DNA damage response (DDR) after radiation exposure and are linked to cell repair, arrest, or apoptosis. We aimed to evaluate whether DDR-markers in peripheral blood lymphocytes (PBLs) may have predictive potential for outcome in metastatic castration-resistant prostate cancer (mCRPC) patients receiving [177Lu]Lu-prostate-specific membrane antigen (PSMA) radioligand therapy (RLT).

METHODS

We prospectively enrolled 20 men with advanced mCRPC scheduled for PSMA-targeted RLT. Prior to the first cycle of [177Lu]Lu-PSMA RLT, all patients underwent [18F]F-PSMA-1007 positron emission tomography (PET)/computed tomography (CT) for assessment of tumor PSMA expression (assessing maximum standardized uptake value (SUVmax) of all tumor lesions). Blood samples were collected prior to, + 1 h after, and + 24 h after administration of [177Lu]Lu-PSMA, and DDR-markers γ-H2AX and 53BP1 were determined in PBLs through immunocytofluorescence. We then tested the predictive performance of DDR-markers relative to clinical and PET-based parameters for progressive disease (PSA-PD) after 2 cycles. In addition, the predictive value for progression-free survival (PSA-PFS, provided as median and 95% confidence interval [CI]) was explored.

RESULTS

Low baseline 53BP1 and γ-H2AX foci (P = 0.17) tended to predict early PSA-PD, whereas low SUVmax was significantly associated with higher risk for PSA-PD (P = 0.04). In Kaplan-Meier analysis, there was a trend towards prolonged PSA-PFS in patients with higher baseline 53BP1 of 6 months (mo; 95%CI, 4-9 mo) compared to 3 mo in patients with low 53BP1 (95% CI, 2-3 mo; P = 0.12). Comparable results were recorded for higher γ-H2AX expression (6 mo [95% CI, 3-9 mo] relative to 3 mo [95% CI, 2-4 mo] in patients with low γ-H2AX; P = 0.12). SUVmax, however, did not demonstrate predictive value (P = 0.29). Consistently, in univariate Cox-regression analysis, baseline 53BP1 foci demonstrated borderline significance for predicting PSA-PFS under [177Lu]Lu-PSMA RLT (P = 0.05).

CONCLUSION

In this prospective study investigating mCRPC patients undergoing [177Lu]Lu-PSMA RLT, low baseline DDR-markers in PBLs tended to predict poor outcome. Although the study group was small and results need further confirmation, these preliminary findings lay the foundation for exploring additive radiosensitizing or treatment intensification in future studies with high-risk individuals scheduled for RLT.

Assessment of γ-H2AX and 53BP1 Foci in Peripheral Blood Lymphocytes to Predict Subclinical Hematotoxicity and Response in Somatostatin Receptor-Targeted Radionuclide Therapy for Advanced Gastroenteropancreatic Neuroendocrine Tumors

BACKGROUND

We aimed to characterize γ-H2AX and 53BP1 foci formation in patients receiving somatostatin receptor-targeted radioligand therapy, and explored its role for predicting treatment-related hematotoxicity, and treatment response.

METHODS

A prospective analysis of double-strand break (DSB) markers was performed in 21 patients with advanced gastroenteropancreatic neuroendocrine tumors. γ-H2AX and 53BP1 foci formation were evaluated in peripheral blood lymphocytes (PBLs) at baseline, +1 h and +24 h after administration of 7.4 GBq (¹⁷⁷Lu)Lu-DOTA-TATE. Hematotoxicity was evaluated using standard hematology. Therapy response was assessed using (⁶⁸Ga)Ga-DOTA-TATE PET/CT before enrollment and after 2 cycles of PRRT according to the volumetric modification of RECIST 1.1.

RESULTS

DSB marker kinetics were heterogeneous among patients. Subclinical hematotoxicity was associated with γ-H2AX and 53BP1 foci formation (e.g., change in platelet count vs change in γ-H2AX⁺ cells between baseline and +1 h (r = -0.6080; p = 0.0045). Patients showing early development of new metastases had less γ-H2AX (p = 0.0125) and less 53BP1 foci per cell at +1 h (p = 0.0289), and demonstrated a distinct kinetic pattern with an absence of DSB marker decrease at +24 h (γ-H2AX: p = 0.0025; 53BP1: p = 0.0008).

CONCLUSIONS

Assessment of γ-H2AX and 53BP1 foci formation in PBLs of patients receiving radioligand therapy may hold promise for predicting subclinical hematotoxicity and early treatment response.

Involvement of Peripheral Opioid Receptors in the Realization of Food Motivation Into Eating Behavior

The involvement of peripheral opioid receptors in the mechanisms of eating behavior is still unclear. The aim of this work was to study the role of peripheral, predominantly gastric mu and delta opioid receptors in the realization of food motivation in conditions of different energy costs for eating behavior. Experiments were performed under a between-sessions progressive ratio schedule of reinforcement in food-deprived rats. The level of food motivation was calculated using a self-developed method. Food intake, motor activity, and metabolic rate were recorded in fed and hungry animals. Results showed that intragastric administration of the mu opioid receptor agonist DAMGO led to an increase in the level of food motivation in the light variant of operant feeding behaviors. Food consumption did not change. At high costs for feeding behavior, the administration of DAMGO did not alter food motivation; however, food consumption and motor activity were reduced. Intragastric administration of the delta opioid receptor agonist DADLE did not lead to changes in the level of food motivation and physical activity, but inhibition of feeding behavior was observed in all reinforcement schedules. Three regulatory pathways of eating behavior in difficult food conditions by peripheral, predominantly gastric opioid receptors are hypothesized: environmental-inhibitory afferentations and suppression of the realization of food motivation into behavior; homeostatic-inhibitory action on food motivation; and rewarding-suppression of the anticipatory reinforcement.

Assessment of the role of translationally controlled tumor protein 1 (TPT1/TCTP) in breast cancer susceptibility and ATM signaling

•

TPT1 sequencing identified one novel, potentially damaging mutation in 200 breast cancer patients.

•

TPT1 is not required for the recognition of radiation-induced DNA damage.

•

Phosphorylation of KAP1 and CHEK2 by ATM is not affected by silencing of TPT1.

•

Nuclear localization and foci formation of TPT1 potentially depends on cell type.

•

TPT1 knockdown might exert a marginally significant effect on residual γH2A.X foci.

Background and purpose

The translationally controlled tumor protein 1 (TPT1/TCTP) has been implicated in the intracellular DNA damage response. We tested the role of TPT1 in breast cancer (BC) predisposition and re-evaluated its function in Ataxia-Telangiectasia mutated (ATM)-mediated damage recognition and DNA repair.

Material and methods

The TPT1 coding sequence was scanned for mutations in genomic DNA from 200 breast cancer patients. TPT1 was down-regulated through siRNA in breast epithelial and fibroblast cell cultures. ATM activation after irradiation (IR) was analyzed by western blotting, and γH2A.X foci were monitored by immunocytochemistry.

Results

The sequencing study identified a novel, potentially damaging missense mutation in a single patient. Silencing of TPT1 did not significantly affect ATM kinase activity and did not impair the initial formation of γH2A.X foci, while we observed a marginally significant effect on residual γH2A.X foci at 6–48 h after IR.

Conclusions

TPT1 does not harbor common mutations as BC susceptibility gene. Consistently, TPT1 protein is not required for the recognition of radiation-induced DNA damage via the ATM-dependent pathway and has only slight impact on timely repair. These results may be important when considering TPT1 as a DNA damage marker.

rs495139 in the TYMS-ENOSF1 Region and Risk of Ovarian Carcinoma of Mucinous Histology

Thymidylate synthase (TYMS) is a crucial enzyme for DNA synthesis. TYMS expression is regulated by its antisense mRNA, ENOSF1. Disrupted regulation may promote uncontrolled DNA synthesis and tumor growth. We sought to replicate our previously reported association between rs495139 in the TYMS-ENOSF1 3' gene region and increased risk of mucinous ovarian carcinoma (MOC) in an independent sample. Genotypes from 24,351 controls to 15,000 women with invasive OC, including 665 MOC, were available. We estimated per-allele odds ratios (OR) and 95% confidence intervals (CI) using unconditional logistic regression, and meta-analysis when combining these data with our previous report. The association between rs495139 and MOC was not significant in the independent sample (OR = 1.09; 95% CI = 0.97⁻1.22; p = 0.15; N = 665 cases). Meta-analysis suggested a weak association (OR = 1.13; 95% CI = 1.03⁻1.24; p = 0.01; N = 1019 cases). No significant association with risk of other OC histologic types was observed (p = 0.05 for tumor heterogeneity). In expression quantitative trait locus (eQTL) analysis, the rs495139 allele was positively associated with ENOSF1 mRNA expression in normal tissues of the gastrointestinal system, particularly esophageal mucosa (r = 0.51, p = 1.7 × 10-28), and nonsignificantly in five MOC tumors. The association results, along with inconclusive tumor eQTL findings, suggest that a true effect of rs495139 might be small.

Variants in genes encoding small GTPases and association with epithelial ovarian cancer susceptibility

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.

Adult height is associated with increased risk of ovarian cancer: a Mendelian randomisation study

BACKGROUND

Observational studies suggest greater height is associated with increased ovarian cancer risk, but cannot exclude bias and/or confounding as explanations for this. Mendelian randomisation (MR) can provide evidence which may be less prone to bias.

METHODS

We pooled data from 39 Ovarian Cancer Association Consortium studies (16,395 cases; 23,003 controls). We applied two-stage predictor-substitution MR, using a weighted genetic risk score combining 609 single-nucleotide polymorphisms. Study-specific odds ratios (OR) and 95% confidence intervals (CI) for the association between genetically predicted height and risk were pooled using random-effects meta-analysis.

RESULTS

Greater genetically predicted height was associated with increased ovarian cancer risk overall (pooled-OR (pOR) = 1.06; 95% CI: 1.01-1.11 per 5 cm increase in height), and separately for invasive (pOR = 1.06; 95% CI: 1.01-1.11) and borderline (pOR = 1.15; 95% CI: 1.02-1.29) tumours.

CONCLUSIONS

Women with a genetic propensity to being taller have increased risk of ovarian cancer. This suggests genes influencing height are involved in pathways promoting ovarian carcinogenesis.

Exome sequencing and case-control analyses identify RCC1 as a candidate breast cancer susceptibility gene

Breast cancer is a genetic disease but the known genes explain a minority of cases. To elucidate the molecular basis of breast cancer in the Tunisian population, we performed exome sequencing on six BRCA1/BRCA2 mutation-negative patients with familial breast cancer and identified a novel frameshift mutation in RCC1, encoding the Regulator of Chromosome Condensation 1. Subsequent genotyping detected the 19-bp deletion in additional 5 out of 153 (3%) breast cancer patients but in none of 400 female controls (p = 0.0015). The deletion was enriched in patients with a positive family history (5%, p = 0.0009) and co-segregated with breast cancer in the initial pedigree. The mutant allele was lost in 4/6 breast tumors from mutation carriers which may be consistent with the hypothesis that RCC1 dysfunction provides a selective disadvantage at the stage of tumor progression. In summary, we propose RCC1 as a likely breast cancer susceptibility gene in the Tunisian population.

Assessment of a FBXW8 frameshift mutation, c.1312_1313delGT, in breast cancer patients and controls from Central Europe

F-box proteins participate in multiple cellular processes through ubiquitylation and subsequent degradation of target proteins, such as cyclin D1 as target of FBXW8. To investigate the spectrum of FBXW8 germ-line mutations in patients with breast cancer and healthy controls, we analyzed the whole FBXW8 coding region and flanking untranslated portions in germ-line DNA samples of 91 breast cancer patients and 277 healthy controls using next-generation amplicon sequencing. Five missense variants, one splice site variant, one frameshift variant, one synonymous variant, and one variant in the 3'-UTR were identified. Frameshift mutation FBXW8 c.1312_1313delGT was considered functionally relevant and was investigated for its potential association with breast cancer risk through subsequent genotyping in two hospital-based breast cancer case-control series from Belarus and Germany, respectively, comprising a total of 2740 breast cancer cases and 2174 controls. The mutation was found in 30 cases and 23 controls with an adjusted Odds Ratio 1.02 (95% CI 0.59-1.77; p = 0.94) in the combined analysis. There was no statistically significant difference when patients were stratified by ER status, PR status, age at diagnosis, ductal histology, contralateral status, family history or tumor grade. Altogether, our data exclude clinically actionable breast cancer risks above two-fold for the FBXW8 c.1312_1313delGT mutation, although larger studies would be needed to exclude low-penetrance associations.

Analysis of a RECQL splicing mutation, c.1667_1667+3delAGTA, in breast cancer patients and controls from Central Europe

RECQL is a DNA helicase required for genomic stability. Two studies have recently identified RECQL as a novel breast cancer susceptibility gene. The most common RECQL mutation, the 4 bp-deletion c.1667_1667+3delAGTA, was five-fold enriched in Polish breast cancer patients, but the exact magnitude of the risk is uncertain. We investigated two hospital-based breast cancer case-control series from Belarus and Germany, respectively, comprising a total of 2596 breast cancer patients and 2132 healthy females. The mutation was found in 9 cases and 6 controls, with an adjusted Odds Ratio 1.23 (95% CI 0.44-3.47; p = 0.69) in the combined analysis. Among the cases, heterozygosity for c.1667_1667+3delAGTA was linked with estrogen-receptor positive breast cancer. There was no significant difference in age at diagnosis between carriers and non-carriers, and only one of the carriers reported a first-degree family history. Meta-analysis with the initial study from Poland suggests an about two-fold increase in risk for this mutation (OR 2.51; 95% CI 1.13-5.57, p = 0.02). Altogether, the data indicate that RECQL* c.1667_1667+3delAGTA is not a high-risk mutation for breast cancer though it could represent a moderate-risk breast cancer susceptibility allele. Further studies will be required to determine the clinical significance of testing for this RECQL mutation.

Germline whole exome sequencing and large-scale replication identifies FANCM as a likely high grade serous ovarian cancer susceptibility gene

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.

Identification of 12 new susceptibility loci for different histotypes of epithelial ovarian cancer

To identify common alleles associated with different histotypes of epithelial ovarian cancer (EOC), we pooled data from multiple genome-wide genotyping projects totaling 25,509 EOC cases and 40,941 controls. We identified nine new susceptibility loci for different EOC histotypes: six for serous EOC histotypes (3q28, 4q32.3, 8q21.11, 10q24.33, 18q11.2 and 22q12.1), two for mucinous EOC (3q22.3 and 9q31.1) and one for endometrioid EOC (5q12.3). We then performed meta-analysis on the results for high-grade serous ovarian cancer with the results from analysis of 31,448 BRCA1 and BRCA2 mutation carriers, including 3,887 mutation carriers with EOC. This identified three additional susceptibility loci at 2q13, 8q24.1 and 12q24.31. Integrated analyses of genes and regulatory biofeatures at each locus predicted candidate susceptibility genes, including OBFC1, a new candidate susceptibility gene for low-grade and borderline serous EOC.

Enrichment of putative PAX8 target genes at serous epithelial ovarian cancer susceptibility loci

BACKGROUND

Genome-wide association studies (GWAS) have identified 18 loci associated with serous ovarian cancer (SOC) susceptibility but the biological mechanisms driving these findings remain poorly characterised. Germline cancer risk loci may be enriched for target genes of transcription factors (TFs) critical to somatic tumorigenesis.

METHODS

All 615 TF-target sets from the Molecular Signatures Database were evaluated using gene set enrichment analysis (GSEA) and three GWAS for SOC risk: discovery (2196 cases/4396 controls), replication (7035 cases/21 693 controls; independent from discovery), and combined (9627 cases/30 845 controls; including additional individuals).

RESULTS

The PAX8-target gene set was ranked 1/615 in the discovery (PGSEA<0.001; FDR=0.21), 7/615 in the replication (PGSEA=0.004; FDR=0.37), and 1/615 in the combined (PGSEA<0.001; FDR=0.21) studies. Adding other genes reported to interact with PAX8 in the literature to the PAX8-target set and applying an alternative to GSEA, interval enrichment, further confirmed this association (P=0.006). Fifteen of the 157 genes from this expanded PAX8 pathway were near eight loci associated with SOC risk at P<10-5 (including six with P<5 × 10-8). The pathway was also associated with differential gene expression after shRNA-mediated silencing of PAX8 in HeyA8 (PGSEA=0.025) and IGROV1 (PGSEA=0.004) SOC cells and several PAX8 targets near SOC risk loci demonstrated in vitro transcriptomic perturbation.

CONCLUSIONS

Putative PAX8 target genes are enriched for common SOC risk variants. This finding from our agnostic evaluation is of particular interest given that PAX8 is well-established as a specific marker for the cell of origin of SOC.

Assessment of an APOBEC3B truncating mutation, c.783delG, in patients with breast cancer

PURPOSE

APOBEC3B belongs to the family of DNA-editing enzymes. A copy number variant targeting the genomic APOBEC3A-APOBEC3B locus has a significant impact on breast cancer risk, but the relative contribution of APOBEC3B is uncertain. In this study, we investigate a loss-of-function mutation that selectively targets APOBEC3B, for its association with breast cancer risk.

METHODS

We performed exome sequencing on genomic DNA samples of 6 Byelorussian patients with familial breast cancer. We then studied through mutation-specific genotyping four hospital-based breast cancer case-control series from Belarus, Russia, Germany, and Iran, respectively, comprising a total of 3070 breast cancer patients and 2878 healthy females. Results were evaluated using fixed-effects meta-analyses.

RESULTS

Exome sequencing uncovered a frameshift mutation, APOBEC3B*c.783delG, that was recurrent in the study populations. Subsequent genotyping identified this mutation in 23 additional breast cancer cases and 9 healthy female controls, with an adjusted Odds Ratio 2.29 (95% CI 1.04; 5.03, P = 0.04) in the combined analysis. There was an enrichment of the c.783delG mutation in patients with breast cancer diagnosed below 50 years of age (OR 3.22, 95% CI 1.37; 7.56, P = 0.007).

CONCLUSIONS

APOBEC3B*c.783delG showed evidence of modest association with breast cancer and seemed to contribute to earlier onset of the disease. These results may need to be reconciled with proposals to consider APOBEC3B as a possible therapeutic target in breast cancer.

Inherited variants affecting RNA editing may contribute to ovarian cancer susceptibility: results from a large-scale collaboration

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.

Assessment of variation in immunosuppressive pathway genes reveals TGFBR2 to be associated with risk of clear cell ovarian cancer

BACKGROUND

Regulatory T (Treg) cells, a subset of CD4+ T lymphocytes, are mediators of immunosuppression in cancer, and, thus, variants in genes encoding Treg cell immune molecules could be associated with ovarian cancer.

METHODS

In a population of 15,596 epithelial ovarian cancer (EOC) cases and 23,236 controls, we measured genetic associations of 1,351 SNPs in Treg cell pathway genes with odds of ovarian cancer and tested pathway and gene-level associations, overall and by histotype, for the 25 genes, using the admixture likelihood (AML) method. The most significant single SNP associations were tested for correlation with expression levels in 44 ovarian cancer patients.

RESULTS

The most significant global associations for all genes in the pathway were seen in endometrioid ( p = 0.082) and clear cell ( p = 0.083), with the most significant gene level association seen with TGFBR2 ( p = 0.001) and clear cell EOC. Gene associations with histotypes at p < 0.05 included: IL12 ( p = 0.005 and p = 0.008, serous and high-grade serous, respectively), IL8RA ( p = 0.035, endometrioid and mucinous), LGALS1 ( p = 0.03, mucinous), STAT5B ( p = 0.022, clear cell), TGFBR1 ( p = 0.021 endometrioid) and TGFBR2 ( p = 0.017 and p = 0.025, endometrioid and mucinous, respectively).

CONCLUSIONS

Common inherited gene variation in Treg cell pathways shows some evidence of germline genetic contribution to odds of EOC that varies by histologic subtype and may be associated with mRNA expression of immune-complex receptor in EOC patients.

Association of vitamin D levels and risk of ovarian cancer: a Mendelian randomization study

BACKGROUND

In vitro and observational epidemiological studies suggest that vitamin D may play a role in cancer prevention. However, the relationship between vitamin D and ovarian cancer is uncertain, with observational studies generating conflicting findings. A potential limitation of observational studies is inadequate control of confounding. To overcome this problem, we used Mendelian randomization (MR) to evaluate the association between single nucleotide polymorphisms (SNPs) associated with circulating 25-hydroxyvitamin D [25(OH)D] concentration and risk of ovarian cancer.

METHODS

We employed SNPs with well-established associations with 25(OH)D concentration as instrumental variables for MR: rs7944926 (DHCR7), rs12794714 (CYP2R1) and rs2282679 (GC). We included 31 719 women of European ancestry (10 065 cases, 21 654 controls) from the Ovarian Cancer Association Consortium, who were genotyped using customized Illumina Infinium iSelect (iCOGS) arrays. A two-sample (summary data) MR approach was used and analyses were performed separately for all ovarian cancer (10 065 cases) and for high-grade serous ovarian cancer (4121 cases).

RESULTS

The odds ratio for epithelial ovarian cancer risk (10 065 cases) estimated by combining the individual SNP associations using inverse variance weighting was 1.27 (95% confidence interval: 1.06 to 1.51) per 20 nmol/L decrease in 25(OH)D concentration. The estimated odds ratio for high-grade serous epithelial ovarian cancer (4121 cases) was 1.54 (1.19, 2.01).

CONCLUSIONS

Genetically lowered 25-hydroxyvitamin D concentrations were associated with higher ovarian cancer susceptibility in Europeans. These findings suggest that increasing plasma vitamin D levels may reduce risk of ovarian cancer.

Functional mechanisms underlying pleiotropic risk alleles at the 19p13.1 breast-ovarian cancer susceptibility locus

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.

Genome-Wide Meta-Analyses of Breast, Ovarian, and Prostate Cancer Association Studies Identify Multiple New Susceptibility Loci Shared by at Least Two Cancer Types

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.

Assessing the genetic architecture of epithelial ovarian cancer histological subtypes

Epithelial ovarian cancer (EOC) is one of the deadliest common cancers. The five most common types of disease are high-grade and low-grade serous, endometrioid, mucinous and clear cell carcinoma. Each of these subtypes present distinct molecular pathogeneses and sensitivities to treatments. Recent studies show that certain genetic variants confer susceptibility to all subtypes while other variants are subtype-specific. Here, we perform an extensive analysis of the genetic architecture of EOC subtypes. To this end, we used data of 10,014 invasive EOC patients and 21,233 controls from the Ovarian Cancer Association Consortium genotyped in the iCOGS array (211,155 SNPs). We estimate the array heritability (attributable to variants tagged on arrays) of each subtype and their genetic correlations. We also look for genetic overlaps with factors such as obesity, smoking behaviors, diabetes, age at menarche and height. We estimated the array heritabilities of high-grade serous disease ([Formula: see text] = 8.8 ± 1.1 %), endometrioid ([Formula: see text] = 3.2 ± 1.6 %), clear cell ([Formula: see text] = 6.7 ± 3.3 %) and all EOC ([Formula: see text] = 5.6 ± 0.6 %). Known associated loci contributed approximately 40 % of the total array heritability for each subtype. The contribution of each chromosome to the total heritability was not proportional to chromosome size. Through bivariate and cross-trait LD score regression, we found evidence of shared genetic backgrounds between the three high-grade subtypes: serous, endometrioid and undifferentiated. Finally, we found significant genetic correlations of all EOC with diabetes and obesity using a polygenic prediction approach.

Adult body mass index and risk of ovarian cancer by subtype: a Mendelian randomization study

BACKGROUND

Observational studies have reported a positive association between body mass index (BMI) and ovarian cancer risk. However, questions remain as to whether this represents a causal effect, or holds for all histological subtypes. The lack of association observed for serous cancers may, for instance, be due to disease-associated weight loss. Mendelian randomization (MR) uses genetic markers as proxies for risk factors to overcome limitations of observational studies. We used MR to elucidate the relationship between BMI and ovarian cancer, hypothesizing that genetically predicted BMI would be associated with increased risk of non-high grade serous ovarian cancers (non-HGSC) but not HGSC.

METHODS

We pooled data from 39 studies (14 047 cases, 23 003 controls) in the Ovarian Cancer Association Consortium. We constructed a weighted genetic risk score (GRS, partial F-statistic = 172), summing alleles at 87 single nucleotide polymorphisms previously associated with BMI, weighting by their published strength of association with BMI. Applying two-stage predictor-substitution MR, we used logistic regression to estimate study-specific odds ratios (OR) and 95% confidence intervals (CI) for the association between genetically predicted BMI and risk, and pooled these using random-effects meta-analysis.

RESULTS

Higher genetically predicted BMI was associated with increased risk of non-HGSC (pooled OR = 1.29, 95% CI 1.03-1.61 per 5 units BMI) but not HGSC (pooled OR = 1.06, 95% CI 0.88-1.27). Secondary analyses stratified by behaviour/subtype suggested that, consistent with observational data, the association was strongest for low-grade/borderline serous cancers (OR = 1.93, 95% CI 1.33-2.81).

CONCLUSIONS

Our data suggest that higher BMI increases risk of non-HGSC, but not the more common and aggressive HGSC subtype, confirming the observational evidence.

Rare ATAD5 missense variants in breast and ovarian cancer patients

ATAD5/ELG1 is a protein crucially involved in replication and maintenance of genome stability. ATAD5 has recently been identified as a genomic risk locus for both breast and ovarian cancer through genome-wide association studies. We aimed to investigate the spectrum of coding ATAD5 germ-line mutations in hospital-based series of patients with triple-negative breast cancer or serous ovarian cancer compared with healthy controls. The ATAD5 coding and adjacent splice site regions were analyzed by targeted next-generation sequencing of DNA samples from 273 cancer patients, including 114 patients with triple-negative breast cancer and 159 patients with serous epithelial ovarian cancer, and from 276 healthy females. Among 42 different variants identified, twenty-two were rare missense substitutions, of which 14 were classified as pathogenic by at least one in silico prediction tool. Three of four novel missense substitutions (p.S354I, p.H974R and p.K1466N) were predicted to be pathogenic and were all identified in ovarian cancer patients. Overall, rare missense variants with predicted pathogenicity tended to be enriched in ovarian cancer patients (14/159) versus controls (11/276) (p = 0.05, 2df). While truncating germ-line variants in ATAD5 were not detected, it remains possible that several rare missense variants contribute to genetic susceptibility toward epithelial ovarian carcinomas.

PPM1D Mosaic Truncating Variants in Ovarian Cancer Cases May Be Treatment-Related Somatic Mutations

Mosaic truncating mutations in the protein phosphatase, Mg(2+)/Mn(2+)-dependent, 1D (PPM1D) gene have recently been reported with a statistically significantly greater frequency in lymphocyte DNA from ovarian cancer case patients compared with unaffected control patients. Using massively parallel sequencing (MPS) we identified truncating PPM1D mutations in 12 of 3236 epithelial ovarian cancer (EOC) case patients (0.37%) but in only one of 3431 unaffected control patients (0.03%) (P = .001). All statistical tests were two-sided. A combination of Sanger sequencing, pyrosequencing, and MPS data suggested that 12 of the 13 mutations were mosaic. All mutations were identified in post-chemotherapy treatment blood samples from case patients (n = 1827) (average 1234 days post-treatment in carriers) rather than from cases collected pretreatment (less than 14 days after diagnosis, n = 1384) (P = .002). These data suggest that PPM1D variants in EOC cases are primarily somatic mosaic mutations caused by treatment and are not associated with germline predisposition to EOC.

Epithelial-Mesenchymal Transition (EMT) Gene Variants and Epithelial Ovarian Cancer (EOC) Risk

Epithelial-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.

Height and Breast Cancer Risk: Evidence From Prospective Studies and Mendelian Randomization

BACKGROUND

Epidemiological studies have linked adult height with breast cancer risk in women. However, the magnitude of the association, particularly by subtypes of breast cancer, has not been established. Furthermore, the mechanisms of the association remain unclear.

METHODS

We performed a meta-analysis to investigate associations between height and breast cancer risk using data from 159 prospective cohorts totaling 5216302 women, including 113178 events. In a consortium with individual-level data from 46325 case patients and 42482 control patients, we conducted a Mendelian randomization analysis using a genetic score that comprised 168 height-associated variants as an instrument. This association was further evaluated in a second consortium using summary statistics data from 16003 case patients and 41335 control patients.

RESULTS

The pooled relative risk of breast cancer was 1.17 (95% confidence interval [CI] = 1.15 to 1.19) per 10cm increase in height in the meta-analysis of prospective studies. In Mendelian randomization analysis, the odds ratio of breast cancer per 10cm increase in genetically predicted height was 1.22 (95% CI = 1.13 to 1.32) in the first consortium and 1.21 (95% CI = 1.05 to 1.39) in the second consortium. The association was found in both premenopausal and postmenopausal women but restricted to hormone receptor-positive breast cancer. Analyses of height-associated variants identified eight new loci associated with breast cancer risk after adjusting for multiple comparisons, including three loci at 1q21.2, DNAJC27, and CCDC91 at genome-wide significance level P < 5×10(-8).

CONCLUSIONS

Our study provides strong evidence that adult height is a risk factor for breast cancer in women and certain genetic factors and biological pathways affecting adult height have an important role in the etiology of breast cancer.

Germline Mutations in the BRIP1, BARD1, PALB2, and NBN Genes in Women With Ovarian Cancer

BACKGROUND

Epithelial ovarian cancer (EOC) is the most lethal gynecological malignancy, responsible for 13 000 deaths per year in the United States. Risk prediction based on identifying germline mutations in ovarian cancer susceptibility genes could have a clinically significant impact on reducing disease mortality.

METHODS

Next generation sequencing was used to identify germline mutations in the coding regions of four candidate susceptibility genes-BRIP1, BARD1, PALB2 and NBN-in 3236 invasive EOC case patients and 3431 control patients of European origin, and in 2000 unaffected high-risk women from a clinical screening trial of ovarian cancer (UKFOCSS). For each gene, we estimated the prevalence and EOC risks and evaluated associations between germline variant status and clinical and epidemiological risk factor information. All statistical tests were two-sided.

RESULTS

We found an increased frequency of deleterious mutations in BRIP1 in case patients (0.9%) and in the UKFOCSS participants (0.6%) compared with control patients (0.09%) (P = 1 x 10(-4) and 8 x 10(-4), respectively), but no differences for BARD1 (P = .39), NBN1 ( P = .61), or PALB2 (P = .08). There was also a difference in the frequency of rare missense variants in BRIP1 between case patients and control patients (P = 5.5 x 10(-4)). The relative risks associated with BRIP1 mutations were 11.22 for invasive EOC (95% confidence interval [CI] = 3.22 to 34.10, P = 1 x 10(-4)) and 14.09 for high-grade serous disease (95% CI = 4.04 to 45.02, P = 2 x 10(-5)). Segregation analysis in families estimated the average relative risks in BRIP1 mutation carriers compared with the general population to be 3.41 (95% CI = 2.12 to 5.54, P = 7×10(-7)).

CONCLUSIONS

Deleterious germline mutations in BRIP1 are associated with a moderate increase in EOC risk. These data have clinical implications for risk prediction and prevention approaches for ovarian cancer and emphasize the critical need for risk estimates based on very large sample sizes before genes of moderate penetrance have clinical utility in cancer prevention.

Fine-scale mapping of the 4q24 locus identifies two independent loci associated with breast cancer risk

BACKGROUND

A recent association study identified a common variant (rs9790517) at 4q24 to be associated with breast cancer risk. Independent association signals and potential functional variants in this locus have not been explored.

METHODS

We conducted a fine-mapping analysis in 55,540 breast cancer cases and 51,168 controls from the Breast Cancer Association Consortium.

RESULTS

Conditional analyses identified two independent association signals among women of European ancestry, represented by rs9790517 [conditional P = 2.51 × 10(-4); OR, 1.04; 95% confidence interval (CI), 1.02-1.07] and rs77928427 (P = 1.86 × 10(-4); OR, 1.04; 95% CI, 1.02-1.07). Functional annotation using data from the Encyclopedia of DNA Elements (ENCODE) project revealed two putative functional variants, rs62331150 and rs73838678 in linkage disequilibrium (LD) with rs9790517 (r(2) ≥ 0.90) residing in the active promoter or enhancer, respectively, of the nearest gene, TET2. Both variants are located in DNase I hypersensitivity and transcription factor-binding sites. Using data from both The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC), we showed that rs62331150 was associated with level of expression of TET2 in breast normal and tumor tissue.

CONCLUSION

Our study identified two independent association signals at 4q24 in relation to breast cancer risk and suggested that observed association in this locus may be mediated through the regulation of TET2.

IMPACT

Fine-mapping study with large sample size warranted for identification of independent loci for breast cancer risk.

Common variants at the CHEK2 gene locus and risk of epithelial ovarian cancer

Genome-wide association studies have identified 20 genomic regions associated with risk of epithelial ovarian cancer (EOC), but many additional risk variants may exist. Here, we evaluated associations between common genetic variants [single nucleotide polymorphisms (SNPs) and indels] in DNA repair genes and EOC risk. We genotyped 2896 common variants at 143 gene loci in DNA samples from 15 397 patients with invasive EOC and controls. We found evidence of associations with EOC risk for variants at FANCA, EXO1, E2F4, E2F2, CREB5 and CHEK2 genes (P ≤ 0.001). The strongest risk association was for CHEK2 SNP rs17507066 with serous EOC (P = 4.74 x 10(-7)). Additional genotyping and imputation of genotypes from the 1000 genomes project identified a slightly more significant association for CHEK2 SNP rs6005807 (r (2) with rs17507066 = 0.84, odds ratio (OR) 1.17, 95% CI 1.11-1.24, P = 1.1×10(-7)). We identified 293 variants in the region with likelihood ratios of less than 1:100 for representing the causal variant. Functional annotation identified 25 candidate SNPs that alter transcription factor binding sites within regulatory elements active in EOC precursor tissues. In The Cancer Genome Atlas dataset, CHEK2 gene expression was significantly higher in primary EOCs compared to normal fallopian tube tissues (P = 3.72×10(-8)). We also identified an association between genotypes of the candidate causal SNP rs12166475 (r (2) = 0.99 with rs6005807) and CHEK2 expression (P = 2.70×10(-8)). These data suggest that common variants at 22q12.1 are associated with risk of serous EOC and CHEK2 as a plausible target susceptibility gene.

Shared genetics underlying epidemiological association between endometriosis and ovarian cancer

Epidemiological studies have demonstrated associations between endometriosis and certain histotypes of ovarian cancer, including clear cell, low-grade serous and endometrioid carcinomas. We aimed to determine whether the observed associations might be due to shared genetic aetiology. To address this, we used two endometriosis datasets genotyped on common arrays with full-genome coverage (3194 cases and 7060 controls) and a large ovarian cancer dataset genotyped on the customized Illumina Infinium iSelect (iCOGS) arrays (10 065 cases and 21 663 controls). Previous work has suggested that a large number of genetic variants contribute to endometriosis and ovarian cancer (all histotypes combined) susceptibility. Here, using the iCOGS data, we confirmed polygenic architecture for most histotypes of ovarian cancer. This led us to evaluate if the polygenic effects are shared across diseases. We found evidence for shared genetic risks between endometriosis and all histotypes of ovarian cancer, except for the intestinal mucinous type. Clear cell carcinoma showed the strongest genetic correlation with endometriosis (0.51, 95% CI = 0.18-0.84). Endometrioid and low-grade serous carcinomas had similar correlation coefficients (0.48, 95% CI = 0.07-0.89 and 0.40, 95% CI = 0.05-0.75, respectively). High-grade serous carcinoma, which often arises from the fallopian tubes, showed a weaker genetic correlation with endometriosis (0.25, 95% CI = 0.11-0.39), despite the absence of a known epidemiological association. These results suggest that the epidemiological association between endometriosis and ovarian adenocarcinoma may be attributable to shared genetic susceptibility loci.

Network-Based Integration of GWAS and Gene Expression Identifies a HOX-Centric Network Associated with Serous Ovarian Cancer Risk

BACKGROUND

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.

METHODS

We selected TF genes within 1 Mb of the top signal at the 12 genome-wide significant risk loci. Mutual information, a form of correlation, was used to build networks of genes strongly coexpressed with each selected TF gene in the unified microarray dataset of 489 serous EOC tumors from The Cancer Genome Atlas. Genes represented in this dataset were subsequently ranked using a gene-level test based on results for germline SNPs from a serous EOC GWAS meta-analysis (2,196 cases/4,396 controls).

RESULTS

Gene set enrichment analysis identified six networks centered on TF genes (HOXB2, HOXB5, HOXB6, HOXB7 at 17q21.32 and HOXD1, HOXD3 at 2q31) that were significantly enriched for genes from the risk-associated end of the ranked list (P < 0.05 and FDR < 0.05). These results were replicated (P < 0.05) using an independent association study (7,035 cases/21,693 controls). Genes underlying enrichment in the six networks were pooled into a combined network.

CONCLUSION

We identified a HOX-centric network associated with serous EOC risk containing several genes with known or emerging roles in serous EOC development.

IMPACT

Network analysis integrating large, context-specific datasets has the potential to offer mechanistic insights into cancer susceptibility and prioritize genes for experimental characterization.

Cis-eQTL analysis and functional validation of candidate susceptibility genes for high-grade serous ovarian cancer

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.

Common Genetic Variation in Circadian Rhythm Genes and Risk of Epithelial Ovarian Cancer (EOC)

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.

Common Genetic Variation In Cellular Transport Genes and Epithelial Ovarian Cancer (EOC) Risk

BACKGROUND

Defective cellular transport processes can lead to aberrant accumulation of trace elements, iron, small molecules and hormones in the cell, which in turn may promote the formation of reactive oxygen species, promoting DNA damage and aberrant expression of key regulatory cancer genes. As DNA damage and uncontrolled proliferation are hallmarks of cancer, including epithelial ovarian cancer (EOC), we hypothesized that inherited variation in the cellular transport genes contributes to EOC risk.

METHODS

In total, DNA samples were obtained from 14,525 case subjects with invasive EOC and from 23,447 controls from 43 sites in the Ovarian Cancer Association Consortium (OCAC). Two hundred seventy nine SNPs, representing 131 genes, were genotyped using an Illumina Infinium iSelect BeadChip as part of the Collaborative Oncological Gene-environment Study (COGS). SNP analyses were conducted using unconditional logistic regression under a log-additive model, and the FDR q<0.2 was applied to adjust for multiple comparisons.

RESULTS

The most significant evidence of an association for all invasive cancers combined and for the serous subtype was observed for SNP rs17216603 in the iron transporter gene HEPH (invasive: OR = 0.85, P = 0.00026; serous: OR = 0.81, P = 0.00020); this SNP was also associated with the borderline/low malignant potential (LMP) tumors (P = 0.021). Other genes significantly associated with EOC histological subtypes (p<0.05) included the UGT1A (endometrioid), SLC25A45 (mucinous), SLC39A11 (low malignant potential), and SERPINA7 (clear cell carcinoma). In addition, 1785 SNPs in six genes (HEPH, MGST1, SERPINA, SLC25A45, SLC39A11 and UGT1A) were imputed from the 1000 Genomes Project and examined for association with INV EOC in white-European subjects. The most significant imputed SNP was rs117729793 in SLC39A11 (per allele, OR = 2.55, 95% CI = 1.5-4.35, p = 5.66x10-4).

CONCLUSION

These results, generated on a large cohort of women, revealed associations between inherited cellular transport gene variants and risk of EOC histologic subtypes.

Fine-mapping identifies two additional breast cancer susceptibility loci at 9q31.2

We recently identified a novel susceptibility variant, rs865686, for estrogen-receptor positive breast cancer at 9q31.2. Here, we report a fine-mapping analysis of the 9q31.2 susceptibility locus using 43 160 cases and 42 600 controls of European ancestry ascertained from 52 studies and a further 5795 cases and 6624 controls of Asian ancestry from nine studies. Single nucleotide polymorphism (SNP) rs676256 was most strongly associated with risk in Europeans (odds ratios [OR] = 0.90 [0.88-0.92]; P-value = 1.58 × 10(-25)). This SNP is one of a cluster of highly correlated variants, including rs865686, that spans ∼14.5 kb. We identified two additional independent association signals demarcated by SNPs rs10816625 (OR = 1.12 [1.08-1.17]; P-value = 7.89 × 10(-09)) and rs13294895 (OR = 1.09 [1.06-1.12]; P-value = 2.97 × 10(-11)). SNP rs10816625, but not rs13294895, was also associated with risk of breast cancer in Asian individuals (OR = 1.12 [1.06-1.18]; P-value = 2.77 × 10(-05)). Functional genomic annotation using data derived from breast cancer cell-line models indicates that these SNPs localise to putative enhancer elements that bind known drivers of hormone-dependent breast cancer, including ER-α, FOXA1 and GATA-3. In vitro analyses indicate that rs10816625 and rs13294895 have allele-specific effects on enhancer activity and suggest chromatin interactions with the KLF4 gene locus. These results demonstrate the power of dense genotyping in large studies to identify independent susceptibility variants. Analysis of associations using subjects with different ancestry, combined with bioinformatic and genomic characterisation, can provide strong evidence for the likely causative alleles and their functional basis.

Evaluating the ovarian cancer gonadotropin hypothesis: a candidate gene study

OBJECTIVE

Ovarian cancer is a hormone-related disease with a strong genetic basis. However, none of its high-penetrance susceptibility genes and GWAS-identified variants to date are known to be involved in hormonal pathways. Given the hypothesized etiologic role of gonadotropins, an assessment of how variability in genes involved in the gonadotropin signaling pathway impacts disease risk is warranted.

METHODS

Genetic data from 41 ovarian cancer study sites were pooled and unconditional logistic regression was used to evaluate whether any of the 2185 SNPs from 11 gonadotropin signaling pathway genes was associated with ovarian cancer risk. A burden test using the admixture likelihood (AML) method was also used to evaluate gene-level associations.

RESULTS

We did not find any genome-wide significant associations between individual SNPs and ovarian cancer risk. However, there was some suggestion of gene-level associations for four gonadotropin signaling pathway genes: INHBB (p=0.045, mucinous), LHCGR (p=0.046, high-grade serous), GNRH (p=0.041, high-grade serous), and FSHB (p=0.036, overall invasive). There was also suggestive evidence for INHA (p=0.060, overall invasive).

CONCLUSIONS

Ovarian cancer studies have limited sample numbers, thus fewer genome-wide susceptibility alleles, with only modest associations, have been identified relative to breast and prostate cancers. We have evaluated the majority of ovarian cancer studies with biological samples, to our knowledge, leaving no opportunity for replication. Using both our understanding of biology and powerful gene-level tests, we have identified four putative ovarian cancer loci near INHBB, LHCGR, GNRH, and FSHB that warrant a second look if larger sample sizes and denser genotype chips become available.

Identification of six new susceptibility loci for invasive epithelial ovarian cancer

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.

Inherited variants in the inner centromere protein (INCENP) gene of the chromosomal passenger complex contribute to the susceptibility of ER-negative breast cancer

The chromosomal passenger complex (CPC) plays a pivotal role in the regulation of cell division. Therefore, inherited CPC variability could influence tumor development. The present candidate gene approach investigates the relationship between single nucleotide polymorphisms (SNPs) in genes encoding key CPC components and breast cancer risk. Fifteen SNPs in four CPC genes (INCENP, AURKB, BIRC5 and CDCA8) were genotyped in 88 911 European women from 39 case-control studies of the Breast Cancer Association Consortium. Possible associations were investigated in fixed-effects meta-analyses. The synonymous SNP rs1675126 in exon 7 of INCENP was associated with overall breast cancer risk [per A allele odds ratio (OR) 0.95, 95% confidence interval (CI) 0.92-0.98, P = 0.007] and particularly with estrogen receptor (ER)-negative breast tumors (per A allele OR 0.89, 95% CI 0.83-0.95, P = 0.0005). SNPs not directly genotyped were imputed based on 1000 Genomes. The SNPs rs1047739 in the 3' untranslated region and rs144045115 downstream of INCENP showed the strongest association signals for overall (per T allele OR 1.03, 95% CI 1.00-1.06, P = 0.0009) and ER-negative breast cancer risk (per A allele OR 1.06, 95% CI 1.02-1.10, P = 0.0002). Two genotyped SNPs in BIRC5 were associated with familial breast cancer risk (top SNP rs2071214: per G allele OR 1.12, 95% CI 1.04-1.21, P = 0.002). The data suggest that INCENP in the CPC pathway contributes to ER-negative breast cancer susceptibility in the European population. In spite of a modest contribution of CPC-inherited variants to the total burden of sporadic and familial breast cancer, their potential as novel targets for breast cancer treatment should be further investigated.

Refined histopathological predictors of BRCA1 and BRCA2 mutation status: a large-scale analysis of breast cancer characteristics from the BCAC, CIMBA, and ENIGMA consortia

INTRODUCTION

The distribution of histopathological features of invasive breast tumors in BRCA1 or BRCA2 germline mutation carriers differs from that of individuals with no known mutation. Histopathological features thus have utility for mutation prediction, including statistical modeling to assess pathogenicity of BRCA1 or BRCA2 variants of uncertain clinical significance. We analyzed large pathology datasets accrued by the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA) and the Breast Cancer Association Consortium (BCAC) to reassess histopathological predictors of BRCA1 and BRCA2 mutation status, and provide robust likelihood ratio (LR) estimates for statistical modeling.

METHODS

Selection criteria for study/center inclusion were estrogen receptor (ER) status or grade data available for invasive breast cancer diagnosed younger than 70 years. The dataset included 4,477 BRCA1 mutation carriers, 2,565 BRCA2 mutation carriers, and 47,565 BCAC breast cancer cases. Country-stratified estimates of the likelihood of mutation status by histopathological markers were derived using a Mantel-Haenszel approach.

RESULTS

ER-positive phenotype negatively predicted BRCA1 mutation status, irrespective of grade (LRs from 0.08 to 0.90). ER-negative grade 3 histopathology was more predictive of positive BRCA1 mutation status in women 50 years or older (LR = 4.13 (3.70 to 4.62)) versus younger than 50 years (LR = 3.16 (2.96 to 3.37)). For BRCA2, ER-positive grade 3 phenotype modestly predicted positive mutation status irrespective of age (LR = 1.7-fold), whereas ER-negative grade 3 features modestly predicted positive mutation status at 50 years or older (LR = 1.54 (1.27 to 1.88)). Triple-negative tumor status was highly predictive of BRCA1 mutation status for women younger than 50 years (LR = 3.73 (3.43 to 4.05)) and 50 years or older (LR = 4.41 (3.86 to 5.04)), and modestly predictive of positive BRCA2 mutation status in women 50 years or older (LR = 1.79 (1.42 to 2.24)).

CONCLUSIONS

These results refine likelihood-ratio estimates for predicting BRCA1 and BRCA2 mutation status by using commonly measured histopathological features. Age at diagnosis is an important variable for most analyses, and grade is more informative than ER status for BRCA2 mutation carrier prediction. The estimates will improve BRCA1 and BRCA2 variant classification and inform patient mutation testing and clinical management.

Common non-synonymous SNPs associated with breast cancer susceptibility: findings from the Breast Cancer Association Consortium

Candidate variant association studies have been largely unsuccessful in identifying common breast cancer susceptibility variants, although most studies have been underpowered to detect associations of a realistic magnitude. We assessed 41 common non-synonymous single-nucleotide polymorphisms (nsSNPs) for which evidence of association with breast cancer risk had been previously reported. Case-control data were combined from 38 studies of white European women (46 450 cases and 42 600 controls) and analyzed using unconditional logistic regression. Strong evidence of association was observed for three nsSNPs: ATXN7-K264R at 3p21 [rs1053338, per allele OR = 1.07, 95% confidence interval (CI) = 1.04-1.10, P = 2.9 × 10(-6)], AKAP9-M463I at 7q21 (rs6964587, OR = 1.05, 95% CI = 1.03-1.07, P = 1.7 × 10(-6)) and NEK10-L513S at 3p24 (rs10510592, OR = 1.10, 95% CI = 1.07-1.12, P = 5.1 × 10(-17)). The first two associations reached genome-wide statistical significance in a combined analysis of available data, including independent data from nine genome-wide association studies (GWASs): for ATXN7-K264R, OR = 1.07 (95% CI = 1.05-1.10, P = 1.0 × 10(-8)); for AKAP9-M463I, OR = 1.05 (95% CI = 1.04-1.07, P = 2.0 × 10(-10)). Further analysis of other common variants in these two regions suggested that intronic SNPs nearby are more strongly associated with disease risk. We have thus identified a novel susceptibility locus at 3p21, and confirmed previous suggestive evidence that rs6964587 at 7q21 is associated with risk. The third locus, rs10510592, is located in an established breast cancer susceptibility region; the association was substantially attenuated after adjustment for the known GWAS hit. Thus, each of the associated nsSNPs is likely to be a marker for another, non-coding, variant causally related to breast cancer risk. Further fine-mapping and functional studies are required to identify the underlying risk-modifying variants and the genes through which they act.

Genetic variation in mitotic regulatory pathway genes is associated with breast tumor grade

Mitotic index is an important component of histologic grade and has an etiologic role in breast tumorigenesis. Several small candidate gene studies have reported associations between variation in mitotic genes and breast cancer risk. We measured associations between 2156 single nucleotide polymorphisms (SNPs) from 194 mitotic genes and breast cancer risk, overall and by histologic grade, in the Breast Cancer Association Consortium (BCAC) iCOGS study (n = 39 067 cases; n = 42 106 controls). SNPs in TACC2 [rs17550038: odds ratio (OR) = 1.24, 95% confidence interval (CI) 1.16-1.33, P = 4.2 × 10(-10)) and EIF3H (rs799890: OR = 1.07, 95% CI 1.04-1.11, P = 8.7 × 10(-6)) were significantly associated with risk of low-grade breast cancer. The TACC2 signal was retained (rs17550038: OR = 1.15, 95% CI 1.07-1.23, P = 7.9 × 10(-5)) after adjustment for breast cancer risk SNPs in the nearby FGFR2 gene, suggesting that TACC2 is a novel, independent genome-wide significant genetic risk locus for low-grade breast cancer. While no SNPs were individually associated with high-grade disease, a pathway-level gene set analysis showed that variation across the 194 mitotic genes was associated with high-grade breast cancer risk (P = 2.1 × 10(-3)). These observations will provide insight into the contribution of mitotic defects to histological grade and the etiology of breast cancer.