==== Front Biomed Res Int Biomed Res Int BMRI BioMed Research International 2314-6133 2314-6141 Hindawi 10.1155/2020/8848919 Research Article A Clinical Model for the Prediction of Acute Exacerbation Risk in Patients with Idiopathic Pulmonary Fibrosis https://orcid.org/0000-0002-0537-6467Wu Qi 1 https://orcid.org/0000-0003-3990-6013Xu Yong 2 https://orcid.org/0000-0002-6970-8571Zhang Ke-jia 3 4 https://orcid.org/0000-0002-1795-8587Jiang Shi-min 3 4 https://orcid.org/0000-0002-2845-6339Zhou Yao zhouyao1226@163.com 3 4 https://orcid.org/0000-0002-1098-2832Zhao Ying njmuzy@126.com 5 1Department of Physiology, Xuzhou Medical University, Xuzhou 221009, China 2Affiliated Hospital of Nanjing University of Chinese Medicine, Nanjing 210029, China 3Department of Pathophysiology, Xuzhou Medical University, Xuzhou 221009, China 4Laboratory of Clinical and Experimental Pathology, Xuzhou Medical University, Xuzhou 221009, China 5Xuzhou Traditional Chinese Medicine Hospital Affiliated to Nanjing University of Chinese Medicine, Xuzhou 221009, China Academic Editor: Bing Niu 2020 7 12 2020 2020 884891917 9 2020 26 11 2020 30 11 2020 Copyright © 2020 Qi Wu et al.2020This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.Objective To develop and validate a risk assessment model for the prediction of the acute exacerbation of idiopathic pulmonary fibrosis (AE-IPF) in patients with idiopathic pulmonary fibrosis (IPF). Methods We enrolled a total of 110 patients with IPF, hospitalized or treated as outpatients at Xuzhou Traditional Chinese Medicine Hospital Affiliated to Nanjing University of Chinese Medicine from July 2012 to July 2020. Of these, 78 and 32 patients were randomly assigned to training and test groups, respectively. The risk factors for AE-IPF were analyzed using logistic regression analysis, and a nomographic model was constructed. The accuracy, degree of calibration, and clinical usefulness of the model were assessed with the consistency index (C-index), calibration diagram, and decision curve analysis (DCA). Finally, the stability of the model was tested using internal validation. Results The results of logistic regression analysis showed that a history of occupational exposure, diabetes mellitus (DM), essential hypertension (EH), and diffusion capacity for carbon monoxide (DLCO)% predicted were independent risk factors for AE-IPF prediction. The nomographic model was constructed based on these independent risk factors, and the C-index was 0.80. The C-index for the internal validation was 0.75, suggesting that the model had good accuracy. The decision curve indicated that for a threshold value of 0.04–0.66, greater clinical benefit was obtained with the AE-IPF risk prediction model. Conclusion A customized AE-IPF prediction model based on a history of occupational exposure, DM, EH, and DLCO% predicted provided a reference for the clinical prediction of AE-IPF. Xuzhou Medical UniversityD2019022D2019005Science and Technology Program of Xuzhou Science and Technology BureauKC17113 ==== Body 1. Introduction Idiopathic pulmonary fibrosis (IPF) is a chronic progressive interstitial lung disease characterized by severe pulmonary tissue fibrosis without a known pathogenesis; the mean survival time after diagnosis is 2–3 years[1]. The occurrence and development of IPF is related to epithelial cell injury, the proliferation of fibroblasts, inflammatory response, and extracellular matrix deposition [2, 3]. The pathologic feature of IPF is usual interstitial pneumonia (UIP). Unexpected respiratory failure, known clinically as the acute exacerbation of IPF (AE-IPF), usually occurs in a proportion of patients with IPF [4–6]. AE-IPF was first reported in 1984. Yoshimura et al. briefly reported the symptoms and mortality rate of AE-IPF and identified the efficacy of a large dosage of glucocorticoids [7]. It was not until 2002 that the terminology of AE-IPF was officially acknowledged in the clinical consensus of interstitial lung diseases. Although multiple retrospective studies have reported on the incidence and mortality rate of AE-IPF, a consensus has not been reached. This can be attributed to the different study designs, inconsistent standards for disease diagnosis, few cases, and loss of patients to follow-up [8]. In a retrospective study of 461 IPF cases, AE-IPF was identified as the major cause for the deterioration of IPF, accounting for 55.2% of the cases. In another study, the incidence of AE-IPF was 14.2% and 20.7% at 1 and 3 years after IPF, respectively [9]. The median overall survival of AE-IPF patients was only 2.2 months; of these patients, 50% died during hospitalization and 40% died in the ICU (mortality rate > 90%) [9]. Currently, there are no effective therapies for the treatment of AE-IPF. Internationally acknowledged suggestions include supportive nursing and glucocorticoid treatment which aim to alleviate symptoms and rectify hyperemia; however, the effects of these strategies are limited [10]. Therefore, it is pertinent to construct a risk assessment model to evaluate the risk factors for AE-IPF in patients with IPF. A wealth of literature has documented that inheritable factors, air pollution, seasonal factors, partial deterioration in lung function (assessed via forced vital capacity (FVC), diffusion capacity for carbon monoxide (DLCO), 6 min walking test, and St. George's respiratory questionnaire), pulmonary arterial hypertension, and coronary artery disease may all increase the risk for AE-IPF [11]. However, a model for the prediction of AE-IPF is lacking. Consequently, in this study, multivariate logistic regression analysis was used to identify high-risk factors for AE-IPF, in order to analyze the degree of risk, construct a simple prediction model, and provide an effective tool for the prediction of AE-IPF. 2. Patients and Methods 2.1. Patients Patients with IPF who were hospitalized in or visited Xuzhou Traditional Chinese Medicine Hospital Affiliated to Nanjing University of Chinese Medicine from July 2012 to July 2020 were enrolled in the present study. In total, 110 patients were included; 78 patients were selected for the construction of the risk analysis model, and the remaining 32 patients were used to test the model. For patients with multiple hospitalizations or visits, the data from the first hospitalization or visit were adopted. The patients provided written informed consent for inclusion in this study, and the study was approved by the Ethics Commission of Xuzhou Hospital of Traditional Chinese Medicine (ID: XZTCM2018LSY-013). 2.2. Diagnostic Criteria Patients with IPF were confirmed based on clinical features, imaging data, and medical history according to the diagnostic criteria for IPF (2011) which included the following: (1) interstitial pneumonia with known etiology (family or occupational exposure, desmosis, and drug dependency) and (2) HRCT-indicated UIP or suspected UIP symptoms [12]. The diagnostic criteria for AE-IPF were as follows: (1) past or present IPF, (2) development or occurrence of deteriorative acute dyspnea within 1 month, (3) newly formed bilateral frosted glass shadow or consolidation shadow at the original reticular shadow or beehive-like shadow UIP via thoracic HRCT, and (4) deterioration of dyspnea that could not be explained by heart failure or liquid overload [13]. 2.3. Data Collection Clinical data were obtained from the medical history of the patients during hospitalization and the follow-up visits. The data included the following: (1) patient characteristics, such as age, sex, occupational history, and history of smoking; (2) complications, such as diabetes mellitus (DM), essential hypertension (EH), coronary heart disease (CHD), cerebral infarction (CI), hypothyroidism, gastroesophageal reflux disease (GERD), and obstructive sleep apnea-hypopnea syndrome (OSAHS); (3) lung function test indices, such as Forced Vital Capacity (FVC)% predicted and diffusing capacity for carbon monoxide (DLCO)% predicted; and (4) medical treatment, including N-acetylcysteine (NAC), pirfenidone (PFD), glucocorticoid (GC), and acid-inhibitory drugs. 2.4. Statistical Analysis All clinical data are shown as count (%) and were analyzed with R version 3.6.3 (R Foundation for Statistical Computing, Vienna, Austria). The Caret package in R was used to randomly assign the 110 patients to the training and test groups in a ratio of 7 : 3 (78 cases and 32 cases, respectively). Information on the training group was subjected to multivariate logistic regression with the odds ratio, 95% confidence interval, andpvalue as characteristics; a p value of <0.05 was considered statistically significant. Factors describing statistical significance were used to construct the nomogram using the “rms” program to create a prediction model. To test the accuracy of the prediction model, we used the “rms” program to calculate the consistency index (C-index). The degree of calibration was evaluated by constructing a calibration curve; the closer the calibration curve of the prediction model to the standard curve, the better the consistency of the prediction model. The area under curve (AUC) value of the model was calculated to predict the degree of discrimination of the prediction model, with a higher AUC value representing a higher degree of discrimination. By quantifying the net benefits of the probability of different threshold values in the array, decision curve analysis (DCA) was used to determine the clinical usefulness of the nomogram. Finally, an internal validation method was used to test the stability of the prediction model; the “rms” program package was used to calculate the C-index and AUC value for the 32 patients with IPF in the test group. 3. Results 3.1. Patient Characteristics In total, 110 patients with IPF, treated as inpatients or outpatients between July 2012 and July 2020, were enrolled in the present study. In the training group (78 cases), there were 16 patients with AE-IPF (20.51%), while in the test group (32 cases), there were four patients with AE-IPF (12.50%). All the patients' clinical data, including demographic and clinical characteristics, complications, lung function, and medication, are presented in Table 1. 3.2. Analysis of Risk Factors for AE-IPF The clinical data of the training group were entered into the multivariate logistic regression model to obtain the coefficients of the corresponding characteristic variables. Among the variables, a history of occupational exposure (p < 0.01), DM (p < 0.01), EH (p < 0.01), and DLCO%-predicted lung function (p < 0.05) were identified as the risk factors for AE-IPF (Table 2). 3.3. Nomography of AE-IPF Risk Based on the results of the logistic regression, we determined the predictive factors for AE-IPF. To facilitate the estimation of the risk for AE-IPF, we screened the four characteristic variables and then constructed a nomographic risk prediction model. Each variable was scored on a scale, and the range of the total score was 0–220. The total score on the AE-IPF risk axis represents the probability of AE-IPF. The higher the score, the higher the risk that a patient with IPF will develop AE-IPF (Figure 1). 3.4. Assessment of the Degree of Calibration for the Risk Prediction Model of AE-IPF We constructed calibration curves to facilitate the assessment of the degree of calibration of the risk prediction model of AE-IPF. As shown in Figure 2, the x-axis represents the predicted risk for AE-IPF, and the y-axis represents the realization of AE-IPF. The bidiagonal dotted line represents the prediction of an ideal model, and the solid line represents the realized prediction capacity. The closer the dotted line is to the bidiagonal line, the greater the prediction capacity. The results of the present study showed a preferable consistency between the nomographic model and the ideal model. 3.5. Assessment of the Accuracy of AE-IPF Risk Prediction To test the accuracy of the model, we first calculated the C-index. An ROC curve indicates that the predictions are close to random guesses. Meanwhile, when an ROC curve is plotted, we can use AUC to measure the performance of the predictor. We then calculated the AUC value of the model and constructed an ROC curve [12]. Based on the calculation, the C-index was 0.80, suggesting good accuracy of the model. As shown in Figure 3, the AUC of the ROC curve was 0.77, suggesting good accuracy of the AE-IPF risk prediction model. 3.6. Internal Validation of the Model for AE-IPF Risk Prediction To further evaluate the stability of the AE-IPF risk prediction model, we calculated the C-index and AUC values of the 32 patients with IPF in the test group using the internal validation group. The C-index value was 0.75, and after constructing the ROC curve to calculate the AUC value for further validation, the AUC value was 0.70 (Figure 4). Therefore, both the C-index and AUC confirmed the good stability of the AE-IPF risk prediction model. 3.7. Clinical Net Benefit DCA was used to evaluate whether the prediction model could improve clinical decision-making. The decision curve of the prediction model is shown in Figure 5. The y-axis indicates the net benefit, and the x-axis represents the probability of the threshold value. The black solid line represents no intervention, at which the net benefit is zero. The black dotted line represents the intervention, and the net benefit is an oblique line with a negative slope. The blue solid line represents the realized profits of the AE-IPF risk prediction model. The decision curve showed that the greatest clinical benefit will be obtained from the AE-IPF risk prediction model at threshold values of 0.04–0.66. 4. Discussion At present, the cause of AE-IPF is unknown. Recurrent local outbreaks of alveolar damage may lead to the progression of IPF [14]. For some patients, recurrent outbreaks of alveolar damage will inflict damage and worsen their condition, eventually leading to acute exacerbation [15]. In addition, it has been suggested that infectious factors and some direct lung injury can also induce the pathophysiological processes of AE-IPF [16]. AE-IPF can occur at any stage during the process of IPF, and AE-IPF may even occur during the relatively stable stage of IPF [17]. Therefore, it is important to predict the occurrence of AE-IPF. However, a reliable model for the prediction of AE-IPF is lacking. In the present study, a model was created to predict the risk of AE-IPF and to provide a reference for the clinical prediction of AE-IPF. Based on the prediction model, occupational exposure is a major risk factor for AE-IPF. Therefore, workers engaged in ore mining, construction work, inorganic dust production, and chemical processing and manufacturing, and those exposed to organic dust production, are at high risk. Owing to a lack of research, the association between occupational exposure and a strong predisposition toward AE-IPF is unknown. It is likely that when productive chemicals with systemic toxicity enter the respiratory tract in the form of aerosols or gas, they readily enter the blood flow through the vascular bed in the lung, leading to respiratory system damage, systemic poisoning, and AE-IPF. In addition to IPF itself, some patients with IPF also have other complications, which may also induce AE-IPF [18]. According to the prediction model, DM and EH are high-risk factors for AE-IPF. It is suggested that DM is an independent factor for the occurrence of IPF [19]. Patients with DM usually experience aberrant lung function, mainly characterized by impairments in lung ventilation, malfunctions in lung diffusion, and abnormalities of bronchomuscular tension and respiratory muscle function [20]. It has also been reported that the malfunction of pulmonary ventilation in DM patients can become compounded over the course of the disease [21]. However, it is unclear why DM leads to IPF exacerbation. It is possible that the final metabolites of DM cause airway endothelial cell injury and abnormalities in the airway defense mechanism, which eventually leads to the progression of IPF. Regarding EH, the association between EH and AE-IPF has received less attention. However, according to the prediction model, EH is another risk factor for the occurrence of AE-IPF. This should be proven by a large-sized clinical trial in the future. Therefore, the treatment of IPF requires a systematic approach, including time discrimination and the treatment of DM, EH, and other complications, in order to optimize treatment efficiency. Lung function measurement is an important method for IPF diagnosis and for assessing the state of illness [22]. According to the consensus for treatment and diagnosis of IPF issued by ATS/ERS in 2011, a ≥10% reduction in FVC compared with the absolute value at baseline and a ≥15% reduction in DLCO% predicted suggest a good prognosis for IPF [23]. This highlights the significance of FVC and DLCO% predicted in diagnosing IPF and assessing the state of illness in IPF. Based on the prediction model, FVC is not a notable risk factor for AE-IPF following logistic regression analysis, whereas DLCO% predicted was an important predictor of AE-IPF. It is suggested in the 2011 guidelines that the lung function of patients with IPF should be examined every 3–6 months. Pulse oximetry saturation during rest and walking in patients with IPF should also be measured to determine the necessity of oxygen treatment. To evaluate the accuracy of the prediction model, we first calculated the C-index and AUC values. The C-index was 0.80, and the AUC value was 0.77, suggesting that the AE-IPF risk prediction model has an acceptable accuracy. To further examine the stability of the prediction model, we used an internal validation method, by evaluating the clinical data of the 32 patients with IPF in the test group, and then calculated the C-index and AUC values. The C-index and AUC values were 0.75 and 0.70, respectively. The results of the internal validation thus suggest acceptable stability of the AE-IPF risk prediction model. 5. Conclusion A history of occupational exposure, DM, EH, and DLCO% predicted were risk factors for AE-IPF as determined by the logistic regression analysis, and the AE-IPF risk prediction model was successfully constructed by creating a nomogram. The nomogram could clearly and intuitively predict the probability of IPF patients' progress to AE-IPF, thus supplying a reference for clinical progression. However, the current study was a single-center retrospective study with a limited sample size. Therefore, the results of this study should be confirmed in a prospective multicenter, large-scale study. In addition, hierarchical studies of patients of different ages or with different illness severities should also be performed. Acknowledgments This study was sponsored by the Science and Technology Program of Xuzhou Science and Technology Bureau (No. KC17113) and the Outstanding Talent Research Funding of Xuzhou Medical University (Grant No. D2019005 and No. D2019022). Data Availability The data used to support the finding of this study are available from the corresponding author upon request. Conflicts of Interest The authors declare that they have no conflict of interest. Authors' Contributions Qi Wu and Yong Xu conceived and designed this study. Yao Zhou, Ke-jia Zhang, and Shi-min Jiang performed the clinical data collection, analysis, and interpretation and wrote the initial draft. Ying Zhao assisted with data interpretation. Qi Wu and Yong Xu contributed equally to this work and should be considered co-first authors. Figure 1 Nomogram to predict the incidence of AE-IPF. The nomogram was constructed with occupational exposure, DM, EH, and DLCO% predicted. Abbreviations: DM: diabetes mellitus; EH: essential hypertension; DLCO: diffusion capacity for carbon monoxide. Figure 2 Calibration curve for the risk prediction model of AE-IPF. Figure 3 ROC curve corresponding to the risk prediction model of AE-IPF. The AUC was 0.77. Figure 4 ROC curve corresponding to the test group. The AUC was 0.70. Figure 5 DCA for the risk prediction model of AE-IPF. The decision curve showed that, at threshold values of 0.04–0.66, using this nomogram to predict the risk of AE-IPF added more benefit. Table 1 Patient characteristics of the training and test groups. Characteristic Training group (n = 78) Verification group (n = 32) AE-IPF 16 (20.51) 4 (12.50) Demographic and clinical characteristics  Gender   Male 48 (61.54) 18 (56.25)   Female 30 (38.46) 14 (43.75)  Age   ≤60 9 (11.54) 4 (12.50)   61-65 11 (14.10) 2 (6.250)   >65 58 (74.36) 26 (81.25)  Occupational 8 (10.26) 3 (9.38)  Smoking 16 (20.51) 4 (12.50) Comorbidities  DM 17 (21.79) 3 (9.38)  EH 7 (8.97) 6 (18.75)  CHD 10 (12.82) 4 (12.50)  CI 11 (14.10) 5 (15.63)  GERD 11 (14.10) 6 (18.75)  Hypothyroidism 4 (5.13) 1 (3.13)  OSAHS 3 (3.85) 0 Lung function  FVC% predicted   >75 16 (20.51) 4 (12.50)   50-75 47 (60.26) 24 (75.00)   <50 15 (19.23) 4 (12.50)  DLco% predicted   >55 19 (24.36) 6 (18.75)   35-55 41 (52.56) 20 (62.50)   <35 18 (23.08) 6 (18.75) Drug therapy  Acetylcysteine 16 (20.51) 3 (9.38)  Pirfenidone 5 (6.41) 0  Glucocorticoid 12 (15.38) 4 (12.50)  Acid-inhibitory drugs 5 (6.41) 2 (6.25) Abbreviations: AE-IPF: acute exacerbation of idiopathic pulmonary fibrosis; DM: diabetes mellitus; EH: essential hypertension; CHD: coronary heart disease; CI: cerebral infarction; GERD: gastroesophageal reflux disease; OSAHS: obstructive sleep apnea-hypopnea syndrome; FVC: forced vital capacity; DLCO: diffusion capacity for carbon monoxide; NAC: N-acetylcysteine; PFD: pirfenidone; GC: glucocorticoid. Table 2 Final regression model for the primary outcome. Intercept β Odds ratio (95% CI) p value Intercept -5.94 0.00 (0.00-0.23) 0.03 Occupational 5.59 267.90 (6.47-40395.40) 0.01 DM 3.76 42.76 (3.84-1227.58) 0.01 EH -4.50 0.01 (0.00-0.27) 0.02 DLCO% predicted  36-55 3.50 32.97 (1.84-2038.88) 0.04  <35 3.97 52.73 (2.22-3825.05) 0.03 Abbreviations: DM: diabetes mellitus; EH: essential hypertension; DLCO: diffusion capacity for carbon monoxide. ==== Refs 1 Richeldi L. Collard H. R. Jones M. G. Idiopathic pulmonary fibrosis The Lancet 2017 389 10082 1941 1952 10.1016/S0140-6736(17)30866-8 2-s2.0-85016515883 2 Raghu G. Ley B. Brown K. K. Risk factors for disease progression in idiopathic pulmonary fibrosis Thorax 2019 75 1 78 80 31611341 3 Wells A. U. Pamrevlumab in idiopathic pulmonary fibrosis The Lancet Respiratory Medicine 2020 8 1 2 3 10.1016/S2213-2600(19)30339-X 31575507 4 Farrand E. Vittinghoff E. Ley B. Butte A. J. Collard H. R. Corticosteroid use is not associated with improved outcomes in acute exacerbation of IPF Respirology 2020 25 6 629 635 10.1111/resp.13753 31846126 5 Marchioni A. Tonelli R. Ball L. Acute exacerbation of idiopathic pulmonary fibrosis: lessons learned from acute respiratory distress syndrome? Critical Care 2018 22 1 p. 80 10.1186/s13054-018-2002-4 2-s2.0-85044324059 29566734 6 Dotan Y. Vaidy A. Shapiro W. B. Effect of acute exacerbation of idiopathic pulmonary fibrosis on lung transplantation outcome Chest 2018 154 4 818 826 10.1016/j.chest.2018.06.027 2-s2.0-85052754528 29966665 7 Yoshimura K. Nakatani T. Nakamori Y. Acute exacerbation in idiopathic interstitial pneumonia Nihon Kyōbu Shikkan Gakkai zasshi 1984 22 11 1012 1020 6530862 8 Kishaba T. Evaluation and management of idiopathic pulmonary fibrosis Respiratory Investigation 2019 57 4 300 311 10.1016/j.resinv.2019.02.003 2-s2.0-85062425671 30853366 9 Song J. W. Hong S.-B. Lim C.-M. Koh Y. Kim D. S. Acute exacerbation of idiopathic pulmonary fibrosis: incidence, risk factors and outcome The European Respiratory Journal 2011 37 2 356 363 10.1183/09031936.00159709 2-s2.0-79953652811 20595144 10 Kolb M. Richeldi L. Behr J. Nintedanib in patients with idiopathic pulmonary fibrosis and preserved lung volume Thorax 2017 72 4 340 346 10.1136/thoraxjnl-2016-208710 2-s2.0-85017017908 27672117 11 Kondoh Y. Taniguchi H. Ebina M. Risk factors for acute exacerbation of idiopathic pulmonary fibrosis – extended analysis of pirfenidone trial in Japan Respiratory Investigation 2015 53 6 271 278 10.1016/j.resinv.2015.04.005 2-s2.0-84954198075 26521104 12 Jiao Y. Du P. Performance measures in evaluating machine learning based bioinformatics predictors for classifications Quantitative Biology 2016 4 4 320 330 10.1007/s40484-016-0081-2 2-s2.0-85031499899 13 Raghu G. Collard H. R. Egan J. J. An official ATS/ERS/JRS/ALAT statement: idiopathic pulmonary fibrosis: evidence-based guidelines for diagnosis and management American Journal of Respiratory and Critical Care Medicine 2011 183 6 788 824 10.1164/rccm.2009-040GL 2-s2.0-79952717349 21471066 14 Tzouvelekis A. Herazo-Maya J. Sakamoto K. Bouros D. Biomarkers in the evaluation and management of idiopathic pulmonary fibrosis Current Topics in Medicinal Chemistry 2016 16 14 1587 1598 10.2174/1568026616666150930120959 2-s2.0-84964226751 26420365 15 Cho S. J. Moon J.-S. Nikahira K. GLUT1-dependent glycolysis regulates exacerbation of fibrosis via AIM2 inflammasome activation Thorax 2020 75 3 227 236 10.1136/thoraxjnl-2019-213571 31822523 16 Collard H. R. Calfee C. S. Wolters P. J. Plasma biomarker profiles in acute exacerbation of idiopathic pulmonary fibrosis American Journal of Physiology-Lung Cellular and Molecular Physiology 2010 299 1 L3 L7 10.1152/ajplung.90637.2008 2-s2.0-77953491551 20418386 17 Ryerson C. J. Cottin V. Brown K. K. Collard H. R. Acute exacerbation of idiopathic pulmonary fibrosis: shifting the paradigm The European Respiratory Journal 2015 46 2 512 520 10.1183/13993003.00419-2015 2-s2.0-84938837190 26232481 18 Cano-Jiménez E. González F. H. Peloche G. Comorbidities and complications in idiopathic pulmonary fibrosis Medical Sciences 2018 6 3 p. 71 10.3390/medsci6030071 30200249 19 Khateeb J. Fuchs E. Khamaisi M. Diabetes and lung disease: a neglected relationship The Review of Diabetic Studies 2019 15 1 15 10.1900/RDS.2019.15.1 2-s2.0-85066411459 30489598 20 Enomoto T. Usuki J. Azuma A. Nakagawa T. Kudoh S. Diabetes mellitus may increase risk for idiopathic pulmonary Fibrosisa Chest 2003 123 6 2007 2011 10.1378/chest.123.6.2007 2-s2.0-0037790724 12796182 21 Wang D. Ma Y. Tong X. Zhang Y. Fan H. Diabetes mellitus contributes to idiopathic pulmonary fibrosis: a review from clinical appearance to possible pathogenesis Frontiers in Public Health 2020 8 p. 196 10.3389/fpubh.2020.00196 22 Abe M. Tsushima K. Sakayori M. Utility of nintedanib for severe idiopathic pulmonary fibrosis: a single-center retrospective study Drug Design, Development and Therapy 2018 12 3369 3375 10.2147/DDDT.S179427 2-s2.0-85055143727 30349191 23 Collard H. R. Ryerson C. J. Corte T. J. Acute exacerbation of idiopathic pulmonary fibrosis. An international working group report American Journal of Respiratory and Critical Care Medicine 2016 194 3 265 275 10.1164/rccm.201604-0801CI 2-s2.0-84988966334 27299520