Three-dimensional printing in radiation oncology: A systematic review of the literature.

PURPOSE/OBJECTIVES: Three-dimensional (3D) printing is recognized as an effective clinical and educational tool in procedurally intensive specialties. However, it has a nascent role in radiation oncology. The goal of this investigation is to clarify the extent to which 3D printing applications are currently being used in radiation oncology through a systematic review of the literature. MATERIALS/METHODS: A search protocol was defined according to preferred reporting items for systematic reviews and meta-analyses (PRISMA) guidelines. Included articles were evaluated using parameters of interest including: year and country of publication, experimental design, sample size for clinical studies, radiation oncology topic, reported outcomes, and implementation barriers or safety concerns. RESULTS: One hundred and three publications from 2012 to 2019 met inclusion criteria. The most commonly described 3D printing applications included quality assurance phantoms (26%), brachytherapy applicators (20%), bolus (17%), preclinical animal irradiation (10%), compensators (7%), and immobilization devices (5%). Most studies were preclinical feasibility studies (63%), with few clinical investigations such as case reports or series (13%) or cohort studies (11%). The most common applications evaluated within clinical settings included brachytherapy applicators (44%) and bolus (28%). Sample sizes for clinical investigations were small (median 10, range 1-42). A minority of articles described basic or translational research (11%) and workflow or cost evaluation studies (3%). The number of articles increased over time (P < 0.0001). While outcomes were heterogeneous, most studies reported successful implementation of accurate and cost-effective 3D printing methods. CONCLUSIONS: Three-dimensional printing is rapidly growing in radiation oncology and has been implemented effectively in a diverse array of applications. Although the number of 3D printing publications has steadily risen, the majority of current reports are preclinical in nature and the few clinical studies that do exist report on small sample sizes. Further dissemination of ongoing investigations describing the clinical application of developed 3D printing technologies in larger cohorts is warranted.

NPIP: A skew line needle configuration optimization system for HDR brachytherapy.

PURPOSE: In this study, the authors introduce skew line needle configurations for high dose rate (HDR) brachytherapy and needle planning by integer program (NPIP), a computational method for generating these configurations. NPIP generates needle configurations that are specific to the anatomy of the patient, avoid critical structures near the penile bulb and other healthy structures, and avoid needle collisions inside the body. METHODS: NPIP consisted of three major components: a method for generating a set of candidate needles, a needle selection component that chose a candidate needle subset to be inserted, and a dose planner for verifying that the final needle configuration could meet dose objectives. NPIP was used to compute needle configurations for prostate cancer data sets from patients previously treated at our clinic. NPIP took two user-parameters: a number of candidate needles, and needle coverage radius, δ. The candidate needle set consisted of 5000 needles, and a range of δ values was used to compute different needle configurations for each patient. Dose plans were computed for each needle configuration. The number of needles generated and dosimetry were analyzed and compared to the physician implant. RESULTS: NPIP computed at least one needle configuration for every patient that met dose objectives, avoided healthy structures and needle collisions, and used as many or fewer needles than standard practice. These needle configurations corresponded to a narrow range of δ values, which could be used as default values if this system is used in practice. The average end-to-end runtime for this implementation of NPIP was 286 s, but there was a wide variation from case to case. CONCLUSIONS: The authors have shown that NPIP can automatically generate skew line needle configurations with the aforementioned properties, and that given the correct input parameters, NPIP can generate needle configurations which meet dose objectives and use as many or fewer needles than the current HDR brachytherapy workflow. Combined with robot assisted brachytherapy, this system has the potential to reduce side effects associated with treatment. A physical trial should be done to test the implant feasibility of NPIP needle configurations.

Feasibility of accelerated image-guided high-dose-rate interstitial brachytherapy with inverse planning simulated annealing (IPSA-HDRBT) for post-operative treatment of pathologically node-negative squamous cell carcinomas of the oral tongue.

PURPOSE: Inverse planning simulated annealing (IPSA) produces highly conformal dose distributions and quick optimizations for high-dose-rate interstitial brachytherapy (HDRBT). We report our dosimetry and overall outcomes using this approach for the accelerated post-operative treatment of pathologically node-negative squamous cell carcinomas of the oral tongue (OTSCC) with high risk of local recurrence. METHODS: Patients with newly diagnosed pN0 OTSCC treated with partial glossectomy, neck dissection, and post-operative HDRBT alone from 2007 to 2021 were retrospectively reviewed. Patients received 30 Gy in 5 fractions over 2.5 days. Target volume and mandible dosimetry are reported. Actuarial rates of local control, regional control, disease-specific survival, and overall survival were estimated using the Kaplan-Meier method. Toxicity was categorized using the Common Terminology Criteria for Adverse Events (CTCAE) version 5.0. RESULTS: 19 consecutive patients were reviewed. Median follow-up was 3.2 years (IQR 1.4-8.2 years) with a 3-year estimated local control rate of 81%. Target volumes were generally small, as the median volume was 12.66 cc. Median V150% and V200% were 52% and 24%, respectively. D1cc and D2cc to the mandible were 17.31 Gy and 14.42 Gy, respectively. CONCLUSIONS: IPSA-HDRBT is feasible and highly efficient for post-operative treatment of the primary tumor bed in patients with pathologically node-negative squamous cell carcinomas of the oral tongue. Further technical optimization and prospective clinical evaluation in a larger patient cohort are planned.

IPIP: A new approach to inverse planning for HDR brachytherapy by directly optimizing dosimetric indices.

PURPOSE: Many planning methods for high dose rate (HDR) brachytherapy require an iterative approach. A set of computational parameters are hypothesized that will give a dose plan that meets dosimetric criteria. A dose plan is computed using these parameters, and if any dosimetric criteria are not met, the process is iterated until a suitable dose plan is found. In this way, the dose distribution is controlled by abstract parameters. The purpose of this study is to develop a new approach for HDR brachytherapy by directly optimizing the dose distribution based on dosimetric criteria. METHODS: The authors developed inverse planning by integer program (IPIP), an optimization model for computing HDR brachytherapy dose plans and a fast heuristic for it. They used their heuristic to compute dose plans for 20 anonymized prostate cancer image data sets from patients previously treated at their clinic database. Dosimetry was evaluated and compared to dosimetric criteria. RESULTS: Dose plans computed from IPIP satisfied all given dosimetric criteria for the target and healthy tissue after a single iteration. The average target coverage was 95%. The average computation time for IPIP was 30.1 s on an Intel(R) Core 2 Duo CPU 1.67 GHz processor with 3 Gib RAM. CONCLUSIONS: IPIP is an HDR brachytherapy planning system that directly incorporates dosimetric criteria. The authors have demonstrated that IPIP has clinically acceptable performance for the prostate cases and dosimetric criteria used in this study, in both dosimetry and runtime. Further study is required to determine if IPIP performs well for a more general group of patients and dosimetric criteria, including other cancer sites such as GYN.

High-dose-rate interstitial brachytherapy for vaginal endometrial cancer recurrence after prior surgery and radiotherapy.

PURPOSE: Characterize the clinical outcomes of endometrial cancer vaginal recurrences after previous surgery and radiation therapy treated with reirradiation including image-guided interstitial high-dose-rate (HDR) brachytherapy. METHODS AND MATERIALS: A single-institution retrospective study identifying women receiving reirradiation for vaginal recurrence of endometrial cancer between 2004 and 2017. RESULTS: Twenty-three women had vaginal recurrences of endometrial cancer, median 13.7 months (range 3.5-104.9) from initial radiation. All received reirradiation with interstitial HDR brachytherapy, and seven also received external beam radiation. Median reirradiation EQD2_10 was 48 Gy (range 24.0-68.81), and median cumulative EQD2_10 was 106.25 Gy (range 62.26-122.0). Median follow-up after reirradiation was 40.2 months (range 4.5-112.7). At 3 years, overall survival was 56%, cancer-specific survival was 61%, and disease-free survival was 46%. 14 patients experienced disease recurrence; 10 including distant sites, one at a regional node only. Three patients experienced local recurrences, two of whom did not complete the prescribed course of reirradiation. The overall crude local control rate was 87%. Three patients experienced Grade 3 vaginal toxicity. There was no bladder or rectal toxicity with Grade >2. CONCLUSIONS: Reirradiation including interstitial HDR brachytherapy is a promising option for vaginal recurrences of endometrial cancer after prior radiation, with high rate of local control and acceptable toxicity. However, distant failure is common. Further studies are needed to determine cumulative radiation dose limits and the role of systemic therapy in this scenario.

Stereotactic Body Radiation Therapy and High-Dose-Rate Brachytherapy Boost in Combination With Intensity Modulated Radiation Therapy for Localized Prostate Cancer: A Single-Institution Propensity Score Matched Analysis.

PURPOSE: To perform a propensity-score matched analysis comparing stereotactic body radiation therapy (SBRT) boost and high-dose-rate (HDR) boost for localized prostate cancer. METHODS AND MATERIALS: A single-institution retrospective chart review was conducted of men treated with pelvic external beam radiation therapy (EBRT) and SBRT boost (21 Gy and 19 Gy in 2 fractions) to the prostate for prostate cancer. A cohort treated at the same institution with HDR brachytherapy boost (19 Gy in 2 fractions) was compared. Propensity-score (PS) matching and multivariable Cox regression were used for analysis. Outcomes were biochemical recurrence freedom (BCRF) and metastasis freedom (MF). RESULTS: One hundred thirty-one men were treated with SBRT boost and 101 with HDR boost with median follow-up of 73.4 and 186.0 months, respectively. In addition, 68.8% of men had high-risk and 26.0% had unfavorable-intermediate disease, and 94.3% received androgen deprivation therapy. Five- and 10-year unadjusted BCRF was 88.8% and 85.3% for SBRT and 91.8% and 74.6% for HDR boost (log-rank P = .3), and 5- and 10-year unadjusted MF was 91.7% and 84.3% for SBRT and 95.8% and 82.0% for HDR (log-rank P = .8). After adjusting for covariates, there was no statistically significant difference in BCRF (hazard ratio [HR] 0.81; 95% confidence interval [CI], 0.37-1.79; P = .6) or MF (HR 1.07; 95% CI, 0.44-2.57; P = .9) between SBRT and HDR boost. Similarly, after PS matching, there was no statistically significant difference between SBRT and HDR (BCRF: HR 0.66, 0.27-1.62, P = .4; MF: HR 0.84, 0.31-2.26, P = .7). Grade 3+ genitourinary and gastrointestinal toxicity in the SBRT cohort were 4.6% and 1.5%, and 3.0% and 0.0% in the HDR cohorts (P = .4, Fisher exact test). CONCLUSIONS: SBRT boost plus pelvic EBRT for prostate cancer resulted in similar BCRF and MF to HDR boost in this single institution, PS matched retrospective analysis. Toxicity was modest. Prospective evaluation of SBRT boost for the treatment of unfavorable-intermediate and high-risk prostate cancer is warranted.

Phase I study of dose escalation to dominant intraprostatic lesions using high-dose-rate brachytherapy.

PURPOSE: Radiation dose escalation for prostate cancer improves biochemical control but is limited by toxicity. Magnetic resonance spectroscopic imaging (MRSI) can define dominant intraprostatic lesions (DIL). This phase I study evaluated dose escalation to MRSI-defined DIL using high-dose-rate (HDR) brachytherapy. MATERIAL AND METHODS: Enrollment was closed early due to low accrual. Ten patients with prostate cancer (T2a-3b, Gleason 6-9, PSA < 20) underwent pre-treatment MRSI, and eight patients had one to three DIL identified. The eight enrolled patients received external beam radiation therapy to 45 Gy and HDR brachytherapy boost to the prostate of 19 Gy in 2 fractions. MRSI images were registered to planning CT images and DIL dose-escalated up to 150% of prescription dose while maintaining normal tissue constraints. The primary endpoint was genitourinary (GU) toxicity. RESULTS: The median total DIL volume was 1.31 ml (range, 0.67-6.33 ml). Median DIL boost was 130% of prescription dose (range, 110-150%). Median urethra V120 was 0.15 ml (range, 0-0.4 ml) and median rectum V75 was 0.74 ml (range, 0.1-1.0 ml). Three patients had acute grade 2 GU toxicity, and two patients had late grade 2 GU toxicity. No patients had grade 2 or higher gastrointestinal toxicity, and no grade 3 or higher toxicities were noted. There were no biochemical failures with median follow-up of 4.9 years (range, 2-8.5 years). CONCLUSIONS: Dose escalation to MRSI-defined DIL is feasible. Toxicity was low but incompletely assessed due to limited patients' enrollment.

Salvage HDR Brachytherapy: Multiple Hypothesis Testing Versus Machine Learning Analysis.

PURPOSE: Salvage high-dose-rate brachytherapy (sHDRB) is a treatment option for recurrences after prior radiation therapy. However, only approximately 50% of patients benefit, with the majority of second recurrences after salvage brachytherapy occurring distantly. Therefore, identification of characteristics that can help select patients who may benefit most from sHDRB is critical. Machine learning may be used to identify characteristics that predict outcome following sHDRB. We aimed to use machine learning to identify patient characteristics associated with biochemical failure (BF) following prostate sHDRB. METHODS AND MATERIALS: We analyzed data for 52 patients treated with sHDRB for locally recurrent prostate cancer after previous definitive radiation therapy between 1998 and 2009. Following pathologic confirmation of locally recurrent disease without evidence of metastatic disease, 36 Gy in 6 fractions was administered to the prostate and seminal vesicles. BF following sHDRB was defined using the Phoenix definition. Sixteen different clinical risk features were collected, and machine learning analysis was executed to identify subpopulations at higher risk of BF. Decision tree-based algorithms including classification and regression trees, MediBoost, and random forests were constructed. RESULTS: Patients were followed up for a minimum of 5 years after sHDRB. Those with a fraction of positive cores ≥0.35 and a disease-free interval <4.12 years after their initial radiation treatment experienced a higher failure rate after sHDRB of 0.75 versus 0.38 for the rest of the population. CONCLUSIONS: Using machine learning, we have identified that patients with a fraction of positive cores ≥0.35 and a disease-free interval <4.1 years might be associated with a high risk of BF following sHDRB.

Magnetic resonance imaging basics for the prostate brachytherapist.

Magnetic resonance imaging (MRI) is increasingly being used in radiation therapy, and integration of MRI into brachytherapy in particular is becoming more common. We present here a systematic review of the basic physics and technical aspects of incorporating MRI into prostate brachytherapy. Terminology and MRI system components are reviewed along with typical work flows in prostate high-dose-rate and low-dose-rate brachytherapy. In general, the brachytherapy workflow consists of five key components: diagnosis, implantation, treatment planning (scan + plan), implant verification, and delivery. MRI integration is discussed for diagnosis; treatment planning; and MRI-guided brachytherapy implants, in which MRI is used to guide the physical insertion of the brachytherapy applicator or needles. Considerations and challenges for establishing an MRI brachytherapy program are also discussed.

A method for restricting intracatheter dwell time variance in high-dose-rate brachytherapy plan optimization.

PURPOSE: To present the algorithm of a modification to the inverse planning simulated annealing (IPSA) optimization engine that allows for restriction of the intracatheter dwell time variance. METHODS AND MATERIALS: IPSA was modified to allow user control of dwell time variance within each catheter through a single parameter, the dwell time deviation constraint (DTDC). The minimum DTDC value (DTDC = 0) does not impose any restriction on dwell time variance, and the maximum value (DTDC = 1) restricts all dwell times within each catheter to take on the same value. The final optimization penalty function value was evaluated as a function of DTDC. RESULTS: The algorithm proposed fully preserves the inverse planning nature of the IPSA algorithm along with the penalty-based dose optimization workflow. Increasing DTDC creates less variance in dwell time between dwell positions in each catheter and may be used to induce a more smooth change in dwell time with dwell position in each catheter. Nonzero DTDC values always increased the optimization penalty function value. CONCLUSIONS: The DTDC was developed as an extension to IPSA to allow restriction of the difference in dwell time between adjacent dwell positions. This results in less variation between neighboring dwell positions which can be clinically desirable. However, the impact of this restriction needs to be considered for its clinical relevance on a case-by-case basis because considerable degradation in dose-volume histogram metrics can result for large DTDC values.

Clinical applications of custom-made vaginal cylinders constructed using three-dimensional printing technology.

PURPOSE: Three-dimensional (3D) printing technology allows physicians to rapidly create customized devices for patients. We report our initial clinical experience using this technology to create custom applicators for vaginal brachytherapy. MATERIAL AND METHODS: Three brachytherapy patients with unique clinical needs were identified as likely to benefit from a customized vaginal applicator. Patient 1 underwent intracavitary vaginal cuff brachytherapy after hysterectomy and chemotherapy for stage IA papillary serous endometrial cancer using a custom printed 2.75 cm diameter segmented vaginal cylinder with a central channel. Patient 2 underwent interstitial brachytherapy for a vaginal cuff recurrence of endometrial cancer after prior hysterectomy, whole pelvis radiotherapy, and brachytherapy boost. We printed a 2 cm diameter vaginal cylinder with one central and six peripheral catheter channels to fit a narrow vaginal canal. Patient 3 underwent interstitial brachytherapy boost for stage IIIA vulvar cancer with vaginal extension. For more secure applicator fit within a wide vaginal canal, we printed a 3.5 cm diameter solid cylinder with one central tandem channel and ten peripheral catheter channels. The applicators were printed in a biocompatible, sterilizable thermoplastic. RESULTS: Patient 1 received 31.5 Gy to the surface in three fractions over two weeks. Patient 2 received 36 Gy to the CTV in six fractions over two implants one week apart, with interstitial hyperthermia once per implant. Patient 3 received 18 Gy in three fractions over one implant after 45 Gy external beam radiotherapy. Brachytherapy was tolerated well with no grade 3 or higher toxicity and no local recurrences. CONCLUSIONS: We established a workflow to rapidly manufacture and implement customized vaginal applicators that can be sterilized and are made of biocompatible material, resulting in high-quality brachytherapy for patients whose anatomy is not ideally suited for standard, commercially available applicators.

Dosimetric analysis of radiation therapy oncology group 0321: the importance of urethral dose.

PURPOSE: Radiation Therapy Oncology Group 0321 is the first multi-institutional cooperative group high-dose-rate (HDR) prostate brachytherapy trial with complete digital brachytherapy dosimetry data. This is a descriptive report of the data and an analysis of toxicity. METHODS AND MATERIALS: Patients are treated with external beam radiation therapy at 45 Gy and 1 HDR implant with 19 Gy in 2 fractions. Implants are done with transrectal ultrasound guidance, and computed tomography (CT)-compatible nonmetallic catheters. HDR planning is done on ≤3-mm-thick CT slices. The "mean DVH" (dose-volume histogram) of the planning target volume (PTV), implanted volume (IP), and organs at risk are calculated. This includes the mean and standard deviation (SD) of the volume at 10-percentage-point intervals from 10% to 200% of the prescribed dose. The conformal index (COIN), homogeneity index (HI), catheters per implant, and patients per institution are calculated. Multivariate analysis and hazard ratios calculation of all the variables against reported grade ≥2 (G2+) genitourinary (GU) adverse events (Common Terminology Criteria for Adverse Events, version 3) are performed. RESULTS: Dosimetry data are based on 122 eligible patients from 14 institutions. The mean of PTV, IP, catheters per implant, and patients per institution are 54 cc, 63 cc, 19 and 9, respectively. The mean of %V100PTV, V80Bladder, V80Rectum, and V120Urethra were 94%, 0.40 cc, 0.15 cc, and 0.25 cc, respectively. There are too few G2+ gastrointestinal adverse event (GI AE) for correlative analysis; thus, the analysis has been performed on the more common G2+ GU AE. There are positive correlations noted between both acute and late G2+ GU AE and urethral dose at multiple levels. Positive correlations with late AE are seen with PTV and IP at high-dose levels. A negative correlation is seen between HI and acute AE. A higher patient accrual rate is associated with a lower rate of G2+ acute and late AE. CONCLUSIONS: Higher urethral dose, larger high-dose volumes, and lower dose homogeneity are associated with greater toxicities. A mean dose-volume histogram comparison at all dose levels should be used for quality control and future research comparison.

Catheter-based ultrasound hyperthermia with HDR brachytherapy for treatment of locally advanced cancer of the prostate and cervix.

A clinical treatment delivery platform has been developed and is being evaluated in a clinical pilot study for providing 3D controlled hyperthermia with catheter-based ultrasound applicators in conjunction with high dose rate (HDR) brachytherapy. Catheter-based ultrasound applicators are capable of 3D spatial control of heating in both angle and length of the devices, with enhanced radial penetration of heating compared to other hyperthermia technologies. Interstitial and endocavity ultrasound devices have been developed specifically for applying hyperthermia within HDR brachytherapy implants during radiation therapy in the treatment of cervix and prostate. A pilot study of the combination of catheter based ultrasound with HDR brachytherapy for locally advanced prostate and cervical cancer has been initiated, and preliminary results of the performance and heating distributions are reported herein. The treatment delivery platform consists of a 32 channel RF amplifier and a 48 channel thermocouple monitoring system. Controlling software can monitor and regulate frequency and power to each transducer section as required during the procedure. Interstitial applicators consist of multiple transducer sections of 2-4 cm length × 180 deg and 3-4 cm × 360 deg. heating patterns to be inserted in specific placed 13g implant catheters. The endocavity device, designed to be inserted within a 6 mm OD plastic tandem catheter within the cervix, consists of 2-3 transducers × dual 180 or 360 deg sectors. 3D temperature based treatment planning and optimization is dovetailed to the HDR optimization based planning to best configure and position the applicators within the catheters, and to determine optimal base power levels to each transducer section. To date we have treated eight cervix implants and six prostate implants. 100 % of treatments achieved a goal of >60 min duration, with therapeutic temperatures achieved in all cases. Thermal dosimetry within the hyperthermia target volume (HTV) and clinical target volume (CTV) are reported. Catheter-based ultrasound hyperthermia with HDR appears feasible with therapeutic temperature coverage of the target volume within the prostate or cervix while sparing surrounding more sensitive regions. (NIHR01CA122276).