A novel multi-network approach reveals tissue-specific cellular modulators of fibrosis in systemic sclerosis
© The Author(s). 2017
Received: 5 August 2016
Accepted: 23 February 2017
Published: 23 March 2017
Systemic sclerosis (SSc) is a multi-organ autoimmune disease characterized by skin fibrosis. Internal organ involvement is heterogeneous. It is unknown whether disease mechanisms are common across all involved affected tissues or if each manifestation has a distinct underlying pathology.
We used consensus clustering to compare gene expression profiles of biopsies from four SSc-affected tissues (skin, lung, esophagus, and peripheral blood) from patients with SSc, and the related conditions pulmonary fibrosis (PF) and pulmonary arterial hypertension, and derived a consensus disease-associate signature across all tissues. We used this signature to query tissue-specific functional genomic networks. We performed novel network analyses to contrast the skin and lung microenvironments and to assess the functional role of the inflammatory and fibrotic genes in each organ. Lastly, we tested the expression of macrophage activation state-associated gene sets for enrichment in skin and lung using a Wilcoxon rank sum test.
We identified a common pathogenic gene expression signature—an immune–fibrotic axis—indicative of pro-fibrotic macrophages (MØs) in multiple tissues (skin, lung, esophagus, and peripheral blood mononuclear cells) affected by SSc. While the co-expression of these genes is common to all tissues, the functional consequences of this upregulation differ by organ. We used this disease-associated signature to query tissue-specific functional genomic networks to identify common and tissue-specific pathologies of SSc and related conditions. In contrast to skin, in the lung-specific functional network we identify a distinct lung-resident MØ signature associated with lipid stimulation and alternative activation. In keeping with our network results, we find distinct MØ alternative activation transcriptional programs in SSc-associated PF lung and in the skin of patients with an “inflammatory” SSc gene expression signature.
Our results suggest that the innate immune system is central to SSc disease processes but that subtle distinctions exist between tissues. Our approach provides a framework for examining molecular signatures of disease in fibrosis and autoimmune diseases and for leveraging publicly available data to understand common and tissue-specific disease processes in complex human diseases.
KeywordsSystemic sclerosis Scleroderma Macrophage Lung disease Functional genomics
Integrative genomics has yielded powerful tissue-specific functional networks that model the interaction of genes in these specialized “microenvironments” . These tools hold promise for understanding how genes may contribute to human diseases  that arise, in part, out of an aberrant interplay of cell types and tissues. Network biology has played a crucial role in our understanding of complex human diseases such as cancer [3, 4] and, more recently, in disorders where the interactions among multiple tissues are dysregulated . Analytical approaches that leverage biological “big data” can be especially fruitful in rare and heterogeneous diseases , for which the risk of mortality is significant and no approved treatments exist. We performed an integrative, multi-tissue analysis for systemic sclerosis (SSc; scleroderma), a disease for which all of these tenets are true, and included samples from patients with pulmonary fibrosis (PF) and pulmonary arterial hypertension (PAH).
SSc is a systemic disease characterized by abnormal vasculature, adaptive immune dysfunction (autoantibody production), and extracellular matrix (ECM) deposition in skin and internal organs. The etiology of SSc is unknown, but it has complex genetic risk  and postulated triggers include immune activation by cancer , infection , or dysbiosis . SSc is clinically heterogeneous, with some patients experiencing rapidly progressive skin and internal organ disease while others have stable disease that is largely limited to skin. Understanding the molecular processes in multiple affected organ systems is critical to understanding the pathogenesis of SSc and other complications, such as PF and PAH, that co-occur in these patients. Here, we ask if deregulated pathways are distinct or common between these tissues affected by SSc and if each organ manifestation has distinct disease signatures at the molecular level.
An integrative genomics study of SSc is of particular importance. Gene expression data from multiple tissues, including skin [11–13], lung [14, 15], and esophagus (ESO) , now exist. However, the rarity of the disease results in studies with small sample sizes and the multi-organ nature makes it difficult to assess molecular changes across organ systems relative to controls. Therefore, analyzing the data from multiple tissues that are more difficult to obtain (e.g., esophagus and lung) in the context of tissues that are more easily assayed (e.g., skin and peripheral blood) is a powerful way to make inferences about pathogenesis in internal organs. In addition, putting SSc disease-specific findings in the context of tools built from biological big data is a way to bolster and refine our findings for this rare disease.
We previously developed mutual information consensus clustering (MICC) to identify gene expression that is conserved across multiple, disparate datasets . Here we expanded MICC to perform an integrative, multi-tissue analysis of SSc and related fibrotic conditions. We included gene expression datasets from ten different cohorts representing four different affected tissues from patients with SSc. Following MICC, we used the Genome-scale Integrated Analysis of gene Networks in Tissues (GIANT) tissue-specific functional genomic networks  to identify gene–gene interactions among those expressed consistently across affected tissues. These big data approaches integrate individual experiments measuring hundreds of disease states and biological perturbations. Integration of these data holds promise for understanding how genes contribute to organ-specific manifestations of human diseases . These GIANT networks are a detailed, genome-scale representation of the functional interactions between genes in different microenvironments.
We identified a pathogenic signature—a common “immune–fibrotic axis”—that is present in all tissues analyzed and is increased in the most severe disease complications, including PF and PAH. Using tissue-specific functional networks , we analyzed the nature of the immune–fibrotic axis to understand the gene–gene interactions that underlie fibrosis across organ systems. Using differential network analysis, we were able to identify skin- and lung-specific gene–gene interactions relevant to macrophage (MØ) plasticity and SSc pathophysiology. We now propose a model that implicates alternatively activated MØs as part of the immune–fibrotic axis that may drive fibrosis in multiple tissues.
Patients and datasets
Datasets included in this study
Phenotypes of interest
Inflammatory subset, proliferative subset
Milano et al. 
Inflammatory subset, proliferative subset
Pendergrass et al. 
Inflammatory subset, proliferative subset
Hinchcliff et al. 
Mahoney et al. 
Christmann et al. 
SSc-PF, IPF, IPAH, SSc-PAH
Hsu et al. 
Inflammatory subset, proliferative subset, SSc-PAH
Taroni et al. 
Pendergrass et al. 
Risbano et al. 
We used the patient disease label (e.g., PAH) as published in the original work for all of these sets. Below, we note some important characteristics (for the purposes of this work) of the included patient populations. As noted in the “Results” section, the two lung datasets contained patients with different histological patterns of lung disease. Some patients included in the PBMC dataset, including those with PAH, also had interstitial lung disease, though exclusion of these patients does not significantly change the interpretation as put forth in Pendergrass et al. . As illustrated in Additional file 1, two datasets (ESO, LSSc) did not contain healthy control samples and three datasets (UCL, LSSc, and PBMC) were comprised entirely of LSSc patients.
Microarray dataset processing
This work contains ten datasets on multiple microarray platforms. Agilent datasets (Pendergrass, PBMC, Milano, Hinchcliff, ESO, UCL, LSSc) used either Agilent Whole Human Genome (4x44K) Microarrays (G4112F) (Pendergrass, PBMC, Milano, Hinchcliff, ESO, UCL) or 8x60K (LSSc). Data were Log2-transformed and lowess normalized and filtered for probes with intensity twofold over local background in Cy3 or Cy5 channels. Data were multiplied by −1 to convert to Log2(Cy3/Cy5) ratios. Probes with >20% missing data were excluded. The Illumina dataset (Bostwick, HumanRef-8 v3.0 BeadChips) was processed using variance-stabilizing transformation xand robust spline normalization using the lumi R package. Dr. Christmann provided the raw data in the form of.CEL files. Dr. Feghali-Bostwick provided Illumina BeadSummary files. Affymetrix datasets (Risbano, HGU133plus2; Christmann, HGU133A_2) were processed using the Robust Multiarray Averaging (RMA) method as implemented in the affy R package. Batch bias was detected in the ESO dataset. To adjust these data, missing values were imputed via k-nearest neighbor algorithm using a GenePattern  module with default parameters and the data were adjusted using ComBat  run as a GenePattern module to eliminate the batch effect.
To compare datasets in our downstream analysis, duplicate genes must not be present in the dataset and must be summarized in some way. First, we annotated each probe with its Entrez gene ID. Agilent 4x44K arrays were annotated using the hgug4112a.db Bioconductor package. LSSc was annotated using UNC Microarray Database with annotations from the manufacturer. Probes annotated to lincRNAs (A19) were removed from the analysis. The Illumina dataset was annotated by converting the gene symbols (provided as part of the BeadSummary file) to Entrez IDs using the org.Hs.eg.db package. The Risbano PBMC dataset was annotated using the hgu133plus2.db package. The Christmann dataset was annotated using an annotation file from the manufacturer. Probes that did not map to any Entrez ID and probes that mapped to multiple Entrez IDs were removed in all cases. Probes that mapped to the same Entrez ID were collapsed to the gene mean using the aggregate function in R, followed by gene median centering.
Clustering of microarray data and statistical tests for phenotype association
The collapsed datasets were used to find coherent coexpression modules. We used Weighted Gene Co-expression Network Analysis (WGCNA), a strong clustering method, which allows us to automatically detect the number of coexpression modules and remove outliers . Each dataset was clustered using the blockwiseModules function in WGCNA R package using the signed network option and power = 12; all other parameters were set to default. WGCNA does not identify large, densely connected coexpression modules in random data  and although changing the soft-thresholding power ultimately changes the resulting modules, we and others find the resulting modules to be stable and concordant across parameter choices .
Using the WGCNA coexpression modules also reduces the dimensionality of the dataset, as it allows us to test for genes’ association with, or differential expression in, a particular pathophenotype of interest on the order of tens, rather than thousands, using the module eigengene. The module eigengene is the first principal component and represents the expression of all genes in a module and an idealized hub of the coexpression module. We used the moduleEigengenes function in the WGCNA R package to extract the eigengenes. A module was considered to be pathophenotype-associated if the module eigengene was significantly differentially expressed in or significantly correlated with a pathophenotype of interest. Only two-class categorical variables were considered using a Mann–Whitney U test (i.e., all pulmonary fibrosis and pulmonary arterial hypertension patients were grouped together regardless of underlying etiology). We used Spearman correlation for continuous values. P values were Bonferroni-corrected on a per-phenotype basis. See Additional files 2, 3, 4, 5, 6, 7, 8, and 9 for complete output of these analyses. In the main text, we discuss categorical pathophenotypes, as these were enriched at the consensus cluster level. We do find instances of coexpression modules that are associated with continuous pathophenotypes, such as pulmonary function test measurements, but these were not apparent at the consensus cluster level of abstraction.
Module overlap network construction and community detection
where N is the total number of genes shared between the two datasets. The overlap scores can be positive, negative, or zero, indicating that the modules overlap more, less, or the same as expected at random, respectively. As shown in Mahoney et al. , the overlap scores can be naturally thresholded using information theory to yield a sparse network of significant overlaps—the module overlap network. We performed a permutation test to test the significance of the mutual information between a pair of partitions (datasets) and found that the true value of the mutual information of partitions was higher than all sampled values of the null distribution (permuted p = 0; see Additional file 10 for permutation test details and Additional file 11 for the results of this test). This is consistent with mutual information being implicitly computed relative to a null model.
The module overlap network is highly structured. For example, a module representing an inflammatory process in skin often significantly overlaps inflammatory modules in other tissues. Thus, the structure of the module overlap network corresponds to the biological processes that are common to multiple datasets. We can identify these processes by clustering the module overlap network itself. Community detection is a procedure used to identify clusters in networks. The type(s) of community detection we employed is based on the concept of modularity (see also Additional file 12; Glossary of terms used in this paper). Networks with high modularity have sets of nodes (here, coexpression modules) that are more densely connected within a set and more sparsely connected outside of that set . Community detection methods based on this concept take into account the expected amount of edges within a set of nodes and detect the sets of nodes that are more densely connected than expected (communities) .
We used two methods of community detection. First, we used fast-greedy modularity maximization (implemented in Matlab) , which yielded large, diffuse communities. (The fast-greedy modularity optimization algorithm has a known bias for the size of communities it selects and is thought to find “low-resolution” clusters in some cases .) We call these “top-level” communities. We tested whether the modularity of the module overlap network as calculated by fast-greedy community detection was significant relative to a network where module labels are randomly permuted. This allowed us to assess whether the module overlap graph had significant community structure; the results were highly significant (permuted p = 0; see Additional file 10 for permutation test details).
Because the above algorithm is greedy, it only finds a local maximum for modularity. To find smaller, more densely connected sub-communities, we used spin-glass community detection (igraph R package implementation, max number of communities = 10, all other parameters were set to default) [26, 27]. This algorithm implements a stochastic algorithm to maximize the modularity function resulting in tighter clusters than the fast-greedy algorithm . We call these “bottom-level” communities. The community/sub-community structure of the module overlap network demonstrates that there is a hierarchy of biological processes that are common across datasets, where large communities contain smaller ones. To display this hierarchical community structure, we first sorted by top-level community label, and then within each community we sorted by bottom-level label. The adjacency matrix of the module overlap network and its node attributes (including dataset of origin and community labels) are supplied in Additional files 13 and 14, respectively.
Bonferroni-corrected p values, Fisher’s exact test pathophenotype-associated modules in top-level communities in the module overlap graph
“In SSc” p value
“In inflammatory” p value
“In proliferative” p value
“In PAH” p value
“In PF” p value
Functional and pathophenotype annotation of the module overlap network
We used biological pathways from the Kyoto Encyclopedia of Genes and Genomes (KEGG) , BioCarta, and Reactome  obtained from Molecular Signatures Database from the Broad Institute (http://software.broadinstitute.org/gsea/msigdb). The Jaccard similarity between the edge and pathway will be equal to one if all of the genes shared between two modules are exactly the same set of genes annotated to the pathway, or zero if no genes are shared between the two sets. To functionally annotate a community in the information graph, we compared the Jaccard similarities of the edges within the community to edges outside of the community using a Mann–Whitney U test (with Bonferroni adjustment). The full results of this analysis are included as Additional files 15, 16, 17, 18, 19, 20, 21, 22, and 23.
Tissue consensus gene sets
As each tissue was considered separately (limited skin and diffuse skin were considered separately), five tissue consensus gene sets were generated; the union of these tissue consensus datasets was used to query the functional genomic networks and is referred to as the “immune–fibrotic axis consensus” gene set or genes throughout the text. For all genes in modules in clusters 4A and 4B, we calculated the Pearson correlation to their respective module eigengene. We compared this correlation of consensus genes to that of non-consensus genes using a Mann–Whitney U test. Additional file 24 contains the tissue consensus genes from 4AB or the immune–fibrotic axis consensus genes.
Querying GIANT functional networks, single tissue network analysis, and network visualization
The GIANT functional genomic networks were obtained as binary (.dab) files and processed using the Sleipnir library for computational functional genomics . We queried all networks (lung, skin, “all tissue”, macrophage) using the immune–fibrotic axis consensus gene sets (as Entrez IDs) and pruned all low probability (<0.5) edges. All networks are available for download from the GIANT webserver (http://giant.princeton.edu/) . For the single tissue analysis (e.g., lung network), we considered only the largest connected component of each network and performed spin-glass community detection as implemented in the igraph R package  to obtain the functional modules. We annotated functional modules using g:Profiler  using all genes in a module as a query. All networks in this work were visualized using Gephi . The network layout was determined by community membership, the strength of connections between communities, and finally the interactions between individual genes. The lung network node attribute file and edge lists are supplied as Additional files 25 and 26.
Differential network analysis
where Alung, Askin, and Aglobal are the lung, skin, and global (all tissues) adjacency matrices from GIANT. The differential lung network is thus the lung network minus the maximum edge weight from the skin and lung networks, where all edges that are stronger in skin or the global network are set to zero. Thus, the differential lung network contains only highly lung-specific interactions. Functional modules in the lung differential network were found using spin-glass community detection (see "Querying GIANT functional networks, single tissue network analysis, and network visualization") within the largest connected component of the network. The differential lung network node attributes and edge list are supplied as Additional files 27 and 28.
To perform the macrophage-specific network analysis in the supplemental material, we subtracted global edge weights from the macrophage network, setting negative edges to zero (as above). We then permuted the order of the adjacency matrix (edges) 1000 times and assessed if the true weight within a community was more than random (red), less than random (blue), or no different from random (white). We performed the same permutation testing on the lung network with global subtraction and found more weight than expected “on-diagonal” and less weight than expected “off-diagonal”; this demonstrates how spin-glass community detection takes into account the expected distribution of edges.
Differential expression and MØ gene set analysis
where x gi is the gene expression value in array/sample i, μ g is the gene mean, and σ g is the gene standard deviation. The average z-score of genes in a set (module from Xue et al. ). computed for an array/sample to summarize gene set expression. Mann–Whitney U tests were used to compare average z-scores between groups. To validate these findings in an independent SSc skin dataset, we used the data from Assassi et al.  as processed by the authors and deposited in NCBI Gene Expression Omnibus (GEO; GSE58095 series matrix). We collapsed duplicated probes to the gene mean, z-scored genes as in Eq. 6, and compared average z-scores as above.
We performed an integrative analysis of ten independent gene expression datasets containing samples from patients with SSc and associated co-morbidities (Table 1). The primary goal of this study was to identify the fundamental processes that occur across end-target and peripheral tissues of patients with SSc and related fibrotic conditions. Secondly, we aimed to identify the presence or absence of common gene expression patterns that underlie the molecular intrinsic subsets of SSc  in different organs. Analysis of multiple tissue biopsies from patients with skin fibrosis, esophageal dysfunction, PF, and PAH allowed us to determine in a data-driven manner whether these tissues were perturbed in a similar manner on a genomic scale.
Number of microarrays and WGCNA coexpression modules in each of the datasets included in this study
Number of arrays
Number of coexpression modules
Severe pathophenotypes share a common immune–fibrotic axis
The module overlap network is agnostic to the clinical phenotypes corresponding to each biopsy. To associate communities in the module overlap network with SSc and fibrotic pathophenotypes, we tested each of the 549 modules for differential expression in relevant pathophenotypes (see “Methods”). For example, every lung module in the PAH cohorts was tested for differential expression in PAH. Clusters 4A and 4B in the module overlap network contain modules with increased expression in all pathophenotypes of interest: the inflammatory and proliferative subsets of SSc, PAH, and PF (Fig. 2b). Thus, the modules in these communities correspond to a common, broad disease signal that is present in every pathophenotype under study. As with our prior studies, we did not find a strong association with autoantibody subtype and the co-expression modules identified here.
Selected pathways that are similar to overlapping coexpression patterns in consensus clusters in the information graph
Summary of selected pathways
Cell–cell junction organization
Metabolism of proteins
T cytotoxic and helper pathway
Antigen processing and presentation
ECM receptor interaction
Signaling by PDGF
G2 M checkpoint
Unwinding of DNA
Nuclear receptors in lipid metabolism and toxicity
Fatty acid metabolism
PPAR signaling pathway
FGFR ligand binding and activation
Edges in 4B encode pro-fibrotic processes, including ECM receptor interaction, collagen formation, and TGF-β signaling (Table 4). Cluster 4B consists of skin inflammatory and fibroproliferative subset-associated modules as well as lung PAH-, late PF- and early PF-associated modules (Fig. 2b; Additional file 30: Figure S1). These results replicate and expand what we have found in our prior meta-analysis of skin data alone : the expression patterns observed in the SSc intrinsic subsets are shared with other tissues and SSc-associated pathophenotypes and indicative of altered immune and fibrotic processes (an immune–fibrotic axis).
To understand how the immune–fibrotic axis and these phenotypes are functionally related, we identified the consensus genes in the combined 4A and 4B clusters (see “Methods”; 2079 unique genes; Additional file 24). Consensus genes are highly central within their respective dataset gene–gene correlation networks and our procedure identifies sets of genes that capture disease-specific variation (Additional file 30: Figure S2). Using a conservative measure, these consensus genes are enriched for genes with increased expression in all disease manifestations (SAM , FDR <5%; PF in both lung datasets p < 2.2e-16; PAH lung, p = 7.88 × 10−5; PAH in both PBMC datasets, p = 3.20 × 10−15, Fisher’s exact test). This demonstrates that the tissue consensus genes are highly relevant to all disease manifestations in this study. The tissue consensus gene sets allow us to rigorously extrapolate from this conservative set a substantially broader, disease-associated signal. This extrapolation is especially important for tissue studies that are underpowered to detect a large number of significantly differentially expressed genes (see “Discussion”). We took the union of the tissue consensus gene sets as a set of “immune–fibrotic axis consensus genes” that are informative about pathology in every tissue.
The lung functional genomic network reveals a coupling of immune and fibrotic processes
Our analysis includes two lung datasets derived from both early SSc-PF (open lung biopsies obtained for diagnostic purposes ) and end-stage or late disease (SSc-PF patients that underwent lung transplantation ). In addition to the differences in disease stage between these two datasets, there is also some difference in the histological patterns of fibrosis in these cohorts. In the Bostwick lung dataset , all patients with SSc-PF had usual interstitial pneumonia (UIP). This study used lung tissues from patients who underwent lung transplantation (late disease). The Christmann lung dataset  contains five patients with non-specific interstitial pneumonia (NSIP) and two patients with centrilobular fibrosis. This study looked at early SSc-PF patients, used open lung biopsies, and specifically avoided honeycombing areas.
Although NSIP and UIP have distinct clinical outcomes, they have been shown to be nearly indistinguishable at the gene expression level . Furthermore, these datasets have overlapping coexpression patterns as demonstrated by their shared community membership in the module overlap network. Comparison of different datasets allows us to determine how genes with increased expression at these different stages and histological subtypes of lung disease are distributed throughout the lung network. Genes overexpressed in SSc-PF (SAM, PF versus Normal comparison, FDR <5%) are distributed throughout the lung network and therefore are predicted to participate in all of the molecular processes identified in the network. Quantification of the distribution of SSc-PF differentially expressed genes throughout the consensus lung network (Fig. 3b) demonstrates that molecular processes can be associated with either a disease stage, histopathological pattern, or both stages/patterns. The cell cycle module contains only early/NSIP SSc-PF genes, the innate immune response/NF-κB/apoptotic processes module contains more late/UIP SSc-PF genes, and the response to TGF-β module contains genes from both histological patterns (Fig. 3a, b).
Hub and bridge genes are associated with the pathogenesis of pulmonary fibrosis
Selected genes in the consensus lung network
Function/potential role in disease
BUB3 mitotic checkpoint protein
Encodes a mitotic cell cycle checkpoint protein that regulates the onset of anaphase
Cell division cycle 7
Regulates MCM complex
Minichromosome maintenance complex component 3
Subunit of minichromosome maintenance (MCM) complex
MutS homolog 6
Participates in DNA mismatch repair.
ECM disassembly/wound healing
CD44 molecule (Indian blood group)
A hyaluronic acid receptor that can interact with many other ligands found in the ECM. Primary idiopathic PF fibroblasts exhibit an invasive phenotype that was abrogated with treatment with anti-CD44 
Has been observed to interact with TIMP1 
Regulates NPC2 secretion, TNF-alpha production, and cholesterol trafficking genes in an animal model of obesity 
Regulates NPC2 secretion, TNF-alpha production, and cholesterol trafficking genes in an animal model of obesity 
Galactosidase, beta 1
Mutations in this gene can lead to GM1-gangliosidosis, a manifestation of which includes foam cell accumulation in the lungs 
Niemann-Pick disease, type C2
Mutations in this gene result in a lipid storage disorder. Functions in the regulation of cholesterol trafficking through the lysosome by binding to cholesterol released from low density lipoproteins taken up by cells
Transforming growth factor, beta-induced
Induced by phagocytosis of apoptotic debris in monocyte-derived MØs and regulates collagen turnover 
TIMP metallopeptidase inhibitor 1
Has been observed to interact with CD63 and overexpression has been noted to inhibit apoptosis in a CD63-dependent manner 
Innate immune response/NFkB signaling/apoptotic process
Baculoviral IAP repeat-containing protein 3
Cysteine-rich, angiogenic inducer, 61
Dual specificity phosphatase 6
Plays a role in the positive regulation of apoptosis 
Fas cell surface death receptor
Cell surface death receptor
Nuclear factor of kappa light polypeptide gene enhancer in B-cells inhibitor, epsilon
Negative regulator of NFkB signaling
Plasminogen activator, urokinase receptor
Also known as uPAR. Contains an SSc risk SNP. Pulmonary fibroblasts from patients with idiopathic PF over express uPAR and that uPAR ligation results in a hypermotile phenotype 
Phospholipid scramblase 1
Regulates phospholipid membrane asymmetry
Tumor necrosis factor, alpha-induced protein 3
Also known as A20. Contains an SSc risk SNP (also associated with other autoimmune conditions). Negative regulator of NFkB signaling
Tumor necrosis factor (ligand) superfamily, member 10
Also known as TRAIL. Elevated in serum of SSc patients 
Tumor necrosis factor receptor superfamily, member 10b
Also known as TRAILR2
Major histocompatibility complex, class I, E
Class I MHC molecule
Major histocompatibility complex, class I, F
Class I MHC molecule
IFN induced transmembrane protein 1
SSc-PF (UIP and NSIP)
IFN induced transmembrane protein 2
IFN induced transmembrane protein 3
IFN regulatory factor 1
Activator of type I IFN signaling
2′-5′-Oligoadenylate synthetase 1, 40/46 kDa
Involved in innate immune response to viral infection
Response to TGF-beta
Contains an SSc risk SNP
Connective tissue growth factor
Dab, mitogen-responsive phosphoprotein, homolog 2 (Drosophila)
SSc-PF (UIP and NSIP)
Extracellular matrix protein.
Laminin gamma1 chain
Expression of this gene is essential for the development of basement membranes 
Mediates cell-to-cell and cell-to-matrix interactions. Putative biomarker of modified Rodnan skin score 
The innate immune response/NF-κB signaling/apoptotic process module contains genes that are highly expressed in late/UIP SSc-PF, including the hub genes CYR61 and TM4SF1 (Fig. 3a, b; Additional file 30: Figure S3). The hub gene TNFAIP3 (A20), which is increased in late SSc-PF (Fig. 3e), is a negative regulator of NF-κB signaling and inhibitor of TNF-mediated apoptosis. The innate immune response/NF-κB signaling/apoptotic process and IFN/antigen presentation modules are bridged by TNFSF10, also known as TRAIL (TNF-related apoptosis inducing ligand; Fig. 3f). These results suggest that the balance of apoptosis is altered in late/UIP SSc-PF. The upregulation of genes with anti-apoptotic function was not reported in the original study , which demonstrates the strength of both the MICC method and the study of functional interactions.
CD44 and PLAUR (uPAR) bridge multiple functional modules in the lung network (Fig. 3f) and have been implicated in IPF [39, 40]. Because these genes link modules important in regulating disease progression, therapeutic targeting of CD44 and uPAR may be an effective strategy in combating SSc-PF. Indeed, anti-CD44 treatment reduces fibroblast invasion and bleomycin-induced lung fibrosis , and inhibition of uPAR ligation significantly reduces motility of pulmonary fibroblasts from patients with idiopathic PF . These results are consistent with our identification of these genes as key genes in the lung network.
The lung microenvironment provides a distinct milieu for pro-fibrotic processes
These highly specific interactions are displayed in Fig. 4b, where a cell is red if it is lung-specific or blue if it is skin-specific (cf. Additional file 30: Figure S4). The number of tissue-specific edges in each functional module is quantified in Fig. 4b, c, which illustrate that most functional modules in lung have fewer interactions than in skin, with the exception of the cell cycle module. Of particular interest is the relationship between the phagolysosome/ECM disassembly genes and response to TGF-β genes, as strong differential connectivity can be observed in this module (Fig. 4b, c). Thus, even though ECM disassembly and TGF-β module genes are coordinately differentially expressed in both lung and skin, they are differentially connected to each other, suggesting that the microenvironment strongly determines the functional consequences of upregulating these pro-fibrotic genes.
To summarize lung-specific biological processes in the immune–fibrotic axis, we clustered the lung-specific interactions (differential lung network) to identify lung-specific pathways (Additional file 30: Figure S5). We identified 23 clusters corresponding to biological processes such as type I IFN signaling (cluster 10), antigen processing and presentation (cluster 4), REACTOME cell surface interactions at the vascular wall (cluster 22), and mitotic cell cycle (cluster 16; shown in Additional file 30: Figure S5b). Taken together, this suggests that within the immune–fibrotic axis we find innate immune and cell proliferation processes that are highly lung-specific. One of the largest of these clusters (cluster 13; Fig. 4d; Additional file 30: Figure S5c) includes NPC2, S100A4, and CTSB, which encode protein products that are highly expressed in normal lung-resident MØs (LR-MØs) [41, 42].
NPC2 is a hub of the ECM disassembly/wound healing module in the full lung network (Fig. 3); many of the genes in cluster 13 also belong to the ECM disassembly/wound healing module in the whole network, including the cathepsins CTSB and CTSL. Alveolar MØs are the main source of cathepsins in bleomycin-induced fibrotic lung tissue . Additional genes associated with development and maintenance of alternative MØ activation include TGFBI , NEU1 , PRCP , and DAB2 . Genes that are specifically associated with alternative activation of lung MØs include PLP2  and IFITM1  (Fig. 4d; Additional file 30: Figure S5c). Based on these genes and the complete lung network in Fig. 3, we identified an LR-MØ signature. These findings are consistent with previous reports of alternative MØ activation in SSc [14, 50].
The LR-MØ cluster in the differential lung network also contains a number of genes implicated in lipid storage disorders, including HEXB, GLB1, and NPC2. Several other LR-MØ cluster genes have been shown to be important for regulating cholesterol trafficking genes in an animal model of obesity, including CTSB, CTSL, and NPC2 . It has been noted that lipid metabolism genes are upregulated in lung MØs relative to other tissue-specific MØs . Furthermore, in the bleomycin injury mouse model of pulmonary fibrosis, lipid-laden MØs have been observed to increase expression of markers associated with alternative MØ activation and to secrete TGF-β .
Distinct MØ gene expression programs are elevated in lung and skin
We hypothesized that early SSc-PF lung samples may have evidence of both alternatively activated and lipid-stimulated MØs and that this may differ from what is observed in skin. The presence of alternatively activated MØs in the inflammatory subset of skin was inferred in our single tissue analysis . To test this hypothesis, we used gene sets associated with classic activation of MØs, alternative activation of MØs, or stimulation of MØs with a variety of activation stimuli, including free fatty acids, taken from Xue et al. . To summarize the expression of each MØ gene set  and compare across tissues in these data, we computed the average expression of all genes in each gene set (see “Methods”; see Additional file 31 for a mapping between Xue et al. modules and our naming scheme). Results are displayed for control and SSc-PF lung, as well as control and SSc-inflammatory skin (Fig. 5b). As shown in Fig. 5b, there is evidence of an increase in alternatively activated and free fatty acid stimulated gene sets in SSc-PF and SSc-inflammatory skin. These data do not show statistically significant differences in expression of gene sets associated with classic MØ activation between controls and SSc-PF or SSc-inflammatory skin (see Additional file 32 for p values of all modules tested).
The discovery of IFN-related genes among the consensus genes indicates that these pathways are increased in pathophenotypes of interest (e.g., SSc-PF and the skin inflammatory subset). Christmann et al. also noted a strong IFN-related gene signature in SSc-PF samples, although the cellular compartment responsible for this signature was not described . Because stimulation with IFN results in classic activation of MØs, we examined the expression of genes from CL 1, as it is most strongly associated with IFN-γ treatment (“classic activation”) in human MØs . However, CL 1 genes’ expression is not different between disease and controls in either skin or lung (Wilcoxon p = 0.76 and 0.80, respectively; Fig. 5b). This result is consistent with our inability to discern differences in classic MØ activation markers between controls and SSc-PF and inflammatory skin and suggests that classically activated MØs are not the source of the reported IFN signature we find.
Modules ALT 1 and ALT 2 are both associated with IL-4 and IL-13 treatment, which are stimuli associated with alternative activation of MØs . These two gene sets are non-overlapping coexpression modules and therefore represent two “parts” of the alternatively activated MØ transcriptional program. We performed functional enrichment analysis for ALT 1 and 2 to understand which biological processes underlie these transcriptional signatures (see “Methods”). Module ALT 1 is enriched for genes involved in oxidative phosphorylation (KEGG, p < 0.0001) and the citric acid cycle (REACTOME, p < 0.0001) pathways. In lung, ALT 1 expression is higher in SSc-PF than in controls (Wilcoxon p = 0.0046). There is no difference between healthy controls and the inflammatory subset in skin (Wilcoxon p = 0.41). Module ALT 2 shows an opposite trend and is enriched for genes implicated in the positive regulation of response to wounding (Gene Ontology (GO) biological process (BP), p = 0.027) and defense response (GO BP, p = 0.00035); this module includes alternatively activated MØ markers such as CD14 and CCL26 [53, 54]. ALT 2 expression is increased in the inflammatory subset in skin (Wilcoxon p = 0.041) and trends toward decreased expression in SSc-PF lung (Wilcoxon p = 0.16). Together, these pathways suggest a metabolic “switch” associated with alternative activation in lung that is not found in skin (for review see ; Fig. 5b).
We also analyzed modules associated with free fatty acid (FFA) stimulation, which are relevant to the question of lipid signaling or exposure in SSc tissues (FFA 1, 2, and 3). We first performed functional enrichment analysis for these modules to gain biological insight into these transcriptional programs. FFA 1 is enriched for genes involved in the unfolded protein response (REACTOME, p = 0.025). FFA 2 is enriched for antigen processing-cross presentation genes (REACTOME; p = 0.00101). FFA 3 is enriched for genes in the ER-phagosome pathway (REACTOME, p = 0.0076). Expression of FFA 1 and 2 is significantly increased in lung (FFA 1 Wilcoxon p = 0.046, p = 0.97 in skin; FFA 2 Wilcoxon p = 0.0013, p = 0.63 in skin), whereas FFA 3 is upregulated in SSc-PF lung (Wilcoxon p = 0.0013) and the SSc inflammatory subset in skin (Wilcoxon p = 0.00056). These results suggest that LR- MØs may have a distinct lipid exposure that strongly diverges from that in skin.
We repeated this analysis in an independent SSc skin dataset (Assassi et al. ) to validate our findings (Additional file 30: Figure S7). Assassi et al. reported that macrophage transcripts are elevated in SSc skin but used a “general” to macrophages signature gene list that does not provide information about activation state. The results from Assassi et al. largely agree with the results from the Hinchcliff dataset: ALT 2, but not ALT1 (Wilcoxon p = 0.0682), is significantly increased in SSc skin (Wilcoxon p = 5.92e-05), and FFA 3 is the only FFA module significantly increased in SSc skin (Wilcoxon p = 3.219e-06; FFA 1 p = 0.928; FFA 2 p = 0.486). The only disparity between the two skin datasets is that we find that CL 1 is significantly increased in the Assassi SSc patients (p = 0.000856). This difference may be because we looked at all SSc patients rather than “fibroinflammatory” patients alone, or due to the increased coverage of the genome on the platform used. Overall, analysis of Assassi et al. supports the differences in MØ alternative activation programs and lipid response in SSc-affected skin and lung.
The differential network analysis (Fig. 4) allowed us to identify highly lung-specific interactions in the immune–fibrotic axis that implicated lipid signaling as a distinct functional process in lung. The higher expression of multiple free fatty acid-associated modules in lung suggests that the role of lipid signaling in MØs may be more important in this tissue than in skin, consistent with what we would predict based on highly lung-specific gene–gene interactions, and based on prior biomedical literature in related conditions [48, 52]. Thus, a major difference between the lung and skin networks can be attributed to the presence of a distinct MØ phenotype in lungs.
SSc is a systemic disease that affects multiple internal organs. Herein, we present the first study of molecular mechanism of disease across multiple affected organ systems in SSc. To our knowledge, we show for the first time that a common set of cell types and pathways are driving disease across these affected organs, and importantly that it can also be found in related fibrotic conditions.
Gene expression data have been collected for multiple tissues in SSc and related conditions. However, these data often have issues that are common to many rare diseases. First, SSc is not prevalent and patients with particular disease manifestations are still rarer, so there is a limit to the amount of biopsy material available for study. Second, for practical and ethical reasons, internal organ biopsies are seldom taken from healthy subjects, making comparisons difficult. Thus, lung, esophagus, and other affected internal organs are more difficult to study than blood and skin tissue. Therefore, there is a critical need to leverage our biological prior knowledge with our understanding of well-studied tissues—like blood and skin—to make plausible inferences about pathogenesis in tissues that are more difficult to study.
The clinical heterogeneity of SSc, particularly the difficulty of predicting internal organ involvement, raises an important question: are the fibrotic processes observed in multiple organs derived from a common disease process, or is each organ manifestation effectively a distinct disease? Our analyses demonstrate that there is a common gene expression signature underlying all severe organ manifestations of SSc—the immune–fibrotic axis—in solid organs. The immune–fibrotic axis underlies both SSc pulmonary manifestations of PF and PAH, and the intrinsic subsets of skin and esophagus. Moreover, coexpression modules from peripheral blood, a mixture of innate and adaptive immune cells, have significant overlap with modules associated with all pathophenotypes studied. Thus, while fibrotic processes were largely associated with solid tissues, the inflammatory component of the immune–fibrotic axis is only found in peripheral blood.
The presence of a common gene expression signature across multiple tissues suggests a common disease driver, but it does not resolve the possible tissue-specific processes that contribute to disease in the internal organs. Indeed, there are many layers of biological regulation between gene expression and whole tissue phenotypes. Resolving the relationship between molecular profiles and phenotypes is a difficult biological problem underlying most biomedical inquiry. However, these relationships have been approximated by integrating high-throughput genomic data into tissue-specific functional networks using big data machine-learning strategies . We addressed tissue specificity in SSc pathology by interpreting the common expression signal—the immune–fibrotic axis—within these tissue-specific functional networks. These networks allowed us to identify critical genes that occupy important positions in molecular pathways in lung. It is clear from this work that the coupling of immune and fibrotic processes is a hallmark of SSc that occurs in SSc-PF and SSc-PAH as well as skin. However, we also find subtle, lung-specific functional differences that we attribute, in part, to the plasticity of the myeloid cell lineage.
The plasticity of the myeloid lineage may drive tissue-specific SSc disease processes
Altered immune function has been implicated in the pathogenesis of SSc [56, 57]. In most prior studies, characterization of macrophage activation has relied on analysis of a very limited number of surface markers and/or a few characteristic mRNAs [56, 57]. Most of these studies have concluded that SSc macrophages bear an M2 activation profile based on CD163 and/or CD206 expression. Macrophage polarization spans a broad spectrum of activation states, ranging from “classically activated” or M1 cells, which largely mediate pro-inflammatory responses to “alternatively activated” or M2 cells, which are predominantly associated with immune suppression and wound healing. While expression of CD206 and CD163 is higher in alternatively activated macrophages compared with “classically activated” macrophages, it is difficult to make global generalizations about macrophage activation based on such limited analysis. While operationally useful, the designation of M1 versus M2 activation has limited utility in vivo as macrophage activation is informed by the local cytokine milieu to which these cells are exposed.
Our study of multiple skin cohorts showed that multiple gene expression markers of activated MØs are elevated in SSc skin across multiple data sets, consistent with gene expression profiling of lung tissue from SSc patients with interstitial lung disease . These data are consistent with elevated levels of IL-4 and IL-13 in SSc sera [58, 59]. Furthermore, CD68+ MØs have been identified as producers of IL-13 in human SSc skin biopsies and genetic deficiency of IL-13 is protective against disease in a mouse model of SSc . IL-13 activates tissue fibrosis  and genetic and observational studies link IL-13 with SSc pathogenesis [62–64].We have further demonstrated that SSc MØs express high levels of profibrotic cytokines, suggesting they play a significant role in mediating fibrosis and in maintaining an inflammatory environment in SSc (unpublished data).
By performing a combined analysis of SSc gene expression in multiple tissues, we are able to observe and infer, in a genome-wide manner, commonalities in the complex mixture of cell types in a tissue at the time of biopsy. Overwhelmingly, we detected a MØ signature associated with severe disease. In the module overlap network, we find that PAH-associated modules from PBMCs [19, 20] have significant overlap with SSc inflammatory subset-associated modules from skin and esophagus (Fig. 2). Indeed, in Pendergrass et al. , we observed that PBMCs from lcSSc patients have significant enrichment in myeloid- and MØ-related gene sets compared to healthy controls. Christmann et al.  expanded on this, showing that highly expressed transcripts in LSSc-PAH CD14+ monocytes were induced in IL-13-stimulated cells, i.e., that PAH monocytes are alternatively activated. We assert that this MØ polarization is a significant part of the immune–fibrotic axis we find in these data and, therefore, is likely a common driver of the complex pathophysiology of SSc. In support of this, an independent study also identified MØs and dendritic cells (DCs) as possible sources of an “inflammatory” signature in lesional SSc skin .
We found evidence for the contribution of LR-MØs to SSc-PF pathobiology, consistent with the alternative activation of MØs and TGF-β production. In our prior analysis of skin, we inferred alternatively activated MØs as modulators of the SSc inflammatory intrinsic subset in skin . Our current study identifies a LR-MØ signature within the functional relationships of immune–fibrotic axis consensus genes in lung (Figs. 4d and 5a). We posit that the differences in fibrotic responses of skin and lung tissue are due, in large part, to innate differences between tissue-resident MØs that have been observed [66, 67], as well as the interactions between infiltrating monocytes and tissue-resident cell types (e.g., alveolar epithelial cells versus keratinocytes). Because MØ phenotype and function are plastic and readily modulated by the local tissue microenvironment, it is likely that differential activation of MØs in these tissues is the result of exposure to distinct cytokine milieu. Indeed, we show that distinct alternative activation gene expression programs have increased expression in SSc-PF lung and inflammatory SSc skin (Fig. 5). In particular, there were multiple lipid-related signatures elevated in SSc-PF lung alone.
We cannot rule out that the MØ changes we observe are a secondary response to the affected organ pathology. Regardless, therapies that target MØ effectors such as IL6R have shown promise in clinical trials  and MØ chemoattractants have been shown to be important in animal models of SSc inflammatory disease, suggesting that MØs play a central role in SSc pathogenesis. We also cannot rule out that DCs contribute to our results, as plasmacytoid DCs are observed to be important in the stiff skin syndrome mouse model . However, some skin-resident DCs have been shown to be transcriptionally similar to peripheral blood monocytes in humans . We speculate that the circulation of peripheral myeloid cells contributes to the multi-organ nature of SSc. Future studies may use in silico and cell-sorting techniques to deconvolve SSc expression data to identify changes in cell proportion and transcriptome throughout the disease course and to finely phenotype myeloid cells from SSc patient tissue samples.
Summary of SSc-PF disease processes
We found that gene signatures that are increased in alternatively activated human MØs and MØs treated with free fatty acids are enriched in early SSc-PF patients and that there is no evidence for enrichment of a pro-inflammatory, IFN-stimulated MØ signature (Fig. 5) .
Christmann et al. had previously identified an increase in IFN- and TGF-β-regulated genes in biopsies from early SSc-PF , but it was unclear which cell types were responsible for the IFN signature or if there was evidence of distinct subpopulations of MØs. Increased CCL18 protein and higher CD163 mRNA were observed in lungs of patients with SSc-associated interstitial lung disease, suggestive of the presence of alternatively activated MØs .
We also find elevated gene expression programs associated with MØ alternative activation (specifically metabolic “reprogramming”) and lipid exposure in Christman et al. (Fig. 5).The LR-MØ signature identified in our differential network analysis consisted of genes with increased expression in early SSc-PF that participate in lipid and cholesterol trafficking (Fig. 4d; Additional file 30: Figure S5). The expression of these genes is correlated with “canonical” MØ genes identified in  (Fig. 5). In the bleomycin injury mouse model of pulmonary fibrosis, lipid-laden MØs, or foam cells, have been observed to upregulate markers associated with alternative MØ activation and to secrete TGF-β . Oxidized phospholipid treatment also causes alternative activation and TGF-β secretion in human MØs . Consistent with this report, recent work demonstrates that foam cell formation in vivo favors the development of a pro-fibrotic MØ activation profile [71, 72]. These studies, along with our results, suggest that lipid exposure or uptake in MØs may be important.
We find genes from both datasets in the response to the TGF-β module of the lung network. TGF-β signaling is a hallmark of SSc and other fibrotic diseases, and was noted in the initial analysis of both SSc lung datasets [14, 15]. However, we also find evidence that the type I IFN signature is present in the Bostwick dataset (Fig. 3). The functional module most strongly associated with late stage disease/UIP is the innate immune, NF-κB, and apoptotic processes module. This module is connected to the TGF-β module through components of the fibrinolysis pathway such as PAI-1 (SERPINE1; Fig. 3). PAI-1 is upregulated in late stage SSc-PF and is known to be important in pulmonary fibrosis [73–75]. One mechanism by which fibrinolysis may contribute to the resolution of fibrosis is through the induction of fibroblast apoptosis . Both TGF-β1 and PAI-1 have been shown to inhibit lung fibroblast apoptosis .
We found evidence for a shift in the balance of apoptosis in the Bostwick dataset, perhaps in myofibroblasts , in our network analyses (Fig. 6). Long-lived myofibroblasts are thought to continually deposit collagen and contribute to persistent fibrosis . This apoptotic-resistance phenotype is related to the stiffness of the matrix , suggesting that a shift in apoptotic processes may occur once the deposition of excess collagen begins. Moreover, impaired phagocytosis of apoptotic cells, or efferocytosis, has been observed in the alveolar MØs of IPF patients . We find genes involved in efferocytosis, specifically in receptors (CD44) and endocytic machinery associated with this process, in the lung network (Figs. 3 and 6) . If the shift in apoptosis and efferocytosis occurs, we speculate that the fibrotic and inflammatory processes in our network will also be altered. Efferocytosis by alveolar MØs plays a key role in the resolution of inflammation in the lung through the subsequent release of TGF-β . We hypothesize that, following initial injury, TGF-β signaling, antifibrinolytic factors, and the disruption of apoptosis and efferocytosis may contribute to progressive fibrosis in SSc-PF (Fig. 6).
Limitations and future directions
A limitation of this study is a lack of post-genomic validation, particularly in lung. This work is in essence hypothesis-generating, but the need for this study is highlighted by the sparseness of biopsy material, and it provides new directions for inquiry into the pathogenesis of the disease.
Our results suggest that alternatively activated MØs likely play a central role in the pathogenesis of SSc by activating fibroblasts. Most importantly, they show for the first time that this is likely to occur across multiple affected organ systems in SSc patients. Future experiments will need to examine these cells functionally to determine if SSc MØs can activate other cell types (e.g., fibroblasts) to produce ECM and to examine the role of these cells in mouse models of fibrosis as well as gene expression in multiple organs from the same patient. Our integrative genomics approach directly compares multiple tissues and manifestations and suggests that there may be subtle differences in the MØ phenotype in SSc-affected skin and lung. This supports the fine phenotyping of these cells from SSc patient tissue samples when possible, and the possibility of targeting these cells therapeutically.
In this study, we have utilized data from multiple tissues to examine the systemic nature of SSc. Our integrative analysis allowed us to leverage well-studied tissues to inform us about SSc manifestations that are under-studied molecularly. This study rigorously tests the notion that patients with severe disease have shared immunological and fibrotic alterations. The common immune–fibrotic axis shows evidence for alternatively activated MØs in multiple SSc tissues. However, there are subtle differences in the MØ gene expression programs detected in skin and lung. Different microenvironments likely provide distinct stimuli to infiltrating MØs that determine the pro-fibrotic character of these cells. The plasticity of this lineage is likely central to the divergence of fibrotic processes in multiple SSc-affected tissues and is a central component of an immune–fibrotic axis driving disease.
False discovery rate
Free fatty acid
Gene Expression Omnibus
Genome-scale Integrated Analysis of gene Networks in Tissues
Idiopathic pulmonary arterial hypertension
Idiopathic pulmonary fibrosis
Lung resident macrophage
Limited cutaneous systemic sclerosis
Mutual information consensus clustering
Non-specific interstitial pneumonia
Pulmonary arterial hypertension
Peripheral blood mononuclear cell
Significance Analysis of Microarrays
scleroderma-associated pulmonary arterial hypertension
scleroderma-associated pulmonary fibrosis
University College London
Usual interstitial pneumonia
Weighted Gene Co-expression Network Analysis
JNT would like to thank members of the Whitfield Laboratory for thoughtful discussion.
This work has been supported by grants from: the Scleroderma Research Foundation (SRF http://www.srfcure.org/) to MLW and MH; grants from the Scleroderma Foundation to PAP and MLW; the Dr. Ralph and Marian Falk Medical Research Trust Catalyst Award to MLW; Dartmouth SYNERGY award to MLW; the Gordon and Betty Moore Foundation (GBMF4552) to CSG; National Institutes of Health (NIH; https://www.nih.gov/) NIAMS grants 2R01AR051089 to RL, P50AR060780 and P30AR061271 to RL and MLW, and K23AR059763 to MH. JNT received support from NIH-NIGMS T32GM008704 and the John H. Copenhaver, Jr. and William H. Thomas, MD 1952 Junior Fellowship from Dartmouth Graduate Studies. CPD received support from the EULAR ODP programme. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
Previously published gene expression datasets used in this study were deposited in the Gene Expression Omnibus repository (http://www.ncbi.nlm.nih.gov/geo/) under the following accession numbers: GSE9285, GSE32413, GSE45485, GSE59785, GSE76808, GSE48149, GSE68698, GSE19617, and GSE22356. The expression datasets that were previously unpublished are available from the Gene Expression Omnibus under the following accession numbers: GSE76806 (LSSc) and GSE76807 (UCL). The functional genomic networks used in this study are publicly available at http://giant.princeton.edu/. The GIANT sub-networks analyzed as part of this study are included in this published article (and its additional files). See Zenodo (https://zenodo.org/record/216614) and figshare (http://doi.org/10.6084/m9.figshare.4490798) for additional expression data files and analysis scripts and intermediate files, respectively.
JNT, JMM, and MLW conceived of the study. JNT, CSG, VM, JMM, and MLW designed data analyses, performed analyses, and interpreted the results. TAW performed the microarray experiments. RBC, HWF, RAL, and CPD designed study cohorts included in this work and contributed samples and/or data. MEH provided clinical expertise and interpreted the results. PAP provided macrophage biology expertise and interpreted the results. JNT, PAP, JMM, and MLW wrote the paper. All authors read, revised, and approved the manuscript.
CPD has been a consultant to Roche, GlaxoSmithKline, Actelion, Inventiva, CSL Behring, Takeda, Merck-Serono, MedImmune, and Biogen. MLW and MH have filed patents for gene expression biomarkers in systemic sclerosis. MLW is a scientific founder of Celdara Medical LLC. MLW has served as consultant to GlaxoSmithKline, Bristol Myers Squib, EMD Serono, Biogen, and Quintiles. RL has received both grants and consulting fees from Genzyme/Sanofi, Shire, Regeneron, Biogen, BMS, Inception, Precision Dermatology, PRISM, UCB, Pfizer, and Roche/Genentech; he received consulting fees from Lycera, Novartis, Celgene, Amira, Celdara, Celltex, Dart Therapeutics, Idera, Intermune, Medimmune, Promedior, Zwitter, Actelion, EMD Serono, Akros, Extera, Reneo, Scholar Rock, and HGS. The remaining authors have no competing interests.
Consent for publication
Ethics approval and consent to participate
The LSSc and UCL studies are previously unpublished. The samples from the LSSc dataset were obtained at Boston University Medical Center (BUMC)/Boston Medical Center (BMC); the BUMC/BMC Institution Review Board approved this study. The samples from the UCL dataset were obtained at University College of London; the London-Hampstead NRES Committee approved this study. The Dartmouth College CPHS approved this work. All subjects gave informed consent. All research conformed to the principles expressed in the Declaration of Helsinki.
Open AccessThis 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.
- Greene CS, Krishnan A, Wong AK, Ricciotti E, Zelaya RA, Himmelstein DS, Zhang R, Hartmann BM, Zaslavsky E, Sealfon SC, Chasman DI, FitzGerald GA, Dolinski K, Grosser T, and Troyanskaya OG. Understanding multicellular function and disease with human tissue-specific networks. Nat Genet. 2015;47:569–576.
- Gross AM, Ideker T. Molecular networks in context. Nat Biotechnol. 2015;33(7):720–1.View ArticlePubMedGoogle Scholar
- Hofree M, Shen JP, Carter H, Gross A, Ideker T. Network-based stratification of tumor mutations. Nat Methods. 2013;10(11):1108–15.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen S, Wang Q, Wu Z, Li Y, Li P, Sun F, Zheng W, Wu Q, Wu C, Deng C, Zhang F, Li Y. Genetic association study of TNFAIP3, IFIH1, IRF5 polymorphisms with polymyositis/dermatomyositis in Chinese Han population. PLoS One. 2014;9(10), e110044.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen JC, Cerise JE, Jabbari A, Clynes R, Christiano AM. Master regulators of infiltrate recruitment in autoimmune disease identified through network-based molecular deconvolution. Cell Syst. 2015;1(5):326–37.View ArticlePubMedPubMed CentralGoogle Scholar
- Sharma A, Menche J, Huang C, Ort T, Zhou X, Kitsak M, Sahni N, Thibault D, Voung L, and Guo F. A disease module in the interactome explains disease heterogeneity, drug response and captures novel pathways and genes. Hum Mol Genet. 2015;24(11):3005–3020.
- Assassi S, Radstake TR, Mayes MD, Martin J. Genetics of scleroderma: implications for personalized medicine? BMC Med. 2013;11:9.View ArticlePubMedPubMed CentralGoogle Scholar
- Joseph CG, Darrah E, Shah AA, Skora AD, Casciola-Rosen LA, Wigley FM, Boin F, Fava A, Thoburn C, Kinde I, Jiao Y, Papadopoulos N, Kinzler KW, Vogelstein B, Rosen A. Association of the autoimmune disease scleroderma with an immunologic response to cancer. Science. 2014;343(6167):152–7.View ArticlePubMedGoogle Scholar
- Farina A, Cirone M, York M, Lenna S, Padilla C, McLaughlin S, Faggioni A, Lafyatis R, Trojanowska M, Farina GA. Epstein-Barr virus infection induces aberrant TLR activation pathway and fibroblast-myofibroblast conversion in scleroderma. J Invest Dermatol. 2014;134(4):954–64.View ArticlePubMedGoogle Scholar
- Arron ST, Dimon MT, Li Z, Johnson ME, A Wood T, Feeney L, G Angeles R, Whitfield ML, Angeles J. High Rhodotorula sequences in skin transcriptome of patients with diffuse systemic sclerosis. J Invest Dermatol. 2014;134(8):2138–45.View ArticlePubMedPubMed CentralGoogle Scholar
- 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(7):e2696.View ArticlePubMedPubMed CentralGoogle Scholar
- Pendergrass SA, Lemaire R, Francis IP, Mahoney JM, Lafyatis R, Whitfield ML. Intrinsic gene expression subsets of diffuse cutaneous systemic sclerosis are stable in serial skin biopsies. J Invest Dermatol. 2012;132(5):1363–73.View ArticlePubMedPubMed CentralGoogle Scholar
- 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 Dermatol. 2013;133(8):1979–89.View ArticlePubMedPubMed CentralGoogle Scholar
- Christmann RB, Sampaio-Barros P, Stifano G, Borges CL, de Carvalho CR, Kairalla R, Parra ER, Spira A, Simms R, Capellozzi VL, Lafyatis R. Association of Interferon- and transforming growth factor β-regulated genes and macrophage activation with systemic sclerosis-related progressive lung fibrosis. Arthritis Rheumatol. 2014;66(3):714–25.View ArticlePubMedPubMed CentralGoogle Scholar
- Hsu E, Shi H, Jordan RM, Lyons-Weiler J, Pilewski JM, Feghali-Bostwick CA. Lung tissues in patients with systemic sclerosis have gene expression patterns unique to pulmonary fibrosis and pulmonary hypertension. Arthritis Rheum. 2011;63(3):783–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Taroni JN, Martyanov V, Huang C, Mahoney JM, Hirano I, Shetuni B, Yang G, Brenner D, Jung B, Wood TA, Bhattacharyya S, Almagor O, Lee J, Sirajuddin A, Varga J, Chang RW, Whitfield ML, and Hinchcliff M. Molecular characterization of systemic sclerosis esophageal pathology identifies inflammatory and proliferative signatures. Arthritis Res Ther. 2015;17(1):194.
- Mahoney JM, Taroni J, Martyanov V, Wood TA, Greene CS, Pioli PA, Hinchcliff ME, Whitfield ML. 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Preliminary criteria for the classification of systemic sclerosis (scleroderma). Subcommittee for scleroderma criteria of the American Rheumatism Association Diagnostic and Therapeutic Criteria Committee. Arthritis Rheum. 1980;23(5):581–90.
- Pendergrass SA, Hayes E, Farina G, Lemaire R, Farber HW, Whitfield ML, Lafyatis R. Limited systemic sclerosis patients with pulmonary arterial hypertension show biomarkers of inflammation and vascular injury. PLoS One. 2010;5(8):e12106.View ArticlePubMedPubMed CentralGoogle Scholar
- Risbano MG, Meadows CA, Coldren CD, Jenkins TJ, Edwards MG, Collier D, Huber W, Mack DG, Fontenot AP, Geraci MW, Bull TM. Altered immune phenotype in peripheral blood cells of patients with scleroderma-associated pulmonary hypertension. Clin Transl Sci. 2010;3(5):210–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Reich M, Liefeld T, Gould J, Lerner J, Tamayo P, Mesirov JP. GenePattern 2.0. Nat Genet. 2006;38(5):500–1.View ArticlePubMedGoogle Scholar
- Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007;8(1):118–27.View ArticlePubMedGoogle Scholar
- Horvath S, Dong J. Geometric interpretation of gene coexpression network analysis. PLoS Comput Biol. 2008;4(8):e1000117.View ArticlePubMedPubMed CentralGoogle Scholar
- Newman ME. Modularity and community structure in networks. Proc Natl Acad Sci U S A. 2006;103(23):8577–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Newman ME. Communities, modules and large-scale structure in networks. Nat Phys. 2012;8(1):25–31.View ArticleGoogle Scholar
- Csardi G, Nepusz T. The igraph software package for complex network research. Inter J Complex Syst. 2006;1695(5):1–9.Google Scholar
- Reichardt J, Bornholdt S. Statistical mechanics of community detection. Phys Rev E. 2006;74(1):016110.View ArticleGoogle Scholar
- Kanehisa M, Goto S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000;28(1):27–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Croft D, Mundo AF, Haw R, Milacic M, Weiser J, Wu G, Caudy M, Garapati P, Gillespie M, Kamdar MR. The Reactome pathway knowledgebase. Nucleic Acids Res. 2014;42(D1):D472–7.View ArticlePubMedGoogle Scholar
- Huttenhower C, Schroeder M, Chikina MD, Troyanskaya OG. The Sleipnir library for computational functional genomics. Bioinformatics. 2008;24(13):1559–61.View ArticlePubMedPubMed CentralGoogle Scholar
- Reimand J, Arak T, and Vilo J. g: Profiler-a web server for functional interpretation of gene lists (2011 update). Nucleic Acids Res. 2011;39(suppl_2):W307–W315.
- Bastian M, Heymann S, Jacomy M. Gephi: an open source software for exploring and manipulating networks. ICWSM. 2009;8:361–2.Google Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Xue J, Schmidt SV, Sander J, Draffehn A, Krebs W, Quester I, De Nardo D, Gohel TD, Emde M, Schmidleithner L, Ganesan H, Nino-Castro A, Mallmann MR, Labzin L, Theis H, Kraut M, Beyer M, Latz E, Freeman TC, Ulas T, Schultze JL. Transcriptome-based network analysis reveals a spectrum model of human macrophage activation. Immunity. 2014;40(2):274–88.View ArticlePubMedPubMed CentralGoogle Scholar
- Assassi S, Swindell WR, Wu M, Tan FD, Khanna D, Furst DE, Tashkin DP, Jahan-Tigh RR, Mayes MD, Gudjonsson JE, and Chang JT. Dissecting the heterogeneity of skin gene expression patterns in systemic sclerosis. Arthritis Rheumatol. 2015;67(11):3016–26.
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Cho JH, Gelinas R, Wang K, Etheridge A, Piper MG, Batte K, Dakhallah D, Price J, Bornman D, Zhang S, Marsh C, Galas D. Systems biology of interstitial lung diseases: integration of mRNA and microRNA expression changes. BMC Med Genomics. 2011;4:8.View ArticlePubMedPubMed CentralGoogle Scholar
- Barabási A-L, Gulbahce N, Loscalzo J. Network medicine: a network-based approach to human disease. Nat Rev Genet. 2011;12(1):56–68.View ArticlePubMedPubMed CentralGoogle Scholar
- Li Y, Jiang D, Liang J, Meltzer EB, Gray A, Miura R, Wogensen L, Yamaguchi Y, Noble PW. Severe lung fibrosis requires an invasive fibroblast phenotype regulated by hyaluronan and CD44. J Exp Med. 2011;208(7):1459–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Grove LM, Southern BD, Jin TH, White KE, Paruchuri S, Harel E, Wei Y, Rahaman SO, Gladson CL, Ding Q, Craik CS, Chapman HA, Olman MA. Urokinase-type plasminogen activator receptor (uPAR) ligation induces a raft-localized integrin signaling switch that mediates the hypermotile phenotype of fibrotic fibroblasts. J Biol Chem. 2014;289(18):12791–804.View ArticlePubMedPubMed CentralGoogle Scholar
- Martinez FO, Helming L, Milde R, Varin A, Melgert BN, Draijer C, Thomas B, Fabbri M, Crawshaw A, Ho LP, Ten Hacken NH, Cobos Jiménez V, Kootstra NA, Hamann J, Greaves DR, Locati M, Mantovani A, Gordon S. Genetic programs expressed in resting and IL-4 alternatively activated mouse and human macrophages: similarities and differences. Blood. 2013;121(9):e57–69.View ArticlePubMedGoogle Scholar
- Uhlen M, Oksvold P, Fagerberg L, Lundberg E, Jonasson K, Forsberg M, Zwahlen M, Kampf C, Wester K, Hober S. Towards a knowledge-based human protein atlas. Nat Biotechnol. 2010;28(12):1248–50.View ArticlePubMedGoogle Scholar
- Koslowski R, Knoch K, Kuhlisch E, Seidel D, Kasper M. Cathepsins in bleomycin-induced lung injury in rat. Eur Respir J. 2003;22(3):427–35.View ArticlePubMedGoogle Scholar
- Nacu N, Luzina IG, Highsmith K, Lockatell V, Pochetuhen K, Cooper ZA, Gillmeister MP, Todd NW, Atamas SP. Macrophages produce TGF-beta-induced (beta-ig-h3) following ingestion of apoptotic cells and regulate MMP14 levels and collagen turnover in fibroblasts. J Immunol. 2008;180(7):5036–44.View ArticlePubMedPubMed CentralGoogle Scholar
- Seyrantepe V, Iannello A, Liang F, Kanshin E, Jayanth P, Samarani S, Szewczuk MR, Ahmad A, Pshezhetsky AV. Regulation of phagocytosis in macrophages by neuraminidase 1. J Biol Chem. 2010;285(1):206–15.View ArticlePubMedGoogle Scholar
- Jackman HL, Tan F, Schraufnagel D, Dragović T, Dezsö B, Becker RP, Erdös EG. Plasma membrane-bound and lysosomal peptidases in human alveolar macrophages. Am J Respir Cell Mol Biol. 1995;13(2):196–204.View ArticlePubMedGoogle Scholar
- Bell-Temin H, Culver-Cochran AE, Chaput D, Carlson CM, Kuehl M, Burkhardt BR, Bickford PC, Liu B, Stevens SM. Novel molecular insights into classical and alternative activation states of microglia as revealed by stable isotope labeling by amino acids in cell culture (SILAC)-based proteomics. Mol Cell Proteomics. 2015;14(12):3173–84.View ArticlePubMedPubMed CentralGoogle Scholar
- Gautier EL, Shay T, Miller J, Greter M, Jakubzick C, Ivanov S, Helft J, Chow A, Elpek KG, Gordonov S, Mazloom AR, Ma'ayan A, Chua WJ, Hansen TH, Turley SJ, Merad M, Randolph GJ, Immunological Genome Consortium. Gene-expression profiles and transcriptional regulatory pathways that underlie the identity and diversity of mouse tissue macrophages. Nat Immunol. 2012;13(11):1118–28.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang J, Nikrad MP, Travanty EA, Zhou B, Phang T, Gao B, Alford T, Ito Y, Nahreini P, Hartshorn K. Innate immune response of human alveolar macrophages during influenza A infection. PLoS One. 2012;7(3):e29879.View ArticlePubMedPubMed CentralGoogle Scholar
- Higashi-Kuwata N, Makino T, Inoue Y, Takeya M, Ihn H. Alternatively activated macrophages (M2 macrophages) in the skin of patient with localized scleroderma. Exp Dermatol. 2009;18(8):727–9.View ArticlePubMedGoogle Scholar
- Hannaford J, Guo H, Chen X. Involvement of cathepsins B and L in inflammation and cholesterol trafficking protein NPC2 secretion in macrophages. Obesity (Silver Spring). 2013;21(8):1586–95.View ArticleGoogle Scholar
- Romero F, Shah D, Duong M, Penn RB, Fessler MB, Madenspacher J, Stafstrom W, Kavuru M, Lu B, Kallen CB, Walsh K, and Summer R. A pneumocyte-macrophage paracrine lipid axis drives the lung toward fibrosis. Am J Respir Cell Mol Biol. 2015;53(1):74–86.
- Clavel C, Ceccato L, Anquetil F, Serre G, and Sebbag M. Among human macrophages polarised to different phenotypes, the M-CSF-oriented cells present the highest pro-inflammatory response to the rheumatoid arthritis-specific immune complexes containing ACPA. Ann Rheum Dis. 2016
- Stubbs VE, Power C, Patel KD. Regulation of eotaxin-3/CCL26 expression in human monocytic cells. Immunology. 2010;130(1):74–82.View ArticlePubMedPubMed CentralGoogle Scholar
- O'Neill LA, Pearce EJ. Immunometabolism governs dendritic cell and macrophage function. J Exp Med. 2016;213(1):15–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Christmann RB, Lafyatis R. The cytokine language of monocytes and macrophages in systemic sclerosis. Arthritis Res Ther. 2010;12(5):1.View ArticleGoogle Scholar
- Johnson ME, Pioli PA, and Whitfield ML. Gene expression profiling offers insights into the role of innate immune signaling in SSc. Semin Immunopathol. 2015;37(5):501–9.
- Scala E, Pallotta S, Frezzolini A, Abeni D, Barbieri C, Sampogna F, De Pità O, Puddu P, Paganelli R, Russo G. Cytokine and chemokine levels in systemic sclerosis: relationship with cutaneous and internal organ involvement. Clin Exp Immunol. 2004;138(3):540–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Riccieri V, Rinaldi T, Spadaro A, Scrivo R, Ceccarelli F, Franco MD, Taccari E, Valesini G. Interleukin-13 in systemic sclerosis: relationship to nailfold capillaroscopy abnormalities. Clin Rheumatol. 2003;22(2):102–6.View ArticlePubMedGoogle Scholar
- Greenblatt MB, Sargent JL, Farina G, Tsang K, Lafyatis R, Glimcher LH, Whitfield ML, Aliprantis AO. Interspecies comparison of human and murine scleroderma reveals IL-13 and CCL2 as disease subset-specific targets. Am J Pathol. 2012;180(3):1080–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Kaviratne M, Hesse M, Leusink M, Cheever AW, Davies SJ, McKerrow JH, Wakefield LM, Letterio JJ, Wynn TA. IL-13 activates a mechanism of tissue fibrosis that is completely TGF-beta independent. J Immunol. 2004;173(6):4020–9.View ArticlePubMedGoogle Scholar
- Fuschiotti P, Medsger TA, Morel PA. Effector CD8+ T cells in systemic sclerosis patients produce abnormally high levels of interleukin-13 associated with increased skin fibrosis. Arthritis Rheum. 2009;60(4):1119–28.View ArticlePubMedGoogle Scholar
- Granel B, Allanore Y, Chevillard C, Arnaud V, Marquet S, Weiller PJ, Durand JM, Harlé JR, Grange C, Frances Y, Berbis P, Gaudart J, de Micco P, Kahan A, Dessein A. IL13RA2 gene polymorphisms are associated with systemic sclerosis. J Rheumatol. 2006;33(10):2015–9.PubMedGoogle Scholar
- Hasegawa M, Sato S, Nagaoka T, Fujimoto M, Takehara K. Serum levels of tumor necrosis factor and interleukin-13 are elevated in patients with localized scleroderma. Dermatology. 2003;207(2):141–7.View ArticlePubMedGoogle Scholar
- Christmann RB, Hayes E, Pendergrass S, Padilla C, Farina G, Affandi AJ, Whitfield ML, Farber HW, Lafyatis R. Interferon and alternative activation of monocyte/macrophages in systemic sclerosis-associated pulmonary arterial hypertension. Arthritis Rheum. 2011;63(6):1718–28.View ArticlePubMedPubMed CentralGoogle Scholar
- Okabe Y, Medzhitov R. Tissue-specific signals control reversible program of localization and functional polarization of macrophages. Cell. 2014;157(4):832–44.View ArticlePubMedPubMed CentralGoogle Scholar
- Gosselin D, Link VM, Romanoski CE, Fonseca GJ, Eichenfield DZ, Spann NJ, Stender JD, Chun HB, Garner H, Geissmann F, Glass CK. Environment drives selection and function of enhancers controlling tissue-specific macrophage identities. Cell. 2014;159(6):1327–40.View ArticlePubMedPubMed CentralGoogle Scholar
- Khanna D, Denton CP, Jahreis A, van Laar JM, Frech TM, Anderson ME, Baron M, Chung L, Fierlbeck G, Lakshminarayanan S, Allanore Y, Pope JE, Riemekasten G, Steen V, Müller-Ladner U, Lafyatis R, Stifano G, Spotswood H, Chen-Harris H, Dziadek S, Morimoto A, Sornasse T, Siegel J, Furst DE. Safety and efficacy of subcutaneous tocilizumab in adults with systemic sclerosis (faSScinate): a phase 2, randomised, controlled trial. Lancet. 2016;387(10038):2630–40.View ArticlePubMedGoogle Scholar
- Gerber EE, Gallo EM, Fontana SC, Davis EC, Wigley FM, Huso DL, Dietz HC. Integrin-modulating therapy prevents fibrosis and autoimmunity in mouse models of scleroderma. Nature. 2013;503(7474):126–30.View ArticlePubMedPubMed CentralGoogle Scholar
- McGovern N, Schlitzer A, Gunawan M, Jardine L, Shin A, Poyner E, Green K, Dickinson R, Wang XN, Low D, Best K, Covins S, Milne P, Pagan S, Aljefri K, Windebank M, Saavedra DM, Larbi A, Wasan PS, Duan K, Poidinger M, Bigley V, Ginhoux F, Collin M, Haniffa M. Human dermal CD14+ cells are a transient population of monocyte-derived macrophages. Immunity. 2014;41(3):465–77.View ArticlePubMedPubMed CentralGoogle Scholar
- Thomas AC, Eijgelaar WJ, Daemen MJ, Newby AC. The pro-fibrotic and anti-inflammatory foam cell macrophage paradox. Genom Data. 2015;6:136–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Thomas AC, Eijgelaar WJ, Daemen MJ, Newby AC. Foam cell formation in vivo converts macrophages to a pro-fibrotic phenotype. PLoS One. 2015;10(7):e0128163.View ArticlePubMedPubMed CentralGoogle Scholar
- Eitzman DT, McCoy RD, Zheng X, Fay WP, Shen T, Ginsburg D, Simon RH. Bleomycin-induced pulmonary fibrosis in transgenic mice that either lack or overexpress the murine plasminogen activator inhibitor-1 gene. J Clin Investig. 1996;97(1):232.View ArticlePubMedPubMed CentralGoogle Scholar
- Günther A, Lübke N, Ermert M, Schermuly RT, Weissmann N, Breithecker A, Markart P, Ruppert C, Quanz K, Ermert L. Prevention of bleomycin-induced lung fibrosis by aerosolization of heparin or urokinase in rabbits. Am J Respir Crit Care Med. 2003;168(11):1358–65.View ArticlePubMedGoogle Scholar
- Chambers RC. Abnormal wound healing responses in pulmonary fibrosis: focus on coagulation signalling. Eur Respir Rev. 2008;17(109):130–7.View ArticleGoogle Scholar
- Horowitz JC, Rogers DS, Simon RH, Sisson TH, Thannickal VJ. Plasminogen activation induced pericellular fibronectin proteolysis promotes fibroblast apoptosis. Am J Respir Cell Mol Biol. 2008;38(1):78–87.View ArticlePubMedGoogle Scholar
- Moodley YP, Caterina P, Scaffidi AK, Misso NL, Papadimitriou JM, McAnulty RJ, Laurent GJ, Thompson PJ, Knight DA. Comparison of the morphological and biochemical changes in normal human lung fibroblasts and fibroblasts derived from lungs of patients with idiopathic pulmonary fibrosis during FasL-induced apoptosis. J Pathol. 2004;202(4):486–95.View ArticlePubMedGoogle Scholar
- Uhal BD. The role of apoptosis in pulmonary fibrosis. Eur Respir Rev. 2008;17(109):138–44.View ArticleGoogle Scholar
- Liu F, Mih JD, Shea BS, Kho AT, Sharif AS, Tager AM, Tschumperlin DJ. Feedback amplification of fibrosis through matrix stiffening and COX-2 suppression. J Cell Biol. 2010;190(4):693–706.View ArticlePubMedPubMed CentralGoogle Scholar
- Morimoto K, Janssen WJ, Terada M. Defective efferocytosis by alveolar macrophages in IPF patients. Respir Med. 2012;106(12):1800–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Vachon E, Martin R, Plumb J, Kwok V, Vandivier RW, Glogauer M, Kapus A, Wang X, Chow C-W, Grinstein S, Downey GP. CD44 is a phagocytic receptor. Blood. 2006;107(10):4149–58.View ArticlePubMedGoogle Scholar
- Noguera A, Gomez C, Faner R, Cosio B, González-Périz A, Clària J, Carvajal A, Agustí A. An investigation of the resolution of inflammation (catabasis) in COPD. Respir Res. 2012;13:101.View ArticlePubMedPubMed CentralGoogle Scholar
- Jung KK, Liu XW, Chirco R, Fridman R, Kim HR. Identification of CD63 as a tissue inhibitor of metalloproteinase-1 interacting cell surface protein. EMBO J. 2006;25(17):3934–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Matsumoto T, Matsumori H, Taki T, Takagi T, Fukuda Y. Infantile GM1-gangliosidosis with marked manifestation of lungs. Acta Pathol Jpn. 1979;29(2):269–76.PubMedGoogle Scholar
- Roy N, Deveraux QL, Takahashi R, Salvesen GS, Reed JC. The c-IAP-1 and c-IAP-2 proteins are direct inhibitors of specific caspases. EMBO J. 1997;16(23):6914–25.View ArticlePubMedPubMed CentralGoogle Scholar
- Deveraux QL, Roy N, Stennicke HR, Van Arsdale T, Zhou Q, Srinivasula SM, Alnemri ES, Salvesen GS, Reed JC. IAPs block apoptotic events induced by caspase-8 and cytochrome c by direct inhibition of distinct caspases. EMBO J. 1998;17(8):2215–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Todorovic V, Chen CC, Hay N, Lau LF. The matrix protein CCN1 (CYR61) induces apoptosis in fibroblasts. J Cell Biol. 2005;171(3):559–68.View ArticlePubMedGoogle Scholar
- Juric V, Chen C-C, Lau LF. Fas-mediated apoptosis is regulated by the extracellular matrix protein CCN1 (CYR61) in vitro and in vivo. Mol Cell Biol. 2009;29(12):3266–79.View ArticlePubMedPubMed CentralGoogle Scholar
- Franzen CA, Chen CC, Todorović V, Juric V, Monzon RI, Lau LF. Matrix protein CCN1 is critical for prostate carcinoma cell proliferation and TRAIL-induced apoptosis. Mol Cancer Res. 2009;7(7):1045–55.View ArticlePubMedPubMed CentralGoogle Scholar
- Rössig L, Haendeler J, Hermann C, Malchow P, Urbich C, Zeiher AM, Dimmeler S. Nitric oxide down-regulates MKP-3 mRNA levels: involvement in endothelial cell protection from apoptosis. J Biol Chem. 2000;275(33):25502–7.View ArticlePubMedGoogle Scholar
- Azab NA, Rady HM, Marzouk SA. Elevated serum TRAIL levels in scleroderma patients and its possible association with pulmonary involvement. Clin Rheumatol. 2012;31(9):1359–64.View ArticlePubMedGoogle Scholar
- Prunier C, Howe PH. Disabled-2 (Dab2) is required for transforming growth factor beta-induced epithelial to mesenchymal transition (EMT). J Biol Chem. 2005;280(17):17540–8.View ArticlePubMedGoogle Scholar
- Penheiter SG, Singh RD, Repellin CE, Wilkes MC, Edens M, Howe PH, Pagano RE, Leof EB. Type II transforming growth factor-beta receptor recycling is dependent upon the clathrin adaptor protein Dab2. Mol Biol Cell. 2010;21(22):4009–19.View ArticlePubMedPubMed CentralGoogle Scholar
- Smyth N, Vatansever HS, Murray P, Meyer M, Frie C, Paulsson M, Edgar D. Absence of basement membranes after targeting the LAMC1 gene results in embryonic lethality due to failure of endoderm differentiation. J Cell Biol. 1999;144(1):151–60.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar