A 3-D CT assisted Monte Carlo evaluation of intracavitary brachytherapy implants

Intracavitary brachytherapy (ICB) combined with external beam irradiation for treatment of cervical cancer is highly successful in achieving local control. The M.D. Anderson Cancer Center employs Fletcher Suit Delclos (FSD) applicators. FSD applicators contain shields to limit dose to critical structures. Dosimetric evaluation of ICB implants is limited to assessing dose at reference points. These points serve as surrogates for treatment intensity and critical structure dose. Several studies have mentioned that the ICRU38 reference points inadequately characterize the dose distribution. Also, the ovoid shields are rarely considered in dosimetry. ^ The goal of this dissertation was to ascertain the influence of the ovoid shields on patient dose distributions. Monte Carlo dosimetry (MCD) was applied to patient computed tomography(CT) scans. These data were analyzed to determine the effect of the shields on dose to standard reference points and the bladder and rectum. The hypothesis of this work is that the ICRU38 bladder and rectal points computed conventionally are not clinically acceptable surrogates for the maximum dose points as determined by MCD. ^ MCD was applied to the tandem and ovoids. The FSD ovoids and tandem were modeled in a single input file that allowed dose to be calculated for any patient. Dose difference surface histograms(DDSH) were computed for the bladder and rectum. Reference point doses were compared between shielded and unshielded ovoids, and a commercial treatment planning system. ^ The results of this work showed the tandem tip screw caused a 33% reduction in dose. The ovoid shields reduced the dose by a maximum of 48.9%. DDSHs revealed on average 5% of the bladder surface area was spared 53 cGy and 5% of the rectal surface area was spared 195 cGy. The ovoid shields on average reduced the dose by 18% for the bladder point and 25% for the rectal point. The Student's t-test revealed the ICRU38 bladder and rectal points do not predict the maximum dose for these organs. ^ It is concluded that modeling the tandem and ovoid internal structures is necessary for accurate dose calculations, the bladder shielding segments may not be necessary, and that the ICRU38 bladder point is irrelevant. ^

A comprehensive assessment of environmental ultraviolet radiation induced DNA damage in marine microorganisms

There is evidence that ultraviolet radiation (UVR) is increasing over certain locations on the Earth's surface. Of primary concern is the annual pattern of ozone depletion over Antarctica and the Southern Ocean. Reduction of ozone concentration selectively limits absorption of solar UV-B (290–320 nm), resulting in higher irradiance at the Earth's surface. The effects of ozone depletion on the human population and natural ecosystems, particularly the marine environment, are a matter of considerable concern. Indeed, marine plankton may serve as sensitive indicators of ozone depletion and UV-B fluctuations. Direct biological effects of UVR result from absorption of UV-B by DNA. Once absorbed, energy is dissipated by a variety of pathways, including covalent chemical reactions leading to the formation of photoproducts. The major types of photoproduct formed are cyclobutyl pyrimidine dimer (CPD) and pyrimidine(6-4)pyrimidone dimer [(6-4)PD]. Marine plankton repair these photoproducts using light-dependent photoenzymatic repair or nucleotide excision repair. The studies here show that fluctuations in CPD concentrations in the marine environment at Palmer Station, Antarctica correlate well with ozone concentration and UV-B irradiance at the Earth's surface. A comparison of photoproduct levels in marine plankton and DNA dosimeters show that bacterioplankton display higher resistance to solar UVR than phytoplankton in an ozone depleted environment. DNA damage in marine microorganisms was investigated during two separate latitudinal transects which covered a total range of 140°. We observed the same pattern of change in DNA damage levels in dosimeters and marine plankton as measured using two distinct quantitative techniques. Results from the transects show that differences in photosensitivity exist in marine plankton collected under varying UVR environments. Laboratory studies of Antarctic bacterial isolates confirm that marine bacterioplankton possess differences in survival, DNA damage induction, and repair following exposure to UVR. Results from DNA damage measurements during ozone season, along a latitudinal gradient, and in marine bacterial isolates suggest that changes in environmental UVR correlate with changes in UV-B induced DNA damage in marine microorganisms. Differences in the ability to tolerate UVR stress under different environmental conditions may determine the composition of the microbial communities inhabiting those environments. ^

Monte Carlo and analytical dose calculations for ocular proton therapy

Uveal melanoma is a rare but life-threatening form of ocular cancer. Contemporary treatment techniques include proton therapy, which enables conservation of the eye and its useful vision. Dose to the proximal structures is widely believed to play a role in treatment side effects, therefore, reliable dose estimates are required for properly evaluating the therapeutic value and complication risk of treatment plans. Unfortunately, current simplistic dose calculation algorithms can result in errors of up to 30% in the proximal region. In addition, they lack predictive methods for absolute dose per monitor unit (D/MU) values. ^ To facilitate more accurate dose predictions, a Monte Carlo model of an ocular proton nozzle was created and benchmarked against measured dose profiles to within ±3% or ±0.5 mm and D/MU values to within ±3%. The benchmarked Monte Carlo model was used to develop and validate a new broad beam dose algorithm that included the influence of edgescattered protons on the cross-field intensity profile, the effect of energy straggling in the distal portion of poly-energetic beams, and the proton fluence loss as a function of residual range. Generally, the analytical algorithm predicted relative dose distributions that were within ±3% or ±0.5 mm and absolute D/MU values that were within ±3% of Monte Carlo calculations. Slightly larger dose differences were observed at depths less than 7 mm, an effect attributed to the dose contributions of edge-scattered protons. Additional comparisons of Monte Carlo and broad beam dose predictions were made in a detailed eye model developed in this work, with generally similar findings. ^ Monte Carlo was shown to be an excellent predictor of the measured dose profiles and D/MU values and a valuable tool for developing and validating a broad beam dose algorithm for ocular proton therapy. The more detailed physics modeling by the Monte Carlo and broad beam dose algorithms represent an improvement in the accuracy of relative dose predictions over current techniques, and they provide absolute dose predictions. It is anticipated these improvements can be used to develop treatment strategies that reduce the incidence or severity of treatment complications by sparing normal tissue. ^

Electron intensity modulation for mixed-beam radiation therapy with an x-ray multi-leaf collimator

The current standard treatment for head and neck cancer at our institution uses intensity-modulated x-ray therapy (IMRT), which improves target coverage and sparing of critical structures by delivering complex fluence patterns from a variety of beam directions to conform dose distributions to the shape of the target volume. The standard treatment for breast patients is field-in-field forward-planned IMRT, with initial tangential fields and additional reduced-weight tangents with blocking to minimize hot spots. For these treatment sites, the addition of electrons has the potential of improving target coverage and sparing of critical structures due to rapid dose falloff with depth and reduced exit dose. In this work, the use of mixed-beam therapy (MBT), i.e., combined intensity-modulated electron and x-ray beams using the x-ray multi-leaf collimator (MLC), was explored. The hypothesis of this study was that addition of intensity-modulated electron beams to existing clinical IMRT plans would produce MBT plans that were superior to the original IMRT plans for at least 50% of selected head and neck and 50% of breast cases. Dose calculations for electron beams collimated by the MLC were performed with Monte Carlo methods. An automation system was created to facilitate communication between the dose calculation engine and the treatment planning system. Energy and intensity modulation of the electron beams was accomplished by dividing the electron beams into 2x2-cm2 beamlets, which were then beam-weight optimized along with intensity-modulated x-ray beams. Treatment plans were optimized to obtain equivalent target dose coverage, and then compared with the original treatment plans. MBT treatment plans were evaluated by participating physicians with respect to target coverage, normal structure dose, and overall plan quality in comparison with original clinical plans. The physician evaluations did not support the hypothesis for either site, with MBT selected as superior in 1 out of the 15 head and neck cases (p=1) and 6 out of 18 breast cases (p=0.95). While MBT was not shown to be superior to IMRT, reductions were observed in doses to critical structures distal to the target along the electron beam direction and to non-target tissues, at the expense of target coverage and dose homogeneity. ^

Quantitative comparison of late effects following photon versus proton external-beam radiation therapies: Toward an evidence-based approach to selecting a treatment modality

Radiation therapy has been used as an effective treatment for malignancies in pediatric patients. However, in many cases, the side effects of radiation diminish these patients’ quality of life. In order to develop strategies to minimize radiogenic complications, one must first quantitatively estimate pediatric patients’ relative risk for radiogenic late effects, which has not become feasible till recently because of the calculational complexity. The goals of this work were to calculate the dose delivered to tissues and organs in pediatric patients during contemporary photon and proton radiotherapies; to estimate the corresponding risk of radiogenic second cancer and cardiac toxicity based on the calculated doses and on dose-risk models from the literature; to test for the statistical significance of the difference between predicted risks after photon versus proton radiotherapies; and to provide a prototype of an evidence-based approach to selecting treatment modalities for pediatric patients, taking second cancer and cardiac toxicity into account. The results showed that proton therapy confers a lower predicted risk of radiogenic second cancer, and lower risks of radiogenic cardiac toxicities, compared to photon therapy. An uncertainty analysis revealed that the qualitative findings of this study are insensitive to changes in a wide variety of host and treatment related factors.

Information bias and lifetime mortality risks of radiation-induced cancer

Additive and multiplicative models of relative risk were used to measure the effect of cancer misclassification and DS86 random errors on lifetime risk projections in the Life Span Study (LSS) of Hiroshima and Nagasaki atomic bomb survivors. The true number of cancer deaths in each stratum of the cancer mortality cross-classification was estimated using sufficient statistics from the EM algorithm. Average survivor doses in the strata were corrected for DS86 random error ($\sigma$ = 0.45) by use of reduction factors. Poisson regression was used to model the corrected and uncorrected mortality rates with covariates for age at-time-of-bombing, age at-time-of-death and gender. Excess risks were in good agreement with risks in RERF Report 11 (Part 2) and the BEIR-V report. Bias due to DS86 random error typically ranged from $-$15% to $-$30% for both sexes, and all sites and models. The total bias, including diagnostic misclassification, of excess risk of nonleukemia for exposure to 1 Sv from age 18 to 65 under the non-constant relative projection model was $-$37.1% for males and $-$23.3% for females. Total excess risks of leukemia under the relative projection model were biased $-$27.1% for males and $-$43.4% for females. Thus, nonleukemia risks for 1 Sv from ages 18 to 85 (DRREF = 2) increased from 1.91%/Sv to 2.68%/Sv among males and from 3.23%/Sv to 4.02%/Sv among females. Leukemia excess risks increased from 0.87%/Sv to 1.10%/Sv among males and from 0.73%/Sv to 1.04%/Sv among females. Bias was dependent on the gender, site, correction method, exposure profile and projection model considered. Future studies that use LSS data for U.S. nuclear workers may be downwardly biased if lifetime risk projections are not adjusted for random and systematic errors. (Supported by U.S. NRC Grant NRC-04-091-02.) ^

A nested case-control study of radiation exposure and brain cancer incidence in the United States Air Force

A nested case-control study design was used to investigate the relationship between radiation exposure and brain cancer risk in the United States Air Force (USAF). The cohort consisted of approximately 880,000 men with at least 1 year of service between 1970 and 1989. Two hundred and thirty cases were identified from hospital discharge records with a diagnosis of primary malignant brain tumor (International Classification of Diseases, 9th revision, code 191). Four controls were exactly matched with each case on year of age and race using incidence density sampling. Potential career summary extremely low frequency (ELF) and microwave-radiofrequency (MWRF) radiation exposures were based upon the duration in each occupation and an intensity score assigned by an expert panel. Ionizing radiation (IR) exposures were obtained from personal dosimetry records.^ Relative to the unexposed, the overall age-race adjusted odds ratio (OR) for ELF exposure was 1.39, 95 percent confidence interval (CI) 1.03-1.88. A dose-response was not evident. The same was true for MWRF, although the OR = 1.59, with 95 percent CI 1.18-2.16. Excess risk was not found for IR exposure (OR = 0.66, 45 percent CI 0.26-1.72).^ Increasing socioeconomic status (SES), as identified by military pay grade, was associated with elevated brain tumor risk (officer vs. enlisted personnel age-race adjusted OR = 2.11, 95 percent CI 1.98-3.01, and senior officers vs. all others age-race adjusted OR = 3.30, 95 percent CI 2.0-5.46). SES proved to be an important confounder of the brain tumor risk associated with ELF and MWRF exposure. For ELF, the age-race-SES adjusted OR = 1.28, 95 percent CI 0.94-1.74, and for MWRF, the age-race-SES adjusted OR = 1.39, 95 percent CI 1.01-1.90.^ These results indicate that employment in Air Force occupations with potential electromagnetic field exposures is weakly, though not significantly, associated with increased risk for brain tumors. SES appeared to be the most consistent brain tumor risk factor in the USAF cohort. Other investigators have suggested that an association between brain tumor risk and SES may arise from differential access to medical care. However, in the USAF cohort health care is universally available. This study suggests that some factor other than access to medical care must underlie the association between SES and brain tumor risk. ^

AN INVESTIGATION OF THE MODIFYING EFFECTS OF VITAMIN A AND HYPOXIC CELL SENSITIZERS IN RADIATION CARCINOGENESIS IN MICE

The effect of vitamin A (retinyl acetate) and three hypoxic cell sensitizers (metronidazole, misonidazole and desmethylmisonidazole) on lung tumor development in strain A mice exposed to radiation was assessed.^ In experiments involving vitamin A, two groups of mice were fed a low vitamin A diet (< 100 IU/100g diet) while the two other groups were fed a high vitamin A diet (800 IU/100g diet). After two weeks one group maintained on the high vitamin A diet and one group maintained on the low vitamin A diet were given an acute dose of 500 rad of gamma radiation to the thoracic region. The circulating level of plasma vitamin A in all four groups of mice was monitored. A difference in circulating vitamin A in the mice maintained on high and low vitamin A diet became evident by 20 weeks and continued for the duration of the experiment. Mice were killed 18, 26, and 40 weeks post irradiation, their lungs were removed and the number of surface adenomas were counted. There was a significant increase in the number of mice bearing lung tumors and the mean number of lung tumors per mouse in the irradiated group maintained on the high vitamin A diet at 40 weeks post irradiation as compared to the irradiated group maintained on a low vitamin A diet (p < 0.05). Under the conditions of this experiment the development of pulmonary adenomas in irradiated strain A mice appears to relate directly to circulating levels of vitamin A.^ In the other experiment two dose levels of the hypoxic cell sensitizers, 0.2mg/g and 0.6mg/g, were used either alone or in combination with 900 rad of gamma radiation in a fractionated dose schedule of twice a week for three weeks. In the groups of mice which received hypoxic cell sensitizers only, the prevalence and the mean number of lung tumors per mouse were somewhat increased (p < 0.10) in the higher dose group (0.6mg/g) of misonidazole but was not significantly different from the control animals in the other two sensitizer groups. The combination of hypoxic cell sensitizer and radiation did not show any significant enhancement of lung tumor response when compared with the group which received radiation only. The dose of radiation used in this study significantly enhanced lung tumor formation in mice when compared with the control group. Thus, under the experimental exposure conditions used in this investigation, which were very similar to the exposure conditions occurring in clinical treatment, all three hypoxic cell sensitizers did not sensitize the mouse to the carcinogenic effects of gamma radiation.^

EVALUATION OF THE EFFECTIVENESS OF ANISOTROPIC ANALYTICAL ALGORITHM IN FLATTENED AND FLATTENING-FILTER-FREE BEAMS FOR HIGH ENERGY LUNG DOSE DELIVERY USING THE RADIOLOGICAL PHYSICS CENTER LUNG PHANTOM

The effectiveness of the Anisotropic Analytical Algorithm (AAA) implemented in the Eclipse treatment planning system (TPS) was evaluated using theRadiologicalPhysicsCenteranthropomorphic lung phantom using both flattened and flattening-filter-free high energy beams. Radiation treatment plans were developed following the Radiation Therapy Oncology Group and theRadiologicalPhysicsCenterguidelines for lung treatment using Stereotactic Radiation Body Therapy. The tumor was covered such that at least 95% of Planning Target Volume (PTV) received 100% of the prescribed dose while ensuring that normal tissue constraints were followed as well. Calculated doses were exported from the Eclipse TPS and compared with the experimental data as measured using thermoluminescence detectors (TLD) and radiochromic films that were placed inside the phantom. The results demonstrate that the AAA superposition-convolution algorithm is able to calculate SBRT treatment plans with all clinically used photon beams in the range from 6 MV to 18 MV. The measured dose distribution showed a good agreement with the calculated distribution using clinically acceptable criteria of ±5% dose or 3mm distance to agreement. These results show that in a heterogeneous environment a 3D pencil beam superposition-convolution algorithms with Monte Carlo pre-calculated scatter kernels, such as AAA, are able to reliably calculate dose, accounting for increased lateral scattering due to the loss of electronic equilibrium in low density medium. The data for high energy plans (15 MV and 18 MV) showed very good tumor coverage in contrast to findings by other investigators for less sophisticated dose calculation algorithms, which demonstrated less than expected tumor doses and generally worse tumor coverage for high energy plans compared to 6MV plans. This demonstrates that the modern superposition-convolution AAA algorithm is a significant improvement over previous algorithms and is able to calculate doses accurately for SBRT treatment plans in the highly heterogeneous environment of the thorax for both lower (≤12 MV) and higher (greater than 12 MV) beam energies.

Improving Cervical Cancer Nodal Boost Radiation Therapy by Quantifying Uncertainties and Exploring Advanced Radiation Therapy Modalities

Radiation therapy for patients with intact cervical cancer is frequently delivered using primary external beam radiation therapy (EBRT) followed by two fractions of intracavitary brachytherapy (ICBT). Although the tumor is the primary radiation target, controlling microscopic disease in the lymph nodes is just as critical to patient treatment outcome. In patients where gross lymphadenopathy is discovered, an extra EBRT boost course is delivered between the two ICBT fractions. Since the nodal boost is an addendum to primary EBRT and ICBT, the prescription and delivery must be performed considering previously delivered dose. This project aims to address the major issues of this complex process for the purpose of improving treatment accuracy while increasing dose sparing to the surrounding normal tissues. Because external beam boosts to involved lymph nodes are given prior to the completion of ICBT, assumptions must be made about dose to positive lymph nodes from future implants. The first aim of this project was to quantify differences in nodal dose contribution between independent ICBT fractions. We retrospectively evaluated differences in the ICBT dose contribution to positive pelvic nodes for ten patients who had previously received external beam nodal boost. Our results indicate that the mean dose to the pelvic nodes differed by up to 1.9 Gy between independent ICBT fractions. The second aim is to develop and validate a volumetric method for summing dose of the normal tissues during prescription of nodal boost. The traditional method of dose summation uses the maximum point dose from each modality, which often only represents the worst case scenario. However, the worst case is often an exaggeration when highly conformal therapy methods such as intensity modulated radiation therapy (IMRT) are used. We used deformable image registration algorithms to volumetrically sum dose for the bladder and rectum and created a voxel-by-voxel validation method. The mean error in deformable image registration results of all voxels within the bladder and rectum were 5 and 6 mm, respectively. Finally, the third aim explored the potential use of proton therapy to reduce normal tissue dose. A major physical advantage of protons over photons is that protons stop after delivering dose in the tumor. Although theoretically superior to photons, proton beams are more sensitive to uncertainties caused by interfractional anatomical variations, and must be accounted for during treatment planning to ensure complete target coverage. We have demonstrated a systematic approach to determine population-based anatomical margin requirements for proton therapy. The observed optimal treatment angles for common iliac nodes were 90° (left lateral) and 180° (posterior-anterior [PA]) with additional 0.8 cm and 0.9 cm margins, respectively. For external iliac nodes, lateral and PA beams required additional 0.4 cm and 0.9 cm margins, respectively. Through this project, we have provided radiation oncologists with additional information about potential differences in nodal dose between independent ICBT insertions and volumetric total dose distribution in the bladder and rectum. We have also determined the margins needed for safe delivery of proton therapy when delivering nodal boosts to patients with cervical cancer.

Pencil -beam redefinition algorithm dose calculations for electron therapy treatment planning

The electron pencil-beam redefinition algorithm (PBRA) of Shiu and Hogstrom has been developed for use in radiotherapy treatment planning (RTP). Earlier studies of Boyd and Hogstrom showed that the PBRA lacked an adequate incident beam model, that PBRA might require improved electron physics, and that no data existed which allowed adequate assessment of the PBRA-calculated dose accuracy in a heterogeneous medium such as one presented by patient anatomy. The hypothesis of this research was that by addressing the above issues the PBRA-calculated dose would be accurate to within 4% or 2 mm in regions of high dose gradients. A secondary electron source was added to the PBRA to account for collimation-scattered electrons in the incident beam. Parameters of the dual-source model were determined from a minimal data set to allow ease of beam commissioning. Comparisons with measured data showed 3% or better dose accuracy in water within the field for cases where 4% accuracy was not previously achievable. A measured data set was developed that allowed an evaluation of PBRA in regions distal to localized heterogeneities. Geometries in the data set included irregular surfaces and high- and low-density internal heterogeneities. The data was estimated to have 1% precision and 2% agreement with accurate, benchmarked Monte Carlo (MC) code. PBRA electron transport was enhanced by modeling local pencil beam divergence. This required fundamental changes to the mathematics of electron transport (divPBRA). Evaluation of divPBRA with the measured data set showed marginal improvement in dose accuracy when compared to PBRA; however, 4% or 2mm accuracy was not achieved by either PBRA version for all data points. Finally, PBRA was evaluated clinically by comparing PBRA- and MC-calculated dose distributions using site-specific patient RTP data. Results show PBRA did not agree with MC to within 4% or 2mm in a small fraction (<3%) of the irradiated volume. Although the hypothesis of the research was shown to be false, the minor dose inaccuracies should have little or no impact on RTP decisions or patient outcome. Therefore, given ease of beam commissioning, documentation of accuracy, and calculational speed, the PBRA should be considered a practical tool for clinical use. ^

Ultraviolet radiation modulates the immune system through the platelet -activating factor receptor

Ultraviolet radiation plays a critical role in the induction of non-melanoma skin cancer. UV radiation is also immune suppressive. Moreover, UV-induced systemic immune suppression is a major risk factor for skin cancer induction. Previous work had shown that UV exposure in vivo activates a cytokine cascade involving PGE2, IL-4, and IL-10 that induces immune suppression. However, the earliest molecular events that occur immediately after UV-exposure, especially those upstream of PGE2, were not well defined. To determine the initial events and mediators that lead to immune suppression after a pathological dose of UV, mouse keratinocytes were analyzed after sunlamp irradiation. It is known that UV-irradiated keratinocytes secrete the phospholipid mediator of inflammation, platelet-activating factor (PAF). Since PAF stimulates the production of immunomodulatory compounds, including PGE2, the hypothesis that UV-induced PAF activates cytokine production and initiates UV-induced immune suppression was tested. Both UV and PAF activated the transcription of cyclooxygenase (COX)-2 and IL-10 reporter gene constructs. A PAF receptor antagonist blocked UV-induced IL, 10 and COX-2 transcription. PAF mimicked the effects of UV in vivo and suppressed delayed-type hypersensitivity (DTH), and immune suppression was blocked when UV-irradiated mice were injected with a PAF receptor antagonist. This work shows that UV generates PAF-like oxidized lipids, that signal through the PAF receptor, activate cytokine transcription, and induce systemic immune suppression. ^