- Department of Pediatrics Department, Dongyang People’s Hospital, Dongyang, Zhejiang, China
Background: This study investigated risk factors for plastic bronchitis (PB) in children with pneumonia and created a nomogram for early detection.
Methods: We retrospectively analyzed data from 487 children with pneumonia who underwent bronchoscopic alveolar lavage between 2018 and 2024. Children were divided into a PB group (n = 65) and a No-PB group (n = 422). Key indicators were identified using regression analysis, and a nomogram prediction model was developed. Its effectiveness was evaluated using receiver operating characteristic (ROC) curves, calibration curves, decision curve analysis (DCA), and the bootstrap (BS) method.
Results: A total of 65 patients (13.3%) out of 487 had PB. Multivariate logistic regression analysis indicated that mycoplasma infection, atelectasis, high fever peak, high neutrophil levels and elevated lactate dehydrogenase levels were risk factors for PB development. Additionally, the BS-ROC curve of the developed nomogram had an area under the curve value of 0.857 (95% confidence interval: 0.808–0.905). The calibration curve analysis revealed a strong agreement between the actual and predicted probabilities of PB development, and DCA highlighted the clinical relevance of the nomogram.
Conclusions: A nomogram for MP infection, fever peak, LDH, maximum neutrophils percentage (N%max), and atelectasis was established to predict the risk of PB in children with pneumonia. The nomogram effectively identifies PB early, and bronchoscopy is advised for high-risk children. However, the nomogram needs external validation before practical use.
1 Introduction
Approximately 344 million patients with lower respiratory tract infections are reported each year, resulting in approximately 502,000 fatalities among children under fiver years of age, making it the second highest cause of death in this demographic (1, 2). Airway mucus hypersecretion is a pathological change associated with pneumonia. If it develops as a bronchial mucus thrombus or leads to plastic bronchitis (PB), it can cause clinical relapse, delayed radiographic resolution, and long-term pulmonary sequelae (3, 4). Therefore, searching for specific clinical features and related markers in children with pneumonia is crucial to determine whether they develop PB at an early stage. Currently, models are related to PB caused by a single pathogen, such as mycoplasma (MP) infection (3, 4). However, increasing evidence also reports models on PB caused by other pathogens, including adenovirus (ADV), respiratory syncytial virus (RSV), bocavirus, and Haemophilus influenzae (5–9). Clinically, children with mixed infections were more prevalent than single infections, occurring in 46.7% of cases vs. 40.7% (10). Consequently, we did not classify the children based on a specific pathogen; instead, we included all children with pneumonia. A risk model for PB in children with pneumonia was constructed by identifying relevant markers based on clinical features and laboratory tests.
2 Methods
2.1 Participants
Patients who met the inclusion criteria at the Dongyang City People's Hospital from April 2018 to August 2024 were retrospectively reviewed. The inclusion criteria for patients were as follows: (1) Hospitalized in our hospital and were ≤14 years old; (2) diagnosed with pneumonia and treated with bronchoscopic alveolar lavage; (3) whose family members agreed and signed an informed consent. The exclusion criteria for patients were as follows: (1) Congenital malformations, including respiratory, circulatory, and urinary malformations; (2) congenital genetic diseases or chromosomal abnormalities; (3) complicated by infections in other parts, such as enteritis and infectious mononucleosis; (4) foreign body aspiration; (5) incomplete clinical data. Figure 1 presents a flowchart of our research.
This retrospective study was approved by the Ethics Committee of Dongyang People's Hospital (Approval No: Dong Ren Yi 2024-YX-156) and required written informed consent from the parents for each case.
2.2 Definition of plastic bronchitis
PB diagnosis was based on BAL condition. PB is characterized by the production of sticky, branching casts within the tracheobronchial tree, which can cause partial or complete airway obstruction. These casts can range from small, segmental casts within the bronchus to large ones that fill the entire airway (11).
2.3 Data collection
Two physicians gathered clinical data, blood tests, and alveolar lavage fluid results from all eligible children. These data included information on gender, age, body mass index (BMI), duration and peak of fever, pulmonary signs, whether corticosteroid was used before bronchoscopy, blood routine examination, coagulation function, pharyngeal pathogenic bacterial polymerase chain reaction (MP, ADV, RSV), alveolar lavage fluid culture, bronchoalveolar lavage fluid (BAL) MP, and PB presence, among others. Pleural effusion and atelectasis were diagnosed by radiologists and pediatricians based on lung CT findings.
2.4 Statistical analysis
Statistical analysis was performed using R (version 4.4.2) software. Continuous data are expressed as mean ± standard deviation or median and interquartile range, with a t-test or rank sum test applied depending on the data distribution. Categorical data are expressed as numbers (%) and analyzed using the chi-square or Fisher's exact test. The least absolute shrinkage and selection operator (LASSO) regression was used to screen for the most variable indices. Multivariate logistic regression was used to identify the primary risk factors, and a two-way method was used to determine independent risk factors. Based on these analyses, a nomogram was constructed. Receiver operating characteristic curves, calibration curves, and decision curve analysis (DCA) were calculated using bootstrapping (BS) (1,000 resamplings) to assess the predictive power and performance of the model. A P < 0.05 was considered statistically significant.
3 Results
3.1 Clinical characteristics of patients with PB
A total of 487 children were included in our study, of whom 65 developed PB and were included in the positive group (PB group), while the remaining 422 were included in the negative group (No-PB group). Statistically significant differences were observed between the two groups in fever duration, peak fever, atelectasis, pleural effusion, MP infection (pharynx and BAL), maximum neutrophils percentage (N%max), lymphocytes (L%max), eosinophilic granulocytes (E%max), platelet minimum, alanine aminotransferase, aspartate aminotransferase, and lactate dehydrogenase (LDH). However, no statistical differences were observed in other factors, such as corticosteroid use before bronchoscopy, age, BMI, and gender (Table 1).
3.2 LASSO regression analysis of PB
Seven indexes were identified by LASSO regression analysis: Peak temperature, atelectasis, pleural effusion, L%max, N%max, MP infection, and LDH (Figures 2A,B). Due to the clear correlation between L%max and N%max, we've opted to include N%max, as it is widely used in clinical practice.

Figure 2. (A) LASSO coefficient profile for 14 variables indicates that more coefficients shrink to zero as penalties increase. (B) Tenfold cross-validation and minimization criteria were used to select the optimal penalty, lambda. The binomial deviance curve was plotted against the log of lambda, with vertical lines indicating 1 standard error. Seven variables with non-zero coefficients were selected with the optimal lambda.
3.3 Logistic regression analysis of PB
Five factors were identified by multivariate logistic regression: Peak fever [odds ratio (OR) = 2.444, 95% confidence interval (CI): 1.527–4.079; P < 0.001], atelectasis (OR = 13.68, 95% CI: 4.225–49.06; P < 0.001), MP infection (OR = 9.363, 95% CI: 2.643–60.81; P = 0.003), N%max (OR = 1.036, 95% CI: 1.007–1.068; P = 0.021), and LDH (OR = 1.003, 95% CI: 1.001–1.006; P = 0.007) were independent risk factors. Pleural effusion was excluded (Figure 3).
3.4 Creation and verification of a predictive model for PB
Independent risk factors from multivariate logistic regression analysis were used to create a nomogram (Figure 4). Each predictor score in the nomogram was summed to obtain a total score, which was then used to determine the likelihood of PB occurrence, aiding in evaluating the clinical necessity of aggressive bronchoscopic intervention. Our analysis demonstrated substantial predictive accuracy and model discrimination, with BS area under the curve values of 0.857 (95% CI: 0.808–0.905) obtained from 1,000 resamplings (Figure 5). Furthermore, the calibration curve of the nomogram demonstrated a reasonably a good match between the actual occurrence of PB and the predicted likelihood, with a P-value of 0.21 from the Hosmer–Lemeshow test confirming this finding (Figure 6). The DCA of the nomogram indicated that net benefits were achieved over a threshold probability range of 0.02–0.98. These findings suggest that the model may be an excellent predictor of PB-complicated pneumonia in children (Figure 7).
4 Discussion
Infection is a major element in PB formation (11, 12). During the acute phase of childhood pneumonia, persistent mucus accumulation can lead to PB, resulting in partial or complete dysfunction of pulmonary ventilation, dyspnea, recurrent high fever, pleural effusion in some cases, atelectasis, and even respiratory failure requiring mechanical ventilation, which can be dangerous to life in severe conditions or result in long-term effects (5, 12, 13). Therefore, early recognition of PB is significant. However, most existing models focus on single pathogens, although clinical cases often present with mixed infections. Consequently, we included all children with pneumonia and searched for relevant specific markers based on clinical characteristics and laboratory tests to construct a risk model for PB in children with pneumonia.
The MP infection independently influenced the PB development in our model, unlike ADV, RSV, or bacterial infections.This process is likely to cause PB for the following reasons: (1) MP infection directly adheres to airway epithelial cells and causes damage via cytotoxic effects or immune mechanisms. (2) MP-induced PB development in mouse models is associated with selective modulation of lymphatic structural changes through the PI3K/AKT/mTOR signaling cascade and pathways involving VEGFR-3 (14). Bronchial mucus plug or PB models based on MP infection in children have been previously reported in the literature (3, 15, 16), whereas MP infection remains one of the primary pathogens associated with PB in mixed infections. The literature indicates that 13 out of 15 cases (86.7%) with bronchoscopic treatment for childhood complex PB were MP-positive (17). In addition, MP infection and other pathogenic bacteria can cause PB. Huang (12) reported 43 cases of PB, of which 14 were MP-positive and 19 were ADV-positive. Chen et al. (18) reported that the top three pathogens in 55 cases of PB were MP (41.8%), ADV (20.0%), and influenza B virus (10.9%). However, our model did not identify ADV infection as a risk factor for PB. This process may be related to the restricted number of ADV infection cases in our group, which comprised only 14 cases, or to the number of ADV vectors. It has been reported that only one-third of children with ADV-infected pneumonia have bronchial mucus embolism, and the ADV number in the alveolar lavage fluid of these children is significantly increased. In mixed infections, particularly those involving ADV combined with MP, the incidence of bronchial mucus embolism is significantly higher (5). The influenza virus is a common cause of PB (18), but we excluded it from our model due to a lack of routine testing. There are limited clinical reports on PB resulting from bacterial or RSV infections, and no models exist for PB caused by a single pathogen, highlighting the need for more research.
Our research indicates that high neutrophil levels increase the risk of developing PB, whereas high lymphocyte levels act as protective factors. Previous studies (16, 19) have found that MP-infected children have more neutrophils and fewer lymphocytes compared with controls. Liu (8) observed that the neutrophil count in the blood or alveolar lavage fluid was elevated in patients with MP infection compared to those with bocavirus infection; however, neutrophils predominated in both groups. This process may be related to the role of neutrophils in immune responses, where they are involved in inflammation and strengthening host defense capabilities. Once their action is complete, neutrophils are rapidly eliminated, which reduces inflammatory responses and promotes tissue healing in pneumonia through the activity of lymphocytes and macrophages (20, 21). Accordingly, increased neutrophil levels indicate a marked inflammatory response.
In this study, we found that elevated LDH levels and fever peaks were independent risk factors for PB, and fever duration was recognized as a risk factor in univariate analysis. This finding is consistent with those reported by Zhao and Zhang (3, 4, 19). They reported higher LDH levels, longer fever durations, and higher fever peaks in patients with PB and MP pneumonia. However, PB caused by ADV infection, as reported by Yuan (22), also exhibited significantly increased LDH levels and recurrent fever. We believe that increased body temperature in children can lead to increased breathing and heart rate without a corresponding increase in water loss. When fluid intake is insufficient, mucus secretion can thicken, making the development of PB easier. In cases of pneumonia or airway obstruction in children, ischemia and hypoxia trigger anaerobic glycolysis, causing a significant release of LDH into the blood. Elevated LDH levels are strongly linked to the severity and prognosis of mycoplasma and adenovirus pneumonia, as shown by meta-analysis studies (23, 24).
In univariate analysis, the risk factors for PB included atelectasis and pleural effusion; however, multivariate analysis suggested that atelectasis was an independent risk factor, whereas pleural effusion was not. This finding is consistent with that of Zhao (19). While Zhang reported that pleural effusion is an independent risk factor for PB, he did not include atelectasis in his analysis (3). Among the 69 cases of PB reported by Lin (25), 35 had pleural effusion (50.72%), whereas 34 did not. In the pleural effusion group, C-reactive protein and LDH levels were significantly higher, suggesting a potential relationship with the degree of the inflammatory reaction (25). Our results conclude that pleural effusion and atelectasis are risk factors for PB, with atelectasis being more specifically linked to this condition. This is related to PB, characterized by infection or other factors that lead to the production of large, gelatinous, or rigid branched airway casts in the trachea. These casts can cause partial or complete airway obstruction, resulting in atelectasis.
Currently, the mechanisms underlying PB remain unclear; however, the most effective and direct treatment method is to remove the plastic foreign body obstructing the airway using bronchoscopy (12, 17). Bronchoscopy is an invasive procedure that involves risks such as mucosal hemorrhage and hypoxemia (26). Besides, it is relatively difficult to identify PB in children with pneumonia at an early stage, making it extremely important to evaluate its presence (27). In this study, a nomogram was created and confirmed to predict PB development in pediatric pneumonia cases using multivariate logistic regression to identify five key factors: MP infection, fever spike, N%max, and atelectasis. The nomogram demonstrated good performance, reflecting its accuracy and discriminatory ability. These five data points, accessible during a child's hospital stay, assist clinicians in early PB detection. Suppose the following two situations: 1. Consider a patient with Mycoplasma positivity, a peak temperature of 39.8°C, 85% neutrophils, atelectasis on CT, and an LDH level of 680 U/L. The scoring is as follows: Mycoplasma: 45 points, Temperature: 60 points, Neutrophils: 60 points, Atelectasis: 55 points, LDH: 40 points. Total score: 260 points, indicating a 0.91 probability and high risk, requiring immediate treatment. 2. For another patient who is Mycoplasma negative, with the same peak temperature of 39.8°C, 70% neutrophils, no atelectasis, and an LDH level of 680 U/L, the scoring is as follows: Mycoplasma: 0 points, Temperature: 60 points, Neutrophils: 50 points, Atelectasis: 0 points, LDH: 40 points. Total score: 150, indicating a 0.05 probability and low risk. For very low PB risk, prioritize drug therapy. For high-risk children, it is strongly recommended to perform early bronchoscopy to prevent complications.
However, our model has certain limitations. Firstly, this investigation was a retrospective analysis conducted at a single center. Despite rigorously adhering to the inclusion and exclusion criteria outlined in the study and employing multiple logistic regression to minimize confounding factors, the potential for selection bias remains. Secondly, our pathogen testing was confined to bacterial culture results from MP, RSV, ADV, and alveolar lavage fluid. Given that clinical presentations often involve mixed infections, there is a possibility of concurrent infections with other undetected viruses or bacteria, thereby imposing limitations on our findings. Lastly, the model is based on data from a single center with a relatively small sample size and lacks external validation, which may constrain its generalizability. Consequently, conducting multicenter and prospective studies is recommended as a subsequent step to enhance and refine the model.
5 Conclusion
A nomogram for MP infection, fever peak, LDH, N%max, and atelectasis was established to predict the risk of PB in children with pneumonia. This model will be helpful for guiding active bronchoscopy in children with pneumonia. But prospective external validation is needed before putting this model into practice.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving humans were approved by Ethics Committee of Dongyang People's Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants' legal guardians/next of kin. Written informed consent was obtained from the minor(s)' legal guardian/next of kin for the publication of any potentially identifiable images or data included in this article.
Author contributions
XF: Data curation, Formal analysis, Methodology, Project administration, Software, Visualization, Writing – original draft, Writing – review & editing. HL: Conceptualization, Data curation, Methodology, Resources, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare no that financial support was received for the research and/or publication of this article.
Acknowledgments
The authors greatly thank home-for-researchers, Song Ge Statistics and other pediatricians for their support.
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.
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.
References
1. Perin J, Mulick A, Yeung D, Villavicencio F, Lopez G, Strong KL, et al. Global, regional, and national causes of under-5 mortality in 2000–19: an updated systematic analysis with implications for the sustainable development goals. Lancet Child Adolesc Health. (2022) 6(2):106–15. doi: 10.1016/s2352-4642(21)00311-4
2. Bender RG, Sirota SB, Swetschinski LR, Dominguez RM, Novotney A, Wool EE, et al. Global, regional, and national incidence and mortality burden of non-COVID-19 lower respiratory infections and aetiologies, 1990–2021: a systematic analysis from the global burden of disease study 2021. Lancet Infect Dis. (2024) 24(9):974–1002. doi: 10.1016/s1473-3099(24)00176-2
3. Zhang H, Yang J, Zhao W, Zhou J, He S, Shang Y, et al. Clinical features and risk factors of plastic bronchitis caused by refractory mycoplasma pneumoniae pneumonia in children: a practical nomogram prediction model. Eur J Pediatr. (2023) 182(3):1239–49. doi: 10.1007/s00431-022-04761-9
4. Yang L, Zhang Y, Shen C, Lu Z, Hou T, Niu F, et al. Clinical features and risk factors of plastic bronchitis caused by mycoplasma pneumoniae pneumonia in children. BMC Pulm Med. (2023) 23(1):468. doi: 10.1186/s12890-023-02766-0
5. Peng L, Liu S, Xie T, Li Y, Yang Z, Chen Y, et al. Predictive value of adenoviral load for bronchial mucus plugs formation in children with adenovirus pneumonia. Can Respir J. (2022) 2022:9595184. doi: 10.1155/2022/9595184
6. Wang W, Zhang L, Ma W-K, He Y-X, Wang W-J, Su G-Y, et al. Plastic bronchitis associated with respiratory syncytial virus infection: a case report. BMC Pediatr. (2023) 23(1):517. doi: 10.1186/s12887-023-04351-0
7. Fei-Zhou Z, Mei-Xia H, Xiaofen T, Lei W, Xuan J, Lan-Fang T. Plastic bronchitis associated with human bocavirus 1 infection in children. Pediatr Pulmonol. (2024) 59(11):2754–60. doi: 10.1002/ppul.27113
8. Liu M, Wei D, Zhang T, Xu Y, Guo W. Distinct clinical characteristics of bocavirus and mycoplasma pneumoniae infection in children plastic bronchitis. Immun Inflamm Dis. (2024) 12(8):e1373. doi: 10.1002/iid3.1373
9. Yamasaki K, Morimoto T, Hashimoto K, Yamaguchi S, Kawamura M, Nemoto K, et al. Plastic bronchitis caused by haemophilus influenzae. Respirol Case Rep. (2023) 11(12):e01248. doi: 10.1002/rcr2.1248
10. Zhao L, Dong X, Cao L, Jiang Y, Zhang H, Mo R, et al. Multicenter evaluation of the acaruiter respiratory panel for diagnosis of respiratory tract infections in Chinese children. Microbiol Spectr. (2023) 11(6):e0058923. doi: 10.1128/spectrum.00589-23
11. Madsen P, Shah SA, Rubin BK. Plastic bronchitis: new insights and a classification scheme. Paediatr Respir Rev. (2005) 6(4):292–300. doi: 10.1016/j.prrv.2005.09.001
12. Huang J, Yang X, Zhuo Z, Yuan L. Clinical characteristics of plastic bronchitis in children: a retrospective analysis of 43 cases. Respir Res. (2022) 23(1):51. doi: 10.1186/s12931-022-01975-1
13. Cheng Q, Zhang H, Shang Y, Zhao Y, Zhang Y, Zhuang D, et al. Clinical features and risk factors analysis of bronchitis obliterans due to refractory mycoplasma pneumoniae pneumonia in children: a nomogram prediction model. BMC Infect Dis. (2021) 21(1):1085. doi: 10.1186/s12879-021-06783-4
14. Jin P, Zhao L-S, Zhang T-Q, Di H, Guo W. Establishment of a mouse model of mycoplasma pneumoniae-induced plastic bronchitis. Microorganisms. (2024) 12(6):1132. doi: 10.3390/microorganisms12061132
15. Xu X, Li H, Sheng Y, Wu L, Wang D, Liu L, et al. Nomogram for prediction of bronchial mucus plugs in children with mycoplasma pneumoniae pneumonia. Sci Rep. (2020) 10(1):4579. doi: 10.1038/s41598-020-61348-w
16. Zhang J, Wang T, Li R, Ji W, Yan Y, Sun Z, et al. Prediction of risk factors of bronchial mucus plugs in children with mycoplasma pneumoniae pneumonia. BMC Infect Dis. (2021) 21(1):67. doi: 10.1186/s12879-021-05765-w
17. Wang L, Wang W, Sun JM, Ni SW, Ding JL, Zhu YL, et al. Efficacy of fiberoptic bronchoscopy and bronchoalveolar lavage in childhood-onset, complicated plastic bronchitis. Pediatr Pulmonol. (2020) 55(11):3088–95. doi: 10.1002/ppul.25016
18. Chen X, Wu S, Huang Z, Lin Y, Xu J, Xu Q, et al. Plastic bronchitis in children: a review of 55 cases over a 10-year period. Int J Pediatr. (2024) 2024:9271324. doi: 10.1155/2024/9271324
19. Zhao L, Zhang T, Cui X, Zhao L, Zheng J, Ning J, et al. Development and validation of a nomogram to predict plastic bronchitis in children with refractory mycoplasma pneumoniae pneumonia. BMC Pulm Med. (2022) 22(1):253. doi: 10.1186/s12890-022-02047-2
20. Moldoveanu B, Otmishi P, Jani P, Walker J, Sarmiento X, Guardiola J, et al. Inflammatory mechanisms in the lung. J Inflamm Res. (2009) 2:1–11.22096348
21. Xiao Y, Cheng Y, Liu W-J, Liu K, Wang Y, Xu F, et al. Effects of neutrophil fate on inflammation. Inflammation Res. (2023) 72(12):2237–48. doi: 10.1007/s00011-023-01811-2
22. Yuan L, Huang JJ, Zhu QG, Li MZ, Zhuo ZQ. Plastic bronchitis associated with adenovirus serotype 7 in children. BMC Pediatr. (2020) 20(1):268. doi: 10.1186/s12887-020-02119-4
23. Zou M, Zhai Y, Mei X, Wei X. Lactate dehydrogenase and the severity of adenoviral pneumonia in children: a meta-analysis. Front Pediatr. (2022) 10:1059728. doi: 10.3389/fped.2022.1059728
24. Wang S, Jiang Z, Li X, Sun C, Zhang Y, Xiao Z. Diagnostic value of serum Ldh in children with refractory mycoplasma pneumoniae pneumoniae: a systematic review and meta-analysis. Front Pediatr. (2023) 11:1094118. doi: 10.3389/fped.2023.1094118
25. Lin X, Xu E, Zhang T, Zhu Q, Zhuang D. Risk factors for pleural effusion in children with plastic bronchitis caused by pneumonia. Front Pediatr. (2024) 12:1412729. doi: 10.3389/fped.2024.1412729
26. Zhou GW, Zhang W, Dong YC, Huang HD, Hu C, Sun J, et al. Flexible bronchoscopy-induced massive bleeding: a 12-year multicentre retrospective cohort study. Respirology. (2016) 21(5):927–31. doi: 10.1111/resp.12784
Keywords: pneumonia, plastic bronchitis, risk factors, nomogram, atelectasis, mycoplasma
Citation: Fang X and Lu H (2025) Nomogram for prediction of plastic bronchitis in Chinese children with pneumonia. Front. Pediatr. 13:1571479. doi: 10.3389/fped.2025.1571479
Received: 5 February 2025; Accepted: 30 April 2025;
Published: 13 May 2025.
Edited by:
Bülent Taner Karadağ, Marmara University, TürkiyeReviewed by:
Vladimir Pohanka, Retired, Liptovská Teolička, SlovakiaTerry Harville, University of Arkansas for Medical Sciences, United States
Copyright: © 2025 Fang and Lu. 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: Hemin Lu, SGVtaW5rZXlAMTYzLmNvbQ==