# Is mean heart dose a relevant surrogate parameter of left ventricle and coronary arteries exposure during breast cancer radiotherapy A dosimetric evaluation based on individually-determined radiation dose (BACCARAT study)

S. Jacob, J. Camilleri, S. Derreumaux, V. Walker, O. Lairez, M. Lapeyre, E.

Bruguière, A. Pathak, M.-O. Bernier, D. Laurier, et al.

► To cite this version:

S. Jacob, J. Camilleri, S. Derreumaux, V. Walker, O. Lairez, et al. Is mean heart dose a relevant surrogate parameter of left ventricle and coronary arteries exposure during breast cancer radiotherapy A dosimetric evaluation based on individually-determined radiation dose (BACCARAT study). *Radiation Oncology, 2019, 14 (1), pp.29.* 10.1186/s13014-019-1234-z. *hal-02521605*

HAL Id: hal-02521605

https://hal.science/hal-02521605v1

Submitted on 27 Mar 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
RESEARCH

Open Access

<figure><img src="image_2.png" /></figure>

CrossMark

# Is mean heart dose a relevant surrogate parameter of left ventricle and coronary arteries exposure during breast cancer radiotherapy: a dosimetric evaluation based on individually-determined radiation dose (BACCARAT study)

Sophie Jacob¹*, Jérémy Camilleri², Sylvie Derreumaux³, Valentin Walker¹, Olivier Lairez⁴, Mathieu Lapeyre⁵, Eric Bruguière⁵, Atul Pathak⁶, Marie-Odile Bernier¹, Dominique Laurier⁷, Jean Ferrieres⁸,⁹, Olivier Gallocher², Igor Latorzeff², Baptiste Pinel², Denis Franck², Christian Chevelle², Gaëlle Jimenez² and David Broggio¹⁰

## Abstract

**Background:** Intra-individual heterogeneity of cardiac exposure is an issue in breast cancer (BC) radiotherapy that was poorly considered in previous cardiotoxicity studies mainly based on mean heart dose (MHD). This dosimetric study analyzes the distribution of individually-determined radiation doses to the heart and its substructures including coronary arteries and evaluate whether MHD is a relevant surrogate parameter of dose for these substructures.

**Methods:** Data were collected from the BACCARAT prospective study that included left or right unilateral BC patients treated with 3D-Conformal Radiotherapy (RT) between 2015 and 2017 and followed-up for 2 years with repeated cardiac imaging examinations. A coronary computed tomography angiography (CCTA) was performed before RT for all patients. Registration of the planning CT and CCTA images allowed delineation of the coronary arteries on the planning CT images. Using the 3D dose matrix generated during treatment planning and the added coronary contours, dose distributions were generated for whole heart and the following substructures: left ventricle (LV), left main coronary artery (LMCA), left anterior descending artery (LAD), left circumflex artery (LCX) and right coronary artery (RCA). A descriptive analysis of the physical doses in Gray (Gy) was performed, Dmean was the volume-weighted mean dose.

**Results:** Dose distributions were generated for 89 left-sided BC patients (MHD = 2.9 ± 1.5 Gy, Dmean\_LAD = 15.7 ± 3.1 Gy) and 15 right-sided BC patients (MHD = 0.5 ± 0.1 Gy; Dmean\_RCA = 1.2 ± 0.4 Gy). For left-sided BC patients, the ratio Dmean\_LAD/MHD was around 5. Pearson correlation coefficients between MHD and Dmean for delineated substructures were all statistically significant. However, for all substructures, the coefficient of determination R² indicated that the proportion of the variance in Dmean of the substructure predictable from MHD was moderate to low (in particular R² = 0.45 for LAD). Among left-sided BC patients with MHD < 3Gy, 56% of patients could nevertheless receive LAD doses above 40Gy (V40 > 0).

(Continued on next page)

<footnote>* Correspondence: sophie.jacob@irsn.fr
¹Institut de Radioprotection et de Sureté Nucléaire (IRSN), PSE-SANTE, SESANE, LEPID, BP17, 92262 Fontenay-aux-roses, Cedex, France
Full list of author information is available at the end of the article</footnote>
(Continued from previous page)

**Conclusion:** Our study illustrates that MHD is not enough to predict with confidence individual patient dose to the LV and coronary arteries, in particular the LAD. For precise radiotherapy-induced cardiotoxicity studies it would be necessary to consider the distribution of doses within these cardiac substructures rather than just the MHD.

**Trial registration:** ClinicalTrials.gov: NCT02605512, Registered 6 November 2015 - Retrospectively registered.

**Keywords:** Breast cancer, Radiotherapy, Heterogeneity in cardiac exposure, Mean heart dose, Coronary arteries doses, Cardiotoxicity

## Introduction

Radiotherapy (RT) for breast cancer is an essential part of adjuvant cancer treatment. RT reduces the risk of local recurrence and the risk of breast cancer mortality [1].

However, left-sided RT, especially, has been shown to induce excess cardiovascular mortality and morbidity [2–6]. The study by Darby et al. [5] found a linear relationship between the mean heart dose (MHD) and the rate of major coronary event, which increased of 7.4% per Gy of the MHD. These results on radiation-induced ischemic heart disease were confirmed in a more recent study of BC patients treated with three dimensional conformal radiation therapy (3D-CRT) [7]. The predominance of ischemic heart disease observed in these studies indicates that the coronary arteries may be the critical structures for the development of radiation induced heart morbidity. Among the three major coronary arteries (left anterior descending, circumflex, and right), the left anterior descending (LAD) supplies a major part of the myocardium. Therefore, occlusion of the LAD, may cause a large area of myocardial necrosis and lead to severe left ventricle impairment and congestive heart failure.

However, in many epidemiological studies on post radiotherapy cardiotoxicity, doses are typically described as those received by the entire heart with the MHD [5, 8]. Therefore, the MHD is used as the reference dose for analyzing dose-response relationship. But, dose distribution in the heart is not homogeneous, highest cardiac radiation doses can be observed in the apex and in the apical-anterior segment and some hot spots>50Gy persist in some parts of the heart [9]. Nevertheless, the anatomic distribution of RT-associated coronary artery diseases has been poorly studied whereas it can be supposed that the increased risk of coronary artery diseases would manifest largely in the coronary arteries that are directly within the radiation field. This was confirmed by a Sweden study [10] where patients with left-sided breast cancer RT had a statistically significant increased rate of stenosis in the LAD when compared with those with right-sided cancer. This observation was concordant with the location of typical RT fields and with the fact that the highest doses are likely to be delivered to the anterior heart, including the LAD. As a consequence, MHD may be a poor

surrogate parameter to reflect the dose to the cardiac sub-structures especially the LAD [11] and MHD may thus be a poorly relevant dose criterion for RT-induced cardiotoxicity studies and dose-response relationship.

However, a major issue is that coronary arteries are difficult structures to delineate because of their small volumes. Some atlas and auto-segmentation methodologies were developed for contouring cardiac substructures [12, 13], but they presented limits due to uncertainties for the smallest heart structures that are the coronary arteries. Another approach for coronary arteries’ radiation doses was presented in a previous work performed by Moignier et al. [14] where radiotherapy simulation CT scans and coronary computed tomography angiography (CCTA) of patients treated for a mediastinal Hodgkin lymphoma were used to merge thoracic and detailed cardiovascular anatomies and allowed personalized coronary arteries dose calculations.

In 2015, we launched the BACCARAT study, a cohort of a hundred of breast cancer patients treated with 3D-CRT and followed prospectively for 2 years with repeated cardiac imaging examinations including echocardiography and CCTA and measurement of circulating biomarkers. BACCARAT study aims to enhance knowledge on detection and prediction of early subclinical cardiac dysfunction and lesions induced by breast RT and on biological mechanisms potentially involved, based on functional and anatomical cardiac imaging combined with simultaneous assessment of multiple circulating biomarkers and accurate heart dosimetry.

With the rare opportunity given by the BACCARAT study to have simultaneous dosimetric and precise anatomical information from the simulation CT and CCTA images for all patients, the aim of this dosimetric study was to analyze the distribution of personalized individually-determined radiation doses to the heart and its substructures, in particular the coronary arteries, and evaluate whether MHD is a relevant surrogate parameter for the dose to cardiac substructures, in particular the LAD.

## Material and methods

The BACCARAT study, a monocentric prospective cohort study further detailed in a previous article [15],
consisted in the inclusion of 114 women treated with adjuvant 3D-CRT for left or right unilateral breast cancer in the Clinic Pasteur in Toulouse, France, without chemotherapy, and followed for 2 years after RT (ClinicalTrials.gov: NCT02605512). Inclusion period lasted from October 2015 to December 2017. Follow-up of patients is still ongoing, the end is foreseen for early 2020.

After the surgical treatment of breast cancer, all patients were treated with 3D-CRT with 6 and 25 MV photon beams by tangential fields. Patients underwent planning CT scan, not contrast enhanced, with CT slices every 1.25 mm.

Patients were planned for RT with or without irradiation of supraclavicular or internal mammary lymph nodes. The planning target volume dose was 50 Gy delivered in 5 weeks with 25 daily doses of 2 Gy or 47 Gy delivered in 5 weeks with 20 daily doses of 2.35Gy. Additional boost of 9–15 Gy could be applied to the tumor site with electron/photons beams with energies ranging from 6 MeV to 18 MeV. Patients were positioned on a breast board with both arms above the head. The treatment planning system (TPS) used to perform dose calculations was Eclipse™ with the Analytical Anisotropic Algorithm (AAA v13.6) (Varian Medical System, Palo Alto, CA, USA). Each patient's radiotherapy was planned such that the dose distribution was optimized and normalized to the International Commission on Radiation Units and Measurements (ICRU) reference point of the breast. Breath-hold gating was used for patients treated for left-sided breast cancer patients with heart very close to the anterior chest wall and for dose constraints achievement (Quantitative Analysis of Normal Tissue Effects in the Clinic - QUANTEC recommends keeping the volume of heart receiving at least 25 Gy (V25) less than 10% to keep the risk of cardiac mortality under 1% [16]; in addition, Clinique Pasteur attempts to keep MHD < 5 Gy). Dose-Volume-Histogram (DVH) for the heart was generated by the Clinic Pasteur radiotherapy department.

Before RT, a coronary computed tomography angiography (CCTA) was performed for all patients as planned in the BACCARAT protocol. One of the challenges of the BACCARAT study was to evaluate precisely absorbed doses to cardiac substructures, in particular focused on coronary arteries in order to enhance in a second step the dose-response analysis when considering early cardiotoxicity. The coronary arteries are narrow and moving structures; consequently the simulation CT scan makes proper contouring difficult or impossible. For dosimetric evaluation, the simulation CT scan, the CCTA, the RT dose and RT structure files in DICOM format were used as done in a previous study [14]. Merging anatomical information from the simulation CT scan and the CCTA was performed. The CCTA provided

an optimal coronary visualization but at a given moment of the cardiac cycle, generally the diastole. On the contrary, the simulation CT scan provided images integrating the breathing and heart beating motion. Consequently, a slight rotation of the heart around the cranio-caudal axis, a homothetic deformation along the three dimensions and a translation to match the coronary arteries origins were carried out for the registration. Once inserted in ISOgray TPS (version 4.2, Dosisoft, Cachan, France; http://www.dosisoft.com/en/radiotherapy/planning-products.html), manual delineation was performed for the left ventricle (LV), the left main coronary artery (LMCA), the left anterior descending artery (LAD), the left circumflex artery (LCX) and the right coronary artery (RCA). Using the 3D dose matrix generated during treatment planning and the new delineated substructures, DVH for LV and coronary arteries were generated with ISOgray TPS by the dosimetric department of IRSN in collaboration with the Clinic Pasteur radiotherapy department (Fig. 1).

From the DVHs, the following absorbed dose metrics for all delineated cardiac structures were calculated: Dmean (in Gy) is the volume-weighted mean dose; D2 (in Gy) is the minimal dose received by the most irradiated 2% of the structure volume, which can be considered as the near maximum dose; V5 (in %) is the relative volume exposed to at least 5 Gy (similar definition for V10, V25 etc.)

Descriptive analysis of the physical doses in Gray (Gy) was performed. Continuous variables are presented with mean, standard deviation, median and range values. Categorical values are presented with percentages. The chi-square test was used to compare categorical variables and Student's t-test or Wilcoxon non-parametric test if necessary was used to compare continuous variables. We defined individual ratio D2/Dmean as a kind of Homogeneity Index (HI) adapted to organ at risk as HI are usually used for planning target volume. Pearson's

<figure><img src="image_3.png" /><figcaption>**Fig. 1** CT dose-planned left tangential breast irradiation showing isodoses and delineated structures: Heart (orange circle), left ventricle (LV), Left anterior descending coronary artery (LAD) and right coronary artery (RCA)</figcaption></figure>
correlation coefficients were used for correlation analysis between MHD and mean doses to the different substructures. The relationship analysis between MHD and mean doses to the different substructures were further investigated based on linear regressions providing the R² value which corresponded to the coefficient of determination indicating the proportion of the variance in Dmean of the substructure predictable from MHD. R² < 0.70 was considered insufficient for prediction and surrogate parameter purpose. P < 0.05 was considered statistically significant. All statistical analysis was performed with SAS software V9.2.

## Results

Retrospective dosimetric evaluation for all cardiac substructures was available for 104 BC patients (89 left-sided and 15 right-sided) from the BACCARAT study (Table 1). Mean age of the population was 58 years old, with no significant difference between right-sided and left-sided BC patients. Most patients (85/104) were diagnosed with an

<table><caption>**Table 1** Patients baseline characteristics</caption><thead><tr><th></th><th>Left-sided<br>BC patients<br>N = 89</th><th>Right-sided<br>BC patients<br>N = 15</th></tr></thead><tbody><tr><td>Age (in years)</td><td>56.0 ± 8.34<br>(40–74)</td><td>58.3 ± 5.8<br>(52–67)</td></tr><tr><td>Type of cancer</td><td></td><td></td></tr><tr><td>In situ</td><td>16 (18%)</td><td>3 (20%)</td></tr><tr><td>Invasive</td><td>73 (82%)</td><td>12 (80%)</td></tr><tr><td>Tumor size (in mm)</td><td>12.7<br>(4–100)</td><td>8.3<br>(4–13)</td></tr><tr><td>Grade</td><td></td><td></td></tr><tr><td>1</td><td>37 (41%)</td><td>7 (47%)</td></tr><tr><td>2</td><td>42 (47%)</td><td>7 (47%)</td></tr><tr><td>3</td><td>10 (11%)</td><td>1 (6%)</td></tr><tr><td>Surgery</td><td></td><td></td></tr><tr><td>Breast conserving</td><td>82 (92%)</td><td>15 (100%)</td></tr><tr><td>Mastectomy</td><td>7 (8%)</td><td>0 (0%)</td></tr><tr><td>Prescribed dose</td><td></td><td></td></tr><tr><td>50 Gy</td><td>68 (76%)</td><td>13 (87%)</td></tr><tr><td>47 Gy</td><td>21 (24%)</td><td>2 (13%)</td></tr><tr><td>Boost</td><td>81 (91%)</td><td>14 (93%)</td></tr><tr><td>9Gy</td><td>0 (0%)</td><td>1 (6%)</td></tr><tr><td>12Gy</td><td>66 (73%)</td><td>4 (27%)</td></tr><tr><td>12.5 Gy</td><td>13 (14%)</td><td>9 (60%)</td></tr><tr><td>15Gy</td><td>2 (2%)</td><td>0</td></tr><tr><td>Regional lymph nodes irradiation</td><td></td><td></td></tr><tr><td>Supraclavicular alone</td><td>1 (1%)</td><td>0</td></tr><tr><td>Internal mammary alone</td><td>3 (3%)</td><td>0</td></tr><tr><td>Both</td><td>27 (30%)</td><td>1 (7%)</td></tr></tbody></table>

invasive ductal carcinoma and underwent breast conserving surgery. The prescribed radiation dose was 50 Gy in 25 sessions for more than three forth of the population. An additional boost of 9–15 Gy was applied to the site of the tumor if clinically indicated. Regional lymph nodes irradiation was performed in 31 left-sided BC patients and 1 right-sided BC patient.

Results for each dose metric for the heart and its sub-structures are shown in Table 2. The mean MHD was 2.95 ± 1.49 Gy for left-sided RT, and 0.46 ± 0.12 Gy for right-sided RT. All patient met the dose constraint of the heart (V25Gy < 10%) with a maximum V25 value of 8.7%. The inter-patient variability in MHD was important especially for left-sided RT with a range of values from 0.87 to 6.72 Gy, further confirmed by the range in near-maximum dose D2 from 3.95 to 48.87 Gy. The intra-patient variability for heart doses was also extremely high, with average individual ratio D2/MHD of 8.4 for left-sided patient and 4.6 for right-sided BC patients illustrating the heterogeneity of doses within whole heart structure. Considering the other cardiac structures, LV and LAD were the most exposed structures with respectively mean doses of 6.2 Gy and 15.7 Gy (as illustrated for a patient in Fig. 1), but intra-patient variability in doses (D2/Dmean) was lower than observed for the whole heart (5.7 and 2.6 respectively). For the RCA, we observed that the mean absorbed dose was higher for right irradiation than left one (0.74 ± 0.53 Gy vs. 1.25 ± 0.51 Gy), and higher than MHD for right-sided BC.

The ratios of Dmean for delineated structure and MHD are presented in Table 3. It allowed observing that for left-sided BC patients, the mean ratio Dmean LAD/MHD was around 5 and around 2 for Dmean LV/MHD. All other ratios were below 1 except for RCA in right-sided BC patients (2.7 ± 0.7). The strongest correlation with MHD was observed for LV (r = 0.78) and LAD (r = 0.67) for left-sided patients, and for LMCA and LAD (r = 0.81 and 0.71 respectively) for right-sided BC patients. Linear regression between Dmean for delineated structure and MHD are presented in Figs. 2 and 3. For every increase in MHD among left-sided patients, the Dmean LV and Dmean LAD increased on average 1.9 Gy and 3.7Gy respectively. However, the coefficients of determination R² values indicated that the proportion of the variance in Dmean LV or LAD predictable from MHD was moderate for LV (R² = 0.63) to low for LAD (R² = 0.45). For every increase in MHD among right-sided BC patients, the Dmean RCA increased on average 2.0 Gy, but again R² remained low (R² = 0.36). For all others structures the relationship was moderate to weak. With no R² value above 0.70, the predictive value of the MHD was not good for cardiac substructures including coronary arteries.
<table><caption>**Table 2** Dosimetric parameters for heart and cardiac substructures including coronary arteries</caption><thead><tr><th rowspan="2"></th><th colspan="2">Left-sided BC patients<br/>N = 89</th><th colspan="2">Right-sided BC patients<br/>N = 15</th></tr><tr><th>Mean ± SD</th><th>Range</th><th>Mean ± SD</th><th>Range</th></tr></thead><tbody><tr><td><strong>Heart</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>MHD (Gy)</td><td>2.95 ± 1.49</td><td>0.87–6.72</td><td>0.46 ± 0.12</td><td>0.25–0.67</td></tr><tr><td>D2 (Gy)</td><td>26.80 ± 17.40</td><td>3.95–48.87</td><td>2.16 ± 0.90</td><td>1.11–3.57</td></tr><tr><td>Individual ratio D2/MHD</td><td>8.4 ± 3.6</td><td>2.9–15.4</td><td>4.6 ± 0.6</td><td>3.6–5.8</td></tr><tr><td><strong>Left Ventricle</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Dmean (Gy)</td><td>6.25 ± 3.53</td><td>1.16–13.42</td><td>0.09 ± 0.04</td><td>0.04–0.17</td></tr><tr><td>D2 (Gy)</td><td>34.01 ± 16.93</td><td>4.49–55.48</td><td>0.34 ± 0.13</td><td>0.17–0.54</td></tr><tr><td>Individual ratio D2/Dmean</td><td>5.7 ± 2.0</td><td>2.6–11.0</td><td>3.9 ± 0.7</td><td>2.6–5.3</td></tr><tr><td><strong>Left Anterior Descending Artery</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Dmean (Gy)</td><td>15.68 ± 8.13</td><td>1.68–37.61</td><td>0.37 ± 0.15</td><td>0.13–0.54</td></tr><tr><td>D2 (Gy)</td><td>37.97 ± 15.93</td><td>2.19–56.42</td><td>0.45 ± 0.17</td><td>0.19–0.61</td></tr><tr><td>Individual ratio D2/Dmean</td><td>2.6 ± 0.9</td><td>1.3–6.0</td><td>2.3 ± 0.6</td><td>1.3–3.2</td></tr><tr><td><strong>Left Circumflex Artery</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Dmean (Gy)</td><td>1.61 ± 0.83</td><td>0.53–4.34</td><td>0.14 ± 0.08</td><td>0.04–0.25</td></tr><tr><td>D2 (Gy)</td><td>1.87 (median)</td><td>0.75–48.66</td><td>0.27 ± 0.13</td><td>0.09–0.45</td></tr><tr><td>Individual ratio D2/Dmean</td><td>1.3 (median)</td><td>0.8–43.4</td><td>2.2 ± 0.8</td><td>1.5–4.0</td></tr><tr><td><strong>Left Main Coronary Artery</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Dmean (Gy)</td><td>1.28 ± 0.69</td><td>0.60–4.72</td><td>0.11 ± 0.05</td><td>0.04–0.21</td></tr><tr><td>D2 (Gy)</td><td>1.39 (median)</td><td>0.69–49.45</td><td>0.25 ± 0.13</td><td>0.06–0.46</td></tr><tr><td>Individual ratio D2/Dmean</td><td>1.2 (median)</td><td>1.0–35.6</td><td>1.2 ± 0.2</td><td>0.8–1.6</td></tr><tr><td><strong>Right Coronary Artery</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Dmean (Gy)</td><td>0.74 ± 0.53</td><td>0.14–4.39</td><td>1.25 ± 0.41</td><td>0.68–2.01</td></tr><tr><td>D2 (Gy)</td><td>1.38 ± 1.04</td><td>0.35–8.61</td><td>1.77 ± 0.74</td><td>0.33–3.54</td></tr><tr><td>Individual ratio D2/Dmean</td><td>1.9 ± 0.4</td><td>1.1–2.9</td><td>1.4 ± 0.4</td><td>0.3–2.2</td></tr></tbody></table>

BC breast cancer, SD standard deviation, MHD Mean heart dose, D2 minimal dose received by the most irradiated 2% of the structure volume, Dmean mean dose to the structure

Considering the left sided BC RT patients, even in the “low exposed” group (MHD < 3Gy) high exposure to the LV and LAD was observed, D2 value that could reach 55.48Gy and 56.42Gy, higher than in the group of patient with MHD > 3Gy (Table 4). Moreover, in the “low exposed” group, 56% of the population could receive doses above 40 Gy to the LAD (75% in the high exposed group) that could reach 51.6% of the LAD volume, and 3 patients even received doses above 50Gy to the LAD with volumes of LAD of 1.2, 5.3 and 24.1% for these 3 patients.

Discussion

Although there was considerable decrease in doses to the heart over past few years [8, 17, 18], radiation induced heart disease is still a concern due to the improvement in breast cancer patients survival. Mean heart dose was often used as the reference measure for cardiotoxicity studies [5]. However, there is increasing

evidence toward the importance to consider individual cardiac substructures as some studies have implicated the left anterior descending artery [10, 19], as well as the left ventricle [7] as important components of the heart associated with radiation induced heart disease.

Our study highlighted that MHD was not representative of order of magnitude of mean dose to the left ventricle or even more to the most exposed coronary, i.e. LAD. The mean LAD dose (15.6Gy) observed among 89 left-sided BC patients was substantially higher than the mean heart dose (2.9 Gy) as previously observed [20]. Inter-individual variability in LAD exposure dose, illustrated by the wide range of LAD dose from 1.68Gy to 37.60 Gy, could be explained by several factors: - variation in heart size (larger heart tend to push the LAD into the radiation field), - variation in breast size (it modifies the extend of the tangential field), - individual coronary topology (coronaries are more or less tortuous [9], and in this study some were delineated longer than others because of
<table><caption>**Table 3** Association parameters between mean dose to the left ventricle and coronary arteries and the mean heart dose</caption><thead><tr><th></th><th>LV</th><th>LAD</th><th>LCX</th><th>LMCA</th><th>RCA</th></tr></thead><tbody><tr><td colspan="6">Left-sided BC patients</td></tr><tr><td>Ratio Dmean struc/MHD (mean ± SD)</td><td>2.1 ± 0.6</td><td>5.6 ± 2.2</td><td>0.7 ± 0.4</td><td>0.5 ± 0.2</td><td>0.3 ± 0.1</td></tr><tr><td>Correlation with MHD (r, p-value)</td><td>0.78, &lt; 0.0001</td><td>0.67, &lt; 0.0001</td><td>0.33, 0.017</td><td>0.47, &lt; 0.0001</td><td>0.60, &lt; 0.0001</td></tr><tr><td colspan="6">Right-sided BC patients</td></tr><tr><td>Ratio Dmean struc/MHD (mean ± SD)</td><td>0.2 ± 0.1</td><td>0.2 ± 0.1</td><td>0.3 ± 0.1</td><td>0.8 ± 0.2</td><td>2.7 ± 0.7</td></tr><tr><td>Correlation with MHD (r, p-value)</td><td>0.55, 0.03</td><td>0.71, 0.0026</td><td>0.64, 0.01</td><td>0.81, 0.0002</td><td>0.60, 0.02</td></tr></tbody></table>

LV Left Ventricle, LAD Left Anterior Descending Artery, LCX Left Circumflex Artery, LMCA Left Main Coronary Artery, RCA Right Coronary Artery, BC breast cancer, SD standard deviation, MHD Mean heart dose, Dmean struc mean dose absorbed by the delineated substructure

intrinsic limitation of the CCTA imaging), – presence or not of boost and regional lymph node irradiation. The difference between LAD dose and MHD is striking, but it can be explained by the fact that the LAD is located in the anterior region of the heart where it is exposed to the tangential fields used in left sided breast RT. Moreover, we observed a mean ratio Dmean LAD/MHD around 5, confirming prior findings [21]. However, the association between the MHD and delineated substructures, in particular the LAD mean dose, was

not precise enough to predict with confidence individual patient dose as shown in Fig. 2 and previously observed [22]. Despite significant correlation as observed in Table 3, this analysis showed that the predictive value of the MHD was not good for cardiac substructures including coronary arteries, even for LAD. More than 55% of patients with mean heart dose<3Gy could still be exposed to 40Gy or more to a part of LAD as illustrated in Table 4. In Darby’s study [5], estimated mean doses to the LAD and MHD were also

<figure><img src="image_4.png" /><figcaption>**Fig. 2** Relationship between MHD and mean doses to left ventricle and coronary arteries for left-sided BC patients</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Fig. 3** Relationship between MHD and mean doses to left ventricle and coronary arteries for right-sided BC patients</figcaption></figure>

correlated (correlation coefficient, 0.76 vs. 0.67 in our study) using a methodology based on a computed tomography of a patient with typical anatomy due to the period of inclusion of patients in a pre-computed tomography era. Such non individually-determined method could have led to inaccurate LAD dose estimations due to difficulties in localizing this structure in contrast with our method based on individual CT and CCTA. This could explain that the MHD was a better predictor of the rate of major coronary events than the mean dose to the LAD and even more surprising that inclusion of the mean LAD dose to MHD failed to improve the prediction of major coronary events rates as coronary arteries, in particular LAD, are critical to the subsequent development of radiation-induced cardiac diseases. In a more recent study [7], dose to the left ventricle (V5 precisely) was a better predictor of major coronary events than MHD. Similarly, in our study, the left ventricle was chosen as a structure of interest since its dose may be more representative of cardiac endpoints than the MHD. Indeed, left

ventricle pumps blood into the systemic circulation to most of the body through the aorta (while the right ventricle fills only the lungs). The most common reason for referral to echocardiography is left ventricular function. Assessment of left ventricular function is extremely important as it correlates with symptoms, prognosis, events, and complications in a large number of conditions and many decisions in cardiology are based on left ventricular function [23]. From a dosimetry point of view, the association between the MHD and the LV dose as illustrated in Fig. 1 showed that MHD also failed in predicting with confidence individual LV dose.

In the cardiotoxicity assessment of breast cancer RT, dosimetric evaluation for RCA was poorly studied. However, RCA, is the second largest artery, after left main coronary artery, supplying the heart. During breast cancer treatment of right breast the proximal part of RCA can be included in the irradiation field. In our study we observed among right sided BC patients that mean heart dose was 0.46Gy but mean RCA dose was 1.25Gy. These
<table><caption>**Table 4** Highest doses to Left Ventricle and Left Anterior Descending Artery according to low or high MHD among left-sided BC patients</caption><thead><tr><th>Left-sided BC patients</th><th>MHD ≤ 3Gy<br>N = 45</th><th>MHD > 3Gy<br>N = 44</th></tr></thead><tbody><tr><td colspan="3"><strong>Left Ventricle</strong></td></tr><tr><td>D2 (Gy): mean (range)</td><td>25.75 (4.49–55.48)</td><td>42.25 (21.30–53.17)</td></tr><tr><td>patients with V40 > 0</td><td>30 (67%)</td><td>34 (77%)</td></tr><tr><td>maximum value of V40</td><td>9.6%</td><td>20.0%</td></tr><tr><td>patients with V50 > 0</td><td>3 (7%)</td><td>14 (32%)</td></tr><tr><td>maximum value of V50</td><td>6.6%</td><td>7.0%</td></tr><tr><td colspan="3"><strong>Left Anterior Descending Artery</strong></td></tr><tr><td>D2 (Gy): mean (range)</td><td>31.99 (2.19–56.42)</td><td>44.08 (18.60–53.32)</td></tr><tr><td>patients with V40 > 0</td><td>25 (56%)</td><td>33 (75%)</td></tr><tr><td>maximum value of V40</td><td>51.6%</td><td>59.5%</td></tr><tr><td>patients with V50 > 0</td><td>3 (7%)</td><td>8 (18%)</td></tr><tr><td>maximum value of V50</td><td>24.1%</td><td>50.7%</td></tr></tbody></table>

MHD mean heart dose, D2 minimal dose received by the most irradiated 2% of the structure volume, V40 (in %) is the relative volume exposed to at least 40 Gy

results were concordant with those previously published [24] with a mean dose to the RCA of 1.92 Gy. The mean ratio Dmean RCA/Dmean Heart was around 3 and to a certain extend the RCA of right irradiated patients may play the same role than LAD for left irradiated patients.

Because of their small volumes and inherent contouring variability, the coronary arteries (in particular the LAD) are difficult structures to delineate. Dreaming of a 'Swiss Army Knife' for cardiac doses, some studies developed methodology based on atlas and auto-segmentation that could help to standardize contouring for these structures [12, 13, 25, 26]. These algorithms were able to generate valid segmentations for major cardiac structures (whole heart, Left ventricle, Right ventricle, ...) but not for the small structures such as coronary arteries. In our study, rather than considering a generic coronary model based on anatomic charts or an atlas, we had the opportunity to perform an individually determined cardiac dosimetry including small cardiac substructures based on registration of the planning CT and CCTA images. However, a limitation of our study is that doses calculated by the TPS might not be reliable enough (out-of-field, lung/heart interface) but currently, it is what is available in the usual clinical context. Further investigations based on our data taking into account some of these dosimetric parameters based on heart DVH parameters and possibly other available clinical parameters, might lead to a better definition of the coronary dose, even if not delineated.

Very few previous studies have clearly defined the anatomic distribution of RT-associated coronary artery disease [14, 19]. One would reasonably hypothesize that

the increased risk of coronary artery disease would be dose dependent and would manifest largely in the coronary arteries that are directly within the radiation field [27] as previously observed [19]. Our study allowed observing that the value of MHD could not help so much in understanding the dose effect relationship for the LAD. More research is needed to determine which indicator of heart dose from breast RT is the most significant for determining cardiac toxicity and morbidity. Doses to the heart has to be as low as possible. Actually, the QUANTEC - Quantitative Analysis of Normal Tissue Effects in the Clinic- recommendations, specify that the heart should always be contoured and that V25Gy < 10% (in 2 Gy per fraction) [16]. No dose recommendation exists for the coronary arteries, except one Danish study proposed a threshold with a maximum of 20 Gy to be given to the LAD [28]. With MHD as low as possible, the coronary arteries exposure could be as low as possible ad gating is an option to decrease the doses [22, 29] but we could not test it in our study as only 4 patients were concerned.

## Conclusion

To our knowledge it is the first study to use patient CCTA combined with patient-specific simulation CT scan to estimate dose to the substructures of the heart, including coronary arteries, in such a large number of patients. It allowed observing that taking into account the MHD only could nevertheless lead to excessive cardiac substructure irradiation and the predictive value of the MHD was not good for cardiac substructures including coronary arteries. The results of the study indicate that for precise radiotherapy-induced cardiotoxicity studies, it would be necessary to assess the dose delivered to the whole heart as well as to the cardiac substructures, in particular the LAD. For our BACCARAT study, this accurate heart dosimetry should be fruitful in the analysis of precise myocardial dysfunction based on cardiac echography and for precise analysis of coronary artery segments based on CCTA.

### Abbreviations

3D-CRT: 3 dimensional conformal radiotherapy; BACCARAT: Breast Cancer and Cardiotoxicity Induced by Radiotherapy; BC: Breast Cancer; CCTA: Coronary computed tomography angiography; CT: Computed tomography; DVH: Dose Volume Histogram; LAD: Left anterior descending coronary artery; LCX: Left circumflex artery; LMCA: Left main coronary artery; LV: Left ventricle; MHD: Mean heart dose; RCA: Right coronary artery; RT: Radiotherapy; TPS: Treatment Planning system

### Funding

The study received funding for collection of data from the Fédération Française de Cardiologie (FFC), from Electricité de France (EDF) and from the H2020 Euratom research and training programme 2014–2018 under grant agreement No 755523 in the frame of the MEDIRAD project.

### Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
**Authors' contributions**

Conceived study: SJ, JC, SD, ML, MOB, JF, GJ and DB. Data analysis: JC, VW, SD, GJ and DB, CM. Statistical analysis: SJ, DB and VW. All authors participated in the writing, and manuscript edition. All authors read and approved the final manuscript.

**Ethics approval and consent to participate**

This study received ethical approval from the French South West Committee for Protection of Persons (ID: CPP2015/66/2015-A00990-69) and from National Agency for Medical and Health products Safety (Reference: 150873B-12). All patients enrolled in the study provided their written informed consent.

**Consent for publication**

Consent for publication of the individual person's image was obtained from that person.

**Competing interests**

The authors declare that they have no competing interests.

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Author details**

¹Institut de Radioprotection et de Sureté Nucléaire (IRSN), PSE-SANTE, SESANE, LEPID, BP17, 92262 Fontenay-aux-roses, Cedex, France. ²Clinique Pasteur, Radiothérapie (Oncorad), Toulouse, France. ³Institut de Radioprotection et de Sureté Nucléaire (IRSN), PSE-SANTE, SER, UEM, Fontenay-aux-Roses, France. ⁴University Hospital Rangueil, Cardiac Imaging Center, Toulouse, France. ⁵Clinique Pasteur, Radiologie, Toulouse, France. ⁶Clinique Pasteur, Unité d'Hypertension Artérielle, Facteurs de Risque et Insuffisance Cardiaque, Toulouse, France. ⁷Institut de Radioprotection et de Sureté Nucléaire (IRSN), PSE-SANTE, SESANE, Fontenay-aux-Roses, France. ⁸Department of Cardiology, Toulouse-Rangueil University Hospital, Toulouse, France. ⁹INSERM, University Paul Sabatier, UMR1027, Epidemiology of Cardiovascular Diseases, Toulouse, France. ¹⁰Institut de Radioprotection et de Sureté Nucléaire (IRSN), PRP-HOM/SDOS/LEDI, Fontenay-aux-Roses, France.

Received: 4 December 2018 Accepted: 28 January 2019

Published online: 07 February 2019

**References**

1. Darby S, McGale P, Correa C, Taylor C, et al. Effect of radiotherapy after breast-conserving surgery on 10-year recurrence and 15-year breast cancer death: meta-analysis of individual patient data for 10,801 women in 17 randomised trials. Lancet. 2011;378(9804):1707-16. https://doi.org/10.1016/S0140-6736(11)61629-2.

2. Clarke M, Collins R, Darby S, Davies C, et al. Effects of radiotherapy and of differences in the extent of surgery for early breast cancer on local recurrence and 15-year survival: an overview of the randomised trials. Lancet. 2005;366(9503):2087-106. https://doi.org/10.1016/S0140-6736(05)67887-7.

3. Darby SC, McGale P, Taylor CW, Peto R. Long-term mortality from heart disease and lung cancer after radiotherapy for early breast cancer: prospective cohort study of about 300,000 women in US SEER cancer registries. Lancet Oncol. 2005;6(8):557-65. https://doi.org/10.1016/S1470-2045(05)70251-5.

4. Giordano SH, Kuo YF, Freeman JL, Buchholz TA, et al. Risk of cardiac death after adjuvant radiotherapy for breast cancer. J Natl Cancer Inst. 2005;97(6): 419-24. https://doi.org/10.1093/jnci/dji067.

5. Darby SC, Ewertz M, McGale P, Bennet AM, et al. Risk of ischemic heart disease in women after radiotherapy for breast cancer. N Engl J Med. 2013; 368(11):987-98. https://doi.org/10.1056/NEJMoa1209825.

6. Jacob S, Ferrières J. Breast cancer radiotherapy: a case of double jeopardy. Arch Cardiovasc Dis. 2016;109(11):587-90. https://doi.org/10.1016/j.acvd. 2016.09.001.

7. van den Bogaard VA, Ta BD, van der Schaaf A, Bouma AB, et al. 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 Off J Am Soc Clin Oncol. 2017;35(11):1171-8. https://doi.org/10.1200/JCO.2016.69.8480.

8. Taylor CW, Nisbet A, McGale P, Darby SC. Cardiac exposures in breast cancer radiotherapy: 1950s-1990s. Int J Radiat Oncol Biol Phys. 2007;69(5):1484-95. https://doi.org/10.1016/j.ijrobp.2007.05.034.

9. Moignier A, Broggio D, Derreumaux S, El Baf F, et al. Dependence of coronary 3-dimensional dose maps on coronary topologies and beam set in breast radiation therapy: a study based on CT angiographies. Int J Radiat Oncol Biol Phys. 2014;89(1):182-90. https://doi.org/10.1016/j.ijrobp.2014.01.055.

10. Correa CR, Litt HI, Hwang WT, Ferrari VA, et al. Coronary artery findings after left-sided compared with right-sided radiation treatment for early-stage breast cancer. J Clin Oncol Off J Am Soc Clin Oncol. 2007;25(21):3031-7. https://doi.org/10.1200/JCO.2006.08.6595.

11. Duma MN, Herr AC, Borm KJ, Trott KR, et al. Tangential field radiotherapy for breast Cancer-the dose to the heart and heart subvolumes: what structures must be contoured in future clinical trials? Front Oncol. 2017;7:130. https:// doi.org/10.3389/fonc.2017.00130.

12. Feng M, Moran JM, Koelling T, Chughtai A, et al. Development and validation of a heart atlas to study cardiac exposure to radiation following treatment for breast cancer. Int J Radiat Oncol Biol Phys. 2011;79(1):10-8. https://doi.org/10.1016/j.ijrobp.2009.10.058.

13. Van Dijk-Peters FBJ, Sijtsema NM, Kierkels RGJ, Vliegenthart R, et al. Validation of a multi-atlas based auto-segmentation of the heart in breast cancer patients. Radiother Oncol. 2015;115(Supp 1):S132-S3.

14. Moignier A, Broggio D, Derreumaux S, Beaudre A, et al. Coronary stenosis risk analysis following Hodgkin lymphoma radiotherapy: a study based on patient specific artery segments dose calculation. Radiother Oncol. 2015; 117(3):467-72. https://doi.org/10.1016/j.radcon.2015.07.043.

15. Jacob S, Pathak A, Franck D, Latorzeff I, et al. Early detection and prediction of cardiotoxicity after radiation therapy for breast cancer: the BACCARAT prospective cohort study. Radiat Oncol. 2016;11:54. https://doi.org/10.1186/ s13014-016-0627-5.

16. Gagliardi G, Constine LS, Moiseenko V, Correa C, et al. Radiation dose-volume effects in the heart. Int J Radiat Oncol Biol Phys. 2010;76(3 Suppl): S77-85. https://doi.org/10.1016/j.ijrobp.2009.04.093.

17. Taylor CW, Povall JM, McGale P, Nisbet A, et al. Cardiac dose from tangential breast cancer radiotherapy in the year 2006. Int J Radiat Oncol Biol Phys. 2008;72(2):501-7. https://doi.org/10.1016/j.ijrobp.2007.12.058.

18. Taylor AM, Wang Z, Macaulay E, Jagsi R, et al. Exposure of the heart in breast cancer radiotherapy 2003-2013: a systematic review. Int J Radiat Oncol, Biol, Phys. 2015;93(4):845-53.

19. Nilsson G, Holmberg L, Garmo H, Duvernoy O, et al. Distribution of coronary artery stenosis after radiation for breast cancer. J Clin Oncol Off J Am Soc Clin Oncol. 2012;30(4):380-6. https://doi.org/10.1200/JCO.2011.34.5900.

20. Drost L, Yee C, Lam H, Zhang L, et al. A systematic review of heart dose in breast radiotherapy. Clin Breast Cancer. 2018;18(5):e819-e92. https://doi.org/ 10.1016/j.clcancer.2018.05.010.

21. Moignier A, Derreumaux S, Broggio D, Beurrier J, et al. Potential of hybrid computational phantoms for retrospective heart dosimetry after breast radiation therapy: a feasibility study. Int J Radiat Oncol Biol Phys. 2013;85(2): 492-9. https://doi.org/10.1016/j.ijrobp.2012.03.037.

22. Becker-Schiebe M, Stockhammer M, Hoffmann W, Wetzel F, et al. Does mean heart dose sufficiently reflect coronary artery exposure in left-sided breast cancer radiotherapy? : Influence of respiratory gating. Strahlenther Onkol. 2016;192(9):624-31. https://doi.org/10.1007/s00066-016-1011-y.

23. Nagueh SF, Smiseth OA, Appleton CP, Byrd BF 3rd, et al. Recommendations for the evaluation of left ventricular diastolic function by echocardiography: an update from the American Society of Echocardiography and the European Association of Cardiovascular Imaging. Eur Heart J Cardiovasc Imaging. 2016;17(12):1321-60. https://doi.org/10.1093/ehjci/jew082.

24. Altinok A, Askeroglu O, Doyuran M, Caglar M, et al. Dosimetric evaluation of right coronary artery in radiotherapy for breast cancer. Med Dosim. 2018. https://doi.org/10.1016/j.meddos.2018.06.006.

25. Duane F, Aznar MC, Bartlett F, Cutter DJ, et al. A cardiac contouring atlas for radiotherapy. Radiother Oncol. 2017;122(3):416-22. https://doi.org/10.1016/j. radcon.2017.01.008.

26. Kaderka R, Gillespie EF, Mundt RC, Bryant AK, et al. Geometric and dosimetric evaluation of atlas based auto-segmentation of cardiac structures in breast cancer patients. Radiother Oncol. 2018. https://doi.org/ 10.1016/j.radcon.2018.07.013.

27. Zagar TM, Marks LB. Breast cancer radiotherapy and coronary artery stenosis: location, location, location. J Clin Oncol Off J Am Soc Clin Oncol. 2012;30(4): 350-2. https://doi.org/10.1200/JCO.2011.38.9304.
28. Offersen B, Hojris I, Overgaard M. Radiation-induced heart morbidity after adjuvant radiotherapy of early breast cancer - is it still an issue? Radiother Oncol. 2011;100(2):157-9. https://doi.org/10.1016/j.radonc.2011.08.007.

29. Edvardsson A, Nilsson MP, Amptoulach S, Ceberg S. Comparison of doses and NTCP to risk organs with enhanced inspiration gating and free breathing for left-sided breast cancer radiotherapy using the AAA algorithm. Radiat Oncol. 2015;10:84. https://doi.org/10.1186/s13014-015-0375-y.

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

<figure><img src="image_6.png" /></figure>