Dose-volume effects in rat thoracolumbar spinal cord: an evaluation of NTCP models

Characterizing Proton-Induced Biological Effects in a Mouse Spinal Cord Model: A Comparison of Bragg Peak and Entrance Beam Response in Single and Fractionated Exposures

PURPOSE

Proton relative biological effectiveness (RBE) is a dynamic variable influenced by factors like linear energy transfer (LET), dose, tissue type, and biological endpoint. The standard fixed proton RBE of 1.1, currently used in clinical planning, may not accurately represent the true biological effects of proton therapy (PT) in all cases. This uncertainty can contribute to radiation-induced normal tissue toxicity in patients. In late-responding tissues such as the spinal cord, toxicity can cause devastating complications. This study investigated spinal cord tolerance in mice subjected to proton irradiation and characterized the influence of fractionation on proton- induced myelopathy at entrance (ENT) and Bragg peak (BP) positions.

METHODS AND MATERIALS

Cervical spinal cords of 8-week-old C57BL/6J female mice were irradiated with single- or multi-fractions (18x) using lateral opposed radiation fields at 1 of 2 positions along the Bragg curve: ENT (dose-mean LET = 1.2 keV/μm) and BP (LET = 6.9 keV/μm). Mice were monitored over 1 year for changes in weight, mobility, and general health, with radiation-induced myelopathy as the primary biological endpoint. Calculations of the RBE of the ENT and BP curve (RBEENT/BP) were performed.

RESULTS

Single-fraction RBEENT/BP for 50% effect probability (tolerance dose (TD50), grade II paresis, determined using log-logistic model fitting) was 1.10 ± 0.06 (95% CI) and for multifraction treatments it was 1.19 ± 0.05 (95% CI). Higher incidence and faster onset of paralysis were seen in mice treated at the BP compared with ENT.

CONCLUSIONS

The findings challenge the universally fixed RBE value in PT, indicating up to a 25% mouse spinal cord RBEENT/BP variation for multifraction treatments. These results highlight the importance of considering fractionation in determining RBE for PT. Robust characterization of proton-induced toxicity, aided by in vivo models, is paramount for refining clinical decision-making and mitigating potential patient side effects.

Primary Hypothyroidism in Childhood Cancer Survivors Treated With Radiation Therapy: A PENTEC Comprehensive Review

PURPOSE

From the Pediatric Normal Tissue Effects in the Clinic (PENTEC) initiative, a systematic review and meta-analysis of publications reporting on radiation dose-volume effects for risk of primary hypothyroidism after radiation therapy for pediatric malignancies was performed.

METHODS AND MATERIALS

All studies included childhood cancer survivors, diagnosed at age <21 years, whose radiation therapy fields exposed the thyroid gland and who were followed for primary hypothyroidism. Children who received pituitary-hypothalamic or total-body irradiation were excluded. PubMed and the Cochrane Library were searched for studies published from 1970 to 2017. Data on age at treatment, patient sex, radiation dose to neck or thyroid gland, specific endpoints for hypothyroidism that were used in the studies, and reported risks of hypothyroidism were collected. Radiation dose-volume effects were modeled using logistic dose response. Relative excess risk of hypothyroidism as a function of age at treatment and sex was assessed by meta-analysis of reported relative risks (RR) and odds ratios.

RESULTS

Fifteen publications (of 1709 identified) were included for systematic review. Eight studies reported data amenable for dose-response analysis. At mean thyroid doses of 10, 20, and 30 Gy, predicted rates of uncompensated (clinical) hypothyroidism were 4%, 7%, and 13%, respectively. Predicted rates of compensated (subclinical) hypothyroidism were 12%, 25%, and 44% after thyroid doses of 10, 20, and 30 Gy, respectively. Female sex (RR = 1.7, P < .0001) and age >15 years at radiation therapy (RR = 1.3, P = .005) were associated with higher risks of hypothyroidism. After a mean thyroid dose of 20 Gy, predicted risks of hypothyroidism were 13% for males <14 years of age, increasing to 29% for females >15 years of age.

CONCLUSION

A radiation dose response for risk of hypothyroidism is evident; a threshold radiation dose associated with no risk is not observed. Thyroid dose exposure should be minimized when feasible. Data on hypothyroidism after radiation therapy should be better reported to facilitate pooled analyses.

Biomedical advances and future prospects of high-precision three-dimensional radiotherapy and four-dimensional radiotherapy

Biomedical advances of external-beam radiotherapy (EBRT) with improvements in physical accuracy are reviewed. High-precision (±1 mm) three-dimensional radiotherapy (3DRT) can utilize respective therapeutic open doors in the tumor control probability curve and in the normal tissue complication probability curve instead of the one single therapeutic window in two-dimensional EBRT. High-precision 3DRT achieved higher tumor control and probable survival rates for patients with small peripheral lung and liver cancers. Four-dimensional radiotherapy (4DRT), which can reduce uncertainties in 3DRT due to organ motion by real-time (every 0.1-1 s) tumor-tracking and immediate (0.1-1 s) irradiation, have achieved reduced adverse effects for prostate and pancreatic tumors near the digestive tract and with similar or better tumor control. Particle beam therapy improved tumor control and probable survival for patients with large liver tumors. The clinical outcomes of locally advanced or multiple tumors located near serial-type organs can theoretically be improved further by integrating the 4DRT concept with particle beams.

Can dose convolution modelling explain bath and shower effects in rat spinal cord?

Objective. ‘Bath and shower’ effects were first seen in proton irradiations of rat spinal cord, where a low dose ‘bath’ reduced the smaller field ‘shower’ dose needed for limb paralysis giving the appearance of sensitisation of the cord or disproportionate response. This was difficult to reconcile with existing tissue complication models. The purpose of this investigation is to explore a different approach using a dose convolution algorithm to model the 50% isoeffect endpoint. Approach. Bath and shower dose distributions were convolved with Gaussian functions with widths specified by the σ parameter. The hypothesis was that the maximum value from the convolved distributions was constant for isoeffect across the modelled scenarios. A simpler field length dependent relative biological effectiveness (FLRBE) approach was also used for a subset of the data which gave results independent of σ. Main results. The maximum values from the convolved distributions were constant within ±17% across the bath and shower experiments for σ = 3.5 mm, whereas the maximum dose varied by a factor of four. The FLRBE results were also within ±14% confirming the validity of the dose convolution approach. Significance. A simple approach using dose convolution modelling of the 50% isotoxicity gave compelling consistency with the full range of bath and shower results, while the FLRBE approach confirmed the results for the symmetric field data. Convolution modelling and the effect of time interval were consistent with a signalling factor diffusion mechanism such as the ‘bystander effect’. The results suggest biological effectiveness is reduced for very small field sizes, requiring a higher isoeffect dose. By implication, the bath dose does not sensitise the cord to the shower dose; when biological effectiveness is accounted for, a small increase in the bath dose requires a significantly larger reduction in shower dose for isoeffect.

Brachial Plexus Tolerance to Single-Session SAbR in a Pig Model

PURPOSE

The single-session dose tolerance of the spinal nerves has been observed to be similar to that of the spinal cord in pigs, counter to the perception that peripheral nerves are more tolerant to radiation. This pilot study aims to obtain a first impression of the single-session dose-response of the brachial plexus using pigs as a model.

METHODS AND MATERIALS

Ten Yucatan minipigs underwent computed tomography and magnetic resonance imaging for treatment planning, followed by single-session stereotactic ablative radiotherapy. A 2.5-cm length of the left-sided brachial plexus cords was irradiated. Pigs were distributed in 3 groups with prescription doses of 16 (n = 3), 19 (n = 4), and 22 Gy (n = 3). Neurologic status was assessed by observation for changes in gait and electrodiagnostic examination. Histopathologic examination was performed with light microscopy of paraffin-embedded sections stained with Luxol fast blue/periodic acid-Schiff and Masson's trichrome.

RESULTS

Seven of the 10 pigs developed motor deficit to the front limb of the irradiated side, with a latency from 5 to 8 weeks after irradiation. Probit analysis of the maximum nerve dose yields an estimated ED₅₀ of 19.3 Gy for neurologic deficit, but the number of animals was insufficient to estimate 95% confidence intervals. No motor deficits were observed at a maximum dose of 17.6 Gy for any pig. Nerve conduction studies showed an absence of sensory response in all responders and absent or low motor response in most of the responders (71%). All symptomatic pigs showed histologic lesions to the left-sided plexus consistent with radiation-induced neuropathy.

CONCLUSIONS

The single-session ED₅₀ for symptomatic plexopathy in Yucatan minipigs after irradiation of a 2.5-cm length of the brachial plexus cords was determined to be 19.3 Gy. The dose-response curve overlaps that of the spinal nerves and the spinal cord in the same animal model. The relationship between the brachial plexus tolerance in pigs and humans is unknown, and caution is warranted when extrapolating for clinical use.

Effects From Nonuniform Dose Distribution in the Spinal Nerves of Pigs: Analysis of Normal Tissue Complication Probability Models

PURPOSE

Our purpose was to evaluate normal tissue complication probability (NTCP) models for their ability to describe the increase in tolerance as the length of irradiated spinal nerve is reduced in a pig.

METHODS AND MATERIALS

Common phenomenological and semimechanistic NTCP models were fit using the maximum likelihood estimate method to dose-response data from spinal nerve irradiation studies in pigs. Statistical analysis was used to compare how well each model fit the data. Model parameters were then applied to a previously published dose distribution used for spinal cord irradiation in rats under the assumption of a similar dose-response.

RESULTS

The Lyman-Kutcher-Burman model, relative seriality, and critical volume model fit the spinal nerve data equally well, but the mean dose logistic and relative seriality models gave the best fit after penalizing for the number of model parameters. The minimum dose logistic regression model was the only model showing a lack of fit. When extrapolated to a 0.5-cm simulated square-wave-like dose distribution, the serial behaving models showed negligible increase in dose-response curve. The Lyman-Kutcher-Burman model and relative seriality models showed significant shifting of NTCP curves due to parallel behaving parameters. The critical volume model gave the closest match to the rat data.

CONCLUSIONS

Several phenomenological and semimechanistic models were observed to adequately describe the increase in the radiation tolerance of the spinal nerves when changing the irradiated length from 1.5 to 0.5 cm. Contrary to common perception, model parameters suggest parallel behaving tissue architecture. Under the assumption that the spinal nerve response to radiation is similar to that of the spinal cord, only the critical volume model was robust when extrapolating to outcome data from a 0.5-cm square-wave-like dose distribution, as was delivered in rodent spinal cord irradiation research.

Existence of a Dose-Length Effect in Spinal Nerves Receiving Single-Session Stereotactic Ablative Radiation Therapy

PURPOSE

The spinal nerves have been observed to have a similar single-session dose tolerance to that of the spinal cord in pigs. Small-animal studies have shown that spinal cord dose tolerance depends on the length irradiated. This work aims to determine whether a dose-length effect exists for spinal nerves.

METHODS AND MATERIALS

Twenty-seven Yucatan minipigs underwent computed tomography and magnetic resonance imaging for treatment planning, followed by single-session stereotactic ablative radiation therapy. A 0.5 cm length of the left-sided C6, C7, and C8 spinal nerves was targeted. The pigs were distributed into 6 groups with prescription doses of 16 Gy (n = 5), 18 Gy (n = 5), 20 Gy (n = 5), 22 Gy (n = 5), 24 Gy (n = 5), or 36 Gy (n = 2) and corresponding maximum doses of 16.7, 19.1, 21.3, 23.1, 25.5, and 38.6 Gy, respectively. Neurologic status was assessed with a serial electrodiagnostic examination and daily observation of gait for approximately 52 weeks. A histopathologic examination of paraffin-embedded sections with Luxol fast blue/periodic acid-Schiff's staining was also performed.

RESULTS

Marked gait change was observed in 8 of 27 irradiated pigs. The latency for responding pigs was 11 to 16 weeks after irradiation. The affected animals presented with a limp in the left front limb, and 62.5% of these pigs had electrodiagnostic evidence of denervation in the C6 and C7 innervated muscles. A probit analysis showed the dose associated with a 50% incidence of gait change is 23.9 Gy (95% confidence interval, 22.5-25.8 Gy), which is 20% higher than that reported in a companion study where a 1.5 cm length was irradiated. All symptomatic pigs had demyelination and fibrosis in the irradiated nerves, but the contralateral nerves and spinal cord were normal.

CONCLUSIONS

A dose-length effect was observed for single-session irradiation of the spinal nerves in a Yucatan minipig model.

Modeling Radiotherapy Induced Normal Tissue Complications: An Overview beyond Phenomenological Models

An overview of radiotherapy (RT) induced normal tissue complication probability (NTCP) models is presented. NTCP models based on empirical and mechanistic approaches that describe a specific radiation induced late effect proposed over time for conventional RT are reviewed with particular emphasis on their basic assumptions and related mathematical translation and their weak and strong points.

Comparison of dose response models for predicting normal tissue complications from cancer radiotherapy: application in rat spinal cord

Seven different radiobiological dose-response models have been compared with regard to their ability to describe experimental data. The first four models, namely the critical volume, the relative seriality, the inverse tumor and the critical element models are mainly based on cell survival biology. The other three models: the Lyman (Gaussian distribution), the parallel architecture and the Weibull distribution models are semi-empirical and rather based on statistical distributions. The maximum likelihood estimation was used to fit the models to experimental data and the χ2-distribution, AIC criterion and F-test were applied to compare the goodness-of-fit of the models. The comparison was performed using experimental data for rat spinal cord injury. Both the shape of the dose-response curve and the ability of handling the volume dependence were separately compared for each model. All the models were found to be acceptable in describing the present experimental dataset (p > 0.05). For the white matter necrosis dataset, the Weibull and Lyman models were clearly superior to the other models, whereas for the vascular damage case, the Relative Seriality model seems to have the best performance although the Critical volume, Inverse tumor, Critical element and Parallel architecture models gave similar results. Although the differences between many of the investigated models are rather small, they still may be of importance in indicating the advantages and limitations of each particular model. It appears that most of the models have favorable properties for describing dose-response data, which indicates that they may be suitable to be used in biologically optimized intensity modulated radiation therapy planning, provided a proper estimation of their radiobiological parameters had been performed for every tissue and clinical endpoint.

Spinal cord tolerance in the age of spinal radiosurgery: lessons from preclinical studies

Clinical implementation of spinal radiosurgery has increased rapidly in recent years, but little is known regarding human spinal cord tolerance to single-fraction irradiation. In contrast, preclinical studies in single-fraction spinal cord tolerance have been ongoing since the 1970s. The influences of field length, dose rate, inhomogeneous dose distributions, and reirradiation have all been investigated. This review summarizes literature regarding single-fraction spinal cord tolerance in preclinical models with an emphasis on practical clinical significance. The outcomes of studies that incorporate uniform irradiation are surprisingly consistent among multiple small- and large-animal models. Extensive investigation of inhomogeneous dose distributions in the rat has demonstrated a significant dose-volume effect while preliminary results from one pig study are contradictory. Preclinical spinal cord dose-volume studies indicate that dose distribution is more critical than the volume irradiated suggesting that neither dose-volume histogram analysis nor absolute volume constraints are effective in predicting complications. Reirradiation data are sparse, but results from guinea pig, rat, and pig studies are consistent with the hypothesis that the spinal cord possesses a large capacity for repair. The mechanisms behind the phenomena observed in spinal cord studies are not readily explained and the ability of dose response models to predict outcomes is variable underscoring the need for further investigation. Animal studies provide insight into the phenomena and mechanisms of radiosensitivity but the true significance of animal studies can only be discovered through clinical trials.

Inter- and intramouse heterogeneity of radiation response for a growing paired organ

An intensive search for predictive markers of individual radiation response of apparently normal tissues in cancer patients is in progress at the genetic and epigenetic levels. However, the relative impact of variability at these levels is not clear. Experimental results obtained in inbred rodents, which have significantly reduced genetic heterogeneity relative to a population of human patients, may help to clarify this issue. We investigated a paired-organ mouse system in a strain of inbred mice to evaluate the intermouse variability of normal tissue radiation response, singled out from measurement errors and stochastic effects. The legs of 5-day-old C3H mice were homogeneously gamma-irradiated with a range of single doses. The lengths of the right and left tibiae were measured in 30 kVp X-ray images taken at the time of irradiation and at 84 days postirradiation. The dose-effect curves were smooth and well defined, with bone growth retardation evident at approximately 14 Gy and higher, and were marginally gender-dependent. The intramouse (left compared to right) variability of the tibia length on day 89, which characterized stochastic effects, was not distinguishable from the measurement error for doses less than 16-18 Gy and slightly exceeded measurement errors only at the largest doses of 20-22 Gy. The corresponding intermouse variability was greater than the measurement error and stochastic effects at all doses used. Interestingly, the total variability, judged by the gamma(50) values of approximately 7 we obtained, was similar to that reported for severe late reactions in human normal tissue. If the variations of response determined by epigenetic events in human patients free of known factors associated with altered radiation sensitivity are comparable to those observed in this mouse model, our results imply a relatively low power of genetic approaches alone to predict individual side effects in radiotherapy.

HI-CHART: a phase I/II study on the feasibility of high-dose continuous hyperfractionated accelerated radiotherapy in patients with inoperable non-small-cell lung cancer

PURPOSE

To determine the feasibility of high-dose continuous hyperfractionated accelerated radiotherapy in patients with inoperable non-small-cell lung cancer (NSCLC).

PATIENTS AND METHODS

In a prospective, Phase I/II study, according to the risk for radiation pneumonitis, three risk groups were defined: V(20) <25%, V(20) 25-37%, and V(20) >37%. The dose was administered in three steps from 61.2 Gy/34 fractions/23 days to 64.8 Gy/36 fractions/24 days to 68.40 Gy/38 fractions/25 days (1.8 Gy b.i.d. with 8-h interval), using a three-dimensional conformal technique. Only the mediastinal lymph node areas that were positive on the pretreatment (18)F-deoxy-D-glucose positron emission tomography scan were included in the target volume. The primary endpoint was toxicity.

RESULTS

A total of 48 Stage I-IIIB patients were included. In all risk groups, 68.40 Gy/38 fractions/25 days could be administered. Maximal toxicity according to the risk groups was as follows: V(20) <25% (n = 35): 1 Grade 4 (G4) lung and 1 G3 reversible esophageal toxicity; V(20) 35-37% (n = 12): 1 G5 lung and 1 G3 reversible esophageal toxicity. For the whole group, local tumor recurrence occurred in 25% (95% confidence interval 14%-40%) of the patients, with 1 of 48 (2.1%; upper one-sided 95% confidence limit 9.5%) having an isolated nodal recurrence. The median actuarial overall survival was 20 months, with a 2-year survival rate of 36%.

CONCLUSIONS

High-dose continuous hyperfractionated accelerated radiotherapy up to a dose of 68.40 Gy/38 fractions/25 days (a biologic equivalent of approximately 80 Gy when delivered in conventional fractionation) in patients with inoperable NSCLC and a V(20) up to 37% is feasible.

Dose-volume effects in rat thoracolumbar spinal cord: the effects of nonuniform dose distribution

PURPOSE

To investigate dose-volume effects in rat spinal cord irradiated with nonuniform dose distributions and to assess regional differences in radiosensitivity.

METHODS AND MATERIALS

A total of 106 rats divided into three groups were irradiated with (192)Ir gamma-rays at a high dose rate. The groups were irradiated with one, two, or six catheters distributed around the thoracolumbar spinal cord to create different dose distributions. After irradiation, the animals were tested for motor function for 9 months. The response was defined as motor dysfunction and WM or nerve root necrosis. Dose-response data were analyzed with a probit analysis as function of the dose level at a percentage of the volume (D(%)) and with different normal tissue complication probability models. Additionally, the histologic responses of the individual dose voxels were analyzed after registration with the histologic sections.

RESULTS

The probit analysis at D(24) (24% of the volume) gave the best fit results. In addition, the Lyman Kutcher Burman model and the relative seriality model showed acceptable fits, with volume parameters of 0.17 and 0.53, respectively. The histology-based analysis revealed a lower radiosensitivity for the dorsal (50% isoeffective dose [ED(50)] = 32.3) and lateral WM (ED(50) = 33.7 Gy) compared with the dorsal (ED(50) = 25.9 Gy) and ventral nerve roots (ED(50) = 24.1 Gy).

CONCLUSIONS

For this nonuniform irradiation, the spinal cord did not show typical serial behavior. No migration terms were needed for an acceptable fit of the dose-response curves. A higher radiosensitivity for the lumbar nerve roots than for the thoracic WM was found.

Toward improving the therapeutic ratio in stereotactic radiosurgery: selective modulation of the radiation responses of both normal tissues and tumor

✓A review of the radiobiological factors that influence the response of the brain to radiation is provided in relation to stereotactic radiosurgery (SRS). The prospects for intervention after radiation treatment to selectively modulate the expression of late central nervous system (CNS) injury is considered, as well as an account of recent interest in the use of radiation enhancers to selectively increase the response of tumors to radiation.

Brain necrosis in humans, after conventional irradiation, indicates that the risk of necrosis increases rapidly after an equivalent single dose of 12 or 13 Gy. When single-dose treatments are extended due to 60Co decay or planned extension of treatment times, account should be taken of the effects of the repair of sublethal radiation damage to DNA on the efficacy of treatment. Both repair capacity and repair kinetics will also influence tumor control, but parameters to quantify this effect have not yet been established.

The volume of CNS tissue that has been irradiated affects the tissue response, but this effect is only significant for volumes less than 0.05 cm3. The gain obtained from irradiation of small volumes is reduced, however, when focal irradiation is given within a wider field of irradiation.

Based on a vascular hypothesis explaining the pathogenesis of late CNS damage, approaches designed to selectively modulate the frequency of late CNS damage have been validated. Given the high intrinsic radioresistance of some tumors, as opposed to the presence of hypoxia, an interest has developed in the use of selective radiation enhancers in the treatment of tumors. The compound presently available has proved to be disappointing clinically due to toxicity at effective doses, when repeated administration is required. However, when given at high single doses it is less toxic and may be more effective. Less toxic radiation enhancers need to be developed.

Physical models and simpler dosimetric descriptors of radiation late toxicity

Predicting radiation damage to specific organs is becoming ever more challenging with the use of intensity-modulated beams, nonuniform dose distributions, partial organ irradiation, and interpatient and even intraorgan variations in radiation sensitivity. Data-based physical models can be of use in summarizing complicated dose-volume data to help describe clinical outcomes and ultimately aid in the prediction of clinical toxicity. This article attempts to provide a brief overview of the use of normal tissue complication probability (NTCP) models and other simple dose/volume metrics to describe a few clinically significant complications (either frequent or serious) associated with radiation therapy of the head and neck, thorax, and abdominal-pelvic regions. Specifically, it reviews the application of these methods for late toxicities of the parotid, lung, heart, spinal cord, liver, and rectum. It focuses on organ-specific NTCP parameters as well as simple dosimetric descriptors that might be used to help treatment plan evaluation in clinical practice.

A normal tissue dose response model of dynamic repair processes

A model is presented for serial, critical element complication mechanisms for irradiated volumes from length scales of a few millimetres up to the entire organ. The central element of the model is the description of radiation complication as the failure of a dynamic repair process. The nature of the repair process is seen as reestablishing the structural organization of the tissue, rather than mere replenishment of lost cells. The interactions between the cells, such as migration, involved in the repair process are assumed to have finite ranges, which limits the repair capacity and is the defining property of a finite-sized reconstruction unit. Since the details of the repair processes are largely unknown, the development aims to make the most general assumptions about them. The model employs analogies and methods from thermodynamics and statistical physics. An explicit analytical form of the dose response of the reconstruction unit for total, partial and inhomogeneous irradiation is derived. The use of the model is demonstrated with data from animal spinal cord experiments and clinical data about heart, lung and rectum. The three-parameter model lends a new perspective to the equivalent uniform dose formalism and the established serial and parallel complication models. Its implications for dose optimization are discussed.