Neonates With Tracheomalacia Generate Auto-Positive End-Expiratory Pressure via Glottis Closure

Comparison of weighting algorithms to mitigate respiratory motion in free-breathing neonatal pulmonary radial UTE-MRI

Background. Thoracoabdominal MRI is limited by respiratory motion, especially in populations who cannot perform breath-holds. One approach for reducing motion blurring in radially-acquired MRI is respiratory gating. Straightforward 'hard-gating' uses only data from a specified respiratory window and suffers from reduced SNR. Proposed 'soft-gating' reconstructions may improve scan efficiency but reduce motion correction by incorporating data with nonzero weight acquired outside the specified window. However, previous studies report conflicting benefits, and importantly the choice of soft-gated weighting algorithm and effect on image quality has not previously been explored. The purpose of this study is to map how variable soft-gated weighting functions and parameters affect signal and motion blurring in respiratory-gated reconstructions of radial lung MRI, using neonates as a model population.Methods. Ten neonatal inpatients with respiratory abnormalities were imaged using a 1.5 T neonatal-sized scanner and 3D radial ultrashort echo-time (UTE) sequence. Images were reconstructed using ungated, hard-gated, and several soft-gating weighting algorithms (exponential, sigmoid, inverse, and linear weighting decay outside the period of interest), with %Nprojrepresenting the relative amount of data included. The apparent SNR (aSNR) and motion blurring (measured by the maximum derivative of image intensity at the diaphragm, MDD) were compared between reconstructions.Results. Soft-gating functions produced higher aSNR and lower MDD than hard-gated images using equivalent %Nproj, as expected. aSNR was not identical between different gating schemes for given %Nproj. While aSNR was approximately linear with %Nprojfor each algorithm, MDD performance diverged between functions as %Nprojdecreased. Algorithm performance was relatively consistent between subjects, except in images with high noise.Conclusion. The algorithm selection for soft-gating has a notable effect on image quality of respiratory-gated MRI; the timing of included data across the respiratory phase, and not simply the amount of data, plays an important role in aSNR. The specific soft-gating function and parameters should be considered for a given imaging application's requirements of signal and sharpness.

Scale resolving simulations of the effect of glottis motion and the laryngeal jet on flow dynamics during respiration

BACKGROUND AND OBJECTIVE

The movement of the respiratory walls has a significant impact on airflow through the respiratory tract. The majority of computational fluid dynamics (CFD) studies assume a static geometry which may not provide a realistic flow field. Furthermore, many studies use Reynolds Averaged Navier-Stokes (RANS) turbulence models that do not resolve turbulence structure. Combining the application of advanced scale-resolving turbulence models with moving respiratory walls using CFD will provide detailed insights into respiratory flow structures.

METHODS

This study simulated a complete breathing cycle involving inhalation and exhalation in a nasal cavity to trachea geometry that incorporated moving glottis walls. A second breathing cycle was simulated with static glottis walls for comparison. A recently developed hybrid RANS-LES turbulence model, the Stress-Blended Eddy Simulation (SBES), was incorporated to resolve turbulent flow structures in fine detail for both transient simulations. Transient results were compared with steady-state RANS simulations for the same respiratory geometry.

RESULTS

Glottis motion caused substantial effects on flow structure through the complete breathing cycle. Significant flow structure and velocity variations were observed due to glottal motion, primarily in the larynx and trachea. Resolved turbulence structures using SBES showed an intense mixing section in the glottis region during inhalation and in the nasopharynx during expiration, which was not present in the RANS simulations.

CONCLUSION

Transient simulations of a realistic breathing cycle uncovered flow structures absent in simulations with a constant flow rate. Furthermore, the incorporation of glottis motion impacted airflow characteristics that suggest rigid respiratory walls do not accurately describe respiratory flow. Future research in respiratory airflow should be conducted using transient scale-resolving models in conjunction with moving respiratory walls to capture flow structures in detail.

Tracheomalacia Reduces Aerosolized Drug Delivery to the Lung

Rationale: Neonates with respiratory issues are frequently treated with aerosolized medications to manage lung disease or facilitate airway clearance. Dynamic tracheal collapse (tracheomalacia [TM]) is a common comorbidity in these patients, but it is unknown whether the presence of TM alters the delivery of aerosolized drugs. Objectives: To quantify the effect of neonatal TM on the delivery of aerosolized drugs. Methods: Fourteen infant subjects with respiratory abnormalities were recruited; seven with TM and seven without TM. Respiratory-gated 3D ultrashort echo time magnetic resonance imaging (MRI) was acquired covering the central airway and lungs. For each subject, a computational fluid dynamics simulation modeled the airflow and particle transport in the central airway based on patient-specific airway anatomy, motion, and airflow rates derived from MRI. Results: Less aerosolized drug reached the distal airways in subjects with TM than in subjects without TM: of the total drug delivered, less particle mass passed through the main bronchi in subjects with TM compared with subjects without TM (33% vs. 47%, p = 0.013). In subjects with TM, more inhaled particles were deposited on the surface of the airway (48% vs. 25%, p = 0.003). This effect becomes greater with larger particle sizes and is significant for particles with a diameter >2 μm (2-5 μm, p ≤ 0.025 and 5-15 μm, p = 0.004). Conclusions: Neonatal patients with TM receive less aerosolized drug delivered to the lungs than subjects without TM. Currently, infants with lung disease and TM may not be receiving adequate and/or expected medication. Particles >2 μm in diameter are likely to deposit on the surface of the airway due to anatomical constrictions such as reduced tracheal and glottal cross-sectional area in neonates with TM. This problem could be alleviated by delivering smaller aerosolized particles.

Pulmonary function of children with tracheomalacia and associated clinical factors

OBJECTIVES

Spirometry is easily accessible yet there is limited data in children with tracheomalacia. Availability of such data may inform clinical practice. We aimed to describe spirometry indices of children with tracheomalacia, including Empey index and flow-volume curve pattern, and determine whether these indices relate with bronchoscopic features.

METHODS

From the database of children with tracheomalacia diagnosed during 2016-2019, we reviewed their flexible bronchoscopy and spirometry data in a blinded manner. We specially evaluated several spirometry indices and tracheomalacia features (cross-sectional lumen reduction, malacic length, and presence of bronchomalacia) and determined their association using multivariable regression.

RESULTS

Of 53 children with tracheomalacia, the mean (SD) peak expiratory flow (PEF) was below the normal range [68.9 percent of predicted value (23.08)]. However, all other spirometry parameters were within normal range [Z-score forced expired volume in 1 s (FEV₁ ) = -1.18 (1.39), forced vital capacity (FVC) = -0.61 (1.46), forced expiratory flow between 25% and 75% of vital capacity​​​​​​ (FEF_25%-75% ) = -1.43 (1.10), FEV₁ /FVC = -1.04 (1.08)], Empey Index = 8.21 (1.59). The most common flow-volume curve pattern was the "knee" pattern (n = 39, 73.6%). Multivariable linear regression identified the presence of bronchomalacia was significantly associated with lower flows: FEV₁ [coefficient (95% CI) -0.78 (-1.54, -0.02)], FEF_25%-75% [-0.61 (-1.22, 0)], and PEF [-12.69 (-21.13, -4.25)], all p ≤ 0.05. Other bronchoscopic-defined tracheomalacia features examined (cross-sectional lumen reduction, malacic length) were not significantly associated with spirometry indices.

CONCLUSION

The "knee" pattern in spirometry flow-volume curve is common in children with tracheomalacia but other indices, including Empey index, cannot be used to characterize tracheomalacia. Spirometry indices were not significantly associated with bronchoscopic tracheomalacia features but children with tracheobronchomalacia have significantly lower flow than those with tracheomalacia alone.

Predicting tracheal work of breathing in neonates based on radiological and pulmonary measurements

Tracheomalacia is an airway condition in which the trachea excessively collapses during breathing. Neonates diagnosed with tracheomalacia require more energy to breathe, and the effect of tracheomalacia can be quantified by assessing flow-resistive work of breathing (WOB) in the trachea using computational fluid dynamics (CFD) modeling of the airway. However, CFD simulations are computationally expensive; the ability to instead predict WOB based on more straightforward measures would provide a clinically useful estimate of tracheal disease severity. The objective of this study is to quantify the WOB in the trachea using CFD and identify simple airway and/or clinical parameters that directly relate to WOB. This study included 30 neonatal intensive care unit subjects (15 with tracheomalacia and 15 without tracheomalacia). All subjects were imaged using ultrashort echo time (UTE) MRI. CFD simulations were performed using patient-specific data obtained from MRI (airway anatomy, dynamic motion, and airflow rates) to calculate the WOB in the trachea. Several airway and clinical measurements were obtained and compared with the tracheal resistive WOB. The maximum percent change in the tracheal cross-sectional area (ρ = 0.560, P = 0.001), average glottis cross-sectional area (ρ = -0.488, P = 0.006), minute ventilation (ρ = 0.613, P < 0.001), and lung tidal volume (ρ = 0.599, P < 0.001) had significant correlations with WOB. A multivariable regression model with three independent variables (minute ventilation, average glottis cross-sectional area, and minimum of the eccentricity index of the trachea) can be used to estimate WOB more accurately (R2 = 0.726). This statistical model may allow clinicians to estimate tracheal resistive WOB based on airway images and clinical data.NEW & NOTEWORTHY The work of breathing due to resistance in the trachea is an important metric for quantifying the effect of tracheal abnormalities such as tracheomalacia, but currently requires complex dynamic imaging and computational fluid dynamics simulation to calculate it. This study produces a method to predict the tracheal work of breathing based on readily available imaging and clinical metrics.