Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Pediatr., 20 November 2025

Sec. Pediatric Pulmonology

Volume 13 - 2025 | https://doi.org/10.3389/fped.2025.1695974

Development and validation of a nomogram for early prediction of macrolide-unresponsive Mycoplasma pneumoniae pneumonia in children


Hui Liu,,&#x;Hui Liu1,2,†Rina Sa,&#x;Rina Sa3,†Yang XinYang Xin4Meilan ZhangMeilan Zhang2Hai Du

Hai Du2*
  • 1Graduate School, Baotou Medical College, Baotou, China
  • 2Department of Radiology, Ordos Central Hospital, Ordos, China
  • 3Department of Pediatrics, Ordos Central Hospital, Ordos, China
  • 4Ordos Clinical Medical College, Inner Mongolia Medical University, Ordos, China

Background: Mycoplasma pneumoniae (MP) is a leading cause of community-acquired pneumonia in children, with a significant increase in incidence following the COVID-19 pandemic. The emergence of macrolide-resistant M. pneumoniae (MRMP) has complicated treatment, leading to the concept of macrolide-unresponsive M. pneumoniae pneumonia (MUMPP), defined as lack of improvement after 72 h of macrolide therapy. Early identification of MUMPP is critical for timely intervention and improved outcomes. This study aimed to develop and validate a nomogram for early prediction of MUMPP in children.

Methods: We conducted a retrospective study involving 278 pediatric patients with MP pneumonia, divided into training (n = 188) and validation (n = 90) sets. Demographic, clinical, laboratory, and chest CT imaging data were collected. Univariate and multivariate logistic regression analyses were used to identify independent predictors of MUMPP. A nomogram was constructed and validated using receiver operating characteristic (ROC) curves, calibration plots, and decision curve analysis (DCA).

Results: Six independent predictors were identified: tree-in-bud pattern, neutrophil-value, lymphocyte-value, creatine kinase (CK), platelet-to-lymphocyte ratio (PLR), and male gender. The nomogram demonstrated strong discriminatory power, with area under the curve (AUC) values of 0.838 (95% CI: 0.779–0.897) in the training set and 0.835 (95% CI: 0.752–0.918) in the validation set. Calibration and DCA confirmed good clinical utility.

Conclusion: We developed and validated a simple-to-use nomogram for predicting MUMPP in early stage. The nomogram demonstrates strong discriminatory power and calibration, and may be a practical tool for clinical practice.

Introduction

Mycoplasma pneumoniae pneumonia (MPP) accounts for approximately 10%–40% of community-acquired pneumonia (CAP) in school-aged children (1). According to reports from the World Health Organization and multiple countries, the burden of respiratory diseases in children has increased significantly following the COVID-19 pandemic, with a notable rise in the incidence of MPP (2). Macrolides are the primary antibiotics used to treat MPP. However, compared to the pre-pandemic period, macrolide-resistant Mycoplasma pneumoniae (MRMP) has become increasingly prevalent worldwide, particularly in Asian countries. Resistance rates have reached approximately 69.67% in South Korea and as high as 90% in Japan and China (35). This trend poses a serious threat to children's health, often leading to prolonged fever, extended hospital stays, increased severity of cases, and challenges in antibiotic selection (6). In response, the 2023 Chinese Guidelines for the Diagnosis and Treatment of MPP in Children introduced the concept of “Macrolide-Unresponsive Mycoplasma pneumoniae Pneumonia (MUMPP), defined as cases showing no clinical or radiological improvement after 72 h of macrolide treatment” (7). MUMPP can serve as an early indicator of MRMP. Therefore, early identification of MUMPP and timely adjustment of antibiotics are crucial for shortening the disease course and reducing the risk of severe complications and sequelae (8).

Imaging plays an increasingly important role in the diagnosis and evaluation of pneumonia. Although chest x-ray remains the first-line imaging method for assessing CAP, its findings in MPP are often non-specific and may be missed due to anatomical overlap. In contrast, computed tomography (CT) provides clearer visualization of parenchymal and interstitial lung abnormalities—such as bronchial wall thickening and bronchiectasis—which correlate well with pathological changes (9, 10). With higher spatial and density resolution, CT offers a more reliable basis for accurately assessing the extent and pattern of lung involvement in children, thereby guiding clinical decision-making.

The nomogram is an intuitive and user-friendly tool for multivariate prediction modeling. It integrates key predictive factors into a visual risk assessment model, facilitating individualized risk evaluation and supporting clinical decision-making (11). Existing prediction models for refractory MPP (RMPP), severe MPP (SMPP), and necrotizing MPP (NMPP) have identified several independent risk factors, including peak body temperature, pleural effusion, extrapulmonary complications, neutrophil ratio, C-reactive protein (CRP), erythrocyte sedimentation rate (ESR), procalcitonin (PCT), D-dimer, lactate dehydrogenase (LDH), and albumin (ALB) (1214). However, few studies have focused on developing clinical prediction models specifically for MUMPP. Current diagnosis relies heavily on symptomatic presentation and isolated laboratory findings, which vary based on physicians' experience and subjective judgment, leading to inconsistent and potentially biased outcomes. Hence, there is an urgent need to develop an effective, economical, and rapid diagnostic tool to support timely and accurate treatment.

This study aimed to identify factors for the early prediction of MUMPP and to develop a simple-to-use nomogram.

Materials and methods

Study subjects and design

This retrospective research was conducted in adherence to the Declaration of Helsinki. Ethical approval, which included a waiver for individual consent, was obtained from the Ethics Committee of Ordos Central Hospital (No. 2025-406).

In this study, 278 children diagnosed with MPP and hospitalized for treatment were enrolled. Among them, 188 cases from the main campus of Ordos Central Hospital between July 2018 and June 2024 were assigned to the training set. An additional 90 cases from the branch campus between April 2019 and June 2024 were collected as the test set, (Figure 1).

Figure 1
Flowchart illustrating the selection process for hospitalized children with MPP from 2018.7 to 2024.6. Out of 3685 children, 278 were included after exclusions for other pathogen infections, prior macrolide treatment, specific CT examination criteria, and incomplete data. The final cohort was divided into a training cohort of 188 and a test cohort of 90, each further separated into MSMMP and MUMPP groups.

Figure 1. Patient screening workflow for this study. MUMPP, macrolide-unresponsive M. pneumonias pneumonia, MSMPP macrolide-sensitive M. pneumonias pneumonia.

Inclusion criteria: (1) were under the age of 14 years; (2) clinical signs and symptoms suggestive of CAP, including fever, cough, adventitious breath sounds (rales or diminished breath sounds), and new infiltrates on chest x-ray or CT; (3) positive serological testing for MP-IgM (serum anti-MP IgM titer ≥1:160 or a ≥ 4-fold increase in antibody levels), or positive MP-RNA or MP-DNA detection in throat swab or bronchoalveolar lavage fluid (BALF); (4) noncontrast chest CT performed within 72 h of admission.

The following tests were used to rule out other respiratory infections and tuberculosis: purified protein derivative (PPD) test, blood culture, pleural fluid culture, nasopharyngeal aspirate/swab culture, viral antigen testing (for respiratory syncytial virus, influenza virus, adenovirus, and parainfluenza virus), and serology for Chlamydia pneumoniae and Legionella pneumophila.

Exclusion criteria included: (1) macrolide antibiotic use for ≥72 h prior to admission; (2) foreign body aspiration; (3) chronic respiratory diseases; (4) malignancy, solid organ transplantation or surgery, immunodeficiency, or use of immunosuppressive drugs; (5) congenital or inherited metabolic disorders and other underlying conditions; (6) incomplete medical records; (7) poor patient cooperation leading to significant imaging artifacts. Furthermore, patients with any missing values in the studied variables were excluded from the final analysis to ensure a complete dataset for model development.

Depending on whether patients had persistent fever, unresolved clinical symptoms, and worsening or stagnant pulmonary imaging following 72 h of macrolide therapy, they were classified into Macrolide-Sensitive Mycoplasma pneumoniae Pneumonia (MSMPP) group and unresponsive (MUMPP) group.

Data collection

We collected data from the electronic health medical record system, including general demographics, clinical features, laboratory results at admission.

(1) Demographic data included gender, age, and ethnicity.

(2) Clinical features encompassed presence of fever, cough, wheezing, pre-admission cough duration (PCD), pre-admission fever duration (PFD), the highest body temperature (T), heart rate (HR) and respiratory rate (R) on admission and length of hospital stay.

(3) Laboratory results consisted of white blood cell count (WBC), neutrophil count (NEUT_value), neutrophil percentage (NEUT), lymphocyte count (LYMPH_value), lymphocyte percentage (LYMPH), monocyte count (MONO), platelet count (PLT), C-reactive protein (CRP), alanine aminotransferase (ALT), aspartate aminotransferase (AST), lactate dehydrogenase (LDH), creatine kinase (CK), creatine kinase isoenzyme (CK-MB), and albumin (ALB). Additionally, ratios of certain serum markers were evaluated, including ratio of neutrophil-value to lymphocyte-value (NLR), ratio of platelet to lymphocyte-value (PLR), ratio of lymphocyte-value to CRP (LCR), and LDH-to-albumin ratio (LAR).

All patients underwent noncontrast chest CT within 72 h of admission.

CT examinations

CT scans were performed while the children were in a calm state. For uncooperative patients, appropriate sedation was administered to ensure successful completion of the examination. Scans were conducted using GE Lightspeed VCT 64-slice, GE Discovery CT 750HD, GE Revolution CT 256-slice, and SIEMENS SOMATOM Force dual-source CT scanners. Scanning parameters were as follows: tube voltage 80–120 kV, tube current 50–350 mA, slice thickness 5 mm, slice interval 5 mm, and reconstructed image thickness 1.25 mm. The CT scan was performed from the lung apex to the lung base. For mediastinal window settings, the window width was 350 HU and window level 40 HU; for lung windows, the window width was 1,500 HU and window level −600 HU.

Imaging analyses

CT images were initially evaluated by a radiologist with 10 years of experience in thoracic imaging and reviewed by another with more than 15 years of experience. In case of disagreement, a consensus was reached through discussion or consultation with a senior attending radiologist. Recorded CT features included lesion location, number of involved lobes, unilateral or bilateral involvement, bronchial wall thickening, tree-in-bud sign, air bronchogram, interlobular septal thickening, consolidation pattern (patchy, segmental, lobar), consolidation with ground-glass opacity (GGO), mosaic pattern, with or without pleural effusion.

Statistical analysis

All statistical analyses in this study were performed using R software (version 4.5.1). Continuous variables following a normal distribution were presented as mean ± standard deviation and compared using the Student's t-test. Non-normally distributed continuous variables were summarized as median with interquartile range [M (P25, P75)] and compared using the Mann–Whitney U test. Categorical variables were expressed as frequency (percentage) and compared using the chi-square test.

For variable selection, potential predictors were first screened by univariate logistic regression, retaining those with a significance level of P < 0.05. Multicollinearity among these selected variables was then assessed using the variance inflation factor (VIF), and variables with VIF > 5 were iteratively removed. The remaining variables with low collinearity were entered into a multivariate logistic regression model. A stepwise backward elimination method was applied, with inclusion and exclusion criteria set at P < 0.05 and P > 0.10, respectively, to identify the final independent predictors and construct the prediction model. In addition to VIF assessment, correlation analysis was conducted to further verify the absence of severe multicollinearity among the final predictors in the model. Furthermore, for continuous variables included in the final model, four-knot restricted cubic spline analysis was employed to assess the presence of nonlinear relationships with the logit-transformed outcome variable.

The predictive performance of the model was evaluated using the area under the receiver operating characteristic (ROC) curve, accuracy, sensitivity, and specificity. Predictive consistency was assessed using calibration curves, and clinical utility was determined via decision curve analysis. All statistical tests were two-sided, P < 0.05 was considered statistically significant.

Results

Clinical characteristics

Among the 278 pediatric cases, there were 111 cases of MUMPP and 167 cases of MSMPP. Including 125 boys (44.96%) and 153 girls (55.04%). The median age was 7.23 ± 2.64 years, and the average length of hospital stay was 8.20 ± 3.10 days. Mycoplasma pneumoniae pneumonia cases were detected throughout the year, with a peak incidence observed in the fourth quarter. Table 1 summarizes the characteristics of children in the training and test sets at admission. No significant differences were observed in gender or ethnicity (P > 0.05). Clinical manifestations, routine blood tests, and biochemical markers within 24 h of admission were compared between the two groups. In the training set, compared with the MSMPP group, the MUMPP group had a longer duration of fever before admission, higher levels of NEUT, NLR, CRP, AST, CK, LDH, PLR, and LAR, as well as lower respiratory rate, LYMPH_value, LYMPH, MONO, ALB, and LCR. In the test set, the MUMPP group exhibited a longer duration of cough before admission, and significant differences were found in NEUT, LYMPH, LYMPH_value, NLR, CK, and PLR (P < 0.05). Notably, in both the training and test sets, children in the MUMPP group were older and had a higher incidence in the fourth quarter. The other results showed no significant differences (P > 0.05).

Table 1
www.frontiersin.org

Table 1. Baseline characteristics of patients in the training and test cohort.

Chest CT imaging features

Figure 6A demonstrates the characteristic tree-in-bud sign on a chest CT image from a child with MUMPP. Chest CT imaging features are summarized in Table 2. Compared with the MSMPP group, the MUMPP group showed a significantly higher prevalence of tree-in-bud signs (P < 0.01). In the training set, mosaic pattern and lobar consolidation were also more common in the MUMPP group (P < 0.05), though these differences were not statistically significant in the test set. Other imaging features, including bronchial wall thickening, air bronchogram, interlobular septal thickening, segmental consolidation, patchy consolidation, bilateral lobar infiltration, and pleural effusion, did not differ significantly between the groups (P > 0.05).

Table 2
www.frontiersin.org

Table 2. Radiological features of patients in the training and test cohort.

Variable selection

Univariate logistic regression analysis was performed on all features. Given that multiple indicators were associated with inflammatory response, 19 independent variables were first assessed for multicollinearity. Following the removal of highly correlated variables (LAR, LYMPH, and NLR), 16 low-collinearity features were retained for multivariate logistic regression. The analysis identified six independent predictors of MUMPP: sex, NEUT_value, LYMPH_value, CK, PLR, and tree-in-bud sign (P < 0.05) (Table 3, which only lists variables with statistical significance in the multivariate analysis).

Table 3
www.frontiersin.org

Table 3. Results of multivariable logistic regression analysis.

A comprehensive correlation analysis confirmed the absence of substantial multicollinearity among the six final predictors (Figure 2). Evaluation of the continuous predictors using a four-knot restricted cubic spline (RCS) analysis revealed no significant nonlinear relationships with the logit of the outcome (P > 0.05) (Supplementary Tables S1, S2). Consequently, these variables were retained in their original linear form in the final model.

Figure 2
Two correlation heatmaps display significant predictors for training and test sets. Both heatmaps show correlations between variables: Sex, NEUT_value, LYMPH_value, CK, PLR, and tree_in_bud_sign. Color ranges from red (positive correlation) to blue (negative correlation) with values shown in each cell.

Figure 2. Heat map analysis of correlations between 6 variables. (A) Training cohort. (B) Test cohort. NEUT-value, neutrophil count; LYMPH-value, lymphocyte count; CK, creatine kinase; PLR, platelet-to-lymphocyte ratio.

Development and evaluation of the predictive nomogram

A nomogram was constructed for the early identification of MUMPP in children containing six independent predictors (Figure 3). A higher total score on the nomogram indicates a greater risk of MUMPP. In the training set, the nomogram showed an area under the curve (AUC) of 0.838 (95% CI: 0.779–0.897) (Figure 4A), while in the test set, the AUC was 0.835 (95% CI: 0.752–0.918) (Figure 4B), demonstrating robust discriminatory ability. Furthermore, the model exhibited robust specificity and PPV in the training cohort, a marked gain in sensitivity was achieved in the test cohort while maintaining strong specificity (Table 4). To rigorously evaluate the internal validity of the model, bootstrap validation with 1,000 resamples was performed, and the corresponding performance metrics are detailed in Supplementary Table S3.

Figure 3
A nomogram with multiple scales labeled for variables: points, sex, NEUT_value, LYMPH_value, CK, PLR, and tree_in_bud_sign. The scales correlate with the total points, linear predictor, and MUMPP probability. Each variable has a corresponding point scale, contributing to overall risk assessment.

Figure 3. The nomogram to predict MUMPP. *Sex: 1 means female; 0 means male; tree_in_bud_sign: lmeans yes; 0 means no.

Figure 4
Graphs showing ROC and calibration curves for model performance. Panel A: Training set ROC, AUC 0.838. Panel B: Test set ROC, AUC 0.835. Panel C: Training set calibration curve, Hosmer-Lemeshow p=0.811. Panel D: Test set calibration curve, Hosmer-Lemeshow p=0.907.

Figure 4. Receiver operating characteristic curves (ROC) of the model in the training set (A) and test set (B); AUC, area under the curve. Calibration curves of the model in the training set (C) and test set (D).

Table 4
www.frontiersin.org

Table 4. Diagnostic performance of model for MUMPP in training and test cohort.

Calibration curves indicated good agreement between the predicted and actual outcomes (Figures 4C,D). The DCA of our nomogram shows that across a wide threshold probability range of approximately 0.1–0.9, its net benefit is consistently higher than the strategies of “intervening on all patients” and “intervening on no patients.” This indicates that in clinical practice, using this model to screen patients requiring early identification of MUMPP can provide greater net clinical benefit compared to simple binary strategies, effectively balancing the risks of missed diagnosis and overdiagnosis (Figure 5).

Figure 5
Two line graphs comparing decision curves. Panel A displays the training data, and Panel B the test data. Both graphs plot net benefit against threshold probability. The red line represents the nomogram, suggesting a higher net benefit than the gray lines, which represent treating all or none.

Figure 5. Decision curve analysis (DCA) of the CT-based nomogram (red line). All (gray solid line, this assumes you intervene on every single patient) and None (thick black line at y = 0, this assumes you do not intervene on any patient). The x-axis indicates the threshold probability. The y-axis indicates the net benefit. DCA comparing the net benefit of the nomogram in the training set (A) and test set (B).

Figure 6 provides an example of how to use the nomogram. Similar to a scoring system, points are assigned for each predictor of MUMPP, which correspond to the risk of MUMPP. A vertical line can be drawn upward from each predictor to determine the points associated with the presence of the tree-in-bud sign on chest CT, gender, neutrophil count, lymphocyte count, creatine kinase, and platelet-to-lymphocyte ratio (PLR). Once points are assigned for all predictors, the total points are calculated. The total points are then converted into the probability of MUMPP by reading the corresponding value on the total points scale.

Figure 6
CT scan showing two images of lungs. Image A includes a red arrow pointing to a specific area within the left lung. Image B displays both lungs without annotations, indicating different internal structures.

Figure 6. The risk of developing MUMPP was calculated for two pediatric patients using a nomogram. (A) A chest CT image from a boy with MPP shows the tree-in-bud sign (red arrow). His laboratory findings were as follows: NEUT-value, 10.22 × 10⁹/L; LYMPH-value, 1.98 × 10⁹/L; CK, 257 U/L; PLR, 267. The total nomogram score was approximately 158 points, corresponding to a probability of developing MUMPP of greater than 90%. (B) A chest CT image from a girl with MPP shows no tree-in-bud sign. Her laboratory findings were as follows: NEUT-value, 6.42 × 10⁹/L; LYMPH-value, 2.76 × 10⁹/L; CK, 197 U/L; PLR, 136. The total nomogram score was approximately 111 points, corresponding to a probability of developing MUMPP of less than 20%. NEUT-value, neutrophil count; LYMPH-value, lymphocyte count; CK, creatine kinase; PLR, platelet-to-lymphocyte ratio.

Consequently, we created an online computing platform (https://dynom.shinyapps.io/DynNomapp/). This platform enables both doctors and patients to carry out calculations directly on the web.

Discussion

In recent years, the prevalence of MPP has increased significantly, with high rates of complications and mortality (15). Macrolides are the primary antibiotics used to treat MPP, however, due to the increased use of macrolides and the emergence of drug-resistant strains, the number of macrolide-resistant cases has been rising. This situation has become particularly complex following the COVID-19 pandemic, further complicating the clinical management of MUMPP in children (16, 17).

We found that a nomogram constructed using logistic regression, incorporating tree-in-bud pattern, lymphocyte-value, neutrophil-value, creatine kinase (CK), platelet-to-lymphocyte ratio (PLR), and sex, can effectively and conveniently predict MUMPP at an early stage. The nomogram demonstrated strong accuracy and discriminative ability, which indicated it may be a practical tool to help pediatricians recognize MUMPP earlier.

In our study, the tree-in-bud sign was identified as a significant independent predictor of MUMPP. The tree-in-bud sign, a CT pattern characterized by centrilobular nodules and branching linear opacities resembling a budding tree, reflects pathological changes in the small airways, including bronchiolar wall thickening, luminal impaction with mucus, pus, or granulation tissue, and peribronchiolar inflammation (18, 19). The pathogenesis of MUMPP involves not only direct microbial damage but also a robust host immune-inflammatory response, leading to a massive release of inflammatory cytokines and subsequent tissue injury (20, 21). A more pronounced immune response is associated with greater organ damage and dictates disease prognosis. In this context, a widespread tree-in-bud sign signifies extensive lung involvement and severe small airway obstruction (22, 23), which can impede pathogen clearance, reduce local antibiotic concentration, and ultimately diminish treatment responsiveness, as supported by previous findings (24). By incorporating this key imaging feature, our model facilitates the early recognition of MUMPP, thereby aiding in timely treatment adjustment, mitigating the risk of complications such as bronchiolitis obliterans, and improving long-term outcomes in children.

Peripheral blood markers provide important insights into the host's immune-inflammatory status. This study revealed significantly elevated neutrophil-value in children with MUMPP, identifying it as an independent risk factor. This aligns with previous studies suggesting that neutrophil-mediated hyperinflammatory responses are closely associated with MPP severity and poor outcomes (2527). A possible explanation is that neutrophils, as the first line of defense against infection, are activated upon MP invasion via granulocyte colony-stimulating factor (G-CSF) produced in large quantities by bronchial epithelial cells, leading to a rapid increase in neutrophil counts and enhanced phagocytic activity (28). Persistent inflammation may result in lymphocyte-value depletion and impaired cellular immune function, which is supported by the decreased lymphocyte-value observed in the MUMPP group. Platelets are not only involved in coagulation but also act as important inflammatory effector cells. MP infection, along with associated cytokine storms and hypoxia, can cause endothelial injury, leading to platelet activation and consumption (29). Previous studies have shown that lower PLT levels are associated with a higher risk of SMPP (30). Although PLT alone did not differ significantly between groups in this study, a lower PLR was identified as a risk factor for MUMPP. As a composite indicator, PLR may more sensitively reflect concurrent platelet consumption and lymphopenia, indicating stronger inflammatory responses and potential immune imbalance, which correlate with disease severity (31).

Additionally, a higher proportion of male patients was observed in the MUMPP group. Although the relationship between sex and MUMPP remains underexplored, this phenomenon could potentially be due to variations in immune function or hormonal regulation between males and females (32). Elevated serum CK was also more common in children with MUMPP. CK is a key enzyme in energy metabolism within muscle cells, predominantly found in skeletal and cardiac muscle, with smaller amounts present in the brain, intestine, liver, spleen, and lungs. Under normal conditions, CK rarely leaks out of cells due to intact cell membranes. Higher CK levels in MUMPP may reflect extrapulmonary involvement (e.g., muscle damage) or a stronger systemic stress response (33).

Furthermore, The COVID-19 pandemic and its subsequent impacts have reshaped the epidemiological landscape and clinical management of MPP. Global studies indicate that public health interventions during the pandemic and shifts in population immunity created favorable conditions for the resurgence and spread of MPP, often accompanied by increased macrolide resistance (34). In terms of technological responses, breakthroughs in COVID-19 vaccine development—particularly highly stable and immunogenic circular RNA platforms—offer potential avenues for future strategies against drug-resistant Mycoplasma strains (35). Meanwhile, national surveillance data on fever of unknown origin (FUO) in China reveal that infectious diseases remain the leading cause of FUO in the post-pandemic era. However, the proportion of undiagnosed cases has risen, reflecting growing diagnostic complexity. Significant regional variations in etiology further underscore the need for individualized diagnosis tailored to local epidemiological patterns (36). Against this backdrop, our study addresses the diagnostic challenges of the post-COVID era by developing an objective prediction tool—the nomogram model presented here—for early identification of MUMPP. This tool aims to support more precise early clinical intervention in complex diagnostic scenarios.

Although the nomogram model established in this study has demonstrated favorable predictive performance, medical research is advancing toward an era driven by artificial intelligence (AI) (37). Future models may leverage deep learning to enable automated analysis of CT images and integrate multimodal data—including clinical features, laboratory indicators (such as the predictors identified in this study), genomic information, and radiomic features—to construct more powerful predictive models (38, 39). Such integrated models would not only predict MUMPP risk as achieved in the current study but also hold the potential for dynamic and individualized assessment of treatment response. By continuously monitoring changes in imaging and clinical data throughout the treatment process, these models could provide data-driven support for real-time adjustments to therapeutic strategies, thereby paving the way for personalized management of MPP.

There are also some limitations in this study. First, this study is retrospective in design, and all data were derived from a single medical center. Although an internal test set was used for validation, its relatively small sample size may lead to an overestimation of model performance. Consequently, the generalizability of our model may be limited when applied to other populations or healthcare institutions with different diagnostic and treatment standards. Second, regarding predictor selection, our model was developed primarily based on routine clinical and imaging indicators. It was thus developed without incorporating potentially influential microbiological factors, such as specific macrolide resistance-associated gene mutations, and does not account for the potential associations with a history of prior COVID-19 infection or vaccination status. The latter is particularly relevant as the post-pandemic immune landscape may influence the presentation and severity of other respiratory infections like MPP. Furthermore, the prediction model is constructed using traditional logistic regression methodology and has not yet incorporated more advanced AI modeling techniques, which might better capture complex, non-linear relationships. Finally, the clinical applicability and robustness of this nomogram require further validation through prospective, multi-center studies. The retrospective, single-center design may introduce selection bias, and the observed differences in baseline characteristics between the training and test cohorts highlight the need for validation in larger, more diverse populations to ensure generalizability.

Conclusion

We developed and validated a simple-to-use nomogram for predicting MUMPP in early stage. The nomogram demonstrates strong discriminatory power and calibration, and may be a practical tool for clinical practice.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by Ethics Committee of Ordos Central Hospital (Approval No. 2025-406). The studies were conducted in accordance with the local legislation and institutional requirements. The ethics committee/institutional review board waived the requirement of written informed consent for participation from the participants or the participants' legal guardians/next of kin due to the retrospective nature of the study.

Author contributions

HL: Data curation, Writing – original draft, Writing – review & editing, Methodology, Validation. RS: Conceptualization, Data curation, Formal analysis, Funding acquisition, Methodology, Resources, Validation, Writing – review & editing. YX: Data curation, Formal analysis, Investigation, Writing – original draft. MZ: Conceptualization, Data curation, Investigation, Writing – original draft. HD: Conceptualization, Formal analysis, Funding acquisition, Methodology, Resources, Supervision, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research and/or publication of this article. This study was Supported by Natural Science Foundation of Inner Mongolia Autonomous Region of China (2023MS08031) and Science and Technology Program of the Joint Fund of Scientific Research for the Public Hospitals of Inner Mongolia Academy of Medical Sciences (2024GLLH1321).

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.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fped.2025.1695974/full#supplementary-material

References

1. Wang YS, Zhou YL, Bai GN, Li SX, Xu D, Chen LN, et al. Expert consensus on the diagnosis and treatment of macrolide-resistant Mycoplasma pneumoniae pneumonia in children. World J Pediatr. (2024) 20:901–14. doi: 10.1007/s12519-024-00831-0

PubMed Abstract | Crossref Full Text | Google Scholar

2. Yan C, Xue GH, Zhao HQ, Feng YL, Cui JH, Yuan J. Current status of Mycoplasma pneumoniae infection in China. World J Pediatr. (2024) 20:1–4. doi: 10.1007/s12519-023-00783-x

PubMed Abstract | Crossref Full Text | Google Scholar

3. Kim K. Global trends in the proportion of macrolide-resistant Mycoplasma pneumoniae infections: a systematic review and meta-analysis. JAMA Netw Open. (2022) 5:e2220949. doi: 10.1001/jamanetworkopen.2022.20949

PubMed Abstract | Crossref Full Text | Google Scholar

4. Yen MH, Yan DC, Wang CJ, Tsao KC, Huang YC, Cheng SW, et al. The clinical significance of and the factors associated with macrolide resistance and poor macrolide response in pediatric Mycoplasma pneumoniae infection: a retrospective study. J Microbiol Immunol Infect. (2023) 56(3):634–40. doi: 10.1016/j.jmii.2023.01.010

PubMed Abstract | Crossref Full Text | Google Scholar

5. Meyer Sauteur PM, Unger WW, Nadal D, Berger C, Vink C, van Rossum AM. Infection with and carriage of Mycoplasma pneumoniae in children. Front Microbiol. (2016) 7:329. doi: 10.3389/fmicb.2016.00329

PubMed Abstract | Crossref Full Text | Google Scholar

6. Chen YC, Hsu WY, Chang TH. Macrolide-resistant Mycoplasma pneumoniae infections in pediatric community-acquired pneumonia. Emerg Infect Dis. (2020) 26(7):1382–91. doi: 10.3201/eid2607.200017

PubMed Abstract | Crossref Full Text | Google Scholar

7. National Health Commission of the People’s Republic of China. Guidelines for the diagnosis and treatment of Mycoplasma pneumoniae pneumonia in children (2023 edition). Electron J Emerg Infect Dis. (2024) 9(1):73–9. doi: 10.19871/j.cnki.xfcrbzz.2024.01.015

Crossref Full Text | Google Scholar

8. Sun Q, Hao J, Zhou Q, Zeng Y. Clinical characteristics and treatment outcomes of pediatric drug-resistant Mycoplasma pneumoniae pneumonia in the post-COVID-19 era. IDR. (2025) 18:957–63. doi: 10.2147/IDR.S502937

Crossref Full Text | Google Scholar

9. Tanaka H. Correlation between radiological and pathological findings in patients with Mycoplasma pneumoniae pneumonia. Front Microbiol. (2016) 7:695. doi: 10.3389/fmicb.2016.00695

PubMed Abstract | Crossref Full Text | Google Scholar

10. Ryu JH. Classification and approach to bronchiolar diseases. Curr Opin PulmMed. (2006) 12(2):145–51. doi: 10.1097/01.mcp.0000208455.80725.2a

PubMed Abstract | Crossref Full Text | Google Scholar

11. Wang X, Lu J, Song Z, Zhou Y, Liu T, Zhang D. From past to future: bibliometric analysis of global research productivity on nomogram (2000–2021). Front Public Health. (2022) 10:997713. doi: 10.3389/fpubh.2022.997713

PubMed Abstract | Crossref Full Text | Google Scholar

12. Cheng S, Lin J, Zheng X, Yan L, Zhang Y, Zeng Q, et al. Development and validation of a simple-to-use nomogram for predicting refractory Mycoplasma pneumoniae pneumonia in children. Pediatr Pulmonol. (2020) 55:968–74. doi: 10.1002/ppul.24684

PubMed Abstract | Crossref Full Text | Google Scholar

13. Luo Y, Wang Y. Risk prediction model for necrotizing pneumonia in children with Mycoplasma pneumoniae pneumonia. J Inflamm Res. (2023) 16:2079–87. doi: 10.2147/JIR.S413161

PubMed Abstract | Crossref Full Text | Google Scholar

14. Zhang X, Sun R, Jia W, Li P, Song C. A new dynamic nomogram for predicting the risk of severe Mycoplasma pneumoniae pneumonia in children. Sci Rep. (2024) 14(1):8260. doi: 10.1038/s41598-024-58784-3

PubMed Abstract | Crossref Full Text | Google Scholar

15. Zhang X, Sun R, Jia W, Li P, Song C. Clinical characteristics of lung consolidation with Mycoplasma pneumoniae pneumonia and risk factors for Mycoplasma pneumoniae necrotizing pneumonia in children. Infect Dis Ther. (2024) 13(2):329–43. doi: 10.1007/s40121-023-00914-x

PubMed Abstract | Crossref Full Text | Google Scholar

16. Ding G, Zhang X, Vinturache A, van Rossum AMC, Yin Y, Zhang Y. Challenges in the treatment of pediatric Mycoplasma pneumoniae pneumonia. Eur J Pediatr. (2024) 183(7):3001–11. doi: 10.1007/s00431-024-05519-1

PubMed Abstract | Crossref Full Text | Google Scholar

17. Du H, Li J, Li X, Zhao J, Lu W, Zhang Q, et al. Impact of the COVID-19 pandemic on epidemiological and clinical characteristics of Mycoplasma pneumoniae pneumonia in children: a multicenter study from Hubei, China. Front Pediatr. (2024) 12:1388132. doi: 10.3389/fped.2024.1388132

PubMed Abstract | Crossref Full Text | Google Scholar

18. Akira M, Kitatani F, Lee YS, Kita N, Yamamoto S, Higashihara T, et al. Diffuse panbronchiolitis: evaluation with high-resolution CT. Radiology. (1988) 168:433–8. doi: 10.1148/radiology.168.2.3393662

PubMed Abstract | Crossref Full Text | Google Scholar

19. Eisenhuber E. The tree-in-bud sign. Radiology. (2002) 222(3):771–2. doi: 10.1148/radiol.2223991980

PubMed Abstract | Crossref Full Text | Google Scholar

20. He J, Liu M, Ye Z, Tan T, Liu X, You X, et al. Insights into the pathogenesis of Mycoplasma pneumoniae. Mol Med Rep. (2016) 14(5):4030–6. doi: 10.3892/mmr.2016.5765

PubMed Abstract | Crossref Full Text | Google Scholar

21. Iannuzo N, Dy ABC, Guerra S, Langlais PR, Ledford JG. The impact of CC16 on pulmonary epithelial-driven host responses during Mycoplasma pneumoniae infection in mouse tracheal epithelial cells. Cells. (2023) 12(15):1984. doi: 10.3390/cells12151984

PubMed Abstract | Crossref Full Text | Google Scholar

22. Rodriguez K, Hariri LP, VanderLaan P, Abbott GF. Imaging of small airways disease. Clin Chest Med. (2024) 45:475–88. doi: 10.1016/j.ccm.2024.02.016

PubMed Abstract | Crossref Full Text | Google Scholar

23. Wen X, Liu J, Li H, Zhao C, Zhao S. Clinicoradiologic features of Mycoplasma pneumoniae bronchiolitis in children. Pediatr Investig. (2019) 2(4):248–52. doi: 10.1002/ped4.12108

PubMed Abstract | Crossref Full Text | Google Scholar

24. Huang Y, Wen H, Liu Y, Xie H, Wu L, Fu Z, et al. Correlation analysis between clinical characteristics of adult patients with Mycoplasma pneumoniae pneumonia and 23S rRNA A2063 gene in alveolar lavage fluid. BMC Pulm Med. (2025) 25:319. doi: 10.1186/s12890-025-03704-y

PubMed Abstract | Crossref Full Text | Google Scholar

25. Zhu Y, Luo Y, Li L, Jiang X, Du Y, Wang J, et al. Immune response plays a role in Mycoplasma pneumoniae pneumonia. Front Immunol. (2023) 14:1189647. doi: 10.3389/fimmu.2023.1189647

PubMed Abstract | Crossref Full Text | Google Scholar

26. Qiu J, Ge J, Cao L. D-dimer: the risk factor of children’s severe Mycoplasma pneumoniae pneumonia. Front Pediatr. (2022) 10:828437. doi: 10.3389/fped.2022.828437

PubMed Abstract | Crossref Full Text | Google Scholar

27. Zhao Q, Zhang T, Zhu B, Bi Y, Jiang SW, Zhu Y, et al. Increasing age affected polymorphonuclear neutrophils in prognosis of Mycoplasma pneumoniae pneumonia. J Inflamm Res. (2021) 14:3933–43. doi: 10.2147/JIR.S321656

PubMed Abstract | Crossref Full Text | Google Scholar

28. Liu Y, Zhang X, Wang Y, Zhu C, Fan M, Dou X, et al. The role of granulocyte macrophage colony stimulating factor in hospitalized children with Mycoplasma pneumoniae pneumonia. J Infect Chemother. (2018) 24(10):789–94. doi: 10.1016/j.jiac.2018.06.003

PubMed Abstract | Crossref Full Text | Google Scholar

29. Sato T, Okamura T, Kojima-Shibata C, Kadowaki H, Suzuki E, Uenishi H, et al. Correlated response of peripheral blood cytokines with selection for reduced mycoplasma pneumonia of swine lesions in Landrace pigs. Anim Sci J. (2016) 87(4):477–83. doi: 10.1111/asj.12462

PubMed Abstract | Crossref Full Text | Google Scholar

30. Li M, Gao J. Correlation and clinical significance of changes in serum soluble P-selectin, D- dimer and platelet levels with the severity of Mycoplasma Pneumoniae infection in children. Altern Ther Health Med. (2025) 31(4):244–9. 38758155

PubMed Abstract | Google Scholar

31. Ye Y, Gao Z, Zhang Z, Chen J, Chu C, Zhou W. A machine learning model for predicting severe mycoplasma pneumoniae pneumonia in school-aged children. BMC Infect Dis. (2025) 25(1):570. doi: 10.1186/s12879-025-10958-8

PubMed Abstract | Crossref Full Text | Google Scholar

32. Koupenova M, Mick E, Mikhalev E, Benjamin EJ, Tanriverdi K, Freedman JE. Sex differences in platelet toll-like receptors and their association with cardiovascular risk factors. Arterioscler Thromb Vasc Biol. (2015) 35(4):1030–7. doi: 10.1161/ATVBAHA.114.304954

PubMed Abstract | Crossref Full Text | Google Scholar

33. Lu Y, Ch Y, Tang X. Investigation on the causes of elevated CK and LD in patients with severe acute respiratory syndrome. Int Med Health Guid News. (2004) 10(18):4–6.

Google Scholar

34. Choi J, Jun M, Jeon J, Kim H, Kim Y, Jeon C, et al. Global lineage evolution pattern of sars-cov-2 in Africa, America, Europe, and Asia: a comparative analysis of variant clusters and their relevance across continents. J Transl Int Med. (2023) 11:410–22. doi: 10.2478/jtim-2023-0118

PubMed Abstract | Crossref Full Text | Google Scholar

35. Xie J, Ye F, Deng X, Tang Y, Liang JY, Huang X, et al. Circular RNA, a promising new star of vaccine. J Transl Int Med. (2023) 11:372–81. doi: 10.2478/jtim-2023-0122

PubMed Abstract | Crossref Full Text | Google Scholar

36. Jiang Y, Zhuo C, Huang J, Tan Y, Guan Y. Etiological distribution and clinical features of fever of unknown origin with pulmonary lesions in South China. Clin Respir J. (2019) 13(11):683–92. doi: 10.1111/crj.13070

PubMed Abstract | Crossref Full Text | Google Scholar

37. Bai HX, Wang R, Xiong Z, Hsieh B, Chang K, Halsey K, et al. Artificial intelligence augmentation of radiologist performance in distinguishing COVID-19 from pneumonia of other origin at chest CT. Radiology. (2021) 299(1):E225. doi: 10.1148/radiol.2021219004

PubMed Abstract | Crossref Full Text | Google Scholar

38. Chu Y, Wang J, Xiong Y, Gao Y, Liu X, Luo G, et al. Point-annotation supervision for robust 3D pulmonary infection segmentation by CT-based cascading deep learning. Comput Biol Med. (2025) 187:109760. doi: 10.1016/j.compbiomed.2025.109760

PubMed Abstract | Crossref Full Text | Google Scholar

39. Wang H, Zhao L, Yang Y, Wang F, Ding H, He L, et al. Chest computed tomography radiomics for determining macrolide resistance-associated gene mutation Status in children with Mycoplasma pneumoniae pneumonia: a two-center study. Acad Radiol. (2024) 31(9):3774–82. doi: 10.1016/j.acra.2024.05.008

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: Mycoplasma pneumoniae pneumonia, macrolide-resistant, children, nomogram, prediction model, tree-in-bud sign, neutrophil, platelet-to-lymphocyte ratio

Citation: Liu H, Sa R, Xin Y, Zhang M and Du H (2025) Development and validation of a nomogram for early prediction of macrolide-unresponsive Mycoplasma pneumoniae pneumonia in children. Front. Pediatr. 13:1695974. doi: 10.3389/fped.2025.1695974

Received: 31 August 2025; Accepted: 5 November 2025;
Published: 20 November 2025.

Edited by:

Zhongjie Shi, Wayne State University, United States

Reviewed by:

Haoru Wang, Children’s Hospital of Chongqing Medical University, China
Jihui Chen, Shanghai Jiao Tong University, China
Dexin Yu, Shandong University, China

Copyright: © 2025 Liu, Sa, Xin, Zhang and Du. 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: Hai Du, bnRoYWl0eUAxMjYuY29t

These authors have contributed equally to this work and share first authorship

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.