Risk of Second Malignant Neoplasms Following Proton Arc Therapy and Volumetric Modulated Arc Therapy for Prostate Cancer

The risk of second malignant neoplasms (SMNs) following prostate radiotherapy is a concern due to the large population of survivors and decreasing age at diagnosis. It is known that parallel-opposed beam proton therapy carries a lower risk than photon IMRT. However, a comparison of SMN risk following proton and photon arc therapies has not previously been reported. The purpose of this study was to predict the ratio of excess relative risk (RRR) of SMN incidence following proton arc therapy to that after volumetric modulated arc therapy (VMAT). Additionally, we investigated the impact of margin size and the effect of risk-minimized proton beam weighting on predicted RRR. Physician-approved treatment plans were created for both modalities for three patients. Therapeutic dose was obtained with differential dose-volume histograms from the treatment planning system, and stray dose was estimated from the literature or calculated with Monte Carlo simulations. Then, various risk models were applied to the total dose. Additional treatment plans were also investigated with varying margin size and risk-minimized proton beam weighting. The mean RRR ranged from 0.74 to 0.99, depending on risk model. The additional treatment plans revealed that the RRR remained approximately constant with varying margin size, and that the predicted RRR was reduced by 12% using a risk-minimized proton arc therapy planning technique. In conclusion, proton arc therapy was found to provide an advantage over VMAT in regard to predicted risk of SMN following prostate radiotherapy. This advantage was independent of margin size and was amplified with risk-optimized proton beam weighting.

Histamine reduces flash sensitivity of on ganglion cells in the primate retina.

PURPOSE. In Old World primates, the retina receives input from histaminergic neurons in the posterior hypothalamus. They are a subset of the neurons that project throughout the central nervous system and fire maximally during the day. The contribution of these neurons to vision, was examined by applying histamine to a dark-adapted, superfused baboon eye cup preparation while making extracellular recordings from peripheral retinal ganglion cells. METHODS. The stimuli were 5-ms, 560-nm, weak, full-field flashes in the low scotopic range. Ganglion cells with sustained and transient ON responses and two cell types with OFF responses were distinguished; their responses were recorded with a 16-channel microelectrode array. RESULTS. Low micromolar doses of histamine decreased the rate of maintained firing and the light sensitivity of ON ganglion cells. Both sustained and transient ON cells responded similarly to histamine. There were no statistically significant effects of histamine in a more limited study of OFF ganglion cells. The response latencies of ON cells were approximately 5 ms slower, on average, when histamine was present. Histamine also reduced the signal-to-noise ratio of ON cells, particularly in those cells with a histamine-induced increase in maintained activity. CONCLUSIONS. A major action of histamine released from retinopetal axons under dark-adapted conditions, when rod signals dominate the response, is to reduce the sensitivity of ON ganglion cells to light flashes. These findings may relate to reports that humans are less sensitive to light stimuli in the scotopic range during the day, when histamine release in the retina is expected to be at its maximum.

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.

Comparison of a traditional and a multilevel Cox Proportional Hazards model

Hierarchically clustered populations are often encountered in public health research, but the traditional methods used in analyzing this type of data are not always adequate. In the case of survival time data, more appropriate methods have only begun to surface in the last couple of decades. Such methods include multilevel statistical techniques which, although more complicated to implement than traditional methods, are more appropriate. ^ One population that is known to exhibit a hierarchical structure is that of patients who utilize the health care system of the Department of Veterans Affairs where patients are grouped not only by hospital, but also by geographic network (VISN). This project analyzes survival time data sets housed at the Houston Veterans Affairs Medical Center Research Department using two different Cox Proportional Hazards regression models, a traditional model and a multilevel model. VISNs that exhibit significantly higher or lower survival rates than the rest are identified separately for each model. ^ In this particular case, although there are differences in the results of the two models, it is not enough to warrant using the more complex multilevel technique. This is shown by the small estimates of variance associated with levels two and three in the multilevel Cox analysis. Much of the differences that are exhibited in identification of VISNs with high or low survival rates is attributable to computer hardware difficulties rather than to any significant improvements in the model. ^

Effects of informative censoring on the hazard function: Application to the proportional hazards regression model

The standard analyses of survival data involve the assumption that survival and censoring are independent. When censoring and survival are related, the phenomenon is known as informative censoring. This paper examines the effects of an informative censoring assumption on the hazard function and the estimated hazard ratio provided by the Cox model.^ The limiting factor in all analyses of informative censoring is the problem of non-identifiability. Non-identifiability implies that it is impossible to distinguish a situation in which censoring and death are independent from one in which there is dependence. However, it is possible that informative censoring occurs. Examination of the literature indicates how others have approached the problem and covers the relevant theoretical background.^ Three models are examined in detail. The first model uses conditionally independent marginal hazards to obtain the unconditional survival function and hazards. The second model is based on the Gumbel Type A method for combining independent marginal distributions into bivariate distributions using a dependency parameter. Finally, a formulation based on a compartmental model is presented and its results described. For the latter two approaches, the resulting hazard is used in the Cox model in a simulation study.^ The unconditional survival distribution formed from the first model involves dependency, but the crude hazard resulting from this unconditional distribution is identical to the marginal hazard, and inferences based on the hazard are valid. The hazard ratios formed from two distributions following the Gumbel Type A model are biased by a factor dependent on the amount of censoring in the two populations and the strength of the dependency of death and censoring in the two populations. The Cox model estimates this biased hazard ratio. In general, the hazard resulting from the compartmental model is not constant, even if the individual marginal hazards are constant, unless censoring is non-informative. The hazard ratio tends to a specific limit.^ Methods of evaluating situations in which informative censoring is present are described, and the relative utility of the three models examined is discussed. ^

Assessing time-by-covariate interactions in Cox proportional hazards regression models using cubic spline functions

This dissertation develops and explores the methodology for the use of cubic spline functions in assessing time-by-covariate interactions in Cox proportional hazards regression models. These interactions indicate violations of the proportional hazards assumption of the Cox model. Use of cubic spline functions allows for the investigation of the shape of a possible covariate time-dependence without having to specify a particular functional form. Cubic spline functions yield both a graphical method and a formal test for the proportional hazards assumption as well as a test of the nonlinearity of the time-by-covariate interaction. Five existing methods for assessing violations of the proportional hazards assumption are reviewed and applied along with cubic splines to three well known two-sample datasets. An additional dataset with three covariates is used to explore the use of cubic spline functions in a more general setting. ^

A SIMULATION STUDY OF THE SIZE AND POWER OF SOME TWO-SAMPLE TESTS FOR UNEQUALLY CENSORED SURVIVAL DATA (RANK TESTS, LOGRANK, PROPORTIONAL HAZARDS, WILCOXON TESTS, EXPONENTIAL DISTRIBUTION)

Sizes and power of selected two-sample tests of the equality of survival distributions are compared by simulation for small samples from unequally, randomly-censored exponential distributions. The tests investigated include parametric tests (F, Score, Likelihood, Asymptotic), logrank tests (Mantel, Peto-Peto), and Wilcoxon-Type tests (Gehan, Prentice). Equal sized samples, n = 18, 16, 32 with 1000 (size) and 500 (power) simulation trials, are compared for 16 combinations of the censoring proportions 0%, 20%, 40%, and 60%. For n = 8 and 16, the Asymptotic, Peto-Peto, and Wilcoxon tests perform at nominal 5% size expectations, but the F, Score and Mantel tests exceeded 5% size confidence limits for 1/3 of the censoring combinations. For n = 32, all tests showed proper size, with the Peto-Peto test most conservative in the presence of unequal censoring. Powers of all tests are compared for exponential hazard ratios of 1.4 and 2.0. There is little difference in power characteristics of the tests within the classes of tests considered. The Mantel test showed 90% to 95% power efficiency relative to parametric tests. Wilcoxon-type tests have the lowest relative power but are robust to differential censoring patterns. A modified Peto-Peto test shows power comparable to the Mantel test. For n = 32, a specific Weibull-exponential comparison of crossing survival curves suggests that the relative powers of logrank and Wilcoxon-type tests are dependent on the scale parameter of the Weibull distribution. Wilcoxon-type tests appear more powerful than logrank tests in the case of late-crossing and less powerful for early-crossing survival curves. Guidelines for the appropriate selection of two-sample tests are given. ^

Volumetric Modulated Arc Therapy Evaluation with the Radiological Physics Center Head and Neck Phantom

Validation of treatment plan quality and dose calculation accuracy is essential for new radiotherapy techniques, including volumetric modulated arc therapy (VMAT). VMAT delivers intensity modulated radiotherapy treatments while simultaneously rotating the gantry, adding an additional level of complexity to both the dose calculation and delivery of VMAT treatments compared to static gantry IMRT. The purpose of this project was to compare two VMAT systems, Elekta VMAT and Varian RapidArc, to the current standard of care, IMRT, in terms of both treatment plan quality and dosimetric delivery accuracy using the Radiological Physics Center (RPC) head and neck (H&N) phantom. Clinically relevant treatment plans were created for the phantom using typical prescription and dose constraints for Elekta VMAT (planned with Pinnacle3 Smart Arc) and RapidArc and IMRT (both planned with Eclipse). The treatment plans were evaluated to determine if they were clinically comparable using several dosimetric criteria, including ability to meet dose objectives, hot spots, conformity index, and homogeneity index. The planned treatments were delivered to the phantom and absolute doses and relative dose distributions were measured with thermoluminescent dosimeters (TLDs) and radiochromic film, respectively. The measured and calculated doses of each treatment were compared to determine if they were clinically acceptable based upon RPC criteria of ±7% dose difference and 4 mm distance-to-agreement. Gamma analysis was used to assess dosimetric accuracy, as well. All treatment plans were able to meet the dosimetric objectives set by the RPC and had similar hot spots in the normal tissue. The Elekta VMAT plan was more homogenous but less conformal than the RapidArc and IMRT plans. When comparing the measured and calculated doses, all plans met the RPC ±7%/4 mm criteria. The percent of points passing the gamma analysis for each treatment delivery was acceptable. Treatment plan quality of the Elekta VMAT, RapidArc and IMRT treatments were comparable for consistent dose prescriptions and constraints. Additionally, the dosimetric accuracy of the Elekta VMAT and RapidArc treatments was verified to be within acceptable tolerances.