Article Text

Neutrophil levels correlate with quantitative extent and progression of fibrosis in IPF: results of a single-centre cohort study
  1. Andrew Achaiah1,
  2. Emily Fraser2,
  3. Peter Saunders2,
  4. Rachel K Hoyles2,
  5. Rachel Benamore3 and
  6. Ling-Pei Ho1,2
  1. 1Translational Immunology Discovery Unit, Weatherall Institute of Molecular Medicine, University of Oxford, Oxford, UK
  2. 2Oxford Interstitial Lung Disease Service, Oxford University Hospitals NHS Foundation Trust, Oxford, UK
  3. 3Thoracic Radiology Department, Oxford University Hospitals NHS Foundation Trust, Oxford, UK
  1. Correspondence to Professor Ling-Pei Ho; ling-pei.ho{at}imm.ox.ac.uk

Abstract

Background Idiopathic pulmonary fibrosis (IPF) is a progressive fibrotic lung disease with poor prognosis. Clinical studies have demonstrated association between different blood leucocytes and mortality and forced vital capacity (FVC) decline. Here, we question which blood leucocyte levels are specifically associated with progression of fibrosis, measured by accumulation of fibrosis on CT scan using a standardised automated method.

Methods Using the Computer-Aided Lung Informatics for Pathology Evaluation and Rating CT algorithm, we determined the correlation between different blood leucocytes (<4 months from CT) and total lung fibrosis (TLF) scores, pulmonary vessel volume (PVV), FVC% and transfer factor of lung for carbon monoxide% at baseline (n=171) and with progression of fibrosis (n=71), the latter using multivariate Cox regression.

Results Neutrophils (but not monocyte or lymphocytes) correlated with extent of lung fibrosis (TLF/litre) (r=0.208, p=0.007), PVV (r=0.259, p=0.001), FVC% (r=−0.127, p=0.029) at baseline. For the 71 cases with repeat CT; median interval between CTs was 25.9 (16.8–39.9) months. Neutrophil but not monocyte levels are associated with increase in TLF/litre (HR 2.66, 95% CI 1.35 to 5.25, p=0.005).

Conclusion Our study shows that neutrophil rather than monocyte levels correlated with quantifiable increase in fibrosis on imaging of the lungs in IPF, suggesting its relative greater contribution to progression of fibrosis in IPF.

  • Innate Immunity
  • Interstitial Fibrosis

Data availability statement

Data are available on reasonable request. N/A.

https://creativecommons.org/licenses/by/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution 4.0 Unported (CC BY 4.0) license, which permits others to copy, redistribute, remix, transform and build upon this work for any purpose, provided the original work is properly cited, a link to the licence is given, and indication of whether changes were made. See: https://creativecommons.org/licenses/by/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

WHAT IS ALREADY KNOWN ON THIS TOPIC

  • Idiopathic pulmonary fibrosis (IPF) is a progressive fibrotic condition. Recently, several human studies have implicated blood leucocyte levels (monocyte, neutrophil and lymphocyte) with forced vital capacity decline and mortality. However, direct association between leucocytes and progression of fibrosis using quantitative CT analysis has not been explored.

WHAT THIS STUDY ADDS

  • This study explored the association between blood monocytes, neutrophils and lymphocytes against increase in fibrosis over time, measured using a quantitative CT algorithm, Computer-Aided Lung Informatics for Pathology Evaluation and Rating. We show that levels of blood neutrophil and lymphocytes but not monocytes were associated with greater risk of progression of fibrosis.

HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY

  • Our study shows that neutrophil rather than monocyte levels correlated with quantifiable increase in fibrosis on imaging of the lungs in IPF, suggesting its relative greater contribution to progression of fibrosis in IPF.

Introduction

Idiopathic pulmonary fibrosis (IPF) is a progressive fibrotic condition. Animal studies have implicated innate immune cells, particularly monocytes and macrophages in the pathogenesis of IPF.1 This is supported by several clinical studies which showed that higher blood monocytes levels are associated with mortality in IPF.2 Blood monocyte levels also correlated with the extent of fibrosis on CT scan,3 and Kreuter et al also found an association with a composite measure of IPF outcome (forced vital capacity (FVC) decline, 6 min walk distance reduction, acute exacerbation and/or mortality).4 Moreover, we (and others) have shown that the neutrophil:lymphocyte ratio (NLR) is also a predictor of mortality and FVC decline,5 6 and patients with higher levels of neutrophils were more likely to progress from indeterminate for Usual Interstitial Pneumonia (UIP) CT pattern to UIP pattern and a clinical diagnosis of IPF.7

However, mortality and FVC decline may not accurately reflect progression of fibrosis. Death can be due to other causes, for example, cardiovascular diseases, especially in the elderly population of IPF. Although FVC decline of >10% in IPF is the validated measure of disease progression predictive of mortality,8 and serial FVC change the recommended monitoring variable in international guidelines,9 this global metric may not accurately reflect regional morphological changes indicative of progression of fibrosis. In fact a recent large study suggests that FVC decline in IPF patients is heterogeneous and may have little correlation with increase in fibrosis in some patients.10

The Computer-Aided Lung Informatics for Pathology Evaluation and Rating algorithm (CALIPER) is an automated quantitative CT application which characterises lung parenchymal features on volumetric high-resolution CT (HRCT), thus quantifying abnormal lung.11 In IPF studies it is highly predictive of mortality,12 but also enhances risk stratification and cohort enrichment for trial end points.13 14 Until now, it has not yet been used to explore leucocyte association with progression of fibrosis.

The aim of the study is twofold—(1) to assess if leucocyte levels are linked to progression in fibrosis per se and (2) to explore which leucocytes demonstrate greatest association with fibrosis. We used CALIPER to provide a standardised automated, method of quantifying lung fibrosis.11

Methods

Study design and patients characteristics

We performed a retrospective analysis of a cohort of patients with IPF, who presented to the Oxford Interstitial Lung Disease (ILD) Service between September 2016 and November 2021. All patients with a multidisciplinary team (MDT) diagnosis of IPF,15 and a CALIPER-compatible HRCT (non-contrast, supine, volumetric HRCT scan) were included.

We asked two main questions—(1) what is the correlation between contemporaneous blood leucocyte levels (within 4 months of the first CT) and baseline amount of lung fibrosis and (2) which blood leucocyte levels (neutrophils, lymphocytes, monocytes and their derived ratios) correlated with progression in amount of lung fibrosis. Association between blood leucocytes and FVC decline and mortality were also examined as comparator end points.

Blood leucocyte measurement

Neutrophil, lymphocyte and monocyte levels captured from standard clinical ‘full blood count’ analysis within 4 months of initial CT (CT1). NLR, monocyte:lymphocyte ratio (MLR) and Systemic Inflammation Response Index (SIRI) (neutrophil×monocyte/lymphocyte) were calculated as described previously.5

CT scans

Non-contrast, supine, volumetric HRCT scan were acquired using a 64-detector row CT scanner. Images were reconstructed using a high spatial resolution algorithm. Non-contrast, volumetric, HRCT scans for appropriate subjects were acquired (0.625 mm slice thickness at an interval of 0.625 mm).

Computer-Aided Lung Informatics for Pathology Evaluation and Rating

Full CALIPER data acquisition and processing are described in online supplemental methods. Briefly, classification of parenchymal features was applied to 15×15×15 voxel volumes of interest (VOI) using texture analysis and computer-based algorithmic interpretation of volumetric histogram signature mapping features.11 CALIPER evaluation included characterisation and quantification of each VOI into one of five radiological parenchymal categories: normal lung, hyperlucent, ground glass opacity (GGO), reticular opacity and honeycombing. Parenchymal features were expressed as a relative percentage of CALIPER-derived total lung volume. Total lung fibrosis (TLF) represented the sum of GGO, reticular and honeycomb percentages.14 We included GGO as these had been reported by our radiologists to be found in areas of reticulation or traction bronchiectasis at MDTs and designated as fine fibrosis in keeping with agreed radiological assessments.16

Assessment of pulmonary vessel volume (PVV) was also performed prior to pulmonary vessel extraction from lung parenchyma using a multiscale tubular structure enhancement filter.17 This assessment excluded large vessels at the hilum of the lung. The PVV was calculated as absolute volume (cm3) and expressed as percentage of total CALIPER-derived lung volume.

The key variables were TLF (the sum of reticulation, GGO and honeycomb, as % of CALIPER-derived lung volume), TLF/litre (TLF expressed as % of lung volume) and PVV. TLF/litre was used to account for difference in lung volume on serial scans with progression in fibrosis.18 %change in fibrosis is (CT2 TLF/litre–CT1 TLF/litre)/CT1 TLF/litre×100.

Lung function tests

Baseline spirometry (forced expiratory volume in 1 s and FVC), transfer factor of lung for carbon monoxide (TLco) and Composite Physiological Index (CPI)19 within 3 months of CT were recorded; all as % predicted for age, height and sex. Annualised change in lung function was calculated as the relative (%) change in absolute value divided by time (years) between tests.

Statistical analysis

Pearson correlation was used to explore association between blood leucocytes with baseline CALIPER variables and pulmonary function tests. Cox proportional hazard modelling was employed to test association between each blood leucocyte level (monocytes, lymphocytes and neutrophils) and outcomes of disease progression—(1)≥7.8% increase in TLF/litre (7.8% being the lower limit of upper tertile) (model A), (2) ≥10% increase in TLF/litre (10% as an arbitrary value) (model B) and (3) ≥10% FVC decline between CT1 and CT2 (model C). Median and upper quartiles for model A were also tested in preliminary modelling but no significant correlation was found. Harrell’s Concordance index (C-index) was used to assess model strength, describing how well a model can discriminate between two survival distributions.20 Kaplan-Meier analysis (log rank test) was used to evaluate time to all-cause mortality from first CT. A censoring time of 1 January 2022 was applied.

All analyses were performed using GraphPad Prism V.9.2 (La Jolla, California, USA) or SPSS V.27 (IBM).

Results

Patients

A total of 171 eligible patients were identified. Of these, n=71 had at least one further follow-on HRCT (the CT with the longest interval from the first was selected). Demographics are shown in table 1. Median interval between CT1 and CT2 was 25.9 (16.8–39.9) months. Median time from first CT to death was 26.2 months (12.4–43.4) vs 45.7 (34.6–61.5) in those surviving to censoring date.

Table 1

Baseline characteristics for patients, at the point of first CT used for CALIPER analysis.

Association between blood leucocyte levels and baseline CALIPER and lung function measures

At baseline (CT1, n=171), there were significant correlations between baseline % predicted FVC, TLCO and CPI, and TLF, TLF/litre and PVV (table 2). In terms of correlation between leucocytes and CALIPER variables and lung function, there was significant correlation between neutrophil levels and TLF/litre (r=0.208, p=0.007), total PVV (r=0.259, p=0.001) and FVC (r=−0.127, p=0.029) (table 3). No other leucocyte levels or their derived measures showed significant corelations with CALIPER or lung function at baseline (table 3 and online supplemental table S1).

Table 2

Association between lung function and CALIPER parameters (Pearson’s correlation) at baseline; n=171 patients.

Table 3

Correlation of blood leucocytes with CALIPER parameters and lung function (n=171 patients, Pearson’s correlation).

Association between blood leucocyte levels and progression of fibrosis

Neutrophil count was significantly higher in cases demonstrating progression of fibrosis as defined by ∆TLF>10%/litre (4.87 (3.86–5.66) vs 4.84 (3.85–6.45), p=0.043). No other differences in leucocyte measure between cases with disease progression and stability was identified (online supplemental table 2).

In model A of Cox proportional hazard analysis, neutrophil levels (HR 1.81, 95% CI 1.10 to 2.99, p=0.020) were significantly associated with increase in TLF/litre≥7.8% (table 4). Lower lymphocyte levels were also significantly associated with TLF/litre≥7.8% (HR 0.26, 95% CI 0.08 to 0.91, p=0.034). Similar findings were observed in model B (where outcome was TLF/litre ≥10%), except that lower lymphocyte levels were not significantly associated with progression of >10%/Litre. In both models A and B, there was a trend of higher monocyte levels with progression of fibrosis, but this was not statistically significant.

Table 4

Cox proportional hazard analysis for progression of lung fibrosis.

In model C (where outcome was change by 10% of FVC), no leucocyte variable significantly associated with FVC decline. A lower baseline FVC was predictive of FVC decline (HR 0.96, 95% CI 0.92 to 0.99, p=0.049).

All leucocyte derived indexes were significantly associated with progression in fibrosis (TLF/litre by 7.8% and 10%) in multivariate analysis (online supplemental table 3).

50.7% of patients who underwent repeat CT were receiving antifibrotics at first CT. When adjusted for antifibrotic use and comorbidities, significance was preserved for neutrophils, NLR and SIRI (online supplemental tables S4,S5).

Association between blood leucocyte levels and mortality

During the study period, 60 all-cause deaths (35.1%) were reported in this cohort. Leucocyte levels were dichotomised by median values (figure 1) or normal reference range limits (online supplemental figure S1). Significantly shorter survival times were observed for cases dichotomised by median monocyte count (p=0.033) and neutrophil levels (p=0.0180) (figure 1), and for cases dichotomised by median NLR, MLR and SIRI values (online supplemental figure S1). Shorter survival times were also observed for those with TLF, PVV and FVC greater than the median value (TLF>3.30%/L, p<0.001; PVV>4.02%, p<0.001 and FVC>76.6%, p<0.001) (figure 1).

Figure 1

Kaplan-Meier curves for time to mortality for monocyte, neutrophil, lymphocyte levels, TLF, PVV and FVC, all at baseline. Leucocyte levels and CALIPER variables were dichotomised by median value. CALIPER, Computer-Aided Lung Informatics for Pathology Evaluation and Rating; FVC, forced vital capacity; PVV, pulmonary vessel volume; TLF, total lung fibrosis.

Discussion

In this study, we showed that when progression of disease is categorised specifically by increase in amount of fibrosis, quantified by an automated quantitative scoring modality, neutrophil levels rather than monocytes were the key immune correlate with progression in amount of fibrosis in the lungs. MLR, NLR and SIRI were also associated with progression. Although monocyte count was not associated with disease progression in this cohort, importantly, as with other studies, higher monocyte levels in our cohort were associated with mortality.4 21 Our findings suggest that neutrophils may be a contributor to active accumulation of fibrosis in IPF, supported by previous findings of high neutrophils in bronchoalveolar lavage (BAL) of IPF patients,22 and more recent findings of higher NLR in IPF patients with greater rates of FVC decline.5

Neutrophils have long been associated with immunopathogenesis of IPF. Neutrophilia in BAL specimens from IPF patients is associated with earlier mortality.22 Neutrophil elastases (NE) are elevated in IPF BAL samples,23 and experimental data using murine models suggest that NE activates transforming growth factor-ß pathway and fibroblast proliferation.24 An intriguing and newly identified fibrosis-promoting function of neutrophils is generation of neutrophil extracellular traps (NETs).25 These proinflammatory collections of chromatin and neutrophils regulate both immune cell function and fibroblast activation.26 While a specific association with IPF has yet to be fully described, enhanced detection of intrapulmonary NETs has been reported in bleomycin models and in non-IPF fibrotic ILD studies.27 28 Here, we are able to link blood neutrophil level measured at baseline CT specifically with progression in amount of fibrosis over time. Further prospective translational studies are required to establish the regulatory role of neutrophils, NET formation over time and cytokine and chemokine activity at different stages of fibrosis.

In addition to TLF change, PVV measures also correlated with neutrophil count. PVV is an intriguing quantitative variable that has proved a consistent predictor of disease severity, progression and mortality in ILD studies.12 14 29–31 Our findings are supportive of the original work performed by Jacob et al who demonstrated that PVV% was independently associated with FVC decline,32 and Chung et al who demonstrated PVV negatively correlated with TLCO%.33 In our study, PVV was independently associated with progression of fibrosis. The pathological mechanisms linking PVV with adverse ILD outcomes are not fully understood and several theories have been suggested. Blood perfusion is reduced in areas of pulmonary fibrosis,34 but increased in adjacent areas of unaffected lung.35 Jacob et al postulated that correlation between ILD extent and vessel calibre may represent regional elevation in pulmonary artery pressures in mildly fibrotic lung or destruction of the capillary bed in more advanced disease leading to neovascularisation and diversion of blood to unaffected lung areas.12 Another possible explanation relates to the negative intrathoracic pressures required of non-compliant fibrotic lungs to generate adequate inspiratory volumes. This could in turn exert additional ‘tractional’ force on the lung vasculature resulting in dilatation in fibrotic regions in a traction-like phenomenon, akin to traction bronchiectasis.33

The trend of lower lymphocytes count with progression of fibrosis mirrors the findings of previous studies.5 36 Lymphocytic aggregates are a recognised pathological feature of IPF lesions.27 37 The association between low blood lymphocyte count and adverse outcomes in IPF is currently unknown but could be explained in part by lymphocyte dysfunction,38 and the sequestering of lymphocytes into sites of inflammation, such as the fibrotic lung. We note that post hoc analysis of the ASCEND and CAPACITY studies Nathan et al reported that serial increase in NLR over 12 months was associated with mortality.39

In our study, we report association between higher neutrophil count and lower lymphocyte count taken from full blood count analysis with progression of fibrosis scoring on HRCT. To our knowledge, this specific association has not been reported elsewhere. Other published studies and post hoc analyses have only reported leucocyte association with mortality, hospitalisation and FVC decline.

Our findings should be framed within the following limitations. First, the retrospective nature of the study meant that we selected patients who had repeat CT scans, which were clinically indicated. This introduced a selection bias towards patients who had a clinical reason to have a CT scan—often, this is worsening in disease either over time or acutely, and CT scans would be performed at different/non-uniform time intervals. Unfortunately, we did not have access to blood leucocyte data measured at the time of the second CT scan. Analysing blood leucocytes/ratios at second CT could may have provided additional information, especially with regard to how any relative change in blood count/ratio may differ between patients with progression of fibrosis or stability at second CT scans. Although the leucocyte levels were measured at the time of the first (presenting) CT scans, this group of patients could have a different disease trajectory compared with those who did not have a second CT scan. There is a possibility that neutrophil levels are more likely a driver in those who are clinically deteriorating. The baseline correlation data (n=171 patients), however, are not limited by these factors. Further analysis of prospective validation cohorts would be required to investigate whether neutrophil count, measured in peripheral blood, could identify patients at greater risk of progression of fibrosis in IPF, and consequently facilitate prognostication of the disease.

From the point of fibrosis scoring, the sum of GGO, reticulation and honeycombing was used as the TLF score.14 Although each scan was reported by a radiologist who ‘called’ the GGO as fine fibrosis and the overall pattern as probable or definite UIP, other causes such as pulmonary oedema, acute exacerbation and infection could have contributed to this CT finding.40

Nevertheless, our study shows that greater neutrophil count was significantly associated with quantifiable increase in fibrosis on imaging of the lungs in IPF. Further studies will be required to validate this finding.

Data availability statement

Data are available on reasonable request. N/A.

Ethics statements

Patient consent for publication

Ethics approval

The study was part of a study examining factors associated with progression of IPF (ethical approval 14/SC/1060 from the Health Research Authority and South-Central National Research Ethics Service).

References

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • Contributors AA conceived project, conducted analysis, interpreted data and wrote the paper. EF, PS and RKH provided guidance with project. RB provided guidance with project and assisted with defining search criteria. L-PH conceived project, interpreted data, wrote paper and supervised the study. L-PH accepts full responsibility for the work and/or the conduct of the study, had access to the data, and controlled the decision to publish.

  • Funding The study was funded by the National Institute for Health Research (NIHR) Oxford Biomedical Research Centre (BRC). L-PH is supported in part by MRC UK (MC_UU_00008/1). Correspondence to L-PH (ling-pei.ho@imm.ox.ac.uk).

  • Competing interests AA, EF, RB, L-PH and RKH report no relevant conflict of interests. PS has received consultancy fees from Trevi Therapeutics, and lecture fees from Boehringer Ingelheim, but no other conflict of interest.

  • Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.