- Open Access
Digital quantification of p16-positive foci in fibrotic interstitial lung disease is associated with a phenotype of idiopathic pulmonary fibrosis with reduced survival
Respiratory Research volume 23, Article number: 147 (2022)
Idiopathic pulmonary fibrosis (IPF) is associated with increased expression of cyclin-dependent kinase inhibitors such as p16 and p21, and subsequent induction of cell cycle arrest, cellular senescence, and pro-fibrotic gene expression. We sought to link p16-expression with a diagnosis of IPF or other fibrotic interstitial lung diseases (ILDs), radiographic pattern, senescent foci-specific gene expression, antifibrotic therapy response, and lung transplant (LTx)-free survival.
Eighty-six cases of fibrosing ILD were identified with surgical lung biopsy. Immunohistochemistry for p16 was performed on sections with the most active fibrosis. p16-positive foci (loose collection of p16-positive fibroblasts with overlying p16-positive epithelium) were identified on digital slides and quantified. Cases were scored as p16-low (≤ 2.1 foci per 100 mm2) or p16-high (> 2.1 foci per 100 mm2). Twenty-four areas including senescent foci, fibrotic and normal areas were characterized using in situ RNA expression analysis with digital spatial profiling (DSP) in selected cases.
The presence of p16-positive foci was specific for the diagnosis of IPF, where 50% of cases expressed any level of p16 and 26% were p16-high. There was no relationship between radiographic pattern and p16 expression. However, there was increased expression of cyclin-dependent kinase inhibitors, collagens and matrix remodeling genes within p16-positive foci, and cases with high p16 expression had shorter LTx-free survival. On the other hand, antifibrotic therapy was significantly protective. DSP demonstrated that fibroblastic foci exhibit transcriptional features clearly distinct from that of normal-looking and even fibrotic areas.
We demonstrated the potential clinical applicability of a standardized quantification of p16-positive fibroblastic foci. This method identifies an IPF phenotype associated with foci-specific upregulation of senescence-associated and matrix remodeling gene expression. While these patients have reduced LTx-free survival, good response to antifibrotic therapies was observed in those who were treated.
Interstitial lung disease (ILD) encompasses a spectrum of conditions characterized by inflammation and fibrosis of the pulmonary interstitium. Although there are many forms of fibrotic ILD, the most common pathologic patterns are usual interstitial pneumonia (UIP) and nonspecific interstitial pneumonia (NSIP). Many cases of fibrotic ILD are diagnosed solely based upon their imaging characteristics , but there are a subset of cases that require surgical lung biopsy (SLB) for definitive characterization of the underlying ILD. Accurate classification is essential due to the significant difference in the treatment and prognosis of different forms of ILD. However, there are often overlapping features of different types of fibrotic ILD, making these cases challenging .
Many cases of UIP can be clinically classified as idiopathic pulmonary fibrosis (IPF) in the absence of any other known cause of a UIP pattern of injury (including autoimmune diseases). IPF carries a poor prognosis in the absence of therapy, with a median survival of 3–5 years from time of diagnosis [3,4,5]. Antifibrotic agents have some limited benefit in delaying the progression of IPF [6, 7], but the only definitive therapy of IPF is lung transplant (LTx) .
Histologically, UIP is characterized by a patchy subpleural fibrosis that is temporally and spatially heterogeneous. Areas of older, dense fibrosis with microscopic honeycombing are intermixed with areas of active fibrosis characterized by fibroblastic foci and uninvolved lung parenchyma. UIP patterns predict a worse prognosis [9, 10] when compared to other patterns on SLB. However, conflicting data have been published regarding the hypothesis that other morphological findings, such as the extent of traction bronchiectasis  or density of fibroblastic foci [12, 13], can predict outcome. The role of fibroblastic foci in the progression of pulmonary fibrosis may instead depend on their activity and current pathobiological role in each individual case. Therefore, none of these findings appears consistently predictive for prognostic outcome and are not routinely assessed in clinical practice.
Abnormal cellular senescence has been identified as a keystone process in the development of pulmonary fibrosis [14,15,16], and a subset of fibrotic ILD has been linked to abnormal cellular senescence. A family of cyclin-dependent kinase (CDK) inhibitors (including p15, p16 and p21) initiate cell cycle arrest in response to context or cell type specific triggers, such as DNA damage and other cellular stressors . The induction of cellular senescence in these cases of pulmonary fibrosis may be mediated by p16 activation [18, 19], arresting the cell cycle through activation of this retinoblastoma-mediated CDK 4/6 pathway. Although this growing body of evidence supports the role of senescence in driving the development of pulmonary fibrosis, the specificity of senescence for IPF and relevance of senescence across the spectrum of fibrotic ILD are both unknown. None of these studies have yet provided a means to quantify these foci in a clinically deployable test, using standardized methodology.
Assays assessing for p16 expression are commonly used in many clinical laboratories, serving as a marker of high-grade squamous lesions associated with human Papilloma virus infection. These well-developed immunohistochemical tools can detect p16 expression, but are unavailable for many of the other CDK inhibitors. By leveraging these assays for p16, we developed a reproducible and clinically deployable method for identifying and quantifying senescent foci defined by p16 expression in epithelial cells and their underlying myofibroblasts. We therefore sought to apply this method on SLBs to identify a UIP pattern, and demonstrate a “pro-senescent” phenotype potentially associated with distinct HRCT pattern, RNA expression profile, prognosis, and response to therapy. Specifically, digital spatial RNA high-resolution profiling with the GeoMX platform was used for the first time in fibrotic ILD to explore transcriptional differences in fibroblastic foci vs. fibrotic and normal-looking areas in IPF cases.
Materials and methods
Eighty-six cases of ILD diagnosed with SLB were identified between 2003 and 2018 at a tertiary-care level hospital (Table 1, Fig. 1). The general diagnosis of fibrotic ILD was first determined by chart review. After excluding all known causes of ILD (including autoimmune disease) and in the absence of a typical UIP radiographic pattern on high resolution CT scan (HRCT), patients were referred to SLB to clarify the specific type of fibrotic ILD. The final diagnosis of each case of fibrotic ILD was based on clinical-radiographic-pathologic multi-disciplinary discussion (MDD) .
All patient characteristics were collected at time of biopsy. Pulmonary function tests (PFTs) were performed according to European Respiratory Society and American Thoracic Society guidelines [21,22,23,24]. The HRCT pattern was evaluated and classified by a chest radiologist (MZ) who was blinded from clinical information, and classified as either probable UIP, indeterminate for UIP or inconsistent for UIP . The overall extent of pulmonary fibrosis, including reticular, ground glass opacities and honeycomb changes, was then evaluated with a quantitative visual score at 4 pre-defined levels (aortic arch, bronchus intermedius, pulmonary veins, lowest scan) in each hemithorax. A 5-point scale representing extent of disease (0 = absence of lesions, 1, 2, 3, 4 = extent of lesions, respectively, 0–25%, 25–50%, 50–75%, 75–100%) was used for the score in each location. The component scores for each scan were then summed, creating a final overall score (100/maximum predicted value, equivalent to 8 times the number of scans performed). The extent of fibrosis was then expressed as a percentage of the total lung volume . Subjects with incomplete follow-up data were excluded.
Formalin-fixed paraffin-embedded tissue sampled from the SLB for each patient was used for immunohistochemical analysis. Haematoxylin and eosin (H&E) stained slides from each case were re-examined and reviewed by a pulmonary pathologist (MJC) blinded from clinical information. Sections with the greatest number of fibroblastic foci or the most active fibrosis were selected. Four micrometer-thick sections were then cut from these tissues for p16 (monoclonal antibody E6H4, Roche, Basel, Switzerland), p21 (monoclonal antibody EPR362, Abcam, Cambridge, United kingdom) and purified retinoblastoma protein (pRB)(monoclonal antibody G3-245, Biosciences, Mississauga, Canada) immunohistochemistry.
The immunohistochemical-stained slides and corresponding H&E slides were digitized using an Aperio ScanScope AT Turbo (Leica Biosystems, Wetzlar, Germany). Using these digital slides, p16-positive foci were defined as loose collections of p16-positive fibroblasts with an overlying p16-positive overlying flat (Fig. 2A), simple cuboidal (Fig. 2B and C) or columnar epithelium (Fig. 2D). Other areas of p16-positivity were disregarded. The total area of lung tissue and location of each p16-positive focus from each case was annotated on the digital slide (Additional file 2: Fig. S1) using the open access QuPath software package . The total surface area was then calculated based on the annotated areas in a semi-automated protocol. The number of p16-positive foci for each case was normalized/100 mm2 lung tissue, and cases were scored as p16-low (≤ 2.1 foci per 100 mm2 lung tissue) or p16-high (> 2.1 foci per 100 mm2 lung tissue), based on receiver operating characteristics analysis (c-statistics) against survival.
Digital spatial RNA profiling
The GeoMX Digital Spatial Profiling (DSP)(NanoString Technologies, Seattle, WA) is a method providing highly multiplex spatial profiling of RNAs on formalin-fixed, paraffin-embedded (FFPE) samples . The DSP platform quantifies the abundance of RNA by counting unique indexing oligonucleotides assigned to each target of interest. The GeoMX DSP allows for up to 800-plex profiling of mRNA or protein using an optical-barcode readout . Slide-mounted FFPE tissue sections underwent antigen retrieval and incubation with indexing oligonucleotides covalently attached to mRNA hybridization probes with a UV-photocleavable (PC) linker . The same sections were stained with fluorescently labeled imaging reagents to identify tissue features and markers of interest. The system is highly sensitive and non-destructive, allowing multiple cycles of high-plex profiling on the same tissue section or subsequent DNA sequencing on the same section.
We selected 24 lung microenvironments from FFPE sections in 2 cases of IPF, both with a high density of p16-positive foci. Tissue sections were stained with standard morphology markers (pan-cytokeratin, alpha-smooth muscle actin, CD45) combined with oligo-tagged ISH probes (RNA) that contained a photocleavable linker. Regions of interest (ROIs) were selected, including 12 fibroblastic foci, 6 areas of pulmonary fibrosis and 6 areas of normal lung tissue in total (Fig. 3). Areas of illumination were consequently identified by morphology markers, then sequentially exposed to UV light to decouple the oligonucleotides from the profiling reagents. Decoupled oligonucleotides were rapidly aspirated using a microcapillary without touching the sample, thereby leaving the sample unaltered. Oligonucleotides were deposited into wells of a microtiter plate, and the information contained within each well were indexed to the ROI on the tissue. The oligonucleotides were hybridized to NanoString® barcodes and sequenced using standard next generation sequencing technologies.
For quality control purposes, counts from each ROI were analyzed mRNAs were targeted by multiple probes. Outlier probes and under-sequenced samples were dropped from downstream data analysis. RNA was detected in a reproducible manner with high correlation to fluorescent in situ hybridization. Semi-quantitative analysis  of these ROIs was then performed by standard next generation oligonucleotide sequencing. In terms of normalization for reproducibility purposes, individual counts were normalized against the signal detected from 31 housekeeping genes across the study. Individual counts were also be normalized against the 75th percentile of signal from their own ROI.
The Cancer Transcriptome Atlas Assay (NanoString®) was used, as considered the most inclusive of genes of interest for pulmonary fibrosis. This atlas included 1825 genes, and principal component analysis mapping was carried out with the GeoMX software (NanoString®). Differentially expressed genes were selected based on both a q-value < 0.05 and a fold change of ≥ 1.5. For hierarchic clustering, Cluster 3.0 and Treeview (Eisen’s Laboratory, Stanford University, Stanford, CA) were used (http://rana.lbl.gov/eisen/?page_id=7). Unsupervised clustering does not take any of the experimental variables such as treatment, phenotype, tissue, etc. into account. For Gene-set Enrichment Analysis (GSEA), C5 GeneOntology was used as gene-set database .
Pearson’s correlation was used to compare p16 expression and HRCT score. The Kolmogorov–Smirnov test assessed variables’ distribution. The Student's t-test, the χ2 method, and Fisher's exact test were used, when appropriate, to compare continuous and categorical variables. Univariate and multivariate Cox proportional hazards regression analyses were used to identify the significant predictors of LTx-free survival. Receiver operating characteristic (ROC) analysis (c-statistics) was used to determine best cut-off of p16 expression to predict LTx-free survival, by examining accuracy (sum of sensitivity and specificity). LTx-free survival from the time of biopsy was also evaluated using Kaplan–Meier curves and the log-rank test. p-values < 0.05 were regarded as significant. For RNA sequencing data, the q-value was adopted as a multiple comparison correction and used to identify differentially expressed genes. The q-value represents the minimal false discovery ratio (FDR) at which an individual hypothesis test may be called significant. Prism-8 (GraphPad, La Jolla, CA), JMP (SAS Institute, Cary, NC) and MedCalc (Mariakerke, Belgium) softwares were used.
Patients characteristics and mortality
After SLB and MDD, 52 cases were diagnosed with IPF and 34 with fibrotic ILDs other than IPF (Table 1, Fig. 1). The average age of patients was 62 years, with a male predominance (56%). On average, these patients had a mild restrictive ventilatory defect (FVC 75%, TLC 70%), with moderate reduction of diffusing lung capacity (48%). HRCT visual score analysis demonstrated an average pre-biopsy extent of fibrosis of 20%.
Post-biopsy follow-up was 34 ± 35 months. Thirty-one patients (36%) were treated with antifibrotic therapy (either pirfenidone or nintedanib), 25 (29%) with immunomodulatory therapy (either prednisone, mycophenolate mofetil, azathioprine or cyclophosphamide) and 30 (35%) did not receive specific treatment (Fig. 1). Of note, several patients were diagnosed with IPF before the approval of pirfenidone (2013) and nintedanib (2015) in Canada. During the post-biopsy observation period, mortality was 49% and 13 patients (15%) underwent LTx (Fig. 1). The combined incidence of death and LTx was therefore 64%.
p16-positive foci specificity to IPF
Positive p16 immunohistochemistry was consistently observed on the epithelium overlying numerous, but not all, fibroblastic foci in IPF cases (Fig. 2).
Any number of p16-positive foci were identified in 40 of the 52 cases of IPF. The presence of any p16-positive foci was sensitive (76.9%) and specific (91.2%) for the final diagnosis of IPF against other types of ILD. There was however a variable expression of p16 in the IPF group, with a range of p16-positive foci between 0 and 33 per slide (median: 3 foci per slide), and the density of senescent foci ranged from 0.43 to 26.32 foci per 100 mm2 lung tissue, with a median density of 1.97 foci per 100 mm2 lung tissue. We selected a cut-off of 2.1 foci per 100 mm2 to maximize the differences in outcome between the two groups, in accordance with the ROC analysis (Additional file 3: Fig. S2). As a result, 19 cases of IPF (47.5%) were classified as p16-low (≤ 2.1 foci per 100 mm2 lung tissue), and 21 cases of IPF (52.5%) were scored as p16-high (> 2.1 foci per 100 mm2 lung tissue) (Table 2).
p16-positive foci were present in only 3 cases of non-IPF fibrotic ILD. Two of the cases had unclassifiable patterns of fibrotic ILD, with 3.82 and 1.08 foci per 100 mm2 lung tissue, respectively, and the remaining one was a case of chronic hypersensitivity pneumonitis with 1 p16-positive focus (1.32 foci per 100 mm2 lung tissue). Senescent foci were absent from all other cases of non-IPF, including 16 cases of NSIP, 5 cases of chronic hypersensitivity pneumonitis, and 11 unclassifiable cases. The specificity of p16-high for the diagnosis of IPF was 97.6%. No relationship was observed between the pre-biopsy HRCT pattern and the density of p16-positive foci (Additional file 4: Fig. S3). No significant associations were identified between PFTs or HRCT fibrosis score and the density of p16-positive foci (data not shown).
Additional markers of senescence were tested with immunohistochemistry. Material was available for 6 cases for p21 and pRB staining. These cases had a range of normalized p16 foci from 0 to 26.32. While the p16 positive foci were found to express pRB (Fig. 4) the expression was not substantially above the background staining seen in the lung and is in keeping with wild-type intact pRB expression. No significant accentuation in of pRB staining was present in the foci. As shown in Fig. 4, p21 is also expressed in p16 positive fibroblastic foci. It was not possible to directly compare all the p16 positive foci to the p21 staining due to variation in the tissue composition on deeper levels. Overall 85% of fibroblastic foci identified on histologic sections have some degree of p21 staining. However, the staining is focal, patchy and challenging to appreciate above the background patchy p21 staining.
Digital spatial RNA profiling
A compact and homogenous unsupervised clustering of fibroblastic foci gene expression was observed, separating them from both the fibrotic and normal areas (Fig. 5). Consistently, principal component analysis (PCA) mapping demonstrated a definite separation across the selected gene array among fibroblastic foci, fibrotic and normal areas (Additional file 5: Fig. S4). Applying strict criteria (q-value < 0.05, fold change ≥ 1.5), 97 differentially expressed genes were identified between the fibroblastic foci and zones of dense fibrosis or normal lung. RNA expression analysis revealed a distinct upregulation in CDK inhibitor p21, extracellular matrix remodelling proteases (MMP-1, MMP-11 and ADAM12), and collagens (Collagen 1A1, Collagen 1A2 and Collagen 3A1) in fibroblastic foci, compared to both areas of dense fibrosis and normal areas (Fig. 6).
GSEA demonstrated that 107 gene sets were upregulated in fibroblastic foci vs. normal-looking areas with an FDR < 0.05. Top 5 upregulated gene sets are shown in Table 3. Not surprisingly, “scarring” was the top upregulated gene set in fibroblastic foci. In contrast, only 4 gene sets were upregulated in fibrotic areas vs. normal-looking areas with an FDR < 0.05. Top 5 upregulated gene sets are shown in Table 4. The top upregulated gene set in fibrotic areas was related to the process of endocytosis. The comparison of fibroblastic foci with fibrotic areas yielded as many as 21 upregulated gene sets with an FDR < 0.05 in the first group (Table 5).
Post-biopsy survival analysis
LTx-free survival was significantly reduced in cases with a high density of p16-positive foci when compared to cases with a low density of senescent foci (Fig. 7A). This remained true considering IPF cases only (Fig. 7B). On regression analysis, a high level of p16 expression was indeed significantly associated with LTx or mortality (HR 2.47, C.I. 1.30–4.62, p = 0.0067). ROC analysis demonstrated that high expression of p16 (> 2.1 foci per 100 mm2 lung tissue) had a 33% sensitivity (95% confidence interval, 21–47%) and 87% specificity (95% confidence interval, 70–96%) in predicting post-biopsy LTx-free survival, resulting in an accuracy (area under the curve) of 60% (p = 0.025) (Additional file 3: Fig. S3).
Considering only IPF patients (and excluding 8 subjects who developed lung cancer post-biopsy), p16 remained a risk factor towards death or LTx, but antifibrotic therapy was significantly protective (HR 2.64, p = 0.0150) on multivariate analysis (Table 6). Importantly, all patients within high p16 expression who were on anti-fibrotic therapy survived (except 1 who died of cancer), while all others either died or required a LTx.
The role of cellular senescence has become increasingly recognized as a major factor in the development of many diseases, with clinical applications for immunohistochemical markers of senescence. In this study, we were able to link immunohistochemistry for p16 to a high likehood of IPF diagnosis and also to a significantly worse outcome, when this expression was high.
Discrete foci of senescence corresponding to senescent epithelial cells and myofibroblasts were identified in recent single cell sequencing studies [14, 15]. The distribution and quantification of p16-positive foci of senescence presented as a challenge because of its low frequency and relative heterogeneity between cases. We felt that simple quantification per high powered field did not adequately standardize nor represent whether senescence appeared to be driving the underlying pathogenesis of disease. We had considered and tried to measure the number of foci in small areas (2 mm2 of lung tissue), but we had found this metric was too variable for these low frequency features. Therefore, we chose to normalize our analysis to the entire lung parenchyma, measuring the entire surface area of fibrosis on the slide using digital image analysis (QuPath). This assessment remains the most accurate and robust measure of the number of p16-positive foci while normalizing for differences in the amount of sampled lung tissue, opening its potential application to smaller cryobiopsy specimens for disease classification and prognostication . We also confirmed the expression of other markers of senescence (p21 and pRB) on p16-positive fibroblastic foci.
One interesting finding of our study is that the presence of p16 overexpression highlights only a subset of fibroblastic foci. Some groups have suggested that the number of fibroblastic foci correlates with poor outcome . However, the definition of a fibroblastic focus remains heterogeneous, with interobserver variability . Therefore, this criterion alone may have limitations as a measure of disease burden. We found that using immunohistochemistry, p16 overexpression highlights a subset of fibroblastic foci, allowing for easier identification of these foci. This technique may be leveraged to clearly define these fibroblastic foci and may make this task more amenable to automation. Importantly, not all fibroblastic foci are positive for p16, suggesting that these collections of p16-positive fibroblasts are driving this disease process through cellular senescence.
We then sought to characterize the features of senescence-associated IPF. A high density of p16-positive foci was associated with a significant worse outcome, despite having no association with any specific HRCT pattern. Our data also demonstrated improved outcomes in patients with a high density of p16-positive foci that were started on these antifibrotics, serving as a marker of increased response from antifibrotic therapy. Validation of these findings in larger cohorts is necessary to determine the value of p16 in guiding management and may support a rationale to obtain the “p16 status” in biopsies for IPF.
Gene expression studies on bulk lung tissues from IPF patients have previously identified increases in extracellular matrix (ECM) remodelling genes, fibroblast proliferation and collagen deposition compared to other types of ILD and normal lung tissue [14, 32,33,34]. Our data builds on this work, assessing 24 different lung microenvironments from 2 IPF lung samples, using the new DSP methodology. A transcriptional increase of CDK inhibitor p21, ECM proteins and matrix metalloproteinases within the fibroblastic foci of senescence-associated IPF was observed. PCA mapping and hierarchical clustering showed a clear demarcation of fibroblastic foci from both normal areas and areas of fibrosis. These data, although only preliminary, highlight the importance of in situ tissue transcriptional analysis and support the hypothesis that fibroblastic foci play an independent role in promoting senescence-associated ECM remodelling in a UIP pattern of lung fibrosis. As a limitation, our data was based on the Cancer Transcriptome Atlas Assay, the largest available panel for FFPE tissues at the time of the study, rather than whole transcriptome. Future work to optimize FFPE and fresh-frozen tissues for transcriptome analysis may highlight pathways that were not identified with these techniques, such as the role of other ECM remodelling proteins or inhibitors that were previously classified in these original studies.
In IPF, cells reported to express senescence markers include fibroblasts and alveolar epithelial type 2 cells . Single-cell RNA studies in IPF identified aberrant, cuboidal epithelial expressing multiple markers of senescence and localized at the edge of fibroblastic foci, exactly where the p16 immunostaining was observed . Importantly, these cells exhibited features of epithelial-to-mesenchymal transition , linking senescence to a fundamental mechanism of the fibrogenic process. Animal model data indeed suggested that senescence and depletion of alveolar epithelial type 2 cells can lead to spontaneous lung fibrosis . Myofibroblasts within fibroblastic foci in IPF display an apoptotic-resistant phenotype , a feature compatible with senescence . However, a diminished proliferative capacity was observed in IPF lung myofibroblasts and is inconsistent with the proliferative process suggested in the fibroblast foci . Overall, it appears that, in IPF, senescence primarily takes place in epithelial cells, but, through SASP, could have a long range of fibrogenic consequences in the microenvironment. Further studies with spatially-resolved transcriptomics on lung biopsies in treatment-naïve patients will likely help clarifying this issue and identify new targets for treatment. From a therapeutic perspective, given consistent transcriptional and immunohistochemical findings, there is an increasing interest in the use of senolytic drugs in IPF. A preliminary experience from Justice et al. with dasatinib plus quercetin demonstrated the feasibility of such trials .
When considering the striking difference in gene sets expression between fibroblastic foci and fibrotic areas, we should consider that the ROI “fibroblastic foci” do include the overlying epithelium. In this regard, the upregulation of the “cadherin binding” gene set in fibroblastic foci (vs. fibrotic areas) is very interesting, as it points again to their central fibrogenic role of epithelial-to-mesenchymal transition, where cadherin modulation is heavily involved . Another important consideration in the interpretation of GSEA results is that “normal-looking areas” in IPF probably do not represent normal controls. A study from Luzina et al. suggested that normal-appearing areas are actually fully engaged in the disease process, with transcriptional expression of connective tissue, immune and inflammatory-related genes .
DSP analysis demonstrated a greater upregulation of p21 than p16 in fibroblastic foci. This is unsurprising, since recent studies have highlighted increased levels of both p16  and p21  in response to critically shortened telomeres, and various drivers of senescence via the p53-dependent cellular senescence pathway . Other groups have been able to highlight p21-positive cells with immunofluorescence in cases of IPF . While p21 appears to have a more significant difference at the gene expression level, it does not appear to translate well at the protein level, given the challenges associated with p21 outlined above. p21 is well known to be regulated at the protein level by the activity of E3 ubiquitin ligases, such as Skp1, which could explain the observed differences . The expression of pRB also appears to be intact on fibroblastic foci and would provide a means for cell cycle blockade in these non-neoplastic cells. We also previously investigated p53 immunostaining, but patchy, wild-type staining was observed on the slides and p53 failed to differentiate UIP from NSIP .
We demonstrated that digital quantification of p16-positive foci is a novel, simple, and inexpensive marker of IPF associated with cellular senescence. We were able to highlight some of the molecular pathways behind senescence-associated IPF in fibroblastic foci, identifying a family of matrix remodelling proteins and CDK inhibitors that may be involved in the rapid progression of fibrosis, leading to worse LTx-free survival. These data suggest that fibroblastic foci-specific activation of cellular senescence is a critical driver for the ongoing fibrosis in IPF. However, encouragingly, the regression analysis in IPF patients demonstrated that this driver is potentially treatable, given the protective effect of anti-fibrotic therapy, when started early post-biopsy.
In conclusion, the identification of p16-positive fibroblastic foci on lung biopsies may be diagnostically helpful and prognostically meaningful, identifying a subgroup of patients who may benefit the most from existing antifibrotic therapies, and potentially from future senolytic therapies.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its additional information files. The original datasets containing patient data that were used and analysed during the current study are available from the corresponding author on reasonable request.
Digital spatial profiler
False discovery ratio
Gene set enrichment analysis
High resolution chest CT scan
Interstitial lung disease
Idiopathic pulmonary fibrosis
Nonspecific interstitial pneumonia
Purified retinoblastoma protein
Pulmonary function test
Receiver operating characteristic
Regions of interest
Senescence-associated signalling protein
Surgical lung biopsy
Usual interstitial pneumonia
Lynch DA, Sverzellati N, Travis WD, et al. Diagnostic criteria for idiopathic pulmonary fibrosis: a Fleischner Society White Paper. Lancet Respir Med. 2018;6:138–53. https://doi.org/10.1016/S2213-2600(17)30433-2.
Adegunsoye A, Vij R, Noth I. Integrating genomics into management of fibrotic interstitial lung disease. Chest. 2019;155:1026–40. https://doi.org/10.1016/j.chest.2018.12.011.
Frankel SK, Schwarz MI. Update in idiopathic pulmonary fibrosis. Curr Opin Pulm Med. 2009;15:463–9. https://doi.org/10.1097/MCP.0b013e32832ea4b3.
Ley B, Collard HR, King TE. Clinical course and prediction of survival in idiopathic pulmonary fibrosis. Am J Respir Crit Care Med. 2011;183:431–40. https://doi.org/10.1164/rccm.201006-0894CI.
Mura M, Porretta MA, Bargagli E, et al. Predicting survival in newly diagnosed idiopathic pulmonary fibrosis: a 3-year prospective study. Eur Respir J. 2012;40:101–9. https://doi.org/10.1183/09031936.00106011.
King TE, Bradford WZ, Castro-Bernardini S, et al. A phase 3 trial of pirfenidone in patients with idiopathic pulmonary fibrosis. N Engl J Med. 2014;370:2083–92. https://doi.org/10.1056/nejmoa1402582.
Richeldi L, du Bois RM, Raghu G, et al. Efficacy and safety of nintedanib in idiopathic pulmonary fibrosis. N Engl J Med. 2014;370:2071–82. https://doi.org/10.1056/nejmoa1402584.
Thabut G. Survival after bilateral versus single-lung transplantation for idiopathic pulmonary fibrosis. Ann Intern Med. 2009;151:767–74. https://doi.org/10.7326/0003-4819-151-11-200912010-00004.
Tiitto L, Heiskanen U, Bloigu R, et al. Thoracoscopic lung biopsy is a safe procedure in diagnosing usual interstitial pneumonia. Chest. 2005;128:2375–80. https://doi.org/10.1378/chest.128.4.2375.
Monaghan H, Wells AU, Colby TV, et al. Prognostic implications of histologic patterns in multiple surgical lung biopsies from patients with idiopathic interstitial pneumonias. Chest. 2004;125:522–6. https://doi.org/10.1378/chest.125.2.522.
Jacob J, Aksman L, Mogulkoc N, et al. Serial CT analysis in idiopathic pulmonary fibrosis: comparison of visual features that determine patient outcome. Thorax. 2020;75:648–54. https://doi.org/10.1136/thoraxjnl-2019-213865.
Flaherty KR, Colby TV, Travis WD, et al. Fibroblastic foci in usual interstitial pneumonia: Idiopathic versus collagen vascular disease. Am J Respir Crit Care Med. 2003;167:1410–5. https://doi.org/10.1164/rccm.200204-373OC.
Hanak V, Ryu JH, de Carvalho E, et al. Profusion of fibroblast foci in patients with idiopathic pulmonary fibrosis does not predict outcome. Respir Med. 2008;102:852–6. https://doi.org/10.1016/j.rmed.2008.01.012.
Adams TS, Schupp JC, Poli S, et al. 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.
Yao C, Guan X, Carraro G, et al. Senescence of alveolar type 2 cells drives progressive pulmonary fibrosis. Am J Respir Crit Care Med. 2021;203:707–17. https://doi.org/10.1164/rccm.202004-1274OC.
Lee JS, La J, Aziz S, et al. Molecular markers of telomere dysfunction and senescence are common findings in the usual interstitial pneumonia pattern of lung fibrosis. Histopathology. 2021;79:67–76. https://doi.org/10.1111/his.14334.
Barnes PJ, Baker J, Donnelly LE. Cellular senescence as a mechanism and target in chronic lung diseases. Am J Respir Crit Care Med. 2019;200:556–64. https://doi.org/10.1164/rccm.201810-1975TR.
Cecchini MJ, Hosein K, Howlett CJ, et al. Comprehensive gene expression profiling identifies distinct and overlapping transcriptional profiles in non-specific interstitial pneumonia and idiopathic pulmonary fibrosis. Respir Res. 2018;19:153. https://doi.org/10.1186/s12931-018-0857-1.
Schafer MJ, White TA, Iijima K, et al. Cellular senescence mediates fibrotic pulmonary disease. Nat Commun. 2017. https://doi.org/10.1038/ncomms14532.
Flaherty KR, King TE, Raghu G, et al. Idiopathic interstitial pneumonia: what is the effect of a multidisciplinary approach to diagnosis? Am J Respir Crit Care Med. 2004;170:904–10. https://doi.org/10.1164/rccm.200402-147OC.
Crapo O, Hankinson JL, Irvin C, et al. Standardization of spirometry: 1994 update. Am J Respir Crit Care Med. 1995;152:1107–36. https://doi.org/10.1164/ajrccm.152.3.7663792.
Miller MR, Hankinson J, Brusasco V, et al. Standardisation of spirometry. Eur Respir J. 2005;26:319–38. https://doi.org/10.1183/09031936.05.00034805.
Wanger J, Clausen JL, Coates A, et al. Standardisation of the measurement of lung volumes. Eur Respir J. 2005;26:511–22. https://doi.org/10.1183/09031936.05.00035005.
Zavorsky GS, Hsia CCW, Hughes JMB, et al. Standardisation and application of the single-breath determination of nitric oxide uptake in the lung. Eur Respir J. 2017. https://doi.org/10.1183/13993003.00962-2016.
Raghu G, Remy-Jardin M, Myers JL, et al. Diagnosis of idiopathic pulmonary fibrosis an official ATS/ERS/JRS/ALAT clinical practice guideline. Am J Respir Crit Care Med. 2018;198:e44-68. https://doi.org/10.1164/rccm.201807-1255ST.
Taha N, D’Amato D, Hosein K, et al. Longitudinal functional changes with clinically significant radiographic progression in idiopathic pulmonary fibrosis: are we following the right parameters? Respir Res. 2020. https://doi.org/10.1186/s12931-020-01371-7.
Bankhead P, Loughrey MB, Fernández JA, et al. QuPath: open source software for digital pathology image analysis. Sci Rep. 2017. https://doi.org/10.1038/s41598-017-17204-5.
Merritt CR, Ong GT, Church SE, et al. Multiplex digital spatial profiling of proteins and RNA in fixed tissue. Nat Biotechnol. 2020;38:586–99. https://doi.org/10.1038/s41587-020-0472-9.
Beechem JM. High-Plex Spatially resolved RNA and protein detection using digital spatial profiling: a technology designed for immuno-oncology biomarker discovery and translational research. In: Methods in Molecular Biology. Humana Press Inc. 2020;563–583. doi:https://doi.org/10.1007/978-1-4939-9773-2_25
Ashburner M, Ball CA, Blake JA, et al. Consortium GO. Gene Ontology: Tool for the unification of biology. Nat Genet. 2000;25:25–29.
Troy LK, Grainge C, Corte TJ, et al. Diagnostic accuracy of transbronchial lung cryobiopsy for interstitial lung disease diagnosis (COLDICE): a prospective, comparative study. Lancet Respir Med. 2020;8:171–81. https://doi.org/10.1016/S2213-2600(19)30342-X.
Selman M, Pardo A, Barrera L, et al. Gene expression profiles distinguish idiopathic pulmonary fibrosis from hypersensitivity pneumonitis. Am J Respir Crit Care Med. 2006;173:188–98. https://doi.org/10.1164/rccm.200504-644OC.
Pardo A, Gibson K, Cisneros J, et al. Up-regulation and profibrotic role of osteopontin in human idiopathic pulmonary fibrosis. PLoS Med. 2005;2:0891–903. https://doi.org/10.1371/journal.pmed.0020251.
Kelly MM, Leigh R, Gilpin SE, et al. Cell-specific gene expression in patients with usual interstitial pneumonia. Am J Respir Crit Care Med. 2006;174:557–65. https://doi.org/10.1164/rccm.200510-1648OC.
Alvarez D, Cardenes N, Sellares J, et al. IPF lung fibroblasts have a senescent phenotype. AmJ Physiol Lung Cell Mol Physiol. 2017;313:L1164-1173.
Justice JJ, Nambiar AM, Tchkonia T, et al. Senolytics in idiopathic pulmonary fibrosis: results from a first-in-human, open-label, pilot study. EBioMedicine. 2019;40:554–63.
Gheldolf A, Berx G. Cadherins and epithelial-to-mesenchymal transition. Prog Mol Biol Transl Sci. 2013;116:317–36.
Luzina IG, Salcedo MV, Rojas-Pena ML, et al. Transcriptomic evidence of immune activation in macroscopically normal-appearing and scarred lung tissues in idiopathic pulmonary fibrosis. Cell Immunol. 2018;325:1–13.
Cecchini MJ, Tarmey T, Ferreira A, et al. Pathology, radiology, and genetics of interstitial lung disease in patients with shortened telomeres. Am J Surg Pathol. 2021;45:871–84. https://doi.org/10.1097/PAS.0000000000001725.
Lu Z, Hunter T. Ubiquitylation and proteasomal degradation of the p21(Cip1), p27(Kip1) and p57(Kip2) CDK inhibitors. Cell Cycle. 2010;9:2342–52.
This research was funded by the Lawson Health Research Institute Internal Research Fund (INN19-006) and by the AMOSO (Academic Medical Organization of Southwestern Ontario) Innovation Fund (recipient Dr. Marco Mura for both).
Ethics approval and consent to participate
This research study was approved by the Western University Research Ethics Board (protocol n.113528) and Western University Clinical Research Impact Committee.
Consent for publication
The authors of this abstract have indicated that they have no conflicts of interest that relate to the content of this publication.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file figure legends.
Sample quantification of p16-positive foci (red) within the total lung parenchyma (yellow) using QuPath.
Recipient operating characteristic analysis of p16-positive foci against LTx-free survival. Area under the curve=0.599, p=0.025.
Stratification of p16-positive foci according to the HRCT pattern (probable, indeterminate or inconsistent for UIP). No association between density of p16-positive senescent foci and HRCT pattern could be identified.
Principal component analysis (PCA) across the analyzed portion of genome. Light brown: fibrotic areas; dark pink: fibroblastic foci; green: normal areas. PCA analysis shows a distinct separation between fibroblastic foci and fibrotic areas.
About this article
Cite this article
Keow, J., Cecchini, M.J., Jayawardena, N. et al. Digital quantification of p16-positive foci in fibrotic interstitial lung disease is associated with a phenotype of idiopathic pulmonary fibrosis with reduced survival. Respir Res 23, 147 (2022). https://doi.org/10.1186/s12931-022-02067-w
- Interstitial lung disease
- Idiopathic pulmonary fibrosis
- Usual interstitial pneumonia