What can computed tomography and magnetic resonance imaging tell us about ventilation?

Imaging the Lung in ARDS: A Primer

Despite periodic changes in the clinical definition of ARDS, imaging of the lung remains a central component of its diagnostic identification. Several imaging modalities are available to the clinician to establish a diagnosis of the syndrome, monitor its clinical course, or assess the impact of treatment and management strategies. Each imaging modality provides unique insight into ARDS from structural and/or functional perspectives. This review will highlight several methods for lung imaging in ARDS, emphasizing basic operational and physical principles for the respiratory therapist. Advantages and disadvantages of each modality will be discussed in the context of their utility for clinical management and decision-making.

Respiratory Invariant Textures From Static Computed Tomography Scans for Explainable Lung Function Characterization

PURPOSE

The inherent characteristics of lung tissue independent of breathing maneuvers may provide fundamental information for function assessment. This paper attempted to correlate textural signatures from computed tomography (CT) with pulmonary function measurements.

MATERIALS AND METHODS

Twenty-one lung cancer patients with thoracic 4-dimensional CT, DTPA-single-photon emission CT ventilation ( VNM ) scans, and available spirometry measurements (forced expiratory volume in 1 s, FEV 1 ; forced vital capacity, FVC; and FEV 1 /FVC) were collected. In subregional feature discovery, function-correlated candidates were identified from 79 radiomic features based on the statistical strength to differentiate defected/nondefected lung regions. Feature maps (FMs) of selected candidates were generated on 4-dimensional CT phases for a voxel-wise feature distribution study. Quantitative metrics were applied for validations, including the Spearman correlation coefficient (SCC) and the Dice similarity coefficient for FM- VNM spatial agreement assessments, intraclass correlation coefficient for FM interphase robustness evaluations, and FM-spirometry comparisons.

RESULTS

At the subregion level, 8 function-correlated features were identified (effect size>0.330). The FMs of candidates yielded moderate-to-strong voxel-wise correlations with the reference VNM . The FMs of gray level dependence matrix dependence nonuniformity showed the highest robust (intraclass correlation coefficient=0.96 and P <0.0001) spatial correlation, with median SCCs ranging from 0.54 to 0.59 throughout the 10 breathing phases. Its phase-averaged FM achieved a median SCC of 0.60, a median Dice similarity coefficient of 0.60 (0.65) for high (low) functional lung volumes, and a correlation of 0.565 (0.646) between the spatially averaged feature values and FEV 1 (FEV 1 /FVC).

CONCLUSIONS

The results provide further insight into the underlying association of specific pulmonary textures with both local ( VNM ) and global (FEV 1 /FVC, FEV 1 ) functions. Further validations of the FM generalizability and the standardization of implementation protocols are warranted before clinically relevant investigations.

Functional imaging for assessing regional lung ventilation in preclinical and clinical research

Dynamic heterogeneity in lung ventilation is an important measure of pulmonary function and may be characteristic of early pulmonary disease. While standard indices like spirometry, body plethysmography, and blood gases have been utilized to assess lung function, they do not provide adequate information on regional ventilatory distribution nor function assessments of ventilation during the respiratory cycle. Emerging technologies such as xenon CT, volumetric CT, functional MRI and X-ray velocimetry can assess regional ventilation using non-invasive radiographic methods that may complement current methods of assessing lung function. As a supplement to current modalities of pulmonary function assessment, functional lung imaging has the potential to identify respiratory disease phenotypes with distinct natural histories. Moreover, these novel technologies may offer an optimal strategy to evaluate the effectiveness of novel therapies and therapies targeting localized small airways disease in preclinical and clinical research. In this review, we aim to discuss the features of functional lung imaging, as well as its potential application and limitations to adoption in research.

Quantification of lung function on CT images based on pulmonary radiomic filtering

PURPOSE

To develop a radiomics filtering technique for characterizing spatial-encoded regional pulmonary ventilation information on lung computed tomography (CT).

METHODS

The lung volume was segmented on 46 CT images, and a 3D sliding window kernel was implemented across the lung volume to capture the spatial-encoded image information. Fifty-three radiomic features were extracted within the kernel, resulting in a fourth-order tensor object. As such, each voxel coordinate of the original lung was represented as a 53-dimensional feature vector, such that radiomic features could be viewed as feature maps within the lungs. To test the technique as a potential pulmonary ventilation biomarker, the radiomic feature maps were compared to paired functional images (Galligas PET or DTPA-SPECT) based on the Spearman correlation (ρ) analysis.

RESULTS

The radiomic feature maps GLRLM-based Run-Length Non-Uniformity and GLCOM-based Sum Average are found to be highly correlated with the functional imaging. The achieved ρ (median [range]) for the two features are 0.46 [0.05, 0.67] and 0.45 [0.21, 0.65] across 46 patients and 2 functional imaging modalities, respectively.

CONCLUSIONS

The results provide evidence that local regions of sparsely encoded heterogeneous lung parenchyma on CT are associated with diminished radiotracer uptake and measured lung ventilation defects on PET/SPECT imaging. These findings demonstrate the potential of radiomics to serve as a complementary tool to the current lung quantification techniques and provide hypothesis-generating data for future studies.

Volumetric multiphase ventilation imaging based on four-dimensional computed tomography for functional lung avoidance radiotherapy

PURPOSE

Current computed tomography (CT)-based lung ventilation imaging (CTVI) techniques derive a static ventilation image without temporal information. This research aims to develop a four-dimensional CT (4DCT)-based multiphase dynamic ventilation imaging framework capable of recovering the entire ventilation process throughout the breathing cycle for functional lung avoidance radiotherapy (FLART).

METHODS

A total of 15 free-breathing thoracic 4DCT scans of lung or esophageal cancer patients were collected from the public datasets. The lung region of each phase image was first delineated, and then the mask-free isotropic total variation image registration algorithm was used to derive the deformation vector fields between the end-expiration (EE) phase and other phases. As a surrogate of ventilation, the voxel-wise local expansion ratio of each phase relative to the EE phase was estimated using the parameterized Integrated Jacobian Formulation method in the EE phase coordinate. Lastly, the dynamic ventilation images were generated by warping these phase-specific local expansion distributions with a same geometry into their respective breathing phases. Quantitative analysis, including interphase Spearman correlation coefficients, voxel-wise, and regional-wise expansion/contraction tracking, were performed to indirectly validate the proposed method.

RESULTS

The proposed method maintains the physiological meaning of ventilation on each phase and enables to recover the dynamic lung ventilation process. The mean interphase Spearman correlations ranged between 0.23 ± 0.20 and 0.93 ± 0.04 and decreased near the EE phase. Only 26.2% (2.59E + 6 out of 9.89E + 6) of lung voxels exhibited the same expansion/contraction pattern as the global lung. Qualitative and quantitative evaluations of the interphase ventilation distribution difference show that ventilation spatiotemporal heterogeneities generally exist during respiration.

CONCLUSIONS

In contrast to conventional CTVI metrics, our method enables to extract additional phase-resolved respiration-correlated information and reflects the generally existed ventilation spatiotemporal heterogeneities. Subsequent studies with quantitative phase-by-phase cross-modality evaluations will further explore its potential to deepen our understanding of lung function and respiration mechanics and also to facilitate more accurate implementation of FLART.

Assessment of Heterogeneity in Lung Structure and Function During Mechanical Ventilation: A Review of Methodologies

The mammalian lung is characterized by heterogeneity in both its structure and function, by incorporating an asymmetric branching airway tree optimized for maintenance of efficient ventilation, perfusion, and gas exchange. Despite potential benefits of naturally occurring heterogeneity in the lungs, there may also be detrimental effects arising from pathologic processes, which may result in deficiencies in gas transport and exchange. Regardless of etiology, pathologic heterogeneity results in the maldistribution of regional ventilation and perfusion, impairments in gas exchange, and increased work of breathing. In extreme situations, heterogeneity may result in respiratory failure, necessitating support with a mechanical ventilator. This review will present a summary of measurement techniques for assessing and quantifying heterogeneity in respiratory system structure and function during mechanical ventilation. These methods have been grouped according to four broad categories: (1) inverse modeling of heterogeneous mechanical function; (2) capnography and washout techniques to measure heterogeneity of gas transport; (3) measurements of heterogeneous deformation on the surface of the lung; and finally (4) imaging techniques used to observe spatially-distributed ventilation or regional deformation. Each technique varies with regard to spatial and temporal resolution, degrees of invasiveness, risks posed to patients, as well as suitability for clinical implementation. Nonetheless, each technique provides a unique perspective on the manifestations and consequences of mechanical heterogeneity in the diseased lung.

Measuring short-term changes in specific ventilation using dynamic specific ventilation imaging

Specific ventilation imaging (SVI) measures the spatial distribution of specific ventilation (SV) in the lung with MRI by using inhaled oxygen as a contrast agent. Because of the inherently low signal to noise ratio in the technique, multiple switches between inspiring air and O₂ are utilized, and the high spatial resolution SV distribution determined as an average over the entire imaging period (~20 minutes). We hypothesized that a trade-off between spatial and temporal resolution could allow imaging at a higher temporal resolution, at the cost of a coarser, yet acceptable, spatial resolution. The appropriate window length and spatial resolution compromise was determined by generating synthetic data with signal- and contrast-to-noise characteristics reflective of that in previously published experimental data, with a known and unchanging distribution of SV, and showed that acceptable results could be obtained in an imaging period of ~7 minutes (80 breaths), with a spatial resolution of ~1cm³. Previously published data were then reanalyzed. The average heterogeneity of the temporally resolved maps of SV were not different to the previous overall analysis, however the temporally resolved maps were less effective at detecting the amount of bronchoconstriction resulting from methacholine administration. The results further indicated that the initial response to inhaled methacholine and subsequent inhalation of albuterol were largely complete within ~22 minutes and ~9 minutes respectively, although there was a tendency for an ongoing developing effect in both cases. These results suggest that it is feasible to use a shortened SVI protocol, with a modest sacrifice in spatial resolution, in order to measure temporally dynamic processes.

Quantitative assessment of pulmonary artery occlusion using lung dynamic perfusion CT

Quantitative measurement of lung perfusion is a promising tool to evaluate lung pathophysiology as well as to assess disease severity and monitor treatment. However, this novel technique has not been adopted clinically due to various technical and physiological challenges; and it is still in the early developmental phase where the correlation between lung pathophysiology and perfusion maps is being explored. The purpose of this research work is to quantify the impact of pulmonary artery occlusion on lung perfusion indices using lung dynamic perfusion CT (DPCT). We performed Lung DPCT in ten anesthetized, mechanically ventilated juvenile pigs (18.6–20.2 kg) with a range of reversible pulmonary artery occlusions (0%, 40–59%, 60–79%, 80–99%, and 100%) created with a balloon catheter. For each arterial occlusion, DPCT data was analyzed using first-pass kinetics to derive blood flow (BF), blood volume (BV) and mean transit time (MTT) perfusion maps. Two radiologists qualitatively assessed perfusion maps for the presence or absence of perfusion defects. Perfusion maps were also analyzed quantitatively using a linear segmented mixed model to determine the thresholds of arterial occlusion associated with perfusion derangement. Inter-observer agreement was assessed using Kappa statistics. Correlation between arterial occlusion and perfusion indices was evaluated using the Spearman-rank correlation coefficient. Our results determined that perfusion defects were detected qualitatively in BF, BV and MTT perfusion maps for occlusions larger than 55%, 80% and 55% respectively. Inter-observer agreement was very good with Kappa scores > 0.92. Quantitative analysis of the perfusion maps determined the arterial occlusion threshold for perfusion defects was 50%, 76% and 44% for BF, BV and MTT respectively. Spearman-rank correlation coefficients between arterial occlusion and normalized perfusion values were strong (− 0.92, − 0.72, and 0.78 for BF, BV and MTT, respectively) and were statically significant (p < 0.01). These findings demonstrate that lung DPCT enables quantification and stratification of pulmonary artery occlusion into three categories: mild, moderate and severe. Severe (occlusion ≥ 80%) alters all perfusion indices; mild (occlusion < 55%) has no detectable effect. Moderate (occlusion 55–80%) impacts BF and MTT but BV is preserved.

Imaging the Injured Lung: Mechanisms of Action and Clinical Use

Acute respiratory distress syndrome (ARDS) consists of acute hypoxemic respiratory failure characterized by massive and heterogeneously distributed loss of lung aeration caused by diffuse inflammation and edema present in interstitial and alveolar spaces. It is defined by consensus criteria, which include diffuse infiltrates on chest imaging-either plain radiography or computed tomography. This review will summarize how imaging sciences can inform modern respiratory management of ARDS and continue to increase the understanding of the acutely injured lung. This review also describes newer imaging methodologies that are likely to inform future clinical decision-making and potentially improve outcome. For each imaging modality, this review systematically describes the underlying principles, technology involved, measurements obtained, insights gained by the technique, emerging approaches, limitations, and future developments. Finally, integrated approaches are considered whereby multimodal imaging may impact management of ARDS.

The VAMPIRE challenge: A multi-institutional validation study of CT ventilation imaging

PURPOSE

CT ventilation imaging (CTVI) is being used to achieve functional avoidance lung cancer radiation therapy in three clinical trials (NCT02528942, NCT02308709, NCT02843568). To address the need for common CTVI validation tools, we have built the Ventilation And Medical Pulmonary Image Registration Evaluation (VAMPIRE) Dataset, and present the results of the first VAMPIRE Challenge to compare relative ventilation distributions between different CTVI algorithms and other established ventilation imaging modalities.

METHODS

The VAMPIRE Dataset includes 50 pairs of 4DCT scans and corresponding clinical or experimental ventilation scans, referred to as reference ventilation images (RefVIs). The dataset includes 25 humans imaged with Galligas 4DPET/CT, 21 humans imaged with DTPA-SPECT, and 4 sheep imaged with Xenon-CT. For the VAMPIRE Challenge, 16 subjects were allocated to a training group (with RefVI provided) and 34 subjects were allocated to a validation group (with RefVI blinded). Seven research groups downloaded the Challenge dataset and uploaded CTVIs based on deformable image registration (DIR) between the 4DCT inhale/exhale phases. Participants used DIR methods broadly classified into B-splines, Free-form, Diffeomorphisms, or Biomechanical modeling, with CT ventilation metrics based on the DIR evaluation of volume change, Hounsfield Unit change, or various hybrid approaches. All CTVIs were evaluated against the corresponding RefVI using the voxel-wise Spearman coefficient r S , and Dice similarity coefficients evaluated for low function lung ( DSC low ) and high function lung ( DSC high ).

RESULTS

A total of 37 unique combinations of DIR method and CT ventilation metric were either submitted by participants directly or derived from participant-submitted DIR motion fields using the in-house software, VESPIR. The r S and DSC results reveal a high degree of inter-algorithm and intersubject variability among the validation subjects, with algorithm rankings changing by up to ten positions depending on the choice of evaluation metric. The algorithm with the highest overall cross-modality correlations used a biomechanical model-based DIR with a hybrid ventilation metric, achieving a median (range) of 0.49 (0.27-0.73) for r S , 0.52 (0.36-0.67) for DSC low , and 0.45 (0.28-0.62) for DSC high . All other algorithms exhibited at least one negative r S value, and/or one DSC value less than 0.5.

CONCLUSIONS

The VAMPIRE Challenge results demonstrate that the cross-modality correlation between CTVIs and the RefVIs varies not only with the choice of CTVI algorithm but also with the choice of RefVI modality, imaging subject, and the evaluation metric used to compare relative ventilation distributions. This variability may arise from the fact that each of the different CTVI algorithms and RefVI modalities provides a distinct physiologic measurement. Ultimately this variability, coupled with the lack of a "gold standard," highlights the ongoing importance of further validation studies before CTVI can be widely translated from academic centers to the clinic. It is hoped that the information gleaned from the VAMPIRE Challenge can help inform future validation efforts.

The spatial pattern of methacholine bronchoconstriction recurs when supine, independently of posture during provocation, but does not recur between postures

The location of lung regions with compromised ventilation (often called ventilation defects) during a bronchoconstriction event may be influenced by posture. We aimed to determine the effect of prone vs. supine posture on the spatial pattern of methacholine-induced bronchoconstriction in six healthy adults (ages 21-41, three females) using specific ventilation imaging. Three postural conditions were chosen to assign the effect of posture to the drug administration and/or imaging phase of the experiment - supine methacholine administration followed by supine imaging, prone methacholine administration followed by supine imaging, and prone methacholine administration followed by prone imaging. The two conditions in which imaging was performed supine had similar spatial patterns of bronchoconstriction despite a change in posture during methacholine administration; the odds ratio for recurrent constriction was mean (SD) = 7.4 (3.9). Conversely, dissimilar spatial patterns of bronchoconstriction emerged when posture during imaging was changed; the odds ratio for recurrent constriction between the prone methacholine/supine imaging condition and the prone methacholine/prone imaging condition was 1.2 (0.9). Logistic regression showed that height above the dependent lung border was a significant negative predictor of constriction in the two supine imaging conditions (p<0.001 for each), but not in the prone imaging condition (p=0.20). These results show that the spatial pattern of methacholine bronchoconstriction is recurrent in the supine posture, regardless of whether methacholine is given prone or supine, but that prone posture during imaging eliminates that recurrent pattern and reduces its dependence on gravitational height.

Assessment of pulmonary structure-function relationships in young children and adolescents with cystic fibrosis by multivolume proton-MRI and CT

BACKGROUND

Lung disease is the most frequent cause of morbidity and mortality in patients with cystic fibrosis (CF), and there is a shortage of sensitive biomarkers able to regionally monitor disease progression and to assess early responses to therapy.

PURPOSE

To determine the feasibility of noncontrast-enhanced multivolume MRI, which assesses intensity changes between expiratory and inspiratory breath-hold images, to detect and quantify regional ventilation abnormalities in CF lung disease, with a focus on the structure-function relationship.

STUDY TYPE

Retrospective.

POPULATION

Twenty-nine subjects, including healthy young children (n = 9, 7-37 months), healthy adolescents (n = 4, 14-22 years), young children with CF lung disease (n = 10, 7-47 months), and adolescents with CF lung disease (n = 6, 8-18 years) were studied.

FIELD STRENGTH/SEQUENCE

3D spoiled gradient-recalled sequence at 1.5T.

ASSESSMENT

Subjects were scanned during breath-hold at functional residual capacity (FRC) and total lung capacity (TLC) through noncontrast-enhanced MRI and CT. Expiratory-inspiratory differences in MR signal-intensity (Δ1 H-MRI) and CT-density (ΔHU) were computed to estimate regional ventilation. MR and CT images were also evaluated using a CF-specific scoring system.

STATISTICAL TESTS

Quadratic regression, Spearman's correlation, one-way analysis of variance (ANOVA).

RESULTS

Δ1 H-MRI maps were sensitive to ventilation heterogeneity related to gravity dependence in healthy lung and to ventilation impairment in CF lung disease. A high correlation was found between MRI and CT ventilation maps (R2  = 0.79, P < 0.001). Globally, Δ1 H-MRI and ΔHU decrease with increasing morphological score (respectively, R2  = 0.56, P < 0.001 and R2  = 0.31, P < 0.001). Locally, Δ1 H-MRI was higher in healthy regions (median 15%) compared to regions with bronchiectasis, air trapping, consolidation, and to segments fed by airways with bronchial wall thickening (P < 0.001).

DATA CONCLUSION

Multivolume noncontrast-enhanced MRI, as a nonionizing imaging modality that can be used on nearly any MRI scanner without specialized equipment or gaseous tracers, may be particularly valuable in CF care, providing a new imaging biomarker to detect early alterations in regional lung structure-function.

LEVEL OF EVIDENCE

3 Technical Efficacy: Stage 3 J. MAGN. RESON. IMAGING 2018;48:531-542.

Spatial persistence of reduced specific ventilation following methacholine challenge in the healthy human lung

Specific ventilation imaging was used to identify regions of the healthy lung (6 supine subjects, ages 21-41 yr, 3 men) that experienced a fall in specific ventilation following inhalation of methacholine. This test was repeated 1 wk later and 3 mo later to test for spatial recurrence. Our data showed that 53% confidence interval (CI; 46%, 59%) of volume elements that constricted during one methacholine challenge did so again in another and that this quantity did not vary with time; 46% CI (28%, 64%) recurred 1 wk later, and 56% CI (51%, 61%) recurred 3 mo later. Previous constriction was a strong predictor for future constriction. Volume elements that constricted during one challenge were 7.7 CI (5.2, 10.2) times more likely than nonconstricted elements to constrict in a second challenge, regardless of whether the second episode was 1 wk [7.7 CI (2.9, 12.4)] or 3 mo [7.7 CI (4.6, 10.8)] later. Furthermore, posterior lung elements were more likely to constrict following methacholine than anterior lung elements (volume fraction 0.43 ± 0.22 posterior vs. 0.10 ± 0.03 anterior; P = 0.005), and basal elements that constricted were more likely than their apical counterparts to do so persistently through all three trials (volume fraction 0.14 ± 0.04 basal vs. 0.04 ± 0.04 apical; P = 0.003). Taken together, this evidence suggests a physiological predisposition toward constriction in some lung elements, especially those located in the posterior and basal lung when the subject is supine. NEW & NOTEWORTHY The spatial pattern of bronchoconstriction following methacholine is persistent over time in healthy individuals, in whom chronic inflammation and airway remodeling are assumed to be absent. This suggests that regional lung inflation and airway structure may play dominant roles in determining the spatial pattern of methacholine bronchoconstriction.

Inter-airway structural heterogeneity interacts with dynamic heterogeneity to determine lung function and flow patterns in both asthmatic and control simulated lungs

Asthma is a disease involving both airway remodelling (e.g. thickening of the airway wall) and acute, reversible airway narrowing driven by airway smooth muscle contraction. Both of these processes are known to be heterogeneous, and in this study we consider a new theoretical model which considers the interactions of both mechanisms: structural heterogeneity (variation in airway remodelling) and dynamic heterogeneity (emergent variation in airway narrowing and flow). By integrating both types of inter-airway heterogeneity in a full human lung geometry, we are able to draw several insights regarding the mechanisms underlying observed ventilation heterogeneity. We show that: (1) bimodal ventilation distributions are driven by paradoxical contraction/dilation patterns for airways of all sizes; (2) structural heterogeneity differences between asthmatic and control subjects significantly influences resulting lung function, and observed ventilation heterogeneity patterns; and (3) individual airway dilation probabilities are uncorrelated with prior airway remodelling of that airway.

Effects of neurally adjusted ventilatory assist on air distribution and dead space in patients with acute exacerbation of chronic obstructive pulmonary disease

Background

Neurally adjusted ventilatory assist (NAVA) could improve patient-ventilator interaction; its effects on ventilation distribution and dead space are still unknown. The aim of this study was to evaluate the effects of varying levels of assist during NAVA and pressure support ventilation (PSV) on ventilation distribution and dead space in patients with acute exacerbation of chronic obstructive pulmonary disease (AECOPD).

Methods

Fifteen mechanically ventilated patients with AECOPD were included in the study. The initial PSV levels were set to 10 cmH₂O for 10 min. Thereafter, the ventilator mode was changed to NAVA for another 10 min with the same electrical activity of the diaphragm as during PSV. Furthermore, the ventilation mode was switched between PSV and NAVA every 10 min in the following order: PSV 5 cmH₂O; NAVA 50%; PSV 15 cmH₂O; and NAVA 150% (relative to the initial NAVA support level). Ventilation distribution in the lung was evaluated in percentages in regions of interest (ROI) of four anteroposterior segments of equal height (ROI1 to ROI4 represents ventral, mid-ventral, mid-dorsal, and dorsal, respectively). Blood gases, ventilation distribution (electrical impedance tomography), diaphragm activity (B-mode ultrasonography), and dead space fraction (PeCO₂ and PaCO₂) were measured.

Results

The trigger and cycle delays were lower during NAVA than during PSV. The work of trigger was significantly lower during NAVA compared to PSV. The diaphragm activities based on ultrasonography were higher during NAVA compared to the same support level during PSV. The ventilation distribution in ROI4 increased significantly (P < 0.05) during NAVA compared to PSV (except for a support level of 50%). Similar results were found in ROI3 + 4. NAVA reduced dead space fraction compared to the corresponding support level of PSV.

Conclusions

NAVA was superior to PSV in AECOPD for increasing ventilation distribution in ROI4 and reducing dead space.

Trial registration

Clinicaltrials.gov, NCT02289573. Registered on 12 November 2014.

Using Hyperpolarized 129Xe MRI to Quantify the Pulmonary Ventilation Distribution

RATIONALE AND OBJECTIVES

Ventilation heterogeneity is impossible to detect with spirometry. Alternatively, pulmonary ventilation can be imaged three-dimensionally using inhaled ¹²⁹Xe magnetic resonance imaging (MRI). To date, such images have been quantified primarily based on ventilation defects. Here, we introduce a robust means to transform ¹²⁹Xe MRI scans such that the underlying ventilation distribution and its heterogeneity can be quantified.

MATERIALS AND METHODS

Quantitative ¹²⁹Xe ventilation MRI was conducted in 12 younger (24.7 ± 5.2 years) and 10 older (62.2 ± 7.2 years) healthy individuals, as well as in 9 younger (25.9 ± 6.4 yrs) and 10 older (63.2 ± 6.1 years) asthmatics. The younger healthy population was used to establish a reference ventilation distribution and thresholds for six intensity bins. These bins were used to display and quantify the ventilation defect region (VDR), the low ventilation region (LVR), and the high ventilation region (HVR).

RESULTS

The ventilation distribution in young subjects was roughly Gaussian with a mean and standard deviation of 0.52 ± 0.18, resulting in VDR = 2.1 ± 1.3%, LVR = 15.6 ± 5.4%, and HVR = 17.4 ± 3.1%. Older healthy volunteers exhibited a significantly right-skewed distribution (0.46 ± 0.20, P = 0.034), resulting in significantly increased VDR (7.0 ± 4.8%, P = 0.008) and LVR (24.5 ± 11.5%, P = 0.025). In the asthmatics, VDR and LVR increased in the older population, and HVR was significantly reduced (13.5 ± 4.6% vs 18.9 ± 4.5%, P = 0.009). Quantitative ¹²⁹Xe MRI also revealed altered ventilation heterogeneity in response to albuterol in two asthmatics with normal spirometry.

CONCLUSIONS

Quantitative ¹²⁹Xe MRI provides a robust and objective means to display and quantify the pulmonary ventilation distribution, even in subjects who have airway function impairment not appreciated by spirometry.

Clustered ventilation defects and bilinear respiratory reactance in asthma

Imaging studies of asthmatics typically reveal clustered ventilation patterns, rather than homogeneous ventilation; furthermore, the variation of these clusters suggests that the causes are at least partially dynamic, rather than structural. Theoretical studies have indicated dynamic mechanisms by which homogeneous ventilation solutions lose stability and clustered solutions emerge. At the same time, it has been demonstrated experimentally that respiratory reactance characteristically has a bilinear relationship with lung volume, and that changes to this relationship are indicative of various aspects of disease progression and control. Moreover, the transition point in the bilinear reactance relationship is thought to relate to reopening/recruitment of airway units, and thus may be connected to the bifurcation via which clustered ventilation solutions emerge. In order to investigate this possibility we develop a new model, including both airway-airway coupling and airway-parenchymal coupling, which exhibits both clustered ventilation defects and also a bilinear reactance relationship. Studying this model reveals that (1) the reactance breakpoint is not coincident with the bifurcation; (2) numerous changes to underlying behaviour can alter the reactance breakpoint in ways which mimic the experimental data; and (3) the location of ventilation defects can be a combination of both structural and dynamic factors.

Systems-level airway models of bronchoconstriction

Understanding lung and airway behavior presents a number of challenges, both experimental and theoretical, but the potential rewards are great in terms of both potential treatments for disease and interesting biophysical phenomena. This presents an opportunity for modeling to contribute to greater understanding, and here, we focus on modeling efforts that work toward understanding the behavior of airways in vivo, with an emphasis on asthma. We look particularly at those models that address not just isolated airways but many of the important ways in which airways are coupled both with each other and with other structures. This includes both interesting phenomena involving the airways and the layer of airway smooth muscle that surrounds them, and also the emergence of spatial ventilation patterns via dynamic airway interaction. WIREs Syst Biol Med 2016, 8:459-467. doi: 10.1002/wsbm.1349 For further resources related to this article, please visit the WIREs website.

Mild loss of lung aeration augments stretch in healthy lung regions

Inspiratory stretch by mechanical ventilation worsens lung injury. However, it is not clear whether and how the ventilator damages lungs in the absence of preexisting injury. We hypothesized that subtle loss of lung aeration during general anesthesia regionally augments ventilation and distension of ventilated air spaces. In eight supine anesthetized and intubated rats, hyperpolarized gas MRI was performed after a recruitment maneuver following 1 h of volume-controlled ventilation with zero positive end-expiratory pressure (ZEEP), FiO2 0.5, and tidal volume 10 ml/kg, and after a second recruitment maneuver. Regional fractional ventilation (FV), apparent diffusion coefficient (ADC) of (3)He (a measurement of ventilated peripheral air space dimensions), and gas volume were measured in lung quadrants of ventral and dorsal regions of the lungs. In six additional rats, computed tomography (CT) images were obtained at each time point. Ventilation with ZEEP decreased total lung gas volume and increased both FV and ADC in all studied regions. Increases in FV were more evident in the dorsal slices. In each lung quadrant, higher ADC was predicted by lower gas volume and by increased mean values (and heterogeneity) of FV distribution. CT scans documented 10% loss of whole-lung aeration and increased density in the dorsal lung, but no macroscopic atelectasis. Loss of pulmonary gas at ZEEP increased fractional ventilation and inspiratory dimensions of ventilated peripheral air spaces. Such regional changes could help explain a propensity for mechanical ventilation to contribute to lung injury in previously uninjured lungs.

Globally optimal co-segmentation of three-dimensional pulmonary ¹H and hyperpolarized ³He MRI with spatial consistence prior

Pulmonary imaging using hyperpolarized (3)He/(129)Xe gas is emerging as a new way to understand the regional nature of pulmonary ventilation abnormalities in obstructive lung diseases. However, the quantitative information derived is completely dependent on robust methods to segment both functional and structural/anatomical data. Here, we propose an approach to jointly segment the lung cavity from (1)H and (3)He pulmonary magnetic resonance images (MRI) by constraining the spatial consistency of the two segmentation regions, which simultaneously employs the image features from both modalities. We formulated the proposed co-segmentation problem as a coupled continuous min-cut model and showed that this combinatorial optimization problem can be solved globally and exactly by means of convex relaxation. In particular, we introduced a dual coupled continuous max-flow model to study the convex relaxed coupled continuous min-cut model under a primal and dual perspective. This gave rise to an efficient duality-based convex optimization algorithm. We implemented the proposed algorithm in parallel using general-purpose programming on graphics processing unit (GPGPU), which substantially increased its computational efficiency. Our experiments explored a clinical dataset of 25 subjects with chronic obstructive pulmonary disease (COPD) across a wide range of disease severity. The results showed that the proposed co-segmentation approach yielded superior performance compared to single-channel image segmentation in terms of precision, accuracy and robustness.

Hyperpolarized gas diffusion MRI for the study of atelectasis and acute respiratory distress syndrome

Considerable uncertainty remains about the best ventilator strategies for the mitigation of atelectasis and associated airspace stretch in patients with acute respiratory distress syndrome (ARDS). In addition to several immediate physiological effects, atelectasis increases the risk of ventilator-associated lung injury, which has been shown to significantly worsen ARDS outcomes. A number of lung imaging techniques have made substantial headway in clarifying the mechanisms of atelectasis. This paper reviews the contributions of computed tomography, positron emission tomography, and conventional MRI to understanding this phenomenon. In doing so, it also reveals several important shortcomings inherent to each of these approaches. Once these shortcomings have been made apparent, we describe how hyperpolarized (HP) gas MRI--a technique that is uniquely able to assess responses to mechanical ventilation and lung injury in peripheral airspaces--is poised to fill several of these knowledge gaps. The HP-MRI-derived apparent diffusion coefficient (ADC) quantifies the restriction of (3) He diffusion by peripheral airspaces, thereby obtaining pulmonary structural information at an extremely small scale. Lastly, this paper reports the results of a series of experiments that measured ADC in mechanically ventilated rats in order to investigate (i) the effect of atelectasis on ventilated airspaces, (ii) the relationship between positive end-expiratory pressure (PEEP), hysteresis, and the dimensions of peripheral airspaces, and (iii) the ability of PEEP and surfactant to reduce airspace dimensions after lung injury. An increase in ADC was found to be a marker of atelectasis-induced overdistension. With recruitment, higher airway pressures were shown to reduce stretch rather than worsen it. Moving forward, HP MRI has significant potential to shed further light on the atelectatic processes that occur during mechanical ventilation.

Validating the distribution of specific ventilation in healthy humans measured using proton MR imaging

Specific ventilation imaging (SVI) uses proton MRI to quantitatively map the distribution of specific ventilation (SV) in the human lung, using inhaled oxygen as a contrast agent. To validate this recent technique, we compared the quantitative measures of heterogeneity of the SV distribution in a 15-mm sagittal slice of lung obtained in 10 healthy supine subjects, (age 37 ± 10 yr, forced expiratory volume in 1 s 97 ± 7% predicted) using SVI to those obtained in the whole lung from multiple-breath nitrogen washout (MBW). Using the analysis of Lewis et al. (Lewis SM, Evans JW, Jalowayski AA. J App Physiol 44: 416-423, 1978), the most likely distribution of SV from the MBW data was computed and compared with the distribution of SV obtained from SVI, after normalizing for the difference in tidal volume. The average SV was 0.30 ± 0.10 MBW, compared with 0.36 ± 0.10 SVI (P = 0.01). The width of the distribution, a measure of the heterogeneity, obtained using both methods was comparable: 0.51 ± 0.06 and 0.47 ± 0.08 in MBW and SVI, respectively (P = 0.15). The MBW estimated width of the SV distribution was 0.05 (10.7%) higher than that estimated using SVI, and smaller than the intertest variability of the MBW estimation [inter-MBW (SD) for the width of the SV distribution was 0.08 (15.8)%]. To assess reliability, SVI was performed twice on 13 subjects showing small differences between measurements of SV heterogeneity (typical error 0.05, 12%). In conclusion, quantitative estimations of SV heterogeneity from SVI are reliable and similar to those obtained using MBW, with SVI providing spatial information that is absent in MBW.

4D CT lung ventilation images are affected by the 4D CT sorting method

PURPOSE

Four-dimensional (4D) computed tomography (CT) ventilation imaging is a novel promising technique for lung functional imaging. The current standard 4D CT technique using phase-based sorting frequently results in artifacts, which may deteriorate the accuracy of ventilation imaging. The purpose of this study was to quantify the variability of 4D CT ventilation imaging due to 4D CT sorting.

METHODS

4D CT image sets from nine lung cancer patients were each sorted by the phase-based method and anatomic similarity-based method, designed to reduce artifacts, with corresponding ventilation images created for each method. Artifacts in the resulting 4D CT images were quantified with the artifact score which was defined based on the difference between the normalized cross correlation for CT slices within a CT data segment and that for CT slices bordering the interface between adjacent CT data segments. The ventilation variation was quantified using voxel-based Spearman rank correlation coefficients for all lung voxels, and Dice similarity coefficients (DSC) for the spatial overlap of low-functional lung volumes. Furthermore, the correlations with matching single-photon emission CT (SPECT) ventilation images (assumed ground truth) were evaluated for three patients to investigate which sorting method provides higher physiologic accuracy.

RESULTS

Anatomic similarity-based sorting reduced 4D CT artifacts compared to phase-based sorting (artifact score, 0.45 ± 0.14 vs 0.58 ± 0.24, p = 0.10 at peak-exhale; 0.63 ± 0.19 vs 0.71 ± 0.31, p = 0.25 at peak-inhale). The voxel-based correlation between the two ventilation images was 0.69 ± 0.26 on average, ranging from 0.03 to 0.85. The DSC was 0.71 ± 0.13 on average. Anatomic similarity-based sorting yielded significantly fewer lung voxels with paradoxical negative ventilation values than phase-based sorting (5.0 ± 2.6% vs 9.7 ± 8.4%, p = 0.05), and improved the correlation with SPECT ventilation regionally.

CONCLUSIONS

The variability of 4D CT ventilation imaging due to 4D CT sorting was moderate overall and substantial in some cases, suggesting that 4D CT artifacts are an important source of variations in 4D CT ventilation imaging. Reduction of 4D CT artifacts provided more physiologically convincing and accurate ventilation estimates. Further studies are needed to confirm this result.

Dynamic multiscale boundary conditions for 4D CT of healthy and emphysematous rats

Changes in the shape of the lung during breathing determine the movement of airways and alveoli, and thus impact airflow dynamics. Modeling airflow dynamics in health and disease is a key goal for predictive multiscale models of respiration. Past efforts to model changes in lung shape during breathing have measured shape at multiple breath-holds. However, breath-holds do not capture hysteretic differences between inspiration and expiration resulting from the additional energy required for inspiration. Alternatively, imaging dynamically--without breath-holds--allows measurement of hysteretic differences. In this study, we acquire multiple micro-CT images per breath (4DCT) in live rats, and from these images we develop, for the first time, dynamic volume maps. These maps show changes in local volume across the entire lung throughout the breathing cycle and accurately predict the global pressure-volume (PV) hysteresis. Male Sprague-Dawley rats were given either a full- or partial-lung dose of elastase or saline as a control. After three weeks, 4DCT images of the mechanically ventilated rats under anesthesia were acquired dynamically over the breathing cycle (11 time points, ≤100 ms temporal resolution, 8 cmH2O peak pressure). Non-rigid image registration was applied to determine the deformation gradient--a numerical description of changes to lung shape--at each time point. The registration accuracy was evaluated by landmark identification. Of 67 landmarks, one was determined misregistered by all three observers, and 11 were determined misregistered by two observers. Volume change maps were calculated on a voxel-by-voxel basis at all time points using both the Jacobian of the deformation gradient and the inhaled air fraction. The calculated lung PV hysteresis agrees with pressure-volume curves measured by the ventilator. Volume maps in diseased rats show increased compliance and ventilation heterogeneity. Future predictive multiscale models of rodent respiration may leverage such volume maps as boundary conditions.