Abstract
Diabetic nephropathy (DN) is a severe diabetic complication of the kidney and is the main cause of end-stage kidney disease worldwide. Cordyceps sinensis (C. sinensis) is not only a traditional Chinese medicine (TCM) but also a healthy food. In China, C. sinensis has been widely used to treat various kidney diseases. Bailing Capsule, which active ingredient is C. sinensis, is approved to treat kidney disease, respiratory disease, and immune disease. However, its underlying mechanism in DN remains unclear. The purpose of the present study was to investigate the underlying mechanism of Bailing Capsule on kidney in diabetic rats. The DN model was established by streptozotocin (STZ) injection. Low and high doses of Bailing Capsule were orally administrated for 12Ā weeks after diabetes induction. Renal function was evaluated by serum creatinine, blood urea nitrogen, 24-h urinary protein, and urinary albumin. Mesangial matrix expansion and renal fibrosis were measured using histopathology staining. We found that the disorder of renal function and pathology in DN rats was significantly modified by Bailing Capsule treatment. Consistently, Bailing Capsule markedly alleviated DN rat glomerulosclerosis, tubulointerstitial injury and renal fibrosis as shown by pathological staining. Moreover, Bailing Capsule significantly reduced the kidney triglyceride content and renal lipid droplet formation in DN rats. The renal transcriptome revealed that Bailing Capsule-treated kidneys had 498 upregulated genes and 448 downregulated genes. These differentially expressed genes were enriched in the peroxisome proliferator activated receptor (PPAR) pathway and fatty acid metabolism function ontology. mRNA and protein expression analyses revealed substantial enhancement of the lipolysis pathway and inhibition of lipogenesis in Bailing Capsule-treated rat kidneys compared to DN rats. Bailing Capsule activated the expression of PPARĪ±, ACOX1 (acyl-CoA oxidase 1), and SCD (stearoyl-CoA desaturase) in diabetic nephropathy while suppressing the expression of FASN (fatty acid synthase). In conclusion, Bailing Capsule could attenuate renal triglyceride accumulation in diabetic rats by moderating PPARĪ± pathway.
Introduction
Diabetes mellitus (DM) is a chronic metabolic disease characterized by hyperglycemia. The main mechanism of DM is insulin deficiency or insulin production and/or action disorder. Currently, DM is trending toward an increasingly severe epidemic worldwide (Zimmet et al., 2016). According to the 10th edition of the International Diabetes Federation (IDF) published in January 2022, the number of patients with diabetes globally reach 537 million in 2021, and will increase to an estimated 783 million by 2045 (Sun et al., 2022).
Poorly controlled chronic hyperglycemia leads to various diabetic complications of the microvasculature and macrovasculature. Diabetic nephropathy (DN) is severe diabetic microvascular complication. Up to 30ā40% of diabetic patients develop DN (Said et al., 2018). The main clinical characteristics are substantial proteinuria and a continuous decrease in glomerular filtration. An augmented accumulation of extracellular matrix (ECM) proteins also occurs (Xu et al., 2018). Eventually, nodular diabetic glomerulosclerosis is induced (Wang et al., 2018). DN has become the main risk factor for chronic kidney disease and end-stage renal disease (ESRD) (Sun et al., 2022) worldwide. Therefore, it is urgent to establish effective treatments to reduce its societal and financial burden.
The pathogenesis of DN is complicated and mainly includes changes in glomerular hypertrophy, mesangial matrix deposition, basement membrane thickening, tubulointerstitial fibrosis and atrophy. The mechanisms of DN are also complex, including the activation of oxidative stress (Baynes, 1991), reactive oxygen species (ROS) (Lapolla et al., 2005), transforming growth factor Ī²1 (TGF-Ī²1) (Aihara et al., 2010), inflammation (Luc et al., 2019), advanced glycation end products (AGEs) (Gross et al., 2020), apoptosis of renal tubular epithelial cells (Ding et al., 2021) and mitochondrial dysfunction (Xiao et al., 2017). However, the treatment of DN is facing challenges, such as drug side effects. Thus, natural product application may be a trend to cure DN. Traditional Chinese medicine (TCM) and healthy food have benefits in treating many diseases (Dai et al., 2020; Wu et al., 2020). In particular, TCMs usually have multiple active components that can affect multiple targets and pathways.
The main ingredient of Bailing Capsule is fermented Cordyceps sinensis (C. sinensis, syn. Hirsutella sinensis) powder. Bailing Capsule is approved by the State Food and Drug Administration of China (state medical license No. Z10910036). C. sinensis belongs to the family Clavicipitaceae, division Ascomycotina. It is a parasitic fungus. C. sinensis, which has multiple active constituents, including cordycepin and its derivatives, polysaccharides, trace elements, and mycelia. These chemical ingredients exert several physiological effects of resistance to oxidation, fibrosis, tumors, viruses, and inflammation. For centuries, C. sinensis has been widely used as a TCM. Recently, C. sinensis has been used in clinical practice to treat chronic kidney disease (CKD) (Deng, 2001). In animal models, C. sinensis reduced fasting blood glucose (Lu et al., 2019) and attenuated renal fibrosis (Du et al., 2015; Dong et al., 2019) and kidney function (Xiao et al., 2018; Lu et al., 2019).
Previous clinical research found that Bailing Capsule could reduce proteinuria and improve renal function in DN (Sheng et al., 2020), nephrotic syndrome (NS) (Xu et al., 2020), lupus nephritis patients (Ren et al., 2019) and renal transplant recipients (Wang et al., 2013). Moreover, Bailing Capsule combined with Western medicine, has better efficacy than Western medicine alone to treat kidney diseases by reducing urine protein and protecting kidney function (Hu et al., 2021; Yu et al., 2021). The mechanism research of Bailing Capsule on moderating kidney function is limited, including anti-inflammatory, anti-oxidative (Niu and Li, 2019), immune moderation (Ma and Yan, 2021).
Thus, the mechanism of the action of Bailing Capsule is still unclear and requires further investigation. Conventional pharmacological methods focused one kind or group of molecules. Whole transcriptome profiling of gene arrays is a useful tool to study complex pathways and drug targets (Logan et al., 2018) and has been widely used in the diabetes field (Zhang et al., 2019). It is important to reveal TCM in regulatory network and build a bridge between TCM and modern pharmacology at whole system level. Thus, we designed this study to determine the effect of Bailing Capsule on the kidney damage of DN rats induced by streptozotocin (STZ) and used a whole transcriptome gene array to identify the differentially expressed genes between the renal tissue of Bailing Capsule-treated group and DN group, screen the pathways and molecular targets, and explore the therapeutic mechanism of Bailing Capsule in delaying the progression of DN.
Materials and methods
Medicine
Bailing Capsule was used in our study. The composition of Bailing Capsule was fermented C. sinensis powder (Cs-C-Q80).
Ultra-performance liquid chromatography analysis of bailing capsule
We performed the chromatographic separation of Bailing Capsule by using a Waters Acquity Ultra-Performance Liquid Chromatography (UPLC, Water Corp., Milford, MA) with a symmetrical C18 column (5Ā Ī¼m, 4.6 Ć 250Ā mm, Water Corp., Milford, MA). The chromatographic parameters were as follows: mobile phase A (0.1% phosphoric acid solution) and mobile phase B (methanol), flow rate of 1.0Ā ml/min, column temperature at 40Ā°C, and detection wavelength of 260Ā nm with an injection volume of 5Ā Ī¼L. Gradient separation was based on the following: 0ā5Ā min, 100% A, 10ā15Ā min, 80% A, 15ā16Ā min, 50% A.
Animal treatments and diets
Thirty-two male Sprague-Dawley (SD) rats, aged 5Ā weeks, were purchased from the Institute of Laboratory Animal Sciences, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China. All rats were housed in a specific pathogen-free facility at a temperature of 25 Ā± 1Ā°C and a standard 12-h light/12-h dark cycle. Sterilized water and food was provided to the rats ad libitum. All of the animal procedures were approval by the Animal Care Committee of the Peking Union Medical Hospital Animal Ethics Committee (Project XHDW-2015-0051, 15 February 2015). DN rat model was developed by a high-fat diet (HFD) feeding combined with low dosage of streptozotocin (STZ) injection, according to previous studies (Danda et al., 2005; Liu et al., 2020; Guo et al., 2021). All rats were randomly divided into a normal control group (n = 8) fed a standard diet (SD, fat content, 10Ā kcal %, D12450B, Research Diets Inc., New Brunswick, NJ) and a DN model group (n = 24) fed a HFD (fat content, 60Ā kcal %, D12492, Research Diets Inc., New Brunswick, NJ), respectively. After 4Ā weeks, the DN model was induced by an intraperitoneal injection of 1% (g/100Ā ml) STZ (S0130, Sigma-Aldrich, St. Louis, MO) at a dose of 40Ā mg/kg in 0.1Ā mol/L citrate buffer, pH 4.5 (Furman, 2015; Nath et al., 2017; Liu et al., 2020). Fasting blood glucose (FBG) was assayed from the tail vein by an automatic glucometer (Contour TS, Bayer, Leverkusen, Germany) after 7Ā days. The rats with FBG >11.1Ā mmol/L were divided into the DN group (n = 8), the low dosage of Bailing Capsule (BC-L, n = 8) group, and the high dosage of Bailing Capsule (BC-H, n = 8) group. After disassembling the capsule shell of Bailing Capsule (Hangzhou Zhongmei East China Pharmaceutical Co., Hangzhou, China), the content was ground. The typical human daily dose of Bailing Capsule is 20Ā g per 60Ā kg of body weight. Thus, according to the formula drat = (37 Ć dhuman)/6, the corresponding dose of Bailing Capsule for rats is 2.06Ā g/kg per day. Therefore, the BC-L group and BC-H group were gavaged with Bailing Capsule at 1 and 2Ā g/kg/day (dissolved in saline), respectively. The gavage volume is 10Ā ml/kg body weight. The rats in the DN group and control group were given an equal volume of saline. At week 12 after treatment, all rats were sacrificed under anesthesia with an intraperitoneal injection of sodium pentobarbital (150Ā mg/kg). The kidneys were immediately collected. Figure 1 illustrates the study procedure.
FIGURE 1
Schematic of group design. SD, standard diet; HFD, high-fat diet; DN, diabetic nephropathy; BC-L, low dosage of Bailing Capsule; BC-H, high dosage of Bailing Capsule.
Measurement of body weight, kidney weight, and fasting blood glucose
At the end of the treatment, the body weight and kidney weight of the SD rats were recorded. Kidney hypertrophy was calculated as the ratio of the kidney weight to the body weight. After 6Ā h of fasting, FBG was determined using an automatic glucometer (Contour TS, Bayer, Leverkusen, Germany).
Measurements of serum creatinine, blood urea nitrogen
After 6Ā h of fasting, blood was collected through the intraorbital retrobulbar plexus. Serum creatinine (OSR6178, Beckman Coulter Inc., Brea, CA), and BUN (OSR6234, Beckman Coulter Inc., Brea, CA) were measured by an automatic biochemical analyzer (AU5800, Beckman Coulter Inc., Brea, CA).
Measurements of 24Ā h urinary protein and urinary albumin
The rats were housed in individual metabolic cages at the end of the 12-week treatment in order to collect 24Ā h urine samples. Urine protein (OSR6132, Beckman Coulter Inc., Brea, CA) and albumin (OSR6202, Beckman Coulter Inc., Brea, CA) levels were measured with a Beckman biochemical analyzer (AU5800, Beckman Coulter Inc., Brea, CA).
Histological observation of the kidney
Kidney samples were fixed in 10% neutral-buffered formalin for 24Ā h and dehydrated with graded ethanol. The specimens were embedded in paraffin, cut into slices (5Ā Ī¼m thick), dewaxed, rehydrated. Then, hematoxylin and eosin (H&E), Periodic Acid Schiff (PAS) and Massonās trichome stained samples were separately used for histopathological and collagen examination. The histological changes were observed using a digital microscope (Nikon, Tokyo, Japan). ImageJ 1.4 software (NIH, Bethesda, MA) was used for image analysis.
Triglyceride assay in the kidney
According to Folch method (Folch et al., 1957), 100Ā mg kidney were ground with PBS and extracted with chloroform and methanol. After centrifugation, the contents of TG (120Ā Ī¼L) were determined at absorption wavelength 420Ā nm by spectrophotometry method (Solarbio, Beijing, China). The TG content was computed as following, TG concentration = ConcentrationStandand Ć (Absorption ValueTestāAbsorption ValueBlank) Ć· (Absorption ValueStandardāAbsorption ValueBlank) Ć· Sample protein concentration.
Oil Red O staining of the kidney
Oil Red O staining was used to detect lipid deposition in the rat renal tissue. Frozen kidney tissue sections (5Ā Ī¼m) were fixed in 10% formalin for 10Ā min, washed in distilled water and then stained in 1% Oil Red O solution (60% isopropanol) for 15Ā min. The sections were rinsed with 60% isopropanol for 2Ā min and counterstained with hematoxylin for 1Ā min. The stained sections were visualized using a digital microscope (Nikon, Tokyo, Japan), and the area of lipid deposition was assessed using ImageJ 1.4 software (NIH, Bethesda, MA) by converting the RGB image to an 8-bit grayscale image and adjusting threshold to analyze (Mehlem et al., 2013).
RNA preparation and transcriptome microarray
Total RNA was collected from kidney tissues with an RNeasy RNA extraction kit (Qiagen, Valencia, CA). The RNA concentration, purity and integrity were determined by a Nanodropā¢ ND-1000 and 1.0% agarose gel electrophoresis. Cy3-UTP-labeled cRNA was fragmented and hybridized to an Agilent Rat gene expression 4 Ć 44Ā K microarray (Agilent Technologies, Santa Clara, CA). The microarray was incubated for 17Ā h at 65Ā°C in an Agilent Hybridization Oven. Then, the microarray was washed, fixed and scanned by using an Agilent DNA Microarray Scanner (G2505C, Agilent Technologies, Santa Clara, CA). Quantile normalization and subsequent data processing were performed using the GeneSpring GX v12.1 software package (Agilent Technologies, Santa Clara, CA).
Data normalization and analysis
Differentially expressed genes with corrected p value <0.05 and fold change >1.5 were judged as statistically significant. The data obtained from the microarray have been deposited in the National Center for Biotechnology Information Gene Expression Omnibus (GEO) database (accession number GSE198021, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE198021).
Gene function enrichment analysis
Gene Ontology (GO) biological process and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed on the above differentially expressed genes by using Database for Annotation, Visualization and Integrated Discovery Software (Huang da et al., 2009) (DAVID, version 6.8, http://david.ncifcrf.gov/). Statistical differences of enriched GO and KEGG terms were defined as p < 0.05. The top 5 and 10 most significant GO and KEGG terms were selected and plotted, respectively.
Quantitative real time polymerase chain reaction analysis
Total RNA and then cDNA were prepared from renal tissue using TRIzol reagent (Invitrogen, Carlsbad, CA) and reverse transcriptase (PrimeSciptā¢ RT reagent Kit, Takara Bio Inc., Shiga, Japan). Acox1 (acyl-CoA oxidase 1), Ppara (peroxisome proliferator activated receptor alpha), Scd (stearoyl-CoA desaturase), and Fasn (fatty acid synthase) were analyzed by qPCR analysis using SYBR Green Mix Kit (Applied Biosystems, Foster City, CA) on the ABI ViiA 7 (Applied Biosystems, Carlsbad, CA). Relative changes in gene expression were calculated using the 2-ā³ā³CT method. Ī²-actin was used as an internal control. The primers used in this study are listed in Table 1.
TABLE 1
| Gene (rat) | Gene bank ID | Forward primer (5ā² to 3ā²) | Reverse primer (3ā² to 5ā²) | Product size (bp) |
|---|---|---|---|---|
| Acox1 | NM_017340 | TGTāCTGāTCAāCTTāCTGāTCGāCC | CGGāACTāGCCāATCāCAAāGATāGT | 135 |
| Ppara | NM_013196 | CTGāTCCāGCTāACTāTCGāAGTāCC | TCAāAGGāGGAāCAAāCCAāGAGāGA | 135 |
| Scd | NM_139192 | GTTāCTTāCATāCGAāCTGāCATāGGC | GAAāCAGāGAAāCTCāAGAāAGCāCCA | 143 |
| Fasn | NM_017332 | GAGāCACāTGAāTGAāGCAāCACāCT | CCAāTCAāGGTāTTCāAGCāCCCāAT | 111 |
Primers used for qPCR.
Acox1, acyl-CoA oxidase 1; Ppara, peroxisome proliferator activated receptor alpha; Scd, stearoyl-CoA desaturase; Fasn, fatty acid synthase.
Immunohistochemistry for ACOX1, PPARĪ±, SCD, and FASN in the kidney
Immunohistochemistry staining was performed. Briefly, the sections were preincubated with 3% hydrogen peroxide and blocked with 2.5% normal horse serum. Then, the slides were incubated with anti-ACOX1 (ab184032, 1:200, Abcam, Cambridge, United Kingdom), PPARĪ± (ab233078, 1:200, Abcam, Cambridge, United Kingdom), SCD (ab19862, 1:300, Abcam, Cambridge, United Kingdom), and FASN (ab22759, 1:200, Abcam, Cambridge, United Kingdom) primary antibodies overnight at 4Ā°C and then with the secondary antibody (GB25301, 1:500, Servicebio Co., Ltd., Wuhan, China) at room temperature for 60Ā min. ImageJ 1.4 software (NIH, Bethesda, MA) was used to quantify the protein expression.
Western blot analysis of ACOX1, PPARĪ±, SCD, and FASN in the kidney
Protein was extracted from renal tissues. Membranes were incubated with primary antibodies against ACOX1 (ab184032,1:1000, Abcam, Cambridge, United Kingdom), PPARĪ± (ab126285, 1:1000, Abcam, Cambridge, United Kingdom), SCD (ab39969, 1:1000, Abcam, Cambridge, United Kingdom), FASN (ab128870, 1:10000, Abcam, Cambridge, United Kingdom), and Ī²-actin (ab8226, 1:1000, Abcam, Cambridge, United Kingdom) at 4Ā°C overnight, and then the membranes were further incubated with horseradish peroxidase-conjugated secondary antibodies (GB25301, 1:5000, Servicebio Co., Ltd., Wuhan, China). Protein band densities were detected by chemiluminescence (Cell Signaling Technology, Danvers, MA), followed by scanning using an Epson V300 scanning system (Epson, Suwa, Japan), analyzed by AlphaEaseFC software (Alpha Innotech, San Leandro, CA), and quantified by ImageJ 1.4 software (NIH, Bethesda, MA).
Statistical analysis
The results are shown as the mean Ā± SD. Comparisons between groups were assessed using one-way ANOVA followed by Tukeyās post hoc test (GraphPad Prism, v 7.0, GraphPad Software, Inc., La Jolla, CA), and the significance level was set at p values <0.05.
Results
Ultra-performance liquid chromatography of bailing capsule
Seven main components in Bailing Capsule were confirmed by UPLC analysis, including 1) adenine (390.85Ā Ī¼g/ml), 2) uracil (73.41Ā Ī¼g/ml), 3) uridine (1342.18Ā Ī¼g/ml), 4) adenosine (2136.96Ā Ī¼g/ml), 5) 2ā²-deoxyadenosine (203.36Ā Ī¼g/ml), 6) guanosine (20.25Ā Ī¼g/ml), and 7) thymidine (1.46Ā Ī¼g/ml, Figure 2).
FIGURE 2
Representative UPLC chromatograms of Bailing Capsule. 1) adenine (390.85Ā Ī¼g/ml), 2) uracil (73.41Ā Ī¼g/ml), 3) uridine (1342.18Ā Ī¼g/ml), 4) adenosine (2136.96Ā Ī¼g/ml), 5) 2ā²-deoxyadenosine (203.36Ā Ī¼g/ml), 6) guanosine (20.25Ā Ī¼g/ml), and 7) thymidine (1.46Ā Ī¼g/ml).
Bailing capsule ameliorates fasting blood glucose in diabetic nephropathy rats
FBG in the DN groups was remarkably elevated compared with the control group, while Bailing Capsule treatment significantly reduced FBG (p < 0.01, Figure 3A).
FIGURE 3
Bailing Capsule improved blood glucose and kidney function. (A) Bailing Capsule significantly reduced fasting blood glucose (FBG). (B) Bailing Capsule significantly reduced kidney hypertrophy. (C) Bailing Capsule significantly reduced 24Ā h urinary protein. (D) Bailing Capsule significantly reduced urinary albuminuria. (E) Bailing Capsule significantly reduced serum creatinine. (F) Bailing Capsule significantly reduced blood urea nitrogen (BUN). DN, diabetic nephropathy; BC-L, low dosage of Bailing Capsule; BC-H, high dosage of Bailing Capsule. Data are expressed as mean Ā± SD (n = 8). *p < 0.05 and **p < 0.01.
Bailing capsule improved kidney hypertrophy and renal function in diabetic nephropathy rats
Kidney hypertrophy is expressed as a ratio of kidney weight to body weight. DN rats showed a significant increase in kidney hypertrophy compared with the control group (p < 0.01, Figure 3B), while the administration of Bailing Capsule to DN rats significantly reduced kidney hypertrophy (p < 0.01, Figure 3B). As shown in Figures 3CāF, DN rats exhibited elevated 24Ā h urinary protein, urinary albuminuria, serum creatinine, and BUN, compared to control rats (p < 0.01, Figures 3CāF). Bailing Capsule treatment significantly reduced 24Ā h urinary protein, urinary albuminuria, serum creatinine, and BUN (p < 0.01, Figures 3CāF).
Bailing capsule moderates histopathological changes in diabetic nephropathy rats
In Figures 4A,B,D, H&E and PAS staining show significant morphological changes in DN rat kidney, including glomerular hypertrophy, glomerular mesangial matrix accumulation, glomerular basement membrane thickening and tubular basement membrane thickening (p < 0.01). Bailing Capsule treatment alleviated these changes (p < 0.01, Figures 4A,B,D). Masson staining showed that glomerular and tubulointerstitial fibrosis occurred in DN rat kidney (p < 0.01, Figures 4C,E). However, the degree of fibrosis was relieved after administration of Bailing Capsule (p < 0.01, Figures 4C,E).
FIGURE 4
Bailing Capsule alleviated morphological changes of kidney. (A) H&E staining (400 X) of kidney showed morphological changes in DN rat kidney, including glomerular basement membrane thickening and tubular basement membrane thickening. Bailing Capsule treatment alleviated these changes. Red arrow showed the basement membrane. (B) PAS staining (400 X) of kidney showed morphological changes in DN rat kidney, including glomerular mesangial matrix accumulation. Bailing Capsule treatment alleviated these changes. Blue arrow showed glomerular mesangial matrix accumulation. (C) Masson staining (400 X) of kidney showed collagen fiber deposition. Bailing Capsule relieved renal fibrosis in DN rats. The blue stain by the black arrow indicated is collagens fiber. Scale bar = 20Ā Ī¼m. (D) The quantification of mesangial area/glomerular area (%) of kidney. (E) The quantification of collagen area (%) of kidney. DN, diabetic nephropathy; BC-L, low dosage of Bailing Capsule; BC-H, high dosage of Bailing Capsule. Data are expressed as mean Ā± SD (n = 8). *p < 0.05 and **p < 0.01.
Bailing capsule administration reduces kidney lipid accumulation in diabetic nephropathy rats
Kidney TG contents were increased significantly in DN rats compared with those in the control group (p < 0.01, Figure 5A), and this change was blocked in Bailing Capsule-treated groups (p < 0.01, Figure 5A). Moreover, enhanced lipid deposition was observed in DN kidneys, compared to control kidneys, as assessed by Oil Red O staining (Figures 5B,C). Bailing Capsule reduced the percentage of lipid deposits in the kidney of diabetic rats (p < 0.01, Figures 5B,C).
FIGURE 5
Bailing Capsule reduced kidney lipid accumulation. (A) Kidney triglyceride (TG) content. Bailing Capsule reduced kidney TG contents in DN rats. (B) Lipid droplets were detected by Oil Red O staining (400 X). Bailing Capsule reduced the percentage of lipid deposits in DN rats. Black arrow showed the lipid droplets. Scale bar = 20Ā Ī¼m. (C) The quantification of percentage of lipid deposits. DN, diabetic nephropathy; BC-L, low dosage of Bailing Capsule; BC-H, high dosage of Bailing Capsule. Data are expressed as mean Ā± SD (n = 8). **p < 0.01.
Renal transcriptome analysis
To clarify the underlying mechanisms of Bailing Capsule on renal protection, we applied transcriptomic analysis to compare the kidneys of the BC-H group and DN group (n = 3). In this study, a total of 946 differentially expressed genes (DEGs) were identified in the BC-H group versus the DN group, 498 of which were upregulated and 448 of which were downregulated (fold change >1.5, p < 0.05).
Gene function enrichment analysis
GO function and KEGG pathway enrichment analyses of differentially expressed genes were performed. The top 5 enriched GO biological processes in each catalog included lipid metabolic process, response to drug, fatty acid metabolic process, response to xenobiotic stimulus, inflammatory response in biological processes, extracellular space, cell surface, apical plasma membrane, extracellular region, basolateral plasma membrane in cellular components, identical protein binding, protein homodimerization activity, integrin binding, receptor binding, and carboxylic ester hydrolase activity in molecular function (p < 0.001, Table 2, Figure 6A). The top 10 enriched KEGG pathways included complement and coagulation cascades, PPAR signaling pathway, cholesterol metabolism, metabolic pathways, bile secretion, biosynthesis of unsaturated fatty acids, peroxisome, phosphoinositide 3-kinase (PI3K)-Akt signaling pathway, chemical carcinogenesis, and metabolism of xenobiotics by cytochrome P450 (p < 0.01, Table 3, Figure 6B). All DEGs in the PPAR pathway are shown in Figure 7.
TABLE 2
| Term ID | Term name | Count | p value | Involved genes | Fold enrichment | FDR | Catalog |
|---|---|---|---|---|---|---|---|
| GO:0006629 | lipid metabolic process | 19 | 2.42 Ć 10ā8 | Slc16a1, Ptprn2, Acsl1, Mogat2, Plcl1, Mttp, Pon1, Cidea, Trpv1, Fabp5, Aadac, Acox1, Apoc2, Dpep2, Hacl1, Apoe, Angptl4, Apob, Enpp6 | 5.16 | 8.14 Ć 10ā5 | biological processes |
| GO:0042493 | response to drug | 40 | 1.72 Ć 10ā6 | Alas2, Cdkn1a, Mgst1, Lpl, Wfdc1, Ptprm, Pdgfa, Slco1a6, Hsd11b2, Mat2a, Oxct1, Erbb4, Cd36, Hmgcs2, Abcb1a, Abcb1b, Tgif1, Ntrk1, Abcc2, Cox8b, Gstm1, Acsl1, Itga3, Mgmt, Apoc3, Mthfr, Dusp6, Gnao1, Hadha, Sfrp2, Ccne1, Gsta3, Apoc2, Acot2, Cyp1a1, Txnip, Pemt, Plin2, Cftr, Tp73 | 2.32 | 0.0027 | biological processes |
| GO:0006631 | fatty acid metabolic process | 15 | 2.49 Ć 10ā6 | Acbd7, Faah, Acot7, Acsl1, Acot12, Lpl, Baat, Acadsb, Cyp1a1, Acot2, Hacl1, Cd36, Gpat2, Ppara, Acot4 | 4.81 | 0.0027 | biological processes |
| GO:0009410 | response to xenobiotic stimulus | 39 | 6.99 Ć 10ā6 | Alas2, Cdkn1a, Mgst1, Lpl, Wfdc1, Ptprm, Pdgfa, Slco1a6, Hsd11b2, Mat2a, Oxct1, Erbb4, Cd36, Hmgcs2, Abcb1a, Abcb1b, Tgif1, Ntrk1, Abcc2, Cox8b, Gstm1, Acsl1, Itga3, Mgmt, Apoc3, Mthfr, Dusp6, Gnao1, Hadha, Sfrp2, Ccne1, Apoc2, Acot2, Cyp1a1, Txnip, Pemt, Plin2, Cftr, Tp73 | 2.22 | 0.0058 | biological processes |
| GO:0006954 | inflammatory response | 27 | 2.17 Ć 10ā5 | Pxk, Calca, Serpina1, Sgms1, Fpr2, Cxcl2, Kng1, C4a, Hyal1, Hrh4, Pdpn, Abcc2, Igfbp4, Rarres2, Apoc3, Trpv1, Il17re, Tirap, Cxcl11, Vnn1, Fasn, Prkd1, Bmpr1b, Klkb1, S100a8, Epha2, Tp73 | 2.56 | 0.014 | biological processes |
| GO:0005615 | extracellular space | 100 | 7.11 Ć 10ā9 | Tfrc, Cyp4a1, Clu, C4a, Efemp2, Efemp1, Tnn, Igfbp1, Cpa1, Igfbp5, Igfbp4, Sema6d, Fibcd1, Defb11, Bcan, Rbp4, Sfrp2, Col4a1, S100a9, S100a8, Calca, Pdgfa, Lpl, Kng1, Klk7, Loc290595, Apoe, Apob, Cga, Fgf21, Ambp, Lrrn2, Rarres2, Inhbb, Igf1, Cxcl11, Bmp1, Tg, Artn, Scgb3a2, Col9a1, Lrrn1, Lgals3bp, Pigr, Serpina1, Wnt2b, Cpxm2, Wfdc1, Cpz, Cxcl2, Areg, Serpina4, Ces2c, Hyal1, Tctn1, Emilin1, Mbl1, Fgb, Fga, Myoc, Vwf, Nrg1, Apoa4, Rt1-Ce5, Nrg2, F2, F3, Serpinb5, Grem2, Ceacam1, Tff3, Angptl4, Pon3, Pon1, Sema3b, Tfpi, Defb42, Cpn1, Il1rl1, Wnt11, Slit3, Slit2, Gc, S100g, Angpt1, Gdf15, Apoc3, Colq, Apln, Kitlg, Epgn, Fabp5, Ces1d, Cfhr2, Ces1e, Apoc2, Ghrl, Ptx3, Ccl28, Klkb1 | 1.80 | 3.52 Ć 10ā6 | cellular components |
| GO:0009986 | cell surface | 53 | 4.98 Ć 10ā8 | Tfrc, Itgb2, Pkd2l1, Slc4a1, Slc2a4, Rt1-Db1, Clu, Areg, Grm7, Tnn, Itgb8, Itgb7, Cd36, Scn5a, Kcnh1, Kcnh2, Fgb, Fga, Gpr37, Entpd2, Abcc2, Fcer1g, Rgd1565355, Itga3, Sema6d, Apoa4, Rt1-Ce5, F3, Bcan, Vcan, Ceacam1, Tnmd, Rt1-Bb, Adtrp, Cryab, Cftr, Epha2, Kcne1, Pdgfa, Lpl, Tfpi, Il1rl1, Slamf9, Gpc3, Mbp, Kcnn2, Apoe, Slit2, Abcb1a, Ntrk1, Ambp, Bst2, Fgfr2 | 2.28 | 1.23 Ć 10ā5 | cellular components |
| GO:0016324 | apical plasma membrane | 35 | 6.01 Ć 10ā7 | Slc26a2, Kcne1, Slc22a2, Slc2a2, Atp12a, Slc7a1, Slc4a5, Slc9a3r2, Gnat3, Pdpn, Ceacam20, Cd36, Abcb1a, Abcb1b, Mpdz, Slc39a4, S100g, Slc16a1, Abcc2, Rgd1565355, Anxa13, Lhfpl5, Nrg1, Pard3b, Bst2, Atp4b, Slc6a6, Ceacam1, Trpv6, Cdhr2, Slc26a4, Rhcg, Atp6v0d2, Slc26a6, Cftr | 2.62 | 9.91 Ć 10ā5 | cellular components |
| GO:0005576 | extracellular region | 57 | 1.74 Ć 10ā6 | Apol9a, Tfrc, Clu, C4a, Esm1, Smpd5, Eva1c, Enpp6, Adamts6, Tgm2, Ccdc39, Vwf, Igfbp5, Il1r2, Nrg1, Nrg2, Bcan, Vcan, Ccdc80, Sfrp2, Col4a1, Tff3, Angptl4, Cryab, Zp4, Calca, Dpt, Itih1, Tfpi, Ltbp1, Kng1, Coch, Loc290595, Wnt11, Spock2, Lcn10, Apoe, Slit2, Cga, Gstm1, Gdf15, Rarres2, Gzmb, Apoc3, Isg15, Inhbb, Apln, Efna1, Igsf10, Kitlg, Nmu, Scgb3a2, Ghrl, Fgf13, Ccl28, Scrg1, Klkb1 | 1.97 | 2.16 Ć 10ā4 | cellular components |
| GO:0016323 | basolateral plasma membrane | 26 | 2.91 Ć 10ā6 | Tfrc, Slc22a2, Mttp, Slc2a2, Slc4a1, Atp12a, Aqp3, Slc7a1, Slc4a5, Nkd2, Erbb4, Pdpn, S100g, Slc16a1, Slc10a1, Itga3, Anxa13, Slc16a12, Cask, Ank2, Slc6a6, Lrrc7, Slco1a1, Rhcg, Slc26a6, Cftr | 2.95 | 2.88 Ć 10ā4 | cellular components |
| GO:0042802 | identical protein binding | 95 | 1.20 Ć 10ā4 | Clic6, Faah, Serpina1, Tfrc, Fmr1, Lect2, Ptprm, Pkd2l1, Grik2, Dnph1, Aqp3, Fgfrl1, Syne1, Gls, Grm7, Creb3l3, Tnn, Tnfaip8l1, Emilin1, Hmgcs2, Sult4a1, Mbl1, Tgm2, Kcnh1, Kcnh2, Acot7, Entpd2, Fcer1g, Vwf, Acot12, Apoa4, Mat1a, Bhlha15, Sdcbp2, Tirap, Grem2, Ceacam1, Prpf3, Padi3, Tff3, Hprt1, Prkd1, Angptl4, Rhcg, Cryab, Dctpp1, Fbp2, Zp4, Prps2, Calca, Mgst1, Pdgfa, Fhit, Cryz, Sav1, Slc9a3r2, Mat2a, Oxct1, Prdx1, Hacl1, Apoe, Slit2, Gbp2, Skil, Cap2, Atp6v1c2, Cd276, E2f8, Ntrk1, Gstm2, Slc16a1, Gstm1, Angpt1, Plk1, Mx1, Trpv1, Acadsb, Esr1, Tm6sf2, Apln, Ski, Dazl, Bst2, Bmp1, Tg, Fabp5, Trpv6, Sp1, Apoc2, Fasn, Ptx3, Mcu, Fgfr2, Slc26a6, Tp73 | 1.46 | 0.053 | molecular function |
| GO:0042803 | protein homodimerization activity | 45 | 1.25 Ć 10ā4 | Prps2, Pon3, Tfrc, Fmr1, Pon1, Lpl, Cib2, Pdgfa, Pdxp, Slc4a1, Dnph1, Syne1, Add2, Efemp2, Creb3l3, Oxct1, Erbb4, Kcnn2, Apoe, Slit2, Ect2, Gbp2, Pdlim4, Mbl1, Kcnh2, Ntrk1, Gstm2, Cadm3, Acot7, Gstm1, Ambp, Gdf15, Glce, Cidea, Inhbb, Apoa4, Sdcbp2, Bst2, Tpst2, Ceacam1, Adrb3, Sp1, Acox1, Cryab, Fgfr2 | 1.83 | 0.053 | molecular function |
| GO:0005178 | integrin binding | 16 | 1.60 Ć 10ā4 | Egfl6, Syk, Vwf, Itga3, Itgb2, Cib2, Isg15, Nrg1, Igf1, Esm1, Sfrp2, Gpnmb, Tnn, Lgals12, Itgb8, Itgb7 | 3.17 | 0.053 | molecular function |
| GO:0005102 | receptor binding | 28 | 3.47 Ć 10ā4 | Calca, Lpl, Clu, Kng1, Syne1, Slc9a3r2, Pdpn, Itgb8, Nefm, Apoe, Cd276, Fgb, Fga, Gstm2, Srms, Syk, Angpt1, Rarres2, Fibcd1, Cask, Nrg1, Rt1-Ce5, Nrg2, F2, Gnao1, Tg, Artn, Angptl4 | 2.12 | 0.0807 | molecular function |
| GO:0052689 | carboxylic ester hydrolase activity | 9 | 4.90 Ć 10ā4 | Ces2c, Acot8, Nceh1, Acot7, Ces1d, Ces1e, Acot2, Acot12, Baat | 4.82 | 0.080 | molecular function |
Top 5 enriched GO terms in each catalog associated with differentially expressed genes in Bailing capsule-treated group vs. DN group (p < 0.001).
FIGURE 6
Bioinformation analysis of enriched GO (gene ontology) term and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways of differentially expressed genes from BC-H (high dosage of Bailing Capsule) group vs. DN (diabetic nephropathy) group. (A) Top 5 significant GO term of DEGs (differentially expressed genes) identified by transcriptional analysis in each catalogue. (B) KEGG enrichment analysis for DEGs identified by transcriptomics analysis.
TABLE 3
| Pathway ID | Pathway term | Count | p value | Involved genes | Fold enrichment | FDR |
|---|---|---|---|---|---|---|
| rno04610 | Complement and coagulation cascades | 15 | 4.62 Ć 10ā6 | Fgb, Fga, Serpina1, Vwf, Itgb2, F2, Tfpi, F3, Clu, Kng1, C4a, Cfhr2, Klkb1, Mbl1, Loc100361907 | 4.5 | 6.84 Ć 10ā4 |
| rno03320 | PPAR signaling pathway | 15 | 4.62 Ć 10ā6 | Acsl1, Rgd1565355, Apoc3, Lpl, Cyp4a1, Acaa1b, Fabp5, Scd, Acox1, Plin2, Cd36, Hmgcs2, Angptl4, Ppara, Slc27a2 | 4.5 | 6.84 Ć 10ā4 |
| rno04979 | Cholesterol metabolism | 11 | 1.40 Ć 10ā5 | Mylip, Nceh1, Rgd1565355, Apoc2, Lpl, Apoc3, Apoa4, Apoe, Cd36, Angptl4, Apob | 5.79 | 0.0013 |
| rno01100 | Metabolic pathways | 89 | 1.87 Ć 10ā4 | Dgkg, Alas2, Gda, Mogat2, Pik3c2g, Cyp4a1, Hnmt, Atp12a, Eno4, Ftcd, Tm7sf2, Gls, Hyal1, Car15, Lrat, Hmgcs2, Enpp6, Lbr, Acot8, Cox8b, Entpd2, Acsl1, Glce, Lancl1, Amt, Acot12, Pipox, Ugt2b17, Mat1a, Baat, Sirt3, Hadha, Acox1, Itpka, Cyp2a1, Rgn, Acot2, Loc688321, Blvrb, Hprt1, L2hgdh, Atp6v0d2, St6galnac4, Dctpp1, Fbp2, Acot4, Blvra, Prps2, Sgms1, Pde1a, Mgst1, Ak2, Pdxp, Acaa1b, Fhit, Papss2, Hsd11b2, Neu3, Cyp3a9, Aanat, Mat2a, Cyp2c11, Oxct1, Aldh3b1, Hsd17b2, Rgd1559459, Hyi, Gpat2, B4galnt3, Atp6v1g3, Atp6v1c2, Cyp2j4, Gstm2, Gstm1, Mthfr, Acadsb, Sgpp2, Atp4b, Bhmt, Cyp24a1, Vnn1, Ggcx, Gsta3, Scd, Fasn, Cyp1a1, Pemt, Dhcr7, Galk1 | 1.43 | 0.0138 |
| rno04976 | Bile secretion | 13 | 2.52 Ć 10ā4 | Abcc2, Slc10a1, Ugt2b17, Slco1a6, Baat, Slc4a5, Nceh1, Slco1a1, Rgd1559459, Kcnn2, Abcb1a, Abcb1b, Cftr | 3.56 | 0.014 |
| rno01040 | Biosynthesis of unsaturated fatty acids | 7 | 0.0017 | Acot7, Scd, Acox1, Acot2, Acaa1b, Baat, Acot4 | 5.31 | 0.084 |
| rno04146 | Peroxisome | 10 | 0.0062 | Acot8, Pex19, Acsl1, Acox1, Prdx1, Pipox, Hacl1, Acaa1b, Baat, Slc27a2 | 2.96 | 0.259 |
| rno04151 | PI3K-Akt signaling pathway | 24 | 0.0075 | Ntrk1, Cdkn1a, Syk, Angpt1, Vwf, Itga3, Pdgfa, Igf1, Efna5, Areg, Efna1, Kitlg, Creb3l3, Ccne1, Tnn, Erbb4, Col4a1, Col9a1, Itgb8, Itgb7, Jak3, Fgfr2, Fgf21, Epha2 | 1.79 | 0.259 |
| rno05204 | Chemical carcinogenesis - DNA adducts | 9 | 0.0087 | Cyp3a9, Gstm2, Gstm1, Cyp2c11, Gsta3, Cyp1a1, Mgst1, Rgd1559459, Ugt2b17 | 3.05 | 0.259 |
| rno00980 | Metabolism of xenobiotics by cytochrome P450 | 9 | 0.0087 | Gstm2, Gstm1, Gsta3, Aldh3b1, Cyp1a1, Loc688321, Mgst1, Rgd1559459, Ugt2b17 | 3.05 | 0.259 |
Top 10 enriched KEGG pathways associated with differentially expressed genes in Bailing Capsule-treated group vs. DN group (p < 0.01).
FIGURE 7
PPAR (peroxisome proliferator activated receptor) signaling pathway affected by Bailing Capsule in DN (diabetic nephropathy) rats. All the different expression genes of BC-H (high dosage of Bailing Capsule) group vs. DN group in PPAR pathway. Red represents up-regulated; green represents down-regulated, white represents no significant change. FC, fold change in BC-H vs. DN group.
Real-time PCR validation of transcriptome analysis
Because the PPAR signaling pathway was among the top-ranked KEGG pathways associated with differentially expressed genes in the BC-H group vs. the DN group, we selected Acox1, Ppara, Scd, and Fasn genes in the PPAR pathway to compare their expression among the four different groups using real-time PCR. As shown in Figure 8, the mRNA levels of typical PPAR pathway genes, such as Acox1, Ppara, Scd, and Fasn, were determined by RT-qPCR. In DN rats, the mRNA levels of Acox1, Ppara, and Scd were significantly reduced while Fasn was obviously elevated compared to the control group, and Bailing Capsule treatment moderated the DN-induced changes in their transcript levels (p < 0.01, Figure 8). The results of the real-time PCR experiments for these genes were consistent with the transcriptome profiling.
FIGURE 8
Real-time PCR analysis for main DEGs (differentially expressed genes) expression. Bailing Capsule increased the mRNA levels of Acox1 (acyl-CoA oxidase 1), Ppara (peroxisome proliferator activated receptor Ī±), Scd (stearoyl-CoA desaturase) and inhibited the mRNA levels of Fasn (fatty acid synthase) in DN (diabetic nephropathy) rats. DN, diabetic nephropathy; BC-L, low dosage of Bailing Capsule; BC-H, high dosage of Bailing Capsule. Data are expressed as mean Ā± SD (n = 8). **p < 0.01.
Effects of bailing capsule on the peroxisome proliferator activated receptor pathway in diabetic nephropathy rats determined by immunohistochemistry staining
Similarly, it was apparent that more FASN-positive cells were present in the DN group, whereas fewer ACOX1-, PPARĪ±-, SCD-positive cells were present in the DN group relative to the control group (p < 0.01, Figure 9). Compared with the DN group, Bailing Capsule group showed a significant reduction in FASN and an increase in ACOX1, PPARĪ±, and SCD (p < 0.01, Figure 9).
FIGURE 9
Effect of Bailing Capsule treatment on the expression levels of PPAR (peroxisome proliferator activated receptor) pathway in kidneys of rats shown by immunohistochemistry staining. (A) Representative immunohistochemistry staining for ACOX1, PPARĪ±, SCD, and FASN (400 X). Scale bar = 20Ā Ī¼m. Black arrows showed ACOX1 (+), PPARĪ± (+), SCD (+) and FASN (+) cells in the kidney. Immunohistochemistry staining showed that Bailing Capsule increased the protein levels of ACOX1, PPARĪ±, SCD and inhibited the protein levels of FASN in DN rats. (B) Intensity of ACOX1 (+), PPARĪ± (+), SCD (+) and FASN (+) cells in kidney of rats. DN, diabetic nephropathy; BC-L, low dosage of Bailing Capsule; BC-H, high dosage of Bailing Capsule. Data are expressed as mean Ā± SD (n = 8). **p < 0.01.
Effects of bailing capsule on peroxisome proliferator activated receptor pathway protein expression levels in diabetic nephropathy rats determined by western blot
Western blotting analysis showed a significant reduction in the expression levels of ACOX1, PPARĪ±, and SCD in the kidneys of DN rats compared to the control group (p < 0.01, Figure 10). The expression levels of ACOX1, PPARĪ±, and SCD were upregulated in DN kidneys by high dosage of Bailing Capsule treatment (p < 0.01, Figure 10). In addition, Bailing Capsule inhibited the overexpression of FASN in DN rats (p < 0.01, Figure 10). Therefore, these findings showed that Bailing Capsule administration has an ameliorative effect on the PPAR pathway protein expression changes induced in the kidney by DN.
FIGURE 10
Effect of Bailing Capsule treatment on the relative protein expression levels of ACOX1, PPARĪ±, SCD and FASN in kidneys of rats shown by western blot. (A) Representative immunoblots of ACOX1, PPARĪ±, SCD, and FASN in kidneys of rats among different groups. Western blot showed that Bailing Capsule increased the protein levels of ACOX1, PPARĪ±, SCD and inhibited the protein levels of FASN in DN rats. (B) Quantification of immunoblots of ACOX1, PPARĪ±, SCD, and FASN. DN, diabetic nephropathy; BC-L, low dosage of Bailing Capsule; BC-H, high dosage of Bailing Capsule. Data are expressed as mean Ā± SD (n = 8). *p < 0.05 and **p < 0.01.
Discussion
In this study, our data show that both low and high doses of Bailing Capsule relieves hyperglycemia. The active ingredient of Bailing Capsule is C. sinensis. Yang et al. found that Cordyceps Cidadae polysaccharides improve glucose tolerance in DN rats (Yang et al., 2020). Hirsutella Sinensis, an anamorph of C. sinensis, reduced the levels of FBG in db/db mice (Lu et al., 2019). However, Kan et al. reported that C. sinensis did not improve hyperglycemia in a KK/HIJ diabetic mouse model (Kan et al., 2012). The differences in these results may be due to the differences in the diabetic animal models. Regarding the renal protection effect, this study found that Bailing Capsule reduced kidney hypertrophy, serum creatinine, BUN, 24Ā h urinary protein, and urinary albumin, and moderated the renal pathology in diabetic rats. C. sinensis was previously found to improve renal function in rodent diabetic models, including STZ-induced diabetic SD rats (Yang et al., 2020), unilateral ureteral obstruction (UUO) rats (Du et al., 2015), and db/db mice (Lu et al., 2019). The mechanism involves the alleviation of the inflammatory response and modulation of gut microbiota dysbiosis via suppression of the toll-like receptor 4 (TLR4)/nuclear factor-ĪŗB (NF-ĪŗB) and TGF-Ī²1/Smad signaling pathways (Yang et al., 2020), inhibiting Bcl-2-associated athanogene 3 (BAG3) and Ī±-smooth muscle actin (Ī±-SMA) in tubular epithelium (Du et al., 2015), inhibiting the p38 and extracellular regulated kinase (ERK) signaling pathways (Dong et al., 2019), reducing oxidative stress (Zhao et al., 2018) and influencing the mitogen-activated protein kinase (MAPK) signaling pathways (Zhao et al., 2018).
Among the pathogenic mechanisms of DN, the disturbance of fatty acid metabolism is an important factor leading to diabetic kidney injury (Herman-Edelstein et al., 2014; Chen et al., 2019). Prolonged hyperglycemia in diabetes aggravates fatty acid synthesis and TG accumulation. In a DN animal model, the effect of renal lipid accumulation was proven to be associated with glomerulosclerosis and tubulointerstitial injury (Jiang et al., 2007; Herman-Edelstein et al., 2014). In this study, we found Bailing Capsule reduced renal lipid deposition in DN rats, which may protect the kidney function.
Renal lipid accumulation in DN kidneys is mainly caused by the activation of fatty acid synthesis or the reduction of fatty acid oxidation. Some key enzymes catalyze fatty acid synthesis, such as FASN and acetyl-CoA carboxylase (ACC). In this study, we found marked upregulation of FASN mRNA and protein expression in DN kidney tissue. Previous studies also showed that the expression of FASN, which promotes fatty acid synthesis, is significantly increased in diabetic kidneys (Wang et al., 2005). High dosage of Bailing Capsule significantly inhibited lipid accumulation and reduced the expression levels of FASN in diabetic kidneys. SCD is a rate-limiting enzyme that converts saturated fatty acids to monounsaturated fatty acids, leading to the formation of neutral lipid droplets. In this study, the reduction in renal SCD expression by DN was alleviated by Bailing Capsule treatment. Iwai et al. found that in HFD-induced diabetic mice, renal SCD expression was significantly reduced. Overexpression of SCD reduced saturated fatty acid (SFA)-induced proximal tubular epithelial cell apoptosis. Enhancing SCD may become a promising therapeutic target to reduce SFA-induced lipotoxicity (Iwai et al., 2016). Thus, our data suggest that Bailing Capsule inhibits lipogenesis by activating SCD and suppressing FASN expression to moderate the progression of DN in rats.
In addition, in this study, we found differentially expressed kidney genes in Bailing Capsule-treated group were enriched in the PPAR pathway. In particular, Bailing Capsule treatment enhanced renal PPARĪ± expression. PPARĪ± is a ligand-dependent nuclear receptor that regulates lipid metabolism (Desvergne and Wahli, 1999). PPARĪ± can be activated by exogenous compounds, such as fatty acids (Desvergne and Wahli, 1999). PPARĪ± is mainly expressed in liver, kidney, and heart tissue, which have abundant mitochondrial and high activities of Ī²-oxidation. The main function of PPARĪ± is the regulation of fatty acid oxidation. In recent decades, several clinical trials demonstrated the beneficial effects of fenofibrate, a PPARĪ± agonist, on type 2 diabetic patients, including decreases in creatinine clearance and estimated glomerular filtration rate (eGFR) (Forsblom et al., 2010; Ismail-Beigi et al., 2010; Wright and Dodson, 2011). These promising results suggest that PPARĪ± is a potential drug target for DN treatment. Activation of PPARĪ± was found to ameliorate kidney function in diabetic animal models through multiple effects, including reduction of renal lipotoxicity (Tanaka et al., 2011), inflammation (Shin et al., 2009), fibrosis (Hou et al., 2010) and oxidative stress (Balakumar et al., 2009). On the other hand, diabetic PPARĪ± knockout mice have more severe kidney dysfunction symptoms, including albuminuria, glomerular sclerosis, and mesangial area expansion (Park et al., 2006). Our results showed that the expression of PPARĪ± was reduced in DN rat kidneys. And Bailing capsule treatment increased renal PPARĪ± expression. In addition, we found that the expression levels of PPARĪ± target genes, such as ACOX1 were also upregulated by Bailing Capsule treatment in DN rats. ACOX1 is a lipolytic enzyme that is regulated by PPARĪ± in lipid oxidation (Morise et al., 2009; Huang et al., 2012). In this study, we found that Bailing Capsule upregulated the expression of PPARĪ± and ACOX1 in diabetic kidneys. These results indicated that Bailing Capsule activates PPARĪ± pathway and downstream ACOX1 to enhance lipolysis, leading to the protection of kidney function.
The current work investigated the mechanism of Bailing Capsule moderates kidney function in DN rats. The transcriptome array is reliable and novel for searching the multiple targets and pathways affected by TCMs which have multiple active components. Furthermore, our study addressed the evidence that Bailing Capsule moderated renal lipid metabolism to moderated kidney function. However, this study has several limitations. Considering that HFD/STZ injection rats is one of DN model, other kidney disease models should be used to validate our findings. Moreover, the beginning time of treatment and duration also need to be explored.
Conclusion
In conclusion, this study provides convincing evidence of the efficacy of Bailing Capsule in preventing the development and progression of DN. The effects of Bailing Capsule are apparently associated with the activation of renal lipolysis and the inhibition of renal lipogenesis, leading to a reduction in kidney lipid accumulation (Figure 11).
FIGURE 11
Mechanism underlying renal-protective effect of Bailing Capsule. Intracellular lipid accumulation is governed by the syntheses and oxidation of triglyceride in kidney. In diabetic condition, lipogenesis is activated by FASN, whereas lipolysis is inhibited, through PPARĪ± and ACOX1, leading to renal lipid accumulation. Bailing capsule inhibited this triglyceride accumulation through PPAR pathway.
Statements
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/geo/query/acc.cgi?acc=GSE198021.
Ethics statement
This study was conducted in strict accordance with the recommendations and with the approval of the Animal Care Committee of the Peking Union medical Hospital Animal Ethics Committee (Project XHDW-2015-0051, 15 February 2015), and all efforts were made to minimize suffering.
Author contributions
XX designed the experiments, contributed reagents and materials; QZ conducted the experiments; MY, ML, and FP analyzed the data; QZ wrote the manuscript.
Funding
This work was supported by the grants from National Natural Science Foundation of China (No. 81870545, 82170854, 81870579, 81570715, 81170736), Beijing Natural Science Foundation (7202163), Beijing Municipal Science & Technology Commission (Z201100005520011), National Key Research and Development Program of China (2021YFC2501700, 2016YFA0101002, 2018YFC2001100), Scientific Activities Foundation for Selected Returned Overseas Professionals of Human Resources and Social Security Ministry, Beijing Dongcheng District Outstanding Talent Funding Project (2019DCT-M-05), CAMS Innovation Fund for Medical Sciences (CIFMS2021-1-I2M-002, CIFMS2017-I2M-1-008).
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.
References
1
AiharaK.IkedaY.YagiS.AkaikeM.MatsumotoT. (2010). Transforming growth factor-Ī²1 as a common target molecule for development of cardiovascular diseases, renal insufficiency and metabolic syndrome. Cardiol. Res. Pract.2011, 175381. 10.4061/2011/175381PubMed Abstract | CrossRef Full Text | Google Scholar
2
BalakumarP.ChakkarwarV. A.SinghM. (2009). Ameliorative effect of combination of benfotiamine and fenofibrate in diabetes-induced vascular endothelial dysfunction and nephropathy in the rat. Mol. Cell. Biochem.320, 149ā162. 10.1007/s11010-008-9917-zPubMed Abstract | CrossRef Full Text | Google Scholar
3
BaynesJ. W. (1991). Role of oxidative stress in development of complications in diabetes. Diabetes40, 405ā412. 10.2337/diab.40.4.405PubMed Abstract | CrossRef Full Text | Google Scholar
4
ChenX.HanY.GaoP.YangM.XiaoL.XiongX.et al (2019). Disulfide-bond A oxidoreductase-like protein protects against ectopic fat deposition and lipid-related kidney damage in diabetic nephropathy. Kidney Int.95, 880ā895. 10.1016/j.kint.2018.10.038PubMed Abstract | CrossRef Full Text | Google Scholar
5
DaiB.WangZ. Z.ZhangH.HanM. X.ZhangG. X.ChenJ. W. (2020). Antihypertensive properties of a traditional Chinese medicine GAO-ZI-YAO in elderly spontaneous hypertensive rats. Biomed. Pharmacother.131, 110739. 10.1016/j.biopha.2020.110739PubMed Abstract | CrossRef Full Text | Google Scholar
6
DandaR. S.HabibaN. M.Rincon-CholesH.BhandariB. K.BarnesJ. L.AbboudH. E.et al (2005). Kidney involvement in a nongenetic rat model of type 2 diabetes. Kidney Int.68, 2562ā2571. 10.1111/j.1523-1755.2005.00727.xPubMed Abstract | CrossRef Full Text | Google Scholar
7
DengY. Y.ChenY. P.HeX. L.LiL. (2001). Study of Cordyceps on mechanism in delaying chronic renal failure. Chin. J. Integr. Traditional West. Nephrol.2, 318ā323. 10.3969/j.issn.1009-587X.2001.07.004CrossRef Full Text | Google Scholar
8
DesvergneB.WahliW. (1999). Peroxisome proliferator-activated receptors: Nuclear control of metabolism. Endocr. Rev.20, 649ā688. 10.1210/edrv.20.5.0380PubMed Abstract | CrossRef Full Text | Google Scholar
9
DingM.TangZ.LiuW.ShaoT.YuanP.ChenK.et al (2021). Burdock fructooligosaccharide attenuates high glucose-induced apoptosis and oxidative stress injury in renal tubular epithelial cells. Front. Pharmacol.12, 784187. 10.3389/fphar.2021.784187PubMed Abstract | CrossRef Full Text | Google Scholar
10
DongZ.SunY.WeiG.LiS.ZhaoZ. (2019). A nucleoside/nucleobase-rich extract from Cordyceps sinensis inhibits the epithelial-mesenchymal transition and protects against renal fibrosis in diabetic nephropathy. Molecules24, E4119. 10.3390/molecules24224119PubMed Abstract | CrossRef Full Text | Google Scholar
11
DuF.LiS.WangT.ZhangH. Y.ZongZ. H.DuZ. X.et al (2015). Cordyceps sinensis attenuates renal fibrosis and suppresses BAG3 induction in obstructed rat kidney. Am. J. Transl. Res.7, 932ā940. PubMed Abstract | Google Scholar
12
FolchJ.LeesM.Sloane StanleyG. H. (1957). A simple method for the isolation and purification of total lipides from animal tissues. J. Biol. Chem.226, 497ā509. 10.1016/s0021-9258(18)64849-5PubMed Abstract | CrossRef Full Text | Google Scholar
13
ForsblomC.HiukkaA.LeinonenE. S.SundvallJ.GroopP. H.TaskinenM. R. (2010). Effects of long-term fenofibrate treatment on markers of renal function in type 2 diabetes: The FIELD helsinki substudy. Diabetes Care33, 215ā220. 10.2337/dc09-0621PubMed Abstract | CrossRef Full Text | Google Scholar
14
FurmanB. L. (2015). Streptozotocin-induced diabetic models in mice and rats. Curr. Protoc. Pharmacol.70, 5.47.1ā5.47.20. 10.1002/0471141755.ph0547s70CrossRef Full Text | Google Scholar
15
GrossC.BelvilleC.LavergneM.CholtusH.JabaudonM.BlondonnetR.et al (2020). Advanced glycation end products and receptor (RAGE) promote wound healing of human corneal epithelial cells. Invest. Ophthalmol. Vis. Sci.61, 14. 10.1167/iovs.61.3.14PubMed Abstract | CrossRef Full Text | Google Scholar
16
GuoZ.JiaJ.TuY.JinC.GuoC.SongF.et al (2021). Altered jagged1-notch1 signaling in enhanced dysfunctional neovascularization and delayed angiogenesis after ischemic stroke in HFD/STZ induced type 2 diabetes rats. Front. Physiol.12, 687947. 10.3389/fphys.2021.687947PubMed Abstract | CrossRef Full Text | Google Scholar
17
Herman-EdelsteinM.ScherzerP.TobarA.LeviM.GafterU. (2014). Altered renal lipid metabolism and renal lipid accumulation in human diabetic nephropathy. J. Lipid Res.55, 561ā572. 10.1194/jlr.P040501PubMed Abstract | CrossRef Full Text | Google Scholar
18
HouX.ShenY. H.LiC.WangF.ZhangC.BuP.et al (2010). PPARalpha agonist fenofibrate protects the kidney from hypertensive injury in spontaneously hypertensive rats via inhibition of oxidative stress and MAPK activity. Biochem. Biophys. Res. Commun.394, 653ā659. 10.1016/j.bbrc.2010.03.043PubMed Abstract | CrossRef Full Text | Google Scholar
19
HuX.WangJ.YangH.JiS.LiY.XuB.et al (2021). Bailing Capsule combined with Ī±-ketoacid tablets for stage 3 chronic kidney disease: Protocol of a double-blinded, randomized, controlled trial. Med. Baltim.100, e25759. 10.1097/MD.0000000000025759CrossRef Full Text | Google Scholar
20
Huang DaW.ShermanB. T.LempickiR. A. (2009). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc.4, 44ā57. 10.1038/nprot.2008.211PubMed Abstract | CrossRef Full Text | Google Scholar
21
HuangJ.JiaY.FuT.ViswakarmaN.BaiL.RaoM. S.et al (2012). Sustained activation of PPARĪ± by endogenous ligands increases hepatic fatty acid oxidation and prevents obesity in ob/ob mice. Faseb J.26, 628ā638. 10.1096/fj.11-194019PubMed Abstract | CrossRef Full Text | Google Scholar
22
Ismail-BeigiF.CravenT.BanerjiM. A.BasileJ.CallesJ.CohenR. M.et al (2010). Effect of intensive treatment of hyperglycaemia on microvascular outcomes in type 2 diabetes: An analysis of the ACCORD randomised trial. Lancet376, 419ā430. 10.1016/S0140-6736(10)60576-4PubMed Abstract | CrossRef Full Text | Google Scholar
23
IwaiT.KumeS.Chin-KanasakiM.KuwagataS.ArakiH.TakedaN.et al (2016). Stearoyl-CoA desaturase-1 protects cells against lipotoxicity-mediated apoptosis in proximal tubular cells. Int. J. Mol. Sci.17, E1868. 10.3390/ijms17111868PubMed Abstract | CrossRef Full Text | Google Scholar
24
JiangT.WangX. X.ScherzerP.WilsonP.TallmanJ.TakahashiH.et al (2007). Farnesoid X receptor modulates renal lipid metabolism, fibrosis, and diabetic nephropathy. Diabetes56, 2485ā2493. 10.2337/db06-1642PubMed Abstract | CrossRef Full Text | Google Scholar
25
KanW. C.WangH. Y.ChienC. C.LiS. L.ChenY. C.ChangL. H.et al (2012). Effects of extract from solid-state fermented Cordyceps sinensis on type 2 diabetes mellitus. Evid. Based. Complement. Altern. Med.2012, 743107. 10.1155/2012/743107PubMed Abstract | CrossRef Full Text | Google Scholar
26
LapollaA.TraldiP.FedeleD. (2005). Importance of measuring products of non-enzymatic glycation of proteins. Clin. Biochem.38, 103ā115. 10.1016/j.clinbiochem.2004.09.007PubMed Abstract | CrossRef Full Text | Google Scholar
27
LiuY.HuangH.GaoR.LiuY. (2020). Dynamic phenotypes and molecular mechanisms to understand the pathogenesis of diabetic nephropathy in two widely used animal models of type 2 diabetes mellitus. Front. Cell Dev. Biol.8, 172. 10.3389/fcell.2020.00172PubMed Abstract | CrossRef Full Text | Google Scholar
28
LoganS.JiangC.YanY.InagakiY.ArzuaT.BaiX. (2018). Propofol alters long non-coding RNA profiles in the neonatal mouse Hippocampus: Implication of novel mechanisms in anesthetic-induced developmental neurotoxicity. Cell. Physiol. biochem.49, 2496ā2510. 10.1159/000493875PubMed Abstract | CrossRef Full Text | Google Scholar
29
LuZ.LiS.SunR.JiaX.XuC.AaJ.et al (2019). Hirsutella sinensis treatment shows protective effects on renal injury and metabolic modulation in db/db mice. Evid. Based. Complement. Altern. Med.2019, 4732858. 10.1155/2019/4732858PubMed Abstract | CrossRef Full Text | Google Scholar
30
LucK.Schramm-LucA.GuzikT. J.MikolajczykT. P. (2019). Oxidative stress and inflammatory markers in prediabetes and diabetes. J. Physiol. Pharmacol.70. 10.26402/jpp.2019.6.01PubMed Abstract | CrossRef Full Text | Google Scholar
31
MaX.YanW. (2021). Effects of bailing capsule combined with valsartan on renal protection and cellular immune function in patients with chronic glomerulonephritis. Chin. Archives Traditional Chin. Med.39, 220ā223. 10.13193/j.issn.1673-7717.2021.09.055CrossRef Full Text | Google Scholar
32
MehlemA.HagbergC. E.MuhlL.ErikssonU.FalkevallA. (2013). Imaging of neutral lipids by oil red O for analyzing the metabolic status in health and disease. Nat. Protoc.8, 1149ā1154. 10.1038/nprot.2013.055PubMed Abstract | CrossRef Full Text | Google Scholar
33
MoriseA.ThomasC.LandrierJ. F.BesnardP.HermierD. (2009). Hepatic lipid metabolism response to dietary fatty acids is differently modulated by PPARalpha in male and female mice. Eur. J. Nutr.48, 465ā473. 10.1007/s00394-009-0037-7PubMed Abstract | CrossRef Full Text | Google Scholar
34
NathS.GhoshS. K.ChoudhuryY. (2017). A murine model of type 2 diabetes mellitus developed using a combination of high fat diet and multiple low doses of streptozotocin treatment mimics the metabolic characteristics of type 2 diabetes mellitus in humans. J. Pharmacol. Toxicol. Methods84, 20ā30. 10.1016/j.vascn.2016.10.007PubMed Abstract | CrossRef Full Text | Google Scholar
35
NiuC.LiH. (2019). Effect of Bailing capsule combined with simvastatin on renal function and oxidative stress in patients with diabeitc nephropathy. Foreine Mecical Sci. Sction Medgeogr.40, 275ā276. 10.3969/j.issn.1001-8883.2019.03.018CrossRef Full Text | Google Scholar
36
ParkC. W.KimH. W.KoS. H.ChungH. W.LimS. W.YangC. W.et al (2006). Accelerated diabetic nephropathy in mice lacking the peroxisome proliferator-activated receptor alpha. Diabetes55, 885ā893. 10.2337/diabetes.55.04.06.db05-1329PubMed Abstract | CrossRef Full Text | Google Scholar
37
RenH. J.SunY. L.YuanB. (2019). Chinese patent medicine bailing capsule for treating lupus nephritis: A protocol for systematic review and meta-analysis. Med. Baltim.98, e17041. 10.1097/MD.0000000000017041PubMed Abstract | CrossRef Full Text | Google Scholar
38
SaidE.ZaitoneS. A.EldosokyM.ElsherbinyN. M. (2018). Nifuroxazide, a STAT3 inhibitor, mitigates inflammatory burden and protects against diabetes-induced nephropathy in rats. Chem. Biol. Interact.281, 111ā120. 10.1016/j.cbi.2017.12.030PubMed Abstract | CrossRef Full Text | Google Scholar
39
ShengX.DongY.ChengD.WangN.GuoY. (2020). Efficacy and safety of bailing capsules in the treatment of type 2 diabetic nephropathy: A meta-analysis. Ann. Palliat. Med.9, 3885ā3898. 10.21037/apm-20-1799PubMed Abstract | CrossRef Full Text | Google Scholar
40
ShinS. J.LimJ. H.ChungS.YounD. Y.ChungH. W.KimH. W.et al (2009). Peroxisome proliferator-activated receptor-alpha activator fenofibrate prevents high-fat diet-induced renal lipotoxicity in spontaneously hypertensive rats. Hypertens. Res.32, 835ā845. 10.1038/hr.2009.107PubMed Abstract | CrossRef Full Text | Google Scholar
41
SunH.SaeediP.KarurangaS.PinkepankM.OgurtsovaK.DuncanB. B.et al (2022). IDF Diabetes Atlas: Global, regional and country-level diabetes prevalence estimates for 2021 and projections for 2045. Diabetes Res. Clin. Pract.183, 109119. 10.1016/j.diabres.2021.109119PubMed Abstract | CrossRef Full Text | Google Scholar
42
TanakaY.KumeS.ArakiS.IsshikiK.Chin-KanasakiM.SakaguchiM.et al (2011). Fenofibrate, a PPARĪ± agonist, has renoprotective effects in mice by enhancing renal lipolysis. Kidney Int.79, 871ā882. 10.1038/ki.2010.530PubMed Abstract | CrossRef Full Text | Google Scholar
43
WangD.ZhangG.ChenX.WeiT.LiuC.ChenC.et al (2018). Sitagliptin ameliorates diabetic nephropathy by blocking TGF-Ī²1/Smad signaling pathway. Int. J. Mol. Med.41, 2784ā2792. 10.3892/ijmm.2018.3504PubMed Abstract | CrossRef Full Text | Google Scholar
44
WangW.ZhangX. N.YinH.LiX. B.HuX. P.LiuH.et al (2013). Effects of bailing capsules for renal transplant recipients: A retrospective clinical study. Chin. Med. J.126, 1895ā1899. 10.3760/cma.j.issn.0366-6999.20130483PubMed Abstract | CrossRef Full Text | Google Scholar
45
WangZ.JiangT.LiJ.ProctorG.McmanamanJ. L.LuciaS.et al (2005). Regulation of renal lipid metabolism, lipid accumulation, and glomerulosclerosis in FVBdb/db mice with type 2 diabetes. Diabetes54, 2328ā2335. 10.2337/diabetes.54.8.2328PubMed Abstract | CrossRef Full Text | Google Scholar
46
WrightA. D.DodsonP. M. (2011). Medical management of diabetic retinopathy: Fenofibrate and ACCORD eye studies. Eye (Lond)25, 843ā849. 10.1038/eye.2011.62PubMed Abstract | CrossRef Full Text | Google Scholar
47
WuJ.WeiZ.ChengP.QianC.XuF.YangY.et al (2020). Rhein modulates host purine metabolism in intestine through gut microbiota and ameliorates experimental colitis. Theranostics10, 10665ā10679. 10.7150/thno.43528PubMed Abstract | CrossRef Full Text | Google Scholar
48
XiaoC.XiaoP.LiX.LiX.LiH.ChenY.et al (2018). Cordyceps sinensis may inhibit Th22 cell chemotaxis to improve kidney function in lgA nephropathy. Am. J. Transl. Res.10, 857ā865. PubMed Abstract | Google Scholar
49
XiaoL.XuX.ZhangF.WangM.XuY.TangD.et al (2017). The mitochondria-targeted antioxidant MitoQ ameliorated tubular injury mediated by mitophagy in diabetic kidney disease via Nrf2/PINK1. Redox Biol.11, 297ā311. 10.1016/j.redox.2016.12.022PubMed Abstract | CrossRef Full Text | Google Scholar
50
XuH.LiX.YuanX.YuanQ.ChenW.PengZ.et al (2020). A meta-analysis of the clinical efficacy and safety of Bailing capsules in the treatment of nephrotic syndrome. Ann. Palliat. Med.9, 3170ā3181. 10.21037/apm-20-1252PubMed Abstract | CrossRef Full Text | Google Scholar
51
XuY.BaiL.ChenX.LiY.QinY.MengX.et al (2018). 6-Shogaol ameliorates diabetic nephropathy through anti-inflammatory, hyperlipidemic, anti-oxidative activity in db/db mice. Biomed. Pharmacother.97, 633ā641. 10.1016/j.biopha.2017.10.084PubMed Abstract | CrossRef Full Text | Google Scholar
52
YangJ.DongH.WangY.JiangY.ZhangW.LuY.et al (2020). Cordyceps cicadae polysaccharides ameliorated renal interstitial fibrosis in diabetic nephropathy rats by repressing inflammation and modulating gut microbiota dysbiosis. Int. J. Biol. Macromol.163, 442ā456. 10.1016/j.ijbiomac.2020.06.153PubMed Abstract | CrossRef Full Text | Google Scholar
53
YuW.DuanS.YuZ. (2021). The effect of Bailing capsules combined with losartan to treat diabetic glomerulosclerosis and the combination's effect on blood and urine biochemistry. Am. J. Transl. Res.13, 6873ā6880. PubMed Abstract | Google Scholar
54
ZhangQ.XiaoX.ZhengJ.LiM.YuM.PingF.et al (2019). Shenqi jiangtang granule ameliorates kidney function by inhibiting apoptosis in a diabetic rat model. Evid. Based. Complement. Altern. Med.2019, 3240618. 10.1155/2019/3240618PubMed Abstract | CrossRef Full Text | Google Scholar
55
ZhaoK.GaoQ.ZongC.GeL.LiuJ. (2018). Cordyceps sinensis prevents contrast-induced nephropathy in diabetic rats: Its underlying mechanism. Int. J. Clin. Exp. Pathol.11, 5571ā5580. PubMed Abstract | Google Scholar
56
ZimmetP.AlbertiK. G.MaglianoD. J.BennettP. H. (2016). Diabetes mellitus statistics on prevalence and mortality: Facts and fallacies. Nat. Rev. Endocrinol.12, 616ā622. 10.1038/nrendo.2016.105PubMed Abstract | CrossRef Full Text | Google Scholar
Summary
Keywords
Cordyceps sinesis, diabetic nephropathy, lipid accumulation, transcriptome analysis, PPARĪ±
Citation
Zhang Q, Xiao X, Li M, Yu M and Ping F (2022) Bailing capsule (Cordyceps sinensis) ameliorates renal triglyceride accumulation through the PPARĪ± pathway in diabetic rats. Front. Pharmacol. 13:915592. doi: 10.3389/fphar.2022.915592
Received
08 April 2022
Accepted
27 June 2022
Published
25 August 2022
Volume
13 - 2022
Edited by
Vincenzo De Feo, University of Salerno, Italy
Reviewed by
Vikram Patial, Institute of Himalayan Bioresource Technology (CSIR), India
Soheb Anwar Mohammed, University of Pittsburgh, United States
Updates
Copyright
Ā© 2022 Zhang, Xiao, Li, Yu and Ping.
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: Xinhua Xiao, xiaoxh2014@vip.163.com
This article was submitted to Ethnopharmacology, a section of the journal Frontiers in Pharmacology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.