Purpose
Radioactive iodine-125 (125I) seed implantation is being increasingly used for local treatment of tumors [1]. However, as the puncture procedure largely depends on a physician’s clinical experience and level of expertise, quality of the procedure and treatment effect vary, limiting promotion of radioactive 125I seed implantation for clinical use. Three-dimensional (3D)-printed template technology is a milestone in the advancement of radioactive 125I seed implantation therapy. With 3D-printed template guidance, accurate control of the seed needle has been realized, greatly improving adherence to pre-operative plan [2, 3]. Recently, various studies have used an image navigation system to assist in the treatment of head and neck tumors using seed implantation, achieving good puncture accuracy [4]. This study applied an image navigation system and 3D-printed template to radioactive 125I seed implantation treatment for tumors in different parts of the body to evaluate puncture and dose accuracy, and provide a theoretical basis and data support for rational selection and optimization of implantation protocols.
Material and methods
Baseline clinical data
A total of 27 patients with malignant tumors, who presented in our department between December 2019 and December 2020 were enrolled in this study. They were treated using an image navigation system, 3D-printed templates, and computed tomography (CT) to guide radioactive 125I seed implantation. All patients had complete imaging data, pre-operative/post-operative plans, and intra-operative data. All patients had been discussed in a multidisciplinary tumor board before receiving treatment, and met the criteria for radioactive seed implantation [1]: 1) patients who showed relapse after surgery or external radiotherapy or refused surgery or external radiotherapy with a tumor diameter ≤ 7 cm, 2) clear pathological diagnosis, 3) a suitable puncture path, 4) no bleeding tendency or hypercoagulable state, 5) good general physical condition (Karnofsky performance status score > 70), so that radioactive seed implantation could be tolerated, and 6) expected survival > 3 months. All patients signed an informed consent prior to treatment. The research protocol was approved by an ethics committee.
Information on systems, materials, and equipment
Brachytherapy treatment planning system KLSIRPS-3D, Beijing University of Aeronautics and Astronautics, and Beijing Astro Technology was used. The source data in the planning system were derived from TG 43 and its updated document issued by the American Association of Physicists in Medicine [5, 6]. Philips’ Brilliance Big Bore CT scanner, and IGS-MO Optical Image Navigation system and Xinbo Medical Technology (Figure 1) were applied. 3D templates were produced from 3D photocuring rapid-prototyping machine printing (Shanghai 3D Union Tech, RS6000) with an accuracy of 0.1 mm, and photocurable resins that comply with the European Communities standards were used as printing materials. 125I seeds, type 6711_1985 (HTA Co., Ltd., China) with a half-life of 59.4 days and dose-rate constant of 0.965 cGy/(h•U) were utilized. For radioactive 125I seed implantation, Mick Radio-Nuclear Instruments and Eckert & Ziegler BEBIG devices were used.
Pre-operative plan
3D template-guided radioactive 125I seed implantation in accordance with the previous research and published professional standards [7, 8] was performed. Image navigation system assistance was applied based on the following aspects:
Positioning and pre-operative plan design:
CT, with a slice thickness of 2.5 mm was performed 2 days before the operation. The position (supine/prone/lateral) was selected according to the location of lesion; the vacuum pad was fixed, and the body surface was marked with a pendulum line.
CT data were transmitted to brachytherapy treatment planning system to develop a pre-operative plan design. The tumor in target area (gross tumor volume – GTV) and organs at risk (OARs) within 2 cm around the target area were outlined to set the prescribed dose and seed activity as well as to identify the direction, distribution, and depth of the seed needle and set up 1-3 seed needles as navigation needles. The number of seeds was calculated and spatial distribution of the seeds was simulated. The design of needle path should follow principles such as: staying within a space of 1 cm, trying to keep parallel arrangement at the same layer, protecting and avoiding important structures, including nerves, large blood vessels, and cavity viscera, and trying to avoid any bones on the puncture path. Based on these, the number of needles should be as little as possible. Seed activity was usually 0.4-0.6 mCi, and seed spacing design was within 0.5-1.0 cm (which may vary with seed activity and prescription dose) and at least 0.5 cm away from OARs. Using computer optimization combined with manual adjustment, GTV D90 (minimum dose to 90% of the volume of GTV) was adjusted to meet the prescribed dose setting requirements. At the same time, the maximum dose of OARs should not exceed the prescribed dose, and the dose of OARs should be as low as possible.
Design and production of 3D templates:
Personalized digital modeling of the treatment area template in brachytherapy treatment planning system, and additional information on alignment axis and puncture characteristics were used to set the printing range of template.
A 3D template was printed using a 3D photocuring rapid-prototyping machine and medical photocuring resin material. The 3D templates included body surface information, alignment marks, and puncture information with respect to the patients’ treatment area. The thickness of the template was 2.5 mm, and the length of the guide hole was 13 mm. Therefore, the total thickness of the part with the guide hole was 15.5 mm, which ensured that the needle insertion was accurately directed. The diameter of the seed needles was 1.27 ±0.01 mm, and the inner wall of the template guide hole was 1.4 mm in diameter, which ensured smooth insertion of the needle and no movement in the guide hole.
Patient reset and image navigation system-assisted 3D-printed template alignment and fixation:
The patient was reset with reference to the positioning mark, and the template was aligned as per the positioning mark, template’s coordinates, and contour of the patient’s body surface.
An optical tracer was placed on CT frame and 1-3 seed needles (navigation needles), so the relative positions of CT, navigation needle, and patient were confirmed. CT was performed with the image information transmitted to the image navigation system, and the virtual navigation needle could be seen on the patient’s CT images in the image navigation system (in an ideal state, the position of the virtual navigation needle in the patient’s image should be consistent with the position of actual navigation needle). Then, the navigation needle was inserted under the guidance of virtual navigation needle.
After the insertion of navigation needle was completed, CT was performed to confirm accuracy of the needles’ position, and fine-tuning was conducted if there was an error. After the navigation needles were correctly calibrated, the 3D template’s position was considered accurate.
Puncture and seed implantation:
After the 3D template was accurately positioned, the rest of seed needles were percutaneously punctured to the predetermined depth through a template guide hole. During the puncture process, CT was performed to monitor the needle insertion path and depth, and the needle was fine-tuned, if necessary, to avoid damage to organs at risk.
After all the needles were in place, the seeds were implanted one by one according to pre-operative plan, and CT was performed to confirm seed distribution.
As there could be differences between the actual intra-operative situation and pre-operative plan, intra-operative optimization was carried out, if necessary, and the needle insertion path and seed distribution were adjusted in real-time to ensure that the dose received by GTV met the requirements of pre-operative plan.
Post-operative dose evaluation: CT was performed post-operatively; the image was transmitted to the brachytherapy treatment planning system, the target area was recontoured on a final CT for dose verification, and the actual doses delivered to GTV and organs at risk were evaluated. Technical process is presented in Figure 2. The post-operative dose quality evaluation complied with the quality evaluation standard defined by the BC Cancer Research Center [9]. According to immediate verification of target areas D90 and V100 (volume when GTV receives 100% of the prescribed dose), the evaluation results were characterized as excellent, good, fair, and poor (excellent: V100 ≥ 90%, 125% ≥ D90 ≥ 100% [prescription dose]; good: 90% > V100 > 85%, 100% > D90 > 90%; fair: 85% ≥ V100 ≥ 75%, 90% ≥ D90 ≥ 80%, or D90 > 125%; poor: V100 < 75%, D90 < 80%).
Collection and comparison of treatment parameters
Puncture information: In the brachytherapy treatment planning system, the images registered after placing intra-operative navigation needles were fused with pre-operative plan images, and rigid registration for the bone was performed. In the fusion image, both the pre-planned and actual puncture characteristics were displayed. The angle and depth of the puncture before and during the procedure were compared, and the absolute value of difference was taken as the error value. Meanwhile, straight-line distances between two punctures’ points on the body surface and between two tips’ positions were recorded.
Dosimetry information: In the brachytherapy treatment planning system, every dosimetry parameter from pre-operative and post-operative plans was recorded and compared, including GTV D90, V100, V150 (volume when GTV receives 150% of the prescribed dose), V200 (volume when GTV receives 200% of the prescribed dose), and minimum peripheral dose (MPD, the minimum dose received at the edge of GTV) in the target area. With regard to conformal index (CI) [10], CI = (VT, ref/VT) × (VT, ref/Vref), where VT, VT, ref, and Vref are the volume of target area, the volume of target area receiving the prescribed dose, and the total volume (cm3) contained in the prescribed dose, respectively. The optimal CI was 1, which indicated that the prescribed dose covered the target area, but the received dose outside the target area was lower than the prescribed dose. A larger CI indicated that the volume inside the target area receiving the prescribed dose was larger, whereas the volume outside the target area receiving the prescribed dose was smaller. In terms of external index (EI) [11], EI = (Vref – VT, ref)/VT × 100%. The optimal EI was 0, which suggested that tissues outside the target area received less than the prescribed dose. A higher EI indicated that the volume of the prescribed dose received outside the target area was larger. Regarding homogeneity index (HI) [11], HI = (VT, ref – VT, 1.5ref)/VT, ref × 100%, where VT, 1.5ref is the volume (cm3) of the target area receiving 150% of the prescribed dose. The ideal optimal HI was 100%. A higher HI suggested a more uniform dose distribution in the target area. As the location of lesions was scattered and adjacent organs at risk varied, this study was not designed to compare doses of organs at risk.
Other treatment information: In the brachytherapy treatment planning system, other pre-operative and post-operative treatment parameters were collected and compared, including GTV and numbers of needles and seeds.
Statistical methods
For comparisons between groups, Shapiro-Wilk test was first used to verify whether the data in each group were normally distributed. A p-value > 0.05 indicated that data corresponded to normal distribution. For normally distributed data, a t-test (including paired and independent sample t-tests) was applied for comparison, and t-value was used to describe the test value. For non-normally distributed data, a non-parametric test was adopted for comparison (Wilcoxon test for correlated samples and Mann-Whitney U test for independent samples), and a z-value was used to describe the test value. P-values ≤ 0.05 were considered statistically significant, and statistical software SPSS version 25 (IBM Corp., Armonk, NY, USA) was used for all calculations.
Results
A total of 27 patients were included in the study, and 52 navigation needles were used (mean, 2 ±1 needles). The percentage of patients who relapsed after surgery or external radiotherapy, and who refused surgery or radiotherapy were 81% (22 cases) and 19% (5 cases), respectively. Navigation-assisted 3D templates combined with CT-guided radioactive 125I seed implantation could be successfully completed according to established technical procedures. Table 1 shows the baseline information of patients, target lesions, and pre-operative plans. The pre-operative and intra-operative needle angles ranged from 72.5-111.1° and 73.9-109.9°, respectively, and the absolute value of difference between the two groups was 0.0-1.7°. The pre-operative and intra-operative needle depths ranged from 20.0-110.0 mm and 20.0-118.0 mm, respectively, and the absolute value of difference between the two groups was 0-8.9 mm. Table 2 describes the results of comparison between the two groups. In addition, the range, median, and mean values of the error of needle insertion point for the two groups were 0-3.9 mm, 1.3 mm, and 1.7 ±1.0 mm, respectively, and those of the needle tip were 0-0.9 mm, 2.9 mm, and 3.1 ±1.8 mm, respectively. The range and mean values of the needle depth in the head & neck and chest wall were 20.0-66.0 mm and 34.0 ±11.0 mm, respectively, and those in the retroperitoneal and paravertebral and pelvic parts were 43.0-118.0 mm and 58.0 ±17.0 mm, respectively. The results of Mann-Whitney U test showed that there was a significant difference in the needle depth between the two groups (p < 0.001). If divided into two groups by position, errors of retroperitoneal and paravertebral and pelvic cavity were significantly larger than those of the head and neck and chest wall (Table 3), and differences in the angle, depth, and insertion point were statistically significant (p < 0.05). Table 4 summarizes the comparison of dosimetry parameters between pre-operative and post-operative plans. With the exception of MPD, other parameters were similar in both groups, and there was no significant difference (p > 0.05). The MPD in the post-operative plan was higher than that in the pre-operative plan (mean, 72.1 Gy and 63.8 Gy, respectively; p = 0.029). The actual number of seeds used during the surgery was higher than that planned before the operation (mean, 61.4 and 59.7, respectively, p = 0.031). The volume of target area (GTV) in the post-operative plan was larger than that in the pre-operative plan (mean, 41.3 cm3 and 40.2 cm3, respectively, p = 0.005). In terms of post-operative dose quality evaluation, the results were excellent in 18 cases (67%), good in 6 cases (22%), fair in 2 cases (7%), and poor in 1 patient (4%). Of the 2 patients with fair evaluation quality, one case involved the head and neck, and the other involved the pelvic area, with a post-operative D90 of 128.4% (166.9 Gy/130 Gy) and 87.7% (131.5 Gy/150 Gy) of the prescribed dose, respectively. One case with poor evaluation quality was a pelvic case, and the post-operative D90 was 79.8% of the prescribed dose (103.7 Gy/130 Gy).
Table 1
Table 2
Table 3
Table 4
Discussion
Although brachytherapy treatment planning system enables doctors to develop more rational treatment plans before procedure, if the position of seed needle is not accurate during surgery, it is often difficult for post-operative dose to meet the expectations of pre-operative plan, thus affecting follow-up treatment. The ultimate goal of both 3D templates and image navigation system technologies is to allow more accurate insertion of the seed needle according to pre-operative plan. In terms of 3D templates, a number of related studies have described their good accuracy, with post-operative dose meeting the requirements of pre-operative plan [2, 7, 8]. Image navigation systems are another solution for puncture intervention technology. This technology fuses and registers CT, magnetic resonance, or positron emission tomography images obtained before the procedure, with intra-operative real-time images. During the operation, the position of puncture needle is tracked in real-time by a tracer (optical or magnetic positioning) to locate the lesion on the fusion image, thereby guiding the operator to perform the puncture [12]. Owing to their real-time, dynamic, and visible characteristics, image navigation systems are widely used in biopsy, ablation, and other puncture-related procedures [13-15]. However, the combination of the two techniques for radioactive 125I seed implantation and its application in the treatment of tumors of many parts of the body have rarely been reported. Current data show that the accuracy of image navigation system-assisted 3D template-guided radioactive 125I seed implantation treatment is good.
According to different positioning principles, image navigation systems mainly include magnetic and optical positioning and navigation [16]: 1. Magnetic positioning and navigation mainly fix the sensor coil on a tracked device (such as puncture device). When the sensor coil moves relative to the magnetic field emitter, it produces different intensity currents and then, locates the tracked device through current signal. 2. Optical positioning and navigation are primarily used to fix the tracer on the CT machine and tracked device, and infrared camera is used to directly detect the position of tracer. Compared with magnetic positioning, optical positioning is faster in data transmission and has higher accuracy. The disadvantage is that there can be no obstruction between the optical camera and tracked device, and the position of the needle tip cannot be tracked. Optical positioning and navigation were adopted in the present study, but the optical positioning and navigation tracer was large and occupied a large space after being fixed to the seed needle. Furthermore, radioactive 125I seed implantation is a multi-needle operation, so it is impossible to equip all seed needles with tracers. Therefore, in the current research, the main function of navigation was to confirm whether template resetting was accurate. We selected 1-3 navigation needles, and if the position of the navigation needle was accurate, it would prove that the position of the template was accurate, but also the precision of the needle under template guidance was good and subsequent needle insertions could be carried out without navigation. We wanted to know the actual needle track error in patients and its impact on dosimetry distribution.
Based on this study, the needle accuracy under the image navigation system with 3D templates was good. The mean errors of the angle, depth, insertion point, and tip were 0.5 ±0.5°, 4.0 ±2.0 mm, 1.7 ±1.0 mm, and 3.1 ±1.8 mm, respectively. Therefore, the mean error of the angle was < 1°, and the mean error of the distance was approximately 3 mm. Compared with previous studies on 3D templates and/or image navigation system-assisted brachytherapy (the average tip error was 0.86-7.00 mm, and the average angle error was 1.9-5.6°), the accuracy of our study was similar or even better [4, 17-19]. Good accuracy reduces the number of CT scan confirmations and shortens the overall operation time. For angle comparison, the pre-operative and intra-operative angles were not significantly different (p = 0.271), but the intra-operative depth was slightly greater than the pre-operative depth (44.0 ±18.0 mm and 43.0 ±18.0 mm, respectively, p = 0.002). The reasons for the errors were considered to be changes in skin contour and thickness caused by anesthesia, template alignment issues, and unsatisfactory fit between template and body surface. In terms of the treatment site, angle, depth, and insertion point, the errors in the pelvic/retroperitoneal region were larger than those in the head, neck, and chest wall (p < 0.05). In addition, these errors might be attributed to longer insertion paths that led to greater errors (the difference in insertion depth between the two groups was statistically significant, p < 0.05). Moreover, the stability of tissues and organs in the head, neck, and chest wall treatment area was better than that in the pelvic/retroperitoneal region. This prompted us to focus on the accuracy of needle insertion and perform multiple CT scans to confirm if it was necessary in areas where the treatment location was deep, and the tissue was soft and prone to change.
By comparing the pre-operative and post-operative dosimetry parameters in this study, we found that most were not significantly different, suggesting that the post-operative dose could better meet the requirements of pre-operative plan. The only dosimetry index with a significant difference was MPD, which was alternatively expressed as D100 (when GTV receives 100% of the dose). MPD in the post-operative plan was higher (p < 0.05), indicating that the error in this study did not result in a reduction in dose to the target area. The reason may be that we replanted seeds in the patients, whose immediate post-operative CT scans showed unsatisfactory seed distribution. This also explains why the number of seeds in the post-operative plan was higher than that in the pre-operative plan, and the difference was statistically significant (p = 0.031). Moreover, GTV in the post-operative plan was larger than that in the pre-operative plan (p = 0.005), which was possibly related to intra-operative and post-operative bleeding, edema, and inflammatory reactions. The post-operative D90 of the three cases with fair and poor evaluation quality did not reach 90% and 80% of the prescribed dose, respectively, suggesting that radioactive 125I seed implantation is a surgery-dependent treatment. Because the dose limit of most normal tissues in the treatment of radioactive 125I seed implantation is not clear, except for rectum and urethra [20], OARs’ dose should not exceed the prescribed dose and should be as low as possible in planning design (but there is no specific value to limit), and then the dose data of OARs for analysis could be recorded if necessary. It is believed that the accuracy of treatment would improve, as physicians become more experienced and proficient with the procedure. Many studies have analyzed the dose accuracy of radioactive 125I seed implantation, guided by 3D templates, and results indicate that dose accuracy was good [2, 8, 21, 22]; however, there was no analysis of needle path error in these studies. Due to large volume of navigation tracer and the need for optical tracking, it can only be placed in the body or tail of the seed needle; therefore, it is impossible to navigate all needles. This study is only a preliminary exploration of the feasibility of navigation combined with 3D templates for seed implantation. Until the instruments are improved, and the combination of navigation systems and 3D-templates is further explored, we cannot be confident whether the combined mode is better than the simple application of template.
Based on above analysis, the solutions we propose to reduce errors of implants include the following: 1. While designing the template, areas with significant body surface features should be included for better accuracy and stability of the alignment of 3D template. 2. Moderate kneading after local infiltration anesthesia, which could not only promote the absorption of anesthetics, but also reduce body surface deformation. 3. The shortest needle track should be selected to reduce the influence of deviation in needle insertion. 4. During seed implantation, the force of pull out the needle and push in the seed should be moderate to avoid changes in needle position and seed shifting. 5. The combination of navigator and needle should be improved to realize multi-needle navigation. 6. Immediate post-operative CT allows for real-time observation of seed distribution. For patients with poor seed distribution, timely replanting of seeds could effectively prevent cold spot of the dose in target area [1]. In addition, high-dose-rate (HDR) brachytherapy can often be used as dose compensation for external beam irradiation [23, 24]; similarly, HDR brachytherapy might be considered to compensate implant inaccuracies with post-operative optimization.
The limitations of this study are as follows: 1. The technical development time was slightly short, and the sample size was small. 2. The results were limited to comparisons between treatment plans, and further observation and follow-up are needed. 3. Due to the small number of cases and scattered treatment sites, effective comparison of the dose distribution to different OARs was not available. 4. The combination and process of the image navigation system with 3D templates remain in the initial stages, and need further exploration and improvement. In a follow-up study, we will include more cases and conduct in-depth and detailed research.
Conclusions
The combination of image navigation system and 3D-printed template technologies for radioactive iodine-125 seed implantation treatment showed good accuracy and feasibility, and exhibited good quality regrading completion of implantation plan. Errors in the actual needle insertion during the surgery were smaller than those before the procedure. The actual D90 of target area, CI, and dose uniformity, all met the design requirements of the pre-operative plan. In a follow-up study, we will expand the number of cases, refine the study, and clarify the efficacy and safety of this approach based on clinical data.