Dosimetry for Quantitative Analysis of the Effects of Low-Dose Ionizing Radiation in Radiation Therapy Patients

Out-of-field scattering from the Integral Quality Monitor® in megavolt photon beams

PURPOSE

To investigate out-of-field scattered doses from the Integral Quality Monitor (IQM) transmission detector in megavoltage photon beams.

MATERIALS AND METHODS

We measured out-of-field point doses for 20 × 20 cm² 6-15 MV photon beams using 10 × 2 cm² Gafchromic EBT3 film strips placed across the surfaces of 5-cm thick water-equivalent RW3 slabs. The films were placed at 10 cm intervals from the central axis (CAX) of each beam, up to 1.0 m away on opposite sides of the CAX. The measurements were conducted at 80 cm and 100 cm source-to-surface distances (SSD) without the IQM, and were repeated with the IQM in the paths of the beams. Measurements were also performed at 90 cm SSD for 20 × 20 and 30 × 30 cm² 15 MV fields. Surface dose profiles were then constructed from the measurements for each beam setup with and without the IQM to examine the differences in scattered dose off-axis. The dose profile for each beam setup was normalised to dose maximum measured on the CAX.

RESULTS

Overall, surface dose profiles acquired with the IQM in the paths of the beams were higher than the corresponding profiles without the IQM. The out-of-field dose increased with increase in photon energy, field size, and shorter SSDs, and decreased with off-axis distance. At 80 cm SSD for the 20 × 20 cm² field, the IQM-induced surface dose ranged from -0.6% ÷ 1.9%, -0.3% ÷ 3.0%, and 0.3% ÷ 6.8% for 6, 10, and 15 MV beams, respectively.

CONCLUSION

The higher surface dose profiles measured with the IQM attached to the linac compared to the profiles without the IQM indicates that the device is acting as an additional source of scattered radiation.

Supplemental computational phantoms to estimate out-of-field absorbed dose in photon radiotherapy

The purpose of this study was to develop a straightforward method of supplementing patient anatomy and estimating out-of-field absorbed dose for a cohort of pediatric radiotherapy patients with limited recorded anatomy. A cohort of nine children, aged 2-14 years, who received 3D conformal radiotherapy for low-grade localized brain tumors (LBTs), were randomly selected for this study. The extent of these patients' computed tomography simulation image sets were cranial only. To approximate their missing anatomy, we supplemented the LBT patients' image sets with computed tomography images of patients in a previous study with larger extents of matched sex, height, and mass and for whom contours of organs at risk for radiogenic cancer had already been delineated. Rigid fusion was performed between the LBT patients' data and that of the supplemental computational phantoms using commercial software and in-house codes. In-field dose was calculated with a clinically commissioned treatment planning system, and out-of-field dose was estimated with a previously developed analytical model that was re-fit with parameters based on new measurements for intracranial radiotherapy. Mean doses greater than 1 Gy were found in the red bone marrow, remainder, thyroid, and skin of the patients in this study. Mean organ doses between 150 mGy and 1 Gy were observed in the breast tissue of the girls and lungs of all patients. Distant organs, i.e. prostate, bladder, uterus, and colon, received mean organ doses less than 150 mGy. The mean organ doses of the younger, smaller LBT patients (0-4 years old) were a factor of 2.4 greater than those of the older, larger patients (8-12 years old). Our findings demonstrated the feasibility of a straightforward method of applying supplemental computational phantoms and dose-calculation models to estimate absorbed dose for a set of children of various ages who received radiotherapy and for whom anatomies were largely missing in their original computed tomography simulations.

Data integration reveals key homeostatic mechanisms following low dose radiation exposure

The goal of this study was to define pathways regulated by low dose radiation to understand how biological systems respond to subtle perturbations in their environment and prioritize pathways for human health assessment. Using an in vitro 3-D human full thickness skin model, we have examined the temporal response of dermal and epidermal layers to 10 cGy X-ray using transcriptomic, proteomic, phosphoproteomic and metabolomic platforms. Bioinformatics analysis of each dataset independently revealed potential signaling mechanisms affected by low dose radiation, and integrating data shed additional insight into the mechanisms regulating low dose responses in human tissue. We examined direct interactions among datasets (top down approach) and defined several hubs as significant regulators, including transcription factors (YY1, MYC and CREB1), kinases (CDK2, PLK1) and a protease (MMP2). These data indicate a shift in response across time - with an increase in DNA repair, tissue remodeling and repression of cell proliferation acutely (24-72h). Pathway-based integration (bottom up approach) identified common molecular and pathway responses to low dose radiation, including oxidative stress, nitric oxide signaling and transcriptional regulation through the SP1 factor that would not have been identified by the individual data sets. Significant regulation of key downstream metabolites of nitrative stress was measured within these pathways. Among the features identified in our study, the regulation of MMP2 and SP1 was experimentally validated. Our results demonstrate the advantage of data integration to broadly define the pathways and networks that represent the mechanisms by which complex biological systems respond to perturbation.

Cell type-dependent gene transcription profile in a three-dimensional human skin tissue model exposed to low doses of ionizing radiation: implications for medical exposures

The concern over possible health risks from exposures to low doses of ionizing radiation has been driven largely by the increase in medical exposures, the routine implementation of X-ray backscatter devices for airport security screening, and, most recently, the nuclear incident in Japan. Because of a paucity of direct epidemiological data at very low doses, cancer risk must be estimated from high dose exposure scenarios. However, there is increasing evidence that low and high dose exposures result in different signaling events and may have different response mechanisms than higher doses. We have examined the radiation-induced temporal response after exposure to 10 cGy of an in vitro three dimensional (3D) human skin tissue model using microarray-based transcriptional profiling. Cell type-specific analysis showed significant changes in gene expression with the levels of >1,400 genes altered in the dermis and >400 genes regulated in the epidermis. The two cell layers rarely exhibited overlapping responses at the mRNA level. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) measurements validated the microarray data in both regulation direction and value. Key pathways identified relate to cell cycle regulation, immune responses, hypoxia, reactive oxygen signaling, and DNA damage repair. The proliferation status as well as the expression of PCNA was examined in histological samples. We discuss in particular the role of proliferation, emphasizing how the disregulation of cellular signaling in normal tissue may impact progression toward radiation-induced secondary diseases.

Monte Carlo modeling of a 6 and 18 MV Varian Clinac medical accelerator for in-field and out-of-field dose calculations: development and validation

There is a serious and growing concern about the increased risk of radiation-induced second cancers and late tissue injuries associated with radiation treatment. To better understand and to more accurately quantify non-target organ doses due to scatter and leakage radiation from medical accelerators, a detailed Monte Carlo model of the medical linear accelerator is needed. This paper describes the development and validation of a detailed accelerator model of the Varian Clinac operating at 6 and 18 MV beam energies. Over 100 accelerator components have been defined and integrated using the Monte Carlo code MCNPX. A series of in-field and out-of-field dose validation studies were performed. In-field dose distributions calculated using the accelerator models were tuned to match measurement data that are considered the de facto 'gold standard' for the Varian Clinac accelerator provided by the manufacturer. Field sizes of 4 cm x 4 cm, 10 cm x 10 cm, 20 cm x 20 cm and 40 cm x 40 cm were considered. The local difference between calculated and measured dose on the percent depth dose curve was less than 2% for all locations. The local difference between calculated and measured dose on the dose profile curve was less than 2% in the plateau region and less than 2 mm in the penumbra region for all locations. Out-of-field dose profiles were calculated and compared to measurement data for both beam energies for field sizes of 4 cm x 4 cm, 10 cm x 10 cm and 20 cm x 20 cm. For all field sizes considered in this study, the average local difference between calculated and measured dose for the 6 and 18 MV beams was 14 and 16%, respectively. In addition, a method for determining neutron contamination in the 18 MV operating model was validated by comparing calculated in-air neutron fluence with reported calculations and measurements. The average difference between calculated and measured neutron fluence was 20%. As one of the most detailed accelerator models for both in-field and out-of-field dose calculations, the model will be combined with anatomically realistic computational patient phantoms into a computational framework to calculate non-target organ doses to patients from various radiation treatment plans.

Assessing probe-specific dye and slide biases in two-color microarray data

BACKGROUND

A primary reason for using two-color microarrays is that the use of two samples labeled with different dyes on the same slide, that bind to probes on the same spot, is supposed to adjust for many factors that introduce noise and errors into the analysis. Most users assume that any differences between the dyes can be adjusted out by standard methods of normalization, so that measures such as log ratios on the same slide are reliable measures of comparative expression. However, even after the normalization, there are still probe specific dye and slide variation among the data. We define a method to quantify the amount of the dye-by-probe and slide-by-probe interaction. This serves as a diagnostic, both visual and numeric, of the existence of probe-specific dye bias. We show how this improved the performance of two-color array analysis for arrays for genomic analysis of biological samples ranging from rice to human tissue.

RESULTS

We develop a procedure for quantifying the extent of probe-specific dye and slide bias in two-color microarrays. The primary output is a graphical diagnostic of the extent of the bias which called ECDF (Empirical Cumulative Distribution Function), though numerical results are also obtained.

CONCLUSION

We show that the dye and slide biases were high for human and rice genomic arrays in two gene expression facilities, even after the standard intensity-based normalization, and describe how this diagnostic allowed the problems causing the probe-specific bias to be addressed, and resulted in important improvements in performance. The R package LMGene which contains the method described in this paper has been available to download from Bioconductor.

A Monte Carlo model for out-of-field dose calculation from high-energy photon therapy

As cancer therapy becomes more efficacious and patients survive longer, the potential for late effects increases, including effects induced by radiation dose delivered away from the treatment site. This out-of-field radiation is of particular concern with high-energy radiotherapy, as neutrons are produced in the accelerator head. We recently developed an accurate Monte Carlo model of a Varian 2100 accelerator using MCNPX for calculating the dose away from the treatment field resulting from low-energy therapy. In this study, we expanded and validated our Monte Carlo model for high-energy (18 MV) photon therapy, including both photons and neutrons. Simulated out-of-field photon doses were compared with measurements made with thermoluminescent dosimeters in an acrylic phantom up to 55 cm from the central axis. Simulated neutron fluences and energy spectra were compared with measurements using moderated gold foil activation in moderators and data from the literature. The average local difference between the calculated and measured photon dose was 17%, including doses as low as 0.01% of the central axis dose. The out-of-field photon dose varied substantially with field size and distance from the edge of the field but varied little with depth in the phantom, except at depths shallower than 3 cm, where the dose sharply increased. On average, the difference between the simulated and measured neutron fluences was 19% and good agreement was observed with the neutron spectra. The neutron dose equivalent varied little with field size or distance from the central axis but decreased with depth in the phantom. Neutrons were the dominant component of the out-of-field dose equivalent for shallow depths and large distances from the edge of the treatment field. This Monte Carlo model is useful to both physicists and clinicians when evaluating out-of-field doses and associated potential risks.

Transient genome-wide transcriptional response to low-dose ionizing radiation in vivo in humans

PURPOSE

The in vivo effects of low-dose low linear energy transfer ionizing radiation on healthy human skin are largely unknown. Using a patient-based tissue acquisition protocol, we have performed a series of genomic analyses on the temporal dynamics over a 24-hour period to determine the radiation response after a single exposure of 10 cGy.

METHODS AND MATERIALS

RNA from each patient tissue sample was hybridized to an Affymetrix Human Genome U133 Plus 2.0 array. Data analysis was performed on selected gene groups and pathways.

RESULTS

Nineteen gene groups and seven gene pathways that had been shown to be radiation responsive were analyzed. Of these, nine gene groups showed significant transient transcriptional changes in the human tissue samples, which returned to baseline by 24 hours postexposure.

CONCLUSIONS

Low doses of ionizing radiation on full-thickness human skin produce a definable temporal response out to 24 hours postexposure. Genes involved in DNA and tissue remodeling, cell cycle transition, and inflammation show statistically significant changes in expression, despite variability between patients. These data serve as a reference for the temporal dynamics of ionizing radiation response following low-dose exposure in healthy full-thickness human skin.

A Monte Carlo model for calculating out-of-field dose from a varian 6 MV beam

Dose to the patient outside of the treatment field is important when evaluating the outcome of radiotherapy treatments. However, determining out-of-field doses for any particular treatment plan currently requires either time-consuming measurements or calculated estimations that may be highly uncertain. A Monte Carlo model may allow these doses to be determined quickly, accurately, and with a great degree of flexibility. MCNPX was used to create a Monte Carlo model of a Varian Clinac 2100 accelerator head operated at 6 MV. Simulations of the dose out-of-field were made and measurements were taken with thermoluminescent dosimeters in an acrylic phantom and with an ion chamber in a water tank to validate the Monte Carlo model. Although local differences between the out-of-field doses calculated by the model and those measured did exceed 50% at some points far from the treatment field, the average local difference was only 16%. This included a range of doses as low as 0.01% of the central axis dose, and at distances in excess of 50 cm from the central axis of the treatment field. The out-of-field dose was found to vary with field size and distance from the central axis, but was almost independent of the depth in the phantom except where the dose increased substantially at depths less than dmax. The relationship between dose and kerma was also investigated, and kerma was found to be a good estimate of dose (within 3% on average) except near the surface and in the field penumbra. Our Monte Carlo model was found to well represent typical Varian 2100 accelerators operated at 6 MV.

Human in vivo dose-response to controlled, low-dose low linear energy transfer ionizing radiation exposure

PURPOSE

The effect of low doses of low-linear energy transfer (photon) ionizing radiation (LDIR, <10 cGy) on human tissue when exposure is under normal physiologic conditions is of significant interest to the medical and scientific community in therapeutic and other contexts. Although, to date, there has been no direct assessment of the response of human tissue to LDIR when exposure is under normal physiologic conditions of intact three-dimensional architecture, vasculature, and cell-cell contacts (between epithelial cells and between epithelial and stromal cells).

EXPERIMENTAL DESIGN

In this article, we present the first data on the response of human tissue exposed in vivo to LDIR with precisely controlled and calibrated doses. We evaluated transcriptomic responses to a single exposure of LDIR in the normal skin of men undergoing therapeutic radiation for prostate cancer (research protocol, Health Insurance Portability and Accountability Act-compliant, Institutional Review Board-approved). Using newly developed biostatistical tools that account for individual splice variants and the expected variability of temporal response between humans even when the outcome is measured at a single time, we show a dose-response pattern in gene expression in a number of pathways and gene groups that are biologically plausible responses to LDIR.

RESULTS

Examining genes and pathways identified as radiation-responsive in cell culture models, we found seven gene groups and five pathways that were altered in men in this experiment. These included the Akt/phosphoinositide-3-kinase pathway, the growth factor pathway, the stress/apoptosis pathway, and the pathway initiated by transforming growth factor-beta signaling, whereas gene groups with altered expression included the keratins, the zinc finger proteins and signaling molecules in the mitogen-activated protein kinase gene group. We show that there is considerable individual variability in radiation response that makes the detection of effects difficult, but still feasible when analyzed according to gene group and pathway.

CONCLUSIONS

These results show for the first time that low doses of radiation have an identifiable biosignature in human tissue, irradiated in vivo with normal intact three-dimensional architecture, vascular supply, and innervation. The genes and pathways show that the tissue (a) does detect the injury, (b) initiates a stress/inflammatory response, (c) undergoes DNA remodeling, as suggested by the significant increase in zinc finger protein gene expression, and (d) initiates a "pro-survival" response. The ability to detect a distinct radiation response pattern following LDIR exposure has important implications for risk assessment in both therapeutic and national defense contexts.