Spontaneous fusion between squamous cell carcinoma and myeloid cells
SCCVII/SF, a squamous cell carcinoma cell line derived from C3H mice, was examined for suitability for tumor-myeloid fusion studies. SCCVII/SF-H2B-GFP cells, which express a nuclear GFP signal, were grown in the dorsum of the syngeneic mice, and the resulting tumor was excised and processed for single-cell suspensions. Since hybrid cells derived from cancer-myeloid cell fusion should carry traits from both parental cells, flow cytometry was used to gate the tumor cells based on their GFP expression and myeloid markers CD45 and CD11b. ~1% of the GFP + cells expressed the myeloid markers, which reflected the fusion rate. An example of a 1.77% fusion rate is shown here (Fig. 1a, left panel). Since phagocytosis of SCCVII/SF-H2B-GFP cells by myeloid cells could produce “pseudo-hybrids,” which would have weak GFP from tumor cell remnants in the phagosomes of the myeloid cells, only cells with high levels of GFP and myeloid markers were gated for detection of true hybrids. The images of the hybrid cells were obtained with an imaging flow cytometry to document nuclear GFP and membranous CD11b (Fig. 1a, right panel). Such specific spatial distribution of the markers provided convincing evidence that fusion between SCCVII/SF and myeloid cells occurred in vivo.
Bone marrow macrophages (BMM) require M-CSF to survive and RANKL to fuse and form multinucleated osteoclasts 24. The properties of these cytokines were exploited to promote fusion between BMM and tumor cells in vitro. When M-CSF and RANKL were added to the co-culture of BMM and SCCVII/SF-H2B-GFP, the resulting multinucleated giant cells expressed not only macrophage marker PU.1, but also GFP from the tumor cells. Notably, while the PU.1 + nuclei may or may not express GFP, all the GFP + nuclei in the giant cells expressed PU.1 (Fig. 1b), indicating the occurrence of nuclear fusion between the two cell types. These multinucleated giant cells were able to survive for more than 1 month with no sign of death as long as the media was replenished. Multinucleated giant cells formed by mouse BMM mono-cultures in the same media conditions usually died after 10 days.
SCCVII/SF-BMM hybrid cell clones exhibit variable tumorigenicity
Most studies show that fusion with macrophages promotes tumor aggression 6,7 4,8−10. However, tumor cells fused with mesenchymal cells or stem cells have been reported to induce anti-tumor events 11 3,25. To study the outcome of fusion between SCCVII/SF and macrophages, we used the same culture conditions described above with some modifications to generate the hybrid cells for subsequent studies. Briefly, BMM were retrovirally transduced to express antibiotic resistance and cytoplasmic GFP before being co-cultured with SCCVII/SF in the presence of M-CSF and RANKL. After the culture formed multinucleated giant cells, the cytokines were withdrawn and antibiotic selection was started to enrich the hybrid cells, which were later cloned. The surviving cells acquired GFP expression and antibiotic resistance from the BMM but had only one nucleus and could survive indefinitely without the supporting cytokines.
Six clones of hybrid cells were injected into the posterior dorsum of mice. SCCVII/SF-GFP, which was SCCVII/SF transduced with the GFP retrovirus used for the hybrid cells, was used as the control. The hybrid cell clones showed variable growth abilities in mice. Clones D2 and D3 were selected for further characterization as they produced tumors that were significantly larger or smaller/not present, compared to those produced by the control cells (Fig. 2a). Most of the mice that received the D2 cells developed significantly larger tumors even though some were underestimated due to loss of tumor mass by ulceration and/or animal premature death (Fig. 2b). Unexpectedly, ~ 16% of the mice injected with the D2 cells failed to develop grossly visible tumors within the study period of 3 weeks. In contrast, the D3 cells either failed to form tumors or produced very small tumors even when the study period was extended. The D3 line never produced large tumors. The growth disparities of these two hybrid cell clones were also seen when grown in the immunocompromised NSG mice, where they grew faster, forcing the study period to be shortened to 16 days to avoid major adverse events (Fig. 2c & d).
SCCVII/SF-BMM hybrid cells have a low growth rate in vitro but exhibit increased ability to recruit stromal cells in vivo
The hybrid cell clones D2 and D3 were analyzed by a growth curve experiment in cell culture. Both cells exhibited slower growth and prolonged doubling times (D2 ~ 23.97 hours, D3 ~ 28.31 hours) than the control cells (~ 16.79 hours) after culturing for four days (Fig. 3a). The D3 cells exhibited the slowest growth, consistent with their diminished abilities to generate tumors in mice. Paradoxically, the D2 cells also showed a slow growth kinetic despite their ability to generate large tumors in mice. Unlike cell culture, the in vivo environment is complex, as non-tumor elements could also influence tumor growth and behavior. The compositions of the D2 and control SCCVII/SF-GFP tumors were examined by flow cytometry performed on tumor cell suspensions. D3 cells were not included as they either failed to form tumors or formed tumors too small to generate sufficient numbers of cells for flow cytometry. Based on the GFP expression of the D2 and control cells, the tumor cells (GFP+) were gated from the non-tumor stromal cells (GFP-) by flow cytometry.
The results showed that the D2 tumors exhibited a higher percentage of GFP- cells (Fig. 3b, left panel) than the control. As a result, a significantly lower ratio between the tumor and stromal cells, expressed as GFP+/GFP- was observed (mean value of control vs. D2: ~0.73 vs. ~0.26) (Fig. 3b, right panel), indicating that only a small portion of the D2 tumor was composed of the D2 hybrid cells. In other words, despite the large size, the bulk mass of the D2 tumor was formed not by the tumor cells but by the stromal cells. Most of the stromal cells were expected to be fibroblasts since they are the most abundant cell type in solid tumors 26.
Flow cytometry was also used to examine the cell cycle distributions of the cells in the tumors. The results showed that the D2 tumor had a significantly higher percentage of tumor cells (GFP+) in the G2/M phase while a lower percentage in the G1/G0 phase compared to the control tumors (mean value of control vs. D2: G2/M: ~24% vs. ~51%. G1/G0: ~63% vs. ~41%) (Fig. 3c & d). The stromal cells (GFP-) from both tumors demonstrated similar cell cycle distributions, with the majority (~ 84%) of the cells were in G0/G1 while only a minority (~ 5%) transitioned through the G2/M phase (Supplementary Fig. S1).
Fusion with BMM increases SCCVII/SF’s ability to recruit immune cells
Tumor-infiltrating immune cells play essential roles in cancer development, and the most abundant immune cells in solid tumors are macrophages, which are present through all stages of tumorigenesis 27. Macrophages can polarize to M1 or M2, a tumoricidal or immunosuppressive phenotype, respectively 28. In established tumors, tumor-associated macrophages (TAM) often polarize to M2 to negatively modulate the tumor immune microenvironment. T regulatory cells (Treg), cytotoxic T cells (Tc), and NK cells are also important tumor-infiltrating immune cells. Similar to macrophages, Treg cells also have an immunosuppressive function.
Flow cytometry was used to gate the tumor cell suspensions with GFP and macrophage markers (CD45 + CD11b + F4/80+) to determine the ability of the tumor cells to recruit macrophages (Fig. 4a), expressed as the TAM/GFP + ratio. The results showed that the D2 tumor had a significantly higher ratio of TAM/GFP+ (mean value of control vs. D2: ~ 0.7 vs. 1.8), indicating that the D2 cells were more capable of recruiting TAM than the control cells (Fig. 4b, upper panel). There were no differences between the two tumors in the percentages of TAM in total leukocytes (CD45+) (Fig. 4b, lower panel). As the D3 tumor was too small to be analyzed by flow cytometry, F4/80 IHC was used to stain TAM in the tumor sections to provide a qualitative assessment. The result showed that the D3 tumor, similar to the D2 tumor, also harbored a substantial number of TAM (Fig. 4c).
Furthermore, TAM polarization was examined by flow cytometry with macrophage markers for M1 (CD11b + F4/80 + iNos+) and M2 (CD11b + F4/80 + Arg1+) (Fig. 4d, 4e). The results showed ~ 62% (mean value) of TAM in the D2 tumor and only ~ 26% of TAM in the control tumors polarized to M2 (Fig. 4e upper panel). In addition, M2 macrophages constituted ~ 14% of the total cell population in the D2 tumor while only 5% was observed in the control tumor (Fig. 4e lower panel). In contrast, M1 macrophages were not significantly present, and their numbers were not different between the two tumors (Fig. 4e lower panel).
The T cell compositions in the tumors were also analyzed by flow cytometry, including Treg (CD4 + CD25 + FoxP3), Tc (CD45 + CD8+), and Th (CD45 + CD4+) (Figs. 5a, 6a). Similar to the TAM results, the D2 tumor also exhibited an increased ability to recruit Treg and Tc, with higher numbers of these immune cells per total cells (Figs. 5b, 6b, top and middle panels). The D2 and the control tumor showed no differences in the percentages of Treg/total helper T cells (CD4+) (Fig. 5b, lower panel). In comparison, the D2 tumor had a slight increase in the percentage of Tc/total leukocytes (CD45+) (Fig. 6b, lower panel). As the D3 tumor was too small to be analyzed by flow cytometry, FoxP3 or CD8 IHC was used to provide a qualitative assessment of Treg or Tc in the tumor. The results showed that the D3 tumor, similar to the D2 tumor, also harbored higher numbers of Treg and Tc compared to the control (Fig. 5c & 6c).
The numbers of tumor-infiltrating Treg and Tc cells were less than 1/10 of the number of M2 macrophages in both tumors (Figs. 5B, 6B, middle panel). There was a scant presence of NK cells (not detectable), expected for SCC-VII/SF-derived tumors, which exhibit NK cell suppression due to high levels of TGF-β 29,30.
Tumor angiogenesis
Tumor growth is highly dependent on neovascularization. In this study, the tumors derived from the hybrid cells recruited high numbers of M2 macrophages. Since M2 macrophages are important pro-angiogenic cells 31, we investigated whether their presence impacted the micro-vessels in the tumors. CD31 IHC staining in endothelial cells showed the micro-vessels in the tumors (Fig. 7a, top panel). The images were segmented (Fig. 7a, lower panel) before being analyzed by the Microvessel-Segmentation MATLAB plugin 32.
The results show that tumors formed by the hybrid cells, D2 and D3, exhibited higher microvessel densities (MVD) and thinner micro-vessel wall thickness (MVWT) than the control tumor (Fig. 7b). There were no significant differences in other angiogenesis parameters, such as blood vessel sizes and percentage of vessel area.
Differential RNA expression
RNA-seq was used to investigate the mechanisms responsible for the behavioral differences between the hybrid and control cells. Ingenuity Pathway Analysis (QIAGEN) was used to generate differential RNA expression profiles, and the unfiltered results are provided in supplementary Excel file “Upstream regulatory analysis of RNA-seq.xlsx”. The file contains: Sheet 1 − D2 hybrid cells vs. SCCVII/SF-GFP; sheet 2 − D3 hybrid cells vs. SCCVII/SF-GFP; sheet 3 − D2 hybrid cells vs. D3 hybrid cells.
Table 1 lists the representative upstream regulatory pathway results relevant to the current study. Compared with the control SCCVII/SF-GFP cells, the D2 cells exhibited increased activation in pathways involving five cytokines - tumor necrosis factor α (TNFα), interleukin (IL)-1α & β, interferon γ (IFNγ), and oncostatin M (OSM), along with two growth factors - epidermal growth factor (EGF) and vascular endothelial growth factor (VEGF)-α. All five cytokines are pro-inflammatory, and upregulation of their activities likely contributed to the increased numbers of inflammatory cells in the D2 tumor. EGF is a mitogen for many cell types, including fibroblasts and epithelial cells 33,34 while TNFα also plays an important role in the maintenance and function of cancer-associated fibroblasts 26,35. Thus, upregulation of TNFα and EGF pathways likely contributed to the significant presence of non-tumor cells that resulted in large sizes of the D2 tumor.
Table 1
Representative upstream regulators that influence the biological behaviors of the hybrid cells. Only pathways of cytokines, growth factors, and major cell cycle molecules are listed. The pathways predicted to be activated are listed on the left side of the z-score column, while those predicted to be inhibited are listed on the right side. A z-score ≥ 2 predicts an activation state, while a z-score ≤ -2 predicts an inhibitory state of the pathway.
RNA-Seq | Category | Upstream Regulator | Activation z-score | p-value of overlap |
D2 vs SCC | Cytokine | TNF | 3.663 | 1.77E-11 |
IL1b | 2.837 | 0.000178 |
IFNg | 2.811 | 0.0312 |
OSM | 2.211 | 0.0156 |
IL1a | 2.007 | 0.00542 |
Growth factor | VEGFa | 2.771 | 0.000839 |
EGF | 2.585 | 0.0374 |
Cell Cycle | TP53 | 2.056 | 0.000185 |
D3 vs SCC | Cytokine | TNF | 2.018 | 1.85E-11 |
Growth factor | VEGFa | 2.244 | 0.00153 |
FGF7 | -2.169 | 0.00779 |
AREG | -3.352 | 9.53E-06 |
Cell Cycle | CDKN1A | 2.532 | 0.0000216 |
TP53 | 2.928 | 1.34E-12 |
CDKN2A | 2.07 | 0.00735 |
D2 vs D3 | Cytokine | IL6 | 2.021 | 0.0126 |
Growth factor | AREG | 2.917 | 8.3E-10 |
FGF7 | 2.028 | 0.0114 |
Cell Cycle | CDKN1A | -2.266 | 0.000000572 |
TP53 | -3.37 | 8.49E-14 |
CCNK | 2.366 | 0.00246 |
The D3 cells also exhibited upregulation of the TNFα pathway, which played important roles in diverse cellular events in tumor, including recruitment and modulation of TAM as well as angiogenesis 35. The D3 cells had decreased signaling from fibroblast growth factor (FGF)-7 and amphiregulin (AREG), the biological implications of which are unclear.
VEGF-α is an essential factor for angiogenesis 36. Upregulation in the VEGF-α pathway increases proliferation and migration of vascular endothelial cells. This process could be enhanced by inflammation, which likely contributed to the increase in tumor vasculature observed in the D2 and D3 tumors.
The RNA-seq results also showed that both D2 and D3 cells exhibited TP53 and cell cycle checkpoint activation. This suggests that during the process of nuclear merging, cells lost chromosomes and developed aneuploidy, leading to TP53 activation 37,38. These events were likely responsible for the altered growth kinetics and cell cycle distribution. When compared with the control, the D3 cells showed induction of cyclin-dependent kinase inhibitor 1A (CDKN1A) and cyclin-dependent kinase inhibitor 2A (CDKN2A). When compared between the D2 and D3 cells, the D3 cells exhibited a higher level of TP53 activation and consequently CDKN1A induction, which was presumably associated with their significant presence in the G0/G1 phase (Fig. S1b) and minimal tumorgenicity (Fig. 2). In contrast, the D2 cells exhibited cyclin K (CCNK) induction, which presumably maintained the genome integrity and allowed for cell cycle progression 39.