Xuehua Xie1,2†
Jiabao Liao3,4†Yuanliang Ai5Jinmei Gao6Jie Zhao2Fei Qu3Chao Xu2
Zhaiyi Zhang7
Weibo Wen2*
Huantian Cui8*
Hongwu Wang9*- 1First College of Clinical Medicine, Nanjing University of Traditional Chinese Medicine, Jiangsu, China
- 2Department of Endocrinology, Yunnan Provincial Hospital of Chinese Medicine, Yunnan, China
- 3Department of Emergency, Jiaxing Hospital of Traditional Chinese Medicine, Zhejiang, China
- 4Jiaxing Key Laboratory of Diabetic Angiopathy Research, Jiaxing Hospital of Traditional Chinese Medicine, Zhejiang, China
- 5Department of Orthopedics, Kunming Municipal Hospital of Traditional Chinese Medicine, Yunnan, China
- 6Department of Rehabilitation, Fujian People’s Hospital of Traditional Chinese Medicine, Fujian, China
- 7College of Integrated Chinese and Western Medicine, Tianjin University of Traditional Chinese Medicine, Tianjin, China
- 8Shandong Provincial Key Laboratory of Animal Cell and Developmental Biology, School of Life Sciences, Shandong University, Shandong, China
- 9College of Traditional Chinese Medicine, Tianjin University of Traditional Chinese Medicine, Tianjin, China
Pi-Dan-Jian-Qing decoction (PDJQ) can been used in the treatment of type 2 diabetes mellitus (T2DM) in clinic. However, the protective mechanisms of PDJQ on T2DM remain unknown. Recent studies have shown that the changes in gut microbiota could affect the host metabolism and contribute to progression of T2DM. In this study, we first investigated the therapeutic effects of PDJQ on T2DM rats. 16S rRNA sequencing and untargeted metabolomics analyses were used to investigate the mechanisms of action of PDJQ in the treatment of T2DM. Our results showed that PDJQ treatment could improve the hyperglycemia, hyperlipidemia, insulin resistance (IR) and pathological changes of liver, pancreas, kidney, and colon in T2DM rats. PDJQ could also decrease the levels of pro-inflammatory cytokines and inhibit the oxidative stress. 16S rRNA sequencing showed that PDJQ could decrease the Firmicutes/Bacteroidetes (F to B) ratio at the phylum level. At the genus level, PDJQ could increase the relative abundances of Lactobacillus, Blautia, Bacteroides, Desulfovibrio and Akkermansia and decrease the relative abundance of Prevotella. Serum untargeted metabolomics analysis showed that PDJQ could regulate tryptophan metabolism, histidine metabolism, tricarboxylic acid (TCA) cycle, phenylalanine, tyrosine and tryptophan biosynthesis and tyrosine metabolism pathways. Correlation analysis indicated that the modulatory effects of PDJQ on the tryptophan metabolism, histidine metabolism and TCA cycle pathways were related to alterations in the abundance of Lactobacillus, Bacteroides and Akkermansia. In conclusion, our study revealed the various ameliorative effects of PDJQ on T2DM, including improving the liver and kidney functions and alleviating the hyperglycemia, hyperlipidemia, IR, pathological changes, oxidative stress and inflammatory response. The mechanisms of PDJQ on T2DM are likely linked to an improvement in the dysbiosis of gut microbiota and modulation of tryptophan metabolism, histamine metabolism, and the TCA cycle.
Introduction
Type 2 diabetes mellitus (T2DM) is a complex endocrine disease caused by a combination of environmental and genetic factors, leading to islet failure or insufficient insulin action. According to data from the World Health Organization, 415 million patients were suffering from T2DM in 2015 and this number will increase to 642 million by 2040 (Rivera-Mancía et al., 2015). The etiology of T2DM has unfortunately not been fully elucidated, and the lack of effective treatment is a major health challenge globally.
The human intestine is a complex micro-ecosystem. The bacteria residing in the intestine are large in number and diverse, and have important roles in regulating host metabolism. Under physiological conditions, the intestinal microbiota can maintain the homeostasis of host glucose and lipid metabolism by digesting and absorbing food and producing metabolites (Ortega et al., 2020). Unhealthy dietary habits can destroy certain species of intestinal bacteria and affect the body’s metabolism, which may eventually lead to the occurrence of T2DM. Compared with the status in healthy people, the abundances of many probiotic bacteria including Lactobacillus and Akkermansia are decreased in the intestine of T2DM patients, while some harmful bacteria such as Enterobacteriaceae are increased (Tanase et al., 2020). Furthermore, the abnormal levels of metabolites caused by intestinal microbiota disorder are closely related to inflammatory reactions, insulin resistance (IR), as well as the disorders of glucose and lipid metabolism (Chávez-Carbajal et al., 2020). Probiotic interventions, including with Roseburia, Bacteroides, Akkermansia, and Lactobacillus, can help improve glucose metabolism disorders and restore insulin sensitivity (Huda et al., 2021). Study have demonstrated that metformin could regulate the abundances of Bacteroidetes, Escherichia, and other microbiota, thus maintaining the integrity of the intestinal barrier, promoting the production of short-chain fatty acids (SCFAs), and regulating bile acid metabolism, thereby showing hypoglycemic effects (Zhang and Hu, 2020).
Traditional Chinese medicine (TCM) has been used to treat T2DM for thousands of years. Meta-analysis indicated that TCM could reduce the risk of progression to T2DM and increase the possibility of regression toward normoglycemia (Pang et al., 2017). It is important to investigate the mechanism of action of TCM in the treatment of T2DM, which can help not only to reveal the implications of its action, but also provide a scientific basis for new drug discovery. Regulation of the intestinal microbiota has become an important target of TCM for the treatment of T2DM. Isomaltulose can increase the abundance of Actinobacteria, Faecalibacterium, Blautia, and Phascolarctobacterium, while it decreases the abundance of Patescibacteria in the gastrointestinal (GI) tract of T2DM rats, thereby regulating the metabolism of lipids, carbohydrates, and bile acids (Yang et al., 2021). Oxymatrine has been shown to help reduce the abundances of Phascolarctobacterium, Roseburia, and Prevotella microbiota in the GI tract and increase the abundances of Bacteroidetes, Ruminococcus, and Streptococcus microbiota to regulate glucose metabolism, lipid metabolism, and carbohydrate metabolism, and to improve pathoglycemia and dyslipidemia in T2DM rats (Shao et al., 2020). Ganoderma lucidum polysaccharides could reduce the abundances of harmful bacteria, such as Ruminococcus, Corynebactrium, and Proteus, in the GI tract of T2DM rats, increase the abundances of Blautia, Dehalobacterium, Parabacteroides, and Bacteroides, and restore the disordered metabolism of amino acids and carbohydrates in T2DM rats (Chen et al., 2020). By increasing the abundance of Parabacteroides, Blautia, and Akkermansia microbiota and decreasing the abundance of Aerococcus, Staphylococcus, and Corynebacterium microbiota in the GI tract, Huang-Lian-Jie-Du decoction helped improve lipid metabolism and bile acid metabolism in T2DM rats, and improve IR (Chen et al., 2018).
Pi-Dan-Jian-Qing decoction (PDJQ), composed of Astragalus mongholicus Bunge, Pseudostellaria heterophylla (Miq.) Pax, Atractylodes lancea (Thunb.) DC., Scrophularia ningpoensis Hemsl., Coptis chinensis Franch., Scutellaria baicalensis Georgi, Pueraria montana var. lobata (Willd.) Maesen & S.M.Almeida ex Sanjappa & Predeep, Potentilla discolor Bunge, Litchi chinensis Sonn. and Salvia miltiorrhiza Bunge, has been used to treat T2DM in Yunnan Provincial Hospital of Chinese Medicine, Kunming Municipal Hospital of Traditional Chinese Medicine and Jiaxing Hospital of Traditional Chinese Medicine for more than 10 years. To date, over 1000 T2DM patients have received the treatment of PDJQ. According to our clinical study, treatment with PDJQ in combination with metformin significantly decreased the fasting blood glucose (FBG) and postprandial blood glucose levels and improved the clinical outcomes such as polydipsia and fatigue in T2DM patients compared with metformin treatment (Xie et al., in-press). However, the mechanism of PDJQ on T2DM remains unclear. In this study, we first investigated the therapeutic effect of PDJQ on T2DM rats. Then, we used 16S rRNA sequencing and untargeted metabolomics analysis to investigate the mechanism of action of PDJQ in the treatment of T2DM.
Materials And Methods
Reagents
High-fat high-carbohydrate (HFHC) diet (65.75% basal chow, 20% sucrose, 10% lard, 3% egg yolk powder, 1% cholesterol, 0.25% pig bile salt) was obtained from Beijing Sibeifu Bioscience Co., Ltd. (Beijing, China). Streptozotocin (STZ) and metformin were purchased from Solarbio Biotechnology Co., Ltd. (Beijing, China). Triglyceride (TG), total cholesterol (TC), low density lipoprotein (LDL), high density lipoprotein (HDL), aspartate aminotransferase (AST), alanine aminotransferase (ALT), creatinine (Cr), blood urea nitrogen (BUN), superoxide dismutase (SOD), methane dicarboxylic aldehyde (MDA), and glutathione peroxidase (GSH-Px) assay kits test kits were obtained from Nanjing Jiancheng Biological Engineering Institute (Nanjing, China). Rat insulin, interleukin (IL)-1β, IL-6, tumor necrosis factor alpha (TNF-α) enzyme-linked immunosorbent assay (ELISA) kit was obtained from Multi Science Biotechnology Co., Ltd. (Hangzhou, China). Reference standards of astragaloside, heterophyllin B, atractylodin, harpagide, berberine, baicalin, puerarin, quercetin, protocatechuic acid and tanshinone II A were obtained from Sichuan Weikeqi Biological Technology Co., Ltd. (Sichuan, China).
Preparation of PDJQ
For the preparation of PDJQ, the following components were used:12 g of Astragalus mongholicus Bunge, 12 g of Peudostellaria heterophylla (Miq.) Pax, 6 g of Atractylodes lancea (Thunb.) DC., 12 g of Scrophularia ningpoensis Hemsl., 12 g of Coptis chinensis Franch., 12 g of Scutellaria baicalensis Georgi, 12 g of Pueraria montana var. lobata (Willd.) Maesen & S.M.Almeida ex Sanjappa & Predeep, 12 g of Potentilla discolor Bunge, 12 g of Litchi chinensis Sonn., and 12 g of Salvia miltiorrhiza Bunge. Eight volumes of water were added to this mixture, followed by decocting for 30 min and concentrating it to 6 g of crude drug/mL.
Quality control of PDJQ was performed using ultra performance liquid chromatography (UPLC; ACQUITY UPLC®, United States) coupled with Xevo G2 quadrupole time-of-flight (Q-TOF) mass spectrometer (MS; Waters Corp., Milford, United States) systems. Briefly, the test solution was injected onto an ACQUITY UPLC BEH C18 column (2.1mm×100mm, 1.7μm) held at 50°C. The flow rate was 0.3 mL/min and the injection volumn was 2 μL. Mobile phase A was 0.1% formic acid aqueous solution and mobile phase B was acetonitrile contained 0.1% formic acid. The mobile phase conditions were: 0 min, 5% B; 1 min, 10% B; 6 min 60% B; 6.5 min 100% B; 10 min 100% B; 10.1 min 5% B; 13 min 5% B.
A Q-TOF MS equipped with an electrospray ionization (ESI) source was used for both positive and negative ionization scan modes (m/z ranges from 50 to 1,200 Da). The scan time was 0.2 s. The capillary voltages were 3,000 V (positive mode) and 2,200 V (negative mode) respectively. The desolvation temperature was 350°C and the source temperature was 100°C. The sample cone voltage was 40 V and the extraction cone voltage was 4V. The cone gas flow was 40 L/h and the desolvation gas flow was 800 L/h (both positive and negative modes).
Grouping and Dosing
Sixty 6–8-week-old specific pathogen free (SPF)-grade Sprague-Dawley (SD) male rats (180–220 g) were purchased from Beijing Huafukang Biotechnology Co., Ltd. (License number: SCXK 2020-0006). The feeding environment was maintained at 21±2°C with relative humidity of 45±10%. During the 12 h light-dark cycle, the rats had free access to water and food. The experiment was approved by the Ethics Committee of Tianjin University of Traditional Chinese Medicine.
After 1 week of acclimation, 10 rats were randomly selected as the control group and fed with normal chow, while the remaining 50 rats were fed with HFHC diet for 8 weeks. All rats were subsequently fasted for 12 h and given a single intraperitoneal injection of STZ at 30 mg/kg (dissolved in 0.1 mol/L citric acid buffer, pH = 4.5). The control group was injected with the same dose of citric acid buffer. Blood was collected from the tail vein 72 h after STZ injection to measure the random blood glucose level, and the model was successfully replicated with blood glucose > 16.7 mmol/L (Yang et al., 2020). After the animal model had been established, the T2DM rats were randomly divided into, T2DM, Metformin, PDJQ low-dose, PDJQ middle-dose and PDJQ high-dose groups with 10 rats in each group. After the normal clinical dose was converted based on the body surface area, a ×6 dose was used in PDJQ middle-dose group, while the low dose was half the middle dose and high dose was double the middle dose. Animals in the control and T2DM groups were gavaged with 2 mL of saline; rats in the Metformin group were gavaged with 0.2 g/kg metformin; while rats in the PDJQ low-, middle-, and high-dose groups were gavaged with 4.9 g, 9.8 g, and 19.6 g of crude drug/kg, respectively, once a day for 4 weeks. The FBG was measured in each group every week after PDJQ treatment (Figure 1).
Figure 1 Overview of the experimental design for all groups.
Oral Glucose Tolerance Test (OGTT)
After 4 weeks of intervention with PDJQ, all rats were fasted without water for 8 h. The blood glucose level was measured with blood glucose testing strips, which was used as the blood glucose value at 0 min. Immediately thereafter, each rat was gavaged with 50% glucose solution (2 g/kg). The blood glucose values were measured after 30, 60, 90 and 120 min of glucose solution administration. The blood glucose-time curve was plotted and the area under the curve (AUC) of OGTT was calculated.
Serum Biochemical Indicator Tests
After 4 weeks of intervention with PDJQ, all rats were fasted for 8 h. The rats were anesthetized by the intraperitoneal injection of pentobarbital sodium (50 mg/kg). Blood was collected from the abdominal aorta using a syringe and the collected blood was centrifuged at 3000 rpm/min for 15 min to collect the serum. The serum levels of lipid-related indicators (TG, TC, HDL, LDL), liver function-related indicators (ALT, AST), kidney function-related indicators (Cr, BUN), oxidative stress-related indicators (SOD, GSH-Px, MDA) were measured using kits based on the instructions provided (Nanjing Jiancheng Biological Engineering Institute).
ELISA
After 4 weeks of intervention with PDJQ, the serum levels of insulin, IL-6, IL-1β, and TNF-α in each group were measured using ELISA, which was performed based on the manufacturer’s instructions (Multi Science Biotechnology Co., Ltd.). Homeostatic model assessment of insulin resistance (HOMA-IR) was calculated using the following formula: HOMA-IR = [fasting blood glucose × fasting insulin (FINS)]/22.5.
Pathology
After 4 weeks of intervention with PDJQ, the rat liver, pancreas, kidney and colon tissues were fixed with formalin solution, paraffin-embedded, and prepared into 3-μm tissue sections, which were routinely stained with hematoxylin and eosin (H&E). The histopathological changes of each were observed under a light microscope.
16S rRNA Sequencing
Extraction of Fecal Genomic DNA
Control, T2DM and PDJQ high-dose groups were selected for intestinal microbiota study. Briefly, six rats were randomly selected from each group and their cecum contents were collected and weighed. The total genomic DNA was extracted by the CTAB/SDS method, and the concentrations of DNA samples were detected using 1% agarose gel. The DNA was diluted to 1 ng/µL with sterile water based on the concentration.
Polymerase Chain Reaction (PCR) Amplification and Sequencing of the 16S rRNA Gene
The V3–V4 region of the 16S rRNA gene was amplified using 338F (5′-ACTCCTACGGGAGGCAGCAG-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′) primers. The PCR amplification system included the following: 10 ng of template DNA, 0.2 µM forward and reverse primers, and 15 µL of Phusion® High-Fidelity PCR Master Mix (New England Biolabs). The reaction conditions were as follows: pre-denaturation at 98°C for 1 min; 15 cycles of denaturation at 95°C for 10 s, annealing at 50°C for 30 s, and extension at 72°C for 30 s; and then maintenance of the reaction at 72°C for 5 min; followed by storage at 4°C. The mixed PCR products were purified using Qiagen Gel Extraction Kit (Qiagen, Germany). Then, 2% agarose gel electrophoresis was used for detection. Sequencing libraries were generated using TruSeq® DNA PCR-Free Sample Preparation Kit (Illumina, USA). The library quality was assessed using a Qubit® 2.0 Fluorometer (Thermo Scientific) and Agilent Bioanalyzer 2100 system. Finally, the libraries were sequenced to by 250 bp of paired-end sequencing on the Illumina NovaSeq platform.
Sequencing Data Analysis
The raw sequencing data were spliced using FLASH (v1.2.7, http://ccb.jhu.edu/software/FLASH/) (Magoč and Salzberg, 2011) and sequences were quality-controlled to obtain the final effective tags. The tags were clustered using Uparse software (Version 7.0.1001, http://drive5.com/uparse/) (Edgar, 2013) to obtain the operational taxonomic units (OTUs). The Mothur-algorithm-based SILVA reference database (Version 138, http://www.arb-silva.de/) (Quast et al., 2013) was used to annotate the OTUs with taxonomic information. Multiple sequence alignments were performed using MUSCLE software (Version 3.8.31, http://www.drive5.com/muscle/) (Edgar, 2004). OTU abundance information was normalized by the sequence number corresponding to the sample with the smallest sequence. Alpha diversity index and beta diversity analysis were subsequently performed. The Wilcoxon rank-sum test was used to test for differences between groups in diversity indices, and the Kruskal–Wallis test was used (the post hoc test was chosen as Games-Howell) combined with the multiple testing method to screen for differential microbiota, where p < 0.05 after false discovery rate (FDR) correction was considered to indicate a significant difference. Finally, Phylogenetic Investigation of Communities by Reconstruction of Unobserved States database (PICRUSt) analysis was applied to predict the relevant gene pathways that may be affected by each strain of differential microbiota.
Metabolomic Analysis
Serum Sample Processing
Control, T2DM and PDJQ high-dose groups were selected for metabolomics study. Briefly, six rats were randomly selected from each group and their serum samples were collected. A total of 100 μL of serum sample was added to 400 μL of 80% methanol solution, vortexed and oscillated, left to stand for 5 min in an ice bath, and then centrifuged for 20 min (15,000 g, 4°C). After centrifugation, the supernatant was diluted with ultrapure water to 53 % methanol and then centrifuged at 15,000 g for 20 min at 4°C. The supernatant was collected for 20 min and used as the testing sample. Quality control (QC) sample was prepared by mixing all samples in equal amounts in order to ensure the stability and accuracy of the measurement throughout the analysis.
Chromatography and Mass Spectrometry Conditions
The chromatography was performed on a Hypersil Goldcolumn (C18) column (2.1 mm × 100 mm, 1.9 μm) with a mobile phase consisting of (A) 0.1% formic acid and (B) methanol. Gradient elution was used, as follows: 0 min, 98% A, 2% B; at 1.5 min, 98% A, 2% B; at 12 min, 0% A, 100% B; at 14 min, 0% A, 100% B; at 14.1 min, 98% A, 2% B; and at 17 min, 98% A, 2% B. The column temperature was set at 40°C, the flow rate was 0.2 mL/min, and the injection volume was 2 μL. The mass spectrometry conditions were as follows: simultaneous detection in positive and negative ion modes using an ESI source. The ESI source settings were as follows: spray voltage: 3.2 kV; sheath gasflow rate: 40 arb; aux gasflow rate: 10 arb; capillary temp: 320°C. QC samples were injected every six samples throughout the analytical run, and the data obtained were used to evaluate stability.
Data Processing and Analysis
Chromatographic characteristic peaks in the samples were detected using high-resolution mass spectrometry. The raw files obtained by MS were imported into Compound Discoverer 3.1 (CD3.1, Thermo Fisher) software for data pre-processing. First, the data were briefly screened by retention time and mass-to-charge ratio parameters, and then the peaks were aligned according to the retention time deviation of 0.2 min and mass deviation (parts per million, ppm) of 5 ppm for different samples to make the identification more accurate; subsequently, the data were matched according to the set values of 5 ppm, signal intensity deviation of 30%, signal-to-noise ratio (S/N) of 3, minimum signal intensity of 100,000, and adduct ions for peak extraction, and also for the quantification of peak area. From molecular ion peaks and fragment ions, the prediction of the molecular formula was performed and compared to the mzCloud (https://www.mzcloud.org/), mzVault, and MassList databases to identify the metabolites. Metabolites with a coefficient of variance of less than 30% (Dai et al., 2017) in QC samples were then retained as final identification results for subsequent analysis. The peaks detected in the samples were integrated using CD3.1, where the peak area of each characteristic peak represented the relative quantitative value of a metabolite, and the quantitative results were normalized using the total peak area, and finally the quantitative results of the metabolites were obtained. The data were then subjected to QC to ensure the accuracy and reliability of the results. Next, the metabolites were subjected to multivariate statistical analysis, including principal component analysis (PCA) and orthogonal partial least squares discriminant analysis (OPLS-DA), to reveal differences in metabolic patterns between different groups. Variable importance of projection (VIP) of each metabolites were obtained based on the OPLS-DA model. The normalized peak areas for each metabolite were also analyzed using two-tailed Student’s t-test. Differential metabolites among each group were screened based on p < 0.05 and, where VIP > 1. Metabolic pathway enrichment analysis was performed for differential metabolites with fold change (FC) > 1.25 or FC < 0.80 based on MetaboAnalyst software (https://www.metaboanalyst.ca/) and Kyoto Encyclopedia of Genes and Genomes (KEGG) database (https://www.kegg.jp/).
Statistics
All data were reported as the mean ± standard deviation (mean ± SD) for the independent experiments. Statistical differences between the experimental groups were accessed by one-way analysis of variance (ANOVA) and post-hoc analysis using SPSS software (version 20.0). A p-value < 0.05 was considered statistically significant. Curve fitting was performed using the GraphPad Prism 5 software.
Results
Identification of Main Compounds in PDJQ by UPLC-MS Analysis
Astragaloside, heterophyllin B, atractylodin, harpagide, berberine, baicalin, puerarin, quercetin, protocatechuic acid and tanshinone II A were used as the reference standards to validate the main compounds in PDJQ. The detailed information of these compounds were shown in Figure S1. The typical based peak intensity (BPI) chromatograms of PDJQ and these reference standards were shown in Figure S2. The characteristic fragment ions of these compounds were shown in Table S2 (Supplementary material). Astragaloside in ragalus mongholicus Bunge, heterophyllin B in Pseudostellaria heterophylla (Miq.) Pax, atractylodin in Atractylodes lancea (Thunb.) DC., harpagide in Scrophularia ningpoensis Hemsl., berberine in Coptis chinensis Franch., baicalin in Scutellaria baicalensis Georgi, puerarin in Pueraria montana var. lobata (Willd.) Maesen & S.M.Almeida ex Sanjappa & Predeep, quercetin in Potentilla discolor Bunge, protocatechuic acid in Litchi chinensis Sonn. and tanshinone II A in Salvia miltiorrhiza Bunge were identified as the preeminent compounds in PDJQ.
Effects of PDJQ on Blood Glucose, Blood Lipids, Liver and Kidney Functions, IR, and Pathological Manifestations in T2DM Rats
The FBG of rats was on average > 16.7 at 72 h after STZ injection. After 4 weeks of PDJQ treatment, the levels of FBG in the Metformin, PDJQ low-dose, PDJQ middle-dose, and PDJQ high-dose groups were significantly decreased compared with those of the T2DM group (p <0.01, respectively, Figure 2A).
Figure 2 PDJQ treatment alleviated the hyperglycemia, IR, pathological changes and inflammatory response in T2DM rats. (A) PDJQ treatment decreased the FBG in T2DM rats. (B, C) The AUC of OGTT was decreased in T2DM rats after metformin and PDJQ treatment. (D) PDJQ treatment significantly improved the pathological changes of liver, pancreas, kidney, and colon in T2DM rats (200×). (E) PDJQ treatment decreased the levels of pro-inflammatory cytokines in serum. Control, T2DM, Metformin, PDJQ low-dose, PDJQ middle-dose and PDJQ high-dose (n = 10 per group) groups. Data are presented as the mean ± SD. ##p < 0.01 as compared to the control group; **p < 0.01 as compared to the T2DM group; △p < 0.05 as compared to the Metformin group; △△p < 0.01 as compared to the Metformin group.
The serum levels of TG, TC, and LDL in T2DM group were significantly increased (p < 0.01, respectively), and the HDL level was significantly decreased, compared with the values in the control group (p < 0.05). The levels of TC, TG and LDL were decreased in Metformin (p < 0.01 respectively), PDJQ middle-dose (p < 0.05, p < 0.05 and p < 0.01, respectively) and PDJQ high-dose groups (p < 0.01, respectively) compared with those in T2DM group. The level of HDL was increased in Metformin (p < 0.05), PDJQ middle-dose(p < 0.05) and PDJQ high-dose (p < 0.01) groups compared with the T2DM group (Table 1).
Table 1 Changes in blood lipid levels after PDJQ treatment.
In addition, the activities of serum ALT and AST were significantly increased in T2DM rats compared with rats in the Control group (p < 0.01, respectively), and the activities of serum ALT and AST were significantly decreased in Metformin (p < 0.01, respectively), PDJQ middle-dose (p < 0.05 and p < 0.01, respectively) and PDJQ high-dose groups (p < 0.01, respectively) compared with those in T2DM group (Table 2). The serum levels of Cr and BUN were higher in the T2DM group as compared to the Control group (p < 0.01, respectively), and the serum levels of Cr and BUN were lower in Metformin (p < 0.01, respectively), PDJQ low-dose (p < 0.05, respectively), PDJQ middle-dose (p < 0.01, respectively) and PDJQ high-dose groups (p < 0.01, respectively) compared with rats in the T2DM group (Table 2).
Table 2 Changes in serum AST, ALT activities and serum Cr, BUN levels after PDJQ treatment.
OGTT showed that the AUC of OGTT was significantly higher in the T2DM group than that in Control group (p < 0.01), and the AUC of OGTT was lower in Metformin, PDJQ low-dose, PDJQ middle-dose and PDJQ high-dose groups compared with rats in the T2DM group (p < 0.01, respectively, Figures 2B, C). In addition, compared with those in the Control group, the FINS level and HOMA-IR indicator were significantly increased in the T2DM group (p < 0.01, respectively). The FINS level was reduced in Metformin, PDJQ middle-dose and PDJQ high-dose groups compared with the T2DM group (p < 0.01, respectively, Table 3). The HOMA-IR was lower in Metformin, PDJQ low-dose, PDJQ middle-dose and PDJQ high-dose groups compared with the T2DM group (p < 0.01, respectively, Table 3).
Table 3 Levels of FINS and HOMA-IR after PDJQ treatment.
H&E staining showed that rats in T2DM group exhibited significant hepatocyte steatosis, hepatocyte necrosis and infiltration of inflammatory cell in liver, islet atrophy and islet cell injury in pancreas, hyperplasia of glomerular basement membrane, renal tubular atrophy, and massive inflammatory cell infiltration in kidney, damage of epithelial cells and incomplete cell arrangement in colon. The intervention of metformin and PDJQ significantly improved the histopathological changes of liver, pancreas, kidney, and colon in T2DM rats (Figure 2D).
Effects of PDJQ on Oxidative Stress and Inflammatory Factors in T2DM Rats
We further evaluated the effects of PDJQ on oxidative stress in T2DM rats by measuring SOD, MDA, and GSH-Px in the serum of each group of rats. Compared with the Control group, the activities of SOD (p < 0.01) and GSH-Px (p < 0.01) in rats’ serum were significantly decreased and the level of MDA (p < 0.01) was significantly increased in the T2DM group. The activities of SOD and GSH-Px were elevated in Metformin (p < 0.01, respectively), PDJQ middle-dose (p < 0.05, respectively) and PDJQ high-dose (p < 0.01, respectively) groups compared with the T2DM group, whereas the level of MDA was reduced in Metformin (p < 0.01), PDJQ middle-dose (p < 0.05) and PDJQ high-dose (p < 0.01) groups compared with the T2DM group (Table 4). Furthermore, the levels of pro-inflammatory factors IL-6, IL-1β, and TNF-α in serum were measured to study the effects of PDJQ on the inflammatory response in T2DM rats. Compared with the control group, the levels of IL-6, IL-1β, and TNF-α in the serum of T2DM rats were significantly increased (p < 0.01, respectively). The level of IL-6 was lower in Metformin, PDJQ low-dose, PDJQ middle-dose and PDJQ high-dose groups compared with the T2DM group (p < 0.01, respectively). The level of IL-1β was lower in Metformin (p < 0.01) and PDJQ high-dose (p < 0.05) groups compared with the T2DM group. The level of TNF-α was lower in Metformin (p < 0.01), PDJQ middle-dose (p < 0.05) and PDJQ high-dose (p < 0.01) groups compared with the T2DM group (Figure 2E). The above results demonstrated that PDJQ had therapeutic effects on T2DM rats, which were most significant at high doses. Therefore, PDJQ high-dose group was selected for the follow-up intestinal microbiota and metabolomic studies. Rats selected for 16S rRNA sequencing and metabolomics analysis showed statistical significances in blood lipid levels (TG, TC, HDL, LDL, Table S2), liver (ALT, AST) and kidney (Cr, BUN) function-related indicators (Table S3), FBG (Figure S3A), AUC of OGTT (Figure S3B, C), FINS and HOMA-IR (Table S4) between Control and T2DM groups and between T2DM and PDJQ high-dose groups.
Table 4 Changes in serum SOD, GSH-Px activities and serum MDA level after PDJQ treatment.
Effects of PDJQ on the Intestinal Microbiota of T2DM Rats
16S rRNA sequencing was used to investigate the changes of intestinal microbiota in T2DM rats after PDJQ treatment. The results showed that the Shannon and Simpson indexes were significantly higher in the T2DM group than in the Control group (p < 0.05) and significantly lower in the PDJQ high-dose group than in the T2DM group (p < 0.05, Figures 3A, B). The principal co-ordinates analysis (PCoA) results showed that sample points of the T2DM group could be completely separated from the Control group, and the sample points of the PDJQ high-dose group were very close to those of the Control group. Likewise, clustering analysis also demonstrated that the distance from the Control group to the PDJQ high-dose group was shorter than the distance from the Control group to the T2DM group (Figures 3C, D). The above results indicated that the overall structure and composition of the intestinal microbiota of T2DM rats had changed significantly, and the high dose of PDJQ could effectively reverse this change.
Figure 3 PDJQ treatment affected the gut microbiota community in T2DM rats. (A, B) Shannon and simpson indexes were lower in PDJQ high-dose group than that in the T2DM group. (C, D) PCoA and system clustering tree indicated more similar beta diversity between PDJQ high-dose and Control groups than that between the T2DM and Control groups (C: Control group; M: T2DM group; P: PDJQ high-dose group). (E, F) At the phylum level, PDJQ treatment decreased the F to B ratio in T2DM rats. (G) At the genus level, PDJQ treatment affected the relative abundances of Lactobacillus, Blautia, Prevotella, Bacteroides, Desulfovibrio and Akkermansia in T2DM rats. (H, I) The differential metabolic pathways (written in red) of PDJQ on T2DM were predicted using PICRUSt analysis based on the 16S rRNA sequencing data. Control, T2DM and PDJQ high-dose (n = 6 per group) groups. Data are presented as the mean ± SD. #p < 0.05 as compared to the control group; ##p < 0.01 as compared to the control group; *p < 0.05 as compared to the T2DM group; **p < 0.01 as compared to the T2DM group.
The composition of the intestinal microbiota in each group of samples at the phylum level is shown in the Figure 3E, showing Firmicutes and Bacteroidetes as the dominant ones. The Firmicutes/Bacteroidetes (F to B) ratio was significantly higher in the T2DM group than in the Control group (p < 0.01), while it was significantly lower in the PDJQ high-dose group than in the T2DM group (p < 0.05, Figure 3F). At the genus level, the abundances of Allobaculum (p < 0.01) and Enterococcus (p < 0.05) were higher, while those of Lactobacillus (p < 0.01), Blautia (p < 0.01), Bacteroides (p < 0.01), and Akkermansia (p < 0.01) were lower in the T2DM group than in the Control group. Compared with the T2DM group, the abundances of Lactobacillus (p < 0.05), Blautia (p < 0.01), Bacteroides (p < 0.05), Desulfovibrio (p < 0.05), and Akkermansia (p < 0.05)were increased, while that of Prevotella (p < 0.01) was decreased in PDJQ high-dose group (Figure 3G).
PICRUSt analysis was used to predict the functional changes of intestinal microbiota in T2DM rats after the intervention with a PDJQ. The results of the prediction based on the KEGG database showed that a total of six categories of biological metabolic pathway functional analyses were obtained at the first functional layer, which were Cellular Processes, Environmental Information Processing, Genetic Information Processing, Human Diseases, Metabolism, and Organismal Systems. The top ten metabolic pathways with the highest proportions and p-values < 0.05 are listed in the Figure 3H (Control group vs. T2DM group) and Figure 3I (T2DM group vs. PDJQ high-dose group). Proportions of metabolic pathways that were increased in the T2DM group compared with Control group but decreased in PDJQ high-dose group compared with T2DM group, or vice versa, were considered as differential pathways. The abundances of glycolysis/gluconeogenesis and tricarboxylic acid (TCA) cycle pathways were all increased in the T2DM group compared with the Control group. The abundances of purine metabolism, tryptophan metabolism and histidine metabolism pathways were lower in T2DM group compared with the Control group (Figure 3H). For the PDJQ high-dose group, the abundances of purine metabolism, trptophan metabolism and histidine metabolism pathways were increased and the abundances of glycolysis/gluconeogenesis and TCA cycle pathways were decreased compared with the T2DM group (Figure 3I).
Effects of PDJQ on the Levels of Serum Metabolites in T2DM Rats
PCA generates new characteristic variables using linear combinations of metabolite variables with certain weights, and groups of data are categorized by major new variables (principal components). The PCA model can reflect the original state of metabolomic data, and the degree of aggregation and dispersion of the samples can be observed from the PCA model. The PCA plot showed that the Control group could be well distinguished from the T2DM group, and the T2DM group and the PDJQ high-dose group were also well distinguished (Figure 4A). To identify the differential metabolites, the OPLS-DA model was used. In agreement with the PCA results, OPLS-DA models showed significant distinctions of metabolomic data between the Control and the T2DM group as well as between the T2DM group and the PDJQ high-dose group (Figures 4B, D). In addition, the explanatory rate (R2) and predictive power (Q2) of OPLA-DA model were evaluated using seven-round cross validation and 200 repetitions of response permutation testing (RPT). Compared with the Control group, the T2DM group had R2 = 0.88 and Q2 = −0.81 (Figure 4C). Compared with the T2DM group, the PDJQ high-dose group had R2 = 0.89 and Q2 = −0.69 (Figure 4E). These results indicated that the OPLS-DA model was stable and had good predictive ability.
Figure 4 PDJQ treatment regulated the serum metabolites in T2DM rats. (A) Scores plots of PCA between the control and T2DM groups and the T2DM and PDJQ high-dose groups. (B, C) Scores plots of OPLS−DA between the Control and T2DM groups and the corresponding coefficient of loading plots. (D, E) Scores plots of OPLS−DA between the T2DM and PDJQ high-dose groups and the corresponding coefficient of loading plots. (F, G) Summary of pathway analysis of serum samples between control and T2DM groups and between T2DM and PDJQ high-dose groups. (A)Tryptophan metabolism; (B) Histidine metabolism; (C) Glycerolipid metabolism; (D) Glyoxylate and dicarboxylate metabolism; (E) TCA cycle; (F) Phenylalanine, tyrosine and tryptophan biosynthesis; (G) Tyrosine metabolism. Control, T2DM and PDJQ high-dose (n = 6 per group) groups.
The detailed trends of the changes of differential metabolites (FC > 1.2 or FC < 0.8, p < 0.05 and VIP > 1.0) are shown in Table 5. We also used MetaboAnalyst software to screen the differential metabolic pathways to screen for significant metabolic pathways (p < 0.05, impact value > 0.10) of the differential metabolites. Compared with the control group, the metabolic pathways that were altered in the T2DM group mainly included the metabolism of tryptophan, histidine, glycerolipid, glyoxylate, and dicarboxylate and the TCA cycle (Figure 4F); while the metabolic pathways affected by PDJQ mainly included tryptophan metabolism, histidine metabolism, the TCA cycle, phenylalanine, tyrosine, and tryptophan biosynthesis, and tyrosine metabolism (Figure 4G). The same metabolic pathway of PICRUSt analysis that had done 16S rRNA sequencing and untargeted metabolomics pathway analysis included tryptophan metabolism, histidine metabolism, and the TCA cycle and were discussed in detail.
Table 5 The differential metabolites in serum after PDJQ treatment.
Correlation Analysis of Physiological Data, Oxidative Stress and Inflammatory Factors, Untargeted Metabolomics and Gut Microbiota
Spearman’s correlation analysis was conducted to analyze the relationship between physiological data, oxidative stress and inflammatory factors, differential serum metabolites and gut microbiota at genus level in the Control, T2DM and PDJQ high-dose groups. As shown in Figure 5A, Lactobacillus, Bacteroides and Akkermansia showed negative correlations and Allobaculum, Prevotella and Enterococcus showed positive correlations with most of the physiological indices. In addition, Blautia, Bacteroides and Akkermansia showed negative correlations and Prevotella, Allobaculum and Enterococcus showed positive correlations with some of the pro-inflamatory cytokines. Desulfovibrio, Prevotella Akkermansia, Lactobacillus, Allobaculum and Enterococcus showed correlations with some of the oxidative stress factors (Figure 5B). Allobaculum, Lactobacillus, Bacteroides and Akkermansia showed correlations with most of the metabolites (Figure 5C).
Figure 5 Correlation analysis of physiological data, oxidative stress and inflammatory factors, untargeted metabolomics and 16S rRNA sequencing. Correlations between physiological data and gut microbiota (A), between oxidative stress and inflammatory factors and gut microbiota (B) and between untargeted metabolomics and gut microbiota (C) were analyzed using spearman’s analysis (heatmap). X-axis represents the physiological indices (A), oxidative stress and inflammatory factors (B) and differential metabolites in the serum (C). Y-axis represents the gut microbiota with differential abundance. The colors of grids represent the correlation analysis value of spearman’s correlation analysis. Grids in red indicate positive correlations (correlation analysis value more than 0.1), while grids in blue indicate negative correlations (correlation analysis value less than -0.1). Color coding scale indicates the correlation analysis value from heatmap, the deeper red or blue indicates higher correlation values. *p < 0.05 between physiological data, oxidative stress and inflammatory factors, differential serum metabolites and gut microbiota. **p < 0.01 between physiological data, oxidative stress and inflammatory factors, differential serum metabolites and gut microbiota.
Discussion
Our results showed that, compared with the control group, the T2DM rats exhibited a significant increase in blood glucose, disorders in lipid levels and elevated ALT, AST, Cr and BUN. In addition, the HOMA-IR and AUC of OGTT were elevated in the T2DM rats, suggesting IR in rats from the T2DM group. The pathological results showed that the hepatocytes in the liver of the T2DM group rats were disordered, enlarged, and had deformation of vacuolated balloon; the pancreatic islets in the pancreas were sparsely distributed and vacuolated; the smooth muscle layer of the colon was thinned and the number of myocytes was reduced; the glomeruli in the kidney tissue were swollen and enlarged; the epithelial cells of the renal tubules were swollen or detached; the renal interstitium was infiltrated with inflammatory cells; and mesangial proliferation and interstitial fibrosis were also observed. The above results are consistent with the pathological manifestations of T2DM (Huang et al., 2014; Fan et al., 2018). The intervention of PDJQ can reduce blood glucose and improve dyslipidemia, IR, and pathological changes of kidney, liver, pancreas, and colon tissues in T2DM rats. This suggests that PDJQ has a therapeutic effect on T2DM, with the most significant findings being found in the high-dose group. In addition, metformin was selected as a positive control for T2DM in our study. Metformin is the first-line regimen for the treatment of T2DM in clinical practice (Inzucchi et al., 2012). Our results showed that there were no significant differences in FBG, FINS, HOMA-IR, and serum levels of biochemical indicators between PDJQ high-dose and Metformin groups. The small differences between the two experimental groups suggested that PDJQ has the potential to substitute for metformin in the treatment of T2DM.
Our results also showed that PDJQ could increase SOD and GSH-Px activities and decrease MDA content in T2DM rats, suggesting that this decoction helped reduce the oxidative stress response in T2DM rats. Oxidative stress is closely related to the development of T2DM (Said et al., 2020) and the development of T2DM is characterized by disorders of glucolipid metabolism, hyperglycemia, and hyperlipidemia, which lead to the production of large amounts of reactive oxygen species (ROS) in mitochondria; these in turn damage mitochondrial functions and cause oxidative stress responses. This then leads to IR by blocking insulin action pathways, and leading to the development of T2DM (Yaribeygi et al., 2020). Oxidative stress also leads to the massive production of free radicals, causing lipid peroxidation, which in turn causes cellular damage (Chen et al., 2013). MDA is a product of lipid peroxidation caused by free radicals or ROS in cells under oxidative stress, and MDA levels can indirectly reflect the degree of oxidative damage in cells (Dai et al., 2011). SOD and GSH-Px are antioxidant enzymes that reflect the strength of antioxidant capacity (Gao et al., 2018). SOD, as an intracellular oxygen radical scavenger, catalyzes the conversion of O2− to O2 and H2O2, and protects the organism from superoxide anions (Ding et al., 2019). GSH-Px is an important peroxide-degrading enzyme that is widespread in the body. It catalyzes the conversion of reduced glutathione to oxidized glutathione and protects cells from peroxide disruption and damage (Tsai et al., 2019).
The results of ELISA showed that PDJQ reduced the levels of the pro-inflammatory factors IL-6, IL-1β, and TNF-α levels in sera and that the inflammatory reaction is a key factor behind the development of T2DM. The production of local and systemic pro-inflammatory cytokines such as IL-6, IL-1β, and TNF-α exacerbates the vicious cycle of T2DM by inducing IR (Baek et al., 2020). IL-6 decreases insulin activity by blocking the expression of insulin receptor signaling, leading to the development of IR (Senn et al., 2002). IL-1β contributes to increase the secretion of intercellular adhesion molecule-1 (Yang et al., 2010), which is responsible for islet beta-cell injury and death, further exacerbating diabetes. TNF-α can induce the phosphorylation of insulin receptor substrate-1 (Rotter et al., 2003), thus disrupt the metabolism of blood glucose.
16S rRNA sequencing was used to study the effects of PDJQ on the composition of intestinal microbiota in T2DM rats. The alpha diversity of intestinal microbiota refers to the diversity of flora in a specific region or ecosystem, and is a comprehensive indicator reflecting the richness and homogeneity of flora. Shannon and Simpson indexes were significantly increased in T2DM model rats, indicating that the alpha diversity of gut microbiota was increased in T2DM model rats. PDJQ treatment decreased the alpha diversity of gut microbiota. The beta diversity of intestinal microbiota of rats was analyzed by PCoA and cluster analysis, and the overall composition of intestinal microbiota in T2DM rats changed significantly. The intestinal microbiota in T2DM rats was restored to the normal level after the intervention of PDJQ. Our results are consistent with previous reports demonstrating that T2DM could alter the alpha and beta diversities of intestinal microbiota (Chen et al., 2018).
Moreover, PDJQ helps to reduce the high F to B ratio caused by T2DM. This ratio is closely linked to many diseases and significantly increased in many models, for instance, T2DM, non-alcoholic fatty liver disease, and depression (Cheung et al., 2019; Jasirwan et al., 2021; Xia et al., 2021). Metabolic disorder and inflammatory reaction can be alleviated by lowering the ratio of F to B ratio (Stojanov et al., 2020). The relative abundances of Allobaculum and Enterococcus were significantly increased, while those of Lactobacillus, Blautia, Bacteroides, and Akkermansia were significantly decreased in the T2DM rat model; PDJQ significantly increased the relative abundances of Lactobacillus, Blautia, Bacteroides, and Desulfovibrio, and Akkermansia and decreased the abundance of Prevotella. Lactobacillus is a probiotic that improves intestinal ecological dysregulation and oxidative stress in T2DM rats (Morshedi et al., 2020). Blautia is associated with anti-inflammatory effects along with many diseases including colorectal carcinoma, hepatic encephalopathy, and graft versus host disease (Jenq et al., 2015). It was found that a decrease in the intestinal probiotic Bacteroides was associated with abnormal glucose tolerance, and had an impact on the sugar and lipid metabolism in the body, which in turn contributed to the development of obesity and the onset of diabetes mellitus (Qiao et al., 2020). Further studies have shown that these processes are mainly related to the production of SCFAs, which are produced by intestinal microbiota breaking down oligosaccharides, polysaccharides, peptides, proteins, and glycoproteins; these in turn directly regulate the number and function of pancreatic beta-cells, and thus participate in hepatic glycogen metabolism (Liu et al., 2020). The increased abundance of Akkermansia reduces the elevated lipopolysaccharide levels caused by a high-fat diet and is important for maintaining glucose homeostasis (Shin et al., 2014). It has been found that metformin treatment could significantly increase the proportion of Akkermansia in obese mice (Lee and Ko, 2014). Likewise, our correlation analysis also showed negative correlations of Lactobacillus, Blautia, Bacteroides, Desulfovibrio and Akkermansia with some of the physiological indicies (FBG, IR indicies, blood lipid levels and blood biochemical indicators), oxidative stress factor (MDA) and pro-inflammatory cytokines(IL-1β, IL-6, TNF-α). Positive correlations were also found between some of the anti-oxidative factors (SOD, GSH-Px) and Lactobacillus, Desulfovibrio. Finally, Prevotella, a pathogenic bacterium specialized in surviving under anaerobic conditions, is usually found in oral infections and is capable of causing wound infection and purulence (Yanagisawa et al., 2006). Correlation analysis also showed positive correlations of Prevotella with physiological indicies (FBG, AUC of OGTT, AST, Cr, BUN, HOMA-IR), pro-inflammatory cytokine (IL-1β) and oxidative stress factor (MDA). However, the role of Prevotella during T2DM progression still need to be further studied. Further studies should also be carried out using fecal transplantation or gut microbiota depletion models to verify whether PDJQ can ameliorate hyperglycemia, hyperlipidemia, IR, oxidative stress and inflammatory responses to treat T2DM through regulating these gut microbiota.
The results of serum untargeted metabolomics showed that PDJQ could affect tryptophan metabolism, histidine metabolism, TCA cycle, phenylalanine, tyrosine and tryptophan biosynthesis, and tyrosine metabolism. After correlating the differential metabolic pathways obtained from metabolomics with those deduced from 16S rRNA sequencing, the pathways of tryptophan metabolism, histidine metabolism, and the TCA cycle were identified as the shared pathways, suggesting that PDJQ may affect tryptophan metabolism, histamine metabolism, and the TCA cycle through the regulation of intestinal microbiota to treat T2DM.
Tryptophan Metabolism
Tryptophan and its related metabolites are closely associated with T2DM and IR (Chen et al., 2016). In our study, we found that the levels of DL-tryptophan and L-kynurenine were reduced in T2DM model rats and significantly increased after treatment with PDJQ. The DL-tryptophan metabolite quinolinic acid inhibits the activity of phosphoenolpyruvate carboxykinase (PEPC), which in turn inhibits gluconeogenesis and lowers blood glucose levels (Dayer et al., 2009). The majority of tryptophan in the body is metabolized by the kynurenine pathway, which is involved in inflammatory and immune reactions (Zhou and Ye, 2018). The key enzymes for the conversion of tryptophan to kynurenine are tryptophan 2,3-dioxygenase (TDO) and indoleamine 2,3-dioxygenase (IDO), which are activated by stress hormones or inflammatory factors and then slowly overproduce downstream metabolites including kynurenine. This ultimately causes the onset of diabetes (Kartika et al., 2020). The ratio of kynurenine to tryptophan is an indicator of IDO activation. When the ratio of tryptophan to kynurenine is high, the organism is more prone to inflammatory reactions and autoimmune disease (Brandacher et al., 2006). Spearman’s analysis showed a positive correlation of Lactobacillus, Bacteroides and Akkermansia with DL-tryptophan and L-kynurenine levels. Additionally, we observed a negative correlation between Allobaculum and L-kynurenine. As 16S r RNA sequencing showed that there were no significant difference in the relative abundance of Allobaculum between T2DM and PDJQ high-dose groups, it’s likely the modulatory effects of PDJQ on tryptophan metabolism may occur through affecting the abundances of Lactobacillus, Bacteroides and Akkermansia.
Histamine Metabolism
Histidine and its metabolites are closely related to the resistance to oxidative stress and inflammatory reactions in the body (Yu et al., 2015). Our results showed that the level of L-histidine was decreased in T2DM rats and increased after the treatment with PDJQ. Histidine can inhibit TNF-α signaling pathway and attenuate the expression of pro-inflammatory cytokines (Lee et al., 2005), and it has been used to treat inflammation in obese patients with metabolic syndrome (Feng et al., 2013). In addition, histidine shows strong binding affinity for Fe2+ ions, thus reducing the amount of ROS produced by the Fenton reaction and protecting cells from damage caused by iron overload (Vera-Aviles et al., 2018; Wake et al., 2020). Histidine enhances the expression and activity of catalase and GSH-Px, thus reducing oxidative stress (Liu et al., 2008; Watanabe et al., 2008). Spearman’s analysis showed a positive correlation of Lactobacillus, Bacteroides and Akkermansia with L-histidine level. Thus, we assumed that the effects of PDJQ on L-histidine level may be associated with modulating the abundances of Lactobacillus, Bacteroides and Akkermansia.
TCA Cycle
The TCA cycle is the final metabolic pathway for carbohydrates, lipids, and amino acids (Chen et al., 2014). The TCA cycle can affect the immunity of the body (Choi et al., 2012). and is closely associated with many different metabolic diseases including T2DM. In our study, we found that D-(+)-malic acid, citric acid, and fumaric acid were significantly elevated in the T2DM group of rats and decreased after administration of the medicine. Citric acid is an important intermediate in the TCA cycle, which has been found to be an important substrate for cellular energy metabolism and is closely associated with T2DM. It is a major source of adenosine triphosphate (ATP) produced by cells and has key roles in inflammation, insulin secretion, and metabolic regulation (Iacobazzi and Infantino, 2014). Citric acid can be cleaved by ATP-citrate synthase (ACLY) to acetyl coenzyme A and oxaloacetate (OAA) (Wang et al., 2021). Acetyl coenzyme A provides a basis for the synthesis of lipids, including arachidonic acid, which produces inflammatory mediators such as prostaglandins, leukotrienes, and thromboxanes (Douzi et al., 2020). OAA is the second product of citric acid and can beconverted to malic acid, which is in turn oxidized by malic enzyme to pyruvic acid to form nicotinamide adenine dinucleotide phosphate (NADPH), which can produce ROS by the action of NADPH oxidase, and the excess ROS can lead to oxidative stress (Munhoz et al., 2016). Metabolic disorder of fumaric acid may be associated with diabetic renal impairment. Fumaric acid is a dicarboxylic acid that is a precursor of L-malate during the TCA cycle. The accumulation of fumaric acid has been shown to lead to oxidative stress (Zheng et al., 2015), and prolonged oxidative stress leads to renal impairment in diabetic nephropathy. Negative correlations between Bacteroides, Akkermansia and D-(+)-malic acid levels and between Bacteroides, Akkermansia Lactobacillus and fumaric acid were found by our spearman’s analysis. Allobaculum showed positive correlations with fumaric acid level. The effect of PDJQ on TCA cycle in T2DM rats could be through affecting the abundances of Bacteroides, Akkermansia Lactobacillus.
Conclusion
In conclusion, our study revealed the various ameliorative effects of PDJQ on T2DM, including improving the liver and kidney functions and alleviating the hyperglycemia, hyperlipidemia, IR, pathological changes, oxidative stress and inflammatory responses. The mechanisms of PDJQ on T2DM are likely linked to an improvement in the dysbiosis of gut microbiota and modulation of tryptophan metabolism, histamine metabolism, and the TCA cycle.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/bioproject/PRJEB48123.
Ethics Statement
The animal study was reviewed and approved by Yunnan Hospital of Traditional Chinese Medicine Ethics Committee.
Author Contributions
XX wrote the manuscript. XX, JL, YA, JG, and JZ conducted animal experiments. XX, JL, FQ, and CX finished molecular bioassays. ZZ, WW, HC, and HW provided technical guidance for the whole work. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by National Science Foundation of China (No. 82104802); Science and Technology Projects of Yunnan Province (No. 202101AT070246); Key scientific project of Yunnan province (No. 2019ZF005).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2021.748872/full#supplementary-material
References
Baek, S. M., Kim, K., Kim, S., Son, Y., Hong, H. S., Yu, S. Y. (2020). SP Prevents T2DM Complications by Immunomodulation. Sci. Rep. 10, 1. doi: 10.1038/s41598-020-73994-1
Brandacher, G., Winkler, C., Aigner, F., Schwelberger, H., Schroecksnadel, K., Margreiter, R., et al. (2006). Bariatric Surgery Cannot Prevent Tryptophan Depletion Due to Chronic Immune Activation in Morbidly Obese Patients. Obes. Surg. 16:5. doi: 10.1381/096089206776945066
Chávez-Carbajal, A., Pizano-Zárate, M. L., Hernández-Quiroz, F., Ortiz-Luna, G. F., Morales-Hernández, R. M., De Sales-Millán, A., et al. (2020). Characterization of the Gut Microbiota of Individuals at Different T2D Stages Reveals a Complex Relationship With the Host. Microorganisms. 8, 1. doi: 10.3390/microorganisms8010094
Chen, S. K., Hsu, C. H., Tsai, M. L., Chen, R. H., Drummen, G. P. (2013). Inhibition of Oxidative Stress by Low-Molecular-Weight Polysaccharides With Various Functional Groups in Skin Fibroblasts. Int. J. Mol. Sci. 14, 10. doi: 10.3390/ijms141019399
Chen, M., Liao, Z., Lu, B., Wang, M., Lin, L., Zhang, S., et al. (2018). Huang-Lian-Jie-Du-Decoction Ameliorates Hyperglycemia and Insulin Resistant in Association With Gut Microbiota Modulation. Front. Microbiol. 9, 2380. doi: 10.3389/fmicb.2018.02380
Chen, L., Liu, T., Zhou, J., Wang, Y., Wang, X., Di, W., et al. (2014). Citrate Synthase Expression Affects Tumor Phenotype and Drug Resistance in Human Ovarian Carcinoma. PloS One 9, 12. doi: 10.1371/journal.pone.0115708
Chen, M., Xiao, D., Liu, W., Song, Y., Zou, B., Li, L., et al. (2020). Intake of Ganoderma Lucidum Polysaccharides Reverses the Disturbed Gut Microbiota and Metabolism in Type 2 Diabetic Rats. Int. J. Biol. Macromol. 155, 890–902. doi: 10.1016/j.ijbiomac.2019.11.047
Chen, T., Zheng, X., Ma, X., Bao, Y., Ni, Y., Hu, C., et al. (2016). Tryptophan Predicts the Risk for Future Type 2 Diabetes. PLoS One 11, e0162192. doi: 10.1371/journal.pone.0162192
Cheung, S. G., Goldenthal, A. R., Uhlemann, A. C., Mann, J. J., Miller, J. M., Sublette, M. E. (2019). Systematic Review of Gut Microbiota and Major Depression. Front. Psychiatry 10, 34. doi: 10.3389/fpsyt.2019.00034
Choi, I., Son, H., Baek, J. H. (2012). Tricarboxylic Acid (TCA) Cycle Intermediates: Regulators of Immune Responses. Life. 11, 1. doi: 10.3390/life11010069
Dai, J., Ma, H., Fan, J., Li, Y., Wang, J., Ni, H., et al. (2011). Crude Polysaccharide From an Anti-UVB Cell Clone of Bupleurum Scorzonerifolium Protect HaCaT Cells Against UVB-Induced Oxidative Stress. Cytotechnology 63, 6. doi: 10.1007/s10616-011-9381-6
Dai, W., Xie, D., Lu, M., Li, P., Lv, H., Yang, C., et al. (2017). Characterization of White Tea Metabolome: Comparison Against Green and Black Tea by a Nontargeted Metabolomics Approach. Food. Res. Int. 96, 40–45. doi: 10.1016/j.foodres.2017.03.028
Dayer, M. R., Safari, I., Dayer, M. S. (2009). New Evidence on Hypoglycemic Effect of Quinolinic Acid in Diabetic Rats. Pak. J. Biol. Sci. 12, 14. doi: 10.3923/pjbs.2009.1025.1030
Ding, X., Jian, T., Wu, Y., Zuo, Y., Li, J., Lv, H., et al. (2019). Ellagic Acid Ameliorates Oxidative Stress and Insulin Resistance in High Glucose-Treated HepG2 Cells via miR-223/Keap1-Nrf2 Pathway. Biomed Pharmacother 110, 85–94. doi: 10.1016/j.biopha.2018.11.018
Douzi, W., Guillot, X., Bon, D., Seguin, F., Boildieu, N., Wendling, D., et al. (2020). 1h-NMR-Based Analysis for Exploring Knee Synovial Fluid Metabolite Changes After Local Cryotherapy in Knee Arthritis Patients. Metabolites. 10, 11. doi: 10.3390/metabo10110460
Edgar, R. C. (2004). MUSCLE: Multiple Sequence Alignment With High Accuracy and High Throughput. Nucleic. Acids Res. 32, 5. doi: 10.1093/nar/gkh340
Edgar, R. C. (2013). UPARSE: Highly Accurate OTU Sequences From Microbial Amplicon Reads. Nat. Methods 10:10. doi: 10.1038/nmeth.2604
Fan, Y. B., He, Z. W., Wang, W., Li, J. J., Hu, A. M., Li, L., et al. (2018). Tangganjian Decoction Ameliorates Type 2 Diabetes Mellitus and Nonalcoholic Fatty Liver Disease in Rats by Activating the IRS/PI3K/AKT Signaling Pathway. Biomed. Pharmacother. 106, 733–737. doi: 10.1016/j.biopha.2018.06.089
Feng, R. N., Niu, Y. C., Sun, X. W., Li, Q., Zhao, C., Wang, C., et al. (2013). Histidine Supplementation Improves Insulin Resistance Through Suppressed Inflammation in Obese Women With the Metabolic Syndrome: A Randomised Controlled Trial. Diabetologia 56, 985–994. doi: 10.1007/s00125-013-2839-7
Gao, Z., Lai, Q., Yang, Q., Xu, N., Liu, W., Zhao, F., et al. (2018). The Characteristic, Antioxidative and Multiple Organ Protective of Acidic-Extractable Mycelium Polysaccharides by Pleurotus Eryngii Var. Tuoliensis on High-Fat Emulsion Induced-Hypertriglyceridemic Mice. Sci. Rep. 8, 1. doi: 10.1038/s41598-018-35706-8
Huang, S., Peng, W. F., Jiang, X. H., Shao, K., Xia, L. L., Tang, Y. B., et al. (2014). The Effect of Chromium Picolinate Supplementation on the Pancreas and Macroangiopathy in Type II Diabetes Mellitus Rats. J. Diabetes. Res. 2014:717219. doi: 10.1155/2014/717219
Huda, M. N., Kim, M., Bennett, B. J. (2021). Modulating the Microbiota as a Therapeutic Intervention for Type 2 Diabetes. Front. Endocrinol. 12, 632335. doi: 10.3389/fendo.2021.632335
Iacobazzi, V., Infantino, V. (2014). Citrate—New Functions for an Old Metabolite. Biol. Chem. 395, 387–399. doi: 10.1515/hsz-2013-0271
Inzucchi, S. E., Bergenstal, R. M., Buse, J. B., Diamant, M., Ferrannini, E., Nauck, M., et al. (2012). Management of Hyperglycemia in Type 2 Diabetes: A Patient-Centered Approach: Position Statement of the American Diabetes Association (ADA) and the European Association for the Study of Diabetes (EASD). Diabetes. Care 35, 6. doi: 10.2337/dc12-0413
Jasirwan, C. O. M., Muradi, A., Hasan, I., Simadibrata, M., Rinaldi, I. (2021). Correlation of Gut Firmicutes/Bacteroidetes Ratio With Fibrosis and Steatosis Stratified by Body Mass Index in Patients With non-Alcoholic Fatty Liver Disease. Biosci. Microbiota Food Health 40, 50–58. doi: 10.12938/bmfh.2020-046
Jenq, R. R., Taur, Y., Devlin, S. M., Ponce, D. M., Goldberg, J. D., Ahr, K. F., et al. (2015). Intestinal Blautia Is Associated With Reduced Death From Graft-Versus-Host Disease. Biol. Blood. Marrow. Transplant. 21, 8. doi: 10.1016/j.bbmt.2015.04.016
Kartika, R., Wibowo, H., Purnamasari, D., Pradipta, S., Larasati, R. A. (2020). Altered Indoleamine 2,3-Dioxygenase Production and Its Association to Inflammatory Cytokines in Peripheral Blood Mononuclear Cells Culture of Type 2 Diabetes Mellitus. Int. J. Tryptophan Res. IJTR 13, 1178646920978236. doi: 10.1177/1178646920978236
Lee, Y. T., Hsu, C. C., Lin, M. H., Liu, K. S., Yin, M. C. (2005). Histidine and Carnosine Delay Diabetic Deterioration in Mice and Protect Human Low Density Lipoprotein Against Oxidation and Glycation. Eur. J. Pharmacol. 513, 145. doi: 10.1016/j.ejphar.2005.02.010
Lee, H., Ko, G. (2014). Effect of Metformin on Metabolic Improvement and Gut Microbiota. Appl. Environ. Microbiol. 80, 19. doi: 10.1128/AEM.01357-14
Liu, W. H., Liu, T. C., Yin, M. C. (2008). Beneficial Effects of Histidine and Carnosine on Ethanol-Induced Chronic Liver Injury. Food. Chem. Toxicol. 46, 5. doi: 10.1016/j.fct.2007.12.013
Liu, J. L., Segovia, I., Yuan, X. L., Gao, Z. H. (2020). Controversial Roles of Gut Microbiota-Derived Short-Chain Fatty Acids (SCFAs) on Pancreatic β-Cell Growth and Insulin Secretion. Int. J. Mol. Sci. 21, 3. doi: 10.3390/ijms21030910
Magoč, T., Salzberg, S. L. (2011). FLASH: Fast Length Adjustment of Short Reads to Improve Genome Assemblies. Bioinformatics. 27, 21. doi: 10.1093/bioinformatics/btr507
Morshedi, M., Saghafi-Asl, M., Hosseinifard, E. S. (2020). The Potential Therapeutic Effects of the Gut Microbiome Manipulation by Synbiotic Containing-Lactobacillus Plantarum on Neuropsychological Performance of Diabetic Rats. J. Transl. Med. 18, 1. doi: 10.1186/s12967-019-02169-y
Munhoz, A. C., Riva, P., Simões, D., Curi, R., Carpinelli, A. R. (2016). Control of Insulin Secretion by Production of Reactive Oxygen Species: Study Performed in Pancreatic Islets From Fed and 48-Hour Fasted Wistar Rats. PloS. One 11, 6. doi: 10.1371/journal.pone.0158166
Ortega, M. A., Fraile-Martínez, O., Naya, I., García-Honduvilla, N., Álvarez-Mon, M., Buján, J., et al. (2020). Type 2 Diabetes Mellitus Associated With Obesity (Diabesity). The Central Role of Gut Microbiota and Its Translational Applications. Nutrients. 12, 9. doi: 10.3390/nu12092749
Pang, B., Lian, F. M., Zhao, X. Y., Zhao, X. M., Jin, D., Lin, Y. Q., et al. (2017). Prevention of Type 2 Diabetes With the Traditional Chinese Patent Medicine: A Systematic Review and Meta-Analysis. Diabetes Res. Clin. Pract. 131, 242–259. doi: 10.1016/j.diabres.2017.07.020
Qiao, S., Bao, L., Wang, K., Sun, S., Liao, M., Liu, C., et al. (2020). Activation of a Specific Gut Bacteroides-Folate-Liver Axis Benefits for the Alleviation of Nonalcoholic Hepatic Steatosis. Cell. Rep. 32, 6. doi: 10.1016/j.celrep.2020.108005
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA Ribosomal RNA Gene Database Project: Improved Data Processing and Web-Based Tools. Nucleic. Acids Res. 41, D590. doi: 10.1093/nar/gks1219
Rivera-Mancía, S., Lozada-García, M. C., Pedraza-Chaverri, J. (2015). Experimental Evidence for Curcumin and its Analogs for Management of Diabetes Mellitus and its Associated Complications. Eur. J. Pharmacol. 5, 756. doi: 10.1016/j.ejphar.2015.02.045
Rotter, V., Nagaev, I., Smith, U. (2003). Interleukin-6 (IL-6) Induces Insulin Resistance in 3T3-L1 Adipocytes and is, Like IL-8 and Tumor Necrosis Factor-Alpha, Overexpressed in Human Fat Cells From Insulin-Resistant Subjects. J. Biol. Chem. 278, 46. doi: 10.1074/jbc.M301977200
Said, E., Mousa, S., Fawzi, M., Sabry, N. A., Farid, S. (2020). Combined Effect of High-Dose Vitamin A, Vitamin E Supplementation, and Zinc on Adult Patients With Diabetes: A Randomized Trial. J. Adv. Res. 2020, 28. doi: 10.1016/j.jare.2020.06.013
Senn, J. J., Klover, P. J., Nowak, I. A., Mooney, R. A. (2002). Interleukin-6 Induces Cellular Insulin Resistance in Hepatocytes. Diabetes 51, 12. doi: 10.2337/diabetes.51.12.3391
Shao, J., Liu, Y., Wang, H., Luo, Y., Chen, L. (2020). An Integrated Fecal Microbiome and Metabolomics in T2DM Rats Reveal Antidiabetes Effects From Host-Microbial Metabolic Axis of EtOAc Extract From Sophora Flavescens. Oxid. Med. Cell. Longev. 2020, 1805418. doi: 10.1155/2020/1805418
Shin, N. R., Lee, J. C., Lee, H. Y., Kim, M. S., Whon, T. W., Lee, M. S., et al. (2014). An Increase in the Akkermansia Spp. Population Induced by Metformin Treatment Improves Glucose Homeostasis in Diet-Induced Obese Mice. Gut. 63, 5. doi: 10.1136/gutjnl-2012-303839
Stojanov, S., Berlec, A., Štrukelj, B. (2020). The Influence of Probiotics on the Firmicutes/Bacteroidetes Ratio in the Treatment of Obesity and Inflammatory Bowel Disease. Microorganisms. 8, 1715. doi: 10.3390/microorganisms8111715
Tanase, D. M., Gosav, E. M., Neculae, E., Costea, C. F., Ciocoiu, M., Hurjui, L. L., et al. (2020). Role of Gut Microbiota on Onset and Progression of Microvascular Complications of Type 2 Diabetes (T2DM). Nutrients 12, 12. doi: 10.3390/nu12123719
Tsai, J. C., Chen, Y. A., Wu, J. T., Cheng, K. C., Lai, P. S., Liu, K. F., et al. (2019). Extracts From Fermented Black Garlic Exhibit a Hepatoprotective Effect on Acute Hepatic Injury. Molecules. 24, 6. doi: 10.3390/molecules24061112
Vera-Aviles, M., Vantana, E., Kardinasari, E., Koh, N. L., Latunde-Dada, G. O. (2018). Protective Role of Histidine Supplementation Against Oxidative Stress Damage in the Management of Anemia of Chronic Kidney Disease. Pharmaceuticals. 11, 4. doi: 10.3390/ph11040111
Wake, H., Takahashi, Y., Yoshii, Y., Gao, S., Mori, S., Wang, D., et al. (2020). Histidine-Rich Glycoprotein Possesses Antioxidant Activity Through Self-Oxidation and Inhibition of Hydroxyl Radical Production via Chelating Divalent Metal Ions in Fenton’s Reaction. Free. Radic. Res. 54, 8. doi: 10.1080/10715762.2020.1825703
Wang, Y., Tang, B., Long, L., Luo, P., Xiang, W., Li, X., et al. (2021). Improvement of Obesity-Associated Disorders by a Small-Molecule Drug Targeting Mitochondria of Adipose Tissue Macrophages. Nat. Commun. 12, 1. doi: 10.1038/s41467-020-20315-9
Watanabe, M., Suliman, M. E., Qureshi, A. R., Garcia-Lopez, E., Barany, P., Heimburger, O., et al. (2008). Consequences of Low Plasma Histidine in Chronic Kidney Disease Patients: Associations With Inflammation, Oxidative Stress, and Mortality. Am. J. Clin. Nutr. 87, 6. doi: 10.1093/ajcn/87.6.1860
Xia, F., Wen, L. P., Ge, B. C., Li, Y. X., Li, F. P., Zhou, B. J. (2021). Gut Microbiota as a Target for Prevention and Treatment of Type 2 Diabetes: Mechanisms and Dietary Natural Products. World J. Diabetes. 12, 1146–1163. doi: 10.4239/wjd.v12.i8.1146
Xie, X. H., Zhao, J., Liao, J. B., Yang, L., Fang, F., Yuan, E. Z., et al. (2021). Clinical Study of Pi-Dan-Jian-Qing Decoction on Type 2 Diabetes Mellitus Patients. TMR Phamacol Res 1, 23. doi: 10.53388/TMRPR20211009023
Yanagisawa, M., Kuriyama, T., Williams, D. W., Nakagawa, K., Karasawa, T. (2006). Proteinase Activity of Prevotella Species Associated With Oral Purulent Infection. Curr. Microbiol. 52, 5. doi: 10.1007/s00284-005-0261-1
Yang, Z. D., Guo, Y. S., Huang, J. S., Gao, Y. F., Peng, F., Xu, R. Y., et al. (2021). Isomaltulose Exhibits Prebiotic Activity, and Modulates Gut Microbiota, the Production of Short Chain Fatty Acids, and Secondary Bile Acids in Rats. Molecules. 26, 9. doi: 10.3390/molecules26092464
Yang, C. M., Luo, S. F., Hsieh, H. L., Chi, P. L., Lin, C. C., Wu, C. C., et al. (2010). Interleukin-1beta Induces ICAM-1 Expression Enhancing Leukocyte Adhesion in Human Rheumatoid Arthritis Synovial Fibroblasts: Involvement of ERK, JNK, AP-1, and NF-kappaB. J. Cell. Physiol. 224, 2. doi: 10.1002/jcp.22153
Yang, L., Lin, H., Lin, W., Xu, S. (2020). Exercise Ameliorates Insulin Resistance of Type 2 Diabetes through Motivating Short-Chain Fatty Acid-Mediated Skeletal Muscle Cell Autophagy. Biology 9, 203. doi: 10.3390/biology9080203
Yaribeygi, H., Sathyapalan, T., Atkin, S. L., Sahebkar, A. (2020). Molecular Mechanisms Linking Oxidative Stress and Diabetes Mellitus. Oxid. Med. Cell. Longevity 2020, 8609213. doi: 10.1155/2020/8609213
Yu, B., Li, A. H., Muzny, D., Veeraraghavan, N., de Vries, P. S., Bis, J. C., et al. (2015). Association of Rare Loss-Of-Function Alleles in HAL, Serum Histidine: Levels and Incident Coronary Heart Disease. Circ. Cardiovasc. Genet. 8, 351–355. doi: 10.1161/CIRCGENETICS.114.000697
Zhang, Q., Hu, N. (2020). Effects of Metformin on the Gut Microbiota in Obesity and Type 2 Diabetes Mellitus. Diabetes. Metab. Syndr. Obes. 13, 5003–5014. doi: 10.2147/DMSO.S286430
Zheng, L., Cardaci, S., Jerby, L., MacKenzie, E. D., Sciacovelli, M., Johnson, T. I., et al. (2015). Fumarate Induces Redox-Dependent Senescence by Modifying Glutathione Metabolism. Nat. Commun. 6, 6001. doi: 10.1038/ncomms7001
Keywords: type 2 diabetes mellitus, Pi-Dan-Jian-Qing decoction, gut microbiota, tryptophan metabolism, histamine metabolism, TCA cycle
Citation: Xie X, Liao J, Ai Y, Gao J, Zhao J, Qu F, Xu C, Zhang Z, Wen W, Cui H and Wang H (2021) Pi-Dan-Jian-Qing Decoction Ameliorates Type 2 Diabetes Mellitus Through Regulating the Gut Microbiota and Serum Metabolism. Front. Cell. Infect. Microbiol. 11:748872. doi: 10.3389/fcimb.2021.748872
Received: 28 July 2021; Accepted: 03 November 2021;
Published: 06 December 2021.
Edited by:
Linhua Zhao, China Academy of Chinese Medical Sciences, ChinaReviewed by:
Palok Aich, National Institute of Science Education and Research (NISER), IndiaJunguk Hur, University of North Dakota, United States
Copyright © 2021 Xie, Liao, Ai, Gao, Zhao, Qu, Xu, Zhang, Wen, Cui and Wang. 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: Weibo Wen, d2Vud2VpYm8yMDIwQDE2My5jb20=; Huantian Cui, MTc2MjMxNjQxMUBxcS5jb20=; Hongwu Wang, d2h3MjAwOUB0anV0Y20uZWR1LmNu
†These authors have contributed equally to this work