Novel Brachytherapy Technology

@article{nokey_29,

title = {Fast DM,M calculation in LDR brachytherapy using deep learning methods},

author = {Francisco Berumen and Shirin A. Enger and Luc Beaulieu},

doi = {10.1088/1361-6560/accd42},

year  = {2023},

date = {2023-05-23},

urldate = {2023-05-23},

journal = {Physics in Medicine & Biology},

abstract = {Objective.The Monte Carlo (MC) method provides a complete solution to the tissue heterogeneity effects in low-energy low-dose rate (LDR) brachytherapy. However, long computation times limit the clinical implementation of MC-based treatment planning solutions. This work aims to apply deep learning (DL) methods, specifically a model trained with MC simulations, to predict accurate dose to medium in medium (DM,M) distributions in LDR prostate brachytherapy.Approach.To train the DL model, 2369 single-seed configurations, corresponding to 44 prostate patient plans, were used. These patients underwent LDR brachytherapy treatments in which125I SelectSeed sources were implanted. For each seed configuration, the patient geometry, the MC dose volume and the single-seed plan volume were used to train a 3D Unet convolutional neural network. Previous knowledge was included in the network as anr2kernel related to the first-order dose dependency in brachytherapy. MC and DL dose distributions were compared through the dose maps, isodose lines, and dose-volume histograms. Features enclosed in the model were visualized.Main results.Model features started from the symmetrical kernel and finalized with an anisotropic representation that considered the patient organs and their interfaces, the source position, and the low- and high-dose regions. For a full prostate patient, small differences were seen below the 20% isodose line. When comparing DL-based and MC-based calculations, the predicted CTVD90metric had an average difference of -0.1%. Average differences for OARs were -1.3%, 0.07%, and 4.9% for the rectumD2cc, the bladderD2cc, and the urethraD0.1cc. The model took 1.8 ms to predict a complete 3DDM,Mvolume (1.18 M voxels).Significance.The proposed DL model stands for a simple and fast engine which includes prior physics knowledge of the problem. Such an engine considers the anisotropy of a brachytherapy source and the patient tissue composition.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{nokey_28,

title = {Applying the column generation method to the intensity modulated high dose rate brachytherapy inverse planning problem},

author = {Majd Antaki and Marc-André Renaud and Marc Morcos and Jan Seuntjens and Shirin A. Enger },

doi = {10.1088/1361-6560/acbc63},

year  = {2023},

date = {2023-03-13},

urldate = {2023-03-13},

journal = {Physics in Medicine & Biology},

abstract = {Objective.Intensity modulated high dose rate brachytherapy (IMBT) is a rapidly developing application of brachytherapy where anisotropic dose distributions can be produced at each source dwell position. This technique is made possible by placing rotating metallic shields inside brachytherapy needles or catheters. By dynamically directing the radiation towards the tumours and away from the healthy tissues, a more conformal dose distribution can be obtained. The resulting treatment planning involves optimizing dwell position and shield angle (DPSA). The aim of this study was to investigate the column generation method for IMBT treatment plan optimization.Approach.A column generation optimization algorithm was developed to optimize the dwell times and shield angles. A retrospective study was performed on 10 prostate cases using RapidBrachyMCTPS. At every iteration, the plan was optimized with the chosen DPSA which would best improve the cost function that was added to the plan. The optimization process was stopped when the remaining DPSAs would not add value to the plan to limit the plan complexity.Main results.The average number of DPSAs and voxels were 2270 and 7997, respectively. The column generation approach yielded near-optimal treatment plans by using only 11% of available DPSAs on average in ten prostate cases. The coverage and organs at risk constraints passed in all ten cases.Significance.The column generation method produced high-quality deliverable prostate IMBT plans. The treatment plan quality reached a plateau, where adding more DPSAs had a minimal effect on dose volume histogram parameters. The iterative nature of the column generation method allows early termination of the treatment plan creation process as soon as the dosimetric indices from dose volume histogram satisfy the clinical requirements or if their values stabilize.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{nokey_30,

title = {Dosimetric impact of a robust optimization approach to mitigate effects from rotational uncertainty in prostate intensity-modulated brachytherapy},

author = {Björn Morén and Majd Antaki and Gabriel Famulari and Marc Morcos and Torbjörn Larsson and Shirin A Enger and Åsa Carlsson Tedgren},

doi = {10.1002/mp.16134},

year  = {2022},

date = {2022-12-22},

journal = {Medical Physics},

abstract = {Background: Intensity-modulated brachytherapy (IMBT) is an emerging technology for cancer treatment, in which radiation sources are shielded to shape the dose distribution. The rotatable shields provide an additional degree of freedom, but also introduce an additional, directional, type of uncertainty, compared to conventional high-dose-rate brachytherapy (HDR BT).



Purpose: We propose and evaluate a robust optimization approach to mitigate the effects of rotational uncertainty in the shields with respect to planning criteria.



Methods: A previously suggested prototype for platinum-shielded prostate 169 Yb-based dynamic IMBT is considered. We study a retrospective patient data set (anatomical contours and catheter placement) from two clinics, consisting of six patients that had previously undergone conventional 192 Ir HDR BT treatment. The Monte Carlo-based treatment planning software RapidBrachyMCTPS is used for dose calculations. In our computational experiments, we investigate systematic rotational shield errors of ±10° and ±20°, and the same systematic error is applied to all dwell positions in each scenario. This gives us three scenarios, one nominal and two with errors. The robust optimization approach finds a compromise between the average and worst-case scenario outcomes.



Results: We compare dose plans obtained from standard models and their robust counterparts. With dwell times obtained from a linear penalty model (LPM), for 10° errors, the dose to urethra (D0.1cc) and rectum (D0.1cc and D1cc) increase with up to 5% and 7%, respectively, in the worst-case scenario, while with the robust counterpart, the corresponding increases were 3% and 3%. For all patients and all evaluated criteria, the worst-case scenario outcome with the robust approach had lower deviation compared to the standard model, without compromising target coverage. We also evaluated shield errors up to 20° and while the deviations increased to a large extent with the standard models, the robust models were capable of handling even such large errors.



Conclusions: We conclude that robust optimization can be used to mitigate the effects from rotational uncertainty and to ensure the treatment plan quality of IMBT.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{weishaupt2022po,

title = {PO-1325 Automated rectal tumor segmentation with inter-observer variability-based uncertainty estimates},

author = {Luca L. Weishaupt and Te Vuong and Alana Thibodeau-Antonacci and A Garant and K Singh and C Miller and A Martin and F Schmitt-Ulms and Shirin A. Enger},

year  = {2022},

date = {2022-01-01},

urldate = {2022-01-01},

journal = {Radiotherapy and Oncology},

volume = {170},

pages = {S1120--S1121},

publisher = {Elsevier},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{weishaupt2022a121,

title = {A121 QUANTIFYING INTER-OBSERVER VARIABILITY IN THE SEGMENTATION OF RECTAL TUMORS IN ENDOSCOPY IMAGES AND ITS EFFECTS ON DEEP LEARNING},

author = {Luca L Weishaupt and Te Vuong and Alana Thibodeau-Antonacci and A Garant and KS Singh and C Miller and A Martin and Shirin A. Enger},

year  = {2022},

date = {2022-01-01},

urldate = {2022-01-01},

journal = {Journal of the Canadian Association of Gastroenterology},

volume = {5},

number = {Supplement_1},

pages = {140--142},

publisher = {Oxford University Press US},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{lecavalier-barsoum_utilization_2021,

title = {Utilization of brachytherapy in Quebec, Canada},

author = {Magali Lecavalier-Barsoum and Farzin Khosrow-Khavar and Krum Asiev and Marija Popovic and Te Vuong and Shirin A. Enger},

doi = {10.1016/j.brachy.2021.07.002},

issn = {1873-1449},

year  = {2021},

date = {2021-08-01},

journal = {Brachytherapy},

pages = {S1538--4721(21)00452--9},

abstract = {BACKGROUND AND PURPOSE: Despite the excellent clinical outcomes from brachytherapy treatments compared with other modalities and the low associated costs, there have been reports of a decline in utilization of brachytherapy. The aim of this study was to investigate in detail the trend in utilization of brachytherapy in the province of Québec, Canada, from 2011 to 2019. 

MATERIALS AND METHODS: All radiotherapy clinics in the province of Quebec, and among these the clinics that provide brachytherapy treatments, were identified. This observational retrospective cohort study involved analysis of data compiled by the Ministère de la Santé et des Services Sociaux du Québec for the period of 2011 to end of 2019 on all brachytherapy procedures performed in the province of Quebec. Time series graphs were used to describe the number of high dose rate (HDR) and low dose rate (LDR) brachytherapy treatments during the studied time period. Statistical analysis was conducted using R statistical software. 

RESULTS: Between 2011 and 2019, 12 hospitals in the province of Québec provided radiotherapy treatments, and all of them offered brachytherapy services. The median annual number of brachytherapy sessions was 4413 (range 3930-4829). HDR brachytherapy represented over 90% of all brachytherapy treatments throughout the study period. Significant changes over time were observed in the number of treatments: at least 5% change was seen only for the two most common subtypes of brachytherapy, HDR interstitial and HDR intracavitary, with an increase of 9.6% and a decrease of 9.2%, respectively. The use of other subtypes of brachytherapy (HDR-plesiotherapy, LDR-interstitial, LDR-intracavitary, LDR-eye plaque) was stable between 2011 and 2019, with ≤ 2.5% variation. 

CONCLUSION: This study demonstrates an overall steady use of brachytherapy between 2011 and 2019 in Quebec. Brachytherapy offers numerous advantages for the treatment of diverse cancer sites. Although more sophisticated external beam radiotherapy treatments have emerged in the last decades, the precision and cost-effectiveness of brachytherapy remain unbeaten. To ensure the continued use and availability of brachytherapy, governments must put in place policies and regulations to that effect. Training and exposure of future health care professionals to brachytherapy within Quebec and Canada is essential to provide all patients the same access to this life saving modality.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{morcos_impact_2021,

title = {On the impact of absorbed dose specification, tissue heterogeneities, and applicator heterogeneities on Monte Carlo-based dosimetry of Ir-192, Se-75, and Yb-169 in conventional and intensity-modulated brachytherapy for the treatment of cervical cancer},

author = {Marc Morcos and Akila N. Viswanathan and Shirin A. Enger},

doi = {10.1002/mp.14802},

issn = {2473-4209},

year  = {2021},

date = {2021-05-01},

journal = {Medical Physics},

volume = {48},

number = {5},

pages = {2604--2613},

abstract = {PURPOSE: The purpose of this study was to evaluate the impact of dose reporting schemes and tissue/applicator heterogeneities for 192 Ir-, 75 Se-, and 169 Yb-based MRI-guided conventional and intensity-modulated brachytherapy. METHODS AND MATERIALS: Treatment plans using a variety of dose reporting and tissue/applicator segmentation schemes were generated for a cohort (n = 10) of cervical cancer patients treated with 192 Ir-based Venezia brachytherapy. Dose calculations were performed using RapidBrachyMCTPS, a Geant4-based research Monte Carlo treatment planning system. Ultimately, five dose calculation scenarios were evaluated: (a) dose to water in water (Dw,w ); (b) Dw,w taking the applicator material into consideration (Dw,wApp ); (c) dose to water in medium (Dw,m ); (d and e) dose to medium in medium with mass densities assigned either nominally per structure (Dm,m (Nom) ) or voxel-by-voxel (Dm,m ). 

RESULTS: Ignoring the plastic Venezia applicator (Dw,wApp ) overestimates Dm,m by up to 1% (average) with high energy source (192 Ir and 75 Se) and up to 2% with 169 Yb. Scoring dose to water (Dw,wApp or Dw,m ) generally overestimates dose and this effect increases with decreasing photon energy. Reporting dose other than Dm,m (or Dm,m Nom ) for 169 Yb-based conventional and intensity-modulated brachytherapy leads to a simultaneous overestimation (up to 4%) of CTVHR D90 and underestimation (up to 2%) of bladder D2cc due to a significant dip in the mass-energy absorption ratios at the depths of nearby targets and OARs. Using a nominal mass-density assignment per structure, rather than a CT-derived voxel-by-voxel assignment for MRI-guided brachytherapy, amounts to a dose error up to 1% for all radionuclides considered. 

CONCLUSIONS: The effects of the considered dose reporting schemes trend correspondingly between conventional and intensity-modulated brachytherapy. In the absence of CT-derived mass densities, MRI-only-based dosimetry can adequately approximate Dm,m by assigning nominal mass densities to structures. Tissue and applicator heterogeneities do not significantly impact dosimetry for 192 Ir and 75 Se, but do for 169 Yb; dose reporting must be explicitly defined since Dw,m and Dw,w may overstate the dosimetric benefits.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{morcos_novel_2021,

title = {A novel minimally invasive dynamic-shield, intensity-modulated brachytherapy system for the treatment of cervical cancer},

author = {Marc Morcos and Majd Antaki and Akila N. Viswanathan and Shirin A. Enger},

doi = {10.1002/mp.14459},

issn = {2473-4209},

year  = {2021},

date = {2021-01-01},

journal = {Medical Physics},

volume = {48},

number = {1},

pages = {71--79},

abstract = {PURPOSE: To present a novel, MRI-compatible dynamicshield intensity modulated brachytherapy (IMBT) applicator and delivery system using 192 Ir, 75 Se, and 169 Yb radioisotopes for the treatment of locally advanced cervical cancer. Needle-free IMBT is a promising technique for improving target coverage and organs at risk (OAR) sparing. 

METHODS AND MATERIALS: The IMBT delivery system dynamically controls the rotation of a novel tungsten shield placed inside an MRI-compatible, 6-mm wide intrauterine tandem. Using 36 cervical cancer cases, conventional intracavitary brachytherapy (IC-BT) and intracavitary/interstitial brachytherapy (IC/IS-BT) (10Ci 192 Ir) plans were compared to IMBT (10Ci 192 Ir; 11.5Ci 75 Se; 44Ci 169 Yb). All plans were generated using the Geant4-based Monte Carlo dose calculation engine, RapidBrachyMC. Treatment plans were optimized then normalized to the same high-risk clinical target volume (HR-CTV) D90 and the D2cc for bladder, rectum, and sigmoid in the research brachytherapy planning system, RapidBrachyMCTPS. Plans were renormalized until either of the three OAR reached dose limits to calculate the maximum achievable HR-CTV D90 and D98 . RESULTS: Compared to IC-BT, IMBT with either of the three radionuclides significantly improves the HR-CTV D90 and D98 by up to 5.2% ± 0.3% (P textless 0.001) and 6.7% ± 0.5% (P textless 0.001), respectively, with the largest dosimetric enhancement when using 169 Yb followed by 75 Se and then 192 Ir. Similarly, D2cc for all OAR improved with IMBT by up to 7.7% ± 0.6% (P textless 0.001). For IC/IS-BT cases, needle-free IMBT achieved clinically acceptable plans with 169 Yb-based IMBT further improving HR-CTV D98 by 1.5% ± 0.2% (P = 0.034) and decreasing sigmoid D2cc by 1.9% ± 0.4% (P = 0.048). Delivery times for IMBT are increased by a factor of 1.7, 3.3, and 2.3 for 192 Ir, 75 Se, and 169 Yb, respectively, relative to conventional 192 Ir BT. 

CONCLUSIONS: Dynamic shield IMBT provides a promising alternative to conventional IC- and IC/IS-BT techniques with significant dosimetric enhancements and even greater improvements with intermediate energy radionuclides. The ability to deliver a highly conformal, OAR-sparing dose without IS needles provides a simplified method for improving the therapeutic ratio less invasively and in a less resource intensive manner.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{morcos_impact_2021b,

title = {On the impact of absorbed dose specification, tissue heterogeneities, and applicator heterogeneities on Monte Carlo-based dosimetry of Ir-192, Se-75, and Yb-169 in conventional and intensity-modulated brachytherapy for the treatment of cervical cancer},

author = {Marc Morcos and Akila N. Viswanathan and Shirin A. Enger},

url = {https://aapm.onlinelibrary.wiley.com/doi/abs/10.1002/mp.14802},

doi = {10.1002/mp.14802},

issn = {2473-4209},

year  = {2021},

date = {2021-01-01},

urldate = {2021-09-08},

journal = {Medical Physics},

volume = {48},

number = {5},

pages = {2604--2613},

abstract = {Purpose The purpose of this study was to evaluate the impact of dose reporting schemes and tissue/applicator heterogeneities for 192Ir-, 75Se-, and 169Yb-based MRI-guided conventional and intensity-modulated brachytherapy. Methods and Materials Treatment plans using a variety of dose reporting and tissue/applicator segmentation schemes were generated for a cohort (n = 10) of cervical cancer patients treated with 192Ir-based Venezia brachytherapy. Dose calculations were performed using RapidBrachyMCTPS, a Geant4-based research Monte Carlo treatment planning system. Ultimately, five dose calculation scenarios were evaluated: (a) dose to water in water (Dw,w); (b) Dw,w taking the applicator material into consideration (Dw,wApp); (c) dose to water in medium (Dw,m); (d and e) dose to medium in medium with mass densities assigned either nominally per structure (Dm,m (Nom)) or voxel-by-voxel (Dm,m). Results Ignoring the plastic Venezia applicator (Dw,wApp) overestimates Dm,m by up to 1% (average) with high energy source (192Ir and 75Se) and up to 2% with 169Yb. Scoring dose to water (Dw,wApp or Dw,m) generally overestimates dose and this effect increases with decreasing photon energy. Reporting dose other than Dm,m (or Dm,m Nom) for 169Yb-based conventional and intensity-modulated brachytherapy leads to a simultaneous overestimation (up to 4%) of CTVHR D90 and underestimation (up to 2%) of bladder D2cc due to a significant dip in the mass-energy absorption ratios at the depths of nearby targets and OARs. Using a nominal mass-density assignment per structure, rather than a CT-derived voxel-by-voxel assignment for MRI-guided brachytherapy, amounts to a dose error up to 1% for all radionuclides considered. Conclusions The effects of the considered dose reporting schemes trend correspondingly between conventional and intensity-modulated brachytherapy. In the absence of CT-derived mass densities, MRI-only-based dosimetry can adequately approximate Dm,m by assigning nominal mass densities to structures. Tissue and applicator heterogeneities do not significantly impact dosimetry for 192Ir and 75Se, but do for 169Yb; dose reporting must be explicitly defined since Dw,m and Dw,w may overstate the dosimetric benefits.},

note = {_eprint: https://aapm.onlinelibrary.wiley.com/doi/pdf/10.1002/mp.14802},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{famulari_monte_2020,

title = {Monte Carlo dosimetric characterization of a new high dose rate 169 Yb brachytherapy source and independent verification using a multipoint plastic scintillator detector},

author = {Gabriel Famulari and Haydee M. Linares Rosales and Justine Dupere and David C. Medich and Luc Beaulieu and Shirin A. Enger},

doi = {10.1002/mp.14336},

issn = {2473-4209},

year  = {2020},

date = {2020-09-01},

journal = {Medical Physics},

volume = {47},

number = {9},

pages = {4563--4573},

abstract = {PURPOSE: A prototype 169 Yb source was developed in combination with a dynamic rotating platinum shield system (AIM-Brachy) to deliver intensity modulated brachytherapy (IMBT). The purpose of this study was to evaluate the dosimetric characteristics of the bare/shielded 169 Yb source using Monte Carlo (MC) simulations and perform an independent dose verification using a dosimetry platform based on a multipoint plastic scintillator detector (mPSD). 

METHODS: The TG-43U1 dosimetric parameters were calculated for the source model using RapidBrachyMCTPS. Real-time dose rate measurements were performed in a water tank for both the bare/shielded source using a custom remote afterloader. For each dwell position, the dose rate was independently measured by the three scintillators (BCF-10, BCF-12, and BCF-60). For the bare source, dose rate was measured at distances up to 3 cm away from the source over a range of 7 cm along the catheter. For the shielded source, measurements were performed with the mPSD placed at 1 cm from the source at four different azimuthal angles ( 0 ∘ , 9 0 ∘ , 18 0 ∘ , and 27 0 ∘ ). 

RESULTS: The dosimetric parameters were tabulated for the source model. For the bare source, differences between measured and calculated along-away dose rates were generally below 5-10%. Along the transverse axis, deviations were, on average (range), 3.3% (0.6-6.2%) for BCF-10, 1.7% (0.9-2.9%) for BCF-12, and 2.2% (0.3-4.4%) for BCF-60. The maximum dose rate reduction due to shielding at a radial distance of 1 cm was 88.8 ± 1.2%, compared to 83.5 ± 0.5% as calculated by MC. 

CONCLUSIONS: The dose distribution for the bare/shielded 169 Yb source was independently verified using mPSD with good agreement in regions close to the source. The 169 Yb source coupled with the partial-shielding system is an effective technique to deliver IMBT.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{famulari_can_2020,

title = {Can intermediate-energy sources lead to elevated bone doses for prostate and head & neck high-dose-rate brachytherapy?},

author = {Gabriel Famulari and Joanne Alfieri and Marie Duclos and Té Vuong and Shirin A. Enger},

doi = {10.1016/j.brachy.2019.12.004},

issn = {1873-1449},

year  = {2020},

date = {2020-04-01},

journal = {Brachytherapy},

volume = {19},

number = {2},

pages = {255--263},

abstract = {PURPOSE: Several radionuclides with high (60Co, 75Se) and intermediate (169Yb, 153Gd) energies have been investigated as alternatives to 192Ir for high-dose-rate brachytherapy. The purpose of this study was to evaluate the impact of tissue heterogeneities for these five high- to intermediate-energy sources in prostate and head & neck brachytherapy. METHODS AND MATERIALS: Treatment plans were generated for a cohort of prostate (n = 10) and oral tongue (n = 10) patients. Dose calculations were performed using RapidBrachyMCTPS, an in-house Geant4-based Monte Carlo treatment planning system. Treatment plans were simulated using 60Co, 192Ir, 75Se, 169Yb, and 153Gd as the active core of the microSelectron v2 source. Two dose calculation scenarios were presented: (1) dose to water in water (Dw,w), and (2) dose to medium in medium (Dm,m). 

RESULTS: Dw,w overestimates planning target volume coverage compared with Dm,m, regardless of photon energy. The average planning target volume D90 reduction was ∼1% for high-energy sources, whereas larger differences were observed for intermediate-energy sources (1%-2% for prostate and 4%-7% for oral tongue). Dose differences were not clinically relevant (textless5%) for soft tissues in general. Going from Dw,w to Dm,m, bone doses were increased two- to three-fold for 169Yb and four- to five-fold for 153Gd, whereas the ratio was close to ∼1 for high-energy sources. 

CONCLUSIONS: Dw,w underestimates the dose to bones and, to a lesser extent, overestimates the dose to soft tissues for radionuclides with average energies lower than 192Ir. Further studies regarding bone toxicities are needed before intermediate-energy sources can be adopted in cases where bones are in close vicinity to the tumor.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{famulari_novel_2020,

title = {A novel 169 Yb-based dynamic-shield intensity modulated brachytherapy delivery system for prostate cancer},

author = {Gabriel Famulari and Marie Duclos and Shirin A. Enger},

doi = {10.1002/mp.13959},

issn = {2473-4209},

year  = {2020},

date = {2020-03-01},

journal = {Medical Physics},

volume = {47},

number = {3},

pages = {859--868},

abstract = {PURPOSE: Intensity modulated brachytherapy (IMBT) is a novel high dose rate brachytherapy (HDR BT) technique which incorporates static or dynamic shielding to increase tumor coverage and/or spare healthy tissues. The purpose of this study is to present a novel delivery system (AIM-Brachy) design that can enable dynamic-shield IMBT for prostate cancer. 

METHODS: The AIM-Brachy system dynamically controls the rotation of platinum shields, placed within interstitial catheters, which partially collimate the radiation emitted from an 169 Yb source. Conventional HDR BT (10 Ci 192 Ir) and IMBT (18 Ci 169 Yb) plans were generated for 12 patients using an in-house column generation-based optimizer, coupled to a Geant4-based dose calculation engine, RapidBrachyMC. Treatment plans were normalized to match the same PTV D90 coverage as the clinical plan. Intershield attenuation effects were taken into account. A sensitivity analysis was performed to evaluate the dosimetric impact of systematic longitudinal source positioning errors ( ± 1 mm, ± 2 mm, and ± 3 mm) and rotational errors ( ± 5 ∘ , ± 10 ∘ , and ± 15 ∘ ) on clinically relevant parameters (PTV D90 and urethra D10 ). 

RESULTS: The platinum shield reduced the dose rate on the shielded side at 1 cm to 18.1% of the dose rate on the unshielded side. For equal PTV D90 coverage, the urethral D10 was reduced by 13.3% ± 4.7%, without change to other plan quality indices (PTV V100 , V150, V200 , bladder V75 , rectum V75 , HI, COIN). Delivery times for HDR BT and IMBT were 9.2 ± 1.6 min and 18.6 ± 4.0 min, respectively. In general, the PTV D90 was more sensitive to source positioning errors than rotational errors, while the urethral D10 was more sensitive to rotational errors than source positioning errors. For a typical range of positioning errors ( ± 1 mm, ± 5 ∘ ), the overall tolerance was textless2%. 

CONCLUSIONS: The AIM-Brachy system was proposed to deliver dynamic-shield IMBT for prostate cancer with the potential to create a low dose tunnel within the urethra. The urethra-sparing properties are desirable to minimize the occurrence and severity of urethral strictures or, alternatively, to provide a method for dose escalation.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{morcos_monte_2020,

title = {Monte Carlo dosimetry study of novel rotating MRI-compatible shielded tandems for intensity modulated cervix brachytherapy},

author = {Marc Morcos and Shirin A. Enger},

doi = {10.1016/j.ejmp.2020.02.014},

issn = {1724-191X},

year  = {2020},

date = {2020-03-01},

journal = {Physica medica: PM: an international journal devoted to the applications of physics to medicine and biology: official journal of the Italian Association of Biomedical Physics (AIFB)},

volume = {71},

pages = {178--184},

abstract = {PURPOSE: Intensity modulated brachytherapy (IMBT) with rotating metal shields enables dose modulation that can better conform to the tumor while reducing OAR doses. In this work, we investigate novel rotating shields, compatible with MRI-compatible tandems used for cervix brachytherapy. Three unique shields were evaluated using the traditional 192Ir source. Additionally, 75Se and 169Yb isotopes were investigated as alternative sources. 

METHODS: Three different IMBT shields were modeled and simulated in RapidBrachyMCTPS. Each tungsten shield was designed to fit inside a 6 mm-wide MRI-compatible tandem. The active core of the source was replaced with 192Ir, 75Se and 169Yb. Transmission factors (TFs) were calculated and defined as the dose ratio at 1 cm on opposite sides of the shielded tandem on the transverse plane. Polar and azimuthal anisotropy plots were extracted from simulations. Dose homogeneities V200%V100% were calculated for all radionuclide-shield combinations. 

RESULTS: TFs are favorable for IMBT and ranged between 12.9% and 32.2% for 192Ir, 4.0%-16.1% for 75Se and 1.2-6.4% for 169Yb for all shield designs. Average beam-widths in the polar and azimuthal directions were reduced to the range of 42°-112° and 27°-107°, respectively, for all shield-radionuclide combinations. Dose homogeneities for all the radionuclide-shield combinations were within 12% of the non-IMBT tandem. 

CONCLUSIONS: This study has quantitatively assessed the influence of various rotating cervical cancer-specific IMBT tandem shields on dosimetry. The dynamic single-channel shields and narrow beam-widths in the polar and azimuthal direction give rise to highly anisotropic distributions. Intermediate-to-high energy radionuclides, 75Se and 169Yb substantially improve the modulation capacity of IMBT and pave the way for treating large and complex cervical cancer without interstitial needle implantation.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{famulari_can_2020b,

title = {Can intermediate-energy sources lead to elevated bone doses for prostate and head & neck high-dose-rate brachytherapy?},

author = {Gabriel Famulari and Joanne Alfieri and Marie Duclos and Té Vuong and Shirin A. Enger},

url = {https://www.sciencedirect.com/science/article/pii/S1538472119306531},

doi = {10.1016/j.brachy.2019.12.004},

issn = {1538-4721},

year  = {2020},

date = {2020-03-01},

urldate = {2021-09-08},

journal = {Brachytherapy},

volume = {19},

number = {2},

pages = {255--263},

abstract = {Purpose 

Several radionuclides with high (60Co, 75Se) and intermediate (169Yb, 153Gd) energies have been investigated as alternatives to 192Ir for high-dose-rate brachytherapy. The purpose of this study was to evaluate the impact of tissue heterogeneities for these five high- to intermediate-energy sources in prostate and head & neck brachytherapy. 

Methods and Materials Treatment plans were generated for a cohort of prostate (n = 10) and oral tongue (n = 10) patients. Dose calculations were performed using RapidBrachyMCTPS, an in-house Geant4-based Monte Carlo treatment planning system. Treatment plans were simulated using 60Co, 192Ir, 75Se, 169Yb, and 153Gd as the active core of the microSelectron v2 source. Two dose calculation scenarios were presented: (1) dose to water in water (Dw,w), and (2) dose to medium in medium (Dm,m). 

Results 

Dw,w overestimates planning target volume coverage compared with Dm,m, regardless of photon energy. The average planning target volume D90 reduction was ∼1% for high-energy sources, whereas larger differences were observed for intermediate-energy sources (1%–2% for prostate and 4%–7% for oral tongue). Dose differences were not clinically relevant (textless5%) for soft tissues in general. Going from Dw,w to Dm,m, bone doses were increased two- to three-fold for 169Yb and four- to five-fold for 153Gd, whereas the ratio was close to ∼1 for high-energy sources. 

Conclusions 

Dw,w underestimates the dose to bones and, to a lesser extent, overestimates the dose to soft tissues for radionuclides with average energies lower than 192Ir. Further studies regarding bone toxicities are needed before intermediate-energy sources can be adopted in cases where bones are in close vicinity to the tumor.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{famulari_monte_2020b,

title = {Monte Carlo dosimetric characterization of a new high dose rate Yb brachytherapy source and independent verification using a multipoint plastic scintillator detector},

author = {Gabriel Famulari and Haydee M. Linares Rosales and Justine Dupere and David C. Medich and Luc Beaulieu and Shirin A. Enger},

url = {https://aapm.onlinelibrary.wiley.com/doi/abs/10.1002/mp.14336},

doi = {10.1002/mp.14336},

issn = {2473-4209},

year  = {2020},

date = {2020-01-01},

urldate = {2021-09-08},

journal = {Medical Physics},

volume = {47},

number = {9},

pages = {4563--4573},

abstract = {Purpose A prototype Yb source was developed in combination with a dynamic rotating platinum shield system (AIM-Brachy) to deliver intensity modulated brachytherapy (IMBT). The purpose of this study was to evaluate the dosimetric characteristics of the bare/shielded Yb source using Monte Carlo (MC) simulations and perform an independent dose verification using a dosimetry platform based on a multipoint plastic scintillator detector (mPSD). Methods The TG-43U1 dosimetric parameters were calculated for the source model using RapidBrachyMCTPS. Real-time dose rate measurements were performed in a water tank for both the bare/shielded source using a custom remote afterloader. For each dwell position, the dose rate was independently measured by the three scintillators (BCF-10, BCF-12, and BCF-60). For the bare source, dose rate was measured at distances up to 3 cm away from the source over a range of 7 cm along the catheter. For the shielded source, measurements were performed with the mPSD placed at 1 cm from the source at four different azimuthal angles ( , 9 , 18 , and 27 ). Results The dosimetric parameters were tabulated for the source model. For the bare source, differences between measured and calculated along-away dose rates were generally below 5–10%. Along the transverse axis, deviations were, on average (range), 3.3% (0.6–6.2%) for BCF-10, 1.7% (0.9–2.9%) for BCF-12, and 2.2% (0.3–4.4%) for BCF-60. The maximum dose rate reduction due to shielding at a radial distance of 1 cm was 88.8 ± 1.2%, compared to 83.5 ± 0.5% as calculated by MC. Conclusions The dose distribution for the bare/shielded Yb source was independently verified using mPSD with good agreement in regions close to the source. The Yb source coupled with the partial-shielding system is an effective technique to deliver IMBT.},

note = {_eprint: https://aapm.onlinelibrary.wiley.com/doi/pdf/10.1002/mp.14336},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{famulari_novel_2020b,

title = {A novel 169Yb-based dynamic-shield intensity modulated brachytherapy delivery system for prostate cancer},

author = {Gabriel Famulari and Marie Duclos and Shirin A. Enger},

url = {https://aapm.onlinelibrary.wiley.com/doi/abs/10.1002/mp.13959},

doi = {10.1002/mp.13959},

issn = {2473-4209},

year  = {2020},

date = {2020-01-01},

urldate = {2021-09-08},

journal = {Medical Physics},

volume = {47},

number = {3},

pages = {859--868},

abstract = {Purpose Intensity modulated brachytherapy (IMBT) is a novel high dose rate brachytherapy (HDR BT) technique which incorporates static or dynamic shielding to increase tumor coverage and/or spare healthy tissues. The purpose of this study is to present a novel delivery system (AIM-Brachy) design that can enable dynamic-shield IMBT for prostate cancer. Methods The AIM-Brachy system dynamically controls the rotation of platinum shields, placed within interstitial catheters, which partially collimate the radiation emitted from an 169Yb source. Conventional HDR BT (10 Ci 192Ir) and IMBT (18 Ci 169Yb) plans were generated for 12 patients using an in-house column generation-based optimizer, coupled to a Geant4-based dose calculation engine, RapidBrachyMC. Treatment plans were normalized to match the same PTV D90 coverage as the clinical plan. Intershield attenuation effects were taken into account. A sensitivity analysis was performed to evaluate the dosimetric impact of systematic longitudinal source positioning errors (1 mm, 2 mm, and 3 mm) and rotational errors (5, 10 and 15) on clinically relevant parameters (PTV D90 and urethra D10). Results The platinum shield reduced the dose rate on the shielded side at 1 cm to 18.1% of the dose rate on the unshielded side. For equal PTV D90 coverage, the urethral D10 was reduced by 13.3% 4.7%, without change to other plan quality indices (PTV V100, V150, V200, bladder V75, rectum V75, HI, COIN). Delivery times for HDR BT and IMBT were 9.2 ± 1.6 min and 18.6 ± 4.0 min, respectively. In general, the PTV D90 was more sensitive to source positioning errors than rotational errors, while the urethral D10 was more sensitive to rotational errors than source positioning errors. For a typical range of positioning errors (1 mm, 5), the overall tolerance was textless2%. Conclusions The AIM-Brachy system was proposed to deliver dynamic-shield IMBT for prostate cancer with the potential to create a low dose tunnel within the urethra. The urethra-sparing properties are desirable to minimize the occurrence and severity of urethral strictures or, alternatively, to provide a method for dose escalation.},

note = {_eprint: https://aapm.onlinelibrary.wiley.com/doi/pdf/10.1002/mp.13959},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{shoemaker_dosimetric_2019,

title = {Dosimetric Considerations for Ytterbium-169, Selenium-75, and Iridium-192 Radioisotopes in High-Dose-Rate Endorectal Brachytherapy},

author = {Tristan Shoemaker and Té Vuong and Harry Glickman and Samar Kaifi and Gabriel Famulari and Shirin A. Enger},

doi = {10.1016/j.ijrobp.2019.07.003},

issn = {1879-355X},

year  = {2019},

date = {2019-11-01},

journal = {International Journal of Radiation Oncology, Biology, Physics},

volume = {105},

number = {4},

pages = {875--883},

abstract = {PURPOSE: To investigate differences between prescribed and postimplant calculated dose in 192Ir high-dose-rate endorectal brachytherapy (HDR-EBT) by evaluating dose to clinical target volume (CTV) and organs at risk (OARs) calculated with a Monte Carlo-based dose calculation software, RapidBrachyMC. In addition, dose coverage, conformity, and homogeneity were compared among the radionuclides 192Ir, 75Se, and 169Yb for use in HDR-EBT. 

METHODS AND MATERIALS: Postimplant dosimetry was evaluated using 23 computed tomography (CT) images from patients treated with HDR-EBT using the 192Ir microSelectron v2 (Elekta AB, Stockholm, Sweden) source and the Intracavitary Mold Applicator Set (Elekta AB, Stockholm, Sweden), which is a flexible applicator capable of fitting a tungsten rod for OAR shielding. Four tissue segmentation schemes were evaluated: (1) TG-43 formalism, (2) materials and nominal densities assigned to contours of foreign objects, (3) materials and nominal densities assigned to contoured organs in addition to foreign objects, and (4) materials specified as in (3) but with voxel mass densities derived from CT Hounsfield units. Clinical plans optimized for 192Ir were used, with the results for 75Se and 169Yb normalized to the D90 of the 192Ir clinical plan. RESULTS: In comparison to segmentation scheme 4, TG-43-based dosimetry overestimates CTV D90 by 6% (P = .00003), rectum D50 by 24% (P = .00003), and pelvic bone D50 by 5% (P = .00003) for 192Ir. For 169Yb, CTV D90 is overestimated by 17% (P = .00003) and rectum D50 by 39% (P = .00003), and pelvic bone D50 is significantly underestimated by 27% (P = .007). Postimplant dosimetry calculations also showed that a 169Yb source would give 20% (P = .00003) lower rectum V60 and 17% (P = .00008) lower rectum D50. 

CONCLUSIONS: Ignoring high-Z materials in dose calculation contributes to inaccuracies that may lead to suboptimal dose optimization and disagreement between prescribed and calculated dose. This is especially important for low-energy radionuclides. Our results also show that with future magnetic resonance imaging-based treatment planning, loss of CT density data will only affect calculated dose in nonbone OARs by 2% or less and bone OARs by 13% or less across all sources if material composition and nominal mass densities are correctly assigned.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}