Polymorphisms in NF-κB Inhibitors and Risk of Epithelial Ovarian Cancer
© White et al; licensee BioMed Central Ltd. 2009
Received: 19 December 2008
Accepted: 06 June 2009
Published: 06 June 2009
The nuclear factor-κB (NF-κB) family is a set of transcription factors with key roles in the induction of the inflammatory response and may be the link between inflammation and cancer development. This pathway has been shown to influence ovarian epithelial tissue repair. Inhibitors of κB (IκB) prevent NF-κB activation by sequestering NF-κB proteins in the cytoplasm until IκB proteins are phosphorylated and degraded.
We used a case-control study to evaluate the association between single nucleotide polymorphisms (SNPs) in NFKBIA and NFKBIB (the genes encoding IκBα and IκBβ, respectively) and risk of epithelial ovarian cancer. We queried 19 tagSNPs and putative-functional SNPs among 930 epithelial ovarian cancer cases and 1,037 controls from two studies.
The minor allele for one synonymous SNP in NFKBIA, rs1957106, was associated with decreased risk (p = 0.03).
Considering the number of single-SNP tests performed and null gene-level results, we conclude that NFKBIA and NFKBIB are not likely to harbor ovarian cancer risk alleles. Due to its biological significance in ovarian cancer, additional genes encoding NF-κB subunits, activating and inhibiting molecules, and signaling molecules warrant interrogation.
Despite estimates of more than 21,000 newly diagnosed cases of ovarian cancer and 15,000 related deaths each year in the United States , the etiology of ovarian cancer remains poorly understood. Known risk factors include increased risk with family history and use of fertility drugs, and decreased risk with oral contraceptive use, parity, and long duration of breast feeding . Rare, high-penetrant mutations in BRCA1 and BRCA2 account for approximately 40% of familial risk, leaving most inherited risk unexplained [3, 4]. The search for additional loci includes thoughtful selection of candidate genes in key biological pathways, an approach which has been successful in identifying new risk alleles for a variety of cancers .
Inflammation has been implicated in ovarian carcinogenesis because of its role in ovulation and post-ovulatory repair. During ovulation the ovarian epithelial surface is damaged, requiring a repair process involving the recruitment of leukocytes and inflammatory cytokines, release of nitrous oxide, DNA repair, and tissue restructuring [6–9]. Over time, this continuous repair of the ovarian epithelial tissue increases the likelihood of errors during replication, potentially leading to carcinogenesis. Nuclear factor-κB (NF-κB) refers to a family of "fast-acting" transcription factors that play a critical role in the inflammatory and innate immune responses . Stimulation by pro-inflammatory cytokines leads to the activation of NF-κB complexes which regulate the expression of key genes controlling apoptosis, angiogenesis, and cell proliferation [10–13]. Aberrant NF-κB functioning can lead to inhibition of apoptosis, constitutive cell replication, and increased angiogenesis, all of which are present in cancer cells . In ovarian cancer, several reports demonstrate the complex relationship between the immune system and established disease, suggesting a role for NF-κB. Immune effectors are thought to assist tumor growth; immunosuppressive regulatory T cells are associated with reduced survival, and the balance of the T cell subsets (regulated by NF-κB) has been shown to be critical to disease outcome . In addition, ovarian tumors acquire aberrant NF-κB functions allowing them to circumvent apoptotic pathways, specifically tumor necrosis factor alpha- (TNFα)-induced apoptosis, and afford protection against environmental insults such as anti-tumor immune effectors or chemotherapy [16–19].
Inhibitors of κB (IκB), IκBα, IκBβ, and IκBε, modulate NF-κB transcription by sequestering complexes of the NF-κB subunits (NF-κB1 [p50/p105], NF-κB2 [p52/p100], RelA [p65], RelB, and c-Rel) in the cytoplasm [10, 20]. In response to stimulation by TNFα, interleukin-1 (IL-1), and toll-like receptor (TLR) and T cell receptor (TCR) ligands, IκB proteins are phosphorylated by IκB kinase (IKK) complexes and degraded by the 26S proteasome, allowing for the release and nuclear localization of NF-κB proteins [11, 12, 21, 22]. Improper functioning of IκB proteins can lead to inhibition or constitutive activation of NF-κB . Because of NF-κB's central role in numerous cancer-related processes and involvement in risk of others cancers [23–26], we hypothesized that inherited variation in the genes encoding the key inhibitors IκBα and IκBβ (NFKBIA and NFKBIB, respectively) is associated with ovarian cancer risk. To examine this hypothesis, we assessed informative single-nucleotide polymorphisms (SNPs) in two case-control study populations.
Participants were recruited at Mayo Clinic in Rochester, MN and at Duke University in Durham, NC. Study protocols were approved by the Mayo Clinic and Duke University Institutional Review Boards, and all study participants provided informed consent. At Mayo Clinic, cases were women over age 20 years with histologically-confirmed epithelial ovarian cancer living in the Upper Midwest and enrolled within one year of diagnosis. Controls without ovarian cancer and without double oophorectomy were recruited from women seen for general medical examinations and frequency-matched to cases on age and region of residence. At Duke University, cases were women between age 20 and 74 years with histologically-confirmed primary epithelial ovarian cancer identified using the North Carolina Central Cancer Registry's rapid case ascertainment system within a 48-county region. Controls without ovarian cancer and who had at least one intact ovary were identified from the same region as the cases using list-assisted random digit dialing and frequency-matched to cases on race and age. Women with borderline and invasive disease were included; cases were 60% serous, 10% mucinous, 14% endometriod, 6% clear cell, and 9% multiple or other histologies. Additional participant details are provided elsewhere .
Data and Biospecimen Collection
Information on known and suspected risk factors were collected through in-person interviews at both sites using similar questionnaires. Mayo Clinic participants had an extra vial of blood drawn during their scheduled medical visit, and DNA was extracted from 10 to 15 mL fresh peripheral blood using the Gentra AutoPure LS Purgene salting out methodology (Gentra, Minneapolis, MN). Duke University participants had venipuncture performed at the conclusion of their interview. DNA samples were transferred to Mayo Clinic and, because of the relatively low quantities of DNA, they were whole-genome amplified (WGA) with the REPLI-G protocol (Qiagen Inc, Valencia CA) which we have shown to yield highly reproducible results with these samples . Genomic and WGA DNA concentrations were adjusted to 50 ηg/μl before genotyping and verified using PicoGreen dsDNA Quantitation kit (Molecular Probes, Inc., Eugene OR). Samples were bar-coded to ensure accurate and reliable sample processing and storage.
The selection of informative tagSNPs from among a larger pool of available SNPs allows for maximal genomic coverage and reduced genotyping redundancy . We identified tagSNPs within five kb of NFKBIA (chromosome 14q13.2, RefSeq NM_020529.1) and NFKBIB (chromosome 19q13.2, RefSeq NM_002503.3) using the algorithm of ldSelect  to bin pairwise-correlated SNPs at r2 ≥ 0.80 with minor allele frequency (MAF) ≥ 0.05 among publicly-available European-American data from the National Heart, Lung, and Blood Institute's Program for Genomic Applications SeattleSNPs gene-resequencing effort . Within bins of SNPs in linkage disequilibrium (LD), tagSNPs with the maximum predicted likelihood of genotype success (Illumina-provided SNP_Score, San Diego, CA) were selected. Within each gene, we binned 26 SNPs resulting in 13 tagSNPs in NFKBIA and eight tagSNPs in NFKBIB; four singleton SNPs in NFKBIA and two singleton SNPs in NFKBIB failed conversion in development of the custom genotype panel and were excluded. The inclusion of additional SNPs with particular suspected functional relevance further increases coverage in a hypothesis-based manner at minimal increased cost; thus, we included all putative-functional SNPs (within 1 kb upstream, 5' UTR, 3' UTR, or non-synonymous) with MAF ≥ 0.05 identified in Ensembl version 34 and Illumina-provided SNP_Score > 0.6, resulting in one additional 3' UTR and three additional 5' upstream SNPs in NFKBIA. A total of 13 NFKBIA SNPs and six NFKBIB SNPs were genotyped (see Additional file 1).
Genotyping of 1,086 genomic and 1,282 WGA DNA samples (total = 2,368 including duplicates and laboratory controls) on 2,051 unique study participants was performed at Mayo Clinic using the Illumina GoldenGate™ BeadArray assay and BeadStudio software for automated genotype clustering and calling separately for genomic and WGA samples according to a standard protocol . A total of 1,536 SNPs in a variety of pathways were attempted (including NFKBIA and NFKBIB), and 57 SNPs failed (poor clustering or call rate < 95%). Of 2,051 participants genotyped, 10 were ineligible and excluded, and 74 samples failed (call rate < 90%). Additional quality control (QC) information on the overall panel is provided elsewhere . In NFKBIA and NFKBIB, 18 of the 19 SNPs were successfully genotyped in both study populations (call rates > 98.9%); NFKBIA rs3138050 was excluded for Duke University samples due to poor clustering. For genotype QC metrics see Additional file 1.
Distributions of demographic and clinical variables were compared across case status using chi-square tests and t-tests as appropriate. Individual SNP associations for ovarian cancer risk were assessed using logistic regression, in which odds ratios (ORs) and 95% confidence intervals (CIs) were estimated. Primary tests for associations assumed an ordinal (log-additive) effect with simple tests for trend, as well as separate comparisons of heterozygous and minor allele homozygous women to major allele homozygous women (referent) using a 2 degree-of-freedom (d.f.) test. In addition, we used a gene-centric principal components analysis to create orthogonal linear combinations of minor allele counts. The component linear combinations that accounted for at least 90% of the variability in the gene were included in a multivariable logistic regression model and simultaneously tested for gene-specific global significance using a likelihood ratio test. Haplotype frequencies were also estimated within each gene and a global haplotype score test of association between haplotypes and ovarian cancer risk was conducted at the gene level using a score test . Individual haplotype tests compared each haplotype to all other haplotypes combined. NFKBIA rs3138050 was excluded from gene-level analyses due to failed genotyping in Duke University participants. All analyses were adjusted for age, race, region of residence, body mass index, hormone therapy use, oral contraceptive use, parity, and age at first birth. We used SAS (SAS Institute, Cary, NC, Version 8, 1999), Haplo.stats http://mayoresearch.mayo.edu/mayo/research/biostat/schaid.cfm, and S-Plus (Insightful Corp, Seattle, WA, Version 7.05, 2005) software systems.
Selected Characteristics of Study Participants
Cases (N = 396)
Controls (N = 469)
Cases (N = 534)
Controls (N = 568)
Mean (S.D.) yrs
Body mass index
< 23 kg/m2
≥ 29 kg/m2
Age at menarche
< 12 yrs
≥ 14 yrs
Oral contraceptive use
≥ 48 months
Postmenopausal hormone use
≥ 60 months
Parity, n/Age at first birth, yrs
1–2/≤ 20 yrs
1–2/> 20 yrs
≥ 3/≤ 20 yrs
≥ 3/> 20 yrs
Ovarian cancer family history
Ovarian or breast cancer family history
Smoking, pack years
High school diploma
Post high school
NFKBIA and NFKBIB Polymorphisms and Adjusted Risk of Epithelial Ovarian Cancer
General Model OR (95%CI)
Ordinal Model OR (95% CI)
bp to next
AB v AA
BB v AA
To assess whether overall variation within each gene was associated with ovarian cancer risk, we performed multiple logistic regression for participants with complete genotype data (N = 1,901 for NFKBIA, N = 1,930 for NFKBIB). Gene-level logistic regression revealed null results (NFKBIA, d.f. = 12, p = 0.23; NFKBIB, d.f. = 6, p = 0.97) as did the potentially more-powerful logistic regression analysis using principal components (NFKBIA, d.f. = 6, p = 0.79; NFKBIB, d.f. = 4, p = 0.89).
Haplotype analysis can reveal hidden associations with alleles at ungenotyped variants. Within NFKBIA, five haplotypes were estimated to have frequencies > 0.05; no associations were observed with any of these. Three rare haplotypes were associated with increased risk (see Additional file 2); however, overall variation among all haplotypes combined was not associated with risk (p = 0.32). Four NFKBIB haplotypes had estimated frequencies > 0.05; no common or rare haplotypes were associated with risk, and overall haplotype associations were null (p = 0.50). In summary, single-SNP, multi-SNP, and haplotype analyses do not indicate that NFKBIA or NFKBIB harbor risk alleles for ovarian cancer.
To our knowledge, this is the first examination of inherited variation in the NF-κB signaling pathway in relation to epithelial ovarian cancer risk. The two genes studied, NFKBIA and NFKBIB, encode IκBs with critical roles in regulating NF-κB transcription by directly binding to NF-κB subunits in the cytoplasm. We assessed a comprehensive set of SNPs in these two genes in a large combined case-control study, and found no evidence of association. Strengths of this study include large sample size, choice of candidate genes, use of multiple study populations, LD-based SNP selection, robust genotyping, control of potential confounding variables, and application of a variety of genetic analysis tools. Limitations of this study include the focus on only two genes in a large pathway, the lack of an independent replication outside of the Mayo Clinic and Duke University datasets, and the lack of functional analyses. This study was designed to detect modest genetic associations with ovarian cancer risk; results suggest that common risk alleles of modest effect size may not reside in NFKBIA or NFKBIB.
Although no association was found here, inherited variation in NFKBIA and NFKBIB have been associated with increased risk of other cancers including melanoma , colorectal cancer , multiple myeloma , and Hodgkin lymphoma . Considering the vast evidence on the importance of NF-κB in carcinogenesis, additional examination of NF-κB including study of inherited variation in the NF-κB pathway and risk of epithelial ovarian cancer is warranted.
Study of inherited variation within the NF-κB pathway has the potential to identify risk alleles accounting for the residual increased familial risk of ovarian cancer . The present analysis is an early epidemiologic assessment which indicates that NFKBIA and NFKBIB are not likely to harbor risk alleles under our statistical assumptions; the key limitation of our study is its focus on only two genes. Other genes to examine are numerous and a more thorough examination of polymorphisms within this pathway is needed to better understand the complexities of ovarian carcinogenesis.
tumor necrosis factor alpha
inhibitors of κB
T cell receptor
minor allele frequency
Grant support: R01 CA88868, R01 CA122443, Fraternal Order of Eagles, Minnesota Ovarian Cancer Alliance.
- Jemal A, Siegel R, Ward E, Hao Y, Xu J, Murray T, Thun MJ: Cancer statistics, 2008. CA Cancer J Clin. 2008, 58 (2): 71-96. 10.3322/CA.2007.0010.View ArticlePubMedGoogle Scholar
- Whittemore AS, Harris R, Itnyre J: Characteristics relating to ovarian cancer risk: collaborative analysis of 12 US case-control studies. II. Invasive epithelial ovarian cancers in white women. Collaborative Ovarian Cancer Group. Am J Epidemiol. 1992, 136 (10): 1184-1203.PubMedGoogle Scholar
- Antoniou AC, Easton DF: Risk prediction models for familial breast cancer. Future oncology (London, England). 2006, 2 (2): 257-274. 10.2217/14796618.104.22.1687.View ArticleGoogle Scholar
- Risch HA, McLaughlin JR, Cole DE, Rosen B, Bradley L, Kwan E, Jack E, Vesprini DJ, Kuperstein G, Abrahamson JL, et al: Prevalence and penetrance of germline BRCA1 and BRCA2 mutations in a population series of 649 women with ovarian cancer. Am J Hum Genet. 2001, 68 (3): 700-710. 10.1086/318787.View ArticlePubMedPubMed CentralGoogle Scholar
- Dong LM, Potter JD, White E, Ulrich CM, Cardon LR, Peters U: Genetic susceptibility to cancer: the role of polymorphisms in candidate genes. Jama. 2008, 299 (20): 2423-2436. 10.1001/jama.299.20.2423.View ArticlePubMedPubMed CentralGoogle Scholar
- Fathalla MF: Incessant ovulation – a factor in ovarian neoplasia?. Lancet. 1971, 2 (7716): 163-10.1016/S0140-6736(71)92335-X.View ArticlePubMedGoogle Scholar
- Bonello N, McKie K, Jasper M, Andrew L, Ross N, Braybon E, Brannstrom M, Norman RJ: Inhibition of nitric oxide: effects on interleukin-1 beta-enhanced ovulation rate, steroid hormones, and ovarian leukocyte distribution at ovulation in the rat. Biology of reproduction. 1996, 54 (2): 436-445. 10.1095/biolreprod54.2.436.View ArticlePubMedGoogle Scholar
- Zackrisson U, Mikuni M, Wallin A, Delbro D, Hedin L, Brannstrom M: Cell-specific localization of nitric oxide synthases (NOS) in the rat ovary during follicular development, ovulation and luteal formation. Human reproduction (Oxford, England). 1996, 11 (12): 2667-2673.View ArticleGoogle Scholar
- Schildkraut JM, Bastos E, Berchuck A: Relationship between lifetime ovulatory cycles and overexpression of mutant p53 in epithelial ovarian cancer. J Natl Cancer Inst. 1997, 89 (13): 932-938. 10.1093/jnci/89.13.932.View ArticlePubMedGoogle Scholar
- Ahn KS, Sethi G, Aggarwal BB: Nuclear factor-kappa B: from clone to clinic. Current molecular medicine. 2007, 7 (7): 619-637. 10.2174/156652407782564363.View ArticlePubMedGoogle Scholar
- Sethi G, Sung B, Aggarwal BB: Nuclear factor-kappaB activation: from bench to bedside. Experimental biology and medicine (Maywood, NJ). 2008, 233 (1): 21-31. 10.3181/0707-MR-196.View ArticleGoogle Scholar
- Inoue J, Gohda J, Akiyama T, Semba K: NF-kappaB activation in development and progression of cancer. Cancer Sci. 2007, 98 (3): 268-274. 10.1111/j.1349-7006.2007.00389.x.View ArticlePubMedGoogle Scholar
- Li H, Lin X: Positive and negative signaling components involved in TNFalpha-induced NF-kappaB activation. Cytokine. 2008, 41 (1): 1-8. 10.1016/j.cyto.2007.09.016.View ArticlePubMedGoogle Scholar
- Naugler WE, Karin M: NF-kappaB and cancer-identifying targets and mechanisms. Curr Opin Genet Dev. 2008, 18 (1): 19-26. 10.1016/j.gde.2008.01.020.View ArticlePubMedPubMed CentralGoogle Scholar
- Curiel TJ, Coukos G, Zou L, Alvarez X, Cheng P, Mottram P, Evdemon-Hogan M, Conejo-Garcia JR, Zhang L, Burow M, et al: Specific recruitment of regulatory T cells in ovarian carcinoma fosters immune privilege and predicts reduced survival. Nat Med. 2004, 10 (9): 942-949. 10.1038/nm1093.View ArticlePubMedGoogle Scholar
- Xiao CW, Yan X, Li Y, Reddy SA, Tsang BK: Resistance of human ovarian cancer cells to tumor necrosis factor alpha is a consequence of nuclear factor kappaB-mediated induction of Fas-associated death domain-like interleukin-1beta-converting enzyme-like inhibitory protein. Endocrinology. 2003, 144 (2): 623-630. 10.1210/en.2001-211024.View ArticlePubMedGoogle Scholar
- Takeyama H, Wakamiya N, O'Hara C, Arthur K, Niloff J, Kufe D, Sakarai K, Spriggs D: Tumor necrosis factor expression by human ovarian carcinoma in vivo. Cancer Res. 1991, 51 (16): 4476-4480.PubMedGoogle Scholar
- Hassan MI, Kassim SK, Saeda L, Laban M, Khalifa A: Ovarian cancer-induced immunosuppression: relationship to tumor necrosis factor-alpha (TNF-alpha) release from ovarian tissue. Anticancer Res. 1999, 19 (6C): 5657-5662.PubMedGoogle Scholar
- Spriggs DR, Imamura K, Rodriguez C, Sariban E, Kufe DW: Tumor necrosis factor expression in human epithelial tumor cell lines. J Clin Invest. 1988, 81 (2): 455-460. 10.1172/JCI113341.View ArticlePubMedPubMed CentralGoogle Scholar
- Tergaonkar V, Correa RG, Ikawa M, Verma IM: Distinct roles of IkappaB proteins in regulating constitutive NF-kappaB activity. Nat Cell Biol. 2005, 7 (9): 921-923. 10.1038/ncb1296.View ArticlePubMedGoogle Scholar
- Weil R, Israel A: Deciphering the pathway from the TCR to NF-kappaB. Cell Death Differ. 2006, 13 (5): 826-833. 10.1038/sj.cdd.4401856.View ArticlePubMedGoogle Scholar
- Braun T, Carvalho G, Fabre C, Grosjean J, Fenaux P, Kroemer G: Targeting NF-kappaB in hematologic malignancies. Cell Death Differ. 2006, 13 (5): 748-758. 10.1038/sj.cdd.4401874.View ArticlePubMedGoogle Scholar
- Osborne J, Lake A, Alexander FE, Taylor GM, Jarrett RF: Germline mutations and polymorphisms in the NFKBIA gene in Hodgkin lymphoma. Int J Cancer. 2005, 116 (4): 646-651. 10.1002/ijc.21036.View ArticlePubMedGoogle Scholar
- Spink CF, Gray LC, Davies FE, Morgan GJ, Bidwell JL: Haplotypic structure across the I kappa B alpha gene (NFKBIA) and association with multiple myeloma. Cancer Lett. 2007, 246 (1–2): 92-99. 10.1016/j.canlet.2006.02.001.View ArticlePubMedGoogle Scholar
- Gao J, Pfeifer D, He LJ, Qiao F, Zhang Z, Arbman G, Wang ZL, Jia CR, Carstensen J, Sun XF: Association of NFKBIA polymorphism with colorectal cancer risk and prognosis in Swedish and Chinese populations. Scandinavian journal of gastroenterology. 2007, 42 (3): 345-350. 10.1080/00365520600880856.View ArticlePubMedGoogle Scholar
- Bu H, Rosdahl I, Sun XF, Zhang H: Importance of polymorphisms in NF-kappaB1 and NF-kappaBIalpha genes for melanoma risk, clinicopathological features and tumor progression in Swedish melanoma patients. J Cancer Res Clin Oncol. 2007, 133 (11): 859-866. 10.1007/s00432-007-0228-7.View ArticlePubMedGoogle Scholar
- Sellers TA, Schildkraut JM, Pankratz VS, Vierkant RA, Fredericksen ZS, Olson JE, Cunningham JM, Taylor W, Liebow M, McPherson CP, et al: Estrogen bioactivation, genetic polymorphisms, and ovarian cancer. Cancer Epidemiol Biomarkers Prev. 2005, 14 (11 Pt 1): 2536-2543. 10.1158/1055-9965.EPI-05-0142.View ArticlePubMedGoogle Scholar
- Cunningham JM, Sellers TA, Schildkraut JM, Fredericksen ZS, Vierkant RA, Kelemen LE, Gadre M, Phelan CM, Huang Y, Meyer JG, et al: Performance of amplified DNA in an Illumina GoldenGate BeadArray assay. Cancer Epidemiol Biomarkers Prev. 2008, 17 (7): 1781-1789. 10.1158/1055-9965.EPI-07-2849.View ArticlePubMedPubMed CentralGoogle Scholar
- Carlson CS, Eberle MA, Rieder MJ, Yi Q, Kruglyak L, Nickerson DA: Selecting a maximally informative set of single-nucleotide polymorphisms for association analyses using linkage disequilibrium. Am J Hum Genet. 2004, 74 (1): 106-120. 10.1086/381000.View ArticlePubMedGoogle Scholar
- Seattle SNPs: Variation Discovery Resource. [http://pga.mbt.washington.edu]
- Oliphant A, Barker DL, Stuelpnagel JR, Chee MS: BeadArray technology: enabling an accurate, cost-effective approach to high-throughput genotyping. Biotechniques. 2002, 56-58. Suppl
- Schaid DJ, Rowland CM, Tines DE, Jacobson RM, Poland GA: Score tests for association between traits and haplotypes when linkage phase is ambiguous. Am J Hum Genet. 2002, 70 (2): 425-434. 10.1086/338688.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/9/170/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.