 Technical Advance
 Open access
 Published:
Modified entropybased procedure detects genegeneinteractions in unconventional genetic models
BMC Medical Genomics volume 13, Article number: 65 (2020)
Abstract
Background
Since it is assumed that genetic interactions play an important role in understanding the mechanisms of complex diseases, different statistical approaches have been suggested in recent years for this task. One interesting approach is the entropybased IGENT method by Kwon et al. that promises an efficient detection of main effects and interaction effects simultaneously. However, a modification is required if the aim is to only detect interaction effects.
Methods
Based on the IGENT method, we present a modification that leads to a conditional mutual information based approach under the condition of linkage equilibrium. The modified estimator is investigated in a comprehensive simulation based on five genetic interaction models and applied to real data from the genomewide association study by the North American Rheumatoid Arthritis Consortium (NARAC).
Results
The presented modification of IGENT controls the type I error in all simulated constellations. Furthermore, it provides high power for detecting pure interactions specifically on unconventional genetic models both in simulation and real data.
Conclusions
The proposed method uses the IGENT software, which is free available, simple and fast, and detects pure interactions on unconventional genetic models. Our results demonstrate that this modification is an attractive complement to established analysis methods.
Background
It is generally assumed that genetic interactions play an important role in understanding the mechanisms of complex diseases such as coronary heart disease, Alzheimer’s disease, breast cancer or diabetes [1]. In the statistical sense, interaction refers to a situation in which the effect of one factor depends on the values of another factor on a given scale. In our case, genetic interactions come in two flavours. Firstly, interactions between genetic loci, usually termed genegene interactions or epistasis, occur, as described for rheumatoid arthritis (RA) [2–4]. Specifically, Liu et al. [2] described interactions between the locations DQA2 and DQB2 in the HLA region on chromosome 6. Secondly, geneenvironment interactions denote interactions between genetic and environmental factors [5]. For example, Chandra et al. reported that the interaction between serum cholesterol levels and the sigma4 genotype [6] plays a role in Alzheimer’s disease.
As a starting point, we will in this paper focus on the detection of genegene interactions in a casecontrol setting, although interactions of higher order can also be of interest. For simplification, our description is restricted to the situation where diallelic genetic markers such as single nucleotide polymorphisms (SNPs) are used leading to three possible genotypes. However, the results are generalizable to geneenvironmentinteractions with categorical environmental factors or indeed interactions of any categorical variables.
Generally, a single locus testing strategy is undertaken as the primary analysis in a genomewide association study (GWAS), but this may be unsuitable to detect loci that interact with other variants since the relevant loci may not display effects on their own [7]. A variety of methods exists to detect or control for the presence of genegene interactions [8]. For a binary phenotype, most of them are based on the saturated logistic regression model for interaction [9, 10] or a simplication of it. The regression parameters are penetrances, odds or log odds. As Cordell [8] notes, this procedure implicitly assumes that the scale used for the regression parameters is the scale of interest. The saturated model has nine twolocus genotypes that are modeled by one intercept, four main effects parameters and four interaction parameters with a dummy coding of the genotypes. Although the saturated model is the best fitting one, a model with fewer parameters might be preferable, e.g. because of greater stability. This can be achieved by assuming a specific genetic model and thus estimating, for example, the additive effect of the number of risk alleles at both genetic loci. In this case, only one parameter is estimated for the interaction effect. In this line, the standard software PLINK [11] provides an overall 4 degreeoffreedom (df) test for interaction or a derived 1df test assuming additive effects at both loci.
However, interaction models have been observed in reality that cannot easily be described by regression models for genegeneinteractions without dummy coding. An example is given by Ziegler and König for sporadic breast cancer [10]. Therefore, the restriction to linear models is not optimal. As an alternative, a number of novel techniques have been developed in the last years that are based on different concepts such as rank building (ANOVA technique) [12], datamining with the multifactor dimensionality reduction (MDR) approach [13, 14], machine learning methods with random forests (RF) [14, 15] or with support vector machines (SVM) [16]. Although promising, it is not always clear which interaction effects can be reliably detected by these methods [17].
Another new idea is borrowed from information theory, the entropybased method. This concept is modelfree and measures the uncertainty or disorder in a system and could therefore lend itself to detect interactions for many genotype constellations. This technique is suggested as particularly powerful and, because of the nonlinearity, as better able to capture nonlinear relationships between genetic variants or other variables [18]. Ferrario et al. reviewed different entropybased measures providing information on suggested test statistics, simulations and implementations [18]. Focusing on second order interaction, there are three important concepts, namely conditional mutual information, information gain, and relative information gain. These are based on the following definitions: First, entropy can be defined as a measure for uncertainty in a random variable [19]. Then, mutual information refers to the reduction of uncertainty of one variable conditional on the knowledge of the other variable [18]. Furthermore, mutual information can also be conditioned on a third variable yielding the conditional mutual information (CMI) [18], which has been used for a test statistic by Zuo et al. [20]. Second, the term information gain is defined in different ways: Fan et al. [21] subtract the mutual information of two genetic variants estimated in the cases from the same quantities estimated in the controls. Alternatively, in the method IGENT Kwon et al. [22] subtract the conditional entropy of the phenotype, given two genetic variants, from the entropy of the phenotype. Kwon et al. [22] also work with the socalled relative information gain, which is given as the relation between the information gain and the entropy of the phenotype.
The advantage of the approach by Kwon et al. lies in the freely available and fast implementation that is also called IGENT [18]. The software has been implemented in C++ and is available at http://statgen.snu.ac.kr/software/igent/(different from the information in the paper from 2014 [22]). In contrast, Zuo et al. [20] work with an individual software.
One characteristic of entropybased procedures is that main effects may also present themselves as a deviation from disorder, i.e. the entropy falls and the test reacts. This is an advantage if we search for main effects and interaction effects simultaneously and do not need to distinguish between either. This might be the case in a first step of an analysis in which the focus is on generally learning about variants having an effect on the outcome. Also, this offers the possibility to reduce the variants for a second computationally more intensive step. However, it should be noted that if we are interested in interaction effects only, a main effect of one or both genetic variants without an interaction will lead to a false positive result. Zuo et al. [20] state that the CMI concept achieves better or comparable control of the false positive error, compared to four previously proposed modelfree metrics [20].
In the following, we are only interested in interaction effects, so it is necessary to eliminate main effects, without diminishing the advantages of the entropy approach. We therefore introduce a modification of IGENT that eliminates the problem of the increased typeIerror in the case of only main effects, but keeps the advantages of the entropy method as far as possible. We illustrate the behavior of the proposed procedure with data simulated for different genotypic models and apply it to the analysis of real data on the genetic background of RA [23, 24]). The same data set was analyzed previously by Liu et al. [2], who utilized a regressionbased approach combined with random forest analyses and Chattopadhyay et al. [3] who worked with three nonparametric scores. Furthermore, most comparisons of interaction methods have so far focused on assessing deviation from additive or multiplicative effects. However, we assume that the strongest advantage of entropybased methods is seen in more unconventional interaction models not following classical genetic models, and we considered these unconventional interaction models in our simulations.
Methods
Entropy and IGENTestimator
Entropy is originally a term from thermodynamics referring to the level of disorder or uncertainty. Information theory has utilized this phenomenon as a measure for the lack of structure in a system [19]. Shannon defines the entropy H of a set of probabilities p_{1},…,p_{n} as \(  \sum p_{i} \log p_{i} \).
In the context of a disease state D depending on the genotypes at two genetic loci, Kwon et al. [22] derive the Information Gain similarly as follows:
First let entropy of the phenotype be written as
with D_{0} and D_{1} denoting the unaffected and affected state, respectively.
The second order entropy (conditional entropy of the disease state on the genotypes G) is then given by the expression
Here, for (i,j=0,1,2), i and j define the genotype at the 1st and 2nd locus. Thus, all possible genetic models are considered.
From that, the Information Gain can be derived as
The estimator (2nd order) then leads to
For this, X_{ijk} are the observations of G_{ij} and D_{k} in N individuals, leading to the mean \(\bar {X} = \sum _{i,j=0}^{2} \sum _{k=0}^{1} X_{ijk}\).
Then, \(\hat P_{ijk} = X_{ijk} / \bar {X}\), \(\hat P_{ij.} = \sum _{k=0}^{1} P_{ijk} \) and \(\hat P_{..k} = \sum _{i,j=0}^{2} {P_{ijk}}\).
Furthermore we define the mutual information as
Finally, the conditional mutual information is given by
Considering two genetic loci with three genotypes each, we construct 3×3 contingency tables which tabulate the values in terms of penetrances or odds for the resulting 9 genotype combinations. In this scenario, the IGENT method estimates an imbalance of the 9 odds, i.e., the deviation for the value of 1 for all odds.
The estimator of IG(DG) (IGENTestimator) asymptotically and approximately follows a gamma distribution under the null hypothesis that genotype combinations and disease states are independent. This null hypothesis is also violated in the case of association with one or both of the genetic variants. Thus, the typeIerror is inflated if this estimator is used as a test for interactions only.
Modification of the iGENTestimator
In the following, we want to utilize the IGENT approach while eliminating the influence of main effects so as to yield a purely interactive effect etimator. For this, we first estimate the main effect of one genetic variant by applying the IGENT approach to a 1st order calculation. Specifically, the Information Gain of 1st order for the first genotype is given by
This can be estimated by
with \(\hat P_{i.k} = \sum _{j=0}^{2} {P_{ijk}}\).
The estimator (1st order) for the second genotype is given analogously.
This leads to the following intuitive modification of IGENT:
All three components can be estimated within the IGENT software.
The estimator of IGmod is
with \( \hat E_{ijk} = \frac { \hat P_{ijk} \hat P_{..k} }{ \hat P_{i.k} \hat P_{.jk} }\) and \( \hat L_{ij.} = \frac { \hat P_{ij.} }{ \hat P_{i..} \hat P_{.j.} }\).
Here, the factor \(\hat E_{ijk}\) takes the value 1 if the two genotypes are conditionally independent, and the factor \(\hat L_{ij.}\) equals 1 under the condition of no correlation between the genetic variants, i.e., no linkage disequilibrium (LD).
Under this condition of no LD, we can simplify the formula for \(\widehat {IGmod}(DG)\) to
with \(\hat P^{k}_{ij} = \frac {\hat P_{ijk}}{\hat P_{..k}}\), \(\hat P^{k}_{i.} = \frac {\hat P_{i.k}}{\hat P_{..k}}\), and \(\hat P^{k}_{.j} = \frac {\hat P_{.jk}}{\hat P_{..k}}\).
This conversion shows that this estimator works with conditional mutual information. Specifically, it estimates the deviation from the conditional independence, and it follows asymptotically and approximately a gamma distribution with shapeparameter 4 and scale parameter \(\frac { 1}{N \ln (2)}\) under the null hypothesis of conditional independence of the genetic variants [25].
Comparison with genoCMI
Recently, Zuo et al. [20] introduced an estimator called GenoCMI that was defined as follows:
Under the condition of no LD, this is identical to our modified IGENTestimator \(\widehat {IGmod}\), except for the basis of the logarithm. GenoCMI follows asymptotically and approximately a χ^{2}(ν)/2N distribution where the degree of freedom ν is 8. This statement is equivalent to the above result of a gamma distribution with shape parameter 4 and scale parameter \(\frac { 1}{N \ln (2)}\). An obvious disadvantage of GenoCMI is that there is no freely available software implementation, whereas IGENT is freely available and efficiently implemented.
Simulation models for genegeneinteractions
The aim of our simulation study was to evaluate the performance of different estimators not only in commonly assumed interaction models but also in more unusual interaction settings. We therefore selected five interaction models, of which the first four models were proposed by Wan et al. [26], the fifth model was based on Ritchie et al. [27] (there model 4) and was generated using the epistasis model discovery method of Moore et al. [28]. Of note, the models display interaction effects but little or no main effects and can be written as 3×3 contingency tables of odds for the first interacting variant with genotypes aa, aA, and AA, and the second interacting variant with genotypes bb, bB, BB, where the minor alleles are denoted by capital letters (see Supplement). The specific values are each determined by a prevalence parameter α and a multiplicative interaction parameter θ and are shown in Table 1 and visualized in the accompanying Figs. 1, 2, 3, 4 and 5 (after conversion to the casecontrol scenario).
Three models (epistasis model, two allele interactionmodel and XOR model) display only two levels of risk, whereas the multiplicative and the no marginmodel a more complicated risk pattern. Specifically, the no marginmodel was selected because it exhibits interaction effects in the absence of any main effects. Furthermore, all marginal odds are equal for this special model.
Data simulation
We simulated data according to the five genotype interaction models using the software GAMETES [29, 30]. Genotypes were generated according to HardyWeinberg proportions, and a range of allele frequencies and patterns of riskgenotype associations was chosen [31].
For the models according to Wan (multiplicative model, epistasis model, two allele interactionmodel and XOR model), three different minor allele frequencies (MAFs) were chosen (0.1, 0.2, and 0.4) leading to twelve odds tables (see Supplement). Setting the prevalence to 0.1 and the heritability to 0.03 for the multiplicative model and to 0.02 for the other models, the prevalence parameter g and the interaction parameter t were determined. The heritability is defined as described by Wan et al. [26]:
For utilization in GAMETES, these odds tables were converted to penetrance tables. The no marginmodel is directly taken from Ritchie et al. [27, 32]. With the MAF for the interacting variants set to 0.25 this model shows no marginal effects.
Each model was simulated with 100 variants, with 2 of them interacting, 1000 replicates (datasets), and 800 cases and 800 controls. The MAFs of the noninteracting SNPs were chosen randomly between 0.05 and 0.5.
For simulations under the null hypothesis, we used the four different interaction models according to Wan et al. [26] with three different MAFs, resulting in 12 models that were simulated with 1000 replicates each. Given that for every model, 100 SNPs were simulated with two interacting SNPs, we had 4949 noninteracting SNP pairs per model. Thus, a total of more than 50 million SNP pairs without interactions were simulated. Simulated data was evaluated at significance thresholds of 5×10^{−2}, 5×10^{−3}, 5×10^{−4}, 5×10^{−5}, and 5×10^{−6}.
In the resulting data sets, we estimated \(\widehat {IGmod}_{0}\) as described above.
Comparison with other approaches
To compare our estimator with previous established approaches, we estimated logistic regression models testing for interaction in an additive genetic model with 1 df as implemented in the module epistasis of the software PLINK [11]. Furthermore, we performed likelihood ratio tests (LRT) comparing a model with 4 parameters (2 additive and 2 dominant terms) with the full model with 8 parameters [33]. Finally, another entropybased approach with the test statistic T_{IG} was used as described by Fan et al. [21].
Essentially, Fan et al. [21] subtract the mutual information of two genetic variants estimated in the cases from the same quantities estimated in the controls.
Submodel classification
To illustrate the underlying genetic models of the simulated data, we utilized one step of the modelbased multifactor dimensionality reduction (MBMDR) algorithm, which is an efficient algorithm to perform multiple testing in epistasis screening [34]. The procedure tabulates the frequencies of cases and controls in the 3×3 genotype combinations and uses a test for association between the trait and the specific genotype combination. The test is performed for every cell of the 3×3 contingency table and denotes the individuals with the genotype combination of the specific cell as having a high risk of being affected (H), or a low risk (L), or not sufficient evidence or information (O). The result is a 3×3 matrix, denoted as HLOmatrix. For illustration, these matrices are shown in Table 1 (right) under the assumption that sample size is large enough to yield sufficient evidence.
Real data
We use the data from the genomewide association study by the North American Rheumatoid Arthritis Consortium (NARAC) that were also analyzed by Liu et al. [2]. The data set comprises genotype data of 2,062 individuals, 868 cases with RA and 1,194 controls, predominantly of Northern European origin. The data had been genotyped on the Illumina 550k platform. After exclusion of monomorphic SNPs and SNPs showing deviation from HWE at p<0.0001, 515,680 SNPs were available for further analysis [24]. Quality control procedures included removing individuals who had a low overall call rate (<95 %) of SNPs [23]. From these data we select the HLAregion on chromosome 6 encompassing 2010 SNPs. This area offers a large number of SNPs of which many are associated with RA, and previous analyses hinted at genegene interactions in this region. To reduce the number of SNPpairs to investigate, we further selected only SNPs overlapping with the 749 SNPs analyzed by Liu et al. [2]. Because of the assumption of no LD, we furthermore eliminated all SNP pairs with a LD of r^{2}>0.01.
Significance thresholds for the modified IGENT interaction evaluation
As described above, the estimator \(\widehat {IGmod}_{0}\) asymptotically and approximately follows a gamma distribution with shape parameter 4 and scale parameter \(\frac { 1}{N \ln (2)}\) under the null hypothesis of conditional independence of the genetic variants [25].
Because of the fact that
the characteristic of the underlying distribution is given only asymptotically and approximately,
the Bonferroni correction is very conservative, and
the test statistic is dependent on allele frequencies and marginal effects,
fixing the shape parameter at 4 is partially very conservative. Thus, we utilize alternative cutoffs to identify relevant pairs of interacting SNPs. For the simulated data, we set the global significance level to α=0.05 and apply a Bonferroni correction to adjust for the number of interactions being tested from the gamma distribution with shape parameter 2 (liberal criterion). For the real data, we also set the global significance level to α=0.05 but apply a Bonferroni correction to adjust for the number of interactions being tested from the gamma distribution with shape parameter 4 (conservative criterion).
Assuming the scale parameter of \(\frac { 1}{N \ln (2)}\), this leads to a cutoff at <0.012836 (see below) for the simulated data and at <0.017331 for the real data.
For the regression analysis we set the global significance level to α=0.05 with Bonferroni correction (for the simulated data based on the number of SNP pairs), which leads to significance levels of 1×10^{−5} for the simulated data and 5×10^{−8} for the real data.
Results
Evaluation of type I error
Firstly, the type I error is evaluated in the simulation data. For the null simulation we took about 59 million SNP pairs like described above. The global significance level of 5% is controlled using the liberal cutoff of 0.012836 (Table 2). This is due to the structure of the null simulation data with different MAFs, but only a small portion of SNPs with main effects. Table 2 shows that the type I error matches the different thresholds of 5×10^{−2}, 5×10^{−3}, 5×10^{−4}, 5×10^{−5}, and 5×10^{−6} under the gamma distribution with shape parameter 2 very well.
Evaluation of power
Figure 6 shows the power to detect the interaction in the simulated data for all interaction models and MAF settings. There are four scenarios that display a limited power of <80%, i.e., the multiplicative model with MAF=0.1, 0.2 and 0.4, and the two allele interactionmodel with MAF=0.1. By comparison, the logistic regression model shows overall lower power to detect the interactions, with scenarios having power <80% including the multiplicative model with MAF=0.4, the epistasis model with MAF=0.1, the two allele interactionmodel as well as the XOR model with all MAFs. The no marginmodel provides satisfying power for both estimators, the proposed estimator and logistic regression. Thus, specifically for the epistasis model at a low MAF and unconventional models at all MAFs, the proposed estimator \(\widehat {IGmod}_{0}\) performs better than the classical logistic regression in detecting interactions. The LRT shows overall comparable power as \(\widehat {IGmod}_{0}\) with better performance in the multiplicative model but lower power in the epistasis and XOR models with low MAF. Finally, the estimator T_{IG} provides good results (power >80%) in the epistasis and XOR models with higher MAFs and for the no marginmodel, but the power is lower than that of the proposed estimator throughout.
Submodel classification
We observed an intermediate power of \(\widehat {IGmod}_{0}\) depending on the specific setting and interaction model. Thus, to get a more detailed impression of the underlying interactions, we classified the frequently occuring identified interaction pairs using HLOmatrices as described above.
Multiplicative model
Table 3 shows the frequencies of the more common submodels occuring in the simulation of the multiplicative model with MAF=0.2 along with the power of \(\widehat {IGmod}_{0}\) and the logistic regression to detect the interaction.
It can be seen that \(\widehat {IGmod}_{0}\) provides high power for mostly complete submodels, in which almost every cell contains enough information to perform the association test. In incomplete submodels, in which not all cells of the HLOmatrix show an effect or the cells contain not enough information to perform the association test, \(\widehat {IGmod}_{0}\) shows more difficulty to detect the interaction. In contrast, the regression model provides a power higher than 80% in all submodels.
Epistasis model
Analogously, Table 4 shows the respective submodel results for the epistasis model with MAF=0.1.
Notably, many results with incomplete submodels were obtained in which one or more cells of the HLOmatrix did not show evidence for association or did not contain enough information for association testing. For example, 41% of the relevant pairs led to HLOmatrices which show association only in four out of nine possible cells. However, \(\widehat {IGmod}_{0}\) provides high power for all submodels, whereas the regression model has more difficulty detecting interactions in incomplete submodels.
Real data
In the analysis of the real data set, we obtain 211 relevant SNP pairs, of which 102 display a model similar to the multiplicative model. Further 90 relevant SNP pairs resemble the epistasis model, and 19 SNP pairs follow more unconventional models. Compared to the regression analysis, our estimator identifies 31 SNP pairs as relevant that are not detected by the regression approach, and these are shown in Tables 5 and 6.
Most of the interaction pairs essentially follow a multiplicative model (multi) or an incomplete multiplicative model (multi incompl). There are seven SNP pairs for which the model resembles the epistatic model (epi), and the remaining are either similar to the XOR model or are difficult to classify (other).
Of note, 39 SNP pairs were detected by the regression approach at a significance level of 1×10^{−9} but not by \(\widehat {IGmod}_{0}\). All of these pairs belong to the categories of multiplicative and epistatic models. They do achieve a relatively high test statistic value of more than 0.01, but do not meet our conservative significance criterion of 0.017331.
Comparing the SNP pairs identified by \(\widehat {IGmod}_{0}\) with the results by Liu et al. [2], we find that all eight SNP pairs reported by Liu et al. [2] are also detected by our proposed estimator (Table 7). Notably, all of these eight pairs belong to either the multiplicative or the epistasis model and have relatively high MAFs (>32%).
Finally, we compared our results with those by Chattopadhyay et al. [3]. Their 20 top ranked 2way interactions within the HLA region contain SNPs with a low MAF, but they were excluded from our evaluation either due to an LD of more than 0.01 or because they were not contained in the set of 749 SNPs analyzed by Liu et al. [2].
Discussion
In this paper we proposed a modification of the IGENT estimator for second order genetic interactions, which exploits the advantages of the entropy methods, but controls the type I error. It is a tool designed specifically to detect various patterns conditional on a baseline model, and foremost to detect interactions. The estimator asymptotically and approximately follows a gamma distribution with a shape parameter that depends on the existence of margin effects. The calculation can be made within the IGENT software, which is simple to apply and very fast.
Our simulations show that the power of our method depends on marginal effects, and the lowest power was observed for the multiplicative model. Instead, the advantages of our modification lie in the possibility to detect interactions for unconventional interaction models (like two allele interactionmodel and XOR model), and for incomplete models (e.g. like last HLOmatrix in Table 4), in which not all genotype combinations are observed with sufficient frequency. Thus, \(\widehat {IGmod}_{0}\) outperforms all other approaches in the XOR and epistasis models especially at low MAFs and is better than all approaches except the LRT in the two allele model, where \(\widehat {IGmod}_{0}\) and the LRT are comparable.
In the analysis of the real data on RA, we find the results of the simulation confirmed. We observe a similar number of SNP pairs detected by one but not the other approach. Again, our modification is more likely to uncover unconventional or incomplete models.
In general, our results confirm the conclusion of Zuo et al. [20] that GenoCMI measures achieve a control of the false positive error in the presence of main effects.
Conclusions
In conclusion, we proposed a modification of the IGENT method, which is a fast and efficient entropybased interaction analysis algorithm [22]. The modification reduces the type I error, so it can easily identify second order genegene interactions on a genomewide scale. The analysis of simulated and real data has shown that, in contrast to classical regression approaches, more unconventional interaction models can be detected with this approach, which makes it an attractive complement to established analysis methods.
Availability of data and materials
Not applicable
Abbreviations
 ANOVA:

Analysis of Variance Technique
 CMI:

Conditional Mutual Information
 GAMETES:

Genetic Architecture Model Emulator for Testing and Evaluating Software
 GenoCMI:

Special Measure for CMI
 GWAS:

Genomewide association study
 HLA:

Human Leukocyte Antigens
 HLO:

High, Low or No Evidence
 HWE:

HardyWeinberg Disequilibrium
 IG:

Information Gain
 IGENT:

Interactions analysis method in Genomewide scale based on ENTtropy
 LD:

Linkage Disequilibrium
 MAF:

Minor Allele Frequency
 MBMDR:

ModelBased Multifactor Dimensionality Reduction
 MDR:

Multifactor Dimensionality Reduction Method
 NARAC:

North American Rheumatoid Arthritis Consortium
 PLINK:

Tool Set for WholeGenome Association and PopulationBased Linkage Analyses
 RA:

Rheumatoid Arthritis
 RF:

Random Forests Method
 SNP:

SingleNucleotide Polymorphism
 SVM:

Support Vector Machine Method
 XOR:

eXclusive OR Logic Disjunction
References
Ritchie MD, Steen KV. The search for genegeneinteractions in genomewide association studies: challenges in abundance of methods, practical considerations, and biological interpretation. Ann Trans Med. 2018; 2018;6(8):157.
Liu C, Ackerman HH. A genomewide screen of genegene interactions for rheumatoid arthritis susceptibility. Hum Genet. 2011; 2011:473–485.
Chattopadhyay AS, et al.Summarizing techniques that combine three nonparametric scores to detect diseaseassociated 2way snpsnp interactions. Gene. 2014; 533:304–312.
Huang CH, et al.Rheumatoid arthritisassociated genegene interaction network for rheumatoid arthritis candidate genes. Int Conf Biomed. 2009; 3(Suppl 7):75.
Aschard H, et al.Challenges and opportunities in genomewide environmental interaction (gwei) studies. Hum Genet. 2012; 131:1591–1613.
Chandra V, Pandav R. Geneenvironment interaction in alzheimer’s disease: A potential role for cholesterol. Neuroepidemiology. 1998; 17:225–232.
Ueki M, Cordell HJ. Improved statistics for genomewide interaction analysys. PLoS Genet. 2012; 8(4):1002625.
Cordell HJ. Epistasis: what it means, what it doesn’t mean, and statistical methods to detect it in humans. Hum Mol Genet. 2002; 11(20):2463–2468.
Gola D, et al.Evaluating predictive biomarkers for a binary outcome with linear logistic regression  practical recommendations for the choice of the model. BioRxiv. 2018; preprint.
Ziegler A, König IR, Pahlke F. A statistical approach to genetic epidemiology: Concepts and applications, with an elearning platform, 2nd edn.Weinheim: Wiley; 2010.
Purcell S, et al.Plink: A tool set for wholegenome association and populationbased linkage analyses. Am J Hum Genet. 2007; 81:559–575.
Brunner E, Puri ML. Nonparametric methods in factorial designs. Stat Pap. 2001; 42:1–52.
Hahn L, Moore J. Ideal discrimination of discrete clinical endpoints using multilocus genotypes. Silico Biol. 2004; 4(2):183–94.
McKinney BA, et al.Machine learning for detecting genegeneinteractions. Appl Bioinform. 2006; 5(2):77–88.
Breiman L. Random forests. Mach Learn. 2001; 45(1):5–32.
Chen SH, et al.A support vector machine approach for detecting genegeneinteraction. Genet Epidemiol. 2008; 32:152–167.
Wright MN, et. al., BMC Bioinformatics. Do little interactions get lost in dark random forests?2016;17–145.
Ferrario PG, König IR. Transferring entropy to the realm of gxg interactions. Brief Bioinform. 2016;1–12.
Shannon CE. A mathematical theory of communication. The Bell System Technical Journal. 1948; 27:379–423623656.
Zuo X, et. al.To control false positives in genegeneinteraction analysis: Two novel conditional entropybased approaches. PloS ONE. 2013; 8(Issue 12):81984.
Fan R, et al.Entropybased information gain approaches to detect and to characterize genegene and geneenvironment interactions/ correlations of complex diseases. Genet Epidemiol. 2011; 35(7):706–721.
Kwon MS, et al.Igent: efficient entropy based algorithm for genomewide genegeneinteraction analysis. BMC Med Genomics. 2014; 7(Suppl 1):6.
MacCluer JW, et al.Genetic analysis workshop 16: Introduction to workshop summaries. Genet Epidemiol. 2009; 33(Suppl 1):1–7.
Kruppa J, et al.Risk estimation and prediction using machinelearning methods. Hum Genet. 2012; 131(19):1639–54.
Goebel B, Dawy Z, Hagenauer J, Mueller JC. An approximation to the distribution of finite sample size mutual information estimates. IEEE International Conference on Communications, 2005. ICC 2005. 2005. 2005; 2:1102–11062.
Wan X, et al.Boost: A fast approach to detecting genegeneinteractions in genomewide casecontrol studies. Am J Hum Genet. 2010; 87:325–340.
Ritchie MD, et al.Power of multifactor dimensionality reduction for detecting genegeneinteractions in the presence of genotyping error, missig data, phenocopy, and genetic heterogeneity. Genet Epidemiol. 2003; 24:150–157.
Moore JH, et al.Application of genetic algorithms to the discovery of complex models for simulation studies in human genetics. Proc Genet Evol Comput Conf. 2002:1150–1155.
Urbanowicz RJ, et al.Predicting the difficulty of pur, strict, epistatic models: metrics for simulate model selection. BioData Min. 2012; 5:15.
Urbanowicz RJ, et al.Gametes: a fast. direct algorithm for generating pure, strict, epistatic models with random archtectures. BioData Min. 2012; 5:16.
Culverhouse RC. A comparison of methods sensitive to interactions with small main effects. Genet Epidemiol. 2012; 36(4):303–311.
Ritchie MD, et al.Genetic heterogeneity is not as threatening as you might think. Genet Epidemiol. 2007; 31:797–800.
Hallgrimsdottir I, et al.A complete classification of epistatic twolocus models. BMC Genetics. 2008; 2008:9–17.
Lishout FV, et al.An efficient algorithm to perform multiple testing in epistasis screening. BMC Bioinformatics. 2013; 14:138.
Acknowledgments
We would like to thank Drs. Christopher I. Amos and Jean W. MacCluer, and Vanessa Olmo for the permission to use the data.
Funding
This work is based on data that was gathered with the support of grants from the National Institutes of Health (NO1AR22263 and RO1AR44422), and the National Arthritis Foundation. The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Author information
Authors and Affiliations
Contributions
The author(s) read and approved the final manuscript
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Additional file 1
Odds table for 13 settings. Values of odds for 4 models with 3 different minor allele frequencies (MAFs) (0.1, 0.2, 0.4 resp. for the interacting SNPs (genotypes)) and 1 model with MAF = 0.25. The minor alleles are denoted by capital letters.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Malten, J., König, I.R. Modified entropybased procedure detects genegeneinteractions in unconventional genetic models. BMC Med Genomics 13, 65 (2020). https://doi.org/10.1186/s1292002007034
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1292002007034