Abstract
Objectives
The peak expiratory flow rate (PEFR) is known to decrease in patients with sarcopenia. However, little is known about the clinical impact of the PEFR in idiopathic pulmonary fibrosis (IPF). This study aimed to confirm whether a decrease in PEFR over 6 months was associated with survival in IPF patients.
Methods
Consecutive IPF patients who had been assessed at a single center were retrospectively analyzed. The relative decline in PEFR over 6 months was assessed. Survival analyses were performed by univariate and multivariate Cox proportional hazard models.
Results
A total of 61 eligible cases (average age 70 years) were examined, and 21 patients (34.4%) died. The univariate Cox regression analysis showed that the body mass index, baseline % predicted forced vital capacity (FVC), baseline % predicted PEFR, % predicted diffusion capacity for carbon monoxide (DLCO), relative decline in FVC, and relative decline in PEFR were prognostic factors. On multivariate analyses, relative decline in PEFR (hazard ratio [HR] 1.037, p < .05) and baseline % predicted FVC (HR 0.932, p < .001) were independent prognostic factors, whereas relative decline in FVC was not.
Conclusion
A decrease in PEFR after 6 months may predict worse survival in patients with IPF.
Keywords
Introduction
Idiopathic pulmonary fibrosis (IPF) is a fibrotic and progressive pulmonary disease that leads to death in most patients. 1 The median survival times of IPF patients were poor at 3 years before the development of anti-fibrotic drugs.1–3 This disease is characterized by decreased lung volumes and reduced gas exchange, and it is associated with symptoms of progressive dyspnea, cough, and reduced exercise capacity. Several known prognostic factors for IPF are: dyspnea score, 4 pulmonary function,5–7 disease severity,2,8 functional exercise capacity,9,10 and fibrotic changes on high-resolution computed tomography (HRCT). 11 Of pulmonary function variables, a decline in forced vital capacity (FVC) is widely known as a useful prognostic factor.5–7
Sarcopenia is a syndrome characterized by a progressive and generalized loss of systemic skeletal muscle mass and strength, and it carries the risk of poor outcomes such as physical disability, poor quality of life, and death. 12 It has been reported that the loss of skeletal muscle mass assessed using computed tomography (CT) images is associated with a poor prognosis in patients with IPF.13,14 Furthermore, we have previously reported that the decrease over the course of 6 months of the cross-sectional area (CSA) of the erector spinae muscle (ESMCSA), evaluated using CT images, was a strong predictor of mortality in patients with IPF. 15 In that study, baseline ESMCSA was not a prognostic factor on multivariate analysis.
The peak expiratory flow rate (PEFR) is obtained by routine spirometry. PEFR is the maximum flow rate generated during a forceful exhalation, starting from full lung inflation. PEFR well reflects large airway flow in patients with asthma and chronic obstructive pulmonary disease, but in patients with IPF, it would strongly depend on the respiratory muscular strength of the patient. Generally, in patients with IPF, radial traction due to scarring of the surrounding parenchyma prevents dynamic compression of the airway and, subsequently, does not reduce PEFR. 16
In previous studies, a decrease in PEFR was correlated with sarcopenia and skeletal muscle mass loss.17,18 Therefore, we hypothesized that a decrease in PEFR over the course of 6 months would be associated with poor survival in patients with IPF.
Methods
This study protocol was approved by the Ethics Committee of Nagoya City University Hospital (approval number 60-21-055) and carried out in accordance with approved guidelines. The need for patient approval and/or informed consent was waived by the Ethics Committee of Nagoya City University Hospital due to the retrospective nature of the study.
Patients
The clinical records of 80 consecutive patients with IPF who were referred to the Nagoya City University Hospital between October 2013 and April 2020 were retrospectively reviewed. IPF was diagnosed through multidisciplinary discussions according to international guidelines.1,19,20 The exclusion criteria were as follows: patients who died sooner than 6 months from baseline; patients whose pulmonary function was not assessed at the initial visit or after 6 months due to any cause; and patients who had lung cancer at baseline. The observation period commenced from the time of the patient’s initial evaluation until death, or 30 April 2021, whichever was sooner. Censored cases were evaluated to confirm their life or death status by telephone, as often as was possible. A pulmonary rehabilitation program was performed during hospitalization by patients who had acute exacerbations of IPF or pneumonia. Outpatients did not perform rehabilitation.
Pulmonary function tests
All patients completed pulmonary function tests using spirometry (CHESTAC-8900; Chest, Tokyo, Japan) according to the ATS/ERS criteria. 21 The diffusion capacity for carbon monoxide (DLco) was also measured (CHESTAC-8900). Percentage of predicted FVC (%FVC), percentage of predicted forced expiratory volume in 1.0 s (FEV1), and percentage of predicted DLco (%DLco) were calculated based on the patients’ height, age, and sex, according to the Japanese standardized methods. 22 Percent (%) predicted PEFR was calculated based on the prediction equation reported by Nunn et al. 23 The percent (%) predicted maximum mid-expiratory flow rate (MMF) was calculated using the prediction equation reported by Schmidt et al. 24
Derivation of ESMCSA by imaging analysis software
SYNAPSE VINCENT (Fujifilm Medical Systems, Tokyo, Japan) CT imaging analysis software was used for the derivation of ESMCSA. ESMCSA was calculated manually according to a previously published method.15,25 Briefly, ESMCSA was measured on a single-slice axial CT image at the level of the spinous process of the 12th thoracic vertebra. For the quantitative analysis of the ESMs, chest HRCT images were reconstructed using the mediastinal window setting (window level, 40 HU; window width, 300 HU). The left and right ESMs were identified and manually shaded, and the ESM area was reported as the sum of the right and left ESMs.
The Charlson comorbidity index
The Charlson comorbidity index was calculated according to a previously reported method. 26
Statistical analyses
Continuous variables are presented as means (±standard deviations) or medians [interquartile ranges]. Categorical variables are presented as numbers and percentages. The relationships between continuous variables were evaluated using Spearman’s rank correlation coefficients. Univariate Cox proportional hazard models were used to examine the associations of selected variables with survival. Stepwise multivariate Cox proportional regression analysis was performed including the significant variables on univariate analysis, with subsequent forward stepwise selection. To avoid multicollinearity, some of the highly correlated variables were excluded on multivariate analysis if they had Spearman’s correlation coefficients greater than 0.7. Survival times were estimated using the Kaplan–Meier method and compared using the log-rank test. To determine the optimal cutoff value for predicting 2-years mortality, receiver operating characteristic (ROC) curves were constructed. Multiple linear regression analysis was performed to identify clinical parameters related to relative decline in PEFR. To assess multicollinearity, variance inflation factor (VIF) values were computed. The differences between patients separated by cutoff values of the relative decline in PEFR were analyzed using Student’s t-test or the Mann-Whitney U test for continuous variables and the Chi-squared test or Fisher’s exact test for categorical variables, as appropriate. p values < .05 were considered significant. Statistical analyses were conducted using IBM SPSS Statistics version 28 (IBM Corp., Armonk, NY, USA).
Results
Patients’ characteristics
Patients’ characteristics.
Data are presented as means (± standard deviation), medians [interquartile range], or numbers (%).
FVC: forced vital capacity; PEFR: peak expiratory flow rate; FEV1: forced expiratory volume in 1.0 s; MMF: maximum mid-expiratory flow rate; DLCO: diffusion capacity of the lung for carbon monoxide; ESMCSA; the cross-sectional area of the erector spinae muscle evaluated using computed tomography images.
aTwo patients could not perform the DLCO maneuver.
bNine cases could not be evaluated.
Correlation between PEFR and ESMCSA
It has been reported that PEFR is an indicator of sarcopenia.
17
On the other hand, assessment of the ESMCSA from chest CT images has been used to evaluate sarcopenia in patients with chronic lung disease.13,25 Therefore, the relationship between PEFR and ESMCSA was examined in patients with IPF. The mean ESMCSA was 29.1 ± 6.4 cm2 in the present cohort. Nine patients' ESMCSA values could not be analyzed because their CT lacked the level of the spinous process of the 12th thoracic vertebra. Baseline PEFR was significantly correlated with ESMCSA (r = 0.431, p < .005) (Figure 1). Correlation between PEFR and ESMCSA. The correlation between baseline PEFR and ESMCSA is shown. PEFR: peak expiratory flow rate; ESMCSA: the cross-sectional area of the erector spinae muscle evaluated on computed tomography images.
Relative decline in PEFR and other respiratory function parameters
The median relative decline in PEFR was 2.90 (−7.27–11.54) %. In addition, the median relative decline in FVC was 2.45 (−2.49–5.54) %, the median relative decline in FEV1 was 3.19 (−2.40–9.03) %, and the median relative decline in MMF was 4.44 (−12.81–15.74) %.
Prognostic survey
Twenty-one of 61 patients (34.4%) died during the study period. The follow-up time of the 61 patients was 1017 days (530–1623). Four patients were lost to follow-up. The median survival estimate of the cohort was 2214 days.
Prognostic factors on univariate and multivariate Cox regression analyses
Prediction of mortality by univariate and multivariate Cox proportional hazards analyses.
HR: hazard ratio; CI: confidence interval; FVC: forced vital capacity; PEFR: peak expiratory flow rate; FEV1: forced expiratory volume in 1.0 s; MMF: maximum mid-expiratory flow rate; DLCO: diffusion capacity of the lung for carbon monoxide; ESMCSA: the cross-sectional area of the erector spinae muscle evaluated using computed tomography images.
HRs and 95% CIs for each parameter on univariate and multivariate Cox regression analyses of only baseline data are shown in Supplemental Table S1.
Kaplan–Meier curves and a log-rank test
The Kaplan–Meier curves of this cohort are shown in Figure 2. To determine the optimal cutoff value for predicting 2-years mortality, an ROC analysis was performed by using the data of 47 patients. A cutoff value of 8.65% (area under the curve = 0.737, specificity: 0.778, sensitivity: 0.763) was identified. The mean survival estimates were as follows: relative decline in PEFR ≥8.65%, 1335 days, and relative decline in PEFR <8.65%, 2169 days. Relative decline in PEFR ≥8.65% had a significantly worse prognosis (p < .01, log-rank test). Kaplan–Meier curves and a log-rank test. Kaplan–Meier survival curves stratified by relative decline in PEFR. The cutoff value was set at 8.65%. Patients with a relative decline in PEFR over 8.65% had significantly worse survival (p < .01, log-rank test). PEFR: peak expiratory flow rate.
Correlations between the relative decline in PEFR and other clinical parameters
Correlations between relative decline in PEFR and clinical parameters.
CI: confidence interval; FVC: forced vital capacity; PEFR: peak expiratory flow rate; FEV1: forced expiratory volume in 1.0 s; MMF: maximum mid-expiratory flow rate; DLCO: diffusion capacity of the lung for carbon monoxide; ESMCSA: the cross-sectional area of the erector spinae muscle evaluated using CT images.
Multiple linear regression analysis for decline in PEFR
Multiple linear regression analysis for decline in PEFR.
CI: confidence interval; VIF: variance inflation factor values; FVC: forced vital capacity; PEFR: peak expiratory flow rate; MMF: maximum mid-expiratory flow rate; DLCO: diffusion capacity of the lung for carbon monoxide.
Characteristics of IPF patients separated by the cutoff value
The clinical characteristics of the IPF patients classified by the cutoff value of 8.65% for relative decline in PEFR are shown in Supplemental Table S2. There was a significant difference in baseline %FVC (p < .05) and baseline %FEV1 (p < .05) between the two groups. In addition, the clinical characteristics with and without anti-fibrotic drug use are shown in Supplemental Table S3.
Discussion
Several studies have reported that a decrease in PEFR was an indicator of sarcopenia and skeletal muscle mass loss in older adults.17,18 Kera et al. 27 created a definition of respiratory sarcopenia based on the PEFR (cut-off values were 4.40 L/s for men, 3.21 L/s for women) and demonstrated that respiratory sarcopenia correlated with conventional sarcopenia. Evaluation of ESMCSA on chest CT images has also been used to assess sarcopenia in patients with chronic lung disease.13,25 Regarding the baseline ESMCSA, Suzuki et al. reported that baseline ESMCSA was related to prognosis, 13 whereas Nakano et al. 15 reported that it was not related to prognosis. In the present study, baseline % predicted PEFR was one of the prognostic factors on univariate Cox regression analysis, whereas baseline ESMCSA was not. Baseline PEFR was correlated with ESMCSA in patients with IPF (Figure 1). The measurement of PEFR seems to be simpler than the analysis of ESMCSA, and we considered it as an indicator of sarcopenia that predicts prognosis in patients with IPF.
The present study demonstrated that a decrease in PEFR over 6 months was a strong prognostic factor in IPF patients. PEFR can be obtained by routine spirometry. However, to the best of our knowledge, this is the first study to examine the use of PEFR for mortality analysis in patients with IPF. Interestingly, on multivariate analysis, relative decline in PEFR (HR 1.037, 95% CI 1.002–1.072, p < .05) was an independent prognostic factor, whereas the relative decline in FVC was not.
Two anti-fibrotic therapies, nintedanib and pirfenidone, have been approved for the treatment of IPF. Anti-fibrotic drugs have been shown to significantly reduce the annual rate of decrease in FVC, and they appear to reduce the risk of all-cause mortality. 28 In the present study, more than 80% of patients were treated with anti-fibrotic drugs. This may cause a reduction in the relative decline in FVC. Even if the effects of anti-fibrotic drugs are taken into account, the decrease in PEFR may have a strong impact on the survival rate in patients with IPF.
Regarding the loss of skeletal muscles in patients with IPF, a low ESMCSA on CT images13–15 and a low fat-free mass index by bioelectrical impedance analysis 29 were demonstrated to be associated with a poor prognosis. In patients with IPF, progressive loss of systemic skeletal muscles can be thought to decrease daily activity, worsen chronic respiratory failure, and increase infectious pneumonia that leads to death. 30
Multiple linear regression analysis showed that age, baseline %FVC, and baseline %PEFR were independent factors contributing to relative decline in PEFR. A decrease in PEFR can be caused by progressive weakness of the respiratory muscles.17,18,27 Loss of respiratory muscle mass may be induced by systemic inflammation, malnutrition, enhanced energy expenditure, and aging. More advanced and elderly IPF patients might have greater losses of systemic skeletal muscle mass and strength. It is possible that a rehabilitation program and nutritional guidance for outpatients with IPF could improve the PEFR decline.
In the present study, baseline % predicted PEFR was also a prognostic factor (HR 0.978, 95% CI 0.956–0.999, p < .05) on univariate analysis, but not on multivariate analysis. The relative decline in PEFR over 6 months was considered to be more useful in predicting survival than the baseline % predicted PEFR.
Recently, 3-month change in FVC has shown prognostic potential in a large IPF meta-analysis. 31 In the present study, the reason for choosing a 6-month period was that most patients had respiratory function tests performed every 6 months. Most patients visited our hospital every 2 months or every 3 months. It would be useful if short-term changes of markers were found to be prognostic factors. Investigation of daily home spirometry 32 might provide more useful markers.
The present study has the following limitations. First, the results were obtained by a retrospective analysis of all Japanese patients from a single center with a small sample size. Second, DLco and ESMCSA after 6 months were not analyzed because there were many cases of missing data. Third, there may be survivor bias and significant lead-time bias. Fourth, an ROC analysis was performed using data of only 47 patients (77% of the cohort) to determine the optimal cutoff values for predicting 2-years mortality. Future studies are necessary to examine the cutoff values in more facilities and in larger samples. The optimal cutoff value may be useful for physicians to note the possibility of weakened respiratory muscles, and to decide starting muscle training with rehabilitation and nutritional guidance. Thus, the present study can be considered a preliminary report. Further studies are warranted to elucidate the importance of relative decline of PEFR in patients with IPF.
In conclusion, a decrease in PEFR after 6 months could be a novel and useful prognostic factor in patients with IPF. Progression of respiratory sarcopenia can be considered the cause of the decrease in PEFR. It is important to pay attention to PEFR in the clinical practice of IPF.
Supplemental Material
Supplemental Material - Decreased peak expiratory flow rate associated with mortality in idiopathic pulmonary fibrosis: A preliminary report
Supplemental Material for Decreased peak expiratory flow rate associated with mortality in idiopathic pulmonary fibrosis: A preliminary report by Kohei Fujita, Hirotsugu Ohkubo, Akiko Nakano, Norihisa Takeda, Kensuke Fukumitsu, Satoshi Fukuda, Yoshihiro Kanemitsu, Takehiro Uemura, Tomoko Tajiri, Ken Maeno, Yutaka Ito, Tetsuya Oguri, Yoshiyuki Ozawa, Takayuki Murase and Akio Niimi in Chronic Respiratory Disease
Footnotes
Acknowledgements
Author contributions
KoF and HO contributed equally to this study. KoF and HO drafted the submitted article and take responsibility for the integrity of the data and the accuracy of the data analysis. ANa, NT, KeF, SF, YK, TU, TT, KM, YI, TO, and ANi contributed to the interpretation of the manuscript. YO contributed as a radiologist. TM contributed as a pathologist.
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
Supplemental Material
Supplemental Material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
