Structural and Functional Features on Quantitative Chest Computed Tomography in the Korean Asian versus the White American Healthy Non-Smokers

Human-airway surface mesh smoothing based on graph convolutional neural networks

BACKGROUND AND OBJECTIVE

A detailed representation of the airway geometry in the respiratory system is critical for predicting precise airflow and pressure behaviors in computed tomography (CT)-image-based computational fluid dynamics (CFD). The CT-image-based geometry often contains artifacts, noise, and discontinuities due to the so-called stair step effect. Hence, an advanced surface smoothing is necessary. The existing smoothing methods based on the Laplacian operator drastically shrink airway geometries, resulting in the loss of information related to smaller branches. This study aims to introduce an unsupervised airway-mesh-smoothing learning (AMSL) method that preserves the original geometry of the three-dimensional (3D) airway for accurate CT-image-based CFD simulations.

METHOD

The AMSL method jointly trains two graph convolutional neural networks (GCNNs) defined on airway meshes to filter vertex positions and face normal vectors. In addition, it regularizes a combination of loss functions such as reproducibility, smoothness and consistency of vertex positions, and normal vectors. The AMSL adopts the concept of a deep mesh prior model, and it determines the self-similarity for mesh restoration without using a large dataset for training. Images of the airways of 20 subjects were smoothed by the AMSL method, and among them, the data of two subjects were used for the CFD simulations to assess the effect of airway smoothing on flow properties.

RESULTS

In 18 of 20 benchmark problems, the proposed smoothing method delivered better results compared with the conventional or state-of-the-art deep learning methods. Unlike the traditional smoothing, the AMSL successfully constructed 20 smoothed airways with airway diameters that were consistent with the original CT images. Besides, CFD simulations with the airways obtained by the AMSL method showed much smaller pressure drop and wall shear stress than the results obtained by the traditional method.

CONCLUSIONS

The airway model constructed by the AMSL method reproduces branch diameters accurately without any shrinkage, especially in the case of smaller airways. The accurate estimation of airway geometry using a smoothing method is critical for estimating flow properties in CFD simulations.

A reference equation for lung volume on computed tomography in Japanese middle-aged and elderly adults

BACKGROUND

Effective use of lung volume data measured on computed tomography (CT) requires reference values for specific populations. This study examined whether an equation previously generated for multiple ethnic groups in the United States, including Asians predominantly composed of Chinese people, in the Multi-Ethnic Study of Atherosclerosis (MESA) could be used for Japanese people and, if necessary, to optimize this equation. Moreover, the equation was used to characterize patients with chronic obstructive pulmonary disease (COPD) and lung hyperexpansion.

METHODS

This study included a lung cancer screening CT cohort of asymptomatic never smokers aged ≥40 years from two institutions (n = 364 and 419) to validate and optimize the MESA equation and a COPD cohort (n = 199) to test its applicability.

RESULTS

In all asymptomatic never smokers, the variance explained by the predicted values (R2) based on the original MESA equation was 0.60. The original equation was optimized to minimize the root mean squared error (RMSE) by adjusting the scaling factor but not the age, sex, height, or body mass index terms of the equation. The RMSE changed from 714 ml in the original equation to 637 ml in the optimized equation. In the COPD cohort, lung hyperexpansion, defined based on the 95th percentile of the ratio of measured lung volume to predicted lung volume in never smokers (122 %), was observed in 60 (30 %) patients and was associated with centrilobular emphysema and air trapping on inspiratory/expiratory CT.

CONCLUSIONS

The MESA equation was optimized for Japanese middle-aged and elderly adults.

Computer-assisted evaluation of small airway disease in CT scans of Iran-Iraq war victims of chemical warfare by a locally developed software: comparison between different quantitative methods

OBJECTIVE

Diagnosis of small airway disease on computed tomography (CT) scans is challenging in patients with a history of chemical warfare exposure. We developed a software package based on different methodologies to identify and quantify small airway disease in CT images. The primary aim was to identify the best automatic methodology for detecting small airway disease in CT scans of Iran-Iraq War victims of chemical warfare.

METHODS

This retrospective case-control study enrolled 46 patients with a history of chemical warfare exposure and 27 controls with inspiratory/expiratory (I/E) CT scans and spirometry tests. Image data were automatically segmented, and inspiratory images were registered into the expiratory images' frame using the locally developed software. Parametric response mapping (PRM) and air trapping index (ATI) mapping were performed on the CT images. Conventional QCT methods, including expiratory/inspiratory mean lung attenuation (E/I MLA) ratio, normal density E/I (ND E/I) MLA ratio, attenuation volume Index (AVI), %low attenuation areas (LAA) < -856 in exhale scans, and %LAA < -950 in inhale scans were also computed. QCT measurements were correlated with spirometry results and compared across the two study groups.

RESULTS

The correlation analysis showed a significant negative relationship between three air trapping (AT) measurements (PRM, ATI, and %LAAExp < -856) and spirometry parameters (Fev1, Fvc, Fev1/Fvc, and MMEF). Moreover, %LAAExp < -856 had the highest significant negative correlation with Fev1/Fvc (r = -0.643, P-value < 0.001). Three AT measurements demonstrated a significant difference between the study groups. The E/I ratio was also significantly different between the two groups (P-value < 0.001). Binary logistic regression models showed PRMFsad, %LAAExp < -856, and ATI as significant and strong predictors of the study outcome. Optimal cut-points for PRMFsad = 19%, %LAAExp < -856 = 23%, and ATI = 27% were identified to classify the participants into two groups with high accuracy.

CONCLUSION

QCT methods, including PRM, ATI, and %LAAExp < -856 can greatly advance the identification and quantification of SAD in chemical warfare victims. The results should be verified in well-designed prospective studies involving a large population.

An unsupervised image registration method employing chest computed tomography images and deep neural networks

BACKGROUND

Deformable image registration is crucial for multiple radiation therapy applications. Fast registration of computed tomography (CT) lung images is challenging because of the large and nonlinear deformation between inspiration and expiration. With advancements in deep learning techniques, learning-based registration methods are considered efficient alternatives to traditional methods in terms of accuracy and computational cost.

METHOD

In this study, an unsupervised lung registration network (LRN) with cycle-consistent training is proposed to align two acquired CT-derived lung datasets during breath-holds at inspiratory and expiratory levels without utilizing any ground-truth registration results. Generally, the LRN model uses three loss functions: image similarity, regularization, and Jacobian determinant. Here, LRN was trained on the CT datasets of 705 subjects and tested using 10 pairs of public CT DIR-Lab datasets. Furthermore, to evaluate the effectiveness of the registration technique, target registration errors (TREs) of the LRN model were compared with those of the conventional algorithm (sum of squared tissue volume difference; SSTVD) and a state-of-the-art unsupervised registration method (VoxelMorph).

RESULTS

The results showed that the LRN with an average TRE of 1.78 ± 1.56 mm outperformed VoxelMorph with an average TRE of 2.43 ± 2.43 mm, which is comparable to that of SSTVD with an average TRE of 1.66 ± 1.49 mm. In addition, estimating the displacement vector field without any folding voxel consumed less than 2 s, demonstrating the superiority of the learning-based method with respect to fiducial marker tracking and the overall soft tissue alignment with a nearly real-time speed.

CONCLUSIONS

Therefore, this proposed method shows significant potential for use in time-sensitive pulmonary studies, such as lung motion tracking and image-guided surgery.

Structural and functional alterations of subjects with cement dust exposure: A longitudinal quantitative computed tomography-based study

Cement dust exposure (CDE) can be a risk factor for pulmonary disease, causing changes in segmental airways and parenchymal lungs. This study investigates longitudinal alterations in quantitative computed tomography (CT)-based metrics due to CDE. We obtained CT-based airway structural and lung functional metrics from CDE subjects with baseline CT and follow-up CT scans performed three years later. From the CT, we extracted wall thickness (WT) and bifurcation angle (θ) at total lung capacity (TLC) and functional residual capacity (FRC), respectively. We also computed air volume (V_air), tissue volume (V_tissue), global lung shape, percentage of emphysema (Emph%), and more. Clinical measures were used to associate with CT-based metrics. Three years after their baseline, the pulmonary function tests of CDE subjects were similar or improved, but there were significant alterations in the CT-based structural and functional metrics. The follow-up CT scans showed changes in θ at most of the central airways; increased WT at the subgroup bronchi; smaller V_air at TLC at all except the right upper and lower lobes; smaller V_tissue at all lobes in TLC and FRC except for the upper lobes in FRC; smaller global lung shape; and greater Emph% at the right upper and lower lobes. CT-based structural and functional variables are more sensitive to the early identification of CDE subjects, while most clinical lung function changes were not noticeable. We speculate that the significant long-term changes in CT are uniquely observed in CDE subjects, different from smoking-induced structural changes.

Bronchial wall parameters on CT in healthy never-smoking, smoking, COPD, and asthma populations: a systematic review and meta-analysis

Objective

Research on computed tomography (CT) bronchial parameter measurements shows that there are conflicting results on the values for bronchial parameters in the never-smoking, smoking, asthma, and chronic obstructive pulmonary disease (COPD) populations. This review assesses the current CT methods for obtaining bronchial wall parameters and their comparison between populations.

Methods

A systematic review of MEDLINE and Embase was conducted following PRISMA guidelines (last search date 25th October 2021). Methodology data was collected and summarised. Values of percentage wall area (WA%), wall thickness (WT), summary airway measure (Pi10), and luminal area (Ai) were pooled and compared between populations.

Results

A total of 169 articles were included for methodologic review; 66 of these were included for meta-analysis. Most measurements were obtained from multiplanar reconstructions of segmented airways (93 of 169 articles), using various tools and algorithms; third generation airways in the upper and lower lobes were most frequently studied. COPD (12,746) and smoking (15,092) populations were largest across studies and mostly consisted of men (median 64.4%, IQR 61.5 – 66.1%). There were significant differences between populations; the largest WA% was found in COPD (mean SD 62.93 ± 7.41%, n = 6,045), and the asthma population had the largest Pi10 (4.03 ± 0.27 mm, n = 442). Ai normalised to body surface area (Ai/BSA) (12.46 ± 4 mm², n = 134) was largest in the never-smoking population.

Conclusions

Studies on CT-derived bronchial parameter measurements are heterogenous in methodology and population, resulting in challenges to compare outcomes between studies. Significant differences between populations exist for several parameters, most notably in the wall area percentage; however, there is a large overlap in their ranges.

Key Points

• Diverse methodology in measuring airways contributes to overlap in ranges of bronchial parameters among the never-smoking, smoking, COPD, and asthma populations.

• The combined number of never-smoking participants in studies is low, limiting insight into this population and the impact of participant characteristics on bronchial parameters.

• Wall area percent of the right upper lobe apical segment is the most studied (87 articles) and differentiates all except smoking vs asthma populations.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00330-022-08600-1.

Quantitative computed tomography imaging-based classification of cement dust-exposed subjects with an artificial neural network technique

BACKGROUND AND OBJECTIVE

Cement dust exposure is likely to affect the structural and functional alterations in segmental airways and parenchymal lungs. This study develops an artificial neural network (ANN) model for identifying cement dust-exposed (CDE) subjects using quantitative computed tomography-based airway structural and functional features.

METHODS

We obtained the airway features in five central and five sub-grouped segmental regions and the lung features in five lobar regions and one total lung region from 311 CDE and 298 non-CDE (NCDE) subjects. The five-fold cross-validation method was adopted to train the following classification models:ANN, support vector machine (SVM), logistic regression (LR), and decision tree (DT). For all the classification models, linear discriminant analysis (LDA) and genetic algorithm (GA) were applied for dimensional reduction and hyperparameterization, respectively. The ANN model without LDA was also optimized by the GA method to observe the effect of the dimensional reduction.

RESULTS

The genetically optimized ANN model without the LDA method was the best in terms of the classification accuracy. The accuracy, sensitivity, and specificity of the GA-ANN model with four layers were greater than those of the other classification models (i.e., ANN, SVM, LR, and DT using LDA and GA methods) in the five-fold cross-validation. The average values of accuracy, sensitivity, and specificity for the five-fold cross-validation were 97.0%, 98.7%, and 98.6%, respectively.

CONCLUSIONS

We demonstrated herein that a quantitative computed tomography-based ANN model could more effectively detect CDE subjects when compared to their counterpart models. By employing the model, the CDE subjects may be identified early for therapeutic intervention.

Generation-based study of airway remodeling in smokers with normal-looking CT with normalization to control inter-subject variability

PURPOSE

With the help of quantitative computed tomography (QCT), it is possible to identify smoking-associated airway remodeling. However, there is currently little information on whether QCT-based airway metrics are sensitive to early airway wall remodeling in subclinical phases of smoking-associated airway disease. This study aimed to evaluate a predictive model that normalized airway parameters and investigate structural airway alterations in smokers with normal-looking CT using the normalization scheme.

METHODS

In this retrospective analysis, 222 non-smokers (male 97, female 125) and 69 smokers (male 66, female 3) from January 2014 to December 2016 were included, and airway parameters were quantitatively analyzed. To control inter-subject variability, multiple linear regressions of tracheal wall thickness (WT), diameter (D), and luminal area (LA) were performed, adjusted for age, sex, and height. Using this normalization scheme, airway parameters with matched generation were compared between smokers and non-smokers.

RESULTS

Using the normalization scheme, it was possible to assess generation-based structural alterations of the airways in subclinical smokers. Smokers showed diffuse luminal narrowing of airways for most generations (P < 0.05, except 3rd generation), no change in wall thickness of the proximal bronchi (1st-3rd generation), and a thinning of distal airways (P <0.05, ≥4th generation).

CONCLUSION

QCT assessment for subclinical smokers can help identify minimal structural changes in airways induced by smoking.

Quantitative CT Analysis in Patients with Pulmonary Emphysema: Do Calculated Differences Between Full Inspiration and Expiration Correlate with Lung Function?

Purpose

The aim of this retrospective study was to evaluate correlations between parameters of quantitative computed tomography (QCT) analysis, especially the 15th percentile of lung attenuation (P15), and parameters of clinical tests in a large group of patients with pulmonary emphysema.

Patients and Methods

One hundred and seventy-two patients with pulmonary emphysema and chronic obstructive pulmonary disease (COPD) global initiative for chronic obstructive lung disease (GOLD) stage 3 or 4 were assessed by nonenhanced thin-section CT scans in full inspiratory and expiratory breath-hold, pulmonary function test (PFT), a 6-minute walk test (6MWT), and quality of life questionnaires (SGRQ and CAT). QCT parameters included total lung volume (TLV), total emphysema score (TES), and P15, all measured at inspiration (IN) and expiration (EX). Differences between inspiration and expiration were calculated for TLV (TLV_Diff), TES (TES_Diff), and P15 (P15_Diff). Spearman correlation analysis was performed.

Results

CT-measured lung volume in inspiration (TLV_IN) correlated strongly with spirometry-measured total lung capacity (TLC) (r=0.81, p<0.001) and moderately to strongly with residual volume (RV), forced vital capacity (FVC), and forced expiratory volume in 1 second (FEV₁)/FVC (r=0.60, 0.56, and −0.49, each p<0.001). Lung volume in expiration (TLV_EX) correlated moderately to strongly with TLC, RV and FEV₁/FVC ratio (r=0.75, 0.66, and −0.43, each p<0.001). TES and P15 showed stronger correlations with the carbon monoxide transfer coefficient (KCO%) (r= −0.42, 0.44, both p<0.001), when measured during expiration. P15_Diff correlated moderately with KCO% and carbon monoxide diffusing capacity (DLCO%) (r= 0.41, 0.40, both p<0.001). The 6MWT and most QCT parameters showed significant differences between COPD GOLD 3 and 4 groups.

Conclusion

Our results suggest that QCT can help predict the severity of lung function decrease in patients with pulmonary emphysema and COPD GOLD 3 or 4. Some QCT parameters, including P15_EX and P15_Diff, correlated moderately to strongly with parameters of pulmonary function tests.

From infancy to adulthood-Developmental changes in pulmonary quantitative computed tomography parameters

Purpose

Quantified computed tomography (qCT) is known for correlations with airflow obstruction and fibrotic changes of the lung. However, as qCT studies often focus on diseased and elderly subjects, current literature lacks physiological qCT values during body development. We evaluated chest CT examinations of a healthy cohort, reaching from infancy to adulthood, to determine physiological qCT values and changes during body development.

Method

Dose-optimized chest CT examinations performed over the last 3 years using a dual-source CT were retrospectively analysed. Exclusion criteria were age >30 years and any known or newly diagnosed lung pathology. Lung volume, mean lung density, full-width-at-half-maximum and low attenuated volume (LAV) were semi-automated quantified in 151 patients. qCT values between different age groups as well as unenhanced (Group 1) and contrast-enhanced (Group 2) protocols were compared. Models for projection of age-dependant changes in qCT values were fitted.

Results

Significant differences in qCT parameters were found between the age groups from 0 to 15 years (p < 0.05). All parameters except LAV merge into a plateau level above this age as shown by polynomial models (r2 between 0.85 and 0.67). In group 2, this plateau phase is shifted back around five years. Except for the volume, significant differences in all qCT values were found between group 1 and 2 (p < 0.01).

Conclusion

qCT parameters underly a specific age-dependant dynamic. Except for LAV, qCT parameters reach a plateau around adolescence. Contrast-enhanced protocols seem to shift this plateau backwards.

Quantitative CT-based structural alterations of segmental airways in cement dust-exposed subjects

Background

Dust exposure has been reported as a risk factor of pulmonary disease, leading to alterations of segmental airways and parenchymal lungs. This study aims to investigate alterations of quantitative computed tomography (QCT)-based airway structural and functional metrics due to cement-dust exposure.

Methods

To reduce confounding factors, subjects with normal spirometry without fibrosis, asthma and pneumonia histories were only selected, and a propensity score matching was applied to match age, sex, height, smoking status, and pack-years. Thus, from a larger data set (N = 609), only 41 cement dust-exposed subjects were compared with 164 non-cement dust-exposed subjects. QCT imaging metrics of airway hydraulic diameter (D_h), wall thickness (WT), and bifurcation angle (θ) were extracted at total lung capacity (TLC) and functional residual capacity (FRC), along with their deformation ratios between TLC and FRC.

Results

In TLC scan, dust-exposed subjects showed a decrease of D_h (airway narrowing) especially at lower-lobes (p < 0.05), an increase of WT (wall thickening) at all segmental airways (p < 0.05), and an alteration of θ at most of the central airways (p < 0.001) compared with non-dust-exposed subjects. Furthermore, dust-exposed subjects had smaller deformation ratios of WT at the segmental airways (p < 0.05) and θ at the right main bronchi and left main bronchi (p < 0.01), indicating airway stiffness.

Conclusions

Dust-exposed subjects with normal spirometry demonstrated airway narrowing at lower-lobes, wall thickening at all segmental airways, a different bifurcation angle at central airways, and a loss of airway wall elasticity at lower-lobes. The airway structural alterations may indicate different airway pathophysiology due to cement dusts.