Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Advertisement
Scientific Reports volume 12, Article number: 17113 (2022)
73
1
Metrics details
Obesity is a global epidemic elevating the risk of various metabolic disorders. As there is a lack of effective drugs to treat obesity, we combined bioinformatics and reverse network pharmacology in this study to identify effective herbs to treat obesity. We identified 1011 differentially expressed genes (DEGs) of adipose tissue after weight loss by analyzing five expression profiles (GSE103766, GSE35411, GSE112307, GSE43471, and GSE35710) from the Gene Expression Omnibus (GEO) database. We identified 27 hub genes from the protein–protein interaction (PPI) network by performing MCODE using the Search Tool for the Retrieval of Interacting Genes (STRING) database. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses revealed that these hub genes have roles in the extracellular matrix–receptor interaction, cholesterol metabolism, PI3K-Akt signaling pathway, etc. Ten herbs (Aloe, Portulacae Herba, Mori Follum, Silybum Marianum, Phyllanthi Fructus, Pollen Typhae, Ginkgo Semen, Leonuri Herba, Eriobotryae Folium, and Litseae Fructus) targeting the nine hub genes (COL1A1, MMP2, MMP9, SPP1, DNMT3B, MMP7, CETP, COL1A2, and MUC1) using six ingredients were identified as the key herbs. Quercetin and (-)-epigallocatechin-3-gallate were determined to be the key ingredients. Lastly, Ingredients-Targets, Herbs-Ingredients-Targets, and Herbs-Taste-Meridian Tropism networks were constructed using Cytoscape to elucidate this complex relationship. This study could help identify promising therapeutic targets and drugs to treat obesity.
Obesity is a global health problem without a definitive cure. Being overweight enhances the risk of chronic diseases, including type 2 diabetes, cardiovascular disease, cancer, and reproductive disorders1,2,3,4. Patients with moderate or severe obesity will also face higher all-factor mortality5. A study revealed that the prevalence of high body mass index (BMI) has increased globally over recent decades. To date, no country has been able to reduce its obesity epidemic by following evidence-based policies6,7,8. Although dietary and lifestyle modifications are considered primary solutions to treat obesity, recent studies recommend medications combined with lifestyle modifications to reduce weight among patients having BMI ≥ 27 kg/m2 and other associated diseases, or BMI ≥ 30 kg/m29. Most drugs used to manage obesity focus on appetite control. These drugs help in reducing food intake by stimulating pro-opiomelanocortin (POMC) neurons to promote satiation, ultimately leading to weight loss10. However, these current anti-obesity drugs are associated with side effects such as insomnia, dry mouth, constipation, adverse gastrointestinal reactions, and acute liver injury 11,12. Therefore, effective alternative strategies are urgently required.
Several studies focus on using Chinese herbs and other natural medications to treat obesity13 owing to the advantages of traditional Chinese medicine (TCM) in enhancing drug safety and reducing complications due to the homologous characteristics of medicine and food14. A systematic review shows that herbal medicines can effectively reduce weight, BMI, waist circumference, hip circumference, and body fat15. They exert anti-obesity effects by suppressing appetite and reducing energy intake, inhibiting pancreatic lipase activity and fat absorption, stimulating thermogenesis and energy expenditure, increasing lipolysis, and reducing lipogenesis15. Moreover, a randomized controlled trial showed that safflower oil intake could modify body weight and shape without lifestyle adjustments and also regulate blood pressure and insulin resistance compared with the placebo16. Drug discovery is associated with challenges such as long development cycles, high costs and turnover rates, and constantly changing regulatory requirements, which are not ideal for inventors17. As opposed to developing brand-new drugs, the extraction of natural ingredients from known herbs has various advantages including considerable saving of both time and economic costs. Therefore, it is necessary to identify novel and natural compounds that are suitable for the treatment of obesity.
Owing to the popularity of gene-microarray and RNA sequencing, integrated bioinformatics has been widely used to analyze high-throughput sequencing data of various diseases. By identifying differentially expressed genes (DEGs) and screening suitable hub genes, disease-related signaling pathways and their mechanisms can be analyzed, becoming the basis for predicting targeted drugs. In addition, network pharmacology has attracted considerable attention for the further exploration of the relationship between drugs, diseases, and targets by using various databases for analysis and simulation18. Many studies have clarified the targets and mechanisms of several Chinese medicines and their components in treating diseases using this technology. However, only a few have focused on mining new targeted drugs. In this study, we used an alternative strategy and reversed this process. Accordingly, the Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/) database and the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP, https://tcmspw.com/tcmsp.php) were used to identify the hub genes associated with obesity, explore the relevant mechanism through enrichment analysis, and predict the herbs and their ingredients with potential anti-obesity effect.
After standardizing gene sets (Fig. 1), 1011 DEGs (|logFC|> 1, p < 0.05) were screened out from GSE103766, GSE35411, GSE112307, GSE43471, and GSE35710 based on the above method. The results included 513 downregulated and 498 upregulated genes, as shown in the volcano plot (Fig. 2 and Supplementary Table S1). The abscissa in the volcano plot is log2 (fold change) value, and the ordinate is log10 (p-value).
Box-plots of the expression profiles after consolidation and standardization. The x-axis label represents the sample symbol and the y-axis label represents gene expression values. The black line in the box-plot represents the median value of gene expression. (a) Standardization of GSE43471, (b) Standardization of GSE35411, (c) Standardization of GSE103766, (d) Standardization of GSE35710, (e) Standardization of GSE112307.
Volcano plot to identify differentially expressed genes (DEGs). (a) GSE43471, (b) GSE35710, (c) GSE35411, (d) GSE103766, (e) GSE112307. The x-axis label represents fold changes and the y-axis label represents the p-values. Red dots represent the 498 upregulated genes and green dots represent the 513 downregulated genes.
As shown in Supplementary Fig. S1, the PPI network of DEGs, based on the Search Tool for the Retrieval of Interacting Genes (STRING) database, includes 584 nodes and 1417 edges. Using the MCODE plugin in Cytoscape software, the most significant modules (score = 6.667) were recognized from the PPI network as comprising 27 hub genes, including ACP5, CETP, COL1A1, COL1A2, CSF1, DNMT3B, EED, HIST1H2AI, HIST1H2BB, HIST1H2BD, HIST1H4B, HIST1H4H, HIST2H3C, HP, LCN2, LIPC, LPA, MMP2, MMP7, MMP9, MSR1, MUC1, PLA2G7, SPP1, THBS1, THBS2, and VLDLR (Table 1 and Fig. 3).
Subnetwork of 27 hub genes from the protein–protein interaction (PPI) network. Node size and temperature color reflect the degree of connectivity (bigger node represents a higher degree and smaller node represents a lower degree; red node represents a higher degree and yellow node represents a lower degree).
An enrichment analysis bubble chart was drawn under GO level 2 classifications using Omicshare tools (Fig. 4 and Supplementary Table S2). As shown in the figure, hub genes were significantly enriched in regulating plasma lipoprotein particle levels, lipid transport, extracellular matrix (ECM) organization, response to reactive oxygen species, and the oxygen-containing compound for biological process (BP). The hub genes were significantly enriched for cell composition (CC) in lipoprotein particles, extracellular regions, ECM, extracellular exosomes, and secretory granules. For molecular function (MF), the hub genes were significantly elevated in lipoprotein particle binding, glycosaminoglycan binding, ECM structural constituents, and peptidase activity.
Biological functions based on Gene Ontology (GO) analysis of obesity-related hub genes. Advanced bubble chart shows significance in GO enrichment items of hub genes in three functional groups: biological process (BP), cell composition (CC), and molecular function (MF). The x-axis label represents the gene ratio (Rich Factor) and the y-axis label represents GO terms.
KEGG pathway enrichment analysis showed that the hub genes were primarily enriched in ECM–receptor interaction, cholesterol metabolism, PI3K-Akt, IL-17, and TNF signaling pathways, endocrine resistance, and leukocyte transendothelial migration (Fig. 5 and Supplementary Table S3).
Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of hub genes. The x-axis label represents the gene ratio (Rich factor) and the y-axis label represents the pathway.
We converted 27 gene names of the hub genes into protein names that could be recognized through the TCMSP database using the Universal Protein Resource (Uniprot). Moreover, the hub genes can be input in the required format to identify potential herbs with anti-obesity effects from the TCMSP database. After excluding the genes that were not present in the databases or those that had no related ingredients, nine were screened for further research, namely, COL1A1, MMP2, MMP9, SPP1, DNMT3B, MMP7, CETP, COL1A2, and MUC1. These genes corresponded to 16 ingredients [(-)-epigallocatechin-3-gallate (EGCG), arachidonic acid, arctiin, baicalein, beta-carotene, capillarisin, deoxypodophyllotoxin, ellagic acid, fisetin, irisolidone, luteolin, matrine, nobiletin, quercetin, rutaecarpine, tanshinone IIa] showing adequate OB and DL values (OB ≥ 30%, DL ≥ 0.18) (Supplementary Table S4).
There were 254 herbs with active ingredients in the databases. The top 10 herbs were Aloe, Portulacae Herba, Mori Follum, Silybum Marianum, Phyllanthi Fructus, Pollen Typhae, Ginkgo Semen, Leonuri Herba, Eriobotryae Folium, and Litseae Fructus. These were associated with more DEGs (related genes = 6) and were, therefore, selected as crucial herbs in our study and annotated using Chinese pharmaceutical properties (CMPs), including characters, tastes, and meridian tropisms (Table 2).
We screened the key ingredients in treating obesity using an Ingredients-Targets network containing 25 nodes and 27 edges (Fig. 6). The nine orange nodes represent the target genes and 16 green nodes represent the active ingredients. As most genes could be linked (degree = 4), quercetin and EGCG were considered the most critical components in the treatment of obesity.
Ingredients-Targets network. Nine orange nodes represent the target genes, whereas the 16 green nodes represent the active compounds. The edges represent the interaction between the compounds and targets.
As shown in Fig. 7a, the Herbs-Ingredients-Targets network containing 24 nodes and 43 edges was constructed to demonstrate the relationship between them: the 10 green nodes represent the key herbs and the six yellow nodes represent the active ingredients in them; the eight blue nodes depict the target genes. By analyzing the network, Phyllanthi Fructus and Portulacae Herba were associated with the most ingredients (degree = 4). Moreover, quercetin was the most frequent active ingredient (degree = 23) found in all herbs. Regarding gene targets, MMP2 was targeted by most ingredients (degree = 5) followed by MMP9 (degree = 4). Other genes were only acted upon by one component (degree = 1).
Herbs-Ingredients-Targets network (a) and Herbs-Taste-Meridian tropism (b) network. (a) Yellow nodes represent the active ingredients and the blue nodes represent the target genes. (b) Yellow nodes represent tastes and purple nodes represent meridian tropisms. In all networks, the light green nodes represent cold-cool herbs, medium green nodes represent calm herbs, and dark green nodes represent warm herbs.
We also established the Herbs-Taste-Meridian tropism network containing 24 nodes and 40 edges to clarify the distribution of CMPs (Fig. 7b). Five yellow nodes represent tastes and eight purple nodes represent meridian tropisms. To indicate different characters, we presented 10 nodes of herbs having different greens (light green, medium green, and dark green). Regarding characters, cold-cool herbs like Mori Follum were the most frequent (nodes = 7), followed by herbs having calm (nodes = 2) and warm (nodes = 1) characters. In terms of taste, herbs were mostly bitter (edges = 6), followed by sweet (edges = 4), acid (edges = 2), symplectic (edges = 2), and astringent (edges = 2). Regarding meridian tropism, most herbs belonged to the liver meridian (edges = 6), followed by the stomach and lung (edges = 4), large intestine (edges = 2), bladder (edges = 2), kidney (edges = 2), pericardium (edges = 2), spleen (edges = 1), and gallbladder (edges = 1) meridians.
Obesity causes energy imbalance due to excess energy intake (intake > expenditure). Controlling energy intake and expenditure is the key to maintaining metabolism and weight loss12. Most anti-obesity drugs lead to weight loss by promoting satiety and reducing food intake and calorie absorption. The TCM theory has different views. Although there is no objection with respect to the influence of energy balance on body weight depending on the law of conservation of energy, TCM theory expounds on obesity from the perspective of body functions other than the dynamic equilibrium of energy and nutrition. Based on TCM, obesity, as a physical state, results from the combined influence of various congenital constitutions and living habits. The obese population can be divided into TCM syndrome types based on their body function status rather than the degree of obesity depending on BMI, abdominal circumference, and other data17. The TCM theory consistently implements this holistic dialectical thought to discuss several biomedical problems and provides new treatments for various diseases18. Owing to the side effects of the current drugs used for weight loss on one hand and the safety of natural components on the other, it is important to predict and further explore potential TCM-based therapies.
Heredity is an important factor contributing to obesity. About 40%–70% of BMI differences among individuals are attributed to genetic differences19. Data from genome-wide association studies indicate more than 140 genetic chromosomal regions associated with obesity, indicating genetic involvement20. About 60% of the cases of obesity are due to genetic differences21. Claussnitzer et al.22 observed that regulating the expression of related genes in adipose tissue could reduce weight and enhance energy consumption in mice without altering their physical activity or appetite. The mechanism is associated with inhibiting adipocyte precursor mitochondria thermogenesis via the obesity-related FTO allele in a tissue-autonomous way. A meta-analysis19 revealed that 97 genetic loci associated with BMI comprised less than 3% (about 2.7%) of BMI variation. In contrast, 21% of BMI variation was likely due to common genetic variation, indicating that any single gene could not explicitly explain BMI variation. Microarray analysis in systems biology is an effective tool to explore the potential genes and pathways in different diseases. It is performed by detecting tens of thousands of gene expression information using a high-throughput pattern23. The GEO database is the most comprehensive public functional genomics database that includes gene expression, mutation, modification, and other related information24. Currently, the GEO database has 65,775 microarray expression profiling by array that has been queried and downloaded by users for free. Therefore, this study was conducted based on the appropriate microarray data from the GEO database. Additionally, the expression data of the same population before and after weight loss was selected from this database in the current study, rather than data from individuals who were obese or of normal weight, to avoid the influence of individual genetic factors during the acquisition and adjustment of obesity.
GO and KEGG enrichment analysis of hub genes revealed that most genes were enriched in ECM-reception interaction, including six DEGs. The ECM is a complex structure with different proteins, proteoglycans, and polysaccharides. It provides a scaffold for cells to control biological processes such as cell adhesion, migration, repair, survival, and development. Obesity is characterized by the extensive expansion of adipose tissue. Thus, ECM remodeling and recombination provide adequate space for adipocyte proliferation and the formation of new fat cells by the adipogenesis of precursor cells25. In addition, this process forms new blood vessels, which are essential for the expansion of healthy adipose tissue. In addition, adipose tissues require a vascular network to receive oxygen and nutrients26. Failure of this process will lead to adipocyte necrosis and hypoxia, which can cause chronic inflammation and fibrosis and induce other chronic diseases. Therefore, the growth and expansion of adipose tissue depend on angiogenesis. Intervention with angiogenesis inhibitors among obesity models can reduce the number of blood vessels in adipose tissues27 and reduce weight without impacting food intake28.
In this study, we constructed an Ingredients-Targets network using Cytoscape. Next, we screened out two compounds as crucial ingredients in treating obesity (quercetin and EGCG), which could each target four hub genes considered key in treating obesity. Quercetin is a plant-derived flavonoid with anti-inflammatory and anti-obesity effects29. Related studies indicate that quercetin also has antioxidant and cardioprotective effects. Its anti-obesity effect is due to the inhibition of adipogenesis by inhibiting the gene expression of peroxisome proliferator-activated receptor-γ (PPAR-γ)30. Quercetin is reported to significantly alleviate high-fat-diet–induced obesity and influence the production of metabolites linked to obesity-related inflammation and oxidative stress by regulating gut microbiota and metabolites31. EGCG is a bioactive polyphenol with antioxidant and anti-inflammatory effects32. Treatment with EGCG leads to reduced visceral adiposity and loss in body weight and is associated with the regulation of Beclin1-dependent autophagy in white adipose tissue33. Moreover, EGCG can inhibit adipocyte differentiation and maturation via the PI3K-AKT-FOXO1 pathway34.
We identified 10 key herbs in treating obesity, many of which have been widely used clinically for weight loss. Eriobotryae Folium can alleviate visceral and central obesity, insulin resistance, dyslipidemia, oxidative stress, and inflammation of metabolic syndrome in high-fat-diet models35. Aloe can inhibit obesity-induced inflammation by activating AMPK expression in muscle and reducing proinflammatory cytokines in white adipose tissue36. Silybum Marianum has lipid-lowering, antihypertensive, antidiabetic, anti-atherosclerotic, hepatoprotective, and anti-obesity effects and is known for its antioxidant, anti-inflammatory, and β-cell regenerative effect; enhancement of insulin sensitivity; and inhibition of gluconeogenesis and Glut4-mediated transport37. However, further investigations, such as randomized controlled trials, are required to verify the clinical effects and elucidate its anti-obesity mechanisms.
One of the innovations of this study was the construction of the Herbs-Taste-Meridian Tropism network and understanding the CMP rules of herbs with promising effects. Although several scholars have begun to study the weight-loss effect of Chinese herbs, research on potential herbs based on CMP is relatively insufficient. CMP is an integral aspect of the TCM theory38 that has guided the clinical practice of TCM practitioners for thousands of years as the primary prescription basis. Some research groups have studied Hedysarum Multijugum Maxim, Coicis Semen, and other Chinese herbs based on the theories of Yin-Yang and Chinese herbal characters. They found that warm-hot herbs could increase cAMP, prostaglandin E2, and other heat-related metabolites in rats with elevated body temperature. Cold-cool herbs could elevate metabolites such as histidine, tyrosine, lipid, and inositol and exert a cooling effect39. A study indicated that40 the characters, tastes, and meridian tropisms of Chinese herbs correlate with human tissues and organs, depicting a nonlinear relationship. Warm herbs have evident effects on the upper part of the liver, heart, and spinal cord, whereas hot herbs affect the stomach, kidneys, and small intestine. According to the theory of TCM, Schisandrae Chinensis Fructus tastes like acid, and acid herbs correspond to the liver. Therefore, Schisandrae Chinensis Fructus is traditionally used in a clinical setting as a liver tonic. Recent studies have revealed that it can alleviate the symptoms of metabolic syndrome by inhibiting oleic acid–induced fatty liver41. From the perspective of meridian tropism, Sanjiao Meridian has a significant impact on the thymus, brain, upper spine, and medulla. In contrast, the Pericardial Meridian has significant effects on the stomach, kidneys, small intestine, and medulla. The above studies could not completely summarize the relationship between CMP and disease targets and the therapeutic mechanism. However, they did explain its regularity and objectivity to a certain extent.
Discrepancies in gender, body weight, age, treatments, platforms, and other factors may have influenced the identification of hub genes, as we used five distinct datasets from the GEO database. Thus, further studies on these genes are required for validation. Next, the two key ingredients and 10 key herbs predicted in this study should be verified using clinical and animal experiments to confirm their effectiveness and heterogeneity in obesity therapy. Lastly, further studies are required to validate the association between drug properties and CMPs and to elucidate potential mechanisms of action.
In this study, we identified 27 hub genes associated with weight loss and analyzed their related signaling pathways using publicly available databases. Additionally, we predicted several potential natural compounds and Chinese herbs targeting these genes using the TCMSP platform. Quercetin and EGCG were screened out as crucial ingredients. Aloe, Portulacae Herba, Mori Follum, Silybum Marianum, Phyllanthi Fructus, Pollen Typhae, Ginkgo Semen, Leonuri Herba, Eriobotryae Folium, and Litseae Fructus were screened out as key herbs. These herbs and compounds are expected to be effective in treating obesity. Supplementary investigations involving randomized controlled trials and molecular biology studies are warranted to verify these findings.
The search strategy included (“obesity”[MeSH terms] OR obesity [All Fields]) AND “Homo sapiens”[porgn] AND “Expression profiling by array”[Filter]. The inclusion criteria were as follows: (i) subcutaneous adipose tissue (SAT) from patients who were obese after lifestyle-based weight-loss interventions; (ii) SAT of patients before weight-loss interventions were considered as controls.
Based on the above retrieval strategy, the five gene expression profiles that were selected for analysis from the GEO database were GSE10376642, GSE3541143, GSE11230744, GSE4347145, and GSE3571046, and their platforms were GPL19109, GPL10335, GPL6947, GPL6947 (same as previous), and GPL96. After the five data sets were retrieved and standardized using the GEOquery47 package in R software (version 4.1.3)48, ggplot2 package was used to design box-plots for gene expression data49.
DEGs with a threshold criterion of |log FC |> 1 and p-value < 0.0550 in SAT samples were screened using the linear models for microarray data package after incorporating weight-loss intervention51. Pheat map, ggplot2, and RColorBrewer packages52,53 were used to formulate the volcano plots of the DEGs.
After summarizing and removing duplicates DEGs, the PPI network of DEGs was determined to identify the most significant clusters of DEGs using STRING (version 11.5; https://cn.string-db.org/) having a combined score > 0.4 as the cutoff point. After STRING analysis, Cytoscape (version 3.9.1), an open-source bioinformatics software platform used for network analysis to visualize their associations, was used to visualize and identify the PPI network. MCODE plugin (version 2.0.0) was used to identify the hub genes, and the parameters of DEG clustering and scoring were as follows: MCODE score = 6.667, Degree Cut-off = 2, Node Score Cutoff = 0.2, k-score = 2, and Max. Depth = 100.
GO enrichment analysis annotates the genetic information based on three aspects: MF, BP, and CC. It is widely used in mining biological function information and the corresponding biological mechanism in microarray results54. KEGG (https://www.genome.jp/kegg/) is a comprehensive database involving pathways, genes, compounds, drugs, and diseases for annotating the biological functions of genes and genomes at the molecular level55. GO and KEGG enrichment analyses were performed on the screened hub genes using Omicshare tools, an online platform for data analysis and gene annotation of gene DENOVO (https://www.omicshare.com/tools/).
When retrieving drug information associated with specific genes (targets) from the TCMSP database, there are requirements on the format of input information. TCMSP can only identify protein names. Nevertheless, the gene expression matrices obtained from the previous steps had been written by the gene names. We needed to transform gene names into protein names that TCMSP could recognize to successfully compare the information of gene targets with Chinese herbs in TCMSP. This step was completed using UniProt (https://www.uniprot.org/help/about), a free comprehensive online database for protein sequences and their annotations56. After entering all protein names for the hub genes, we selected the active ingredients of Chinese herbs based on their oral bioavailability (OB) and drug-likeness (DL) values (OB ≥ 30%, DL ≥ 0.18).
After screening the herbs containing active ingredients that could treat obesity, the relationship between the herbs and target genes was studied. If different components in a particular herb were associated with the same target, only one association was considered between the herb and the target gene. Herbs associated with more target genes were considered more significant in treating obesity. We defined a part of the herbs acting on more target genes as the key herbs in treating obesity. Then, referring to the Chinese Pharmacopoeia 2020 Edition (I), key herbs were annotated with CMPs, such as characters, tastes, and meridian tropisms.
We established the Ingredients-Targets, Herbs-Ingredients-Targets, and Herbs-Taste-Meridian Tropism networks using Cytoscape to determine the association between the above nodes. The importance of each node in the network was determined based on the degree value of the network topology parameter. Thus, the importance of the node changed significantly with the increase in degree value, and the degree value was represented by the edge connecting the node57.
The data that support the findings of this study are available in GEO database (http://www.ncbi.nlm.nih.gov/geo), reference number [GSE103766, GSE35411, GSE112307, GSE43471, and GSE35710].
Tay, J. et al. Effects of an energy-restricted low-carbohydrate, high unsaturated fat/low saturated fat diet versus a high-carbohydrate, low-fat diet in type 2 diabetes: A 2-year randomized clinical trial. Diabetes Obes. Metab. 20(4), 858–871 (2018).
CAS PubMed Article Google Scholar
Hu, L. et al. Prevalence of overweight, obesity, abdominal obesity and obesity-related risk factors in southern China. PLoS ONE 12(9), e0183934 (2017).
PubMed PubMed Central Article Google Scholar
Ringel, A. E. et al. Obesity shapes metabolism in the tumor microenvironment to suppress anti-tumor immunity. Cell 183(7), 1848–1866 (2020).
CAS PubMed PubMed Central Article Google Scholar
Lainez, N. M. et al. Diet-induced obesity elicits macrophage infiltration and reduction in spine density in the hypothalami of male but not female mice. Front. Immunol. 9, 1992 (2018).
PubMed PubMed Central Article Google Scholar
Flegal, K. M. et al. Association of all-cause mortality with overweight and obesity using standard body mass index categories: a systematic review and meta-analysis. JAMA 309(1), 71–82 (2013).
CAS PubMed PubMed Central Article Google Scholar
Roberto, C. A. et al. Patchy progress on obesity prevention: Emerging examples, entrenched barriers, and new thinking. Lancet 385(9985), 2400–2409 (2015).
PubMed Article Google Scholar
Dai, H. et al. The global burden of disease attributable to high body mass index in 195 countries and territories, 1990–2017: An analysis of the global burden of disease study. PLoS Med 17(7), e1003198 (2020).
PubMed PubMed Central Article Google Scholar
Collaborators, G. B. D. O. et al. Health Effects of Overweight and Obesity in 195 Countries over 25 Years. N. Engl. J. Med. 377(1), 13–27 (2017).
Article Google Scholar
Jensen, M. D. et al. 2013 AHA/ACC/TOS guideline for the management of overweight and obesity in adults: A report of the American College of Cardiology/American heart association task force on practice guidelines and the obesity society. J. Am. Coll. Cardiol. 63(25 Pt B), 2985–3023 (2014).
PubMed Article Google Scholar
Keith, J. N. Pharmacotherapy in treatment of obesity. Gastroenterol. Clin. North Am. 45(4), 663–672 (2016).
PubMed Article Google Scholar
Singh, A. K. & Singh, R. Pharmacotherapy in obesity: A systematic review and meta-analysis of randomized controlled trials of anti-obesity drugs. Expert Rev. Clin. Pharmacol. 13(1), 53–64 (2020).
CAS PubMed Article Google Scholar
Gadde, K. M. et al. Obesity: Pathophysiology and management. J. Am. Coll. Cardiol. 71(1), 69–84 (2018).
PubMed PubMed Central Article Google Scholar
Sun, N. N., Wu, T. Y. & Chau, C. F. Natural dietary and herbal products in anti-obesity treatment. Molecules 21(10), 1351 (2016).
PubMed Central Article Google Scholar
Xu, L. et al. Chinese medicine in the battle against obesity and metabolic diseases. Front. Physiol. 9, 850 (2018).
PubMed PubMed Central Article Google Scholar
Payab, M. et al. Effect of the herbal medicines in obesity and metabolic syndrome: A systematic review and meta-analysis of clinical trials. Phytother. Res. 34(3), 526–545 (2020).
CAS PubMed Article Google Scholar
Benjamini, Y. et al. Controlling the false discovery rate in behavior genetics research. Behav. Brain Res. 125(1–2), 279–284 (2001).
CAS PubMed Article Google Scholar
Li, C., Zhang, H. & Li, X. The mechanism of traditional chinese medicine for the treatment of obesity. Diabetes Metab. Syndr. Obes. 13, 3371–3381 (2020).
PubMed PubMed Central Article Google Scholar
Yang, Y. et al. Traditional Chinese medicine in the treatment of patients infected with 2019-new coronavirus (SARS-CoV-2): A review and perspective. Int. J. Biol. Sci. 16(10), 1708–1717 (2020).
CAS PubMed PubMed Central Article Google Scholar
Locke, A. E. et al. Genetic studies of body mass index yield new insights for obesity biology. Nature 518(7538), 197–206 (2015).
CAS PubMed PubMed Central Article Google Scholar
Fall, T., Mendelson, M. & Speliotes, E. K. Recent advances in human genetics and epigenetics of adiposity: Pathway to precision medicine?. Gastroenterology 152(7), 1695–1706 (2017).
CAS PubMed Article Google Scholar
Elks, C. E. et al. Variability in the heritability of body mass index: A systematic review and meta-regression. Front. Endocrinol. (Lausanne) 3, 29 (2012).
Article Google Scholar
Claussnitzer, M. et al. FTO obesity variant circuitry and adipocyte browning in humans. N. Engl. J. Med. 373(10), 895–907 (2015).
CAS PubMed PubMed Central Article Google Scholar
Tao, J. et al. Critical roles of ELVOL4 and IL-33 in the progression of obesity-related cardiomyopathy via integrated bioinformatics analysis. Front. Physiol. 11, 542 (2020).
PubMed PubMed Central Article Google Scholar
Song, D. et al. Identification of biomarkers associated with diagnosis and prognosis of gastroesophageal junction adenocarcinoma-a study based on integrated bioinformatics analysis in GEO and TCGA database. Med. (Baltim.) 99(51), e23605 (2020).
CAS Article Google Scholar
Schoettl, T., Fischer, I. P. & Ussar, S. Heterogeneity of adipose tissue in development and metabolic function. J. Exp. Biol. 221(Pt Suppl 1), 162958 (2018).
Article Google Scholar
Cao, Y. Angiogenesis modulates adipogenesis and obesity. J. Clin. Invest. 117(9), 2362–2368 (2007).
CAS PubMed PubMed Central Article Google Scholar
Brakenhielm, E. et al. Angiogenesis inhibitor, TNP-470, prevents diet-induced and genetic obesity in mice. Circ. Res. 94(12), 1579–1588 (2004).
CAS PubMed Article Google Scholar
Rupnick, M. A. et al. Adipose tissue mass can be regulated through the vasculature. Proc. Natl. Acad. Sci. U S A 99(16), 10730–10735 (2002).
ADS CAS PubMed PubMed Central Article Google Scholar
Ohmae, S. et al. Quercetin attenuates adipogenesis and fibrosis in human skeletal muscle. Biochem. Biophys. Res. Commun. 615, 24–30 (2022).
CAS PubMed Article Google Scholar
Ulusoy, H. G. & Sanlier, N. A minireview of quercetin: From its metabolism to possible mechanisms of its biological activities. Crit. Rev. Food Sci. Nutr 60(19), 3290–3303 (2020).
CAS PubMed Article Google Scholar
Su, L. et al. Quercetin improves high-fat diet-induced obesity by modulating gut microbiota and metabolites in C57BL/6J mice. Phytother. Res. (2022).
Du, W. Q. & Xue, T. Progress of researches on anti-infective properties of epigallocatechin-3-gallate. Zhongguo Xue Xi Chong Bing Fang Zhi Za Zhi 34(1), 102–108 (2021).
CAS PubMed Google Scholar
Choi, C. et al. Epigallocatechin-3-gallate reduces visceral adiposity partly through the regulation of beclin1-dependent autophagy in white adipose tissues. Nutrients 12(10), 3072 (2020).
CAS PubMed Central Article Google Scholar
Lu, Y. et al. Epigallocatechin-3-gallate suppresses differentiation of adipocytes via regulating the phosphorylation of FOXO1 mediated by PI3K-AKT signaling in 3T3-L1 cells. Oncotarget 9(7), 7411–7423 (2018).
PubMed Article Google Scholar
Li, F. et al. Eriobotrya japonica leaf triterpenoid acids ameliorate metabolic syndrome in C57BL/6J mice fed with high-fat diet. Biomed. Pharmacother. 132, 110866 (2020).
CAS PubMed Article Google Scholar
Shin, E. et al. Dietary aloe reduces adipogenesis via the activation of AMPK and suppresses obesity-related inflammation in obese mice. Immune Netw. 11(2), 107–113 (2011).
PubMed PubMed Central Article Google Scholar
Tajmohammadi, A., Razavi, B. M. & Hosseinzadeh, H. Silybum marianum (milk thistle) and its main constituent, silymarin, as a potential therapeutic plant in metabolic syndrome: A review. Phytother. Res. 32(10), 1933–1949 (2018).
CAS PubMed Article Google Scholar
Wang, X. R. et al. Quantification of “cold-hot” medicinal properties of chinese medicines based on primary metabolites and Fisher’s analysis. Comput. Math. Methods Med. 2022, 5790893 (2022).
PubMed PubMed Central Google Scholar
Huang, Y. et al. The Chinese medicinal herbs of spleen-meridian property regulate body temperature in yeast-induced fever rats. Phytomedicine 74, 152815 (2020).
CAS PubMed Article Google Scholar
Yang, X. et al. Information integration research on cumulative effect of “Siqi, Wuwei, and Guijing” in Traditional Chinese Medicine. J. Tradit. Chin. Med. 36(4), 538–546 (2016).
PubMed Article Google Scholar
Wat, E. et al. An in vitro and in vivo study of a 4-herb formula on the management of diet-induced metabolic syndrome. Phytomedicine 42, 112–125 (2018).
PubMed Article Google Scholar
Bollepalli, S. et al. Subcutaneous adipose tissue gene expression and DNA methylation respond to both short- and long-term weight loss. Int. J. Obes. (Lond.) 42(3), 412–423 (2018).
CAS Article Google Scholar
Johansson, L. E. et al. Differential gene expression in adipose tissue from obese human subjects during weight loss and weight maintenance. Am. J. Clin. Nutr. 96(1), 196–207 (2012).
CAS PubMed Article Google Scholar
Dao, M. C. et al. A data integration multi-omics approach to study calorie restriction-induced changes in insulin sensitivity. Front. Physiol. 9, 1958 (2018).
PubMed Article Google Scholar
Campbell, K. L. et al. Gene expression changes in adipose tissue with diet- and/or exercise-induced weight loss. Cancer Prev. Res. (Phila.) 6(3), 217–231 (2013).
CAS Article Google Scholar
Nookaew, I. et al. Adipose tissue resting energy expenditure and expression of genes involved in mitochondrial function are higher in women than in men. J. Clin. Endocrinol. Metab. 98(2), E370–E378 (2013).
CAS PubMed Article Google Scholar
Davis, S. & Meltzer, P. S. GEOquery: A bridge between the gene expression omnibus (GEO) and BioConductor. Bioinformatics 23(14), 1846–1847 (2007).
PubMed Article Google Scholar
Zhao, C. Identifying the hub gene and immune infiltration of osteoarthritis by bioinformatical methods. Clin. Rheumato.l 40(3), 1027–1037 (2021).
Article Google Scholar
Maag, J. L. V. gganatogram: An R package for modular visualisation of anatograms and tissues based on ggplot2. F1000Res 7, 1576 (2018).
PubMed PubMed Central Article Google Scholar
Liu, Y. et al. Integrative analyses of biomarkers and pathways for adipose tissue after bariatric surgery. Adipocyte 9(1), 384–400 (2020).
CAS PubMed PubMed Central Article Google Scholar
Ritchie, M. E. et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucl. Acids Res. 43(7), e47 (2015).
PubMed PubMed Central Article Google Scholar
Khomtchouk, B. B., Van Booven, D. J. & Wahlestedt, C. HeatmapGenerator: high performance RNAseq and microarray visualization software suite to examine differential gene expression levels using an R and C++ hybrid computational pipeline. Source Code Biol. Med. 9(1), 30 (2014).
PubMed PubMed Central Article Google Scholar
Ramasamy, R. K. et al. STRUCTURE PLOT: A program for drawing elegant STRUCTURE bar plots in user friendly interface. Springerplus 3, 431 (2014).
PubMed PubMed Central Article Google Scholar
Gene Ontology, C. The gene ontology (GO) project in. Nucl. Acids Res. 34(Database issue), D322–D326 (2006).
Article Google Scholar
Kanehisa, M. et al. KEGG: New perspectives on genomes, pathways, diseases and drugs. Nucl. Acids Res. 45(D1), D353–D361 (2017).
CAS PubMed Article Google Scholar
UniProt, C. UniProt: The universal protein knowledgebase in 2021. Nucl. Acids Res. 49(D1), D480–D489 (2021).
Article Google Scholar
Wang, M. L. et al. Network pharmacology-based approach uncovers the mechanism of GuanXinNing tablet for treating thrombus by MAPKs signal pathway. Front. Pharmacol. 11, 652 (2020).
CAS PubMed PubMed Central Article Google Scholar
Download references
This research was supported by the National Natural Science Foundation of China (82174525) and Jilin Natural Science Foundation (YDZJ202201ZYTS195).
Changchun University of Chinese Medicine, Changchun, 130117, China
Yuxing Tai, Hongying Tian & Tianjiao Gao
Jilin Ginseng Academy, Changchun University of Chinese Medicine, Changchun, 130117, China
Xiaoqian Yang
Beijing University of Chinese Medicine, Beijing, 100029, China
Shixing Feng
Department of Acupuncture and Tuina, Changchun University of Chinese Medicine, Changchun, 130117, China
Shaotao Chen, Chongwen Zhong, Xiaochao Gang & Mingjun Liu
Acupuncture and Massage Center of the Third Affiliated Clinical Hospital, Changchun University of Chinese Medicine, Changchun, 130117, China
Shaotao Chen & Mingjun Liu
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
Y.T. and H.T.: Study Design, Writing – original Draft and Statistical Analysis. X.Y. and S.F.: Software, Visualization. S.C.: Manuscript Revision. C.Z., T.G., and X.G.: Data Collection, Data Interpretation. M.L.: Project Administration, Investigation and Manuscript Revision. All authors reviewed the manuscript.
Correspondence to Mingjun Liu.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Tai, Y., Tian, H., Yang, X. et al. Identification of hub genes and candidate herbal treatment in obesity through integrated bioinformatic analysis and reverse network pharmacology. Sci Rep 12, 17113 (2022). https://doi.org/10.1038/s41598-022-22112-4
Download citation
Received: 21 April 2022
Accepted: 10 October 2022
Published: 12 October 2022
DOI: https://doi.org/10.1038/s41598-022-22112-4
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Advertisement
© 2022 Springer Nature Limited
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.