Copyright ©ERS Journals Ltd 2007 Relationship between airway narrowing, patchy ventilation and lung mechanics in asthmatics1 Depts of Anesthesia and Critical Care, 3 Medicine (Pulmonary and Critical Care Unit and Clinical Epidemiology Unit), Massachusetts General Hospital, 2 Dept of Biomedical Engineering, Boston University, and 4 Harvard Medical School, Boston, MA, USA. CORRESPONDENCE: J. G. Venegas, Dept of Anesthesia and Critical Care (CLN-237F), Massachusetts General Hospital, Boston, 02114 MA, USA. E-mail: jvenegas{at}partners.org Keywords: Airway obstruction, constriction, elastance, heterogeneous, oscillatory, resistance
Received: August 29, 2006
Bronchoconstriction in asthma results in patchy ventilation forming ventilation defects (VDefs). Patchy ventilation is clinically important because it affects obstructive symptoms and impairs both gas exchange and the distribution of inhaled medications. The current study combined functional imaging, oscillatory mechanics and theoretical modelling to test whether the degrees of constriction of airways feeding those units outside VDefs were related to the extent of VDefs in bronchoconstricted asthmatic subjects. Positron emission tomography was used to quantify the regional distribution of ventilation and oscillatory mechanics were measured in asthmatic subjects before and after bronchoconstriction. For each subject, ventilation data was mapped into an anatomically based lung model that was used to evaluate whether airway constriction patterns, consistent with the imaging data, were capable of matching the measured changes in airflow obstruction. The degree and heterogeneity of constriction of the airways feeding alveolar units outside VDefs was similar among the subjects studied despite large inter-subject variability in airflow obstruction and the extent of the ventilation defects. Analysis of the data amongst the subjects showed an inverse relationship between the reduction in mean airway conductance, measured in the breathing frequency range during bronchoconstriction, and the fraction of lung involved in ventilation defects. The current data supports the concept that patchy ventilation is an expression of the integrated system and not just the sum of independent responses of individual airways. In asthma, ventilation distribution and lung mechanics are extensively affected by heterogeneous airway constriction. Detailed imaging studies in the bronchoconstricted lung using magnetic resonance imaging 1, 2 and positron emission tomography (PET) 3, 4 have demonstrated the presence of patchy ventilation characterised by large ventilation defects (VDefs). The extent of these VDefs initially led investigators to speculate that they were likely to be caused by severe constriction of central airways feeding the constricted regions. However, experimental evidence and theoretical modelling suggest that most of the functional impairment in asthma is caused by heterogeneous constriction of peripheral airways 5. This mechanical evidence is thus inconsistent with the concept that constriction of large airways could be exclusively responsible for VDefs. However, the current authors have recently demonstrated that even a symmetric airway tree undergoing uniform smooth muscle activation can develop highly heterogeneous ventilation as a result of severely constricted peripheral airways clustered in VDefs similar to those observed with imaging 3, 4. Furthermore, this model predicted that, for constrictive stimuli beyond a critical level, airway narrowing did not progress throughout the entire airway tree but was localised, increasing the extent of the lung covered by VDefs. The latter prediction cannot be directly tested experimentally because limited spatial resolution of current imaging techniques does not allow visualisation of all the affected airways in a breathing lung. However, in a recent study the present authors showed that combining PET imaging with measurements and analysis of the frequency dependence of dynamic resistance (RL) and elastance (EL), could be used to test hypothetical airway constriction patterns imposed on an anatomically consistent model of the airway tree 6. These image-guided modelling tools are used in the present study to evaluate the relationship between obstructive changes of the airway tree and relative extent (size) of the VDefs in bronchoconstricted asthmatics. For this purpose, VDefs defined from PET images of each subject were mapped into the corresponding acini of an anatomically consistent three-dimensional (3-D) model of the airway tree 7. Bronchoconstriction was simulated by systematically varying the lumen of specific airways of the tree. This allowed the definition of spatial patterns of airway obstruction that best approached the changes in oscillatory mechanics after bronchoprovocation in each of the subjects. The spatial patterns were evaluated to test whether constriction of airways feeding those units outside the VDefs could be associated with the differences, observed amongst subjects, in the extent of the lung covered by VDefs.
A group of 10 mild-to-moderate asthmatic subjects were studied before and after bronchoconstriction with a protocol approved by the Massachusetts General Hospital review board (Boston, MA, USA; table 1 20% of the highest residual tracer concentration. A 20% threshold for residual tracer concentration was selected because at thresholds between 20 and 30%, the change in size of VDefs was found to be rather insensitive to the threshold level selected. In addition, a constant 20% threshold for all subjects yielded the largest possible VDefs size within that range of thresholds.
Oscillatory mechanics were assessed by delivering, through a mouthpiece, seven superimposed sinusoidal waveforms with frequencies ranging 0.158 Hz. Transpulmonary pressure (Ptp), measured as the difference between airway opening pressure and oesophageal pressure, and flow (Qao) at the airway opening measured with a pneuomotachograph, were electronically acquired and monitored for 23 min. Signal processing on Qao and Ptp was used to estimate values of RL and EL as a function of frequency 8.
Computational modelling
For bronchoconstriction the level of airway narrowing applied to the model was defined based on whether or not the individual airway fed VDefs. First, a fixed fraction (40%) of the terminal bronchioles that led to VDefs were reduced in diameter to a fixed fraction (20%) of that at baseline while the rest of the airway tree was left at the baseline size (pattern A). Next, in addition to the constriction conditions in pattern A, all bronchi exclusively feeding the VDefs were constricted to a fixed fraction (30%) of baseline diameter (pattern B). The rationale for choosing the specific degrees of airway constriction in patterns A and B comes from a previous modelling study 5. Finally, in addition to the constriction conditions in pattern B, heterogeneous constriction was imposed on the airways that did not lead to VDefs, by random sampling of Gaussian distributions covering a grid of means and SDs (pattern C).
The simulated lung mechanics best approaching those measured, for a given subject was taken as those yielding the lowest simulation index (SI) defined as:
where Re and Im are the real and imaginary part of lung impedance, respectively, i denotes each of the nf frequencies studied, d refers to experimental and m to modelled data. Therefore, the pair of mean and SD yielding a minimum value of SI (fig. 2
Statistical analysis Comparisons were performed using paired t-tests with significance level of p<0.05. Linear regressions were conducted for baseline-normalised airway diameters (D/Do) and baseline-normalised airway conductance (G/Go) versus the extent of the model's VDefs. Linear regressions were also conducted for baseline-normalised airway conductance at 0.15 Hz during bronchoconstriction (GMch/GBase) versus the fraction of imaged lung volume covered by VDefs.
The sets of airway diameters of the model that simulated measured lung mechanics at baseline for all subjects, involved a substantial constriction with a small level of heterogeneity (mean±SD 0.42±0.08). These adjusted airway diameters defined baseline conditions for each subject.
During bronchoconstriction all subjects showed a considerable increase in the frequency dependence and mean levels of RL and EL compared with baseline (fig. 3
It is therefore clear that, in addition to severe constriction of terminal bronchioles within VDefs and the airways leading to them, a certain degree of heterogeneous constriction needs to be applied to the rest of the airway tree to explain the changes in oscillatory mechanics measured in the subjects during bronchoconstriction. There was substantial variability in the extent of VDefs amongst subjects, ranging 1980% of the imaged lung (fig. 1
where TB is the fraction of terminal bronchioles encompassed within the VDefs (fig. 6a
The assumption of an average G/Go for the model as proportional to D/Do to the fourth power (D/Do4), also yielded a negative correlation (r = 0.91) between G/Go and TB which can be approximated by linear regression as:
This relationship (fig. 6
where RMch and RBase are oscillatory lung resistance measured at 0.15 Hz during bronchoconstriction and at baseline, respectively, and VDef is the fraction of voxels of imaged lung encompassed within the VDefs (fig. 6b
The main findings of the current study were as follows. 1) Severe constriction of terminal bronchioles and airways leading to VDefs was not enough to explain the mechanical obstruction measured in bronchoconstriction. 2) A significant degree of heterogeneous constriction of the rest of the airway tree was necessary to account for the measured changes in oscillatory mechanics of each subject. 3) Inter-subject differences in the average degree of constriction of airways not leading to VDefs were not correlated with the extent of VDefs. Before discussing these findings, it is important to point out methodological assumptions and limitations of the current methodology. Baseline airway diameters of the model were initially adjusted to match measured baseline RL and EL. Initially, a previously published model of airway dimensions 13 derived from a generic normal lung in the upright position for baseline conditions was used, as a result, the airway diameters had to be adjusted to account for the airway size differences between that generic model and the lung of each asthmatic subject studied. The actual amount of adjustment is also justified because the subjects were studied when lying down, a condition expected to reduce lung volume at FRC. Thus, it can be expected that to match the baseline oscillatory mechanics, the generic model required a significant mean level of constriction (mean 0.42). Interestingly, the SD of the adjustment was low (0.08) and thus represented a relatively homogenous constriction pattern. Regions of VDefs, identified from PET images, were mapped on to an anatomically consistent 3-D model of the lungs. In this model, severe constriction of the terminal bronchioles of as much as 20% of the baseline diameter was defined because previous modelling analysis 5 showed that at least that much constriction was required to simulate the elevation of elastance at normal breathing frequencies observed in asthmatics during bronchoconstriction. Terminal bronchioles in the model were not fully obstructed because analysis of the PET images showed that units within VDefs ventilated slowly but were usually not totally unventilated. Preliminary analysis of the data also showed that such a reduction in diameter resulted in regional ventilation of VDefs, which was consistent with that observed in poorly ventilated units with PET 3, 4, 14.
Previous analysis of the PET images has shown that voxels within the VDefs were not uniformly hypoventilated but instead had a mixture of ventilating units and hypoventilating units 4, 15. To account for this it was assumed that 40% of terminal bronchioles within the VDefs of the model were severely constricted while the diameters of the rest were kept at baseline levels. No significant differences in the impedance spectrum predicted by the model were observed when as much as 100% of the terminal bronchioles within the VDefs of the model were severely constricted. Diameters of the airways feeding VDefs were reduced to 30% of baseline consistent with previous modelling results 6. The simulations showed that exclusive constriction to
Integration of ventilation imaging and lung mechanics
Both the mean and SD of constriction in the best match simulations were similar amongst subjects and uncorrelated with the extent of their respective VDefs (fig. 5
The best match simulations yielded a negative correlation between the average baseline-normalised airway diameters and the extent of the VDefs among subjects. This finding implies that part of the differences in the global mechanical response among subjects (fig. 6 The inverse relationship between mean baseline-normalised diameter of all airways of the model and the fraction of terminal bronchi involved in the VDefs was consistent with the equivalent relationship for experimental data obtained between the baseline-normalised airway conductance, measured at 0.15 Hz, and the fraction of the imaged lung involved in VDefs.
Both modelling and experimental correlations can be expected given the concentration of severely constricted terminal airways in the VDefs from the modelled and imaged lungs. A less expected finding was that both modelling and experimental relationships (fig. 6
G calculated as the inverse of RL is usually regarded as a surrogate for changes of the airway lumen. Both G and RL are typically measured at relatively high frequencies (58 Hz) because of the minimal contribution of tissue resistance and parallel heterogeneity exhibited at these frequencies as previously depicted by network models of lung mechanics using rigid airways. When using models with rigid airways, heterogeneity is in fact more of a contributor to the low frequency data (0.15 Hz) compared with the high frequency data (8 Hz) 20. To test whether in realistic conditions the presence of airway compliance could affect that conclusion, G/Go was also estimated at 8 Hz assuming rigid walls with the constriction conditions found to best match the measured changes in mechanics during bronchoconstriction for each subject. If G/Go at 8 Hz for the rigid airways is taken as the best indicator of changes in average airway lumen during bronchoconstriction, then the deviation between this and G/Go derived with compliant airways, at any frequency, could be considered as the error in assessing the change in airway lumen under realistic conditions at that frequency. It was found that, in the models with compliant airways, G/Go at 8 Hz (mean±SD 0.15±0.13 L1·cmH20·s1) deviated three times more than G/Go at 0.15 Hz (0.053±0.0408 L1·cmH20·s1) relative to G/Go at 8 Hz obtained assuming rigid airways: this was true in the models of nine out of 10 subjects. It is concluded, due to the presence of airway compliance, that the change in global G measured at a breathing frequency of 0.15 Hz can better approximate the change in airway calibre than when measured at 8 Hz. This may explain the better correlation obtained in figure 6 The suggestion that a critical reduction of conductance takes place before a patchy pattern of ventilation distribution develops is consistent with the behaviour of the network model of the airway tree 4. This network model is predictive of the formation of the VDefs and change in airway diameter as a result of dynamic interactions among different airways of the tree unlike the 3-D airway tree model, which is not predictive of the formation of the VDefs, and requires airway diameters and the location of ventilation defects as inputs. More specifically, the network model incorporated the effects of tidal stretch on smooth muscle 21, 22, the interdependence between airway wall and parenchymal forces 23, 24 and the interdependent behaviour between airways within a branching tracheobronchial tree 25. This model was assumed to be virtually symmetrical and with uniform structural and functional properties to demonstrate that heterogeneity in those properties was not a requirement to develop the heterogeneity in constriction observed in asthma. Indeed, the model showed that progressive increases in smooth muscle tone, up to a critical level, resulted in progressive but uniform airway constriction. However, smooth muscle activation beyond that critical level led to a very heterogeneous pattern of ventilation and airway obstruction with large clusters of constricted terminal bronchioles 4 equivalent to large VDefs observed with PET 15. Further activation of smooth muscle in the model expanded the extent of the VDefs as well as the number of severely constricted terminal bronchioles within them, while the rest of the lung remained uniformly ventilated. The findings of the present study are consistent with the following predictions of that theoretical model 4. 1) Exclusive constriction of airways leading to ventilation defects, whether peripheral or more central was not enough to explain the obstructive changes in lung mechanics measured in the asthmatic subjects. To explain those changes, a certain degree of heterogeneous airway constriction to the rest of the airway tree was also required. 2) Although the differences in mechanical response to bronchoconstriction amongst different subjects were associated with differences in the extent of their VDefs, the constriction of airways leading to the rest of the lung was not associated with the inter-subject differences in the extent of the VDefs. In other words, it appears that once VDefs start to form, further smooth muscle constriction increases the relative extent of the VDefs but not the degree of constriction of the airways leading to the rest of the lung. 3) A non-unity y-axis intercept between the baseline-normalised conductance and the relative extent of VDefs is consistent with the prediction that a critical level of airway narrowing needs to be reached before sizable ventilation defects are developed. In contrast to the highly uniform distribution of ventilation predicted by the theoretically symmetrical network model for regions outside the VDefs 4, the mechanical model in the present study showed that to explain the measured changes in oscillatory mechanics heterogeneity had to be present in airways that did not lead to VDefs. This is also consistent with results of regional analysis of the ventilation PET images 15. This result is to be expected because the asymmetry in structural and functional properties intrinsic to the asthmatic lung. The present study demonstrates that in addition to the large-scale ventilation heterogeneity caused by clustering of severely constricted distal airways in VDefs 4, an average degree and heterogeneity of constriction in the rest of the airway tree was required to explain the measured changes in mechanical obstruction. Such heterogeneity was similar among the subjects studied and was not related to the extent of the VDefs. However, differences in mechanical response amongst patients could be attributed, in part, to the extent of the lung involved in ventilatory defects for each patient. These findings are consistent with predictions of a recent integrative computational model of a constricting lung 4 in which smooth muscle activation above a critical level resulted in clusters of highly constricted distal airways and a certain level of constriction of the rest of the tree, albeit to a much lesser extent than that of airways leading to VDefs. Both experiments and modelling suggest that for smooth muscle activation exceeding that critical level, VDefs increase in their extent while the rest of the airway tree does not constrict any further.
Conclusion
The present study was supported by a National Institutes of Health Grant HL-68011 and the American Association of University Women Selected Professions Fellowship.
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||