Skip to main content

Epigenetic approach in obesity: DNA methylation in a prepubertal population which underwent a lifestyle modification

Abstract

Background

Metabolically healthy obesity (MHO) is a considerably controversial concept as it is considered a transitory condition towards the development of different pathologies (type 2 diabetes, insulin resistance, or cardiovascular disease). MHO is closely related to lifestyle and environmental factors. Epigenetics has become an essential biological tool to analyze the link between obesity and metabolic status. The aim of this study was to determine whether MHO status is conditioned by the DNA methylation (DNAm) of several genes related to lipid metabolism (lipoprotein lipase, retinoid X receptor alpha, liver X receptor, stearoyl-CoA desaturase, sterol regulatory element binding factor 1), and inflammation (LEP) in peripheral blood mononuclear cells (PBMCs) from 131 prepubertal subjects with MHO phenotype after lifestyle modifications with personalized Mediterranean diet (MedDiet) combined with a physical activity (PA) program.

Results

The DNAm of all studied genes were significantly modified in the population after 12 months of lifestyle modifications (MedDiet and PA). In addition, associations were found between the DNAm studies and BMI, homeostatic model assessment of insulin resistance, monounsaturated fatty acid and polyunsaturated fatty acid, moderate-vigorous PA, fat mass, and adherence to MedDiet.

Conclusions

It was found that DNAm of genes related to lipid metabolism and inflammation are also present in childhood and that this methylation profile can be modified by interventions based on MedDiet and PA.

Introduction

Obesity is a global health problem that has long-term health repercussions and involves a chronic state of inflammation and increased risk of developing cardiovascular disease, type 2 diabetes, and some types of cancer [1,2,3]. The World Health Organization (WHO) defines obesity as the abnormal or excessive accumulation of fat that can be harmful to health.

In 2016, the WHO estimated that 41 million children under the age of 5 were overweight or obese (http://www.who.int/dietphysicalactivity/childhood/es/), making obesity the most prevalent nutritional disorder in childhood.

It is known that obesity is a multifactorial entity influenced by environmental (including dietary habits and physical activity (PA)), psychosocial, and neuroendocrine status as well as genetics and epigenetics [4]. The increase in the prevalence of this pathology has led to an increase in the number of studies on its treatment and prevention, including some that study the influence of epigenetics in response to diet and predisposition to weight gain [5]. Childhood and adolescence are critical periods for the development of obesity. Many studies have demonstrated that obesity developed during infancy is associated with a high risk of being overweight in adult life [6]. However, obesity is phenotypically a very heterogeneous pathology and, concretely, includes the metabolically healthy obese (MHO) phenotype. This phenotype describes individuals with obesity who do not meet the criteria of the metabolic syndrome [7]. They may have an epigenetic profile different than that of other phenotypes associated with obesity.

Epigenetics is defined as inheritable and reversible phenomena that affect gene expression without altering the DNA sequence through DNA and histone protein modifications, including DNA methylation (DNAm, the most stable epigenetic modification), covalent histone modifications, chromatin folding, and the regulatory noncoding miRNAs. Epigenetic changes can modify gene transcription by altering the accessibility of gene transcription machinery [8, 9]. These changes are influenced by environmental factors such as nutrition, inflammation, hypoxia, PA, sex, and age, with these epigenetic changes acting as a regulator of the environment-gene interaction and the way the genome responds to lifestyle changes [10, 11].

Obesity is a multifactorial disease involving interactions between an individual’s genetic makeup and an unhealthy lifestyle (unhealthy food intake and sedentarism). It is known that epigenetic modifications affect the expression of different metabolic genes, including lipid metabolism and inflammation genes involved in obesity pathogenesis. Identifying subjects who present with changes in DNAm could help to predict their susceptibility to gain or lose weight, thus helping to prevent obesity via the implementation of new therapeutic approaches [12,13,14,15].

A variety of genes have been proposed as candidates for methylation studies. Variations in the methylation profiles of various genes related to the metabolic pathway in obesity are associated with BMI, adiposity, and waist circumference (WC) [16].

The objective of this study was to analyze the methylation profile of genes involved in lipid metabolism after a personalized lifestyle modification based on the Mediterranean diet (MedDiet) and physical activity (PA), during 4-month and 12-month periods, in a prepubertal population of children with metabolically healthy obesity.

Results

Anthropometric and clinical changes

The characteristics of the population at baseline, 4 months, and 12 months after the intervention and by gender and for all participants are summarized in Table 1.

Table 1 Anthropometric (A) and analytical (B) parameters at baseline, after 4 months, and 12 months of intervention in the total study population and by gender (mean ± SD)

Regarding anthropometric (Table 1A) variables, participants’ weight and height showed significant differences at 12 months compared to baseline, both in the general population and by sex. After 12 months of intervention, weight and height showed a significant increase with respect to baseline conditions in the total population (+ 4.0 ± 0.9 kg and + 5.9 ± 1.5 cm). When data were analyzed according to sex, significant increases were also observed (boys + 2.5 ± 1.2 kg and + 4.8 ± 1.1 cm; girls + 5.7 ± 1.0 kg and + 7.0 ± 0.3 cm). Of the total population, 91.8% of participants gained weight. In the female group, 46% gained weight compared to 44% in the male group. On the other hand, only 8.2% of the total population lost weight, with boys losing more weight than girls (− 3.0 ± 3.0 kg vs − 0.3 ± 0.2 kg). All participants gained height (100%), with girls gaining slightly more height than boys. No significant changes in BMI were observed in the total population (0.0 ± 0.4 kg/m2). When studying the genders separately, BMI decreased slightly in boys (− 0.4 ± 0.2 kg/m2) but increased in girls (+ 0.5 ± 0.4 kg/m2). Also, WC increased significantly in the total population (1.4 ± 0.7 cm), although this effect was only observed in male sex (1.0 ± 0.6 cm). Finally, SBP increased in the population as a whole (3.0 ± 0.8 mmHg), in this case due to a significant increase in SBP in girls (6.0 ± 0.6 mmHg).

In regard to the analytical parameters (Table 1B), after 12 months of the intervention, the entire population showed a statistically significant increase in fasting blood glucose (within normal range). Insulin levels and homeostatic model assessment of insulin resistance (HOMA-IR) decreased after 12 months of intervention. The lipid profile also underwent several changes, although they were not statistically significant. Both interleukin 6 (IL-6) and tumor necrosis factor alpha (TNFa) levels increased significantly after the intervention.

Analysis of Mediterranean diet adherence in the population showed that it improved by almost 2 points in the total population after 4 and 12 months of intervention (baseline 7.2 ± 1.7; 4 months 8.8 ± 2.0, p (B vs 4 m) < 0.0001; 12 months: 8.4 ± 2.0, p (B vs 12 m) = 0.01).

Lifestyle modification

In Table 2, energy intake is summarized. The energy intake significantly decreased after the intervention. In addition, all components of the diet decreased significantly except for the consumption of fiber and vitamin D. Data on energy expenditure reflect significant variations in both levels of sedentarism as well as in light, moderate, and vigorous PA (Table 3) after the intervention. Sedentary time per day increased in the total population in a significant manner (p < 0.0001) while moderate and vigorous exercise time both increased as well. Lastly, changes in participants’ body composition (Table 4) showed that lean mass and total mass were significantly increased and total fat decreased in the entire population after 12 months of intervention.

Table 2 Energy and food intake at baseline, after 4 months and 12 months of intervention in the total study population and by gender (mean ± SD). All, n = 131 participants (n = 70 boys and n = 61 girls; p = 0.65)
Table 3 Energy expenditure measured using a GENEActiv Actigraph GT3X+ accelerometer at baseline, after 4 months and 12 months of intervention in the total study population and by gender (mean ± SD)
Table 4 Body composition as measured by a bone densitometer-DXA at baseline, after 4 months and 12 months of intervention in the total study population and by gender (mean ± SD)

DNAm profiles

Results on the methylation status of the studied genes after 4 and 12 months are summarized in Table 5. Methylation of the lipoprotein lipase (LPL) gene increased in the female group after 4 months (p = 0.02) and was maintained at 12 months (p = 0.04). Also, methylation of leptin (LEP) gene increased in the female group (p = 0.02) after 12 months of intervention. No significant modifications were observed in retinoid X receptor alpha (RXRa) methylation. The methylation level of the liver X receptor (LXR) gene decreased slightly at 4 months in the total population (p = 0.03) and in the female group (p = 0.04), but returned to baseline levels at 12 months. Regarding stearoyl-CoA desaturase (SCD) gene methylation, a slight decrease was found at 4 months that only continued at 12 months in the male group (p = 0.05). For the sterol regulatory element binding factor 1 (SREBF1) gene, methylation levels increased slightly in the total population after 4 months but returned to baseline levels after 12 months. In the male group, methylation levels increased at both times (p = 0.03 and p = 0.04, respectively) while, in girls, methylation levels decreased at 12 (p = 0.03) months.

Table 5 Variation in percentage of methylation of genes studied after 4 months and 1 year of a lifestyle intervention

Correlations between DNAm levels and anthropometric parameters, clinical parameters, energy intake and consumption, and body composition

The methylation levels of every gene along with anthropometric variables, analytical parameters, food intake, PA, and body composition all showed statistically significant correlations.

There were associations between anthropometric variables, analytical parameters, energy and food intake, energy expenditure, body composition, and adherence to MedDiet in the MHO prepubertal population. At baseline, the methylation of the LPL gene correlated negatively with abdominal visceral fat (r = − 0.24, p = 0.03) and positively with adherence to the MedDiet (r = 0.22, p = 0.04) and moderate (r = 0.27, p = 0.02) and vigorous (r = 0.26, p = 0.02) PA. Methylation level of the LXR gene correlated positively with weight (r = 0.33, p = 0.03), abdominal visceral fat mass (r = 0.39, p = 0.01), and abdominal total mass (r = 0.36, p = 0.02) in males. SCD gene methylation correlated positively with energy (r = 0.27, p = 0.02), total fat (r = 0.32, p = 0.01), monounsaturated fatty acid (MUFA) (r = 0.38, p = 0.001), polyunsaturated fatty acid (PUFA) (r = 0.25, p = 0.04), cholesterol (r = 0.27, p = 0.02), and vitamin D (r = 0.35, p = 0.003) intake. SREBF1 gene methylation correlated negatively with moderate (r = − 0.42, p = 0.001) and vigorous (r = 0.34, p = 0.002) PA. Methylation of the LEP gene correlated negatively with sedentarism (r = − 0.33, p = 0.03) and light PA (r = − 0.25, p = 0.02) and positively with HDL-c levels (r = 0.21, p = 0.04). In girls, RXRa methylation correlated negatively with energy intake (r = − 0.45, p = 0.01), fat intake (r = − 0.58, p < 0.001), and saturated fatty acid (SFA) and MUFA intake (r = − 0.45, p = 0.01; r = − 0.58, p < 0.001, respectively) and positively with moderate (r = 0.36, p = 0.03) and vigorous (r = 0.40, p = 0.02) PA. In all populations, this gene correlated positively with adherence to the MedDiet (r = 0.30, p = 0.01).

After 4 months of intervention, the methylation of LPL gene correlated negatively to BMI (r = − 0.22, p = 0.04), WC (r = − 0.23, p = 0.03), and abdominal viscera fat (r = − 0.23, p = 0.04) and positively with adherence to the MedDiet (r = 0.33, p = 0.04). The methylation level of the LXR gene correlated positively with BMI (r = 0.48, p = 0.001), WC (r = 0.38, p = 0.01), abdominal visceral fat mass (r = 0.41, p = 0.01), abdominal total mass (r = 0.39, p = 0.01), and fiber intake (r = 0.40, p = 0.03) in males. The SCD gene methylation correlated negatively with PUFA intake (r = − 0.26, p = 0.04), Furthermore, in females, it correlated negatively with glucose levels (r = − 0.33, p = 0.04) and HOMA-IR (r = 0.34, p = 0.05) and positively with adherence to MedDiet (r = 0.24, p = 0.03). SREBF1 gene methylation correlated negatively with light PA (r = − 0.31, p = 0.01) and MUFA intake (r = − 0.30, p = 0.02) and positively with sedentarism (r = 0.28, p = 0.02). Methylation of the LEP gene correlated negatively with LDL-c levels (r = −0.31, p = 0.01). In females, RXRa methylation correlated negatively with SBP (r = − 0.38, p = 0.02) and abdominal viscera fat (r = − 0.24, p = 0.03) and positively with HDL-c levels (r = 0.49, p = 0.01).

Finally, after 12 m of intervention, the methylation of the LPL gene correlated negatively with abdominal visceral fat (r = − 0.24, p = 0.04) and positively with adherence to the MedDiet (r = 0.22, p = 0.04) and moderate (r = 0.27, p = 0.02) and vigorous (r = 0.26, p = 0.02) PA. In males, RXRa methylation correlated negatively with SBP (r = − 0.38, p = 0.01) and DBP (r = − 0.35, p = 0.02) and positively with HDL-c levels (r = 0.38, p = 0.02) and with adherence to the MedDiet (r = 0.38, p = 0.02). In males, methylation levels of the LXR gene correlated negatively with cholesterol intake (r = − 0.42, p = 0.01). SCD gene methylation correlated positively with insulin levels (r = 0.33, p = 0.005) and HOMA-IR (r = 0.27, p = 0.02). Methylation of the LEP gene correlated negatively with light (r = − 0.38, p = 0.02) and vigorous (r = − 0.36, p = 0.07) PA and positively with insulin levels (r = 0.36, p = 0.03) in females.

To study the strength of the association observed in the correlation analyses, we performed linear regression analyses. We observed that the nutrition, PA, and analytical parameters could explain DNAm levels of our target genes (Tables 6), a regression that was corrected for BMI. Baseline DNAm levels of the LPL gene correlated negatively with energy intake (p = 0.03) and hsPCR levels (p = 0.01) and positively with carbohydrate intake (p = 0.02), protein (p = 0.04), MUFA (p = 0.01), PUFA (p = 0.02), MedDiet adherence (p = 0.03), and resistin levels (p = 0.03). In addition, DNAm levels of the LEP gene at baseline correlated negatively with BMI (p = 0.02) and vigorous PA (p = 0.05) and positively with moderate PA (p = 0.05). After 12 months of intervention, this gene correlated negatively with light PA (p = 0.02). On the other hand, DNAm levels of the SREBF gene correlated negatively at baseline with moderate PA (p = 0.01) and positively with glucose levels (p = 0.002), HOMA-IR (p = 0.01), and hsPCR levels (p = 0.01). DNAm levels of the RXRa gene correlated negatively with hsPCR levels (p = 0.03) at baseline and negatively with both light and vigorous PA (p = 0.01, p = 0.03, respectively) after 12 months. DNAm levels of the LXR gene correlated negatively with glucose levels (p = 0.02) and positively with HOMA-IR index (p = 0.03) at 12 months. Lastly, DNAm levels of the SCD gene correlated negatively with HOMA-IR index (p = 0.03), high-density lipoprotein cholesterol (HDL-c) (p = 0.01), resistin (p = 0.01), and IL-6 (p = 0.001) levels and positively with insulin (p = 0.01) and TNFa (p = 0.02) levels.

Table 6 Linear regression analysis with DNA methylation (DNAm) as the dependent variable and nutrition components, physical activity, or analytical parameters as independent variables and corrected for BMI

Discussion

Personalized lifestyle modification is a great tool for analyzing the pattern in DNA modification from peripheral blood mononuclear cells (PBMCs). Epigenetic marks can change due to metabolic status, dietary intake or PA [17]. In this work, we studied these variations in DNAm in prepubertal population with obesity that was metabolically healthy. After 12 months of intervention, our population showed an improvement in insulin resistance and hyperinsulinemia compared to baseline conditions, components that play a fundamental role in the development of the obesity-related cardiometabolic risk factors [18]. This improvement in clinical parameters was produced as a result of a decrease in intake of different nutrients and total fat and an increase in moderate-vigorous PA. Both lean mass and total mass decreased after 12 months of intervention in the entire population.

Epigenetic marks are used in studies as predictors of the response to weight loss programs. As in our intervention, the DNA has different methylation patterns at the beginning of the study in comparison to the end of follow-up, data supported by Moleres et al. [19]. Some studies have also proposed epigenetics as a predictor of the predisposition towards developing obesity in both childhood and adulthood [16]. Herrera et al. [20] demonstrated the importance of an obesogenic environment. Early environmental influences induce epigenetic variation, which permanently affects metabolism and risk of chronic disease. Evidence suggests that the establishment of the epigenome can be affected by environmental factors during critical developmental periods [21]. There are potential interactions between lifestyle and epigenetic mechanisms that mediate the expression of genes associated with increased BMI and adiposity, as is suggested by the effect of diet on the methylation of genes involved in lipid metabolism and inflammation such as leptin [22]. On the other hand, dietary intake of methyl groups (choline, methionine, genistein, and folate) during critical periods of developmental alters DNA promoter and histone methylation, thereby resulting in lifelong changes in gene expression and alteration of the epigenome that inclines the individual towards obesity in adulthood [23]. Some cross-sectional studies have reported a significant association between obesity or adiposity status, and DNA methylation [24]. These studies demonstrate that some regions of different genes studied were hypomethylated in obese children and located in the gene body region, and revealed a unique cluster of obese individuals that was differentiated from the normal-weight children [25]. Some of these genes are implicated in lipid and glucose metabolism, differential body size, and body composition in children [26].

The genes studied are involved in lipid metabolism and inflammation. The implication of these genes in the obesity process is still unknown. In this context, we observed that the methylation pattern is modified after different periods of follow-up (4 and 12 months) with lifestyle modification (MedDiet and PA) which, according to other studies, showed the relationship between obesity and epigenetic marks. In our study, BMI, an indicator of body adiposity, influenced DNAm. After 12 months of intervention, we verified that the DNAm is located in genes involved in lipid metabolism and inflammatory pathways, results in concordance with Wahl et al. [27]. As far as is known, LPL mediates the used of blood triglycerides. It is a hydrolytic enzyme located in the capillary endothelium that catalyzes the cleavage of triglycerides into fatty acids. Therefore, LPL activity is critical for plasma triglyceride (TG) clearance and tissue uptake of fatty acids [28]. We found a significant correlation between the LPL DNAm levels and adherence to the MedDiet, moderate and vigorous PA, BMI, WC, and abdominal visceral fat at all points of follow-up, in line with previous results [29, 30]. LPL is expressed in different tissues, including adipose tissue and foam cells, and its expression is regulated by LXR/RXR heterodimers. LPL activity is closely linked with RXR [31]. RXRa is a member of the steroid/thyroid hormone superfamily of nuclear receptors that function as transcription factors binding to promoter regions of genes. These receptors exert their action by binding, as homodimers or as heterodimers, to specific sequences in the promoters of the target genes and the regulation of their transcription. In our population, we found DNAm of RXRa had a significant correlation with energy, fat, SFA and MUFA intake, moderate and vigorous PA, and adherence to the MedDiet at baseline. During the follow-up, the DNAm of RXRa showed a significant correlation with HDL-c levels and fat mass at the end of intervention [32]. It is known that RXR heterodimerizes with subfamily 1 nuclear receptors, including LXR [33]. LXRs are members of the nuclear hormone receptor superfamily of ligand-activated transcription factors that regulate cholesterol and lipid metabolism. Their regulation of gene expression in response to changes in intracellular cholesterol levels link lipid metabolism, inflammation, and immune cell function [34].

LXR can be activated by increased intracellular cholesterol concentrations [35]. According to these data, in our MHO prepubertal population, the DNAm of LXR significantly correlated with weight, BMI, WC, cholesterol and fiber intake, and body composition during the intervention’s follow-up. In addition, we found a significant correlation with the DNAm of SCD. SCD is a key enzyme in the conversion of polyunsaturated fatty acids (PUFAs) to monounsaturated fatty acids (MUFAs), and it has been shown to be overexpressed in adipose tissue in obesity. The DNAm of SCD correlated in our population with energy, total fat, PUFA, MUFA, and cholesterol and vitamin D intake at baseline, in accordance with other data [29]. Recent studies suggest that decreased functional SCD promotes lipid oxidation over storage, decreased lipogenesis, increased beta-oxidation, and increased glucose utilization [36]. In this sense, during the follow-up on our participants, the DNAm of SCD correlated with glycemic metabolism parameters (glucose, insulin levels, and HOMA-IR) [37] and adherence to the MedDiet.

Within lipid metabolism, SRE sites are fundamental in the regulation of LPL gene expression by sterols. On the other hand, sterols contribute to the control of triglyceride metabolism via binding of SREBF1 to the LPL regulatory sequences [38]. SREBPs are loop-helix leucine zipper transcription factors that regulate the synthesis and cellular uptake of cholesterol and fatty acids [39]. According to the relation of this gene to fat metabolism, we found correlations with the DNAm of LPL, SCD, and RXRa. In addition, DNAm of SREBF1 significantly correlates with PA and MUFA intake.

Finally, leptin is a hormone secreted by adipocytes whose function is to inhibit intake and stimulate energy expenditure, thus allowing for body weight maintenance. Leptin resistance is related to diet-induced obesity (high-fat diets), which is the main cause of obesity in humans. The decrease in leptin sensitivity leads to deregulation of satiety, which increases intake and affects nutrient absorption, metabolism, insulin sensitivity, and energy balance [40]. In our population, DNAm of leptin correlated with sedentarism and HDL-c levels at baseline. It is known that leptin can regulate cholesterol-ester metabolism and the atherosclerotic process associated with obesity [41]. A significant association has also been found between the DNAm of leptin and insulin and LDL-c levels as well as vigorous PA at the end of follow-up. In accordance with these data, some studies about dietary intervention identify DNAm as biomarkers of a weight loss response [42]. These genomic loci are located in proximity to obesity-associated genes such as LEP. These epigenetic signatures may potentially be used in precision medicine as baseline biomarkers to predict the effectiveness of a weight loss intervention [5].

Limitations

The limitations of this study include the low weight loss and lack of big changes in the BMI of our participants. However, the period of follow-up of the lifestyle modification of our prepubertal population, which took place over a year, is a strength of this study, since it provides an evaluation of the impact of nutrition and PA over the epigenetic modifications of different genes involved in the inflammation process and lipid metabolism. In addition, MedDiet adherence questionnaire used in this study is validated for Mediterranean populations, but not specifically for a prepubertal population. Another limitation is that there was no randomization or control group (normal-weight subjects) to safely attribute changes, though the results and the longitudinal design are strengths.

Conclusion

This study shows the results of the DNAm of specific genes related to lipid metabolism and inflammation in metabolically healthy obesity and in a prepubertal population. The DNAm of the genes studied are closely related to different anthropometric and analytical parameters, components of food intake, type of PA, and body composition as well as between those genes. In addition, these results also suggest that the influence of methylation profile is present even in childhood and that this methylation profile can be modified by interventions based on MedDiet and PA.

In terms of public health, a healthy lifestyle in the entire population from childhood to old age is essential so that the epigenetic modifications that occur in human DNA do not lead to the development of different molecular process and pathologies associated with obesity, such as insulin resistance, dyslipidemia, inflammatory processes, hypertension, alterations of the vascular structure, and cardiovascular diseases.

Material and methods

Study population

In this study, a prepubescent population aged 4–9 years, of both sexes, was included. Inclusion criteria were boys (from 4 years to testicular volume < 3 ml) and girls (from 4 years to Tanner S2, breast bud elevation) with obesity (≥ 95 percentile) [43]. Participants were metabolically healthy, meaning they presented ≤ 1 of the following criteria: WC and blood pressure ≥ 90 percentile, triglycerides > 150 mg/dl, HDL-c < 40 mg/dl, or impaired fasting glucose [44]. Participants who met more than one criterion of metabolic syndrome, did not meet the age criterion at the start of the study, were diabetic, or had any metabolic pathology were excluded from the study.

The recruitment period was from November 2016 to May 2017 and took place in different schools in the city of Málaga (Andalusia, Spain). Recruitment was carried out via visits to pre-school and elementary school classrooms. Once the possible participants were selected, families were contacted to inform them about the study and invite them to participate. The selected subjects and their parents/guardians were summoned to the Civil Hospital of Malaga in order to inform them about the design and objectives of the study as well as to request their written informed consent for the voluntary participation of their children in the study.

Study design

The subjects who wished to participate in the study had an initial visit with a nurse. Weight, height, BMI, WC, hip circumference, waist-hip ratio, and blood pressure measurements were taken. Next, blood samples were taken to analyze analytical parameters (triglyceride levels, HDL-c, total cholesterol, low-density lipoprotein cholesterol (LDL-c), glycosylated hemoglobin A1c (HbA1c) (%), creatinine, uric acid, microalbuminuria, insulin, HOMA-IR, and glucose after 12 h of fasting) to confirm that they met the inclusion criteria as well as to determine baseline levels before the intervention. The samples were tested in the Clinical Analysis Laboratory of the Regional Hospital of Malaga. Lastly, lifestyle questionnaires (nutritional questionnaires on the frequency of food consumption, a PA questionnaire, and questionnaires on adherence to the Mediterranean diet) were completed.

Once the parents/guardians granted permission for the inclusion of their child in the study and it was verified that the child met the inclusion criteria, they had an interview with a nutritionist at the beginning of the study. The nutritionist informed the parents/guardians and the children about the characteristics and recommendations of the Mediterranean diet and the healthy lifestyle they should follow during the study period. The participants and their parents/guardians participated in an interview with the nutritionist in which they completed lifestyle questionnaires (nutritional questionnaires about the frequency of food consumption, PA questionnaires, and questionnaires about adherence to the Mediterranean diet) in order to know their lifestyle prior to the intervention. Dietary records completed for three non-consecutive days and a food frequency questionnaire for the previous year (number of times/day, number of days/week, number of days/14 days, number of days/month, rarely, or never) were completed at the baseline visit and 4 and 12 months after the intervention [45] (https://www.healthychildren.org/Spanish/healthy-living/nutrition/Paginas/Energy-In-Recommended-Food-Drink-Amounts-for-Children.aspx). A food frequency questionnaire divided into nine sections showing food groups was used. These sections are diary, eggs, meat and fish, fruits, vegetables, legumes and cereals, oils and fats, pastries, miscellaneous, and beverages. The FFQ includes 9 possible responses for frequency of consumption ranging from never/rarely to > 6 servings/day. In addition, a self-reported record of food questionnaire was used for three consecutive days, 2 days of daily, and another day of the weekend.

The Mediterranean diet recommended to the study subjects included extra virgin olive oil and nuts. The recommended calorie intake was 1500 kcal/day (https://www.healthychildren.org/Spanish/healthy-living/nutrition/Paginas/Energy-In-Recommended-Food-Drink-Amounts-for-Children.aspx) with the following macronutrient distribution [46]: 30% fats (5–8% saturated fatty acids, 15–18% monounsaturated fatty acids, 5–8% polyunsaturated fatty acids, and less than 300 mg of cholesterol per day), 50% carbohydrates (less than 10% simple sugars, 40% complex sugars, and low glycemic index), and 20% protein. Adherence to the Mediterranean diet was determined according to the instructions of Trichopoulou et al. [47].

They were also informed that they should perform regular PA with two trained supervisors with a degree in the physical activity and sports sciences field. International PA guidelines recommend that children exercise daily (http://www.health.gov/paguidelines/chapterthree). To do so, participants were given access to sports programs. We offered the possibility of attending an exercise program daily from Monday to Friday (5 sessions/week for 120 min/session). Participants were required to attend a minimum of 3 sessions per week. Attendance was recorded at the beginning of each PA session. The sessions consisted of aerobic exercises, resistance training, and exercises to improve flexibility and balance. Sessions were held for the duration of the study. In addition, participants kept a PA record using a GENEActiv Actigraph GT3X+ accelerometer to measure energy expenditure. It was used for one full week at the baseline visit and again at 4 and 12 months.

In addition, parents/guardians were provided with the telephone number of those in charge of the study so that they could ask questions about nutrition and PA that may have arisen before the next visit.

All patients participating in the study gave their informed consent, and protocols were approved by the institutional ethics committee (Comité de Ética de la Investigación Provincial de Málaga, belonging to the Andalusian Health Service).

Epigenetic study

For methylation analysis, PBMCs were used to extract genomic DNA using TRIzol™ Reagent. The epigenetic profiling of genomic DNA was analyzed in 131 MHOCh using the methylation assays (Qiagen). Bisulfite conversion of DNA samples was performed using the EZ DNAm kit (Zymo Research, Orange, CA, USA) on 500 ng of DNA. A PCR of methylated DNA was carried out to amplify the genes of interest using PyroMark PCR Kit (Qiagen) and Primers PyroMark CpG Assays (Qiagen). Pyrosequencing was carried out using PyroMark Gold Q96 and primers PyroMark CpG Assay (Qiagen). The proportion of methylation (%) for each subject at each CpG site was computed by first subtracting the background signal intensity of negative controls from both the methylated and unmethylated signals then dividing the ratio of the methylated signal intensity by the sum of both methylated and unmethylated signals.

Table 7 shows 6 amplicons for the validated CpGs. These gene candidates are lipoprotein lipase (LPL), retinoid X receptor alpha (RXRa), liver X receptor (LXR), stearoyl-CoA desaturase (SCD), sterol regulatory element binding protein (SREBF1), and leptin (LEP).

Table 7 Gene characteristics of differentially methylated CpGs

Statistical analysis

Quantitative variables were expressed as mean ± standard deviation (SD), and qualitative variables were expressed as percentages. Student’s t test was used to compare quantitative variables whereas the chi-square test was used for qualitative variables in order to contrast variables measured within each group at different time periods. Bivariate correlations were determined using the Pearson correlation coefficient analysis. To calculate the sample size (calculated using SISA, Simple Interactive Statistical Analysis), we relied on other clinical studies that have shown the metabolic benefits of weight loss in the pediatric population with obesity [48]. Assuming a 95% confidence level (α error of 0.5%), a statistical power of 80%, and a 5% loss rate, a sample of 110 MHO subjects was required. To account for non-participation and possible loss to follow-up, recruitment of at least 130 participants was planned. Statistical analysis was performed using SPSS for Windows, version 22.0 (IBM Corporation INC. Somers, NY, USA).

Abbreviations

WHO:

World Health Organization

MHO:

Metabolically healthy obese

DNAm:

Deoxyribonucleic acid methylation

miRNAs:

Micro ribonucleic acid

MedDiet:

Mediterranean diet

PA:

Physical activity

BMI:

Body mass index

WC:

Waist circumference

SBP/DBP:

Systolic blood pressure/diastolic blood pressure

HbA1c:

Glycosylated hemoglobin A1c

HOMA-IR:

Homeostatic model assessment of insulin resistance

LDL-c:

Low-density lipoprotein cholesterol

HDL-c:

High-density lipoprotein cholesterol

TG:

Triglycerides

IL-6:

Interleukin 6

TNFa:

Tumor necrosis factor alpha

LPL:

Lipoprotein lipase

RXRa:

Retinoid X receptor alpha

LXR:

Liver X receptor

SCD:

Stearoyl-CoA desaturase

SREBF1:

Sterol regulatory element binding factor 1

LEP:

Leptin

SFA:

Saturated fatty acid

MUFA:

Monounsaturated fatty acid

PUFA:

Polyunsaturated fatty acid

PBMCs:

Peripheral blood mononuclear cells

References

  1. Assunção SNF, Boa Sorte NCA, Alves CAD, Mendes PSA, Alves CRB, Silva LR. Glucose and insulin resistance alteration in asymptomatic obese children and adolescents. J Pediatr. 2017;(16):30322–9.

  2. Genoni G, Menegon V, Secco GG, Sonzini M, Martelli M, Castagno M, Ricotti R, Monzani A, Aronici M, Grossini E, Di Mario F, Bona G, Bellone S, Prodama F. Insulin resistance, serumuric acid and metabolic syndrome are linked to cardiovascular dysfunction in pediatric obesity. Int J Cardiol. 2017;(17):33656–2.

  3. Williams EP, Mesidor M, Winters K, Dubbert PM, Wyatt SB. Overweight and obesity: prevalence, consequences, and causes of a growing public health problem. Curr Obes Rep. 2015;4:363–70.

    Article  PubMed  Google Scholar 

  4. Pigeyre M, Yazdi FT, Kaur Y, Meyre D. Recent progress in genetics, and epigenetics Metagenomics unveils the pathophysiology of human obesity. Clin Sci. 2016;130:943–86.

    Article  CAS  Google Scholar 

  5. Aronica L, Levine AJ, Brennan K, Mi J, Gardner C, Haile RW, Hitchins MP. A systematic review of studies of DNA methylation in the context of a weight loss intervention. Epigenomics. 2017;9(5):769–87.

    Article  CAS  PubMed  Google Scholar 

  6. Druet C, et al. Prediction of childhood obesity by infancy weight gain: an individual-level meta-analysis. Paediatr. Perinat. Epidemiol. 2012;26:19–26.

    Google Scholar 

  7. Teixeira TF, Alves RD, Moreira AP, Peluzio Mdo C. Main characteristics of metabolically obese normal weight and metabolically healthy obese phenotypes. Nutr Rev. 2015;73:175–90.

    Article  PubMed  Google Scholar 

  8. Kaushik P. and. Anderson J.T. Obesity: epigenetic aspects. Biomol Concepts 2016; 7(3):145-55.

  9. Goldberg AD, Allis CD, Bernstein E. Epigenetics: a landscape takes shape. Cell. 2007;128(4):635–8.

    Article  CAS  PubMed  Google Scholar 

  10. Seisenberger S, Peat JR, Hore TA, Santos F, Dean W, Reik W. Reprogramming DNA methylation in the mammalian life cycle: building and breaking epigenetic barriers. Philos Trans R Soc Lond Ser B Biol Sci. 2013;368:20110330.

    Article  CAS  Google Scholar 

  11. Mansego ML, Milagro FI, Zulet MA, Moreno-Aliaga MJ, Martínez JA. Differential DNA methylation in relation to age and health risks of obesity. Int J Mol Sci. 2015;16:16816–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Campion J, Milagro F, Martinez JA. Epigenetics and obesity. Prog Mol Biol Transl Sci. 2010;10:291–347.

    Article  Google Scholar 

  13. Martinez J, Milagro F, Claycombe K, Schalinske K. Epigenetics in adipose tissue, obesity, weight loss, and diabetes. Adv Nutr. 2014;5:71–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. McAllister E, Dhurandhar N, Keith S, Aronne L, Barger J, Baskin M, Benca RM, Biggio J, Boggiano MM, Eisenmann JC, Elobeid M, Fontaine KR, Gluckman P, Hanlon EC, Katzmarzyk P, Pietrobelli A, Redden DT, Ruden DM, Wang C, Waterland RA, Wright SM, Allison DB. Ten putative contributors to the obesity epidemic. Crc Cr Rev Food Sci. 2009;49:868–913.

    Article  Google Scholar 

  15. Milagro F, Alfredo Martinez J. Epigenetics of obesity and weight loss. Endocrinol Nutr. 2013;60:12–4.

    Article  PubMed  Google Scholar 

  16. Van Dijk SJ, Molloy PL, Varinli H, Morrison JL, Muhlhausler BS. Epigenetics and human obesity. Int J Obes. 2014;39(1):85–97.

    Article  CAS  Google Scholar 

  17. Chamberlain JA, Dugué PA, Bassett JK, Hodge AM, Brinkman MT, Joo JE, Jung CH, Makalic E, Schmidt DF, Hopper JL, Buchanan DD, English DR, Southey MC, Giles GG, Milne RL. Dietary intake of one-carbon metabolism nutrients and DNA methylation in peripheral blood. Am J Clin Nutr. 2018;108(3):611–21.

    Article  PubMed  Google Scholar 

  18. Damanhoury S, Newton AS, Rashid M, Hartling L, Byrne JLS, Ball GDC. Defining metabolically healthy obesity in children: a scoping review. Obes Rev. 2018;19:1476–91.

    Article  CAS  PubMed  Google Scholar 

  19. Moleres A, Campión J, Milagro FI, Marcos A, Campoy C, Garagorri JM, Gómez-Martínez S, Martínez JA, Azcona-Sanjulián MC, Martí A. Differential DNA methylation patterns between high and low responders to a weight loss intervention in overweight or obese adolescents: the EVASYON study. FASEB J. 2013;27(6):2504–12.

    Article  CAS  PubMed  Google Scholar 

  20. Herrera BM, Keildson S, Lindgren CM. Genetics and epigenetics of obesity. Maturitas. 2011;69(1):41–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Waterland RA, Jirtle RL. Early nutrition, epigenetic changes at transposons and imprinted genes, and enhanced susceptibility to adult chronic diseases. Nutrition. 2004;20(1):63–8.

    Article  CAS  PubMed  Google Scholar 

  22. Milagro FI, Campion J, Garcia-Diaz DF, Goyenechea E, Paternain L, Martinez JA. High fat diet-induced obesity modifies the methylation pattern of leptin promoter in rats. J Physiol Biochem. 2009;65(1):1–9.

    Article  CAS  PubMed  Google Scholar 

  23. Zeisel SH. Epigenetic mechanisms for nutrition determinants of later health outcomes. Am J Clin Nutr. 2009;89:1488S–93S.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Samblas M, Milagro FI, Martínez A. DNA methylation markers in obesity, metabolic syndrome, and weight loss. Epigenetics. 2019;14(5):421–44.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Rhee JK, Lee JH, Yang HK, Kim TM, Yoon KH. DNA Methylation profiles of blood cells are distinct between early-onset obese and control individuals. Genomics Inform. 2017;15(1):28–37.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Rzehak P, Covic M, Saffery R, Reischl E, Wahl S, Grote V, Weber M, Xhonneux A, Langhendries JP, Ferre N, Closa-Monasterolo R, Escribano J, Verduci E, Riva E, Socha P, Gruszfeld D, Koletzko B. DNA methylation and body composition in preschool children: epigenome- wide-analysis in the European Childhood Obesity Project (ChOP)-study. Sci Rep. 2017;7:14349.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Wahl S, Drong A, Lehne B, et al. Epigenome-wide association study of body mass index, and the adverse outcomes of adiposity. Nature. 2017;541(7635):81–6.

    Article  CAS  PubMed  Google Scholar 

  28. Kersten S. Physiological regulation of lipoprotein lipase. Biochim Biophys Acta. 2014;1841(7):919–33.

    Article  CAS  PubMed  Google Scholar 

  29. Castellano-Castillo D, Moreno-Indias I, Fernández-García JC, Alcaide-Torres J, Moreno-Santos I, Ocaña L, Gluckman E, Tinahones F, Queipo-Ortuño MI, Cardona F. Adipose tissue LPL methylation is associated with triglyceride concentrations in the metabolic syndrome. Clin Chem. 2018;64(1):210–8.

    Article  CAS  PubMed  Google Scholar 

  30. Perfilyev A, Dahlman I, Gillberg L, Rosqvist F, Iggman D, Volkov P, Nilsson E, Risérus U, Ling C. Impact of polyunsaturated and saturated fat overfeeding on the DNA-methylation pattern in human adipose tissue: a randomized controlled trial. Am J Clin Nutr. 2017;105(4):991–1000.

    Article  CAS  PubMed  Google Scholar 

  31. Davies PJ, Berry SA, Shipley GL, Eckel RH, Hennuyer N, Crombie DL, Ogilvie KM, Peinado-Onsurbe J, Fievet C, Leibowitz MD, Heyman RA, Auwerx J. Metabolic effects of rexinoids: tissue-specific regulation of lipoprotein lipase activity. Mol Pharmacol. 2001;59(2):170–6.

    Article  CAS  PubMed  Google Scholar 

  32. Godfrey KM, Sheppard A, Gluckman PD, Lillycrop KA, Burdge GC, McLean C, Rodford J, Slater-Jefferies JL, Garratt E, Crozier SR, Emerald BS, Gale CR, Inskip HM, Cooper C, Hanson MA. Epigenetic gene promoter methylation at birth is associated with child's later adiposity. Diabetes. 2011;60(5):1528–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Germain P, Chambon P, Eichele G, Evans RM, Lazar MA, Leid M, De Lera AR, Lotan R, Mangelsdorf DJ, Gronemeyer H. International Union of Pharmacology. LXIII. Retinoid X receptors. Pharmacol Rev. 2006;58(4):760–72.

    Article  PubMed  CAS  Google Scholar 

  34. Schulman IG. Liver X receptors link lipid metabolism and inflammation. FEBS Lett. 2017;591(19):2978–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Janowski BA, Grogan MJ, Jones SA, Wisely GB, Kliewer SA, Corey EJ, Mangelsdorf DJ. Structural requirements of ligands for the oxysterol liver X receptors LXRalpha and LXRbeta. Proc Natl Acad Sci U S A. 1999;96(1):266–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Poudyal H, Brown L. Stearoyl-CoA desaturase a vital checkpoint in the development and progression of obesity. Endocr Metab Immune Disord Drug Targets. 2011;11(3):217–31.

    Article  CAS  PubMed  Google Scholar 

  37. Yamakawa T, Ogihara K, Utsunomiya H, Muraoka T, Kadonosono K, Terauchi Y. Colestimide improves glycemic control via hepatic glucose production in db/db mice. Endocr J. 2014;61(5):425–36.

    Article  CAS  PubMed  Google Scholar 

  38. Schoonjans K, Gelman L, Haby C, Briggs M, Auwerx J. Induction of LPL gene expression by sterols is mediated by a sterol regulatory element and is independent of the presence of multiple E boxes. J Mol Biol. 2000;304(3):323–34.

    Article  CAS  PubMed  Google Scholar 

  39. Madison BB. Srebp2: a master regulator of sterol and fatty acid synthesis. J Lipid Res. 2016;57(3):333–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Crujeiras AB, Carreira MC, Cabia B, Andrade S, Amil M, Casanueva FF. Leptin resistance in obesity: an epigenetic landscape. Life Sci. 2015;140:57–63.

    Article  CAS  PubMed  Google Scholar 

  41. O'Rourke L, Yeaman SJ, Shepherd PR. Insulin and leptin acutely regulate cholesterol ester metabolism in macrophages by novel signaling pathways. Diabetes. 2001;50(5):955–61.

    Article  CAS  PubMed  Google Scholar 

  42. Milagro FI, Campión J, Cordero P, Goyenechea E, Gómez-Uriz AM, Abete I, Zulet MA, Martínez JA. A dual epigenomic approach for the search of obesity biomarkers: DNA methylation in relation to diet-induced weight loss. FASEB J. 2011;25(4):1378–89.

    Article  CAS  PubMed  Google Scholar 

  43. Cole TJ, Bellizzi MC, Flegal KM, Dietz WH. Establishing a standard definition for child overweight and obesity worldwide: international survey. BMJ. 2000;320(7244):1240–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Cook S, Kavey RE. Dyslipidemia and pediatric obesity. Pediatr Clin N Am. 2011;58(6):1363–73.

    Article  Google Scholar 

  45. Martin-Moreno JM, Weevil L. Assessment of dietary intake at the population level through individual questionnaires: methodological shadows and lights. Rev Esp Public Health. 2007;81:507–18.

    Google Scholar 

  46. Davis C, Bryan J, Hodgson J, Murphy K. Definition of the Mediterranean Diet: a literature review. Nutrients. 2015;7(11):9139–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Trichopoulou A, Costacou T, Bamia C, Trichopoulos D. Adherent to a Mediterranean diet and survival in a Greek population. N Engl J Med. 2003;348(26):2599–608.

    Article  PubMed  Google Scholar 

  48. Reinehr T, Elfers C, Lass N, Roth CL. Irisin and its relation to insulin resistance and puberty in obese children: a longitudinal analysis. J Clin Endocrinol Metab. 2015;100(5):2123–30.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to thank the “Miraflores de los Angeles” school for allowing us to use its sports facilities for PA and the rest of schools in Malaga city for facilitating contact with the parents of MHOCh participants in this study and to Claire Alexandra Conrad for her help with the final English language version.

Funding

This work was supported by grants from the Instituto de Salud Carlos III, cofinanced by the Fondo Europeo de Desarrollo Regional-FEDER (“Centros de Investigación En Red” (CIBER, CB06/03/0018)). M Rosa Bernal-Lopez was supported by “Miguel Servet Type I” program (CP15/00028) from the ISCIII-Madrid (Spain), cofinanced by the Fondo Europeo de Desarrollo Regional-FEDER.

Author information

Authors and Affiliations

Authors

Contributions

CG-E analyzed the methylation experiments and written and revised the manuscript. VB and AV-C performed the methylation experiments. MIR-M performed the data collection and clinical interviews and analyzed the anthropometric and analytical parameters. AV-P performed the nutrition program and analyzed the nutrition questionnaires. JB-P and ARR-G performed the physical activity program and analyzed the accelerometry results. AR-M performed the data collection and analyzed the anthropometric and analytical parameters. RG-H supervised the clinical intervention and revised the manuscript. MRB-L recruited the population, designed the study, analyzed the results, and wrote and reviewed manuscript. The authors read and approved the final manuscript.

Corresponding authors

Correspondence to Ricardo Gomez-Huelgas or M. Rosa Bernal-Lopez.

Ethics declarations

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gallardo-Escribano, C., Buonaiuto, V., Ruiz-Moreno, M.I. et al. Epigenetic approach in obesity: DNA methylation in a prepubertal population which underwent a lifestyle modification. Clin Epigenet 12, 144 (2020). https://doi.org/10.1186/s13148-020-00935-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13148-020-00935-0

Keywords