3D computed tomography airway geometry for predicting bronchoscopic accessibility in peripheral pulmonary nodules: a prospective study
Highlight box
Key findings
• Quantitative computed tomography-derived airway geometric features—specifically bifurcation angle, curvature, and luminal narrowing—are independent predictors of bronchoscopic accessibility for peripheral pulmonary nodules.
• A composite model integrating these 3D geometric metrics with clinical variables demonstrated improved predictive performance (area under the curve =0.84), offering incremental value beyond models based on conventional clinical or radiologic features alone.
What is known and what is new?
• Bronchoscopic accessibility is a prerequisite for successful radial probe endobronchial ultrasound-guided biopsy, yet its prediction has traditionally relied on subjective operator intuition or simple, non-quantitative bronchus sign assessments.
• This study establishes a data-driven framework that converts 3D airway morphology into quantifiable geometric biomarkers. By integrating metrics like bifurcation angles and curvature, this approach provides reproducible navigation mapping that replaces subjective qualitative assessments.
What is the implication, and what should change now?
• Preprocedural planning using quantitative airway analysis enables early identification of challenging pathways, allowing clinicians to optimize bronchoscopic strategies and select appropriate navigational tools.
• Integrating these objective geometric assessments into routine clinical practice can reduce operator-dependent bias and improve the overall efficiency and safety of peripheral lung nodule biopsies.
Introduction
Diagnosis of peripheral pulmonary nodules (PPNs) through bronchoscopy has been extensively investigated in recent literature (1). Bronchoscopic interventions, such as radial probe endobronchial ultrasound (rEBUS)-guided transbronchial lung biopsy (TBLB) based on computed tomography (CT)-guided approaches, have been extensively explored (2). Preprocedural chest CT imaging enables precise procedural planning by delineating key anatomical structures, such as lesion location and bronchial pathways, whereas rEBUS provides intraprocedural real-time localization of the lesion relative to the probe, enhancing the accuracy of imaging-guided bronchoscopic interventions when other forms of intraprocedural imaging are unavailable (3,4).
However, PPN identification with rEBUS depends heavily on operator knowledge of bronchial anatomy and procedural difficulty may vary according to the anatomical course of the target bronchus (5,6). Certain bronchial pathways may present greater mechanical challenges during scope advancement, underscoring the importance of additional preprocedural information for procedural planning. When a bronchoscope is advanced into an upper lobar bronchus, the direction of pushing force applied by the operator is opposite to the direction in which the bronchoscope needs to proceed (6). Consequently, accessibility to PPNs in the upper lobes can be lower than in other lobes. These anatomical constraints underscore the need for additional preprocedural information to support bronchoscopic navigation.
Lung airway morphometry using chest CT has emerged as a powerful tool for quantifying airway geometry and assessing structural abnormalities (7). Prior research has proposed CT-derived quantitative airway measures as predictors of lung function decline and disease progression (8-10). However, those findings have largely been limited to observational studies, and their utility in guiding real-world clinical procedures remains underexplored. In our previous retrospective study, we demonstrated the feasibility of CT-based airway analysis in preoperative planning (11). Although our results suggested potential clinical value, their application in prospective clinical settings had not been validated.
Therefore, in this study, we extend our previous work by implementing a prospective design to develop and clinically validate a predictive model for estimating the accessibility of PPNs during rEBUS-TBLB. We further assessed the utility of our predictive model as a preprocedural guide in routine bronchoscopic workflows. We present this article in accordance with the TRIPOD reporting checklist (available at https://tlcr.amegroups.com/article/view/10.21037/tlcr-2026-1-0133/rc).
Methods
Study participants
The study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. The study was approved by the Institutional Review Board of Samsung Medical Center (No. 2018-03-021). Informed consent was taken from all the patients (or a statement that it was not required and why). Patient information was de-identified and anonymized prior to the analysis.
We collected data from 212 patients whose rEBUS-TBLB procedures were planned during a prospective trial at Samsung Medical Center, a single tertiary referral center, between December 2021 and February 2023. Since the procedures were performed using rEBUS alone—without additional guidance from fluoroscopy, cone-beam CT, guide sheaths, or navigation systems—patients were primarily referred for the procedure when a bronchus leading to the lesion was identifiable on pre-procedural CT to ensure procedural feasibility. During the initial clinical screening, and prior to any quantitative airway analysis, 10 patients were excluded for the following reasons: the operator determined that the lesion was clinically unsuitable for biopsy (n=6), the lesions were small pure ground-glass opacities (n=3), or the lesion improved on follow-up imaging before the procedure (n=1). Quantitative airway analysis was then attempted in the remaining 202 patients. However, three patients were further excluded because of processing failures due to technical or software issues. Consequently, a total of 219 PPNs from 199 patients were included in the final analysis. Because several patients presented with multiple PPNs, each lesion was prospectively assessed and included in the analysis to reflect the full clinical spectrum encountered in real-world practice. The overall study flow is illustrated in Figure 1.
CT acquisition and preoperative evaluation
All patients underwent chest CT scanning prior to rEBUS-TBLB as part of the prospective study protocol. CT scans were acquired with standard thin slices (≤1 mm) and high-resolution reconstruction; detailed scanning parameters are available in a supplementary appendix online (Appendix 1).
Several radiologic characteristics were visually assessed on chest CT: lesion size, lesion depth from the pleural surface, number of bronchial branches from the trachea to the lesion, and lesion solidity (solid, sub-solid, or cavitary). The lesion location was categorized into five anatomical regions: right upper, right middle, right lower, left upper, and left lower lobes. Furthermore, to account for anatomical and physiological factors affecting procedural difficulty, we categorized the segments into three bronchial clusters: (I) segments requiring acute angulation for access, specifically the apical and superior segments (RB1, LB1, RB6, LB6); (II) segments with high respiratory-induced mobility, including the basal segments (RB8–10, LB8–10); and (III) all other segments. This classification reflects established challenges in bronchoscopy, such as the acute angles of the upper lobes (6) and the significant respiratory displacement common in the lower lobes (12). All CT image assessments were performed by a board-certified thoracic radiologist (H.Y.L.) with 18 years of experience in thoracic imaging.
To assess preoperative airway information, we performed three-dimensional airway segmentation and lesion localization using chest CT data and commercial deep-learning software (AVIEW; Coreline Soft, Seoul, Korea). The segmented airway and lesion volumes were reviewed by an expert radiologist (H.Y.L.). Quantitative airway geometrical features were extracted using in-house software (AirGeo) that was developed and examined in a previous study (11). The extracted features included cross-sectional dimensions (e.g., minimum diameter, perimeter, circularity), spatial orientation (e.g., curvature, torsion), and branching characteristics (e.g., bifurcation angle). The detailed quantitative analysis process and a comprehensive list of airway geometry variables are provided in Appendix 1 and Table S1.
The extracted CT-based airway information was provided to the bronchoscopist prior to the procedure, enabling preoperative planning. It is important to note that the AirGeo-based analysis provided only raw quantitative data (e.g., subsegmental routes and bronchial angles) for preoperative planning. It did not offer any qualitative assessment of predicted difficulty. Thus, it served as an informative tool to guide the anatomical approach, rather than as a behavioral intervention. The proceduralist’s effort was standardized across all cases, and the study’s objective was to determine which of these objective preoperative parameters were significantly associated with the actual technical difficulty encountered during the procedure. Figure 2 shows the overall workflow of our study.
Procedure for rEBUS-TBLB
All rEBUS-TBLB procedures were performed on patients under conscious sedation using midazolam and fentanyl, as previously reported (13). Briefly, a 4.1-mm bronchoscope (BF P260F; Olympus, Tokyo, Japan) was advanced to the sub-subsegmental bronchus nearest the lesion after reviewing the CT. Then, the rEBUS probe (1.4-mm, 20-MHz, UM S20-17S; Olympus) was inserted through the bronchoscope’s working channel. Upon confirming lesion visualization via ultrasound, a 1.8-mm biopsy forceps was inserted through the working channel for TBLB after the probe was removed. All procedures were performed by one operator (B.H.J.), a bronchoscopist with 11 years of clinical experience, without a guide sheath, fluoroscopy, or navigation system.
When the target was identified on ultrasound, the findings were classified as “within” (the radial probe was within the target lesion, surrounded by lesion material) or “adjacent” (the probe did not pierce the lesion but was nearby). When the nodule was not able to be visualized by the probe, it was defined as “invisible”, indicating failure of endobronchial navigation using thin bronchoscopy.
Independent of the relationship between the probe and the lesion, the difficulty in reaching the lesion was defined by the operator as “easily accessible”, “difficult but accessible”, or “inaccessible”. To provide a more objective assessment, “easily accessible” was defined as reaching the target lesion within one or two attempts after identifying the bronchial path on preprocedural CT. “Difficult” access (including “difficult but accessible” and “inaccessible”) was defined as cases requiring more than three attempts (>3 attempts) to successfully navigate the probe to the target.
We analyzed 219 procedures from 199 patients (some with multiple lesions) to validate the prediction model. Additionally, the procedure time was recorded from radial probe insertion to biopsy completion to assess procedural efficiency.
Statistical analysis
In this study, we focused on airway geometric characteristics as three levels of feature sets. First, six global level features were defined using five first-order statistics about the bifurcation angles of the total branches and the total length of each sectional branch. Next, 64 features were considered from the last two branches leading to the target lesion because they are likely most relevant for accessing the lesion. All cases were grouped by final accessibility (easily accessible vs. difficult or inaccessible).
For the prediction model, multicollinearity among the extracted features was assessed using the variance inflation factor (VIF), excluding features with VIF >4 (14). Univariate logistic regression identified variables associated with accessibility, and redundant variables (P<0.05) were removed. Feature selection and predictive scoring were performed using logistic-least absolute shrinkage and selection operator (LASSO) with 10-fold cross-validation to determine the optimal penalty term. The predictive score for each patient was calculated as a weighted sum of the selected features.
We developed three prediction models: clinical, airway geometrical, and composite (combining both feature sets). Model performance was evaluated using receiver operating characteristic (ROC) curves and corresponding area under the curve (AUC) values. DeLong’s test was applied to compare model classification performance (15). Calibration was assessed by plotting observed versus predicted probabilities, where an ideal nomogram would closely follow the 45-degree reference line (16).
Normality of continuous variables was assessed using the Shapiro-Wilk test. Because most continuous variables were not normally distributed, they are presented as median [interquartile range (IQR)] and were compared using the Mann-Whitney U test. Categorical variables were compared using the Chi-squared test. All statistical analyses were performed using both SPSS version 20.0.0.0 (IBM Corp., Armonk, NY, USA) and R 3.2.5 software (R Foundation for Statistical Computing, Vienna, Austria). A P value less than 0.05 was considered statistically significant.
Results
Characteristics of cohorts
Demographics for the enrolled patients are shown in Table 1. The patients [median age, 70.0 years (IQR, 62.0–77.0 years); 118 male and 81 female] had lesions with three types of solidity: solid (n=150, 68.5%), sub-solid (n=47, 21.5%), and cavitary (n=22, 10.0%). Compared with easily accessible cases, difficult or inaccessible cases had significantly smaller target lesions [19.2 mm (IQR, 15.6–22.8 mm) vs. 22.8 mm (IQR, 17.7–29.5 mm), P=0.008]. There was no significant difference between groups in the number of branch levels (P=0.56). Although lesion depth from the pleural surface did not reach statistical significance, it was numerically smaller in difficult or inaccessible cases [4.3 mm (IQR, 0–10.4 mm) vs. 6.9 mm (IQR, 0–16.9 mm), P=0.08]. The distribution of the PPNs in the five lobes did not differ statistically between the groups. However, in terms of bronchus clustering, more difficult or inaccessible cases than easily accessible cases tended to be found in the RB1/LB1/RB6/LB6 cluster (45.9% vs. 25.8%, P=0.059). The ‘invisible’ bronchus sign was significantly more frequent in the difficult or inaccessible group than the easily accessible group (16.2% vs. 3.3%; P=0.007). While the “within” was more frequently observed in the easily accessible group (155/182, 85.2%) than in the difficult or inaccessible group (27/37, 73.0%), its presence in a substantial proportion of the difficult or inaccessible group cases was notable.
Table 1
| Variables | Total cases (n=219) | Easily accessible (n=182) | Difficult or inaccessible (n=37) | P |
|---|---|---|---|---|
| Age (years)† | 70.0 [62.0–77.0] | 70.0 [62.0–77.0] | 68.0 [62.0–75.0] | 0.71 |
| Male† | 118/199 (59.3) | 98/166 (59.0) | 20/33 (60.6) | 0.58 |
| Lesion size (mm) | 22.4 [17.5–28.0] | 22.8 [17.7–29.5] | 19.2 [15.6–22.8] | 0.008 |
| Lesion depth from pleural surface (mm) | 6.8 [0.0–15.8] | 6.9 [0–16.9] | 4.3 [0–10.4] | 0.08 |
| Solidity | 0.91 | |||
| Solid | 150 (68.5) | 124 (68.1) | 26 (70.3) | |
| Sub-solid | 47 (21.5) | 39 (21.4) | 8 (21.6) | |
| Cavitary | 22 (10.0) | 19 (10.4) | 3 (8.1) | |
| Number of branch levels | 7.0 [6.0–8.0] | 7.0 [6.0–8.0] | 7.0 [6.0–8.0] | 0.56 |
| Lesion location | 0.12 | |||
| Right upper | 47 (21.5) | 37 (20.3) | 10 (27.0) | |
| Right middle | 16 (7.3) | 16 (8.8) | 0 | |
| Right lower | 47 (21.5) | 35 (19.2) | 12 (32.4) | |
| Left upper | 69 (31.5) | 59 (32.4) | 10 (27.0) | |
| Left lower | 40 (18.3) | 35 (19.2) | 5 (13.5) | |
| Bronchus clustering | 0.059 | |||
| RB1, LB1, RB6, LB6 | 64 (29.2) | 47 (25.8) | 17 (45.9) | |
| RB8–10, LB8–10 | 57 (26.0) | 50 (27.5) | 7 (18.9) | |
| Others | 98 (44.7) | 85 (46.7) | 13 (35.1) | |
| Bronchus sign on CT | 0.007 | |||
| Within | 182 (83.1) | 155 (85.2) | 27 (73.0) | |
| Adjacent | 25 (11.4) | 21 (11.5) | 4 (10.8) | |
| Invisible | 12 (5.5) | 6 (3.3) | 6 (16.2) |
Data are reported as median [interquartile range] or number (%). P values for continuous variables were calculated using the Mann-Whitney U test, and P values for categorical variables were calculated using the Chi-squared test. We performed rEBUS-TBLB for 219 PPNs in 199 patients (19 patients had multiple PPNs). †, age and sex were calculated for 199 patients. All remaining variables were analyzed for 219 PPNs. CT, computed tomography; LB1, apical segmental bronchus of the left upper lobar bronchus; LB6, superior segmental bronchus of the left lower lobar bronchus; LB8–10, anterior, lateral, and posterior segmental bronchi of the left lower lobar bronchus; PPN, peripheral pulmonary nodule; RB1, apical segmental bronchus of the right upper lobar bronchus; RB6, superior segmental bronchus of the right lower lobar bronchus; RB8–10, anterior, lateral, and posterior segmental bronchi of the right lower lobar bronchus; rEBUS-TBLB, radial probe endobronchial ultrasound-guided transbronchial lung biopsy.
Procedure results are shown in Table S2. The procedure time was longer in the difficult or inaccessible cases than in the easily accessible cases (15.8±8.7 vs. 9.4±5.0 min, P<0.001). Radial probe positioning within the target lesion during the procedure was more common in the easily accessible cases than in the difficult or inaccessible cases (89.0% vs. 45.9%, P<0.001). Diagnostic results from pathology were obtained for 136/182 (74.7%) easily accessible cases and only 10/37 (27.0%) difficult or inaccessible cases (P<0.001).
Risk factors associated with difficult or inaccessible PPNs
We performed univariable logistic regression analyses to identify clinical variables (Table 2) associated with difficult or inaccessible PPNs. Larger lesion size [per 1 mm, odds ratio (OR) =0.931, P=0.007] and centrality (1 mm farther from pleura, OR =0.963, P=0.046) are likely to increase target lesion accessibility. For bronchus clustering, lower accessibility was observed if the targeted lesion was located at RB1, RB6, LB1, or LB6 (OR =2.288, P=0.044). When the bronchus sign on preprocedural CT was ‘invisible’ (OR =5.741, P=0.004) or the procedure time was longer (per 1 min, OR =1.156, P<0.001), bronchoscopic access to the target lesion was more likely to be difficult.
Table 2
| Variables | OR (95% CI) | P |
|---|---|---|
| Age | 0.995 (0.961–1.031) | 0.79 |
| Male (vs. female) | 1.219 (0.596–2.494) | 0.58 |
| Lesion size | 0.931 (0.883–0.981) | 0.007 |
| Lesion depth from pleural surface | 0.963 (0.927–0.999) | 0.046 |
| Solidity | ||
| Sub-solid (vs. solid) | 0.978 (0.041–2.336) | 0.96 |
| Cavitary (vs. solid) | 0.753 (0.208–2.733) | 0.66 |
| Number of branch levels | 1.040 (0.853–1.269) | 0.69 |
| Lesion location | ||
| Left lower (vs. left upper) | 0.843 (0.266–2.667) | 0.77 |
| Right upper and middle (vs. left upper) | 1.113 (0.430–2.884) | 0.82 |
| Right lower (vs. left upper) | 2.023 (0.792–5.166) | 0.14 |
| Bronchus clustering | ||
| RB8–10, LB8–10 (vs. others) | 0.905 (0.338–2.418) | 0.84 |
| RB1, LB1, RB6, LB6 (vs. others) | 2.288 (1.024–5.116) | 0.044 |
| Bronchus sign on CT before the procedure | ||
| Adjacent (vs. within) | 1.093 (0.348–3.435) | 0.87 |
| Invisible (vs. within) | 5.741 (1.724–19.120) | 0.004 |
| Procedure time | 1.156 (1.090–1.226) | <0.001 |
CI, confidence interval; CT, computed tomography; LB1, apical segmental bronchus of the left upper lobar bronchus; LB6, superior segmental bronchus of the left lower lobar bronchus; LB8–10, anterior, lateral, and posterior segmental bronchi of the left lower lobar bronchus; OR, odds ratio; PPNs, peripheral pulmonary nodules; RB1, apical segmental bronchus of the right upper lobar bronchus; RB6, superior segmental bronchus of the right lower lobar bronchus; RB8–10, anterior, lateral, and posterior segmental bronchi of the right lower lobar bronchus.
Several airway geometrical features were significantly associated with difficult or inaccessible PPNs (Table S3). Among bifurcation angle metrics at the total branch level, higher sum (OR =1.505, P=0.02), maximum (OR =1.595, P=0.01), and average (OR =1.636, P=0.006) values were linked to reduced accessibility. At the branch preceding the final segment, a higher minimum of the min-max diameter ratio indicated greater accessibility (OR =0.641, P=0.02). At the final branch level, increased values of the average max inscribed sphere R (OR =0.623, P=0.01), maximum hydraulic luminal diameter (OR =0.349, P=0.040), and average hydraulic luminal diameter (OR =0.496, P=0.008) were associated with better accessibility. Conversely, greater maximum and average curvature were linked to inaccessibility (OR =1.512, P=0.01; OR =1.474, P=0.03, respectively). Representative cases illustrating these associations are shown in Figure 3.
Performance of the prediction models
To identify patients with PPNs that were difficult to access with rEBUS, we applied the logistic-LASSO method to develop a prediction model and assessed the performance metrics using 10-fold cross-validation (Table 3). In the final composite model, smaller lesion size; closeness to the pleural surface; location in RB1, RB6, LB1, or LB6; invisible bronchus sign; larger bifurcation angle-in; smaller min-max diameter ratio; larger curvature; and smaller hydraulic luminal diameter all predicted difficult PPN access for rEBUS.
Table 3
| Proposed model | Variables | OR (95% CI) |
|---|---|---|
| Clinical model | ||
| Clinical variables | Lesion size | 0.280 (0.231–0.340) |
| Lesion depth from pleural surface | 0.467 (0.439–0.498) | |
| Bronchus clustering (RB1, RB6, LB1, LB6) | 1.372 (1.294–1.454) | |
| Bronchus sign (invisible) | 1.347 (1.308–1.387) | |
| Airway geometrical model | ||
| Total branches | Sum of bifurcation angle-in | 1.195 (1.118–1.277) |
| Max of bifurcation angle-in | 1.621 (1.490–1.763) | |
| Average of bifurcation angle-in | 1.073 (1.022–1.127) | |
| Previous branch to final branch | Min of min–max diameter ratio | 0.634 (0.591–0.680) |
| Final branch | Average of max inscribed sphere R | 0.789 (0.732–0.850) |
| Max curvature | 1.585 (1.474–1.704) | |
| Max hydraulic luminal diameter | 0.528 (0.486–0.575) | |
| Composite model | ||
| Clinical variables | Lesion size | 0.363 (0.299–0.440) |
| Lesion depth from pleural surface | 0.471 (0.415–0.533) | |
| Bronchus clustering (RB1, RB6, LB1, LB6) | 1.278 (1.183–1.381) | |
| Bronchus sign (invisible) | 1.355 (1.288–1.425) | |
| Total branches | Sum of bifurcation angle-in | 1.214 (1.131–1.302) |
| Max of bifurcation angle-in | 1.436 (1.313–1.570) | |
| Previous branch to final branch | Min of min–max diameter ratio | 0.563 (0.518–0.611) |
| Final branch | Max curvature | 1.377 (1.294–1.464) |
| Average curvature | 1.105 (1.044–1.169) | |
| Max hydraulic luminal diameter | 0.481 (0.417–0.554) | |
CI, confidence interval; LASSO, least absolute shrinkage and selection operator; LB1, apical segmental bronchus of the left upper lobar bronchus; LB6, superior segmental bronchus of the left lower lobar bronchus; OR, odds ratio; PPNs, peripheral pulmonary nodules; RB1, apical segmental bronchus of the right upper lobar bronchus; RB6, superior segmental bronchus of the right lower lobar bronchus.
The composite model, which combines clinical variables with airway geometrical variables, performed better than either model alone [AUC (95% CI); composite model =0.84 (0.76–0.92), clinical model =0.68 (0.60–0.81), airway geometrical model =0.79 (0.74–0.84)] (Table 4, Figure 4). The Delong test showed significant model differences (P=0.05). Calibration curves showed good agreement between predicted and actual accessibility, indicating well-developed models (Figure S1). The composite model showed greater consistency between observed and predicted probabilities.
Table 4
| Models | Performance and statistics | |||||
|---|---|---|---|---|---|---|
| AUC | Sensitivity | Specificity | Delong test | |||
| Clinical | Airway geometrical | |||||
| Clinical | 0.68 (0.60–0.81) | 0.15 (0.02–0.28) | 0.99 (0.98–1.00) | Reference | ||
| Airway geometrical | 0.79 (0.74–0.84) | 0.36 (0.22–0.50) | 0.99 (0.97–1.00) | 0.43 | Reference | |
| Composite | 0.84 (0.76–0.92) | 0.51 (0.28–0.74) | 0.98 (0.96–1.00) | 0.05 | 0.28 | |
AUC, sensitivity, and specificity are represented as mean (95% confidence interval). Model performance was evaluated using 10-fold cross validation. The Delong test was performed to compare the model performance (e.g., AUC). AUC, area under the curve.
Discussion
In this study, we prospectively evaluated the accessibility of 219 PPNs (median size, 22.4 mm; median distance from the pleura, 6.8 mm) to rEBUS. To create a model for predicting accessibility, we analyzed both clinical and airway geometrical factors. The composite model showed superior predictive performance (AUC =0.84) compared with the clinical model (0.68). These results suggest that objective quantitative assessment using CT-based airway geometrical data may be more useful than clinical information alone.
Several meta-analyses have reported pooled diagnostic yields of approximately 70% for PPNs measuring 15–40 mm when using rEBUS (17-19). In our study, the median lesion size was 22.4 mm and the diagnostic yield was 66.7%, suggesting comparable lesion characteristics. Unlike previous studies focusing on diagnostic yield, we analyzed the accessibility of PPNs to rEBUS. Because diagnostic yield varies substantially depending on malignancy status and biopsy technique (20), we focused on anatomical accessibility using CT-based airway information. Clinical factor analysis showed that larger lesion size, farther distance from the pleura, and presence of the bronchus sign were associated with easier access—findings consistent with previous meta-analyses (18,21,22). Moreover, we assessed accessibility by bronchus clustering rather than lobe-wise locations and found that apical segments of the upper lobes (RB1, LB1) and posterior segments of the lower lobes (RB6, LB6) were significantly less accessible, likely due to their greater curvature (5).
In the composite model, six airway geometrical features were selected as significant risk factors for difficult or inaccessible cases (Table S4). These included larger bifurcation angles, lower min-max diameter ratio, greater curvature, and smaller hydraulic luminal diameter. Clinically, these features indicate that more acutely angled, irregularly shaped, and narrowed distal airways present barriers to successful rEBUS navigation—findings that are consistent with our previous pilot study (11). The present results not only validate these risk factors in a prospective cohort but also extend prior work by integrating both clinical and quantitative airway data within a unified predictive model.
Unlike earlier studies, our analysis directly compared the predictive value of clinical variables, airway geometrical features, and their combination. The composite model achieved a higher AUC (0.84) than either the clinical (0.68) or airway geometrical (0.79) models alone, with notably improved sensitivity. These findings support the clinical utility of combining quantitative CT-based airway measurements with conventional factors to better predict bronchoscopic accessibility.
Notably, 73.0% of the lesions classified as difficult or inaccessible exhibited a ‘within’ bronchus sign on preprocedural CT, highlighting the limitations of conventional visual assessment. The bronchus sign can often fail to capture underlying anatomical complexities such as sharp bifurcation angles, severe curvature, or airway narrowing. In contrast, our quantitative CT-based airway analysis could provide a more objective and comprehensive assessment, enabling early identification of difficult cases—even when CT appearances seemed favorable. This approach has the potential to enhance preprocedural planning and procedural success, supporting its integration into real-world bronchoscopic workflows.
The proposed CT-based airway geometric analysis is not intended to replace existing bronchoscopic navigation systems, but rather to complement them by providing objective preprocedural information. Unlike many navigation technologies that are applied primarily during bronchoscopy, our approach uses quantitative CT-derived airway measurements to estimate procedural difficulty before the intervention, thereby supporting earlier and more structured decision-making. Because this framework relies on reproducible geometric parameters rather than visual bronchial assessment alone, it may reduce subjectivity in procedural planning. Furthermore, by identifying lesions likely to be difficult to access in advance, this method may help clinicians selectively allocate adjunctive navigation or robotic systems to cases in which they are most likely to provide incremental benefit, thereby improving procedural efficiency and potential cost-effectiveness.
The procedures in this study were performed under conscious sedation without the use of a guide sheath, fluoroscopy, or navigation systems. While these adjuncts have been shown to enhance diagnostic yield in several studies (5,23-25), their routine use across all cases may be limited by medical costs and resource availability (26). Recent meta-analyses of high-cost robotic-assisted bronchoscopy (RAB) reported diagnostic yields approaching 80% (27,28). However, as noted in previous literature, RAB studies often involve a higher proportion of challenging targets, including smaller nodules, bronchus sign-negative lesions, and subsolid opacities, which are inherently more difficult to sample than the lesions typically included in standalone rEBUS studies. Therefore, rather than a direct comparison of diagnostic yields, we propose that our CT-driven airway geometrical analysis (AirGeo) could serve as a tool for procedural triage. By preoperatively identifying lesions predicted to be “inaccessible” via standalone rEBUS, clinicians could strategically refer these high-difficulty cases for advanced platforms like RAB or navigation-guided bronchoscopy under general anesthesia. Such a risk-stratified approach could optimize both clinical outcomes and cost-effectiveness by matching the complexity of the lesion with the appropriate level of technological intervention.
Our study has several limitations. First, although the study was conducted prospectively, it was performed at a single center, and all bronchoscopic procedures were carried out by one experienced bronchoscopist. This design may limit the generalizability of the findings to other institutions, procedural settings, and operators with different levels of experience. However, the use of a single operator also reduced procedural heterogeneity and enabled a more consistent assessment of the relationship between CT-derived airway geometrical features and bronchoscopic accessibility. Future multicenter studies including bronchoscopists with varying levels of expertise are warranted to confirm the external validity and reproducibility of the proposed approach. Another limitation is that the operational definition of bronchoscopic accessibility—although based on the number of procedural attempts, and therefore more structured than a purely subjective assessment—may still be affected by operator-dependent factors such as procedural technique, persistence, and intraprocedural judgment. Accordingly, this definition should be interpreted as a pragmatic procedural surrogate rather than a fully anatomy-specific measure of lesion accessibility. In this study, the use of a single experienced bronchoscopist helped reduce inter-operator variability; however, residual operator dependence cannot be excluded. Future studies are needed to validate this concept using more standardized and less operator-sensitive criteria. Second, our proposed model was not externally validated due to the limited dataset. However, to minimize overfitting and improve model robustness, we conducted internal validation using cross-validation during logistic-LASSO regression modeling. Third, we evaluated only accessibility and not the diagnostic yield of rEBUS for PPNs. As noted in recent literature, diagnostic yield is heavily influenced by confounding variables such as the type of biopsy instruments used (e.g., forceps vs. cryobiopsy), the prevalence of malignancy, and the pathological nature of the lesion (19,20). By focusing strictly on accessibility—defined as the ability to navigate the probe to the target—we aimed to isolate the impact of airway geometric factors from these variables and provide a purer assessment of navigational predictability. Despite those limitations, our study identified the potential of CT-based quantitative assessment for enhancing bronchoscopic procedures and improving PPN accessibility to rEBUS.
Conclusions
In conclusion, this study prospectively evaluated the bronchoscopic accessibility of PPNs by constructing predictive models that incorporate CT-derived airway geometrical features. By integrating these imaging-derived variables with clinical information, we demonstrated their clinical significance in preprocedural planning. Our findings validated the utility of a quantitative airway analysis as an objective tool for assessing procedural difficulty and support its potential role in real-world bronchoscopic workflows. The prospective study design further strengthens the evidence that CT-based airway geometry holds practical relevance and promising potential for improving bronchoscopic strategy and decision-making in clinical practice.
Acknowledgments
None.
Footnote
Reporting Checklist: The authors have completed the TRIPOD reporting checklist. Available at https://tlcr.amegroups.com/article/view/10.21037/tlcr-2026-1-0133/rc
Data Sharing Statement: Available at https://tlcr.amegroups.com/article/view/10.21037/tlcr-2026-1-0133/dss
Peer Review File: Available at https://tlcr.amegroups.com/article/view/10.21037/tlcr-2026-1-0133/prf
Funding: This work was supported by
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://tlcr.amegroups.com/article/view/10.21037/tlcr-2026-1-0133/coif). The authors have no conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. The study was approved by the Institutional Review Board of Samsung Medical Center (No. 2018-03-021). Informed consent was taken from all the patients (or a statement that it was not required and why).
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- Shen YC, Chen CH, Tu CY. Advances in Diagnostic Bronchoscopy. Diagnostics (Basel) 2021;11:1984. [Crossref] [PubMed]
- Steinfort DP, Khor YH, Manser RL, et al. Radial probe endobronchial ultrasound for the diagnosis of peripheral lung cancer: systematic review and meta-analysis. Eur Respir J 2011;37:902-10. [Crossref] [PubMed]
- Chen A, Chenna P, Loiselle A, et al. Radial probe endobronchial ultrasound for peripheral pulmonary lesions. A 5-year institutional experience. Ann Am Thorac Soc 2014;11:578-82.
- Pritchett MA, Bhadra K, Calcutt M, et al. Virtual or reality: divergence between preprocedural computed tomography scans and lung anatomy during guided bronchoscopy. J Thorac Dis 2020;12:1595-611. Erratum in: J Thorac Dis 2020;12:4593-5.
- Kurimoto N, Miyazawa T, Okimasa S, et al. Endobronchial ultrasonography using a guide sheath increases the ability to diagnose peripheral pulmonary lesions endoscopically. Chest 2004;126:959-65. [Crossref] [PubMed]
- Oki M, Saka H, Imabayashi T, et al. Guide sheath versus non-guide sheath method for endobronchial ultrasound-guided biopsy of peripheral pulmonary lesions: a multicentre randomised trial. Eur Respir J 2022;59:2101678. [Crossref] [PubMed]
- Kirby M, Smith BM, Quantitative CT. Scan Imaging of the Airways for Diagnosis and Management of Lung Disease. Chest 2023;164:1150-8. [Crossref] [PubMed]
- Diaz AA, Rahaghi FN, Ross JC, et al. Understanding the contribution of native tracheobronchial structure to lung function: CT assessment of airway morphology in never smokers. Respir Res 2015;16:23. [Crossref] [PubMed]
- Nakano Y, Muro S, Sakai H, et al. Computed tomographic measurements of airway dimensions and emphysema in smokers. Correlation with lung function. Am J Respir Crit Care Med 2000;162:1102-8.
- Washko GR, Diaz AA, Kim V, et al. Computed tomographic measures of airway morphology in smokers and never-smoking normals. J Appl Physiol 1985;2014:668-73. [Crossref] [PubMed]
- Cho HH, Choe J, Kim J, et al. 3D airway geometry analysis of factors in airway navigation failure for lung nodules. Cancer Imaging 2024;24:84. [Crossref] [PubMed]
- Salahuddin M, Sarkiss M, Sagar AS, et al. Ventilatory Strategy to Prevent Atelectasis During Bronchoscopy Under General Anesthesia: A Multicenter Randomized Controlled Trial (Ventilatory Strategy to Prevent Atelectasis -VESPA- Trial). Chest 2022;162:1393-401. [Crossref] [PubMed]
- Moon SM, Choe J, Jeong BH, et al. Diagnostic Performance of Radial Probe Endobronchial Ultrasound without a Guide-Sheath and the Feasibility of Molecular Analysis. Tuberc Respir Dis (Seoul) 2019;82:319-27. [Crossref] [PubMed]
- Miles J. Tolerance and variance inflation factor. In: Everitt BS, Howell DC. editors. Encyclopedia of Statistics in Behavioral Science. John Wiley & Sons, Ltd.; 2005. Available online: https://onlinelibrary.wiley.com/doi/epdf/10.1002/9781118445112.stat06593
- DeLong ER, DeLong DM, Clarke-Pearson DL. Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics 1988;44:837-45.
- Lagendijk J, Steyerberg EW, Daalderop LA, et al. Validation of a prognostic model for adverse perinatal health outcomes. Sci Rep 2020;10:11243. [Crossref] [PubMed]
- McGuire AL, Myers R, Grant K, et al. The Diagnostic Accuracy and Sensitivity for Malignancy of Radial-Endobronchial Ultrasound and Electromagnetic Navigation Bronchoscopy for Sampling of Peripheral Pulmonary Lesions: Systematic Review and Meta-analysis. J Bronchology Interv Pulmonol 2020;27:106-21. [Crossref] [PubMed]
- Lee J, Song JU. Diagnostic yield of radial probe endobronchial ultrasonography-guided transbronchial biopsy without fluoroscopy in peripheral pulmonary lesions: A systematic review and meta-analysis. Thorac Cancer 2023;14:195-205. [Crossref] [PubMed]
- Sryma PB, Mittal S, Madan NK, et al. Efficacy of Radial Endobronchial Ultrasound (R-EBUS) guided transbronchial cryobiopsy for peripheral pulmonary lesions (PPL...s): A systematic review and meta-analysis. Pulmonology 2023;29:50-64. [Crossref] [PubMed]
- Ali MS, Trick W, Mba BI, et al. Radial endobronchial ultrasound for the diagnosis of peripheral pulmonary lesions: A systematic review and meta-analysis. Respirology 2017;22:443-53. [Crossref] [PubMed]
- Tian S, Li X, Liu J, et al. Radial endobronchial ultrasound - guided bronchoscopy for the diagnosis of peripheral pulmonary lesions: A systematic review and meta-analysis of prospective trials. Heliyon 2024;10:e29446. [Crossref] [PubMed]
- Kim SH, Chung HS, Kim J, et al. Development of the Korean Association for Lung Cancer Clinical Practice Guidelines: Recommendations on Radial Probe Endobronchial Ultrasound for Diagnosing Lung Cancer - An Updated Meta-Analysis. Cancer Res Treat 2024;56:464-83. [Crossref] [PubMed]
- Jiang S, Xie F, Mao X, et al. The value of navigation bronchoscopy in the diagnosis of peripheral pulmonary lesions: A meta-analysis. Thorac Cancer 2020;11:1191-201. [Crossref] [PubMed]
- Tamiya M, Okamoto N, Sasada S, et al. Diagnostic yield of combined bronchoscopy and endobronchial ultrasonography, under LungPoint guidance for small peripheral pulmonary lesions. Respirology 2013;18:834-9. [Crossref] [PubMed]
- Eberhardt R, Anantham D, Ernst A, et al. Multimodality bronchoscopic diagnosis of peripheral lung lesions: a randomized controlled trial. Am J Respir Crit Care Med 2007;176:36-41. [Crossref] [PubMed]
- Vakil E, Jackson N, Sainz-Zuñega PV, et al. Optimizing Diagnostic and Staging Pathways for Suspected Lung Cancer: A Decision Analysis. Chest 2021;160:2304-23. [Crossref] [PubMed]
- Zhang C, Xie F, Li R, et al. Robotic-assisted bronchoscopy for the diagnosis of peripheral pulmonary lesions: A systematic review and meta-analysis. Thorac Cancer 2024;15:505-12. [Crossref] [PubMed]
- Ali MS, Ghori UK, Wayne MT, et al. Diagnostic Performance and Safety Profile of Robotic-assisted Bronchoscopy: A Systematic Review and Meta-Analysis. Ann Am Thorac Soc 2023;20:1801-12. [Crossref] [PubMed]

