PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of hmgLink to Publisher's site
 
Hum Mol Genet. 2009 June 15; 18(12): 2297–2304.
Published online 2009 March 20. doi:  10.1093/hmg/ddp138
PMCID: PMC2685754
EMSID: UKMS4378

Association between invasive ovarian cancer susceptibility and 11 best candidate SNPs from breast cancer genome-wide association study

Abstract

Because both ovarian and breast cancer are hormone-related and are known to have some predisposition genes in common, we evaluated 11 of the most significant hits (six with confirmed associations with breast cancer) from the breast cancer genome-wide association study for association with invasive ovarian cancer. Eleven SNPs were initially genotyped in 2927 invasive ovarian cancer cases and 4143 controls from six ovarian cancer case–control studies. Genotype frequencies in cases and controls were compared using a likelihood ratio test in a logistic regression model stratified by study. Initially, three SNPs (rs2107425 in MRPL23, rs7313833 in PTHLH, rs3803662 in TNRC9) were weakly associated with ovarian cancer risk and one SNP (rs4954956 in NXPH2) was associated with serous ovarian cancer in non-Hispanic white subjects (P-trend < 0.1). These four SNPs were then genotyped in an additional 4060 cases and 6308 controls from eight independent studies. Only rs4954956 was significantly associated with ovarian cancer risk both in the replication study and in combined analyses. This association was stronger for the serous histological subtype [per minor allele odds ratio (OR) 1.07 95% CI 1.01–1.13, P-trend = 0.02 for all types of ovarian cancer and OR 1.14 95% CI 1.07–1.22, P-trend = 0.00017 for serous ovarian cancer]. In conclusion, we found that rs4954956 was associated with increased ovarian cancer risk, particularly for serous ovarian cancer. However, none of the six confirmed breast cancer susceptibility variants we tested was associated with ovarian cancer risk. Further work will be needed to identify the causal variant associated with rs4954956 or elucidate its function.

INTRODUCTION

Recent advances in high-throughput genotyping technologies have enabled rapid and efficient genotyping to be performed for hundreds of thousands of genetic variants without prior knowledge of gene function as part of genome-wide association studies (GWAS). Several GWAS have led to the identification of novel loci for many different common complex diseases, including diabetes, Crohn's disease, rheumatoid arthritis and breast, prostate and colorectal cancers, confirming that susceptibility to these diseases has a polygenic component (14). The known ovarian cancer susceptibility genes, such as BRCA1 and BRCA2, appear to explain <40% of the excess familial risk of this disease (5). It is likely that a combination of multiple low or moderate penetrance genetic variants contribute to the remaining unexplained excess ovarian cancer risks. It is evident that some loci regulate carcinogenic pathways common to multiple cancers, for example, the variant rs6983267 in 8q24 is associated with colorectal, prostate and ovarian cancer risk (6). Also, BRCAC1 and BRCA2 gene mutations are associated with the variation in the risks of both breast and ovarian cancers. The development of breast and ovarian cancers both have a hormonal basis of female cancers, therefore, it seems logical to hypothesize that these malignancies share some common genetic risks and that breast cancer susceptibility loci may also be associated with the risk of ovarian cancer.

The aim of this study is to evaluate whether the 11 SNPs that were most strongly associated with breast cancer risk in our breast cancer GWAS (7) were associated with epithelial ovarian cancer (EOC) risk in 14 case–control studies which comprised 6987 invasive EOC cases and 10451 controls. This represents the largest ovarian cancer case–control analysis conducted to date. These 11 SNPs included six breast cancer susceptibility variants that reach genome-wide significance (P < 10−8) [i.e. rs2981582 (FGFR2), rs12443621 (TNRC9), rs13281615 (8q), rs3817198 (LSP1), rs3803662 (TNRC9) and rs889312 (MAP3K1)]. The remaining five top hits from breast cancer GWAS, namely rs4666451 (located on chromosome 2p), rs2107425 (MRPL23), rs7313833 (PTHLH), rs981782 (on chromosome 5p) and rs4954956 (NXPH2) are also strong candidate for breast cancer associations (P < 10−5).

RESULTS

We have genotyped the 11 SNPs identified through our GWAS for breast cancer in a set of six ovarian cancer case–control studies. Genotype distributions in controls were consistent with Hardy–Weinberg equilibrium (HWE) except for rs4666451 in the USC (P = 0.02) and UKO study

(P = 0.01) and rs7313833 in the AUS study (P = 0.04). These deviations are likely to be due to chance, rather than a reflection of poor genotyping, because inspection of the cluster plots indicated good discrimination between genotype; furthermore deviation from HWE was not observed in cases. In addition, the genotyping for rs981782 failed in the UKO study and rs7313833 failed in the GER study.

Genotype-specific odds ratios (ORs) and tests of association are presented in Table 1. The genotype-specific risks for serous ovarian cancer, estimated from the combined data, are also presented in Table 1. The observed genotype frequencies for each of the data sets are presented in Supplementary Material, Table S1. There was no association in controls between age and genotype frequency for any of the SNPs and age-adjusted genotype-specific risks were similar to the unadjusted ORs (data not shown). Two SNPs (rs7313833 and rs210742) showed some evidence of association with all types of invasive ovarian cancer, whereas rs4954956 and rs210742 showed some evidence of association with serous type ovarian. The association of rs3803662 with ovarian cancer risk was of borderline significance (P = 0.07) with ovarian cancer risks. There was no association for the remaining seven SNPs (P > 0.1).

Table 1.
Genotype-specific risks (95% CI) for all types of invasive ovarian cancer and serous type ovarian cancer

Carriers of the minor allele of rs7313833 were at increased risk of ovarian cancer: per minor allele OR = 1.09, 95% CI 1.01–1.18, P-trend = 0.027. Carriers of the minor allele of rs2107425 were at decreased risk of ovarian cancer overall and serous type ovarian cancer: per minor allele OR = 0.91 95% CI 0.84–0.98 (P = 0.027) and OR = 0.91 95% CI 0.82–1.00, (P-trend = 0.048), respectively. Carriers of the minor allele of rs4954956 were at increased risk of serous type ovarian cancer: per minor allele OR = 1.10 95% CI 1.02–1.19, (P-trend = 0.01).

The four SNPs (rs2107425, rs3803662, rs4954956 and rs7313833) with some evidence of association in either all cases or serous type ovarian cancer were then genotyped in the validation set. The results of the validation component alone and in combination with the initial set are presented in Table 2. SNP rs4954956 (NXPH2) was associated with overall ovarian cancer risk and with the risk of serous subtype both in the validation set alone and in the combined analysis, although the risk was stronger for the serous subtype (Table 2). The risk for serous ovarian cancer per minor allele of rs4954956 was 1.17 (95% CI 1.06–1.28, P = 0.0011) for the validation set and 1.14 (95% CI 1.07–1.22, P = 0.00017) for the combined analyses. For all types of ovarian cancer, per minor allele risk for rs4954956 was 1.10 (95% CI 1.02–1.19; P = 0.01) for the validation set and 1.07 (95% CI 1.01–1.13; P = 0.02) for the combined analyses. Figure 1 shows the genotype-specific ORs for each ovarian study and for the combined analysis for rs4954956. The effect of rs4954956 was slightly attenuated after adjusting for a first degree family history of breast cancer [per rare allele OR 1.06, 95% CI 0.97–1.16 (P = 0.17) for all types of ovarian cancer and 1.10, 95% CI 0.99–1.22 (P = 0.066) for serous type of ovarian cancer, respectively]. The remaining three SNPs were not validated. There was no evidence for between-study heterogeneity (P > 0.05) for all the SNPs tested except rs3817198 (P = 0.0002) in the initial set.

Figure 1.
Genotype-specific risks of SNP rs4954956 for ovarian cancer by study in White non-Hispanic subjects. (A) All ovarian cancer subtypes included. (B) Analysis restricted to serous type ovarian cancers.
Table 2.
Genotype-specific risks 95% CI for all types and serous type of invasive epithelial ovarian cancer in the validation (stage 2) studies for non-Hispanic White subjects

DISCUSSION

This is the largest ovarian cancer association study conducted to date involving 14 studies from Ovarian Cancer Association Consortium (OCAC) and comprising 5876 invasive EOC cases and 9273 controls of non-Hispanic origin. We observed an association for the minor allele of SNP rs4954956 with increased risks of EOC both in the initial set and in the replication set with the strongest gene-dose effect for serous type EOC. We urge caution in the interpretation of these results as a number of reported positive associations in the literature have not been replicated by the subsequent studies. Indeed, the proportion of studies with false-positive findings can be as high as 95% in association studies between genetic variants and disease risks (810). To estimate the likelihood that our results represent a true association with ovarian cancer risks, we calculated the false-positive report probability (FPRP) under different prior probability scenarios (11). The FPRP depends on the prior probability that a true association exists, the observed level of significance (α) and the statistical power to detect the OR of the alternative hypothesis at the given α. As there are a large number of common SNPs in the genome, the overall prior probability of association is very low (<1 in 106). However, the prior probability that rs4954956 is associated with ovarian cancer is more favourable as this is one of the best candidates from our breast cancer GWAS and therefore a good candidate for ovarian cancer susceptibility. The FPRPs for rs4954956 under various prior probabilities and the power to detect the association at our observed significance level α (assuming the true effect size is equal to that observed) are presented in Table 3. For example, assuming the prior probability to be 1 in 100 or 1 in 1000, the FPRP for association of rs4954956 with serous ovarian cancer would be 0.03 and 0.23, respectively. This, along with the fact that results were indeed independently replicated in the case–control validation studies (Table 2), suggests that this association is robust (Table 3). The evidence, however, is weaker for its association with all types of ovarian cancer under the same prior probabilities.

Table 3.
False-positive report probability values for rs4954956

Underlying population stratification is another explanation for a spurious association. This occurs when allele frequencies differ between population subgroups and cases and controls are drawn differentially from those subgroups. To minimize the impact of population stratification, analyses were restricted to White subjects with non-Hispanic origin. If population stratification were present, it is unlikely that the same degree of stratification would be found in all 14 studies. We did not observe any heterogeneity between different studies in the initial or replication studies or the combined analysis for rs4954956, thus providing evidence against substantial population stratification or other study-specific biases.

If the observed association is confirmed, the SNP may be directly causal or an indirect marker in linkage disequilibrium with the real cause of malignancy. SNP rs4954956 is in an intergenic region situated <7 kb upstream of the gene NXPH2, which encodes the protein neurexophilin 2. NXPH2 is expressed in kidney and brain and acts as a signalling molecule; it has also been shown to be expressed in the ovary (http://www.genecards.org/). It is a signalling molecule that resembles neuropeptides and acts by binding to alpha-neurexins and possibly other receptors (12). However, it is not known whether rs4954956 directly affects NXPH2 gene expression: rs4954956 is not in a highly conserved region (http://genome.ucsc.edu/cgi-bin/hgGateway) and various bioinformatics tools such as SNAP (http://www.broad.mit.edu/mpg/snap) did not reveal any variants associated with this SNP (r2 > 0.8) that have a putative function. It is possible that the functional effects of rs4954956 are due to other, as yet unidentified variants that are strongly correlated with this SNP.

Although breast cancer and ovarian cancer are both hormonal related female cancers and share some common genetic risk factors such as BRCA1/BRCA2 mutations, our data suggest that the overlap between ovarian and breast cancer susceptibility alleles is limited. We found no evidence of an association of ovarian cancer risk with the remaining 10 breast GWAS hits tested. Gates et al. (13) recently reported a null association for seven breast cancer susceptibility alleles in two ovarian cancer case–control populations, five of these alleles (i.e. rs2981582, rs3803662, rs889312, rs3817198 and rs13281615) were also genotyped in our study. Our meta-analysis pooling Gates et al.'s data together with ours for these five SNPs appears to confirm that they are not associated with ovarian cancer risks (data not shown). By combined data from 14 ovarian cancer case–control studies (5876 cases/9273 controls of non-Hispanic origin), we were able to provide at least 90% power to detect a co-dominant allele with a minor allele frequency of 0.27 that confers a relative risk of 1.1 at a Type 1 error of 0.05. However, our power to detect the alleles associated with smaller ovarian cancer risks is low, and we cannot exclude the possibility that the alleles investigated are associated with smaller risks for ovarian cancer.

In conclusion, we have found that rs4954956 is associated with an increased ovarian cancer risk, but none of the top six confirmed breast cancer susceptibility variants tested is associated with ovarian cancer. Further work will be needed to determine a functional rational for rs4954956 or any correlated variants in causing ovarian cancer.

MATERIALS AND METHODS

Study subjects

Initial set

Six case–control studies contributed data to the initial (stage1) analysis, including three studies from Europe (SEA, MAL and UKO), two from the USA (STA and USC), and one from Australia (AUS). Table 4 provides details for each of the studies which followed a population-based case–control design. Participation rates for cases and controls were generally excellent, and included largely White non-Hispanic women. In total, stage 1 comprised 2927 invasive ovarian cases and 4143 controls.

Table 4.
Study description

Validation set

Eight case–control studies contributed data to the validation set (stage 2) which included six studies from the USA (DOV, HOP, MAY, NCO, UCI, HAW) and two studies from Europe (GER, POL). Stage 2 also included additional samples from Australia (AUS) (413 cases/448controls), UKO (180 cases/333 controls) and USC (323 cases/479 controls) studies described above. Thus, the total validation set comprised 4060 invasive ovarian cancer cases and 6308 controls (Table 4).

These 14 case–control studies from OCAC contained a total of 6987 invasive ovarian cancer cases and 10451 controls when combined. Details of all these case–control studies have been published (14,15).

To reduce the possibility of population stratification, the analyses were limited to the 5876 cases and 9273 controls who were of non-Hispanic White origin for whom genotype information was available. All studies were approved by the review boards and Ethics Committees of their parent institutions and written informed consent was obtained from all participants.

Genotyping

Genotyping was performed at 11 different centres in 384-well plate formats and all but one study (AUS) used Taqman™ 7900HT Sequence Detection System according to the manufacturer's instructions. The AUS study used iPlex technology (Sequenom) for genotyping according to the manufacturer's instructions. Genotypes were determined using Allelic Discrimination Sequence Detection Software (Applied Biosystems, Warrington, UK). Assays were carried out in 384-well plates and included at least 3% duplicate samples in each plate for quality control. The six studies in the initial set were genotyped either at the Department of Oncology, University of Cambridge or at the Gynaecological Oncology Unit, University College London. The validation studies were genotyped by the individual study centres. Each assay was carried out using 10 ng DNA in a 5 µl reaction using TaqMan universal PCR master mix, forward and reverse primers and FAM and VIC labelled probes designed by Applied Biosystems (ABI Assay-by-design). Details of primer and probe sequences and assay conditions used for each polymorphism analysed are available upon request.

Genotyping quality control

We compared genotype call rates and concordance by study and overall. We used the following criteria as a measure of acceptable genotyping: (1) >3% sample duplicates included; (2) concordance rate for the duplicates ≥98%; (3) overall call rate (by study) >95% and (4) call rates >90% for each individual 384-well plate. The data for any SNP failing these criteria in any study were excluded from the final analyses. The HWE among White non-Hispanic controls was used to examine the quality of genotyping. For any SNP that was out of HWE (P < 0.05), the genotyping call rate was reviewed and the data excluded if the genotype clusters was found to be suboptimal. However, some studies with genotypes out of HWE were included if their genotypes, based on clusters, were of excellent quality. Genotyping consistency across labs was also evaluated by genotyping a common panel of CEPH-Utah trios including 90 individual DNA samples, five duplicate samples and one negative control (http://ccr.coriell.org/Sections/Search/Panel_Detail.aspx?PgId=202&Ref=HAPMAPPT01). The concordance of genotyping results between the centres was required to be >98% in order for the genotype data to be included. No attempt was made to repeat genotyping in DNA samples that did not provide a clear genotype at the first attempt resulting in variations in the number of studies/samples that were successfully genotyped for each polymorphism.

Statistics

Deviation of genotype frequencies from those expected under HWE was assessed by χ2 tests with one degree of freedom (1df) for each study of controls as part of the genotyping quality control. The primary test of association was the comparison of genotype frequencies in cases and controls using a test for gene-dose effect for each SNP through an interval variable with three levels: 0, 1, 2; one assigned to each genotype. This was done using unconditional logistic regression stratified by study. OR for allele dosage and associated 95% CI were also estimated by unconditional logistic regression. We tested for heterogeneity between study strata by comparing logistic regression models with and without a genotype-stratum interaction term using likelihood ratio tests. A subgroup analysis was used to compare genotype-specific risks by disease subgroup with the controls. We limited subgroup analysis to the serous histology type as the number of cases diagnosed with other subtypes was low.

FUNDING

This work was supported by Cancer Research UK, The Roswell Park Alliance, The Danish Cancer Society and The National Cancer Institute (CA71766, CA16056, RO1 CA61107 and RO1 CA122443). We are grateful to the family and friends of Kathryn Sladek Smith for their generous support of OCAC through their donations to the Ovarian Cancer Research Fund. US Army Medical Research and Materiel Command under DAMD17-01-1-0729, the Cancer Council Tasmania and Cancer Foundation of Western Australia (AOCS study), The National Health and Medical Research Council of Australia (199600) (ACS study). DOV study was supported by the US national Cancer Institute grants R01 CA87538 and R01 CA112523. The German Ovarian Cancer Study was supported by the German Federal Ministry of Education and Research of Germany, Programme of Clinical Biomedical Research grant 01 GB 9401 and the genotyping in part by the state of Baden-Württemberg through Medical Faculty of the University of Ulm (P.685). HAW study was supported by US Public Health Service grant R01-CA- 58598, and contracts N01-CN-55424 and N01-PC-35137 from the National Cancer Institute, NIH, Department of Health and Human Services. UCI study was supported by the National Institutes of Health, National Cancer Institute grants CA-58860, CA-92044 and the Lon V Smith Foundation grant LVS-39420. The UKOPS study is funded by the OAK Foundation. Some of this work was undertaken at UCLH/UCL who received a proportion of funding from the Department of Health's NIHR Biomedical Research Centre funding scheme.

Supplementary Material

[Supplementary Data]

ACKNOWLEDGEMENT

We thank all the individuals who took part in this study. We thank: Hannah Munday, Barbara Perkins, Mitul Shah, Clare Jordan, Judy West, Anabel Simpson, Sue Irvine, Anne Stafford, the SEARCH team: the local general practices and nurses and the Eastern Cancer Registry for recruitment of the UK cases and the EPIC-Norfolk investigators for recruitment of the UK controls for SEA study; We thank all members of the research team, including research nurses, research scientists, data entry personnel and consultant gynaecological oncologists for their help in establishing the UKOPS case–control collection, particularly we thank Eva Wozniak for the efforts in genotyping and Andy Ryan and Jeremy Ford for data and sample management for UKO study. We thank Ursula Eilber and Tanja Koehler for competent technical assistance for German Ovarian Cancer study. The AOCS Management Group (D. Bowtell, G. Chenevix-Trench, A. deFazio, D. Gertig, A. Green, P. Webb) gratefully acknowledges the contribution of all the clinical and scientific collaborators (see http://www.aocstudy.org/). The AOCS and ACS Management Group (A. Green, P. Parsons, N. Hayward, P. Webb, D. Whiteman) thank all of the project staff, collaborating institutions and study participants. We thank all the funding agencies and institutions for their support. D.F.E. is a Principal Research Fellow of Cancer Research UK, P.D.P.P. is CRUK Senior Clinical Research Fellow. S.J.R. is supported by the Mermaid/Eve Appeal, G.C.T. and P.W. are supported by the NHMRC.

Conflict of Interest statement. None declared.

REFERENCES

1. The Wellcome Trust Case Control Consortium. Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 2007;447:661–678. [PMC free article] [PubMed]
2. Eeles R.A., Kote-Jarai Z., Giles G.G., Olama A.A., Guy M., Jugurnauth S.K., Mulholland S., Leongamornlert D.A., Edwards S.M., Morrison J., et al. Multiple newly identified loci associated with prostate cancer susceptibility. Nat. Genet. 2008;40:316–321. [PubMed]
3. Gudmundsson J., Sulem P., Steinthorsdottir V., Bergthorsson J.T., Thorleifsson G., Manolescu A., Rafnar T., Gudbjartsson D., Agnarsson B.A., Baker A., et al. Two variants on chromosome 17 confer prostate cancer risk, and the one in TCF2 protects against type 2 diabetes. Nat. Genet. 2007;39:977–983. [PubMed]
4. Easton D.F., Eeles R.A. Genome-wide association studies in cancer. Hum. Mol. Genet. 2008;17:R109–R115. [PubMed]
5. Pharoah P.D., Ponder B.A. The genetics of ovarian cancer. Best Pract. Res. Clin. Obstet. Gynaecol. 2002;16:449–468. [PubMed]
6. Ghoussaini M., Song H., Koessler T., Al Olama A.A., Kote-Jarai Z., Driver K.E., Pooley K.A., Ramus S.J., Kjaer S.K., Hogdall E., et al. Multiple Loci with different cancer specificities within the 8q24 gene desert. J. Natl Cancer Inst. 2008;100:962–966. [PMC free article] [PubMed]
7. Easton D.F., Pooley K.A., Dunning A.M., Pharoah P.D., Thompson D., Ballinger D.G., Struewing J.P., Morrison J., Field H., Luben R., et al. Genome-wide association study identifies novel breast cancer susceptibility loci. Nature. 2007;447:1087–1093. [PMC free article] [PubMed]
8. Ioannidis J.P., Ntzani E.E., Trikalinos T.A., Contopoulos-Ioannidis D.G. Replication validity of genetic association studies. Nat. Genet. 2001;29:306–309. [PubMed]
9. Lohmueller K.E., Pearce C.L., Pike M., Lander E.S., Hirschhorn J.N. Meta-analysis of genetic association studies supports a contribution of common variants to susceptibility to common disease. Nat. Genet. 2003;33:177–182. [PubMed]
10. Colhoun H.M., McKeigue P.M., Davey S.G. Problems of reporting genetic associations with complex outcomes. Lancet. 2003;361:865–872. [PubMed]
11. Wacholder S., Chanock S., Garcia-Closas M., El Ghormli L., Rothman N. Assessing the probability that a positive report is false: an approach for molecular epidemiology studies. J. Natl Cancer Inst. 2004;96:434–442. [PubMed]
12. Missler M., Sudhof T.C. Neurexophilins form a conserved family of neuropeptide-like glycoproteins. J. Neurosci. 1998;18:3630–3638. [PubMed]
13. Gates M.A., Tworoger S.S., Terry K.L., De Vivo I., Hunter D.J., Hankinson S.E., Cramer D.W. Breast cancer susceptibility alleles and ovarian cancer risk in 2 study populations. Int. J. Cancer. 2008;124:729–733. [PMC free article] [PubMed]
14. Gayther S.A., Song H., Ramus S.J., Kjaer S.K., Whittemore A.S., Quaye L., Tyrer J., Shadforth D., Hogdall E., Hogdall C., et al. Tagging single nucleotide polymorphisms in cell cycle control genes and susceptibility to invasive epithelial ovarian cancer. Cancer Res. 2007;67:3027–3035. [PubMed]
15. Ramus S.J., Vierkant R.A., Johnatty S.E., Pike M.C., Van Den Berg D.J., Wu A.H., Pearce C.L., Menon U., Gentry-Maharaj A., Gayther S.A., et al. Consortium analysis of 7 candidate SNPs for ovarian cancer. Int. J. Cancer. 2008;123:380–388. [PMC free article] [PubMed]
16. Rossing M.A., Cushing-Haugen K.L., Wicklund K.G., Doherty J.A., Weiss N.S. Menopausal hormone therapy and risk of epithelial ovarian cancer. Cancer Epidemiol. Biomarkers Prev. 2007;16:2548–2556. [PubMed]

Articles from Human Molecular Genetics are provided here courtesy of Oxford University Press