Modeling Genetic and Environmental Factors in Biological Systems Using Structural Equation Modeling: An Application to Energy Balance

To improve our understanding of the role(s) that genes and environmental factors play in a complex disease, we need statistical approaches that model multiple factors simultaneously in a hierarchical manner that aims to reflect the underlying biological system(s). We present an approach that models genes as latent constructs, defined by multiple variants (single nucleotide polymorphisms, SNPs) within each gene, using the multivariate statistical framework of structural equation modeling (SEM) to model multiple, putative genetic and environmental factors involved in energy imbalance (dasiaobesitypsila) using subjects from a colon polyp case-control study. We found that modeling constructs for the leptin receptor (LEPR) gene (defined by SNPs rs1137100, rs1137101, rs1805096, rs6588147) and the fat mass-and-obesity-associated (FTO) gene (defined by SNPs rs9939609, rs1421085, rs8044769) together with demographic (age, race, gender), physical activity, diet and sleep variables increased the strength of the association (beta_std=-0.13 plusmn 0.06; p=0.03) between the FTO and obesity constructs compared to that observed in a reduced model with only the FTO and LEPR constructs and demographic variables (beta_std=-0.05 plusmn 0.03; p=0.08). Several indirect paths, including an association between the LEPR and physical activity constructs (beta_std=-0.15 plusmn 0.04; p=0.01), were found. Interestingly, removing FTO revealed a marginal association between the LEPR and obesity constructs (beta_std=0.24 plusmn 0.14; p=0.09), which was not present when FTO was in the model. These results illustrate the importance of modeling multiple relevant genes and other factors in the same model, which is a major strength of this approach. Moreover, our latent gene construct approach exploits the correlation structure between SNPs while capturing overall effects of variation in that gene, which will enable better utilization of candidate gene and genome- -wide SNP array data.