Validation of dual-energy CT-based composition analysis using fresh animal tissues and composition-optimized tissue equivalent samples.

Objective Proton therapy allows for highly conformal dose deposition, but is sensitive to range uncertainties. Several approaches currently under development measure composition-dependent secondary radiation to monitor the delivered proton range in-vivo. To fully utilize these methods, an estimate of the elemental composition of the patient's tissue is often needed. Approach A published dual-energy computed tomography (DECT)-based composition-extraction algorithm was validated against reference compositions obtained with two independent methods. For this purpose, a set of phantoms containing either fresh porcine tissue or tissue-mimicking samples with known, realistic compositions were imaged with a CT scanner at two different energies. Then, the prompt gamma-ray (PG) signal during proton irradiation was measured with a PG detector prototype. The PG workflow used pre-calculated Monte Carlo simulations to obtain an optimized estimate of the sample's carbon and oxygen contents. The compositions were also assessed with combustion analysis (CCA), and the stopping-power ratio (SPR) was measured with a multi-layer ionization chamber. The DECT images were used to calculate SPR-, density- and elemental composition maps, and to assign voxel-wise compositions from a selection of human tissues. For a more comprehensive set of reference compositions, the original selection was extended by 135 additional tissues, corresponding to spongiosa, high-density bones and low-density tissues. Results The root-mean-square error for the soft tissue carbon and oxygen content was 8.5wt% and 9.5wt% relative to the CCA result and 2.1wt% and 10.3wt% relative to the PG result. The phosphorous and calcium content were predicted within 0.4wt% and 1.1wt% of the CCA results, respectively. The largest discrepancies were encountered in samples whose composition deviated the most from tabulated compositions or that were more inhomogeneous. Significance Overall, DECT-based composition estimations of relevant elements were in equal or better agreement with the ground truth than the established SECT-approach and could contribute to in-vivo dose verification measurements.

Breast imaging with an ultra-low field MRI scanner: a pilot study.

Breast cancer screening is necessary to reduce mortality due to undetected breast cancer. Current methods have limitations, and as a result many women forego regular screening. Magnetic resonance imaging (MRI) can overcome most of these limitations, but access to conventional MRI is not widely available for routine annual screening. Here, we used an MRI scanner operating at ultra-low field (ULF) to image the left breasts of 11 women (mean age, 35 years ±13 years) in the prone position. Three breast radiologists reviewed the imaging and were able to discern the breast outline and distinguish fibroglandular tissue (FGT) from intramammary adipose tissue. Additionally, the expert readers agreed on their assessment of the breast tissue pattern including fatty, scattered FGT, heterogeneous FGT, and extreme FGT. This preliminary work demonstrates that ULF breast MRI is feasible and may be a potential option for comfortable, widely deployable, and low-cost breast cancer diagnosis and screening.

Global democratisation of proton radiotherapy.

Proton radiotherapy is an advanced treatment option compared with conventional x-ray treatment, delivering much lower doses of radiation to healthy tissues surrounding the tumour. However, proton therapy is currently not widely available. In this Review, we summarise the evolution of proton therapy to date, together with the benefits to patients and society. These developments have led to an exponential growth in the number of hospitals using proton radiotherapy worldwide. However, the gap between the number of patients who should be treated with proton radiotherapy and those who have access to it remains large. We summarise the ongoing research and development that is contributing to closing this gap, including the improvement of treatment efficiency and efficacy, and advances in fixed-beam treatments that do not require an enormously large, heavy, and costly gantry. The ultimate goal of decreasing the size of proton therapy machines to fit into standard treatment rooms appears to be within reach, and we discuss future research and development opportunities to achieve this goal.

Combined clinical and research training in medical physics in a multi-institutional setting: 13-year experience of Harvard Medical Physics Residency Program.

PURPOSE: This manuscript describes the structure, management and outcomes of a multi-institutional clinical and research medical physics residency program (Harvard Medical Physics Residency Program, or HMPRP) to provide potentially useful information to the centers considering a multi-institutional approach for their training programs. METHODS: Data from the program documents and public records was used to describe HMPRP and obtain statistics about participating faculty, enrolled residents, and graduates. Challenges associated with forming and managing a multi-institutional program and developed solutions for effective coordination between several clinical centers are described. RESULTS: HMPRP was formed in 2009 and was accredited by the Commission on Accreditation of Medical Physics Education Programs (CAMPEP) in 2011. It is a 3-year therapy program, with a dedicated year of research and the 2 years of clinical training at three academic hospitals. A CAMPEP-accredited Certificate Program is embedded in HMPRP to allow enrolled residents to complete a formal didactic training in medical physics if necessary. The clinical training covers the material required by CAMPEP. In addition, training in protons, CyberKnife, MR-linac, and at network locations is included. The clinical training and academic record of the residents is outstanding. All graduates have found employment within clinical medical physics, mostly at large academic centers and graduates had a 100% pass rate at the oral American Board of Radiology exams. On average, three manuscripts per resident are published during residency, and multiple abstracts are presented at conferences. CONCLUSIONS: A multi-institutional medical physics residency program can be successfully formed and managed. With a collaborative administrative structure, the program creates an environment for high-quality clinical training of the residents and high productivity in research. The main advantage of such program is access to a wide variety of resources. The main challenge is creating a structure for efficient management of multiple resources at different locations. This report may provide valuable information to centers considering starting a multi-institutional residency program.

Automated clinical target volume delineation using deep 3D neural networks in radiation therapy of Non-small Cell Lung Cancer.

BACKGROUND AND PURPOSE: Clinical targeted volume (CTV) delineation accounting for the patient-specific microscopic tumor spread can be a difficult step in defining the treatment volume. We developed an intelligent and automated CTV delineation system for locally advanced non-small cell lung carcinoma (NSCLC) to cover the microscopic tumor spread while avoiding organs-at-risk (OAR). MATERIALS AND METHODS: A 3D UNet with a customized loss function was used, which takes both the patients' respiration-correlated ("4D") CT scan and the physician contoured internal gross target volume (iGTV) as inputs, and outputs the CTV delineation. Among the 84 identified patients, 60 were randomly selected to train the network, and the remaining as testing. The model performance was evaluated and compared with cropped expansions using the shape similarities to the physicians' contours (the ground-truth) and the avoidance of critical OARs. RESULTS: On the testing datasets, all model-predicted CTV contours followed closely to the ground truth, and were acceptable by physicians. The average dice score was 0.86. Our model-generated contours demonstrated better agreement with the ground-truth than the cropped 5 mm/8 mm expansion method (median of median surface distance of 1.0 mm vs 1.9 mm/2.0 mm), with a small overlap volume with OARs (0.4 cm3 for the esophagus and 1.2 cm3 for the heart). CONCLUSIONS: The CTVs generated by our CTV delineation system agree with the physician's contours. This approach demonstrates the capability of intelligent volumetric expansions with the potential to be used in clinical practice.

Toward MR-integrated proton therapy: modeling the potential benefits for liver tumors.

Magnetic resonance imaging (MRI)-integrated proton therapy (MRiPT) is envisioned to improve treatment quality for many cancer patients. However, given the availability of alternative image-guided strategies, its clinical need is yet to be justified. This study aims to compare the expected clinical outcomes of MRiPT with standard of practice cone-beam CT (CBCT)-guided PT, and other MR-guided methods, i.e. offline MR-guided PT and MR-linac, for treatment of liver tumors. Clinical outcomes were assessed by quantifying the dosimetric and biological impact of target margin reduction enabled by each image-guided approach. Planning target volume (PTV) margins were calculated using random and systematic setup, delineation and motion uncertainties, which were quantified by analyzing longitudinal MRI data for 10 patients with liver tumors. Proton treatment plans were created using appropriate PTV margins for each image-guided PT method. Photon plans with margins equivalent to MRiPT were generated to represent MR-linac. Normal tissue complication probabilities (NTCP) of the uninvolved liver were compared. We found that PTV margin can be reduced by 20% and 40% for offline MR-guided PT and MRiPT, respectively, compared with CBCT-guided PT. Furthermore, clinical target volume expansion could be largely alleviated when delineating on MRI rather than CT. Dosimetric implications included decreased equivalent mean dose of the uninvolved liver, i.e. up to 24.4 Gy and 27.3 Gy for offline MR-guided PT and MRiPT compared to CBCT-guided PT, respectively. Considering Child-Pugh score increase as endpoint, NTCP of the uninvolved liver was significantly decreased for MRiPT compared to CBCT-guided PT (up to 48.4%,p < 0.01), offline MR-guided PT (up to 12.9%,p < 0.01) and MR-linac (up to 30.8%,p < 0.05). Target underdose was possible in the absence of MRI-guidance (D90 reduction up to 4.2 Gy in 20% of cases). In conclusion, MRiPT has the potential to significantly reduce healthy liver toxicities in patients with liver tumors. It is superior to other image-guided techniques currently available.

Inter-patient variations of radiation-induced normal-tissue changes in Gd-EOB-DTPA-enhanced hepatic MRI scans during fractionated proton therapy.

BACKGROUND AND PURPOSE: Previous MRI studies have shown a substantial decrease in normal-tissue uptake of a hepatobiliary-directed contrast agent 6-9 weeks after liver irradiation. In this prospective clinical study, we investigated whether this effect is detectable during the course of proton therapy. MATERIAL AND METHODS: Gd-EOB-DTPA enhanced MRI was performed twice during hypo-fractionated proton therapy of liver lesions in 9 patients (plus two patients with only one scan available). Dose-correlated signal changes were qualitatively scored based on difference images from the two scans. We evaluated the correlation between the MRI signal change with the planned dose map. The GTV was excluded from all analyses. In addition, were examined timing, irradiated liver volume, changes in liver function parameters as well as circulating biomarkers of inflammation. RESULTS: Strong, moderate or no dose-correlated signal changes were detected for 2, 3 and 5 patients, respectively. Qualitative scoring was consistent with the quantitative dose to signal change correlation. In an exploratory analysis, the strongest correlation was found between the qualitative scoring and pretreatment IL-6 concentration. For all patients, a clear dose-correlated signal decrease was seen in late follow-up scans. CONCLUSION: Radiation-induced effects can be detected with Gd-EOB-DTPA enhanced MRI in a subgroup of patients within a few days after proton irradiation. The reason for the large inter-patient variations is not yet understood and will require validation in larger studies.

Feasibility study of in vivo MRI based dosimetric verification of proton end-of-range for liver cancer patients.

PURPOSE: To investigate the feasibility of using MRI to verify proton beam distal range for liver tumor treatment in a retrospective study. METHODS AND MATERIALS: Because the follow-up hepatocyte-specific functional MR imaging can detect the radiobiological change of liver tissue after radiation, we firstly registered the contrast-enhanced MR images to the planning CT images from 5 liver patients, then overlaid the prescribed dose distribution on the MR images. Since dose calculation is most accurate at the penumbra dose region, we correlated the MR signal intensity (SI) to the radiation dose at the superior/inferior penumbra region. This dose-SI correlation was finally employed on registered MR images to estimate the proton end-of-range. RESULTS: Statistically significant correlations between radiation dose and MR SI were observed in superior/inferior penumbra regions, with correlation coefficient ranging from 0.93 to 0.99. By applying the dose-SI correlation to the distal region of each proton beam, the mean difference between MR-estimated and the planned dose range was -2.18 ± 4.89 mm for anterior-posterior beams and -3.90 ± 5.87 mm for lateral beams. CONCLUSIONS: This feasibility study proved the principle that proton dose range can be verified in vivo by follow-up MR images after proton liver treatment.

Clinical application of in-room positron emission tomography for in vivo treatment monitoring in proton radiation therapy.

PURPOSE: The purpose of this study is to evaluate the potential of using in-room positron emission tomography (PET) for treatment verification in proton therapy and for deriving suitable PET scan times. METHODS AND MATERIALS: Nine patients undergoing passive scattering proton therapy underwent scanning immediately after treatment with an in-room PET scanner. The scanner was positioned next to the treatment head after treatment. The Monte Carlo (MC) method was used to reproduce PET activities for each patient. To assess the proton beam range uncertainty, we designed a novel concept in which the measured PET activity surface distal to the target at the end of range was compared with MC predictions. The repositioning of patients for the PET scan took, on average, approximately 2 minutes. The PET images were reconstructed considering varying scan times to test the scan time dependency of the method. RESULTS: The measured PET images show overall good spatial correlations with MC predictions. Some discrepancies could be attributed to uncertainties in the local elemental composition and biological washout. For 8 patients treated with a single field, the average range differences between PET measurements and computed tomography (CT) image-based MC results were <5 mm (<3 mm for 6 of 8 patients) and root-mean-square deviations were 4 to 11 mm with PET-CT image co-registration errors of approximately 2 mm. Our results also show that a short-length PET scan of 5 minutes can yield results similar to those of a 20-minute PET scan. CONCLUSIONS: Our first clinical trials in 9 patients using an in-room PET system demonstrated its potential for in vivo treatment monitoring in proton therapy. For a quantitative range prediction with arbitrary shape of target volume, we suggest using the distal PET activity surface.

Improved planning time and plan quality through multicriteria optimization for intensity-modulated radiotherapy.

PURPOSE: To test whether multicriteria optimization (MCO) can reduce treatment planning time and improve plan quality in intensity-modulated radiotherapy (IMRT). METHODS AND MATERIALS: Ten IMRT patients (5 with glioblastoma and 5 with locally advanced pancreatic cancers) were logged during the standard treatment planning procedure currently in use at Massachusetts General Hospital (MGH). Planning durations and other relevant planning information were recorded. In parallel, the patients were planned using an MCO planning system, and similar planning time data were collected. The patients were treated with the standard plan, but each MCO plan was also approved by the physicians. Plans were then blindly reviewed 3 weeks after planning by the treating physician. RESULTS: In all cases, the treatment planning time was vastly shorter for the MCO planning (average MCO treatment planning time was 12 min; average standard planning time was 135 min). The physician involvement time in the planning process increased from an average of 4.8 min for the standard process to 8.6 min for the MCO process. In all cases, the MCO plan was blindly identified as the superior plan. CONCLUSIONS: This provides the first concrete evidence that MCO-based planning is superior in terms of both planning efficiency and dose distribution quality compared with the current trial and error-based IMRT planning approach.

Monitoring proton radiation therapy with in-room PET imaging.

We used a mobile positron emission tomography (PET) scanner positioned within the proton therapy treatment room to study the feasibility of proton range verification with an in-room, stand-alone PET system, and compared with off-line equivalent studies. Two subjects with adenoid cystic carcinoma were enrolled into a pilot study in which in-room PET scans were acquired in list-mode after a routine fractionated treatment session. The list-mode PET data were reconstructed with different time schemes to generate in-room short, in-room long and off-line equivalent (by skipping coincidences from the first 15 min during the list-mode reconstruction) PET images for comparison in activity distribution patterns. A phantom study was followed to evaluate the accuracy of range verification for different reconstruction time schemes quantitatively. The in-room PET has a higher sensitivity compared to the off-line modality so that the PET acquisition time can be greatly reduced from 30 to <5 min. Features in deep-site, soft-tissue regions were better retained with in-room short PET acquisitions because of the collection of (15)O component and lower biological washout. For soft tissue-equivalent material, the distal fall-off edge of an in-room short acquisition is deeper compared to an off-line equivalent scan, indicating a better coverage of the high-dose end of the beam. In-room PET is a promising low cost, high sensitivity modality for the in vivo verification of proton therapy. Better accuracy in Monte Carlo predictions, especially for biological decay modeling, is necessary.

The role of medical physicists and the AAPM in the development of treatment planning and optimization.

Developments in radiotherapy treatment planning and optimization by medical physicists and the American Association of Physicists in Medicine are reviewed, with emphasis on recent work in optimization. It is shown that medical physicists have played a vital role in the creation of innovative treatment planning techniques throughout the past century, most significantly since the advent of computerized tomography for three-dimensional (3D) imaging and high-powered computers capable of 3D planning and optimization. Some early advances in 3D planning made by physicists include development of novel planning algorithms, beam's-eye-view, virtual simulation, dose-volume histogram analysis tools, and bioeffect modeling. Most of the recent developments have been driven by the need to develop treatment planning for conformal radiotherapy, especially intensity modulated radiation therapy. These advances include inverse planning, handling the effects of motion and uncertainty, biological planning, and multicriteria optimization.

Multicriteria optimization in intensity-modulated radiation therapy treatment planning for locally advanced cancer of the pancreatic head.

PURPOSE: Intensity-modulated radiation therapy (IMRT) affords the potential to decrease radiation therapy-associated toxicity by creating highly conformal dose distributions. However, the inverse planning process can create a suboptimal plan despite meeting all constraints. Multicriteria optimization (MCO) may reduce the time-consuming iteration loop necessary to develop a satisfactory plan while providing information regarding trade-offs between different treatment planning goals. In this exploratory study, we examine the feasibility and utility of MCO in physician plan selection in patients with locally advanced pancreatic cancer (LAPC). METHODS AND MATERIALS: The first 10 consecutive patients with LAPC treated with IMRT were evaluated. A database of plans (Pareto surface) was created that met the inverse planning goals. The physician then navigated to an "optimal" plan from the point on the Pareto surface at which kidney dose was minimized. RESULTS: Pareto surfaces were created for all 10 patients. A physician was able to select a plan from the Pareto surface within 10 minutes for all cases. Compared with the original (treated) IMRT plans, the plan selected from the Pareto surface had a lower stomach mean dose in 9 of 10 patients, although often at the expense of higher kidney dose than with the treated plan. CONCLUSION: The MCO is feasible in patients with LAPC and allows the physician to choose a satisfactory plan quickly. Generally, when given the opportunity, the physician will choose a plan with a lower stomach dose. The MCO enables a physician to provide greater active clinical input into the IMRT planning process.

Accelerating IMRT optimization by voxel sampling.

This paper presents a new method for accelerating intensity-modulated radiation therapy (IMRT) optimization using voxel sampling. Rather than calculating the dose to the entire patient at each step in the optimization, the dose is only calculated for some randomly selected voxels. Those voxels are then used to calculate estimates of the objective and gradient which are used in a randomized version of a steepest descent algorithm. By selecting different voxels on each step, we are able to find an optimal solution to the full problem. We also present an algorithm to automatically choose the best sampling rate for each structure within the patient during the optimization. Seeking further improvements, we experimented with several other gradient-based optimization algorithms and found that the delta-bar-delta algorithm performs well despite the randomness. Overall, we were able to achieve approximately an order of magnitude speedup on our test case as compared to steepest descent.

Radiotherapy treatment of early-stage prostate cancer with IMRT and protons: a treatment planning comparison.

PURPOSE: To compare intensity-modulated photon radiotherapy (IMRT) with three-dimensional conformal proton therapy (3D-CPT) for early-stage prostate cancer, and explore the potential utility of intensity-modulated proton therapy (IMPT). METHODS AND MATERIALS: Ten patients were planned with both 3D-CPT (two parallel-opposed lateral fields) and IMRT (seven equally spaced coplanar fields). Prescribed dose was 79.2 Gy (or cobalt Gray-equivalent, [CGE] for protons) to the prostate gland. Dose-volume histograms, dose conformity, and equivalent uniform dose (EUD) were compared. Additionally, plans were optimized for 3D-CPT with nonstandard beam configuration, and for IMPT assuming delivery with beam scanning. RESULTS: At least 98% of the planning target volume received the prescription dose. IMRT plans yielded better dose conformity to the target, whereas proton plans achieved higher dose homogeneity and better sparing of rectum and bladder in the range below 30 Gy/CGE. Bladder volumes receiving more than 70 Gy/CGE (V70) were reduced, on average, by 34% with IMRT vs. 3D-CPT, whereas rectal V70 were equivalent. EUD from 3D-CPT and IMRT plans were indistinguishable within uncertainties for both bladder and rectum. With the use of small-angle lateral-oblique fields in 3D-CPT and IMPT, the rectal V70 was reduced by up to 35% compared with the standard lateral configuration, whereas the bladder V70 increased by less than 10%. CONCLUSIONS: In the range higher than 60 Gy/CGE, IMRT achieved significantly better sparing of the bladder, whereas rectal sparing was similar with 3D-CPT and IMRT. Dose to healthy tissues in the range lower than 50% of the target prescription was substantially lower with proton therapy.

Approximating convex pareto surfaces in multiobjective radiotherapy planning.

Radiotherapy planning involves inherent tradeoffs: the primary mission, to treat the tumor with a high, uniform dose, is in conflict with normal tissue sparing. We seek to understand these tradeoffs on a case-to-case basis, by computing for each patient a database of Pareto optimal plans. A treatment plan is Pareto optimal if there does not exist another plan which is better in every measurable dimension. The set of all such plans is called the Pareto optimal surface. This article presents an algorithm for computing well distributed points on the (convex) Pareto optimal surface of a multiobjective programming problem. The algorithm is applied to intensity-modulated radiation therapy inverse planning problems, and results of a prostate case and a skull base case are presented, in three and four dimensions, investigating tradeoffs between tumor coverage and critical organ sparing.