Abstract
Systemic sclerosis (SSc) is a severe multi-organ disease with interstitial lung disease (ILD) being the major cause of death. While targeted therapies are emerging, biomarkers for sub-stratifying patients based on individual profiles are lacking. Herein, we investigated how levels of serum metabolites correlated with different stages of SSc and SSc-ILD. Serum samples of patients with SSc without ILD, stable and progressive SSc-ILD as well as of healthy controls (HC) were analysed using liquid targeted tandem mass spectrometry. The best discriminating profile consisted of 4 amino acids (AA) and 3 purine metabolites. l-tyrosine, l-tryptophan, and 1-methyl-adenosine distinguished HC from SSc patients. l-leucine, l-isoleucine, xanthosine, and adenosine monophosphate differentiated between progressing and stable SSc-ILD. In SSc-ILD, both, l-leucine and xanthosine negatively correlated with changes in FVC% predicted. Additionally, xanthosine was negatively correlated with changes in DLco% predicted and positively with the prognostic GAP index. Validation of l-leucine and l-isoleucine by an enzymatic assay confirmed both the sub-stratification of SSc-ILD patients and correlation with lung function and prognosis score. Serum metabolites may have potential as biomarkers for discriminating SSc patients based on the presence and severity of ILD. Confirmation in larger cohorts will be needed to appreciate their value for routine clinical care.
Similar content being viewed by others
Introduction
Systemic sclerosis (SSc) is a multi-system fibrotic disease with high morbidity and mortality1. Fibrosing interstitial lung disease (ILD) affects 60% of patients and is the leading cause of disease-related deaths2,3. Disease severity at time point of ILD diagnosis and progression over time determine the prognosis2,4,5. Currently available measures such as imaging, pulmonary function tests, or patient reported outcome scores have been shown to be of limited predictive value for the identification of SSc-ILD patients with poor prognosis6,7. Circulating biomarkers as readily available tools have been studied extensively with focus on serum proteins and auto-antibodies, but have not yet been established for daily use8.
A more recent research trend focuses on metabolic changes in diseased states9. Metabolic pathways represent central hubs within cells that have significant impact on complex signalling networks as well as cell functions. Metabolites are the final products of pathophysiologic processes caused by diseased states or by exposure to environmental factors or drugs and as such offer an integrated view of the biochemical milieu of bodily fluids or tissues10. The profiling of the metabolites present in biological samples—called metabolomics—has the potential to provide novel insights into pathophysiological processes11. Disturbed metabolism in the setting of (lung) fibrosis is increasingly acknowledged12. In general, the balance shifts from a catabolic to an anabolic state. In accordance, corresponding changes in the respective energy pathways occur, including e.g. an altered amino acid (AA) metabolism and enhanced protein production, an increased glycolysis with depletion of glucose, oxidative stress with mitochondrial dysfunction, decreased ATP levels, and reduced beta-oxidation of fatty acids13,14,15. The biologic relevance is supported by recent data showing that metabolic reprogramming successfully retransforms myofibroblasts to fibroblasts and activates degradation of extracellular matrix16,17,18, which offers novel therapeutic perspectives.
Metabolomics studies in ILD have so far mainly analysed lung tissue, bronchoalveolar lavage fluid (BALF) or exhaled breath condensate (EBC) from patients with idiopathic pulmonary fibrosis (IPF), the most prevalent form of ILD14,19. One serum-based study identified LysoPC, a potential precursor of lysophosphatic acid as discriminator between patients and healthy controls20. In SSc, data are scarce. Two recent studies showed some diagnostic potential since SSc patients could be distinguished from healthy controls based on their serum and/or urine metabolome with alterations in hexose, glycerolipid, and amino acid metabolism21,22.
Given the high unmet need of sub-stratifying biomarkers for ILD in SSc, the aim of our study was to evaluate the potential of serum metabolites to identify SSc patients according to their disease state with special focus on differentiating those with progressive from non-progressive ILD as base of therapeutic decision making.
Results
Patients’ characteristics
To decrease potential confounding factors for metabolomics studies, we matched the participating individuals (n = 12 per group, n = 48 in total) for age, sex, and time point of blood withdrawal. Demographic, clinical, and laboratory data are summarized in Table 1. The mean age was 61.5 ± 9.9 years, the majority of patients was female (n = 30, 83.3%), 25.0% (n = 9) had diffuse cutaneous and 27.8% (n = 10) had limited cutaneous disease, whereas 22.2% (n = 8) and 25.0% (n = 9) had only sclerodactyly or no skin involvement, respectively. The mean disease duration of SSc defined as onset of first non-Raynaud’s phenomenon was 9.0 ± 7.7 years. The mean disease duration of ILD was 3.2 ± 2.5 years. No differences in non-lung organ involvement were observed.
All SSc patients included in the study fulfilled the 2013 ACR/EULAR classification criteria23 at the time of blood withdrawal.
Patients with progressing ILD show different clinical features compared with patients with stable ILD
With respect to autoantibody profiles, a higher prevalence of anti-Centromere antibodies occurred in non-ILD patients (n = 9, 75%; p = 0.0045), of which the majority had only sclerodactyly or no skin involvement (n = 7, 58.3%) (Table 1). Anti-topoisomerase1-positivity, previously reported in the context of severe skin and lung involvement, was predominantly found in progressive SSc-ILD patients (n = 7, 58.3%, p = 0.0045, Fig. 1a) (Table 1, Fig. 1b).
While no difference was seen in time from onset from first non-Raynaud’s symptoms between the different patients’ groups, progressive patients had a significantly longer average duration of ILD than stable patients (mean ± SD = 4.7 ± 2.7 vs. 1.7 ± 1.0 years, p = 0.0028) and a more severe skin fibrosis as assessed by modified Rodnan skin score (mRSS) than both stable SSc-ILD and non-ILD patients (Table 1; Fig. 1c, mean ± SD = 10.3 ± 7.7 vs. 4.5 ± 5.8 and 3.7 ± 4.1, respectively; p = 0.024 and 0.011).
In addition, progressive SSc-ILD patients had a significantly lower forced vital capacity (FVC% predicted) than stable (mean ± SD = 80.7 ± 22.8 versus 108.1 ± 18.5 versus 102.1 ± 18.8, p = 0.0059 and p = 0.036, respectively) and non-ILD patients, while their diffusion capacity for carbon monoxide (DLco% predicted) was significantly lower than that of SSc patients without ILD (mean ± SD = 58.6 ± 25.1 versus 79.2 ± 13.1, p = 0.035) (Table 1, Fig. 1d,e). During the observation period, patients with progressive ILD experienced a significant decline in lung function compared with stable ILD patients (change FVC% predicted mean ± SD − 14.5 ± 7.1 versus 4.4 ± 7.1, p < 0.0001 and change DLco% predicted − 16.1 ± 21.6 versus 2.1 ± 14.6, p = 0.026; Table 1). Moreover, both the gender-age-physiology (GAP)24 and composite physiologic index (CPI)25 values, developed for the prediction of mortality in IPF patients, were higher in progressive than in stable SSc-ILD patients indicating a worse prognosis (mean ± SD = 2.6 ± 2.1 versus 1.2 ± 1.1 and 38.2 ± 20.8 versus 22.4 ± 11.5, p = 0.031 and 0.048, respectively, Table 1; Fig. 1f,g).
While there was a tendency towards a higher prevalence of immunosuppressive treatment, higher erythrocyte sedimentation rate (ESR) and higher C-reactive protein (CRP) levels in progressive SSc-ILD patients, changes were not statistically significant (Table 1). In addition, there was a non-significant tendency towards a lower body mass index (BMI) in progressive SSc-ILD compared to stable SSc-ILD and non-ILD patients (Table 1) (p = 0.35).
Metabolic serum profiling detects differences between disease subtypes
To identify differentially regulated serum metabolites as potential discriminators between healthy individuals and different SSc subgroups, targeted metabolic profiling for 110 metabolites (Supplementary Table S1) using targeted LC–MS/MS analysis was performed. After data processing and filtration, a total of 85 serum metabolites was detected, 56 in ESI (electron spray ionization) positive, 24 in ESI negative mode, and 5 in both modes (Supplementary Table S2). To test our hypothesis of distinct and discriminating metabolite patterns we performed multivariate analysis (hierarchical clustering and PLS-DA), followed by univariate analysis. Performance of significant metabolites from both analyses was then assessed by ROC analysis.
Hierarchical clustering of the targeted serum profiles suggested differences between the four groups (Fig. 2), which became more apparent in the subsequent PLS-DA (Fig. 3). Interestingly, the clearest separation was observed when comparing patients with progressing and stable ILD (Fig. 3b).
Serum metabolites show potential as biomarkers
To define the best discriminating metabolites we used variable importance in projection (VIP) scores (Supplementary Tables S3a,b) with a rather strict threshold set at ≥ 2. Based on this approach, we found distinct metabolite signatures for the different group comparisons (Table 2). Healthy individuals and SSc patients were best distinguished by changes in the amino acids (AAs) l-tyrosine and l-tryptophan, whereas SSc patients with and without ILD were best classified by dysregulation of l-threonine, xanthosine, 3-aminoisobutyric acid, and adenosine monophosphate. Progressors compared to stable ILD patients were characterized by alterations of l-leucine, l-isoleucine, xanthosine, and adenosine monophosphate.
In accordance with the results from the multivariate analysis, ANOVA-based univariate analysis with FDR correction for multiple testing identified l-leucine (p = 0.028), l-tyrosine (p = 0.077), xanthosine (p = 0.032), l-tryptophan (p = 0.028 for ESI+ and ESI− modes), and 1-methyladenosine (p = 0.077) as significantly altered between groups.
Analysis of peak areas showed that levels of l-leucine and l-isoleucine were highest in healthy individuals and gradually decreased from SSc patients without ILD to those with stable ILD. Interestingly, patients with progressing ILD had increased levels compared with stable SSc-ILD patients (Fig. 4a,b). For xanthosine, we observed significantly lower levels in healthy individuals compared with non-ILD SSc patients as well as for stable SSc-ILD patients compared with progressing SSc-ILD patients (Fig. 4c).
l-Tryptophan, l-tyrosine, l-threonine, and adenosine monophosphate showed a gradual decline from healthy individuals, SSc without ILD, stable ILD to progressive ILD (Fig. 4d–h), whereas the opposite occurred for 3-aminoisobutyric acid and 1-methyl-adenosine (Fig. 4i,j).
Amino acid profiles discriminate between groups and correlate with changes in lung function and prognostic scores
To select the best discriminators between two groups, we performed ROC curve analysis for the previously identified candidate serum metabolites.
l-Tryptophan, detected in both ESI+ and ESI− modes, best differentiated between healthy individuals and SSc patients (AUC 0.859, 95% CI 0.745–0.973; AUC 0.884, 95% CI 0.788–0.981), followed by 1-methyl-adenosine (AUC 0.822, 95% CI 0.705–0.939) and l-tyrosine (AUC 0.812, 95% CI 0.667–0.958) (Fig. 5a).
In the group of ILD patients, l-leucine and l-isoleucine performed best to distinguish progressing from stable SSc-ILD patients (AUC 0.847, 95% CI 0.695–1.00; AUC 0.826, 95% CI 0.656–0.997), followed by adenosine monophosphate (AUC 0.785, 95% CI 0.598–0.971), and xanthosine (AUC 0.771, 95% CI 0.551–0.990) (Fig. 5b). No significant performance of metabolites in separation of non-ILD SSc and SSc-ILD patients was observed.
Notably, in patients with SSc-ILD, both, l-leucine and xanthosine negatively correlated with changes in the lung function parameter FVC% predicted (r = − 0.48 and − 0.51; p = 0.016 and 0.011, respectively), while xanthosine also negatively correlated with changes in DLco% predicted (r = − 0.58; p = 0.0038) as well as with the GAP index (r = 0.67, p = 0.0005) (Fig. 6). The correlations with changes in absolute values of FVC and DLco% predicted are shown in Supplementary Fig. S1a,b.
Validation of higher values of the BCAAs l-leucine and l-isoleucine in progressive SSc-ILD compared to stable disease with the use of an enzymatic assay resulted in similar results as LC–MS/MS, with significantly higher values detected in progressive patients (mean = 286.5 and 235.5 µM, for progressive and stable patients, respectively; p = 0.005) (Fig. 7a). In ROC analysis (AUC 0.818, 95% CI 0.631–1.032), a cut-off value of 250.3 µM separated stable from progressive patients with a sensitivity of 72.7% and a specificity of 83.3% (Fig. 7b). Furthermore, BCAA levels negatively correlated with changes in FVC (r = − 0.55; p = 0.0063) and DLco% predicted (r = − 0.56; p = 0.0064) and positively with the GAP index (r = 0.66, p = 0.0005) (Fig. 7c–e). The correlations with changes in absolute values of FVC and DLco% predicted are shown in Supplementary Fig. S1c.
For external validation, we assessed BCAA levels in an independent, prospectively followed cohort of SSc-ILD patients from Paris. Herein, we found a similar trend towards higher BCAA levels in patients with future progression of ILD as compared with patients, who remained stable during follow-up (mean = 311.8 ± 26.5 and 289.6 ± 49.4 µM, respectively; p = 0.26). In ROC analysis, BCAA levels correlated negatively with DLco% predicted (r = − 0.38, p = 0.027), and positively with the mortality-predicting GAP index (r = 0.39, p = 0.022). The results are shown in Supplementary Fig. S2, for details including patients’ characteristics refer to data supplement.
In order to assess whether elevated BCAA levels specifically reflected fibrotic and/or pulmonary processes, we measured BCAA serum levels in 29 patients with non-fibrotic, primary myositis26 without concomitant ILD. Interestingly, BCAA levels significantly correlated with disease activity (Spearman’s rho = 0.52, p = 0.0037). Patients with CK levels exceeding the defined upper limit of normal (170 U/l) had significantly higher BCAA levels than patients with low CK levels (mean = 309.5 ± 88.2 and 253.4 ± 41.6 U/l for active and inactive disease, respectively; p = 0.032). No correlation with systemic inflammation, measured by CRP as a surrogate marker, was observed (Pearson’s r = − 0.15, p = 0.46). All results are visualized in Supplementary Fig. S3, for details including patients’ characteristics refer to data supplement.
Discussion
Our study assessed the potential of serum metabolites as circulating biomarkers for disease stage and severity of SSc(-ILD).
Serum metabolite profiling yielded a final set of 4 amino acids and 3 purine metabolites. Changes in the levels of l-tyrosine, l-tryptophan, and 1-methyl-adenosine distinguished HC from SSc and alterations in l-leucine, l-isoleucine, xanthosine, and adenosine monophosphate profiles differentiated between progressing and stable SSc-ILD with l-tryptophan and l-leucine being the best performing discriminators in the respective groups. Increased serum levels of l-leucine and l-isoleucine in progressing compared with stable SSc-ILD patients were confirmed with an independent enzymatic assay with definition of a critical threshold.
Our results are in accordance with the previously reported changes in energy metabolism in fibrotic conditions13,14,15. Similar to our findings, previous studies in SSc and/or fibrosing ILD reported decreased serum levels of l-tryptophan27,28 or upregulation in IPF lungs14. Serotonin29, a downstream product of l-tryptophan, was upregulated in fibrotic conditions30 including SSc31. Inhibition of the serotonin receptors had anti-fibrotic effects in experimental conditions and in a first in human proof-of-concept study31,32,33,34. Elevated levels of 1-methyl-adenosine29 have so far been associated with proliferative and/or metastatic tumors35,36. Non-modified adenosine, generated extracellularly from ATP/ADP breakdown29, is a well-known pro-fibrotic mediator37,38,39,40,41 and its inhibition was beneficial in in fibrotic animal models37,39,40. As previously observed in lung fibrosis and SSc14, there was a tendency towards lower adenosine monophosphate (AMP) levels in progressing SSc-ILD patients. Interestingly, AMP-activated protein kinase (AMPK), a critical sensor of energy sufficiency, acts as central metabolic switch in cell metabolism and thereby opposes mTOR (mechanistic target of rapamycin) signaling42. Activation of AMPK causes a shift from anabolism to catabolism to generate ATP to restore energy homeostasis. In fibrotic conditions, AMPK activity was decreased43. Activation of AMPK reversed established lung fibrosis in the bleomycin-induced lung fibrosis model16,44. Patients with progressing compared with stable ILD displayed increased levels of l-leucine and l-isoleucine. Both essential AA are highly abundant in elastin45, a major extracellular matrix (ECM) protein, and could thus reflect increased ECM degradation 46. Both AAs were upregulated in the lung tissue and the exhaled breath condensate of patients with IPF14,19. BCAAs, particularly l-leucine, stimulate protein synthesis and reduce protein breakdown via the phosphorylation of mTOR47. mTOR plays an important role in anabolic processes by causing cells to switch from oxidative phosphorylation to aerobic glycolysis48. In SSc und pulmonary fibrosis, mTOR activity is increased49. Its inhibition by rapamycin prevented experimental fibrosis and showed some benefit in diffuse cutaneous SSc patients in a small, randomized, phase 1 study48. Finally, levels of xanthosine, a purine nucleoside, were higher in progressing compared with stable SSc-ILD. Purine receptors were suggested to play a role in fibrosis50.
Clinically, in our study, progressive SSc-ILD patients were characterized by a higher prevalence of anti-topoisomerase 1 antibodies, more severe skin fibrosis, worse lung function and worse prognostic scores. In addition, they showed a substantial decline of lung function in the observation period. Of note, both l-leucine and xanthosine negatively correlated with pulmonary function (changes in FVC% and DLco% predicted respectively) and one ILD-mortality prediction score (GAP index), which underlines their relevance as biomarker candidates. Most importantly, we could validate l-leucine and l-isoleucine, the best discriminators of progressing vs. stable SSc-ILD, in an independent experiment. In the BCAA assay, a defined cut-off value of 250.3 µM distinguished stable from progressive patients with good sensitivity and specificity and again BCAA levels negatively correlated with changes of lung function parameters and positively with the GAP index. Similar results were obtained by analysing an independent, external validation cohort of SSc-ILD patients from Paris.
The limitations of our study mainly arise from the limited number of patients. Validation in external multi-centre cohorts will be needed to assess the future usefulness in clinical routine. Furthermore, potential correlations with other protein biomarkers should be assessed. Prediction modelling for progression of ILD with circulating biomarkers, clinical, functional, and imaging parameters would be ideal to test the performance of circulating biomarkers only models versus mixed models. This, however, again warrants large datasets and multi-centre cohorts.
Disturbances in AA metabolism have been reported in other studies on SSc and IPF. Differences in identified metabolite profiles might arise from different analysis and detection methodologies (i.e. mass spectrometry or ion exchange chromatography) or differences in sample collection, storage and processing. We decided on a large-scale targeted analysis on a triple quadrupole mass spectrometer rather than an untargeted full-scan approach applying high-resolution mass spectrometry due to the increased sensitivity, linearity, reproducibility and straight-forward metabolite identification of targeted LC–MS/MS acquisition51. Although the applied assay covered multiple differentially regulated pathways, it is limited to the tested 110 metabolites. Further research could be carried out using untargeted metabolic profiling or a targeted assay covering a wider range of metabolites such as carbohydrates or phospholipids.
Additionally, given the real-life scenario, we cannot exclude that the fasting state and/or the diet might have some influence especially on the measured AA levels. Supplementation with micronutrients and/or vitamins did not occur in our patients’ cohort and we draw blood at approximately the same time in the morning to eliminate potentially confounding factors as well as possible. The fact, however, that changes in (BC)AA levels were reported consistently in other studies of lung fibrosis argues against a strong or exclusively dietary influence.
Furthermore, we have to take into account that SSc is a multi-organ disease. In our study, patients with progressive ILD also had more extensive skin disease pointing towards a more severe disease state. We can therefore not assume that changes in BCAA serum levels exclusively reflect lung pathology. Taking into account the correlation between CK and BCAA levels in primary myositis patients, it seems likely that in both myositis and SSc(-ILD), high serum BCAA levels reflect overall disease activity characterized by a switch to an anabolic state with subsequent changes in AA metabolism. Thus, in these complex diseases, changes in BCAA levels can probably not be attributed to either pro-fibrotic or immune processes since they may rather reflect the disturbed metabolism that arises from global tissue remodelling with varying contributions of different cell types. This argues, however, not against the usefulness of BCAA as progression markers in a given disease context.
In conclusion, our study suggests that serum metabolites might have potential as circulating biomarkers for discriminating stable and progressive SSc patients. Confirmation in larger multi-cohorts will be needed to fully appreciate their value for routine clinical care.
Methods
Patients and controls
For this study, SSc patients from the University Hospital Zurich’s prospective SSc patients cohort were divided in the following three groups: patients without ILD (non-ILD), patients with stable ILD and patients with progressive ILD (n = 12 per group).
Progressive ILD was defined as either a relative decrease in FVC% predicted of ≥ 10% independent of changes in DLco% predicted, a decrease in FVC% predicted of 5–9% combined with a decrease in DLco% predicted of ≥ 15%, or an increase of pathologic lung involvement in high resolution computed tomography (HRCT) from < 20% to > 20% compared to the previous visit [mean follow-up interval = 14 months (range = 9–26)]52,53. Stable ILD was defined as the absence of the above-mentioned criteria for progression in any of the visits recorded in the EUSTAR database54. Non-ILD patients were defined as having no evidence of lung involvement on HRCT scans.
Progressive SSc-ILD patients were matched with stable SSc-ILD and non-ILD SSc patients as well as healthy controls (HC, n = 12) for age, sex, and time point of blood withdrawal (morning).
In addition, serum BCAA levels were analyzed in two additional cohorts of SSc-ILD and primary myositis patients. Detailed information on these patients can be found in the Supplementary Methods, data supplement pp. 15.
Serum collection and processing was performed following a standardized protocol in accordance with international guidelines55. Aliquots of serum were stored at − 80 °C until further processing.
Written informed consent was obtained from all enrolled individuals.
The study was approved by the canton of Zurich’s ethics committee (approval numbers: pre-BASEC-EK-839 (KEK-no.-2016-01515), KEK-ZH-no. 2010-158/5, BASEC-no. 2018-01873, BASEC-no. 2018-02165, BASEC-no. 2017-01298, KEK-ZH-Nr. 2012-0419) and all experiments were conducted in accordance with Swiss legislation and regulation.
Ultra-high performance liquid chromatography coupled to tandem mass spectrometry (UHPLC-MS/MS)
Serum profiling of 110 metabolites of SSc patients and HC was performed using a targeted ultra-high performance liquid chromatography coupled to tandem mass spectrometry (LC–MS/MS) assay as described previously56.
For metabolite extraction, frozen serum samples were thawed at room temperature and 300 μL of ice-cold acetonitrile:ethanol (1:1) were added to 100 μL of each sample. Samples as well as an equally treated water-only control were vortexed, incubated at − 80 °C for 30 min, centrifuged at 14,000g and 4 °C for 15 min in order to pellet precipitated proteins. 350 μL of supernatant of each sample were transferred to fresh tubes and dried using a speed vacuum centrifuge at room temperature and minimum pressure 5.1 Torr (Savant SPD1010, Thermo). Dried samples were stored at − 80 °C until further processing. Before data acquisition, samples were reconstituted in 300 μL 10% methanol, resulting in a final dilution of 1:3 in respect to the initial serum volume, sonicated for 1 min in a water bath in order to ensure complete reconstitution, and centrifuged at 14,000g and 4 °C for 20 min. Supernatants were transferred into glass vials and a quality control (QC) sample was generated by pooling 10 µL of all extracted samples. After, samples were stored at 10 °C until analysis.
Extracted samples were analysed on a Xevo TQ-S triple quadrupole mass spectrometer interfaced with an electrospray ionization source and coupled to an ACQUITY UPLC I-Class system (both Waters, USA). Chromatographic separation of 1 µL extract was performed using reversed-phase chromatography (ACQUITY UPLC HSS T3, 2.1 × 100 mm, 1.8 µm column, Waters, USA) with mobile phases composed of (A) 0.1% formic acid in H2O and (B) 0.1% formic acid in methanol. Further LC conditions and MS-specific parameters were previously described56. LC–MS/MS analysis was performed in negative electrospray ionization (ESI−) and positive electrospray ionization (ESI+) mode by two independent injections.
Before analysis of the first sample, the instrument was equilibrated by injecting the QC sample 10 times. Samples were then injected in a randomized block design order with intermittent analysis of the QC sample after every fifth sample in order to observe instrumental fluctuations.
Branched chain amino acid assay
For quantification of the branched chain amino acids (BCAAs) leucine, isoleucine and valine, a commercial colorimetric analysis kit (ab83374, Abcam, United Kingdom) was used following the manufacturer’s instructions. Briefly, leucine standards and 2.5 times diluted serum samples were incubated for 30 min at room temperature with an equal volume of enzyme and substrate mix, initiating the colorimetric reaction. Absorption at 450 nm was measured using a GloMax-Multi Detection System microplate reader (Promega, USA) and sample BCAA concentration was calculated in relation to the standard curve.
Data analysis
Raw metabolomics data were processed in Skyline 4.2 (MacCoss Lab Software, USA). Metabolites were selected for statistical analysis if the peak area decreased linearly in diluted samples, no detection of background noise in extracted blank sample was observed, the coefficient of variation was ≤ 20% in QC samples, as well as, a proper peak shape was detected. Peak areas were then used for further statistical analysis as specified below.
Statistical data analysis and graphical visualization was performed using R 3.6 with the mixOmics package and GraphPad Prism 8.0.0 (GraphPad Software, USA).
For univariate analysis and hierarchical clustering, metabolite peak areas were Z-score transformed for normalization.
Univariate analysis of metabolomics data was performed by applying one-way ANOVA as well as Tukey’s post-hoc test for multi-group data. Two-group data was analyzed by Student’s t-test or Mann–Whitney U test for parametric and nonparametric data, respectively. Categorical data were analysed using a Chi-Square test. Performance of potential biomarkers was assessed by Receiver Operating Characteristic (ROC) curve analysis. Pearson’s (parametric) and Spearman’s (non-parametric) correlation was used to assess linear relationships between metabolites and clinical parameters. For multivariate analysis, data were subjected to partial least-square discriminant analysis (PLS-DA) with variable importance in projection (VIP) scores of ≥ 2 being considered statistically significant.
Data are presented either as medians with interquartile range (boxplots; horizontal line = median, boxes = interquartile range) or as means with standard deviation (SD; tables). For false discovery rate (FDR)-corrected univariate analysis of metabolites (excluding post-hoc testing), the significance level was set to 0.1, while for all other analyses p-values < 0.05 were considered statistically significant.
Data availability
All data are presented either in the main text or in the data supplement. The datasets generated and/or analysed during the current study are available from the corresponding author on reasonable request.
References
Elhai, M. et al. Mapping and predicting mortality from systemic sclerosis. Ann. Rheum. Dis. 76, 1897–1905. https://doi.org/10.1136/annrheumdis-2017-211448 (2017).
Hoffmann-Vold, A.-M. et al. Tracking impact of interstitial lung disease in systemic sclerosis in a complete nationwide cohort. Am. J. Respir. Crit. Care Med. https://doi.org/10.1164/rccm.201903-0486OC (2019).
Steen, V. D. & Medsger, T. A. Changes in causes of death in systemic sclerosis, 1972–2002. Ann. Rheum. Dis. 66, 940–944. https://doi.org/10.1136/ard.2006.066068 (2007).
Goh, N. S. et al. Interstitial lung disease in systemic sclerosis: A simple staging system. Am. J. Respir. Crit. Care Med. 177, 1248–1254. https://doi.org/10.1164/rccm.200706-877OC (2008).
Goh, N. S. et al. Short-term pulmonary function trends are predictive of mortality in interstitial lung disease associated with systemic sclerosis. Arthritis Rheumatol.. 69, 1670–1678. https://doi.org/10.1002/art.40130 (2017).
Wu, W. et al. Prediction of progression of interstitial lung disease in patients with systemic sclerosis: The SPAR model. Ann. Rheum. Dis. 77, 1326–1332. https://doi.org/10.1136/annrheumdis-2018-213201 (2018).
Salaffi, F. et al. Computer-aided tomographic analysis of interstitial lung disease (ILD) in patients with systemic sclerosis (SSc). Correlation with pulmonary physiologic tests and patient-centred measures of perceived dyspnea and functional disability. PLoS ONE 11, e0149240. https://doi.org/10.1371/journal.pone.0149240 (2016).
Hoffmann-Vold, A.-M., Fretheim, H., Meier, C. & Maurer, B. Circulating biomarkers of systemic sclerosis—interstitial lung disease. J. Scleroderma Relat. Disord. https://doi.org/10.1177/2397198319894851 (2020).
Eming, S. A., Wynn, T. A. & Martin, P. Inflammation and metabolism in tissue repair and regeneration. Science 356, 1026–1030. https://doi.org/10.1126/science.aam7928 (2017).
Athersuch, T. J. & Keun, H. C. Metabolic profiling in human exposome studies. Mutagenesis 30, 755–762. https://doi.org/10.1093/mutage/gev060 (2015).
Rinschen, M. M., Ivanisevic, J., Giera, M. & Siuzdak, G. Identification of bioactive metabolites using activity metabolomics. Nat. Rev. Mol. Cell Biol. 20, 353–367. https://doi.org/10.1038/s41580-019-0108-4 (2019).
Reisdorph, N. A., Cruickshank-Quinn, C., Nkrumah-Elie, Y. & Reisdorph, R. Application of metabolomics in lung research. Methods Mol. Biol. (Clifton, N.J.) 263–288, 2018. https://doi.org/10.1007/978-1-4939-8570-8_17 (1809).
Zhao, Y. D. et al. Metabolic heterogeneity of idiopathic pulmonary fibrosis: A metabolomic study. BMJ Open Respir. Res. 4, e000183. https://doi.org/10.1136/bmjresp-2017-000183 (2017).
Kang, Y. P. et al. Metabolic profiling regarding pathogenesis of idiopathic pulmonary fibrosis. J. Proteome Res. 15, 1717–1724. https://doi.org/10.1021/acs.jproteome.6b00156 (2016).
Rabinowitz, J. D. & Mutlu, G. M. A metabolic strategy to reverse fibrosis?. Nat. Metab. 1, 12–13. https://doi.org/10.1038/s42255-018-0013-8 (2019).
Rangarajan, S. et al. Metformin reverses established lung fibrosis in a bleomycin model. Nat. Med. 24, 1121–1127. https://doi.org/10.1038/s41591-018-0087-6 (2018).
Zhao, X. et al. Metabolic regulation of dermal fibroblasts contributes to skin extracellular matrix homeostasis and fibrosis. Nat. Metab. 1, 147–157. https://doi.org/10.1038/s42255-018-0008-5 (2019).
Xie, N. et al. Glycolytic reprogramming in myofibroblast differentiation and lung fibrosis. Am. J. Respir. Crit. Care Med. 192, 1462–1474. https://doi.org/10.1164/rccm.201504-0780OC (2015).
Gaugg, M. T. et al. Molecular breath analysis supports altered amino acid metabolism in idiopathic pulmonary fibrosis. Respirology 24, 437–444. https://doi.org/10.1111/resp.13465 (2019).
Rindlisbacher, B., Schmid, C., Geiser, T., Bovet, C. & Funke-Chambour, M. Serum metabolic profiling identified a distinct metabolic signature in patients with idiopathic pulmonary fibrosis—A potential biomarker role for LysoPC. Respir. Res. 19, 7. https://doi.org/10.1186/s12931-018-0714-2 (2018).
Murgia, F. et al. Metabolomic profile of systemic sclerosis patients. Sci. Rep. 8, 7626. https://doi.org/10.1038/s41598-018-25992-7 (2018).
Fernandez-Ochoa, A. et al. Urinary and plasma metabolite differences detected by HPLC-ESI-QTOF-MS in systemic sclerosis patients. J. Pharm. Biomed. Anal. 162, 82–90. https://doi.org/10.1016/j.jpba.2018.09.021 (2019).
van den Hoogen, F. et al. 2013 classification criteria for systemic sclerosis: An American College of Rheumatology/European League against Rheumatism collaborative initiative. Arthritis Rheum. 65, 2737–2747. https://doi.org/10.1002/art.38098 (2013).
Ryerson, C. J. et al. Predicting survival across chronic interstitial lung disease: The ILD-GAP model. Chest 145, 723–728. https://doi.org/10.1378/chest.13-1474 (2014).
Wells, A. U. et al. Idiopathic pulmonary fibrosis: A composite physiologic index derived from disease extent observed by computed tomography. Am. J. Respir. Crit. Care Med. 167, 962–969. https://doi.org/10.1164/rccm.2111053 (2003).
Lundberg, I. E. et al. European League Against Rheumatism/American College of Rheumatology classification criteria for adult and juvenile idiopathic inflammatory myopathies and their major subgroups. J. Ann. Rheum. Dis. 76(1955–1964), 2017. https://doi.org/10.1136/annrheumdis-2017-211468 (2017).
Csipo, I. et al. Decreased serum tryptophan and elevated neopterin levels in systemic sclerosis. Clin. Exp. Rheumatol. 13, 269–270 (1995).
Meyer, K. C. et al. Tryptophan metabolism in chronic inflammatory lung disease. J Lab Clin Med 126, 530–540 (1995).
Wishart, D. S. et al. HMDB 4.0: The human metabolome database for 2018. Nucleic Acids Res. 46, D608–D617. https://doi.org/10.1093/nar/gkx1089 (2017).
Dolivo, D. M., Larson, S. A. & Dominko, T. Tryptophan metabolites kynurenine and serotonin regulate fibroblast activation and fibrosis. Cell Mol. Life Sci. 75, 3663–3681. https://doi.org/10.1007/s00018-018-2880-2 (2018).
Dees, C. et al. Platelet-derived serotonin links vascular disease and tissue fibrosis. J. Exp. Med. 208, 961–972. https://doi.org/10.1084/jem.20101629 (2019).
Beretta, L. et al. A polymorphism in the human serotonin 5-HT2A receptor gene may protect against systemic sclerosis by reducing platelet aggregation. Arthritis Res. Ther. 10, R103. https://doi.org/10.1186/ar2495 (2008).
Chaturvedi, S. et al. 5-HT2 and 5-HT2B antagonists attenuate pro-fibrotic phenotype in human adult dermal fibroblasts by blocking TGF-beta1 induced non-canonical signaling pathways including STAT3: Implications for fibrotic diseases like scleroderma. Int. J. Rheum. Dis. 21, 2128–2138. https://doi.org/10.1111/1756-185X.13386 (2018).
Distler O. M. B., Vettori, S., Blumhardt, S., Frey, D., Distler, A., Beyer, C., Distler, J. H. The serotonin receptor 2 inhibitor terguride has beneficial effects on skin fibrosis: Results from a phase 2 proof of concept study. Ann. Rheum. Dis. 66 (2016).
Vreken, P. & Tavenier, P. Urinary excretion of six modified nucleosides by patients with breast carcinoma. Ann. Clin. Biochem. 24(Pt 6), 598–603. https://doi.org/10.1177/000456328702400608 (1987).
Langridge, J. I. et al. Gas chromatography/mass spectrometric analysis of urinary nucleosides in cancer patients; potential of modified nucleosides as tumour markers. Rapid Commun. Mass Spectrom. 7, 427–434. https://doi.org/10.1002/rcm.1290070605 (1993).
Karmouty-Quintana, H. et al. The antifibrotic effect of A2B adenosine receptor antagonism in a mouse model of dermal fibrosis. Arthritis Rheumatol. 70, 1673–1684. https://doi.org/10.1002/art.40554 (2018).
Della Latta, V., Cabiati, M., Rocchiccioli, S., Del Ry, S. & Morales, M. A. The role of the adenosinergic system in lung fibrosis. Pharmacol. Res. 76, 182–189. https://doi.org/10.1016/j.phrs.2013.08.004 (2013).
Chan, E. S. et al. Adenosine A2A receptors in diffuse dermal fibrosis: Pathogenic role in human dermal fibroblasts and in a murine model of scleroderma. Arthritis Rheum. 54, 2632–2642. https://doi.org/10.1002/art.21974 (2006).
Sun, C. X. et al. Role of A2B adenosine receptor signaling in adenosine-dependent pulmonary inflammation and injury. J. Clin. Investig. 116, 2173–2182. https://doi.org/10.1172/JCI27303 (2006).
Luo, F. et al. Extracellular adenosine levels are associated with the progression and exacerbation of pulmonary fibrosis. FASEB J. 30, 874–883. https://doi.org/10.1096/fj.15-274845 (2016).
Mukhopadhyay, S. et al. Reciprocal regulation of AMP-activated protein kinase and phospholipase D. J. Biol. Chem. 290, 6986–6993. https://doi.org/10.1074/jbc.M114.622571 (2015).
Wang, J. et al. Exacerbated cardiac fibrosis induced by beta-adrenergic activation in old mice due to decreased AMPK activity. Clin. Exp. Pharmacol. Physiol. 43, 1029–1037. https://doi.org/10.1111/1440-1681.12622 (2016).
Lu, J. et al. Activation of AMPK by metformin inhibits TGF-beta-induced collagen production in mouse renal fibroblasts. Life Sci. 127, 59–65. https://doi.org/10.1016/j.lfs.2015.01.042 (2015).
Debelle, L. & Alix, A. J. The structures of elastins and their function. Biochimie 81, 981–994. https://doi.org/10.1016/s0300-9084(99)00221-7 (1999).
Jenkins, R. G. et al. Longitudinal change in collagen degradation biomarkers in idiopathic pulmonary fibrosis: An analysis from the prospective, multicentre PROFILE study. Lancet Respir. Med. 3, 462–472. https://doi.org/10.1016/S2213-2600(15)00048-X (2015).
Wilkinson, D. J. et al. Effects of leucine and its metabolite beta-hydroxy-beta-methylbutyrate on human skeletal muscle protein metabolism. J. Physiol. 591, 2911–2923. https://doi.org/10.1113/jphysiol.2013.253203 (2013).
Suto, T. & Karonitsch, T. The immunobiology of mTOR in autoimmunity. J. Autoimmun. https://doi.org/10.1016/j.jaut.2019.102373 (2019).
Lawrence, J. & Nho, R. The role of the mammalian target of rapamycin (mTOR) in pulmonary fibrosis. Int. J. Mol. Sci. https://doi.org/10.3390/ijms19030778 (2018).
Distler, O. & Cozzio, A. Systemic sclerosis and localized scleroderma-current concepts and novel targets for therapy. Semin. Immunopathol. 38, 87–95. https://doi.org/10.1007/s00281-015-0551-z (2016).
Chen, S. et al. Pseudotargeted metabolomics method and its application in serum biomarker discovery for hepatocellular carcinoma based on ultra high-performance liquid chromatography/triple quadrupole mass spectrometry. Anal. Chem. 85, 8326–8333. https://doi.org/10.1021/ac4016787 (2013).
Distler, O. et al. Nintedanib for systemic sclerosis-associated interstitial lung disease. N. Engl. J. Med. 380, 2518–2528. https://doi.org/10.1056/NEJMoa1903076 (2019).
Tashkin, D. P. et al. Mycophenolate mofetil versus oral cyclophosphamide in scleroderma-related interstitial lung disease (SLS II): A randomised controlled, double-blind, parallel group trial. Lancet Respir. Med. 4, 708–719. https://doi.org/10.1016/S2213-2600(16)30152-7 (2016).
Meier, F. M. et al. Update on the profile of the EUSTAR cohort: An analysis of the EULAR Scleroderma Trials and Research group database. Ann. Rheum. Dis. 71, 1355–1360. https://doi.org/10.1136/annrheumdis-2011-200742 (2012).
Beyer, C. et al. EUSTAR biobanking: Recommendations for the collection, storage and distribution of biospecimens in scleroderma research. Ann. Rheum. Dis. 70, 1178–1182. https://doi.org/10.1136/ard.2010.142489 (2011).
Freiburghaus, K. et al. Metabolomics by UHPLC–MS: Benefits provided by complementary use of Q-TOF and QQQ for pathway profiling. Metabolomics 15, 120. https://doi.org/10.1007/s11306-019-1585-3 (2019).
Funding
This work was supported by the Prof. Max Cloetta Foundation (to BM) and by the Forschungskredit UZH Postdoc (FK-19-046 to JS).
Author information
Authors and Affiliations
Contributions
B.M., C.M., and C.B. were centrally involved in the conception and the planning of the study. Experiments were performed by C.M., K.F., and C.B. C.M., B.M., C.N., K.F., C.B., J.S., Y.A. and O.D. participated in the data analysis and interpretation. C.M., B.M., and C.N. prepared the figures. B.M. and C.M. wrote the manuscript. All authors revised and edited the manuscript.
Corresponding author
Ethics declarations
Competing interests
CM, KF, CB, JS, CN: none to declare. OD had speaker fees from Actelion, Bayer, Boehringer Ingelheim, Medscape, Novartis, Roche, Menarini, Mepha, MSD, iQone, Pfizer, consultancy fees from Abbvie, Actelion, Acceleron Pharma, Amgen, AnaMar, Bayer, Baecon Discovery, Blade Therapeutics, Boehringer, CSL Behring, ChemomAb, Curzion Pharmaceuticals, Ergonex, Galapagos NV, GSK, Glenmark Pharmaceuticals, Inventiva, Italfarmaco, iQvia, medac, Medscape, Mitsubishi Tanabe Pharma, MSD, Roche, Sanofi, UCB, Amgen, Lilly, Target BioScience, Pfizer, project scoring fees from Abbvie, interview fees from Catenion, travel support from Pfizer, research grants from Actelion, Bayer, Boehringer Ingelheim, Competitive Drug Development International Ltd, Mitsubishi Tanabe, patent issued: mir-29 for the treatment of systemic sclerosis (US8247389, EP2331143). BM had grant/research support from AbbVie, Protagen, Novartis, and congress support from Pfizer, Roche, Actelion, and MSD, speakers honorary from Novartis. In addition, B. Maurer has a patent mir-29 for the treatment of systemic sclerosis registered (US8247389, EP2331143). The real or perceived potential conflicts listed above are accurately stated. YA had consultancy relationships and/or has received research funding from Actelion, Bayer, Boehringer Ingelheim, ChemomAb, Genentech/Roche, Inventiva, Pfizer, Sanofi, Servier, in the area of potential treatments of scleroderma and its complications. AH received consulting fees and/or research funding from Actelion, ARXX, Bayer, Boehringer Ingelheim, Medscape, MSD and Roche.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Meier, C., Freiburghaus, K., Bovet, C. et al. Serum metabolites as biomarkers in systemic sclerosis-associated interstitial lung disease. Sci Rep 10, 21912 (2020). https://doi.org/10.1038/s41598-020-78951-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-020-78951-6
This article is cited by
-
The interplay between rheumatic diseases and pulmonary health
Rheumatology International (2024)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.