Skip to main content

Radiobiological and dosimetric comparison of 60Co versus 192Ir high-dose-rate intracavitary-interstitial brachytherapy for cervical cancer

Abstract

Background

High-dose-rate (HDR) intracavitary-interstitial brachytherapy (IC-ISBT) is an effective treatment for bulky, middle, and advanced cervical cancer. In this study, we compared the differences between 60Co and 192Ir HDR IC-ISBT plans in terms of radiobiological and dosimetric parameters, providing a reference for clinical workers in brachytherapy.

Methods

A total of 30 patients with cervical cancer receiving HDR IC-ISBT were included in this study, and IC-ISBT plans for each individual were designed with both 60Co and 192Ir at a prescribed dose of CTV D90 = 6 Gy while keeping the dose to OARs as low as possible. Physical dose and dose–volume parameters of CTV and OARs were extracted from TPS. The EQD2, EUBED, EUD, TCP, and NTCP were calculated using corresponding formulas. The differences between the 60Co and 192Ir IC-ISBT plans were compared using the paired t-test.

Results

In each patient's 60Co and 192Ir IC-ISBT plan, the average physical dose and EQD2 of 60Co were lower than those of 192Ir, and there were statistically significant differences in D2cc and D1cc for the OARs (p < 0.05); there were statistically significant differences in D0.1 cc for the bladder (p < 0.05) and no significant differences in D0.1 cc for the rectum or intestines (p > 0.05). The EUBED ratio (60Co/192Ir) at the CTV was mostly close to 1 when neither 60Co or 192Ir passed their half-lives or when both passed two half-lives, and the difference between them was not significant; at the OARs, the mean value of 60Co was lower than that of 192Ir. There was no statistical difference between 60Co and 192Ir in the EUD (93.93 versus 93.92 Gy, p > 0.05) and TCP (97.07% versus 97.08%, p > 0.05) of the tumors. The mean NTCP value of 60Co was lower than that of 192Ir.

Conclusions

Considering the CTV and OARs, the dosimetric parameters of 60Co and 192Ir are comparable. Compared with 192Ir, the use of 60Co for HDR IC-ISBT can ensure a similar tumor control probability while providing better protection to the OARs. In addition, 60Co has obvious economic advantages and can be promoted as a good alternative to 192Ir.

Introduction

Cervical cancer is the fourth leading cause of cancer death among women worldwide, and morbidity and mortality continue to increase each year in regions with poor access to prevention and screening measures [1, 2]. Cervical cancer has high radiosensitivity and is often treated with external beam radiotherapy (EBRT) and intracavitary/interstitial brachytherapy (IC/ISBT). For bulky, middle, and advanced cervical cancer, the treatment coverage of EBRT combined with ICBT has certain limitations [3,4,5], while in ISBT, an insertion needle can be placed according to the tumor shape and then used in combination with an intracavitary applicator to ensure that the central part of the cervix receives a high dose, covering the target area accurately and conformally so as to achieve more effective treatment. At present, cervical cancer brachytherapy is mainly based on a high dose rate (HDR). The two radioactive sources recommended in ICRU report 89 are 192Ir and 60Co. Previously, due to the large geometric size of 60Co and its limited clinical use, 192Ir occupied most of the market share. After continuous technical improvements, the geometry of miniaturized 60Co sources has become comparable to that of 192Ir sources. The 60Co (1.25 MeV) radiation source has higher gamma energy than 192 Ir (0.38 MeV), which indicates that a higher standard of medical protection is required for the use of 60Co. There are concerns that 60Co may have some damaging consequences distinct from those of 192Ir, such as excessive doses to normal organs and tissues and potentially increased toxicity. However, because 60Co has a long half-life (about 5 years) compared with 192Ir (74 days), it has obvious advantages in terms of human resources, logistics, and the economy. For developing countries, it is undoubtedly more cost-effective to choose a 60Co source [6, 7]. There are few reports on the differences between using 60Co and 192Ir for intracavitary-interstitial brachytherapy (IC-ISBT). Therefore, this study compares the differences in the physical dose and radiobiological effects on tumors and organs at risk (OARs) using these two radiation sources in IC-ISBT, conducting an analysis of the findings to provide a reference for clinical workers in brachytherapy.

Methods and materials

Patient selection and applicator insertion

Thirty patients diagnosed with cervical cancer (stage IB2 to IVA according to FIG.O 2018, squamous cell carcinoma) in Sichuan Cancer Hospital from January 1, 2021, to October 31, 2021, were selected for retrospective analysis. All patients underwent pelvic EBRT (45 Gy) and IC/IC-ISBT (6 Gy × 5 fractions). Before ISBT, the first CT scan was performed in a molded fixed position, ranging from the plane of the iliac spine to 2 cm below the vaginal opening, and the slice thickness was 3 mm. The insertion was immediately performed as a routine operation by three experienced clinicians, using a ProGuide sharp needle (Part# 189.601) and a Henschke titanium applicator (Part#110.437), and comprehensive evaluation of the patient's tumor location and size was performed to determine the required number of needles and the depth of insertion into the tumor tissue. A second CT scan was performed immediately after insertion, and the image data were transmitted to the Oncentra brachytherapy planning system. Before each CT scan, the catheter should be opened to empty the urine, and then 100 mL normal saline should be injected into the bladder to ensure that the filling state of the bladder remains unchanged. All patients were transported using a special transfer bed to prevent both position changes and movement of the implantation needle and metal applicator.

Contouring and treatment planning

We used the Oncentra treatment planning system (Elekta AB, Stockholm, Sweden, version 4.3) to design two brachytherapy plans on the CT for each patient, one using 60Co and one using 192Ir. The tumor and the OARs (including the bladder, rectum, sigmoid, and intestines) were outlined by oncologists according to the recommendations of the GEC-ESTRO working group. Then, a physicist performed applicator reconstruction, referring to the ICRU recommendations. The step length of the radiation source was 5 mm, the length of the metal applicator was 150 cm with an OFFSET value (the distance from the top of the applicator to the first dwelling position) of -0.6 cm, and the needle length was 124 cm with an OFFSET value of − 0.4 cm. All plans were optimized using the hybrid inverse planning optimization (HIPO) algorithm. When more needles are used, the HIPO algorithm can provide a better dose distribution in the high-dose region [8]. The dwell time gradient ratio (DTGR) was 0.6. The target parameters were continuously optimized so that the CTV D90 in the 60Co and 192Ir treatment plans of all patients reached 600 cGy (whereby the difference did not exceed 1 cGy), and the OAR D2cc value was as low as possible. The total dose of EBRT and HDR BT should meet the requirements of HR-CTV ≥ 85–90 Gy EQD2. For patients with severe disease, these requirements are HR-CTV D90 ≥ 87 Gy EQD2, rectum D2cc ≤ 65–75 Gy EQD2, sigmoid D2cc ≤ 70–75 Gy EQD2, and bladder D2cc ≤ 80–90 Gy EQD2 [9].

Parameters calculation

The use of the linear quadratic (LQ) model has been validated, and it has become the dominant model for radiobiology. The model proposes the use of α/β values in order to evaluate the radiotherapy sensitivity of tumors and normal tissues to fractionated irradiation doses and dose rate changes; increasing the dose of a single fraction of irradiation and reducing the number of fractions will significantly increase the radiation damage to late-reacting tissues [10]. The design of HDR BT plans based on the LQ model helps clinicians to conduct comprehensive evaluation of applying a single irradiation dose and determining the required number of fractions for patients. We extracted dose–volume data (D0.1 cc, D1cc, and D2cc) of HR-CTV and the OARs (the bladder, rectum, and intestines) from each patient's 60Co and 192Ir ISBT plans and strictly normalized the CTV D90 of all patients' ISBT plans to 600 cGy. The physical dose and equivalent dose in 2 Gy/f (EQD2) of the OARs were calculated for D0.1 cc, D1cc, and D2cc using Eq. (1).

$$EQD_{2} = D\frac{d + \alpha /\beta }{{2 + \alpha /\beta }}$$
(1)

where D is the total irradiation dose (fraction n × single irradiation dose d). The Groupe Europeen de Curietherapie and the European Society for Therapeutic Radiology and Oncology (GEC-ESTRO) working group recommends an α⁄β of 10 for cervical cancer (CTV) and 3 for normal tissues (OARs).

The biologically effective dose (BED) can be used for comparisons between different fraction schemes. In addition, in the LQ model, the G factor takes into account the time and dose changes, which can better describe the dose rate and cell sublethal damage repair ability through the effect on the β component [11]. The expressions of BED and the G factor are as follows:

$$BED = D\left( {1 + \frac{D}{\alpha /G\beta }} \right)$$
(2)
$$G = \left( {\frac{2}{{d^{2} }}} \right)\mathop \smallint \limits_{0}^{T} {\text{du}} I\left( u \right) \mathop \smallint \limits_{0}^{T} {\text{dw}} I\left( w \right) H\left( {u - w} \right) \times {\text{exp}}\left( { - \left( {u - w} \right)/T_{1/2} } \right)$$
(3)

where d is the dose of a single fraction; I(t) is the dose rate at time t; T1/2 is the repair half-time; and H(u-w) is the Heaviside function, defined as unity for u ≥ w and zero otherwise.

Both BED and EQD2 calculations assume a homogeneous dose distribution, but within the range of HDR BT treatment, there is anatomical heterogeneity in the irradiated tissue. Because each reference point is different from the radiation source, the dose rate and dose deposition differ, the dose distribution has obvious inhomogeneity, and the biological effects on the different parts also vary. Therefore, BED and EQD2 do not describe the radiobiological effects of HDR BT with sufficient precision. The concept of equivalent uniform biologically effective dose (EUBED) was thus introduced for improved description of the role of radiobiological effects arising from dose inhomogeneity in clinical outcomes [12,13,14,15].

$${\text{EUBED}} = - \frac{1}{{\alpha_{j} }}\ln \left( {\mathop \sum \limits_{i = 1}^{N} v_{i} e^{{ - \alpha \cdot {\text{BED}}}} } \right)$$
(4)

The parameter range of \({\alpha }_{j}\) is 0.05–0.5 Gy−1.

Mathematical models of normal tissue complication probability (NTCP) and tumor control probability (TCP) have been proposed for better optimization and evaluation of HDR BT plans. To address inhomogeneous dose distribution, Niemierko proposed the concept of an equivalent uniform dose (EUD) [16, 17], which refers to a uniform dose with the same probability of radiobiological damage as that produced by an inhomogeneous dose distribution. These can be calculated from the dose–volume data of the tumors and OARs, which are extracted from the TPS using the following expressions:

$$EUD = \left( {\mathop \sum \limits_{i = 1} \left( {V_{i} D_{i}^{a} } \right)} \right)^{\frac{1}{a}}$$
(5)

where a is a specific unitless model parameter, Vi is unitless, Di is the dose received by the voxel volume at i, and the sum of all partial volumes Vi is equal to 1. The local control rate of the tumor is likely to depend on the volume receiving the smallest dose, which represents the place with the highest clonal survival of tumor cells. Therefore, the EUD of the tumor will be close to the minimum dose, and parameter a should be a large negative number. In normal tissues with tandem structures, the EUD will be close to the maximum dose, and parameter a is usually a large positive number [18,19,20]. atumor = − 10,abladder = 2,arectum = 8.33,aintestines = 6。Based on the concept of EUD, Niemierko proposed the TCP/NTCP model:

$${\text{TCP}} = \frac{1}{{1 + \left( {\frac{{{\text{TCD}}_{50} }}{{{\text{EUD}}}}} \right)^{{4\gamma_{50} }} }}$$
(6)
$${\text{NTCP}} = \frac{1}{{1 + \left( {\frac{{{\text{TD}}_{50} }}{{{\text{EUD}}}}} \right)^{{4\gamma_{50} }} }}$$
(7)

We used the published code of Gay, H. A. et al. for calculating TCP/NTCP[20], which takes into account the total dose of radiation therapy (EBRT 45 Gy + IC-ISBT 6 Gy × 5 fractions). TCD50 is the irradiation dose received at a tumor control rate of 50%, \({\gamma }_{50}\)=3; TD50 is the dose tolerated by normal tissue at a complication rate of 50% for OARs (the bladder, rectum, and intestines), \({\gamma }_{50}\)=4[21, 22].

Statistics

We extracted the dose–volume data from the TPS of the CTV and OARs (the bladder, rectum, and small bowel) at D2cc, D1cc, and D0.1 cc for all patients with 60Co versus 192Ir BT plans. The differences in brachytherapy between the two sources, 60Co and 192Ir, were evaluated by comparing their dosimetric and radiobiological parameters. We used the paired-sample t-test to evaluate these differences on the premise that the samples met an approximately normal distribution, and the difference was statistically significant at p < 0.05. All statistical analyses were performed using IBM SPSS Statistics, version 22. (Table 1).

Table 1 Constraints for the HIPO algorithm

Results

Table 2 shows the mean ± SD of the dose values of D2cc, D1cc, and D0.1 cc in the OARs (the bladder, rectum, and intestines) for 60Co and 192Ir for a single ISBT irradiation dose in 30 patients. As can be seen in Table 2, as the cumulative dose values of each OAR gradually increased from D2cc to D0.1 cc, the difference between the mean values of the 60Co and 192Ir plans gradually decreased, and the dose values of the 60Co plans were all lower than those of the 192Ir plans. At D2cc and D1cc, the differences between the physical doses of 60Co versus 192Ir were statistically different in the bladder, rectum, and intestines (p < 0.05); at D0.1 cc, the differences when comparing 60Co and 192Ir in the bladder remained statistically significant (p < 0.05), whereas there was no statistically significant difference (p > 0.05) in the rectum (p = 0.240) or intestines (p = 0.298). In the 60Co and 192Ir plans, the difference in the mean physical dose to the bladder was the largest among the OARs: at D2cc, the difference was 0.11 Gy, and 60Co was 2.29% lower than 192Ir; at D1cc, the difference was 0.11 Gy, and 60Co was 2.16% lower than 192Ir; and at D0.1 cc, the difference was 0.1 Gy, and 60Co was 1.73% lower than 192Ir.

Table 2 Differences between 60Co and 192Ir according to physical dose to OARs (the bladder, rectum, and intestines)

Table 3 presents the mean ± SD of the EQD2 of the 30 patients, and they were calculated using formula (1). The calculation results are based on the physical dose, so the data distribution trend in Table 3 is the same as that in Table 2. At D2cc, D1cc, and D0.1 cc of each OAR, EQD2 was lower in the 60Co than in the 192Ir plans. At D2cc and D1cc, for 60Co versus 192Ir, there was a statistically significant difference in the physical dose for each OAR (p < 0.05); at D0.1 cc, there was a statistically significant difference in the bladder (p < 0.05), while there was no statistically significant difference in the rectum (p = 0.265) or the intestines (p = 0.354) (p > 0.05). In the plans of the two radioactive sources, the difference in the mean EQD2 of the bladder was the largest among the OARs.

Table 3 Differences between 60Co and 192Ir according to EQD2 of OARs (the bladder, rectum, and intestines)

Figure 1 displays a box plot of the physical dose and EQD2 received by the bladder, rectum, and intestines of the 30 patients. The boxes show the median value and 5%-95% range, the cross symbols indicate the maximum and minimum values, and the plus sign indicates the mean value. Figure 1 shows that, in general, the intestines receive a lower dose than the bladder and rectum, and there are large individual differences among the different patients. The average value of 60Co is lower than that of 192Ir, and the difference between them is small.

Fig. 1
figure 1

A box plot of physical dose a and EQD2 b received by the bladder, rectum, and intestines for 30 patients. Boxes show median value and 5–95% range, cross symbols indicate maximum and minimum values, and the plus sign indicates the mean value

The half-lives of 60Co and 192Ir are quite different. Considering the decay of the radioactive source, we calculated the EUBED values generated by the IC-ISBT plans of the 30 patients at different activity levels. Figure 2 shows the results of the EUBED calculation for all patients, and Table 4 shows the mean ± SD EUBED values for the CTVs and OARs. As can be seen in the above graphs, the bioeffective doses considering the CTV and OARs gradually decrease with the decay of the radioactive sources. In addition, we calculated the EUBED ratios (60Co/192Ir) at different activity levels (Fig. 3). Regarding CTV, when neither 60Co or 192Ir pass their half-life (as in Fig. 3. (a), (b), (d), (e)), or when both pass two half-lives (as in Fig. 3. (h)), their EUBED ratios are mostly close to 1, and the difference between them is not obvious; when 60Co is 1 or 2 Ci and 192Ir is 2 Ci, as in Fig. 3. (c) and (f), their EUBED ratios are mostly greater than 1, and 60Co produces a higher EUBED; when 192Ir is 8 or 4 Ci and 60Co is 0.5 Ci, as in Fig. 3. (g) and (h), their EUBED ratios are mostly less than 1, and 192Ir produces a higher EUBED.

Fig. 2
figure 2

A box plot of EUBED values for the CTV, bladder, rectum, and intestines at different activity levels for 60Co and 192Ir. Boxes show median value and 5–95% range, cross symbols indicate maximum and minimum values, and the plus sign indicates the mean value

Table 4 EUBED for the CTV, bladder, rectum, and intestines at different activity levels for 60Co and 192Ir
Fig. 3
figure 3

The EUBED ratio produced by 60Co and 192Ir at different activity levels

Table 5 shows the mean EUD for the CTV and OARs for the 30 patients. The statistical results are based on the total dose of EBRT (45 Gy) and IC-ISBT (6 Gy × 5 fractions). As can be seen in the Table 5, there was no statistically significant difference in EUD between the 60Co and 192Ir plans in tumors (p > 0.05), while there was a statistically significant difference in the OARs (p < 0.05), with 60Co being 1.26% lower than 192Ir in the bladder, 1.31% lower in the rectum, and 0.82% lower in the intestines.

Table 5 EUD values of 60Co and 192Ir for the CTV, bladder, rectum, and intestines

Table 6 shows the mean TCP and NTCP values in the 60Co and 192Ir plans for the 30 patients. The tumor control probability values of 60Co and 192Ir both reached 97% and were not statistically significantly different (p > 0.05). The NTCP of the bladder (0.27 for 60Co and 0.33 for 192Ir) was the lowest of all the OARs; the use of 60Co was 18.18% lower than that of 192Ir, and there was a statistical difference between the two sources (p < 0.05). The NTCP values of 60Co and 192Ir were significantly different in the rectum and intestines (p < 0.05), where 60Co was 16.87% and 8.23% lower than 192Ir, respectively.

Table 6 Mean tumor control probability and normal tissue complication probability in 60Co and 192Ir plans for 30 patients

Discussion

HDR BT is currently an important component of cervical cancer treatment. For bulky, middle, and advanced cervical cancer, IC-ISBT can better conform to the tumor area and has excellent therapeutic effects [23]. The dose distribution of brachytherapy is characterized by a high dose close to the center of the radiation source, rapid drop in the dose to the surrounding area, and uneven distribution of the actual brachytherapy dose due to heterogeneity of the human anatomical structure. The ICRU recommends the use of rectal and bladder reference points for the dose assessment of normal organ tissues in HDR BT; however, with the continuous development of 3D BT technology, studies have confirmed that a DVH dose assessment based on CTV and OARs often differs from those based on ICRU reference points. DVH-based dose assessment results are more conformal and reliable, while TCP and NTCP calculations can be performed based on the results, which is conducive to individually designed radiotherapy plans and improved tumor control rates [24,25,26].

Radiation sources 60Co and 192Ir have similar dose distribution trends. Many previous studies have demonstrated that the dose parameters of the HR-CTV, bladder, and rectum in brachytherapy plans of 60Co are comparable to those of 192Ir. Although the dose line distribution in the cephalocaudal direction is more prominent for 60Co than for 192Ir, they essentially become the same after optimization of the radiotherapy plan by a physicist, with the difference reduced to the extent that it is no longer clinically significant [27,28,29,30]. In this study, we selected 30 patients and designed two radiotherapy plans of 60Co and 192Ir for each individual patient, and all plans were given the prescribed dose of CTV D90 6 Gy. The data extracted from TPS were processed using calculations, and the physical dose values of CTV D90 and EQD2 were equal for both the 60Co and 192Ir plans. Table 2 shows that there is a statistically significant difference between the two sources in the OAR D2cc and D1cc doses, with 60Co producing a lower mean physical dose than 192Ir. This may be due to the steeper dose gradient of 60Co. At distances less than 1 cm from the source, the 60Co source produced higher dose rates and, therefore, higher cumulative physical dose values, while at distances greater than 1 cm from the source, the 192Ir source produced higher dose rates [31]. Moreover, the bladder and rectum had higher doses than the intestines due to their location being closer to the source of radiation. It has previously been demonstrated that 60Co produces a larger hot spot than 192Ir in the CTV of tumors and that the dose rate falls more rapidly as the distance from the source increases [32,33,34]. Compared to the published results, the internal, external, and obturator lymph nodes etc., (locations far away from the center of the source) receive a lower dose when using 60Co for ICBT [35]. These OARs are located relatively far away from the source of radiation and, therefore, receive less exposure. The results for D2cc, D1cc, and D0.1 cc not only represent dose–volume data in treatment plans but also allow for an overall dose assessment and complication prediction for OARs [36,37,38]. Based on the analysis of the experimental data, we speculate that the use of 60Co in IC-ISBT will provide better protection to OARs than 192Ir.

Cervical cancer has a strong self-healing ability. Although a large number of cells may proliferate in the interval between two treatments, hypoxic cells may also increase, resulting in increased radioresistance. However, because the dose distribution of interstitial implants may be more adequate, the dose of brachytherapy is high enough to overcome the increased radiation resistance caused by hypoxia, thereby maintaining a good tumor control rate [39, 40]. High dose rates of brachytherapy above 12 Gy/h are very beneficial for tumor tissue regression. It can be seen that, in order to achieve better therapeutic effects, HDR BT has certain requirements in terms of the activity of the radioactive source. It can be seen in Fig. 3 that the dose rates of the two radioactive sources, whether 60Co or 192Ir, show the same trend, and a decrease in activity leads to a decrease in EUBED. In the case of higher activity of the radioactive source, there was no significant difference in the EUBED values produced by the two sources at the CTV. In the OARs, the EUBED of 60Co was significantly lower than that of 192Ir, which may be due to the fact that the OARs received a lower radiation dose and were farther away from the radioactive source in response to 60Co, thus resulting in a lower dose rate. The half-life of 60Co is about 5.27 years, and the half-life of 192Ir is about 74 days. In order to ensure the therapeutic effect, 192Ir should be replaced about every six months, and 60Co thus has obvious advantages from the perspective of economic efficiency.

When designing a radiotherapy plan, in addition to optimization based on the distribution of physical doses, another important evaluation method is based on TCP and NTCP. In the treatment of cervical cancer, TCP is positively correlated with irradiation dose. Increasing the dose rate and shortening the total treatment time can improve TCP, but at the same time, the risks of normal tissue damage and late complications increase with dose [41]. In this study, we chose the model proposed by Niemierko, which is calculated based on EUD (Eq. 5). The calculation for EQD2 in the previous paper was based on D90 and did not take into account the non-uniform dose distribution in brachytherapy, so it is possible that the calculated value was lower than the actual value. The EUD weighting factor was introduced in Eq. 5, which can better describe the tumor area with a non-uniform dose distribution and has a better fitting ability. In Table 5, the EUD values in the tumors for 60Co and 192Ir are not statistically different, indicating that the choice of radiation source has no significant effect on the equivalent uniform dose distribution at the CTV; for the OARs, the use of 60Co compared with 192Ir produces lower mean EUD values. However, since the actual physical dose distribution of the OARs is more heterogeneous than the tumor area in HDR BT, routine EUD calculation and evaluation of the OARs are not recommended, and D2cc, D1cc, or D0.1 cc is recommended instead [36]. Table 2 shows that for TCP, the use of 60Co does not cause results that are significantly different from those of 192Ir; in the OARs with statistical differences, the mean NTCP is slightly lower for 60Co than for 192Ir, especially regarding the NTCP of the bladder, which is reduced by 18.18%. Intestinal cells have a higher ionizing radiation sensitivity and a higher probability of complications (enterocolitis, intestinal fistula, etc.) after radiotherapy, and the predicted outcomes of 60Co and 192Ir have the same trend of distribution, so there is no disadvantage of 60Co compared to 192Ir in terms of disease control. In published clinical trials, there are results showing that patients with cervical cancer treated with 60Co HDR ISBT have similar survival and toxicity outcomes to those treated with 192Ir [42,43,44]. However, there is still a lack of large sample sizes and randomized clinical trials using 60Co and 192Ir for IC-ISBT in patients with cervical cancer to establish definitive conclusions. Therefore, further research is needed to investigate the differences in clinical outcomes, such as the local control rate and survival time, between the two radioactive sources.

In this study, data from 30 patients treated with a single ISBT irradiation dose were selected to investigate the radiobiological differences between 60Co and 192Ir, but in the actual clinical treatment process, the biological characteristics of tumors and normal tissues are not only affected by a single irradiation. The following five major biological factors underlie the rationale for fractionated radiation therapy: radiosensitivity, repair, regeneration, redistribution, and reoxygenation. After each radiotherapy, the tumor tissue regression and patient's physical tolerance will influence the next IC-ISBT. An increase in tumor hypoxic cells after treatment also increases radioresistance to the next treatment. Residual cells after irradiation may “repopulate”, and if a sufficiently large number of tumor cells repopulate in the interval between fractions, the treatment may fail. The repair of sublethal damage to normal tissues may not be completed in time, so the control of the tumor cannot be separated from the time of the entire radiotherapy [45]. In the study conducted by Dayyani, M. et al., considering the influence of RBE correction and a patient's full treatment regimen, the BED produced by 60Co was lower, so it was recommended that the prescribed dose be increased by 4% when using a 60Co source; the EQD2 of OARs after dose escalation was close to or even lower than that of 192Ir [33]. However, our study was mainly conducted to compare the physical dose and radiobiological differences in a single IC-ISBT irradiation dose with regard to the choice of radioactive source, so the other treatments, including the EBRT that was performed and the remaining four ISBTs, were treated uniformly and were simplified. The statistical results can provide a reference for brachytherapy workers. In order to obtain higher level clinical evidence, this comparative trial needs to be further extended to address the aforementioned issues.

Conclusion

By comparing the results of 192Ir and 60Co IC-ISBT plans, it can be concluded that, after constant optimization and given the same constraints, there is no significant difference in the radiobiological effect parameters and TCP at the CTV between the two radioactive sources. For the OARs, the average physical dose and EQD2 of 60Co in D2cc were lower than those of 192Ir, the average NTCP of the rectum and intestines was also lower, and the NTCP of the bladder was not significantly different. Therefore, it is speculated that 60Co may have a better protective effect on OARs without reducing the dose at the CTV. In addition, 60Co is cost effective, so the use of 60Co as an alternative to 192Ir for IC-ISBT it feasible and advantageous.

Availability of data and materials

The data used during this study are available from the corresponding author on reasonable request.

Abbreviations

HDR:

High dose rate

IC-ISBT:

Intracavitary-interstitial brachytherapy

OARs:

Organs at risk

CTV:

Clinical target volume

EBRT:

External beam radiotherapy

TPS:

Treatment planning system

HIPO:

Hybrid inverse planning optimization

DTGR:

Dwell time gradient ratio

BED:

Biologically effective dose

EQD2:

Equivalent dose in 2 Gy/f

EUBED:

Equivalent uniform biologically effective dose

EUD:

Equivalent uniform dose

TCP:

Tumor control probability

NTCP:

Normal tissue complication probability

References

  1. Xia C, Dong X, Li H, Cao M, Sun D, He S, et al. Cancer statistics in China and United States, 2022: profiles, trends, and determinants. Chin Med J (Engl). 2022;135(5):584–90.

    Article  Google Scholar 

  2. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.

    Article  Google Scholar 

  3. Tanderup K, Eifel PJ, Yashar CM, Potter R, Grigsby PW. Curative radiation therapy for locally advanced cervical cancer: brachytherapy is NOT optional. Int J Radiat Oncol Biol Phys. 2014;88(3):537–9.

    Article  Google Scholar 

  4. Barraclough LH, Swindell R, Livsey JE, Hunter RD, Davidson SE. External beam boost for cancer of the cervix uteri when intracavitary therapy cannot be performed. Int J Radiat Oncol Biol Phys. 2008;71(3):772–8.

    Article  Google Scholar 

  5. Han K, Milosevic M, Fyles A, Pintilie M, Viswanathan AN. Trends in the utilization of brachytherapy in cervical cancer in the United States. Int J Radiat Oncol Biol Phys. 2013;87(1):111–9.

    Article  Google Scholar 

  6. Strohmaier S, Zwierzchowski G. Comparison of (60)Co and (192)Ir sources in HDR brachytherapy. J Contemp Brachytherapy. 2011;3(4):199–208.

    Article  Google Scholar 

  7. Granero D, Perez-Calatayud J, Ballester F. Technical note: dosimetric study of a new Co-60 source used in brachytherapy. Med Phys. 2007;34(9):3485–8.

    Article  CAS  Google Scholar 

  8. Frohlich G, Geszti G, Vizkeleti J, Agoston P, Polgar C, Major T. Dosimetric comparison of inverse optimisation methods versus forward optimisation in HDR brachytherapy of breast, cervical and prostate cancer. Strahlenther Onkol. 2019;195(11):991–1000.

    Article  Google Scholar 

  9. Abu-Rustum NR, Yashar CM, Bean S, Bradley K, Campos SM, Chon HS, et al. NCCN guidelines insights: cervical cancer version 1.2000. J Natl Compr Canc Netw. 2020;18(6):660–6.

    Article  Google Scholar 

  10. van Leeuwen CM, Oei AL, Crezee J, Bel A, Franken NAP, Stalpers LJA, et al. The alfa and beta of tumours: a review of parameters of the linear-quadratic model, derived from clinical radiotherapy studies. Radiat Oncol. 2018;13(1):96.

    Article  Google Scholar 

  11. Manning MA, Zwicker RD, Arthur DW, Arnfield M. Biologic treatment planning for high-dose-rate brachytherapy. Int J Radiat Oncol Biol Phys. 2001;49(3):839–45.

    Article  CAS  Google Scholar 

  12. Barendsen GW. Dose fractionation, dose rate and iso-effect relationships for normal tissue responses. Int J Radiat Oncol Biol Phys. 1982;8(11):1981–97.

    Article  CAS  Google Scholar 

  13. Niemierko A. Reporting and analyzing dose distributions: a concept of equivalent uniform dose. Med Phys. 1997;24(1):103–10.

    Article  CAS  Google Scholar 

  14. Jones LC, Hoban PW. Treatment plan comparison using equivalent uniform biologically effective dose (EUBED). Phys Med Biol. 2000;45(1):159–70.

    Article  CAS  Google Scholar 

  15. Fields EC, Melvani R, Hajdok G, D’Souza D, Jones B, Stuhr K, et al. A multi-institution, retrospective analysis of cervix intracavitary brachytherapy treatments part 1: Is EQD2 good enough for reporting radiobiological effects? Int J Radiat Oncol Biol Phys. 2017;99(1):219–26.

    Article  Google Scholar 

  16. Niemierko A. A generalized concept of equivalent uniform dose (EUD). Med Phys. 1999;26(6):1100.

    Google Scholar 

  17. Niemierko A, editor A unified model of tissue response to radiation. Proceedings of the 41th AAPM annual meeting; 1999: Nashville: Wikipedia.

  18. Niemierko A, Goitein M. Modeling of normal tissue response to radiation: the critical volume model. Int J Radiat Oncol Biol Phys. 1993;25(1):135–45.

    Article  CAS  Google Scholar 

  19. Niemierko A, Goitein M. Implementation of a model for estimating tumor control probability for an inhomogeneously irradiated tumor. Radiother Oncol. 1993;29(2):140–7.

    Article  CAS  Google Scholar 

  20. Gay HA, Niemierko A. A free program for calculating EUD-based NTCP and TCP in external beam radiotherapy. Phys Med. 2007;23(3–4):115–25.

    Article  Google Scholar 

  21. Oinam AS, Singh L, Shukla A, Ghoshal S, Kapoor R, Sharma SC. Dose volume histogram analysis and comparison of different radiobiological models using in-house developed software. J Med Phys. 2011;36(4):220–9.

    Article  Google Scholar 

  22. Sukhikh ES, Sukhikh LG, Lushnikova PA, Tatarchenko MA, Abdelrahman AR. Dosimetric and radiobiological comparison of simultaneous integrated boost and sequential boost of locally advanced cervical cancer. Phys Med. 2020;73:83–8.

    Article  CAS  Google Scholar 

  23. Tiwari R, Narayanan GS, Narayanan S, Suresh KP. Long-term effectiveness and safety of image-based, transperineal combined intracavitary and interstitial brachytherapy in treatment of locally advanced cervical cancer. Brachytherapy. 2020;19(1):73–80.

    Article  Google Scholar 

  24. Hashim N, Jamalludin Z, Ung NM, Ho GF, Malik RA, Phua VC. CT based 3-dimensional treatment planning of intracavitary brachytherapy for cancer of the cervix: comparison between dose-volume histograms and ICRU point doses to the rectum and bladder. Asian Pac J Cancer Prev. 2014;15(13):5259–64.

    Article  Google Scholar 

  25. Paul AG, Nalichowski A, Abrams J, Paximadis P, Zhuang L, Miller S. Dosimetric evaluation of point A and volume-based high-dose-rate plans: a single institution study on adaptive brachytherapy planning for cervical cancer. J Contemp Brachytherapy. 2018;10(3):202–10.

    Article  Google Scholar 

  26. Srivastava S, Painuly NK, Mishra SP, Srivastava K, Singh N, Singh S, et al. Evaluation of volumetric doses of organs at risk in carcinoma cervix patients with HDR intracavitary brachytherapy and comparison of CT-based and conventional plans. J Biomed Phys Eng. 2019;9(6):603–12.

    CAS  Google Scholar 

  27. Shukla AK, Jangid PK, Rajpurohit VS, Verma A, Dangayach SK, Gagrani V, et al. Dosimetric comparison of (60)Co and (192)Ir high dose rate source used in brachytherapy treatment of cervical cancer. J Cancer Res Ther. 2019;15(6):1212–5.

    Article  Google Scholar 

  28. Sinnatamby M, Nagarajan V, Kanipakam Sathyanarayana R, Karunanidhi G, Singhavajala V. Study of the dosimetric differences between (192)Ir and (60)Co sources of high dose rate brachytherapy for breast interstitial implant. Rep Pract Oncol Radiother. 2016;21(5):453–9.

    Article  Google Scholar 

  29. Sinnatamby M, Kandasamy S, Karunanidhi G, Neelakandan V, Ramapandian S, Kannan M, et al. Image-guided brachytherapy a comparison between 192Ir and 60Co sources in carcinoma uterine cervix. Gulf J Oncolog. 2022;1(38):7–14.

    Google Scholar 

  30. Richter J, Baier K, Flentje M. Comparison of 60cobalt and 192iridium sources in high dose rate afterloading brachytherapy. Strahlenther Onkol. 2008;184(4):187–92.

    Article  Google Scholar 

  31. Perez-Calatayud J, Ballester F, Das RK, Dewerd LA, Ibbott GS, Meigooni AS, et al. Dose calculation for photon-emitting brachytherapy sources with average energy higher than 50 keV: report of the AAPM and ESTRO. Med Phys. 2012;39(5):2904–29.

    Article  Google Scholar 

  32. Safigholi H, Meigooni AS, Song WY. Comparison of (192) Ir, (169) Yb, and (60) Co high-dose rate brachytherapy sources for skin cancer treatment. Med Phys. 2017;44(9):4426–36.

    Article  CAS  Google Scholar 

  33. Dayyani M, Hoseinian-Azghadi E, Miri-Hakimabad H, Rafat-Motavalli L, Abdollahi S, Mohammadi N. Radiobiological comparison between Cobalt-60 and Iridium-192 high-dose-rate brachytherapy sources: Part I-cervical cancer. Med Phys. 2021;48(10):6213–25.

    Article  CAS  Google Scholar 

  34. Safigholi H, Han DY, Mashouf S, Soliman A, Meigooni AS, Owrangi A, et al. Direction modulated brachytherapy (DMBT) for treatment of cervical cancer: a planning study with (192) Ir, (60) Co, and (169) Yb HDR sources. Med Phys. 2017;44(12):6538–47.

    Article  CAS  Google Scholar 

  35. Siavashpour Z, Aghili M, Anjidani S, Zayeri F, Moghani MM, Maleki M, et al. Dose to pelvic lymph nodes during brachytherapy of locally advanced cervical cancer with 60Co HDR source. Brachytherapy. 2021;21(2):158.

    Article  Google Scholar 

  36. Potter R, Haie-Meder C, Van Limbergen E, Barillot I, De Brabandere M, Dimopoulos J, et al. Recommendations from gynaecological (GYN) GEC ESTRO working group (II): concepts and terms in 3D image-based treatment planning in cervix cancer brachytherapy-3D dose volume parameters and aspects of 3D image-based anatomy, radiation physics, radiobiology. Radiother Oncol. 2006;78(1):67–77.

    Article  Google Scholar 

  37. Haie-Meder C, Potter R, Van Limbergen E, Briot E, De Brabandere M, Dimopoulos J, et al. Recommendations from gynaecological (GYN) GEC-ESTRO working group (I): concepts and terms in 3D image based 3D treatment planning in cervix cancer brachytherapy with emphasis on MRI assessment of GTV and CTV. Radiother Oncol. 2005;74(3):235–45.

    Article  Google Scholar 

  38. Zhou YC, Zhao LN, Wang N, Hu J, Sun XH, Zhang Y, et al. Late rectal toxicity determined by dose-volume parameters in computed tomography-based brachytherapy for locally advanced cervical cancer. Cancer Med. 2016;5(3):434–41.

    Article  Google Scholar 

  39. Liu ZS, Guo J, Zhao YZ, Lin X, Chen B, Zhang M, et al. Salvage interstitial brachytherapy based on computed tomography for recurrent cervical cancer after radical hysterectomy and adjuvant radiation therapy: case presentations and introduction of the technique. J Contemp Brachytherapy. 2016;8(5):415–21.

    Article  Google Scholar 

  40. Lindblom E, Dasu A, Beskow C, Toma-Dasu I. High brachytherapy doses can counteract hypoxia in cervical cancer-a modelling study. Phys Med Biol. 2017;62(2):560–72.

    Article  CAS  Google Scholar 

  41. Lambin P, Gerbaulet A, Kramar A, Scalliet P, Haie-Meder C, Malaise EP, et al. Phase III trial comparing two low dose rates in brachytherapy of cervix carcinoma: report at two years. Int J Radiat Oncol Biol Phys. 1993;25(3):405–12.

    Article  CAS  Google Scholar 

  42. Tantivatana T, Rongsriyam K. Treatment outcomes of high-dose-rate intracavitary brachytherapy for cervical cancer: a comparison of Ir-192 versus Co-60 sources. J Gynecol Oncol. 2018;29(5): e86.

    Article  Google Scholar 

  43. Abtahi M, Gholami S, Nashi HHJ. High dose rate(192)Ir versus high dose rate(60)Co brachytherapy: an overview of systematic reviews of clinical responses of gynecological cancers from 1984 to 2020. Biomed Phys Eng Express. 2021;7(5):55022.

    Article  Google Scholar 

  44. Kumar M, Thangaraj R, Alva RC, Koushik K, Ponni A, Janaki MG. Interstitial high-dose-rate brachytherapy using cobalt-60 source for cervical cancer: dosimetric and clinical outcomes from a single institute. J Contemp Brachytherapy. 2020;12(4):351–5.

    Article  Google Scholar 

  45. Mazeron R, Castelnau-Marchand P, Dumas I, del Campo ER, Kom LK, Martinetti F, et al. Impact of treatment time and dose escalation on local control in locally advanced cervical cancer treated by chemoradiation and image-guided pulsed-dose rate adaptive brachytherapy. Radiother Oncol. 2015;114(2):257–63.

    Article  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This study was in part funded by a research grant (No: 2021-DF-016) of Elekta AB (Stockholm, Sweden). The funders had no role in study design, data collection and analysis, and decisions on preparation of the manuscript. This study was partially supported by Sichuan Science and Technology Program (No. 2022YFG0194, No. 2022YFS0047, No. 2021YFG0320) and Medical Engineering Innovation Fund for Cancer (No. ZYGX2021YGCX002).

Author information

Authors and Affiliations

Authors

Contributions

AW and XW designed the study and collected the data with support from PW and JL. AW, CY, BW and JL designed all the brachytherapy plans and did statistical analysis. AW wrote the manuscript. All co-authors read and approved the final manuscript.

Corresponding author

Correspondence to Pei Wang.

Ethics declarations

Ethics declarations consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wen, A., Wang, X., Wang, B. et al. Radiobiological and dosimetric comparison of 60Co versus 192Ir high-dose-rate intracavitary-interstitial brachytherapy for cervical cancer. Radiat Oncol 17, 206 (2022). https://doi.org/10.1186/s13014-022-02170-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13014-022-02170-8

Keywords