ORIGINAL ARTICLE: GENETICS
Evidence of a genetic link between endometriosis and ovarian cancer Alice W. Lee, M.P.H.,a Claire Templeman, M.D.,b Douglas A. Stram, B.S.,a Jonathan Beesley, Ph.D.,c Jonathan Tyrer, Ph.D.,d Andrew Berchuck, M.D.,e Paul P. Pharoah, Ph.D.,d,f Georgia Chenevix-Trench, Ph.D.,c and Celeste Leigh Pearce, Ph.D.,a,g on behalf of the Ovarian Cancer Association Consortium a Department of Preventive Medicine, Keck School of Medicine, University of Southern California, Los Angeles, California; b Department of Obstetrics and Gynecology, Keck School of Medicine, University of Southern California, Los Angeles, California; c Cancer Division, QIMR Berghofer Medical Research Institute, Brisbane, Queensland, Australia; d Department of Oncology and Strangeways Research Laboratory, University of Cambridge, Cambridge, United Kingdom; e Department of Obstetrics and Gynecology, Duke University Medical Center, Durham, North Carolina; f Department of Public Health and Primary Care and Strangeways Research Laboratory, University of Cambridge, Cambridge, United Kingdom; and g Department of Epidemiology, School of Public Health, University of Michigan, Ann Arbor, Michigan
Objective: To evaluate whether endometriosis-associated genetic variation affects risk of ovarian cancer. Design: Pooled genetic analysis. Setting: University hospital. Patient(s): Genetic data from 46,176 participants (15,361 ovarian cancer cases and 30,815 controls) from 41 ovarian cancer studies. Intervention(s): None. Main Outcome Measure(s): Endometriosis-associated genetic variation and ovarian cancer. Result(s): There was significant evidence of an association between endometriosis-related genetic variation and ovarian cancer risk, especially for the high-grade serous and clear cell histotypes. Overall we observed 15 significant burden statistics, which was three times more than expected.
Received June 8, 2015; revised August 17, 2015; accepted September 16, 2015. P.P.P. reports grants from the National Institutes of Health and Cancer Research UK during the conduct of the study. A.W.L. has nothing to disclose. C.T. has nothing to disclose. D.A.S. has nothing to disclose. J.B. has nothing to disclose. J.T. has nothing to disclose. A.B. has nothing to disclose. G.C.-T. has nothing to disclose. C.L.P. has nothing to disclose. This work was supported by the National Cancer Institute (K07-CA095666, K07-CA80668, K07-CA143047, K22-CA138563, N01-CN55424, N01-PC67001, N01PC067010, N01-PC035137, P01-CA017054, P01-CA087696, P30-CA072720, P30-CA15083, P30-CA008748, P50-CA159981, P50-CA105009, P50-CA136393, R01-CA149429, R01-CA014089, R01-CA016056, R01-CA017054, R01-CA049449, R01-CA050385, R01-CA054419, R01-CA058598, R01-CA058860, R01CA061107, R01-CA061132, R01-CA063678, R01-CA063682, R01-CA067262, R01-CA071766, R01-CA074850, R01-CA076016, R01-CA080978, R01CA083918, R01-CA087538, R01-CA092044, R01-CA095023, R01-CA122443, R01-CA112523, R01-CA114343, R01-CA126841, R01-CA136924, R03CA113148, R03-CA115195, U01-CA069417, U01-CA071966, and Intramural research funds); the European Commission's Seventh Framework Programme (agreement number 223175 HEALTH F2 2009-223175); a National Cancer Institute, Cancer Post-GWAS Initiative (U19-CA148112); the Genetic Associations and Mechanisms in Oncology (GAME-ON); the Ovarian Cancer Research Fund (thanks to donations by the family and friends of Kathryn Sladek Smith) (PPD/RPCI.07); the National Institute of Environmental Health Sciences (T32ES013678 to A.W.L.); National Health and Medical Research Council (to G.-C.T.); American Cancer Society Early Detection Professorship (SIOP-06-258-01-COUN to B.K.); National Center for Advancing Translational Sciences (NCATS) (UL1TR000124 to B.K.); California Cancer Research Program (00-01389V-20170, N01-CN25403, 2II0200); the Canadian Institutes of Health Research (MOP-86727); Cancer Australia; Cancer Council Victoria; Cancer Council Queensland; Cancer Council New South Wales; Cancer Council South Australia; Cancer Council Tasmania; Cancer Foundation of Western Australia; the Cancer Institute of New Jersey; Cancer Research UK (C490/A6187, C490/A10119, C490/A10124); the Danish Cancer Society (94-222-52); the ELAN Program of the University of Erlangen-Nuremberg; the Eve Appeal; the Helsinki University Central Hospital Research Fund; Helse Vest; the Norwegian Cancer Society; the Norwegian Research Council; the Ovarian Cancer Research Fund; Nationaal Kankerplan of Belgium; the L & S Milken Foundation; the Polish Ministry of Science and Higher Education (4 PO5C 028 14, 2 PO5A 068 27); the Roswell Park Cancer Institute Alliance Foundation; the National Institutes of Health/National Center for Research Resources/General Clinical Research Center (MO1-RR000056); the US Army Medical Research and Material Command (DAMD17-01-1-0729, DAMD1702-1-0666, DAMD17-02-1-0669, W81XWH-07-0449, W81XWH-10-1-02802); the U.S. Public Health Service (PSA-042205); the National Health and Medical Research Council of Australia (199600, 400281); the German Federal Ministry of Education and Research of Germany Programme of Clinical Biomedical Research (01 GB 9401); the State of Baden-Wurttemberg through Medical Faculty of the University of Ulm (P.685); the German Cancer Research Center; the Minnesota Ovarian Cancer Alliance; the Mayo Foundation; the Fred C. and Katherine B. Andersen Foundation; the Lon V. Smith Foundation (LVS-39420); the Oak Foundation; the Oregon Health & Science University Foundation; the Mermaid I project; the Rudolf-Bartling Foundation; the UK National Institute for Health Research Biomedical Research Centres at the University of Cambridge, Imperial College London, University College Hospital ‘‘Womens Health Theme’’ and the Royal Marsden Hospital; and WorkSafeBC 14. A.W.L., C.T., G.C.-T., and C.L.P. should be considered similar in author order. Reprint requests: Celeste Leigh Pearce, Ph.D., Department of Epidemiology, University of Michigan School of Public Health, 1415 Washington Heights, SPH Tower, Office #4642, Ann Arbor, MI 48109-2029 (E-mail:
[email protected]). Fertility and Sterility® Vol. -, No. -, - 2015 0015-0282/$36.00 Copyright ©2015 American Society for Reproductive Medicine, Published by Elsevier Inc. http://dx.doi.org/10.1016/j.fertnstert.2015.09.023 VOL. - NO. - / - 2015
1
FLA 5.4.0 DTD FNS29891_proof 21 October 2015 9:02 am ce JS
ORIGINAL ARTICLE: GENETICS Conclusion(s): By focusing on candidate regions from a phenotype associated with ovarian cancer, we have shown a clear genetic link between endometriosis and ovarian cancer that warrants further follow-up. The functional significance of the identified regions and SNPs is presently uncertain, though future fine mapping and histotype-specific functional analyses may shed light on the etiologies of both gynecologic conditions. (Fertil SterilÒ 2015;-: -–-. Ó2015 by American Society for Reproductive Medicine.) Key Words: Endometriosis, ovarian cancer, genetic variation, SNPs Discuss: You can discuss this article with its authors and with other ASRM members at http:// fertstertforum.com/leea-endometriosis-ovarian-cancer/
O
varian carcinoma (ovarian cancer) is the most fatal malignancy in the female reproductive system, accounting for more than 140,000 deaths annually worldwide (1). Endometriosis, the presence of ectopic endometrial glands and tissue in the peritoneum, is a common gynecologic condition, occurring in 6%–10% of the general female population (2). Studies that have controlled for parity have shown that endometriosis is a well-established ovarian cancer risk factor, especially for the endometrioid and clear cell histotypes (3). Although the etiology of endometriosis remains enigmatic, it is influenced by genetic factors, with an estimated heritability of 51% and an incidence that is approximately seven times higher in relatives of women with endometriosis than in women without such family history (4, 5). To date, seven variants reaching genome-wide significance have been identified in association with risk of endometriosis (6–8). In addition, the most recent meta-analysis by Nyholt et al. (9) has identified multiple additional variants associated with risk, albeit some at a sub-genome-wide significance level (P%1 105). Genetics also plays a role in the etiology of ovarian cancer: women with first-degree family histories of the disease have more than a twofold increased risk (10). High-penetrance susceptibility genes, such as BRCA1 and BRCA2, as well as 18 published common variants identified through genome-wide association studies (GWASs) account for a substantial portion of ovarian cancer's familial risk, but at least 60% remains unexplained (6–8,11–16). Lee et al. (17) have previously shown that the cumulative effect of many risk variants contributes to disease heritability, and hence it is likely that the germline genetic contributions to ovarian cancer are not limited to GWAS-identified variants. Given the important role genetics plays in the etiologies of both endometriosis and ovarian cancer and the consistent epidemiologic evidence of their association with one another, the two gynecologic conditions may also have a similar genetic profile. Hence, on the basis of the results presented by Nyholt et al. of their endometriosis GWAS meta-analysis, we present the first report that evaluates whether variation in the 18 regions harboring the top 38 endometriosis-associated singlenucleotide polymorphisms (SNPs) is associated with risk of ovarian cancer.
Use your smartphone to scan this QR code and connect to the discussion forum for this article now.* * Download a free QR code scanner by searching for “QR scanner” in your smartphone’s app store or app marketplace.
MATERIALS AND METHODS All studies included in this report obtained institutional ethics committee approval, and all participating subjects provided written, informed consent.
Study Populations Our analysis included 41 studies participating in the Ovarian Cancer Association Consortium (OCAC), an international collaboration of ovarian cancer studies founded in 2005. In total, 20 studies were conducted in Europe, 19 in North America, and 2 in Australia. Only participants of European ancestry were included, as was determined from the program LAMP (Local Ancestry in Admixed Populations) (see ‘‘Statistical Analysis,’’ below). We used a combined total of 46,176 participants (15,361 ovarian cancer cases and 30,815 controls) in our analyses; borderline tumors were excluded. Details regarding sample quality control have been previously published (15). Supplemental Table 1 (available online) provides an overview of each study's characteristics and the numbers of subjects included.
Genotyping and Imputation Analyses The genetic data for our analyses came from three population-based ovarian cancer GWASs, which comprised 2,162 cases and 2,564 controls from a GWAS in North America (‘‘US GWAS’’) (18), 1,763 cases and 6,118 controls from a GWAS in the United Kingdom (‘‘UK GWAS’’) (11), and 443 cases and 441 controls from a second GWAS in North America. In addition, 11,030 cases and 21,693 controls were genotyped using the iCOGS array, which was a large-scale genotyping project by the Collaborative Oncological Geneenvironment Study (COGS) (15). The US and UK GWASs included several independent case–control studies, and samples from these studies were also genotyped using the iCOGS array. All duplicates were removed from the analyses, resulting in genetic data for 15,361 women diagnosed with ovarian cancer and 30,815 controls. Study sets were created on the basis of the scope of genotyping information (GWAS vs. COGS) available for imputation. These sets are indicated in Supplemental Table 1. Details regarding the genotyping platform for each dataset have been published previously (15). Imputation of the entire scope of genetic variation in the genome was carried out separately for iCOGS samples and VOL. - NO. - / - 2015
2
FLA 5.4.0 DTD FNS29891_proof 21 October 2015 9:02 am ce JS
Fertility and Sterility® each of the GWASs. We imputed variants by combining all available genotype data with information from the April 2012 release of the 1000 Genomes Project using the program IMPUTE2 (19). In addition, all data were pre-phased using the software SHAPEIT, to improve computation efficiency (20).
SNP Selection The SNPs evaluated here were based on the results of Nyholt et al.'s genome-wide association meta-analysis (9), which included two large endometriosis GWASs: one conducted in a Japanese sample obtained from BioBank Japan (21) and the other in a European sample from Australia and the UK by the International Endogene Consortium (22). Nyholt et al. found six SNPs (rs7521902, rs12700667, rs10859871, rs4141819, rs1537377, and rs7739264) that reached genomewide significance (P%5 108) in their analysis, as well as an additional SNP (rs13394619) that reached genome-wide significance when combined with the results from a previous meta-analysis of two Japanese case–control cohorts (23). These seven SNPs, as well as an additional 31 SNPs that showed suggestive evidence of an association with endometriosis (P%1 105), were included in our analyses for a total of 38 index SNPs. A list of these SNPs is provided in Table 1. In addition, we evaluated genetic variation in the surrounding regions of each of the 38 index SNPs. Initially these regions were defined as the areas 25 kb up- and downstream of each index SNP. All variants in these 50-kb regions were analyzed, and the SNP with the strongest ovarian cancer association in each region was identified. The final regions, which we have labeled as regions A to R, were defined as including all SNPs with an r2 R 0.2 with the most significantly associated ovarian cancer SNPs. On the basis of each SNP's linkage disequilibrium pattern, a total of 18 regions of varying sizes was identified because some of the 38 SNPs were strongly correlated with each other. Only SNPs with an imputation r2 R 0.5 and a minor allele frequency (MAF) R0.05 were considered. In total, 6,981 SNPs (including the original 38 SNPs) were assessed.
Statistical Analysis We used LAMP to assign intercontinental ancestry according to genotype frequencies in European populations according to HapMap (24). Subjects were classified as European if they had 90% or more European ancestry. Principal components analysis to control for population substructure was also performed using a set of 37,000 unlinked markers, as well as an in-house program written in Cþþ that used the Intel MKL library for eigenvectors (http://ccge.medschl.cam.ac.uk/software/pccalc). We carried out unconditional logistic regression analyses, adjusting for the first five eigenvalues from the principal components analysis for European ancestry, to determine the association between each SNP and risk of ovarian cancer. For all ovarian cancer cases combined and for high-grade serous cases, the analyses were carried out within each study set, and the set-specific results were summarized using a fixed-effects meta-analysis approach. For
the less-common histotypes (mucinous, endometrioid, clear cell), the cases and controls were pooled across all study sets, and study set was included as a term in the model. A logadditive mode of inheritance was used, with each SNP modeled as an ordinal variable. Hence, the effect estimates reflect per-allele odds ratios (ORs). All P values reported herein are two-sided. To evaluate whether overall genetic variation implicated in risk of endometriosis also plays a role in risk of ovarian cancer, we calculated burden statistics using the admixture likelihood (AML) method (25). We used the AMLcalc program to accomplish this using 1,000 simulations with the maximum proportion of associated SNPs set to 0.2 on the genotyped and imputed data and adjusting for the first five ancestry principal components. P values for the AML trend test are provided. We calculated burden statistics across the 38 SNPs identified by Nyholt et al., as well as across each of the 18 regions for ovarian cancer overall risk and for the four main histotypes. This allowed us to take a global approach to assess whether the combined variation that exists among the 38 endometriosis-associated SNPs or across the 18 regions plays a role in ovarian cancer risk after accounting for the correlation between SNPs. Manhattan plots were generated by plotting the -log P vs. the chromosomal position using GraphPad Prism 6 to provide a picture of the distribution of P values by histotype across each region. This plotting was done for all regions except region A because it was presented in a recent meta-analysis by OCAC and the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA) (16). Linkage disequilibrium plots depicting pairwise correlation data from the 1000 Genomes CEU (Utah residents [CEPH] with Northern and Western European Ancestry) population were generated using HaploView. Epigenomic data made available from the ENCODE and Roadmap Epigenomics Consortia were obtained and visualized in the University of California, Santa Cruz Genome Browser.
RESULTS In total, 38 SNPs with a P%1 105 spanning 18 unique, uncorrelated regions were identified from the endometriosis GWAS meta-analysis carried out by Nyholt et al. (9) (Table 1). Table 1 presents the association between each of those 38 SNPs and risk of ovarian cancer. Eight SNPs (rs7515106, rs7521902, rs742356, rs4858692, rs1603995, rs4241991, rs6907340, and rs10777670) from five different regions (regions A, F, I, J, and R) showed statistically significant associations with ovarian cancer risk (P%.05), with the most significantly associated SNP being rs7515106 from region A on chromosome 1 (P¼4.9 106). When considering all 38 SNPs together, the calculated AML burden statistic showed significant evidence of an association with ovarian cancer risk (P¼ .001); it slightly attenuated when the SNPs from region A were excluded (P¼ .052). Given that each region was defined by the area containing SNPs with an r2 R 0.2 with the SNP most significantly associated with ovarian cancer risk, Table 2 presents the coordinates as well as the number of SNPs with genotyped or imputed data available for each region. Overall, a total of
VOL. - NO. - / - 2015
3
FLA 5.4.0 DTD FNS29891_proof 21 October 2015 9:02 am ce JS
ORIGINAL ARTICLE: GENETICS
TABLE 1 Association between the top 38 SNPs identified through endometriosis genome-wide scans and ovarian cancer risk. Endometriosisc a
Region A A A B C C D E F F G H I J J J J J J J K L M N N N N O O P Q Q Q Q R R R R
Chromosome 1p36 1p36 1p36 1p31 2p25 2p25 2p14 3p24 3p24 3p24 3p14 3q13 4q12 6p22 6p22 6p22 6p22 6p22 6p22 6p22 6q21 6q25 6q25 7p15 7p15 7p15 7p15 7q33 7q33 9p21 10q26 10q26 10q26 10q26 12q22 12q22 12q22 12q22
b
Position
22473410 22490724 22501846 80278733 11725241 11727507 67864675 24927203 25075167 25075530 55216779 104575864 55931246 19706761 19708481 19729003 19761215 19785588 19790809 19803768 113111808 152638209 158144402 25860812 25871109 25873110 25901639 134607676 134618710 22169700 117325021 117393524 117396288 117398195 95574831 95631276 95702385 95711876
SNP
OR e
Ovarian cancerd
P value 6
2.1 10 4.6 108 3.1 106 9.6 106 8.1 106 6.1 108 4.0 107 9.6 106 1.3 106 1.6 106 7.4 106 8.6 106 4.6 106 6.3 106 4.3 106 3.1 107 3.4 106 1.3 107 7.3 107 5.9 107 8.2 106 8.1 106 4.4 106 2.1 106 5.6 107 7.4 107 9.3 1010 2.3 106 4.7 106 2.5 106 1.5 106 8.2 106 7.6 106 6.9 106 1.6 107 4.8 107 2.1 106 5.5 109 Burden statistic (all SNPs): Burden statistic (excluding Region A SNPs):
rs7515106 rs7521902e rs742356 rs12744944e rs12470971 rs13394619 rs4141819e rs10510551 rs4858692 rs1603995 rs9311552 rs12635480 rs4241991 rs1016251 rs9366312 rs9356708e rs6916251 rs7739264e rs2223361 rs6907340 rs9487982 rs1890100 rs16900375 rs7798431e rs1055144e rs10282436e rs12700667e rs7809057 rs6973420 rs1537377e rs1572396 rs2769422e rs2804250 rs2769417 rs10777670e rs10859856 rs11107973 rs10859871e
1.18 1.18 1.15 1.15 1.14 1.15 1.16 1.14 1.16 1.16 1.13 1.34 1.16 1.14 1.14 1.16 1.14 1.16 1.15 1.15 1.22 1.16 1.29 1.16 1.18 1.18 1.22 1.15 1.14 1.14 1.16 1.14 1.15 1.15 1.19 1.15 1.14 1.18
OR
P value
1.09 1.08 1.05 1.00 0.97 0.98 1.00 1.01 0.96 0.96 1.00 0.99 0.95 1.02 1.02 1.02 1.02 0.98 1.03 1.03 0.99 1.01 1.01 0.98 0.97 0.97 0.97 1.00 1.00 1.03 1.02 1.01 1.01 1.01 1.05 1.01 1.01 1.03
4.9 106* 7.3 106* .004* .80 .062 .23 .86 .54 .025* .017 .85 .78 .018* .154 .14 .132 .13 .13 .093 .045* .68 .70 .84 .28 .19 .18 .072 .81 .87 .071 .24 .55 .54 .55 .021* .39 .36 .088 .001* .052
a
The 38 SNPs represent 18 distinct regions of linkage disequilibrium (Regions A-R). Based on build 37. Odds ratios and P values reported in Nyholt et al. (9) for the association between the SNP and risk of endometriosis. d Odds ratios and P-values for the association between the SNP and risk of ovarian cancer. e SNP was genotyped. * Significant P values at P%.05 for ovarian cancer. b c
Lee. Endometriosis SNPs and ovarian cancer. Fertil Steril 2015.
6,981 SNPs across 18 regions were evaluated. Region O was the largest, spanning almost 500 kb, but region R included the largest number of SNPs with an MAF R0.05 and an imputation r2 R 0.8 (n ¼ 1,309). Table 2 also presents the burden statistics for each region for all ovarian cancers combined and for the four main histotypes. In total there were 15 significant burden statistics at a P%.05 level, covering eight different regions. Of these, the most significant burden statistic was for region A (P¼5.0 105 for overall). However, there was a total of six regions (A, C, D, I, P, and R) that showed a significant association with ovarian cancer overall, and among them, five regions also had significant histotype-specific burden statistics. High-grade serous and clear cell cancers had the greatest
number of significant burden statistics, each with four and both sharing region A. When burden statistics were calculated globally across all 18 regions, significant evidence of an association was found between overall genetic variation and risk of all ovarian cancers combined (P¼ .002), as well as risk of the high-grade serous and clear cell histotypes (P¼ .002 and P¼ .039, respectively). Odds ratios and 95% confidence intervals (CIs) for the most significant SNPs associated with ovarian cancer are presented in Table 3. Region A contained the most significant SNP, rs10917151 (OR 1.11, P¼4.0 107), but regions L, P, and R also had SNPs of notable significance (P¼3.3 104 for rs73007780, P¼2.5 104 for rs1333052, P¼1.2 105 for rs7397212, respectively). Supplemental Tables 2 VOL. - NO. - / - 2015
4
FLA 5.4.0 DTD FNS29891_proof 21 October 2015 9:02 am ce JS
Fertility and Sterility®
TABLE 2 Burden statistics for the association between the 18 genomic regions, harboring SNPs associated with endometriosis, and risk of ovarian cancer by histotype. Burden statisticsc Region
Chr.
A 1p36 B 1p31 C 2p25 D 2p14 E 3p24 F 3p24 G 3p14 H 3q13 I 4q12 J 6p22 K 6q21 L 6q25 M 6q25 N 7p15 O 7q33 P 9p21 Q 10q26 R 12q22 Global –
Region coordinatesa 22347396–22531542 80152740–80338250 11661140–11712075 67983706–68019054 24940850–25060867 25035804–25075645 55181646–55221542 104584856–104828081 55988082–56027573 19567942–19632693 112966971–113226112 152544371–152677815 158184000–158192932 25888067–25897388 134455276–134936517 22140648–22201586 117345263–117407423 95304634–95709315 –
Region No. of SNPs size (kb) Coverage in regionb 184.146 185.510 50.935 35.348 120.017 39.841 39.896 243.225 39.491 64.751 259.141 133.444 8.932 9.321 481.241 60.938 62.160 404.681 –
0.99 0.99 0.98 0.97 0.96 0.96 0.94 0.99 0.96 0.97 0.97 1.00 0.94 0.95 0.98 0.96 0.99 0.98 –
Overall
High-grade serous
5.0 105* 1.0 104* .71 .057 .041* .70 .044* .21 .67 .85 .15 .32 .30 .32 .84 .95 .011* .010* .61 .54 .15 .63 .26 .39 .23 .50 .16 .38 .54 .96 .002* .026* .82 .56 .011* .025* .002* .002*
594 607 107 77 390 168 113 846 135 194 550 370 19 15 1,082 184 204 1,309 –
Mucinous Endometrioid Clear cell .78 .26 .11 .16 .28 .17 .79 .86 .23 .20 .19 .75 .075 .79 .70 .63 .13 .53 .61
.35 .93 .23 .43 .57 .036* .54 .34 .19 .23 .67 .26 .48 .22 .12 .17 .70 .15 .20
.007* .015* .014* .54 .65 .28 .68 .48 .018* .29 .82 .070 .61 .12 .94 .14 .27 .28 .039*
Note: ‘‘Global’’ refers to when all SNPs across all regions are considered. Chr. ¼ chromosome. a Defined as including all SNPs with an r2 R 0.2 with the most significant ovarian cancer-associated SNP (when considering the areas 25 kb up- and downstream of the index SNP as presented in Table 1). b Based on r2 R 0.5 and MAF R 0.05. c Took into account the first five ancestry principal components; P values for trend reported. * Significant P values at P%.05. Lee. Endometriosis SNPs and ovarian cancer. Fertil Steril 2015.
and 3 present the most significant SNPs for each region by histotype (high-grade serous and mucinous in Supplemental Table 2; endometrioid and clear cell in Supplemental Table 3). Similar to the results for ovarian cancer overall,
region A for high-grade serous cancer contained the most significant SNP, rs3754496 (P¼5.3 108), whereas region R for high-grade serous cancer harbored the largest number of significant SNPs (n ¼ 327) when looking across the four
TABLE 3 Odds ratios and 95% CIs for the most significant SNP in each of the 18 genomic regions for ovarian cancer overall. Region
No. of significant SNPs
Most significant SNP
ORa
95% CI
P value
192 2 53 14 5 22 10 6 87 1 59 24 2 1 0 68 2 256
rs10917151 rs28556754 rs142034631 rs17034382 rs62228500 rs1846555 rs17054879 rs114063826 rs12331507b rs9717730 rs4437516 rs73007780 rs7742463 rs73093677 N/A rs1333052 rs2804241 rs7397212
1.11 0.96 1.05 0.94 0.94 1.04 0.94 1.08 0.95 1.07 1.04 0.88 0.94 0.93 – 1.06 1.05 1.13
1.07–1.15 0.93–0.99 1.01–1.10 0.91–0.98 0.89–1.00 1.01–1.08 0.90–0.98 1.02–1.14 0.92–0.98 1.02–1.13 1.01–1.07 0.83–0.95 0.90–0.98 0.88–0.99 – 1.03–1.09 1.00–1.11 1.07–1.20
4.0 107 .017 .004 .004 .033 .015 .005 .007 1.0 103 .012 .018 3.3 104 .002 .023 – 2.5 104 .041 1.2 105
A B C D E F G H I J K L M N O P Q R
Note: ‘‘Significant’’ reflects a P%.05. N/A means there are no significant SNPs in that region. a OR ¼ per-allele odds ratio, adjusted for the first five principal components; based on study set-specific ORs that were meta-analyzed. b SNP was genotyped. Lee. Endometriosis SNPs and ovarian cancer. Fertil Steril 2015.
VOL. - NO. - / - 2015
5
FLA 5.4.0 DTD FNS29891_proof 21 October 2015 9:02 am ce JS
ORIGINAL ARTICLE: GENETICS
print & web 4C=FPO
FIGURE 1
Manhattan and linkage disequilibrium plots for region P. These plots depict the results for all 184 SNPs in region P (chromosome 9p21, position 22140648–22201586). The Manhattan plot includes ovarian cancer overall (‘‘invasive’’) and its four histotypes, with the x axis corresponding to the chromosomal position (in Mb), the y axis to the -log P, and the line to P¼5.0 108. Rs1333052, the most significant SNP in the region (P¼2.5 104 for ovarian cancer overall), is indicated. The linkage disequilibrium plot depicts pairwise correlation data from the 1000 Genomes CEU (Utah residents [CEPH] with Northern and Western European Ancestry) population. Epigenomic data from the ENCODE and Roadmap Epigenomics Consortia were obtained and visualized in the University of California, Santa Cruz Genome Browser. Lee. Endometriosis SNPs and ovarian cancer. Fertil Steril 2015.
histotypes. Other SNPs of notable significance were found in region R for high-grade serous cancer (rs6538605, P¼1.2 104) and endometrioid cancer (rs11107893, P¼4.3 104) as well as in regions A (rs4654785, P¼1.3 104) and L (rs71575922, P¼1.4 104) for clear cell cancer. Figure 1 presents the Manhattan and linkage disequilibrium plots for region P, the region showing the most significant burden statistic across all 18 regions, excluding region A (P¼ .002 for ovarian cancer overall). A clear elevation of P values is present, especially when looking at ovarian cancer overall, with the most significant SNP rs1333052 (P¼2.5 104) indicated. In addition, Supplemental Figure 1 presents the Manhattan plots for the remaining regions with significant burden statistics, whereas Supplemental Figure 2 presents these plots for the regions that did not have significant burden statistics. With the exception of region L showing some elevation of P values for clear cell cancer, reasonable given its borderline significant burden statistic (P¼ .07), all plots in Supplemental Figure 2 look relatively flat.
DISCUSSION We took a novel approach toward evaluating endometriosisassociated SNPs with ovarian cancer risk and found appreciable support for a shared genetic etiology between these two gynecologic conditions. Across the 18 regions harboring
putative endometriosis SNPs, we calculated a total of 15 significant burden statistics for ovarian cancer risk, compared with approximately 5 expected by chance at the P%.05 level (i.e., 18 regions 5 types of ovarian cancer (overall, highgrade serous, mucinous, endometrioid, clear cell) ¼ 90 total), a conservative estimate owing to the strong correlation between overall and high-grade serous. Endometriosis and ovarian cancer were first linked because of their frequent co-occurrence in surgical specimens. Most recently a pooled analysis of 13 case–control studies by Pearce et al. (3) showed that after adjusting for oral contraceptive use and parity, women with a history of endometriosis were 46% more likely to develop ovarian cancer, with the association primarily restricted to the endometrioid and clear cell histotypes. Endometriosis may be a precursor lesion for endometrioid and clear cell ovarian cancer, but the process of its malignant transformation is not well understood (26, 27). The two gynecologic conditions are likely to have a shared pathophysiology, given evidence from pathology case series reporting endometrioid and clear cell ovarian cancers arising from endometriotic foci as well as epidemiologic studies that highlight their similar hormonerelated risk factors, such as nulliparity (27). In addition, both traits seem to thrive in similar hormonal and immune environments, highlighting a possible shared inflammatory etiology. VOL. - NO. - / - 2015
6
FLA 5.4.0 DTD FNS29891_proof 21 October 2015 9:02 am ce JS
Fertility and Sterility® Some of these results are contrary to what we expected. Consistent with our hypothesis, we calculated a significant global burden statistic for clear cell ovarian cancer when considering all endometriosis SNPs together. However, we did not see this association for the endometrioid histotype despite its well-established association with endometriosis. In addition, we observed a link with the high-grade serous histotype, which previous epidemiologic studies have not found to be associated with endometriosis; we did not observe any associations with the low-grade serous histotype (data not shown owing to small numbers). Our findings suggest that although endometriosis is not a precursor lesion for high-grade serous ovarian cancer, the genetic pathways related to risk of these two diseases are shared. This finding paves the way for interesting lines of inquiry related to shared pathways. Perhaps women with endometriosis are more likely to have endosalpingiosis, which is common but infrequently noted in pathology reports, and which may predispose to high-grade serous ovarian cancer. Our ovarian cancer susceptibility GWASs have identified thousands of nominally significant associations, but most are due to chance. Deciphering the noise from true signals is difficult despite our large sample size. In this analysis we took a novel, two-pronged approach, first selecting candidate regions from a phenotype linked to ovarian cancer and then assessing whether the burden of associations in these regions was significant. The burden statistics implicate some of these regions in risk of ovarian cancer, although the exact relevance of the regions remains unknown. Region A includes the most significantly associated SNP among all 6,981 considered, rs3754496 (P¼5.4 108 for high-grade serous), an imputed SNP located near WNT4, a gene involved in steroidogenesis and in the development of the ovarian follicle and the female reproductive tract, biological functions that make its role in the development of both endometriosis and ovarian cancer compelling (28). Our data from OCAC for region A have previously been combined with those of CIMBA, and a genome-wide significant association was observed (16). In addition, Table 1 shows that an association with the remaining regions may exist even after excluding region A (P¼ .052). Only 3 of the 15 significant burden statistics presented in Table 2 are found in region A, with the other 12 spanning seven other regions. Region P had the most significant burden statistic (P¼ .001 for ovarian cancer overall), with rs1333052 as its most significant SNP (P¼2.5 104). This SNP is located adjacent to CDKN2B, a tumor suppressor gene whose methylation has been commonly seen in ovarian carcinogenesis (29). In addition, rs1333052 lies within a binding site for GATA2, one of six factors that constitute the GATA family of transcriptional regulatory proteins, which has been shown to play a role in ovarian function; it was identified by genome-wide ChIP sequencing (30). Rs12331507 in region I, which, other than region A, had the greatest number of significant burden statistics, lies in an intergenic region approximately 20 kb upstream of KDR, a gene that encodes one of the two receptors for
vascular endothelial growth factor (VEGF); VEGF has been shown to be expressed by tumor cells in ovarian cancer (31). Interestingly, the most significant SNP in region C, rs142034631, is located in GREB1, an estrogen-responsive gene that modulates tumor progression in models of ovarian cancer (32). Although the specific relevance of these genes and regions as well as the functional significance of their most significantly associated SNPs cannot be determined at this time, these results suggest that additional real associations exist between ovarian cancer risk and variation in the endometriosis-related regions. More importantly, they highlight the presence of a genetic commonality that high-grade serous and clear cell ovarian cancers and endometriosis are likely to share. In conclusion, we have shown significant evidence of a genetic link between endometriosis and ovarian cancer that warrants further follow-up. Whether the association between these two gynecologic conditions is causal remains unknown. However, fine mapping studies of the regions we have identified will greatly contribute to our knowledge regarding the etiologies of both diseases, as well as shed some light on their likely shared pathophysiology. The clinical implications of these results remain to be established, but next steps would include fine mapping of the regions with significant burden statistics (regions A, B, C, D, F, I, P, and R) and functional analyses of the most likely causal SNP(s). Acknowledgments: Ovarian Cancer Association Consortium: Roberta B. Ness, Agnieszka DansonkaMieszkowska, Aleksandra Gentry-Maharaj, Alexander Hein, Alice S. Whittemore, Allan Jensen, Andreas du Bois, Angela Brooks-Wilson, Anja Rudolph, Anna Jakubowska, Anna H. Wu, Argyrios Ziogas, Arif B. Ekici, Arto Leminen, Australian Cancer Study (Ovarian Cancer), Australian Ovarian Cancer Study Group, Barry Rosen, Beata Spiewankiewicz, Beth Y. Karlan, Britton Trabert, Brooke L. Fridley, C. Blake Gilks, Camilla Krakstad, Catherine M. Phelan, Cezary Cybulski, Christine Walsh, Claus Hogdall, Daniel W. Cramer, David G. Huntsman, Diana Eccles, Diether Lambrechts, Dong Liang, Douglas A. Levine, Edwin S. Iversen, Elisa V. Bandera, Elizabeth M. Poole, Ellen L. Goode, Els Van Nieuwenhuysen, Estrid Hogdall, Fiona Bruinsma, Florian Heitz, Francesmary Modugno, Graham G. Giles, Harvey A. Risch, Helen Baker, Helga B. Salvesen, Heli Nevanlinna, Hoda Anton-Culver, Honglin Song, Iain McNeish, Ian G. Campbell, Ignace Vergote, Ingo B. Runnebaum, Ingvild L. Tangen, Ira Schwaab, Jacek Gronwald, James Paul, Jan Lubinski, Jennifer A. Doherty, Jenny Chang-Claude, Jenny Lester, Joellen M. Schildkraut, John R. McLaughlin, Jolanta Lissowska, Jolanta Kupryjanczyk, Jonathan Tyrer, Joseph L. Kelley, Joseph H. Rothstein, Julie M. Cunningham, Karen Lu, Karen Carty, Kathryn L. Terry, Katja K.H. Aben, Kirsten B. Moysich, Kristine G. Wicklund, Kunle Odunsi, Lambertus A. Kiemeney, Lara Sucheston-Campbell, Lene Lundvall, Leon F.A.G. Massuger, Liisa M. Pelttari, Linda E. Kelemen, Linda S. Cook, Line Bjorge, Lotte Nedergaard, Louise A. Brinton, Lynne R. Wilkens, Malcolm C. Pike, Marc T. Goodman, Maria Bisogna, Mary Anne Rossing, Matthias
VOL. - NO. - / - 2015
7
FLA 5.4.0 DTD FNS29891_proof 21 October 2015 9:02 am ce JS
ORIGINAL ARTICLE: GENETICS W. Beckmann, Matthias D€ urst, Melissa C. Southey, Melissa Kellar, Michelle A.T. Hildebrandt, Nadeem Siddiqui, Natalia Antonenkova, Natalia Bogdanova, Nhu D. Le, Nicolas Wentzensen, Pamela J. Thompson, Patricia Harrington, Penelope M. Webb, Peter A. Fasching, Peter Hillemanns, Philipp Harter, Piotr Sobiczewski, Rachel Palmieri Weber, Ralf Butzow, Robert P. Edwards, Robert A. Vierkant, Rosalind Glasspool, Sandra Orsulic, Sandrina Lambrechts, Sara H. Olson, Shan Wang-Gohrke, Shashi Lele, Shelley S. Tworoger, Simon A. Gayther, Stacey A. Missmer, Steven A. Narod, Susan J. Ramus, Susanne K. Kjaer, Tanja Pejovic, Thilo D€ ork, Ursula Eilber, Usha Menon, Valerie McGuire, Weiva Sieh, Xifeng Wu, Yukie Bean, Yurii B Shvetsov. The authors thank all the individuals who took part in this study and all the researchers, clinicians, and technical and administrative staff who have made possible the many studies contributing to this work. In particular, the authors thank D. Bowtell, A. deFazio, D. Gertig, A. Green, P. Parsons, N. Hayward, and D. Whiteman (AUS) (see Supplemental Table 1 for site abbreviations); G. Peuteman, T. Van Brussel, and D. Smeets (BEL); the staff of the genotyping unit, S. LaBoissiere and F. Robidoux (Genome Quebec); L. Gacucova (HMO); P. Schurmann, F. Kramer, W. Zheng, T. W. Park, Simon, K. Beer-Grondke, and D. Schmidt (HJO); S. Windebank, C. Hilker, and J. Vollenweider (MAY); the state cancer registries of AL, AZ, AR, CA, CO, CT, DE, FL, GA, HI, ID, IL, IN, IA, KY, LA, ME, MD, MA, MI, NE, NH, NJ, NY, NC, ND, OH, OK, OR, PA, RI, SC, TN, TX, VA, WA, and WYL (NHS); L. Paddock, M. King, L. Rodriguez-Rodriguez, A. Samoila, and Y. Bensman (NJO); M. Sherman, A. Hutchinson, N. SzeszeniaDabrowska, B. Peplonska, W. Zatonski, A. Soni, P. Chao, and M. Stagner (POL); C. Luccarini, P. Harrington the SEARCH team and ECRIC (SEA); I. Jacobs, M. Widschwendter, E. Wozniak, N. Balogun, A. Ryan, C. Karpinskyj, and J. Ford (UKO); Carole Pye (UKR); A. Amin Al Olama, J. Dennis, E. Dicks, K. Michilaidou, K. Kuchenbaker (COGS); and G. Montgomery.
REFERENCES 1. 2. 3.
4.
5.
6.
7.
Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global cancer statistics. CA Cancer J Clin 2011;61:69–90. Bulletti C, Coccia ME, Battistoni S, Borini A. Endometriosis and infertility. J Assist Reprod Genet 2010;27:441–7. Pearce CL, Templeman C, Rossing MA, Lee A, Near AM, Webb PM, et al. Association between endometriosis and risk of histological subtypes of ovarian cancer: a pooled analysis of case-control studies. Lancet Oncol 2012;13: 385–94. Treloar SA, O'Connor DT, O'Connor VM, Martin NG. Genetic influences on endometriosis in an Australian twin sample. Fertil Steril 1999;71: 701–10. Matalliotakis IM, Arici A, Cakmak H, Goumenou AG, Koumantakis G, Mahutte NG. Familial aggregation of endometriosis in the Yale Series. Arch Gynecol Obstet 2008;278:507–11. Uno S, Zembutsu H, Hirasawa A, Takahashi A, Kubo M, Akahane T, et al. A genome-wide association study identifies genetic variants in the CDKN2BAS locus associated with endometriosis in Japanese. Nat Genet 2010;42:707–10. Painter JN, Anderson CA, Nyholt DR, Macgregor S, Lin J, Lee SH, et al. Genome-wide association study identifies a locus at 7p15.2 associated with endometriosis. Nat Genet 2011;43:51–4.
8.
Adachi S, Tajima A, Quan J, Haino K, Yoshihara K, Masuzaki H, et al. Metaanalysis of genome-wide association scans for genetic susceptibility to endometriosis in Japanese population. J Hum Genet 2010;55:816–21. 9. Bojesen SE, Pooley KA, Johnatty SE, Beesley J, Michailidou K, Tyrer JP, et al. Multiple independent variants at the TERT locus are associated with telomere length and risks of breast and ovarian cancer. Nat Genet 2013;45: 371–84. 384e1–2. 10. Bolton KL, Ganda C, Berchuck A, Pharaoh PD, Gayther SA. Role of common genetic variants in ovarian cancer susceptibility and outcome: progress to date from the Ovarian Cancer Association Consortium (OCAC). J Intern Med 2012;271:366–78. 11. Kuchenbaecker KB, Ramus SJ, Tyrer J, Lee A, Shen HC, Beesley J, et al. Identification of six new susceptibility loci for invasive epithelial ovarian cancer. Nat Genet 2015;47:164–71. 12. Nyholt DR, Low SK, Anderson CA, Painter JN, Uno S, Morris AP, et al. Genome-wide association meta-analysis identifies new endometriosis risk loci. Nat Genet 2012;44:1355–9. 13. Pearce CL, Rossing MA, Lee AW, Ness RB, Webb PM, for Australian Cancer Study (Ovarian Cancer), et al. Combined and interactive effects of environmental and GWAS-identified risk factors in ovarian cancer. Cancer Epidemiol Biomarkers Prev 2013;22:880–90. 14. Song H, Ramus SJ, Tyrer J, Bolton KL, Gentry-Maharaj A, Wozniak E, et al. A genome-wide association study identifies a new ovarian cancer susceptibility locus on 9p22.2. Nat Genet 2009;41:996–1000. 15. Goode EL, Chenevix-Trench G, Song H, Ramus SJ, Notaridou M, Lawrenson K, et al. A genome-wide association study identifies susceptibility loci for ovarian cancer at 2q31 and 8q24. Nat Genet 2010;42:874–9. 16. Bolton KL, Tyrer J, Song H, Ramus SJ, Notaridou M, Jones C, et al. Common variants at 19p13 are associated with susceptibility to ovarian cancer. Nat Genet 2010;42:880–4. 17. Permuth-Wey J, Lawrenson K, Shen HC, Velkova A, Tyrer JP, Chen Z, et al. Identification and molecular characterization of a new ovarian cancer susceptibility locus at 17q21.31. Nat Commun 2013;4:1627. 18. Pharoah PD, Tsai YY, Ramus SJ, Phelan CM, Goode EL, Lawrenson K, et al. GWAS meta-analysis and replication identifies three new susceptibility loci for ovarian cancer. Nat Genet 2013;45:362–70. 370e1–2. 19. Shen H, Fridley BL, Song H, Lawrenson K, Cunningham JM, Ramus SJ, et al. Epigenetic analysis leads to identification of HNF1B as a subtype-specific susceptibility gene for ovarian cancer. Nat Commun 2013;4:1628. 20. Lee SH, Harold D, Nyholt DR, ANZGene Consortium, International Endogene Consortium, Genetic and Environmental Risk for Alzheimer's Disease Consortium, et al. Estimation and partitioning of polygenic variation captured by common SNPs for Alzheimer's disease, multiple sclerosis and endometriosis. Hum Mol Genet 2013;22:832–41. 21. Permuth-Wey J, Kim D, Tsai YY, Lin HY, Chen YA, Barnholtz-Sloan J, et al. LIN28B polymorphisms influence susceptibility to epithelial ovarian cancer. Cancer Res 2011;71:3896–903. 22. Howie B, Marchini J, Stephens M. Genotype imputation with thousands of genomes. G3 (Bethesda) 2011;1:457–70. 23. Howie B, Fuchsberger C, Stephens M, Marchini J, Abecasis GR. Fast and accurate genotype imputation in genome-wide association studies through pre-phasing. Nat Genet 2012;44:955–9. 24. Sankararaman S, Sridhar S, Kimmel G, Halperin E. Estimating local ancestry in admixed populations. Am J Hum Genet 2008;82:290–303. 25. Tyrer J, Pharoah PD, Easton DF. The admixture maximum likelihood test: a novel experiment-wise test of association between disease and multiple SNPs. Genet Epidemiol 2006;30:636–43. 26. Nezhat F, Datta MS, Hanson V, Pejovic T, Nezhat C, Nezhat C. The relationship of endometriosis and ovarian malignancy: a review. Fertil Steril 2008; 90:1559–70. 27. Ness RB. Endometriosis and ovarian cancer: thoughts on shared pathophysiology. Am J Obstet Gynecol 2003;189:280–94. 28. Franco HL, Dai D, Lee KY, Rubel CA, Roop D, Boerboom D, et al. WNT4 is a key regulator of normal postnatal uterine development and progesterone signaling during embryo implantation and decidualization in the mouse. FASEB J 2011;25:1176–87.
VOL. - NO. - / - 2015
8
FLA 5.4.0 DTD FNS29891_proof 21 October 2015 9:02 am ce JS
Fertility and Sterility® 29.
30.
Ozdemir F, Altinisik J, Karateke A, Coksuer H, Buyru N. Methylation of tumor suppressor genes in ovarian cancer. Exp Ther Med 2012;4: 1092–6. Viger RS, Guittot SM, Anttonen M, Wilson DB, Heikinheimo M. Role of the GATA family of transcription factors in endocrine development, function, and disease. Mol Endocrinol 2008;22:781–98.
31.
32.
Boocock CA, Charnock-Jones DS, Sharkey AM, McLaren J, Barker PJ, Wright KA, et al. Expression of vascular endothelial growth factor and its receptors flt and KDR in ovarian carcinoma. J Natl Cancer Inst 1995;87:506–16. Laviolette LA, Hodgkinson KM, Minhas N, Perez-Iratxeta C, Vanderhyden BC. 17beta-estradiol upregulates GREB1 and accelerates ovarian tumor progression in vivo. Int J Cancer 2014;135:1072–84.
VOL. - NO. - / - 2015
9
FLA 5.4.0 DTD FNS29891_proof 21 October 2015 9:02 am ce JS
ORIGINAL ARTICLE: GENETICS
SUPPLEMENTAL FIGURE 1
(A–F) Manhattan plots for the regions that resulted in a significant burden statistic at the P%.05 level. These plots include all SNPs in each region for ovarian cancer overall (‘‘invasive’’) and its four histotypes, with the x axis corresponding to the chromosomal position (in Mb), the y axis to the -log P, and the line to P¼5.0 108. Region A, which is presented in a previous publication (16), and region P, which is shown in the article, are not included. Lee. Endometriosis SNPs and ovarian cancer. Fertil Steril 2015.
VOL. - NO. - / - 2015
9.e1
FLA 5.4.0 DTD FNS29891_proof 21 October 2015 9:02 am ce JS
Fertility and Sterility®
web 4C=FPO
SUPPLEMENTAL FIGURE 1 Continued
Lee. Endometriosis SNPs and ovarian cancer. Fertil Steril 2015.
VOL. - NO. - / - 2015
9.e2
FLA 5.4.0 DTD FNS29891_proof 21 October 2015 9:02 am ce JS
ORIGINAL ARTICLE: GENETICS
SUPPLEMENTAL FIGURE 2
(A–J) Manhattan plots for the regions that did not result in a significant burden statistic at the P%.05 level. These plots include all SNPs in each region for ovarian cancer overall (‘‘invasive’’) and its four histotypes, with the x axis corresponding to the chromosomal position (in Mb), the y axis to the -log P, and the line to P¼5.0 108. Lee. Endometriosis SNPs and ovarian cancer. Fertil Steril 2015.
VOL. - NO. - / - 2015
9.e3
FLA 5.4.0 DTD FNS29891_proof 21 October 2015 9:02 am ce JS
Fertility and Sterility®
SUPPLEMENTAL FIGURE 2 Continued
Lee. Endometriosis SNPs and ovarian cancer. Fertil Steril 2015.
VOL. - NO. - / - 2015
9.e4
FLA 5.4.0 DTD FNS29891_proof 21 October 2015 9:02 am ce JS
ORIGINAL ARTICLE: GENETICS
web 4C=FPO
SUPPLEMENTAL FIGURE 2 Continued
Lee. Endometriosis SNPs and ovarian cancer. Fertil Steril 2015.
VOL. - NO. - / - 2015
9.e5
FLA 5.4.0 DTD FNS29891_proof 21 October 2015 9:02 am ce JS
VOL. - NO. - / - 2015
SUPPLEMENTAL TABLE 1 Description of OCAC study sites included in analyses. Study set abbreviation AUS
BAV FLA 5.4.0 DTD FNS29891_proof 21 October 2015 9:02 am ce JS
BEL BWH1 DAN DOV
GER HAW HJO HMO HOC HPE
LA2
MSK
Site abbreviation
Australian Ovarian Cancer Study/Australian Cancer Study (Ovarian Cancer) Melbourne Collaborative Cohort Study Bavarian Ovarian Cancer Cases and Controls Belgian Ovarian Cancer Study New England Case–Control Study Danish Malignant Ovarian Tumor Study Danish Pelvic Mass Study Diseases of the Ovary and their Evaluation Oregon Ovarian Cancer Registry German Ovarian Cancer Study Dr. Horst Schmidt Kliniken Hawaii Ovarian Cancer Case–Control Study Hannover-Jena Ovarian Cancer Study Hannover-Minsk Ovarian Cancer Study Helsinki Ovarian Cancer Study Gilda Radner Familial Ovarian Cancer Registry Hormones and Ovarian Cancer Prediction
AUS
Australia
Population-based/case– control
977
880
498
42
106
52
MCC
Melbourne, Australia
66
66
19
7
7
6
BAV
Southeast Germany
143
96
42
8
13
6
BEL
Belgium
1,349
276
183
23
22
23
NEC
142
104
103
0
0
0
MAL
USA: NH and Eastern MA Denmark
828
440
183
42
54
33
PVD DOV
Denmark USA: Western WA
0 1,487
169 904
119 435
10 26
14 148
8 65
ORE
USA: Portland, OR
Cohort/Nested case– control Population-based/case– control Hospital-based/case– control Population-based/case– control Population-based/case– control Case-only Population-based/case– control Case-only
0
55
37
2
5
2
GER
Southern Germany
413
189
68
21
21
6
HSK HAW
0 157
147 60
103 36
1 3
16 12
0 5
HJO
Germany USA: HI and Southern CA Germany
273
261
107
9
26
4
HMO
Belarus
Population-based/case– control Case-only Population-based/case– control Hospital-based/case– control Case-control
138
142
1
7
12
1
HOC
Helsinki, Finland
Case-control
447
210
0
39
26
11
GRR
USA
Familial cancer/case-only
0
134
0
5
18
11
HOP
USA: Western PA, Northeastern OH, Western NY USA: Southern CA
Population-based/case– control
1,464
568
275
29
83
42
Case-only
0
277
201
12
18
9
USC
USA: Southern CA (Los Angeles County)
Population-based/case– control
1,047
693
341
44
79
35
MDA
USA: TX
384
375
135
27
28
4
MSK
USA: NY
Hospital-based/case– control Case–control
593
467
343
0
20
18
Women's Cancer Research Institute - Cedars-Sinai Medical Center Los Angeles County Case– Control Studies of Ovarian Cancer MD Anderson Ovarian Cancer Study Memorial Sloan Kettering Cancer Center Gynecology Tissue Bank
9.e6
Lee. Endometriosis SNPs and ovarian cancer. Fertil Steril 2015.
LAX
Location
Site type
Controls Cases
Fertility and Sterility®
MDA
High-grade Clear serous Mucinous Endometrioid cell
Site namea
Continued. Study set abbreviation MYO1 NCO NEC NHS FLA 5.4.0 DTD FNS29891_proof 21 October 2015 9:02 am ce JS
NOR NTH OVA PLD
POL1 SEA
STA
UCI
Site namea
Site abbreviation
Mayo Clinic Ovarian Cancer Case–Control Study North Carolina Ovarian Cancer Study New England Case–Control Study Nurses' Health Study
MAY
NHS
USA: NH and Eastern MA USA
University of Bergen Norway Study Nijmegen Polygene Study & Nijmegen Biomedical Study Ovarian Cancer in Alberta and British Columbia Study Poland Ovarian Cancer Study
NOR
NCI Ovarian Case–Control Study in Poland Warsaw Ovarian Cancer Study NCI Ovarian Case–Control Study in Poland UK Studies of Epidemiology and Risk Factors in Cancer Heredity (SEARCH) Ovarian Cancer Study Familial Registry for Ovarian Cancer and Genetic Epidemiology of Ovarian Cancer UC Irvine Ovarian Cancer Study
Lee. Endometriosis SNPs and ovarian cancer. Fertil Steril 2015.
Location
Controls Cases
High-grade Clear serous Mucinous Endometrioid cell
Clinic-based/case–control
441
443
343
11
34
20
172
266
133
16
35
24
980
618
331
41
140
33
455
138
0
7
14
6
Norway
Population-based/case– control Population-based/case– control Population-based/nested case–control Case-control
371
237
85
15
27
11
NTH
Eastern Netherlands
Case-control
323
295
64
33
64
20
OVA
Canada: Alberta and British Columbia
Case-control
748
640
0
26
103
57
POC
Poland: Szczecin, Poznan, Opole, Rzeszow Poland: Warsaw, Lodz
Case-control
417
422
0
33
39
9
Population-based/case– control Case-control
186
42
15
4
10
2
204
202
131
8
20
17
Population-based/case– control Population-based/case– control
555
212
57
13
28
9
6023
316
108
43
29
28
NCO NEC
POL WOC POL
USA: MN, SD, ND, IL, IA, WI USA: NC
Site type
Poland: Warsaw, Central Poland Poland, Warsaw, Lodz
SEA
UK: East Anglia and West Midlands
STA
USA: San Francisco Bay
Population-based/case– control
313
251
135
16
32
20
UCI
USA: Southern CA (Orange, San Diego, Imperial Counties)
Population-based/case– control
367
277
143
19
48
23
ORIGINAL ARTICLE: GENETICS
9.e7
SUPPLEMENTAL TABLE 1
VOL. - NO. - / - 2015
VOL. - NO. - / - 2015
SUPPLEMENTAL TABLE 1 Continued. Study set abbreviation UK
FLA 5.4.0 DTD FNS29891_proof 21 October 2015 9:02 am ce JS
UK2
US
USN
b
Site abbreviation
Royal Marsden Hospital Case Series UK Studies of Epidemiology and Risk Factors in Cancer Heredity (SEARCH) Ovarian Cancer Study UK Ovarian Cancer Population Study UK Familial Ovarian Cancer Registry Welcome Trust Case– Control Consortium Royal Marsden Hospital Case Series Southampton Ovarian Cancer Study Scottish Randomised Trial in Ovarian Cancer UK Ovarian Cancer Population Study UK Familial Ovarian Cancer Registry Mayo Clinic Hospital
RMH
UK: London
Hospital-based/case-only
0
138
0
15
26
16
SEA
UK: East Anglia, West Midlands
Population-based/case– control
1
1,130
346
111
205
127
UKO
Population-based/case– control Case-only/familial register
1
494
186
53
82
42
UKR
UK: England, Wales, Northern Ireland UK: National
0
44
0
3
6
2
WTC
UK
6,116
0
0
0
0
0
RMH
UK: London
Blood donors, 1958 Birth Cohort Controls Hospital-based/case-only
0
10
0
1
1
1
SOC
UK: Wessex
Hospital-based/case-only
0
274
72
33
62
11
SRO
UK: Glasgow, patients recruited worldwide UK: England, Wales, Northern Ireland UK: National
Case-only from clinical trial
0
158
87
3
17
8
Population-based/case– control Case-only/familial register
1,103
238
93
24
32
25
0
3
0
0
0
0
Clinic-based/case–control
519
358
233
11
71
25
Population-based/case– control Population-based/case– control Population-based/case– control Clinic-based/case–control
654
487
249
18
78
56
169
227
129
13
28
12
524
729
0
56
154
46
10
10
9
0
1
0
Case-control
181
169
80
7
27
20
74
20
0
1
7
2
North Carolina Ovarian Cancer Study Tampa Bay Ovarian Cancer Study Toronto Ovarian Cancer Study Mayo Clinic Ovarian Cancer Case–Control Study New Jersey Ovarian Cancer Study Toronto Ovarian Cancer Study
UKO UKR MAY
Location
NCO
USA: MN, SD, ND, IL, IA, WI USA: NC
TBO
USA: FL
TOR
Canada: Ontario
MAY NJO
USA: MN, SD, ND, IL, IA, WI USA: NJ
TOR
Canada: Ontario
Some sites appear in more than one study set owing to differences in genotyping platform within each site. Sum of numbers do not equal total number of cases because some cases were not classified as one of the four histotypes.
Lee. Endometriosis SNPs and ovarian cancer. Fertil Steril 2015.
Site type
Population-based/case– control Total
Controls Cases
30,815 15,361
b
6,258
991
b
b
2,148
1,013b
9.e8
Fertility and Sterility®
a
High-grade Clear serous Mucinous Endometrioid cell
Site namea
ORIGINAL ARTICLE: GENETICS
SUPPLEMENTAL TABLE 2 Odds ratios and 95% CIs for the most significant SNP in each of the 18 genomic regions for high-grade serous and mucinous ovarian cancers. High-grade serous (n [ 6,258) No. of significant Region SNPs A B C D E F G H I J K L M N O P Q R
221 50 0 0 0 3 6 1 101 0 1 13 0 0 1 52 0 327
Most significant SNP rs3754496 rs1524169 N/A N/A N/A rs1846555 rs73063613 chr3:104626317:I rs1547512 N/A chr6:113076251:D chr6:152616764:D N/A N/A rs13311403 rs10965263 N/A rs6538605
ORa 1.16 1.05 – – – 1.05 1.07 1.08 1.06 – 1.08 0.87 – – 1.07 1.07 – 0.92
95% CI
Mucinous (n [ 991)
P value 8
1.10–1.23 5.6 10 1.01–1.10 9.7 103 – – – – – – 1.00–1.10 .031 1.01–1.13 .015 1.01–1.15 .031 1.02–1.10 .006 – – 1.00–1.15 .037 0.80–0.96 .003 – – – – 1.01–1.14 .024 1.02–1.12 .005 – – 0.88–0.96 1.2 104
No. of significant SNPs
Most s ignificant SNP
ORb
4 2 27 3 14 23 0 5 2 1 23 6 3 0 15 1 35 14
rs2501276 chr1:80243322:I rs142034631 rs71413109 rs12487626 rs4858684 N/A rs116203658 rs62304730 chr6:19617536:I rs7751257 rs12206761 rs7742463 N/A rs6954307 rs1930591c rs2769419 rs7397485
0.80 0.90 1.15 0.74 0.89 0.85 – 0.76 1.14 1.11 1.21 1.18 0.81 – 0.79 0.91 1.13 0.86
95% CI
P value
0.67–0.95 .012 0.82–0.99 .028 1.02–1.29 .020 0.60–0.91 .005 0.80–0.98 .024 0.76–0.94 .003 – – 0.60–0.97 .028 1.02–1.28 .024 1.01–1.22 .034 1.08–1.34 5.8 104 1.04–1.34 .013 0.72–0.92 7.1 104 – – 0.67–0.93 .005 0.83–1.00 .040 1.03–1.24 .009 0.78–0.96 .006
Note: ‘‘Significant’’ reflects a P%.05. N/A means there are no significant SNPs in that region. a OR ¼ Per-allele odds ratio, adjusted for the first five principal components; based on study set–specific ORs that were meta-analyzed. b OR ¼ Per-allele odds ratio, adjusted for the first five principal components and study set. c SNP was genotyped. Lee. Endometriosis SNPs and ovarian cancer. Fertil Steril 2015.
VOL. - NO. - / - 2015
9.e9
FLA 5.4.0 DTD FNS29891_proof 21 October 2015 9:02 am ce JS
Fertility and Sterility®
SUPPLEMENTAL TABLE 3 Odds ratios and 95% CIs for the most significant SNP in each of the 18 genomic regions for endometrioid and clear cell ovarian cancers. Endometrioid (n [ 2,148)
Region A B C D E F G H I J K L M N O P Q R
No. of significant SNPs
Most significant SNP
ORa
95% CI
49 0 47 3 13 53 1 4 9 1 5 18 0 1 72 6 0 106
rs12024016 N/A rs1898003 rs17034338 rs7639387 rs1603996 rs17054879 rs75281964 rs2412627 rs10434838 rs34384052 rs9479288 N/A rs73093677 rs12530891b rs6475610b N/A rs11107893
1.17 – 1.09 0.86 0.87 1.15 0.91 1.14 0.86 1.08 1.14 0.82 – 0.84 1.10 1.09 – 1.18
1.02–1.33 – 1.01–1.17 0.75–0.99 0.77–0.99 1.06–1.24 0.84–1.00 1.01–1.28 0.79–0.94 1.01–1.16 1.03–1.25 0.72–0.94 – 0.73–0.95 1.03–1.18 1.02–1.16 – 1.08–1.30
Clear cell (n [ 1,013)
P value
No. of significant SNPs
Most significant SNP
ORa
95% CI
P value
.023 – .020 .038 .030 .001 .050 .033 .001 .029 .009 .004 – .008 .006 .008 – 4.3 104
176 210 51 0 0 3 1 25 52 40 1 57 0 2 5 27 0 93
rs4654785 chr1:80301656:D rs10929757 N/A N/A rs62228504 rs6767715 rs77484561 rs12331507b rs2743617 rs72942081 rs71575922 N/A rs73093677 rs3800729b rs10965256b N/A chr12:95380417:D
1.22 0.87 1.18 – – 0.83 1.13 1.22 0.85 1.24 1.26 1.25 – 0.80 0.87 0.73 – 0.86
1.10–1.35 0.79–0.95 1.07–1.31 – – 0.70–0.99 1.02–1.26 1.04–1.43 0.77–0.93 1.04–1.49 1.03–1.55 1.12–1.41 – 0.66–0.97 0.77–0.98 0.60–0.89 – 0.78–0.96
1.3 104 .003 .001 – – .033 .020 .013 8.3 104 .017 .028 1.4 104 – .023 .021 .002 – .005
Note: ‘‘Significant’’ reflects a P%.05. N/A means there are no significant SNPs in that region. a OR ¼ Per-allele odds ratio, adjusted for the first five principal components and study set. b SNP was genotyped. Lee. Endometriosis SNPs and ovarian cancer. Fertil Steril 2015.
VOL. - NO. - / - 2015
9.e10
FLA 5.4.0 DTD FNS29891_proof 21 October 2015 9:02 am ce JS