Skip to main content

ORIGINAL RESEARCH article

Front. Physiol., 01 February 2024
Sec. Clinical and Translational Physiology

Tracheal intubation in patients with Pierre Robin sequence: development, application, and clinical value based on a 3-dimensional printed simulator

Updated
Yu Mao&#x;Yu Mao1Lu Liu&#x;Lu Liu2John ZhongJohn Zhong3Pei QinPei Qin4Rui MaRui Ma4Mingzhang ZuoMingzhang Zuo5Li Zhang
Li Zhang2*Lifang Yang
Lifang Yang4*
  • 1Department of Cardiovascular Surgery, Xijing Hospital, Air Force Medical University, Xi’an, China
  • 2Department of Anesthesiology, Children's Hospital of Nanjing Medical University, Nanjing, China
  • 3Department of Anesthesiology and Pain Management, University of Texas Southwestern Medical Center, Dallas, TX, United States
  • 4Department of Anesthesiology, Xi’an Children Hospital, Xi’an, China
  • 5Department of Anesthesia, Beijing Hospital, National Center of Gerontology, Institute of Geriatric Medicine, Chinese Academy of Medical Sciences, Beijing, China

Background: The main clinical manifestations of patients with Pierre Robin sequence (PRS) include micrognathia, the glossoptosis and dyspnoea. The difficulty of tracheal intubation (TI) in such patients is increased.

Objective: The purpose of the study was to evaluate the reliability and efficacy of the PRS simulator.

Methods: A PRS simulator was developed by using 3-dimensional (3D) printing technology under computer-aided design. A total of 12 anaesthesiologists each trained 5 times for TI on the PRS Training Simulator-1 and recorded the simulation time. After the training, they were randomly divided into three groups with a total of 12 nontrained anaesthesiologists, and the simulation was completed on PRS Simulator-2, 3 and 4. The simulation time was recorded, and the performance was evaluated by three chief anaesthesiologists. Then, all 24 anaesthesiologists completed the questionnaire.

Results: A PRS simulator developed by 3D printing was used to simulate the important aspects of TI. The average number of years worked was 6.3 ± 3.1 years, and 66.7% were female. The time for the 12 anaesthesiologists to complete the training gradually decreased (p < 0.01). Compared with the trained anaesthesiologists, the simulation time of TI in the nontrained anaesthesiologists was much longer (all p < 0.01). In addition, the simulation performance of the trained anaesthesiologists was relatively better (all p < 0.01).

Conclusion: The reliability and efficacy of the PRS simulator is herein preliminarily validated, and it has potential to become a teaching and training tool for anaesthesiologists.

Introduction

Pierre Robin sequence (PRS) is a rare congenital defect caused by mandibular dysplasia in new-borns characterised by micrognathia, glossoptosis and airway obstruction (Robin, 1923). PRS, if not handled in time, will lead to a series of serious complications (Cole et al., 2008). The treatment of PRS generally involves postural therapy, continuous positive pressure ventilation, nasopharyngeal ventilation, orthodontic treatment and other nonsurgical treatments. When nonsurgical treatment is ineffective or difficult to maintain, surgical treatment is needed. In recent years, mandibular distraction osteogenesis (MDO) has become the main surgical method for PRS treatment and has achieved good therapeutic effects, replacing tracheotomy and tongue and labial adhesions (TLAs) (Gómez et al., 2018).

Successful airway assessment and management is important in pre-MDO anaesthesia. However, the unique anatomy of PRS undoubtedly increases the challenge of tracheal intubation (TI), which in turn can lead to fatal anaesthesia-related complications (Artime and Hagberg, 2015). It has been reported that one in every 22,000 patients undergoing anaesthesia develops serious airway-related complications (Cook et al., 2011). Among them, the major associated adverse outcomes include death, hypoxic brain injury, cardiopulmonary arrest, airway trauma, aspiration of stomach contents, pulmonary oedema, and dental damage (Cook et al., 2011).

In recent years, 3-dimensional (3D) printing has become an easy technology to use and versatile. It is a technology that is based on the 3D reconstruction data formed by medical imaging scans and printed layer by layer of the corresponding model exactly consistent with the reconstructed model (MacDonald and Wicker, 2016). This technology has been widely applied in orthopaedics, neurosurgery and orofacial surgery, but its use in anaesthesia has not yet been full-scale (Sinn et al., 2006; Suomalainen et al., 2015; Mulford et al., 2016). Previous studies have shown that the use of 3D printing to create a patient’s airway model can help anaesthesiologists more fully understand the morphology of the trachea, choose the appropriate catheter, and determine a more secure intubation approach in advance (Fiadjoe et al., 2015; Park et al., 2021; Shaylor et al., 2023). Furthermore, in real-world procedures, the intubation time has been shortened, the success rate of intubation has been increased, and complications caused by repeated intubation have mainly been avoided (Han et al., 2016).

Therefore, to enhance familiarity of anaesthesiologists’s handling difficult airway such as PRS and potentially lower complications caused by mishandling, we designed an innovative, highly realistic PRS simulator to guide anaesthesiologists to simulate and practice key parts of procedures. The purpose of this study was to evaluate the reliability and efficacy of the simulator for preprocedural evaluation and training in such patients.

Methods

Study design

Four patients with PRS who received MDO treatment in Nanjing Children’s Hospital and Xi’an Children’s Hospital in May 2023 were enrolled. Based on their computed tomography (CT) data, 3D printed models were built and printed by professional engineers under computer-aided design, which were named PRS Training Simulator (Patient #1) and PRS Testing Simulator-1 (Patient #2), 2 (Patient #3), and 3 (Patient #4). Twenty-four anaesthesiologists from Xi’an Children’s Hospital and Nanjing Children’s Hospital were randomly assigned to participate in this study and to confirm the effects of the simulator (Figure 1). PRS Testing Simulator-2,3 and 4 were used in Groups A, B and C, respectively. This study protocol was approved by the Ethics Committee of the two centres (202301005-1), and all patient data were desensitised.

Figure 1
www.frontiersin.org

Figure 1. Flow chart. Red figures represent the trained anaesthesiologists. Orange, green and blue figures represent the nontrained anaesthesiologists in Groups A, B, and C, respectively.

Development of PRS simulator for tracheal intubation

The Digital Imaging and Communications in Medicine (DICOM) format of CT data of PRS patients was imported into Materialise Mimics Version 21.0 software (Leuven, Belgium), and the 3D reconstructed model of the head and airway using the threshold segmentation function was obtained. In Materialise 3-Matic software (Leuven, Belgium), the obtained 3D model was extracted, clipped, smoothed and repaired, and the above anatomical structures were completely restored 1:1. Furthermore, the Standard Triangle Language (STL) file of the 3D reconstructed model was exported to a Stratasys Polyjet 850 multimaterial full colour 3D printer (Supplementary Figure S1, Supplementary Table S1). Different tissues matched hard and soft materials with different colours for editing and printing to obtain a specific PRS simulator (Supplementary Table S1, Figures 2A–E).

Figure 2
www.frontiersin.org

Figure 2. Pierre Robin sequence simulator and its simulation procedures. (A–E) The simulator components. The simulator consisted of the head, neck and half of the upper body of the specific patient’s model. The mandible and the upper airway were detachable and replaceable. The simulator was used to simulate different types of laryngoscopes (F, G) to complete tracheal intubation (H, I). (J) Finally, the anaesthesiologist examined the lung function to verify the simulation effect.

Training and simulation

Twelve anaesthesiologists participated in the training performed 5 times each on the PRS Training Simulator, recording the time of each training (from the beginning use of the laryngoscope to the observance of thoracic undulating) (Figures 2F–J). After the training, they were randomly divided into three groups with a total of 12 nontrained anaesthesiologists, and the simulation was completed on PRS Testing Simulator-2, 3, 4. The training time was recorded, and the performance was evaluated by three chief anaesthesiologists (10-point scale). After the evaluation, all 24 anaesthesiologists filled out the Likert scale scorecard to evaluate the simulator’s performance and function in real-world procedures. The Likert scale contents include the overall perspective, the simulator material and operability, and the performance of simulation steps (Table 1).

Table 1
www.frontiersin.org

Table 1. Questionnaire for the evaluation of specific training by the Pierre Robin sequence simulator.

Statistical analysis

The measurement data conforming to a normal distribution are presented as the mean ± standard deviation. The rank sum test was used for rating, and the paired Student’s t-test was used for analysis. Bilateral p values <0.05 were considered statistically significant. All statistical analyses were performed using SPSS software version 26.0 (Armonk, United States).

Results

Study population

The demographic characteristics and preoperative CT measurement results of 4 patients with PRS are shown in Table 2. We measured the distance between the inner margins of the mandible (L1), the vertical distance from the inner margin of the mandible to L1 (L2), the distance of the narrowest level of the oropharyngeal region in the upper airway (L3), and the distance of the upper airway perpendicular to L3 (L4). Furthermore, we assessed the cross-sectional area of the glossopharyngeal region of the upper airway (CSA) and the angle between the two inner margins of the mandible and the inner margin of the mandibular body (α). Importantly, the patients’ clinical grades of Pierre Robin sequence were all up to Grade 3 (Lee et al., 2015).

Table 2
www.frontiersin.org

Table 2. Baseline characteristics of the patients.

Baseline characteristics of anaesthesiologists and PRS simulator assembly costs

The baseline characteristics of all 24 anaesthesiologists are shown in Table 3. The average age was 31.8 ± 4.3 years, 66.7% (n = 16) were female, and the average number of working years was 6.3 ± 3.1. There were no significant difference among the three groups in terms of average age and working years. Notably, 4 anaesthesiologists in each group were trained on the PRS Training Simulator, and there was no significant difference in the average age or working years between the trained anaesthesiologists and the nontrained anaesthesiologists. As an assembler simulator, the consumables for the detachable part cost 32.76 dollars (of which the mandible model cost 14.08 dollars and the upper airway model cost 18.68 dollars), and other fixed components cost 346.48 dollars (Supplementary Table S1).

Table 3
www.frontiersin.org

Table 3. Baseline characteristics of the anaesthesiologists.

Training, simulation and questionnaire outcomes

The simulated time results on the PRS Training Simulator are shown in Figure 3. As expected, the time for the 12 anaesthesiologists to complete the training using videoscope gradually decreased (4.3 ± 1.6 min vs. 2.8 ± 1.2 min vs. 2.2 ± 0.9 min vs. 1.4 ± 0.5 min vs. 0.8 ± 0.3 min, p < 0.01). The anaesthesiologists of Groups A, B and C completed a simulated operation of TI on PRS Simulator-1, 2 and 3, respectively. Table 4 shows the results of the simulation time. Compared with the nontrained anaesthesiologists, the trained anaesthesiologists had a significantly reduced simulation time (1.2 ± 0.7 vs. 4.3 ± 2.1, p < 0.01). Importantly, the simulated performance of the trained anaesthesiologists was significantly better than that of the nontrained anaesthesiologists (8.6 ± 2.2 vs. 5.4 ± 2.0, p < 0.01). Notably, all 24 anaesthesiologists approved the reliability and efficacy of the PRS simulator (Figure 4).

Figure 3
www.frontiersin.org

Figure 3. The duration of training on the Pierre Robin Sequence Training Simulator.

Table 4
www.frontiersin.org

Table 4. Simulation performance among the three groups.

Figure 4
www.frontiersin.org

Figure 4. The Likert scale results.

Discussion

This study presents an optimised TI simulator for anaesthesiologists. The results show that the PRS simulator can help significantly shorten the successful TI time defined by beginning of laryngoscope and spotting of chest undulating. The trained group was able to cut average time from 4.3 min at first attempt to 0.8 min at 5th one. When facing a different PRS Train Simulator, their successful TI time was also significantly shorter compared to the groups of not being simulated.

The main clinical features of PRS are micrognathia, glossoptosis, and airway obstruction (Robin, 1923), and approximately 58%–90% of patients have a cleft palate (Maas and Poets, 2014). The aetiology of PRS is not completely clear, but it may be related to inclusion body virus infection of giant cells in the uterus (Jakobsen et al., 2007). Previous studies reported that the pathogenesis of PRS may be based on mandibular dysplasia and bilateral symmetric retraction, resulting in a small oral capacity, an abnormal tongue position and secondary damage of palate closure. All those features plus rare occurrence of PRS make their airway handling challenging, a lot of anesthesiologists are not comfortable handing their airways due to lack of training or exposure. Commercially available mannequin are expensive, not readily available and provide only one dimension of airway difficulty (Parotto et al., 2017; Park et al., 2021). In our study, due to the anaesthesiologists could practice various 3D printed models of PRS, and the procedural plan may be formulated by using the simulator to improve the success rate of TI and reduce the potential complications from repeated airway instrumentation.

Difficult airways are one of the most urgent clinical challenges in anaesthesia, and it is very important to fully understand and prepare them before surgery. Based on the accurate design, 3D printing can be used to create highly simulated models, which is one of the essential factors that improve the quality of simulations and reduce the learning curve of anaesthesiologists. First, the accurate estimation of 3D spatial structures are the unique advantages of 3D printing. Anaesthesiologists often face problems with 3D anatomical structures, and 3D printing has shown its important role in dealing with these challenges (Gauger et al., 2018). In a previous study of cricothyroid incisions, 3D reconstruction may highly restore real anatomical structures and distinguish thyroid cartilage, cricoid cartilage and trachea (Gauger et al., 2018). Similarly, the evaluation of the chief anaesthesiologists in this study also endorsed the fidelity of the PRS simulator. The detachable part of the PRS simulator cost 32.76 dollars, and it was able to be used not only to train anaesthesiologists to grasp procedural skills quickly but also to simulate specific challenging cases before procedures.

Furthermore, 3D-printed models may not only be employed to visually identify anatomical structures but also help anaesthesiologists develop preprocedural airway management plans to reduce the incidence of intubation complications. Bustamante et al. (2014) and Parotto et al. (2017) produced a simulated trachea-bronchial tree model with 3D printing for the clinical training of residents to promote their proficiency in trachea-bronchial anatomy, fibreoptic bronchoscopy operation and lung isolation. Han et al. (2016) reported a case of 3D printing used for airway evaluation. A preoperative CT scan of this patient revealed a soft tissue mass at the tracheotomy position. The anaesthesiologist used 3D printing to develop a simulated airway and found scar contraction around the tracheal incision and slight narrowing below the incision. With the guidance of the model, they successfully achieved airway control. In addition, Park et al. (2021) reported a case of using a 3D printed model to simulate TI before anaesthesia to improve the quality of airway management. Meanwhile, Shaylor et al. (2023) carried out a randomized trial of the video laryngoscope to the flexible fiberoptic bronchoscope in a Pierre Robin manikin, the first-attempt success rates were 88.3% and 85.0%, respectively. In this study, the trained anaesthesiologists made significant improvements in simulation performance after training, and the training time on the training simulator was gradually reduced. Importantly, the training performance of the trained anaesthesiologists was also significantly better than that of their nontrained counterparts.

There are some limitations in this study. First, the sample size of this study is limited. We conducted the study and analysis on the models of only 4 PRS patients, and more patients should be included for evaluation in the future. Second, none of the anaesthesiologists performed TI on real PRS patients to further confirm the reliability of the simulator. Third, because the manufacturing of 3D-printed models is relatively time-consuming (9.5 days of average production time were used in the study), which needs to be reduced in the next step of simulator development. Finally, the physical properties of each anatomical structure of the simulator in terms of elongation at break of and tensile strength of various tissues are still different from those of the real human body, which needs to be improved in the next step.

Conclusion

In this study, a new potential training tool to optimise teaching and training for PRS cases was proposed. The simulator can help anaesthesiologists significantly shortening successful TI time without causing any harm.to real patients. The anaesthesiologists considered our PRS simulator is of high value in application and training. The performance of the trained anaesthesiologists was not only improved after training but also significantly better than that of nontrained anaesthesiologists. Although the simulator is not yet able to give anaesthesiologists a completely realistic simulation experience, it may greatly improve the success rate and efficacy of TI in PRS patients. Therefore, the reliable simulation quality of the simulator makes it promising for clinical applications in the future, and future studies will investigate the effects of using this simulator as a training tool and the potential impact on difficult airway management.

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 authors.

Ethics statement

Written informed consent was obtained from the individual(s), and minor(s)’ legal guardian/next of kin, for the publication of any potentially identifiable images or data included in this article.

Author contributions

YM: Writing–original draft. LL: Writing–original draft. JZ: Writing–review and editing. PQ: Writing–original draft, Data curation. RM: Writing–original draft, Data curation. MZ: Writing–review and editing. LZ: Writing–review and editing. LY: Writing–review and editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was supported by National Natural Science Foundation of China, Award Numbers: (No. 81774415, No. 82174493); Outstanding Youth Foundation of Shaanxi Province, Award Number: (No. 2021JC-49).

Acknowledgments

We would like to thank Make Medical Technology Co., LTD. (Xi’an, China) for supplying the 3D printed models.

Conflict of interest

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

Publisher’s note

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

Supplementary material

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

Abbreviations

CT, computed tomography; CSA, cross-sectional area; MDO, mandibular distraction osteogenesis; PRS, Pierre Robin sequence; 3D, 3-dimensional.

References

Artime C. A., Hagberg C. A. (2015). Is there a gold standard for management of the difficult airway? Anesthesiol. Clin. 33 (2), 233–240. doi:10.1016/j.anclin.2015.02.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Bustamante S., Bose S., Bishop P., Klatte R., Norris F. (2014). Novel application of rapid prototyping for simulation of bronchoscopic anatomy. J. Cardio⁃ Thorac. Vasc. Anesth. 28 (4), 1122–1125. doi:10.1053/j.jvca.2013.08.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Cole A., Lynch P., Slator R. (2008). A new grading of Pierre Robin sequence. Cleft Palate Craniofac J. 45 (6), 603–606. doi:10.1597/07-129.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Cook T. M., Woodall N., Frerk C.Fourth National Audit Project (2011). Major complications of airway management in the UK: results of the fourth national audit project of the royal college of anaesthetists and the difficult airway society. Part 1: anaesthesia. Br. J. Anaesth. 106 (5), 617–631. doi:10.1093/bja/aer058

PubMed Abstract | CrossRef Full Text | Google Scholar

Fiadjoe J. E., Hirschfeld M., Wu S., Markley J., Gurnaney H., Jawad A. F., et al. (2015). A randomized multi-institutional crossover comparison of the GlideScope® Cobalt Video laryngoscope to the flexible fiberoptic bronchoscope in a Pierre Robin manikin. Paediatr. Anaesth. 25 (8), 801–806. doi:10.1111/pan.12668

PubMed Abstract | CrossRef Full Text | Google Scholar

Gauger V. T., Rooney D., Kovatch K. J., Richey L., Powell A., Berhe H., et al. (2018). A multidisciplinary in⁃ ternational collaborative implementing low-cost, high-fidelity 3D printed airway models to enhance Ethiopian anesthesia resident emergency cricothyroidotomy skills. Int. J. Pediatr. Otorhinolar⁃ yngol 114, 124–128. doi:10.1016/j.ijporl.2018.08.040

PubMed Abstract | CrossRef Full Text | Google Scholar

Gómez O. J., Barón O. I., Peñarredonda M. L. (2018). Pierre Robin Sequence: an evidence based treatment proposal. J. Craniofac Surg. 29 (2), 332–338. doi:10.1097/SCS.0000000000004178

PubMed Abstract | CrossRef Full Text | Google Scholar

Han B., Liu Y., Zhang X., Wang J. (2016). Three-dimensional printing as an aid to airway evaluation after tracheotomy in a patient with laryngeal carcinoma. BMC Anesthesiol. 16, 6. doi:10.1186/s12871-015-0170-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Jakobsen L. P., Ullmann R., Christensen S. B., Jensen K. E., Mølsted K., Henriksen K. F., et al. (2007). Pierre Robin sequence may be caused by dysregulation of SOX9 and KCNJ2. J. Med. Genet. 44 (6), 381–386. doi:10.1136/jmg.2006.046177

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee S. L., Hosford C., Lee Q. T., Parnes S. M., Shapshay S. M. (2015). Mallampati class, obesity, and a novel airway trajectory measurement to predict difficult laryngoscopy. Laryngoscope 125 (1), 161–166. doi:10.1002/lary.24829

PubMed Abstract | CrossRef Full Text | Google Scholar

Maas C., Poets C. F. (2014). Initial treatment and early weight gain of children with Robin Sequence in Germany: a prospective epidemiological study. Arch. Dis. Child. Fetal Neonatal 99 (6), 491–494. doi:10.1136/archdischild-2014-306472

PubMed Abstract | CrossRef Full Text | Google Scholar

MacDonald E., Wicker R. (2016). Multiprocess 3D printing for increasing component functionality. Science 353 (6307), aaf2093. doi:10.1126/science.aaf2093

PubMed Abstract | CrossRef Full Text | Google Scholar

Mulford J. S., Babazadeh S., Mackay N. (2016). Three-dimensional printing in orthopaedic surgery: review of current and future applications. ANZ J. Surg. 86 (9), 648–653. doi:10.1111/ans.13533

PubMed Abstract | CrossRef Full Text | Google Scholar

Park S., Ahn J., Kim H. J., Choi E. J. (2021). Endotracheal intubation using a three-dimensional printed airway model in a patient with Pierre Robin sequence and a history of tracheostomy -a case report. Korean J. Anesthesiol. 74 (3), 262–265. doi:10.4097/kja.20430

PubMed Abstract | CrossRef Full Text | Google Scholar

Parotto M., Jiansen J. Q., AboTaiban A., Ioukhova S., Agzamov A., Cooper R., et al. (2017). Evaluation of a low⁃cost, 3D⁃printed model for bronchoscopy training. Anaes⁃ thesiol Intensive Ther. 49 (3), 189–197. doi:10.5603/AIT.a2017.0035

PubMed Abstract | CrossRef Full Text | Google Scholar

Robin P. (1923). Backword lowering of the root of the tongue causing respiratory disturbances. Bull. Acad. Med. 89 (2), 37–41.

Google Scholar

Shaylor R., Weiniger C. F., Rachman E., Sela Y., Kohn A., Lahat S., et al. (2023). A prospective observational crossover study comparing intubation by pediatric residents using video laryngoscopy and direct laryngoscopy on a Pierre Robin simulation manikin. Pediatr. Emerg. Care 39 (3), 159–161. doi:10.1097/PEC.0000000000002923

PubMed Abstract | CrossRef Full Text | Google Scholar

Sinn D. P., Cillo J. E., Miles B. A. (2006). Stereolithography for craniofacial surgery. J. Craniofac Surg. 17 (5), 869–875. doi:10.1097/01.scs.0000230618.95012.1d

PubMed Abstract | CrossRef Full Text | Google Scholar

Suomalainen A., Stoor P., Mesimäki K., Kontio R. K. (2015). Rapid prototyping modelling in oral and maxillofacial surgery: a two-year retrospective study. J. Clin. Exp. Dent. 7 (5), e605–e612. doi:10.4317/jced.52556

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Pierre Robin sequence, tracheal intubation, 3-dimensional printing, simulator, difficult airway

Citation: Mao Y, Liu L, Zhong J, Qin P, Ma R, Zuo M, Zhang L and Yang L (2024) Tracheal intubation in patients with Pierre Robin sequence: development, application, and clinical value based on a 3-dimensional printed simulator. Front. Physiol. 14:1292523. doi: 10.3389/fphys.2023.1292523

Received: 11 September 2023; Accepted: 11 December 2023;
Published: 01 February 2024.

Edited by:

Christina Maria Pabelick, Mayo Clinic, United States

Reviewed by:

Yary Volpe, University of Florence, Italy
Gunpreet Oberoi, Medical University of Vienna, Austria
Lindsay Guevara, Mayo Clinic, United States

Copyright © 2024 Mao, Liu, Zhong, Qin, Ma, Zuo, Zhang and Yang. 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: Li Zhang, drzhangli@njmu.edu.cn; Lifang Yang, 27891533@qq.com

These authors have contributed equally to this work

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