Article Text

Biomarkers of haemodynamic severity of systemic sclerosis-associated pulmonary arterial hypertension by serum proteome analysis
Free
  1. Sébastien Sanges1,2,3,4,5,6,
  2. Lisa Rice1,
  3. Ly Tu7,8,
  4. Eleanor Valenzi9,
  5. Jean-Luc Cracowski10,
  6. David Montani7,8,11,
  7. Julio C Mantero1,
  8. Camille Ternynck12,
  9. Guillemette Marot12,13,14,
  10. Andreea M Bujor1,
  11. Eric Hachulla2,3,4,5,6,
  12. David Launay2,3,4,5,6,
  13. Marc Humbert7,8,11,
  14. Christophe Guignabert7,8,
  15. Robert Lafyatis15
  1. 1 Boston University School of Medicine, E5 Arthritis Center, Boston, Massachusetts, USA
  2. 2 Univ. Lille, U1286 – INFINITE – Institute for Translational Research in Inflammation, Lille, France
  3. 3 INSERM, Lille, France
  4. 4 CHU Lille, Département de Médecine Interne et Immunologie Clinique, Lille, France
  5. 5 Centre National de Référence Maladies Auto-immunes Systémiques Rares du Nord et Nord-Ouest de France (CeRAINO), Lille, France
  6. 6 Health Care Provider of the European Reference Network on Rare Connective Tissue and Musculoskeletal Diseases Network (ReCONNET), Lille, France
  7. 7 Université Paris Saclay, School of Medicine, Le Kremlin-Bicêtre, France
  8. 8 INSERM UMR_S 999, Hôpital Marie-Lannelongue, Le Plessis-Robinson, France
  9. 9 Division of Pulmonary, Allergy and Critical Care Medicine, Department of Medicine, University of Pittsburgh Medical Center, Pittsburgh, Pennsylvania, USA
  10. 10 Univ. Grenoble Alpes, INSERM, HP2, Grenoble, France
  11. 11 AP-HP, Department of Respiratory and Intensive Care Medicine, Hôpital Bicêtre, Kremlin-Bicêtre, France
  12. 12 Univ. Lille, CHU Lille, ULR 2694 - METRICS: Évaluation des technologies de santé et des pratiques médicales, Lille, France
  13. 13 Inria, MODAL: MOdels for Data Analysis and Learning, Lille, France
  14. 14 Univ. Lille, CNRS, Inserm, CHU Lille, Institut Pasteur de Lille, UAR 2014 - US 41 - PLBS, bilille, Lille, France
  15. 15 Division of Rheumatology, Department of Medicine, University of Pittsburgh Medical Center, Pittsburgh, Pennsylvania, USA
  1. Correspondence to Dr Sébastien Sanges, Département de Médecine Interne et Immunologie Clinique, CHU Lille, F-59037 Lille Cedex, France; sebastien.sanges{at}univ-lille.fr

Abstract

Objectives To mine the serum proteome of patients with systemic sclerosis-associated pulmonary arterial hypertension (SSc-PAH) and to detect biomarkers that may assist in earlier and more effective diagnosis and treatment.

Methods Patients with limited cutaneous SSc, no extensive interstitial lung disease and no PAH-specific therapy were included. They were classified as cases if they had PAH confirmed by right heart catheterisation (RHC) and serum collected on the same day as RHC; and as controls if they had no clinical evidence of PAH.

Results Patients were mostly middle-aged females with anticentromere-associated SSc. Among 1129 proteins assessed by a high-throughput proteomic assay (SOMAscan), only 2 were differentially expressed and correlated significantly with pulmonary vascular resistance (PVR) in SSc-PAH patients (n=15): chemerin (ρ=0.62, p=0.01) and SET (ρ=0.62, p=0.01). To validate these results, serum levels of chemerin were measured by ELISA in an independent cohort. Chemerin levels were confirmed to be significantly higher (p=0.01) and correlate with PVR (ρ=0.42, p=0.04) in SSc-PAH patients (n=24). Chemerin mRNA expression was detected in fibroblasts, pulmonary artery smooth muscle cells (PA-SMCs)/pericytes and mesothelial cells in SSc-PAH lungs by single-cell RNA-sequencing. Confocal immunofluorescence revealed increased expression of a chemerin receptor, CMKLR1, on SSc-PAH PA-SMCs. SSc-PAH serum seemed to induce higher PA-SMC proliferation than serum from SSc patients without PAH. This difference appeared neutralised when adding the CMKLR1 inhibitor α-NETA.

Conclusion Chemerin seems an interesting surrogate biomarker for PVR in SSc-PAH. Increased chemerin serum levels and CMKLR1 expression by PA-SMCs may contribute to SSc-PAH pathogenesis by inducing PA-SMC proliferation.

  • systemic sclerosis
  • autoimmune diseases
  • cytokines

Data availability statement

Data are available on reasonable request.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

WHAT IS ALREADY KNOWN ON THIS TOPIC

  • Invasive assessment of pulmonary vascular resistance (PVR) by right heart catheterisation is essential in the therapeutic and risk stratification strategies for patients with systemic sclerosis-associated pulmonary arterial hypertension (SSc-PAH).

WHAT THIS STUDY ADDS

  • Using a wide-scale proteomic approach, we identified chemerin as a potential non-invasive surrogate for PVR in SSc-PAH patients.

  • Chemerin seems to induce pulmonary arterial smooth muscle cell proliferation and contribute to the progression of SSc-PAH through binding with its cognate receptor chemokine-like receptor 1 (CMKLR1).

HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY

  • Chemerin could be used as a non-invasive assessment of haemodynamic severity in SSc-PAH patients.

  • Further studies could evaluate the ability of chemerin to assess treatment response during follow-up, as well as the efficacy of therapeutic strategies targeting the chemerin-CMKLR1 axis.

Introduction

Pulmonary arterial hypertension (PAH) is the most devastating complication of systemic sclerosis (SSc).1 A common pathological hallmark of PAH is remodelling of precapillary pulmonary arteries (PAs), manifested by pulmonary vasoconstriction and medial thickening due to increased expansion of PA smooth muscle cells (SMCs).2 This alteration of the small/medium-sized PAs contributes to the progressive and rapid increase in pulmonary vascular resistance (PVR) that eventually lead to right ventricular failure.3 Right heart catheterisation (RHC) is required to confirm the diagnosis of PAH, to assess the severity of the haemodynamic impairment and the response to treatment.4

PAH has become a leading cause of death in SSc, with a standardised mortality ratio of 5.27 and a median survival of 3 years.1 5 Improving its prognosis is thus a major challenge when managing SSc patients,3 requiring early diagnosis to allow for timely initiation of treatment, as well as rapid detection of treatment failure to allow for the immediate adjustment of medications. However, since diagnosis and disease progression in PAH have haemodynamic definitions, this supposes performing iterative RHCs, highlighting the need to identify biomarkers for the accurate and non-invasive prediction of PAH in SSc patients.

Several works have previously tried to address this issue.6–12 Most of them focused on identifying surrogates for haemodynamic diagnostic parameters, with biomarkers that could accurately discriminate SSc patients with and without PAH.7–10 Few studies, however, have tried to find surrogate markers for severity parameters, such as PVR.6 11 This is a matter of importance, since PVR reflects the ongoing vascular remodelling underlying the disease progression and can be used to guide treatment initiation and assess therapeutic efficacy.11 13

This work was an exploratory study that aimed to identify proteins that correlate with haemodynamic severity in SSc-PAH patients and to serve as a base on which future hypotheses could be tested. To achieve this, we used a wide-scale approach by investigating alterations in the patients’ serum proteome with a high-throughput assay. In a second step, we assessed whether the identified candidate biomarkers were involved in disease pathogenesis.

Methods

Serum proteome signature of SSc-PAH

Study population

In a first step, to identify candidate biomarkers, a discovery cohort was recruited from Boston University Arthritis Centre. Patients were included if they met the following criteria: (1) diagnosis of SSc according to the 2013 American College of Rheumatology (ACR)/European Alliance of Associations for Rheumatology (EULAR) criteria;14 (2) limited cutaneous (lc) subset of SSc according to LeRoy’s criteria;15 (3) absence of extensive interstitial lung disease (ILD) according to Goh’s criteria; and16 (4) absence of PAH-specific therapy. In that first step, patients with diffuse cutaneous (dc)SSc and/or extensive ILD were excluded to avoid our biomarker screening to be interfered with by any active organ involvement other than PAH, and our results to be mediated by any pathogenic process other than pulmonary microangiopathy. In addition, extensive ILD was also considered an exclusion criterion to avoid including patients with secondary (group 3) pulmonary hypertension that probably did not have significant pulmonary microangiopathy. Patients fulfilling these inclusion criteria were classified as cases if they had an RHC-proven diagnosis of PAH according to the 2015 European Society of Cardiology (ESC)/European Respiratory Society (ERS) guidelines4; and as controls if they had no evidence of PAH. PAH screening modalities, assessment of PAH probability and referral for RHC followed ESC/ERS guidelines in effect at the time of sample collection. For cases, all serum samples were collected on the same day as RHC.

In a second step, in order to confirm the validity of the identified candidate biomarkers, an independent validation cohort was recruited from Boston University, Grenoble University Hospital and Lille national reference centre for SSc. Patients were included if they met the same criteria as the discovery cohort. Cases and controls were also similarly defined. For cases, serum samples were collected on the same day as RHC (except for four samples collected within 6 weeks before).

In a third step, in order to determine if our candidate biomarkers were influenced by other SSc manifestations apart from PAH, three other patient groups were recruited from Boston University Arthritis Centre. The first group included patients with dcSSc, no extensive ILD and no evidence of PAH. A second group included patients with lcSSc, extensive ILD and no evidence of PAH. A third group included healthy controls.

All patients were recruited consecutively by screening the databases of participating centres in 2015. Serum samples had been collected between 2004 and 2014.

SOMAscan assay

The relative expression levels of 1129 serum proteins were assessed in the discovery cohort using the SOMAscan platform (SomaLogic, Boulder, Colorado, USA), a highly multiplexed aptamer microarray, as previously described.17 18 Briefly, samples were deposited in a 96-well plate and incubated with a mixture of the 1129 SOMAmer reagents. Two sequential bead-based immobilisation and washing steps eliminated unbound proteins, non-specifically bound proteins and unbound SOMAmer reagents. The remaining specifically bound SOMAmer reagents were isolated, and each reagent was quantified simultaneously on a custom Agilent hybridisation array. The amount of each SOMAmer measured was quantitatively proportional to the protein concentration in the original sample. Results were expressed in relative fluorescence units.

ELISA assays

To confirm result reproducibility and cross-platform consistency, serum concentrations of candidate proteins were measured in duplicate in the validation cohort using commercial ELISA assays (Quantikine ELISA Human Chemerin Immunoassay, RD systems; SET Nuclear Oncogene (SET) ELISA Kit (Human), cat. #ABIN1152447, antibodies-online.com) at appropriate dilutions (1:100 for chemerin; 1:1 for SET), according to the manufacturer’s protocol.

Statistical analyses

For the description of study populations, quantitative variables were expressed as means±SD or medians (IQR) for non-normal distributions, and categorical variables were expressed as numbers (percentage). Normality of distributions was assessed using histograms and tested using the Shapiro-Wilk test.

In a first step, data from the discovery cohort were log2 transformed and a differential analysis was performed between cases and controls with the Bioconductor R (V.3.6.1) package limma (V.3.38.3).19 Limma uses an empirical Bayesian approach to estimate variances in moderated t-tests, which has proven to improve results on standard t-tests, especially when the number of replicates is small. Raw p values were adjusted with the Benjamini Hochberg method20 and proteins with adjusted p values <0.05 were considered differentially expressed. For cases, the Spearman correlations between PVR and expressions of these differentially expressed proteins were tested and were considered significant if raw p values were <0.05.

In a second step, differential expression of candidate biomarkers identified in the discovery cohort was assessed between cases and controls from the validation cohort using Mann-Whitney test; and correlations with PVR in cases were analysed using Spearman test.

In a third step, serum levels of the candidate biomarkers were compared between the five patient groups using the Kruskal-Wallis test. Pairwise comparison was conducted by post-hoc Dunn test followed by Bonferroni correction.

As protein production could always be quantified in all datasets, no imputation for missing data was performed. A heatmap was drawn for differentially expressed proteins with the R package pheatmap (V.1.0.12) after standardisation of the expression data. Other figures were created using GraphPad Prism V.9.1.0 (GraphPad Software, California, USA).

Chemerin/CMKLR1 axis and pulmonary vascular remodelling

Study population

Patients were recruited from Le Kremlin-Bicêtre and Pittsburgh reference centres for PAH. Lung specimens were obtained during transplantation in patients diagnosed with SSc-PAH and/or ILD or idiopathic PAH (iPAH),4 and during lobectomy or pneumonectomy for localised lung cancer in control subjects (non-SSc controls). Lung specimens from controls were collected at a distance from the tumour foci. Preoperative echocardiography was performed in controls to rule out PAH.

Single-cell RNA sequencing in explanted human lungs

Data generation

Lung samples from four SSc-PAH patients and four non-SSc controls were transported in Perfadex and processed within 30 min of explantation as previously described.21 Of note, three SSc patients had extensive ILD. Following mechanical and enzymatic digestion by Liberase DL and DNase I, resulting single-cell suspensions were loaded into 10 x Genomics Chromium instrument (Pleasanton, California, USA), with scRNA-seq library preparation per manufacturer’s protocol. Sequencing was performed on an Illumina NovaSeq 6000 instrument through the UPMC Genome Centre.

Data analysis

Data analysis was performed with R (V.4.0.4) package Seurat (V.4.1.0). Samples were merged with batch correction by Harmony for individual sample, followed by normalisation, identification and visualisation of cell populations, and differential expression testing. Cell populations were identified by gene markers and visualised by uniform manifold and approximation projection plot.

Confocal immunofluorescence analyses in explanted human lungs

Immunofluorescent staining for chemerin and chemokine-like receptor 1 (CMKLR1) was performed on lung paraffin sections from SSc-PAH patients, SSc-no PAH patients and non-SSc controls as previously described.22 Of note, both SSc patient groups had extensive ILD. Briefly, lung sections were deparaffined and incubated with retrieval buffer. Then, sections were saturated with blocking buffer and incubated overnight with specific antibodies (sc-398769, CliniSciences, France), followed by corresponding secondary fluorescent-labelled antibodies (Thermo-Fisher Scientific, France). Nuclei were labelled using 4',6-diamidino-2-phénylindole (DAPI) (Thermo-Fisher Scientific). Mounting was done using ProLong Gold antifade reagent (Thermo-Fisher Scientific). Images were taken using LSM700 confocal microscope (Zeiss, France).

Role of CMKLR1 in the proliferation of human PA-SMCs

PA-SMC proliferation experiments

PA-SMCs were isolated from distal PA of lung explants from iPAH patients and non-SSc controls, and cultured as previously described.23 24 The isolated PA-SMCs were strongly positive for α-smooth muscle actin (α-SMA), smooth muscle-specific SM22 protein and calponin, and negative for von Willebrand factor and CD31. Cells were used at passage <5.

The mRNA expression of CMKLR1 was measured by real-time quantitative PCR using TaqMan gene expression assay (assay ID: Hs01081979_s1) as previously described.22 25 Relative quantification was calculated by normalising the Ct (threshold cycle) of the gene of interest to the Ct of 18S in the same sample, according to the comparative CT method (ΔΔCT method).

PA-SMCs were cultured with 5% serum from 5 SSc-PAH and 5 SSc-no PAH patients included in the proteomic validation cohort, in the presence or absence of 2-(alpha-naphthoyl)ethyltrimethylammonium iodide (α-NETA) at 1 µM. Proliferation was assessed by 5-bromo-2-deoxyuridine (BrdU) incorporation and direct cell counting.

Descriptive analysis

Figures were created using GraphPad Prism V.9.1.0 (GraphPad Software, California, USA). Differences were assessed visually without statistical analysis due to small sample size.

Results

Serum levels of chemerin are elevated and correlate with PVR in SSc-PAH patients

In a first step, serum expression of 1129 biomarkers was assessed by SOMAscan in 15 cases (SSc-PAH patients) and 16 controls (SSc-no PAH patients) from a discovery cohort (table 1; for additional details, see also online supplemental table 1). We identified 53 proteins differentially expressed between the 2 groups (figure 1). Among them, two analytes showed a significant correlation with PVR in cases: chemerin (ρ=0.62, p=0.01), an adipokine (figure 2); and SET nuclear protooncogene (ρ=0.62, p=0.01). Similar results were found when removing an outlier patient with markedly high PVR value from the analysis.

Supplemental material

Figure 1

Heatmap depicting the differential expression of the 53 candidate biomarkers identified in the discovery cohort: standardised expression values of the 53 differentially expressed proteins in the discovery cohort. Red values indicate overexpression and blue values underexpression compared with the mean expression level.

Figure 2

Differential expression of chemerin and correlation with PVR in the discovery and the validation cohorts. PAH, pulmonary arterial hypertension; PVR, pulmonary vascular resistance; RFU, relative fluorescence units; SSc, systemic sclerosis.

Table 1

Characteristics of SSc-PAH patients from the discovery and the validation cohorts

In a second step, to confirm these results, serum levels of chemerin and SET were then measured by ELISA in 24 cases and 17 controls from an independent validation cohort (table 1; online supplemental table 1). Consistently, serum levels of chemerin were significantly higher (p=0.01) and correlated significantly with PVR levels (ρ=0.42, p=0.04) in cases (figure 2). Serum SET levels were undetectable in both cases and controls when measured by ELISA, consistent with its intracellular nature.

In a third step, in order to determine if other SSc manifestations influenced biomarker concentrations, chemerin serum levels from our discovery cases and controls were compared with other patient groups included in the SOMAscan dataset: patients with dcSSc, no extensive ILD and no PAH (n=11); patients with lcSSc, extensive ILD and no PAH (n=10); and healthy controls (n=11). Circulating concentrations of chemerin were only increased in the PAH group, and similar to healthy controls in all others (figure 3).

Figure 3

Differential expression of chemerin in different SSc patient groups. DC, diffuse cutaneous; ext ILD, extensive interstitial lung disease; HC, healthy controls; lc, limited cutaneous; PAH, pulmonary arterial hypertension; RFU, relative fluorescence units; SSc, systemic sclerosis.

Overall, these results could suggest that chemerin seems a potential surrogate biomarker for PVR in SSc-PAH. Since results with SET were not reproduced, and since literature data suggested chemerin to be a more promising pathophysiological lead, we chose to focus our investigations on this latter protein.

Expression of chemerin receptor CMKLR1 is increased on PA-SMCs from SSc-PAH patients

Because the binding of chemerin to its receptor CMKLR1 stimulates the proliferation and migration of PA-SMCs,26 we first performed a single cell RNA-sequencing (scRNAseq) analysis in lungs from four SSc-PAH and four non-SSc controls to obtain a global view of chemerin and CMKLR1 expression patterns. Our scRNAseq data indicate that CMKLR1 was predominantly expressed by a subpopulation of cells expressing α-SMA and clustering with PA-SMCs/pericytes and myofibroblasts; as well as by endothelial cells and macrophages (figure 4). However, no significant difference between SSc-PAH patients and non-SSc controls was observed in CMKLR1 mRNA expression levels in these different cell populations. In addition, our scRNAseq data indicated that chemerin is mostly expressed by fibroblasts, PA-SMCs/pericytes and mesothelial cells, with a significant increase in SSc-PAH patients as compared with non-SSc controls for this last cell population (figure 4).

Figure 4

Expression of CMKLR1, chemerin and α-smooth muscle actin mRNA in cell populations from SSc-PAH lungs assessed by single-cell RNA-sequencing. Top panels provide information from all cell populations; Bottom panels provide information on mesenchymal cell populations. α-SMA, α-smooth muscle actin; CMKLR1, chemokine-Like Receptor 1; PAH, pulmonary arterial hypertension; SSc, systemic sclerosis.

To confirm our scRNAseq data, confocal microscopic analyses were next performed on lung specimens dually labelled with CMLKR1 and a specific PA-SMC marker, α-SMA. Examination of CMLKR1 protein expression patterns showed CMLKR1 staining to be more pronounced in the smooth muscle layer in both patients with SSc and SSc-PAH groups relative to non-SSc controls (figure 5). No difference was observed in chemerin staining between groups (online supplemental figure S1).

Supplemental material

Figure 5

Representative images of lung sections immunostained with DAPI (blue), α-SMA (green) and CMKLR1 (red) from non-SSc controls (top row), SSc-no PAH (middle row) and SSc-PAH (bottom row). Scale bar=50 µm in all sections. α-SMA, a-smooth muscle actin; DAPI, 4',6-diamidino-2-phénylindole; PAH, pulmonary arterial hypertension; SSc, systemic sclerosis.

Taken altogether, these data could suggest an upregulation in the expression pattern of the chemerin/CMRLK1 axis in SSc-PAH pulmonary vessels.

The SSc-PAH serum-induced PA-SMC proliferation seems inhibited by a chemerin-CMKLR1 inhibitor

To further examine the functional consequences of these changes on PA-SMC proliferation, we next determined the efficacy of the CMKLR1 antagonist α-NETA to attenuate the proliferation of PA-SMCs derived from iPAH patients. First, we confirmed an increased CMKLR1 mRNA production in PA-SMCs from iPAH vs control (online supplemental figure S2).

Supplemental material

Then, PA-SMCs from iPAH patients were stimulated with serum obtained from five SSc-PAH and five SSc-no PAH patients in the presence or absence of α-NETA (figure 6). Serum from SSc-PAH cases seemed to induce higher PA-SMC proliferation than serum from SSc-no PAH controls. This difference seemed neutralised in the presence of α-NETA.

Figure 6

Proliferation of PA-SMCs from idiopathic PAH patients, after stimulation with serum from SSc-PAH and SSc-no PAH patients, in the presence or absence of α-NETA. α-NETA: 2-(anaphthoyl)ethyltrimethylammonium iodide; BrdU, 5-bromo-2-deoxyuridine; i, idiopathic; PAH, pulmonary arterial hypertension; PA-SMC, pulmonary arterial smooth muscle cells; SSc, systemic sclerosis.

Overall, these findings suggest that serum from SSc-PAH patients seems to induce proliferation of PA-SMCs and that inhibition of CMKLR1 activation could partly abolish this phenomenon.

Discussion

To our knowledge, this is the first study that focused on identifying a surrogate marker for haemodynamic severity in SSc-PAH using a wide-scale approach. Our results can be summarised as follows: (1) chemerin seems to be a potential surrogate marker for haemodynamic severity in SSc-PAH, as it showed robust correlations with PVR; (2) in lungs, chemerin mRNA was detected in fibroblasts, PA-SMCs/pericytes and mesothelial cell populations; (3) elevated chemerin serum levels and increased expression of its receptor CMKLR1 by PA-SMCs could contribute to SSc-PAH pathogenesis by inducing PA-SMC proliferation.

Chemerin as a surrogate marker for haemodynamic severity in SSc-PAH

Since PAH is characterised by progressive pulmonary vascular remodelling leading to increased PVR, PVR is thought to reflect the severity of this process.13 As such, it can be used to monitor treatment efficacy and failure. Developing non-invasive methods to assess PVR is thus a major unmet need in the field of SSc-PAH management, that has only been scarcely investigated so far.6 11 Although there is some dispersion on the scatterplots that could limit future works, the correlation between PVR and chemerin serum levels was reproducible and can be described as moderate to strong based on Spearman coefficient values. In a recent work,11 Bauer et al measured 313 proteins on sera from the DETECT cohort and identified a panel of 8 analytes capable of accurately discriminating SSc-PAH from SSc-no PAH patients. When testing the ability of these biomarkers to predict PVR, they found moderate correlations for a model comprised of 5 proteins: RAGE, NT-pro-BNP, IGFBP-7, SP-D and VCAM-1. Since this study focused primarily on identifying diagnostic biomarkers, and not PVR surrogates, further comparison with our results is challenging.

Interestingly, the increase in chemerin concentrations appeared restricted to SSc-PAH patients in our study. Previous publications have reported on chemerin serum levels in SSc, although with conflicting results.27–29 Akamata et al found similar values between SSc patients and healthy controls, and no difference between cutaneous subsets of the disease.27 Chemerin levels correlated with modified Rodnan skin score (mRSS) and disease duration in dcSSc patients, as well as with digital ulcers in the whole cohort; but not with DLCO nor elevated systolic pulmonary arterial pressure (sPAP) on echocardiography.27 According to Chighizola et al, chemerin levels were lower in dcSSc patients but similar in lcSSc patients when compared with healthy controls; and correlated with disease duration (especially in dcSSc) but not with capillaroscopic pattern, ILD or digital ulcers.28 Sawicka et al observed higher chemerin values in SSc patients, especially among dcSSc, than in controls.29 They were also associated with acute phase reactants in the whole cohort and with mRSS in dcSSc patients, but not with disease duration, forced vital capacity nor DLCO.29 In summary, while no association was identified with ILD, previous studies have reported chemerin levels in dcSSc patients as similar, lower and higher that lcSSc and healthy controls. This discrepancy is challenging to explain. It should be noted however that our work used a different design, as we deliberately constituted homogeneous patient subgroups and compared chemerin levels between them, rather than including a broad heterogeneous patient population and assessed correlations. We believe that this allowed us to better ascertain the disease parameter responsible for elevated chemerin levels. Moreover, it is possible that the associations observed in dcSSc patients were in fact mediated by undiagnosed PAH (as chemerin levels were reported to correlate with sPAP among dcSSc patients29). Finally, this may also be explained by measurement of different chemerin isoforms.30 In any case, further studies including well-defined homogeneous patient groups are needed to determine whether abnormal serum levels of chemerin are specific to PAH in SSc.

Of note, serum chemerin levels were also studied in in patients with iPAH and were significantly increased compared with controls, with a satisfactory sensitivity (85.7%) and specificity (100%) at a concentration of 471.76 pg/mL.31

Chemerin as a contributor to the pathogenesis of SSc-PAH

Chemerin is pleiotropic protein that exerts multiple effects (acting to varying degrees as a chemokine, an adipokine and/or a growth factor) on different tissues (notably the immune system, the adipose tissue and the vasculature).30 32 In blood vessels, chemerin can act on endothelial cells (regulating angiogenesis, cell adhesion and nitric oxide production)33–36 as well as SMCs (inducing contraction, migration and proliferation).37–40 As such, it has been implicated in several cardiovascular diseases, such as arterial hypertension, atherosclerosis, diabetic microangiopathy and pre-eclampsia.30 32 Chemerin also yields proinflammatory effects by attracting and stimulating cytokine production from macrophages, dendritic cells and natural killer cells41–43; and thus has also been involved in various immunologic conditions, such as rheumatoid arthritis and psoriasis.30

In our work, chemerin mRNA was expressed primarily by fibroblasts and smooth muscle/pericytes, though not increased compared with non-SSc controls; and its protein was not detected by immunostaining. This discrepancy may reflect difficulties in detecting this secreted protein, limitations of the staining antibody or lack of mRNA translation, possibly suggesting extrapulmonary production. In our previous work detailing scRNAseq data in iPAH patients, chemerin was also expressed primarily by fibroblasts and upregulated 2.075-fold.44 These data suggest that upregulated secretion of chemerin by adventitial pulmonary fibroblasts and/or peripheral tissues may contribute to the pathogenesis of SSc-PAH.

Interestingly, previous studies have suggested a role of chemerin in both SSc and PAH pathogenesis.26 27 45 In SSc patients, chemerin expression is increased in dermal endothelial cells; and circulating levels of chemerin correlated with occurrence of digital ulcers.27 Chemerin expression was also decreased in dermal fibroblasts, due to autocrine TGF-β secretion and Fli1 deficiency; and serum chemerin levels correlated with mRSS.27 Of note, CMKLR1 expression in SSc patients’ skin was similar to healthy controls’.27 In PA from healthy rats, chemerin and CMKLR1 expression is detected in both endothelial cells and SMCs.46 Exposure to chemerin potentiated vascular responses to vasoconstrictors (phenylephrine, endothelin-1 and serotonin) in rat PA and impaired acetylcholine-induced PA vasodilatation, by mechanisms involving at least partly NO signalling and oxidative stress.46 In primary cultured PASMCs from healthy rats, exposure to chemerin induced a dose-dependent proliferation and migration (potentiated by combination with endothelin-1),26 31 and reduced staurosporine-induced apoptosis.26 Chemerin treatment did not alter gene expression of IL-6, IL-6R and IL-1β, but increased TNF-α at high doses.26 Healthy rat PASMCs exposed to chemerin up-regulated the ERK1/2 pathways; and treatment by an ERK inhibitor annulled the chemerin-induced proliferation.31 In the monocrotaline (MCT) rat model of PAH, chemerin-induced contraction of intrapulmonary arteries is increased compared with controls.45 Lungs from MCT rats showed increased expression of CMKLR1 and production of various chemerin isoforms (possibly due to glycosylation) compared with healthy animals.45 Expression of CMKLR1 in MCT rats was increased in SMCs but decreased expression in endothelial cells compared with controls.31 45 Finally, in studies performed in other vascular diseases, chemerin was able to induce angiogenesis in vitro47; and knockdown of chemerin significantly inhibited SMC proliferation and neointimal hyperplasia in a rat model of vascular injury induced by balloon angioplasty.48

Overall, these data suggest that the chemerin-CMKLR1 axis could contribute to the pathogenesis of SSc-PAH through vascular effects but could also be involved in the fibrotic and inflammatory events occurring during the disease.

Strengths and limitations

Our study draws strength from the sample collection on the same day as RHC (which allowed to perform accurate haemodynamic correlations), its wide-scale proteomic approach, the validation of our results on an independent cohort, and evidences for the pathophysiological relevance of the candidate biomarkers.

It also has limitations. First, we limited our study population to untreated patients with lcSSc and non-extensive ILD, so that our biomarker screening would not be interfered with by an active cutaneous or interstitial lung involvement, nor by the effect of PAH-specific therapy. Although it allowed us to generate reproducible results in this specific subgroup, it also limited the generalisability of our findings. Further studies should try to investigate the relevance of chemerin as a PVR surrogate in other causes of PH (especially groups 2 and 3), as well as its PAH-specificity in other patient subgroups (notably dcSSc and extensive ILD). Second, although we screened a large number of candidate biomarkers, more recent panels now include several thousand proteins, that were not tested here. Third, as we wanted to obtain reliable hits, we used statistical methods that limited the risk of false discoveries. Due to this stringent approach, some other relevant biomarkers may not have reached statistical significance and thus been missed. Fourth, due to the small sample size and missing data, we could not assess the added value of chemerin serum levels when integrated in PAH diagnostic algorithm and/or risk assessment. This should be addressed in future studies by performing multilevel analysis on larger cohorts. Fifth, due to sample availability considerations, lung experiments included different patient populations, such as iPAH, SSc-ILD and patients with participation of group 3 PH. Finally, as our work had a cross-sectional design, we could not assess the capacity of chemerin level variations to reflect changes in PVR during follow-up. Longitudinal studies are warranted to test its value as a surrogate for treatment efficacy and failure.

Conclusions

Using a wide-scale proteomic approach, we identified chemerin as a potential surrogate for PVR in SSc-PAH patients, which could be interesting as a non-invasive assessment of haemodynamic severity. Chemerin and its receptor CMKLR1 could contribute to the pathogenesis of the disease through PA-SMC proliferation, but could also be involved in inflammatory and fibrotic processes occurring during SSc. Further studies should investigate the potential of this pathway as a marker of disease progression and as a therapeutic target.

Data availability statement

Data are available on reasonable request.

Ethics statements

Patient consent for publication

Ethics approval

This study involves human participants and was approved by institutional review boards from Boston University (IRB# H-31479), Lille University (Comité de Protection des Personnes CPP# 2019-87), Grenoble University (IRB# 6705), Le Kremlin-Bicetre centre (Comité de Protection des Personnes CPP# 8.06.06) and University of Pittsburgh (IRB# PRO14010265, CORID #415 and 718). Written informed consent was obtained from all patients before sample collection. Participants gave informed consent to participate in the study before taking part.

Acknowledgments

The authors wish to acknowledge the contributions of Giuseppina Stifano (Boston University), Eric A. Stratton (Boston University), Barbara Girerd (Le Kremlin-Bicêtre), Nicolas Gonnet (Grenoble), Alain Duhamel (Lille) and Sylvain Dubucquoi (Lille).

References

Supplementary materials

Footnotes

  • Handling editor Josef S Smolen

  • Twitter @Montanidavid

  • SS and LR contributed equally.

  • Contributors All individuals listed as authors met the ICMJE guidelines for determining authorship. SS and RL conceived, designed and supervised the overall study. LR, JCM and RL screened and included patients in the discovery cohort, and contributed samples for the serum experiments. SS, DM and J-LC screened and included patients in the validation cohort, and contributed samples for the serum experiments. SS, LR and JCM performed experiments on serum samples. AMB contributed additional clinical data. SS and LR created and managed the database. CT and GM performed statistical analyses and created Figure 1. SS created figures 2 and 3. EV and RL provided lung samples for the single-cell RNA sequencing. EV performed sequencing experiments, analysed the generated dataset and created Figure 4. LT, DM, MH and CG provided lung samples for the confocal and PASMCs for the proliferation experiments. EH and DL contributed serum samples for the proliferation experiments. LT performed the confocal stainings and PASMC proliferation experiments, and created figure 5 and supplemental figure S1. CT and GM performed statistical analyses. SS created figure 6 and supplemental figure S2. SS wrote the first draft of the manuscript and created the Tables. LT, CG and RL made major revisions to the manuscript. MH provided his PAH expertise. EH and DL provided their SSc expertise. All authors read and approved the submitted version. SS acts as the guarantor of this study.

  • Funding This work was supported by NIH National Institute of Arthritis Musculoskeletal and Skin Disease grants to RL: Scleroderma Core Centers (5P30AR061271), Scleroderma Center of Research Translation (1P50AR060780) and 2R01AR051089.

  • Competing interests SS reports travel grants from Shire, Sanofi-Genzyme, SOBI and Novartis; consulting fees from Novartis; outside the submitted work. J-LC reports grants from United Therapeutic, Bioprojet and Topadur, outside the submitted work. DM reports grants and personal fees from Actelion, grants and personal fees from Bayer, personal fees from GSK, personal fees from Pfizer, personal fees from MSD, personal fees from Chiesi, outside the submitted work. EH reports consulting fees from Actelion, Boehringer Ingelheim, GSK, Roche-Chugai, Sanofi-Genzyme; speaking fees from Actelion, GSK, Roche-Chugai; and research funding from CSL Behring, GSK, Roche-Chugai, outside the submitted work. DL reports personal fees from Actelion, grants and personal fees from Takeda, grants and personal fees from CSL Behring, outside the submitted work. MH reports personal fees from Actelion, grants and personal fees from Bayer, grants and personal fees from GSK, personal fees from Merck, personal fees from United Therapeutics, personal fees from Acceleron, outside the submitted work. RL reports research grants from Formation, Elpidera and Kiniksa, outside the submitted work. RL has served as a consultant for Bristol Myers Squibb, Boehringer-Mannheim, Merck, Magenta and Genentech/Roche, outside the submitted work. Over the last three years, CG, reports grants from Acceleron, ShouTi, and Janssen and grants and personal fees from Merck, outside the submitted work.

  • Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.