Interaction among apoptosis-associated sequence variants and joint effects on aggressive prostate cancer

doi:10.1186/1755-8794-5-11

BMC Med Genomics. 2012; 5: 11.

Published online 2012 April 30. doi: 10.1186/1755-8794-5-11

PMCID: PMC3355002

Interaction among apoptosis-associated sequence variants and joint effects on aggressive prostate cancer

Nicole A Lavender,¹ Erica N Rogers,¹ Susan Yeyeodu,² James Rudd,² Ting Hu,⁵ Jie Zhang,⁴ Guy N Brock,⁴ Kevin S Kimbro,² Jason H Moore,³ David W Hein,¹ and La Creis R Kidd^1,⁵

¹Department of Pharmacology & Toxicology, School of Medicine, University of Louisville (UofL), Louisville, KY, USA

²Department of Biology, The Julius L. Chambers Biomedical/Biotechnology Research Institute, North Carolina Central University, Durham, NC, USA

³Departments of Genetics and Community and Family Medicine, Institute for Quantitative Biomedical Sciences, Norris Cotton Cancer Center, Dartmouth Medical School, Lebanon, NH, USA

⁴Department of Bioinformatics & Biostatistics, School of Public Health and Information Sciences, UofL, Louisville, KY, USA

⁵505 South Hancock Street, Clinical & Translational Research Building, Room 306, Louisville, KY 40202, USA

Corresponding author.

Nicole A Lavender: nlavender/at/tgen.org; Erica N Rogers: enroge01/at/louisville.edu; Susan Yeyeodu: syeyeodu/at/nccu.edu; James Rudd: jrudd/at/nccu.edu; Ting Hu: ting.hu/at/dartmouth.edu; Jie Zhang: j0Zhang/at/louisville.edu; Guy N Brock: g0broc01/at/lousiville.edu; Kevin S Kimbro: kkimbro/at/nccu.edu; Jason H Moore: Jason.H.Moore/at/dartmouth.edu; David W Hein: dwhein01/at/louisville.edu; La Creis R Kidd: lrkidd01/at/louisville.edu

Received May 10, 2011; Accepted April 30, 2012.

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.

This article has been cited by other articles in PMC.

Abstract

Background

Molecular and epidemiological evidence demonstrate that altered gene expression and single nucleotide polymorphisms in the apoptotic pathway are linked to many cancers. Yet, few studies emphasize the interaction of variant apoptotic genes and their joint modifying effects on prostate cancer (PCA) outcomes. An exhaustive assessment of all the possible two-, three- and four-way gene-gene interactions is computationally burdensome. This statistical conundrum stems from the prohibitive amount of data needed to account for multiple hypothesis testing.

Methods

To address this issue, we systematically prioritized and evaluated individual effects and complex interactions among 172 apoptotic SNPs in relation to PCA risk and aggressive disease (i.e., Gleason score ≥ 7 and tumor stages III/IV). Single and joint modifying effects on PCA outcomes among European-American men were analyzed using statistical epistasis networks coupled with multi-factor dimensionality reduction (SEN-guided MDR). The case-control study design included 1,175 incident PCA cases and 1,111 controls from the prostate, lung, colo-rectal, and ovarian (PLCO) cancer screening trial. Moreover, a subset analysis of PCA cases consisted of 688 aggressive and 488 non-aggressive PCA cases. SNP profiles were obtained using the NCI Cancer Genetic Markers of Susceptibility (CGEMS) data portal. Main effects were assessed using logistic regression (LR) models. Prior to modeling interactions, SEN was used to pre-process our genetic data. SEN used network science to reduce our analysis from > 36 million to < 13,000 SNP interactions. Interactions were visualized, evaluated, and validated using entropy-based MDR. All parametric and non-parametric models were adjusted for age, family history of PCA, and multiple hypothesis testing.

Results

Following LR modeling, eleven and thirteen sequence variants were associated with PCA risk and aggressive disease, respectively. However, none of these markers remained significant after we adjusted for multiple comparisons. Nevertheless, we detected a modest synergistic interaction between AKT3 rs2125230-PRKCQ rs571715 and disease aggressiveness using SEN-guided MDR (p = 0.011).

Conclusions

In summary, entropy-based SEN-guided MDR facilitated the logical prioritization and evaluation of apoptotic SNPs in relation to aggressive PCA. The suggestive interaction between AKT3-PRKCQ and aggressive PCA requires further validation using independent observational studies.

Keywords: Prostate cancer, Apoptosis, Single nucleotide polymorphisms, Gene-gene interactions, Multifactor dimensionality reduction (MDR), Statistical epistasis networks (SEN)

Background

Prostate cancer (PCA) is the most frequently diagnosed cancer and the 2^ndleading cause of cancer-related deaths among men in the United States [1]. The American Cancer Society estimates that 26-29% of all new cancer cases and cancer-related deaths are attributed to PCA cancer. Well established PCA risk factors include older age, black race, and family history of PCA. However, other potential contributors of this disease may include lifestyle and genetic factors as well as imbalances within important biological pathways.

Apoptosis or programmed cell death is one such biological process that moderates cell differentiation, proliferation, death, whole body homeostasis and tumorigenesis [2-4]. This process is controlled by cell death (e.g. BAD, CASP, BIK) and cell survival proteins (e.g. BCL2, NFκB, AKT3) that induce or block apoptosis, respectively, as summarized in Table Table11[2,3,5]. Decreased apoptotic cell death and increased cell proliferation may lead to clonal expansion and tumor growth [2]. Failure to undergo apoptosis permits survival of transformed cells, leading to subsequent genetic alteration, genomic instability, and ultimately a more invasive cancer phenotype [3]. Imbalances in apoptosis-related genes may play an important role in PCA susceptibility as well as disease progression. For example, several independent studies have shown that overexpression of cell survival indicators (e.g., BCL-2, CARD8, IKBKE, PRKCQ, and PIK3CB, AKT3) or down-regulation of cell death markers (e.g., BCL2L14) are associated with more aggressive phenotypes, higher Gleason grade, increased tumor progression, and poor PCA prognosis [6-19].

Table 1

Selected genes involved in the regulation of apoptosis

There is mounting epidemiological evidence that genetic alterations in apoptosis-related genes play an important role in tumorigenesis. Apoptosis-associated sequence variants, when considered individually, may minimally influence the risk of developing numerous cancers, such as multiple myeloma, squamous cell carcinoma, chronic lymphocytic leukemia, non-Hodgkin lymphoma, colorectal, ovarian, breast, pancreatic, and non-small cell lung [20-37]. However, the impact of individual apoptosis-related single nucleotide polymorphisms (SNPs) and their interactions on PCA outcomes remains understudied.

Genome wide association studies (GWAS), involving the evaluation of millions of SNPs within various biological pathways, has resulted in the detection of numerous PCA susceptibility loci [38]. However, most GWAS and PCA epidemiology studies place emphasis on individual SNP effects. Consequently, researchers often ignore the fact that complex diseases such as PCA are governed by complex gene-gene and gene-environment interactions within distinct biological pathways. Consequently, we sought to evaluate millions of interactions among apoptosis SNPs and their joint modifying impact on PCA risk and disease progression.

This report focuses on the impact of complex interactions among 172 apoptosis-related SNPs on PCA among 2,286 European-American male participants of the Cancer Genetic Markers of Susceptibility Study (CGEMS). The exhaustive assessment of all possible two-, three-, and four-way interactions among 172 SNPs, totaling more than 36 million SNP combinations is computationally burdensome. The statistical challenge stems from the prohibitive amount of data needed for multiple hypothesis testing. To address this issue, for the first time, we use statistical epistasis networks (SEN) guided multi-factor dimensionality reduction (MDR) to efficiently pre-process our genetic data, prior to modeling higher order interactions. As previously reported by Hu and co-workers (2011), SEN uses network science to generate a large connected cluster of pairwise interactions embedded within a genetic dataset [39]. The resulting network, consisting of a subset of high susceptibility SNPs, is used to guide SNP interactions using MDR. MDR is a rigorous statistical tool designed, in part, to evaluate main effects and complex interactions in relation to discrete outcomes [40-48]. This report may serve as a foundation for researchers interested in a state-of the art bioinformatics technique, namely SEN-guided MDR, which facilitates the logical prioritization of SNPs for gene-gene interaction analyses in relation to PCA susceptibility and prognosis.

Results

Cancer Genetic Markers of Susceptibility (CGEMS) Population Description

CGEMS study participants consisted of middle aged non-Hispanic Caucasian men, ranging in age from 55 to 81. Relative to control subjects, men diagnosed with PCA were more likely to have a family history of the disease (11.4% versus 6.3%), PSA levels ≥ 4 ng/ml (48.5% versus 6.5%), and abnormal DRE exams (60.7% versus 46.2%) [data not shown]. Men with aggressive PCA had both a Gleason score ≥ 7 and tumor stage III/IV. Although non-aggressive PCA cases primarily had a Gleason score < 7 and tumor stage I/II, a small percentage (1.8%) of them had a Gleason score > 6 and a tumor stage I/II. Among cases, there were no significant differences in family history of disease or PSA levels between men diagnosed comparing men with aggressive and non-aggressive disease (data not shown). The vast majority of subjects with aggressive PCA had higher Gleason scores than subjects with non-aggressive PCA, although 9 of the 488 non-aggressive cases (1.08%) had a Gleason score greater than 6.

Single SNP Effects

We examined 172 apoptosis-related sequence variants in relation to PCA risk and disease aggressiveness, as shown in Tables Tables22 and and3.3. Age and family history of PCA did not modify the relationship between the sequence variants and the risk of developing PCA or aggressive disease. Inheritance of one or more TNFRSF10B rs1001793 or TNFSF10 rs4894559 A alleles was nominally associated with a 1.23-1.38 fold increase in PCA risk (p ≤ 0.039), as shown in Table Table2.2. Whereas, eight apoptosis-related SNPs were associated with a moderate 18-50% reduction in disease susceptibility [BIK rs4988366 AG & AG + GG; BNIP3L rs10503786 TT & TC + TT; CASP9 rs1052571 TT; IKBKE rs1539243 TC & TC + TT; PRKCE rs608139 CC; PRKCE rs935673 GG & AG + GG; TNFRSF1A rs4149576 AG, GG, & AG + GG; and TNFRSF1A rs4149577 TC & TC + CC] (p ≤ 0.043). However, after applying Bonferonni correction, none of these findings remained statistically significant under the per-allele, dominant and recessive genetic models.

Table 2

Association of apoptosis SNPs with prostate cancer

Table 3

Association of apoptosis SNPs with aggressive prostate cancer

In terms of disease progression, eleven apoptosis SNPs were modestly associated with a 1.28-1.81 fold increase in the risk of developing aggressive PCA [AKT3 rs10157763 TC & TC + TT; AKT3 rs10927067 AG & AG + AA; AKT3 rs12031994 TC & TC + CC;AKT3 rs2125230 AG & AG + AA; AKT3 rs2125231 AG + AA; AKT3 rs2345994 TC & TC + CC; AKT3 rs4132509 AC & AC + AA; BCL2L14 rs885720 AG & AG + AA; CARD8 rs10405717 TC & TC + TT; IKBKE rs11578093 GG & TG + GG; and RAF1 rs13060691 TG] (p ≤ 0.031), as summarized in Table Table3.3. In addition, possession of IKBKE rs1539243 TT and PIK3CB rs500687 CC genotypes were moderately linked with a 44-56% reduction in aggressive PCA susceptibility (p ≤ 0.019). However, none of these loci remained significant after adjusting for multiple comparison bias.

Combined Gene Effects

We investigated the ability of gene-gene interactions to predict PCA outcomes using statistical epistasis networks (SEN) guided multifactor dimensionality reduction (MDR), information gain (IG) measure, and hierarchical interactions graphs (hIG). SEN was used to generate a network to visualize the genetic architecture of PCA outcomes, prior to an exhaustive search of SNP interactions. This approach uses information theory to build an epistasis network from the strongest pairwise SNP interactions. A network, consisting of vertices (i.e., SNPs) joined by edges (i.e., SNP pairs), enabled us to limit our MDR analysis to a subset of the most informative SNPs in relation to PCA outcomes. We were unable to build a statistically significant epistasis network for PCA susceptibility. However, the network for aggressive PCA was topographically significant (p < 0.05), given the results from 1000 permutated data sets. This network consisted of 91 vertices (i.e., individual SNP effects), 80 edges (i.e., pairwise interactions), and 18 total components (i.e., "sub-networks" of vertices and edges). The largest connected component of the aggressive PCA SEN diagram involved 24 vertices and 34 edges, as shown in Figure Figure11 andAdditional file 1. Consequently, MDR analysis was limited to < 13,000 interactions among 24 SNPs, instead of more than 36 million interactions for 172 SNPs. SEN-guided MDR revealed complex interactions between AKT3 rs12031994-PRKCQ rs571715 as well as AKT3 rs12031994-BID rs366542-PRKCQ rs571715 that were significantly associated with disease aggressiveness (p ≤ 0.009). However, the network and entropy graphs indicated these findings were mainly driven by AKT3 rs12031994. This is evident by the marker's large vertex size (i.e., large main effect) in the SEN graph (Figure (Figure1)1) and a lack of a strong pairwise interaction between AKT3 rs12031994-PRKCQ rs571715 in the entropy graph. The AKT3-PRKCQ interaction resulted in joint mutual information (I) of only 0.54% relative to the mutual information (I) of 0.87% for AKT3 alone. Therefore, prior to repeating MDR, we removed the AKT3rs12031994 marker. This resulted in the detection of a modest synergistic interaction between AKT3 rs2125230-PRKCQ rs571715 (p = 0.011), as shown in Table Table4.4. Based on MDR, this model was the best two-factor model. It was selected 9 out of 10 times by MDR based on a cross-validation consistency (CVC) score of 90% and predictive accuracy of 56.3%. Furthermore, we observed a higher mutual information for the combined effect of AKT3 rs2125230-PRKCQ rs571715 (I = 0.66%), in contrast to AKT3 rs2125230 alone (I = 0.45%) or PRKCQrs571715 alone (I = 0.21%), as revealed in Figure Figure22.

Figure 1

Visualization of the interaction among apoptosis-related sequence variants using Statistical Epistasis Network Modeling (SEN)-guided MDR. The Statistical Epistasis Network (SEN) was generated using a pairwise interaction strength cut-off or 1%. Overall, (more ...)

Table 4

Multi-dimensionality reduction models for apoptosis-related markers in relation to aggressive prostate cancer (PCA)

Figure 2

Interaction Entropy model. This graphical model, describes the percent entropy that is explained by each apoptosis-related SNP or pairwise combination within our study population. Positive percent entropy indicates information gain (IG) or synergy; whereas, (more ...)

To confirm the validity of the SEN results, MDR analysis was also conducted on all 172 SNPs (Figure (Figure11 andAdditional file 1) as well as 148 loci that were not included in the main component of the network. As expected, we did not detect any significant two-, three-, or four- way interaction models (p ≥ 0.184).

Discussion

Previous genome-wide association studies (GWAS) have identified linkages about 40 PCA loci, including between PCA genetic alterations detected in the 8q24 region, β-microseminoprotein (MSMB), and allele -8 of the microsatellite DG8S737 [38]. These studies have limited their scope to individual SNPs across the entire genome. Such assessments tend to ignore the genetic architecture of PCA that potentially involves complex interactions along key regulatory pathways. For the first time, the current study evaluates complex interactions among 172 apoptosis-related SNPs in relation to PCA risk and disease aggressiveness among 2,286 European-American men using SEN-guided MDR. Specifically, SEN was used to build a topographically significant aggressive PCA epistasis network, prior to evaluating complex interactions. This inferred epistasis network consisted of 24 SNPs and 34 pairwise interactions, and reduced MDR analysis from > 36 million to < 13,000 SNP interactions. Consequently, we observed a non-linear and modest interaction between AKT3 rs2125230-PRKCQ rs571715 in relation to aggressive PCA. This state-of-the-art bioinformatics technique facilitates the logical prioritization of SNPs for gene-gene interaction analyses in relation to complex diseases.

Unfortunately, there are no published reports on the functional consequence of these two intronic SNPs in AKT3 and PRKCQ in relation to mRNA stability/expression protein expression/structure/function or PCA outcomes. However, we speculate that the AKT3 rs2125230 and PRKCQ rs571715 sequence variants, with minor allele frequencies ranging from 14.4-22.9 among men of European descent, may alter transcription regulation, leading to increased mRNA expression. Increased mRNA/protein expression AKT3 rs2125230 and PRKCQ rs571715 may cause: decreased apoptosis, an escape of transformed cells from programmed cell death, increased accumulation of genetic alterations, genomic instability, and ultimately an invasive PCA phenotype. Thus, in vitro and in vivo assays using (short hairpin RNAs) shRNAs or small interfering RNAs (siRNAs) are needed to elucidate the impact of AKT3-PRKCQ genetic alterations on protein expression, apoptosis capacity, and prostate tumorigenesis.

The impact of a non-linear interaction along the AKT3 rs2125230-PRKCQ rs571715 axis in relation to aggressive PCA may be attributed to markers involved in the apoptosis signaling pathway. Overexpression of PRKCQ and AKT3 are associated with invasive cancer phenotypes [6-19]. In fact, AKT3 is responds to insulin and growth factors, transduces signals including cell death, and is upregulated in androgen-independent PCA cell lines [49]. PRKCQ is also associated with apoptosis. In particular PRKCQ, a protein kinase C (PRKC) family member, promotes cell survival by inactivating BAD (BCL2-associated agonist of cell death), which subsequently results in NFκB activation.

Although AKT3 and PRKCQ are involved in pro-survival pathways, their interaction is not fully understood. However, their interactions with other related protein kinases may offer biological clues on the mechanism by which AKT3 and PRKCQ synergistically influence aggressive PCA. For example, PRKCQ interacts with another AKT family member, AKT1, to activate NFκB [50]. If the AKT3-PRKCQ axis has a similar function as other protein kinases, namely AKT1 and PRKCQ, then these pro-survival markers may synergistically activate NFκB. As a result, activated NFκB may enable the tumor to escape programmed cell death and progress toward an aggressive PCA phenotype.

Numerous observational studies evaluated the impact of one or more apoptosis-related SNPs in relation to cancer outcomes [20-37]. Among 15 case-control studies, less than 10% of the sequence variants evaluated in the current study were significantly associated with various tumors, including of the colon, rectal, ovarian, breast, pancreatic, and non-small cell lung cancers [20-37]. In particular, modest cancer risk estimates were observed among 8 apoptosis-related SNPs detected in CASP3, CASP8, CASP9, TP53, NFKB2, and NFKBIA. However, some of these studies were limited by a small sample size, small number of analyzed SNPs, or failure to consider the impact of multiple SNPs on disease susceptibility. In the current study, there were 24 SNPs detected in 10 apoptosis-related genes [i.e., AKT3, BIK, BNIP3L, CARD8, CASP9, IKBKE, PRKCE, TNFSF10, and TNFRSF10 (B, D)]. These apoptosis-related SNPs, after adjusting for confounders, were modestly associated with PCA risk and/or aggressive disease. Yet, these findings lost statistical significance after adjusting for multiple hypothesis testing. Main effects were not observed for AKT3 rs2125230 and PRKCQ rs571715 in relation to PCA risk or disease progression in our study set. Wang and co-workers (2009) evaluated interactions among 5 apoptosis-related SNPs, including death receptor 4 (DR4), and pack-years of smoking in relation to bladder cancer using entropy-based MDR [51]. MDR analysis revealed a significant additive interaction between DR4 -397 G > T and smoking on bladder cancer. Unfortunately, this study analyzed a relatively small number of sequence variants in the DR4 apoptosis-related gene; hence, making it difficult to compare study findings. In a post-hoc analysis, we did not observe any interactions among the selected apoptosis-related markers and sources of reactive oxygen species, antioxidants, and anti-inflammatory agents, including cigarette smoking, dietary supplements, aspirin, ibuprofen, and meat-derived carcinogens (data not shown).

We considered the strengths, limitations and future directions of the current study. Although SEN-guided MDR only identified a nominally significant network for disease aggressiveness in the current study, this approach overcomes the computational challenge of detecting all possible two-, three- and four-way SNP combinations involved in PCA progression. SEN, in the current study, was used to prioritize > 10 million possible interactions by focusing on a "sub-network" of informative SNPs in relation to aggressive PCA. Reduction of our genetic data set to the most informative markers improved the feasibility to detect interactions that may have otherwise remained undetected. Unfortunately, the genome-wide association studies (GWAS) database used for the current study did not include apoptosis-related sequence variants (e.g., TNF-308 rs1800629, TNFSF10 rs1131532, BCL2 -936 rs2279115) previously reported in published cancer epidemiology studies [52-54]. Future studies in our laboratory will focus on high-throughput targeted sequencing to evaluate the impact of novel and commonly reported sequence variants on PCA susceptibility and disease prognosis. In light of recent GWAS reports, it is tempting to assume that extremely large case-control study sets, involving tens of thousands subjects are required to evaluate millions of SNP interactions in relation to PCA outcomes. However, the current study had adequate statistical power to evaluate individual and SNP combination effects in relation to prostate cancer. In particular, MDR has 80% statistical power to evaluate all possible two-, three-, and four-way gene-gene interactions with as low as 200 cases and 200 controls [55]. MDR remains effective even in the presence of 5% genotyping errors and/or 5% missing data [55]. It is anticipated that 5% of the about 13,000 possible interactions among 24 apoptosis-related SNPs will result in approximately 650 significant relationships due to chance alone. However, MDR coupled with permutation testing adjust for multiple comparison bias. Given the low prediction accuracy affiliated with the interaction between AKT3 and PRKCQ, our study findings require replication within independent study sets. However, recent simulation studies demonstrate that even modest disparities in genotype allele frequencies among study participants of independent study sets may interfere with the capacity to replicate complex interactions [56]. Thus, to ensure reproducibility within future studies, we plan to select study sets with similar genetic architecture (i.e., ancestry identification markers and SNP genotypes) as CGEMS PCA case-control subjects.

Conclusion

In summary, we have identified a marginal interaction between two apoptosis-related SNPs linked to PCA disease aggressiveness using SEN-guided MDR. Future molecular and pre-clinical studies may help to clarify the functionality of these genetic variants and their role in PCA disease progression. Our state-of-the art bioinformatics approach will enable researchers to pre-process millions of SNP interactions using GWAS or cancer consortia genotype data. In order to replicate findings within independent study sets, researchers will have to ensure that comparative sub-populations have the same genetic architecture and ancestry backgrounds. Therefore, the application of SEN-guided MDR may ultimately help researchers identify and validate gene-gene and/or gene-environment interactions to serve as effective cancer prognostication tools.

Methods

Study population

The Prostate, Lung, Colon, and Ovarian (PLCO) Cancer Screening Trial is a randomized, well-designed, and multi-center investigation sponsored and coordinated by the National Cancer Institute (NCI) [57,58]. Between 1993 and 2001, the PLCO Trial recruited study participants aged 55 to 81 years old to evaluate the effect of screening on disease specific mortality, relative to standard care. Men randomized to the PCA arm of the trial, received annual PSA and DRE exams. For the current study, 2,286 European-American men were included in our study if they had GWAS data available through the Cancer Genetic Markers of Susceptibility (CGEMS) data portal. (http://cgems.cancer.gov/) [38]. The CGEMS study population consists of nationally available genetic data for nearly 500,000 sequence variants. Incident PCA cases (488 non-aggressive and 687 aggressive) were identified through various sources including: screening exams; reports from patients, physicians, or relatives; or linkage with the National Death Index or state cancer registries. Incident PCA cases were pathologically confirmed with aggressive (Gleason score ≥ 7 and tumor stage III/IV) or non-aggressive [Gleason score (< 7) and tumor stage I/II] disease. Controls (n = 1,111) were matched to cases identified based on age, time since initial screening, and year of blood draw using incidence density sampling. All participants signed informed consent documents approved by both the NCI and local institutional review boards. Access to clinical and background data collected through examinations and questionnaires was approved for use by the PLCO.

Gene selection

A panel of 73 apoptosis-associated genes was selected from published cancer epidemiology or molecular biology studies as well as pathway databases and tools, including Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.genome.jp/kegg), Gene Ontology (GO), BioCarta (http://www.biocarta.com), ProteinLounge (http://www.proteinlounge.com), and Ingenuity (http://www.ingenuity.com) [59-63]. We searched PubMed for articles using the following keywords: [(prostate OR prostatic) AND (cancer OR neoplasms) AND apoptosis (gene variants OR single nucleotide polymorphisms OR targets) AND epidemiology]. Pathway tools, such as KEGG, BioCarta, ProteinLounge, and Ingenuity were used to visualize gene-gene and protein-protein interactions essential to regulating apoptosis [59-61,63]. The Entrez Gene database of the National Center for Biotechnology Information (NCBI, http://www.ncbi.nlm.nih.gov) was used to ensure that selected targets were involved in the apoptotic pathway [64].

Criteria for SNP selection and data management

Prior to uploading our initial list (i.e., 73 apoptosis genes) into the CGEMS data portal, we secured the HUGO gene name equivalents for the targets of interest using NCBI Entrez Gene. SNP profile data was available for more than 1197 SNPs [38]. However, upon further investigation, many of these SNPs were eliminated if they were: not located within 2.5 kb of the 5' or 3' end of the gene; or detected within a gene that was not related to apoptosis. Among SNPs that were associated with apoptosis, we focused on gene regions located in coding, promoter, or "near gene" regions. "Near gene" regions were defined as 2.5 kb up- and downstream of the 5' or 3' ends of the selected genes, respectively.

We excluded any sequence variants that had a minor allele frequency (MAF) < 5% as reported in the NCBI Entrez SNP [64]. In addition, one SNP was removed since its genotype frequency distribution among controls deviated substantially from the Hardy-Weinberg Equilibrium, p ≤ 0.005 (n = 1). Following data cleanup, 172 apoptotic SNPs were analyzed among 2,286 men of European descent (687 aggressive cases, 488 non-aggressive cases, and 1,111 controls).

Predicted function of selected SNPs

The SNPinfo (http://snpinfo.niehs.nih.gov/) webserver enabled us to annotate and/or predict the functional consequence of selected apoptosis variants, as summarized in Additional file 2[65]. This server consists of several pipelines to predict whether alternative alleles of a SNP may alter one or more of the following: transcriptional regulation via transcription factor binding site (TFBS) activity; premature termination of amino-acid sequence (stop codons); the splicing pattern or efficiency at mRNA splice sites, exonic splicing enhancers (ESE) or silencers (ESS); protein structures or properties by changing single amino acids (i.e., non-synonymous SNPs); or mRNA transcription or protein translation by altering microRNA (miRNA) binding sites.

Statistical analysis for single gene effects

Univariate and multivariate analyses were used to evaluate apoptosis associated SNPs among men of European descent in relation to PCA outcomes. To assess whether possession of one or more apoptotic alleles influence the risk of developing PCA, we tested for significant differences in the distribution of homozygous major, heterozygous, or homozygous minor genotypes between cases and controls using the chi-square test of heterogeneity. A case-only analysis was used to examine the relationship between apoptosis-related alleles and aggressive PCA. We evaluated differences in the distribution and inheritance of apoptosis-related genes comparing men with aggressive disease (Gleason score ≥ 7 and tumor stage > 2) to those with non-aggressive disease (i.e., other PCA cases that did not have both Gleason score ≥ 7 and tumor stage > 2). The associations between PCA outcomes and selected polymorphic genes, expressed as odds ratios (ORs) and corresponding 95% confidence intervals (CIs), were estimated using unconditional multivariate logistic regression (LR) models adjusted for age and family history of PCA. LR analyses for genetic variants and PCA outcomes were conducted using the major or common genotype as the referent category. All chi-square test and LR analyses were conducted using SAS 9.2 (SAS Institute Inc., Cary, NC) and SVS software (Golden Helix, Inc., Bozeman, MT). Statistical significance was assessed using a p-value < 0.05. Adjustments for multiple comparisons were made using Bonferroni correction and permutation testing.

Statistical power

We performed calculations to determine the statistical power of our sample to detect significant relationships between apoptosis-related polymorphisms and PCA development. The expected risk estimates of our study were determined by specifying values for a number of parameters, including: an average minor allele frequency of at least 26.0%, NCI's estimate of PCA disease prevalence (19%); statistical power (80%); sample size (1,175 PCA cases and 1,111 disease-free individuals or 687 aggressive PCA cases and 488 non-aggressive cases). We assumed the outcome was in complete linkage disequilibrium with an apoptosis-predisposing variant (r²= 1.0). Based on the aforementioned parameters, we have > 80% power to detect genetic markers with odds ratios (ORs) of ≥ 1.4 (or 0.71 for protective effects) for PCA risk and ≥ 1.6 (or 0.62 for protective effects) for disease aggressiveness, assuming a codominant model with 1 degree of freedom (df). Power calculations were performed using Power for Genetic Association Version 2 Software [66].

Analysis of gene interactions using multi-factor dimensionality reduction (MDR)

To evaluate the single- and joint- modifying effects of 172 candidate apoptosis-related SNPs within a large dataset, such as CGEMS, is computationally challenging [67,68]. In order to overcome this problem, open-source and freely available MDR 2.0 software (htpp://http://www.epistasis.org) was used to analyze interactions among apoptotic sequence variants in relation to PCA outcomes [69]. MDR is a data reduction approach, designed to detect and characterize high-order interactions in case-control studies. With MDR, information from various genetic loci is categorized into high-risk or low-risk groups of disease. The resulting one-dimensional multi-locus genotype variable is then evaluated for its ability to classify and predict a disease outcome through cross-validation and permutation testing. MDR uses a 10-fold cross validation to estimate the testing accuracy of a model by leaving out 1/10^thof the data as an independent test set. The model is developed on 9/10^thof the data and then evaluated on the remaining test set. This process is repeated for each 1/10^thof the data and the resulting prediction accuracies are averaged.

In the current study, the model with the greatest cross validation consistency (i.e., CVC ≥ 8/10) and highest prediction accuracy [i.e., Average Testing Accuracy (ATA)] was selected as the best predictor of disease outcome. Accuracy is a function of the percentage of true positives (TP), true negatives (TN), false positives (FP), and false negatives (FN) as defined as (TP + TN)/(TP + TN + FP + FN). ATAs are averaged across all 10 pieces of the data, in order to provide an estimate of the predictive ability of the loci in relation to the outcome of interest. We used cross-validation consistency (CVC) to measure the degree to which the same best MDR model was selected across the 10 divisions of the data. Models with a CVC of ≥ 8/10 using a 10-fold cross-validation were considered more carefully. If the MDR model met the CVC criteria, we selected models that had the highest ATAs. The combination of CVC and permutation testing were used to control for multiple hypothesis testing. Permutation testing results were statistically significant at the ≤ 0.05 level.

The current version of MDR used in this study enabled the incorporation and adjustment of multiple covariates [70]. To remove the covariate effects [i.e., age-group and family history of PCA (yes or no)], we integrated two sampling methods (i.e., over- and under-sampling). This approach is computationally efficient, since it allows for the adjustment of multiple covariates without significantly increasing computational burden.

Visualization of interaction models using interaction entropy algorithms, hierarchical graphs and statistical epistasis network (SEN)

The interaction entropy algorithm, based on information theory, is a method to verify, visualize, and interpret combination effects identified by parametric (e.g., LR) and non-parametric (e.g., MDR) statistical test[43,71-74]. Jakulin and Brakto (2003) have developed a metric to gauge whether the gain in information (i.e., information gain) about a class variable (i.e., ability to predict disease status) from a combination of two variables provides more information than each variable considered independently [72,73]. In the current study, measures of interaction were used to build interaction entropy graphs to visualize and interpret interactions among selected apoptosis-related markers in relation to PCA risk and disease progression. Individual and all possible pairwise loci were assigned a mutual information (I) percentage score in relation to disease risk or aggressiveness; whereby typical scores for genetic loci are < 5%. Pairwise SNP combinations were deemed important if the joint mutual information [I = (SNP₁, SNP₂; disease status)] was greater than the total mutual information of each individual locus considered separately [I (SNP₁; disease status) + I (SNP₂; disease status)]. Interactions were further visualized using an interaction entropy graph, which uses a color-coding scheme to interpret interactions. Strong interacting factors, coded either red or orange, indicated high and medium levels of synergistic effects on outcomes, respectively. Weak interacting factors, coded either blue or green, denoted high or modest levels of redundancy between markers, respectively. Gold depicted independence and a midway point between synergy and redundancy. Design of interaction entropy graphs was accomplished using the Orange software [75].

An important limitation of any gene-gene interaction analysis method is the combinatorial nature of the problem. Exhaustive analysis of all possible two-, three-, four-way combinations of SNPs creates a computational burden and significant multiple testing problems. We addressed this issue by pre-processing the data using statistical epistasis network (SEN) modeling, which has been described elsewhere [39]. Briefly, this tool uses information theory, as previously described, to develop a network based on the amount of mutual information from SNP₁, SNP₂, and a class variable (i.e., disease status). Thus, it describes the pairwise interactive effect on a discrete outcome. An epistasis network, or graph, consists of a set of vertices (individual SNPs) and edges (pairwise interactions) that attach the vertices. In the current study, the statistical epistasis networks (SEN) were built by incrementally adding edges if their strengths were greater than a given threshold. Each vertex and each edge connecting two vertices were assigned a weight commensurate with the strength of main and interactive effects, respectively. The weight of a vertex is depicted by its size; whereby, a larger size is indicative of a higher main effect. Similarly, the thickness of the line between SNP pairs is directly proportional to the strength of their interaction.

We used SEN to visually summarize the strongest individual SNP effects and SNP-pairwise interactions. The topology of the graph, consisting of SNPs and their pairwise interactions, was used to guide MDR modeling. This resulted in a reduction in the total number of SNP combinations. For instance, we reduced the gene-gene interaction analysis from 172 to 24 SNPs that resided within the major connected component with the largest number of vertices (i.e., main effects). The derived major component plot was validated by comparing: (1) the number of vertices and edges among the various connected "sub-networks"; and (2) the size of the major component plot of real data versus randomized data that underwent 1000 permutations. The latter compared data using a pairwise interaction strength cut-off of 1% and significance level of 0.05. To ensure that no significant interactions were identified, non-informative SNPs outside the major component plot were also analyzed by MDR. In essence, SEN-guided MDR facilitated logical prioritization of SNPs for gene-gene interaction analysis in relation to PCA outcomes.

Abbreviations

PCA: (PCA); SNP: (single nucleotide polymorphism); MDR: (multifactor dimensionality reduction); SEN: (Statistical Epistasis Network).

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

NAL: conceptualized the project, conducted the statistical analysis and initial MDR modeling, interpreted study results, and composed the initial manuscript draft. ENR, SY, LRK: Intellectually contributed toward the introduction and/or discussion. TH: performed statistical epistasis network modeling. JZ and JR: conducted MDR modeling and/or SEN-guided MDR analysis. GNB & LRK: co-supervised data-management and statistical/bioinformatics analyses. LRK, DWH, GNB, JHM, KSK: served as mentors in the study design and implementation as well as data interpretation. All authors contributed to draft reviewing, editing, and approving the final manuscript draft.

Pre-publication history

The pre-publication history for this paper can be accessed here:

http://www.biomedcentral.com/1755-8794/5/11/prepub

Supplementary Material

Additional file 1

Table S1. Corresponding Gene and dbSNP IDs for apoptosis-related sequence variants depicted in the Statistical Epistasis Network Modeling Graph, presented in Figure Figure11.

Click here for file^{(18K, DOCX)}

Additional file 2

Table S2. Functional consequence of selected apoptosis-related Polymorphisms.

Click here for file^{(270K, DOC)}

Acknowledgements

The authors thank CGEMS and PLCO projects for making their data available for use in this study. This work was partially supported by the James Graham Brown Cancer Center (JGBCC) Pilot Project Initiative Grant to LRK, and the JGBCC Bucks for Brains "Our Highest Potential" in Cancer Research Endowment to LRK, National Cancer Institute/National Institute of Health grants (R03 CA128028, 3R01 CA034627-19S) to LRK, and the National Institute of Health/National Institute of Environmental Health Sciences grants P30-ES014443 to GNB, National Institute of Environmental Health Sciences training grant T32 ES011564 to DWH, and P20-MD000175 NIH NCMHD to KSK.

References

American Cancer S. Cancer Facts and Figures 2012. Atlanta, Georgia: American Cancer Society; 2012.
Tapia-Vieyra JV, Mas-Oliva J. Apoptosis and cell death channels in prostate cancer. ArchMedRes. 2001;32(3):175–185. [PubMed]
Zhivotovsky B, Orrenius S. Carcinogenesis and apoptosis: paradigms and paradoxes. Carcinogenesis. 2006;27(10):1939–1945. [PubMed]
Limoli CL, Hartmann A, Shephard L, Yang CR, Boothman DA, Bartholomew J, Morgan WF. Apoptosis, reproductive failure, and oxidative stress in Chinese hamster ovary cells with compromised genomic integrity. Cancer Res. 1998;58(16):3712–3718. [PubMed]
Danial NN, Korsmeyer SJ. Cell death: critical control points. Cell. 2004;116(2):205–219. [PubMed]
Furuya Y, Krajewski S, Epstein JI, Reed JC, Isaacs JT. Expression of bcl-2 and the progression of human and rodent prostatic cancers. Clinical cancer research: an official journal of the American Association for Cancer Research. 1996;2(2):389–398. [PubMed]
Colombel M, Symmans F, Gil S, O'Toole KM, Chopin D, Benson M, Olsson CA, Korsmeyer S, Buttyan R. Detection of the apoptosis-suppressing oncoprotein bc1-2 in hormone-refractory human prostate cancers. Am J Pathol. 1993;143(2):390–400. [PubMed]
McDonnell TJ, Troncoso P, Brisbay SM, Logothetis C, Chung LW, Hsieh JT, Tu SM, Campbell ML. Expression of the protooncogene bcl-2 in the prostate and its association with emergence of androgen-independent prostate cancer. Cancer Res. 1992;52(24):6940–6944. [PubMed]
Fleischmann A, Huland H, Mirlacher M, Wilczak W, Simon R, Erbersdobler A, Sauter G, Schlomm T. Prognostic relevance of Bcl-2 overexpression in surgically treated prostate cancer is not caused by increased copy number or translocation of the gene. The Prostate. 2011. in press doi: 10.1002/pros.21504. [PubMed]
Abate-Shen C, Shen MM. Molecular genetics of prostate cancer. Genes Dev. 2000;14(19):2410–2434. [PubMed]
Thomas DJ, Robinson M, King P, Hasan T, Charlton R, Martin J, Carr TW, Neal DE. p53 expression and clinical outcome in prostate cancer. Br J Urol. 1993;72(5 Pt 2):778–781. [PubMed]
Shurbaji MS, Kalbfleisch JH, Thurmond TS. Immunohistochemical detection of p53 protein as a prognostic indicator in prostate cancer. Hum Pathol. 1995;26(1):106–109. [PubMed]
Stackhouse GB, Sesterhenn IA, Bauer JJ, Mostofi FK, Connelly RR, Srivastava SK, Moul JW. p53 and bcl-2 immunohistochemistry in pretreatment prostate needle biopsies to predict recurrence of prostate cancer after radical prostatectomy. J Urol. 1999;162(6):2040–2045. [PubMed]
Pickard MR, Edwards SE, Cooper CS, Williams GT. Apoptosis regulators Fau and Bcl-G are down-regulated in prostate cancer. Prostate. 2010;70(14):1513–1523. [PubMed]
Boehm JS, Zhao JJ, Yao J, Kim SY, Firestein R, Dunn IF, Sjostrom SK, Garraway LA, Weremowicz S, Richardson AL. et al. Integrative genomic approaches identify IKBKE as a breast cancer oncogene. Cell. 2007;129(6):1065–1079. [PubMed]
Guo JP, Shu SK, He L, Lee YC, Kruk PA, Grenman S, Nicosia SV, Mor G, Schell MJ, Coppola D. et al. Deregulation of IKBKE is associated with tumor progression, poor prognosis, and cisplatin resistance in ovarian cancer. Am J Pathol. 2009;175(1):324–333. [PubMed]
Hill KM, Kalifa S, Das JR, Bhatti T, Gay M, Williams D, Taliferro-Smith L, De Marzo AM. The role of PI 3-kinase p110beta in AKT signally, cell survival, and proliferation in human prostate cancer cells. Prostate. 2010;70(7):755–764. [PubMed]
Gasparian AV, Yao YJ, Kowalczyk D, Lyakh LA, Karseladze A, Slaga TJ, Budunova IV. The role of IKK in constitutive activation of NF-kappaB transcription factor in prostate carcinoma cells. J Cell Sci. 2002;115(Pt 1):141–151. [PubMed]
Jia S, Liu Z, Zhang S, Liu P, Zhang L, Lee SH, Zhang J, Signoretti S, Loda M, Roberts TM. et al. Essential roles of PI(3)K-p110beta in cell growth, metabolism and tumorigenesis. Nature. 2008;454(7205):776–779. [PMC free article] [PubMed]
Ricks-Santi L, Mason T, Apprey V, Ahaghotu C, McLauchlin A, Josey D, Bonney G, Dunston GM. p53 Pro72Arg polymorphism and prostate cancer in men of African descent. Prostate. 2010;70(16):1739–1745. [PMC free article] [PubMed]
Srivastava K, Srivastava A, Sharma KL, Mittal B. Candidate gene studies in gallbladder cancer: a systematic review and meta-analysis. Mutat Res. 2011;728(1-2):67–79. [PMC free article] [PubMed]
Lee EB, Jeon HS, Yoo SS, Choi YY, Kang HG, Cho S, Cha SI, Choi JE, Park TI, Lee BH. et al. Polymorphisms in apoptosis-related genes and survival of patients with early-stage non-small-cell lung cancer. Ann Surg Oncol. 2010;17(10):2608–2618. [PubMed]
Bye H, Prescott NJ, Matejcic M, Rose E, Lewis CM, Parker MI, Mathew CG. Population-specific genetic associations with oesophageal squamous cell carcinoma in South Africa. Carcinogenesis. 2011;32(12):1855–1861. [PMC free article] [PubMed]
Campa D, Kaaks R, Le Marchand L, Haiman CA, Travis RC, Berg CD, Buring JE, Chanock SJ, Diver WR, Dostal L. et al. Interactions between genetic variants and breast cancer risk factors in the breast and prostate cancer cohort consortium. J National Cancer Inst. 2011;103(16):1252–1263. [PMC free article] [PubMed]
Broeks A, Schmidt MK, Sherman ME, Couch FJ, Hopper JL, Dite GS, Apicella C, Smith LD, Hammet F, Southey MC. et al. Low penetrance breast cancer susceptibility loci are associated with specific breast tumor subtypes: findings from the Breast Cancer Association Consortium. Hum Mol Genet. 2011;20(16):3289–3303. [PMC free article] [PubMed]
Couch FJ, Wang X, McWilliams RR, Bamlet WR, de Andrade M, Petersen GM. Association of breast cancer susceptibility variants with risk of pancreatic cancer. Cancer epidemiology, biomarkers & prevention: a publication of the American Association for Cancer Research, cosponsored by the American Society of Preventive Oncology. 2009;18(11):3044–3048. [PMC free article] [PubMed]
Sergentanis TN, Economopoulos KP. Association of two CASP8 polymorphisms with breast cancer risk: a meta-analysis. Breast Cancer Res Treat. 2010;120(1):229–234. [PubMed]
Enjuanes A, Benavente Y, Bosch F, Martin-Guerrero I, Colomer D, Perez-Alvarez S, Reina O, Ardanaz MT, Jares P, Garcia-Orad A. et al. Genetic variants in apoptosis and immunoregulation-related genes are associated with risk of chronic lymphocytic leukemia. Cancer Res. 2008;68(24):10178–10186. [PubMed]
Sigurdson AJ, Bhatti P, Doody MM, Hauptmann M, Bowen L, Simon SL, Weinstock RM, Linet MS, Rosenstein M, Stovall M. et al. Polymorphisms in apoptosis- and proliferation-related genes, ionizing radiation exposure, and risk of breast cancer among U.S. Radiologic Technologists. Cancer epidemiology, biomarkers & prevention: a publication of the American Association for Cancer Research, cosponsored by the American Society of Preventive Oncology. 2007;16(10):2000–2007. [PubMed]
Cox A, Dunning AM, Garcia-Closas M, Balasubramanian S, Reed MW, Pooley KA, Scollen S, Baynes C, Ponder BA, Chanock S. et al. A common coding variant in CASP8 is associated with breast cancer risk. Nat Genet. 2007;39(3):352–358. [PubMed]
Lan Q, Zheng T, Chanock S, Zhang Y, Shen M, Wang SS, Berndt SI, Zahm SH, Holford TR, Leaderer B. et al. Genetic variants in caspase genes and susceptibility to non-Hodgkin lymphoma. Carcinogenesis. 2007;28(4):823–827. [PubMed]
Lou Y, Fang CQ, Li JH. A study on the expression of CASP9 gene and its polymorphism distribution in non-small cell lung cancer. Zhonghua Yi Xue Yi Chuan Xue Za Zhi. 2007;24(1):59–62. [PubMed]
Consortium BCA. Commonly studied single-nucleotide polymorphisms and breast cancer: results from the Breast Cancer Association Consortium. J National Cancer Inst. 2006;98(19):1382–1396. [PubMed]
Du J, Huo J, Shi J, Yuan Z, Zhang C, Fu W, Jiang H, Yi Q, Hou J. Polymorphisms of nuclear factor-kappaB family genes are associated with development of multiple myeloma and treatment outcome in patients receiving bortezomib-based regimens. Haematologica. 2011;96(5):729–737. [PubMed]
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. Scand J Gastroenterol. 2007;42(3):345–350. [PubMed]
Schildkraut JM, Iversen ES, Wilson MA, Clyde MA, Moorman PG, Palmieri RT, Whitaker R, Bentley RC, Marks JR, Berchuck A. Association between DNA damage response and repair genes and risk of invasive serous ovarian cancer. PLoS One. 2010;5(4):e10061. [PMC free article] [PubMed]
Sprague BL, Trentham-Dietz A, Garcia-Closas M, Newcomb PA, Titus-Ernstoff L, Hampton JM, Chanock SJ, Haines JL, Egan KM. Genetic variation in TP53 and risk of breast cancer in a population-based case control study. Carcinogenesis. 2007;28(8):1680–1686. [PubMed]
Yeager M, Orr N, Hayes RB, Jacobs KB, Kraft P, Wacholder S, Minichiello MJ, Fearnhead P, Yu K, Chatterjee N. et al. Genome-wide association study of prostate cancer identifies a second risk locus at 8q24. Nat Genet. 2007;39(5):645–649. [PubMed]
Hu T, Sinnott-Armstrong NA, Kiralis JW, Andrew AS, Karagas MR, Moore JH. Characterizing genetic interactions in human disease association studies using statistical epistasis networks. BMC Bioinforma. 2011;12:364. [PMC free article] [PubMed]
Duell EJ, Bracci PM, Moore JH, Burk RD, Kelsey KT, Holly EA. Detecting pathway-based gene-gene and gene-environment interactions in pancreatic cancer. Cancer epidemiology, biomarkers & prevention: a publication of the American Association for Cancer Research, cosponsored by the American Society of Preventive Oncology. 2008;17(6):1470–1479. [PubMed]
Mason RA, Morlock EV, Karagas MR, Kelsey KT, Marsit CJ, Schned AR, Andrew AS. EGFR pathway polymorphisms and bladder cancer susceptibility and prognosis. Carcinogenesis. 2009;30(7):1155–1160. [PMC free article] [PubMed]
Andrew AS, Karagas MR, Nelson HH, Guarrera S, Polidoro S, Gamberini S, Sacerdote C, Moore JH, Kelsey KT, Demidenko E. et al. DNA repair polymorphisms modify bladder cancer risk: a multi-factor analytic strategy. HumHered. 2008;65(2):105–118. [PMC free article] [PubMed]
Andrew AS, Nelson HH, Kelsey KT, Moore JH, Meng AC, Casella DP, Tosteson TD, Schned AR, Karagas MR. Concordance of multiple analytical approaches demonstrates a complex relationship between DNA repair gene SNPs, smoking and bladder cancer susceptibility. Carcinogenesis. 2006;27(5):1030–1037. [PubMed]
Lavender NA, Benford ML, VanCleave TT, Brock GN, Kittles RA, Moore JH, Hein DW, Kidd LC. Examination of polymorphic glutathione S-transferase (GST) genes, tobacco smoking and prostate cancer risk among men of African descent: a case-control study. BMC Cancer. 2009;9:397. [PMC free article] [PubMed]
Lavender NA, Komolafe OO, Benford M, Brock G, Moore JH, Vancleave TT, States JC, Kittles RA, Kidd LC. No association between variant DNA repair genes and prostate cancer risk among men of African descent. Prostate. 2010;70(2):113–119. [PMC free article] [PubMed]
Chen M, Kamat AM, Huang M, Grossman HB, Dinney CP, Lerner SP, Wu X, Gu J. High-order interactions among genetic polymorphisms in nucleotide excision repair pathway genes and smoking in modulating bladder cancer risk. Carcinogenesis. 2007;28(10):2160–2165. [PubMed]
Benford ML, VanCleave TT, Lavender NA, Kittles RA, Kidd LR. 8q24 sequence variants in relation to prostate cancer risk among men of African descent: a case-control study. BMC Cancer. 2010;10:334. [PMC free article] [PubMed]
VanCleave TT, Moore JH, Benford ML, Brock GN, Kalbfleisch T, Baumgartner RN, Lillard JW Jr, Kittles RA, Kidd LC. Interaction among variant vascular endothelial growth factor (VEGF) and its receptor in relation to prostate cancer risk. Prostate. 2010;70(4):341–352. [PubMed]
Nakatani K, Thompson DA, Barthel A, Sakaue H, Liu W, Weigel RJ, Roth RA. Up-regulation of Akt3 in estrogen receptor-deficient breast cancers and androgen-independent prostate cancer lines. J Biol Chem. 1999;274(31):21528–21532. [PubMed]
Bauer B, Krumbock N, Fresser F, Hochholdinger F, Spitaler M, Simm A, Uberall F, Schraven B, Baier G. Complex formation and cooperation of protein kinase C theta and Akt1/protein kinase B alpha in the NF-kappa B transactivation cascade in Jurkat T cells. J Biol Chem. 2001;276(34):31627–31634. [PubMed]
Wang M, Cheng G, Zhang Z, Fu G. Genetic variants in the death receptor 4 gene contribute to susceptibility to bladder cancer. Mutat Res. 2009;661(1-2):85–92. [PubMed]
Nam RK, Zhang WW, Jewett MA, Trachtenberg J, Klotz LH, Emami M, Sugar L, Sweet J, Toi A, Narod SA. The use of genetic markers to determine risk for prostate cancer at prostate biopsy. Clinical cancer research: an official journal of the American Association for Cancer Research. 2005;11(23):8391–8397. [PubMed]
Jung JH, Chae YS, Moon JH, Kang BW, Kim JG, Sohn SK, Park JY, Lee MH, Park HY. TNF superfamily gene polymorphism as prognostic factor in early breast cancer. J Cancer Res Clin Oncol. 2010;136(5):685–694. [PubMed]
Kidd LR, Coulibaly A, Templeton TM, Chen W, Long LO, Mason T, Bonilla C, Akereyeni F, Freeman V, Isaacs W. et al. Germline BCL-2 sequence variants and inherited predisposition to prostate cancer. Prostate Cancer ProstaticDis. 2006;9(3):284–292. [PubMed]
Ritchie MD, Hahn LW, Moore JH. Power of multifactor dimensionality reduction for detecting gene-gene interactions in the presence of genotyping error, missing data, phenocopy, and genetic heterogeneity. Genetic epidemiology. 2003;24(2):150–157. [PubMed]
Greene CS, Penrod NM, Williams SM, Moore JH. Failure to replicate a genetic association may provide important clues about genetic architecture. PLoS One. 2009;4(6):e5639. [PMC free article] [PubMed]
Gohagan JK, Prorok PC, Hayes RB, Kramer BS. The Prostate, Lung, Colorectal and Ovarian (PLCO) Cancer Screening Trial of the National Cancer Institute: history, organization, and status. Control ClinTrials. 2000;21(6 Suppl):251S–272S. [PubMed]
Hayes RB, Sigurdson A, Moore L, Peters U, Huang WY, Pinsky P, Reding D, Gelmann EP, Rothman N, Pfeiffer RM. et al. Methods for etiologic and early marker investigations in the PLCO trial. MutatRes. 2005;592(1-2):147–154. [PubMed]
Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008. pp. D480–D484. [PMC free article] [PubMed]
Kanehisa M, Goto S, Hattori M, oki-Kinoshita KF, Itoh M, Kawashima S, Katayama T, Araki M, Hirakawa M. From genomics to chemical genomics: new developments in KEGG. Nucleic Acids Res. 2006. pp. D354–D357. [PMC free article] [PubMed]
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30. [PMC free article] [PubMed]
Ingenuity Systems. 2008.
Biocarta LLC. 2000.
National Center for Biotechnology Information. 2011.
Xu J, Kibel AS, Hu JJ, Turner AR, Pruett K, Zheng SL, Sun J, Isaacs SD, Wiley KE, Kim ST. et al. Prostate cancer risk associated loci in African Americans. Cancer EpidemiolBiomarkers Prev. 2009;18(7):2145–2149. [PMC free article] [PubMed]
Menashe I, Rosenberg PS, Chen BE. PGA: power calculator for case-control genetic association analyses. BMC Genet. 2008;9:36. [PMC free article] [PubMed]
Greene CS, Penrod NM, Kiralis J, Moore JH. Spatially uniform relieff (SURF) for computationally-efficient filtering of gene-gene interactions. BioData mining. 2009;2(1):5. [PMC free article] [PubMed]
Moore JH, Asselbergs FW, Williams SM. Bioinformatics challenges for genome-wide association studies. Bioinformatics. 2010;26(4):445–455. [PMC free article] [PubMed]
Moore JH, Gilbert JC, Tsai CT, Chiang FT, Holden T, Barney N, White BC. A flexible computational framework for detecting, characterizing, and interpreting statistical patterns of epistasis in genetic studies of human disease susceptibility. J Theor Biol. 2006;241(2):252–261. [PubMed]
Gui J, Andrew AS, Andrews P, Nelson HM, Kelsey KT, Karagas MR, Moore JH. A simple and computationally efficient sampling approach to covariate adjustment for multifactor dimensionality reduction analysis of epistasis. Hum Hered. 2010;70(3):219–225. [PMC free article] [PubMed]
McGill WL. Multivariate information transmission. Psychometrika vol. 1954;19:97–116.
Jakulin A, Bratko I. Analyzing attribute interations. Lecture Notes in Artificial Intelligence. 2003;2838:229.
Jakulin A, Bratko I, Smrke D, Demsar J, Supan B. Attribute interactions in medical data analysis. Lecture Notes in Artificial Intelligence. 2003;2780:229.
Moore JH, Gilbert JC, Tsai CT, Chiang FT, Holden T, Barney N, White BC. A flexible computational framework for detecting, characterizing, and interpreting statistical patterns of epistasis in genetic studies of human disease susceptibility. JTheorBiol. 2006;241(2):252–261. [PubMed]
Mramor M, Leban G, Demsar J, Zupan B. Visualization-based cancer microarray data classification analysis. Bioinformatics. 2007;23(16):2147–2154. [PubMed]

Articles from BMC Medical Genomics are provided here courtesy of
BioMed Central