Thoracic target volume delineation using various maximum-intensity projection computed tomography image sets for stereotactic body radiation therapy

The motion of lung tumors during respiration makes the accurate delivery of radiation therapy to the thorax difficult because it increases the uncertainty of target position. The adoption of four-dimensional computed tomography (4D-CT) has allowed us to determine how a tumor moves with respiration for each individual patient. Using information acquired during a 4D-CT scan, we can define the target, visualize motion, and calculate dose during the planning phase of the radiotherapy process. One image data set that can be created from the 4D-CT acquisition is the maximum-intensity projection (MIP). The MIP can be used as a starting point to define the volume that encompasses the motion envelope of the moving gross target volume (GTV). Because of the close relationship that exists between the MIP and the final target volume, we investigated four MIP data sets created with different methodologies (3 using various 4D-CT sorting implementations, and one using all available cine CT images) to compare target delineation. It has been observed that changing the 4D-CT sorting method will lead to the selection of a different collection of images; however, the clinical implications of changing the constituent images on the resultant MIP data set are not clear. There has not been a comprehensive study that compares target delineation based on different 4D-CT sorting methodologies in a patient population. We selected a collection of patients who had previously undergone thoracic 4D-CT scans at our institution, and who had lung tumors that moved at least 1 cm. We then generated the four MIP data sets and automatically contoured the target volumes. In doing so, we identified cases in which the MIP generated from a 4D-CT sorting process under-represented the motion envelope of the target volume by more than 10% than when measured on the MIP generated from all of the cine CT images. The 4D-CT methods suffered from duplicate image selection and might not choose maximum extent images. Based on our results, we suggest utilization of a MIP generated from the full cine CT data set to ensure a representative inclusive tumor extent, and to avoid geometric miss.

EVALUATION OF THE QUANTITATIVE ACCURACY OF A COMMERCIALLY-AVAILABLE POSITRON EMISSION MAMMOGRAPHY SCANNER

Objective: The PEM Flex Solo II (Naviscan, Inc., San Diego, CA) is currently the only commercially-available positron emission mammography (PEM) scanner. This scanner does not apply corrections for count rate effects, attenuation or scatter during image reconstruction, potentially affecting the quantitative accuracy of images. This work measures the overall quantitative accuracy of the PEM Flex system, and determines the contributions of error due to count rate effects, attenuation and scatter. Materials and Methods: Gelatin phantoms were designed to simulate breasts of different sizes (4 – 12 cm thick) with varying uniform background activity concentration (0.007 – 0.5 μCi/cc), cysts and lesions (2:1, 5:1, 10:1 lesion-to-background ratios). The overall error was calculated from ROI measurements in the phantoms with a clinically relevant background activity concentration (0.065 μCi/cc). The error due to count rate effects was determined by comparing the overall error at multiple background activity concentrations to the error at 0.007 μCi/cc. A point source and cold gelatin phantoms were used to assess the errors due to attenuation and scatter. The maximum pixel values in gelatin and in air were compared to determine the effect of attenuation. Scatter was evaluated by comparing the sum of all pixel values in gelatin and in air. Results: The overall error in the background was found to be negative in phantoms of all thicknesses, with the exception of the 4-cm thick phantoms (0%±7%), and it increased with thickness (-34%±6% for the 12-cm phantoms). All lesions exhibited large negative error (-22% for the 2:1 lesions in the 4-cm phantom) which increased with thickness and with lesion-to-background ratio (-85% for the 10:1 lesions in the 12-cm phantoms). The error due to count rate in phantoms with 0.065 μCi/cc background was negative (-23%±6% for 4-cm thickness) and decreased with thickness (-7%±7% for 12 cm). Attenuation was a substantial source of negative error and increased with thickness (-51%±10% to -77% ±4% in 4 to 12 cm phantoms, respectively). Scatter contributed a relatively constant amount of positive error (+23%±11%) for all thicknesses. Conclusion: Applying corrections for count rate, attenuation and scatter will be essential for the PEM Flex Solo II to be able to produce quantitatively accurate images.

AN IMPLANTABLE MOSFET DOSIMETER MODIFIED TO ACT AS A FIDUCIAL MARKER

In external beam radiation therapy, it is imperative that the prescribed dose is administered to the correct location and in the correct amount. Though several ex vivo methods of quality assurance are currently employed to achieve this goal, verifying that the correct dose is received within the patient in situ is impossible without the capability of measuring dose inside the patient. Recently, a method of measuring dose delivered within the patient has been developed, an implantable MOSFET dosimeter. This dosimeter is implanted within the patient and records the dose received. Since the dosimeter is implanted in the patient, it could serve a dual function as a fiducial marker for image guided radiation therapy (IGRT) treatment if it could be modified to be visible on x-rays. In this study, modifications to the MOSFET dosimeter were made to increase its visibility for IGRT treatment. To test whether the modifications hindered the dosimeter’s ability to accurately measure and transmit dose information, the energy dependence, angular dependence and wireless read range of the modified dosimeter were measured and compared to unmodified dosimeters. It was found that the modified dosimeter performed as well as the unmodified dosimeter while also being suitable for use as a fiducial marker for IGRT treatment.

Clinical Impact of Couch Top and Rails on IMRT and Arc Therapy

Purpose: To evaluate the clinical impact of the Varian Exact Couch on dose and volume coverage to targets and critical structures and tumor control probability (TCP) for 6-MV IMRT and Arc Therapy. Methods: Five clinical prostate patients were planned with both, 6-MV 8-field IMRT and 6-MV 2-field RapidArc using the Eclipse treatment planning system (TPS). These plans neglected treatment couch attenuation, as is standard clinical practice. Dose distributions were then recalculated in Eclipse with the inclusion of the Varian Exact Couch (imaging couch top) and the rails in varying configurations. The changes in dose and coverage were evaluated using the DVHs from each plan iteration. We used a tumor control probability (TCP) model to calculate losses in tumor control resulting from not accounting for the couch top and rails. We also verified dose measurements in a phantom. Results: Failure to account for the treatment couch and rails resulted in clinically unacceptable dose and volume coverage losses to the target for both IMRT and RapidArc. The couch caused average dose losses (relative to plans that ignored the couch) to the prostate of 4.2% and 2.0% for IMRT with the rails out and in, respectively, and 3.2% and 2.9% for RapidArc with the rails out and in, respectively. On average, the percentage of the target covered by the prescribed dose dropped to 35% and 84% for IMRT (rails out and in, respectively) and to 18% and 17% for RapidArc (rails out and in, respectively). The TCP was also reduced by as much as 10.5% (6.3% on average). Dose and volume coverage losses for IMRT plans were primarily due to the rails, while the imaging couch top contributed most to losses for RapidArc. Both the couch top and rails contribute to dose and coverage losses that can render plans clinically unacceptable. A follow-up study we performed found that the less attenuating unipanel mesh couch top available with the Varian Exact couch does not cause a clinically impactful loss of dose or coverage for IMRT but still causes an unacceptable loss for RapidArc. Conclusions: Both the imaging couch top and rails contribute to dose and coverage loss to a degree that, if included, would prevent the plan from meeting clinical planning criteria. Therefore, the imaging and mesh couch tops and rails should be accounted for in Arc Therapy and the imaging couch and rails only in IMRT treatment planning.

Assessment of Collimator Jaw Optimization in Reducing Normal Tissue Irradiation with Intensity Modulated Radiation Therapy

Purpose: To evaluate normal tissue dose reduction in step-and-shoot intensity-modulated radiation therapy (IMRT) on the Varian 2100 platform by tracking the multileaf collimator (MLC) apertures with the accelerator jaws. Methods: Clinical radiation treatment plans for 10 thoracic, 3 pediatric and 3 head and neck patients were converted to plans with the jaws tracking each segment’s MLC apertures. Each segment was then renormalized to account for the change in collimator scatter to obtain target coverage within 1% of that in the original plan. The new plans were compared to the original plans in a commercial radiation treatment planning system (TPS). Reduction in normal tissue dose was evaluated in the new plan by using the parameters V5, V10, and V20 in the cumulative dose-volume histogram for the following structures: total lung minus GTV (gross target volume), heart, esophagus, spinal cord, liver, parotids, and brainstem. In order to validate the accuracy of our beam model, MLC transmission measurements were made and compared to those predicted by the TPS. Results: The greatest change between the original plan and new plan occurred at lower dose levels. The reduction in V20 was never more than 6.3% and was typically less than 1% for all patients. The reduction in V5 was 16.7% maximum and was typically less than 3% for all patients. The variation in normal tissue dose reduction was not predictable, and we found no clear parameters that indicated which patients would benefit most from jaw tracking. Our TPS model of MLC transmission agreed with measurements with absolute transmission differences of less than 0.1 % and thus uncertainties in the model did not contribute significantly to the uncertainty in the dose determination. Conclusion: The amount of dose reduction achieved by collimating the jaws around each MLC aperture in step-and-shoot IMRT does not appear to be clinically significant.

IN-VIVO CT DOSIMETRY DURING VIRTUAL COLONOSCOPY

Virtual colonoscopy (VC) is a minimally invasive means for identifying colorectal polyps and colorectal lesions by insufflating a patient’s bowel, applying contrast agent via rectal catheter, and performing multi-detector computed tomography (MDCT) scans. The technique is recommended for colonic health screening by the American Cancer Society but not funded by the Centers for Medicare and Medicaid Services (CMS) partially because of potential risks from radiation exposure. To date, no in‐vivo organ dose measurements have been performed for MDCT scans; thus, the accuracy of any current dose estimates is currently unknown. In this study, two TLDs were affixed to the inner lumen of standard rectal catheters used in VC, and in-vivo rectal dose measurements were obtained within 6 VC patients. In order to calculate rectal dose, TLD-100 powder response was characterized at diagnostic doses such that appropriate correction factors could be determined for VC. A third-order polynomial regression with a goodness of fit factor of R2=0.992 was constructed from this data. Rectal dose measurements were acquired with TLDs during simulated VC within a modified anthropomorphic phantom configured to represent three sizes of patients undergoing VC. The measured rectal doses decreased in an exponential manner with increasing phantom effective diameter, with R2=0.993 for the exponential regression model and a maximum percent coefficient of variation (%CoV) of 4.33%. In-vivo measurements yielded rectal doses ranged from that decreased exponentially with increasing patient effective diameter, in a manner that was also favorably predicted by the size specific dose estimate (SSDE) model for all VC patients that were of similar age, body composition, and TLD placement. The measured rectal dose within a younger patient was favorably predicted by the anthropomorphic phantom dose regression model due to similarities in the percentages of highly attenuating material at the respective measurement locations and in the placement of the TLDs. The in-vivo TLD response did not increase in %CoV with decreasing dose, and the largest %CoV was 10.0%.

THE MECHANISM OF TUMORIGENESIS IN THE IMMORTALIZED HUMAN PANCREATIC CELL LINES: CELL CULTURE MODELS OF HUMAN PANCREATIC CANCER

The mechanism of tumorigenesis in the immortalized human pancreatic cell lines: cell culture models of human pancreatic cancer Pancreatic ductal adenocarcinoma (PDAC) is the most lethal cancer in the world. The most common genetic lesions identified in PDAC include activation of K-ras (90%) and Her2 (70%), loss of p16 (95%) and p14 (40%), inactivation p53 (50-75%) and Smad4 (55%). However, the role of these signature gene alterations in PDAC is still not well understood, especially, how these genetic lesions individually or in combination contribute mechanistically to human pancreatic oncogenesis is still elusive. Moreover, a cell culture transformation model with sequential accumulation of signature genetic alterations in human pancreatic ductal cells that resembles the multiple-step human pancreatic carcinogenesis is still not established. In the present study, through the stepwise introduction of the signature genetic alterations in PDAC into the HPV16-E6E7 immortalized human pancreatic duct epithelial (HPDE) cell line and the hTERT immortalized human pancreatic ductal HPNE cell line, we developed the novel experimental cell culture transformation models with the most frequent gene alterations in PDAC and further dissected the molecular mechanism of transformation. We demonstrated that the combination of activation of K-ras and Her2, inactivation of p16/p14 and Smad4, or K-ras mutation plus p16 inactivation, was sufficient for the tumorigenic transformation of HPDE or HPNE cells respectively. We found that these transformed cells exhibited enhanced cell proliferation, anchorage-independent growth in soft agar, and grew tumors with PDAC histopathological features in orthotopic mouse model. Molecular analysis showed that the activation of K-ras and Her2 downstream effector pathways –MAPK, RalA, FAK, together with upregulation of cyclins and c-myc were involved in the malignant transformation. We discovered that MDM2, BMP7 and Bmi-1 were overexpressed in the tumorigenic HPDE cells, and that Smad4 played important roles in regulation of BMP7 and Bmi-1 gene expression and the tumorigenic transformation of HPDE cells. IPA signaling pathway analysis of microarray data revealed that abnormal signaling pathways are involved in transformation. This study is the first complete transformation model of human pancreatic ductal cells with the most common gene alterations in PDAC. Altogether, these novel transformation models more closely recapitulate the human pancreatic carcinogenesis from the cell origin, gene lesion, and activation of specific signaling pathway and histopathological features.

DEVELOPMENT AND IMPLEMENTATION OF A REMOTE AUDIT TOOL FOR HIGH DOSE RATE (HDR) 192IR BRACHYTHERAPY USING OPTICALLY STIMULATED LUMINESCENCE DOSIMETRY

This work aimed to create a mailable and OSLD-based phantom with accuracy suitable for RPC audits of HDR brachytherapy sources at institutions participating in NCI-funded cooperative clinical trials. An 8 × 8 × 10 cm3 prototype with two slots capable of holding nanoDot Al2O3:C OSL dosimeters (Landauer, Glenwood, IL) was designed and built. The phantom has a single channel capable of accepting all 192Ir HDR brachytherapy sources in current clinical use in the United States. Irradiations were performed with an 192Ir HDR source to determine correction factors for linearity with dose, dose rate, and the combined effect of irradiation energy and phantom construction. The uncertainties introduced by source positioning in the phantom and timer resolution limitations were also investigated. It was found that the linearity correction factor was where dose is in cGy, which differed from that determined by the RPC for the same batch of dosimeters under 60Co irradiation. There was no significant dose rate effect. Separate energy+block correction factors were determined for both models of 192Ir sources currently in clinical use and these vendor-specific correction factors differed by almost 2.6%. For Nucletron sources, this correction factor was 1.026±0.004 (99% Confidence Interval) and for Varian sources it was 1.000±0.007 (99% CI). Reasonable deviations in source positioning within the phantom and the limited resolution of the source timer had insignificant effects on the ability to measure dose. Overall measurement uncertainty of the system was estimated to be ±2.5% for both Nucletron and Varian source audits (95% CI). This uncertainty was sufficient to establish a ±5% acceptance criterion for source strength audits under a formal RPC audit program. Trial audits of eight participating institutions resulted in an average RPC-to-institution dose ratio of 1.000 with a standard deviation of 0.011.

Development and Implementation of the Use of Optically Stimulated Luminescent Detectors in the Radiological Physics Center Anthropomorphic Quality Assurance Phantoms

The Radiological Physics Center (RPC) uses both on-site and remote reviews to credential institutions for participation in clinical trials. Anthropomorphic quality assurance (QA) phantoms are one tool the RPC uses to remotely audit institutions, which include thermoluminescent dosimeters (TLDs) and radiochromic film. The RPC desires to switch from TLD as the absolute dosimeter in the phantoms, to optically stimulated luminescent dosimeters (OSLDs), but a problem lies in the angular dependence exhibited by the OSLD. The purpose of this study was to characterize the angular dependence of OSLD and establish a correction factor if necessary, to provide accurate dosimetric measurements as a replacement for TLD in the QA phantoms. A 10 cm diameter high-impact polystyrene spherical phantom was designed and constructed to hold an OSLD to study the angular response of the dosimeter under the simplest of circumstances for both coplanar and non-coplanar treatment deliveries. OSLD were irradiated in the spherical phantom, and the responses of the dosimeter from edge-on angles were normalized to the response when irradiated with the beam incident normally on the surface of the dosimeter. The average normalized response was used to establish an angular correction factor for 6 MV and 18 coplanar treatments, and for 6 MV non-coplanar treatments specific to CyberKnife. The RPC pelvic phantom dosimetry insert was modified to hold OSLD, in addition to the TLD, adjacent to the planes of film. Treatment plans of increasing angular beam delivery were developed, three in Pinnacle v9.0 (4-field box, IMRT, and VMAT) and one in Accuray’s MultiPlan v3.5.3 (CyberKnife). The plans were delivered to the pelvic phantom containing both TLD and OSLD in the target volume. The pelvic phantom was also sent to two institutions to be irradiated as trials, one delivering IMRT, and the other a CyberKnife treatment. For the IMRT deliveries and the two institution trials, the phantom also included film in the sagittal and coronal planes. The doses measured from the TLD and OSLD were calculated for each irradiation, and the angular correction factors established from the spherical phantom irradiations were applied to the OSLD dose. The ratio of the TLD dose to the angular corrected OSLD dose was calculated for each irradiation. The corrected OSLD dose was found to be within 1% of the TLD measured dose for all irradiations, with the exception of the in-house CyberKnife deliveries. The films were normalized to both TLD measured dose and the corrected OSLD dose. Dose profiles were obtained and gamma analysis was performed using a 7%/4 mm criteria, to compare the ability of the OSLD, when corrected for the angular dependence, to provide equivalent results to TLD. The results of this study indicate that the OSLD can effectively be used as a replacement for TLD in the RPC’s anthropomorphic QA phantoms for coplanar treatment deliveries when a correction is applied for the dosimeter’s angular dependence.

Evaluation of Polymer Gel Dosimeters for Measurements of Dose and LET in Proton Beams

This project assessed the effectiveness of polymer gel dosimeters as tools for measuring the dose deposited by and LET of a proton beam. A total of three BANG® dosimeter formulations were evaluated: BANG®-3-Pro-2 BANGkits™ for dose measurement and two BANG®-3 variants, the LET-Baseline and LET-Meter dosimeters, for LET measurement. All dosimeters were read out using an OCT scanner. The basic characteristics of the BANGkits™ were assessed in a series of photon and electron irradiations. The dose-response relationship was found to be sigmoidal with a threshold for response of approximately 15 cGy. The active region of the dosimeter, the volume in which dosimeter response is not inhibited by oxygen, was found to make up roughly one fourth of the total dosimeter volume. Delivering a dose across multiple fractions was found to yield a greater response than delivering the same dose in a single irradiation. The dosimeter was found to accurately measure a dose distribution produced by overlapping photon fields, yielding gamma pass rates of 95.4% and 93.1% from two planar gamma analyses. Proton irradiations were performed for measurements of proton dose and LET. Initial irradiations performed through the side of a dosimeter led to OCT artifacts. Gamma pass rates of 85.7% and 89.9% were observed in two planar gamma analyses. In irradiations performed through the base of a dosimeter, gel response was found to increase with height in the dosimeter, even in areas of constant dose. After a correction was applied, gamma pass rates of 94.6% and 99.3% were observed in two planar gamma analyses. Absolute dose measurements were substantially higher (33%-100%) than the delivered doses for proton irradiations. Issues encountered while calibrating the LET-Meter gel restricted analysis of the LET measurement data to the SOBP of a proton beam. LET-Meter overresponse was found to increase linearly with track-average LET across the LET range that could be investigated (1.5 keV/micron – 3.5 keV/micron).

Computer simulation of glioma growth and morphology.

Despite major advances in the study of glioma, the quantitative links between intra-tumor molecular/cellular properties, clinically observable properties such as morphology, and critical tumor behaviors such as growth and invasiveness remain unclear, hampering more effective coupling of tumor physical characteristics with implications for prognosis and therapy. Although molecular biology, histopathology, and radiological imaging are employed in this endeavor, studies are severely challenged by the multitude of different physical scales involved in tumor growth, i.e., from molecular nanoscale to cell microscale and finally to tissue centimeter scale. Consequently, it is often difficult to determine the underlying dynamics across dimensions. New techniques are needed to tackle these issues. Here, we address this multi-scalar problem by employing a novel predictive three-dimensional mathematical and computational model based on first-principle equations (conservation laws of physics) that describe mathematically the diffusion of cell substrates and other processes determining tumor mass growth and invasion. The model uses conserved variables to represent known determinants of glioma behavior, e.g., cell density and oxygen concentration, as well as biological functional relationships and parameters linking phenomena at different scales whose specific forms and values are hypothesized and calculated based on in vitro and in vivo experiments and from histopathology of tissue specimens from human gliomas. This model enables correlation of glioma morphology to tumor growth by quantifying interdependence of tumor mass on the microenvironment (e.g., hypoxia, tissue disruption) and on the cellular phenotypes (e.g., mitosis and apoptosis rates, cell adhesion strength). Once functional relationships between variables and associated parameter values have been informed, e.g., from histopathology or intra-operative analysis, this model can be used for disease diagnosis/prognosis, hypothesis testing, and to guide surgery and therapy. In particular, this tool identifies and quantifies the effects of vascularization and other cell-scale glioma morphological characteristics as predictors of tumor-scale growth and invasion.

Computer simulation of glioma growth and morphology.

Despite major advances in the study of glioma, the quantitative links between intra-tumor molecular/cellular properties, clinically observable properties such as morphology, and critical tumor behaviors such as growth and invasiveness remain unclear, hampering more effective coupling of tumor physical characteristics with implications for prognosis and therapy. Although molecular biology, histopathology, and radiological imaging are employed in this endeavor, studies are severely challenged by the multitude of different physical scales involved in tumor growth, i.e., from molecular nanoscale to cell microscale and finally to tissue centimeter scale. Consequently, it is often difficult to determine the underlying dynamics across dimensions. New techniques are needed to tackle these issues. Here, we address this multi-scalar problem by employing a novel predictive three-dimensional mathematical and computational model based on first-principle equations (conservation laws of physics) that describe mathematically the diffusion of cell substrates and other processes determining tumor mass growth and invasion. The model uses conserved variables to represent known determinants of glioma behavior, e.g., cell density and oxygen concentration, as well as biological functional relationships and parameters linking phenomena at different scales whose specific forms and values are hypothesized and calculated based on in vitro and in vivo experiments and from histopathology of tissue specimens from human gliomas. This model enables correlation of glioma morphology to tumor growth by quantifying interdependence of tumor mass on the microenvironment (e.g., hypoxia, tissue disruption) and on the cellular phenotypes (e.g., mitosis and apoptosis rates, cell adhesion strength). Once functional relationships between variables and associated parameter values have been informed, e.g., from histopathology or intra-operative analysis, this model can be used for disease diagnosis/prognosis, hypothesis testing, and to guide surgery and therapy. In particular, this tool identifies and quantifies the effects of vascularization and other cell-scale glioma morphological characteristics as predictors of tumor-scale growth and invasion.

Trends in utilization of computer tomography (CT) scans for victims of motor vehicle collision (MVCs) in a level I trauma center from 1996 to 2006

Purpose. To evaluate trends in the utilization of head, abdominal, thoracic and other body regions CTs in the management of victims of MVC at a level I trauma center from 1996 to 2006.^ Method. From the trauma registry, I identified patients involved in MVC's in a level I trauma center and categorized them into three age groups of 13-18, 19-55 and ≥56. I used International Classification of Disease (ICD-9-CM) codes to find the type and number of CTs examinations performed for each patient. I plotted the mean number of CTs per patient against year of admission to find the crude estimate of change in utilization pattern for each type of CT. I used logistic regression to assess whether repetitive CTs (≥ 2) for head, abdomen, thorax and other body regions were associated with age group and year of admission for MVC patients. I adjusted the estimates for gender, ethnicity, insurance status, mechanism and severity of injury, intensive care unit admission status, patient disposition (dead or alive) and year of admission.^ Results. Utilization of head, abdominal, thoracic and other body regions CTs significantly increased over 11-year period. Utilization of head CT was greatest in the 13-18 age group, and increased from 0.58 CT/patient in 1996 to 1.37 CT/patient in 2006. Abdominal CTs were more common in the ≥56+ age group, and increased from 0.33 CT/patient in 1996 to 0.72 CT/patient in 2006. Utilization of thoracic CTs was higher in the 56+ age group, and increased from 0.01 CT/patient in 1996 to 0.42 CT/patient in 2006. Utilization of other CTs did not change materially during the study period for adolescents, adults or older adults. In the multivariable analysis, after adjustment for potential confounders, repetitive head CTs significantly increased in the 13-18 age group (95% CI: 1.29-1.87, p=<0.001) relative to the 19-55 age group. Repetitive thoracic CT use was lower in adolescents (95% CI: 0.22-0.70, p=<0.001) relative to the 19-55 age group.^ Conclusion. There has been a substantial increase in the utilization of head, abdominal, thoracic and other CTs in the management of MVC patients. Future studies need to identify if increased utilization of CTs have resulted in better health outcome for these patients. ^