Skip to main content
  • Research article
  • Open access
  • Published:

Metabolomic alterations associated with Behçet’s disease



The diagnosis of Behçet’s disease (BD) remains challenging due to the lack of diagnostic biomarkers. This study aims to identify potential serum metabolites associated with BD and its disease activity.


Medical records and serum samples of 24 pretreated BD patients, 12 post-treated BD patients, and age-matched healthy controls (HC) were collected for metabolomics and lipidomics profiling using UPLC-QTOF-MS and UPLC-QTOF-MSE approaches. Additionally, serum samples from an independent cohort of BD patients, disease controls including rheumatoid arthritis (RA), systemic lupus erythematosus (SLE), Takayasu’s arteritis (TA), Crohn’s disease (CD) patients, and HC were collected for further validation of two potential biomarkers using UPLC-QTOFMS analysis.


Unsupervised principal component analysis (PCA) showed a clear separation of metabolomics profiles of BD patients from HC. Statistical analysis of the data revealed differential metabolites between BD patients and HC. The serum levels of some phosphatidylcholines (PCs) were found to be significantly lower in BD patients, while the levels of several polyunsaturated fatty acids (PUFAs) were increased markedly in the BD group compared with HC. Furthermore, the serum level of two omega-6 PUFAs, linoleic acid (LA) and arachidonic acid (AA), were dramatically decreased in patients with remission. A validation cohort confirmed that the serum LA and AA levels in BD patients were significantly higher than those in HC and patients with RA, SLE, TA, and CD. In addition, receiver operating characteristic (ROC) analysis indicated good sensitivity and specificity.


The serum metabolomics profiles in BD patients are altered. Serum LA and AA are promising diagnostic biomarkers for BD.


Behçet’s disease (BD) is a chronic multisystem inflammatory disorder characterized by recurrent oral and genital ulceration, uveitis, and skin lesions. Currently, the diagnosis of BD is primarily based on clinical manifestations, and no diagnostic biomarkers are available. BD is a multisystem vasculitis, so both arteries and veins of all sizes may be involved [1]. Since the clinical symptoms of BD are variable, it is sometimes difficult to distinguish it from other diseases such as inflammatory bowel disease and Reiter’s syndrome. As a result, early diagnosis remains a challenge in clinical practice.

Metabolomics, an emerging “omics” science, uses state-of-the-art quantitative analysis approaches and advanced bioinformatic methods to characterize the metabolome. It reflects both physiological and pathological states, and it may detect the alterations of affected metabolites at the early stages of disease due to its great sensitivity [2]. Metabolomic methods have been used for evaluating clinical diagnosis and therapeutic treatment in a variety of diseases, such as cancer, diabetes, multiple sclerosis, primary biliary cirrhosis, and autoimmune hepatitis [3,4,5,6,7]. Metabolic abnormalities in BD remain elusive. Given that serum is an accessible and informative biofluid, this study aims to identify serum metabolites in BD and to elucidate the metabolites responsive to treatment using a metabolomics approach.


Patients and controls

For metabolomics and lipidomics profiling, 24 BD patients and 25 gender- and age-matched healthy controls (HC) (without a personal or family history of autoimmune diseases) were enrolled from Peking Union Medical College Hospital (PUMCH) between March 2014 and November 2014. For further validation, an independent cohort of BD (n = 25), rheumatoid arthritis (RA) (n = 12), systemic lupus erythematosus (SLE) (n = 12), Takayasu’s arteritis (TA) (n = 15), and Crohn’s disease (CD) (n = 15) patients, and 19 HC were enrolled from March 2014 to July 2018. All BD patients fulfilled the 1990 International Study Group BD criteria or the new International Criteria for Behçet’s Disease (ICBD) [8, 9]. RA, SLE, TA and CD patients fulfilled their respective diagnostic and classification criteria [10,11,12,13]. All participants underwent a clinical evaluation, and hospital records were reviewed. The following data were collected: disease duration, clinical manifestations, erythrocyte sedimentation rate (ESR)/C-reactive protein (CRP) level, and treatment. BD disease activities were evaluated according to the BD Current Activity Form 2006 (BDCAF 2006;

This study was carried out in accordance with the recommendations of the institutional committee for the Protection of Human Subjects from PUMCH. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocol was approved by the institutional committee for the Protection of Human Subjects from PUMCH. All methods were performed in accordance with the relevant guidelines and regulations.

Sample preparation for metabolomics/lipidomics profiling

Sterile siliconized 0.6-mL Eppendorf tubes were used for sample preparation, and 25 μL serum was added to the tubes followed by 100 μL cold chloroform/methanol (2/1) containing lipid standards at predetermined concentrations as described previously [14]. The mixture was vortexed for 30 s at room temperature and then centrifuged at 13,000×g for 5 min to separate the polar and nonpolar species. Upper and lower phases were collected separately and transferred to new tubes. The white interphase was discarded. The collected samples were dried using a speed vacuum. The pellet of the upper phase, which primarily contained polar metabolites, was resuspended in 200 μL 50% methanol for metabolomic profiling. The lower phase was resuspended in 200 μL isopropanol/acetonitrile/water (50/25/25) for lipidomic analysis.


The chromatographic and mass spectrometric parameters were used as described previously [15]. The ultra-performance liquid chromatography (UPLC) column eluent was introduced directly into the mass spectrometer by electrospray. For metabolomic profiling, a 2-μL sample was injected onto a reverse-phase ACQUITY BEH C18 50 × 2.1 mm 1.7-μm column (Waters Corp., Milford, MA) using an ACQUITY UPLC system (Waters Corp., Milford, MA). For lipidomics profiling, an Acquity CSH C18 50 × 2.1 mm 1.7-μm column (Waters Corp., Milford, MA) was used in the UPLC-quadrupole time-of-flight mass spectrometry (UPLC-QTOF-MSE) analysis. MSE is a technique by which both the precursor and fragment mass spectra are acquired by alternating between high and low collision energy during a single chromatographic run. Mass spectrometric analysis was performed on a XEVO G2 QTOF (Waters) operating in both positive and negative modes. Accurate mass was maintained by introducing the LockSpray interface of sulfadimethoxine (311.0814 [M + H]+ or 309.0658 [M − H]−) at a concentration of 250 pg/μL in 50% aqueous ACN at a rate of 150 μL/min. For biomarker validation, a 2-μL sample was injected into a reverse-phase ACQUITY HSS T3 C18 100 × 2.1 mm 1.7-μm column (Waters Corp., Milford, MA) and analyzed using the consistent UPLC-QTOF-MS system. The mobile phase consisted of acetonitrile (A) and water containing 0.1% (v/v) formic acid (B), while the gradient elution program (0.0–18.0 min, 10%–95% A; 18.1–20.0 min, 100% A) was applied for favorable separation. The flow rate was set at 0.4 mL/min. The column temperature was 40 °C. All chromatograms and mass spectrometric data were acquired in centroid mode using the MassLynx software (Waters Corp., Milford, MA).

Data processing and multivariate data analysis

Raw mass spectrometric data were processed using Progenesis QI software (Nonlinear Dynamics, Durham, NC) to generate a data matrix that consisted of the retention time, m/z value, and the normalized peak area. Statistical analysis and putative ion identification on the postprocessed data were conducted using MetaboLyzer [16]. Statistically significant ions were putatively identified in MetaboLyzer, which utilizes the Human Metabolome Database (HMDB), LipidMaps, and the Kyoto Encyclopaedia of Genes and Genomes (KEGG) database [8] while accounting for possible adducts, H+, Na+, and NH4+ in the ESI+ mode, and H and Cl in the ESI mode. The m/z values were compared with the exact mass of small molecules in the databases, from which putative metabolites were identified with a mass error of 10 ppm or less. KEGG annotated pathways associated with these putative metabolites were also identified. Lipid ions were validated with the fragmentation in MSE results based on their identifier fragments and retention time with the help of nonendogenous lipid standards from each lipid class and/or the comparison of tandem mass spectrometry (MS/MS) fragments with reference spectra provided in METLIN and LipidMaps databases. SIMCA-P+ (Umetrics, Umea, Sweden) was used for principal component analysis (PCA). The heat map displaying the relative levels of differential metabolites was generated by the Random Forests (RF) algorithm as explained in detail in previous studies [17]. MS data acquired in negative ion mode were employed for quantitation of arachidonic acid (AA) and linoleic acid (LA).

Sample preparation for biomarker verification

For the preparation of calibration samples, reference substances of rosmarinic acid, AA, and LA were purchased from Sigma-Aldrich Company (MO, USA). Rosmarinic acid was dissolved in methanol to produce the internal standard (IS) solution at the concentration of 1 μg/mL. The reference solutions of the two targeted metabolites, AA and LA, at the stock concentrations of 1/1000 (v/v) for each were serially diluted with the IS solution to produce a series of calibration standard solutions. All calibration standard solutions were sealed and stored at 4 °C until use.

For serum samples, 100 μL serum was added to the tube, followed by 400 μL cold methanol. The mixture was vortexed for 30 s at room temperature and then centrifuged at 13,000×g at 4 °C for 5 min to precipitate the protein. The supernatant was transferred to the Eppendorf tube, dried using a speed vacuum at 25 °C, and then dissolved in 100 μL of the IS solution for LC/MS analysis.

Statistical analysis

Experimental values are presented as mean ± SD. Statistical analysis was performed using GraphPad Prism (San Diego, CA). The significance of the metabolite level was determined using a two-tailed student t test. P values less than 0.05 were considered significant.


Baseline clinical characteristics of BD patients

24 BD patients (15 men and 9 women) were all of the Han Chinese population (100%). Their mean age was, on average, 35.83 ± 11.96 years old. The median disease duration of BD was 120 months (range 13–379 months). The median age at diagnosis of BD was 28.5 years (range 17–60 years). All cases initially presented with oral ulcers. The interval between the onset of an oral ulcer and the BD diagnosis ranged from 1 to 379 months, with a median time of 71 months. The baseline clinical characteristics and medications of the BD patients are shown in Table 1. Twelve cases had been reviewed with a mean follow-up time of 7.67 ± 2.06 months. After being treated with glucocorticoids or immunosuppressants, 11 patients (91.7%) improved as measured by decreased ESR and CRP levels.

Table 1 Demographics and clinical characteristics of Behçet’s disease (BD) patients

Serum metabolomics

To characterize metabolomic alterations associated with BD, we analyzed the metabolomics of BD patients and HC using a UPLC-QTOF-MS approach. Unsupervised PCA plots were generated by SIMCA-P software, and differential analysis was carried out using Metabolyzer. The volcano plot showing the differential ions between HC and BD patients is shown in Fig. 1. Individuals and variables with similar profiles are grouped together in the plot. In this PCA plot, the purple and green spots, which represent individual pretreated BD patients and HC, respectively, form two segregated clusters. Statistical analysis of the metabolomics data from pretreated BD patients and the control group revealed differential ions between these two groups with statistical significance as shown in the volcano plot (Fig. 1b). In the volcano plot, the red dots represent ions that show significantly different (p value less than 0.05) levels between BD patients and HC. Putative molecules of these differential ions were designated by screening the accurate mass in metabolite databases, as stated in the Methods.

Fig. 1
figure 1

Multivariate data analyses of serum metabolomic profiles segregate Behçet’s disease (BD) patients from the healthy cohort (HC; Nor). a Principle components analysis (PCA) unsupervised clustering plots for healthy (green triangle) and BD (purple triangle) are shown. b The volcano plot displays t test results of samples from BD and HC. Ions marked in red show a significant difference in intensity between BD and HC. c KEGG pathway enrichment analysis results of the differential ions identified in the multivariate analysis of serum metabolomic profiles from HC and BD patients. FDR false discovery rate

KEGG pathway analysis results indicate the metabolic pathways associated with the differential metabolites. Prominent pathways with a false discovery rate (FDR)-corrected p value less than 0.25 are shown in Fig. 1c.

Serum lipidomics

Since the major pathways in Fig. 1c pointed to lipid metabolism, we performed the lipidomic profiling analysis by UPLC-QTOF-MSE. To address whether the treatments have an effect on the disease-associated metabolomics, serum lipidomic profiling was performed using samples from twelve diagnosed patients before and after treatment, as well as a healthy cohort of the same number. Statistically significant differential metabolites between HC and pretreatment BD patients were determined using Metabolyzer. PCA analysis results for the healthy, pretreatment, and post-treatment groups based on these differential metabolites are shown in Fig. 2a. This PCA plot shows that the post-treatment cluster (red circle with a dotted line) is located in the middle, between the pretreatment cluster (dark blue circle with a dotted line) and the healthy cluster (green circle with a dotted line), suggesting a drift of diseased data points towards the direction of the healthy group after treatment. RF algorithms were used to generate a heat map of these differential metabolites for these three groups (Fig. 2b).

Fig. 2
figure 2

Serum lipidomic profile phenotype of treated Behçet’s disease (BD) patients moves towards that of the healthy cohort (HC). a PCA clustering plot of samples from these three groups with differential metabolites identified between pretreatment BD patients and HC. b Heat map of the metabolites with significantly different abundance between BD and HC

From the list of the differential ions, lower levels of several ions with the putative identification of phosphatidylcholine (PC) were found in the BD patient group compared with the HC group. Verification of the three selected ions was performed through MS/MS, which confirmed the identity of these ions as PCs (Additional file 1). The scatter plots in Fig. 3a show the lower level of the three PCs detected in the serum of the pretreatment BD patients compared with those in the serum of HC, but no significant difference was seen between the pretreatment and post-treatment groups.

Fig. 3
figure 3

Levels of phosphatidylcholines (PCs) and polyunsaturated fatty acids (PUFAs) are affected by treatments in a different way. Comparison of the level of three PCs (a) and two PUFAs (b) in healthy controls (HC), pretreatment Behçet’s disease (BD) (Pre-) patients, and post-treatment BD (Post-) patients. The significance of metabolite levels was determined using a two-tailed student t test. *p < 0.05, ***p < 0.001. ns not significant

In addition to the above phospholipids, levels of several free fatty acids from pretreatment BD patients were found to be significantly different from those of HC. We observed markedly lower levels of several polyunsaturated fatty acids (PUFAs) in the HC group compared with the pretreatment BD group, including two omega-6 (n-6) fatty acids (Fig. 3b) LA (18:2n-6) and AA (20:4n-6), and oleic acid (OA), an n-9 PUFA (Additional file 2). In contrast to the small difference between pretreatment and post-treatment BD groups for PCs, a significantly lower level of PUFAs was found in the post-treatment BD group compared with the pretreatment BD group, indicating that the treatment effectively corrected the abnormal increases of these PUFAs in BD patients. Validations of these PUFAs by MS/MS are shown in Additional file 2 and Additional file 3.

Receiver operating characteristic (ROC) analysis along with sensitivities and specificities of the area under the curve (AUC) > 0.85, is shown in Table 2 and in Additional file 4. ROC curves showed that AA was the most efficient diagnostic performance (AUC = 0.9495), compared with PCs and LA (Table 2 and Additional file 4). The sensitivity of PC(34:3), PC(40:8), LA, and AA in diagnosis of BD were comparable (0.96%, 0.88%, 0.9474%, and 0.9474%, respectively), but the specificity of PC(40:8) and AA were higher than PC(34:3) and LA (Table 2).

Table 2 The diagnostic values of phosphatidylcholines (PCs), arachidonic acid (AA), and linoleic acid (LA) in Behçet’s disease (BD)

To verify these findings, the concentrations of LA and AA were further determined using reference standards in an independent cohort containing BD, RA, SLE, TA, and CD patients and HC. As shown in Fig. 4 and Additional file 5, the serum levels of LA and AA in BD patients were significantly higher than those in HC (p = 2.35 × 10−3 and p = 6.1 × 10−6). ROC curves were also made to indicate their diagnostic efficiency (Table 2 and Additional file 4). The serum levels of LA and AA in patients with RA, SLE, and TA were comparable with those in HC (Fig. 4). In addition, the levels of LA and AA in BD patients were significantly higher than they were in disease controls, suggesting that LA and AA might serve as specific markers for BD. Intriguingly, the serum levels of LA and AA in patients with CD were significantly lower than those in HC (p = 0.016 and p = 0.002, respectively) which need further investigation.

Fig. 4
figure 4

Serum level of arachidonic acid (AA) and linoleic acid (LA) in the validation cohort. Levels of AA (a) and LA (b) in patients with Behçet’s disease (BD) (n = 25), systemic lupus erythematosus (SLE) (n = 12), rheumatoid arthritis (RA) (n = 12), Takayasu’s arteritis (TA) (n = 15), and Crohn’s disease (CD) (n = 15), and healthy controls (HC) (n = 19). Statistical significance between disease groups and HC was determined using a two-tailed student t test. *p < 0.05, **p < 0.01, ***p < 0.001. n.s. not significant


BD is a chronic autoimmune disease characterized by various clinical manifestations that may be similar to other diseases. Given the lack of specific serological markers, it is difficult to diagnose the disease early and to treat it. This study is the first to use a metabolomics approach for exploring the potential diagnostic markers of BD. Our study suggests that the altered levels of PCs and PUFAs may be indicative of the diagnosis of BD. Two n-6 fatty acids, LA and AA, may provide insight into therapeutic effects.

PCs, the major structural components of cell membranes, serve as fatty acid carriers and play an important role in metabolism and signaling [18, 19]. Table 3 lists PC and lysophosphatidylcholine (LPC) species that showed significantly different serum levels in BD patients compared with HC. PCs have been studied as potential metabolic biomarkers for the diagnosis of several diseases, such as calcific coronary artery disease [20] and endometriosis [21]. Since PCs are involved in pathogenic processes such as chronic inflammation, autoimmunity, and allergy [22,23,24], it has been suggested that PCs act as predictive metabolites corresponding with the activation of inflammatory, oxidant, and fibrotic pathways in progressive nephropathy [22]. Increased levels of various polyunsaturated PCs were positively associated with asthma [23]. In addition, PC/LPC ratios in plasma may be indicators of the early stages of RA, and they may be a reliable measure of inflammation. PC/LPC ratios could increase on therapy with tumor necrosis factor (TNF)α inhibitors [24]. In our study, the levels of PCs, PC(35:2), PC(36:6), and PC(P-40:6) in the pretreated BD group were lower than in the HC group. This may be associated with the hyperinflammatory status of BD. Intriguingly, however, the decreased level of PCs did not recover after glucocorticoid or immunosuppressant treatment. The implications of PCs in the pathophysiology of BD need to be further studied.

Table 3 The levels of ions with putative phosphatidylcholine (PC) identification in Behçet’s disease (BD) patients and healthy controls (HC)

The n-6 and n-3 PUFAs play an important role in the regulation of biological functions, inflammation, and immunity. Eicosanoids derived from n-6 PUFAs have a proinflammatory role, while those derived from n-3 PUFAs have an anti-inflammatory role [25]. It has been suggested that inflammatory and autoimmune diseases can be managed by regulating the intake of n-3 and n-6 PUFAs in the diet. In fact, modulation of the n-6/n-3 PUFA proportion is beneficial in several diseases, such as RA, ulcerative colitis, and cardiovascular diseases. It can decrease disease activity and minimize the requirements for anti-inflammatory drugs [26,27,28,29].

LA, one of the n-6 PUFAs, is an essential fatty acid because it cannot be synthesized in the human body. LA can be converted to the metabolically active AA, an n-6 PUFA that is present in the phospholipids of biomembranes. AA can be metabolized to several proinflammatory eicosanoids via multiple metabolic pathways, including the cyclooxygenase, lipoxygenase, and cytochrome P450 monooxygenases pathways [30]. AA can be involved in the regulation of inflammation through its eicosanoid metabolites, such as prostaglandin E2, thromboxane A2, and leukotriene B4 [31]. It is reported that AA-derived eicosanoids can reduce inflammatory Th17 and Th1 cell-mediated inflammation and improve colitis-associated immunopathology [32]. In our study, increased levels of LA and AA were found in pretreated BD patients compared with HC. This may reflect enhanced inflammation and relate to the occurrence and development of the disease. Our results showed a reduced level of two n-6 PUFAs in post-treatment BD patients, which indicated that these PUFAs, as indicators of inflammatory symptoms, may be useful for treatment assessment.

OA, an n-9 PUFA, is present in human plasma, cell membranes, and adipose tissue. OA can regulate physiological and pathological changes in cells through cell surface receptors or nuclear receptors [33, 34]. OA has been linked with metabolic and inflammatory diseases, and OA induces neutrophil accumulation and the release of inflammatory cytokines [35]. OA can also sensitize dendritic cells, resulting in augmented secretion of Th1/17 cytokines upon proinflammatory stimulation, and it can further promote an inflammatory response [36]. Our study suggests that OA may provide insights for the diagnosis and therapeutic effects of BD.

Ahn et al. [37] recently reported that the serum metabolite profiles of BD patients were distinctively separate from those of HC using gas chromatography with time-of-flight mass spectrometry (GC/TOF-MS). Five metabolic biomarkers, namely decanoic acid, fructose, tagatose, LA, and OA, were selected and validated as potential metabolite biomarkers for diagnosing BD. While GC/MS and LC/MS can be complimentary in terms of detecting different metabolites, the application of GC is limited to those who are volatile before or after derivatization. Our metabolomics profiling pointed to a different lipid metabolism in BD patients, so we designed UPLC-QTOF-MSE methods for lipidomics, which was not included in the study of Ahn et al. These differences in analytical methods may lead to different biomarkers from the previous study [37]. In addition to identifying the differential metabolites between BD patients and HC, we have also compared lipidomic profiles before and after treatment to search for potential biomarkers with therapeutic effects.

To further assess the diagnostic efficiency of these biomarkers, an independent validation cohort was employed. Since all patients were enrolled from a single center with relatively small sample sizes, we cannot exclude the possibility that our conclusions may have some specific limitations to the Chinese population. A multicenter study with a large sample size would therefore strengthen this study. In addition, we found that serum levels of LA and AA could distinguish BD patients from HC and other inflammatory or autoimmune diseases, including RA, SLE, TA, and CD, suggesting that these serum biomarkers might be specific markers for BD diagnosis.


In conclusion, our study supports the importance of PCs, LA, AA, and OA in the diagnosis and therapeutic effects of BD. This study is the first to use a metabolomic approach in the study of BD. Further investigations are required to explore the implication of these metabolomics alterations in the pathophysiology of BD.



Arachidonic acid


Area under the curve


Behçet’s disease

BDCAF 2006:

BD Current Activity Form 2006


Crohn’s disease


C-reactive protein


Erythrocyte sedimentation rate


Gas chromatography


Healthy controls


International Criteria for Behçet’s Disease


Internal standard


Kyoto Encyclopaedia of Genes and Genomes


Linoleic acid


Liquid chromatography




Tandem mass spectrometry


Oleic acid




Principal component analysis


Polyunsaturated fatty acid


Peking Union Medical College Hospital


Rheumatoid arthritis


Random Forests


Receiver operating characteristic


Systemic lupus erythematosus


Takayasu’s arteritis


Ultra-performance liquid chromatography-quadrupole time-of-flight mass spectrometry


  1. O'Duffy JD. Vasculitis in Behcet’s disease. Rheum Dis Clin N Am. 1990;16(2):423–31.

    CAS  Google Scholar 

  2. Nicholson JK, Lindon JC. Systems biology: metabonomics. Nature. 2008;455(7216):1054–6.

    Article  CAS  Google Scholar 

  3. Zhou L, Wang Q, Yin P, Xing W, Wu Z, Chen S, et al. Serum metabolomics reveals the deregulation of fatty acids metabolism in hepatocellular carcinoma and chronic liver diseases. Anal Bioanal Chem. 2012;403(1):203–13.

    Article  CAS  Google Scholar 

  4. Wu H, Xue R, Tang Z, Deng C, Liu T, Zeng H, et al. Metabolomic investigation of gastric cancer tissue using gas chromatography/mass spectrometry. Anal Bioanal Chem. 2010;396(4):1385–95.

    Article  CAS  Google Scholar 

  5. Xu J, Zhang J, Cai S, Dong J, Yang JY, Chen Z. Metabonomics studies of intact hepatic and renal cortical tissues from diabetic db/db mice using high-resolution magic-angle spinning 1H NMR spectroscopy. Anal Bioanal Chem. 2009;393(6–7):1657–68.

    Article  CAS  Google Scholar 

  6. Lian J-S, Liu W, Hao S-R, Chen D-Y, Wang Y-Y, Yang J-L, et al. A serum metabolomic analysis for diagnosis and biomarker discovery of primary biliary cirrhosis and autoimmune hepatitis. Hepatobiliary Pancreat Dis Int. 2015;14(4):413–21.

    Article  Google Scholar 

  7. Reinke S, Broadhurst D, Sykes B, Baker G, Catz I, Warren K, et al. Metabolomic profiling in multiple sclerosis: insights into biomarkers and pathogenesis. Mult Scler J. 2014;20(10):1396–400.

    Article  CAS  Google Scholar 

  8. Criteria for diagnosis of Behcet’s disease. International study group for Behcet’s disease. Lancet. 1990;335(8697):1078–80.

    Google Scholar 

  9. International Team for the Revision of the International Criteria for Behcet’s D. The international criteria for Behcet's disease (ICBD): a collaborative study of 27 countries on the sensitivity and specificity of the new criteria. J Eur Acad Dermatol Venereol. 2014;28(3):338–47.

    Article  Google Scholar 

  10. Aletaha D, Neogi T, Silman AJ, Funovits J, Felson DT, Bingham CO 3rd, 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  Google Scholar 

  11. Petri M, Orbai AM, Alarcón GS Gordon C, Merrill JT, Fortin PR, et al. Derivation and validation of the systemic lupus international collaborating clinics classification criteria for systemic lupus erythematosus. Arthritis Rheum. 2012;64(8):2677–86.

    Article  Google Scholar 

  12. Arend WP, Michel BA, Bloch DA, Hunder GG, Calabrese LH, Edworthy SM, et al. The American College of Rheumatology 1990 criteria for the classification of Takayasu arteritis. Arthritis Rheum. 1990;33:1129–34.

    Article  CAS  Google Scholar 

  13. Gomollón F, Dignass A, Annese V, Tilg H, Van Assche G, Lindsay JO, et al. 3rd European evidence-based consensus on the diagnosis and management of crohn's disease 2016: part 1: diagnosis and medical management. J Crohns Colitis. 2017;11(1):3–25.

    Article  Google Scholar 

  14. Goudarzi M, Weber WM, Chung J, Doyle-Eisele M, Melo DR, Mak TD, et al. Serum dyslipidemia is induced by internal exposure to strontium-90 in mice, lipidomic profiling using a data-independent liquid chromatography-mass spectrometry approach. J Proteome Res. 2015;14(9):4039–49.

    Article  CAS  Google Scholar 

  15. Li HH, Tyburski JB, Wang YW, Strawn S, Moon BH, Kallakury BV, et al. Modulation of fatty acid and bile acid metabolism by peroxisome proliferator-activated receptor alpha protects against alcoholic liver disease. Alcohol Clin Exp Res. 2014;38(6):1520–31.

    Article  CAS  Google Scholar 

  16. Mak TD, Laiakis EC, Goudarzi M, Fornace AJ. MetaboLyzer: a novel statistical workflow for analyzing postprocessed LC–MS metabolomics data. Anal Chem. 2014;86(1):506–13.

    Article  CAS  Google Scholar 

  17. Laiakis EC, Hyduke DR, Fornace AJ. Comparison of mouse urinary metabolic profiles after exposure to the inflammatory stressors gamma radiation and lipopolysaccharide. Radiat Res. 2012;177(2):187–99.

    Article  CAS  Google Scholar 

  18. Pietiläinen KH, Róg T, Seppänen-Laakso T, Virtue S, Gopalacharyulu P, Tang J, et al. Association of lipidome remodeling in the adipocyte membrane with acquired obesity in humans. PLoS Biol. 2011;9(6):e1000623.

    Article  Google Scholar 

  19. Wymann MP, Schneiter R. Lipid signalling in disease. Nat Rev Mol Cell Biol. 2008;9(2):162–76.

    Article  CAS  Google Scholar 

  20. Vorkas PA, Isaac G, Holmgren A, Want EJ, Shockcor JP, Holmes E, et al. Perturbations in fatty acid metabolism and apoptosis are manifested in calcific coronary artery disease: an exploratory lipidomic study. Int J Cardiol. 2015;197:192–9.

    Article  Google Scholar 

  21. Vouk K, Ribič-Pucelj M, Adamski J, Rižner TL. Altered levels of acylcarnitines, phosphatidylcholines, and sphingomyelins in peritoneal fluid from ovarian endometriosis patients. J Steroid Biochem Mol Biol. 2016;159:60–9.

    Article  CAS  Google Scholar 

  22. Zhao YY, Wang HL, Cheng XL, Wei F, Bai X, Lin RC, et al. Metabolomics analysis reveals the association between lipid abnormalities and oxidative stress, inflammation, fibrosis, and Nrf2 dysfunction in aristolochic acid-induced nephropathy. Sci Rep. 2015;5:12936.

    Article  CAS  Google Scholar 

  23. Ried JS, Baurecht H, Stuckler F, Krumsiek J, Gieger C, Heinrich J, et al. Integrative genetic and metabolite profiling analysis suggests altered phosphatidylcholine metabolism in asthma. Allergy. 2013;68(5):629–36.

    Article  CAS  Google Scholar 

  24. Fuchs B, Schiller J, Wagner U, Häntzschel H, Arnold K. The phosphatidylcholine/lysophosphatidylcholine ratio in human plasma is an indicator of the severity of rheumatoid arthritis: investigations by 31P NMR and MALDI-TOF MS. Clin Biochem. 2005;38(10):925–33.

    Article  CAS  Google Scholar 

  25. Calder PC. Polyunsaturated fatty acids, inflammation, and immunity. Lipids. 2001;36(9):1007–24.

    Article  CAS  Google Scholar 

  26. Wall R, Ross R, Fitzgerald G, Stanton C. Fatty acids from fish: the anti-inflammatory potential of long-chain omega-3 fatty acids. Nutr Rev. 2010;68(5):280–9.

    Article  Google Scholar 

  27. Simopoulos AP. Omega-3 fatty acids in inflammation and autoimmune diseases. J Am Coll Nutr. 2002;21(6):495–505.

    Article  CAS  Google Scholar 

  28. Simopoulos AP. The importance of the omega-6/omega-3 fatty acid ratio in cardiovascular disease and other chronic diseases. Exp Biol Med. 2008;233(6):674–88.

    Article  CAS  Google Scholar 

  29. Investigators IBDiES, Tjonneland A, Overvad K, Bergmann MM, Nagel G, Linseisen J, et al. Linoleic acid, a dietary n-6 polyunsaturated fatty acid, and the aetiology of ulcerative colitis: a nested case–control study within a European prospective cohort study. Gut 2009;58(12):1606.

  30. Subhash PK, David SG, David RJ, Letts LG. Eicosanoids in inflammation: biosynthesis, pharmacology, and therapeutic frontiers. Curr Top Med Chem. 2007;7(3):311–40.

    Article  Google Scholar 

  31. Bogatcheva NV, Sergeeva MG, Dudek SM, Verin AD. Arachidonic acid cascade in endothelial pathobiology. Microvasc Res. 2005;69(3):107–27.

    Article  CAS  Google Scholar 

  32. Monk JM, Turk HF, Fan Y-Y, Callaway E, Weeks B, Yang P, et al. Antagonizing arachidonic acid-derived eicosanoids reduces inflammatory Th17 and Th1 cell-mediated inflammation and colitis severity. Mediat Inflamm. 2014;2014:14.

    Article  Google Scholar 

  33. Lopez S, Bermudez B, Montserrat-de la Paz S, Jaramillo S, Varela LM, Ortega-Gomez A, et al. Membrane composition and dynamics: a target of bioactive virgin olive oil constituents. Biochim Biophys Acta. 2014;1838(6):1638–56.

    Article  CAS  Google Scholar 

  34. Catalá A. Five decades with polyunsaturated fatty acids: chemical synthesis, enzymatic formation, lipid peroxidation and its biological effects. J Lipids. 2013;2013:19.

    Article  Google Scholar 

  35. Gonçalves-de-Albuquerque CF1, Silva AR, Burth P, de Moraes IM, Oliveira FM, Younes-Ibrahim M, et al. Oleic acid induces lung injury in mice through activation of the ERK pathway. Mediat Inflamm. 2012;2012:11.

    Google Scholar 

  36. Stelzner K, Herbert D, Popkova Y, Lorz A, Schiller J, Gericke M, et al. Free fatty acids sensitize dendritic cells to amplify TH1/TH17-immune responses. Eur J Immunol. 2016;46(8):2043–53.

    Article  CAS  Google Scholar 

  37. Ahn JK, Kim J, Hwang J, Song J, Kim KH, Cha H-S. Potential metabolomic biomarkers for reliable diagnosis of Behcet's disease using gas chromatography/time-of-flight-mass spectrometry. Joint Bone Spine. 2017.

    Article  CAS  Google Scholar 

Download references


We thank the health professional staff from the Department of Rheumatology & Clinical Immunology, Peking Union Medical College Hospital, and the patients for their participation in this study. We are especially thankful to the healthy volunteers for the donation of blood. We also thank the Proteomic and Metabolomics Shared Resources at Georgetown University, NIH P30CA51008, for providing UPLC-QTOF-MS service.


This study was supported by the National Institute of Allergy and Infectious Diseases (NIAID grant number U19 A1067773–09), the National Natural Science Foundation of China (grant number 81571598), the National Key Research and Development Program: “Precise Medical Research” (grant number 2016YFC0906201), and CAMS Innovation Fund for Medical Sciences (CIFMS; grant number 2016-I2M-1-013). The funders had no role in study design, data collection, data analysis, data interpretation, or writing of the manuscript.

Availability of data and materials

All data generated and/or analyzed during this study are included in this published article and its supplementary information files.

Author information

Authors and Affiliations



WZ and H-HL conceptualized, designed, and supervised the project. XW, WZ, and H-HL performed the experiments and wrote the manuscript. H-HL analyzed the metabolomics and lipidomics data. MG performed validations for PCs, AA, and LA. JS, CL, and JL participated in the sample collection and clinical analysis. WS performed further validation tests. HC, XZ, and XfZ critically reviewed the manuscript and provided valuable input. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Wenjie Zheng or Heng-Hong Li.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the institutional committee for the Protection of Human Subjects from PUMCH. All subjects gave written informed consent in accordance with the Declaration of Helsinki.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional files

Additional file 1:

Verification of PCs by multiple reaction monitoring. These panels show MS/MS spectra of the indicated ions. Multiple reaction monitoring transitions were monitored for PC signature fragmentation (m/z 184). (PDF 46 kb)

Additional file 2:

Verification of PUFAs by MS/MS. Retention time of two n-6 PUFAs, linoleic acid and arachidonic acid were compared with that of the pure chemicals. MS/MS spectra are shown. (PDF 46 kb)

Additional file 3:

Treatment reverses the increased level of oleic acid in serum. (A) Abundance of oleic acid in healthy volunteers, pretreatment BD (Pre-) patients, and post-treatment BD (Post-) patients. ***p < 0.001. (B) Verification of oleic acid by MS/MS. (PDF 59 kb)

Additional file 4:

The ROC curve of PCs, AA, and LA in BD patients. (A) The ROC curve of PCs with area under the curve (AUC) > 0.85 in BD patients. (B) The ROC curve of AA and LA in BD patients. (PDF 18 kb)

Additional file 5:

Comparison of the content of AA and LA in an independent cohort. (A) Identification of AA and LA in serum samples by comparison with reference standards. (B) Calibration curves of AA and LA with rosmarinic acid as the internal standard. (C) Contents of LA and AA in BD and HC samples. (D) Correlation of ESR and CRP with the serum levels of LA and AA in BD. (PDF 253 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zheng, W., Wu, X., Goudarzi, M. et al. Metabolomic alterations associated with Behçet’s disease. Arthritis Res Ther 20, 214 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: