Tyrosine kinase inhibitors (TKI) are medications of interest in the treatment of Systemic Sclerosis (SSc) because of their ability to inhibit pathways involved in fibrosis. In this open-label pilot trial, our objectives were to assess the safety, efficacy, and molecular change associated with treatment of patients with diffuse cutaneous (dc)SSc with the TKI nilotinib (Tasigna™).
Ten adult patients with early dcSSc were treated with nilotinib. Primary endpoints were safety and change in modified Rodnan Skin Score (MRSS) after 6 months. Lesional skin biopsies at baseline, 6 and 12 months of treatment were assessed by histopathology, immunohistochemistry, and DNA microarray.
Patients had early and active dcSSc with median disease duration of 0.7 years (range 0.5, 1.7) and increasing MRSS in the month prior to baseline (mean +2.9, p=0.02). Seven out of ten patients completed 6 and 12 months of treatment. Seventy-one adverse events (AEs) including 2 serious AEs were observed, and 92 % of AEs were grade 1-2. Two patients discontinued the medication due to mild QTc prolongation. MRSS improved by a mean of 4.2 points (16 %) at 6 months and by 6.3 points (23 %) at 12 months in the 7 completers, p=0.02 and 0.01, respectively. Patients with a decrease in MRSS >20 % from baseline at 12 months (classified as improvers) had significantly higher expression of transforming growth factor beta receptor (TGFBR) and platelet-derived growth factor receptor beta (PDGFRB) signaling genes at baseline than non-improvers, and the expression of these genes significantly decreased in improvers post-treatment.
Nilotinib was well tolerated by the majority of patients in this study, with tolerability limited primarily by mild QTc-prolongation. Significant MRSS improvement was observed in these early, active patients, but is not conclusive of treatment effect given the open-label study-design and small number of patients in this pilot study. Improvers had higher levels of expression of genes associated with TGFBR and PDGFRB signaling at baseline, and a significant decrease in the expression of these genes occurred only in patients with higher MRSS improvement. The findings of this pilot study warrant more conclusive evaluation.
Systemic sclerosis (SSc; scleroderma) is a multisystem disorder characterized by vasculopathy, autoimmunity, inflammation, and fibrosis . Patients with diffuse cutaneous SSc (dcSSc) have increased morbidity and mortality when compared to patients with other rheumatic diseases . Although several medications are used to treat the skin disease associated with dcSSc, there are no universally effective therapies, and the treatment of scleroderma skin disease remains an area of unmet need .
Transforming growth factor beta (TGFB) and the platelet-derived growth factor (PDGF) are cytokines implicated in the pathological fibrosis of dcSSc [4, 5]. Nilotinib (Tasigna™; Novartis, Basel, Switzerland) is a tyrosine kinase inhibitor (TKI) with antagonistic activity against Abelson tyrosine kinase (c-Abl), the PDGF receptor (PDGFR), and other tyrosine kinases. It is approved in the USA for the treatment of chronic myelogenous leukemia (CML) . It is a therapy of interest for dcSSc because of its ability to interfere with both TGFB and PDGF signaling. Nilotinib has been shown to decrease fibrosis in vitro and in bleomycin models of SSc similarly to imatinib . However, these models have shortcomings in their ability to predict clinical impact in SSc .
Several groups have studied imatinib for the treatment of dcSSc with variable experiences . Unfortunately, none of the studies have been definitive due to the open-label study design, inclusion of patients with limited cutaneous SSc (lcSSc) and morphea, or inadequate power [10–13]. Adverse events (AE), in particular fluid retention, were prominent in these studies, but may be less frequent when imatinib is used at a low dose . Although fluid retention has been observed in other populations, it has been particularly problematic in dcSSc, even leading to the early termination of one trial. Subcutaneous edema may also elevate the modified Rodnan skin score (MRSS), as edema can be difficult to distinguish from dermal thickening.
Nilotinib is a second generation TKI that blocks c-abl and PDGFR (more potently c-abl than imatinib and less potently PDGFR) . In populations with CML, edema has been seen in only 10 % of patients, and this represents an advantage of nilotinib over imatinib as a candidate therapy for SSc. We report here the results of our pilot trial.
Patients were recruited from November 2010 until December 2011. Patients fulfilled the 1980 American College of Rheumatology classification criteria for SSc  and had the diffuse subtype . Patients were over 18 years old, had a disease duration <3 years since the first SSc-related symptom other than Raynaud’s phenomenon, and had a baseline MRSS ≥16. Patients were excluded if they had a baseline corrected QT (QTc) interval on electrocardiogram (EKG) >450 msec. Additional exclusion criteria included treatment with immunosuppressive therapies within 3 months before baseline (including prednisone equivalent >10 mg), pregnancy, serious medical conditions, diffusion capacity of carbon monoxide (DLCO) <30 % predicted, or ejection fraction (EF) <50 %.
This was an investigator-initiated, single-center, open-label pilot study. The primary objective was to assess the safety and tolerability of nilotinib in patients with dcSSc as assessed by the number of AE and serious adverse events (SAE). The primary efficacy endpoint was change in MRSS after 6 months of treatment. Secondary efficacy endpoints included change in MRSS at 12 months, forced vital capacity (FVC) and DLCO on pulmonary function testing (PFT) as well as change in the short form 36 (SF-36) mental (MC) and physical components (PC) and scleroderma health assessment questionnaire disability index (SHAQ-DI). Skin biopsies were assessed using histopathologic and immunohistochemical analysis, and gene expression profiling with DNA microarray to assess change with treatment and explore the biologic basis of the clinical changes observed.
The protocol was approved by the institutional review board at the Hospital for Special Surgery. Patients provided written informed consent before enrollment. An independent Data and Safety Monitoring Board (DSMB) regularly reviewed safety data. The trial was registered at ClinicalTrials.gov (NCT01166139).
Patients were assessed at visits occurring every month for AE ascertainment, interval history, physical examination, clinical laboratory measurements, and 12-lead EKG, and called our center about issues between visits. AE were listed according to the common terminology of the National Cancer Institute . AE were graded as follows: 1 - mild, not requiring intervention; 2 - moderate, requiring minimal local or noninvasive intervention; 3 - severe, 4 - life-threatening; 5 - death. AE were attributed as unrelated, unlikely related, possibly related, probably related, or definitely related to the study medication in the opinion of the investigators and as adjudicated by the Data Safety Monitoring Board (DSMB). The MRSS was measured at screening, baseline, and every 3 months by the same physician (RS or JG). PFT with measurement of FVC and DLCO was performed at baseline and after 6 and 12 months.
Patients started nilotinib 200 mg by mouth daily for 1 week, which was up-titrated to 200 mg twice daily for 3 weeks and then to 400 mg twice daily. A 12-lead EKG with measurement of QTc was checked 1 week after any dosing change. This titration scheme is not required but was chosen by the investigators to carefully observe any dose-related side effects .
Two 3-mm punch biopsies of lesional, extensor-surface, forearm skin were performed at baseline, and after 6 and 12 months of treatment. The post-treatment biopsies were taken 1 cm adjacent to the previous biopsy. At each time point one specimen was formalin-fixed and paraffin-embedded, and the other was stored in RNAlater (Thermo Fisher Scientific, Waltham, MA, USA). Sections for histopathology were stained with hematoxylin and eosin (H&E), anti-α-smooth muscle actin (αSMA), and Masson trichrome (TC), CD-34, pro-Collagen (proCol), phosphorylated cAbl (p-cAbl), and phosphorylated PDGFR (p-PDGFR) using standard techniques. A dermatopathologist (CM), blinded to treatment status, compared each case. Slides were scored semiquantitatively using a scale (0, 1, 2, 3) based on collagen density and degree of infiltrate on H&E, and the intensity of immunohistochemical stains. Skin thickness was measured by a micrometer from the epidermis to the subcutis. Eccrine coils and hair follicles were counted per section.
Gene expression by DNA microarray
Tissue samples stored in RNAlater were homogenized and RNA was purified as previously described . RNA integrity was assessed using the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) and all samples had RNA integrity numbers (RIN). Total RNA (25 ng) was amplified and labeled with Agilent Low Input Quick Amp Labeling Kit. cRNA was hybridized to Agilent SurePrint G3 Human Gene Expression 8x60K Microarrays (G4851A). Agilent Feature Extraction Image Analysis Software (Version 10.7.3) was used to extract data from raw microarray image files. Microarray data were log2-lowess normalized and filtered for probes with intensity ≥1.5-fold over local background in Cy3 or Cy5 channels. Expression values were multiplied by −1 to convert them to log2(Cy3/Cy5) ratios. Probes with >20 % missing data were excluded resulting in 31,762 probes that passed the filtering criteria. Paired samples (baseline and post-treatment biopsies at 12 months) were available for six patients and were used for analyses. We categorized patients as improvers based on a decrease in MRSS of >20 % at 12 months compared to baseline, and 4/6 patients with paired 12-month biopsy specimens met this criterion. The expression data from the study are deposited to NCBI GEO [GEO:GSE65405].
Microarray data pre-processing
Microarray data were pre-processed and analyzed using GenePattern modules  with default parameters unless stated otherwise. Expression data in Stanford preclustering (PCL) file format were converted to Gene Cluster Text (GCT) file format via PclToGct module. Missing expression values were imputed via ImputeMissingValuesKNN module. In total, 31,762 probes were collapsed to 18,398 unique gene symbols via CollapseDataset module using annotation file for the Agilent 8x60K microarray platform. Class label (CLS) files were created via ClsFileCreator module to define phenotypic classes for microarray samples. Expression data were median-centered gene-wise in Cluster 3.0 .
Differential expression analysis
For both baseline (improvers vs non-improvers) and improver (baseline vs post treatment) comparisons, the ComparativeMarkerSelection  module from GenePattern was used to identify differentially expressed genes. The ‘number of permutations’ parameter was set to 0 and the ‘log transformed data’ parameter was set to ‘yes’. For improver comparison, the ‘test statistic’ parameter was set to ‘Paired T-Test’. Expression data for significant genes (p <0.05, not corrected for multiple hypothesis testing) were extracted via the ExtractComparativeMarkerResults module and converted to PCL file format using the GctToPcl module. Expression data were then hierarchically clustered gene-wise and array-wise in Cluster 3.0 using the uncentered correlation similarity metric and average linkage clustering method, and were visualized in TreeView .
Pathway enrichment analysis
For baseline and improver comparisons, pathways with significant changes in expression were identified by gene set enrichment analysis (GSEA) [25, 26] and single-sample GSEA (ssGSEA)  using corresponding GenePattern modules. All GSEA analyses were corrected for multiple hypothesis testing. GSEA and ssGSEA were run against the Canonical Pathways database comprising gene sets from several pathway databases. For GSEA, the ‘permutation type’ parameter was set to ‘gene set’. ssGSEA enrichment scores were normalized by dividing by the maximum ssGSEA enrichment score for this expression dataset. Normalized ssGSEA enrichment scores for significant pathways (false discovery rate (FDR) <5 %) were extracted, clustered and visualized as described above for the expression data.
Intrinsic subset assignment
Intrinsic probes (n = 995) from Milano et al.  were collapsed to 793 unique genes. Separately the entire nilotinib dataset comprising expression data for 24 samples (including all baseline, 6-month and 12-month biopsies) was combined with 4 healthy control samples analyzed on the same DNA microarray platform to provide the proper data distribution across groups. From these samples 27,276 probes passed quality filters and were collapsed to 16,580 unique genes. Overlap with the 793 unique genes from Milano et al. resulted in 651/793 genes (82.1 %) in common between the two datasets. These 651 genes were used to organize the gene expression data from nilotinib and healthy control samples by unsupervised hierarchical clustering.
Intrinsic subset assignment was performed using the 651 intrinsic genes to calculate Spearman non-parametric statistics (correlation coefficients and p values) between each sample from the study and three centroids corresponding to fibroproliferative, inflammatory and normal-like samples from Milano et al. Limited was excluded because no limited SSc samples were included in this study. Centroids were created by averaging expression values for each gene across all samples assigned to a given intrinsic subset in Milano et al. The intrinsic subset assignment for each nilotinib sample was to the subset centroid with the highest Spearman correlation coefficient and the lowest p value.
Descriptive analysis of AE and SAE was performed for assessment of clinical data. Continuous clinical variables were analyzed using the unpaired (baseline vs baseline) and paired (baseline vs post treatment) t test.
For gene expression analysis, the relationship between baseline intrinsic subset and response status was analyzed using the chi-square test. In cases where three or more groups were compared, one-way analysis of variance (ANOVA) was used and p values were corrected for multiple testing using Tukey’s multiple comparisons test with a single pooled variance. Statistical analyses were performed using SPSS software version 17.0 (IBM, Armonk, NY, USA) and GraphPad Prism Windows 6.05 (GraphPad Software, San Diego, CA, USA).
Nineteen patients were screened for study inclusion. Five were excluded due to baseline QTc >450 msec and four were excluded due to meeting other exclusion criteria. Baseline characteristics of all patients are listed in Table 1. The patients in this study represented a group of patients with early and active dcSSc. The median disease duration based on the time since the first non-Raynaud’s symptom of SSc was 0.7 years (range 0.5−1.7). In the 1-month screening period prior to starting medication, 70 % patients had an increasing MRSS, and the group had a mean MRSS increase of 2.9 points, p = 0.02. Fifty percent of patients were positive for the RNA Polymerase 3 antibody. Seven patients tolerated 6 months of nilotinib and continued treatment for at least 12 months. Two patients had to discontinue nilotinib within the first month due to grade 1 (QTc of 453 msec) or 2 (QTc of 483 msec) QTc prolongation. One patient discontinued due to progression of preexisting coronary artery disease after 3 months of treatment.
Seventy-one AE including two SAE were observed during the 12-month treatment period (Additional file 1); 75 % were considered to be at least possibly related to nilotinib (Table 2). Of the AE, 92 % were grade 1 or 2 and no unexpected events were observed. One common AE was grade 1 (<2.5 × upper limit of normal) liver function test (LFT) abnormality, which occurred in 50 % of patients. All LFT abnormalities are shown in Additional file 1, and all were asymptomatic. No medication adjustment was necessary and the abnormal values resolved when the laboratories were rechecked. There were no LFT abnormalities of grade 2 or higher. QTc prolongation occurred in 60 % of patients. Although QTc prolongation was grade 1 or 2, our protocol required discontinuation of two patients due to this. The other four patients had normalized QTc on repeat EKG or adjustment of medication dosage.
Two SAEs occurred both of which were considered possibly related to the medication. One patient experienced a syncopal episode and was hospitalized. The medication was held and the patient’s evaluation, including EKG with QTc measurement and telemetry, was normal with the syncopal episode considered to be a vasovagal episode. However, when the patient followed up and off medication for 3 weeks, her baseline was QTc >450 msec and she could not restart medication as she now met an exclusion criterion. Another patient with known prior history of coronary artery disease (CAD) was hospitalized for coronary artery bypass graft (CABG) surgery after 3 months of treatment. The patient had multiple risk factors for vascular disease including age, diabetes mellitus, longstanding essential hypertension, hyperlipidemia, and obesity. She had atypical chest pain, which she described as heartburn and which had been attributed to her gastroesophageal reflux disease predating the trial. She underwent cardiac testing during the trial, which in turn led to indication for CABG. It is possible that worsening of her vascular disease was related to nilotinib use, although with her multiple risk factors and symptoms that predated the use of nilotinib, the attribution of this SAE was not clear cut.
In the seven patients who completed 12 months of treatment, the MRSS statistically significantly improved from a mean of 26.9 ± 5.4 to 22.7 ± 8 (by 16 %) at 6 months and to 20.6 ± 7.7 (by 23 %) at 12 months, p = 0.02 and 0.0112, respectively (Table 3). The mean change in the MRSS in the completers with the RNA Polymerase3 antibody (n = 4) was −7.0 ± 1.8 versus −5.5 ± 7.5 in those without (n = 3). This was not statistically significant, p = 0.68.
Two patients who discontinued study medication had follow-up MRSS at 6 months. One patient’s MRSS decreased from 43 to 21 on 4 months of mycophenolate, after stopping nilotinib after less than 1 month due to QTc prolongation. One patient’s MRSS was relatively unchanged from 44 to 46 after stopping nilotinib after 3 months. The third patient continued her care locally.
Out of 10 patients, 3 had ILD at baseline and 2 of those patients did not complete the study. The mean FVC was 77.4 ± 12.9 % predicted at baseline, 75.4 ± 12.6 % at 6 months and 71.7 ± 11.7 at 12 months, p = 0.17 and 0.054, respectively. The DLCO was 72.0 ± 9.9 % predicted at baseline and 69.9 ± 18.1 % at 6 months, and 69.3 ± 12.4 % at 12 months, p = 0.56 and 0.25, respectively. The completer with ILD had stable PFT parameters. The trend to decline in FVC was driven by one patient without evidence of ILD at entry, who went on to develop mild reticular opacities after one year.
The physician global assessment (PGA) significantly improved from 62.5 ± 10.8 at baseline to 40.2 ± 12.7 at 6 months, and to 32.3 ±15.3 at 12 months, p = 0.0033 and 0.0013, respectively. There were no significant differences in ESR, SHAQ-DI, SF-36 mental or physical components, oral aperture, hand extension, or finger-to-palm distance.
Eight patients had biopsies at baseline and 6 months and six patients had biopsies at 12 months. The mean skin thickness at baseline was 2.3 ± 0.6 mm, 2.7 ± 0.6 mm at 6 months, and 2.6 ± 0.3 mm at 12 months, with no statistically significant difference (p = 0.06 at 6 months and p = 0.36 at 12 months). One patient sample is shown in Fig. 1. There was no significant change in the group in collagen density, degree of infiltrate on H&E, number of follicles and eccrine structures, and staining intensity of α-SMA, trichrome, CD-34, pro-Collagen, p-PDGFR, or p-cAbl; these data are summarized in Additional file 2. There was no difference in qualitatively assessed staining of pPDGFR or p-cAbl staining between those patients with MRSS improvement >20 % and those who did not meet this level of improvement.
Gene expression changes in skin
Intrinsic gene expression subset assignment before and after treatment
We used intrinsic subset genes from Milano et al.  (see Methods and Additional file 3) to assign the samples in this trial to intrinsic gene expression subsets (Fig. 2a). Subsets were assigned based on Spearman correlation statistics to a gene expression centroid calculated from Milano et al. samples (see Methods and Additional file 4). Clustering was primarily driven by a strong inflammatory signature evident in the skin of these patients (Fig. 2b). Baseline intrinsic subset assignment showed that three out of four improvers were classified as non-fibroproliferative at baseline (two inflammatory and one normal-like), whereas one improver was classified as fibroproliferative. The patients who improved generally lost their major subset signature and three out of four were classified as normal-like post treatment (Table 4). The two non-improvers were assigned to the fibroproliferative subset (p = 0.0833, chi-square test for enrichment of fibroproliferative patients as non-improvers) and remained in the fibroproliferative subset post treatment, consistent with the results of Pendergrass et al. . Slight discrepancies between the Spearman correlation intrinsic subset assignments and the hierarchical clustering tree occured mainly in post-treatment samples.
TGFBR and PDGFRB signaling pathways are highly expressed in improvers at baseline
We performed differential gene expression analysis between baseline improver (n = 4) and non-improver (n = 2) samples (Fig. 3a): 2,242 genes were significantly differentially expressed at baseline (p <0.05, unpaired t test) (Fig. 3b).
We used GSEA to identify 18 pathways with significantly differential expression between improvers and non-improvers (FDR <5 %; corrected for multiple testing). Pathways that increased in improvers included TGFB receptor signaling (TGFBR) and nuclear factor kappa B (NFKB) signaling, whereas pathways with increased expression in non-improvers included telomere maintenance and meiosis, suggesting enrichment in terms associated with the cell cycle and cell proliferation, consistent with the assignment of the two non-improvers to the fibroproliferative subset (Fig. 3c).
The four improvers with baseline and 12-month biopsies were included in the microarray data analysis (Fig. 4a). We identified 666 genes significantly differentially expressed in improvers between baseline and post treatment (p <0.05, paired t test) that comprised an improver gene signature (Additional file 5). Genes that showed decreased expression post treatment included TGFB-regulated (SMAD2, ACTB, COL15A1, ABI1, EGR2 and EGR3) and inflammatory (SIGLEC7, IL23A, ICAM1, CCL2 and MMP14) gene signatures (Fig. 4b, Additional file 5). These genes had significantly higher expression in improvers at baseline relative to non-improvers, and were significantly decreased in improvers post treatment (p <0.0001 for both comparisons). In contrast, their expression remained stable in non-improvers (p = 0.1057) (Fig. 4c). Non-improvers had 74 genes with significant differential expression between baseline and post-treatment samples, and no significant overlap with the improver signature. There were only 2/74 genes in common with the improver gene signature but their directionality was opposite to that for improvers.
We identified 25 pathways that were highly significantly decreased in improvers post-treatment (FDR <1 %), including multiple immune response pathways such as interleukin and interferon signaling (Fig. 4d). A broader list of 106 pathways that were significantly decreased in improvers post treatment (FDR <5 %) is shown in Additional files 6 and 7. We were interested in TGFBR and PDGFRB pathways because of their relevance to fibrosis in SSc and their convergence on c-Abl, a target of nilotinib. We identified 19 genes involved in TGFBR signaling (including TGFB2, TGFB3, TGFBR1 and TGFBR2) and 50 genes involved in PDGFRB signaling (including PDGFRB and PDGFB) that comprised the core enrichment groups for these pathways in improvers based on GSEA results (Additional file 8).
We examined the expression of TGFBR and PDGFRB signaling gene sets across six completers (combined improvers and non-improvers), four improvers and two non-improvers at baseline and post treatment (Fig. 5). Both pathways were significantly higher at baseline in improvers compared to non-improvers (p <0.0001 for both gene sets). Both pathways significantly decreased post treatment in improvers (p <0.0001 for both gene sets). While genes associated with TGFBR signaling showed stable expression in non-improvers (p = 0.1092) (Fig. 5a), PDGFRB signaling genes were significantly increased in non-improvers post treatment (p = 0.0058; Fig. 5b). No pathways were significantly differentially expressed in samples from non-improvers at baseline and post-treatment.
The initial positive in vitro and murine studies on the use of TKIs with c-Abl and PDGFR inhibition (e.g., imatinib, nilotinib, and dasatinib) for dcSSc presented the hope for a profound treatment effect for the cutaneous fibrosis seen in dcSSc. Several studies have attempted to address the safety and effectiveness of imatinib in this population, but design limitations make it impossible to make definitive conclusions. The collective imatinib experience has been mixed. A modest improvement in skin thickening was seen in some open-label studies, but side effects, especially fluid retention, were prevalent and prevented some patients from continuing the medication. This paper is the first published report on the use of nilotinib for the treatment of dcSSc.
In this pilot study of nilotinib, a different but related TKI with fewer fluid-related side effects, our goals were to assess safety in a population with dcSSc, generate pilot data on the MRSS to help us adequately power future studies, and to use histopathology and gene expression profiling to further understand our clinical observations.
Nilotinib was well-tolerated by the majority of patients in this study. Although the number of AE that were at least possibly related to nilotinib was high, many of the AE represented asymptomatic laboratory abnormalities that resolved without intervention. The side effects observed were similar to what has been seen in patients with CML [30, 31]. However, nilotinib is known to have important cardiovascular side effects, which were observed and are particularly relevant to SSc patients. Monitoring the QTc with periodic EKG is required. Conduction abnormalities including QTc prolongation can be seen in dcSSc patients  or as a side effect of other medications they may use. In our study QTc prolongation was seen in five SSc patients at screening, which limited their eligibility. Persistent mild QTc prolongation without arrhythmia was the cause of two patients discontinuing treatment. Acceleration of peripheral arterial vascular disease has been reported with nilotinib in CML patients both with and without other risk factors for CAD [33, 34]. In this small study we observed exacerbation of preexisting CAD in one patient with other risk factors. Careful consideration of the use of this medicine in patients with CAD or multiple risk factors would be needed in patients with SSc, as it would in CML.
We observed that the MRSS improved in a statistically significant manner with 6 and 12 months of treatment in this group of patients with very early and active SSc. The range of the minimally clinically important difference in MRSS has been estimated to be 3.2−5.3 points . Therefore, the patients in this study were observed to change in a clinically meaningful manner. At the time of study entry the median disease duration was 0.7 years, 70 % of patients had a worsening MRSS, 50 % had the RNA Polymerase3 antibody and 40 % had tendon friction rubs. The combination of these clinical attributes put these patients at high risk of exacerbation of disease during the study period. However, the MRSS improvement is inconclusive given the uncontrolled study design and the observation that the MRSS tends to improve both with time  and in the context of clinical trials . Nonetheless, this clinical finding of improved MRSS in group of patients with such early and active dcSSc is worthy of further and more definitive study.
The FVC and DLCO were not significantly changed, but our trial was not designed to evaluate pulmonary endpoints, and most patients in this study did not have ILD. Recently another TKI, nintedanib, which has some overlapping spectrum of TK inhibition, has been shown to slow the decline of FVC in idiopathic pulmonary fibrosis, a related fibrotic condition . The additional outcome measures, including the SF-36 and SHAQ, at the 6- and 12-month time points also remained stable, and this has been shown to occur typically in SSc clinical trials .
In this study we used gene expression profiling and the concept of SSc intrinsic subsets to characterize biological response and distinguish between nilotinib improvers and non-improvers. In terms of the intrinsic subset assignment, three of four improvers were classified as non-fibroproliferative (with two out of four classified as inflammatory), whereas non-improvers were classified as fibroproliferative at baseline. Improvers tended to lose their intrinsic subset signature and become normal-like post treatment, while non-improvers remained fibroproliferative (Table 4). Genes from specific pathways targeted by nilotinib (e.g., TGFBR and PDGFRB signaling) showed significantly higher expression in improvers at baseline and were significantly downregulated in those improvers by treatment, whereas their expression in non-improvers was stable or displayed directionality opposite to improvers, with low expression at baseline. Interestingly, a recent analysis of imatinib and nilotinib effects in SSc mouse models suggested that high activation status and expression pattern of TKI targets was predictive of the potential response to therapy . It is not clear based on this single-group and uncontrolled clinical trial whether our results represent the natural history of the disease or a treatment effect, although an apparent change in the post-treatment intrinsic subset assignment to the one with less severe disease in improvers might reflect beneficial effects of the therapy. As the use of gene expression in the context of clinical trials continues, the significance of these findings will be more conclusive.
In this study, we found that TGFBR signaling was significantly increased in improvers at baseline and that this signature spanned the inflammatory and fibroproliferative subsets. In our original study of TGFB signaling in SSc skin  we observed its enrichment primarily in the fibroproliferative subset. However, a recent meta-analysis of three independent gene expression datasets from SSc skin, which dramatically increases the number of samples analyzed, shows that TGFBR signaling appears to span the inflammatory and fibroproliferative subsets [41, 42].
The results here are consistent with those data, as our improvers, who are primarily defined by high expression of TGFBR signaling at the initiation of treatment, appear to span across the fibroproliferative and inflammatory intrinsic gene expression subsets. This indicates that a biomarker of TGFBR signaling will be a useful stratification method in addition to biomarkers of the intrinsic gene expression subsets.
The histopathology of the forearm skin biopsy did not show a significant morphological improvement even in cases where a clinical improvement was seen, and phospho-staining of known targets of nilotinib in the skin was not differentially activated before and after treatment. It is possible that these dermatopathological investigations are less sensitive than microarray in picking up this degree of change, and this would be a topic for future investigation.
There were several limitations to this study, including its open-label design and small sample size. Because of these factors conclusions about efficacy cannot be definitively drawn. As noted, the MRSS tends to improve due to regression to the mean  and in the context of clinical trials . The lack of a control group additionally makes the changes in gene expression observed in the microarray studies inconclusive with respect to whether they represent a change due to treatment or a change based on the natural history of the disease. In terms of the differential gene expression analyses, we had to use uncorrected p values in the unpaired and paired t tests as the small sample size precluded the adjustment for multiple comparisons. We controlled for this by focusing on pathways found to be deregulated using GSEA, which were corrected for multiple hypothesis testing. This was a pilot study designed to provide preliminary data that could drive more conclusive hypothesis testing.
This was an open-label, single-group pilot trial, which included histopathologic and microarray analyses of skin to further delineate biological response. In this pilot clinical trial a significant improvement in the MRSS was observed in patients who would have been expected to have progression of disease. Seventy percent of patients were able to tolerate the medication for 12 months, and side effects were predictable. Improvers were most clearly defined by having high expression of genes associated with TGFBR/PDGFRB signaling at the time of treatment initiation, while those that did not improve did not have evidence of activation. TGFBR and PDGFRB signaling pathways significantly decreased in improvers post treatment and this was not observed in non-improvers. Further evaluation of nilotinib or related TKIs in a randomized and controlled setting is warranted.
Abelson murine leukemia viral oncogene homolog
coronary artery disease
chemokine (C-C motif) ligand 2
chronic myelogenous leukemia
diffuse cutaneous systemic sclerosis
diffusion lung capacity of carbon monoxide
early growth response protein
erythrocyte sedimentation rate
false discovery rate
forced vital capacity
gene set enrichment analysis
hematoxylin and eosin
intercellular adhesion molecule 1
interstitial lung disease
limited cutaneous systemic sclerosis
matrix metalloproteinase 14
modified Rodnan skin score
nuclear factor kappa-light-chain-enhancer of activated B cells
platelet-derived growth factor
platelet-derived growth factor receptor
pulmonary function test
physician global assessment
corrected QT interval - referring to an interval on an electrocardiogram
serious adverse event
short form 36-item health survey
short form 36-item health survey mental component
short form 36-item health survey physical component
scleroderma health assessment questionnaire - disability index
sialic acid-binding immunoglobulin-type lectin 7
transforming growth factor beta
transforming growth factor beta receptor
tyrosine kinase inhibitor
alpha smooth muscle actin
Varga J, Abraham D. Systemic sclerosis: a prototypic multisystem fibrotic disorder. J Clin Invest. 2007;117:557–67.
Thomas E, Symmons DP, Brewster DH, Black RJ, Macfarlane GJ. National study of cause-specific mortality in rheumatoid arthritis, juvenile chronic arthritis, and other rheumatic conditions: a 20 year followup study. J Rheumatol. 2003;30:958.
Kantarjian HM, Giles F, Gattermann N, Bhalia K, Alimena G, Palandri F, et al. Nilotinib (formerly AMN107), a highly selective Bcr-Abl tyrosine kinase inhibitor, is effective in patients with Philadelphia chromosome-positive chronic myelogenous leukemia in chronic phase following imatinib resistance and intolerance. Blood. 2007;110:3540–6.
Akhmetshina A, Dees C, Pileckyte M, Maurer B, Axmann R, Jüngel A, et al. Dual inhibition of c-abl and PDGF receptor signaling by dasatinib and nilotinib for the treatment of dermal fibrosis. FASEB J. 2008;22:2214–22.
Spiera RF, Gordon JK, Mersten JN, Magro CM, Mehta M, Wildman HF, et al. Imatinib mesylate (Gleevec) in the treatment of diffuse cutaneous systemic sclerosis: results of a 1-year, phase IIa, single-arm, open-label clinical trial. Ann Rheum Dis. 2011;70:1003–9.
Khanna D, Saggar R, Mayes MD, Abtin F, Clements PJ, Maranian P, et al. A one-year, phase I/IIa, open-label pilot trial of imatinib mesylate in the treatment of systemic sclerosis-associated active interstitial lung disease. Arthritis Rheum. 2011;63:3540–6.
Pope J, McBain D, Petrlich L, Watson S, Vanderhoek L, de Leon F, et al. Imatinib in active diffuse cutaneous systemic sclerosis: Results of a six-month, randomized, double-blind, placebo-controlled, proof-of-concept pilot study at a single center. Arthritis Rheum. 2011;63:3547–51.
Fraticelli P, Gabrielli B, Pomponio G, Valentini G, Bosello S, Riboldi P, et al. Low-dose oral imatinib in the treatment of systemic sclerosis interstitial lung disease unresponsive to cyclophosphamide: a phase II pilot study. Arthritis Res Ther. 2014;16:R144.
Subcommittee for Scleroderma Criteria of the American Rheumatism Association Diagnostic and Therapeutic Criteria Committee Preliminary criteria for the classification of systemic sclerosis (scleroderma). Arthritis Rheum 1980;23:581–90.
LeRoy EC, Medsger Jr TA. Criteria for the classification of early systemic sclerosis. J Rheumatol. 2001;28:1573–6.
Hinchcliff M, Huang CC, Wood TA, Matthew Mahoney J, Martyanov V, Bhattacharyya S, et al. Molecular signatures in skin associated with clinical improvement during mycophenolate treatment in systemic sclerosis. J Invest Dermatol. 2013;133:1979–89.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–50. Epub 2005 Sep 30.
Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, et al. PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003;34:267–73.
Kantarjian HM, Giles FJ, Bhalla KN, Pinilla-Ibarz J, Larson RA, Gattermann N, et al. Nilotinib is effective in patients with chronic myeloid leukemia in chronic phase after imatinib resistance or intolerance: 24-month follow-up results. Blood. 2011;117:1141–5.
Medsger T. Natural History of systemic sclerosis and the assessment of disease activity, severity, functional status, and psychologic well-being. Rheum Dis Clin N Am. 2003;255–273.
Amjadi S, Maranian P, Furst DE, Clements PJ, Wong WK, Postlethwaite AE, et al. Course of the modified Rodnan skin thickness score in systemic sclerosis clinical trials: analysis of three large multicenter, double-blind, randomized controlled trials. Arthritis Rheum. 2009;60:2490–8.
Merkel PA, Silliman NP, Clements PJ, Denton CP, Furst DE, Mayes MD, et al. Patterns and predictors of change in outcome measures in clinical trials in scleroderma: an individual patient meta-analysis of 629 subjects with diffuse cutaneous systemic sclerosis. Arthritis Rheum. 2012;64:3420–9.
Sargent JL, Milano A, Bhattacharyya S, Varga J, Connolly MK, Chang HY, et al. A TGFbeta-responsive gene signature is associated with a subset of diffuse scleroderma with increased disease severity. J Invest Dermatol. 2010;130:694–705.
Mahoney JM, Taroni J, Martyanov V, Wood TA, Greene CS, Pioli PA, et al. Systems level analysis of systemic sclerosis shows a network of immune and profibrotic pathways connected with genetic polymorphisms. PLoS Comput Biol. 2015;11(1):e1004005.
The authors would like to thank Yifang Li for processing dermatopathology specimens, Christopher Hatzis, Morgana Davids, and Uzunma Udeh for their role as study coordinators. This study was supported by an investigator-initiated grant from Novartis and by the Rudolph Rupert Scleroderma Program at the Hospital for Special Surgery.
Authors and Affiliations
Department of Rheumatology, Hospital for Special Surgery, 535 East 70th St, New York, NY, 10021, USA
Jessica K. Gordon, Wei-Ti Huang, Mary K. Crow & Robert F. Spiera
Geisel School of Medicine at Dartmouth, Department of Genetics, Remsen 7400, Hanover, NH, 03755, USA
Viktor Martyanov, Tammara A. Wood & Michael L. Whitfield
Weill Cornell Medical Center, Department of Dermatopathology, 525 East 68th St, New York, NY, 10065, USA
Weill Cornell Medical Center, Department of Dermatology, 1305 York Ave, New York, NY, 10021, USA
This study was supported by an investigator-initiated grant from Novartis and by the Rudolph Rupert Scleroderma Program at the Hospital for Special Surgery. JKG has been supported by a New Investigator Grant from the Scleroderma Foundation and by the Clinician Scientist Development Award through the Kellen Foundation at the Hospital for Special Surgery. TAW, VM and MLW were supported by grants from the Scleroderma Research Foundation and by NIH P30AR061271 from the National Institute of Arthritis Musculoskeletal and Skin Diseases (NIAMS) the National Institutes of Health. RFS has received grant funding from Novartis, Human Genome Sciences, Roche, Actelion, GlaxoSmithKline, Genentech, InterMune, and NIH 5U19AI056363-10 from NIAMS. None of the authors report any non-financial competing interests.
JKG participated in the study design, acquisition of clinical data, and data analysis, and drafted the manuscript. RFS conceived of the study, participated in the study design, acquisition of clinical data, and data analysis, and revised the manuscript. TAW extracted RNA from the skin samples, ran the microarray studies, participated in the microarray analysis and participated in drafting the manuscript. VM performed the microarray analysis and participated in drafting the manuscript. MLW performed the microarray studies and analysis and participated in drafting the manuscript. WH performed the clinical statistical analysis and helped to revise the manuscript. HFW contributed to study design and acquisition of skin samples and participated in drafting the manuscript. CM contributed to the dermatopathologic aspect of study design, performed all of the dermatopathologic readings and participated in drafting the manuscript. MKC contributed to study design and interpretation of the data and participated in drafting the manuscript. All authors read and approved the final manuscript.
All adverse events. This table lists all adverse events recorded in order of frequency. Grades listed are based on definitions set forth in CTCAE v2.0. Attribution as per the investigators and depending on the clinical scenario is listed as well. 1 signifies not related to the medication, 2 signifies unlikely related to the medication, 3 signifies possibly related to the medication, 4 signifies probably related to the medication, and 5 signifies definitely related to the medication. Some side effects have different attribution for different events given clinical context. QTC corrected QT interval, ESR erythrocyte sedimentation rate, CABG coronary artery bypass graft. (DOCX 17 kb)
List of 651 intrinsic genes used to cluster nilotinib samples in Fig. 2. Gene symbols are color-coded based on their association with inflammatory (purple), normal-like (green) or fibroproliferative (red) intrinsic subsets. (XLSX 32 kb)
Detailed output for the 666 gene signature in improvers from Fig. 4b. For each gene in the improver gene signature, the following information is provided: where it is upregulated (baseline or post treatment), gene symbol, gene name, raw p value, false discovery rate (FDR)(BH), q value, fold change, baseline mean expression value and post-treatment mean expression value. (XLSX 73 kb)
Heat maps for core enrichment genes from transforming growth factor beta receptor (TGFBR) (A) and platelet-derived growth factor receptor beta (PDGFRB) (B) signaling pathways corresponding to Fig. 5. (PDF 316 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Gordon, J.K., Martyanov, V., Magro, C. et al. Nilotinib (Tasigna™) in the treatment of early diffuse systemic sclerosis: an open-label, pilot clinical trial.
Arthritis Res Ther17, 213 (2015). https://doi.org/10.1186/s13075-015-0721-3