- Open Access
Microarray profiling reveals suppressed interferon stimulated gene program in fibroblasts from scleroderma-associated interstitial lung disease
Respiratory Researchvolume 14, Article number: 80 (2013)
Interstitial lung disease is a major cause of morbidity and mortality in systemic sclerosis (SSc), with insufficiently effective treatment options. Progression of pulmonary fibrosis involves expanding populations of fibroblasts, and the accumulation of extracellular matrix proteins. Characterisation of SSc lung fibroblast gene expression profiles underlying the fibrotic cell phenotype could enable a better understanding of the processes leading to the progressive build-up of scar tissue in the lungs. In this study we evaluate the transcriptomes of fibroblasts isolated from SSc lung biopsies at the time of diagnosis, compared with those from control lungs.
We used Affymetrix oligonucleotide microarrays to compare the gene expression profile of pulmonary fibroblasts cultured from 8 patients with pulmonary fibrosis associated with SSc (SSc-ILD), with those from control lung tissue peripheral to resected cancer (n=10). Fibroblast cultures from 3 patients with idiopathic pulmonary fibrosis (IPF) were included as a further comparison. Genes differentially expressed were identified using two separate analysis programs following a set of pre-determined criteria: only genes significant in both analyses were considered. Microarray expression data was verified by qRT-PCR and/or western blot analysis.
A total of 843 genes were identified as differentially expressed in pulmonary fibroblasts from SSc-ILD and/or IPF compared to control lung, with a large overlap in the expression profiles of both diseases. We observed increased expression of a TGF-β response signature including fibrosis associated genes and myofibroblast markers, with marked heterogeneity across samples. Strongly suppressed expression of interferon stimulated genes, including antiviral, chemokine, and MHC class 1 genes, was uniformly observed in fibrotic fibroblasts. This expression profile includes key regulators and mediators of the interferon response, such as STAT1, and CXCL10, and was also independent of disease group.
This study identified a strongly suppressed interferon-stimulated gene program in fibroblasts from fibrotic lung. The data suggests that the repressed expression of interferon-stimulated genes may underpin critical aspects of the profibrotic fibroblast phenotype, identifying an area in pulmonary fibrosis that requires further investigation.
Pulmonary fibrosis, characterised by the destruction of lung architecture leading to organ failure, is, together with pulmonary hypertension (PH), the leading cause of death in patients with systemic sclerosis (SSc) . Interstitial lung disease is more common in SSc (SSc-ILD) than in any other connective tissue disease, occurring in > 70% of patients , and is most frequently associated with a pattern of non-specific interstitial pneumonia (NSIP) . Despite intense research efforts, the underlying mechanisms of SSc-ILD remain largely unknown , and there are currently limited therapeutic options for this serious complication .
While a large number of hypothesis-driven studies have identified potential profibrotic mediators [4, 5], translation of these into therapeutic targets has so far been largely disappointing . The search for more effective targets in lung fibrosis is now being addressed using hypothesis generating microarray-based strategies [7, 8]. The majority of these studies have investigated gene expression in whole lung tissue samples, mostly in the idiopathic setting . Matrix metalloproteinase (MMP) 7 , osteopontin , Twist1 , and MMP19 , are among suggested mediators identified using this strategy in idiopathic pulmonary fibrosis (IPF), a disease characterised by a histological pattern of usual interstitial pneumonia (UIP) [13, 14]. In SSc, most microarray studies have been performed on skin biopsies and dermal fibroblasts . However, a recent study compared whole lung tissue and fibroblasts isolated, at the time of transplant, from SSc-ILD lungs with a histological pattern of UIP, with those from IPF and idiopathic PH. The investigators reported gene profiles of SSc-ILD/UIP, with either predominant fibrosis or PH, overlapping with profiles of IPF and idiopathic PH, respectively .
While the initiating factors for fibrosis development may vary between diseases, the progressive accumulation of scar tissue in the lung is characterised by common themes, including expanding populations of activated fibroblasts, and excessive accumulation of extracellular matrix (ECM) proteins . An important strategy to identify potential therapeutic targets, therefore, is to define fibrotic fibroblast phenotypes so as to delineate underlying key mechanisms for fibrosis progression.
Here we report analysis of the transcriptome of fibroblasts isolated from surgical lung biopsies at the time of diagnosis, from patients with well defined SSc-ILD and the histopathological pattern of NSIP. Although the main aim of this study was to compare SSc-ILD/NSIP fibroblast gene expression profiles with those of control lung fibroblasts, we also included a small number of IPF-derived fibroblast lines, as a separate fibrotic group. Our study confirms, with a robust signature in both diseases, the aberrant expression of previously reported myofibroblast markers and fibrosis mediators, and identifies a number of novel, co-expressed putative disease targets. We also observed the suppression of a large gene program, the interferon stimulated genes (ISGs), reported here for the first time. From the known function of some of these genes , it is possible to hypothesise that this suppressed program underlies key fibrotic fibroblast properties, such as hyper-proliferation, and apoptosis resistance. This study therefore identifies a potential new area for investigation and possible intervention in pulmonary fibrosis.
Patients and primary lung fibroblasts
Primary adult pulmonary fibroblasts were cultured from control tissue samples of unaffected lung from patients undergoing cancer-resection surgery (n=10), and from surgical lung biopsy samples of 11 patients with pulmonary fibrosis (SSc-ILD n=8 and IPF n=3). Independent reviews of the clinical and histopathologic diagnoses were performed and conformed to established criteria [17, 18]. All of the SSc-ILD biopsies were characterised by a fibrotic NSIP pattern, and the IPF biopsies by a UIP pattern, based on current consensus criteria for these histological patterns . The control tissue was histologically normal. Median age (range) was 60 (52–78) in controls, 48 (38–69) in SSc-ILD, and 61 (44–67) in IPF. The gender distribution (M/F) was as follows: controls 6/4; SSc-ILD 2/6; IPF 2/1. Four of the SSc-ILD and two of the IPF patients were ex-smokers. Smoking status was not available for all control cases. In SSc-ILD patients, median (range) percent predicted FVC was 72.5% (61–106), median FEV1 was 79% (58–92) and median DLCO was 50% (24–58). In IPF patients, median FVC was 70% (64–75), median FEV1 was 66% (55–79), and median DLCO was 50% (35–53). Patients had not been on corticosteroids or other immunosuppressants prior to surgical biopsy, as the biopsy was performed at the time of diagnosis of the ILD pattern, prior to initiation of treatment. Informed written consent was given by all subjects, and authorisation given by the Royal Brompton Hospital Ethics Committee. Fibroblasts were obtained from the biopsies by explant culture, and cell cultures maintained, as previously described [20, 21]. Experiments were performed on fibroblasts at passage 2–5. Only one sample (S1) was used at passage 2. There was no difference in the median passage number between the control (median: 4.5; range: 3–5), SSc-ILD (median: 4; range: 2–5) and IPF (median: 4; range 3–5) groups.
Microarray gene expression and analysis
At confluence, lung fibroblasts were serum-deprived for 42 hours (media changed at 18 hours) in the presence of 0.1% bovine serum albumin (Sigma). Total RNA was harvested (Trizol, Life Technologies), quantified, and the integrity verified by denaturing gel electrophoresis. Samples with a 28S:18S ratios of approximately 2:1 were accepted for further analysis by the Genomics Laboratory, CSC-MRC, Imperial College London, Hammersmith. RNA samples were prepared for chip hybridisation following manufacturer’s guidelines (Affymetrix). Hybridisation of cRNA to Affymetrix human U133Av2 chips, containing approximately 14,500 well characterised human genes, signal amplification, and data collection were performed using an Affymetrix fluidics station and chip reader, according to manufacturer’s protocol. Array normalisation, using the invariant set normalisation method, and subsequent calculation of model-based expression values, was performed using DNA-Chip Analyzer (dChip) . A list of differentially expressed genes was generated in dChip using fold change ≥2, difference in means ≥100, and p<0.05. Significance analysis of microarrays (SAM) v 4.0  was also used to determine significantly differentially expressed genes with fold change ≥2, difference in means ≥100, delta =1, and a false discovery rate <0.01. Only genes identified as differentially expressed according to both programs were considered as different between groups. Microarray data has been deposited in the Gene Expression Omnibus database , accession number GSE40839. Although the main aim of this study was to assess global gene expression profiles in SSc-ILD compared to controls, for completeness we present the comparison between both SSc-ILD and IPF and controls separately in tables. dChip software was used for data visualisation and hierarchical average linkage clustering using Pearson’s correlation .
Functional category analysis
Functional categories enriched in the differentially expressed genes were identified using the functional annotation and clustering tool of the Database for Annotation, Visualisation, and Integrated Discovery (DAVID) v 6.7 [25, 26]. The probability that a Gene Ontology (GO) biological process term  is overrepresented was determined by a modified Fisher’s exact test, comparing the proportion of genes in the whole genome which are part of that GO term, to the proportion of the differentially expressed genes which are part of the same GO term, and was expressed as an EASE score. Clusters of overrepresented GO terms were then generated based on the similarity of differentially expressed genes assigned to each functional GO term. Clusters were considered significantly overrepresented if they contained a minimum of five GO terms with an EASE score of ≤ 0.01. A summary description of each cluster was generated based on the constitutive GO term names of that cluster which achieved an EASE score <0.05 following Benjamini-Hotchberg correction of multiple comparisons. Only clusters with enrichment scores >3 (minus log transformed geometric mean of the EASE scores of the constitutive terms, equivalent to average EASE score=0.001) were selected. The open access database INTERFEROME  was used to identify differentially expressed genes which have been shown experimentally to be regulated by interferons.
RNA was extracted using the RNeasy® Mini kit (Qiagen) according to manufacturer’s instructions. Samples were quantified and quality assessed using the NanoDrop spectrophotometry system (Thermo Scientific). Complementary DNA (cDNA) was synthesised from 500 ng RNA in a 20 μl reaction using the QuanTitect® reverse transcription kit (Qiagen). Expression levels were measured using a Rotor Gene 6000 (Corbett) in 10 μl reactions containing 2 μl cDNA (five-fold dilution), 1 × SensiMix™ SYBR NO-ROX (Bioline), and 0.5 μM of each forward and reverse primer (Table 1). PCR conditions were: 10 minutes at 95°C, followed by 40 cycles of 10 seconds at 95°C, 15 seconds at 57°C, and 5 seconds at 72°C. All reactions were performed in duplicate and non-template controls were included for each gene. Standard curves were generated for each gene studied using seven two-fold serial dilutions, high standard of 1×107 copies/μl, of primer set amplicons generated from cDNA. Threshold cycle was manually determined at a fixed value of 10-0.5 and the template quantity calculated using Rotor Gene 6000 Series Software 1.7 (Corbett). Expression levels were normalised to YWHAZ and HPRT1[29, 30].
Western blot analysis
Following pre-incubation for 24 hours in serum-free media (DMEM, 0.1% BSA, penicillin/streptomycin), pulmonary fibroblasts from healthy controls, SSc-ILD, and IPF (n=3 for each phenotype), were cultured for a further 24 hours in fresh serum-free media. Cells were lysed and western blot analysis was performed using the following primary antibodies: CTGF and STAT1 (Santa Cruz Biotechnology); αSMA (Dako); IFITM1-3, ISG15 and GAPDH (Abcam); IRF-1 (Cell Signaling Technology); horseradish peroxidise conjugated secondary antibodies (Dako and Cell Signaling Technology); and ECL detection (Amersham).
Gene expression profiles of fibrotic lung fibroblasts: approximately two-thirds of differentially expressed genes are down-regulated
Using an Affymetrix platform (U133Av2), we determined basal (serum free) global gene expression levels in fibroblasts prepared from lung tissue of 8 patients with SSc-ILD and 10 control lungs. As a further comparison we also included 3 fibroblast cultures from lung tissue of IPF patients. Unsupervised hierarchical cluster analysis of samples and genes resulted in an overall separation of fibrotic samples from controls (Figure 1). Two of the SSc-ILD samples clustered among the normal controls, demonstrating recognised fibrotic fibroblast sample heterogeneity.
When stringent criteria were applied, as described in the methods section, 478 and 744 probe sets (probes in future), equivalent to 360 and 547 genes, displayed differential expression by at least two-fold in SSc-ILD and IPF fibroblasts, respectively, compared to control fibroblasts. In SSc-ILD fibroblasts, 125 probes (99 genes) showed significantly higher, and 353 probes (261 genes) lower expression levels compared with controls (Additional file 1). In IPF fibroblasts, 239 and 505 probes (181 and 366 genes) had significantly higher and lower expression, respectively, compared with control fibroblasts (Additional file 2). Thus, approximately two thirds of differentially expressed genes were expressed at lower levels in fibroblasts from the two disease groups compared with controls. The sets of underexpressed genes in the disease groups, compared to normal controls, also contained the most significantly differentially expressed genes. Table 2 lists the 20 most significantly differentially expressed, and highest ranking (fold change) genes in SSc-ILD and IPF samples, separately, versus controls. It should be noted here that, while 379 out of a total of 843 probes were differentially expressed in both SSc-ILD and IPF vs. controls, most genes falling outside this overlap were very nearly significant in the other disease group, i.e. followed similar trends without meeting cut-off criteria. When the two disease groups were compared directly, only 8 probes (7 genes) were significantly differentially expressed (Additional file 3). The 843 probes differentially expressed in SSc-ILD and/or IPF, which also include the genes differentially expressed between these two groups, are listed in Additional file 4. Within the pooled fibrotic samples (IPF and SSc-ILD), no significant difference was observed in the 843 probes according to smoking status (SAM analysis, data not shown).
Expression levels of a subset of genes identified by the present microarray analysis were verified by qRT-PCR in the microarray RNA samples demonstrating good correlation between these two techniques (Figure 2). Protein levels for six differentially expressed genes, fibrosis related genes; connective tissue growth factor (CTGF), and alpha-smooth muscle actin (αSMA), and interferon stimulated genes (ISG); signal transducer and activator of transcription 1 (STAT1), IFITM1-3 (all three isoforms are detected by this antibody), ISG15 and IRF-1, were determined by western blot analysis in independent preparations of additional SSc-ILD, IPF, and control fibroblasts (Figure 3).
Functional groups over-represented among differentially expressed genes
As predicted from the similarities between the two disease groups with regards to differentially expressed genes, gene ontology analysis revealed a major overlap in overrepresented functional groups. Enriched functional groups representing six broad categories (GO clusters) were identified among the genes with higher expression in disease fibroblasts compared with normal controls (see Table 3 for a summary of enriched clusters of GO terms, and Additional file 5 for full GO analysis): anatomical structure development, regulation of cell cycle, response to stress and wounding, regulation of apoptosis, cell migration, and smooth muscle contraction. Among underexpressed/downregulated genes, GO clusters included: inflammatory and immune response, response to biotic stimulus, regulation of apoptosis, regulation of cell migration, regulation of cell proliferation, and regulation of I-κB/NF-κB cascade. By far the most enriched functional groups, compared with control cells, in both SSc-ILD and IPF lung fibroblasts, are genes involved in immune system processes and in microbial/viral defence, which are strongly suppressed in both disease groups. These genes are also among the most significantly differentially expressed genes in this study.
Cluster analysis of differentially expressed genes – ISGs; a major group suppressed in fibrotic lung fibroblasts
To visualise the differential expression across samples and to identify co-expressed genes, average linkage cluster analysis was performed using expression data for the 843 probes sets which displayed differential expression in at least one of the two comparisons: SSc-ILD vs. control, and IPF vs. control (Figure 4A). Within this set are also the eight probes differentially expressed when SSc-ILD and IPF samples were compared directly. Parts of identified gene clusters were selected to illustrate co-expression among upregulated (Figure 4, Panels B-F) and downregulated (Panels H-J) genes, and also to highlight different patterns of sample heterogeneity. Overall, a heterogeneous expression pattern was observed in the upregulated genes, whereas downregulated genes had more uniform expression patterns across samples for the majority of genes in both disease groups. Among upregulated genes, a TGF-β response signature [31, 32] including genes encoding for fibrosis mediators and myofibroblast markers, such as SERPINE1 (PAI1), connective tissue growth factor, smoothelin, and transgelin (SM22) is prominent (Panel B). Co-expressed with these are strongly upregulated genes: growth arrest and DNA-damage inducible β (GADD45), xylosyltransferase 1 (XYLT1), N-cadherin, and elastin, with potential roles in the fibrotic disease process. Groups B, C and D all contain genes involved in contraction and migration, however, the degree of heterogeneity between samples differ between these groups: in group B, the majority of fibrotic samples have elevated expression of smoothelin and transgelin compared with controls; in group C, fewer fibrotic samples, 7 out of 11, have enhanced levels of α2 smooth muscle actin (ACTA2) expression; and in group D, all three IPF samples, but only 3 out of 8 SSc-ILD samples, have elevated expression of calponin 1 and actin gamma 2 smooth muscle (ACTG2). This may indicate different degrees of contractile/migratory phenotypes among these fibrotic cell preparations. Panel E contains ID1 and ID3, which are in the top 20 differentially regulated genes in both disease groups, and are upregulated in most of the fibrotic samples. Group F depicts a cluster of co-expressed cell-cycle associated genes, including cyclins and TOPO2, which exhibit heterogeneous expression in both disease groups. Panel G illustrates an area with less clustering, which however includes possible disease specific genes, e.g. Secreted protein, acidic, cysteine-rich (SPARC) (IPF, Table 2B) and desmoplakin (SSc-ILD) (Additional file 4). Desmoplakin is among the top 20 most upregulated genes in SSc-ILD with an elevated expression in 7 out of the 8 SSc-ILD fibroblast lines, but with low expression in the three IPF cell preparations and in all controls. Desmoplakin is part of the desmosome complex which forms tight cell-cell contacts , and its enhanced expression in SSc-ILD fibroblasts may define a different pathogenesis and cell origin. Panels H-J shows the marked suppressed expression of interferon stimulated genes (ISGs), such as antiviral genes (Group H), chemokine (Group I) and MHC class I genes (Group J). This cluster also includes key regulators of the interferon gene program, STAT1, interferon regulatory factor 1 (IRF1) and interferon regulatory factor 7 (IRF7), as well as chemokine (C-X-C motif) ligand 10 (CXCL10/IP10), one of the most strongly suppressed genes in the study, and the most strongly repressed chemokine in both disease groups (Panel H and I).
Comparing the SSc-ILD and IFP fibroblast gene expression profiles with the interferome
Since we observe a clear fibrosis/TGF-β signature together with a strongly suppressed ISG program, and there is a well-documented antagonistic relationship between TGF-β and interferon signalling in the fibrosis literature , we next interrogated the Interferome, a database of interferon regulated genes (IRGs) reported in the literature . This database includes 1996 human IRGs, of which 1581 are induced and 415 repressed by interferons. In our study, out of the 99 and 181 overexpressed genes in SSc-ILD and IPF fibroblasts, respectively, 40 (40.4%) and 55 (30.1%) genes were in the Interferome database. Among our underexpressed genes, out of a total of 261 and 366 in SSc-ILD and IPF fibroblasts, respectively, 134 (51.3%) and 173 (74.3%) genes were in the Interferome database. The genes overlapping with IRGs in this database are listed in Additional file 6. The comparison revealed that many of the TGF-β responsive genes upregulated in our microarray data set are indeed genes repressed by IFNs.
Previous studies have shown that fibroblasts isolated from SSc-ILD  and IPF  lungs, while displaying substantial heterogeneity, are generally more proliferative, migratory, resistant to apoptosis, and ECM producing, than control lung fibroblasts . A greater proportion display elevated αSMA expression, also enhancing their contractility. These features are all consistent with the fibroblast accumulation and scar tissue formation observed in fibrotic lung. The so called myofibroblast phenotype is maintained over several passages in culture , suggesting an underlying epigenetic regulation, whether established and maintained in local cells within the chronic disease setting, or supporting the phenotype of a specialised infiltrating wound healing cell type . Regardless of cell origin, this feature enables in vitro studies of mechanisms underlying the fibrotic fibroblast phenotype.
In this study we observed a high number of differentially expressed genes between SSc-ILD and/or IPF derived fibroblasts, compared with controls. There was a large overlap between expression profiles of SSc-ILD and IPF fibroblasts, suggesting several common pathways at this stage of the two diseases. Indeed, a direct comparison demonstrated only seven genes with significantly different expression between the two disease groups. Caution should be applied, however, when interpreting data from this direct comparison since only three IPF samples were included here. It is possible that more genes with differential expression between the two disease groups would be identified with a larger numbers of IPF fibroblasts samples. Such studies would also be required to verify the sporadic observations made in this study, such as the elevated expression levels of desmoplakin in SSc-ILD, but not in IPF fibroblasts. Therefore, we stress that the main objective of the study presented here was to gain an overview of potential SSc-ILD target genes, and the study was not designed to detect differences between different fibrotic entities.
Among significantly upregulated genes in both SSc-ILD and IPF/UIP fibroblasts, we identify a recognised fibrosis signature, including smooth muscle actin (ACTA2), CTGF, and PAI1 (SERPINE1), along with genes more recently associated with lung fibrogenesis, including ID1, ID3, IL11, and NOX4. Inhibitor of DNA binding 1 and 3 (ID1 and ID3), are target genes of bone morphogenetic proteins, and control cell differentiation by dominant negative inhibition of helix-loop-helix transcription factors . Chambers et al. described upregulation of ID1 in lung fibroblasts in response to TGF-β, and linkage to the smooth muscle phenotypic switch . Increased ID1 and ID3 gene expression was also observed in lung tissue and fibroblasts from patients with SSc-ILD by Hsu and colleagues, . The ID proteins are overexpressed in many cancers, controlling cell growth and apoptosis, and have been suggested as a therapeutic target [39, 40]. In fibroproliferative stages of fibrosis, high ID1/ID3 expression could maintain fibroblasts in a dedifferentiated, hyperproliferative, and apoptosis resistant state. Another significantly upregulated gene, though somewhat variably in SSc-ILD fibroblasts, is NOX4, encoding a member of the NADPH oxidase (NOX) proteins which generate superoxide by electron transfer to oxygen . Through its involvement in TGF-β-induced fibroblast differentiation into myofibroblasts, two recent studies have suggested a role for NOX4 in IPF, where it is found to be overexpressed . Specific NOX4 inhibitors are now being developed as possible antifibrotic agents . Another potential novel target gene, GADD45B, highly induced and coexpressed with the fibrosis related genes in this study, is a pro-survival factor associated with stress-resistant tumours , and has been found to be upregulated in SSc skin biopsies .
The most striking of our observations is the strongly and uniformly repressed ISG profile in fibrotic fibroblasts. Among these are genes coding for: antiviral protein myxovirus (influenza virus) resistance 1 (MX1), interferon gamma inducible protein p16 (IFI16), 2’,5’-oligoadenylate synthetase 1 (OAS1), the chemokines chemokine (C-X-C motif) ligand 10 (CXCL10/IP10) and chemokine (C-C motif) ligand 11 (CCL11), and antigen presenting MHC I molecules. Key transcriptional regulators of this program, including IRF1, IRF7, and STAT1 are also suppressed. To our knowledge, the reduction in expression of a large set of immune response/interferon related genes has not previously been described in either SSc-ILD or IPF derived fibroblasts, or other fibrotic lung diseases. Global gene expression was recently evaluated in scleroderma whole skin biopsies and matched dermal fibroblasts. While 26 genes were differentially expressed in both scleroderma whole skin and fibroblasts compared to controls, nine were found to be discordant . Interestingly, the majority of the discordant genes were upregulated in skin biopsies but downregulated in fibroblasts, including the ISG genes we identify as repressed in lung fibroblasts, such as MX1, IFI16, intercellular adhesion molecule 1 (ICAM1), and OAS1. While the authors argue that the discordant genes are an indication of a scarce representativeness of skin fibroblast gene expression in vivo, in light of our data, it is likely that they were in fact observing the same phenomenon reported in the current study. As these are mainly interferon regulated genes, the discordance between fibroblasts and whole tissue gene expression could be explained by the known increase of immune cell populations in SSc skin, likely to be overexpressing immunoregulatory genes. The upregulation of MX1, IRF7 and STAT1 in PBMC from SSc patients would support this notion [46, 47]. The downregulation of the ISG program both in SSc-ILD and in IPF fibroblasts suggests that this phenomenon relates to a local fibroblast specific, rather than systemic, profibrotic process, perhaps underpinned by a general susceptibility for tissue fibrosis, common to both diseases.
Additional support for our finding of an aberrantly regulated ISG program in fibrotic lung fibroblasts, comes from work on two of the signature genes from this group, CXCL10 (IP10) and STAT1. IP10 levels were found to be downregulated in IPF lung fibroblasts by Keane et al. . More recently, Coward et al. have shown, again in IPF fibroblasts, that epigenetic dysregulation involving both histone deacetylation and hypermethylation is responsible for targeted repression of IP10. By contrast, Hsu et al. did not report a significant difference in expression in isolated lung fibroblasts , a difference which may relate to experimental design, as discussed below. Further support of a role of suppressed ISGs in pulmonary fibrosis comes from two animal models. Mice deficient in IP10, and in STAT1, displayed enhanced susceptibility to pulmonary fibrosis. Based on these studies, the suppressed ISG program indentified in fibroblasts isolated from SSc-ILD and IPF lung, as presented here, would support enhanced lung fibrosis progression through promoting fibroblast proliferation, migration, and apoptosis resistance. Interestingly, IFN-γ treatment in IPF has failed to show a benefit , and there was a suggestion of worse pulmonary outcomes in a study investigating treatment of SSc patients with IFN-α compared with placebo , with the latter providing indirect support for IFN-related mechanisms involved in organ-specific SSc complications. There could be several possible explanations for these disappointing results, including unexpected adverse effects through circulating cell populations. Activation of pathways downstream of systemically administered interferons is likely to have different direct and indirect effects depending on the cell type and tissue location. The findings shown here add important information to this complexity, and need to be investigated in future detailed mechanistic cell and animal studies.
The signatures observed in our study are remarkably strong, both in terms of fold difference and statistical significance. One possible reason for this is that lung tissue samples were obtained from biopsies at the time of diagnosis, when the disease may be at a relatively early or active wound healing stage, and when perhaps fibroblast proliferation (accumulation) and elastin synthesis, rather than contraction and collagen remodelling, dominate. This is in contrast to the study by Hsu et al., in which gene expression profiles were investigated in SSc lung tissue and fibroblasts from transplant, and therefore possible end-stage, material, also noted by the authors as a potential limitation . Another possibility relates to differences in in vitro culture condition as we employed serum free media before harvest, similarly to others, including Coward et al. , as opposed to in low serum (0.5%) as applied by e.g. Hsu et al. . Serum withdrawal, a form of cellular stress, may evoke the clear differential expression profiles observed in our study. Whereas in the fibrotic fibroblasts an anti-apoptotic survival gene program is maintained, which may be a result of the suppressed ISGs, this gene program may not be subject to repression in the normal fibroblasts.
While an accepted source of control tissue in studies of ILD [10, 12], the use of control fibroblasts from cancer resected specimens, rather than from healthy control subjects (not available for this study), represents a potential limitation. Although obtained from areas of lung with normal histological appearance, differential gene expression in fibroblasts derived from lungs in which cancer has developed cannot be excluded. However, significant gene expression differences were observed in lung cancer associated fibroblasts compared to matched fibroblasts from areas of normal lung from the same patient, suggesting that the cancer associated phenotype of lung fibroblasts is regionally limited to the cancer stroma . A further possible source of bias in gene expression is smoking history. Smoking has been shown to be associated with interstitial fibrotic changes , and is itself likely to cause changes in the expression of certain genes. We observed no significant differences among differentially expressed genes according to smoking status in the pooled fibrotic samples, suggesting that the observed changes were related to the fibrotic lung disease itself rather than to smoking. However, as subgroup numbers were small, and it was not possible to separately analyse the two fibrotic lung diseases, further studies are needed to carefully assess the contribution of smoking to gene expression changes in the context of fibrotic lung disease.
While the general hypothesis, that a repressed interferon stimulated gene program at least in part underpins the fibrotic fibroblast phenotype, will be tested in future studies, it is interesting to note the same phenomenon in several other clinical settings where hyperplasia and apoptosis resistance are key features; certain viruses, including high-risk human papillomaviruses (HPV), have evolved a mechanism to down-regulate ISGs in host cells as an immunoevasive strategy , and persistent HPV infection may lead to cervical cancer development; breast cancer metastasis is promoted by IRF7 silencing ; fibroblasts from patients with Li-Fraumeni syndrome become spontaneously immortalised through the downregulation of interferon pathway genes . Conversely, IRF1 expression reverts the phenotype of oncogenically transformed fibroblasts , and IRF-1 enhancing drugs with tumour suppressing properties are currently being developed . Many similar examples in the literature lead to questions about whether fibrosis is a pre-cancerous state . The repressed, or aberrantly regulated, fibroblast specific interferon response network may therefore be a common necessary determinant allowing lung fibrosis progression to occur.
In summary, in this study comparing gene expression profiles of fibroblasts explanted from fibrotic lung tissue (SSc-ILD and IPF), with control fibroblasts from areas of normal lung, we observe: an overall elevated expression of previously reported fibrosis associated genes, with marked heterogeneity across samples; differentially regulated myofibroblast markers which correlate with the expression heterogeneity between samples; and a strongly suppressed interferon stimulated gene program, uniformly present across fibrotic samples. This suppressed gene program displays both the greatest significance and largest fold differences in expression in our data set. Similarly to functional findings in parallel fields, particularly cancer, this group of genes, and the suppression of their expression, could explain essential aspects of the profibrotic fibroblast phenotype. This hypothesis will need to be tested by future studies, with particular focus on epigenetic silencing as a potential underlying mechanism.
Database for annotation, visualisation, and integrated discovery
Idiopathic pulmonary fibrosis
Interferon regulated genes
Interferon stimulated gene
Nonspecific interstitial pneumonia
Significance analysis of microarrays
Interstitial lung disease in SSc
Usual interstitial pneumonia.
Steen VD, Medsger TA: Changes in causes of death in systemic sclerosis, 1972–2002. Ann Rheum Dis. 2007, 66: 940-944. 10.1136/ard.2006.066068.
Tan A, Denton CP, Mikhailidis DP, Seifalian AM: Recent advances in the diagnosis and treatment of interstitial lung disease in systemic sclerosis (scleroderma): a review. Clin Exp Rheumatol. 2011, 29: S66-S74.
Bouros D, Wells AU, Nicholson AG, Colby TV, Polychronopoulos V, Pantelidis P, Haslam PL, Vassilakis DA, Black CM, du Bois RM: Histopathologic subsets of fibrosing alveolitis in patients with systemic sclerosis and their relationship to outcome. Am J Respir Crit Care Med. 2002, 165: 1581-1586. 10.1164/rccm.2106012.
du Bois RM: Mechanisms of scleroderma-induced lung disease. Proc Am Thorac Soc. 2007, 4: 434-438. 10.1513/pats.200608-152MS.
Scotton CJ, Chambers RC: Molecular targets in pulmonary fibrosis: the myofibroblast in focus. Chest. 2007, 132: 1311-1321. 10.1378/chest.06-2568.
Datta A, Scotton CJ, Chambers RC: Novel therapeutic approaches for pulmonary fibrosis. Br J Pharmacol. 2011, 163: 141-172. 10.1111/j.1476-5381.2011.01247.x.
Sargent JL, Whitfield ML: Capturing the heterogeneity in systemic sclerosis with genome-wide expression profiling. Expert Rev Clin Immunol. 2011, 7: 463-473. 10.1586/eci.11.41.
Kaminski N, Rosas IO: Gene expression profiling as a window into idiopathic pulmonary fibrosis pathogenesis: can we identify the right target genes?. Proc Am Thorac Soc. 2006, 3: 339-344. 10.1513/pats.200601-011TK.
Zuo F, Kaminski N, Eugui E, Allard J, Yakhini Z, Ben-Dor A, Lollini L, Morris D, Kim Y, DeLustro B, Sheppard D, Pardo A, Selman M, Heller RA: Gene expression analysis reveals matrilysin as a key regulator of pulmonary fibrosis in mice and humans. Proc Natl Acad Sci U S A. 2002, 99: 6292-6297. 10.1073/pnas.092134099.
Pardo A, Gibson K, Cisneros J, Richards TJ, Yang Y, Becerril C, Yousem S, Herrera I, Ruiz V, Selman M, Kaminski N: Up-regulation and profibrotic role of osteopontin in human idiopathic pulmonary fibrosis. PLoS Med. 2005, 2: e251-10.1371/journal.pmed.0020251.
Bridges RS, Kass D, Loh K, Glackin C, Borczuk AC, Greenberg S: Gene expression profiling of pulmonary fibrosis identifies Twist1 as an antiapoptotic molecular “rectifier” of growth factor signaling. Am J Pathol. 2009, 175: 2351-2361. 10.2353/ajpath.2009.080954.
Yu G, Kovkarova-Naumovski E, Jara P, Parwani A, Kass D, Ruiz V, Lopez-Otin C, Rosas IO, Gibson KF, Cabrera S, Ramirez R, Yousem SA, Richards TJ, Chensny LJ, Selman M, Kaminski N, Pardo A: MMP19 is a key regulator of lung fibrosis in mice and humans. Am J Respir Crit Care Med. 2012, 186: 752-762. 10.1164/rccm.201202-0302OC.
Maher TM, Wells AU, Laurent GJ: Idiopathic pulmonary fibrosis: multiple causes and multiple mechanisms?. Eur Respir J. 2007, 30: 835-839. 10.1183/09031936.00069307.
Lynch JP, Saggar R, Weigt SS, Zisman DA, White ES: Usual interstitial pneumonia. Semin Respir Crit Care Med. 2006, 27: 634-651. 10.1055/s-2006-957335.
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: 783-794. 10.1002/art.30159.
Saha B, Jyothi PS, Chandrasekar B, Nandi D: Gene modulation and immunoregulatory roles of interferon gamma. Cytokine. 2010, 50: 1-14. 10.1016/j.cyto.2009.11.021.
Subcommittee for scleroderma criteria of the American Rheumatism Association Diagnostic and Therapeutic Criteria Committee: Preliminary criteria for the classification of systemic sclerosis (scleroderma). Arthritis Rheum. 1980, 23: 581-590. 10.1002/art.1780230510.
American Thoracic Society: Idiopathic pulmonary fibrosis: diagnosis and treatment. International consensus statement. American Thoracic Society (ATS), and the European Respiratory Society (ERS). Am J Respir Crit Care Med. 2000, 161: 646-664.
American Thoracic Society/European Respiratory Society International Multidisciplinary Consensus Classification of the Idiopathic Interstitial Pneumonias: This joint statement of the American Thoracic Society (ATS), and the European Respiratory Society (ERS) was adopted by the ATS board of directors, June 2001 and by the ERS Executive Committee, June 2001. Am J Respir Crit Care Med. 2002, 165: 277-304.
Abraham DJ, Shiwen X, Black CM, Sa S, Xu Y, Leask A: Tumor necrosis factor alpha suppresses the induction of connective tissue growth factor by transforming growth factor-beta in normal and scleroderma fibroblasts. J Biol Chem. 2000, 275: 15220-15225. 10.1074/jbc.275.20.15220.
Shi-wen X, Pennington D, Holmes A, Leask A, Bradham D, Beauchamp JR, Fonseca C, du Bois RM, Martin GR, Black CM, Abraham DJ: Autocrine overexpression of CTGF maintains fibrosis: RDA analysis of fibrosis genes in systemic sclerosis. Exp Cell Res. 2000, 259: 213-224. 10.1006/excr.2000.4972.
Li C, Wong WH: Model-based analysis of oligonucleotide arrays: expression index computation and outlier detection. Proc Natl Acad Sci U S A. 2001, 98: 31-36. 10.1073/pnas.98.1.31.
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: 5116-5121. 10.1073/pnas.091062498.
NCBI: GEO Gene Expression Omnibus. [http://www.ncbi.nlm.nih.gov/geo/]
Huang dW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-57.
Huang dW, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37: 1-13. 10.1093/nar/gkn923.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.
Samarajiwa SA, Forster S, Auchettl K, Hertzog PJ: INTERFEROME: the database of interferon regulated genes. Nucleic Acids Res. 2009, 37: D852-D857. 10.1093/nar/gkn732.
Vandesompele J, De PK, Pattyn F, Poppe B, Van RN, De PA, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: RESEARCH0034-
Stock CJ, Leoni P, Shi-wen X, Abraham DJ, Nicholson AG, Wells AU, Renzoni EA, Lindahl GE: Identification of stable housekeeping genes for real-time PCR in human pulmonary fibroblasts [abstract]. Eur Respir J. 2011, 38 (Suppl. 55): 3805-
Chambers RC, Leoni P, Kaminski N, Laurent GJ, Heller RA: Global expression profiling of fibroblast responses to transforming growth factor-beta1 reveals the induction of inhibitor of differentiation-1 and provides evidence of smooth muscle cell phenotypic switching. Am J Pathol. 2003, 162: 533-546. 10.1016/S0002-9440(10)63847-3.
Renzoni EA, Abraham DJ, Howat S, Shi-wen X, Sestini P, Bou-Gharios G, Wells AU, Veeraraghavan S, Nicholson AG, Denton CP, Leask A, Pearson JD, Black CM, Welsh KI, du Bois RM: Gene expression profiling reveals novel TGFbeta targets in adult lung fibroblasts. Respir Res. 2004, 5: 24-10.1186/1465-9921-5-24.
Pieperhoff S, Barth M, Rickelt S, Franke WW: Desmosomal molecules in and out of adhering junctions: normal and diseased States of epidermal, cardiac and mesenchymally derived cells. Dermatol Res Pract. 2010, 2010: 139167-
Leask A, Abraham DJ: TGF-beta signaling and the fibrotic response. FASEB J. 2004, 18: 816-827. 10.1096/fj.03-1273rev.
Shi-wen X, Kennedy L, Renzoni EA, Bou-Gharios G, du Bois RM, Black CM, Denton CP, Abraham DJ, Leask A: Endothelin is a downstream mediator of profibrotic responses to transforming growth factor beta in human lung fibroblasts. Arthritis Rheum. 2007, 56: 4189-4194. 10.1002/art.23134.
Hinz B, Phan SH, Thannickal VJ, Prunotto M, Desmouliere A, Varga J, De WO, Mareel M, Gabbiani G: Recent developments in myofibroblast biology: paradigms for connective tissue remodeling. Am J Pathol. 2012, 180: 1340-1355. 10.1016/j.ajpath.2012.02.004.
Balestrini JL, Chaudhry S, Sarrazy V, Koehler A, Hinz B: The mechanical memory of lung myofibroblasts. Integr Biol (Camb). 2012, 4: 410-421. 10.1039/c2ib00149g.
Phan SH: Genesis of the myofibroblast in lung injury and fibrosis. Proc Am Thorac Soc. 2012, 9: 148-152. 10.1513/pats.201201-011AW.
Perk J, Iavarone A, Benezra R: Id family of helix-loop-helix proteins in cancer. Nat Rev Cancer. 2005, 5: 603-614. 10.1038/nrc1673.
Mern DS, Hasskarl J, Burwinkel B: Inhibition of Id proteins by a peptide aptamer induces cell-cycle arrest and apoptosis in ovarian cancer cells. Br J Cancer. 2010, 103: 1237-1244. 10.1038/sj.bjc.6605897.
Crestani B, Besnard V, Boczkowski J: Signalling pathways from NADPH oxidase-4 to idiopathic pulmonary fibrosis. Int J Biochem Cell Biol. 2011, 43: 1086-1089. 10.1016/j.biocel.2011.04.003.
Laleu B, Gaggini F, Orchard M, Fioraso-Cartier L, Cagnon L, Houngninou-Molango S, Gradia A, Duboux G, Merlot C, Heitz F, Szyndralewiez C, Page P: First in class, potent, and orally bioavailable NADPH oxidase isoform 4 (Nox4) inhibitors for the treatment of idiopathic pulmonary fibrosis. J Med Chem. 2010, 53: 7715-7730. 10.1021/jm100773e.
Engelmann A, Speidel D, Bornkamm GW, Deppert W, Stocking C: Gadd45 beta is a pro-survival factor associated with stress-resistant tumors. Oncogene. 2008, 27: 1429-1438. 10.1038/sj.onc.1210772.
Whitfield ML, Finlay DR, Murray JI, Troyanskaya OG, Chi JT, Pergamenschikov A, McCalmont TH, Brown PO, Botstein D, Connolly MK: Systemic and cell type-specific gene expression patterns in scleroderma skin. Proc Natl Acad Sci U S A. 2003, 100: 12319-12324. 10.1073/pnas.1635114100.
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-1973. 10.1002/art.21894.
Tan FK, Zhou X, Mayes MD, Gourh P, Guo X, Marcum C, Jin L, Arnett FC: Signatures of differentially regulated interferon gene expression and vasculotrophism in the peripheral blood cells of systemic sclerosis patients. Rheumatology (Oxford). 2006, 45: 694-702. 10.1093/rheumatology/kei244.
York MR, Nagai T, Mangini AJ, Lemaire R, van Seventer JM, Lafyatis R: A macrophage marker, Siglec-1, is increased on circulating monocytes in patients with systemic sclerosis and induced by type I interferons and toll-like receptor agonists. Arthritis Rheum. 2007, 56: 1010-1020. 10.1002/art.22382.
Keane MP, Arenberg DA, Lynch JP, Whyte RI, Iannettoni MD, Burdick MD, Wilke CA, Morris SB, Glass MC, DiGiovine B, Kunkel SL, Strieter RM: The CXC chemokines, IL-8 and IP-10, regulate angiogenic activity in idiopathic pulmonary fibrosis. J Immunol. 1997, 159: 1437-1443.
Coward WR, Watts K, Feghali-Bostwick CA, Jenkins G, Pang L: Repression of IP-10 by interactions between histone deacetylation and hypermethylation in idiopathic pulmonary fibrosis. Mol Cell Biol. 2010, 30: 2874-2886. 10.1128/MCB.01527-09.
Tager AM, Kradin RL, LaCamera P, Bercury SD, Campanella GS, Leary CP, Polosukhin V, Zhao LH, Sakamoto H, Blackwell TS, Luster AD: Inhibition of pulmonary fibrosis by the chemokine IP-10/CXCL10. Am J Respir Cell Mol Biol. 2004, 31: 395-404. 10.1165/rcmb.2004-0175OC.
Walters DM, Antao-Menezes A, Ingram JL, Rice AB, Nyska A, Tani Y, Kleeberger SR, Bonner JC: Susceptibility of signal transducer and activator of transcription-1-deficient mice to pulmonary fibrogenesis. Am J Pathol. 2005, 167: 1221-1229. 10.1016/S0002-9440(10)61210-2.
Spagnolo P, Del GC, Luppi F, Cerri S, Balduzzi S, Walters EH, D’Amico R, Richeldi L: Non-steroid agents for idiopathic pulmonary fibrosis. Cochrane Database Syst Rev. 2010, 9: CD003134-
Black CM, Silman AJ, Herrick AI, Denton CP, Wilson H, Newman J, Pompon L, Shi-wen X: Interferon-alpha does not improve outcome at one year in patients with diffuse cutaneous scleroderma: results of a randomized, double-blind, placebo-controlled trial. Arthritis Rheum. 1999, 42: 299-305. 10.1002/1529-0131(199902)42:2<299::AID-ANR12>3.0.CO;2-R.
Navab R, Strumpf D, Bandarchi B, Zhu CQ, Pintilie M, Ramnarine VR, Ibrahimov E, Radulovich N, Leung L, Barczyk M, Panchal D, To C, Yun JJ, Der S, Shepherd FA, Jurisica I, Tsao MS: Prognostic gene-expression signature of carcinoma-associated fibroblasts in non-small cell lung cancer. Proc Natl Acad Sci U S A. 2011, 108: 7160-7165. 10.1073/pnas.1014506108.
Katzenstein AL, Mukhopadhyay S, Zanardi C, Dexter E: Clinically occult interstitial fibrosis in smokers: classification and significance of a surprisingly common finding in lobectomy specimens. Hum Pathol. 2010, 41: 316-325. 10.1016/j.humpath.2009.09.003.
Reiser J, Hurst J, Voges M, Krauss P, Munch P, Iftner T, Stubenrauch F: High-risk human papillomaviruses repress constitutive kappa interferon transcription via E6 to prevent pathogen recognition receptor and antiviral-gene expression. J Virol. 2011, 85: 11372-11380. 10.1128/JVI.05279-11.
Bidwell BN, Slaney CY, Withana NP, Forster S, Cao Y, Loi S, Andrews D, Mikeska T, Mangan NE, Samarajiwa SA, de Weerd NA, Gould J, Argani P, Moller A, Smyth MJ, Anderson RL, Hertzog PJ, Parker BS: Silencing of Irf7 pathways in breast cancer cells promotes bone metastasis through immune escape. Nat Med. 2012, 18: 1224-1231. 10.1038/nm.2830.
Kulaeva OI, Draghici S, Tang L, Kraniak JM, Land SJ, Tainsky MA: Epigenetic silencing of multiple interferon pathway genes after cellular immortalization. Oncogene. 2003, 22: 4118-4127. 10.1038/sj.onc.1206594.
Kroger A, Dallugge A, Kirchhoff S, Hauser H: IRF-1 reverts the transformed phenotype of oncogenically transformed cells in vitro and in vivo. Oncogene. 2003, 22: 1045-1056. 10.1038/sj.onc.1206260.
Gao J, Wang Y, Xing Q, Yan J, Senthil M, Akmal Y, Kowolik CM, Kang J, Lu DM, Zhao M, Lin Z, Cheng CH, Yip ML, Yim JH: Identification of a natural compound by cell-based screening that enhances interferon regulatory factor-1 activity and causes tumor suppression. Mol Cancer Ther. 2011, 10: 1774-1783. 10.1158/1535-7163.MCT-11-0304.
Vancheri C: Idiopathic pulmonary fibrosis: an altered fibroblast proliferation linked to cancer biology. Proc Am Thorac Soc. 2012, 9: 153-157. 10.1513/pats.201203-025AW.
This work was supported by the Royal Brompton & Harefield Clinical Research Fund (RBHCRF), the Raynaud’s and Scleroderma Association (RSA), the Scleroderma Society, and by the Arthritis Research Council (now AR UK). EAR was funded by a project grant from RSA. CS and PL were funded by a program grant from RSA. The study was also supported by the Respiratory Biomedical Research Unit at Royal Brompton & Harefield NHS Foundation Trust. We would like to thank the Microarray Centre, Genomics Laboratory (CSC, MRC, Imperial College London and Hammersmith Hospital) for support with the microarray analysis.
The authors declare that they have no competing interests.
GEL CJWS contributed to design, data generation, and data analysis of the study, and wrote and prepared the manuscript. XS performed all explant cell preparations, cell culture experiments and western blot analyses. PS contributed to data analysis and interpretation, and critical review of the manuscript. PL SLH contributed to the verification of the microarray data and valuable discussions. GG participated in the microarray work. AGN contributed human lung tissue samples and expert histopathological typing. CPD and JCG contributed patient material and information. DJA, CPD AUW TMM contributed to design and critical supervision of the study. EAR designed the study, performed experimental work and wrote the manuscript. All authors contributed to critical revision of the manuscript. All authors read and approved the final manuscript.
Gisela E Lindahl, Carmel JW Stock contributed equally to this work.