Correlation of chemokines and growth factors with radiation-induced liver injury after interstitial high dose rate (HDR) brachytherapy of liver metastases.

BACKGROUND: Locoregional therapies, as imaging-guided tumor-directed procedures, are emerging treatment strategies in the management of primary and secondary liver malignancies such as e.g. colorectal cancer liver metastases. As one of those, irradiation-based interstitial high dose rate brachytherapy (iBT) of liver metastases bears a risk of developing focal radiation-induced liver injury (fRILI). Since little is known about biological factors involved in hepatic dysfunction after irradiation, the aim of this study was to identify factors, that may play a role in the underlying mechanism of fRILI, and that potentially may serve as biomarkers for post-therapeutic fRILI to improve specific management and treatment of patients.
METHODS: Twenty-two patients with hepatic malignancies (tumor patients, TP) underwent iBT with total ablative doses of radiation to the target volume ranging from e.g. 15 to 25 Gy. Hepatobiliary magnetic resonance imaging (MRI) was performed 6 weeks after iBT to quanitify fRILI. Blood samples were taken before (pre) and 6 weeks after (post) iBT from TP, and from ten healthy volunteers (HV controls) for the analyses of humoral mediators: monocyte chemoattractant protein-1 (MCP-1), chemokine (C-X3-C motif) ligand 1 (CX3CL1), vascular endothelial growth factor (VEGF) and beta-nerve growth factor (beta-NGF) using the Multi-Analyte Flow Assay via flow cytometry. Correlation analyses between the humoral mediators (pre and post iBT) with the tumor volume and fRILI were performed.
RESULTS: While MCP-1 and CX3CL1 tended to decrease in TP vs. HV, VEGF was significantly decreased in TP vs. HV pre and post iBT (p < 0.05). Beta-NGF levels were significantly increased in TP vs. HV pre and post iBT (p < 0.05). Baseline circulating levels of MCP-1, VEGF and beta-NGF have shown significant positive correlations with the hepatic tumor volume (p < 0.05). Circulating levels of humoral mediators before treatment did not correlate with fRILI, while CX3CL1 and VEGF after iBT have shown significant positive correlations with fRILI (p < 0.05).
CONCLUSION: Tumor volume and threshold dose of irradiation damage correlated positively with MCP-1 and VEGF as well as NGF and CX3CL, respectively. Thus, investigation of biological mediators in blood samples from tumor patients may provide an appropriate tool to predict fRILI after interstitial HDR brachytherapy of liver metastases.

Introduction

Local-ablative therapies are emerging treatment options for primary and secondary liver malignancies such as hepatocellular carcinoma or colorectal cancer liver metastases (Vogel et al. 2021; Cutsem et al. 2016). In the last two decades, several techniques have been developed including thermal and irradiation-based ablation therapy, such as e.g. Yttrium90 radioembolization (Y90 RE), stereotactic-ablative body radiotherapy (SABR) and high dose rate (HDR) interstitial brachytherapy (iBT) (Garza-Ramos and Toskich 2021; Ricke and Wust 2011; Manchec et al. 2021; Romesser et al. 2021). A common limitation of the liver-directed radiotherapy is radiation-induced liver disease (RILD), often originating from sinusoidal obstruction syndrome (SOS), formerly known as hepatic veno-occlusive disease (VOD) (Fan and Crawford 2014; Bayraktar et al. 2017). Histopathological findings of RILD include deposits of extracellular matrix leading to a congestion of sinusoids and central veins with subsequent necrosis of hepatocytes (Lawrence et al. 1995; Fajardo and Colby 1980). In case that liver dysfunction induced by irradiation exceeds a tolerable limit, clinical symptoms of RILD such as weight gain and jaundice by hyperbilirubinemia occur (Sangro et al. 2008; Doi et al. 2017). RILD is stratified in two types: classic RILD and non-classic RILD. Patients with classic RILD usually exert fatigue, increased abdominal girth, abdominal pain, hepatomegaly and anicteric ascites approximately 2 weeks to 4 months after hepatic irradiation (Lawrence et al. 1995). The pathological hallmark is VOD. While the levels of alkaline phosphatase increase, the levels of transaminases and bilirubin remain normal (Liang et al. 2011). Non-classic RILD develops in patients who suffer from chronic hepatic diseases, e.g. cirrhosis and viral hepatitis (Cheng et al. 2004). Patients with non-classic RILD exert more dysregulated hepatic functions with jaundice and/or, in contrast the classic RILD, markedly increased levels of transaminases rather than alkaline phosphatase (Cheng et al. 2004; Pan et al. 2010; Kim and Jung 2017). Compared with loco-regional treatment by Y90 RE, SABR or iBT provide high-conformal irradiation of liver malignancies and offer the possibility to quantify focal radiation-induced liver injury (fRILI), which is visualized in hepatobiliary magnetic resonance imaging (MRI) by fusion with the 3D irradiation plan (Hass et al. 2019). This technique facilitates the correct assessment of the threshold isodose of hepatic irradiation tolerance. Regarding molecular mechanisms of hepatic dysfunction after irradiation only little is known. Although the knowledge derives from analyses of VOD/SOS after bone-marrow transplantation (BMT), this is not necessarily representative for conditions induced by SABR and iBT (Shulman and Hinterberger 1992).

Although our knowledge of fRILI pathogenesis has improved in the last years, the underlying molecular pathomechanisms remain unclear. The irradiation causes early DNA damage, oxidative stress and the production of reactive oxygen species, which lead to hepatocellular apoptosis and acute inflammatory responses in irradiated regions (Yahyapour et al. 2018; Baran et al. 2022; Robbins and Zhao 2004). During the past decade considerable developments of liquid biopsy technologies as noninvasive and real-time approached to represent the tumor burden and comprehensively reflect the genetic profile of HCC, e.g. molecular analysis of the tumor-derived fraction of the cell-free DNA (cfDNA) known as circulating tumor DNA (ctDNA), have emerged (Zhang et al. 2021; Roy et al. 2021; Bettegowda et al. 2014). Postoperative ctDNA sequencing has been shown of a great prognostic value in patients with liver cancer (Zhang et al. 2021; Ye et al. 2022). Thus, analysis of ctDNA has many potential applications, including noninvasive tumor genotyping, early detection of tumor recurrence but also detection of minimal residual disease following radiotherapy (Chaudhuri et al. 2015). Adjacent to ctDNA techniques, other biological structures including proteins could lead to personalized therapy based on their identification and alterations in their presence before and after treatment As such, several regulators of the repair responses to liver damage e.g. growth factors and cytokines, including tumor necrosis factor alpha (TNF-α), transforming growth factor beta (TGF-beta) and hedgehog have been associated with fRILI (Lee and Friedman 2011). TGF-beta has been identified to mediate sinusoidal obstruction by deposition of collagen matrix after radiation exposure (Pihusch et al. 2005). Cytokines are highly related to cell-to-cell communication reflecting the local or systemic inflammatory milieu (O’Shea et al. 2002). Since their systemic levels may correlate with organ injury, they can be applied as biomarkers for diseases. Furthermore, they play a key role in tumor escape mechanisms, which are responsible for metastatic outgrowth (Chen et al. 2020; Kany et al. 2019). Although pro-inflammatory cytokines may display anti-cancer activity, they are also responsible for tumor growth. Their local expression levels have been associated with liver integrity, but regarding their systematic levels, cytokine release during therapy or pathology itself can be used for fRILI quantification (Llovet et al. 2021; Conlon et al. 2019; Dranoff 2004). Local therapies are known to induce antigen and cytokine release, as reported in case of inhibitors of the vascular endothelium growth factor (VEGF) and tyrosine kinase that boost immunity and prime tumors for checkpoint inhibition. However, the knowledge about their precise role in fRILI is limited (Llovet et al. 2021; Herranz-Itúrbide et al. 2021). Thus, local therapies can provoke various biological effect caused by the different release pattern of humoral mediators including cytokines and growth factors, which on the one hand enhance the immune response or induce organ damage (e.g. fRILI), and promote tumor progression on the other hand (Llovet et al. 2021; Kim et al. 2013; Fornaro et al. 2014). While VEGF is already in clinical use, several other potential candidates such as nerve growth factor (NGF) have been suggested as markers of tumor progression and potential targets for future innovative therapies (Demir et al. 2016). Thus, identifying further clinical and biological factors that play an important role in the timely diagnose and fRILI progression will improve the understanding of the pathogenesis and potentially maximize the iBT treatment efficacy by reducing the fRILI. Therefore, investigating the correlation between systemic cytokine levels and patient’s outcome, may be useful in understanding their etiology in the clinical course of the pathology. As example, pre- and posttreatment “liquid biopsy” can distinguish between hepatocellular carcinoma patients with progression versus non-progression after radiotherapy (Reiss et al. 2021). With this background, it remains fundamental to understand the role of humoral mediators e.g. cytokines, as potential indicators and/or mediators of local radio sensitivity, radio resistance and distant abscopal effects (Powerski et al. 2020). This latter, radiation-induced effect on the neighboring (bystander) cells as well as fRILI are possibly related to the release of soluble mediators e.g. chemokines and growth factors (Bentzen 2006; Rodemann and Blaese 2007). The aim of this exploratory study was to identify circulating factors in patients undergoing iBT of liver metastases and correlate their levels with both volume and threshold doses of irradiated liver tissue, to potentially propose biomarkers for post-therapeutic fRILI.

Materials and methods

Patient cohort

Twenty-two patients (male n = 11, female n = 11, age 39–80 years) undergoing interstitial HDR brachytherapy of liver metastases (colorectal carcinoma n = 16, breast cancer n = 4, neuroendocrine tumor n = 1, lymphangiosarcoma n = 1) were prospectively included in this study. The study was conducted in accordance with the Declaration of Helsinki. Exclusion criteria were clinically ongoing infections and advanced liver cirrhosis with diminished uptake of hepatobiliary contrast agents in MRI. All included patients gave their written informed consent. Prospective data collection and analyses were approved by the local ethics committee. Inclusion criteria were: (i) irresectable liver metastases scheduled for iBT, (ii) no pre-existing liver cirrhosis, (iii) no contraindication to undergo hepatobiliary MRI, (iv) age 18–85 years, (v) ECOG performance status 0 or 1; (vi) no prior irradiation therapy of the same liver lobe.

Treatment by interstitial brachytherapy

The procedure of interstitial HDR brachytherapy was previously described (Ricke and Wust 2011). Briefly, the 6F catheter sheaths (Terumo Radifocus® Introducer II, Terumo Europe, Leuven, Belgium) and 6F irradiation catheters (afterloading catheter, Primed® medical GmbH, Halberstadt, Germany) were inserted into the liver metastases under computed tomography (CT) fluroscopy utilizing Seldinger’s technique and concomitant conscious sedation with Fentanyl and Midazolam. In complex or larger lesions, multiple catheters may be placed to achieve a sufficient dose distribution and to lower the detrimental radiation effect onto adjacent organs. Then, a single fraction irradiation with an Iridium 192 source was applied in afterloading technique according to a 3D treatment plan (Oncentra® Brachy, Elekta Instrument AB, Stockholm, Sweden). Ablative doses to the planning target volume (PTV) were selected according to the tumor entity (e.g. 25 Gy in colorectal cancer metastases, 15 Gy in metastases of breast cancer or neuroendocrine tumors) (Wieners et al. 2011; Ricke et al. 2010; Schippers et al. 2017). Finally, catheters were removed leaving gelatin sponge in the catheter path to minimize the risk of postinterventional bleeding events. The interstitial brachytherapy data set is described in Fig. 1.

Fig. 1

Brachytherapy of a colorectal cancer metastasis. A Pre-treatment hepatobiliary magnetic resonance imaging (MRI) depicting a large lesion in liver segment 7. B computed tomography-guided insertion of irradiation catheters into the target volume. C Follow up hepatobiliary MRI after 6 weeks with shrinkage of the metastasis and additional rim of focal radiation-induced liver injury (arrow heads)

Data acquisition and follow-up

Regular schedule of interstitial brachytherapy in the Department of Radiology and Nuclear Medicine (Otto-von-Guericke University, Magdeburg, Germany) includes baseline imaging with hepatobiliary MRI, a body CT as well as the laboratory evaluation of liver specific parameters and coagulation status one day prior to irradiation and three days after irradiation. Routine follow-up is assigned every 3 months including imaging and assessment of the standard laboratory parameters. Within this study, hepatobiliary MRI was additionally performed 6 weeks after iBT to detect the maximal extent of fRILI at this time, while the definition of fRILI followed Seidensticker et al. (2014). Blood samples for cytokine analyses were taken after 6 weeks correspondingly. Imaging data of all patients were entirely acquired in a 1.5 T MRI scanner (Achieva, Philips, Best, The Netherlands) including 3D axial THRIVE sequences (T1 high resolution isotropic volume excitation) 20 min after intravenous application of hepatocyte-specific contrast agent Gd-EOB-DTPA (Primovist, Bayer Healthcare, Leverkusen, Germany). MRI data sets where then semi-automatically fused with the CT-based 3D irradiation plan using a point-to-point registration. Afterward, the threshold of fRILI was correlated with the specific isodose. The method was described aby Seidensticker et al. (2014).

Laboratory analysis of cytokines

Blood samples were obtained one day prior to (baseline) and 6 weeks after irradiation in prechilled citrate tubes (BD Vacutainer, Becton Dickinson Diagnostics, Aalst, Belgium) and were kept on ice. Blood was centrifuged at 2000×g for 15 min at 4 °C. Then, plasma was stored at − 80 °C until the sample analyses. Blinded specimens were used for measurements of MCP-1, CX3CL1, VEGF, beta-NGF and TGF-beta by the laboratory of the Experimental Radiology, Department of Radiology and Nuclear Medicine at the University Hospital of the Otto-von-Guericke University Magdeburg using a highly specific commercially available using LEGENDplex Multi-Analyte Flow Assay Kit (BioLegend, San Diego, California) according to manufacturer’s instructions. The measurements and analyses were performed with the BD FACSCanto II™ (BD Heidelberg, Germany). Twenty-two patients were included, and ten healthy volunteers served as controls.

Statistical analysis

Statistical analyses were performed by IBM SPSS Statistics 22.0®. Besides descriptive statistics, changes of cytokine levels in healthy volunteer controls (HV) and tumor patients (TP) before interstitial brachytherapy (pre iBT) and after interstitial brachytherapy (post iBT) were analyzed by the non-parametric Kruskal–Wallis test. Furthermore, spearman correlation analysis was performed to explore correlations between the levels of the assessed humoral mediators and clinical characteristics of the treatment, including tumor and fRILI volumes as well as the isodose. Data are presented as mean ± standard deviation (SD) and in the results given as mean ± standard error of mean (SEM). Statistical significance was assumed for p < 0.05.

Results

Data description

Treatment characteristics

In 22 patients, a total of 30 lesions (1–3 per patient) were ablated by interstitial brachytherapy. Mean liver volume was 1416 ± 382 ml and mean clinical target volume (equalling tumor volume) was 28 ± 35 ml. A D100 of 21.62 ± 4.24 Gy and D90 of 30.46 ± 6.57 Gy was achieved. All patients underwent initial follow-up after 6 weeks including hepatobiliary MRI as well as laboratory evaluation of liver specific and inflammatory parameters. In 19 out of 22 patients (86%), imaging revealed a fRILI with a mean volume of 92 ± 97 ml. Image fusion with the 3D irradiation plan correlated fRILI to a mean hepatic tolerance dose of 14.2 ± 5.7 Gy. A summary of patient and treatment characteristics is given in Table 1.

Table 1

Patient and treatment characteristics (n = 22)

Variables	n (%) or mean ± standard deviation
Age	62.0 ± 10.9 years
Sex (m/w)	n = 11 (50%)/n = 11 (50%)
Tumor entities
Colorectal cancer	n = 16 (72%)
Breast cancer	n = 4 (18%)
Neuroendocrine tumor	n = 1 (5%)
Lymphangiosarcoma	n = 1 (5%)
Number of lesions	n = 30
Liver volume	1416 ± 382 mL
Clinical target volume (CTV)	28 ± 35 mL
Dose to CTV
D100	21.6 ± 4.2 Gy
D90	30.5 ± 6.6 Gy
Focal radiation-induced liver injury (fRILI)	n = 19 (86%)
Threshold isodose	14.2 ± 5.7 Gy
Volume	92 ± 97 mL

Humoral mediators in tumor patients pre and post iBT versus healthy volunteer controls

The data have shown a trend to decreased MCP-1 levels in TP both pre and post treatment compared to HV; however, no statistical changes were found between the groups (HV: 145.30 ± 33.17 vs. TP pre iBT: 94.30 ± 12.40 vs. TP post iBT: 88.21 ± 9.68 pg/mL; Fig. 2A). Similarly, the CX3CL1 plasma levels demonstrated no significant changes between the three groups (HV: 429.70 ± 165.30 vs. TP pre iBT: 65.84 ± 14.36 vs. TP post iBT: 68.78 ± 14.16 pg/mL; Fig. 2B). Interestingly, the VEGF levels have shown a significant decrease in TP pre and post iBT compared to HV (p < 0.05, HV: 118.60 ± 27.26 vs. TP pre iBT: 35.52 ± 3.68 vs. TP post iBT: 35.19 ± 3.81 pg/mL; Fig. 2C). The beta-NGF levels were significantly increased in both TP pre and post iBT compared to HV (p < 0.05, HV: 2.16 ± 0.41 vs. TP pre iBT: 7.08 ± 0.93 vs. TP post iBT: 6.20 ± 0.94 pg/mL; Fig. 2D).

Fig. 2

Humoral mediators in tumor patients (TP) pre and post interstitial high dose rate brachytherapy (iBT) versus healthy volunteer controls (HV). Systemic levels of humoral mediators A monocyte chemoattractant protein-1 (MCP-1), B chemokine (C-X3-C motif) ligand 1 (CX3CL1), C vascular endothelial growth factor and D beta-nerve growth factor (beta-NGF) are shown. Plasma levels of all humoral mediators were compared from pre iBT (baseline) and post iBT (after 6 weeks) in TP and HV (n = 10). Data are presented as mean ± standard error of the mean. *p < 0.05 between the indicated group

Correlation analyses of humoral mediators in tumor patients pre and post iBT with the tumor volume

Baseline circulating levels of MCP-1, VEGF and beta-NGF have shown significant positive correlations with the hepatic tumor volume (Fig. 3A MCP-1: r = 0.6640, p = 0.0008. Figure 3C VEGF: r = 0.5804, p = 0.0046, and Fig. 3D, beta-NGF: r = 0.4802, p = 0.0237, respectively). CX3CL1 did not correlate with the tumor volume (Fig. 3B). There were no significant correlations between MCP-1, VEGF, beta-NGF and CX3CL1 with the hepatic tumor volume index 6 weeks after interstitial brachytherapy (Fig. 4).

Fig. 3

Correlation analyses of humoral mediators in tumor patients (TP) before interstitial high dose rate brachytherapy (pre iBT) with the tumor volume. Spearman correlation analyses between the tumor volume with the systemic levels of A monocyte chemoattractant protein-1 (MCP-1), B chemokine (C-X3-C motif) ligand 1 (CX3CL1), C vascular endothelial growth factor or D beta-nerve growth factor (beta-NGF) pre iBT are shown

Fig. 4

Correlation analyses of humoral mediators in tumor patients (TP) 6 weeks after the interstitial high dose rate brachytherapy (iBT) with the tumor volume. Spearman correlation analyses between the tumor volume with the systemic levels of A monocyte chemoattractant protein-1 (MCP-1), B chemokine (C-X3-C motif) ligand 1 (CX3CL1), C vascular endothelial growth factor or (D) beta-nerve growth factor (beta-NGF) post iBT are shown

Correlation analyses of humoral mediators in tumor patients pre and post iBT with fRILI

Circulating levels of humoral mediators before the treatment did not correlate with fRILI after irradiation (Fig. 5). CX3CL1 and VEGF have shown significant positive correlations with fRILI after treatment (Fig. 5B, CX3CL1: r = 0.6246, p = 0.0333 and Fig. 5C, VEGF: r = 0.7793, p = 0.0024, respectively). MCP-1 and beta-NGF levels assessed after interstitial brachytherapy did not correlate with fRILI doses (Fig. 6).

Fig. 5

Correlation analyses of humoral mediators in tumor patients (TP) before interstitial high dose rate brachytherapy (pre iBT) with focal radiation-induced liver injury (fRILI). Spearman correlation analyses between fRILI with the systemic levels of A monocyte chemoattractant protein-1 (MCP-1), B chemokine (C-X3-C motif) ligand 1 (CX3CL1), C vascular endothelial growth factor or (D) beta-nerve growth factor (beta-NGF) pre iBT are shown

Fig. 6

Correlation analyses of humoral mediators in tumor patients (TP) 6 weeks after the interstitial high dose rate brachytherapy (iBT) with focal radiation-induced liver injury (fRILI). Spearman correlation analyses between fRILI with the systemic levels of A monocyte chemoattractant protein-1 (MCP-1), B chemokine (C-X3-C motif) ligand 1 (CX3CL1), C vascular endothelial growth factor or D beta-nerve growth factor (beta-NGF) post iBT are shown

Discussion

SABR and iBT are irradiation techniques for the therapy of primary and secondary liver malignancies that deliver high-conformal irradiation to target volumes while preserving surrounding healthy liver tissue. However, RILD is still a limiting factor for the therapy application in large tumors and/or multifocal disease and might lead to a dose reduction with subsequent risk of non-ablative treatment. Although the histopathology of RILD is well described and related to VOD or SOS, little is known about the underlying molecular mechanisms, which on the one hand might allow to identifying patients at risk, or on the other hand, may uncover potential targets for preventive medication (Lawrence et al. 1995; Fajardo and Colby 1980). Several findings identified TGF-beta to play a major role in sinusoidal obstruction by promoting proliferation of myofibroblasts, and thereby the deposition of extracellular matrix in hepatic central veins and sinusoids during RILD (Pihusch et al. 2005).

Although hepatobiliary MRI after iBT enables a precise determination of the irradiation damage, we did not find any relationship between circulating TGF-beta levels and either the volume or the threshold isodose of fRILI as the TGF-beta levels were below the detection limit. Thus, our results might suggest that the long-time storage of plasma samples may cause a TGF-beta degradation. Moreover, this could affect the levels of cytokines, since activated platelet fraction, which are the major enhancer of TGF-beta, is reduced using the citrate tube (Assoian and Sporn 1986). This fact should be critically considered since we cannot avoid the inaccurate reflection of circulating levels of TGF-beta ex vivo (Kropf et al. 1997). Interestingly, the lack of TGF-beta in plasma samples might be related to radiation effects directly and locally affecting the liver tissue also. Consequently, the local ablative therapy may have impacted also the local TGF-beta levels, where it exerts its autocrine and paracrine effects (Braun et al. 1988). Importantly, TGF-beta is involved also in liver repair processes (Yang et al. 2016). It has been shown in preclinical studies that the active form of TGF-beta appears to be rapidly sequestered and degraded by the liver (Wakefield et al. 1990). However, we were unable to discern between its active and latent form, and thus, further investigations are required.

In contrast to TGF-beta, significant correlations were found between both MCP-1 and VEGF with the hepatic tumor volume. Furthermore, the data have shown that plasma levels of CXCL1 and VEGF positively correlated with the threshold isodose fRILI. It is known that the tumor growth, metastasis, and veno-occlusive disease (SOD), may depend on the degree of tissue vascularization mediated by VEGF (Kong et al. 2021; Iguchi et al. 2001). However, there is no much information regarding the response to iBT. In our study, we observed that those patients had decreased circulating levels of CX3CL1, VEGF and MCP-1 as compared to healthy donors. Moreover, the correlation between those factors and the tumor volume suggests that the lower baseline circulating levels may indicate the lower tumor volume. A reasonable explanation may be provided by the fact that during the iBT induced tumor cyto-reduction, the bioavailability of VEGF may also have been reduced. This is underlined by the ability of free VEGF to bind to vacant sites, either on its receptors (VEGFR) or binding proteins, which may arise mostly from the tissue damaged by the therapy (Tsai et al. 2021). Taken together, baseline levels of VEGF correlated with the tumor volume, while the post therapy VEGF levels correlated with the fRILI, indicating at a potentially important role of VEGF as a biomarker for therapy outcome.

The current literature indicates that either healthy hepatic or malignant cells are sources of MCP-1 (Wyler et al. 2016; Czaja et al. 1994). Therefore, the positive correlation between MCP-1 and tumor volume reduction may be associated with the iBT, which might consequentially also reduce the circulating protein levels. Similar scenario may also be supported by other chemotactic cytokines such as CX3CL1 (Bazan et al. 1997). Interestingly, CX3CL1 and its receptor can be upregulated in injured human liver (Efsen et al. 2002). Moreover, in hepatocellular carcinoma patients, a good prognosis and low recurrence rates are related to high expression levels of CX3CL1/CX3R1 (Matsubara et al. 2007).

The neuronal tracking can be mediated by various molecules including chemokines and NGF secreted by the tumor microenvironment. Conversely, the latter was found to be associated with wound healing, combined with radiation in vitro (Shi et al. 2003). Therefore, the iBT may provoke specific effects of NGF induction with the aim to protect from radiation and to induce local tissues regeneration. Taken together, further investigations are needed to understand the association between humoral mediators and iBT/fRILI.

Limitations

Our study has several limitations. One major critical point is the time frame between the clinical therapy application and the laboratory measurements. The stability of cytokines may be affected by long-term storage (de Jager et al. 2009a, b). Therefore, the biomaterial has been processed as soon as possible, and stored at − 80 °C avoiding freezing and thawing processes. This issue needs to be considered in the overall assessment. Moreover, the inclusion criteria did not account whether the tumor patients had natural tolerance or sensitivity to radiation. Therefore, correlating the therapy outcome with biomarkers of liver radiation sensitivity would be advantageous in future. Another limitation is the analysis of the patient cohort with metastases deriving from different primary tumors (e.g. colorectal carcinoma and breast cancer) which might pose a bias as different tumor entities may cause specific microenvironments in involved tissues. This limits the conclusions about tumor volume and circulating cytokines and should be addressed in future studies. Also more frequent time course for the assessment of biomarkers after therapy should be included in a future study. Same applies for the small number of patients included in this study which prohibits detailed statistical tests and/or subgroup analyses. Furthermore, detection of fRILI in hepatobiliary MRI and correlation of with a threshold isodose by 3D image fusion was only reported in a few studies and relies on the diminished uptake of the hepatocyte-specific contrast agent after affecting a single transporting polypeptide (Organic Anion Transporting Polypeptide 1, OATP1). Thus, not all aspects of liver cell functions are reflected by this imaging technique.

Conclusions

This exploratory study describes several humoral mediators to correlate with the focal hepatic dysfunction after interstitial brachytherapy of liver metastases. Positive correlations of tumor volume and threshold dose of irradiation damage were found for MCP-1 and VEGF as well as NGF and CX3CL, respectively. Therefore, the results from this study suggest that the investigation of humoral mediators in blood samples as “liquid biopsy” from patients may be an appropriate tool to predict side effects such as fRILI after interstitial HDR brachytherapy of liver metastases. However, as a perspective, the patient cohort should be expanded and adapted to homogeneous criteria. We likewise require examining tissue materials from tumor biopsies to characterize the status of the pathology.