On the Use of Variance per Genotype as a Tool to Identify Quantitative Trait Interaction Effects: A Report from the Women's Genome Health Study
Testing for genetic effects on mean values of a quantitative trait has been a very successful strategy. However, most studies to date have not explored genetic effects on the variance of quantitative traits as a relevant consequence of genetic variation. In this report, we demonstrate that, under plausible scenarios of genetic interaction, the variance of a quantitative trait is expected to differ among the three possible genotypes of a biallelic SNP. Leveraging this observation with Levene's test of equality of variance, we propose a novel method to prioritize SNPs for subsequent gene–gene and gene–environment testing. This method has the advantageous characteristic that the interacting covariate need not be known or measured for a SNP to be prioritized. Using simulations, we show that this method has increased power over exhaustive search under certain conditions. We further investigate the utility of variance per genotype by examining data from the Women's Genome Health Study. Using this dataset, we identify new interactions between the LEPR SNP rs12753193 and body mass index in the prediction of Creactive protein levels, between the ICAM1 SNP rs1799969 and smoking in the prediction of soluble ICAM1 levels, and between the PNPLA3 SNP rs738409 and body mass index in the prediction of soluble ICAM1 levels. These results demonstrate the utility of our approach and provide novel genetic insight into the relationship among obesity, smoking, and inflammation.
Published in the journal:
. PLoS Genet 6(6): e32767. doi:10.1371/journal.pgen.1000981
Category:
Research Article
doi: 10.1371/journal.pgen.1000981
Summary
Testing for genetic effects on mean values of a quantitative trait has been a very successful strategy. However, most studies to date have not explored genetic effects on the variance of quantitative traits as a relevant consequence of genetic variation. In this report, we demonstrate that, under plausible scenarios of genetic interaction, the variance of a quantitative trait is expected to differ among the three possible genotypes of a biallelic SNP. Leveraging this observation with Levene's test of equality of variance, we propose a novel method to prioritize SNPs for subsequent gene–gene and gene–environment testing. This method has the advantageous characteristic that the interacting covariate need not be known or measured for a SNP to be prioritized. Using simulations, we show that this method has increased power over exhaustive search under certain conditions. We further investigate the utility of variance per genotype by examining data from the Women's Genome Health Study. Using this dataset, we identify new interactions between the LEPR SNP rs12753193 and body mass index in the prediction of Creactive protein levels, between the ICAM1 SNP rs1799969 and smoking in the prediction of soluble ICAM1 levels, and between the PNPLA3 SNP rs738409 and body mass index in the prediction of soluble ICAM1 levels. These results demonstrate the utility of our approach and provide novel genetic insight into the relationship among obesity, smoking, and inflammation.
Introduction
Most genetic association studies have so far emphasized the effect of SNPs on mean values of a trait of interest. While this approach has been widely successful in identifying genetic determinants of quantitative traits, it does not make use of differences in the variance of these traits per genotype, which may reflect an alternative mode of genetic influence and signify biological relevance. Several scenarios can result in the variance of a quantitative trait being different when individuals are stratified by SNP genotype (or any other genetic variation) in crosssectional designs. For examples, such difference could be ascribed to increased variation over time among individuals of a given genotype, increased variation in the homeostatic setpoint of a trait (without increased variation over time), or an interaction effect with another SNP or exposure. This latter possibility is particularly attractive since information on variance per genotype could thus be leveraged to reduce the dimensionality of interaction analysis by prioritizing SNPs for further interaction testing.
Several methods have been developed to address the statistical challenge of assessing the degree and nature of interaction between two or more genetic loci [1]–[4]. While these offer invaluable insight, most do not address the main problem pertaining to interaction testing in the context of whole genome association studies: the multiple hypothesis testing incurred when searching for interactions between hundreds of thousands of genetic variants. In theory, excluding unlikely interactions on an a priori basis will result in a lower number of false positives tests and increased power. Current data reduction methods use evidence of coadaptation between loci to restrict the search space [5] or specify a minimum marginal effect to select a subset of SNPs for interaction testing [6]–[8]. This latter method is based on the assumption that interacting SNPs will show at least a moderate effect on the trait of interest. While this premise might be valid in most cases, it need not be, as has been observed [2], [9], [10]. Furthermore, even when restricting the search for interaction to SNPs with marginal effect Pvalues lower than 0.1, as has been suggested [7], a sizeable number of hypotheses will have to be tested. Thus there is a need for novel methods to identify SNPs likely to be involved in genegene and geneenvironment interactions. Murcray et al. [11] recently proposed a twostep approach to test for geneenvironment interactions in genomewide association studies. In the context of dichotomous outcomes, this method first tests for geneenvironment association, followed by conventional interaction testing of the subset of markers with significant association in the combined casecontrol sample.
In this report, we first derive expressions for the variance according to SNP genotype under models of genetic interactions and demonstrate that the variance of a quantitative trait will differ between the three possible genotypes of a biallelic SNP in the presence of genetic interactions. Building on these derivations, we then simulate data to explore contexts where testing for inequality of variance between genotypes could be used to prioritize SNPs for further interaction testing. Leveraging these observations, we propose a novel method to prioritize SNPs for genegene and geneenvironment interaction, “variance prioritization”. Variance prioritization entails two steps. In the first step, Levene's test of equality of variance is used to prioritize SNPs for further interaction testing. In the second step, prioritized SNPs are tested for interaction effects against environmental covariates or other SNPs using linear regression. We use data from the Women's Genome Health Study (WGHS) on two inflammation markers, Creactive protein and soluble ICAM1, to test whether (1) there is empirical confirmation of SNPs demonstrating inequality of variance between genotypes and (2) whether such SNPs show evidence of either geneenvironment or genegene interactions.
Methods
Ethics statement
All analyses were performed with approval of the institutional review board of the Brigham and Women's Hospital. All members of the WGHS cohort were participants in the WHS who provided an adequate baseline blood sample for plasma and DNA analysis and who gave consent for bloodbased analyses and longterm followup.
Variance of a quantitative trait by genotype in the presence of an interaction effect
Consider a biallelic SNP G with minor allele frequency p and in HardyWeinberg equilibrium. Without loss of generality and assuming an additive genetic effect of G on a quantitative trait of interest, G can be encoded as 2p, 12p, or 22p such that mean genotype equals zero. Assume that there is an additive interaction between G and a covariate C (with mean equal to zero, variance equal to one, and uncorrelated to G) such that the quantitative trait Y for the ith individual is represented as:(1)Where β_{0}, β_{1}, β_{2}, and β_{3} represent constants for the intercept, the marginal effect of G, the marginal effect of C, and the interaction of G and C, respectively. It is assumed the outcome is standardized so it has unit variance conditional on G and C, that is, the error term ε_{i}∼N(0,1). The variance in Y conditional on genotype of G can be written as:(2)Where coefficients represent true rather than estimated effects. Thus, the variance of Y will differ by genotype group (g) when β_{3}≠0. Importantly, the variance does not depend on the marginal effect of G on Y (i.e. β_{1}). Furthermore, as we did not specify the distribution form of the interacting covariate C, Equation 3 can be generalized to different types of interaction such as dichotomous environmental exposures, continuous environmental exposures or the three genotype classes of a biallelic SNP. The total variance of Y is given by:(3)Where the three covariance terms are zero because of the assumption of G and C independence. It results from equation (3) that the proportion of total variance attributable to the marginal effect of the covariate (i.e. β_{2}) can be defined as:(4)Similarly, the proportion of total variance attributable to the SNPcovariate interaction (i.e. β_{3}) can be defined as:(5)
Levene's test for equality of variance
Levene's test [12] is used to test if k samples have equal variances and is less sensitive to deviation from normality than the alternative Bartlett test [13]. Indeed, use of Bartlett test with empirical data led to marked inflation of type I error and was not pursued further (data not shown). The Levene test assesses the null hypothesis that the population variances are equal:

H_{0}:

H_{a}: for at least one pair (i,j).
Given a variable Y with sample of size N divided into k subgroups, where N_{i} is the sample size of the ith subgroup, the Levene test statistic is defined as:Where with the value of Y for the jth observation of the ith subgroup, the mean of the ith subgroup, are the group means of the and is the overall mean of the .
The Levene test rejects the hypothesis that the variances are equal if:Where is the upper critical value of the F distribution with k−1 and N−k degrees of freedom at a significance level of .
Simulations of interaction effects
We used simulations to estimate the power of Levene's test to detect differences in variance of a quantitative trait between individuals of differing genotypes. We first simulated the case of a SNP (denoted SNP_{1}) interacting with a covariate (denoted C_{1}) in determining the value of a quantitative trait (denoted Y; see Figure 1). We assumed the main effect of SNP_{1} to be null in all simulations because it is expected to affect neither the within genotype variance (see Equation 2) nor the power to detect the interaction. For simplicity, only results from the additive genetic model (assuming an additive contribution of each minor allele; see Equation 1) are shown but similar conclusions can be derived from more general models of association. Under a range of models for the main effect of C_{1} (on the quantitative trait), allelic frequency of SNP_{1}, and interaction between SNP_{1} and C_{1}, we tested the power of Levene's test to reject the null hypothesis that the variance of Y does not vary between subsets of individuals with differing genotypes (of SNP_{1}). In all simulations, the quantitative trait was assumed to have a normally distributed error term with mean equal to zero and variance equal to 1. The covariate was normally distributed with mean equal to zero and variance equal to 1. Interactions are expected to be modest at best and simulations were therefore restricted to interaction effects explaining up to 1% of the quantitative trait's total variance. Each condition was simulated 1,000 times using 15,000 individuals unless otherwise specified. In a second set of simulations, we tested the power of Levene's test of equality of variance to detect SNPSNP interactions. These simulations are based on a model similar to the one previously introduced, with the exception that the normally distributed covariate (C_{1}) is replaced by a second, independent, SNP (denoted SNP_{2}). To simplify simulations, both SNPs were assigned the same allelic frequency as well as the same main effect (i.e. SNPY association). Varying the strength of the SNP_{1} SNP_{2} interaction, we simulated each condition 1,000 times using 15,000 individuals unless otherwise specified.
Variance prioritization
Leveraging two observations that (1) the signature of a SNPcovariate interaction may be inequality of variance between SNP genotypes and (2) a test for inequality of variance is independent of a subsequent direct test of interaction, we propose a novel method to prioritize SNPs for genegene and geneenvironment interaction, termed “variance prioritization”. This prioritization procedure first applies Levene's test of variance equality to selects SNPs below a prespecified threshold and then tests these “prioritized” SNPs for geneenvironment or genegene interaction using standard linear regression models (see Figure 2), with Bonferroni correction only for the number statistical tests prioritized in the first step. Using simulation models, we explored the power of variance prioritization and compared it to an exhaustive search for interaction, without prior prioritization, followed by Bonferroni correction for multiple hypothesis testing.
Independence of variance per genotype from estimates of interaction effects under the null of no interaction is important because the validity of the method rests on having correct type I error. A proof of independence between variance per genotype and estimation of interaction effect under the null of no interaction is presented in Text S1 and Figure S1.
WGHS participants
The study population was derived from the Women's Genome Health Study (WGHS) [14]. Briefly, participants in the WGHS include American women from the Women's Health Study (WHS) with no prior history of cardiovascular disease, cancer, or other major chronic illness who also provided a baseline blood sample at the time of study enrollment from which genomic DNA was extracted [15]. Individuals with selfreported diabetes at enrollment (or with glycated hemoglobin >7%) or taking lipid lowering medication were excluded from analysis.
All baseline blood samples underwent measurement for highsensitivity CRP (hsCRP) via a validated immunoturbidometric method (Denka Seiken, Tokyo, Japan). CRP is an inflammatory biomarker with strong prognostic value for development of metabolic syndrome, diabetes, myocardial infarction, and stroke. Logtransformed CRP levels were used in all analyses. Soluble ICAM1 is an inflammatory marker that has been associated with several common diseases such as diabetes, heart disease, stroke and malaria. Circulating plasma sICAM1 concentrations were determined using a commercial ELISA assay (R&D Systems, Minneapolis, Minn.) The environmental factors tested for interaction included selfreported body mass index, menopause status, hormone replacement therapy use, smoking, age and exercise, expressed as quintiles of total energy expenditure spent from all recreational activities.
Genotyping
DNA samples were genotyped with the Infinium II technology from Illumina. Either the HumanHap300 DuoPlus chip or the combination of the HumanHap300 Duo and ISelect chips were used. In either case, the custom content was identical and consisted of candidate SNPs chosen without regard to allele frequency to increase coverage of genetic variation with impact on biological function including metabolism, inflammation or cardiovascular diseases. Genotyping at 318,237 HumanHap300 Duo SNPs and 45,571 custom content SNPs was attempted, for a total of 363,808 SNPs. Genetic context for all annotations are derived from human genome build 36.1 and dbSNP build 126.
SNPs with call rates <90% were excluded from further analysis. Likewise, all samples with percentage of missing genotypes higher than 2% were removed. Among retained samples, SNPs were further evaluated for deviation from HardyWeinberg equilibrium using an exact method and were excluded when the Pvalue was lower than 10^{−6}. SNPs with minor allele frequency >1% in Caucasians were used for analysis. 339,596 SNPs were left after quality control. All statistical analyses were done with R, unless specified otherwise. All analyses were performed on individuals of selfreported Caucasian ancestry, with confirmation using a principal component analysis. A total of 21,799 individuals were included in our final dataset.
Results
Building on the theoretical observation that the within genotype variance of a quantitative trait will vary when a genetic interaction is present, we explored the power of Levene's test of equality of variance to detect interactions. We first tested the power of Levene's test to detect a SNP (denoted SNP_{1}) interacting with a single covariate, C_{1} (see Figure 1). As shown in Figure 3, the power of Levene's test to identify a SNP as interacting with a covariate was modest in most situations. For example, when the interaction effect explained 0.4% of the variance of a quantitative trait, the allelic frequency was 0.2 and the covariate was associated with the quantitative trait with β_{2} = 0.3 (explaining 8.2% of the variance; see Figure 3B), the power of Levene's test to identify a SNP as “interacting” at P<0.05 was 82.5% whereas the power of linear regression to identify the interaction at P<1.5×10^{−7} was 99.6%. In other words, even after accounting for 340,000 tests using Bonferroni correction, testing for geneenvironment interactions with linear regression was more powerful than testing for inequality of variance between genotypes at nominal significance.
Nevertheless, as can be observed in Figure 3, there were situations involving small interaction effects explaining less than 0.4% of the variance for which a nominally significant Levene's test was more powerful than linear regression testing (after multiple hypothesis correction). Because SNP effects on variance per genotype is independent of the estimates of interaction effects, variance prioritization could provide a method for decreasing the multiple testing burden in these situations, with ensuing increased power. To test the power of variance prioritization for small interaction effects, we performed further simulations in which we varied Levene's test Pvalue threshold used to prioritize SNPs for interaction testing (see Figure 4). For these simulations, we used the same parameters as in Figure 3B (Figure 4A and 4B) and Figure 3D (Figure 4C and 4D), fixing β _{3} such that the interaction explained 0.1% or 0.15% of total variance. We subsequently interrogated the prioritized SNPs for interaction effects using linear regression and applied a Bonferroni correction assuming 340,000 SNPs were initially tested (ex.: if Levene's test Pvalue cutoff is 0.01 then upon further interaction testing a SNP will be deemed significant if its interaction Pvalue is <0.05/(0.01×340,000)). We then compared the results of variance prioritization with a “conventional” method correcting for 340,000 tests. Each condition was simulated 2,000 times. As illustrated in Figure 4A–4D, variance prioritization was more powerful than conventional methods to detect geneenvironment interactions, with the increase in power dependent on the Levene's test Pvalue threshold used. In almost all cases, the power was at least equal to performing an exhaustive search for geneenvironment interactions with all 340,000 SNPs.
Similar simulations were performed for SNPSNP interactions (see Figure 5, and Figure 4E–4H). For simplicity, the allelic frequency of both SNPs and their main effect on the quantitative trait was set to be equal. To account for all possible pairwise SNPSNP interactions between 340,000 SNPs, the “conventional” Pvalue threshold was set at 4.3×10^{−13}. Figure 5 illustrates the power that either SNP will have a Levene's test Pvalue lower than the prespecified cutoff, and compares it to the power to detect the interaction by exhaustive linear regression testing. Once again, the power to detect either SNP with a nominally significant Levene's test was higher than the power to identify the SNPSNP interaction using exhaustive linear regression testing when the interaction effect was weak, providing the basis for prioritization of SNPs with Levene's test. Results of variance prioritization are shown in Figure 4E and 4F and Figure 4G and 4H, matching the conditions of Figure 5A and 5B, respectively, and fixing B_{3} such that 0.2% or 0.25% of the variance was explained by the interaction. In these simulations, all SNPs with Levene's test Pvalue lower than the specified threshold were prioritized and tested for interaction against all other SNPs, prioritized or not (with appropriate correction for multiple hypothesis testing). This approach was shown (by simulations) to be more powerful than testing for pairwise interactions between Levene's test prioritized SNPs only, which was underpowered as compared to an exhaustive search (data not shown). As illustrated in Figure 4E, the power of variance prioritization was superior to exhaustive search under specific conditions. For example, when a Levene's Pvalue lower than 0.1 was chosen as the prioritization threshold, the power of variance prioritization was 9.05% as compared to 7.30% when testing all possible pairwise interactions between 340,000 SNPs. However, variance prioritization was associated with decreased power in certain conditions, as illustrated in Figure 4G and 4H.
To empirically test whether some SNPs exhibit inequality of variance between genotypes and whether these effects might be due to interactions, we applied Levene's test of equality of variance to ongoing genome scans of two inflammatory markers, Creactive protein (CRP) and soluble ICAM1 (sICAM1). Only SNPs with genomewide significant Levene's Pvalue were selected for subsequent study in order to ensure they represent true positives. From the analysis of our sample of 21,799 women, 23 SNPs had a Levene's test Pvalue lower than 1.5×10^{−7} for CRP (all clustering near the LEPR gene on chromosome 1). Because the LEPR SNPs are in linkage disequilibrium, only rs12753193 will be further considered since it had the most significant Levene's test Pvalue. Two SNPs, the PNPLA3 SNP rs738409 and the ICAM1 SNP rs1799969, had genomewide significant Levene's test Pvalue for sICAM1 (see Figure 6 and Table S1). To further characterize these SNPs, we then tested all pairwise interactions with the 339,595 remaining SNPs as well as environmental covariates, using linear regression. No SNPSNP interaction was observed after correction for multiple hypothesis testing. Three SNPenvironment interactions were significant after correction for multiple hypothesis testing (P<0.05/(3 SNPs×6 traits); see Table 1). The LEPR SNP rs12753193 interacted with BMI in the prediction of CRP levels (P = 7.2×10^{−10}), the ICAM1 SNP rs1799969 interacted with smoking status in the prediction of sICAM1 levels (P = 4.8×10^{−9}) and the PNPLA3 SNP rs738409 interacted with BMI in the prediction of sICAM1 levels (P = 1.6×10^{−7}).
We sought to determine the contribution of the three interactions to Levene's test of inequality of variance. While Levene's test Pvalue was 1.6×10^{−29} for the LEPR SNP rs12753193 in the CRP analysis, the Pvalue was 9.1×10^{−21} after adjustment for the interaction for BMI. Similarly, Levene's test Pvalues were 2.1×10^{−9} for the ICAM1 SNP rs1799969 and 1.9×10^{−10} for the PNPLA3 SNP rs738409 in the sICAM1 analysis but were slightly less significant after adjustment for the intraction with BMI and smoking status, respectively, at 3.3×10^{−9} and 1.6×10^{−7}. Overall, these interactions explained 0.18% of CRP total variance, and 0.16% and 0.13% of sICAM1 total variance for the ICAM1 and PNPLA3 interactions, respectively. Neither the LEPR SNP rs12753193 nor the PNPLA3 SNP rs738409 was associated with BMI (P>0.05). The ICAM1 SNP rs1799969 was not associated with smoking status (P>0.05). These data collectively suggest that individual interactions are likely to have modest effects, and that additional covariates and/or SNPs might interact with the validated SNP.
Discussion
Genetic effects on the variance of quantitative traits have been largely ignored so far in the genetic analysis. As shown in this report, differences in the variance per genotype can reflect interaction effects and these differences can be leveraged to prioritize SNPs for interaction testing. The proposed method, variance prioritization, presents the distinct advantage that it does not rely on the marginal effect of SNPs on the means to prioritize them. Importantly, it can be used in combination with prioritization based on marginal effects (on the means) if it is assumed all interacting SNPs have a main effect on the trait of interest. In such case, combining variance prioritization to prioritization based on a predefined marginal effect Pvalue provides the greatest power as compared to each procedure individually, if the assumption that interacting SNPs have main effects is true (data not shown). Another interesting property of this method is the fact that the interacting covariate need not be measured or even known to arouse suspicion of an interaction effect, perhaps leading to novel biological insight.
The goal of any genetic prioritization method is to form a subset of genetic variants enriched for the effect of interest. Variance prioritization results in substantial reduction in multiple hypothesis correction and hence increases power by using the information embedded within the variance per genotype, which we showed is independent from subsequent interaction testing. In this case, a frequentist approach was adopted since prior knowledge on the nature of the contribution of genetic interactions to complex traits is currently scarce. Our simulations have shown that the optimal Levene's test Pvalue threshold for variance prioritization will depend upon the allelic frequency, sample size, association of the covariate with the quantitative trait, and expected strength of the interaction. In most cases, all of these parameters can be directly observed with the exception of the interaction strength. The optimal threshold can therefore be empirically selected through simulations when the potentially interacting covariate (or SNP) is characterized and the expected interaction effect size defined.
Inequality of variance between genotypes could be the result of factors other than interactions, such as increased variability. Furthermore, the heteroscedacity of variances also depends on an assumption of homogeneity of variances within genotype, which could be violated if variances tend to be proportional to means. This might reflect an inappropriate choice of scale for the phenotype and a variancestabilizing transform like the log or square root might then be indicated. Linkage disequilibrium with variants of large effect sizes could also theoretically lead to inequality of variance between genotypes. Conversely, nonadditive interactions or SNPs with a multiplicity of interactions may not result in deviation of Levene's test. Failure to prioritize “interacting” SNPs or prioritization of “noninteracting” SNPs can theoretically lead to decreased statistical power but should not lead to inflated type I error since the first step is a screening test followed by confirmatory analysis. While care was taken to simulate models with plausible effect sizes (i.e. explaining less than 1% of total variance), the scenarios described in this report undoubtedly represent simplified situations. Nevertheless, these simulations suggest that variance prioritization has increased power to detect genetic interactions of very modest strength in situations where conventional methods perform poorly.
Genetic effects on variance in an empirical dataset led to the identification of three novel interactions. Given the strength of the observed interactions, they would have also been identified with traditional methods. The significance of these results resides in (1) the existence of SNPs with inequality of variance between genotypes at genomewide significance and (2) the presence of biologically plausible interactions between these SNPs and environmental factors. The first interaction involved the leptin receptor SNP rs12753193 in the prediction of CRP. Leptin is a protein synthesized and secreted by the adipose tissue whose main function is to regulate appetite trough interaction with the leptin receptor in the hypothalamus [16]. As such, leptin levels are directly proportional to BMI and patients deficient in either leptin or leptin receptor are characterized by marked obesity [17]. Interestingly, LEPR is also expressed by the liver, the source of circulating plasma CRP, and physiological concentrations of leptin stimulate expression of CRP in primary hepatocytes [18]. This might explain why only a small fraction of the inequality of variance between the LEPR SNP rs12753193 genotypes was due to the interaction with BMI: the statistical interaction may only capture a small fraction of the underlying biological interaction between leptin levels and its receptor.
The LEPR SNP rs12753193 is known to be associated with CRP concentration at a genomewide level [19]. The minor allele of rs12753193 (i.e. “G”) is associated with lower levels of CRP and was shown in this report to be associated with a greater effect of BMI on CRP concentrations. This confirms the soundness of using marginal effects of SNPs to select a subset for interaction testing. No other SNP with genomewide significance for CRP showed evidence of inequality of variance between genotypes (excluding SNPs in linkage disequilibrium with rs12753193).
The second interaction observed involved the ICAM1 SNP rs1799969 in the prediction of sICAM1 levels. The latter SNP is a coding nonsynonymous SNP (G241R) known to be associated with sICAM1 levels [20]–[22] and was shown in this report to interact with smoking status. Smoking is a very strong determinant of sICAM1 levels (explaining up to 13% of sICAM1 total variance) and was associated with a stronger effect of rs1799969 in WGHS. Interestingly, an interaction between another ICAM1 nonsynonymous SNP (rs5498, r^{2} = 0.17 with rs1799969) and smoking was recently described with regards to susceptibility to coronary artery disease [23].
The third interaction identified involved the PNPLA3 SNP rs738409 in the prediction of sICAM1. PNPLA3 encodes a protein of unknown function that belongs to the patatinlike phospholipase family. The methionine allele of the missense PNPLA3 SNP rs738409 (Ile148Met) has recently been associated with increased hepatic fat levels, hepatic inflammation and plasma levels of liver enzymes (traits often linked to insulin resistance and obesity) [24]. Interestingly, that same variant has been associated with insulin secretion in response to oral glucose tolerance test [25]. Levels of the inflammatory marker sICAM1 are known to be correlated with insulin resistance and obesity [26]. Consistent with this literature, the methionine allele was associated with sICAM1 itself (P = 2.2×10^{−7}; Beta = 4.9 ng/mL per minor allele) while being associated with a greater effect of BMI on sICAM1 levels.
In this report, we demonstrated that the presence of genegene or geneenvironment interactions can result in changes in variance of a quantitative trait per genotype and that these changes can be exploited to prioritize SNPs for interaction testing. We used this property of quantitative trait interactions on a real dataset and discovered three new, biologically relevant, interactions. Variance prioritization presents two advantages over other methods. First, it does not depend on the marginal effect of SNPs on the means. Second, it can identify SNPs likely to be involved in an interaction even when the interacting genetic or environmental covariate(s) is not measured. For instance, had BMI or smoking status not been assessed in the WGHS, the LEPR SNP rs12753193, the ICAM1 SNP rs1799969 and the PNPLA3 SNP rs738409 would still have been selected as candidates for interaction effects, perhaps leading to novel biological hypotheses.
Supporting Information
Zdroje
1. CordellHJ
2002 Epistasis: what it means, what it doesn't mean, and statistical methods to detect it in humans. Hum Mol Genet 11 2463 2468
2. NelsonMR
KardiaSL
FerrellRE
SingCF
2001 A combinatorial partitioning method to identify multilocus genotypic partitions that predict quantitative trait variation. Genome Res 11 458 470
3. RitchieMD
HahnLW
MooreJH
2003 Power of multifactor dimensionality reduction for detecting genegene interactions in the presence of genotyping error, missing data, phenocopy, and genetic heterogeneity. Genet Epidemiol 24 150 157
4. RitchieMD
WhiteBC
ParkerJS
HahnLW
MooreJH
2003 Optimization of neural network architecture using genetic programming improves detection and modeling of genegene interactions in studies of human diseases. BMC Bioinformatics 4 28
5. BochdanovitsZ
SondervanD
PerillousS
van BeijsterveldtT
BoomsmaD
2008 Genomewide prediction of functional genegene interactions inferred from patterns of genetic differentiation in mice and men. PLoS ONE 3 e1593 doi:10.1371/journal.pone.0001593
6. MillsteinJ
ContiDV
GillilandFD
GaudermanWJ
2006 A testing framework for identifying susceptibility genes in the presence of epistasis. Am J Hum Genet 78 15 27
7. MarchiniJ
DonnellyP
CardonLR
2005 Genomewide strategies for detecting multiple loci that influence complex diseases. Nat Genet 37 413 417
8. EvansDM
MarchiniJ
MorrisAP
CardonLR
2006 Twostage twolocus models in genomewide association. PLoS Genet 2 e157 doi:10.1371/journal.pgen.0020157
9. PuttW
PalmenJ
NicaudV
TregouetDA
TahriDaizadehN
2004 Variation in USF1 shows haplotype effects, gene : gene and gene : environment associations with glucose and lipid parameters in the European Atherosclerosis Research Study II. Hum Mol Genet 13 1587 1597
10. ChoYM
RitchieMD
MooreJH
ParkJY
LeeKU
2004 Multifactordimensionality reduction shows a twolocus interaction associated with Type 2 diabetes mellitus. Diabetologia 47 549 554
11. MurcrayCE
LewingerJP
GaudermanWJ
2009 Geneenvironment interaction in genomewide association studies. Am J Epidemiol 169 219 226
12. LeveneH
1960 Contributions to Probability and Statistics: Essays in Honor of Harold Hotelling; Olkin et al., editor: Stanford University Press
13. NIST/SEMATECH eHandbook of Statistical Methods
14. RidkerPM
ChasmanDI
ZeeRY
ParkerA
RoseL
2008 Rationale, design, and methodology of the Women's Genome Health Study: a genomewide association study of more than 25,000 initially healthy american women. Clin Chem 54 249 255
15. RidkerPM
CookNR
LeeIM
GordonD
GazianoJM
2005 A randomized trial of lowdose aspirin in the primary prevention of cardiovascular disease in women. N Engl J Med 352 1293 1304
16. LouisGW
MyersMGJr
2007 The role of leptin in the regulation of neuroendocrine function and CNS development. Rev Endocr Metab Disord 8 85 94
17. FarooqiS
O'RahillyS
2006 Genetics of obesity in humans. Endocr Rev 27 710 718
18. ChenK
LiF
LiJ
CaiH
StromS
2006 Induction of leptin resistance through direct interaction of Creactive protein with leptin. Nat Med 12 425 432
19. RidkerPM
PareG
ParkerA
ZeeRY
DanikJS
2008 Loci related to metabolicsyndrome pathways including LEPR,HNF1A, IL6R, and GCKR associate with plasma Creactive protein: the Women's Genome Health Study. Am J Hum Genet 82 1185 1192
20. PareG
ChasmanDI
KelloggM
ZeeRY
RifaiN
2008 Novel association of ABO histoblood group antigen with soluble ICAM1: results of a genomewide association study of 6,578 women. PLoS Genet 4 e1000118 doi:10.1371/journal.pgen.1000118
21. PonthieuxA
LambertD
HerbethB
DroeschS
PfisterM
2003 Association between Gly241Arg ICAM1 gene polymorphism and serum sICAM1 concentration in the Stanislas cohort. Eur J Hum Genet 11 679 686
22. PuthothuB
KruegerM
BernhardtM
HeinzmannA
2006 ICAM1 aminoacid variant K469E is associated with paediatric bronchial asthma and elevated sICAM1 levels. Genes Immun 7 322 326
23. SareckaHujarB
ZakI
KrauzeJ
2009 Interactions between rs5498 polymorphism in the ICAM1 gene and traditional risk factors influence susceptibility to coronary artery disease. Clin Exp Med 9 117 124
24. RomeoS
KozlitinaJ
XingC
PertsemlidisA
CoxD
2008 Genetic variation in PNPLA3 confers susceptibility to nonalcoholic fatty liver disease. Nat Genet 40 1461 1465
25. JohanssonLE
LindbladU
LarssonCA
RastamL
RidderstraleM
2008 Polymorphisms in the adiponutrin gene are associated with increased insulin secretion and obesity. Eur J Endocrinol 159 577 583
26. SongY
MansonJE
TinkerL
RifaiN
CookNR
2007 Circulating levels of endothelial adhesion molecules and risk of diabetes in an ethnically diverse cohort of women. Diabetes 56 1898 1904
Štítky
Genetika Reprodukční medicínaČlánek vyšel v časopise
PLOS Genetics
2010 Číslo 6
Nejčtenější v tomto čísle
Tomuto tématu se dále věnují…
 Copy Number Variation and Transposable Elements Feature in Recent, Ongoing Adaptation at the Locus
 On the Use of Variance per Genotype as a Tool to Identify Quantitative Trait Interaction Effects: A Report from the Women's Genome Health Study
 Amplification of a Cytochrome P450 Gene Is Associated with Resistance to Neonicotinoid Insecticides in the Aphid
 Use of GenomeWide Expression Data to Mine the “Gray Zone” of GWA Studies Leads to Novel Candidate Obesity Genes