Systematic metabolomic studies identified adult adiposity biomarkers with acetylglycine associated with fat loss in vivo

1 Introduction

Obesity, defined as an excess of adiposity, is a global epidemic widely recognized as a leading cause of various metabolic disorders and cancers (Cirulli et al., 2019; Koenen et al., 2021; He et al., 2023). Excess central fat accumulation causing metabolic disturbance can exacerbate the risk of diabetes and cardiovascular diseases (Vasan et al., 2018; Marsh et al., 2023). Biomarkers reported to influence abdominal obesity have also been shown to be prognostic of early cardiometabolic risk independent of general obesity (Supriya et al., 2018). However, not all abdominal areas in which fat accumulates are associated with harm. It is known that android fat deposited around the waist increases the risk of metabolic disease, while gynoid fat deposited around the hips does not, and some studies suggest gynoid fat may have a beneficial effect (Vasan et al., 2018; Marsh et al., 2023). Therefore, uncovering novel biomarkers for fat distribution can provide more precise insight into fat accumulation and obesity-related diseases.

Metabolomics provides comprehensive profiling of distinct exogenous/endogenous small molecules functioning as intermediates or end products of cellular metabolism. Since metabolites represent the downstream expression of genomic, transcriptomic, and proteomic factors, their study can reveal biomarkers and pathways that link genotypes to phenotypes (Cirulli et al., 2019). Multiple studies have revealed crucial metabolites for obesity or regulating BF (Chen et al., 2015; Bogl et al., 2016; Ahmad et al., 2022; O’Keeffe et al., 2022), some of which were used as a basis for developing treatments and preventive interventions for obesity-related disorders (Okosun et al., 2015; Chen and Gerszten, 2020) and prediction models for visceral adipose tissue, an important prognostic factor that is difficult to measure in practice (Boone et al., 2022). Although there have been numerous efforts to identify metabolic biomarkers of obesity, these studies were largely limited by small sample sizes, the number of metabolites considered, and lack of replication/functional validation. Moreover, they primarily relied on proxy measures of body composition such as BMI and/or waist-to-hip ratio. These proxy measures may miss valuable information provided by the regional BF, and therefore hinder assessment of metabolic alterations and differential risk of metabolic disorders among individuals presenting adiposity differently throughout the body (Rangel-Huerta et al., 2019).

In the present study, we aim to identify shared and specific metabolites and their associated pathways related to various BF traits. The study has three main components: discovery, validation, and in vivo experiments. First, we analyzed untargeted liquid chromatography-mass spectrometry (LC-MS) metabolomics profiling on a sample of 517 Chinese women for discovery. Next, we verified the identified BF-related metabolites in two independent cohorts with various ethnicities and genders. Both the validated metabolites and those associated with all the studied BF traits were selected for more detailed examination of their biological functions and relevance as revealed in prior studies. One such metabolite, meeting both criteria, was chosen for in vivo experiments for further validation and functional significance in mice. Our study design and analytical approach are outlined in Supplementary Figure S1.

2 Materials and methods

2.1 Subjects and study design

A total of 517 unrelated peri-/post-menopausal Chinese women in the discovery cohort were recruited from the Third Affiliated Hospital of Southern Medical University (Guang Zhou, Guang Dong Province, China) between June 2014 and January 2018. The mean age was 52.9 years (standard deviation (SD) = 2.9). Two independent samples for validation were sampled from the Louisiana Osteoporosis Study (LOS) (Zhao et al., 2018). The first validation cohort consisted of 136 Caucasian women (mean aged 31.48 years (SD = 5.08)) sampled by a discordant phenotype design based on the top/bottom 20% of the hip bone mineral density Z-score (Zhao et al., 2018). The second validation cohort comprised 700 males (295 African Americans and 405 Caucasians) with mean age 37.8 (SD = 8.45) years overall. The study designs and sample characteristics are detailed in Supporting Material S1 and Supplementary Tables S1-S2, respectively. Signed consent forms were obtained from all participants. All clinical measurements and experiment procedures were subject to the Helsinki Declaration II and regulations of the Institutional Review Boards.

2.2 Clinical measurements

All respondents in the discovery and validation cohorts completed a questionnaire regarding demographics, lifestyle, dietary factors, and reproductive and medical history. Trained research staff collected clinical measurements with participants. The BF distribution phenotypes considered in this study were android fat/whole BF mass ratio (A/W ratio), gynoid fat/whole BF mass ratio (G/W ratio), android fat to gynoid fat ratio (A/G ratio), and whole BF percentage (W%) (Bogl et al., 2016). Daily calibrated DXA machines (discovery cohort: Lunar, GE Healthcare, Madison, WI, United States; validation cohorts: Hologic Inc., Bedford, MA, United States) were used to measure BF (Supplementary Figure S2). The detail of the region of interest is provided in Supporting Material S1.

2.3 Metabolite analysis

Untargeted LC-MS metabolite profiling was conducted on serum samples from all cohorts. The details of experimental protocols, which are in accordance with the suggestions of the Metabolomics Standards Initiative reporting standards (Sumner et al., 2007), have been previously described (Gong et al., 2021) and are provided in Supporting Material S2. A high-resolution tandem mass spectrometer TripleTOF5600plus (SCIEX, United Kingdom) at Lian-Chuan Biotechnology Co., Ltd (Hangzhou, China) was used to detect metabolites eluted from the column in the discovery cohort. LC-MS metabolomics platforms were also used to perform the metabolomic analyses on serum samples for both validation cohorts. The experimental protocols for the validation cohorts have been previously described (Zhao et al., 2018) and are listed in Supporting Material S2. All metabolites profiling with a unique or high-confidence annotation were log-transformed and auto-scaled prior to the subsequent analyses.

2.4 In vivo validation study in mice for acetylglycine

Acetylglycine consistently demonstrated an anti-obesity association with observed BF traits across sex and ethnicity. We performed an in vivo study that extended the work of Harper et al. (Harper et al., 2010) by exploring the effect of acetylglycine on high-fat diet-induced obese mice for further validation. A total of 61 female C57BL/6 J mice were randomized into five groups: 1) standard chow diet with vehicle (control; n = 12), 2) high fat (60% kcal from fat; n = 12) diet with vehicle (HFD; n = 12), 3) HFD with low-dosage acetylglycine (HFD + ACE500; 500 mg/kg; n = 12), 4) HFD with medium-dosage acetylglycine (HFD + ACE1000; 1,000 mg/kg; n = 12), and 5) HFD with high-dosage acetylglycine (HFD + ACE1500; 1,500 mg/kg; n = 13). We administered drinking water or acetylglycine (treatment) to corresponding mice by oral gavage every day starting at 8 weeks of age. Experimental measurements were recorded weekly during the 16-week intervention. We used micro-computed tomography (µ-CT) to measure total abdominal fat (between L1 and L5 vertebrae), visceral BF, and subcutaneous BF (Luu et al., 2009). The biochemical analysis involved the measurement of fasting plasma glucose (FPG) after 6-h and 12-h fasting, total cholesterol (TC), triglyceride (TG), high-density lipoprotein cholesterol (HDL), low-density lipoprotein cholesterol (LDL), alanine transaminase (ALT), and aspartate aminotransferase (AST). All experimental procedures are provided in Supporting Material S3.

2.5 Statistical analysis

We investigated the relationship between the metabolites and BF traits using two complementary approaches, seemingly unrelated regression (SUR) and sparse partial least squares (sPLS) regression. The most important BF-related metabolites were identified as those with variable importance in projection (VIP) score >1 in sPLS and SUR false discovery rate (FDR) q-value ≤0.2. The selection of cutoffs for the FDR q-value and VIP score usually depends on the study objective. In this study, we aimed to provide a robust list of potential BF-related metabolites for future replication, and therefore chose not to apply highly stringent significance criteria since the analysis is exploratory in nature. A pathway analysis was performed to identify biological functions with p-score <0.05. The Spearman and partial Spearman correlation tests were conducted in the validation cohorts to examine the associations between the metabolites and traits after adjustment for age, height, physical activity for both cohorts, and ethnicity in the second cohort. We also performed the tests stratified by ethnicity. p-value ≤0.05 was considered statistically significant. All analyses were detailed in Supporting Material S4. The metabolites associated with all BF traits in the discovery cohort or replicated the associations in the validation cohorts will be further discussed in our study (Supporting Material S5).

The Generalized Estimating Equations (GEE) model was applied to investigate the associations between acetylglycine and the measure of weight and weight gain in the animal models. The covariates included baseline weight, treatment (water, low-, medium-, or high-dose acetylglycine), and diet (normal/HFD) in the weight model, and treatment and diet in the weight gain model. Tukey method was used to examine whether there were significant pairwise differences (p-value ≤0.05) in weight and weight gain. The Kruskal–Wallis test was used as a global test for a one-time measurement to compare biomedical indexes across different treatment groups. The null hypothesis of this test is that there is no significant difference in the biomedical measurement across all treatment groups, while the alternative hypothesis is that at least one group differs significantly from at least one other group. Following the Kruskal–Wallis test, Wilcoxon tests were conducted to examine whether acetylglycine had significant differences (FDR q-value ≤0.05) in pairwise comparisons for both the µ-CT and biochemical measurements.

3 Results

3.1 Relationships between BF traits and serum metabolites in discovery

Among 3,075 annotated metabolites with MS1 and MS2 confidence levels in our discovery cohort (Supplementary Table S3 and Supporting Material S3), we identified 18 statistically significant BF-metabolite associations with FDR q-values <0.2 in the SUR models and VIP scores >1 in sPLS (Figure 1 and Supplementary Table S4). Figure 1A illustrates how the metabolites relate to multiple traits: none of the 18 metabolites exclusively related to a single BF trait. All were significantly associated with G/W ratio. Overall, two metabolites had significant relationships with three traits (G/W, A/W, and A/G ratios) representing central BF and five with all four traits representing general and central BF.

FIGURE 1

FIGURE 1. Summary results for the joint analysis in the discovery cohort. Panel (A) Venn diagram depicting metabolite sets that were significantly associated with body fat traits (FDR q-value <0.2 and VIP >1). Traits are represented by ovals of different colors. The number in each area represents the count of identified BF-associated metabolites associated with the given traits. Panel (B). Important metabolites selected on the basis of VIP, score >1 and FDR q-value ≤0.2. Each tile represents the effect size of the metabolite on each trait. VIP: the variable importance in projection scores, FDR: false discovery rate, A/W ratio: Android fat/Whole body fat ratio, G/W ratio: Gynoid fat/Whole body fat ratio, A/G ratio: Android fat to gynoid fat ratio, and W: Whole body fat %. *: p-value <0.05; **: p-value <0.01; ***: p-value <0.001.

Directions of coefficient effects on the four traits were consistent across the 18 metabolites; A/G ratio, A/W ratio, and W% were always affected in opposite directions as G/W ratio (Figure 1B). Those having positive relationships with A/G ratio, A/W ratio, and W% but a negative relation with G/W ratio were defined as obesity-risk metabolites, the others as protective metabolites. Among the six metabolites with VIP scores >2, the top three (austalide J, pyridinoline, and deoxycoformycin) are obesity-risk metabolites, while acetylglycine, creatinine, and cortisone are protective metabolites. The top six in the VIP score also had relatively large effect sizes on central BF compared with general BF (Supplementary Figure S3). Acetylglycine, for example, had stronger negative effect sizes on A/G ratio (β = −0.19, p-value = 1.09 × 10⁻⁵) and A/W ratio (β = −0.17, p-value = 6.23 × 10⁻⁵) than on W% (β = −0.09, p-value = 3.24 × 10⁻³).

3.2 Functional topological pathway analysis in discovery

Six metabolites out of 18 could be mapped to the Kyoto Encyclopedia of Genes and Genomes (KEGG) database and used in the pathway analysis. We identified 14 significant pathways, three modules, 28 enzymes, 39 reactions, and 21 compounds (refer to Table 1 for an overview and Supplementary Table S6 for a comprehensive list). The topological results showed that thiamine had the highest connectivity at the compound level (degree = 10) and thiamine metabolism had the highest connectivity at the pathway level (degree = 9) in Supplementary Figure S4. The connectivity degree indicates the importance of thiamine for essential metabolic processes.

TABLE 1

TABLE 1. Summary of topological attributes of the top 5% significant biological components at each biological level.

3.3 Metabolite and BF trait relations validation in LOS cohorts

We examined the correlations of four metabolites in the female validation cohort and five in the male validation cohort, which resulted in 22 significant relations between the metabolites and traits (Table 2). The correlations between acetylglycine and the traits were significant in both cohorts, with the strongest negative association with W% and A/W ratio respectively in Caucasian women (Spearman correlation ρ) = −0.34, partial Spearman correlation (ρ′) = −0.29, both p-value < 0.001) and African American and Caucasian men (ρ = −0.28, ρ' = −0.26, both p-value < 0.001). Glycocholic acid was negatively associated with A/G ratio, A/W ratio, and W% in African American and Caucasian men only, which may be attributed to the relatively small sample size of Caucasian women (Supplementary Tables S7-S14).

TABLE 2

TABLE 2. Summary of body-fat associated metabolites in the validation cohorts.

In the ethnicity-stratified analyses, we replicated the correlations of five metabolites related to at least one BF trait in African American or Caucasian men (Supplementary Table S5). Acetylglycine had the highest absolute correlation coefficients on the BF traits (except for the G/W ratio in the African American men). Cortisone had the significant negative correlations with A/G ratio, A/W ratio, and W% exclusively in the African American men; creatinine was also negatively related to W% in Caucasian men.

3.4 Relationship of acetylglycine with body weight, body weight gain, and abdominal fat mass in mice

The acetylglycine treatments with equivalent baseline weight and diet were significantly associated with an average decrease in body weight of 0.51 (p-value = 0.01) and 1.03 (p-value = 2.30 × 10⁻⁷) grams, respectively, in ACE1000 and ACE1500 groups, compared with groups without the treatment. Similarly, we found negative marginal effect sizes on weight gain while adjusting for diet (Table 3). Body weight and weight gain were lower in all the acetylglycine treatment groups and the control group than the HFD-only group (Supplementary Figure S5). Overall food intake presented no differences between the HFD-only and HFD with treatment groups during the intervention period.

TABLE 3

TABLE 3. Results of GEE models for the associations of acetylglycine on mouse weight and weight change.

Among the acetylglycine treatment groups, we observed significant differences in the µ-CT measures: abdominal fat, visceral fat, subcutaneous fat, and fat/weight ratio (Figure 2 and Supplementary Tables S1-S15). The HFD groups had significantly higher measurements than the control group, indicating the success of the obesity model. Despite similar HFD consumption, we observed a decreasing trend in the four µ-CT measures as the treatment dosage increased, and statistically significant differences between the HFD-only group and the HFD with the medium and high-dose groups.

FIGURE 2

FIGURE 2. In vivo study with acetylglycine. Panels (A–D) show the pairwise comparisons for the µ-CT measures among different groups. The abdominal region was defined as between L1–L5 lumbar vertebrae. Panels (E, F) are the pairwise comparisons for the fasting plasma glucose test. FDR: false discovery rate; *: FDR q-value <0.05; **: FDR q-value <0.01; ***: FDR q-value <0.001; ****: FDR q-value <0.0001.

3.5 Association of acetylglycine with biochemical measurements in mice

The overall tests indicated significant differences in the levels of glucose, CHO, HDL, and LDL among the groups (Supplementary Table S16). The HFD-only group in the 6-h FPG test had significantly higher glucose level (mean = 8.57 mmol/L, SD = 1.49) than the control, HFD + ACE500, and HFD + ACE1000 groups (Figure 2). In the 12-h FPG test, there were significant differences between the HFD + ACE1000 and control/HFD-only groups. ALT and AST, two enzymes for testing liver toxicity or damage, were not significantly different across groups (Supplementary Figure S6). Lipid profiles were significantly higher in the treatment groups with HFD compared to those in the control group, and there were no differences in the cholesterol levels among the intervention groups (Supplementary Figure S7). Taken together, the HFD-only diet successfully induced obesity in the mice, as reflected by the higher weight, weight gain, abdominal fat mass, and glucose level observed in the study.

4 Discussion

We performed a joint analysis that revealed the relationships between the metabolites and four DXA-derived BF measures among Chinese women. In contrast to proxy phenotypes, such as BMI, BF measures provide a direct assessment of body adiposity distribution, enabling the identification of novel metabolites related to general or central BF. We identified 18 metabolites that demonstrated consistent associations with A/G ratio, A/W ratio, and W%, and opposing associations with G/W ratio. We further partially validated the consistent relationships of the metabolites on the BF traits in two independent samples with different ethnicity and sex characteristics. Lastly, we validated the protective effect of acetylglycine in the obese mouse experiments.

Pyridinoline, deoxycoformycin, acetylglycine, creatinine, and cortisone were robustly associated with general and central BF. Austalide J and 28-homobrassinolide were only associated with central BF, and austalide J had the highest VIP score, reflecting a strong contribution in distinguishing abdominal obesity. Interestingly, no metabolite was exclusively related to any single BF trait, implying that the identified metabolites simultaneously influence total and central region fat deposits with different effect sizes. As gynoid fat may reduce the risk of metabolic diseases (Fu et al., 2013; Okosun et al., 2015), our findings of 12 metabolites positively associated with G/W ratio suggest further replication and merit for future study in preventing obesity.

Pyridinoline and deoxycoformycin were observed as positive associations with general and central BF, except gynoid fat, in our study. These two metabolites have been recognized from exogenous sources and have close relationships with complex diseases. Pyridinoline, a well-known biomarker of bone resorption and formation (Delmas et al., 1991; Kuo and Chen, 2017), has been reported to be positively associated with weight, BMI, and bone mineral density at different skeletal sites (New et al., 2000; Wang et al., 2016). It is a cross-linked compound that provides chemical stabilization in cartilage, bone, ligaments, and blood vessels (Kuo and Chen, 2017). Next, deoxycoformycin is an inhibitor of adenosine deaminase and a common anticancer chemotherapeutic drug for leukemia and lymphoma treatment (Tusup et al., 2022). In addition, deoxycoformycin can be produced by fungi such as C. militaris, which is regarded as a beneficial food supplement (Chen et al., 2020). However, it increases the susceptibility to toxic accumulation in kidney, liver, and central nervous system (Demain and Sanchez, 2009; Chen et al., 2020). The effect of the impaired tissue/organs by deoxycoformycin on the development of BF and its dose-dependent toxicity in humans is still unclear. Further research on the underlying mechanism of this relationship is warranted.

We demonstrated the potential anti-obesity effect of acetylglycine on both general and central body fat, which is consistent with previous findings of its association with lower BMI levels and fat percentages in population-based studies (Moore et al., 2014; Cirulli et al., 2019). However, a longitudinal study of Mexican American women aged 20–72 years reported that increased levels of acetylglycine elevated the risk of weight gain (Zhao et al., 2016). Further investigations into the long-term effect are needed to clarify these contradictory results. In addition to obesity-related traits, a study on the causal effect of metabolites on cardiovascular diseases found that acetylglycine is associated with a lower risk of diastolic blood pressure and is considered a potential cause (Qiao et al., 2021). The glycine N-acyltransferase (GLYAT) enzyme, linked with phenylalanine metabolism (Supplementary Figure S4), forms N-acylglycines (Westhuizen et al., 2000) and catalyzes glycine conjugation of acyl-Coa-species, producing intermediate products such as acetylglycine in amino acid and fatty acid metabolism (Badenhorst et al., 2013). A mice study showed that lower mRNA expression levels of GLYAT in adipose tissue were observed in a fat-susceptible group compared to a fat-resistant group (Fedry et al., 2016). The perturbation of GLYAT induces harmful disruption to CoA homeostasis (Badenhorst et al., 2013), musculoskeletal development (Badenhorst et al., 2013), and obesity-related metabolic disturbances (Fedry et al., 2016; Alves et al., 2019), suggesting that acetylglycine plays a role in the development of obesity-related diseases.

Creatinine, a breakdown product of creatine, reflects muscle metabolism and kidney function and plays a significant role in muscle energy production (Patel et al., 2013). Individuals with low muscle mass and body weight have low creatinine levels (Patel et al., 2013) and high glomerular hyperfiltration rates, which may increase the risk of metabolic diseases and diabetes (Hjelmesæth et al., 2010). The association between cortisone and general/central BF was consistent with a recent study reporting a negative association on BMI in European women (Cirulli et al., 2019). Interestingly, when assessing the relationship of cortisone in childhood obesity, cortisol and cortisone levels were significantly positively correlated with BMI and A/G ratio (Noppe et al., 2016).

The age-dependent differences in activities of 11β-hydroxysteroid dehydrogenases, an enzyme catalyzing the interconversion between cortisol and cortisone, may reverse the associations on BF deposition (Vierhapper et al., 2007; Chapman and Seckl, 2008). Cortisone participating in cortisol metabolism for generating steroid hormone is a known signature for obesity, type 2 diabetes (Gawlik et al., 2020), and other diseases (Supplementary Figure S4).

In our validation analyses, we observed the negative relationships for glycocholic acid, a key a bile acid regulator of fat absorption, cholesterol level, and energy homeostasis (Marco-Ramell et al., 2018). Glycocholic acid was only replicated in African American or Caucasian men but not in Caucasian women. This discrepancy could be attributed to aging-related decline in physiological functions, such as a decrease in female estrogen and hormone levels. A study reported a large difference in glycocholic acid levels between obese and lean women aged between 50 and 70 years old, and a relatively small difference in women aged 30–40 years old (Xie et al., 2015). This is consistent with the findings from our study samples. Further work is needed to understand the complex relationship between glycocholic acid and BF which may be mediated by age-related factors, estrogen, and hormone level.

The negative associations for N-Phenylacetyl-L-Glutamine and thiamine were validated in the men’s and women’s replication samples, respectively. A recent study reported that there are no sex differences in excretion of the amino acid N-Phenylacetyl-L-Glutamine (Zheng et al., 2014). The abnormal activity of N-Phenylacetyl-L-Glutamine could be attributed to the difference in phenylalanine metabolism and further reflect the development of obesity (Bogl et al., 2016). However, some studies also observed that the concentration was lower in obese men compared with normal-weight men (Yu et al., 2018). The mechanism through which N-Phenylacetyl-L-Glutamine may influence BF requires further investigation. Moreover, many studies support thiamine as an essential micronutrient in glucose metabolism that is negatively associated with obesity (Maguire et al., 2018). To the best of our knowledge, our study provided the first evidence of the relationship between thiamine and BF.

The arginine and proline metabolism pathway had a high closeness and betweenness centralities, which implies that the pathway has considerable influence on neighboring molecules/pathways. The pathway is known to be associated with obesity (Cirulli et al., 2019). Elevated arginine levels can reduce BF accretion in humans and animals (Wu et al., 2009). Further investigation of the identified biological components involved in the arginine and proline metabolism may reveal an important influence on obesity.

The ATP-binding cassette (ABC) transporter pathway, a pivotal connector pathway, regulates the import/export of membrane proteins, such as thiamine, phthalate, L-glutamine, and glycine via different enzymes or reactions for the overall BF-related network (Rees et al., 2009). The pathway may delineate associations between genes and physiologic changes before, during, and after the onset of obesity in mice (Donepudi et al., 2016) and beneficial influences on obese women (Teixeira et al., 2020). However, we observed that phthalate in ABC transporters was negatively associated with A/G, A/W, and W%, in accordance with the findings from Hatch et al.’s study of BMI and waist circumference (Hatch et al., 2008). Other studies found contradictory results, demonstrating that phthalate was an endocrine-disrupting chemical (Apau et al., 2020) and linked to adverse health outcomes such as obesity and diabetes (Hatch et al., 2008; Apau et al., 2020). The various characteristic of phthalates on obesity may depend upon endogenous hormone levels and vary across sex and age groups (Hatch et al., 2008). The influence on BF and other pathways via the ABC transporter pathway may warrant further investigation.

Additionally, our study suggested that impaired cholesterol metabolism may cause adverse health outcomes. The cholesterol metabolism pathway was connected to pathways for aldosterone-regulated sodium reabsorption and prostate cancer. Cortisone and cortisol were closely connected to the disease pathway and may be diagnostic biomarkers for obesity. The impaired conversion of cortisone and cortisol has been reported as an important factor for insulin sensitivity and central obesity (Weaver et al., 1998).

Harper et al. observed no statistical differences in body weights and feed consumption between the acetylglycine treatment groups and the control group (Harper et al., 2010). Our animal model extended the prior research by exploring the effect of acetylglycine on HFD-induced obese mice. For the first time, we demonstrated a protective effect of mid and high doses of acetylglycine on body weight and fat in mice. The evidence has also been seen in a case of smoking-cessation-induced weight gain mice study that showed acetylglycine ameliorated weight gain rate compared with HFD control mice (Fluhr et al., 2021). The single-cell transcriptomics analysis of epididymal-adipose immune cells exemplified that acetylglycine is a potent signaling molecule that modulates modulated multiple adipose-tissue gene expressions in obesity-associated pathways such as immune response, lysosome function, and tissue remodeling (Fluhr et al., 2021).

Biochemical tests showed significant increases in the fasting glucose, total cholesterol, HDL, and LDL of the HFD groups, and unaltered triglycerides. This is consistent with what was reported in a previous HFD-induced obesity study in mice (Eisinger et al., 2014; Fluhr et al., 2021). The results of the FPG tests for the mid-dose group provided evidence for a reduction in the glucose levels, and the total and abdominal fat mass. Acetylglycine may ameliorate fat mass accumulation and further improve glucose metabolism since the decreased capacity for adipocyte differentiation and angiogenesis is reported to alleviate lipogenesis and lipolysis activities as well as insulin resistance (Patel and Abate, 2013). Elevated acetylglycine levels have been associated with decreased risk of impaired fasting glucose and onset of diabetes (Menni et al., 2013) and improved glucose tolerance (Fluhr et al., 2021). However, this evidence may not be applicable to all situations, such as high variations in glucose reduction in our high-dose group. Further investigations are needed to estimate the dose-response effects of acetylglycine on glucose.

Our study has several strengths. We had precise measures of BF mass from DXA scans to distinguish the association of metabolites with obesity. Our study systematically examined the BF-metabolites associations and partially validated them in two independent samples with different sex and ethnicities. We confirmed our novel finding with in vivo experiments showing that acetylglycine can significantly affect adiposity. One inherent limitation is the difficulty harmonizing metabolites from different metabolomics platforms. Although we partially replicated our findings in other sex and ethnicity groups, more research on obesity is still needed since metabolic responses affect or are affected by diverse lifestyles and diets. Second, we cannot draw causality in our findings due to the cross-sectional design. However, our validations from the independent replication cohort and in vivo mice studies may support the robustness of our findings to some extent.

Our investigation of adiposity phenotypes systematically identified BF metabolomic signatures and several relevant pathways. Our study provides evidence that acetylglycine, creatinine, and cortisone may have a protective role against body fat accumulation. External validation replicated six metabolites associated with BF in either Caucasian women or African American and Caucasian men. The protective effects of acetylglycine were consistent across different samples and were further validated in vivo in mice. The findings open new possibilities for utilizing acetylglycine as a potential diagnostic biomarker and therapeutic target of obesity or obesity-related diseases.

Data availability statement

The data that support the findings of this study are available from the corresponding author upon request and approval of the team and respective institutions. The Kyoto Encyclopedia of Genes and Genomes database and network documents for the topological pathway analysis have been deposited on the Mendeley data website: https://data.mendeley.com/datasets/rhm2b8hz75

Ethics statement

The studies involving human participants were reviewed and approved by Tulane University Institutional Review Board and Third Affiliated Hospital of Southern Medical University Institutional Review Board. The patients/participants provided their written informed consent to participate in this study.

Author contributions

The authors’ responsibilities were as follows—H-WD: Supervision, Funding acquisition. JiS: Supervision, Funding acquisition. H-MX: Project administration, Funding acquisition. K-JS: Conceptualization, Methodology, Writing—Original Draft, Writing—Review and Editing, Visualization, Formal analysis, Investigation. X-YC: Validation, Resources, Investigation for the animal study. RG and XL: Data Curation, Conceptualization, Resources. QZ, JG, and HuS: Conceptualization, Writing—Review and Editing. S-DH, M-CF, YL, and Y-HZ: Resources. QT: Data Curation for the validation data. All authors provided critical feedback and helped shape the research, analysis, and manuscript.

Funding

H-WD and HuS were partially supported by grants from National Institutes of Health [R01AR069055, U19AG055373]. JiS was partially supported by grants from the Science and Technology Program of Guangzhou, China [201604020007], and the National Natural Science Foundation of China [81770878]. H-MX was partially supported by the National Key R&D Program of China (2016YFC1201805 and 2017YFC1001100). This research was supported in part using high-performance computing (HPC) resources and services provided by Technology Services at Tulane University, New Orleans, LA.

Acknowledgments

The authors would like to acknowledge the English editing service from Loula Burton, Tulane Office of Research.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2023.1166333/full#supplementary-material

References

Ahmad, S., Hammar, U., Kennedy, B., Salihovic, S., Ganna, A., Lind, L., et al. (2022). Effect of general adiposity and central body fat distribution on the circulating metabolome: A multicohort nontargeted metabolomics observational and mendelian randomization study. Diabetes 71, 329–339. doi:10.2337/db20-1120

PubMed Abstract | CrossRef Full Text | Google Scholar

Alves, A., Bassot, A., Bulteau, A.-L. L., Pirola, L., and Morio, B. (2019). Glycine metabolism and its alterations in obesity and metabolic diseases. doi:10.3390/nu11061356

CrossRef Full Text | Google Scholar

Apau, J., Sefah, W., and Adua, E. (2020). Human contact with phthalates during early life stages leads to weight gain and obesity. doi:10.1080/23312009.2020.1815273

CrossRef Full Text | Google Scholar

Badenhorst, C. P. S., Van Der Sluis, R., Erasmus, E., and Van Dijk, A. A. (2013). Glycine conjugation: Importance in metabolism, the role of glycine N-acyltransferase, and factors that influence interindividual variation. Expert Opin. Drug Metab. Toxicol. 9, 1139–1153. doi:10.1517/17425255.2013.796929

PubMed Abstract | CrossRef Full Text | Google Scholar

Bogl, L. H., Kaye, S. M., Rämö, J. T., Kangas, A. J., Soininen, P., Hakkarainen, A., et al. (2016). Abdominal obesity and circulating metabolites: A twin study approach. Metabolism 65, 111–121. doi:10.1016/j.metabol.2015.10.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Boone, S. C., van Smeden, M., Rosendaal, F. R., le Cessie, S., Groenwold, R. H. H., Jukema, J. W., et al. (2022). Evaluation of the value of waist circumference and metabolomics in the estimation of visceral adipose tissue. Am. J. Epidemiol. 191, 886–899. doi:10.1093/AJE/KWAB298

PubMed Abstract | CrossRef Full Text | Google Scholar

Chapman, K. E., and Seckl, J. R. (2008). 11beta-HSD1, inflammation, metabolic disease and age-related cognitive (dys)function. Neurochem. Res. 33, 624–636. doi:10.1007/s11064-007-9504-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, B., Sun, Y., Luo, F., and Wang, C. (2020). Bioactive metabolites and potential mycotoxins produced by cordyceps fungi: A review of safety. Toxins (Basel) 12, 410. doi:10.3390/TOXINS12060410

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, H. H., Tseng, Y. J., Wang, S. Y., Tsai, Y. S., Chang, C. S., Kuo, T. C., et al. (2015). The metabolome profiling and pathway analysis in metabolic healthy and abnormal obesity. Int. J. Obes. 39, 1241–1248. doi:10.1038/ijo.2015.65

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Z.-Z., and Gerszten, R. E. (2020). Metabolomics and proteomics in type 2 diabetes. Circ. Res. 126, 1613–1627. doi:10.1161/CIRCRESAHA.120.315898

PubMed Abstract | CrossRef Full Text | Google Scholar

Cirulli, E. T., Guo, L., Leon Swisher, C., Shah, N., Huang, L., Napier, L. A., et al. (2019). Profound perturbation of the metabolome in obesity is associated with health risk. Cell Metab. 29, 488–500.e2. doi:10.1016/J.CMET.2018.09.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Delmas, P. D., Schlemmer, A., Gineyts, E., Riis, B., and Christiansen, C. (1991). Urinary excretion of pyridinoline crosslinks correlates with bone turnover measured on iliac crest biopsy in patients with vertebral osteoporosis. J. Bone Min. Res. 6, 639–644. doi:10.1002/JBMR.5650060615

PubMed Abstract | CrossRef Full Text | Google Scholar

Demain, A. L., and Sanchez, S. (2009). Microbial drug discovery: 80 years of progress. J. Antibiot. 621 62, 5–16. doi:10.1038/ja.2008.16

PubMed Abstract | CrossRef Full Text | Google Scholar

Donepudi, A. C., Cheng, Q., Lu, Z. J., Cherrington, N. J., and Slitt, A. L. (2016). Hepatic transporter expression in metabolic syndrome: Phenotype, serum metabolic hormones, and transcription factor expression. Drug Metab. Dispos. 44, 518–526. doi:10.1124/dmd.115.066779

PubMed Abstract | CrossRef Full Text | Google Scholar

Eisinger, K., Liebisch, G., Schmitz, G., Aslanidis, C., Krautbauer, S., and Buechler, C. (2014). Lipidomic analysis of serum from high fat diet induced obese mice. Int. J. Mol. Sci. 15, 2991–3002. doi:10.3390/ijms15022991

PubMed Abstract | CrossRef Full Text | Google Scholar

Fedry, J., Blais, A., Even, P. C., Piedcoq, J., Fromentin, G., Gaudichon, C., et al. (2016). Urinary metabolic profile predicts high-fat diet sensitivity in the C57Bl6/J mouse. J. Nutr. Biochem. 31, 88–97. doi:10.1016/J.JNUTBIO.2015.12.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Fluhr, L., Mor, U., Kolodziejczyk, A. A., Dori-Bachash, M., Leshem, A., Itav, S., et al. (2021). Gut microbiota modulates weight gain in mice after discontinued smoke exposure. Nature 600, 713–719. doi:10.1038/s41586-021-04194-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, X., Song, A., Zhou, Y., Xiaoguang, M., Jiao, J., Yang, M., et al. (2013). Association of regional body fat with metabolic risks in Chinese women. Public Health Nutr. 17, 2316–2324. doi:10.1017/S1368980013002668

PubMed Abstract | CrossRef Full Text | Google Scholar

Gawlik, A. M., Shmoish, M., Hartmann, M. F., Wudy, S. A., and Hochberg, Z. (2020). Steroid metabolomic signature of insulin resistance in childhood obesity. Diabetes Care 43, 405–410. doi:10.2337/DC19-1189

PubMed Abstract | CrossRef Full Text | Google Scholar

Gong, R., Xiao, H.-M., Zhang, Y.-H., Zhao, Q., Su, K.-J., Lin, X., et al. (2021). Identification and functional characterization of metabolites for bone mass in peri- and postmenopausal Chinese women. J. Clin. Endocrinol. Metab. 106, e3159–e3177. doi:10.1210/clinem/dgab146

PubMed Abstract | CrossRef Full Text | Google Scholar

Harper, M. S., Amanda Shen, Z., Barnett, J. F., Krsmanovic, L., Dakoulas, E. W., and Delaney, B. (2010). Toxicology studies with N-acetylglycine. Food Chem. Toxicol. 48, 1321–1327. doi:10.1016/j.fct.2010.02.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Hatch, E. E., Nelson, J. W., Qureshi, M. M., Weinberg, J., Moore, L. L., Singer, M., et al. (2008). Association of urinary phthalate metabolite concentrations with body mass index and waist circumference: A cross-sectional study of nhanes data, 1999-2002. Environ. Heal. A Glob. Access Sci. Source 7, 27–15. doi:10.1186/1476-069X-7-27

PubMed Abstract | CrossRef Full Text | Google Scholar

He, X., Zhu, Z., Zang, J., Wang, Z., Liao, P., Wang, W., et al. (2023). Percent body fat, but not body mass index, is associated with cardiometabolic risk factors in children and adolescents. Chronic Dis. Transl. Med. doi:10.1002/CDT3.54

CrossRef Full Text | Google Scholar

Hjelmesæth, J., Røislien, J., Nordstrand, N., Hofsø, D., Hager, H., and Hartmann, A. (2010). Low serum creatinine is associated with type 2 diabetes in morbidly obese women and men: A cross-sectional study. BMC Endocr. Disord. 10, 6. doi:10.1186/1472-6823-10-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Koenen, M., Hill, M. A., Cohen, P., and Sowers, J. R. (2021). Obesity, adipose tissue and vascular dysfunction. Circ. Res. 128, 951–968. doi:10.1161/CIRCRESAHA.121.318093

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuo, T. R., and Chen, C. H. (2017). Bone biomarker for the clinical assessment of osteoporosis: Recent developments and future perspectives. Biomark. Res. 5, 18–19. doi:10.1186/s40364-017-0097-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Luu, Y. K., Lublinsky, S., Ozcivici, E., Capilla, E., Pessin, J. E., Rubin, C. T., et al. (2009). In vivo quantification of subcutaneous and visceral adiposity by micro-computed tomography in a small animal model. Med. Eng. Phys. 31, 34–41. doi:10.1016/j.medengphy.2008.03.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Maguire, D., Talwar, D., Shiels, P. G., and McMillan, D. (2018). The role of thiamine dependent enzymes in obesity and obesity related chronic disease states: A systematic review. Clin. Nutr. ESPEN 25, 8–17. doi:10.1016/j.clnesp.2018.02.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Marco-Ramell, A., Tulipani, S., Palau-Rodriguez, M., Gonzalez-Dominguez, R., Miñarro, A., Jauregui, O., et al. (2018). Untargeted profiling of concordant/discordant phenotypes of high insulin resistance and obesity to predict the risk of developing diabetes. J. Proteome Res. 17, 2307–2317. doi:10.1021/acs.jproteome.7b00855

PubMed Abstract | CrossRef Full Text | Google Scholar

Marsh, M. L., Oliveira, M. N., and Vieira-Potter, V. J. (2023). Adipocyte metabolism and health after the menopause: The role of exercise. Nutr 15, Page 444 15, 444. doi:10.3390/NU15020444

CrossRef Full Text | Google Scholar

Menni, C., Fauman, E., Erte, I., Perry, J. R. B., Kastenmüller, G., Shin, S. Y., et al. (2013). Biomarkers for type 2 diabetes and impaired fasting glucose using a nontargeted metabolomics approach. Diabetes 62, 4270–4276. doi:10.2337/db13-0570

PubMed Abstract | CrossRef Full Text | Google Scholar

Moore, S. C., Matthews, C. E., Sampson, J. N., Stolzenberg-Solomon, R. Z., Zheng, W., Cai, Q., et al. (2014). Human metabolic correlates of body mass index. Metabolomics 10, 259–269. doi:10.1007/s11306-013-0574-1

PubMed Abstract | CrossRef Full Text | Google Scholar

New, S. A., Robins, S. P., Campbell, M. K., Martin, J. C., Garton, M. J., Bolton-Smith, C., et al. (2000). Dietary influences on bone mass and bone metabolism: Further evidence of a positive link between fruit and vegetable consumption and bone health? Am. J. Clin. Nutr. 71, 142–151. doi:10.1093/AJCN/71.1.142

PubMed Abstract | CrossRef Full Text | Google Scholar

Noppe, G., Van Den Akker, E. L. T., De Rijke, Y. B., Koper, J. W., Jaddoe, V. W., and Van Rossum, E. F. C. (2016). Long-term glucocorticoid concentrations as a risk factor for childhood obesity and adverse body-fat distribution. Int. J. Obes. 40, 1503–1509. doi:10.1038/ijo.2016.113

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Keeffe, L. M., Bell, J. A., O’Neill, K. N., Lee, M. A., Woodward, M., Peters, S. A. E., et al. (2022). Sex-specific associations of adiposity with cardiometabolic traits in the UK: A multi-life stage cohort study with repeat metabolomics. Plos Med. 19, e1003636. doi:10.1371/JOURNAL.PMED.1003636

PubMed Abstract | CrossRef Full Text | Google Scholar

Okosun, I. S., Seale, J. P., and Lyn, R. (2015). Commingling effect of gynoid and android fat patterns on cardiometabolic dysregulation in normal weight American adults. Nutr. Diabetes 5, e155. doi:10.1038/nutd.2015.5

PubMed Abstract | CrossRef Full Text | Google Scholar

Patel, P., and Abate, N. (2013). Body fat distribution and insulin resistance. Nutrients 5, 2019–2027. doi:10.3390/nu5062019

PubMed Abstract | CrossRef Full Text | Google Scholar

Patel, S. S., Molnar, M. Z., Tayek, J. A., Ix, J. H., Noori, N., Benner, D., et al. (2013). Serum creatinine as a marker of muscle mass in chronic kidney disease: Results of a cross-sectional study and review of literature. J. Cachexia. Sarcopenia Muscle 4, 19–29. doi:10.1007/s13539-012-0079-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiao, J., Zhang, M., Wang, T., Huang, S., Zeng, P., Tang, J., et al. (2021). Evaluating causal relationship between metabolites and six cardiovascular diseases based on GWAS summary statistics 12. doi:10.3389/fgene.2021.746677

PubMed Abstract | CrossRef Full Text | Google Scholar

Rangel-Huerta, O. D., Pastor-Villaescusa, B., and Gil, A. (2019). Are we close to defining a metabolomic signature of human obesity? A systematic review of metabolomics studies. Metabolomics 15, 93. doi:10.1007/s11306-019-1553-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Rees, D. C., Johnson, E., and Lewinson, O. (2009). ABC transporters: The power to change. Nat. Rev. Mol. Cell Biol. 10, 218–227. doi:10.1038/nrm2646

PubMed Abstract | CrossRef Full Text | Google Scholar

Sumner, L. W., Amberg, A., Barrett, D., Beale, M. H., Beger, R., Daykin, C. A., et al. (2007). Proposed minimum reporting standards for chemical analysis: Chemical analysis working group (CAWG) metabolomics standards initiative (MSI). Metabolomics 3, 211–221. doi:10.1007/s11306-007-0082-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Supriya, R., Tam, B. T., Yu, A. P., Lee, P. H., Lai, C. W., Cheng, K. K., et al. (2018). Adipokines demonstrate the interacting influence of central obesity with other cardiometabolic risk factors of metabolic syndrome in Hong Kong Chinese adults. PLoS One 13, e0201585. doi:10.1371/JOURNAL.PONE.0201585

PubMed Abstract | CrossRef Full Text | Google Scholar

Teixeira, M. D., Tureck, L. V., Do Nascimento, G. A., de Souza, R. L. R., and Furtado-Alle, L. (2020). Is it possible abc transporters genetic variants influence the outcomes of a weight-loss diet in obese women? Genet. Mol. Biol. 43, e20190326. doi:10.1590/1678-4685-GMB-2019-0326

PubMed Abstract | CrossRef Full Text | Google Scholar

Tusup, M., Kündig, T. M., and Pascolo, S. (2022). Epitranscriptomics modifier pentostatin indirectly triggers Toll-like receptor 3 and can enhance immune infiltration in tumors. Mol. Ther. 30, 1163–1170. doi:10.1016/J.YMTHE.2021.09.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Vasan, S. K., Osmond, C., Canoy, D., Christodoulides, C., Neville, M. J., Di Gravio, C., et al. (2018). Comparison of regional fat measurements by dual-energy X-ray absorptiometry and conventional anthropometry and their association with markers of diabetes and cardiovascular disease risk. Int. J. Obes. 42, 850–857. doi:10.1038/ijo.2017.289

PubMed Abstract | CrossRef Full Text | Google Scholar

Vierhapper, H., Heinze, G., and Nowotny, P. (2007). Sex-specific difference in the interconversion of cortisol and cortisone in men and women. Obesity 15, 820–824. doi:10.1038/oby.2007.592

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Liu, D., Li, Y., Guo, L., Cui, Y., Zhang, X., et al. (2016). Metabolomic analysis of serum from obese adults with hyperlipemia by UHPLC-Q-TOF MS/MS. Biomed. Chromatogr. 30, 48–54. doi:10.1002/bmc.3491

PubMed Abstract | CrossRef Full Text | Google Scholar

Weaver, J. U., Taylor, N. F., Monson, J. P., Wood, P. J., and Kelly, W. F. (1998). Sexual dimorphism in 11 β hydroxysteroid dehydrogenase activity and its relation to fat distribution and insulin sensitivity; a study in hypopituitary subjects. Clin. Endocrinol. (Oxf). 49, 13–20. doi:10.1046/j.1365-2265.1998.00494.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Westhuizen, F. H. van der, Pretorius, P. J., and Erasmus, E. (2000). The utilization of alanine, glutamic acid, and serine as amino acid substrates for glycine N-acyltransferase. J. Biochem. Mol. Toxicol. 14, 102–109. Available at: https://onlinelibrary.wiley.com/doi/abs/10.1002/%28SICI%291099-0461%282000%2914%3A2%3C102%3A%3AAID-JBT6%3E3.0.CO%3B2-H?sid=nlm%3Apubmed (Accessed June 26, 2021).

PubMed Abstract | Google Scholar

Wu, G., Bazer, F. W., Davis, T. A., Kim, S. W., Li, P., Marc Rhoads, J., et al. (2009). Arginine metabolism and nutrition in growth, health and disease. Amino Acids 37, 153–168. doi:10.1007/s00726-008-0210-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, G., Wang, Y., Wang, X., Zhao, A., Chen, T., Ni, Y., et al. (2015). Profiling of serum bile acids in a healthy Chinese population using UPLC–MS/MS. J. Proteome Res. 14, 850–859. doi:10.1021/pr500920q

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, H. T., Fu, X. Y., Xu, B., Zuo, L. L., Ma, H. B., and Wang, S. R. (2018). Untargeted metabolomics approach (UPLC-Q-TOF-MS) explores the biomarkers of serum and urine in overweight/obese young men. Asia Pac. J. Clin. Nutr. 27, 1067–1076. doi:10.6133/apjcn.052018.07

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, H., Shen, J., Djukovic, D., Daniel-MacDougall, C., Gu, H., Wu, X., et al. (2016). Metabolomics-identified metabolites associated with body mass index and prospective weight gain among Mexican American women. Obes. Sci. Pract. 2, 309–317. doi:10.1002/osp4.63

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Q., Shen, H., Su, K.-J. J., Zhang, J.-G. G., Tian, Q., Zhao, L.-J. J., et al. (2018). Metabolomic profiles associated with bone mineral density in US Caucasian women. Nutr. Metab. 15, 57. doi:10.1186/s12986-018-0296-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, H., Yde, C. C., Arnberg, K., Mølgaard, C., Michaelsen, K. F., Larnkjær, A., et al. (2014). NMR-based metabolomic profiling of overweight adolescents: An elucidation of the effects of inter-/intraindividual differences, gender, and pubertal development. Biomed. Res. Int. 2014, 537157. doi:10.1155/2014/537157

PubMed Abstract | CrossRef Full Text | Google Scholar