Quantitative assessment of lung microstructure in healthy mice using an MR-based 3He lung morphometry technique

Acinar micromechanics in health and lung injury: what we have learned from quantitative morphology

Within the pulmonary acini ventilation and blood perfusion are brought together on a huge surface area separated by a very thin blood-gas barrier of tissue components to allow efficient gas exchange. During ventilation pulmonary acini are cyclically subjected to deformations which become manifest in changes of the dimensions of both alveolar and ductal airspaces as well as the interalveolar septa, composed of a dense capillary network and the delicate tissue layer forming the blood-gas barrier. These ventilation-related changes are referred to as micromechanics. In lung diseases, abnormalities in acinar micromechanics can be linked with injurious stresses and strains acting on the blood-gas barrier. The mechanisms by which interalveolar septa and the blood-gas barrier adapt to an increase in alveolar volume have been suggested to include unfolding, stretching, or changes in shape other than stretching and unfolding. Folding results in the formation of pleats in which alveolar epithelium is not exposed to air and parts of the blood-gas barrier are folded on each other. The opening of a collapsed alveolus (recruitment) can be considered as an extreme variant of septal wall unfolding. Alveolar recruitment can be detected with imaging techniques which achieve light microscopic resolution. Unfolding of pleats and stretching of the blood-gas barrier, however, require electron microscopic resolution to identify the basement membrane. While stretching results in an increase of the area of the basement membrane, unfolding of pleats and shape changes do not. Real time visualization of these processes, however, is currently not possible. In this review we provide an overview of septal wall micromechanics with focus on unfolding/folding as well as stretching. At the same time we provide a state-of-the-art design-based stereology methodology to quantify microarchitecture of alveoli and interalveolar septa based on different imaging techniques and design-based stereology.

Preclinical MRI Using Hyperpolarized 129Xe

Although critical for development of novel therapies, understanding altered lung function in disease models is challenging because the transport and diffusion of gases over short distances, on which proper function relies, is not readily visualized. In this review we summarize progress introducing hyperpolarized ¹²⁹Xe imaging as a method to follow these processes in vivo. The work is organized in sections highlighting methods to observe the gas replacement effects of breathing (Gas Dynamics during the Breathing Cycle) and gas diffusion throughout the parenchymal airspaces (3). We then describe the spectral signatures indicative of gas dissolution and uptake (4), and how these features can be used to follow the gas as it enters the tissue and capillary bed, is taken up by hemoglobin in the red blood cells (5), re-enters the gas phase prior to exhalation (6), or is carried via the vasculature to other organs and body structures (7). We conclude with a discussion of practical imaging and spectroscopy techniques that deliver quantifiable metrics despite the small size, rapid motion and decay of signal and coherence characteristic of the magnetically inhomogeneous lung in preclinical models (8).

Pulmonary acini exhibit complex changes during postnatal rat lung development

Pulmonary acini represent the functional gas-exchanging units of the lung. Due to technical limitations, individual acini cannot be identified on microscopic lung sections. To overcome these limitations, we imaged the right lower lobes of instillation-fixed rat lungs from postnatal days P4, P10, P21, and P60 at the TOMCAT beamline of the Swiss Light Source synchrotron facility at a voxel size of 1.48 μm. Individual acini were segmented from the three-dimensional data by closing the airways at the transition from conducting to gas exchanging airways. For a subset of acini (N = 268), we followed the acinar development by stereologically assessing their volume and their number of alveoli. We found that the mean volume of the acini increases 23 times during the observed time-frame. The coefficients of variation dropped from 1.26 to 0.49 and the difference between the mean volumes of the fraction of the 20% smallest to the 20% largest acini decreased from a factor of 27.26 (day 4) to a factor of 4.07 (day 60), i.e. shows a smaller dispersion at later time points. The acinar volumes show a large variation early in lung development and homogenize during maturation of the lung by reducing their size distribution by a factor of 7 until adulthood. The homogenization of the acinar sizes hints at an optimization of the gas-exchange region in the lungs of adult animals and that acini of different size are not evenly distributed in the lungs. This likely leads to more homogeneous ventilation at later stages in lung development.

Validating in vivo hyperpolarized 129 Xe diffusion MRI and diffusion morphometry in the mouse lung

PURPOSE

Diffusion and lung morphometry imaging using hyperpolarized gases are promising tools to quantify pulmonary microstructure noninvasively in humans and in animal models. These techniques assume the motion encoded is exclusively diffusive gas displacement, but the impact of cardiac motion on measurements has never been explored. Furthermore, although diffusion morphometry has been validated against histology in humans and mice using 3 He, it has never been validated in mice for 129 Xe. Here, we examine the effect of cardiac motion on diffusion imaging and validate 129 Xe diffusion morphometry in mice.

THEORY AND METHODS

Mice were imaged using gradient-echo-based diffusion imaging, and apparent diffusion-coefficient (ADC) maps were generated with and without cardiac gating. Diffusion-weighted images were fit to a previously developed theoretical model using Bayesian probability theory, producing morphometric parameters that were compared with conventional histology.

RESULTS

Cardiac gating had no significant impact on ADC measurements (dual-gating: ADC = 0.020 cm2 /s, single-gating: ADC = 0.020 cm2 /s; P = .38). Diffusion-morphometry-generated maps of ADC (mean, 0.0165 ± 0.0001 cm2 /s) and acinar dimensions (alveolar sleeve depth [h] = 44 µm, acinar duct radii [R] = 99 µm, mean linear intercept [Lm ] = 74 µm) that agreed well with conventional histology (h = 45 µm, R = 108 µm, Lm = 63 µm).

CONCLUSION

Cardiac motion has negligible impact on 129 Xe ADC measurements in mice, arguing its impact will be similarly minimal in humans, where relative cardiac motion is reduced. Hyperpolarized 129 Xe diffusion morphometry accurately and noninvasively maps the dimensions of lung microstructure, suggesting it can quantify the pulmonary microstructure in mouse models of lung disease.

Preclinical hyperpolarized 129 Xe MRI: ventilation and T2 * mapping in mouse lungs at 7 T using multi-echo flyback UTE

Fast apparent transverse relaxation (short T₂ *) is a common obstacle when attempting to perform quantitative ¹ H MRI of the lungs. While T₂ * times are longer for pulmonary hyperpolarized (HP) gas functional imaging (in particular for gaseous ¹²⁹ Xe), T₂ * can still lead to quantitative inaccuracies for sequences requiring longer echo times (such as diffusion weighted images) or longer readout duration (such as spiral sequences). This is especially true in preclinical studies, where high magnetic fields lead to shorter relaxation times than are typically seen in human studies. However, the T₂ * of HP ¹²⁹ Xe in the most common animal model of human disease (mice) has not been reported. Herein, we present a multi-echo radial flyback imaging sequence and use it to measure HP ¹²⁹ Xe T₂ * at 7 T under a variety of respiratory conditions. This sequence mitigates the impact of T₁ relaxation outside the animal by using multiple gradient-refocused echoes to acquire images at a number of effective echo times for each RF excitation. After validating the sequence using a phantom containing water doped with superparamagnetic iron oxide nanoparticles, we measured the ¹²⁹ Xe T₂ * in vivo for 10 healthy C57Bl/6 J mice and found T₂ * ~ 5 ms in the lung airspaces. Interestingly, T₂ * was relatively constant over all experimental conditions, and varied significantly with sex, but not age, mass, or the O₂ content of the inhaled gas mixture. These results are discussed in the context of T₂ * relaxation within porous media.

Alveolar dynamics during mechanical ventilation in the healthy and injured lung

Mechanical ventilation is a life-saving therapy in patients with acute respiratory distress syndrome (ARDS). However, mechanical ventilation itself causes severe co-morbidities in that it can trigger ventilator-associated lung injury (VALI) in humans or ventilator-induced lung injury (VILI) in experimental animal models. Therefore, optimization of ventilation strategies is paramount for the effective therapy of critical care patients. A major problem in the stratification of critical care patients for personalized ventilation settings, but even more so for our overall understanding of VILI, lies in our limited insight into the effects of mechanical ventilation at the actual site of injury, i.e., the alveolar unit. Unfortunately, global lung mechanics provide for a poor surrogate of alveolar dynamics and methods for the in-depth analysis of alveolar dynamics on the level of individual alveoli are sparse and afflicted by important limitations. With alveolar dynamics in the intact lung remaining largely a "black box," our insight into the mechanisms of VALI and VILI and the effectiveness of optimized ventilation strategies is confined to indirect parameters and endpoints of lung injury and mortality.In the present review, we discuss emerging concepts of alveolar dynamics including alveolar expansion/contraction, stability/instability, and opening/collapse. Many of these concepts remain still controversial, in part due to limitations of the different methodologies applied. We therefore preface our review with an overview of existing technologies and approaches for the analysis of alveolar dynamics, highlighting their individual strengths and limitations which may provide for a better appreciation of the sometimes diverging findings and interpretations. Joint efforts combining key technologies in identical models to overcome the limitations inherent to individual methodologies are needed not only to provide conclusive insights into lung physiology and alveolar dynamics, but ultimately to guide critical care patient therapy.

The role of three-dimensionality and alveolar pressure in the distribution and amplification of alveolar stresses

Alveolar stresses are fundamental to enable the respiration process in mammalians and have recently gained increasing attention due to their mechanobiological role in the pathogenesis and development of respiratory diseases. Despite the fundamental physiological role of stresses in the alveolar wall, the determination of alveolar stresses remains challenging, and our current knowledge is largely drawn from 2D studies that idealize the alveolar septal wall as a spring or a planar continuum. Here we study the 3D stress distribution in alveolar walls of normal lungs by combining ex-vivo micro-computed tomography and 3D finite-element analysis. Our results show that alveolar walls are subject to a fully 3D state of stresses rather than to a pure axial stress state. To understand the contributions of the different components and deformation modes, we decompose the stress tensor field into hydrostatic and deviatoric components, which are associated with isotropic and distortional stresses, respectively. Stress concentrations arise in localized regions of the alveolar microstructure, with magnitudes that can be up to 27 times the applied alveolar pressure. Interestingly, we show that the stress amplification factor strongly depends on the level of alveolar pressure, i.e, stresses do not scale proportional to the applied alveolar pressure. In addition, we show that 2D techniques to assess alveolar stresses consistently overestimate the stress magnitude in alveolar walls, particularly for lungs under high transpulmonary pressure. These findings take particular relevance in the study of stress-induced remodeling of the emphysematous lung and in ventilator-induced lung injury, where the relation between transpulmonary pressure and alveolar wall stress is key to understand mechanotransduction processes in pneumocytes.

Comparative Assessment of Aspergillosis by Virtual Infection Modeling in Murine and Human Lung

Aspergillus fumigatus is a ubiquitous opportunistic fungal pathogen that can cause severe infections in immunocompromised patients. Conidia that reach the lower respiratory tract are confronted with alveolar macrophages, which are the resident phagocytic cells, constituting the first line of defense. If not efficiently removed in time, A. fumigatus conidia can germinate causing severe infections associated with high mortality rates. Mice are the most extensively used model organism in research on A. fumigatus infections. However, in addition to structural differences in the lung physiology of mice and the human host, applied infection doses in animal experiments are typically orders of magnitude larger compared to the daily inhalation doses of humans. The influence of these factors, which must be taken into account in a quantitative comparison and knowledge transfer from mice to humans, is difficult to measure since in vivo live cell imaging of the infection dynamics under physiological conditions is currently not possible. In the present study, we compare A. fumigatus infection in mice and humans by virtual infection modeling using a hybrid agent-based model that accounts for the respective lung physiology and the impact of a wide range of infection doses on the spatial infection dynamics. Our computer simulations enable comparative quantification of A. fumigatus infection clearance in the two hosts to elucidate (i) the complex interplay between alveolar morphometry and the fungal burden and (ii) the dynamics of infection clearance, which for realistic fungal burdens is found to be more efficiently realized in mice compared to humans.

Advanced pulmonary MRI to quantify alveolar and acinar duct abnormalities: Current status and future clinical applications

There are serious clinical gaps in our understanding of chronic lung disease that require novel, sensitive, and noninvasive in vivo measurements of the lung parenchyma to measure disease pathogenesis and progressive changes over time as well as response to treatment. Until recently, our knowledge and appreciation of the tissue changes that accompany lung disease has depended on ex vivo biopsy and concomitant histological and stereological measurements. These measurements have revealed the underlying pathologies that drive lung disease and have provided important observations about airway occlusion, obliteration of the terminal bronchioles and airspace enlargement, or fibrosis and their roles in disease initiation and progression. ex vivo tissue stereology and histology are the established gold standards and, more recently, micro-computed tomography (CT) measurements of ex vivo tissue samples has also been employed to reveal new mechanistic findings about the progression of obstructive lung disease in patients. While these approaches have provided important understandings using ex vivo analysis of excised samples, recently developed hyperpolarized noble gas MRI methods provide an opportunity to noninvasively measure acinar duct and terminal airway dimensions and geometry in vivo, and, without radiation burden. Therefore, in this review we summarize emerging pulmonary MRI morphometry methods that provide noninvasive in vivo measurements of the lung in patients with bronchopulmonary dysplasia and chronic obstructive pulmonary disease, among others. We discuss new findings, future research directions, as well as clinical opportunities to address current gaps in patient care and for testing of new therapies. Level of Evidence: 5 Technical Efficacy: Stage 5 J. Magn. Reson. Imaging 2019;50:28-40.

Design of biomimetic substrates for long-term maintenance of alveolar epithelial cells

There is a need to establish in vitro lung alveolar epithelial culture models to better understand the fundamental biological mechanisms that drive lung diseases. While primary alveolar epithelial cells (AEC) are a useful option to study mature lung biology, they have limited utility in vitro. Cells that survive demonstrate limited proliferative capacity and loss of phenotype over the first 3-5 days in traditional culture conditions. To address this limitation, we generated a novel physiologically relevant cell culture system for enhanced viability and maintenance of phenotype. Here we describe a method utilizing e-beam lithography, reactive ion etching, and replica molding to generate poly-dimethylsiloxane (PDMS) substrates containing hemispherical cavities that mimic the architecture and size of mouse and human alveoli. Primary AECs grown on these cavity-containing substrates form a monolayer that conforms to the substrate enabling precise control over cell sheet architecture. AECs grown in cavity culture conditions remain viable and maintain their phenotype over one week. Specifically, cells grown on substrates consisting of 50 μm diameter cavities remained 96 ± 4% viable and maintained expression of surfactant protein C (SPC), a marker of type 2 AEC over 7 days. While this report focuses on primary lung alveolar epithelial cells, our culture platform is potentially relevant and useful for growing primary cells from other tissues with similar cavity-like architecture and could be further adapted to other biomimetic shapes or contours.

Optical model of the murine lung to optimize pulmonary illumination

We describe a Monte Carlo model of the mouse torso to optimize illumination of the mouse lung for fluorescence detection of low levels of pulmonary pathogens, specifically Mycobacterium tuberculosis. After validation of the simulation with an internally illuminated optical phantom, the entire mouse torso was simulated to compare external and internal illumination techniques. Measured optical properties of deflated mouse lungs were scaled to mimic the diffusive properties of inflated lungs in vivo. Using the full-torso model, a 2  ×   to 3  ×   improvement in average fluence rate in the lung was seen for dorsal compared with ventral positioning of the mouse with external illumination. The enhancement in average fluence rate in the lung using internal excitation was 40  ×   to 60  ×   over external illumination in the dorsal position. Parameters of the internal fiber optic source were manipulated in the model to guide optimization of the physical system and experimental protocol for internal illumination and whole-body detection of fluorescent mycobacteria in a mouse model of infection.

Hyperpolarized gas diffusion MRI of biphasic lung inflation in short- and long-term emphysema models

During lung inflation, airspace dimensions are affected nonlinearly by both alveolar expansion and recruitment, potentially confounding the identification of emphysematous lung by hyperpolarized helium-3 diffusion magnetic resonance imaging (HP MRI). This study aimed to characterize lung inflation over a broad range of inflation volume and pressure values in two different models of emphysema, as well as in normal lungs. Elastase-treated rats (n = 7) and healthy controls (n = 7) were imaged with HP MRI. Gradual inflation was achieved by incremental changes to both inflation volume and airway pressure. The apparent diffusion coefficient (ADC) was measured at each level of inflation and fitted to the corresponding airway pressures as the second-order response equation, with minimizing residue (χ² < 0.001). A biphasic ADC response was detected, with an initial ADC increase followed by a decrease at airway pressures >18 cmH₂O. Discrimination between treated and control rats was optimal when airway pressure was intermediate (between 10 and 11 cmH₂O). Similar findings were confirmed in mice following long-term exposure to cigarette smoke, where optimal discrimination between treated and healthy mice occurred at a similar airway pressure as in the rats. We subsequently explored the evolution of ADC measured at the intermediate inflation level in mice after prolonged smoke exposure and found a significant increase (P < 0.01) in ADC over time. Our results demonstrate that measuring ADC at intermediate inflation enhances the distinction between healthy and diseased lungs, thereby establishing a model that may improve the diagnostic accuracy of future HP gas diffusion studies.

Diffusion lung imaging with hyperpolarized gas MRI

Lung imaging using conventional 1 H MRI presents great challenges because of the low density of lung tissue, lung motion and very fast lung tissue transverse relaxation (typical T2 * is about 1-2 ms). MRI with hyperpolarized gases (3 He and 129 Xe) provides a valuable alternative because of the very strong signal originating from inhaled gas residing in the lung airspaces and relatively slow gas T2 * relaxation (typical T2 * is about 20-30 ms). However, in vivo human experiments should be performed very rapidly - usually during a single breath-hold. In this review, we describe the recent developments in diffusion lung MRI with hyperpolarized gases. We show that a combination of the results of modeling of gas diffusion in lung airspaces and diffusion measurements with variable diffusion-sensitizing gradients allows the extraction of quantitative information on the lung microstructure at the alveolar level. From an MRI scan of less than 15 s, this approach, called in vivo lung morphometry, allows the provision of quantitative values and spatial distributions of the same physiological parameters as measured by means of 'standard' invasive stereology (mean linear intercept, surface-to-volume ratio, density of alveoli, etc.). In addition, the approach makes it possible to evaluate some advanced Weibel parameters characterizing lung microstructure: average radii of alveolar sacs and ducts, as well as the depth of their alveolar sleeves. Such measurements, providing in vivo information on the integrity of pulmonary acinar airways and their changes in different diseases, are of great importance and interest to a broad range of physiologists and clinicians. We also discuss a new type of experiment based on the in vivo lung morphometry technique combined with quantitative computed tomography measurements, as well as with gradient echo MRI measurements of hyperpolarized gas transverse relaxation in the lung airspaces. Such experiments provide additional information on the blood vessel volume fraction, specific gas volume and length of the acinar airways, and allow the evaluation of lung parenchymal and non-parenchymal tissue. Copyright © 2015 John Wiley & Sons, Ltd.

Detection of the mild emphysema by quantification of lung respiratory airways with hyperpolarized xenon diffusion MRI

PURPOSE

To demonstrate the feasibility to quantify the lung respiratory airway in vivo with hyperpolarized xenon diffusion magnetic resonance imaging (MRI), which is able to detect mild emphysema in the rat model.

MATERIALS AND METHODS

The lung respiratory airways were quantified in vivo using hyperpolarized xenon diffusion MRI (7T) with eight b values (5, 10, 15, 20, 25, 30, 35, 40 s/cm² ) in five control rats and five mild emphysematous rats, which were induced by elastase. The morphological results from histology were acquired and used for comparison.

RESULTS

The parameters D_L (longitudinal diffusion coefficient), r (internal radius), h (alveolar sleeve depth), Lm (mean linear intercept), and S/V (surface area to lung volume ratio) derived from the hyperpolarized xenon diffusion MRI in the emphysematous group showed significant differences from those in the control group (P < 0.05). Additionally, these parameters correlated well with the Lm obtained by the traditional histological sections (Pearson's correlation coefficients >0.8).

CONCLUSION

The lung respiratory airways can be quantified by hyperpolarized xenon diffusion MRI, showing the potential for mild emphysema diagnosis. Also, the study suggested that the hyperpolarized xenon D_L is more sensitive than D_T (transverse diffusion coefficient) to detect mild emphysema.

LEVEL OF EVIDENCE

1 J. Magn. Reson. Imaging 2017;45:879-888.

Mechanical Forces Accelerate Collagen Digestion by Bacterial Collagenase in Lung Tissue Strips

Most tissues in the body are under mechanical tension, and while enzymes mediate many cellular and extracellular processes, the effects of mechanical forces on enzyme reactions in the native extracellular matrix (ECM) are not fully understood. We hypothesized that physiological levels of mechanical forces are capable of modifying the activity of collagenase, a key remodeling enzyme of the ECM. To test this, lung tissue Young's modulus and a nonlinearity index characterizing the shape of the stress-strain curve were measured in the presence of bacterial collagenase under static uniaxial strain of 0, 20, 40, and 80%, as well as during cyclic mechanical loading with strain amplitudes of ±10 or ±20% superimposed on 40% static strain, and frequencies of 0.1 or 1 Hz. Confocal and electron microscopy was used to determine and quantify changes in ECM structure. Generally, mechanical loading increased the effects of enzyme activity characterized by an irreversible decline in stiffness and tissue deterioration seen on both confocal and electron microscopic images. However, a static strain of 20% provided protection against digestion compared to both higher and lower strains. The decline in stiffness during digestion positively correlated with the increase in equivalent alveolar diameters and negatively correlated with the nonlinearity index. These results suggest that the decline in stiffness results from rupture of collagen followed by load transfer and subsequent rupture of alveolar walls. This study may provide new understanding of the role of collagen degradation in general tissue remodeling and disease progression.

Probing lung microstructure with hyperpolarized noble gas diffusion MRI: theoretical models and experimental results

The introduction of hyperpolarized gases ((3)He and (129)Xe) has opened the door to applications for which gaseous agents are uniquely suited-lung MRI. One of the pulmonary applications, diffusion MRI, relies on measuring Brownian motion of inhaled hyperpolarized gas atoms diffusing in lung airspaces. In this article we provide an overview of the theoretical ideas behind hyperpolarized gas diffusion MRI and the results obtained over the decade-long research. We describe a simple technique based on measuring gas apparent diffusion coefficient (ADC) and an advanced technique, in vivo lung morphometry, that quantifies lung microstructure both in terms of Weibel parameters (acinar airways radii and alveolar depth) and standard metrics (mean linear intercept, surface-to-volume ratio, and alveolar density) that are widely used by lung researchers but were previously available only from invasive lung biopsy. This technique has the ability to provide unique three-dimensional tomographic information on lung microstructure from a less than 15 s MRI scan with results that are in good agreement with direct histological measurements. These safe and sensitive diffusion measurements improve our understanding of lung structure and functioning in health and disease, providing a platform for monitoring the efficacy of therapeutic interventions in clinical trials.

Experimental evidence of age-related adaptive changes in human acinar airways

The progressive decline of lung function with aging is associated with changes in lung structure at all levels, from conducting airways to acinar airways (alveolar ducts and sacs). While information on conducting airways is becoming available from computed tomography, in vivo information on the acinar airways is not conventionally available, even though acini occupy 95% of lung volume and serve as major gas exchange units of the lung. The objectives of this study are to measure morphometric parameters of lung acinar airways in living adult humans over a broad range of ages by using an innovative MRI-based technique, in vivo lung morphometry with hyperpolarized (3)He gas, and to determine the influence of age-related differences in acinar airway morphometry on lung function. Pulmonary function tests and MRI with hyperpolarized (3)He gas were performed on 24 healthy nonsmokers aged 19-71 years. The most significant age-related difference across this population was a 27% loss of alveolar depth, h, leading to a 46% increased acinar airway lumen radius, hence, decreased resistance to acinar air transport. Importantly, the data show a negative correlation between h and the pulmonary function measures forced expiratory volume in 1 s and forced vital capacity. In vivo lung morphometry provides unique information on age-related changes in lung microstructure and their influence on lung function. We hypothesize that the observed reduction of alveolar depth in subjects with advanced aging represents a remodeling process that might be a compensatory mechanism, without which the pulmonary functional decline due to other biological factors with advancing age would be significantly larger.

Early stage radiation-induced lung injury detected using hyperpolarized (129) Xe Morphometry: Proof-of-concept demonstration in a rat model

PURPOSE

Radiation-induced lung injury (RILI) is still the major dose-limiting toxicity related to lung cancer radiation therapy, and it is difficult to predict and detect patients who are at early risk of severe pneumonitis and fibrosis. The goal of this proof-of-concept preclinical demonstration was to investigate the potential of hyperpolarized (129) Xe diffusion-weighted MRI to detect the lung morphological changes associated with early stage RILI.

METHODS

Hyperpolarized (129) Xe MRI was performed using eight different diffusion sensitizations (0.0-115 s/cm(2) ) in a small group of control rats (n = 4) and rats 2 wk after radiation exposure (n = 5). The diffusion-weighted images were used to obtain morphological estimates of the pulmonary parenchyma including external radius (R), internal radius (r), alveolar sleeve depth (h), and mean airspace chord length (Lm ). The histological mean linear intercept (MLI) were obtained for five control and five irradiated animals.

RESULTS

Mean R, r, and Lm were both significantly different (P < 0.02) in the irradiated rats (74 ± 17 µm, 43 ± 12 µm, and 54 ± 17 µm, respectively) compared with the control rats (100 ± 12 µm, 67 ± 10 µm, and 79 ± 12 µm, respectively). Changes in measured Lm values were consistent with changes in MLI values observed by histology.

CONCLUSIONS

Hyperpolarized (129) Xe MRI provides a way to detect and measure regional microanatomical changes in lung parenchyma in a preclinical model of RILI. Magn Reson Med 75:2421-2431, 2016. © 2015 Wiley Periodicals, Inc.

Interactions of Francisella tularensis with Alveolar Type II Epithelial Cells and the Murine Respiratory Epithelium

Francisella tularensis is classified as a Tier 1 select agent by the CDC due to its low infectious dose and the possibility that the organism can be used as a bioweapon. The low dose of infection suggests that Francisella is unusually efficient at evading host defenses. Although ~50 cfu are necessary to cause human respiratory infection, the early interactions of virulent Francisella with the lung environment are not well understood. To provide additional insights into these interactions during early Francisella infection of mice, we performed TEM analysis on mouse lungs infected with F. tularensis strains Schu S4, LVS and the O-antigen mutant Schu S4 waaY::TrgTn. For all three strains, the majority of the bacteria that we could detect were observed within alveolar type II epithelial cells at 16 hours post infection. Although there were no detectable differences in the amount of bacteria within an infected cell between the three strains, there was a significant increase in the amount of cellular debris observed in the air spaces of the lungs in the Schu S4 waaY::TrgTn mutant compared to either the Schu S4 or LVS strain. We also studied the interactions of Francisella strains with human AT-II cells in vitro by characterizing the ability of these three strains to invade and replicate within these cells. Gentamicin assay and confocal microscopy both confirmed that F. tularensis Schu S4 replicated robustly within these cells while F. tularensis LVS displayed significantly lower levels of growth over 24 hours, although the strain was able to enter these cells at about the same level as Schu S4 (1 organism per cell), as determined by confocal imaging. The Schu S4 waaY::TrgTn mutant that we have previously described as attenuated for growth in macrophages and mouse virulence displayed interesting properties as well. This mutant induced significant airway inflammation (cell debris) and had an attenuated growth phenotype in the human AT-II cells. These data extend our understanding of early Francisella infection by demonstrating that Francisella enter significant numbers of AT-II cells within the lung and that the capsule and LPS of wild type Schu S4 helps prevent murine lung damage during infection. Furthermore, our data identified that human AT-II cells allow growth of Schu S4, but these same cells supported poor growth of the attenuated LVS strain in vitro. Collectively, these data further our understanding of the role of AT-II cells in Francisella infections.

In vivo lung morphometry with hyperpolarized (3) He diffusion MRI: reproducibility and the role of diffusion-sensitizing gradient direction

PURPOSE

Lung morphometry with hyperpolarized gas diffusion MRI is a highly sensitive technique for the noninvasive measurement of acinar microstructural parameters traditionally only accessible by histology. The goal of this work is to establish the reproducibility of these measurements in healthy volunteers and their dependence on the direction of the applied diffusion-sensitizing gradient.

METHODS

Hyperpolarized helium-3 ((3) He) lung morphometry MRI was performed on a total of five healthy subjects. Two subjects received duplicate imaging on the same day and three subjects received duplicate imaging after a 4-month or 27-month delay to assess reproducibility. Four subjects repeated the measurement during the same session with different diffusion-sensitizing gradient directions to determine the effect on the parameter estimates.

RESULTS

The (3) He lung morphometry measurements were reproducible over the short term and long term (e.g., % coefficient of variation [CV] of mean chord length, Lm = 2.1% and 2.9%, respectively) and across different diffusion gradient directions (Lm % CV = 2.6%). Results also show independence of field inhomogeneity effects at 1.5T.

CONCLUSION

(3) He lung morphometry is a reproducible technique for measuring acinar microstructure and is effectively independent of the choice of diffusion gradient direction. This provides confidence for the use of this technique to compare populations and treatment efficacy.

Animal models of bronchopulmonary dysplasia. The preterm baboon models

Much of the progress in improved neonatal care, particularly management of underdeveloped preterm lungs, has been aided by investigations of multiple animal models, including the neonatal baboon (Papio species). In this article we highlight how the preterm baboon model at both 140 and 125 days gestation (term equivalent 185 days) has advanced our understanding and management of the immature human infant with neonatal lung disease. Not only is the 125-day baboon model extremely relevant to the condition of bronchopulmonary dysplasia but there are also critical neurodevelopmental and other end-organ pathological features associated with this model not fully discussed in this limited forum. We also describe efforts to incorporate perinatal infection into these preterm models, both fetal and neonatal, and particularly associated with Ureaplasma/Mycoplasma organisms. Efforts to rekindle the preterm primate model for future evaluations of therapies such as stem cell replacement, early lung recruitment interventions coupled with noninvasive surfactant and high-frequency nasal ventilation, and surfactant therapy coupled with antioxidant or anti-inflammatory medications, to name a few, should be undertaken.

Efficient estimation of the total number of acini in adult rat lung

Pulmonary airways are subdivided into conducting and gas‐exchanging airways. An acinus is defined as the small tree of gas‐exchanging airways, which is fed by the most distal purely conducting airway. Until now a dissector of five consecutive sections or airway casts were used to count acini. We developed a faster method to estimate the number of acini in young adult rats. Right middle lung lobes were critical point dried or paraffin embedded after heavy metal staining and imaged by X‐ray micro‐CT or synchrotron radiation‐based X‐rays tomographic microscopy. The entrances of the acini were counted in three‐dimensional (3D) stacks of images by scrolling through them and using morphological criteria (airway wall thickness and appearance of alveoli). Segmentation stopper were placed at the acinar entrances for 3D visualizations of the conducting airways. We observed that acinar airways start at various generations and that one transitional bronchiole may serve more than one acinus. A mean of 5612 (±547) acini per lung and a mean airspace volume of 0.907 (±0.108) μL per acinus were estimated. In 60‐day‐old rats neither the number of acini nor the mean acinar volume did correlate with the body weight or the lung volume.

An efficient method to estimate the number of acini in young adult rats has been developed. All entrances of the acini were counted, labeled, and visualized in three‐dimensional stacks of X‐ray tomographic images. We observed that acinar airways start at various generations and that one transitional bronchiole may serve more than one acinus. A mean of 5612 (±547) acini per lung and a mean acinar airspace volume of 0.907 (±0.108) µL were estimated.

Commentary on "The influence of lung airways branching structure and diffusion time on measurements and models of short-range 3He gas MR diffusion"

In a recently published paper by Parra-Robles and Wild, the authors challenge the in vivo lung morphometry technique (based on hyperpolarized gas diffusion MRI) developed by our Washington University research group. In this Commentary we demonstrate that the main conclusion of Parra-Robles and Wild, that our MRI-based lung morphometry technique "produces inaccurate estimates of the airway dimensions", does not have any scientific basis and is not in agreement with the considerable body of peer-reviewed scientific reports as well as with Parra-Robles and Wild's own data. On the contrary, our technique has a strong theoretical background, is validated, and provides accurate 3D tomographic information on lung microstructural parameters previously available only from invasive biopsy specimens. This technique has already produced a number of results related to lung morphology and function that were not previously available. In our Commentary we also discuss a number of other incorrect statements in and shortcomings of Parra-Robles and Wild's paper.

Response to Commentary on "The influence of lung airways branching structure and diffusion time on measurements and models of short-range 3He gas MR diffusion"

Our extensive investigation of the cylinder model theory through numerical modelling and purpose-designed experiments has demonstrated that it does produce inaccurate estimates of airway dimensions at all diffusion times currently used. This is due to a variety of effects: incomplete treatment of non-Gaussian effects, finite airway size, branching geometry, background susceptibility gradients and diffusion time dependence of the (3)He MR diffusion behaviour in acinar airways. The cylinder model is a good starting point for the development of a lung morphometry technique from (3)He diffusion MR but its limitations need to be understood and documented in the interest of reliable clinical interpretation.

Longitudinal, noninvasive monitoring of compensatory lung growth in mice after pneumonectomy via (3)He and (1)H magnetic resonance imaging

In rodents and some other mammals, partial pneumonectomy (PNX) of adult lungs results in rapid compensatory lung growth. In the past, quantification of compensatory lung growth and realveolarization could only be accomplished after killing the animal, removal of lungs, and histologic analysis of lungs at single time points. Hyperpolarized (3)He diffusion magnetic resonance imaging (MRI) allows in vivo morphometry of human lungs; our group has adapted this technique for application to mouse lungs. Through imaging, we can obtain maps of lung microstructural parameters that allow quantification of morphometric and physiologic measurements. In this study, we employed our (3)He MRI technique to image in vivo morphometry at baseline and to serially assess compensatory growth after left PNX of mice. (1)H and hyperpolarized (3)He diffusion MRI were performed at baseline (pre-PNX), 3-days, and 30-days after PNX. Compared with the individual mouse's own baseline, MRI was able to detect and serially quantify changes in lung volume, alveolar surface area, alveolar number, and regional changes in alveolar size that occurred during the course of post-PNX lung growth. These results are consistent with morphometry measurements reported in the literature for mouse post-PNX compensatory lung growth. In addition, we were also able to serially assess and quantify changes in the physiologic parameter of lung compliance during the course of compensatory lung growth; this was consistent with flexiVent data. With these techniques, we now have a noninvasive, in vivo method to serially assess the effectiveness of therapeutic interventions on post-PNX lung growth in the same mouse.

Multiscale imaging and registration-driven model for pulmonary acinar mechanics in the mouse

A registration-based multiscale method to obtain a deforming geometric model of mouse acinus is presented. An intact mouse lung was fixed by means of vascular perfusion at a hydrostatic inflation pressure of 20 cmH(2)O. Microcomputed tomography (μCT) scans were obtained at multiple resolutions. Substructural morphometric analysis of a complete acinus was performed by computing a surface-to-volume (S/V) ratio directly from the 3D reconstruction of the acinar geometry. A geometric similarity is observed to exist in the acinus where S/V is approximately preserved anywhere in the model. Using multiscale registration, the shape of the acinus at an elevated inflation pressure of 25 cmH(2)O is estimated. Changes in the alveolar geometry suggest that the deformation within the acinus is not isotropic. In particular, the expansion of the acinus (from 20 to 25 cmH(2)O) is accompanied by an increase in both surface area and volume in such a way that the S/V ratio is not significantly altered. The developed method forms a useful tool in registration-driven fluid and solid mechanics studies as displacement of the alveolar wall becomes available in a discrete sense.

Imaging the interaction of atelectasis and overdistension in surfactant-depleted lungs

OBJECTIVE

Atelectasis and surfactant depletion may contribute to greater distension-and thereby injury-of aerated lung regions; recruitment of atelectatic lung may protect these regions by attenuating such overdistension. However, the effects of atelectasis (and recruitment) on aerated airspaces remain elusive. We tested the hypothesis that during mechanical ventilation, surfactant depletion increases the dimensions of aerated airspaces and that lung recruitment reverses these changes.

DESIGN

Prospective imaging study in an animal model.

SETTING

Research imaging facility.

SUBJECTS

Twenty-seven healthy Sprague Dawley rats.

INTERVENTIONS

Surfactant depletion was obtained by saline lavage in anesthetized, ventilated rats. Alveolar recruitment was accomplished using positive end-expiratory pressure and exogenous surfactant administration.

MEASUREMENTS AND MAIN RESULTS

Airspace dimensions were estimated by measuring the apparent diffusion coefficient of He, using diffusion-weighted hyperpolarized gas magnetic resonance imaging. Atelectasis was demonstrated using computerized tomography and by measuring oxygenation. Saline lavage increased atelectasis (increase in nonaerated tissue from 1.2% to 13.8% of imaged area, p < 0.001), and produced a concomitant increase in mean apparent diffusion coefficient (~33%, p < 0.001) vs. baseline; the heterogeneity of the computerized tomography signal and the variance of apparent diffusion coefficient were also increased. Application of positive end-expiratory pressure and surfactant reduced the mean apparent diffusion coefficient (~23%, p < 0.001), and its variance, in parallel to alveolar recruitment (i.e., less computerized tomography densities and heterogeneity, increased oxygenation).

CONCLUSIONS

Overdistension of aerated lung occurs during atelectasis is detectable using clinically relevant magnetic resonance imaging technology, and could be a key factor in the generation of lung injury during mechanical ventilation. Lung recruitment by higher positive end-expiratory pressure and surfactant administration reduces airspace distension.

Imaging of the mouse lung with scanning laser optical tomography (SLOT)

The current study focuses on the use of scanning laser optical tomography (SLOT) in imaging of the mouse lung ex vivo. SLOT is a highly efficient fluorescence microscopy technique allowing rapid scanning of samples of a size of several millimeters, thus enabling volumetric visualization by using intrinsic contrast mechanisms of previously fixed lung lobes. Here, we demonstrate the imaging of airways, blood vessels, and parenchyma from whole, optically cleared mouse lung lobes with a resolution down to the level of single alveoli using absorption and autofluorescence scan modes. The internal structure of the lung can then be analyzed nondestructively and quantitatively in three-dimensional datasets in any preferred planar orientation. Moreover, the procedure preserves the microscopic structure of the lung and allows for subsequent correlative histologic studies. In summary, the current study has shown that SLOT is a valuable technique to study the internal structure of the mouse lung.

Evidence for adult lung growth in humans

A 33-year-old woman underwent a right-sided pneumonectomy in 1995 for treatment of a lung adenocarcinoma. As expected, there was an abrupt decrease in her vital capacity, but unexpectedly, it increased during the subsequent 15 years. Serial computed tomographic (CT) scans showed progressive enlargement of the remaining left lung and an increase in tissue density. Magnetic resonance imaging (MRI) with the use of hyperpolarized helium-3 gas showed overall acinar-airway dimensions that were consistent with an increase in the alveolar number rather than the enlargement of existing alveoli, but the alveoli in the growing lung were shallower than in normal lungs. This study provides evidence that new lung growth can occur in an adult human.

Monitoring in vivo changes in lung microstructure with ³He MRI in Sendai virus-infected mice

Recently, a Sendai virus (SeV) model of chronic obstructive lung disease has demonstrated an innate immune response in mouse airways that exhibits similarities to the chronic airway inflammation in human chronic obstructive pulmonary disease (COPD) and asthma, but the effect on distal lung parenchyma has not been investigated. The aim of our study is to image the time course and regional distribution of mouse lung microstructural changes in vivo after SeV infection. (1)H and (3)He diffusion magnetic resonance imaging (MRI) were successfully performed on five groups of C57BL/6J mice. (1)H MR images provided precise anatomical localization and lung volume measurements. (3)He lung morphometry was implemented to image and quantify mouse lung geometric microstructural parameters at different time points after SeV infection. (1)H MR images detected the SeV-induced pulmonary inflammation in vivo; spatially resolved maps of acinar airway radius R, alveolar depth h, and mean linear intercept Lm were generated from (3)He diffusion images. The morphometric parameters R and Lm in the infected group were indistinguishable from PBS-treated mice at day 21, increased slightly at day 49, and were increased with statistical significance at day 77 (p = 0.02). Increases in R and Lm of infected mice imply that there is a modest increase in alveolar duct radius distal to airway inflammation, particularly in the lung periphery, indicating airspace enlargement after virus infection. Our results indicate that (3)He lung morphometry has good sensitivity in quantifying small microstructural changes in the mouse lung and that the Sendai mouse model has the potential to be a valid murine model of COPD.

Morphometric changes in the human pulmonary acinus during inflation

Despite decades of research into the mechanisms of lung inflation and deflation, there is little consensus about whether lung inflation occurs due to the recruitment of new alveoli or by changes in the size and/or shape of alveoli and alveolar ducts. In this study we use in vivo (3)He lung morphometry via MRI to measure the average alveolar depth and alveolar duct radius at three levels of inspiration in five healthy human subjects and calculate the average alveolar volume, surface area, and the total number of alveoli at each level of inflation. Our results indicate that during a 143 ± 18% increase in lung gas volume, the average alveolar depth decreases 21 ±5%, the average alveolar duct radius increases 7 ± 3%, and the total number of alveoli increases by 96 ± 9% (results are means ± SD between subjects; P < 0.001, P < 0.01, and P < 0.00001, respectively, via paired t-tests). Thus our results indicate that in healthy human subjects the lung inflates primarily by alveolar recruitment and, to a lesser extent, by anisotropic expansion of alveolar ducts.

Mapping of (3) He apparent diffusion coefficient anisotropy at sub-millisecond diffusion times in an elastase-instilled rat model of emphysema

Hyperpolarized (3) He gas can provide detailed anatomical maps of the macroscopic airways in the lungs (i.e., ventilation) as well as insight into the lung microstructure through the apparent diffusion coefficient. In particular, the apparent diffusion coefficient of (3) He in the lung exhibits anisotropic effects that depend on diffusion time (δ), and it has been shown to be extraordinarily sensitive to enlargement in terminal airways and alveoli associated with emphysema. In this study, the anisotropic nature of the (3) He apparent diffusion coefficient is studied in a rat model of emphysema, based on elastase instillation, specifically for δ values less than one millisecond. Longitudinal (D(L) ) and transverse (D(T) ) diffusion coefficients were mapped at δ = 360 μs and δ = 800 μs based on a cylinder model of lung structure and correlated with histological measurement of alveolar damage based on mean linear intercept (L(m) ). Whole-lung mean D(T) measured at δ = 360 μs in the elastase-instilled rat lungs (0.14 ± 0.09 cm(2) /s) demonstrated the most significant increase (p = 0.00195) compared to the sham-instilled cohort (0.06 ± 0.06 cm(2) /s) and had a strong linear correlation with L(m) (Pearson's correlation coefficient of 0.9). These results suggest that measurement of (3) He apparent diffusion coefficient anisotropy, specifically D(T) , can provide a sensitive indicator of emphysema, particularly at very short diffusion times (δ = 360 μs).

In vivo detection of acinar microstructural changes in early emphysema with (3)He lung morphometry

PURPOSE

To quantitatively characterize early emphysematous changes in the lung microstructure of current and former smokers with noninvasive helium 3 ((3)He) lung morphometry and to compare these results with the clinical standards, pulmonary function testing (PFT) and low-dose computed tomography (CT).

MATERIALS AND METHODS

This study was approved by the local institutional review board, and all subjects provided informed consent. Thirty current and former smokers, each with a minimum 30-pack-year smoking history and mild or no abnormalities at PFT, underwent (3)He lung morphometry. This technique is based on diffusion MR imaging with hyperpolarized (3)He gas and yields quantitative localized in vivo measurements of acinar airway geometric parameters, such as airway radii, alveolar depth, and number of alveoli per unit lung volume. These measurements enable calculation of standard morphometric characteristics, such as mean linear intercept and surface-to-volume ratio.

RESULTS

Noninvasive (3)He lung morphometry was used to detect alterations in acinar structure in smokers with normal PFT findings. When compared with smokers with the largest forced expiratory volume in 1 second (FEV(1)) to forced vital capacity (FVC) ratio, those with chronic obstructive pulmonary disease had significantly reduced alveolar depth (0.07 mm vs 0.13 mm) and enlarged acinar ducts (0.36 mm vs 0.3 mm). The mean alveolar geometry measurements in the healthiest subjects were in excellent quantitative agreement with literature values obtained by using invasive techniques (acinar duct radius, 0.3 mm; alveolar depth, 0.14 mm at 1 L above functional residual capacity). (3)He lung morphometry depicted greater abnormalities than did PFT and CT. No adverse events were associated with inhalation of (3)He gas.

CONCLUSION

(3)He lung morphometry yields valuable noninvasive insight into early emphysematous changes in alveolar geometry with increased sensitivity compared with conventional techniques.

Lung morphometry with hyperpolarized 129Xe: theoretical background

The (3) He lung morphometry technique, based on MRI measurements of hyperpolarized (3) He gas diffusion in lung airspaces, provides unique information on the lung microstructure at the alveolar level. In vivo 3D tomographic images of standard morphological parameters (airspace chord length, lung parenchyma surface-to-volume ratio, and number of alveoli per unit volume) can be generated from a rather short (several seconds) MRI scan. The technique is based on a theory of gas diffusion in lung acinar airways and experimental measurements of diffusion-attenuated MRI signal. The present work aims at developing the theoretical background of a similar technique based on hyperpolarized (129) Xe gas. As the diffusion coefficient and gyromagnetic ratio of (129) Xe gas are substantially different from those of (3) He gas, the specific details of the theory and experimental measurements with (129) Xe should be amended. We establish phenomenological relationships between acinar airway geometrical parameters and the diffusion-attenuated MR signal for human and small animal lungs, both normal lungs and lungs with mild emphysema. Optimal diffusion times are shown to be about 5 ms for human and 1.3 ms for small animals. The expected uncertainties in measuring main morphometrical parameters of the lungs are estimated in the framework of Bayesian probability theory.

Microscopy-based quantitative analysis of lung structure: application in diagnosis

INTRODUCTION

Stereology provides a set of methods that are appropriate for a microscopy-based quantitative assessment of lung structure. In general, the aim of stereology is to obtain information on three-dimensional structures from two-dimensional sections of these structures. The inherent impartiality of stereological principles is critical in order to meet the requirements of 'good laboratory practice'.

AREAS COVERED

This article is a systematic review of the applications of stereology to characterize pathological alterations of emphysema, fibrosis, acute lung injury and tumor grading. The reader is provided with a general overview of unbiased or design-based stereology and is provided with some examples of how these methods could be integrated into a diagnostic work-up of lung diseases in humans and animal models. The article also reviews the implications of a published statement, which defines standards for quantitative assessment of lung structure based on stereology, by the American Thoracic Society and the European Respiratory Society.

EXPERT OPINION

In view of the recently published standards for quantitative assessment of lung structure, unbiased stereological methods are strongly recommended, particularly as they provide valuable information in diagnosing lung diseases and allow a statistically valid quantitative comparison between different groups. Future developments will make the application of stereology in lung biology and pathology even more efficient. Moreover, there is also the potential for combing the principles of stereology with other imaging modalities (e.g., radiological), which will allow for non-invasive lung stereology.

Imaging lung microstructure in mice with hyperpolarized 3He diffusion MRI

Quantitative measurement of lung microstructure is of great significance in assessment of pulmonary disease, particularly in the earliest stages. The technique for MRI-based 3He lung morphometry was previously developed and validated for human lungs, and was recently extended to ex vivo mouse lungs. The technique yields accurate, quantitative information about the microstructure and geometry of acinar airways. In this study, the 3He lung morphometry technique is successfully implemented for in vivo studies of mice. Results indicate excellent agreement between in vivo morphometry via 3He MRI and microscopic morphometry after sacrifice. This opens up new avenues for application of the technique as a precise, noninvasive, in vivo biomarker of changes in lung microstructure, within various mouse models of lung disease.

³He lung morphometry technique: accuracy analysis and pulse sequence optimization

The (3)He lung morphometry technique (Yablonskiy et al., JAP, 2009), based on MRI measurements of hyperpolarized gas diffusion in lung airspaces, provides unique information on the lung microstructure at the alveolar level. 3D tomographic images of standard morphological parameters (mean airspace chord length, lung parenchyma surface-to-volume ratio, and the number of alveoli per unit lung volume) can be created from a rather short (several seconds) MRI scan. These parameters are most commonly used to characterize lung morphometry but were not previously available from in vivo studies. A background of the (3)He lung morphometry technique is based on a previously proposed model of lung acinar airways, treated as cylindrical passages of external radius R covered by alveolar sleeves of depth h, and on a theory of gas diffusion in these airways. The initial works approximated the acinar airways as very long cylinders, all with the same R and h. The present work aims at analyzing effects of realistic acinar airway structures, incorporating airway branching, physiological airway lengths, a physiological ratio of airway ducts and sacs, and distributions of R and h. By means of Monte-Carlo computer simulations, we demonstrate that our technique allows rather accurate measurements of geometrical and morphological parameters of acinar airways. In particular, the accuracy of determining one of the most important physiological parameter of lung parenchyma - surface-to-volume ratio - does not exceed several percent. Second, we analyze the effect of the susceptibility induced inhomogeneous magnetic field on the parameter estimate and demonstrate that this effect is rather negligible at B(0) ≤ 3T and becomes substantial only at higher B(0) Third, we theoretically derive an optimal choice of MR pulse sequence parameters, which should be used to acquire a series of diffusion-attenuated MR signals, allowing a substantial decrease in the acquisition time and improvement in accuracy of the results. It is demonstrated that the optimal choice represents three not equidistant b-values: b(1)=0, b(2)∼2 s/cm(2), b(3)∼8 s/cm(2).