Log-linear model-based multifactor dimensionality reduction method to detect gene gene interactions

MOTIVATION: The identification and characterization of susceptibility genes that influence the risk of common and complex diseases remains a statistical and computational challenge in genetic association studies. This is partly because the effect of any single genetic variant for a common and complex disease may be dependent on other genetic variants (gene-gene interaction) and environmental factors (gene-environment interaction). To address this problem, the multifactor dimensionality reduction (MDR) method has been proposed by Ritchie et al. to detect gene-gene interactions or gene-environment interactions. The MDR method identifies polymorphism combinations associated with the common and complex multifactorial diseases by collapsing high-dimensional genetic factors into a single dimension. That is, the MDR method classifies the combination of multilocus genotypes into high-risk and low-risk groups based on a comparison of the ratios of the numbers of cases and controls. When a high-order interaction model is considered with multi-dimensional factors, however, there may be many sparse or empty cells in the contingency tables. The MDR method cannot classify an empty cell as high risk or low risk and leaves it as undetermined. RESULTS: In this article, we propose the log-linear model-based multifactor dimensionality reduction (LM MDR) method to improve the MDR in classifying sparse or empty cells. The LM MDR method estimates frequencies for empty cells from a parsimonious log-linear model so that they can be assigned to high-and low-risk groups. In addition, LM MDR includes MDR as a special case when the saturated log-linear model is fitted. Simulation studies show that the LM MDR method has greater power and smaller error rates than the MDR method. The LM MDR method is also compared with the MDR method using as an example sporadic Alzheimer's disease.     

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

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.     

A flexible computational framework for detecting, characterizing, and interpreting statistical patterns of epistasis in genetic studies of human disease susceptibility

Detecting, characterizing, and interpreting gene-gene interactions or epistasis in studies of human disease susceptibility is both a mathematical and a computational challenge. To address this problem, we have previously developed a multifactor dimensionality reduction (MDR) method for collapsing high-dimensional genetic data into a single dimension (i.e. constructive induction) thus permitting interactions to be detected in relatively small sample sizes. In this paper, we describe a comprehensive and flexible framework for detecting and interpreting gene-gene interactions that utilizes advances in information theory for selecting interesting single-nucleotide polymorphisms (SNPs), MDR for constructive induction, machine learning methods for classification, and finally graphical models for interpretation. We illustrate the usefulness of this strategy using artificial datasets simulated from several different two-locus and three-locus epistasis models. We show that the accuracy, sensitivity, specificity, and precision of a na?ve Bayes classifier are significantly improved when SNPs are selected based on their information gain (i.e. class entropy removed) and reduced to a single attribute using MDR. We then apply this strategy to detecting, characterizing, and interpreting epistatic models in a genetic study (n = 500) of atrial fibrillation and show that both classification and model interpretation are significantly improved.     

Disease-driven detection of differential inherited SNP modules from SNP network

Detection of the synergetic effects between variants, such as single-nucleotide polymorphisms (SNPs), is crucial for understanding the genetic characters of complex diseases. Here, we proposed a two-step approach to detect differentially inherited SNP modules (synergetic SNP units) from a SNP network. First, SNP-SNP interactions are identified based on prior biological knowledge, such as their adjacency on the chromosome or degree of relatedness between the functional relationships of their genes. These interactions form SNP networks. Second, disease-risk SNP modules (or sub-networks) are prioritised by their differentially inherited properties in IBD (Identity by Descent) profiles of affected and unaffected sibpairs. The search process is driven by the disease information and follows the structure of a SNP network. Simulation studies have indicated that this approach achieves high accuracy and a low false-positive rate in the identification of known disease-susceptible SNPs. Applying this method to an alcoholism dataset, we found that flexible patterns of susceptible SNP combinations do play a role in complex diseases, and some known genes were detected through these risk SNP modules. One example is GRM7, a known alcoholism gene successfully detected by a SNP module comprised of two SNPs, but neither of the two SNPs was significantly associated with the disease in single-locus analysis. These identified genes are also enriched in some pathways associated with alcoholism, including the calcium signalling pathway, axon guidance and neuroactive ligand-receptor interaction. The integration of network biology and genetic analysis provides putative functional bridges between genetic variants and candidate genes or pathways, thereby providing new insight into the aetiology of complex diseases.     

Alternative contingency table measures improve the power and detection of multifactor dimensionality reduction

Background  

Multifactor Dimensionality Reduction (MDR) has been introduced previously as a non-parametric statistical method for detecting gene-gene interactions. MDR performs a dimensional reduction by assigning multi-locus genotypes to either high- or low-risk groups and measuring the percentage of cases and controls incorrectly labelled by this classification – the classification error. The combination of variables that produces the lowest classification error is selected as the best or most fit model. The correctly and incorrectly labelled cases and controls can be expressed as a two-way contingency table. We sought to improve the ability of MDR to detect gene-gene interactions by replacing classification error with a different measure to score model quality.     

High-order SNP combinations associated with complex diseases: efficient discovery, statistical power and functional interactions

There has been increased interest in discovering combinations of single-nucleotide polymorphisms (SNPs) that are strongly associated with a phenotype even if each SNP has little individual effect. Efficient approaches have been proposed for searching two-locus combinations from genome-wide datasets. However, for high-order combinations, existing methods either adopt a brute-force search which only handles a small number of SNPs (up to few hundreds), or use heuristic search that may miss informative combinations. In addition, existing approaches lack statistical power because of the use of statistics with high degrees-of-freedom and the huge number of hypotheses tested during combinatorial search. Due to these challenges, functional interactions in high-order combinations have not been systematically explored. We leverage discriminative-pattern-mining algorithms from the data-mining community to search for high-order combinations in case-control datasets. The substantially improved efficiency and scalability demonstrated on synthetic and real datasets with several thousands of SNPs allows the study of several important mathematical and statistical properties of SNP combinations with order as high as eleven. We further explore functional interactions in high-order combinations and reveal a general connection between the increase in discriminative power of a combination over its subsets and the functional coherence among the genes comprising the combination, supported by multiple datasets. Finally, we study several significant high-order combinations discovered from a lung-cancer dataset and a kidney-transplant-rejection dataset in detail to provide novel insights on the complex diseases. Interestingly, many of these associations involve combinations of common variations that occur in small fractions of population. Thus, our approach is an alternative methodology for exploring the genetics of rare diseases for which the current focus is on individually rare variations.     

GPNN: Power studies and applications of a neural network method for detecting gene-gene interactions in studies of human disease

Background  

The identification and characterization of genes that influence the risk of common, complex multifactorial disease primarily through interactions with other genes and environmental factors remains a statistical and computational challenge in genetic epidemiology. We have previously introduced a genetic programming optimized neural network (GPNN) as a method for optimizing the architecture of a neural network to improve the identification of gene combinations associated with disease risk. The goal of this study was to evaluate the power of GPNN for identifying high-order gene-gene interactions. We were also interested in applying GPNN to a real data analysis in Parkinson's disease.     

Lower-order effects adjustment in quantitative traits model-based multifactor dimensionality reduction

Identifying gene-gene interactions or gene-environment interactions in studies of human complex diseases remains a big challenge in genetic epidemiology. An additional challenge, often forgotten, is to account for important lower-order genetic effects. These may hamper the identification of genuine epistasis. If lower-order genetic effects contribute to the genetic variance of a trait, identified statistical interactions may simply be due to a signal boost of these effects. In this study, we restrict attention to quantitative traits and bi-allelic SNPs as genetic markers. Moreover, our interaction study focuses on 2-way SNP-SNP interactions. Via simulations, we assess the performance of different corrective measures for lower-order genetic effects in Model-Based Multifactor Dimensionality Reduction epistasis detection, using additive and co-dominant coding schemes. Performance is evaluated in terms of power and familywise error rate. Our simulations indicate that empirical power estimates are reduced with correction of lower-order effects, likewise familywise error rates. Easy-to-use automatic SNP selection procedures, SNP selection based on "top" findings, or SNP selection based on p-value criterion for interesting main effects result in reduced power but also almost zero false positive rates. Always accounting for main effects in the SNP-SNP pair under investigation during Model-Based Multifactor Dimensionality Reduction analysis adequately controls false positive epistasis findings. This is particularly true when adopting a co-dominant corrective coding scheme. In conclusion, automatic search procedures to identify lower-order effects to correct for during epistasis screening should be avoided. The same is true for procedures that adjust for lower-order effects prior to Model-Based Multifactor Dimensionality Reduction and involve using residuals as the new trait. We advocate using "on-the-fly" lower-order effects adjusting when screening for SNP-SNP interactions using Model-Based Multifactor Dimensionality Reduction analysis.     

On the Use of Variance per Genotype as a Tool to Identify Quantitative Trait Interaction Effects: A Report from the Women's Genome Health Study

Testing for genetic effects on mean values of a quantitative trait has been a very successful strategy. However, most studies to date have not explored genetic effects on the variance of quantitative traits as a relevant consequence of genetic variation. In this report, we demonstrate that, under plausible scenarios of genetic interaction, the variance of a quantitative trait is expected to differ among the three possible genotypes of a biallelic SNP. Leveraging this observation with Levene''s test of equality of variance, we propose a novel method to prioritize SNPs for subsequent gene–gene and gene–environment testing. This method has the advantageous characteristic that the interacting covariate need not be known or measured for a SNP to be prioritized. Using simulations, we show that this method has increased power over exhaustive search under certain conditions. We further investigate the utility of variance per genotype by examining data from the Women''s Genome Health Study. Using this dataset, we identify new interactions between the LEPR SNP rs12753193 and body mass index in the prediction of C-reactive protein levels, between the ICAM1 SNP rs1799969 and smoking in the prediction of soluble ICAM-1 levels, and between the PNPLA3 SNP rs738409 and body mass index in the prediction of soluble ICAM-1 levels. These results demonstrate the utility of our approach and provide novel genetic insight into the relationship among obesity, smoking, and inflammation.     

Gene-based Genomewide Association Analysis: A Comparison Study

The study of gene-based genetic associations has gained conceptual popularity recently. Biologic insight into the etiology of a complex disease can be gained by focusing on genes as testing units. Several gene-based methods (e.g., minimum p-value (or maximum test statistic) or entropy-based method) have been developed and have more power than a single nucleotide polymorphism (SNP)-based analysis. The objective of this study is to compare the performance of the entropy-based method with the minimum p-value and single SNP–based analysis and to explore their strengths and weaknesses. Simulation studies show that: 1) all three methods can reasonably control the false-positive rate; 2) the minimum p-value method outperforms the entropy-based and the single SNP–based method when only one disease-related SNP occurs within the gene; 3) the entropy-based method outperforms the other methods when there are more than two disease-related SNPs in the gene; and 4) the entropy-based method is computationally more efficient than the minimum p-value method. Application to a real data set shows that more significant genes were identified by the entropy-based method than by the other two methods.     

Genetic variants in TGF-β pathway are associated with ovarian cancer risk

The transforming growth factor-β (TGF-β) signaling pathway is involved in a diverse array of cellular processes responsible for tumorigenesis. In this case-control study, we applied a pathway-based approach to evaluate single-nucleotide polymorphisms (SNPs) in the TGF-β signaling pathway as predictors of ovarian cancer risk. We systematically genotyped 218 SNPs from 21 genes in the TGF-β signaling pathway in 417 ovarian cancer cases and 417 matched control subjects. We analyzed the associations of these SNPs with ovarian cancer risk, performed haplotype analysis and identified potential cumulative effects of genetic variants. We also performed analysis to identify higher-order gene-gene interactions influencing ovarian cancer risk. Individual SNP analysis showed that the most significant SNP was SMAD6: rs4147407, with an adjusted odds ratio (OR) of 1.60 (95% confidence interval [CI], 1.14–2.24, P = 0.0066). Cumulative genotype analysis of 13 SNPs with significant main effects exhibited a clear dose-response trend of escalating risk with increasing number of unfavorable genotypes. In gene-based analysis, SMAD6 was identified as the most significant gene associated with ovarian cancer risk. Haplotype analysis further revealed that two haplotype blocks within SMAD6 were significantly associated with decreased ovarian cancer risk, as compared to the most common haplotype. Gene-gene interaction analysis further categorized the study population into subgroups with different ovarian cancer risk. Our findings suggest that genetic variants in the TGF-β signaling pathway are associated with ovarian cancer risk and may facilitate the identification of high-risk subgroups in the general population.     

Multiple pathway-based genetic variations associated with tobacco related multiple primary neoplasms

Background

In order to elucidate a combination of genetic alterations that drive tobacco carcinogenesis we have explored a unique model system and analytical method for an unbiased qualitative and quantitative assessment of gene-gene and gene-environment interactions. The objective of this case control study was to assess genetic predisposition in a biologically enriched clinical model system of tobacco related cancers (TRC), occurring as Multiple Primary Neoplasms (MPN).

Methods

Genotyping of 21 candidate Single Nucleotide Polymorphisms (SNP) from major metabolic pathways was performed in a cohort of 151 MPN cases and 210 cancer-free controls. Statistical analysis using logistic regression and Multifactor Dimensionality Reduction (MDR) analysis was performed for studying higher order interactions among various SNPs and tobacco habit.

Results

Increased risk association was observed for patients with at least one TRC in the upper aero digestive tract (UADT) for variations in SULT1A1 Arg²¹³His, mEH Tyr¹¹³His, hOGG1 Ser³²⁶Cys, XRCC1 Arg²⁸⁰His and BRCA2 Asn³⁷²His. Gene - environment interactions were assessed using MDR analysis. The overall best model by MDR was tobacco habit/p53(Arg/Arg)/XRCC1(Arg³⁹⁹His)/mEH(Tyr¹¹³His) that had highest Cross Validation Consistency (8.3) and test accuracy (0.69). This model also showed significant association using logistic regression analysis.

Conclusion

This is the first Indian study on a multipathway based approach to study genetic susceptibility to cancer in tobacco associated MPN. This approach could assist in planning additional studies for comprehensive understanding of tobacco carcinogenesis.     

Weighted Interaction SNP Hub (WISH) network method for building genetic networks for complex diseases and traits using whole genome genotype data

Background

High-throughput genotype (HTG) data has been used primarily in genome-wide association (GWA) studies; however, GWA results explain only a limited part of the complete genetic variation of traits. In systems genetics, network approaches have been shown to be able to identify pathways and their underlying causal genes to unravel the biological and genetic background of complex diseases and traits, e.g., the Weighted Gene Co-expression Network Analysis (WGCNA) method based on microarray gene expression data. The main objective of this study was to develop a scale-free weighted genetic interaction network method using whole genome HTG data in order to detect biologically relevant pathways and potential genetic biomarkers for complex diseases and traits.

Results

We developed the Weighted Interaction SNP Hub (WISH) network method that uses HTG data to detect genome-wide interactions between single nucleotide polymorphism (SNPs) and its relationship with complex traits. Data dimensionality reduction was achieved by selecting SNPs based on its: 1) degree of genome-wide significance and 2) degree of genetic variation in a population. Network construction was based on pairwise Pearson's correlation between SNP genotypes or the epistatic interaction effect between SNP pairs. To identify modules the Topological Overlap Measure (TOM) was calculated, reflecting the degree of overlap in shared neighbours between SNP pairs. Modules, clusters of highly interconnected SNPs, were defined using a tree-cutting algorithm on the SNP dendrogram created from the dissimilarity TOM (1-TOM). Modules were selected for functional annotation based on their association with the trait of interest, defined by the Genome-wide Module Association Test (GMAT). We successfully tested the established WISH network method using simulated and real SNP interaction data and GWA study results for carcass weight in a pig resource population; this resulted in detecting modules and key functional and biological pathways related to carcass weight.

Conclusions

We developed the WISH network method which is a novel 'systems genetics' approach to study genetic networks underlying complex trait variation. The WISH network method reduces data dimensionality and statistical complexity in associating genotypes with phenotypes in GWA studies and enables researchers to identify biologically relevant pathways and potential genetic biomarkers for any complex trait of interest.

     

Association of single nucleotide polymorphisms with form traits in three New Zealand populations of radiata pine in the presence of genotype by environment interactions

Associations of single nucleotide polymorphisms (SNP) with the form traits, branch cluster frequency and stem straightness, were studied in three radiata pine (Pinus radiata D. Don) breeding populations. For branch cluster frequency, high genotype by environment (G×E) interactions were found between two sites of the POP1 trial series (Tarawera and Glenledi) and between two sites of the POP2 trial series (Tarawera and Woodhill). For stem straightness, high G×E interactions were found between two sites (Tarawera and Woodhill) in both the POP2 and POP3 trial series. Thirty-two and 26 SNPs were associated with branch cluster frequency and stem straightness respectively, across a number of sites and trial series in a genetic model assuming heterogeneous additive genetic variances across sites. These SNPs had significant associations with branch cluster frequency or stem straightness either in one site only, across sites within a trial series or across different trial series. Seven SNPs had significant associations with branch cluster frequency and seven SNPs with stem straightness across two of three trial series. One SNP had significant association with stem straightness across all three trial series. Strong G×E interactions resulted in fewer significant SNPs in common across the sites among which there were large effect differences. The percentage of genetic variance explained by SNPs showing significant associations ranged from 0.23 to 8.76 % for branch cluster frequency and from 0.37 to 12.75 % for stem straightness. This study showed that SNP effects for a trait change across environments if G×E interactions exist for that trait.     

Preventive SNP–SNP interactions in the mitochondrial displacement loop (D-loop) from chronic dialysis patients

Chronic dialysis association study involving individual single nucleotide polymorphisms (SNPs) in the mitochondrial displacement loop (D-loop) has previously been reported. However, possible SNP–SNP interactions for SNPs in the D-loop which could be associated with a reduced risk for chronic dialysis were not investigated. The purpose of this study was to propose an effective algorithm to identify protective SNP–SNP interactions in the D-loop from chronic dialysis patients. We introduce ISGA that uses an initialization strategy for genetic algorithms (GA) to improve the computational analysis for protective SNP–SNP interactions. ISGA generates genotype patterns with combined SNPs (SNP barcodes) for chronic dialysis. Using our previously reported 77 SNPs in the D-loop, the algorithm-generated protective SNP barcodes for chronic dialysis were evaluated. ISGA provides the SNP barcodes with the maximum frequency differences of occurrence between the cases and controls. The identified SNP barcodes with the lowest odds ratio (OR) values were regarded as the best preventive SNP barcodes against chronic dialysis. The best ISGA-generated SNP barcodes (two to nine SNPs) are more closely associated with the prevention of chronic dialysis when more SNPs are chosen (OR = 0.64 to 0.32; 95% confidence interval = 0.882 to 0.198). The cumulative effects of SNP–SNP interactions were more dominant in ISGA rather than in GA without the initialization strategy. We provide a fast identification of chronic dialysis-associated protective SNP barcodes and demonstrate that the SNP–SNP interactions may have a cumulative effect on prediction for chronic dialysis.     

Research on Single Nucleotide Polymorphisms Interaction Detection from Network Perspective

Single Nucleotide Polymorphisms (SNPs) found in Genome-Wide Association Study (GWAS) mainly influence the susceptibility of complex diseases, but they still could not comprehensively explain the relationships between mutations and diseases. Interactions between SNPs are considered so important for deeply understanding of those relationships that several strategies have been proposed to explore such interactions. However, part of those methods perform poorly when marginal effects of disease loci are weak or absent, others may lack of considering high-order SNPs interactions, few methods have achieved the requirements in both performance and accuracy. Considering the above reasons, not only low-order, but also high-order SNP interactions as well as main-effect SNPs, should be taken into account in detection methods under an acceptable computational complexity. In this paper, a new pairwise (or low-order) interaction detection method IG (Interaction Gain) is introduced, in which disease models are not required and parallel computing is utilized. Furthermore, high-order SNP interactions were proposed to be detected by finding closely connected function modules of the network constructed from IG detection results. Tested by a wide range of simulated datasets and four WTCCC real datasets, the proposed methods accurately detected both low-order and high-order SNP interactions as well as disease-associated main-effect SNPS and it surpasses all competitors in performances. The research will advance complex diseases research by providing more reliable SNP interactions.     

A simple and computationally efficient sampling approach to covariate adjustment for multifactor dimensionality reduction analysis of epistasis

Epistasis or gene-gene interaction is a fundamental component of the genetic architecture of complex traits such as disease susceptibility. Multifactor dimensionality reduction (MDR) was developed as a nonparametric and model-free method to detect epistasis when there are no significant marginal genetic effects. However, in many studies of complex disease, other covariates like age of onset and smoking status could have a strong main effect and may potentially interfere with MDR's ability to achieve its goal. In this paper, we present a simple and computationally efficient sampling method to adjust for covariate effects in MDR. We use simulation to show that after adjustment, MDR has sufficient power to detect true gene-gene interactions. We also compare our method with the state-of-art technique in covariate adjustment. The results suggest that our proposed method performs similarly, but is more computationally efficient. We then apply this new method to an analysis of a population-based bladder cancer study in New Hampshire.     

Interaction between SNPs in the RXRA and near ANGPTL3 gene region inhibits apoB reduction after statin-fenofibric acid therapy in individuals with mixed dyslipidemia

The mixed dyslipidemia phenotype is characterized by elevated triglycerides (TG), low HDL cholesterol (HDL-C), increased ApoB levels, and premature coronary atherosclerosis. Fibrate-statin combination therapy reduces ApoB levels and coronary events in the mixed dyslipidemia population. We sought to identify gene-gene interactions that affect ApoB response to statin-fenofibric acid therapy in the mixed dyslipidemia population. Using a predefined subset of single-nucleotide polymorphisms (SNPs) that were previously associated with TG, VLDL, or HDL-C, we applied gene-gene interaction testing in a randomized, double-blind, clinical trial examining the response to fenofibric acid (FNA) and its combination with statin in 1,865 individuals with mixed dyslipidemia. Of 11,783 possible SNP pairs examined, we detected a single significant interaction between rs12130333, located within the ANGPTL3 gene region, and rs4240705, within the RXRA gene, on ApoB reduction after statin-FNA therapy (P = 4.0 × 10⁻⁶). ApoB response to therapy gradually reduced with the increasing number of T alleles in the rs12130333 but only in the presence of the GG genotype of rs4240705. Individuals doubly homozygous for the minor alleles at rs12130333 and rs4240705 showed a paradoxical increase of 1.8% in ApoB levels after FNA-statin combination therapy. No gene-gene interaction was identified other than an interaction between SNPs in the ANGPTL3 and RXRA regions, which results in the inhibition of ApoB reduction in response to statin-FNA therapy. Further study is required to examine the clinical applicability of this genetic interaction and its effect on coronary events.     

Powerful multilocus tests of genetic association in the presence of gene-gene and gene-environment interactions

In modern genetic epidemiology studies, the association between the disease and a genomic region, such as a candidate gene, is often investigated using multiple SNPs. We propose a multilocus test of genetic association that can account for genetic effects that might be modified by variants in other genes or by environmental factors. We consider use of the venerable and parsimonious Tukey's 1-degree-of-freedom model of interaction, which is natural when individual SNPs within a gene are associated with disease through a common biological mechanism; in contrast, many standard regression models are designed as if each SNP has unique functional significance. On the basis of Tukey's model, we propose a novel but computationally simple generalized test of association that can simultaneously capture both the main effects of the variants within a genomic region and their interactions with the variants in another region or with an environmental exposure. We compared performance of our method with that of two standard tests of association, one ignoring gene-gene/gene-environment interactions and the other based on a saturated model of interactions. We demonstrate major power advantages of our method both in analysis of data from a case-control study of the association between colorectal adenoma and DNA variants in the NAT2 genomic region, which are well known to be related to a common biological phenotype, and under different models of gene-gene interactions with use of simulated data.     

Gene–Gene and Gene-Sex Epistatic Interactions of MiR146a,IRF5, IKZF1, ETS1 and IL21 in Systemic Lupus Erythematosus

Several confirmed genetic susceptibility loci involved in the interferon signaling and Th17/B cell response for SLE in Chinese Han populations have been described. Available data also indicate that sex-specific genetic differences contribute to SLE susceptibility. The aim of this study was to test for gene–gene/gene-sex epistasis (interactions) in these known lupus susceptibility loci. Six single-nucleotide polymorphisms (SNPs) in MiR146a, IRF5, IKZF1, ETS1 and IL21 were genotyped by Sequenom MassArray system. A total of 1,825 subjects (858 SLE patients and 967 controls) were included in the final analysis. Epistasis was tested by additive model, multiplicative model and multifactor dimensionality reduction (MDR) method. Additive interaction analysis revealed interactions between IRF5 and IKZF1 (OR 2.26, 95% CI 1.48–3.44 [P = 1.21×10⁴]). A similar tendency was also observed between IL21 and ETS1 by parametric methods. In addition, multiple high dimensional gene-gene or gene-sex interactions (three-and four-way) were identified by MDR analysis. Our study identified novel gene–gene/gene-sex interactions in lupus. Furthermore, these findings highlight sex, interferon pathway, and Th17/B cells as important contributors to the pathogenesis of SLE.