Pregnancy is a valuable model to study the association between DNA methylation and several cardiometabolic traits, due to its direct potential to influence mother’s and child’s health. Epigenetics in Pregnancy (EPIPREG) is a population-based sample with the aim to study associations between DNA-methylation in pregnancy and cardiometabolic traits in South Asian and European pregnant women and their offspring. This cohort profile paper aims to present our sample with genetic and epigenetic data and invite researchers with similar cohorts to collaborative projects, such as replication of ours or their results and meta-analysis. In EPIPREG we have quantified epigenome-wide DNA methylation in maternal peripheral blood leukocytes in gestational week 28±1 in Europeans (n = 312) and South Asians (n = 168) that participated in the population-based cohort STORK Groruddalen, in Norway. DNA methylation was measured with Infinium MethylationEPIC BeadChip (850k sites), with technical validation of four CpG sites using bisulphite pyrosequencing in a subset (n = 30). The sample is well characterized with few missing data on e.g. genotype, universal screening for gestational diabetes, objectively measured physical activity, bioelectrical impedance, anthropometrics, biochemical measurements, and a biobank with maternal serum and plasma, urine, placenta tissue. In the offspring, we have repeated ultrasounds during pregnancy, cord blood, and anthropometrics up to 4 years of age. We have quantified DNA methylation in peripheral blood leukocytes in nearly all eligible women from the STORK Groruddalen study, to minimize the risk of selection bias. Genetic principal components distinctly separated Europeans and South Asian women, which fully corresponded with the self-reported ethnicity. Technical validation of 4 CpG sites from the methylation bead chip showed good agreement with bisulfite pyrosequencing. We plan to study associations between DNA methylation and cardiometabolic traits and outcomes.

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.

Although there are some epigenome-wide association studies (EWAS) of insulin resistance, most of them did not replicate their findings and are focused in populations of European ancestry limiting the generalizability. In EPIPREG (294 Europeans and 162 South Asians), we conducted an EWAS of insulin resistance in maternal peripheral blood leukocytes, with replication in Born in Bradford (n=879; 430 Europeans and 449 South Asians), MENA (n=320) and Botnia (n=56) cohorts. In EPIPREG, we identified six CpG sites inversely associated with insulin resistance across ancestry, whereof five were replicated in independent cohorts (cg02988288, cg19693031, and cg26974062 in TXNIP, cg06690548 in SLC7A11, cg04861640 in ZSCAN26). From methylation quantitative trait loci analysis in EPIPREG, we identified gene variants related to all five replicated cross-ancestry CpG sites, which were associated with several cardiometabolic phenotypes. Mediation analyses suggested that the gene variants regulate insulin resistance through DNA methylation. To conclude, our cross-ancestry EWAS identified five CpG sites related with lower insulin resistance.