Tumor radiomic heterogeneity: Multiparametric functional imaging to characterize variability and predict response following cervical cancer radiation therapy.

BACKGROUND: Robust approaches to quantify tumor heterogeneity are needed to provide early decision support for precise individualized therapy. PURPOSE: To conduct a technical exploration of longitudinal changes in tumor heterogeneity patterns on dynamic contrast-enhanced (DCE) magnetic resonance imaging (MRI), diffusion-weighted imaging (DWI) and FDG positron emission tomography / computed tomography (PET/CT), and their association to radiation therapy (RT) response in cervical cancer. STUDY TYPE: Prospective observational study with longitudinal MRI and PET/CT pre-RT, early-RT (2 weeks), and mid-RT (5 weeks). POPULATION: Twenty-one FIGO IB2 -IVA cervical cancer patients receiving definitive external beam RT and brachytherapy. FIELD STRENGTH/SEQUENCE: 1.5T, precontrast axial T1 -weighted, axial and sagittal T2 -weighted, sagittal DWI (multi-b values), sagittal DCE MRI (<10 sec temporal resolution), postcontrast axial T1 -weighted. ASSESSMENT: Response assessment 1 month after completion of treatment by a board-certified radiation oncologist from manually delineated tumor volume changes. STATISTICAL TESTS: Intensity histogram (IH) quantiles (DCE SI10% and DWI ADC10% , FDG-PET SUVmax ) and distribution moments (mean, variance, skewness, kurtosis) were extracted. Differences in IH features between timepoints and modalities were evaluated by Skillings-Mack tests with Holm's correction. Area under receiver-operating characteristic curve (AUC) and Mann-Whitney testing was performed to discriminate treatment response using IH features. RESULTS: Tumor IH means and quantiles varied significantly during RT (SUVmean : ↓28-47%, SUVmax : ↓30-59%, SImean : ↑8-30%, SI10% : ↑8-19%, ADCmean : ↑16%, P < 0.02 for each). Among IH heterogeneity features, FDG-PET SUVCoV (↓16-30%, P = 0.011) and DW-MRI ADCskewness decreased (P = 0.001). FDG-PET SUVCoV was higher than DCE-MRI SICoV and DW-MRI ADCCoV at baseline (P < 0.001) and 2 weeks (P = 0.010). FDG-PET SUVkurtosis was lower than DCE-MRI SIkurtosis and DW-MRI ADCkurtosis at baseline (P = 0.001). Some IH features appeared to associate with favorable tumor response, including large early RT changes in DW-MRI ADCskewness (AUC = 0.86). DATA CONCLUSION: Preliminary findings show tumor heterogeneity was variable between patients, modalities, and timepoints. Radiomic assessment of changing tumor heterogeneity has the potential to personalize treatment and power outcome prediction. LEVEL OF EVIDENCE: 2 Technical Efficacy: Stage 3 J. Magn. Reson. Imaging 2018;47:1388-1396.

Electron beam energy QA - a note on measurement tolerances.

Monthly QA is recommended to verify the constancy of high-energy electron beams generated for clinical use by linear accelerators. The tolerances are defined as 2%/2 mm in beam penetration according to AAPM task group report 142. The practical implementation is typically achieved by measuring the ratio of readings at two different depths, preferably near the depth of maximum dose and at the depth corresponding to half the dose maximum. Based on beam commissioning data, we show that the relationship between the ranges of energy ratios for different electron energies is highly nonlinear. We provide a formalism that translates measurement deviations in the reference ratios into change in beam penetration for electron energies for six Elekta (6-18 MeV) and eight Varian (6-22 MeV) electron beams. Experimental checks were conducted for each Elekta energy to compare calculated values with measurements, and it was shown that they are in agreement. For example, for a 6 MeV beam a deviation in the measured ionization ratio of ± 15% might still be acceptable (i.e., be within ± 2 mm), whereas for an 18 MeV beam the corresponding tolerance might be ± 6%. These values strongly depend on the initial ratio chosen. In summary, the relationship between differences of the ionization ratio and the corresponding beam energy are derived. The findings can be translated into acceptable tolerance values for monthly QA of electron beam energies.

Essentials and guidelines for clinical medical physics residency training programs: executive summary of AAPM Report Number 249.

There is a clear need for established standards for medical physics residency training. The complexity of techniques in imaging, nuclear medicine, and radiation oncology continues to increase with each passing year. It is therefore imperative that training requirements and competencies are routinely reviewed and updated to reflect the changing environment in hospitals and clinics across the country. In 2010, the AAPM Work Group on Periodic Review of Medical Physics Residency Training was formed and charged with updating AAPM Report Number 90. This work group includes AAPM members with extensive experience in clinical, professional, and educational aspects of medical physics. The resulting report, AAPM Report Number 249, concentrates on the clinical and professional knowledge needed to function independently as a practicing medical physicist in the areas of radiation oncology, imaging, and nuclear medicine, and constitutes a revision to AAPM Report Number 90. This manuscript presents an executive summary of AAPM Report Number 249.

Preclinical Ultra-High Dose Rate (FLASH) Proton Radiation Therapy System for Small Animal Studies.

Purpose: Animal studies with ultrahigh dose-rate radiation therapy (FLASH, >40 Gy/s) preferentially spare normal tissues without sacrificing antitumor efficacy compared with conventional dose-rate radiation therapy (CONV). At the University of Washington, we developed a cyclotron-generated preclinical scattered proton beam with FLASH dose rates. We present the technical details of our FLASH radiation system and preliminary biologic results from whole pelvis radiation. Methods and Materials: A Scanditronix MC50 compact cyclotron beamline has been modified to produce a 48.7 MeV proton beam at dose rates between 0.1 and 150 Gy/s. The system produces a 6 cm diameter scattered proton beam (flat to ± 3%) at the target location. Female C57BL/6 mice 5 to 6 weeks old were used for all experiments. To study normal tissue effects in the distal colon, mice were irradiated using the entrance region of the proton beam to the whole pelvis, 18.5 Gy at different dose rates: control, CONV (0.6-1 Gy/s) and FLASH (50-80 Gy/s). Survival was monitored daily and EdU (5-ethynyl-2´-deoxyuridine) staining was performed at 24- and 96-hours postradiation. Cleaved caspase-3 staining was performed 24-hours postradiation. To study tumor control, allograft B16F10 tumors were implanted in the right flank and received 18 Gy CONV or FLASH proton radiation. Tumor growth and survival were monitored. Results: After 18.5 Gy whole pelvis radiation, survival was 100% in the control group, 0% in the CONV group, and 44% in the FLASH group (P < .01). EdU staining showed cell proliferation was significantly higher in the FLASH versus CONV group at both 24-hours and 96-hours postradiation in the distal colon, although both radiation groups showed decreased proliferation compared with controls (P < .05). Lower cleaved caspase-3 staining was seen in the FLASH versus conventional group postradiation (P < .05). Comparable flank tumor control was observed in the CONV and FLASH groups. Conclusions: We present our preclinical FLASH proton radiation system and biologic results showing improved survival after whole pelvis radiation, with equivalent tumor control.

Validation of the Pinnacle³ photon convolution-superposition algorithm applied to fast neutron beams.

We evaluate a photon convolution-superposition algorithm used to model a fast neutron therapy beam in a commercial treatment planning system (TPS). The neutron beam modeled was the Clinical Neutron Therapy System (CNTS) fast neutron beam produced by 50 MeV protons on a Be target at our facility, and we implemented the Pinnacle3 dose calculation model for computing neutron doses. Measured neutron data were acquired by an IC30 ion chamber flowing 5 cc/min of tissue equivalent gas. Output factors and profile scans for open and wedged fields were measured according to the Pinnacle physics reference guide recommendations for photon beams in a Wellhofer water tank scanning system. Following the construction of a neutron beam model, computed doses were then generated using 100 monitor units (MUs) beams incident on a water-equivalent phantom for open and wedged square fields, as well as multileaf collimator (MLC)-shaped irregular fields. We compared Pinnacle dose profiles, central axis doses, and off-axis doses (in irregular fields) with 1) doses computed using the Prism treatment planning system, and 2) doses measured in a water phantom and having matching geometry to the computation setup. We found that the Pinnacle photon model may be used to model most of the important dosimetric features of the CNTS fast neutron beam. Pinnacle-calculated dose points among open and wedged square fields exhibit dose differences within 3.9 cGy of both Prism and measured doses along the central axis, and within 5 cGy difference of measurement in the penumbra region. Pinnacle dose point calculations using irregular treatment type fields showed a dose difference up to 9 cGy from measured dose points, although most points of comparison were below 5 cGy. Comparisons of dose points that were chosen from cases planned in both Pinnacle and Prism show an average dose difference less than 0.6%, except in certain fields which incorporate both wedges and heavy blocking of the central axis. All clinical cases planned in both Prism and Pinnacle were found to be comparable in terms of dose-volume histograms and spatial dose distribution following review by the treating clinicians. Variations were considered minor and within clinically acceptable limits by the treating clinicians. The Pinnacle TPS has sufficient computational modeling ability to adequately produce a viable neutron model for clinical use in treatment planning.

Report on the American Association of Medical Physics Undergraduate Fellowship Programs.

The American Association of Physicists in Medicine (AAPM) sponsors two summer undergraduate research programs to attract top performing undergraduate students into graduate studies in medical physics: the Summer Undergraduate Fellowship Program (SUFP) and the Minority Undergraduate Summer Experience (MUSE). Undergraduate research experience (URE) is an effective tool to encourage students to pursue graduate degrees. The SUFP and MUSE are the only medical physics URE programs. From 2001 to 2012, 148 fellowships have been awarded and a total of $608,000 has been dispersed to fellows. This paper reports on the history, participation, and status of the programs. A review of surveys of past fellows is presented. Overall, the fellows and mentors are very satisfied with the program. The efficacy of the programs is assessed by four metrics: entry into a medical physics graduate program, board certification, publications, and AAPM involvement. Sixty-five percent of past fellow respondents decided to pursue a graduate degree in medical physics as a result of their participation in the program. Seventy percent of respondents are currently involved in some educational or professional aspect of medical physics. Suggestions for future enhancements to better track and maintain contact with past fellows, expand funding sources, and potentially combine the programs are presented.

Effects of tissue heterogeneity and comparisons of collapsed cone and Monte Carlo fast neutron patient dosimetry using the University of Washington clinical neutron therapy system (CNTS).

Fast neutron therapy is a high linear energy transfer (LET) radiation treatment modality offering advantages over low LET radiations. Multileaf collimator technology reduces normal-tissue dose (toxicity) and makes neutron therapy more comparable to MV x-ray treatments. Published clinical-trial and other experiences with fast neutron therapy are reported. Early comparative studies failed to consider differences in target-dose spatial conformality between x-ray and neutron treatments, which is especially important for organs-at-risk close to tumor targets. Treatments planning systems (TPS) for high-energy neutrons lag behind TPS tools for MV x-rays, creating challenges for comparative studies of clinical outcomes. A previously published Monte Carlo model of the University of Washington (UW) Clinical Neutron Therapy System (CNTS) is refined and integrated with the RayStation TPS as an external dose planning/verification tool. The collapsed cone (CC) dose calculations in the TPS are based on measured dose profiles and output factors in water, with the absolute dose determined using a tissue-equivalent ionization chamber. For comparison, independent (external) Monte Carlo simulation computes dose on a voxel-by-voxel basis using an atlas that maps Hounsfield Unit (HU) numbers to elemental composition and density. Although the CC algorithm in the TPS accurately computes neutron dose to water compared to Monte Carlo calculations, calculated dose to water differs from bone or tissue depending largely on hydrogen content. Therefore, the elemental composition of tissue and bone, rather than the material or electron density, affects fast neutron dose. While the CC algorithm suffices for reproducible patient dosimetry in fast neutron therapy, adopting methods that consider tissue heterogeneity would enhance patient-specific neutron dose accuracy relative to national standards for other types of ionizing radiation. Corrections for tissue composition have a significant impact on absolute dose and the relative biological effectiveness (RBE) of neutron treatments compared to other radiation types (MV x-rays, protons, and carbon ions).

Scattering kernels for fast neutron therapy treatment planning.

The University of Washington (UW) Clinical Neutron Therapy System (CNTS) has been used to treat over 3300 patients. Treatment planning for these patients is currently performed using an MV x-ray model in Pinnacle® adapted to fit measurements of fast neutron output factors, wedge factors, depth-dose and lateral profiles. While this model has provided an adequate representation of the CNTS for 3D conformal treatment planning, later versions of Pinnacle did not allow for isocentric gantry rotation machines with a source-to-axis distance of 150 cm. This restriction limited the neutron model to version 9.0 while the photon and electron treatment planning at the UW had moved on to newer versions. Also, intensity modulated neutron therapy (IMNT) is underdevelopment at the UW, and the Pinnacle neutron model developed cannot be used for inverse treatment planning. We have used the MCNP6 Monte Carlo code system to develop Collapsed Cone (CC) and Singular Value Decomposition (SVD) neutron scattering kernels suitable for integration into the RayStation treatment planning system. Kernels were generated for monoenergetic neutrons with energies ranging from 1 keV to 51 MeV, i.e. the energy range most relevant to the CNTS. Percent depth dose (PDD) profiles computed in RayStation for the CC and SVD kernels are in excellent agreement with each other for depths beyond the beam's dose build-up region (depths greater than about 1.7 cm) for open 2.8 × 2.8 cm2, 10.3 × 10.3 cm2, and 28.8 × 32.8 cm2 fields. Lateral profiles at several depths, as well as the PDD, calculated using the CC kernels in RayStation for the full CNTS energy spectrum pass a 3%/3 mm gamma test for field sizes of 2.8 × 2.8 cm2, 10.0 × 10.3 cm2, and 28.8 × 32.8 cm2.

Reducing Cardiac Radiation Dose From Breast Cancer Radiation Therapy With Breath Hold Training and Cognitive Behavioral Therapy.

The delivery of radiation therapy shares many of the challenges encountered in imaging procedures. As in imaging, such as MRI, organ motion must be reduced to a minimum, often for lengthy time periods, to effectively target the tumor during imaging-guided therapy while reducing radiation dose to nearby normal tissues. For patients, radiation therapy is frequently a stress- and anxiety-provoking medical procedure, evoking fear from negative perceptions about irradiation, confinement from immobilization devices, claustrophobia, unease with equipment, physical discomfort, and overall cancer fear. Such stress can be a profound challenge for cancer patients' emotional coping and tolerance to treatment, and particularly interferes with advanced radiation therapy procedures where active, complex and repetitive high-level cooperation is often required from the patient.In breast cancer, the most common cancer in women worldwide, radiation therapy is an indispensable component of treatment to improve tumor control and outcome in both breast-conserving therapy for early-stage disease and in advanced-stage patients. High technological complexity and high patient cooperation is required to mitigate the known cardiac toxicity and mortality from breast cancer radiation by reducing the unintended radiation dose to the heart from left breast or left chest wall irradiation. To address this, radiation treatment in daily deep inspiration breath hold (DIBH), to create greater distance between the treatment target and the heart, is increasingly practiced. While holding the promise to decrease cardiac toxicity, DIBH procedures often augment patients' baseline stress and anxiety reaction toward radiation treatment. Patients are often overwhelmed by the physical and mental demands of daily DIBH, including the nonintuitive timed and sustained coordination of abdominal thoracic muscles for prolonged breath holding.While technologies, such as DIBH, have advanced to millimeter-precision in treatment delivery and motion tracking, the "human factor" of patients' ability to cooperate and perform has been addressed much less. Both are needed to optimally deliver advanced radiation therapy with minimized normal tissue effects, while alleviating physical and cognitive distress during this challenging phase of breast cancer therapy.This article discusses physical training and psychotherapeutic integrative health approaches, applied to radiation oncology, to leverage and augment the gains enabled by advanced technology-based high-precision radiation treatment in breast cancer. Such combinations of advanced technologies with training and cognitive integrative health interventions hold the promise to provide simple feasible and low-cost means to improve patient experience, emotional outcomes and quality of life, while optimizing patient performance for advanced imaging-guided treatment procedures - paving the way to improve cardiac outcomes in breast cancer survivors.

Combined use of Monte Carlo DNA damage simulations and deterministic repair models to examine putative mechanisms of cell killing.

A kinetic repair-misrepair-fixation (RMF) model is developed to better link double-strand break (DSB) induction to reproductive cell death. Formulas linking linear-quadratic (LQ) model radiosensitivity parameters to DSB induction and repair explicitly account for the contribution to cell killing of unrejoinable DSBs, misrepaired and fixed DSBs, and exchanges formed through intra- and intertrack DSB interactions. Information from Monte Carlo simulations is used to determine the initial yields and complexity of DSBs formed by low- and high-LET radiations. Our analysis of published survival data for human kidney cells suggests that intratrack DSB interactions are negligible for low-LET radiations but increase rapidly with increasing LET. The analysis suggests that no class of DSB is intrinsically unrejoinable or that DSB reparability is not strictly determined by the number of lesions forming the DSB. For radiations with LET >110 keV/mum, the model predicts that the relative cell killing efficiency, per unit absorbed dose, should continue to increase, whereas data from published experiments indicate a reduced cell killing efficiency. This observation suggests that the Monte Carlo simulation overestimates the DSB yield beyond 110 keV/microm or that other biological phenomena not included in the model, such as proximity effects, are important. For 200-250 kVp X rays ( approximately 1.9 keV/microm), only about 1% of the one-track killing is attributed to intratrack binary misrepair interactions. The analysis indicates that the remaining 99% of the lethal damage is due to other types of one-track damage, including possible unrepairable, misrepaired and fixed damage. Compared to the analysis of the X-ray results, 48% of the one-track lethal damage caused by 5.1 MeV alpha particles (approximately 88 keV/microm) is due to intratrack DSB interactions while the remainder is due to other forms of one-track damage.

Dosimetric impact of intrafraction motion for compensator-based proton therapy of lung cancer.

Compensator-based proton therapy of lung cancer using an un-gated treatment while allowing the patient to breathe freely requires a compensator design that ensures tumor coverage throughout respiration. Our investigation had two purposes: one is to investigate the dosimetric impact when a composite compensator correction is applied, or is not, and the other one is to evaluate the significance of using different respiratory phases as the reference computed tomography (CT) for treatment planning dose calculations. A 4D-CT-based phantom study and a real patient treatment planning study were performed. A 3D MIP dataset generated over all phases of the acquired 4D-CT scans was adopted to design the field-specific composite aperture and compensator. In the phantom study, the MIP-based compensator design plan named plan D was compared to the other three plans, in which average intensity projection (AIP) images in conjunction with the composite target volume contour copied from the MIP images were used. Relative electron densities within the target envelope were assigned either to original values from the AIP image dataset (plan A) or to predetermined values, 0.8 (plan B) and 0.9 (plan C). In the patient study, the dosimetric impact of a compensator design based on the MIP images (plan ITV(MIP)) was compared to designs based on end-of-inhale (EOI) (plan ITV(EOI)) and middle-of-exhale (MOE) CT images (plan ITV(MOE)). The dose distributions were recalculated for each phase. Throughout the ten phases, it shows that D(GTV)(min) changed slightly from 86% to 89% (SD = 0.9%) of prescribed dose (PD) in the MIP plan, while varying greatly from 10% to 79% (SD = 26.7%) in plan A, 17% to 73% (SD = 22.5%) in plan B and 53% to 73% (SD = 6.8%) in plan C. The same trend was observed for D(GTV)(mean) and V95 with less amplitude. In the MIP-based plan ITV(MIP), D(GTV)(mean) was almost identically equal to 95% in each phase (SD = 0.5%). The patient study verified that the MIP approach increased the minimum value of D99 of the clinical target volume (CTV) by 58.8% compared to plan ITV(EOI) and 12.9% compared to plan ITV(MOE). Minimum values of D99 were 37.60%, 83.50% and 96.40% for plan ITV(EOI), plan ITV(MOE) and plan ITV(MIP), respectively. Standard deviations of D99 were significantly decreased (SD = 0.5%) in the MIP plan as compared to plan ITV(EOI) (SD = 18.9%) or plan ITV(MOE) (SD = 4.0%). These studies demonstrate that the use of MIP images to design the patient-specific composite compensators provide superior and consistent tumor coverage throughout the entire respiratory cycle whilst maintaining a low average normal lung dose. The additional benefit of the MIP-based design approach is that the dose calculation can be implemented on any single phase as long as it uses the aperture and compensator optimized from the MIP images. This also reduces the requirement for contouring on all breathing phases down to just one.

Dosimetric characteristics of the University of Washington Clinical Neutron Therapy System.

The University of Washington (UW) Clinical Neutron Therapy System (CNTS), which generates high linear energy transfer fast neutrons through interactions of 50.5 MeV protons incident on a Be target, has depth-dose characteristics similar to 6 MV x-rays. In contrast to the fixed beam angles and primitive blocking used in early clinical trials of neutron therapy, the CNTS has a gantry with a full 360° of rotation, internal wedges, and a multi-leaf collimator (MLC). Since October of 1984, over 3178 patients have received conformal neutron therapy treatments using the UW CNTS. In this work, the physical and dosimetric characteristics of the CNTS are documented through comparisons of measurements and Monte Carlo simulations. A high resolution computed tomography scan of the model 17 ionization chamber (IC-17) has also been used to improve the accuracy of simulations of the absolute calibration geometry. The response of the IC-17 approximates well the kinetic energy released per unit mass (KERMA) in water for neutrons and photons for energies from a few tens of keV up to about 20 MeV. Above 20 MeV, the simulated model 17 ion chamber response is 20%-30% higher than the neutron KERMA in water. For CNTS neutrons, simulated on- and off-axis output factors in water match measured values within ~2% ± 2% for rectangular and irregularly shaped field with equivalent square areas ranging in a side dimension from 2.8 cm to 30.7 cm. Wedge factors vary by less than 1.9% of the measured dose in water for clinically relevant field sizes. Simulated tissue maximum ratios in water match measured values within 3.3% at depths up to 20 cm. Although the absorbed dose for water and adipose tissue are within 2% at a depth of 1.7 cm, the absorbed dose in muscle and bone can be as much as 12 to 40% lower than the absorbed dose in water. The reported studies are significant from a historical perspective and as additional validation of a new tool for patient quality assurance and as an aid in ongoing efforts to clinically implement advanced treatment techniques, such as intensity modulated neutron therapy, at the UW.

Does Neutron Radiation Therapy Potentiate an Immune Response to Merkel Cell Carcinoma?

BACKGROUND: Merkel cell carcinoma (MCC) is a rare and aggressive cutaneous malignancy. In the advanced setting, MCC is often treated with immune checkpoint inhibitors such as anti-PD-1/PD-L1 antibodies. X-ray radiation therapy (XRT) is commonly used for palliation. There is an unmet need for new treatment options in patients progressing on immunotherapy and XRT. We present 2 patients with progressive MCC who were successfully treated with high linear energy transfer neutron radiation therapy (NRT). CLINICAL OBSERVATIONS: Patient A, an 85-year-old white male with chronic lymphocytic leukemia had progressive MCC with multiple tumors on the face despite prior XRT and ongoing treatment with pembrolizumab. The 5 most symptomatic lesions were treated with a short course of NRT (2 × 3 Gy) while continuing pembrolizumab. All irradiated facial lesions demonstrated a complete response 2 weeks after NRT. Remarkably, an additional 4 lesions located outside the NRT fields also completely resolved. Patient B, a 78-year-old white male with no immunosuppressive condition had recurrent MCC in the scalp and bilateral cervical nodes. The painful, ulcerative tumors on his scalp were progressing despite multiple courses of XRT and multiple immunotherapy regimens, including pembrolizumab. He was treated with NRT (16-18 Gy) to the scalp and had a complete response with successful palliation. While his disease subsequently progressed outside the NRT fields, the response to NRT bridged him to receive further investigational immunotherapies, and he remains disease free 3 years later. CONCLUSION: Short courses of high linear energy transfer particle therapy deserve consideration as a promising modality for local tumor control in XRT refractory tumors. The out-of-field response suggests that NRT has potential for synergizing with immunotherapy. While more data are required to identify optimal NRT parameters, the NRT dose that potentiates an antitumor immune response appears to be well below organ-at-risk tolerance.

Scattered neutron dose equivalent from an active scanning proton beam delivery system.

A study of neutron production from a novel active scanning proton beam delivery system at the Midwest Proton Radiotherapy Institute (MPRI) has been performed. The neutron dose equivalent was determined using a neutron rem (roentgen equivalent in man) detector which has an upper energy limit of 10 MeV. Measurement were taken at 0, 45, and 90 degrees from the proton beam central axis and for various proton beam energies (127-208 MeV) and scanned field sizes (25-144 cm2). The maximum neutron dose observed was 0.43 mSv / (proton treatment Gy) at 90 degrees from the beam axis for a beam energy of 208.4 MeV and a scanned field size of 144 cm2. It is still possible to further mitigate this secondary neutron dose during treatment by optimizing parameters within the treatment nozzle and using shielding.

Functional lung avoidance and response-adaptive escalation (FLARE) RT: Multimodality plan dosimetry of a precision radiation oncology strategy.

PURPOSE: Nonsmall cell lung cancer (NSCLC) patient radiation therapy (RT) is planned without consideration of spatial heterogeneity in lung function or tumor response. We assessed the dosimetric and clinical feasibility of functional lung avoidance and response-adaptive escalation (FLARE) RT to reduce dose to [99m Tc]MAA-SPECT/CT perfused lung while redistributing an escalated boost dose within [18 F]FDG-PET/CT-defined biological target volumes (BTV). METHODS: Eight stage IIB-IIIB NSCLC patients underwent FDG-PET/CT and MAA-SPECT/CT treatment planning scans. Perfused lung objectives were derived from scatter/collimator/attenuation-corrected MAA-SPECT uptake relative to ITV-subtracted lung to maintain < 20 Gy mean lung dose (MLD). Prescriptions included 60 Gy to the planning target volume (PTV) and concomitant boost of 74 Gy mean to biological target volumes (BTV = GTV + PET gradient segmentation) scaled to each BTV voxel by relative FDG-PET SUV. Dose-painting-by-numbers prescriptions were integrated into commercial treatment planning systems via uptake threshold discretization. Dose constraints for lung, heart, cord, and esophagus were defined. FLARE RT plans were optimized with volumetric modulated arc therapy (VMAT), proton pencil beam scanning (PBS) with 3%-3 mm robust optimization, and combination of PBS (avoidance) plus VMAT (escalation). The high boost dose region was evaluated within a standardized SUVpeak structure. FLARE RT plans were compared to reference VMAT plans. Linear regression between radiation dose to BTV and normalized FDG PET SUV at every voxel was conducted, from which Pearson linear correlation coefficients and regression slopes were extracted. Spearman rank correlation coefficients were estimated between radiation dose to lung and normalized SPECT uptake. Dosimetric differences between treatment modalities were evaluated by Friedman nonparametric paired test with multiple sampling correction. RESULTS: No unacceptable violations of PTV and normal tissue objectives were observed in 24 FLARE RT plans. Compared to reference VMAT plans, FLARE VMAT plans achieved a higher mean dose to BTV (73.7 Gy 98195. 61.3 Gy), higher mean dose to SUVpeak (89.7 Gy vs. 60.8 Gy), and lower mean dose to highly perfused lung (7.3 Gy vs. 14.9 Gy). These dosimetric gains came at the expense of higher mean heart dose (9.4 Gy vs. 5.8 Gy) and higher maximum cord dose (50.1 Gy vs. 44.6 Gy) relative to the reference VMAT plans. Between FLARE plans, FLARE VMAT achieved higher dose to the SUVpeak ROI than FLARE PBS (89.7 Gy vs. 79.2 Gy, P = 0.01), while FLARE PBS delivered lower dose to lung than FLARE VMAT (11.9 Gy vs. 15.6 Gy, P < 0.001). Voxelwise linear dose redistribution slope between BTV dose and FDG PET uptake was higher in magnitude for FLARE PBS + VMAT (0.36 Gy per %SUVmax ) compared to FLARE VMAT (0.27 Gy per %SUVmax ) or FLARE PBS alone (0.17 Gy per %SUVmax ). CONCLUSIONS: FLARE RT is clinically feasible with VMAT and PBS. A combination of PBS for functional lung avoidance and VMAT for FDG PET dose escalation balanced target and normal tissue objective tradeoffs. These results provide a technical platform for testing of FLARE RT safety and efficacy within a precision radiation oncology trial.

Reconstruction of electron spectra from depth doses with adaptive regularization.

Electron spectral reconstruction of medical accelerators from measured depth doses is a practical method for providing the input initial phase space distribution at the patient surface that is required by Monte Carlo treatment planning systems. The posed inverse problem of spectral reconstruction is ill conditioned and this may lead to nonphysical oscillations in the reconstructed spectra. Use of a variational method of solution with a regularization technique removes the oscillations but tends to smooth the sharp (deltalike) energy peak that is a common feature in electron spectra generated by medical accelerators. Because the sharp peak contains a large percentage of the electrons in the spectrum, an accurate estimate of the peak width, height and position is critical to the success of the technique for spectrum reconstruction with regularization. We propose use of an adaptive regularization term as a special form of the general Tichonov regularization function. The variational method with the adaptive regularization term is applied to reconstruct electron spectra for the 6, 9, and 18 MeV electron beams of a Varian Clinac 2100C accelerator and proves to be a very simple, effective and accurate approach. Results using this variational method with adaptive regularization almost perfectly reconstruct electron spectra from depth dose distributions.

Neutron scattered dose equivalent to a fetus from proton radiotherapy of the mother.

Scattered neutron dose equivalent to a representative point for a fetus is evaluated in an anthropomorphic phantom of the mother undergoing proton radiotherapy. The effect on scattered neutron dose equivalent to the fetus of changing the incident proton beam energy, aperture size, beam location, and air gap between the beam delivery snout and skin was studied for both a small field snout and a large field snout. Measurements of the fetus scattered neutron dose equivalent were made by placing a neutron bubble detector 10 cm below the umbilicus of an anthropomorphic Rando phantom enhanced by a wax bolus to simulate a second trimester pregnancy. The neutron dose equivalent in milliSieverts (mSv) per proton treatment Gray increased with incident proton energy and decreased with aperture size, distance of the fetus representative point from the field edge, and increasing air gap. Neutron dose equivalent to the fetus varied from 0.025 to 0.450 mSv per proton Gray for the small field snout and from 0.097 to 0.871 mSv per proton Gray for the large field snout. There is likely to be no excess risk to the fetus of severe mental retardation for a typical proton treatment of 80 Gray to the mother since the scattered neutron dose to the fetus of 69.7 mSv is well below the lower confidence limit for the threshold of 300 mGy observed for the occurrence of severe mental retardation in prenatally exposed Japanese atomic bomb survivors. However, based on the linear no threshold hypothesis, and this same typical treatment for the mother, the excess risk to the fetus of radiation induced cancer death in the first 10 years of life is 17.4 per 10,000 children.

Dosimetric impact of a CT metal artefact suppression algorithm for proton, electron and photon therapies.

Accurate dose calculation is essential to precision radiation treatment planning and this accuracy depends upon anatomic and tissue electron density information. Modern treatment planning inhomogeneity corrections use x-ray CT images and calibrated scales of tissue CT number to electron density to provide this information. The presence of metal in the volume scanned by an x-ray CT scanner causes metal induced image artefacts that influence CT numbers and thereby introduce errors in the radiation dose distribution calculated. This paper investigates the dosimetric improvement achieved by a previously proposed x-ray CT metal artefact suppression technique when the suppressed images of a patient with bilateral hip prostheses are used in commercial treatment planning systems for proton, electron or photon therapies. For all these beam types, this clinical image and treatment planning study reveals that the target may be severely underdosed if a metal artefact-contaminated image is used for dose calculations instead of the artefact suppressed one. Of the three beam types studied, the metal artefact suppression is most important for proton therapy dose calculations, intermediate for electron therapy and least important for x-ray therapy but still significant. The study of a water phantom having a metal rod simulating a hip prosthesis indicates that CT numbers generated after image processing for metal artefact suppression are accurate and thus dose calculations based on the metal artefact suppressed images will be of high fidelity.

Energy response of an aluminium oxide detector in kilovoltage and megavoltage photon beams: an EGSnrc Monte Carlo simulation study.

A Monte Carlo study of the energy response of an aluminium oxide (Al(2)O(3)) detector in kilovoltage and megavoltage photon beams relative to (60)Co gamma rays has been performed using EGSnrc Monte Carlo simulations. The sensitive volume of the Al(2)O(3) detector was simulated as a disc of diameter 2.85 mm and thickness 1 mm. The phantom material was water and the irradiation depth chosen was 2.0 cm in kilovoltage photon beams and 5.0 cm in megavoltage photon beams. The results show that the energy response of the Al(2)O(3) detector is constant within 3% for photon beam energies in the energy range of (60)Co gamma rays to 25 MV X rays. However, the Al(2)O(3) detector shows an enhanced energy response for kilovoltage photon beams, which in the case of 50 kV X rays is 3.2 times higher than that for (60)Co gamma rays. There is essentially no difference in the energy responses of LiF and Al(2)O(3) detectors irradiated in megavoltage photon beams when these Al(2)O(3) results are compared with literature data for LiF thermoluminescence detectors. However, the Al(2)O(3) detector has a much higher enhanced response compared with LiF detectors in kilovoltage X-ray beams, more than twice as much for the case of 50 kV X rays.

MCNP6 model of the University of Washington clinical neutron therapy system (CNTS).

A MCNP6 dosimetry model is presented for the Clinical Neutron Therapy System (CNTS) at the University of Washington. In the CNTS, fast neutrons are generated by a 50.5 MeV proton beam incident on a 10.5 mm thick Be target. The production, scattering and absorption of neutrons, photons, and other particles are explicitly tracked throughout the key components of the CNTS, including the target, primary collimator, flattening filter, monitor unit ionization chamber, and multi-leaf collimator. Simulations of the open field tissue maximum ratio (TMR), percentage depth dose profiles, and lateral dose profiles in a 40 cm × 40 cm × 40 cm water phantom are in good agreement with ionization chamber measurements. For a nominal 10 × 10 field, the measured and calculated TMR values for depths of 1.5 cm, 5 cm, 10 cm, and 20 cm (compared to the dose at 1.7 cm) are within 0.22%, 2.23%, 4.30%, and 6.27%, respectively. For the three field sizes studied, 2.8 cm × 2.8 cm, 10.4 cm × 10.3 cm, and 28.8 cm × 28.8 cm, a gamma test comparing the measured and simulated percent depth dose curves have pass rates of 96.4%, 100.0%, and 78.6% (depth from 1.5 to 15 cm), respectively, using a 3% or 3 mm agreement criterion. At a representative depth of 10 cm, simulated lateral dose profiles have in-field (⩾ 10% of central axis dose) pass rates of 89.7% (2.8 cm × 2.8 cm), 89.6% (10.4 cm × 10.3 cm), and 100.0% (28.8 cm × 28.8 cm) using a 3% and 3 mm criterion. The MCNP6 model of the CNTS meets the minimum requirements for use as a quality assurance tool for treatment planning and provides useful insights and information to aid in the advancement of fast neutron therapy.