Label-free Quantitative Proteomics for Investigating the Therapeutic Mechanism of Sijunzi Decoction on Spleen Deficiency Syndrome

Background: Sijunzi decoction is a traditional Chinese medicine prescription which is used to treat spleen deficiency syndrome. But the therapeutic mechanism of Sijunzi decoction on spleen deficiency syndrome is unclear. Methods: A rat model of reserpine-induced spleen deficiency syndrome was established. The therapeutic effects of Sijunzi decoction were assessed by histopathological examination and enzyme activity analysis. A label-free proteomics method was used to assess the changes in the endogenous protein. The differentially expressed proteins (DEPs) were analyzed using bioinformatics methods. Western blot was used to validate the proteomic data. Results: Sijunzi decoction could significantly restore damaged liver tissue and enzymatic activities. DEPs of Sijunzi decoction-intervened spleen deficiency syndrome mainly regulated various metabolic processes involved in cellular responses to enzyme activity for catalysing, transforming and transferring. The four pathways (bile secretion, pyruvate metabolism, cysteine and methionine metabolism, and fatty acid metabolism) might be closely related with the intervention mechanism of Sijunzi decoction. Protein H2afx, interferon stimulated gene 15 protein, and signal transducers and activators of transcription1 proteins were predicted to play the important roles in Sijunzi decoction action networks. Conclusions : Sijunzi decoction could exert intervention action by restoring DEPs involved in potential target pathways. Current results wish to provide novel perspectives for revealing the potential protein targets of Sijunzi decoction on spleen deficiency syndrome and facilitating its clinical usage.

4 recent years, label-free proteomics methods have been considered as a promising alternative because of the reliability and validity. There are relatively few reports on the proteomics research of spleen deficiency syndrome. Li et al first reported the differentially expressed proteins and pathways of spleen-deficiency syndrome rats using isotope labelling proteomics methods [19]. Liu et al investigated the effect of traditional Chinese medicine Yi Tang Kang on metabolic syndrome of spleen deficiency with iTRAQ-based proteomics anslysis [20]. At present, there are almost no studies on spleen deficiency syndrome using label-free proteomics approaches.
It is unclear whether the pharmacological mechanisms of SJZD on spleen deficiency syndrome target multiple pathways. In this study, under the guidance of TCM theory of the spleen, and by reflecting on the concept of "wholeness", a characteristic subject of TCM, we integrated data from clinical symptoms monitoring, physiological and biochemical indices determinations, histopathological observations, and protein expression level quantification. To achieve the differentially expressed spectra of proteins using label-free proteomic technology was the primary focus of this study. Furthermore, the mechanism of action of Sijunzi decoction in spleen deficiency syndrome was further developed using bioinformatic analysis. This study takes a holistic approach to develop our understanding of the effect of Sijunzi decoction on spleen deficiency syndrome. Medicine, Xinxiang Medical University, (Xinxiang, China). The preparation and manufacturing process of all herbs followed standard operating procedures according to a previous study [21]. The raw materials were socked in distilled water and boiled twice for 40 min, respectively. The drug solution was filtered and concentrated in a water bath at 70 o C to 80 o C. The final concentration was 2 g of crude drug materials/mL, and the solution was stored at 4 o C. For the delivery of Sijunzi decoction using the drug gavage method, 1 mL/100 g of body weight was delivered. The dosage of Sijunzi decoction used was the clinical dose converted into an experimental equivalent dose based on the standard body shape coefficient [22].

Animals and treatments
Thirty male Sprague Dawley 3-month-old rats with a body mass of 220 ± 10 g were provided by the Laboratory Animal Center of Henan Province (Zhengzhou, China) (Rodent license No. SCXK 20100002). After the cages were separated, the animals were kept in sterilized rooms at a constant temperature (20 ± 2 o C) and humidity (50 ± 20%), as well as natural circadian rhythm illumination. The rats were observed for 3 days to ensure they adapted to the environment and were healthy.
These rats were then randomly divided into three groups: the normal control group (NC, n = 10), the disease model group (M, n = 10), and the Sijunzi decoction-treated group (SJZD, n = 10). The preparation of the animal model followed standard operating procedures according to a previous study [21].NC rats were given a subcutaneous injection of 0.2 mL normal saline daily at 9 am and an intragastric dose of 2 mL distilled water daily at 4 pm.M rats were given a daily subcutaneous injection of 0.6 mg/kg of reserpineat 9 am and an intragastric dose of 2 mL distilled water at 4 pm. Rats in the Sijunzi decoction-treated group were given a subcutaneous injection of 0.6 mg/kg of 6 reserpinedaily at 9 am and were given S Sijunzi decoction via gavage (1 mL/100 g of body weight) at 4 pm.The rats were treated daily for 14 days and weighed every morning.
Successful modeling of spleen deficiency syndrome was confirmed by diarrhea, loss of appetite, weight loss, wilting, dull coat, and collapsing [23].
Animal welfare and experimental procedures were performed in accordance with the Guide for the Care and Use of Laboratory Animals (Ministry of Science and Technology of China, 2006) and were approved by the Animal Ethics Committee of Xinxiang Medical University.

Sample collection
After the model was finished, the rats were anesthetized using 10% chloral hydrate and dissected. Liver tissues were resected and divided into three parts. One part of the samples was fixed in 3.7% (v/v) formalin solution for histopathological examination. The others were frozen in liquid nitrogen immediately and stored at -80 o C for further analysis (enzyme activity determination and proteomics analysis).

Histopathological examination
The same parts from resected liver tissues of each group were fixed in 3.7% (v/v) formalin solution immediately and imbedded into paraffin wax. Sections were then stained with hematoxylin-eosin stain and examined under a biological microscope to observe any morphological differences.

Enzyme activity analysis
Rat livers were comminuted using liquid nitrogen after removal from the refrigerator. The activities of SDH and LDH in the liver tissues were analyzed using the SDH and LDH assay kits (Nanjing Jiancheng Biological Engineering Co., Ltd, Nanjing, China), respectively.

Protein extraction and Protein digestion with the FASP method
The proteins from same parts of liver tissues were cut into pieces with clean scissors and lysed in an appropriate amount of SDT lysis buffer (SDS, dithiothreitol, Tris). Then, the liver lysate was ground repeatedly with a full-automatic sample frozen grinding machine (JXFSTPRP-III, Shanghai Jingxin Industrial Development CO., LTD, Shanghai, China) until completely pulverized. Samples were homogenized with a Vortex mixer (Digital Vortex-Genie 2; Scientific Industries, Inc., Bohemia, NY) and Ultrasonic treatments (100 W, 10s with 10s rest, 10 times in total) were utilized. After 5min of boiling in a water bath, the mixture was cooled to room temperature. After centrifugation, the supernatant of each sample was quantified using the fluorescence protein quantification method and analysed by SDS-PAGE.
A modified protocol based on the previous method was used for protein digestion [24].
Approximately 300 μg of protein extract from each sample was added into 200 μL of urea (UA) buffer (8 M Urea, 150 mM Tris-HCl pH 8.5). After mixing, the mixture was transferred into Microcon 10 kD ultra-fraction tube (Merck-Millipore, Germany) and then centrifuged at 14,000 g for 10min. The fluid in the bottom of the collecting pipe was discarded. 200 μL of UA buffer were put into an ultra-fraction tube, rotated at 14000 rpm (4 °C) for 15min. With the eluates discarded, 100 μL of iodoacetamide (IAA, Bio-Rad) buffer (50 mM IAA in UA buffer) were pipetted to the ultra-fraction tubes, the ultra-fraction tubes were shaken at 600 rpm for 1min, after which they were left at room temperature without light for 30 min, and then centrifuged at 14,000 × g for 10 min. To wipe off IAA, 100 μL UA buffer was added to the ultra-fraction tubes, and the spun in the above conditions; this was repeated two times. After that, the ultra-fraction tubes were washed twice with 100 μL of 50 mM NH 4 HCO 3 and centrifuged at 14,000 g for 10min at room temperature. Then, 40 μL of trypsin buffer (5 μL of trypsin diluted in 40μL of 50 mM NH 4 HCO 3 ; trypsin from Promega, Madison, WI,) was added into the ultra-fraction tubes. The samples were oscillated at 600 rpm for 1min. Then digestion was achieved at 37 °C for 16 to 18 h. The eluates were collected with new collecting tubes after centrifuging at 14,000 g for 10min. 25μL 25 mM 8 NH 4 HCO 3 was injected into each ultra-fraction tube. The ultra-fraction tubes were centrifuged at 14,000 g for 10min after being oscillated at 600 rpm for 1 min.This was repeated one times. Then, the filtrate was incorporated and lyophilized. The freezingdrying peptides were dissolved with 50 μL of a 0.1% trifluoroacetic acid (TFA, Sigma-Aldrich, USA). The fluorescence method (excitation wavelength of 295 nm, emission wavelength of 350 nm) was carried out to quantify the peptides and then analyzed by liquid chromatography-mass spectrometry/mass spectrometry (LC/MS/MS).

LC-MS/MS analysis
A 5 μg sample of each group was analyzed in triplicate using LC/MS/MS. Analysis was conducted as previously described with few modifications [25]. The separation was done using an EASY-nLC1000 nanoflow HPLC (Thermo Scientific, USA). The mass spectrometry proteomics data have been deposited to the ProteomeXchange 9 Consortium (http://proteomecentral. proteomexchange.org) via the iProX partner repository with the dataset identifier PXD015449 [26].

Identification of DEPs
DEPs were identified with a fold-change (FC) greater than 1.5 or less than 0.667. These selected DEPs were used in a random variance model t-test performed using SPSS version 20.0 (IBM Corp., Armonk, NY, USA). Only proteins with a p value < 0.05 were considered to be significantly differentially expressed as previously reported [28,29].

Bioinformatics analysis
A clustering analysis map was built using cluster 3.0 version 2.3 (Bio Fly Bioscience; www.bangfeibio.com/company) to identify DEPs with similar expression. Gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation for each protein were retrieved from the GO (http://www.geneontology.org/) and KEGG pathway database (http://www.genome.jp/Pathway) [30,31], respectively. Proteins in the GO analysis were related at three ontologies-biological process (BP), cellular components (CC), and molecular functions (MF) [32]. Protein-protein interaction networks of the DEPs were analyzed by ingenuity pathway analysis (IPA) version 9.0 (Redwood City, CA, http://www. ingenuity.com ).

Western blot analysis
Total proteins were extracted from samples of liver tissues (100 mg), and 40 μg proteins were separated by 12% SDS-PAGE and then transferred to PVDF membranes (Merck-Millipore, Germany). The membranes were blocked with 5% (w/v) non-fat milk in TBST for 1 h and then incubated with the primary antibodies overnight at 4 °C. After washing three times, the membranes were incubated with HRP-conjugated secondary antibody for 1 h at room temperature. Subsequently, the membranes were washed and were visualized with the enhanced chemiluminescence reagent (Sigma Chemical Co, USA) and the images were generated using an AzureC600 instrument (Azure Biosystems, USA) and were calculated by Gray analysis. All the primary and secondary antibodies used in this western blot analysis came from Cell Signaling Technology (Cell Signaling Technology, USA).

Statistical analysis
In the proteomics experiments, statistically significant differences (p value < 0.05) were analyzed between the M and SJZD or M and NC group. In the enzyme activity assays, data were reported as the mean ± standard error of mean (SEM) and statistically significant differences (p value < 0.05) were analyzed between the M and SJZD or M and NC group by one-way analysis of variance (ANOVA), followed by Tukey's multiple comparisons test. All statistical analysis was performed with SPSS 20.0 (SPSS, Armonk, New York, USA).

Body weight and viscera indices
Body weight and viscera indices results are shown in Table 1. Compared with the NC group, the body weight and viscera (spleen and thymus) index of rats in the M group were significantly reduced (p value < 0.05) by the 14th day. Meanwhile, these indicators in the SJZD group were significantly increased (p value < 0.05) compared to the M group. These results indicate that SJZD significantly reduces organ damage and increases the body weight and immunity of rats with spleen deficiency syndrome.

Histopathological examination of liver tissues
Liver histology was applied to assess the effect of SJZD on spleen deficiency syndrome

Enzyme activity analysis
The SDH and LDH assay kits (Nanjing Jiancheng Biological Engineering Co., Ltd, China) were used for enzyme activity analysis. The SDH and LDH activity in the liver of each group can be seen in Fig. 2. These results indicated that the activity of SDH was significantly decreased, while LDH activity was significantly increased (p value < 0.05) in the M group compared to the NC group. However, the activity of SDH was significantly increased, and LDH activity was significantly decreased (p value < .05) in the SJZD-treated group compared with the M group rats (p value < 0.05). This demonstrated that the activity of SDH was significantly elevated and LDH activity reduced after SJZD treatment.
Interestingly, the activities of SDH and LDH showed no significant difference between the SJZD group and the NC group (p value > 0.05).

Label-free proteomics analysis
In this study, proteomics techniques were used to analyze the effect of SJZD on spleen deficiency syndrome and a total of 20,759 peptides and 3,021 proteins were identified.  (Table S1-S4). As seen in Fig. 3, a Venn diagram analysis highlighted that there were 98 DEPs in common between the two comparisons.

Cluster analysis of DEPs
In order to further investigate the accuracy of the selected DEPs, clustering analysis was used to evaluate the DEPs. As seen in Fig.4, the DEPs identified in the M group and the SJZD-treated group can distinguish between the two groups. This indicates that the DEPs identified in this study may represent key differences between the two groups of samples.

GO analysis
The results of functional annotations of the DEPs using GO analysis are displayed in Fig. 5.
When the CC of these DEPs was analyzed, the majority of the DEPs were associated with the cytoplasm, cell membrane, and extracellular region (Fig.5A). In the MF analysis, the majority of the DEPs were associated with enzyme activity in catalyzing, transforming and 13 transferring (Fig. 5B). The DEPs analyzed for their BP were found to be enriched in various metabolic processes (Fig.5C).

KEGG pathway analysis
To identify the potential target pathway of SJZD in spleen deficiency syndrome, we analysed the target pathway of the DEPs in rats with spleen deficiency syndrome and the target pathway of the DEPs in SJZD-treated rats with spleen deficiency syndrome. Pathway enrichment analysis of the DEPs was shown in Table 2 and Table 3. As seen in Table 2, the DEPs in rats with spleen deficiency syndrome were enriched in 11 biological pathways, namely, cysteine and methionine metabolism, bile secretion, retinol metabolism, drug metabolism-cytochrome P450, fatty acid metabolism, peroxisome, pyruvate metabolism, butanoate metabolism, phenylalanine metabolism, terpenoid backbone biosynthesis, and the AMPK signaling pathway. As seen in Table 3, the DEPs in SJZD-treated rats with spleen deficiency syndrome were significantly enriched in five biological pathways. These were bile secretion, pyruvate metabolism, cysteine and methionine metabolism, complement and coagulation cascades, and fatty acid metabolism. By comparing these two results, it was found that four metabolic pathways were common to both groups and these were bile secretion, pyruvate metabolism, cysteine and methionine metabolism, and fatty acid metabolism. It can be speculated that the four metabolic pathways are potential target pathways of SJZD treatment spleen deficiency syndrome.

Protein interaction analysis
The protein-protein interaction (PPI) network of the DEPs was shown in Fig. 6. In the PPI networks, nodes represent the identified DEPs and links represent correlation among different nodes. Proteins with high connectivity may be key points that affect metabolism or signal transduction pathways throughout the system. Data analysis indicated that some proteins in PPI networks could play key roles in SJZD intervention spleen deficiency 14 syndrome and these included H2AFX, ISG15 and STAT1.

Western blot verification of the differentially expressed proteins
Three proteins H2AFX, ISG15 and STAT1 were selected to validate the results of label-free proteomics by performing western blot analyses. The results showed the expression levels of three proteins were significantly different between the M and NC groups or M and SJZD groups (Fig. 7A). We further analysed the bands based on the gray-level values of the three proteins. H2AFX was significantly upregulated in M group compared with the NC and SJZD groups, while ISG15 and STAT1 were significantly downregulated, consistent with the proteomics results (Fig. 7B).

Discussion
In this study, the results of body weight and viscera indices, enzyme activity assays, histopathological observations and proteomic analysis suggest that Sijunzi decoction has a protective effect on rats with spleen deficiency syndrome. Most important of all, the proteomic analysis revealed the potential therapeutic mechanism of Sijunzi decoction in rats with spleen deficiency syndrome.
In addition to being closely related to the digestion and absorption of nutrients, the function of the spleen is also closely related to material and energy metabolism. Studies have reported that the activity of enzymes related to energy metabolism in patients with spleen deficiency syndrome differ from that of healthy individuals [33]. SDH is one of the key enzymes that catalyzes oxidative phosphorylation reactions in the inner mitochondria of cells. The first dehydrogenase in the tricarboxylic acid cycle in the central metabolism of the body reflects the energy supply of oxidative phosphorylation in cells. Xue et al.
found that SDH activity in spleen deficiency rats was significantly decreased, and was one of the important causes of digestive dysfunction [34]. LDH is the most active glycolytic enzyme in the liver, and also exists in the cytoplasm of all tissues in the body. LDH can catalyse lactic acid transformation into pyruvate, and therefore it can be used as a standard to predict the anaerobic metabolic capacity of the body. In this study it was found that compared with the NC group, the SDH activity in the liver of the M group was significantly decreased (p value < 0.05), while the activity of LDH was significantly increased (p value < 0.05). This indicates that aerobic oxidation metabolism was decreased, and anaerobic oxidation reaction was increased in the liver tissues of the rats with spleen deficiency syndrome. This is consistent with previous studies where increased LDH activity was observed in rats with spleen deficiency syndrome [9]. After treatment with Sijunzi decoction, the enzyme activities recovered back to normal levels. This suggests that Sijunzi decoction can significantly aid in restoring mitochondrial health caused by spleen deficiency syndrome in rats and thereby improve the energy supply level in the body.
In this study, pathological observations were carried out on liver tissues of model rats, and it was found that the liver tissues of rats with spleen deficiency syndrome had varying degrees of pathological changes, consistent with previous studies [35]. Additionally, these results are consistent with previous research on the pathological changes in spleen deficiency syndrome such as congestion, cell degeneration, atrophy and inflammatory cell infiltration in rat liver tissues [36].
Label free proteomics technology was employed to identify differentially expressed proteins, to study the mechanism of action of Sijunzi decoction in rats with spleen deficiency syndrome. In this study, 272 DEPs (128 proteins up-regulated, 144 proteins down-regulated) in rats with spleen deficiency syndrome were found when comparing the protein expression spectrum of the M group with the NC group. When comparing the M group to the SJZD group, 247 DEPs (81 proteins upregulated, 166 proteins downregulated) were found. Venn diagram showed there were 98 common DEPs. At the same time, we found that 81 proteins showed a reversal of the trend in the SJZD group, indicating Sijunzi decoction can adjust most of the DEPs and adjust their expression back to normal. It can be speculated that the common DEPs are potential target proteins of Sijunzi decoction when treating spleen deficiency syndrome. GO enrichment analysis showed that the DEPs in Sijunzi decoction intervention rats were mainly distributed in the cytoplasm, cell membrane and extracellular region. Additionally, the main function of the proteins was found to be in the activity of enzyme catalysis, transformation and transfer as well as other various metabolic processes. By comparing the target pathways in rats with spleen deficiency syndrome and Sijunzi decoction intervention rats with spleen deficiency syndrome, it was revealed that the mechanism of action of Sijunzi decoction was mainly through the regulation of four pathways (bile secretion, pyruvate metabolism, cysteine and methionine metabolism, and fatty acid metabolism).
The bile secretion pathway is an important secretion pathway in the body. The secretion and excretion of bile are regulated by a variety of bodily fluids, with gastrin promoting the secretion of bile [37]. Our pharmacodynamic studies showed that gastrin content in the blood of rats with spleen deficiency syndrome decreased significantly compared with that of the NC group (p value < 0.05). This suggests that the decrease of gastrin inhibits bile secretion, hindering the digestion and absorption of fats and vitamins. Fat metabolism in rats with spleen deficiency syndrome was restrained, leading to a shortage of nutrients and energy for the body, and consequently, rats with spleen deficient syndrome appeared thin and tired. This is consistent with the results of pharmacodynamics studies [38]. In the bile secretion pathway, the expression of solute carrier organic anion transporter family member 1A4 and bile salt export pump were up-regulated. The expression of ATP-binding cassette and neutral cholesterol ester hydrolase 1 were down-regulated. This indicates that bile secretion is dysfunctional in the liver of rats with spleen deficiency syndrome, leading to the inhibition of digestion and absorption functions. After Sijunzi decoction intervention intervention, gastrin content was found to be increased and bile secretion returned to normal. Additionally, DEPs (ATP-binding cassette, neutral cholesterol ester hydrolase 1) returned to normal, so that bile secretion in the liver was closer to normal and digestion and absorption ability were stable.
The cysteine and methionine metabolism pathway is a very important amino acid metabolism pathway in the body. Methionine and its metabolites have antioxidant capacities and maintain the redox state in cells. Cysteine and methionine metabolism can produce pyruvate in the body through a series of reactions. Pyruvate, the metabolic hub, enters many different metabolic pathways. There were four upregulated DEPs in cysteine and methionine metabolic pathways of the liver of rats with spleen deficiency syndrome (p value < 0.01), including L-serine dehydratase /L-threonine deaminase, aspartate aminotransferase, tyrosine aminotransferase, and S-methyl methionine-homocysteine Smethyltransferase BHMT2. This suggested that the metabolism of cysteine and methionine in liver tissues of rats with spleen deficiency syndrome was enhanced, indicating that the body of rats with spleen deficiency syndrome was under oxidative stress. Oxidative stress damages cells and causes their dysfunction, which leads to viscera damage and all kinds of metabolic disturbances. The expression of these proteins was significantly downregulated after intervention with Sijunzi decoction (p value < 0.05), therefore highlighting that the metabolism of cysteine and methionine was gradually restored to normal in the liver.
Pyruvate acid is an important intermediate product in metabolic processes. Accordingly, pyruvate metabolic disorders disrupt a variety of metabolic processes in the body. In this study, the expression of L-lactate dehydrogenase B chain and acetyl-CoA carboxylase 1 were upregulated in the M group, but the glyoxylate reductase/hydroxypyruvate reductase (GRHPR) protein was downregulated, indicating that the rats with spleen deficiency syndrome were suffering from pyruvate metabolic disorders. It can be speculated that there were three major nutrient disorders in spleen deficiency rats due to the changes in pyruvate metabolism, leading to further deficiency of nutrients and body emaciation. This is also consistent with the results of the physical indicators in this study. In addition, pyruvate is the energy source of mitochondria, and changes to pyruvate metabolism leads to functional changes in the mitochondria and a reduction in the energy supply to the body [39]. This was also consistent with the results of the enzymatic indicators of energy metabolism in this study. After treatment with Sijunzi decoction, pyruvate metabolism was restored to normal, and this then meant that all metabolic processes related to pyruvate metabolism were also restored to normal.
One of the functions of fatty acids in animals is to produce energy in the form of ATP.
Compared with other nutrient carbohydrates and proteins, complete oxidation of fatty acids produces the most ATP energy [40]. In the mitochondria, fatty acids are degraded into acetyl-CoA by β-oxidation. On the other hand, acetyl-CoA can also synthesize fatty acids through the synthesis pathway of fatty acids in the cytoplasm. The two pathways are mutually inhibitory. Studies have showed that acetyl-CoA carboxylase can prevent acetyl-CoA produced from fatty acid β-oxidation entering the fatty acid synthesis pathway [40] .
The expression of acetyl-CoA carboxylase was found to be up-regulated in the tissues of liver from the M group in this study, and the synthesis of fatty acids was inhibited.
Therefore, due to the body's demand for energy, the breakdown of fatty acids was increased and fat was consumed in much larger quantities, resulting in the emaciation of the body. After treatment with Sijunzi decoction intervention, the expression of acetyl-CoA carboxylase decreased and fatty acid metabolism was returned to normal.
By constructing protein interaction networks represented by nodes and links, effective information of proteins can be extracted from different system scales and comprehensive information that cannot be obtained by a single protein is obtained. Potential protein targets of Sijunzi decoction intervention treatment spleen deficiency syndrome were predicted by PPI networks. ISG15, one of our potential protein targets, is known to be involved in the immune response, and has been proved to exert antiviral activity against a variety of viruses [41]. STAT1 is necessary for congenital immunity of the body. In addition, it can also act as a cell growth inhibitor and apoptosis activator. In this study, the expression of ISG15 and STAT1 were found to be downregulated in the liver tissue of M group rats, resulting in immune system dysregulation. Protein H2AFX levels are associated with damage to DNA in cells [42]. The expression levels of H2AFX were found to be upregulated in M group rats, which suggested the relevant DNA in the cell may be damaged in rats with spleen deficiency syndrome. After treatment with Sijunzi decoction, the expression of ISG15, STAT1 and H2AFX were called back. ISG15, STAT1 and H2AFX were detected by western blots and further verified the results of proteomics.

Conclusions
In conclusion, Current study demonstrated that Sijunzi decoction could remarkably improve spleen deficiency syndrome induced by reserpine. Liver proteomics speculated that bile secretion, pyruvate metabolism, cysteine and methionine metabolism, and fatty

Acknowledgements
We thank International Science Editing for providing language editing services for this manuscript.

Publication statement
The contents of this manuscript have not been published or submitted for publication elsewhere.

Consent for publication
Not applicable.

Authors' contributions
Qinghua Yan designed study, performed experiments, conducted data analysis and prepared manuscript. Shijiang Ma, contributed to experimental design and assisted with data interpretation. Qiqing Zhang helped with study design, data interpretation and manuscript preparation. Hongli Chen contributed to data interpretation and manuscript preparation. All authors read and approved the finally manuscript.

Institutions of Henan Province (20A180002) and the Doctoral Scientific Research
Foundation of Xinxiang Medical University (505304).

Availability of data
The mass spectrometry proteomics data have been deposited to the ProteomeXchange 21 Consortium (http://proteomecentral. proteomexchange.org) via the iProX partner repository with the dataset identifier PXD015449. The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.