- Open Access
The dose-response relationship for cardiovascular disease is not necessarily linear
Radiation Oncology volume 12, Article number: 74 (2017)
The probability for a complication after radiotherapy is usually a function of dose and volume in the organ or tissue of interest. In most epidemiological studies the risk for a complication is stratified in terms of dose, but not irradiated volume. We show that the obtained risk cannot generally be applied to radiotherapy patients.The epidemiological data of Darby et al. (N Engl J Med 368:2527, 2013) who found a linear relationship between the excess relative risk of major coronary events as function of mean heart dose in patients treated with tangential breast irradiation are analyzed. We have used the relative seriality model for a partly irradiated heart (“a lot to a little”) which models radiation therapy using two tangential fields. The relative seriality model was then used to predict NTCP of cardiovascular disease for a homogenously irradiated heart (“a little to a lot”). The relative seriality model was fitted to the data of Darby et al. (N Engl J Med 368:2527, 2013) for tangential breast irradiation. For the situation “a little to a lot” it was found that the dose-response relationship is sigmoidal and contradicts the findings of Darby et al. (N Engl J Med 368:2527, 2013). It was shown in this work that epidemiological studies which predict a linear dose-response relationship for cardiovascular disease can be reproduced by bio-physical models for normal tissue complication. For irradiation situations which were not included in the epidemiological studies, e.g. a homogenous irradiation of the heart (“a little to a lot”) the dose-response curve can be different. This could have consequences whether or not IMRT should be used for treating breast cancer. We believe that the results of epidemiological studies should not be generally used to predict normal tissue complications. It is better to use such data to optimize bio-physical models which can then be applied (with caution) to general treatment situations.
The results from epidemiological studies of side effects from radiation therapy can provide valuable information for the decision-making process in clinical treatment planning and can ultimately help to treat a patient successfully with radiation while minimizing adverse reactions. However, epidemiologically obtained data describe a side effect for a well-defined treatment situation and cannot be generally applied to all radiotherapy patients. The application of specifically obtained dose-response relationships to other treatment techniques/modalities requires the extrapolation of the obtained risks. This extrapolation can be achieved by the use of biophysical models.
The probability for a complication is usually a function of dose and volume in the organ or tissue of interest. In most epidemiological studies the risk for a complication is stratified in terms of dose, but not irradiated volume. We are concerned that this data, investigated for a very specific dose-volume situation, is used in clinical practice in a generalized way for completely different clinical settings.
In this report we want to point out the pits and pearls by using as an example the rate of major coronary events after radiation therapy. In particular we analyze the epidemiological data of Darby et al.  who found a linear relationship between the excess relative risk (ERR) of major coronary events as function of mean heart dose in patients treated with tangential breast irradiation. They obtained a risk factor of 7.4% per Gy of mean dose to the heart which is currently widely used in clinical practice. We believe that the general usage of this factor is dangerous, since it is a greatly simplified representation of complex biological processes which can only be applied safely to exactly the same situation for which the data were obtained (tangential breast irradiation). For this purpose we show that if the epidemiological data are used as input into a bio-physical model describing cardiovascular complications, the shape of the resulting dose-response relationship is highly dependent on the specific dose-volume situation. We believe further that the data of Darby et al.  should be better used for verifying and optimizing existing NTCP models of coronary events which can then be used in a more general setting.
Darby et al.  conducted a population-based case–control study of major coronary events (i.e., myocardial infarction, coronary revascularization, or death from ischemic heart disease) in 2168 women who underwent radiotherapy for breast cancer between 1958 and 2001 in Sweden and Denmark. The study included 963 women with major coronary events and 1205 controls. For each woman, the mean radiation doses to the whole heart were estimated from her radiotherapy chart. They found that excess rates of major coronary events increased linearly with the mean dose to the heart by 7.4% per Gy. This obtained dose response-relationship for major coronary events is illustrated by the symbols in Fig. 2 with the corresponding 95% confidence interval as a function of mean heart dose. The results of Darby et al. were validate by van den Bogaard et al.  in an independent cohort using individual 3D CT planning data.
Probability of a complication is function of dose and volume
The probability of a complication in radiation therapy is usually a function of dose and irradiated volume. The complex relationship between heart complication, dose and volume is illustrated in Fig. 1. Here the normal tissue complication probability (NTCP) for cardiac mortality from ischemic heart disease and myocardial infarction was taken from Gagliardi et al.  and modeled with the relative seriality model  resulting in the model parameters D50 = 52.3 Gy, γ = 1.28 and s = 1. The tangential breast treatments, which were analyzed by Darby et al. , were dosimetrically reproduced by assuming a simplified geometry of the heart. If it is assumed that the heart as a spherical volume which is irradiated partly with the treated dose D T = 60 Gy, then the % of irradiated heart volume V H is related to mean heart dose: V H = D mean /D T . The relative seriality model yields then NTCP as a function of mean heart dose:
Usually NTCP is given in absolute quantities. However, to be able to compare NTCP with the relative results of Darby et al. , we divide in Eq. 1 the original NTCP by R base = 0.38% , which is the base line risk for major coronary events and was taken from Table S11 of  for the age at exposure (59.6 years) and attained age (67.1 years) of the cohort analyzed by Darby et al. .
The ERR from Darby et al.  were then fitted to Eqs. 1 and 2. The data fits were produced with a least-squares algorithm with the software package PVWave (PV-Wave Advantage, PV-Wave Command Language, Version 9.01 - Numerics, Inc - 2008) using the epidemiological errors as weights. The model parameters D 50 , γ and s were obtained.
Fit of the Darby data to “a lot to a little”
The dose-response curve from the relative seriality model based on the model parameters of Gagliardi et al.  is shown as the solid line in Fig. 2 (D 50 = 52.3 Gy, γ = 1.28 and s = 1.0). This dose-response curve represents the situation “a lot to a little”, where the parts of the heart in the tangential treatment beams receive the treatment dose and the rest of the heart more or less no dose. The NTCP model shows qualitatively a linear dependence on mean heart dose as it was observed by Darby et al. . Quantitatively the model can be fitted to the Darby data which results in the dotted line with the model parameters D 50 = 36.5 Gy, γ = 1.29 and s = 0.75.
Model results for “a little to a lot”
We used the NTCP model based on the Gagliardi parameters (solid line) and the parameters obtained from the Darby data (dotted line) to calculate NTCP for a homogenous dose distribution as shown in Fig. 3. It is obvious that risk as a function of mean dose for a homogenously irradiated heart is not linear. In fact the dose-response curve exhibits the typical sigmoidal behavior with a negligible risk for major coronary events up to a dose of 15 Gy. For this reason the risk of 7.4% per Gy determined by Darby et al.  which is also shown in Fig. 3 as the symbols cannot be used in situations where the heart is irradiated “a little to a lot” and would strongly overestimate the risk of major cardiac events.
It must be noted that the epidemiological life span study of the A-bomb survivors shows a significant quadratic dose-response relationships for hypertension and for myocardial infarction among survivors exposed at less than 40 years of age . As the A-bomb survivors were irradiated with a more or less homogenous doses smaller than 3 Gy it corresponds to the scenario (“a little to a lot”). At a first glance these findings contradict our results illustrated in Fig. 3. However, the A-bomb survivors received a whole body radiation exposure and, since radiation dose is also associated with increased hypertension incidence, it can be suggested that radiation dose may be associated with chronic renal failure, thus explaining part of the mechanism for increased cardiovascular disease. Indeed Adams et al.  have shown that part of the risk of cardiovascular disease, particularly myocardial infarction risk, is mediated by renal dysfunction which is a result of the whole body irradiation (see Fig. 4). This relationship between renal dysfunction and cardiovascular disease is usually not present in radiotherapy, which is typically a local treatment.
It should also be noted that Ghobadi et al.  have shown that irradiation of heart, lung, or both independently induces specific cardiac dysfunction. These results show that treatment of the lung with radiation therapy can enhance cardiac toxicity through an unknown mechanism which is not yet included in bio-physical models.
It was shown in this work that epidemiological studies which predict a linear dose-response relationship for cardiovascular disease can be reproduced by bio-physical models for normal tissue complication. We have used the relative seriality model for a partly irradiated heart (“a lot to a little”) which models radiation therapy using two tangential fields. It was also shown that the model predicts a sigmoidal dose-response curve - and not a linear one - for cardiovascular disease in the situation of a homogenously irradiated heart (“a little to a lot”). This has direct consequences for the application of new treatment techniques, e.g. the use of the IMRT technique for treating the breast. IMRT is reducing the volume irradiated at high doses and may therefore result in more favorable cardiac outcomes than treating with tangents if the heart cannot be removed from the tangent by, e.g., deep inspiration breath hold techniques (DIBH). This is supported by a recent publication of Zagar et al.  who found in a prospective study that tangential breast treatment with DIBH for left-sided breast cancer is an effective means to avoid early RT-associated cardiac perfusion defects.
We believe that the results of epidemiological studies should not generally be used to predict normal tissue complications. It is better to use such data to optimize bio-physical models which can then be applied (with caution) to general treatment situations.
Darby SC, Ewertz M, Hall P. Ischemic heart disease after breast cancer radiotherapy. N Engl J Med. 2013;368(26):2527.
van den Bogaard VA, Ta BD, van der Schaaf A, Bouma AB, Middag AM, Bantema-Joppe EJ, van Dijk LV, van Dijk-Peters FB, Marteijn LA, de Bock GH, Burgerhof JG, Gietema JA, Langendijk JA, Maduro JH, Crijns AP. Validation and modification of a prediction model for acute cardiac events in patients with breast cancer treated with radiotherapy based on three-dimensional dose distributions to cardiac substructures. J Clin Oncol. 2017;35(11):1171–8.
Gagliardi G, Lax I, Ottolenghi A, et al. Long-term cardiac mortality after radiotherapy of breast cancer – application of the relative seriality model. Br J Radiol. 1996;69:839–46.
Källman P, Agren A, Brahme A. Tumour and normal tissue responses to fractionated non-uniform dose delivery. Int J Radiat Biol. 1992;62:249–62.
Yamada M, Wong FL, Fujiwara S, Akahoshi M, Suzuki G. Noncancer disease incidence in atomic bomb survivors, 1958-1998. Radiat Res. 2004;161(6):622–32.
Adams MJ, Grant EJ, Kodama K, Shimizu Y, Kasagi F, Suyama A, Sakata R, Akahoshi M. Radiation dose associated with renal failure mortality: a potential pathway to partially explain increased cardiovascular disease mortality observed after whole-body irradiation. Radiat Res. 2012;177(2):220–8.
Ghobadi G, van der Veen S, Bartelds B, de Boer RA, Dickinson MG, de Jong JR, Faber H, Niemantsverdriet M, Brandenburg S, Berger RM, Langendijk JA, Coppes RP, van Luijk P. Physiological interaction of heart and lung in thoracic irradiation. Int J Radiat Oncol Biol Phys. 2012;84(5):e639–46.
Zagar TM, Kaidar-Person O, Tang X, Jones EE, Matney J, Das SK, Green RL, Sheikh A, Khandani AH, McCartney WH, Oldan JD, Wong TZ, Marks LB. Utility of deep inspiration breath hold for left-sided breast radiation therapy in preventing early cardiac perfusion defects: a prospective study. Int J Radiat Oncol Biol Phys. 2017;97(5):903–9.
Availability of data and materials
All data are kept in the institute’s data reposity and are made available upon request.
US designed the study and wrote the manuscript. ME and MH helped with literature search and argumentation. All authors reviewed and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Schneider, U., Ernst, M. & Hartmann, M. The dose-response relationship for cardiovascular disease is not necessarily linear. Radiat Oncol 12, 74 (2017) doi:10.1186/s13014-017-0811-2
- Dose-response curves
- Cardiovascular disease
- Bio-physical modeling