Vascular response patterns to targeted therapies in murine breast cancer models with divergent degrees of malignancy
Breast Cancer Research volume 25, Article number: 56 (2023)
Response assessment of targeted cancer therapies is becoming increasingly challenging, as it is not adequately assessable with conventional morphological and volumetric analyses of tumor lesions. The tumor microenvironment is particularly constituted by tumor vasculature which is altered by various targeted therapies. The aim of this study was to noninvasively assess changes in tumor perfusion and vessel permeability after targeted therapy in murine models of breast cancer with divergent degrees of malignancy.
Low malignant 67NR or highly malignant 4T1 tumor-bearing mice were treated with either the multi-kinase inhibitor sorafenib or immune checkpoint inhibitors (ICI, combination of anti-PD1 and anti-CTLA4). Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) with i.v. injection of albumin-binding gadofosveset was conducted on a 9.4 T small animal MRI. Ex vivo validation of MRI results was achieved by transmission electron microscopy, immunohistochemistry and laser ablation-inductively coupled plasma-mass spectrometry.
Therapy-induced changes in tumor vasculature differed between low and highly malignant tumors. Sorafenib treatment led to decreased tumor perfusion and endothelial permeability in low malignant 67NR tumors. In contrast, highly malignant 4T1 tumors demonstrated characteristics of a transient window of vascular normalization with an increase in tumor perfusion and permeability early after therapy initiation, followed by decreased perfusion and permeability parameters. In the low malignant 67NR model, ICI treatment also mediated vessel-stabilizing effects with decreased tumor perfusion and permeability, while ICI-treated 4T1 tumors exhibited increasing tumor perfusion with excessive vascular leakage.
DCE-MRI enables noninvasive assessment of early changes in tumor vasculature after targeted therapies, revealing different response patterns between tumors with divergent degrees of malignancy. DCE-derived tumor perfusion and permeability parameters may serve as vascular biomarkers that allow for repetitive examination of response to antiangiogenic treatment or immunotherapy.
Since angiogenesis has been identified as one hallmark of cancer , the role of tumor vasculature in tumorigenesis is of prime interest. Angiogenesis is defined as the formation of new blood vessels from a pre-existing vascular bed, driven by predominance of proangiogenic factors and corresponding signaling pathways [2, 3]. In contrast to physiological blood vessel formation, which is an organized and controlled process, rapid tumor growth causes angiogenesis in tumor tissue to be dysregulated, leading to immature tumor vasculature  with several structural and functional abnormalities. Tumor-associated blood vessels are characterized by variable diameters and distorted shapes [3,4,5]. The endothelial layer is often tortuous and disrupted, leading to enhanced permeability with increased interstitial pressure [2, 5]. While these characteristics promote tumor dissemination to distant organs, they also lead to treatment limitations by restricted drug delivery into the tumor lesion and consequently reduced therapeutic effectiveness [5,6,7].
Many targeted cancer therapies closely interact with the tumor vasculature. To specifically target the aberrant properties of tumor vasculature itself, a variety of antiangiogenic therapies have been developed, mainly addressing the vascular endothelial growth factor pathway as one of the most significant proangiogenic factors . By restoring the balance between pro- and antiangiogenic pathways, antiangiogenic therapies aim for structural stabilization of the blood vessel wall, mediated by establishment of inter-endothelial contacts, re-development of regularly structured basement membranes and association and integration of vascular mural cells [9, 10]. This structural stabilization with subsequently reduced vascular permeability contributes to vascular normalization, the overall process of remodeling tumor blood vessels to restore their structure and function, and thereby leading to improved blood supply and reduction in intratumoral hypoxia .
In addition to antiangiogenic therapies, other classes of targeted therapies, especially immunotherapies, are also closely coupled to the vascular system . Since tumor blood vessels constitute a component of the tumor microenvironment (TME), there is a complex bidirectional interplay between tumor vasculature and immune cells. On the one hand, characteristics of tumor vasculature, e.g., endothelial permeability and expression of adhesion molecules, influence trafficking and migration of immune cells into the tumor lesion [13, 14]. On the other hand, activated immune cells within the TME, especially tumor-associated macrophages, neutrophils and T-cells, can produce additional angiogenic factors with subsequent effects on tumor perfusion and permeability .
Even though these targeted therapies are approved for a wide range of solid tumors [5, 16], including breast cancer , their effectiveness is limited to only a subgroup of patients while others show innate or acquired resistances [18, 19]. This variable response is presumably based on pronounced intra- and intertumoral heterogeneity of the TME. In this context, varying distribution in tumor vasculature within a single tumor lesion itself or between different cancer patients might lead to variable therapeutic drug delivery [20, 21]. Thus, as there is a significant number of patients who do not benefit from targeted therapies, early and repeated examination of therapy response is needed.
A noninvasive imaging approach to investigate the characteristics of tumor vasculature, a main component of the TME that is altered by targeted therapies, is constituted by dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI). By obtaining fast MRI sequences before, during and after rapid intravenous administration of gadolinium-based contrast agents and subsequent pharmacokinetic modeling, DCE-MRI allows for quantification of tumor perfusion and permeability . This MRI technique has already been applied to assess response to antiangiogenic treatment both in preclinical [23, 24] and clinical settings .
Here, DCE-MRI was performed to compare the vascular response patterns between low malignant 67NR and highly malignant 4T1 murine breast cancer models to different classes of targeted therapies. First, the sensitivity of the applied MR imaging approach to dynamically capture specific modulation of tumor vasculature characteristics was tested. To this end, dedicated manipulation of vasculature features was conducted by treating tumors with angiopoietin-1, which induces stabilization of blood vessel walls . Then, treatment with targeted therapies, including either a tyrosine kinase inhibitor or immune checkpoint inhibitors (ICI), was performed. While sorafenib is a multi-kinase inhibitor directly influencing tumor vasculature by inhibiting angiogenesis , immune checkpoint inhibitors primarily address immune cells and tumor-promoting inflammation with secondary effects on tumor vasculature . MRI results were confirmed by ex vivo analysis of tumor vasculature, including histology, immunohistochemistry and transmission electron microscopy.
Materials and methods
All animal experiments were approved by the responsible authorities (Landesamt für Natur, Umwelt und Verbraucherschutz NRW, Protocol No. 81-02.04.2018.A010). A total of 111 female BALB/c mice (Charles River Laboratories, Sulzberg, Germany) were used at the age of eight to twelve weeks and were kept under a 12 h light–dark-cycle, provided with food and water ad libitum. Two syngeneic murine breast cancer models have been used: highly malignant 4T1 tumors, metastasizing in regional lymph nodes and distant organs, primarily the lung, liver and bones, and 67NR tumors that grow noninvasively and do not develop metastases [29, 30]. Primarily based on their different growth kinetics, both tumor models contain significant differences in tumor vasculature. Due to rapid tumor growth, highly malignant 4T1 tumors show distorted blood vessels with a thin and disrupted endothelial layer. In contrast, slower growth of low malignant 67NR tumors enables intratumoral neoangiogenesis to proceed in a more controlled manner, resulting in increased maturity and integrity of intratumoral blood vessels .
4T1 and 67NR cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM, Thermo Fisher Scientific, Waltham, Massachusetts, USA), supplemented with 10% fetal bovine serum in standard conditions (37 °C, 5% CO2).
For tumor implantation, 106 4T1 or 67NR cells, resuspended in 25 µL cell culture medium, were implanted orthotopically into the lower left mammary fat pad. For daily scoring of the mice, tumor sizes were measured daily using a digital caliper.
Angiopoietin-1-induced vessel stabilization
To test whether the MR imaging approach is sensitive enough to detect targeted manipulation of tumor vasculature, 4T1 tumor-bearing mice were treated with angiopoieitin-1 that is known to promote vessel stabilization . Starting the day after tumor implantation, angiopoietin-1 treatment was administered by tail vein injection (5 µg/kg) every day . MR imaging was conducted on day six after tumor implantation. Afterward, mice were killed and tumors were prepared for ex vivo analysis.
Treatment initiation of both targeted therapies (sorafenib or immune checkpoint inhibitors) was started on day three after tumor implantation. Sorafenib treatment (Selleck Chemicals, Houston, Texas, USA) was given as a daily oral gavage (30 mg/kg) . Besides a weight loss of up to 5% during therapy, mice tolerated treatment well. As combination treatment enhances efficient anticancer activity, especially in 4T1 tumors , immune checkpoint blockade was induced by a combination of anti-CTLA4 and anti-PD1 inhibition (BioXCell, Lebanon, USA) via i.p. injection (10 mg/kg of each drug) every second day . Mice tolerated treatment well without any observable side effects. Data were compared to untreated control groups that received no therapy after tumor implantation. MR imaging was conducted on two sets of animals on either day six or day nine, and after MRI scans, mice were killed and tumors were prepared for ex vivo analysis.
MR imaging and analysis
MR imaging was performed on a 9.4 T small animal BioSpec system (Bruker BioSpin GmbH, Ettlingen, Germany), equipped with a 1H quadrature volume resonator for signal excitation and a 10 mm surface coil for signal reception (Bruker), using ParaVision 6.0.1. Mice were anesthetized with 1.5% isoflurane in 1.5 L/min of oxygen and compressed air (1:4) under continuous respiratory and temperature monitoring. Mice were placed in supine position and for reduction in susceptibility artifacts, the area between tumor and surface coil was covered in alginate (Johannes Weithas, Lütjenburg, Germany). For anatomical information, a T2-weighted multi-slice RARE (Rapid Acquisition with Relaxation Enhancement) sequence (TR = 2500 ms, TE = 11 ms, effective TE = 55 ms, 2 averages, 1 mm slice thickness, 9 slices, 20 × 20 mm2 FOV, 256 × 256 matrix, RARE factor = 12, acquisition time = 1:45 min:s) was acquired. The slice with the greatest mean tumor diameter was chosen for the subsequent single-slice sequences. T1 mapping was performed before and after injection of the contrast agent using a single-slice RARE sequence with variable TR (RARE-VTR, TR = 7500, 5000, 3000, 1500, 800, 400, 311, 123 ms, TE = 50 ms, 1 average, 1 mm slice thickness, 18 × 15 mm2 FOV, 108 × 96 matrix, acquisition time = 14:39 min:s) with consecutive quantification of T1 relaxation times using the image sequence analysis tool within ParaVision. For DCE imaging, the blood pool contrast agent gadofosveset trisodium (Lantheus Medical Imaging, North Billerica, Massachusetts, USA) was injected via a tail vein catheter (Klinika Medical GmbH, Usingen, Germany) with a concentration of 0.6 mmol/kg, using a perfusion pump (World Precision Instruments, Sarasota, Florida, USA) at a rate of 240 µL/min. The injection was initiated one minute after starting a FLASH (Fast Low Angle Shot) scan (TR = 24.6 ms, TE = 1.5 ms, 15° flip angle, 1 average, 610 repetitions, 18 × 15 mm2 FOV, 96 × 96 matrix, acquisition time = 20:00 min:s). As gadofosveset is a macromolecular contrast agent that binds reversibly to serum albumin and can thus only extravasate in case of enhanced vascular leakage, it enables not only assessment of tumor perfusion but also of endothelial permeability. Overall acquisition time of the entire MRI protocol was 51:03 min:s.
T2-weighted images were used to evaluate macroscopic changes after therapy and enabled assessment of tumor volume, analyzed using 3D Slicer (version 4.11.2021, https://www.slicer.org/). For subsequent analysis of tumor perfusion and permeability parameters, regions of interest (ROIs) were drawn around the viable tumor periphery with exclusion of necrotic tumor areas, using T2-weighted images as guidance. Three representative, non-overlapping ROIs were placed in the viable tumor periphery, and their mean yielded the average value per animal. To eliminate interindividual differences in pre-contrast T1 relaxation times, ROIs were copied from pre- to post-contrast T1 maps and delta of T1 values (∆T1) before and after contrast agent injection were calculated. Dynamic assessment of contrast enhancement after gadofosveset injection was used to derive the parameters area under the curve (AUC), maximum slope (slopemax) and volume transfer constant Ktrans. While the area under the signal intensity curve and its maximum slope describe characteristics of tumor perfusion, the transfer constant Ktrans assesses capillary permeability . Furthermore, plasma volume fraction vp, a marker for vascularization, was calculated. Calculation of these parameters was performed with an in-house developed software based on the Pk modeling extension for 3D Slicer (https://github.com/millerjv/PkModeling). First, the acquired pre-contrast T1 map was resampled to the matrix and pixel size of the DCE-MRI using a linear interpolator. Subsequently, pharmacokinetic modeling of the DCE measurement was performed using a three-parameter Tofts model (also referred to as extended Tofts model)  and the parameters (AUC, slopemax and Ktrans) were calculated pixel by pixel. A population-based arterial input function was used . The longitudinal relaxivity of the applied contrast agent was set to 0.0045 ms−1 mM−1, determined by extrapolation from literature values , and the hematocrit value was assumed to be 0.45.
Histology and immunohistochemistry
For histological analysis, tumors were paraffin-embedded and sectioned in 5 µm tissue slices using a rotary microtome (Leica Microsystems, Nussloch, Germany). Hematoxylin and eosin (H&E) staining was performed according to standard protocols and enabled identification of intratumoral hemorrhage and necrosis. Immunohistochemical staining of CD31, a marker for vessel density and angiogenesis , was performed as previously described . In brief, it was started with dewaxing and rehydration, followed by incubation in 3% H2O2 for 10 min. Afterward, sections were incubated with primary antibody CD31 (ab124432, Abcam, Cambridge, UK) for one hour at 4 °C in a dilution of 1:200 in blocking buffer, followed by horseradish peroxidase/diaminobenzidine detection. To analyze changes of the intratumoral immune cell infiltrate after immune checkpoint blockade, further immunohistochemical staining of CD3 (T-cells), F4/80 (macrophages), Ly6G (neutrophils), CD49b (natural killer cells) and CD45R (B-cells) was performed (Additional file 1: Meth. 1). For each group, n = 3 samples have been analyzed qualitatively.
Transmission electron microscopy
Transmission electron microscopy enabled assessment of endothelial junction integrity and intactness of endothelial layers and basal laminae. Tumors were fixed in 2% (v/v) formaldehyde and 2.5% (v/v) glutaraldehyde in 100 mM cacodylate buffer. After washing in PBS, post-fixation was conducted with 0.5% (v/v) osmium tetroxide and 1% (w/v) potassium hexacyanoferrate (III) in 0.1 M cacodylate buffer. After washing in distilled water and dehydration in an ascending ethanol series, tumors were incubated in propylenoxide twice for 15 min each. Next, small tissue pieces (1–3 mm3) of the tumor border were embedded in Epon (Sigma-Aldrich, St. Louis, MO, USA) using flat embedding molds. 80 nm sections were sectioned using a ultramicrotome, collected on copper grids and negatively stained with 2% uranyl acetate for 15 min . Electron micrographs were taken at a Philips EM-410 electron microscope (Philips, Amsterdam, The Netherlands). For each group, two samples have been analyzed qualitatively.
Shapiro–Wilk test was used to check for normal distribution of data (α = 0.05). Depending on p-values, treatment groups were compared to untreated control groups using a two-tailed t-test or Mann–Whitney U test. Statistics were performed using SAS (version 9.4, Copyright SAS Institute Inc., Cary, NC, USA), and graphs were created using GraphPad Prism (version 9.2.0, GraphPad Software Inc., San Diego, California, USA). Group sizes were as followed: n = 12 for 67NR 6d untreated, n = 8 for 67NR 6d sorafenib, n = 8 for 67NR 6d ICI, n = 10 for 67NR 9d untreated, n = 8 for 67NR 9d sorafenib, n = 7 for 67NR 9d ICI, n = 12 for 4T1 6d untreated, n = 6 for 4T1 6d angiopoietin-1, n = 8 for 4T1 6d sorafenib, n = 8 for 4T1 6d ICI, n = 9 for 4T1 9d untreated, n = 8 for 4T1 9d sorafenib, n = 7 for 4T1 9d ICI. Descriptive statistics and further details for analysis of significance are provided in Additional file 1: Tables S1–S4.
Angiopoietin-1-induced vessel stabilization
To evaluate if the MRI protocol allows accurate assessment of vessel stabilization, 4T1 tumor-bearing mice were treated with angiopoietin-1.
Angiopoietin-1 treatment had no significant influence on tumor volume compared to untreated control tumors on day six. However, sequential T1-weighted images and corresponding contrast enhancement curves of angiopoietin-1 treated tumors showed a much more pronounced wash-out of the contrast agent compared to untreated controls (Fig. 1a). All calculated tumor perfusion and permeability parameters, including area under the curve (AUC), maximum slope (slopemax) and Ktrans, decreased significantly after treatment. In line, angiopoietin-1 treatment also reduced the intratumoral retention of gadofosveset, assessed by decreasing ∆T1 compared to untreated controls (Fig. 1b, c). Furthermore, plasma volume fraction vp, a marker for overall vascularization, decreased significantly after treatment of 4T1 tumors (Additional file 1: Fig. S1). Ex vivo analysis of tumor vasculature confirmed the in vivo imaging results qualitatively: While H&E staining showed intratumoral hemorrhage in untreated tumors, angiopoietin-1 treated tumors exhibited almost no intratumoral hemorrhage and the expression of CD31 was reduced. In addition, ultrastructural analysis of blood vessels using transmission electron microscopy confirmed improved integrity of the endothelial layer after angiopoietin-1 treatment: Compared to the thin and disrupted endothelium of untreated tumors, the angiopoietin-1 treated endothelium was overall thicker, with almost no endothelial gaps and an intact basal lamina (Fig. 1d).
Impact of sorafenib treatment on tumor growth and vasculature
To assess the changes in tumor vasculature after antiangiogenic treatment, 4T1 or 67NR tumor-bearing mice were treated with the multi-kinase inhibitor sorafenib. Compared to untreated controls, sorafenib treatment resulted in a significantly reduced volume of 67NR tumors on both day six and day nine (Fig. 2a). However, sorafenib treatment had no significant effect on the volume of 4T1 tumors at any observed time point (Fig. 2b), with tumors showing pronounced hemorrhage in the tumor center.
Anatomical T2-weighted images revealed small areas of necrosis in sorafenib-treated 67NR tumors compared to homogenous untreated controls. Sorafenib treatment of 67NR tumors resulted in significantly decreased tumor perfusion (AUC, slopemax) and permeability (Ktrans) parameters with reduced retention of intratumoral gadofosveset, assessed by lower ∆T1 compared to untreated control group (Fig. 3a, b). In addition, vp was reduced after treatment (Additional file 1: Fig. S1). In line with T2-weighted images, H&E staining of sorafenib-treated 67NR tumors showed small areas of necrosis, increasing from day six to day nine, which were not visible in untreated 67NR tumors. Again, qualitative ex vivo analysis of tumor vasculature validated the in vivo imaging results of reduced tumor perfusion and permeability. Compared to high CD31 expression of untreated 67NR tumors, CD31 expression decreased after treatment (Fig. 3c). Transmission electron microscopy of sorafenib-treated 67NR tumor vessels showed intact endothelial layers with a continuous basal lamina on both time points (Fig. 3d). In addition, in a separate set of experiments, mean gadolinium concentrations within the ex vivo tumor slices were quantified using laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS, Additional file 1: Meth. 2). The mean gadolinium concentrations showed a significant negative correlation with T1 relaxation times after gadofosveset injection (Additional file 1: Fig. S2). Since gadofosveset is a macromolecular, albumin-binding contrast agent, gadolinium extravasation and subsequent detection in LA-ICP-MS analysis occurs in tissues with increased endothelial permeability, confirming that the applied MRI approach can accurately assess vascular permeability.
In sorafenib-treated 4T1 tumors, T2-weighted images revealed increasing intratumoral hemorrhage and necrosis, specifically arranged as a round structure in the tumor center. Tumor vascularization marker vp was reduced on both day six and day nine compared to control groups (Additional file 1: Fig. S1). However, tumors exhibited significantly increased tumor perfusion and permeability parameters early on day six (about 1.5-fold), with significantly elevated ∆T1 compared to untreated control group. From day six to nine, all DCE-derived perfusion and permeability parameters decreased. On day nine, tumor perfusion parameters (AUC, slopemax) were still slightly elevated compared to control groups, while tumor permeability parameter Ktrans as well as ∆T1 showed no significant differences as compared to control groups anymore (Fig. 4a, b). The specific tumor architecture of sorafenib-treated 4T1 tumors, including the round area of necrosis and hemorrhage, was also visible in H&E staining and LA-ICP-MS (Additional file 1: Fig. S3). In contrast to DCE-derived perfusion parameters being elevated compared to control groups on day six, CD31 expression of 4T1 tumors was reduced after treatment on both day six and day nine, as assessed qualitatively (Fig. 4c). While in vivo DCE-MRI assessed increasing permeability on day six, sorafenib-treated 4T1 tumors showed less interruptions of the endothelial layer compared to control groups on both time points (Fig. 4d).
Impact of immune checkpoint blockade on tumor growth and vasculature
To further evaluate the effects of immunotherapy on tumor vasculature, 4T1 or 67NR tumor-bearing mice were treated with a combination of anti-PD1 and anti-CTLA4 immune checkpoint inhibitors. While size of ICI-treated 67NR tumors was significantly reduced on day six and day nine (Fig. 2a), ICI combination therapy led to a significant increase in size of 4T1 tumors as compared to control groups (Fig. 2b). Both tumor models exhibited substantial infiltration of T-cells after ICI therapy, as assessed by immunohistochemistry (Additional file 1: Fig. S4).
ICI-treated 67NR tumors showed only little intratumoral necrosis compared to untreated controls in T2-weighted imaging. Similar to sorafenib treatment, all DCE-derived perfusion and permeability parameters decreased, paralleled by lower intratumoral gadolinium retention (Fig. 5a, b) and reduction in vascularization marker vp (Additional file 1: Fig. S1). In line with in vivo imaging results, H&E staining of ICI-treated 67NR tumors showed small areas of necrosis with little intratumoral hemorrhage, while qualitatively assessed CD31 expression of 67NR tumors decreased after treatment (Fig. 5c). Transmission electron microscopy revealed intact and thick endothelium with a continuous surrounding basal lamina, again similar to the characteristics after antiangiogenic sorafenib treatment (Fig. 5d).
In contrast, ICI-treated 4T1 tumors exhibited elevated DCE-derived perfusion and permeability parameters with elevated gadolinium retention after treatment (Fig. 6a, b), accompanied by decreased vascularization (Additional file 1: Fig. S1). Even though qualitatively assessed CD31 expression decreased, H&E staining showed excessive intratumoral hemorrhage and increasing tumor necrosis. Transmission electron micrographs confirmed enhanced vascular leakage, by capturing very thin and tortuous endothelial layers with markedly increasing interruptions and endothelial cells being detached from the basal lamina after treatment (Fig. 6d).
The aim of this study was to noninvasively assess and compare the vascular response patterns to two different targeted therapies (sorafenib or immune checkpoint inhibitors) between highly malignant 4T1 and low malignant 67NR tumors using dynamic contrast-enhanced MRI.
Conventional chemotherapy targets proliferating cells unselectively, and its successful antitumor effects can be reliably detected by decreasing lesion size. In comparison, targeted therapies alter specific characteristics of the tumor microenvironment, including tumor perfusion and permeability. These changes frequently occur without an immediate effect on lesion size; thus, response assessment becomes more challenging . Indeed, tumor volumetry, performed using T2-weighted images, was not suitable to reliably indicate a therapy response in this study either. Although low malignant 67NR tumors demonstrated a reduction in tumor size after both targeted therapies, this was not the case for the highly malignant 4T1 model, which showed no significant changes in tumor size after sorafenib treatment and even increasing size after immune checkpoint blockade, compared to control groups.
Instead, this study clearly demonstrates that DCE-MRI with an albumin-binding contrast agent is able to capture early changes in tumor vasculature after both antiangiogenic treatment and immune checkpoint blockade, while revealing different response patterns between tumors with varying grade of malignancy. Low malignant 67NR tumors showed decreasing tumor perfusion and permeability after sorafenib treatment, while 4T1 tumors were characterized by an initial increase in DCE-derived parameters on day six, with a subsequent decrease to day nine. While ICI therapy also led to decreased tumor perfusion and permeability in 67NR tumors, parameters strongly increased in ICI-treated 4T1 tumors. The suitability of the DCE-MRI approach with albumin-binding gadofosveset to evaluate both tumor perfusion and vessel permeability was indicated by reduced perfusion and permeability parameters after angiopoietin-1 treatment, which is known for its vessel-stabilizing effects . Furthermore, the vascular changes implied by noninvasive DCE-MRI were confirmed by ex vivo structural analysis of tumor blood vessels using histology and transmission electron microscopy, providing evidence that the calculated DCE-MRI parameters truly reflect tumor vasculature characteristics.
In low malignant 67NR tumors, sorafenib treatment was able to inhibit angiogenesis and stabilize intratumoral blood vessels. In line with immunohistochemistry showing a reduced CD31 expression and electron microscopy revealing intact and continuous endothelial layers of tumor capillaries, DCE imaging parameters confirmed decreasing tumor perfusion and vessel permeability after treatment. These vascular changes after antiangiogenic therapy have already been recorded in other preclinical and clinical studies [24, 25, 44, 45]. They are in accordance with the initial rationale that sorafenib inhibits blood vessel formation, which in turn will cause profound vascular regression and perfusion, leading to tumors ‘starving to death’ [46, 47]. Using DCE-MRI, these therapy-induced effects can be monitored noninvasively and longitudinally over time. Compared to 67NR tumors, the effects of sorafenib treatment on vasculature of highly malignant 4T1 tumors differed strongly. On day three after treatment initiation, sorafenib-treated 4T1 tumors exhibited increasing tumor perfusion and permeability compared to control groups. Three days later, perfusion and permeability decreased to values almost comparable with control groups again. These early vascular changes after antiangiogenic treatment indicate a stabilization of the vessel barrier . By reverting tumor vasculature from a grossly abnormal structure toward a more physiological state, antiangiogenic therapies can create a transient window of enhanced tumor perfusion and oxygenation for a few days, followed by subsequent reduction in perfusion [47,48,49]. While we were able to capture vascular normalization on treatment day three, followed by decreasing parameters on day six, we are not able to draw conclusions about the potential long-term changes beyond this period. Increases of Ktrans values in the 4T1 model have already been shown after bevacizumab, another antiangiogenic treatment . In recent studies, the vascular normalization window after antiangiogenic therapy and its subsequent enhanced tumor perfusion was used to improve the drug delivery of other therapies, e.g., chemotherapies, into the tumor lesion . However, as it is only a transient stage of improved tumor perfusion, it is important to be able to precisely detect these vascular changes, which is possible using DCE-MRI as presented. In this study, the vascular normalization phenomenon was only visible in the highly malignant 4T1 model, as only these tumors were initially characterized by a tortuous and abnormal tumor vasculature, while the vascular system of low malignant 67NR tumors was generally more intact at baseline before treatment .
The response patterns of both tumor models did not only differ after sorafenib treatment but also after immune checkpoint blockade. In the low malignant 67NR model, ICI treatment also mediated stabilizing effects on intratumoral blood vessels with decreasing tumor perfusion and permeability, in line with reduction in CD31 expression and microstructurally intact blood vessels. These therapeutic effects after ICI treatment have also been observed in other preclinical studies [51, 52], applied to different murine breast cancer models and a fibrosarcoma model, as well as first clinical studies in patients with metastatic melanoma . They may be attributed to vascular remodeling after cytokine secretion of tumor-infiltrating immune cells, such as TNFα . Contrary to that, highly malignant 4T1 tumors exhibited enhanced DCE-derived perfusion and permeability parameters with overall slightly increasing CD31 expression. In this context, therapy-induced hyperperfusion and enhanced endothelial permeability were recently found to be correlated with successful activation of antitumor T-cell immunity [15, 53], which again underlines that the response patterns after targeted therapies are marked by huge intra- and intertumoral heterogeneity.
The divergent response patterns to both targeted therapies between the two tumor models used within this study might be related to differences in a variety of tumor characteristics determined by their respective malignant potential. While low malignant 67NR tumors demonstrated overall maturity and integrity of their vascular system, rapid tumor growth of highly malignant 4T1 tumors led to abnormal and tortuous intratumoral blood vessels, hindering drug delivery into the tumor lesion. The dominating intratumoral immune cell populations and their secretion of different angiogenic factors as well as PD-L1 expression of the tumors might additionally contribute to the divergent vascular response patterns [15, 54]. Furthermore, different response patterns might be attributable to differences in mutational status, epigenetic modification and cancer cell plasticity between tumors with varying malignant potential .
Despite the presented accuracy for detecting different vascular response patterns between both tumor models, the presented study is limited in some aspects of methodology. Due to the high temporal resolution, only single-slice MR imaging was performed, which led to loss of information on slice-adjacent tumor areas. However, multi-slice T2-weighted imaging was used to place DCE slices across the largest tumor diameter, suggesting representative tumor areas and reducing partial volume effects. ROIs were placed around the viable tumor border to exclude tumor necrosis, where perfusion parameters do not lead to meaningful conclusions about tumor vasculature. For future studies, reduction in temporal in favor of spatial resolution could be performed to further assess intratumoral heterogeneity, a main contributor to malignant tumor progression, metastasis formation and therapy resistance . For precise analysis of not only tumor perfusion but also permeability, numerous studies have demonstrated that macromolecular contrast agents (MMCA) are more favorable than small molecular contrast agents (SMCA) [57,58,59]. While SMCAs have an unselective extracellular extravasation profile and already extravasate in physiological tissue, MMCAs follow an intravascular distribution profile and are only able to extravasate when the endothelial permeability is pathologically enhanced, as known in tumor vasculature. Due to the use of MMCAs, DCE-MRI analysis using two-compartment, unidirectional pharmacokinetic models, which neglect the efflux of the contrast agent from the interstitial space back to intravascular space, could potentially lead to altered results when compared to the two-compartment, bidirectional extended Tofts model used within this study . Instead of longitudinal measurements, mice were killed after each time point to enable ex vivo analysis of tumor tissue. Since it is very difficult to compare the exact same ROIs in histology and MRI, ex vivo analysis is able to validate the key findings of in vivo DCE-MRI, but no direct correlations between histology and MRI can be drawn. While this study performed qualitative ex vivo analysis of tumor tissue, further quantitative ex vivo analysis would enable stronger validation of in vivo MRI data. In this regard, immunohistochemical analysis of tumor vasculature could be extended to αSMA for more detailed analysis of vessel maturity. The initial increase in size of 4T1 tumors after ICI therapy may suggest pseudoprogression. However, to validate a subsequent decrease in tumor size and to assess long-term therapy response, tumors would have had to be observed for a longer period of time including overall survival analyses, which was not possible due to the animal protocol.
Divergent trends in treatment responses were observed between the highly malignant 4T1 and low malignant 67NR tumor model. To further validate that the differing response patterns observed in this study are caused by the divergent degrees of malignancy, the study should be expanded to other tumor entities. It should ideally also include a direct comparison of responding and non-responding lesions.
In summary, DCE-MRI was able to noninvasively assess different vascular response patterns to treatment with the multi-kinase inhibitor sorafenib or immune checkpoint blockade between low malignant and highly malignant tumors. By capturing early changes of tumor vasculature, a key component of the tumor microenvironment, it may allow for personalized treatment decisions, with the possibility to guide early therapeutic decisions in patients that fail fast.
Availability of data and materials
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
Area under the curve
Fast low angle shot
Field of view
Hematoxylin and eosin
Immune checkpoint inhibitors
- K trans :
Volume transfer constant
Laser ablation-inductively coupled plasma-mass spectrometry
Macromolecular contrast agents
Magnetic resonance imaging
Rapid acquisition with relaxation enhancement
Variable repetition time
Region of interest
- slopemax :
Small molecular contrast agents
Transmission electron microscopy
Hanahan D, Weinberg RA. The hallmarks of cancer. Cell. 2000;100:57–70.
Lugano R, Ramachandran M, Dimberg A. Tumor angiogenesis: causes, consequences, challenges and opportunities. Cell Mol Life Sci. 2020;77:1745–70.
Barlow KD, Sanders AM, Soker S, et al. Pericytes on the tumor vasculature: jekyll or hyde? Cancer Microenviron. 2013;6:1–17.
Siemann DW. The unique characteristics of tumor vasculature and preclinical evidence for its selective disruption by tumor-vascular disrupting agents. Cancer Treat Rev. 2011;37:63–74.
Treps L, Ager A, Hida K. Editorial: tumor vessels as directors of the tumor microenvironment: new findings, current challenges & perspectives. Front Cell Dev Biol. 2022;10: 885670.
Bhise NS, Shmueli RB, Sunshine JC, et al. Drug delivery strategies for therapeutic angiogenesis and antiangiogenesis. Expert Opin Drug Deliv. 2011;8:485–504.
Azzi S, Hebda JK, Gavard J. Vascular permeability and drug delivery in cancers. Front Oncol. 2013;3:211.
Teleanu RI, Chircov C, Grumezescu AM, et al. Tumor angiogenesis and anti-angiogenic strategies for cancer treatment. J Clin Med. 2019;9:E84.
Klein D. The tumor vascular endothelium as decision maker in cancer therapy. Front Oncol. 2018;8:367.
Ergün S, Tilki D, Kilic N. Vascular stabilization. In: Schwab M, editor. Encyclopedia of cancer. Berlin: Springer; 2018. p. 1–8.
Yang T, Xiao H, Liu X, et al. Vascular normalization: a new window opened for cancer therapies. Front Oncol. 2021;11: 719836.
Schaaf MB, Garg AD, Agostinis P. Defining the role of the tumor vasculature in antitumor immunity and immunotherapy. Cell Death Dis. 2018;9:115.
Lau D, McLean MA, Priest AN, et al. Multiparametric MRI of early tumor response to immune checkpoint blockade in metastatic melanoma. J Immunother Cancer. 2021;9: e003125.
Zhao Y, Ting KK, Coleman P, et al. The tumour vasculature as a target to modulate leucocyte trafficking. Cancers. 2021;13:1724.
Saida Y, Brender JR, Yamamoto K, et al. Multimodal molecular imaging detects early responses to immune checkpoint blockade. Cancer Res. 2021;81:3693–705.
Ghasemali S, Farajnia S, Barzegar A, et al. New developments in anti-angiogenic therapy of cancer, review and update. Anticancer Agents Med Chem. 2021;21:3–19.
Aalders KC, Tryfonidis K, Senkus E, et al. Anti-angiogenic treatment in breast cancer: facts, successes, failures and future perspectives. Cancer Treat Rev. 2017;53:98–110.
Haibe Y, Kreidieh M, El Hajj H, et al. Resistance mechanisms to anti-angiogenic therapies in cancer. Front Oncol. 2020;10:221.
Pezzella F. Mechanisms of resistance to anti-angiogenic treatments. Cancer Drug Resist. 2019;2:595–607.
Nagy JA, Chang S-H, Shih S-C, et al. Heterogeneity of the tumor vasculature. Semin Thromb Hemost. 2010;36:321–31.
Hida K, Ohga N, Hida Y, et al. Significance of anti-angiogenic therapy in head and neck cancer—heterogeneity of tumor endothelium. Jpn Dent Sci Rev. 2010;46:26–32.
Türkbey B, Thomasson D, Pang Y, et al. The role of dynamic contrast-enhanced MRI in cancer diagnosis and treatment. Diagn Interv Radiol. 2010;16:186–92.
Merz M, Komljenovic D, Zwick S, et al. Sorafenib tosylate and paclitaxel induce anti-angiogenic, anti-tumour and anti-resorptive effects in experimental breast cancer bone metastases. Eur J Cancer. 2011;47:277–86.
Jeon TY, Kim CK, Kim J-H, et al. Assessment of early therapeutic response to sorafenib in renal cell carcinoma xenografts by dynamic contrast-enhanced and diffusion-weighted MR imaging. Br J Radiol. 2015;88:20150163.
Campos M, Candelária I, Papanikolaou N, et al. Perfusion magnetic resonance as a biomarker for sorafenib-treated advanced hepatocellular carcinoma: a pilot study. GE - Port J Gastroenterol. 2019;26:260–7.
Satoh N, Yamada Y, Kinugasa Y, et al. Angiopoietin-1 alters tumor growth by stabilizing blood vessels or by promoting angiogenesis. Cancer Sci. 2008;99:2373–9.
Keating GM. Sorafenib: a review in hepatocellular carcinoma. Target Oncol. 2017;12:243–53.
Petitprez F, Meylan M, de Reyniès A, et al. The tumor microenvironment in the response to immune checkpoint blockade therapies. Front Immunol. 2020;11:784.
Helfen A, Große Hokamp N, Geyer C, et al. Target-specific imaging of cathepsin and S100A8/A9 reflects specific features of malignancy and enables estimation of tumor malignancy. Mol Imaging Biol. 2019;22(1):66–72.
Becker A, Große Hokamp N, Zenker S, et al. Optical in vivo imaging of the alarmin S100A9 in tumor lesions allows for estimation of the individual malignant potential by evaluation of tumor-host cell interaction. J Nucl Med. 2015;56:450–6.
Gerwing M, Hoffmann E, Kronenberg K, et al. Multiparametric MRI enables for differentiation of different degrees of malignancy in two murine models of breast cancer. Front Oncol. 2022;12:1000036.
David S, Park J-K, van Meurs M, et al. Acute administration of recombinant angiopoietin-1 ameliorates multiple-organ dysfunction syndrome and improves survival in murine sepsis. Cytokine. 2011;55:251–9.
Kissel M, Berndt S, Fiebig L, et al. Antitumor effects of regorafenib and sorafenib in preclinical models of hepatocellular carcinoma. Oncotarget. 2017;8:107096–108.
Rupp T, Genest L, Babin D, et al. Anti-CTLA-4 and anti-PD-1 immunotherapies repress tumor progression in preclinical breast and colon model with independent regulatory T cells response. Transl Oncol. 2022;20: 101405.
Kim K, Skora AD, Li Z, et al. Eradication of metastatic mouse cancers resistant to immune checkpoint blockade by suppression of myeloid-derived cells. Proc Natl Acad Sci U S A. 2014;111:11774–9.
Jahng G-H, Li K-L, Ostergaard L, et al. Perfusion magnetic resonance imaging: a comprehensive update on principles and techniques. Korean J Radiol. 2014;15:554–77.
Sourbron SP, Buckley DL. On the scope and interpretation of the tofts models for DCE-MRI. Magn Reson Med. 2011;66:735–45.
Theis D, Keil B, Heverhagen JT, Béhé M, Klose KJ, Fiebich M. Determination of the arterial input function in mouse-models using clinical MRI. Z F Med Phys. 2008;18(2):85–90.
Blockley NP, Jiang L, Gardener AG, et al. Field strength dependence of R1 and R2* relaxivities of human whole blood to prohance, vasovist, and deoxyhemoglobin. Magn Reson Med. 2008;60:1313–20.
Vanchinathan V, Mizramani N, Kantipudi R, et al. the vascular marker CD31 also highlights histiocytes and histiocyte-like cells within cutaneous tumors. Am J Clin Pathol. 2015;143:177–85.
Kimm MA, Haas H, Stölting M, et al. Targeting endothelin receptors in a murine model of myocardial infarction using a small molecular fluorescent probe. Mol Pharm. 2020;17:109–17.
Becker F, Kebschull L, Rieger C, et al. Bryostatin-1 attenuates ischemia-elicited neutrophil transmigration and ameliorates graft injury after kidney transplantation. Cells. 2022;11:948.
Gerwing M, Herrmann K, Helfen A, et al. The beginning of the end for conventional RECIST — novel therapies require novel imaging approaches. Nat Rev Clin Oncol. 2019;16:442–58.
Choi J-I, Imagawa DK, Bhosale P, et al. Magnetic resonance imaging following treatment of advanced hepatocellular carcinoma with sorafenib. Clin Mol Hepatol. 2014;20:218–22.
Ravoori MK, Nishimura M, Singh SP, et al. Tumor T1 relaxation time for assessing response to bevacizumab anti-angiogenic therapy in a mouse ovarian cancer model. PLoS ONE. 2015;10: e0131095.
Folkman J. Tumor angiogenesis: therapeutic implications. N Engl J Med. 1971;285:1182–6.
Goel S, Duda DG, Xu L, et al. Normalization of the vasculature for treatment of cancer and other diseases. Physiol Rev. 2011;91:1071–121.
Yang J, Liao C, Liu Y, et al. MR imaging biomarkers evaluating vascular normalization window after anti-vessel treatment. Oncotarget. 2018;9:11964–76.
Liu Z, Zhao Q, Zheng Z, et al. Vascular normalization in immunotherapy: a promising mechanisms combined with radiotherapy. Biomed Pharmacother. 2021;139:111607.
Zhang J, Winters K, Kiser K, et al. Assessment of tumor treatment response using active contrast encoding (ACE)-MRI: comparison with conventional DCE-MRI. PLoS ONE. 2020;15:e0234520.
Zheng X, Zhang N, Qian L, et al. CTLA4 blockade promotes vessel normalization in breast tumors via the accumulation of eosinophils. Int J Cancer. 2020;146:1730–40.
Johansson A, Hamzah J, Payne CJ, et al. Tumor-targeted TNFα stabilizes tumor vessels and enhances active immunotherapy. Proc Natl Acad Sci USA. 2012;109:7841–6.
Zheng X, Fang Z, Liu X, Deng S, Zhou P, Wang X, Zhang C, Yin R, Hu H, Chen X, Han Y. Increased vessel perfusion predicts the efficacy of immune checkpoint blockade. J Clinic Invest. 2018;128(5):2104–15.
Hubert P, Roncarati P, Demoulin S, et al. Extracellular HMGB1 blockade inhibits tumor growth through profoundly remodeling immune microenvironment and enhances checkpoint inhibitor-based immunotherapy. J Immunother Cancer. 2021;9: e001966.
Marrari A, Wagner AJ, Hornick JL. Predictors of response to targeted therapies for gastrointestinal stromal tumors. Arch Pathol Lab Med. 2012;136:483–9.
Lüönd F, Tiede S, Christofori G. Breast cancer as an example of tumour heterogeneity and tumour cell plasticity during malignant progression. Br J Cancer. 2021;125:164–75.
Cyran CC, Schwarz B, Paprottka PM, et al. In vivo monitoring of sorafenib therapy effects on experimental prostate carcinomas using dynamic contrast-enhanced MRI and macromolecular contrast media. Cancer Imaging. 2013;13:557–66.
Cyran CC, Sennino B, Fu Y, et al. Permeability to macromolecular contrast media quantified by dynamic MRI correlates with tumor tissue assays of vascular endothelial growth factor (VEGF). Eur J Radiol. 2012;81:891–6.
Raatschen H-J, Simon GH, Fu Y, et al. Vascular permeability during antiangiogenesis treatment: MR imaging assay results as biomarker for subsequent tumor growth in rats. Radiology. 2008;247:391–9.
Li K-L, Wilmes LJ, Henry RG, et al. Heterogeneity in the angiogenic response of a BT474 human breast cancer to a novel vascular endothelial growth factor-receptor tyrosine kinase inhibitor: assessment by voxel analysis of dynamic contrast-enhanced MRI. J Magn Reson Imaging. 2005;22:511–9.
We thank Claudia Terwesten-Solé, Karin Gäher, Klaudia Niepagenkemper, Richard Holtmeier and Dr. Stephan Niland for excellent technical support.
Open Access funding enabled and organized by Projekt DEAL. This study was supported by the German Research Foundation (DFG, GE 3336/1-1, grant no. 446302350; SFB1009-A09 and -Z02, grant no. 194468054, CRC1450, grant no. 431460824), the Interdisciplinary Center for Clinical Research Münster (PIX) and the Medical Faculty of the University of Münster (fellowships to EH, MG).
Ethics approval and consent to participate
All animal experiments were carried out in accordance with the ARRIVE guidelines and were approved by the responsible authorities (Landesamt für Natur, Umwelt und Verbraucherschutz NRW, Protocol No. 81-02.04.2018.A010).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Immunohistochemical analysis of the intratumoral immune cell infiltrate. Method 2: Laser ablation-inductively coupled plasma-mass spectrometry. Fig. S1: Analysis of plasma volume fraction vp. Fig. S2: Correlation between T1 relaxation times after contrast agent injection and LA-ICP-MS analysis. Fig. S3: Architecture of sorafenib-treated 4T1 tumors. Fig. S4: Immunohistochemical analysis of the intratumoral immune cell infiltrate after immune checkpoint blockade. Table S1: Descriptive statistics of 67NR tumors. Table S2: Analysis of significance of 67NR tumors. Table S3: Descriptive statistics of 4T1 tumors. Table S4: Analysis of significance of 4T1 tumors.
About this article
Cite this article
Hoffmann, E., Gerwing, M., Krähling, T. et al. Vascular response patterns to targeted therapies in murine breast cancer models with divergent degrees of malignancy. Breast Cancer Res 25, 56 (2023). https://doi.org/10.1186/s13058-023-01658-9