The clinical impact of removing rectal gas on high-dose-rate brachytherapy dose distributions for gynecologic cancers.

PURPOSE: To evaluate the impact of gas removal on bladder and rectal doses during intracavitary and interstitial high-dose-rate brachytherapy (HDRB) for gynecologic cancers. MATERIAL AND METHODS: Fifteen patients treated with definitive external beam radiation followed by HDRB for gynecologic cancers for a total of 21 fractions, presented with a significant amount of rectal gas at initial CT imaging (CTGAS ) after implantation. The gas was removed via rectal tubing followed by subsequent scan acquisition (CTCLINICAL ), which was used for planning and treatment delivery. To assess the effect of gas removal on dosimetry, both bladder and rectum volumes were recontoured on CTGAS . In order to evaluate the clinical impact on the total Equivalent-Dose-in-2Gy-fraction (EQD2 ), each fraction was also replanned to maintain clinically delivered target coverage (HRCTV D90). EQD2 D2cm3 for bladder and rectum were compared between plans. The Wilcoxon signed rank test was performed to evaluate statistically significant differences for all comparisons (P < 0.05). RESULTS: Mean rectum and bladder Dmax , D0.1cm3 , D1cm3 , D2cm3 , and D5cm3 were significantly different between CTGAS and CTCLINICAL . The mean percent increases on CTGAS for bladder were 12.3, 8.4, 9.9, 10.2, and 9.5% respectively and for rectum were 27.0, 19.6, 18.1, 18.5, and 19.4%, respectively. After replanning with CTGAS to maintain HRCTV D90 EQD2 , bladder and rectum EQD2 D2 cm3 resulted in significantly higher doses. The mean EQD2 D2 cm3 difference was 2.4 and 4.1 Gy for bladder and rectum, revealing a higher impact of gas removal on rectal DVH. CONCLUSION: Rectal gas removal resulted in statistically significant differences for both bladder and rectum. The resulting larger EQD2 D2 cm3 for bladder and rectum demonstrates that if patients were treated without removing gas, target coverage would need to be sacrificed to satisfy the rectum constraints and prevent toxicities. Therefore, this study demonstrates the importance of gas removal for gynecologic HDRB patients.

The combined use of 2D scout and 3D axial CT images to accurately determine the catheter tips for high-dose-rate brachytherapy plans.

PURPOSE: To develop a method combining CT scout images with axial images to improve the localization accuracy of catheter tips in high-dose-rate (HDR) brachytherapy treatments. MATERIALS AND METHODS: CT scout images were utilized along with conventionally reconstructed axial images to aid the localization of catheter tips used during HDR treatment planning. A method was developed to take advantage of the finer image resolution of the scout images to more precisely identify the tip coordinates. The accuracies of this method were compared with the conventional method based on the axial CT images alone, for various slice thicknesses, in a computed tomography dose index (CTDI) head phantom. A clinical case which involved multiple interstitial catheters was also selected for the evaluation of this method. Locations of the catheter tips were reconstructed with the conventional CT-based method and this newly developed method, respectively. Location coordinates obtained via both methods were quantitatively compared. RESULTS: Combination of the scout and axial CT images improved the accuracy of identification and reconstruction of catheter tips along the longitudinal direction (i.e., head-to-foot direction, more or less parallel to the catheter tracks), compared to relying on the axial CT images alone. The degree of improvement was dependent on CT slice thickness. For the clinical patient case, the coordinate differences of the reconstructed catheter tips were 2.6 mm ± 0.9 mm in the head-to-foot direction, 0.4 mm ± 0.2 mm in the left-to-right direction, and 0.6 mm ± 0.2 mm in the anterior-to-posterior direction, respectively. CONCLUSION: Combining CT scout and axial images demonstrates the ability to provide a more accurate identification and reconstruction of the interstitial catheter tips for HDR brachytherapy treatment, especially in the longitudinal direction. The method developed in this work has the potential to be implemented clinically together with automatic segmentation in modern brachytherapy treatment planning systems, in order to improve the reconstruction accuracy of HDR catheters.

Using gEUD based plan analysis method to evaluate proton vs. photon plans for lung cancer radiation therapy.

The goal of this study was to exam the efficacy of current DVH based clinical guidelines draw from photon experience for lung cancer radiation therapy on proton therapy. Comparison proton plans and IMRT plans were generated for 10 lung patients treated in our proton facility. A gEUD based plan evaluation method was developed for plan evaluation. This evaluation method used normal lung gEUD(a) curve in which the model parameter "a" was sampled from the literature reported value. For all patients, the proton plans delivered lower normal lung V5 Gy with similar V20 Gy and similar target coverage. Based on current clinical guidelines, proton plans were ranked superior to IMRT plans for all 10 patients. However, the proton and IMRT normal lung gEUD(a) curves crossed for 8 patients within the tested range of "a", which means there was a possibility that proton plan would be worse than IMRT plan for lung sparing. A concept of deficiency index (DI) was introduced to quantify the probability of proton plans doing worse than IMRT plans. By applying threshold on DI, four patients' proton plan was ranked inferior to the IMRT plan. Meanwhile if a threshold to the location of curve crossing was applied, 6 patients' proton plan was ranked inferior to the IMRT plan. The contradictory ranking results between the current clinical guidelines and the gEUD(a) curve analysis demonstrated there is potential pitfalls by applying photon experience directly to the proton world. A comprehensive plan evaluation based on radio-biological models should be carried out to decide if a lung patient would really be benefit from proton therapy.

Investigation on using high-energy proton beam for total body irradiation (TBI).

This work investigated the possibility of using proton beam for total body irradia-tion (TBI). We hypothesized the broad-slow-rising entrance dose from a monoen-ergetic proton beam can deliver a uniform dose to patient with varied thickness. Comparing to photon-based TBI, it would not require any patient-specific com-pensator or beam spoiler. The hypothesis was first tested by simulating 250 MeV, 275 MeV, and 300 MeV protons irradiating a wedge-shaped water phantom in a paired opposing arrangement using Monte Carlo (MC) method. To allow ± 7.5% dose variation, the maximum water equivalent thickness (WET) of a treatable patient separation was 29 cm for 250 MeV proton, and > 40 cm for 275 MeV and 300 MeV proton. The compared 6 MV photon can only treat patients with up to 15.5 cm water-equivalent separation. In the second step, we simulated the dose deposition from the same beams on a patient's whole-body CT scan. The maximum patient separation in WET was 23 cm. The calculated whole-body dose variations were ± 8.9%, ± 9.0%, ± 9.6%, and ± 14% for 250 MeV proton, 275 MeV proton, 300 MeV proton, and 6 MV photon. At last, we tested the current machine capability to deliver a monoenergetic proton beam with a large uniform field. Experiments were performed on a compact double scattering single-gantry proton system. With its C-shaped gantry design, the source-to-surface distance (SSD) reached 7 m. The measured dose deposition curve had 22 cm relatively flat entrance region. The full width half maximum field size was measured 105 cm. The current scatter filter had to be redesigned to produce a uniform intensity at such treatment distance. In con-clusion, this work demonstrated the possibility of using proton beam for TBI. The current commercially available proton machines would soon be ready for such task.

ETV6-NTRK3 is a common chromosomal rearrangement in radiation-associated thyroid cancer.

BACKGROUND: In their previous analysis of papillary thyroid carcinomas (PTCs) from an Ukrainian-American cohort that was exposed to iodine-131 ((131) I) from the Chernobyl accident, the authors identified RET/PTC rearrangements and other driver mutations in 60% of tumors. METHODS: In this study, the remaining mutation-negative tumors from that cohort were analyzed using RNA sequencing (RNA-Seq) and reverse transcriptase-polymerase chain reaction to identify novel chromosomal rearrangements and to characterize their relation with radiation dose. RESULTS: The ETS variant gene 6 (ETV6)-neurotrophin receptor 3 (NTRK3) rearrangement (ETV6-NTRK3) was identified by RNA-Seq in a tumor from a patient who received a high (131) I dose. Overall, the rearrangement was detected in 9 of 62 (14.5%) post-Chernobyl PTCs and in 3 of 151 (2%) sporadic PTCs (P = .019). The most common fusion type was between exon 4 of ETV6 and exon 14 of NTRK3. The prevalence of ETV6-NTRK3 rearrangement in post-Chernobyl PTCs was associated with increasing (131) I dose, albeit at borderline significance (P = .126). The group of rearrangement-positive PTCs (ETV6-NTRK3, RET/PTC, PAX8-PPARγ) was associated with significantly higher dose response compared with the group of PTCs with point mutations (BRAF, RAS; P < .001). In vitro exposure of human thyroid cells to 1 gray of (131) I and γ-radiation resulted in the formation of ETV6-NTRK3 rearrangement at a rate of 7.9 × 10(-6) cells and 3.0 × 10(-6) cells, respectively. CONCLUSIONS: The authors report the occurrence of ETV6-NTRK3 rearrangements in thyroid cancer and demonstrate that this rearrangement is significantly more common in tumors associated with exposure to (131) I and has a borderline significant dose response. Moreover, ETV6-NTRK3 rearrangement can be directly induced in thyroid cells by ionizing radiation in vitro and, thus, may represent a novel mechanism of radiation-induced carcinogenesis.

Radiomics for the Prediction of Pathological Complete Response to Neoadjuvant Chemoradiation in Locally Advanced Rectal Cancer: A Prospective Observational Trial.

(1) Background: An increasing amount of research has supported the role of radiomics for predicting pathological complete response (pCR) to neoadjuvant chemoradiation treatment (nCRT) in order to provide better management of locally advanced rectal cancer (LARC) patients. However, the lack of validation from prospective trials has hindered the clinical adoption of such studies. The purpose of this study is to validate a radiomics model for pCR assessment in a prospective trial to provide informative insight into radiomics validation. (2) Methods: This study involved a retrospective cohort of 147 consecutive patients for the development/validation of a radiomics model, and a prospective cohort of 77 patients from two institutions to test its generalization. The model was constructed using T2-weighted, diffusion-weighted, and dynamic contrast-enhanced MRI to understand the associations with pCR. The consistency of physicians' evaluations and agreement on pathological complete response prediction were also evaluated, with and without the aid of the radiomics model. (3) Results: The radiomics model outperformed both physicians' visual assessments in the prospective test cohort, with an area under the curve (AUC) of 0.84 (95% confidence interval of 0.70-0.94). With the aid of the radiomics model, a junior physician could achieve comparable performance as a senior oncologist. (4) Conclusion: We have built and validated a radiomics model with pretreatment MRI for pCR prediction of LARC patients undergoing nCRT.

Determination of optimal fiducial marker across image-guided radiation therapy (IGRT) modalities: visibility and artifact analysis of gold, carbon, and polymer fiducial markers.

The purpose of this study was to evaluate the visibility and artifact created by gold, carbon, and polymer fiducial markers in a simple phantom across computed tomography (CT), kilovoltage (kV), and megavoltage (MV) linear accelerator imaging and MV tomotherapy imaging. Three types of fiducial markers (gold, carbon, and polymer) were investigated for their visibility and artifacts in images acquired with various modalities and with different imaging parameters (kV, mAs, slice thickness). The imaging modalities include kV CT, 2D linac-based kilovoltage and megavoltage X-ray imaging systems, kV cone-beam CT, and normal and fine tomotherapy imaging. The images were acquired on a phantom constructed using Superflab bolus in which markers of each type were inserted into the center layer. The visibility and artifacts produced by each marker were assessed qualitatively and quantitatively. All tested markers could be identified clearly on the acquired CT and linac-based kV images; gold markers demonstrated the highest contrast. On the CT images, gold markers produced a significant artifact, while no artifacts were observed with polymer markers. Only gold markers were visible when using linac-based MV and tomotherapy imaging. For linac-based kV images, the contrast increased with kV and mAs values for all the markers, with the gold being the most pronounced. On CT images, the contrast increased with kV for the gold markers, while decreasing for the polymer and carbon marker. With the bolus phantom used, we found that when kV imaging-based treatment verification equipment is available, polymer and carbon markers may be the preferred choice for target localization and patient treatment positioning verification due to less image artifacts. If MV imaging will be the sole modality for positioning verification, it may be necessary to use gold markers despite the artifacts they create on the simulation CT images.

The impact of phantom design and material-dependence on repeatability and reproducibility of CT-based radiomics features.

PURPOSE: To understand the design of radiomics phantom and material-dependence on repeatability and reproducibility of computed tomography (CT) radiomics features. METHODS: A radiomics phantom consisting of various materials with uniformity, textural, and biological components, was constructed. The phantom was scanned with different manufacturer CT scanners and the scans were repeated multiple times on the same scanner with different acquisition settings as kVp, mAs, orientation, field of view (FOV), slice thickness, pitch, reconstruction kernels, and acquisition mode. A total of 72 phantom scans were included. For each scan, 18 different regions of interest (ROI) were contoured and 708 radiomics features were extracted from each ROI via an open source radiomics tool, IBEX. To relate the phantom data to patient data, the radiomics features from different phantom materials were compared with those extracted from 50 patients' images of five disease sites as brain, head-and-neck, breast, liver, and lung cases using box-plots comparison and principal component analysis (PCA). The temporal stability of imaging features was then evaluated with respect to a controlled scenario (test-retest) via the intra-class correlation coefficient (ICC). The reproducibility of radiomics features with respect to different scanners or acquisition settings were further evaluated with concordance correlation coefficients (CCC). RESULTS: Among all phantom materials, the biological component had feature values closest to human tissues, especially for tumors in brain and liver. The textural component showed similar ranges of variation to lung lesions, particularly for cartridges of rice, cereal, and the 3D-printed textural phantom with fine and rough grid. It also showed that certain materials, such as polystyrene foam, plaster, and peanuts, did not have comparable values to human tissue and could be excluded for future phantom design. High repeatability was observed in the test-retest study as indicated by an ICC value of 0.998 ± 0.020. All materials were used for feature stability analysis. For the inter-scanner study, shape-related features were the most-reliable category with 94% of features having CCC ≥ 0.9, while gradient orientation histogram (GOH) were the least-reliable with only 14.6% meeting the criteria. For the intra-scanner study, the reproducibility of CT-based radiomics features showed material-dependence. In general, the instability of radiomics features introduced by kVp, mAs, pitch, acquisition mode, and orientation were relatively mild. However, the homogeneous materials were more vulnerable to those changes compared to materials with textural patterns. Regardless of material compositions, resolution parameters like FOV and slice thickness, could have large impact on feature stability. Switching between standard and bone reconstruction kernels could also result significant changes to feature reproducibility. CONCLUSION: We have built a radiomics phantom using materials that cover a wide span of tumor textures seen in oncological CT images. The designed phantom presents a preliminary opportunity for investigating reproducibility of radiomics features and the reproducibility can be material dependent. Thus, in the radiomics quality assurance design, it is important to choose appropriate materials that can provide a close range of radiomics features to patients with specific disease sites dependency taken into consideration.

Optimization of treatment isocenter location in single-isocenter LINAC-based stereotactic radiosurgery for management of multiple brain metastases.

PURPOSE: Single-isocenter linear accelerator (LINAC)-based stereotactic radiosurgery (SRS) has become a promising treatment technique for the management of multiple brain metastases. Because of the high prescription dose and steep dose gradient, SRS plans are sensitive to geometric errors, resulting in loss of target coverage and suboptimal local tumor control. Current planning techniques rely on adding a uniform and isotropic setup margin to all gross tumor volumes (GTVs) to account for rotational uncertainties. However, this setup margin may be insufficient, since the magnitude of rotational uncertainties varies and is dependent upon the distance between a GTV and the isocenter. In this study, we designed a framework to determine the optimal isocenter of a single-isocenter SRS plan for multiple brain metastases using stochastic optimization to mitigate potential errors resulting from rotational uncertainties. METHODS: Planning target volumes (PTVs), defined as GTVs plus a 1-mm margin following common SRS planning convention, were assumed to be originally treated with a prescription dose and therefore covered by the prescription isodose cloud. The dose distribution, including the prescription isodose, was considered invariant assuming small rotations throughout the study. A stochastic optimization scheme was developed to determine the location of the optimal isocenter, so that the prescription dose coverage of rotated GTVs, equivalent to the intersecting volumes between the rotated GTVs and original PTVs, was maximized for any random small rotations about the isocenter. To evaluate the coverage of GTVs, the expected V 100 % undergoing random rotations was approximated as the sample average V 100 % undergoing a predetermined number of rotations. The expected V 100 % of each individual GTV and total GTVs was then compared between the plans using the optimal isocenter and the center-of-mass (CoM), respectively. RESULTS: Twenty-two patients previously treated for multiple brain metastases in a single institute were included in this retrospective study. Each patient was initially treated for more than three brain metastases (mean: 7.6; range: 3-15) with the average GTV volume of 0.89 cc (range: 0.03-11.78 cc). The optimal isocenter found for each patient was significantly different from the CoM, with the average Euclidean distance between the optimal isocenter and the CoM being 4.36 ± 2.59 cm. The dose coverage to GTVs was also significantly improved (paired t-test; p < 0.001) when the optimal isocenter was used, with the average V 100 % of total GTVs increasing from 87.1% (standard deviation as std: 11.7%; range: 39.9-98.2%) to 94.2% (std: 5.4%; range: 77.7-99.4%). The volume of a GTV was positively correlated with the expected V 100 % regardless of the isocenter used (Spearman coefficient: ρ = 0.66 ; p < 0.001). The distance between a GTV and the isocenter was negatively correlated with the expected V 100 % when the CoM was used ( ρ = - 0.21 ; p = 0.004), however no significant correlation was found when the optimal isocenter was used ( ρ = - 0.11 ; p = 0.137). CONCLUSION: The proposed framework provides an effective approach to determine the optimal isocenter of single-isocenter LINAC-based SRS plans for multiple brain metastases. The implementation of the optimal isocenter results in SRS plans with consistently higher target coverage despite potential rotational uncertainties, and therefore significantly improves SRS plan robustness against random rotational uncertainties.

SWOT Analysis and Preliminary Study on Prevention and Control Management of Temporary Integrated Isolation Ward During COVID-19 Outbreak.

Background: The world faced crises of prevention and control and shortage of medical resources during the COVID-19 (Corona Virus Disease 2019) outbreak. The establishment of temporary integrated isolation wards in hospitals, which is universal and representative in China, is one of the most-effective strategies in solving these problems according to China's experiences. Aim: To conduct a preliminary study on the establishment of a temporary integrated isolation ward during the outbreak of COVID-19 and to evaluate related impact. Methods: SWOT analysis was used to analyze the advantages, disadvantages, opportunities, and risks in the establishment of the temporary integrated isolation ward, and corresponding corrective measures were made according to the analysis results. Findings: The ward has formulated more than 10 related work procedures and prevention and control measures. A total of 93 patients with 18 critically ill patients were admitted for treatment and isolation. They were all evaluated based on established procedures and protocols. Twenty-four supplementary nucleic acid tests were ordered and conducted. One new patient with COVID-19 was confirmed and was successfully transferred to the designated COVID-19 infectious control hospital. There were no missed diagnosis or misdiagnosis, no cross-infection of patients, no cluster outbreak, and no infection of medical workers during the entire process. Conclusion: SWOT analysis is helpful in guiding the establishment of a temporary integrated isolation ward and the formulation of prevention and control measures in Hebei General Hospital during the COVID-19 outbreak. It provides the guidance and reference of significance for the establishment of similar types of wards in the future.

Proton beam therapy can achieve lower vertebral bone marrow dose than photon beam therapy during chemoradiation therapy of esophageal cancer.

Chemoradiation therapy plays an important role in both the neoadjuvant and definitive management of esophageal cancer (EC). Prior studies have suggested that advanced planning techniques can better spare organs at risk including the heart. Although multiple toxicities can result from esophageal radiotherapy, one less studied acute toxicity is that of myelosuppression, which can result, in part, from the combination of chemotherapy and incidental radiotherapy administration to the vertebral bodies (VBs), which abut the posterior aspect of the esophagus, especially in the lower thoracic esophagus. Traditionally, VB bone marrow doses are not accounted during EC radiation therapy planning. We sought to compare the doses to VBs between proton and photon radiation therapy as part of chemoradiation therapy for EC treatment. By reducing doses to the vertebrae, radiation therapy can decrease treatment-related myelosuppression, which can avoid delays or chemotherapy dose reductions in therapy, which likely affect long-term patient survival. Dose constraints are not routinely employed for bone marrow in radiation treatment planning. In our previous work, we identified thresholds to avoid grade ≥3 leukopenia, including VB V10Gy, VB V20Gy, and a mean VB dose (MVD) of 18.8 Gy. Herein we perform a retrospective dosimetric planning study comparing passive- or double-scattering proton beam therapy (PS-PBT), volumetric-modulated arc therapy (VMAT) (photon-based), and intensity-modulated radiation therapy (IMRT) (photon-based) in 25 patients with locally advanced EC who were treated originally with photon RT at our institution between 2011 and 2016. The aforementioned dose constraints were included in the retrospective planning process for PS-PBT, VMAT, and IMRT to determine the feasibility of achieving these VB constraints while maintaining reasonable target coverage and planned, consistent constraints to other organs at risk including lungs, spinal cord, and stomach. PS-PBT plans were found to achieve lower doses for VB V10Gy, V20Gy, and MVD than VMAT and static IMRT plans while achieving the same target coverage. PS-PBT resulted in lower organs at risk dosimetric parameters than the photon plans, with p < 0.0001. Student's paired t-test p-values in favor of proton therapy's ability to spare organs were as follows: for PS-PBT vs VMAT and PS-PBT vs IMRT in mean doses for lung, liver, and VB and VB V10Gy and VB V20Gy were all <0.001 (Bonferroni corrected α=0.017). One-way ANOVA found that VB doses (VB V10Gy, VB V20Gy, and MVD) were significantly lower for proton therapy (p < 0.006) among the 3 planning techniques.

A novel approach to embed eye shields in customized bolus on nasal dorsum treatment for electron radiotherapy.

We aim to demonstrate the unique use of embedded lead eye shields in an electron wax bolus when treating the nasal dorsum. A patient presented to the clinic with squamous cell carcinoma of the nasal dorsum requiring treatment with en face electrons. A 3D customized wax bolus was designed and imported into the treatment planning system (TPS) to calculate the dose distribution. Due to high lens dose, the bolus was customized further to create 2 milled open slots in the wax, over the lens of eye, to allow lead sheets totaling 4 mm to be slid into the wax. The patient was brought back to the clinic to be scanned with the wax bolus fitting snugly over the nose, eyes, and cheek regions. The 3D milled insert holes were contoured on the CT in the TPS, assigned HU of 2758, to mimic the lead insertion. The lens dose with lead inserts was compared to the plan without lead insert. To further confirm the lens dose, EBT3 films were placed on the right and left eye under the bolus, and nose dorsum on the first day of treatment. The maximum dose of right lens, as calculated in the TPS with the simulated lead shields in place, decreased from 989.5cGy to 457cGy. The maximum dose of left lens decreased from 1085.4cGy to 501cGy. The dose readings from EBT3 films were in good agreement with the TPS, with deviation of 3.32%, 0.26%, and 3.44% for right lens, left lens, and nose, respectively. Daily positioning deviations compared to the plan were 0.65 ± 0.16cm and 0.63 ± 0.29cm for right eye and left eye, respectively. This novel device demonstrated the feasibility, in terms of dose calculation accuracy in the TPS and fabrication, of using customized bolus with lead inserts to conveniently shield the lens of the eyes in an electron treatment for the nose, enabling a streamlined daily setup.

Adaptive Radiation Therapy in the Treatment of Lung Cancer: An Overview of the Current State of the Field.

Lung cancer treatment is constantly evolving due to technological advances in the delivery of radiation therapy. Adaptive radiation therapy (ART) allows for modification of a treatment plan with the goal of improving the dose distribution to the patient due to anatomic or physiologic deviations from the initial simulation. The implementation of ART for lung cancer is widely varied with limited consensus on who to adapt, when to adapt, how to adapt, and what the actual benefits of adaptation are. ART for lung cancer presents significant challenges due to the nature of the moving target, tumor shrinkage, and complex dose accumulation because of plan adaptation. This article presents an overview of the current state of the field in ART for lung cancer, specifically, probing topics of: patient selection for the greatest benefit from adaptation, models which predict who and when to adapt plans, best timing for plan adaptation, optimized workflows for implementing ART including alternatives to re-simulation, the best radiation techniques for ART including magnetic resonance guided treatment, algorithms and quality assurance, and challenges and techniques for dose reconstruction. To date, the clinical workflow burden of ART is one of the major reasons limiting its widespread acceptance. However, the growing body of evidence demonstrates overwhelming support for reduced toxicity while improving tumor dose coverage by adapting plans mid-treatment, but this is offset by the limited knowledge about tumor control. Progress made in predictive modeling of on-treatment tumor shrinkage and toxicity, optimizing the timing of adaptation of the plan during the course of treatment, creating optimal workflows to minimize staffing burden, and utilizing deformable image registration represent ways the field is moving toward a more uniform implementation of ART.

Dosimetric comparison and biological evaluation of fixed-jaw intensity-modulated radiation therapy for T-shaped esophageal cancer.

BACKGROUND: To evaluate the dosimetric and biological benefits of the fixed-jaw (FJ) intensity-modulated radiation therapy (IMRT) technique for patients with T-shaped esophageal cancer. METHODS: FJ IMRT plans were generated for thirty-five patients and compared with jaw tracking (JT) IMRT, static jaw (SJ) IMRT and JT volumetric modulated arc therapy (VMAT). Dosimetric parameters, tumor control probability (TCP) and normal tissue complication probability (NTCP), monitor units (MUs), delivery time and gamma passing rate, as a measure of dosimetric verification, were compared. The correlation between the length of PTV-C below the upper boundary of lung tissue (PTV-Cinferior) and dosimetric parameters and NTCP of the lung tissue were analyzed. RESULTS: The homogeneity and conformity of the target in the four plans were basically equivalent. When compared to the JT IMRT and SJ IMRT plans, FJ IMRT plan led to a statistically significant improvement in the NTCP and low-middle dosimetric parameters of the lung, and the improvement had a moderately positive correlation with the length of PTV-Cinferior, with a correlation coefficient ranging from 0.523 to 0.797; the FJ IMRT plan exhibited better lung sparing in low-dose volumes than the JT VMAT plan. The FJ IMRT plan had similar MUs (888 ± 99) and delivery times (516.1 ± 54.7 s) as the JT IMRT plan (937 ± 194, 522 ± 5.6 s) but higher than SJ IMRT (713 ± 137, 488.8 ± 45.2 s) and JT VMAT plan (517 ± 59, 263.7 ± 43.3 s). CONCLUSIONS: The FJ IMRT technique is superior in reducing the low-dose volumes of lung tissues for patients with T-shaped esophageal cancer.

Comprehensive Output Estimation of Double Scattering Proton System With Analytical and Machine Learning Models.

OBJECTIVES: The beam output of a double scattering proton system varies for each combination of beam option, range, and modulation and therefore is difficult to be accurately modeled by the treatment planning system (TPS). This study aims to design an empirical method using the analytical and machine learning (ML) models to estimate proton output in a double scattering proton system. MATERIALS AND METHODS: Three analytical models using polynomial, linear, and logarithm-polynomial equations were generated on a training dataset consisting of 1,544 clinical measurements to estimate proton output for each option. Meanwhile, three ML models using Gaussian process regression (GPR) with exponential kernel, squared exponential kernel, and rational quadratic kernel were also created for all options combined. The accuracy of each model was validated against 241 additional clinical measurements as the testing dataset. Two most robust models were selected, and the minimum number of samples needed for either model to achieve sufficient accuracy ( ± 3%) was determined by evaluating the mean average percentage error (MAPE) with increasing sample number. The differences between the estimated outputs using the two models were also compared for 1,000 proton beams with a randomly generated range, and modulation for each option. RESULTS: The polynomial model and the ML GPR model with exponential kernel yielded the most accurate estimations with less than 3% deviation from the measured outputs. At least 20 samples of each option were needed to build the polynomial model with less than 1% MAPE, whereas at least a total of 400 samples were needed for all beam options to build the ML GPR model with exponential kernel to achieve comparable accuracy. The two independent models agreed with less than 2% deviation using the testing dataset. CONCLUSION: The polynomial model and the ML GPR model with exponential kernel were built for proton output estimation with less than 3% deviations from the measurements. They can be used as an independent output prediction tool for a double scattering proton beam and a secondary output check tool for a cross check between themselves.

The Impact of COVID-19 on Brachytherapy During the Pandemic: A Rutgers-Robert Wood Johnson Barnabas Health Multisite Experience.

PURPOSE: This study aimed to evaluate whether the coronavirus disease of 2019 (COVID-19) pandemic resulted in treatment delays in patients scheduled for or undergoing brachytherapy. METHODS AND MATERIALS: A retrospective cohort study was conducted across 4 affiliated sites after local institutional review board approval. The eligibility criteria were defined as all patients with cancer whose treatment plan included brachytherapy during the COVID-19 pandemic from February 24, 2020 to June 30, 2020. Treatment delays, cancellations, alterations of fractionation regimens, and treatment paradigm changes were evaluated. RESULTS: A total of 47 patients were eligible for the analysis. Median patient age at the time of treatment was 62 years (interquartile range, 56-70 years). Endometrial, cervical, and prostate cancers were the most common sites included in this analysis. Three patients (6.4%) with cervical cancer were diagnosed with COVID-19 during the course of their treatment. Interruptions of external beam radiation therapy (EBRT), cancellations of EBRT, cancellations of brachytherapy, and treatment delays due to COVID occurred in 5 (10.6%), 3 (6.4%), 8 (17%), and 9 (19%) patients, respectively. The mean and median number of days delayed for patients who experienced treatment interruptions were 16.3 days (standard deviation: 13.9 days) and 14 days (interquartile range, 5.75-23.75 days), respectively. For patients with cervical cancer, the mean and median overall treatment times defined as the time from the start of EBRT to the end of brachytherapy were 56 and 49 days, respectively. CONCLUSIONS: Despite the challenges the health care system faced during the pandemic, most patients with cancer were safely treated with minor treatment delays and interruptions. Long-term follow up is needed to assess the impact of COVID-19 and treatment interruptions on oncologic outcomes.

A 3D global-to-local deformable mesh model based registration and anatomy-constrained segmentation method for image guided prostate radiotherapy.

PURPOSE: In the external beam radiation treatment of prostate cancers, successful implementation of adaptive radiotherapy and conformal radiation dose delivery is highly dependent on precise and expeditious segmentation-and registration of the prostate volume between the simulation and the treatment images. The purpose of this study is to develop a novel, fast, and accurate segmentation and registration method to increase the computational efficiency to meet the restricted clinical treatment time requirement in image guided radiotherapy. METHODS: The method developed in this study used soft tissues to capture the transformation between the 3D planning CT (pCT) images and 3D cone-beam CT (CBCT) treatment images. The method incorporated a global-to-local deformable mesh model based registration framework as well as an automatic anatomy-constrained robust active shape model (ACRASM) based segmentation algorithm in the 3D CBCT images. The global registration was based on the mutual information method, and the local registration was to minimize the Euclidian distance of the corresponding nodal points from the global transformation of deformable mesh models, which implicitly used the information of the segmented target volume. The method was applied on six data sets of prostate cancer patients. Target volumes delineated by the same radiation oncologist on the pCT and CBCT were chosen as the benchmarks and were compared to the segmented and registered results. The distance-based and the volume-based estimators were used to quantitatively evaluate the results of segmentation and registration. RESULTS: The ACRASM segmentation algorithm was compared to the original active shape mode (ASM) algorithm by evaluating the values of the distance-based estimators. With respect to the corresponding benchmarks, the mean distance ranged from -0.85 to 0.84 mm for ACRASM and from -1.44 to 1.17 mm for ASM. The mean absolute distance ranged from 1.77 to 3.07 mm for ACRASM and from 2.45 to 6.54 mm for ASM. The volume overlap ratio ranged from 79% to 91% for ACRASM and from 44% to 80% for ASM. These data demonstrated that the segmentation results of ACRASM were in better agreement with the corresponding benchmarks than those of ASM. The developed registration algorithm was quantitatively evaluated by comparing the registered target volumes from the pCT to the benchmarks on the CBCT. The mean distance and the root mean square error ranged from 0.38 to 2.2 mm and from 0.45 to 2.36 mm, respectively, between the CBCT images and the registered pCT. The mean overlap ratio of the prostate volumes ranged from 85.2% to 95% after registration. The average time of the ACRASM-based segmentation was under 1 min. The average time of the global transformation was from 2 to 4 min on two 3D volumes and the average time of the local transformation was from 20 to 34 s on two deformable superquadrics mesh models. CONCLUSIONS: A novel and fast segmentation and deformable registration method was developed to capture the transformation between the planning and treatment images for external beam radiotherapy of prostate cancers. This method increases the computational efficiency and may provide foundation to achieve real time adaptive radiotherapy.

Development and Validation of a Radiomics Nomogram Model for Predicting Postoperative Recurrence in Patients With Esophageal Squamous Cell Cancer Who Achieved pCR After Neoadjuvant Chemoradiotherapy Followed by Surgery.

Background and purpose: Although patients with esophageal squamous cell carcinoma (ESCC) can achieve a pathological complete response (pCR) after neoadjuvant chemoradiotherapy (nCRT) followed by surgery, one-third of these patients with a pCR may still experience recurrence. The aim of this study is to develop and validate a predictive model to estimate recurrence-free survival (RFS) in those patients who achieved pCR. Materials and methods: Two hundred six patients with ESCC were enrolled and divided into a training cohort (n = 146) and a validation cohort (n = 60). Radiomic features were extracted from contrast-enhanced computed tomography (CT) images of each patient. Feature reduction was then implemented in two steps, including a multiple segmentation test and least absolute shrinkage and selection operator (LASSO) Cox proportional hazards regression method. A radiomics signature was subsequently constructed and evaluated. For better prediction performance, a clinical nomogram based on clinical risk factors and a nomogram incorporating the radiomics signature and clinical risk factors was built. Finally, the prediction models were further validated by calibration and the clinical usefulness was examined in the validation cohort to determine the optimal prediction model. Results: The radiomics signature was constructed using eight radiomic features and displayed a significant correlation with RFS. The nomogram incorporating the radiomics signature with clinical risk factors achieved optimal performance compared with the radiomics signature (P < 0.001) and clinical nomogram (P < 0.001) in both the training cohort [C-index (95% confidence interval [CI]), 0.746 (0.680-0.812) vs. 0.685 (0.620-0.750) vs. 0.614 (0.538-0.690), respectively] and validation cohort [C-index (95% CI), 0.724 (0.696-0.752) vs. 0.671 (0.624-0.718) vs. 0.629 (0.597-0.661), respectively]. The calibration curve and decision curve analysis revealed that the radiomics nomogram outperformed the other two models. Conclusions: A radiomics nomogram model incorporating radiomics features and clinical factors has been developed and has the improved ability to predict the postoperative recurrence risk in patients with ESCC who achieved pCR after nCRT followed by surgery.

Time dependence of energy spectra of brachytherapy sources and its impact on the half and the tenth value layers.

PURPOSE: Several factors including radionuclide purity influence the photon energy spectra from sealed brachytherapy sources. The existence of impurities and trace elements in radioactive materials as well as the substrate and encapsulation may not only alter the spectrum at a given time but also cause change in the spectra as a function of time. The purpose of this study is to utilize a semiempirical formalism, which quantitatively incorporates this time dependence, to calculate and evaluate the shielding requirement impacts introduced by this time dependence for a 103Pd source. METHODS: The formalism was used to calculate the NthVL thicknesses in lead for a 103Pd model 200 seed. Prior to 2005, the 103Pd in this source was purified to a level better than 0.006% of the total 103Pd activity, the key trace impurity consisting of 65Zn. Because 65Zn emits higher energy photons and has a much longer half-life of 244 days compared to 103Pd, its presence in 103Pd seeds led to a time dependence of the photon spectrum and other related physical quantities. This study focuses on the time dependence of the NthVL and the analysis of the corresponding shielding requirements. RESULTS: The results indicate that the first HVL and the first TVL in lead steadily increased with time for about 200 days and then reached a plateau. The increases at plateau were more than 1000 times compared to the corresponding values on the zeroth day. The second and third TVLs in lead reached their plateaus in about 100 and 60 days, respectively, and the increases were about 19 and 2.33 times the corresponding values on the zeroth day, respectively. All the TVLs demonstrated a similar time dependence pattern, with substantial increases and eventual approach to a plateau. CONCLUSIONS: The authors conclude that the time dependence of the emitted photon spectra from brachytherapy sources can introduce substantial variations in the values of the NthVL with time if certain impurities are present. The contribution of 65Zn to the dose rate constant was less than 0.03% in the earlier 103Pd seeds, and because of the use of new processing technologies since 2005, this impurity has been essentially eliminated, as demonstrated in the measured spectra of current 103Pd model 200 seeds. This study illustrates the importance of performing photon spectroscopy of the manufactured radioactive sources as a quality assurance test for an assessment over time of both the radiation protection and the dosimetric properties.

Dosimetric evaluations of the interplay effect in respiratory-gated intensity-modulated radiation therapy.

The interplay between a mobile target and a dynamic multileaf collimator can compromise the accuracy of intensity-modulated radiation therapy (IMRT). Our goal in this study is to investigate the dosimetric effects caused by the respiratory motion during IMRT. A moving phantom was built to simulate the typical breathing motion. Different sizes of the gating windows were selected for gated deliveries. The residual motions during the beam-on period ranged from 0.5 to 3 cm. An IMRT plan with five treatment fields from different gantry angles were delivered to the moving phantom for three irradiation conditions: Stationary condition, moving with the use of gating system, and moving without the use of gating system. When the residual motion was 3 cm, the results showed significant differences in dose distributions between the stationary condition and the moving phantom without gating beam control. The overdosed or underdosed areas enclosed about 33% of the treatment area. In contrast, the dose distribution on the moving phantom with gating window set to 0.5 cm showed no significant differences from the stationary phantom. With the appropriate setting of the gating window, the deviation of dose from the respiratory motion can be minimized. It appeals that limiting the residual motion to less than 0.5 cm is critical for the treatments of mobile structures.