Fractal analysis reveals functional unit of ventilation in the lung.

Ventilation is inhomogeneous in the lungs across species. It has been hypothesized that ventilation inhomogeneity is largely determined by the design of the airway branching network. Because exchange of gases at the alveolar barrier is more efficient when gas concentrations are evenly distributed at subacinar length scales, it is assumed that a 'functional unit' of ventilation exists within the lung periphery, where gas concentration becomes uniform. On the other hand, because the morphology of pulmonary airways and alveoli, and the distribution of inhaled fluorescent particles show self-similar fractal properties over a wide range of length scales, it has been predicted that fractal dimension of ventilation approaches unity within an internally homogeneous functional unit of ventilation. However, the existence of such a functional unit has never been demonstrated experimentally due to lack of in situ gas concentration measurements of sufficient spatial resolution in the periphery of a complex bifurcating network. Here, using energy-subtractive synchrotron radiation tomography, we measured the distribution of an inert gas (Xe) in the in vivo rabbit lung during Xe wash-in breathing manoeuvres. The effects of convective flow rate, diffusion and cardiac motion were also assessed. Fractal analysis of resulting gas concentration and tissue density maps revealed that fractal dimension was always smaller for Xe than for tissue density, and that only for the gas, a length scale existed where fractal dimension approached unity. The length scale where this occurred was seen to correspond to that of a rabbit acinus, the terminal structure comprising only alveolated airways. KEY POINTS: Gas ventilation is inhomogeneous in the lung of many species. However, it is not known down to what length scales this inhomogeneity persists. It is generally assumed that ventilation becomes homogeneous at subacinar length scales, beyond the spatial resolution of commonly available imaging techniques, hence this has not been demonstrated experimentally. Here we measured the distribution of inhaled Xe gas in the rabbit lung using synchrotron radiation energy-subtractive imaging and used fractal analysis to show that ventilation becomes internally uniform within regions about the size of rabbit lung acini.

STRUCTURE-FUNCTION IN SMOKERS : WHEN A SMALL AIRWAYS TEST REALLY REFLECTS THE SMALL AIRWAYS.

INTRODUCTION: If Multiple Breath Washout (MBW) derived acinar ventilation heterogeneity (Sacin) really represents peripheral units, the N2 phase-III of the first MBW exhalation should be curvilinear. This is due to the superposed effect of gas diffusion and convection resulting in an equilibration of N2 concentrations between neighbouring lung units throughout exhalation. We investigated this in smokers with CT-proven functional small airway disease. METHODS: Instantaneous N2slopes were computed over 40ms intervals throughout phase-III and normalized by mean phase-III N2 concentration. N2phase-III (concave) curvilinearity was quantified as the rate at which the instantaneous N2 slope decreases past the phase-II peak over a 1s interval; for a linear N2phase-III unaffected by diffusion, this rate would amount to 0L-1/s. N2phase-III curvilinearity was obtained on the experimental curves and on existing model simulations of N2curves from a normal peripheral lung model and one with missing terminal bronchioles (either 50% or 30% TB left). RESULTS: In forty-six smokers (66 (+8) years; 49 (+26) packyears) with CT-evidence of peripheral lung destruction, instantaneous N2 slope decrease was compared between those with (fSAD+fEmphys)>20% (-0.26+0.14(SD) L-1/s;n=24) and those with (fSAD+fEmphys)<20% (-0.16+0.12(SD)L-1/s;n=22) (P=0.014). Experimental values fell in the range predicted by a realistic peripheral lung model with progressive reduction of terminal bronchioles: values of instantaneous N2-slope decrease obtained from model simulations were -0.09L-1/s (normal lung;100%TB left), -0.17L-1/s (normal lung 50%TB left) and -0.29L-1/s (30%TB left). DISCUSSION : In smokers with CT-based evidence of functional small airways alterations, it is possible to demonstrate that Sacin really does represent the most peripheral airspaces.

Lung clearance index in adult cystic fibrosis patients: the role of convection-dependent lung units.

We investigated the extent of convective ventilation heterogeneity contributing to the observed lung clearance index values in adult cystic fibrosis patients, as the result of two simulated scenarios, either 1) a fixed part of the lungs becoming increasingly less compliant, or 2) an increasingly greater part of the lung being less compliant. In 25 cystic fibrosis patients and 25 matched controls, we computed the lung clearance index and also quantified curvilinearity of the washout concentration curve, where curvilinearity is equal to 0 (linear in semilog plot) reflects homogeneous ventilation and curvilinearity equal to 1 corresponds to the presence of an infinitesimally slowly emptying lung unit. In the cystic fibrosis group (forced expiratory volume in 1 s 27-100% predicted), lung clearance index and curvilinearity average±sd values were 10.3±2.3 and 0.57±0.13, respectively, and were significantly different from control values (6.2±0.4 and 0.18±0.07; both p<0.001); lung clearance index and curvilinearity were also correlated (R = 0.67; p<0.001). The average curvilinearity value in the cystic fibrosis group was found to be compatible with a cumulative volume of underventilated lung of 40-50%, depending on the simulation scenario considered. The degree of washout curvilinearity observed here indicates that a major determinant of the abnormal lung clearance index values observed in adult cystic fibrosis patients is ventilation heterogeneity generated between convection-dependent lung units, while the remainder of lung clearance index abnormality with respect to normal controls potentially represents the small airways within these lung zones.

Consensus statement for inert gas washout measurement using multiple- and single- breath tests.

Inert gas washout tests, performed using the single- or multiple-breath washout technique, were first described over 60 years ago. As measures of ventilation distribution inhomogeneity, they offer complementary information to standard lung function tests, such as spirometry, as well as improved feasibility across wider age ranges and improved sensitivity in the detection of early lung damage. These benefits have led to a resurgence of interest in these techniques from manufacturers, clinicians and researchers, yet detailed guidelines for washout equipment specifications, test performance and analysis are lacking. This manuscript provides recommendations about these aspects, applicable to both the paediatric and adult testing environment, whilst outlining the important principles that are essential for the reader to understand. These recommendations are evidence based, where possible, but in many places represent expert opinion from a working group with a large collective experience in the techniques discussed. Finally, the important issues that remain unanswered are highlighted. By addressing these important issues and directing future research, the hope is to facilitate the incorporation of these promising tests into routine clinical practice.

Correction: Vector mapping and bloodmeal metabarcoding demonstrate risk of urban Chagas disease transmission in Caracas, Venezuela.

[This corrects the article DOI: 10.1371/journal.pntd.0010613.].

Vector mapping and bloodmeal metabarcoding demonstrate risk of urban Chagas disease transmission in Caracas, Venezuela.

Chagas disease is a significant public health risk in rural and semi-rural areas of Venezuela. Triatomine infection by the aetiological agent Trypanosoma cruzi is also observed in the Metropolitan District of Caracas (MDC), where foodborne T. cruzi outbreaks occasionally occur but active vector-to-human transmission (infection during triatomine bloodmeal) is considered absent. Citizen science-based domiciliary triatomine collection carried out between 2007 and 2013 in the MDC has advanced understanding of urban T. cruzi prevalence patterns and represents an important public awareness-building tool. The present study reports on the extension of this triatomine collection program from 2014 to 2019 and uses mitochondrial metabarcoding to assess feeding behavior in a subset of specimens. The combined, thirteen-year dataset (n = 4872) shows a high rate of T. cruzi infection (75.2%) and a predominance of Panstrongylus geniculatus (99.01%) among triatomines collected in domiciliary areas by MDC inhabitants. Collection also involved nymphal stages of P. geniculatus in 18 of 32 MDC parishes. Other collected species included Triatoma nigromaculata, Triatoma maculata, Rhodnius prolixus, and Panstrongylus rufotuberculatus. Liquid intestinal content indicative of bloodmeal was observed in 53.4% of analyzed specimens. Dissection pools representing 108 such visually blooded P. geniculatus specimens predominantly tested positive for human cytochrome b DNA (22 of 24 pools). Additional bloodmeal sources detected via metabarcoding analysis included key sylvatic T. cruzi reservoirs (opossum and armadillo), rodents, and various other synanthropic and domesticated animals. Results suggest a porous sylvatic-domiciliary transmission interface and ongoing adaptation of P. geniculatus to the urban ecotope. Although P. geniculatus defecation traits greatly limit the possibility of active T. cruzi transmission for any individual biting event, the cumulation of this low risk across a vast metropolitan population warrants further investigation. Efforts to prevent triatomine contact with human food sources also clearly require greater attention to protect Venezuela's capital from Chagas disease.

Ventilation heterogeneity in the acinar and conductive zones of the normal ageing lung.

RATIONALE: Small airways function studies in lung disease have used three promising multiple breath washout (MBW) derived indices: indices of ventilation heterogeneity in the acinar (S(acin)) and conductive (S(cond)) lung zones, and the lung clearance index (LCI). Since peripheral lung structure is known to change with age, ventilation heterogeneity is expected to be affected too. However, the age dependence of the MBW indices of ventilation heterogeneity in the normal lung is unknown. OBJECTIVES: The authors systematically investigated S(acin), S(cond) or LCI as a function of age, testing also the robustness of these relationships across two laboratories. METHODS: MBW tests were performed by never-smokers (50% men) in the age range 25-65 years, with data gathered across two laboratories (n=120 and n=60). For comparison with the literature, the phase III slopes from classical single breath washout tests were also acquired in one group (n=120). MEASUREMENTS AND MAIN RESULTS: All three MBW indices consistently increased with age, representing a steady worsening of ventilation heterogeneity in the age range 25-65. Age explained 7-16% of the variability in S(acin) and S(cond) and 36% of the variability in LCI. There was a small but significant gender difference only for S(acin). Classical single breath washout phase III slopes also showed age dependencies, with gender effects depending on the normalisation method used. CONCLUSIONS: With respect to the clinical response, age is a small but consistent effect that needs to be factored in when using the MBW indices for the detection of small airways abnormality in disease.

Large urban outbreak of orally acquired acute Chagas disease at a school in Caracas, Venezuela.

BACKGROUND: Trypanosoma cruzi oral transmission is possible through food contamination by vector's feces. Little is known about the epidemiology and clinical features of microepidemics of orally acquired acute Chagas disease (CD). METHODS: A case-control, cohort-nested, epidemiological study was conducted during an outbreak of acute CD that affected a school community. Structured interviews were designed to identify symptoms and sources of infection. Electrocardiograms were obtained for all patients. Specific serum antibodies were assessed by immunoenzimatic and indirect hemagglutination tests. In some cases, parasitemia was tested directly or by culture, animal inoculation, and/or a polymerase chain reaction technique. RESULTS: Infection was confirmed in 103 of 1000 exposed individuals. Of those infected, 75% were symptomatic, 20.3% required hospitalization, 59% showed ECG abnormalities, parasitemia was documented in 44, and 1 child died. Clinical features differed from those seen in vectorial transmission. The infection rate was significantly higher among younger children. An epidemiological investigation incriminated contaminated fresh guava juice as the sole source of infection. CONCLUSIONS: This outbreak was unique, because it affected a large, urban, predominantly young, middle-class, otherwise healthy population and resulted in an unprecedented public health emergency. Rapid diagnosis and treatment avoided higher lethality. Food-borne transmission of T. cruzi may occur more often than is currently recognized.

Effect of airways constriction on exhaled nitric oxide.

While airway constriction has been shown to affect exhaled nitric oxide (NO), the mechanisms and location of constricted airways most likely to affect exhaled NO remain obscure. We studied the effects of histamine-induced airway constriction and ventilation heterogeneity on exhaled NO at 50 ml/s (Fe(NO,50)) and combined this with model simulations of Fe(NO,50) changes due to constriction of airways at various depths of the lung model. In 20 normal subjects, histamine induced a 26 +/- 15(SD)% Fe(NO,50) decrease, a 9 +/- 6% forced expiratory volume in 1 s (FEV(1)) decrease, a 19 +/- 9% mean forced midexpiratory flow between 25% and 75% forced vital capacity (FEF(25-75)) decrease, and a 94 +/- 119% increase in conductive ventilation heterogeneity. There was a significant correlation of Fe(NO,50) decrease with FEF(25-75) decrease (P = 0.006) but not with FEV(1) decrease or with increased ventilation heterogeneity. Simulations confirmed the negligible effect of ventilation heterogeneity on Fe(NO,50) and showed that the histamine-induced Fe(NO,50) decrease was due to constriction, with associated reduction in NO flux, of airways located proximal to generation 15. The model also indicated that the most marked effect of airways constriction on Fe(NO,50) is situated in generations 10-15 and that airway constriction beyond generation 15 markedly increases Fe(NO,50) due to interference with the NO backdiffusion effect. These mechanical factors should be considered when interpreting exhaled NO in lung disease.

Contribution of upper airway geometry to convective mixing.

We investigated the axial dispersive effect of the upper airway structure (comprising mouth cavity, oropharynx, and trachea) on a traversing aerosol bolus. This was done by means of aerosol bolus experiments on a hollow cast of a realistic upper airway model (UAM) and three-dimensional computational fluid dynamics (CFD) simulations in the same UAM geometry. The experiments showed that 50-ml boluses injected into the UAM dispersed to boluses with a half-width ranging from 80 to 90 ml at the UAM exit, across both flow rates (250, 500 ml/s) and both flow directions (inspiration, expiration). These experimental results imply that the net half-width induced by the UAM typically was 69 ml. Comparison of experimental bolus traces with a one-dimensional Gaussian-derived analytical solution resulted in an axial dispersion coefficient of 200-250 cm(2)/s, depending on whether the bolus peak and its half-width or the bolus tail needed to be fully accounted for. CFD simulations agreed well with experimental results for inspiratory boluses and were compatible with an axial dispersion of 200 cm(2)/s. However, for expiratory boluses the CFD simulations showed a very tight bolus peak followed by an elongated tail, in sharp contrast to the expiratory bolus experiments. This indicates that CFD methods that are widely used to predict the fate of aerosols in the human upper airway, where flow is transitional, need to be critically assessed, possibly via aerosol bolus simulations. We conclude that, with all its geometric complexity, the upper airway introduces a relatively mild dispersion on a traversing aerosol bolus for normal breathing flow rates in inspiratory and expiratory flow directions.

A simulation study of diffusion-convection interaction and its effect on multiple breath washout indices.

Experimental studies of acinar ventilation heterogeneity (Sacin) derived from the multiple breath washout have shown the potential of Sacin to pick up structural change in lung disease. Recent Sacin data suggest that even when intra-acinar structure is unaltered, the combination of convection, diffusion and number of acini fed by patent terminal bronchioles can modify Sacin. We show here how Sacin is affected by structural features such as the secondary alveolar septa, intra-acinar ramification and number of ventilated acini. The simulations also predict relationships between respective alterations in Sacin and washout indices such as lung clearance index (LCI) and alveolar mixing efficiency (AME). This was verified experimentally, with highly significant correlations between Sacin and LCI (r = +0.85;p < 0.001) and between Sacin and AME (r = - 0.92; p < 0.001). We have shown how acinar ventilation heterogeneity can be affected by a reduction of number of ventilated acini, a change in overall alveolization or in intra-acinar alveolization pattern, via their impact on the balance between convection and diffusion at acinar level.

Simulation of the apparent diffusion of helium-3 in the human acinus.

Functional MRI of the lungs with hyperpolarized helium provides an index of apparent diffusion measured over several seconds (ADCsec) that is only 2% of its free diffusion in air (0.88 cm2/s). The potential of ADCsec to noninvasively assess in vivo lung structure of diseased lungs at the length scales corresponding to several seconds is critically dependent on the exact link between ADCsec and lung peripheral structure. To understand the intruigingly small ADCsec, numerical simulations of gas transport were performed in 1) a trumpet model, 2) a symmetrical, and 3) an asymmetrical multiple-branch-point model of the human acinus. For initial gas boluses in different locations of the acinar models, ADCsec was quantified as follows. At different time intervals, we computed a coefficient of variation (CoV) of the concentration distributions within each acinar model. The slope in the semilog plot of log(CoV) vs. time was proportional to the ADCsec generated by the internal model structure, provided that the outer model boundaries were similar across all models (i.e., similar cumulative cross section vs. average path length). The simulations revealed an ADCsec that amounted to approximately 1% of free diffusion in the trumpet model of the acinus, i.e., corresponding to free diffusion within the acinar geometric boundaries. Our simulations show that for initial conditions corresponding to those used in MRI experiments, intra-acinar branching introduces a dramatic diffusion delay, comparable to what is observed experimentally.

Single-breath test in lateral decubitus reflects function of single lungs grafted for emphysema.

The slope of alveolar plateau for nitrogen derived from the single-breath test is useful to assess the function of bilateral lung grafts, but this technique is not applicable to patients with single-lung grafts due to the confounding influence of the native lung. We tested the hypothesis that the nitrogen slope measured in lateral decubitus with the graft in nondependent position may primarily reflect the distribution of ventilation in this lung. Fifteen patients with single-lung transplantation for emphysema, 10 healthy controls, and 7 patients with advanced emphysema performed single-breath washouts in right and left lateral decubitus; nitrogen slope was measured between 75 and 100% of expired volume. In 10 transplant recipients, the volume of each lung was measured in the two postures by computerized tomography. Nitrogen slope was unaffected by posture in normal controls and emphysema patients. On the other hand, nitrogen slope in transplant recipients was invariably smaller, with the graft in nondependent vs. in dependent position. Values of nitrogen slope with the graft in nondependent position were similar to those obtained in normal controls but significantly smaller than those obtained in emphysema patients. Computerized tomography studies in this position indicated that the volume expired below functional residual capacity was exclusively contributed by the graft. We conclude that, in patients with single-lung transplantation for emphysema, 1) measuring nitrogen slope in lateral decubitus allows to distinguish between the graft and the native lung, and 2) nitrogen slope obtained with the graft in nondependent position reflects ventilation distribution in this lung.

Anatomically based three-dimensional model of airways to simulate flow and particle transport using computational fluid dynamics.

We have studied gas flow and particle deposition in a realistic three-dimensional (3D) model of the bronchial tree, extending from the trachea to the segmental bronchi (7th airway generation for the most distal ones) using computational fluid dynamics. The model is based on the morphometrical data of Horsfield et al. (Horsfield K, Dart G, Olson DE, Filley GF, and Cumming G. J Appl Physiol 31: 207-217, 1971) and on bronchoscopic and computerized tomography images, which give the spatial 3D orientation of the curved ducts. It incorporates realistic angles of successive branching planes. Steady inspiratory flow varying between 50 and 500 cm(3)/s was simulated, as well as deposition of spherical aerosol particles (1-7 microm diameter, 1 g/cm(3) density). Flow simulations indicated nonfully developed flows in the branches due to their relative short lengths. Velocity flow profiles in the segmental bronchi, taken one diameter downstream of the bifurcation, were distorted compared with the flow in a simple curved tube, and wide patterns of secondary flow fields were observed. Both were due to the asymmetrical 3D configuration of the bifurcating network. Viscous pressure drop in the model was compared with results obtained by Pedley et al. (Pedley TJ, Schroter RC, and Sudlow MF. Respir Physiol 9: 387-405, 1970), which are shown to be a good first approximation. Particle deposition increased with particle size and was minimal for approximately 200 cm(3)/s inspiratory flow, but it was highly heterogeneous for branches of the same generation.