Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Development and validation of the predictive score for pediatric COVID-19 pneumonia: A nationwide, multicenter study

  • Araya Satdhabudha,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Writing – original draft, Writing – review & editing

    Affiliation Department of Pediatrics, Faculty of Medicine, Thammasat University Hospital, Thammasat University, Pathumthani, Thailand

  • Chanapai Chaiyakulsil ,

    Roles Data curation, Formal analysis, Investigation, Methodology, Project administration, Writing – original draft, Writing – review & editing

    chanapai.chai@hotmail.com

    Affiliation Department of Pediatrics, Faculty of Medicine, Thammasat University Hospital, Thammasat University, Pathumthani, Thailand

  • Rattapon Uppala ,

    Contributed equally to this work with: Rattapon Uppala, Watit Niyomkarn, Prakarn Tovichien, Vasinee Norasettekul, Kanokpan Ruangnapa, Chutima Smathakanee, Bararee Choursamran, Aunya Kulbun, Rasintra Jaroenying, Harutai Kamalaporn, Tidarat Sriboonyong, Koonkoaw Roekworachai, Kanokkarn Sunkonkit, Auchara Tangsathapornpong, Pornumpa Bunjoungmanee, Patcharapa Thaweekul, Thanyarat Ratanavongkosol, Chutima Thongnual

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Pediatrics, Faculty of Medicine, Khon Kaen University, Khon Kaen, Thailand

  • Watit Niyomkarn ,

    Contributed equally to this work with: Rattapon Uppala, Watit Niyomkarn, Prakarn Tovichien, Vasinee Norasettekul, Kanokpan Ruangnapa, Chutima Smathakanee, Bararee Choursamran, Aunya Kulbun, Rasintra Jaroenying, Harutai Kamalaporn, Tidarat Sriboonyong, Koonkoaw Roekworachai, Kanokkarn Sunkonkit, Auchara Tangsathapornpong, Pornumpa Bunjoungmanee, Patcharapa Thaweekul, Thanyarat Ratanavongkosol, Chutima Thongnual

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Pediatrics, Faculty of Medicine, Chulalongkorn University, Bangkok, Thailand

  • Prakarn Tovichien ,

    Contributed equally to this work with: Rattapon Uppala, Watit Niyomkarn, Prakarn Tovichien, Vasinee Norasettekul, Kanokpan Ruangnapa, Chutima Smathakanee, Bararee Choursamran, Aunya Kulbun, Rasintra Jaroenying, Harutai Kamalaporn, Tidarat Sriboonyong, Koonkoaw Roekworachai, Kanokkarn Sunkonkit, Auchara Tangsathapornpong, Pornumpa Bunjoungmanee, Patcharapa Thaweekul, Thanyarat Ratanavongkosol, Chutima Thongnual

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Pediatrics, Faculty of Medicine Siriraj Hospital, Mahidol University, Bangkok, Thailand

  • Vasinee Norasettekul ,

    Contributed equally to this work with: Rattapon Uppala, Watit Niyomkarn, Prakarn Tovichien, Vasinee Norasettekul, Kanokpan Ruangnapa, Chutima Smathakanee, Bararee Choursamran, Aunya Kulbun, Rasintra Jaroenying, Harutai Kamalaporn, Tidarat Sriboonyong, Koonkoaw Roekworachai, Kanokkarn Sunkonkit, Auchara Tangsathapornpong, Pornumpa Bunjoungmanee, Patcharapa Thaweekul, Thanyarat Ratanavongkosol, Chutima Thongnual

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Pediatrics, Charoenkrung Pracharak Hospital, Bangkok, Thailand

  • Kanokpan Ruangnapa ,

    Contributed equally to this work with: Rattapon Uppala, Watit Niyomkarn, Prakarn Tovichien, Vasinee Norasettekul, Kanokpan Ruangnapa, Chutima Smathakanee, Bararee Choursamran, Aunya Kulbun, Rasintra Jaroenying, Harutai Kamalaporn, Tidarat Sriboonyong, Koonkoaw Roekworachai, Kanokkarn Sunkonkit, Auchara Tangsathapornpong, Pornumpa Bunjoungmanee, Patcharapa Thaweekul, Thanyarat Ratanavongkosol, Chutima Thongnual

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Pediatrics, Faculty of Medicine, Prince of Songkla University, Songkhla, Thailand

  • Chutima Smathakanee ,

    Contributed equally to this work with: Rattapon Uppala, Watit Niyomkarn, Prakarn Tovichien, Vasinee Norasettekul, Kanokpan Ruangnapa, Chutima Smathakanee, Bararee Choursamran, Aunya Kulbun, Rasintra Jaroenying, Harutai Kamalaporn, Tidarat Sriboonyong, Koonkoaw Roekworachai, Kanokkarn Sunkonkit, Auchara Tangsathapornpong, Pornumpa Bunjoungmanee, Patcharapa Thaweekul, Thanyarat Ratanavongkosol, Chutima Thongnual

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Pediatrics, Hat Yai Hospital, Songkhla, Thailand

  • Bararee Choursamran ,

    Contributed equally to this work with: Rattapon Uppala, Watit Niyomkarn, Prakarn Tovichien, Vasinee Norasettekul, Kanokpan Ruangnapa, Chutima Smathakanee, Bararee Choursamran, Aunya Kulbun, Rasintra Jaroenying, Harutai Kamalaporn, Tidarat Sriboonyong, Koonkoaw Roekworachai, Kanokkarn Sunkonkit, Auchara Tangsathapornpong, Pornumpa Bunjoungmanee, Patcharapa Thaweekul, Thanyarat Ratanavongkosol, Chutima Thongnual

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Pediatrics, Sawanpracharak Hospital, Nakhon Sawan, Thailand

  • Aunya Kulbun ,

    Contributed equally to this work with: Rattapon Uppala, Watit Niyomkarn, Prakarn Tovichien, Vasinee Norasettekul, Kanokpan Ruangnapa, Chutima Smathakanee, Bararee Choursamran, Aunya Kulbun, Rasintra Jaroenying, Harutai Kamalaporn, Tidarat Sriboonyong, Koonkoaw Roekworachai, Kanokkarn Sunkonkit, Auchara Tangsathapornpong, Pornumpa Bunjoungmanee, Patcharapa Thaweekul, Thanyarat Ratanavongkosol, Chutima Thongnual

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Pediatrics, Her Royal Highness Maha Chakri Sirindhorn Medical Center, Faculty of Medicine, Srinakharinwirot University, Nakhon Nayok, Thailand

  • Rasintra Jaroenying ,

    Contributed equally to this work with: Rattapon Uppala, Watit Niyomkarn, Prakarn Tovichien, Vasinee Norasettekul, Kanokpan Ruangnapa, Chutima Smathakanee, Bararee Choursamran, Aunya Kulbun, Rasintra Jaroenying, Harutai Kamalaporn, Tidarat Sriboonyong, Koonkoaw Roekworachai, Kanokkarn Sunkonkit, Auchara Tangsathapornpong, Pornumpa Bunjoungmanee, Patcharapa Thaweekul, Thanyarat Ratanavongkosol, Chutima Thongnual

    Roles Data curation, Investigation, Writing – review & editing

    Affiliations Department of Pediatrics, Faculty of Medicine, Phramongkutklao Hospital, Bangkok, Thailand, Phramongkutklao Hospital Sleep Disorder Center, Faculty of Medicine, Phramongkutklao Hospital, Bangkok, Thailand

  • Harutai Kamalaporn ,

    Contributed equally to this work with: Rattapon Uppala, Watit Niyomkarn, Prakarn Tovichien, Vasinee Norasettekul, Kanokpan Ruangnapa, Chutima Smathakanee, Bararee Choursamran, Aunya Kulbun, Rasintra Jaroenying, Harutai Kamalaporn, Tidarat Sriboonyong, Koonkoaw Roekworachai, Kanokkarn Sunkonkit, Auchara Tangsathapornpong, Pornumpa Bunjoungmanee, Patcharapa Thaweekul, Thanyarat Ratanavongkosol, Chutima Thongnual

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Pediatrics, Faculty of Medicine, Ramathibodi Hospital, Mahidol University, Bangkok, Thailand

  • Tidarat Sriboonyong ,

    Contributed equally to this work with: Rattapon Uppala, Watit Niyomkarn, Prakarn Tovichien, Vasinee Norasettekul, Kanokpan Ruangnapa, Chutima Smathakanee, Bararee Choursamran, Aunya Kulbun, Rasintra Jaroenying, Harutai Kamalaporn, Tidarat Sriboonyong, Koonkoaw Roekworachai, Kanokkarn Sunkonkit, Auchara Tangsathapornpong, Pornumpa Bunjoungmanee, Patcharapa Thaweekul, Thanyarat Ratanavongkosol, Chutima Thongnual

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Pediatrics, Faculty of Medicine, Ramathibodi Hospital, Mahidol University, Bangkok, Thailand

  • Koonkoaw Roekworachai ,

    Contributed equally to this work with: Rattapon Uppala, Watit Niyomkarn, Prakarn Tovichien, Vasinee Norasettekul, Kanokpan Ruangnapa, Chutima Smathakanee, Bararee Choursamran, Aunya Kulbun, Rasintra Jaroenying, Harutai Kamalaporn, Tidarat Sriboonyong, Koonkoaw Roekworachai, Kanokkarn Sunkonkit, Auchara Tangsathapornpong, Pornumpa Bunjoungmanee, Patcharapa Thaweekul, Thanyarat Ratanavongkosol, Chutima Thongnual

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Pediatrics, Nakornping Hospital, Chiangmai, Thailand

  • Kanokkarn Sunkonkit ,

    Contributed equally to this work with: Rattapon Uppala, Watit Niyomkarn, Prakarn Tovichien, Vasinee Norasettekul, Kanokpan Ruangnapa, Chutima Smathakanee, Bararee Choursamran, Aunya Kulbun, Rasintra Jaroenying, Harutai Kamalaporn, Tidarat Sriboonyong, Koonkoaw Roekworachai, Kanokkarn Sunkonkit, Auchara Tangsathapornpong, Pornumpa Bunjoungmanee, Patcharapa Thaweekul, Thanyarat Ratanavongkosol, Chutima Thongnual

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Pediatrics, Faculty of Medicine, Chiang Mai University, Chiang Mai, Thailand

  • Auchara Tangsathapornpong ,

    Contributed equally to this work with: Rattapon Uppala, Watit Niyomkarn, Prakarn Tovichien, Vasinee Norasettekul, Kanokpan Ruangnapa, Chutima Smathakanee, Bararee Choursamran, Aunya Kulbun, Rasintra Jaroenying, Harutai Kamalaporn, Tidarat Sriboonyong, Koonkoaw Roekworachai, Kanokkarn Sunkonkit, Auchara Tangsathapornpong, Pornumpa Bunjoungmanee, Patcharapa Thaweekul, Thanyarat Ratanavongkosol, Chutima Thongnual

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Pediatrics, Faculty of Medicine, Thammasat University Hospital, Thammasat University, Pathumthani, Thailand

  • Pornumpa Bunjoungmanee ,

    Contributed equally to this work with: Rattapon Uppala, Watit Niyomkarn, Prakarn Tovichien, Vasinee Norasettekul, Kanokpan Ruangnapa, Chutima Smathakanee, Bararee Choursamran, Aunya Kulbun, Rasintra Jaroenying, Harutai Kamalaporn, Tidarat Sriboonyong, Koonkoaw Roekworachai, Kanokkarn Sunkonkit, Auchara Tangsathapornpong, Pornumpa Bunjoungmanee, Patcharapa Thaweekul, Thanyarat Ratanavongkosol, Chutima Thongnual

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Pediatrics, Faculty of Medicine, Thammasat University Hospital, Thammasat University, Pathumthani, Thailand

  • Wanida Pao-in,

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Pediatrics, Faculty of Medicine, Thammasat University Hospital, Thammasat University, Pathumthani, Thailand

  • Patcharapa Thaweekul ,

    Contributed equally to this work with: Rattapon Uppala, Watit Niyomkarn, Prakarn Tovichien, Vasinee Norasettekul, Kanokpan Ruangnapa, Chutima Smathakanee, Bararee Choursamran, Aunya Kulbun, Rasintra Jaroenying, Harutai Kamalaporn, Tidarat Sriboonyong, Koonkoaw Roekworachai, Kanokkarn Sunkonkit, Auchara Tangsathapornpong, Pornumpa Bunjoungmanee, Patcharapa Thaweekul, Thanyarat Ratanavongkosol, Chutima Thongnual

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Pediatrics, Faculty of Medicine, Thammasat University Hospital, Thammasat University, Pathumthani, Thailand

  • Pichaya Tantiyavarong,

    Roles Formal analysis, Investigation, Methodology, Writing – review & editing

    Affiliation Department of Clinical Epidemiology, Faculty of Medicine, Thammasat University, Pathumthani, Thailand

  • Thanyarat Ratanavongkosol ,

    Contributed equally to this work with: Rattapon Uppala, Watit Niyomkarn, Prakarn Tovichien, Vasinee Norasettekul, Kanokpan Ruangnapa, Chutima Smathakanee, Bararee Choursamran, Aunya Kulbun, Rasintra Jaroenying, Harutai Kamalaporn, Tidarat Sriboonyong, Koonkoaw Roekworachai, Kanokkarn Sunkonkit, Auchara Tangsathapornpong, Pornumpa Bunjoungmanee, Patcharapa Thaweekul, Thanyarat Ratanavongkosol, Chutima Thongnual

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Pediatrics, Charoenkrung Pracharak Hospital, Bangkok, Thailand

  • Chutima Thongnual ,

    Contributed equally to this work with: Rattapon Uppala, Watit Niyomkarn, Prakarn Tovichien, Vasinee Norasettekul, Kanokpan Ruangnapa, Chutima Smathakanee, Bararee Choursamran, Aunya Kulbun, Rasintra Jaroenying, Harutai Kamalaporn, Tidarat Sriboonyong, Koonkoaw Roekworachai, Kanokkarn Sunkonkit, Auchara Tangsathapornpong, Pornumpa Bunjoungmanee, Patcharapa Thaweekul, Thanyarat Ratanavongkosol, Chutima Thongnual

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Pediatrics, Hat Yai Hospital, Songkhla, Thailand

  • Paskorn Sritipsukho,

    Roles Formal analysis, Methodology, Project administration, Writing – review & editing

    Affiliations Department of Pediatrics, Faculty of Medicine, Thammasat University Hospital, Thammasat University, Pathumthani, Thailand, Center of Excellence in Applied Epidemiology, Thammasat University, Pathumthani, Thailand

  •  [ ... ],
  • Jitladda Deerojanawong

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Pediatrics, Faculty of Medicine, Chulalongkorn University, Bangkok, Thailand

  • [ view all ]
  • [ view less ]

Abstract

Background

Due to the possibility of asymptomatic pneumonia in children with COVID-19 leading to overexposure to radiation and problems in limited-resource settings, we conducted a nationwide, multi-center study to determine the risk factors of pneumonia in children with COVID-19 in order to create a pediatric pneumonia predictive score, with score validation.

Methods

This was a retrospective cohort study done by chart review of all children aged 0–15 years admitted to 13 medical centers across Thailand during the study period. Univariate and multivariate analyses as well as backward and forward stepwise logistic regression were used to generate a final prediction model of the pneumonia score. Data during the pre-Delta era was used to create a prediction model whilst data from the Delta one was used as a validation cohort.

Results

The score development cohort consisted of 1,076 patients in the pre-Delta era, and the validation cohort included 2,856 patients in the Delta one. Four predictors remained after backward and forward stepwise logistic regression: age < 5 years, number of comorbidities, fever, and dyspnea symptoms. The predictive ability of the novel pneumonia score was acceptable with the area under the receiver operating characteristics curve of 0.677 and a well-calibrated goodness-of-fit test (p = 0.098). The positive likelihood ratio for pneumonia was 0.544 (95% confidence interval (CI): 0.491–0.602) in the low-risk category, 1.563 (95% CI: 1.454–1.679) in the moderate, and 4.339 (95% CI: 2.527–7.449) in the high-risk.

Conclusion

This study created an acceptable clinical prediction model which can aid clinicians in performing an appropriate triage for children with COVID-19.

Introduction

Coronavirus disease 2019 (COVID-19) is a pandemic caused by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-COV2) that emerged in 2019. It rapidly spread from Wuhan, China, to all corners of the planet. More than 494 million cases of COVID-19 and 6 million deaths have been reported. In the pediatric context, COVID-19 cases were initially estimated to range from 0.1–15% of total cases. With new variants identified in the second half of 2021, however, the pediatric COVID-19 infection rate grew to 28% of total cases including the adult population. Severe pediatric cases and deaths have been approximately 0.6–5% and 0.3%, respectively [112]. Although infection rate and symptom severity are disproportionately lower in children versus adults, children are still at risk for severe disease and death due to the limitations of vaccination in young children. Children with comorbidities are also at higher risk for the development of a more severe clinical course of infection.

All confirmed COVID-19 cases in Thailand are ostensibly admitted to the healthcare system. Initially, most patients were admitted to the hospital in the pre-Delta era. Due to the increase in patients during the Delta wave, the government added external health care facilities such as field hospitals and home isolation as options for patients with mild disease. A triage system also began for patients: green (asymptomatic or mild symptoms), yellow (mild dyspnea/high-risk group with mild symptoms), red (severe pneumonia with an oxygen saturation < 96% or patients with severe symptoms requiring ventilatory or inotropic support). Green patients were admitted to the field hospitals or home isolation system, whilst yellow and red patients were hospitalized. Bedside physical examination in these patients was limited due to the enhanced personal protective equipment (PPE) and the isolation process which impeded bedside evaluation and chest auscultation. As protocol, all hospitalized patients underwent chest X-rays (CXR), but patients in home isolation underwent CXR only if exhibiting signs of lower respiratory tract infection such as dyspnea or shortness of breath.

A meta-analysis in children by Ding et al. reported that approximately 20% of asymptomatic patients had abnormal chest X-rays (CXR) [8] and another study by Liu et al. revealed that abnormal computed tomography (CT) of the chest compatible with COVID-19 could be found in asymptomatic children [13]. In this regard, this protocol might lead to excessive CXR to rule out lower respiratory tract infections in hospitalized patients but a possible under-detection of pneumonia in home isolation patients. A new predictive score derived from basic patient demographic data and clinical clues could reduce the risk of unnecessary exposure to radiation and also reduce the risk of pneumonia misdiagnosis.

Some data on risk factors for the development of severe diseases in pediatric COVID-19 has come to light. While almost all reports have been from China, Europe, and the United States, it has demonstrated that risk factors for severe diseases include underlying comorbidities, hypoalbuminemia, lymphopenia, and high acute phase reactant [13, 1027]. However, much of this was determined by laboratory investigations that may not be practical for real-life clinical practice in limited-resource areas. Data has also been lacking on pediatric COVID-19 from Southeast and South Asia [4, 28]. In Thailand, there is only a single-center retrospective study looking at the epidemiology and risk factors in pneumonia development and severe pediatric COVID-19 [28]; nonetheless, the information on the pertinent risk factors was still insufficient to create any predictive score. Thus, the main objective of our study was to conduct a nationwide, multi-center study to determine the risk factors associated with the development of pneumonia in children with COVID-19, create a predictive score system, and validate this score.

Materials and methods

Study design

We conducted a retrospective cohort study of all pediatric patients aged from 0–15 years with PCR-confirmed COVID-19 infections between January 13th, 2020, and October 31st, 2021, who were admitted to the healthcare systems which include the hospital, field hospital, and the home isolation system at 13 centers in Thailand: 8 sites in central Thailand, 2 sites in the north, 2 sites in the south, and 1 site in the northeast. Children diagnosed with perinatal COVID-19 infection and those without CXR were excluded from the analysis. Approval for the study was granted by the Central Research Ethics Committee (CREC) of Thailand: CREC0016/2565. An online standardized database was set up using REDCap (Research Electronic Data Capture) with the main coordinating center at Thammasat University. Demographic data and clinical information were collected through a manual chart review and entered through REDCap by experienced pediatricians. Due to the study’s retrospective nature, the informed consent was waived by the CREC.

Model development of pneumonia score

A clinical prediction model for predicting pneumonia in children with COVID-19 was created using patient data in the pre-Delta period and then using this model to externally validate the Delta population. Confirmed COVID-19 cases were classified into two groups based on the Ministry of Public Health of Thailand (MoPH) nationwide surveillance data of SARS-CoV2 variants: a pre-Delta group (diagnosed before July 1st, 2021) and a Delta dominant group (from July 1st, 2021, to October 31st, 2021).

Univariate analyses were used to determine potential risk factors for radiographic pneumonia at a statistical significance level of 5%. Whenever univariate statistical analyses showed the same significance level in two or more parameters that correlated with each other, the variable that provided better discriminative ability was selected. Afterward, backward and forward stepwise logistic regression was performed to acquire a final prediction model. Multivariate analyses were then performed after adjusting for statistically significant confounding factors from the univariate analysis. Clinical relevance was determined by calculating the odds ratios (ORs) with 95% confidence intervals (95% CI) and coefficient. The regression coefficients of each clinical predictor were divided by the smallest coefficient from the multivariate model and rounded to the nearest integer to transform into an item risk score. Scores from each clinical predictor were added up to obtain a total risk score. We called the new scoring system the “pneumonia score”.

Data collection

Clinical information from both the score development cohort and validation one was collected: demographic data, underlying medical conditions, nutritional status, clinical history, initial vital signs, clinical symptoms, CXR findings, disease severity, and clinical outcome.

Operational definition

For children < 5 years of age, obesity was defined as weight-for-height > 3 standard deviations (SD) above the World Health Organization (WHO) Child Growth Standard median [29, 30]. Children and adolescents aged between 5–15 years were defined as obese if BMI-for-age was > 2 SD above the WHO Growth Reference Median [29, 31]. Tachycardia and tachypnea for age were defined using WHO criteria [32]. Desaturation was defined as oxygen saturation < 95%. Pneumonia was defined as patients with clinical signs and symptoms of lower respiratory tract infection such as dyspnea, tachypnea, and desaturation with abnormal CXR or asymptomatic patients with abnormal CXR. All CXR were reviewed by radiologists and/or pediatric pulmonologists to ensure accuracy. Disease severity categorization was based on the National Institutes of Health (NIH) [33]. Five categories were defined: asymptomatic, mild (mild symptoms without pneumonia), moderate (pneumonia without desaturation or no symptoms but abnormal CXR), severe (pneumonia with SpO2 < 94%), or critically ill (patients requiring mechanical ventilation or vasoactive medication; acute respiratory distress syndrome; or septic shock). Immunocompromised children were defined as children who underwent bone marrow or solid organ transplantation, children who are currently taking immunosuppressives medications, children who are undergoing cancer treatment, children with primary immunodeficiencies, and children with advanced or untreated HIV infection [34].

Statistical analyses

Data were analyzed using STATA for Windows v14.0 (StataCorp LLC, Texas, USA). Clinical characteristics in terms of continuous data were reported as mean and standard deviation or median with interquartile range (IQR) and categorical variables were reported as the frequency with percentage. One-way ANOVA, Wilcoxon rank-sum test, and Kruskal-Wallis test were used to compare continuous data. Categorical data analysis was conducted using a Chi-square test adjusted for multiple comparisons.

The external validity of the developed model was evaluated in terms of both discrimination and calibration. The discriminative ability of the model was determined using Mann–Whitney U test. Predictive accuracy was presented as the area under the receiver operating characteristic (AuROC) curve, along with 95% CI. The calibration plot was illustrated, and the Hosmer–Lemeshow test was conducted to demonstrate the results of model calibration. Pneumonia scores were categorized into three risk levels: high, moderate, and low risk. The predictive ability of each risk score level was calculated and presented as a positive likelihood ratio and 95% CI.

Results

Overall patient characteristics

A total of 4,628 confirmed pediatric COVID-19 cases were enrolled during the study period. Approximately 3,932 patients (85%) had CXR and were included in the development and validation analysis. The score development cohort consisted of 1,076 patients in the pre-Delta era, and the validation cohort included 2,856 patients from the Delta era. Weight-for-height was obtained in 2844 cases due to the difficulty in obtaining height with the isolation process at some centers. The proportion of children with obesity did not differ among groups (18.5% VS 20.7%; p = 0.174).

The demographic data of children within all cohorts were described in Table 1. There were significantly fewer patients who were asymptomatic in the validation cohort (14.6%) when compared to the score development cohort (26.6%) (p < 0.001). There was a significantly higher proportion of pneumonia in the validation cohort (32.7% VS 28.8%; p = 0.020). Nevertheless, there were no significant differences in the utilization of respiratory support, intensive care admission, and mortality rates among the cohorts.

thumbnail
Table 1. Demographic data and clinical characteristics in different periods of 3,932 confirmed pediatric COVID-19 cases.

https://doi.org/10.1371/journal.pone.0273842.t001

Model generation

In the development cohort, pneumonia was diagnosed in 310 patients (28.8%); 87 children (28.1%) were asymptomatic. There was a significantly higher proportion of children < 5 years old who developed pneumonia versus older ones (52.9% VS 33.5%; p < 0.001). Children with comorbidities (12.2% VS 7.8%; p = 0.011) and children exhibiting dyspnea symptoms (4.2% VS 0.9%; p < 0.001) were more likely to have pneumonia. There were a significantly higher proportion of children with pneumonia in children who exhibited the symptoms of fever (38.7% VS 25.3%; p < 0.001) and cough (46.8% VS 38.8%; p = 0.016). Children with pneumonia were more likely to be tachycardic for age and have desaturation when compared to children without pneumonia. Gender, obesity status, and immunocompromised status were not associated with pneumonia. Associated risk factors for the development of radiographic pneumonia as determined by univariate analyses are given in Table 2.

thumbnail
Table 2. Univariate analysis of risk factors in the development of radiographic pneumonia for 1,076 confirmed pediatric COVID-19 cases in the pre-Delta period.

https://doi.org/10.1371/journal.pone.0273842.t002

Multivariate analyses were then performed after adjusting for statistically significant confounding factors from the univariate analysis. Four predictors remained after backward and forward stepwise logistic regression: age < 5 years, number of comorbidities, fever, and dyspnea symptoms. Coefficients in terms of the point estimate, odds ratios, and 95% CI were created for the statistically significant clinical variables to generate a score for each parameter. The scores from each parameter were then added to create a sum for the predictive score with a maximum of 16. Details are in Table 3.

Model validation

During the validation period, a total of 933 patients (32.7%) were diagnosed with pneumonia. A total of 113 patients with pneumonia (12.1%) were asymptomatic. The discriminative ability of the pneumonia score, which ranged from 0 to 16, was presented with the ROC curve as shown in S1 Fig. The score predicted abnormal CXR in pediatric patients with an AuROC curve of 0.677 (95% CI: 0.656–0.686; p < 0.001). The pneumonia scores were considered a well-calibrated scoring system, given that the Hosmer–Lemeshow test demonstrated the goodness-of-fit test with the p-value of 0.098. The calibration plot is illustrated in S2 Fig. The pneumonia score was categorized into three risk groups: low (< 3), moderate (3 to 7), and high (> 7) to facilitate clinical interpretation. The positive likelihood ratio (LHR+) for pneumonia was 0.544 (95%CI: 0.491–0.602) in the low-risk category, 1.563 (95% CI: 1.454–1.679) in the moderate group, and 4.339 (95% CI: 2.527–7.449) in the high-risk category (Table 4).

thumbnail
Table 4. Probability categories and the predictive ability.

https://doi.org/10.1371/journal.pone.0273842.t004

Discussion

This study was the first in Thailand to develop and validate a predictive score for pneumonia prediction in children with COVID-19. By using this new predictive score, the aim was to reduce both overuse and underuse of CXR in pneumonia diagnosis. An appropriate triage system is ideal for both patient care and economic/budget considerations, especially where there are limited resources and medical access.

As mentioned, bedside evaluation and physical examination could be extremely limited during patient isolation. Based on the risk score, we were able to categorize children into low risk, intermediate-risk, and high risk for having pneumonia. With the low number of patients with pneumonia in the high-risk group, overuse of CXR was deemed to be as low as 0.9%. Furthermore, a high LHR+ suggested that CXR should be performed in children in high-risk groups. Our study revealed that 20.8% of asymptomatic children during the pre-Delta period and 22.1% during the Delta one had abnormal CXR. These results were comparable with the previous study by Ding et al. which demonstrated pneumonia in 20% of asymptomatic children [8]. In the low-risk group, the results demonstrated that approximately 31.5% of patients had pneumonia. Nevertheless, with the low pneumonia score, most children were either asymptomatic or presented with only mild symptoms, which did not require medical intervention. The decision to perform CXR on low-risk patients might not be cost-effective. Nevertheless, patient monitoring remains essential, and CXR should be performed when the patients exhibit more clinical characteristics which constitute a higher pneumonia score.

To date, only one large, observational cohort study in Brazil was conducted to create a clinical prediction model in children with COVID-19. Despite the difference in the endpoint, which aimed to predict death rather than pneumonia, the study by Oliveira et al. revealed similar results: children with more respiratory symptoms and comorbidities were more likely to have negative outcomes [35]; this coincides with a systematic review by Tsankov et al. [36]. Similar to an Asian multinational study by Wong et al. and Bellino et al. in Italy, our study also found younger children were more likely to have a more severe clinical course [4, 22].

Our study was the first large, nationwide, multicenter study of pediatric COVID-19 from 13 different centers across Thailand, which is a major strength. Unlike other nations, Thailand attempted to admit all patients into the healthcare system, which led to the founding of field hospitals and our home isolation system. Due to this national policy, we were able to obtain large amounts of epidemiological data not only from hospitalized children but also from children in field hospitals and home isolation. With this large and robust clinical dataset, we created a prediction score with a well-fitted calibration model and validated cohort to be used across the country, regardless of region. Different regions of Thailand have certain disparities in terms of access to medical attention. There is also a limited number of pediatricians in rural areas. Therefore, by using our model with its acceptable discrimination power and accuracy in predicting pneumonia, a triage system in these less privileged regions can be set up and provide helpful clinical assistance; this could also aid in public policy decision-making [37].

Nevertheless, several limitations must be outlined. Firstly, due to the retrospective nature of the study and limitation in obtaining the height of children in isolation, weight-for-height was obtained in 2844 cases (72% of the cohort). Despite the missing data, the remaining population was adequately powered for risk factor determination for pneumonia. In these patients, obesity was not associated with pneumonia, unlike for adults, which is in accordance with several other pediatric studies [1, 1012, 14, 1617].

Secondly, due to the low prevalence of severe disease and deaths in children with COVID-19 in Thailand, we were unable to create a prediction model for severe disease. Thus, we used a diagnosis of pneumonia, classified as moderate severity, as a surrogate. Another limitation worth mentioning is that this study was conducted when the pediatric vaccination rollout was still limited; thus, we were unable to include vaccination in the predictive score calculation. Another validation study might be warranted when vaccination coverage in children is higher. Finally, our study was conducted during the pre-Delta and the Delta dominant eras of COVID-19, so we were unable to make any conclusions on the effectiveness of its discrimination power on the future variants of COVID-19.

Conclusion

This study created an acceptable clinical prediction model to predict pneumonia for children with COVID-19. Thus, by using this clinical prediction score, children can be triaged more appropriately and received improved medical care.

Supporting information

S1 Fig. Receiver-operating characteristic curve of the risk score to predict pneumonia.

https://doi.org/10.1371/journal.pone.0273842.s001

(TIF)

S2 Fig. Calibration plot of pneumonia score.

Pr (obs_pneumonia) = predicted risk; obsrisk = observed risk.

https://doi.org/10.1371/journal.pone.0273842.s002

(TIF)

Acknowledgments

The authors are grateful to The Royal College of Pediatricians of Thailand and Pediatric Society of Thailand, The Thai Society of Pediatric Respiratory and Critical Care Medicine, and the Research Group of Pediatrics, from the Faculty of Medicine, Thammasat University for the great support and advice for making this project possible. We are also grateful to the pediatric department staff and healthcare workers responsible for patient care at all centers. We would also like to show great appreciation to Ms. Debra Kim Liwiski for her great comprehensive English review and correction.

References

  1. 1. Dong Y, Mo X, Hu Y, Qi X, Jiang F, Jiang Z, et al. Epidemiology of COVID-19 among children in China. Pediatrics. 2020;145: e20200702. pmid:32179660
  2. 2. COVID-19 CDC Response Team, Bialek S, Gierke R, Hughes M, McNamara LA, Pilishvili T, Skoff T. Coronavirus Disease 2019 in children—United States, February 12-April 2, 2020. MMWR Morb Mortal Wkly Rep. 2020;69: 422–6. pmid:32271728
  3. 3. Guo CX, He L, Yin JY, Meng XG, Tan W, Yang GP, et al. Epidemiological and clinical features of pediatric COVID-19. BMC Med. 2020;18: 250. pmid:32762696
  4. 4. Wong JJM, Abbas Q, Chuah SL, Malisie RF, Pon KM, Katsuta T, et al. Comparative analysis of pediatric COVID-19 infection in Southeast Asia, South Asia, Japan, and China. Am J Trop Med Hyg. 2021;105: 413–20. pmid:34129517
  5. 5. Hoang A, Chorath K, Moreira A, Evans M, Burmeister-Morton F, Burmeister F, et al. COVID-19 in 7780 pediatric patients: a systematic review. EClinicalMedicine. 2020;24: 100433. pmid:32766542
  6. 6. Badal S, Bajgain KP, Badal S, Thapa R, Bajgain BB, Santana MJ. Prevalence, clinical characteristics, and outcomes of pediatric COVID-19: A systematic review and meta-analysis. J Clin Virol. 2021;135: 104715. pmid:33348220
  7. 7. Irfan O, Muttalib F, Tang K, Jiang L, Lassi ZS, Bhutta Z. Clinical characteristics, treatment and outcomes of paediatric COVID-19: a systematic review and meta-analysis. Arch Dis Child. 2021;106: 440–8.
  8. 8. Ding Y, Yan H, Guo W. Clinical characteristics of children with COVID-19: A meta-analysis. Front Pediatr. 2020;8: 431. pmid:32719759
  9. 9. Wang Z, Zhou Q, Wang C, Shi Q, Lu S, Ma Y, et al. COVID-19 Evidence and recommendations working group. Clinical characteristics of children with COVID-19: a rapid review and meta-analysis. Ann Transl Med. 2020;8: 620.
  10. 10. Lu X, Zhang L, Du H, Zhang J, Li YY, Qu J, et al. SARS-CoV-2 infection in children. N Engl J Med. 2020;382: 1663–5. pmid:32187458
  11. 11. DeBiasi RL, Song X, Delaney M. Bell M, Smith K, Pershad J, et al. Severe Coronavirus Disease-2019 in children and young adults in the Washington, DC, Metropolitan Region. J Pediatr. 2020;223: 199–203e1. pmid:32405091
  12. 12. Parri N, Lenge M. Children with Covid-19 in pediatric emergency departments in Italy. N Engl J Med. 2020;383: 187–90. pmid:32356945
  13. 13. Liu M, Song Z, Xiao K. High resolution computed tomography manifestations of 5 pediatric patients with 2019 novel Coronavirus. J Comp Assist Tomogr. 2020;44: 311–3. pmid:32217900
  14. 14. Chao JY, Derespina KR, Herold BC, Goldman DL, Aldrich M, Weingarten J, et al. Clinical characteristics and outcomes of hospitalized and critically ill children and adolescents with Coronavirus Disease 2019 at a tertiary care medical center in New York City. J Pediatr. 2020;223: 14–9. e2. pmid:32407719
  15. 15. Whittaker E, Bamford A, Kenny J, Kaforou M, Jones CE, Shah P, et al. Clinical characteristics of 58 children with a pediatric inflammatory multisystem syndrome temporally associated with SARS-CoV-2. JAMA. 2020;324: 259–69. pmid:32511692
  16. 16. Shekerdemian LS, Mahmood NR, Wolfe KK, Riggs BJ, Ross CE, McKiernan CA, et al. Characteristics and outcomes of children with Coronavirus Disease 2019 (COVID-19) infection admitted to US and Canadian pediatric intensive care units. JAMA Pediatr. 2020;174: 868–73. pmid:32392288
  17. 17. Tagarro A, Epalza C, Santos M, Sanz-Santaeufemia FJ, Otheo E, Moraleda C, et al. Screening and severity of coronavirus disease 2019 (COVID-19) in children in Madrid, Spain. JAMA Pediatr. 2020: e201346. pmid:32267485
  18. 18. Qiu H, Wu J, Hong L, Luo Y, Song Q, Chen D, et al. Clinical and epidemiological features of 36 children with coronavirus disease 2019 (COVID-19) in Zhejiang, China: an observational cohort study. Lancet Infect Dis. 2020;20: 689–96. pmid:32220650
  19. 19. Belhadjer Z, Meot M, Bajolle F, Khraiche D, Legendre A, Abakka S, et al. Acute heart failure in multisystem inflammatory syndrome in children (MIS-C) in the context of global SARS-CoV-2 pandemic. Circulation. 2020;142: 429–36.
  20. 20. Bandi S, Nevid MZ, Mahdavinia M. African American children are at higher risk for COVID-19 infection. Pediatr Allergy Immunol. 2020;31: 861–4.
  21. 21. Zheng F, Liao C, Fan QH, Chen HB, Zhao XG, Xie ZG, et al. Clinical characteristics of children with Coronavirus Disease 2019 in Hubei, China. Curr Med Sci. 2020;40: 275–80. pmid:32207032
  22. 22. Bellino S, Punzo O, Rota MC, Manso MD, Urdiales AM, Andrianou X, et al. COVID-19 disease severity risk factors for pediatric patients in Italy. Pediatrics. 2020;146: e2020009399. pmid:32665373
  23. 23. Zhou B, Yuan Y, Wang S, Zhang Z, Yang M, Deng X, et al. Risk profiles of severe illness in children with COVID-19: a meta-analysis of individual patients. Pediatr Res. 2021;90: 347–52. pmid:33753892
  24. 24. Graff K, Smith C, Silveira L, Jung S, Curran-Hays S, Jarjour J, et al. Risk factors for severe COVID-19 in children. Pediatr Infect Dis J. 2021;40: e137–e145. pmid:33538539
  25. 25. Moreno-Noguez M, Rivas-Ruiz R, Roy-García IA, Pacheco-Rosas DO, Moreno-Espinosa S, Flores-Pulido AA. Risk factors associated with SARS-CoV-2 pneumonia in the pediatric population. Bol Med Hosp Infant Mex. 2021;78: 251–8.
  26. 26. Tsabouri S, Makis A, Kosmeri C, Ekaterini Siomou E. Risk factors for severity in children with Coronavirus Disease 2019: A comprehensive literature review. Pediatr Clin North Am. 2021;68: 321–38. pmid:33228941
  27. 27. Mania A, Faltin K, Mazur-Melewska K, Małecki P, Jończyk-Potoczna K, Lubarski K, et al. Clinical picture and risk factors of severe respiratory symptoms in COVID-19 in children. Viruses. 2021;13: 2366. pmid:34960635
  28. 28. Anugulruengkitt S, Teeraananchai S, Chantasrisawad N, Promsena P, Jantarabenjakul W, Puthanakit T. Clinical outcomes of pediatric COVID-19 in a tertiary care center in Bangkok, Thailand. IJID Regions. 2021;1: 159–62. pmid:35721777
  29. 29. World Health Organization. Obesity and overweight. 2021. Available from: https://www.who.int/news-room/fact-sheets/detail/obesity-and-overweight.
  30. 30. World Health Organization. World Health Organization Child Growth Standards. 2006. Available from: https://www.who.int/tools/child-growth-standards.
  31. 31. de Onis M, Onyango AW, Borghi E, Siyam A, Nishida C, Siekmann J. Development of a WHO growth reference for school-aged children and adolescents. Bull World Health Organ. 2007;85: 660–7 pmid:18026621
  32. 32. World Health Organization. The management of acute respiratory infection in children In: Practical guideline for outpatient care. Geneva, WHO 1995
  33. 33. Coronavirus Disease 2019 (COVID-19) Treatment Guidelines. Special consideration in children. 2022. Available from: https://www.covid19treatmentguidelines.nih.gov/special-populations/children/.
  34. 34. National Center for Immunization and Respiratory Diseases (NCIRD), Division of Viral Diseases. COVID-19 Vaccines for People who are Moderately or Severely Immunocompromised. 2022. Available from: https://www.cdc.gov/coronavirus/2019-ncov/vaccines/recommendations/immuno.html.
  35. 35. Oliveira EA, Simoes e Silva A, Oliveira MCL, Colosimo EA, Mak RH, Vasconcelos MA, et al. Comparison of the first and second waves of the Coronavirus Disease 2019 (COVID-19) pandemic in children and adolescents in a middle-income country: Clinical impact associated with Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Gamma lineage. J Pediatr. 2022; S0022-3476(22)00001-4. pmid:35031347
  36. 36. Tsankov BK, Allaire JM, Irvine MA, Lopez AA, Sauve LJ, Vallance BA et al. Severe COVID-19 infection and pediatric comorbidities: A systematic review and meta-analysis. Int J Infect Dis 2021;103: 246–256. pmid:33227520
  37. 37. Wynants L, Sotgiu G. Improving clinical management of COVID-19: the role of prediction model. Lancet Respir Med 2021;9; 320–1.