
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39251615
68477
10.1038/s41598-024-68477-6
Article
Red blood cell metabolomics identify ergothioneine as a key metabolite in DMARD-naïve rheumatoid arthritis and response to methotrexate
Sigaux Johanna 12
Junot Christophe 3
Boissier Marie-Christophe marie-christophe.boissier@aphp.fr

12
Petit Mylène 2
Breckler Magali 1
Castelli Florence 3
Fenaille François 3
Roméo Paul-Henri 4
Semerano Luca 12
1 https://ror.org/0199hds37 grid.11318.3a 0000 0001 2149 6883 Inserm U1125, Université Sorbonne Paris Nord, Bobigny, France
2 https://ror.org/00pg5jh14 grid.50550.35 0000 0001 2175 4109 Rheumatology Department, Assistance Publique-Hôpitaux de Paris, GH HUPSSD, Bobigny, France
3 grid.457334.2 0000 0001 0667 2738 Service de Pharmacologie et Immuno-Analyse (SPI), Laboratoire d’Etude du Métabolisme des Médicaments, CEA, INRA, Université Paris Saclay, MetaboHUB, 91191 Gif-Sur-Yvette, France
4 https://ror.org/010j2gw05 grid.457349.8 0000 0004 0623 0579 Institut de Biologie François Jacob (IBFJ), CEA, 92260 Fontenay-Aux-Roses, France
9 9 2024
9 9 2024
2024
14 2100522 3 2024
24 7 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Using a new red blood cell (RBC) metabolite extraction protocol, we performed a metabolomic analysis on RBCs in rheumatoid arthritis (RA) patients treated or not with methotrexate (MTX), with the two following objectives: to compare the RBC metabolic profiles of MTX-naïve RA patients and healthy controls (HC), and to investigate whether RBC profiles before and after MTX treatment in RA differed between responders and non-responders. Plasma analysis was performed in parallel. Metabolites were extracted and identified in RBCs and plasma by liquid chromatography-mass spectrometry. We compared the metabolomic fingerprints of 31 DMARD-naïve RA patients and 39 HCs. We also compared the RBC and plasma metabolomes of 25 RA patients who responded or not to MTX therapy before (M0) and after a 3-month treatment period (M3). Significance was determined by Storey’s false discovery rate (FDR) q-values to correct for multiple testing. RA patients and HCs differed in the metabolomic signature of RBCs. The signature mainly contained amino acids (AA). Eleven metabolites, including 4 metabolites belonging to the carbohydrate subclass and 2 amino acids (creatine and valine) showed accumulation in RBCs from RA patients. Conversely, citrulline (fold change = 0.83; q = 0.025), histidine (fold change = 0.86; q = 0.014) and ergothioneine (EGT) (fold change = 0.66; q = 0.024), were lower in RBC of RA patients. Five plasma metabolites, including succinic acid and hydroxyproline, were higher in RA patients, and 7 metabolites, including DHEA sulfate, alanine, threonine and ornithine, were lower. Among RA patients undergoing MTX treatment pre-treatment (M0), EGT values were significantly lower in non-responders. In conclusion, low RBC levels of EGT, a food-derived AA barely detectable in plasma, characterize DMARD naïve RA patients and lack of response to MTX treatment.

Keywords

Rheumatoid arthritis
Metabolomic
Ergothioneine
Diet
Environment
Subject terms

Metabolomics
Rheumatology
Société Française de RhumatologieBourse "Mobilité" 2021 Sigaux Johanna issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Omics sciences are increasingly used to study complex multifactorial diseases. Among these omics, metabolomics study the whole set of metabolites in a biological system that define its metabolome. Given that metabolites are the final products of cellular regulatory pathways, the metabolome integrates and reflects all the interactions between biochemical entities in cells or tissues and their environment. Thus, metabolomics enable the qualitative and quantitative description of a biological system’s ultimate response to environmental changes1. With regard to rheumatoid arthritis (RA), metabolomics can identify specific signatures of the disease (e.g. RA vs. healthy controls (HC)), trace putative metabolicpathways relevant for disease pathogenesis and treatment, and provide additional information on how patients with RA respond to the environment, including diet and exposure to xenobiotic substances. Several metabolomics studies have been published in the field of RA and have identified a large number of metabolites as putative markers of RA2, response to treatment3, and disease activity4. Most of them are involved in glycolysis (lactic acid, glucose), citrate cycle (citric acid), amino acid (leucine, isoleucine, valine, threonine, serine, proline, alanine, histidine) and lipid (cholesterol, glycerol) metabolisms5. All the RA metabolomic studies were performed on serum, plasma, urine and synovial fluid, but none of these studies examined red blood cells (RBC) as an alternative matrix.

RBCs play a major role in interorgan transport of metabolites during their life span, thereby reflecting exogenous metabolic processes. A case in point is the measure of glycated hemoglobin A1c, which reflects mean glycemia over the past 3 months. RBCs are basic cells devoid of nucleus, mitochondria and ribosome. They therefore have an elementary metabolism.

In this proof-of-concept study, we developed a metabolite extraction protocol to perform metabolomic analysis on RBCs with three objectives. The first was to characterize any differences in the RBC metabolic profiles of Disease-Modifying antirheumatic drugs (DMARD)-naïve RA patients and HCs. The second was to determine whether RBC profiles before and after methotrexate (MTX) treatment in RA patients differed between responders and non-responders. The third was to identify relevant RBC discriminant metabolites. In parallel, we performed metabolomic analysis of the plasma to demonstrate the diagnostic value of RBC metabolome and to study any association between RBC and plasma metabolomes.

Results

Clinical parameters

RA patients were predominantly female (71.0%), with a mean age of 50.0 years, and mostly ACPA (92.3%) and RF (95.2%) positive. Disease activity before MTX initiation was moderate to high (mean DAS28-CRP: 4.7 ± 1.4) (Table 1). Of the 31 initial RA patients, 25 were treated with MTX for at least 3 months. After a 3-month MTX treatment period, 12 patients were classified as responders and 13 as non-responders. The two groups before treatment did not significantly differ in baseline DAS28-CRP (4.6 ± 1.1 vs. 3.8 ± 1.5, p = 0.17 for responders and non-responders, respectively) or in tender and swollen joint counts (Supplementary Table S1). Table 1 Baseline RA patient characteristics (M0).

	RA (N = 31)	
Women, n (%)	22 (71.0)	
Age (years), mean (sd)	50.0 (18.8)	
Disease duration (years), mean (sd)	4.9 (7.0)	
ACPA + , n (%)	24 (92.3)	
RF + , n (%)	20 (95.2)	
History or actual MTX or other DMARDs, n (%)	0 (0)	
History or actual bDMARDs, n (%)	0 (0)	
Oral steroid, n (%)	14 (45.2)	
Oral steroid dose* (mg), mean (sd)	4.9 (6.2)	
TJC (/28), mean (sd)	8.3 (7.0)	
SJC (/28), mean (sd)	6.3 (4.3)	
VAS pain (/100), mean (sd)	58.4 (29.3)	
ESR (mm), mean (sd)	36.7 (33.8)	
CRP (mg/L), mean(sd)	22.1 (26.6)	
DAS28-CRP, mean (sd)	4.7 (1.4)	
Type II diabetes, n (%)	2 (6.7)	
Dyslipidaemia, n (%)	5 (16.7)	
High blood pressure, n (%)	2 (6.7)	
*prednisone equivalent.

ACPA anti-citrullinated protein antibodies, RF rheumatoid factor, TJC tender joint count, SJC swollen joint count, VAS visual analogue scale, ESR erythrocyte sedimentation rate, CRP C-reactive protein.

RBC metabolites associated with RA (vs. healthy controls)

Under our conditions, 33 metabolite features from the C18( +) and 74 from the HILIC(-) analysis matched the accurate masses and retention times of the metabolites included in our chemical database, yielding 105 unique annotated metabolites (See supplementary file). Univariate analysis of the differences in mean expression between the two groups showed a significant difference in HILIC(−) analysis (p-value < 0.05) for 14 metabolites (Fig. 1A, Table 2) and all of them remained significant after correction for multiple testing (FDR q-value < 0.05) (Supplementary Figure S1). Three of them – histidine, ergothioneine (EGT), and citrulline—were found to be less expressed in RBCs from RA patients (RA/HC ratio: 0.86, 0.66 and 0.83, respectively, with all three FDR qvalue < 0.02). Meanwhile, eleven—including 4 metabolites in the carbohydrate subclass and 2 amino acids (creatine and valine)—showed accumulation. Although the results obtained under C18( +) conditions were statistically less robust (Fig. 1B, Supplementary Table S2), the same trend was observed for EGT (RA/HC ratio: 0.70, p = 0.040) and citrulline (RA/HC ratio: 0.83, p = 0.027), showing lower concentrations in the RA group compared to HCs. These differences did not reach FDR-adjusted significance. In this dataset, 5 metabolites involved in the arginine and proline pathway showed accumulation in RA RBCs compared to HCs (aminovaleric acid, creatine, spermidine, spermine and stachydrine). Metabolites were also studied in plasma samples from RA patients and HCs (See supplementary Data, tables S3 and S4). We found lower concentrations of several amino acids in the plasma of RA patients compared to HCs (DHEA, alanine, threonine, ornithine, theophylline) and higher hydroxyproline.Figure 1 Red blood cell metabolite levels in RA patients compared to HCs Volcano plots showing the results of comparisons of RBC metabolite levels (using the HILIC(−) dataset) in RA patients relative to HCs. The vertical dashed lines indicate the threshold for the 0.2Log2 fold abundance difference. The horizontal dashed line indicates the p = 0.05 threshold. Between-group comparisons as described in statistical methods. Figure 1A shows the results obtained with the HILIC(−) dataset. Figure 1B shows the results obtained with the C18( +) dataset.

Table 2 RBC metabolites associated with RA (HILIC(-) dataset), ranked lowest to highest RA/HC ratio.

Metabolite name	KEGG ID¤	Class¤	Subclass¤	Pathway¤	Pvalue*	FDR-Qvalue*	RA/HC ratio	
Ergothioneine	C05570	Carboxylic acids and derivatives	Amino acids, peptides, and analogues	Histidine metabolism	0.016	0.024	0.66	
Citrulline	C00327	Carboxylic acids and derivatives	Amino acids, peptides, and analogues	Arginine and proline metabolism	0.018	0.025	0.83	
Histidine	C00768	Carboxylic acids and derivatives	Amino acids, peptides, and analogues	Histidine metabolism	0.005	0.014	0.86	
NADP	C00006	Phenols	Benzenediols	Histidine metabolism	0.046	0.051	1.14	
Hexoses-bisphosphate	C01231	Organooxygen compounds	Carbohydrates and carbohydrate conjugates	Amino sugar and nucleotide sugar metabolism	0.011	0.020	1.16	
Creatine	C00300	Carboxylic acids and derivatives	Amino acids, peptides, and analogues	Arginine and proline metabolism	0.019	0.025	1.17	
Uridine 5’-diphosphogalactose	C00052	Pyrimidine nucleotides	Pyrimidine nucleotide sugars	Pyrimidine metabolism	0.005	0.014	1.19	
Gluconic acid	C00257	Organooxygen compounds	Carbohydrates and carbohydrate conjugates	Ketogluconate metabolim	 < 0.001	0.002	1.22	
Hydroxyglutaric acid	C02630	Hydroxy acids and derivatives	Short-chain hydroxy acids and derivatives	Butanoate metabolism	0.011	0.020	1.22	
Phosphoenolpyruvic acid	C00074	Organic phosphoric acids and derivatives	Phosphate esters	Pyruvate metabolism	0.007	0.016	1.23	
Diphosphoglyceric acid	C01159	Organooxygen compounds	Carbohydrates and carbohydrate conjugates	Gluconeogenesis	 < 0.001	0.002	1.29	
Glyceraldehyde -phosphate	C00118	Organooxygen compounds	Carbohydrates and carbohydrate conjugates	Gluconeogenesis	 < 0.001	0.002	1.29	
Valine or betaine	C00719	Carboxylic acids and derivatives	Amino acids, peptides, and analogues	Glycolysis	0.040	0.048	1.30	
Guanosine – diphospho-mannose/Guanosine diphosphoglucose (GDP-Man/GDP-Gluc)	C00096/C00394	Purine nucleotides	Purine ribonucleotides	Purine Metabolism	 < 0.001	 < 0.001	1.86	
*Significant changes in metabolite levels were tested using the t-test, and false discovery rate (FDR)-q value were determined according to Storey and Tibshirani.

¤Metabolite ID, class, subclass and pathway were provided by the human metabolome database (HMDB).

High level of RBC ergothioneine is associated with favorable methotrexate (MTX) response

Since RBC analysis identified EGT as a differentiating metabolite between RA patients and HCs, we investigated whether the concentration of EGT differed between responders and non-responders before (M0) and/or after (M3) MTX initiation. Of the 31 included patients, 25 completed the follow-up and were available for analysis. Among the 6 patients that did not complete the study, 4 missed their 3 month follow-up visit and 2 did not tolerate MTX. Those 6 patients, had similar values of DAS28-CRP (p = 0.47), mean corticosteroid dose (p = 0.21) and EGT levels (p = 0.10) as those that completed their follow-up.

A total of 12 patients responded to MTX, while 13 were classified as non-responders. We found significantly lower EGT values in non-responders at M0 (p = 0.036) when using the HILIC(−) dataset (Fig. 2). A similar trend was observed using C18( +) (p = 0.07, Wilcoxon test). Responders and non-responders did not differ in baseline corticosteroid dose (p = 0.84) or in baseline disease activity expressed by DAS28CRP values (p = 0.17) (supplementary table 1). Higher that the median EGT was the only variable associated with clinical response to methotrexate (OR 6.00; 95% CI 1.02 to 35.38) in multivariable logistics regression model taking into account these two confounders.Figure 2 Boxplots illustrating differential ergothioneine level (using the HILIC(−) dataset) between methotrexate non-responder and responder RA patients before MTX onset (M0). *Between-group comparisons were performed using Mann–Whitney test (two-sided), p = 0.036.

Discussion

This study is the first RBC metabolomics characterization in RA. This proof-of-concept study identified EGT—an amino acid barely detectable in plasma—as a key metabolite that is reduced in untreated RA and particularly in future non-responders to MTX. In alignment with several RA metabolomic studies, our RBC and plasma analysis shows a predominantly amino acidic RA signature. Since the RBC metabolome is not limited to pathways essential for RBC functions, but also includes metabolites from external sources such as plasma and other cells, the observed variations in AA concentration may reflect impaired systemic metabolism.

We found a significant decrease in EGT, histidine (EGT precursor), and citrulline in RA patient RBCs vs. HCs. This decrease may be related to the predominance of the ammonia recycling pathway and/or the nitric oxide (NO) pathway involved in RA, and may reflect increased catabolism and resulting muscle loss6. Interestingly, the plasma metabolome also revealed the potential role of the ammonia-recycling pathway, notably through the lower ornithine in RA patients. Furthermore, the fatigue frequently experienced by RA patients has been associated with downregulation of the urea cycle that involves citrulline7. Citrulline concentration was lower in RA patients’ RBC metabolome compared to HCs. Since RBCs do not have any NO synthase activity, the citrulline may be coming from plasma. The presence of endothelial dysfunction in RA patients, even at an early stage, is well established and it has been suggested that patients may exhibit impaired NO reactivity8. Thus, we can hypothesize that the citrulline decrease reflects the endothelial dysfunction observed in untreated RA and may be a consequence of decreased NO synthase activity and/or lack of arginine (NO synthase substrate) intake.

The metabolomic analysis performed on plasma in our study is consistent with the relevant literature. In a meta-analysis of 10 publications reporting metabolomic analyses, Li et al. identified 17 metabolites that could differentiate RA patients and controls5. This study confirmed two of them: threonine and alanine. Of note, these metabolites were significantly reduced in collagen-induced arthritis mice compared to controls9. We also identified DHEA sulfate as the metabolite with the lowest RA/HC ratio, in alignment with the literature10.

The data obtained from the RBC metabolome analysis is therefore partly consistent with the plasma analysis (involvement of the ammonia recycling pathway), but also identified other metabolites and pathways that appear to be dysregulated in RA patients, especially EGT. Our analyses showed a lower EGT concentration in RA RBCs compared to HC RBCs, and a higher EGT concentration in MTX responders compared to non-responders before MTX onset. EGT is a natural AA that derives from histidine. It is strongly associated with a “health conscious food pattern” that was linked to a reduced risk of cardiometabolic disease in a 3,236-participant study11. EGT cannot be synthesized by humans and is provided through diet. It can be found in actinobacteria, cyanobacteria, some fungi (such as porcini), oysters and garlic. EGT is a powerful antioxidant and can induce the expression of SIRT112, which appears lower in RA endothelial and synovial vessels13.

EGT uses OCTN1—a dependent organic cation transporter primarily expressed in blood cells including erythrocytes—to enter cells14. EGT is therefore up to 100 times more concentrated in RBCs than in plasma, where it is barely detectable. Interestingly, a study highlighted a polymorphism of SLC22A4 (coding for OCTN1) strongly associated with RA in East Asian15 but not Caucasian populations16. The functional consequences of this polymorphism are unknown, as is the actual biological role of EGT.

Our result differs from a study showing that the erythrocyte concentration of ergothioneine in RA patients was higher than in HCs17. This discrepancy may be related to the cohorts used, since the patients included in this previous study were nearly all treated with MTX and showed low disease activity (i.e. responders). Conversely, our study focused on DMARD-naïve RA patients. If, as our results suggest, patients with higher levels of EGT have a higher probability of responding to methotrexate treatment, it is conceivable that these patients are also more likely to remain on treatment over the long term. Therefore, by selecting patients in remission on methotrexate treatment, the other study might have focused on RA patients with the highest levels of EGT.

Our work has some limitations, mainly the absence of matching of RA patients and HCs on potential confounders like smoking status and BMI.

Moreover, like most proof-of-concept studies, this work is underpowered to predict the response to methotrexate. Based on our ratio of responders to non-responders, an ideal number of 59 subjects should have been included to achieve 80% power with an alpha error of 0.05. Also, our results were not replicated in a validation cohort. However, identifying a metabolic signature that may predict response to methotrexate treatment in early RA is beyond the scope of this work and would require larger cohorts and adjustments for multiple confounders.

Another limitation is the fact that the biological role of EGT and its transporter in the pathogenesis of RA is still unknown. The observation that RA patients have lower levels of EGT compared to healthy controls, and that those with the lowest levels are also less likely to respond to methotrexate treatment, is intriguing. This finding warrants mechanistic studies to understand the biological role of this metabolite. Whether EGT may be a marker of refractory RA, a potential adjuvant to methotrexate treatment, or merely a random bystander cannot be established by this work.

However, this completely novel method in rheumatology preliminarily identifies EGT, a food-derived AA that has antioxidant properties and accumulates in tissues, as a potential metabolite of interest in untreated RA and MTX response. Larger studies are warranted to bolster these original data.

Methods

Patients

The study group consisted of 31 prospectively recruited consecutive patients that met the 2010 ACR-EULAR Classification Criteria for RA and were DMARD-naïve. All patients were recruited from the Avicenne Rheumatology Department (Bobigny, France). Following an overnight fast, RA patients underwent blood sampling and clinical assessment before (M0) and 3 months after (M3) MTX initiation. At M3, all patients with a good or moderate DAS28-CRP-based EULAR response or in DAS28-CRP-based remission (i.e. DAS28-CRP < 2.6)) were classified as MTX responders. Only patients who had completed 3 months of MTX treatment and completed their 3-month follow-up visit were included in the longitudinal analysis (i.e. responders vs. non-responders). Blood samples of age- and sex-matched blood donors (Etablissement Français du Sang) were used as healthy controls (HC). This study was approved by the local review board of Avicenne Hospital (NI-2016-11-01) and performed in accordance with the Declaration of Helsinki and the institutional guidelines. Informed consent was obtained from all subjects.

Preparation of RBC lysates and RBC and plasma metabolite extraction

Preparation of RBC lysates and metabolite extraction were performed as previously described (see supplementary methods).

Untargeted metabolomics by ultra-high performance liquid chromatography (UHPLC) coupled with high-resolution mass spectrometry (HRMS) and data processing

The ultra-high performance liquid chromatographic (UHPLC) separation was performed on a Hypersil GOLD C18 1.9 µm, 2.1 mm × 150 mm column (RP) at 30 °C (Thermo Fisher Scientific, les Ulis, France), and HPLC chromatographic separations were performed on a Sequant ZICpHILIC 5 µm, 2.1 × 150 mm (HILIC) at 15 °C (Merck, Darmstadt, Germany). All chromatographic systems were equipped with an on-line prefilter (Thermo Fisher Scientific, Courtaboeuf, France). The data processing is described elsewhere18 and reported in the supplementary methods. LC–MS analyses were performed using a U3000 liquid chromatography system coupled to an Exactive mass spectrometer from Thermo Fisher Scientific (Courtaboeuf, France) fitted with an electrospray source. The software interface was Xcalibur (version 2.1) (Thermo Fisher Scientific, Courtaboeuf, France) (see supplementary methods).

Data processing and metabolite annotation procedures are described elsewhere (see supplementary methods).

Statistical analysis

Metabolite concentrations expressed as chromatographic peak areas (HILIC (-) or C18 ( +)) were analyzed using R software (version 4.2.1)19. For a given metabolite, the normality of the data distributions for each of the two groups—RA and HC—was evaluated by the Shapiro–Wilk test and visual examination of the qqplot graphs. The equality of variance was tested using the var.test function (R stats library). If the hypothesis of normality of at least one of the two distributions was rejected, the comparison of the values of the two groups was done by the row_Wilcoxon_two-sample test (R MatrixTest library). Otherwise, this comparison was done by Welch t test (unequal variances) or Student t test. Significant changes in metabolite levels were tested using the t-test, and False Discovery Rate (FDR) were determined according to Storey and Tibshirani procedure to control false discovery rate in multiple statistical tests. For the comparison between responders and non-responders, baseline corticosteroid doses and DAS28CRP values were compared with independent sample T test after checking for normality. A multivariable logistic regression model was used to check the independent association between EGT baseline values and 3-month clinical response to methotrexate adjusting on baseline corticosteroids and DAS28CRP values. To check for a potential attrition bias, DAS28-CRP, mean corticosteroid dose and EGT level were compared between patients lost-to-follow-up and those who completed the study by Mann–Whitney test to the sample size.

Supplementary Information

Supplementary Information 1.

Supplementary Information 2.

Supplementary Information 3.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-68477-6.

Author contributions

PHR, MCB and LS conceptualized the study JS, PHR and LS created the methodology for the study JS, MP, MB, FC, FF carried out the data collection JS and LS did the data analysis and verified the data JS and CJ supervised the data collection JS and LS wrote the original draft All authors reviewed and edited the manuscript All authors had full access to all of the data and the final responsibility to submit for publication.

Funding

This work was supported by the French Rheumatology Society (Societé Française de Rhumatologie, SFR) and Entrepreneurs&Go.

Data availability

The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Competing interests

The authors declare no competing interests.

Publisher's note

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

1. Semerano L Roméo P-H Boissier M-C Metabolomics for rheumatic diseases: Has the time come? Ann. Rheum. Dis. 2015 74 1325 1326 10.1136/annrheumdis-2014-206618 25609411
Semerano, L., Roméo, P.-H. & Boissier, M.-C. Metabolomics for rheumatic diseases: Has the time come?. Ann. Rheum. Dis. 74, 1325–1326 (2015).25609411 10.1136/annrheumdis-2014-206618
2. Chu SH Circulating plasma metabolites and risk of rheumatoid arthritis in the Nurses’ Health Study Rheumatology (Oxford) 2020 59 3369 3379 10.1093/rheumatology/keaa125 32310291
Chu, S. H. et al. Circulating plasma metabolites and risk of rheumatoid arthritis in the Nurses’ Health Study. Rheumatology (Oxford) 59, 3369–3379 (2020).32310291 10.1093/rheumatology/keaa125
3. Priori R 1H-NMR-based metabolomic study for identifying serum profiles associated with the response to etanercept in patients with rheumatoid arthritis PLoS One 2015 10 e0138537 10.1371/journal.pone.0138537 26558759
Priori, R. et al. 1H-NMR-based metabolomic study for identifying serum profiles associated with the response to etanercept in patients with rheumatoid arthritis. PLoS One 10, e0138537 (2015).26558759 10.1371/journal.pone.0138537
4. Hur B Plasma metabolomic profiling in patients with rheumatoid arthritis identifies biochemical features predictive of quantitative disease activity Arthritis Res. Ther. 2021 23 164 10.1186/s13075-021-02537-4 34103083
Hur, B. et al. Plasma metabolomic profiling in patients with rheumatoid arthritis identifies biochemical features predictive of quantitative disease activity. Arthritis Res. Ther. 23, 164 (2021).34103083 10.1186/s13075-021-02537-4
5. Li C Metabolomics in the development and progression of rheumatoid arthritis: A systematic review Joint Bone Spine 2020 87 425 430 10.1016/j.jbspin.2020.05.005 32473419
Li, C. et al. Metabolomics in the development and progression of rheumatoid arthritis: A systematic review. Joint Bone Spine 87, 425–430 (2020).32473419 10.1016/j.jbspin.2020.05.005
6. Farrow M Muscle deterioration due to rheumatoid arthritis: Assessment by quantitative MRI and strength testing Rheumatology (Oxford) 2021 60 1216 1225 10.1093/rheumatology/keaa364 32910153
Farrow, M. et al. Muscle deterioration due to rheumatoid arthritis: Assessment by quantitative MRI and strength testing. Rheumatology (Oxford) 60, 1216–1225 (2021).32910153 10.1093/rheumatology/keaa364
7. Surowiec I Metabolomics study of fatigue in patients with rheumatoid arthritis naïve to biological treatment Rheumatol. Int. 2016 36 703 711 10.1007/s00296-016-3426-2 26803313
Surowiec, I. et al. Metabolomics study of fatigue in patients with rheumatoid arthritis naïve to biological treatment. Rheumatol. Int. 36, 703–711 (2016).26803313 10.1007/s00296-016-3426-2
8. Bergholm R Impaired responsiveness to NO in newly diagnosed patients with rheumatoid arthritis Arterioscler. Thromb. Vasc. Biol. 2002 22 1637 1641 10.1161/01.ATV.0000033516.73864.4E 12377742
Bergholm, R. et al. Impaired responsiveness to NO in newly diagnosed patients with rheumatoid arthritis. Arterioscler. Thromb. Vasc. Biol. 22, 1637–1641 (2002).12377742 10.1161/01.ATV.0000033516.73864.4E
9. He M Spontaneous ultra-weak photon emission in correlation to inflammatory metabolism and oxidative stress in a mouse model of collagen-induced arthritis J. Photochem. Photobiol. B 2017 168 98 106 10.1016/j.jphotobiol.2016.12.036 28199905
He, M. et al. Spontaneous ultra-weak photon emission in correlation to inflammatory metabolism and oxidative stress in a mouse model of collagen-induced arthritis. J. Photochem. Photobiol. B 168, 98–106 (2017).28199905 10.1016/j.jphotobiol.2016.12.036
10. Lee Y-J A discovery of screening markers for rheumatoid arthritis by liquid chromatography mass spectrometry: A metabolomic approach Int. J. Rheum. Dis. 2020 23 1353 1362 10.1111/1756-185X.13935 32845094
Lee, Y.-J. et al. A discovery of screening markers for rheumatoid arthritis by liquid chromatography mass spectrometry: A metabolomic approach. Int. J. Rheum. Dis. 23, 1353–1362 (2020).32845094 10.1111/1756-185X.13935
11. Smith E Ergothioneine is associated with reduced mortality and decreased risk of cardiovascular disease Heart 2020 106 691 697 10.1136/heartjnl-2019-315485 31672783
Smith, E. et al. Ergothioneine is associated with reduced mortality and decreased risk of cardiovascular disease. Heart 106, 691–697 (2020).31672783 10.1136/heartjnl-2019-315485
12. D’Onofrio N Ergothioneine oxidation in the protection against high-glucose induced endothelial senescence: Involvement of SIRT1 and SIRT6 Free Radic. Biol. Med. 2016 96 211 222 10.1016/j.freeradbiomed.2016.04.013 27101740
D’Onofrio, N. et al. Ergothioneine oxidation in the protection against high-glucose induced endothelial senescence: Involvement of SIRT1 and SIRT6. Free Radic. Biol. Med. 96, 211–222 (2016).27101740 10.1016/j.freeradbiomed.2016.04.013
13. Leblond A Implication of the deacetylase sirtuin-1 on synovial angiogenesis and persistence of experimental arthritis Ann. Rheum. Dis. 2020 79 891 900 10.1136/annrheumdis-2020-217377 32381568
Leblond, A. et al. Implication of the deacetylase sirtuin-1 on synovial angiogenesis and persistence of experimental arthritis. Ann. Rheum. Dis. 79, 891–900 (2020).32381568 10.1136/annrheumdis-2020-217377
14. Gründemann D Discovery of the ergothioneine transporter Proc. Natl. Acad. Sci. U.S.A. 2005 102 5256 5261 10.1073/pnas.0408624102 15795384
Gründemann, D. et al. Discovery of the ergothioneine transporter. Proc. Natl. Acad. Sci. U.S.A. 102, 5256–5261 (2005).15795384 10.1073/pnas.0408624102
15. Tokuhiro S An intronic SNP in a RUNX1 binding site of SLC22A4, encoding an organic cation transporter, is associated with rheumatoid arthritis Nat. Genet. 2003 35 341 348 10.1038/ng1267 14608356
Tokuhiro, S. et al. An intronic SNP in a RUNX1 binding site of SLC22A4, encoding an organic cation transporter, is associated with rheumatoid arthritis. Nat. Genet. 35, 341–348 (2003).14608356 10.1038/ng1267
16. Lee YH Meta-analysis of SLC22A4 and RUNX1 polymorphisms : Associations with rheumatoid arthritis susceptibility Z. Rheumatol. 2015 74 351 358 10.1007/s00393-014-1447-3 24972750
Lee, Y. H. et al. Meta-analysis of SLC22A4 and RUNX1 polymorphisms : Associations with rheumatoid arthritis susceptibility. Z. Rheumatol. 74, 351–358 (2015).24972750 10.1007/s00393-014-1447-3
17. Taubert D Association of rheumatoid arthritis with ergothioneine levels in red blood cells: A case control study J. Rheumatol. 2006 33 2139 2145 17086603
Taubert, D. et al. Association of rheumatoid arthritis with ergothioneine levels in red blood cells: A case control study. J. Rheumatol. 33, 2139–2145 (2006).17086603
18. Darghouth D Pathophysiology of sickle cell disease is mirrored by the red blood cell metabolome Blood 2011 117 e57 66 10.1182/blood-2010-07-299636 21135259
Darghouth, D. et al. Pathophysiology of sickle cell disease is mirrored by the red blood cell metabolome. Blood 117, e57-66 (2011).21135259 10.1182/blood-2010-07-299636
19. R Core Team (2021). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/
