Skip to main content

Traditional Chinese medicine Biqi capsule compared with leflunomide in combination with methotrexate in patients with rheumatoid arthritis: a randomized controlled trial



Biqi capsule is a traditional Chinese medicine widely used as a complementary and alternative treatment for rheumatoid arthritis (RA). The objective is to understand the efficacy, safety and mechanism of Biqi combined with methotrexate (MTX) in RA.


We present a randomized, controlled pilot trial on Biqi combined with MTX against patients with active RA. Seventy patients were randomized 1:1 to receive Biqi + MTX or Leflunomide (LEF) + MTX for 24 weeks, and were assessed at baseline, 4, 12 and 24 weeks. Serum and urine samples were collected for metabolomics.


Overall, 81.2% patients in Biqi group achieved ACR20 at 24 weeks. No statistically significant differences were observed in primary or secondary outcomes between the two groups. A better safety profile was observed for Biqi with significantly fewer adverse effects reported (11.4%) compared to LEF group (40%, P < 0.05). Comparison between treatment responders and non-responders indicated a unique urine metabolic profile of enriched fatty acids and decreased acylcarnitines associated with Biqi responders, indicating a restored energy homeostasis in response to Biqi. The gene targets of these metabolites were significantly enriched in interleukin-4 and interleukin-13 pathways, implying that Biqi could ameliorate Th2-derived inflammatory response. Multivariate network analysis indicated that patient morning stiffness and SJC were key factors associated with metabolomics in Biqi-treated patients, whereas CRP was the main factor in LEF group. Therefore, Biqi and LEF likely work by influencing different patient clinical parameters.


Our study suggests that Biqi capsule can be a promising alternative option in combination with MTX for RA treatment, and demonstrates the capability of using metabolomics to interrogate mechanism of action for traditional Chinese medicine.

Trial registration This trial is registered with ChiCTR, No. ChiCTR-IPR-16009029. Registered August 15, 2016.


Rheumatoid arthritis (RA) is an autoimmune disease characterized by progressive articular destruction and associated comorbidities in vascular, metabolic, bone, and psychological domains [1, 2]. Disease-modifying antirheumatic drugs (DMARDs) are recommended by the American College of Rheumatology (ACR) for the treatment of RA [3]. One of the most widely used DMARDs, Methotrexate (MTX) is the anchor drug for treatment of RA [4]. However, some patients failed to achieve low disease activity or clinical remission state after receiving MTX monotherapy [5, 6]. Leflunomide (LEF) is another DMARD with different mode of action from MTX [7]. It was considered as one of the fundamental therapeutics for RA by the European League Against Rheumatism (EULAR) in 2010 [3]. The combination of LEF and MTX (MTX + LEF) has been widely applied for treatment of refractory RA [8]. While it shows improved efficacy, there are also increasing risks of hepatotoxicity, bone marrow suppression, tuberculosis and infection [9,10,11,12], which has led an increasing number of RA patients to seek complementary and integrative medicine [13].

Traditional Chinese medicine (TCM) has been practiced in China for thousands of years and it offers a holistic approach to patient management. A variety of TCMs have been used for the treatment of RA in China for centuries, and their efficacy and safety were inferred from clinical experience and tested in clinical trials [14, 15]. Biqi capsule is a well-recognized traditional TCM formula that has been approved by the State Food and Drug Administration as an add-on therapy for RA due to its clinical benefits and low side effects [16]. Biqi capsule consists of ten TCM medicines: Strychni Semen (Maqianzi), Pheretima (Dilong), Codonopsis Radix (Dangshen), Poria (Fuling), Atractylodis Macrocephalae Rhizoma (Baizhu), Glycyrrhizae Radix Et Rhizoma (Gancao), Chuanxiong Rhizoma (Chuangxiong), Salviae Miltiorrhizae Radix Et Rhizoma (Danshen), Notoginseng Radix Et Rhizoma (Sanqi) and Cyathulae Radix (Chuanniuxi). Previous studies showed that Biqi capsule improved clinical symptoms and reduced inflammation in RA patients [17, 18]. A meta-analysis of ten randomized controlled trials (RCTs) revealed that combination treatment with MTX and Biqi (MTX + Biqi) was more effective than MTX monotherapy for RA. However, the majority of the included RCTs showed low methodological quality [19]. Despite its efficacy, little is known regarding the underlying mechanism of action of Biqi capsule for RA due to the complexity of its bioactive ingredients and paucity of knowledge on their therapeutic targets in human body.

Previously, we have explored potential mechanisms of action (MoA) for TCMs in treating RA and other rheumatic disorders through a number of pre-clinical studies including in silico network modeling [20], in vitro cellular assays [20, 21] and in vivo animal experiment [22]. Although insightful, results of these studies warrant further validation in the clinical settings. On the other hand, while a few RCTs have been conducted to assess the efficacy of TCMs in RA [23], a comprehensive understanding of MoA for TCMs by directly measuring the metabolic perturbation and immune response of patients has been lacking. Metabolomics emerge as a useful tool to provide novel insights into the therapeutic effects and underlying mechanisms of multi-biochemical component medicines. In addition, it endows the power to identify biomarkers indicative of treatment response. Therefore, a systems biology approach integrating measurements of clinical parameters and metabolomics (known as pharmacometabolomics [24]) in the RCTs offers the promise to interrogate MoA for TCMs in RA and identify novel biomarkers toward precision medicine.

The aim of this randomized, controlled pilot trial is to evaluate the efficacy and safety of MTX + Biqi compared to MTX + LEF for RA. We also employed a systems biology approach by exploring the differential MoAs for biqi and LEF in treating RA through collective serum and urine metabolomics. Results in this study will serve as a basis for future large-scale studies on this combination therapy and other TCMs in general.


Trial design

This study was a multi-centric, open-label, randomized trial conducted in four hospitals in Guangdong province in China, the Second Affiliated Hospital of Guangzhou University of Chinese Medicine, Shenzhen Hospital of Traditional Chinese Medicine, Dongguan Hospital of Traditional Chinese Medicine and Guangzhou Hospital of Integrated Traditional Chinese and Western Medicine, between September 2016 and July 2018. The study was conducted in accordance with the Declaration of Helsinki. All participants provided written informed consent and the protocol was approved by the Medical Ethics Committee of the Second Affiliated Hospital of Guangzhou University of Chinese Medicine (B2016-073-01) and registered with the World Health Organization clinical trial registry (No. ChiCTR-IPR-16009029).

Setting and participants

Eligible patients had to meet the following criteria for this trial: (1) aged 18–65 years; (2) RA was confirmed by 2010 ACR/EULAR classification criteria [25]; (3) imaging results suggested Class I, II or III disease (according to the 1987 American Rheumatism Association (ARA) classification standard [26]); (4) Chinese medicine inclusion criteria: according to the syndrome of wind and damp stagnation syndrome, or cold and damp stagnation, or phlegm and stasis stagnation, or deficiency of Qi and Blood [27]; (5) provided written informed consent.

Patients were excluded from the study if they: (1) used in the last month or are using glucocorticoids, MTX, hydroxychloroquine, willow nitrogen sulfanilamide pyridine, cyclophosphamide, penicillamine and gold preparations and other immunosuppressive drugs or slow-acting drugs; (2) had a history of cardiac, hepatic, renal, or mental diseases, other rheumatic autoimmune diseases, any current infection, or any cancer; (3) currently in pregnancy, or were planning on being pregnant during the study period; (4) were unwilling or unable to comply with treatment or assessment regimen; (5) experienced an allergic reaction to the medicine; (6) complicated with active gastrointestinal diseases or diagnosed with esophagus or digestive ulcer in the last month; (7) was participating in another clinical trial within 4 weeks prior to screening.

Randomization and interventions

Patients were randomly assigned 1:1 to receive either Biqi capsule at a dose of 1.2 g twice daily, or LEF at a dose of 20 mg once daily, for a period of 24 weeks. All included patients received MTX at a dose of 10–15 mg once weekly with 5 mg of folic acid two or three times a day. Participants were allowed to take glucocorticoids at 2.5–10 mg/day if their patient’s assessment of pain were more than 40 mm on a 100-mm visual analogue scale. In addition, nonsteroidal anti-inflammatory drugs (NSAIDs), calcitriol/calcium carbonate and antacids were allowed during the study. Patients could discontinue the assigned treatment at any time for lack of effectiveness or adverse effects or by their own choice.

Randomization was performed by an independent statistician using SAS v9.2. Allocation concealment was achieved using a centralized service. The statistician who prepared the list had no further role in the study. We did not use stratification or blocking. Blinding was not practical for this investigator-initiated clinical trial, and the treatments were open to both the researchers and participants.

Outcomes and measurements

Outcomes were assessed at baseline and after 4, 12 and 24 weeks. The primary efficacy endpoint was the proportion of patients who achieved a 20% improvement in ACR criteria (ACR20) at 24 weeks, defined as at least a 20% reduction from baseline in the number of both tender and swollen joints (TJC and SJC) and at least a 20% improvement in three or more of the following: the physician’s or patient’s global assessment of disease activity (PhGADA or PaGADA), the patient’s assessment of pain on a visual analogue scale (VAS), the patient’s assessment of function using the Health Assessment Questionnaire (HAQ), and the serum erythrocyte sedimentation rate (ESR) or C-reactive protein (CRP) level [28]. Secondary endpoints included the ACR50 and ACR70, EULAR response, and other clinical parameters (supplementary methods). The EULAR response criteria were assessed based on the 28-joint disease activity score (DAS28). A moderate EULAR response was a decrease of > 0.6 and ≤ 1.2, and a good response is a decrease of > 1.2 [29]. All of the outcome measures were assessed and recorded at baseline, weeks 4, 12 and 24 by trained evaluators who were unaware of the specific therapeutic regimen.

Safety endpoints included adverse events (AEs), serious AEs, and laboratory abnormalities. Safety was monitored by physical examination, chest radiography, electrocardiography, blood pressure, pulse rate and body temperature. Standard hematological and biochemical tests and urinalysis were also performed. The occurrence of adverse events was documented and included those spontaneously reported by patients, as well as responses elicited by general questioning.

Statistical analysis

The sample size was calculated based on the primary outcome ACR20 [30, 31]. It was estimated that 98 participants (49 per arm) is sufficient to detect a 20% difference in ACR20, using a 5% significance level with 80% power from baseline to end of treatment, allowing for a 20% loss to follow-up. Full analysis set was used to test differences in baseline characteristics of patients between two groups. All efficacy endpoints of the two group were analyzed using the intention-to-treat (ITT) and per-protocol set (PP). The ITT analysis included all patients who received at least one dose of the drug, and was performed as the primary outcome analysis. The PP analysis included patients who strictly followed the protocol and completed the study at each time point.

The mean ± standard deviation were used for the description of continuous variables. The proportion was used for the description of categorical variables. Dichotomous variables were analyzed using Fisher’s exact test for the primary efficacy endpoints (ACR), the EULAR response, as well as the safety endpoints. The secondary efficacy endpoints, including the TJC28, SJC28, patient’s assessment of pain, PaGADA, PhGADA, CRP and ESR, were assessed by one-way repeated measures ANOVA of the mean values from baseline to weeks 4, 12 and 24 for each group, followed by Benjamin-Hochberg (BH) post hoc correction to address the multiple comparisons on clinical measures and in between timepoints. The differences between the two treatment groups in these secondary efficacy endpoints were analyzed by the non-parametric Mann–Whitney U test to accommodate their non-normal distributions. All statistical tests were 2-sided. P-values less than 0.05 were considered significant. An independent statistician who was blinded to the group allocation conducted the data analysis using SPSS software (version 17.0).

Serum and urine metabolomic analysis

Serum and urine samples were collected from patients at baseline, weeks 4, 12 and 24, and were analyzed using Agilent 1290 ultra-performance liquid chromatography mass spectrometry (UHPLC-MS). Identification of metabolites was performed by searching HMDB [32], METLIN [33] and KEGG databases [34]. Metabolomic data analysis was performed using MetaboAnalyst [35]. The data were normalized by sum, log2 transformed and auto-scaled. Analysis of similarities (ANOSIM) was performed to test the dissimilarity between metabolic profiles among groups. Differentially abundant metabolites were identified for paired samples in between different timepoints using one-way repeated measures ANOVA. Differential metabolites at 24 weeks between treatment responders and non-responders were identified using analysis of covariance (ANCOVA), in which the metabolite abundances at baseline were controlled as a covariate. Human gene targets of the metabolites were identified using the STITCH database [36]. Metabolite-target interactions with activation or inhibition effects and with confidence score > 0.8 were retained. The fold-changes of targets were inferred as the additive fold-change of all metabolites targeting the genes, adjusted by activation or inhibition effects of the metabolites to the genes.

Correlation was performed between serum and urine metabolites and patient clinical variables. Each metabolite and clinical variable was first residualized using a linear mixed-effect model with age, gender and timepoint as fixed-effect variables and patient identity as random-effect variable. A hierarchical spearman correlation was performed on the residues of metabolites and clinical variables using HAllA [37]. Correlation networks were established using the top 50 correlations between metabolites and clinical variables.


Patient characteristics

The main procedure of the clinical trial was illustrated in Fig. 1. Of 100 screened patients, 70 eligible patients were enrolled in the study and 1:1 randomly assigned to one of the two treatment arms (Fig. 2). Fifty-nine (84.3%) patients completed 24 weeks of evaluation (32 in Biqi and 27 in LEF group). There were no statistically significant differences between two groups in demographics, disease characteristics or concomitant medications at baseline (Table 1, Additional file 1: Table S1).

Fig. 1
figure 1

The main procedure for the clinical trial

Fig. 2
figure 2

Consort diagram for patient participation flow throughout the trial

Table 1 The demographic and clinical characteristics of RA patients at baseline in the full analysis set

Comparable clinical efficacy for Biqi and LEF treatment groups

Overall similar clinical efficacy was observed for the two treatment arms. 81.2% of patients in Biqi group and 81.5% of patients in LEF group achieved ACR20 response as the primary outcome at 24 weeks (Fig. 3a). Similar ACR50 and ACR70 rates were observed at 24 weeks for the two groups. The ACR20, ACR50 and ACR70 rates were also comparable between two groups during 4 and 12 weeks, except for a notable higher ACR50 rate for Biqi compared to LEF group at 12 weeks (51.5% versus 35.4%, P = 0.17). There were slightly higher patient EULAR response rates in LEF compared to Biqi group (Fig. 3b). Compared to baseline, most clinical measurements showed significant and continuous improvement in both two groups over time (Fig. 3d, FDR-adjusted P < 0.05). No significant differences were found between the two groups. PP analysis on the 59 patients that completed 24 weeks of treatment showed consistent results with those in the ITT analysis (Additional file 1: Fig. S1).

Fig. 3
figure 3

Patient clinical outcomes for Biqi and LEF group in the ITT analysis. These include a patient response rates for ACR20, ACR50, ACR70, b patient response rates for EULAR good or moderate response, c proportion of patients with reported adverse events, and d patient clinical measurements at baseline, 4, 12 and 24 weeks. Paired Wilcoxon test was performed followed by Benjamin-Hochberg (BH) post hoc correction for multiple comparisons on clinical measures and in between timepoints. The FDR-adjusted P-values between 24 weeks and baseline were reported. *** FDR P < 0.001, ** P < 0.01, * P < 0.05

A better safety profile for Biqi treatment

Of all 70 patients, 18 (25.7%) experienced one or more adverse events (AEs), including 4 (11.4%) and 14 (40%) patients receiving Biqi and LEF, respectively (Fig. 3c, P = 0.006). In Biqi group, all AEs were regarding to elevated liver enzymes above the normal range. In LEF group, the AEs included hepatic side effects (9 cases), abdominal discomfort (2 cases), hypertension (1 case), rash (1 case), and herpes zoster (1 case). In particular, 5 patients in LEF group suffered from the increasing ALT/AST and discontinued the study. Moreover, 12 of the 14 LEF-treated patients with AEs showed symptoms at week 0 or week 4, whereas all AEs were for Biqi-treated patients were observed at week 12 or 24 (week 4: P = 2e−3, Fig. 3c). These results indicate a milder adverse effect of Biqi possibly due to its slow-acting nature, as opposed to LEF which likely provoked an acute host response at the onset of treatment.

Serum and urine metabolomics over time for Biqi and LEF treatment arms

To investigate potential mechanism of action for Biqi capsule, we carried out metabolomic analysis for patient serum and urine samples. A total of 106 serum and 103 urine samples were collected from 31 to 27 patients respectively (Additional file 1: Table S2). Quality assessment indicated that all samples had good repeatability and stability (Additional file 1: Fig. S2). A total of 14,956 serum and 18,775 urine metabolites were detected in UHPLC-MS, of which 275 and 474 metabolites were successfully resolved to identity. There were significant serum metabolomic shifts at 24 weeks in both Biqi and LEF groups as indicated in PCA plots (ANOSIM 24w vs baseline: P < 0.05, Fig. 4a, b). Similar trend was observed for urine samples of Biqi-treated patients, whereas for LEF-treated patients both week 12 and 24 samples together formed a separate cluster from other samples. Consistently, the majority of differentially abundant serum and urine metabolites were identified at 24 weeks for Biqi-treated patients, while for LEF-treated patients the number of differential urine metabolites peaked at 12 weeks (Additional file 1: Fig. S3). The different timing for patient metabolomic shifts supported the slow-acting nature of Biqi compared to LEF.

Fig. 4
figure 4

Serum and urine metabolic shifts over time for Biqi and LEF groups. a PCA plots for serum and urine samples at baseline, 4, 12 and 24 weeks for patients receiving Biqi or LEF treatment respectively. Samples were colored by timepoints. The ellipses represent 95% confidence interval for samples within each timepoint. b Volcano plots showing differentially abundant serum and urine metabolites comparing 24 weeks versus baseline for Biqi and LEF-treated patients. c Venn diagrams comparing serum and urine metabolites up and down-regulated at 24 weeks versus baseline for Biqi and LEF groups. d Top serum and urine metabolites that were commonly downregulated in Biqi and LEF groups

Next, we identified differentially abundant metabolites at 24 weeks versus baseline for Biqi and LEF groups. Comparison of the differential metabolites indicates a large proportion of serum and urine metabolites commonly altered in both groups. For example, 136 out of 195 (69.7%) differential serum metabolites were commonly altered in both Biqi and LEF groups at 24 weeks (FDR P < 0.05, log2 fold-change (FC) > 1). Of them, 134 metabolites were downregulated (Fig. 4c, Additional file 1: Table S3). Likewise, 77 urine metabolites were commonly altered in both groups, of which 66 metabolites were downregulated (Fig. 4c, Additional file 1: Table S3). Amino acid and sphingolipid metabolisms were most significantly enriched for the downregulated metabolites in both serum and urine (Fig. Additional file 1: S4a, b). P-Chlorophenylalanine was most significantly decreased in serum and urine for both groups (Fig. 4d). Of note, for the metabolites that were deemed uniquely up or downregulated in each treatment group, they mostly showed similar trend of changes in the other group, despite at a smaller scale that didn’t reach statistical significance (Additional file 1: Fig. S5). Therefore, overall patients in Biqi and LEF groups appeared to have similar serum and urine metabolomic shifts over time.

Urine metabolomic profiles associated with patient responsiveness to Biqi treatment

To further disentangle patient metabolomic shifts as a result of Biqi or LEF treatment alone, we divide both patient groups into responders (R) and non-responders (NR) based on their clinical outcomes at 24 weeks. We chose to use the ACR50 as the criteria for treatment responsiveness which divided patients into two subgroups with roughly even sample size. As expected, PCA indicated significant serum and urine metabolic shifts for both responders and non-responders at week 24 versus baseline for both treatment groups. There was no clear separation between R and NR groups for serum metabolites of both Biqi and LEF-treated patients, and for urine metabolites of LEF-treated patients (Fig. 5a). However, for urine metabolites of Biqi-treated patients, the R and NR were well separated in PCA both at baseline and 24 weeks, indicating a significant association of urine metabolomic profiles with Biqi treatment responsiveness (ANOSIM P = 0.039, Fig. 5a, b).

Fig. 5
figure 5

Serum and urine metabolomic profiles for responders and non-responders to Biqi and LEF treatment. a PCA plots for paired serum and urine samples at baseline and 24 weeks for patients receiving Biqi or LEF treatment respectively. Samples were colored by timepoints. The paired baseline and 24 weeks samples from the same patient were linked by solid lines colored by responders or non-responders. The dotted and solid ellipses represent 95% confidence interval for samples at baseline and 24 weeks, respectively. b Volcano plots showing differential abundant serum and urine metabolites comparing responders versus non-responders at 24 weeks in Biqi and LEF groups using ANCOVA. c Venn diagrams comparing serum and urine metabolites up and down-regulated in responders versus non-responders at 24 weeks for Biqi and LEF groups. d Top urine carboxylic acids and acylcarnitines that were uniquely up and downregulated in Biqi responders versus non-responders at 24 weeks

We then used an ANCOVA approach to identify metabolites that significantly differed between R and NR for Biqi and LEF treatments at 24 weeks, controlling for their initial abundances at baseline. A total of 152 urine metabolites were significant different in abundance between Biqi R and NR, 111 uniquely upregulated in Biqi responders (Fig. 5c, Additional file 1: Table S4). In comparison, only 18 metabolites were significantly altered between LEF R and NR at 24 weeks. The top downregulated metabolite in LEF responders was uracil, which was congruent with the established mechanism of action for LEF in inhibiting synthesis of pyrimidine ribonucleotide uridine monophosphate [38]. The upregulated metabolites in Biqi responders were most enriched in propionate and glutathione metabolisms (Additional file 1: Fig. S4c). Meanwhile, multiple medium or long-chain carboxylic acids such as suberic acid, decanoic acid, docosatrienoic acid and hexadecanedioic acid were significantly upregulated (Fig. 5d). The 37 downregulated metabolites in Biqi responders included a wide range of acylcarnitines, including butenylcarnitine, pimelylcarnitine, hexanoycarnitine and propionylcarnitine (Fig. 5d, Additional file 1: Table S4).

To understand how perturbations of metabolites might impact host immune and inflammatory response, we performed an in silico analysis to identify human gene targets for these 152 metabolites using the STITCH database [36]. A total of 167 human gene targets were identified as potentially activated or inhibited by these metabolites (confidence score > 0.8, Additional file 1: Fig. S6). Of them, 91 genes were inferred to be downregulated in Biqi R versus NR at 24 weeks. These genes were most significantly enriched in interleukin (IL)-4 and IL-13 signaling pathways, key Th2-derived anti-inflammatory pathways (Additional file 1: Fig. S7). Conversely, no pathway was significantly enriched for the 70 upregulated genes.

Distinct associations between metabolites and clinical parameters for Biqi and LEF treatments

To further link metabolic alternations to patient clinical outcomes, we performed a correlation network analysis between serum and urine metabolites and patient clinical parameters. We used residualized correlation to identify correlations independent of longitudinal measurements and other demographic co-factors. We observed distinct correlation network modules between serum and urine metabolites and clinical variables for Biqi and LEF groups (Fig. 6). For Biqi-treated patients, SJC and morning stiffness were key clinical factors associated with both serum and urine metabolites (Fig. 6a, c). For example, serum l-glutamine, vitamin-D3, decanoic acid and prostaglandin F1a showed significant negative correlations with patient morning stiffness levels (Fig. 6a). Notably, these metabolites were uniquely upregulated in Biqi responders at 24 weeks, suggesting such correlations were also related to Biqi responsiveness. Conversely, multiple urine acylcarnitines, all downregulated in Biqi responders, showed positive correlations with morning stiffness (Fig. 6c). For LEF-treated patients, the main clinical factor associated with metabolites was CRP (Fig. 6b, d). In particular, there was a mutual positive association between multiple serum lysophosphatidylcholines (lysoPCs) that together were positively correlated with CRP (Fig. 6b). Urine metabolites such as ethanolamide and octadecanoid acid were also negatively associated with CRP in LEF group (Fig. 6d).

Fig. 6
figure 6

Interaction networks between metabolites and patient clinical parameters in Biqi and Lef group. This includes a serum of Biqi-treated patients, b serum of LEF-treated patients, c urine of Biqi-treated patients, and d urine of LEF-treated patients. Each node in round shape represents a metabolite colored based on its fold-change in responders versus non-responders at 24 weeks. Each node in diamond shape represents one clinical measurement. Each edge represents a significant correlation between the two nodes (FDR P < 0.05, residual correlation), colored by positive or negative correlations. For visualization purpose, only top 50 significant correlations between metabolites and patient clinical variables, along with any significant correlations within these metabolites are shown


This 24-week pilot trial confirmed both efficacy and tolerability for Biqi, one of the commonly prescribed TCM formulas, in combination with MTX for RA treatment. Intergroup comparisons suggested an approximately equivalent efficacy of Biqi when compared with LEF that had demonstrated therapeutic efficacy for patients with refractory or active RA [39]. The proportion of patients who attained the primary and secondary endpoints were comparable between LEF and Biqi groups at 4, 12 and 24 weeks. A slightly higher proportion of patients achieved EULAR response in LEF than Biqi group, whereas higher ACR50 was observed for Biqi group in 12 weeks. All clinical measurements were significantly improved at 24 weeks in both groups, except for the non-significant changes in CRP, ESR and rheumatoid factor (RF) in Biqi group. Together these results indicated a clear clinical efficacy for MTX + Biqi comparable to MTX + LEF combination therapy.

Importantly, our results indicated that MTX + Biqi had a better safety profile. The frequency of total adverse events in Biqi group was significantly lower than that of LEF group. Furthermore, all adverse events in Biqi group occurred after 4 weeks, as opposed to those in LEF that occurred at the onset of the treatment, indicating a slow-acting nature for Biqi capsule. Liver abnormality was the most frequently reported adverse events in both treatments. All patients in Biqi group could be immediately relieved after treatment with hematinics, while 14.3% patients in the LEF group discontinued the trial. Reported adverse effects suggested that there were concerns regarding the toxicity of Biqi owing to its active ingredients such as strychnine and brucine [40]. Future studies focused on the possible adverse effects of Biqi in vital organs and systems are therefore highly essential.

The serum and urine metabolomic profiling provided many insights into the potential mechanisms of action for Biqi that might be different from those of LEF. Overall there were similar serum and urine metabolic changes with decreased amino acid and sphingolipid metabolism in both groups when comparing 24 weeks versus baseline. This might be explained as the effect of concomitant MTX treatment that could have masked that of Biqi and LEF. Further patient stratification revealed a unique association of urine metabolomics with Biqi responsiveness, which was characterized by enriched carboxylic acids and decreased acylcarnitines. As acylcarnitines constitute the rate-limiting step for fatty acid beta-oxidation, the increased fatty acids coupled with decreased carnitine together indicates a restored energy homeostasis accompanied with Biqi responsiveness, which may also reflect suppressed host immune responses by Biqi [41]. Accordingly, in silico analysis showed that these metabolites might inhibit cytokines in Th2-derived (IL-4 and IL-13) inflammatory pathways, including IL-1, IL-4, IL-6, IL-8 and IL-13. Bian et al. showed that Biqi inhibited the progression of collagen-induced arthritis by decreasing JAK3, STAT3, IL–4, IL–1 and other inflammatory cytokines in RA rat models [42]. Other clinical studies also confirmed that Biqi markedly decreases IL–17, IL–1 and TNF–a in RA patients [43, 44]. Our results extended these findings in the clinical setting and further supported the anti-inflammatory properties of Biqi capsule.

Using a systems biology approach integrating metabolomics and patient clinical parameters, we showed that there were different types of interactions between metabolites and disease activities in Biqi and LEF groups. In Biqi-treated patients, patient morning stiffness levels were positively correlated with serum anti-inflammatory agents such as prostaglandin F1a [45], vitamin K3 [46] and l-glutamine [47] while negatively correlated with urine acylcarnitines. Importantly, many of these metabolites were also associated with Biqi responsiveness (such as urine acylcarnitines), which further supported that they were key components mediating patient clinical response to Biqi. On the other hand, CRP was the main factor associated with metabolites in the LEF group. The anti-CRP effect of LEF in RA is well established [48] and our interaction network implies that LEF could potentially reduce CRP levels via decreasing lysoPCs in serum, and increasing creatine, guanidinoacetic acid and thiodiacetic acid in urine. Our analysis suggests that Biqi and LEF likely achieved anti-RA effects through impacting different types of patient clinical characteristics. The integrated patient stratification and clinical correlation analyses allowed us to navigate the important metabolites that most likely played mechanistic roles underlying Biqi treatment.

Metabolomics has been increasingly applied to animal experiments and clinical studies on TCMs to characterize their therapeutic efficacies and mechanisms of action in treating RA [49]. For instance, by applying GC–MS based metabolomic approach in a rat model, Wang Y. et al. investigated the effect of Simiaowan and showed that it achieved the anti-arthritic effect by decreasing amino acid metabolism and restoring energy metabolism [50]. Some of their findings were also observed in our results. Comparative studies with additional TCMs in RA are therefore of interest to further determine whether there might be shared or specific metabolic effects between different TCMs in the treatment of RA. In another study, using combined untargeted and targeted metabolomic methods, Wang M. et al. examined the plasma metabolic perturbation in patients receiving a combo therapy of MTX and Chinese herbal medicine tripterygium glycosides, and showed that it modified a wide range of human metabolism including amino acid, nucleic acid, lipid metabolisms and oxidative stress [51]. By comparing the combo-therapy group with mono-therapy groups, they further showed the combo treatment may have specific effects on the levels of alanine, adenosine, and lactic acid. Their and our studies both demonstrated the importance of a rigorous clinical trial design with proper control arm, longitudinal follow-up, and a finer-scale patient stratification analysis, in order to distinguish the metabolic effects specific to TCMs from those non-specific effects from other concomitant treatments or interpersonal variations.

This pilot study has several limitations. Firstly, our study was an open-label clinical trial. In order to achieve more objective outcomes, a double-blind RCT would be necessary. The second one is the duration of the study. According to the EULAR guiding principle for RA treatment, it is usually admitted that a 24 weeks duration study is sufficient to evaluate the statistically significant efficacy of a therapeutic regimen [52]. However, the evaluation of radiological progression of the affected joints in X-ray is not available and requires probably studies of longer duration. The third one is the relatively small cohort size, in particular with respect to safety evaluation and sub-analysis of treatment responsiveness. Future studies with double-blind design, larger group of patients and a prolonged period of time are needed to validate findings of this pilot study and to identify the subgroup of patients best responsive to Biqi capsule.


Our study showed that Biqi capsule in combination of MTX achieved equivalent efficacy with MTX + LEF in the treatment of RA with a significantly better safety profile. We showed that Biqi capsule might ameliorate RA symptoms through restoring urine energy homeostasis and suppressing Th2 inflammatory response. The improvement in symptoms and patient quality of life makes this combination therapy an attractive alternative treatment to currently available agents for RA patient management. However, the long-term effects and toxicities of Biqi capsule in combination with MTX or other antirheumatic therapies need to be addressed in future studies.

Availability of data and materials

The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.



American College of Rheumatology


Analysis of covariance


Analysis of similarities


C-reactive protein


28-joint disease activity score


Disease-modifying antirheumatic drugs


European league against rheumatism


Erythrocyte sedimentation rate


Health Assessment Questionnaire










Patient’s global assessment of disease activity


Physician’s global assessment of disease activity




Rheumatoid arthritis


Rheumatoid factor


Visual analogue scale


  1. Scott DL, Wolfe F, Huizinga TW. Rheumatoid arthritis. Lancet. 2010;376(9746):1094–108.

    Article  PubMed  Google Scholar 

  2. McMurdie PJ, Holmes S. Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput Biol. 2014;10(4):e1003531.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  3. Smolen JS, Landewe R, Bijlsma J, Burmester G, Chatzidionysiou K, Dougados M, Nam J, Ramiro S, Voshaar M, van Vollenhoven R, et al. EULAR recommendations for the management of rheumatoid arthritis with synthetic and biological disease-modifying antirheumatic drugs: 2016 update. Ann Rheum Dis. 2017;76(6):960–77.

    Article  PubMed  Google Scholar 

  4. Brown PM, Pratt AG, Isaacs JD. Mechanism of action of methotrexate in rheumatoid arthritis, and the search for biomarkers. Nat Rev Rheumatol. 2016;12(12):731–42.

    Article  CAS  PubMed  Google Scholar 

  5. Aletaha D, Smolen JS. Joint damage in rheumatoid arthritis progresses in remission according to the Disease Activity Score in 28 joints and is driven by residual swollen joints. Arthritis Rheum. 2011;63(12):3702–11.

    Article  PubMed  Google Scholar 

  6. Gandjbakhch F, Conaghan PG, Ejbjerg B, Haavardsholm EA, Foltz V, Brown AK, Moller Dohn U, Lassere M, Freeston J, Boyesen P, et al. Synovitis and osteitis are very frequent in rheumatoid arthritis clinical remission: results from an MRI study of 294 patients in clinical remission or low disease activity state. J Rheumatol. 2011;38(9):2039–44.

    Article  PubMed  Google Scholar 

  7. Kremer JM. Methotrexate and leflunomide: biochemical basis for combination therapy in the treatment of rheumatoid arthritis. Semin Arthritis Rheum. 1999;29(1):14–26.

    Article  CAS  PubMed  Google Scholar 

  8. Scott DL, Ibrahim F, Farewell V, O’Keeffe AG, Walker D, Kelly C, Birrell F, Chakravarty K, Maddison P, Heslin M, et al. Tumour necrosis factor inhibitors versus combination intensive therapy with conventional disease modifying anti-rheumatic drugs in established rheumatoid arthritis: TACIT non-inferiority randomised controlled trial. BMJ. 2015;350:h1046.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Toth P, Bernd R. Severe leukopenia in a rheumatoid arthritis patient treated with a methotrexate/leflunomide combination. Rev Bras Reumatol. 2014;54(2):152–4.

    Article  PubMed  Google Scholar 

  10. Rutanen J, Kononoff A, Arstila L, Elfving P, Koskela H, Kaipiainen-Seppanen O. Five cases of interstitial lung disease after leflunomide was combined with methotrexate therapy. Scand J Rheumatol. 2014;43(3):254–6.

    Article  CAS  PubMed  Google Scholar 

  11. Criner GJ, Connett JE, Aaron SD, Albert RK, Bailey WC, Casaburi R, Cooper JA Jr, Curtis JL, Dransfield MT, Han MK, et al. Simvastatin for the prevention of exacerbations in moderate-to-severe COPD. N Engl J Med. 2014;370(23):2201–10.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Jenks KA, Stamp LK, O’Donnell JL, Savage RL, Chapman PT. Leflunomide-associated infections in rheumatoid arthritis. J Rheumatol. 2007;34(11):2201–3.

    CAS  PubMed  Google Scholar 

  13. Moudgil KD, Berman BM. Traditional Chinese medicine: potential for clinical treatment of rheumatoid arthritis. Expert Rev Clin Immunol. 2014;10(7):819–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. He YT, Ou AH, Yang XB, Chen W, Fu LY, Lu AP, Yan XP, Feng XH, Su L, Song YJ, et al. Traditional Chinese medicine versus western medicine as used in China in the management of rheumatoid arthritis: a randomized, single-blind, 24-week study. Rheumatol Int. 2014;34(12):1647–55.

    Article  PubMed  Google Scholar 

  15. Lee S, Cho Y, Kim J, Kang JW, Yoon GY, Lee JH, Jung SY, Kwon O, Shin KM, Lee JD. The efficacy and safety of the herbal medicine geonchildan for patients with active rheumatoid arthritis: study protocol for a randomized, double-blind, placebo-controlled, parallel pilot trial. Trials. 2018;19(1):471.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. Liu W, Zhang L, Xu Z. Clinical observation on treatment of rheumatoid arthritis with biqi capsule. Zhongguo Zhong Xi Yi Jie He Za Zhi. 2006;26(2):157–9.

    PubMed  Google Scholar 

  17. Wu Q, Ding Z, Lu Y, Sun D, Feng L, Jie H, Su X, Chen W, Tang H. Clinical study of Biqi capsule for treatment of rheumatism. In: Proceedings of the 8th national conference on rheumatic diseases in Traditional Chinese Medicine and Western Medicine. 2010; p. 104–109. (In Chinese).

  18. Yun C. Treatment of rheumatoid arthritis with Biqi capsule in 36 cases. J Pract Tradit Chin Med. 2004;20:509 (In Chinese).

    Google Scholar 

  19. Chen XM, Wu JQ, Huang QC, Zhang JY, Pen JH, Huang ZS, Chu YL, He XH, Wang MJ, Huang RY. Systematic review and meta-analysis of the efficacy and safety of Biqi capsule in rheumatoid arthritis patients. Exp Ther Med. 2018;15(6):5221–30.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. Wang Z, Linghu KG, Hu Y, Zuo H, Yi H, Xiong SH, Lu J, Chan G, Yu H, Huang RY. Deciphering the Pharmacological Mechanisms of the Huayu-Qiangshen-Tongbi Formula Through Integrating Network Pharmacology and In Vitro Pharmacological Investigation. Front Pharmacol. 2019;10:1065.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Huang RY, Lin JH, He XH, Li X, Lu CL, Zhou YY, Cai J, He YT. Anti-inflammatory activity of extracts of Bushen-Qiangdu-Zhilv decoction, a Chinese medicinal formula, in M1-polarized RAW264.7. BMC Complement Altern Med. 2014;14:268.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Zhou YY, Huang RY, Lin JH, Xu YY, He XH, He YT. Bushen-Qiangdu-Zhilv decoction inhibits osteogenic differentiation of rat fibroblasts by regulating connexin 43. Exp Ther Med. 2016;12(1):347–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Huang RY, Pan HD, Wu JQ, Zhou H, Li ZG, Qiu P, Zhou YY, Chen XM, Xie ZX, Xiao Y, et al. Comparison of combination therapy with methotrexate and sinomenine or leflunomide for active rheumatoid arthritis: a randomized controlled clinical trial. Phytomedicine. 2019;57:403–10.

    Article  CAS  PubMed  Google Scholar 

  24. Kaddurah-Daouk R, Weinshilboum RM, Pharmacometabolomics Research N. Pharmacometabolomics: implications for clinical pharmacology and systems pharmacology. Clin Pharmacol Ther. 2014;95(2):154–67.

    Article  CAS  PubMed  Google Scholar 

  25. Aletaha D, Neogi T, Silman AJ, Funovits J, Felson DT, Bingham CO 3rd, Birnbaum NS, Burmester GR, Bykerk VP, Cohen MD, et al. 2010 Rheumatoid arthritis classification criteria: an American College of Rheumatology/European League Against Rheumatism collaborative initiative. Arthritis Rheum. 2010;62(9):2569–81.

    Article  PubMed  Google Scholar 

  26. Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, Cooper NS, Healey LA, Kaplan SR, Liang MH, Luthra HS, et al. The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum. 1988;31(3):315–24.

    Article  CAS  PubMed  Google Scholar 

  27. Liu T. Chinese medicine industry standard of the People’s Republic of China. Criteria for the diagnosis and treatment of diseases and syndromes in traditional Chinese medicine 2006.

  28. Felson DT, Anderson JJ, Boers M, Bombardier C, Furst D, Goldsmith C, Katz LM, Jr. Lightfoot R, Paulus H, Strand V, et al. American College of Rheumatology Preliminary definition of improvement in rheumatoid arthritis. Arthritis Rheum. 1995;38(6):727–35.

    Article  CAS  PubMed  Google Scholar 

  29. Fransen J, van Riel PL. The Disease Activity Score and the EULAR response criteria. Rheum Dis Clin North Am. 2009;35(4):745–57.

    Article  PubMed  Google Scholar 

  30. Strand V, Cohen S, Schiff M, Weaver A, Fleischmann R, Cannon G, Fox R, Moreland L, Olsen N, Furst D, Leflunomide Rheumatoid Arthritis Investigators Group, et al. Treatment of active rheumatoid arthritis with leflunomide compared with placebo and methotrexate. Arch Intern Med. 1999;159(21):2542–50.

    Article  CAS  PubMed  Google Scholar 

  31. Brien S, Lachance L, Prescott P, McDermott C, Lewith G. Homeopathy has clinical benefits in rheumatoid arthritis patients that are attributable to the consultation process but not the homeopathic remedy: a randomized controlled clinical trial. Rheumatology. 2011;50(6):1070–82.

    Article  PubMed  Google Scholar 

  32. Wishart DS, Feunang YD, Marcu A, Guo AC, Liang K, Vazquez-Fresno R, Sajed T, Johnson D, Li C, Karu N, et al. HMDB 4.0: the human metabolome database for 2018. Nucleic Acids Res. 2018;46(D1):D608–17.

    Article  CAS  PubMed  Google Scholar 

  33. Guijas C, Montenegro-Burke JR, Domingo-Almenara X, Palermo A, Warth B, Hermann G, Koellensperger G, Huan T, Uritboonthai W, Aisporna AE, et al. METLIN: a technology platform for identifying knowns and unknowns. Anal Chem. 2018;90(5):3156–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 1999;27(1):29–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Chong J, Soufan O, Li C, Caraus I, Li S, Bourque G, Wishart DS, Xia J. MetaboAnalyst 4.0: towards more transparent and integrative metabolomics analysis. Nucleic Acids Res. 2018;46(W1):W486–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Szklarczyk D, Santos A, von Mering C, Jensen LJ, Bork P, Kuhn M. STITCH 5: augmenting protein-chemical interaction networks with tissue and affinity data. Nucleic Acids Res. 2016;44(D1):D380–4.

    Article  CAS  PubMed  Google Scholar 

  37. Rahnavard G, Franzosa EA, Mclver LJ, Schwager E, Lloyd-Price J, Weingart G, Moon YS, Morgan XC, Waldron L, Huttenhower C. High-sensitivity pattern discovery in large multi omic datasets. huttenhowersphharvardedu/halla 2019.

  38. Fox RI, Herrmann ML, Frangou CG, Wahl GM, Morris RE, Strand V, Kirschbaum BJ. Mechanism of action for leflunomide in rheumatoid arthritis. Clin Immunol. 1999;93(3):198–208.

    Article  CAS  PubMed  Google Scholar 

  39. Antony T, Jose VM, Paul BJ, Thomas T. Efficacy and safety of leflunomide alone and in combination with methotrexate in the treatment of refractory rheumatoid arthritis. Indian J Med Sci. 2006;60(8):318–26.

    Article  CAS  PubMed  Google Scholar 

  40. Zheng H, Wang Z, Liu W, Jin H, Zhang J. Toxicokinetics of strychnine and brucine after the oral administration of Biqi capsule to rats by RRLC-MS/MS. Biomed Chromatogr. 2018;32(3):e4117.

    Article  CAS  Google Scholar 

  41. Seca S, Franconi G. Understanding Chinese medicine patterns of rheumatoid arthritis and related biomarkers. Medicines. 2018;5(1):17.

    Article  PubMed Central  CAS  Google Scholar 

  42. Bian XQ, Liu W. Experimental study on pharmacodynamics of Biqi capsule for treatment of rheumatoid arthritis. China J Tradit Chin Med Pharm. 2006;21:773–4 (In Chinese).

    Google Scholar 

  43. Liu W, Zhang L, Zhou YL, Liu XY. Influence of Biqi capsule on cytokine in rheumatoid arthritis. Tianjin J Tradit Chin Med. 2005;22:290–1 (In Chinese).

    Google Scholar 

  44. Wan HJ, Liu J, Xu JY, Fang HY. The efficacy of Biqi capsules combined with low dose methotrexate in the treatment of rheumatoid arthritis. Chin J Biochem Pharm. 2012;33:669–71 (In Chinese).

    Google Scholar 

  45. Scher JU, Pillinger MH. The anti-inflammatory effects of prostaglandins. J Investig Med. 2009;57(6):703–8.

    Article  CAS  PubMed  Google Scholar 

  46. Checker R, Sharma D, Sandur SK, Khan NM, Patwardhan RS, Kohli V, Sainis KB. Vitamin K3 suppressed inflammatory and immune responses in a redox-dependent manner. Free Radic Res. 2011;45(8):975–85.

    Article  CAS  PubMed  Google Scholar 

  47. Jain P, Khanna NK. Evaluation of anti-inflammatory and analgesic properties of l-glutamine. Agents Actions. 1981;11(3):243–9.

    Article  CAS  PubMed  Google Scholar 

  48. Kraan MC, Smeets TJ, van Loon MJ, Breedveld FC, Dijkmans BA, Tak PP. Differential effects of leflunomide and methotrexate on cytokine production in rheumatoid arthritis. Ann Rheum Dis. 2004;63(9):1056–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Wang P, Wang Q, Yang B, Zhao S, Kuang H. The progress of metabolomics study in Traditional Chinese Medicine Research. Am J Chin Med. 2015;43(7):1281–310.

    Article  CAS  PubMed  Google Scholar 

  50. Wang Y, Guo X, Xie J, Hou Z, Li Y. A GC-MS based metabonomics study of rheumatoid arthritis and the interventional effects of the Simiaowan in rats. Molecules. 2015;20(12):21364–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Wang M, Huang J, Fan H, He D, Zhao S, Shu Y, Li H, Liu L, Lu S, Xiao C, et al. Treatment of rheumatoid arthritis using combination of methotrexate and tripterygium glycosides tablets—a quantitative plasma pharmacochemical and pseudotargeted metabolomic approach. Front Pharmacol. 2018;9:1051.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Smolen JS, Breedveld FC, Burmester GR, Bykerk V, Dougados M, Emery P, Kvien TK, Navarro-Compan MV, Oliver S, Schoels M, et al. Treating rheumatoid arthritis to target: 2014 update of the recommendations of an international task force. Ann Rheum Dis. 2016;75(1):3–15.

    Article  PubMed  Google Scholar 

Download references


Not applicable.


This study was supported by National Natural Science Foundation of China (No. 81774218), the Science and Technology Foundation of Guangdong Province (No. 2016A020226041), the Key Research Project of Guangzhou University of Chinese Medicine (No. XK2019021), the Science and Technology Foundation of Guangzhou city (No. 201710010076), as well as Chinese Medical Science and Technology research funding from Guangdong Provincial Hospital of Chinese Medicine (No. YN2018ZD06). Funders had no role in study design, collection, analysis and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations



ZW, JW and RH contributed to the study conception and design. JW, DL, XT, YZ, XC, XHC, XMC coordinated the clinical trial and the collection of samples and clinical data. ZW and JW performed statistical analyses on clinical and metabolomic data. QH and RH supervised the study. ZW, JW and RH wrote the draft of the manuscript with additional content provided by all authors. All authors approved the final manuscript.

Corresponding author

Correspondence to Runyue Huang.

Ethics declarations

Ethics approval and consent to participate

All participants provided written informed consent and the protocol was approved by the Medical Ethics Committee of the Second Affiliated Hospital of Guangzhou University of Chinese Medicine (B2016-073-01).

Consent for publication

Not applicable.

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.

Supplementary information

Additional file 1.

Patient clinical outcomes for Biqi and LEF treatment arms in the PP analysis.

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 The Creative Commons Public Domain Dedication waiver ( 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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, Z., Wu, J., Li, D. et al. Traditional Chinese medicine Biqi capsule compared with leflunomide in combination with methotrexate in patients with rheumatoid arthritis: a randomized controlled trial. Chin Med 15, 36 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: