Skip to main content

Quantitative proteomics reveals the effect of Yigu decoction (YGD) on protein expression in bone tissue

Abstract

Background

Osteoporosis (OP) is a systemic bone disease characterized by decreased bone mass, destruction of the bone tissue microstructure, increased bone brittleness and an increased risk of fracture. OP has a high incidence rate and long disease course and is associated with serious complications. Yigu decoction (YGD) is a compound prescription in traditional Chinese medicine that is used to treat OP. However, its mechanism in OP is not clear. This study used a tandem mass tag (TMT)quantitative proteomics method to explore the potential bone-protective mechanism of YGD in an osteoporotic rat model.

Materials and methods

A rat model of OP was established by ovariectomy. Eighteen 12-week-old specific-pathogen-free female Wistar rats weighing 220 ± 10 g were selected. The eighteen rats were randomly divided into 3 groups (n = 6 in each group): the normal, model and YGD groups. The right femurs from each group were subjected to quantitative biological analysis. TMT quantitative proteomics was used to analyze the proteins extracted from the bone tissue of rats in the model and YGD groups, and the differentially expressed proteins after intervention with YGD were identified as biologically relevant proteins of interest. Functional annotation correlation analysis was also performed to explore the biological function and mechanism of YGD.

Result

Compared with the model group, the YGD group showed significant upregulation of 26 proteins (FC > 1.2, P < 0.05) and significant downregulation of 39 proteins (FC < 0.833, P < 0.05). Four important targets involved in OP and 5 important signaling pathways involved in bone metabolism were identified.

Conclusions

YGD can significantly increase the bone mineral density (BMD) of osteoporotic rats and may play a therapeutic role by regulating target proteins involved in multiple signaling pathways. Therefore, these results improve the understanding of the OP mechanism and provide an experimental basis for the clinical application of YGD in OP treatment.

Background

Osteoporosis (OP) is a bone disease characterized by low bone mass, bone tissue destruction and pain [1]. The main manifestations are low bone mass, bone tissue metamorphism, destruction of the bone structure, reduced bone strength and an increased risk of fracture [2]. Low bone mass is characterized by decreased bone mineral density (BMD) and increased bone brittleness. Bone tissue metamorphism and bone structure destruction are bone homeostasis changes caused by decreased osteoblast production and increased osteoclast proliferation, and fracture, mainly vertebral fracture, hip fracture or forearm fracture, is the end point event of OP. Fracture caused by OP increases the medical burden [3, 4], and refracture [5,6,7] and death after fracture are also harmful to public health, which necessitates studies on the etiology of OP. Some studies have attempted to explain the pathogenesis of OP via different basic research methods and identify therapeutic targets by proteomic and metabonomic investigation of the pathogenic mechanism of OP. Previous studies [8,9,10,11,12] have shown that a variety of signaling pathways, such as the ER pathway, OPG/RANKL pathway, BMP-2/Smad pathway, classical Wnt/β-catenin signaling pathway, and PI3K/akt signaling pathway, can affect osteoblast proliferation and inhibit osteoclast differentiation.

During the development of modern medicine, due to the influence of cultural exchange and integration, traditional Chinese medicine (TCM) has attracted increasing attention because of its unique role. Researchers have extracted active ingredients from various herbs to treat different diseases. Yigu decoction (YGD) consists of Fructus psoraleae, Rhizoma Drynariae, Radix Rehmanniae, Salvia miltiorrhiza, Epimedii Folium and Dioscorea opposita. Fructus psoraleae could reportedly improve the damage caused by bone resorption at the level of bone metabolism [13] and increase BMD in mice [14]. Studies have confirmed that Rhizoma Drynariae has antioxidant effects on osteoblasts [15]. The use of Radix Rehmannia extract increased the thickness of cortical bone and the number of trabeculae in the bone marrow cavity of osteoporotic rats [16]. According to previous studies, the extract of Salvia miltiorrhiza could improve BMD [17] and increase the levels of alkaline phosphatase and tartrate resistant acid phosphatase (TRACP) [18]. The extract of Epimedii Folium can significantly increase the expression of transforming growth factor-beta (TGF-β), reduce the movement and bone resorption activity of isolated osteoclasts [19], accelerate bone formation, and contribute to bone plastic surgery and remodeling [20]. A literature search revealed that Dioscorea opposita is beneficial for the proliferation of osteoblasts and inhibits the differentiation of osteoclasts [21]. In basic research, our team found that YGD could increase the proliferation and differentiation of osteoblasts [22, 23], increase the serum levels of bone Gla protein (BGP) and Ca2+ [24], improve the biomechanical properties of bone [25, 26], and promote the expression of osteoblast BMP [25, 27] to ameliorate the inhibitory effect on bone resorption and achieve the prevention and treatment of OP. The effects of YGD in the treatment of OP are mainly achieved by promoting osteoblast proliferation. Our previous research [26,27,28,29,30] mainly focused on the effect of YGD on the classical Wnt/β-catenin and BMP-2/Smad signaling pathways in the bone tissue of osteoporotic rats. However, whether YGD is also involved in the regulation of other proteins and the expression of signaling pathway components warrants further study. Microarray studies have been performed to explore the pathogenesis of OP. A larger-sample size expression profile can be obtained by integrated microarray analysis [37], which is helpful for identifying the differentially expressed genes (DEGs) in OP more accurately than a single microarray. Based on the above research, we established an OP model in ovariectomized rats and treated them with YGD. The use of TMT quantitative proteomics technology to study protein expression differences, mass spectrometry experiments and data analysis led to possible conclusions, which guided subsequent research.

Materials and methods

Drugs and reagents

YGD contains six kinds of traditional Chinese medicine: Rhizoma Drynariae, Radix Rehmanniae, Fructus psoraleae, Salvia miltiorrhiza, Epimedii Folium, and Dioscorea opposita, which were purchased from The Third Affiliated Hospital of Zhejiang Chinese Medical University.

Animal culture

Eighteen12-week-old specific-pathogen-free female Wistar rats weighing 220 ± 10 g were used in this study. All the experimental animals were provided by Shanghai Shrek Experimental Animal Co., Ltd. (license number: SCXK, Shanghai 2013–0016, certificate code: 2013001806499). The rats were raised in the Animal Experimental Center of Zhejiang University of Traditional Chinese Medicine (laboratory facility license number: SYXK [Zhejiang] 2013–0184). The feeding room was well ventilated, the room temperature was controlled at 21 ± 1 °C, the humidity was 63%, the noise was less than 55 decibels, food and water were provided ad libitum, and light/dark conditions were alternated every 12 h. All animals received meticulous and humanitarian care.

Eighteen rats were randomly divided into 3 groups (6 rats in each group): the normal group, model group and YGD group. Under anesthesia via an intraperitoneal injection of ketamine (5 mg 100 g−1), rats in the model group and YGD group were ovariectomized. Three days after the operation, each rat was injected with 40,000 U/d penicillin intramuscularly to prevent infection. The rats in the normal group did not receive any treatment. There was no death of rats after the operation, and the success rate of the model was 100%. After the operation, the rats were fed a routine diet for 12 weeks, and the drug was given once a day from the 13th week. The rats in the normal group and the model group were perfused with 10 ml·kg−1 double distilled water, and the rats in the YGD group were fed 10 ml kg−1 water extract of YGD. The whole process lasted for 12 weeks.

After 12 weeks, the rats were anesthetized via injection of 10% pentobarbital sodium in the right femur, and the left ventricle was intubated and perfused with normal saline. All bone tissues were quickly frozen and preserved at − 80 °C.

Determination of BMD in rats

A dual-energy X-ray bone densitometer was used to measure the BMD of the small animals (produced by GE, USA, Lunar Prodigy). The femur was placed on the measuring table for automatic measurement.

HE staining

The bone tissue of the rats was removed, fixed with 4% paraformaldehyde (in 0.1 MPBS containing 0.1% DEPC, pH 7.0–7.6), and decalcified with 10% EDTA. The distal femoral intercondylar fossa was cut lengthwise into a 1.5 cm bone segment. Tissue samples were processed by conventional paraffin embedding, sectioning, processing, and HE staining and examined by light microscopy.

Protein extraction and peptide enzymatic hydrolysis

The rats were anesthetized with 10% pentobarbital sodium at 0.3 ml 100 g−1. The blood was collected from the heart in vitro, placed at room temperature for 40 min at 4 °C, and centrifuged at 3600 rpm for 15 min. The upper yellowish serum was transferred into a new centrifuge tube. The sample was filtered through a 0.22 µm membrane, sterilized on an ultraclean worktable and then placed in a freezer at -80 °C for long-term storage.

SDT (4% SDS, 100 mM Tris–HCl, pH 7.6, 0.1 M DTT) buffer was used for sample lysis and protein extraction. The amount of protein was quantified with the BCA Protein Assay Kit (Bio-Rad, USA). Protein digestion by trypsin was performed according to the filter-aided sample preparation (FASP) procedure described by Matthias Mann. The digested peptides of each sample were desalted on C18 Cartridges (Empore SPE Cartridges C18 [standard density], bed I.D. 7 mm, volume 3 ml, Sigma), concentrated by vacuum centrifugation and reconstituted in 40 µl of 0.1% (v/v) formic acid (FA).

TMT labeling

For each sample, a 100 μg of peptide mixture was labeled using TMT reagent according to the manufacturer’s instructions (Thermo Scientific). The tagging information is shown in the table below. There were 3 groups designed for the project, and each group contained 3 or 2 biological repetitive samples, for a total of 8 samples (Table 1).

Table 1 TMT labeling scheme for the 8 samples

Reversed-phase (RP) grading

The labeled peptides in each group were mixed in equal amounts and graded by a High pH Reversed-Phase Peptide Fractionation Kit. First, acetonitrile and 0.1% trifluoroacetic acid (TFA) were used for column equilibration, and then, the mixed labeled peptide samples were added to the column with pure water and desalted at low speed by centrifugation. Finally, the column-bound peptides were gradient eluted with high-pH acetonitrile solutions of progressively increasing concentrations. After vacuum drying, 12 μl of 0.1% FA was used to dissolve the freeze-dried sample, and the peptide concentration was determined by measuring the OD280.

LC–MS/MS data acquisition

LC–MS/MS analysis was performed on a Q Exactive mass spectrometer (Thermo Scientific) that was coupled to an Easy nLC (Proxeon Biosystems, now Thermo Fisher Scientific) for 90 min. The peptides were loaded onto a reversed-phase trap column (Thermo Scientific Acclaim PepMap100, 100 μm * 2 cm, nanoViper C18) connected to a C18 reversed-phase analytical column (Thermo Scientific Easy Column, 10 cm long, 75 μm inner diameter, 3 μm resin) in buffer A (0.1% FA) and separated with a linear gradient of buffer B (84% acetonitrile and 0.1% FA) at a flow rate of 300 nl/min controlled by IntelliFlow technology. The mass spectrometer was operated in positive ion mode. MS data were acquired using a data-dependent top10 method by dynamically choosing the most abundant precursor ions from the survey scan (300–1800 m/z) for HCD fragmentation. The automatic gain control (AGC) target was set to 3e6, and the maximum injection time was set to 10 ms. The dynamic exclusion duration was 40.0 s. Survey scans were acquired at a resolution of 70,000 at 200 m/z, the resolution for the HCD spectra was set to 17,500 at 200 m/z, and the isolation width was 2 m/z. The normalized collision energy was 30 eV, and the underfill ratio, which specifies the minimum percentage of the target value likely to be reached at the maximum fill time, was defined as 0.1%. The instrument was run with peptide recognition mode enabled.

Data analysis

The MS raw data for each sample were searched using the MASCOT engine (Matrix Science, London, UK; version 2.2) embedded into Proteome Discoverer 1.4 for identification and quantitation analysis.

Bioinformatics analysis

Cluster 3.0 (http://bonsai.hgc.jp/~mdehoon/software/cluster/software.htm) and Java Treeview (http://jtreeview.sourceforge.net) were used to perform hierarchical clustering analysis; the Euclidean distance algorithm was selected for similarity measurement and the average linkage clustering algorithm for clustering. A heat map is often presented as a visual aid in addition to the dendrogram.

Blast2GO (https://www.blast2go.com/) was used to annotate the target protein set with GO, and the process can be summarized into four steps: sequence alignment (Blast), GO entry extraction (mapping), GO annotation (annotation) and InterProScan supplementary annotation (annotation augmentation).

The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway of the target protein set was annotated by KAAS (KEGG Automatic Annotation Server) software. All differentially expressed proteins were compared to all of the experimentally identified proteins with KEGG annotation results to reveal the enriched pathways as determined by Fisher's exact test.

Based on the PPI relationship in the STRING database, a PPI network map was constructed for the differentially expressed proteins in the comparison group by using CytoScape software. The direct and indirect interactions between the target proteins were found based on the information in the STRING (http://string-db.org/) database, and the interaction network was generated and analyzed by using CytoScape software (version number: 3.2.1).

Statistical analysis

All the data were from at least three independent experiments and are expressed as the mean ± standard deviation (x ± s). The univariate comparison between the two groups was analyzed by Student's T-test in GraphPad Prism software, and the multivariate data were compared by analysis of variance. The difference was statistically significant (P < 0.05).

Results

The ability of YGD to alleviate OP

To explore the ability of YGD to alleviate OP, we first constructed a rat model of OP and measured bone density. The results showed that the bone density value of the model group was significantly lower than that of the normal group (P < 0.01), which indicated that the OP rat model was successfully established. Compared with that of the model group, the BMD of the YGD group increased significantly, and the difference was statistically significant (P < 0.01) (Table 2). In parallel, our team observed the pathological histomorphology of the left femur of the rats. The visible area of the trabecular bone was smaller, the width was narrower, the cortex was thinner, the number was reduced, the morphology and structure were irregular, the arrangement was disordered, and the marrow was expanded in the model group compared with the normal group. After YGD administration for 3 months, the rats in the YGD group were similar to those in the blank group. The bone trabecula was slightly thinner, the shape remained regular, the arrangement remained orderly, and the gap was slightly enlarged in some areas (Fig. 1).

Table 2 Comparison of bone mineral density results of rats in each group (x ± s) (unit: g/cm2)
Fig. 1
figure 1

Pathological observation of the right femur of rats (HE staining, 100×)

TMT quantitative proteomic analysis of differentially expressed proteins after YGD treatment

To further explore the mechanism of action of YGD in the treatment of OP, through TMT quantitative proteomics analysis, a total of 340,987 secondary spectra were obtained from bone tissues from the control group (control), OP model group (model) and OP model intervention group (YGD). The total number of spectra with matches in the database was 29,695, including 3382 proteins, and 3363 matches showed high similarity (Additional file 1, Fig. 2).

Fig. 2
figure 2

Statistical picture of identification and quantitative results

Bioinformatics analysis of differentially expressed proteins after YGD treatment

To analyze the differentially expressed proteins between different groups, the experimental data were further screened. Biological analysis of the differentially expressed proteins showed that a total of 65 proteins were significantly affected after treatment. Compared with the model group, the YGD group showed significant upregulation of 26 proteins and downregulation of 39 proteins (FC ≥ 1.2, P < 0.05). To determine whether the differential protein expression changes represent significant effects caused by biological treatment of the samples, a hierarchical clustering algorithm was used to group the differentially expressed proteins of the compared groups and present them as a heat map. The heat map of the results shows that the significantly differentially expressed proteins can effectively separate the comparison groups, indicating that the screening of differentially expressed proteins can represent the impact of biological treatment on the samples. The detailed analysis is shown in Fig. 3. (Please refer to Additional file 2 for the specific proteins.)

Fig. 3
figure 3

TMT quantitative proteomic analysis of differentially expressed proteins after YGD treatment. a The volcano plot was drawn using two factors, log2 fold change (log2FC) and P value (Logarithmic transformation based on 10, Student's t-test). The significantly downregulated proteins are annotated in blue (FC < 0.833 and P < 0.05), the significantly upregulated proteins are annotated in red (FC > 1.2 and P < 0.05), and the proteins without differences are indicated in gray. b The hierarchical clustering results are presented as a tree-type heatmap with the abscissa showing the sample information and the ordinate showing the significantly differentially expressed proteins, and the expression amounts of the significantly different proteins in different samples were exhibited in the Heatmap by different colors after normalization using the log2 method. Where red represents the significantly upregulated proteins, blue represents the significantly downregulated proteins and grey represents proteins with no quantitative information

Functional analysis of differentially expressed proteins

To better understand the biological functions of differentially expressed proteins after YGD treatment, we carried out functional annotation via GO and KEGG analysis. A total of 1656 GO terms were obtained by GO functional enrichment analysis, including 1195 biological process (BP) entries, 218 cellular component (CC) entries and 243 molecular function (MF) entries. The first 20 items are shown in Fig. 4a–c). Significant differences are observed in biological processes related to OP, such as protein regulation and muscle regulation, the molecular function of enzyme regulation, and the localization protein collagen trimer.

Fig. 4
figure 4figure 4figure 4

a Go functional enrichment bubble plot under biological process classification. b Bubble plot of go functional enrichment under cell component classification. c Bubble plot of go functional enrichment under molecular function classification

To analyze the biological process and drug action mechanism more systematically and comprehensively, we performed KEGG pathway enrichment analysis of the 65 differential proteins and found that 93 signaling pathways were changed significantly. The main signaling pathways involved included the protein digestion and absorption signaling, vitamin digestion and absorption signaling, Thyroid hormone synthesis (TH signaling) and TGF-beta (TGF-β) signaling pathways. The protein digestion and absorption pathways were visualized and analyzed, and the top 20 signaling pathways are shown in Fig. 5.

Fig. 5
figure 5

Bubble diagram of the KEGG pathway (top20) enrichment analysis of differentially expressed proteins

The abscissa in the graph is the enrichment factor (rich factor ≤ 1), which represents the proportion of the number of differentially expressed proteins annotated to a certain go functional category over the number of all identified proteins annotated to that go functional category.The ordinate represents the differential protein statistical outcome under each go functional classification; Where bubble color indicates the significance of enriched go functional categories, i.e., the p value was calculated based on Fisher's exact test, color gradient represents the magnitude of the p value (taken—log 10), and closer color to red represents smaller p value, corresponding to higher significance level of enrichment of go functional categories.

The abscissa of the plot shows the enrichment factor (richness factor ≤ 1). The enrichment factor indicates the number of differentially expressed proteins annotated to a GO functional category as a proportion of the number of all identified proteins annotated to the GO functional category. The ordinate represents the differential protein statistics under each KEGG pathway, and the bubble color represents the significance of the enriched KEGG pathways. Colors closer to red represent smaller p values, corresponding to a higher level of significance for metabolic pathway enrichment. The size of the circles indicates the number of proteins represented.

Protein–protein interaction (PPI) analysis

An important mechanism by which proteins function is via interactions with other proteins and regulation of biological processes through protein-mediated pathways or the formation of complexes. According to PPI network analysis, a total of 21 among the 65 differentially expressed proteins may directly interact. Alb, DCN, Col6a2, and Smad2 were involved in many ways. Interestingly, we found that the differentially expressed protein DPP4 is also in figure, which is a protein closely related to diabetes. This protein may also be involved in and affect bone metabolic pathways. (see Fig. 6).

Fig. 6
figure 6

Interaction network analysis

The nodes are proteins, the lines represent functional associations between the proteins, and the different line colors represent the types of evidence for the predicted functional associations.

Discussion

OP is defined as a hip or lumbar BMD less than or equal to 2.5 standard deviations of the average BMD of the young reference population [31]. According to the guidelines of the International Osteoporosis Foundation [32], patients with osteopenia or OP should be treated with vitamin D combined with calcium [33], and patients with severe cases need supplementary drug intervention [2]. The use of calcium agents such as bisphosphate can improve the condition of patients with OP to some extent. However, these drugs have limited patient compliance and cause long-term adverse reactions [34], such as increased bone mineral density without improved bone quality or bone strength and effects on calcium and phosphorus metabolism [35, 36]. YGD treatment was shown to be effective, but the mechanism remains unclear. Quantitative mass spectrometry of proteins showed a possible relationship between pathways such as protein digestion and absorption signaling pathway obtained from GO enrichment analysis and bone metabolism, and protein network interaction showed that proteins such as albumin (Alb) have a relatively close relationship with OP. Thus, YGD may function through the process described above.

The normal digestion and absorption of proteins is fundamental for maintaining human life [38]. A variety of trace elements needed by the human body, such as vitamin B and cobalamin, must be absorbed through the intestine, and the absorption of cobalamin needs to be completed by binding proteins [39]. The protein digestion and absorption signaling pathway is mainly involved in the effective digestion and absorption of intestinal protein and fat [40]. Studies have shown that the effective digestion and absorption of proteins may be beneficial to cell metabolism, and high-quality cell status can help maintain immunity and functional stability in the body [41]. Dietary proteins participate in the physiological functions of the intestine to control amino acid, glucose and lipid metabolism so that bone cells can maintain normal physiological functions [42]. The protein digestion and absorption signaling pathway may play a role mainly through intestinal flora [43, 44], thus affecting bone calcium metabolism. By downregulating the key proteins in this signaling pathway, namely, Col6a2 and Col4a1, the protein DPP4 may be upregulated to promote bone formation and inhibit bone resorption. This confirmed the accuracy of our study that col6a2 and Col4a1 proteins were suppressed and that DPP4 protein was upregulated. There is a positive correlation between the vitamin digestion and absorption signaling pathways and aging [45]. Its relationship with protein digestion and absorption signal pathway is synergistic. With increasing age, the absorption of vitamins through the intestinal flora decreases significantly, the most important change is the decrease in vitamin D absorption, and this decrease in vitamin D absorption in turn leads to the occurrence of OP. In the vitamin digestion and absorption signaling pathway, we found that vitamin B was also involved [46]. The existing conclusions also confirm our findings that vitamins B1, B12 (cobalamin) and B9 (folic acid) have bone protective effects [47, 48]. Vitamin B12 deficiency can increase serum homocysteine levels, leading to OP and osteoporotic fractures [49]. In this study, we believe that the use of TCM may regulate the content of B vitamins in rats and protect bone tissue through this pathway. TH is directly involved in the development of osteocytes and controls the linear growth and maturation of bones [50], which requires protein to provide enough power. TGF-β is the main transforming growth factor, and the disturbance of TGF-β signaling is the basis of inflammatory diseases [51], which play an important role in OP [52,53,54]. This signaling pathway was also mentioned in our previous research [18]. Therefore, we believe that YGD may directly and indirectly participate in the regulation of this signaling pathway, thus affecting bone metabolism, and thereby controlling the occurrence and development of OP to a certain extent.

Interactive network analysis showed that Alb participates in a variety of processes, and related studies [55,56,57] have shown that the protein is closely related to energy and bone metabolism. Alb is abundant in human blood and is mainly involved in the transport of various substances. Alb also could enhance bone formation [58]. It has been reported that there is a direct relationship between serum Alb content and OP, and a low protein diet can lead to a decrease in bone density and eventually trigger OP [71]. According to our experimental findings, downregulation of Alb protein may be through the TH signaling pathway and thereby inhibit bone destruction [59, 69]. Through downregulating Alb thereby inhibited the expression of osteoclastogenic proteins and stimulated osteoclast apoptosis [72]. For preventing bone resorption. Decorin (DCN) is one of several glycosaminoglycan (GAG) no-collagen proteins secreted by bones. It is characterized as a leucine-rich repetitive core protein in bone [60]. DCN can directly interact with TGF-β in a dynamic environment to affect cell migration [61]. According to our study, the downregulation of DCN protein may be through TGF-β Signaling pathway inhibited the fibrotic response in skeletal muscle [70]. By decreasing the content of DCN, the acceleration of bone mineralization can be achieved [63]. Bone loss can be reduced by inhibiting the DCN protein. It can also promote the activity of osteoblasts [62]. Smad is the main transducer of bone morphogenetic protein (BMP) and the TGF-β signaling pathway that regulates mitochondrial function and apoptosis [64]. According to our results, the downregulation of Smad proteins similarly confirms that YGD may through TGF-β signaling pathway inhibits bone remodeling and reduces osteoclast activity [65]. Similarly, Smad inhibits the apoptosis of mature osteoblasts through BMP2 [66]. Interestingly, we found that DPP4 expression was significantly increased among the proteins regulated by YGD, and this protein was closely related to diabetes [67, 68]. The high expression of DPP4 in OP may indicate that different signaling pathways have beneficial effects on bone formation and bone resorption.

Conclusions

Using TMT-based proteomics technology and proteome microarray technology, the effect of YGD on protein expression in bone tissue was studied. The results showed that YGD participated in several signaling pathways related to OP and some proteins involved in the regulation of signaling pathways. This confirmed the conclusions reached in some reports and proved our research hypothesis regarding YGD. That is, that YGD participates in the regulation of different target proteins, is involved in signal transduction via multiple signaling pathways and plays a therapeutic role. At the same time, we also found that there is a certain relationship between the differentially expressed protein DPP4 and OP, but how DPP4 participates in signal pathway transduction and affects the occurrence and development of OP is worthy of further study.

Availability of data and materials

The data are included in the article as figures and tables. All data generated during this study are included in this published article.

Abbreviations

OP:

Osteoporosis

YGD:

Yigu decoction

TCM:

Traditional Chinese medicine

TMT:

Tandem mass tag

TGF-β:

Transforming growth factor-beta

TRACP:

Tartrate resistant acid phosphates

BGP:

Bone Gla protein

DEGs:

Differentially expressed genes

FA:

Formic acid

FASP:

Filter-aided sample preparation

BP:

Biological process

CC:

Cellular

MF:

Molecular function

PPIs:

Protein–protein interactions

TH:

Thyroid hormone synthesis

ECM:

Extracellular matrix

DCN:

Decorin

GAG:

Glycosaminoglycan

KEGG:

Kyoto Encyclopedia of Genes and Genomes

KAAS:

KEGG Automatic Annotation Server

BMD:

Bone mineral density

Alb:

Albumin

References

  1. Consensus development conference. prophylaxis and treatment of osteoporosis. Osteoporos Int. 1991;1(2):114–7.

    Article  Google Scholar 

  2. Cosman F, de Beur SJ, LeBoff MS, Lewiecki EM, Tanner B, Randall S, et al. Clinician’s guide to prevention and treatment of osteoporosis. Osteoporos Int. 2014;25(10):2359–81.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  3. Burge R, Dawson-Hughes B, Solomon DH, Wong JB, King A, Tosteson A. Incidence and economic burden of osteoporosis-related fractures in the United States, 2005–2025. J Bone Miner Res. 2007;22(3):465–75.

    Article  Google Scholar 

  4. Lin H, Sohn J, Shen H, Langhans MT, Tuan RS. Bone marrow mesenchymal stem cells: Aging and tissue engineering applications to enhance bone healing. Biomaterials. 2019;203:96–110.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  5. Lewiecki EM, Laster AJ. Clinical review: Clinical applications of vertebral fracture assessment by dual-energy x-ray absorptiometry. J Clin Endocrinol Metab. 2006;91(11):4215–22.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  6. Abrahamsen B, van Staa T, Ariely R, Olson M, Cooper C. Excess mortality following hip fracture: a systematic epidemiological review. Osteoporos Int. 2009;20(10):1633–50.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  7. Colón-Emeric C, Kuchibhatla M, Pieper C, William H, Lisa F, Jay M, et al. The contribution of hip fracture to risk of subsequent fractures: data from two longitudinal studies. Osteoporos Int. 2003;14(11):879–83.

    PubMed  Article  PubMed Central  Google Scholar 

  8. Pang X, Li Y, Bao N, Cong Y. Research progress of main signal pathways and signal molecules of bone metabolism. Basic Medicine and Clinic. 2018;38(12):1799–2180.

    Google Scholar 

  9. Xiong Y, Zhou Y, Lei T. Research progress of bone metabolic signal pathway. Chin J of Ost. 2014;20(02):200–4.

    CAS  Google Scholar 

  10. Wang K, Chen Y, Gao S, Maosi W, Mengmeng G, Qian Y, et al. Norlichexanthone purified from plant endophyte prevents postmenopausal osteoporosis by targeting ER α to inhibit RANKL signaling. Acta Pharm Sin B. 2021;11(2):442–55.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  11. Mo H, Zhang N, Li H, Li F, Pu R. Beneficial effects of Cuscuta chinensis extract on glucocorticoid-induced osteoporosis through modulation of RANKL/OPG signals. Braz J Med Biol Res. 2019;52(12):e8754.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  12. Wang Y, Han X, Zang T, Kang P, Jiang W, Niu N. miR-29b enhances the proliferation and migration of bone marrow mesenchymal stem cells in rats with castration-induced osteoporosis through the PI3K/AKT and TGF-β/Smad signaling pathways. Exp Ther Med. 2020;20(4):3185–95.

    CAS  PubMed  PubMed Central  Google Scholar 

  13. He J, Li X, Wang Z, Samuel B, Kai C, Zhifeng X, et al. Therapeutic anabolic and anticatabolic benefits of natural chinese medicines for the treatment of osteoporosis. Front Pharmacol. 2019;10:1344.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  14. Tsai MH, Huang GS, Hung YC, Bin L, Liao LT, Lin LW. Psoralea corylifolia extract ameliorates experimental osteoporosis in ovariectomized rats. Am J Chin Med. 2007;35(4):669–80.

    PubMed  Article  PubMed Central  Google Scholar 

  15. Sun JS, Lin CY, Dong GC, Sheu SY, Lin FH, Chen LT, et al. The effect of Gu-Sui-Bu (Drynaria fortunei J. Sm) on bone cell activities. Biomaterials. 2002;23(16):3377–85.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  16. Oh KO, Kim SW, Kim JY, Ko SY, Kim HM, Baek JH, et al. Effect of Rehmannia glutinosa Libosch extracts on bone metabolism. Clin Chim Acta. 2003;334(1–2):185–95.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  17. Guo Y, Li Y, Xue L, Richele PS, Sihua G, Jianzhao N, et al. Salvia miltiorrhiza: an ancient Chinese herbal medicine as a source for anti-osteoporotic drugs. J Ethnopharmacol. 2014;155(3):1401–16.

    PubMed  Article  PubMed Central  Google Scholar 

  18. Miao B, Wang J, Zhu Y, Yue C, Chen M. Experimental study on effect of Salvia miltiorrhiza on alveolar bone metabolism and variation in bone mass in diabetic rats. China J Chin Mat Med. 2012;37(11):1659–62.

    Google Scholar 

  19. Huang J, Yuan L, Wang X, Zhang TL, Wang K. Icaritin and its glycosides enhance osteoblastic, but suppress osteoclastic, differentiation and activity in vitro. Life Sci. 2007;81(10):832–40.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  20. Ding Q, Sun N, Wang X. Effect of glycoside of short-horned epimedium Herb on the expressions of transforming growth factor β (1) and basic fibroblast growth factor induced membrane by Masquelet technology. Chin J pro rec sur. 2020;34(11):1405–9.

    Google Scholar 

  21. Yin J, Tezuka Y, Kouda K, Quan LT, Tatsuro M, Yingjie C, et al. Antiosteoporotic activity of the water extract of Dioscorea spongiosa. Biol Pharm Bull. 2004;27(4):583–6.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  22. Bang JH, Yin SZ, Jian WY. Yigu decoction containing serum promotes osteoblast proliferation and differentiation via canonical Wnt signaling pathway. New Chin med. 2017;49(03):10–3.

    Google Scholar 

  23. He BJ, Zhu YS, Ying JW, Li GJ, Lu Y, Yao XM, et al. Study on the effect of Yigu decoction-containing serum on promoting osteoblast proliferation and differentiation through classical Wnt signal pathway. New trad Chin Med. 2017;49(03):10–3.

    CAS  Google Scholar 

  24. Li GJ, Yao XM. Experimental study on the effect of Yigu decoction on the expression of related genes in osteoblasts of osteoporotic rats. Zhejiang J trad Chin Med. 2018;53(10):727–9.

    Google Scholar 

  25. Li W, Yao XM, Li XY, Yu B, Dong ZC. Effect of Yigu decoction on bone mineral density and BMP-2 signal pathway in osteoporotic rats. J Yunnan Col TCM. 2016;39(06):5–9.

    Google Scholar 

  26. Lin XF, Yao XM, Li W, Chen XB, Fang F. Effect of Yigu decoction on classical signal pathway of Wnt/ β-catenin in bone tissue of ovariectomized rats. J Zhejiang Uni TCM. 2018;42(02):97–104.

    Google Scholar 

  27. Yang L, Yao XM, Huang J, Chen YD, Wang J. Experimental study on the effect of Yigu decoction on bone mineral density, biomechanics and BMP-2 expression in ovariectomized rats. J Zhejiang Uni TCM. 2005;(03):53–55.

    Google Scholar 

  28. Yang L, Yao XM, Huang J, Sun QR, Lu JH. Effect of Yigu decoction on serum hormone and osteoblast proliferation in ovariectomized rats. TCM Zhenggu. 2006;8–9(12):79–80.

    Google Scholar 

  29. Yao XM, Chen YD, Fang F. Experimental study on the effect of Yigu decoction containing serum on osteoblast proliferation and ALP. J Zhejiang Uni trad Chin Med. 2007;2:158–9.

    Google Scholar 

  30. Yao XM, Chen YD, Fang F, Shao XM. Effect of medicated serum of Yigu decoction on the expression of osteoprotegerin mRNA in rat osteoblasts. Chin J Ort Trau of TCM. 2007;05:30–2.

    Google Scholar 

  31. Kanis JA, Cooper C, Rizzoli R, Reginster JY, Scientific Advisory Board of the European Society for Clinical and Economic Aspects of Osteoporosis (ESCEO) and the Committees of Scientific Advisors and National Societies of the International Osteoporosis Foundation (IOF). European guidance for the diagnosis and management of osteoporosis in postmenopausal women. Osteoporos Int. 2019;30(1):3–44.

    CAS  Article  Google Scholar 

  32. Baran DT, Faulkner KG, Genant HK, Miller PD, Pacifici R. Diagnosis and management of osteoporosis: guidelines for the utilization of bone densitometry. Calcif Tissue Int. 1997;61(6):433–40.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  33. Liao EY, Zhang ZL, Xia WB, Lin H, Cheng Q, Wang L, et al. Calcifediol (25-hydroxyvitamin D) improvement and calcium-phosphate metabolism of alendronate sodium/vitamin D3 combination in Chinese women with postmenopausal osteoporosis: a post hoc efficacy analysis and safety reappraisal. BMC Musculoskelet Disord. 2018;19(1):210.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  34. Novince CM, Ward BB, McCauley LK. Osteonecrosis of the jaw: an update and review of recommendations. Cells Tissues Organs. 2009;189(1–4):275–83.

    PubMed  Article  PubMed Central  Google Scholar 

  35. Hiligsmann M, Rabenda V, Gathon HJ, Ethgen O, Reginster JY. Potential clinical and economic impact of nonadherence with osteoporosis medications. Calcif Tissue Int. 2010;86(3):202–10.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  36. Ma B, Li X, Zhang Q, Wu D, Wang G, et al. Metabonomic profiling in studying anti-osteoporosis effects of strontium fructose 1,6-diphosphate on estrogen deficiency-induced osteoporosis in rats by GC/TOF-MS. Eur J Pharmacol. 2013;718(1–3):524–32.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  37. Capurso C, Bellanti F, Lo Buglio A, Vendemiale G. The mediterranean diet slows down the progression of aging and helps to prevent the onset of frailty: a narrative review. Nutrients. 2019;12(1):35.

    PubMed Central  Article  CAS  Google Scholar 

  38. Gong R, Ren S, Chen M, Wang Y, Zhang G, Shi L, et al. Corrigendum to “Bioinformatics analysis reveals the altered gene expression of patients with postmenopausal osteoporosis using liuweidihuang pills treatment.” Biomed Res Int. 2019;2019:6897187.

    PubMed  PubMed Central  Google Scholar 

  39. Alpers DH. Absorption and blood/cellular transport of folate and cobalamin: Pharmacokinetic and physiological considerations. Biochimie. 2016;126:52–6.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  40. Jiang Y, Bian Y, Lian N, Wang Y, Xie K, Qin C, et al. iTRAQ-based quantitative proteomic analysis of intestines in murine polymicrobial sepsis with hydrogen gas treatment. Drug Des Devel Ther. 2020;14:4885–900.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  41. Watson MD, Cross BL, Grosicki GJ. Evidence for the contribution of gut microbiota to age-related anabolic resistance. Nutrients. 2021;13(2):706.

    PubMed  PubMed Central  Article  Google Scholar 

  42. Jahan-Mihan A, Luhovyy BL, El Khoury D, Anderson GH. Dietary proteins as determinants of metabolic and physiologic functions of the gastrointestinal tract. Nutrients. 2011;3(5):574–603.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  43. Freeman HJ, Kim YS, Sleisenger MH. Protein digestion and absorption in man. Normal mechanisms and protein-energy malnutrition. Am J Med. 1979;67(6):1030–6.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  44. Wu YR, Kuang GY, Lu FG, Wang HX, Lu M, Zhou Q. Pathological relationship between intestinal flora and osteoarthritis and intervention mechanism of chinese medicine. Chin J Integr Med. 2019;25(9):716–20.

    PubMed  Article  PubMed Central  Google Scholar 

  45. Gallagher JC. Vitamin D and aging. Endocrinol Metab Clin North Am. 2013;42(2):319–32.

    PubMed  PubMed Central  Article  Google Scholar 

  46. Khera S, Pramanik SK, Patnaik SK. Transcobalamin deficiency: vitamin B12 deficiency with normal serum B12 levels. BMJ Case Rep. 2019;12(10):e232319.

    PubMed  Article  PubMed Central  Google Scholar 

  47. Bucciarelli P, Martini G, Martinelli I, Ceccarelli E, Gennari L, Bader R, et al. The relationship between plasma homocysteine levels and bone mineral density in post-menopausal women. Eur J Intern Med. 2010;21(4):301–5.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  48. Ma Q, Liang M, Wang Y, Ding N, Wu Y, Duan L, et al. Non-coenzyme role of vitamin B1 in RANKL-induced osteoclastogenesis and ovariectomy induced osteoporosis. J Cell Biochem. 2020;121(7):3526–36.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  49. Bloise FF, Cordeiro A, Ortiga-Carvalho TM. Role of thyroid hormone in skeletal muscle physiology. J Endocrinol. 2018;236(1):R57–68.

    PubMed  Article  PubMed Central  Google Scholar 

  50. Gouveia CHA, Miranda-Rodrigues M, Martins GM, Neofiti-Papi B. Thyroid hormone and skeletal development. Vitam Horm. 2018;106:383–472.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  51. Batlle E, Massagué J. Transforming growth factor-β signaling in immunity and cancer. Immunity. 2019;50(4):924–40.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  52. Che N, Qiu W, Wang JK, Sun XX, Xu LX, Liu R, et al. MOTS-c improves osteoporosis by promoting the synthesis of type I collagen in osteoblasts via TGF-β/SMAD signaling pathway. Eur Rev Med Pharmacol Sci. 2019;23(8):3183–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  53. Qiu J, Huang G, Na N, Chen L. MicroRNA-214-5p/TGF-β/Smad2 signaling alters adipogenic differentiation of bone marrow stem cells in postmenopausal osteoporosis. Mol Med Rep. 2018;17(5):6301–10.

    CAS  PubMed  PubMed Central  Google Scholar 

  54. Chen X, Wang Z, Duan N, Zhu G, Schwarz EM, Xie C. Osteoblast-osteoclast interactions. Connect Tissue Res. 2018;59(2):99–107.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  55. Pan XY, Zhang ZW. MFGE8, ALB, APOB, APOE, SAA1, A2M, and C3 as novel biomarkers for stress cardiomyopathy. Cardiovasc Ther. 2020;11:1615826.

    Google Scholar 

  56. Shi JM, Yang Z, Li FQ, Wang GJ. Preliminary study of human serum albumin level in early warning onset of preeclampsia. Chin J Obstet Gynecol. 2020;55(1):29–35.

    CAS  Google Scholar 

  57. Luo HC, Fu YQ, You CY, Liu CJ, Xu F. Comparison of admission serum albumin and hemoglobin as predictors of outcome in children with moderate to severe traumatic brain injury: a retrospective study. Medicine (Baltimore). 2019;98(44):e17806.

    Article  Google Scholar 

  58. Goetsch KP, Niesler CU. The extracellular matrix regulates the effect of decorin and transforming growth factor beta-2 (TGF-β2) on myoblast migration. Biochem Biophys Res Commun. 2016;479(2):351–7.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  59. Schandl K, Horváthy DB, Doros A, Majzik E, Schwarz CM, Csönge L, et al. Bone-Albumin filling decreases donor site morbidity and enhances bone formation after anterior cruciate ligament reconstruction with bone-patellar tendon-bone autografts. Int Orthop. 2016;40(10):2097–104.

    PubMed  Article  PubMed Central  Google Scholar 

  60. Verde Z, Giaquinta A, Sainz CM, Ondina MD, Araque AF. Bone mineral metabolism status, quality of life, and muscle strength in older people. Nutrients. 2019;11(11):2748.

    CAS  PubMed Central  Article  Google Scholar 

  61. Licini C, Vitale-Brovarone C, Mattioli-Belmonte M. Collagen and non-collagenous proteins molecular crosstalk in the pathophysiology of osteoporosis. Cytokine Growth Factor Rev. 2019;49:59–69.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  62. He R, Lu Y, Ren J, Wang Z, Huang J, Zhu L, et al. Decreased fibrous encapsulation and enhanced osseointegration in vitro by decorin-modified titanium surface. Colloids Surf B Biointerfaces. 2017;155:17–24.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  63. Mochida Y, Parisuthiman D, Pornprasertsuk-Damrongsri S, Atsawasuwan P, Sricholpech M, Boskey AL, et al. Decorin modulates collagen matrix assembly and mineralization. Matrix Biol. 2009;28(1):44–52.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  64. Schmierer B, Hill CS. TGFbeta-SMAD signal transduction: molecular specificity and functional flexibility. Nat Rev Mol Cell Biol. 2007;8(12):970–82.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  65. Moon YJ, Yun CY, Choi H, Ka SO, Kim JR, Park BH, et al. Smad4 controls bone homeostasis through regulation of osteoblast/osteocyte viability. Exp Mol Med. 2016;48(9):e256.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  66. Hyzy SL, Olivares-Navarrete R, Schwartz Z, Boyan BD. BMP2 induces osteoblast apoptosis in a maturation state and noggin-dependent manner. J Cell Biochem. 2012;113(10):3236–45.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  67. Zheng T, Yang L, Liu Y, Liu H, Yu J, Zhang X, et al. Plasma DPP4 activities are associated with osteoporosis in postmenopausal women with normal glucose tolerance. J Clin Endocrinol Metab. 2015;100(10):3862–70.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  68. Wittrant Y, Gorin Y, Woodruff K, Horn D, Abboud HE, Mohan S, et al. High d (+) glucose concentration inhibits RANKL-induced osteoclastogenesis. Bone. 2008;42(6):1122–30.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  69. Kim HY, Mohan S. Role and mechanisms of actions of thyroid hormone on the skeletal development. Bone Res. 2013;1(2):146–61.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  70. Cabello-Verrugio C, Santander C, Cofré C, Acuña MJ, Melo F, Brandan E. The internal region leucine-rich repeat 6 of decorin interacts with low density lipoprotein receptor-related protein-1, modulates transforming growth factor (TGF)-β-dependent signaling, and inhibits TGF-β-dependent fibrotic response in skeletal muscles. J Biol Chem. 2012;287(9):6773–87.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  71. Afshinnia F, Pennathur S. Association of hypoalbuminemia with osteoporosis: analysis of the national health and nutrition examination survey. J Clin Endocrinol Metab. 2016;101(6):2468–74.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  72. Shang N, Wu J. Egg white ovotransferrin attenuates RANKL-induced osteoclastogenesis and bone resorption. Nutrients. 2019;11(9):2254. https://doi.org/10.3390/nu11092254.

    CAS  Article  PubMed Central  Google Scholar 

Download references

Acknowledgements

Technical assistance with LC–MS/MS was provided by APPLIED PROTEIN TECHNOLOGY (Shanghai) Co., Ltd.

Funding

This study was supported by the Science and Technology Traditional Chinese Medicine Plan of Zhejiang Province, China (No. 2021ZB149); Major Disease Prevention and Control of a Public Relations Program of Chinese Medicine of Zhejiang Province, China (No. 2012ZCG004); Science and Technology Traditional Chinese Medicine Plan of Zhejiang Province, China (No. 2016ZA105); and the support of the Native Chinese medicine studio of Dr. Xinmiao Yao.

Author information

Authors and Affiliations

Authors

Contributions

RZ and KY made substantial contributions to the concept and design of the study and wrote the manuscript. YW performed most of the experiments. XY, GL and LG assisted with research and protocol design. ZC helped with the data analysis and revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Zhineng Chen.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the laboratory animal management and ethics committee of Zhejiang Chinese Medical University (Ethics Application No. 2015018).

Consent for publication

All authors consent to the publication of this manuscript.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Ruikun Zhang and Kun Yan are CO-author

Supplementary Information

Additional file 1. 

Identification versus quantification results statistical table.

Additional file 2. 

Table of differentially expressed protein statistics.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhang, R., Yan, K., Wu, Y. et al. Quantitative proteomics reveals the effect of Yigu decoction (YGD) on protein expression in bone tissue. Clin Proteom 18, 24 (2021). https://doi.org/10.1186/s12014-021-09330-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12014-021-09330-0

Keywords

  • Yigu decoction
  • Osteoporosis
  • TMT proteomics
  • Bone metabolism