Next Article in Journal
Altered Phenotype and Enhanced Antibody-Producing Ability of Peripheral B Cells in Mice with Cd19-Driven Cre Expression
Next Article in Special Issue
A Focused Review on Primary Graft Dysfunction after Clinical Lung Transplantation: A Multilevel Syndrome
Previous Article in Journal
Going Micro in Leptospirosis Kidney Disease
Previous Article in Special Issue
A Comprehensive Review on the Surgical Aspect of Lung Transplant Models in Mice and Rats
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantitative CT Correlates with Local Inflammation in Lung of Patients with Subtypes of Chronic Lung Allograft Dysfunction

1
Department of Radiology, University of Michigan, Ann Arbor, MI 48109, USA
2
Antwerp Surgical Training, Anatomy and Research Centre (ASTARC), Faculty of Medicine and Health Sciences, University of Antwerp, Wilrijk, 2610 Antwerp, Belgium
3
Department of Chronic Diseases and Metabolism (CHROMETA), Katholieke Universiteit Leuven, 3000 Leuven, Belgium
4
Division of Pulmonary and Critical Care Medicine, University of Michigan, Ann Arbor, MI 48109, USA
5
Department of Biostatistics, School of Public Health, University of Michigan, Ann Arbor, MI 48109, USA
6
Respiratory Division, University Hospital Leuven, 3000 Leuven, Belgium
7
Imbio, LLC, Minneapolis, MN 55405, USA
*
Author to whom correspondence should be addressed.
Cells 2022, 11(4), 699; https://doi.org/10.3390/cells11040699
Submission received: 1 December 2021 / Revised: 8 February 2022 / Accepted: 11 February 2022 / Published: 16 February 2022
(This article belongs to the Special Issue Advances in Lung Transplantation)

Abstract

:
Chronic rejection of lung allografts has two major subtypes, bronchiolitis obliterans syndrome (BOS) and restrictive allograft syndrome (RAS), which present radiologically either as air trapping with small airways disease or with persistent pleuroparenchymal opacities. Parametric response mapping (PRM), a computed tomography (CT) methodology, has been demonstrated as an objective readout of BOS and RAS and bears prognostic importance, but has yet to be correlated to biological measures. Using a topological technique, we evaluate the distribution and arrangement of PRM-derived classifications of pulmonary abnormalities from lung transplant recipients undergoing redo-transplantation for end-stage BOS (N = 6) or RAS (N = 6). Topological metrics were determined from each PRM classification and compared to structural and biological markers determined from microCT and histopathology of lung core samples. Whole-lung measurements of PRM-defined functional small airways disease (fSAD), which serves as a readout of BOS, were significantly elevated in BOS versus RAS patients (p = 0.01). At the core-level, PRM-defined parenchymal disease, a potential readout of RAS, was found to correlate to neutrophil and collagen I levels (p < 0.05). We demonstrate the relationship of structural and biological markers to the CT-based distribution and arrangement of PRM-derived readouts of BOS and RAS.

1. Introduction

Transplanted lungs have one of the highest rejection rates among all solid organ transplantations, and approximately 50% of lung transplant recipients suffer from chronic lung allograft dysfunction (CLAD) within five years of transplantation [1,2]. The main subtypes of CLAD are (a) bronchiolitis obliterans syndrome (BOS), which is characterized by persistent spirometric airflow obstruction and mostly air trapping on expiratory X-ray computed tomography (CT), with or without bronchiectasis, and (b) restrictive allograft syndrome (RAS), which is defined by decreased total lung capacity ≤90% of baseline and persistent radiographic opacities [3,4,5]. The underlying pathology of both subtypes of CLAD is irreversible fibroproliferation, either strictly airway-related in BOS, or involving both the airways and the parenchyma in RAS. In RAS, fibroproliferation leads to fibrotic remodeling of the lung, which is preceded by immune-mediated graft injury. Once persistent lung function loss occurs, prognosis is poor, with most patients demonstrating ongoing decline [6]. Of note, prognosis is significantly better for BOS than for RAS patients (3–5 years post-diagnosis versus 6–18 months).
We previously demonstrated using a combination of CT and microCT that obstruction of pre-terminal airways is common in both BOS and RAS (albeit significantly higher in more distal airways in BOS), while in RAS there is also parenchymal destruction leading to further distortion and remodeling of terminal bronchioles [7,8]. Lungs with end-stage RAS were characterized by a more pronounced inflammatory response with increased presence of neutrophils, eosinophils, mast cells, B cells and cytotoxic T cells compared to controls. In contrast, lungs with end-stage BOS only showed increased T and B cells, along with increased neutrophilia around the airways, compared to controls [9].
Parametric response mapping (PRM) is a quantitative CT-based algorithm that, when applied to paired inspiratory and expiratory high-resolution computed tomography (HRCT), has been shown to simultaneously detect parenchymal abnormalities associated with emphysema and small airways disease (SAD) in patients with chronic obstructive pulmonary disease (COPD) [10,11]. The ability of this technique to quantify SAD in COPD has been validated using microCT analysis of core samples from explanted lungs [12]. This study was the first to demonstrate the correlation of PRM-derived SAD to areas of lung tissue with terminal bronchial loss, luminal narrowing and obstruction as determined by microCT. In addition, PRM has been extended as a readout for other lung diseases, offering insights into pulmonary complications associated with chronic rejection from hematopoietic stem cell and lung transplantations [13,14,15,16]. This CT-based technique has recently been demonstrated to also offer valuable prognostic information for potential CLAD, identifying patients most at risk of developing CLAD or dying [17].
In these studies, PRM classifications, inherently three-dimensional (3D) datasets of pulmonary abnormalities, have been reported simply as percentages of lung volume. We have expanded the clinical application of PRM by demonstrating how topological analysis of PRM (tPRM) provides information on the distribution and arrangement of PRM classifications that correlates to clinically meaningful measures of COPD severity [18]. In this study, we extend tPRM to analysis of CLAD, as we investigate the relationship between topological features of PRM-derived classifications and structural and biological markers obtained from core samples extracted from explants of lung transplant recipients diagnosed with BOS and RAS who underwent re-transplantation or autopsy.

2. Materials and Methods

2.1. Patient Cohort

All patients who underwent double lung transplant reached a best postoperative forced expiratory volume at 1 s percent predicted (FEV1pp) > 80% and received azithromycin treatment for potential CLAD but were found to be nonresponsive. BOS and RAS were diagnosed according to the most recent International Society for Heart and Lung Transplantation guidelines. In brief, BOS patients showed a decline in FEV1pp > 20% compared to baseline but did not show evidence of concomitant restriction (no total lung capacity pp (TLCpp) < 90% or forced vital capacity pp (FVCpp) < 80% of baseline), without evidence of persistent opacities on CT. In contrast, RAS patients showed a simultaneous decrease in FEV1pp and TLCpp and persistent pleuroparenchymal opacities. Inspiratory and expiratory CT scans were acquired from all patients and the last CT scan prior to redo transplantation was used for further analysis. Clinical diagnosis was confirmed with pathological analysis of the contralateral lung with typical interstitial and (sub)pleural fibrosis in RAS and mostly inconspicuous parenchyma and airway obliteration in BOS. This study was approved by the local hospital’s ethical committee (S57752).

2.2. Computed Tomography

CT data were obtained as whole lung volumetric CT scans at full inspiration (TLC) and at relaxed expiration (functional residual capacity (FRC)) on Somatom scanner (Siemens, Munich, Germany) and reconstructed using a b60 or b70 reconstruction kernel. Slice thickness was 1 mm for all scans. All CT scans were checked for Hounsfield unit (HU) drift and if necessary corrected based on aortic blood (50 HU) and central air (1000 HU) as previously described [19].

2.3. Parametric Response Map

PRM was applied to all paired inspiratory and expiratory CT scans from both study groups. Briefly, lungs from paired scans were segmented from the thoracic cavity using an in-house algorithm written in MATLAB (The MathWorks, Inc., Natick, MA, USA). The whole-lung inspiration CT scan was spatially registered and aligned to the CT scan obtained at expiration using a deformable image registration algorithm provided in Elastix, an open-source image registration software package [20,21]. Once complete, the images share the same geometric space where each voxel, the smallest unit of volume in a 3D image dataset, consists of HU values at inspiration and expiration. To minimize the effects of noise on the PRM analysis, a 3 × 3 × 3 median filter was applied to both spatially aligned CT scans. Individual voxels were then classified based on predetermined thresholds as normal (color coded green), functional small airways disease (fSAD; color coded yellow), emphysema (color coded red), and parenchymal disease (color coded magenta). Voxels with values ≥ −950 HU and <−810 HU at inspiration and ≥−856 HU at expiration were classified normal, ≥−950 HU and <−810 HU at inspiration and <−856 HU at expiration were fSAD, <−950 HU at inspiration and <856 HU at expiration were emphysema, and ≥−810 HU at inspiration were parenchymal disease (PD) [13]. Global PRM measures were presented as relative lung volumes.

2.4. Topological Analysis of PRM (tPRM)

Topological properties of each PRM classification map were defined in this study through the Minkowski measures (local estimates of the Minkowski functionals) associated with 3D distributions: Volume (V, in mm3), Surface Area (S, in mm2), Mean Breadth (B, in mm), and the Euler–Poincaré characteristic (χ) [22]. A detailed description of these parameters is provided in the Supplementary Materials in [18]. Maps of Minkowski measures (i.e., V, S, B, and χ) were computed using a moving window of size 213 voxels evaluated on a grid of every 5th voxel. V, S, and B were normalized by the Minkowski estimate of the mask within the same local window volume (rather than a direct calculation of the mask volume in the window as previously described) and χ by the masked window voxel count. Minkowski measures were quantified per subject as the mean local normalized value over the entire lung volume for group comparisons and regression. Final displayed representations of spatially resolved indices have been linearly interpolated back to original dimensions. It is important to note that the volume density of a given PRM classification is proportional to relative volumes of PRM classification. All image processing was performed using in-house algorithms developed in a technical computing language (MATLAB, The MathWorks Inc., Natick, MA, USA).

2.5. Sampling of the Ex Vivo Lung

Since 2009, explant lungs have been prospectively collected from all lung transplant patients undergoing lung redo transplantation or at autopsy (IRB approval UZ/KU Leuven S52174). Permission for CT analysis was requested separately from the UZ/KU Leuven ethics committee via approval number S57752. Explant lungs were cannulated and inflated to total lung capacity by applying a transpulmonary pressure of 30 cm H2O. They were then deflated to 10 cm H2O transpulmonary pressure, followed by freezing in the fumes of liquid nitrogen vapor. Lungs were scanned using HRCT (Siemens Definition Flash 120 kV; 110 mAs; Siemens, Munich, Germany). Following CT, the lungs were cut with a band saw, from apex to base, into slices of 2 cm and systematically sampled. Photographs were taken before and after sampling to ensure proper matching. Four cores per lung were selected for further analysis. These cores were chosen to present the different spectra of disease in the lungs. Specifically, for RAS patients, we included both fibrotic and non-fibrotic areas.

2.6. Topological PRM Analysis of Lung Cores

Topological features of individual cores were determined through spatial alignment of the in vivo paired CT scans to a photograph of the explanted cored lung section. This method was previously reported [12]. The process involved the following five steps. First, the pre-transplant inspiration CT scan was spatially aligned to the paired expiration CT scans using the Elastix deformable registration algorithm. Second, the lung used for the microCT analysis was segmented from the thoracic cavity in the pre-transplant expiration CT scan and spatially aligned to the CT scan of the explanted inflated frozen lung using the same deformable registration algorithm. The third step, registering the inflated lung to the uncored section RGB photograph, required additional post-processing. The photograph of the uncored section was converted to gray scale and subsampled by a factor of 10. A histogram equalization algorithm was applied to improve image contrast. To obtain the 3D topological maps for each PRM classification within the core, the image of the uncored lung section was replicated such that a 3D data set was constructed with the same number of slices as the inflated frozen lung 3D CT scan. The uncored lung section image on the center slice was then segmented. A similar transformation, i.e., rotate-translate-isotropic scaling, was performed to spatially align the inflated lung CT scan to the segmented uncored section in the dataset. Fourth, the uncored lung section was spatially aligned to the cored lung section using post-processing steps described in step three. Finally, the transformation matrices from each step were applied to the topological maps for each PRM classification so that each map (i.e., V, S, B and χ) was aligned to the cored section. Mean values for the topological features were calculated for the individual cores. The core volume of interest was determined by contouring the cored volume in the 3D cored lung section dataset. The number of contoured slices was determined by dividing the section thickness by the transformed voxel size (35 slices = 20 mm section (0.82 scaling factor × 0.7 mm slice thickness for inflated lung CT scan)).

2.7. MicroCT Imaging and Analysis of the Ex Vivo Lung Cores

The frozen cores were scanned in frozen condition with microCT (Skyscan 1172; Bruker, Kontich, Belgium) using a cooling stage (resolution 10 µm; 40 kV, 226 mA, 0.5° rotation steps). After scanning, the images were reconstructed with NRecon (Bruker, Kontich, Belgium) and processed using CT for assessment of the parenchymal morphometry (tissue percentage of core and surface density) and by manually counting the number of terminal bronchioles (TB/mL) [23]. A portion of the core was vacuum embedded in OCT (Sakura, Tokyo, Japan), and cryosections (8 µm) were made. Sections were stained for hematoxylin and picrosirius red (total collagen). Antibodies were used to detect collagen type I (ab6308, Abcam, Cambridge, UK) and collagen type III (2150-0100, Bio-Rad, Hercules, CA, USA). Detailed methods for histology labeling and processing steps are available in [24]. To assess the presence of relevant immune cells, immunohistochemical stains for T-lymphocytes (CD4, CD8), B-lymphocytes (CD20), neutrophils (MPO), eosinophils (ECP), mast cells (CD117), and macrophages (CD68) were made using an AEC chromogen. Slides were digitalized using an Aperio digital pathology slide scanner (Leica Biosystems Inc., Wetzlar, Germany) and analyzed using Qupath (version v0.3.2) and ImageJ (version 1.53o). Details on the methods used for immunohistochemical staining protocols are provided in Supplementary Table S1 and the analysis is described in [9]. The relative number of positively stained cells was calculated by normalizing stained cells to all identified cells.

2.8. Statistical Analysis

All data are presented as the mean ± the standard deviation unless stated otherwise. Differences in metrics between CLAD subtypes were assessed by Mann–Whitney U Test for all continuous variables and Fisher Exact test for categorical variables. Correlation analysis was determined using Spearman rho. All statistical computations were performed using a statistical software package (SPSS Software Products, version 28.0.1). Results were considered statistically significant at the two-sided 5% comparison-wise significance level (p < 0.05).

3. Results

3.1. Study Population

Patient characteristics are presented in Table 1. There were no significant differences in recipient age, gender, height, and weight. Subtle differences were observed in underlying disease between groups. FEV1pp (20.3 ± 4.0 and 30.2 ± 10.7 for BOS and RAS, respectively) and FEV1/FVC (0.42 ± 0.08 and 0.72 ± 0.20 for BOS and RAS, respectively) were found to be significantly lower in BOS than in RAS patients. Mean lung density as measured on CT at expiration was significantly lower in BOS patients (−806 ± 55 HU) compared with RAS patients (−678 ± 112 HU).

3.2. Representative Cases

Provided in Figure 1 are CT coronal views of the thoracic region with lungs at full inspiration and expiration for a patient with BOS and one with RAS. The BOS case is a male diagnosed with BOS 206 days post-transplant. This patient had a FEV1pp of 17% and FEV1/FVC of 0.4 when undergoing second transplantation. The RAS case is a female diagnosed with RAS 2649 days post-transplant. This patient had a FEV1pp of 50% and FEV1/FVC of 0.88 when undergoing second transplantation. The corresponding volume density maps for PRM-derived normal parenchyma (green), fSAD (yellow) and PD (magenta) are also presented. The lung transplant recipient who developed BOS had more PRM-derived fSAD distributed throughout the lungs, whereas the patient who developed RAS had more PRM-derived PD (Figure 1).

3.3. Group Comparison of tPRM Whole Lung Analysis

Volume density, surface area, and Euler–Poincaré characteristic metrics associated with PRM-derived fSAD were found to differ significantly between BOS and RAS (Figure 2). The volume density of fSAD, a readout of disease extent (VfSAD), was on average more than double in BOS patients (0.48 ± 0.14) than in RAS patients (0.17 ± 0.13). The surface area, a readout of disease distribution (SfSAD), was also found to be elevated in the BOS group (0.66 ± 0.08 mm−1) compared to RAS (0.37 ± 0.25 mm−1), suggesting most patients have a diffuse spread of BOS throughout the lungs. Euler–Poincaré characteristic (χfSAD), a readout of disease pockets or holes, was negative for the BOS group (−0.007 ± 0.004). Negative values in χfSAD occur when PRM-derived fSAD is the dominant classification within the lungs with holes that may comprise of other PRM classes, or structures such as airways and blood vessels. In contrast, the observed mean value of χfSAD for RAS patients was positive (0.003 ± 0.001), which is associated with small pockets of PRM-derived fSAD. Although the mean curvature of PRM-derived fSAD, a readout of tunneling within this classification, was higher in RAS (0.04 ± 0.01 mm−2) than in BOS patients (0.02 ± 0.02 mm−2), this was not found to be significant.
In contrast to PRM-derived fSAD, the topology of PD did not result in significant differences between groups for any of the metrics. Notwithstanding, the volume density and surface area of PRM-derived PD generated expected trends. For VPD and SPD, patients with RAS (0.34 ± 0.19 and 0.53 ± 0.18 mm−1, respectively) had higher values compared to patients with BOS (0.20 ± 0.03 (p = 0.145) and 0.36 ± 0.12 mm−1 (p = 0.091), respectively). Near identical values in mean curvature were observed between groups (BOS at 0.04 ± 0.01 mm−2 and RAS at 0.04 ± 0.03 mm−2). Although the mean value of χPD was lower in RAS (0.0003 ± 0.006) than in BOS (0.003 ± 0.001) patients, no significant difference (p = 0.132) was observed due to the large spread in values for RAS patients.
PRM-derived normal and emphysema provided mixed results. Only mean curvature for PRM-derived normal parenchyma showed a significant difference between groups (Figure 2). PRM-derived emphysema was found to generate significant differences between BOS and RAS for surface area, mean curvature and Euler–Poincaré characteristic. The volume density of PRM-derived emphysema on average was around 0.05 or less for all cases (0.05 is equivalent to 5% emphysema of total lung volume). This is partly attributed to scanner noise and airways that were not adequately segmented from the lung volume.

3.4. Correlation of Whole-Lung tPRM Metrics and PFT

Based on whole-lung analysis of individual tPRM metrics, the volume density of PRM-derived fSAD and PD were elevated for BOS and RAS, respectively. Using the total study population (N = 12), correlations to pulmonary function measurements were assessed for all tPRM metrics of these two PRM classifications (Table 2). FVC and FEV1/FVC were strongly correlated to multiple tPRM metrics, except mean curvature (B) for fSAD. As PRM-derived fSAD serves as a readout of BOS, patients with BOS typically manifest higher levels of FVC and fSAD than their RAS counterparts. The strongest correlation was observed for FVC against χfSAD (R = −0.848, p-value < 0.001). Here, patients with BOS had higher FVC and negative χfSAD, whereas RAS patients, except for one, had low FVC and positive χfSAD. VfSAD and SfSAD demonstrated significant correlations for FVC and FEV1/FVC (Table 2). For PRM-derived PD, only the surface area (S) was found to result in significant correlates (FVC: R = −0.734, p-value = 0.006 and FEV1/FVC: R = 0.615, p-value = 0.033). Nevertheless, VPD and χPD demonstrated mild correlations to FVC with p-values equal to 0.07 and 0.08, respectively. This trend was not observed when correlated to FEV1/FVC (p-values = 0.175 for VPD and 0.152 for χPD).

3.5. tPRM Metrics from Cored Samples

Figure 3 illustrates how tPRM metrics for individual cores depend on the location of the lung sample (core location provided on explant (blue circle) and CT scan (yellow circle)). This case is a male diagnosed with BOS 447 days post-transplant. This patient had a FEV1pp of 18% and FEV1/FVC of 0.41 when undergoing second transplantation. The distribution and arrangement of each PRM-derived classification generates topological measures that are spatially dependent on the underlying condition of the lung parenchyma. This lung section, approximately midway along the vertical axis (head to foot), comprised primarily of fSAD (yellow on volume density map) and normal parenchyma (green on volume density map). Regions with extreme values in volume density (low and high) resulted in low values of surface area (blue on surface area maps). Lung regions with a diffuse distribution of PRM classes resulted in higher values in S (red on surface area maps). This is clear for SNorm where normal lung parenchyma regionally coexists with local parenchyma classified as fSAD and PD (Figure 3 SNorm). The low values in SNorm were a consequence of that entire lung region being classified as fSAD. The Euler–Poincaré characteristic provides an indication of the dominance of any PRM class, as well as its transition to or from another PRM class. As seen for PRM-derived fSAD, χfSAD generated negative values where VfSAD was highest. Within that same lung region, χNorm had positive values, suggesting an inferior contribution to this space (red regions). These trends are consistent with the whole lung analysis for χNorm and χfSAD in the BOS group (Figure 2).

3.6. tPRM Metrics of Individual Cores

Topological values averaged over each core are presented for all cores sampled from each patient (Figure 4). As PRM-derived fSAD and PD are readouts of BOS and RAS, respectively, only these PRM classes are provided. Five out of the six patients diagnosed with BOS were found to have cores with VfSAD > VPD, with median values of VfSAD around or greater than 0.5 (i.e., 50% of the core volume was identified as fSAD). Only one case (231) was found to have similar contributions of PRM-derived fSAD (median VfSAD = 0.24) and PD (median VPD = 0.34) within the core sample. The surface area (S) and Euler-Poincaré characteristic (χ) were also consistent among BOS cases, except for 231. Unlike the BOS patients, no consistent topology of fSAD and PD were identified from the core samples obtained from lung sections of RAS patients. Two of the six RAS patients had core samples with elevated levels of VPD with median values above 0.6. Three of the six (cases 133, 193, and 234) had similar volume densities for fSAD and PD. The sixth case (143) had higher values of VfSAD than VPD, like trends observed in BOS patients.

3.7. MicroCT and Histopathology Analysis

Figure 5 presents microCT and immunohistochemistry stains of representative cores from patients diagnosed with BOS (Figure 5A,B) and RAS (Figure 5C,D). Morphologically inconspicuous lung parenchyma can be seen on microCT with segmental airway scarring and obstruction (Figure 5A). Representative slides for all stains for each case and a control are provided in Supplementary Figure S1. In contrast, the microCT of the RAS patient demonstrated loss of alveolar space, opacities, and extensive scarring of the interstitium (Figure 5C). Infiltration of MPO-expressing cells was evident in the airways for the BOS patient (Figure 5B). Collagen I (COL I) deposits were identified in the interstitial space of the RAS patient (Figure 5D). Group comparison of the microCT measures found significant differences in the tissue percentage of the cores and surface density (Table 3). The BOS group was found to have significantly lower levels of CD4+ T cells than the RAS group (Table 3). Although not significant (p = 0.06), neutrophil levels were more elevated in BOS than RAS subjects.

3.8. Correlation of tPRM Metrics to Structural and Biological Markers

For PRM classes fSAD and PD, correlations were determined for topology metrics and structural and biological markers (Table 4). This analysis was repeated using pooled core data (N = 45) and core data from each group (N = 23 for BOS and N = 22 for RAS). Combining all cores from each group resulted in significant correlations (yellow text in Table 4 indicates p-value < 0.05) of tissue percentage of cores (T%) and surface density (SD) to nearly all the topology measures for PRM-derived fSAD. The strongest correlation was observed between VfSAD and T% (R = −0.599, p < 0.001). PRM-derived PD was also found to demonstrate significant correlations with the volume density (VPD) to T% (R = 0.457, p = 0.002) and the number of terminal bronchioles per mL (TB/mL) (R = 0.380, p = 0.010). MPO levels were found to correlate significantly to the volume density and surface area (S) of both PRM classes, as well as to the Euler–Poincaré characteristic of PD (χPD). CD117 was found to correlate with V and S of fSAD and COL I with χPD.
Correlations for individual groups generated different results from the pooled data. No significant correlations were observed within the BOS group. Nevertheless, non-significant correlations for BOS subjects were seen for VPD and TB/mL (R = 0.363, p = 0.089). Cores obtained from RAS patients showed significant correlations between tPRM metrics and structural and biological markers. Like the pooled data, TB/mL, MPOs, and COL I were found to significantly correlate with various topology metrics from both PRM classes. The strongest correlation, irrespective of pooling, was observed between VPD and MPO levels, with a R of −0.684 and p-value < 0.001.

4. Discussion

In this study, we confirmed the role of PRM-derived fSAD and PD as CT-based readouts of BOS and RAS, respectively, in lung transplant recipients. We also demonstrated the relationship between the extent, distribution, and arrangement of our PRM-derived classifications of pulmonary abnormalities, and the structural and biological microenvironment of lung parenchyma in these same CLAD subtypes. We provide further support for our findings by demonstrating elevated levels of fSAD and PD in BOS and RAS patients, respectively, and a correlation between topological measures of fSAD and PD and microCT-measured surface density and MPO levels in cored lung samples.
An important find in our study was confirmation that PRM-derived classifications are strongly associated with clinically relevant measures of CLAD subtypes. Consistent with prior publications by both our groups [13,14,17,25] and researchers from Stanford University [15], the extent of PRM-derived fSAD, as measured by its volume density, was observed to be a readout of BOS. In our study, VfSAD was found to be on average around 0.48 ± 0.18 of the total lung volume. This is equivalent to saying 48% of the lungs are PRM-derived fSAD. This value was higher than the ~30% reported previously, which has been demonstrated as a link to prognosis [14,15,16]. As described in the Methods section, the mean of VfSAD over the lungs is equal to the percentage of PRM-derived fSAD normalized to 100. Taking this into account, we attribute the differences in the extent of fSAD reported here and in earlier studies to the severity of the disease and not to differences in calculation. Unlike the previous studies where PRM was determined at the time of diagnosis, patients in our cohort were scheduled for second lung transplants due to end-stage CLAD. With respect to the extent of PRM-derived PD, elevated values found in RAS patients (Figure 2) were consistent with previously published work [14,15,17]. Despite this consistency, wide variability in measurements occurs across the literature, likely due to differences in the time of measurement as part of the study designs, as well as to the inherent variability in the degree and location of fibrosis which has been previously discussed [26].
Case-specific and spatial variability in topological features within cores were evident in our study population. Between the cohorts, the lung transplant recipients diagnosed with BOS had similar topological values of volume density (V), surface area (S) and Euler–Poincaré characteristic (χ) for PRM-derived fSAD and PD. Out of the 6 BOS cases, 5 were found to have large differences in the median values over all cores per case in fSAD and PD. The cores from these cases had extensive small airways disease (VfSAD > VPD) that was highly diffuse (elevated SfSAD) and was the dominant pathology within the cores (i.e., fSAD) (negative χfSAD) (Figure 4). This is consistent with our observation over whole lungs (Figure 2), where significant differences were only observed between BOS and RAS for V, S and χ. A large range of topology values between case-specific cores was observed, as depicted by the whiskers in Figure 4. These findings suggest local variability in the distribution and arrangement of disease subtypes within individual cases.
A unique trend in topological features could not be identified in lung transplant recipients diagnosed with RAS. Only cases 184 and 233, as seen in Figure 4, generated elevated VPD and SPD and negative χPD, which is like the topology of PRM-derived fSAD in BOS patients. Further, case 143, who had radiographically-identified scarring, presented several similarities consistent with the BOS group. Among these were trends in V, S and χ of the core samples, and volume densities of 0.37 for fSAD and 0.11 for PD when averaged over the lungs. This patient also had FVC (L) and FEV1/FVC of 1.95 L and 0.68, respectively, again consistent with BOS. It is important to note that PRM is defined within a range of −1000 HU to −250 HU on both the inspiration and expiration CT scans. Unlike air trapping due to BOS, pulmonary fibrosis, a hallmark of RAS, may have densities well outside this range. As such, PRM only detects and quantifies parenchyma with mild opacities that may be attributed to pneumonitis [15] and fibrosis [14], which is the likely explanation for the observed discrepancy.
Topology measures for PRM-derived fSAD and PD demonstrated strong correlations to microCT and biological measures when CLAD subtypes were pooled (N = 45). The most pronounced correlations for microCT and biology to topologies were identified in percentage of core tissue and MPOs (Table 4). These findings changed when analyzing cores from specific CLAD subtypes. All significant correlations observed in the pooled analysis were lost in the BOS cohort. We attribute the lack of correlation in our topology measures to microCT and biological measures to the limited range in disease severity between the BOS cases. As shown in Figure 4 and discussed above, 5 out of 6 cases had very similar values in topological features for both PRM-derived disease classifications (i.e., fSAD and PD). In contrast, the RAS group presented significant correlations. Most notable was surface area of PRM-derived fSAD, which correlated with terminal bronchiole count/mL and MPO levels, and volume density of PRM-derived PD, which was found to correlate with MPO and COL I levels (Table 4). The retention of these significant finds from the pooled analysis is attributed to wide variability in disease subtypes found in the RAS cohort, which provided more dynamic range in CT, microCT and biological measures. We suspect that infiltration of MPO-expressing cells occurs prior to collagen deposits associated with interstitial scarring. This would explain not only the negative correlation between VPD and MPOs (R = −0.684, p < 0.001) but also the positive correlation between VPD and COL I levels (R = 0.478, p = 0.028).
Our study was subject to some limitations worth addressing. First, the number of cases (N = 12) accrued as part of this study limited our statistical power. Nevertheless, whole-lung analysis of our cohort generated volume densities of PRM-derived fSAD and PD (Figure 2) consistent with PRM results from earlier studies, specifically elevated fSAD and PD for BOS and RAS, respectively [14,15,17]. Another limitation of the small number of cases was that core measurements, around 3–4 cores per case, were treated as independent measures. Though variations between cases were observed, especially within the RAS group, large intra-case variability occurred as identified by large whiskers in Figure 4. We also acknowledge that we could not mark 100% of all cells for proper quantification as could be done by FACS analysis. However, the current study does resemble and very closely approaches the global picture of the presence and amount of the different immune cells and collagen. Finally, we recognize the potential errors associated with spatially aligning clinical volumetric CT scans to images of explanted lung sections. The process requires three transformations, each producing some degree of misalignment. Nevertheless, we have previously demonstrated this approach validates PRM-derived fSAD as a measure of bronchiolitis in COPD patients [12]. Like all experimental techniques, great care was taken to minimize variability in final measurements while providing a consistent and reproducible technique to correlate clinical CT scans to microCT and biological core measurements.

5. Conclusions

We have offered a view of the PRM landscape in a group of lung transplant recipients who developed specific CLAD subtypes. We identified a connection between the structural and biological makeup of the lung microenvironment and the topology of a clinically relevant CT-based readout of BOS and RAS. Our investigative group is working on developing additional techniques to elucidate the nature of small airway and parenchymal abnormality identified by PRM topology. These steps could further confirm the potential use and benefit of PRM for monitoring CLAD progression in lung transplant recipients.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cells11040699/s1, Figure S1: Representative Slides for Each Stain, Table S1: Antibodies and Staining Protocols.

Author Contributions

Conceptualization, C.J.G. and S.E.V.; methodology, C.J.G., B.A.H. and S.E.V.; software, B.A.H., S.R. and A.J.B.; validation, C.J.G., A.J.B. and S.R.; formal analysis, S.E.V., S.R., C.J.G. and S.M.; investigation, S.R., C.J.G. and S.E.V.; resources, C.J.G. and S.E.V.; data curation, S.E.V., S.R., A.J.B., B.M.V., R.V. and G.M.V.; writing—original draft preparation, S.R., C.J.G. and S.E.V.; writing—review and editing, M.K.H., S.G., W.W.L., V.N.L., C.R.H., E.A.K., S.M., C.J.G., S.R. and S.E.V.; visualization, C.J.G. and S.E.V.; supervision, C.J.G., S.E.V., B.M.V., R.V. and G.M.V.; project administration, C.J.G., S.E.V. and B.M.V.; funding acquisition, C.J.G. and V.N.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Heart, Lung, and Blood Institute/NIH, grant number R01HL139690 and the Cystic Fibrosis Foundation, grant number LAMA20AB0. The APC was funded by R01HL139690.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Ethics Committee of UZ/KU Leuven (approval number S57752). The acquisition of explant lungs for this study was approved by the Institutional Review Board of UZ/KU Leuven (approval number S52174). CT analysis was approved by the UZ/KU Leuven Ethics Committee (approval number S57752).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available in order to protect the privacy of the subjects of this study.

Acknowledgments

The authors thank Lee Olsen for assisting with manuscript preparation and editing.

Conflicts of Interest

C.J.G. is co-inventor of PRM, which the University of Michigan has licensed to Imbio, LLC. B.A.H. and C.J.G. are co-inventors of tPRM, which the University of Michigan has licensed to Imbio, LLC. C.J.G. and B.A.H. have a financial interest in Imbio, LLC. C.R.H. is employed by Imbio, LLC. W.W.L. reports personal fees from Konica Minolta and Continuing Education Alliance. M.K.H. reports personal fees from GlaxoSmithKline, AstraZeneca, Boehringer Ingelheim, Cipla, Chiesi, Novartis, Pulmonx, Teva, Verona, Merck, Mylan, Sanofi, DevPro, Aerogen, Polarian, Regeneron, United Therapeutics, Medscape, and Integrity. She has received either in kind research support or funds paid to the institution from the National Institutes of Health (NIH), Novartis, Sunovion, Nuvaira, Sanofi, AstraZeneca, Boehringer Ingelheim, Gala Therapeutics, Biodesix, the COPD Foundation, and the American Lung Association. She has participated in Data Safety Monitoring Boards for Novartis and Medtronic with funds paid to the institution. She has received stock options from Meissa Vaccines. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Verleden, G.M.; Raghu, G.; Meyer, K.C.; Glanville, A.R.; Corris, P. A new classification system for chronic lung allograft dysfunction. J. Heart Lung Transplant. 2014, 33, 127–133. [Google Scholar] [CrossRef] [PubMed]
  2. Yusen, R.D.; Edwards, L.B.; Kucheryavaya, A.Y.; Benden, C.; Dipchand, A.I.; Goldfarb, S.B.; Levvey, B.J.; Lund, L.H.; Meiser, B.; Rossano, J.W.; et al. The Registry of the International Society for Heart and Lung Transplantation: Thirty-second Official Adult Lung and Heart-Lung Transplantation Report--2015; Focus Theme: Early Graft Failure. J. Heart Lung Transplant. 2015, 34, 1264–1277. [Google Scholar] [CrossRef] [PubMed]
  3. Verleden, G.M.; Glanville, A.R.; Lease, E.D.; Fisher, A.J.; Calabrese, F.; Corris, P.A.; Ensor, C.R.; Gottlieb, J.; Hachem, R.R.; Lama, V.; et al. Chronic lung allograft dysfunction: Definition, diagnostic criteria, and approaches to treatment-A consensus report from the Pulmonary Council of the ISHLT. J. Heart Lung Transplant. 2019, 38, 493–503. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Meyer, K.C.; Raghu, G.; Verleden, G.M.; Corris, P.A.; Aurora, P.; Wilson, K.C.; Brozek, J.; Glanville, A.R.; ISHLT/ATS/ERS BOS Task Force Committee. An international ISHLT/ATS/ERS clinical practice guideline: Diagnosis and management of bronchiolitis obliterans syndrome. Eur. Respir. J. 2014, 44, 1479–1503. [Google Scholar] [CrossRef] [PubMed]
  5. Glanville, A.R.; Verleden, G.M.; Todd, J.L.; Benden, C.; Calabrese, F.; Gottlieb, J.; Hachem, R.R.; Levine, D.; Meloni, F.; Palmer, S.M.; et al. Chronic lung allograft dysfunction: Definition and update of restrictive allograft syndrome—A consensus report from the Pulmonary Council of the ISHLT. J. Heart Lung Transplant. 2019, 38, 483–492. [Google Scholar] [CrossRef] [Green Version]
  6. Lama, V.N.; Murray, S.; Mumford, J.A.; Flaherty, K.R.; Chang, A.; Toews, G.B.; Peters-Golden, M.; Martinez, F.J. Prognostic value of bronchiolitis obliterans syndrome stage 0-p in single-lung transplant recipients. Am. J. Respir. Cri.t Care Med. 2005, 172, 379–383. [Google Scholar] [CrossRef]
  7. Verleden, S.E.; Vasilescu, D.M.; McDonough, J.E.; Ruttens, D.; Vos, R.; Vandermeulen, E.; Bellon, H.; Geenens, R.; Verbeken, E.K.; Verschakelen, J.; et al. Linking clinical phenotypes of chronic lung allograft dysfunction to changes in lung structure. Eur. Respir. J. 2015, 46, 1430–1439. [Google Scholar] [CrossRef] [Green Version]
  8. Verleden, S.E.; Vasilescu, D.M.; Willems, S.; Ruttens, D.; Vos, R.; Vandermeulen, E.; Hostens, J.; McDonough, J.E.; Verbeken, E.K.; Verschakelen, J.; et al. The site and nature of airway obstruction after lung transplantation. Am. J. Respir. Crit. Care Med. 2014, 189, 292–300. [Google Scholar] [CrossRef]
  9. Vandermeulen, E.; Lammertyn, E.; Verleden, S.E.; Ruttens, D.; Bellon, H.; Ricciardi, M.; Somers, J.; Bracke, K.R.; Van Den Eynde, K.; Tousseyn, T.; et al. Immunological diversity in phenotypes of chronic lung allograft dysfunction: A comprehensive immunohistochemical analysis. Transpl. Int. 2017, 30, 134–143. [Google Scholar] [CrossRef] [Green Version]
  10. Galbán, C.J.; Han, M.K.; Boes, J.L.; Chughtai, K.A.; Meyer, C.R.; Johnson, T.D.; Galbán, S.; Rehemtulla, A.; Kazerooni, E.A.; Martinez, F.J.; et al. Computed tomography-based biomarker provides unique signature for diagnosis of COPD phenotypes and disease progression. Nat. Med. 2012, 18, 1711–1715. [Google Scholar] [CrossRef]
  11. Bhatt, S.P.; Soler, X.; Wang, X.; Murray, S.; Anzueto, A.R.; Beaty, T.H.; Boriek, A.M.; Casaburi, R.; Criner, G.J.; Diaz, A.A.; et al. Association between Functional Small Airway Disease and FEV1 Decline in Chronic Obstructive Pulmonary Disease. Am. J. Respir. Crit. Care Med. 2016, 194, 178–184. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Vasilescu, D.M.; Martinez, F.J.; Marchetti, N.; Galban, C.J.; Hatt, C.; Meldrum, C.A.; Dass, C.; Tanabe, N.; Reddy, R.M.; Lagstein, A.; et al. Noninvasive Imaging Biomarker Identifies Small Airway Damage in Severe Chronic Obstructive Pulmonary Disease. Am. J. Respir. Crit. Care Med. 2019, 200, 575–581. [Google Scholar] [CrossRef] [PubMed]
  13. Galbán, C.J.; Boes, J.L.; Bule, M.; Kitko, C.L.; Couriel, D.R.; Johnson, T.D.; Lama, V.; Telenga, E.D.; van den Berge, M.; Rehemtulla, A.; et al. Parametric response mapping as an indicator of bronchiolitis obliterans syndrome after hematopoietic stem cell transplantation. Biol. Blood Marrow Transplant. 2014, 20, 1592–1598. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Belloli, E.A.; Degtiar, I.; Wang, X.; Yanik, G.A.; Stuckey, L.J.; Verleden, S.E.; Kazerooni, E.A.; Ross, B.D.; Murray, S.; Galban, C.J.; et al. Parametric Response Mapping as an Imaging Biomarker in Lung Transplant Recipients. Am. J. Respir. Crit. Care Med. 2017, 195, 942–952. [Google Scholar] [CrossRef] [Green Version]
  15. Sharifi, H.; Lai, Y.K.; Guo, H.; Hoppenfeld, M.; Guenther, Z.D.; Johnston, L.; Brondstetter, T.; Chhatwani, L.; Nicolls, M.R.; Hsu, J.L. Machine Learning Algorithms to Differentiate Among Pulmonary Complications after Hematopoietic Cell Transplant. Chest 2020, 158, 1090–1103. [Google Scholar] [CrossRef]
  16. Verleden, S.E.; Vos, R.; Vandermeulen, E.; Ruttens, D.; Bellon, H.; Heigl, T.; Van Raemdonck, D.E.; Verleden, G.M.; Lama, V.; Ross, B.D.; et al. Parametric Response Mapping of Bronchiolitis Obliterans Syndrome Progression After Lung Transplantation. Am. J. Transplant. 2016, 16, 3262–3269. [Google Scholar] [CrossRef] [Green Version]
  17. Belloli, E.A.; Gu, T.; Wang, Y.; Vummidi, D.; Lyu, D.M.; Combs, M.P.; Chughtai, A.; Murray, S.; Galban, C.J.; Lama, V.N. Radiographic Graft Surveillance in Lung Transplantation: Prognostic Role of Parametric Response Mapping. Am. J. Respir. Crit. Care Med. 2021, 204, 967–976. [Google Scholar] [CrossRef]
  18. Hoff, B.A.; Pompe, E.; Galban, S.; Postma, D.S.; Lammers, J.J.; Ten Hacken, N.H.T.; Koenderman, L.; Johnson, T.D.; Verleden, S.E.; de Jong, P.A.; et al. CT-Based Local Distribution Metric Improves Characterization of COPD. Sci. Rep. 2017, 7, 2999. [Google Scholar] [CrossRef]
  19. Stoel, B.C.; Stolk, J. Optimization and standardization of lung densitometry in the assessment of pulmonary emphysema. Investig. Radiol. 2004, 39, 681–688. [Google Scholar] [CrossRef] [Green Version]
  20. Klein, S.; Staring, M.; Murphy, K.; Viergever, M.A.; Pluim, J.P. elastix: A toolbox for intensity-based medical image registration. IEEE Trans. Med. Imaging 2010, 29, 196–205. [Google Scholar] [CrossRef]
  21. Shamonin, D.P.; Bron, E.E.; Lelieveldt, B.P.; Smits, M.; Klein, S.; Staring, M.; Alzheimer’s Disease Neuroimaging Initiative. Fast parallel image registration on CPU and GPU for diagnostic classification of Alzheimer’s disease. Front. Neuroinform. 2013, 7, 50. [Google Scholar] [CrossRef] [PubMed]
  22. Legland, D.; Kiêu, K.; Devaux, M. Computation of Minkowski Measures on 2D and 3D Binary Images. Image Anal. Stereol. 2011, 26, 83–92. [Google Scholar] [CrossRef] [Green Version]
  23. Verleden, S.E.; Kirby, M.; Everaerts, S.; Vanstapel, A.; McDonough, J.E.; Verbeken, E.K.; Braubach, P.; Boone, M.N.; Aslam, D.; Verschakelen, J.; et al. Small airway loss in the physiologically ageing lung: A cross-sectional study in unused donor lungs. Lancet Respir. Med. 2021, 9, 167–174. [Google Scholar] [CrossRef]
  24. McDonough, J.E.; Ahangari, F.; Li, Q.; Jain, S.; Verleden, S.E.; Herazo-Maya, J.; Vukmirovic, M.; DeIuliis, G.; Tzouvelekis, A.; Tanabe, N.; et al. Transcriptional regulatory model of fibrosis progression in the human lung. JCI Insight 2019, 4, e131597. [Google Scholar] [CrossRef] [PubMed]
  25. Cheng, G.S.; Selwa, K.E.; Hatt, C.; Ram, S.; Fortuna, A.B.; Guerriero, M.; Himelhoch, B.; McAree, D.; Hoffman, T.C.; Brisson, J.; et al. Multicenter evaluation of parametric response mapping as an indicator of bronchiolitis obliterans syndrome after hematopoietic stem cell transplantation. Am. J. Transplant. 2020, 20, 2198–2205. [Google Scholar] [CrossRef]
  26. Verleden, S.E.; Ruttens, D.; Vandermeulen, E.; Bellon, H.; Dubbeldam, A.; De Wever, W.; Dupont, L.J.; Van Raemdonck, D.E.; Vanaudenaerde, B.M.; Verleden, G.M.; et al. Predictors of survival in restrictive chronic lung allograft dysfunction after lung transplantation. J. Heart Lung Transplant. 2016, 35, 1078–1084. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Representative cases diagnosed with BOS and RAS. Presented are the coronal view of CT scans acquired at inspiration and expiration with corresponding topological parametric response mapping (tPRM)-derived volume density maps for normal parenchyma (green), functional small airways disease (fSAD; yellow) and parenchymal disease (PD; magenta).
Figure 1. Representative cases diagnosed with BOS and RAS. Presented are the coronal view of CT scans acquired at inspiration and expiration with corresponding topological parametric response mapping (tPRM)-derived volume density maps for normal parenchyma (green), functional small airways disease (fSAD; yellow) and parenchymal disease (PD; magenta).
Cells 11 00699 g001
Figure 2. Chronic lung allograft dysfunction (CLAD) subtype group comparison of all tPRM metrics. Barplots are presented as the mean and standard error of the mean for topological measures across both BOS (N = 6) and RAS (N = 6), for PRMNorm (green), PRMfSAD (yellow), PRMEmph (red) and PRMPD (magenta). BOS = Bronchiolitis Obliterans Syndrome, RAS = Restrictive Allograft Syndrome, Norm = normal parenchyma, fSAD = functional small airways disease, Emph = emphysema, PD = parenchymal disease. (A) Volume density, describing class magnitude (relative amounts of voxels). (B) Surface area, describing class exposure (exposed faces of voxels). (C) mean breadth, which can be interpreted as average curvature (convex/concave) of class surfaces. (D) Euler–Poincaré characteristic, describing class homology, determined by number and type of holes within class volumes. The legend indicates measures from BOS and RAS as light and dark shades, respectively, for a given PRM color-code.
Figure 2. Chronic lung allograft dysfunction (CLAD) subtype group comparison of all tPRM metrics. Barplots are presented as the mean and standard error of the mean for topological measures across both BOS (N = 6) and RAS (N = 6), for PRMNorm (green), PRMfSAD (yellow), PRMEmph (red) and PRMPD (magenta). BOS = Bronchiolitis Obliterans Syndrome, RAS = Restrictive Allograft Syndrome, Norm = normal parenchyma, fSAD = functional small airways disease, Emph = emphysema, PD = parenchymal disease. (A) Volume density, describing class magnitude (relative amounts of voxels). (B) Surface area, describing class exposure (exposed faces of voxels). (C) mean breadth, which can be interpreted as average curvature (convex/concave) of class surfaces. (D) Euler–Poincaré characteristic, describing class homology, determined by number and type of holes within class volumes. The legend indicates measures from BOS and RAS as light and dark shades, respectively, for a given PRM color-code.
Cells 11 00699 g002
Figure 3. Representative case illustrating spatial heterogeneity in tPRM measures. Presented is a photo of an explanted lung section, with corresponding aligned inspiratory CT scan and tPRM-derived volume density map, surface area and Euler–Poincaré characteristic. Volume density map depicts the local extent of normal parenchyma (green), fSAD (yellow) and PD (magenta). Surface area and Euler–Poincaré characteristic are provided for each PRM classification. The circles on explant (blue) and CT (yellow) provide approximate location and size of an example core sample. This case is a male diagnosed with BOS 447 days post-transplant. This patient had a FEV1pp of 18% and FEV1/FVC of 0.41 when undergoing redo transplantation.
Figure 3. Representative case illustrating spatial heterogeneity in tPRM measures. Presented is a photo of an explanted lung section, with corresponding aligned inspiratory CT scan and tPRM-derived volume density map, surface area and Euler–Poincaré characteristic. Volume density map depicts the local extent of normal parenchyma (green), fSAD (yellow) and PD (magenta). Surface area and Euler–Poincaré characteristic are provided for each PRM classification. The circles on explant (blue) and CT (yellow) provide approximate location and size of an example core sample. This case is a male diagnosed with BOS 447 days post-transplant. This patient had a FEV1pp of 18% and FEV1/FVC of 0.41 when undergoing redo transplantation.
Cells 11 00699 g003
Figure 4. Boxplots for all four tPRM metrics of the dominant readouts from each patient. Presented is the (A) volume density, (B) surface area, (C) mean breadth, and (D) Euler–Poincaré characteristic, for fSAD (yellow) and PD (purple). Approximately four cores were acquired from each patient. Box lines represent the lower quartile, median, and upper quartile, and whiskers represent minimum and maximum values. Cases were separated based on disease subtype.
Figure 4. Boxplots for all four tPRM metrics of the dominant readouts from each patient. Presented is the (A) volume density, (B) surface area, (C) mean breadth, and (D) Euler–Poincaré characteristic, for fSAD (yellow) and PD (purple). Approximately four cores were acquired from each patient. Box lines represent the lower quartile, median, and upper quartile, and whiskers represent minimum and maximum values. Cases were separated based on disease subtype.
Cells 11 00699 g004
Figure 5. Representative microCT and immunohistochemistry for BOS ((A,B), respectively) and RAS ((C,D), respectively) cases. (A) MicroCT image of BOS core where normal airway parenchyma is seen. The arrow indicates an obliterated airway. (B) Neutrophil staining (pink) corresponds to the site of the arrow on (A). (C) MicroCT image of RAS core where some patchy fibroblast foci can be seen. (D) Collagen I stain (pink) is shown at the site of the arrow on (C). Scale bar (red line in (A)) represents 1 mm on both microCT (A,C). (B,C) were counter stained with hematoxylin (blue).
Figure 5. Representative microCT and immunohistochemistry for BOS ((A,B), respectively) and RAS ((C,D), respectively) cases. (A) MicroCT image of BOS core where normal airway parenchyma is seen. The arrow indicates an obliterated airway. (B) Neutrophil staining (pink) corresponds to the site of the arrow on (A). (C) MicroCT image of RAS core where some patchy fibroblast foci can be seen. (D) Collagen I stain (pink) is shown at the site of the arrow on (C). Scale bar (red line in (A)) represents 1 mm on both microCT (A,C). (B,C) were counter stained with hematoxylin (blue).
Cells 11 00699 g005
Table 1. Subject characteristics.
Table 1. Subject characteristics.
MeasurementsBOSRASp-Value
Number of Patients, n66
Donor age, years34 (11)45 (13)0.09
Recipient age at explant, years47 (14)49 (16)
Recipient gender, n (%)
Male2 (33%)3 (50%)
Female4 (67%)3 (50%)
Recipient height, cm170 (5)172 (13)
Recipient weight, kg59 (13)59 (11)
Underlying disease, n (%)
A1ATD0 (0%)1 (17%)
CF3 (49%)1 (17%)
COPD1 (17%)2 (32%)
CHP1 (17%)0 (0%)
IPF1 (17%)0 (0%)
PHT0 (0%)1 (17%)
Systemic Sclerosis + PHT0 (0%)1 (17%)
FEV1 at CT, L0.64 (0.11)0.87 (0.36)
FEV1pp at CT, (%) 20.3 (4.0)30.2 (10.7)0.04
FVC at CT, L1.6 (0.6)1.2 (0.4)
FVCpp at CT, (%) 41.7 (11.7)31.3 (10.0)
FEV1/FVC at CT0.42 (0.08)0.72 (0.20)0.01
CT lung volume, L
Expiration3.7 (1.3)2.3 (0.9)0.07
Inspiration4.1 (1.4)3.1 (1.1)
CT mean lung density, HU
Expiration−806 (55)−678 (112)0.03
Inspiration−826 (34)−754 (73)0.09
Continuous variables presented as mean (standard deviation). Fisher’s Exact test was used to compare sex. Mann–Whitney U Test was used for continuous data. p-values less than 0.1 are provided with significant values bolded (p < 0.05). BOS, bronchiolitis obliterans syndrome; RAS, restrictive allograft syndrome; A1ATD, alpha-1-antitrypsin deficiency; CF, cystic fibrosis; COPD, chronic obstructive pulmonary disease; CHP, chronic hypersensitivity pneumonitis; IPF, idiopathic pulmonary fibrosis; PHT, pulmonary hypertension; FEV1, forced expiratory volume at 1 s; FEV1pp, FEV1 percent predicted; FVC, forced vital capacity; CT, computed tomography; HU, Hounsfield Unit.
Table 2. Correlations of whole lung tPRM metrics of dominant phenotypes to pulmonary function measurements.
Table 2. Correlations of whole lung tPRM metrics of dominant phenotypes to pulmonary function measurements.
MetricFEV1FEV1ppFVCFVCppFEV1/FVC
VfSAD−0.196−0.5040.6660.301−0.825
SfSAD−0.007−0.4900.6800.217−0.615
BfSAD0.0000.198−0.392−0.4620.413
χfSAD−0.1120.374−0.848−0.5940.727
VPD−0.0840.413−0.543−0.2310.420
SPD−0.1400.307−0.743−0.4340.615
BPD−0.098−0.367−0.168−0.3850.196
χPD0.007−0.4900.5220.098−0.441
Notes: Presented are the correlation coefficients (R) obtained using Spearman rho on pooled data (N = 6 BOS and 6 RAS). Conditional formatting from red (min value) to blue (max value) through white (median value) has been applied over the entire table. Statistical significance was assessed at p < 0.05 and is indicated with yellow font.
Table 3. Structural and Biological Measures from Core Samples.
Table 3. Structural and Biological Measures from Core Samples.
MeasuresBOSRASp-Value
Number of Cores, n2322
Mean number of cores per patient, n3.83.7
Core Volume, mL1.9 × 100 (2.9 × 10−1)2.0 × 100 (4.0 × 10−1)
StructuralTissue Percentage of Core, %2.9 × 101 (5.3 × 100)5.2 × 101 (1.9 × 101)<0.001
Surface density, µm−11.3 × 10−2 (1.8 × 10−3)8.9 × 10−3 (3.1 × 10−3)<0.001
Number of terminal bronchioles per volume, (n/mL)4.2 × 100 (1.9 × 100)4.9 × 100 (2.2 × 100)
BiologicalB cells (CD20)3.1 × 10−3 (8.1 × 10−3)3.2 × 10−3 (9.0 × 10−3)
Neutrophils (MPO)4.3 × 10−2 (3.1 × 10−2)3.0 × 10−2 (3.1 × 10−2)0.06
Mast cells (CD117)6.1 × 10−5 (1.1 × 10−4)1.3 × 10−4 (2.0 × 10−4)
Macrophages (CD68)3.4 × 10−2 (2.7 × 10−2)3.9 × 10−2 (2.9 × 10−2)
CD4+ T cells5.5 × 10−4 (1.7 × 10−3)1.3 × 10−3 (2.1 × 10−3)0.02
CD8+ T cells4.1 × 10−3 (7.3 × 10−3)3.5 × 10−3 (4.6 × 10−3)
Eosinophils (ECP)6.3 × 10−3 (1.2 × 10−2)4.3 × 10−3 (7.7 × 10−3)
Collagen I1.6 × 10−3 (1.3 × 10−3)2.2 × 10−3 (2.1 × 10−3)
Collagen III1.2 × 10−2 (1.8 × 10−2)1.9 × 10−2 (4.1 × 10−2)
Continuous variables presented as mean (standard deviation). Group comparisons were performed using Mann–Whitney U test. Statistical significance was assessed at p < 0.05. p-values less than 0.1 are provided with significant values bolded. Tissue percentage of Core (%) was determined by normalizing the tissue volume to the core volume times 100. Biological measurements were determined by normalizing the number of positively stained cells to all identified cells.
Table 4. Correlations of tPRM metrics of dominant subtypes to structural and biological measurements.
Table 4. Correlations of tPRM metrics of dominant subtypes to structural and biological measurements.
MetricT%SDTB/mLCD20MPOCD117CD68CD4CD8ECPCOL ICOL III
Pooled (N = 45)VfSAD−0.5990.448−0.282−0.1150.342−0.310−0.037−0.171−0.225−0.036−0.1980.207
SfSAD−0.4450.309−0.270−0.0620.450−0.3130.059−0.060−0.1350.031−0.2870.167
BfSAD0.262−0.315−0.0900.1210.1230.0130.0370.1810.148−0.004−0.104−0.101
χfSAD0.479−0.4530.0990.118−0.2650.184−0.0490.1710.119−0.0480.146−0.161
VPD0.457−0.2360.380−0.062−0.4700.2620.0740.0690.0880.0500.292−0.076
SPD0.425−0.2490.287−0.020−0.3300.2330.0220.0880.2170.0270.151−0.081
BPD−0.063−0.040−0.0950.0880.274−0.135−0.1040.1590.190−0.093−0.260−0.018
χPD−0.209−0.023−0.2310.1110.374−0.2410.0600.1540.035−0.043−0.3520.020
BOS (N = 23)VfSAD−0.304−0.070−0.132−0.045−0.106−0.1810.1190.122−0.225−0.008−0.0700.307
SfSAD0.030−0.121−0.002−0.0840.289−0.2060.2420.157−0.0540.063−0.2610.157
BfSAD0.292−0.0060.1250.0110.1250.075−0.145−0.1190.106−0.1000.109−0.299
χfSAD0.036−0.0150.053−0.052−0.0600.091−0.242−0.248−0.103−0.2070.278−0.310
VPD0.1480.2830.363−0.0130.0180.150−0.112−0.1030.1120.0130.124−0.236
SPD0.2270.2090.3200.0310.0400.108−0.093−0.0540.137−0.0010.126−0.252
BPD0.1300.0500.169−0.0530.104−0.018−0.113−0.0530.147−0.1060.063−0.268
χPD−0.045−0.128−0.125−0.1150.1160.0100.0980.0170.221−0.092−0.296−0.039
RAS (N = 22)VfSAD−0.119−0.078−0.3980.1820.399−0.3120.0200.169−0.1490.119−0.2750.135
SfSAD−0.129−0.115−0.4600.1620.460−0.2970.0540.180−0.0850.098−0.3270.133
BfSAD0.053−0.394−0.4590.1150.319−0.2470.1060.3140.083−0.019−0.3800.212
χfSAD0.224−0.4030.029−0.004−0.141−0.054−0.0180.177−0.020−0.189−0.0670.072
VPD0.217−0.1560.328−0.225−0.6840.2700.130−0.112−0.057−0.0120.4780.124
SPD0.188−0.2660.136−0.117−0.3460.2520.0220.0360.337−0.1110.0920.097
BPD−0.097−0.224−0.2730.2590.420−0.206−0.0580.4040.283−0.063−0.4910.119
χPD−0.238−0.049−0.2470.3110.390−0.4190.0310.305−0.1580.060−0.2890.106
Presented are the correlation coefficients (R) obtained using Spearman rho on pooled data from cores (N = 23 BOS and 22 RAS). T%, Tissue percentage of core; SD, Surface density; TB/mL, Number of terminal bronchioles per volume; CD20, B cells; MPO, Neutrophils; CD117, Mast cells; CD68, Macrophages; CD4, CD4+ T cells; CD8, CD8+ T cells; ECP, Eosinophils; COL I, Collagen I; COL III, Collagen III. Conditional formatting from red (min value) to blue (max value) through white (median value) has been applied over the entire table. Statistical significance was assessed at p < 0.05 and is indicated with yellow font.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ram, S.; Verleden, S.E.; Bell, A.J.; Hoff, B.A.; Labaki, W.W.; Murray, S.; Vanaudenaerde, B.M.; Vos, R.; Verleden, G.M.; Kazerooni, E.A.; et al. Quantitative CT Correlates with Local Inflammation in Lung of Patients with Subtypes of Chronic Lung Allograft Dysfunction. Cells 2022, 11, 699. https://doi.org/10.3390/cells11040699

AMA Style

Ram S, Verleden SE, Bell AJ, Hoff BA, Labaki WW, Murray S, Vanaudenaerde BM, Vos R, Verleden GM, Kazerooni EA, et al. Quantitative CT Correlates with Local Inflammation in Lung of Patients with Subtypes of Chronic Lung Allograft Dysfunction. Cells. 2022; 11(4):699. https://doi.org/10.3390/cells11040699

Chicago/Turabian Style

Ram, Sundaresh, Stijn E. Verleden, Alexander J. Bell, Benjamin A. Hoff, Wassim W. Labaki, Susan Murray, Bart M. Vanaudenaerde, Robin Vos, Geert M. Verleden, Ella A. Kazerooni, and et al. 2022. "Quantitative CT Correlates with Local Inflammation in Lung of Patients with Subtypes of Chronic Lung Allograft Dysfunction" Cells 11, no. 4: 699. https://doi.org/10.3390/cells11040699

APA Style

Ram, S., Verleden, S. E., Bell, A. J., Hoff, B. A., Labaki, W. W., Murray, S., Vanaudenaerde, B. M., Vos, R., Verleden, G. M., Kazerooni, E. A., Galbán, S., Hatt, C. R., Han, M. K., Lama, V. N., & Galbán, C. J. (2022). Quantitative CT Correlates with Local Inflammation in Lung of Patients with Subtypes of Chronic Lung Allograft Dysfunction. Cells, 11(4), 699. https://doi.org/10.3390/cells11040699

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop