Immune modulation resulting from MR-guided high intensity focused ultrasound in a model of murine breast cancer

Fite, Brett Z.; Wang, James; Kare, Aris J.; Ilovitsh, Asaf; Chavez, Michael; Ilovitsh, Tali; Zhang, Nisi; Chen, Weiyu; Robinson, Elise; Zhang, Hua; Kheirolomoom, Azadeh; Silvestrini, Matthew T.; Ingham, Elizabeth S.; Mahakian, Lisa M.; Tam, Sarah M.; Davis, Ryan R.; Tepper, Clifford G.; Borowsky, Alexander D.; Ferrara, Katherine W.

doi:10.1038/s41598-020-80135-1

Download PDF

Article
Open access
Published: 13 January 2021

Immune modulation resulting from MR-guided high intensity focused ultrasound in a model of murine breast cancer

Brett Z. Fite¹,
James Wang¹,
Aris J. Kare^1,2,
Asaf Ilovitsh¹,
Michael Chavez²,
Tali Ilovitsh¹,
Nisi Zhang¹,
Weiyu Chen¹,
Elise Robinson¹,
Hua Zhang³,
Azadeh Kheirolomoom¹,
Matthew T. Silvestrini³,
Elizabeth S. Ingham³,
Lisa M. Mahakian³,
Sarah M. Tam³,
Ryan R. Davis⁴,
Clifford G. Tepper⁵,
Alexander D. Borowsky⁴ &
…
Katherine W. Ferrara¹

Scientific Reports volume 11, Article number: 927 (2021) Cite this article

4599 Accesses
46 Citations
1 Altmetric
Metrics details

Subjects

Abstract

High intensity focused ultrasound (HIFU) rapidly and non-invasively destroys tumor tissue. Here, we sought to assess the immunomodulatory effects of MR-guided HIFU and its combination with the innate immune agonist CpG and checkpoint inhibitor anti-PD-1. Mice with multi-focal breast cancer underwent ablation with a parameter set designed to achieve mechanical disruption with minimal thermal dose or a protocol in which tumor temperature reached 65 °C. Mice received either HIFU alone or were primed with the toll-like receptor 9 agonist CpG and the checkpoint modulator anti-PD-1. Both mechanical HIFU and thermal ablation induced a potent inflammatory response with increased expression of Nlrp3, Jun, Mefv, Il6 and Il1β and alterations in macrophage polarization compared to control. Furthermore, HIFU upregulated multiple innate immune receptors and immune pathways, including Nod1, Nlrp3, Aim2, Ctsb, Tlr1/2/4/7/8/9, Oas2, and RhoA. The inflammatory response was largely sterile and consistent with wound-healing. Priming with CpG attenuated Il6 and Nlrp3 expression, further upregulated expression of Nod2, Oas2, RhoA, Pycard, Tlr1/2 and Il12, and enhanced T-cell number and activation while polarizing macrophages to an anti-tumor phenotype. The tumor-specific antigen, cytokines and cell debris liberated by HIFU enhance response to innate immune agonists.

Low Intensity Focused Ultrasound (LOFU)-mediated Acoustic Immune Priming and Ablative Radiation Therapy for in situ Tumor Vaccines

Article Open access 29 October 2019

Karin A. Skalina, Saurabh Singh, … Chandan Guha

Microwave ablation of primary breast cancer inhibits metastatic progression in model mice via activation of natural killer cells

Article Open access 08 May 2020

Muxin Yu, Hong Pan, … Shui Wang

Molecularly targeted photothermal ablation improves tumor specificity and immune modulation in a rat model of hepatocellular carcinoma

Article Open access 17 December 2020

Nina M. Muñoz, Crystal Dupuis, … Rahul A. Sheth

Introduction

High intensity focused ultrasound (HIFU) is a non-invasive method for the treatment of solid tumors, including those grown in pancreas, liver, kidney, and prostate^1,2,3,4. HIFU is capable of delivering intense acoustic energy to a small focus with high spatial precision, preventing damage to surrounding tissue while resulting in the rapid destruction of tissue within the focal volume through mechanical disruption, thermal damage, or a combination of the two. Furthermore, magnetic resonance guided focused ultrasound (MRgFUS) affords the ability to non-invasively measure temperature in real time, facilitating precise control of temperature and thermal dose^5,6,7,8.

HIFU ablation can exert both thermal and mechanical effects on the targeted tissue. Thermal effects are produced when tissue absorbs acoustic energy, resulting in an increase in temperature; the treatment time and acoustic intensity can be varied to achieve a wide variety of temperature elevations, from only 1–2 °C up to the boiling point of tissue. Mechanical effects primarily result from the production of cavitation bubbles (histotripsy) or microscopic boiling bubbles (boiling histotripsy)^9,10. Mechanical disruption typically results in homogenization and eventual resorption of the liquefied tissue with a narrow transition region^11,12. Conversely, thermal ablation has a broader transition region due to heat diffusion and results in the formation of a scar, which is not as easily resorbed and is resistant to reinvasion. Both mechanical and thermal HIFU rapidly and effectively destroy tissue and have been employed for local control of tumors.

Pre-clinical results suggest local ablative techniques can cause immunogenic cell death (ICD), increased tumor infiltration of macrophages, CD4⁺ and CD8⁺ lymphocytes¹³, release of endogenous danger signals¹⁴ and tumor specific antigens¹⁵, and can result in an antitumor response^16,17,18,19. Clinical results also demonstrate modulation of the immune system following HIFU^13,20. However, ablation alone generates only a weak abscopal effect^16,21; the initial inflammatory response, mediated by HIFU’s release of damage associated molecular patterns (DAMPs) from damaged or destroyed cells is not accompanied by simultaneous pathogen associated molecular pattern (PAMP) recognition, which is likely one major requirement for generating a strong immune response²². Thus, without additional stimuli to innate immune receptors²³, HIFU alone is insufficient to overcome tumor-mediated immune suppression in the vast majority of cases. Nevertheless, the cytokines and cell debris liberated by HIFU can be harnessed for immune stimulation when combined with immunotherapies.

Previously, we have combined HIFU with checkpoint modulation (anti-PD-1) and toll-like receptor (TLR) 9 ligation (using CpG) to generate a durable cure in a majority of mice^15,21 in a “priming” strategy. Our work has demonstrated the importance of timing and sequence of multiple immunomodulatory therapeutics^24,25,26. Aside from debulking, HIFU’s ability to release tumor-specific antigen and its effect on the tumor immune environment can synergize with immunotherapeutic agents. However, the dynamism of the tumor cytokine and immune cell population profile in response to immunotherapeutics and HIFU necessitates further investigation to optimize their combination. Furthermore, the variability in tissue effects generated through different HIFU treatment schemes (e.g. thermal ablation, hyperthermia, histotripsy, boiling histotripsy, low intensity ultrasound) may necessitate different therapeutic combinations to achieve optimal therapeutic effects. Herein, we aim to elucidate the immune effects of HIFU and its combination with immunotherapy. Specifically, we use a combination of anti-PD-1 therapy to increase T-cell response and the TLR9 agonist CPG to activate the myeloid compartment and increase T-cell infiltration within the tumor. Without a robust T-cell infiltrate in the tumor, PD-1 treatment alone may prove less efficacious. We first start with dimensional reduction analysis of transcriptomes from various treatments to understand the differences across treatment combinations. We then demonstrate that while all therapies reduce tumor burden and downregulate cancer and proliferation-related genes, HIFU induces a potent inflammatory response, which is mitigated by priming with an innate immune agonist. Further, thermal ablation upregulates innate immune receptors, and this effect is enhanced by addition of an agonist. Activation of the innate immune system, especially in concert with agonists, generates a systemic, adaptive response, mediated through T cells, B220⁺ cells, and cytokines.

Results

Overall summary of the magnitude of the effect

We hypothesized that different treatment protocols would elicit distinct biological effects observable from transcriptomic expression data. Our choice of time points was based on the ultrasound parameters we chose to employ and our experience with their biological effects and immunohistochemistry studies; namely, we found cell death/destruction was significantly more rapid with mechanical ablation than with thermal ablation. Previously reported results with histotripsy²⁷ have also suggested cell death immediately following treatment. Here, matching immunohistochemistry (IHC) datasets at 72 h and 7 days and flow cytometry at 72 h were acquired for both treatments. These datasets were augmented by an earlier time point for mechanical ablation (24 h for transcriptomics and 48 h for IHC) and a later dataset for thermal ablation (7 day transcriptomics).

To test this hypothesis, we first performed principal component analysis on whole transcriptomes from our NDL tumor-bearing mouse cohorts and found that different treatments (Fig. 1A) have distinct clustering patterns. Both the directly-treated (AI-T) and distant (AI-C) tumors of animals that were treated with immunotherapy and thermal ablation clustered together (Fig. 1B). Tumors treated with thermal ablation alone (A-T) clustered together and were distinct from control tumors and tumors from animals that received immunotherapy prior to thermal ablation. AI-T tumors were 2.2 times further from control tumors than A-T tumors (calculated as the mean Euclidean distance between each cohort and control across all principal components), and AI-C tumors were 5.1 times further from control tumors than A-C tumors. Tumors treated with mechanical HIFU (M-T) and distant tumors in animals that received thermal ablation alone (A-C) clustered with control tumors (NTC), consistent with the lack of response to therapy in these cohorts. A hierarchical clustering of the entire transcriptome (Fig. 1C) similarly demonstrated separate clustering of AI-T, AI-C, and A-T, while M-T, A-C, and NTC tumors clustered together.

To further understand the differences in gene expression across treatment cohorts, we utilized volcano plots of log2 fold change of gene expression (compared to control) with respect to the negative log2 of the adjusted p value of the change in expression (Fig. 1D) to illustrate the differences between the treatment cohorts. In the absence of immunotherapy, gene expression in treated tumors was altered with both mechanical (4673 genes upregulated versus 4011 genes downregulated) and thermal (6406 genes upregulated versus 5553 genes downregulated) ablation. Changes in distant tumors were small; A-C tumors exhibit no change from control (one gene had p_adj < 0.05) (Fig. 1B), and a similar result was observed with mechanical HIFU. After priming with immunotherapy prior to ablation, similar numbers of genes were downregulated and upregulated as compared to HIFU alone (4560 genes upregulated and 4755 genes downregulated for AI-T) in the treated tumor. However, the magnitude of the effect was greater (Fig. 1E) in AI-T tumors; 41% of AI-T genes were differentially expressed with |FC|> 2 (compared to control) as compared to 33% for A-T and 11.8% for M-T. In distant tumors treated with the combined protocol, 5152 genes were upregulated and 4396 downregulated compared to a single gene in A-C tumors. Specific genes associated with adaptive immune response that were upregulated in AI-T compared to control and A-T included Prf1, Gzmb, Batf3, Pdcd1, Ifng, Cd244a, Tigit. A greater fraction of upregulated genes overlap between A-T and AI-T than between A-T and M-T (Fig. 1F) suggesting mechanistic differences between thermal ablation and mechanical HIFU therapy.

HIFU reduces viable tumor

Since HIFU is a non-invasive method for tumor debulking, we sought to characterize the cellular and biomolecular differences of thermal ablation, mechanical HIFU, and thermal ablation combined with immunotherapy. Cancer and proliferation-related genes such as Wnt7b, S100a14, and Erbb2, were downregulated by thermal ablation in directly-treated tumors (A-T) 1 week after ablation, and to a lesser extent, 24 h after mechanical HIFU (M-T) compared to control (NTC) (Fig. 2A). However, this effect did not extend to distant tumors (A-C), consistent with HIFU being a local ablative therapy. Addition of immunotherapy into the ablation protocol further reduced expression of proliferative and cancer-related genes in directly-treated (AI-T) tumors and distant (AI-C) tumors, demonstrating an abscopal effect. We have previously shown that such a combined protocol more than doubles survival in animals with multi-focal disease compared to immunotherapy alone²¹. Similarly, both thermal ablation (Fig. 2B) and mechanical HIFU (Fig. 2C) reduced viable tumor cells only in directly-treated, but not in distant, tumors. Both mechanical HIFU alone and mechanical HIFU pretreated with immunotherapy resulted in increased survival (Fig. 2D) compared to control, but were less efficacious than immunotherapy alone (IT). H&E stained sections of thermally-ablated tumors (Fig. 2E) exhibited a small rim of viable tumor cells while the combination of thermal ablation and immunotherapy eliminated the majority of viable tumor cells (Fig. 2F) and reduced tumor size. Mechanical HIFU caused substantial hemorrhage (Fig. 2G) but remaining viable tumor volume was larger compared to thermally-ablated tumors. However, the addition of immunotherapy substantially reduced viable tumor (Fig. 2H) compared to control tumors (Fig. 2I). Furthermore, the addition of immunotherapy to mechanical ablation reduced viable tumor at 24 h, 72 h, and 7 days post HIFU (Supplementary Fig. S2).

HIFU induces a potent inflammatory response

To understand the potential biological implications of HIFU induced tumor debulking, we performed gene ontology analysis with DAVID. Extracellular matrix (ECM) remodeling and wound healing pathways were induced by thermal ablation and mechanical HIFU (Fig. 3A), but the effect was greatest in thermally-ablated tumors, in which 165 cell adhesion related genes and 100 angiogenesis related genes were upregulated (Supplementary Table S1) compared to 64 and 45 genes, respectively for mechanical ablation.

We further reasoned that wound healing events would be correlated to an inflammatory response and should result in an up-regulation of pro-inflammatory genetic expression. Pro-inflammatory cytokines Il6 (Fig. 3B) and Il1β (Fig. 3C) had increased expression in thermally-ablated tumors, while mechanical HIFU primarily increased expression of the upstream Il1β. Nevertheless, both mechanical HIFU and thermal ablation cohorts were enriched in the Il6 protein production pathway (GO:0032755) with 13 and 22 genes enhanced, respectively. Addition of immunotherapy prior to ablation decreased Il6 gene expression as compared to thermal ablation alone, but remained higher than no treatment control Il6 expression levels. Furthermore, for the combined protocol, the magnitude of the increases in expression of genes associated with the adaptive immune response was greater than that of Il6 expression. In the immunotherapy-ablation cohort, distant tumors (AI-C) had a Cd8a:Il6 expression ratio of 1708 (Fig. 3D) with Il6 expression significantly attenuated. In addition, the immunotherapy-ablation protocol enriched the negative regulatory Il6 pathway (GO:0032715).

Type 3 collagen (Col3) is found in connective tissues, and upregulation of this type of collagen occurs during wound healing processes²⁸. Furthermore, Col3A expression may “play a key role promoting myofibroblast proliferation and progression of fibroepithelial lesions in the breast”²⁹. Thermal ablation increased Col3a protein within the treated tumor (Fig. 3E); however, priming with immunotherapy prior to thermal ablation reduced tumoral Col3a deposition (Fig. 3F). Mechanical ablation (Fig. 3G) did not alter Col3a levels as compared to untreated control tumors, with sparse, scattered Col3a observed throughout the tumor on IHC (Fig. 3H).

Innate immune receptors and signaling pathways are locally upregulated by HIFU and systemically upregulated by adding immunotherapy

Activation of innate immune receptors, especially on myeloid cells, is a key element facilitating a systemic immune response after focal therapy. Pattern recognition receptors (PRRs) are innate immune receptors which recognize endogenous (DAMPs) or exogenous (PAMPs) molecular patterns associated with tissue damage and pathogen infection³⁰. DAMPs, such as heat-shock proteins and HMGB1, are released following tissue damage. In the absence of additional stimuli, DAMP detection triggers a sterile inflammatory response that promotes wound healing and tissue repair. Detection of PAMPs, such as lipopolysaccharides and bacterial DNA, polarize the inflammatory response towards pathogen elimination. Innate immune PAMP sensors must be activated to achieve an anti-tumor immune response.

Ablation alone releases DAMPs and upregulates multiple families of innate receptors, including Toll-like receptors (TLRs), NOD-like receptors (NLRs), cytosolic DNA sensors^30,31,32 (Supplementary Fig. S3), and RIG-I-like receptors (RLRs) (Fig. 4A), which sensitize the innate compartment to detection of PAMPs. Specifically, expression of the innate sensors Nod1, Nlrp3, Aim2, Ctsb, Tlr1/2/4/7/8/9, Oas2, and RhoA increased following thermal ablation (Fig. 4B). However, thermal ablation also upregulates sterile inflammatory pathway genes including, Nlrp3, Jun, Mefv, and Il6. Il1β expression is enhanced with ablation and is associated with Th1 based cellular immunity via the caspase-1 pathway, Il1α can act as a DAMP and stimulate TLR pathways³³. Notably, ablation alone did not alter innate immune expression profiles of distant tumors, which remained similar to control tumors.

The combined ablation-immunotherapy protocol produces a broad effect that is not limited to a particular receptor or cell-type; it upregulates innate receptors across multiple families and alters the polarization of infiltrating macrophages, thereby facilitating a strong systemic immune response. Addition of the innate agonist CpG to the ablation protocol further upregulates Nod2, Oas2, RhoA, Pycard, Tlr1/2 and Il12 compared to ablation alone, augmenting the innate immune response and further enriching the following KEGG^34,35 (developed by Kanehisa Laboratories) innate immune sensing and signaling pathways: NLR (KEGG: mmu04621), RIG-I (KEGG: mmu04622), cGAS cytosolic DNA sensing (KEGG: mmu04623), and TLR (KEGG: mmu04620). Innate immune response (GO:0045087) was the second most enriched pathway for AI-T, with expression of 111 genes altered, while response to lipopolysaccharide (LPS) (GO:0032496) was the second most enriched for both thermal and mechanical ablation, including 70 and 32 genes, respectively (Supplementary Table S2). Moreover, the combined immunotherapy-ablation protocol did not increase expression of Nlrp3 and other inflammation-associated genes. Additionally, ligation of TLR9 by CpG potentiates macrophage phagocytosis³⁶ and may synergize with ablation-induced DAMP release. For example, increased Adgre1 expression in AI-T tumors is consistent with enhanced macrophage infiltration into the treated tumors.

To investigate the cellular effects of ablation, with or without immunotherapy, we homogenized tumor tissues and analyzed the cellular components via flow cytometry. Both thermal (Fig. 4C) and mechanical (Fig. 4D) HIFU increased dendritic cell (DC) infiltration when added to an immunotherapy protocol. However, infiltration of DCs was not significantly altered in tumors treated with either mechanical HIFU (Fig. 4E) or thermal ablation (Fig. 4F) in the absence of immunotherapy. Although we expect increased infiltration of antigen presenting cells following an ablative treatment, it is possible our choices of time points for flow cytometry did not capture the infiltration. Expression of MHCII was consistent with increased DC infiltration at distant sites in animals that received both thermal ablation and immunotherapy (Fig. 4G).

At the 1-week timepoint, although macrophage numbers (as assessed by flow cytometry) were similar across treatments, the phenotype was distinct. Siglec1 is expressed on macrophages and DCs associated with the presentation of dead cell antigens. HIFU upregulated Siglec1 (Cd169) expression within directly-treated tumors, with AI-T > A-T > M-T (Fig. 4H). Expression of Cd43, the T-cell counterreceptor for Siglec1, increased only with the combined protocol as compared to control (Fig. 4I). Tumor macrophage infiltration was unchanged following either mechanical HIFU (Fig. 4J) or thermal ablation (Fig. 4K) as quantified by flow cytometry. Macrophage polarization was altered by both immunotherapy and mechanical HIFU, with distinct populations developing from 24 h to 7 days post HIFU (Supplemental Fig. S4). Macrophage polarization shifted towards M2 in tumors treated with mechanical HIFU (Supplementary Fig. S5A). With the addition of immunotherapy to thermal ablation, M2 macrophages were reduced compared to control or immunotherapy alone. Furthermore, the combined mechanical HIFU-immunotherapy protocol altered macrophage polarization to predominantly non-M1/M2, compared to control or mechanical HIFU alone (Supplementary Fig. S5B). IHC sections stained for F4/80 were concordant with flow cytometry results: macrophage infiltration was sparse for thermally-ablated tumors at the 1-week time point (Supplementary Fig. S5C); however, addition of immunotherapy increased macrophage recruitment and infiltration (Supplementary Fig. S5D). Mechanical HIFU did not alter macrophage infiltration (Supplementary Fig. S5E), and neither form of ablation was sufficient to enhance macrophage populations compared to control tumors (Supplementary Fig. S5F).

Adaptive immune response is augmented by the addition of immunotherapy

In addition to macrophages, we hypothesized that tumor-infiltrating T cells increase as a result of ablative treatment. To test this hypothesis, we first investigated T cell activation markers within the transcriptome and found that T cell activation markers and adaptive immune genes were upregulated by thermal ablation in treated, but not distant tumors (Fig. 5A), which was expected of a localized treatment modality. The combined immunotherapy-ablation protocol further upregulated T cell activation markers such as Gzma, Gzmb, and Batf3 in the directly-treated tumor compared to ablation alone. Additionally, the combined protocol extended the effect to distant sites, upregulating expression of adaptive immune genes such as Eomes, Cd8a, Prf1, and Ifng in distant tumors, extending the effect of localized ablation treatment to systemic effects. T cell and adaptive response pathways were likewise upregulated by the combination of ablation and immunotherapy (Supplementary Table S3) based on gene ontology enrichment analysis by DAVID. Mechanical HIFU did not alter adaptive immune gene expression compared to control. Based on the RNA-seq signature, macrophages were the predominant immune cell infiltrating tumors (Fig. 5B) but were underestimated by flow cytometry (Supplementary Fig. S6). Thermal ablation with immunotherapy increased CD8⁺ T cell infiltration into distant tumors (AI-C) (Fig. 5C).

CD8⁺ T cell infiltration at distant sites increased under the combined immunotherapy-HIFU protocols compared to the no treatment control cohort (Fig. 5D). However, only thermal ablation primed with immunotherapy improved CD8⁺ T cell recruitment compared to immunotherapy alone, indicating the possibility of synergistic effects from simultaneous DAMP and PAMP signaling. The thermal ablation-immunotherapy protocol also increased CD4⁺ T cell infiltration at both directly-treated and distant sites (Fig. 5E). The impact of immunotherapy on the directly-treated tumor is further shown in Supplementary Fig. S7.

In addition to effector tumor infiltrating lymphocytes, immunosuppressive cell populations were induced by HIFU but in distinct patterns for thermal ablation compared to mechanical HIFU. Thermal ablation expanded immunosuppressive populations at distant sites while mechanical HIFU primarily enhanced these populations within directly-treated tumors. Tregs were enhanced in both directly-treated and distant tumors following thermally ablative therapies, while mechanical HIFU increased Tregs in directly-treated tumors (Fig. 5C). Both thermal ablation (Fig. 5F) and mechanical HIFU (Fig. 5G) increased the number of myeloid-derived suppressor cells (MDSCs) within the treated tumor compared to control (1.9-fold greater for thermal and 2.3-fold for mechanical), but the effect was only significant for mechanical ablation. Thermal ablation generated an ~ threefold increase in MDSC infiltration at distant tumor sites, while mechanical ablation did not enhance distant MDSC infiltration.

Pdcd1 (PD-1) gene expression was decreased in HIFU monotherapy cohorts compared to control, but the combined immunotherapy-ablation protocol increased Pdcd1 expression at both directly-treated and distant sites (Fig. 5A). On IHC, control tumors exhibited relatively scattered expression of Pdcd1 (Fig. 5H); thermal ablation reduced PD-1⁺ cells in both the directly-treated (Fig. 5I) and distant (Fig. 5J) tumors. Administration of CpG and anti-PD-1 followed by thermal ablation increased the number and infiltration of PD-1⁺ cells in both the directly-treated (Fig. 5K) and distant tumors (Fig. 5L). Pdcd1 expression was higher in AI-C as compared to AI-T tumors. Tumors treated with mechanical HIFU had similar PD-1 expression (Fig. 5M) as thermally-ablated tumors.

Foxp3 is primarily expressed on immunosuppressive CD4⁺ T cells, which suppress the activity of immune effector cells while simultaneously attenuating the inflammatory response. Treg infiltration has been associated with worse outcomes in some types of cancers (e.g. melanoma, breast, and ovarian) but paradoxically with improved tumor control in other cancers (e.g. head and neck and bladder cancers). Foxp3 expression, as observed on IHC, was induced by thermal ablation but reduced under the combined immunotherapy-ablation protocol. Control tumors had minimal Foxp3 expression (Fig. 5N) which was increased after thermal ablation in the directly-treated (Fig. 5O) and distant tumors (Fig. 5P). Priming with immunotherapy reduced the number of Foxp3⁺ cells in the directly-treated tumor (Fig. 5Q) compared to thermal ablation alone but marginally increased Foxp3⁺ cells in the distant tumor (Fig. 5R). Few Foxp3⁺ cells were observed following mechanical HIFU (Fig. 5S).

The combination of ablation and immunotherapy highly enriched multiple T cell activation, signaling and proliferation pathways (Supplementary Table S3). Thermal ablation enriched similar pathways to a lesser degree.

Inflammatory cytokines and chemokines are upregulated by HIFU

To further understand how the combination of HIFU and immunotherapy extended the overall immune stimulatory effect towards systemic lymphocyte activation and trafficking, we evaluated the genetic expression levels of inflammatory cytokines and chemokines. Across the combined immunotherapy-ablation treatment, mechanical HIFU and thermal ablation, T-cell chemotactic signaling chemokines and cytokines were enhanced in the treated tumors (Fig. 6A). The most enriched pathway across all treatment groups was neutrophil chemotaxis (GO:0030593) (Supplementary Table S4), which was reflected by strong upregulation of Cxcl9, the leukocyte chemotaxis signaling molecule, in the distant tumors of animals treated with the combined ablation-immunotherapy protocol (Fig. 6B). Chemokines related to leukocyte recruitment were upregulated in A-T tumors, which were largely distinct from those upregulated in M-T tumors. Thermal ablation upregulated Cxcl10, preferentially at distant tumor sites (Fig. 6C).

Combining TLR9 ligation with ablation enhances the recruitment of B220⁺ cells

Since Il1β is related to the induction of humoral immunity³⁷, we hypothesized that plasma-based immune cell gene expressions varied across treatment cohorts. To test this hypothesis, we first interrogated our mouse transcriptome for humoral-related gene expression and found that it was enhanced by thermal ablation as compared to control (Fig. 7A) in the directly-treated tumors. The addition of immunotherapy further increased genes associated with plasma-based immune cell receptor signaling in both the directly-treated and distant tumors (Supplementary Table S5). Thermal ablation upregulated Siglech in directly-treated tumors (Supplementary Fig. S8) compared to both control and AI-T tumors, while Cd20 expression was significantly higher in AI-C tumors (Supplementary Fig. S8). IHC confirmed B220⁺ cell infiltration resulting from thermal ablation (Fig. 7B,F) and the further enhancement from priming with immunotherapy (Fig. 7C,G). Mechanical HIFU did not enhance B220⁺ cell infiltration (Fig. 7D,H) as compared with control tumors (Fig. 7E,I). On magnified sections, B220⁺ cells were localized to the rim of thermally-ablated tumors (Fig. 7F), but infiltrate within AI-T tumors (Fig. 7G). B220⁺ cells were not detected in tumors treated with mechanical HIFU (Fig. 7H) and control tumors (Fig. 7I) on magnified sections.

Discussion

We have previously reported on the efficacy of a combined ablation-immunotherapy protocol in generating an adaptive immune response and reducing systemic tumor burden¹⁵. Here, we probe the basis of that response by examining innate immune receptors and cytokines that initiate the adaptive response and quantifying the immunomodulatory effects of thermal and mechanical HIFU alone. Thus, we can directly compare across two HIFU parameter sets and the effects of incorporating an immune agonist.

Thermal ablation has a broad effect on innate immune sensors and pathways. However, by itself, thermal ablation primarily releases DAMPs and produces a strong, sterile inflammatory response. Priming with an immune agonist attenuates Il6 gene expression and enhances upregulation of innate receptors across multiple families, sensitizing the innate immune system to danger signals. When combined with innate immune agonists, tumor specific antigens released by ablation are processed by antigen presenting cells (APCs) such as macrophages and dendritic cells and cross-presented to T-cells, generating a systemic anti-tumor response. The strength and breadth of the innate immune response to HIFU offers multiple opportunities for synergy with agonists.

HIFU ablation of solid tumors is attractive due to its non-invasive nature and limited side-effect profile. Moreover, HIFU can be repeated as needed without regard to lifetime dose limits, unlike chemotherapy or radiation therapy. However, while HIFU can rapidly kill all viable tumor cells within a treated volume, as a monotherapy, HIFU is insufficient to generate a robust systemic anti-tumor response. In the mouse models we have studied, ablation alone does not result in any significant increase in survival²¹. Indeed, ablation releases DAMPs from damaged cells, recruiting neutrophils and monocytes to the ablated site, but in the absence of additional stimuli for pattern recognition receptors, DAMPs will be cleared and the leukocytes recruited by the initial inflammatory insult transition from acute inflammatory to reparative, or even generate a pro-tumor chronic inflammatory state via alteration in macrophage metabolism³⁸. Nevertheless, both mechanical HIFU and thermal ablation can synergize with immunotherapy due to their immunomodulatory effects, both stimulatory and immunosuppressive at varying stages, and their ability to release substantial tumor-specific antigen. Ultimately, the efficacy of a combined ablation-immunotherapy protocol is highly dependent on the immunomodulatory agents and the timing of their administration relative to ablation.

We have previously reported on the importance of treatment sequence when incorporating ablation into an immunotherapy protocol²¹, where we found that a combined ablation-immunotherapy regimen was less effective than immunotherapy alone when the two treatments began simultaneously. However, the combined protocol significantly enhanced survival compared to immunotherapy alone when immunotherapy was initiated prior to ablation. In this study we found ablation broadly increases PRR expression across multiple sensor types. This would seem to suggest addition of innate agonists after ablation would result in increased survival. However, we hypothesize that the large increase in Il6 expression and the resulting sterile inflammatory response following ablation alters expression in PRR pathways. When pre-treated with immunotherapy, the increase in Il6 expression is smaller following ablation (compared to ablation alone). We hypothesize this attenuation of Il6 expression is one important component of generating an optimal anti-tumor immune response in the model studied here. Priming with immunotherapy substantially reduces the ablation-induced spike in Il6 expression, perhaps via reducing expression of upstream Il1β. Furthermore, some DAMPs released by ablation, such as oxidized phosphorylcholine derivatives, can drive dendritic cells pre-exposed to TLR ligands into a hyperactive activation state, enhancing their antigen presentation capabilities^39,40. Conversely, exposure to the same DAMPs prior to TLR ligation attenuated the inflammatory activity of DCs upon subsequent TLR ligation⁴¹, desensitized DCs to CD40-CD40L interaction and may limit migration of DCs to lymph nodes⁴².

Pre-treatment with innate immune agonists may slow the clearance of cellular debris generated by both thermal ablation or mechanical HIFU via generation of an inflammatory environment prior to and during the clearance process⁴³. Thus, tumor specific antigen released by ablation may remain in situ and available for uptake by APCs longer when previously primed with immunotherapy. Aside from providing a reservoir of tumor antigen, ablation may also free amino and nucleic acids capable of inducing TLR-independent interferon production^44,45,46,47, which may provide an additional mechanism for synergy. TLR tolerance can result in reduced type-1 IFN release after repeated ligation⁴⁸, thus an alternative mechanism for type-1 IFN release, such as via HIFU, may be beneficial in those cases.

Thermal ablation increased MDSCs as a fraction of leukocytes within both the treated and distant tumors. Although Il1β can generate an adaptive anti-tumor immune response, following ablation it may produce a pro-tumor inflammatory response that ultimately results in MDSC recruitment and induction^49,50. Furthermore, Il1β induces Il6 production, which promotes sustained, excessive inflammation, fostering a pro-tumor immune environment. Mechanical HIFU caused a similar increase in MDSCs in the directly-treated tumor, but mechanical HIFU did not increase MDSC expansion at distant sites. It has been suggested that release of tumor-derived microRNAs can induce MDSCs⁵¹, which may be one mechanism for the increase in MDSCs following destruction of tumor cells by mechanical HIFU. Furthermore, in the absence of immune adjuvants or co-stimulatory factors such as CD40 agonists, ablation alone may result in sub-optimal T cell priming and decreased T cell expansion^52,53, and has been shown to induce resistance to checkpoint modulation⁵⁴.

In addition to increased MDSC infiltration, thermal ablation increased the number and infiltration of FoxP3⁺ cells at both directly-treated and distant sites. Mechanical HIFU caused a slight decrease in FoxP3⁺ cells compared to control. Thermal and mechanical HIFU both reduced the number of PD-1⁺ cells on IHC compared to control, while pre-treatment with immunotherapy substantially increased PD-1⁺ cells in both directly-treated and distant tumors. The increase in PD-1 expression may relate to decreased Il6 expression in cohorts that received immunotherapy given that Il6 blockade has been demonstrated to upregulate PD-L1⁵⁵. We hypothesize that increased Pdcd1 expression increases sensitivity to anti-PD-1 therapy, which was a component of the ablation-immunotherapy combined treatment protocol employed in this study.

Our study has multiple limitations. Based on flow cytometry, we chose to evaluate gene expression at 7 days following thermal ablation and 1 day after mechanical ablation. Our reasoning is that cell death from thermal ablation occurred over several days. Alternatively, the mechanical ablation resulted in immediate cell death with healing and tumor regrowth observed within a few days.

In our work, the same system (an MRgFUS system with a center frequency of 3 MHz) was used for all thermal and mechanical protocols. For thermal ablation, continuous wave insonation was applied at 3.1 MPa peak negative pressure (PNP), whereas for mechanical ablation a 16.9 MPa PNP was applied and the duty cycle was minimized such that any temperature increase was less than 2 °C. The resulting studies provide a starting point for the comparison of regimes; however, in the future, comparisons to the histotripsy protocols used in^56,57,58,59 must also be made directly with the dedicated lower frequency equipment required for those studies. A disadvantage of such equipment (limiting a comparative study of thermal and mechanical mechanisms) is the larger insonified volume (relative to mouse tumors) and the incompatibility with MR guidance. The use of MRgFUS for the thermal ablation protocols is essential in order to quantify the temperature increase, assuring that a minimum temperature of 60 °C is achieved and the entire tumor is not exposed to this temperature for an extended period of time. Such an exposure can heat fix the tissue and reduce the immune effect. Such parameters will be the subject of future study. Finally, with bulk RNA-seq, we are unable to determine if enhanced/reduced transcript levels are the result of upregulation of expression of resident cells within the tumor or greater infiltration of cells constitutively expressing the examined genes at differing levels than tumor cells.

Overall, we found that both mechanical HIFU and thermal ablation rapidly kills viable tumor, with corresponding reductions in cancer and proliferation-related genes and pathways. Although HIFU does not appreciably activate an adaptive immune response, it induced inflammation via upregulation of Il1β and Il6, which leads to recruitment of MDSCs and subsequent tumor regrowth. Pre-treatment with an innate immune agonist, CpG in this case, and modulation of the PD-1/PD-L1 axis attenuates the excessive inflammatory response induced by ablation and overcomes rebound immune suppression.

Methods

Reagents

The TLR9 agonist, CpG 1826 (class B), was purchased from InvivoGen (San Diego, CA). Anti-mouse PD-1 (CD279) mAb (rat IgG2a, clone RMP1-14), was purchased from Bio X Cell (West Lebanon, NH).

Ethics statement

All experiments and methods were performed in accordance with relevant guidelines and regulations. Specifically, all animal experiments were conducted with approval from the Stanford University Administrative Panel on Laboratory Animal Care (APLAC) and the University of California, Davis, (Davis, California) Institutional Animal Care and Use Committee (IACUC).

Animal model

The neu exon deletion line (NDL), a syngeneic model of mammary adenocarcinoma, was obtained from the Alexander Borowsky Laboratory (UC Davis, Davis, CA). Four-week-old FVB/n female mice purchased from Charles River (Wilmington, MA) were transplanted with NDL tumor biopsies (~ 1 mm³) bilaterally into the fourth and ninth inguinal mammary fat pads. Approximately 21 days later, when tumors reached ~ 4 mm in longest dimension, mice were randomized into treatment groups.

Therapeutic and sequencing protocols

A total of 167 mice were studied. For all studies, one tumor was directly-treated with either: immunotherapy alone (I-T), thermal ablation alone (A-T), mechanical ablation alone (M-T), thermal ablation primed with immunotherapy (AI-T), or mechanical ablation primed with immunotherapy (MI-T). RNA-seq data for the combination of thermal ablation and immunotherapy were reported previously in¹⁵ and are included for comparison. Distant (contralateral) tumors were also examined for each cohort (I-C, A-C, M-C, AI-C, and MI-C respectively). For flow cytometry quantification of immune cell populations, tumors were collected at 72 h and 7 days post thermal ablation (with and without immunotherapy), and at 24 h, 72 h, and 7 days post mechanical ablation. These studies provided the basis for the timing of the following sequencing studies. For RNA-seq studies, tumors were collected and sequenced 7 days following thermal ablation (with and without immunotherapy), or 24 h after mechanical HIFU. Our choice of the 24-h timepoint for mechanical HIFU was based on the ultrasound parameters we chose to employ and our experience with their biological effects and immunohistochemistry studies that we performed at 2 h, 24 h, 48 h, 72 h, and 7 days post mechanical ablation (Supplementary Fig. S1 showing tissue morphology on H&E stained sections); namely, cell death/destruction was significantly more rapid with mechanical ablation than with thermal ablation. Following thermal ablation, cells slightly distanced from the focal zone receive a lethal thermal dose due to heat diffusion but may take tens of hours to die. Consequently, cell debris is more rapidly released after mechanical ablation than after thermal ablation.

Animals receiving ablation and immunotherapy were primed with immunotherapy prior to ablation as follows: 100 μg CpG injected intratumorally (i.t.) on days 21, 24, 28, and 30; 200 μg anti-PD-1 injected intraperitoneally (i.p.) on days 21, 28, and 34, and ablation was performed on day 30.

MRgFUS ablation protocol

All ablations were performed under MR guidance on a Bruker BioSpec 7T small animal MR system (Bruker Biospin) with core body temperature monitoring using a 16-element annular array transducer operating at 3 MHz (Imasonic SAS)⁶⁰. Acoustic pressure was calibrated with a fiber optic hydrophone (HFO690, Onda Corp.) in a degassed water tank under free-field conditions. Prior to ablation, mice were given 0.05–0.1 mg/kg buprenorphine subcutaneously (s.c.) and 0.05 mmol/kg gadoteridol (Bracco Imaging) i.p. and imaged with a T1w RARE (TE/TR = 12.5/750 ms, FOV = 4.8 cm × 4.8 cm, MTX = 256 × 256, ST/SI = 1/1 mm, 17 slices) sequence for tumor localization and treatment planning. Tumors were then ablated and temperature was monitored in real time via the MR proton resonance frequency shift using Thermoguide Software (Image Guided Therapy), with α = − 0.0101 ppm/°C, TE/TR = 4.5/21 ms⁵.

For thermal ablation, continuous wave (CW) insonation was employed at 3.1 MPa in a circular pattern (diameter of 2 mm, scan speed of 1 revolution per second) until the targeted volume reached at least 60 °C and a thermal dose in cumulative equivalent minutes at 43 degrees (CEM43) of more than 5000 was achieved¹⁵.

For mechanical ablation, tumors were ablated in a grid pattern with individual points separated by 0.5 mm, treating the entirety of the tumor with the exception of a ~ 0.5–0.75 mm rim. Each point of the grid pattern was insonated for 5 ms at 16.9 MPa (peak negative pressure), 95.3 MPa (peak positive pressure), with a duty cycle of 0.5%. The duty cycle was chosen to minimize thermal effects and maintain temperature elevations ≤ 2 °C. The grid pattern was repeated 10 times.

Histology and immunohistochemistry

Tissues for microscopic analysis were prepared as previously described¹⁵. Briefly, tissues were fixed in 10% buffered formalin for 24 h or longer, followed with 70% ethanol overnight before further processing. A Tissue-Tek VIP autoprocessor (Sakura, Torrance, CA) was used to process samples for paraffin-embedding. Tissue blocks were then sectioned to 4 μm, and the sections were mounted on glass slides. Hematoxylin and Eosin (H&E) staining was performed at the University of California, Davis, Department of Pathology and Laboratory Medicine. For staining of cytotoxic T cells, the rat anti-mouse CD8a primary antibody was used (1:500; 14-0808, eBiosciences). All IHC was performed manually without the use of an automated immunostainer and using the ABC method. Antigen retrieval was performed using a Decloaking Chamber (Biocare Medical, Concord, CA) with citrate buffer at pH 6.0, 125 °C and pressurized to 15 psi. The total time slides were in the chamber was 45 min. Incubation with the primary antibody was performed at ambient temperature overnight in a humidified chamber. Normal goat serum was used for blocking. Biotinylated goat anti-rat (1:500; Vector Labs, Burlingame, CA) was the secondary antibody used with a Vectastain ABC Kit Elite and a Peroxidase Substrate Kit DAB (both from Vector Labs) used for the amplification and visualization of signal, respectively. Mouse spleen was used as a positive control. Stained slides were scanned on an AT2 Scanscope (Leica Biosystems) and digital images viewed using the ImageScope program (Leica Biosystems).

Flow cytometry antibodies

Flow cytometry was performed with mouse-specific fluorochrome-conjugated monoclonal antibodies (mAbs) as previously described¹⁵. Briefly, Pacific blue (PB)-anti-CD45 (30-F11), fluorescein isothiocyanate (FITC)-anti-F4/80 (BM8), phycoerythrin (PE)-Cy7-anti-CD11c (N418), PE-Cy7-anti-CD3 (145-2C11), allophycocyanin (APC)-Cy7-anti-CD11b (M1/70), Alexa Fluor (AF)-700-anti-Ly6G/Ly6C (Gr-1, RB6-8C5), and AF700-anti-CD8a (53–6.7) were purchased from BioLegend (San Diego, CA). FITC-anti-CD4 (GK1.5) and APC-anti-CD206 (C068C2) were purchased from BD Biosciences (San Jose, CA). PE-Cy5-anti-MHCII (M5/114.15.2) and PE-anti-CD86 (GL1) were purchased from eBioscience (San Diego, CA). Negative fluorescence-minus-one (FMO) control staining was performed with isotype-matched mouse, rat or hamster IgG mAbs, and nonspecific binding was blocked with the FcγIII/II receptor-mediated anti-CD16/CD32 antibody (2.4G2) from BD Biosciences.

Flow cytometry

Mice were sacrificed at indicated time points following treatment, and tissues (tumors, tumor draining lymph nodes, spleen, and blood) were collected and processed by mechanical and enzymatic disruption to single-cell suspensions for immune cell profiling via flow cytometry. Briefly, after euthanizing tumor-burdened mice, bilateral tumors were dissected and stored in 1 mL of DMEM on ice. For mechanical disruption, tumors were then finely minced using dissection scissors. Following this, the tumor solution was enzymatically digested by adding 40 μL of Liberase DL (28 U/mL stock), 80 μL of Liberase TL (14 U/mL stock), and 40 μL of DNase I (15 mg/mL stock) and incubated at 37 °C under continuous rotation⁶¹. After incubation, the digested tumor solution was triturated, passed through a 70 μm cell strainer, and washed with several mL of DMEM + 10% FBS. The solution was spun at 300G for 5 min at 4 °C and resuspended in 1 mL of fresh DMEM + 10% FBS for subsequent cell counting and flow cytometry antibody staining.

Live-dead cell staining with the LIVE/DEAD Fixable Aqua Dead Cell Stain Kit (Invitrogen, Carlsbad, CA) was performed according to the manufacturer’s instructions prior to all other antibody staining in order to exclude dead cells from analysis. Antibody panel combinations used to distinguish immune cell populations were CD45⁺ (leukocytes) plus the following: CD11b⁺, F4/80⁺, Gr-1⁻ (macrophages); CD11b⁺, F4/80⁺, Gr-1⁻, CD86⁺, CD206⁻ (M1 macrophages); CD11b⁺, F4/80⁺, Gr-1⁻, CD86⁻, CD206⁺ (M2 macrophages); CD11b⁺, F4/80⁺, Gr-1⁻, CD86⁺, CD206⁺ (M1/M2 macrophages); CD11b⁺, F4/80⁺, Gr-1⁻, CD86⁻, CD206⁻ (non-M1/M2 macrophages); CD11c⁺, MHCII⁺, F4/80⁻ (dendritic cells); CD11b⁺, Ly6G/Ly6C⁺ (myeloid derived suppressor cells); CD3⁺ (T-cells); CD3⁺, CD4⁺ (CD4⁺ T-cells); CD3⁺, CD8⁺ (CD8⁺ T-cells). Cell samples were fixed in Cytofix buffer (BD Biosciences), diluted to 1% paraformaldehyde (PFA) in PBS−/−, and run within 24 h on an LSRII flow cytometer (BD Biosciences). All data were analyzed using FlowJo v10 software (TreeStar).

RNA isolation

Snap-frozen NDL tumors were submitted to the UC Davis Comprehensive Cancer Center’s Genomics Shared Resource for isolation of total cellular RNA and subsequent RNA-seq analysis. Total cellular RNA was isolated using the TRIzol Reagent (Invitrogen) using a modified protocol that incorporates an additional extraction with phenol/chloroform/isoamyl alcohol (25:24:1, pH 4.3) followed by clean-up with an RNeasy spin column (Qiagen). RNA concentration and purity were assessed with a NanoDrop 1000 Spectrophotometer (Thermo Scientific) and quality assessments were made using an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA).

Directional RNA-seq library preparation and next-generation sequencing

Indexed RNA-seq libraries were prepared using the NEBNext Ultra Directional RNA Library Prep Kit (New England BioLabs, Ipswich, MA), according to the manufacturer’s standard protocol as previously described¹⁵. Briefly, poly(A) mRNA was captured from total RNA (100 ng) by binding to magnetic oligo(dT)₂₅ beads, fragmented, and then double-stranded cDNA generated by random-primed first-strand synthesis and second strand synthesis in the presence of dUTP for strand marking^62,63. The double-stranded cDNA was then end repaired and 3′-dA tailed, followed by ligation of an Illumina-compatible adaptor and subsequent uridine excision with the USER (Uracil-Specific Excision Reagent) enzyme. The libraries were then indexed and enriched by high-fidelity PCR amplification (15 cycles) with i7 primer/index primers and Q5 High-Fidelity DNA Polymerase. Subsequently, libraries were pooled for multiplex sequencing on an Illumina HiSeq 4000 System (150-bp, paired-end, ~ 25–30 × 10⁶ reads per sample).

NGS data processing and analysis

RNA-seq data was analyzed using a STAR-StringTie-Cufflinks pipeline. Raw sequence reads (FASTQ format) were mapped to the reference mouse genome assembly (GENCODE, GRCm38, release 05/2017) using STAR (Spliced Transcripts Alignment to a Reference) software with a 2-pass alignment approach with an initial alignment to detect novel junctions and to insert these junctions into the genome index, followed by a second pass to re-align reads using both annotated (GENCODE, Release M14, GRCm38.p5) and novel junctions⁶⁴. Mapped reads were then passed onto StringTie for transcript assembly⁶⁵. Subsequently, transcript-level expression was quantified with Cufflinks tools (e.g. Cuffquant, Cuffnorm) to yield normalized expression as FPKM (fragments per kilobase of transcript per million fragments mapped reads) and test for differential expression (Cuffdiff)⁶⁶.

Cuffnorm files were imported to MATLAB for processing and analysis. Z-scores were calculated using raw FPKM values across all presented samples for a single gene in R. Principal component analysis was performed in R. Gene set enrichment analysis was performed on the set of all differentially expressed genes for each Cuffdiff analysis (i.e. genes with a differential expression at p < 0.01) using DAVID’s Functional Annotation Tool. Specifically, the GO Biological Processes and the KEGG Pathways were investigated⁶⁷.

CIBERSORTx

Bulk RNA-seq data was input to CIBERSORTx⁶⁸ at https://cibersortx.stanford.edu, which imputed immune cell fractions.

Statistical analysis

Statistical analyses were performed using Prism 8 software (GraphPad Software Inc). Results are presented as mean ± SD, unless otherwise indicated. One-way ANOVA was performed for all analyses of three or more groups followed by a with Tukey correction for multiple hypotheses in GraphPad Prism. Analysis of differences between two groups was performed using an unpaired t-test assuming unequal variance. p values less than 0.05 were considered significant.

References

Ning, Z. et al. HIFU is safe, effective, and feasible in pancreatic cancer patients: A monocentric retrospective study among 523 patients. Onco Targets Ther. 12, 1021–1029. https://doi.org/10.2147/OTT.S185424 (2019).
Article CAS PubMed PubMed Central Google Scholar
Strunk, H. M. et al. Clinical use of high-intensity focused ultrasound (HIFU) for tumor and pain reduction in advanced pancreatic cancer. Rofo 188, 662–670. https://doi.org/10.1055/s-0042-105517 (2016).
Article CAS PubMed Google Scholar
Khokhlova, T. D. & Hwang, J. H. HIFU for palliative treatment of pancreatic cancer. Adv Exp. Med. Biol. 880, 83–95. https://doi.org/10.1007/978-3-319-22536-4_5 (2016).
Article CAS PubMed Google Scholar
de Senneville, B. D., Moonen, C. & Ries, M. MRI-guided HIFU methods for the ablation of liver and renal cancers. Adv. Exp. Med. Biol. 880, 43–63. https://doi.org/10.1007/978-3-319-22536-4_3 (2016).
Article CAS PubMed Google Scholar
Fite, B. Z. et al. Magnetic resonance thermometry at 7T for real-time monitoring and correction of ultrasound induced mild hyperthermia. PLoS ONE 7, e35509. https://doi.org/10.1371/journal.pone.0035509 (2012).
Article ADS CAS PubMed PubMed Central Google Scholar
Kohler, M. O. et al. Volumetric HIFU ablation under 3D guidance of rapid MRI thermometry. Med. Phys. 36, 3521–3535. https://doi.org/10.1118/1.3152112 (2009).
Article PubMed Google Scholar
Gaur, P. & Grissom, W. A. Accelerated MRI thermometry by direct estimation of temperature from undersampled k-space data. Magn. Reson. Med. 73, 1914–1925. https://doi.org/10.1002/mrm.25327 (2015).
Article PubMed Google Scholar
Rieke, V. & Butts Pauly, K. MR thermometry. J. Magn. Reson. Imaging 27, 376–390. https://doi.org/10.1002/jmri.21265 (2008).
Article PubMed PubMed Central Google Scholar
Maxwell, A. D. et al. Cavitation clouds created by shock scattering from bubbles during histotripsy. J. Acoust. Soc. Am. 130, 1888–1898. https://doi.org/10.1121/1.3625239 (2011).
Article ADS PubMed PubMed Central Google Scholar
Simon, J. C. et al. Miniature acoustic fountain mechanism for tissue emulsification during millisecond boiling in high intensity focused ultrasound fields. J. Acoust. Soc. Am. 129, 2478–2478. https://doi.org/10.1121/1.3588151 (2011).
Article ADS Google Scholar
Hall, T. L. et al. Histotripsy of rabbit renal tissue in vivo: Temporal histologic trends. J. Endourol. 21, 1159–1166. https://doi.org/10.1089/end.2007.9915 (2007).
Article PubMed Google Scholar
Hoogenboom, M. et al. In vivo MR guided boiling histotripsy in a mouse tumor model evaluated by MRI and histopathology. NMR Biomed. 29, 721–731. https://doi.org/10.1002/nbm.3520 (2016).
Article PubMed Google Scholar
Lu, P. et al. Increased infiltration of activated tumor-infiltrating lymphocytes after high intensity focused ultrasound ablation of human breast cancer. Surgery 145, 286–293. https://doi.org/10.1016/j.surg.2008.10.010 (2009).
Article PubMed Google Scholar
Hu, Z. et al. Release of endogenous danger signals from HIFU-treated tumor cells and their stimulatory effects on APCs. Biochem. Biophys. Res. Commun. 335, 124–131. https://doi.org/10.1016/j.bbrc.2005.07.071 (2005).
Article CAS PubMed PubMed Central Google Scholar
Chavez, M. et al. Distinct immune signatures in directly treated and distant tumors result from TLR adjuvants and focal ablation. Theranostics 8, 3611–3628. https://doi.org/10.7150/thno.25613 (2018).
Article CAS PubMed PubMed Central Google Scholar
van den Bijgaart, R. J. et al. Thermal and mechanical high-intensity focused ultrasound: Perspectives on tumor ablation, immune effects and combination strategies. Cancer Immunol. Immunother. 66, 247–258. https://doi.org/10.1007/s00262-016-1891-9 (2017).
Article PubMed Google Scholar
Reits, E. A. et al. Radiation modulates the peptide repertoire, enhances MHC class I expression, and induces successful antitumor immunotherapy. J. Exp. Med. 203, 1259–1271. https://doi.org/10.1084/jem.20052494 (2006).
Article CAS PubMed PubMed Central Google Scholar
Takahashi, Y. et al. Optimized magnitude of cryosurgery facilitating anti-tumor immunoreaction in a mouse model of Lewis lung cancer. Cancer Immunol. Immunother. 65, 973–982. https://doi.org/10.1007/s00262-016-1858-x (2016).
Article CAS PubMed Google Scholar
Kroemer, G., Galluzzi, L., Kepp, O. & Zitvogel, L. Immunogenic cell death in cancer therapy. Annu. Rev. Immunol. 31, 51–72. https://doi.org/10.1146/annurev-immunol-032712-100008 (2013).
Article CAS PubMed Google Scholar
Wu, F. et al. Activated anti-tumor immunity in cancer patients after high intensity focused ultrasound ablation. Ultrasound Med. Biol. 30, 1217–1222. https://doi.org/10.1016/j.ultrasmedbio.2004.08.003 (2004).
Article PubMed Google Scholar
Silvestrini, M. T. et al. Priming is key to effective incorporation of image-guided thermal ablation into immunotherapy protocols. JCI Insight 2, e90521. https://doi.org/10.1172/jci.insight.90521 (2017).
Article PubMed PubMed Central Google Scholar
Iwasaki, A. & Medzhitov, R. Regulation of adaptive immunity by the innate immune system. Science 327, 291–295. https://doi.org/10.1126/science.1183021 (2010).
Article ADS CAS PubMed PubMed Central Google Scholar
Nathan, C. Points of control in inflammation. Nature 420, 846–852. https://doi.org/10.1038/nature01320 (2002).
Article ADS CAS PubMed Google Scholar
Messenheimer, D. J. et al. Timing of PD-1 blockade is critical to effective combination immunotherapy with anti-OX40. Clin. Cancer Res. 23, 6165–6177. https://doi.org/10.1158/1078-0432.CCR-16-2677 (2017).
Article CAS PubMed PubMed Central Google Scholar
Hotz, C. & Bourquin, C. Systemic cancer immunotherapy with Toll-like receptor 7 agonists: Timing is everything. Oncoimmunology 1, 227–228. https://doi.org/10.4161/onci.1.2.18169 (2012).
Article CAS PubMed PubMed Central Google Scholar
Sato-Kaneko, F. et al. Combination immunotherapy with TLR agonists and checkpoint inhibitors suppresses head and neck cancer. JCI Insight. https://doi.org/10.1172/jci.insight.93397 (2017).
Article PubMed PubMed Central Google Scholar
Lundt, J. E. et al. Non-invasive, rapid ablation of tissue volume using histotripsy. Ultrasound Med. Biol. 43, 2834–2847. https://doi.org/10.1016/j.ultrasmedbio.2017.08.006 (2017).
Article PubMed PubMed Central Google Scholar
Hurme, T., Kalimo, H., Sandberg, M., Lehto, M. & Vuorio, E. Localization of type I and III collagen and fibronectin production in injured gastrocnemius muscle. Lab. Investig. 64, 76–84 (1991).
CAS PubMed Google Scholar
Wang, Y. et al. Collagen type III alpha1 as a useful diagnostic immunohistochemical marker for fibroepithelial lesions of the breast. Hum. Pathol. 57, 176–181. https://doi.org/10.1016/j.humpath.2016.07.017 (2016).
Article CAS PubMed Google Scholar
Kanneganti, T. D. Central roles of NLRs and inflammasomes in viral infection. Nat. Rev. Immunol. 10, 688–698. https://doi.org/10.1038/nri2851 (2010).
Article CAS PubMed PubMed Central Google Scholar
Barton, G. M. & Medzhitov, R. Toll-like receptors and their ligands. Curr. Top. Microbiol. Immunol. 270, 81–92 (2002).
CAS PubMed Google Scholar
He, Y., Hara, H. & Nunez, G. Mechanism and regulation of NLRP3 inflammasome activation. Trends Biochem. Sci. 41, 1012–1021. https://doi.org/10.1016/j.tibs.2016.09.002 (2016).
Article CAS PubMed PubMed Central Google Scholar
Dinarello, C. A. Overview of the IL-1 family in innate inflammation and acquired immunity. Immunol. Rev. 281, 8–27. https://doi.org/10.1111/imr.12621 (2018).
Article CAS PubMed PubMed Central Google Scholar
Kanehisa, M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 28, 1947–1951. https://doi.org/10.1002/pro.3715 (2019).
Article CAS PubMed PubMed Central Google Scholar
Kanehisa, M. & Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28, 27–30. https://doi.org/10.1093/nar/28.1.27 (2000).
Article CAS PubMed PubMed Central Google Scholar
Liu, M. et al. Metabolic rewiring of macrophages by CpG potentiates clearance of cancer cells and overcomes tumor-expressed CD47-mediated “don’t-eat-me” signal. Nat. Immunol. https://doi.org/10.1038/s41590-018-0292-y (2019).
Article PubMed PubMed Central Google Scholar
Nakae, S., Asano, M., Horai, R. & Iwakura, Y. Interleukin-1 beta, but not interleukin-1 alpha, is required for T-cell-dependent antibody production. Immunology 104, 402–409. https://doi.org/10.1046/j.1365-2567.2001.01337.x (2001).
Article CAS PubMed PubMed Central Google Scholar
Di Gioia, M. et al. Endogenous oxidized phospholipids reprogram cellular metabolism and boost hyperinflammation. Nat. Immunol. 21, 42–53. https://doi.org/10.1038/s41590-019-0539-2 (2020).
Article CAS PubMed Google Scholar
Zanoni, I. et al. An endogenous caspase-11 ligand elicits interleukin-1 release from living dendritic cells. Science 352, 1232–1236. https://doi.org/10.1126/science.aaf3036 (2016).
Article ADS CAS PubMed PubMed Central Google Scholar
Zanoni, I., Tan, Y., Di Gioia, M., Springstead, J. R. & Kagan, J. C. By capturing Inflammatory lipids released from dying cells, the receptor CD14 induces inflammasome-dependent phagocyte hyperactivation. Immunity 47, 697–709. https://doi.org/10.1016/j.immuni.2017.09.010 (2017).
Article CAS PubMed PubMed Central Google Scholar
Bochkov, V. N. et al. Protective role of phospholipid oxidation products in endotoxin-induced tissue damage. Nature 419, 77–81. https://doi.org/10.1038/nature01023 (2002).
Article ADS CAS PubMed Google Scholar
Bluml, S. et al. Oxidized phospholipids negatively regulate dendritic cell maturation induced by TLRs and CD40. J. Immunol. 175, 501–508. https://doi.org/10.4049/jimmunol.175.1.501 (2005).
Article PubMed Google Scholar
Munoz, L. E., Lauber, K., Schiller, M., Manfredi, A. A. & Herrmann, M. The role of defective clearance of apoptotic cells in systemic autoimmunity. Nat. Rev. Rheumatol. 6, 280–289. https://doi.org/10.1038/nrrheum.2010.46 (2010).
Article PubMed Google Scholar
Baccala, R., Hoebe, K., Kono, D. H., Beutler, B. & Theofilopoulos, A. N. TLR-dependent and TLR-independent pathways of type I interferon induction in systemic autoimmunity. Nat. Med. 13, 543–551. https://doi.org/10.1038/nm1590 (2007).
Article CAS PubMed Google Scholar
Andrade, F., Casciola-Rosen, L. & Rosen, A. Apoptosis in systemic lupus erythematosus. Clinical implications. Rheum. Dis. Clin. N. Am. 26, 215–227. https://doi.org/10.1016/s0889-857x(05)70136-8 (2000).
Article CAS Google Scholar
Janssen, E. et al. Efficient T cell activation via a Toll-Interleukin 1 receptor-independent pathway. Immunity 24, 787–799. https://doi.org/10.1016/j.immuni.2006.03.024 (2006).
Article CAS PubMed Google Scholar
Fritz, J. H., Ferrero, R. L., Philpott, D. J. & Girardin, S. E. Nod-like proteins in immunity, inflammation and disease. Nat. Immunol. 7, 1250–1257. https://doi.org/10.1038/ni1412 (2006).
Article CAS PubMed Google Scholar
Koga-Yamakawa, E. et al. TLR7 tolerance is independent of the type I IFN pathway and leads to loss of anti-tumor efficacy in mice. Cancer Immunol. Immunother. 64, 1229–1239. https://doi.org/10.1007/s00262-015-1730-4 (2015).
Article CAS PubMed Google Scholar
Kaplanov, I. et al. Blocking IL-1beta reverses the immunosuppression in mouse breast cancer and synergizes with anti-PD-1 for tumor abrogation. Proc. Natl. Acad. Sci. U.S.A. 116, 1361–1369. https://doi.org/10.1073/pnas.1812266115 (2019).
Article CAS PubMed Google Scholar
Voronov, E. & Apte, R. N. Targeting the tumor microenvironment by intervention in interleukin-1 biology. Curr. Pharm. Des. 23, 4893–4905. https://doi.org/10.2174/1381612823666170613080919 (2017).
Article CAS PubMed Google Scholar
Huber, V. et al. Tumor-derived microRNAs induce myeloid suppressor cells and predict immunotherapy resistance in melanoma. J. Clin. Investig. 128, 5505–5516. https://doi.org/10.1172/JCI98060 (2018).
Article PubMed PubMed Central Google Scholar
Sanchez, P. J., McWilliams, J. A., Haluszczak, C., Yagita, H. & Kedl, R. M. Combined TLR/CD40 stimulation mediates potent cellular immunity by regulating dendritic cell expression of CD70 in vivo. J. Immunol. 178, 1564–1572. https://doi.org/10.4049/jimmunol.178.3.1564 (2007).
Article CAS PubMed Google Scholar
Byrne, K. T. & Vonderheide, R. H. CD40 stimulation obviates innate sensors and drives T cell immunity in cancer. Cell Rep. 15, 2719–2732. https://doi.org/10.1016/j.celrep.2016.05.058 (2016).
Article CAS PubMed PubMed Central Google Scholar
Verma, V. et al. PD-1 blockade in subprimed CD8 cells induces dysfunctional PD-1(+)CD38(hi) cells and anti-PD-1 resistance. Nat. Immunol. https://doi.org/10.1038/s41590-019-0441-y (2019).
Article PubMed PubMed Central Google Scholar
Tsukamoto, H. et al. Combined blockade of IL6 and PD-1/PD-L1 signaling abrogates mutual regulation of their immunosuppressive effects in the tumor microenvironment. Cancer Res. 78, 5011–5022. https://doi.org/10.1158/0008-5472.CAN-18-0118 (2018).
Article CAS PubMed Google Scholar
Pahk, K. J. et al. Boiling histotripsy-induced partial mechanical ablation modulates tumour microenvironment by promoting immunogenic cell death of cancers. Sci. Rep. 9, 9050. https://doi.org/10.1038/s41598-019-45542-z (2019).
Article ADS CAS PubMed PubMed Central Google Scholar
Hendricks, A. D. et al. Histotripsy initiates local and systemic immunological response and reduces tumor burden in breast cancer. J. Immunol. 202, 130 (2019).
Google Scholar
Qu, S. et al. Non-thermal histotripsy tumor ablation promotes abscopal immune responses that enhance cancer immunotherapy. J. Immunother. Cancer https://doi.org/10.1136/jitc-2019-000200 (2020).
Article PubMed PubMed Central Google Scholar
Schade, G. R. et al. Boiling histotripsy ablation of renal cell carcinoma in the eker rat promotes a systemic inflammatory response. Ultrasound Med. Biol. 45, 137–147. https://doi.org/10.1016/j.ultrasmedbio.2018.09.006 (2019).
Article PubMed Google Scholar
Wong, A. W. et al. Ultrasound ablation enhances drug accumulation and survival in mammary carcinoma models. J. Clin. Investig. 126, 99–111. https://doi.org/10.1172/JCI83312 (2016).
Article PubMed Google Scholar
Cassetta, L. et al. Isolation of mouse and human tumor-associated macrophages. Adv. Exp Med Biol 899, 211–229. https://doi.org/10.1007/978-3-319-26666-4_12 (2016).
Article CAS PubMed PubMed Central Google Scholar
Borodina, T., Adjaye, J. & Sultan, M. A strand-specific library preparation protocol for RNA sequencing. Methods Enzymol. 500, 79–98. https://doi.org/10.1016/B978-0-12-385118-5.00005-0 (2011).
Article CAS PubMed Google Scholar
Levin, J. Z. et al. Comprehensive comparative analysis of strand-specific RNA sequencing methods. Nat Methods 7, 709–715. https://doi.org/10.1038/nmeth.1491 (2010).
Article CAS PubMed PubMed Central Google Scholar
Dobin, A. & Gingeras, T. R. Mapping RNA-seq reads with STAR. Curr. Protoc. Bioinform. 51, 11–19. https://doi.org/10.1002/0471250953.bi1114s51 (2015).
Article Google Scholar
Pertea, M. et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 33, 290–295. https://doi.org/10.1038/nbt.3122 (2015).
Article CAS PubMed PubMed Central Google Scholar
Trapnell, C. et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 7, 562–578. https://doi.org/10.1038/nprot.2012.016 (2012).
Article CAS PubMed PubMed Central Google Scholar
da Huang, W., Sherman, B. T. & Lempicki, R. A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4, 44–57. https://doi.org/10.1038/nprot.2008.211 (2009).
Article CAS Google Scholar
Newman, A. M. et al. Determining cell type abundance and expression from bulk tissues with digital cytometry. Nat. Biotechnol. 37, 773–782. https://doi.org/10.1038/s41587-019-0114-2 (2019).
Article CAS PubMed PubMed Central Google Scholar

Download references

Acknowledgements

The authors wish to thank Ms. Stephenie Y. Liu of the UC Davis Comprehensive Cancer Center (UCDCCC) Genomics Shared Resource (GSR) for her expert technical assistance. The UC Davis Comprehensive Cancer Center Genomics Shared Resource is supported by Cancer Center Support Grant P30CA093373 awarded by the NCI. The authors also wish to thank Meredith Weglarz, Lisa Nichols, and Ricardo Zermeno of the Stanford Shared FACS Facility for their guidance in panel construction and analysis. We acknowledge the support of R01CA210553, R01CA253316, and T32GM007276.

Author information

Authors and Affiliations

Department of Radiology, Stanford University, 3165 Porter Dr, Palo Alto, CA, 94305, USA
Brett Z. Fite, James Wang, Aris J. Kare, Asaf Ilovitsh, Tali Ilovitsh, Nisi Zhang, Weiyu Chen, Elise Robinson, Azadeh Kheirolomoom & Katherine W. Ferrara
Department of Biomedical Engineering, Stanford University, Palo Alto, CA, 94305, USA
Aris J. Kare & Michael Chavez
Department of Biomedical Engineering, University of California Davis, Davis, CA, 95616, USA
Hua Zhang, Matthew T. Silvestrini, Elizabeth S. Ingham, Lisa M. Mahakian & Sarah M. Tam
Department of Pathology and Laboratory Medicine, School of Medicine, University of California Davis, Sacramento, CA, 95817, USA
Ryan R. Davis & Alexander D. Borowsky
Department of Biochemistry and Molecular Medicine, School of Medicine, University of California Davis, Sacramento, CA, 95817, USA
Clifford G. Tepper

Authors

Brett Z. Fite
View author publications
You can also search for this author in PubMed Google Scholar
James Wang
View author publications
You can also search for this author in PubMed Google Scholar
Aris J. Kare
View author publications
You can also search for this author in PubMed Google Scholar
Asaf Ilovitsh
View author publications
You can also search for this author in PubMed Google Scholar
Michael Chavez
View author publications
You can also search for this author in PubMed Google Scholar
Tali Ilovitsh
View author publications
You can also search for this author in PubMed Google Scholar
Nisi Zhang
View author publications
You can also search for this author in PubMed Google Scholar
Weiyu Chen
View author publications
You can also search for this author in PubMed Google Scholar
Elise Robinson
View author publications
You can also search for this author in PubMed Google Scholar
Hua Zhang
View author publications
You can also search for this author in PubMed Google Scholar
Azadeh Kheirolomoom
View author publications
You can also search for this author in PubMed Google Scholar
Matthew T. Silvestrini
View author publications
You can also search for this author in PubMed Google Scholar
Elizabeth S. Ingham
View author publications
You can also search for this author in PubMed Google Scholar
Lisa M. Mahakian
View author publications
You can also search for this author in PubMed Google Scholar
Sarah M. Tam
View author publications
You can also search for this author in PubMed Google Scholar
Ryan R. Davis
View author publications
You can also search for this author in PubMed Google Scholar
Clifford G. Tepper
View author publications
You can also search for this author in PubMed Google Scholar
Alexander D. Borowsky
View author publications
You can also search for this author in PubMed Google Scholar
Katherine W. Ferrara
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

K.W.F. and B.Z.F. conceived the experiments, B.Z.F., A.I., T.I., and M.T.S carried out MRgFUS experimental studies. A.Kare., N.Z., W.C., E.R., H.Z., A.K., E.S.I., L.M.M, and S.M.T., performed flow cytometry. A.D.B. performed immnohistochemistry and analysis. R.R.D and C.G.T. carried out RNA sequencing and analysis. B.Z.F, K.W.F., J.W., M.C., and A.Kare. performed data analysis, B.Z.F. and K.W.F wrote the manuscript, all authors read the paper and participated in the discussions.

Corresponding author

Correspondence to Katherine W. Ferrara.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Publisher's note

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

Supplementary Information

Supplementary Information.

Rights and permissions

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

Reprints and permissions

About this article

Cite this article

Fite, B.Z., Wang, J., Kare, A.J. et al. Immune modulation resulting from MR-guided high intensity focused ultrasound in a model of murine breast cancer. Sci Rep 11, 927 (2021). https://doi.org/10.1038/s41598-020-80135-1

Download citation

Received: 15 July 2020
Accepted: 16 December 2020
Published: 13 January 2021
DOI: https://doi.org/10.1038/s41598-020-80135-1

This article is cited by

A glycoside hydrolase 12 protein from Cytospora chrysosperma triggers plant immunity but is not essential to virulence
- Zhiye Xu
- Zheng Luo
- Chengming Tian
Phytopathology Research (2023)
Thermal immuno-nanomedicine in cancer
- Zhe Yang
- Di Gao
- Xingcai Zhang
Nature Reviews Clinical Oncology (2023)
Values of OAS gene family in the expression signature, immune cell infiltration and prognosis of human bladder cancer
- Lijuan Gao
- Ruimin Ren
- Jimin Cao
BMC Cancer (2022)
Development of an ultrasound guided focused ultrasound system for 3D volumetric low energy nanodroplet-mediated histotripsy
- Bar Glickstein
- Ramona Aronovich
- Tali Ilovitsh
Scientific Reports (2022)
Development of a Systematic Review Protocol and a Scoping Review of Ultrasound-Induced Immune Effects in Peripheral Tumors
- Anne Rix
- Renée Girbig
- Fabian Kiessling
Molecular Imaging and Biology (2022)

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.