TY - JOUR
AU1 - Rivera, Judith N.
AU2 - Kierski, Thomas M.
AU3 - Kasoji, Sandeep K.
AU4 - Abrantes, Anthony S.
AU5 - Dayton, Paul A.
AU6 - Chang, Sha X.
AB - Purpose To identify key dosimetric parameters that have close associations with tumor treatment response and body weight change in SFRT treatments with a large range of spatial-fractionation scale at dose rates of several Gy/min. Methods Six study arms using uniform tumor radiation, half-tumor radiation, 2mm beam array radiation, 0.3mm minibeam radiation, and an untreated arm were used. All treatments were delivered on a 320kV x-ray irradiator. Forty-two female Fischer 344 rats with fibrosarcoma tumor allografts were used. Dosimetric parameters studied are peak dose and width, valley dose and width, peak-to-valley-dose-ratio (PVDR), volumetric average dose, percentage volume directly irradiated, and tumor- and normal-tissue EUD. Animal survival, tumor volume change, and body weight change (indicative of treatment toxicity) are tested for association with the dosimetric parameters using linear regression and Cox Proportional Hazards models. Results The dosimetric parameters most closely associated with tumor response are tumor EUD (R2 = 0.7923, F-stat = 15.26*; z-test = -4.07***), valley (minimum) dose (R2 = 0.7636, F-stat = 12.92*; z-test = -4.338***), and percentage tumor directly irradiated (R2 = 0.7153, F-stat = 10.05*; z-test = -3.837***) per the linear regression and Cox Proportional Hazards models, respectively. Tumor response is linearly proportional to valley (minimum) doses and tumor EUD. Average dose (R2 = 0.2745, F-stat = 1.514 (no sig.); z-test = -2.811**) and peak dose (R2 = 0.04472, F-stat = 0.6874 (not sig.); z-test = -0.786 (not sig.)) show the weakest associations to tumor response. Only the uniform radiation arm did not gain body weight post-radiation, indicative of treatment toxicity; however, body weight change in general shows weak association with all dosimetric parameters except for valley (minimum) dose (R2 = 0.3814, F-stat = 13.56**), valley width (R2 = 0.2853, F-stat = 8.783**), and peak width (R2 = 0.2759, F-stat = 8.382**). Conclusions For a single-fraction SFRT at conventional dose rates, valley, not peak, dose is closely associated with tumor treatment response and thus should be used for treatment prescription. Tumor EUD, valley (minimum) dose, and percentage tumor directly irradiated are the top three dosimetric parameters that exhibited close associations with tumor response. Introduction Spatially-fractionated radiation therapy (SFRT) is a nonconventional radiation therapy that is characterized by intentionally-created high dose inhomogeneities, ultra-high maximum doses, and single fraction treatments [1, 2]. The dose inhomogeneity consists of many small sub-regions that vary in pattern, shape, size and with alternating high and low doses throughout the treatment volume. For example, in clinical GRID Therapy, which is part of the SFRT family, [1, 3] the high dose sub-regions are 1cm-wide spaced ~2cm apart in a hexagonal pattern, while in preclinical microbeam radiation therapy (MRT) [4] the sub-regions are 10’s to 100s microns wide micro-planar arrays of radiation, spaced 100–4000 microns apart. Both the clinical and preclinical forms of SFRT have a decades-long history demonstrating a superior therapeutic ratio compared to conventional radiation therapy, especially in terms of normal organ sparing. Detailed summaries can be found in two recent reviews by Billena and Khan [5] for GRID therapy and by Eling et al. [6, 7] for MRT. Today, there are a number of modern treatment delivery technologies available for clinical SFRT including multi-leaf collimator generated [8–11], Tomotherapy [5], and particle GRID therapy [12, 13]. For preclinical SFRT, newer technologies include “minibeams” with larger spatial fractionation scales (on the order of millimeter instead of the tens of microns used in classical MRT [14, 15] and with conventional dose-rates [16, 17], [4]. The conventional dose rate SFRT radiations, such as the ones used in this study, are highly relevant to translational research for LINAC-based SFRT clinical applications, where conventional dose rates are also used. Despite the long history and well demonstrated therapeutic ratio advantage over conventional uniform dose radiation therapy, SFRT remains an experimental therapy. There are several reasons attributed to the sluggish clinical translation progress including a lack of understanding of SFRT working mechanisms and of the association between SFRT treatment response and dosimetry. While we have verified treatment dosimetry and tumor control outcome correlations for conventional radiation therapy (i.e., tumor minimum dose and Equivalent Uniform Dose (EUD) are closely correlated with tumor control) [18] we do not yet have such understanding for SFRT, which has significantly more complex dosimetry than that of conventional radiation therapy. Unique SFRT dosimetric parameters that describe the dosimetry include peak dose, valley dose, peak-to-valley-dose-ratio (PVDR), peak width, valley width, and percentage tumor volume directly irradiated. It is reasonable to assume that not all these dosimetric parameters have the same clinical significance. To effectively advance SFRT clinical translation it is critically important to identify which parameters have strong/weak associations with a given treatment response. The goal of this study is to identify key dosimetric parameters that are most closely associated with treatment response using a preclinical animal model. We hypothesize that while peak dose has always been used to prescribe SFRT treatment for both clinical and preclinical applications, peak dose may not be the dosimetric parameter most closely associated with SFRT tumor control or treatment toxicity. If it is not, which SFRT dosimetric parameters are? Further, we ask that, for a given pattern of SFRT treatment, what is its conventional radiation therapy equivalence for a given treatment response? These questions are crucial to advance clinical translation of SFRT; however, decades of synchrotron-based MRT studies may not be well-suited to answer these questions due to the use of ultrahigh dose rates (1000s Gy/sec) [19]. Recent research on ultra-high dose-rate FLASH radiation has shown that radiation with dose-rates of 100 Gy/s or higher selectively spares normal tissue not tumor [4, 20, 21]. This finding revealed that the ultrahigh dose-rate alone is likely to be partially responsible for the observed high therapeutic-ratio demonstrated in the majority of SFRT research published so far [6]. This study is designed to discern the impact of radiation spatial fractionation on treatment response and the association between dosimetric parameters and treatment response at dose rates relevant to clinical SFRT treatments. Today, SFRT is receiving much deserved renewed attention and enthusiasm in the field of radiation oncology. In 2018 National Cancer Institute and Radiosurgery Society jointly held the first workshop on Understanding High-Dose, Ultra-Dose-Rate and Spatially Fractionated Radiotherapy and created three standing working groups (clinical, biology, and physics) aiming to provide guidelines on SFRT research and clinical application [22]. We hope this work will assist in this endeavor by shedding light on the clinical impact of different SFRT dosimetry parameters. Materials and methods Study design The secret of SFRT lays in its radiation dose spatial fractionation. Although this work does not address the very much needed understanding of working mechanism it addresses another important matter for SFRT application—the association of SFRT dosimetric parameters with treatment response at conventional dose rates (dose rate ranges from 4.27 to 5.25 Gy/min was used). Fig 1 shows a six-arm study design using a very large span of radiation spatial fractionation, constructed to explore the impact of radiation spatial fractionation. Table 1 summarizes the dosimetric parameters of each of the six arms. To study the effect of radiation spatial fractionation under the condition of equal volume-averaged dose we used the following four study arms: 20GyUniformRT (entire tumor directly irradiated), 20GyHalfRT (only one-half of tumor directly irradiated), 20Gy2mmSFRT (50% of tumor directly irradiated by 2mm-wide planar beam array), and 20GySFRT (20% of tumor directly irradiated with 0.3mm-wide planar beam array). Note that the doses are volume-averaged doses computed for the entire tumor volume. A 50GySFRT arm (50 Gy volume-averaged dose, beam width 0.31mm) is added as it has a peak dose of 225 Gy, which is within the known minibeam peak dose range showing tumor control. To account for unavoidable variations in tumor position under the 20Gy2mmSFRT treatment beams during animal irradiations, we computed the maximum and minimum beam coverage positions and calculated their corresponding dosimetric specifications. The 20Gy2mmSFRT treatment arm dosimetric values reported in Table 1 correspond to the average at these positions for a 10mm diameter tumor. For example, a 10mm sized tumor is irradiated by at most three 2mm-peaks and at minimum two 2mm-peaks and the average dosimetric parameter at these two positions was calculated. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 1. Illustration of the SFRT spatial fractionation study design. A very large range of radiation spatial fractionation scale was used to derive the impact of radiation spatial fractionation. Four arms share the same 20 Gy volume-average dose. The high dose 50GySFRT arm is added because 20GySFRT is not known to have tumor control. The dosimetric parameters studied and number of animals per study arm are listed in Table 1. https://doi.org/10.1371/journal.pone.0229053.g001 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 1. Summary of nine SFRT dosimetric parameter specifications in the six-arm study. https://doi.org/10.1371/journal.pone.0229053.t001 Custom-made radiation blocks and collimators made of Cerrobend or tungsten were used to define the 2cmx2cm field for 20GyUniformRT arm treatment, the 2cmx1cm for 20GyHalfSFRT treatment, and the beamlet array 2cmx2cm fields for both the 20Gy2mmSFRT and 20GySFRT/50GySFRT treatments. The 2cm field size in the direction of the uniform dose within each of SFRT planar beams is made possible by the very large focal spot size (8mm2) of the XRad irradiator (Precision X-ray Inc., North Branford, CT USA). All irradiations in this study used the same irradiator. Animal tumor model This study was carried out in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health (NIH). The University of North Carolina- Chapel Hill Institutional Animal Care and Use Committee (IACUC) reviewed and approved the animal protocol (IACUC ID: 15–366.0) in accordance with NIH standards. All animal surgical, radiation, and imaging procedures were performed under general anesthesia and all efforts were made to minimize suffering. Forty-two eight-week-old female Fischer 344 rats from Charles River Labs and rat fibrosarcoma tumor allografts were used [23]. The rat fibrosarcoma (FSA) allograft model has been well characterized in several radiotherapy response studies by our and collaborator labs [23–25]. Rat FSA is characterized as a local, non-metastasizing tumor that is highly vascular and oxygen dependent. It is an appropriate tumor model for our long-term study goal that investigates the association of SFRT dosimetric parameters with treatment responses, which is reported here, and the association between SFRT treatment response and tumor vascular change post radiation using 3D acoustic angiography. The latter is ongoing research for future publication. All surgical, radiation, and imaging procedures were performed under general anesthesia, induced in the animals initially using 5% vaporized isoflurane mixed with pure oxygen as the carrier gas and then maintained at 2.5% isoflurane mixed with pure oxygen throughout each procedure. Depth of anesthesia was monitored by toe pinch reflex and breathing rate. Opthalmic ointment was placed on the animal’s eyes during anesthesia to provide lubrication and body temperature under anesthesia was maintained via electronically controlled heating pad. Tumors were grown in each rat by implanting freshly resected tumor tissue (1mm3) that was harvested from tumor-bearing donor rats into the subcutaneous space of the rodent flank using blunt dissection. Postoperative care included daily incision surveillance, body temperature monitoring, and a water bottle containing 6mg/mL cherry-flavored, dye-free children’s Tylenol diluted in water for a minimum of 24-hrs post-surgery to alleviate any associated pain from the implantation procedure. Animals were used for experiments 2–3 weeks post-implantation, when the tumors reached the target RT treatment size of approximately 5-10mm. In preclinical studies the pre-treatment tumor volume is known to be strongly correlated with treatment tumor control [23]. We minimize this unwanted effect by controlling the pre-treatment tumor volume in a randomized, matched group study design. We binned animals according to their pre-treatment tumor volume and then randomly assigned these matched bins of animals such that at least one animal from each bin is assigned to each treatment arm. This technique resulted in an average initial tumor volume across arms of 566 +/- 47 mm3 on RT treatment day (pre-treatment tumor volumes for each treatment groups were: 511.84mm3 for the Untreated arm; 602.73mm3 for 20GyUniformRT arm; 608.77mm3 for the 20GyHalfSFRT arm; 489.2mm3 for the 20Gy2mmSFRT arm; 599.81mm3 for the 20GySFRT arm; and 584.06mm3 for the 50GySFRT arm). Biological variability was minimized by ordering animals from the same vendor and of the same age (6 weeks old), implanting tumor on the same day and from the same donor animal, treating with radiation on the same day, and housing animals in the same Vivarium location with identical husbandry conditions. All animals (mixed caged) were provided identical standard laboratory rodent diets of (23%> crude protein) and water ad libitum throughout the study. In addition, all animal diets were supplemented with high-calorie, nutritionally fortified water-based gel cups to help mitigate any potential significant weight loss and dehydration post-radiation. The animals body weight and tumor volumes are monitored prior to radiation and every third day thereafter for up to 30 days. Study endpoints are maximum tumor burden (2.5cm or larger in any dimension), weight loss in excess of 15%, body condition scores [26] less than or equal to 2, or other signs of pain, discomfort, or moribundity as recommended by University of North Carolina- Chapel Hill Division of Comparative Medicine veterinary staff. Animals that met study end-point criteria will be ethically euthanized primarily via compressed carbon dioxide gas or vaporized isoflurane overdose followed by thoracotomy as a secondary means of physical euthanasia per the approved animal study protocol. Animal radiation dosimetry XRad Irradiator and 320kV x-rays were used in this study. Surface dose rates ranging from 4.27 to 5.25 Gy/min were used for all study arms. Fig 2 shows the treatment setup, the radiation light field on animal seen by the camera, and treatment verification films. Dosimetry was measured via EBT-3 film calibrated by an ADCL-calibrated ion chamber under large field conditions. Acrylic phantom measurement setup and beam profile and percentage depth dose (PDD) dosimetry are shown in Figs 3 and 4. The volume-averaged tumor dose was approximated by computing the film average dose within an area of 1 cm x 1 cm (depth) of the PDD film. Tumor and normal tissue EUD were computed using 2D dose distributions measured in a phantom using PDD films under the assumptions that 1) there is no variation in dose distribution along the long dimension of the minibeam array and 2) there is no tissue inhomogeneity. Both assumptions, especially the latter, are not correct. However, we expect that their impact on dosimetry are similar for all study arms and thus have minimum impact on the results of the inter-comparative study among the arms. The widely used EUD calculation has been described by Niemierko et al [27]. EUD dose calculation is a model to compute the impact of organ dose inhomogeneity on radiobiological effect. Parameter ‘a’ used in EUD calculation is estimated per tissue type by Wu et al. for tumors and normal tissues, where a large negative ‘a’ value is characteristic of tumors and a large positive value of ‘a’ is consistent with normal tissues. In this study we use the reported values of a = -10 for tumors and a = 5 for normal tissue such as the bladder or rectum [28]. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 2. Animal irradiation setup and treatment alignment and verification. (A—B) The treatment setup components include (1) X-ray source, (2) endoscopic camera (lens shielded), (3) field shaping collimator for all treated arms (20GySFRT shown), (4) animal and tumor, and the (5) 3-axial heated animal positioning stage. (C) Photo of the built-in irradiator light shines through the 50GySFRT collimator and onto the outlined tumor as seen from the beams-eye view camera (live feed used to position tumor within the treatment fields). (D) EBT-3 treatment verification films with a cutout in the tumor region. The films were reviewed for all treated animals for treatment targeting verification. https://doi.org/10.1371/journal.pone.0229053.g002 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 3. Phantom dosimetry measurement. EBT-3 films were calibrated by ion chamber under large field conditions. All beam profiles and corresponding percentage depth dose were measured using two films as shown: one is on the surface perpendicular to radiation beam (A) and one sandwiched between two small phantom blocks parallel to radiation beam (B). The circles indicate the film areas used for volume-average dose calculation estimates. The following assumption was made for volume-averaged tumor dose and EUD calculations: dose value does not vary +/-1cm along the direction parallel to the same valleys/peaks. https://doi.org/10.1371/journal.pone.0229053.g003 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 4. Measured dose beam profiles and percentage depth doses for all treatment arms. (A-D) Figures display the percentage depth doses for each of the 20 Gy volume-averaged treatment arms. (E-H) Figures display the corresponding SFRT beam profiles for each of the 20 Gy volume-averaged treatment arms. Note that the 20GySFRT and 50GySFRT arms share the same SFRT collimator and thus the same relative dosimetry. The large non-uniformity of the peak doses in the SFRT radiation is due to the finite x-ray target size and the nondivergence of the SFRT collimator. However, the actual peak dose non-uniformity in the treated tumor (diameter of ~10mm) is within 10%. https://doi.org/10.1371/journal.pone.0229053.g004 Animal radiation delivery and verification All of the RT collimators were aligned with x-ray target of the irradiator using film dosimetry. Animals were anesthetized with vaporized isoflurane mixed with oxygen carrier gas and positioned on an electronically controlled heating pad (Fig 2, panels A and B). For radiation tumor targeting we used the light field and a PC-linked camera before radiation and verified it with film dosimetry during each irradiation (Fig 2, panel C). Live video-feed from the camera was used for animal tumor-radiation alignment and for animal monitoring during treatment. Radiation targeting is achieved by (a) delineating the tumor boundary on animal skin using marker pre-treatment, (b) transferring the marking onto the verification EBT-3 film taped on skin and cutting out the tumor portion of the film, (c) taping the film back with the tumor inside the cutout, (d) placing the animal in the irradiator and align the tumor with the radiation, and (e) animal monitoring throughout irradiation. The treatment verification films were reviewed post-radiation for radiation targeting documentation (Fig 2, panel D). Tumor volume imaging and body weight monitoring Three-dimensional B-mode ultrasound imaging of the tumors was performed using a Vevo 770 preclinical ultrasound scanner (Vevo 770, VisualSonics, Toronto, ON, Canada) and the resulting images used to calculate tumor volume, as described in a previous publication [23]. Imaging was performed on the day before treatment as well as every third day post-treatment for approximately 30 days, or when maximal tumor burden was met, at which point the animals were humanely sacrificed per IACUC-approved animal protocol. Fig 5 shows an illustration of the 3D ultrasound tumor imaging setup and acquisition. Three-dimensional imaging is performed by mechanically stepping the ultrasound probe in the elevational dimension and acquired a two-dimensional image at each step (100um step size, 2cm elevational scan length). The reconstructed 3D ultrasound images were used to calculate tumor volume. The longest orthogonal tumor dimensions in each 3D image were measured using the digital caliper feature on the Vevo 770 imaging software and tumor volume was approximated using the volume formula for an ellipsoid, , where V is the calculated tumor volume, and a, b, and c are each the half lengths of the principal axes of the tumor [29]. A sample tumor volume change post radiation from a 20GyHalfSFRT arm animal shows no tumor control (Fig 5, panel D). Animal body weight was measured using the same schedule. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 5. Illustration of 3D ultrasound imaging-based tumor volume measurement. Figure (A) is an illustration of the 3D ultrasound imaging setup with anesthetized animal (23). Two-dimensional transverse image slices (B) are acquired along the elevational direction and are then reconstructed into 3D images [30] (C). Tumors are visually identified on the ultrasound images. Resulting 3D images (C) are used to measure the tumor dimensions and calculate tumor volume. Imaging data is acquired pre-treatment (D) and every ~third day thereafter (E-G). In images D-G the tumor (yellow dotted line) and corresponding tumor volume grow over time following a 20GyHalfSFRT treatment. https://doi.org/10.1371/journal.pone.0229053.g005 Association between SFRT dosimetry and treatment response We analyzed the associations between animal treatment responses and each of the nine dosimetric parameters, listed in Table 1. The treatment responses are time-to-euthanasia, proportion of animals surviving to Day 17, and change in animal body weight on Day 17. In this study we present data on tumor volume change with time; however, we did not use it in determining the association of tumor response with dosimetric parameters for the following reason. When tumors reach the maximum tumor mass defined by the IACUC-approved animal protocol ethical animal euthanasia is performed. As a result, the non-responders with uncontrolled tumor growth are removed from the study and the tumor volume change data only represent the remaining responder animals. The survival data, on the other hand, considers all animals in the study arms, responders and non-responders alike. Day 17 was chosen for the linear regression association studies between treatment response and dosimetry because at this timepoint there is a good compromise between the number of animals available for statistical consideration and the magnitude of radiation effects (20GyUniformRT n = 8, 20GyHalfSFRT n = 3, 20Gy2mmSFRT n = 4, 20GySFRT n = 4, 50GySFRT n = 5, Untreated n = 0). We also fit a more robust Cox Proportional Hazards (CoxPH) model to the full data set that includes all animals. Animal body weight change on Day 17 is used as an indicator of treatment toxicity. Animal body weight change is a gross assessment on treatment toxicity, especially in this study where tumors were implanted in the rodent flank, near the lower gastro-intestinal tract (including the rodent anus, rectum, colon, and cecum) and parts of the upper gastro-intestinal tract (including portions of the small bowel). We speculate that some treatment arms may induce more GI toxicity that others. We subtracted the tumor weight from the measured body weight and regard this “net” animal body weight change as an indication, not confirmation, of treatment toxicity. To confirm any lower GI toxicity, additional tissue histological staining or organ function examination studies would be necessary, both of which are beyond the scope of this work. Statistical methods We computed Product-Limit (Kaplan-Meier) Estimator and Logrank (Mantel-Haenszel) test for statistical significance of survival difference between each pair of treatment arms [31]. Multiple simple linear regression models [32] were used to study the association between dosimetric parameters with animal body weight and percentage survival within treatment arm on Day 17. R2 (square of the Pearson correlation) coefficient is computed to estimate the proportion of variance explained in each of the linear regression models. In general, the greater the magnitude of the test statistic (t or F), the more closely associated the dosimetric parameter studied is with the treatment response (survival or body weight). In addition to linear regressions, we fit Cox Proportional Hazard (CoxPH) models with individual animal survival as the time-to-event outcome, which used data from all dates including Day 17. This allowed us to calculate the hazard ratio associated with the impact of dosimetric parameters on treatment response. We also used a Pearson Correlation matrix to show the cross-correlation between each pair of the dosimetric parameters. All data collected were analyzed using R (version 3.5.3) statistical software available from R Core Team. The associated p-values for all statistics are computed and reported for the significance levels *p < 0.05, **p < 0.01, ***p <0.001. Study design The secret of SFRT lays in its radiation dose spatial fractionation. Although this work does not address the very much needed understanding of working mechanism it addresses another important matter for SFRT application—the association of SFRT dosimetric parameters with treatment response at conventional dose rates (dose rate ranges from 4.27 to 5.25 Gy/min was used). Fig 1 shows a six-arm study design using a very large span of radiation spatial fractionation, constructed to explore the impact of radiation spatial fractionation. Table 1 summarizes the dosimetric parameters of each of the six arms. To study the effect of radiation spatial fractionation under the condition of equal volume-averaged dose we used the following four study arms: 20GyUniformRT (entire tumor directly irradiated), 20GyHalfRT (only one-half of tumor directly irradiated), 20Gy2mmSFRT (50% of tumor directly irradiated by 2mm-wide planar beam array), and 20GySFRT (20% of tumor directly irradiated with 0.3mm-wide planar beam array). Note that the doses are volume-averaged doses computed for the entire tumor volume. A 50GySFRT arm (50 Gy volume-averaged dose, beam width 0.31mm) is added as it has a peak dose of 225 Gy, which is within the known minibeam peak dose range showing tumor control. To account for unavoidable variations in tumor position under the 20Gy2mmSFRT treatment beams during animal irradiations, we computed the maximum and minimum beam coverage positions and calculated their corresponding dosimetric specifications. The 20Gy2mmSFRT treatment arm dosimetric values reported in Table 1 correspond to the average at these positions for a 10mm diameter tumor. For example, a 10mm sized tumor is irradiated by at most three 2mm-peaks and at minimum two 2mm-peaks and the average dosimetric parameter at these two positions was calculated. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 1. Illustration of the SFRT spatial fractionation study design. A very large range of radiation spatial fractionation scale was used to derive the impact of radiation spatial fractionation. Four arms share the same 20 Gy volume-average dose. The high dose 50GySFRT arm is added because 20GySFRT is not known to have tumor control. The dosimetric parameters studied and number of animals per study arm are listed in Table 1. https://doi.org/10.1371/journal.pone.0229053.g001 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 1. Summary of nine SFRT dosimetric parameter specifications in the six-arm study. https://doi.org/10.1371/journal.pone.0229053.t001 Custom-made radiation blocks and collimators made of Cerrobend or tungsten were used to define the 2cmx2cm field for 20GyUniformRT arm treatment, the 2cmx1cm for 20GyHalfSFRT treatment, and the beamlet array 2cmx2cm fields for both the 20Gy2mmSFRT and 20GySFRT/50GySFRT treatments. The 2cm field size in the direction of the uniform dose within each of SFRT planar beams is made possible by the very large focal spot size (8mm2) of the XRad irradiator (Precision X-ray Inc., North Branford, CT USA). All irradiations in this study used the same irradiator. Animal tumor model This study was carried out in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health (NIH). The University of North Carolina- Chapel Hill Institutional Animal Care and Use Committee (IACUC) reviewed and approved the animal protocol (IACUC ID: 15–366.0) in accordance with NIH standards. All animal surgical, radiation, and imaging procedures were performed under general anesthesia and all efforts were made to minimize suffering. Forty-two eight-week-old female Fischer 344 rats from Charles River Labs and rat fibrosarcoma tumor allografts were used [23]. The rat fibrosarcoma (FSA) allograft model has been well characterized in several radiotherapy response studies by our and collaborator labs [23–25]. Rat FSA is characterized as a local, non-metastasizing tumor that is highly vascular and oxygen dependent. It is an appropriate tumor model for our long-term study goal that investigates the association of SFRT dosimetric parameters with treatment responses, which is reported here, and the association between SFRT treatment response and tumor vascular change post radiation using 3D acoustic angiography. The latter is ongoing research for future publication. All surgical, radiation, and imaging procedures were performed under general anesthesia, induced in the animals initially using 5% vaporized isoflurane mixed with pure oxygen as the carrier gas and then maintained at 2.5% isoflurane mixed with pure oxygen throughout each procedure. Depth of anesthesia was monitored by toe pinch reflex and breathing rate. Opthalmic ointment was placed on the animal’s eyes during anesthesia to provide lubrication and body temperature under anesthesia was maintained via electronically controlled heating pad. Tumors were grown in each rat by implanting freshly resected tumor tissue (1mm3) that was harvested from tumor-bearing donor rats into the subcutaneous space of the rodent flank using blunt dissection. Postoperative care included daily incision surveillance, body temperature monitoring, and a water bottle containing 6mg/mL cherry-flavored, dye-free children’s Tylenol diluted in water for a minimum of 24-hrs post-surgery to alleviate any associated pain from the implantation procedure. Animals were used for experiments 2–3 weeks post-implantation, when the tumors reached the target RT treatment size of approximately 5-10mm. In preclinical studies the pre-treatment tumor volume is known to be strongly correlated with treatment tumor control [23]. We minimize this unwanted effect by controlling the pre-treatment tumor volume in a randomized, matched group study design. We binned animals according to their pre-treatment tumor volume and then randomly assigned these matched bins of animals such that at least one animal from each bin is assigned to each treatment arm. This technique resulted in an average initial tumor volume across arms of 566 +/- 47 mm3 on RT treatment day (pre-treatment tumor volumes for each treatment groups were: 511.84mm3 for the Untreated arm; 602.73mm3 for 20GyUniformRT arm; 608.77mm3 for the 20GyHalfSFRT arm; 489.2mm3 for the 20Gy2mmSFRT arm; 599.81mm3 for the 20GySFRT arm; and 584.06mm3 for the 50GySFRT arm). Biological variability was minimized by ordering animals from the same vendor and of the same age (6 weeks old), implanting tumor on the same day and from the same donor animal, treating with radiation on the same day, and housing animals in the same Vivarium location with identical husbandry conditions. All animals (mixed caged) were provided identical standard laboratory rodent diets of (23%> crude protein) and water ad libitum throughout the study. In addition, all animal diets were supplemented with high-calorie, nutritionally fortified water-based gel cups to help mitigate any potential significant weight loss and dehydration post-radiation. The animals body weight and tumor volumes are monitored prior to radiation and every third day thereafter for up to 30 days. Study endpoints are maximum tumor burden (2.5cm or larger in any dimension), weight loss in excess of 15%, body condition scores [26] less than or equal to 2, or other signs of pain, discomfort, or moribundity as recommended by University of North Carolina- Chapel Hill Division of Comparative Medicine veterinary staff. Animals that met study end-point criteria will be ethically euthanized primarily via compressed carbon dioxide gas or vaporized isoflurane overdose followed by thoracotomy as a secondary means of physical euthanasia per the approved animal study protocol. Animal radiation dosimetry XRad Irradiator and 320kV x-rays were used in this study. Surface dose rates ranging from 4.27 to 5.25 Gy/min were used for all study arms. Fig 2 shows the treatment setup, the radiation light field on animal seen by the camera, and treatment verification films. Dosimetry was measured via EBT-3 film calibrated by an ADCL-calibrated ion chamber under large field conditions. Acrylic phantom measurement setup and beam profile and percentage depth dose (PDD) dosimetry are shown in Figs 3 and 4. The volume-averaged tumor dose was approximated by computing the film average dose within an area of 1 cm x 1 cm (depth) of the PDD film. Tumor and normal tissue EUD were computed using 2D dose distributions measured in a phantom using PDD films under the assumptions that 1) there is no variation in dose distribution along the long dimension of the minibeam array and 2) there is no tissue inhomogeneity. Both assumptions, especially the latter, are not correct. However, we expect that their impact on dosimetry are similar for all study arms and thus have minimum impact on the results of the inter-comparative study among the arms. The widely used EUD calculation has been described by Niemierko et al [27]. EUD dose calculation is a model to compute the impact of organ dose inhomogeneity on radiobiological effect. Parameter ‘a’ used in EUD calculation is estimated per tissue type by Wu et al. for tumors and normal tissues, where a large negative ‘a’ value is characteristic of tumors and a large positive value of ‘a’ is consistent with normal tissues. In this study we use the reported values of a = -10 for tumors and a = 5 for normal tissue such as the bladder or rectum [28]. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 2. Animal irradiation setup and treatment alignment and verification. (A—B) The treatment setup components include (1) X-ray source, (2) endoscopic camera (lens shielded), (3) field shaping collimator for all treated arms (20GySFRT shown), (4) animal and tumor, and the (5) 3-axial heated animal positioning stage. (C) Photo of the built-in irradiator light shines through the 50GySFRT collimator and onto the outlined tumor as seen from the beams-eye view camera (live feed used to position tumor within the treatment fields). (D) EBT-3 treatment verification films with a cutout in the tumor region. The films were reviewed for all treated animals for treatment targeting verification. https://doi.org/10.1371/journal.pone.0229053.g002 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 3. Phantom dosimetry measurement. EBT-3 films were calibrated by ion chamber under large field conditions. All beam profiles and corresponding percentage depth dose were measured using two films as shown: one is on the surface perpendicular to radiation beam (A) and one sandwiched between two small phantom blocks parallel to radiation beam (B). The circles indicate the film areas used for volume-average dose calculation estimates. The following assumption was made for volume-averaged tumor dose and EUD calculations: dose value does not vary +/-1cm along the direction parallel to the same valleys/peaks. https://doi.org/10.1371/journal.pone.0229053.g003 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 4. Measured dose beam profiles and percentage depth doses for all treatment arms. (A-D) Figures display the percentage depth doses for each of the 20 Gy volume-averaged treatment arms. (E-H) Figures display the corresponding SFRT beam profiles for each of the 20 Gy volume-averaged treatment arms. Note that the 20GySFRT and 50GySFRT arms share the same SFRT collimator and thus the same relative dosimetry. The large non-uniformity of the peak doses in the SFRT radiation is due to the finite x-ray target size and the nondivergence of the SFRT collimator. However, the actual peak dose non-uniformity in the treated tumor (diameter of ~10mm) is within 10%. https://doi.org/10.1371/journal.pone.0229053.g004 Animal radiation delivery and verification All of the RT collimators were aligned with x-ray target of the irradiator using film dosimetry. Animals were anesthetized with vaporized isoflurane mixed with oxygen carrier gas and positioned on an electronically controlled heating pad (Fig 2, panels A and B). For radiation tumor targeting we used the light field and a PC-linked camera before radiation and verified it with film dosimetry during each irradiation (Fig 2, panel C). Live video-feed from the camera was used for animal tumor-radiation alignment and for animal monitoring during treatment. Radiation targeting is achieved by (a) delineating the tumor boundary on animal skin using marker pre-treatment, (b) transferring the marking onto the verification EBT-3 film taped on skin and cutting out the tumor portion of the film, (c) taping the film back with the tumor inside the cutout, (d) placing the animal in the irradiator and align the tumor with the radiation, and (e) animal monitoring throughout irradiation. The treatment verification films were reviewed post-radiation for radiation targeting documentation (Fig 2, panel D). Tumor volume imaging and body weight monitoring Three-dimensional B-mode ultrasound imaging of the tumors was performed using a Vevo 770 preclinical ultrasound scanner (Vevo 770, VisualSonics, Toronto, ON, Canada) and the resulting images used to calculate tumor volume, as described in a previous publication [23]. Imaging was performed on the day before treatment as well as every third day post-treatment for approximately 30 days, or when maximal tumor burden was met, at which point the animals were humanely sacrificed per IACUC-approved animal protocol. Fig 5 shows an illustration of the 3D ultrasound tumor imaging setup and acquisition. Three-dimensional imaging is performed by mechanically stepping the ultrasound probe in the elevational dimension and acquired a two-dimensional image at each step (100um step size, 2cm elevational scan length). The reconstructed 3D ultrasound images were used to calculate tumor volume. The longest orthogonal tumor dimensions in each 3D image were measured using the digital caliper feature on the Vevo 770 imaging software and tumor volume was approximated using the volume formula for an ellipsoid, , where V is the calculated tumor volume, and a, b, and c are each the half lengths of the principal axes of the tumor [29]. A sample tumor volume change post radiation from a 20GyHalfSFRT arm animal shows no tumor control (Fig 5, panel D). Animal body weight was measured using the same schedule. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 5. Illustration of 3D ultrasound imaging-based tumor volume measurement. Figure (A) is an illustration of the 3D ultrasound imaging setup with anesthetized animal (23). Two-dimensional transverse image slices (B) are acquired along the elevational direction and are then reconstructed into 3D images [30] (C). Tumors are visually identified on the ultrasound images. Resulting 3D images (C) are used to measure the tumor dimensions and calculate tumor volume. Imaging data is acquired pre-treatment (D) and every ~third day thereafter (E-G). In images D-G the tumor (yellow dotted line) and corresponding tumor volume grow over time following a 20GyHalfSFRT treatment. https://doi.org/10.1371/journal.pone.0229053.g005 Association between SFRT dosimetry and treatment response We analyzed the associations between animal treatment responses and each of the nine dosimetric parameters, listed in Table 1. The treatment responses are time-to-euthanasia, proportion of animals surviving to Day 17, and change in animal body weight on Day 17. In this study we present data on tumor volume change with time; however, we did not use it in determining the association of tumor response with dosimetric parameters for the following reason. When tumors reach the maximum tumor mass defined by the IACUC-approved animal protocol ethical animal euthanasia is performed. As a result, the non-responders with uncontrolled tumor growth are removed from the study and the tumor volume change data only represent the remaining responder animals. The survival data, on the other hand, considers all animals in the study arms, responders and non-responders alike. Day 17 was chosen for the linear regression association studies between treatment response and dosimetry because at this timepoint there is a good compromise between the number of animals available for statistical consideration and the magnitude of radiation effects (20GyUniformRT n = 8, 20GyHalfSFRT n = 3, 20Gy2mmSFRT n = 4, 20GySFRT n = 4, 50GySFRT n = 5, Untreated n = 0). We also fit a more robust Cox Proportional Hazards (CoxPH) model to the full data set that includes all animals. Animal body weight change on Day 17 is used as an indicator of treatment toxicity. Animal body weight change is a gross assessment on treatment toxicity, especially in this study where tumors were implanted in the rodent flank, near the lower gastro-intestinal tract (including the rodent anus, rectum, colon, and cecum) and parts of the upper gastro-intestinal tract (including portions of the small bowel). We speculate that some treatment arms may induce more GI toxicity that others. We subtracted the tumor weight from the measured body weight and regard this “net” animal body weight change as an indication, not confirmation, of treatment toxicity. To confirm any lower GI toxicity, additional tissue histological staining or organ function examination studies would be necessary, both of which are beyond the scope of this work. Statistical methods We computed Product-Limit (Kaplan-Meier) Estimator and Logrank (Mantel-Haenszel) test for statistical significance of survival difference between each pair of treatment arms [31]. Multiple simple linear regression models [32] were used to study the association between dosimetric parameters with animal body weight and percentage survival within treatment arm on Day 17. R2 (square of the Pearson correlation) coefficient is computed to estimate the proportion of variance explained in each of the linear regression models. In general, the greater the magnitude of the test statistic (t or F), the more closely associated the dosimetric parameter studied is with the treatment response (survival or body weight). In addition to linear regressions, we fit Cox Proportional Hazard (CoxPH) models with individual animal survival as the time-to-event outcome, which used data from all dates including Day 17. This allowed us to calculate the hazard ratio associated with the impact of dosimetric parameters on treatment response. We also used a Pearson Correlation matrix to show the cross-correlation between each pair of the dosimetric parameters. All data collected were analyzed using R (version 3.5.3) statistical software available from R Core Team. The associated p-values for all statistics are computed and reported for the significance levels *p < 0.05, **p < 0.01, ***p <0.001. Results Overall treatment response Fig 6 shows (A) percentage animal survival, (B) normalized averaged tumor volume, and (C) normalized averaged “net” body weight post treatment for all 6 study arms. In addition, the differences between survival curves pairs are shown in S2 Table. No animal died of body condition deterioration in this study. All endpoints were due to ethical animal euthanasia triggered by tumors exceeding the maximum allowable burden per IACUC-approved animal protocol limitations. Our data shows that the 20GyUnformRT arm has the best tumor control followed by the 50GySFRT and 20Gy2mmSFRT arms. Note that among the four arms sharing similar volume-averaged dose (20 Gy or 18 Gy) survival varies greatly, from 33% to 100% on Day 17, which is a strong indication that volume-averaged dose is poorly associated with tumor treatment response. The tumor volume data indicate that although 50GySFRT arm and 20Gy2mmSFRT arm have similar survival the former has a better tumor volume reduction than the latter arm. Only the 20GyUniformRT arm experienced weight loss post-treatment and then recovered back to pre-treatment weight after week three. The 20GySFRT and 20Gy2mmSFRT arms experienced similar body weight gains as the untreated arm, indicating little treatment toxicity from the two SFRT treatments. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 6. Animal survival, tumor volume change, and body weight change post-treatment. Animal survival (A), normalized averaged tumor volume (B), and normalized averaged body weight (C) are shown for all six study arms. In plots (B) and (C), the error bars at each time point represent one standard deviation from the mean for each treatment arm. The statistical significance values for each pair of survival curves are shown in S2 Table. https://doi.org/10.1371/journal.pone.0229053.g006 Cross-correlation in the SFRT dosimetry parameters In this study not all dosimetric parameters studied are independent variables and their cross-correlations are shown in the table of Pearson Correlation coefficients (Table 2). The larger the magnitude of the coefficient, the more co-linear and correlated the pair of dosimetric parameters. In this study, peak width and valley width are perfectly co-linear (correlation of 1.0) by study design. Valley (minimum) dose, a parameter used in tumor EUD calculation, is also highly correlated with tumor EUD (correlation of 0.99). These strong correlations explain the similar statistical associations of these parameters with treatment responses. These correlations also limited the study’s ability to better exam the association between a given treatment response with each of the dosimetric parameters. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 2. Pearson Correlation coefficient matrix for the eight SFRT dosimetric parameters relevant for tumor treatment response. https://doi.org/10.1371/journal.pone.0229053.t002 Association between tumor response and SFRT dosimetry We associated eight dosimetric parameters with percentage of animals surviving to Day 17 and with the survival curves shown in Fig 6. Fig 7 shows scatter plots of eight tumor-related dosimetric parameters vs. percentage survival at Day 17, each fitted with a corresponding regression line, R2 (Fig 7). Tumor EUD (R2 = 0.7923, F-stat = 15.26*), Valley dose (R2 = 0.7636, F-stat = 12.92*), and percentage volume directly irradiated (R2 = 0.7153, F-stat = 10.05*) are the top three most statistically significant dosimetric parameters in terms of association with the animal survival at Day 17 (see S1 Table). Peak dose (R2 = 0.04472, F-stat = 0.6874 (not sig.)) and AVG Dose (R2 = 0.2745, F-stat = 1.514 (not sig.)) showed little association with survival. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 7. Associations between Percentage Survival (Day 17) and eight dosimetric parameters. The percentage survival at Day 17 for each treatment arm is plotted against their corresponding (A) Tumor EUD, (B) valley dose, (C) percentage volume irradiated, (D) valley width, (E) peak width, (F) volume-averaged dose, (G) peak dose, (H) and PVDR. In addition, their corresponding regression lines and R2 values are presented. Eight linear regression models with single covariates, one for each dosimetric parameter, were used to calculate the R2 value and corresponding statistics. https://doi.org/10.1371/journal.pone.0229053.g007 To validate the above finding in Fig 7 we used data from the entire survival curves in Univariate Cox Proportional Hazards analysis and the results are shown in Table 3. The results from the Univariate Cox Proportional Hazards analysis confirms the results from the linear regression analysis—among the eight dosimetric parameters analyzed tumor EUD (z-stat = -4.07***), valley (minimum) dose (z-stat = -4.338***), and percentage tumor volume directly irradiated (z-stat = -3.837***) have the closest associations with animal survival. Compared to the linear regression analysis (Fig 7) the improved p-values in the CoxPH model analysis is likely due to the increased sample size. The Hazard Ratio shows the impact of change in each of the dosimetric parameters to the hazard rate (risk of death). For instance, when valley (minimum) dose parameter changes by 1 Gy, the hazard rate (risk of death) changes by 19% (95% CI, 26% - 11%) with p-value of 1.44x10-5. For a 1 Gy change in peak dose, the corresponding change in hazard rate is 0.2% (95% CI, 0.7% - 0.3%) with p-value of 0.432. Three additional statistical tests were used to validate the CoxPH z-test statistics results for each model (Likelihood Ratio Test, Wald Test, and Logrank Test) and all three tests largely agree with the results presented in Table 3. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 3. Table of coefficients for univariate Cox Proportional Hazards analysis of survival. https://doi.org/10.1371/journal.pone.0229053.t003 Association between body weight change and SFRT dosimetry Eight dosimetric parameters are associated with the body weight change on Day-17. Note that the body weight is the measured body weight subtracted the measured tumor weight to remove the influence of tumor size on the analysis. Fig 8 is a scatter plot of the dosimetric parameters vs. the body weight at Day 17 for each individual animal in the study, as well as the average body weight for each of the treatment arms. This time point was chosen for both the tumor and body weight study because it is a good compromise between data statistics and magnitude of treatment response. Table 4 is a table of coefficients for the corresponding linear regression models used in Fig 8 In general, the greater the magnitude of the t statistic, the greater the individual parameter association with Body Weight (Day 17). For the F-statistic, the greater the statistic value, the more closely associated the model is with Body Weight (Day 17). Based on the t statistics and F-statistics, among the eight dosimetric parameters studied the Valley Dose has the greatest, yet modest, association with Body Weight (Day 17). The order of decreasing association with the body weight change are: valley dose (R2 = 0.3814, F-stat = 13.45**), valley width (R2 = 0.2853, F-stat = 8.783*), peak width (R2 = 0.2759, F-stat = 8.382*), percentage volume irradiated (R2 = 0.1985, F-stat = 5.448*), PVDR (R2 = 0.1203, F-stat = 3.009 (not sig.)), volume-averaged dose (R2 = 0.03308, F-stat = 0.7526 (not sig.)), normal tissue EUD (R2 = 1.022x10-03, F-stat = 0.882 (not sig.)), and peak dose (R2 = 5.99x10-06, F-stat = 1.32x10-04 (not sig.)). A strong similarity between the peak width and valley width association with body weight is expected (see discussion on Cross-Correlations, Table 2). No significant association is observed between body weight change post radiation and PVDR, average dose, normal tissue EUD, and peak dose. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 8. Associations between body weight change (Day 17) and eight dosimetric parameters. Individual animal body weight change are plotted against each of the 8 treatment dosimetric parameters: valley dose (A), valley width (B), peak width (C), percentage volume irradiated (D), normal tissue EUD (E), PVDR (F), volume-averaged dose (G), and peak dose (H) vs % Body Weight at Day 17 and their corresponding regression lines and R2 values are shown. Eight linear regression models with single covariates, one for each dosimetric parameter, were used to calculate the R2 value and corresponding statistics. The body weight change averages for each of the treatment arms are also shown (black bars). https://doi.org/10.1371/journal.pone.0229053.g008 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 4. Table of coefficients for univariate linear regression analysis of Body Weight (Day 17). https://doi.org/10.1371/journal.pone.0229053.t004 Overall treatment response Fig 6 shows (A) percentage animal survival, (B) normalized averaged tumor volume, and (C) normalized averaged “net” body weight post treatment for all 6 study arms. In addition, the differences between survival curves pairs are shown in S2 Table. No animal died of body condition deterioration in this study. All endpoints were due to ethical animal euthanasia triggered by tumors exceeding the maximum allowable burden per IACUC-approved animal protocol limitations. Our data shows that the 20GyUnformRT arm has the best tumor control followed by the 50GySFRT and 20Gy2mmSFRT arms. Note that among the four arms sharing similar volume-averaged dose (20 Gy or 18 Gy) survival varies greatly, from 33% to 100% on Day 17, which is a strong indication that volume-averaged dose is poorly associated with tumor treatment response. The tumor volume data indicate that although 50GySFRT arm and 20Gy2mmSFRT arm have similar survival the former has a better tumor volume reduction than the latter arm. Only the 20GyUniformRT arm experienced weight loss post-treatment and then recovered back to pre-treatment weight after week three. The 20GySFRT and 20Gy2mmSFRT arms experienced similar body weight gains as the untreated arm, indicating little treatment toxicity from the two SFRT treatments. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 6. Animal survival, tumor volume change, and body weight change post-treatment. Animal survival (A), normalized averaged tumor volume (B), and normalized averaged body weight (C) are shown for all six study arms. In plots (B) and (C), the error bars at each time point represent one standard deviation from the mean for each treatment arm. The statistical significance values for each pair of survival curves are shown in S2 Table. https://doi.org/10.1371/journal.pone.0229053.g006 Cross-correlation in the SFRT dosimetry parameters In this study not all dosimetric parameters studied are independent variables and their cross-correlations are shown in the table of Pearson Correlation coefficients (Table 2). The larger the magnitude of the coefficient, the more co-linear and correlated the pair of dosimetric parameters. In this study, peak width and valley width are perfectly co-linear (correlation of 1.0) by study design. Valley (minimum) dose, a parameter used in tumor EUD calculation, is also highly correlated with tumor EUD (correlation of 0.99). These strong correlations explain the similar statistical associations of these parameters with treatment responses. These correlations also limited the study’s ability to better exam the association between a given treatment response with each of the dosimetric parameters. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 2. Pearson Correlation coefficient matrix for the eight SFRT dosimetric parameters relevant for tumor treatment response. https://doi.org/10.1371/journal.pone.0229053.t002 Association between tumor response and SFRT dosimetry We associated eight dosimetric parameters with percentage of animals surviving to Day 17 and with the survival curves shown in Fig 6. Fig 7 shows scatter plots of eight tumor-related dosimetric parameters vs. percentage survival at Day 17, each fitted with a corresponding regression line, R2 (Fig 7). Tumor EUD (R2 = 0.7923, F-stat = 15.26*), Valley dose (R2 = 0.7636, F-stat = 12.92*), and percentage volume directly irradiated (R2 = 0.7153, F-stat = 10.05*) are the top three most statistically significant dosimetric parameters in terms of association with the animal survival at Day 17 (see S1 Table). Peak dose (R2 = 0.04472, F-stat = 0.6874 (not sig.)) and AVG Dose (R2 = 0.2745, F-stat = 1.514 (not sig.)) showed little association with survival. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 7. Associations between Percentage Survival (Day 17) and eight dosimetric parameters. The percentage survival at Day 17 for each treatment arm is plotted against their corresponding (A) Tumor EUD, (B) valley dose, (C) percentage volume irradiated, (D) valley width, (E) peak width, (F) volume-averaged dose, (G) peak dose, (H) and PVDR. In addition, their corresponding regression lines and R2 values are presented. Eight linear regression models with single covariates, one for each dosimetric parameter, were used to calculate the R2 value and corresponding statistics. https://doi.org/10.1371/journal.pone.0229053.g007 To validate the above finding in Fig 7 we used data from the entire survival curves in Univariate Cox Proportional Hazards analysis and the results are shown in Table 3. The results from the Univariate Cox Proportional Hazards analysis confirms the results from the linear regression analysis—among the eight dosimetric parameters analyzed tumor EUD (z-stat = -4.07***), valley (minimum) dose (z-stat = -4.338***), and percentage tumor volume directly irradiated (z-stat = -3.837***) have the closest associations with animal survival. Compared to the linear regression analysis (Fig 7) the improved p-values in the CoxPH model analysis is likely due to the increased sample size. The Hazard Ratio shows the impact of change in each of the dosimetric parameters to the hazard rate (risk of death). For instance, when valley (minimum) dose parameter changes by 1 Gy, the hazard rate (risk of death) changes by 19% (95% CI, 26% - 11%) with p-value of 1.44x10-5. For a 1 Gy change in peak dose, the corresponding change in hazard rate is 0.2% (95% CI, 0.7% - 0.3%) with p-value of 0.432. Three additional statistical tests were used to validate the CoxPH z-test statistics results for each model (Likelihood Ratio Test, Wald Test, and Logrank Test) and all three tests largely agree with the results presented in Table 3. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 3. Table of coefficients for univariate Cox Proportional Hazards analysis of survival. https://doi.org/10.1371/journal.pone.0229053.t003 Association between body weight change and SFRT dosimetry Eight dosimetric parameters are associated with the body weight change on Day-17. Note that the body weight is the measured body weight subtracted the measured tumor weight to remove the influence of tumor size on the analysis. Fig 8 is a scatter plot of the dosimetric parameters vs. the body weight at Day 17 for each individual animal in the study, as well as the average body weight for each of the treatment arms. This time point was chosen for both the tumor and body weight study because it is a good compromise between data statistics and magnitude of treatment response. Table 4 is a table of coefficients for the corresponding linear regression models used in Fig 8 In general, the greater the magnitude of the t statistic, the greater the individual parameter association with Body Weight (Day 17). For the F-statistic, the greater the statistic value, the more closely associated the model is with Body Weight (Day 17). Based on the t statistics and F-statistics, among the eight dosimetric parameters studied the Valley Dose has the greatest, yet modest, association with Body Weight (Day 17). The order of decreasing association with the body weight change are: valley dose (R2 = 0.3814, F-stat = 13.45**), valley width (R2 = 0.2853, F-stat = 8.783*), peak width (R2 = 0.2759, F-stat = 8.382*), percentage volume irradiated (R2 = 0.1985, F-stat = 5.448*), PVDR (R2 = 0.1203, F-stat = 3.009 (not sig.)), volume-averaged dose (R2 = 0.03308, F-stat = 0.7526 (not sig.)), normal tissue EUD (R2 = 1.022x10-03, F-stat = 0.882 (not sig.)), and peak dose (R2 = 5.99x10-06, F-stat = 1.32x10-04 (not sig.)). A strong similarity between the peak width and valley width association with body weight is expected (see discussion on Cross-Correlations, Table 2). No significant association is observed between body weight change post radiation and PVDR, average dose, normal tissue EUD, and peak dose. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 8. Associations between body weight change (Day 17) and eight dosimetric parameters. Individual animal body weight change are plotted against each of the 8 treatment dosimetric parameters: valley dose (A), valley width (B), peak width (C), percentage volume irradiated (D), normal tissue EUD (E), PVDR (F), volume-averaged dose (G), and peak dose (H) vs % Body Weight at Day 17 and their corresponding regression lines and R2 values are shown. Eight linear regression models with single covariates, one for each dosimetric parameter, were used to calculate the R2 value and corresponding statistics. The body weight change averages for each of the treatment arms are also shown (black bars). https://doi.org/10.1371/journal.pone.0229053.g008 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 4. Table of coefficients for univariate linear regression analysis of Body Weight (Day 17). https://doi.org/10.1371/journal.pone.0229053.t004 Discussion Study limitations There are several limitations in this study, many of which are discussed below. (i) There was no image-guidance used in the irradiation study. Our remedy for the lack of online imaging technology included the use of light field and video-based animal alignment, of treatment verification film, and lastly, removal of 2 treatment-misaligned animals from the study. This was judged from reviewing the treatment verification film for each animal. Our remedy worked well resulting in a 20GyHalfRT arm % volume irradiated of 47.8% (±2.2) near the target value of 50% (S1 Fig) and the 2 animals were removed in cases where the magnitude of the % Area Treated differed by more than 5%. (ii) No CT-based treatment planning. Based on the anatomical location of the implanted tumor (rodent flank) we believe a portion of the rodent GI tract may have been irradiated but the actual volume irradiated is unknown. Because all animals were randomized across study arms such that all arms have the same average pre-treatment tumor size and similar tumor location distribution, it is reasonable to assume that any variations in portions of GI track irradiated do not bias any particular study arm. (iii) Only a single tumor model used. The FSA rat tumor model does not represent tumors with low vascularity, which may have different treatment responses. The study should be repeated using different tumor and animal models. (iv) The dosimetric parameters have strong cross correlations in this study, which is discussed in the Cross-Correlations section and in Table 2. The potential impact of spatial fractionation pattern (lines vs. dots, for instance) on treatment response is beyond the scope of this work. However, it is a very important question that deserves methodical investigations as some spatial fractionation patterns are easier to achieve than others in practical application. Our data shows that valley (minimum) dose has the closest association with treatment response for tumor and body weight. However, different spatial fractionation patterns with the same valley dose may not lead to the same treatment response when a different endpoint is used. In our study the 20Gy2mmSFRT arm and the 20GySFRT arm have similar valley doses but dissimilar survival fraction on Day 17. To investigate the impact of radiation spatial fractionation pattern alone on given treatment responses, carefully designed new studies are needed. The exciting noncytotoxic effects of SFRT, such as induction therapy to sensitize tumor to increase therapeutic ratio of the following therapy including anti-tumor immunotherapy, remain largely underexplored [33]; however, they are also beyond of the scope of this work. Our own and others’ work have demonstrated that SFRT radiation impacts tumor microenvironment and modulates immune system very differently than uniform radiation therapy [34–36]. We intend to conduct similar studies to identify associations between dosimetric parameters and these indirect effects of SFRT in the future. SFRT dosimetric association with treatment tumor response Valley dose and tumor EUD. The importance of tumor minimum dose to tumor control has long been established in conventional radiation therapy [37]. Does the same association between tumor control and minimum/valley dose hold for SFRT? For some the answer is yes and sophisticated techniques have been developed to “fill up” the dose valleys in an MRT beam by interlacing the microbeams from MRT from different irradiation angles. As a result, a uniform dose distribution inside the tumor is reached [38] while the surrounding normal tissue out of the “cross-firing” range still receive largely MRT radiation pattern of peaks and valleys. In a synchrotron-MRT study Ibahim et al. [39] reported that valley dose is closely correlated with cell survival, but valley dose alone does not determine the observed radiobiological effects. Our study shows that the tumor EUD (a = -10) and valley (minimum) tumor dose have the highest linear associations (R2 = 0.7923, F-stat = 15.26*; R2 = 0.7636, F-stat = 12.92*, respectively) with tumor treatment response (Fig 7 and S1 Table). This observed association between tumor treatment response with tumor valley minimum) dose and tumor EUD dose in this preclinical study is consistent with their known association in tumor treatment response seen in clinical conventional uniform dose radiation therapy. Our data suggest that valley (minimum) dose or Tumor EUD are more appropriate than peak dose for SFRT treatment prescription. When tumor control is the endpoint, we suggest that equal valley or minimum dose be used for comparative study between a uniform radiation and SFRT therapy or among different SFRT treatments. PVDR. Our data showed that PVDR has a consistent but not statistically significant association with tumor treatment response (R2 = 0.7194, F-stat = 7.691) (Fig 7, S1 Table). The linear regression analysis on day 17 was not statistically significant. The CoxPH analysis using the entire survival data set show a modest association with survival. Although not statistically significant, an inverse association is observed between PVDR value and survival fraction on Day 17—the higher PVDR value the less survival fraction. The inverse association is largely determined by the uniform radiation arm where PVDR value is 1.0. If this data point is removed, the PVDR association with survival for all SFRT arms is inconclusive (Fig 7). We believe this result of inverse association is likely biased by the study design that has very limited PVDR values (4 values) and strong cross-correlations between PVDR and other SFRT parameters (see Cross-Correlations discussion and Table 2, above). A better understanding of PVDR’s association with a given treatment response requires a carefully designed new study that focuses on the impact of PVDR value on treatment response. Percentage volume irradiated, peak width, and valley width. It seems logical that tumor treatment response is closely associated with the tumor volume irradiated. However, this is not supported by a clinical GRID therapy study by Neuner et al. [2] where both MLC-based and collimator-based GRID treatments showed similar response rates for pain, mass effect, other patient complaints, and have similar adverse reactions. The collimator-generated GRID had 50% of the radiation field open while the MLC-generated GRID had only 31% open. In our study the 20GyHalfSFRT and 20Gy2mmSFRT arms have similar percentage-volume-irradiated (as well as PDD curves) but there is a difference of 5 days in the 50% survival time (Fig 1 and Fig 6). Nonetheless, our data shows that percentage-volume-irradiated has the 3rd highest linear association (R2 = 0.7153, F-stat = 10.05*) with tumor treatment response (Fig 7 and Table 3). Since percentage-volume-irradiated is jointly determined by peak width and valley width it is understandable to see moderate associations between tumor treatment response and peak width (R2 = 0.4201, F-stat = 2.898 (not sig.)) and valley width (R2 = 0.4296, F-stat = 3.012 (not sig.)). In a synchrotron microbeam brain study using multiple beams Serduc et al. kept valley dose constant while varying peak width and peak dose. They concluded that the latter two parameters have strong influence therapeutic ratio [40]. Volume-averaged dose and peak dose. This study is designed to scrutinize the association of volume-averaged dose with tumor treatment response (Fig 1). The four study arms sharing very similar volume-averaged doses (20 or 18 Gy) exhibited very different tumor treatment responses (Figs 6 and 7) showing the survival rate at day 17 varied from 100% to 33%. Therefore, the association between volume-average dose and tumor treatment response is weak. We found that peak dose has little to no association with tumor treatment response (R2 = 0.04472, F-stat = 0.6874 (not sig.)) (Fig 7, S1 Table, Table 3). This finding is significant because peak dose has been used for treatment prescription in practically all SFRT treatments (8) (9). Although the linear regression analysis on day 17 showed a weak association between peak dose and survival that was not statistically significant, the CoxPH analysis using the entire survival data set did show a modest association with survival. SFRT dosimetric association with normal tissue toxicity We did not study treatment induced normal tissue toxicity directly in this study. We used body weight change post radiation (targeted to the flank, lower abdominal region of the animal) as an indicator, not an evidence of normal tissue toxicity. We did not see a strong association between animal body weight change and any of the eight dosimetric parameters studied, except a modest association with valley (minimum) dose. Valley dose. The strongest association we observed is a weak one between body weight change and valley (minimum) dose (R2 = 0.3814, F-stat = 13.56**) (Table 4). Note that valley (minimum) dose is also strongly associated with tumor treatment response (R2 = 0.7636, F-stat = 12.92*). Our finding is consistent with a normal mouse brain MRT study Nakayma et al. reported that valley dose is one of the important factors to determine normal brain dose tolerance [41]. Our data suggests that among the nine dosimetric parameters we have studied, valley dose may have a close correlation with both tumor control and toxicity, and thus is a crucial dosimetric parameter in SFRT treatment. Valley width, peak width, percentage volume irradiated. The valley width, peak width, and percentage volume of the tumor that is irradiated parameters were only weakly associated with animal body weight change post radiation (R2 = 0.2853, F-stat = 8.783**; R2 = 0.2759, F-stat = 8.382**; and R2 = 0.1985, F-stat = 5.448*, respectively) (Fig 8 and Table 4). Note that in this study peak width and valley width are closely correlated (see discussion on cross-correlations, Table 2). Percentage volume directly irradiated showed no statistically important association with body weight change. Neuner et al. reported that they observed similar treatment responses from clinical GRID treatments of different percentages of volume directly irradiated [2]. Normal tissue EUD, PVDR, volume-averaged dose, peak dose. The normal tissue EUD, PVDR, volume-averaged dose, and peak dose parameters showed little to no association with body weight change post radiation (R2 = 1.022x10-03, F-stat = 0.882 (not sig.); R2 = 0.1203, F-stat = 3.009 (not sig.); R2 = 0.03308, F-stat = .7526 (not sig.); and R2 = 5.99x10-06, F-stat = 1.32x10-04 (not sig.), respectively). Our finding is consistent with a rat normal brain minibeam study by Prezado et al. showing arms with similar volume-average-doses have drastic differences in survival (14) and inconsistent with a MRT study on normal mouse skin by Priyadarshika et al. concluded that integrated dose (i.e., volume-averaged dose) rather than peak or valley dose, may dictate the acute skin toxicity [42]. 2mm wide beam array SFRT. Our data indicates that the 20Gy2mmSFRT arm is not only the most relevant to clinical application because of its millimeter scale, but it also has the potential for superior therapeutic ratio. The 20Gy2mmSFRT arm showed similar survival with the 50GySFRT arm but has significantly lower valley dose (6.2 Gy vs. 17 Gy). At the same time, it showed the least, if any, body weight change compared to the untreated arm while the 50GySFRT arm with 0.31mm beam width exhibited significant body weight growth deficit (Table 1 and Fig 6). The 20GyUniform arm has the best tumor treatment response, but also had a large drop in body weight. Our data indicated that for KV x-rays the 2mm wide beam array has the potential to be the SFRT pattern for high therapeutic ratio and deserves further investigation. Study limitations There are several limitations in this study, many of which are discussed below. (i) There was no image-guidance used in the irradiation study. Our remedy for the lack of online imaging technology included the use of light field and video-based animal alignment, of treatment verification film, and lastly, removal of 2 treatment-misaligned animals from the study. This was judged from reviewing the treatment verification film for each animal. Our remedy worked well resulting in a 20GyHalfRT arm % volume irradiated of 47.8% (±2.2) near the target value of 50% (S1 Fig) and the 2 animals were removed in cases where the magnitude of the % Area Treated differed by more than 5%. (ii) No CT-based treatment planning. Based on the anatomical location of the implanted tumor (rodent flank) we believe a portion of the rodent GI tract may have been irradiated but the actual volume irradiated is unknown. Because all animals were randomized across study arms such that all arms have the same average pre-treatment tumor size and similar tumor location distribution, it is reasonable to assume that any variations in portions of GI track irradiated do not bias any particular study arm. (iii) Only a single tumor model used. The FSA rat tumor model does not represent tumors with low vascularity, which may have different treatment responses. The study should be repeated using different tumor and animal models. (iv) The dosimetric parameters have strong cross correlations in this study, which is discussed in the Cross-Correlations section and in Table 2. The potential impact of spatial fractionation pattern (lines vs. dots, for instance) on treatment response is beyond the scope of this work. However, it is a very important question that deserves methodical investigations as some spatial fractionation patterns are easier to achieve than others in practical application. Our data shows that valley (minimum) dose has the closest association with treatment response for tumor and body weight. However, different spatial fractionation patterns with the same valley dose may not lead to the same treatment response when a different endpoint is used. In our study the 20Gy2mmSFRT arm and the 20GySFRT arm have similar valley doses but dissimilar survival fraction on Day 17. To investigate the impact of radiation spatial fractionation pattern alone on given treatment responses, carefully designed new studies are needed. The exciting noncytotoxic effects of SFRT, such as induction therapy to sensitize tumor to increase therapeutic ratio of the following therapy including anti-tumor immunotherapy, remain largely underexplored [33]; however, they are also beyond of the scope of this work. Our own and others’ work have demonstrated that SFRT radiation impacts tumor microenvironment and modulates immune system very differently than uniform radiation therapy [34–36]. We intend to conduct similar studies to identify associations between dosimetric parameters and these indirect effects of SFRT in the future. SFRT dosimetric association with treatment tumor response Valley dose and tumor EUD. The importance of tumor minimum dose to tumor control has long been established in conventional radiation therapy [37]. Does the same association between tumor control and minimum/valley dose hold for SFRT? For some the answer is yes and sophisticated techniques have been developed to “fill up” the dose valleys in an MRT beam by interlacing the microbeams from MRT from different irradiation angles. As a result, a uniform dose distribution inside the tumor is reached [38] while the surrounding normal tissue out of the “cross-firing” range still receive largely MRT radiation pattern of peaks and valleys. In a synchrotron-MRT study Ibahim et al. [39] reported that valley dose is closely correlated with cell survival, but valley dose alone does not determine the observed radiobiological effects. Our study shows that the tumor EUD (a = -10) and valley (minimum) tumor dose have the highest linear associations (R2 = 0.7923, F-stat = 15.26*; R2 = 0.7636, F-stat = 12.92*, respectively) with tumor treatment response (Fig 7 and S1 Table). This observed association between tumor treatment response with tumor valley minimum) dose and tumor EUD dose in this preclinical study is consistent with their known association in tumor treatment response seen in clinical conventional uniform dose radiation therapy. Our data suggest that valley (minimum) dose or Tumor EUD are more appropriate than peak dose for SFRT treatment prescription. When tumor control is the endpoint, we suggest that equal valley or minimum dose be used for comparative study between a uniform radiation and SFRT therapy or among different SFRT treatments. PVDR. Our data showed that PVDR has a consistent but not statistically significant association with tumor treatment response (R2 = 0.7194, F-stat = 7.691) (Fig 7, S1 Table). The linear regression analysis on day 17 was not statistically significant. The CoxPH analysis using the entire survival data set show a modest association with survival. Although not statistically significant, an inverse association is observed between PVDR value and survival fraction on Day 17—the higher PVDR value the less survival fraction. The inverse association is largely determined by the uniform radiation arm where PVDR value is 1.0. If this data point is removed, the PVDR association with survival for all SFRT arms is inconclusive (Fig 7). We believe this result of inverse association is likely biased by the study design that has very limited PVDR values (4 values) and strong cross-correlations between PVDR and other SFRT parameters (see Cross-Correlations discussion and Table 2, above). A better understanding of PVDR’s association with a given treatment response requires a carefully designed new study that focuses on the impact of PVDR value on treatment response. Percentage volume irradiated, peak width, and valley width. It seems logical that tumor treatment response is closely associated with the tumor volume irradiated. However, this is not supported by a clinical GRID therapy study by Neuner et al. [2] where both MLC-based and collimator-based GRID treatments showed similar response rates for pain, mass effect, other patient complaints, and have similar adverse reactions. The collimator-generated GRID had 50% of the radiation field open while the MLC-generated GRID had only 31% open. In our study the 20GyHalfSFRT and 20Gy2mmSFRT arms have similar percentage-volume-irradiated (as well as PDD curves) but there is a difference of 5 days in the 50% survival time (Fig 1 and Fig 6). Nonetheless, our data shows that percentage-volume-irradiated has the 3rd highest linear association (R2 = 0.7153, F-stat = 10.05*) with tumor treatment response (Fig 7 and Table 3). Since percentage-volume-irradiated is jointly determined by peak width and valley width it is understandable to see moderate associations between tumor treatment response and peak width (R2 = 0.4201, F-stat = 2.898 (not sig.)) and valley width (R2 = 0.4296, F-stat = 3.012 (not sig.)). In a synchrotron microbeam brain study using multiple beams Serduc et al. kept valley dose constant while varying peak width and peak dose. They concluded that the latter two parameters have strong influence therapeutic ratio [40]. Volume-averaged dose and peak dose. This study is designed to scrutinize the association of volume-averaged dose with tumor treatment response (Fig 1). The four study arms sharing very similar volume-averaged doses (20 or 18 Gy) exhibited very different tumor treatment responses (Figs 6 and 7) showing the survival rate at day 17 varied from 100% to 33%. Therefore, the association between volume-average dose and tumor treatment response is weak. We found that peak dose has little to no association with tumor treatment response (R2 = 0.04472, F-stat = 0.6874 (not sig.)) (Fig 7, S1 Table, Table 3). This finding is significant because peak dose has been used for treatment prescription in practically all SFRT treatments (8) (9). Although the linear regression analysis on day 17 showed a weak association between peak dose and survival that was not statistically significant, the CoxPH analysis using the entire survival data set did show a modest association with survival. Valley dose and tumor EUD. The importance of tumor minimum dose to tumor control has long been established in conventional radiation therapy [37]. Does the same association between tumor control and minimum/valley dose hold for SFRT? For some the answer is yes and sophisticated techniques have been developed to “fill up” the dose valleys in an MRT beam by interlacing the microbeams from MRT from different irradiation angles. As a result, a uniform dose distribution inside the tumor is reached [38] while the surrounding normal tissue out of the “cross-firing” range still receive largely MRT radiation pattern of peaks and valleys. In a synchrotron-MRT study Ibahim et al. [39] reported that valley dose is closely correlated with cell survival, but valley dose alone does not determine the observed radiobiological effects. Our study shows that the tumor EUD (a = -10) and valley (minimum) tumor dose have the highest linear associations (R2 = 0.7923, F-stat = 15.26*; R2 = 0.7636, F-stat = 12.92*, respectively) with tumor treatment response (Fig 7 and S1 Table). This observed association between tumor treatment response with tumor valley minimum) dose and tumor EUD dose in this preclinical study is consistent with their known association in tumor treatment response seen in clinical conventional uniform dose radiation therapy. Our data suggest that valley (minimum) dose or Tumor EUD are more appropriate than peak dose for SFRT treatment prescription. When tumor control is the endpoint, we suggest that equal valley or minimum dose be used for comparative study between a uniform radiation and SFRT therapy or among different SFRT treatments. PVDR. Our data showed that PVDR has a consistent but not statistically significant association with tumor treatment response (R2 = 0.7194, F-stat = 7.691) (Fig 7, S1 Table). The linear regression analysis on day 17 was not statistically significant. The CoxPH analysis using the entire survival data set show a modest association with survival. Although not statistically significant, an inverse association is observed between PVDR value and survival fraction on Day 17—the higher PVDR value the less survival fraction. The inverse association is largely determined by the uniform radiation arm where PVDR value is 1.0. If this data point is removed, the PVDR association with survival for all SFRT arms is inconclusive (Fig 7). We believe this result of inverse association is likely biased by the study design that has very limited PVDR values (4 values) and strong cross-correlations between PVDR and other SFRT parameters (see Cross-Correlations discussion and Table 2, above). A better understanding of PVDR’s association with a given treatment response requires a carefully designed new study that focuses on the impact of PVDR value on treatment response. Percentage volume irradiated, peak width, and valley width. It seems logical that tumor treatment response is closely associated with the tumor volume irradiated. However, this is not supported by a clinical GRID therapy study by Neuner et al. [2] where both MLC-based and collimator-based GRID treatments showed similar response rates for pain, mass effect, other patient complaints, and have similar adverse reactions. The collimator-generated GRID had 50% of the radiation field open while the MLC-generated GRID had only 31% open. In our study the 20GyHalfSFRT and 20Gy2mmSFRT arms have similar percentage-volume-irradiated (as well as PDD curves) but there is a difference of 5 days in the 50% survival time (Fig 1 and Fig 6). Nonetheless, our data shows that percentage-volume-irradiated has the 3rd highest linear association (R2 = 0.7153, F-stat = 10.05*) with tumor treatment response (Fig 7 and Table 3). Since percentage-volume-irradiated is jointly determined by peak width and valley width it is understandable to see moderate associations between tumor treatment response and peak width (R2 = 0.4201, F-stat = 2.898 (not sig.)) and valley width (R2 = 0.4296, F-stat = 3.012 (not sig.)). In a synchrotron microbeam brain study using multiple beams Serduc et al. kept valley dose constant while varying peak width and peak dose. They concluded that the latter two parameters have strong influence therapeutic ratio [40]. Volume-averaged dose and peak dose. This study is designed to scrutinize the association of volume-averaged dose with tumor treatment response (Fig 1). The four study arms sharing very similar volume-averaged doses (20 or 18 Gy) exhibited very different tumor treatment responses (Figs 6 and 7) showing the survival rate at day 17 varied from 100% to 33%. Therefore, the association between volume-average dose and tumor treatment response is weak. We found that peak dose has little to no association with tumor treatment response (R2 = 0.04472, F-stat = 0.6874 (not sig.)) (Fig 7, S1 Table, Table 3). This finding is significant because peak dose has been used for treatment prescription in practically all SFRT treatments (8) (9). Although the linear regression analysis on day 17 showed a weak association between peak dose and survival that was not statistically significant, the CoxPH analysis using the entire survival data set did show a modest association with survival. SFRT dosimetric association with normal tissue toxicity We did not study treatment induced normal tissue toxicity directly in this study. We used body weight change post radiation (targeted to the flank, lower abdominal region of the animal) as an indicator, not an evidence of normal tissue toxicity. We did not see a strong association between animal body weight change and any of the eight dosimetric parameters studied, except a modest association with valley (minimum) dose. Valley dose. The strongest association we observed is a weak one between body weight change and valley (minimum) dose (R2 = 0.3814, F-stat = 13.56**) (Table 4). Note that valley (minimum) dose is also strongly associated with tumor treatment response (R2 = 0.7636, F-stat = 12.92*). Our finding is consistent with a normal mouse brain MRT study Nakayma et al. reported that valley dose is one of the important factors to determine normal brain dose tolerance [41]. Our data suggests that among the nine dosimetric parameters we have studied, valley dose may have a close correlation with both tumor control and toxicity, and thus is a crucial dosimetric parameter in SFRT treatment. Valley width, peak width, percentage volume irradiated. The valley width, peak width, and percentage volume of the tumor that is irradiated parameters were only weakly associated with animal body weight change post radiation (R2 = 0.2853, F-stat = 8.783**; R2 = 0.2759, F-stat = 8.382**; and R2 = 0.1985, F-stat = 5.448*, respectively) (Fig 8 and Table 4). Note that in this study peak width and valley width are closely correlated (see discussion on cross-correlations, Table 2). Percentage volume directly irradiated showed no statistically important association with body weight change. Neuner et al. reported that they observed similar treatment responses from clinical GRID treatments of different percentages of volume directly irradiated [2]. Normal tissue EUD, PVDR, volume-averaged dose, peak dose. The normal tissue EUD, PVDR, volume-averaged dose, and peak dose parameters showed little to no association with body weight change post radiation (R2 = 1.022x10-03, F-stat = 0.882 (not sig.); R2 = 0.1203, F-stat = 3.009 (not sig.); R2 = 0.03308, F-stat = .7526 (not sig.); and R2 = 5.99x10-06, F-stat = 1.32x10-04 (not sig.), respectively). Our finding is consistent with a rat normal brain minibeam study by Prezado et al. showing arms with similar volume-average-doses have drastic differences in survival (14) and inconsistent with a MRT study on normal mouse skin by Priyadarshika et al. concluded that integrated dose (i.e., volume-averaged dose) rather than peak or valley dose, may dictate the acute skin toxicity [42]. 2mm wide beam array SFRT. Our data indicates that the 20Gy2mmSFRT arm is not only the most relevant to clinical application because of its millimeter scale, but it also has the potential for superior therapeutic ratio. The 20Gy2mmSFRT arm showed similar survival with the 50GySFRT arm but has significantly lower valley dose (6.2 Gy vs. 17 Gy). At the same time, it showed the least, if any, body weight change compared to the untreated arm while the 50GySFRT arm with 0.31mm beam width exhibited significant body weight growth deficit (Table 1 and Fig 6). The 20GyUniform arm has the best tumor treatment response, but also had a large drop in body weight. Our data indicated that for KV x-rays the 2mm wide beam array has the potential to be the SFRT pattern for high therapeutic ratio and deserves further investigation. Valley dose. The strongest association we observed is a weak one between body weight change and valley (minimum) dose (R2 = 0.3814, F-stat = 13.56**) (Table 4). Note that valley (minimum) dose is also strongly associated with tumor treatment response (R2 = 0.7636, F-stat = 12.92*). Our finding is consistent with a normal mouse brain MRT study Nakayma et al. reported that valley dose is one of the important factors to determine normal brain dose tolerance [41]. Our data suggests that among the nine dosimetric parameters we have studied, valley dose may have a close correlation with both tumor control and toxicity, and thus is a crucial dosimetric parameter in SFRT treatment. Valley width, peak width, percentage volume irradiated. The valley width, peak width, and percentage volume of the tumor that is irradiated parameters were only weakly associated with animal body weight change post radiation (R2 = 0.2853, F-stat = 8.783**; R2 = 0.2759, F-stat = 8.382**; and R2 = 0.1985, F-stat = 5.448*, respectively) (Fig 8 and Table 4). Note that in this study peak width and valley width are closely correlated (see discussion on cross-correlations, Table 2). Percentage volume directly irradiated showed no statistically important association with body weight change. Neuner et al. reported that they observed similar treatment responses from clinical GRID treatments of different percentages of volume directly irradiated [2]. Normal tissue EUD, PVDR, volume-averaged dose, peak dose. The normal tissue EUD, PVDR, volume-averaged dose, and peak dose parameters showed little to no association with body weight change post radiation (R2 = 1.022x10-03, F-stat = 0.882 (not sig.); R2 = 0.1203, F-stat = 3.009 (not sig.); R2 = 0.03308, F-stat = .7526 (not sig.); and R2 = 5.99x10-06, F-stat = 1.32x10-04 (not sig.), respectively). Our finding is consistent with a rat normal brain minibeam study by Prezado et al. showing arms with similar volume-average-doses have drastic differences in survival (14) and inconsistent with a MRT study on normal mouse skin by Priyadarshika et al. concluded that integrated dose (i.e., volume-averaged dose) rather than peak or valley dose, may dictate the acute skin toxicity [42]. 2mm wide beam array SFRT. Our data indicates that the 20Gy2mmSFRT arm is not only the most relevant to clinical application because of its millimeter scale, but it also has the potential for superior therapeutic ratio. The 20Gy2mmSFRT arm showed similar survival with the 50GySFRT arm but has significantly lower valley dose (6.2 Gy vs. 17 Gy). At the same time, it showed the least, if any, body weight change compared to the untreated arm while the 50GySFRT arm with 0.31mm beam width exhibited significant body weight growth deficit (Table 1 and Fig 6). The 20GyUniform arm has the best tumor treatment response, but also had a large drop in body weight. Our data indicated that for KV x-rays the 2mm wide beam array has the potential to be the SFRT pattern for high therapeutic ratio and deserves further investigation. Summary In this conventional dose rate small animal SFRT study we used a large range of radiation spatial fractionation scales to study the association of dosimetric parameters with treatment response. Among the SFRT dosimetric parameters studied we found that valley (minimum) dose, tumor EUD and percentage tumor irradiated have strong and proportional associations with tumor treatment response while peak dose, which has been routinely used to prescribe SFRT treatments, exhibited little association. We found weak associations between body weight change and the dosimetric parameters studied. SFRT dosimetry is significantly more complex than that of conventional radiation therapy. Therefore, identification of key SFRT dosimetric parameters that are closely associated with treatment responses is crucial to the advancement of SFRT research, clinical application, and treatment technology development. Supporting information S1 Fig. Illustration of treatment delivery verification analysis by film. (A) The post-treatment verification film for a 20GyHalfSFRT treated tumor shows that only one-half the tumor was treated as intended. The black dashed line in the photograph was drawn to illustrate which half of the tumor was irradiated. (B) The verification films for all 5 animals included in the study arm were analyzed by calculating the percentage area of the tumor irradiated. https://doi.org/10.1371/journal.pone.0229053.s001 (TIF) S2 Fig. Tumor volumes changes post-radiation for all six study arms. Tumor volumes for all animals within each treatment arm are averaged and then plotted over time in units of mm3. Error bars represent one standard deviation from the mean. https://doi.org/10.1371/journal.pone.0229053.s002 (TIF) S1 Table. Univariate linear regression analysis of survival day 17. This is the full table of coefficients for the corresponding linear regression models used in Fig 7. We analyze 8 models with single covariates, one for each dosimetric parameter and list their corresponding statistics. Tumor EUD and Valley Dose have the largest magnitude of effect on Survival (Day 17) and together with % Volume Irradiated are statistically significant. However, analyzing data for a single timepoint (Day17) is limited by animal losses at Day 17 (missing data), so we include a more robust statistical model that utilizes all the data in Table 4. https://doi.org/10.1371/journal.pone.0229053.s003 (TIF) S2 Table. Table of statistical significance values for the percentage animal surviving. This is a table of statistical significance values for the animal survival curves, shown in Fig 6. https://doi.org/10.1371/journal.pone.0229053.s004 (TIF) S1 File. The ARRIVE guidelines checklist. https://doi.org/10.1371/journal.pone.0229053.s005 (PDF) S1 Data. https://doi.org/10.1371/journal.pone.0229053.s006 (CSV) Acknowledgments One of the authors (Chang) acknowledges Dr. Mark W. Dewhirst for his decade long unwavering encouragements and expert radiobiology advices, which are invaluable for this (and other) original work exploring the magic of radiation spatial fractionation. One of authors (Rivera) expresses her appreciation to Leith Rankine, MS for his kind help on EUD calculation.
TI - Conventional dose rate spatially-fractionated radiation therapy (SFRT) treatment response and its association with dosimetric parameters—A preclinical study in a Fischer 344 rat model
JF - PLoS ONE
DO - 10.1371/journal.pone.0229053
DA - 2020-06-22
UR - https://www.deepdyve.com/lp/public-library-of-science-plos-journal/conventional-dose-rate-spatially-fractionated-radiation-therapy-sfrt-300WSZhdg0
SP - e0229053
VL - 15
IS - 6
DP - DeepDyve
ER -