Early events marking lung fibroblast transition to profibrotic state in idiopathic pulmonary fibrosis
Respiratory Research volume 24, Article number: 116 (2023)
Idiopathic Pulmonary Fibrosis (IPF) is an age-associated progressive lung disease with accumulation of scar tissue impairing gas exchange. Previous high-throughput studies elucidated the role of cellular heterogeneity and molecular pathways in advanced disease. However, critical pathogenic pathways occurring in the transition of fibroblasts from normal to profibrotic have been largely overlooked.
We used single cell transcriptomics (scRNA-seq) from lungs of healthy controls and IPF patients (lower and upper lobes). We identified fibroblast subclusters, genes and pathways associated with early disease. Immunofluorescence assays validated the role of MOXD1 early in fibrosis.
We identified four distinct fibroblast subgroups, including one marking the normal-to-profibrotic state transition. Our results show for the first time that global downregulation of ribosomal proteins and significant upregulation of the majority of copper-binding proteins, including MOXD1, mark the IPF transition. We find no significant differences in gene expression in IPF upper and lower lobe samples, which were selected to have low and high degree of fibrosis, respectively.
Early events during IPF onset in fibroblasts include dysregulation of ribosomal and copper-binding proteins. Fibroblasts in early stage IPF may have already acquired a profibrotic phenotype while hallmarks of advanced disease, including fibroblast foci and honeycomb formation, are still not evident. The new transitional fibroblasts we discover could prove very important for studying the role of fibroblast plasticity in disease progression and help develop early diagnosis tools and therapeutic interventions targeting earlier disease states.
Idiopathic Pulmonary Fibrosis (IPF) is an age-dependent chronic lung disease affecting individuals generally over 60 years old . The mechanisms driving the disease development and progression are still not fully understood . Gene-by-gene analysis and high-throughput studies have promoted the field over the years, offering valuable insights into the pathophysiology of the disease [3,4,5]. Similar to other diseases, single cell approaches have the potential to advance this knowledge even further by fully dissecting the disease milieu and allowing investigators to assess individualized cell contributions and complex cellular dynamics during disease emergence and progression [6,7,8].
Extensive research has provided valuable insight into the origin and contribution of fibroblasts and especially myofibroblasts in pulmonary fibrosis [9, 10]. All evidence supports the hypothesis that these are heterogeneous groups of cells, undergoing very distinct transition processes dictated by the disease microenvironments [9,10,11]. They share common characteristics and express a set of biomolecules that drive fibrosis . Single cell transcriptomics allow dissection of this heterogeneous population. Murine models have been used to better understand how fibroblast subtypes contribute to fibrosis . More recently, single cell RNA sequencing (scRNA-seq) of fresh human tissue revealed cell-specific differences between normal and fibrotic tissue [8, 13].
To study the disease onset and progression, our group has performed scRNA-seq on upper and lower lobes of fresh human explanted IPF lungs and on healthy controls, identifying a subpopulation of proliferating SPP1Hi macrophages with a potential role in lung fibrosis . Here, we reanalyze the raw sequencing data after performing an imputation step and focus on the early events that could drive transition of normal lung fibroblasts to profibrotic. Our analyses reveal four major fibroblast clusters showing unique characteristics regarding their gene expression and related to their tissue of origin (control or IPF). A closer look at each cluster confirms previously reported changes associated with IPF and—more importantly—reveals a profibrotic state of the upper (unaffected) lung. Gene expression patterns in each cluster reveal an expected dysregulation of fibrosis-associated genes. A novel finding is the dysregulation of genes coding for copper-binding proteins during both early and late stages of the disease. Immunohistochemistry assays in IPF lung fibroblasts show high levels of one of the top differentially expressed copper-binding genes, MOXD1. Furthermore, pseudotime analysis identifies a distinct group of fibroblasts in the process of acquiring a profibrotic phenotype while undergoing a global downregulation of genes coding for ribosomal proteins. Both copper-binding  and ribosomal protein pathways [15,16,17] have been previously associated to senescence. The discovery of this new type of transitional fibroblasts is important as it provides insights into the early events of fibrosis and may prove more helpful in finding suitable targets for early diagnosis and possible therapeutic interventions.
Single cell RNA sequencing (scRNA-seq)
Raw scRNA-seq data were derived from our previous publication . Briefly, normal and IPF lung tissue (Additional file 1: Table S1, Additional file 2: Fig. S1) was obtained and processed to obtain single cell suspension. Single cell libraries were prepared using the 10X Genomics Chromium instrument and V2 chemistry. Sequencing was performed on an Illumina NextSeq-500 instrument.
scRNA-seq data analysis
scRNA-seq raw count and cell-UMI (Unique Molecular Identifier) count matrix were generated by Cell Ranger . Single-cell Analysis Via Expression Recovery (Saver) was used to impute dropout events in gene expression . Seurat (version 2.3.4) was used to normalize gene expression, perform differentially expressed gene analysis, identify distinct cell populations and visualize clusters graphically [19, 20]. The cell-UMI matrix was filtered and only cells expressing at least 200 genes were further analyzed. Cells containing greater than 35% of mitochondrial genes were also excluded from the analysis. Highly variable genes were identified, based on their average expression and dispersion, and would be used in the downstream analysis. Data were scaled and the number of UMIs per cell as well as the percentage of mitochondrial gene content were regressed out. In this study, we further removed the effect from technical or biological confounders using Harmony, which integrates multiple scRNA seq datasets by projecting cells into a shared embedding in which cells are grouped by cell type, not the specific conditions related to the datasets . t-Distributed stochastic neighbor embedding (t-SNE) plots based on Harmony embeddings were generated to assign clusters and each cluster was identified by differentially expressed genomic signatures.
Analysis of fibroblast cluster from scRNA-seq data
The fibroblast cells from scRNA-seq data were processed by Seurat and Harmony to identify subclusters of fibroblast cells. Destiny was used for pseudotime analysis of the subclusters . Velocyto was used to estimate the time derivative of the gene expression state . RNA velocity was estimated using gene-relative model with k-nearest neighbor cell pooling (k = 20) and velocity fields were projected into a UMAP-based embedding through SeuratWrappers and Seurat (version 3.1.0). We also used AddModuleScore of Seurat to calculate average expression level of (i) copper binding , (ii) senescence  and (iii) ribosome biogenesis (Gene Ontology Browser) gene sets on a single cell level, normalized by randomly selected control feature set.
Fresh lung tissue (Additional file 1: Table S1) was fixed with 4% paraformaldehyde and embedded in OCT (Tissue-Tek® Sakura® Finetek, US). Cultured fibroblasts (Additional file 1: Table S1) were plated in chamber-slides and fixed with 2% paraformaldehyde. OCT sections (5 μM thick) and fixed fibroblasts were stained with the primary/secondary antibodies listed in Additional file 1: Table S2. Slides were mounted using ProLong Gold Antifade Mountant with DAPI (Life technologies, USA). Images were obtained using an Olympus Fluoview 1000-3 Confocal Microscope (20× objective).
Briefly, formalin-fixed paraffin-embebbed tissue sections were deparaffinised, rehydrated and subjected to an antigen retrieval step, blocked with 10% Normal Serum (PK-6200, Vector Laboratories, Burlingame, USA) in DAKO diluent (S3022, DAKO, Glostrup, Denmark), washed and incubated overnight at 4 °C with 1:300 of anti-MOXD1 (bs-17733R, Bioss, Massachusetts, USA) antibody. Specific staining was detected with secondary biotinylated horse anti-rabbit antibody (PK-6200, Vector Laboratories, Burlingame, USA) followed by Vectorstain ABC Elite reagent and DAB peroxidase substrate (SK-4105, Vector laboratories, Burlingame, USA) and counterstained with haematoxylin. Slides were mounted with permount mounting medium (SP15-100, Fisher Scientific, Pittsburgh, PA, USA) and analyzed in a light microscopy (Nikon Eclipse 55i, Melville, NY, USA).
Cell types in IPF and healthy lungs
A mixed population of cells (48,023 cells) from the lungs of three healthy control (21,485 cells) and six IPF (26,538 cells) samples (from three donors) was sequenced using next generation single cell sequencing. An important study design includes tissue from three upper lobes (IPF_UP; 12,577 cells) and fibrotic tissue from three lower lobes (IPF_LOW; 13,961 cells) (Fig. 1A). In general, samples from the upper lobes presented a relatively unaffected lung histology, while all samples from the lower lobes showed advanced fibrosis . Here, we performed imputation of dropouts and gene expression profiles were clustered and visualized (t-SNE). Our analysis resulted in 20 clusters (Fig. 1B) that were identified as distinct cell types (Fig. 1C), using previously described cell markers . Each cluster had variable cell counts coming from healthy and IPF upper and lower lung lobes (Fig. 1D, E and Additional file 2: Fig. S1).
Fibroblast (FB) clusters’ characteristics
Fibroblasts (Fig. 1C, green box) were re-clustered and visualized in two dimensions using t-SNE (Fig. 2A). Overall, the clustering of fibroblasts from 449 control, 404 IPF_UP and 1085 IPF_LOW samples revealed four different groups. A fifth cluster, neighboring B-cells, was found to contain cells expressing HLA type II family genes and was excluded from this analysis (Additional file 2: Fig. S2). Despite our best efforts, one of the IPF upper lung lobe samples (SC154) showed extensive fibrosis compared to normal controls (Additional file 2: Fig. S1A), posing a risk of masking any differences between the IPF_UP and IPF_LOW fibroblasts. We decided to include this sample in our analysis since it clustered similarly to the other fibroblasts in IPF_UP samples (Additional file 2: Fig. S3A). Similarly, we checked the age distribution of the fibroblast cells. Despite the fact that IPF samples generally come from older people (since IPF usually presents itself > 55 yrs), we did not see any other age bias (Additional file 2: Fig. S3B).
Cluster-1 (608 cells, 31% of all FB) contained mainly cells from the healthy control donors (80% of the control FBs were in Cluster-1) as well as a smaller percentage of the IPF_UP and IPF_LOW cells (24% and 14% respectively). Cluster-2, a much smaller cluster (197 cells, 10% of all FB), included cells from each sample origin. Cluster-3 (920 cells, 48% of all FB) contained mainly IPF_UP and IPF_LOW (56% and 64% respectively). Only 1% of the healthy control FBs classified as Cluster-3. Cluster-4 (213 cells, 11% of all FB) contained a slightly higher percentage of the IPF_UP and IPF_LOW (12% and 13% respectively) than the healthy control fibroblasts (6%). The breakdown of each cluster to sample types is summarized in Fig. 2B, C while the breakdown of each sample to different clusters is summarized in Fig. 2D.
Since fibroblasts in Cluster-1 came predominantly from control samples we consider this to represent the “normal state”; while those in Cluster-3, which came almost exclusively form IPF_UP and IPF_LOW, represent the “disease state”. Cluster-1, Cluster-2 and Cluster-4 had distinct expression profiles, but they were all present in healthy and IPF lungs (Fig. 2C). Overall, IPF_UP and IPF_LOW fibroblasts showed similar clustering and compositional patterns (Fig. 2B, D). These results showed that fibroblasts from the relatively non-fibrotic upper lobes had similar molecular signatures to those from the highly fibrotic lower lobes, suggesting that molecular changes in lung fibroblasts precede the morphological changes identified by histological examination.
We also noticed that both Cluster-1 and CLuster-3 contained some alveolar fibroblasts (SPINT2+FGFR4+GPC3+), while cluster 2 was predominately composed of adventitial fibroblasts (SFRP2+PI16+SERPINF1+) (Additional file 2: Fig. S4), which were previously reported . Additionally, some fibroblast cells in Cluster-3 highly expressed LGR5, and LGR5+ fibroblasts serve as a mesenchymal niche for distal airway epithelial cells, as it was recently reported . From the UMAP plots, LGR5+ fibroblasts are next adventitial fibroblasts (Additional file 2: Fig. S4A). Furthermore, part of Cluster-3 and Cluster-4 are CTHRC1+ fibroblasts, which are uniquely found in fibrotic lungs (Additional file 2: Fig. S4B). In a previous RNA velocity analysis in mouse fibroblasts, vectors from alveolar fibroblasts and adventitial fibroblasts pointed to CTHRC1+ fibroblasts . Consistently, our velocity analysis showed that vectors from myofibroblast Cluster-4, alveolar fibroblasts, and LGR5+ fibroblasts, point towards CTHRC1+ fibroblasts.
Differentially expressed genes in each cluster are listed in Additional file 1: Table S3 and expression of the top 16 genes in these lists is further described per cluster and sample type (Figs. 3, 4, 5, 6). Briefly, fibroblasts in Cluster-1 had higher expression in 263 genes (FDR < 0.05) with the highest being: ICAM1, CXCL1, CXCL2, CXCL3, CXCL8, CCL2, IL6, PTX3, IER3, GADD45B, THBS1, SOD2 and NFKBIA which are associated with immune response and regulation of inflammation (Fig. 3 and Additional file 1: Table S4).
Among the 118 genes that were significantly upregulated in Cluster-2, CXCL14, SFRP2 and SFRP4 had higher expression levels in IPF samples. Interestingly, the highest levels of CXCL14 were detected in IPF_LOW fibroblasts, being one of the most significant differences observed overall between IPF_UP and IPF_LOW samples. Profibrotic gene CXCL14 is shown to play a regulatory role in immune response and inflammation while SFRP2 and SFRP4 are regulators (inhibitors) of Wnt signaling (Fig. 4 and Additional file 1: Table S4).
The majority of the top 16 upregulated genes in Cluster-3 (out of the total 219 differentially expressed) have been previously shown to increase in lung fibrosis (Fig. 5 and Additional file 1: Table S4). Our results corfirmed that cells expressing myofibroblasts markers like POSTN and ASPN also expressed LTBP2, LTBP1, BGN, DPT and the highest levels of COL3A1 and COL8A1. Interestingly, 68–93% of the cells on this cluster that had low expression of myofibroblast-associated genes also had high expression levels of IPF-associated profibrotic markers such as MFAP4, LTBP1, BGN, COMP, MMP2, COL3A1 and COL8A1. All 16 genes had similar expression in IPF_UP and IPF_LOW samples. CXCL14, which was upregulated in Cluster-2, was also upregulated in the Cluster-3 IPF_LOW compared to IPF_UP.
Cells in Cluster-4 expressed profibrotic and myofibroblast markers like TNC, collagens, SPARC, POSTN, FN1 and TPM2 in higher levels than cells in Cluster-3. All 16 top upregulated genes in this group had similar expression levels in IPF_UP and IPF_LOW samples (Fig. 6 and Additional file 1: Table S4).
Full names and relevant references for all genes mentioned above are listed in Additional file 1: Table S4.
Trajectory and RNA velocity analysis reveals early disease-associated events
Single cell trajectory inference and pseudotime analysis are often used to study cellular dynamics or transitional states and special organization of cells in tissues . Each cell is assigned a numeric value (pseudotime) which indicates where in the underlying dynamic biological process that cell falls into. These methods allow the visualization of intermediate stages that connect distinct cell states, which are often overlooked during cell clustering. We used diffusion map (destiny) to perform pseudotime analysis on the four fibroblast subclusters (Fig. 7A). The diffusion components, as measures of pseudotime , revealed temporal ordering and cellular decision on the single-cell transcriptome level. Interestingly, cells in Cluster-1 had two distinct subpopulations (Fig. 7A, ovals) which corresponded to cells from Control and IPF samples with the latter being closer to Cluster-2 and 3 (Fig. 7B). This suggests that IPF fibroblasts in Cluster-1, despite their overall “healthy” phenotype, have been influenced by the disease environment.
To further our understanding of the mechanisms driving disease onset, we calculated RNA velocity , which is based on the balance of unspliced (nascent) and spliced (mature) mRNA. This high-dimensional vector can act as an indicator of the future state of mature mRNA, driving cell state. The directionality of the plot arrows depicts the direction of cell state progression. RNA velocity analysis of all fibroblasts (Fig. 7C) showed that the majority of the fibroblasts in Cluster-1 (composed mainly of non-IPF fibroblasts) trended towards a steady state (SS1). Even fibroblasts originating from IPF samples in this cluster appeared to be able to express genes driving them towards SS1. Similarly, cluster-3 fibroblasts are attracted to another steady state (SS2) which is likely disease related since almost all of these cells originated from IPF lungs. Although not as clear, Cluster-4 also showed a directional flow towards SS2 suggesting that these cells may transition to a profibrotic phenotype but have a cellular transcriptome that differentiates them from fibroblasts in Cluster-3 (Fig. 7C). TIMP1, COL1A1, CTHRC1, TUBA1B, SH3BGRL3 and S100A16 were the top six differentially expressed genes between Cluster-3 and Cluster-4 (upregulated in Cluster-4).
Cluster-2 fibroblasts showed a more complex pattern (Fig. 7C) with some cells going towards Cluster-1 (subgroup C2G1; Fig. 7D, gold) and others towards Cluster-3 (subgroup C2G2; Fig. 7D, purple) suggesting that cells in Cluster-2 were captured in a transition state that could either be resolved (SS1) or progress to a profibrotic state (SS2). We found 74 of the 103 genes coding for ribosomal proteins to be significantly downregulated in C2G2 vs C2G1 (Additional file 1: Table S5). The remaining 29 ribosomal protein genes were unaffected. Based on these 74 genes, C2G1 and C2G2 cells were assigned a ribosome score which was significantly decreased in C2G2 (p-value ≪ 0.001). In addition, we calculated a ribosome biogenesis score based on 267 other genes associated with this process (Gene Ontology Browser). Again, C2G2 showed a significant decrease for this score (p-value ≪ 0.001). Both these observations (Fig. 7E) are consistent with the hypothesis that there is a dysregulation of ribosomal function at the early stages of IPF.
To validate this finding, we used a second publicly available scRNA-seq dataset , that had similar tissue collection and processing procedures as ours. This dataset consisted of eight transplant donors and four IPF ex-plants, but it did not have distinct upper/lower lobe samples. We processed and reanalyzed the raw scRNA-seq data in the same way as our dataset. Due to the low number of IPF fibroblasts (see Additional file 2: Fig. S5B) we couldn’t confidently identify a subcluster similar to our Cluster-2 in the fibroblast cluster, but we observed three areas containing distinct cell populations originating from Controls, IPF and a mixture of the two, respectively (Additional file 2: Fig. S5A, green). Analysis of those cells showed 32 ribosomal genes to be differentially expressed between control and IPF with 29 of them downregulated. Notably, all 29 were also included in the previous list of 74 differentially expressed genes in our Cluster-2 (p-value < 0.01).
Monooxygenase DBH like 1 (MOXD1): a novel early IPF biomarker
To identify early IPF biomarkers specific to fibroblasts, we examined the Cluster-3 top upregulated genes that also showed low expression in the rest of the cell types (Fig. 5). The five most significant of those included known IPF players (metalloprotease, collagen and two -LTBP1 LTBP2) and MOXD1 (Additional file 1: Table S3). MOXD1 was not previously associated with IPF, but it is expressed significantly higher in Cluster 3 and it is significantly downregulated in controls (Fig. 8A–C). Furthermore, it had been described in replicative senescent fibroblasts . Given that IPF fibroblasts can be resistant to apoptosis  and have a senescent phenotype , we selected MOXD1 for validation as a new possible fibroblast biomarker in early and advanced IPF.
Immunofluorescence and immunohistochemistry staining of human lung tissues (Fig. 8D, Additional file 2: Figs. S6 and S7) validated the mRNA findings. Cells from healthy donors’ tissues expressed minimal MOXD1 protein. IPF lung upper lobes, showing less advanced fibrotic disease, also showed increased MOXD1 compared to control lungs. IPF lower lobes, showing advanced fibrotic disease, expressed even higher levels of MOXD1. We then proceed to characterize the fibroblast cells expressing MOXD1 by co-staining the sections with PDGFRa (marker of fibroblasts) and alpha smooth muscle actin (α-SMA, a marker of myofibroblasts, which also marks pericytes and smooth muscle cells). As expected, α-SMA staining was increased in upper lobes and even more highly increased in lower lobes of IPF lungs compared to control lungs. Most MOXD1 staining in upper lobes was found in cells not staining with α-SMA, while these markers co-stained many of the same cells in the lower lobe suggesting that MOXD1 is a marker of cells that are maturing into myofibroblasts. This is consistent with increased levels of myofibroblasts in IPF. PDGFRa staining showed that fibroblasts in control lungs are mostly MOXD1/SMA negative, while fibroblasts in lower lobes of IPF lungs are mostly MOXD1/SMA positive, co-staining with these markers indicated that most of the MOXD1 staining cells in IPF are myofibroblasts. To assess location of MOXD1 fibroblasts, we performed immunostaining of whole lung tissues. Lung sections from age matched donor control were negative for MOXD1, whereas sections from IPF lungs showed numerous positive cells (Fig. 8D). Immunohistochemistry staining also showed increased MOXD1 expression in fibroblast foci (Additional file 2: Fig. S7).
Velocity analysis of fibroblasts using their spliced (Fig. 8E, left) and unspliced (Fig. 8E, center) MOXD1 mRNA context indicated transcriptional kinetics towards the intersection between Cluster-3 and Cluster-4. A closer inspection of this region showed high expression of the myofibroblast associated genes POSTN, ASPN, ACTA2 and TAGLN (Fig. 8G) [11, 33]. It was clear from the UMAP plots that the cells carrying these myofibroblast markers did not segregate in one distinct cluster (Fig. 7C) suggesting that myofibroblast specific genes are having similar signatures to the original cells that transitioned to a myofibroblast.
Copper-binding proteins are upregulated early in disease
MOXD1 is a copper-binding protein and copper itself has been implicated in pulmonary fibrosis . Blockhuys et al.  have identified 54 human Cu-binding proteins, divided in 9 groups according to their cellular localization. We found 39 of them to be significantly up- and 10 down-regulated in Cluster-3 (“fibrotic”) vs Cluster-1 (“control”) (FDR < 0.05; Additional file 1: Table S6A). Furthermore, 20 of them were already significantly upregulated in C2G2 compared to C2G1 subcluster (Additional file 1: Table S6B). Among the top differentially expressed genes we found SPARC, LOXL1, and APP which have been implicated in lung fibrosis [35,36,37]. In similar analysis to ribosome genes, we found 29 genes to be significantly upregulated in the validation cohort; and 24 of them were included in the 39 of Cluster-3. Next, we assigned a copper-binding score to all fibroblasts based on the expression of these proteins in our dataset. IPF samples, regardless of lung topology, showed higher scores compared to the control samples (Fig. 8F, left). This difference was also reflected in the fibroblast clusters where Cluster-1 cells had the lowest Cu-binding score and Cluster-3 and 4 had the highest. Cluster-2, consistent with an intermediate cell state, showed an intermediate score (Fig. 8F, right).
Fibroblasts have been previously studied for their role in fibrosis in general and lung fibrosis in particular . Research has been focused on the lower subpleural lung regions where fibrotic scar tissue is primarily localized, assuming that fibroblasts from these regions contribute the most to disease development and progression. In the present work we are extending research to generally less disease-affected areas, the upper lobes of the lungs, to investigate early events in disease onset.
Our results show that the majority of IPF lung fibroblasts have a distinct molecular phenotype, since they form a distinct cluster. Interestingly, fibroblasts from the upper and lower lobes have similar molecular profiles, which is consistent with the hypothesis that gene expression changes in fibroblasts happen at an early disease stage. The observed early shift in gene expression in IPF fibroblasts, especially in extracellular matrix associated genes (collagens, MMP2, etc.), agrees with computationally deconvoluted bulk RNA gene expression signatures . These models indicated that deciding regulatory molecular events can happen early or late in lung fibrosis, depending on the cell type. Our scRNA-seq results confirm this prediction. However, we note as a potential limitation of this analysis that the control samples come from younger individuals compared to the IPF samples.
Our results indicate that underlying cell state differences (specifically ribosome and copper binding proteins) influence the differentiation of fibroblasts to myofibroblasts. One of the top differentially expressed genes in the IPF fibroblasts (Cluster-3) is MOXD1, a copper-binding enzyme, which has not been previously studied in IPF. MOXD1 shows a tight association with the ER membrane and is not secreted . Currently its substrate is unknown. It is found to be upregulated in senescent human fibroblasts and human vascular endothelial cells (HUVEC) . Most importantly, IPF fibroblasts, under both in vitro and in vivo conditions and regardless of their lung localization, show high expression of MOXD1. Our velocity analysis results of MOXD1 nascent and mature mRNAs suggests that it start expressing earlier than known myofibroblast markers and could be potentially used for early diagnosis.
We find that MOXD1 upregulation in fibrotic cells is part of the more general upregulation of copper binding enzymes, including Lysyl Oxidase Like 1 (LOXL1) and Secreted Protein Acidic and Cystein Rich (SPARC), which have been studied for their pro-fibrotic role in multiple organs [37, 40]. The dysregulation of copper-binding proteins along with the recent demonstration of copper induced lung fibrosis in mice  and copper-induced premature senescence of human fibroblasts  supports a more comprehensive role of copper and fibroblast senescence in the pathogenesis of IPF. It would be interesting to see whether IPF fibroblasts accumulate intracellular copper similar to the progressive accumulation of intracellular copper in aging, kidney fibrosis and other diseases [41,42,43].
Depending on their localization (IPF_LOW, IPF_UP), approximately 15% to 25% of fibroblasts from IPF samples retain a transcriptional profile that resembles that of normal control lung fibroblasts. However, they form a subgroup which indicates that their “normal” phenotype is influenced by the disease environment. Another 25% of IPF fibroblasts along with 15% of control fibroblasts seems to be in “perturbed” states that are characterized by increased expression of ECM genes. It is reasonable to assume that these states may represent normal wound healing processes occurring in both control and IPF lungs. In the case of IPF lungs, these processes would eventually run uncontrolled resulting in severe fibrosis .
One particular group of cells (Cluster-2) seem to capture this transition. It shows two distinct cell sub-populations with tendency towards opposing states. One subgroup shows a tendency to revert to a normal state while the other seems to be transitioning towards a profibrotic state. A closer examination of the profibrotic sub-population reveals a downregulation of almost all ribosomal proteins and dysfunctional ribosomal biosynthesis; and upregulation of the majority of copper-binding genes. We validate these findings in a similar cohort from a previous study . Posttranslational regulation has been shown to play an important role in heart fibrosis through the largely underexplored roles of RNA binding proteins and variation in ribosome occupancy, which affect protein expression levels independent of mRNA levels . Importantly, dysregulation of translational control on the level of polyribosome formation has been shown to play a significant role in the emergence of IPF myofibroblasts . Although we observe a decrease in ribosomal gene expression, we wouldn’t necessarily associate it with decreased genome-wide translational activity in these cells. Disproportionate expression of ribosomal proteins could affect ribosomal specificity [47, 48], favoring the translation of a profibrotic proteome in fibroblasts. Moreover, quantitative change in ribosome biogenesis could increase competition for ribosome binding and translation initiation, causing a variable effect on the translation of cellular mRNAs .
Impairment of ribosomal biogenesis not only could have an impact on the ribosomes’ cellular housekeeping role in protein synthesis but could also affect cell cycle and proliferation [15, 50]. Ribosome biogenesis can be impaired at multiple steps by a variety of stress factors, resulting in ribosomal stress and cellular senescence [15,16,17]. Our results support the hypothesis that ribosomal impairment may be a general mechanism activated during switching between two cell fates. In the future, experimental manipulation of ribosomal proteins will be necessary to decipher whether ribosome dysregulation is a driver of IPF or one of the many phenotypes associated with the dysregulated homeostasis caused by the disease.
In summary, we examine the events that take place in fibroblasts when cells transition from healthy state to disease. We demonstrate that, at the single cell level, changes in the transcriptome of fibroblasts appear early in the relatively unaffected upper lobes of the IPF lungs. We also identify two focal points (steady states) towards where the normal and the fibrotic cells coalesce. More importantly, we discover a point of transition between these two states, which is characterized by a significant decrease in ribosomal protein genes and a significant increase in copper-binding proteins. We examine the second most upregulated copper-binding protein, MOXD1, and find it to be fibroblast- and IPF-specific. Furthermore, we find that MOXD1 is also expressed at early stages of IPF suggesting that it can be used as an early biomarker of fibroblast pulmonary fibrosis.
Increasing evidence supports the idea that IPF fibroblasts have acquired resistance to apoptosis  and have a senescent phenotype . Furthermore, senescence can mediate the fibrotic phenotype , perhaps by secreting profibrotic factors . Thus, early targeting of senescent fibroblasts with senolytic cocktails may help manage or even reverse the disease [51, 52]. To that extend, in future studies the transitional fibroblasts we discovered can farther our mechanistic understanding of the transition process and can be utilized for the development of potential diagnostic and therapeutic strategies.
At present, a major limitation of this study is the lack of available biomarkers for these transitioning fibroblasts. Future studies can reveal how senescence and alteration of ribosomal protein and copper-binding pathways may affect the switch of fibroblasts from the normal steady state to a profibrotic state. Furthermore, in-depth studies of these fibroblasts may provide a useful diagnostic tool for early detection of IPF and possibly a therapeutic window before extensive tissue damage occurs.
Availability of data and materials
All data presented in this manuscript are available from GEO (GSE128033, GSE122960).
Wolters PJ, Collard HR, Jones KD. Pathogenesis of idiopathic pulmonary fibrosis. Annu Rev Pathol. 2014;9:157–79.
du Bois RM. Strategies for treating idiopathic pulmonary fibrosis. Nat Rev Drug Discov. 2010;9:129–40.
Vukmirovic M, Kaminski N. Impact of transcriptomics on our understanding of pulmonary fibrosis. Front Med (Lausanne). 2018;5:87–87.
Habermann AC, Gutierrez AJ, Bui LT, Yahn SL, Winters NI, Calvi CL, Peter L, Chung M-I, Taylor CJ, Jetter C, Raju L, Roberson J, Ding G, Wood L, Sucre JM, Richmond BW, Serezani AP, McDonnell WJ, Mallal SB, Bacchetta MJ, Loyd JE, Shaver CM, Ware LB, Bremner R, Walia R, Blackwell TS, Banovich NE, Kropski JA. Single-cell RNA-sequencing reveals profibrotic roles of distinct epithelial and mesenchymal lineages in pulmonary fibrosis. Sci Adv. 2020; 6:eaba1972. https://doi.org/10.1126/sciadv.aba1972.
McDonough JE, Ahangari F, Li Q, Jain S, Verleden SE, Herazo-Maya J, Vukmirovic M, DeIuliis G, Tzouvelekis A, Tanabe N, Chu F, Yan X, Verschakelen J, Homer RJ, Manatakis DV, Zhang J, Ding J, Maes K, De Sadeleer L, Vos R, Neyrinck A, Benos PV, Bar-Joseph Z, Tantin D, Hogg JC, Vanaudenaerde BM, Wuyts WA, Kaminski N. Transcriptional regulatory model of fibrosis progression in the human lung. JCI Insight. 2019;4: e131597.
Reyfman PA, Walter JM, Joshi N, Anekalla KR, McQuattie-Pimentel AC, Chiu S, Fernandez R, Akbarpour M, Chen C-I, Ren Z, Verma R, Abdala-Valencia H, Nam K, Chi M, Han S, Gonzalez-Gonzalez FJ, Soberanes S, Watanabe S, Williams KJN, Flozak AS, Nicholson TT, Morgan VK, Winter DR, Hinchcliff M, Hrusch CL, Guzy RD, Bonham CA, Sperling AI, Bag R, Hamanaka RB, Mutlu GM, Yeldandi AV, Marshall SA, Shilatifard A, Amaral LAN, Perlman H, Sznajder JI, Argento AC, Gillespie CT, Dematte J, Jain M, Singer BD, Ridge KM, Lam AP, Bharat A, Bhorade SM, Gottardi CJ, Budinger GRS, Misharin AV. Single-cell transcriptomic analysis of human lung provides insights into the pathobiology of pulmonary fibrosis. Am J Respir Crit Care Med. 2019;199:1517–36.
Adams TS, Schupp JC, Poli S, Ayaub EA, Neumark N, Ahangari F, Chu SG, Raby BA, DeIuliis G, Januszyk M, Duan Q, Arnett HA, Siddiqui A, Washko GR, Homer R, Yan X, Rosas IO, Kaminski N. Single Cell RNA-seq reveals ectopic and aberrant lung resident cell populations in idiopathic pulmonary fibrosis. Sci Adv. 2020;6:eaba1983. https://doi.org/10.1126/sciadv.aba1983.
Morse C, Tabib T, Sembrat J, Buschur KL, Bittar HT, Valenzi E, Jiang Y, Kass DJ, Gibson K, Chen W, Mora A, Benos PV, Rojas M, Lafyatis R. Proliferating SPP1/MERTK-expressing macrophages in idiopathic pulmonary fibrosis. Eur Respir J. 2019;54:1802441.
Xie T, Wang Y, Deng N, Huang G, Taghavifar F, Geng Y, Liu N, Kulur V, Yao C, Chen P, Liu Z, Stripp B, Tang J, Liang J, Noble PW, Jiang D. Single-cell deconvolution of fibroblast heterogeneity in mouse pulmonary fibrosis. Cell Rep. 2018;22:3625–40.
Habiel DM, Hogaboam CM. Heterogeneity of fibroblasts and myofibroblasts in pulmonary fibrosis. Curr Pathobiol Rep. 2017;5:101–10.
Guerrero-Juarez CF, Dedhia PH, Jin S, Ruiz-Vega R, Ma D, Liu Y, Yamaga K, Shestova O, Gay DL, Yang Z, Kessenbrock K, Nie Q, Pear WS, Cotsarelis G, Plikus MV. Single-cell analysis reveals fibroblast heterogeneity and myeloid-derived adipocyte progenitors in murine skin wounds. Nat Commun. 2019;10:650.
Kendall RT, Feghali-Bostwick CA. Fibroblasts in fibrosis: novel roles and mediators. Front Pharmacol. 2014;5:123.
Valenzi E, Bulik M, Tabib T, Morse C, Sembrat J, Trejo Bittar H, Rojas M, Lafyatis R. Single-cell analysis reveals fibroblast heterogeneity and myofibroblasts in systemic sclerosis-associated interstitial lung disease. Ann Rheum Dis. 2019;78:1379–87.
Matos L, Gouveia A, Almeida H. Copper ability to induce premature senescence in human fibroblasts. Age (Dordr). 2012;34:783–94.
Zhou X, Liao W-J, Liao J-M, Liao P, Lu H. Ribosomal proteins: functions beyond the ribosome. J Mol Cell Biol. 2015;7:92–104.
Nishimura K, Kumazawa T, Kuroda T, Katagiri N, Tsuchiya M, Goto N, Furumai R, Murayama A, Yanagisawa J, Kimura K. Perturbation of ribosome biogenesis drives cells into senescence through 5S RNP-mediated p53 activation. Cell Rep. 2015;10:1310–23.
Moore M, Gammon L, Dreger S, Koh J, Garbe J, Stampfer M, Philpott M, Jones L, Bishop C. Ribosomal stress-induced senescence as a novel pro-senescence strategy for p16 positive basal-like breast cancer. bioRxiv. 2018. https://doi.org/10.1101/469445.
Huang M, Wang J, Torre E, Dueck H, Shaffer S, Bonasio R, Murray JI, Raj A, Li M, Zhang NR. SAVER: gene expression recovery for single-cell RNA sequencing. Nat Methods. 2018;15:539–42.
Satija R, Farrell JA, Gennert D, Schier AF, Regev A. Spatial reconstruction of single-cell gene expression data. Nat Biotechnol. 2015;33:495–502.
Macosko EZ, Basu A, Satija R, Nemesh J, Shekhar K, Goldman M, Tirosh I, Bialas AR, Kamitaki N, Martersteck EM, Trombetta JJ, Weitz DA, Sanes JR, Shalek AK, Regev A, McCarroll SA. Highly parallel genome-wide expression profiling of individual cells using nanoliter droplets. Cell. 2015;161:1202–14.
Korsunsky I, Millard N, Fan J, Slowikowski K, Zhang F, Wei K, Baglaenko Y, Brenner M, Loh P-R, Raychaudhuri S. Fast, sensitive and accurate integration of single-cell data with Harmony. Nat Methods. 2019;16:1289–96.
La Manno G, Soldatov R, Zeisel A, Braun E, Hochgerner H, Petukhov V, Lidschreiber K, Kastriti ME, Lonnerberg P, Furlan A, Fan J, Borm LE, Liu Z, van Bruggen D, Guo J, He X, Barker R, Sundstrom E, Castelo-Branco G, Cramer P, Adameyko I, Linnarsson S, Kharchenko PV. RNA velocity of single cells. Nature. 2018;560:494–8.
Blockhuys S, Celauro E, Hildesjö C, Feizi A, Stål O, Fierro-González JC, Wittung-Stafshede P. Defining the human copper proteome and analysis of its expression variation in cancers. Metallomics. 2017;9:112–23.
Hernandez-Segura A, de Jong TV, Melov S, Guryev V, Campisi J, Demaria M. Unmasking transcriptional heterogeneity in senescent cells. Curr Biol. 2017;27:2652-2660.e2654.
Travaglini KJ, Nabhan AN, Penland L, Sinha R, Gillich A, Sit RV, Chang S, Conley SD, Mori Y, Seita J, Berry GJ, Shrager JB, Metzger RJ, Kuo CS, Neff N, Weissman IL, Quake SR, Krasnow MA. A molecular cell atlas of the human lung from single-cell RNA sequencing. Nature. 2020;587:619–25.
KadurLakshminarasimha Murthy P, Sontake V, Tata A, Kobayashi Y, Macadlo L, Okuda K, Conchola AS, Nakano S, Gregory S, Miller LA, Spence JR, Engelhardt JF, Boucher RC, Rock JR, Randell SH, Tata PR. Human distal lung maps and lineage hierarchies reveal a bipotent progenitor. Nature. 2022;604:111–9.
Tsukui T, Sun KH, Wetter JB, Wilson-Kanamori JR, Hazelwood LA, Henderson NC, Adams TS, Schupp JC, Poli SD, Rosas IO, Kaminski N, Matthay MA, Wolters PJ, Sheppard D. Collagen-producing lung cell atlas identifies multiple subsets with distinct localization and relevance to fibrosis. Nat Commun. 2020;11:1920.
Saelens W, Cannoodt R, Todorov H, Saeys Y. A comparison of single-cell trajectory inference methods. Nat Biotechnol. 2019;37:547–54.
Haghverdi L, Buttner M, Wolf FA, Buettner F, Theis FJ. Diffusion pseudotime robustly reconstructs lineage branching. Nat Methods. 2016;13:845–8.
Chambers KJ, Tonkin LA, Chang E, Shelton DN, Linskens MH, Funk WD. Identification and cloning of a sequence homologue of dopamine beta-hydroxylase. Gene. 1998;218:111–20.
Thannickal VJ, Horowitz JC. Evolving concepts of apoptosis in idiopathic pulmonary fibrosis. Proc Am Thorac Soc. 2006;3:350–6.
Álvarez D, Cárdenes N, Sellarés J, Bueno M, Corey C, Hanumanthu VS, Peng Y, D’Cunha H, Sembrat J, Nouraie M, Shanker S, Caufield C, Shiva S, Armanios M, Mora AL, Rojas M. IPF lung fibroblasts have a senescent phenotype. Am J Physiol Lung Cell Mol Physiol. 2017;313:L1164–73.
Xie X, Percipalle P. Elevated transforming growth factor beta signaling activation in beta-actin-knockout mouse embryonic fibroblasts enhances myofibroblast features. J Cell Physiol. 2018;233:8884–95.
Lai X, Zhao H, Zhang Y, Guo K, Xu Y, Chen S, Zhang J. Intranasal delivery of copper oxide nanoparticles induces pulmonary toxicity and fibrosis in C57BL/6 mice. Sci Rep. 2018;8:4499.
Lahiri DK, Maloney B, Greig NH. Are pulmonary fibrosis and Alzheimer’s disease linked? Shared dysregulation of two miRNA species and downstream pathways accompany both disorders. J Biol Chem. 2017;292:20353–20353.
Wong SL, Sukkar MB. The SPARC protein: an overview of its role in lung cancer and pulmonary fibrosis and its potential role in chronic airways disease. Br J Pharmacol. 2017;174:3–14.
Aumiller V, Strobel B, Romeike M, Schuler M, Stierstorfer BE, Kreuz S. Comparative analysis of lysyl oxidase (like) family members in pulmonary fibrosis. Sci Rep. 2017;7:149.
Pardo A, Selman M. Lung fibroblasts, aging, and idiopathic pulmonary fibrosis. Ann Am Thorac Soc. 2016;13(Suppl 5):S417-s421.
Xin X, Mains RE, Eipper BA. Monooxygenase X, a member of the copper-dependent monooxygenase family localized to the endoplasmic reticulum. J Biol Chem. 2004;279:48159–67.
Kuhn C, Mason RJ. Immunolocalization of SPARC, tenascin, and thrombospondin in pulmonary fibrosis. Am J Pathol. 1995;147:1759–69.
Brewer GJ. Risks of copper and iron toxicity during aging in humans. Chem Res Toxicol. 2010;23:319–26.
Barnham KJ, Bush AI. Metals in Alzheimer’s and Parkinson’s diseases. Curr Opin Chem Biol. 2008;12:222–8.
Niu Y-Y, Zhang Y-Y, Zhu Z, Zhang X-Q, Liu X, Zhu S-Y, Song Y, Jin X, Lindholm B, Yu C. Elevated intracellular copper contributes a unique role to kidney fibrosis by lysyl oxidase mediated matrix crosslinking. Cell Death Dis. 2020;11:211.
Martinez FJ, Collard HR, Pardo A, Raghu G, Richeldi L, Selman M, Swigris JJ, Taniguchi H, Wells AU. Idiopathic pulmonary fibrosis. Nat Rev Dis Prim. 2017;3:17074.
Chothani S, Schafer S, Adami E, Viswanathan S, Widjaja AA, Langley SR, Tan J, Wang M, Quaife NM, Jian Pua C, D’Agostino G, Guna Shekeran S, George BL, Lim S, Yiqun Cao E, van Heesch S, Witte F, Felkin LE, Christodoulou EG, Dong J, Blachut S, Patone G, Barton PJR, Hubner N, Cook SA, Rackham OJL. Widespread translational control of fibrosis in the human heart by RNA-binding proteins. Circulation. 2019;140:937–51.
Larsson O, Diebold D, Fan D, Peterson M, Nho RS, Bitterman PB, Henke CA. Fibrotic myofibroblasts manifest genome-wide derangements of translational control. PLoS ONE. 2008;3: e3220.
Xue S, Barna M. Specialized ribosomes: a new frontier in gene regulation and organismal biology. Nat Rev Mol Cell Biol. 2012;13:355–69.
Shi Z, Fujii K, Kovary KM, Genuth NR, Röst HL, Teruel MN, Barna M. Heterogeneous ribosomes preferentially translate distinct subpools of mRNAs genome-wide. Mol Cell. 2017;67:71-83.e77.
Bastide A, David A. The ribosome, (slow) beating heart of cancer (stem) cell. Oncogenesis. 2018;7:34.
Kulkarni S, Dolezal JM, Wang H, Jackson L, Lu J, Frodey BP, Dosunmu-Ogunbi A, Li Y, Fromherz M, Kang A, Santana-Santos L, Benos PV, Prochownik EV. Ribosomopathy-like properties of murine and human cancers. PLoS ONE. 2017;12: e0182705.
Schafer MJ, White TA, Iijima K, Haak AJ, Ligresti G, Atkinson EJ, Oberg AL, Birch J, Salmonowicz H, Zhu Y, Mazula DL, Brooks RW, Fuhrmann-Stroissnigg H, Pirtskhalava T, Prakash YS, Tchkonia T, Robbins PD, Aubry MC, Passos JF, Kirkland JL, Tschumperlin DJ, Kita H, LeBrasseur NK. Cellular senescence mediates fibrotic pulmonary disease. Nat Commun. 2017;8:14532.
Merkt W, Bueno M, Mora AL, Lagares D. Senotherapeutics: targeting senescence in idiopathic pulmonary fibrosis. Semin Cell Dev Biol. 2020;101:104–10.
This work was funded by the National Institutes of Health grant U01HL145550 to PVB, RL, ALM, MR; and R01HL157879, R01HL159805, R01DK130294, R01HL127349 to PVB.
Ethics approval and consent to participate
All data obtained from previously published samples from our work or work of others.
Consent for publication
RL served as a consultant for Bristol Myers Squibb, Formation, Sanofi, Biocon, Boehringer-Mannheim, Merck and Genentech/Roche, and holds or recently had research grants from Corbus, Formation, Elpidera, Regeneron, Pfizer and Kiniksa.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Jia, M., Rosas, L., Kapetanaki, M.G. et al. Early events marking lung fibroblast transition to profibrotic state in idiopathic pulmonary fibrosis. Respir Res 24, 116 (2023). https://doi.org/10.1186/s12931-023-02419-0