Distinct Associations of BMI and Fatty Acids With DNA Methylation in Fasting and Postprandial States in Men

We have previously shown that blood global DNA methylation (DNAm) differs between postprandial state (PS) and fasting state (FS) and is associated with BMI and polyunsaturated fatty acid (PUFA) (negatively and positively, respectively) in 12 metabolically healthy adult Mexican men (AMM cohort) equally distributed among conventional BMI classes. Here, we detailed those associations at CpG dinucleotide level by exploiting the Infinium methylation EPIC array (Illumina). We sought differentially methylated CpG (dmCpG) that were (1) associated with BMI (BMI-dmCpG) and/or fatty acids (FA) (FA-dmCpG) in FS or PS and (2) different across FS and PS within a BMI class. BMI-dmCpG and FA-dmCpG were more numerous in FS compared to PS and largely prandial state-specific. For saturated and monounsaturated FA, dmCpG overlap was higher across than within the respective saturation group. Several BMI- and FA-dmCpG mapped to genes involved in metabolic disease and in some cases matched published experimental data sets. Notably, SETDB1 and MTHFS promoter dmCpG could explain the previously observed associations between global DNAm, PUFA content, and BMI in FS. Surprisingly, overlap between BMI-dmCpG and FA-dmCpG was limited and the respective dmCpG were differentially distributed across functional genomic elements. BMI-dmCpG showed the highest overlap with dmCpG of the saturated FA palmitate, monounsaturated C20:1 and PUFA C20:2. Of these, selected promoter BMI-dmCpG showed opposite associations with palmitate compared to C20:1 and C20:2. As for the comparison between FS and PS within BMI classes, dmCpG were strikingly more abundant and variably methylated in overweight relative to normoweight or obese subjects (∼70–139-fold, respectively). Overweight-associated dmCpG-hosting genes were significantly enriched in targets for E47, SREBP1, and RREB1 transcription factors, which are known players in obesity and lipid homeostasis, but none overlapped with BMI-dmCpG. We show for the first time that the association of BMI and FA with methylation of disease-related genes is distinct in FS and PS and that limited overlap exists between BMI- and FA-dmCpG within and across prandial states. Our study also identifies a transcriptional regulation circuitry in overweight that might contribute to adaptation to that condition or to transition to obesity. Further work is necessary to define the pathophysiological implications of these findings.