Construction and Evaluation of a Model for the Prediction of Risk Factors Associated with Severe Pneumonia in Patients with Positive Hepatitis B Core Antibody and COVID-19 Infection in China

authors:

avatar Sun Zhenmin 1 , avatar Yang Nan 1 , avatar Zhou Jiansuo 1 , 2 , avatar Wang Jun ORCID 1 , 3 , *

Department of Transfusion, Peking University Third Hospital, Beijing, China
Department of Laboratory, Peking University Third Hospital, Beijing, China
Department of Anesthesiology, Peking University Third Hospital, Beijing, China

How To Cite Zhenmin S, Nan Y, Jiansuo Z, Jun W. Construction and Evaluation of a Model for the Prediction of Risk Factors Associated with Severe Pneumonia in Patients with Positive Hepatitis B Core Antibody and COVID-19 Infection in China. Jundishapur J Microbiol. 2024;17(8):e148377. https://doi.org/10.5812/jjm-148377.

Abstract

Background:

Hepatitis B virus (HBV) infection is a high-risk factor for severe COVID-19 cases, leading clinicians to enhance the monitoring of patients with concurrent HBV and COVID-19 infections. However, this focus is typically on patients with active HBV infection. Patients who are HBV surface antigen (HBsAg) negative and HBV core antibody (anti-HBc) positive are often considered to have either a past infection that has naturally cleared or to be in a low-replication state, and thus receive less clinical attention. However, in immunosuppressed states, such as during immunosuppressive treatment following COVID-19 infection, the virus in these patients may reactivate. In such cases, this reactivation can increase the risk of COVID-19 progressing to a severe illness. Therefore, being HBsAg negative and anti-HBc positive could be a potential risk factor for severe COVID-19. Studying the combination of anti-HBc positivity and COVID-19 infection can help identify high-risk populations, allowing for the implementation of targeted prevention and management measures, thereby reducing the occurrence of severe COVID-19 cases.

Objectives:

To establish and evaluate a model for predicting severe pneumonia in patients with positive anti-HBc combined with COVID-19 infection.

Methods:

We retrospectively enrolled 380 patients who tested positive for anti-HBc, negative for HBsAg, and HBV e-antigen (HBeAg), combined with COVID-19 infection, in our hospital from December 2022 to May 2023. Based on the inclusion criteria, 163 patients were included in the study. We applied the Lasso binary logistic regression model to optimize feature selection, identifying eight non-zero coefficients using a minimum of one standard error. Using the multiple logistic regression method with backward selection, we screened six factors from the eight selected by the Lasso binary logistic regression model. These six factors were used to construct the predictive model and a nomogram. The validity of our nomogram was assessed using the area under the receiver operating characteristic curve (AUC), calibration curve, Hosmer-Lemeshow (HL) test, and decision-curve analysis (DCA).

Results:

Hypertension, diabetes, decreased absolute lymphocyte count, prolonged prothrombin time, elevated aspartate aminotransferase, and decreased albumin are high-risk factors for severe pneumonia in patients with positive anti-HBc combined with COVID-19 infection. The AUC of the predictive model constructed using these six factors is 0.785, with a 95% confidence interval of (0.709-0.862). The HL test, performed using the calibration curve, yielded a p-value of 0.868. The application of this diagnostic curve will increase the net benefit when the threshold probability is between 5% and 75%.

Conclusions:

The constructed nomogram can be used to predict the risk of patients with positive anti-HBc combined with COVID-19 infection progressing to severe pneumonia, based on routine blood parameters, liver function, coagulation, lactate dehydrogenase levels, and the patient's underlying disease. The predictive model demonstrates good discrimination, calibration, and clinical utility.

1. Background

The COVID-19 pandemic has impacted over 770 million people and resulted in approximately 7 million deaths (1). Although COVID-19 is no longer considered a Public Health Emergency of International Concern, the entire world should reflect on the lessons learned and be prepared for the next public health emergency (2). Predicting the disease course and improving patient prognosis are important for COVID-19 patients, as those with severe COVID-19 experience worse outcomes, including high in-hospital mortality, exacerbation of underlying conditions, reactivation of latent pathogens, and long-term post-acute sequelae of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Therefore, identifying high-risk patients is crucial for facilitating better clinical management (3). Although there are models predicting COVID-19 severity (4-6), it is necessary to develop generalizable severity stratification models for different high-risk groups (such as those with chronic lung diseases, chronic liver diseases, immunodeficiencies, etc.) to guide more personalized treatment plans. This would provide a valuable reference for clinical management during the next pandemic.

Globally, approximately 296 million people are infected with Hepatitis B virus (HBV) (7). Hepatitis B has consistently been the leading cause of morbidity and mortality in mainland China, with the reported number of cases ranking first among all notifiable infectious diseases (8, 9). Individuals who have been previously infected with HBV will have positive HBV core antibody (anti-HBc) in their serum, even if the HBV surface antigen (HBsAg) becomes negative after infection. Patients with positive anti-HBc remain at risk of HBV reactivation or acute liver failure when receiving immunosuppressive therapy (10). Since immunotherapy is one of the treatment regimens for severe COVID-19 patients, it is necessary to develop personalized disease severity prediction models for individuals with positive anti-HBc who are also infected with COVID-19. This would help assess the risk of severe illness in COVID-19 patients with positive anti-HBc and assist in formulating appropriate treatment plans and preventive measures against HBV reactivation.

2. Objectives

In this study, we conducted a retrospective analysis of patients with COVID-19 who were positive for anti-HBc, negative for HBsAg and HBV e-antigen (HBeAg), and were admitted to our hospital from December 2022 to May 2023. The aim was to provide data-supported clinical management plans for patients with COVID-19 and positive anti-HBc to improve the prognosis of this population.

3. Methods

3.1. Study Design

This study is a single-center, retrospective analysis focusing on patients with COVID-19 who were anti-HBc positive, and negative for HBsAg and HBeAg, admitted to Peking University Third Hospital from December 2022 to May 2023. By analyzing patients' routine blood tests, liver function, coagulation parameters, lactate dehydrogenase levels, and underlying diseases, we aimed to determine the risk factors for severe pneumonia and construct a predictive model. The flowchart outlining the study is presented in Figure 1.

Flowchart outlining the study protocol
Flowchart outlining the study protocol

3.2. Case Selection

We enrolled a total of 1,275 cases diagnosed with "novel coronavirus infection" from January 2022 to May 2023, and ultimately selected 163 cases based on the inclusion and exclusion criteria.

Inclusion criteria: Patients who met the diagnostic and classification criteria for novel coronavirus infection outlined in the "Diagnosis and Treatment Protocol for COVID-19" (trial version 10) (11), jointly issued by the General Office of the National Health Commission and the Comprehensive Department of the National Administration of Traditional Chinese Medicine on January 5, 2023; Patients diagnosed with "novel coronavirus infection" at the time of discharge; Patients who had received the COVID-19 vaccination; patients aged over 18 years. For patients with multiple test results during hospitalization, we selected the test result closest to the first day of hospitalization for analysis. For patients without test results during hospitalization, we selected the outpatient test result taken within 3 days closest to the first day of hospitalization for analysis.

Exclusion criteria: Pregnant patients were excluded. The diagnostic criteria for severe COVID-19 (11) are as follows: An adult meets any of the following criteria, with no explanation other than COVID-19 infection: (1) presents with dyspnea, with a respiratory rate ≥ 30 breaths/min; (2) resting state Oxygen saturation ≤ 93% when breathing room air; (3) arterial oxygen partial pressure (PaO2)/fraction of inspired oxygen (FiO2) ≤ 300 mmHg (1 mmHg = 0.133 kPa). For high-altitude areas (above 1000 meters), PaO2/FiO2 should be corrected using the formula: PaO2/FiO2 × [760/atmospheric pressure (mmHg)]; (4) progressive worsening of clinical symptoms, with lung imaging showing a clear progression of the lesion by more than 50% within 24 - 48 hours.

3.3. Data Collection

The following test parameters were selected: Medical history (cardiovascular disease, renal disease, hypertension, diabetes), demographic characteristics (age, gender), white blood cell (WBC) count, red blood cell (RBC) count, absolute lymphocyte count (ALC), absolute neutrophil count (ANC), absolute eosinophil count (AEC), absolute basophil count (ABC), absolute monocyte count (AMC), prothrombin time (PT), activated partial thromboplastin time (APTT), alanine aminotransferase (ALT), aspartate aminotransferase (AST), total bilirubin (T.Bil), alkaline phosphatase (ALP), total protein (TP), albumin (ALB), gamma-glutamyl transferase (γ.GT), globulin (GLB), and lactate dehydrogenase (LDH). The study also included the anti-HBc and HBsAg test results for each patient from 1 year prior to admission up to the time of discharge.

The acquisition of other experimental data was based on the following principles: For patients with multiple test results during hospitalization, we selected the test result closest to the first day of hospitalization for analysis. For patients without any test results during hospitalization, we selected the outpatient test result from within 3 days closest to the first day of hospitalization for analysis. The study was conducted retrospectively using hospital data.

3.4. Statistical Analysis

Data analysis, model building, and evaluation were performed using R version 4.2.1. Continuous variables that followed a normal distribution were expressed as mean ± standard deviation, while non-normally distributed continuous variables were expressed as median (interquartile range), and categorical variables were expressed as frequencies and percentages. The compare groups package was used to perform baseline analysis, with P < 0.05 considered statistically significant. The glmnet package was employed to perform Lasso binary logistic regression, and the minimum criterion with one standard error was used to select non-zero coefficients (12-14). The factors selected by Lasso binary logistic regression were included in a multivariable Logit regression model, with the glm package used for multivariable logistic regression analysis. The minimum Akaike information criterion (AIC) was utilized for model factor selection (15, 16).

The pROC, ggROC, and fbroc packages were used to perform discrimination analysis, drawing receiver operating characteristic (ROC) curves and calculating areas under the curve (AUC). The AUC was used to evaluate the discriminative power of the model. Calibration was performed using the val.prob function and calibrate from the rms package. The Hosmer-Lemeshow (HL) test was performed using the ResourceSelection package. The calibration curve, along with the HL test, was used to compare the observed probabilities and predicted probabilities. Decision curve analysis (DCA) and cross-validated DCA were performed using the rmda package, and DCA was used to evaluate clinical utility. The nomogram was created using the rms package, and the nomogram was used to present the research results.

4. Results

4.1. Statistical Analysis of Study Data

Based on the inclusion and exclusion criteria of the study, a total of 163 patients were included: Forty-four patients with severe COVID-19 pneumonia and 119 patients with non-severe pneumonia. The demographic and clinical characteristics of patients with severe and non-severe pneumonia are shown in Table 1. A total of 24 variables were included in the statistical analysis, and significant differences (P ≤ 0.05) were found between the two groups of patients in experimental parameters such as hypertension, WBC count, ALC, ANC, AEC, PT, AST, TP, ALB, and LDH.

Table 1.

Differences in Demographic and Clinical Characteristics of Patients with Severe and Non-Severe COVID-19 Combined with Positive Anti-HBc and Negative HBsAg and HBeAg a

CharacteristicsAll; (N = 166)Non-severe; (N = 119)Severe; (N = 44)P-Value
Cardiovascular disease0.353
No118 (72.4)89 (74.8)29 (65.9)
Yes45 (27.6)30 (25.2)15 (34.1)
Renal disease0.290
No153 (93.9)110 (92.4)43 (97.7)
Yes10 (6.13)9 (7.56)1 (2.27)
Hypertension0.046
No63 (38.7)52 (43.7)11 (25.0)
Yes100 (61.3)67 (56.3)33 (75.0)
Diabetes0.127
No99 (60.7)77 (64.7)22 (50.0)
Yes64 (39.3)42 (35.3)22 (50.0)
Gender0.473
Male102 (62.6)72 (60.5)30 (68.2)
Female61 (37.4)47 (39.5)14 (31.8)
Age (y)78.0 [69.0; 85.0]77.0 [69.0; 84.5]82.0 [71.8; 87.0]0.117
WBC count (× 109/L)7.53 [5.15; 9.71]6.88 [5.01; 9.04]8.27 [5.94; 10.6]0.046
RBC count (× 1012/L)4.02 (0.74)4.04 (0.74)3.98 (0.75)0.642
ALC (× 109/L)0.92 [0.64; 1.21]0.95 [0.71; 1.26]0.66 [0.34; 1.02]< 0.001
ANC (× 109/L)5.92 [3.67;7.92]5.24 [3.46;7.31]6.94 [5.46;9.20]0.002
AEC (× 109/L)0.02 [0.00; 0.08]0.02 [0.01; 0.09]0.00 [0.00; 0.03]0.001
ABC (× 109/L)0.01 [0.01; 0.02]0.01 [0.01; 0.02]0.01 [0.00; 0.02]0.288
AMC (× 109/L)0.46 [0.31; 0.60]0.47 [0.32; 0.61]0.42 [0.29; 0.52]0.170
PT(s)12.0 [11.2; 12.8]11.7 [11.0; 12.6]12.1 [11.8; 13.1]0.015
APTT (s)30.5 [27.8; 33.6]30.3 [28.0; 33.6]31.0 [26.8; 33.0]0.707
ALT (U/L)20.0 [14.9; 30.5]20.0 [14.7; 30.5]21.5 [15.0; 31.0]0.599
AST (U/L)27.0 [20.0; 39.0]27.0 [20.0; 35.5]33.0 [23.8; 53.2]0.020
T.Bil (μmol/L)9.80 [7.55; 13.8]9.80 [7.50; 13.9]9.60 [7.75; 13.3]0.914
ALP (U/L)73.0 [60.0; 92.0]72.0 [60.5; 97.0]76.4 [59.0; 89.0]0.952
TP (g/L)60.2 [56.6; 67.0]61.8 [57.7; 67.8]57.6 [53.9; 60.7]< 0.001
ALB (g/L)32.5 (4.93)33.3 (4.96)30.3 (4.13)< 0.001
γ.GT (U/L)29.0 [19.0; 43.5]27.0 [19.0; 43.5]32.0 [20.5; 43.2]0.619
GLB (g/L)29.0 [25.0; 31.0]29.0 [26.0; 32.0]28.2 [24.0; 30.9]0.151
LDH (U/L)273 [229; 382]258 [218; 324]326 [274; 468]< 0.001

4.2. Construction of Predictive Model

Lasso binary logistic regression was used for feature selection. Selecting a threshold of 1 standard error allows for choosing the smallest number of feature variables while still maintaining the model's predictive performance, thereby achieving model simplicity and interpretability. All factors were included in the Lasso binary logistic regression for screening, and eight non-zero coefficients were selected using a minimum of one standard error, including hypertension, diabetes, ALC, ANC, PT, AST, TP, and ALB (Figure 2).

Feature selection using LASSO binary logistic regression model with ten-fold cross-validation. A, a coefficient profile plot of the logarithmic (Lambda) sequence, displaying the Log (Lambda) values of 24 features in the LASSO model; B, a ten-fold cross-validation curve of the LASSO model. Using the minimum criterion (left dotted line) and the minimum criterion with one standard error (right dotted line) a vertical dotted line was drawn at the optimal value. The minimum criterion with one standard error was used to select 8 non-zero coefficients. LASSO is the abbreviation for Least Absolute Shrinkage and Selection Operator.
Feature selection using LASSO binary logistic regression model with ten-fold cross-validation. A, a coefficient profile plot of the logarithmic (Lambda) sequence, displaying the Log (Lambda) values of 24 features in the LASSO model; B, a ten-fold cross-validation curve of the LASSO model. Using the minimum criterion (left dotted line) and the minimum criterion with one standard error (right dotted line) a vertical dotted line was drawn at the optimal value. The minimum criterion with one standard error was used to select 8 non-zero coefficients. LASSO is the abbreviation for Least Absolute Shrinkage and Selection Operator.

Multiple-factor backward logistic regression was used to construct a predictive model for patients with COVID-19 who were positive for anti-HBc and negative for HBsAg and HBeAg. The eight LASSO-selected factors were incorporated into a multivariable regression model. Using the minimum AIC, six factors were selected to construct the final model. Hypertension, diabetes, ALC, PT, AST, and ALB were finally included in the model (Table 2). A nomogram of the predictive model was plotted (Figure 3), and the total score of each patient was calculated based on the corresponding scores of different risk factors. Finally, the risk of a patient developing severe COVID-19 was determined based on the total score.

Table 2.

Multivariable Backward Logistic Regression Analysis of Patients with Severe and Non-severe COVID-19 Testing Positive for Anti-HBc and Negative for HBsAg and HBeAg

CharacteristicsBSEORCIZP
Intercept-0.7112.0680.4910.491 (0.008 - 27.46)-0.3440.731
Hypertension0.7280.4412.0712.071 (0.891 - 5.084)1.6520.099
Diabetes0.7140.4192.0412.041 (0.902 - 4.703)1.7030.089
ALC-1.1400.5430.3200.320 (0.105 - 0.886)-2.0970.036
PT0.2130.1071.2381.238 (1.018 - 1.559)1.9960.046
AST0.0140.0081.0141.014 (0.998 - 1.032)1.7120.087
ALB-0.1030.0450.9020.902 (0.823 - 0.985)-2.2750.023
Nomogram for risk prediction of severe COVID-19 infection among patients testing positive for anti-HBc and negative for HBsAg and HBeAg. ALC, absolute lymphocyte count; PT, prothrombin time; AST, aspartate aminotransferase; ALB, albumin.
Nomogram for risk prediction of severe COVID-19 infection among patients testing positive for anti-HBc and negative for HBsAg and HBeAg. ALC, absolute lymphocyte count; PT, prothrombin time; AST, aspartate aminotransferase; ALB, albumin.

4.3. Analysis and Validation of the Constructed Model

First, using the ROC curve, we assessed the model's ability to predict whether a patient will develop severe COVID-19 infection. The AUC of the predictive model was plotted, with the AUC of the model being 0.785 (95% CI: 0.709 - 0.862), and the AUC obtained from the 500-Bootstrap method for internal validation was 0.785 (95% CI: 0.717 - 0.864), indicating that the model has good discrimination (Figure 4). We evaluated and visualized the disparity between the model’s predicted values and actual values. The goodness of fit was assessed using the calibration curve, and the calibrated curve, based on the results of the 500-Bootstrap method for internal validation, was consistent with both the original and the ideal calibration curves (Figure 5).

The areas under the curve (AUC) values of the predictive model and the AUC value of the internal validation using the bootstrap method. A, the AUC value of the predictive model was 0.785, with a 95% confidence interval of (0.709 - 0.862); B, the AUC value of the model validated internally via 500-bootstrap method was 0.785, with a 95% confidence interval of (0.717 - 0.864).
The areas under the curve (AUC) values of the predictive model and the AUC value of the internal validation using the bootstrap method. A, the AUC value of the predictive model was 0.785, with a 95% confidence interval of (0.709 - 0.862); B, the AUC value of the model validated internally via 500-bootstrap method was 0.785, with a 95% confidence interval of (0.717 - 0.864).
The calibration curve of the predictive model, which was internally validated by the bootstrap sampling method. The black dashed line represents the ideal calibration curve, the blue solid line represents the original calibration curve, and the red solid line represents the calibration curve using 500-bootstrap method.
The calibration curve of the predictive model, which was internally validated by the bootstrap sampling method. The black dashed line represents the ideal calibration curve, the blue solid line represents the original calibration curve, and the red solid line represents the calibration curve using 500-bootstrap method.

The horizontal axis represents predicted values, while the vertical axis represents actual values. The black dashed line represents the ideal calibration curve, the blue solid line labeled “Apparent” represents the fit between the model’s predicted values and actual values, and the red solid line labeled “Bias-corrected” indicates the fit between the model’s predicted values and actual values after calibration. The shape of the curve suggests good consistency between the model's predicted values and actual values. The calibration curve was further tested using the HL test, and the P-value was 0.868 (P > 0.05), indicating a high goodness of fit. The clinical applicability of the predictive model was evaluated via DCA (Figure 6).

The decision-curve analysis (DCA) and the 5-fold cross-validated DCA of the predictive model. A, the DCA of the predictive model using the 500-bootstrap method. The red line represents the diagnostic curve, the light solid line represents the overall positives, and the dark solid line represents the overall negatives; B, a five-fold cross-validation of the DCA. The red dashed line represents the DCA of the predictive model; the blue solid line is the DCA of cross-validation; the light solid line indicates overall positives, and the dark solid line represents overall negatives.
The decision-curve analysis (DCA) and the 5-fold cross-validated DCA of the predictive model. A, the DCA of the predictive model using the 500-bootstrap method. The red line represents the diagnostic curve, the light solid line represents the overall positives, and the dark solid line represents the overall negatives; B, a five-fold cross-validation of the DCA. The red dashed line represents the DCA of the predictive model; the blue solid line is the DCA of cross-validation; the light solid line indicates overall positives, and the dark solid line represents overall negatives.

When the threshold probability ranged between 5% and 75%, the net benefit was enhanced using this diagnostic curve. The DCA was also validated via 5-fold cross-validation, showing a high degree of consistency with the predictive model, indicating high model performance. Finally, a rationality analysis of the predictive model (Figure 7) showed that the AUC of the predictive model was higher than that of any single predictive factor. The constructed model demonstrated superior performance in predicting whether a patient will develop severe COVID-19 infection compared to a single factor alone. Additionally, the DCA of the predictive model was higher than that of a single factor. The results demonstrated that, within the model’s threshold probability range, it provides a better net benefit to the clinic compared to a single-factor model.

The goodness-of-fit analysis of the predictive model. A, comparison of AUC of the predictive model and individual factors; B, comparison of decision curves using the predictive model and individual factors
The goodness-of-fit analysis of the predictive model. A, comparison of AUC of the predictive model and individual factors; B, comparison of decision curves using the predictive model and individual factors

5. Discussion

This study developed a personalized nomogram for predicting COVID-19 infection in patients testing positive for anti-HBc and negative for HBsAg and HBeAg antibodies. We carried out internal validation of the model. We identified six easily assessable variables, including hypertension, diabetes, ALC, PT, AST, and ALB, which can be used for clinical prediction of the risk of COVID-19 in patients testing positive for anti-HBc and negative for HBsAg and HBeAg developing into severe pneumonia. The model demonstrated good discrimination, calibration, and clinical applicability in predicting progression to severe pneumonia. In addition, the DCA showed that the model can be used to make significant clinical decisions across a wide range of probability thresholds.

According to the European Association for the Study of the Liver (EASL), HBV infection can progress in five stages that are not necessarily consecutive. In some patients, an HBsAg-negative phase may occur, known as "occult HBV infection," which is defined as the absence of serum HBsAg but the presence of anti-HBc, with or without anti-HBs, and with or without detectable serum HBV DNA. However, patients in this stage may have covalently closed circular DNA (cccDNA) in their liver. These patients may experience reactivation of HBV infection if they receive immunosuppressive therapy (17).

The American Gastroenterological Association (AGA) Institute has published guidelines for the prevention and management of HBV reactivation during immunosuppressive therapy. According to the dosage and duration of immunosuppressive medication use, the risk of HBV reactivation is classified into high, moderate, and low levels. For moderate-risk patients, the AGA suggests that antiviral HBV treatment should be continued for 6 months after discontinuation of immunosuppressive therapy. For low-risk patients, the AGA recommends against routinely using prophylactic anti-HBV treatment but advises monitoring HBV DNA levels during immunosuppressive therapy to detect HBV reactivation and identify liver function abnormalities as early as possible (18). However, due to the strain the pandemic has placed on the healthcare system, it cannot be ensured that low-risk patients will receive proper follow-up care after discharge, especially elderly patients returning to long-term care facilities. Moderate-risk patients who do not receive prophylactic treatment are at risk of HBV reactivation. Table 1 shows that the majority of admitted COVID-19 patients are elderly males, which is also a risk factor for HBV reactivation (19). Severe acute respiratory syndrome coronavirus 2 infection can cause significant lymphocyte depletion, increasing the likelihood of HBV reactivation (20). Therefore, evaluating the risk factors for severe COVID-19 in these patients can help develop a reasonable immunotherapy plan, as well as HBV prophylaxis and monitoring regimens, which can improve the prognosis of patients with anti-HBc positivity concurrent with COVID-19 infection. At the same time, this reminds clinicians to pay more attention to patients with positive anti-HBc.

Therefore, the management of anti-HBc-positive patients should be strengthened to reduce the impact of HBV on liver function. Abnormal liver function is a risk factor for mortality in patients with COVID-19 (21). The results of univariate analysis showed that AST levels were elevated in severe COVID-19 cases compared to non-severe cases (P ≤ 0.05), and TP and ALB levels were decreased in severe patients compared to non-severe cases (P ≤ 0.05), which is consistent with the results of other studies (21, 22). This suggests that further quantitative testing of anti-HBc-positive patients for Hepatitis B core antibody is necessary to determine whether the patient has occult Hepatitis B. Some studies have suggested that the use of angiotensin-converting enzyme inhibitors may increase the risk of severe COVID-19. Our model identified hypertension and diabetes as risk factors for severe COVID-19, which is consistent with the results of other studies (23, 24). This suggests that alternative treatment options, such as calcium channel blockers, may be considered for these patients (25). A decrease in ALC is a high-risk factor for severe COVID-19, which is consistent with other research findings (26, 27).

The diagnosis and treatment plan for novel coronavirus infection (10th edition, trial implementation), jointly issued by the General Office of the National Health Commission and the Comprehensive Department of the National Administration of Traditional Chinese Medicine on January 5, 2023, pointed out that a decrease in lymphocyte count is an early warning sign of severe/critical COVID-19. Thrombocytopenia, moderately prolonged PT, and elevated D-dimer are coagulation abnormalities associated with COVID-19 infection. This may be due to elevated von Willebrand factor levels in patients with COVID-19, leading to a deficiency of the von Willebrand factor-cleaving protease ADAMTS13 (A disintegrin and metalloproteinase with thrombospondin motifs 13), resulting in platelet aggregation and microthrombus formation, followed by activation of the fibrinolysis system and elevation of D-dimer levels (28). This study also evaluated factors such as occupation, education level, marital status, and alcohol consumption, and found that alcohol consumption increased the risk of severe COVID-19 by 3.52-fold (29). In the future, we may include additional parameters and data to construct a more comprehensive model for clinical prediction and conduct multicenter studies to externally validate the model for improved predictive accuracy.

5.1. Conclusions

Our model, constructed using routine clinical laboratory tests, achieved an AUC of 0.785 and demonstrated good calibration, clinical applicability, and rationality. The nomogram we developed using only six variables achieved high accuracy, facilitating rapid clinical detection of patient risk factors and the development of individualized treatment plans. However, this study has limitations. It is a retrospective single-center analysis, and the sample selection only included patients from our hospital, resulting in a small sample size. In future studies, it would be beneficial to design multicenter experiments to further validate the aforementioned research findings. In the subsequent multicenter study, external validation of the model results can be conducted using data from other hospitals. This will help to reduce the model's bias and enhance its generalizability.

References

  • 1.

    Falzon D, Zignol M, Bastard M, Floyd K, Kasaeva T. The impact of the COVID-19 pandemic on the global tuberculosis epidemic. Front Immunol. 2023;14:1234785. [PubMed ID: 37795102]. [PubMed Central ID: PMC10546619]. https://doi.org/10.3389/fimmu.2023.1234785.

  • 2.

    Frutos R, Gavotte L, Serra-Cobo J, Chen T, Devaux C. COVID-19 and emerging infectious diseases: The society is still unprepared for the next pandemic. Environ Res. 2021;202:111676. [PubMed ID: 34252435]. [PubMed Central ID: PMC8268624]. https://doi.org/10.1016/j.envres.2021.111676.

  • 3.

    Davis HE, McCorkell L, Vogel JM, Topol EJ. Long COVID: Major findings, mechanisms and recommendations. Nat Rev Microbiol. 2023;21(3):133-46. [PubMed ID: 36639608]. [PubMed Central ID: PMC9839201]. https://doi.org/10.1038/s41579-022-00846-2.

  • 4.

    Feyaerts D, Hedou J, Gillard J, Chen H, Tsai ES, Peterson LS, et al. Integrated plasma proteomic and single-cell immune signaling network signatures demarcate mild, moderate, and severe COVID-19. Cell Rep Med. 2022;3(7):100680. [PubMed ID: 35839768]. [PubMed Central ID: PMC9238057]. https://doi.org/10.1016/j.xcrm.2022.100680.

  • 5.

    Liu A, Hammond R, Donnelly PD, Kaski JC, Coates ARM. Effective prognostic and clinical risk stratification in COVID-19 using multimodality biomarkers. J Intern Med. 2023;294(1):21-46. [PubMed ID: 37106509]. https://doi.org/10.1111/joim.13646.

  • 6.

    Knight SR, Ho A, Pius R, Buchan I, Carson G, Drake TM, et al. Risk stratification of patients admitted to hospital with covid-19 using the ISARIC WHO Clinical Characterisation Protocol: development and validation of the 4C Mortality Score. BMJ. 2020;370:m3339. [PubMed ID: 32907855]. [PubMed Central ID: PMC7116472]. https://doi.org/10.1136/bmj.m3339.

  • 7.

    Hsu YC, Huang DQ, Nguyen MH. Global burden of hepatitis B virus: Current status, missed opportunities and a call for action. Nat Rev Gastroenterol Hepatol. 2023;20(8):524-37. [PubMed ID: 37024566]. https://doi.org/10.1038/s41575-023-00760-9.

  • 8.

    Su S, Wong WC, Zou Z, Cheng DD, Ong JJ, Chan P, et al. Cost-effectiveness of universal screening for chronic hepatitis B virus infection in China: An economic evaluation. Lancet Glob Health. 2022;10(2):e278-87. [PubMed ID: 35063115]. [PubMed Central ID: PMC8789560]. https://doi.org/10.1016/S2214-109X(21)00517-9.

  • 9.

    Liu J, Liang W, Jing W, Liu M. Countdown to 2030: Eliminating Hepatitis B disease, China. Bull World Health Organ. 2019;97(3):230-8. [PubMed ID: 30992636]. [PubMed Central ID: PMC6453311]. https://doi.org/10.2471/BLT.18.219469.

  • 10.

    Rodriguez-Tajes S, Miralpeix A, Costa J, Lopez-Sune E, Laguno M, Pocurull A, et al. Low risk of Hepatitis B reactivation in patients with severe COVID-19 who receive immunosuppressive therapy. J Viral Hepat. 2021;28(1):89-94. [PubMed ID: 32969557]. [PubMed Central ID: PMC7537127]. https://doi.org/10.1111/jvh.13410.

  • 11.

    Wang SX, Wang Y, Lu YB, Li JY, Song YJ, Nyamgerelt M, et al. Diagnosis and treatment of novel coronavirus pneumonia based on the theory of traditional Chinese medicine. J Integr Med. 2020;18(4):275-83. [PubMed ID: 32446813]. [PubMed Central ID: PMC7159848]. https://doi.org/10.1016/j.joim.2020.04.001.

  • 12.

    Shang Y, Liu T, Wei Y, Li J, Shao L, Liu M, et al. Scoring systems for predicting mortality for severe patients with COVID-19. E Clinic Med. 2020;24:100426. [PubMed ID: 32766541]. [PubMed Central ID: PMC7332889]. https://doi.org/10.1016/j.eclinm.2020.100426.

  • 13.

    Chen Y, Qi Y, Li T, Lin A, Ni Y, Pu R, et al. A more objective PD diagnostic model: Integrating texture feature markers of cerebellar gray matter and white matter through machine learning. Front Aging Neurosci. 2024;16:1393841. [PubMed ID: 38912523]. [PubMed Central ID: PMC11190310]. https://doi.org/10.3389/fnagi.2024.1393841.

  • 14.

    Kang H, Zhao D, Xiang H, Li J, Zhao G, Li H. Large-scale transcriptome sequencing in broiler chickens to identify candidate genes for breast muscle weight and intramuscular fat content. Genet Sel Evol. 2021;53(1):66. [PubMed ID: 34399688]. [PubMed Central ID: PMC8369645]. https://doi.org/10.1186/s12711-021-00656-9.

  • 15.

    Krull KR, Brinkman TM, Li C, Armstrong GT, Ness KK, Srivastava DK, et al. Neurocognitive outcomes decades after treatment for childhood acute lymphoblastic leukemia: A report from the St Jude lifetime cohort study. J Clin Oncol. 2013;31(35):4407-15. [PubMed ID: 24190124]. [PubMed Central ID: PMC3842908]. https://doi.org/10.1200/JCO.2012.48.2315.

  • 16.

    Domb BG, Ouyang VW, Go CC, Gornbein JA, Shapira J, Meghpara MB, et al. Personalized medicine using predictive analytics: A machine learning-based prognostic model for patients undergoing hip arthroscopy. Am J Sports Med. 2022;50(7):1900-8. [PubMed ID: 35536218]. https://doi.org/10.1177/03635465221091847.

  • 17.

    European Association for the Study of the Liver. EASL 2017 clinical practice guidelines on the management of Hepatitis B virus infection. J Hepatol. 2017;67(2):370-98. [PubMed ID: 28427875]. https://doi.org/10.1016/j.jhep.2017.03.021.

  • 18.

    Reddy KR, Beavers KL, Hammond SP, Lim JK, Falck-Ytter YT, American Gastroenterological Association I. American Gastroenterological Association Institute guideline on the prevention and treatment of hepatitis B virus reactivation during immunosuppressive drug therapy. Gastroenterol. 2015;148(1):215-9. quiz e16-7. [PubMed ID: 25447850]. https://doi.org/10.1053/j.gastro.2014.10.039.

  • 19.

    Loomba R, Rowley A, Wesley R, Liang TJ, Hoofnagle JH, Pucino F, et al. Systematic review: the effect of preventive lamivudine on hepatitis B reactivation during chemotherapy. Ann Intern Med. 2008;148(7):519-28. [PubMed ID: 18378948]. [PubMed Central ID: PMC2729097]. https://doi.org/10.7326/0003-4819-148-7-200804010-00008.

  • 20.

    Maini MK, Gehring AJ. The role of innate immunity in the immunopathology and treatment of HBV infection. J Hepatol. 2016;64(1):S60-70. https://doi.org/10.1016/j.jhep.2016.01.028.

  • 21.

    Lei F, Liu YM, Zhou F, Qin JJ, Zhang P, Zhu L, et al. Longitudinal association between markers of liver injury and mortality in COVID-19 in China. Hepatol. 2020;72(2):389-98. [PubMed ID: 32359177]. [PubMed Central ID: PMC7267515]. https://doi.org/10.1002/hep.31301.

  • 22.

    Wang J, Lu Z, Jin M, Wang Y, Tian K, Xiao J, et al. Clinical characteristics and risk factors of COVID-19 patients with chronic hepatitis B: A multi-center retrospective cohort study. Front Med. 2022;16(1):111-25. [PubMed ID: 34387851]. [PubMed Central ID: PMC8362646]. https://doi.org/10.1007/s11684-021-0854-5.

  • 23.

    Mahamat-Saleh Y, Fiolet T, Rebeaud ME, Mulot M, Guihur A, El Fatouhi D, et al. Diabetes, hypertension, body mass index, smoking and COVID-19-related mortality: a systematic review and meta-analysis of observational studies. BMJ Open. 2021;11(10). e052777. [PubMed ID: 34697120]. [PubMed Central ID: PMC8557249]. https://doi.org/10.1136/bmjopen-2021-052777.

  • 24.

    Schlesinger S, Neuenschwander M, Lang A, Pafili K, Kuss O, Herder C, et al. Risk phenotypes of diabetes and association with COVID-19 severity and death: A living systematic review and meta-analysis. Diabetol. 2021;64(7):1480-91. [PubMed ID: 33907860]. [PubMed Central ID: PMC8079163]. https://doi.org/10.1007/s00125-021-05458-8.

  • 25.

    Fang L, Karakiulakis G, Roth M. Are patients with hypertension and diabetes mellitus at increased risk for COVID-19 infection? Lancet Respir Med. 2020;8(4). e21. [PubMed ID: 32171062]. [PubMed Central ID: PMC7118626]. https://doi.org/10.1016/S2213-2600(20)30116-8.

  • 26.

    Illg Z, Muller G, Mueller M, Nippert J, Allen B. Analysis of absolute lymphocyte count in patients with COVID-19. Am J Emerg Med. 2021;46:16-9. [PubMed ID: 33706251]. [PubMed Central ID: PMC7923864]. https://doi.org/10.1016/j.ajem.2021.02.054.

  • 27.

    Wagner J, DuPont A, Larson S, Cash B, Farooq A. Absolute lymphocyte count is a prognostic marker in Covid-19: A retrospective cohort review. Int J Lab Hematol. 2020;42(6):761-5. [PubMed ID: 32779838]. [PubMed Central ID: PMC7405282]. https://doi.org/10.1111/ijlh.13288.

  • 28.

    D'Ardes D, Boccatonda A, Cocco G, Fabiani S, Rossi I, Bucci M, et al. Impaired coagulation, liver dysfunction and COVID-19: Discovering an intriguing relationship. World J Gastroenterol. 2022;28(11):1102-12. [PubMed ID: 35431501]. [PubMed Central ID: PMC8985482]. https://doi.org/10.3748/wjg.v28.i11.1102.

  • 29.

    Tekle E, Gelaw Y, Dagnew M, Gelaw A, Negash M, Kassa E, et al. Risk stratification and prognostic value of prothrombin time and activated partial thromboplastin time among COVID-19 patients. PLoS One. 2022;17(8). e0272216. [PubMed ID: 35951632]. [PubMed Central ID: PMC9371343]. https://doi.org/10.1371/journal.pone.0272216.