Skip to main content


We’d like to understand how you use our websites in order to improve them. Register your interest.

Limited cutaneous systemic sclerosis skin demonstrates distinct molecular subsets separated by a cardiovascular development gene expression signature



Systemic sclerosis (SSc; scleroderma) is an uncommon autoimmune rheumatic disease characterised by autoimmunity, vasculopathy and fibrosis. Gene expression profiling distinguishes scleroderma from normal skin, and can detect different subsets of disease, with potential to identify prognostic biomarkers of organ involvement or response to therapy. We have performed gene expression profiling in skin samples from patients with limited cutaneous SSc (lcSSc).


Total RNA was extracted from clinically uninvolved skin biopsies of 15 patients with lcSSc and 8 healthy controls (HC). Gene expression profiling was performed on a DNA oligonucleotide microarray chip. Differentially expressed genes (DEG) were identified using significance analysis of microarrays (SAM). Functional enrichment analysis of gene signatures was done via g:Profiler.


There were 218 DEG between lcSSc and HC samples (false discovery rate <10%): 181/218 DEG were upregulated in lcSSc samples. Hierarchical clustering of DEG suggested the presence of two separate groups of lcSSc samples: “limited 1” and “limited 2”. The limited-1 group (13 samples, 10 unique patients) showed upregulation of genes involved in cell adhesion, cardiovascular system (CVS) development, extracellular matrix and immune and inflammatory response. The CVS development signature was of particular interest as its genes showed very strong enrichment in response to wounding, response to transforming growth factor (TGF)-β and kinase cascade. Neither limited-2 samples (six samples, five unique patients) nor HC samples showed functional enrichment. There were no significant differences in demographic or clinical parameters between these two groups. These results were confirmed using a second independent cohort.


Our study suggests the presence of molecular subsets in lcSSc based on gene expression profiling of biopsies from uninvolved skin. This may reflect important differences in pathogenesis within these patient groups. We identify differential expression of a subset of genes that relate to CVS and are enriched in fibrotic signalling. This may shed light on mechanisms of vascular disease in SSc. The enrichment in profibrotic profile suggests that dysregulated gene expression may contribute to vasculopathy and fibrosis in different disease subsets.


Systemic sclerosis (SSc; scleroderma) represents a major clinical challenge and offers insight into fundamental processes relating to autoimmunity, fibrosis and vascular injury and pathology. It is an uncommon immune-mediated rheumatic disease with a very high clinical burden, high mortality and limited treatment options [1, 2] and provides a model for more common forms of organ-based fibrosis in the lung, liver and kidney. Recent observational cohort studies have highlighted remarkable clinical diversity in terms of pattern and extent of skin and internal organ involvement, clinical outcome and response to therapy [3,4,5]. Most current and emerging treatment strategies focus on intensive immunosuppression though greater understanding of the biology of disease outcomes, in particular mechanisms that determine improvement in skin and organ-based disease, may help to personalize treatment strategies more effectively [6].

The limited cutaneous subset of SSc (lcSSc) is characterised by less severe and extensive skin fibrosis but patients can develop major internal organ complications and the vascular manifestations of SSc, particularly pulmonary arterial hypertension and digital ulceration, are prominent in this subset [3]. Detailed gene expression analysis in SSc biopsies has recently been used to define molecular intrinsic subsets of the disease and to provide mechanistic insight into the pathobiology [7,8,9,10,11,12,13,14]. Interestingly, clinically uninvolved skin in the more extensive diffuse cutaneous subset of SSc has closely replicated gene expression signatures compared with biopsies from clinically involved skin [7] and a number of vasoactive genes have been identified in this way in skin from patients with the limited subset, who cluster separately from healthy controls and from patients with diffuse disease [7]. It has been shown that lesional and non-lesional lcSSc biopsies consistently cluster together and show concordance in their deregulated pathways [7, 14]. The numbers of patients with limited disease in those studies are small but the results are consistent. On this basis we sought to develop greater understanding about the gene expression abnormalities in the more prevalent lcSSc subset through detailed transcriptional analysis of skin biopsies taken from uninvolved forearm skin. We hypothesise that this may give valuable insight into the key pathogenetic processes underlying the disease and also provide potential for defining and characterising molecular subsets of lcSSc. Since the molecular subsets of SSc may also inform clinical decision-making and treatment selection, there may be additional value from extending this concept more broadly into lcSSc.

In this study, we demonstrate differential gene expression in a cohort of lcSSc patients and healthy controls and describe a distinct lcSSc subgroup not identifiable by clinical or serological assessment showing enrichment in cell adhesion, cardiovascular system (CVS) development and extracellular matrix genes. We confirm our findings in a second independent cohort of samples. The CVS development signature was significantly different between both subgroups of limited patients and a control group.


Inclusion criteria and study participants

Demographic information, clinical history including organ involvement, other diagnoses and extent of skin involvement, autoantibody profiles and treatments were retrieved from medical notes or obtained at the time of skin biopsy. Diagnosis of SSc was made according to the 2013 American College of Rheumatologists (ACR)/European League Against Rheumatism (EULAR) classification criteria and assignment to the limited cutaneous subset reflected the distribution and extent of skin thickening at the time of biopsy together with other typical disease characteristics [15]. Autoantibodies were measured in an accredited institutional autoimmune serology laboratory using validated commercial tests with appropriate quality control and blinded assessment of the results at the time of biopsy. In brief, antinuclear antibody (ANA) pattern was screened by indirect immunofluorescence using an HEp-2 cell substrate and further characterisation of defined extractable nuclear antigen (ENA) was by counterimmunoelectrophoresis for anti-ENA using soluble extracts from human spleen and rabbit thymus acetone powder (Pelfreez Biologicals, Rogers, AR, USA) as antigen. All patients included in this study signed informed consent for their clinical and laboratory data to be used in this clinical research project.

Sample collection

Forearm 4-mm skin punch biopsies were performed both in patients with lcSSc and in healthy controls, and were stored for extraction in RNAlaterTM RNA stabilisation reagent (Ambion, Austin, TX, USA) at 4 °C overnight followed by longer-term storage at -70 °C. To avoid a potential confounding effect of different degrees of skin thickening or fibrosis, biopsy sites were selected from skin that was not clinically involved. There were no changes in pigmentation or dryness.

Microarray processing

Tissue homogenization was performed using Qiagen TissueLyser II. RNA purification was carried out in a QIAcube with Qiagen RNeasy Fibrous Tissue Mini Kit (Qiagen, Gaithersburg, MD, USA). The Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA, USA) was used to assess the RNA integrity of samples with numbers >7. RNA concentration was measured with Thermo Scientific NanoDrop 2000 Spectrophotometer (Wilmington, DE, USA). A measure of 200 ng of total RNA was amplified and labelled with Agilent's Quick-Amp Labelling Kit. Cy3-labelled samples and Cy5-labelled Universal Human Reference RNA (Stratagene, La Jolla, CA, USA) were co-hybridized to Agilent's Human Genome (4 × 44 K) Microarrays (G4112F). Data were log2 lowess-normalised and filtered for probes with intensity ≥1.5-fold over local background in Cy3 or Cy5 channels. Data were multiplied by -1 to convert them to log2(Cy3/Cy5) ratios. Probes with >20% missing data were excluded.

Gene expression data pre-processing

Expression data were pre-processed using GenePattern [16] modules with default settings unless stated otherwise. Missing values were imputed using the ImputeMissingValuesKNN module. Expression data were collapsed from probes to unique genes using the CollapseDataset module with the Agilent 4 × 44 K chip platform. As microarrays were processed in three separate batches, batch bias was adjusted for using GenePattern implementation of ComBat [17] by means of parametric prior method and information about lcSSc samples and controls as a covariate. Batch bias was assessed before and after ComBat using guided principal component analysis (gPCA) [18]. Expression data were adjusted by median-centering genes in Cluster 3.0 [19].

Differential expression and functional enrichment analysis

Differentially expressed genes were identified using significance analysis of microarrays (SAM) [20]. Two-class unpaired response type was used for comparing two groups and multiclass response type was used for comparing three groups. The number of permutations was set to 500. Expression data for significantly differentially expressed genes were hierarchically clustered in Cluster 3.0 and visualized in Java TreeView [21], version 1.1.6r4. Functional enrichment analysis of gene signatures was done using g:Profiler [22] with the following settings: maximum size of functional category was set to 3500, default multiple testing correction method (g:SCS significance threshold) was used and regulatory motif and protein-protein interaction databases were excluded from the analyses.


Subject selection and clinical characteristics

There were 23 subjects included in the discovery cohort: 15 patients with lcSSc and 8 healthy controls. The study samples were collected and analysed in two stages, first in 10 biopsies from patients with lcSSc and from 5 healthy controls that generated the initial CVS signature, and then in samples from an extended cohort with an additional 6 patients with lcSSc and 3 healthy controls. These later samples were collected independently of the first set of samples.

Demographics, clinical characteristics, organ-based disease and therapies are listed in Table 1. The demographics of the control group were broadly similar to patients in terms of age, sex and ethnicity and these are also representative of the single-centre cohort that included the sample population. All patients had established lcSSc with minimum disease duration of 2 years, and most had had the condition for more than 10 years at the time of biopsy. Patients included in the study retained the subset designation according to their early stage disease and therefore no patients who would previously have been designated as having ‘diffuse’ SSc or those who had higher skin scores in the past were included. Serological findings were again representative of a cohort of lcSSc patients but kept purposefully broad: six patients had centromere pattern staining on indirect immunofluorescence (IIF) for ANA, and four had scl-70 reactivity. Two patients had overlap Ro antibodies and two were ANA-negative; one had serological features of lupus with high double-stranded (ds) DNA antibodies and low C4. There were 6/15 patients with digital ulcers at the time of biopsy and the majority took therapy for a vascular complication, either for significant Raynaud’s phenomenon or for digital ulceration. No patients from this cohort had a diagnosis of pulmonary arterial hypertension and none had had a scleroderma renal crisis or myositis. Four had interstitial lung disease and three (more than may be expected) had an overlap diagnosis of inflammatory arthritis: these two latter diagnoses resulted in the significant use of immunosuppressive medications compared with most cohorts of lcSSc patients. The broad but characteristic range of clinical and serological features demonstrated in this cohort allowed us to identify pathogenic factors that may exist at the gene expression level but that cannot be explained by standard outpatient assessment.

Table 1 Demographic, clinical and serological characteristics of lcSSc patients and healthy controls

Overview of gene expression profiles

There were 218 differentially expressed genes (DEG) between lcSSc and healthy control samples (false discovery rate (FDR) <10%). Of these 218 genes, 181 (83%) were upregulated in lcSSc samples and were significantly enriched in several terms related to the extracellular matrix (ECM), e.g. ECM organisation and ECM component and response to growth factor, tissue development and regulation of the serine/threonine kinase signalling pathway. This gene signature included genes previously implicated in the pathogenesis of SSc such as ANGPT2, CD163, COMP, CTGF and TIMP2, among others (Fig. 1a). The full list of significantly differentially expressed genes is available in Additional file 1.

Fig. 1

Gene expression in whole skin samples from patients with limited cutaneous systemic sclerosis (lcSSc) and healthy controls (Controls). a Heatmap: red indicates upregulation, green indicates downregulation. Gene expression in the limited 1 subgroup of samples from patients with lcSSc differs markedly from both the limited 2 subgroup and from healthy controls, which shared similar gene expression patterns. Sample genes of interest are listed. b Hierarchical clustering of genes distinguishing lcSSc samples from healthy control samples identifies three distinct clusters, termed limited 1, limited 2 and healthy control. Norm Normal. c Pairwise comparisons between limited 1, limited 2 and controls. Data are plotted as mean with SEM values. P values were derived from the Kruskal-Wallis test followed by Dunn's multiple comparisons test

Even though we specifically asked for genes differentially expressed between lcSSc samples and healthy controls, 6 lcSSc samples (5 unique patients) clustered with controls, whereas the remaining 13 lcSSc samples (10 unique patients) formed a distinct cluster (Fig. 1b). We designated 13 lcSSc samples as the “limited 1” group and 6 lcSSc samples as the “limited 2” group. The expression of 181 genes significantly upregulated in lcSSc vs. healthy controls was also significantly different in all pairwise comparisons between limited 1, limited 2 and control samples (Fig. 1c). Mean ± standard error of the mean expression values for 181 gene signatures were as follows: 0.23 ± 0.01 for limited 1, -0.02 ± 0.03 for limited 2 and -0.46 ± 0.03 for healthy control samples.

Since comparison between lcSSc and controls was suggestive of the presence of two lcSSc groups, we performed multiclass SAM to identify DEG between limited-1, limited-2 and healthy control samples: 807 genes were differentially expressed between these three groups (FDR <10%). Again, limited-2 samples clustered with controls and separately from the limited-1 group (Fig. 2; see Additional file 2 for the full gene list).

Fig. 2

Multiclass significance analysis of microarrays (SAM) between the three identified clusters. Multiclass SAM was used to identify 807 differentially expressed genes (false discovery rate <10%) between limited -1, limited -2 and healthy control samples. Sample genes are listed. The majority of differentially expressed genes (475/807, 58.9%) had increased expression in the limited-1 group

Genes with significantly increased expression in either the limited-2 group or in healthy controls were not functionally enriched. Genes with significantly increased expression in samples from the limited-1 group (475/807, 58.9%) displayed very strong enrichment in functional terms related to ECM and vasculature development (p < 10-10) and strong enrichment in cell adhesion, response to growth factor, mitogen-activated protein kinase (MAPK) cascade and response to wounding (p < 10-5), among others. Enrichment in immune signalling was comparatively weak (e.g. p = 0.0183 for immune system process and p = 0.0208 for inflammatory response). The entire g:Profiler output is listed in Additional file 3.

For the limited-1 group, the term with the most significant functional enrichment was cardiovascular system development (p = 3.01 × 10-14): 70 out of 475 genes with increased expression in limited-1 samples were involved in this biological process. This set of 70 genes essentially recaptured many of the biological themes observed in the bigger limited-1 gene expression signature. For example, there was very strong enrichment in angiogenesis, response to growth factor and cell proliferation (p < 10-10) and strong enrichment in ECM, cartilage development, MAPK cascade, cell adhesion and response to wounding (p < 10-5), among others. It was also strongly enriched in response to transforming growth factor (TGF)-β (p = 2.87 × 10-7, 12/70 genes). Additional file 4 contains the entire g:Profiler output for this analysis. There were no significant differences in demographic, clinical or serological parameters that correlated with limited-1 and limited-2 groups, in particular, the presence of significant microvascular involvement and hallmark antibody reactivities were similar across each group as shown in Table 2 below.

Table 2 Clinical and laboratory parameters of the limited-1 and limited-2 subgroups

Validation of microarray gene expression profiles in an independent lcSSc cohort

We applied the same approach described above to validate our findings in an independent cohort of patients with lcSSc derived from an American population (the Boston University cohort), with samples comprising 24 lcSSc and 4 control samples. Similar to this study cohort, analysis of differential expression between lcSSc and control samples suggested the presence of two subgroups of patients with lcSSc (group 1 with 14 lcSSc samples and group 2 with 10 lcSSc samples). We then compared the gene expression of the CVS development gene signature across subgroups of lcSSc patients and controls from the validation study. This gene expression signature was significantly different among both subgroups of patients with LcSSc and a control group (Fig. 3b), paralleling the results from this study (Fig. 3a).

Fig. 3

Cardiovascular development trends in discovery and validation cohorts. Pairwise comparisons between limited 1, limited 2 and control samples in the discovery cohort (a) and an independent validation cohort (b). Data are plotted as mean with SEM values. P values were derived from the Kruskal-Wallis test followed by Dunn's multiple comparisons test. CVS cardiovascular system


In this study we have assessed molecular heterogeneity in gene expression profiles from skin biopsies taken from uninvolved skin of lcSSc patients compared with matched healthy controls. We have identified two subgroups within the lcSSc patient population in this study, which we termed limited 1 and limited 2. Patients from the limited-2 subgroup correspond to a third of the lcSSc cases and are generally characterised by a subtle alteration in gene expression that resembles but is distinct from the control samples. The majority of lcSSc cases cluster into a limited-1 subgroup associated with a substantial number of genes with significant differential expression, which are involved in multiple functional terms with known or potential relevance to SSc pathogenesis. We have confirmed our findings in a second cohort.

Most compelling is the identification of a cardiovascular system (CVS) development gene expression signature. This is interesting because of the well-recognised vascular abnormalities that are seen in lcSSc [3]. The idea that these pathways may be upregulated or altered in lcSSc even at sites that are not clinically affected is notable and may point towards an inherent susceptibility to vasculopathy that may be a hallmark of SSc. It is plausible that using uninvolved lcSSc skin we were able to minimise the impact of pathways reflecting inflammation and that this may explain the novel findings in our study. It was reassuring to find that several genes identified within this group were well-recognised in scleroderma pathogenesis, including, for instance, COMP, THBS1 and CTGF. These have traditionally been considered as pro-fibrotic markers, but all have a role in the regulation of vascular function [23,24,25]. More traditional “vascular” markers such as MMP19 and COL4A1 were also included. It is also attractive to consider that recapitulation of developmental gene expression programmes may be central to susceptibility to development of SSc since it can be envisaged as a disease of perturbed connective tissue repair. There is a strong precedent for the same genetic pathways and programmes that are used in embryonic and post-natal growth and development to be recapitulated in acquired disease [26]. This may also fit with the model of SSc being a susceptibility genotype that we have recently postulated for the key complication of pulmonary arterial hypertension [27].

It is also notable that there were no clear clinical or serological associations with the molecular subgroups of lcSSc. This is relevant since it suggests that gene-expression-based subsets may indeed add to the clinical and serological factors that are already used in the clinic to subgroup patients. This would be analogous to the intrinsic subsets of diffuse cutaneous SSc (dcSSc) that have been reported and that are already of clinical utility [7]. Of note, in the cohort from which the study samples were taken, 40% of patients with limited disease were anti-topoisomerase (Scl-70)-positive and while this is traditionally thought to represent diffuse scleroderma, the patients in this cohort are categorised according to the distribution of skin disease at an early stage of disease and retain that subset over the course of their follow-up period. There were, therefore, no patients included in this study who would have experienced remission to a low skin score with previous diffuse or more severe skin disease.

It is striking that many of the genes involved in CVS development are also seen in other relevant pathways such as ECM and response to TGF-β. Not only does this give some key candidate genes and pathways that could be perturbed in SSc but it also is a reminder that these factors that have been focussed on for a role in fibrosis may also have other roles in the biology or the disease. Thus it is plausible that the patterns of gene expression are altered in SSc and that this may have distinct consequences depending upon the stage and subset of disease. Further work will be required to assess whether these differences in expression reflect susceptibility to SSc that depends on genetic or epigenetic factors or whether they are reflective of pathogenetic processes that are occurring within the clinically uninvolved skin in limited SSc.

Strengths of the study include the well-characterised patients from a single centre, which may reduce the variability that can be a hallmark of multicentre cohorts and ensures standard classification and clinical evaluation. In addition, the use of well-established gene expression analysis platforms and semi-automated sample processing minimises technical variation that has been a limitation in some studies. Inclusion of a matched control group of healthy individuals was important to allow reliable interpretation of the “limited” gene expression signature.

Weaknesses of the study are the limitations of a relatively small study cohort that may not allow generalisability of our findings. However, the results are in line with those previously reported for subjects with lcSSc included in earlier gene expression studies [7]. We have validated our results in an independent cohort of patients with consistent findings across all groups. There is also a risk, given the small size of the study that there may be confounding effects on gene expression due to differences in immunosuppressive therapies, as an example. There were no statistically significant differences or strong trends between the limited-1 and limited-2 groups in terms of demographics, clinical history, treatments or serology in our study.

We have not included patients with clinically involved forearm skin in this study. There are no adequately powered studies that confirm concordance between clinically involved and uninvolved skin but previous data do suggest this; these confirmatory studies are required to extend and validate our findings further as part of a larger future study of gene expression in lcSSc. Involved distal skin is not biopsied in our patients due to concerns about wound healing. Finally, both further post-transcriptomic functional studies, and modern techniques that allow structural analysis of mRNA expression most notably to examine the dermal microcirculation and perivascular space would localise gene expression changes and verify that differential gene expression, particularly from the CVS development cluster, are reflected by differential protein expression.


We showed that gene expression profiling of biopsies from uninvolved skin in lcSSc differentiates two potential subgroups that overlap with other clinical and serological features. This may reflect important differences in pathogenesis within these patient groups. In addition, we identified differential expression of a subset of genes that relate to CVS development. Since the lcSSc subset is characterised by vasculopathy in the skin and internal organs, this may shed light on underlying mechanisms of vascular disease in SSc. The clinical implications of our findings will need to be analysed in future larger studies.



Antinuclear antibody




Cardiovascular system


Diffuse cutaneous systemic sclerosis


Differentially expressed genes


Extracellular matrix


extractable nuclear antigen


European League Against Rheumatism


False discovery rate


Healthy controls


Indirect immunofluorescence


Limited cutaneous systemic sclerosis


Mitogen-activated protein kinase


Modified Rodnan skin score




Significance analysis of microarrays


Systemic sclerosis


Transforming growth factor


  1. 1.

    Denton C and Moinzadeh P. Systemic sclerosis. In: Oxford Textbook of Rheumatology. Fourth Edition. Edited by Richard A. Watts, Philip G. Conaghan, Christopher Denton, Helen Foster, John Isaacs and Ulf Müller-Ladner. Oxford: Oxford University Press; 2013.

  2. 2.

    Varga J. Overview: Pathogenesis Integrated. In: John V, Denton CP, Wigley FM, editors. Scleroderma, from Pathogenesis to Comprehensive Management. New York: Springer; 2012. p. 163–4.

  3. 3.

    Nihtyanova SI, Schreiber BE, Ong VH, Rosenberg D, Moinzadeh P, Coghlan JG, Wells AU, Denton CP. Prediction of pulmonary complications and long term survival in systemic sclerosis. Arthritis Rheum. 2014;66(6):1625–35.

  4. 4.

    Herrick A, Lunt M, Whidby N, Ennis H, Silman A, McHugh N, Denton C. Observational study of treatment outcome in early diffuse cutaneous systemic sclerosis. J Rheumatol. 2010;37:116–24.

  5. 5.

    Villalta D, Imbastaro T, Di Giovanni S, Lauriti C, Gabini M, Turi MC, Bizzaro N. Diagnostic accuracy and predictive value of extended autoantibody profile in systemic sclerosis. Autoimmun Rev. 2012;12:114–20.

  6. 6.

    Chakravarty EF, Martyanov V, Fiorentino D, Wood TA, Haddon DJ, Jarrell JA, Utz PJ, Genovese MC, Whitfield ML, Chung L. Gene expression changes reflect clinical response in a placebo-controlled randomized trial of abatacept in patients with diffuse cutaneous systemic sclerosis. Arthritis Res Ther. 2015;17:159.

  7. 7.

    Milano A, Pendergrass SA, Sargent JL, George LK, McCalmont TH, Connolly MK, Whitfield ML. Molecular subsets in the gene expression signatures of scleroderma skin. PLoS One. 2008;3:e2696.

  8. 8.

    Hinchcliff M, Huang CC, Wood TA, Matthew Mahoney J, Martyanov V, Bhattacharyya S, Tamaki Z, Lee J, Carns M, Podlusky S, Sirajuddin A, Shah SJ, Chang RW, Lafyatis R, Varga J, Whitfield ML. Molecular signatures in skin associated with clinical improvement during mycophenolate treatment in systemic sclerosis. J Invest Dermato. 2013;133:1979–89.

  9. 9.

    Gardner H, Shearstone JR, Bandaru R, Crowell T, Lynes M, Trojanowska M, Pannu J, Smith E, Jablonska S, Blaszczyk M, Tan FK, Mayes MD. Gene profiling of scleroderma skin reveals robust signatures of disease that are imperfectly reflected in the transcript profiles of explanted fibroblasts. Arthritis Rheum. 2006;54:1961–73.

  10. 10.

    Rice LM, Stifano G, Ziemek J, Lafyatis R. Local skin gene expression reflects both local and systemic skin disease in patients with systemic sclerosis. Rheumatology (Oxford). 2016;55(2):377–9.

  11. 11.

    Rice LM, Ziemek J, Stratton EA, McLaughlin SR, Padilla CM, Mathes AL, Christmann RB, Stifano G, Browning JL, Whitfield ML, Spiera RF, Gordon JK, Simms RW, Zhang Y, Lafyatis R. A longitudinal biomarker for the extent of skin disease in patients with diffuse cutaneous systemic sclerosis. Arthritis Rheumatol. 2015;67(11):3004–15.

  12. 12.

    Assassi S, Swindell WR, Wu M, Tan FD, Khanna D, Furst DE, Tashkin DP, Jahan-Tigh RR, Mayes MD, Gudjonsson JE, Chang JT. Dissecting the heterogeneity of skin gene expression patterns in systemic sclerosis. Arthritis Rheumatol. 2015;67(11):3016–26.

  13. 13.

    Farina G, Lafyatis D, Lemaire R, Lafyatis R. A four-gene biomarker predicts skin disease in patients with diffuse cutaneous systemic sclerosis. Arthritis Rheum. 2010;62:580–8.

  14. 14.

    Johnson ME, Mahoney JM, Taroni J, Sargent JL, Marmarelis E, Wu MR, Varga J, Hinchcliff ME, Whitfield ML. Experimentally-derived fibroblast gene signatures identify molecular pathways associated with distinct subsets of systemic sclerosis patients in three independent cohorts. PLoS One. 2015;10(1):e0114017.

  15. 15.

    van den Hoogen F, Khanna D, Fransen J, Johnson SR, Baron M, Tyndall A, Matucci-Cerinic M, Naden RP, Medsger Jr TA, Carreira PE, Riemekasten G, Clements PJ, Denton CP, Distler O, Allanore Y, Furst DE, Gabrielli A, Mayes MD, van Laar JM, Seibold JR, Czirjak L, Steen VD, Inanc M, Kowal-Bielecka O, Müller-Ladner U, Valentini G, Veale DJ, Vonk MC, Walker UA, Chung L, Collier DH, Csuka ME, Fessler BJ, Guiducci S, Herrick A, Hsu VM, Jimenez S, Kahaleh B, Merkel PA, Sierakowski S, Silver RM, Simms RW, Varga J, Pope JE. 2013 classification criteria for systemic sclerosis: an American College of Rheumatology/European League against Rheumatism collaborative initiative. Arthritis Rheum. 2013;65:2737–47.

  16. 16.

    Reich M, Liefeld T, Gould J, Lerner J, Tamayo P, Mesirov JP. GenePattern 2.0. Nat Genet. 2006;38(5):500–1.

  17. 17.

    Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007;8(1):118–27.

  18. 18.

    Reese SE, Archer KJ, Therneau TM, Atkinson EJ, Vachon CM, de Andrade M, Kocher JP, Eckel-Passow JE. A new statistic for identifying batch effects in high-throughput genomic data that uses guided principal component analysis. Bioinformatics. 2013;29(22):2877–83.

  19. 19.

    de Hoon MJ, Imoto S, Nolan J, Miyano S. Open source clustering software. Bioinformatics. 2004;20(9):1453–4.

  20. 20.

    Tusher VG, Tibshirani R, Chu G. Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001;98(9):5116–21.

  21. 21.

    Saldanha AJ. Java Treeview − extensible visualization of microarray data. Bioinformatics. 2004;20(17):3246–8.

  22. 22.

    Reimand J, Kull M, Peterson H, Hansen J, Vilo J. g:Profiler–a web-based toolset for functional profiling of gene lists from large-scale experiments. Nucleic Acids Res. 2007;35(Web Server issue):W193–200.

  23. 23.

    Fu Y, Kong W. Cartilage Oligomeric Matrix Protein: Matricellular and Matricrine Signaling in Cardiovascular Homeostasis and Disease. Curr Vasc Pharmacol. 2017 [epub ahead of print].

  24. 24.

    Lawler PR, Lawler J. Molecular basis for the regulation of angiogenesis by thrombospondin-1 and -2. Cold Spring Harb Perspect Med. 2012;2:a006627.

  25. 25.

    Sonnylal S, Shi-Wen X, Leoni P, Naff K, Van Pelt CS, Nakamura H, Leask A, Abraham D, Bou-Gharios G, de Crombrugghe B. Selective expression of connective tissue growth factor in fibroblasts in vivo promotes systemic tissue fibrosis. Arthritis Rheum. 2010;62(5):1523–32.

  26. 26.

    Rayner KJ, Liu PP. Long Noncoding RNAs in the Heart: The regulatory roadmap of cardiovascular development and disease. Circ Cardiovasc Genet. 2016;9(2):101–3.

  27. 27.

    Gilbane AJ, Derrett-Smith E, Trinder SL, Good RB, Pearce A, Denton CP, Holmes AM. Impaired bone morphogenetic protein receptor II signaling in a transforming growth factor-β-dependent mouse model of pulmonary hypertension and in systemic sclerosis. Am J Respir Crit Care Med. 2015;191(6):665–77.

Download references


Not applicable.


This work was supported by a grant from EULAR Orphan Diseases Programme. The funding body was not involved in design of the study and collection, analysis or interpretation of data or writing the manuscript. In addition, our research activity benefits from support from grants awarded by Arthritis Research UK, Scleroderma & Raynaud’s UK and the British Heart Foundation. PM’s work was supported by a grant from Koeln Fortune 224/2011, 131/2013 and 155/2014.

Availability of data and materials

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

Authors’ contributions

CD and MW conceived and designed the study; ED-S, CCh, PM, KK and CCa obtained the data; ED-S, VM, CCh, PM and TW analyzed the data. ED-S, VM, VO, DA, PLM, RL, MW and CD interpreted the data. ED-S, VM, MW and CD drafted the manuscript and CCh, PM, KK, TW, CCa, VO, DA, PLM and RL critically appraised it. All authors read and approved the final manuscript and are accountable.

Authors’ information

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

All patients included in this study signed informed consent for their clinical and laboratory data to be used in this clinical research project. Ethics approval was obtained from the local ethics committee: NRES Committee London-Hampstead HRA reference 6398.

Publisher’s Note

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

Author information



Corresponding author

Correspondence to Christopher P. Denton.

Additional files

Additional file 1:

Genes significantly differentially expressed between lcSSc and control groups. Green cells: genes with increased expression in controls. Yellow cells: genes with increased expression in lcSSc samples. (XLSX 16 kb)

Additional file 2:

Genes significantly differentially expressed between the limited-1, limited-2 and control groups. Orange cells: genes with increased expression in limited 1 group. Green cells: genes with increased expression in controls. Grey cells: genes with increased expression in limited 2 group. (XLSX 38 kb)

Additional file 3:

g:Profiler output for 475 genes with increased expression in the limited-1 group (multiclass SAM, FDR <10%). Additional files 3 contain the following column headers. Descriptions are as follows: p value: significance of enrichment in a given term corrected for multiple hypothesis testing using default g:GOSt method g:SCS. Q&T: overlap between genes in the query (Q) and genes in the genome annotated to a given term (T). t type: term category - Gene Ontology: Biological Process (BP), Molecular Function (MF), Cellular Component (CC); KEGG (keg); Reactome (rea). t name: term name. Q&T list: list of genes forming the overlap between query (Q) and genome lists for a given term (T). (XLSX 25 kb)

Additional file 4:

g:Profiler output for 70 genes with increased expression in the limited-1 group annotated to cardiovascular system development. Additional file 4 contains the following column headers. Descriptions are as follows: p value: significance of enrichment in a given term corrected for multiple hypothesis testing using default g:GOSt method g:SCS. Q&T: overlap between genes in the query (Q) and genes in the genome annotated to a given term (T). t type: term category - Gene Ontology: Biological Process (BP), Molecular Function (MF), Cellular Component (CC); KEGG (keg); Reactome (rea). t name: term name. Q&T list: list of genes forming the overlap between query (Q) and genome lists for a given term (T). (XLSX 29 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Derrett-Smith, E.C., Martyanov, V., Chighizola, C.B. et al. Limited cutaneous systemic sclerosis skin demonstrates distinct molecular subsets separated by a cardiovascular development gene expression signature. Arthritis Res Ther 19, 156 (2017).

Download citation


  • Scleroderma
  • Genetics
  • Systemic sclerosis
  • Vasculopathy
  • Microarray