A clinical trial to compare a 3D‐printed bolus with a conventional bolus with the aim of reducing cardiopulmonary exposure in postmastectomy patients with volumetric modulated arc therapy

Abstract Background We compared the dosimetry, application, and acute toxicity of a 3D‐printed and a conventional bolus for postmastectomy radiotherapy (PMRT) with volumetric modulated arc therapy (VMAT). Materials and Methods Eligible patients (n = 75) with PMRT breast cancer were randomly selected to receive VMAT with a conventional bolus or a 3D‐printed bolus. The primary endpoint was a 10% decrease in the mean heart dose to left‐sided breast cancer patients. The secondary endpoint was a 5% decrease in the mean ipsilateral lung dose to all patients. A comparative analysis was carried out of the dosimetry, normal tissue complication probability (NTCP), acute skin toxicity, and radiation pneumonitis. Results Compared to a conventional bolus, the mean heart dose in left‐sided breast cancer was reduced by an average of 0.8 Gy (5.5 ± 1.3 Gy vs. 4.7 ± 0.8 Gy, p = 0.035) and the mean dose to the ipsilateral lung was also reduced by an average of 0.8 Gy (12.4 ± 1.0 Gy vs. 11.6 ± 0.8 Gy, p < 0.001). The values for V50Gy of the PTV of the chest wall for the 3D‐printed and conventional boluses were 95.4 ± 0.6% and 94.8 ± 0.8% (p = 0.026) and the values for the CI of the entire PTV were 0.83 ± 0.02 and 0.80 ± 0.03 (p < 0.001), respectively. The NTCP for the 3D‐printed bolus was also reduced to an average of 0.14% (0.32 ± 0.19% vs. 0.18 ± 0.11%, p = 0.017) for the heart and 0.45% (3.70 ± 0.67% vs. 3.25 ± 0.18%, p < 0.001) for the ipsilateral lung. Grade 2 and Grade 1 radiation pneumonitis were 0.0% versus 7.5% and 14.3% versus 20.0%, respectively (p = 0.184). Conclusions The 3D‐printed bolus may reduce cardiopulmonary exposure in postmastectomy patients with volumetric modulated arc therapy.


| INTRODUCTION
With advances in treatment technology, survival rates for breast cancer patients are increasing. 1 Cardiovascular disease has become the leading cause of non-tumor-related death in breast cancer patients with long-term survival, 2 especially when anthracycline, trastuzumab, and radiotherapy are involved. 3,4 Protection of cardiopulmonary function is therefore treated as paramount by radiation oncologists. Although there is no clear upper limit for a safe cardiac exposure dose, the mean heart dose (MHD) is an important practical parameter that has been used to predict radiotherapy-induced heart disease (RIHD) in past decades. [4][5][6][7] Modified radical surgery continues to play a fundamental role for breast cancer patients in China, and reaching up to 80% of breast cancer patients undergo this surgery in some regions. Radiotherapy for breast cancer after mastectomy contributes to a large degree to a reduction in the risk of local relapse. [6][7][8] The megavoltage photon beam creates a dose build-up effect, 9,10 and to ensure that the dose distribution conforms to the target volume, a bolus is frequently used on the patient's skin. 11 Recently, several studies have demonstrated that patients benefit from the application of a three-dimensional (3D) printed bolus. [12][13][14][15][16][17] However, studies that have suggested a dose reduction for the heart and lung refer only to electron beam radiotherapy for breast cancer patients. 15,16 Volumetric modulated arc therapy (VMAT) can significantly shorten the length of radiotherapy and improve the efficiency of equipment, and previous research studies have confirmed that VMAT technology can reduce the dose to the heart and ipsilateral lung. [17][18][19] However, there are no clinical practice data so far on whether VMAT has similar advantages in terms of dosimetry as a 3D-printed bolus on the chest wall. In this study, we focus on a comparison of the dosimetry characteristics, radiobiology, and acute toxicity of VMAT treatment using a 3D-printed bolus and a conventional bolus.

Study design and participants
A total of 75 eligible T1-4 and N1-3 breast cancer patients receiving postmastectomy chest wall radiation were recruited between March and December 2020 after approval was granted by the ethics board of our institution. Clinical staging was carried out based on the eighth edition of the TNM classification. The selection criteria included: (i) Karnofsky score ≥80; (ii) an intention to undergo VMAT treatment; and (iii) a requirement for a bolus during treatment. No patient was disenrolled from the clinical trial.
Each patient received a series of medical evaluations, including a medical history, a physical examination, a magnetic resonance imaging scan of the chest, blood testing, abdominal ultrasound, chest computed tomography (CT) scan, and a whole-body bone scan (if necessary) before enrollment.

| Randomization and treatment procedures
The random assignment of patients was performed at Nanchang University using a computer-generated random number code. This was an exploratory open-label study based on clinical feasibility. There was no formal sample size calculation, as no formal hypothesis test was applied. Depending on clinical feasibility, we planned to include about 70 cases. Patients were randomly distributed in a 1:1 ratio to the 3D-printed bolus and conventional bolus groups, and the left and right treatment sides formed a stratification factor. All patients understood the treatment assignment, and were fully informed of the available treatment option for bolus. While part of the patients enrolled with 3D printed bolus still opted for conventional bolus due to concerns about the 3D technology, and the appropriate treatment consent forms were signed and kept in the medical record. Thirty-five patients (19 right-sided and 16 left-sided) and 40 patients (20 right-sided and 20 left-sided) were enrolled and were assigned either a 3Dprinted bolus or a 5.0-mm thick conventional bolus.
All patients were placed in a supine position on a custom mold fixed onto a Solo Align Full Body System (Klarity) with the head tilted to the opposite side, and both upper arms were lifted in abduction. A conventional (standard vinyl gel sheet) bolus or a 3D-printed bolus was placed on the chest wall as uniformly as possible, to minimize the extent of any air gaps, and tape was applied to assist in conforming to the bolus to the concavities of the patient's surface as necessary. A thermoplastic mask was affixed to the patient to immobilize the neck and chest. The details of setting up the immobilization devices are shown in Figure 1.
The patients were given a CT scan under free breathing with a slice thickness of 5 mm (Siemens Medical Systems). Each clinical target volume (CTV) was delineated according to the breast cancer atlas for radiation therapy planning consensus definitions of the Radiation Therapy Oncology Group (RTOG). Planning target volumes (PTVs) were obtained from each CTV by expanding the 5 mm margin in three dimensions. PTV1 (the ipsilateral chest wall and internal mammary node regions) and PTV2 (the supraclavicular node regions) were restricted to 0 and 4 mm under the skin surface, respectively, and excluded the build-up region. Organs at risk (OARs) such as the ipsilateral lung, heart, contralateral breast and lung, liver, spinal cord, and trachea were outlined in axial CT sections. CBCT images were designed to verify the positioning error and the geometrical accuracy of the fit of the bolus. Age, tumor classification, and volumes of PTVs for all patients are given in Table 1.
The boluses in this study were designed based on CT images and manufactured using a 3D printer (Rehearsal), involving fused deposition modeling of polylactic acid (PLA) filament and 100% infill. 20 The design of the bolus for 3D printing was based on a series of CT images with a 1.5 mm slice, and the bolus area was marked with lead wires on the surface of the patient's body according to the extent of surgical clearance during CT scanning. The CT images were then imported into Mimics 19.0 (Materialise, Belgium) to create the 3D skin model. This model was imported into Geomagic Studio 2013 (Geomagi) in STL file format. The skin contours within the target area were extracted with 5 mm, and the thickness of the thin-walled area was increased appropriately based on cross-sectional CT images. The individualized 3D printed bolus created in this way was imported into the slicing software for layerby-layer slicing, and the resulting slicing file was input into  the 3D printer for printing. The electron density of the 3Dprinted bolus relative to water was set as 1.11 ± 0.01 (range 1.07-1.13) in our study, which is the same as in previous reports. 12 For each patient, VMAT plans were created using the Pinnacle treatment planning system (version 9.10) with 6 MV photon energy beams from a Varian TrueBeam linear accelerator. The prescription dose was 50 Gy in 25 fractions, and at least 95% of the PTV received 50 Gy. Dose-volume histograms (DVHs) were used for a dosimetry analysis. The homogeneity index (HI) of PTV1 and PTV2 was calculated 21 as HI = D 2% − D 98% ∕D 50% . The conformity index (CI) of the PTV as proposed by Paddick et al. 22 was evaluated as CI = (VPTV 50 ) 2 ∕VPTV × V 50 , where VPTV is the target volume, V 50 is the volume of the prescribed isodose value, and VPTV 50 is the volume of the target covered by the prescribed isodose value. The aim, starting objective, and constraints of OARs for two groups planning optimization are specified in Table 2.

| Outcomes and follow-up
MHD is an important practical parameter used to predict radiotherapy-induced heart disease. In a systematic review of heart doses, 23 the average MHD was round 8 Gy when the IMC was irradiated in left-sided breast cancer and 3.3 Gy for right-sided breast cancer. We previously reported that the MHD was 5.2 ± 0.9 Gy in left breast cancer using a virtual bolus with the VMAT technique. 19 Since a safe dose for the heart in clinical practice is always controversial, the lower the dose to the heart, the better, in principle. The primary endpoint was therefore a 10% decrease in MHD for left-sided breast cancer patients, defined as D mean (3D − heart) − D mean (conventional − heart)∕ D mean (conventional − heart) × 100% . This would mean that the average MHD for left-sided breast cancer would drop below 5 Gy. The secondary endpoint was a 5% decrease in the mean lung dose (MLD) to the ipsilateral lung, defined as The normal tissue complication probability (NTCP) for radiation-induced pneumonitis and mortality was computed for the lung and heart using the Lyman-Kutcher-Berman 24 and relative seriality models, respectively. 25 Toxicity was scored using the Radiation Therapy Oncology Group (RTOG) radiation morbidity scoring system. 26 No data on long-term toxicity were collected. Radiation pneumonitis toxicity was evaluated within 12 weeks after radiotherapy.

| Statistical analysis
An independent two-sample t-test was used to compare the air gaps between the bolus and skin and the dosimetry parameters. A chi-squared test was used for a comparison of radiation toxicity. Values of p ≤ 0.05 were considered statistically significant. All statistical analyses were twosided and were performed using SPSS version 18.0 (IBM Corporation).  Figure 2 shows CT and CBCT images of four patients with a conventional bolus (images A, C, E, and G) and a 3D-printed bolus (images B, D, F, and H), respectively. The 3D-printed bolus was superior, as smaller air gaps between the bolus and skin can be seen from the images. A frequency histogram for the maximum dimension of the air gap in the CT measurements is given in Figure 3. These distributions have different frequencies, where the frequency is 3.2 mm for the conventional bolus and 1.8 mm for the 3D-printed bolus. For all patients, the maximum mean air gap was 3.9 ± 1.4 mm for the conventional bolus and only 1.9 ± 0.9 mm for the 3D-printed bolus (p < 0.001). In addition, the air gap volume differed significantly (p < 0.001), with the 3D printed bolus resulting in lower values. Over all observations, the mean total air gap volumes were 6.5 ± 5.1 and 14.1 ± 8.8 ml for 3D printed bolus and for conventional bolus, respectively. These larger air gaps occurred at regions of surface complexity and, in particular, locations featuring large convexities or concavities (see Figure 2). Figure 4 shows average dose-volume histograms (DVHs) for OARs with all conventional bolus plans and all 3D-printed bolus plans. The relevant dosimetry parameters for the OARs (ipsilateral lung, heart, contralateral lung and breast, liver, and all normal tissue [Body-PTV]) are presented in Table 3. Compared to the conventional bolus, the MHD was reduced by an average of 14.55% (0.8 Gy) for the left-sided breast, while the V20 and mean dose for the ipsilateral lung were decreased by an average of 7.56% (1.7%) and 6.45% (0.8 Gy), respectively. No substantial differences were observed for the other OARs. PTV dosimetry parameters for the conventional and 3D-printed bolus plans for all patients are presented in Figure 5. For PTV1, the coverage of the 100% prescription dose was higher, with a score of 95.4 ± 0.6% versus 94.8 ± 0.8% (p = 0.001), and the HI in the 3D-printed bolus was better than for the conventional bolus (0.10 ± 0.01 vs. 0.11 ± 0.01, p = 0.019). No substantial differences were found for PTV2. For the overall PTV, the mean CI for the 3D-printed bolus was higher than for the conventional bolus (0.83 ± 0.02 vs. 0.80 ± 0.03, p < 0.001).
Compared with the conventional bolus, lower values for the NTCP of the heart and ipsilateral lung were also found for the 3D-printed bolus; these were reduced by an average of 0.14% (0.32 ± 0.19% vs. 0.18 ± 0.11%, p = 0.017) and 0.45% (3.70 ± 0.67% vs. 3.25 ± 0.18%, p = 0.001). No ≥Grade 2 radiation pneumonitis (RP) occurred with the 3D-printed bolus, whereas three patients from the conventional bolus group were identified to have ≥Grade 2 RP. With or without a 3D-printed bolus, Grade 1 RP occurred in 14.3% versus 20.0% of the cases (X 2 = 3.39, p = 0.184). Twenty-nine (38.7%) patients in total were identified as having ≥Grade 2 radiation dermatitis (RD) which occurred where the bolus was applied, and the difference in the occurrence rates of ≥Grade 2 RD between the conventional bolus (17/40) and the 3D-printed bolus (12/35) was not statistically significant (X 2 = 0.531, p = 0.466).

| DISCUSSION
This study aimed to show that the use of an individualized 3D-printed bolus combined with VMAT technology can reduce the irradiated dose to heart and lung in PMRT, resulting in a reduced incidence of acute toxicity. Our results revealed a 14.5% (0.8 Gy) decrease in the MHD for left-sided breast cancer treated with a 3D-printed bolus, compared to a conventional standard super-flat bolus. Our results are similar to those of previous studies on electronic electron beam technology for breast cancer patients. 15,16 In an earlier evaluation of clinical radiation treatment planning, the recommended dose limit for the heart was set at V 25 <10%. However, Beaton found that even if V 25 did not exceed 10% in all patients with left-sided breast cancer, the 10-year incidence of cardiovascular disease was still 1.7%. 27 He therefore suggested that the MHD should be used as a criterion for considering cardiac limits. The left anterior descending (LAD) artery is the branch of the coronary artery that is most likely to be damaged by radiotherapy, so some scholars have also proposed using the mean dose of the LAD to predict the incidence of major coronary events, although Shah found that this index was less predictive than MHD. 28 Darby proposed that the incidence of major coronary events increased linearly with increasing MHD, with a 7.4% increase per gray level. 4 In other words, compared to a conventional bolus, left-sided breast cancer patients treated with a 3D-printed bolus may be reduced by 5.92% in MHD. By looking at a meta-analysis of data from 75 trials in which patients were randomly assigned to breast cancer radiotherapy versus no radiotherapy, Taylor found a cardiac mortality ratio of 1.3 per 1 Gy of MHD. 29 Although Van believed that MHD alone might be obsolete as a risk assessment index for cardiac injury, and proposed a new index based on the volume of the left ventricle receiving 5 Gy as a predictor for an acute coronary event, 30 this index has not been widely used in clinical practice. Trott suggested that NTCP based on anatomical dose distribution can also be predictive of cardiopulmonary radiation risk. 31 In this study, the results for the 3D-printed bolus showed a lower NTCP on the heart (0.14%). For right-sided breast cancer, although the dose to the heart did not change between the 3D-printed

F I G U R E 2 Examples of conventional bolus and 3D-printed bolus CT (A-D) and CBCT images (E-H).
A and E were right-sided breast with the conventional bolus, B and F were right-sided breast with the 3D-printed bolus, C and G were left-sided breast with the conventional bolus, and D and H were left-sided breast with 3D printed bolus bolus and the conventional bolus, this value was lower than the results reported in the literature. 23 In addition to the lower dose to the heart, the dose to the ipsilateral lung was also reduced (6.45%, 0.8 Gy), which resulted in a decrease in NTCP. Although patients with a 3D-printed bolus had a lower incidence of Grade 1 and 2 RP, the difference in the outcomes between the two groups was not statistically significant. The frequency of symptomatic RP has been reported as 1%-7% after treating local and regional nodes with a total dose of 45-50 Gy. 32,33 We reduced the volume of irradiation dose to 5-30 Gy in the ipsilateral lung using the VMAT technique, especially in the 3D-printed bolus. In the current study, no patients in the 3D-printed bolus group developed Grade 2 RP. Furthermore, for clinical postmastectomy plans, the dose distribution for the PTV of the chest wall with a 3D-printed bolus did not suffer from the problems experienced by the conventional bolus. In our study, the PTV of the chest wall for the V 100% range was 91.7%-95.6% for the conventional bolus and 94.5%-97.8% for the 3Dprinted bolus. The results were also considerably more homogeneous than for a conventional bolus, and this improvement in the homogeneity of irradiation will reduce the rates of acute complications and long-term fibrosis. 4  Our results are inconsistent with previous results for the electronic irradiation of the chest wall, 34 which had a similar target dose coverage. The 3D-printed bolus reduces the dose to the heart and lung, and provides higher coverage of the chest wall, mainly because it gives a better fit and the bolus is thicker on a thinner chest wall. In terms of the accuracy of fit between the skin and the bolus, the 3D-printed bolus was superior to the conventional bolus as its outline was more consistent with the chest wall ( Figure 2). It is encouraging to see that the average dimension of the maximum air gap for both the conventional and the 3D-printed bolus was lower than previously reported, 7 and there were no larger air gaps (>10 mm) for either bolus. This could be ascribed to the use of a thermoplastic mold for immobilization. The radiotherapist pressed the thermoplastic mold above the bolus to improve the fit between the bolus and skin, and continued to press as the mold cooled. The average MHD and MLD using a 5-mm virtual bolus were 5.2 and 11.5 Gy, respectively, in our previous study, 19 which were lower than the values for the real conventional bolus of 5.5 and 12.4 Gy. This further demonstrates that air gaps increase the exposure dose to the heart and ipsilateral lungs.

T A B L E 3 Dosimetric difference between conventional bolus and 3D printed bolus on OARs
In addition to the effect of the air gap on the dose distribution of the OAR and targets, the thickness of the bolus will also change the dose to the skin. [35][36][37] Monajemi et al. showed that average skin doses in postmastectomy VMAT treatments with a 3, 5, and 10 mm bolus were 0.96, 1.03, and 1.04, respectively. 32 It should be noted that the thickness of the chest wall after mastectomy is uneven. The differences in the chest wall thickness for the majority of the patients who took part in this study were between 6 and 8 mm, and the largest exceeded 14 mm. It is therefore not appropriate to use boluses of uniform thickness. The results of previous studies 36 confirm that for a 3D conformal electron beam therapy of the chest wall after radical mastectomy, modifications to the patient's skin surface with a variable bolus can improve the coverage and homogeneity of the target dose. This method is also suitable for x-rays. Yoon et al. showed 37 that the dose reduction due to the bolus was more significant at shorter distances from the beam, and most of the dose reductions occurred in the first 2 cm of depth and stopped at a 4 cm of depth for IMRT or VMAT. We therefore designed the thickness of the bolus based on the inner profile of the patient's chest wall, so that the outline of the bolus would be similar to that of the chest wall, thus ensuring that there was sufficient bolus over the thinner chest wall for dose build-up. The average MHD using a 5 mm virtual bolus in our previous report 19 was also lower than for the 3Dprinted bolus (5.2 vs. 4.7 Gy), while the MLD was lower for the virtual bolus (115 vs. 11.6 Gy). It may be that the effect of the air gap is greater than the effect of the thickness of the bolus.
To further confirm that the air gap and the thickness of the bolus can modify the dose to the heart and ipsilateral lung, three patients with left-sided breast cancer for conventional bolus who had a high heart dose were selected. We modified the bolus by filling the air gap with virtual water material and adding the 5 mm thickness of the virtual bolus to the conventional bolus, and then re-optimizing it based on the treatment planning system. The new results showed that the MLD was reduced by 0.6, 1.2, and 0.7 Gy, respectively, and the ipsilateral MLD was reduced by 1.0, 0.9, and 0.9 Gy, respectively. More specifically, the coverage of the prescription dose for the PTV of the chest wall was increased by 0.1%, 2.5%, and 0.4%, respectively. The conformity index was increased by 0.04, 0.03, and 0.06, respectively, and there were no significant changes in the other OARs. These results mean that the completely fitted bolus and the thicker bolus over a thin chest wall can reduce the radiation dose to the heart and lung, and improve the target dose coverage of normal tissue on these areas. A custom 3D-printed bolus was therefore more suitable for radiotherapy following mastectomy, and also offered clinical benefits due to the improvement in PTV dosimetry and the decrease in the OAR dosimetry.
In this study, we found that there was no difference in the radiation dose to the liver between the conventional bolus and 3D-printed bolus for the right-sided breast. This is likely to be because only a small part of the liver was located behind the PTV of the chest wall, meaning that it had little effect on the overall radiation dose to the liver.
The RP was decreased for the 3D-printed bolus compared to the conventional bolus, but the difference was not statistically significant, which may have been due to the small sample used and the absolute value of differences in mean the ipsilateral lung between 3D-printed bolus and conventional bolus was small in this case. In addition, the probability of RP is relatively low for breast cancer, and it may not only be related to the dose of ipsilateral lung. The RD ≥Grade 2 (38.7%) was lower than that reported in other studies. [38][39][40] As in previous reports, the frequency of ≥Grade 2 RD was 40%-90% for patients treated with a conventional technique, 35 and 42% for patients treated with IMRT. 39 The results of Arsenault's study (38%) are similar to our findings. 40 Chen et al. 41 found that the risk factor for acute skin toxicity in adjuvant breast radiotherapy was mainly related to V 107% and V 110% . This corresponds to the dosimetry results of this study.
There are several potential limitations of this study that merit consideration. The primary limitations are the small sample size and the fact that fewer left-sided patients received the 3D-printed bolus. A further limitation relates to missing data for long-term clinical follow-up. Additional studies are therefore required to examine these issues.

| CONCLUSIONS
The use of a 3D-printed bolus in VMAT plans achieved better dose coverage and homogeneity in the PTV of the chest wall and the ipsilateral lung/heart irradiated dose, which resulted in reductions in the NTCP for those OARs and the rates of RD and RP. In conclusion, applications of new technologies such as 3D printing may enable more patients to reap the benefits of adjuvant radiation therapy for breast cancer treatment.