Exploratory Study on COPD Phenotypes and their Progression: Integrating SPECT and qCT Imaging Analysis

Frank Li, Xuan Zhang, Alejandro P. Comellas, Eric A. Hoffman, Michael M. Graham, Ching-Long Lin
doi: https://doi.org/10.1101/2024.04.10.24305577

ABSTRACT

Background The objective of this study is to understand chronic obstructive pulmonary disease (COPD) phenotypes and their progressions by quantifying heterogeneities of lung ventilation from the single photon emission computed tomography (SPECT) images and establishing associations with the quantitative computed tomography (qCT) imaging-based clusters and variables.

Methods Eight COPD patients completed a longitudinal study of three visits with intervals of about a year. CT scans of these subjects at residual volume, functional residual capacity, and total lung capacity were taken for all visits. The functional and structural qCT-based variables were derived, and the subjects were classified into the qCT-based clusters. In addition, the SPECT variables were derived to quantify the heterogeneity of lung ventilation. The correlations between the key qCT-based variables and SPECT-based variables were examined.

Results The SPECT-based coefficient of variation (CVTotal), a measure of ventilation heterogeneity, showed strong correlations (|r| ≥ 0.7) with the qCT-based functional small airway disease percentage (fSAD%Total) and emphysematous tissue percentage (Emph%Total) in the total lung on cross-sectional data. As for the two-year changes, the SPECT-based maximum tracer concentration (TCmax), a measure of hot spots, exhibited strong negative correlations with fSAD%Total, Emph%Total, average airway diameter in the left upper lobe, and airflow distribution in the middle and lower lobes.

Conclusion Small airway disease is highly associated with the heterogeneity of ventilation in COPD lungs. TCmax is a more sensitive functional biomarker for COPD progression than CVTotal. Besides fSAD%Total and Emph%Total, segmental airways narrowing and imbalanced ventilation between upper and lower lobes may contribute to the development of hot spots over time.

INTRODUCTION

Chronic obstructive pulmonary disease (COPD) is a prevalent respiratory disease that imposes a significant burden on healthcare systems worldwide (1). The Global Initiative for Chronic Obstructive Lung Disease (GOLD) report (2) has been continuously updated to better characterize COPD heterogeneities and improve treatment outcomes (3,4). In addition, imaging modalities, such as computed tomography (CT) and single-photon emission computed tomography (SPECT), have been utilized to visualize and quantify regional functional and structural changes associated with various pathological processes in the lung diseases such as emphysema (5).

Innovative biomarkers based on quantitative CT (qCT) imaging have been developed to investigate the underlying causes of COPD. These biomarkers assess a wide range of risk factors and associated defects, such as airway-branch variation (6), dysanapsis (7), and pulmonary vascular dysfunction (8,9). The introduction of the qCT-based parametric response map (PRM) (10) has also significantly advanced research on functional small airway disease (fSAD) and emphysema (Emph). Both fSAD and Emph are critical phenotypes in COPD, with the former considered a precursor to the latter (11). It is suggested that targeting small airways with appropriate treatments may potentially control the progression of both airway and parenchymal diseases in COPD (12). Novel multiscale qCT biomarkers have been developed to capture a wide range of COPD phenotypes at different disease stages. These imaging-based biomarkers allow for identification of clusters (also known as subtypes, subgroups, or subpopulations) of COPD patients that show strong associations with clinical characteristics. For instance, Haghighi et al. applied unsupervised clustering techniques to cross-sectional qCT biomarkers and identified four clinically-relevant imaging-based clusters among former and current smokers, respectively, to eliminate the effect of inflammation on the emphysema index (13,14). These clusters were sorted in order of severity, with cluster 1 (C1) being the mildest and cluster 4 (C4) being the most severe in terms of airflow obstruction. Classification of COPD patients into imaging-based clusters can aid in comprehending the heterogeneous nature of the disease more quantitatively and locally. In addition, Zou et al. utilized the multiscale qCT phenotyping and clustering mentioned above to analyze longitudinal data from two visits to identify COPD progression clusters (15). Cross-sectional clustering is based on static severity stages or phenotypes, while longitudinal clustering aims to capture disease activity mechanisms or endotypes.

Assessment of global and regional lung ventilation can be achieved through SPECT ventilation imaging, which measures the concentration of a radioactive tracer serving as a biomarker for ventilation if the tracer aerosol size is small enough to reach alveoli. De Backer et al. demonstrated a strong correlation between CT-based lobar air volume change and lobar SPECT tracer aerosol concentration with aerosol diameter less than 2µm (16). They further demonstrated that the hot spots observed on SPECT images corresponded to airway constrictions on CT images as well as the hot spots determined by computational fluid dynamics simulation. Studies have also demonstrated a correlation between heterogeneities observed on SPECT images and impaired lung functions. Xu et al. found that the coefficient of variation (CV) observed in SPECT images, which is a measure of ventilation heterogeneity, not only differentiated patients with emphysema from non-emphysematous smokers and non-smokers, but also exhibited a correlation with pulmonary functions. (17). Thus, the use of combined SPECT and CT imaging enables the evaluation of the functional and structural relationships within different regions of the lungs (5). Recently approved by the FDA, hyperpolarized xenon-129 magnetic resonance imaging (XeMRI) can serve as an alternative approach to evaluate lung ventilation distribution without exposing patients to ionizing radiation (18–21).

In this study, we aimed to evaluate COPD progression by quantifying structural and functional changes in CT and SPECT images of COPD subjects acquired at three visits, with approximately one year between each visit. Our objective was to establish the links between qCT variables and SPECT biomarkers within the previously established COPD qCT-based clusters (13,14). Given the sample size in this study, the associations observed with the imaging clusters may assist with our understanding of COPD progression. To the best of our knowledge, this study was the first attempt to conduct a longitudinal CT/SPECT analysis, shedding light on COPD progression.

METHOD

Fourteen COPD patients were recruited to participate in a longitudinal study of three visits at 0, 12.84 ± 1.68, and 26.20 ± 3.06 months, denoted by V0, V1, and V2, respectively. Five subjects dropped out of the study due to other health issues or personal reasons, and one subject was excluded due to an abnormal airway structure having an accessary bronchus connected to the right main bronchus. Thus, eight subjects who completed the three visits were analyzed in this study. Pulmonary function tests (PFT), CT scans, and SPECT scans were acquired at all visits. This study was approved by the institutional review board of the University of Iowa and the informed consent was obtained from all the patients before the study.

CT Images and qCT Variables

Three static three-dimensional (3D) CT scans were acquired using the Siemens SOMATOM Force Scanner at the lung volumes of total lung capacity (TLC), functional residual capacity (FRC) and residual volume (RV), using the lung volume control system proposed in the previous study (22). The TLC scan protocol used dose modulation with 36 reference mAs, 120 kV, pitch of 1.0, and rotation time of 0.25 s. The FRC and RV scans used dose modulation with 15 reference mAs, 120 kV, pitch of 1.0, and rotation time of 0.25 s The CT images have a spatial resolution of 0.5 mm in the z axis, and 0.5-0.7 mm in the x and y axes. The subjects were positioned on the CT table in the supine position and were instructed to breathe normally through the mouthpiece for a few breathing cycles. While the CT scans were performed, they were instructed to take in a deepest breath and hold it at the desired lung volume of TLC, FRC, or RV. The subject breathed normally for a few breathing cycles before the next CT scan was performed.

The CT images were processed by VIDA Vision (VIDA Diagnostics Inc., Coralville, Iowa) to segment the lung, lobes, and central airways. A mass preserving image registration technique was performed to match FRC (or RV) with TLC images (23). TLC and RV images were registered to derive functional variables, while TLC and FRC images were registered to match SPECT images. A total of 69 structural and functional multiscale qCT variables presented in previous studies were derived (13–15) as shown in Figure 1. The structural and functional variables describe the regional alterations of lung structures and capture the regional alterations of lung functions, respectively. The qCT variables most relevant to this study include normalized airway wall thickness (WT*), normalized airway hydraulic diameter (Dh*), fractional air volume change (Embedded Image), determinant of Jacobian matrix (J), anisotropic deformation index (ADI), fraction-based small airways disease (fSAD%), fraction-based emphysema (Emph%), and tissue fraction at TLC (βtissue). Please refer to Supplementary Material for detailed explanation of these variables and a list of abbreviations of qCT variables. Subsequently, the subjects were classified into pre-established imaging-based clusters (13,14) using these qCT variables, aiming to enhance our understanding of COPD progression within a broader population context. The assignment of cluster membership was determined based on the distances to the geometric centers of the clusters within the qCT feature space.

Figure 1.
Figure 1.

The multiscale structural and functional qCT variables.

The region, where the qCT variable was measured, was specified as the subscript of the variable, namely {Variable}{Region}. The structural variables were measured from the lung shape, central airways, and segmental airways. The lung shape was measured at TLC as the ratio of apical-basal distance and ventral-dorsal distance. The central airways include trachea, left main bronchus (LMB), right main bronchus (RMB), right intermediate bronchus (BronInt), and trifurcation of left lower lobe (TriLLB). The segmental airways were grouped by lobes, such as sub-grouped segmental airways of left upper lobe (sLUL), sub-grouped segmental airways of left lower lobe (sLLL), sub-grouped segmental airways of right upper lobe (sRUL), sub-grouped segmental airways of right middle lobe (sRML), and sub-grouped segmental airways of right lower lobe (sRLL). The functional variables were measured at lobar level, including left upper lobe (LUL), left lower lobe (LLL), right upper lobe (RUL), right middle lobe (RML), and right lower lobe (RLL), as well as total lung level (Total).

SPECT Images and Variables

The subject was positioned in the supine position on the SPECT scanner to duplicate the same position used in the above CT scans. The patient inhaled through a mouthpiece, with the nose occluded, 99mTc sulfur colloid as a radioactive tracer aerosol that was generated by a specialized nebulizer. Then a 3D SPECT scan of the whole lung was performed. The particle size of the sulfur colloid is below 1 µm (24), enabling distribution of the tracer aerosol deep into the lung (25). Since 99mTc sulfur colloid is non-absorbable, it is not cleared from the lung rapidly. The 3D SPECT scan took approximately 20 minutes to complete and was acquired during continuous tidal breathing. The SPECT images have a spatial resolution of 3.895 mm in the x, y, and z axes. It was estimated that the subject retained a dose of 0.17 mSv, which was lower than the standard of the U.S. Nuclear Regulatory Commission (50 mSv/year). It is noted that the original SPECT imaging protocol was to use 99mTc-HMPAO labeled neutrophils to evaluate the distribution of inflammation, and was adopted for the first two subjects at V0. However, the resulting SPECT images did not provide the expected results. We then modified the protocol. Thus, the first two subjects did not have SPECT ventilation scans at V0.

SPECT image is a dynamic scan which is a time-averaged image acquired during tidal breathing. On the other hand, the CT images are static scans acquired at TLC, FRC, and RV. It is preferable to match the SPECT images with the FRC CT images, which were taken at a lung volume closer to tidal breathing, than with the RV or TLC CT images. The lung deformation from FRC to the peak of tidal volume was assumed relatively small and quasi-linear. Thus, the SPECT images were first matched with the FRC CT images by affine transforms, using mutual information as the cost function. Displacement fields of transformation from FRC to TLC images were obtained by image registration. Subsequently, the transformed SPECT images were deformed by applying the displacement fields to match the TLC CT images, which provide detailed anatomical information of the lungs.

After matching the SPECT images with TLC CT images, the tracer concentrations (TC) within the lung and each lobe were quantified at the domain of TLC CT images. The ratio of TC in a lobe over the total lung (TC%), TC coefficient of variation of the total lung (CVTotal), and standardized maximum TC (TCMax) in the total lung were used to describe the distribution, heterogeneity, and hot spot of TC. Specifically, TC% was used to measure the lobar ventilation distribution. CVTotal was a standardized measure of the heterogeneity of TC distribution in the total lung, calculated by normalizing the standard deviation with the average of TC (17). TCMax was used to quantify the local hot spots magnitude as a measure of airway constriction. The normalized CVTotal and TCMax aimed to factor out the effect of different amount of TC inhaled. Moreover, the correlation between TC% and Embedded Image was used to assess the efficacy of the image registration of SPECT and CT images. The Pearson’s correlation coefficients between CVTotal and PFT values were calculated and validated against previous studies (17,26).

Analysis of SPECT and qCT Variables

The associations between SPECT and qCT variables were evaluated to assist the interpretation of the SPECT variables. Exploratory factor analysis (EFA) was applied to reduce the large number of qCT variables derived from the 799 current and former smokers (13,14) to a lower number of factors that preserved the variability of qCT-captured features. The number of factors was determined by the parallel analysis (27). The factors were extracted by principal component analysis (PCA) with Varimax rotation (28). Varimax rotation is an orthogonal rotation, so the factors are independent.

We then selected key qCT variables contributing significantly to each factor as the surrogates of the factors for interpretability of disease phenotypes. The qCT variables were regarded as key variables if their loading values on their corresponding factors were greater than 0.7. We then examined the correlations between key qCT variables and SPECT variables, based on the cross-sectional data at respective V0, V1, and V2 as well as the longitudinal data with changes between any two visits of V1-V0, V2-V1, and V2-V0. Since this study was exploratory and the sample size was small, we only focused on the large effect size, specifically the Pearson’s correlations with coefficient greater than 0.7 (29). A large effect size indicates that even with a small sample size, the observed effect is likely to be substantial and meaningful.

RESULT

Five former smokers (Subject 1-5) and three current smokers (Subject 6-8) were analyzed (age: 63.1 ± 11.5 years; range: 51-84 years). The cluster characteristics of former smokers and current smokers were summarized in Table 1. The cluster membership, clinical data, and PFTs of these subjects were presented in Table 2 (see Table S1 for more data). The subjects had an average of Tiffeneau ratio (FEV1/FVC%) of 52.9 ± 16.6 %, and an average of predicted forced expiratory volume ratio (FEV1% predicted) of 67.9 ± 16.3% at V0. The mean changes between V1 and V0 in terms of FEV1/FVC% and FEV1% predicted were 0.0 ± 1.9% and −2.5 ± 5.4%, respectively. Moreover, the mean changes between V2 and V1 in terms of FEV1/FVC% and FEV1% predicted were 0.8 ± 1.9% and 3.5 ± 4.6%, respectively.

View this table:
Table 1. Characteristics of clusters of the current and former smokers, defined by Haghighi et al. (13,14).
View this table:
Table 2. Clinical data, PFT at V0, and changes in terms of PFT between visits. Subject 1-5, former smokers; Subject 6-8, current smokers.

SPECT Quantification

TC% and Embedded Image of each lobe at each visit were listed in Table S2 of the supplementary material. The correlation coefficient between TC% and Embedded Image was 0.73, indicating a strong correlation between the two variables (Figure 2 (a)). In addition, the correlations between CVTotal and PFT values showed a strong negative correlation with FEV1% predicted at −0.74 and an even stronger negative correlation with FEV1/FVC% at −0.80 (Figure 2 (b)).

Figure 2.
Figure 2.

(a) A scatter plot showing the strong correlation (r=0.73) between TC% and ΔVairF of the lobes. (b) A scatter plot showing the strong correlations between CVTotal and PFT results (FEV1% predicted: r=-0.74, FEV1/FVC%: r=-0.80). The data points were from both three visits.

The SPECT images for all the subjects were displayed in Figure 3. CVTotal and TCMax at each visit together with their changes between visits were shown in Table 3. At all the visits, Subject 2 and 5 were characterized by greater CVTotal and extremely high TCMax, while Subject 1, 6, 7, and 8 exhibited more homogeneous TC distributions. Subject 3 and 4 had large TCMax at V1 and V2, respectively. Hot spots revealed by concentrated red areas were notable in the SPECT images of Subject 2 and 5 for all the visits. A distinguished hot spot was also observed in the SPECT image of Subject 3 at V2. The changes of CVTotal and TCMax between visits for all the subjects were small, being around 10%.

Figure 3.
Figure 3.

Normalized SPECT images transformed to TLC domain at the three visits. Subject 1 and 2 did not have ventilation scans at V0.

View this table:
Table 3. CVTotal and TCMax at each visit, and their changes between visits.

Subject Classification based on qCT Variables

At V0, Subject 1-5 were classified into respective cluster 3, 4, 1, 3, and 3 (C3, C4, C1, C3, and C3), of former-smoker clusters (14), respectively. The memberships of these subjects did not change from V0 to V1. However, from V1 to V2, Subject 1 changed from C3 to C1, Subject 5 moved from C3 to C4, and Subject 2-4 remained the same clusters. On the other hand, Subject 6, 7, and 8 were classified into C1, C3, and C1 of current-smoker clusters (13), respectively, at both V0 and V1. However, from V1 to V2, Subject 7 changed from C3 to C1 and Subject 6 and 8 remained the same clusters. In summary, among these 24 imaging-based classifications (for 8 subjects at 3 visits), only three changed at V2, indicating the stability of cluster membership over two years. These classifications were placed on the space formed by the first and the second principal components (PC0 and PC1) to visualize the progressive loci of these subjects in the cluster landscape depicted by the cohorts of former and current smokers studied previously (13,14) (Figure 4). The imaging metrics that contributed the most to PC0 were fSAD%Total and Emph%Total, while the imaging metrics that contributed the most to PC1 were βtissue,Total, JTotal, and ADITotal (Figure S1). Subject 1 (former smoker) and Subject 7 (current smoker) displayed the similar loci, moving downward along PC1 and changing from C3 to C1 at V2, while Subject 5 moving rightward along PC0 and changing from C3 to C4 at V2. Subject 2, 3, 4, 6, and 8 remained within their respective cluster domains during the two-year duration.

Figure 4.
Figure 4.

(a) Former smokers in this study placed on the two-dimensional space formed by the first and the second principal components derived from previous study (14). (b) Current smokers in this study placed on the two-dimensional space formed by the first and the second principal components derived from previous study (13). C1, C2, C3, and C4 in the figure legend denote all of the subjects analyzed in previous studies (13,14). Those enclosed by a circle denote the change of cluster membership.

Factors Extracted from Former and Current Smokers

Seven factors were extracted from the qCT variables of the former and current smokers (13,14). The factors were able to explain 78.4% of the variance of the original features. The important qCT variables which had moderate-to-high and high contribution (factor loading > 0.6) to the factors were listed in Table 4. ADITotal and JTotal were selected to represent Factor 0 (F0). βtissue,Total and fSAD%Total were selected to represent Factor 1 (F1). Dh*sLUL was selected to represent Factor 2 (F2). Emph%Total was selected to represent Factor 3 (F3). Embedded Image was selected to represent Factor 4 (F4). WT*sLLL was selected to represent Factor 5 (F5). Dh*LMB was selected to represent Factor 6 (F6).

View this table:
Table 4. The important qCT variables which have large contribution to the factors. The variables shown in this table had the magnitude of factor loadings greater than 0.6. A “(-)” indicates a negative factor loading.

Association between SPECT and qCT Variables

Association tests were conducted using both cross-sectional and longitudinal data. The cross-sectional tests captured inter-subject variation, while the longitudinal tests captured intra-subject progression. On cross-sectional data, both SPECT variables of CVTotal and TCMax had the strongest correlations with fSAD%Total and Emph%Total among all qCT variables. The correlations of CVTotal with fSAD%Total and Emph%Total were 0.90 and 0.71, respectively. Similarly, the correlations of TCMax with fSAD%Total and Emph%Total were 0.86 and 0.77, respectively. The strength of correlations between the one-year changes in SPECT variables and those in qCT variables was found to be weak. As for the two-year changes, TCMax exhibited strong correlations with several qCT variables, including fSAD%Total (r=-0.70), Dh,sLUL* (r=-0.74), Emph%Total (r=- 0.75), and Embedded Image (r=0.75) (Figure 5). This finding suggests that TCmax is a more sensitive functional biomarker for COPD progression than CVTotal.

Figure 5.
Figure 5.

The correlations between the SPECT variables and the qCT variables for cross-sectional and longitudinal data. Magnitudes of correlations greater than 0.7 were bounded by blue boxes.

On longitudinal data, the selected qCT variables at each visit and their changes between visits were shown in Table 5. For all the visits, Subject 2 and 5 exhibited higher fSAD%Total and Emph%Total compared to the other subjects. Moreover, Subject 1, 3, and 4 had moderate fSAD%Total and Emph%Total, while Subject 6-8 had low fSAD%Total and Emph%Total. However, the changes of the selected key qCT variables between visits for all subjects were small. The fSAD maps at V1, which marked the fSAD-pixels on CT coronal slices of the subjects, were compared with the SPECT images in Figure 6. Visual inspection supported a correlation between qCT-based fSAD maps and SPECT-based ventilation patterns by groups of Subject (1,3,4), (2,5), and (6,7,8) with moderate, high, and low fSAD%Total, respectively.

Figure 6.
Figure 6.

The fSAD-voxel maps for each subject (left column) along with their SPECT images (right column). These images were plotted in coronal planes of the CT images

View this table:
Table 5. The selected qCT variables at each visit, and their changes between visits.

DISCUSSION

In this work, five former smokers (Subject 1-5) and three current smokers (Subject 6-8) with COPD were studied using combined CT and SPECT imaging biomarkers. At V0, Subject 2 and 5 had moderate-to-severe (GOLD 2/3) and severe (GOLD 3) airflow obstruction, respectively. Subject 1, 3, 4, and 6 had either mild (GOLD 1) or moderate (GOLD 2) airflow obstruction. Subject 7 and 8 were regarded as subjects at risk of COPD (GOLD 0) according to their PFT results. The lung functions of all the subjects were well maintained over two years (Table 2).

The robust correlation observed between TC% and Embedded Image of the lobes (Figure 2: r=0.73) suggests that qCT-based Embedded Image can be utilized as an alternative measure for lobar ventilation. The observed correlation also serves as evidence supporting the effectiveness of the registration process utilized to align SPECT and CT images. Furthermore, the notable negative correlations were found between TC heterogeneity (CVTotal) and PFT results (Figure 2: r=-0.74 with FEV1% predicted; r=- 0.80 with FEV1/FVC%), agreeing with previous studies (r=-0.84 to -0.89 with FEV1% predicted) (17,18).

Subject Classification based on qCT Variables

Subject 2 and 5 were patients with severe COPD. They had the most heterogeneous TC distribution with higher CVTotal and TCMax compared to other subjects (Figure 3). The former-smoker clusters of C3 and C4 were characterized by elevated small airway disease and emphysema (Table 1). Subject 2 was classified into C4 for all the visits. Subject 5 was classified into C3 at V0, but progressed to C4 at V2, remaining near the boundary of C3 and C4 in the 2D latent space for all the visits (Figure 4 (a)). Thus, as expected, Subject 2 and 5 had greater fSAD%Total and Emph%Total than others.

Subject 1, 3, and 4 were COPD patients at moderate stages. Subject 1 and 4 were classified as C3 at V0 and V1. Although Subject 1 changed to C1 at V2, this subject remained near the boundary of C1 and C3 (Figure 4 (a)). Therefore, compared to C1 subjects, greater fSAD%Total was observed, yet elevated Emph%Total was not found for Subject 1 and 4 at all visits. This may explain why they were close to the boundaries of C1 and C2 and shared the characteristics of the two clusters, having marginal emphysema with relatively preserved lung function (14). It is noteworthy that the lung functions may not be sensitive enough to detect the elevated fSAD% (30). Due to the greater JTotal and less βtissue, Total found in Subject 1 at V2, the cluster membership shifted from C3 to C1. Subject 1 and 4 showed moderate dispersion of the TC distribution, exhibiting more heterogeneous TC distribution than Subject 6-8 (Figure 3). Subject 3 was classified as C1 for all the visits due to large lung deformation from RV to TLC lung (as measured by JTotal) regardless of the fact that Subject 3 had increased fSAD%Total and Emph%Total. The large deformation might be attributable to the hyperinflated lung caused by emphysematous regions with less elastic air sacs (31). Since this subject was close to the boundary of C1 and C3 (Figure 4 (a)) and had elevated fSAD%Total and Emph%Total, Subject 3 may be regarded as a C3 subject. The TC distributions for V0 and V1 were similar to those of Subject 1 and 4 at all visits. Furthermore, a clear hot spot was observed on the SPECT image at V2 (Figure 3).

Subject 6-8 were patients with mild COPD. They can be treated as C1 of current smokers at all visits given that they had marginal emphysema and fSAD as well as preserved lung function (13), although Subject 7 was located at the boundary of C1 and C3 at V0 and V1. They showed homogeneous TC distributions and no clear hot spot on the SPECT images for all visits (Figure 3).

Our findings indicate that fSAD was the primary factor determining both disease stage and the heterogeneity of TC distribution, and it was the most significant imaging variable distinguishing the clusters. Furthermore, the clusters appeared to provide a more sensitive representation of lung ventilation patterns compared to GOLD stages.

Association between SPECT and qCT Variables

CVTotal and TCMax demonstrated strong correlations with fSAD%Total (CVTotal: r=0.90; TCMax: r=0.86) and Emph%Total (CVTotal: r=0.71; TCMax: r=0.77) on cross-sectional data, meaning that fSAD and emphysema were associated with heterogeneity of lung ventilation. Compared to emphysema, fSAD was found to have a stronger association with the heterogeneity of lung ventilation. To further interpret and complement the aforementioned findings, a cross-lagged panel analysis study (see Supplementary Material) was conducted to examine the causal relationship between fSAD and heterogeneity of lung ventilation. The analysis indicated that fSAD is the cause of heterogeneity of lung ventilation. In COPD, narrowing of small airways leads to higher airflow resistance and greater particle deposition in the lungs (32–34). As a result, subjects with greater fSAD%Total were observed to have localized areas of increased TC activity, or ‘hot spots’. Regarding the longitudinal data collected over two years, it was observed that there were strong correlations between TCMax and Dh,sLUL* (r=-0.74) as well as Embedded Image (r=0.75) (Figure 5). These findings suggest that reduced diameters of segmental airways and imbalanced ventilation between the upper and lower lobes may also play a role in the development of TC hot spots. Upon visualizing the segmental airways in the LUL of Subject 3 and 5, there existed constricted branches at V2 and V0, respectively (Figure 7). This finding reveals the negative correlations between the changes of Dh,sLUL* and the changes of SPECT ventilation images. This observation is in line with a previous study that identified Dh,sLUL* as a significant discriminant qCT variable for clusters of former smokers (14). Since changes of TCMax over two years were negatively correlated with those of fSAD% (r=-0.70) and Emph% (r=-0.75), changes in Dh,sLUL* and Embedded Image could moderate the correlation between changes in TCMax and changes in fSAD%Total and Emph%Total when the disease stages remained relatively stable over time. This suggests that the formation of hot spots in individuals with COPD could be influenced by alterations in the segmental airways located in the left upper lobe, as well as changes in the distribution of airflow between the upper and lower lobes, in addition to factors such as fSAD and emphysema. Because of the small sample size and the minimal changes observed in SPECT and qCT variables between visits, further investigation with a larger sample size is necessary to validate the correlations found in the longitudinal data.

Figure 7.
Figure 7.

A constricted segmental airway in LUL (red arrow) of Subject 3 at V2, which may contribute to the TC hot spot found in LUL. In addition, a constricted segmental airway in LUL (red arrow) of Subject 5 at V0 resolved at V2, which may increase the homogeneity of ventilation.

This study was subject to several limitations. Firstly, the spatial resolution of SPECT images, which was 3.985 mm, may not be as precise as that of CT scans, which had a resolution of 0.5-0.7 mm. Consequently, there may be some inaccuracies in the quantification of SPECT images in certain regions (such as instances where the tracer may appear outside of the lung on CT images). Furthermore, the radiation exposure associated with SPECT is not insignificant, and its widespread availability is limited. An alternative avenue for this study involves lung ventilation derived from the XeMRI technique. This approach has been demonstrated to exhibit a high correlation with ventilation SPECT (18). Secondly, the sample size was relatively small, and thus increasing the size of the sample may be necessary to obtain more robust and reliable results. Furthermore, the fluctuations of JTotal between visits observed in Subject 1 and 7 may cause the unstable cluster memberships when the subjects were located near the boundaries of two clusters, even though there was no clear disease progression or recovery. These fluctuations of JTotal might be related to the variations of lung volume control during the CT scans. Lastly, it should be noted that the subjects in this study effectively managed their disease for a period of two years in terms of severity, and no significant changes in lung function or imaging were observed during this time. Thus, in order to obtain a clearer understanding of the progression of COPD, it may be necessary to extend the interval between visits beyond the two-year period evaluated in this study.

CONCLUSION

To summarize, our analysis revealed three distinct categories of subjects: those with severe COPD (Subjects 2 and 5), those with moderate COPD (Subjects 1, 3, and 4), and those who exhibited mild symptoms or were at-risk for developing COPD (Subjects 6, 7, and 8). Interestingly, all subjects within each category displayed similar patterns of ventilation as observed in the SPECT images. Our findings suggest that fSAD could be a crucial factor that contributes to the heterogeneous ventilation observed in patients with COPD. Furthermore, the formation of hot spots in patients with COPD could potentially be influenced by alterations in the segmental airways located in the left upper lobe, along with modifications in the distribution of airflow between the upper and lower lobes, as well as the degree of fSAD and emphysema. Further investigation of fSAD could lead to a better understanding of the diverse characteristics of COPD, improved identification of COPD phenotypes, and enhanced knowledge of the progression of the disease.

Data Availability

The data generated or analyzed during this study are included in this published article and its supplementary information files. The imaging datasets analyzed during the current study are not publicly available due privacy concern but are available from the corresponding author on reasonable request.

DECLARATIONS

Ethics Approval and Consent to Participate

This study (Title: Integrative lung model Y3-Y5) complies with the Declaration of Helsinki. It was approved by the institutional review board of the University of Iowa and informed consent was obtained before the study to have individual patient data published. The manuscript and the results were presented in a way that the patients cannot be identified.

Consent for Publication

All authors approve the version to be published.

Availability of Data and Material

The data generated or analyzed during this study are included in this published article and its supplementary information files. The imaging datasets analyzed during the current study are not publicly available due privacy concern but are available from the corresponding author on reasonable request.

Competing Interests

There is no conflict of interest for all authors.

Funding/Support

Supports for this study were provided, in part, by NIH grants R01-HL168116, U01-HL114494, R01-HL112986, S10-RR022421, T32-HL-144461, and ED grant P116S210005.

Author’s Contributions

All authors made a significant contribution to the work reported, whether that is in the conception, study design, execution, acquisition of data, analysis and interpretation, or in all these areas; took part in drafting, revising or critically reviewing the article; gave final approval of the version to be published; have agreed on the journal to which the article has been submitted; and agree to be accountable for all aspects of the work.

REFERENCES

  1. 1.↵
    Halpin DMG, Vogelmeier CF, Agusti A. Lung Health for All: Chronic Obstructive Lung Disease and World Lung Day 2022. Am J Respir Crit Care Med. 2022;206(6):669–71.
    OpenUrl
  2. 2.↵
    GLOBAL STRATEGY FOR PREVENTION, DIAGNOSIS AND MANAGEMENT OF COPD: 2023 Report [Internet]. 2023. Available from: https://goldcopd.org/2023-gold-report-2/
  3. 3.↵
    Ferrera MC, Labaki WW, Han MK. Advances in Chronic Obstructive Pulmonary Disease. Annu Rev Med. 2021;72:119–34.
    OpenUrl
  4. 4.↵
    Vogelmeier CF, Román-Rodríguez M, Singh D, Han MLK, Rodríguez-Roisin R, Ferguson GT. Goals of COPD treatment: Focus on symptoms and exacerbations. Respir Med. 2020;166(March).
  5. 5.↵
    Israel O, Pellet O, Biassoni L, De Palma D, Estrada-Lobato E, Gnanasegaran G, et al. Two decades of SPECT/CT – the coming of age of a technology: An updated review of literature evidence. Eur J Nucl Med Mol Imaging. 2019;46(10):1990–2012.
    OpenUrl
  6. 6.↵
    Smith BM, Traboulsi H, Austin JHM, Manichaikul A, Hoffman EA, Bleecker ER, et al. Human airway branch variation and chronic obstructive pulmonary disease. Proc Natl Acad Sci U S A. 2018;115(5):E974–81.
    OpenUrlAbstract/FREE Full Text
  7. 7.↵
    Smith BM, Kirby M, Hoffman EA, Kronmal RA, Aaron SD, Allen NB, et al. Association of Dysanapsis with Chronic Obstructive Pulmonary Disease among Older Adults. JAMA - J Am Med Assoc. 2020;323(22):2268–80.
    OpenUrl
  8. 8.↵
    Alford SK, Van Beek EJR, McLennan G, Hoffman EA. Heterogeneity of pulmonary perfusion as a mechanistic image-based phenotype in emphysema susceptible smokers. Proc Natl Acad Sci U S A. 2010;107(16):7485–90.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    Iyer KS, Newell JD, Jin D, Fuld MK, Saha PK, Hansdottir S, et al. Quantitative dual-energy computed tomography supports a vascular etiology of smoking-induced inflammatory lung disease. Am J Respir Crit Care Med. 2016;193(6):652–61.
    OpenUrlCrossRef
  10. 10.↵
    Galbán CJ, Han MK, Boes JL, Chughtai KA, Meyer CR, Johnson TD, et al. Computed tomography-based biomarker provides unique signature for diagnosis of COPD phenotypes and disease progression. Nat Med. 2012;18(11):1711–5.
    OpenUrlCrossRefPubMed
  11. 11.↵
    Ostridge K, Williams N, Kim V, Bennett M, Harden S, Welch L, et al. Relationship between pulmonary matrix metalloproteinases and quantitative CT markers of small airways disease and emphysema in COPD. Thorax. 2016 Feb 1;71(2):126–32.
    OpenUrlAbstract/FREE Full Text
  12. 12.↵
    Singh D. Small airway disease in patients with chronic obstructive pulmonary disease. Tuberc Respir Dis (Seoul). 2017;80(4):317–24.
    OpenUrl
  13. 13.↵
    Haghighi B, Choi S, Choi J, Hoffman EA, Comellas AP, Newell JD, et al. Imaging-based clusters in current smokers of the COPD cohort associate with clinical characteristics: The SubPopulations and Intermediate Outcome Measures in COPD Study (SPIROMICS) 11 Medical and Health Sciences 1102 Cardiorespiratory Medicine and Haemato. Respir Res. 2018 Sep 18;19(1).
  14. 14.↵
    Haghighi B, Choi S, Choi J, Hoffman EA, Comellas AP, Newell JD, et al. Imaging-based clusters in former smokers of the COPD cohort associate with clinical characteristics: The SubPopulations and intermediate outcome measures in COPD study (SPIROMICS). Respir Res. 2019 Jul 15;20(1).
  15. 15.↵
    Zou C, Li F, Choi J, Haghighi B, Choi S, Rajaraman PK, et al. Longitudinal imaging-based clusters in former smokers of the copd cohort associate with clinical characteristics: The subpopulations and intermediate outcome measures in copd study (spiromics). Int J COPD. 2021;16(April):1477–96.
    OpenUrl
  16. 16.↵
    De Backer JW, Vos WG, Vinchurkar SC, Claes R, Drollmann A, Wulfrank D, et al. Validation of computational fluid dynamics in CT-based airway models with SPECT/CT. Radiology. 2010;257(3):854–62.
    OpenUrlCrossRefPubMedWeb of Science
  17. 17.↵
    Xu J, Moonen M, Johansson Å, Gustafsson A, Bake B. Quantitative analysis of inhomogeneity in ventilation SPET. Eur J Nucl Med. 2001;28(12):1795–800.
    OpenUrlCrossRefPubMedWeb of Science
  18. 18.↵
    Doganay O, Matin T, Chen M, Kim M, McIntyre A, McGowan DR, et al. Time-series hyperpolarized xenon-129 MRI of lobar lung ventilation of COPD in comparison to V/Q-SPECT/CT and CT. Eur Radiol. 2019;29(8):4058–67.
    OpenUrl
  19. 19.
    Peiffer JD, Altes T, Ruset IC, Hersman FW, Mugler JP, Meyer CH, et al. Hyperpolarized 129Xe MRI, 99mTc scintigraphy, and SPECT in lung ventilation imaging: a quantitative comparison. Acad Radiol [Internet]. 2023;1–10. Available from: doi:10.1016/j.acra.2023.10.038
    OpenUrlCrossRef
  20. 20.
    Kim M, Doganay O, Matin TN, Povey T, Gleeson F V. CT-based Airway Flow Model to Assess Ventilation in Chronic Obstructive Pulmonary Disease: A pilot study. Radiology. 2019;293(3):666–73.
    OpenUrl
  21. 21.↵
    Bayat S, Wild J, Winkler T. Lung functional imaging. Breathe [Internet]. 2023;19(3). Available from: doi:10.1183/20734735.0272-2022
    OpenUrlAbstract/FREE Full Text
  22. 22.↵
    Fuld MK, Grout RW, Guo J, Morgan JH, Hoffman EA. Systems for Lung Volume Standardization during Static and Dynamic MDCT-based Quantitative Assessment of Pulmonary Structure and Function. Acad Radiol [Internet]. 2012;19(8):930–40. Available from: doi:10.1016/j.acra.2012.03.017
    OpenUrlCrossRefPubMed
  23. 23.↵
    Haghighi B, D. Ellingwood N, Yin Y, Hoffman EA, Lin CL. A GPU-based symmetric non-rigid image registration method in human lung. Med Biol Eng Comput. 2018;56(3):355–71.
    OpenUrl
  24. 24.↵
    Krogsgaard OW. Technetium-99m-Sulfur Colloid. 1976;5.
  25. 25.↵
    Sirr SA, Juenemann PJ, Tom H, Boudreau RJ, Chandler RP, Loken MK. Effect of ethanol on droplet size, efficiency of delivery, and clearance characteristics of technetium-99m DTPA aerosol. J Nucl Med. 1985;26(6):643–6.
    OpenUrlAbstract/FREE Full Text
  26. 26.↵
    Nagao M, Murase K. Measurement of heterogeneous distribution on Technegas SPECT images by three-dimensional fractal analysis. Ann Nucl Med. 2002;16(6):369–76.
    OpenUrlPubMed
  27. 27.↵
    Horn J. Factors in Factor Analysis. Psychometrika. 1965;30(2):179–85.
    OpenUrlCrossRefPubMedWeb of Science
  28. 28.↵
    Kaiser HF. The varimax criterion for analytic rotation in factor analysis. Psychometrika. 1958;23(3):187–200.
    OpenUrlCrossRefWeb of Science
  29. 29.↵
    Akoglu H. User’s guide to correlation coefficients. Turkish J Emerg Med. 2018;18(3):91–3.
    OpenUrlCrossRef
  30. 30.↵
    Han MLK, Agusti A, Calverley PM, Celli BR, Criner G, Curtis JL, et al. Chronic obstructive pulmonary disease phenotypes: The future of COPD. Am J Respir Crit Care Med. 2010;182(5):598–604.
    OpenUrlCrossRefPubMedWeb of Science
  31. 31.↵
    O’Donnell DE, Laveneziana P. Physiology and consequences of lung hyperinflation in COPD. Eur Respir Rev. 2006;15(100):61–7.
    OpenUrlAbstract/FREE Full Text
  32. 32.↵
    Marsh S, Aldington S, Williams M V., Nowitz M, Kingzett-Taylor A, Weatherall M, et al. Physiological associations of computerized tomography lung density: a factor analysis. Int J Chron Obstruct Pulmon Dis. 2006;1(2):181–7.
    OpenUrlPubMed
  33. 33.
    Kim CS, Kang TC. Comparative measurement of lung deposition of inhaled fine particles in normal subjects and patients with obstructive airway disease. Am J Respir Crit Care Med [Internet]. 1997 Mar;155(3):899–905. Available from: https://www.atsjournals.org/doi/10.1164/ajrccm.155.3.9117024
    OpenUrl
  34. 34.↵
    Zhang X, Li F, Rajaraman PK, Choi J, Comellas AP, Hoffman EA, et al. A computed tomography imaging-based subject-specific whole-lung deposition model. Eur J Pharm Sci [Internet]. 2022;177(July):106272. Available from: doi:10.1016/j.ejps.2022.106272
    OpenUrlCrossRef
Back to top
Posted April 13, 2024.
Download PDF

Supplementary Material

Data/Code
Email
Exploratory Study on COPD Phenotypes and their Progression: Integrating SPECT and qCT Imaging Analysis
Frank Li, Xuan Zhang, Alejandro P. Comellas, Eric A. Hoffman, Michael M. Graham, Ching-Long Lin
medRxiv 2024.04.10.24305577; doi: https://doi.org/10.1101/2024.04.10.24305577
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area