Hamed Mahdikia1,2, Fariba Saadati2, Eric Freund2,3, Udo S Gaipl4, Keivan Majidzadeh-A5, Babak Shokri1,6, Sander Bekeschus2. 1. Laser and Plasma Research Institute, Shahid Beheshti University, Tehran, Iran. 2. Center for Innovation Competence (ZIK) Plasmatis, Leibniz Institute for Plasma Science and Technology (INP), Greifswald, Germany. 3. Department of General, Visceral, Thoracic and Vascular Surgery, Greifswald University Medical Center, Greifswald, Germany. 4. Department of Radiation Oncology, Universitätsklinikum Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany. 5. Recombinant Proteins Department, Breast Cancer Research Center, Motamed Cancer Institute, ACECR, Tehran, Iran. 6. Department of Physics, Shahid Beheshti University,Tehran, Iran.
Abstract
While many new and emerging therapeutic concepts have appeared throughout the last decades, cancer still is fatal in many patients. At the same time, the importance of immunology in oncotherapy is increasingly recognized, not only since the advent of checkpoint therapy. Among the many types of tumors, also breast cancer has an immunological dimension that might be exploited best by increasing the immunogenicity of the tumors in the microenvironment. To this end, we tested a novel therapeutic concept, gas plasma irradiation, for its ability to promote the immunogenicity and increase the toxicity of breast cancer cells in vitro and in vivo. Mechanistically, this emerging medical technology is employing a plethora of reactive oxygen species being deposited on the target cells and tissues. Using 2D cultures and 3D tumor spheroids, we found gas plasma-irradiation to drive apoptosis and immunogenic cancer cell death (ICD) in vitro, as evidenced by an increased expression of calreticulin, heat-shock proteins 70 and 90, and MHC-I. In 4T1 breast cancer-bearing mice, the gas plasma irradiation markedly decreased tumor burden and increased survival. Interestingly, non-treated tumors injected in the opposite flank of mice exposed to our novel treatment also exhibited reduced growth, arguing for an abscopal effect. This was concomitant with an increase of apoptosis and tumor-infiltrating CD4+ and CD8+ T-cells as well as dendritic cells in the tissues. In summary, we found gas plasma-irradiated murine breast cancers to induce toxicity and augmented immunogenicity, leading to reduced tumor growth at a site remote to the treatment area.
While many new and emerging therapeutic concepts have appeared throughout the last decades, cancer still is fatal in many patients. At the same time, the importance of immunology in oncotherapy is increasingly recognized, not only since the advent of checkpoint therapy. Among the many types of tumors, also breast cancer has an immunological dimension that might be exploited best by increasing the immunogenicity of the tumors in the microenvironment. To this end, we tested a novel therapeutic concept, gas plasma irradiation, for its ability to promote the immunogenicity and increase the toxicity of breast cancer cells in vitro and in vivo. Mechanistically, this emerging medical technology is employing a plethora of reactive oxygen species being deposited on the target cells and tissues. Using 2D cultures and 3D tumor spheroids, we found gas plasma-irradiation to drive apoptosis and immunogenic cancer cell death (ICD) in vitro, as evidenced by an increased expression of calreticulin, heat-shock proteins 70 and 90, and MHC-I. In 4T1breast cancer-bearing mice, the gas plasma irradiation markedly decreased tumor burden and increased survival. Interestingly, non-treated tumors injected in the opposite flank of mice exposed to our novel treatment also exhibited reduced growth, arguing for an abscopal effect. This was concomitant with an increase of apoptosis and tumor-infiltrating CD4+ and CD8+ T-cells as well as dendritic cells in the tissues. In summary, we found gas plasma-irradiated murinebreast cancers to induce toxicity and augmented immunogenicity, leading to reduced tumor growth at a site remote to the treatment area.
Breast cancer is the most frequent type of cancer in women worldwide, with many of the cases being fatal.[1-3] The development of efficient and new treatment methods is difficult because of the complex nature of that type of cancer.[4] Most breast cancer-related deaths occur if metastatic cells begin to detach from the primary tumor and start to spread to surrounding tissue or secondary sites. During this process, cancer cells detach from the primary tumor by invading the basement membrane, crossing connective tissue, and penetrating other organs through blood vessels or lymphatic vessels.[5,6] In general, novel cancer therapies should be effective, safe, and associated with a minimum of side effects.[7,8]There is a variety of different traditional and novel treatment options for cancer, such as surgery, chemotherapy, radiotherapy, photodynamic therapy, and immunotherapy.[9-11] Especially the latter received increased attention during the recent years.[12] Clinical trials using immune checkpoint inhibitors in patients with triple-negative breast cancer, such as atezolizumab targeting PD-L1, showed promising results.[13,14] This builds the case for the power of unleashing antitumor immunity in breast cancerpatients, which was observed already for other tumor entities.[15,16] Tumor cells mutate at high frequencies, leading to an elevated mutational burden that correlates to responses to immunotherapy in patients.[17] For breast cancer, on average, somatic mutations are regularly observed,[18] facilitating antitumor T-cells responses for attacking metastatic lesions.However, for antitumor immunity to form, tumor antigens not only need to be present but also accessible and processable for antigen-presenting cells to promote T-cell mediated immunity. Radiotherapy has been proven to exert immune modulatory effects.[19] Locally applied ionizing radiation can thereby potentiate antitumor immunity, leading to reduced cancer growth also of distant, untreated sites (abscopal effect).[20] It is noteworthy that radiotherapy, to some extent, relies on the local generation of reactive oxygen species (ROS).[21] Further, radiotherapy synergizes with other ROS-inducing treatments such as gas plasma in the induction of cancer cell death with immune stimulatory properties.[22] Therapeutic levels of ROS are increasingly appreciated as novel paradigm for cancer treatment.[23]Gas plasma technology is an emerging innovation in oncology.[24] It is recognized that the main mode of action of these partially ionized gases operated at body temperature relies on the generation of ROS.[25] The unique feature of gas plasma systems is that a plethora of different types of ROS is being produced simultaneously,[26] facilitating oxidative damage of tumor cells and, ultimately, leading to their demise. While the concept of ROS-mediated killing using gas plasmas has been explored to some extent,[27] the immunological consequences of gas plasma irradiation of breast cancer cells are unknown. To this end, we investigated gas plasma-irradiated breast cancer cells in terms of toxicity and immunogenicity in vitro and in vivo. Gas plasma irradiation not only effectively controlled tumor growth but also was concomitant with an upregulation of several markers of the immunogenic cancer cell death (ICD), such as calreticulin (CRT) and the damage-associated molecular pattern (DAMP) molecules heat-shock protein (HSP) 70 and HSP90.[28,29] Finally, we found gas plasma irradiation to mediate an abscopal effect in breast-cancer bearing mice, mounting evidence of a promising role of therapeutic ROS derived from gas plasma systems in promoting anticancer immunity.
Materials and methods
Cell culture and in vitro gas plasma irradiation
Human breast cancer cell lines MCF-7 (ATCC: HTB-22) and MDA-MB (ATCC: HTB-26) were cultivated in Dulbecco’s Modified Eagles’ Medium (DMEM; Corning). The murinebreast cancer cell line 4T1 (ATCC: CRL-2539) was cultured in Roswell Park Memorial Institute (RPMI) 1640 medium (Corning). All culture media were supplemented with 10% fetal bovine serum, 2% L-glutamine, and 1% penicillin (all Sigma Aldrich), and were kept under standard culture conditions at 37°C, 5% CO2, and 95% humidity in a cell culture incubator (Binder). Twenty-four hours before experiments, the cells were seeded at 5 × 104 cells per ml of fully supplemented medium in 24-well tc-treated cell culture plates (Eppendorf). Three-dimensional tumor spheroids were formed as described before.[30] Briefly, cells were seeded at 4 × 103 cells per 200 µl fully supplemented medium in 96-well ultra-low attachment plates (Corning). After centrifugation at 1000 x g for 10 min, the cells were cultured for 72 h before experiments. For the gas plasma irradiation of the cells in vitro, a helium-driven plasma jet was utilized. Its technical parameters and physical properties have been described in detail before.[26] It was operated with a gas flow four standard liters helium (purity: 99.999%; Air Liquide). The distance of the target to plasma as well as the treatment time was kept constant using a software-controlled xyz-table (CNC Step). Vehicle controls received helium gas only, through the same system.
Reactive species detection
For the detection of reactive species introduced into the cell culture medium following gas plasma irradiation, several different assays were applied. Hydrogen peroxide (H2O2) was quantified utilizing the Amplex Ultra Red assay kit (ThermoFisher) according to the manufacturers’ instructions. For the detection of nitrite (NO2−) inside the medium, a deterioration product of nitric oxide (NO), the Griess assay (Cayman Chemical) was carried out. Hypochlorous acid (HOCl) was detected by the colorimetric taurinechloramine assay, as previously described.[31] The measurements were performed using a multiplate reader (Tecan Infinite F200 PRO) at λex/em = 535 nm/590 nm for fluorescence detection (H2O2), or the (Tecan Infinite M200 PRO) at λ = 540 nm (NO2−) and at λ = 645 nm (HOCl) for absorbance detection. The pH levels of gas plasma-irradiated cell culture medium were quantified using a pH-meter (Mettler Toledo).
Metabolic activity and cytotoxicity
Resazurin (Alfa Aesar), a non-fluorescent dye becoming fluorescent after reduction via cell-derived NADPH, was added at 22 h or 46 h after gas plasma irradiation at a final concentration of 100 µM and incubated for 2 h. Fluorescence was detected by a multimode plate reader (Tecan Infinite F200 PRO) at λex 535 nm and λem 590 nm. For assessing cytotoxicity via microscopy, 4,6-Diamidin-2-phenylindol (DAPI, final concentration 1 µM; BioLegend) was added to each well. Subsequently, the cells were imaged using a high-content imaging analysis system (Operetta CLS; PerkinElmer) with a 20x (NA = 0.4) air objective (Zeiss) equipped with laser-based autofocus (785 nm). Imaging was performed using the brightfield channel (with a 740 nm transmission light source), the digital phase contrast (DPC) channel (generating a single image from differential contrasts that reflect the cytosolic region of a cell), and the DAPI channel (λex 365 nm and λem 465 ± 35 nm). The data were analyzed using quantitative imaging and segmentation software (Harmony 4.9; PerkinElmer). For analysis of tumor spheroids, they were stained with either MitoTracker Red (MTR) CMXRos (ThermoFisher), sytox green (ThermoFisher), or anti-calreticulin (CRT) Alexa Fluor (AF) 647 (Novus) to assess mitochondrial activity, cell death, and immunogenic cell death, respectively. Imaging was done by using the high content imaging device with a 5x air (NA = 0.16) objective and four z-stacks per spheroid in non-confocal mode. These z-stacks of each spheroid were merged using maximum intensity projections. The spheroids area was segmented and evaluated using Harmony software for the extent of cell death (sytox green intensity per spheroid), the oxidation (MTR intensity per spheroid), and CRT expression (AF647 intensity per spheroid).
Flow cytometry
To investigate the extent of apoptosis following gas plasma irradiation, the cells were incubated with caspase 3/7 detection reagent (ThermoFisher) and DAPI. Cells were detached using accutase (BioLegend) and analyzed using flow cytometry (CytoFLEX S and LX; Beckman-Coulter). For analysis of immunogenic cell death (ICD) markers, the cells were incubated with monoclonal antibodies targeted against CRT and labeled with AF647 (Novus), heath-shock protein (HSP) 70 AF488 (BioLegend), HSP90 PE (Enzo), MHC-I PE-Cy7 (BioLegend), CD274 BV650 (BioLegend), as well as with iFluor 860 (AAT Bioquest). For the detection of intracellular antigens, the cells were fixed and permeabilized (fixation and permeabilization wash buffer; BioLegend) and stained with antibodies targeted against LC3 AF488 (Cell Signaling), or phoshpho-eIF2α labeled with AF488+ (Invitrogen; conjugate ThermoFisher), phospho-ATM labeled with PE (BioLegend), and ERp75 labeled with AF647+ (Santa Cruz; conjugate ThermoFisher). The cells were acquired using flow cytometry, and data analysis was performed using Kaluza 2.1.1 software (Beckman-Coulter).
Supernatant analysis
Cellular supernatants were collected 1 h and 48 h after in vitro treatment procedures. Residual cells were removed via centrifugation. ATP was detected utilizing the Luminescent ATP Detection Kit (Abcam), and luminescence was measured using a multimode plate reader (Tecan). Enzyme-linked immunosorbent assays were utilized according to the manufacturers’ procedures for the quantification of HSP70 (ThermoFisher), IFN2α, IFNγ, and IL-6 (all BioLegend). The absorbance was measured using the same multiplate reader (Tecan).
In vivo tumor model
The in vivo experiments were approved by the Shahid Beheshti University of Medical Science Animal Ethics Committee (approval number IR.SBMU.RETECH.REC.1397.1212) and were conducted at the Shahid Beheshti University. The animals were housed under standard laboratory conditions (room temperature = 23°C), with an ad libitum supply of nutrition and water, and circadian environmental control (12 h light shifts). All interventions at the animals were performed under anesthesia through the injection of ketamine and xylazine. Subcutaneous breast cancers were induced by subcutaneous injection of 5 × 105 4T1 cells (ATCC CRL-2539) in each flank of 6–7 weeks old female Balb/CfC3H mice (Royan Institute). Daily gas plasma irradiations were started after the tumor diameters reached 5 mm at day 12 and were performed for 21 days in total. A helium gas-driven plasma multi-jet was utilized that was operated at 20 kHz frequency with a total input power of 1 W. The gas plasma irradiation time was 300 s per tumor and was only applied to the left tumors in the irradiation group. In order to investigate tumor growth and immune cell infiltration, a portion of the animals was sacrificed at day 36 for further downstream analysis (untreated tumor: n = 16; irradiation left and right: n = 8 each). The survival of untreated and gas plasma-irradiated animals was investigated until day 72 (untreated: n = 10; irradiation: n = 10).
Tumor growth monitoring in vivo
The tumors’ greatest longitudinal and transverse diameter were measured every three days. The tumor volume was calculated using a modified ellipsoidal formula ()). The tumor weight at day 36 was measured after tumor explantation via a precision balance. The tumor density was calculated to be 1.056 g/ml (). For all tumor growth evaluations, the calculated tumor volume was normalized to that of the corresponding initial tumor volume at day 12. Moreover, sonography (ultrasound) was performed to measure the greatest tumor length, width, and height, as well as to identify hypo or hyper-dense structures within the tissue. Representative animals were additionally analyzed via computer tomography (CT) at the Faculty of Veterinary Medicine of the University of Tehran. The calculation of three-dimensional animal reconstructions, surface rendering, and tumor volume quantification was performed with the Horos 3.3.6 medical image viewer (Horos Project).
Tissue sectioning and staining
Tissue sections of 5 µm thickness were cut and immersed in xylene and different percentages of ethanol to remove the paraffin and perform rehydration. For antigen retrieval, tissue slides were placed in citrate buffer (10 mM, PH 6.0) and were washed three times with PBS containing 1% Tween (PBS-T). Peroxidase blockade was done by adding a H2O2 and methanol mixture in PBS-TTween for 30 min at room temperature, followed by a similar procedure using 3% Bovine Serum Albumin (BSA) in PBS-T. After the blocking buffer was removed, 70 µl of primary antibodies solution was added. These contained anti-mouse antibodies targeted against active caspase 3, CRT, CD4, CD8, or CD11 (all cell signaling technology). After overnight incubated at 4°C, and followed by incubation with an AF555 secondary antibody for 1 h at room temperature, Antifade Mounting Medium (Vectashield) with DAPI (LSBio) was used for nuclei staining and mounting. Slide imaging was done using the high content imaging system (Operetta CLS; PerkinElmer) in 10 z- stuck with 5 µm thickness with 20x air objective to capture DAPI and AF555 emission (λex 505 nm and λem 553 ± 23 nm). Quantitative image analysis was performed using Harmony 4.9 Software (PerkinElmer).Terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) technology labels DNA fragmentation in the last phase of apoptosis (In Situ Cell Death Detection Kit, Fluorescein, Sigma-Aldrich). Paraffin-embedded tumor tissue sections were deparaffinized and rehydrated by immersing in xylene and different percentages of ethanol. Then, the tissue sections were pretreated with proteinase k working solution (15 µg/ml in 10 mM Tris/HCl) as a permeabilization solution for 15 min at 37°C. Slides were rinsed three times with PBS and pretreated with 50 µl of TUNEL reaction mixture (50 µl enzyme solution mixed with 450 µl label solution) in a 37°C humidified atmosphere in the dark for 60 min. The label solution also was used as a negative control for accurate target detection. In the following, the slides were incubated with 50 µl DAPI (20 µM) for 15 min and were covered by a drop of mounting medium and cover slide. TUNEL-stained and negative control tissues were imaged and quantified, as described above. For the detection of intracellular FOXP3 and IL-17 in CD4+ T-cells, the tissue slides were processed as described above and later fixed and permeabilized with TBS-T, as well as Triton X-100 and FOXP3 staining buffer (BioLegend). Afterward, the nuclei were counterstained with DAPI, and the slides were incubated with antibodies targeting CD4, FOXP3, IL-17A, and IL17F (all BioLegend). The outer tissue border region was quantified by the Harmony software by counting double-positive cells.
Results
Gas plasma irradiation generated ROS that impaired breast cancer cell growth
A helium (He) gas plasma jet was used for the treatment of human breast cancer cells (MCF-7 and MDA-MB) in the in vitro experiments (Figure 1a) that expels reactive oxygen and nitrogen species. In gas plasma-irradiated liquids, these species deteriorate then further to oxidants that can be quantified more conveniently with redox chemical assays. Accordingly, a gas plasma irradiation time-dependent increase of hydrogen peroxide (H2O2) as a product of OH (Figure 1b), nitrite (NO2−) as a product of nitric oxide (Figure 1c), and hypochlorous acid (HOCl) as a product of O (Figure 1d) was observed in the gas plasma-irradiated cell culture medium. Concomitantly, a modest increase in the pH of the gas plasma-irradiated liquid was observed (Figure 1e). Next, gas plasma irradiation was applied to breast cancer cells in vitro (Figure 1f). The treatment led to a significant increase in terminally dead MCF-7 cells (Figure 1g), and a subsequent significant decrease in the metabolic activity of the cells was observed at 24 h and 48 h post-treatment measured in the entire well containing the cells (Figure 1h). In MDA-MB breast cancer cells, gas plasma-induced terminal cell death was more pronounced as compared to MCF-7 cells (Figure 1i), while the extent of reduction of metabolic activity at 24 h and 48 h post-irradiation measured in the entire well containing the cells was comparable (Figure 1j). Altogether, gas plasma irradiation generated reactive species in the gas phase and treated liquids, which subsequently reduced breast cancer cell viability and metabolic activity significantly.
Figure 1.
Gas plasma irradiation-generated redox chemistry and metabolic activity in breast cancer cells. (a) scheme of the helium (He) gas plasma irradiation of cells in vitro; (b-e) quantification of hydrogen peroxide (H2O2, b), nitrite (NO2−, c), hypochlorous acid (HOCl, d), and pH (e) in gas plasma-irradiated cell culture medium; (f) representative brightfield and DAPI (terminally dead cells, blue) images of MCF-7 and MDA-MB breast cancer cells; (g) quantitative image analysis of dead cells in MCF-7 cultures; (h) metabolic activity in MCF-7 cells at 24 h and 48 h post gas plasma irradiation measured in the entire well containing the cells; (i) quantitative image analysis of dead cells in MDA-MB cultures; (j) metabolic activity in MDA-MB cells at 24 h and 48 h post gas plasma irradiation measured in the entire well containing the cells. Cell data are from three independent experiments. Data are presented as mean (+SD). Statistical analysis was performed using one-way analysis of variances with p < .05 (*) and p < .001 (***); ns = non-significant; vehicle = helium gas treatment alone (plasma ignition off); scale bar is 100 µm
Gas plasma irradiation-generated redox chemistry and metabolic activity in breast cancer cells. (a) scheme of the helium (He) gas plasma irradiation of cells in vitro; (b-e) quantification of hydrogen peroxide (H2O2, b), nitrite (NO2−, c), hypochlorous acid (HOCl, d), and pH (e) in gas plasma-irradiated cell culture medium; (f) representative brightfield and DAPI (terminally dead cells, blue) images of MCF-7 and MDA-MB breast cancer cells; (g) quantitative image analysis of dead cells in MCF-7 cultures; (h) metabolic activity in MCF-7 cells at 24 h and 48 h post gas plasma irradiation measured in the entire well containing the cells; (i) quantitative image analysis of dead cells in MDA-MB cultures; (j) metabolic activity in MDA-MB cells at 24 h and 48 h post gas plasma irradiation measured in the entire well containing the cells. Cell data are from three independent experiments. Data are presented as mean (+SD). Statistical analysis was performed using one-way analysis of variances with p < .05 (*) and p < .001 (***); ns = non-significant; vehicle = helium gas treatment alone (plasma ignition off); scale bar is 100 µm
Gas plasma irradiation promoted toxicity and immunogenicity in 2D and 3D cell models
To identify the mode and immunogenic consequences of gas plasma irradiation in breast cancer cells, flow cytometry was subsequently employed. In MCF-7 cells (Figure 2a), gas plasma irradiation led to a significant decline in the percentage of cells negative for active caspases 3 and 7 (Figure 2b) being a marker of apoptosis. In MDA-MB cells (Figure 2c), a much more pronounced induction of apoptosis was observed (Figure 2d) as compared to the MCF-7 cells. The latter were subsequently analyzed for the expression of CRT, a molecule known to be increasingly expressed in cells succumbing to immunogenic cancer cell death (ICD), heat-shock-proteins (HSP) 70 and HSP90, major histocompatibility complex (MHC)-1, as well as death-ligand 1 (PD-L1) (Figure 2e). All molecules are reported to be associated with ICD and anticancer immunity,[32-36] and a significant elevation was observed in gas plasma-irradiated but viable (caspase 3/7−) cells across all markers investigated (Figure 2f-j). When analyzing these markers in MDA-MB cells (Figure 2k), a significant increase of CRT was observed (Figure 2l). Also for HSP70 (Figure 2m), HSP90 (Figure 2n), MHC-I (Figure 2o), and PD-L1 (Figure 2p), significantly elevated levels were observed in gas plasma-irradiated tumor cells. In increase in the autophagy-related marker LC3, however, was observed only in MCF-7 but not MDA-MB cells (Figure 2q). As additional markers for cell death signaling, also the endoplasmic reticulum protein 57 (ERp57), the phosphorylated ATM kinase (pATM), as well as the phosphorylated eukaryotic initiation factor 2a (peIF2a) were significantly upregulated in MCF-7, but not MDA-MB cells post-treatment Supplementary (Figure 1a-c). Besides the surface-upregulation of the ICD marker, the secretion of several inflammation and DAMP-factors cytokines was investigated in the cells’ supernatants (Figure 2r-s). The levels differed when analyzed at early (1 h) versus late (48 h) collection time points. Both cell lines showed a significant increase of HSP70 as well as interferon (IFN) α2 in the cell culture supernatant (Figure 2r-s). ATP was increased early (1 h) but not late (48 h) after gas plasma irradiation. IFNγ was only found to be increased after the treatment of MDA-MB cells but not MCF-7(Figure 2s). In contrast, interleukin (IL)-6, known as a general pro-inflammatory signal molecule, was elevated 48 h after exposure of MCF-7 cells to gas plasma irradiation (Figure 2r).
Figure 2.
Cell viability and markers of immunogenic cancer cell death in breast cancer cells. (a) representative flow cytometry forward scatter (FSC) and caspase 3/7 dot-plots of MCF-7 cells; (b) quantification of viability data from flow cytometry of MCF-7 cells; (c) representative flow cytometry forward scatter (FSC) and caspase 3/7 dot-plots of MDA-MB cells; (d) quantification of viability data from flow cytometry of MDA-MB cells; (e) overlay histogram from flow cytometry of calreticulin (CRT), heat-shock protein (HSP) 70, HSP 90, major histocompatibility complex (MHC)-I, and programmed death-ligand 1 (PD-L1) in MCF-7 cells; (f-j) quantification of CRT (f), HSP70 (g), HSP90 (h), MHC-I (i), and PD-L1 (j) in viable (caspase 3/7−) MCF-7 cells; (k) overlay histogram from flow cytometry of CRT, HSP70, HSP90, MHC-I, and PD-L1 in MDA-MB cells; (l-p) quantification of CRT (l), HSP70 (m), HSP90 (n), MHC-I (o), and PD-L1 (p) in viable (caspase 3/7−) MDA-MB cells; (q) quantification of LC3 expression 6 h post treatment in MCF-7 and MDA-MB cells; (r-s) early (1 h) and late (48 h) cytokine expression of adenosine triphosphate (ATP), HSP70, interferon (IFN) α2, IFNγ, and interleukin (IL)-6 in MCF-7 (r), and MDA-MB cells (s).Data are representative of at least three independent experiments and presented as mean (min-max boxplot; b,d), individual values (e-j, l-p), mean + SD (q), or mean +SEM (r-s). Statistical analysis was performed using one-way analysis of variances with p < .05 (*) and p < .001 (***)
Cell viability and markers of immunogenic cancer cell death in breast cancer cells. (a) representative flow cytometry forward scatter (FSC) and caspase 3/7 dot-plots of MCF-7 cells; (b) quantification of viability data from flow cytometry of MCF-7 cells; (c) representative flow cytometry forward scatter (FSC) and caspase 3/7 dot-plots of MDA-MB cells; (d) quantification of viability data from flow cytometry of MDA-MB cells; (e) overlay histogram from flow cytometry of calreticulin (CRT), heat-shock protein (HSP) 70, HSP 90, major histocompatibility complex (MHC)-I, and programmed death-ligand 1 (PD-L1) in MCF-7 cells; (f-j) quantification of CRT (f), HSP70 (g), HSP90 (h), MHC-I (i), and PD-L1 (j) in viable (caspase 3/7−) MCF-7 cells; (k) overlay histogram from flow cytometry of CRT, HSP70, HSP90, MHC-I, and PD-L1 in MDA-MB cells; (l-p) quantification of CRT (l), HSP70 (m), HSP90 (n), MHC-I (o), and PD-L1 (p) in viable (caspase 3/7−) MDA-MB cells; (q) quantification of LC3 expression 6 h post treatment in MCF-7 and MDA-MB cells; (r-s) early (1 h) and late (48 h) cytokine expression of adenosine triphosphate (ATP), HSP70, interferon (IFN) α2, IFNγ, and interleukin (IL)-6 in MCF-7 (r), and MDA-MB cells (s).Data are representative of at least three independent experiments and presented as mean (min-max boxplot; b,d), individual values (e-j, l-p), mean + SD (q), or mean +SEM (r-s). Statistical analysis was performed using one-way analysis of variances with p < .05 (*) and p < .001 (***)To confirm terminal cell death and ICD in a more realistic in vitro model, 3D tumor spheroids were generated and exposed to gas plasma irradiation (Figure 3a). In MCF-7 cells, the effect on terminal cell death (Figure 3b) and cellular oxidation (Figure 3c) was modest, while CRT expression was significantly elevated with the treatment (Figure 3d). For MDA-MB cells, a significantly increased amount of terminally dead cells was observed (Figure 3e), which was concomitant with a significant increase in intracellular oxidation (Figure 3f) and CRT (Figure 3g). The differences between both cell lines in terms of cytotoxic effects might be related to the MCF-7 cells having a more epithelial-like phenotype at the spheroids’ vicinity, leading to denser 3D tumors with confined boarders. The data analysis was based on the fluorescence signals from the segmented spheroids, not taking into account the viability of individual cells due to the technical limitation of unambiguously identifying all individual cells within the tumor spheroids. In summary, gas plasma irradiation induced apoptosis in human breast cancer cell lines in 2D monolayer and 3D tumor spheroid cultures, and led to a significant elevation of the molecules CRT, HSP70, HSP90, MHC-I, as well as the secretion of ATP and interferons known to be involved in ICD.
Figure 3.
Viability, oxidation, and CRT expression in 3D breast cancer spheroids. (a) representative maximum intensity projection images of MCF-7 and MDA-MB 3D multicellular breast cancer spheroids generated from an overlay digital phase contrast and a sytox green image; (b-d) quantitative image analysis of the dead cell marker sytox green (b), oxidation marker mitotracker red (c), and CRT (d) of the segmented spheroid image area of MCF-7 cells at 24 h and 48 h post gas plasma irradiation exposure; (e-g) quantitative image analysis of sytox green (e), mitochondrial oxidation (f), and CRT (g) of the segmented spheroid image area of MDA-MB cells at 24 h and 48 h post gas plasma irradiation exposure. The data analysis was based on the fluorescence signals from the segmented spheroids, not taking into account the individual viability of the cells (except for sytox green). Data are from three to eight spheroids and presented as mean (min-max boxplot; ±SD). Statistical analysis was performed using one-way analysis of variances with p < .001 (***); scale bar is 300 µm
Viability, oxidation, and CRT expression in 3D breast cancer spheroids. (a) representative maximum intensity projection images of MCF-7 and MDA-MB 3D multicellular breast cancer spheroids generated from an overlay digital phase contrast and a sytox green image; (b-d) quantitative image analysis of the dead cell marker sytox green (b), oxidation marker mitotracker red (c), and CRT (d) of the segmented spheroid image area of MCF-7 cells at 24 h and 48 h post gas plasma irradiation exposure; (e-g) quantitative image analysis of sytox green (e), mitochondrial oxidation (f), and CRT (g) of the segmented spheroid image area of MDA-MB cells at 24 h and 48 h post gas plasma irradiation exposure. The data analysis was based on the fluorescence signals from the segmented spheroids, not taking into account the individual viability of the cells (except for sytox green). Data are from three to eight spheroids and presented as mean (min-max boxplot; ±SD). Statistical analysis was performed using one-way analysis of variances with p < .001 (***); scale bar is 300 µm
Gas plasma irradiation reduced tumor growth in vivo at the treated and remote side
The next question was whether gas plasma irradiation conferred cytotoxic and pro-immunogenic effects not only in vitro but also in vivo. For this, 4T1 syngeneic tumor cells were injected in both flanks of Balb/c mice (Figure 4a). In the control group, both tumors of each mouse remained untreated. In the gas plasma irradiation group, only the left tumor in each mouse was treated, while the right tumor of that same animal remained untreated. This way, we aimed at investigating remote, immune-mediated effects of the treatment on distant tumor sites following exposure to the gas plasma irradiation scheme (Figure 4b). The treatment led to a significantly increased survival in tumor-bearing mice (Figure 4c), which was overall prolonged compared to the control group (Figure 4c-d). This was also appreciated when comparing the tumor growth of the individual animals (Figure 4e) and as well as the groups (Figure 4f). Remarkably, the gas plasma irradiation not only decreased the tumor growth and weight of the directly treated breast tumor (left flank) but simultaneously, also the remote breast tumor at the opposite, non-treated side (right flank) showed a significantly impaired growth as compared to tumors of the untreated animals (Figure 4g), suggesting distant, immune-mediated effects. The principal differences between left and right tumors in the gas plasma-irradiated group, and tumors of the control group, were confirmed using Doppler sonography (Figure 4h) and computer tomography (CT; Figure 4i). All in all, these results implicated irradiation via gas plasma technology not only to decelerate tumor growth, leading to increased survival of breast cancer-bearing mice but also to stimulate possibly immune-mediated effects resulting in impaired tumor growth at distant, untreated sites.
Figure 4.
Syngeneic breast cancer animal model and tumor growth. (a-b) scheme of the animal treatment sides (a) and the timeline of the treatment of the breast cancer animal model (b); (c-d) Kaplan-Meier survival curve (c) and area under the curve (d) of the curves; (e) tumor sizes of individual mice of the control group (black) and the gas plasma-irradiated group carrying one tumor on each flank, one being subjected to gas plasma irradiation treatment (magenta) and the tumor on the opposite side remaining untreated (turquoise); (f-g) mean tumor growth (f) and weight (g) of the untreated tumors, gas plasma irradiation-treated tumors as well as untreated opposite flank tumors of the same mice; (h-i) representative images of ultrasound (sonography) measurements of the tumors (h), and of a representative mouse of each group using computer tomography (CT). Data are from eight mice per group. Statistical analysis was performed using one-way analysis of variances with p < .05 (*) and p < .01 (**)
Syngeneic breast cancer animal model and tumor growth. (a-b) scheme of the animal treatment sides (a) and the timeline of the treatment of the breast cancer animal model (b); (c-d) Kaplan-Meier survival curve (c) and area under the curve (d) of the curves; (e) tumor sizes of individual mice of the control group (black) and the gas plasma-irradiated group carrying one tumor on each flank, one being subjected to gas plasma irradiation treatment (magenta) and the tumor on the opposite side remaining untreated (turquoise); (f-g) mean tumor growth (f) and weight (g) of the untreated tumors, gas plasma irradiation-treated tumors as well as untreated opposite flank tumors of the same mice; (h-i) representative images of ultrasound (sonography) measurements of the tumors (h), and of a representative mouse of each group using computer tomography (CT). Data are from eight mice per group. Statistical analysis was performed using one-way analysis of variances with p < .05 (*) and p < .01 (**)
Gas plasma irradiation augmented immune infiltration in the tumor microenvironment
To next understand the composition of the tumor microenvironment that might be related to the antitumor effects observed with gas plasma irradiation, tissue sections were done followed by immunofluorescence staining. The quantification of positively-stained cells was done using algorithm-driven image segmentation. The TUNEL staining is indicative of apoptosis (Figure 5a), and a significant increase of TUNEL+ cells was observed not only in gas plasma-irradiated tumors but also in opposite-flank tumors (Figure 5b). Along similar lines, the number of cells staining positive for active caspase 3, an executioner caspase of the apoptotic pathway, was significantly increased in both left and right flank tumors of gas plasma-irradiated mice when compared to untreated control tumors (Figure 5c). Subsequently, the immune infiltrate in tumor tissues was quantified (Figure 5d). For CD8+ cytotoxic T-cells, a significant increase was observed in the gas plasma-irradiated tumors (Figure 5e). The change in the untreated tumors of irradiated mice was modest but not significant compared to the control tumors. By contrast, the proportion of CD4+ T-helper cells was significantly elevated in both the gas plasma-irradiated and the opposite-flank tumor within the same animal when compared to the numbers identified in control mice (Figure 5f). A similar trend was observed for CD11+ cells, indicative of dendritic cells, which were significantly elevated in the gas plasma-treated mice (Figure 5g). Investigating CRT as an indicator ICD, a significantly pronounced staining was observed in both the gas plasma-treated tumors as well as the distant non-treated tumors within the same animal, when compared to those levels found in the tumors of untreated animals. Finally, when staining for subpopulations of CD4+ cells and quantifying their numbers in the tissues, a significantly lower amount of FOXP3+ CD4+ cells was found at the irradiation site of the animals (Figure 5i-j). IL-17 expression showed a pronounced upregulation in CD4+ T-cells of the left irradiated tumor, in contrast to the opposite sided tumor and the untreated group (Figure 5k). The analysis of CD4, CD8, CD11c, CRT, FOXP3, and IL-17 was done for all cells staining positive for these markers, not taking into account their individual viability as technical limitations did not allow for such degree of multiplexing. In sum, analysis of the tumor microenvironment revealed a significant elevation of apoptosis and cell death, which was concomitant with the promotion of immunogenicity and immune cell infiltration following irradiation with gas plasma technology.
Figure 5.
Quantitative imaging analysis of apoptosis, immune cell infiltrates, and CRT in breast cancer tissues. (a) representative images of tumor tissue sections stained with DAPI (nuclei, blue) for all cells and TUNEL (green) for apoptotic cells; (b-c) quantitative image analysis of TUNEL+ (b) and caspase 3+ (c) cells in tissue sections; (d) representative whole tissue section image that was stained (left) with DAPI and (blue) and anti-CD8 antibodies (orange), and segmented for different image regions (right) based on algorithm-driven image processing for unsupervised detection of cells; (e-h) quantitative image analysis results obtained for objects staining positive for CD8 (e), CD4 (f), CD11 c (g), and CRT (h); (i) software detection of tissue border region with the quantification of (j) FOXP3+, and (k) IL-17+ subpopulations of CD4+ cells. The latter analysis was done for all cells staining positive for these markers, not taking into account their viability as technical limitations did not allow for such degree of multiplexing. Data show individual values and mean or min-max boxplots, or as bar graphs with mean +SEM of at least seven tumor sections per group. Statistical analysis was performed using one-way analysis of variances, or t-tests with p < .05 (*), p < .01 (**), and p < .001 (***); scale bars are 100 µm, 800 µm, and 50 µm, respectively
Quantitative imaging analysis of apoptosis, immune cell infiltrates, and CRT in breast cancer tissues. (a) representative images of tumor tissue sections stained with DAPI (nuclei, blue) for all cells and TUNEL (green) for apoptotic cells; (b-c) quantitative image analysis of TUNEL+ (b) and caspase 3+ (c) cells in tissue sections; (d) representative whole tissue section image that was stained (left) with DAPI and (blue) and anti-CD8 antibodies (orange), and segmented for different image regions (right) based on algorithm-driven image processing for unsupervised detection of cells; (e-h) quantitative image analysis results obtained for objects staining positive for CD8 (e), CD4 (f), CD11 c (g), and CRT (h); (i) software detection of tissue border region with the quantification of (j) FOXP3+, and (k) IL-17+ subpopulations of CD4+ cells. The latter analysis was done for all cells staining positive for these markers, not taking into account their viability as technical limitations did not allow for such degree of multiplexing. Data show individual values and mean or min-max boxplots, or as bar graphs with mean +SEM of at least seven tumor sections per group. Statistical analysis was performed using one-way analysis of variances, or t-tests with p < .05 (*), p < .01 (**), and p < .001 (***); scale bars are 100 µm, 800 µm, and 50 µm, respectively
Discussion
Despite recent advances in therapy, breast cancer still is fatal in many patients. In thriving for finding novel treatment avenues, we employed an innovative medical gas plasma irradiation approach for tackling breast cancer. Treatment with this technology not only limited tumor growth in vitro and in vivo but also had an immunological dimension. Besides an upregulation of markers of the immunogenic cell death (ICD), we found in vivo untreated tumors showing decelerated growth if present in mice having tumors that had received the treatment, arguing for an abscopal effect (Figure 6).
Figure 6.
Scheme of gas plasma irradiation in breast cancer treatment. Gas plasma irradiation generates multiple types of ROS and RNS simultaneously, which subsequently induce cytotoxicity in breast cancer cells in vitro and in vivo. Concurrently, cell death occurs in an immunogenic fashion, leading to immune-mediated tumor reduction at a distant site not exposed to gas plasma irradiation (abscopal effect) and improved survival
Scheme of gas plasma irradiation in breast cancer treatment. Gas plasma irradiation generates multiple types of ROS and RNS simultaneously, which subsequently induce cytotoxicity in breast cancer cells in vitro and in vivo. Concurrently, cell death occurs in an immunogenic fashion, leading to immune-mediated tumor reduction at a distant site not exposed to gas plasma irradiation (abscopal effect) and improved survivalThe term abscopal effect was introduced in 1953,[37] describing radiotherapy of one tumor site evoking hampered growth of another, distant, untreated metastatic tumor sites. This phenomenon has occasionally been reported in the medical literature and speculated to be immune-mediated.[38] Today, it is well recognized that the immune system makes a significant contribution to tumor control, as exemplified by the clinical application of checkpoint immunotherapy.[39] This understanding also sheds new light on the abscopal effect observed with predominantly local treatment, such as radiotherapy.[40] Many studies using experimental tumor models have reported abscopal effects, for instance, in pancreatic cancer,[41] malignant melanoma,[42] colon cancer,[43] lung cancer,[44] and breast cancer .[45] The data of our experimental model suggest gas plasma irradiation to mediate an abscopal effect, too. The distant tumor lesion responding to the local therapy usually is characterized not only by retarded growth but also an enhanced immune infiltrate, such as T-cells.[46] Especially cytotoxic CD8+ T-cells were in focus in recent years when it came to harnessing the power of antitumor T-cell responses.[47] This subtype, however, was not significantly enhanced in our study when investigating the tumor microenvironment (TME). By contrast, we found significantly elevated numbers of CD4+ T-cells in the TME, also on the gas plasma irradiated tumor site. This is of note as we previously determined T-cells to be highly susceptible to gas plasma irradiation-induced oxidative stress[48-50] in vitro, which was predominantly found in antigen-experienced memory T-cells.[51] Hence, gas plasma irradiation-induced T-cell death either is of minor relevance in vivo or not resolved kinetically as new T-cells constantly migrate into the TME from the circulation. In the current studies, the CD4+ T-cells expressed pro-inflammatory IL-17, while the immunosuppressive subtype was less pronounced, as FOXP3 expression suggested. CD4+ T cells have long been underappreciated in their ability to contribute antitumor effects in the TME,[52] albeit their positive prognostic value in cancers of, for instance, the head and neck,[53,54] the lung,[55] and the breast.[56] The value of CD4+ T-helper cells lies in their supportive interplay with other types of immune cells, such as cytotoxic CD8+ T-cells that subsequently exert the antitumor effects in experimental models.[57,58] The polyfunctional role of CD4+ T-cells is even extended to direct cytotoxic action as premature evidence points into mechanism related to MHCII-expression on tumor cells.[59,60] CD4+ T-cells, moreover, promote the antitumor efficacy of dendritic cells (DCs),[61] a cell type we also found to be increased in the gas plasma-irradiated breast cancers, as compared to controls.DCs are professional antigen-presenting cells that respond to damage-associated molecular patterns (DAMPs).[62] DAMPs are intimately linked to the immunogenic cancer cell death (ICD),[35] a pro-inflammatory cell death program that overrides the predominantly immunosuppressive nature of many types of regulated cell death.[63] In this regard, we found the DAMPs and ICD-associated molecules HSP70, HSP90, calreticulin (CRT), and MHC-I, as well as ATP secretion, IFNα2, IFNγ, and the interleukin 6 (IL-6) to be upregulated following gas plasma irradiation. Also, the programmed death-ligand 1 (PD-L1) which contributes to a muted immune response was found to be upregulated as it is also known for radiation therapy.[64,65] However, this might be outbalanced by the pro-inflammatory molecules that were increasingly present as well, such as CRT. CRT levels depend on the type of chemotherapeutic agent used and are associated with elevated rejection rates of tumors in vivo.[28] Intriguingly, also localized therapies such as ionizing radiation and photodynamic therapy (PDT) led to increased CRT levels in tumor cells and subsequently hampered tumor growth in vivo that was associated with increased immune infiltrates into the TME.[66-69] The effects of both PDT and radiotherapy are known to be mediated significantly via reactive oxygen species (ROS).[70-72] ROS, in turn, have significant implications in the release and appearance of DAMPs and ICD in terms of not only oxidative damage but also redox regulation.[73-76]Medical gas plasma systems are a novel technology for oncotherapy. Gas plasmas expel a multitude of different ROS and RNS onto the irradiated target in a localized fashion.[25] In contrast to radiotherapy and PDT, these reactive species are not generated inside the target cells but are being deposited from the extracellular gas plasma phase. Our results showing gas plasma irradiation-mediated cytotoxicity in breast cancer cells are in line with a previous study reporting increased apoptosis and cell damage related to antioxidant profiles and degradation of PKD2.[77-79] Using a microcapillary plasma jet, a limitation of breast cancer growth in vivo and an increase in the number of TUNEL+ cells was also noted before.[80] Interestingly, gas plasma irradiation was shown to induce hallmarks of immunogenic cell death before in malignant melanoma cells,[81,82] lung cancer,[83] and colorectal cancer.[84,85] Some reports also observed ICD-induction in tumor cells exposed to cell culture medium[86] or phosphate-buffered saline[87] previously irradiated with gas plasma, implicating that some of the long-lived ROS may contribute to effects observed.[88] These treatments were also associated with higher tumor immuno-infiltrates in vivo when applied as intraperitoneal lavage to target peritoneal carcinomatosis of colorectal or pancreatic cancer origin.[87,89] In vivo, direct gas plasma-induced ICD was not only found to correlate with anticancer efficacy and immune infiltration[85,90,91] but also showed protection of subsequent re-challenge using the syngeneic vaccination model with gas plasma-inactivated tumor cells.[31,81] Along with our findings, these results suggest gas plasma irradiation to be a novel and potent tool for promoting antitumor immunity by eliciting ICD in tumor cells. An open question, however, remains, which single or composite of reactive species in terms of quantity and quality contribute to gas plasma irradiation-induced ICD. Among the species being generated, especially in our system, is, for instance, hydroxyl radicals, superoxide anion, nitrite, singlet delta oxygen, hydrogen peroxide, nitrate, atomic oxygen, nitric oxide, ozone, peroxynitrite, and several types of metastables, each having distinct spatio-temporal distributions and reaction kinetics with other molecules.[92-96] While it is known that the multi-ROS nature of gas plasmas exceeds the extent of ICD when compared to individual species,[81] it is also clear the individual distribution within this ROS cocktail affects ICD-related outcomes as well.[31] The proof-of-concept of ROS-induced ICD remains exciting as a recent study had suggested increased tumor cell immunogenicity and more potent therapeutic effects of autologous DCs in patients suffering from metastatic ovarian cancer.[97,98]Novel cancer therapies not only need to be effective but also safe. While exact technical parameters and the type of feed gas differ between gas plasma systems, their unifying theme is the generation of ROS together with other plasma components, such as electrons and ions, electric fields, and UV-emission. A large body of evidence, however, points to ROS and RNS to be the primary mode of action.[25] For some gas plasma jets, extensive safety testing has been performed before.[99] This includes the absence of mutagenic effects using the OECD-accredited assays for HRPT[100] and micronucleus in vitro[101] and in vivo.[102] Moreover, a one-year follow-up study in mice[103] and human volunteers[104] did not report on any tumorigenic effects of repeated gas plasma irradiation. Moreover, gas plasma irradiation seems to target malignant over nonmalignant cells selectively.[24,105] It is suggested that this is due to faster proliferation,[106] the different lipid composition of the cell plasma membrane,[107] lower concentration of cholesterol,[108] lower levels of antioxidants,[109] and/or higher ROS concentration in cancer cells.[110,111] For breast cancer cells, it was speculated that their accelerated genome mutation rate, as well as hyper-activated MAPK/JNK and NF-κB pathways, promote their vulnerability toward gas plasma irradiation when compared to nonmalignant cells.[112]A limitation of our study was the lack of the gold standard prophylactic vaccination. The non-orthotopic model of breast cancer is only partially realistic in terms of microenvironment and cellular differentiation. The blockage of effector cell subtypes and its effect on tumor growth was not investigated but cellular subtypes and tumor immuno-infiltration was analyzed.
Authors: Sander Bekeschus; Anke Schmidt; Axel Kramer; Hans-Robert Metelmann; Frank Adler; Thomas von Woedtke; Felix Niessner; Klaus-Dieter Weltmann; Kristian Wende Journal: Environ Mol Mutagen Date: 2018-02-08 Impact factor: 3.216
Authors: Cheryl Lai-Lai Chiang; Lana E Kandalaft; Janos Tanyi; Andrea R Hagemann; Gregory T Motz; Nikolaos Svoronos; Kathleen Montone; Gina M Mantia-Smaldone; Lori Smith; Harvey L Nisenbaum; Bruce L Levine; Michael Kalos; Brian J Czerniecki; Drew A Torigian; Daniel J Powell; Rosemarie Mick; George Coukos Journal: Clin Cancer Res Date: 2013-07-09 Impact factor: 12.531
Authors: Michel Obeid; Antoine Tesniere; François Ghiringhelli; Gian Maria Fimia; Lionel Apetoh; Jean-Luc Perfettini; Maria Castedo; Grégoire Mignot; Theoharis Panaretakis; Noelia Casares; Didier Métivier; Nathanael Larochette; Peter van Endert; Fabiola Ciccosanti; Mauro Piacentini; Laurence Zitvogel; Guido Kroemer Journal: Nat Med Date: 2006-12-24 Impact factor: 53.440
Authors: Michael Rückert; Lisa Deloch; Rainer Fietkau; Benjamin Frey; Markus Hecht; Udo S Gaipl Journal: Strahlenther Onkol Date: 2018-03-02 Impact factor: 3.621
Authors: Sander Bekeschus; Kristian Wende; Mohamed Mokhtar Hefny; Katrin Rödder; Helena Jablonowski; Anke Schmidt; Thomas von Woedtke; Klaus-Dieter Weltmann; Jan Benedikt Journal: Sci Rep Date: 2017-06-05 Impact factor: 4.379
Authors: Anke Schmidt; Thomas von Woedtke; Jan Stenzel; Tobias Lindner; Stefan Polei; Brigitte Vollmar; Sander Bekeschus Journal: Int J Mol Sci Date: 2017-04-19 Impact factor: 5.923
Authors: Abraham Lin; Joey De Backer; Delphine Quatannens; Bart Cuypers; Hanne Verswyvel; Edgar Cardenas De La Hoz; Bart Ribbens; Vasiliki Siozopoulou; Jonas Van Audenaerde; Elly Marcq; Filip Lardon; Kris Laukens; Steve Vanlanduit; Evelien Smits; Annemie Bogaerts Journal: Bioeng Transl Med Date: 2022-04-21
Authors: Hager Mohamed; Eric Gebski; Rufranshell Reyes; Samuel Beane; Brian Wigdahl; Fred C Krebs; Katharina Stapelmann; Vandana Miller Journal: Cancers (Basel) Date: 2021-05-18 Impact factor: 6.639