Abstract
Accumulating evidence suggests a relationship between endometrial cancer and epithelial ovarian cancer. For example, endometrial cancer and epithelial ovarian cancer share epidemiological risk factors and molecular features observed across histotypes are held in common (e.g. serous, endometrioid and clear cell). Independent genome-wide association studies (GWAS) for endometrial cancer and epithelial ovarian cancer have identified 16 and 27 risk regions, respectively, four of which overlap between the two cancers. Using GWAS summary statistics, we explored the shared genetic etiology between endometrial cancer and epithelial ovarian cancer. Genetic correlation analysis using LD Score regression revealed significant genetic correlation between the two cancers (rG = 0.43, P = 2.66 × 10−5). To identify loci associated with the risk of both cancers, we implemented a pipeline of statistical genetic analyses (i.e. inverse-variance meta-analysis, co-localization, and M-values), and performed analyses by stratified by subtype. We found seven loci associated with risk for both cancers (PBonferroni < 2.4 × 10−9). In addition, four novel regions at 7p22.2, 7q22.1, 9p12 and 11q13.3 were identified at a sub-genome wide threshold (P < 5 × 10−7). Integration with promoter-associated HiChIP chromatin loops from immortalized endometrium and epithelial ovarian cell lines, and expression quantitative trait loci (eQTL) data highlighted candidate target genes for further investigation.
Introduction
Ovarian cancer is the eighth most commonly diagnosed cancer in women with 295,000 new cases annually1. Epithelial ovarian cancer accounts for ~90% of ovarian tumors and is commonly divided into five major histotypes: high-grade serous, low-grade serous, mucinous, clear cell and endometrioid2. Herein, “ovarian cancer” refers to epithelial types of this disease. On both histological and molecular levels, it is evident that ovarian cancer is a highly heterogeneous disease. Endometrial cancer (cancer of the uterine lining) is a comparatively understudied gynecological cancer, although it ranks fifth for cancer incidence in women globally, with 380,000 new cases diagnosed annually1. Endometrial cancer also has several histotypes, the most common being endometrioid (~80% of cases) but also includes serous, mucinous and clear cell.
Comparison of the epidemiology and histopathology of endometrial cancer and ovarian cancer has identified a number of similarities suggesting that shared molecular mechanisms underlie the pathology of these two diseases. Both cancers are hormone related, with epidemiological studies showing concordant direction of effect in relation to exposure to estrogen and progesterone (reviewed by Cramer 3). Protective factors for both types of cancer include early menopause4,5, late age of menarche6,7, longer periods of breastfeeding8,9, and longer use of contraceptives that include progesterone10,11 (i.e. factors that decrease exposure to unopposed estrogen). Although more strongly associated with endometrial cancer risk, higher body mass index (BMI) has been reported to be associated with increased risk of both cancers12,13.
The histotypes of endometrial cancer mirror those of ovarian cancer, albeit with varied frequencies observed across the two cancers. For example, serous histology is found in ~70% of ovarian tumors, compared with 10% of endometrial tumors, while endometrioid histology is found in ~10% of ovarian tumors and 80% of endometrial tumors. Clear cell and mucinous histologies are found in a relatively low frequency in both ovarian and endometrial tumors. Common features have been observed in similar histotypes regardless of the organ of origin. Tumors with serous histology from both the endometrium and ovary are characterized by somatic defects in the tumor suppressor gene, TP5314,15. Endometrioid endometrial and endometrioid ovarian tumors have both been found to contain somatic alterations in PTEN, PIK3CA, ARID1A, PPP2R1A and CTNNB1, although the frequencies of these mutations vary by tissue type (reviewed by McConechy, et al. 16). Methylation profiling has found endometrioid endometrial and endometrioid ovarian tumors cluster together17, and similar gene expression patterns have been observed for clear cell endometrial and clear cell ovarian tumors18. Further, there is increasing evidence that clear cell and endometrioid ovarian tumors arise in part from endometriosis (reviewed by King, et al. 19). Endometriosis is a chronic disease affecting reproductive aged women, in which endometrium grows outside of the uterus, suggesting these ovarian cancer subtypes and endometrial cancer develop from similar precursor endometrial epithelial cells.
Some, but not all germline cancer risk variants are also shared between endometrial cancer and ovarian cancer. Lynch Syndrome, characterized by germline pathogenic variants in the mismatch repair genes (i.e. MLH1, MSH2 and MSH6), is associated with 40-60% and 8-15% lifetime risks of endometrial cancer and ovarian cancer, respectively20. Additionally, separate genome-wide association studies (GWAS) of the two cancer types have identified four genetic risk regions common to both cancers21,22.
Meta-analyses of GWAS datasets across etiologically-related diseases have successfully been used to increase statistical power and identify novel genetic risk regions23,24. Hence, in the current study, we have performed a joint meta-analysis of the largest endometrial cancer and ovarian cancer GWAS datasets to identify novel genetic loci associated with risk of both cancers, including risk variation specific to less common ovarian cancer subtypes. To identify candidate target genes at such loci, we have intersected risk variation with chromatin looping data enriched for promoter-enhancer interactions. We have also assessed associations between risk variation and gene expression to provide evidence of candidate target gene regulation and reveal further candidate genes.
Methods
GWAS Datasets
GWAS summary statistics were obtained from the latest meta-analyses performed by the Endometrial Cancer Association Consortium (ECAC)21 and the Ovarian Cancer Association Consortium (OCAC)22. Because of the low number of non-endometrioid endometrial cancer available in ECAC, summary statistics were provided for all endometrial cancer risk (including all endometrial cancer cases) and analyses restricted to endometrioid cases only. OCAC summary statistics were available for all ovarian cancer risk (including all ovarian cancer cases), as well analyses restricted to eight different subtypes: endometrioid histology, serous (including borderline, high- and low-grade serous cases), serous high-grade histology, serous low-grade histology, serous borderline histology, serous low-grade and borderline cases combined, clear cell histology and mucinous histology. Sample sizes for each study and subgroups analyzed are provided in Table 1. Details on genotyping, quality control and imputation have been previously described21,22. Data for approximately 10 million genetic variants (imputation quality score > 0.4 and minor allele frequency > 0.01) were available for both cancers for the present study.
Genetic Correlation Analyses
Genetic correlation (i.e. the estimated proportion of variance shared between two traits due to genetic factors) between endometrial cancer and ovarian cancer was assessed using linkage disequilibrium (LD) Score Regression25. Genetic correlation was also assessed between each of the ovarian cancer subtypes analyzed by OCAC and all endometrial cancer as well as restricted to endometrioid endometrial cancer. For this analysis, the complete set of GWAS variants were pruned to the HapMap3 variant list (~1 million variants) to provide variants with high confidence imputation scores. The major histocompatibility complex (MHC) region was removed from this analysis because of its complex LD structure.
Cross-cancer GWAS meta-analyses
To identify joint endometrial and ovarian cancer genetic risk variants, summary statistics from ECAC and OCAC were combined by inverse-variance meta-analysis, adjusting for unknown sample overlap using MTAG26. Because of the significant heterogeneity in risk estimates observed for genetic variants across ovarian cancer subtypes22, we additionally performed meta-analysis combining results from ECAC (all endometrial cancer or endometrioid endometrial cancers) with summary statistics from each of the nine ovarian cancer subtypes analyzed by OCAC (listed in Table 1). To minimize false positives, following inverse-variance meta-analysis, output variants were restricted to those meeting the following criteria: (i) concordant direction of effect on risk of both cancers; (ii) no significant heterogeneity in risk estimates between the two cancers (Phet > 0.05); and (iii) associated with each cancer at nominal significance (P < 0.05). Counts of variants meeting these criteria are provided in Supplementary Table S1. M-values27 were generated for variants reaching with suggestive evidence of association (P < 5 x 10−7) using METASOFT28. Variants with a posterior probability for an effect in each study (M-value > 0.9) were retained for further consideration.
Loci containing variants that were statistically significant in the meta-analysis were further evaluated for co-localization by GWAS-PW29, using all genetic variants at the query locus. GWAS-PW estimates Bayes factors and posterior probabilities of association (PPA) for four models: (i) a locus associates with risk of endometrial cancer only; (ii) a locus associates with risk of ovarian cancer only; (iii) a locus contains a risk signal that associates with risk of both endometrial and ovarian cancers; or (iv) a locus contains two risk signals that associate independently with risk of either endometrial or ovarian cancer. Risk signals located in loci that were classified as meeting model (iii) were considered to be joint endometrial and ovarian cancer signals (PPA > 0.5).
Cell culture
IOSE11 (immortalized ovarian surface epithelial)30 cells were gifted from Prof S Gayther (Cedars-Sinai Medical Center). Cells were authenticated using STR profiling and confirmed to be negative for Mycoplasma contamination. For routine culture, IOSE11 were grown in 1:1 MCDB105:Medium 199 with 15% FBS and antibiotics (100 IU/ml penicillin and 100 μg/ml streptomycin).
Cell fixation
For fixation, cells (~80% confluent on 10 cm tissue culture plates) were washed with PBS and fixed at room temperature in 1% formaldehyde in PBS. After 10 min, the reaction wasquenched by washing with 125 mM glycine in PBS and then adding fresh glycine-PBS. Cells were removed from the dish with a cell scraper and washed with PBS before storing cell pellets at -80°C.
HiChIP library generation
HiChIP libraries were generated as previously31. Briefly, cell nuclei were extracted from fixed cell pellets and digested overnight with DpnII. After digestion, restriction fragment overhangs were filled in with biotin-dATP using the DNA polymerase I, large Klenow fragment. Proximity ligation was then performed, nuclei lysed and chromatin sheared. Sheared chromatin was incubated overnight with H3K27Ac antibody (Abcam, EP16602) to enrich for chromatin associated with promoters or enhancers. The next day Protein A beads were used to capture H3K27Ac-associated chromatin which was eluted and purified with Zymo Research concentrator columns. The DNA concentration of the purified chromatin was used to estimate the amount of TDE1 enzyme (Illumina) needed for tagmentation which was performed with biotin-labelled chromatin captured on streptavidin beads. Sequencing libraries were then generated using tagmented samples and the Nextera DNA preparation kit (Illumina). Size selection was performed using Ampure XP beads to capture 300–700 bp fragments. Two independent sequencing libraries were pooled together to provide 25 μl of library at ≥ 10 nM for Illumina HiSeq4000 (AGRF, Brisbane, QLD, Australia) paired-end sequencing with read lengths of 75 bp.
HiChIP bioinformatics analyses
HiChIP reads (fastq files) were aligned to the human reference genome (hg19) using HiC-Pro v2.9.032 and default settings were used to remove duplicate reads, assign reads to DpnII restriction fragments and filter for valid interactions as previously31 All valid reads from Hi-Pro were processed by the hichipper pipeline v0.7.033 as previously31. Chromatin interactions were filtered using a minimum distance of 5 kb and a maximum of 2 Mb. The final set of chromatin loops used for further investigation were interactions that were supported by a minimum of two unique paired end tags and with a Mango34 q-value < 5%. Promoter- associated chromatin loops were defined as HiChIP loops with anchors within ± 3 kb of a transcription start site. Promoter-associated chromatin looping data was also available from our previous analysis of a normal immortalized endometrial cell line (E6E7hTERT)31.
Credible candidate risk variants
Using 100:1 log likelihood ratios, “credible variants” (CVs) were identified at each of the joint endometrial and ovarian cancer risk regions. To identify genes that could be distally regulated by a CV, intersections of CVs with promoter-associated chromatin loops were performed using bedtools v2.28.0. Identification of genes whose expression is associated with a CV was performed by lookup of publicly available eQTL databases, including precomputed eQTL results from 336 endometrial and 318 ovarian tumors from the Cancer Genome Atlas (https://albertlab.shinyapps.io/tcga_eqtl)35, and from 101 non-cancerous uterus samples and 122 ovarian tissue samples from GTEx (data release v7; http://gtexportal.org)36. Additionally, due to the substantially increased power the sample size provided over solid tissue analyses, we accessed eQTL results from 31,684 whole blood samples (http://eqtlgen.org)37. Genes were considered potential targets if their expression associated with CVs that had a p-values within two orders of magnitude of the best eQTL variant in any of these eQTL datasets.
Significant genetic correlation was observed between all endometrial cancer and all ovarian cancer (rG = 0.43, P = 2.66 × 10−5; Table 2). When broken down by ovarian cancer subtype, we observed significant correlation between endometrial cancer and the following subgroups; endometrioid (rG = 0.53, P = 7.0 × 10−3), serous (rG = 0.42, P = 1.0 × 10−4) and high-grade serous ovarian cancers (rG = 0.44, P = 1.0 × 10−4). These correlations remained significant, although attenuated, when using endometrioid endometrial cancers only (Table 2).
Seven genetic loci displaying evidence of a joint association with risk of both endometrial cancer (all or endometrioid histology) and ovarian cancer (all or one of the subtypes) (i.e. PPA > 0.5 for GWAS-PW model iii), passed Bonferroni-correction for multiple testing (5 × 10−8/17 tests = 2.9 × 10−9; Table 3). Three of these loci belong to regions that have previously been reported as being associated with risk of both cancers (8q24, 17q12 and 17q21.32), although the 17q21.32 region had not been reported to be associated with the specific subtypes of ovarian cancer found in this meta-analysis (Table 3). One of the seven loci (2p16.1) has been previously reported as being associated with risk of endometrial cancer, but not with ovarian cancer risk. The three remaining loci (5p15.33, 9q34.2 and 10p12.31) have been previously reported as associated with risk of all ovarian cancer and serous ovarian cancer but not with endometrial cancer risk below GWAS significance levels; however, associations between endometrial cancer and variants in the 5p15.33 (TERT) region have been reported in a candidate-region study38. Additionally, we identified four novel loci with sub-GWAS significance levels (P < 5 × 10−7) that had not been previously reported as being associated with risk of either cancer at genome-wide levels of significance (7p22.2, 7q22.1, 9p12 and 11q13.3, Figure 1).
We identified a total of 22 candidate target genes at the 11 identified joint endometrial and ovarian cancer risk loci using a number of approaches (Table 4, Supplementary Table S2). Log likelihood ratios identified a median of 20 CVs per locus (range 1–73, Supplementary Table S3). Using H3K27Ac-associated chromatin looping data from normal immortalized ovarian surface epithelial cells and the same data previously generated from a normal immortalized endometrium cell line31, we intersected CVs coincident with putative enhancers (marked by H3K27Ac) belonging to promoter-associated loops. We found looping between such enhancers and the promoters of 14 genes (at five of the 11 loci) to be common to both immortalized endometrium and ovarian surface epithelial cell lines (e.g. Figure 1). These included genes which encode proteins involved in relevant processes such as steroid hormone metabolism (CYP3A43), estrogen response (GPER) and oncogenesis (MYC, CCDN1). Four of the 14 candidate target genes identified by chromatin looping also had CV located in the promoter, indicating potential to regulate expression (Table 4). An additional five genes were identified as candidate targets with CVs located in the corresponding promoters (Table 4). Interrogation of five relevant public eQTL databases revealed CVs to be associated with the expression of four genes (ABO, BCL11A, HOXB2 and SNX11), highlighting them as candidate targets. One of these, SNX11, had also been identified through the chromatin looping analyses and a CV was located in its promoter. Notably, we observed that increased expression of ABO associated with risk allele of CVs at the 9q34.2 locus in all five eQTL datasets: blood, non-cancerous uterine and ovarian tissues, and endometrial and ovarian tumors.
Discussion
In this study, we have performed the first cross-cancer GWAS analysis of endometrial cancer and ovarian cancer. Genetic correlation analyses found significant correlation between the two cancers, particularly between all endometrial cancer (and its endometrioid subtype) and the serous (high- and low-grade combined) or endometrioid ovarian cancer subtypes. Our pipeline of genetic analyses, stratifying by subtype, allowed us to identify seven joint endometrial cancer and ovarian cancer genetic risk loci. Three of these loci were located in regions that had been previously associated with both cancers, one was located in a known endometrial cancer risk region and the remaining three were located in known ovarian cancer risk regions. Four novel genetic risk loci for these two cancers did not reach the statistical threshold for significance but were highlighted as of potential interest, requiring further study to confirm their status.
Joint endometrial and ovarian cancer risk loci are located in the 8q24.21 and 5p15.33 regions, previously described as “cancer GWAS nexus regions”39 since genetic variation at these regions has been associated with many different types of cancer. 8q24.21 has been previously identified as a genetic risk region for both endometrial cancer and ovarian cancer21,22. CVs in a putative enhancer at the 8q24.21 joint endometrial and ovarian cancer risk locus showed evidence of chromatin looping to the promoter of the pan-cancer MYC oncogene in immortalized endometrial epithelial and ovarian surface epithelial cell lines. A previous study of the 5p15.33 multi-cancer risk region, containing the TERT gene, identified two independent signals for ovarian cancer risk: one (lead variant rs7705526) associated with serous borderline ovarian cancer risk and the other (lead variant rs10069690) associated with serous invasive ovarian cancer risk40. Although not previously associated with risk of endometrial cancer at genome-wide significance, a candidate fine-mapping study of 5p15.33 did highlight three independent endometrial cancer risk signals at this locus at study-wide significance38, one of which was shared with the serous borderline ovarian cancer risk signal. The present analysis identified this signal as a joint endometrial and ovarian cancer risk signal, with CVs in the TERT promoter highlighting this gene as a likely target. Moreover, the TERT protein has been heavily implicated in cancer development (reviewed in Yuan, et al. 41) and has oncogenic interactions with MYC (reviewed in Pestana, et al. 42).
Our results suggest, at a sub-genome wide significance level, a potential joint endometrial and ovarian cancer risk signal at another cancer GWAS nexus region, 11q13.3. Originally identified as a prostate cancer risk locus, 11q13.3 also contains risk signals for melanoma, breast cancer and renal cancer (https://www.ebi.ac.uk/gwas/). Although the results from the present study require validation, the identification of a shared endometrial and ovarian cancer risk signal at 11q13.3 would provide further evidence that this region is important for cancer development. At this locus, chromatin looping data showed that CVs in a putative enhancer looped to the promoters of MYEOV and CCND1, in immortalized endometrial epithelial and ovarian surface epithelial cell lines. CCND1 (encoding cyclin D1) is of particular interest as it is frequently amplified in human cancers and has been identified as a pan-cancer driver gene43. Cyclin D1 is considered an oncogene due to its central role in cell cycle regulation, and ability to promote cell proliferation44. CCND1 has been found to be significantly mutated in gynecological (endometrial, ovarian and cervical cancer, and uterine carcinosarcoma) and breast cancers45. The results of our genetic association analyses and integrative analyses of chromatin interactions results provide additional support that CCND1 is important in the development of endometrial cancer and ovarian cancer.
Our analysis identified the 17q12 region as a joint endometrial and ovarian cancer risk region, associating with clear cell ovarian cancer. The 17q12 region, containing HNF1B, has been previously associated with risk of endometrial cancer and ovarian cancer46–49. Significant heterogeneity in risk estimates has been observed across ovarian cancer histotypes at this locus. The minor allele of the lead ovarian cancer risk variant previously identified at this region associated with increased serous (high- and low-grade combined) ovarian cancer risk but decreased clear cell ovarian cancer risk48,49. Further genotyping had resolved this region into two risk signals for ovarian cancer risk: one in intron 1 of HNF1B for clear cell ovarian cancer risk (rs11651775; the same signal for endometrial cancer risk) and another in intron 3 for serous ovarian cancer risk (rs7405776)49. Our results confirm that joint endometrial and ovarian cancer risk variants at 17q12 map to the same signal as that for that previously reported for endometrial cancer and the clear cell ovarian subtype. HNF1B is a likely target of endometrial and ovarian cancer risk variation, with CVs located in the promoter region of this gene. We have previously demonstrated that these variants affect activity of the HNF1B promoter46, which may lead to increased secretion of insulin, a risk factor for endometrial cancer50.
The 17q21.32 region is a known shared endometrial21 and ovarian cancer22 risk region. The joint endometrial and ovarian cancer signal found in the present study (lead SNP rs882380) is the same as that previously identified for endometrial cancer, but is independent of the signal previously found for all invasive and high-grade serous ovarian cancer risk (lead SNP rs7207826, r2 = 0.06 with rs882380). The joint endometrial and ovarian cancer signal associates specifically with risk of clear cell, endometrioid, serous low-grade, serous low-grade and borderline combined, and serous borderline ovarian cancer subtypes. Clear cell, endometrioid and serous low-grade ovarian cancers are often referred to as endometriosis-associated ovarian cancers due to the increased risk of these ovarian cancer subtypes with endometriosis51. Epidemiological and molecular data provide strong evidence that clear cell and endometrioid ovarian cancer arise in part from endometriosis (reviewed by King, et al. 19). The joint endometrial and ovarian cancer signal identified in the present study at 17q21.32 was also found in a joint GWAS analysis of endometrial cancer and endometriosis52, and subsequently found to be associated with endometriosis risk independently53. Five candidate target genes were identified at this locus, all of which we had previously found to be candidate targets of the original endometrial cancer signal through chromatin looping studies of endometrial cancer cell lines31.
Another potential joint endometrial and ovarian cancer signal, 9p12, associated with risk of serous low-grade ovarian cancer, has also been previously identified as a joint endometrial cancer and endometriosis risk locus52. These findings at 17q21.32 and 9p12, add to the body of evidence for the relationship between endometriosis and specific ovarian cancer subtypes19,51, and provide further support for shared genetic etiology between endometriosis and endometrial cancer52. CVs at the 9p12 joint risk locus were located intronic to PTPRD, but no candidate target genes were identified through chromatin looping or eQTL analyses. PTPRD protein is involved in the STAT3 pathway which has been implicated as a potential target for both endometrial cancer54 and ovarian cancer55.
The 2p16.1 region is a known endometrial cancer risk locus and was found to associate with the risk of clear cell ovarian cancer only. Interestingly, we previously found evidence that this locus may have a stronger association with risk of non-endometrioid endometrial cancer, with the strongest effect observed for clear cell endometrial cancer subtype (128 cases and 26,638 controls; rs148261157 OR 2.36; 95% CI 1.07 - 5.19)21. BCL11A was identified as a candidate target gene through eQTL analysis of endometrial tumors. We had previously found that BCL11A was a candidate target gene at the endometrial cancer risk locus through chromatin looping studies in endometrial cancer cells31. The eQTL finding suggested that reduced expression of BCL11A may increase endometrial/clear cell ovarian cancer risk. Indeed, some studies have shown that BCL11A acts as a proto-oncogene56,57; however, others suggest that overexpression of BCL11A results in anti-cancer effects58. Notably, BCL11A has been found to be mutated in clear cell ovarian cancer59,60, providing further evidence that the expression of BCL11A explains, at least in part, the mechanism underlying the risk association with both endometrial cancer and clear cell ovarian cancer.
The 9q34.2 region is a known ovarian cancer risk locus that is highly pleiotropic, having been previously associated with gastric and pancreatic cancers, in addition to a wide range of traits including blood cell counts, the tumor marker CEA (carcinoembryonic antigen), circulating cholesterol, bone mineral density and levels of proteins related to angiogenesis (e.g. VEGFR-2 and angiopoietin)(https://www.ebi.ac.uk/gwas/). eQTL data from normal, tumor endometrial and ovarian tissue, as well as blood, provide evidence that ABO is a regulatory target of CVs at this locus. ABO encodes an enzyme with glycosyltransferase activity and determines human ABO blood group antigens. It is not immediately apparent how ABO may mediate cancer risk but its encoded glycosyltransferase can affect cell recognition and adhesion, and activation of T and natural killer cells (reviewed by Arend 61).
The 10p12.31 region is another known ovarian cancer risk locus that is also pleiotropic, having been previously associated with breast cancer as well as with traits related to obesity such as BMI, body fat percentage and physical activity (https://www.ebi.ac.uk/gwas/). MLLT10 was identified as a candidate target gene at this locus, through chromatin looping analysis and localization of a CV to its promoter, and has been found to be a partner gene for chromosomal rearrangements that result in leukaemia62. Another biologically relevant candidate target gene at this locus is MIR1915 whose expression is upregulated by p53 in response to DNA damage, subsequently leading to increased apoptosis63.
Two of the sub-genome wide significant endometrial/ovarian cancer risk regions (7q22.1 and 7p22.2) may relate to circulating hormone levels or regulation. At 7q22.1, GWAS have previously revealed associations with androgen and progesterone levels64. The sole candidate target gene at this locus, CYP3A43, encodes a cytochrome P450 enzyme that may be involved in androgen metabolism65 and is upregulated in ovarian tumors66. At 7p22.2, the candidate target gene GPER1, identified through chromatin looping, encodes an estrogen receptor that induces endometrial and ovarian cancer cell proliferation in response to estrogen (reviewed in Prossnitz and Barton 67). Further, it appears that androgen can also bind to GPER1 to stimulate cancer cell growth68.
Despite these findings, the present study does have some limitations. The low numbers of non-endometrioid endometrial cancers meant we could not explore the relationship of these endometrioid histotypes and their relationship with ovarian cancer. Another limitation was the use of cell lines to model chromatin looping that occurs in endometrial and ovarian tissue, with chromatin looping potentially impacted by the immortalization and 2D-culturing processes of cell lines, or mutations gained through routine passaging of cells. Only one endometrial and one ovarian cell line were used to identify chromatin looping events. These experiments should be repeated in additional endometrial and ovarian cell lines, representing tumor subtypes. One of the four regions previously identified to be associated with both cancers, located at 1p34, was not identified in the present analysis. This locus was originally found in a combined analysis of the OCAC with a cohort of BRCA1/2 carriers with ovarian cancer, i.e. the Consortium of Investigators of Modifiers of BRCA1 and BRCA2 (CIMBA)69. The CIMBA study was not included in the present study, perhaps explaining why it was not identified in our analysis as a joint endometrial and ovarian cancer locus. Future analysis of this region, in the context of BRCA1/2 carrier status will be required to explore how this region affects endometrial cancer and ovarian cancer risk.
In summary, using endometrial and ovarian cancer GWAS summary statistics we have been able to identify seven joint risk loci for these cancers, with an additional four novel potential risk regions at a sub-GWAS significance level. Further studies are required to validate these findings in larger sample sets. Notably, we also found significant genetic correlation between the two cancers, supported by the observed epidemiological and histopathological similarities. These findings support the need for larger GWAS of endometrial and ovarian cancer, in particular focusing on their minor subtypes to further explore shared genetic etiology. Integration of CVs with chromatin looping and eQTL data has identified several plausible candidate target genes, including those at potentially novel risk loci. Although the role of these genes in endometrial and ovarian cancer development should be explored in future studies, the current findings provide insights into the shared biology of endometrial and ovarian cancer.
Data Availability
OCAC and ECAC summary statistics are available from the GWAS catalog. The E6E7hTERT HiChIP data are available from GEO and IOSE11 HiChIP data will be made available upon publication.
Acknowledgements
We thank ECAC and OCAC for provision of summary statistics to perform this study. We thank Siddhartha Kar for his helpful discussions and advice for designing the genetic analysis approaches. Full acknowledgements and funding for ECAC and OCAC can be found in the Supplementary Note.
Footnotes
Financial Support: TAO’M is supported by a National Health and Medical Research Council (NHMRC) Early Career Fellowship (APP1111246) and Investigator Fellowship (APP1173170), ABS is supported by an NHMRC Senior Research Fellowship (APP1061779) and Investigator Fellowship (APP1177524). This work was supported by a Cancer Australia PdCCRS Project Grant, funded by Cure Cancer Australia and the CanToo Foundation (#1138084), NHMRC Project Grants (APP1158083 and APP1109286), QIMR Berghofer Medical Research Institute Near Miss Funding and a special purpose donation from gratefully received from Sarah Stork.
Conflict of Interest Statement: UM has stocks in Abcodia awarded to her by UCL. All other authors declare no potential conflicts of interest.