- Open Access
Simulating gene-gene and gene-environment interactions in complex diseases: Gene-Environment iNteraction Simulator 2
- Michele Pinelli†1, 2Email author,
- Giovanni Scala†1, 3,
- Roberto Amato1, 3, 4,
- Sergio Cocozza1, 2 and
- Gennaro Miele1, 3, 4
© Pinelli et al.; licensee BioMed Central Ltd. 2012
- Received: 20 December 2011
- Accepted: 10 May 2012
- Published: 14 June 2012
The analysis of complex diseases is an important problem in human genetics. Because multifactoriality is expected to play a pivotal role, many studies are currently focused on collecting information on the genetic and environmental factors that potentially influence these diseases. However, there is still a lack of efficient and thoroughly tested statistical models that can be used to identify implicated features and their interactions. Simulations using large biologically realistic data sets with known gene-gene and gene-environment interactions that influence the risk of a complex disease are a convenient and useful way to assess the performance of statistical methods.
The Gene-Environment iNteraction Simulator 2 (GENS2) simulates interactions among two genetic and one environmental factor and also allows for epistatic interactions. GENS2 is based on data with realistic patterns of linkage disequilibrium, and imposes no limitations either on the number of individuals to be simulated or on number of non-predisposing genetic/environmental factors to be considered. The GENS2 tool is able to simulate gene-environment and gene-gene interactions. To make the Simulator more intuitive, the input parameters are expressed as standard epidemiological quantities. GENS2 is written in Python language and takes advantage of operators and modules provided by the simuPOP simulation environment. It can be used through a graphical or a command-line interface and is freely available from http://sourceforge.net/projects/gensim. The software is released under the GNU General Public License version 3.0.
Data produced by GENS2 can be used as a benchmark for evaluating statistical tools designed for the identification of gene-gene and gene-environment interactions.
- Gene-environment interaction
- Computer simulation
- Complex disease
- Genome-wide association study
- Linkage disequilibrium
Most of the common human diseases with high mortality rates (such as cancer, heart disease, obesity, diabetes, and several common psychiatric and neurological conditions) are classified as complex diseases [1, 2]. By definition, a complex disease is a multifactorial complex trait generally caused by multiple predisposing loci (possibly interacting) and by the exposure to particular environmental factors . Although several genetic and environmental factors have been shown to affect susceptibility to particular complex diseases, the intricate sets of relationships between these factors and disease susceptibility are not yet exhaustively understood. For this reason, typically, the proportion of risk accountable to genetics and environment remains mostly unpredictable . Explanations for their unpredictability might include the occurrence of still unidentified factors and/or the presence of non-linear interactions among already identified factors; for example, some combinations of genetic and environmental factors could have disease risks that are consistently higher than those predicted by a single component.
Gene-environment interactions (G×E) are expected to influence complex phenotypes, for example, disease risk. Hence individuals with predisposing genetics are more likely to develop a disease when exposed to a damaging environment than individuals, exposed to the same environment, without predisposing genetics [5, 6]. The role of G×E is so relevant that it is generally accepted that neglecting them can lead to an underestimation of disease risk, and may explain some of the inconsistencies in replications in different studies .
Complex phenotypes are regulated by pathways and biochemical mechanisms that involve many genetic products. Hence, in addition to interactions among genes and environment, interactions among different genetic loci (G×G) can also influence disease risk. In particular, G×G are defined as epistatic when the allelic variations of one gene alters the effect of variations of another gene . Epistasis has been identified in human diseases [9, 10], and its role in public health has been highlighted .
Surprisingly, despite the general agreement on the relevance of G×E and G×G for correct disease risk estimations, only a few epidemiological studies have attempted to identify them. Indeed, studying the complex interactions among risk factors is a daunting task that requires large data sets and specific research designs. Furthermore, the best statistical method for the identification of G×G and G×E in case-control data sets [11, 12] is still widely debated. The performance of statistical methods that are used for the identification of G×G and G×E are typically influenced by many factors: sample size, number of involved factors, type of interaction, model of inheritance, allelic frequencies, distributions of the environmental factors, and relative strength of different factors affecting disease risk. Unfortunately, only a few of these features are generally assessable in real populations. A further limitation of the epidemiological studies that have been performed so far, is the limited knowledge about the impact of linkage disequilibrium (LD) on association statistics in the presence of G×G and G×E. It has been demonstrated that when G×G occur and the assayed SNPs are not the functional ones but SNPs that are in LD with them, common statistics like r2 are generally inappropriate and tend to lead to an over/underestimation of disease risk .
a Multi-Logistic Model (MLM) that can model any type of G×G and G×E,
a mathematical approach (Knowledge Aided Parameterization System, KAPS) that can translate biological and epidemiological information to MLM parameters, and
GENS (Gene Environment iNteraction Simulator), a software that produces simulated data sets.
Using that approach interactions between one genetic and one environmental factor only could be simulated; therefore, it was not possible to account for epistatic G×G. Moreover, all simulated loci were considered to be independent and thus it was not possible to account for LD .
In the present paper, we describe an extension of the previous model that overcomes such limitations using a new strategy that simulates up to two-genes×one-environment interactions with the possible inclusion of epistasis. Importantly, the present algorithm can be easily extended to manage more than two genetic and one environmental factors. However, to simplify the design of biologically meaningful interactions, we limited the number of features (see the Discussion section for details). Furthermore, the inclusion of two genetic factors (with epistatic interaction) that in turn interact with a continuous environmental factor heavily increased of the complexity of the model. Indeed, statistical methods that can deal with even two genetic factors are still far from being functionally useful for real, large data sets . To provide a realistic genetic background to the simulated populations, we implemented our extended model as a module which can be integrated with simuPOP, a forward-time populations simulator that reproduces realistic demographic and evolutionary features .
The epidemiologic parameters that were used for the sample simulation
1) Starting data (Hap Map)
Chromosomes, or chromosome regions, or markers and marker distance.
The genomic regions containing the loci that will be simulated
The starting frequency and linkage data to be used in the sumulation
2) Simulation of sample’s genetic data
DPLs (Disease Predisposing Loci)
Loci that will influence the disease risk.
Target allelic frequency
Final allelic frequencies at the end of simuPOP simulation
Final sample size
Number of individuals that compose the population by simuPOP
simuPOP generated sample
Sample data generated with simuPOP
The expected disease prevalence in the whole sample
Environmental factor distribution
Distribution of the environmental exposure in the whole sample
Environmental factor OR
Odds ratio associated with one-unit-increase of the environmental exposure
Noisy Environmental variables
As many as desired confounding environmental exposures not associated with the disease risk (gaussian, binomial or uniform distributed)
These are the same DPLs as selected in the simuPOP simulation
High risk alleles
The allele, for each DPL, associated with the highest disease risk
DPLs genotypic RR
The relative risk of the high risk homozygote versus low risk homozygote, for each DPL
The relationship of the risk associated with the heterozygote with that associated with the homozygotes (recessive, dominant, codominant)
Epistasis model (G×G)
Percent increase of the risk associated with each combined genotype
Gene Environment interaction
One of the four predefined interaction models
Generation of the synthetic data set
The generation of the starting sample is carried out by a series of simuPOP scripts  that
download phased genomic data from the HapMap public database ,
select a subset of SNPs or entire genomic regions, and
let the population evolve until it reaches the desired size and frequencies for some disease predisposing loci (DPLs).
To obtain a synthetic data set, simuPOP drives a forward-time simulation to obtain a population that is composed of the desired number of individuals and genotypic frequencies for all the markers. The use of this simulator helps to retain genetic realism in the final population, in particular with respect to the patterns of LD (for a detailed description of this process, please see ). When the genetic information for each individual has been obtained, GENS2 assigns environmental exposures following a user defined Gaussian distribution for the disease environmental variable, and several other user defined distributions (Gaussian, Uniform or Binomial) for a number of other environmental variables not related with the disease (environmental-confusing variables).
Definition of the penetrance model
The second branch of the simulation procedure (right side of Figure 1) is implemented in the Knowledge Aided Parameterization System 2 (KAPS2) subsystem which accepts the input of some standard epidemiological measures for the relevant features listed here:
the expected prevalence of the disease in the sample,
the id in the input data set of one or two DPLs,
the allelic frequencies of DPLs (calculated automatically from the input population),
the effect on disease risk of each DPL in terms of the relative risk of the high risk homozygote compared with the other homozygote,
the dominance relation of each DPL (W), expressed as a number in the interval [0−1], representing the dominance relation (W=0 dominant, W=1 recessive, 0<W<1 co-dominant), and
the distribution parameters and the effect of the environmental factor on disease risk, expressed as odds ratio (OR) of the risk related to one-unit increase in the exposure.
G×E models are chosen from a set of four predefined models, two models of interaction between DPLs and the environment, and two special models in which there is no gene-environment interaction but in which only one genetic or environmental factor contributes to the disease risk (see Table 2).
Predefined gene-environment interaction models in GENS2
Disease risk depends only on the genetics of an individual
Environmental Model (ENV)
Disease risk depends only on environmental exposure of an individual
Gene Environment interaction Model (GEM)
The genetics modifies the effect of the environment in modulating the disease risk
The effects of environment and genetics are independent and sums in modulating the disease risk
G×G models (epistasis) are accepted in the form of percentage variations on the risk associated with a maximum number of three (out of the possible nine) combined genotypes.
Consequently, for each combined genotype, KAPS2 computes the coefficients of a penetrance function of the environmental exposure that is associated with the combined genotype in the MLM. In this step G×E are also modeled; the Additive model (ADD) assumes that combined genotypes with higher penetrance have a higher basal disease risk, while the risk associated with the environmental factor is just added. On the other hand, in the modulative model (GEM), combined genotypes with higher penetrance have the same basal risk although they are more ’sensitive’ to the effect of the environment (see the Methods section).
Disease risk of an individual
In the final step the two branches of the procedure (Figure 1) join. Once the genetics and the environmental exposure for each individual are given, its disease risk can be obtained by evaluating the penetrance function associated with its combined genotype. The risk is then used to assign a disease status using a random process.
To create simulated populations, we employed an existing tool, SimuPOP, and the implementation of the above described algorithm. Using SimuPOP it is possible to drive a forward-time simulation that results in a population composed by the desired number of individuals and having specified genotype frequencies for a set of markers. To be usable in GENS2, populations should be created in SimuPOP as described previously . This procedure generates a data set that contains genotypic data as a set of bi-allelic loci, for each individual in the population.
GENS2 accepts as input a population produced by SimuPOP and the ids of the DPLs chosen from those present in the data set. For each DPL, the penetrance values are provided as relative risk (RR) and dominance model (see ). RR is calculated as the risk ratio of the high over the low risk homozygotes. For each combination of DPL genotypes, the percentage variation of the penetrance value if epistasis occurs can be assigned. GENS2 assigns environmental factor exposure on the basis of a user-defined Gaussian distribution. Any remaining non-implicated loci can be considered as background noise. Similarly, GENS2 can produce an arbitrary number of non-associated random environmental exposures, following a Gaussian, uniform or binomial distribution, that act as further background noise.
On the basis of the selected type of G×G and G×E, GENS2 calculates the coefficients of the MLM as described in the Method section.
For each individual, GENS2 assigns the disease status (affected or unaffected) on the basis of its disease risk by applying the MLM and using a random process.
The main output of the software can be either a single file or several files for a set of subpopulations of a given size produced by means of a subsampling procedure. Subsampling allows bootstrapping procedures to be executed on data sets produced with the same features. The output of GENS2 is in the form of a table in which each row represents an individual and the columns contain, from the left to the right, disease status, gender, environmental exposures and genotypes for each individual.
Two possible formats for the genetics output are available: phased haplotypes or genotypes. In both output formats the initial columns are identical to those described above; however, they differ in the way the genetic information for each individual is represented. In the phased haplotype format, there are two columns for each SNP that report the allele status (either A, C, T or G) on each chromosome. In the genotype format, each SNP is represented by one number (1, 2 or 3), where 2 represents the heterozygote and 3 represents, for DLPs, the high risk homozygote or, for all the other SNPs, the lower frequency homozygote.
In addition to the main output file, GENS2 outputs a log file that contains an extensive report of all the intermediate steps in the procedure and the values used to obtain the populations. Optionally, a file containing the ID and genomic position of the SNPs in the data set can be returned.
GENS2 is designed to be used either from the command line as a Python script, or through a graphical user interface, similar to a wizard, that prompts the user in the specification of all required parameters [see Additional file 1].
Overall, the computational time complexity of the simulation procedure depends by both simuPOP and GENS2. For GENS2, the procedure is dominated by the assignment of the disease status to all individuals in the population. Indeed, after the KAPS2 module has performed the translation of user provided parameters into MLM parameters in bounded constant time, the time complexity becomes linear in the number of individuals and the number of represented variables (genotypes and environmental exposures) for each individual in the simulated population. On the other hand, the amount of time required to perform a simulation with simuPOP depends on the size of the simulation and scales roughly linearly with the number of markers and individuals that are used . GENS2 is written in Python and is completely open-source. The software is freely available from .
Here we describe a method based on the MLM to simulate two genetic and one environmental factors interacting in the determination of a disease risk. The method is implemented in GENS2, a software that is freely available.
To test populations produced by GENS2, we performed a set of analyses on some representative populations. The aim was to emulate a case in which GENS2 was used to assess the performances of a feature selection method. In particular, all the analyses were performed using a logistic regression (glm function in R) with a different model for each test and considering the status as the dependent variable.
The first test was a single-marker analysis on a population of 1,000 cases and 1,000 controls with two DPLs in two distinct genomic regions, with no epistasis and an additive G×E model. The association of each marker with the status was tested using logistic regression analyses with model: disease risk = genetic factor + environmental factor. As expected, the most significant associations were those of DPLs [see Additional file 2]. The result showed that the environmental variable was also associated with the disease (p<10−6). Furthermore, non-causative markers in LD with the two DPLs also showed a significant association that was roughly proportional to the value of r2 with the DPLs.
The second test was similar to the first, except that 10,000 cases and 10,000 controls and a modulative G×E model for the DPLs were used. For this test, the logistic regression was used by considering both an additive model (disease risk = genetic factor + environmental factor) and a multiplicative model (disease risk = genetic factor * environmental factor). None of the markers, when tested by additive model, reached a Bonferroni corrected significance level [see Additional file 3, middle panel]. Conversely, DPLs were found to be significant when the multiplicative model was explicitly considered. Non-causative markers in LD with them were also found to be significant. Notably, this more complex model required a 10-fold increase in the sample size to achieve the same significance level as the previous test.
The model described here can handle, in principle, any number of DPLs and environmental variables. However, we chose to limit the implementation to a relatively small number of factors (two genetic and one environmental) so that setting up the model does not become too complicated for the user. In this way, we reached a balance between the complexity of the represented phenomena and simplicity in the definition of the model. Moreover, the best strategy to identify even simple interactions as single G×G and G×E with binary environmental variables it is still debated (for an example of the debate, see the report on the 2009 Genetic Analysis Workshop [11, 12, 17]). For this reason, we believe that a set of simulated populations in which all features are known provides an important tool for the identification of the best strategy to identify and study G×G and G×E.
Several methods simulating genetic data have been proposed, many of them also handle complex LD patterns and polygenic traits [22, 23]. For example, HAPGEN2  simulates multiple polymorphic loci that are in LD and in the same chromosomal region. HAPGEN2 can simulate G×G, including epistasis, between two loci; however, the available models are limited to a predefined set. Moreover, this program cannot handle G×E and the number of cases and controls that are produced cannot be controlled. Another tool, gs, similarly simulates multiple loci with a realistic pattern of LD; it can modulate a quantitative trait (as disease risk) and can also manage G×G and G×E. In gs the user can define G×G between two loci in two ways: one, by providing a penetrance matrix for combined genotypes or, two, by selecting a penetrance matrix from a predefined set of G×G models. However, both these approaches have some limitations. Although the first allows great control over the penetrance of each genotype, it easily leads to a loosening of control on marginal effects, making the replication of real populations difficult. The second approach, on the contrary, is too restrictive and does not allow any possible interaction to be simulated. For simulation of G×G, gs allows the user to input a list of rules regarding specific combinations of genotypes and levels of environmental values, and the corresponding risk levels. Again, this approach makes it very difficult to control the overall characteristics of populations in terms of marginal effects of genetic and environmental factors. The approaches described above can simulate complex interactions by loosening control on overall population characteristics or, alternatively, can keep the overall population characteristics under control by limiting the freedom to specify the interactions. Another strategy to simulate G×G and G×E is to manually write functions or sets of rules that associate each combination of genes and environmental factors to a risk value. Although this strategy provides more freedom, it is very difficult to set up when the control of marginal effects of single factors is desired. We believe that GENS2 provides a better balance between the freedom to define possible interactions among factors and the control of the overall population characteristics.
The simulated populations produced with GENS2 can be thought of as a sampling of an ideal infinite population that has the characteristics specified by the user. From this point of view, it is easy to understand that fluctuations of observed values around the expected ones can occur. Among the elements that mostly affect these fluctuations, are sample size, allele frequencies, and penetrance values. In particular, small sample sizes increase the effect of sampling error and thus, as expected, these fluctuations tend to vanish as the sample size is increased [see Additional file 4]. This property can be used to test statistical methods because it allows the user to assess how well the methods perform as population and sample features change; conversely, it helps to assess the statistical power of methods in the presence of population features that can only be presumed.
Although the GENS2 part of the simulation process is reasonably fast, the procedure to simulate large populations using simuPOP takes time to complete. It would be difficult to simulate a large number of samples without a cluster system, unless multiple (small) samples are drawn from the same large population.
GENS2 allows the simulation of gene-gene and gene-environment interactions among two genetic and one environmental factor in relation to the risk to develop a complex disease. It is based on data with a realistic pattern of LD and it has no limitations either on the number of individuals that can be simulated or on the number of genetic and environmental factors within a simulated data set. Furthermore, a large amount of effort has been channeled into allowing the input of parameters as standard epidemiological quantities so that the software is immediately usable by the biomedical community.
GENS2 provides large biologically realistic data sets with known features that can be used to challenge, and eventually improve, the statistical tools that are designed to identify those interactions.
the genetics can influence the disease risk either directly or by modifying the effect of the environment.
the genetic loci can have independent effects (no epistasis) or can interact in an epistatic manner, and
the DPLs are not in LD.
The Multi-Logistic Model
where and are free parameters determined by the genetic factors that determine the shape of the function.
To simulate a population, the coefficients and of the logistic functions that produce the desired population features have to be determined. This task is performed by the KAPS2 module that finds the corresponding MLM coefficients by considering all the desired population features within a set of biological constraints.
Determination of MLM parameters
where X is the domain of the environmental variable.
Because every logistic function in MLM is characterized by its own parameters, the 3×3 pairs of values ( that satisfy the constraints expressed by equation (2) need to be found.
In general, equation (2) admits infinite solutions. However, the G×E model imposes some constraints on the coefficients. Thus, by fixing the value of one of the coefficients , hereafter denoted as β AB , the number of degrees of freedom of the system can be reduced, drawing one solution from the equation system. By convention, we chose to associate β AB to the genotype with highest risk; it is easy to show that this value corresponds to the natural logarithm of the odds ratio of the risk which is related to the increase of one unit of the environmental exposure. Constraints imposed on the system by each one of the proposed gene environment interactions model are summarized below:
Genetic effect (GEN): and and .
Environmental effect (ENV): and .
Modulative effect (GEM): and .
Additive effect (ADD): and .
When the interaction model, the matrix containing total risk values for each combination of genotypes, namely , and the value for the coefficient β AB have been defined, a set of six transcendent equations can be written with the coefficients of the logistic functions (except β AB ) as the unknown variables; these equations admit exactly one solution .
Notice that the superscript “I” is a reminder that the independent polygenic model has been assumed.
Using an independent polygenic model and a deformation procedure, epistatic interactions among DPLs can be modeled to obtain a matrix (where superscript “E” stands for epistatic) that still complies with constraints (4) and (5). In this approach, epistasis is modeled as a departure from the independent polygenic model via a change (positive or negative) in one or more entries of .
Once the quantities , , and are given, the constraints (9) and (10) define a convex region in R3×3 in which the elements are assignments for the entries of matrix . The specification of an epistatic model is, therefore, made through the definition of an increments matrix that complies with the constraints (9) and (10) and that also satisfies the positivity condition for . It can be difficult for a user to specify such a matrix in a way that does not violate the above constraints. However, the number of entries of that the user has to provide (paying attention to avoiding extreme or off-range values) can be reduced by letting the system find the remaining entries.
More precisely (in the two variables case), given the constraints of Eq.s (9) and (10) from one up to three entries for can be provided following the rule that any pair must lie on the same row or in the same column. If the user correctly provides three values, the system admits only one assignment for unspecified values of ; however, if the user provides less than three values, there are an infinite number of ways to choose the remaining entries of . In such a case, instead of randomly choosing a solution, a solution that maximizes an “objective function” is chosen. The problem of fixing the remaining values of can be represented as a continuous mathematical programming problem with decision variables that are the non-user-provided entries of and whose admissible region can be determined by Eq.s (9), (10) and (6).
An objective function can be used to minimize the variance of the set of ratios computed ∀(g a ,g b ) corresponding to non-user-assigned . Such a function is suitable for use in all situations in which the relationships between existing variables for which the user does not provide increments are to be maintained as far as is possible.
Establishing the disease status
Once the coefficients of the MLM are fixed, the disease risk for each individual in a population can be established by substituting the coefficients associated with the carried genotype into Eq. (1) and then by evaluating the resulting logistic function forthe exposure level of the environmental disease factor. Finally, to assign the disease status to each individual, the disease risk is compared with a random number drawn from a uniform distribution.
Project name: Gene-Environment iNteraction Simulator 2
Project home page: http://sourceforge.net/projects/gensim/
Operating system(s): Platform independent
Programming language: Python
Other requirements: SimuPop, OpenOpt, wxPython (optional)
License: GNU GPLv3
RA is the recipient of a fellowship from the Doctorate of Computational Biology and Bioinformatics, University “Federico II”, Naples, Italy. The funders had no role in the study design, data collection and analysis, decision to publish, or in preparation of the manuscript.
- Group TGCR: New models of collaboration in genome-wide association studies: the Genetic Association Information Network. Nat Genet. 2007, 39 (9): 1045-51. 10.1038/ng2127. [http://view.ncbi.nlm.nih.gov/pubmed/17728769]View ArticleGoogle Scholar
- Lohmueller KE, Pearce CL, Pike M, Lander ES, Hirschhorn JN: Meta-analysis of genetic association studies supports a contribution of common variants to susceptibility to common disease. Nat Genet. 2003, 33: 177-182. 10.1038/ng1071.View ArticlePubMedGoogle Scholar
- Weeks DE, Lathrop GM: Polygenic disease: methods for mapping complex disease traits. Trends Genet. 1995, 11 (12): 513-9. 10.1016/S0168-9525(00)89163-5. [http://view.ncbi.nlm.nih.gov/pubmed/8533170]View ArticlePubMedGoogle Scholar
- Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, Hunter DJ, McCarthy MI, Ramos EM, Cardon LR, Chakravarti A, Cho JH, Guttmacher AE, Kong A, Kruglyak L, Mardis E, Rotimi CN, Slatkin M, Valle D, Whittemore AS, Boehnke M, Clark AG, Eichler EE, Gibson G, Haines JL, Mackay TF, McCarroll SA, Visscher PM: Finding the missing heritability of complex diseases. Nature. 2009, 461: 747-753. 10.1038/nature08494.PubMed CentralView ArticlePubMedGoogle Scholar
- Cambien F, Poirier O, Lecerf L, Evans A, Cambou JP, Arveiler D, Luc G, Bard JM, Bara L, Ricard S: Deletion polymorphism in the gene for angiotensin-converting enzyme is a potent risk factor for myocardial infarction. Nature. 1992, 359: 641-644. 10.1038/359641a0.View ArticlePubMedGoogle Scholar
- Jacques PF, Bostom AG, Williams RR, Ellison RC, Eckfeldt JH, Rosenberg IH, Selhub J, Rozen R: Relation between folate status, a common mutation in methylenetetrahydrofolate reductase, and plasma homocysteine concentrations. Circulation. 1996, 93: 7-9. 10.1161/01.CIR.93.1.7.View ArticlePubMedGoogle Scholar
- Hunter DJ: Gene-environment interactions in human diseases. Nat Rev Genet. 2005, 6 (4): 287-98. [http://view.ncbi.nlm.nih.gov/pubmed/15803198]View ArticlePubMedGoogle Scholar
- Musani SK, Shriner D, Liu N, Feng R, Coffey CS, Yi N, Tiwari HK, Allison DB: Detection of gene x gene interactions in genome-wide association studies of human population data. Hum Hered. 2007, 63: 67-84. 10.1159/000099179.View ArticlePubMedGoogle Scholar
- Small KM, Wagoner LE, Levin AM, Kardia SL, Liggett SB: Synergistic polymorphisms of beta1- and alpha2C-adrenergic receptors and the risk of congestive heart failure. N Engl J Med. 2002, 347: 1135-1142. 10.1056/NEJMoa020803.View ArticlePubMedGoogle Scholar
- Howard TD, Koppelman GH, Xu J, Zheng SL, Postma DS, Meyers DA, Bleecker ER: Gene-gene interaction in asthma: IL4RA and IL13 in a Dutch population with asthma. Am J Hum Genet. 2002, 70: 230-236. 10.1086/338242.PubMed CentralView ArticlePubMedGoogle Scholar
- Engelman CD, Baurley JW, Chiu YF, Joubert BR, Lewinger JP, Maenner MJ, Murcray CE, Shi G, Gauderman WJ: Detecting gene-environment interactions in genome-wide association data. Genet Epidemiol. 2009, 33 (Suppl 1): 68-73.View ArticleGoogle Scholar
- An P, Mukherjee O, Chanda P, Yao L, Engelman CD, Huang CH, Zheng T, Kovac IP, Dube MP, Liang X, Li J, de Andrade M, Culverhouse R, Malzahn D, Manning AK, Clarke GM, Jung J, Province MA: The challenge of detecting epistasis (G x G interactions): Genetic Analysis Workshop 16. Genet Epidemiol. 2009, 33 (Suppl 1): 58-67.View ArticleGoogle Scholar
- Reimherr M, Nicolae DL: You’ve gotta be lucky: Coverage and the elusive gene-gene interaction. Ann Hum Genet. 2011, 75: 105-111. 10.1111/j.1469-1809.2010.00615.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Amato R, Pinelli M, D’Andrea D, Miele G, Nicodemi M, Raiconi G, Cocozza S: A novel approach to simulate gene-environment interactions in complex diseases. BMC Bioinformatics. 2010, 11: 8-10.1186/1471-2105-11-8.PubMed CentralView ArticlePubMedGoogle Scholar
- North Shore LIJ Research Institute List of Genetic Analysis Software page. [http://www.nslij-genetics.org/soft/]
- Peng B: Simulating gene-environment interactions in complex human diseases. Genome Med. 2010, 2: 21-10.1186/gm142.PubMed CentralView ArticlePubMedGoogle Scholar
- Moore JH, Ritchie MD: STUDENTJAMA. The challenges of whole-genome approaches to common diseases. JAMA. 2004, 291: 1642-1643. 10.1001/jama.291.13.1642.View ArticlePubMedGoogle Scholar
- Peng B, Amos CI, Kimmel M: Forward-Time Simulations of Human Populations with Complex Diseases. PLoS Genetics. 2007, 3 (3): e47+-[http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1829403&rendertype=abstract]PubMed CentralView ArticlePubMedGoogle Scholar
- Peng B, Amos CI: Forward-time simulation of realistic samples for genome-wide association studies. BMC Bioinformatics. 2010, 11: 442-10.1186/1471-2105-11-442.PubMed CentralView ArticlePubMedGoogle Scholar
- HapMap project home page. [http://hapmap.ncbi.nlm.nih.gov/]
- Gene-Environment iNteraction Simulator 2 Download Page. [http://sourceforge.net/projects/gensim]
- Su Z, Marchini J, Donnelly P: HAPGEN2: simulation of multiple disease SNPs. Bioinformatics. 2011, 27: 2304-2305. 10.1093/bioinformatics/btr341.PubMed CentralView ArticlePubMedGoogle Scholar
- Li J, Chen Y: Generating samples for association studies based on HapMap data. BMC Bioinformatics. 2008, 9: 44-10.1186/1471-2105-9-44.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.