Preprint
Essay

Pueraria thomsonii Radix Water Extract Alleviates Type 2 Diabetes Mellitus in db/db Mice through the Integrated Regulation of Metabolism and Gut Microbiota

Altmetrics

Downloads

221

Views

64

Comments

0

  † Jiarong Li and Yuhui Liu have contributed equally to this work.

This version is not peer-reviewed

Submitted:

05 March 2023

Posted:

06 March 2023

You are already at the latest version

Alerts
Abstract
Type 2 diabetes mellitus(T2DM) is an increasingly prevalent and serious health problem. Its onset is typically associated with metabolic disorders and disturbances in the gut microbiota. Previous studies have reported the anti-T2DM effects of Pueraria thomsonii Radix as a functional food. However, the mechanism of action is still unknown. In this study, the effects of Pueraria thomsonii Radix water extract(PTR) on db/db mice were evaluated by pharmacology, metabolomics, and 16S rRNA gene sequencing. The results showed that PTR could improve the body weight of mice, reduce fasting blood glucose(FBG), urinary glucose(UGLU), homeostasis model assessment insulin resistance(HOMA-IR), urinary albumin/creatinine ratio(UACR), and reduce pancreatic tissue damage. Metabolomics showed that the Model group produced 109 differential metabolites, of which 74 could be regulated by PTR. 16S rRNA sequencing was done in fecal samples and results showed that PTR could reduce the Firmicutes/Bacteroidetes ratio(F/B) associated with obesity and regulate 3 beneficial bacteria and 1 harmful bacteria. In conclusion, the results showed that PTR could ameliorate T2DM symptoms, metabolic disorder, and gut microbiota imbalance of db/db mice, and it was superior to metformin in some aspects. We suggested for the first time that γ-aminobutyric(GABA) may be involved in the regulation of the microbiota-gut-brain axis(MGB) and thus affect the metabolic disorders associated with T2DM. This study will provide a scientific basis for the development of functional food with PTR.
Keywords: 
Subject: Medicine and Pharmacology  -   Endocrinology and Metabolism

Introduction

Diabetes is a complex metabolic disease associated with disorders of glucose or lipid metabolism, with chronic hyperglycemia as the main feature. The number of adults with diabetes worldwide reached 537 million in 2021, and it is estimated that this number will reach 642 million by 2040, with type 2 diabetes mellitus (T2DM) accounting for the vast majority[1]. Diabetes has become a major challenge for people’s health and brought economic burden worldwide. Currently, dietary intervention, intensive exercise programs, insulin and oral hypoglycemic agents are the most common approaches to combat T2DM, while dietary intervention and lifestyle modification are considered to be more conducive to alleviate T2DM[2]. The theory of “medicine and food homology” was formally proposed in the 1920s and 1930s[3]. Medicinal and food homologous plants are often developed into functional foods due to being rich in bioactive compounds which have health-promoting effects with minimal side effects. As a low glycemic index product, Pueraria thomsonii.
Radix has been used for food and medicine for thousands of years in Asian countries such as China, Thailand, Vietnam, Japan, and others countries. Furthermore, modern studies have emphasized that the bioactive components derived from Pueraria thomsonii Radix can help to improve T2DM in db/db mice, with limited side effects[4,5]. In previous studies, many monomer compounds extracted from Pueraria thomsonii Radix, such as puerarin, daidzein, and genistein, have been shown to regulate animal insulin resistance [6,7], but there is no evidence that the effect of a single monomer component is better than the overall effect of Pueraria thomsonii Radix. Consequently, there is still a limitation to a comprehensive understanding of the integrative mechanism underlying the intake of Pueraria thomsonii Radix to prevent the development of T2DM.
Mounting evidence has demonstrated the combination of metabolism and gut microbiota may serve as a potentially crucial factor in T2DM diagnosis, pharmaceutical discovery, as well as therapeutic response monitoring[8,9]. In the present study, we compared PTR with metformin on meliorating the clinical pathological feature of db/db mice, and then explored the gut microbiota and host metabolomic reactions through UPLC-IM-Q-TOF-MS[10,11] and 16S rRNA gene sequencing to understand the comprehensive antidiabetic mechanism of PTR. We further proposed the potential mechanism that PTR participates in metabolic regulation through the microbiota-gut-brain axis (MGB).

Materials and methods

Preparation of PTR. Pueraria thomsonii Radix, purchased from Jiangzhong Traditional Chinese Medicine Co., Ltd. (Nanchang, Jiangxi) was identified as the dried root of the legume plant Pueraria Thomsonii Benth by Professor Fei Ge and Professor Ronghua Liu of Jiangxi University of Chinese Medicine.PTR was extracted by Jiangxi Xinglin Baima Pharmaceutical Co., Ltd. Raw material of 200 kg was extracted three times continuously with 15 times the amount of water each time, filtered in a basket, combined with the filtrate, and concentrated under reduced pressure. After vacuum drying, 29.4 kg of PTR(yield: 14.70%) was obtained by crushing, which was used for subsequent animal experiments and preparation of test solution. In our previous experiments, UPLC-Q-TOF/MS technique was used to analyze the absorbed components of PTR in rat blood, mainly including 15 polyphenols(Supplementary Table S1)[12], and Elisa was used to analyze the content of total polysaccharides in PTR, these polyphenols and polysaccharides are the main functional components of PTR.
Animals and sample collection. Db/db mice and db/m mice with half female at 7-8 weeks old in an SPF environment were obtained from Nanjing Cavins Biotech, LTD. [Animal Certificate Number: SCXK (Su) 2016-0013] and were raised in the protective environment of the Animal Science and Technology Experimental Center of Jiangxi University of Chinese Medicine (humidity: 50±5%, temperature: 22±2 °C, and light-dark cycle for 12 hours). The treatment of animals during the experiment was in accordance with the eighth edition of the Regulations on the Management of Laboratory Animals of China and was approved by the Ethics Committee of the Laboratory Animal Science and Technology Center of Jiangxi University of Chinese Medicine (Approval No. JZLLSC20210075). After 3 days of adaptive feeding, 11 db /m mice were used as the normal control group (Normal group), 39 db /db mice were modeled for T2DM and randomly divided into 3 groups with 13 mice in each: Model group, PTR group, and Metformin group.
The clinical dosage of Pueraria thomsonii Radix ranges from 10 to 60 g, with 40 g being the most commonly used dosage for human administration. According to the formula of human and mouse body surface area and the yield of extract powder prepared from raw drug, the calculated dose of PTR group was 0.89 g/kg/d.
Each group was respectively given the corresponding drugs by intragastric administration, namely, the mice in Normal group and Model group were intragastrically administered with distilled water, the mice in PTR group were intragastrically administered with PTR 0.89 g/kg/d, and Metformin group was intragastrically administered with the original drug of metformin 0.364 g/kg/d (Sino American Shanghai Squibb Pharmaceutical Co., Ltd.). The mice were given intragastric administration strictly according to their body weight, once a day for 8 consecutive weeks, and the changes of body weight and blood sugar were recorded weekly.
The mice were fasting without water after 18:00 on the night before sampling. The next morning, about 50 uL of blood was collected from the orbital venous plexus of mice and put into common EP tubes for the determination of fasting blood glucose and fasting insulin levels of mice. One hour after the last administration, mice were anesthetized by intraperitoneal injection of 1% Pentobarbital sodium. Whole blood was collected from abdominal aorta into common EP tubes and centrifuged at 8000 r/min for 15 min in a cryogenic centrifuge at 4 ℃. After centrifugation, the supernatant was collected and stored at -80 ℃ for serum metabolomics study. Under aseptic operation, tissues and organs such as pancreas, colon and intestinal contents were carefully collected, and then flash-frozen in liquid nitrogen and stored at -80 ℃ for subsequent experimental analysis. Figure 1 summarizes the schematic diagram of the study.
Analysis of physiological and Biochemical Indexes. The weekly body weight of mice was measured using an electronic analytical balance (Beijing Sartorius Scientific Instruments Co., LTD.). The levels of fasting blood glucose (FBG) were measured weekly using a blood glucose meter and a blood glucose strip (Roche Blood Glucose Health Care). Fasting insulin levels were measured using Mouse Insulin (INS) Elisa kit (Nanjing Jiancheng Bioengineering Institute). Homeostasis model assessment insulin resistance(HOMA-IR) was calculated using the following formula: fasting blood glucose× fasting insulin / 22.5. Urinary glucose(UGLU) was measured using a 7180 automatic biochemical analyzer (Hitachi). Urinary microalbumin was measured using an A15 automatic protein analyzer (Biosystems) and the urinary albumin/creatinine ratio(UACR) was calculated.
Observation of pancreatic hematoxylin-eosin (HE) staining. After tissue sampling, paraffin embedding machine and microtome (Leica) were used to embed and slice the tissue with a thickness of 4 μm. The sections were routinely stained with HE, sealed with neutral gum, and observed by light microscopy and the corresponding images were taken.
Untargeted metabolomics analysis. Thawed the frozen serum samples naturally, taken 50 μL serum, added 200 μL methanol working solution (containing 10.06 µg/mL 2-chloro-L-phenylalanine), mixed them with MTV-100 vortex mixer (Hangzhou Aosheng Group Co., Ltd.) for 30 s, and after standing for 10 min, centrifuged at 12000 r/min for 15 min at 4 °C. The supernatant was taken for untargeted metabolomic analysis in UPLC-IM-Q-TOF-MS (UPLC I-CLASS liquid chromatography system and SYNAPT G2-Si mass spectrometer, Waters, USA). In addition, quality control (QC) samples were prepared by taking an equal amount from each frozen serum sample, vortexed and mixed according to the above method.
Liquid chromatographic separation of the samples was achieved using an ACQUITY UPLC BEN C18 column (2.1ITY UPLC BEm, Waters Corporation). The injection volume was 2 µL, the column temperature was 40 °C, the flow rate was 0.35 mL/min, the eluent A was water containing 0.1% formic acid, and the eluent B was acetonitrile. Solvent gradient setting in positive ion mode: 0-0.2 min, 5% B; 0.2-2 min, 5-20% B; 2-7 min, 20-50% B; 7-17 min, 50-65% B; 17-22 min, 65-80% B; 22-23 min, 80-95% B; 23-26 min, 95-5% B; negative ion mode: 0-0.2 min, 5% B; 0.2-3 min, 5-20% B; 3- 5 min, 20-45% B; 5-7 min, 45-55% B; 7-13 min, 55-65% B; 13-16 min, 65% B; 16-21 min, 65-80% B; 21-23 min, 80-95% B; 23-26 min, 95-5% B.
ESI ion source was used to carry out data acquisition in positive and negative ion mode respectively in HDMSE mode. The mass range (m/z) was 50-1200 Da, ion source temperature was 120°C, curtain rate cone gas flow rate was 50 L/h, desolvation gas flow rate was 800 L/h, desolvation temperature was 400°C, cone voltage was 40 V, collision energy was 20-40 V, scan time and inter scan delay were 0.3 s and 0.015 s, capillary ESI+ was 3 kV, and ESI- was 2.3 kV. In order to ensure the accuracy and reproducibility of the experimental data, the standard product sodium formate was used to establish the mass axis standard curve, and at the same time, leucine enkephalin was used for real-time mass correction, and polyalanine (purchased from Sigma p/n P9003) was used to perform CCS Correction.
In this experiment, in order to ensure the stability of the analysis system as a whole, method validation was performed using quality control (QC) samples. Before collecting each group of samples, QC samples were advanced 5 times to balance the instrument. After stabilization, 1 QC sample was injected for every 5 experimental samples to monitor the operation of the system in real-time. The real-time monitoring system should run to analyze the data in real-time and evaluate the data quality. Finally, the peak intensities of 10 typical mass spectrometry peaks were analyzed, and the value of relative standard deviation (RSD) was used to confirm that the instrument injection was stable and the method was reproducible.
The mass spectrometry raw data obtained by UPLC-IM-QTOF-MS was imported into Progenesis QI V2.0 software (QI, Waters, USA) for normalization processing, and then the data matrix was exported and saved. After grouping the data, exported the data to EZinfo3.0 (Waters, USA) software to establish labels (VIP>1 and p<0.05) to screen differential metabolites. The SIMCA14.1 (Swedish Umetrics Company) software was used to perform principal component analysis (PCA) and orthogonal partial least squares discriminant analysis (OPLS-DA) on the data obtained in the previous step. According to the data of mass-to-charge ratio, retention time, CCS value, and secondary mass spectrometry fragment of the differential metabolites, they were matched and identified in the HMDB database, and then were analyzed by Pathway Analysis in the Metabo Analyst 5.0 (https://www.metaboanalyst.ca) online database. The identified differential metabolites were enriched for metabolic pathways, and the differential metabolites were finally determined to obtain metabolic pathways and conduct pathway analysis.
Gut microbiota analysis. Gut microbiota analysis was performed by Shanghai Magi Biomedical Technology Co., LTD. (Shanghai, China). DNA was extracted from differential intestinal content samples using E.Z.N.A.® Soil DNA Kit (Omega Biotek, Norcross, GA, U.S.). The quality of DNA extraction was examined using 1% agarose gel electrophoresis, and the concentration and purity of DNA were determined using a NanoDrop2000 ultramicro spectrophotometer (Thermo Scientific, Wilmington, USA). The hypervariable region V3-V4 of the bacterial 16S rRNA gene was amplified with primer pairs 338F (5’-ACTCCTACGGGAGGCAGCAG-3’) and 806R (5’-GGACTACHVGGGTWTCTAAT-3’) by an ABI GeneAmp® 9700 PCR thermocycler (ABI, CA, USA). There were three replicates for each sample. The PCR products of each sample were mixed in triplicate and then the PCR products were recovered on 2% agarose gel. AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA) was used to purify the recovered product and detect it by 2% agarose gel electrophoresis. Quantus™ Fluorometer (Promega, USA) was used to detect and quantify the recovered products. The NEXTflexTM Rapid DNA-Seq Kit (Bioo Scientific, USA) was used for library construction. Sequencing was performed using Illumina’s Miseq PE300 platform. The fastp software (https://github.com/OpenGene/fastp, version 0.20.0) was used to quality control of original sequencing sequence. FLASH software (http://www.cbcb.umd.edu/software/flash, version 1.2.7) was used to pick up the fight. UPARSE software (http://drive5.com/uparse/,Version 7.1) was used to cluster OTU sequences according to 97% similarity. RDP classifier (http://rdp.cme.msu.edu/, version 2.2) was used to annotate each sequence for species classification. The comparison threshold was set to 70% for Silva 16S rRNA database (V138).
Statistical analyses. These values were defined as mean ± standard deviation (SD). The difference between the two groups of metabolites was calculated using student’s t-test by GraphPad prism 9 (GraphPad Software Inc., La Jolla, CA, USA). The difference between the two groups of gut microbiota was calculated using the Wilcox rank sum test. Compared with Normal group or Model group, the significance of the difference was described as * P < 0.05, ** P < 0.01, and *** P < 0.001.

Results

PTR improves major physiological and biochemical Indexes in db/db mice. Obesity provokes the impairment of insulin sensitivity, ultimately developing into diabetes. To study the modulatory effects of PTR on diabetes, the db/db mice were applied in this study. Due to deficient leptin receptor gene, db/db mice can spontaneously develop the obesity and chronic hyperglycemia. As shown in Figure 3a, the body weight of Model group was significantly higher than that of the group taken normal diet, indicating db/db mice have developed obesity. Although the body weight of mice in PTR group was higher than that of Normal group, while PTR group was lower than (Figure 2a). Compared with Normal group, Model group significantly increased the level of FBG and HOMA-IR of the mice, indicating that db/db mice developed T2DM in the present study. However, FBG and HOMA-IR were significantly decreased in mice after administration of PTR and metformin (Figure 2b, c). Compared with Normal group, UGLU and UACR of the mice in Model group were significantly increased, and both UGLU and UACR were significantly decreased after administration of PTR and metformin (Figure 2d, e). In summary, the role of PTR in improving key physiological and biochemical parameters in db/db mice was comparable to that of metformin.
PTR alleviates pancreatic tissue damage in db/db mice. The observation of mouse pancreas HE staining (Figure 3a–d) showed that compared with Normal group, the islets in Model group were atrophied, the diameter was reduced, the number of islet cells was reduced, and the boundary with pancreatic acinar cells was unclear. Moreover, a large number of pancreatic acinar cells were necrotic and degenerative. After administration of PTR and metformin, a clear boundary of islet cells and increased numbers from the treatment groups were observed in the present study. It suggested that the pancreatic tissue damage of db/db mice can be potentially alleviated by the administration of PTR and metformin.
The modulatory effects of PTR on serum metabolites in db/db mice. A significant separation was observed among Normal group, Model group, PTR, or Metformin groups in the PCA and OPLS-DA score graphs. The R2 and Q2 values in OPLS-DA were smaller than the original values, indicating that the established discriminant model was reliable. The above results indicated that PTR and metformin could significantly alter the composition of serum metabolites in db/db mice (Figure 4). VIP > 1 and p < 0.05 have been set up to determine the potential biomarkers. In this study, 13 kinds of substances were mainly assessed, including lipids, lipid-like molecules, bile acids, alcohols and their derivatives, glycerol phospholipids, nucleosides, nucleotides and analogues, steroids and steroid derivatives, organic oxygen compounds, and organic heterocyclic compounds. There was a total of 109 differential metabolites among the above categories of compounds identified in this study (Supplementary Table S2). Heat map analysis of these relative metabolites were shown in Figure 5a. The results of analysis of differential metabolites with callback after administration of db/db mice showed that there were 74 callback differential metabolites after administration of PTR, and their abundance values were shown in Supplementary Table S3. In addition, there were 49 callback differential metabolites after administration of metformin, and their abundance values were shown in Supplementary Table S4. To determine the difference between PTR group and Metformin group, the concentration table of differential metabolites introduced into PTR and metformin was used to analyze the enrichment path of their differential metabolites in MetaboAnalyst Path analysis and KEGG (https://www.metaboanalyst.ca/). Finally, 40 metabolic pathways (Supplementary Table S5) were particularly enriched by PTR, and 26 metabolic pathways are enriched by metformin (Supplementary Table S6). There are 24 main metabolic pathways overlapped between PTR and metformin, as shown in the Bubble Diagram (Figure 5b,c).
PTR improved the gut microbiota in db/db mice. A total of 2,866,188 colon content sequences were detected by 16S rRNA gene sequencing in this study. The dilution curve was close to the saturation plateau and the sequencing was reasonable (Figure 6a). Based on the 16S rRNA gene sequencing assessment, Chao index reflects the richness of gut microbiota, Shannon index reflects community diversity, and Coverage index refers to the coverage of the sample library. The higher the value, the higher the probability of detecting sequences in the sample. The indexes of each group of gut microbiota were analyzed as shown in Figure 6c–e. Compared with Normal group, there was no significant difference in the Chao index and Shannon index of Model group, indicating that there was no significant change in the richness and diversity of gut microbiota in db/db mice. Compared with Model group, the Chao index and Shannon index of PTR group were increased but no significant difference being detected. Furthermore, a significant decrease of Chao index and a moderate decrease of Shannon index were identified in Metformin group. Compared with Metformin group, a significant increase of Chao index, Shannon index and Coverage index were identified in PTR group, suggesting that intake of PTR contributes to restoring the gut microbiota structure.
In the Beta diversity analysis, the Bray Curtis distance matrix PCoA analysis based on the OTU level showed a clear separation between Model group and Normal group. Furthermore, we observed that all groups tend to be a distinct cluster, suggesting that each group in our study had a distinct gut microbiota composition at the OTU levels (Figure 6b). Based on OTU sequence analysis at 97% similarity level, 778 OTU sequences were found in Normal group, 872 in Model group, 861 in PTR group, and 850 in Metformin group (Figure 7a). At the phylum level, the dominant bacteria in the gut microbiota of the four groups of samples were Bacteroidota, Firmicutes, Proteobacteria, Desulfobacterota, Actinobacteriota, Verrucomicrobiota, and Campilobacterota (Figure 7b). According to the histogram of the gut microbiota at the phylum level, the relative abundance of Firmicutes increased and that of Bacteroidetes decreased in Model group, leading to the increase of the Firmicutes/Bacteroidetesratio(F/B) without a significant difference. Compared with Model group, no significant difference of F/B was found in PTR group. However, the F/B of Metformin group was significantly increased (P≤0.01) compared to Model group. Likely, the F/B of Metformin group was significantly increased (P≤0.01) compared with PTR group.
As shown in Figure 7c, the dominant bacteria in the gut microbiota of the four groups of samples at the genus level were norank_f__Muribaculaceae, Escherichia-Shigella, Lactobacillus, unclassified_f__Lachnospiraceae, Akkermansia, Alistipes. The Wilcox rank sum test was performed at the genus level to obtain the corresponding different bacterium. The different bacterium between Normal group and Model group included norank_f__Muribaculaceae, Escherichia-Shigella,Alistipes, Dubosiella, Lactobacillus, Ileibacterium, unclassified_f__Lachnospiraceae, Bifidobacterium, Alloprevotella, norank_f__Lachnospiraceae, Muribaculum,
Odoribacter, Rikenellaceae_RC9_gut_group, Lachnoclostridium, norank_f__norank_o__Clostridia_UCG-014, etc. (Figure 7d). The differential bacterium between PTR and Model groups included Parabacteroides, Akkermansia, norank_f__Desulfovibrionaceae, norank_f__Ruminococcaceae, Bifidobacterium,
unclassified_f__Oscillospiraceae, Anaerotruncus, Oscillibacter, Colidextribacter,
norank_f__norank_o__Clostridia_UCG-014, Muribaculum, Intestinimonas, Marvinbryantia,
Rikenellagca-900066575 etc. (Figure 7e). The differential bacterium between Meltormin group and Model group included norank_f__Muribaculaceae, Lactobacillus, Escherichia-Shigella, Akkermansia, Lachnospiraceae_NK4A136_group, Desulfovibrio, Bacteroides, norank_f__Lachnospiraceae,Prevotellaceae_NK3B31_group, Enterorhabdus, norank_f__Desulfovibrionaceae, norank_f__Ruminococcaceae, Parabacteroides,unclassified_o__Bacteroidales, Bifidobacterium etc. (Figure 7f). The differential bacterium between Metformin group and PTR group included norank_f__Muribaculaceae, Lactobacillus, Escherichia-Shigella, Akkermansia, Bacteroides, Lachnospiraceae_NK4A136_group, Desulfovibrio, norank_f__Eubacterium_coprostanoligenes_group, Parabacteroides, unclassified_o__Bacteroidales, Blautia, Streptococcus, Alloprevotella, norank_f__Oscillospiraceae, Proteus, muribaculum, etc. (Figure 7g).
Correlation analysis of gut microbiota and differential metabolites. Spearman correlation analysis was used to calculate the correlation between gut microbiota and differential metabolites in serum collected from db/db mice after treated with PTR (Figure 8). Among them, Muribaculum was positively correlated with 3 - O-Sulfogalactosylceramide ( d18: 1 / 18 : 1 ( 9Z ) ) and Taurocholic acid, while negatively correlated with 5 - Acetylamino-6 - formylamino-3 - methyluracil and PE ( 20 : 3 ( 8Z, 11Z, 14Z ) / P-18 : 0 ) ( | r | > 0.5, p < 0.05, 0.01, 0.001 ). Additionally, Desulfovibrio was positively correlated with 3beta-Hydroxypregn-5-en-20-one sulfate, 11b,21-Dihydroxy-3,20-oxo-5b-pregnan-18-al,13-L-Hydroperoxylinoleic acid, Leukotriene D4, Deoxyguanosine, Oxidized glutathione, Melibiitol, Cortisol, Inosine, PC (22:4(7Z,10Z,13Z,16Z)/P-16:0) (|r|>0.5,p<0.05,0.01,0.001). Lastly, Blautia was positively correlated with 5a-Pregnane-3, 20-dione, Riboflavin reduced, 2-Methoxy-estradiol-17b 3-glucuronide, Leukotriene D4, SM(d18:1/26:1(17Z)), Orotidylic acid, PC (22:4(7Z,10Z,13Z,16Z)/P-16:0) (|r|>0.5,p<0.05,0.01,0.001). The above results indicated that the changes of serum metabolites were closely linked to the changes of gut microbiota after PTR supplementation.

Discussion

Since ancient times, Pueraria thomsonii Radix is known as “Asian Ginseng”, and is recognised as a medicinal and food homologous plant. In this study, the therapeutic effects on db/db mice were firstly demonstrated, and then the gut microbiota and host metabolomic responses were explored to understand the integrative mechanism underlying PTR ameliorating T2DM and the interplay between the host and microbial.
Metformin is often used in T2DM patients who are unsatisfied with dietary control alone, especially those who are obese and hyperinsulinemic, not only for its hypoglycaemic effect but also for its effect on weight loss and hyperinsulinemia. Metformin is also commonly used in the development of new anti-diabetic drugs to evaluate the efficacy of new drugs against diabetes. In experimental diabetes models, metformin was effective in downregulating blood glucose levels and improving insulin resistance in diabetic mice and rats. In this study, we demonstrated for the first
time that PTR has therapeutic effects in the treatment of T2DM, improving body weight, reducing FBG, UGLU, HOMA-IR and UACR in mice, and reducing pancreatic tissue damage. We found that PTR was comparable to metformin in down-regulating blood glucose levels and improving insulin resistance. This further suggests that the anti-diabetic potential of PTR is in urgent need of further clinical studies, or the use of dietary therapy to control diabetes.
Regulating metabolism is an effective measure for the treatment of diabetes. In this study, PTR reversed 74 of the 109 differential metabolites in db/db mice, mainly including purines, pyrimidines, glycerophospholipids, and steroids. Purine and pyrimidine metabolism are central to the adenosine system and play a key role in the regulation of glucose homeostasis as well as the pathophysiology of diabetes[13,14].In this study, PTR affected the purine metabolic pathway by regulating the levels of Inosine, Inosinic acid and Deoxyguanosine in db/db mice, and the pyrimidine metabolic pathway by regulating the concentrations of Orotidylic acid, dUMP and 5-Thymidylic acid, thereby improving T2DM. Triglycerides are a clinically important feature of adult-onset T2DM[15].In this study, PTR could modulate PE (20:3(8Z,11Z,14Z)/P-18:0), PC (22:4(7Z,10Z,13Z,16Z)/P-16:0) and LysoPA (18:1(9Z)/0:0) concentrations in db/db mice. Two biomarkers, cortisol and cortisone, which play an important role in steroid hormone regulation, affect the energy metabolism associated with T2DM[16].In this experiment, PTR effectively regulated 5a-Pregnane-3,20-dione, Corticosterone, 11b,21-Dihydroxy-3, 20-oxo-5b-pregnan-18-al, Dehydroepiandrosterone sulfate Cortisol, Estrone, and 2-Methoxy-estradiol-17b 3-glucuronide concentrations. In conclusion, PTR was effective in alleviating T2DM by regulating metabolism.
There is growing evidence that dysbiosis of the intestinal flora is a typical feature of the pathogenesis of diabetes. Clinical studies have found that Bifidobacterium, Bacillus spp., Clostridium pretense and Ackermannia spp. are negatively associated with T2DM, and that Clostridium tumefaciens, Clostridium perfringens, and Braunschweiger spp. are positively associated with T2DM. In animal models of experimental diabetes, the species and number of beneficial bacteria were significantly lower, while pathogenic bacteria were significantly higher[17].In this study,we found that PTR significantly upregulated the abundance of Muribaculum, norank_f__Muribaculaceae, and Parabacteroides in db/db mice while significantly reducing the abundance of the harmful bacterium Desulfovibrio. A study has demonstrated that Plukenetia volubilis L. leaf significantly increased Muribaculum abundance in db/db mice, a beneficial bacterium that can improve inflammation, dyslipidemia, and glucose intolerance[18,19].Norank_f_Muribaculaceae was a beneficial dominant bacterium in the mouse intestine and was negatively correlated with the glucose levels of postprandial blood[20]. Parabacteroides is a beneficial bacterium that significantly improves insulin resistance and has been shown in vitro to activate the intestinal gluconeogenic pathway through the conversion of succinate, thereby exerting a hypoglycaemic effect[21]. Desulfovibrio is a harmful bacterium whose increased abundance can cause inflammation and insulin resistance[22]. Moreover, it is well known that an increase of the F/B has been considered as a marker of obesogenic flora[23,24].In this study, we also found an increased F/B in db/db mice. In a word, this study provides a basis for further investigation of the anti-obesity role of key PTR-dependent bacteria.
In addition, Spearman’s correlation analysis further revealed a positive correlation between Desulfovibrio and Cortisol in db/db mice. Cortisol is an important product of the steroid hormone biosynthetic pathway and is considered an active indicator of the Hypothalamic-pituitary-adrenal(HPA)axis. A recent study has shown that cortisol is associated with dyslipidemia and impaired glucose metabolism[25]. Desulfovibrio has been reported to be negatively correlated with γ-aminobutyric(GABA)[26], which can inhibit glucagon secretion and enhance insulin secretion[27]. Notably, hypothalamic Cortisol concentrations are influenced by GABA, and reduced GABA levels lead to over-activation of the HPA axis, resulting in increased Cortisol levels[28,29]. Recent research hotspots have reported a stronger bidirectional communication between the gut and brain through the nervous, endocrine, and immune systems, i.e., MGB[30,31]. In this experiment, PTR reduced the abundance of Desulfovibrio and the levels of Cortisol metabolism. Therefore, we speculate that PTR administration resulted in a decrease in Desulfovibrio abundance, leading to an increase in GABA levels in the gut, which in turn inhibited HPA axis activation via the enteric nervous system (ENS) via the MGB pathway, contributing to a decrease in Cortisol metabolism levels (Figure 9). This is the focus of our subsequent studies, and perhaps related studies will be reported in the near future.

Supplementary Materials

The following supporting information can be downloaded at the website of this paper posted on Preprints.org. Supplementary Table S1:LC-MS data of the absorbed components of PTR in rat blood; Supplementary Table S2: Differential metabolites; Supplementary Table S3: Differential metabolites for PTR callbacks; Supplementary Table S4: Differential metabolites for Metformin callbacks; Supplementary Table S5: Analysis of PTR metabolic pathway enrichment; Supplementary Table S6: Analysis of Metformin metabolic pathway enrichment.

Author Contributions

Jiarong Li and Yuhui Liu have contributed equally to this work. Conceptualization, Ronghua Liu, Hua Zhang and Jiarong Li; Methodology, Yuhui Liu, Hui Ouyang and Yong Sun; Software, Lifang Wang and Jiaxing Huang; Data curation, Jiarong Li, Junchang Chen and Lifang Wang; Formal Analysis, Jiarong Li; Investigation, Mingyao Li; Resources, Yuhui Liu; Writing—original draft preparation, Jiarong Li, Yuhui Liu and Lifang Wang; Writing—review and editing, Youbao Zhong and Hua Zhang; Supervision, Xiangyan Tan; Visualization, Yue Hu. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (2017YFC1702902).

Data Availability Statement

Data is contained within the article or supplementary material.

Conflicts of Interest

There are no conflicts of interest to declare.

Abbreviations

PTR Pueraria thomsonii Radix water extract
T2DM Type 2 diabetes mellitus
FBG Fasting blood glucose
UGLU Urinary glucose
HOMA-IR Homeostasis model assessment insulin resistance
UACR Urinary albumin/creatinine ratio
CCS Collision cross section
OPLS-DA Orthogonal partial least squares discriminate analysis
GABA γ-aminobutyric acid
HPA Hypothalamic-pituitary-adrenal
MGB The microbiota-gut-brain axis
F/B The Firmicutes/Bacteroidetes ratio
HE Hematoxylin-eosin
QC Quality control
RSD Relative standard deviation

References

  1. Pan, L.; Li, Z.; Wang, Y.; Zhang, B.; Liu, G.; Liu, J. Network pharmacology and metabolomics study on the intervention of traditional Chinese medicine Huanglian Decoction in rats with type 2 diabetes mellitus. J. Ethnopharmacol. 2020, 258, 112842. [Google Scholar] [CrossRef] [PubMed]
  2. Yang, Z.; Wang, D.; Li, Y.; Zhou, X.; Liu, T.; Shi, C.; Li, R.; Zhang, Y.; Zhang, J.; Yan, J.; et al. Untargeted metabolomics analysis of the anti-diabetic effect of Red ginseng extract in Type 2 diabetes Mellitus rats based on UHPLC-MS/MS. Biomed. Pharmacother. 2021, 146, 112495. [Google Scholar] [CrossRef]
  3. Ji, M.-Y.; Bo, A.; Yang, M.; Xu, J.-F.; Jiang, L.-L.; Zhou, B.-C.; Li, M.-H. The Pharmacological Effects and Health Benefits of Platycodon grandiflorus—A Medicine Food Homology Species. Foods 2020, 9, 142. [Google Scholar] [CrossRef]
  4. Z. Zhang, T. N. Lam and Z. Zuo, Radix Puerariae: an overview of its chemistry, pharmacology, pharmacokinetics, and clinical use, The Journal of Clinical Pharmacology, 2013, 53, 787–811.
  5. Lee, J.-S. Effects of soy protein and genistein on blood glucose, antioxidant enzyme activities, and lipid profile in streptozotocin-induced diabetic rats. Life Sci. 2006, 79, 1578–1584. [Google Scholar] [CrossRef] [PubMed]
  6. Duru, K.C.; Mukhlynina, E.A.; Moroz, G.A.; Gette, I.F.; Danilova, I.G.; Kovaleva, E.G. Anti-diabetic effect of isoflavone rich kudzu root extract in experimentally induced diabetic rats. J. Funct. Foods 2020, 68, 103922. [Google Scholar] [CrossRef]
  7. Carlson, S.; University of Alabama at Birmingham; Prasain, J. K.; Peng, N.; Dai, Y.; Wyss, J.M. Acute and Chronic Kudzu Improves Plasma Glucose Tolerance in Non-Diabetic CD-1 Mice. J. Endocrinol. Diabetes Mellit. 2014, 2, 70–77. [Google Scholar] [CrossRef]
  8. Luo, D.; Dong, X.; Huang, J.; Huang, C.; Fang, G.; Huang, Y. Pueraria lobata root polysaccharide alleviates glucose and lipid metabolic dysfunction in diabetic db/db mice. Pharm. Biol. 2021, 59, 380–388. [Google Scholar] [CrossRef] [PubMed]
  9. Song, X.; Dong, H.; Zang, Z.; Wu, W.; Zhu, W.; Zhang, H.; Guan, Y. Kudzu Resistant Starch: An Effective Regulator of Type 2 Diabetes Mellitus. Oxidative Med. Cell. Longev. 2021, 2021, 1–15. [Google Scholar] [CrossRef]
  10. Hernández-Mesa, M.; D’atri, V.; Barknowitz, G.; Fanuel, M.; Pezzatti, J.; Dreolin, N.; Ropartz, D.; Monteau, F.; Vigneau, E.; Rudaz, S.; et al. Interlaboratory and Interplatform Study of Steroids Collision Cross Section by Traveling Wave Ion Mobility Spectrometry. Anal. Chem. 2020, 92, 5013–5022. [Google Scholar] [CrossRef]
  11. Paglia, G.; Williams, J.P.; Menikarachchi, L.; Thompson, J.W.; Tyldesley-Worster, R.; Halldórsson, S.; Rolfsson, O.; Moseley, A.; Grant, D.; Langridge, J.; et al. Ion Mobility Derived Collision Cross Sections to Support Metabolomics Applications. Anal. Chem. 2014, 86, 3985–3993. [Google Scholar] [CrossRef] [PubMed]
  12. Yu, H.-H.; Meng, X.-W.; Li, J.-R.; Ouyang, H.; Yin, Q.-Q.; Lin, L.-H.; Zhu, W.-F.; Liu, R.-H. [Comparative study of rat serum pharmacochemistry between Puerariae Lobatae Radix and Puerariae Thomsonii Radix based on UPLC-Q-TOF-MS]. . 2022, 47, 528–536. [Google Scholar] [CrossRef]
  13. Mousavi, S.M.; Jalilpiran, Y.; Karimi, E.; Aune, D.; Larijani, B.; Mozaffarian, D.; Willett, W.C.; Esmaillzadeh, A. Dietary Intake of Linoleic Acid, Its Concentrations, and the Risk of Type 2 Diabetes: A Systematic Review and Dose-Response Meta-analysis of Prospective Cohort Studies. Diabetes Care 2021, 44, 2173–2181. [Google Scholar] [CrossRef] [PubMed]
  14. Yuan, E.-F.; Yang, Y.; Cheng, L.; Deng, X.; Chen, S.-M.; Zhou, X.; Liu, S.-M. Hyperglycemia affects global 5-methylcytosine and 5-hydroxymethylcytosine in blood genomic DNA through upregulation of SIRT6 and TETs. Clin. Epigenetics 2019, 11, 63. [Google Scholar] [CrossRef]
  15. Shi, Z.-Q.; Wang, L.-Y.; Zheng, J.-Y.; Xin, G.-Z.; Chen, L. Lipidomics characterization of the mechanism of Cynomorium songaricum polysaccharide on treating type 2 diabetes. J. Chromatogr. B 2021, 1176, 122737. [Google Scholar] [CrossRef]
  16. Ji, Y.; Yao, Y.; Duan, Y.; Zhao, H.; Hong, Y.; Cai, Z.; Sun, H. Association between urinary organophosphate flame retardant diesters and steroid hormones: A metabolomic study on type 2 diabetes mellitus cases and controls. Sci. Total. Environ. 2020, 756, 143836. [Google Scholar] [CrossRef] [PubMed]
  17. Tao, S.; Zheng, W.; Liu, Y.; Li, L.; Li, L.; Ren, Q.; Shi, M.; Liu, J.; Jiang, J.; Ma, H.; et al. Analysis of serum metabolomics among biopsy-proven diabetic nephropathy, type 2 diabetes mellitus and healthy controls. RSC Adv. 2019, 9, 18713–18719. [Google Scholar] [CrossRef]
  18. Yuan, Y.; Zhou, J.; Zheng, Y.; Xu, Z.; Li, Y.; Zhou, S.; Zhang, C. Beneficial effects of polysaccharide-rich extracts from Apocynum venetum leaves on hypoglycemic and gut microbiota in type 2 diabetic mice. Biomed. Pharmacother. 2020, 127, 110182. [Google Scholar] [CrossRef]
  19. J. Lin, J. Wen, N. Xiao, Y.-T. Cai, J. Xiao, W. Dai, J.-P. Chen, K.-W. Zeng, F. Liu, B. Du and P. Li, Anti-diabetic and gut microbiota modulation effects of sacha inchi (Plukenetia volubilis L.) leaf extract in streptozotocin-induced type 1 diabetic mice, Journal of the Science of Food and Agriculture, 2022, 102, 4304–4312.
  20. S. Li, M. S. Li, M. Wang, C. Li, Q. Meng, Y. Meng, J. Ying, S. Bai, Q. Shen and Y. Xue, Beneficial Effects of Partly Milled Highland Barley on the Prevention of High-Fat Diet-Induced Glycometabolic Disorder and the Modulation of Gut Microbiota in Mice, Nutrients, 2022, 14.
  21. Wang, K.; Liao, M.; Zhou, N.; Bao, L.; Ma, K.; Zheng, Z.; Wang, Y.; Liu, C.; Wang, W.; Wang, J.; et al. Parabacteroides distasonis Alleviates Obesity and Metabolic Dysfunctions via Production of Succinate and Secondary Bile Acids. Cell Rep. 2019, 26, 222–235e5. [Google Scholar] [CrossRef]
  22. Zhang, L.-Y.; Zhou, T.; Zhang, Y.-M.; Xu, X.-M.; Li, Y.-Y.; Wei, K.-X.; He, J.-P.; Ding, N.; Zhou, H.; Liu, Y.-Q. Guiqi Baizhu Decoction Alleviates Radiation Inflammation in Rats by Modulating the Composition of the Gut Microbiota. Evidence-Based Complement. Altern. Med. 2020, 2020, 1–13. [Google Scholar] [CrossRef] [PubMed]
  23. Ussar, S.; Fujisaka, S.; Kahn, C.R. Interactions between host genetics and gut microbiome in diabetes and metabolic syndrome. Mol. Metab. 2016, 5, 795–803. [Google Scholar] [CrossRef] [PubMed]
  24. F. Alberto Arguelles-Tello, A. Kammar-Garcia, C. Adolfo Trejo-Jasso, J. Carlos Huerta-Cruz, L. Marcela Barranco-Garduno, H. Isaac Rocha-Gonzalez and J. Gerardo Reyes-Garcia, Metformin improves the weight reduction effect of mazindol in obese Mexican, International Journal of Clinical Pharmacology and Therapeutics, 2022, 60, 336–345.
  25. Jin, Y.; Wei, D.; Liu, P.; Chen, F.; Li, R.; Zhang, J.; Zhang, R.; Liu, Z.; Huo, W.; Li, L.; et al. Serum Cortisol, 25 (OH)D, and Cardiovascular Risk Factors in Patients with Type 2 Diabetes Mellitus. Int. J. Endocrinol. 2022, 2022, 1–9. [Google Scholar] [CrossRef] [PubMed]
  26. Tian, X.-Y.; Xing, J.-W.; Zheng, Q.-Q.; Gao, P.-F. 919 Syrup Alleviates Postpartum Depression by Modulating the Structure and Metabolism of Gut Microbes and Affecting the Function of the Hippocampal GABA/Glutamate System. Front. Cell. Infect. Microbiol. 2021, 11. [Google Scholar] [CrossRef] [PubMed]
  27. Wang, Q.; Prud'Homme, G.; Wan, Y. GABAergic system in the endocrine pancreas: a new target for diabetes treatment. Diabetes, Metab. Syndr. Obesity: Targets Ther. 2015, 8, 79–87. [Google Scholar] [CrossRef] [PubMed]
  28. Colmers, P.L.; Bains, J.S. Balancing tonic and phasic inhibition in hypothalamic corticotropin-releasing hormone neurons. J. Physiol. 2018, 596, 1919–1929. [Google Scholar] [CrossRef] [PubMed]
  29. Gao, Y.; Zhou, J.-J.; Zhu, Y.; Kosten, T.; Li, D.-P. Chronic Unpredictable Mild Stress Induces Loss of GABA Inhibition in Corticotrophin-Releasing Hormone-Expressing Neurons through NKCC1 Upregulation. Neuroendocrinology 2016, 104, 194–208. [Google Scholar] [CrossRef] [PubMed]
  30. Luan, H.; Wang, X.; Cai, Z. Mass spectrometry-based metabolomics: Targeting the crosstalk between gut microbiota and brain in neurodegenerative disorders. Mass Spectrom. Rev. 2017, 38, 22–33. [Google Scholar] [CrossRef]
  31. Rosario, D.; Boren, J.; Uhlen, M.; Proctor, G.; Aarsland, D.; Mardinoglu, A.; Shoaie, S. Systems Biology Approaches to Understand the Host–Microbiome Interactions in Neurodegenerative Diseases. Front. Neurosci. 2020, 14, 716. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the study.
Figure 1. Schematic diagram of the study.
Preprints 69444 g001
Figure 2. Effect of Pueraria thomsonii Radix water extract(PTR) on type 2 diabetes mellitus (T2DM) symptoms in db/db mice: (a)Weight; (b)Fasting blood glucose (FBG); (c)Homeostasis model assessment insulin resistance(HOMA-IR); (d)Urinary glucose (UGLU); (e)Urinary albumin/creatinine ratio(UACR). Compared with Normal group or Model group, the significance of the difference was described as *P<0.05, **P<0.01, ***P<0.001.
Figure 2. Effect of Pueraria thomsonii Radix water extract(PTR) on type 2 diabetes mellitus (T2DM) symptoms in db/db mice: (a)Weight; (b)Fasting blood glucose (FBG); (c)Homeostasis model assessment insulin resistance(HOMA-IR); (d)Urinary glucose (UGLU); (e)Urinary albumin/creatinine ratio(UACR). Compared with Normal group or Model group, the significance of the difference was described as *P<0.05, **P<0.01, ***P<0.001.
Preprints 69444 g002
Figure 3. Effect of PTR on histopathology of pancreas in db/db mice:(a) Normal group; (b) Model group; (c) PTR group; (d) Metformin group.Pancreas tissues are stained with H&E (100×).
Figure 3. Effect of PTR on histopathology of pancreas in db/db mice:(a) Normal group; (b) Model group; (c) PTR group; (d) Metformin group.Pancreas tissues are stained with H&E (100×).
Preprints 69444 g003
Figure 4. Untargeted metabolomics analysis of fecal samples by using SYNAPT G2-Si, including PCA, OPLS-DA and Permutation test in the positive(POS) and in negative(NEG).
Figure 4. Untargeted metabolomics analysis of fecal samples by using SYNAPT G2-Si, including PCA, OPLS-DA and Permutation test in the positive(POS) and in negative(NEG).
Preprints 69444 g004
Figure 5. Effect of PTR on metabolomics of serum in db/db mice:(a) Heat map analysis of potential biomarker relative content;(b) Summary plot for pathway analysis of PTR (1–15,18-19 metabolic pathway satisfies the conditions of p < 0.05 or impact value > 0.1);(c) Summary plot for pathway analysis of Metformin (1-18 metabolic pathway satisfies the conditions of p < 0.05 or impact value > 0.1).
Figure 5. Effect of PTR on metabolomics of serum in db/db mice:(a) Heat map analysis of potential biomarker relative content;(b) Summary plot for pathway analysis of PTR (1–15,18-19 metabolic pathway satisfies the conditions of p < 0.05 or impact value > 0.1);(c) Summary plot for pathway analysis of Metformin (1-18 metabolic pathway satisfies the conditions of p < 0.05 or impact value > 0.1).
Preprints 69444 g005
Figure 6. Effects of PTR on the overall structure of gut microbiota in db/db mice:(a) Rarefaction curves among different groups; (b)Sobs index; (c)Chao index; (d)Shannon index; (e)Unweighted unifrac-based principal coordinates analysis.
Figure 6. Effects of PTR on the overall structure of gut microbiota in db/db mice:(a) Rarefaction curves among different groups; (b)Sobs index; (c)Chao index; (d)Shannon index; (e)Unweighted unifrac-based principal coordinates analysis.
Preprints 69444 g006
Figure 7. Effects of PTR on the overall composition of gut microbiota in db/db mice:(a) Venn diagram analysis at OUT level; (b)The relative abundance of gut microflora at phyla level; (c) The relative abundance of gut microflora at genus level. (d) Wilcoxon rank-sum test bar plot on Genus level (P<0.05), the difference between Normal and Model group; (e)the difference between Model and PTR group; (f) the difference between Model and Metformin group; (g) the difference between PTR and Metformin group.Compared with Normal group or Model group, the significance of the difference was described as *P<0.05, **P<0.01, ***P<0.001.
Figure 7. Effects of PTR on the overall composition of gut microbiota in db/db mice:(a) Venn diagram analysis at OUT level; (b)The relative abundance of gut microflora at phyla level; (c) The relative abundance of gut microflora at genus level. (d) Wilcoxon rank-sum test bar plot on Genus level (P<0.05), the difference between Normal and Model group; (e)the difference between Model and PTR group; (f) the difference between Model and Metformin group; (g) the difference between PTR and Metformin group.Compared with Normal group or Model group, the significance of the difference was described as *P<0.05, **P<0.01, ***P<0.001.
Preprints 69444 g007
Figure 8. Spearman analysis between bacterial genera and microbial metabolites. The R values are represented by gradient colors, where pink and blue cells indicate positive and negative correlations, respectively.* p < 0.05, ** p < 0.01 and *** p < 0.001.
Figure 8. Spearman analysis between bacterial genera and microbial metabolites. The R values are represented by gradient colors, where pink and blue cells indicate positive and negative correlations, respectively.* p < 0.05, ** p < 0.01 and *** p < 0.001.
Preprints 69444 g008
Figure 9. The schematic of mechanism of action of PTR.
Figure 9. The schematic of mechanism of action of PTR.
Preprints 69444 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated