Second cancer incidence risk estimates using BEIR VII models for standard and complex external beam radiotherapy for early breast cancer.

PURPOSE: To compare organ specific cancer incidence risks for standard and complex external beam radiotherapy (including cone beam CT verification) following breast conservation surgery for early breast cancer. METHOD: Doses from breast radiotherapy and kilovoltage cone beam CT (CBCT) exposures were obtained from thermoluminescent dosimeter measurements in an anthropomorphic phantom in which the positions of radiosensitive organs were delineated. Five treatment deliveries were investigated: (i) conventional tangential field whole breast radiotherapy (WBRT), (ii) noncoplanar conformal delivery applicable to accelerated partial beast irradiation (APBI), (iii) two-volume simultaneous integrated boost (SIB) treatment, (iv) forward planned three-volume SIB, and (v) inverse-planned three volume SIB. Conformal and intensity modulated radiotherapy methods were used to plan the complex treatments. Techniques spanned the range from simple methods appropriate for patient cohorts with a low cancer recurrence risk to complex plans relevant to cohorts with high recurrence risk. Delineated organs at risk included brain, salivary glands, thyroid, contralateral breast, left and right lung, esophagus, stomach, liver, colon, and bladder. Biological Effects of Ionizing Radiation (BEIR) VII cancer incidence models were applied to the measured mean organ doses to determine lifetime attributable risk (LAR) for ages at exposure from 35 to 80 yr according to radiotherapy techniques, and included dose from the CBCT imaging. RESULTS: All LAR decreased with age at exposure and were lowest for brain, thyroid, liver, and bladder (<0.1%). There was little dependence of LAR on radiotherapy technique for these organs and for colon and stomach. LAR values for the lungs for the three SIB techniques were two to three times those from WBRT and APBI. Uncertainties in the LAR models outweigh any differences in lung LAR between the SIB methods. Constraints in the planning of the SIB methods ensured that contralateral breast doses and LAR were comparable to WBRT, despite their added complexity. The smaller irradiated volume of the ABPI plan contributed to a halving of LAR for contralateral breast compared with the other plan types. Daily image guided radiotherapy (IGRT) for a left breast protocol using kilovoltage CBCT contributed <10% to LAR for the majority of organs, and did not exceed 22% of total organ dose. CONCLUSIONS: Phantom measurements and calculations of LAR from the BEIR VII models predict that complex breast radiotherapy techniques do not increase the theoretical risk of second cancer incidence for organs distant from the treated breast, or the contralateral breast where appropriate plan constraints are applied. Complex SIB treatments are predicted to increase the risk of second cancer incidence in the lungs compared to standard whole breast radiotherapy; this is outweighed by the threefold reduction in 5 yr local recurrence risk for patients of high risk of recurrence, and young age, from the use of radiotherapy. APBI may have a favorable impact on risk of second cancer in the contralateral breast and lung for older patients at low risk of recurrence. Intensive use of IGRT increased the estimated values of LAR but these are dominated by the effect of the dose from the radiotherapy, and any increase in LAR from IGRT is much lower than the models' uncertainties.

Kilovoltage energy imaging with a radiotherapy linac with a continuously variable energy range.

PURPOSE: In this paper, the effect on image quality of significantly reducing the primary electron energy of a radiotherapy accelerator is investigated using a novel waveguide test piece. The waveguide contains a novel variable coupling device (rotovane), allowing for a wide continuously variable energy range of between 1.4 and 9 MeV suitable for both imaging and therapy. METHOD: Imaging at linac accelerating potentials close to 1 MV was investigated experimentally and via Monte Carlo simulations. An imaging beam line was designed, and planar and cone beam computed tomography images were obtained to enable qualitative and quantitative comparisons with kilovoltage and megavoltage imaging systems. The imaging beam had an electron energy of 1.4 MeV, which was incident on a water cooled electron window consisting of stainless steel, a 5 mm carbon electron absorber and 2.5 mm aluminium filtration. Images were acquired with an amorphous silicon detector sensitive to diagnostic x-ray energies. RESULTS: The x-ray beam had an average energy of 220 keV and half value layer of 5.9 mm of copper. Cone beam CT images with the same contrast to noise ratio as a gantry mounted kilovoltage imaging system were obtained with doses as low as 2 cGy. This dose is equivalent to a single 6 MV portal image. While 12 times higher than a 100 kVp CBCT system (Elekta XVI), this dose is 140 times lower than a 6 MV cone beam imaging system and 6 times lower than previously published LowZ imaging beams operating at higher (4-5 MeV) energies. CONCLUSIONS: The novel coupling device provides for a wide range of electron energies that are suitable for kilovoltage quality imaging and therapy. The imaging system provides high contrast images from the therapy portal at low dose, approaching that of gantry mounted kilovoltage x-ray systems. Additionally, the system provides low dose imaging directly from the therapy portal, potentially allowing for target tracking during radiotherapy treatment. There is the scope with such a tuneable system for further energy reduction and subsequent improvement in image quality.

An end-to-end assessment on the accuracy of adaptive radiotherapy in an MR-linac.

PURPOSE: To develop and demonstrate an end-to-end assessment procedure for adaptive radiotherapy (ART) within an MR-guided system. METHODS AND MATERIALS: A 3D printed pelvic phantom was designed and constructed for use in this study. The phantom was put through the complete radiotherapy treatment chain, with planned internal changes made to model prostate translations and shape changes, allowing an investigation into three ART techniques commonly used. Absolute dosimetry measurements were made within the phantom using both gafchromic film and alanine. Comparisons between treatment planning system (TPS) calculations and measured dose values were made using the gamma evaluation with criteria of 3 mm/3% and 2 mm/2%. RESULTS: Gamma analysis evaluations for each type of treatment plan adaptation investigated showed a very high agreement with pass rates for each experiment ranging from 98.10% to 99.70% and 92.60% to 97.55%, for criteria of 3%/3 mm and 2%/2 mm respectively. These pass rates were consistent for both shape and position changes. Alanine measurements further supported the results, showing an average difference of 1.98% from the TPS. CONCLUSION: The end-to-end assessment procedure provided demanding challenges for treatment plan adaptations to demonstrate the capabilities and achieved high consistency in all findings.

The use of PET images for radiotherapy treatment planning: an error analysis using radiobiological endpoints.

PURPOSE: There is significant current interest in the use of biological image guidance in radiotherapy planning. In lung-cancer treatment, tumor motion due to respiration is known to be a limitation. This is particularly true in PET, where image data are collected over a number of minutes. An in-house-developed 4D PET acquisition mode is described and an analysis of the effects of acquisition parameters on the reconstructed image quality is presented. The potential impact of the resulting biological image quality on radiotherapy planning is then quantified in terms of tumor control probability (TCP). METHODS: Data were acquired using a human torso phantom comprised of a hot '8F-filled spheroidal "tumor" (40 mm in diameter) suspended in an air-filled "lung" cylinder and surrounded by a warm 18F-filled background. Two different sphere-to-background (S/B) ratios were used. The tumor was connected to a 3-axis computer-controlled motion stage and could be moved during PET data acquisition. Images were acquired with a range of count statistics, motion blurring, and CT attenuation correction (CTAC) misalignment. Four simple models were proposed for the assignment of clonogenic cell density according to the voxel value. The impact of image artifacts was then assessed by calculating the TCP, which is the probability that no clonogenic tumor cell remains after a given dose of radiation. TCP was calculated for a uniform dose distribution in the tumor. RESULTS: Reduced count statistics and misaligned CTAC images had the most detrimental impact on the image fidelity. It was found that in both cases the images became less intense, demonstrated by smaller number of voxels at the maximum values. The maximum TCP difference between images with the least and most noise was 3.4% (S/B=3), and with weakest and strongest CT misalignment artifacts, it was 3.2% (S/B = 10). Motion blurring only contributed weakly to the TCP imprecision at 1.7% (S/B=10) between best- and worst-case images. However, the model-calculated TCP showed increasing differences from the ground truth as the complexity of the model increased [maximum difference of approximately 8% (model 3)], which could be attributed to the partial volume effect. CONCLUSIONS: Based on the results of this study, it is believed that simple techniques of biologically guided radiotherapy planning for lung cancer should be feasible at intermediate contrast levels (tumor-to-background ratio of approximately 10) with the clinically achievable image quality.

Dose prescription complexity versus tumor control probability in biologically conformal radiotherapy.

The technical feasibility and potential benefits of voxel-based nonuniform dose prescriptions for biologically heterogeneous tumors have been widely demonstrated. In some cases, an "ideal" dose prescription has been generated by individualizing the dose to every voxel within the target, but often this voxel-based prescription has been discretized into a small number of compartments. The number of dose levels utilized and the methods used for prescribing doses and assigning tumor voxels to different dose compartments have varied significantly. The authors present an investigation into the relationship between the complexity of the dose prescription and the tumor control probability (TCP) for a number of these methods. The linear quadratic model of cell killing was used in conjunction with a number of modeled tumors heterogeneous in clonogen density, oxygenation, or proliferation. Models based on simple mathematical functions, published biological data, and biological image data were investigated. Target voxels were assigned to dose compartments using (i) simple rules based on the initial biological distribution, (ii) iterative methods designed to maximize the achievable TCP, or (iii) methods based on an ideal dose prescription. The relative performance of the simple rules was found to depend on the form of heterogeneity of the tumor, while the iterative and ideal dose methods performed comparably for all models investigated. In all cases the maximum achievable TCP was approached within the first few (typically two to five) compartments. Results suggest that irrespective of the pattern of heterogeneity, the optimal dose prescription can be well approximated using only a few dose levels but only if both the compartment boundaries and prescribed dose levels are well chosen.

Electron beam quality control using an amorphous silicon EPID.

An amorphous silicon EPID has been investigated to determine whether it is capable of quality control constancy measurements for linear accelerator electron beams. The EPID grayscale response was found to be extremely linear with dose over a wide dose range and, more specifically, for exposures of 95-100 MU. Small discrepancies of up to 0.8% in linearity were found at 6 MeV (8-15 MeV showed better agreement). The shape of the beam profile was found to be significantly altered by scatter in air over the approximately 60 cm gap between the end of the applicator and the EPID. Nevertheless, relative changes in EPID-measured profile flatness and symmetry were linearly related to changes in these parameters at 95 cm focus to surface distance (FSD) measured using a 2D diode array. Similar results were obtained at 90 degrees and 270 degrees gantry angles. Six months of daily images were acquired and analyzed to determine whether the device is suitable as a constancy checker. EPID output measurements agreed well with daily ion chamber measurements, with a 0.8% standard deviation in the difference between the two measurement sets. When compared to weekly parallel plate chamber measurements, this figure dropped to 0.5%. A Monte Carlo (MC) model of the EPID was created and demonstrated excellent agreement between MC-calculated profiles in water and the EPID at 95 and 157 cm FSD. Good agreement was also found with measured EPID profiles, demonstrating that the EPID provides an accurate measurement of electron profiles. The EPID was thus shown to be an effective method for performing electron beam daily constancy checks.

An efficient Monte Carlo-based algorithm for scatter correction in keV cone-beam CT.

A new method is proposed for scatter-correction of cone-beam CT images. A coarse reconstruction is used in initial iteration steps. Modelling of the x-ray tube spectra and detector response are included in the algorithm. Photon diffusion inside the imaging subject is calculated using the Monte Carlo method. Photon scoring at the detector is calculated using forced detection to a fixed set of node points. The scatter profiles are then obtained by linear interpolation. The algorithm is referred to as the coarse reconstruction and fixed detection (CRFD) technique. Scatter predictions are quantitatively validated against a widely used general-purpose Monte Carlo code: BEAMnrc/EGSnrc (NRCC, Canada). Agreement is excellent. The CRFD algorithm was applied to projection data acquired with a Synergy XVI CBCT unit (Elekta Limited, Crawley, UK), using RANDO and Catphan phantoms (The Phantom Laboratory, Salem NY, USA). The algorithm was shown to be effective in removing scatter-induced artefacts from CBCT images, and took as little as 2 min on a desktop PC. Image uniformity was greatly improved as was CT-number accuracy in reconstructions. This latter improvement was less marked where the expected CT-number of a material was very different to the background material in which it was embedded.

Target-tracking deliveries on an Elekta linac: a feasibility study.

A target-tracking, intensity-modulated delivery on an Elekta MLCi system was assessed by film measurement with a simulated target-motion trajectory. A toroidally shaped idealized target surrounding an organ at risk necessitating multiple field segments to irradiate the target and spare the organ at risk was defined in a solid-water phantom. The phantom was programmed to move following a reproducible 2D elliptical trajectory in the beam's-eye view with a period of 10 s. Static and target-tracking treatments were planned for delivery on a standard Elekta Precise series linac with integrated MLCi system. Dose was delivered in three ways: (i) a static treatment to a static phantom, (ii) a static treatment to a moving phantom and (iii) a target-tracking treatment to a moving phantom. The dose delivered was assessed by film measurement on the central plane through the target and organ at risk. The target dose blurring was quantified by the standard deviation of the dose to the target which was evaluated as 2.8% for the static treatment to the static phantom, 7.2% for the static treatment to the moving phantom and 2.6% for the tracking treatment to the moving phantom. The mean organ-at-risk dose was 38.2%, 54.0% and 38.2% of the prescription dose for each delivery case. We have therefore shown that the linac is capable of delivering target-tracking fields with MLCs for the target trajectories tested.

The stability of imaging biomarkers in radiomics: a framework for evaluation.

This paper studies the sensitivity of a range of image texture parameters used in radiomics to: (i) the number of intensity levels, (ii) the method of quantisation to select the intensity levels and (iii) the use of an intensity threshold. 43 commonly used texture features were studied for the gross target volume outlined on the CT component of PET/CT scans of 50 patients with non-small cell lung carcinoma (NSCLC). All cases were quantised for all values between 4 and 128 intensity levels using four commonly used quantisation methods. All results were analysed with and without a threshold range of -200 HU to 300 HU. Cases were ranked for each texture feature and for all quantisation methods with the Spearman's rank correlation coefficient determined to evaluate stability. Results showed large fluctuations in ranking, particularly for low numbers of levels, differences between quantisation methods and with the use of a threshold, with values Spearman's Rank Correlation for many parameters below 0.2. Our results demonstrated the sensitivity of radiomics features to the parameters used during analysis and highlight the risk of low reproducibility comparing studies with slightly different parameters. In terms of the lung cancer CT datasets, this study supports the use of 128 intensity levels, the same uniform quantiser applied to all scans and thresholding of the data. It also supports several of the features recommended in the literature for such studies such as skewness and kurtosis. A recommended framework is presented for curation of the data analysis process to ensure stability of results.

A theoretical framework for prescribing radiotherapy dose distributions using patient-specific biological information.

We present a formalism for using functional imaging both to derive patient-specific radiobiological properties and consequently to prescribe optimal nonuniform radiotherapy dose distributions. The ability to quantitatively assess the response to an initial course of radiotherapy would allow the derivation of radiobiological parameters for individual patients. Both an iterative optimization and an analytical approach to this problem were investigated and illustrated by application to the linear-quadratic model of cell killing using simulated parametric data for a modeled tumor. Potential gains in local control were assessed by comparing uniform dose distributions with optimized dose distributions of equal integral dose. The effect on local prescribed dose of variations in effective radiosensitivity, tumor burden, and proliferation rate was investigated, with results suggesting that dose variations would be significant but clinically achievable. The sensitivity of derived parameters to image noise and the effect of varying the initial fractionation and imaging schedule were assessed. The analytical approach proved remarkably robust, with 10% image noise resulting in dose errors of approximately 1% for a clinically relevant set of parameters. Potential benefits were demonstrated by using this formalism to prescribe nonuniform dose distributions for model tumors using a range of literature-derived parameters. The redistribution of dose improved tumor control probability by factors between 1.03 and 4.27 for a range of model tumors.

An investigation into the use of CMOS active pixel technology in image-guided radiotherapy.

The increased intelligence, read-out speed, radiation hardness and potential large size of CMOS active pixel sensors (APS) gives them a potential advantage over systems currently used for verification of complex treatments such as IMRT and the tracking of moving tumours. The aim of this work is to investigate the feasibility of using an APS-based system to image the megavoltage treatment beam produced by a linear accelerator (Linac), and to demonstrate the logic which may ultimately be incorporated into future sensor and FPGA design to evaluate treatment and track motion. A CMOS APS was developed by the MI(3) consortium and incorporated into a megavoltage imaging system using the standard lens and mirror configuration employed in camera-based EPIDs. The ability to resolve anatomical structure was evaluated using an Alderson RANDO head phantom, resolution evaluated using a quality control (QC3) phantom and contrast using an in-house developed phantom. A complex intensity-modulated radiotherapy (IMRT) treatment was imaged and two algorithms were used to determine the field-area and delivered dose, and the position of multi-leaf collimator (MLC) leaves off-line. Results were compared with prediction from the prescription and found to agree within a single image frame time for dose delivery and 0.02-0.03 cm for the position of collimator leaves. Such a system therefore shows potential as the basis for an on-line verification system capable of treatment verification and monitoring patient motion.

A low Z linac and flat panel imager: comparison with the conventional imaging approach.

Experimental and Monte Carlo simulations were conducted for an Elekta Ltd Precise Treatment System linac fitted with a low Z insert of sufficient thickness to remove all primary electrons. A variety of amorphous silicon based panels employing different scintillators were modelled to determine their response to a variety of x-ray spectra and produce an optimized portal imaging system. This study has shown that in a low Z configuration the vast majority of x-rays are produced in the nickel electron window, and with a combination of a carbon insert and caesium iodide based XVI-panel, significant improvement in the object contrast was achieved. For thin, head and neck-type geometries, contrast is 4.62 times greater for 1.6 cm bone in 5.8 cm water than the standard 6 MV/iViewGT system. For thicker, pelvis-type geometries contrast increases by a factor of 1.3 for 1.6 cm of bone in 25.8 cm water. To obtain images with the same signal-to-noise ratio as the 6 MV/iViewGT system, dose reductions of a factor of 15 and 4.2 are possible for 5.8 cm and 25.8 cm phantoms respectively. This design has the advantage of being easily implemented on a standard linac and provides a portal image directly from the therapy beam aperture.

Planning lung radiotherapy using 4D CT data and a motion model.

This work is a feasibility study to use a four-dimensional computed tomography (4D CT) dataset generated by a continuous motion model for treatment planning in lung radiotherapy. The model-based 4D CT data were derived from multiple breathing cycles. Four patients were included in this retrospective study. Treatment plans were optimized at end-exhale for each patient and the effect of respiratory motion on the dose delivery investigated. The accuracy of the delivered dose as determined by the number of intermediate respiratory phases used for the calculation was considered. The time-averaged geometry of the anatomy representing the mid-ventilation phase of the breathing cycle was generated using the motion model and a treatment plan was optimized for this phase for one patient. With respiratory motion included, the mid-ventilation plan achieved better target coverage than the plan optimized at end-exhale when standard margins were used to expand the clinical target volume (CTV) to planning target volume (PTV). Using a margin to account for set-up uncertainty only, resulted in poorer target coverage and healthy tissue sparing. For this patient cohort, the results suggest that conventional three-dimensional treatment planning was sufficient to maintain target coverage despite respiratory motion. The motion model has proved a useful tool in 4D treatment planning.

Amorphous silicon EPID calibration for dosimetric applications: comparison of a method based on Monte Carlo prediction of response with existing techniques.

For EPID dosimetry, the calibration should ensure that all pixels have a similar response to a given irradiation. A calibration method (MC), using an analytical fit of a Monte Carlo simulated flood field EPID image to correct for the flood field image pixel intensity shape, was proposed. It was compared with the standard flood field calibration (FF), with the use of a water slab placed in the beam to flatten the flood field (WS) and with a multiple field calibration where the EPID was irradiated with a fixed 10x10 field for 16 different positions (MF). The EPID was used in its normal configuration (clinical setup) and with an additional 3 mm copper slab (modified setup). Beam asymmetry measured with a diode array was taken into account in MC and WS methods. For both setups, the MC method provided pixel sensitivity values within 3% of those obtained with the MF and WS methods (mean difference<1%, standard deviation<2%). The difference of pixel sensitivity between MC and FF methods was up to 12.2% (clinical setup) and 11.8% (modified setup). MC calibration provided images of open fields (5x5 to 20x20 cm2) and IMRT fields to within 3% of that obtained with WS and MF calibrations while differences with images calibrated with the FF method for fields larger than 10x10 cm2 were up to 8%. MC, WS and MF methods all provided a major improvement on the FF method. Advantages and drawbacks of each method were reviewed.

Thresholds of biodiversity and ecosystem function in a forest ecosystem undergoing dieback.

Ecological thresholds, which represent points of rapid change in ecological properties, are of major scientific and societal concern. However, very little research has focused on empirically testing the occurrence of thresholds in temperate terrestrial ecosystems. To address this knowledge gap, we tested whether a number of biodiversity, ecosystem functions and ecosystem condition metrics exhibited thresholds in response to a gradient of forest dieback, measured as changes in basal area of living trees relative to areas that lacked recent dieback. The gradient of dieback was sampled using 12 replicate study areas in a temperate forest ecosystem. Our results provide novel evidence of several thresholds in biodiversity (namely species richness of ectomycorrhizal fungi, epiphytic lichen and ground flora); for ecological condition (e.g. sward height, palatable seedling abundance) and a single threshold for ecosystem function (i.e. soil respiration rate). Mechanisms for these thresholds are explored. As climate-induced forest dieback is increasing worldwide, both in scale and speed, these results imply that threshold responses may become increasingly widespread.

The UK HeartSpare Study (Stage II): Multicentre Evaluation of a Voluntary Breath-hold Technique in Patients Receiving Breast Radiotherapy.

AIMS: To evaluate the feasibility and heart-sparing ability of the voluntary breath-hold (VBH) technique in a multicentre setting. MATERIALS AND METHODS: Patients were recruited from 10 UK centres. Following surgery for early left breast cancer, patients with any heart inside the 50% isodose from a standard free-breathing tangential field treatment plan underwent a second planning computed tomography (CT) scan using the VBH technique. A separate treatment plan was prepared on the VBH CT scan and used for treatment. The mean heart, left anterior descending coronary artery (LAD) and lung doses were calculated. Daily electronic portal imaging (EPI) was carried out and scanning/treatment times were recorded. The primary end point was the percentage of patients achieving a reduction in mean heart dose with VBH. Population systematic (Σ) and random errors (σ) were estimated. Within-patient comparisons between techniques used Wilcoxon signed-rank tests. RESULTS: In total, 101 patients were recruited during 2014. Primary end point data were available for 93 patients, 88 (95%) of whom achieved a reduction in mean heart dose with VBH. Mean cardiac doses (Gy) for free-breathing and VBH techniques, respectively, were: heart 1.8 and 1.1, LAD 12.1 and 5.4, maximum LAD 35.4 and 24.1 (all P<0.001). Population EPI-based displacement data showed Σ =+1.3-1.9 mm and σ=1.4-1.8 mm. Median CT and treatment session times were 21 and 22 min, respectively. CONCLUSIONS: The VBH technique is confirmed as effective in sparing heart tissue and is feasible in a multicentre setting.

Assessing the effect of electron density in photon dose calculations.

Photon dose calculation algorithms (such as the pencil beam and collapsed cone, CC) model the attenuation of a primary photon beam in media other than water, by using pathlength scaling based on the relative mass density of the media to water. In this study, we assess if differences in the electron density between the water and media, with different atomic composition, can influence the accuracy of conventional photon dose calculations algorithms. A comparison is performed between an electron-density scaling method and the standard mass-density scaling method for (i) tissues present in the human body (such as bone, muscle, etc.), and for (ii) water-equivalent plastics, used in radiotherapy dosimetry and quality assurance. We demonstrate that the important material property that should be taken into account by photon dose algorithms is the electron density, and not the mass density. The mass-density scaling method is shown to overestimate, relative to electrondensity predictions, the primary photon fluence for tissues in the human body and water-equivalent plastics, where 6%-7% and 10% differences were observed respectively for bone and air. However, in the case of patients, differences are expected to be smaller due to the large complexity of a treatment plan and of the patient anatomy and atomic composition and of the smaller thickness of bone/air that incident photon beams of a treatment plan may have to traverse. Differences have also been observed for conventional dose algorithms, such as CC, where an overestimate of the lung dose occurs, when irradiating lung tumors. The incorrect lung dose can be attributed to the.incorrect modeling of the photon beam attenuation through the rib cage (thickness of 2-3 cm in bone upstream of the lung tumor) and through the lung and the oversimplified modeling of electron transport in convolution algorithms. In the present study, the overestimation of the primary photon fluence, using the mass-density scaling method, was shown to be a consequence of the differences in the hydrogen content between the various media studied and water. On the other hand, the electron-density scaling method was shown to predict primary photon fluence in media other than water to within 1%-2% for all the materials studied and for energies up to 5 MeV. For energies above 5 MeV, the accuracy of the electron-density scaling method was shown to depend on the photon energy, where for materials with a high content of calcium (such as bone, cortical bone) or for primary photon energies above 10 MeV, the pair-production process could no longer be neglected. The electron-density scaling method was extended to account for pair-production attenuation of the primary photons. Therefore the scaling of the dose distributions in media other than water became dependent on the photon energy. The extended electron-scaling method was shown to estimate the photon range to within 1% for all materials studied and for energies from 100 keV to 20 MeV, allowing it to be used to scale dose distributions to media other than water and generated by clinical radiotherapy photon beams with accelerator energies from 4 to 20 MV.

A quantitative study of IMRT delivery effects in commercial planning systems for the case of oesophagus and prostate tumours.

This study focuses on understanding the impact of intensity-modulated radiotherapy (IMRT) delivery effects when applied to plans generated by commercial treatment-planning systems such as Pinnacle (ADAC Laboratories Inc.) and CadPlan/Helios (Varian Medical Systems). These commercial planning systems have had several version upgrades (with improvements in the optimization algorithm), but the IMRT delivery effects have not been incorporated into the optimization process. IMRT delivery effects include head-scatter fluence from IMRT fields, transmission through leaves and the effect of the rounded shape of the leaf ends. They are usually accounted for after optimization when leaf sequencing the "optimal" fluence profiles, to derive the delivered fluence profile. The study was divided into two main parts: (a) analysing the dose distribution within the planning-target volume (PTV), produced by each of the commercial treatment-planning systems, after the delivered fluence had been renormalized to deliver the correct dose to the PTV; and (b) studying the impact of the IMRT delivery technique on the surrounding critical organs such as the spinal cord, lungs, rectum, bladder etc. The study was performed for tumours of (i) the oesophagus and (ii) the prostate and pelvic nodes. An oesophagus case was planned with the Pinnacle planning system for IMRT delivery, via multiple-static fields (MSF) and compensators, using the Elekta SL25 with a multileaf collimator (MLC) component. A prostate and pelvic nodes IMRT plan was performed with the Cadplan/Helios system for a dynamic delivery (DMLC) using the Varian 120-leaf Millennium MLC. In these commercial planning systems, since IMRT delivery effects are not included into the optimization process, fluence renormalization is required such that the median delivered PTV dose equals the initial prescribed PTV dose. In preparing the optimum fluence profile for delivery, the PTV dose has been "smeared" by the IMRT delivery techniques. In the case of the oesophagus, the critical organ, spinal cord, received a greater dose than initially planned, due to the delivery effects. The increase in the spinal cord dose is of the order of 2-3 Gy. In the case of the prostate and pelvic nodes, the IMRT delivery effects led to an increase of approximately 2 Gy in the dose delivered to the secondary PTV, the pelvic nodes. In addition to this, the small bowel, rectum and bladder received an increased dose of the order of 2-3 Gy to 50% of their total volume. IMRT delivery techniques strongly influence the delivered dose distributions for the oesophagus and prostate/pelvic nodes tumour sites and these effects are not yet accounted for in the Pinnacle and the CadPlan/Helios planning systems. Currently, they must be taken into account during the optimization stage by altering the dose limits accepted during optimization so that the final (sequenced) dose is within the constraints.

An experimental demonstration of a new type of proton computed tomography using a novel silicon tracking detector.

PURPOSE: Radiography and tomography using proton beams promise benefit to image guidance and treatment planning for proton therapy. A novel proton tracking detector is described and experimental demonstrations at a therapy facility are reported. A new type of proton CT reconstructing relative "scattering power" rather than "stopping power" is also demonstrated. Notably, this new type of imaging does not require the measurement of the residual energies of the protons. METHODS: A large area, silicon microstrip tracker with high spatial and temporal resolution has been developed by the Proton Radiotherapy Verification and Dosimetry Applications consortium and commissioned using beams of protons at iThemba LABS, Medical Radiation Department, South Africa. The tracker comprises twelve planes of silicon developed using technology from high energy physics with each plane having an active area of ∼10 × 10 cm segmented into 2048 microstrips. The tracker is organized into four separate units each containing three detectors at 60° to one another creating an x-u-v coordinate system. Pairs of tracking units are used to reconstruct vertices for protons entering and exiting a phantom containing tissue equivalent inserts. By measuring the position and direction of each proton before and after the phantom, the nonlinear path for each proton through an object can be reconstructed. RESULTS: Experimental results are reported for tracking the path of protons with initial energies of 125 and 191 MeV. A spherical phantom of 75 mm diameter was imaged by positioning it between the entrance and exit detectors of the tracker. Positions and directions of individual protons were used to create angular distributions and 2D fluence maps of the beam. These results were acquired for 36 equally spaced projections spanning 180°, allowing, for the first time, an experimental CT image based upon the relative scattering power of protons to be reconstructed. CONCLUSIONS: Successful tracking of protons through a thick target (phantom) has demonstrated that the tracker discussed in this paper can provide the precise directional information needed to perform proton radiography and tomography. When synchronized with a range telescope, this could enable the reconstruction of proton CT images of stopping power. Furthermore, by measuring the deflection of many protons through a phantom, it was demonstrated that it is possible to reconstruct a new kind of CT image (scattering power) based upon this tracking information alone.

On the susceptibility of human platelets to aggregation by concanavalin A and the effect of this lectin on their response to ADP.

Concanavalin A aggregated gel-filtered platelets in 0.9% NaCl solution signifying cross-bridging by the lectin. Aggregation of these platelets by concanavalin A was temperature dependent; it did not occur at 0-4 degrees C unless the platelets were previously trypsinized. The level of aggregation of trypsinized platelets by concanavalin A at 0-4 degrees C was similar to that of untreated platelets at 37 degrees C. It is suggested that trypsin facilitates platelet aggregation by concanavalin A at 0-4 degrees C by causing a configurational change in membrane glycoproteins which orientates concanavalin A receptor sites into positions that favour lectin cross-bridging. Concanavalin A failed to aggregate platelets in plasma. Radioisotope studies showed that the amount of [3H]concanavalin A which combined with platelets in plasma was extremely low compared with gel-filtered platelets in saline. The aggregation of Ehrlich ascites cells by concanavalin A was considerably reduced when platelet-free plasma was added to the medium suggesting that it was due to the presence of concanavalin A-reactive components in the plasma. Concanavalin A inhibited the ADP-induced aggregation of platelets suspended in plasma or in a salts solution supplemented with calcium and fibrinogen, although the inhibitory effect was more conspicuous in the latter case. The results suggests that concanavalin A produces its inhibitory effect on ADP-induced platelet aggregation by interacting with membrane glycoproteins, and this further suggests their involvement in aggregation.