Serum metabolomics analysis of deficiency pattern and excess pattern in patients with rheumatoid arthritis
Chinese Medicine volume 17, Article number: 71 (2022)
Rheumatoid arthritis (RA) is a chronic and refractory autoimmune disease. Deficiency pattern (DP) and excess pattern (EP), as crucial types of Chinese medicine pattern diagnoses published by International Classification of Diseases 11th Revision (ICD-11), could provide new strategies for RA diagnosis. However, the biological basis of DP and EP of RA is not explicit.
19 female RA DP patients, 41 female RA EP patients and 30 female healthy participants were included in the study. The serums of participants were collected and analyzed by metabolomics based on ultra-performance liquid chromatography coupled with quadrupole time-of-flight mass spectrometry to profile metabolic characteristics of RA DP and EP. Furthermore, bioinformatics analysis results were obtained by using Ingenuity Pathway Analysis (IPA) and statistical analysis was performed by SAS version 9.4 for further identification of potential biomarkers.
Serum metabolic profiling revealed 25 and 24 differential metabolites in RA DP and EP respectively, and 19 metabolites were common to RA DP and EP. Compared with DP group, L-Homocysteic acid, LysoPE(P-16:0/0:0), N(omega)-Hydroxyarginine and LysoPC(16:0/0:0) decreased (P < 0.05), and Pyruvic acid, D-Ribose, Gamma-Glutamylserine, PE(22:0/24:1(15Z)), Inosinic acid increased (P < 0.05) in EP group. Menawhile, S-Nitrosoglutathione, 5-Thymidylic acid, SN38 glucuronide, PE(22:0/24:0), PC(24:0/24:1(15Z)) and Bisdiphosphoinositol tetrakisphosphate increased significantly in DP group compared to EP group (P < 0.05). For the unique metabolites, bioinformatics analysis results showed that 5-Methoxytryptamine involved in Melatonin Degradation II and Superpathway of Melatonin Degradation is the key metabolite to RA DP. Meanwhile, GABA is the key metabolite in EP group, which involved in Glutamate Dependent Acid Resistance, GABA Receptor Signaling, Glutamate Degradation III (via 4-aminobutyrate) and 4-aminobutyrate Degradation I. Bioinformatics analysis between unique metabolites of RA DP and EP groups with human target genes for RA showed that 5-methoxytryptamine and LysoPC(18:1(9Z)/0:0), the unique metabolites of RA DP, might participate in colorectal cancer metastasis signaling, tumor microenvironment pathway, apoptosis signaling, MYC mediated apoptosis signaling, erythropoietin signaling pathway and LXR/RXR activation. Simultaneously, GABA, LysoPA(18:1(9Z)/0:0) and L-Targinine, the unique metabolites of RA EP, might participate in neuroinflammation signaling pathway, osteoarthritis pathway, glucocorticoid receptor signaling, ILK signaling, IL-17 signaling and HIF1α signaling.
The study indicates that serum metabolomics preliminarily revealed the biological basis of RA DP and EP. 5-methoxytryptamine, LysoPC(18:1(9Z)/0:0) and GABA, LysoPA(18:1(9Z)/0:0), L-Targinine might be the predictors to distinguish the DP and EP of RA respectively. These interesting results provide thoughts for further study of traditional medicine patterns of ICD-11. It also contributes to provide strategy for personalized precision treatment of RA and further validation is needed.
Rheumatoid arthritis (RA) is a heterogeneous and refractory disease, which primarily involves the joints and cause disability . The etiology of RA is thought to be associated with genetic and epigenetic components, environmental factors and autoimmunity . Accordingly, new strategies and therapies are urgently required for the reason that many patients still do not respond to current therapies . Chinese medicine (CM) is a system that uses an individualized approach based on the unique combination of symptoms and signs of each patient over 3000 years [3, 4]. CM pattern diagnosis is an approach that takes into account a broad spectrum of symptoms and signs and will lead more scientific new strategies and therapies for RA [3, 4].
Deficiency pattern (DP) and excess pattern (EP) are two crucial types of CM pattern diagnoses, which belongs to the principle-based patterns in the category of traditional medicine patterns published by International Classification of Diseases 11th Revision (ICD-11) [5,6,7]. Pattern differentiation in CM can help define specific indications for biomedical therapy in the disease treatment . Based on the traditional medicine theory, principle-based patterns are an approach that take into account a wide range of signs and symptoms . DP is manifested as pale or sallow complexion, bradypsychia, fatigue, shortness of breath, loose stools, uracratia, tongue with little or no coating, weak pulse, and others. Whereas EP is shown in red complexion, dysphoria, delirium, stentorian and raucous, constipation, dysuria, tongue coating thick and greasy, forceful pulse, and so on . In CM clinical practice, patients with the same disease can be divided into different groups according to traditional medicine patterns . In the near future, CM pattern diagnoses will facilitate further individualized medicine.
Metabolomics, the analysis of concentration profiles of low molecular weight metabolites present in biological fluids, has immense potential as a novel diagnostic tool in complex disorders such as RA, systemic lupus erythematosus [11, 12]. Simultaneously, metabolomics has highly discriminatory for biological perturbations or disease states in identifying new biomarkers . Thus, metabolomics can be used to identify biological characteristics of CM pattern, which is an important approach to distinguish disease states. The urine metabolomics of diabetes mellitus (DM) indicated that xylose and C4 sugar 2 were higher in the EP than in the DP. These potential biomarkers reflected the deregulation of glucose metabolism in diabetic individuals that might be helpful for the diagnosis and therapy of DM . Xylopyranoside for the diagnosis of DP and ribonic acid, uric acid, d-Ribose, and cyclohexanone for the diagnosis of EP in patients with chronic hepatitis B by urinary metabolomics. Meanwhile, the immune function of EP was demonstrated higher than that of DP . The difference could help us transform the concept of DP and EP to modern therapeutic approaches and prompts us to further precision medicine.
Previous studies have studied the biological basis of cold pattern and heat pattern of RA by using plasma metabolomics [15, 16]. Meanwhile, the characteristic gene expression profile of DP in RA patients was explored by bioinformatics analysis approach , and the DP-associated autoantibodies was identified by using protein chips . The above findings established the foundation to study CM pattern of RA. Therefore, this study will use serum metabolomics to explore the biological characteristics of the DP and EP in RA.
Materials and methods
The participants including female RA DP patients, female RA EP patients and healthy females were recruited from November 2018 to June 2019 in The First Affiliated Hospital of Henan University of CM. The sample size was based on previous studies [14, 19, 27]. The inclusion criteria include the following: (1) Patients with diagnosed of the RA and the diagnostic criteria followed 1987 or 2010 RA classification criteria [20, 21]; (2) DP and EP were classified according to the guidelines for the diagnosis and treatment of TCM patterns in RA , manifestations and other diagnostic information were determined independently by two experienced traditional Chinese physicians to ensure an objective evaluation; (3) The participants ranged in age from 18 to 80; (4) Participants signed the informed consent. Exclusion criteria: (1) RA patients with severe diseases requiring glucocorticoids; (2) RA patients who have been taking drugs for a long time or have not stopped taking drugs within 1 week; (3) Patients with circulatory, respiratory, digestive, urogenital, hematopoietic system, central nervous system and other serious diseases; (4) Pregnant or breast-feeding patients; (5) Patients with unclear prior medication records; (6) RA patients with other TCM patterns. Meanwhile, this study was approved by the Ethics Committee at the Institute of Basic Research in Clinical Medicine, China Academy of Chinese Medical Sciences (No.201,608).
Each blood sample collected in a fasting condition was left standing at room temperature for 30 min. Then blood samples were centrifuged for 20 min at 3000 g at 4 °C to obtain the serum and transferred into clean tubes. All serum samples were stored at − 80 °C until LC-MS analysis.
Chromatographic separation was performed on a 1.8 μm, 2.1 mm × 100 mm ACQUITY UPLC®HSS T3 column (Waters, USA) using an ACQUITY Ultra Performance LC system (Waters corp., Milford, MA, USA) equipped with a binary solvent delivery system, an auto-sampler, and high temperature column oven. The column was maintained at 30 °C. The flow rate was set at 0.4 ml/min. The sample injection volume was 1.2 µl. Solvent A was water mixed with 0.1% formic acid, and solvent B was acetonitrile mixed with 0.1% formic acid. The gradient elution programmes were: 0–0.5 min 5% B, 0.5–7 min 100% B, 7–8.6 min 100% B, 8.6–11.5 min 100% B, 11.5–11.6 min 5% B, 11.6–16 min 5% B.
Mass spectrometry detection was acquired on a quadrupole time-of-flight mass spectrometer (Waters Corp., Milford, MA, USA) equipped with an electrospray ion source. The positive ion mode detection: The desolvation gas was set at 600 L/h at a temperature of 350 °C. The cone gas was set at 37 L/h. The source temperature was set at 100 °C. The capillary and cone voltage were optimized at 1.5 kV and 40 V, respectively. Full scan spectra were performed between 50 and 1000 m/z with a 0.2 s scan time and a 0.1 s interscan delay. The negative ion mode detection: The desolvation gas was set at 800 L/h at a temperature of 500 °C. The cone gas was set at 50 L/h. The source temperature was set at 120 °C. The capillary and cone voltage were optimized at 2 kV and 40 V, respectively. Full scan spectra were performed between 50 and 1200 m/z with a 0.2 s scan time and a 0.1 s interscan delay.
Data pre-processing, multivariate statistical analysis and metabolite identification
The raw serum LC-MS data was pre-processed using Waters Progenesis QI 2.0 software (Nonlinear Dynamics, Newcastle, U.K.). Progenesis QI can provide an accurate measurement of the compounds in the bio-samples and facilitate flexible comparisons across runs. It included steps of importing data, reviewing alignment, experiment design setup, picking peaks, identifying and reviewing compounds and performing compound statistical analysis. Then, the data were exported into SIMCA-P software (v13.0) (Umetrics, Malmö, Sweden) for multivariable statistical analysis (MVA). The MVA of principal component (PCA), partial least squares discriminant (PLS-DA) and orthogonal partial least square-discriminant (OPLS-DA) models were both applied to observe the classifications for different groups in score plots. Next, to find significant metabolites contributing to the classifications, ions with variable importance in the projection (VIP) > 1 were picked out from the OPLS-DA loading plots. Then univariate analysis of t-test and fold-change were also applied to these ions. Ions with P-values < 0.05, fold-changes > 2.0 (or < 0.5) and relative standard deviation (RSD) < 0.3 were finally regarded as differentiated metabolite ions. Then they were structurally identified and interpreted based on searches of their accurate masses in HMDB (www.hmdb.ca) database. Lastly, the isotopic distribution, retention time and fragments of commercial standards were further confirmed for the metabolites of interest.
The analysis of canonical pathways, bio-functions and the networks were conducted by using the Ingenuity Pathway Analysis system (IPA, http://www.ingenuity.com) for the candidate metabolites and proteins to gain further insights into the typical metabolic alterations.
Continuous variables were described using mean ± standard deviation. The count data were expressed as the number of cases or percentages. Differences between groups were assessed using analysis of t-tests for continuous variables and the Chi-square test or Fisher’s exact test for categorical variables. SAS version 9.4 was used for all statistical analyses. P value was set as 0.05 for statistical significance.
Clinical characteristics of participants
A total of 90 participants including 19 female RA DP patients, 41 female RA EP patients (Fig. 1) and 30 healthy females were recruited. The characteristics of the enrolled subjects, including age, duration of RA, body mass index (BMI), white blood cell (WBC), red blood cell (RBC), hemoglobin (HGB), platelet (PLT), albumin (ALB), erythrocyte sedimentation rate (ESR), C-reactive protein (CRP), rheumatoid factor (RF)-IgG, RF-IgA, RF-IgM, smoked, hypertension, diabetes and coronary heart disease were shown in Table 1. In terms of age, there was no statistical difference between DP, EP and healthy participants, respectively (P > 0.05), but there was statistical difference between DP and EP (P = 0.003). There were no significant differences in BMI, WBC, RBC, HGB, PLT, ALB, between the three groups (P > 0.05). In terms of duration of RA, ESR, CRP, RF-IgG, RF-IgA, RF-IgM, smoked, hypertension, diabetes and coronary heart disease, there were no statistically significant differences between DP and EP (P > 0.05).
PCA and PLS-DA analysis of metabolomics profiles in DP, EP, and healthy groups
PCA was used to determine the presence of inherent similarities in spectral profiles and the corresponding PLS-DA analysis was used to identify discriminating metabolites and differentiate the three groups. Although there was some overlap between DP and EP group, the PCA (Fig. 2A, B) and PLS-DA (Fig. 2C, D) results revealed a clear separation between the serum samples of participants from the DP, EP and healthy groups which demonstrating specific metabolic profiles of different groups. In PLS-DA positive ion mode: R2X = 0.471 (cum), R2Y = 0.862 (cum), Q2 = 0.739(cum). And in PLS-DA negative ion mode: R2X = 0.323 (cum), R2Y = 0.808 (cum), Q2 = 0.772 (cum).
PCA and OPLS-DA analysis of metabolomics profiles in DP and healthy groups
PCA (Fig. 3A, B) and OPLS-DA (Fig. 3C, D) analysis results showed that DP group was clearly separated from those in the healthy group, indicating that the RA DP patients have totally different metabolic serum characteristics from healthy individuals. In OPLS-DA positive ion mode: R2X = 0.689 (cum), R2Y = 0.992 (cum), Q2 = 0.972 (cum). And in OPLS-DA negative ion mode: R2X = 0.446 (cum), R2Y = 0.99 (cum), Q2 = 0.972 (cum). By metabolite identification, 25 differential metabolites including 12 in positive ion mode and 13 in negative ion mode were obtained between DP and healthy groups (Table 2). Compared with the healthy group, 10 metabolites in DP group showed a downward trend, and 15 metabolites showed an upward trend (P < 0.01).
PCA and OPLS-DA analysis of metabolomics profiles in EP and healthy groups
PCA (Fig. 4A, B) and OPLS-DA (Fig. 4C, D) analysis results showed that EP group was clearly separated from those in the healthy group, indicating that the RA EP patients have totally different metabolic serum characteristics from healthy individuals. In OPLS-DA positive ion mode: R2X = 0.667 (cum), R2Y = 0.985 (cum), Q2 = 0.968 (cum). And in OPLS-DA negative ion mode: R2X = 0.496 (cum), R2Y = 0.983 (cum), Q2 = 0.965 (cum). By metabolite identification, 24 differential metabolites including 11 in positive ion mode and 13 in negative ion mode were obtained between DP and healthy groups (Table 3). Compared with the healthy group, 7 metabolites in DP group showed a downward trend, and 17 metabolites showed an upward trend (P < 0.01).
Comparative analysis of common differential metabolites between DP and EP groups
19 metabolites were common to DP and EP groups, including 8 in positive ion mode: L-Homocysteic acid, S-Nitrosoglutathione, LysoPE(P-16:0/0:0), LysoPC(16:0/0:0), LysoPC(P-18:0/0:0), Bisdiphosphoinositol tetrakisphosphate, PE(22:0/24:0), PC(24:0/24:1(15Z)), and 11 in negative ion mode: Pyruvic acid, D-Ribose, N(omega)-Hydroxyarginine, Citric acid, Gamma-Glutamylserine, 5-Thymidylic acid, Inosinic acid, SN38 glucuronide, PC(O-22:1(13Z)/22:3(10Z,13Z,16Z)), PE(22:0/24:1(15Z)), Octanoyl-CoA.
Compared with DP group, L-Homocysteic acid (Fig. 5A), LysoPE(P-16:0/0:0) (Fig. 5B), N(omega)-Hydroxyarginine (Fig. 5C) and LysoPC(16:0/0:0) (Fig. 5D) decreased (P < 0.05), and Pyruvic acid (Fig. 5E), D-Ribose (Fig. 5F), Gamma-Glutamylserine (Fig. 5G), PE(22:0/24:1(15Z)) (Fig. 5H), Inosinic acid (Fig. 5I) increased (P < 0.05) in EP group. Menawhile, S-Nitrosoglutathione (Fig. 5J), 5-Thymidylic acid (Fig. 5K), SN38 glucuronide (Fig. 5L), PE(22:0/24:0) (Fig. 5M), PC(24:0/24:1(15Z)) (Fig. 5N) and Bisdiphosphoinositol tetrakisphosphate (Fig. 5O) increased significantly in DP group compared to EP group (P < 0.05). However, there was no statistical difference in Citric acid (Fig. 6A), Octanoyl-CoA (Fig. 6B), LysoPC(P-18:0/0:0) (Fig. 6C) and PC(O-22:1(13Z)/22:3(10Z,13Z,16Z)) (Fig. 6D) between DP and EP group (P > 0.05).
Bioinformatics analysis of unique metabolites of DP and EP groups
There were 6 unique metabolites in DP group: 5-Methoxytryptamine, 4-Fumarylacetoacetic acid, LysoPC(18:1(9Z)/0:0), LysoPC(20:1(11Z)/0:0), Glucosylceramide (d18:1/12:0), PE(20:3 (5Z,8Z,11Z) /22:6 (4Z,7Z,10Z,13Z,16Z,19Z)), and 5 unique metabolites in EP group: Gamma-Aminobutyric acid (GABA), LysoPA(18:1(9Z)/0:0), DG(18:3 (6Z,9Z,12Z)/22:6 (4Z,7Z,10Z,13Z,16Z,19Z) /0:0), L-Targinine, PC(24:0/24:0). IPA core function analysis results show that Melatonin Degradation II, Tyrosine Degradation I, Superpathway of Melatonin Degradation are biological pathways closely related to DP group (Fig. 7A). And Glutamate Dependent Acid Resistance, GABA Receptor Signaling, Glutamate Degradation III (via 4-aminobutyrate), 4-aminobutyrate Degradation I, Putrescine Degradation III, Circadian Rhythm Signaling, Neuroinflammation Signaling Pathway, Neurovascular Coupling Signaling Pathway are related biological pathways in EP group (Fig. 7B). In network analysis results, 5-Methoxytryptamine involved in Melatonin Degradation II and Superpathway of Melatonin Degradation is the key metabolite to DP group (Fig. 7C). Meanwhile, GABA is the key metabolite in EP group (Fig. 7D), which involved in Glutamate Dependent Acid Resistance, GABA Receptor Signaling, Glutamate Degradation III (via 4-aminobutyrate) and 4-aminobutyrate Degradation I.
Bioinformatics analysis between unique metabolites of RA DP and EP groups with human target genes for RA
GeneCards, the human gene database, was used to retrieve RA human target genes. Through the IPA analysis between unique metabolites of RA DP and EP groups with human target genes for RA, the network relationship was obtained (Fig. 8A). Among the 6 unique metabolites of RA DP, 5-methoxytryptamine and LysoPC(18:1(9Z)/0:0) were closely related to RA target genes. In particular, melatonin can link 5-methoxytryptamine to target genes in the extracellular space, plasma membrane, cytoplasm, and nucleus that are related to RA. Furthermore, 5-methoxytryptamine might play a role in cancer, tumor morphology, gastrointestinal disease, hepatic system disease, inflammatory disease, cardiovascular disease (Fig. 8D) etc. through colorectal cancer metastasis signaling, tumor microenvironment pathway, apoptosis signaling, MYC mediated apoptosis signaling, erythropoietin signaling pathway, LXR/RXR activation (Fig. 8B) etc. Meanwhile, GABA, LysoPA(18:1(9Z)/0:0) and L-Targinine, the unique metabolites of RA EP, were closely related to RA target genes. GABA, in particular, has been linked to RA-related target genes found in the extracellular space, plasma membrane, cytoplasm, and nucleus. Through neuroinflammation signaling pathway, osteoarthritis pathway, glucocorticoid receptor signaling, ILK signaling, IL-17 signaling and HIF1α signaling (Fig. 8C) etc. GABA might participate in inflammatory response, organismal injury and abnormalities, cancer, infectious diseases, skeletal and muscular disorders, neurological disease (Fig. 8E) etc.
CM is a medical system with at least 3000 years of uninterrupted clinical practice. With the release of ICD-11, traditional medicine, commonly used across many countries including China, Japan, and Korea, was brought within the perimeter of mainstream medicine to be considered right . From global extension of traditional medicine to the integration of traditional medicine with multiple Western medicine-based disciplines, people all over the world will benefit . Subgrouping of patients with complex diseases, such as RA, according to CM pattern diagnosis has the potential to provide opportunities for new strategies and better treatment outcomes [8, 25]. Moreover, metabolomics was used to explore the biological characteristics of DP and EP can provide scientific evidence for CM pattern diagnosis [5, 13, 14, 26].
RA is a heterogeneous and refractory inflammatory joint disease, which still require new strategies and therapies urgently . CM pattern differentiation based on clinical signs and symptoms showed a diverse range of biomolecules, proteins and genes from RA patients correlated well with cold pattern and heat pattern . Particularly, the cold pattern was related to Toll-like receptor signaling pathway, cell proliferation and the Jak-STAT cascade. The following related pathways in heat pattern were included: Calcium signaling pathway, cell adhesion molecules, PPAR signaling pathway, fatty acid metabolism and I-κB kinase/NF-κB cascade [27, 28]. Previous studies have shown that traditional Chinese medicine would be more effective in RA patients with DP, and western medicine would be more effective in RA patients with cold pattern. Based on different symptoms and signs, superior efficacy can be achieved in RA patients with different CM patterns . These results suggest that a better understanding of CM patterns might help to choose the most appropriate therapies .
In the present study, metabonomic technology was used to stratify the DP and EP in RA patients. Although the age of patients with DP and EP has a statistical difference, RA is a chronic autoimmune disease characterized by immunomodulatory disorder and a failure of spontaneous resolution of inflammation . Previous studies have shown that the human immune system deteriorates with age and immunosenescence is a series of age-related changes that affect the immune system . As an autoimmune disease, the immune function of patients with RA may also be affected by age. In CM clinical practice, the age of patients with DP was higher than that of patients with EP . This is consistent with our research. The results illustrate that substance fundaments and biological characteristics are separated clearly between DP and EP groups by metabolic profiles (Fig. 9). 19 differential metabolites in common were identified between DP and EP groups. There was no statistical difference between citric acid and octanoyl-COA in RA DP and EP groups. These two metabolites and their participation in the tricarboxylic acid cycle  are the potential biological characteristics of RA. Among the 15 statistical differences metabolites, L-Homocysteic acid and N(omega)-Hydroxyarginine decreased in EP is more obvious than that in DP. This suggests that amino acid metabolism is more significant in EP than in DP. Meanwhile, pyruvic acid, D-Ribose, inosinic acid and Gamma-Glutamylserine increased in EP is more obvious than that in DP. This reflects that energy metabolism is more significant in EP than in DP. In addition, the changes of lipid metabolites were different in DP and EP groups which indicates there was also lipid metabolism disorder in RA. These common metabolites have potential value for the identification of RA DP and EP.
Among the specific metabolites of RA DP, 5-methoxytryptamine as the intermediate of melatonin showed a decreasing trend. Melatonin secreted by the pineal gland widely exists in the body, which has the function of scavenging reactive oxygen species and active nitrogen and plays a certain role in hormone secretion, biological rhythm, anti-free radical damage and anti-aging [35, 36]. Previous study has found that the level of melatonin with Heart-spleen Deficiency Pattern, one type of CM DPs, is lower than that of normal people, reflecting that melatonin has a certain correlation with DP . In our study, the level of 5-Methoxytryptamine is lower than normal in RA DP. Simultaneously, Melatonin Degradation II and the Superpathway of Melatonin Degradation have been confirmed to be closely related to RA DP by bioinformatics analysis.
GABA, one of the specific metabolites of RA EP, showed a decreasing trend. GABA is the main inhibitory neurotransmitter in the brain, mainly involved in cell proliferation, differentiation and migration . GABA is mainly synthesized by glutamic acid via glutamic acid decarboxylase and acts through electrophilic GABA (A) receptors and/or hyper electrophilic GABA (B) receptors . In neurological diseases, GABA has a good anticonvulsant effect and can be used to fight the occurrence of epilepsy, by inhibiting the sudden abnormal firing of neurons in the brain, thereby normalizing brain dysfunction . The increase of GABA can alleviate the empirical pathological state, in other words, the decrease of GABA will aggravate the empirical pathological state in the body. The level of GABA decreased in RA EP, which reflected that the inhibition of neurotransmitters was reduced and bodily functions are more active than normal state in patients with RA EP.
The results of this study indicate that metabolomics can provide a basis for the classification of RA DP and EP. In other words, the classification of CM patterns can be distinguished by the difference of metabolites. The biological basis of the classification of RA DP and EP provides a new insight into the modern scientific connotation of RA, which will be helpful for the personalized treatment of RA patients. However, there are still some limitations in this study: (1) The sample size of this study was small, but representative RA DP and EP samples were still collected based on the rigorous diagnostic criteria, inclusion criteria and exclusion criteria, and meaningful results had been obtained; (2) The participants included in the study were all from Henan Province, which might reduce the influence of environmental factors on traditional medicine patterns of RA to some extent. However, the results could also reflect the traditional medicine patterns of RA in central China. Future studies under different regions and environments are still needed; (3) Since there is no animal and cell models of DP and EP, the study did not verify the findings based on cell biology experiments. However, metabolomics and bioinformatics analysis found some interesting results, which can provide a scientific basis for the biological basis of RA DP and EP and thoughts for further study of traditional medicine patterns of ICD-11; (4) DP and EP are the basic steps of the principle-based patterns for the diagnosis of RA. Subgroups can still be classified under each CM pattern diagnosis. For example, liver and kidney deficiency syndrome is one of the important subgroups of DP, while cold and dampness syndrome is one of the important subgroups of EP, which needs to be further studied.
In conclusion, serum metabolomics preliminarily revealed the biological basis of RA DP and EP. 5-methoxytryptamine, LysoPC(18:1(9Z)/0:0) and GABA, LysoPA(18:1(9Z)/0:0), L-Targinine might be the predictors to distinguish the DP and EP of RA respectively, which provides a potential basis for the understanding of CM pattern diagnoses. These interesting results provide thoughts for further study of traditional medicine patterns of ICD-11. It also contributes to provide strategy for personalized precision treatment of RA and further validation is needed.
The datasets generated and/or analyzed during the current study are available from the corresponding author upon reasonable request.
International Classification of Diseases 11th Revision
Multivariable statistical analysis
Partial least squares discriminant
Orthogonal partial least square-discriminan.
Variable importance in the projection
Relative standard deviation
Ingenuity Pathway Analysis system
Body mass index
White blood cell
Red blood cell
Erythrocyte sedimentation rate
Smolen JS, Aletaha D, McInnes IB. Rheumatoid arthritis. Lancet. 2016;388(10055):2023–38.
Fan D, Xia Y, Lu C, Ye Q, Xi X, Wang Q, Wang Z, Wang C, Xiao C. Regulatory role of the RNA N6-methyladenosine modification in immunoregulatory cells and immune-related bone homeostasis associated with rheumatoid arthritis. Front cell Dev Biol. 2020;8:627893.
Lu A, Chen K. Chinese medicine pattern diagnosis could lead to innovation in medical sciences. Chin J Integr Med. 2011;17(11):811–17.
Seca S, Franconi G. Understanding Chinese medicine patterns of rheumatoid arthritis and related biomarkers. Medicines (Basel). 2018;5(1):17.
Maeda-Minami A, Yoshino T, Katayama K, Horiba Y, Hikiami H, Shimada Y, Namiki T, Tahara E, Minamizawa K, Muramatsu S, Yamaguchi R, Imoto S, Miyano S, Mima H, Mimura M, Nakamura T, Watanabe K. Discrimination of prediction models between cold-heat and deficiency-excess patterns. Complement Ther Med. 2020;49:102353.
Choi S. A proposed revision of the international classification of diseases, 11th revision, chap 26. Integr Cancer Ther. 2020;19:1475702865.
World HO. International Classification of Diseases 11th Revision-The global standard for diagnostic health information. 2019.
Lu C, Zha Q, Chang A, He Y, Lu A. Pattern differentiation in traditional Chinese medicine can help define specific indications for biomedical therapy in the treatment of rheumatoid arthritis. J Altern Complement Med. 2009;15(9):1021–5.
Zhang Y, Li JT, Liu T, Zhao KP, Zhang T. Identification of deficiency pattern and excess pattern between Shang Han Lun and principle-based patterns. Traditional Chinese Medicinal Research. 2021;34(7):13–5.
Jiang M, Lu C, Zhang C, Yang J, Tan Y, Lu A, Chan K. Syndrome differentiation in modern research of traditional Chinese medicine. J Ethnopharmacol. 2012;140(3):634–42.
Li J, Che N, Xu L, Zhang Q, Wang Q, Tan W, Zhang M. LC-MS-based serum metabolomics reveals a distinctive signature in patients with rheumatoid arthritis. Clin Rheumatol. 2018;37(6):1493–502.
Scrivo R, Casadei L, Valerio M, Priori R, Valesini G, Manetti C. Metabolomics approach in allergic and rheumatic diseases. Curr Allergy Asthma R. 2014;14(6):445.
Wu T, Yang M, Wei HF, He SH, Wang SC, Ji G. Application of metabolomics in traditional chinese medicine differentiation of deficiency and excess syndromes in patients with diabetes mellitus. Evid Based Complement Alternat Med. 2012;2012:968083.
Sun S, Dai J, Fang J, Gou X, Cao H, Zheng N, Wang Y, Zhang W, Zhang Y, Jia W, Hu Y. Differences of excess and deficiency zheng in patients with chronic hepatitis B by urinary metabonomics. Evid Based Complement Altern Med. 2013;2013:738245.
Gu Y, Lu C, Zha Q, Kong H, Lu X, Lu A, Xu G. Plasma metabonomics study of rheumatoid arthritis and its Chinese medicine subtypes by using liquid chromatography and gas chromatography coupled with mass spectrometry. Mol Biosyst. 2012;8(5):1535.
Guo H, Niu X, Gu Y, Lu C, Xiao C, Yue K, Zhang G, Pan X, Jiang M, Tan Y, Kong H, Liu Z, Xu G, Lu A. Differential Amino Acid, carbohydrate and lipid metabolism perpetuations involved in a subtype of rheumatoid arthritis with Chinese medicine cold pattern. Int J Mol Sci. 2016;17(10):1757.
Wang M, Chen G, Lu C, Xiao C, Li L, Niu X, He X, Jiang M, Lu A. Rheumatoid arthritis with deficiency pattern in traditional chinese medicine shows correlation with cold and hot patterns in gene expression profiles. Evid Based Complement Altern Med. 2013;2013:248650.
Zhao H, Zhang Y, Liu B, Li L, Zhang L, Bao M, Guo H, Xu H, Feng H, Xiao L, Yi W, Yi J, Chen P, Lu C, Lu A. Identification of characteristic autoantibodies associated with deficiency pattern in traditional Chinese medicine of rheumatoid arthritis using protein chips. Front Pharmacol. 2019;10:755.
Zhao C, Chen L, Cai H, Yao W, Zhou Q, Zhu H, Gao Y, Liu P, Gou X, Zhang H. Classification of Gan Dan Shi Re pattern and Gan Shen Yin Xu pattern in patients with hepatitis B cirrhosis using metabonomics. Evid Based Complement Altern. 2018;2018:1–13.
Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, Cooper NS, Healey LA, Kaplan SR, Liang MH, Luthra HS, Et A. The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum. 1988;31(3):315–24.
Aletaha D, Neogi T, Silman AJ, Funovits J, Felson DT, Bingham CO, Birnbaum NS, Burmester GR, Bykerk VP, Cohen MD, Combe B, Costenbader KH, Dougados M, Emery P, Ferraccioli G, Hazes JMW, Hobbs K, Huizinga TWJ, Kavanaugh A, Kay J, Kvien TK, Laing T, Mease P, Ménard HA, Moreland LW, Naden RL, Pincus T, Smolen JS, Stanislawska-Biernat E, Symmons D, Tak PP, Upchurch KS, Vencovský J, Wolfe F, Hawker G. 2010 Rheumatoid arthritis classification criteria: an American College of Rheumatology/European League Against Rheumatism collaborative initiative. Arthritis Rhuematol. 2010;62(9):2569–81.
Jiang Q, Wang H, Gong X, Luo C. Guidelines for the diagnosis and treatment of TCM syndromes in rheumatoid arthritis. J Tradit Chin Med. 2018;20(59):1794–800.
The L. ICD-11. Lancet. 2019;393(10188):2275.
Lam WC, Lyu A, Bian Z. ICD-11: impact on traditional Chinese medicine and world healthcare systems. Pharm Med. 2019;33(5):373–7.
van Wietmarschen H, Yuan K, Lu C, Gao P, Wang J, Xiao C, Yan X, Wang M, Schroen J, Lu A, Xu G, van der Greef J. Systems biology guided by Chinese medicine reveals new markers for sub-typing rheumatoid arthritis patients. J Clin Rheumatol. 2009;15(7):330–37.
Jo J, Lee HL, Lee SH, Baek Y, Lee Y, Kim H. Exploratory analysis of cold, heat, deficiency, or excess pattern distribution in women with dysmenorrhea. Explore (NY). 2020;16(3):185–88.
Lu C, Xiao C, Chen G, Jiang M, Zha Q, Yan X, Kong W, Lu A. Cold and heat pattern of rheumatoid arthritis in traditional Chinese medicine: distinct molecular signatures indentified by microarray expression profiles in CD4-positive T cell. Rheumatol Int. 2012;32(1):61–8.
Lu C, Niu X, Xiao C, Chen G, Zha Q, Guo H, Jiang M, Lu A. Network-based gene expression biomarkers for cold and heat patterns of rheumatoid arthritis in traditional Chinese medicine. Evid-Based Compl Alt. 2012;2012:1–17.
He Y, Lu A, Lu C, Zha Y, Yan X, Song Y, Zeng S, Liu W, Zhu W, Su L, Feng X, Qian X, Tsang I. Symptom combinations assessed in traditional Chinese medicine and its predictive role in ACR20 efficacy response in rheumatoid arthritis. Am J Chin Med. 2008;36(4):675–83.
Jiang M, Zha Q, Zhang C, Lu C, Yan X, Zhu W, Liu W, Tu S, Hou L, Wang C, Zhang W, Liang Q, Fan B, Yu J, Zhang W, Liu X, Yang J, He X, Li L, Niu X, Liu Y, Guo H, He B, Zhang G, Bian Z, Lu A. Predicting and verifying outcome of Tripterygium wilfordii Hook F. based therapy in rheumatoid arthritis: from open to double-blinded randomized trial. Sci Rep. 2015;5:9700.
Chen Z, Bozec A, Ramming A, Schett G. Anti-inflammatory and immune-regulatory cytokines in rheumatoid arthritis. Nat Rev Rheumatol. 2019;15(1):9–17.
Nikolich-Žugich J. The twilight of immunity: emerging concepts in aging of the immune system. Nat Immunol. 2018;19(1):10–9.
Zhu Y, Pan Z, Wei Y, Du J, Wang C, Shuai X, Zhang M. Study of the relationship among TCM excess and defi ciency syndromes, heart rate, HRV in the hypertension population. China J Traditional Chin Med Pharm. 2013;28(12):3671–74.
Zheng T, Su S, Dai X, Zhang L, Duan J, Ou-Yang Z. Metabolomic analysis of biochemical changes in the serum and urine of freund’s adjuvant-induced arthritis in rats after treatment with silkworm excrement. Molecules. 2018;23(6):1490.
Maffei ME. 5-Hydroxytryptophan (5-HTP): natural occurrence, analysis, biosynthesis, biotechnology, physiology and toxicology. Int J Mol Sci. 2021;22(1):181.
Baburina Y, Lomovsky A, Krestinina O. Melatonin as a potential multitherapeutic agent. J Person Med. 2021;11(4):274.
Shi J, Chai R, Huang Y, Jiang T, Zhu Y, Zhang Z. Effect of Yishen Ningxin Formula (益神宁心方) on melatonin and cortisol levels in saliva of primary insomnia patients with heart-spleen deficiency pattern. J Tradit Chin Med. 2021;62(24):2157–63.
Tang X, Jaenisch R, Sur M. The role of GABAergic signalling in neurodevelopmental disorders. Nat Rev Neurosci. 2021;22(5):290–307.
Watanabe M, Maemura K, Oki K, Shiraishi N, Shibayama Y, Katsu K. Gamma-aminobutyric acid (GABA) and cell proliferation: focus on cancer cells. Histol Histopathol. 2006;21(10):1135–41.
Castro-Torres RD, Ureña-Guerrero ME, Morales-Chacón LM, Lorigados-Pedre L, Estupiñan-Díaz B, Rocha L, Orozco-Suárez S, Rivera-Cervantes MC, Alonso-Vanegas M, Beas-Zárate C. New Aspects of VEGF, GABA, and glutamate signaling in the neocortex of human temporal lobe pharmacoresistant epilepsy revealed by RT-qPCR arrays. J Mol Neurosci. 2020;70(6):916–29.
The study was supported by the Fundamental Research Funds for the Central public welfare research institutes (ZZ11-056), the National Natural Science Foundation of China (82074269) and the Innovation Team and Talents Cultivation Program of National Administration of Traditional Chinese Medicine (ZYYCXTD-D-202005).
Ethics approval and consent to participate
The informed consent for inclusion was taken consciously from all the participants included in this study. The study was conducted under the Declaration of Helsinki, and the protocol was approved by the Ethics Committee of the Institute of Basic Research in Clinical Medicine (No. 201608).
Consent for publication
These authors declare that there are no conflicts of interest regarding this work.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Liu, B., Guo, H., Li, L. et al. Serum metabolomics analysis of deficiency pattern and excess pattern in patients with rheumatoid arthritis. Chin Med 17, 71 (2022). https://doi.org/10.1186/s13020-022-00632-5