Correlation between clinical risk factors and tracheal intubation difficulty in infants with Pierre-Robin syndrome: A retrospective study

DOI: https://doi.org/10.21203/rs.2.20445/v2

Abstract

Background:  Difficult tracheal intubation is a common problem encountered by anesthesiologists in the clinic. This study was conducted to assess the difficulty of tracheal intubation in infants with Pierre Robin syndrome (PRS) by incorporating computed tomography (CT) to guide airway management for anesthesia. Methods:  In this retrospective study, we analyzed case-level clinical data and CT images of 96 infants with PRS. First, a clinically experienced physician labeled CT images, after which the color space conversion, binarization, contour acquisition, and area calculation processing were performed on the annotated files. Finally, the correlation coefficient between the seven clinical factors and tracheal intubation difficulty, as well as the differences in each risk factor under tracheal intubation difficulty were calculated. Results:  The absolute value of the correlation coefficient between the throat area and tracheal intubation difficulty was 0.54; the observed difference was statistically significant. Body surface area, weight, and gender also showed significant difference under tracheal intubation difficulty. Conclusions: There is a significant correlation between throat area and tracheal intubation difficulty in infants with PRS. Body surface area, weight and gender may have an impact on tracheal intubation difficulty in infants with PRS.

Background

Difficult tracheal intubation is common in clinical practice, and it mostly refers to tracheal intubation that cannot be successfully completed by an ordinary indirect laryngoscope [1]. It represents the most difficult problem encountered by anesthesiologists in their daily work and is mainly caused by anatomical deformities, restricted back tilting activities, obesity and limited mouth opening [2]. These factors have an adverse effect on treatment. In practice, the level of difficulty is evaluated before the formal implementation of tracheal intubation. For patients with different levels of difficulty, preparations should be done in advance to avoid local mucosal damage caused by multiple intubation or complications such as dislocation of the circular cartilage [3].

In 2016, Münster et al [4] have reported that the position of vocal cords is related to laryngeal exposure and that difficult laryngoscopy is more likely to occur when vocal cords are closer to the head. From 2016 to 2018, many studies have utilized ultrasound for the clinical diagnosis of difficult tracheal intubation [5][6][7][8][9][10]. Ultrasound provides not only real-time images but also reveals dynamic structural changes of the airway. In 2019, Lee et al [11] found that the distance from the mandibular groove to the  hyoid bone and the distance from the inner edge of the mandible to the hyoid bone on X-ray images of the lateral neck were important for predicting difficult tracheal intubation in patients with acromegaly. However, there are only a few available  methods for infant airway assessment and their accuracy is relatively poor[12].

Pierre Robin syndrome [13][14] is the triad of micrognathia, glossoptosis, and cleft palate. These conditions could easily lead to difficult tracheal intubation which is the most significant risk factor for intubation anesthesia. Accurate preoperative prediction of intubation difficulty and adequate preparations are essential for ensuring successful airway management in infants with PRS[15][16]. There are many methods for assessing the difficulty of tracheal intubation [3][17]; yet, no existing method is suitable for infants, especially infants with PRS. Moreover, few reports have focused on the application of CT on tracheal intubation difficulty assessment in infants with PRS [18][19]. Therefore, this study was conducted to assess the difficulty of tracheal intubation in infants with PRS by incorporating CT to guide airway management for anesthesia[20].

Methods

Dataset

Clinical information and CT images were collected from infants with PRS who underwent intubation anesthesia in 2018 at Children’s Hospital of Nanjing Medical University. This retrospective study was approved by the local institutional review board that waived the requirement for informed consent based on minimal harm to the patient.

Seven clinical risk factors [21] that may have an impact on tracheal intubation difficulty were provided by experienced clinicians, including gender, height, weight, body surface area (BSA), throat area, age, and pneumonia (Table 1). The calculation of the throat area was elaborated below, and the remaining indicators could be directly obtained or simply calculated. Tracheal intubation difficulty is divided into three levels based on whether glottis can be completely observed under visual laryngoscope, where level Ⅰ refers to complete observation, level Ⅱ refers to partial observation, and level Ⅲ refers to the case when the only epiglottis can be observed.

Labeling criteria

To assess the impact of the throat area on tracheal intubation difficulty, the collected CT images (Figure 1a) were labeled according to the irregularity of the area being labeled using Labelme, an annotation tool which is based on the Python language and allows for irregular area annotation [22]. A radiologist with 20 years of clinical experience, who was blinded to the infants’ difficulty level, was responsible for labeling. Through a three-dimensional reconstruction technique, the median sagittal image of the upper airway of the infants was obtained, after which then the area of the oropharyngeal cavity (ie, the pharyngeal area between the plane of the tongue and the glottis) was labeled.

Annotation file processing and area calculation

The overall workflow is shown in Figure 2. The annotation file generated by Labelme is in the format of .json (Figure 1b) [23]. To calculate the throat area, the annotation file was first converted to a single-channel image in .png format (Figure 1c).

OpenCV performed subsequent processing in the Python environment. First, the single-channel image that was obtained during the previous step underwent color space conversion using the cvtColor function of OpenCV and was converted into a grayscale image (Figure 1d) [24][26]. The grayscale image was then thresholded (the threshold was set to 1) using the threshold function and becoming a binary image (Figure 1e) [25][26]. The throat contour information of the marker was then obtained by the findContours function, with pixel position difference between two adjacent points in all contour points no larger than 1 [26][27]. Finally, the contour information obtained in the previous step in the form of a point set was input into the contourArea function of OpenCV to calculate the area [26][28].

Correlation analysis

Correlation coefficients were used to assess the impact of each risk factor on tracheal intubation difficulty. Clinical risk factors highly correlated with difficulty level had better predicative effects in the clinic.

Statistical analysis

Since clinical risk factors include numerical and categorical variables and tracheal intubation difficulty is categorical, the correlation was measured by the Spearman rank correlation coefficient. Besides, to analyze whether there is a significant difference in each clinical risk factor under tracheal intubation difficulty, the Kruskal-Wallis test was used for numerical factors, and Pearson's Chi-squared test was used for categorical factors.

Results

The flow chart of the study is shown in Figure 2. Eight infants were excluded due to censored data (4 cases of censored pneumonia data and 4 cases of censored throat area data). Finally, 96 infants were included in the study, among whom 29 were level Ⅰ difficulty, 43 were level Ⅱ difficulty, and 24 were level Ⅲ difficulty of tracheal intubation. Additional data with sufficient clinical information were collected.

The correlation coefficients are integrated in Figure 3, where darker color indicates stronger correlations, while the lighter color represents weaker correlations. The correlation was strongest between the throat area and tracheal intubation difficulty with the correlation coefficient of -0.54. Risk factors that were moderately correlated with tracheal intubation difficulty were BSA, weight, and gender with correlation coefficients of -0.29, -0.29 and 0.26, respectively. All numerical risk factors were negatively correlated with tracheal intubation difficulty. Among categorical risk factors, males were more difficult to intubate than females, and infants with pneumonia had a lower level of difficulty in intubation than infants without pneumonia.

The results of the internal difference analysis of risk factors are shown in Table 2. The difference in throat area under tracheal intubation difficulty was significant, with P < 0.0001 (Level Ⅰ vs. Ⅱ: P = 0.0022, Level Ⅱ vs. Ⅲ: P = 0.0002, Level Ⅰ vs.Ⅲ: P < 0.0001). The differences in BSA, weight, and gender under tracheal intubation difficulty were also significant, and their corresponding P values were 0.0117, 0.0117 and 0.0043, respectively. BSA, weight, and gender were significantly different when comparing level Ⅱ to level Ⅲ and level Ⅰ to level Ⅲ. Height, age, and pneumonia showed no significant difference under tracheal intubation difficulty.

Discussion

In this study, we used clinical data from 96 PRS infants who underwent intubation anesthesia to perform correlation analysis, which demonstrated that the throat area had a significant effect on tracheal intubation difficulty. Our results revealed that a larger  throat area was associated with a lower level of tracheal intubation difficulty, which is consistent with the clinician's subjective perception. Besides, we found that high BSA and weight corresponded to low tracheal intubation difficulty, which may be related to the better physical development of these infants. Moreover, male infants had a higher tracheal intubation difficulty than females. Pneumonia, age, and height were slightly correlated with the difficulty of tracheal intubation, which may be due to the small amount of collected data and thus needs to be further analyzed.

After further P-value analysis, we found that four factors, namely throat area, gender, weight, and BSA, were internally different under the difficulty of tracheal intubation. Among them, the difference in the throat area was significant between all levels of tracheal intubation difficulty. Gender, weight, and BSA were only significantly different between level Ⅱ and level Ⅲ, level Ⅰ, and level Ⅲ. We speculate that it may be because the sample size of the level Ⅰ tracheal intubation difficulty is too small. In addition height, age, and pneumonia under tracheal intubation difficulty were not statistically significant, which may be related to the small sample size.

Attention should be paid to some of the limitations of our research. First, we studied the correlation between risk factors and tracheal intubation difficulty without building a predictive model, because the limited number of cases obtained in this study could not meet the requirements for modelling. Second, in order to facilitate the drawing of the correlation coefficient map, the correlation measure was based on the Spearman rank correlation coefficient. In addition, this was a single-center study. Finally, the annotation of the region of interest in the throat was done by one experienced doctor, which may be subjectively biased.

This study has few limitations: first, future studies should expand the number of cases collected and construct a predictive model of intubation difficulty. Secondly, the regional annotation should be performed by multiple physicians, and artificial intelligence annotation tools should be constructed. Finally, the integration of labeling and difficulty prediction should be performed.

Conclusion

The throat area may be helpful for predicting the difficulty of tracheal intubation in infants with PRS. Besides, gender, weight and BSA may also affect the prediction of the difficulty of airway intubation to some extent.

List of Abbreviations

PRS: Pierre Robin Syndrome

CT: Computed Tomography

BSA: Body Surface Area

Declarations

Ethics approval and consent to participate

The local institutional review committee (Nanjing, China) ethically approved the study. The committee waived the requirement for written informed consent.

Consent for publication

Not Applicable.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Competing interests

The authors declare that they have no competing interests.

Funding

The study was funded by departmental resources.

Authors' contributions

YlL is the main contributor in writing the manuscript. SZ is responsible for the collection and annotation of CT images. JSW processes the image and calculates the area, and performs statistical analysis. All authors read and approved the final manuscript.

Acknowledgements

Not applicable

References

  1. Xu Z, Ma W, Hester DL, et al. Anticipated and unanticipated difficult airway management. Curr Opin Anaesthesiol, 2018, 31 (1) :96-103.
  2. Cook TM, Woodall N, Frerk C. 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, 2011, 106 (5) :617-631.
  3. Rosenblatt WH. Preoperative planning of airway management in critical carepatients. Critical Care Medicine, 2004, 32 (4) :186-192.
  4. Münster T, Hoffmann M, Schlaffer S, et al. Anatomical location of the vocal cords in relation to cervical vertebrae. Eur J Anaesthesiol, 2016, 33 (4) :257-262.
  5. Guttman J, Nelson BP. Diagnostic emergency ultrasound: assessment techniques inthe pediatric patient. Pediatr Emerg Med Pract, 2016, 13 (1) :1-27.
  6. Erer OF, Erol S, Anar C, et al. Contribution of cell block obtained by endobronchial ultrasound-guided transbronchial needle aspiration in the diagnosis of malignant diseases and sarcoidosis. Endosc Ultrasound, 2017, 6 (4) :265-268.
  7. Leversedge FJ, Cotterell IH, Nickel B, et al. Ultrasonography guided de Quervain injection: accuracy and anatomic considerations in a cadaver model. J Am Acad Orthop Surg, 2016, 24 (6) :399-404.
  8. Li Y, Wang W, Yang T, et al. Incorporating uterine artery embolization in the treatment of cesarean scar pregnancy following diagnostic ultrasonography. Int J Gynaecol Obstet, 2016, 134 (2) :202-207.
  9. Osman A, Sum KM. Role of upper airway ultrasound in airway management. J Intensive Care, 2016, 4 (1) :52.
  10. Falcetta S, Cavallo S, Gabbanelli V, Pelaia P, Sorbello M, Zdravkovic I, Donati A. (2018). Evaluation of two neck ultrasound measurements as predictors of difficult direct laryngoscopy: A prospective study. Eur J Anaesthesiol. 2018; 35: 605-612. 
  11. Lee HC, Kim MK, Kim YH, et al. Radiographic predictors of difficult laryngoscopy in acromegaly patients. J Neurosurg Anesthesiol, 2019, 31 (1) :50-56.
  12. Karišik M, Janjević D, Sorbello M. Fiberoptic bronchoscopy versus video laryngoscopy in pediatric airway management. Acta Clin Croat. 2016 Mar;55 Suppl 1:51-4. 
  13. Benko S, Fantes JA, Amiel J, et al. Highly conserved non-coding elements on either side of SOX9 associated with Pierre Robin sequence[J].Nat Genet, 2009, 41 (3) :359-364.
  14. 2011 May; 127(5): 936-948. Robin Sequence: From Diagnosis to Development of an Effective Management Plan. Kelly N. Evans et al.
  15. J Clin Anesth. 2020 Jan 14;62:109726. doi: 10.1016/j.jclinane.2020.109726. A proposal for pediatric specific difficult airway management guidelines . Komasawa N, Matsunami S, Minami T.
  16. Int J Pediatr Otorhinolaryngol. 2019 Dec 28;130:109855. doi: 10.1016/j.ijporl.2019.109855. [Epub ahead of print] International Pediatric ORL Group (IPOG) Robin Sequence consensus recommendations. Fayoux P, et al.
  17. DAS Guidelines
  18. Adoración Martínez Plaza, Ricardo Fernández Valadés, Antonio España López, et al. Changes in airway dimensions after mandibular distraction in patients with Pierre-Robin sequence associated with malformation syndromes[J]. Revista Española De Cirugía Oral Y Maxilofacial, 2015, 37(2):71-79.
  19. Frova G, Guarino A, Petrini F, et Al. Recommendations for airway control and difficult airway management in paediatric patients. Minerva Anestesiol. 2006 Sep;72(9):723-48
  20. J Pediatr Intensive Care. 2018 Sep; 7(3): 115-125. Management of the Difficult Airway in the Pediatric PatientSenthil G. Krishna, Jason F. Bryant, and Joseph D. Tobias 
  21. Loftus PA , Ow TJ , Siegel B , et al. Risk Factors for Perioperative Airway Difficulty and Evaluation of Intubation Approaches Among Patients With Benign Goiter[J]. Annals of Otology, Rhinology & Laryngology, 2014, 123(4):279-285.
  22. Xue FS , Yuan YJ , Wang Q , et al. Difficulties and possible solutions for tracheal intubation with the Airway Scope[J]. American Journal of Emergency Medicine, 2011, 29(1):0-124.
  23. Hong L , Jin Q , Li X , et al. Image and medical annotations using non-homogeneous 2D ruler learning models[J]. Computers & Electrical Engineering, 2016, 50:102-110.
  24. Domínguez, César, Heras, Jónathan, Pascual V. IJ-OpenCV: Combining ImageJ and OpenCV for processing images in biomedicine[J]. Computers in Biology and Medicine, 2017, 84:189-194.
  25. Chernov V, Alander J, Bochko V. Integer-based accurate conversion between RGB and HSV color spaces[J]. Computers & Electrical Engineering, 2015, 46:328-337.
  26. Culjak I, Abram D, Pribanic T, et al. A brief introduction to OpenCV[C]// MIPRO, 2012 Proceedings of the 35th International Convention. IEEE, 2012.
  27. Shin JW. High-accuracy skin lesion segmentation and size determination[J]. Dissertations & Theses - Gradworks, 2011.
  28. Raymond WH, Garder A. A Spatial Filter for Use in Finite Area Calculations[J]. Monthly Weather Review, 2009, 116(1):209-222.

Tables

Table 1. Clinical information for children with PRS.

Gender

Male: Female

48: 48

Height (Unit: m)

Median (1st Qu., 3rd Qu )

0.5000 (0.5000, 0.5300)

Weight (Unit: kg)

Median (1st Qu., 3rd Qu )

3.400 (3.000, 3.800)

BSA (Unit: m2)

Median (1st Qu., 3rd Qu )

0.2190 (0.2050, 0.2330)

Throat area (Unit: pixel)

Median (1st Qu., 3rd Qu )

1440.5 (1237.2, 2034.4)

Age (Unit: day)

Median (1st Qu., 3rd Qu )

33.00 (13.67, 50.00)

Pneumonia

Yes: No

32: 64

Descriptive statistics of the seven clinical risk factors for 96 infants enrolled in the study. For categorical variables, the frequency of each category is listed. For numerical variables, the first quartile, median, and third quartiles are calculated.

 



Table 2. Difference analysis results of various factors.

 

Level 1 vs. 2

Level 2 vs. 3

Level 1 vs. 3

Total

Gender

1

0.0042 **

0.0125 *

0.0043**

Height

0.2473

0.4621

0.0526

0.1772

Weight

0.476

0.0264*

0.0025**

0.0117*

BSA

0.476

0.0264*

0.0025**

0.0117*

Throat area

0.0022**

0.0002 ***

<0.0001***

<0.0001***

Age

0.4694

0.2924

0.0503

0.1949

Pneumonia

0.4703

1

0.5253

0.5438

P-values for each risk factor under tracheal intubation difficulty. Among them, P values for a numerical variable were calculated by the Kruskal-Wallis test and for the categorical variable by Pearson's Chi-squared test.


*   P<0.05

**   P<0.01

***   P<0.001