Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Endocrinol., 31 July 2025

Sec. Obesity

Volume 16 - 2025 | https://doi.org/10.3389/fendo.2025.1625501

This article is part of the Research TopicTargeting Adipose Tissue for the Treatment of Metabolic AlterationsView all 7 articles

Integrated proteomic and metabolomic profiling identifies distinct molecular signatures and metabolic pathways associated with obesity and potential targets for anti-obesity therapies

Yi Li,&#x;Yi Li1,2†Huawu Yang&#x;Huawu Yang3†Xinpeng Zhang&#x;Xinpeng Zhang4†Xingyu HeXingyu He5Anke LiuliAnke Liuli5Rui LiRui Li6Xingyu HanXingyu Han5Yongmei Li*Yongmei Li2*Pan Gao*Pan Gao7*
  • 1Department of Radiology, The Third People’s Hospital of Chengdu, Chengdu, Sichuan, China
  • 2Department of Radiology, The First Affiliated Hospital of Chongqing Medical University, Chongqing, China
  • 3The Center of Obesity and Metabolic Diseases, Department of General Surgery, The Third People’s Hospital of Chengdu, Chengdu, Sichuan, China
  • 4General Surgery Day Ward, Department of General Surgery, The Third People’s Hospital of Chengdu, Chengdu, Sichuan, China
  • 5College of Life Science and Engineering, Southwest Jiaotong University, Chengdu, Sichuan, China
  • 6School of Basic Medical Sciences, Southwest Medical University, Luzhou, Sichuan, China
  • 7Obesity and Metabolism Medicine-Engineering Integration Laboratory, Department of General Surgery, The Third People’s Hospital of Chengdu, Chengdu, Sichuan, China

Background: Adipose tissue remodeling induced by bariatric surgery plays a pivotal role in promoting weight loss and metabolic improvement. However, the underlying molecular mechanisms, particularly protein-metabolite regulatory networks, remain poorly understood. This integrative proteomic and metabolomic study identifies key pathway alterations and molecular signatures associated with metabolic phenotypes, offering novel mechanistic insights into the therapeutic efficacy of bariatric surgery.

Methods: Visceral adipose tissue samples were analyzed using label-free DIA quantitative proteomics and LC-MS/MS metabolomics. Proteomic and metabolomic data were processed with MaxQuant software and XCMS R package, respectively.

Results: Proteomic and metabolomic analyses were performed on visceral adipose tissue from 10 obese patients undergoing sleeve gastrectomy and 10 controls. Proteomic profiling quantified identified 135 differentially expressed proteins (57 upregulated, 78 downregulated), with PHACTR2 and PLIN2 upregulated in obesity and ADAR down-regulated in obesity. Enrichment analyses indicated disruptions in lipid droplet formation, muscle processes, and protein autophosphorylation, with KRT1/MYH9 and NF1/ATR identified as hub proteins. Metabolomics revealed 191 differential metabolites (110 upregulated, 81 downregulated), with 4-Vinylcyclohexene positively correlated with BMI and asparagine-betaxanthin negatively correlated. KEGG analysis showed disturbances in purine/pyrimidine metabolism, AMPK signaling, and cortisol biosynthesis. Integrated protein-metabolite network analysis identified OSBPL10, CUL2, and PRTN3 as potential regulators of lipid metabolism and insulin resistance, offering insights into obesity-associated metabolic dysfunction.

Conclusions: This study integrated proteomic and metabolomic data from visceral adipose tissue obtained through sleeve gastrectomy, identifying obesity-related functional pathways and molecular signatures linked to metabolic phenotypes, highlighting the value of multi-omics in understanding adipose tissue remodeling and postoperative metabolic improvement.

Introduction

Obesity has emerged as a major global public health challenge, currently affecting more than 878 million adults worldwide and closely associated with the development of numerous chronic diseases (1), including type 2 diabetes, cancer, and cardiovascular disorders (24). Bariatric surgery, particularly laparoscopic sleeve gastrectomy (SG), is recognized as one of the most effective clinical interventions for severe obesity, capable of inducing substantial weight loss and improving metabolic health (5). Although the clinical benefits of SG (including weight reduction, improvements in glycemic and lipid profiles) have been widely reported (68), the molecular adaptation mechanisms within adipose tissue that mediate these long-term metabolic improvements remain poorly understood.

As a critical organ regulating metabolic homeostasis, adipose tissue plays a central role in the development of obesity and its associated metabolic dysfunctions (9). In obesity, adipose tissue function becomes impaired, characterized by insulin resistance and dysregulated lipid metabolism (10, 11). Weight loss induced by bariatric surgery partially restores adipose tissue function and enhances insulin sensitivity (12). While these phenomena have been observed, a comprehensive understanding of the proteomic and metabolomic changes in adipose tissue following SG, as well as their associations with long-term metabolic outcomes, remains limited. Most existing multi-omics studies have focused primarily on blood samples or subcutaneous adipose tissue, while relatively little attention has been paid to visceral adipose tissue (VAT), which is closely associated with insulin resistance and metabolic risk.

Recent advances in proteomics and metabolomics have begun to uncover molecular pathway alterations associated with obesity and weight loss (13, 14). For example, studies have reported changes in sphingolipid metabolism and inflammatory signaling pathways following bariatric surgery (15). However, these studies often lack longitudinal clinical data and rarely integrate proteomic and metabolomic profiles with detailed postoperative phenotypic outcomes. In this study, we performed an integrated proteomic and metabolomic analysis of visceral adipose tissue from 10 obese patients undergoing SG and matched 10 healthy controls. By combining molecular data with one-year clinical follow-up information, we aimed to systematically characterize adipose tissue-specific pathway alterations associated with weight loss and metabolic improvement, and to identify potential molecular signatures linked to long-term outcomes. The uniqueness of this study lies in its first demonstration of the associations between adipose tissue proteomic and metabolomic signatures and clinical outcomes such as body mass index (BMI) reduction and ​metabolic dysfunction-associated steatotic liver disease (MASLD) remission after SG. These findings provide novel insights into the mechanisms of adipose tissue remodeling in the context of obesity intervention. Moreover, they deepen our understanding of the molecular basis of bariatric surgery and offer a foundation for individualized postoperative prognosis assessment and the development of new therapeutic targets.

Materials and methods

Subjects and study design

This study enrolled 30 subjects from the Third People’s Hospital of Chengdu, including 20 obese patients (8 males and 12 females) who underwent SG as part of the Longitudinal Study of Bariatric Surgery in Western China (registration number: ChiCTR2300073353), and another 10 healthy control participants (5 males and 5 females). All subjects met the following inclusion criteria: 1) obesity group (BMI > 28 kg/m²), control group (BMI < 24 kg/m²); 2) age between 18 and 65 years. Besides, we systematically excluded obese subjects diagnosed with polycystic ovary syndrome (PCOS). Within three days prior to surgery, several clinical phenotypic information was systematically collected for each subject, including BMI, total cholesterol (TC), low-density lipoprotein (LDL), high-density lipoprotein (HDL), triglycerides (TG), fasting blood glucose (Glu), aspartate aminotransferase (AST), alanine aminotransferase (ALT), uric acid, urea nitrogen, creatinine, and other metabolic indices, as well as clinical information such as history of MASLD and diabetes.

For postoperative follow-up, the 20 SG patients were monitored longitudinally for 1 year, with regular assessments of body weight. Venous blood samples were drawn at 1-year post-surgery, and several metabolic-related clinical parameters were measured to evaluate the metabolic improvement resulting from the surgical intervention. For adipose tissue sample selection, we prioritized obese patients with complete 1-year follow-up data to ensure prognostic reliability. To maintain balanced group sizes for comparative analysis, we ultimately selected 10 SG patients from the original 20 obese participants for VAT collection, matched with VAT samples from 10 controls. All tissue specimens were collected, processed, and stored under standardized protocols to ensure sample integrity and consistency. This study was conducted in strict accordance with the ethical principles of the Declaration of Helsinki. The study protocol was reviewed and approved by the Ethics Committee of Chengdu Third People’s Hospital.

Samples collection

Standardized biospecimen collection and processing protocols were followed. All subjects fasted for 12 hours and abstained from medication for 24 hours prior to blood sampling. In the morning, trained medical personnel collected 5 mL of venous blood from each subject, and the samples were immediately stored at 4 °C. During the sleeve gastrectomy procedure, the surgical team collected omental adipose tissue samples according to standard procedures. Excised tissues were immediately transferred on dry ice and, under biosafety conditions, divided into two portions: one for DIA proteomic analysis (including library preparation and sample analysis) and one for untargeted metabolomic analysis. Any remaining tissue was processed and stored at -80 °C to ensure availability for future supplementary assays and validation studies. All procedures strictly followed standardized operating protocols to ensure sample quality and data reliability.

Proteome analysis

Peptide and protein identification was performed using MaxQuant (v2.6.7) against the UniProt homo sapiens reference proteome database (16, 17). Key parameters were defined in the configuration file mqpar.xml, with the following settings: DIA and MS2-based quantification modes; spectral library-based quantification; label-free quantification enabled; enzyme specificity set according to the digestion protocol; and “Match Between Runs” enabled to reduce missing values across samples. Intensity-based absolute quantification (iBAQ) was used for protein abundance normalization. The false discovery rate (FDR) was controlled at 1% at both the peptide and protein levels. Identifications based on single peptides and low-quality spectra were excluded. Following identification, raw protein intensities were normalized using median centering and batch correction. Differential expression analysis was initially performed using unadjusted P < 0.05 and |log>FC| ≥ 1 as cutoffs. Proteins meeting these criteria were considered high-confidence differentially expressed proteins (DEPs).

Metabolome analysis

The raw mass spectrometry data were first converted into mzXML format using ProteoWizard MSConvert (18). Subsequent data processing, including peak detection, retention time correction, peak alignment, and peak area extraction, was performed using the XCMS package (19). Specifically, peak detection was conducted using the centWave algorithm with the following parameters: ppm = 15 and peak width = 5–10 seconds. Peak alignment was carried out with a bandwidth (bw) of 5 and a binSize of 0.025. Ion features were retained for downstream analysis if more than 50% of samples in at least one group exhibited nonzero intensity measurements (20). Metabolites identification was performed according to the standards proposed by the Metabolomics Standards Initiative. Metabolite annotation was conducted by integrating information from the Human Metabolome Database (HMDB) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) (21, 22), based on accurate mass matching with a mass error threshold of less than 5 ppm. Significantly altered metabolites were identified based on the following criteria: P < 0.05, and variable importance in projection (VIP) than 1.

Pathway analysis and protein–protein interaction network

Functional enrichment analysis was performed to investigate the biological significance of the DEPs and metabolites. Gene Ontology (GO) enrichment analysis, covering Biological Process (BP), Cellular Component (CC), and Molecular Function (MF) categories, was conducted using the clusterProfiler R package (23). KEGG pathway analysis was similarly performed to identify significantly enriched pathways, with a P < 0.05 considered statistically significant. For protein-protein interaction (PPI) network construction, DEPs were input into the STRING database (v11.5) with a minimum interaction score threshold of 0.4 (24). The resulting network was visualized using Cytoscape software (v3.9.1) (25). Key hub proteins were identified based on network topology parameters such as degree centrality and betweenness centrality.

Statistical analysis

All statistical analyses were conducted using R software (v4.3.2) in this study. Continuous variables were expressed as mean ± sd or median with interquartile range (IQR) as appropriate. Paired-sample t-tests or Wilcoxon signed-rank tests were used for pre- and post-surgical comparisons within the same group. Categorical variables were analyzed using the chi-square test. Pearson correlation coefficients were calculated to assess associations between molecular features (proteins or metabolites) and clinical phenotypes. For multiple comparisons, p-values were adjusted using the Benjamini-Hochberg FDR method.

Results

Clinical characteristics of participants

This study analyzed comprehensive clinical data from a cohort of 30 participants, comprising 20 obese patients undergoing SG and 10 healthy controls. The clinical data encompassing key metabolic indicators such as BMI, TC, LDL, HDL, TG, ALT, AST, Glu, uric acid, creatinine, urea nitrogen, and documented histories of MASLD and diabetes (Supplementary Table S1). Longitudinal follow-up data were obtained over 12 months for all surgical patients, capturing metabolic changes (n = 12) and weight trajectories (n = 20) (Table 1). Comparative analysis revealed several interesting findings. No significant differences were observed in age between the obese and control groups, minimizing potential demographic bias. As expected, preoperative BMI in the obese group was significantly higher than in controls, while the BMI values declined markedly postoperatively, approaching those of healthy individuals, demonstrating the efficacy of SG in inducing substantial weight loss.

Table 1
www.frontiersin.org

Table 1. Clinical characteristics of participants before and after bariatric surgery and healthy controls.

Compared to controls, obese patients exhibited significantly reduced HDL levels preoperatively, a finding that improved significantly after surgery. TG levels were elevated in the obesity group, although the difference did not reach statistical significance, and a marked reduction in TG levels was observed one year after SG. Moreover, postoperative improvements were evident across several biochemical parameters, including significant reductions in Glu, ALT, AST, and uric acid, suggesting systemic metabolic amelioration following weight loss surgery. Notably, all 20 obese participants presented with MASLD prior to SG. One year postoperatively, only 6 of the 12 followed-up patients continued to exhibit signs of hepatic steatosis, indicating partial reversal of MASLD in a subset of patients. These results underscore the therapeutic impact of SG not only on weight but also on broader metabolic health, including liver function and lipid regulation.

Comprehensive proteomic changes in visceral adipose tissue between obesity and controls

To systematically investigate proteomic alterations in VAT associated with obesity, we performed label-free DIA quantitative proteomic profiling on VAT samples obtained from 10 obese patients undergoing bariatric surgery and 10 control individuals undergoing abdominal surgery for non-metabolic indications. A total of 4813 proteins were identified and quantified across all samples based on database search and stringent filtering criteria. Principal component analysis (PCA) based on the global proteomic profiles revealed a clear separation between the obesity and controls, with principal components 1 (PC1) and principal components 2 (PC2) explaining 13.49% and 13.2% of the total variance, respectively (Figure 1A). This suggests a robust divergence in protein expression patterns associated with the obese state. A total of 135 proteins were found to be differentially expressed between the two groups, meeting the criteria of statistical significance (P < 0.05) and a fold change threshold (|log<FC| > 1) (Supplementary Table S2). Of these, 57 proteins were upregulated and 78 were downregulated in the obese group relative to the controls (Figure 1B). Hierarchical clustering of these DEPs further demonstrated distinct expression signatures between obesity and controls (Figure 1C). Next, violin plots highlighted several representative proteins with significant expression differences. For example, PHACTR2, PLIN2, and SERPING1 were markedly upregulated in the obese group, whereas ADAR, PLTP and RBMS1 exhibited higher expression in controls (Figure 1D). These protein-level alterations may reflect key molecular mechanisms underlying VAT dysfunction in obesity and provide potential targets for further mechanistic studies.

Figure 1
Panel A shows a PCA plot differentiating samples by type; Ob (purple) and Ctl (green). Panel B is a volcano plot with log2 fold change versus -log10 FDR, highlighting significant genes. Panel C presents a heatmap indicating expression levels across samples, categorized by diagnosis and BMI. Panel D displays violin plots comparing expression for specific genes (PHACTR2, PLIN2, SERPING1, ADAR, PLTP, and RBMS1) between Ob and Ctl groups, with significant differences marked by asterisks.

Figure 1. Identification of differentially expressed proteins in visceral adipose tissue between obesity post-bariatric surgery and healthy-weight control groups (Ob, Obesity; Ctr, Control). (A) Principal component analysis (PCA) of protein expression profiles across samples. (B) Volcano plot of DEPs with thresholds set at |log>FC| > 1 and P < 0.05. (C) Hierarchical clustering heatmap of DEPs. (D) Violin plots illustrating expression levels of representative DEPs between two groups. * means P<0.05.

Protein functional pathways alterations associated with obesity

To further elucidate the biological implications of proteomic alterations in VAT of obese individuals, we conducted GO and KEGG pathway enrichment analyses based on the set of DEPs. We showed the 15 representative GO biological processes that were significantly enriched (P < 0.05), indicating functional convergence of obesity-related proteins (Figure 2A, Supplementary Table S3). Of particular note was the lipid droplet organization pathway, indicating that lipid processing and remodeling in VAT were significantly affected in the obese state. Besides, our results also show that obesity-related DEPs are related to muscle system processes and protein autophosphorylation functions. In obesity, skeletal muscle undergoes tissue remodeling through ectopic lipid deposition and inflammatory response, leading to metabolic disorders, and is closely related to the development of insulin resistance and MASLD (26). In addition, the role of protein phosphatases is essential for a variety of physiological responses, and protein kinase-phosphoprotein interactions also play an important role in obesity (27). These findings are consistent with the central role of adipose tissue in systemic energy homeostasis and lipid storage. Similarly, we displayed the top 15 significantly KEGG pathways (Figure 2B, Supplementary Table S4), among which the top three pathways, adherens junction, endocrine system, and salivary secretion were upregulated in obesity compared to healthy controls. These pathways are functionally related to epithelial barrier dysfunction, appetite regulation, and salivary organ secretion in obesity (2830), all of which are known to be altered in obesity-induced adipose tissue dysfunction.

Figure 2
Panel A shows a density plot of GO pathways with Log2FoldChange values ranging from -6 to 4, highlighting processes like muscle contraction. Panel B displays a similar plot for KEGG pathways with a range of -4 to 2, emphasizing pathways like adherens junction. Panels C and D present network diagrams illustrating gene interactions. Nodes are color-coded and connected by lines indicating relationships, with central genes highlighted.

Figure 2. Functional enrichment analysis of DEPs between obesity and control groups. (A) Bubble chart of significantly enriched Gene Ontology terms based on DEPs. (B) Bubble chart of significantly enriched KEGG pathways. (C) Protein-protein interaction (PPI) network of upregulated proteins in obesity (nodes: proteins; edges: interactions; node color: expression fold change). (D) PPI network of downregulated proteins in obesity, highlighting key hub proteins.

To investigate the functional interplay among the DEPs, we constructed protein-protein interaction (PPI) networks for both the upregulated and downregulated protein subsets. The network of upregulated proteins demonstrated a dense interaction architecture with KRT1 and MYH9 emerging as key hub nodes, suggesting their potential core roles in the functional transformation of VAT in obesity (Figure 2C). In contrast, the PPI network of downregulated proteins highlighted NF1, ATR, and LRRK2 as major interaction hubs (Figure 2D), indicating proteins that play a key role in the lean state. These results collectively underscore the complexity and specificity of proteomic remodeling in VAT during obesity and reveal functional clusters that may contribute to adipose tissue pathology and systemic metabolic dysregulation.

Untargeted metabolomic profiling of visceral adipose tissue with obesity

To investigate the metabolic reprogramming of VAT in obesity, we conducted untargeted metabolomic profiling using liquid chromatography-tandem mass spectrometry (LC-MS/MS) in both positive and negative ionization modes. A total of 6303 metabolites were identified, including 2572 metabolites in positive mode and 3731 metabolites in negative mode, which were merged for downstream analysis. Metabolites annotation and classification were carried out based on the Human Metabolome Database, with the chemical taxonomy analysis revealing diverse classes of metabolites, including lipids and lipid-like molecules (41.36%), organic acids and derivatives (17.88%), organoheterocyclic compounds (10.71%) and etc (Figure 3A). We next identified 191 significantly altered metabolites through differential expression analysis (VIP > 1.0, P < 0.05), including 110 upregulated and 81 downregulated metabolites in the obese group compared to the controls (Figure 3B, Supplementary Table S5). PCA analysis based on metabolite intensities revealed a clear separation between obese and control groups, with PC1 and PC2 accounting for 22.1% and 12.8% of total variance, respectively (Figure 3C). Further supervised multivariate analyses using PLS-DA and OPLS-DA confirmed robust group discrimination. In PLS-DA, component 1 and component 2 explained 21.7% and 9.2% of the variance, respectively (Figure 3D). In OPLS-DA, component 1 and component 2 explained 24.8% and 10.4% (Figure 3E), further supporting metabolic divergence between the two groups.

Figure 3
A composite image with seven panels. A: Pie chart showing distribution of chemical classes, with lipids and lipid-like molecules as the largest group at 41.36%. B: Volcano plot of differentially expressed compounds, with red dots indicating upregulation, gray dots non-significant, and blue dots downregulation. C: PCA plot showing distinct clusters for Ob and Ctl groups. D: PLS-DA plot depicting separation of Ob and Ctl groups. E: OPLS-DA plot with clear group differentiation. F: Bar chart of VIP scores for top compounds influencing group separation. G: Box plots of normalized expression for six highlighted compounds in Ob and Ctl groups.

Figure 3. Liquid chromatography-mass spectrometry-based metabolomics analysis of obesity and controls groups. (A) Pie chart of chemical taxonomy super classes annotated using the Human Metabolome Database. (B) Volcano plot of differential metabolites (VIP > 1.0, P < 0.05). (C) PCA score plot of metabolite expression profiles. (D) Partial least squares-discriminant analysis (PLS-DA) score plot. (E) Orthogonal PLS-DA (OPLS-DA) score plot. (F) Bubble chart of top 25 metabolites ranked by VIP scores. (G) Box plots of representative differentially expressed metabolites.

Ranking by VIP scores highlighted the top 25 metabolites (Figure 3F), with Diniconazole, 4-Vinylcyclohexene, (R)-(+)-Ropivacaine, Octanamide and Ropivacaine emerging as the top 5 metabolites. Among the representative metabolites, we showed that the levels of 5 metabolites, including the above ones, were increased in the obese group, whereas Asparagine-betaxanthin was notably downregulated (Figure 3G). These findings reveal significant alterations in VAT metabolomic profiles associated with obesity, indicating potential metabolic biomarkers and dysregulated pathways involved in adipose tissue dysfunction.

Functional pathways analysis of visceral adipose tissue metabolites associated with obesity

To clarify the biological implications of the metabolomic changes observed in VAT from obese patients, we conducted functional annotation and pathway enrichment analyses using the KEGG database. Our classification of all identified differential metabolites via KEGG pathway analysis demonstrated their participation in a wide array of biological processes (Figure 4A). Specifically, the majority of metabolites were predominantly enriched in five major KEGG functional categories: Metabolism, Human Diseases, Environmental Information Processing, Organismal Systems, and Cellular Processes. Notably, the “Metabolism” category was predominant, with a particularly high representation in pathways related to amino acid metabolism (n = 20), nucleotide metabolism (n = 19), and lipid metabolism (n = 18) (Supplementary Table S6). These results underscore the profound metabolic reprogramming that occurs in obese VAT. The top 15 significantly enriched pathways, visualized by bubble plot (Figure 4B), further highlight key metabolic and disease-related processes. Purine metabolism, nucleotide metabolism, and pyrimidine metabolism were the most significantly enriched pathways, with the highest number of associated metabolites and the greatest statistical significance. Additional pathways of interest included tryptophan metabolism, AMPK signaling, and cortisol biosynthesis, suggesting potential disruptions in amino acid turnover, cellular energy sensing, and endocrine regulation in obesity. Circular visualization of pathway enrichment offers an integrated overview of pathway statistics, including the number of involved metabolites (inner ring), enrichment significance (−log10 P values; middle ring), and specific pathway IDs (outer ring) (Figure 4C). Consistent with the bubble plot, several metabolic pathways (e.g., ko00230: purine metabolism, ko00240: pyrimidine metabolism, ko00380: tryptophan metabolism) and human disease pathways (e.g., ko04924: cortisol synthesis and secretion) were highly enriched and statistically significant. The presence of signaling-related pathways (e.g., AMPK signaling) further supports a multifaceted metabolic dysfunction in obese VAT. Together, these analyses suggest that the visceral adipose tissue of obese individuals exhibits extensive alterations in core metabolic pathways, nucleotide turnover, and hormonal biosynthesis, which may underlie systemic metabolic disturbances associated with obesity.

Figure 4
Diagram consisting of three panels: A shows a bar chart of KEGG pathway annotations divided into categories like metabolism, human diseases, and cellular processes. B is a dot plot showing the top 15 pathway enrichments with varying dot sizes and colors based on gene number and significance. C displays a circular diagram indicating the number of genes, selection count, and rich factors across various pathways, with color coding for different categories.

Figure 4. Functional annotation of adipose tissue metabolites in obesity vs. controls. (A) KEGG pathway classification of metabolites. (B) Bubble chart of top 15 enriched KEGG pathways. (C) Circular visualization of KEGG pathway enrichment (inner ring: numbers of KEGG terms, middle ring: −log10Pvalue, outer ring: KEGG pathway id).

Association analysis between representative metabolites, proteins, and clinical phenotypes

We performed an integrative correlation analysis between representative metabolites, proteins, and a panel of clinical phenotypes across VAT samples collected from bariatric surgery patients and healthy controls for investigating the biological relevance of the identified differential molecules. As shown in Figure 5A, a Pearson correlation heatmap was generated, where the columns represent selected proteins and metabolites with significant differential expression, and the rows correspond to key clinical indices, including BMI, blood lipid profiles (TG, HDL, LDL, TC), liver enzymes (ALT, AST), renal function markers (creatinine, urea nitrogen), glucose metabolism, and clinical histories (e.g., MASLD, diabetes) (Supplementary Table S7). The color scale reflects the direction and magnitude of correlation coefficients, while the asterisk indicates that the correlation coefficient is greater than 0.4. Notably, 4-Vinylcyclohexene showed a strong positive correlation with BMI (r = 0.85, P < 0.01), whereas asparagine-betaxanthin displayed a significant negative correlation with BMI (r = -0.76, P < 0.01). Other significant correlations included KRT1 with TG (r = 0.40, P < 0.05), NF1 with TG (r = -0.43, P < 0.05), ATR with HDL (r = 0.40, P < 0.05), and (R)-(+)-ropivacaine with HDL (r = -0.41, P < 0.05).

Figure 5
Heatmap and scatter plots showing relationships between various metabolic markers and compounds. Panel A displays a clustered heatmap with color gradients from purple to green, indicating correlation strengths among biomarkers like BMI, HDL, and various genes and compounds such as ATR and NF1. Asterisks denote significant correlations. Panel B contains six scatter plots with linear regression lines, depicting correlations between BMI, triglycerides, HDL, and compounds like 4-vinylcyclohexene and asparagine-betaxanthin, with correlation coefficients and P-values indicated.

Figure 5. Correlation analysis between molecules and clinical phenotypes. (A) Heatmap of spearman correlations between proteins/metabolites and clinical parameters. (B) Scatter plots of representative correlations. * means P<0.05.

To enhance the visualization of these findings, scatter plots of six representative molecule-phenotype pairs with moderate to strong correlations are presented in Figure 5B. These plots reveal the direction and strength of the associations, underscoring the potential physiological roles of these biomolecules in metabolic regulation. Overall, the results demonstrate that specific metabolites and proteins identified in adipose tissue are closely linked to obesity-related clinical characteristics. These molecular correlations may indicate underlying changes in adipose tissue function and provide candidate targets for further mechanistic studies or therapeutic interventions.

Integrated proteomic and metabolomic analysis reveals its core regulatory network

To explore the potential molecular interactions underlying metabolic alterations in visceral adipose tissue, we performed an integrative analysis of proteomic and metabolomic data based on Pearson correlation. Representative DEPs and metabolites were selected for pairwise correlation analysis across all samples. As shown in Figure 6A, we constructed a heatmap displaying all protein-metabolite pairs with correlation coefficients above 0.7, indicating strong associations. The heatmap highlights coordinated changes between specific proteins and metabolites, suggesting their potential co-regulation interplay in adipose tissue metabolism. Subsequently, we extracted these highly correlated pairs to construct a protein–metabolite interaction network, revealing a core regulatory module (Figure 6B). Within this network, several key proteins emerged as hub nodes, including OSBPL10 (Oxysterol binding protein like 10), CUL2 (Cullin 2), LRCH3 (Leucine rich repeats and calponin homology domain containing 3), and PRTN3 (Proteinase 3). These proteins were closely associated with multiple metabolites, which may imply that they play a central role in mediating the proteomic-metabolomic crosstalk in the VAT microenvironment. The identification of these protein hubs may reflect underlying pathways involved in dyslipidemia (OSBPL10), beige fat biogenesis (CUL2), fatty liver disease and insulin resistance (PRTN3) (3133). This integrated network provides novel insights into the molecular architecture governing adipose tissue dysfunction in obesity and may guide further exploration of regulatory mechanisms and therapeutic targets.

Figure 6
Panel A shows a heatmap with hierarchical clustering of gene expression data, displaying patterns in red and blue. Panel B is a network diagram illustrating interactions among molecules and genes with nodes and connecting lines, highlighting relationships.

Figure 6. Integrative analysis of proteins and metabolites. (A) Heatmap of pearson correlations between DEPs and differential metabolites (|cor| > 0.7). (B) Interaction network of strongly correlated protein-metabolite pairs (nodes: proteins/metabolites, node color: yellow represents positive correlation, and blue represents negative correlation). * means P<0.05, ** means P<0.01, *** means P<0.001.

Discussions

In this study, we observed significant alterations in both the proteomic and metabolomic profiles of VAT in obese patients undergoing SG compared to individuals with normal body weight. These molecular signatures reflect a profound remodeling of adipose tissue biology. Previous studies have reported that weight loss can reduce adipocyte size and reverse obesity-associated inflammation and insulin resistance (34, 35). Consistent with these findings, our results revealed changes in proteins involved in lipid handling (e.g., Lipid droplet organization pathway) and energy sensing (e.g., Muscle system process pathway). Lipid droplets are key organelles for triglyceride storage in adipocytes, impaired degradation of PLIN2 in obesity can lead to abnormally enlarged lipid droplets, contributing to lipotoxicity (36, 37). In addition, muscle contraction could regulate energy expenditure and adipogenesis by AMPK signaling (38, 39). Additionally, the protein autophosphorylation pathway was found to be suppressed in post-surgical patients. Dysregulated protein phosphorylation profile underlies muscle insulin resistance in type 2 diabetes (40). Collectively, compared to controls, the VAT of obese patients following SG exhibits distinct molecular signatures in terms of lipid droplet dynamics, inflammatory status, energy metabolism, and insulin sensitivity.

Among the pathway alterations, the most striking were those related to lipid metabolism signaling. Functional annotation analysis of DEPs and metabolites revealed significant disturbances in pathways governing lipid droplet biogenesis and turnover. Enlarged lipid droplets, as a hallmark of adipocyte hypertrophy, are one of the key contributor to obesity-induced insulin resistance (41). In fact, evidences have shown that weight loss can reduce adipose tissue mass, alleviate lipid overload, and partially reverse type 2 diabetes (42, 43). Our findings suggest that SG may induce a coordinated metabolic shift away from lipid accumulation, characterized by reduced lipid droplet formation and enhanced lipid mobilization or oxidation. This is further supported by the observed upregulation of the muscle contraction signaling pathway, through which skeletal muscle activity activates the AMPK/PGC-1α axis, promoting mitochondrial biogenesis, increasing fatty acid oxidation, and mitigating obesity-associated metabolic dysfunction (44).

Inflammatory activity also underwent marked changes following obesity. Many of the DEPs identified in our study are involved in innate immune responses and protein-metabolite interaction networks centered on inflammatory mediators. Notably, the neutrophil-derived serine protease PRTN3 (proteinase 3) emerged as a key hub. PRTN3 promotes adipose inflammation by activating proinflammatory cytokines such as IL-1β (45). Importantly, experimental models have shown that PRTN3 is upregulated in the steatotic liver of obese mice, while PRTN3 knockout prevents hepatic steatosis and insulin resistance (46). In light of our findings, the differential level of PRTN3 suggests that weight-loss interventions may suppress neutrophil-driven inflammation in visceral adipose tissue. This suppression may help restore adipokine homeostasis and improve systemic insulin sensitivity. In addition, we observed that several differentially expressed metabolites were enriched in purine and pyrimidine metabolism pathways, which are known to reflect shifts in inflammatory signaling (47). Collectively, the reduction in proinflammatory proteins and metabolites likely reflects an attenuation of chronic adipose immune stress following weight loss.

Our findings also highlight key regulatory factors associated with adipocyte phenotypes. For example, OSBPL10 was identified as a critical node. OSBPL10 binds sterols and has been implicated in dyslipidemia (48). Its presence in our network may suggest that alterations in cholesterol sensing accompany VAT remodeling. Similarly, CUL2 emerged as another central hub. Recent studies have shown that silencing CUL2 in adipocytes stabilizes PRDM16, thereby strongly promoting fatty acid oxidation while simultaneously suppressing pro-inflammatory and pro-fibrotic pathways (49). The expression level of CUL2 in VAT following SG may reflect the mechanistic basis of weight-loss efficacy. Clinically, these molecular features may possess predictive or therapeutic relevance. For example, given that genetic variants in OSBPL10 affect serum lipid profiles, patients who fail to downregulate OSBPL10 postoperatively may be predisposed to persistent dyslipidemia or delayed hepatic lipid clearance. Likewise, the role of CUL2 in thermogenic regulation suggests that elevated CUL2 expression in VAT could indicate attenuated weight loss or reduced metabolic benefit. Equally important, the core gene PRTN3 is directly linked to the pathophysiology of MASLD. If patients with SG who show greater remission of MASLD also exhibit a larger reduction in adipose PRTN3 levels, it may indicate a mechanistic link between adipose immune signaling and hepatic outcomes. As a result, these proteins and their associated metabolites represent promising candidate biomarkers for tracking and potentially predicting postoperative trajectories, including BMI changes, insulin sensitivity, and MASLD progression.

Our integrative proteomic-metabolomic analysis reveals key drivers of adipose remodeling after bariatric surgery. Central proteins such as OSBPL10, CUL2, and PRTN3 link lipid metabolism, energy balance, and immune signaling. While these findings provide mechanistic insights and reveal potential therapeutic targets, several limitations warrant consideration. The absence of longitudinal VAT sampling due to ethical constraints necessitates cautious interpretation of metabolic improvements, requiring further validation through animal models and in vitro systems. Notably, the dynamic nature of VAT remodeling may be influenced by interindividual variability and postoperative lifestyle modifications, underscoring the need for: (1) extended clinical follow-up with serial tissue profiling; (2) functional validation of candidate biomarkers using organoid models; and (3) rigorous preclinical evaluation of identified targets. These investigations will be critical for translating our discoveries into clinical treatment targets and biomarkers for obesity.

Conclusions

In summary, this study integrates proteomic and metabolomic analyses to elucidate the molecular remodeling of visceral adipose tissue (VAT) following sleeve gastrectomy (SG). Compared to normal-weight controls, obese patients undergoing SG exhibited significant alterations in protein and metabolite profiles within VAT. These changes predominantly involved coordinated regulation of key pathways, including lipid metabolism, inflammatory signaling, and insulin sensitivity. Notably, pathways related to lipid metabolism, lipid droplet organization, and protein autophosphorylation were markedly perturbed. Further analysis identified several hub proteins, such as OSBPL10, CUL2, and PRTN3, which may play pivotal roles in linking adipose tissue remodeling to postoperative clinical benefits, including improvements in MASLD. These hub proteins also hold potential as biomarkers for monitoring surgical outcomes. Collectively, the findings provide a theoretical framework and potential therapeutic targets for the development of non-surgical interventions that mimic the metabolic benefits of bariatric surgery. Future studies should explore how modulating these factors might replicate surgical benefits through less invasive strategies.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving humans were approved by The Third People’s Hospital of Chengdu ethical committee. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

YL: Writing – review & editing, Formal Analysis, Data curation, Investigation, Resources. HY: Investigation, Resources, Writing – review & editing, Data curation. XZ: Resources, Data curation, Writing – review & editing, Investigation. XHe: Investigation, Writing – review & editing. AL: Investigation, Writing – review & editing. RL: Investigation, Writing – review & editing. XHa: Investigation, Writing – review & editing. YML: Writing – review & editing, Supervision, Investigation, Resources. PG: Investigation, Writing – original draft, Visualization, Supervision, Formal Analysis, Methodology, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. The authors declare that the study and/or publication of this article is supported by the Fundamental Research Funds for the Central Universities (No: 2682025ZTPY057) and the Third People's Hospital of Chengdu Clinical Research Program (No: CSY-YN-01- 2023-054).

Acknowledgments

We sincerely thank all the volunteers who participated in this study. We also gratefully acknowledge the Longitudinal Study of Bariatric Surgery in Western China (registration number: ChiCTR2300073353) for providing comprehensive and high-quality clinical data, which were essential for the successful conduct of this study.

Conflict of interest

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

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

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/fendo.2025.1625501/full#supplementary-material

Supplementary Table 1 | Clinical phenotypic characteristics of obese patients and controls.

Supplementary Table 2 | The results of differential protein analysis.

Supplementary Table 3 | GO annotation analysis results of proteins.

Supplementary Table 4 | KEGG annotation analysis results of proteins.

Supplementary Table 5 | The results of differential metabolites analysis.

Supplementary Table 6 | KEGG annotation analysis results of metabolites.

Supplementary Table 7 | Correlation analysis results between key molecules and clinical information.

References

1. Phelps NH, Singleton RK, Zhou B, Heap RA, Mishra A, James EB, et al. Worldwide trends in underweight and obesity from 1990 to 2022: a pooled analysis of 3663 population-representative studies with 222 million children, adolescents, and adults. Lancet. (2024) 403:1027–50. doi: 10.1016/S0140-6736(23)02750-2

PubMed Abstract | Crossref Full Text | Google Scholar

2. Ruze R, Liu T, Zou X, Song J, Chen Y, Xu R, et al. Obesity and type 2 diabetes mellitus: connections in epidemiology, pathogenesis, and treatments. Front Endocrinol (Lausanne). (2023) 14:1161521. doi: 10.3389/fendo.2023.1161521

PubMed Abstract | Crossref Full Text | Google Scholar

3. Avgerinos KI, Spyrou N, Mantzoros CS, and Dalamaga M. Obesity and cancer risk: Emerging biological mechanisms and perspectives. Metabolism. (2019) 92:121–35. doi: 10.1016/j.metabol.2018.11.001

PubMed Abstract | Crossref Full Text | Google Scholar

4. Koskinas KC, Van Craenenbroeck EM, Antoniades C, Blüher M, Gorter TM, Henner H, et al. Obesity and cardiovascular disease: an ESC clinical consensus statement. Eur Heart J. (2024) 45:4063–98. doi: 10.1093/eurheartj/ehae508

PubMed Abstract | Crossref Full Text | Google Scholar

5. Salminen P, Helmiö M, Ovaska J, Juuti A, Leivonen M, Peromaa-Haavisto P, et al. Effect of Laparoscopic Sleeve Gastrectomy vs Laparoscopic Roux-en-Y Gastric Bypass on Weight Loss at 5 Years Among Patients With Morbid Obesity: The SLEEVEPASS Randomized Clinical Trial. JAMA. (2018) 319:241–54. doi: 10.1001/jama.2017.20313

PubMed Abstract | Crossref Full Text | Google Scholar

6. Svanevik M, Lorentzen J, Borgeraas H, Sandbu R, Seip B, Medhus AW, et al. Patient-reported outcomes, weight loss, and remission of type 2 diabetes 3 years after gastric bypass and sleeve gastrectomy (Oseberg); a single-centre, randomised controlled trial. Lancet Diabetes Endocrinol. (2023) 11:555–66. doi: 10.1016/S2213-8587(23)00127-4

PubMed Abstract | Crossref Full Text | Google Scholar

7. Han Y, Jia Y, Wang H, Cao L, and Zhao Y. Comparative analysis of weight loss and resolution of comorbidities between laparoscopic sleeve gastrectomy and Roux-en-Y gastric bypass: A systematic review and meta-analysis based on 18 studies. Int J Surg. (2020) 76:101–10. doi: 10.1016/j.ijsu.2020.02.035

PubMed Abstract | Crossref Full Text | Google Scholar

8. Biter LU, Hart JW, Noordman BJ, Smulders JF, Nienhuijs S, Dunkelgrün M, et al. Long-term effect of sleeve gastrectomy vs Roux-en-Y gastric bypass in people living with severe obesity: a phase III multicentre randomised controlled trial (SleeveBypass). Lancet Reg Health Eur. (2024) 38:100836. doi: 10.1016/j.lanepe.2024.100836

PubMed Abstract | Crossref Full Text | Google Scholar

9. Kawai T, Autieri MV, and Scalia R. Adipose tissue inflammation and metabolic dysfunction in obesity. Am J Physiol Cell Physiol. (2021) 320:C375–91. doi: 10.1152/ajpcell.00379.2020

PubMed Abstract | Crossref Full Text | Google Scholar

10. Zatterale F, Longo M, Naderi J, Raciti GA, Desiderio A, Miele C, et al. Chronic adipose tissue inflammation linking obesity to insulin resistance and type 2 diabetes. Front Physiol. (2020) 10:1607. doi: 10.3389/fphys.2019.01607

PubMed Abstract | Crossref Full Text | Google Scholar

11. Cho YK, Lee S, Lee J, Doh J, Park JH, Jung YS, et al. Lipid remodeling of adipose tissue in metabolic health and disease. Exp Mol Med. (2023) 55:1955–73. doi: 10.1038/s12276-023-01071-4

PubMed Abstract | Crossref Full Text | Google Scholar

12. Douros JD, Tong J, and D’Alessio DA. The effects of bariatric surgery on islet function, insulin secretion, and glucose control. Endocr Rev. (2019) 40:1394–423. doi: 10.1210/er.2018-00183

PubMed Abstract | Crossref Full Text | Google Scholar

13. Rodriguez-Muñoz A, Motahari-Rad H, Martin-Chaves L, Benitez-Porres J, Rodriguez-Capitan J, Gonzalez-Jimenez A, et al. A systematic review of proteomics in obesity: unpacking the molecular puzzle. Curr Obes Rep. (2024) 13:403–38. doi: 10.1007/s13679-024-00561-4

PubMed Abstract | Crossref Full Text | Google Scholar

14. Beyoğlu D, Popov YV, and Idle JR. Metabolomic hallmarks of obesity and metabolic dysfunction-associated steatotic liver disease. Int J Mol Sci. (2024) 25:12809. doi: 10.3390/ijms252312809

PubMed Abstract | Crossref Full Text | Google Scholar

15. Carlsson ER, Grundtvig JLG, Madsbad S, and Fenger M. Changes in serum sphingomyelin after Roux-en-Y gastric bypass surgery are related to diabetes status. Front Endocrinol (Lausanne). (2018) 9:172. doi: 10.3389/fendo.2018.00172

PubMed Abstract | Crossref Full Text | Google Scholar

16. Cox J and Mann M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotechnol. (2008) 26:1367–72. doi: 10.1038/nbt.1511

PubMed Abstract | Crossref Full Text | Google Scholar

17. UniProt Consortium. UniProt: the universal protein knowledgebase in 2021. Nucleic Acids Res. (2021) 49:D480–9. doi: 10.1093/nar/gkaa1100

PubMed Abstract | Crossref Full Text | Google Scholar

18. Adusumilli R and Mallick P. Data conversion with ProteoWizard msConvert. Methods Mol Biol. (2017) 1550:339–68. doi: 10.1007/978-1-4939-6747-6_23

PubMed Abstract | Crossref Full Text | Google Scholar

19. Domingo-Almenara X and Siuzdak G. Metabolomics data processing using XCMS. Methods Mol Biol. (2020) 2104:11–24. doi: 10.1007/978-1-0716-0239-3_2

PubMed Abstract | Crossref Full Text | Google Scholar

20. Chen Q, Liang X, Wu T, Jiang J, Jiang Y, Zhang S, et al. Integrative analysis of metabolomics and proteomics reveals amino acid metabolism disorder in sepsis. J Transl Med. (2022) 20:123. doi: 10.1186/s12967-022-03320-y

PubMed Abstract | Crossref Full Text | Google Scholar

21. Wishart DS, Guo A, Oler E, Wang F, Anjum A, Peters H, et al. HMDB 5.0: the human metabolome database for 2022. Nucleic Acids Res. (2022) 50:D622–31. doi: 10.1093/nar/gkab1062

PubMed Abstract | Crossref Full Text | Google Scholar

22. Kanehisa M, Furumichi M, Sato Y, Kawashima M, and Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. (2023) 51:D587–92. doi: 10.1093/nar/gkac963

PubMed Abstract | Crossref Full Text | Google Scholar

23. Wu T, Hu E, Xu S, Chen M, Guo P, Dai Z, et al. clusterProfiler 4.0: A universal enrichment tool for interpreting omics data. Innovation (Camb). (2021) 2:100141. doi: 10.1016/j.xinn.2021.100141

PubMed Abstract | Crossref Full Text | Google Scholar

24. Szklarczyk D, Gable AL, Nastou KC, Lyon D, Kirsch R, Pyysalo S, et al. The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. (2021) 49:D605–12. doi: 10.1093/nar/gkaa1074

PubMed Abstract | Crossref Full Text | Google Scholar

25. Doncheva NT, Morris JH, Gorodkin J, and Jensen LJ. Cytoscape StringApp: network analysis and visualization of proteomics data. J Proteome Res. (2019) 18:623–32. doi: 10.1021/acs.jproteome.8b00702

PubMed Abstract | Crossref Full Text | Google Scholar

26. Li Y, Guo W, Li H, Wang Y, Liu X, Kong W, et al. The change of skeletal muscle caused by inflammation in obesity as the key path to fibrosis: thoughts on mechanisms and intervention strategies. Biomolecules. (2024) 15:20. doi: 10.3390/biom15010020

PubMed Abstract | Crossref Full Text | Google Scholar

27. Engin A. Human protein kinases and obesity. Adv Exp Med Biol. (2017) 960:111–34. doi: 10.1007/978-3-319-48382-5_5

PubMed Abstract | Crossref Full Text | Google Scholar

28. AlMarzooqi SK, Almarzooqi F, Sadida HQ, Jerobin J, Ahmed I, Abou-Samra AB, et al. Deciphering the complex interplay of obesity, epithelial barrier dysfunction, and tight junction remodeling: Unraveling potential therapeutic avenues. Obes Rev. (2024) 25:e13766. doi: 10.1111/obr.13766

PubMed Abstract | Crossref Full Text | Google Scholar

29. Ignatieva EV, Afonnikov DA, Saik OV, Rogaev EI, and Kolchanov NA. A compendium of human genes regulating feeding behavior and body weight, its functional characterization and identification of GWAS genes involved in brain-specific PPI network. BMC Genet. (2016) 17:158. doi: 10.1186/s12863-016-0466-2

PubMed Abstract | Crossref Full Text | Google Scholar

30. Roa I and Del Sol M. Obesity, salivary glands and oral pathology. Colomb Med (Cali). (2018) 49:280–7. doi: 10.25100/cm.v49i3.3919

PubMed Abstract | Crossref Full Text | Google Scholar

31. Koriyama H, Nakagami H, Katsuya T, Akasaka H, Saitoh S, Shimamoto H, et al. Variation in OSBPL10 is associated with dyslipidemia. Hypertens Res. (2010) 33:511–4. doi: 10.1038/hr.2010.28

PubMed Abstract | Crossref Full Text | Google Scholar

32. Wang Q, Li H, Tajima K, Verkerke ARP, Taxin ZH, Hou Z, et al. Post-translational control of beige fat biogenesis by PRDM16 stabilization. Nature. (2022) 609:151–8. doi: 10.1038/s41586-022-05067-4

PubMed Abstract | Crossref Full Text | Google Scholar

33. Toonen EJ, Mirea AM, Tack CJ, Stienstra R, Ballak DB, Diepen JA, et al. Activation of proteinase 3 contributes to Non-alcoholic Fatty Liver Disease (NAFLD) and insulin resistance. Mol Med. (2016) 22:202–14. doi: 10.2119/molmed.2016.00033

PubMed Abstract | Crossref Full Text | Google Scholar

34. Pasarica M, Tchoukalova YD, Heilbronn LK, Fang X, Albu JB, Kelley DE, et al. Differential effect of weight loss on adipocyte size subfractions in patients with type 2 diabetes. Obes (Silver Spring). (2014) 17:1976–8. doi: 10.1038/oby.2009.219

PubMed Abstract | Crossref Full Text | Google Scholar

35. Villarreal-Calderón JR, Cuéllar RX, Ramos-González MR, Rubio-Infante N, Castillo EC, Elizondo-Montemayor L, et al. Interplay between the adaptive immune system and insulin resistance in weight loss induced by bariatric surgery. Oxid Med Cell Longev. (2019) 2019:3940739. doi: 10.1155/2019/3940739

PubMed Abstract | Crossref Full Text | Google Scholar

36. Cohen S. Lipid droplets as organelles. Int Rev Cell Mol Biol. (2018) 337:83–110. doi: 10.1016/bs.ircmb.2017.12.007

PubMed Abstract | Crossref Full Text | Google Scholar

37. Konige M, Wang H, and Sztalryd C. Role of adipose specific lipid droplet proteins in maintaining whole body energy homeostasis. Biochim Biophys Acta. (2014) 1842:393–401. doi: 10.1016/j.bbadis.2013.05.007

PubMed Abstract | Crossref Full Text | Google Scholar

38. Mengeste AM, Rustan AC, and Lund J. Skeletal muscle energy metabolism in obesity. Obes (Silver Spring). (2021) 29:1582–95. doi: 10.1002/oby.23227

PubMed Abstract | Crossref Full Text | Google Scholar

39. Kjøbsted R, Hingst JR, Fentz J, Foretz M, Sanz MN, Pehmoller C, et al. AMPK in skeletal muscle function and metabolism. FASEB J. (2018) 32:1741–77. doi: 10.1096/fj.201700442R

PubMed Abstract | Crossref Full Text | Google Scholar

40. Batista TM, Jayavelu AK, Wewer Albrechtsen NJ, Iovino S, Lebastchi J, Pan H, et al. A cell-autonomous signature of dysregulated protein phosphorylation underlies muscle insulin resistance in type 2 diabetes. Cell Metab. (2020) 32:844–859.e5. doi: 10.1016/j.cmet.2020.08.007

PubMed Abstract | Crossref Full Text | Google Scholar

41. Sanjabi B, Dashty M, Özcan B, Akbarkhanzadeh V, Rahimi M, Vinciguerra M, et al. Lipid droplets hypertrophy: a crucial determining factor in insulin regulation by adipocytes. Sci Rep. (2015) 5:8816. doi: 10.1038/srep08816

PubMed Abstract | Crossref Full Text | Google Scholar

42. Taylor R, Al-Mrabeh A, and Sattar N. Understanding the mechanisms of reversal of type 2 diabetes. Lancet Diabetes Endocrinol. (2019) 7:726–36. doi: 10.1016/S2213-8587(19)30076-2

PubMed Abstract | Crossref Full Text | Google Scholar

43. Magkos F, Fraterrigo G, Yoshino J, Luecking C, Kirbach K, Kelly SC, et al. Effects of moderate and subsequent progressive weight loss on metabolic function and adipose tissue biology in humans with obesity. Cell Metab. (2016) 23:591–601. doi: 10.1016/j.cmet.2016.02.005

PubMed Abstract | Crossref Full Text | Google Scholar

44. Lira VA, Benton CR, Yan Z, and Bonen A. PGC-1alpha regulation by exercise training and its influences on muscle function and insulin sensitivity. Am J Physiol Endocrinol Metab. (2010) 299:E145–161. doi: 10.1152/ajpendo.00755.2009

PubMed Abstract | Crossref Full Text | Google Scholar

45. Toonen EJM, Mirea AM, Tack CJ, Stienstra R, Ballak DB, Diepenet JAV, et al. Activation of proteinase 3 contributes to nonalcoholic fatty liver disease and insulin resistance. Mol Med. (2016) 22:202–14. doi: 10.2119/molmed.2016.00033

PubMed Abstract | Crossref Full Text | Google Scholar

46. Mirea AM, Stienstra R, Kanneganti TD, Tack CJ, Chavakis T, Toonen EJM, et al. Mice deficient in the IL-1β Activation genes Prtn3, Elane, and Casp1 are protected against the development of obesity-induced NAFLD. Inflammation. (2020) 43:1054–64. doi: 10.1007/s10753-020-01190-4

PubMed Abstract | Crossref Full Text | Google Scholar

47. Linden J, Koch-Nolte F, and Dahl G. Purine release, metabolism, and signaling in the inflammatory response. Annu Rev Immunol. (2019) 37:325–47. doi: 10.1146/annurev-immunol-051116-052406

PubMed Abstract | Crossref Full Text | Google Scholar

48. Perttilä J, Merikanto K, Naukkarinen J, Surakka I, Martin NW, Tanhuanpää K, et al. OSBPL10, a novel candidate gene for high triglyceride trait in dyslipidemic Finnish subjects, regulates cellular lipid metabolism. J Mol Med (Berl). (2009) 87:825–35. doi: 10.1007/s00109-009-0490-z

PubMed Abstract | Crossref Full Text | Google Scholar

49. Mao L, Lu J, Hou Y, and Nie T. Directly targeting PRDM16 in thermogenic adipose tissue to treat obesity and its related metabolic diseases. Front Endocrinol (Lausanne). (2024) 15:1458848. doi: 10.3389/fendo.2024.1458848

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: obesity, adipose tissue, signature, metabolomic, proteomic

Citation: Li Y, Yang H, Zhang X, He X, Liuli A, Li R, Han X, Li Y and Gao P (2025) Integrated proteomic and metabolomic profiling identifies distinct molecular signatures and metabolic pathways associated with obesity and potential targets for anti-obesity therapies. Front. Endocrinol. 16:1625501. doi: 10.3389/fendo.2025.1625501

Received: 09 May 2025; Accepted: 15 July 2025;
Published: 31 July 2025.

Edited by:

Alma Martelli, University of Pisa, Italy

Reviewed by:

George Grant, Independent Researcher, Aberdeen, United Kingdom
Anna Citko, Akademia Medyczna Nauk Stosowanych i Holistycznych, Poland

Copyright © 2025 Li, Yang, Zhang, He, Liuli, Li, Han, Li and Gao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yongmei Li, bHltemhhbmc3MEAxNjMuY29t; Pan Gao, Z3Bmb2N1c0AxNjMuY29t

These authors have contributed equally to this work

Disclaimer: 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.