Hu, Yanle

PurposeA magnetic resonance (MR) biologic marker (biomarker) is a measurable quantitative characteristic that is an indicator of normal biological and pathogenetic processes or a response to therapeutic intervention derived from the MR imaging process. There is significant potential for MR biomarkers to facilitate personalized approaches to cancer care through more precise disease targeting by quantifying normal versus pathologic tissue function as well as toxicity to both radiation and chemotherapy. Both of which have the potential to increase the therapeutic ratio and provide earlier, more accurate monitoring of treatment response. The ongoing integration of MR into routine clinical radiation therapy (RT) planning and the development of MR guided radiation therapy systems is providing new opportunities for MR biomarkers to personalize and improve clinical outcomes. Their appropriate use, however, must be based on knowledge of the physical origin of the biomarker signal, the relationship to the underlying biological processes, and their strengths and limitations. The purpose of this report is to provide an educational resource describing MR biomarkers, the techniques used to quantify them, their strengths and weakness within the context of their application to radiation oncology so as to ensure their appropriate use and application within this field.

AbstractPurposeTo demonstrate that variation in chemical composition has a negligible effect on the mapping curve from relative electron density (RED) to proton stopping power ratio (SPR), and to establish the theoretical framework of using Megavoltage (MV) computed tomography (CT), instead of kilovoltage (kV) dual energy CT, to accurately estimate proton SPR.MethodsA simulation study was performed to evaluate the effect of chemical composition variation on kVCT number and proton SPR. The simulation study involved both reference and simulated human tissues. The reference human tissues, together with their physical densities and chemical compositions, came from the ICRP publication 23. The simulated human tissues were created from the reference human tissues assuming that elemental percentage weight followed a Gaussian distribution. For all tissues, kVCT number and proton SPR were obtained through (a) theoretical calculation from tissue’s physical density and chemical composition which served as the ground truth, and (b) estimation from RED using the calibration curves established from the stoichiometric method. Deviations of the estimated values from the calculated values were quantified as errors in using RED to estimate kVCT number and proton SPR.ResultsGiven a chemical composition variation of 5% (1σ) of the nominal percentage weights, the total estimation error of using RED to estimate kVCT number was 0.34%, 0.62%, and 0.77% and the total estimation error of using RED to estimate proton SPR was 0.30%, 0.22%, and 0.16% for fat tissues, non‐fat soft tissues and bone tissues, respectively.ConclusionChemical composition had a negligible effect on the method of using RED to determine proton SPR. RED itself is sufficient to accurately determine proton SPR. MVCT number maintains a superb linear relationship with RED because it is highly dominated by Compton scattering. Therefore, MVCT has great potential in reducing the proton range uncertainty.

The use of dedicated magnetic resonance simulation (MR‐SIM) platforms in Radiation Oncology has expanded rapidly, introducing new equipment and functionality with the overall goal of improving the accuracy of radiation treatment planning. However, this emerging technology presents a new set of challenges that need to be addressed for safe and effective MR‐SIM implementation. The major objectives of this report are to provide recommendations for commercially available MR simulators, including initial equipment selection, siting, acceptance testing, quality assurance, optimization of dedicated radiation therapy specific MR‐SIM workflows, patient‐specific considerations, safety, and staffing. Major contributions include guidance on motion and distortion management as well as MRI coil configurations to accommodate patients immobilized in the treatment position. Examples of optimized protocols and checklists for QA programs are provided. While the recommendations provided here are minimum requirements, emerging areas and unmet needs are also highlighted for future development.

PurposeTo clinically implement and comprehensively evaluate two independent methods for beam monitor calibration of scanning proton beam.MethodsSeven proton energies that represent the lowest to highest energy proton beams were selected. Single energy layer circular fields of diameter 15 cm with 2.5 mm spot spacing and 10 times of repainting (FS15cm) were designed for all energies. The effective measurement points of Bragg peak chamber (BPC), advanced Markus chamber (AMC) and farmer chamber (FC) were all aligned to 2 cm depth in water using SSD setup. The BPC and AMC were cross‐calibrated with farmer chamber (FC) using the field FS15cm. In order to evaluate BPC's lateral response uniformity, a collimated narrow proton beam (5.8 mm diameter) was delivered to the active area and edge of the BPC. The dose area product (DAP) was measured using two methods by two BPCs, one AMC and one FC. For method 1, a single spot proton beam was delivered to the geometric center of the BPC. For method 2, the fields FS15cm were delivered to FC and AMC, respectively. Accumulated charges by these chambers were converted to DAPs, and the quantitative difference of DAPs between both methods was calculated. The causes of the uncertainties were discussed, and the advantages of the two methods were compared.ResultsThe two BPCs showed different lateral response uniformity. BPC1 has a uniform response from the center up to a radius of 3.5 cm. BPC2 has a uniform response only to 2 cm and the response dropped 1% to 2% at 3.5 cm from center. BPC2 also has significant over‐response compared to BPC1. A 2.2% systematic error would be transferred to DAP if the over‐response from BPC2 was not considered. The DAPs measured by method 1 with two BPCs and by method 2 with FC and AMC were consistent to 0.5%. The major uncertainty component of method 1 is from the cross‐calibration of the BPC.ConclusionsThe two independent methods for DAP were shown to give consistent results, given the sources of uncertainties were carefully addressed in the measurements. Direct measurement of DAP with BPC is very efficient, but it may be subject to more than 2% systematic error if the BPC lateral response is not carefully evaluated.

PurposeApproximate dose calculation methods were used in the nominal dose distribution and the perturbed dose distributions due to uncertainties in a commercial treatment planning system (CTPS) for robust optimization in intensity‐modulated proton therapy (IMPT). We aimed to investigate whether the approximations influence plan quality, robustness, and interplay effect of the resulting IMPT plans for the treatment of locally advanced lung cancer patients.Materials and methodsTen consecutively treated locally advanced nonsmall cell lung cancer (NSCLC) patients were selected. Two IMPT plans were created for each patient using our in‐house developed TPS, named “Solo,” and also the CTPS, EclipseTM (Varian Medical Systems, Palo Alto, CA, USA), respectively. The plans were designed to deliver prescription doses to internal target volumes (ITV) drawn by a physician on averaged four‐dimensional computed tomography (4D‐CT). Solo plans were imported back to CTPS, and recalculated in CTPS for fair comparison. Both plans were further verified for each patient by recalculating doses in the inhalation and exhalation phases to ensure that all plans met clinical requirements. Plan robustness was quantified on all phases using dose‐volume‐histograms (DVH) indices in the worst‐case scenario. The interplay effect was evaluated for every plan using an in‐house developed software, which randomized starting phases of each field per fraction and accumulated dose in the exhalation phase based on the patient's breathing motion pattern and the proton spot delivery in a time‐dependent fashion. DVH indices were compared using Wilcoxon rank‐sum test.ResultsCompared to the plans generated using CTPS on the averaged CT, Solo plans had significantly better target dose coverage and homogeneity (normalized by the prescription dose) in the worst‐case scenario [ITV D95%: 98.04% vs 96.28%, Solo vs CTPS, P = 0.020; ITV D5%–D95%: 7.20% vs 9.03%, P = 0.049] while all DVH indices were comparable in the nominal scenario. On the inhalation phase, Solo plans had better target dose coverage and cord Dmax in the nominal scenario [ITV D95%: 99.36% vs 98.45%, Solo vs CTPS, P = 0.014; cord Dmax: 20.07 vs 23.71 Gy(RBE), P = 0.027] with better target coverage and cord Dmax in the worst‐case scenario [ITV D95%: 97.89% vs 96.47%, Solo vs CTPS, P = 0.037; cord Dmax: 24.57 vs 28.14 Gy(RBE), P = 0.037]. On the exhalation phase, similar phenomena were observed in the nominal scenario [ITV D95%: 99.63% vs 98.87%, Solo vs CTPS, P = 0.037; cord Dmax: 19.67 vs 23.66 Gy(RBE), P = 0.039] and in the worst‐case scenario [ITV D95%: 98.20% vs 96.74%, Solo vs CTPS, P = 0.027; cord Dmax: 23.47 vs 27.93 Gy(RBE), P = 0.027]. In terms of interplay effect, plans generated by Solo had significantly better target dose coverage and homogeneity, less hot spots, and lower esophageal Dmean, and cord Dmax [ITV D95%: 101.81% vs 98.68%, Solo vs CTPS, P = 0.002; ITV D5%–D95%: 2.94% vs 7.51%, P = 0.002; cord Dmax: 18.87 vs 22.29 Gy(RBE), P = 0.014].ConclusionsSolo‐generated IMPT plans provide improved cord sparing, better target robustness in all considered phases, and reduced interplay effect compared with CTPS. Consequently, the approximation methods currently used in commercial TPS programs may have space for improvement in generating optimal IMPT plans for patient cases with locally advanced lung cancer.

BackgroundEsophageal carcinoma is the eighth most common cancer in the world. Volumetric‐modulated arc therapy (VMAT) is widely used to treat distal esophageal carcinoma due to high conformality to the target and good sparing of organs at risk (OAR). It is not clear if small‐spot intensity‐modulated proton therapy (IMPT) demonstrates a dosimetric advantage over VMAT. In this study, we compared dosimetric performance of VMAT and small‐spot IMPT for distal esophageal carcinoma in terms of plan quality, plan robustness, and interplay effects.Methods35 distal esophageal carcinoma patients were retrospectively reviewed; 19 patients received small‐spot IMPT and the remaining 16 of them received VMAT. Both plans were generated by delivering prescription doses to clinical target volumes (CTVs) on phase‐averaged 4D‐CT's. The dose‐volume‐histogram (DVH) band method was used to quantify plan robustness. Software was developed to evaluate interplay effects with randomized starting phases for each field per fraction. DVH indices were compared using Wilcoxon rank‐sum test. For fair comparison, all the treatment plans were normalized to have the same CTVhigh D95% in the nominal scenario relative to the prescription dose.ResultsIn the nominal scenario, small‐spot IMPT delivered statistically significantly lower liver Dmean and V30Gy[RBE], lung Dmean, heart Dmean compared with VMAT. CTVhigh dose homogeneity and protection of other OARs were comparable between the two treatments. In terms of plan robustness, the IMPT and VMAT plans were comparable for kidney V18Gy[RBE], liver V30Gy[RBE], stomach V45Gy[RBE], lung Dmean, V5Gy[RBE], and V20Gy[RBE], cord Dmax and , liver Dmean, heart V20Gy[RBE], and V30Gy[RBE], but IMPT was significantly worse for CTVhigh D95%, , and D5%‐D95%, CTVlow D95%, heart Dmean, and V40Gy[RBE], requiring careful and experienced adjustments during the planning process and robustness considerations. The small‐spot IMPT plans still met the standard clinical requirements after interplay effects were considered.ConclusionsSmall‐spot IMPT decreases doses to heart, liver, and total lung compared to VMAT as well as achieves clinically acceptable plan robustness. Our study supports the use of small‐spot IMPT for the treatment of distal esophageal carcinoma.

PurposeBiomechanical properties can be used as biomarkers to diagnose tumors, monitor tumor development, and evaluate treatment efficacy. The purpose of this preliminary study is to characterize the biomechanical environment of two typical liver tumors, hemangiomas (HEMs) and hepatocellular carcinomas (HCCs), and to investigate the potential of using strain metrics as biomarkers for tumor diagnosis, based on a limited clinical dataset.MethodsMagnetic resonance (MR) tagging was used to quantify the motion and deformation of the two types of liver tumors. Displacements of the tumors arising from a heartbeat were measured over one cardiac cycle. Local biomechanical conditions of the tumors were characterized by estimating two principal strains (ε1 and ε2) and an octahedral shear strain (εsoct) of the tumor and its peripheral region. Biomechanical conditions of the tumors were compared with those of the arbitrarily selected regions from healthy volunteers.ResultsWe observed that the HCCs had significantly smaller strain values compared to their peripheral tissues. However, the HEMs did not have significantly different strains from those of the peripheral tissues, and were similar to healthy liver regions. The sensitivity of using ε1, ε2, and εsoct to diagnose HCC were all 1, while the sensitivity of using ε1, ε2, and εsoct to diagnose HEM were 0.67, 0.17, and 0.67, respectively.ConclusionsLagrangian strain metrics provide insight into the biomechanical conditions of certain liver tumors in the human body and may provide another perspective for tumor characterization and diagnosis.

PurposeTo establish a mathematical model to guide prospective T2‐weighted four‐dimensional magnetic resonance imaging (4DMRI) acquisition and to propose an efficient solution to expedite prospective T2‐weighted 4DMRI acquisition.MethodsProspective T2‐weighted 4DMRI acquisition was characterized by a mathematical model with 4DMRI acquisition time as the objective function and completeness of the image set, acquisition timing, image contrast, and image artifacts as constraints. Given the irregular nature of human respiration, an efficient solution based on the greedy strategy (ESGS) was proposed. The efficiency of the ESGS method was validated using healthy human subjects. Comparisons were made with the previous 4DMRI method incorporating the prefixed‐order respiratory state splitting (PO‐RSS) technique.Results4DMRI image sets acquired using the ESGS and PO‐RSS methods had similar image quality. The average time to acquire a 4DMRI image set covering 60 slices at 10 respiratory states was reduced by 30%, from 13.1 min using the PO‐RSS method to 9.0 min using the ESGS method. It was demonstrated that high‐quality T2‐weighted 4DMRI could be obtained within a reasonable amount of time and all slices within each of the three‐dimensional volumes were indeed acquired at the same respiratory state.ConclusionsThe ESGS method substantially reduces the acquisition time for T2‐weighted 4DMRI, making it ready for clinical evaluation to obtain abdominal tumor motion for radiotherapy treatment planning.

AbstractMonte Carlo (MC) simulation has been used to generate commissioning data for the beam modeling of treatment planning system (TPS). We have developed a method called radial projection (RP) for postprocessing of MC‐simulation‐generated data. We used the RP method to reduce the statistical uncertainty of the lateral profile of proton pencil beams with axial symmetry. The RP method takes advantage of the axial symmetry of dose distribution to use the mean value of multiple independent scores as the representative score. Using the mean as the representative value rather than any individual score results in substantial reduction in statistical uncertainty. Herein, we present the concept and step‐by‐step implementation of the RP method, as well as show the advantage of the RP method over conventional measurement methods for generating lateral profile. Lateral profiles generated by both methods were compared to demonstrate the uncertainty reduction qualitatively, and standard error comparison was performed to demonstrate the reduction quantitatively. The comparisons showed that statistical uncertainty was reduced substantially by the RP method. Using the RP method to postprocess MC data, the corresponding MC simulation time was reduced by a factor of 10 without quality reduction in the generated result from the MC data. We concluded that the RP method is an effective technique to increase MC simulation efficiency for generating lateral profiles for axially symmetric pencil beams.