1. Introduction
Globally, there were 2.5 billion adults living with overweight and obesity [
1]. In the last 12 years, Malaysia has seen a 10% rise in the prevalence of overweight and obesity [
2] and this trend is not unique to Malaysia but is evident in many countries globally. The persistent rise in global obesity prevalence remains a significant public health concern worldwide, given its well-established association with an increased risk of developing chronic conditions, such as type 2 diabetes (T2D), hypertension (HPT) and cardiovascular diseases (CVD) [
3,
4]. Efforts to combat obesity through various weight loss interventions are widespread. However, the outcomes associated with weight loss intervention vary significantly among individuals.
Variations in the clinical manifestations of obesity could be one of the attributes that lead to differences in the responses towards obesity intervention regimens. Weight loss resulting from lifestyle interventions among overweight/obese women categorized as metabolically healthy has been linked to changes in lipid metabolism, activation of sulfation processes, and modulation of microbiota metabolism, potentially indicating a metabolically protective effect [
5]. However, gaps persist in understanding disparities in metabolite profiles among obesity phenotypes and the underlying mechanisms that could contribute to differential responses to interventions. Comprehensive investigation is required to clarify the regulatory mechanisms that may vary across different obesity phenotypes.
Metabolomics provides a platform that allows us to capture the dynamic physiological conditions corresponding to current health conditions by analyzing low-molecular-weight metabolites present in tissues or biological fluids, such as lipids, amino acids, peptides, organic acids, and carbohydrates [
6,
7]. This comprehensive profiling of metabolites offers valuable insights into the etiology of obesity and other diseases, facilitating the discovery of potential biomarkers that could enhance our current knowledge of obesity pathophysiology and its related comorbidities.
Hence, the aims of our study were to explore the differences in serum metabolite profiles between women classified as metabolically healthy obese/overweight (MHO) and metabolically unhealthy obese/overweight (MUO) and examine the changes in these metabolite profiles following a 6-month lifestyle intervention for weight loss.
2. Materials and Methods
2.1. Study Design and Study Participants
This study utilized archived samples collected from an extension study of a lifestyle intervention for weight loss known as My Body is Fit and Fabulous at Home (MyBFF@home). The study was conducted in 2015 and the samples were stored at -80°C freezer. Information on the study design and recruitment details has been published elsewhere [
8]. Participants in this study consisted of obese and overweight housewives aged 18 - 59 years old, living in the low-cost flats in Klang Valley, Malaysia. All the participants underwent six months lifestyle intervention for weight loss program consisting of dietary counselling, physical activity (PA) and self-monitoring tools (PA diary, food diary and pedometer). Fasting blood samples, anthropometric and clinical measurements were assessed at baseline and after 6 months of intervention. Biochemical profiles were analyzed in the laboratory according to methods described previously [
9]. Dietary intakes were recorded using 3-day food diary and nutrient intakes were calculated using Nutritionist Pro TM version 2.4 (First Data Bank, The Hearst Corp, NY USA) as per described elsewhere [
10]. A total of 70 (MHO=36, MUO=34) serum were randomly selected from the archived follow-up samples of MyBFF@home. Participants were considered as metabolically healthy obese (MHO) based on the following criteria: HbA1c <6.5% and systolic or diastolic blood pressure of <140mmHg, <90mmHg respectively, while participants were considered as metabolically unhealthy obese (MUO) by the following criteria: HbA1c ≥6.5%, systolic or diastolic blood pressure of ≥140mmHg, ≥90mmHg respectively, and self-reported of being diagnosed with T2D or HPT [
11].
2.2. Sample Preparation for Metabolite Profiling
All samples were kept at -80°C prior to the analysis. Serum samples were thawed on ice to minimize metabolite degradation. Serum samples were vortexed and centrifuged (20,000g x 5min at 4°C). 200μL of the serum supernatant were mixed with 400μL of phosphate buffer (KH2PO4), pH 7.4 in deuterium oxide (D2O) and 0.1% TSP (Merck, Darmstadt, Germany) and 0.1% imidazole (Sigma-Aldrich, St. Louis, MO, USA) in a 1:1 ratio.
2.3. NMR Metabolomic Analysis
Untargeted metabolomic approach were used to analyze the metabolites present in the study. The NMR analysis method was adapted and optimized based on the work of Maulidani et al. [
12]. 1D 1H-NMR spectra were collected at a temperature of 26°C using a 600 MHz Jeol NMR. The combination of PRESAT with CPMG pulse sequence were used to suppress water signals and wide protein resonances. NMR spectra with a spectral width of 10ppm were acquired using a total of 128 scans and a 660-second acquisition time. Spectra were processed using the Chenomx NMR suite version 9.0 software (Chenomx Inc., Edmonton, AB, Canada) with the following settings: 0.50 Hz line broadening, autophasing, baseline correction (Whittaker spline), referenced to TSP as an internal standard, and referenced to imidazole as a pH indicator. The spectral band between 0.50 and 10.00 ppm were divided into equal bins using intelligent binning (0.04 ppm). The peak ppm readings were calibrated against the 0 ppm TSP signal. The area of the spectrum associated with residual water and imidazole were eliminated prior to analysis. Subsequently, the corresponding spectra were transformed to a table of common integrals that has a non-negative value for multivariate data analysis (MVDA).
2.4. Statistical Analysis
Prior to MVDA, the spectral data were mean-centered and pareto-scaled to improve normality. Principal component analysis (PCA) and partial least squares discriminant analysis (OPLS-DA) were performed to detect outliers and observe trends and separation of metabolites between groups using the standard algorithm as implemented in the SIMCA® software version 17.0.2 (Sartorius, Göttingen, Germany). The validation of the OPLS-DA model was carried out using cross-validated analysis of variance (CV-ANOVA), with the results expressed as p-values for the model. Prior to univariate data analysis the metabolomics data were log-transformed and metabolites not found in at least 20% of the samples were removed. Statistical analysis was conducted using IBM SPSS Statistics for Windows, Version 27.0 (Armonk, NY, USA). The normality of continuous data was determined using Shapiro-Wilk test. Differences in the baseline parameters of the participants between the two groups were analyzed using independent t-tests. Generalized estimated equation (GEE) was used to assess significantly changed metabolites between groups (time X group). The analysis was adjusted for sociodemographic characteristics (age, education level and household income), baseline value of anthropometry parameters (BMI and WC), FPG, HbA1c, systolic blood pressure, the used of medication and baseline value of dietary intake (energy, carbohydrate, protein, total fat, saturated fat, sodium, potassium, and dietary fiber. The Benjamini–Hochberg (B-H) method was applied to correct for multiple testing with false discovery rate (FDR) set at 5%. To assess the extent of changes experienced after the intervention for all participants, delta values were calculated. These values were derived from the ratio of the sixth month measurements to the baseline values of significant metabolites, anthropometric measures, blood pressure, and biochemical variables. Pearson correlation analysis was then performed using the delta values to evaluate the relationship between the significant metabolites and anthropometry (weight, BMI and WC), biochemical (FPG, cholesterol, HDL, LDL and TG) and clinical data (systolic and diastolic blood pressure).
2.5. Pathway Analysis
To explore the metabolic pathways affected by changes in metabolites following the intervention, pathway analysis was conducted in the MetaboAnalyst 6.0 [
13] web application, which utilizes databases for Homo sapiens from the Kyoto Encyclopedia of Genes and Genomes (KEGG) and the Human Metabolome Database (HMDB) as the pathway library. The Materials
4. Discussion
Our study demonstrated that prior to the intervention, the spectral regions 3.42ppm and 3.46ppm were highly predictive of differentiating between MHO and MUO. Following the intervention, regions 3.42ppm and 3.74ppm were highly predictive of the separation between MHO and MUO. Based on the spectral data generated and profiled in this study, these regions corresponded to glucose which was found to be significantly different between MHO and MUO both at baseline and sixth month of the intervention. Initial investigation of metabolite differences at baseline identified three metabolites were significantly different between the MHO and MUO groups. However, after six months of intervention a greater number of metabolites (15 metabolites) showed significant differences between the groups. Among these were the branched-chain amino acids (BCAAs): leucine, isoleucine, and valine which have been shown to be associated with insulin resistance, metabolic syndrome and T2D in various cohorts [
14,
15,
16]. This aligns with our findings, where the levels of BCAAs were significantly higher in MUO women compared to MHO women. The proposed mechanisms linking BCAAs and T2D involve elevated BCAAs levels triggering activation of the mTOR pathway. This activation increases the production of lipid intermediates such as diacylglycerol (DAG), which subsequently results in the phosphorylation of insulin receptor substrate (IRS) proteins, impairing downstream insulin signalling [
17,
18]. The frequent association between BCAAs and obesity and T2D has led to suggestions that MHO and MUO could be characterized based on these metabolites [
19].
Further analysis revealed that seven metabolites significantly changed following the intervention between MHO and MUO women namely TMAO, arginine, ribose, aspartate, carnitine, choline and tyrosine. All these metabolites displayed a significant decreasing trend in the MHO group. The substantial decrease in circulating choline, carnitine, and TMAO were likely be attributed to reduced intake of protein-rich foods such as red meat, which are high in choline and carnitine during the intervention. TMAO is produced by the gut microbiota from these precursors [
20], hence a lower intake of choline- and carnitine-rich foods putatively leads to decreased TMAO levels. This direct relationship was demonstrated previously, showing that the consumption of fish and meat is associated with plasma concentrations of TMAO [
21]. Additionally, the marked reduction in arginine levels in MHO could also be linked to the decreased protein intake as studies have shown that dietary intake of arginine is significantly associated with serum and plasma arginine levels, suggesting that a lower intake of arginine-rich foods, during the intervention may contribute to this observed decrease [
22,
23].
Tyrosine is an aromatic amino acid categorized as non-essential since it can be synthesized in our body from its precursor, phenylalanine. Our study does not only demonstrate that tyrosine levels were markedly elevated in MUO post-intervention but was also found to be higher in the MUO compared to MHO at baseline. Similar to BCAAs, elevated levels of tyrosine have been widely linked to insulin resistance and an increased risk of T2D [
16,
24,
25]. This aligns with our findings, as some participants in the MUO group either had diagnosed T2D or exhibited elevated HbA1c levels beyond the normal range. Additionally, tyrosine has been linked to aging [
26], which was also evident in our study as we found a positive correlation between tyrosine levels and age.
Despites the substantial changed observed in the metabolomic profile of MHO and MUO women in the present study, only minimal change was observed in the anthropometry parameter while no significant changes were observed in the clinical parameters. Although the changes in body weight and WC across the groups were minimal but the change were significant. In this study, carnitine levels decreased in MHO, but increased in MUO following the intervention which was accompanied by significant reduction in body weight and WC in the latter. Additionally, a significant correlation was found between carnitine and both body weight and BMI. A meta-analysis previously has revealed that supplementation with L-carnitine significantly reduced body weight [
27]. This could pertain to the role of carnitine in fatty acid metabolism. Carnitine is crucial for transporting fatty acids into the mitochondria, where they undergo β-oxidation and are broken down into short and medium chain fatty acids to produce energy [
28]. The increased energy expenditure through fatty acid metabolism is presumed to be the factor that promoted weight loss in the MUO women. Additionally, during periods of reduced calorie intake and weight loss, it is possible that more carnitine was released from the tissues [
29]. This could potentially account for the high levels of carnitine observed in MUO individuals despite their reduced protein intake. The same pattern of changes and correlation with weight loss was also observed with tyrosine. Previous research has shown weight reduction is associated with decreased tyrosine levels in overweight/obese [
30] and those with metabolic syndrome [
31]. However, in the present study tyrosine levels was elevated in the MUO after the intervention and was significantly correlated with weight reduction. A study on 3-week weight loss program also reported significant weight loss with increase tyrosine levels following intervention. While elevated levels of tyrosine in the MUO have often been reported [
32,
33,
34] and are associated with insulin resistance and T2D [
16,
24,
25], the discrepancy in the findings could potentially be attributed to several factors including differences in the study population, participants’ dietary intake and the methods used to measure the metabolites.
We conducted a pooled correlation analysis to investigate the association between the significant metabolites and various anthropometric and clinical parameters using the delta data and found TMAO, arginine, carnitine and ribose were correlated with lipid parameters namely TG, LDL, TC and HDL respectively. These metabolites are mostly involved in the energy metabolism, lipid metabolism and vascular function. TMAO were found to be negatively correlated with both BMI and TG. These correlations were also reported in other studies involving overweight individuals, though they found a positive correlation [
35,
36]. The correlation might be related to the role of the gut microbiota, as it has been shown to also affect changes in TG levels and BMI [
37]. As previously mentioned, the circulating levels of TMAO are influenced by the dietary intake of protein-rich food and its biosynthesis by the gut microbiota [
20]. In this study, participants across all groups exhibited a significant reduction in protein intake and corresponding decreased in TMAO levels. This suggests that changes in dietary protein intake could influence circulating TMAO levels by modulating the metabolism of the gut microbiota [
38,
39]. Another factor that could also influence the circulating levels of TMAO is the activity of the flavin monooxygenase 3 (FMO3) enzyme that may link the relationship between TMAO and lipid regulation. FMO3 does not only responsible for the conversion of TMA to TMAO in the liver but it also involves in the regulation of lipid via the farnesoid X receptor (FXR) and liver X receptor (LXR) which resulted in reduced reverse cholesterol transport to the intestine and leads to accumulation of cholesterol [
40,
41].
In the context of the positive correlation between arginine and LDL, it's plausible that changes in arginine could be a consequence of obesity rather than the effect on arginine itself on LDL changes as it has been reported in a meta-analysis that supplementation with arginine did not have any significant impact on LDL [
42]. Apart from excessive body weight, participants in the present study exhibited both impaired fasting glucose and hyperlipidemia. This aligns with the well-documented coexistence of obesity with insulin resistance, metabolic syndrome, inflammation, and oxidative stress [
43]. Chronic inflammation and metabolic abnormalities associated with obesity are thought to increase the activity of arginase 1, an enzyme that competes with nitric oxide (NO) synthase for arginine [
44,
45]. This competition reduces NO production, which is crucial for maintaining healthy vascular function. The decrease in NO leads to increased production of reactive oxygen species (ROS), further promoting oxidative stress and endothelial dysfunction [
44,
45]. Consequently, this could result in lipid dysregulation, including elevated LDL levels [
46,
47], as observed in our study. As previously suggested, the significant reduction in carnitine may be attributed to decreased protein intake during the intervention. Therefore, the observed association between carnitine and total cholesterol (TC) might also be due to metabolic disturbances associated with obesity. It is worth noting that studies have shown that supplementation with L-carnitine significantly reduces TC levels [
48,
49] which contradicted with our findings.
We also managed to demonstrate correlation between ribose-FPG and ribose-HDL. Unfortunately, no previous human study has reported such correlation. While small amount of ribose can be sourced from dietary intake, our body primary source of ribose is synthesized through pentose phosphate pathway (PPP) [
50]. Obesity is known to be associated with elevated activity of glucose-6-phosphate dehydrogenase (G6PD), the rate-limiting enzyme in the pentose phosphate pathway (PPP). The increase in G6PD activity led to enhanced oxidative PPP activity, resulting in higher NADPH production. The surplus NADPH fueled NADPH oxidase (NOX)-mediated reactive oxygen species (ROS) generation, which exacerbated inflammatory response and induced DNA damage [
51,
52,
53,
54,
55]. The non-oxidative phase of the pentose phosphate pathway (PPP) produces ribose-5-phosphate (R5P), a crucial building block for nucleotide synthesis [
56]. Obesity-induced DNA damage could putatively increase the non-oxidative PPP activity to boost nucleotide synthesis necessary for DNA repair leading to reduced circulating levels of ribose as observed in the present study.
This study has certain limitations. We did not measure the inflammatory cytokines associated with obesity; thus, we could not verify the association between obesity-related inflammation and the metabolites. Additionally, since the samples were obtained from a community-based intervention program, we lacked access to the medical history of participants with comorbidities (T2D and HPT), particularly the duration since their diagnosis. This is a notable limitation, as existing evidence indicates that metabolite levels significantly change as these diseases progress [
57]. Consequently, we were unable to control for this factor in our analysis. Furthermore, dietary intake data were recorded using self-reported 3-day food diaries, which may lead to underreporting. However, to minimize this issue, each participant's food diary was reviewed one-on-one with a professional dietician or nutritionist at each visit.