High-dose-rate brachytherapy for cervical cancer after subtotal hysterectomy using standardized applicators cannot achieve a good absorbed-dose coverage of the target volume in special tumor morphologies and topographies due to the steep dose gradient. The aim of this pictorial essay is to present an individualized cylindrical vaginal applicator with oblique guide holes using 3D modeling and printing technologies used at the China-Japan Union Hospital of Jilin University for cervical cancer patients. We use images to describe the steps of this method. Copyright:
High-dose-rate brachytherapy for cervical cancer after subtotal hysterectomy using standardized applicators cannot achieve a good absorbed-dose coverage of the target volume in special tumor morphologies and topographies due to the steep dose gradient. The aim of this pictorial essay is to present an individualized cylindrical vaginal applicator with oblique guide holes using 3D modeling and printing technologies used at the China-Japan Union Hospital of Jilin University for cervical cancerpatients. We use images to describe the steps of this method. Copyright:
External beam radiotherapy (EBRT), concomitant chemotherapy, and brachytherapy is the modern treatment modality for advanced cervical cancer. Brachytherapy is an integral part of radiation treatment of cervical cancer. Studies have shown a clear dose-volume effect between modern dosimetric parameters (D90 of high-risk clinical target volume [HR-CTV] and intermediate-risk clinical target volume [IR-CTV]) and the probability of achieving local control [1,2]. For cervical cancer brachytherapy, common standardized applicators include tandem and ring applicators, tandem and ovoid applicators, vaginal cylinders, combined intracavitary and interstitial (IC/IS) applicators, and the Martinez universal perineal interstitial template [3,4,5,6]. Each applicator has its own unique structural characteristics and is applicable to a specific situation.For special tumor morphology and topography, brachytherapy cannot achieve a good absorbed-dose coverage of the target volume for using standardized applicators due to a steep dose gradient, thereby not achieving good local control. Many studies have made many attempts in the field of individualized applicators and have achieved good results [7,8,9,10]. The aim of this pictorial essay is to demonstrate the workflow of an individualized cylindrical vaginal applicator with oblique guide holes using 3D modeling and printing technologies.
Material and methods
A 47-year-old woman was diagnosed with stage IIIA cervical cancer after having a subtotal hysterectomy for uterine fibroids 10 years prior. A gynecological examination and magnetic resonance imaging (MRI) showed that the tumor was located at the anterior lip of the cervix, the upper one-third of the anterior wall of the vagina, and the left one-third of the parauterine organ, with a size of 5.6 × 3 × 3.8 cm3. The dose of EBRT was 45 Gy/25 fractions delivered by intensity-modulated radiation therapy. After EBRT, a gynecological examination showed that the tumor involved the upper one-third of the anterior vaginal wall with a 2 cm maximum extension in depth. In consideration of her stenotic vagina and tumor location, we applied the individualized vaginal cylinder applicator with guide holes. The design of the applicator was based on the preplanning steps: computed tomography (CT) simulation and contour (Figure 1), establishment of coordinate system of the applicator (Figure 2), optimized design of needle path (Figures 3-5), and dose distribution optimization (Figure 6). The delineation of the target volumes was contoured according to recommendations from the GYN GEC-ESTRO Working Group [11]. Data on the direction for the needles based on the coordinate system are shown in Table 1. According to the data from the preplanning phase, the 3D modeling of the applicator (Figures 7-10) was implemented in 3ds Max (2015 version, Autodesk Inc., San Rafael, CA, USA). The 3D printing process of the applicator is shown in Figures 11-13. The implantation of the applicator and the needles was performed under the guidance of transrectal ultrasound in every application. We used the following brachytherapy treatment schema: 4 fractions, 7 Gy/fraction; 1 week apart for total treatment time of 57 days (including EBRT).
Fig. 1
I. Preplanning: CT simulation was performed after a vaginal model (a half sphere 3 cm in diameter and a cylinder 6 cm in length and 3 cm in diameter) was implanted into the vagina. The target volumes were delineated on reconstructed CT images with reference MR images before implantation. The red dotted line is the HR-CTV, and the green dotted line is the IR-CTV: A) para-axial plane (according to dummy applicator), B) sagittal plane, C) para-coronal plane (according to dummy applicator), D) 3D view
Fig. 2
I. Preplanning: To design convenient directions of the implanted needles, we established the coordinate system of the applicator as follows: the x-axis is the left and right direction, where the left of the patient is positive; the z-axis is the anterior and posterior direction, where the anterior of the patient is positive; the y-axis is the head and foot direction, and the head of the patient is positive. We designed an ideal implanted needle according to the relationship between the tumor (red coil) and the vaginal cylinder (blue) while considering the needle distribution. Two points determine a straight line (the implant needle). The superior point is called the “tip end”, and the inferior point is called the “base end”
Fig. 3
I. Preplanning: The needle directions parallel to the axis of the applicator: A) the transverse section with labels for the other subgraphs: B) the first, C) the second, D) the third to the fifth, E) the sixth, and F) the seventh needles
Fig. 5
I. Preplanning: We obtained the directional data of the needles, which are shown in Table 1, according to the relationship between the needles and the coordinate system: A) base end of applicator, B) 3.3 cm upper from the base end of applicator
Fig. 6
I. Preplanning: The relationship between the preplanned isodose curves and the target volume: A) para-axial plane (according to applicator), B) sagittal plane, C) para-coronal plane (according to applicator), D) 3D view
Table 1
The coordinate data for two points of each implant needle (cm)
Needle umber
Base end
Tip end
x
y
z
x
y
z
1
–1.0
0.0
0.1
–1.0
6.0
0.1
2
–0.6
0.0
1.0
–0.6
6.0
1.0
3
0.0
0.0
0.3
0.0
7.0
2.7
4
0.0
0.0
–0.3
0.0
7.0
1.8
5
0.0
0.0
–1.0
0.0
7.0
0.1
6
0.6
0.0
1.0
0.6
6.0
1.0
7
1.0
0.0
0.1
1.0
6.0
0.1
8
0.0
0.0
1.0
0.5
3.3
1.5
Fig. 7
II. 3D modeling: Modeling of the main body of the vaginal applicator without guide holes: the Boolean union of a cylinder and a sphere
Fig. 10
II. 3D modeling: A mark (red arrow) was included on the main body of the applicator, which just pointed to the anterior of the patient; the Foley catheter passed through this mark
Fig. 11
III. 3D printing: The modeling data were exported as an “*.stl” file to MakerBot Print (version 3.0, MakerBot Industries, USA). In MakerBot Print, the file was processed and exported as a “*. makerbot” file to the 3D printer
Fig. 13
III. 3D printing: The 3D graphics in the preplanning phase (A) and the cylindrical applicator with oblique needles after 3D printing (B) were compared
The coordinate data for two points of each implant needle (cm)I. Preplanning: CT simulation was performed after a vaginal model (a half sphere 3 cm in diameter and a cylinder 6 cm in length and 3 cm in diameter) was implanted into the vagina. The target volumes were delineated on reconstructed CT images with reference MR images before implantation. The red dotted line is the HR-CTV, and the green dotted line is the IR-CTV: A) para-axial plane (according to dummy applicator), B) sagittal plane, C) para-coronal plane (according to dummy applicator), D) 3D viewI. Preplanning: To design convenient directions of the implanted needles, we established the coordinate system of the applicator as follows: the x-axis is the left and right direction, where the left of the patient is positive; the z-axis is the anterior and posterior direction, where the anterior of the patient is positive; the y-axis is the head and foot direction, and the head of the patient is positive. We designed an ideal implanted needle according to the relationship between the tumor (red coil) and the vaginal cylinder (blue) while considering the needle distribution. Two points determine a straight line (the implant needle). The superior point is called the “tip end”, and the inferior point is called the “base end”I. Preplanning: The needle directions parallel to the axis of the applicator: A) the transverse section with labels for the other subgraphs: B) the first, C) the second, D) the third to the fifth, E) the sixth, and F) the seventh needlesI. Preplanning: To obtain better target coverage, we set an oblique needle, which was both forward and left with respect to the patient: A) para-axial plane (according to needle No. 8), B) para-sagittal plane (according to needle No. 8), C) para-coronal plane (according to needle No. 8), D) 3D viewI. Preplanning: We obtained the directional data of the needles, which are shown in Table 1, according to the relationship between the needles and the coordinate system: A) base end of applicator, B) 3.3 cm upper from the base end of applicatorI. Preplanning: The relationship between the preplanned isodose curves and the target volume: A) para-axial plane (according to applicator), B) sagittal plane, C) para-coronal plane (according to applicator), D) 3D viewII. 3D modeling: Modeling of the main body of the vaginal applicator without guide holes: the Boolean union of a cylinder and a sphereII. 3D modeling: The needles were constructed as cylinders (shown within the main body of the vaginal applicator), and the distribution and direction were adjusted according to the data shown in Table 1II. 3D modeling: A cylindrical applicator with oblique guide holes was obtained using Boolean subtractionII. 3D modeling: A mark (red arrow) was included on the main body of the applicator, which just pointed to the anterior of the patient; the Foley catheter passed through this markIII. 3D printing: The modeling data were exported as an “*.stl” file to MakerBot Print (version 3.0, MakerBot Industries, USA). In MakerBot Print, the file was processed and exported as a “*. makerbot” file to the 3D printerIII. 3D printing: The model was printed using a 3D printer (Replicator +, MakerBot USA) with polylactic acid (PLA) materialIII. 3D printing: The 3D graphics in the preplanning phase (A) and the cylindrical applicator with oblique needles after 3D printing (B) were compared
Results
An example (the first fraction) of the relationship between the isodose curves and the target volumes or organs at risk (OARs) is shown in Figure 14. The DVH parameters for all the treatment plans, including the preplanning phase, are shown in Table 2. Brachytherapy was well tolerated and achieved a satisfactory total equivalent dose in 2 Gy/fraction (EQD2) for HR-CTV. At the last follow-up at 24 months after brachytherapy, there was no evidence of recurrence in the pelvic region and no signs or symptoms of complications (see Figure 15).
Fig. 14
IV. Treatment plan: An example (the first fraction) of the relationship between the isodose curves and the target volumes or OARs: A) para-axial plane (according to applicator), B) sagittal plane, C) para-coronal plane (according to applicator), D) 3D view
Table 2
Dose parameters of target volumes and OARs (cGy)
Parameter
Preplanning
BT1
BT2
BT3
BT4
Total*
Planning aims**
Limits dose**
HR-CTV D98
588
618
660
548
591
7662
> 7500
–
HR-CTV D90
714
724
756
732
732
8684
> 9000
> 8500
HR-CTV D50
1177
1157
1112
1200
1173
12786
–
–
IR-CTV D98
415
415
501
318
387
6338
> 6000
–
Bladder D0.1cc
659
691
638
570
539
8783
–
–
Bladder D2cc
571
531
550
474
443
7530
< 8000
< 9000
Rectum D0.1cc
493
445
544
542
383
7337
–
–
Rectum D2cc
383
361
403
381
294
6232
< 6500
< 7500
Sigmoid D0.1cc
677
633
392
617
557
8130
–
–
Sigmoid D2cc
490
467
292
436
417
6622
< 7000
< 7500
Bowel D2cc
459
480
446
413
471
7049
< 7000
< 7500
Total dose was EQD2, which is the dose using a linear-quadratic (LQ) equation with α/β = 10 for a tumor and α/β = 3 for OARs, including 45 Gy/25 f EBRT.
Objective recommendation from EMBRACE II
Fig. 15
IV. Clinical outcome: MRI showed no evidence of recurrence in the pelvic region at A) 9 months, B) 12 months, C) 16 months, and D) 24 months after brachytherapy
Dose parameters of target volumes and OARs (cGy)Total dose was EQD2, which is the dose using a linear-quadratic (LQ) equation with α/β = 10 for a tumor and α/β = 3 for OARs, including 45 Gy/25 f EBRT.Objective recommendation from EMBRACE IIIV. Treatment plan: An example (the first fraction) of the relationship between the isodose curves and the target volumes or OARs: A) para-axial plane (according to applicator), B) sagittal plane, C) para-coronal plane (according to applicator), D) 3D viewIV. Clinical outcome: MRI showed no evidence of recurrence in the pelvic region at A) 9 months, B) 12 months, C) 16 months, and D) 24 months after brachytherapy
Discussion
In conventional 2D brachytherapy, individualized treatment cannot be achieved. The use of 3D brachytherapy shows a trend towards increased local control and improved overall survival with reduced toxicity compared to 2D brachytherapy [12]. At present, individualized treatment is a recommended strategy for cervical cancer brachytherapy in the modern management of cervical cancer [8,13]. This paper introduced the design and application of an individualized cylindrical vaginal applicator using 3D modeling and printing technologies and achieved good clinical results. This paper introduces only the production process of the individualized vaginal cylindrical applicators. The specific 3D modeling software and 3D printer used in this study are replaceable. As long as other software and hardware have the corresponding functions, they can be used to replace the software and hardware introduced in this paper. In the preplanning phase, the delineation of target volumes and OARs is not directly based on MRI but on CT images cross-referenced with an MRI image. Although this method is not the best strategy, it is an acceptable secondary option to some extent. Readers can follow the workflow of this study and perform the preplanning phase based on MRI, which will bring more optimized results. In this study, the design of an individualized cylindrical vaginal applicator was based on the data acquired in the preplanning phase, which is an effective method, and preplanning is also recommended in all brachytherapy plans [14].Four fractions of brachytherapy achieved the ideal dose distribution. Figure 6 shows that the isodose line of the prescribed dose covered the HR-CTV well and avoided OARs well. From Table 2, it can be seen that the total EQD2 of the target volumes and the OARs all achieved the limit dose, and 5 of the 7 items achieved the planning aim dose. This result fully shows that the individualized cylindrical vaginal applicator using 3D modeling and printing technologies has good availability, practicability and repeatability.For cervical cancer brachytherapy after subtotal hysterectomy, cylindrical vaginal applicators are commonly used applicators. Because of the high dose gradient in brachytherapy, when the depth of vaginal infiltration is far greater than 5 mm, it is impossible to achieve both a high dose coverage of the prescribed dose for the target volumes and a low dose to OARs [6]. The individualized applicator introduced in this paper can be used for deep vaginal invasion, vaginal tumors with asymmetric growth, vaginal recurrences, and vaginal cuff recurrence after subtotal hysterectomy.
Conclusions
When standardized applicators cannot achieve a good dose distribution for cervical cancer after subtotal hysterectomy, one can achieve a higher dose coverage of the target and better sparing of the OARs by using a 3D-printed, individualized cylindrical vaginal applicator.
Authors: Jacob Christian Lindegaard; Mikkel Lænsø Madsen; Anders Traberg; Bjarne Meisner; Søren Kynde Nielsen; Kari Tanderup; Harald Spejlborg; Lars Ulrik Fokdal; Ole Nørrevang Journal: Radiother Oncol Date: 2015-12-29 Impact factor: 6.280
Authors: Christian Kirisits; Stefan Lang; Johannes Dimopoulos; Daniel Berger; Dietmar Georg; Richard Pötter Journal: Int J Radiat Oncol Biol Phys Date: 2006-06-01 Impact factor: 7.038
Authors: Christel N Nomden; Astrid A C de Leeuw; Marinus A Moerland; Judith M Roesink; Robbert J H A Tersteeg; Ina Maria Jürgenliemk-Schulz Journal: Int J Radiat Oncol Biol Phys Date: 2011-06-12 Impact factor: 7.038
Authors: Ericka Wiebe; Harry Easton; Gillian Thomas; Lisa Barbera; Laura D'Alimonte; Ananth Ravi Journal: Brachytherapy Date: 2015-01-24 Impact factor: 2.362
Authors: Christine Haie-Meder; Richard Pötter; Erik Van Limbergen; Edith Briot; Marisol De Brabandere; Johannes Dimopoulos; Isabelle Dumas; Taran Paulsen Hellebust; Christian Kirisits; Stefan Lang; Sabine Muschitz; Juliana Nevinson; An Nulens; Peter Petrow; Natascha Wachter-Gerstner Journal: Radiother Oncol Date: 2005-03 Impact factor: 6.280
Authors: Kris Derks; Jacco L G Steenhuijsen; Hetty A van den Berg; Saskia Houterman; Jeltsje Cnossen; Paul van Haaren; Katrien De Jaeger Journal: J Contemp Brachytherapy Date: 2018-02-22