- Research
- Open Access
- Published:

# A model of gene-gene and gene-environment interactions and its implications for targeting environmental interventions by genotype

*Theoretical Biology and Medical Modelling*
**volume 3**, Article number: 35 (2006)

## Abstract

### Background

The potential public health benefits of targeting environmental interventions by genotype depend on the environmental and genetic contributions to the variance of common diseases, and the magnitude of any gene-environment interaction. In the absence of prior knowledge of all risk factors, twin, family and environmental data may help to define the potential limits of these benefits in a given population. However, a general methodology to analyze twin data is required because of the potential importance of gene-gene interactions (epistasis), gene-environment interactions, and conditions that break the 'equal environments' assumption for monozygotic and dizygotic twins.

### Method

A new model for gene-gene and gene-environment interactions is developed that abandons the assumptions of the classical twin study, including Fisher's (1918) assumption that genes act as risk factors for common traits in a manner necessarily dominated by an additive polygenic term. Provided there are no confounders, the model can be used to implement a top-down approach to quantifying the potential utility of genetic prediction and prevention, using twin, family and environmental data. The results describe a solution space for each disease or trait, which may or may not include the classical twin study result. Each point in the solution space corresponds to a different model of genotypic risk and gene-environment interaction.

### Conclusion

The results show that the potential for reducing the incidence of common diseases using environmental interventions targeted by genotype may be limited, except in special cases. The model also confirms that the importance of an individual's genotype in determining their risk of complex diseases tends to be exaggerated by the classical twin studies method, owing to the 'equal environments' assumption and the assumption of no gene-environment interaction. In addition, if phenotypes are genetically robust, because of epistasis, a largely environmental explanation for shared sibling risk is plausible, even if the classical heritability is high. The results therefore highlight the possibility – previously rejected on the basis of twin study results – that inherited genetic variants are important in determining risk only for the relatively rare familial forms of diseases such as breast cancer. If so, genetic models of familial aggregation may be incorrect and the hunt for additional susceptibility genes could be largely fruitless.

## Background

Some geneticists have predicted a genetic revolution in healthcare: involving a future in which individuals take a battery of genetic tests, at birth or later in life, to determine their individual 'genetic susceptibility' to disease [1, 2]. In theory, once the risk of particular combinations of genotype and environmental exposure is known, medical interventions (including lifestyle advice, screening or medication) could then be targeted at high-risk groups or individuals, with the aim of preventing disease [3].

However, there are also many critics of this strategy, who argue that it is likely to be of limited benefit to health [4–8]. One area of debate concerns the proportion of cases of a given common disease that might be avoided by targeting environmental or lifestyle interventions to those at high genotypic risk. Known genetic risk factors have to date shown limited utility in this respect [9]. However, some argue that combinations of multiple genetic risk factors may prove more useful in the future [10].

There are two possible approaches to considering this issue. The 'bottom-up' approach seeks to identify individual genetic and environmental risk factors and their interactions and quantify the risks. However, this approach is limited by the difficulties in establishing the statistical validity of genetic association studies and of quantifying gene-gene and gene-environment interactions: see, for example, [11–14].

A 'top-down' approach instead considers risks at the population level using twin and family studies and data on the importance of environmental factors in determining a trait. However, analysis of twin data is usually limited by the assumptions made in the classical twin study [15], including that: (i) there are no gene-gene interactions (epistasis); (ii) there are no gene-environment interactions; (iii) the effects of environmental factors shared by twins are independent of zygosity (the 'equal environments' assumption). These assumptions have all been individually explored and shown to be important in influencing the conclusions drawn from twin and family data [16–18]. In addition, the magnitude of any gene-environment interaction is critically important in determining the utility of targeting environmental interventions by genotype [19]. Although a general methodology to analyze twin data without making these assumptions has been developed, the algebra becomes intractable once multiple loci are involved [17]. This is problematic because, for common diseases, the impacts of multiple genetic variants, and potentially the whole genetic sequence, on disease susceptibility (here called 'genotypic risk') may be important.

The four-category model of population risks developed by Khoury and others [19] is a useful starting point for a top-down analysis of genetic prediction and prevention. It allows the merits of a targeted intervention strategy (which seeks to reduce the exposure of the high-risk genotype group only) to be explored, and can readily be extended to include more than four risk categories [10]. However, this model's use to date has been limited to bottom-up consideration of single genetic variants or to studying hypothetical examples of multiple variants. The four-category model is limited by the assumption of no confounders, which means it is applicable to only a subset of possible models of gene-gene and gene-environment interaction. However, situations where the 'no confounders' assumption is valid are arguably most likely to be of relevance to public health.

The aim of this paper is to combine the four-category model with population level data from twin, family and environmental studies, without adopting the classical twin model assumptions. This model of gene-gene and gene-environment interactions is then used to implement a 'top-down' approach to quantifying the utility of genetic 'prediction and prevention'.

## Method

### The four-category model

Consider a population divided into genotypic or environmental risk categories for a given trait (Figure 1a and 1b). The fraction of the population in the 'high environmental risk group' (designated by subscript e) is ε, and this subpopulation is at risk r_{e}. The remainder of the population is at risk r_{oe}. The fraction of the population in the 'high genotypic risk' group (designated by the subscript g) is γ, and this subpopulation is at risk r_{g}, with the remainder of the population at risk r_{og}. The total risk r_{t} for this trait in this population is then given by:

*r*_{
t
}= γ*r*_{
g
}+ (1-γ)*r*_{
og
} (1)

or by:

*r*_{
t
}= ε*r*_{e} + (1-ε)*r*_{
oe
} (2)

The same population can alternatively be divided into four categories, making a four-category model (Figure 1c)) with risks R_{oo}, R_{oe}, R_{go} and R_{ge}. Table 1 shows the risk categories in this model.

The risks are related to the previous definitions by:

*r*_{
g
}= *ε* *R*_{
ge
}+ (1-*ε*) *R*_{
go
} (3)

*r*_{
og
}= *ε* *R*_{
oe
}+ (1-*ε*) *R*_{
oo
} (4)

*r*_{e} = *γ* *R*_{
ge
}+ (1-*γ*) *R*_{
oe
} (5)

*r*_{
oe
}= *γ* *R*_{
og
}+ (1-*γ*) *R*_{
oo
} (6)

The category risks R remain constant in different populations (i.e. as ε and γ vary), provided there are no confounders. This assumption restricts the model to special cases of gene-gene and gene-environment interaction. Note that for a single genetic variant, r_{g} corresponds to the penetrance of the variant, and that in general (provided R_{ge} ≠ R_{go}) this varies with the proportion of the population in the high exposure group, ε, as has been observed [20, 21].

The total risk for the given trait is given by:

*r*_{
t
}= *γ* *ε* *R*_{
ge
}+ *γ*(1-*ε*)*R*_{
go
}+ *ε*(1-*γ*)*R*_{
oe
}+ (1-*ε*)(1-*γ*)*R*_{
oo
} (7)

The subpopulation of cases has different characteristics from the general population: for example, it contains a higher proportion of people from the 'ge' subgroup. The relative risk for a person drawn randomly from a subpopulation with the same genotypic and environmental characteristics as the cases, RR^{cases}, is given by the sum of the relative risks for each category shown in Table 1:

R{R}^{cases}=\frac{\gamma \epsilon {R}_{ge}^{2}+\gamma (1-\epsilon ){R}_{go}^{2}+\epsilon (1-\gamma ){R}_{oe}^{2}+(1-\epsilon )(1-\gamma ){R}_{oo}^{2}}{{r}_{t}^{2}}\left(8\right)

Similarly, the relative risk for a person drawn randomly from a subpopulation with the same genotypic characteristics as the cases (but with the environmental characteristics of the general population) is:

R{R}_{gen}^{cases}=\frac{\gamma {r}_{g}^{2}+(1-\gamma ){r}_{og}^{2}}{{r}_{t}^{2}}\left(9\right)

The relative risk for a person drawn randomly from a subpopulation with the same environmental characteristics as the cases (but with the genotypic characteristics of the general population) is:

R{R}_{env}^{cases}=\frac{\epsilon {r}_{e}^{2}+(1-\epsilon ){r}_{oe}^{2}}{{r}_{t}^{2}}\left(10\right)

### Population attributable fractions

Provided there are no confounders, the population attributable fraction (PAF^{E}_{e}) due to the presence of the high exposure (E) in the high exposure population subgroup (e) may be defined as:

PA{F}_{e}^{E}=\frac{\epsilon ({r}_{e}-{r}_{oe})}{{r}_{t}}=\epsilon \left\{\gamma ({R}_{ge}-{R}_{go})+(1-\gamma )({R}_{oe}-{R}_{oo})\right\}/{r}_{t}\left(11\right)

If the trait is a disease, PAF^{E}_{e} is the proportion of cases that could be avoided if an environmental intervention (such as a lifestyle change or reduction in exposure) succeeds in moving everyone in the 'high environmental risk group' to the 'low environmental risk' category, as shown in Figure 1b.

The targeted population attributable fraction (PAF^{E}_{ge}) may be defined as the proportion of cases that could be avoided by targeting the same environmental intervention at the 'high genotypic + high environmental risk' subgroup only (the 'ge' subgroup), as shown in Figure 1c. Again assuming no confounders, it is given by:

PA{F}_{ge}^{E}=\epsilon \gamma ({R}_{ge}-{R}_{go})/{r}_{t}\left(12\right)

Note that PAF^{E}_{ge} differs from PAF_{ge} as defined by Khoury & Wagener [19]. The latter implicitly assumes that both environmental and genetic risk factors are reduced and thus is inappropriate for assessing the merits of a targeted environmental intervention. PAF^{E}_{ge} as defined here is instead equivalent to the targeted attributable fraction (AF_{T}) defined by Khoury et al. [10]. To avoid confusion, the notation adopted here specifies both the nature of the intervention (environmental, denoted by superscript E) and the target subpopulation (the 'ge' subgroup, at both high genotypic and high environmental risk). Thus, the proportion of cases that would be avoided were it possible to move the 'high genotypic risk' subgroup to 'low genotypic risk' (as shown in Figure 1a) is written as PAF^{G}_{g}, given by:

PA{F}_{g}^{G}=\frac{\gamma ({r}_{g}-{r}_{og})}{{r}_{t}}=\gamma \left\{\epsilon ({R}_{ge}-{R}_{oe})+(1-\epsilon )({R}_{go}-{R}_{oo})\right\}/{r}_{t}\left(13\right)

Although in practice it is not possible to change the genotype of the population, the parameter PAF^{G}_{g} is nevertheless useful in the calculations that follow.

### Measures of utility

Khoury et al. [10] define the Population Impact (PI) as:

PI=PA{F}_{ge}^{E}/PA{F}_{e}^{E}\left(14\right)

PI is one possible measure of the usefulness of targeting the environmental intervention (E) at the 'ge' subgroup. It measures the proportion of cases avoided by targeting the 'high genotypic + high environmental risk' subgroup (the 'ge' subgroup), compared to the proportion avoided by applying the environmental intervention to the whole 'high environmental risk' group. PI has the property:

0 ≤ *PI* ≤ 1 (15)

and has its maximum value when PAF^{E}_{ge} = PAF^{E}_{e}. However, as a measure of the utility of genotyping, PI has the disadvantage that it takes no account of the proportion of the population γ in the high genotypic risk group. This means PI = 1 when γ = 1 simply because the whole population is then in the high genotypic risk group, although using genotyping to target environmental interventions is more likely to be useful if PI = 1 and γ is also small.

Therefore, consider an alternative utility parameter U_{ge}, defined by:

{U}_{ge}=\frac{PA{F}_{ge}^{E}}{PA{F}_{e}^{E}}-\gamma =\frac{\gamma (1-\gamma )\left[({R}_{ge}-{R}_{go})-({R}_{oe}-{R}_{oo})\right]}{\left[\gamma ({R}_{ge}-{R}_{go})+(1-\gamma )({R}_{oe}-{R}_{oo})\right]}\left(16\right)

which has the property

*γ* ≤ *U*_{
ge
}≤ (1-*γ*) (17)

U_{ge} tends to 1 only if PI = 1 and γ is also small. It is a measure of the utility of using genotyping to target the environmental intervention at the 'ge' subgroup, compared to randomly selecting the same proportion γ of the population to receive the intervention. U_{ge} is positive if those at high genotypic risk have *more to gain* than those at low genotypic risk from the intervention ((R_{ge}-R_{go}) ≥ (R_{oe}-R_{oo})) and negative if they have *less to gain* from the intervention. This reflects the fact that targeting those who have least to gain through an intervention is worse than using random selection in terms of its impact on population health.

Note that even if genotyping is better than random selection, other types of test that are more useful may be available [22]; a population-based approach still has the potential to reduce more cases of disease [9, 19, 23]; and such targeting also has broader psychological and social implications. Therefore a positive U_{ge} does not necessarily imply that genotyping is the best means of selecting a subpopulation to target, or that a targeted approach is necessarily effective or socially acceptable. Note also that the measure U_{ge} applies only to interventions that are considered applicable to the whole population (such as smoking cessation) and neglects other relevant issues such as cost-effectiveness and the burden of disease [24]. In addition, it is necessary to consider the magnitude of the Population Attributable Fraction, PAF^{E}_{e} before proposing this approach. This is because both PI and U_{ge} may tend to unity even if only a small proportion of cases can be avoided by means of environmental interventions.

### Limits on parameters

Consider only populations where r_{g} ≥ r_{og} and r_{e} ≥ r_{oe} for all values of ε and γ. Then the risks in the four box model must be ordered such that:

1 ≥ *R*_{
ge
}≥ *R*_{
oe
}≥ *R*_{
oo
}≥ 0 (18)

and

*R*_{
ge
}≥ *R*_{
go
}≥ *R*_{
oo
} (19)

Using the known relationships (Equations (11), (13) and (16)) between PAF^{E}_{e}, PAF^{G}_{g}, U_{ge} and the risks R_{oo}, R_{go}, R_{oe} and R_{ge}, leads to the limits on the utility parameter U_{ge} shown in Table 2. These conditions also ensure that PAF^{E}_{e}, PAF^{G}_{g} and PAF^{E}_{ge} are all positive. The two remaining inequalities (R_{ge} ≤ 1 and R_{oo} ≥ 0) are considered later, where they are used to derive limits on the proportion of the population in the 'high genotypic risk' group, γ. This step is not possible at this stage because PAF^{E}_{e}, PAF^{G}_{g} and PAF^{E}_{ge} are themselves dependent on γ.

### The twin and familial risks model

Data from studies of monozygotic and dizygotic twins are commonly used to estimate the genetic and environmental variances V_{g} and V_{e} of a trait. Here, the aim is to use twin and other data to estimate the possible magnitudes of the population attributable fractions and measures of utility defined above. To do this it is necessary to estimate V_{g}, V_{e} and the variance due to gene-environment interaction, V_{ge}. The standard methodology for twin data analysis is inappropriate because it assumes V_{ge} = 0.

First note that we are interested in the extent to which relatives share *risk categories* (which may be either environmental or genotypic, or both), rather than a particular genetic variant. The probability that a relative of a proband is also a case depends on the extent to which their environmental and genotypic risks are correlated with those of the proband. Rather than adopting a specific form for the genetic model, define p^{rel}_{g} as the correlation in genotypic risk category (g) between relatives of type denoted by the superscript 'rel'. The parameter p^{rel}_{g} is the probability that the genotypic risk category (high or low) is identical by descent.

For monozygotic (MZ) twins, assumed to share their entire genome, p^{MZ}_{g} = 1. For dizygotic (DZ) twins and other siblings, who share half their genome, p^{DZ}_{g} = p^{sib}_{g} = 1/2 for a single allele model (dominant Mendelian disorder) or an additive polygenic model. For a two allele model (recessive Mendelian disorder) or the dominance term of a polygenic model (in which multiple pairs of alleles interact), p^{DZ}_{g} = p^{sib}_{g} = 1/4. Here, allowing for the possibility of multiple gene-gene interactions (epistasis), require only that:

1/2\ge {p}_{g}^{DZ}\ge 0\left(20\right)

The meaning of p^{DZ}_{g} and its relationship to the polygenic risk model first adopted by Ronald Fisher in 1918 is discussed further below.

Similarly, define p^{rel}_{e} as the correlation in environmental risk category (e) between relatives of type "rel", requiring only that:

1\ge {p}_{e}^{rel}\ge 0\left(21\right)

Assume that p^{rel}_{g} and p^{rel}_{e} are independent (so that there is no genotype-environment correlation) and that risks within a category are randomly distributed. The relative risk for a relative of type "rel" may then be written:

{\lambda}_{rel}=(1-{p}_{g}^{rel})(1-{p}_{e}^{rel})+{p}_{g}^{rel}(1-{p}_{e}^{rel})R{R}_{gen}^{cases}+(1-{p}_{g}^{rel}){p}_{e}^{rel}R{R}_{env}^{cases}+{p}_{g}^{rel}{p}_{e}^{rel}R{R}^{cases}\left(22\right)

Substituting for the relative risks RR^{cases}_{gen}, RR^{cases}_{env} and RR^{cases} using Equations (8), (9) and (10) leads (after some algebra) to:

{\lambda}_{rel}-1={p}_{g}^{rel}\frac{{V}_{g}}{{r}_{t}^{2}}+{p}_{e}^{rel}\frac{{V}_{e}}{{r}_{t}^{2}}+{p}_{g}^{rel}{p}_{e}^{rel}\frac{{V}_{ge}}{{r}_{t}^{2}}\left(23\right)

where

\frac{{V}_{e}}{{r}_{t}^{2}}=\frac{(1-\epsilon )}{\epsilon}{\left[PA{F}_{e}^{E}\right]}^{2}\left(24\right)

\frac{{V}_{g}}{{r}_{t}^{2}}=\frac{(1-\gamma )}{\gamma}{\left[PA{F}_{g}^{G}\right]}^{2}\left(25\right)

\frac{{V}_{ge}}{{r}_{t}^{2}}=\frac{(1-\epsilon )}{\epsilon \gamma (1-\gamma )}{\left[{U}_{ge}PA{F}_{e}^{E}\right]}^{2}\left(26\right)

Note that if the G-E interaction component of the variance, V_{ge}, is zero, the utility of targeting the environmental intervention by genotype, U_{ge}, is also zero (Equation (26)), because those at high genotypic risk have no more to gain from the intervention than those at low genotypic risk (R_{ge}-R_{go} = R_{oe}-R_{oo}).

Equation (23) can also be derived more formally using matrix methods (Appendix A).

### The gene-environment interaction factor and remaining inequalities

Without loss of generality, define the gene-environment interaction factor f_{ge} such that:

\frac{{V}_{ge}}{{r}_{t}^{2}}={f}_{ge}^{2}\frac{{V}_{g}}{{r}_{t}^{2}}.\frac{{V}_{e}}{{r}_{t}^{2}}\left(27\right)

and choose its sign so that (combining Equations (24), (25) and (26)):

{U}_{ge}={f}_{ge}\sqrt{\gamma (1-\gamma )\frac{{V}_{g}}{{r}_{t}^{2}}}\left(28\right)

U_{ge} is zero if f_{ge} = 0 (i.e. for an additive G-E model, with no G-E interaction), but for a given γ and V_{g}, U_{ge} increases with increasing gene-environment interaction factor, f_{ge}. For a fixed f_{ge} and genetic variance component V_{g}, U_{ge} is maximum when γ = 1/2, i.e. when half the population is in the high genotypic risk group, provided solutions with γ = 1/2 exist (see also below: *cases where γ*_{maxge} *< 1/2*).

Using the definitions of V_{e}, V_{g} and V_{ge} (Equations (24), (25) and (26)) and the remaining inequalities, R_{ge} ≤ 1 and R_{oo} ≥ 0, two limits can be derived on the proportion of the population in the 'high genotypic risk' group, γ (see Table 2).

### Scoping studies

The general system of equations represented by Equation (23) may be simplified where data exist from monozygotic twins, dizygotic twins and other siblings, such that λ_{DZ} > λ_{sib}. This implies that environmental risks are more strongly correlated in dizygotic twins than in other siblings, p^{e}_{DZ} > p^{e}_{sib}. Remembering that p^{MZ}_{g} = 1 and p^{sib}_{g} = p^{DZ}_{g}, three independent equations for the relative risk in monozygotic, dizygotic twins and siblings may then be written:

{\lambda}_{MZ}-1=\frac{{V}_{g}}{{r}_{t}^{2}}+{p}_{e}^{MZ}\frac{{V}_{e}}{{r}_{t}^{2}}+{p}_{e}^{MZ}\frac{{V}_{ge}}{{r}_{t}^{2}}\left(29\right)

{\lambda}_{DZ}-1={p}_{g}^{DZ}\frac{{V}_{g}}{{r}_{t}^{2}}+{p}_{e}^{DZ}\frac{{V}_{e}}{{r}_{t}^{2}}+{p}_{g}^{DZ}{p}_{e}^{DZ}\frac{{V}_{ge}}{{r}_{t}^{2}}\left(30\right)

{\lambda}_{sib}-1={p}_{g}^{DZ}\frac{{V}_{g}}{{r}_{t}^{2}}+{p}_{e}^{sib}\frac{{V}_{e}}{{r}_{t}^{2}}+{p}_{g}^{DZ}{p}_{e}^{sib}\frac{{V}_{ge}}{{r}_{t}^{2}}\left(31\right)

To solve, assume the recurrence risks λ are known (see Appendix B and [25]) and define:

{R}_{MD}=\frac{{\lambda}_{MZ}-1}{{\lambda}_{DZ}-1}\left(32\right)

{R}_{SD}=\frac{{\lambda}_{sib}-1}{{\lambda}_{DZ}-1}\left(33\right)

with

*R*_{
MD
}≥ 1 (34)

and

0 ≤ *R*_{
SD
}≤ 1. (35)

Note that if R_{SD} = 1, Equations (30) and (31) are identical, p^{e}_{DZ} = p^{e}_{sib}, and more relatives are needed to obtain solutions, except in the special case where there is no environmental variance (see below: *no environmental variance*).

In addition, define the variable parameters (assumed unknown):

{c}_{MD}=\frac{{p}_{e}^{MZ}}{{p}_{e}^{DZ}}\left(36\right)

{c}_{SD}=\frac{{p}_{e}^{sib}}{{p}_{e}^{DZ}}\left(37\right)

with

*c*_{
MD
}≥ 1 (38)

and

0 ≤ *c*_{
SD
}≤ 1. (39)

For λ_{DZ} > 1 and R_{SD} < 1 the simultaneous Equations (29), (30) and (31) can then be solved to give:

\frac{{V}_{g}}{{r}_{t}^{2}}=\frac{\left({\lambda}_{DZ}-1\right)}{{p}_{g}^{DZ}}.\frac{\left({R}_{SD}-{c}_{SD}\right)}{\left(1-{c}_{SD}\right)}\left(40\right)

\frac{{V}_{e}}{{r}_{t}^{2}}=\frac{({\lambda}_{DZ}-1)}{{p}_{e}^{DZ}{c}_{MD}(1-{p}_{g}^{DZ})}\left[\frac{({c}_{MD}-1)(1-{R}_{SD})}{(1-{c}_{SD})}+(1-{p}_{g}^{DZ}{R}_{MD})\right]\left(41\right)

\frac{{V}_{ge}}{{r}_{t}^{2}}=\frac{({\lambda}_{DZ}-1)}{{p}_{e}^{DZ}{p}_{g}^{DZ}{c}_{MD}(1-{p}_{g}^{DZ})}\left[\frac{(1-{c}_{MD}{p}_{g}^{DZ})(1-{R}_{SD})}{(1-{c}_{SD})}+(1-{p}_{g}^{DZ}{R}_{MD})\right]\left(42\right)

provided {p}_{g}^{DZ} ≠ 0, {p}_{e}^{DZ} ≠ 0 and *c*_{
SD
}≠ 1 (see also below).

For situations in which a targeted intervention is under consideration, the population attributable fraction PAF^{E}_{e} and exposure ε are likely to be known, allowing V_{e} to be treated as an input variable. However, p^{DZ}_{e} is usually unknown, since environmental correlations are often difficult to measure. Therefore, it is useful to eliminate p^{DZ}_{e} from Equations (41) and (42), leading to:

\frac{{V}_{ge}}{{V}_{e}}=\frac{\left\{\frac{{p}_{g}^{DZ}}{{p}_{g\mathrm{min}}^{DZ}}-1\right\}\frac{({R}_{SD}-{c}_{SD})}{(1-{c}_{SD})}}{{p}_{g}^{DZ}{R}_{MD}({p}_{gtop}^{DZ}-{p}_{g\mathrm{min}}^{DZ})}\left(43\right)

where

{p}_{gtop}^{DZ}=\frac{1}{{R}_{MD}}\left\{1+\frac{({c}_{MD}-1)(1-{R}_{SD})}{(1-{c}_{SD})}\right\}\left(44\right)

and

{p}_{g\mathrm{min}}^{DZ}=\frac{({R}_{SD}-{c}_{SD})}{\left\{{R}_{MD}(1-{c}_{SD})-{c}_{MD}(1-{R}_{SD})\right\}}\left(45\right).

Equations (27), (40) and (43) allow the gene-environment interaction factor f_{ge} to be written as:

{f}_{ge}^{2}=\frac{\left\{\frac{{p}_{g}^{DZ}}{{p}_{g\mathrm{min}}^{DZ}}-1\right\}}{({\lambda}_{DZ}-1){R}_{MD}({p}_{gtop}^{DZ}-{p}_{g}^{DZ})}\left(46\right).

The parameter p^{DZ}_{g}, which defines the form of the genetic model, is then given by:

\frac{{p}_{g}^{DZ}}{{p}_{g\mathrm{min}}^{DZ}}=\frac{1+{f}_{ge}^{2}({\lambda}_{DZ}-1){R}_{MD}{p}_{gtop}^{DZ}}{1+{f}_{ge}^{2}({\lambda}_{DZ}-1){R}_{MD}{p}_{g\mathrm{min}}^{DZ}}\left(47\right).

For known R_{MD}, R_{SD} and λ_{DZ} a solution space can now be mapped, which includes all possible variances consistent with the data and with the inequalities derived above.

Requiring the variances to be positive leads to the additional conditions on p^{DZ}_{g} and c_{SD} shown in Table 3.

The limits on U_{ge} shown in Table 2 set limits on the range of gene-environment interaction models such that:

-\frac{\epsilon}{(1-\epsilon )PA{F}_{e}^{E}}\le {f}_{ge}\le \frac{1}{PA{F}_{e}^{E}}\left(48\right)

Noting that f_{ge} = 0 corresponds to p^{DZ}_{g} = p^{DZ}_{gmin} (Equation (64)), this implies that, for U_{ge} ≥ 0, the solution space may be defined by:

{p}_{g\mathrm{min}}^{DZ}\le {p}_{g}^{DZ}\le {p}_{g\mathrm{max}}^{DZ}\left(49\right)

where p^{DZ}_{gmax} is given by Equation (47) with f_{ge} = 1/PAF^{E}_{e}.

For U_{ge} ≤ 0, the solution space may be defined by:

{p}_{g\mathrm{min}}^{DZ}\le {p}_{g}^{DZ}\le {p}_{gneg}^{DZ}\left(50\right)

where p^{DZ}_{gneg} is given by Equation (47) with f_{ge} = -ε/(1-ε)PAF^{E}_{e}.

The remaining limits on U_{ge} lead to the additional conditions on the range of γ values (the proportion of the population in the high risk group) shown in Table 2. These conditions on γ may be written:

*γ*_{min} ≤ *γ* ≤ *γ*_{max} (51)

where (noting that γ_{maxge} = γ_{o} when f_{ge} = 1):

{\gamma}_{\mathrm{max}}=\{\begin{array}{c}{\gamma}_{\mathrm{max}ge}\text{for}{f}_{ge}\ge 1\\ {\gamma}_{0}\text{for}{f}_{ge}\le 1\end{array}\left(52\right)

and (noting that γ_{minge} = γ_{neg} when f_{ge} = -r_{t}/(1-r_{t})):

{\gamma}_{\mathrm{min}}=\{\begin{array}{c}{\gamma}_{\mathrm{min}ge}\text{for}{f}_{ge}\ge -{r}_{t}/\left(1-{r}_{t}\right)\\ {\gamma}_{neg}\text{for}{f}_{ge}\le -{r}_{t}/\left(1-{r}_{t}\right)\end{array}\left(53\right)

Two transition lines can therefore be defined such that p^{DZ}_{g} = p^{DZ}_{gt} when f_{ge} = 1 and p^{DZ}_{g} = p^{DZ}_{gnegt} when f_{ge} = -r_{t}/(1-r_{t}). The values of p^{DZ}_{gt} and p^{DZ}_{gnegt} may be calculated using Equation (47).

The full range of gene-environment interaction models specified by f_{ge} (within the limits given by Equation (48)) and the corresponding range of γ values are summarized in Table 4. Note that the risk distribution associated with f_{ge} = 1 corresponds to a multiplicative model of gene-environment interaction. If f_{ge} ≥ 1 solutions with population impact PI = 1 may exist (i.e. with PAF^{E}_{ge} = PAF^{E}_{e}), provided the proportion of the population in the high risk genotypic group takes the maximum value consistent with the data (γ = γ_{maxge}). For lower values of f_{ge}, solutions with PI = 1 cannot exist.

One additional condition is necessary for solutions to exist, namely:

*γ*_{max} ≥ *γ*_{min} (54)

This condition is always met if

*λ*_{
MD
}≤ *y*_{
e
}+ 1 (55)

where

{y}_{e}=\{\begin{array}{c}{F}_{1}/{f}_{ge}\text{for}{f}_{ge}\ge 1\\ {F}_{1}/{F}_{2}\text{for}1\ge {f}_{ge}\ge -{r}_{t}/\left(1-{r}_{t}\right)\\ -{F}_{2}/{f}_{ge}\text{for}{f}_{ge}\le -{r}_{t}/\left(1-{r}_{t}\right)\end{array}\left(56\right)

and F_{1} and F_{2} are given by:

{F}_{1}=\frac{\left[\left(\frac{1-{r}_{t}}{{r}_{t}}\right)-\left(\frac{1-\epsilon}{\epsilon}\right)PA{F}_{e}^{E}\right]}{\left[1+{f}_{ge}\left(\frac{1-\epsilon}{\epsilon}\right)PA{F}_{e}^{E}\right]}={\scriptscriptstyle \frac{(1-{r}_{e})}{\left[{r}_{t}+{f}_{ge}({r}_{e}-{r}_{t})\right]}}\left(57\right)

{F}_{2}=\frac{\left(1-PA{F}_{e}^{E}\right)}{\left(1-{f}_{ge}PA{F}_{e}^{E}\right)}\left(58\right).

However, if λ_{MD} is greater than this, the requirement γ_{max} ≥ γ_{min} further restricts the values of c_{SD} that lie within the solution space (Table 3).

If V_{e} and ε are known, a solution space can be now be mapped for p^{DZ}_{g} and f_{ge} with known input data from twin and sibling studies (λ_{MZ}, λ_{DZ} and λ_{sib}), for a given c_{MD} and all values of c_{SD} within the assumed range. The boundaries of the solution space are determined by the limits on f_{ge} given by Equation (48), the condition γ_{max} ≥ γ_{min} (Equation (54)), and the requirement that p^{DZ}_{g} is less than or equal to 1/2 (Equation (20)) – no other condition on the genetic model is specified a priori. For each genetic risk model and gene-environment interaction model in the solution space, defined by p^{DZ}_{g} and f_{ge} respectively, the variances V_{g} and V_{ge} can then be calculated, as can γ_{max} and γ_{min}. For a chosen γ value in the allowed range, U_{ge} can then be calculated from Equation (28).

The model code is available as [Additional file 1: heritability12.xls].

Note that the condition on p^{DZ}_{g} ≤ 1/2 may also be rewritten using Equation (47), so that:

{p}_{g}^{DZ}\le 1/2\Rightarrow \frac{\left({p}_{g\mathrm{min}}^{DZ}-\frac{1}{2}\right)}{{p}_{g\mathrm{min}}^{DZ}}\le {R}_{MD}{f}_{e}^{2}\left({\lambda}_{DZ}-1\right)\left(1/2-{p}_{gtop}^{DZ}\right)\left(59\right)

which is always met if

{p}_{gtop}^{DZ}\le 1/2\left(60\right).

Before mapping the solution space, first consider some special cases and a comparison of the model with the classical twin studies approach.

### Special cases

#### 1. No genetic variance

If V_{g} = 0, Equation (27) implies that V_{ge} = 0 also. Equations (29), (30) and (31) then give:

*R*_{
SD
}= *c*_{
SD
} (61)

and

*R*_{
MD
}= *c*_{
MD
} (62)

Under the usual assumption that c_{MD} = 1 (the 'equal environments' assumption), this is the well-known result that genetic variance can be zero only when the concordance in monozygotic and dizygotic twins is the same (leading to R_{MD} = 1). However, if the equal environments assumption is not met (c_{MD} > 1), values of R_{MD} greater than 1 do not necessarily imply that a genetic component to the variance exists (see, for example, [18]).

#### 2. No environmental variance

If V_{e} = 0, Equation (27) implies that V_{ge} = 0 also. Equations (29), (30) and (31) then give:

*R*_{
SD
}= 1 (63)

and

{R}_{MD}=1/{p}_{g}^{DZ}\left(64\right)

For a purely genetic model with no environmental variance, Equation (64) implies that if R_{MD} > 2, p^{DZ}_{g} < 1/2. This is consistent with Risch's finding [16] that neither an additive genetic model nor a single dominant gene model (both with p^{DZ}_{g} = 1/2) can fit the data for conditions such as schizophrenia (which has an R_{MD} value significantly greater than 2).

#### 3. Classical twin study assumptions

Assuming no gene-environment interaction (V_{ge} = 0); an additive genetic risk model (p^{DZ}_{g} = 1/2); and the 'equal environments' assumption (c_{MD} = 1) in Equations (29), (30) and (31) gives:

\frac{{V}_{g}}{{r}_{t}^{2}}=2\left({\lambda}_{MZ}-{\lambda}_{DZ}\right)\left(65\right)

This is the classical twin study result, assuming the dominance term of the genetic variance is negligible. Note that, if R_{MD} = 2, the classical solution implies that the environmental variance terms in Equations (29) to (31) are zero and shared sibling risk is due to entirely to shared genes.

#### 4. No correlation in genotypic risk in siblings (p^{DZ}_{g} = 0)

Equation (20) allows p^{DZ}_{g} to tend to zero. Substituting p^{DZ}_{g} = 0 in Equations (29), (30) and (31) and using the definition of the gene-environment interaction factor (Equation (28)) gives:

*R*_{
SD
}= *c*_{
SD
} (66)

and

\frac{{V}_{g}}{{r}_{t}^{2}}=\frac{\left({\lambda}_{DZ}-1\right)\left({R}_{MD}-{c}_{MD}\right)}{\left[1+{f}_{ge}^{2}{c}_{MD}\left({\lambda}_{DZ}-1\right)\right]}\left(67\right)

Note that, from Equations (30) and (31), p^{DZ}_{g} = 0 corresponds to a purely environmental explanation for shared sibling risks (although there may remain a genetic component to shared risks in monozygotic twins, from Equation (29)). The solution p^{DZ}_{g} = 0 may not exist in reality; however, the solution at this limit is of interest because low values of p^{DZ}_{g} are plausible.

Also, note that if f_{ge} = 0 (no gene-environment interaction) and c_{MD} = 1 (the 'equal environments' assumption), the genetic variance V_{g} given by Equation (67) is half the classical twin study result (Equation (65)).

#### 5. Cases where γ_{max} = γ_{min}

If the line γ_{max} = γ_{min} exists within the solution space, some special cases may arise with risk distributions of particular interest (including, for example, a solution with R_{ge} = 1 and all other risks zero). These special cases and the conditions that they meet are shown in Table 5.

#### 6. Cases where γ_{maxge} < 1/2

Equation (27) shows that for a fixed gene-environment interaction factor f_{ge} and genetic variance component V_{g}, the utility U_{ge} is maximum when γ = 1/2, i.e. when half the population is in the high genotypic risk group, provided this solution exists. However, if γ_{max} < 1/2, utility is maximum when γ = γ_{max}. As a smaller proportion of the population is then targeted, these solutions are of particular interest. Because solutions with population impact PI = 1 may exist when 1 ≤ f_{ge} ≤ 1/PAF^{E}_{e} if γ = γ_{maxge} (Table 4), it is of interest to identify the area of the solution space with γ_{maxge} < 1/2. Maximum utility is then obtained when γ = γ_{maxge} (where PI = 1 and U_{ge} = 1-γ_{maxge}). For the condition

{\gamma}_{\mathrm{max}ge}<1/2\Rightarrow {p}_{g}^{DZ}>{p}_{gx}^{DZ}\left(68\right)

where p^{DZ}_{gx} is given by:

{R}_{MD}(1-{c}_{SD}){({p}_{gx}^{DZ})}^{2}+\left[(1-{c}_{SD})({R}_{MD}-1)-(2{c}_{MD}-1)(1-{R}_{SD})\right]{p}_{gx}^{DZ}-({R}_{SD}-{c}_{SD})=0\left(69\right)

solving for p^{DZ}_{gx} allows the region of the solution space where γ_{maxge} < 1/2 to be defined.

#### 7. Cases where the 'equal environments' assumption holds (c_{MD} = 1)

In the special case where the 'equal environments' assumption holds (c_{MD} = 1, and hence p^{DZ}_{gtop} = 1/R_{MD}), Equation (63) simplifies to give R_{MD} ≥ 2. Equation (62) also simplifies to give:

{p}_{g}^{DZ}\le 1/2\Rightarrow {c}_{SD}\ge {c}_{1}\left(70\right)

where

{c}_{1}=1-\frac{(1-{R}_{SD})\left[1+{f}_{ge}^{2}({\lambda}_{DZ}-1)(2-{R}_{MD})\right]}{(2-{R}_{MD})\left[1+{f}_{ge}^{2}({\lambda}_{DZ}-1)\right]}\left(71\right)

Meeting the condition p^{DZ}_{g} ≤ 1/2 at c_{SD} = 0 then requires:

{R}_{MD}\ge 2-\frac{(1-{R}_{SD})}{\left[1+{f}_{ge}^{2}({\lambda}_{DZ}-1){R}_{SD}\right]}\left(72\right).

It follows that if c_{MD} = 1, solutions with p^{DZ}_{g} = 1/2 (an additive genetic model) and positive utility exist only when the following condition holds for R_{MD}:

{R}_{MD}\le 2-\frac{(1-{R}_{SD})}{\left[1+({\lambda}_{DZ}-1){R}_{SD}/{(PA{F}_{e}^{E})}^{2}\right]}\left(73\right).

Further, all three classical twin study assumptions (c_{MD} = 1, p^{DZ}_{g} = 1/2 and f_{ge} = 0) can be met only for values of R_{MD} that are low enough to satisfy:

1 + *R*_{
SD
}≥ *R*_{
MD
}> 1 (74).

If R_{MD} lies within this range, the classical twin study gives one possible solution; however, other solutions also exist. All alternative solutions favour a less 'genetic' and more 'environmental' explanation for shared sibling risks (i.e. they have higher values of c_{SD}). If R_{MD} is greater than 1+R_{SD}, all three assumptions of the classical twin study cannot be met simultaneously.

### Comparison with the classical twins approach

Table 6 summarizes the differences between the classical twin studies approach and the method adopted here.

A central feature of the model is that it abandons Fisher's assumption [26] that genes act as risk factors for common traits in a manner necessarily dominated by an additive polygenic term. In his historic 1918 paper, Fisher synthesized Mendelian inheritance with Darwin's theory of evolution by showing that the genetic variance of a continuous trait could be decomposed into additive and non-additive components [26, 27]. Following Fisher, the classical twin study analysis depends on writing the genetic component of a trait as a convergent series of terms, consisting of an additive term (the sum of contributions of individual alleles at each locus) plus a smaller dominance term (the sum of contributions from pairs of alleles at each locus) and – usually neglected – epistatic terms (involving potentially multiple interactions between alleles at multiple loci) [15]. Often the additive term is assumed to dominate the series (equivalent to assuming p^{DZ}_{g} = 1/2).

Fisher saw his polygenic model as "*abandon* [ing] *the strictly Mendelian mode of inheritance, and treat* [ing] *Galton's 'particulate inheritance' in almost its full generality*" [26]. However, it can be argued that Fisher's model is flawed in so far as it fails to distinguish between the function of alleles and the properties of traits [4, 28]. In particular, epistasis (although referred to here as 'gene-gene interaction') is not strictly an interaction between genes, but can be shown to depend on the structure and interdependence of metabolic pathways [28].

The alternative model adopted here is based on correlations in *risk categories* for a trait (which may be either environmental or genetic, or both), rather than single or multiple genetic variants. Adopting Porteous' critique [28], there is no *a priori* biological reason why the parameter p^{DZ}_{g} (the probability that the genotypic risk category of a dizygotic twin pair is identical by descent) cannot take any value between 1/2 (its value if the additive model holds) and zero. Low p^{DZ}_{g} can then be understood to mean either a situation in which Fisher's polygenic model [26] is dominated by negative (synergistic) epistatic terms (for example, p^{DZ}_{g} = 1/2^{n} implies that interactions between n deleterious alleles are necessary to produce a phenotypic effect), or, more meaningfully, a situation in which human phenotypes are *biologically robust* to individual genetic variants [29]. Thus, in the extreme case where numerous genetic variants combine to influence a trait through the interdependence of metabolic pathways, the trait may be highly correlated in monozygotic twins (who share all the genetic variants) but not correlated at all (p^{DZ}_{g} = 0) in dizygotic twins or siblings (who share only half the relevant variants by descent). Although p^{DZ}_{g} = 0 may not be realistic, low values of p^{DZ}_{g} are plausible, and may even be typical of complex diseases.

The classical twin study assumptions (see above) allow a single solution to be calculated from the under-determined system of simultaneous Equations (29), (30) and (31). However, in the absence of prior knowledge about the form of the genetic model, the presence or absence of gene-environment interactions, and the validity of the 'equal environments' assumption, the approach adopted here is more rigorous.

## Results

### General model solutions

First consider the behaviour of the model when the 'equal environments' assumption holds and hence c_{MD} = 1 (as described above).

Figures 2, 3 and 4 show the possible solution spaces for an arbitrary set of plausible input parameters satisfying the requirement R_{MD} > 1+R_{SD} necessary for the classical twin study solution to exist. In Figure 2 the gene-environment interaction factor f_{ge} and hence utility, U_{ge}, are both positive and in Figure 3 they are negative. The horizontal axis shows c_{SD}/R_{SD}, which is zero if shared sibling risk is due to shared genetic factors only and 1 if shared sibling risk is due to shared environmental factors only. The vertical axis shows p^{DZ}_{g}, which is 1/2 if the additive genetic model holds, but may reduce to zero if epistasis dominates and the phenotype is robust to genetic variation. The three curved solid lines represent three models of gene-environment (G-E) interaction: an additive G-E model (i.e. no gene-environment interaction, f_{ge} = 0); a multiplicative G-E model (f_{ge} = 1); and maximum G-E interaction (f_{ge} = 1/PAF^{E}_{e}). The possible solution spaces are shaded grey. Each point in each shaded solution space corresponds to a given genetic model (defined by p^{DZ}_{g}) and a given G-E interaction model (defined by f_{ge}). Figure 4 plots the entire solution space (including both negative and positive utility) by transforming the horizontal axis to represent the G-E interaction parameter, f_{ge}. Although the classical twin model can fit the data, an infinite number of other solutions corresponding to different genetic and gene-environment interaction models also exist. In this example, the line γ_{max} = γ_{min} lies outside the solution space and no solutions exist with γ_{maxge} < 1/2.

For lower values of R_{MD}, the curves defining the solution space are shifted downwards [see Additional files 2 to 9], so that the line f_{ge} = 0 (corresponding to no gene-environment interaction) lies entirely below the line p^{DZ}_{g} = 1/2 (corresponding to an additive genetic model). The classical twin study solution does not exist, but many other combinations of genetic and gene-environment interaction models may fit the data.

When c_{MD} > 1, lines of constant f_{ge} no longer decrease monotonically to zero, and are also shifted upwards, so that solutions with strong G-E interactions are no longer possible [see Additional files 10 to 12].

### Example applications using twin, sibling and environmental data

#### Input values

Consider example applications of the model for male lung cancer, female breast cancer and schizophrenia. The model input variables used are shown in Table 7.

The recurrence risks, λ, and total risks, r_{t}, for breast and lung cancer are those calculated by Risch [30], based on Scandinavian twin data reported by Lichtenstein et al. [31] (involving more than 44,000 twin pairs) and Swedish familial data reported by Doug and Hemminki [32] (involving more than 2 million families). The proportion of the population exposed, ε, and population attributable fraction, PAF^{E}_{e}, for breast cancer are taken from those reported by Rockhill et al. [33] for a US population. Although strictly speaking these values may not be appropriate for a Scandinavian population, and include a component due to family history that may be (at least partly) genetic, they give a low V_{e}, consistent with the known environmental risk factors for breast cancer, and results are not sensitive to these input values (because V_{e} is so small). For lung cancer, it is assumed that 15% of the Scandinavian population smokes and that 86% of lung cancer cases could be avoided if they did not (giving a risk of lung cancer in smokers of 10%).

The recurrence risks λ, and total risk, r_{t}, for schizophrenia are those used by Risch [16], based on European data summarized by McGue et al. [34]. More recent twin studies for schizophrenia have given variable results and this example should be treated as illustrative only. Further, environmental exposures and population attributable fractions are unknown for schizophrenia. Two exploratory sets of results are therefore reported, using data consistent with a low environmental variance (based on the values used for breast cancer), and high environmental variance (based on the values used for smoking and lung cancer).

Detailed results for the three diseases are shown in [Additional file 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33]. The key findings are outlined below.

#### Breast cancer results

For breast cancer, the PAF^{E}_{e} associated with known environmental factors is low. The value of the model is therefore less in calculating the utility of targeted environmental interventions than in exploring the solution space for a complex disease with R_{MD} close to 2.

Although strictly speaking the classical twin study solution (with an additive genetic model, p^{DZ}_{g} = 1/2, and an additive G-E model, f_{ge} = 0) does not exist as a solution, it might lie within the margin of error of the data. However, an infinite number of other models also could also fit the data. The classical twin model result always overestimates the genetic component of the variance, which reduces as the gene-environment interaction factor f_{ge} increases, and also as p^{DZ}_{g} decreases (i.e. as epistatic terms begin to dominate the genetic model). These alternative models imply that shared environmental factors may partially explain familial aggregation of breast cancer. This contrasts with the classical twin method result (see earlier), which for R_{MD} = 2 leads to the inevitable conclusion that shared sibling risk must be due solely to shared genes [35].

In theory, a model with p^{DZ}_{g} = 0 (where shared sibling risk is due entirely to shared environmental factors) could fit the data. However, for breast cancer the existence of known mutations that significantly increase risk (particularly mutations in the BRCA1 and BRCA2 genes, which are relatively common) rules out this solution. Although it is not possible to subtract out the effect of these mutations from the model, it is possible to show that they could be sufficient to explain the twin data if a G-E interaction also exists. For example, one possible solution consistent with the data could involve one or more dominant genes (p^{DZ}_{g} = 1/2), a strong G-E interaction (f_{ge} = 1/PAF^{E} _{e}), but a largely environmental explanation for shared sibling risk (say c_{SD}/R_{SD} = 0.9). This solution implies that the genetic component of the variance is less than a fifth of the classical twin study result, which could be low enough to be explained by mutations in the BRCA1 and BRCA2 genes alone [35]. If this model were correct it would have important implications for women with such mutations, but would not contribute significantly to reducing the incidence of breast cancer in the population as a whole, because the affected proportion of the population γ would be rather small. Other solutions, involving different genetic models with lower p^{DZ}_{g}, and/or less gene-environment interaction, are also possible.

The line γ_{max} = γ_{min} does not occur within the solution space for breast cancer; however, in some circumstances the lines γ_{max} and γ_{min} may be rather close together. This suggests that, although as expected there is always a trade-off between selecting a small proportion (γ_{min}) of the population with a high Positive Predictive Value (PPV), or a larger proportion of the population (γ_{max}) with a higher Population Impact (PI) [19], some possible solutions could exist for breast cancer where the PPV and PI are both relatively high. Further, γ_{max} is often less than 1/2, so that, in these regions of the possible solution space, maximum utility might be obtained by targeting less than 50% of the population. However, known environmental factors for breast cancer are often not amenable to intervention and other possible solutions, with low, zero or negative utility, also exist.

#### Lung cancer results

For lung cancer, all the possible solutions imply that shared sibling risk is largely due to shared environmental factors (smoking) because solutions occur only when c_{SD}/R_{SD} is close to 1. Unlike for breast cancer, the line γ_{max} = γ_{min} lies outside the solution space, even for negative f_{ge}, as does the area of solutions with γ_{maxge} < 1/2. However, the classical twin study solution, with f_{ge} = 0 and p^{DZ}_{g} = 1/2, clearly lies within the solution space.

Although the classical twin model again provides an upper limit to the genetic component of the variance, even the classical result indicates that the risk of lung cancer is dominated by smoking in this population and the variance has at most a small genetic component.

Unlike the breast cancer example, γ_{max} and γ_{min} are always far apart, suggesting a strong trade off between high Postive Predictive Value (R_{ge}) for a genotypic test and a high Population Impact (PI) for a targeted intervention. This means that a genotypic test that predicts which smokers will get lung cancer cannot exist. To predict all cases of lung cancer in smokers (i.e. to obtain PI = 1), 95% or more of the population would have to be in the high genotypic risk group, and the predictive value of such a test would be very low.

Because the genetic component of the variance is so small, it follows that the utility of genetic 'prediction and prevention' (measured by U_{ge}) is also small (from Equation (28)). Utility is maximum when γ = 1/2, but even then values are low. The maximum utility of genotyping occurs when about 60% of cases could be prevented by targeting the 50% of smokers at high genotypic risk. However, other possible solutions have zero or negative utility.

#### Schizophrenia results

For schizophrenia, the classical twin study solution (with f_{ge} = 0 and p^{DZ}_{g} = 1/2 and c_{MD} = 1) cannot not fit the data. If the 'equal environments' assumption holds, neither a single dominant gene (p^{DZ}_{g} = 1/2), nor additive polygenic model (also with p^{DZ}_{g} = 1/2), nor single recessive gene (p^{DZ}_{g} = 1/4) can explain the twin and family data, consistent with Risch's 1990 findings [16]. This may suggest that the genetic model for schizophrenia is likely to be dominated by epistatic terms. However, if gene-environment interactions are important, it is also possible that a recessive gene, combined with at least multiplicative G-E interaction (p^{DZ}_{g} = 1/4 and f_{ge} = 1 or higher), could explain the data.

The possible solution spaces include purely genetic explanations for shared sibling risk (at c_{SD}/R_{SD} = 0), or purely environmental ones (at c_{SD}/R_{SD} = 1, applicable if p^{DZ}_{g} = 0).

Assuming a small environmental component to the variance, there is no region of the solution space for which γ_{maxge} < 1/2, suggesting that the utility of targeted environmental interventions under these assumptions is likely to be low. However, if the environmental component of the variance is assumed to be much larger, the available solution space changes dramatically, because the line γ_{max} = γ_{min} now constrains the solution space to a much smaller area, which excludes solutions with no G-E interaction (f_{ge} = 0). Special solutions may exist along the line γ_{max} = γ_{min}, as shown in Table 5. Because the environmental factors contributing to schizophrenia are unknown, it is impossible to draw any conclusions about the potential benefits of targeting environmental interventions at those at high genotypic risk.

Because prenatal development is thought to be important in schizophrenia, it is plausible that monozygotic twins are more likely to share environmental risk factors than dizygotic twins are. Breaking the 'equal environments' assumption changes the shape of the solution space significantly, and, assuming a small environmental component to the variance, only limited G-E interactions are now possible (the multiplicative G-E model, f_{ge} = 1, lies largely outside the solution space). The utility of targeting environmental interventions by genotype is then likely to be low. However, in these circumstances it is possible that an additive genetic model (p^{DZ}_{g} = 1/2) with some G-E interaction, or a recessive gene (p^{DZ}_{g} = 1/4) with no G-E interaction, could explain the data.

## Discussion

If Fisher's polygenic model [26] is abandoned, along with the usual twin study assumption that there are no gene-environment interactions, the four-category model developed by Khoury and others can be combined with twin, family and environmental data to implement a 'top down' approach to assessing the utility of targeting environmental/lifestyle interventions by genotype. Scoping studies, valid when R_{SD} ≠ 1, provide a first step to modelling the health of populations [23].

Abandoning Fisher's assumption that the polygenic model is necessarily dominated by an additive term can be justified by the growing evidence that phenotypic effects can result from the synergistic action of alleles in many genes [36]. For example, Bardet-Biedl Syndrome, historically assumed to be a recessive trait, has been shown to involve three interacting mutations at two loci in some patients (implying that p^{DZ}_{g} = 1/8) and, more recently, an additional locus has been identified that can also interact to change disease severity and symptoms [37]. Both positive and negative gene-environment interactions have also been observed in human diseases, although there are difficulties in confirming their statistical validity [38, 39].

The model also allows the impact of the much criticised 'equal environments' assumption to be explored.

A number of conclusions can be drawn about the merits of the classical twin study and the utility of genetic 'prediction and prevention'.

Firstly, the model confirms that the classical twin study solution is not always valid and gives at best an upper limit to the genetic component of the variance of a trait. The importance of the 'equal environments' assumption and of gene-environment interactions have previously been recognised [17, 18]; however, less attention has been paid to the potential role of gene-gene interactions (epistasis). For larger values of R_{MD} (greater than 1+R_{SD}), observed for conditions such as schizophrenia, the model generalizes Risch's findings [16] to show that the three assumptions of the classical twin model cannot all be satisfied simultaneously. For intermediate R_{MD} values, observed for conditions such as breast cancer (for which R_{MD} is approximately 2), the model illustrates that the conclusion drawn from classical twin studies, that familial aggregation is due entirely to shared genetic factors, may be erroneous. This raises the possibility – previously rejected on the basis of twin study results [35] – that genetic variants are important in determining risk only for the relatively rare familial forms of cancer. If so, genetic models of familial aggregation (for example [40]) may be incorrect and the hunt for additional susceptibility genes could be largely fruitless. Existing published findings might then reflect prevailing bias, rather than true associations [14].

Secondly, the model confirms that the potential for reducing the incidence of common diseases using environmental/lifestyle interventions targeted by genotype may be limited [7] by:

(i) the low importance of genetic differences in determining the risk of some conditions (for example, lung cancer);

(ii) the complexity of gene-gene and gene-environment interactions and/or lack of knowledge of environmental factors (for example, schizophrenia).

Targeting environmental/lifestyle interventions at those at 'high genotypic risk' can be of high utility only in specific circumstances. The utility of targeting environmental interactions by genotype (compared to randomly selecting the same number of people from the population) is zero if there is no gene-environment interaction. Utility can also be negative in the presence of a negative interaction (i.e. if the people at high genotypic risk have *less to gain* by the intervention than people at low genotypic risk). The finding that utility increases with gene-environment interaction is consistent with Khoury and Wagener [19] but the relationship is considerably clarified by the adoption here of different measures of the population attributable fraction associated with a targeted intervention (PAF^{E}_{ge}) and of utility (U_{ge}). Further, by formally introducing constraints on the model (for example, that risks are positive and do not exceed 100%), it is possible to demonstrate that both the gene-environment interaction factor and utility have maximum values, which cannot be exceeded for a given data set.

The lung cancer example is apparently trivial but also of critical importance. The R_{MD} value for lung cancer is close to 1, and neither the Scandinavian data used here [31], nor earlier US studies [41], have identified a significant heritable component. It follows from Equation (27) that if the genetic component of the variance, V_{g}, is zero, V_{ge} (the G-E component of the variance) is also zero and using genotyping to target an intervention such as smoking cessation is therefore of zero utility (no better than randomly selecting the same number of individuals). This approximate conclusion is confirmed by the results presented for lung cancer, which show extremely low utility. The detailed calculations may at first sight seem unnecessary, particularly because smoking causes multiple diseases and targeting smoking cessation on the basis of lung cancer risk alone is therefore ill-advised. However, the idea that a genetic test will one day predict which smokers get lung cancer has been widely promoted in the literature and has driven much research aimed at identifying the supposed 'genes for lung cancer' [42]. The results presented here strongly suggest that there will never be a genetic or genotypic test that predicts which smokers will get lung cancer, because the genetic component of the variance is not high enough.

Finally, the model illustrates the argument of Terwilliger and Weiss [11] that the potential for population biobanks to quantify risks for complex disease is limited by a 'multiple testing' problem caused by the large number of genetic and gene-environment interaction models that could fit existing data. Each point in each solution space described above represents a different combination of a genetic risk model (defined by p^{DZ}_{g}) and a G-E interaction model (defined by f_{ge}). Further, any given value of p^{DZ}_{g} may be obtained by an infinite number of different combinations of different alleles acting through multiple biological pathways. Because the number of hypotheses that could be tested is essentially infinite, sample sizes necessary to quantify the risks (R_{oo}, R_{go}, R_{oe} and R_{ge}) could *"plausibly be larger than the number of people that have ever lived*" [11].

The model has several limitations. Measurements of shared sibling risk (λ_{sib}) are needed from the same population as twin data, and the scoping studies are only valid for λ_{DZ} > λ_{sib}, implying that environmental risks are more strongly correlated in dizygotic twins than other siblings. Some data exist to support this assumption for smoking [43] but for other exposures its validity is usually unknown. However, the model does not reduce to the classical twin study solution if this condition is not met: instead, data from more relatives are needed. In principle the model could, and should, be expanded to include data from more relatives, other data (such as migration study data), more risk categories and error terms. However, the number of unknown parameters will then increase, unless more data are available to quantify exposures (which change from generation to generation) and to estimate the extent to which environments are correlated between different types of relative.

Treating exposure and environmental variance (or population attributable fraction) as input data is also problematic when the effects of environmental factors on risk are often unknown. Further, the simple nature of the model (with one environmental axis) cannot adequately represent the complexity of environmental (including socio-economic) causes of disease. However, if targeting environmental interventions by genotype is to be considered, this implies that at least something is known (or expected to be learned) about environmental factors, such as particular exposures, that are amenable to intervention.

The assumption of no gene-environment correlation will often hold (for example it is rather implausible that the same genes strongly influence both lung cancer risk and nicotine addiction), but is not necessarily always true. Adult lactose intolerance is an example of a condition with a strong gene-environment interaction where targeted intervention to avoid drinking milk may be of high utility. However, the model is invalid for lactose intolerance unless exposures are applied equally to the population studied because, in general, people who are lactose intolerant may be less likely to drink milk (a gene-environment correlation) owing to the unpleasant symptoms.

A more fundamental problem is caused by the assumptions that: (i) the risks R_{oo}, R_{go}, R_{oe} and R_{ge} are inherent properties of a given trait within a given population (with a given γ and ε) and that there are therefore no confounders; and (ii) risks are randomly distributed within these categories.

These assumptions, although often made, are implausible in many situations. The assumption of no confounders means that the model can only represent a subset of the potential models of gene-gene and gene-environment interaction described by more complex models (for example [17]). It is unlikely to be met if multiple genetic factors interact with multiple environmental ones [44]. Although this may well render the results presented here invalid, such complexity is likely to reduce the utility of targeting by genotype, rather than enhance it. Hence, situations where the 'no confounders' assumption at least approximately holds are those most likely to be of relevance to public health.

The second assumption neglects the fact that for most exposures there is a gradient in risk, with higher exposure meaning higher risk, and that the same may also be true of genetic factors. This means that increasing the number of categories in the model will increase V_{e} (see [45]) and perhaps V_{g}. Further, these subcategories may be differently correlated between relatives (for example, the twin of a heavy smoker may be more likely to be a heavy smoker than a light one). If so, a relative of a proband may not be representative of their allocated risk category in the four-category model and Equation (22) then becomes invalid.

More broadly, these assumptions make the model, like the classical twin model, essentially deterministic: it assumes that all the factors contributing to correlations in risk between relatives are perfectly known and are either environmental or genetic. Retention of these assumptions here may be problematic and could limit the applicability of the results. Nevertheless, all the other questionable assumptions of the classical twin model have been simultaneously removed.

## Conclusion

The model shows that the potential for reducing the incidence of common diseases using environmental interventions targeted by genotype may be limited, except in special cases. The model also confirms that the importance of an individual's genotype in determining their risk of complex diseases tends to be exaggerated by the classical twin studies method, owing to the 'equal environments' assumption and the assumption of no gene-environment interaction. In addition, if phenotypes are genetically robust, because of epistasis, a largely environmental explanation for shared sibling risk is plausible, even if the classical heritability is high. The model therefore highlights the possibility – previously rejected on the basis of twin study results – that inherited genetic variants are important in determining risk only for the relatively rare familial forms of diseases such as breast cancer. If so, genetic models of familial aggregation may be incorrect and the hunt for additional susceptibility genes could be largely fruitless.

## Appendix A: formal derivation of equation (31)

Equation (23) may be derived more formally by extending the matrix method of Li and Sacks [46].

Define the probability that an affected proband is in genotypic risk category z and environmental risk category w as P_{zw} and assume that risks are randomly distributed within categories. Using the definitions of the four category model given in Table 1, a vector **P** may be defined:

P=\left(\begin{array}{c}{P}_{oo}\\ {P}_{oe}\\ {P}_{go}\\ {P}_{ge}\end{array}\right)=\left(\begin{array}{c}(1-\epsilon )(1-\gamma ){R}_{oo}/{r}_{t}\\ \epsilon (1-\gamma ){R}_{oe}/{r}_{t}\\ \gamma (1-\epsilon ){R}_{go}/{r}_{t}\\ \gamma \epsilon {R}_{ge}/{r}_{t}\end{array}\right)\left(A1\right)

A risk vector **R** may also be defined:

R=\left(\begin{array}{c}{R}_{oo}\\ {R}_{oe}\\ {R}_{go}\\ {R}_{ge}\end{array}\right)\left(A2\right)

Now define G_{xy} as the conditional probability P(relative is in genotypic risk category y|proband is in genotypic risk category x). Similarly, define E_{xy} as the conditional probability P(relative is in environmental risk category y|proband is in environmental risk category x). Using the definitions of p^{rel}_{g} and p^{rel}_{e} given in Section 2.5, matrices **G** and **E** may be written such that:

{G}^{rel}=\left(\begin{array}{cc}{G}_{oo}& {G}_{og}\\ {G}_{go}& {G}_{gg}\end{array}\right)=\left(\begin{array}{cc}{p}_{g}^{rel}+(1-\gamma )(1-{p}_{g}^{rel})& \gamma (1-{p}_{g}^{rel})\\ (1-\gamma )(1-{p}_{g}^{rel})& {p}_{g}^{rel}+\gamma (1-\gamma ){p}_{g}^{rel}\end{array}\right)\left(A3\right)

{E}^{rel}=\left(\begin{array}{cc}{E}_{oo}& {E}_{oe}\\ {E}_{eo}& {E}_{ee}\end{array}\right)=\left(\begin{array}{cc}{p}_{e}^{rel}+(1-\epsilon )(1-{p}_{e}^{rel})& \epsilon (1-{p}_{e}^{rel})\\ (1-\epsilon )(1-{p}_{e}^{rel})& {p}_{e}^{rel}+\epsilon (1-\epsilon ){p}_{e}^{rel}\end{array}\right)\left(A4\right)

Finally, define X_{ab-cd} as the conditional probability P(relative is in risk category cd|proband is in risk category ab), where the risk categories are as defined in Table 1 (for example risk categorgy 'ge' implies high-genotypic and high-environmental risk). Provided p^{rel}_{g} and p^{rel}_{e} are independent (there are no gene-environment correlations), the gene-environment interaction matrix **M**^{rel}_{ge} may be written as:

{M}_{ge}^{rel}=\left(\begin{array}{cccc}{X}_{oo-oo}& {X}_{oo-oe}& {X}_{oo-go}& {X}_{oo-ge}\\ {X}_{oe-oo}& {X}_{oe-oe}& {X}_{oe-go}& {X}_{oe-ge}\\ {X}_{go-oo}& {X}_{go-oe}& {X}_{go-go}& {X}_{go-ge}\\ {X}_{ge-oo}& {X}_{ge-oe}& {X}_{ge-go}& {X}_{ge-ge}\end{array}\right)=\left(\begin{array}{cccc}{G}_{oo}{E}_{oo}& {G}_{oo}{E}_{oe}& {G}_{og}{E}_{oo}& {G}_{og}{E}_{oe}\\ {G}_{oo}{E}_{eo}& {G}_{oo}{E}_{ee}& {G}_{og}{E}_{eo}& {G}_{og}{E}_{ee}\\ {G}_{go}{E}_{oo}& {G}_{go}{E}_{oe}& {G}_{gg}{E}_{oo}& {G}_{gg}{E}_{oe}\\ {G}_{go}{E}_{eo}& {G}_{go}{E}_{ee}& {G}_{gg}{E}_{eo}& {G}_{gg}{E}_{ee}\end{array}\right)\left(A5\right)

Then the risk in a relative of the proband is given by:

{\lambda}_{rel}{r}_{t}=P.\left({M}_{ge}^{rel}R\right)\left(A6\right)

After some algebra, this yields equation (23).

## Appendix B: calculating recurrence risks for twins

The sibling recurrence risk λ_{sib} is often available directly from familial studies. For twins the recurrence risks, if not reported, may be calculated from the case-wise concordance (Cc):

*λ*_{
MZ
}= *Cc*_{
MZ
}/*r*_{
t
} (B1)

*λ*_{
DZ
}= *Cc*_{
DZ
}/*r*_{
t
} (B2)

where, if there is complete ascertainment of all affected twins in a population,

*Cc* = 2*C*/(2*C* + *D*) (B3)

and C is the number of concordant and D the number of discordant pairs [25].

## References

Collins FS: Shattuck Lecture – medical and societal consequences of the Human Genome Project. New Engl J Med. 1999, 341: 28-37. 10.1056/NEJM199907013410106.

Bell J: The new genetics in clinical practice. BMJ. 1998, 316 (7131): 618-620.

Collins FS, McKusick VA: Implications of the Human Genome Project for medical science. J Am Med Assoc. 2001, 285: 540-544. 10.1001/jama.285.5.540.

Strohman RC: The coming Kuhnian revolution in biology. Nat Biotechnol. 1997, 15: 194-200. 10.1038/nbt0397-194.

Holtzman NA, Marteau TM: Will genetics revolutionize medicine?. New Engl J Med. 2000, 343: 141-144. 10.1056/NEJM200007133430213.

Vineis P, Schulte P, McMichael AJ: Misconceptions about the use of genetic tests in populations. Lancet. 2001, 357: 709-712. 10.1016/S0140-6736(00)04136-2.

Baird P: The Human Genome Project, genetics and health. Community Genet. 2001, 4: 77-80. 10.1159/000051161.

Cooper RS, Psaty BM: Genetics and medicine: distraction, incremental progress, or the dawn of a new age?. Ann Intern Med. 2003, 138: 576-580.

Vineis P, Ahsan H, Parker M: Genetic screening and occupational and environmental exposures. Occup Environ Med. 2004, 62: 657-662. 10.1136/oem.2004.019190.

Khoury MJ, Yang Q, Gwinn M, Little J, Flanders WD: An epidemiologic assessment of genetic profiling for measuring susceptibility to common diseases and targeting interventions. Genet Med. 2004, 6 (1): 38-47.

Terwilliger JD, Weiss KM: Confounding, ascertainment bias, and the blind quest for a genetic 'fountain of youth'. Ann Med. 2003, 35: 532-544. 10.1080/07853890310015181.

Ioannidis JPA, Ntzani EE, Trikalinos TA, Contopoulos-Ionnidis DG: Replication validity of genetic association studies. Nat Genet. 2001, 29: 306-309. 10.1038/ng749.

Cordell HJ, Clayton DG: Genetic association studies. Lancet. 2005, 366: 1121-1131. 10.1016/S0140-6736(05)67424-7.

Ioannidis J: Why most published research findings are false. PloS Med. 2005, 2 (8): e124-10.1371/journal.pmed.0020124. DOI: 10.137/journal.pmed.0020124

Layzer D: Heritability analyses of IQ scores: science or numerology?. Science. 1974, 183: 1259-1266. 10.1126/science.183.4131.1259.

Risch N: Linkage strategies for genetically complex traits. I. Multilocus models. Am J Hum Genet. 1990, 46: 222-228.

Guo S-W: Gene-environment interaction and the mapping of complex traits: some statistical models and their interpretation. Hum Hered. 2000, 50: 286-303. 10.1159/000022931.

Hopper JL: Why 'common environmental effects' are so uncommon in the literature. Advances in twin and sib-pair analysis. Edited by: Spector TD, Sneider H, MacGregor AJ. 2000, London: Greenwich Medical Media Ltd

Khoury MJ, Wagener DK: Epidemiological evaluation of the use of genetics to improve the predictive value of disease risk factors. Am J Hum Genet. 1995, 56: 835-844.

Lewis SJ, Brunner EJ: Methodological problems in genetic association studies of longevity – the apolipoprotein E gene as an example. Int J Epidemiol. 2004, 33: 962-970. 10.1093/ije/dyh214.

Tryggvadottir L, Sigvaldason H, Olafsdottir GH, Jonasson JG, Jonsson T, Tulinius H, Eyfjord JE: Population-based study of changing breast cancer risk in Icelandic BRCA2 mutation carriers, 1920–2000. J Natl Cancer Inst. 2006, 98 (2): 116-122.

Humphries S, Ridker PM, Talmud PJ: Genetic testing for cardiovascular disease susceptibility: a useful clinical management tool or possible misinformation?. Arterioscler Thromb Vasc Biol. 2004, 24: 628-636. 10.1161/01.ATV.0000116216.56511.39.

Rose G: Sick individuals and sick populations. Int J Epidemiol. 1985, 14: 32-38.

Khoury MJ, Jones K, Grosse SD: Quantifying the health benefits of genetic tests: The importance of a population perspective. Genet Med. 2006, 8 (3): 191-195.

MacGregor AJ: Practical approaches to account for bias and confounding in twin data. Advances in twin and sib-pair analysis. Edited by: Spector TD, Sneider H, MacGregor AJ. 2000, London: Greenwich Medical Media Ltd

Fisher RA: The correlation between relatives on the supposition of Mendelian inheritance. Trans R Soc Edinb. 1918, 52: 399-433.

Hopper JL: Variance components for statistical genetics: applications in medical research to characteristics related to human diseases and health. Stat Methods Med Res. 1993, 2: 199-223.

Porteous JW: A rational treatment of Mendelian genetics. Theor Biol Med Model. 2004, 1: 6-10.1186/1742-4682-1-6. DOI: 10.1186/1742-4682-1-6

Azevedo RBR, Lohaus R, Srinivasan S, Dang KK, Burch CL: Sexual reproduction selects for robustness and negative epistasis in artificial gene networks. Nature. 2006, 440: 87-90. 10.1038/nature04488.

Risch N: The genetic epidemiology of cancer: interpreting family and twin studies and their implications for molecular genetic approaches. Cancer Epidemiol Biomarkers Prev. 2001, 10: 733-741.

Lichtenstein P, Holm NV, Verkasalo PK, Iliadou A, Kaprio J, Koskenvuo M, Pukkala E, Skytthe A, Hemminki K: Environmental and heritable factors in the causation of cancer. New Engl J Med. 2000, 343: 78-85. 10.1056/NEJM200007133430201.

Dong C, Hemminki K: Modification of cancer risks in offspring by sibling and parental cancers from 2,112,616 nuclear families. Int J Cancer. 2001, 92: 144-150. 10.1002/1097-0215(200102)9999:9999<::AID-IJC1147>3.0.CO;2-C.

Rockhill B, Weinberg CR, Newman B: Population attributable fraction estimation for established breast cancer risk factors: considering the issues of high prevalence and unmodifiability. Am J Epidemiol. 1998, 147 (9): 826-833.

McGue M, Gottesman II, Rao DC: The transmission of schizophrenia under a multifactorial threshold model. Am J Hum Genet. 1983, 35: 1161-1178.

Easton DF: How many more breast cancer predisposition genes are there?. Breast Cancer Res. 1999, 1 (1): 14-17. 10.1186/bcr6.

Badano JL, Katsanis N: Beyond Mendel: an evolving view of human genetic disease transmission. Nat Rev Genet. 2002, 3: 779-789. 10.1038/nrg910.

Badano JL, Leitch CC, Ansley SJ, May-Simera H, Lawson S, Lewis RA, Beales PL, Dietz HC, Fisher S, Katsanis N: Dissection of epistasis in oligogenic Bardet-Biedl syndrome. Nature. 2006, 439: 326-330. 10.1038/nature04370.

Taioli E, Zocchetti C, Garte S: Models of interaction between metabolic genes and environmental exposure in cancer susceptibility. Environ Health Perspect. 1998, 106 (2): 67-70.

Hunter DJ: Gene-environment interactions in human diseases. Nat Rev Genet. 2005, 6: 287-298. 10.1038/nrg1578.

Antoniou AC, Pharoah PDP, McMullan G, Day NE, Stratton MR, Peto J, Ponder BJ, Easton DF: A comprehensive model for familial breast cancer incorporating BRCA1, BRCA2 and other genes. Br J Cancer. 2002, 86: 76-83. 10.1038/sj.bjc.6600008.

Braun MM, Caporaso NE, Page WF, Hoover RN: A cohort study of twins and cancer. Cancer Epidemiol Biomarkers Prev. 1995, 4 (5): 469-473.

Hall W, Madden P, Lynskey M: The genetics of tobacco use: methods, findings and policy implications. Tob Control. 2002, 11: 119-124. 10.1136/tc.11.2.119.

Vink JM, Willemsen G, Boomsma DI: The association of current smoking behavior with the smoking behavior of parents, siblings, friends and spouses. Addiction. 2003, 98: 923-931. 10.1046/j.1360-0443.2003.00405.x.

Taioli E, Garte S: Covariates and confounding in epidemiologic studies using metabolic gene polymorphisms. Int J Cancer. 2002, 100: 97-100. 10.1002/ijc.10448.

Guo S: The behaviors of some heritability estimators in the complete absence of genetic factors. Hum Hered. 1999, 49 (4): 215-228. 10.1159/000022878.

Li CC, Sacks L: The derivation of joint distribution and correlation between relatives by the use of stochastic matrices. Biometrics. 1954, 10: 347-360. 10.2307/3001590.

## Acknowledgements

The author is grateful to the Joseph Rowntree Charitable Trust for funding the completion of this work.

## Author information

### Authors and Affiliations

### Corresponding author

## Additional information

### Competing interests

The author(s) declare that they have no competing interests.

## Electronic supplementary material

### 12976_2006_98_MOESM1_ESM.xls

Additional File 1: Gene-gene and gene-environment interaction model. Contains the Visual Basic macro (Twincal), input and output datasheets and charts used to calculate the solutions described in the text. The program is run by entering parameters in the 'Inputs' sheet and clicking on the 'Run' button. Note that for the final chart ('fe') the number of categories on the horizontal axis changes depending on the environmental input parameters ε and PAF^{E}_{e}. If these parameters are changed it is therefore necessary to delete the lower part of the output sheet prior to running the model and, after the run, to redraw the chart using the source data option from the chart. All other charts are drawn automatically. The line γ_{max} = γ_{min} is calculated exactly for the chart 'fe' but is approximated in the charts 'pgdz' and 'pgdzneg' using Newton's method and an initial guess for f_{ge} (f0) and step (fet). For some input parameters it may be necessary to change these values by editing the Visual Basic code (Twincal) to obtain a valid solution. (XLS 814 KB)

### 12976_2006_98_MOESM2_ESM.bmp

Additional File 2: Supplementary Figure 1: Example model solution space with R_{MD} = 1.7 and U_{ge} ≥ 0. Model solution space with U_{ge} ≥ 0 for the same input parameters as Figure 2, apart from λ_{MZ} = 4.4. (BMP 73 KB)

### 12976_2006_98_MOESM3_ESM.bmp

Additional File 3: Supplementary Figure 2: Example model solution space with R_{MD} = 1.8 and U_{ge} ≥ 0. Model solution space with U_{ge} ≥ 0 for the same input parameters as Figure 2, apart from λ_{MZ} = 4.6. (BMP 73 KB)

### 12976_2006_98_MOESM4_ESM.bmp

Additional File 4: Supplementary Figure 3: Example model solution space with R_{MD} = 1.95 and U_{ge} ≥ 0. Model solution space with U_{ge} ≥ 0 for the same input parameters as Figure 2, apart from λ_{MZ} = 4.9. (BMP 73 KB)

### 12976_2006_98_MOESM5_ESM.bmp

Additional File 5: Supplementary Figure 4: Example model solution space with R_{MD} = 2.1 and U_{ge} ≥ 0. Model solution space with U_{ge} ≥ 0 for the same input parameters as Figure 2, apart from λ_{MZ} = 5.2. (BMP 73 KB)

### 12976_2006_98_MOESM6_ESM.bmp

Additional File 6: Supplementary Figure 5: Example full solution space with R_{MD} = 1.7. Full model solution space for the same input parameters as Figure 5, transformed so that f_{ge} is on the horizontal axis. (BMP 71 KB)

### 12976_2006_98_MOESM7_ESM.bmp

Additional File 7: Supplementary Figure 6: Example full solution space with R_{MD} = 1.8. Full model solution space for the same input parameters as Figure 6, transformed so that f_{ge} is on the horizontal axis. (BMP 71 KB)

### 12976_2006_98_MOESM8_ESM.bmp

Additional File 8: Supplementary Figure 7: Example full solution space with R_{MD} = 1.95. Full model solution space for the same input parameters as Figure 7, transformed so that f_{ge} is on the horizontal axis. (BMP 71 KB)

### 12976_2006_98_MOESM9_ESM.bmp

Additional File 9: Supplementary Figure 8: Example full solution space with R_{MD} = 2.1. Full model solution space for the same input parameters as Figure 8, transformed so that f_{ge} is on the horizontal axis. (BMP 71 KB)

### 12976_2006_98_MOESM10_ESM.bmp

Additional File 10: Supplementary Figure 9: Example model solution with c_{MD} > 1 and U_{ge} ≥ 0. Input parameters: λ_{MZ} = 5.2, λ_{DZ} = 3, λ_{sib} = 2, ε = 0.2, PAF^{E}_{e} = 0.5, c_{MD} = 2, r_{t} = 0.1. (BMP 2 MB)

### 12976_2006_98_MOESM11_ESM.bmp

Additional File 11: Supplementary Figure 10: Example model solution with c_{MD} > 1 and U_{ge} ≥ 0. Input parameters as for Figure 13. (BMP 88 KB)

### 12976_2006_98_MOESM12_ESM.bmp

Additional File 12: Supplementary Figure 11: Example full solution space with c_{MD} > 1. Full model solution space for the same parameters as Figure 13, transformed so that f_{ge} is on the horizontal axis. (BMP 71 KB)

### 12976_2006_98_MOESM13_ESM.bmp

Additional File 13: Supplementary Figure 12: Breast cancer solution space with U_{ge} ≥ 0. Input parameters are as shown in Table 5, with c_{MD} = 1. The solution space is shown (shaded) for positive f_{ge}, assuming the 'equal environments' assumption holds (c_{MD} = 1). The darker shaded area shows the part of the solution space for which γ_{maxge} < 1/2. Utility U_{ge} is at its maximum when γ = 1/2 except within this darker shaded area. (BMP 75 KB)

### 12976_2006_98_MOESM14_ESM.bmp

Additional File 14: Supplementary Figure 13: Breast cancer variances with f_{ge} = 0. Input parameters as for Figure 16. Additive model of G-E interaction (f_{ge} = 0). Variance components are genetic (V_{g}) or environmental (V_{e}). (BMP 66 KB)

### 12976_2006_98_MOESM15_ESM.bmp

Additional File 15: Supplementary Figure 14: Breast cancer variances with f_{ge} = 1. Input parameters as for Figure 16. Multiplicative G-E interaction model (f_{ge} = 1). Variance components are genetic (V_{g}), environmental (V_{e}) or due to gene-environment interaction (V_{ge}). (BMP 66 KB)

### 12976_2006_98_MOESM16_ESM.bmp

Additional File 16: Supplementary Figure 15: Breast cancer variances with f_{ge} = 1/PAF^{E}_{e}. Input parameters as for Figure 16. Maximum G-E interaction model (f_{ge} = 1/PAF^{E}_{e}). Variance components are genetic (V_{g}), environmental (V_{e}) or due to gene-environment interaction (V_{ge}). (BMP 66 KB)

### 12976_2006_98_MOESM17_ESM.bmp

Additional File 17: Supplementary Figure 16: Breast cancer γ values with f_{ge} = 0. Input parameters as for Figure 16. The proportion of the population in the 'high genotypic risk' group, γ, may take any value in the shaded area. γ_{min} occurs when R_{ge} = 1, i.e. when the Positive Predictive Value (PPV) of being in the 'ge' subgroup is 100%. γ_{max} occurs when R_{oo} = 1 for an additive G-E model and solutions with a Population Impact of 100% (PI = 1) cannot exist. (BMP 67 KB)

### 12976_2006_98_MOESM18_ESM.bmp

Additional File 18: Supplementary Figure 17: Breast cancer γ values with f_{ge} = 1. Input parameters as for Figure 16. The proportion of the population in the 'high genotypic risk' group, γ, may take any value in the shaded area. A solution with a Population Impact of 100% (PI = 1) may exist if γ = γ_{max}. (BMP 62 KB)

### 12976_2006_98_MOESM19_ESM.bmp

Additional File 19: Supplementary Figure 18: Breast cancer γ values with f_{ge} = 1/PAF^{E}_{e}. Input parameters as for Figure 16. The proportion of the population in the 'high genotypic risk' group, γ, may take any value in the shaded area. A solution with a Population Impact of 100% (PI = 1) may exist if γ = γ_{max}. (BMP 64 KB)

### 12976_2006_98_MOESM20_ESM.bmp

Additional File 20: Supplementary Figure 19: Breast cancer solution space with U_{ge} ≤ 0. Input parameters are as for Figure 16. The solution space is shown for negative f_{ge} (where the utility of targeting environmental interventions at the high genotypic risk group is negative, U_{ge} ≤ 0). Solutions exist only in the shaded area where γ_{max} ≥ γ_{min}. (BMP 59 KB)

### 12976_2006_98_MOESM21_ESM.bmp

Additional File 21: Supplementary Figure 20: Breast cancer: full solution space. Input parameters are as for Figure 16. The same solution space as Figures 16 and 23 is shown (shaded), transformed so that the G-E interaction factor is plotted on the horizontal axis. Again, each point in the shaded solution space represents a genetic model defined by p^{DZ}_{g} and a G-E interaction model defined by f_{ge}. The area of solutions with γ_{maxge} < 1/2 is highlighted with darker shading. The classical twin study solution lies on the vertical axis (f_{ge} = 0) at the point p^{DZ}_{g} = 1/2, and is slightly outside the solution space. (BMP 73 KB)

### 12976_2006_98_MOESM22_ESM.bmp

Additional File 22: Supplementary Figure 21: Lung cancer solution space with U_{ge} ≥ 0. Input parameters are as shown in Table 5, with c_{MD} = 1. (BMP 59 KB)

### 12976_2006_98_MOESM23_ESM.bmp

Additional File 23: Supplementary Figure 22: Lung cancer variances with f_{ge} = 0. Input parameters as for Figure 25. Note that the horizontal axis has been expanded to show high values of c_{SD}/R_{SD} only. (BMP 66 KB)

### 12976_2006_98_MOESM24_ESM.bmp

Additional File 24: Supplementary Figure 23: Lung cancer variances with f_{ge} = 1. Input parameters as for Figure 25. Note that the horizontal axis has been expanded to show high values of c_{SD}/R_{SD} only. (BMP 63 KB)

### 12976_2006_98_MOESM25_ESM.bmp

Additional File 25: Supplementary Figure 24: Lung cancer variances with f_{ge} = 1/PAF^{E}_{e}. Input parameters as for Figure 25. Note that the horizontal axis has been expanded to show high values of c_{SD}/R_{SD} only. (BMP 63 KB)

### 12976_2006_98_MOESM26_ESM.bmp

Additional File 26: Supplementary Figure 25: Lung cancer γ values for f_{ge} = 1. Input parameters as for Figure 25. The proportion of the population in the 'high genotypic risk' group, γ, may take any value in the shaded area. (BMP 59 KB)

### 12976_2006_98_MOESM27_ESM.bmp

Additional File 27: Supplementary Figure 26: Lung cancer γ values for f_{ge} = 1/PAF^{E}_{e}. Input parameters as for Figure 25. The proportion of the population in the 'high genotypic risk' group, γ, may take any value in the shaded area. (BMP 63 KB)

### 12976_2006_98_MOESM28_ESM.bmp

Additional File 28: Supplementary Figure 27: Lung cancer U_{ge} values for f_{ge} = 1. Input parameters as for Figure 25. The utility parameter, U_{ge}, may take any value in the shaded area, but is maximum when γ = 1/2. (BMP 63 KB)

### 12976_2006_98_MOESM29_ESM.bmp

Additional File 29: Supplementary Figure 28: Lung cancer U_{ge} values for f_{ge} = 1/PAF^{E}_{e}. Input parameters as for Figure 25. The utility parameter, U_{ge}, may take any value in the shaded area, but is maximum when γ = 1/2. (BMP 59 KB)

### 12976_2006_98_MOESM30_ESM.bmp

Additional File 30: Supplementary Figure 29: Lung cancer: full solution space. Input parameters as for Figure 25. (BMP 2 MB)

### 12976_2006_98_MOESM31_ESM.bmp

Additional File 31: Supplementary Figure 30: Schizophrenia U_{ge} ≥ 0, small environmental variance and c_{MD} ≥ 1. Input parameters are as shown in Table 5, with ε = 0.62, PAF^{E}_{e} = 0.15 and c_{MD} = 1. (BMP 59 KB)

### 12976_2006_98_MOESM32_ESM.bmp

Additional File 32: Supplementary Figure 31: Schizophrenia U_{ge} ≥ 0, small environmental variance and c_{MD} > 1. Input parameters are as shown in Table 5, with ε = 0.62, PAF^{E}_{e} = 0.15 and c_{MD} = 3.8. (BMP 59 KB)

### 12976_2006_98_MOESM33_ESM.bmp

Additional File 33: Supplementary Figure 32: Schizophrenia U_{ge} ≥ 0, large environmental variance and c_{MD} = 1. Input parameters are as shown in Table 5, with ε = 0.15, PAF^{E}_{e} = 0.86 and c_{MD} = 1. (BMP 59 KB)

## Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

## Rights and permissions

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.

## About this article

### Cite this article

Wallace, H.M. A model of gene-gene and gene-environment interactions and its implications for targeting environmental interventions by genotype.
*Theor Biol Med Model* **3**, 35 (2006). https://doi.org/10.1186/1742-4682-3-35

Received:

Accepted:

Published:

DOI: https://doi.org/10.1186/1742-4682-3-35

### Keywords

- Environmental Intervention
- Dizygotic Twin
- Population Attributable Fraction
- Genotypic Risk
- Population Impact