Abstract

Obesity is increasing worldwide in prepubertal children, reducing the age of onset of associated comorbidities, including type 2 diabetes. Sulfur-containing amino acids, methionine, cysteine, and their derivatives play important roles in the transmethylation and transsulfuration pathways. Dysregulation of these pathways leads to alterations in the cellular methylation patterns and an imbalanced redox state. Therefore, we tested the hypothesis that one-carbon metabolism is already dysregulated in prepubertal children with obesity. Peripheral blood was collected from 64 children, and the plasma metabolites from transmethylation and transsulfuration pathways were quantified by HPLC. The cohort was stratified by BMI z-scores and HOMA-IR indices into healthy lean (HL), healthy obese (HO), and unhealthy obese (UHO). Fasting insulin levels were higher in the HO group compared to the HL, while the UHO had the highest. All groups presented normal fasting glycemia. Furthermore, high-density lipoprotein (HDL) was lower while triglycerides and lactate levels were higher in the UHO compared to HO subjects. S-adenosylhomocysteine (SAH) and total homocysteine levels were increased in the HO group compared to HL. Additionally, glutathione metabolism was also altered. Free cystine and oxidized glutathione (GSSG) were increased in the HO as compared to HL subjects. Importantly, the adipocyte secretory function was already compromised at this young age. Elevated circulating leptin and decreased adiponectin levels were observed in the UHO as compared to the HO subjects. Some of these alterations were concomitant with alterations in the DNA methylation patterns in the obese group, independent of the impaired insulin levels. In conclusion, our study informs on novel and important metabolic alterations in the transmethylation and the transsulfuration pathways in the early stages of obesity. Moreover, the altered secretory function of the adipocyte very early in life may be relevant in identifying early metabolic markers of disease that may inform on the increased risk for specific future comorbidities in this population.

1. Introduction

Obesity is a rapidly growing epidemic that is contributing to the significant increase in metabolic diseases worldwide. It is characterized by excess adipose tissue expansion and is associated with low-grade inflammation and metabolic dysfunction [1]. The continuous release of proinflammatory cytokines [1] and adipokines (e.g., leptin) by dysregulated adipose tissue may contribute to the obesity-associated inflammation [2]. It is thought that chronic low-grade inflammation induces chronic oxidative stress, and that both contribute to the obesity-related insulin resistance (IR) and type 2 diabetes (T2D) development [3]. Due to the drastic increase in early childhood obesity, the journey to T2D development is starting earlier in life. This in turn increases the risk for other severe health complications over the lifespan, such as hypertension, cardiovascular diseases (CVD), retinopathy, and neuropathy, that appear to increase as the age of T2D onset decreases [46]. Further, there are several important differences in the pathophysiology of obesity-associated comorbidities in adults compared to children, including early β-cell decline and time to T2D treatment failure, as well as the lack of appropriate pharmacological medications approved for earlier ages, and longer duration of the disease [710].

Early alterations in the redox and methylation status that are associated with obesity may play a significant role in the early onset of metabolic disturbances in children with obesity. The thiol group plays an important role in biological systems [11]. It appears in the sulfur-containing amino acids methionine and cysteine, and their derivatives, such as glutathione (GSH) and other low molecular weight intermediates in the transmethylation and transsulfuration pathways [11, 12], also known as aminothiols. Thiols are responsible for scavenging reactive oxygen species (ROS) and maintaining redox homeostasis [13]. In particular, cysteine is primarily responsible for maintaining the redox state in plasma [14], while GSH maintains intracellular redox homeostasis, acting directly or indirectly through enzymatic activity [13]. Recent studies postulate dysfunction in the redox homeostasis in obese children [15, 16]. Lechuga-Sancho et al. [15] have identified an altered oxidative status in erythrocytes from obesity-associated insulin resistant children, even before those changes occurred in plasma. Besides, Zalewska et al. [16] reported alterations in the saliva redox status followed by higher oxidative damage in obese when compared to overweight children.

Transmethylation and transsulfuration intermediates are also critically important for methylation of DNA, proteins, and lipids [14] with methionine-derived S-adenosylmethionine (SAM), being the primary methyl group donor [17]. Imbalance in the transmethylation and transsulfuration pathways is linked with obesity-related inflammation [18]. It has also been shown that high levels of circulating homocysteine, resulting from the S-adenosylhomocysteine (SAH) degradation, are linked to an increased oxidation status in circulation [19]. Moreover, the hyperhomocysteinemia resulting from the imbalance of transsulfuration and transmethylation metabolites has been linked to further risk of obesity-associated CVD, such as atherosclerosis [18, 19]. Interestingly, a study conducted in mice suggested that high levels of SAH in the circulation could be involved with alterations at the epigenetic levels by inhibiting the DNA methyltransferase enzymes. The same study also indicated a possible relation between high levels of SAH and endothelial dysfunction [20].

One-carbon metabolism pathways, including those described above, have been implicated in important metabolic processes that include redox defenses and epigenetic alterations, which are both altered in obesity [11]. However, it is not known how soon this can happen in life and whether these pathways, if becoming altered in prepubertal children with obesity, can facilitate the early onset of obesity-related comorbidities.

Therefore, the present study is mainly aimed at testing the hypothesis that one-carbon metabolism perturbation is already present in the early stages of obesity development, in prepubertal children. Therefore, transsulfuration and transmethylation metabolite levels were quantified and related with their systemic oxidative stress, genomic methylation status, and inflammatory marker levels in children of normal weight or with overweight/obesity.

2. Material and Methods

2.1. Study Cohort

A group of 64 prepubertal children (5-9 years old, Tables 1 and 2 from Results) were recruited after approval of the study by the Institutional Review Board (IRB) (protocol number 206164) at the University of Arkansas for Medical Science and following the guidelines of Declaration of Helsinki (1964). This clinical study was registered at ClinicalTrials.gov (NCT03323294). The inclusion criteria were age 5-9 years at the date of the visit (i.e., 5- <10 years), and the exclusion criteria were the presence of known chronic illnesses/disorders that might affect study outcome measures, such as type 1 diabetes mellitus, neurologic, developmental, endocrine, hepatic, autoimmune, cardiac, and renal disorders; use of any medication could affect study outcomes, e.g., antipsychotics, thyroid hormone replacement therapy, inhalation/oral steroids, insulin, anabolic drugs and stimulants, or being classified as underweight based on the CDC growth charts (http://www.cdc.gov/growthcharts).

Anthropometric variables were collected for all study participants and sex (male/female), age (years), weight (kg), height (cm), and waist circumference (cm) were included (Tables 1 and 2). For data analyses, children with an years and 6 months, but <10 years, were considered 10 years old (Tables 1 and 2). The weight was measured using a calibrated Avery Berkel, HL122 Series Platform Scale (Dynamic Scales, Terre Haute, IN, USA) wearing minimal clothing, while height was obtained using a stadiometer (Novel Products, Rockton, IL, USA). The waist circumference was measured as reported previously [21].

Body mass index (BMI) was calculated from body mass and height as kg/m2 and adjusted for age and sex according to the Centers for Diseases Control and Prevention (http://www.cdc.gov/growthcharts). The participants were considered overweight or obese if their age- and sex-adjusted BMI was above the 85th percentile (i.e., ). Although during the stratification present in Statistical Analysis, all the participants with overweight and obesity were included in obese groups. Clinical outcomes such as systolic and diastolic blood pressure (mmHg), as well as heart rate (bpm), were also measured using a digital sphygmomanometer (Tables 1 and 2). The measurement was performed on an arm rested at heart level, and the cuff was placed two fingers above the brachial artery. These measurements were performed at the Pediatric Clinical Research Unit from Arkansas Children’s Hospital using a GE Carescape V100 Dinamap vital sign monitor following the standard procedures for this unit. The instrument is calibrated for children.

Fat-free mass, fat mass, and total body water were also measured using the Tanita Body Composition Analyzer (Model TBF-300A; Tanita Corporation of America, Inc., Arlington Heights, IL, USA).

2.2. Blood Collection and Processing

Fasting venous blood samples were collected in EDTA tubes to isolate peripheral blood mononuclear cells (PBMCs), as previously described [21]. Plasma samples were collected after whole blood centrifugation (1,500 × g for 30 min at 4°C). Thereafter, samples were stored up to 1-2 years at -80°C, until the study was concluded, so that all samples could be measured together to reduce batch effects. Then, the plasma volume was replaced with wash buffer consisting of Ca2+/Mg2+-free PBS supplemented with 2 mM EDTA and 0.1% BSA (Sigma Aldrich, St. Louis, MO). To perform the gradient separation, Histopaque-1077 (Sigma Aldrich) was used. The diluted blood was layered on histopaque and centrifuged at 400 × g for 30 min at room temperature. The white cloudy layer of PBMCs was collected and washed two times with ~20 ml of room temperature wash buffer. PBMCs were counted using a hemocytometer (Bright-Line; Hausser Scientific, Horsham, PA), and 2-5 million PBMCs were pelleted, snap frozen on dry ice, and stored up to two years at -80°C, until the study was concluded [21].

2.3. Biochemical Measures

Fasting insulin concentration was measured in plasma using the Mesoscale Discovery Platform (MSD Multi-Array Assay System, Gaithersburg, MD, USA) according to the manufacturer’s protocol. Fasting plasma glucose was measured using YSI 2900 biochemistry analyzer (YSI Life Sciences, Yellow Springs, OH, USA). The lipid profile was quantified in plasma using a RX Daytona clinical analyzer accordingly to the manufacturer’s instructions (Randox Laboratories-IS Limited, Kearneysville, WV, USA)—nonesterified fatty acids (NEFA: mmol/L), glycerol (μmol/L), high-density lipoprotein (HDL: mmol/L), low-density lipoprotein (LDL: mmol/L), and triglycerides (TGs: mmol/L). Additionally, plasma lactate (mmol/L) and CRP were measured using the same methodology.

Fasting insulin (μU/mL) and glucose concentrations (mmol/L) were used to calculate HOMA-IR and HOMA-β using the following equations:

HOMA-IR was used to determine the insulin sensitivity status for each participant. When , the participant was considered insulin resistant [2224].

2.4. Sample Preparation for Aminothiol Analysis

Plasma was prepared for analysis as previously described by Melnyk et al. [25] in order to determine free reduced and oxidized or total reduced aminothiols. Briefly, to assess the total concentration of aminothiols, 50 μL of a solution containing 1.43 M of sodium borohydride, 66 mM sodium hydroxide, 1.5 μM EDTA, and 10 μL n-amyl alcohol was added to 200 μL of plasma and incubated for 30 min at 40°C. Thereafter, the proteins were precipitated by incubation for 10 min with cold 10% meta-phosphoric acid, the samples were centrifuged for 15 min at 14,000 RPM, and 20 μL of supernatant was measured by HPLC. To assess the free and oxidized aminothiol concentration, an equal volume of 10% meta-phosphoric acid was added to the plasma samples and treated as previously described [25].

2.5. Aminothiols and Oxidative Damage Marker Identification

Total and free aminothiols were separated using a Shimadzu HPLC with a Shimadzu pump model 580 on a 5 μm,  mm i.d. reverse-phase C18 column (MCM, Inc., Tokyo, Japan) with the thermostat at 25°C. An isocratic mobile phase composed of 50 mM sodium phosphate, 1.0 mM of reagent OSA, and 2% acetonitrile (v/v) at pH 2.7 was used. The detection of all compounds was carried out using a Coulochem II EC detector, model 5200A (ESA, Inc.). The identification was carried out using external standards for each compound: methionine, homocysteine, cysteine, cystine, cysteinylglycine, reduced and oxidized glutathione, gamma-glutamylcysteine, 3-nitro-tyrosine, and 3-chloro-tyrosine, as previously described [25].

The percentage of oxidized GSH was obtained using the following equation [14]:

2.6. Plasma Adipokines and Cytokine Quantification

A plasma adipokine and cytokine kit was used to measure leptin, IL-1β, IL-6, IL-8, MCP-1, and TNF-α by multiplexing using a Milliplex® Map Human Adipokine Panel (Millipore®, MA, USA). Adiponectin was also measured using a human Adiponectin ELISA (Millipore®, MA, USA). All procedures were performed according to the manufacturer’s instructions.

2.7. DNA Methylation Profile

DNA methylation was assessed in the isolated PBMCs. The Puregene Blood Kit (Gentra Systems, Inc., Minneapolis, MN, USA) was used to extract the DNA, and it was further bisulfite-converted and purified using an EZ DNA Methylation-Gold kit (Zymo Research, Irvine, CA, USA) according to the manufacturer’s protocol [26]. After bisulfite-conversion, the methylation was determined using the Infinium MethylationEPIC bead chip from Illumina®. The acquired data was followed by a quality control analysis of samples and probes, followed by further normalization using the Bioconductor packages minfi v1.34.0 and watermelon version 1.32.0 [27, 28] in R version 4.0.2 [29, 30]. In order to reduce the bias within-array, the data was normalized combining Noob+BMIQ (β-mixture quantile normalization) in order to improve signal intensities [27, 31]. After normalization, the data was filtered and probes that failed ( value >0.01) were removed. All probes mapped to the X and Y chromosomes were also removed to avoid sex chromosome bias. Finally, cross-reactive probes [32] and probes including known SNPs were also removed, according to Illumina recommendations, before final statistical analysis [33].

2.8. Statistical Analysis

The original statistical power for the present study was computed with a total of 110 children based on a one-factor ANOVA, 80% power and were assumed, and a minimum detectable Cohen’s effect size of 0.33 was used. To test the differences between HL and HO and between HO and UHO, a -test was performed when data fulfilled all the assumptions—the normal distribution was tested by the Shapiro-Wilk test, and the variance homogeneity was tested by the Levene’s test. Otherwise, a Wilcoxon signed-rank test was performed. Results are presented as and median (Q1–Q3) according to the respective test. The correlation between continuous variables was assessed using the Spearman’s Rank-Order correlation, and the coefficient (rho) is shown for each correlation. A value <0.05 was considered statistically significant. These tests were performed using the R version 4.0.2 [29, 30].

For statistical analysis of methylation results, differentially methylated positions (DMPs) were tested among groups using limma v3.44.3 [34] R package, after converting values into values. Covariates such as sex, age, and race were adjusted to the linear model. The values were adjusted by the Benjamini-Hochberg method (false discovery rate [FDR]) [35]. DMPs were considered significative for . Gene set enrichment analyses were performed using the webtool STRING database v11.0b [36], for Gene Ontology (GO), KEGG (Kyoto Encyclopedia of Genes and Genomes), and Reactome pathways. Significant results were defined as .

3. Results

3.1. Physiologic and Biochemical Characterization of Study Population

The physiologic and biochemical characteristics of the study population were stratified according to the BMIz and HOMA-IR, as shown in Tables 1 and 2, respectively. The HO showed a significantly higher BMIz () and waist circumference (), as compared with the HL subjects. In addition, the HO displayed a significantly higher fat-free mass () and fat mass () compared with HL. Interestingly, the HO showed significantly elevated plasma insulin levels () when compared with the HL, despite normal fasting glucose levels and a . The HL presented a better β-cell insulin secretory function (HOMA-β) compared to the HO group (Table 1). Although obesity is normally characterized by dyslipidemia, this was not observed in the HO group, as their lipid profile was similar to that of HL. The systolic blood pressure was elevated in the HO as compared to HL subjects ().

Differences between the HO and UHO groups are presented in Table 2. The data indicate that UHO participants had higher BMIz as compared to the HO. Their higher BMIz was caused by a significant higher fat mass, as well as fat-free mass, that was accompanied by a significantly higher total body water in UHO compared to the HO. Insulin was significantly higher in the UHO group, in comparison to the HO group. This was also accompanied by a significant insulin secretory dysfunction, as represented by the HOMA-β index (). Interestingly, these metabolic defects are already present in this prepubertal cohort of obese children, even in the presence of normal fasting plasma glycemia. The fasting glucose levels are similar among all the groups. While there were no differences in plasma cholesterol levels in the HL vs. HO, the UHO presented dyslipidemia which was characterized by a decrease in HDL cholesterol () and significant increase in triglycerides levels (). The LDL cholesterol levels were slightly increased in the UHO compared to the HO participants, but the difference did not reach statistical significance () in the UHO vs. HO participants. In spite of the young age of this cohort of UHO subjects, they already presented significantly elevated systolic blood pressure, as compared to the HO. Furthermore, the plasma lactate concentration was also significantly higher in the UHO compared to the HO. The lactate concentration was fairly well correlated with high levels of plasma insulin (Spearman’s correlation, , value <0.01) and with HOMA-IR (Spearman’s correlation, , value ≤0.01) (Supplementary Table 1B).

3.2. Transmethylation Metabolites

Transmethylation metabolites, such as methionine, SAH, SAM, adenosine, and homocysteine were quantified in the three groups of participants (Table 3). Of the evaluated metabolites, SAH levels were significantly increased in the HO group compared with the HL subjects. Interestingly, similar results were also observed for homocysteine levels. However, there was no difference in the SAM/SAH ratio. This ratio is frequently used to predict the methylation capacity of the cells. When UHO were compared with HO, no significant differences were found between groups.

3.3. Transsulfuration Metabolites

In parallel to the measurements of plasma transmethylation metabolites, Table 4 shows the levels of plasma metabolites related with the redox state as well as the transsulfuration pathway. When comparing the HL with the HO prepubertal children, several metabolites were significantly altered. The HO children exhibited increased levels of total cysteine (tCysteine; free circulating+protein-bound), while free cysteine (fCysteine) was not different as compared to the HL. Similarly, increased levels of cystine were observed. The cysteine oxidation ratio (fCysteine/cystine) is an important redox buffer responsible for maintaining the plasma redox state. However, the ratio was not different between groups. The homocysteine levels described above play an important role in the transsulfuration pathway, since it acts as an intermediary metabolite between both the transmethylation and the transsulfuration pathways. In fact, homocysteine is the main source of cysteine that is used to synthetize glutathione (GSH). Total reduced GSH (tGSH) and free reduced GSH (fGSH) were also measured but showed no significant differences between groups. On the other hand, the oxidized glutathione (GSSG) was significantly elevated in the HO compared with the HL subjects. The evaluation of tGSH/GSSG and fGSH/GSSG ratios showed a significant reduction in the antioxidant capacity and consequently an increase in oxidative stress in plasma of the HO children. The percentage of oxidized GSH was also higher in HO subjects compared to HL. Additionally, the levels of cysteinylglycine in HO group were also significantly higher, when compared with the HL.

Biomarkers of nitrosative stress, 3-chloro-tyrosine and 3-nitro-tyrosine, showed no significant differences between groups.

3.4. Inflammation Patterns

Oxidative stress is normally accompanied by an increase in the systemic inflammatory status. Therefore, we measured inflammatory cytokines as well as the inflammatory marker, C-reactive protein (CRP) in the plasma (Table 5). While CRP was elevated in HO as compared with the HL, no significant alterations were observed for inflammatory cytokines except for the unexpected finding of reduced TNF alpha levels in UHO as compared to HO subjects. Also, presented in Table 5 are the circulating levels of leptin and adiponectin, which are important cytokines secreted by adipose tissue. Adiponectin was reduced in the UHO as compared to the HO subjects while leptin levels and the leptin/adiponectin ratio were significantly elevated in the UHO as compared to the HO as well as in HO when compared to the HL groups. Increased levels of CRP were correlated with the leptin levels (Spearman’s correlation, , ) (Supplementary Table 1A).

3.5. DNA Methylation Pattern

Since the transmethylation and transsulfuration pathways have a complementary loop as shown in Figure 1, and to further support our described findings of metabolic perturbations with possible effect at the DNA methylation mechanism, we analyzed the DNA methylation profile in PBMCs from a subset of the participants (, , and ). Since no significant differences were observed between the HO and the UHO groups, both groups were merged (OverallObese) to achieve higher statistical power. Therefore, the HL () were compared with the OverallObese subjects () adjusting the model for HOMA-IR, thus enabling us to better isolate the effect of obesity.

From the cytosine-phosphate-guanines (CpGs) analyzed, 4677 were differentially methylated between the two groups () (Supplementary Table 2). Furthermore, 35% of the significant DMPs presented a reduction in the methylation status (hypomethylation) in the OverallObese group, while 65% of the DMPs were hypermethylated in the same group. Moreover, 24 DMPs were selected from the list of 4677 DMPs based on their association with genes that are directly or indirectly related with one-carbon metabolism and consequently associated with oxidative stress and/or methylation processes [3739], as explained before (Table 6).

Some of these DMPs are associated with genes involved in the expression of important enzymes, such as methionine synthase reductase (MTRR gene), methylenetetrahydrofolate dehydrogenase (MTHFD1 gene), methylenetrahydrofolate reductase (MTHFR gene), and glycine-N-methyltransferase (GNMT gene). These enzymes are key in one-carbon metabolism and are responsible for methionine regeneration through the homocysteine conversion and the transmethylation pathway [40]. From these important results, it is possible to predict a downregulation in these enzymes since their genes are hypermethylated (positive value or fold change in Table 6). The alteration in the methylation status of these enzymes could indeed explain the observed increase in SAH and homocysteine plasma levels observed in the HO (Table 3). Moreover, the methylation results also showed alteration in the regulation of genes involved in oxidative stress, which include glutathione peroxidase 1 and 7 (GPX1 and GPX7), oxidation resistance 1 (OXR1), and superoxide dismutase 1 (SOD1). All these genes were hypermethylated in the OverallObese group, suggesting epigenetic-mediated downregulation. Interestingly, the enzymes responsible for DNA methylation, i.e., DNA methyltransferases (transcribed by the gene DNMT1), seem to be downregulated in the OverallObese group as well.

In order to disclose the biological meaning of these results, gene enrichment analysis was performed for the 24 DMP set, using the STRING database [36]. Importantly, the enrichment analysis showed 207 significant biological processes (Supplementary Table 3). Furthermore, 19 of them are terms involved in aminothiol metabolism, such as “response to oxidative stress,” “methylation,” “S-adenosylhomocysteine metabolic process,” and “homocysteine metabolic process” (Figure 2(a)). In addition, 16 molecular functions, including SAM-dependent methyltransferase and peroxidase activity and 10 cellular components, including mitochondrial components, appear to be affected through the methylation pattern (Figures 2(b) and 2(c)). Interestingly, mitochondria seem to be of the most affected cellular component (Figure 2(c)). The KEGG pathway analysis, shown in Figure 3(a), indicates alterations in cysteine, methionine, and glutathione metabolism, as well as “one-carbon pool by folate.” The Reactome pathway enrichment corroborates previous findings by showing alterations in metabolism and their association to epigenetic regulation (Figure 3(b)).

4. Discussion

Overnutrition and poor-quality diets are triggering a severe increase in obesity worldwide, with resulting metabolic disorders starting early in life, in particular during childhood [41]. Our results show that obesity has already caused profound changes in several aspects of metabolism in a cohort of prepubertal children, particularly in the one-carbon metabolism, of particular importance in the folate, transmethylation, and transsulfuration pathways. These studies, evaluating early childhood obesity, can inform on the potential origin of the related comorbidities that start plaguing many, already in young adulthood. We first chose to stratify children with obesity into two groups, the healthy and unhealthy obese, in order to understand the main differences between both conditions, where the unhealthy obese already presented insulin resistance. Our main criteria for this stratification were based on the , although different phenotypes could be used to characterize and differentiate the metabolic state of this population, as reviewed by Phillips [41], who described various valid ways to characterize and distinguish between metabolically healthy and metabolically unhealthy obese children [41]. This separation between obese groups may give some important insights into adulthood comorbidities linked to obesity [41], and our analyses have revealed essential differences in several metabolic processes already emerging in this pediatric cohort. All pediatric subjects in this study presented normoglycemia, while the state of insulin resistance driven by the high circulating insulin concentrations observed in the UHO group likely reflects an attempt to maintain their euglycemia [42]. This is in agreement with other studies indicating that the insulin-resistant state could be present years before any alteration in circulating glucose are detected [43, 44]. Chronic periods of insulin resistance, even in the absence of elevated fasting glucose, may be an important contributing factor to the early blood pressure alterations also observed in our cohort [45, 46]. The increase in insulin production is reflected by the increased HOMA-β observed in the obese unhealthy subjects, reflecting increased β-cell activity. Similar to insulin, lactate levels were also elevated in the insulin-resistant obese group. This phenomenon has already been observed by Hosking et al. [47] in children with insulin resistance, indicating a positive correlation between insulin resistance and lactate levels during childhood and adolescence, even after controlling for BMIz as covariate BMIz. Furthermore, Berhane et al. [48] have shown that plasma lactate was increased in adults during a hyperinsulinemic euglycemic clamp, a method that mimics the hyperinsulinemic state. Besides, they also indicated that high levels of lactate were present even before the insulin resistance was clinically detected [48]. Our data show that children as young as 5 to 10 years of age already show significantly elevated levels of lactate in the circulation, when comparing the HO to the UHO subjects. Moreover, there was a positive correlation between insulin and lactate levels as previously demonstrated [47]. Obesity-related insulin resistance has been involved in the development of other metabolic conditions, including CVD [49, 50] and lipid dysregulation [51, 52]. The insulin-resistant UHO subjects also presented lipid dysregulation and increased systolic blood pressure, risk factors for CVD. Furthermore, chronic low-grade inflammation and increased oxidative stress are characteristics of obesity [15]. The quantification of thiol molecules during prepubertal obesity and insulin resistance has not been well described. Some of these molecules are intermediaries of one-carbon metabolism and are important for the maintenance of redox homeostasis and the methylation capacity of cells by acting on transsulfuration and transmethylation pathways. Our data clearly show important differences in the transmethylation pathway, especially when comparing the HL and HO subjects. The HO showed higher levels of SAH and total homocysteine, similarity to the previously reported findings by Kumar et al. [53]. These molecules have already been implicated in the development of atherosclerosis and CVD [20, 54]. Similar alterations have been identified in diabetic patients with renal dysfunction, and the homocysteine levels were also positively correlated with insulin levels [40]. In addition, Chiang et al. [40] have observed the effect of insulin (1 μM) in HepG2 cell lines and showed that mRNA expression of different enzymes related with one-carbon metabolism, such as MTRR, was reduced during the treatment with insulin. In agreement, the HO subjects in the present study had higher fasting insulin levels that are accompanied with high levels of homocysteine in plasma, compared to their lean counterparts. Importantly, DNA methylation analysis, from PBMCs isolated from the prepubertal children with obesity, showed hypermethylation of different CpGs that are localized in the vicinity of important metabolic genes, including PEMT, MTRR, and MTHFR, suggesting that an epigenetic effect could be the cause for this downregulation. Gene Ontology enrichment analysis also showed alteration in different key biological processes, molecular functions, and cellular components related with one-carbon metabolism, as shown in Figure 2. Even though the SAM/SAH ratio was not altered in our cohort, there were already alterations in the methylation profile of PBMCs in the obese subjects. Yi et al. [55] correlated high levels of homocysteine and SAH with a decrease in DNA methylation of lymphocytes. High SAH levels are also described as inhibitors of methyltransferase processes [20]. Furthermore, alterations in the transsulfuration pathway were also identified, in particular when comparing the HL with the HO subjects. This pathway is important in the redox state maintenance, since it is responsible for the synthesis of molecules, such as cysteine and GSH, that act as important ROS scavengers [14]. Interestingly, increased insulin levels have been shown to disrupt the redox homeostasis and consequently alter the defense mechanisms against excess ROS production [56]. Elshorbagy et al. [57] indicated that alterations in plasma total cysteine are associated with an increase in BMI, in adults. Similarly, we also observed higher levels of total cysteine in the HO subjects compared to their lean counterparts. On the other hand, the levels of fCysteine remain unchanged between groups even in the presence of higher levels of cystine, especially in the obese groups. Furthermore, the fCysteine/cystine ratio is an important plasma marker that defines an imbalance in redox homeostasis [58]. However, no differences were identified in that ratio between groups in our cohort, although the levels of GSSG were increased in the obese groups. Similar results were postulated by Choromańska et al. [59] who reported high levels of GSSG in plasma of obese adults with hypertension. While the fCysteine/cystine ratio was not altered, we observed significant alterations in the GSH/GSSG ratio in the HO group, with an increase in the percentage of GSH oxidation compared to HL. Few studies of this kind have been performed in prepubertal children with obesity; however, it is known that in adults, insulin resistance and dyslipidemia have a significant impact on oxidative stress [3, 60]. In line with the alterations detected in the redox homeostasis, increased plasma inflammation was also observed, in the children with obesity, in particular CRP and leptin levels. Pedersen et al. [61] demonstrated that hyperinsulinemia during obesity induced the expression of inflammation-associated genes. Importantly, the role leptin plays in inflammatory exacerbation has been reviewed [2]. Interestingly, high levels of homocysteine have also previously been implicated with inflammation [62].

These metabolic alterations are key factors for future development of obesity-related comorbidities, particularly in people with insulin resistance, and have been described, particularly in adults, as important predictive markers of CVD, including atherosclerosis [20, 54, 63]. Additionally, our data show significant differences in adipocyte secretion of adipokines. While leptin levels were increased in both obese groups, adiponectin levels were decreased in the UHO group. Our results are in agreement with other studies indicating low adiponectin levels, as a marker of adipocyte secretory dysfunction, together with elevated leptin secretion [64, 65]. Our data indicate an evident alteration in adipocyte secretory patterns, especially in the obese- and insulin-resistant subjects. Landgraf et al. [66] corroborate the presence of adipose tissue dysfunction during childhood obesity, making the correlation between adipocyte expansion and inflammation. Interestingly, the high leptin and low adiponectin levels together with the increased plasma lactate are early important markers that differentiate obese subjects with insulin resistance from their metabolic healthy counterpart. Our results show a deregulated adipocyte secretory function, even under normoglycemia in prepubertal children with obesity.

Moreover, our data also show significant changes in the methylation patterns of PBMCs. Circulating cells have been considered important surrogate markers for different type of diseases, since they are in contact with the continuous changing circulating molecules from surrounding tissues. Therefore, small metabolic changes that may occur in the body will induce key modifications in the biology of these cells [67]. In agreement, our results reflect the impact of obesity on the epigenetic patterns of PBMC DNA. In particular, we noticed alterations of the methylation degree in the vicinity of genes that are involved in one-carbon metabolism, such as the transmethylation and transsulfuration pathways. Previous studies have shown alterations in methylation patterns associated with obesity and obesity-related insulin resistance [6870]. The lack of significant differences in these metabolic pathways and in the DNA methylation results when comparing the HO and the UHO subjects might likely be due to the low study power when we split the obese subjects into two different groups. Our study is limited by relatively small sample sizes as we were unable to fully reach our recruitment goals. Unfortunately, we do not have the DNA methylation data for the entire cohort. In future studies, we would like to increase the in each group and invite the subjects back two years later for a follow-up visit.

Finally, the physiological and biochemical characteristics of the young prepubertal children show that among the obese subjects, there are indeed a group of HO that are still metabolically healthy as the HL, while the UHO already presented metabolic dysfunction, even under normoglycemia. This suggests that the adult obesity phenotype is starting to set in at an early age, raising important questions about the reversibility of this condition and the future health of these subjects.

In conclusion, our study showed the presence of deep alteration in one-carbon metabolism and related pathways in children with obesity (Figure 1), as young as 5-9 years of age. These alterations are mainly driven by obesity and elevated insulin levels. Most importantly, even the metabolic healthy obese subjects show increased levels of important metabolic markers, such as SAH, that is related with developing future comorbidities. This may reflect an important transition phase between different stages of obesity and insulin resistance. Specific markers that could identify this transition would be extremely important in identifying populations early at high risk of developing CVD and T2D. The data presented provide insight on important metabolic changes that occur in obesity early in life (summarized in Figure 1). Importantly, the metabolic healthy obese children seem to have some compensatory mechanisms to maintain important features unchanged, including the lipid profile and glucose levels. Also, some important alterations are already established at the epigenetic level, although some of these alterations could still be reversed by possible changing lifestyle habits [39].

Data Availability

Additional data can be found in supplementary tables.

Conflicts of Interest

The authors have no conflicts of interests to declare.

Authors’ Contributions

Shannon Rose and Eugenia Carvalho contributed equally to this work and are considered senior authors.

Acknowledgments

Research reported in this publication was financed by the National Institute of General Medical Sciences of the National Institutes of Health under a COBRE Award Number P20GM109096 and the Arkansas Biosciences Institute/Arkansas Children’ Research Institute, through the Discovery Acceleration Initiative/Program Project Planning Grants and bridging funds. Further funds were obtained by the European Regional Development Fund, through the Centro 2020 Regional Operational Programme Healthy Aging2020-CENTRO-01-0145-FEDER-000012 and through COMPETE 2020–Operational Programme for Competitiveness and Internationalisation and Portuguese national funds via FCT—Fundação para a Ciência e a Tecnologia, projects POCI-01-0145-FEDER-007440, POCI-01-0145-FEDER-022184 (GenomePT), UIDB/04539/2020, UIDP/04539/2020, UIDP/04501/2020 (iBiMED-UA), and PB PhD grant SFRH/BD/143849/2019. EB was partly funded by NIH/NIGMS UL1 TR003107 and KL2 TR003108 and USDA/ARS 6026-51000-012-06-S. The content is solely the responsibility of the authors and does not necessarily represent the official views of the funding agencies. We thank all the participants in the study. We acknowledge the technical assistance of Matthew Cotter, Oleksandra Pavliv, Stewart Macleod, Stefan Graw, and Stephanie Byrum.

Supplementary Materials

Supplementary 1. Supplementary Table 1: correlation of clinical and biochemical parameters (A) in healthy lean and healthy obese participants and (B) in healthy obese and unhealthy obese participants.

Supplementary 2. Supplementary Table 2: complete list of differentially methylated positions (DMPs).

Supplementary 3. Supplementary Table 3: GO term enrichment results for biological processes.