Patients
We retrospectively analyzed the records of 111 patients who underwent thoracoscopic segmentectomy for small subsolid pulmonary nodules at the Second Hospital of Shandong University from September 2019 to September 2021. The patients were all clinically diagnosed with peripheral T1aN0M0 stage IA lung cancer. This study was approved by the Ethics Committee and Medical Administration Division of the Second Hospital of Shandong University. Written informed consent was obtained from each of the enrolled patients. The ClinicalTrials.gov ID is NCT04842578. All methods performed in our study were conducted in accordance with the relevant guidelines and regulations.
Our selection criteria for VATS segmentectomy in the present study were according to the NCCN guideline as follows: (I) peripheral GGO suspicious for malignancy, (II) lesion <2 cm in diameter and meeting at least one of the following criteria: (i) pure AIS histology, (ii) ≥50% ground-glass appearance on CT, and (iii) a long doubling time (≥400 days) confirmed by radiologic surveillance.
Reading of thin-slice CT images
Nonenhanced CT or contrast-enhanced CT images could be used, and the slice thickness could be 0.625 mm, 1.0 mm or 1.25 mm. The nomenclature adopted is an accordance with that in a previous publication[10].
The reading, decision-making and planning methods are summarized as follows.
1. The bronchus, pulmonary artery and pulmonary vein should be identified. For the right upper lobe, the reading direction of CT images was first from caudad to cephalad (Supplementary Figures 1). After identification of the bronchus of the right middle lobe, we scanned the thin-slice CT (1.25 mm) images cephalad. The bronchus ventralis was identified as B3, and it could be divided into B3a (Rm. lateralis) and B3b (Rm. medialis), B3b (Rm. medialis) could be divided into B3bi and B3bii (Figure 1a), as presented in the 3D simulation shown in Figure 1b. Then, scanning cephalad was continued uninterrupted. The bronchus dorsalis was identified as B2, and it could be divided into B2a (Rm. dorsalis) and B2b (Rm. horizontalis) (Figure 1c), as presented in the 3D simulation in Figure 1d. The bronchus apicalis was identified as B1, and it could be divided into B1a (Rm. apicalis proprius) and B1b (Rm. ventralis) (Figure 1e), as presented in the 3D simulation in Figure 1f. In our opinion, the definition of the bronchus should follow the "common trunk" principle for the purposes of unification and simplification, for example, B3bii was not defined as B1b in this case. We followed this principle in all cases. It should be noted that the nomenclature does not concern the resection range, which is based on the location of the nodule and on achieving a sufficient surgical margin. Then, the segmental artery, subsegmental artery and vein are identified. Discrimination of the pulmonary artery and pulmonary vein is a key point that can be difficult during the reading of CT images or reconstruction procedure of 3D simulations. The discrimination method applied was as follows (nonenhanced CT): (1) the most important characteristic of the pulmonary artery is that it follows the bronchus tightly, especially in the relatively peripheral part of the lung, while the vein follows a different course, going between or among the bronchi, especially in the relatively peripheral part of the lung, (2) sometimes, two vessels can be seen beside the bronchus, and continuous thin-slice CT images should be read carefully to find the root of the vessel and recognize the origin (Figure 2a-2d, vessel 2 is the pulmonary vein, and vessel 1 is the pulmonary artery. Consecutive illustrations can be found in Supplementary Figures 2). Alternatively, peripheral part of the vessels on CT images should be read, the vein will gradually leave the bronchus, while the artery will continue to follow the bronchus. In this book[10], veins are defined between segments, however, we cannot identify segments through direct observation during clinical practice, and the bronchus is always used instead of segments to define veins. For example, if a vein runs between B1a and B1b (on CT or 3D simulation), it is recognized as V1a. The branches and the vein that runs between them can be observed on the same CT image (V1+2a is between B1+2a and B3c, Figure 3a-b). The branches and the vein that runs between them can be also observed on the longitudinal axis of the human body (V1+2c is between B1+2b and B1+2c), when the continuous reading direction is caudad, is that the main body of B1+2b can be observed (Figure 3c), followed by the main body of V1+2c (Figure 3d) and then the main body of B1+2c (Figure 3e). The 3D anatomic relationship can also be observed on 3D simulation (Figure 3f), and consecutive illustrations of this case can be found in Supplementary Figures 3.
Figure 1 (Figure 1a-1f): The definition of the bronchus.
Figure 2 (Figure 2a-2d): One method to discriminate the pulmonary artery and vein.
Figure 3 (Figure 3a-3f): The identification of pulmonary vein.
2. Then the location of the nodule and the resection range should be evaluated. Reading thin-slice CT images has advantages. Taking this case as an example, interestingly, there was confusion about the location of the nodule on 3D simulation (Figure 4a), which could be resolved well by reading thin-slice CT images. There are two methods to determine the location: (1) if there are branches of intrasegmental pulmonary structures penetrating the nodule, the nodule is located in the corresponding segment, and (2) the location can be determined according to the intersegmental vein. Here, a tiny branch was observed penetrating the nodule (Figure 4b), this branch was identified as a distal branch of V2b (intrasegmental vein of S2) (Figure 4c), the nodule was confirmed to be located in S2. Regarding the resection range, the distance between the upper pole of the nodule and the V2a branch (intersegmental vein between S1 and S2) was 16 mm (Figure 4d). If the nodule and the intersegmental vein are not present in the same section, the most vertical distance between the nodule and the peripheral branch of the intersegmental vein should be measured by counting the CT layers and multiplying by the slice thickness to evaluate resection margin.
Figure 4 (Figure 4a-4d): One method to determine the location of the nodule and evaluate the resection range.
3. Adjacent relationship of structures related to the target segment should be identified for surgical planning (Supplementary Figures 2). The position of surrounding structures can be easily determined on the same CT image, while that of structures superior and inferior to the lesion can be recognized by reading the CT images continuously. For example, in this case, reading the CT images cephalad, first, V2b and V2c are observed (Figure 5a), it can be determined that V2b runs posterior and toward the head (when the reading direction is cephalad, one vessel can be observed running peripherally, it is toward the head) and that V2c is located anterior to V2b, running toward the lateral chest wall and cephalad (Figure 5b). The next important structure is B2, which also runs posterior and cephalad, and V2a is located anterior to B2 (Figure 5c). As the reading procedure continues, A2a and A2b can be seen (no ascending A2 in this case), with A2a running posterior and cephalad and A2b running toward the lateral chest wall and caudad (when the reading direction is cephalad, one vessel can be observed running toward the pulmonary hilum, it is toward caudad). V2a and B1 are located anterior to A2b (Figure 5d).
Figure 5 Figure 5a-5d: To identify the adjacent relationship of pulmonary structures related to the target segment for surgical planning.
3D CT simulation
At our institution, digital imaging and communications in medicine (DICOM) data of thin-slice (0.625 mm or 1.25 mm) CT images were imported into Mimics 21.0 (developed by Materialise Nv Co., Materialise's interactive medical image control system, Kingdom of Belgium) or the InterOperation Thorax Planning system (Beijing Infervision Technology Co., Ltd.) for the 3D reconstruction of pulmonary structures. Contrast-enhanced CT was used for 3D simulation in the 3D simulation group. A sphere extending 2 cm outside the primary tumor can be used to evaluate the safe resection margin.
Nonenhanced CT can be used to finish the 3D reconstruction in only 15-30 minutes by another doctor in case the surgeon of CT group requires structural verification during surgery, meanwhile, a doctor (not the surgeon) can review the 3D simulation to check whether the surgical decision and plan are consistent with those made according to thin-slice CT.
Surgical technique
VATS segmentectomy was performed in all enrolled cases. Computed tomography- guided localization for lung nodules was performed before surgery for all cases. The target segmental artery and bronchus were isolated and divided according to preoperative plan and intraoperative anatomy. The hilar lymph nodes were then routinely obtained for cryosectioning and examination. The venous branches running into the target segment were divided according to preoperative plan or after the distal stump of the target segmental bronchus was lifted and denuded in the peripheral direction. The inflation-deflation line was indicated after the lung was reinflated with pure oxygen. Electrocautery was used to cut toward the inflation-deflation line, the peripheral parenchyma was divided with a stapler, and intersegmental veins were preserved. The location of the nodule was also taken in full consideration to guarantee enough margin. If the lesion was proven to be malignant, N1 and N2 lymph nodes were sampled or dissected.
Statistical analysis
Statistical analysis and graph drawing were performed with Stata 12.0 (StataCorp LP, College Station, TX, USA). Variables were analyzed as the mean and standard deviation. ANOVA analysis was used for comparison when groups are more than two. Mean values were compared using Student's t test, and frequency distributions were compared using the chi-squared test or Fisher's exact test. P values of <0.05 were regarded as significant.