Skip to main content

ORIGINAL RESEARCH article

Front. Cell. Infect. Microbiol., 10 March 2022
Sec. Clinical Microbiology

A Clinically Applicable Nomogram for Predicting the Risk of Invasive Mechanical Ventilation in Pneumocystis jirovecii Pneumonia

Rongjun Wan,,,,&#x;Rongjun Wan1,2,3,4,5†Lu Bai,,,,&#x;Lu Bai1,2,3,4,5†Yusheng Yan&#x;Yusheng Yan6†Jianmin Li&#x;Jianmin Li7†Qingkai LuoQingkai Luo8Hua HuangHua Huang9Lingmei HuangLingmei Huang10Zhi XiangZhi Xiang11Qing LuoQing Luo12Zi GuZi Gu13Qing GuoQing Guo14Pinhua Pan,,,,Pinhua Pan1,2,3,4,5Rongli Lu,,,,Rongli Lu1,2,3,4,5Yimin Fang,,,,Yimin Fang1,2,3,4,5Chengping Hu,,,,Chengping Hu1,2,3,4,5Juan Jiang,,,,*Juan Jiang1,2,3,4,5*Yuanyuan Li,,,,*Yuanyuan Li1,2,3,4,5*
  • 1Department of Respiratory Medicine, National Key Clinical Specialty, Branch of National Clinical Research Center for Respiratory Disease, Xiangya Hospital, Central South University, Changsha, China
  • 2Center of Respiratory Medicine, Xiangya Hospital, Central South University, Changsha, China
  • 3Clinical Research Center for Respiratory Diseases in Hunan Province, Changsha, China
  • 4Hunan Engineering Research Center for Intelligent Diagnosis and Treatment of Respiratory Disease, Changsha, China
  • 5National Clinical Research Center for Geriatric Disorders, Xiangya Hospital, Changsha, China
  • 6Department of Pulmonary and Critical Care Medicine, First Hospital of Changsha, Changsha, China
  • 7Department of Pulmonary and Critical Care Medicine, Hunan Provincial People’s Hospital, First Affiliated Hospital of Hunan Normal University, Changsha, China
  • 8Department of Pulmonary and Critical Care Medicine, First People’s Hospital of Chenzhou, Chenzhou, China
  • 9Medical Center of Tuberculosis, Second People’s Hospital of Chenzhou, Chenzhou, China
  • 10Department of Pulmonary and Critical Care Medicine, Yueyang Central Hospital, Yueyang, China
  • 11Department of Respiratory Medicine, First People’s Hospital of Huaihua, Huaihua, China
  • 12Department of Pulmonary and Critical Care Medicine, Second Affiliated Hospital, Hengyang Medical School, University of South China, Hengyang, China
  • 13Department of Pulmonary and Critical Care Medicine, Xiangtan Central Hospital, Xiangtan, China
  • 14Department of Pulmonary and Critical Care Medicine, Yiyang Central Hospital, Yiyang, China

Objective: Pneumocystis jirovecii pneumonia (PCP) is a life-threatening disease associated with a high mortality rate among immunocompromised patient populations. Invasive mechanical ventilation (IMV) is a crucial component of treatment for PCP patients with progressive hypoxemia. This study explored the risk factors for IMV and established a model for early predicting the risk of IMV among patients with PCP.

Methods: A multicenter, observational cohort study was conducted in 10 hospitals in China. Patients diagnosed with PCP were included, and their baseline clinical characteristics were collected. A Boruta analysis was performed to identify potentially important clinical features associated with the use of IMV during hospitalization. Selected variables were further analyzed using univariate and multivariable logistic regression. A logistic regression model was established based on independent risk factors for IMV and visualized using a nomogram.

Results: In total, 103 patients comprised the training cohort for model development, and 45 comprised the validation cohort to confirm the model’s performance. No significant differences were observed in baseline clinical characteristics between the training and validation cohorts. Boruta analysis identified eight clinical features associated with IMV, three of which were further confirmed to be independent risk factors for IMV, including age (odds ratio [OR] 2.615 [95% confidence interval (CI) 1.110–6.159]; p = 0.028), oxygenation index (OR 0.217 [95% CI 0.078–0.604]; p = 0.003), and serum lactate dehydrogenase level (OR 1.864 [95% CI 1.040–3.341]; p = 0.037). Incorporating these three variables, the nomogram achieved good concordance indices of 0.829 (95% CI 0.752–0.906) and 0.818 (95% CI 0.686–0.950) in predicting IMV in the training and validation cohorts, respectively, and had well-fitted calibration curves.

Conclusions: The nomogram demonstrated accurate prediction of IMV in patients with PCP. Clinical application of this model enables early identification of patients with PCP who require IMV, which, in turn, may lead to rational therapeutic choices and improved clinical outcomes.

Introduction

Pneumocystis jirovecii pneumonia (PCP) is an opportunistic fungal infection that commonly affects immunocompromised patient populations. It remains a severe disease associated with a high mortality rate among patients either positive or negative for human immunodeficiency virus (HIV) infection (Roux et al., 2014; Alanio et al., 2016). While the incidence of PCP among HIV-positive patients has declined with the widespread use of antiretroviral and prophylactic therapy (Morris et al., 2004; Salzer et al., 2018), PCP has become more frequent in HIV-negative populations in recent years due to the increasing prevalence of immunosuppressive conditions such as hematological malignancies, solid tumors, systemic corticosteroid and immunosuppressive therapies, organ transplantation, and hematopoietic stem cell transplantation (Jin et al., 2021; Wang et al., 2021). PCP has been reported to be associated with mortality rates as high as 10%–20% in HIV and 30%–60% in non-HIV patients (Cilloniz et al., 2019).

Patients with PCP are usually characterized by acute or subacute hypoxemia and exhibit symptoms of progressive dyspnea. Adherence of Pneumocystis to the alveoli triggers a severe inflammatory response in hosts that causes diffuse lung injury and impairs gas exchange (Hahn and Limper, 2003). In a substantial proportion of cases, hypoxemia can progress to respiratory failure and even acute respiratory distress syndrome within a short period of time. Mechanical ventilation, including invasive mechanical ventilation (IMV), is a crucial component of the clinical management of patients with PCP. According to data from previous studies, 16% of HIV-infected and 50%–64% of non-HIV-infected PCP patients require mechanical ventilation during hospitalization (Nüesch et al., 1999; Kim et al., 2014). A large cohort study also reported that >40% of patients required admission to the intensive care unit, and approximately 37% required mechanical ventilation (Schmidt et al., 2018). Approximately 28% of patients with PCP require IMV. In particular, the use of IMV is more frequent (56.6%) in patients with severe PCP (Gaborit et al., 2019). For patients who require mechanical ventilation, mortality can be as high as 62% (Monnet et al., 2008). Late intubation has been recognized to be associated with increased mortality in those with acute respiratory distress syndrome (Kangelaris et al., 2016). Furthermore, it has been found that failure of non-invasive mechanical ventilation is significantly associated with hospital survival among PCP patients, with a mortality rate as high as 88.9% (Wang et al., 2020). Therefore, early identification of patients who require IMV can provide guidance for prompt initiation of IMV and optimize the allocation of medical resources. However, while prognostic factors for mortality have been well established (Walzer et al., 2008; Fei et al., 2009; Wang et al., 2011; Weng et al., 2016; Gaborit et al., 2019), little is known about the early prediction of requirements for IMV in patients with PCP so far.

Herein, we developed a visual and easy-to-use nomogram to predict the risk for IMV in PCP patients on admission to the hospital using a machine learning method.

Methods

Study Design and Patients

The present multicenter, observational cohort study collected clinical data from all immunosuppressed patients who were diagnosed with PCP between January 1, 2018, and October 31, 2021, in 10 tertiary care hospitals in China. Clinical data included demographic information, clinical manifestations, comorbidities, laboratory indicators, vital signs, anti-Pneumocystis jirovecii treatment and respiratory support approaches. All baseline laboratory tests were performed within 24 h of patient admission. Adult patients (age ≥ 18 years) were eligible for enrollment if they fulfilled all of the following criteria: immunosuppressed conditions, including HIV-positive and HIV-negative patients (Ramirez et al., 2020); typical clinical manifestations of PCP, including fever, cough (usually dry cough), and progressive dyspnea; radiological findings suggestive of PCP in the bilateral lungs that newly emerged on chest computed tomography; and a positive result for Pneumocystis jirovecii by Gomori methenamine silver staining, real-time polymerase chain reaction assay, and/or metagenomic next-generation sequencing (Jiang et al., 2021). Diagnosis of PCP was confirmed after consultation with ≥ 2 senior pulmonary specialists at each research center, based on clinical symptoms, laboratory findings, chest radiology, microbiological tests, and treatment response. Patients < 18 years of age, those who were pregnant, and those with incomplete medical records were excluded from the study. Eligible patients, who were diagnosed with PCP between January 1, 2018, and November 30, 2020, were included in the training cohort. Clinical data from these patients were retrospectively collected for the development of a nomogram. Patients who were diagnosed with PCP between January 1, 2021, and October 31, 2021, were included in the validation cohort. Clinical data from these patients were prospectively collected to validate the nomogram.

Relevant Variable Identification Using Random Forest

To identify clinical features that had a larger than random association with IMV use, Boruta analysis was performed (Kursa and Rudnicki, 2010) to cover all data regarding selected variables obtained when patients were admitted to hospital. The Boruta analysis is a random forest-based feature selection algorithm (Breiman, 2001) that performs multiple runs of random forest to compare shuffled random variables to the original variables. At the same time, scores representing importance are assigned to each variable. All analyzed variables were divided into rejected, tentative, and accepted groups based on their importance scores. Briefly, variables with importance scores below the shadowMax value are rejected by Boruta, while those above are accepted; other variables with importance scores around shadowMax value, which cannot be rejected or accepted, are divided into a tentative group. All Boruta analyses were performed using the Boruta package (Kursa and Rudnicki, 2010) version 6.0.0 using the default parameters (ntree = 500, maxRuns = 100, p = 0.01). After the Boruta analysis, the accepted variables were further analyzed using univariable and multivariable logistic regression to calculate odds ratio (OR) and corresponding 95% confidence interval (CI) to confirm whether the selected variables are risk factors for IMV in patients with PCP.

Development and Validation of a Logistic Regression Model

To develop a practical and easy-to-use predictive model for IMV using clinical data acquired on admission, clinical data from all enrolled patients in the training cohort (n = 103) were analyzed to start the establishment procedure (Table 1). Important variables identified by the Boruta algorithm were further analyzed using univariable logistic regression to calculate the OR and to assess significance in the training cohort. Variables that were statistically significant were further included in the multivariable logistic regression analysis to confirm whether they were independent risk factors. Subsequently, a logistic regression model was established based on factors with statistical significance in multivariable logistic regression and visualized using a nomogram. Briefly, the nomogram is based on proportionally converting each regression coefficient in multivariate logistic regression analysis to a 0- to 100-point scale. The effect of the variable with the highest β coefficient (absolute value) was assigned 100 points. The points are added across independent variables to obtain the total points, which are converted to predicted probabilities (Lei et al., 2016). Furthermore, the model was evaluated in a validation cohort (n = 45), for which the number of patients was determined based on the ratio of the samples in the training and validation sets recommended in previous studies (Frizzell et al., 2017; Huang et al., 2018). Predictive performance was assessed by discrimination and calibration.

TABLE 1
www.frontiersin.org

Table 1 Baseline characteristics of patients.

Statistical Analysis

All continuous variables were tested for normal distribution using the Lilliefors test (Lilliefors, 1967) before statistical expression and comparison. Normally distributed continuous variables are expressed as mean ± standard deviation (SD) and compared using Student’s t-test, while non-normally distributed continuous variables are expressed as median (interquartile range) and were compared using the Mann–Whitney U test. Categorical variables are expressed as count and percentage and were compared using the chi-squared or Fisher’s exact tests.

Important variables selected by the Boruta algorithm were successively included in the univariate and multivariate logistic regression analyses, and a clinical predictive model was constructed using multivariable logistic regression. Furthermore, the discriminative ability of the model was assessed using receiver operator characteristic (ROC) curve and area under the ROC curve (AUROC), which was assessed using a bootstrap-based calibration process (Steyerberg and Vergouwe, 2014; Lindhiem et al., 2020) and Spiegelhalter’s z test (Spiegelhalter, 1986) to evaluate the relationship between the actual and predicted probability of IMV. A non-significant result in the Spiegelhalter’s z-test indicated that the model was well fitted. These analyses were performed using the rms version 5.1 and pROC version 1.16.2 R packages.

All statistical analyses were performed using R software version 3.6.3 (www.r-project.org). Except for the nomogram, figures were plotted using ggplot2 version 3.3.2 R package. For all analyses, differences with p < 0.05 were considered to be statistically significant.

Ethics Statement

The present study was approved by the Ethics Committee of Xiangya Hospital, Central South University (Changsha, Hunan, China; 202112655). Due to the use of de-identified and anonymized research data, requirements for informed written consent were waived. The study was conducted in accordance with the Declaration of Helsinki.

Results

Baseline Characteristics of Patients

A total of 148 PCP patients who met the inclusion criteria were enrolled in this study. Of these, 103 patients were included into the training cohort, and 45 were in validation cohort. Baseline clinical characteristics are shown in Table 1. The average age of all patients was 51 years old, and 66.2% of them were men. The most common comorbidities included AIDS (30.4%), rheumatic diseases (25.7%), and hematological malignancies (10.8%). Marked laboratory abnormalities were observed, including decreased oxygenation index (the ratio of PaO2 to FiO2), lymphocytopenia, elevated serum levels of (1,3)-β-D-glucan, and lactate dehydrogenase (LDH). Co-infections were common in patients with PCP. Bacterial and viral co-infections were found in 56.8% and 30.1% of patients, respectively. Among all patients, 48.0% received TMP-SMZ only, while 45.3% received TMP-SMZ combined with echinocandins. The median time from admission to initiation of anti-Pneumocystis jirovecii treatment was 1 day. Systemic glucocorticoids was used in 79.1% of patients. Overall, there was no significant difference in baseline clinical characteristics between training and validation cohorts. IMV was used in 45 (43.7%) and 23 (51.1%) patients in these two cohorts, respectively.

Identification of Important Variables by Boruta Analysis

Based on the baseline characteristics of all patients, 38 clinical features were evaluated by using Boruta feature selection method. The rank plot of feature importance is shown in Figure 1. There were 8 clinical features confirmed as important, including oxygenation index, AIDS, combination with echinocandins, age, serum LDH, procalcitonin, serum creatinine, and blood urea nitrogen. Unsurprisingly, baseline oxygenation index had the highest importance associated with the risk of IMV. The other variables with an importance score lower than shadowMax were all confirmed as unimportant.

FIGURE 1
www.frontiersin.org

Figure 1 The feature importance in the Boruta feature selection process. The red box showed the features that are confirmed as important, the green box showed the features confirmed as unimportant, the yellow box showed the features that are tentative, and the blue box showed Boruta parameters. AIDS, acquired immune deficiency syndrome; PTB, pulmonary tuberculosis; COPD, chronic obstructive pulmonary disease; ILD, interstitial lung disease; OI, oxygenation index; WBC, white blood cell; Neu, neutrophil; Lym, lymphocyte; Hb, hemoglobin; PLT, platelet; PCT, procalcitonin; CRP, C-reaction protein; GM, galactomannan; LDH, lactate dehydrogenase; TBil, total bilirubin; SCr, serum creatine; BUN, blood urea nitrogen; HR, heart rate; RR, respiratory rate; SpO2, the percent saturation of oxygen in the blood; SBP, systolic blood pressure; DBP, diastolic blood pressure.

Development and Validation of a Predictive Nomogram of IMV

Next, these 8 important variables selected by Boruta were further analyzed by univariate logistic regression, and the results are presented in Table 2. There were 5 variables associated with the risk of IMV, including age, AIDS, oxygenation index, serum LDH, and blood urea nitrogen. In subsequent multivariable logistic regression, age (OR, 2.615; 95% CI, 1.110–6.159; p = 0.028), oxygenation index (OR, 0.217; 95% CI, 0.078–0.604; p = 0.003), and serum LDH (OR, 1.864; 95% CI, 1.040-3.341; p = 0.037) were confirmed as independently associated with the risk of IMV in PCP patients (Table 3). These factors were used to construct a predictive model for IMV risk by multivariable logistic regression, which was further visualized by nomogram as shown in Figure 2.

TABLE 2
www.frontiersin.org

Table 2 Univariate logistic regression analysis of invasive mechanical ventilation based on baseline characteristics in the training cohort.

TABLE 3
www.frontiersin.org

Table 3 Multivariate logistic regression analysis of invasive mechanical ventilation based on baseline characteristics in the training cohort.

FIGURE 2
www.frontiersin.org

Figure 2 Nomogram for the estimation of invasive mechanical ventilation risk in PCP patients. To use the nomogram, find the position of each variable on the corresponding axis, and draw a line to the points axis (top) for the number of points. Then, sum up the points from all of the variables, and draw a line from the total points axis to the lower line of the nomogram to determine the predicted risk of invasive mechanical ventilation. OI, oxygenation index; LDH, lactate dehydrogenase.

The predictive model was internally validated using the bootstrap validation method. Receiver operating characteristic curves for the predictive model of training and validation cohorts are shown in Figures 3A, B. The nomogram demonstrated good accuracy in estimating the risk of IMV, with a bootstrap-corrected C index of 0.829 (95% CI, 0.752–0.906). Furthermore, calibration plots graphically showed good agreement between the estimated risk by the nomogram and the actual use of IMV in PCP patients (Figure 3C). Next, we applied our model to the validation cohort and performed relative analysis to evaluate the performance of nomogram. the results displayed a C index of 0.818 (95% CI, 0.686–0.950) for the estimation of IMV risk. There was also a good calibration curve for the risk estimation (Figure 3D). The sensitivity, specificity, positive predictive value, and negative predictive value of this nomogram in identifying the patients who require IMV were 88.9%, 65.5%, 66.7%, and 88.4% in the training cohort, and 82.6%, 72.7%, 76.0%, and 80.0% in the validation cohort, respectively.

FIGURE 3
www.frontiersin.org

Figure 3 Receiver operating characteristic curves and predictive performance of the nomogram for invasive mechanical ventilation in PCP patients. (A) ROC curve for the predictive model of the training cohort. Area under the curve was 0.829. (B) ROC curve for the predictive model of the validation cohort. Area under the curve was 0.818. (C) Validity of the predictive performance of the nomogram in estimating the risk of invasive mechanical ventilation in the training cohort. (D) Validity of the predictive performance of the nomogram in estimating the risk of invasive mechanical ventilation in the validation cohort. C index, concordance index; ROC, receiver operating characteristic.

Discussion

PCP is a life-threatening disease associated with the frequent use of IMV in immunocompromised patient populations. In this multicenter study, we designed and developed an easy-to-use visual nomogram to predict the risk for IMV in patients with PCP. The nomogram demonstrated good predictive ability for IMV using three common patient variables collected on admission; more specifically, age, oxygenation index, and serum LDH level. To the best of our knowledge, this is the first study to report a clinically applicable predictive model for the early estimation of IMV risk in patients with PCP.

In this study, our data demonstrated that oxygenation index was the most important risk factor for IMV. It is well recognized that the severity of PCP is determined based on the partial pressure of oxygen (PaO2) while breathing room air and at rest. Patients were classified as having mild (PaO2 > 70 mmHg), moderate (60 mmHg ≤ PaO2 ≤ 70 mmHg), or severe (PaO2 < 60 mmHg) PCP (Moon et al., 2011). A low oxygenation index reflects severe hypoxemia, which conceivably increases the risk for disease deterioration and requirement for IMV. Previous studies have demonstrated that poor oxygenation on admission to the hospital is a risk factor for mortality among HIV-infected patients (Walzer et al., 2008; Fei et al., 2009). Moreover, it has also been reported that age is independently associated with disease severity and mortality in patients with PCP. In a prospective observational study that enrolled 107 patients with PCP, Gaborit et al. reported that older age (OR 3.36 [95% CI 1.4–8.5]) was independently associated with disease severity (Gaborit et al., 2019). Another bicentric retrospective study also suggested that age (OR 1.051 [95% CI 1.007–1.097]) was an independent risk factor for hospital mortality for non-HIV patients in the intensive care unit, although its effect was only slight. Old age has been consistently correlated with worse prognosis in patients with connective tissue disease-associated PCP (Ishikawa et al., 2021). In the present study, multivariate logistic regression analysis confirmed that age was independently associated with the risk for IMV (OR 2.615 [95% CI 1.110–6.159]), indicating that older PCP patients are more likely to require IMV during hospitalization. Together with previous studies, our data support that close attention should be devoted to elderly patients with PCP.

LDH has long been recognized as a useful but non-specific serological marker in the diagnosis of PCP (White et al., 2017). Serum LDH levels can reflect the severity of radiographic abnormalities in PCP (Boldt and Bai, 1997). Schmidt et al. found that baseline serum LDH level was a predictor of in-hospital mortality (OR 1.17 [95% CI 1.09–1.27]) (Schmidt et al., 2018). However, there is no established association between serum LDH level and the risk for IMV in PCP. Interestingly, our data revealed that serum LDH level was an independent risk factor for IMV during hospitalization (OR 1.864 [95% CI 1.040–3.341]). This can be explained by the hypothesis that extracellular LDH levels indicate cell damage or cell death, and elevated serum LDH level is correlated with lung tissue damage. In fact, serum LDH level has been reported to be significantly higher (9.79 ± 7.0 versus 5.93 ± 3.9 μkat/L) among patients with severe PCP compared to those with non-severe PCP (Gaborit et al., 2019), which supports the hypothesis proposed above.

Based on our data, it is noteworthy that HIV infection was not independently associated with the risk for IMV among patients with PCP. It is well known that HIV-negative PCP patients usually experience a more progressive disease course, more severe hypoxemia, and a higher mortality rate (Roux et al., 2014; Cilloniz et al., 2019; Ghembaza et al., 2020). In the present study, HIV infection appeared to be associated with IMV use according to univariate logistic regression analysis. However, no statistical significance was observed in the multivariate logistic regression analysis. The effect of confounding factors can be an important reason, including oxygenation index, age, serum level of LDH, and blood urea nitrogen. It has been known that clinical characteristics of HIV-positive and HIV-negative PCP patients are markedly different. Overall, HIV-positive PCP patients tend to be younger and have higher oxygenation index, lower serum LDH, and lower blood urea nitrogen levels (Mundo et al., 2020; Feng et al., 2021). Therefore, these variables, which reflect the severity of disease, can be confounding factors for AIDS/HIV infection in multivariate logistic regression analysis. Our results indicate that HIV infection itself did not directly affect the risk for IMV, even though disease severity is more prominent in HIV-negative PCP patients.

The use of a nomogram for estimating the risk for IMV in patients with PCP is a new concept. Our predictive model is visual and easy to use, enabling early and accurate identification of patients with PCP who require IMV before clinical decompensation. By calculating the total points based on three easily accessible clinical variables (i.e., age, oxygenation index, and serum LDH level), physicians can predict the risk for IMV for HIV-positive and HIV-negative PCP patients on admission to hospital in either urban or rural settings. More importantly, this model demonstrated good predictive performance, with C indexes of 0.829 (95% CI 0.752–0.906) and 0.818 (95% CI 0.686–0.950) in predicting IMV in the training and validation cohorts. In addition, the calibration curves were well fitted. The use of this predictive model could provide important guidance for clinical decision-making and optimizing the allocation of medical resources, with the hope of improving patient outcomes and quality of care. Meanwhile, it is of great significance to apply this predictive model in primary hospitals, where IMV is usually not possible. Physicians can easily estimate the risk of IMV once patients are admitted, which is useful in deciding whether to transfer them to higher-level hospitals. By using this predictive model, early and accurate identification of patients who are likely to require IMV will enable closer monitoring for signs and symptoms of clinical deterioration and optimize allocation of medical resources, which can facilitate clinical decision-making and may potentially improve the clinical outcomes of patients with PCP.

There were several limitations to the present study. First, clinical data from patients in the training cohort were retrospectively collected; thus, inherent bias was unavoidable. Second, although 10 tertiary care hospitals participated in this study, the sample size was relatively small, which may have affected the validity of our predictive model. Finally, whether this predictive model will improve the clinical outcomes of patients with PCP by enabling early identification of patients who require IMV warrants further evaluation in large-scale prospective studies. Nevertheless, we believe that this visual predictive model represents a useful tool for clinical physicians who are focused on PCP—especially severe PCP—but needs confirmation in future clinical trials.

Conclusion

We developed a clinically applicable nomogram for predicting the risk for IMV in patients with PCP. This nomogram enables early and accurate estimation of IMV risk when PCP patients are admitted to the hospital, which may provide guidance for clinical decision-making.

Data Availability Statement

Data are available related to this study can be accessible from the corresponding authors with reasonable requests.

Ethics Statement

The present study was approved by the Ethics Committee of Xiangya Hospital, Central South University (Changsha, Hunan, China; 202112655). Due to the use of deidentified and anonymized research data, requirements for informed written consent were waived. The study was conducted in accordance with the Declaration of Helsinki.

Author Contributions

Study concept and design: YL and JJ. Acquisition of data: YY, JL, LB, RW, QKL, HH, LH, ZX, QL, ZG, QG, PP, RL, YF, and CH. Statistical analysis: RW and LB. Analysis and interpretation of data: YL, JJ, RW, LB, YY, and JL. Drafting of the manuscript: JJ, LB, and RW. Critical revision of the manuscript: YL and JJ. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by grants from the National Natural Science Foundation of China (82170041, 82100099, and 81873406), the Innovative Research Platform of Hunan Development and Reform Commission (2021-212), and the Youth Research Foundation of Xiangya Hospital (2018Q015).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We appreciate the professionalism and compassion demonstrated by all the healthcare workers involved in patient care. We also acknowledge all the patients for their involvement in this study.

References

Alanio, A., Hauser, P. M., Lagrou, K., Melchers, W. J., Helweg-Larsen, J., Matos, O., et al. (2016). ECIL Guidelines for the Diagnosis of Pneumocystis Jirovecii Pneumonia in Patients With Haematological Malignancies and Stem Cell Transplant Recipients. J. Antimicrob. Chemother. 71 (9), 2386–2396. doi: 10.1093/jac/dkw156

PubMed Abstract | CrossRef Full Text | Google Scholar

Boldt, M. J., Bai, T. R. (1997). Utility of Lactate Dehydrogenase vs Radiographic Severity in the Differential Diagnosis of Pneumocystis Carinii Pneumonia. Chest 111 (5), 1187–1192. doi: 10.1378/chest.111.5.1187

PubMed Abstract | CrossRef Full Text | Google Scholar

Breiman, L. (2001). Random Forests. Mach. Learn. 45 (1), 5–32. doi: 10.1023/A:1010933404324

CrossRef Full Text | Google Scholar

Cilloniz, C., Dominedo, C., Alvarez-Martinez, M. J., Moreno, A., Garcia, F., Torres, A., et al. (2019). Pneumocystis Pneumonia in the Twenty-First Century: HIV-Infected Versus HIV-Uninfected Patients. Expert Rev. Anti Infect. Ther. 17 (10), 787–801. doi: 10.1080/14787210.2019.1671823

PubMed Abstract | CrossRef Full Text | Google Scholar

Fei, M. W., Kim, E. J., Sant, C. A., Jarlsberg, L. G., Davis, J. L., Swartzman, A., et al. (2009). Predicting Mortality From HIV-Associated Pneumocystis Pneumonia at Illness Presentation: An Observational Cohort Study. Thorax 64 (12), 1070–1076. doi: 10.1136/thx.2009.117846

PubMed Abstract | CrossRef Full Text | Google Scholar

Feng, Q., Hao, J., Li, A., Tong, Z. (2021). Prognostic Nomogram for Death From Pneumocystis Pneumonia in Non-HIV-And HIV-Infected Patients. doi: 10.21203/rs.3.rs-487120/v2

CrossRef Full Text | Google Scholar

Frizzell, J. D., Liang, L., Schulte, P. J., Yancy, C. W., Heidenreich, P. A., Hernandez, A. F., et al. (2017). Prediction of 30-Day All-Cause Readmissions in Patients Hospitalized for Heart Failure: Comparison of Machine Learning and Other Statistical Approaches. JAMA Cardiol. 2 (2), 204–209. doi: 10.1001/jamacardio.2016.3956

PubMed Abstract | CrossRef Full Text | Google Scholar

Gaborit, B. J., Tessoulin, B., Lavergne, R.-A., Morio, F., Sagan, C., Canet, E., et al. (2019). Outcome and Prognostic Factors of Pneumocystis Jirovecii Pneumonia in Immunocompromised Adults: A Prospective Observational Study. Ann. Intensive Care 9 (1), 1–10. doi: 10.1186/s13613-019-0604-x

CrossRef Full Text | Google Scholar

Ghembaza, A., Vautier, M., Cacoub, P., Pourcher, V., Saadoun, D. (2020). Risk Factors and Prevention of Pneumocystis Jirovecii Pneumonia in Patients With Autoimmune and Inflammatory Diseases. Chest 158 (6), 2323–2332. doi: 10.1016/j.chest.2020.05.558

PubMed Abstract | CrossRef Full Text | Google Scholar

Hahn, P. Y., Limper, A. H. (2003). The Role of Inflammation in Respiratory Impairment During Pneumocystis carinii Pneumonia. Semin. Respir. Infect. 18, 40–47. doi: 10.1053/srin.2003.50004

CrossRef Full Text | Google Scholar

Huang, C., Murugiah, K., Mahajan, S., Li, S.-X., Dhruva, S. S., Haimovich, J. S., et al. (2018). Enhancing the Prediction of Acute Kidney Injury Risk After Percutaneous Coronary Intervention Using Machine Learning Techniques: A Retrospective Cohort Study. PloS Med. 15 (11), e1002703. doi: 10.1371/journal.pmed.1002703

PubMed Abstract | CrossRef Full Text | Google Scholar

Ishikawa, Y., Nakano, K., Tokutsu, K., Miyata, H., Fujino, Y., Matsuda, S., et al. (2021). Estimation of Treatment and Prognostic Factors of Pneumocystis Pneumonia in Patients With Connective Tissue Diseases. RMD Open 7 (1), e001508. doi: 10.1136/rmdopen-2020-001508

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, J., Bai, L., Yang, W., Peng, W., An, J., Wu, Y., et al. (2021). Metagenomic Next-Generation Sequencing for the Diagnosis of Pneumocystis Jirovecii Pneumonia in Non-HIV-Infected Patients: A Retrospective Study. Infect. Dis. Ther. 10 (3), 1733–1745. doi: 10.1007/s40121-021-00482-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, F., Xie, J., Wang, H.-L. (2021). Lymphocyte Subset Analysis to Evaluate the Prognosis of HIV-Negative Patients With Pneumocystis Pneumonia. BMC Infect. Dis. 21 (1), 1–9. doi: 10.1186/s12879-021-06124-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Kangelaris, K. N., Ware, L. B., Wang, C. Y., Janz, D. R., Hanjing, Z., Matthay, M. A., et al. (2016). Timing of Intubation and Clinical Outcomes in Adults With ARDS. Crit. Care Med. 44 (1), 120. doi: 10.1097/CCM.0000000000001359

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, S. J., Lee, J., Cho, Y.-J., Park, Y. S., Lee, C.-H., Yoon, H. I., et al. (2014). Prognostic Factors of Pneumocystis Jirovecii Pneumonia in Patients Without HIV Infection. J. Infect. 69 (1), 88–95. doi: 10.1016/j.jinf.2014.02.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Kursa, M. B., Rudnicki, W. R. (2010). Feature Selection With the Boruta Package. J. Stat. Softw. 36 (11), 1 – 13. doi: 10.18637/jss.v036.i11

CrossRef Full Text | Google Scholar

Lei, Z., Li, J., Wu, D., Xia, Y., Wang, Q., Si, A., et al. (2016). Nomogram for Preoperative Estimation of Microvascular Invasion Risk in Hepatitis B Virus–Related Hepatocellular Carcinoma Within the Milan Criteria. JAMA Surg. 151 (4), 356–363. doi: 10.1001/jamasurg.2015.4257

PubMed Abstract | CrossRef Full Text | Google Scholar

Lilliefors, H. W. (1967). On the Kolmogorov-Smirnov Test for Normality With Mean and Variance Unknown. J. Am. Stat. Assoc. 62, 399–402. doi: 10.1080/01621459.1967.10482916

CrossRef Full Text | Google Scholar

Lindhiem, O., Petersen, I. T., Mentch, L. K., Youngstrom, E. A. (2020). The Importance of Calibration in Clinical Psychology. Assessment 27 (4), 840–854. doi: 10.1177/1073191117752055

PubMed Abstract | CrossRef Full Text | Google Scholar

Monnet, X., Vidal-Petiot, E., Osman, D., Hamzaoui, O., Durrbach, A., Goujard, C., et al. (2008). Critical Care Management and Outcome of Severe Pneumocystis Pneumonia in Patients With and Without HIV Infection. Crit. Care 12 (1), 1–9. doi: 10.1186/cc6806

CrossRef Full Text | Google Scholar

Moon, S. M., Kim, T., Sung, H., Kim, M. N., Kim, S. H., Choi, S. H., et al. (2011). Outcomes of Moderate-to-Severe Pneumocystis Pneumonia Treated With Adjunctive Steroid in Non-HIV-Infected Patients. Antimicrob. Agents Chemother. 55 (10), 4613–4618. doi: 10.1128/aac.00669-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Morris, A., Lundgren, J. D., Masur, H., Walzer, P. D., Hanson, D. L., Frederick, T., et al. (2004). Current Epidemiology of Pneumocystis Pneumonia. Emerg. Infect. Dis. 10 (10), 1713–1720. doi: 10.3201/eid1010.030985

PubMed Abstract | CrossRef Full Text | Google Scholar

Mundo, W., Morales-Shnaider, L., Tewahade, S., Wagner, E., Archuleta, S., Bandali, M., et al. (2020). Lower Mortality Associated With Adjuvant Corticosteroid Therapy in Non-HIV-Infected Patients With Pneumocystis jirovecii Pneumonia: A Single-Institution Retrospective US Cohort Study. Open Forum Infect. Dis. 7 (9), 1–7. doi: 10.1093/ofid/ofaa354

CrossRef Full Text | Google Scholar

Nüesch, R., Bellini, C., Zimmerli, W. (1999). Pneumocystis Carinii Pneumonia in Human Immunodeficiency Virus (HIV)—Positive and HIV-Negative Immunocompromised Patients. Clin. Infect. Dis. 29 (6), 1519–1523. doi: 10.1086/313534

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramirez, J. A., Musher, D. M., Evans, S. E., Cruz, C. D., Crothers, K. A., Hage, C. A., et al. (2020). Treatment of Community-Acquired Pneumonia in Immunocompromised Adults: A Consensus Statement Regarding Initial Strategies. Chest 158 (5), 1896–1911. doi: 10.1016/j.chest.2020.05.598

PubMed Abstract | CrossRef Full Text | Google Scholar

Roux, A., Canet, E., Valade, S., Gangneux-Robert, F., Hamane, S., Lafabrie, A., et al. (2014). Pneumocystis Jirovecii Pneumonia in Patients With or Without AIDS, France. Emerg. Infect. Dis. 20 (9), 1490–1497. doi: 10.3201/eid2009.131668

PubMed Abstract | CrossRef Full Text | Google Scholar

Salzer, H. J., Schäfer, G., Hoenigl, M., Günther, G., Hoffmann, C., Kalsdorf, B., et al. (2018). Clinical, Diagnostic, and Treatment Disparities Between HIV-Infected and Non-HIV-Infected Immunocompromised Patients With Pneumocystis Jirovecii Pneumonia. Respiration 96 (1), 52–65. doi: 10.1159/000487713

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmidt, J. J., Lueck, C., Ziesing, S., Stoll, M., Haller, H., Gottlieb, J., et al. (2018). Clinical Course, Treatment and Outcome of Pneumocystis Pneumonia in Immunocompromised Adults: A Retrospective Analysis Over 17 Years. Crit. Care 22 (1), 1–9. doi: 10.1186/s13054-018-2221-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Spiegelhalter, D. J. (1986). Probabilistic Prediction in Patient Management and Clinical Trials. Stat. Med. 5 (5), 421–433. doi: 10.1002/sim.4780050506

PubMed Abstract | CrossRef Full Text | Google Scholar

Steyerberg, E. W., Vergouwe, Y. (2014). Towards Better Clinical Prediction Models: Seven Steps for Development and an ABCD for Validation. Eur. Heart J. 35 (29), 1925–1931. doi: 10.1093/eurheartj/ehu207

PubMed Abstract | CrossRef Full Text | Google Scholar

Walzer, P. D., Evans, H. E., Copas, A. J., Edwards, S. G., Grant, A. D., Miller, R. F. (2008). Early Predictors of Mortality From Pneumocystis Jirovecii Pneumonia in HIV-Infected Patients: 1985-2006. Clin. Infect. Dis. 46 (4), 625–633. doi: 10.1086/526778

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, H.-W., Lin, C.-C., Kuo, C.-F., Liu, C.-P., Lee, C.-M. (2011). Mortality Predictors of Pneumocystis Jirovecii Pneumonia in Human Immunodeficiency Virus-Infected Patients at Presentation: Experience in a Tertiary Care Hospital of Northern Taiwan. J. Microbiol. Immunol. Infect. 44 (4), 274–281. doi: 10.1016/j.jmii.2010.08.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, M., Wu, Y., Qu, S., Yang, X., Song, L. (2020). Risk Factors for Noninvasive Ventilation Failure in Acute Respiratory Failure With Pneumocystis Jirovecii Pneumonia in Patients Without Human Immunodeficiency Virus Infection. Eur. Respir. Soc. 56 (Suppl. 64). doi: 10.1183/13993003.congress-2020.2020

CrossRef Full Text | Google Scholar

Wang, Y., Zhou, X., Saimi, M., Huang, X., Sun, T., Fan, G., et al. (2021). Risk Factors of Mortality From Pneumocystis Pneumonia in Non-HIV Patients: A Meta-Analysis. Front. Public Health 9, 704. doi: 10.3389/fpubh.2021.680108

CrossRef Full Text | Google Scholar

Weng, L., Huang, X., Chen, L., Feng, L. Q., Jiang, W., Hu, X. Y., et al. (2016). Prognostic Factors for Severe Pneumocystis Jiroveci Pneumonia of Non-HIV Patients in Intensive Care Unit: A Bicentric Retrospective Study. BMC Infect. Dis. 16 (1), 528. doi: 10.1186/s12879-016-1855-x

PubMed Abstract | CrossRef Full Text | Google Scholar

White, P. L., Backx, M., Barnes, R. A. (2017). Diagnosis and Management of Pneumocystis Jirovecii Infection. Expert Rev. Anti Infect. Ther. 15 (5), 435–447. doi: 10.1080/14787210.2017.1305887

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Pneumocystis jirovecii pneumonia (PCP), invasive mechanical ventilation (IMV), predictive model, nomogram, machine learning

Citation: Wan R, Bai L, Yan Y, Li J, Luo Q, Huang H, Huang L, Xiang Z, Luo Q, Gu Z, Guo Q, Pan P, Lu R, Fang Y, Hu C, Jiang J and Li Y (2022) A Clinically Applicable Nomogram for Predicting the Risk of Invasive Mechanical Ventilation in Pneumocystis jirovecii Pneumonia. Front. Cell. Infect. Microbiol. 12:850741. doi: 10.3389/fcimb.2022.850741

Received: 08 January 2022; Accepted: 07 February 2022;
Published: 10 March 2022.

Edited by:

Francisco Antunes, Universidade de Lisboa, Portugal

Reviewed by:

Francisco Esteves, New University of Lisbon, Portugal
Enrique Calderón, Spanish National Research Council (CSIC), Spain

Copyright © 2022 Wan, Bai, Yan, Li, Luo, Huang, Huang, Xiang, Luo, Gu, Guo, Pan, Lu, Fang, Hu, Jiang and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yuanyuan Li, bGVlcm91bmRAY3N1LmVkdS5jbg==; Juan Jiang, anVhbmppYW5nQGNzdS5lZHUuY24=

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.