Morphometric differences between central vs. surface acini in A/J mice using high-resolution micro-computed tomography

The fix is not yet in: recommendation for fixation of lungs within physiological/pathophysiological volume range in preclinical pulmonary structure-function studies

Quantitative characterization of lung structures by morphometrical or stereological analysis of histological sections is a powerful means of elucidating pulmonary structure-function relations. The overwhelming majority of studies, however, fix lungs for histology at pressures outside the physiological/pathophysiological respiratory volume range. Thus, valuable information is being lost. In this perspective article, we argue that investigators performing pulmonary histological studies should consider whether the aims of their studies would benefit from fixation at functional transpulmonary pressures, particularly those of end-inspiration and end-expiration. We survey the pressures at which lungs are typically fixed in preclinical structure-function studies, provide examples of conditions that would benefit from histological evaluation at functional lung volumes, summarize available fixation methods, discuss alternative imaging modalities, and discuss challenges to implementing the suggested approach and means of addressing those challenges. We aim to persuade investigators that modifying or complementing the traditional histological approach by fixing lungs at minimal and maximal functional volumes could enable new understanding of pulmonary structure-function relations.

Mapping the strain-stiffening behavior of the lung and lung cancer at microscale resolution using the crystal ribcage

Lung diseases such as cancer substantially alter the mechanical properties of the organ with direct impact on the development, progression, diagnosis, and treatment response of diseases. Despite significant interest in the lung's material properties, measuring the stiffness of intact lungs at sub-alveolar resolution has not been possible. Recently, we developed the crystal ribcage to image functioning lungs at optical resolution while controlling physiological parameters such as air pressure. Here, we introduce a data-driven, multiscale network model that takes images of the lung at different distending pressures, acquired via the crystal ribcage, and produces corresponding absolute stiffness maps. Following validation, we report absolute stiffness maps of the functioning lung at microscale resolution in health and disease. For representative images of a healthy lung and a lung with primary cancer, we find that while the lung exhibits significant stiffness heterogeneity at the microscale, primary tumors introduce even greater heterogeneity into the lung's microenvironment. Additionally, we observe that while the healthy alveoli exhibit strain-stiffening of ∼1.75 times, the tumor's stiffness increases by a factor of six across the range of measured transpulmonary pressures. While the tumor stiffness is 1.4 times the lung stiffness at a transpulmonary pressure of three cmH2O, the tumor's mean stiffness is nearly five times greater than that of the surrounding tissue at a transpulmonary pressure of 18 cmH2O. Finally, we report that the variance in both strain and stiffness increases with transpulmonary pressure in both the healthy and cancerous lungs. Our new method allows quantitative assessment of disease-induced stiffness changes in the alveoli with implications for mechanotransduction.

The pathogenesis of influenza in intact alveoli: virion endocytosis and its effects on the lung's air-blood barrier

Lung infection by influenza A virus (IAV) is a major cause of global mortality from lung injury, a disease defined by widespread dysfunction of the lung's air-blood barrier. Endocytosis of IAV virions by the alveolar epithelium - the cells that determine barrier function - is central to barrier loss mechanisms. Here, we address the current understanding of the mechanistic steps that lead to endocytosis in the alveolar epithelium, with an eye to how the unique structure of lung alveoli shapes endocytic mechanisms. We highlight where future studies of alveolar interactions with IAV virions may lead to new therapeutic approaches for IAV-induced lung injury.

Crystal ribcage: a platform for probing real-time lung function at cellular resolution

Understanding the dynamic pathogenesis and treatment response in pulmonary diseases requires probing the lung at cellular resolution in real time. Despite advances in intravital imaging, optical imaging of the lung during active respiration and circulation has remained challenging. Here, we introduce the crystal ribcage: a transparent ribcage that allows multiscale optical imaging of the functioning lung from whole-organ to single-cell level. It enables the modulation of lung biophysics and immunity through intravascular, intrapulmonary, intraparenchymal and optogenetic interventions, and it preserves the three-dimensional architecture, air-liquid interface, cellular diversity and respiratory-circulatory functions of the lung. Utilizing these capabilities on murine models of pulmonary pathologies we probed remodeling of respiratory-circulatory functions at the single-alveolus and capillary levels during disease progression. The crystal ribcage and its broad applications presented here will facilitate further studies of nearly any pulmonary disease as well as lead to the identification of new targets for treatment strategies.

Quantitative Imaging Metrics for the Assessment of Pulmonary Pathophysiology: An Official American Thoracic Society and Fleischner Society Joint Workshop Report

Multiple thoracic imaging modalities have been developed to link structure to function in the diagnosis and monitoring of lung disease. Volumetric computed tomography (CT) renders three-dimensional maps of lung structures and may be combined with positron emission tomography (PET) to obtain dynamic physiological data. Magnetic resonance imaging (MRI) using ultrashort-echo time (UTE) sequences has improved signal detection from lung parenchyma; contrast agents are used to deduce airway function, ventilation-perfusion-diffusion, and mechanics. Proton MRI can measure regional ventilation-perfusion ratio. Quantitative imaging (QI)-derived endpoints have been developed to identify structure-function phenotypes, including air-blood-tissue volume partition, bronchovascular remodeling, emphysema, fibrosis, and textural patterns indicating architectural alteration. Coregistered landmarks on paired images obtained at different lung volumes are used to infer airway caliber, air trapping, gas and blood transport, compliance, and deformation. This document summarizes fundamental "good practice" stereological principles in QI study design and analysis; evaluates technical capabilities and limitations of common imaging modalities; and assesses major QI endpoints regarding underlying assumptions and limitations, ability to detect and stratify heterogeneous, overlapping pathophysiology, and monitor disease progression and therapeutic response, correlated with and complementary to, functional indices. The goal is to promote unbiased quantification and interpretation of in vivo imaging data, compare metrics obtained using different QI modalities to ensure accurate and reproducible metric derivation, and avoid misrepresentation of inferred physiological processes. The role of imaging-based computational modeling in advancing these goals is emphasized. Fundamental principles outlined herein are critical for all forms of QI irrespective of acquisition modality or disease entity.

Mathematical modeling of pulmonary acinus structure: Verification of acinar shape effects on pathway structure using rat lungs

The pulmonary acinus is the gas exchange unit in the lung and has a very complex microstructure. The structure model is essential to understand the relationship between structural heterogeneity and mechanical phenomena at the acinus level with computational approaches. We propose an acinus structure model represented by a cluster of truncated octahedra in conical, double-conical, inverted conical, or chestnut-like conical confinement to accommodate recent experimental information of rodent acinar shapes. The basis of the model is the combined use of Voronoi and Delaunay tessellations and the optimization of the ductal tree assuming the number of alveoli and the mean path length as quantities related to gas exchange. Before applying the Voronoi tessellation, controlling the seed coordinates enables us to model acinus with arbitrary shapes. Depending on the acinar shape, the distribution of path length varies. The lengths are more widely spread for the cone acinus, with a bias toward higher values, while most of the lengths for the inverted cone acinus primarily take a similar value. Longer pathways have smaller tortuosity and more generations, and duct length per generation is almost constant irrespective of generation, which agrees well with available experimental data. The pathway structure of cone and chestnut-like cone acini is similar to the surface acini's features reported in experiments. According to space-filling requirements in the lung, other conical acini may also be acceptable. The mathematical acinus structure model with various conical shapes can be a platform for computational studies on regional differences in lung functions along the lung surface, underlying respiratory physiology and pathophysiology.

Origins of and lessons from quantitative functional X-ray computed tomography of the lung

Functional CT of the lung has emerged from quantitative CT (qCT). Structural details extracted at multiple lung volumes offer indices of function. Additionally, single volumetric images, if acquired at standardized lung volumes and body posture, can be used to model function by employing such engineering techniques as computational fluid dynamics. With the emergence of multispectral CT imaging including dual energy from energy integrating CT scanners and multienergy binning using the newly released photon counting CT technology, function is tagged via use of contrast agents. Lung disease phenotypes have previously been lumped together by the limitations of spirometry and plethysmography. QCT and its functional embodiment have been imbedded into studies seeking to characterize chronic obstructive pulmonary disease, severe asthma, interstitial lung disease and more. Reductions in radiation dose by an order of magnitude or more have been achieved. At the same time, we have seen significant increases in spatial and density resolution along with methodologic validations of extracted metrics. Together, these have allowed attention to turn towards more mild forms of disease and younger populations. In early applications, clinical CT offered anatomic details of the lung. Functional CT offers regional measures of lung mechanics, the assessment of functional small airways disease, as well as regional ventilation-perfusion matching (V/Q) and more. This paper will focus on the use of quantitative/functional CT for the non-invasive exploration of dynamic three-dimensional functioning of the breathing lung and beating heart within the unique negative pressure intrathoracic environment of the closed chest.

Recent advances in the understanding of alveolar flow

Understanding the dynamics of airflow in alveoli and its effect on the behavior of particle transport and deposition is important for understanding lung functions and the cause of many lung diseases. The studies on these areas have drawn substantial attention over the last few decades. This Review discusses the recent progress in the investigation of behavior of airflow in alveoli. The information obtained from studies on the structure of the lung airway tree and alveolar topology is provided first. The current research progress on the modeling of alveoli is then reviewed. The alveolar cell parameters at different generation of branches, issues to model real alveolar flow, and the current numerical and experimental approaches are discussed. The findings on flow behavior, in particular, flow patterns and the mechanism of chaotic flow generation in the alveoli are reviewed next. The different flow patterns under different geometrical and flow conditions are discussed. Finally, developments on microfluidic devices such as lung-on-a-chip devices are reviewed. The issues of current devices are discussed.

Novel Subtypes of Pulmonary Emphysema Based on Spatially-Informed Lung Texture Learning: The Multi-Ethnic Study of Atherosclerosis (MESA) COPD Study

Pulmonary emphysema overlaps considerably with chronic obstructive pulmonary disease (COPD), and is traditionally subcategorized into three subtypes previously identified on autopsy. Unsupervised learning of emphysema subtypes on computed tomography (CT) opens the way to new definitions of emphysema subtypes and eliminates the need of thorough manual labeling. However, CT-based emphysema subtypes have been limited to texture-based patterns without considering spatial location. In this work, we introduce a standardized spatial mapping of the lung for quantitative study of lung texture location and propose a novel framework for combining spatial and texture information to discover spatially-informed lung texture patterns (sLTPs) that represent novel emphysema subtype candidates. Exploiting two cohorts of full-lung CT scans from the MESA COPD (n = 317) and EMCAP (n = 22) studies, we first show that our spatial mapping enables population-wide study of emphysema spatial location. We then evaluate the characteristics of the sLTPs discovered on MESA COPD, and show that they are reproducible, able to encode standard emphysema subtypes, and associated with physiological symptoms.

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.

Ultra-high-resolution computed tomography shows changes in the lungs related with airway hyperresponsiveness in a murine asthma model

In vivo presentation of airway hyper-responsiveness (AHR) at the different time points of the allergic reaction is not clearly understood. The purpose of this study was to investigate how AHR manifests in the airway and the lung parenchyma in vivo following exposure to different stimuli and in the early and late phases of asthma after allergen exposure. Ovalbumin (OVA)-induced allergic asthma model was established using 6-week female BALB/c mice. Enhanced pause was measured with a non-invasive method to assess AHR. The dynamic changes of the airway and lung parenchyma were evaluated with ultra-high-resolution computed tomography (128 multi-detector, 1024 × 1024 matrix) for 10 h. While the methacholine challenge showed no grossly visible changes in the proximal airway and lung parenchyma despite provoking AHR, the OVA challenge induced significant immediate changes manifesting as peribronchial ground glass opacities, consolidations, air-trapping, and paradoxical proximal airway dilatations. After resolution of immediate response, multiple episodes of AHRs occurred with paradoxical proximal airway dilatation and peripheral air-trapping in late phase over a prolonged time period in vivo. Understanding of airflow limitation based on the structural changes of asthmatic airway would be helpful to make an appropriate drug delivery strategy for the treatment of asthma.

A semi-automated method for unbiased alveolar morphometry: Validation in a bronchopulmonary dysplasia model

Reproducible and unbiased methods to quantify alveolar structure are important for research on many lung diseases. However, manually estimating alveolar structure through stereology is time consuming and inter-observer variability is high. The objective of this work was to develop and validate a fast, reproducible and accurate (semi-)automatic alternative. A FIJI-macro was designed that automatically segments lung images to binary masks, and counts the number of test points falling on tissue and the number of intersections of the air-tissue interface with a set of test lines. Manual selection remains necessary for the recognition of non-parenchymal tissue and alveolar exudates. Volume density of alveolar septa (VVsep) and mean linear intercept of the airspaces (Lm) as measured by the macro were compared to theoretical values for 11 artificial test images and to manually counted values for 17 lungs slides using linear regression and Bland-Altman plots. Inter-observer agreement between 3 observers, measuring 8 lungs both manually and automatically, was assessed using intraclass correlation coefficients (ICC). VVsep and Lm measured by the macro closely approached theoretical values for artificial test images (R² of 0.9750 and 0.9573 and bias of 0.34% and 8.7%). The macro data in lungs were slightly higher for VVsep and slightly lower for Lm in comparison to manually counted values (R² of 0.8262 and 0.8288 and bias of -6.0% and 12.1%). Visually, semi-automatic segmentation was accurate. Most importantly, manually counted VVsep and Lm had only moderate to good inter-observer agreement (ICC 0.859 and 0.643), but agreements were excellent for semi-automatically counted values (ICC 0.956 and 0.900). This semi-automatic method provides accurate and highly reproducible alveolar morphometry results. Future efforts should focus on refining methods for automatic detection of non-parenchymal tissue or exudates, and for assessment of lung structure on 3D reconstructions of lungs scanned with microCT.

Micrometer-resolution X-ray tomographic full-volume reconstruction of an intact post-mortem juvenile rat lung

In this article, we present an X-ray tomographic imaging method that is well suited for pulmonary disease studies in animal models to resolve the full pathway from gas intake to gas exchange. Current state-of-the-art synchrotron-based tomographic phase-contrast imaging methods allow for three-dimensional microscopic imaging data to be acquired non-destructively in scan times of the order of seconds with good soft tissue contrast. However, when studying multi-scale hierarchically structured objects, such as the mammalian lung, the overall sample size typically exceeds the field of view illuminated by the X-rays in a single scan and the necessity for achieving a high spatial resolution conflicts with the need to image the whole sample. Several image stitching and calibration techniques to achieve extended high-resolution fields of view have been reported, but those approaches tend to fail when imaging non-stable samples, thus precluding tomographic measurements of large biological samples, which are prone to degradation and motion during extended scan times. In this work, we demonstrate a full-volume three-dimensional reconstruction of an intact rat lung under immediate post-mortem conditions and at an isotropic voxel size of (2.75 µm)³. We present the methodology for collecting multiple local tomographies with 360° extended field of view scans followed by locally non-rigid volumetric stitching. Applied to the lung, it allows to resolve the entire pulmonary structure from the trachea down to the parenchyma in a single dataset. The complete dataset is available online (https://doi.org/10.16907/7eb141d3-11f1-47a6-9d0e-76f8832ed1b2).

Heterogeneous structure and surface tension effects on mechanical response in pulmonary acinus: A finite element analysis

BACKGROUND

The pulmonary acinus is a dead-end microstructure that consists of ducts and alveoli. High-resolution micro-CT imaging has recently provided detailed anatomical information of a complete in vivo acinus, but relating its mechanical response with its detailed acinar structure remains challenging. This study aimed to investigate the mechanical response of acinar tissue in a whole acinus for static inflation using computational approaches.

METHODS

We performed finite element analysis of a whole acinus for static inflation. The acinar structure model was generated based on micro-CT images of an intact acinus. A continuum mechanics model of the lung parenchyma was used for acinar tissue material model, and surface tension effects were explicitly included. An anisotropic mechanical field analysis based on a stretch tensor was combined with a curvature-based local structure analysis.

FINDINGS

The airspace of the acinus exhibited nonspherical deformation as a result of the anisotropic deformation of acinar tissue. A strain hotspot occurred at the ridge-shaped region caused by a rod-like deformation of acinar tissue on the ridge. The local structure becomes bowl-shaped for inflation and, without surface tension effects, the surface of the bowl-shaped region primarily experiences isotropic deformation. Surface tension effects suppressed the increase in airspace volume and inner surface area, while facilitating anisotropic deformation on the alveolar surface.

INTERPRETATION

In the lungs, the heterogeneous acinar structure and surface tension induce anisotropic deformation at the acinar and alveolar scales. Further research is needed on structural variation of acini, inter-acini connectivity, or dynamic behavior to understand multiscale lung mechanics.

MRI and CT lung biomarkers: Towards an in vivo understanding of lung biomechanics

BACKGROUND

The biomechanical properties of the lung are necessarily dependent on its structure and function, both of which are complex and change over time and space. This makes in vivo evaluation of lung biomechanics and a deep understanding of lung biomarkers, very challenging. In patients and animal models of lung disease, in vivo evaluations of lung structure and function are typically made at the mouth and include spirometry, multiple-breath gas washout tests and the forced oscillation technique. These techniques, and the biomarkers they provide, incorporate the properties of the whole organ system including the parenchyma, large and small airways, mouth, diaphragm and intercostal muscles. Unfortunately, these well-established measurements mask regional differences, limiting their ability to probe the lung's gross and micro-biomechanical properties which vary widely throughout the organ and its subcompartments. Pulmonary imaging has the advantage in providing regional, non-invasive measurements of healthy and diseased lung, in vivo. Here we summarize well-established and emerging lung imaging tools and biomarkers and how they may be used to generate lung biomechanical measurements.

METHODS

We review well-established and emerging lung anatomical, microstructural and functional imaging biomarkers generated using synchrotron x-ray tomographic-microscopy (SRXTM), micro-x-ray computed-tomography (micro-CT), clinical CT as well as magnetic resonance imaging (MRI).

FINDINGS

Pulmonary imaging provides measurements of lung structure, function and biomechanics with high spatial and temporal resolution. Imaging biomarkers that reflect the biomechanical properties of the lung are now being validated to provide a deeper understanding of the lung that cannot be achieved using measurements made at the mouth.

How high resolution 3-dimensional imaging changes our understanding of postnatal lung development

During the last 10 + years biologically and clinically significant questions about postnatal lung development could be answered due to the application of modern cutting-edge microscopic and quantitative histological techniques. These are in particular synchrotron radiation based X-ray tomographic microscopy (SRXTM), but also ³Helium Magnetic Resonance Imaging, as well as the stereological estimation of the number of alveoli and the length of the free septal edge. First, the most important new finding may be the following: alveolarization of the lung does not cease after the maturation of the alveolar microvasculature but continues until young adulthood and, even more important, maybe reactivated lifelong if needed to rescue structural damages of the lungs. Second, the pulmonary acinus represents the functional unit of the lung. Because the borders of the acini could not be detected in classical histological sections, any investigation of the acini requires 3-dimensional (imaging) methods. Based on SRXTM it was shown that in rat lungs the number of acini stays constant, meaning that their volume increases by a factor of ~ 11 after birth. The latter is very important for acinar ventilation and particle deposition.

Biomimetics of the pulmonary environment in vitro: A microfluidics perspective

The entire luminal surface of the lungs is populated with a complex yet confluent, uninterrupted airway epithelium in conjunction with an extracellular liquid lining layer that creates the air-liquid interface (ALI), a critical feature of healthy lungs. Motivated by lung disease modelling, cytotoxicity studies, and drug delivery assessments amongst other, in vitro setups have been traditionally conducted using macroscopic cultures of isolated airway cells under submerged conditions or instead using transwell inserts with permeable membranes to model the ALI architecture. Yet, such strategies continue to fall short of delivering a sufficiently realistic physiological in vitro airway environment that cohesively integrates at true-scale three essential pillars: morphological constraints (i.e., airway anatomy), physiological conditions (e.g., respiratory airflows), and biological functionality (e.g., cellular makeup). With the advent of microfluidic lung-on-chips, there have been tremendous efforts towards designing biomimetic airway models of the epithelial barrier, including the ALI, and leveraging such in vitro scaffolds as a gateway for pulmonary disease modelling and drug screening assays. Here, we review in vitro platforms mimicking the pulmonary environment and identify ongoing challenges in reconstituting accurate biological airway barriers that still widely prevent microfluidic systems from delivering mainstream assays for the end-user, as compared to macroscale in vitro cell cultures. We further discuss existing hurdles in scaling up current lung-on-chip designs, from single airway models to more physiologically realistic airway environments that are anticipated to deliver increasingly meaningful whole-organ functions, with an outlook on translational and precision medicine.

Capturing complexity in pulmonary system modelling

Respiratory disease is a significant problem worldwide, and it is a problem with increasing prevalence. Pathology in the upper airways and lung is very difficult to diagnose and treat, as response to disease is often heterogeneous across patients. Computational models have long been used to help understand respiratory function, and these models have evolved alongside increases in the resolution of medical imaging and increased capability of functional imaging, advances in biological knowledge, mathematical techniques and computational power. The benefits of increasingly complex and realistic geometric and biophysical models of the respiratory system are that they are able to capture heterogeneity in patient response to disease and predict emergent function across spatial scales from the delicate alveolar structures to the whole organ level. However, with increasing complexity, models become harder to solve and in some cases harder to validate, which can reduce their impact clinically. Here, we review the evolution of complexity in computational models of the respiratory system, including successes in translation of models into the clinical arena. We also highlight major challenges in modelling the respiratory system, while making use of the evolving functional data that are available for model parameterisation and testing.

Development of the lung

To fulfill the task of gas exchange, the lung possesses a huge inner surface and a tree-like system of conducting airways ventilating the gas exchange area. During lung development, the conducting airways are formed first, followed by the formation and enlargement of the gas exchange area. The latter (alveolarization) continues until young adulthood. During organogenesis, the left and right lungs have their own anlage, an outpouching of the foregut. Each lung bud starts a repetitive process of outgrowth and branching (branching morphogenesis) that forms all of the future airways mainly during the pseudoglandular stage. During the canalicular stage, the differentiation of the epithelia becomes visible and the bronchioalveolar duct junction is formed. The location of this junction stays constant throughout life. Towards the end of the canalicular stage, the first gas exchange may take place and survival of prematurely born babies becomes possible. Ninety percent of the gas exchange surface area will be formed by alveolarization, a process where existing airspaces are subdivided by the formation of new walls (septa). This process requires a double-layered capillary network at the basis of the newly forming septum. However, in parallel to alveolarization, the double-layered capillary network of the immature septa fuses to a single-layered network resulting in an optimized setup for gas exchange. Alveolarization still continues, because, at sites where new septa are lifting off preexisting mature septa, the required second capillary layer will be formed instantly by angiogenesis. The latter confirms a lifelong ability of alveolarization, which is important for any kind of lung regeneration.