Skip to main content

ORIGINAL RESEARCH article

Front. Med., 10 November 2022
Sec. Rheumatology
Volume 9 - 2022 | https://doi.org/10.3389/fmed.2022.1011734

Effects of tofacitinib therapy on arginine and methionine metabolites in association with vascular pathophysiology in rheumatoid arthritis: A metabolomic approach

Boglárka Soós1† Attila Hamar1† Anita Pusztai1 Monika Czókolyová1 Edit Végh1 Szilvia Szamosi1 Zsófia Pethő1 Katalin Gulyás1 György Kerekes2 Sándor Szántó1,3 Gabriella Szűcs1 Uwe Christians4 Jelena Klawitter4 Tamás Seres4‡ Zoltán Szekanecz1*‡
  • 1Department of Rheumatology, Faculty of Medicine, University of Debrecen, Debrecen, Hungary
  • 2Intensive Care Unit, Department of Internal Medicine, Faculty of Medicine, University of Debrecen, Debrecen, Hungary
  • 3Department of Sports Medicine, Faculty of Medicine, University of Debrecen, Debrecen, Hungary
  • 4Department of Anesthesiology, University of Colorado Anschutz Medical Campus, Aurora, CO, United States

Introduction: Rheumatoid arthritis (RA) has been associated with changes in lipid, arginine and NO metabolism with increased cardiovascular (CV) risk. The aim of this study is to evaluate the effect of tofacitinib, a Janus kinase (JAK) inhibitor, on arginine and methionine metabolism in correlation with inflammation, functional and pathological vascular changes during one-year treatment of patients with RA.

Materials and methods: Thirty RA patients with active disease were treated with either 5 mg bid or 10 mg bid tofacitinib for 12 months. We determined DAS28, CRP, IgM rheumatoid factor (RF) and anti-cyclic citrullinated peptide (CCP) levels. We assessed brachial artery flow-mediated vasodilation (FMD), carotid intima-media thickness (IMT) and pulse-wave velocity (PWV) by ultrasound at baseline and after 6 and 12 months. We also determined plasma L-arginine, L-citrulline, L-ornithine, inducible nitric oxide synthase (iNOS), asymmetric (ADMA) and symmetric dimethylarginine (SDMA), L-N-monomethyl-arginine (L-NMMA), cysteine, homocysteine, and methionine levels at these time points.

Results: Twenty-six patients (13 on each arm) completed the study. CRP, ESR and DAS28 decreased significantly during one-year treatment with tofacitinib. Arginine and ADMA showed a negative univariate correlation with CRP but not with FMD, PWV or IMT. Tofacitinib at 10 mg bid significantly increased L-arginine, L-ornithine, iNOS and methionine levels after 12 months. ADMA and SDMA levels did not change in our study. Methionine showed negative correlation with FMD at baseline and positive correlation with PWV after 12 months. No change was observed in FMD and PWV but a significant increase was measured in IMT at 6 and 12 months. Multivariate analysis indicated variable correlations of L-arginine, L-citrulline, ADMA, L-NMMA, homocysteine and methionine with DAS28, CRP, ESR and RF but not with anti-CCP after one-year treatment. With respect to vascular pathophysiology, only PWV and methionine correlated with each other.

Conclusion: One-year tofacitinib treatment suppressed systemic inflammation and improved functional status in RA. FMD, PWV have not been affected by one-year tofacitinib treatment., while IMT increased further despite treatment. Increased arginine and methionine might contribute to the anti-inflammatory effects of tofacitinib. Increased arginine availability with no changing ADMA may protect FMD and PWV from deterioration. The increase of IMT in the anti-inflammatory environment cannot be explained by arginine or methionine metabolism in this study.

Introduction

Rheumatoid arthritis (RA) is an autoimmune inflammatory rheumatic and musculoskeletal disease (RMD), which causes progressive deformation of different joints (1). Protein citrullination in the RA synovium is the first connection of the disease to arginine metabolism (25). At the joint level, there are known interactions between fibroblast-like synoviocytes (FLS), macrophages and helper T (TH)-, B- and plasma cells (6). In this environment, pro-inflammatory cytokines including tumor necrosis factor α (TNF-α), interleukin 1β (IL-1β) and IL-6 are abundantly produced (6, 7). This pro-inflammatory milieu promotes monocytes to differentiate to M1 type macrophages (6, 8). These macrophages express inducible nitric oxide synthase (iNOS) that stimulates the production of pro-inflammatory cytokines, NO and ONOO radicals at cytotoxic levels (9, 10). High expression of iNOS convert arginine into NO and citrulline (Figure 1) (10, 11). Synovial tissue damage initiates the differentiation of monocytes into M2 type macrophages (8), which cells produce large amount of arginase 1 and convert arginine to ornithine and urea (Figure 1) (10, 12, 13). Ornithine is converted to polyamines and proline promoting cell division and growth as well as collagen synthesis (11, 13). M2 macrophages produce anti-inflammatory cytokines and are involved in wound healing and tissue regeneration at the site of inflammation (8). The chronic inflammatory process is changing in time depending on which macrophages and cytokines are present and how the fibroblast-like synoviocytes (FLS) are activated (68). Pro-inflammatory cytokines can induce high amount of NO production with decreased arginase activation in human FLS (1416). At the same time synovial M2 type macrophages show no increase in NO production but significant increase in arginase activity in the presence of pro-inflammatory cytokines (10, 12). The shift in the balance between iNOS and arginase 1 in the inflammatory process might shift the arginine metabolism either to citrulline or ornithine over-production, respectively (Figure 1) (10, 12).

FIGURE 1
www.frontiersin.org

Figure 1. Metabolites of the urea cycle. ADMA, asymmetric dimethylarginine; L-NMMA, L-N-monomethyl-arginine; NO, nitric oxide; NOS, nitric oxide synthase; SDMA, symmetric dimethylarginine.

Rheumatoid arthritis (RA), similarly to diabetes, obesity and hypertension, is associated with chronic systemic inflammation leading to endothelial dysfunction, accelerated inflammatory atherosclerosis manifesting in high risk of cardiovascular diseases (CVD) (1720). In addition to coronary artery disease, endothelial dysfunction may cause myocardial hypertrophy and interstitial fibrosis manifesting in diastolic dysfunction and heart failure with preserved ejection fraction (EF) as also observed in diabetes, obesity and hypertension (21). The unique characteristic of CVD in RA is that there is higher risk for CVD at lower levels of cholesterol but high CRP levels. This lipid paradox may represent a special form of endothelial dysfunction where pro-inflammatory cytokines maintain the process of atherosclerosis (22, 23). Arginine is methylated during post-translational modification of different kind of proteins especially in the nucleus by protein methyl transferase enzymes (PRMT) (11). Three types of modifications are present on the proteins: N-monomethyl-L-arginine (L-NMMA), asymmetric (ADMA) and symmetric dimethyl arginine (SDMA) (11, 24). During protein turnover in the proteolysis phase the free L-NMMA and ADMA are metabolized to citrulline and dimethylamine by the enzyme dimethylarginine dimethylaminohydrolase (DDAH) (Figure 1) (11). Free SDMA is eliminated by the kidneys (25). L-NMMA and ADMA are competitive inhibitors of arginine at the NOS molecules causing decreased NO production and NOS uncoupling which turns NOS to produce superoxide anion instead of NO (11). ADMA and SDMA production has been associated with systemic inflammatory states (2628). In arthritides, serum ADMA level is higher than normal and is associated with endothelial dysfunction (2830). ADMA may cause endothelial dysfunction by disruption of the physiological endothelial NO production (11, 2830). Physiological NO production inhibits platelet aggregation, tissue factor expression, expression of adhesion molecules, as well as smooth muscle cell (SMC) contraction, proliferation and migration (11). These physiological roles keep the endothelial layer intact and inhibit atherosclerosis or thrombus formation (11).

Non-invasive, ultrasound-based methods are able to evaluate endothelial function, such as flow mediated dilation of brachial artery (FMD) (19, 31). Endothelial dysfunction may enhance media thickening via SMC hypertrophy and proliferation, as well as media fibrosis leading to increasing arterial stiffness. This can be evaluated by aortic pulse wave velocity (PWV) changes (31). Similarly, overt atherosclerosis is reflected by increasing carotid intima media thickness (IMT) and the development of carotid plaques, which can be detected by B-mode ultrasound (19, 31).

Methionine is the source of methylation by PMRT enzymes (Figure 2) (32). Homocysteine level is suppressed by cysteine synthesis, which enhance glutathione (GSH) synthesis and protection against oxidative stress (33). Increased methionine intake was associated with increased plasma concentration of homocysteine and reduced FMD in healthy human subjects (34, 35). The mechanism of inhibition of endothelial function by methionine load is not clear. Although the role of elevated ADMA level or oxidative stress were hypothesized they were not supported uniformly by various studies (3436).

FIGURE 2
www.frontiersin.org

Figure 2. Metabolites of the methionine and folate cycles. DHF, dihydrofolate; DHFR, dihydrofolate reductase; MTHF, methylene-tetradhydrofolate; MTX, methotrexate; THF, tetrahydrofolate.

Methionine load with increasing plasma homocysteine concentration did not alter PWV although it changes aortic distensibility in healthy individuals (37). Population based interventional study using vitamin B12 and folic acid in hyperhomocysteinemic patients did not show any effect of decreasing homocysteine level on PWV or IMT (38).

Targeted therapies including biologic and targeted synthetic disease-modifying drugs dampen systemic inflammation, RA disease activity and they may also have beneficial effects on CV outcomes (17, 39, 40). As of today, four Janus kinase (JAK) inhibitors, tofacitinib, baricitinib, upadacitinib and filgotinib have been approved for the treatment of RA (41). JAK inhibition has been associated with elevation of lipids, possibly due to the lipid paradox described above (17, 42, 43). We have recently reported that tofacitinib dampened aortic wall inflammation by PET/CT (44). There has been only one study assessing IMT in tofacitinib-treated patients (45).

To our best knowledge there have been no other studies on the effects of JAK inhibitors on vascular pathophysiology measured by FMD, PWV and IMT in parallel with arginine and methionine metabolism. No previous studies evaluated arginine metabolism as an indicator of change in inflammation or endothelial dysfunction.

We conducted a prospective one-year study in order to assess the effects of tofacitinib on inflammation and functional status as well as FMD, PWV and IMT in RA patients. We propose that one-year tofacitinib treatment will, as a JAK-STAT inhibitor decrease pro-inflammatory cytokine levels, decrease inflammation and improve functional status with signals from methionine and arginine metabolism such as citrulline and ornithine levels. Improving inflammation may improve vascular status parallel with indicators of endothelial dysfunction in arginine metabolism like L-NMMA or ADMA levels.

Patients and methods

Patients and study design

Thirty patients with active RA were recruited for this tofacitinib interventional study. Patient characteristics are presented in Table 1. Inclusion criteria included definitive diagnosis of RA according to the 2010 European League Against Rheumatism (EULAR)/American College of Rheumatology (ACR) classification criteria for RA (1); moderate-high disease activity (DAS28 > 3.2) at baseline and clinical indication of targeted therapy. Patients were either naïve to any targeted therapies (n = 16) or initiated tofacitinib after stopping a biologic followed by an appropriate washout period (n = 14). Exclusion criteria included inflammatory diseases other than RA, acute/recent infection, standard contraindications to JAK inhibition, uncontrolled CV disease or hypertension, chronic renal or liver failure and malignancy within 10 years.

TABLE 1
www.frontiersin.org

Table 1. Patient characteristics.

The 30 enrolled patients were randomly assigned in a 1:1 ratio to either 5mg or 10mg tofacitinib twice daily (bid) treatment arms. All patients received tofacitinib in combination with either methotrexate (MTX) (n = 23) or leflunomide (n = 7). MTX and leflunomide had been taken in stable dose at least one year prior to the present study. No dose changes of these DMARDs were allowed throughout the course of the study. Although most patients may have received corticosteroids prior to the study, none of the patients had been on corticosteroids for at least 3 months prior to and during the study.

Clinical assessments were performed at baseline, and after 6 and 12 months of therapy. Four patients (2 on each arm) completed the 6-month follow-up but did not complete the one-year treatment. Twenty-six patients completed the one-year treatment period and were included in the data analysis.

The study was approved by the Hungarian Scientific Research Council Ethical Committee (approval No. 56953-0/2015-EKL). Written informed consent was obtained from each patient and assessments were carried out according to the Declaration of Helsinki and its amendments.

Clinical assessment

First, a detailed medical history was taken. We inquired about history of CVD, as well as current smoking, experience of chest pain resembling angina pectoris, hypertension and diabetes mellitus during the last 2 years prior to the start of this study by a questionnaire (Table 1). Further clinical assessments including physical examination were performed at baseline, and after 3, 6 and 12 months of tofacitinib therapy.

Laboratory measurements and assessment of disease activity

Blood samples were drawn from fasting patients in the morning into ethylene-diamine-tetraacetate (EDTA)-treated tubes and were immediately processed, aliquoted and stored at –70°C until use. Blood samples were taken at baseline, as well as after 6 and 12 months of tofacitinib treatment.

Serum high sensitivity C reactive protein (hsCRP; normal: ≤5 mg/l) and IgM rheumatoid factor (RF; normal: ≤50 IU/ml) were measured by quantitative nephelometry (Cobas Mira Plus, Roche Diagnostics, Basel, Switzerland), using CRP and RF reagents (both Dialab Ltd, Budapest, Hungary). ACPA (CCP) autoantibodies were detected in serum samples using a second generation Immunoscan-RA CCP2 ELISA test (Euro Diagnostica, Malmö, Sweden; normal: ≤25 IU/ml). The assay was performed according to the manufacturer’s instructions.

Disease activity of RA was calculated as DAS28-CRP (3 variables) (46).

Assessment of vascular physiology by ultrasound

Brachial artery FMD was assessed as described before (47). In brief, ultrasound examination was performed on the right arm using 10 MHz linear array transducer (ultrasound system: HP Sonos 5500, Hewlett Packard, Palo Alto, CA, USA) by a single trained sonographer after 30 min resting in a temperature-controlled room (basal value for FMD). A B-mode longitudinal section was obtained of the brachial artery above the antecubital fossa. In order to assess FMD, reactive hyperemia was induced by release of a pneumatic cuff around the forearm inflated to suprasystolic pressure for 4.5 min. After deflation the maximal flow velocity and the arterial diameter was continuously recorded for 90 s. Flow velocities, the baseline diameter, as well as FMD were ECG gated and detected offline. FMD values were expressed as% change from baseline (resting) value.

The IMT measurements were carried out as described before (47). Briefly, a duplex ultrasound system (HP Sonos 5500, 10 MHz linear array transducer) was used to assess the common carotid arteries by a single observer. Longitudinal high-resolution B-mode ultrasound scan were employed over both right and left common carotid arteries and were R-synchronized and recorded. The offline measurements were performed 1 cm proximal to the carotid bulb in the far wall. IMT was defined as the distance between the first and second echogenic lines from the lumen taking the average of 10 measurements on both sides. IMT values were expressed in mm.

With respect to arterial stiffness, PWV was calculated automatically by a TensioClinic arteriograph system (Tensiomed Ltd, Budapest, Hungary) as the quotient of the distance between the jugular fossa and symphysis as described before (47). If an artery is elastic, PWV is low. With decreased arterial elasticity, PWV rises. The arteriograph assesses this parameter from the oscillometric data obtained from the 35 mmHg suprasystolic pressure of the brachial artery. In order to obtain reproducible results, the patient had to rest in a supine position for at least 10 min before the assessment in a quiet room. PWV is expressed in m/s. Reproducibility of the three techniques expressed in intraclass correlation is included in the Statistical analysis section.

Assessment of metabolites in the urea and methionine cycles

L-arginine, L-citrulline, L-NMMA, ADMA, SDMA, cysteine, homocysteine, methionine and ornithine were quantified using a validated HPLC-MS/MS assay (48, 49). In brief, 50 μl of internal standard containing solution (50 μM d7-ADMA, d7-arginine, d4-cysteine, d6-citrulline, d8-homocystine, d3-methionine and d6-ornithine, all in HPLC water) and 40 μl of 500 mM DTT solution were added to 100 μl serum. For protein precipitation, 400 μl of acetonitrile containing 0.1% trifluoric acid was added to the sample. The sample was vortexed for 5 min, centrifuged for 20 min at 16,000 g, and transferred into a HPLC vial. Ten μl of the supernatant was injected onto a 4.6 × 12.5-mm guard column (Eclipse XDB-C8, 5 μm, Agilent Technologies, Palo Alto, CA, USA) in line with a 3.0 × 150-mm analytical column (RP-Amide, 3.5 μm, Supelco, St. Louis, MO, USA). The API5000 mass spectrometer (AB Sciex, Concord, ON, Canada) was run in the positive electrospray ionization mode (ESI) using multiple reaction monitoring (MRM) (48, 49).

Measurement of human inducible nitric oxide synthase

Inducible nitric oxide synthase (iNOS) was measured in human plasma by sandwich-enzyme-linked immunosorbent assay (ELISA) according to the manufacturer’s instructions (Novus Biologicals, CO, USA). The micro ELISA plate provided in the kit had been pre-coated with an antibody specific to human iNOS. Fifty micro liter of standard (serially diluted from 1000 to 15.6 pg/ml) or test sample was added to the plate wells and incubated for 2 h. This was followed by addition of 100 μl biotinylated detection antibody specific for human iNOS for 1 h and 100 μl of avidin-horseradish peroxidase (HRP) conjugate for 0.5 h. The substrate solution (90 μl) was added to each well for 20 min. The enzyme-substrate reaction was terminated by the addition of stop solution. Optical density (OD) was measured at a wavelength of 450 nm, Concentration of iNOS in the samples was calculated by comparing the OD of the samples to the standard curve. All incubations were done on a shaker at room temperature and plates were washed with wash buffer between incubation periods.

Statistical analysis

Statistical analysis was performed using SPSS version 22.0 (IBM, Armonk, NY, USA) software. Data are expressed as the mean ± SD for continuous variables and percentages for categorical variables. The distribution of continuous variables was evaluated by Kolmogorov-Smirnov test. Continuous variables were evaluated by paired two-tailed t-test and Wilcoxon test. Nominal variables were compared between groups using the chi-squared or Fisher’s exact test, as appropriate. Correlations were determined by Pearson’s analysis. Univariate and multivariate regression analysis using the stepwise method were applied to investigate independent associations between metabolites of the urea or methionine cycle (dependent variables) and clinical (disease activity), laboratory (CRP, ESR) or vascular physiology parameters (independent variables). The β standardized linear coefficients showing linear correlations between two parameters were determined. The B (+95% CI) regression coefficient indicated independent associations between dependent and independent variables during changes. General linear model (GLM) repeated measures analysis of variance (RM-ANOVA) was performed in order to determine the additional effects of multiple parameters including therapy on 6- or 12-month changes of metabolite levels. In this analysis, partial η2 is given as indicator of effect size, with values of 0.01 suggesting small, 0.06 medium and 0.14 large effects. The power was estimated using the G*-Power 3 software (50). P < 0.05 were considered significant.

The reliability of the vascular ultrasound measurements was tested by inter-item correlation and intraclass correlation (ICC). With respect to the FMD, IMT and PWV tests, ICC = 0.470; F-test value: 1.887; p = 0.001. The power was estimated using the G*-Power software (50). P < 0.05 were considered significant.

Results

Characteristics of patients

Patient characteristics are seen in Table 1. Altogether 6 patients (3-3 on each arm) had a positive CV history. All these patients had coronary artery disease, which was properly treated and stable. A total of 14 patients had hypertension, 2 had diabetes mellitus and 7 had been current smokers at the time of inclusion. Hypertension and diabetes mellitus were also well-controlled in those patients (Table 1).

Clinical response to tofacitinib therapy

Eventually a total of 4 patients, 2-2 each treatment arms, dropped out after 6 months of treatment but before the end of the study. Out of the 4 patients, 2 had inefficacy, one had significantly elevated transaminases and one moved abroad. Thus, 13-13 patients on each arm completed the study and were eligible for further data analysis (Table 1).

Tofacitinib treatment significantly decreased DAS28 after 6 months (3.31 ± 0.91; p < 0.001) and 12 months of treatment (3.32 ± 1.12; p < 0.001) compared to baseline (5.05 ± 0.77). CRP decreased from 14.8 ± 14.9 mg/l at baseline to 5.3 ± 5.3 mg/l after 6 months (p < 0.001) and 7.4 ± 7.7 mg/l after 12 months (p = 0.001). Similar observations were made in the 5 mg bid and 10 mg bid subsets (data not shown).

Effects of tofacitinib on vascular pathophysiology

Carotid IMT significantly increased after 6 months (0.56 ± 0.12 mm; p = 0.05) and 12 months (0.59 ± 0.14 mm; p = 0.002) compared to baseline (0.53 ± 0.11 mm) in the full cohort. In the 5mg bid subset, there was no difference in IMT between 6 months and baseline. However, after 12 months, IMT significantly increased vs. baseline (p = 0.007). In the 10 mg bid subset, there were no significant differences in IMT after 6 or 12 months compared to baseline (Figure 3). In the total cohort, as well as in the 5 mg bid and 10 mg bid subsets, neither FMD nor PWV showed any significant changes over time (Figure 3).

FIGURE 3
www.frontiersin.org

Figure 3. Effects of one-year tofacitinib therapy on (A) flow-mediated vasodilation (FMD). (B) carotid intima-media thickness (IMT) and (C) arterial pulse-wave velocity (PWV). In the full cohort and in the 5 mg bid group there was a progression in IMT after 12 months. In contrast, 10 mg bid tofacitinib was able to halt IMT progression. Both doses of tofacitinib slowed down FMD and PWV progression.

Effects of tofacitinib therapy on l-arginine metabolism

In the full cohort, tofacitinib significantly increased L-arginine levels by 6 (p = 0.004) and 12 months (p = 0.043). Similar pattern was observed in the 10 mg bid subset (p = 0.004 and p = 0.013, respectively) but not in the 5 mg bid subset (Figure 4A).

FIGURE 4
www.frontiersin.org

Figure 4. Effects of one-year tofacitinib therapy on metabolites in the urea cycle. (A) L-arginine levels increased in the full cohort and in the 10 mg bid subset after 6 and 12 months. (B) L-citrulline production transiently increased in the full cohort and the 10 mg bid subset after 6 months. (C) L-ornithine levels increased overtime in the full cohort and in the 10 mg bid subset after 6 and 12 months. (D) iNOS release increased in the full cohort and in the 10 mg bid subset after 12 months. (E) ADMA and (F) SDMA levels did not change overtime. (G) L-NMMA production transiently increased in the full cohort and the 10 mg bid subset after 6 months.

L-citrulline levels showed a transient increase in the full cohort. L-citrulline increased after 6 months vs. baseline (p = 0.006) but significantly decreased after 12 months compared to 6 months (p = 0.023). Tofacitinib 10 mg bid also transiently increased L-citrulline after 6 months (p = 0.018). No change was observed in the 5 mg bid subset (Figure 4B).

In the full cohort, L-ornithine also increased after 6 months vs. baseline (p = 0.025). Then, L-ornithine levels showed a tendency to remain higher after 12 months (p = 0.119). There was no significant difference between L-ornithine levels after 6 and 12 months. In the 10 mg bid subset, L-ornithine levels were significantly higher after both 6 (p = 0.018) and 12 months (p = 0.020) vs. baseline. Again, L-ornithine production did not change over time in the 5 mg bid subset (Figure 4C).

We also determined L-arginine/L-citrulline and L-arginine/L-ornithine ratios at 6 and 12 months. These ratios reflect the above described balances between pro- and anti-inflammatory mechanisms that may, at least in part, lead to clinical improvement upon tofacitinib treatment. Both L-arginine/L-citrulline (p = 0.007) and L-arginine/L-ornithine ratio (p = 0.035) significantly but transiently decreased after 6 months. Both ratios returned to the baseline level after 12 months (data not shown).

Inducible nitric oxide synthase (iNOS) levels after 12 months were significantly higher compared to both baseline (p = 0.045) and 6 months (p = 0.020) in the full cohort. There was no difference between 6 months and baseline. Similar pattern was observed in the 10 mg bid subset (p = 0.047 and p = 0.043, respectively), but not in the 5 mg bid subset (Figure 4D).

Tofacitinib treatment did not change ADMA (Figure 4E) or SDMA levels (Figure 4F) over time. No changes were seen in the full cohort, 5 mg or 10 mg bid subsets.

Finally, L-NMMA levels also showed a transient increase upon tofacitinib therapy. In the full cohort, L-NMMA was significantly increased after 6 months (p = 0.008), but then levels dropped after 12 months compared to 6 months (p = 0.048). Similar pattern was observed in the 10 mg bid subset (p = 0.044 and p = 0.047, respectively). There was no change in L-NMMA in the 5 mg bid subset (Figure 4G).

Effect of tofacitinib on metabolites in the methionine cycle

In the full cohort, cysteine showed a “late” increase after 12 months compared to baseline (p = 0.028) and 6 months (p = 0.005). There was no difference between baseline and 6 months. This increase was not observed in the 5 mg bid or 10 mg bid subset (Figure 5A).

FIGURE 5
www.frontiersin.org

Figure 5. Effects of one-year tofacitinib therapy on metabolites in the methionine cycle. (A) Cysteine levels increased overtime in the full cohort after 12 months. (B) Homocysteine production decreased in the full cohort and the 10 mg bid subset after 12 months. (C) Methionine levels increased overtime in the full cohort and in the 10 mg bid subset after 6 and 12 months and also in the 5 mg bid subset after 12 months.

Homocysteine levels did not change between 6 months and baseline, however, later significant decreases were observed between 6 and 12 months in the full cohort (p = 0.047), as well as in the 10 mg bid subset (p = 0.049). No changes in homocysteine levels were seen in the 5 mg bid subset (Figure 5B).

Finally, methionine levels increased after 6 months versus baseline in the full cohort (p = 0.005) and in the 10 mg bid subset (p < 0.001). In addition, methionine levels increased after 12 months compared to baseline in the full cohort (p < 0.001), as well as in the 5 mg bid (p = 0.002) and the 10 mg bid arms (p < 0.001) (Figure 5C).

Correlations between levels of metabolites with each other and with other parameters

Table 2 shows uni- and multivariate regression analysis of metabolites as dependent variables. In general, the univariate analysis suggests that baseline and 12-month L-arginine, L-citrulline, L-NMMA, ADMA, homocysteine and methionine levels variably correlate with disease activity/inflammation (DAS28, CRP, ESR) or seropositivity (RF, CCP). The correlations with DAS28 are positive, while those with CRP, ESR and CCP are negative. RF positively correlated with L-arginine and L-citrulline, but negatively with ADMA, L-NMMA and homocysteine. Baseline disease duration may also determine L-NMMA after 12 months.

TABLE 2
www.frontiersin.org

Table 2. Univariable and multivariable analysis of determinants of metabolites as dependent variables.

There were sporadic correlations between vascular parameters and arginine or methionine metabolites. FMD negatively correlated with arginine and methionine at baseline. IMT positively correlated with ADMA and L-NMMA at 6 months. PWV positively correlated with homocysteine at baseline and methionine at 12 months.

The multivariate analysis confirmed variable correlations of L-arginine, L-citrulline, ADMA, L-NMMA, homocysteine and methionine after one-year treatment with DAS28, CRP, ESR and RF, but not with anti-CCP. Interestingly, most correlations, except for two, were observed regarding 12-month levels of the urea and methionine cycle metabolites. With respect to vascular pathophysiology, only PWV and methionine after 12 months correlated with each other (Table 2).

Repeated measures analysis of variance (RM-ANOVA) was performed in order to determine the effects of tofacitinib treatment in combination with other parameters on changes in the levels of metabolites over time. Six-month change (increase) in L-citrulline levels was significantly determined by tofacitinib treatment together with either lower ESR or higher RF at baseline. After 6 months, L-citrulline decreased to the level of baseline. Moreover, 12-month changes in homocysteine levels correlated with treatment in combination with lower ESR at baseline (Table 3).

TABLE 3
www.frontiersin.org

Table 3. Significant results of general linear model (GLM) repeated measures analysis of variance (RM-ANOVA) test determining the effects of treatment and other independent variables on 12-month changes in the levels of metabolites as dependent variables.

Discussion

Rheumatoid arthritis (RA) is characterized by systemic inflammation and high risk for CVD. In the present study, we have followed endothelial function (FMD), vascular stiffness (PWV) and intima-media thickness (IMT) during one-year treatment with tofacitinib. We applied a metabolomic approach in order to evaluate the dose dependent effect of tofacitinib treatment on arginine and methionine metabolism in association with the level of inflammation and vascular pathology. The two amino acids are interacting with each other as methionine is the source of methylation of arginine on proteins and alters homocysteine and cysteine synthesis (5153).

One-year tofacitinib treatment decreased CRP, ESR and DAS28 showing significant antiinflammatory effects and improving functional status.

Tofacitinib has significantly increased arginine and methionine levels during the one-year treatment. Inflammation with increased Il-6 level decreases muscle protein synthesis and increases amino-acid uptake of other organs decreasing the amino-acid level (54). As tofacitinib, among other cytokines, suppresses IL-6 production it might potentially reverse the effect described above, promoting an increase in plasma level of amino acids.

Citrulline and ornithine level may represent the balance between inflammatory and anti-inflammatory macrophages or other cell types. Citrulline and ornithine levels increased significantly at 6 months and returned to baseline in the whole group. The increasing citrulline may represent flaring inflammation with producing destructive NO production by M1 type macrophages or other inflammatory cells. The decrease of CRP, ESR and DAS28 at 6 months represented inhibition of inflammation, thus the increased NO production most likely belonged to endothelial cells. The increased ornithine level, however, may represent increased arginase activity in M2 type macrophages in parallel with the anti-inflammatory effect of tofacitinib. We do not have enough data to explain why citrulline and ornithine levels returned to baseline at 12 months and what is the clinical significance of the significantly higher level of ornithine in the 10mg BID tofacitinib group at 12 months. There has been only one cross-sectional study assessing metabolites of arginine in RA (27). In that study, decreased levels of plasma arginine and citrulline, elevated levels of L-ornithine and elevated arginase activity was found in RA patients compared to controls. In that patient group the CRP and ESR levels were almost normal so those patients represented marked anti-inflammatory macrophage activity with increased arginase level. In our patient group the CRP and ESR values were high even after tofacitinib treatment, which might explain we did not see decreasing arginine and major shift in the ornithine and citrulline ratio.

Inducible nitric oxide synthase (iNOS) level was significantly higher at 12 months in the whole ant the 10mg BID groups. In an in vitro observation of the effect of tofacitinib on human dendritic cells and macrophages, tofacitinib induced M1 macrophage phenotypes with iNOS expression and IL-12 and IL-23 production (55). However, it should be considered that tofacitinib was combined with methotrexate in our study. Methotrexate inhibits M1 and promotes M2 macrophage differentiation through adenosine receptors. At the same time it is JAK1/JAK2 inhibitor also (56). Our study suggests that the interaction between methotrexate and tofacitinib might shift to M2 macrophage activity with slight increase in arginase activity and ornithine level and inhibition of the cells with iNOS expression especially at higher dose of tofacitinib.

Similarly to arginine, methionine level significantly increased at 6 and 12 months. Methionine, as an essential amino acid, is highly involved in the methylation of different proteins and nucleic acid (32). Methylations could modulate T cell functions, impairing Th1/Th2 cytokines release, and decreasing T cell proliferation and activation. Furthermore, epigenetic regulations play a significant role in gene expression affecting immunocyte function and signaling pathways (32). Global DNA hypomethylation in RA synovial fibroblasts contributes to their intrinsic activation (57). In our study the increased methionine availability during tofacitinib treatment may provide at least appropriate methylation status mediating anti-inflammatory effects. However, homocysteine level did not change significantly from the baseline level with increasing methionine indicating the same level of methylation reactions. Following methionine and homocysteine levels are not conclusive to describe the change in inflammation during tofacitinib treatment. There was an increase in cysteine and a decrease in homocysteine synthesis at 12 months. Increased cysteine synthesis initiates more glutathione production, which enhance cellular protection against oxidative stress (58).

Tofacitinib therapy did not alter endothelial function and vascular stiffness measured by FMD and PWV, respectively. In most studies, impairment of FMD (19, 20, 31, 47, 59, 60) and increasing PWV (31, 47, 6163) were observed in RA. Our study suggests that one year tofacitinib therapy was able to prevent further worsening of endothelial function and arterial stiffness, respectively. In contrast, carotid atherosclerosis (IMT) further deteriorated overtime despite tofacitinib therapy. In numerous studies, IMT was between 0.63 and 0.73 mm in RA compared to 0.54-0.62 mm in healthy controls (19, 20, 31, 6466). In our present study, the mean baseline IMT was only 0.53 mm (range: 0.40–0.66 mm), which is similar to healthy subjects in other studies. Although, IMT progressed despite tofacitinib treatment, it does not mean that tofacitinib would promote atherosclerosis. One year of tofacitinib treatment may not be long enough to suppress atherosclerosis. Very recently, in the ORAL Surveillance study, tofacitinib treatment in comparison to anti-TNF therapy was associated with increase risk of major cardiovascular events (67). That study, on contrast to ours, included RA patients with at least one cardiovascular risk at baseline and the average age of that patient group was 10-year older, than ours. The latest integrated safety analysis did not find increased cardiovascular risk upon tofacitinib therapy (68).

One-year tofacitinib treatment decreased CRP, ESR and DAS28 showing significant antinflammatory effects but it did not show major effect on vascular function. There has been only one study where IMT did not significantly change during one year tofacitinib treatment (10mg/day) in a methotrexate resistant patient group. In this study cardio-ankle vascular index and augmentation index decreased significantly after one-year treatment (45). The notable difference between their study and ours, that the CRP level was twice higher in our study at 1 year. This might explain why only 20mg/day tofacitinib were effective in our study to inhibit significant increase in IMT. We did not find any other studies where tofacitinib or any other JAK inhibitors were studied in relation to FMD or PWV. Other investigators also did not find association between vascular function or morphology and ADMA or SDMA in RA (69, 70).

Arginine, citrulline, ornithine and global arginine bioavailability ratio (GABR) [arginine/(citrulline + ornithine)] were measured in patients with and without coronary artery disease (CAD). The prevalence of CAD and risk of death or major adverse cardiovascular events (MACE) were evaluated at 3 years with decreasing quartiles of arginine and GABR and increasing quartiles of citrulline and ornithine. The best predictor for CAD, MACE and death was GABR and >1.46 characterized the smallest risk. This study suggested that arginine availability is an important predictor in vascular pathology (71). In our study, the calculated GABR was 1.47 ± 0.41, 1.31 ± 0.36, 1.49 ± 0.38 at 0, 6, and 12 months, respectively. These numbers put our patient group in low-risk categories for CAD, death and MACE. The individual concentrations of arginine, citrulline and ornithine represented low risk categories also during tofacitinib treatment. This may explain that tofacitinib has not altered endothelial function (FMD) and vascular stiffness (PWV) during one-year treatment. At higher concentration tofacitinib inhibited carotid atherosclerosis maybe because that patient group had the highest arginine bioavailability. We could not find any correlation between the vascular parameters and arginine or methionine metabolism, most likely because the measured parameters were in the normal range. This study was not controlled, however, as a historical control, previous studies indicated that impaired FMD and increased IMT and PWV have been associated with RA (17, 18, 31, 39, 47, 72) and biologics including TNF blockers and others may, at least transiently, improve endothelial function, atherosclerosis and arterial stiffness or at least slow down the progression of these parameters [reviewed in Szekanecz et al. (39)].

Methylated arginine derivatives are enriched in the plasma after proteolysis (7375) and the DDAH enzyme transforms ADMA and L-NMMA to L-citrulline (73). Increased circulating ADMA levels have been associated with RA and spondyloarthritides (28, 29). Similarly to increasing citrulline and ornithine, increasing ADMA level increased the risk for MACE in a patient groups of peripheral artery disease or in CAD (30, 76, 77). In patients with peripheral artery disease ADMA level < 0.68 μmol/l represented the lowest quartile for hazard ratio for MACE. In our study the ADMA level was between 0.6 and 0.64 μmol/l range representing a normal level and low risk for MACE. This ADMA level did not change significantly during the study period. The effects of targeted therapies on ADMA may be controversial (24, 78). TNF-α inhibits ADMA degradation by inhibition of DDAH and thus increases ADMA levels (78). Anti-TNF therapy decreased ADMA levels in a few (78) but, similarly to ours, not in most studies (7981). The L-arginine/ADMA ratio is important for modulation of NOS activity and it is a risk factor for atherosclerosis (82, 83). TNF-α inhibitors may increase L-arginine/ADMA ratios at 3 and 12 months of treatment (84). We observed a similar pattern in our present study applying tofacitinib therapy. Although the L-arginine/ADMA ratio did not change significantly in our study, the higher L-arginine levels along with unchanged levels of ADMA might contribute to the inhibition of carotid wall thickening at the higher dose of tofacitinib therapy. We have recently reported that tofacitinib attenuated vascular inflammation as determined by PET/CT (44). Based on our study the increased arginine level, high arginine availability for NO synthesis and high arginine/ADMA ratio may explain the attenuation of vascular inflammation during tofacitinib treatment.

Although, methionine showed a negative correlation with FMD at baseline, increased methionine level at 6 and 12 months did not change FMD significantly. The effect of chronic elevation of methionine on FMD is different from the results of previous studies where methionine load increased homocysteine level and decreased FMD in healthy individuals and patients with various diseases (3436). Elevated homocysteine and cysteine levels were associated with decreased FMD across menopausal stages in healthy women. In this study methionine showed a positive correlation with FMD (49). Homocysteine level is higher in RA compared to normal population (85). In our longitudinal study homocysteine level was higher than normal and there was no significant difference between baseline and 12-month values although the methionine level increased significantly. In conclusion, tofacitinib treatment increase methionine level without altering homocysteine level and endothelial function. There was a positive correlation between methionine and PWV at 12 month in our study. On the other hand, there was no significant change in PWV at 12 months compared to baseline values. In comparison with FMD, PWV is not changing with methionine load (37). However, multiple studies showed association between homocysteine and PWV in various diseases (8688). One year tofacitinib treatment did not alter PWV significantly in the presence of high homocysteine and methionine level suggesting not progressing vascular stiffness in our patient population.

Arginine, citrulline, ADMA, L-NMMA, homocysteine and methionine variably correlated with DAS28, CRP and ESR. Interestingly, these metabolites showed inverse correlations with CRP, a marker of systemic inflammation suggesting that the levels of some of these metabolites increase in parallel with attenuation of CRP levels by JAK inhibition. Some reports also found low levels of arginine and citrulline in severe infectious and inflammatory states associated with elevated CRP release (89, 90).

Our study has certain advantages and limitations. There have been no similar prospective studies assessing metabolites of the urea and methionine cycle in association with disease activity, autoantibodies and vascular pathophysiology in RA patients undergoing any targeted treatments. Possible limitations may include the relatively low number of RA patients and the lack of control group due to the fact that this is a self-controlled prospective study.

In conclusion, we propose that in our RA cohort, tofacitinib treatment, in parallel with its anti-inflammatory action, normalized the balance of amino acid metabolism and protected vascular function at least for one year. Although we did not assess M1 and M2 macrophages our results suggested that tofacitinib might shift the arginine metabolism to anti-inflammatory direction. The higher methionine level might also contribute to anti-inflammatory immune modulation. Endothelial function and arterial stiffness have not changed in one-year tofacitinib treatment partly due to the high arginine availability and not increasing ADMA level. Increasing methionine level did not associated with higher homocysteine level and did not alter FMD and PWV. The mechanism of increasing IMT during tofacitinib treatment needs further investigation.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving human participants were reviewed and approved by Hungarian Scientific Research Council Ethical Committee (approval No. 56953-0/2015-EKL). The patients/participants provided their written informed consent to participate in this study.

Author contributions

BS: study concept, patient recruitment and examination, and manuscript drafting. AH, EV, SSzam, ZP, KG, SSzán, and GS: patient recruitment and examination, and data collection. AP and MC: laboratory tests and data analysis. GK: vascular assessments and data analysis. UC and JK: metabolic assessments, data analysis, and interpretation. TS and ZS: senior PIs, supervisors, study conceptualization, and manuscript finalization. All authors contributed to the article and approved the submitted version.

Funding

This research was supported by the European Union and the State of Hungary and co-financed by the European Social Fund in the framework of TAMOP-4.2.4.A/2-11/1-2012-0001 “National Excellence Program” (ZS); by the European Union grant GINOP-2.3.2-15-2016-00050 (ZS) and by the WI188341 investigator-initiated research (IIR) grant obtained from Pfizer US (ZS). The funder (Pfizer US) was not involved in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Aletaha D, Neogi T, Silman AJ, Funovits J, Felson DT, Bingham CO III, et al. 2010 rheumatoid arthritis classification criteria: an American College of Rheumatology/European League Against Rheumatism collaborative initiative. Ann Rheum Dis. (2010) 69:1580–8.

Google Scholar

2. Wegner N, Lundberg K, Kinloch A, Fisher B, Malmstrom V, Feldmann M, et al. Autoimmunity to specific citrullinated proteins gives the first clues to the etiology of rheumatoid arthritis. Immunol Rev. (2010) 233:34–54. doi: 10.1111/j.0105-2896.2009.00850.x

PubMed Abstract | CrossRef Full Text | Google Scholar

3. De Rycke L, Nicholas AP, Cantaert T, Kruithof E, Echols JD, Vandekerckhove B, et al. Synovial intracellular citrullinated proteins colocalizing with peptidyl arginine deiminase as pathophysiologically relevant antigenic determinants of rheumatoid arthritis-specific humoral autoimmunity. Arthritis Rheum. (2005) 52:2323–30. doi: 10.1002/art.21220

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Lundberg K, Bengtsson C, Kharlamova N, Reed E, Jiang X, Kallberg H, et al. Genetic and environmental determinants for disease risk in subsets of rheumatoid arthritis defined by the anticitrullinated protein/peptide antibody fine specificity profile. Ann Rheum Dis. (2013) 72:652–8. doi: 10.1136/annrheumdis-2012-201484

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Dekkers J, Toes RE, Huizinga TW, van der Woude D. The role of anticitrullinated protein antibodies in the early stages of rheumatoid arthritis. Curr Opin Rheumatol. (2016) 28:275–81. doi: 10.1097/BOR.0000000000000277

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Szekanecz Z, Koch AE. Update on synovitis. Curr Rheumatol Rep. (2001) 3:53–63. doi: 10.1007/s11926-001-0051-0

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Szekanecz Z, Koch AE, Kunkel SL, Strieter RM. Cytokines in rheumatoid arthritis. Potential targets for pharmacological intervention. Drugs Aging. (1998) 12:377–90. doi: 10.2165/00002512-199812050-00004

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Tardito S, Martinelli G, Soldano S, Paolino S, Pacini G, Patane M, et al. Macrophage M1/M2 polarization and rheumatoid arthritis: A systematic review. Autoimmun Rev. (2019) 18:102397. doi: 10.1016/j.autrev.2019.102397

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Forstermann U, Sessa WC. Nitric oxide synthases: regulation and function. Eur Heart J. (2012) 33:829–37. doi: 10.1093/eurheartj/ehr304

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Rodriguez PC, Ochoa AC, Al-Khami AA. Arginine metabolism in myeloid cells shapes innate and adaptive immunity. Front Immunol. (2017) 8:93. doi: 10.3389/fimmu.2017.00093

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Wu G, Morris SM Jr. Arginine metabolism: nitric oxide and beyond. Biochem J. (1998) 336(Pt 1):1–17. doi: 10.1042/bj3360001

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Yang Z, Ming XF. Functions of arginase isoforms in macrophage inflammatory responses: impact on cardiovascular diseases and metabolic disorders. Front Immunol. (2014) 5:533. doi: 10.3389/fimmu.2014.00533

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Panfili E, Gerli R, Grohmann U, Pallotta MT. Amino acid metabolism in rheumatoid arthritis: Friend or Foe? Biomolecules. (2020) 10:1280. doi: 10.3390/biom10091280

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Caldwell RB, Toque HA, Narayanan SP, Caldwell RW. Arginase: an old enzyme with new tricks. Trends Pharmacol Sci. (2015) 36:395–405. doi: 10.1016/j.tips.2015.03.006

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Corraliza I, Moncada S. Increased expression of arginase II in patients with different forms of arthritis. Implications of the regulation of nitric oxide. J Rheumatol. (2002) 29:2261–5.

Google Scholar

16. Miyoshi F, Sato K, Mimura T. Changes in the pattern of cytokine production from peripheral blood mononuclear cells in patients with rheumatoid arthritis treated with infliximab and their relation to plasma arginase activity. Int J Rheum Dis. (2018) 21:1907–14. doi: 10.1111/1756-185X.12864

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Agca R, Heslinga SC, Rollefstad S, Heslinga M, McInnes IB, Peters MJ, et al. EULAR recommendations for cardiovascular disease risk management in patients with rheumatoid arthritis and other forms of inflammatory joint disorders: 2015/2016 update. Ann Rheum Dis. (2017) 76:17–28. doi: 10.1136/annrheumdis-2016-209775

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Nurmohamed MT, Heslinga M, Kitas GD. Cardiovascular comorbidity in rheumatic diseases. Nat Rev Rheumatol. (2015) 11:693–704. doi: 10.1038/nrrheum.2015.112

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Kerekes G, Szekanecz Z, Der H, Sandor Z, Lakos G, Muszbek L, et al. Endothelial dysfunction and atherosclerosis in rheumatoid arthritis: a multiparametric analysis using imaging techniques and laboratory markers of inflammation and autoimmunity. J Rheumatol. (2008) 35:398–406.

Google Scholar

20. Gonzalez-Gay MA, Gonzalez-Juanatey C, Vazquez-Rodriguez TR, Martin J, Llorca J. Endothelial dysfunction, carotid intima-media thickness, and accelerated atherosclerosis in rheumatoid arthritis. Semin Arthritis Rheum. (2008) 38:67–70. doi: 10.1016/j.semarthrit.2008.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Blyszczuk P, Szekanecz Z. Pathogenesis of ischaemic and non-ischaemic heart diseases in rheumatoid arthritis. RMD Open. (2020) 6:e001032. doi: 10.1136/rmdopen-2019-001032

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Myasoedova E, Crowson CS, Kremers HM, Roger VL, Fitz-Gibbon PD, Therneau TM, et al. Lipid paradox in rheumatoid arthritis: the impact of serum lipid measures and systemic inflammation on the risk of cardiovascular disease. Ann Rheum Dis. (2011) 70:482–7. doi: 10.1136/ard.2010.135871

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Kerekes G, Nurmohamed MT, Gonzalez-Gay MA, Seres I, Paragh G, Kardos Z, et al. Rheumatoid arthritis and metabolic syndrome. Nat Rev Rheumatol. (2014) 10:691–6. doi: 10.1038/nrrheum.2014.121

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Dowsett L, Higgins E, Alanazi S, Alshuwayer NA, Leiper FC, Leiper J. ADMA: A key player in the relationship between vascular dysfunction and inflammation in atherosclerosis. J Clin Med. (2020) 9:3026. doi: 10.3390/jcm9093026

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Baylis C. Arginine, arginine analogs and nitric oxide production in chronic kidney disease. Nat Clin Pract Nephrol. (2006) 2:209–20. doi: 10.1038/ncpneph0143

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Dimitroulas T, Hodson J, Sandoo A, Smith J, Kitas GD. Endothelial injury in rheumatoid arthritis: a crosstalk between dimethylarginines and systemic inflammation. Arthritis Res Ther. (2017) 19:32. doi: 10.1186/s13075-017-1232-1

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Chandrasekharan UM, Wang Z, Wu Y, Wilson Tang WH, Hazen SL, Wang S, et al. Elevated levels of plasma symmetric dimethylarginine and increased arginase activity as potential indicators of cardiovascular comorbidity in rheumatoid arthritis. Arthritis Res Ther. (2018) 20:123. doi: 10.1186/s13075-018-1616-x

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Kemeny-Beke A, Gesztelyi R, Bodnar N, Zsuga J, Kerekes G, Zsuga M, et al. Increased production of asymmetric dimethylarginine (ADMA) in ankylosing spondylitis: association with other clinical and laboratory parameters. Joint Bone Spine. (2011) 78:184–7. doi: 10.1016/j.jbspin.2010.05.009

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Surdacki A, Martens-Lobenhoffer J, Wloch A, Marewicz E, Rakowski T, Wieczorek-Surdacka E, et al. Elevated plasma asymmetric dimethyl-L-arginine levels are linked to endothelial progenitor cell depletion and carotid atherosclerosis in rheumatoid arthritis. Arthritis Rheum. (2007) 56:809–19. doi: 10.1002/art.22424

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Zafari P, Zarifian A, Alizadeh-Navaei R, Taghadosi M, Rafiei A, Samimi Z, et al. Asymmetric and symmetric dimethylarginine concentration as an indicator of cardiovascular diseases in rheumatoid arthritis patients: a systematic review and meta-analysis of case-control studies. Clin Rheumatol. (2020) 39:127–34. doi: 10.1007/s10067-019-04713-z

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Kerekes G, Soltesz P, Nurmohamed MT, Gonzalez-Gay MA, Turiel M, Vegh E, et al. Validated methods for assessment of subclinical atherosclerosis in rheumatology. Nat Rev Rheumatol. (2012) 8:224–34. doi: 10.1038/nrrheum.2012.16

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Li, H, Lu H, Tang W, Zuo J. Targeting methionine cycle as a potential therapeutic strategy for immune disorders. Expert Opin Ther Targets. (2017) 1–17. doi: 10.1080/14728222.2017.1370454

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Djuric D, Jakovljevic V, Zivkovic V, Srejovic I. Homocysteine and homocysteine-related compounds: an overview of the roles in the pathology of the cardiovascular and nervous systems. Can J Physiol Pharmacol. (2018) 96:991–1003. doi: 10.1139/cjpp-2018-0112

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Boger RH, Lentz SR, Bode-Boger SM, Knapp HR, Haynes WG. Elevation of asymmetrical dimethylarginine may mediate endothelial dysfunction during experimental hyperhomocyst(e)inaemia in humans. Clin Sci (Lond). (2001) 100:161–7. doi: 10.1042/cs1000161

CrossRef Full Text | Google Scholar

35. Doshi S, McDowell I, Goodfellow J, Stabler S, Boger R, Allen R, et al. Relationship between S-adenosylmethionine, S-adenosylhomocysteine, asymmetric dimethylarginine, and endothelial function in healthy human subjects during experimental hyper- and hypohomocysteinemia. Metabolism. (2005) 54:351–60. doi: 10.1016/j.metabol.2004.09.015

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Nightingale AK, James PP, Morris-Thurgood J, Harrold F, Tong R, Jackson SK, et al. Evidence against oxidative stress as mechanism of endothelial dysfunction in methionine loading model. Am J Physiol Heart Circ Physiol. (2001) 280:H1334–9. doi: 10.1152/ajpheart.2001.280.3.H1334

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Eleftheriadou I, Grigoropoulou P, Moyssakis I, Kokkinos A, Perrea D, Toutouzas K, et al. The effect of hyperhomocysteinemia on aortic distensibility in healthy individuals. Nutrition. (2013) 29:876–80. doi: 10.1016/j.nut.2012.12.026

PubMed Abstract | CrossRef Full Text | Google Scholar

38. van Dijk SC, Enneman AW, Swart KM, van Wijngaarden JP, Ham AC, Brouwer-Brolsma EM, et al. Effects of 2-year vitamin B12 and folic acid supplementation in hyperhomocysteinemic elderly on arterial stiffness and cardiovascular outcomes within the B-PROOF trial. J Hypertens. (2015) 33:1897–906;discussion1906. doi: 10.1097/HJH.0000000000000647

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Szekanecz Z, Kerekes G, Soltesz P. Vascular effects of biologic agents in RA and spondyloarthropathies. Nat Rev Rheumatol. (2009) 5:677–84. doi: 10.1038/nrrheum.2009.219

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Greenberg JD, Furer V, Farkouh ME. Cardiovascular safety of biologic therapies for the treatment of RA. Nat Rev Rheumatol. (2011) 8:13–21. doi: 10.1038/nrrheum.2011.168

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Smolen JS, Landewe RBM, Bijlsma JWJ, Burmester GR, Dougados M, Kerschbaumer A, et al. EULAR recommendations for the management of rheumatoid arthritis with synthetic and biological disease-modifying antirheumatic drugs: 2019 update. Ann Rheum Dis. (2020) 79:685–99.

Google Scholar

42. Xie W, Huang Y, Xiao S, Sun X, Fan Y, Zhang Z. Impact of Janus kinase inhibitors on risk of cardiovascular events in patients with rheumatoid arthritis: systematic review and meta-analysis of randomised controlled trials. Ann Rheum Dis. (2019) 78:1048–54. doi: 10.1136/annrheumdis-2018-214846

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Souto A, Salgado E, Maneiro JR, Mera A, Carmona L, Gomez-Reino JJ. Lipid profile changes in patients with chronic inflammatory arthritis treated with biologic agents and tofacitinib in randomized clinical trials: a systematic review and meta-analysis. Arthritis Rheumatol. (2015) 67:117–27. doi: 10.1002/art.38894

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Hamar A, Hascsi Z, Pusztai A, Czokolyova M, Vegh E, Petho Z, et al. Prospective, simultaneous assessment of joint and vascular inflammation by PET/CT in tofacitinib-treated patients with rheumatoid arthritis: associations with vascular and bone status. RMD Open. (2021) 7:e001804. doi: 10.1136/rmdopen-2021-001804

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Kume K, Amano K, Yamada S, Kanazawa T, Ohta H, Hatta K, et al. Tofacitinib improves atherosclerosis despite up-regulating serum cholesterol in patients with active rheumatoid arthritis: a cohort study. Rheumatol Int. (2017) 37:2079–85. doi: 10.1007/s00296-017-3844-9

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Fransen J, Welsing PMJ, De Keijzer RMH, van Riel PLCM. Development and validation of the DAS28 using CRP. Ann Rheum Dis. (2003) 62:10.

Google Scholar

47. Vegh E, Kerekes G, Pusztai A, Hamar A, Szamosi S, Vancsa A, et al. Effects of 1-year anti-TNF-alpha therapy on vascular function in rheumatoid arthritis and ankylosing spondylitis. Rheumatol Int. (2019) 40:427–36. doi: 10.1007/s00296-019-04497-0

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Klawitter J, Reed-Gitomer BY, McFann K, Pennington A, Klawitter J, Abebe KZ, et al. Endothelial dysfunction and oxidative stress in polycystic kidney disease. Am J Physiol Renal Physiol. (2014) 307:F1198–206. doi: 10.1152/ajprenal.00327.2014

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Keller AC, Klawitter J, Hildreth KL, Christians U, Putnam K, Kohrt WM, et al. Elevated plasma homocysteine and cysteine are associated with endothelial dysfunction across menopausal stages in healthy women. J Appl Physiol. (2019) 126:1533–40. doi: 10.1152/japplphysiol.00819.2018

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Faul F, Erdfelder E, Lang AG, Buchner A. G*Power 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav Res Methods. (2007) 39:175–91. doi: 10.3758/BF03193146

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Lorin J, Zeller M, Guilland JC, Cottin Y, Vergely C, Rochette L. Arginine and nitric oxide synthase: regulatory mechanisms and cardiovascular aspects. Mol Nutr Food Res. (2014) 58:101–16. doi: 10.1002/mnfr.201300033

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Morris SM Jr. Arginine metabolism: boundaries of our knowledge. J Nutr. (2007) 137:1602S–9S. doi: 10.1093/jn/137.6.1602S

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Sanderson SM, Gao X, Dai Z, Locasale JW. Methionine metabolism in health and cancer: a nexus of diet and precision medicine. Nat Rev Cancer. (2019) 19:625–37. doi: 10.1038/s41568-019-0187-8

PubMed Abstract | CrossRef Full Text | Google Scholar

54. van Hall G, Steensberg A, Fischer C, Keller C, Moller K, Moseley P, et al. Interleukin-6 markedly decreases skeletal muscle protein turnover and increases nonmuscle amino acid utilization in healthy individuals. J Clin Endocrinol Metab. (2008) 93:2851–8. doi: 10.1210/jc.2007-2223

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Stalder R, Zhang B, Jean Wrobel L, Boehncke WH, Brembilla NC. The Janus Kinase inhibitor tofacitinib impacts human dendritic cell differentiation and favours M1 macrophage development. Exp Dermatol. (2020) 29:71–8. doi: 10.1111/exd.14059

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Gremese E, Alivernini S, Tolusso B, Zeidler MP, Ferraccioli G. JAK inhibition by methotrexate (and csDMARDs) may explain clinical efficacy as monotherapy and combination therapy. J Leukoc Biol. (2019) 106:1063–8. doi: 10.1002/JLB.5RU0519-145R

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Karouzakis E, Gay RE, Gay S, Neidhart M. Increased recycling of polyamines is associated with global DNA hypomethylation in rheumatoid arthritis synovial fibroblasts. Arthritis Rheum. (2012) 64:1809–17. doi: 10.1002/art.34340

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Zhao G, He F, Wu C, Li P, Li N, Deng J, et al. Betaine in inflammation: Mechanistic aspects and applications. Front Immunol. (2018) 9:1070. doi: 10.3389/fimmu.2018.01070

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Bosello S, Santoliquido A, Zoli A, Di Campli C, Flore R, Tondi P, et al. TNF-alpha blockade induces a reversible but transient effect on endothelial dysfunction in patients with long-standing severe rheumatoid arthritis. Clin Rheumatol. (2008) 27:833–9. doi: 10.1007/s10067-007-0803-y

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Chatterjee Adhikari M, Guin A, Chakraborty S, Sinhamahapatra P, Ghosh A. Subclinical atherosclerosis and endothelial dysfunction in patients with early rheumatoid arthritis as evidenced by measurement of carotid intima-media thickness and flow-mediated vasodilatation: an observational study. Semin Arthritis Rheum. (2012) 41:669–75. doi: 10.1016/j.semarthrit.2011.08.003

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Protogerou AD, Zampeli E, Fragiadaki K, Stamatelopoulos K, Papamichael C, Sfikakis PP. A pilot study of endothelial dysfunction and aortic stiffness after interleukin-6 receptor inhibition in rheumatoid arthritis. Atherosclerosis. (2011) 219:734–6. doi: 10.1016/j.atherosclerosis.2011.09.015

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Maki-Petaja KM, Hall FC, Booth AD, Wallace SM, Yasmin PW, Bearcroft S, et al. Rheumatoid arthritis is associated with increased aortic pulse-wave velocity, which is reduced by anti-tumor necrosis factor-alpha therapy. Circulation. (2006) 114:1185–92. doi: 10.1161/CIRCULATIONAHA.105.601641

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Pieringer H, Stuby U, Pohanka E, Biesenbach G. Arterial stiffness in a muscular artery in women with longstanding rheumatoid arthritis compared with healthy controls and patients with traditional cardiovascular risk factors. Rheumatol Int. (2010) 30:1335–9. doi: 10.1007/s00296-009-1155-5

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Fan CY, Zhang ZY, Mei YF, Wu CJ, Shen BZ. Impaired brachial artery flow-mediated dilation and increased carotid intima-media thickness in rheumatoid arthritis patients. Chin Med J (Engl). (2012) 125:832–7.

Google Scholar

65. Alkaabi JK, Ho M, Levison R, Pullar T, Belch JJ. Rheumatoid arthritis and macrovascular disease. Rheumatology (Oxford). (2003) 42:292–7. doi: 10.1093/rheumatology/keg083

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Daza L, Aguirre M, Jimenez M, Herrera R, Bollain JJ. Common carotid intima-media thickness and von Willebrand factor serum levels in rheumatoid arthritis female patients without cardiovascular risk factors. Clin Rheumatol. (2007) 26:533–7. doi: 10.1007/s10067-006-0338-7

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Ytterberg SR, Bhatt DL, Mikuls TR, Koch GG, Fleischmann R, Rivas JL, et al. Cardiovascular and cancer risk with tofacitinib in rheumatoid arthritis. N Engl J Med. (2022) 386:316–26. doi: 10.1056/NEJMoa2109927

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Cohen SB, Tanaka Y, Mariette X, Curtis JR, Lee EB, Nash P, et al. Long-term safety of tofacitinib up to 9.5 years: a comprehensive integrated analysis of the rheumatoid arthritis clinical development programme. RMD Open. (2020) 6:e001395. doi: 10.1136/rmdopen-2020-001395

PubMed Abstract | CrossRef Full Text | Google Scholar

69. Sandoo A, Dimitroulas T, Veldhuijzen van Zanten JJ, Smith JP, Metsios GS, Nightingale P, et al. Lack of association between asymmetric dimethylarginine and in vivo microvascular and macrovascular endothelial function in patients with rheumatoid arthritis. Clin Exp Rheumatol. (2012) 30:388–96.

Google Scholar

70. Dimitroulas T, Sandoo A, Hodson J, Smith JP, Kitas GD. In vivo microvascular and macrovascular endothelial function is not associated with circulating dimethylarginines in patients with rheumatoid arthritis: a prospective analysis of the DRACCO cohort. Scand J Clin Lab Invest. (2016) 76:331–7. doi: 10.1080/00365513.2016.1177852

PubMed Abstract | CrossRef Full Text | Google Scholar

71. Tang WH, Wang Z, Cho L, Brennan DM, Hazen SL. Diminished global arginine bioavailability and increased arginine catabolism as metabolic profile of increased cardiovascular risk. J Am Coll Cardiol. (2009) 53:2061–7. doi: 10.1016/j.jacc.2009.02.036

PubMed Abstract | CrossRef Full Text | Google Scholar

72. Szekanecz Z, Kerekes G, Vegh E, Kardos Z, Barath Z, Tamasi L, et al. Autoimmune atherosclerosis in 3D: How it develops, how to diagnose and what to do. Autoimmun Rev. (2016) 15:756–69. doi: 10.1016/j.autrev.2016.03.014

PubMed Abstract | CrossRef Full Text | Google Scholar

73. Boisvert FM, Chenard CA, Richard S. Protein interfaces in signaling regulated by arginine methylation. Sci STKE. (2005) 2005:re2. doi: 10.1126/stke.2712005re2

PubMed Abstract | CrossRef Full Text | Google Scholar

74. Zsuga J, Torok J, Magyar MT, Valikovics A, Gesztelyi R, Lenkei A, et al. Dimethylarginines at the crossroad of insulin resistance and atherosclerosis. Metabolism. (2007) 56:394–9. doi: 10.1016/j.metabol.2006.10.023

PubMed Abstract | CrossRef Full Text | Google Scholar

75. Boger RH. Asymmetric dimethylarginine, an endogenous inhibitor of nitric oxide synthase, explains the “L-arginine paradox” and acts as a novel cardiovascular risk factor. J Nutr. (2004) 134(Suppl 10):2842S–7S. doi: 10.1093/jn/134.10.2842S

PubMed Abstract | CrossRef Full Text | Google Scholar

76. Wilson AM, Shin DS, Weatherby C, Harada RK, Ng MK, Nair N, et al. Asymmetric dimethylarginine correlates with measures of disease severity, major adverse cardiovascular events and all-cause mortality in patients with peripheral arterial disease. Vasc Med. (2010) 15:267–74. doi: 10.1177/1358863X10364552

PubMed Abstract | CrossRef Full Text | Google Scholar

77. Ye J, Dai Y, Mao H, Zheng W, Zhang J. Prognostic value of asymmetric dimethylarginine in patients with coronary artery disease: A meta-analysis. Nitric Oxide. (2021) 110:50–6. doi: 10.1016/j.niox.2021.03.002

PubMed Abstract | CrossRef Full Text | Google Scholar

78. Spinelli FR, Di Franco M, Metere A, Conti F, Iannuccelli C, Agati L, et al. Decrease of asymmetric dimethyl arginine after anti-TNF therapy in patients with rheumatoid arthritis. Drug Dev Res. (2014) 75(Suppl 1):S67–9. doi: 10.1002/ddr.21200

PubMed Abstract | CrossRef Full Text | Google Scholar

79. Korkosz M, Gasowski J, Surdacki A, Leszczynski P, Pawlak-Bus K, Jeka S, et al. Disparate effects of anti-TNF-alpha therapies on measures of disease activity and mediators of endothelial damage in ankylosing spondylitis. Pharmacol Rep. (2013) 65:891–7. doi: 10.1016/S1734-1140(13)71070-3

CrossRef Full Text | Google Scholar

80. Turiel M, Tomasoni L, Sitia S, Cicala S, Gianturco L, Ricci C, et al. Effects of long-term disease-modifying antirheumatic drugs on endothelial function in patients with early rheumatoid arthritis. Cardiovasc Ther. (2010) 28:e53–64. doi: 10.1111/j.1755-5922.2009.00119.x

PubMed Abstract | CrossRef Full Text | Google Scholar

81. Genre F, Lopez-Mejias R, Miranda-Filloy JA, Carnero-Lopez B, Gomez-Acebo I, Blanco R, et al. Asymmetric dimethylarginine serum levels in non-diabetic ankylosing spondylitis patients undergoing TNF-alpha antagonist therapy. Clin Exp Rheumatol. (2013) 31:749–55.

Google Scholar

82. Notsu Y, Yano S, Shibata H, Nagai A, Nabika T. Plasma arginine/ADMA ratio as a sensitive risk marker for atherosclerosis: Shimane CoHRE study. Atherosclerosis. (2015) 239:61–6. doi: 10.1016/j.atherosclerosis.2014.12.030

PubMed Abstract | CrossRef Full Text | Google Scholar

83. Mammedova JT, Sokolov AV, Freidlin IS, Starikova EA. The mechanisms of l-arginine metabolism disorder in endothelial cells. Biochemistry (Mosc). (2021) 86:146–55. doi: 10.1134/S0006297921020036

PubMed Abstract | CrossRef Full Text | Google Scholar

84. Angel K, Provan SA, Mowinckel P, Seljeflot I, Kvien TK, Atar D. The L-arginine/asymmetric dimethylarginine ratio is improved by anti-tumor necrosis factor-alpha therapy in inflammatory arthropathies. Associations with aortic stiffness. Atherosclerosis. (2012) 225:160–5. doi: 10.1016/j.atherosclerosis.2012.08.033

PubMed Abstract | CrossRef Full Text | Google Scholar

85. Roubenoff R, Dellaripa P, Nadeau MR, Abad LW, Muldoon BA, Selhub J, et al. Abnormal homocysteine metabolism in rheumatoid arthritis. Arthritis Rheum. (1997) 40:718–22. doi: 10.1002/art.1780400418

PubMed Abstract | CrossRef Full Text | Google Scholar

86. Alvarez-Cienfuegos A, Cantero-Nieto L, Garcia-Gomez JA, Sabio JM, Gonzalez-Gay MA, Ortego-Centeno N. Association between homocysteine levels and arterial stiffness in women with systemic sclerosis. Clin Exp Rheumatol. (2020) 125(38 Suppl):92–7.

Google Scholar

87. Wang XN, Ye P, Cao RH, Yang X, Xiao WK, Zhang Y, et al. Plasma homocysteine is a predictive factor for arterial stiffness: A community-based 4.8-year prospective study. J Clin Hypertens (Greenwich). (2015) 17:594–600. doi: 10.1111/jch.12555

PubMed Abstract | CrossRef Full Text | Google Scholar

88. Elias MF, Crichton GE, Abhayaratna WP. Interactions between plasma homocysteine and arterial stiffness in chronic kidney disease in community-dwelling individuals: The Maine-Syracuse Study. J Hum Hypertens. (2015) 29:726–31. doi: 10.1038/jhh.2015.17

PubMed Abstract | CrossRef Full Text | Google Scholar

89. van Waardenburg DA, de Betue CT, Luiking YC, Engel M, Deutz NE. Plasma arginine and citrulline concentrations in critically ill children: strong relation with inflammation. Am J Clin Nutr. (2007) 86:1438–44. doi: 10.1093/ajcn/86.5.1438

PubMed Abstract | CrossRef Full Text | Google Scholar

90. De Pietri S, Frandsen TL, Christensen M, Grell K, Rathe M, Muller K. Citrulline as a biomarker of bacteraemia during induction treatment for childhood acute lymphoblastic leukaemia. Pediatr Blood Cancer. (2020) 68:e28793. doi: 10.1002/pbc.28793

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: rheumatoid arthritis, tofacitinib, atherosclerosis, metabolomics, arginine, methionine

Citation: Soós B, Hamar A, Pusztai A, Czókolyová M, Végh E, Szamosi S, Pethő Z, Gulyás K, Kerekes G, Szántó S, Szűcs G, Christians U, Klawitter J, Seres T and Szekanecz Z (2022) Effects of tofacitinib therapy on arginine and methionine metabolites in association with vascular pathophysiology in rheumatoid arthritis: A metabolomic approach. Front. Med. 9:1011734. doi: 10.3389/fmed.2022.1011734

Received: 04 August 2022; Accepted: 26 October 2022;
Published: 10 November 2022.

Edited by:

Tue Wenzel Kragstrup, Aarhus University, Denmark

Reviewed by:

Patrick Dessein, University of the Witwatersrand, South Africa
Theodoros Dimitroulas, Aristotle University of Thessaloniki, Greece

Copyright © 2022 Soós, Hamar, Pusztai, Czókolyová, Végh, Szamosi, Pethő, Gulyás, Kerekes, Szántó, Szűcs, Christians, Klawitter, Seres and Szekanecz. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Zoltán Szekanecz, szekanecz.zoltan@med.unideb.hu

These authors have contributed equally to this work and share first authorship

These authors have contributed equally to this work and share last authorship

Download