University Of Tasmania
141388 - in epigenomic studies.pdf (1.07 MB)
Download file

In epigenomic studies, including cell-type adjustments in regression models can introduce multicollinearity, resulting in apparent reversal of direction of association

Download (1.07 MB)
journal contribution
posted on 2023-05-20, 18:41 authored by Barton, SJ, Phillip MeltonPhillip Melton, Titcombe, P, Murray, R, Rauschert, S, Lillycrop, KA, Huang, R-C, Holbrook, JD, Godfrey, KM

Background: Association studies of epigenome-wide DNA methylation and disease can inform biological mechanisms. DNA methylation is often measured in peripheral blood, with heterogeneous cell types with different methylation profiles. Influences such as adiposity-associated inflammation can change cell-type proportions, altering measured blood methylation levels. To determine whether associations between loci-specific methylation and outcomes result from cellular heterogeneity, many studies adjust for estimated blood cell proportions, but high correlations between methylation and cell-type proportions could violate the statistical assumption of no multicollinearity. We examined these assumptions in a population-based study.

Methods: CDKN2A promoter CpG methylation was measured in peripheral blood from 812 adolescents aged 17 years (Western Australian Pregnancy Cohort Study). Loge adolescent BMI was used as the outcome in a regression analysis with DNA methylation as predictor, adjusting for age/sex. Further regression analyses additionally adjusted for estimated cell-type proportions using the reference-based Houseman method, and simulations modeled the effects of varying levels of correlation between cell proportions and methylation. Correlations between estimated cell proportions and CpG methylation from Illumina 450K were measured.

Results: Lower DNA methylation was associated with higher BMI when cell-type adjustment was not included; for CpG4, β = −0.004 logeBMI/%methylation (95% CI −0.0065, −0.001; p = 0.003). The direction of association reversed when adjustment for six cell types was made; for CpG4, β = 0.004 logeBMI/%methylation (−0.0002, 0.0089; p = 0.06). Correlations between CpG methylation and cell-type proportions were high, and variance inflation factors (VIFs) were extremely high (25 to 113.7). Granulocyte count was correlated with BMI, and removing granulocytes from the regression model reduced all VIFs to <3.1, with persistence of a positive association between methylation and BMI [CpG4 β = 0.004 logeBMI/%methylation (−0.0002, 0.0088; p = 0.06)]. Simulations supported major effects of multicollinearity on regression results.

Conclusions: Where cell types are highly correlated with other covariates in regression models, the statistical assumption of no multicollinearity may be violated. This can result in reversal of direction of association, particularly when examining associations with phenotypes related to inflammation, as CpG methylation may associate with changes in cell-type proportions. Removing predictors with high correlations from regression models may remove the multicollinearity. However, this might hinder biological interpretability.


Publication title

Frontiers in Genetics



Article number









Menzies Institute for Medical Research


Frontiers Research Foundation

Place of publication


Rights statement

Copyright © 2019 Barton, Melton, Titcombe, Murray, Rauschert, Lillycrop, Huang, Holbrook and Godfrey. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms

Repository Status

  • Open

Socio-economic Objectives

Expanding knowledge in the biological sciences; Expanding knowledge in the mathematical sciences