Theranostics 2019; 9(24):7370-7383. doi:10.7150/thno.37750
Normalization of Tumor Vasculature by Oxygen Microbubbles with Ultrasound
1. Department of Biomedical Engineering and Environmental Sciences, National Tsing Hua University, Hsinchu, Taiwan.
2. Institute of Bioinformatics and Structural Biology & Department of Medical Sciences, National Tsing Hua University, Hsinchu, Taiwan.
3. Department of Neurosurgery, Chang Gung Memorial Hospital, Taoyuan, Taiwan.
Ho YJ, Chu SW, Liao EC, Fan CH, Chan HL, Wei KC, Yeh CK. Normalization of Tumor Vasculature by Oxygen Microbubbles with Ultrasound. Theranostics 2019; 9(24):7370-7383. doi:10.7150/thno.37750. Available from http://www.thno.org/v09p7370.htm
Tumor microenvironment influences the efficacy of anti-cancer therapies. The dysfunctional tumor vasculature limits the efficiency of oxygenation and drug delivery to reduce treatment outcome. A concept of tumor vascular normalization (VN), which inhibits angiogenesis to improve vessel maturity, blood perfusion, and oxygenation, has been demonstrated under the anti-angiogenic therapy. The efficiency of drug delivery and penetration is increased by enhancing perfusion and reducing interstitial fluid pressure during the time window of VN. However, anti-angiogenic agents only induce transient VN and then prune vessels to aggravate tumor hypoxia. To repair tumor vessels without altering vessel density, we proposed to induce tumor VN by local oxygen release via oxygen microbubbles with ultrasound. With tumor perfusion enhancement under ultrasound contrast imaging tracing, the time window of VN was defined as 2-8 days after a single oxygen microbubble treatment. The enhanced tumor oxygenation after oxygen microbubble treatment inhibited hypoxia inducible factor-1 alpha (HIF-1α)/vascular endothelial growth factor (VEGF) pathway to improve the morphology and function of tumor vasculature. The pericyte coverage and Hoechst penetration of tumor vessels increased without any changes to the vessel density. Finally, the intratumoral accumulation of anti-cancer drug doxorubicin could be increased 3-4 folds during tumor VN. These findings demonstrate that regulating tumor oxygenation by oxygen microbubbles could normalize dysfunctional vessels to enhance vascular maturity, blood perfusion, and drug penetration. Furthermore, ultrasound perfusion imaging provides a simple and non-invasive way to detect the VN time window, which increases the feasibility of VN in clinical cancer applications.
Keywords: oxygen microbubbles, tumor perfusion, vascular normalization, oxygenation, ultrasound
Tumor oxygenation regulates cellular metabolism to influence the malignancy of the tumor microenvironment. Since the proliferation rate of tumor cells is faster than angiogenesis, the inadequate oxygen supply from blood circulation causes the tumor microenvironment to become hypoxic. In order to compensate for adequate oxygen delivery, increased expression of hypoxia inducible factor-1 alpha (HIF-1α) activates tumor cells to secrete vascular endothelial growth factor (VEGF) for angiogenesis . However, the overexpression of VEGF not only promotes neovascularization, but also aggravates vascular abnormity . The immature tumor vessels with leaky arrangement of endothelial cells, absent pericyte coverage, and incomplete basement membrane contribute to structural and functional abnormalities of the vascular network. Although tumor vascular density is promoted by angiogenesis, the irregular distribution of vasculature and heterogeneous blood perfusion restricts the efficiency of oxygen transport and impedes tumor reoxygenation . Under hypoxic conditions, overexpression of HIF-1α changes the proteome and genome of tumor cells, causing the tumor microenvironment to contribute to treatment resistance and metastasis . Hence, the HIF-1α/VEGF pathway regulates the deterioration of tumor microenvironments including oxygenation, vessel maturity, blood perfusion, and hypoxia to restrict the efficacy of tumor therapy.
The concept of tumor VN induced by US with O2-MBs. The local oxygen release within tumors during O2-MBs treatment enhanced oxygenation and inhibited hypoxia/angiogenesis pathway to accomplish VN.(Click on the image to enlarge.)
In the last two decades, there is a growing interest in regulating the tumor microenvironment to decrease malignancy and improve treatment outcome . The concept is that “vascular normalization (VN)” remodels tumor vasculature into normal and mature phenotype under the anti-angiogenic status, which repairs vessel function to enhance blood perfusion and oxygenation . Bevacizumab (Avastin®), a clinical anti-angiogenic agent, has been used to promote endothelial cell arrangement tightening, pericyte coverage, and blood perfusion for normalizing tumor vasculature at 1-4 days after treatment [7, 8]. Since anti-angiogenic agents only induce transient VN and then prune vessels to aggravate tumor hypoxia, the tumor oxygenation, vessel morphology, or blood perfusion should be traced over time to define the normalized window for combined therapies . Previous studies have combined chemo- or radio-therapy in the VN time window to improve the treatment outcome by the increased efficiency of drug delivery and tumor oxygenation because of the enhanced tumor perfusion [7, 8, 10]. Furthermore, regulating the upstream or downstream genes of VEGF, like prolyl hydroxylase domain-containing protein 2 (PHD2), angiopoietin-1, and epidermal growth factor receptor, also provides various pathways for accomplishing VN .
Oxygen plays an essential role in cellular function and metabolism. Microbubbles comprising a lipid shell and gaseous core can enhance the contrast of ultrasound (US) imaging and deliver various therapeutic gases or drugs for tumor theranosis applications [12, 13]. Previous studies have used oxygen microbubbles (O2-MBs) to locally release oxygen into tumors by external US stimulation, and then improve therapeutic efficacy by reducing hypoxic treatment resistance [14-18]. Since the elevated oxygenation decreases HIF-1α expression, the downstream factor VEGF would be inhibited to restrict tumor angiogenesis and deterioration . Therefore, we speculated that the promotion of tumor oxygenation by US-mediated O2-MBs destruction might regulate the expression of HIF-1α/VEGF pathway to accomplish VN (Figure 1). The feasibility of VN for murine prostate cancer induced by local oxygen release was evaluated to define the VN inducible dose of O2-MBs and corresponding normalized time window. In contrast to anti-angiogenic agents, VN induced by O2-MBs treatment might repair the morphology and function of tumor vessels without decreasing vascular density. In addition, tumor perfusion traced by ultrasound contrast imaging provides a potential and convenient way to evaluate the time window of VN for assisting clinical cancer therapy.
Materials and Methods
Characteristics of O2-MBs
To modify the optimal fabrication of O2-MBs, the gas mixture with various volume ratios of perfluoropropane (C3F8) and oxygen was infused into the phospholipid solution. The phospholipid solution in a sealed vial was composed of 2.5 mg of 1,2-Distearoyl-sn-glycero-3-phosphorylcholine (DSPC) and 1 mg of 1,2-distearoyl-sn-glycero-3-phosphoethanolamine-N-[10-(trimethoxysilyl)undecanamide(polyethylene glycol)-2000] (DSPE-PEG2000; Avanti Lipids Polar, Alabaster, AL, USA) in 0.8 mL of glycerol-phosphate-buffered saline (PBS, 5 wt%). The phospholipid solution was degassed and infused with the gas mixture with various volume ratios of perfluoropropane (C3F8; ABCR GmbH, Karlsruhe, Germany) and oxygen. Then, O2-MBs were fabricated via 45 seconds of violent shaking by an agitator (VIALMIX, Bristol-Myers Squibb Medical Imaging, New York, NY, USA). The additional C3F8 was used to improve the stability of O2-MBs due to the lower water solubility. In contrast, C3F8-MBs were made by the infusion with pure C3F8 gas to be the non-O2 MBs.
The concentration and size distribution of O2-MBs were detected by a coulter counter (Multisizer 3, Beckman Coulter, Miami, FL, USA). In order to estimate the amount of oxygen encapsulated within O2-MBs, the oxygen partial pressure (pO2) within the O2-MBs emulsion (1×107 O2-MBs) was measured by a fiberoptic pO2 probe with an OxyLite 2000 system (Oxford Optronics, Oxford, UK). The oxygen release from O2-MBs was triggered by a sonication (2510, Branson, Danbury, CT, USA) for 5 min at 37 °C to totally destruct O2-MBs. To detect the stability of O2 within O2-MBs, the pO2 levels of O2-MBs in the PBS, degas PBS, and O2-PBS without US sonication were measured . The PBS solution was degassed for 3 min and infused O2 for 1 min to prepare O2-PBS. During O2 infusion, the needle was immersed into PBS to observe the bubble production.
Acoustic characteristics of O2-MBs
The acoustic characteristics were measured to evaluate the difference between C3F8-MBs and O2-MBs. The contrast enhancement of US imaging and changes of pO2 levels were calculated to estimate the in vitro stability of MBs. Diluted MBs emulsion (concentration of 5×106 MBs/mL) was infused into a cylindrical hollow chamber within the 3% agarose phantom at 37 °C. A fiberoptic pO2 probe with an OxyLite 2000 system was inserted into a cylindrical hollow chamber to detect the changes of pO2 levels. The US images were recorded over time with a US imaging system (7-MHz, Model t3000, Terason, Burlington, MA, USA) and the contrast-to-noise ratio produced from the MBs was analyzed. Since the violent oscillation and destruction of MBs might cause vessel damage, the in vitro destruction threshold and inertial cavitation dose (ICD) of MBs were evaluated to establish the safety US parameters for the subsequent local oxygen release. A 2-MHz high-intensity focused US (HIFU; model SU-101, Sonic Concepts, Bothell, WA, USA; diameter of 35 mm, focus length of 55 mm, focal size of 1.2 mm width × 13.3 mm depth) sonication system was used to stimulate O2-MBs (concentration of 5×106 MBs/mL) in a PE-50 tube (Sunpoint scientific Instrument Co, Ltd, Taipei, Taiwan) with the flow rate of 6 mL/h. The 3-cycle US pulses with 2.25 Hz pulse repetition frequency (PRF) were used to stimulate MBs for measuring the destruction threshold and ICD under various peak negative pressures. The destruction threshold of MBs was evaluated by relative contrast enhancement of US images after US sonication. Moreover, a 15-MHz US transducer was used to receive the broadband noises, which were analyzed to calculate the ICD.
All animal experiments were approved by the animal experiment committee at National Tsing Hua University (approval number: 10633). The animal procedures were performed following the guidelines of the Institutional Animal Care Committee. C57BL/6JNarl mice (male, 6-10 weeks old) were intraperitoneally anesthetized with a 1:1 mixture of Rompun 2% (Bayer HealthCare, Leverkusen, Germany) and Zoletil 50 (Virbac, Carros, France).
In vivo vascular bioeffects
Since the ICD produced by O2-MBs violent oscillation and destruction might cause vessel damage, a mouse dorsal skinfold window chamber model (SM100, APJ Trading, Ventura, CA, USA) was used to observe the in vivo vascular bioeffects induced by US-stimulated O2-MBs destruction [20, 21]. An acousto-optical system comprising a 2-MHz HIFU sonication system and an inverted microscope (IX71, Olympus, Tokyo, Japan) was used to recorded intravital images. Window chamber-mice were intravenously injected with 2×107 O2-MBs and then stimulated by a 3-cycle single pulsed US with various peak negative pressures to observe the vessel response. The intravital images before and after US stimulation were collected to evaluate the safety parameters of the HIFU sonication system for O2-MBs destruction, which could locally release oxygen without vascular damage.
Cell line and tumor model
The transgenic adenocarcinoma mouse prostate (TRAMP) cell line was cultured in Dulbecco's modified Eagle's medium supplemented with 1% penicillin-streptomycin and 10% fetal bovine serum (Gibco, Grand Island, NY, USA). Approximately 1×106 TRAMP cells in 100 μL PBS were subcutaneously injected into the right leg of each mouse for tumor growth. Tumors were traced until they became established and reached a mean tumor volume of approximately 100-200 mm3.
Local oxygen release within tumors
The integrated US image-guided treatment system was used to stimulate O2-MBs for local oxygen release within tumors . When the tumor volume reached 100-200 mm3 (day 0), the treatment dose of 2×107 O2-MBs (50 μL) was retro-orbitally injected into mice. A 2-MHz HIFU transducer (2 MPa, 1000-cycle, PRF of 2 Hz) was moved manually by a triaxial platform at 0.5 mm intervals to finish the whole tumor scanning without heating . The acoustic pulse duration was prolonged to decrease the intensity attenuation in solid tumors. The total treatment time for whole tumor sonication was about 30 min. A fiberoptic pO2 probe was inserted into tumors to detect oxygen release kinetics within tumors during O2-MBs treatment by an OxyLite 2000 system. In order to confirm the possible biological damage after O2-MBs treatment, tumors (without pO2 probe insertion) were removed after 24 h and stained by hematoxylin and eosin (H&E) to evaluate tissue hemorrhage and necrosis.
Tumor perfusion tracing
Since tumor perfusion is a potential biomarker for the occurrence of VN [23, 24], US contrast imaging with C3F8-MBs infusion was recorded to trace the relative perfusion intensity of whole tumors over time after O2-MBs treatment. To maintain the concentration of C3F8-MBs within the vasculature during whole tumor scanning, C3F8-MBs (2×109 /mL) were infused by an injection pump with a velocity of 0.3 mL/h. The whole tumor perfusion images were collected with an interval of 0.5 mm. The enhanced signal between pre- and post-C3F8-MBs injection images was quantified to evaluate the tumor perfusion intensity and then normalized to day 0. The VN time window was defined when the tumor perfusion intensity was both higher than the control group and the baseline at day 0. Perfusion density was calculated as the area with enhanced signal in tumors to evaluate the density of the functional vessels. Moreover, the length, width, and depth of tumors were measured on two orthogonal US images to calculate and trace tumor volume.
In vivo microdialysis analysis
Intratumoral oxygen release by US-mediated O2-MBs destruction could improve tumor oxygenation and reduce hypoxia. To verify the correlated gene pathway between oxygenation, angiogenesis, and cell proliferation, the variability of VEGF and transforming growth factor-beta (TGF-β) expression within tumors were traced. The in vivo microdialysis analysis was used to collect intratumoral tissue fluid at day 0, 2, 4, 6, and 8 after O2-MBs treatment. Microdialysis provided a non-sacrificed long-time protein collection in the same mice, which could trace the changes of microenvironment in tumors. The microdialysis probe (MAB 5, Microbiotech, Stockholm, Sweden) was implanted into the center of the tumor and collection was delayed for 30 min to reduce the interference from potential tumor damage. During dialysate collection, a probe was perfused with saline by a microinjection pump (CMA 402, CMA, Stockholm, Sweden) for 30 min (flow rate of 2.0 mL/h). Then, the protein concentrations of VEGF and TGF-β were analyzed in the microdialysis samples using the commercial ELISA kits (Rockland, Pennsylvania, USA) according to the manufacturer's protocols. After reaction, samples were measured using a plate reader system (Tecan Infinite M200, Tecan Trading AG, Switzerland) at an optical density of 450 nm to quantify the concentration of VEGF and TGF-β. The concentration of VEGF and TGF-β at each time point was normalized to that measured at day 0.
Intratumoral protein expression during VN time window
Western blot analysis was used to detect the protein expression in whole tumors during VN time window, the time point was based on the significant perfusion enhancement and VEGF reduction. Tumors (N = 15) were removed at day 4 after O2-MBs treatment to confirm the protein expression during the period of perfusion enhancement. The collection obtained more proteins by homogeneous tissue than tissue fluid to suitably analyze various proteins expression. Tumor tissue was mechanically disrupted with RIPA lysis buffer (APOLLO, CA, USA). The samples were centrifuged at 13,000 rpm for 30 min at 4 °C. Proteins in the supernatants were quantified by Coomassie Protein Assay Reagent (BioRad Laboratories, Hercules, CA, USA). We diluted 30 μg of proteins in Laemmli sample buffer (pH 6.8, 50 mM Tris, 10% (v/v) glycerol, 2% SDS (w/v), 0.01% (w/v) bromophenol blue) and separated proteins by 1D-SDS-PAGE. The proteins were separated onto 0.45 μm Immobilon P membranes (Millipore, Bedford, MA, USA) and the membranes were immersed in a blocking solution with 5% w/v skim milk in TBS-T (50 mM Tris pH 8.0, 0.1% Tween-20 (v/v), and 150 mM NaCl) for 1 h. After that, membranes were incubated with primary antibodies (1: 2000; GeneTex, Irvine, CA, USA) of HIF-1α, VEGF, TGF-β, and PHD2 overnight at 4 °C and washed 4 times with TBST (150 mM NaCl and 0.1% Tween 20 in 10 mM Tris-HCL, pH 7.4) for 10 min. The horseradish peroxidase-conjugated AffiniPure Goat Anti-Rabbit IgG (H+L) (1:10000; Jackson ImmunoResearch, West Grove, PA, USA) was added at room temperature for 1 h, and then membranes were washed 6 times with TBST for 10 min. Finally, membranes were visualized by an enhanced chemiluminescence method (Visual Protein Biotechnology, Taipei, Taiwan) and the expression levels of the proteins were measured by an Image-Quant TL software (GE Healthcare Bio-Sciences, Uppsala, Sweden).
The modulation of tumor microenvironment was confirmed by assessing the function and morphology of vasculature by immunohistochemistry. At day 4 after O2-MBs treatment, mice were intravenously injected with Hoechst 33342 (15 mg/kg; Invitrogen, Carlsbad, CA, USA) and the dye was allowed to circulate for 1 min before the mouse was sacrificed and perfused with 0.9% normal saline to remove the Hoechst from the circulation. The Hoechst stain was used to evaluate the diffusion function of vessels. Furthermore, 20-μm-thick tumor sections were stained for endothelial cells and pericytes using rat anti-mouse CD31 (1: 100; GeneTex, Irvine, CA, USA) and mouse monoclonal alpha-smooth muscle actin (α-SMA) directly conjugated with fluorescein isothiocyanate secondary antibody (1: 100; GeneTex, Irvine, CA, USA). Tumors were separated into three equal parts, with each part defined as central (inner half-radii) or peripheral regions (outer half-radii); the three random fields per tumor region were analyzed. Each tumor had a total of 18 fields for histological quantification: (1) vessel density: density of endothelial cells; (2) vessel diameter: size of tumor vessels; (3) vascular maturity: ratio of pericyte coverage on endothelial cells (α-SMA/CD31); (4) Hoechst intensity: the fluorescence intensity of Hoechst per field; (5) Hoechst leakage: quantification of Hoechst intensity normalized to vessel density [25, 26]; (6) Hoechst penetration: the fluorescence intensity curve of Hoechst from the vessel wall to the tumor tissue.
Clinical anti-cancer drug accumulation
In this study, we used chemotherapeutic drug doxorubicin (DOX) to confirm the improvement of drug accumulation within tumors after enhancing tumor perfusion by VN. Tumor-bearing mice were separated into the control, C3F8-MBs+US, and O2-MBs+US groups. After MBs treatment 4 days, mice were intravenously injected DOX (7 mg/kg) without MBs infusion. After 24 h, tumors were removed to evaluate the intratumoral drug distribution by histological assessment. The red-fluorescence intensity of DOX in histological images was analyzed using MATLAB (MathWorks, Natick, MA, USA) software to quantify the accumulation of DOX in solid tumors. The quantification of microscopic images was followed the above procedure of immunohistochemistry analysis.
Results were presented as mean ± standard deviation. In each graph, the error bars indicated the standard deviation. For in vivo experiments, four to eight mice were recruited to each group. The multiple comparisons were made by one-way ANOVA (SPSS 13.0, SPSS Inc., IBM, Armonk, NY, USA), and the statistical significance was defined as p-value less than 0.05.
Characteristics of O2-MBs
Different volume ratios of C3F8 and O2 were encapsulated in the lipid-shelled O2-MBs. Based on the smaller size and higher concentration than other groups, the optimal fabrication of O2-MBs was made by a C3F8:O2 volume ratio of 7:5 (Figure 2A and Figure S1A). The mean size and concentration of O2-MBs (C3F8:O2=7:5) were 1.02±0.03 μm (mean±standard deviation) and (20.54±3.46) ×109 MBs/mL, respectively. The pO2 levels within the O2-MBs emulsion increased 21±7 mmHg relative to the PBS, and then O2-MBs released oxygen in a sonication bath to promote pO2 levels from 166±7 to 186±8 mmHg (p<0.05; Figure 2B). In order to demonstrate that the enhanced pO2 level was from O2-MBs not buffer, O2-MBs was centrifuged 2000 rcf for 1 min and washed by PBS, called washed O2-MBs. The pO2 levels showed no significant difference between O2-MBs and washed O2-MBs which indicated that O2 was mostly encapsulated within O2-MBs. Besides, the pO2 levels of 1×107 O2-MBs in the PBS and O2-PBS showed more stable than that in the degas PBS, which demonstrated the loading stability of O2 within O2-MBs (Figure S2).
To evaluate the acoustic features of the O2-MBs, the in vitro stability, destruction threshold, and ICD were measured. The contrast enhancement of US images and pO2 levels by C3F8-MBs and O2-MBs revealed no significant reduction after 60 min at 37 °C (Figure S1B). A 2-MHz HIFU with an acoustic pressure of 2 MPa showed 100% MBs destruction, which was used to stimulate O2-MBs destruction and trigger local oxygen release (Figure S1C). In order to evaluate the possible bioeffects induced by the mechanical force from the MBs destruction, the ICD was measured, providing a direct proportion to the acoustic pressures in the C3F8-MBs and O2-MBs (p>0.05; Figure S1D). Since ICD might have an anti-vascular effect that reduces tumor perfusion for anti-vascular therapy , we used a mouse dorsal skinfold window chamber model to confirm the vascular bioeffects during O2-MBs destruction. The intravital images revealed vasoconstriction after US stimulation with 2 MPa and showed vascular disruption under 3 to 5 MPa (Figure 2C). According to the above-mentioned results, US-mediated O2-MBs destruction could increase the pO2 without inducing vascular bioeffects under the acoustic pressure of 2 MPa.
Tumor perfusion enhancement during VN
To study the impact of local oxygen release from O2-MBs on tumors, we used TRAMP tumors to evaluate the feasibility of VN. The normalized intratumoral pO2 was significantly increased and maintained 132±13% over 30 min after O2-MBs treatment (Figure 3A). The intratumoral pO2 levels were increased from 16±4 to 22±6 mmHg at 0 to 60 min (Figure S3). In Figure 3B, the perfusion images from the maximum section of the tumor show higher perfusion intensity in the O2-MBs+US group than in the control and C3F8-MBs+US groups, especially at day 4 and 6 after treatment. Quantification of tumor perfusion showed an enhanced signal of C3F8-MBs on US images. The perfusion intensity ratio in the O2-MBs+US group was 1.00±00, 1.35±0.20, 1.10±0.27, and 1.16±0.24 at day 0, 2, 4, and 6, respectively (Figure 3C). In comparison to the control group, the perfusion density showed no significant reduction by the O2-MBs treatment (Figure 3D). Although angiogenesis is promoted during tumor growth, the vessel density in the control group showed no significant difference due to the tumor volume increasing [7, 27-29]. Since the O2-MBs treatment did not influence the vessel density of the tumors, the subsequent hypoxia aggravation would be mitigated to improve the treatment outcome. These results further confirm that tumor perfusion is enhanced through the repair of vessel function during VN instead of angiogenesis. The trace of tumor volume showed no significant difference in growth rate between the groups, demonstrating that increased tumor oxygenation does not accelerate tumor cell proliferation (Figure 3E).
(A) Volume distribution of O2-MBs with various volume ratios of C3F8 and O2. (B) The pO2 levels within O2-MBs emulsion (C3F8:O2=7:5). Oxygen release was triggered by sonication for 5 min at 37 °C. The legends of w/o and w/ mean without US and with US sonication, respectively. (C) The intravital images were collected to evaluate the in vivo vascular bioeffects induced by O2-MBs destruction. A mouse dorsal skinfold window chamber model was placed on an acousto-optical system. Intravital images revealed the vascular constriction (→) and disruption (➝) during O2-MBs destruction under various acoustic pressures. Scale bar, 100 μm. Quantitative results are presented as mean ± standard deviation and were analyzed by one-way ANOVA (* p<0.05; † p<0.01).(Click on the image to enlarge.)
(A) The normalized intratumoral pO2 was significantly increased by O2 release from O2-MBs during US sonication. (B) The blood perfusion of TRAMP tumors was traced over time by US contrast imaging. The infusion of C3F8-MBs (2×109 /mL) displayed the enhanced signal on US imaging to estimate the perfusion ability. (C) The perfusion intensity ratios were relative to day 0 of each group, which were traced over time after various treatment protocols. The treatment dose of C3F8- or O2-MBs was 2×107 per mouse. Data are presented as box-and-whiskers plots, each dot represents a sample (N=5 to 8 mice per group). (D) Perfusion density ratios showed no difference in the O2-MBs+US group. (E) Tumor volume was traced to compare the growth rate between each group (p>0.05). Data are presented as mean ± standard deviation and were analyzed by one-way ANOVA (* p<0.05; † p<0.01).(Click on the image to enlarge.)
The dose dependence of tumor perfusion enhancement by O2-MBs treatment was also investigated in the present study. In the control group, the end point of perfusion tracing was at day 6 because the subsequent tumor perfusion was always lower than day 0 and decreased over time. In a comparison of the perfusion intensity ratio of the control tumors (w/o O2-MBs injection), the dose of 0.5×107 O2-MBs showed no significant difference from day 2 to day 6 (each p>0.05; Figure 4A). The dose of 1.0×107 O2-MBs displayed perfusion enhancement at day 6 relative to the control group (p<0.05), which was lower than that at day 0. Notably, tumor perfusion relative to day 0 in the 2.0 and 4.0×107 O2-MBs groups maintained the enhancement until day 8 after treatment, but the difference between doses is not significant. At day 4 after treatment, the doses of 0, 0.5, 1.0, 2.0, and 4.0×107 O2-MBs per mouse showed perfusion intensity ratios of 0.70±0.12, 0.71±0.10, 0.97±0.27, 1.10±0.27, and 1.43±0.14, respectively. The correlation coefficient for the relationship between the O2-MBs dose and perfusion intensity ratio was 0.659, 0.795, and 0.719 at day 2, 4, and 6 after treatment, respectively (Figure 4B). Therefore, the VN inducible dose of O2-MBs was 2.0×107 O2-MBs per mouse, with a normalized time window of 2-8 days.
Further, blood perfusion can be enhanced by MBs cavitation, where the increased shear stress on endothelial cells promotes nitric oxide generation and induces vasodilation . The transient effect of perfusion enhancement by MBs cavitation is shown in Figure 4C. However, this vascular enhancing effect by MBs cavitation only presents after 1-30 min and then returns to the baseline at 60 min. Thus, the long-term enhancement of tumor perfusion after O2-MBs treatment was induced by VN instead of MBs cavitation. The biological effects induced by MBs cavitation were also evaluated by tumor perfusion and histological images. Pervious study used non-recovery perfusion reduction to demonstrate the efficiency of anti-vascular effect induced by MBs cavitation [20, 31]. In the present study, intravital images of vessel network on the mouse dorsal skin showed vasoconstriction at small vessels after US stimulation but there was no vessel disruption and hemorrhage (Figure 2C). The enhanced tumor perfusion at 1-30 min after O2-MBs treatment indicated the absence of anti-vascular effect (Figure 4C). Moreover, the histological images stained by H&E revealed intact tumor structure without hemorrhage and necrosis (Figure S4). So that, tumor vessel density showed no significant reduction by the O2-MBs treatment in comparison to the control group (Figure 3D). These results demonstrated that the slight inertial cavitation was produced by O2-MBs under 2 MPa US sonication and biological damage did not be observed.
(A) The threshold of tumor perfusion enhancement at various O2-MBs doses was evaluated to determine the VN inducible dose of O2-MBs. The legend of w/o means without O2-MBs injection. (B) Linear correlation between perfusion intensity ratio and O2-MBs dose. The correlation coefficients are 0.659, 0.795, and 0.719 at day 2, 4, and 6 after O2-MBs treatment, respectively. Data are presented as mean ± standard deviation. (C) The transient effect of perfusion enhancement by MBs cavitation. Data are presented as box-and-whiskers plots and analyzed by one-way ANOVA (* p<0.05; † p<0.01). Each dot represents a sample (N=5 to 8 mice per group).(Click on the image to enlarge.)
Regulation of protein expression after O2-MBs treatment
To estimate the potential gene pathway after tumor VN induced by O2-MBs treatment, the intratumoral protein expression was measured. In the control tumors, the concentration of VEGF showed a non-significant increase of 0.16±0.03 pg/mL at day 0 to 0.22±0.08 pg/mL at day 6 (p>0.05; Figure S5A). The level of VEGF expression normalized to day 0 maintained consistent ratios in the control and C3F8-MBs+US groups; but started to significantly decrease to 0.68±0.24 at day 4 in the O2-MBs+US group (Figure 5A). The TGF-β expression was measured to evaluate the variability of cell proliferation under well oxygenation condition. The concentration and relative TGF-β expression presented no significant difference between each group, revealing consistent results of tumor volume tracing (p>0.05; Figure S5B and Figure 5B). Moreover, the protein expression in whole tumors by Western blot analysis was evaluated at day 4 after treatment because of the significant perfusion enhancement and VEGF reduction. There was a significant reduction of 0.43±0.21, 0.43±0.24, and 0.41±0.26 for HIF-1α, VEGF, and PHD2 expression, respectively, relative to control tumors (Figure 5C). The variability of intratumoral protein expression after O2-MBs treatment demonstrated that promoting tumor oxygenation can reduce tumor hypoxia to inhibit the expression of VEGF for VN.
Variability in protein expression after O2-MBs treatment. The relative expressions of (A) VEGF and (B) TGF-β were traced over time by in vivo microdialysis and measured by ELISA assay. (C) Whole tumoral HIF1-α, VEGF, TGF-β, and PHD2 expression was determined by Western blot. Data are presented as mean ± standard deviation and were analyzed by one-way ANOVA (* p<0.05; † p<0.01).(Click on the image to enlarge.)
Repairing morphology and function of tumor vessels
To demonstrate the occurrence of tumor VN, the morphology and function of tumor vessels were evaluated at day 4 after local oxygen release. Histological images revealed more pericyte (αSMA) coverage on the endothelial cells (CD31) in the central and peripheral regions after O2-MBs treatment (Figure 6A). The extravascular leakage of Hoechst dye was presented to identify the repaired vessel function of perfusion and permeability in the O2-MBs+US group. Quantification showed no significant difference in the vessel density and diameter between each group (each p>0.05; Figure 6B and 6C). The vessel maturity of O2-MBs+US group increased 1.99±0.35- and 1.53±0.39-fold in the tumor central and peripheral regions, respectively, relative to the control tumors (each p<0.01; Figure 6D). Hoechst intensity was analyzed to evaluate the vascular perfusion of tumors, and the O2-MBs+US group revealed significant enhancement with the control and C3F8-MBs+US groups in the tumor center and periphery (each p<0.01; Figure 6E). The Hoechst leakage of the O2-MBs+US group increased 2.20±3.33- and 2.54±4.02-fold in the tumor central and peripheral regions relative to the control tumors (each p<0.01; Figure 6F). The increase in functional tumor vessels was caused by promoting the vascular maturity. Moreover, the distance of Hoechst penetration was significantly increased after O2-MBs treatment. The Hoechst intensity at a penetration distance of 80 μm was 11±8, 12±11, and 41±81% in the control, C3F8-MBs+US, and O2-MBs+US groups, respectively at the tumor center (Figure 6G), and 8±6, 25±20, and 43±19% at the tumor periphery (Figure 6H).
In order to demonstrate the efficiency of intratumoral drug accumulation improved by VN, the chemotherapeutic drug DOX was administrated at day 4 after O2-MBs treatment. Histological images showed more DOX distribution within the tumor core and periphery in the O2-MBs+US group (Figure 7A). The drug accumulation in the control, C3F8-MBs+US, and O2-MBs+US group was 1.96±1.12, 1.87±1.72, and 5.50±5.86, respectively, in the tumor center, and 2.01±1.75, 2.03±2.01, and 9.34±8.82, respectively, in the tumor periphery (Figure 7B). In the histological assessment, the enhancement of vessel maturity, Hoechst intensity, and Hoechst leakage demonstrated that dysfunctional tumor vasculature was restored by normalizing the morphology of the vessels. Moreover, the drug accumulation within tumors can also be improved by enhancing the efficiency of drug delivery during the VN time window.
Histological assessment of tumor vascular morphology and function. (A) CD31 and αSMA staining indicated endothelial cells (red) and pericytes (green). The Hoechst dye (blue) was used to display the ability of vascular perfusion and permeability. Scale bar, 100 μm. Quantification results from histological images included (B) vessel density, (C) vessel diameter, (D) vessel maturity, (E) Hoechst intensity, and (F) Hoechst leakage. The distances of Hoechst penetration into the (G) tumor center and (H) tumor periphery were measured by fluorescence intensity curve from the vessel wall to the tumor tissue. Each group contained 5 mice (N=5). Images were analyzed from nine random fields per region in the tumor center or periphery. Data are presented as box-and-whiskers plots and analyzed by one-way ANOVA (* p<0.05; † p<0.01). Each dot represents an image (n=45).(Click on the image to enlarge.)
The tumor vasculature provides gas exchange and molecular transport to influence the tumor microenvironment. Excessive angiogenesis causes distorted and leaky tumor vessels that reduce the efficiency of blood perfusion. Due to the porous vessel wall, intravascular plasma, molecular, and nano-sized particles passively penetrate into tumor tissue and are retained due to the dysfunctional lymphatic system, called the enhanced permeability and retention (EPR) effect . Although the EPR effect can increase the nanodrug accumulation within tumors, the elevated interstitial fluid pressure (IFP) restricts the penetration distance of nanodrugs, resulting in non-uniform drug distribution . Tumor VN can remodel abnormal tumor vessels into normal and mature phenotypes under anti-angiogenic status . The normalized tumor vasculature has a mature morphology to enhance blood perfusion through repair of vessel function, reduction of IFP by vessel pore shrinking, and prevention of hypoxia by promoting oxygenation, which improves the malignant tumor microenvironment and treatment outcome . Chauhan et al. proved that the intratumoral penetration was enhanced for the 12 nm nanoparticles after VN, but showed no improvement for the 60 and 125 nm nanoparticles . Although tumor VN might decrease the EPR effect by shrinking the vascular pore size, the decreased IFP could enhance the penetration of small-nanosized or molecular drugs to improve therapeutic efficacy [7, 28]. The balance between tumor VN and the EPR effect was mainly associated with the size of drugs in drug penetration. Thus, we used small molecular drug DOX to exclude the interference of drug sizes and demonstrated the enhancement of drug penetration after tumor VN in the present study. Besides, the efficiency of tumor VN induced by O2-MBs treatment was associated with the initial tumor size. The initial tumor diameter of 8 mm reduced perfusion intensity over time after O2-MBs treatment, but tumor necrosis occurred according to the histological images staining by H&E (data was not shown). Thus, the initial tumor diameter to 6 mm (100-200 mm3) was used to investigate tumor VN in this study.
Evaluation of intratumoral drug accumulation. DOX (7 mg/kg) was injected into mice at day 4 after O2-MBs treatment. (A) Histological images for DOX (red) and CD31 staining (endothelial cells, green) show the distribution of drugs within the whole tumor, tumor center (white dotted line), and tumor periphery (yellow dashed line). (B) Quantification of intratumoral DOX accumulation. Each group contained 5 mice (N=5). Images were analyzed from nine random fields per region in the tumor center or periphery. Data are presented as box-and-whiskers plots and analyzed by one-way ANOVA (* p<0.05; † p<0.01). Each dot represents an image (n=45).(Click on the image to enlarge.)
Tumor oxygenation influences HIF-1α expression to regulate the balance between pro- and anti-angiogenesis in neovasculature progression . In tumor hypoxic conditions, tumor cells secrete HIF-1α to activate VEGF expression for angiogenesis. Thus, improving tumor oxygenation could inhibit HIF-1α secretion, and then reduce the downstream VEGF protein level. In the present study, we proposed a potential pathway to induce VN by improving tumor oxygenation. The tumor perfusion enhancement, HIF-1α and VEGF expression reduction, and mature morphology of vessels after O2-MBs treatment demonstrated the feasibility of VN induced by local oxygen release. The increased penetration distance of Hoechst dye indicates the possibility of IFP reduction after VN. Since the improvement of intratumoral drug accumulation after VN shows high correlation with the size of nanodrugs, the small clinical anti-cancer drug DOX was used in our study and was shown to have more drug accumulation within tumors after O2-MBs treatment. Furthermore, US-mediated MBs cavitation has been widely applied to promote cellular and vascular permeability for drug and gene delivery [36, 37]. However, these vascular enhancing effects only present short-term changes to assist with immediate drug delivery. The intratumoral distribution of drugs might be non-uniform and only accumulate around the vessel wall because the dysfunctional vessel network and high IFP have yet to be solved.
Previous studies used commercial US imaging system (mechanical index= 0.3-1.3) with MBs infusion to evaluate the possible vascular bioeffects for contrast enhanced sonography [38-40]. MBs cavitation not only increases the intravascular shear stress to induce vasodilation, the intracellular hydrogen peroxide, calcium ion flux, and ATP generation are also increased to regulate the metabolism of endothelial cells [30, 41, 42]. Besides, the augmentation of blood flow by MBs cavitation has been applied to provide a prevention way for ischemia-reperfusion injury after stroke, thrombolysis, surgery, and organ transplantation [30, 43, 44]. In the present study, this transient enhancement of blood perfusion induced by C3F8-MBs treatment with therapeutic US (mechanical index=1.4) was also observed in the tumor vasculature. The enhanced tumor perfusion only maintained 30 min after C3F8-MBs treatment but revealed no significant difference in the tumor vascular morphology. The intratumoral drug penetration in the present study was investigated at day 4 after C3F8-MBs treatment. The penetration distance of hoechst within tumor periphery showed significant increment in the C3F8-MBs +US group relative to the control tumors. Notably, the time point of hoechst administration not accompanied with MBs injection. Thus, the enhanced penetration of hoechst was not caused by the increased vascular permeability during MBs cavitation. Although the mechanism of penetration improvement after C3F8-MBs treatment was not clear, the change in endothelial cell metabolism induced by the increased shear stress might be the potential reason.
An oxygen-sensor PHD2 regulates endothelial cell arrangement and pericyte coverage in angiogenesis to influence the maturity of vessels . Inhibition of PHD2 expression can reset oxygen sensing in tumors to adapt to hypoxia and induce VN [46-48]. Besides, TGF-β expression is associated with tumor cell proliferation, which was used to evaluate whether the improvement of tumor oxygenation would affect tumor growth or not in our study. The expression of TGF-β also inhibited the degradation of HIF-1α by decreasing the protein levels of PHD2 . Liu et al. blocked TGF-β to decrease HIF-1α and VEGF for inducing tumor VN in breast carcinoma . However, since TGF-β is not the downstream factors of HIF-1α, improving tumor oxygenation showed no significant difference in the protein levels of TGF-β. According to the upstream gene regulation, O2-MBs treatment reduces the level of HIF-1α to affect downstream VEGF and PHD2 expression. The mechanism of VN induced by O2-MBs treatment involves inhibition of HIF-1α, VEGF, and PHD2 expression to regulate multiple hypoxic and angiogenic genes. A previous study reported that dual inhibition of VEGF receptor and angiopoietin-2 could prolong the VN time window compared with VEGF receptor inhibition alone . Since the anti-angiogenic agent is used to inhibit vessel growth, the balance between pro- and anti-angiogenesis for inducing VN is transient [7, 8]. The excessive pruning tumor vessels after VN time window reduced perfusion and increased hypoxia to cause probable treatment resistance and poor prognosis. On the other hand, O2-MBs-induced VN was regulated by hypoxic and angiogenic gene expression under tumor reoxygenation. The non-significant difference in vascular density after O2-MBs treatment might palliate the aggravation of tumor microenvironment, which might be another reason to prolong the VN time window.
Tumor VN can also reverse the immunosuppressive microenvironment because the enhanced tumor perfusion promotes the infiltration and activity of immune cells . The tumor microenvironment usually accumulates abundant myeloid-derived suppressor cells, tumor-associated macrophages, and regulatory T cells to evade host immunosurveillance . The hypoxic condition in tumors reduces the activity of immune cells and increases the recruitment of immune suppressors to contribute to immune suppression by activating HIF-1 and VEGF pathways . The low dose of anti-angiogenic treatment promotes the polarization of tumor-inhibitory M1-like tumor-associated macrophages and the infiltration of CD4+ and CD8+ T cells . Since the high IFP would compress the intratumoral lymphatic vessels, the reduction of IFP by VN might help antigen-presenting cells return to home to lymph nodes for immune activation. Moreover, VN induced by PHD2 inhibition prevented tumor cell migration through the porous vessel wall into the circulation [56, 57]. The shrinking pore size of tumor vessels reduced the probability of tumor metastasis. Therefore, VN not only remodels defective tumor vasculature, but also repairs malignant tumor microenvironment to reverse immunosuppression and impede metastasis. In the future work, the possibility of prolonging the VN time window should be investigated by multiple injections of O2-MBs or by combination of anti-angiogenic agents and O2-MBs. The regulations of the tumor microenvironment for immunosuppression and metastasis should be discussed under the O2-MBs treatment.
In summary, we proposed that local oxygen release by US-mediated O2-MBs destruction regulates the tumor microenvironment to inhibit hypoxia and angiogenesis for accomplishing VN. Tumor vascular maturity and function were improved to enhance tumor perfusion during the VN time window of 2 to 8 days after O2-MBs treatment. And then, enhanced tumor perfusion promotes the efficiency of drug delivery for increasing intratumoral drug accumulation. In light of this view, US perfusion imaging could be applied to trace tumor perfusion in vivo, allowing us to evaluate whether VN is induced or not after O2-MBs or anti-angiogenic treatment. When tumor perfusion was higher than the first-time measurement, the time window of VN for subsequent chemo- or radio-therapy could be defined to establish personal medicine information. This convenient and easy approach of US perfusion imaging can be used as part of clinical diagnosis for precision tumor therapy.
Supplementary figures and tables.
The authors gratefully acknowledge the support of the Ministry of Science and Technology, Taiwan under Grant No. MOST 107-2221-E-007-002, 107-2627-M-007-005, 107-2119-M-182-001, 108-2221-E-007-040-MY3, 108-2221-E-007-041-MY3, and 106-2218-E-007-022-MY3; National Tsing Hua University under Grant No. 107Q2511E1 and 108Q2717E1; Chang Gung Memorial Hospital and National Tsing Hua University Joint Research Program under Grant No. 107-CMRPG3H0221.
The authors have declared that no competing interest exists.
1. Vaupel P, Harrison L. Tumor hypoxia: causative factors, compensatory mechanisms, and cellular response. Oncologist. 2004;9(Suppl 5):4-9
2. Azzi S, Hebda JK, Gavard J. Vascular permeability and drug delivery in cancers. Front Oncol. 2013;3:211
3. Shannon AM, Bouchier-Hayes DJ, Condron CM, Toomey D. Tumour hypoxia, chemotherapeutic resistance and hypoxia-related therapies. Cancer Treat Rev. 2003;29:297-307
4. Vaupel P. The role of hypoxia-induced factors in tumor progression. Oncologist. 2004;9(Suppl 5):10-7
5. Danhier F. To exploit the tumor microenvironment: Since the EPR effect fails in the clinic, what is the future of nanomedicine?. J Control Release. 2016;244:108-21
6. Jain RK. Normalizing tumor vasculature with anti-angiogenic therapy: a new paradigm for combination therapy. Nat Med. 2001;7:987-9
7. Dickson PV, Hamner JB, Sims TL, Fraga CH, Ng CY, Rajasekeran S. et al. Bevacizumab-induced transient remodeling of the vasculature in neuroblastoma xenografts results in improved delivery and efficacy of systemically administered chemotherapy. Clin Cancer Res. 2007;13:3942-50
8. McGee MC, Hamner JB, Williams RF, Rosati SF, Sims TL, Ng CY. et al. Improved intratumoral oxygenation through vascular normalization increases glioma sensitivity to ionizing radiation. Int J Radiat Oncol Biol Phys. 2010;76:1537-45
9. Jain RK. Normalization of tumor vasculature: an emerging concept in antiangiogenic therapy. Science. 2005;307:58-62
10. Tolaney SM, Boucher Y, Duda DG, Martin JD, Seano G, Ancukiewicz M. et al. Role of vascular density and normalization in response to neoadjuvant bevacizumab and chemotherapy in breast cancer patients. Proc Natl Acad Sci U S A. 2015;112:14325-30
11. Carmeliet P, Jain RK. Principles and mechanisms of vessel normalization for cancer and other angiogenic diseases. Nat Rev Drug Discov. 2011;10:417-27
12. Fix SM, Borden MA, Dayton PA. Therapeutic gas delivery via microbubbles and liposomes. J Control Release. 2015;209:139-49
13. Rojas JD, Lin F, Chiang YC, Chytil A, Chong DC, Bautch VL. et al. Ultrasound Molecular Imaging of VEGFR-2 in Clear-Cell Renal Cell Carcinoma Tracks Disease Response to Antiangiogenic and Notch-Inhibition Therapy. Theranostics. 2018;8:141-55
14. Sun J, Yin M, Zhu S, Liu L, Zhu Y, Wang Z. et al. Ultrasound-mediated destruction of oxygen and paclitaxel loaded lipid microbubbles for combination therapy in hypoxic ovarian cancer cells. Ultrason Sonochem. 2016;28:319-26
15. Eisenbrey JR, Shraim R, Liu JB, Li JZ, Stanczak M, Oeffinger B. et al. Sensitization of Hypoxic Tumors to Radiation Therapy Using Ultrasound-Sensitive Oxygen Microbubbles. International Journal of Radiation Oncology Biology Physics. 2018;101:88-96
16. Szabo T. Diagnostic Ultrasound Imaging. Insite Out Burlington: Elsevier Academic. 2004
17. Khan MS, Hwang J, Seo Y, Shin K, Lee K, Park C. et al. Engineering oxygen nanobubbles for the effective reversal of hypoxia. Artif Cells Nanomed Biotechnol. 2018:1-10
18. Kwan JJ, Kaya M, Borden MA, Dayton PA. Theranostic Oxygen Delivery Using Ultrasound and Microbubbles. Theranostics. 2012;2:1174-84
19. D'Angelo G, Duplan E, Boyer N, Vigne P, Frelin C. Hypoxia up-regulates prolyl hydroxylase activity - A feedback mechansim that limits HIF-1 responses during reoxygenation. J Biol Chem. 2003;278:38183-7
20. Ho YJ, Wang TC, Fan CH, Yeh CK. Spatially Uniform Tumor Treatment and Drug Penetration by Regulating Ultrasound with Microbubbles. ACS Appl Mater Interfaces. 2018;10:17784-91
21. ter Haar G. Safety and bio-effects of ultrasound contrast agents. Med Biol Eng Comput. 2009;47:893-900
22. Ho YJ, Yeh CK. Concurrent anti-vascular therapy and chemotherapy in solid tumors using drug-loaded acoustic nanodroplet vaporization. Acta Biomater. 2017;49:472-85
23. Goel S, Duda DG, Xu L, Munn LL, Boucher Y, Fukumura D. et al. Normalization of the vasculature for treatment of cancer and other diseases. Physiol Rev. 2011;91:1071-121
24. Arjaans M, Schroder CP, Oosting SF, Dafni U, Kleibeuker JE, de Vries EG. VEGF pathway targeting agents, vessel normalization and tumor drug uptake: from bench to bedside. Oncotarget. 2016;7:21247-58
25. Viski C, Konig C, Kijewska M, Mogler C, Isacke CM, Augustin HG. Endosialin-Expressing Pericytes Promote Metastatic Dissemination. Cancer Res. 2016;76:5313-25
26. Franco M, Man S, Chen L, Emmenegger U, Shaked Y, Cheung AM. et al. Targeted anti-vascular endothelial growth factor receptor-2 therapy leads to short-term and long-term impairment of vascular function and increase in tumor hypoxia. Cancer Res. 2006;66:3639-48
27. Winkler F, Kozin SV, Tong RT, Chae SS, Booth MF, Garkavtsev I. et al. Kinetics of vascular normalization by VEGFR2 blockade governs brain tumor response to radiation: role of oxygenation, angiopoietin-1, and matrix metalloproteinases. Cancer Cell. 2004;6:553-63
28. Tong RT, Boucher Y, Kozin SV, Winkler F, Hicklin DJ, Jain RK. Vascular normalization by vascular endothelial growth factor receptor 2 blockade induces a pressure gradient across the vasculature and improves drug penetration in tumors. Cancer Res. 2004;64:3731-6
29. El Alaoui-Lasmaili K, Djermoune EH, Tylcz JB, Meng D, Plenat F, Thomas N. et al. A new algorithm for a better characterization and timing of the anti-VEGF vascular effect named "normalization". Angiogenesis. 2017;20:149-62
30. Belcik JT, Mott BH, Xie A, Zhao Y, Kim S, Lindner NJ. et al. Augmentation of limb perfusion and reversal of tissue ischemia produced by ultrasound-mediated microbubble cavitation. Circ Cardiovasc Imaging. 2015:8
31. Liu Z, Gao S, Zhao Y, Li P, Liu J, Li P. et al. Disruption of tumor neovasculature by microbubble enhanced ultrasound: a potential new physical therapy of anti-angiogenesis. Ultrasound Med Biol. 2012;38:253-61
32. Fang J, Nakamura H, Maeda H. The EPR effect: Unique features of tumor blood vessels for drug delivery, factors involved, and limitations and augmentation of the effect. Adv Drug Deliv Rev. 2011;63:136-51
33. Li Y, Wang J, Wientjes MG, Au JL. Delivery of nanomedicines to extracellular and intracellular compartments of a solid tumor. Adv Drug Deliv Rev. 2012;64:29-39
34. Jain RK. Normalizing tumor microenvironment to treat cancer: bench to bedside to biomarkers. J Clin Oncol. 2013;31:2205-18
35. Chauhan VP, Stylianopoulos T, Martin JD, Popovic Z, Chen O, Kamoun WS. et al. Normalization of tumour blood vessels improves the delivery of nanomedicines in a size-dependent manner. Nat Nanotechnol. 2012;7:383-8
36. Hernot S, Klibanov AL. Microbubbles in ultrasound-triggered drug and gene delivery. Adv Drug Deliv Rev. 2008;60:1153-66
37. Lentacker I, De Cock I, Deckers R, De Smedt SC, Moonen CT. Understanding ultrasound induced sonoporation: definitions and underlying mechanisms. Adv Drug Deliv Rev. 2014;72:49-64
38. Miller DL. Overview of experimental studies of biological effects of medical ultrasound caused by gas body activation and inertial cavitation. Prog Biophys Mol Biol. 2007;93:314-30
39. Gao Y, Gao S, Zhao B, Zhao Y, Hua X, Tan K. et al. Vascular effects of microbubble-enhanced, pulsed, focused ultrasound on liver blood perfusion. Ultrasound Med Biol. 2012;38:91-8
40. Bertuglia S. Mechanisms by which low-intensity ultrasound improve tolerance to ischemia-reperfusion injury. Ultrasound Med Biol. 2007;33:663-71
41. Juffermans LJ, van Dijk A, Jongenelen CA, Drukarch B, Reijerkerk A, de Vries HE. et al. Ultrasound and microbubble-induced intra- and intercellular bioeffects in primary endothelial cells. Ultrasound Med Biol. 2009;35:1917-27
42. Belcik JT, Davidson BP, Xie A, Wu MD, Yadava M, Qi Y. et al. Augmentation of Muscle Blood Flow by Ultrasound Cavitation Is Mediated by ATP and Purinergic Signaling. Circulation. 2017;135:1240-52
43. Chen G, Yang L, Zhong L, Kutty S, Wang Y, Cui K. et al. Delivery of Hydrogen Sulfide by Ultrasound Targeted Microbubble Destruction Attenuates Myocardial Ischemia-reperfusion Injury. Sci Rep. 2016;6:30643
44. Cui H, Zhu Q, Gao YH, Xia HM, Tan KB, He Y. et al. Ultrasound Mediated Microbubbles Destruction Augmented Sonolysis: An In Vitro and In Vivo Study. BioMed Res Int. 2017
45. Bautch VL. Endothelial cells form a phalanx to block tumor metastasis. Cell. 2009;136:810-2
46. Mazzone M, Dettori D, Leite de Oliveira R, Loges S, Schmidt T, Jonckx B. et al. Heterozygous deficiency of PHD2 restores tumor oxygenation and inhibits metastasis via endothelial normalization. Cell. 2009;136:839-51
47. De Bock K, De Smet F, De Oliveira RL, Anthonis K, Carmeliet P. Endothelial oxygen sensors regulate tumor vessel abnormalization by instructing phalanx endothelial cells. Journal of Molecular Medicine-Jmm. 2009;87:561-9
48. Koyama S, Matsunaga S, Imanishi M, Maekawa Y, Kitano H, Takeuchi H. et al. Tumour blood vessel normalisation by prolyl hydroxylase inhibitor repaired sensitivity to chemotherapy in a tumour mouse model. Sci Rep. 2017;7:45621
49. McMahon S, Charbonneau M, Grandmont S, Richard DE, Dubois CM. Transforming growth factor beta 1 induces hypoxia-inducible factor-1 stabilization through selective inhibition of PHD2 expression. J Biol Chem. 2006;281:24171-81
50. Liu J, Liao S, Diop-Frimpong B, Chen W, Goel S, Naxerova K. et al. TGF-beta blockade improves the distribution and efficacy of therapeutics in breast carcinoma by normalizing the tumor stroma. Proc Natl Acad Sci U S A. 2012;109:16618-23
51. Peterson TE, Kirkpatrick ND, Huang Y, Farrar CT, Marijt KA, Kloepper J. et al. Dual inhibition of Ang-2 and VEGF receptors normalizes tumor vasculature and prolongs survival in glioblastoma by altering macrophages. Proc Natl Acad Sci U S A. 2016;113:4470-5
52. Huang Y, Kim BYS, Chan CK, Hahn SM, Weissman IL, Jiang W. Improving immune-vascular crosstalk for cancer immunotherapy. Nat Rev Immunol. 2018
53. Vasievich EA, Huang L. The suppressive tumor microenvironment: a challenge in cancer immunotherapy. Mol Pharm. 2011;8:635-41
54. Chouaib S, Messai Y, Couve S, Escudier B, Hasmim M, Noman MZ. Hypoxia promotes tumor growth in linking angiogenesis to immune escape. Front Immunol. 2012;3:21
55. Huang Y, Yuan J, Righi E, Kamoun WS, Ancukiewicz M, Nezivar J. et al. Vascular normalizing doses of antiangiogenic treatment reprogram the immunosuppressive tumor microenvironment and enhance immunotherapy. Proc Natl Acad Sci U S A. 2012;109:17561-6
56. Peinado H, Zhang H, Matei IR, Costa-Silva B, Hoshino A, Rodrigues G. et al. Pre-metastatic niches: organ-specific homes for metastases. Nat Rev Cancer. 2017
57. Kuchnio A, Moens S, Bruning U, Kuchnio K, Cruys B, Thienpont B. et al. The Cancer Cell Oxygen Sensor PHD2 Promotes Metastasis via Activation of Cancer-Associated Fibroblasts. Cell Rep. 2015;12:992-1005
Corresponding author: Department of Biomedical Engineering and Environmental Sciences, National Tsing Hua University, No. 101, Section 2, Kuang-Fu Road, Hsinchu 30013, Taiwan. Tel: +886-3-571-5131 ext. 34240; Fax: +886-3-571-8649; E-mail address: ckyehnthu.edu.tw