Imaging for lung physiology: what do we wish we could measure?

Imaging drug delivery to the lungs: Methods and applications in oncology

Direct delivery to the lung via inhalation is arguably one of the most logical approaches to treat lung cancer using drugs. However, despite significant efforts and investment in this area, this strategy has not progressed in clinical trials. Imaging drug delivery is a powerful tool to understand and develop novel drug delivery strategies. In this review we focus on imaging studies of drug delivery by the inhalation route, to provide a broad overview of the field to date and attempt to better understand the complexities of this route of administration and the significant barriers that it faces, as well as its advantages. We start with a discussion of the specific challenges for drug delivery to the lung via inhalation. We focus on the barriers that have prevented progress of this approach in oncology, as well as the most recent developments in this area. This is followed by a comprehensive overview of the different imaging modalities that are relevant to lung drug delivery, including nuclear imaging, X-ray imaging, magnetic resonance imaging, optical imaging and mass spectrometry imaging. For each of these modalities, examples from the literature where these techniques have been explored are provided. Finally the different applications of these technologies in oncology are discussed, focusing separately on small molecules and nanomedicines. We hope that this comprehensive review will be informative to the field and will guide the future preclinical and clinical development of this promising drug delivery strategy to maximise its therapeutic potential.

An Inhalable Theranostic System for Local Tuberculosis Treatment Containing an Isoniazid Loaded Metal Organic Framework Fe-MIL-101-NH2-From Raw MOF to Drug Delivery System

The theranostic approach to local tuberculosis treatment allows drug delivery and imaging of the lungs for a better control and personalization of antibiotic therapy. Metal-organic framework (MOF) Fe-MIL-101-NH2 nanoparticles were loaded with isoniazid. To optimize their functionality a 2³ factorial design of spray-drying with poly(lactide-co-glycolide) and leucine was employed. Powder aerodynamic properties were assessed using a twin stage impinger based on the dose emitted and the fine particle fraction. Magnetic resonance imaging (MRI) contrast capabilities were tested on porous lung tissue phantom and ex vivo rat lungs. Cell viability and uptake studies were conducted on murine macrophages RAW 246.9. The final product showed good aerodynamic properties, modified drug release, easier uptake by macrophages in relation to raw isoniazid-MOF, and MRI contrast capabilities. Starting from raw MOF, a fully functional inhalable theranostic system with a potential application in personalized tuberculosis pulmonary therapy was developed.

Pulmonary Vascular Dynamics

The pulmonary circulation carries deoxygenated blood from the systemic veins through the pulmonary arteries to be oxygenated in the capillaries that line the walls of the pulmonary alveoli. The pulmonary circulation carries the cardiac output with a relatively low driving pressure, and so differs considerably in structure and function from the systemic circulation to maintain a low-resistance vascular system. The pulmonary circulation is often considered to be a quasi-static system in both experimental and computational studies of pulmonary perfusion and its matching to ventilation (air flow) for exchange. However, the system is highly dynamic, with cardiac output and regional perfusion changing with posture, exercise, and over time. Here we review this dynamic system, with a focus on understanding the physiology of pulmonary vascular dynamics across spatial and temporal scales, and the changes to these dynamics that are reflective of disease. © 2019 American Physiological Society. Compr Physiol 9:1081-1100, 2019.

Correlation of Regional Lung Ventilation and Gas Transfer to Red Blood Cells: Implications for Functional-Avoidance Radiation Therapy Planning

PURPOSE

To investigate the degree to which lung ventilation and gas exchange are regionally correlated, using the emerging technology of hyperpolarized (HP)-129Xe magnetic resonance imaging (MRI).

METHODS AND MATERIALS

Hyperpolarized-129Xe MRI studies were performed on 17 institutional review board-approved human subjects, including 13 healthy volunteers, 1 emphysema patient, and 3 non-small cell lung cancer patients imaged before and approximately 11 weeks after radiation therapy (RT). Subjects inhaled 1 L of HP-129Xe mixture, followed by the acquisition of interleaved ventilation and gas exchange images, from which maps were obtained of the relative HP-129Xe distribution in three states: (1) gaseous, in lung airspaces; (2) dissolved interstitially, in alveolar barrier tissue; and (3) transferred to red blood cells (RBCs), in the capillary vasculature. The relative spatial distributions of HP-129Xe in airspaces (regional ventilation) and RBCs (regional gas transfer) were compared. Further, we investigated the degree to which ventilation and RBC transfer images identified similar functional regions of interest (ROIs) suitable for functionally guided RT. For the RT patients, both ventilation and RBC functional images were used to calculate differences in the lung dose-function histogram and functional effective uniform dose.

RESULTS

The correlation of ventilation and RBC transfer was ρ = 0.39 ± 0.15 in healthy volunteers. For the RT patients, this correlation was ρ = 0.53 ± 0.02 before treatment and ρ = 0.39 ± 0.07 after treatment; for the emphysema patient it was ρ = 0.24. Comparing functional ROIs, ventilation and RBC transfer demonstrated poor spatial agreement: Dice similarity coefficient = 0.50 ± 0.07 and 0.26 ± 0.12 for the highest-33%- and highest-10%-function ROIs in healthy volunteers, and in RT patients (before treatment) these were 0.58 ± 0.04 and 0.40 ± 0.04. The average magnitude of the differences between RBC- and ventilation-derived functional effective uniform dose, fV20Gy, fV10Gy, and fV5Gy were 1.5 ± 1.4 Gy, 4.1% ± 3.8%, 5.0% ± 3.8%, and 5.3% ± 3.9%, respectively.

CONCLUSION

Ventilation may not be an effective surrogate for true regional lung function for all patients.

Regional Distribution of Pulmonary Blood Volume with Dual-Energy Computed Tomography: Results in 42 Subjects

RATIONALE AND OBJECTIVES

The noninvasive approach of lung perfusion generated from dual-energy computed tomography acquisitions has entered clinical practice. The purpose of this study was to analyze the regional distribution of iodine within distal portions of the pulmonary arterial bed on dual-source, dual-energy computed tomography examinations in a cohort of subjects without cardiopulmonary pathologies.

MATERIALS AND METHODS

The study population included 42 patients without cardiorespiratory disease, enabling quantitative and qualitative analysis of pulmonary blood volume after administration of a 40% contrast agent. Qualitative analysis was based on visual assessment. Quantitative analysis was obtained after semiautomatic division of each lung into 18 areas.

RESULTS

The iodine concentration did not significantly differ between the right (R) and left (L) lungs (P = .49), with a mean attenuation of 41.35 Hounsfield units (HU) and 41.14 HU, respectively. Three regional gradients of attenuation were observed between: (a) lung bases and apices (P < .001), linked to the conditions of examination (mean Δ: 6.23 in the R lung; 5.96 in the L lung); (b) posterior and anterior parts of the lung (P < .001) due to gravity (mean Δ: 11.92 in the R lung ; 15.93 in the L lung); and (c) medullary and cortical lung zones (P < .001) (mean Δ: 9.35 in the R lung ; 8.37 in the L lung). The intensity of dependent-nondependent (r = 0.42; P < .001) and corticomedullary (r = 0.58; P < .0001) gradients was correlated to the overall iodine concentration.

CONCLUSION

Distribution of pulmonary blood volume is influenced by physiological gradients and scanning conditions.

Three-dimensional oxygen-enhanced MRI of the human lung at 1.5T with ultra-fast balanced steady-state free precession

PURPOSE

To assess the feasibility of 3D oxygen-enhanced (OE) MRI of the lung at 1.5T using multi-volumetric ultra-fast balanced steady-state free precession (ufSSFP) acquisitions.

METHODS

Isotropic imaging of the lung for OE-MRI was performed with an adapted 3D ufSSFP sequence using five breath-hold acquisitions ranging from functional residual capacity to tidal inspiration under both normoxic (room air) and hyperoxic (100% O₂ ) gas conditions. For each O₂ concentration, a sponge model (which captures the parenchymal signal intensity variation as a function of the lung volume) was fitted to the acquired multi-volumetric datasets after semiautomatic lung segmentation and deformable image registration. From the retrieved model parameters, 3D oxygen-enhancement maps were calculated.

RESULTS

For OE ufSSFP imaging, the maximum parenchymal signal is observed for flip angles around 23° under both normoxic and hyperoxic conditions. It is found that the sponge model accurately describes parenchymal signal at different breathing positions, thereby mitigating the confounding bias in the estimated oxygen enhancement from residual density modulations. From the model, an average lung oxygen enhancement of 7.0% ± 0.3% was found in the healthy volunteers, and the oxygen-enhancement maps indicate a ventral to dorsal gravitation-related gradient.

CONCLUSION

The study demonstrates the feasibility of whole-lung OE-MRI from multi-volumetric ufSSFP in healthy volunteers. Magn Reson Med 79:246-255, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Quantitative assessment of lung stiffness in patients with interstitial lung disease using MR elastography

PURPOSE

To investigate the use of magnetic resonance elastography (MRE) in the quantitative assessment of pulmonary fibrosis by comparing quantitative shear stiffness measurements of lung parenchyma in patients diagnosed with fibrotic interstitial lung disease (ILD) and healthy controls.

MATERIALS AND METHODS

A 1.5T spin-echo, echo planar imaging MRE (SE-EPI MRE) pulse sequence was utilized to assess absolute lung shear stiffness in 15 patients with diagnosed ILD and in 11 healthy controls. Data were collected at residual volume (RV) and total lung capacity (TLC). Spirometry data were obtained immediately prior to scanning. To test for statistical significance between RV and TLC shear stiffness estimates a two-sample t-test was performed. To assess variability within individual subject shear stiffness estimates, the intraclass correlation coefficient (ICC) and Krippendorff's alpha were calculated.

RESULTS

Patients with ILD exhibited an average (±1 standard deviation) shear stiffness of 2.74 (±0.896) kPa at TLC and 1.32 (±0.300) kPa at RV. The corresponding values for healthy individuals were 1.33 (±0.195) kPa and 0.849 (±0.250) kPa, respectively. The difference in shear stiffness between RV and TLC was statistically significant (P < 0.001). At TLC, the ICC and alpha values were 0.909 and 0.887, respectively. At RV, the ICC and alpha values were 0.852 and 0.862, respectively.

CONCLUSION

In subjects with known fibrotic interstitial lung disease, parenchymal shear stiffness is increased when compared to normal controls at both RV and TLC, with TLC demonstrating the most significant difference. MRE-derived parenchymal shear stiffness is a promising new noninvasive imaging-based biomarker of interstitial lung disease.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 2 J. MAGN. RESON. IMAGING 2017;46:365-374.

Correlation of 68Ga Ventilation-Perfusion PET/CT with Pulmonary Function Test Indices for Assessing Lung Function

Pulmonary function tests (PFTs) are routinely used to assess lung function, but they do not provide information about regional pulmonary dysfunction. We aimed to assess correlation of quantitative ventilation-perfusion (V/Q) PET/CT with PFT indices.

METHODS

Thirty patients underwent V/Q PET/CT and PFT. Respiration-gated images were acquired after inhalation of (68)Ga-carbon nanoparticles and administration of (68)Ga-macroaggregated albumin. Functional volumes were calculated by dividing the volume of normal ventilated and perfused (%NVQ), unmatched and matched defects by the total lung volume. These functional volumes were correlated with forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC), FEV1/FVC, and diffusing capacity for carbon monoxide (DLCO).

RESULTS

All functional volumes were significantly different in patients with chronic obstructive pulmonary disease (P < 0.05). FEV1/FVC and %NVQ had the highest correlation (r = 0.82). FEV1 was also best correlated with %NVQ (r = 0.64). DLCO was best correlated with the volume of unmatched defects (r = -0.55). Considering %NVQ only, a cutoff value of 90% correctly categorized 28 of 30 patients with or without significant pulmonary function impairment.

CONCLUSION

Our study demonstrates strong correlations between V/Q PET/CT functional volumes and PFT parameters. Because V/Q PET/CT is able to assess regional lung function, these data support the feasibility of its use in radiation therapy and preoperative planning and assessing pulmonary dysfunction in a variety of respiratory diseases.

Imaging regional lung function: a critical tool for developing inhaled antimicrobial therapies

Alterations in regional lung function due to respiratory infection have a significant effect on the deposition of inhaled treatments. This has consequences for treatment effectiveness and hence recovery of lung function. In order to advance our understanding of respiratory infection and inhaled treatment delivery, we must develop imaging techniques that can provide regional functional measurements of the lung. In this review, we explore the role of functional imaging for the assessment of respiratory infection and development of inhaled treatments. We describe established and emerging functional lung imaging methods. The effect of infection on lung function is described, and the link between regional disease, function, and inhaled treatments is discussed. The potential for lung function imaging to provide unique insights into the functional consequences of infection, and its treatment, is also discussed.

3D MRI of impaired hyperpolarized 129Xe uptake in a rat model of pulmonary fibrosis

A variety of pulmonary pathologies, in particular interstitial lung diseases, are characterized by thickening of the pulmonary blood-gas barrier, and this thickening results in reduced gas exchange. Such diffusive impairment is challenging to quantify spatially, because the distributions of the metabolically relevant gases (CO2 and O2) cannot be detected directly within the lungs. Hyperpolarized (HP) (129)Xe is a promising surrogate for these metabolic gases, because MR spectroscopy and imaging allow gaseous alveolar (129)Xe to be detected separately from (129)Xe dissolved in the red blood cells (RBCs) and the adjacent tissues, which comprise blood plasma and lung interstitium. Because (129)Xe reaches the RBCs by diffusing across the same barrier tissues (blood plasma and interstitium) as O2, barrier thickening will delay (129)Xe transit and, thus, reduce RBC-specific (129)Xe MR signal. Here we have exploited these properties to generate 3D, MR images of (129)Xe uptake by the RBCs in two groups of rats. In the experimental group, unilateral fibrotic injury was generated prior to imaging by instilling bleomycin into one lung. In the control group, a unilateral sham instillation of saline was performed. Uptake of (129)Xe by the RBCs, quantified as the fraction of RBC signal relative to total dissolved (129)Xe signal, was significantly reduced (P = 0.03) in the injured lungs of bleomycin-treated animals. In contrast, no significant difference (P = 0.56) was observed between the saline-treated and untreated lungs of control animals. Together, these results indicate that 3D MRI of HP (129)Xe dissolved in the pulmonary tissues can provide useful biomarkers of impaired diffusive gas exchange resulting from fibrotic thickening.

Automated measurement of heterogeneity in CT images of healthy and diseased rat lungs using variogram analysis of an octree decomposition

Background

Assessing heterogeneity in lung images can be an important diagnosis tool. We present a novel and objective method for assessing lung damage in a rat model of emphysema. We combined a three-dimensional (3D) computer graphics method–octree decomposition–with a geostatistics-based approach for assessing spatial relationships–the variogram–to evaluate disease in 3D computed tomography (CT) image volumes.

Methods

Male, Sprague-Dawley rats were dosed intratracheally with saline (control), or with elastase dissolved in saline to either the whole lung (for mild, global disease) or a single lobe (for severe, local disease). Gated 3D micro-CT images were acquired on the lungs of all rats at end expiration. Images were masked, and octree decomposition was performed on the images to reduce the lungs to homogeneous blocks of 2 × 2 × 2, 4 × 4 × 4, and 8 × 8 × 8 voxels. To focus on lung parenchyma, small blocks were ignored because they primarily defined boundaries and vascular features, and the spatial variance between all pairs of the 8 × 8 × 8 blocks was calculated as the square of the difference of signal intensity. Variograms–graphs of distance vs. variance–were constructed, and results of a least-squares-fit were compared. The robustness of the approach was tested on images prepared with various filtering protocols. Statistical assessment of the similarity of the three control rats was made with a Kruskal-Wallis rank sum test. A Mann-Whitney-Wilcoxon rank sum test was used to measure statistical distinction between individuals. For comparison with the variogram results, the coefficient of variation and the emphysema index were also calculated for all rats.

Results

Variogram analysis showed that the control rats were statistically indistinct (p = 0.12), but there were significant differences between control, mild global disease, and severe local disease groups (p < 0.0001). A heterogeneity index was calculated to describe the difference of an individual variogram from the control average. This metric also showed clear separation between dose groups. The coefficient of variation and the emphysema index, on the other hand, did not separate groups.

Conclusion

These results suggest the octree decomposition and variogram analysis approach may be a rapid, non-subjective, and sensitive imaging-based biomarker for characterizing lung disease.

Estimation of the absolute shear stiffness of human lung parenchyma using (1) H spin echo, echo planar MR elastography

PURPOSE

To develop a rapid proton MR elastography (MRE) technique that can quantify the absolute shear stiffness of lung parenchyma, to investigate the ability to differentiate respiration-dependent stiffness variations of the lung, and to demonstrate clinical feasibility.

MATERIALS AND METHODS

A spin-echo echo planar imaging MRE sequence (SE-EPI MRE) with a very short echo time was developed and tested in a series of five healthy volunteers at three different lung volumes: (i) residual volume (RV), (ii) total lung capacity (TLC), (iii) and midway between RV and TLC (MID). At each volume, lung density was quantified using a MR-based density mapping sequence. For reference, data were acquired using the previously described spin-echo lung MRE sequence (SE-MRE). MRE data were also acquired in a patient with proven Idiopathic Pulmonary Fibrosis (IPF) to test clinical feasibility.

RESULTS

The SE-EPIMRE sequence reduced total acquisition time by a factor of 2 compared with the SE-MRE sequence. Lung parenchyma median shear stiffness for the 5 volunteers quantified with the SE-EPI MRE sequence was 0.9 kPa, 1.1 kPa, and 1.6 kPa at RV, MID, and TLC, respectively. The corresponding values obtained with the SE-MRE sequence were 0.9 kPa, 1.1 kPa, and 1.5 kPa. Absolute shear stiffness was also successfully measured in the IPF patient.

CONCLUSION

The results indicate that stiffness variations due to respiration could be measured with the SE-EPIMRE technique and were equivalent to values generated by the previously described SE-MRE approach. Preliminary data obtained from the patient demonstrate clinical feasibility.