- Open Access
Transcriptomics‐based analysis of the mechanism by which Wang-Bi capsule alleviates joint destruction in rats with collagen‐induced arthritis
Chinese Medicine volume 16, Article number: 31 (2021)
Rheumatoid arthritis (RA) is a chronic autoimmune disease accompanied with joint destruction that often leads to disability. Wang-Bi capsule (WB), a traditional Chinese medicine-based herbs formula, has exhibited inhibition effect on joint destruction of collagen-induced arthritis (CIA) animal model in our previous study. But its molecular mechanisms are still obscure.
CIA rats were treated intragastrical with WB for eight weeks, and the effect of joints protection were evaluated by hematoxylin and eosin (H&E) staining, safranin O fast green staining, tartrate-resistant acid phosphatase (TRAP) staining and micro‑CT scanning analysis. The transcriptomic of tarsal joints were used to investigate how WB alleviated joint destruction.
The histological examination of ankle joints showed WB alleviated both cartilage damage and bone destruction of CIA rats. This protective effect on joints were further evidenced by micro-CT analysis. The transcriptomic analysis showed that WB prominently changed 12 KEGG signaling pathways (“calcium signaling pathway”, “cAMP signaling pathway”, “cell adhesion molecules”, “chemokine signaling pathway”, “complement and coagulation cascades”, “MAPK signaling pathway”, “NF-kappa B signaling pathway”, “osteoclast differentiation”, “PI3K-Akt signaling pathway”, “focal adhesion”, “Gap junction” and “Rap1 signaling pathway”) associated with bone or cartilage. Several genes (including Il6, Tnfsf11, Ffar2, Plg, Tnfrsf11b, Fgf4, Fpr1, Siglec1, Vegfd, Cldn1, Cxcl13, Chad, Arrb2, Fgf9, Egfr) regulating bone resorption, bone formation and cartilage development were identified by further analysis. Meanwhile, these differentially expressed genes were validated by real-time quantitative PCR.
Overall, the protective effect of WB treatment on joint were confirmed in CIA rats, and its basic molecular mechanisms may be associated with regulating some genes (including Il6, Tnfsf11, Ffar2 and Plg etc.) involved in bone resorption, bone formation and cartilage development.
Rheumatoid arthritis (RA) is a chronic, autoimmune disease that is accompanied by persistent inflammation of synovial tissue as well as joint destruction. Joint destruction often interferes with physical function and even results in disability, which further affects productivity and quality of life. Approximately 80% of patients have deformed joints, and 40% lose work due to disability within 10 years of disease onset . Many earlier studies indicated that joint destruction was correlated with signs and symptoms of inflammation, and some patients with RA have benefited from anti-inflammatory therapies. However, joint destruction still occurs with a high incidence after anti-inflammatory treatment. Accumulating evidence in recent years has revealed that inflammatory activity and joint destruction may be dissociated from one another [2, 3]. Thus, anti-inflammatory drugs are not sufficient to protect joints, and agents that can directly inhibit joint destruction are essential for RA patients. Although several biologic agents have exhibited articular protection to some extent, how to effectively block joint destruction and protect joint function in RA is still a great challenge for clinical doctors because of the lack of sufficient joint protection agents to date [4,5,6].
Wang-Bi capsule (WB), approved by the China Food and Drug Administration (approval ID of CFDA: Z20080096), is a traditional Chinese medicine-based herbal formula. It has been used for RA treatment in China for many years. WB displayed joint protective effects in a RA mouse model, and it could regulate osteoclast-osteoblast functions and inhibit bone destruction in joints in our previous study . Icariin and paeoniflorin are both active ingredients of WB. Icariin abrogated osteoclast differentiation in vitro and inhibited bone loss by regulating the coupling process of osteogenic and osteoclastic activity in rats with experimental osteoporosis [8,9,10,11]. Paeoniflorin suppresses osteoclast differentiation by controlling the receptor activator of nuclear factor-κ B ligand (RANKL)/osteoprotegerin (OPG) ratio . Furthermore, icariin and paeoniflorin were also proven to regulate the activity of chondrocytes and inhibit articular cartilage degradation [13, 14]. These studies hinted that WB was a promising candidate drug for joint protection in RA patients. However, it is still unclear how WB inhibits joint destruction in RA. In this study, we aimed to demonstrate the joint protective effects of WB in a CIA rat model and explore the molecular mechanism of action via transcriptomic analysis.
Male 6- to 8-week-old Sprague Dawley rats (Vital River Laboratory Animal Technology Co. Ltd., Beijing, China) were used for this study, and all rats were housed in a specific pathogen-free facility with free access to sterilized water and chow. In addition, the experiments were reviewed and approved by the Research Ethics Committee of the Institute of Basic Theory of Chinese Medicine, China Academy of Chinese Medical Sciences.
Induction of the CIA model
The CIA model was established in line with the method reported in our previous study . The rats were given an intradermal injection of 200 µl of an emulsified mixture of equal proportions bovine collagen type II and incomplete Freund’s adjuvant (Chondrex, Redmond, WA, USA) on day 1. In addition, the rats were immunized on day 7 by another injection of 100 µl of the emulsified mixture.
WB was provided by Liaoning China Resources Benxi Sanyao Co., Ltd. (Liaoning, China, No. 20180205). The rats were randomly divided into three groups after successful induction of the CIA model, with eight rats in each group: the model group, WB group (WB 0.74 g/kg/d, equal to that used in RA patients) and normal group. The rats in the WB group were intragastrically administered WB for 8 weeks after successful induction of the CIA model. The rats in the normal and model groups were administered an equal volume of double-distilled water in the same way. All animals were sacrificed at the end of the experiment. The sera were harvested for enzyme-linked immunosorbent assay (ELISA). The hind limbs were collected for micro-CT analysis, and the ankle joints were removed for histological analysis. The tarsal joints were separated from the hind limbs and used for transcriptomics analysis and real-time quantitative PCR (RT-qPCR).
The serum levels of type I collagen N-terminal propeptide (PINP) and C-terminal telopeptide of type I collagen (CTX-1) were detected using commercially available ELISA kits (Sino-uk bio, Beijing, China). The procedures were performed according to the manufacturer’s instructions.
The ankle joints of rats were fixed in 10% phosphate-buffered formalin for 3 days and then decalcified in 10% EDTA (pH = 7.2). Section (4 μm) were stained with hematoxylin and eosin (H&E) for histopathological assessment, and the histopathological characteristics were scored on a scale of 0–4 using the criteria described in a previous study . In addition, sections were stained with safranin O fast green (Solarbio, Beijing, China) for cartilage observation, and cartilage damage was evaluated with the following criteria: 0 = normal, 1 = mild loss of safranin O staining with no obvious chondrocyte loss, 2 = moderate loss of staining with focal mild chondrocyte loss, 3 = marked loss of staining with multifocal marked chondrocyte loss, and 4 = severe diffuse loss of staining with multifocal severe chondrocyte loss . TRAP staining was performed using a TRAP staining kit (Sigma, St. Louis, MO, USA) to visualize osteoclast formation, and the number of osteoclasts per bone surface around the taluses was determined .
Micro‑CT scanning analysis
The hind paws and femurs were scanned for reconstruction into a three-dimensional (3-D) structure with a SKYSCAN 1174 micro-CT (Bruker, Belgium). The hind paws were analyzed to test the bone volume (BV)/tissue volume (TV) ratio, and the femurs were analyzed to determine the bone mineral density (BMD) of the periarticular bone, trabecular number (Tb.N), trabecular pattern factor (Tb.Pf) and trabecular separation (Tb.S) .
Transcriptome of the tarsal joints
Three tarsal joints were sampled per group (normal group, model group and WB group) at the end of the experiment. The skin, muscle, tendon and ligament were removed, but bone tissue was included in each sample . The samples were stored at − 80 °C. Total RNA was extracted from the tissues using TRIzol (Invitrogen, Carlsbad, CA, USA) according to the instructions in the product manual. Subsequently, total RNA was qualified and quantified using a NanoDrop (NanoDrop, Madison, USA) and an Agilent 2100 bioanalyzer (Agilent, Santa Clara, USA). The products were enriched with PCR to create the final cDNA library for mRNA sequencing. The constructed library was amplified to make DNA nanoballs (DNBs) with more than 200 copies of DNA per molecule. The DNBs were loaded into the patterned nanoarray, and single-end 100-base reads were generated by sequence-based synthesis. The qualified libraries were subjected to paired-end sequencing on the BGISEQ-500 System (BGI-Shenzhen, China).
The sequencing data were filtered with SOAPnuke (v1.5.2), and clean reads were obtained and stored in FASTQ format. The clean reads were mapped to the reference genome using HISAT2 (v2.0.4). Bowtie2 (v2.2.5) was applied to align the clean reads to the reference coding gene set, and the expression levels of genes were calculated with RSEM (v1.2.12).
Differential expression analysis was performed using DESeq2 (v1.4.5) with a fold change > 2.0 and a P-value < 0.01. To gain insight into the change in phenotype, KEGG (https://www.kegg.jp/) enrichment analysis of annotated differentially expressed genes was performed with Phyper (https://en.wikipedia.org/wiki/Hypergeometric_distribution) based on a hypergeometric test. The significance levels of pathways were filtered by P-value, with a rigorous threshold (P-value < 0.01) in the Bonferroni test. Network analysis was performed by ingenuity pathway analysis (IPA, Ingenuity Systems).
RNA (2 µg per reaction) was reverse transcribed to cDNA using TRUEscript RT MasterMix (Aibosen Biotechnologies Co., Ltd., Beijing, China). qPCR was performed with 2× SYBR Green qPCR Master Mix (Biomake, Houston, TX, USA). PCR was initiated at 95 °C for 5 min, and the cycling conditions were 95 °C for 10 s, 60 °C for 15 s and 75 °C 15 s (40 cycles). Samples were analyzed in triplicate and normalized by subtraction of the sample mean Ct value from the mean Ct value of the housekeeping gene GAPDH. The primer information is shown in Additional file 1: Table S1.
GraphPad Prism v.6 (GraphPad Software, San Diego, CA, USA) was used for statistical analyses. Ranked data including histological scoring were assessed with the Kruskal-Wallis test followed by a multiple comparisons test. Measurement data including RT-qPCR were analyzed by ANOVA, followed by Bonferroni’s multiple comparison test. P < 0.05 was identified significant.
WB inhibited joint destruction in rats with CIA
The joints of rats in the model group displayed severe swelling and deformation, and WB treatment remarkably improved these symptoms compared with those of model group rats (Fig. 1a). To further evaluate the effect of WB on joint destruction in rats with CIA, H&E staining was performed. The results showed no detectable pathological changes in the rats in the normal group. In contrast, the rats in the model group exhibited severe inflammatory cell infiltration and bone and cartilage damage, whereas these pathological changes in the WB group were significantly relieved compared with those in the model group (Fig. 1b, e). Next, safranin O fast green staining was used to further investigate the changes in articular cartilage. The results showed severe cartilage damage in the model group compared with the normal group, represented by a markedly decreased amount of cartilage matrix. However, the cartilage damage in the WB group exhibited obvious mitigation compared with that of the model group (Fig. 1c, f). Finally, TRAP staining was performed to detect the number of osteoclasts. We did not detect positive staining for osteoclasts in the normal group, but a large number of osteoclasts were observed in the model group. In contrast, the number of osteoclasts with positive staining in the WB group was significantly decreased compared with that in the model group (Fig. 1d, g). These results indicated that WB could effectively relieve bone and cartilage destruction in rats with CIA.
WB increased the PINP/CTX-1 ratio of rats with CIA
PINP and CTX-1 are important biochemical markers of bone formation and bone resorption, respectively. The serum levels of PINP and CTX-1 were all significantly increased in the model group compared with the normal group (Fig. 2a, b). However, the PINP/CTX-1 ratio was higher in the WB group than in the model group (Fig. 2c), which provided evidence that WB could promote bone formation in rats with CIA.
WB inhibited bone erosion and bone loss in rats with CIA
Micro-CT and 3D image analysis were used to examine the effect of WB on bone erosion and bone microstructure. The CT analysis showed that bone erosion was substantially more severe in the model group than in the normal group. After WB treatment, bone erosion was inhibited (Fig. 3a). Furthermore, the BV/TV ratio in the WB group was significantly increased compared with that in the model group (Fig. 3b). In addition, the distal femur microstructure was reconstituted (Fig. 3c). We found that BMD and Tb.N in the model group were obviously lower than those in the normal group, whereas these two parameters in the WB group were significantly increased compared with those in the model group (Fig. 3d, e). Both the Tb.S and Tb.Pf of the model group increased compared to those of the normal group, and WB treatment significantly decreased the Tb.S and Tb.Pf (Fig. 3f, g). Overall, these results indicated that WB could inhibit bone erosion and bone loss in rats with CIA.
WB regulated signaling pathways associated with bone or cartilage
Although WB exhibited a protective effect on joints, the potential mechanism was still unclear. Therefore, we wanted to reveal the mechanism via transcriptomic analysis. A total of 3790 differentially expressed genes (DEGs) were identified in Mod-VS-Nor; 2988 DEGs showed downregulated expression, and 802 DEGs showed upregulated expression. A total of 4139 DEGs were identified in WB-VS-Mod; 3362 DEGs showed upregulated expression, and 777 DEGs showed downregulated expression (Fig. 4a). By KEGG analysis of DEGs in Mod-VS-Nor and WB-VS-Mod, we found 12 obviously enriched (P value < 0.05) signaling pathways associated with bone or cartilage, including “calcium signaling pathway”, “cell adhesion molecules (CAMs)”, “complement and coagulation cascades”, “osteoclast differentiation”, “PI3K-Akt signaling pathway”, “cAMP signaling pathway”, “chemokine signaling pathway”, “NF-kappa B signaling pathway”, “MAPK signaling pathway”, “focal adhesion”, “Rap1 signaling pathway” and “gap junction” (Fig. 4b).
WB regulated the levels of some genes that controlled bone and cartilage function in rats with CIA
Having identified the 12 signaling pathways that regulated bone and cartilage function, we next focused on the DEGs belonging to these signaling pathways. A total of 373 DEGs from Mod-VS-Nor and 425 DEGs from WB-VS-Mod belonged to the 12 pathways associated with bone or cartilage identified by KEGG analysis, and there were 238 overlapping DEGs between these two groups. These genes were input into IPA, and three functional networks related to bone or cartilage were generated: the bone resorption network, bone formation network and cartilage development network (Fig. 5). Moreover, the levels of genes promoting bone resorption (including Il6, Tnfsf11, Cxcl1, Acp5 and Mef2c) were elevated in the model group but were suppressed by WB treatment, while the levels of genes inhibiting bone resorption (including Ffar2, Plg, Tnfrsf11b, Hcn4, Adrb1, Prkaa2 and Gipr) were decreased in the model group but were increased by WB treatment (Additional file 1: Table S2). The transcript levels of genes promoting bone formation (including Fgf4, Fpr1, Siglec1, Vegfd, Cldn1, Taok3, Grin3a, Tshr, Plcb1 and Ntf3) were decreased in the model group but were increased by WB treatment (Additional file 1: Table S3). Genes involved with cartilage development (including Cxcl13, Chad, Arrb2, Fgf9, Egfr and Clu) exhibited low levels in the model group but higher levels after WB treatment (Additional file 1: Table S4).
DEG validation by RT-qPCR
The top 5 DEGs (according to fold change) associated with bone resorption, bone formation and cartilage development were chosen for validation by RT-qPCR. The mRNA expression levels of genes promoting bone formation (Fgf4, Fpr1, Siglec1, Vegfd, and Cldn1) were remarkably decreased in the model group compared with the normal group, but they were all significantly increased in the WB group compared with the model group (Fig. 6a). The levels of Il6 and Tnfsf11 were much higher in the model group than in the normal group; however, they were significantly downregulated by WB treatment. In contrast, the levels of genes suppressing osteoclast differentiation, such as Ffar2, Plg, and Tnfrsf11b, were markedly elevated by WB treatment (Fig. 6b). Furthermore, the transcript levels of genes controlling cartilage development, such as Cxcl13, Chad, Arrb2, Fgf9 and Egfr, were decreased in the model group, whereas they were markedly upregulated in the WB group (Fig. 6c). The changes in the transcript levels of these genes were consistent with the tendencies observed in the transcriptome analysis.
Joint destruction is one of the predominant types of damage underlying RA disease. Accumulating evidence has proven a positive correlation between joint destruction and disability, therefore, it is essential to consider whether therapeutic agents have a direct impact on progressive joint damage. Our previous studies revealed that WB not only inhibited the inflammatory response but also alleviated joint destruction . In this study, we mainly focused on WB’s protective effects on joints. WB treatment obviously improved the histopathological manifestations of CIA in rats, including inhibition of bone destruction and cartilage damage. Bone loss often occurs in RA, which increases fracture risk. Rats with CIA displayed high bone turnover and osteoporosis in previous studies . In our study, the rats with CIA also exhibited high bone turnover and presented elevated serum levels of both PINP (a bone formation marker) and CTX-1 (a bone resorption marker). WB treatment increased the PINP/CTX-1 ratio, which indicated that WB could inhibit bone loss in rats with CIA. Furthermore, the protective effect of WB treatment on bone was also evidenced by the micro-CT examination results.
The mechanisms by which WB protects joints are complicated because of the intricate ingredients of WB. Accordingly, transcriptomic analysis of the tarsal joints was performed to evaluate the impacts of WB on cartilage and bone. The KEGG analysis of DEGs identified several signaling pathways associated with bone and cartilage. We know that the balance of osteoclasts and osteoblasts is important for bone function, and “osteoclast differentiation” is a core factor that controls osteoclast formation . The KEGG database shows that pathways including the “calcium signaling pathway”, “MAPK signaling pathway”, “PI3K-Akt signaling pathway” and “NF-kappa B signaling pathway” participate in the induction of “osteoclast differentiation”. Some studies have also revealed that these pathways are deeply involved in regulating osteoclast differentiation or activation, which further impacts bone resorption [22,23,24,25]. “Focal adhesion” and the “MAPK signaling pathway” are important parts of the “Wnt signaling pathway”, which plays an essential role in osteoblast differentiation. Furthermore, osteochondrogenitor cells express both Runx2 and Sox9, and “cAMP signaling” modifies these transcription factors and drives differentiation to osteoblasts or chondrocytes [26, 27]. In addition, cAMP induces RAP1 activation, which promotes integrin clustering and cell adhesion to the bone matrix [28, 29]. The migration and maturation of osteoclasts and osteoblasts are important processes in bone remodeling. Pathways including “cell adhesion molecules”, “chemokine signaling pathway” and “gap junctions” all help to facilitate these processes [30,31,32,33,34]. “Complement and coagulation cascades” are seldom reported in the context of regulating bone remodeling in RA. However, several studies have hinted at a role for these pathways. For example, FVIIIKO and PAR1KO mice displayed a lower bone/tissue volume ratio and a smaller number of bone trabeculae than wild-type mice [35, 36]. Overall, we found that WB could regulate these signaling pathways, which helped us understand how WB influenced the joints of rats with CIA.
To further dissect the mechanism by which WB protects bone and cartilage, we analyzed these DEGs (belonging to the signaling pathways mentioned above) from several profiles. Several DEGs (Il6, Tnfsf11, Cxcl1, Acp5, and Mef2c) promoting bone resorption were first identified. Their expression levels were all decreased by WB treatment. IL-6 is a proinflammatory cytokine, and multiple studies in recent years have shown that it is also an important factor for promoting osteoclast differentiation. The inhibition of IL-6 retarded the progression of joint damage in RA independent of its anti-inflammatory effects, which suggested disassociation of the link between inflammation and joint destruction [37, 38]. RANKL (encoded by Tnfsf11) is absolutely necessary for osteoclast differentiation, and both Cxcl1 and Mef2c can enhance RANKL-induced osteoclast differentiation [39,40,41]. Furthermore, Acp5, encoding TRAP (a marker of osteoclasts), was also suppressed by WB treatment. In addition, the mRNA levels of genes (Ffar2, Plg, Tnfrsf11b, Hcn4, Adrb1, Prkaa2 and Gipr) that suppressed bone resorption were elevated by WB treatment. Plasminogen (Plg) is an important regulator of bone metabolism that has not been described in RA. Exogenous Plg suppressed RAW264.7 cell differentiation to osteoclasts. Moreover, Plg-/- mice displayed a low BMD . OPG (encoded by Tnfrsf11b), a soluble decoy receptor for RANKL, can inhibit bone loss by controlling osteoclast differentiation. AMPKα2 (encoded by Prkaa2) could inhibit RANKL signaling . Ffar2 and Hcn4 are both seldom reported in RA studies, but gene knockdown experiments have shown decreased formation of osteoblasts or BMD [44, 45]. Stimulation of the sympathetic nervous system by beta-2 adrenergic receptor (β2AR) plays an important role in mediating bone remodeling . β1AR (encoded by Adrb1) has also been proven to mitigate negative changes in cancellous bone microarchitecture and increase the bone mass of mice . Glucose-dependent insulinotropic peptide (GIP), an intestinally secreted hormone, could decrease bone resorption. Excess signaling during GIPR (encoded by Gipr) activation resulted in increased bone mass . In this study, we found that rats with CIA displayed low transcript levels of Ffar2, Plg, Tnfrsf11b, Hcn4, Adrb1, Prkaa2 and Gipr compared with those of normal rats, whereas the levels of all these genes transcripts were significantly increased by WB treatment, which demonstrated that WB treatment could effectively inhibit bone resorption in rats with CIA.
In addition to influencing genes regulating bone resorption, WB also regulated several genes regulating bone formation. Bone morphogenetic protein family members (BMPs) are important signaling molecules regulating bone formation. Fgf4, with a very high transcript level after WB treatment, could cooperate with BMP-2-induced bone formation, and administration of FGF4 could increase bone mass in rats [49, 50]. Similar to FGF4, PLCβ1 can also cooperate with BMPs. Increased expression of PLCβ1 induced by BMP-2 could promote osteoblast differentiation in vitro . Activation of the N-formyl peptide receptor (FPR) can trigger multiple biochemical cascades and eventually lead to cellular activation. Its role during the differentiation of osteocytes has received attention in recent years. Activation of FPR1 could promote osteogenesis and mineralization, but activation of FPR2 or FPR3 could not . The anabolic process is an important process of bone remodeling. For example, CD169 (encoded by Siglec1) could provide pro-anabolic support during the process of bone repair . VEGF-D could promote the formation of mineralized nodules and accelerate bone formation, which provides another option to promote bone formation . Claudins (Cldns) are major components of tight junctions. In addition to functioning as tight junctions, they have been found to play roles in cell proliferation and differentiation. Cldn1 was downregulated during osteoclast differentiation but upregulated during osteoblast differentiation, and Cldn1 knockdown decreased the expression of Runx2 . Proteins encoded by genes such as Taok3, Tshr, Grin3a and Ntf3 could also promote bone formation by enhancing osteoblast differentiation or mineralization [56,57,58,59,60]. The transcript levels of these genes were elevated by WB treatment, which might contribute to bone formation through multiple pathways.
Articular cartilage damage is a very important pathological change during joint destruction in RA. Cartilage tissue is composed of extracellular matrix (ECM) that is synthesized by chondrocytes. The degradation of cartilage ECM also provides an opportunity for the hyperplastic synovium to directly invade subchondral bone. Although cartilage tissue has limited self-repair ability, recent studies have shown that cartilage contains a population of stem cells or progenitor cells that facilitate cartilage development and repair . Thus, cartilage may benefit from WB treatment through several genes, such as Cxcl13, Chad, Arrb2, Fgf9, Egfr and Clu, that regulate cartilage development. Cartilage homeostasis is disrupted when damage occurs. Chemokines and chondroadherin (CHAD) could contribute to the maintenance of cartilage homeostasis by promoting ECM production or assembly [57, 61]. FGF9 is a member of the fibroblast growth factor (FGF) family, and studies have shown that it could enhance chondrogenesis in dental pulp stem cells . However, the role of FGF9 with respect to the articular cartilage damage underlying RA has not been well described. The epidermal growth factor receptor (EGFR) system plays important roles in multiple processes, such as the differentiation and proliferation of osteoclasts, osteoblasts and chondrocytes . Both beta-arrestin-2 (ARRB2) and clusterin (CLU) could inhibit the inflammatory response of chondrocytes, but CLU could also influence the proliferation and differentiation of chondrocytes [64, 65]. In this study, we found that the expression of Cxcl13, Chad, Arrb2, Fgf9, Egfr and Clu was obviously suppressed in rats with CIA, but WB treatment increased it to a certain degree.
Overall, the protective effect of WB on joints was confirmed in rats with CIA, and its basic molecular mechanisms may be associated with regulation of the expression levels of some genes (including Il6, Tnfsf11, Ffar2, Plg, Tnfrsf11b, Fgf4, Fpr1, Siglec1, Vegfd, Cldn1, Cxcl13, Chad, Arrb2, Fgf9, and Egfr) involved in bone resorption, bone formation and cartilage development.
Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.
Receptor activator of nuclear factor-κ B ligand
Type I collagen N-terminal propeptide
C-terminal telopeptide of type I collagen
Bone volume/tissue volume
Bone mineral density
Trabecular pattern factor
Differentially expressed genes
Aletaha D, Smolen JS. Diagnosis and management of rheumatoid arthritis: a review. Electronic. 2018;25:1538–3598.
Drossaers-Bakker KW, et al. Long-term course and outcome of functional capacity in rheumatoid arthritis: the effect of disease activity and radiologic damage over time. Arthritis Rheum. 1999;42(9):1854–60.
Scott DL, et al. The links between joint damage and disability in rheumatoid arthritis. Rheumatology. 2000;39(2):122–32.
Binder NB, et al. Tumor necrosis factor-inhibiting therapy preferentially targets bone destruction but not synovial inflammation in a tumor necrosis factor-driven model of rheumatoid arthritis. Arthritis Rheum. 2013;65(3):608–17.
Smolen JS, et al. Rheumatoid arthritis. Nat Rev Dis Primers. 2018;4:18001.
Smolen JS, et al. Evidence of radiographic benefit of treatment with infliximab plus methotrexate in rheumatoid arthritis patients who had no clinical improvement: a detailed subanalysis of data from the anti-tumor necrosis factor trial in rheumatoid arthritis with concomitant therapy study. Arthritis Rheum. 2005;52(4):1020–30.
Cui H, et al. Wang-Bi Capsule Alleviates the Joint Inflammation and Bone Destruction in Mice with Collagen-Induced Arthritis. Evid Based Complement Alternat Med. 2020;2020:1015083.
Jing X, et al. Icariin protects against iron overload-induced bone loss via suppressing oxidative stress. J Cell Physiol. 2019;234(7):10123–37.
Kim B, Lee KY, Park B. Icariin abrogates osteoclast formation through the regulation of the RANKL-mediated TRAF6/NF-kappaB/ERK signaling pathway in Raw264.7 cells. Phytomedicine. 2018;51:181–90.
Xie Y, et al. Icariin-loaded porous scaffolds for bone regeneration through the regulation of the coupling process of osteogenesis and osteoclastic activity. Int J Nanomedicine. 2019;14:6019–33.
Zhang S, et al. Icariin influences adipogenic differentiation of stem cells affected by osteoblast-osteoclast co-culture and clinical research adipogenic. Biomed Pharmacother. 2017;88:436–42.
Xu H, et al. Paeoniflorin ameliorates collagen-induced arthritis via suppressing nuclear factor-kappaB signalling pathway in osteoclast differentiation. Immunology. 2018;12:58.
Zuo S, et al. Icariin alleviates IL-1beta-induced matrix degradation by activating the Nrf2/ARE pathway in human chondrocytes. Drug Des Devel Ther. 2019;13:3949–61.
Hu PF, et al. Paeoniflorin inhibits IL-1beta-induced MMP secretion via the NF-kappaB pathway in chondrocytes. Exp Ther Med. 2018;16(2):1513–9.
Guo Q, et al. Wu-Tou decoction in rheumatoid arthritis: integrating network pharmacology and in vivo pharmacological evaluation. Front Pharmacol. 2017;8:230.
Mukai T, et al. Loss of SH3 domain-binding protein 2 function suppresses bone destruction in tumor necrosis factor-driven and collagen-induced arthritis in mice. Arthritis Rheumatol. 2015;67(3):656–67.
Akagi T, et al., Effect of angiotensin II on bone erosion and systemic bone loss in mice with tumor necrosis factor-mediated arthritis. Int J Mol Sci. 2020;21:11.
Zhao H, et al. Yi Shen Juan Bi pill ameliorates bone loss and destruction induced by arthritis through modulating the balance of cytokines released by different subpopulations of T cells. Front Pharmacol. 2018;9:262–2.
Denninger KC, et al. Kinetics of gene expression and bone remodelling in the clinical phase of collagen-induced arthritis. Arthritis Res Ther. 2015;17:43.
Zhang X, et al. Effect of prednisone treatment for 30 and 90 days on bone metabolism in collagen-induced arthritis (CIA) rats. J Bone Miner Metab. 2018;36(6):628–39.
Cho E, et al. PSTP-3,5-Me inhibits osteoclast differentiation and bone resorption. Molecules. 2019;24:18.
Lee K, et al. Roles of mitogen-activated protein kinases in osteoclast biology. Int J Mol Sci. 2018;19:10.
Mishra R, et al. NF-kappaB signaling negatively regulates osteoblast dedifferentiation during zebrafish bone regeneration. Dev Cell. 2020;52(2):167-82 e7.
Ding N, et al. Physalin D inhibits RANKL-induced osteoclastogenesis and bone loss via regulating calcium signaling. BMB Rep. 2020;53(3):154–9.
Chen LL, et al. PI3K/AKT pathway involvement in the osteogenic effects of osteoclast culture supernatants on preosteoblast cells. Tissue Eng Part A. 2013;19(19–20):2226–32.
Topol L, et al. Sox9 inhibits Wnt signaling by promoting beta-catenin phosphorylation in the nucleus. J Biol Chem. 2009;284(5):3323–33.
Carroll SH, Ravid K. Differentiation of mesenchymal stem cells to osteoblasts and chondrocytes: a focus on adenosine receptors. Expert Rev Mol Med. 2013;15:e1.
Jeevaratnam K, et al. Regulatory actions of 3’,5’-cyclic adenosine monophosphate on osteoclast function: possible roles of Epac-mediated signaling. Ann N Y Acad Sci. 2018;1433(1):18–28.
Su X, et al. The ADP receptor P2RY12 regulates osteoclast function and pathologic bone remodeling. J Clin Invest. 2012;122(10):3579–92.
Kim H, et al. CD44 can compensate for IgSF11 deficiency by associating with the scaffold protein PSD-95 during osteoclast differentiation. Int J Mol Sci. 2020;21:7.
Xuan W, et al. Osteoclast differentiation gene expression profiling reveals chemokine CCL4 mediates RANKL-induced osteoclast migration and invasion via PI3K pathway. Cell Biochem Funct. 2017;35(3):171–7.
Mathieu PS, Loboa EG. Cytoskeletal and focal adhesion influences on mesenchymal stem cell shape, mechanical properties, and differentiation down osteogenic, adipogenic, and chondrogenic pathways. Tissue Eng Part B Rev. 2012;18(6):436–44.
Matemba SF, Lie A, Ransjo M. Regulation of osteoclastogenesis by gap junction communication. J Cell Biochem. 2006;99(2):528–37.
Ilvesaro J, Tuukkanen J. Gap-junctional regulation of osteoclast function. Crit Rev Eukaryot Gene Expr. 2003;13(2–4):133–46.
Aronovich A, et al. A novel role for factor VIII and thrombin/PAR1 in regulating hematopoiesis and its interplay with the bone structure. Blood. 2013;122(15):2562–71.
Nguyen TS, Lapidot T, Ruf W. Extravascular coagulation in hematopoietic stem and progenitor cell regulation. Blood. 2018;132(2):123–31.
Smolen JS, Avila JC, Aletaha D. Tocilizumab inhibits progression of joint damage in rheumatoid arthritis irrespective of its anti-inflammatory effects: disassociation of the link between inflammation and destruction. Ann Rheum Dis. 2012;71(5):687–93.
Sims NA. Cell-specific paracrine actions of IL-6 family cytokines from bone, marrow and muscle that control bone formation and resorption. Int J Biochem Cell Biol. 2016;79:14–23.
Boyce BF, Xing L. Functions of RANKL/RANK/OPG in bone modeling and remodeling. Arch Biochem Biophys. 2008;473(2):139–46.
Valerio MS, et al. Critical role of MKP-1 in lipopolysaccharide-induced osteoclast formation through CXCL1 and CXCL2. Cytokine. 2015;71(1):71–80.
Kramer I, et al. Mef2c deletion in osteocytes results in increased bone mass. J Bone Miner Res. 2012;27(2):360–73.
Kanno Y, et al. Plasminogen/plasmin modulates bone metabolism by regulating the osteoblast and osteoclast function. J Biol Chem. 2011;286(11):8952–60.
Kang H, Viollet B, Wu D. Genetic deletion of catalytic subunits of AMP-activated protein kinase increases osteoclasts and reduces bone mass in young adult mice. J Biol Chem. 2013;288(17):12187–96.
Montalvany-Antonucci CC, et al. Short-chain fatty acids and FFAR2 as suppressors of bone resorption. Bone. 2019;125:112–21.
Notomi T, et al. Zinc-Induced Effects on Osteoclastogenesis Involves Activation of Hyperpolarization-Activated Cyclic Nucleotide Modulated Channels via Changes in Membrane Potential. J Bone Miner Res. 2015;30(9):1618–26.
Bonnet N, et al. Alteration of trabecular bone under chronic beta2 agonists treatment. Med Sci Sports Exerc. 2005;37(9):1493–501.
Swift JM, et al. Beta-1 adrenergic agonist treatment mitigates negative changes in cancellous bone microarchitecture and inhibits osteocyte apoptosis during disuse. PLoS One. 2014;9(9):e106904.
Xie D, et al. Glucose-dependent insulinotropic peptide-overexpressing transgenic mice have increased bone mass. Bone. 2007;40(5):1352–60.
Kubota K, et al. Synergistic effect of fibroblast growth factor-4 in ectopic bone formation induced by bone morphogenetic protein-2. Bone. 2002;31(4):465–71.
Kuroda S, et al. Anabolic effect of aminoterminally truncated fibroblast growth factor 4 (FGF4) on bone. Bone. 1999;25(4):431–7.
Ramazzotti G, et al. BMP-2 induced expression of PLCbeta1 that is a positive regulator of osteoblast differentiation. J Cell Physiol. 2016;231(3):623–9.
Shin MK, et al. N-formyl-methionyl-leucyl-phenylalanine (fMLP) promotes osteoblast differentiation via the N-formyl peptide receptor 1-mediated signaling pathway in human mesenchymal stem cells from bone marrow. J Biol Chem. 2011;286(19):17133–43.
Batoon L, et al. CD169(+) macrophages are critical for osteoblast maintenance and promote intramembranous and endochondral ossification during bone repair. Biomaterials. 2019;196:51–66.
Orlandini M, et al. Vascular endothelial growth factor-D activates VEGFR-3 expressed in osteoblasts inducing their differentiation. J Biol Chem. 2006;281(26):17961–7.
Alshbool FZ, Mohan S. Differential expression of claudin family members during osteoblast and osteoclast differentiation: Cldn-1 is a novel positive regulator of osteoblastogenesis. PLoS One. 2014;9(12):e114357.
Su YW, et al. Roles of neurotrophins in skeletal tissue formation and healing. J Cell Physiol. 2018;233(3):2133–45.
Su YW, et al. Neurotrophin-3 induces BMP-2 and VEGF activities and promotes the bony repair of injured growth plate cartilage and bone in rats. J Bone Miner Res. 2016;31(6):1258–74.
Boutin A, et al., Beta-Arrestin-1 mediates thyrotropin-enhanced osteoblast differentiation. FASEB J. 2014;28(8):3446-55.
Li JL, et al. NMDA enhances stretching-induced differentiation of osteoblasts through the ERK1/2 signaling pathway. Bone. 2008;43(3):469–75.
Li Z, et al. TAOK3 is a MAP3K contributing to osteoblast differentiation and skeletal mineralization. Biochem Biophys Res Commun. 2020;531(4):497–502.
Batista MA, et al. Nanomechanical phenotype of chondroadherin-null murine articular cartilage. Matrix Biol. 2014;38:84–90.
Dai J, et al. The effect of co-culturing costal chondrocytes and dental pulp stem cells combined with exogenous FGF9 protein on chondrogenesis and ossification in engineered cartilage. Biomaterials. 2012;33(31):7699–711.
Schneider MR, Sibilia M, Erben RG. The EGFR network in bone biology and pathology. Trends Endocrinol Metab. 2009;20(10):517–24.
Campo GM, et al. Beta-arrestin-2 negatively modulates inflammation response in mouse chondrocytes induced by 4-mer hyaluronan oligosaccharide. Mol Cell Biochem. 2015;399(1–2):201–8.
Tarquini C, et al. Clusterin exerts a cytoprotective and antioxidant effect in human osteoarthritic cartilage. Aging. 2020;12(11):10129–46.
This research is supported by the National Key R&D Program of China (2018YFC1705205).
Ethics approval and consent to participate
This study was reviewed and approved by the Research Ethics Committee of Institute of Basic Theory of Chinese Medicine, China Academy of Chinese Medical Sciences.
Consent for publication
The authors declare that there are no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Primers used for RT-qPCR. Table S2. Fold change of DEGs related to bone resorption regulated by WB. Table S3. Fold change of DEGs related to bone formation regulated by WB. Table S4. Fold change of DEGs related to cartilage development regulated by WB.
About this article
Cite this article
Shu, H., Zhao, H., Shi, Y. et al. Transcriptomics‐based analysis of the mechanism by which Wang-Bi capsule alleviates joint destruction in rats with collagen‐induced arthritis. Chin Med 16, 31 (2021). https://doi.org/10.1186/s13020-021-00439-w
- Wang-Bi capsule
- Rheumatoid arthritis
- Joint destruction
- Molecular mechanism