Chronic inflammation induced by persistent microbial infection plays an essential role in tumor progression. Although it is well documented that Epstein–Barr virus (EBV) infection is closely associated with nasopharyngeal carcinoma (NPC), how EBV-induced inflammation promotes NPC progression remains largely unknown. Here, we report that tumor infiltration of tumor-associated macrophages (TAM) and expression of CCL18, the cytokine preferentially secreted by TAM, closely correlate with serum EBV infection titers and tumor progression in two cohorts of NPC patients. In vitro, compared with EBV NPC cell lines, EBV+ NPC cell lines exhibited superior capacity to attract monocytes and skew them to differentiate to a TAM-like phenotype. Cytokine profiling analysis revealed that NPC cells with active EBV replications recruited monocytes by VEGF and induced TAM by GM-CSF in an NF-κB–dependent manner. Reciprocally, TAM induced epithelial–mesenchymal transition and furthered NF-κB activation of tumor cells by CCL18. In humanized mice, NPC cells with active EBV replications exhibited increased metastasis, and neutralization of CCL18, GM-CSF, and VEGF significantly reduced metastasis. Collectively, our work defines a feed-forward loop between tumor cells and macrophages in NPC, which shows how metastatic potential can evolve concurrently with virus-induced chronic inflammation. Cancer Res; 77(13); 3591–604. ©2017 AACR.

Microbial infection causally accounts for more than 20% of malignancies worldwide (1). Epstein–Barr virus (EBV) is a human gammaherpesvirus that infects approximately 95% of the global population and linked to various types of epithelial and lymphoid cancer, including nasopharyngeal carcinoma (NPC; ref. 2). Elevated EBV infection titers in serum, induced by environmental stress like smoking and other chemicals like phorbol-12-myristate-13-acetate (TPA; ref. 3) and sodium butyrate (4), is closely associated with NPC progression (5).

Tumor initiation and progression caused by viral infections reflect properties that are intrinsic to tumor cells, including expression of viral oncoprotein and nontranslated viral RNAs that directly contribute to tumor growth (6), metastasis (7), and drug resistance (8). Viral carcinogenesis can also occur indirectly. In many cases, although the immune response to pathogens that establish persistent infections is designed to promote host defense followed by tissue repair, it can also stimulate chronic inflammation and tumor growth (9). Mounting evidence has indicated that persistent viral infections can foster chronic inflammatory microenvironments that extrinsically promote tumor development (9, 10). For example, Hepatitis B and C viruses can lead to the chronic production of protumor cytokines via the induction of cell death (11). In addition, inflammation can also be induced by virus-encoding oncoproteins. For example, the oncoprotein encoded by human papillomavirus (HPV) in cervical epithelial cells can induce a hyperinflammatory response by promoting IL17A production (12). Although it has been well documented that EBV plays a central role in the pathogenesis of NPC, the interaction between EBV activity and the tumor-promoting inflammatory microenvironment in NPC remains elusive.

Macrophages are the most abundant recruited inflammatory cells in the tumor microenvironment, and they are closely associated with poor prognosis in various tumor types (13–16). Chemokine (C-C motif) ligand (CCL)18 has been shown to be exclusively secreted by myeloid-derived cells and is one of the most abundantly produced cytokines of tumor-associated macrophages (TAM) in various types of cancer, including breast, gastric, and ovarian cancer (17–19). TAM-derived CCL18 has been reported to promote cancer cell invasion by inducing the epithelial–mesenchymal transition (EMT; refs. 17, 20). However, the clinical significance of CCL18 and its underlying mechanism of action in NPC remain unknown.

In this study, we investigated the contribution of EBV replication in NPC cells to macrophage activation and the interaction between NPC cells and TAMs. The results demonstrated the underlying mechanism linking EBV replication and tumor metastasis.

Clinical specimens and study design

We obtained 580 pathologically proven and nondistant metastatic paraffin-embedded nasopharyngeal carcinoma specimens for this study. A total of 321 specimens were obtained from the Sun Yat-sen University Cancer Center (Guangzhou, China) between Jan 1, 2010 and Dec 31, 2011, and 259 specimens were collected from the Sun Yat-Sen Memorial Hospital (Guangzhou, China) between Jan 1, 2009 and Dec 31, 2011. All samples were pathologically reassessed by two pathologists. None of the patients had received radiotherapy or chemotherapy before biopsy sampling. All MRI and CT scans were reassessed separately by two radiologists, and any disagreements were resolved by consensus. Clinical staging was reclassified according to the criteria of the American Joint Committee on Cancer (AJCC)/Union for International Cancer Control (UICC)-staging manual (seventh edition). All patients were treated with 2.0 Gy per fraction with five daily fractions per week for 6–7 weeks using either two-dimensional (2D) radiotherapy (2DRT), IMRT, or three-dimensional (3D) conformal radio therapy (3DCRT), in accordance with the treatment policy adopted by each center. The cumulative radiation doses were 66 Gy or greater to the primary tumor and 60–66 Gy to the involved neck area. All potential sites of local infiltration and bilateral cervical lymphatics were irradiated to 50 Gy or greater. Patients with stage III to IVB disease also received concurrent platinum-based chemotherapy. This study was approved by the institutional ethical review boards of both hospitals, and written informed consent was obtained from all patients. The studies were conducted in accordance with International Ethical Guidelines for Biomedical Research Involving Human Subjects (CIOMS).

We used the 321 specimens from Sun Yat-sen University Cancer Center as a training set to obtain the optimal cut-off point of TAM count or CCL18+ TAM count per field for disease-free survival (DFS), overall survival (OS), and disease-specific survival (DSS). To assess whether TAM count or CCL18+TAM count also had prognostic value in different patient population, we used the 259 samples from Sun Yat-sen Memorial Hospital as an independent set to externally validate the results.

DSS was defined as the time from surgery to death from NPC or death due to treatment; patients known to be alive or lost to follow-up at the time of analysis were censored at their last follow-up and deaths from other causes without disease were treated as a competing risk. Overall survival (OS) was defined as the time from surgery to death from any cause; patients known to be still alive at the time of analysis were censored at the time of their last follow-up. DFS was defined as the time from surgery to first evidence of recurrent/metastatic disease or death from any cause; patients known to be alive without recurrent/metastatic disease or lost to follow-up at the time of analysis were censored at the time of their last follow-up.

Examination of EBV VCA-IgA antibody titers in blood sample of NPC patients

A total of 5–10 mL of each patient's blood was collected. The VCA-IgA antibody titers of EBV were measured using a commercial kit (Zhongshan Bio-tech Co Ltd) based on standard immunoenzymatic techniques, which were described previously (5). To minimize bias, samples were blinded to laboratory personnel. Serial dilutions of quality control sera (1:10, 1:20, 1:40, and 1:80) were applied to every assay for the evaluation of intraset variability. To minimize experimental error, all tests were conducted in the same laboratory by the same technician staff.

Double immunostaining for patients' specimens

Paraffin-embedded samples were sectioned at 4-μm thickness. Antigen retrieval was performed using a pressure cooker for 30 minutes in 0.01 mol/L citrate buffer (pH 6.0). For immunofluorescence staining, the specimens were incubated with antibodies specific for CD68 (catalog no. sc7083, Santa Cruz Biotechnology) and CCL18 (catalog no. MAB394, R&D Systems) overnight at 4°C, followed by incubation with Alexa Fluor 555 donkey anti-goat IgG (H+L) and Alexa Fluor 488 donkey anti-mouse IgG (H+L; Life Technologies). Cells stained with the indicated antibody were imaged using a confocal laser-scanning microscope (Carl Zeiss) with a core data acquisition system (Applied Precision) and calculated per field of view, and at least 10 fields of view per section were evaluated at 400× magnification.

Cell culture

All the NPC cell lines used in this study were kindly provided by Dr. Musheng Zeng from the Sun Yat-Sen University Cancer Center in 2014. All cell lines were passaged less than 10 times or 6 months after initial revival from frozen stocks. The cell lines were tested negative for mycoplasma contamination (catalog no. 201203, Qiagen). All cell lines were authenticated by short tandem repeat profiling prior to use. Two poorly differentiated NPC cell lines (CNE2 and HNE1) and a highly differentiated NPC cell line (HK1) were grown in RPMI1640 medium (Life Technologies) containing 10% FCS (Life Technologies) in a humidified 5% CO2 incubator at 37°C. EBV+ NPC cells were generated from parental cell lines by coculturing of the surface IgG crosslinked EGFP-neo'EBV–infected Akata cells according to the previous report's protocol (21).

Generation of humanized mice and transplantation of NPC cells

Fresh human cord blood was obtained from Sun Yat-Sen Memorial Hospital, according to guidelines approved by the ethics boards and the Clinical Research Committee at Sun Yat-Sen Memorial Hospital. The hematopoietic stem cells (HSC) were isolated by using direct CD34 Progenitor Cell Isolation Kit (Miltenyi Biotec) as described previously (20). All animal work was conducted in accordance with a protocol approved by the Institutional Animal Care and Use Committee at the medical college of Sun Yat-Sen University. Humanized mice were generated as described previously (20). Briefly, NOD/SCID mice (3–4 weeks old) were subjected to 200 cGy total body irradiation 12 hours before being injected with 2 × 105 HSCs in 0.2 mL of medium via the tail vein. At 8 weeks after HSC transplantation, 5 × 106 CNE2-EBV, CNE2-EBV+, or CSE-treated CNE2-EBV+ NPC cells were subcutaneously injected into the back of mice. For antibody treatment, mice were injected with CCL18-specific neutralizing antibody (catalog no. ab9849; Abcam) and/or GM-CSF–specific neutralizing antibody (catalog no. MAB215; R&D Systems) and/or VEGF-specific neutralizing antibody (catalog no. ab9570; Abcam) via the tail vein at 10 mg/mouse twice a week after the xenografts became palpable (about 0.4 cm in diameter). Tumor growth was evaluated by monitoring tumor volume (length × width2 × 0.5) every 3 days. The animals were sacrificed when the xenografts reached 1.5 cm in diameter. The tumor xenografts and lungs of the sacrificed mice were harvested for further investigation. The xenografts were collected to exam the state of EMT and macrophages infiltration, while the lungs were harvested to exam the metastasis by IHC staining as described before (22, 23). Similar animal experiments were performed with HNE1 cell line.

Statistical analysis

All statistical analyses were performed using SPSS for Windows version 13.0 (SPSS). We selected the optimum cut-off score for the count of CCL18+ TAM in nasopharyngeal carcinoma using X-tile plots based on the association with the patients' DSS. X-tile plots provide a single and intuitive method to assess the association between variables and survival, which can automatically select the optimum data cut-off point according to the highest χ2 value (minimum P value). We did the X-tile plots using the X-tile software version 3.6.1 (Yale University School of Medicine, New Haven, CT). Spearman correlation and regression analysis was performed to assess the relationship between EBV titres in the serum of patients and CCL18+ counts, CD68+ counts, and CCL18+TAM counts in human nasopharyngeal tissue. Kaplan–Meier survival curves were plotted, and the log-rank test was performed. We used the Cox regression model to do the multivariable survival analysis. The mean level of serum EBV VCA-IgA was calculated by geometric mean. Groups of discrete variables were compared using the Mann–Whitney U test or Kruskal–Wallis nonparametric ANOVA. All cell culture experiments were performed independently in triplicate at least three times. In all cases, P values <0.05 were considered statistically significant.

TAM infiltration and their CCL18 expression are correlated with EBV infection and poor survival in NPC patients

A total of 580 patients with NPC were included in this study, including 321 patients from Sun Yat-sen Cancer Center as the training cohort and 259 patients from Sun Yat-sen Memorial Hospital as the validation cohort. The median follow-up time of the training cohort was 43 months (range, 8–58 months), among which 98 patients (30.5%) developed local recurrences and/or distant metastases. The median follow-up time of the validation cohort was 67 months (range, 7–79 months), among which 84 patients (32.4%) had tumor relapse. The baseline demographic and clinical characteristics among two groups showed no significant statistical difference (Supplementary Table S1).

To investigate the clinical significance of macrophage infiltration and CCL18 expression in NPC, we performed double immunostaining of CCL18 and CD68 in the tumor samples. We observed that CCL18 expression was predominantly in CD68+ TAMs (Fig. 1A), which is consistent with the previous reports in other cancer types (17, 24) and 63.49% ± 12.21% of TAMs in the NPC samples were positive for CCL18. Next, we correlated the total number of macrophages (CD68+ cells) and the number of CCL18+ TAMs (CCL18+CD68+ cells) with the prognosis of NPC patients. The cut-off points of total TAM count and CCL18+ TAM count for patient outcome were determined by X-tile analysis in the training set. Total TAM count cut-off points for DFS, OS, and DSS were 22.20, 23.20, and 22.00 cells per field, while CCL18+ TAM count cut-off points for DFS, OS, and DSS were 20.70, 22.30, and 20.70 cells per field, respectively (Supplementary Table S2). We selected the cut-off points for DSS (≤ versus > 22.00 cells/field for total TAM number; ≤ versus > 20.70 cells/field for CCL18+TAM count;) as the uniform cutoff points for DFS, OS, and DSS in both training cohort and validation cohort.

Figure 1.

TAM infiltration and their CCL18 expression are correlated with EBV infection and poor survival in NPC patients. A, Representative double immunofluorescence staining of CD68 (red) and CCL18 (green) in NPC samples with low (top) or high (bottom) CCL18 expression. The isotype antibodies for CCL18 and CD68 are shown as negative controls. Scale bar, 20 μm. B, Kaplan–Meier survival curve of DFS, OS, and DSS in two cohorts. NPC patients with low (≤20.70 cells/field; training set: n = 180; validation set: n = 130) and high (>20.70 cells/field; training set: n = 141; validation set: n = 129) CCL18+ TAM counts (P value by log-rank test; HR was calculated by univariate Cox regression analysis). The numbers of surviving patients in each follow-up period are indicated below the graph. C, The count of CD68+ cell, CCL18+ cell, and CCL18+CD68+ cell in NPC samples was positively correlated with EBV titers in serum in the training cohort [left, correlation analysis: n = 321, Spearman R and P < 0.001; right, the CD68+ cell, CCL18+ cell, and CCL18+CD68+ cell count in the EBV+ (n = 288)/EBV (n = 33) group]. D, Smokers (n = 309) have higher serum concentration of total EBV viral capsid antigen-IgA than nonsmokers (n = 271) in entire cohort (geometric mean ± SE; **, P < 0.01 using Mann–Whitney U test nonparametric analysis). In the training cohort, smokers (n = 173) have higher CD68+ cell, CCL18+ cell, and CCL18+CD68+ cell count compared with nonsmokers (n = 148; mean ± SEM; **, P < 0.01; ***, P < 0.001; two-tailed Student t test).

Figure 1.

TAM infiltration and their CCL18 expression are correlated with EBV infection and poor survival in NPC patients. A, Representative double immunofluorescence staining of CD68 (red) and CCL18 (green) in NPC samples with low (top) or high (bottom) CCL18 expression. The isotype antibodies for CCL18 and CD68 are shown as negative controls. Scale bar, 20 μm. B, Kaplan–Meier survival curve of DFS, OS, and DSS in two cohorts. NPC patients with low (≤20.70 cells/field; training set: n = 180; validation set: n = 130) and high (>20.70 cells/field; training set: n = 141; validation set: n = 129) CCL18+ TAM counts (P value by log-rank test; HR was calculated by univariate Cox regression analysis). The numbers of surviving patients in each follow-up period are indicated below the graph. C, The count of CD68+ cell, CCL18+ cell, and CCL18+CD68+ cell in NPC samples was positively correlated with EBV titers in serum in the training cohort [left, correlation analysis: n = 321, Spearman R and P < 0.001; right, the CD68+ cell, CCL18+ cell, and CCL18+CD68+ cell count in the EBV+ (n = 288)/EBV (n = 33) group]. D, Smokers (n = 309) have higher serum concentration of total EBV viral capsid antigen-IgA than nonsmokers (n = 271) in entire cohort (geometric mean ± SE; **, P < 0.01 using Mann–Whitney U test nonparametric analysis). In the training cohort, smokers (n = 173) have higher CD68+ cell, CCL18+ cell, and CCL18+CD68+ cell count compared with nonsmokers (n = 148; mean ± SEM; **, P < 0.01; ***, P < 0.001; two-tailed Student t test).

Close modal

In the training set, both high macrophage number and CCL18+ TAM count were correlated with early relapse, poor long-term overall survival and disease-specific survival (Fig. 1B; Supplementary Fig. S1A, top). Among them, the CCL18+ TAM counts showed higher prognostic value than total macrophage count, as indicated by the P value and HR (DFS: CCL18+ TAM: log-rank P < 0.001; HR: 2.037; 95% CI, 1.366–3.037 versus total TAM: log-rank P = 0.042; HR: 1.517; 95% CI, 1.008–2.283; OS: CCL18+ TAM: log-rank P = 0.001; HR: 2.062; 95% CI, 1.345–3.159 versus total TAM: log-rank P = 0.040; HR: 1.582; 95% CI, 1.021–2.452; DSS: CCL18+ TAM: log-rank P < 0.001; HR: 3.099; 95% CI, 1.904–5.043 versus total TAM: log-rank P = 0.022; HR: 1.738; 95% CI, 1.069–2.826). In addition, univariate Cox regression analysis showed that DSS was associated with T stage (tumor size), N stage (lymph node metastasis status), UICC stage, and CCL18+ TAM count (Table 1). Furthermore, adjusting for tumor size, node status, and UICC stage, CCL18+TAM counts remained an independent prognostic factor of survival in multivariate analysis (HR: 2.770; 95% CI, 1.669–4.595; P = 0.001, n = 321, Table 1).

Table 1.

Univariate and multivariate Cox Regression analyses of DSS in relation to clinicopathologic features

Training set (n = 321)Validation set (n = 259)
HR (95% CI)PHR (95% CI)P
Univariate analysis 
 Age (years) 
  ≤45 (ref.) vs. >45 0.930 (0.581–1.491) 0.765 1.461 (0.928–2.297) 0.101 
 Histologic grade 
  I (ref.) vs. II–III 2.561 (0.400–10.013) 0.485 0.775 (0.356–1.688) 0.521 
 T stage 
  T1–T2 (ref.) vs. T3–T4 1.695 (1.026–2.799) 0.039a 1.804 (1.121–2.902) 0.015a 
 N stage 
  N0 (ref.) vs. N1–N3 2.994 (1.488–6.022) 0.002a 5.832 (2.994–11.361) <0.001a 
 UICC Stage 
  I–II (ref.) vs. III–IV 2.195 (1.344–3.585) 0.002a 4.818 (2.539–9.144) <0.001a 
 CCL18+TAM count 
  low (ref.) vs. high 3.079 (1.892–5.009) <0.001a 3.115 (1.867–5.198) <0.001a 
 Smoke 
  Nonsmoker (ref.) vs. smoker 1.167 (0.734–1.853) 0.514 1.228 (0.777–1.940) 0.380 
Multivariate analysis 
 T stage 
  T1–T2 (ref.) vs. T3–T4 0.844 (0.487–1.461) 0.544 1.857 (1.149–3.002) 0.011a 
 N stage 
  N0 (ref.) vs. N1–N3 2.109 (1.002–4.437) 0.049a 4.217 (2.129–8.353) <0.001a 
 UICC Stage 
  I–II (ref.) vs. III–IV 1.761 (1.033–3.003) 0.038a 2.766 (1.439–5.319) 0.002a 
 CCL18+TAM count 
  Low (ref.) vs. high 2.770 (1.669–4.595) 0.001a 2.177 (1.292–3.668) 0.003a 
Training set (n = 321)Validation set (n = 259)
HR (95% CI)PHR (95% CI)P
Univariate analysis 
 Age (years) 
  ≤45 (ref.) vs. >45 0.930 (0.581–1.491) 0.765 1.461 (0.928–2.297) 0.101 
 Histologic grade 
  I (ref.) vs. II–III 2.561 (0.400–10.013) 0.485 0.775 (0.356–1.688) 0.521 
 T stage 
  T1–T2 (ref.) vs. T3–T4 1.695 (1.026–2.799) 0.039a 1.804 (1.121–2.902) 0.015a 
 N stage 
  N0 (ref.) vs. N1–N3 2.994 (1.488–6.022) 0.002a 5.832 (2.994–11.361) <0.001a 
 UICC Stage 
  I–II (ref.) vs. III–IV 2.195 (1.344–3.585) 0.002a 4.818 (2.539–9.144) <0.001a 
 CCL18+TAM count 
  low (ref.) vs. high 3.079 (1.892–5.009) <0.001a 3.115 (1.867–5.198) <0.001a 
 Smoke 
  Nonsmoker (ref.) vs. smoker 1.167 (0.734–1.853) 0.514 1.228 (0.777–1.940) 0.380 
Multivariate analysis 
 T stage 
  T1–T2 (ref.) vs. T3–T4 0.844 (0.487–1.461) 0.544 1.857 (1.149–3.002) 0.011a 
 N stage 
  N0 (ref.) vs. N1–N3 2.109 (1.002–4.437) 0.049a 4.217 (2.129–8.353) <0.001a 
 UICC Stage 
  I–II (ref.) vs. III–IV 1.761 (1.033–3.003) 0.038a 2.766 (1.439–5.319) 0.002a 
 CCL18+TAM count 
  Low (ref.) vs. high 2.770 (1.669–4.595) 0.001a 2.177 (1.292–3.668) 0.003a 

Abbreviations: DSS, disease-specific survival; CI, confidence interval.

aP < 0.05, statistical difference.

To further validate the prognostic value of the TAM number and CCL18+ TAM count, we applied the above cut-off setting to an independent validation group with 259 cases. Similarly, patients with higher TAM number or CCL18+ TAM count had shorter DFS, OS, and DSS than those with lower TAM number or CCL18+ TAM count (Fig. 1B; Supplementary Fig. S1A, bottom). Again, the CCL18+TAM counts showed higher prognostic value than total macrophage count (DFS: CCL18+ TAM: log-rank P < 0.001; HR: 2.172; 95% CI, 1.414–3.337; versus total TAM: log-rank P = 0.018; HR: 1.928; 95% CI, 1.109–3.350; OS: CCL18+ TAM: log-rank P = 0.003; HR: 1.814; 95% CI, 1.209–2.720; versus total TAM: log-rank P = 0.022; HR: 1.708; 95% CI, 1.071–2.725; DSS: CCL18+ TAM: log-rank P < 0.001; HR: 3.115; 95% CI, 1.867–5.198; versus total TAM: log-rank P = 0.017; HR: 2.067; 95% CI, 1.247–3.425). In consistent to the training cohort, Cox regression multivariate analysis showed that CCL18+TAM counts were an independent prognostic factor for patient DSS in the validation cohort (HR: 2.177; 95% CI, 1.292–3.668; P = 0.003, n = 259, Table 1).

In addition, we investigated the association between the CCL18+ TAM count and the clinicopathological features in both cohorts. A higher CCL18+ TAM count was significantly associated with a higher histologic grade, more aggressive lymph node, distant metastases, and a more advanced cancer stage in both cohorts (Supplementary Table S3). Collectively, these data suggest that the TAM-derived-CCL18 plays an important role in the progression of human NPC.

Elevated EBV infection titer in serum is closely associated with NPC progression. Therefore, we investigated whether EBV activity is correlated with macrophage infiltration and their CCL18 expression in NPC patients. Our data showed that the total TAM count, CCL18+ cell count, and CCL18+ TAM count all significantly increased in the EBV+ group compared with the EBV group in both training cohort (Fig. 1C) and validation cohort (Supplementary Fig. S1B). Moreover, EBV infection titer in serum is positively correlated with the total TAM count, CCL18+ cell count, and CCL18+ TAM count in NPC patients (Fig. 1C; Supplementary Table S3). In addition, the percentages of CCL18+ TAMs (CCL18+CD68+) in total TAMs significantly increase in the EBV+ group compared with the EBV group in both cohorts (Supplementary Fig. S1C).

It has been reported that smoking is associated with an increased EBV seropositivity among NPC patients (5). In agreement, our data showed that the serum concentration of EBV VCA IgA in smokers (n = 309, geometric mean = 675.12, SE = 37.55) was higher than that in nonsmokers (n = 271, geometric mean = 536.78, SE = 24.14, P = 0.001 by using Mann–Whitney U test, Fig. 1D). Importantly, we observed that the total TAM count, CCL18+ cell count, and CCL18+ TAM count all significantly increased in the smokers compared with the nonsmokers in both training cohort (Fig. 1D) and validation cohort (Supplementary Fig. S1D). Collectively, these data indicated that tumor infiltration of macrophages and their CCL18 expression are closely associated with EBV infection, smoking, and poor prognosis of NPC patients.

T-cell infiltration has been reported to be associated with a better prognosis in many tumors (25, 26). We evaluated T-cell infiltration in 580 NPC samples by CD3 immunostaining and observed both TAM count and CCL18+ TAM count were negatively correlated with T-cell counts (TAM count: R2 = 0.199, P < 0.001; CCL18+TAM count: R2 = 0.206, P < 0.001; Supplementary Fig. S1E), suggesting macrophage infiltration and their CCL18 expression may be associated with NPC immunosuppression.

NPC cells with different EBV replications exhibit distinctive capacity to recruit and activate macrophages

We next investigated the interaction between macrophages and NPC cells with different EBV replication in vitro by using a Transwell coculture system. Freshly isolated human monocytes were added to the top chamber, and condition media (CM) of three GFP-EBV–infected NPC cell lines (CNE2-EBV+, HNE1-EBV+, and HK1-EBV+) or that of parental EBV cells was added to the bottom chamber of the transwell apparatus. We observed significantly more monocytes migrated toward CM of EBV+ NPC cells compared with that of EBV NPC cells and control medium (Fig. 2A). To investigate whether EBV+ NPC cells can activate monocytes to become TAM-like macrophages, we placed monocytes into the bottom transwell chamber to be cocultured with EBV NPC cells or EBV+ NPC cells for 5 days and evaluated surface markers (CD206high/HLA-DRlow; ref. 20) and protumor cytokines associated with TAMs (CCL18, CCL17, CCL22, and IL10; refs. 17, 20, 27) in the cocultured macrophages. We found that only macrophages cocultured with EBV+ NPC cells, rather than those cocultured with EBV NPC cells exhibited a CD206high/HLA-DRlow phenotype compared with untreated macrophages (Fig. 2B and C; Supplementary Fig. S2A). In consistent with the changes in surface markers, macrophages cocultured with EBV+ NPC cells produced significantly higher amounts of protumor cytokines compared with untreated macrophages and macrophages cocultured with EBV NPC cells (Fig. 2D).

Figure 2.

NPC cells with different EBV replications exhibit distinctive capacity to recruit and activate macrophages. A, Monocyte chemotaxis in response to CM from EBV NPC cells or EBV+ ones with or without CSE pretreatment (CNE2, HNE1, and HK1) was determined using a transwell assay (mean ± SEM; n = 5; **, P < 0.01; ***, P < 0.001; two-tailed Student t test). B, The mean fluorescence intensity (MFI) of CD206/HLA-DR in macrophages, with indicated treatments evaluated by flow cytometry (mean ± SEM; n = 5; **, P < 0.01; ***, P < 0.001; two-tailed Student t test). C, The representative plots of CD206/HLA-DR levels in macrophages, with indicated treatments. The histograms are representative of five independent experiments of macrophages from five different donors. The values of mean fluorescence intensity (×103, mean ± SEM, n = 5) are indicated on the top-right corners of the representative plots. D, Cytokine levels in the media of macrophages, with indicated treatments (mean ± SEM; n = 5; **, P < 0.01; ***, P < 0.001; two-tailed Student t test).

Figure 2.

NPC cells with different EBV replications exhibit distinctive capacity to recruit and activate macrophages. A, Monocyte chemotaxis in response to CM from EBV NPC cells or EBV+ ones with or without CSE pretreatment (CNE2, HNE1, and HK1) was determined using a transwell assay (mean ± SEM; n = 5; **, P < 0.01; ***, P < 0.001; two-tailed Student t test). B, The mean fluorescence intensity (MFI) of CD206/HLA-DR in macrophages, with indicated treatments evaluated by flow cytometry (mean ± SEM; n = 5; **, P < 0.01; ***, P < 0.001; two-tailed Student t test). C, The representative plots of CD206/HLA-DR levels in macrophages, with indicated treatments. The histograms are representative of five independent experiments of macrophages from five different donors. The values of mean fluorescence intensity (×103, mean ± SEM, n = 5) are indicated on the top-right corners of the representative plots. D, Cytokine levels in the media of macrophages, with indicated treatments (mean ± SEM; n = 5; **, P < 0.01; ***, P < 0.001; two-tailed Student t test).

Close modal

Environmental stress like smoking can increase EBV replication. Indeed, we observed fluorescence intensity in GFP-EBV–infected NPC cell lines was enhanced in a dose-dependent manner following treatment with proportions of cigarette smoke extract (CSE) ranging from 5% to 20% (Supplementary Fig. S2B–S2D), which is consistent with a previous report (5). To further investigate whether EBV replication was associated with NPC cell's capacity to recruit and activate macrophages, we treated the monocyte/macrophages with the CM of EBVNPC cells and EBV+NPC cells with/without CSE pretreatment. The chemotaxis assay results showed that the CM of CSE-treated EBV+NPC cells attracted more monocytes compared with that of untreated EBV+NPC cells. In contrast, there is no difference between the EBVNPC cells with and without CSE treatment (Fig. 2A). In addition, macrophages cocultured with CSE-treated EBV+ cells exhibited a more TAM-like phenotype than those cocultured with untreated EBV+ cells, which was indicated by higher CD206 and lower HLA-DR expression and increased secretion of protumor cytokines. Similar to the chemotaxis assay, the activation state of macrophages has no difference between EBVcells with and without the CSE treatment (Fig. 2B–D; Supplementary Fig. S2A). Collectively, these data suggest that EBV replication in NPC cells plays a crucial role in their capacity to recruit and activate macrophages.

EBV-replicating NPC cells recruit monocytes and activate TAMs via VEGF and GM-CSF respectively

To identify the EBV+ NPC-derived cytokines that recruit and activate macrophages, the cytokine profiles of the CM from CNE2-EBV, CNE2-EBV+, and CSE-treated CNE2-EBV+ cells were analyzed using a RayBio Human Cytokine Antibody Array. The gray intensity analysis showed that four cytokines were increased in the CM of CSE-treated EBV+ cells compared with the CM of EBV cells and that of EBV+ NPC cells; including GM-CSF, VEGF, GRO, and IL1α (Fig. 3A). The ELISA assay further confirmed the increases of those cytokines in EBV+ NPC cell lines compared with EBV ones, and further increase in CSE-treated EBV+ NPC cell lines, rather than CSE-treated EBV ones (Fig. 3B; Supplementary Fig. S3A). We then investigated which of those cytokines was responsible for monocyte recruitment and TAM activation. Using different concentrations of recombinant cytokines, we found that among these four cytokines, only VEGF recruited monocytes to a level comparable with that of CSE-treated EBV+ cell CM (Fig. 3C). Moreover, only GM-CSF activated macrophages to produce protumor cytokines to levels comparable with that of CSE-treated EBV+ cell CM (Fig. 3D). Furthermore, the addition of neutralizing GM-CSF antibody to the CM of EBV+ cells and CSE-treated EBV+ cells abrogated the induction of TAM-related surface marker expression (Fig. 3E; Supplementary Fig. S3B) and cytokine secretion (Fig. 3F; Supplementary Fig. S3C), whereas the addition of neutralizing VEGF antibody had no appreciable effect on surface marker expression or cytokine production (Fig. 3E and F). In comparison, neutralization of VEGF, rather than GM-CSF, markedly inhibited the monocyte chemotaxis induced by EBV+ cells and CSE-treated EBV+ cells (Fig. 3G; Supplementary Fig. S3D). These findings indicate that EBV-replicating NPC cells recruit monocytes and activate TAMs via VEGF and GM-CSF, respectively.

Figure 3.

EBV-replicating NPC cells recruit monocytes and activate TAMs via VEGF and GM-CSF, respectively. A, Cytokine array of the CM of CNE2-EBV, CNE2-EBV+, and CSE-treated CNE2-EBV+ NPC cells. A table summarizing the relative signal intensity of the indicated cytokines is presented on the right. B, The levels of GM-CSF, IL1α, GROα, and VEGF in the CM of CNE2-EBV, CSE-treated CNE2-EBV, CNE2-EBV+, and CSE-treated CNE2-EBV+ NPC cells were measured by ELISA (mean ± SEM; n = 5; **, P < 0.01 and ***, P < 0.001 by two-tailed Student t test). C, Monocyte chemotaxis in response to the indicated recombinant human cytokines at the indicated concentrations (ng/mL; mean ± SEM; n = 3 independent experiments of macrophages from three different donors; **, P < 0.01; ***, P < 0.001 compared with control medium without cytokine treatment by two-tailed Student t test). D, Cytokine levels in the media of macrophages treated with the indicated recombinant human cytokines at the indicated concentrations (ng/mL; mean ± SEM; n = 3 independent experiments of macrophages from three different donors; **, P < 0.01; ***, P < 0.001 compared with control medium without cytokine treatment by two-tailed Student t test). E, Expression of CD206/HLA-DR in macrophages treated with 30% of the indicated tumor-CM in the presence of GM-CSF and/or VEGF neutralizing antibodies. The histograms are representative of four independent experiments of macrophages from four different donors. The values of mean fluorescence intensity (×103, mean ± SEM, n = 5) are indicated on the top-right corners of the representative plots. F, Cytokine levels in the media of macrophages treated as described in E (mean ± SEM; n = 4 independent experiments; N.S., no statistical significance; **, P < 0.01; ***, P < 0.001 by two-tailed Student t test). G, Monocyte chemotaxis in response to the indicated tumor-CM in the presence of GM-CSF and/or VEGF neutralizing antibodies (mean ± SEM; n = 4; NS, no statistical significance; **, P < 0.01; ***, P < 0.001; two-tailed Student t test).

Figure 3.

EBV-replicating NPC cells recruit monocytes and activate TAMs via VEGF and GM-CSF, respectively. A, Cytokine array of the CM of CNE2-EBV, CNE2-EBV+, and CSE-treated CNE2-EBV+ NPC cells. A table summarizing the relative signal intensity of the indicated cytokines is presented on the right. B, The levels of GM-CSF, IL1α, GROα, and VEGF in the CM of CNE2-EBV, CSE-treated CNE2-EBV, CNE2-EBV+, and CSE-treated CNE2-EBV+ NPC cells were measured by ELISA (mean ± SEM; n = 5; **, P < 0.01 and ***, P < 0.001 by two-tailed Student t test). C, Monocyte chemotaxis in response to the indicated recombinant human cytokines at the indicated concentrations (ng/mL; mean ± SEM; n = 3 independent experiments of macrophages from three different donors; **, P < 0.01; ***, P < 0.001 compared with control medium without cytokine treatment by two-tailed Student t test). D, Cytokine levels in the media of macrophages treated with the indicated recombinant human cytokines at the indicated concentrations (ng/mL; mean ± SEM; n = 3 independent experiments of macrophages from three different donors; **, P < 0.01; ***, P < 0.001 compared with control medium without cytokine treatment by two-tailed Student t test). E, Expression of CD206/HLA-DR in macrophages treated with 30% of the indicated tumor-CM in the presence of GM-CSF and/or VEGF neutralizing antibodies. The histograms are representative of four independent experiments of macrophages from four different donors. The values of mean fluorescence intensity (×103, mean ± SEM, n = 5) are indicated on the top-right corners of the representative plots. F, Cytokine levels in the media of macrophages treated as described in E (mean ± SEM; n = 4 independent experiments; N.S., no statistical significance; **, P < 0.01; ***, P < 0.001 by two-tailed Student t test). G, Monocyte chemotaxis in response to the indicated tumor-CM in the presence of GM-CSF and/or VEGF neutralizing antibodies (mean ± SEM; n = 4; NS, no statistical significance; **, P < 0.01; ***, P < 0.001; two-tailed Student t test).

Close modal

NF-κB mediates cytokine production in NPC cells with active EBV replications

Recent studies show that EBV replication activates NF-κB in NPC cells via the interaction between latent membrane protein 1 and tumor necrosis factor receptor-associated factors (28). In addition, both VEGF and GM-CSF have been identified as target genes of NF-κB (20, 29). Thus, we examined NF-κB activation in NPC cells with different EBV replication. We found that EBV+ NPC cells exhibited enhanced NF-κB activation compared with EBV NPC cells, including IκB kinase (IKK) phosphorylation, IκB phosphorylation and degradation (Fig. 4A; Supplementary Fig. S4A), p65 nuclear translocation (Fig. 4C; Supplementary Fig. S4B); and NF-κB binding to its DNA consensus sequence (Fig. 4B). CSE treatment, which increased EBV replication, further promoted NF-κB activation in EBV+ cells. In contrast, CSE treatment could not activate the NF-κB pathway in EBV NPC cells (Fig. 4A–C; Supplementary Fig. S4A and S4B). To determine whether NF-κB activation in responsible for cytokine production in EBV+ cells, we suppressed NF-κB using the pharmacologic inhibitors 4-methyl-N1-(3-phenyl-propyl)-benzene-1,2-diamine (JSH-23; an inhibitor of NF-κB nuclear translocation) or BAY-117082 (an IKK inhibitor) or two small interfering RNAs (siRNAs) targeting p65. We found that NF-κB inhibition significantly decreased GM-CSF and VEGF production from CSE-treated EBV+ NPC cells (Fig. 4D; Supplementary Fig. S4C). These above data indicate that NF-κB is essential for cytokine production in NPC cells with active EBV replications.

Figure 4.

NF-κB mediates cytokine production in NPC cells with active EBV replications. A, Western blot results showing the total and phosphorylated IKK and IκBα levels in untreated or CSE-treated CNE2-EBV and CNE2-EBV+ NPC cells (n = 3 independent experiments). B, NF-κB activities of untreated or CSE-treated CNE2-EBV and CNE2-EBV+ NPC cells, as determined using an electrophoretic mobility shift assay (n = 3 independent experiments). C, Confocal images showing p65 translocation in NPC cells treated as described in A. P65 antibody isotype and quantification of p65 are shown on the right side (mean ± SEM; n = 4; ***, P < 0.001 by two-tailed Student t test). D, Cytokine levels in the media of CSE-treated CNE2-EBV+cells pretreated with dimethyl sulfoxide (DMSO), BAY-117082, or JSH-23 or pretransfected with GFP siRNA or p65 siRNAs (mean ± SEM; n = 4; **, P < 0.01 and ***, P < 0.001 by two-tailed Student t test).

Figure 4.

NF-κB mediates cytokine production in NPC cells with active EBV replications. A, Western blot results showing the total and phosphorylated IKK and IκBα levels in untreated or CSE-treated CNE2-EBV and CNE2-EBV+ NPC cells (n = 3 independent experiments). B, NF-κB activities of untreated or CSE-treated CNE2-EBV and CNE2-EBV+ NPC cells, as determined using an electrophoretic mobility shift assay (n = 3 independent experiments). C, Confocal images showing p65 translocation in NPC cells treated as described in A. P65 antibody isotype and quantification of p65 are shown on the right side (mean ± SEM; n = 4; ***, P < 0.001 by two-tailed Student t test). D, Cytokine levels in the media of CSE-treated CNE2-EBV+cells pretreated with dimethyl sulfoxide (DMSO), BAY-117082, or JSH-23 or pretransfected with GFP siRNA or p65 siRNAs (mean ± SEM; n = 4; **, P < 0.01 and ***, P < 0.001 by two-tailed Student t test).

Close modal

TAMs induce EMT and enhance the invasiveness of NPC cells through CCL18

Previous studies show that TAMs promote metastasis of various types of tumors by inducing EMT (17). Thus, we examined whether activated macrophages can induce EMT and increase the motility of NPC cells. We found that CNE2 cells cocultured with macrophages activated by GM-CSF changed from a rounded to an elongated morphology, decreased E-cadherin expression, increased vimentin expression and EMT-transcription factors (snail, twist, and ZEB1; refs. 30–34), indicating that NPC cells underwent EMT (Fig. 5A–C). As EMT is an essential step of metastasis, we tested the motility of the NPC cells in the absence or presence of macrophages in coculture. Consistently, the abilities of NPC cells to migrate and invade were significantly enhanced when cocultured with activated macrophages (Fig. 5D). Furthermore, compared with the negative control IgG antibody, the addition of neutralizing CCL18 antibody to the coculture system suppressed EMT and the motility induced by GM-CSF-activated macrophages. Moreover, CCL18 alone can induce EMT and enhance the migration and invasiveness of NPC cells to a comparable level as activated macrophages (Fig. 5A–E). Similar results were observed in other two NPC cell lines (Supplementary Fig. S5A–S5C). To investigate whether TAM-derived CCL18 induces EMT in clinical samples, we investigated the correlation of CCL18+ TAM count and ZEB1 expression of NPC cells in human tumor samples (Fig. 5E; Supplementary Fig. S5D). The numbers of nuclear ZEB1+ tumor cells, determined by double immunostaining for ZEB1 and cytokeratin increased in the high CCL18+ TAM group compared with the low one (Fig. 5F), and were positively correlated to CCL18+ TAM count (Pearson R2 = 0.720, P < 0.001, n = 580, Fig. 5G). The above data indicated that CCL18 derived from TAMs promotes NPC EMT and aggressiveness.

Figure 5.

TAMs induce EMT and enhance the invasiveness of NPC cells through CCL18. A, Morphology, E-cadherin, and vimentin immunofluorescent staining of CNE2 cells alone, cocultured with either control macrophages (MΦ) or GM-CSF–activated macrophages in the presence or absence of control IgG or anti-CCL18–neutralizing antibody, or treated with recombinant human CCL18. Scale bar, 20 μm. B, Fluorescence intensity quantification of E-cadherin and vimentin in CNE2 cells treated as described in A. Fluorescence intensity was normalized to control levels (mean ± SEM; n = 3; **, P < 0.01; ***, P < 0.001 by two-tailed Student t test). C, The levels of EMT markers (E-cadherin and vimentin) and EMT transcription factors (Snail, Twist, and ZEB1) in CNE2 cells treated as described in A, as detected using Western blotting (n = 3 independent experiments). D, Migration and invasion assays of CNE2 cells treated as described in A. Scale bar, 50 μm (mean ± SEM; n = 3; **, P < 0.01 by two-tailed Student t test). E, The representative images of CCL18 and CD68 double staining or ZEB1 and cytokeratin (CK) in the human NPC samples with low or high CCL18+ TAM count. Scale bar, 20 μm. F, The number of nuclear ZEB1+ tumor cells significantly increased in the high CCL18+ TAM group compared with the low one (entire cohorts: n = 270 in the high CCL18+ TAM group; n = 310 in the low CCL18+ TAM group; ***, P < 0.001 by two-tailed Student t test). G, The numbers of nuclear ZEB1+ tumor cells were significantly correlated to the CCL18+ TAM counts (n = 580, Pearson R and P < 0.001).

Figure 5.

TAMs induce EMT and enhance the invasiveness of NPC cells through CCL18. A, Morphology, E-cadherin, and vimentin immunofluorescent staining of CNE2 cells alone, cocultured with either control macrophages (MΦ) or GM-CSF–activated macrophages in the presence or absence of control IgG or anti-CCL18–neutralizing antibody, or treated with recombinant human CCL18. Scale bar, 20 μm. B, Fluorescence intensity quantification of E-cadherin and vimentin in CNE2 cells treated as described in A. Fluorescence intensity was normalized to control levels (mean ± SEM; n = 3; **, P < 0.01; ***, P < 0.001 by two-tailed Student t test). C, The levels of EMT markers (E-cadherin and vimentin) and EMT transcription factors (Snail, Twist, and ZEB1) in CNE2 cells treated as described in A, as detected using Western blotting (n = 3 independent experiments). D, Migration and invasion assays of CNE2 cells treated as described in A. Scale bar, 50 μm (mean ± SEM; n = 3; **, P < 0.01 by two-tailed Student t test). E, The representative images of CCL18 and CD68 double staining or ZEB1 and cytokeratin (CK) in the human NPC samples with low or high CCL18+ TAM count. Scale bar, 20 μm. F, The number of nuclear ZEB1+ tumor cells significantly increased in the high CCL18+ TAM group compared with the low one (entire cohorts: n = 270 in the high CCL18+ TAM group; n = 310 in the low CCL18+ TAM group; ***, P < 0.001 by two-tailed Student t test). G, The numbers of nuclear ZEB1+ tumor cells were significantly correlated to the CCL18+ TAM counts (n = 580, Pearson R and P < 0.001).

Close modal

The interaction between EBV+NPC cells and macrophages promotes metastasis in a humanized mouse model

CCL18, the key cytokine released by human TAMs, does not have a mouse homolog. Therefore, to mimic the human tumor microenvironment, we evaluated the role of the interaction between EBV+ NPC cells and macrophages in a humanized mouse model as previously described (20, 35). CNE2-EBV cells, CNE2-EBV+ or CSE-treated CNE2-EBV+ cells were injected into the subcutaneous tissue of irradiated NOD/SCID mice without HSC transplantation and humanized mice, respectively. Consistent with the findings in vitro, injection of CNE-EBV+ cells into the subcutaneous tissue of humanized mice led to more CCL18+ macrophage infiltration in the xenografts, increased GM-CSF and VEGF production and exhibited EMT phenotype (decrease of E-cadherin, increase of vimentin, and ZEB1) compared with injection of EBV cells into the humanized mice. In addition, injection of CSE-treated CNE2-EBV+ cells further promoted the above effects. In sharp contrast with humanized mice, the xenografts of nonhumanized mice had no CCL18+ macrophage infiltration, low GM-CSF and VEGF production and absence of EMT phenotype regardless of EBV+ or not (Fig. 6A; Supplementary Fig. S6A–S6F). In addition, EBV replication in NPC cells led to increased metastasis in the lungs of the humanized mice, but not in those of nonhumanized mice, as determined by counting of metastatic colonies in the lungs (Fig. 6B) and quantification of human hypoxanthine phosphor-ribosyl transferase (HPRT) mRNA expression (Fig. 6C). Similar results were observed in the other NPC cell line HNE1 (Fig. 6B and C). To study the potential therapeutic values of targeting inflammatory cytokines in NPC, we treated the humanized mice bearing CSE-treated EBV+NPC xenografts with anti-GM-CSF, anti-CCL18, and anti-VEGF–neutralizing antibody. We found that the addition of either antibody could reduce lung metastases, while the combination use of these three neutralizing antibodies had the strongest effect on metastases inhibition (Fig. 6D and E; Supplementary Fig. S6G). The above data suggest that GM-CSF, VEGF, and CCL18 could represent therapeutic targets for inhibiting EBV-infected NPC metastasis.

Figure 6.

The interaction between EBV+NPC cells and macrophages promotes metastasis in a humanized mouse model. A, Representative images of IHC staining of E-cadherin, vimentin, GM-CSF, CCL18, VEGF, and ZEB1 in sections from the indicated tumor xenografts. The isotype antibodies for indicated antibodies are shown as negative controls. Scale bar, 100 μm. B, The metastasis in lungs of humanized mice or NOD/SCID mice bearing different tumors was determined by quantification of metastatic tumor nodules (n = 9 for each group; *, P < 0.05; ***, P < 0.001 by two-tailed Student t test). C, The metastasis in lungs of humanized mice or NOD/SCID mice bearing different tumors was determined by expression of human HPRT mRNA relative to mouse 18S rRNA (n = 9 for each group; **, P < 0.01; ***, P < 0.001 by two-tailed student t test). D, Representative hematoxylin and eosin (H&E)-stained and human cytokeratin (CK)-immunostained lung sections from humanized mice bearing EBV or CSE-treated EBV+ xenografts injected with control IgG, anti-CCL18, anti-GM-CSF, and/or anti-VEGF–neutralizing antibody via the tail vein. Scale bars, 100 μm. E, The metastasis in lungs of humanized mice as described in D was determined by quantification of metastatic tumor nodules (top) and expression of human HPRT mRNA relative to mouse 18S rRNA (bottom; n = 9 for each group; **, P < 0.01; ***, P < 0.001 by two-tailed Student t test).

Figure 6.

The interaction between EBV+NPC cells and macrophages promotes metastasis in a humanized mouse model. A, Representative images of IHC staining of E-cadherin, vimentin, GM-CSF, CCL18, VEGF, and ZEB1 in sections from the indicated tumor xenografts. The isotype antibodies for indicated antibodies are shown as negative controls. Scale bar, 100 μm. B, The metastasis in lungs of humanized mice or NOD/SCID mice bearing different tumors was determined by quantification of metastatic tumor nodules (n = 9 for each group; *, P < 0.05; ***, P < 0.001 by two-tailed Student t test). C, The metastasis in lungs of humanized mice or NOD/SCID mice bearing different tumors was determined by expression of human HPRT mRNA relative to mouse 18S rRNA (n = 9 for each group; **, P < 0.01; ***, P < 0.001 by two-tailed student t test). D, Representative hematoxylin and eosin (H&E)-stained and human cytokeratin (CK)-immunostained lung sections from humanized mice bearing EBV or CSE-treated EBV+ xenografts injected with control IgG, anti-CCL18, anti-GM-CSF, and/or anti-VEGF–neutralizing antibody via the tail vein. Scale bars, 100 μm. E, The metastasis in lungs of humanized mice as described in D was determined by quantification of metastatic tumor nodules (top) and expression of human HPRT mRNA relative to mouse 18S rRNA (bottom; n = 9 for each group; **, P < 0.01; ***, P < 0.001 by two-tailed Student t test).

Close modal

More than 20% of human cancers are attributed to chronic microbial infection. However, how chronic inflammation promotes tumor progression remains incompletely understood. Our data indicate that NPC with different EBV replications has distinctive capacity to foster a tumor-promoting microenvironment. The NPC cells with active EBV replications recruit monocytes by VEGF and activate monocyte to TAMs by GM-CSF in an NF-κB-dependent manner. Reciprocally, TAMs promote tumor metastasis and further NF-κB activation by CCL18. Collectively, our study demonstrated a feedforward loop between tumor cells and macrophages, suggesting that the metastatic potential may evolve hand in hand with virus-induced chronic inflammation.

Pathogen-induced chronic inflammation plays a central role in various types of tumors. Persistent microbial infections elicit inflammatory cytokine storm. This overactive immune response causes tissue damage and DNA mutation, which subsequently promote tumor progression in various malignancies, such as hepatocellular, cervical, and gastric carcinoma (36, 37). Similarly, it has been reported TAM infiltration is tightly associated with poor prognosis in NPC with EBV infection (38), suggesting the interaction between EBV activity and macrophages may have a role in NPC development. However, despite it is well-documented that EBV is associated with NPC, how EBV-induced inflammation promotes NPC progression remains largely unknown. Our data showed the EBV replication in NPC cells can induce an inflammatory microenvironment by recruiting and activating macrophages. This is consistent with the previous reports on other tumor types. However, the mechanism here seems to be different—through the NF-κB pathway directly activated by EBV replications in tumor cells. Indeed, a number of oncoproteins encoded by EBV can activate NF-κB. For example, Latent membrane protein 1 (LMP1), a major EBV oncoprotein, interacts with TRAFs, and activates both noncanonical and canonical NF-κB through its functional domains termed CTAR1/2/3 (39, 40). Another viral protein termed EBNA1, which enables the replication and segregation of EBV genomes in proliferating infected cells, can increase cell migration via recruitment of the transcription factors Ap1, Sp1, and NF-κB to upregulate matrix metalloproteinase-9 (MMP-9; ref. 41). It is possible that upregulation of viral oncoproteins like LMP1 and EBNA by EBV replications induces NF-κB activation and sequence cytokine production.

GM-CSF and VEGF have individually been identified to be involved in TAM formation (20, 42). Moreover, GM-CSF has also been proved to play a role in the formation of tumor immunosuppressive microenvironment (43). However, GM-CSF and VEGF have not previously been considered for dual targeting. Emerging evidence has validated the concept that inhibition of key signaling pathways of TAM could elicit potent antitumor activities in preclinical tumor models and cancer patients (16). For example, the RG7155 antibody, a humanized anti-human CSF1 receptor mAb, specifically depletes macrophages and exhibits antitumor efficacy in patients of various malignancies (44, 45). However, the presence of tumor-derived GM-CSF–induced RG7155 resistance of TAMs, indicating the therapeutic potential for targeting GM-CSF (44). GM-CSF neutralization is a rapidly expanding market with multiple agents in phase I–II development. These agents are well tolerated and demonstrate efficacy in various chronic inflammatory diseases, including rheumatic arthritis (46). However, the role of GM-CSF neutralization in virus-induced tumor progression has not been reported. In addition, increased VEGF-A expression was observed in EBV+ NPC tumors and related to lymph node metastasis and poor long-term survival (47), which is in consistent with our data. Bevacizumab, VEGF-neutralizing antibody, has already become the standard treatment for various cancer types (48, 49). Our study demonstrated that neutralization of GM-CSF and VEGF significantly reduced NPC metastasis in vivo, highlighting the rationales to investigate the therapeutic values of GM-CSF and VEGF-neutralizing antibodies in NPC patients in future clinical trials.

In addition, EBV is also etiologically linked to a wide range of other human tumors, including several B-cell malignancies, extranodal lymphomas of T/NK-cell origin and a distinct subset of gastric carcinomas (50). Our study demonstrated a positive feed-forward loop between EBV-infected tumor cells and macrophages play a crucial role in NPC progression. Whether this vicious cycle in tumor microenvironment also mediates the progression of other malignancies related to EBV infection warrants further investigation.

In summary, our study suggests an important role for the interaction between NPC cells and any step of this interaction, including GM-CSF, CCL18, and VEGF, could be a potential therapeutic target for NPC treatment.

No potential conflicts of interest were disclosed.

Conception and design: S.-J. Song, S.-C. Su

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): D. Huang, S.-J. Song, Z.-Z. Wu, W. Wu, J.-N. Chen

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): D. Huang, S.-J. Song, M.-S. Zeng

Writing, review, and/or revision of the manuscript: D. Huang, S.-J. Song

Study supervision: X.-Y. Cui

This work was supported by grants from the Natural Science Foundation of China (81472468, 81672614, 81622036, 81490750, 81372816, 81621004, 81230060, and 81442009), the National Key Research and Development Program of China (2016YFC1302300), Science Foundation of Guangdong Province (2016A030306023, 2014A0303130940, S2012030006287, and 2016B030229004), China Postdoctoral Science Foundation (2014M550449), Guangzhou Science and Technology Project (201710010083), Translational medicine public platform of Guangdong Province (4202037), Tip-top Scientific and Technical Innovative Youth Talents of Guangdong special support program (No. 2016TQ03R553), Guangdong Department of Science & Technology Translational Medicine Center grant (2011A080300002) and Sun Yat-sen University (16ykpy07).

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Parkin
DM
. 
The global health burden of infection-associated cancers in the year 2002
.
Int J Cancer
2006
;
118
:
3030
44
.
2.
Cohen
JI
. 
Epstein–Barr virus infection
.
N Engl J Med
2000
;
343
:
481
92
.
3.
zur Hausen
H
,
O'Neill
FJ
,
Freese
UK
,
Hecker
E
. 
Persisting oncogenic herpesvirus induced by the tumour promotor TPA
.
Nature
1978
;
272
:
373
5
.
4.
Fang
CY
,
Lee
CH
,
Wu
CC
,
Chang
YT
,
Yu
SL
,
Chou
SP
, et al
Recurrent chemical reactivations of EBV promotes genome instability and enhances tumor progression of nasopharyngeal carcinoma cells
.
Int J Cancer
2009
;
124
:
2016
25
.
5.
Xu
FH
,
Xiong
D
,
Xu
YF
,
Cao
SM
,
Xue
WQ
,
Qin
HD
, et al
An epidemiological and molecular study of the relationship between smoking, risk of nasopharyngeal carcinoma, and Epstein–Barr virus activation
.
J Natl Cancer Inst
2012
;
104
:
1396
410
.
6.
Qiu
J
,
Smith
P
,
Leahy
L
,
Thorley-Lawson
DA
. 
The Epstein–Barr virus encoded BART miRNAs potentiate tumor growth in vivo
.
PLoS Pathog
2015
;
11
:
e1004561
.
7.
Kong
QL
,
Hu
LJ
,
Cao
JY
,
Huang
YJ
,
Xu
LH
,
Liang
Y
, et al
Epstein–Barr virus-encoded LMP2A induces an epithelial-mesenchymal transition and increases the number of side population stem-like cancer cells in nasopharyngeal carcinoma
.
PLoS Pathog
2010
;
6
:
e1000940
.
8.
Banerjee
AS
,
Pal
AD
,
Banerjee
S
. 
Epstein Barr virus-encoded small non-coding RNAs induce cancer cell chemoresistance and migration
.
Virology
2013
;
443
:
294
305
.
9.
Kuraishy
A
,
Karin
M
,
Grivennikov
SI
. 
Tumor promotion via injury- and death-induced inflammation
.
Immunity
2011
;
35
:
467
77
.
10.
Moore
PS
,
Chang
Y
. 
Why do viruses cause cancer? Highlights of the first century of human tumour virology
.
Nat Rev Cancer
2010
;
10
:
878
89
.
11.
Grivennikov
SI
,
Greten
FR
,
Karin
M
. 
Immunity, inflammation, and cancer
.
Cell
2010
;
140
:
883
99
.
12.
Tran
le S
,
Mittal
D
,
Mattarollo
SR
,
Frazer
IH
. 
Interleukin-17A promotes arginase-1 production and 2,4-dinitrochlorobenzene-induced acute hyperinflammation in human papillomavirus E7 oncoprotein-expressing skin
.
J Innate Immun
2015
;
7
:
392
404
.
13.
Condeelis
J
,
Pollard
JW
. 
Macrophages: obligate partners for tumor cell migration, invasion, and metastasis
.
Cell
2006
;
124
:
263
6
.
14.
Pollard
JW
. 
Tumour-educated macrophages promote tumour progression and metastasis
.
Nat Rev Cancer
2004
;
4
:
71
8
.
15.
Allavena
P
,
Sica
A
,
Solinas
G
,
Porta
C
,
Mantovani
A
. 
The inflammatory micro-environment in tumor progression: the role of tumor-associated macrophages
.
Crit Rev Oncol Hematol
2008
;
66
:
1
9
.
16.
Ruffell
B
,
Coussens
LM
. 
Macrophages and therapeutic resistance in cancer
.
Cancer Cell
2015
;
27
:
462
72
.
17.
Chen
J
,
Yao
Y
,
Gong
C
,
Yu
F
,
Su
S
,
Liu
B
, et al
CCL18 from tumor-associated macrophages promotes breast cancer metastasis via PITPNM3
.
Cancer Cell
2011
;
19
:
541
55
.
18.
Schutyser
E
,
Struyf
S
,
Proost
P
,
Opdenakker
G
,
Laureys
G
,
Verhasselt
B
, et al
Identification of biologically active chemokine isoforms from ascitic fluid and elevated levels of CCL18/pulmonary and activation-regulated chemokine in ovarian carcinoma
.
J Biol Chem
2002
;
277
:
24584
93
.
19.
Hou
X
,
Zhang
Y
,
Qiao
H
. 
CCL18 promotes the invasion and migration of gastric cancer cells via ERK1/2/NF-kappaB signaling pathway
.
Tumour Biol
2016
;
37
:
641
51
.
20.
Su
S
,
Liu
Q
,
Chen
J
,
Chen
F
,
He
C
,
Huang
D
, et al
A positive feedback loop between mesenchymal-like cancer cells and macrophages is essential to breast cancer metastasis
.
Cancer Cell
2014
;
25
:
605
20
.
21.
Lo
AKF
,
Lo
KW
,
Tsao
SW
,
Wong
HL
,
Hui
JWY
,
To
KF
, et al
Epstein–Barr virus infection alters cellular signal cascades in human nasopharyngeal epithelial cells
.
Neoplasia
2006
;
8
:
173
80
.
22.
Zhao
QY
,
Wang
PP
,
Huang
ZL
,
Peng
L
,
Lin
CS
,
Gao
ZL
, et al
Tumoral indoleamine 2, 3-dioxygenase 1 is regulated by monocytes and T lymphocytes collaboration in hepatocellular carcinoma
.
Oncotarget
2016
;
7
:
14781
90
.
23.
Huang
D
,
Su
S
,
Cui
X
,
Shen
X
,
Zeng
Y
,
Wu
W
, et al
Nerve fibers in breast cancer tissues indicate aggressive tumor progression
.
Medicine
2014
;
93
:
e172
.
24.
Pettersen
JS
,
Fuentes-Duculan
J
,
Suarez-Farinas
M
,
Pierson
KC
,
Pitts-Kiefer
A
,
Fan
L
, et al
Tumor-associated macrophages in the cutaneous SCC microenvironment are heterogeneously activated
.
J Invest Dermatol
2011
;
131
:
1322
30
.
25.
Galon
J
,
Costes
A
,
Sanchez-Cabo
F
,
Kirilovsky
A
,
Mlecnik
B
,
Lagorce-Pages
C
, et al
Type, density, and location of immune cells within human colorectal tumors predict clinical outcome
.
Science
2006
;
313
:
1960
4
.
26.
Ansell
SM
,
Stenson
M
,
Habermann
TM
,
Jelinek
DF
,
Witzig
TE
. 
Cd4+ T-cell immune response to large B-cell non-Hodgkin's lymphoma predicts patient outcome
.
J Clin Oncol
2001
;
19
:
720
6
.
27.
Su
SC
,
Zhao
QY
,
He
CH
,
Huang
D
,
Liu
J
,
Chen
F
, et al
miR-142–5p and miR-130a-3p are regulated by IL-4 and IL-13 and control profibrogenic macrophage program
.
Nat Commun
2015
;
6
:
8523. doi: 10.1038/ncomms9523
.
28.
Shair
KH
,
Bendt
KM
,
Edwards
RH
,
Bedford
EC
,
Nielsen
JN
,
Raab-Traub
N
. 
EBV latent membrane protein 1 activates Akt, NFkappaB, and Stat3 in B cell lymphomas
.
PLoS Pathog
2007
;
3
:
e166
.
29.
Ko
HM
,
Seo
KH
,
Han
SJ
,
Ahn
KY
,
Choi
IH
,
Koh
GY
, et al
Nuclear factor kappaB dependency of platelet-activating factor-induced angiogenesis
.
Cancer Res
2002
;
62
:
1809
14
.
30.
Sosic
D
,
Richardson
JA
,
Yu
K
,
Ornitz
DM
,
Olson
EN
. 
Twist regulates cytokine gene expression through a negative feedback loop that represses NF-kappa B activity
.
Cell
2003
;
112
:
169
80
.
31.
Peinado
H
,
Olmeda
D
,
Cano
A
. 
Snail, Zeb and bHLH factors in tumour progression: an alliance against the epithelial phenotype?
Nat Rev Cancer
2007
;
7
:
415
28
.
32.
Gregory
PA
,
Bert
AG
,
Paterson
EL
,
Barry
SC
,
Tsykin
A
,
Farshid
G
, et al
The miR-200 family and miR-205 regulate epithelial to mesenchymal transition by targeting ZEB1 and SIP1
.
Nat Cell Biol
2008
;
10
:
593
601
.
33.
Wellner
U
,
Schubert
J
,
Burk
UC
,
Schmalhofer
O
,
Zhu
F
,
Sonntag
A
, et al
The EMT-activator ZEB1 promotes tumorigenicity by repressing stemness-inhibiting microRNAs
.
Nat Cell Biol
2009
;
11
:
1487
95
.
34.
Kim
T
,
Veronese
A
,
Pichiorri
F
,
Lee
TJ
,
Jeon
YJ
,
Volinia
S
, et al
p53 regulates epithelial-mesenchymal transition through microRNAs targeting ZEB1 and ZEB2
.
J Exp Med
2011
;
208
:
875
83
.
35.
Zitvogel
L
,
Pitt
JM
,
Daillere
R
,
Smyth
MJ
,
Kroemer
G
. 
Mouse models in oncoimmunology
.
Nat Rev Cancer
2016
;
16
:
759
73
.
36.
Wang
Y
,
Kato
N
,
Hoshida
Y
,
Yoshida
H
,
Taniguchi
H
,
Goto
T
, et al
Interleukin-1beta gene polymorphisms associated with hepatocellular carcinoma in hepatitis C virus infection
.
Hepatology
2003
;
37
:
65
71
.
37.
Guidotti
LG
,
Ishikawa
T
,
Hobbs
MV
,
Matzke
B
,
Schreiber
R
,
Chisari
FV
. 
Intracellular inactivation of the hepatitis B virus by cytotoxic T lymphocytes
.
Immunity
1996
;
4
:
25
36
.
38.
Liao
Q
,
Zeng
Z
,
Guo
X
,
Li
X
,
Wei
F
,
Zhang
W
, et al
LPLUNC1 suppresses IL-6-induced nasopharyngeal carcinoma cell proliferation via inhibiting the Stat3 activation
.
Oncogene
2014
;
33
:
2098
109
.
39.
Eliopoulos
AG
,
Caamano
JH
,
Flavell
J
,
Reynolds
GM
,
Murray
PG
,
Poyet
JL
, et al
Epstein–Barr virus-encoded latent infection membrane protein 1 regulates the processing of p100 NF-kappaB2 to p52 via an IKKgamma/NEMO-independent signalling pathway
.
Oncogene
2003
;
22
:
7557
69
.
40.
Luftig
M
,
Yasui
T
,
Soni
V
,
Kang
MS
,
Jacobson
N
,
Cahir-McFarland
E
, et al
Epstein–Barr virus latent infection membrane protein 1 TRAF-binding site induces NIK/IKK alpha-dependent noncanonical NF-kappaB activation
.
Proc Natl Acad Sci U S A
2004
;
101
:
141
6
.
41.
Ning
S
,
Hahn
AM
,
Huye
LE
,
Pagano
JS
. 
Interferon regulatory factor 7 regulates expression of Epstein–Barr virus latent membrane protein 1: a regulatory circuit
.
J Virol
2003
;
77
:
9359
68
.
42.
Qian
BZ
,
Zhang
H
,
Li
JF
,
He
TF
,
Yeo
EJ
,
Soong
DYH
, et al
FLT1 signaling in metastasis-associated macrophages activates an inflammatory signature that promotes breast cancer metastasis
.
J Exp Med
2015
;
212
:
1433
48
.
43.
Gallina
G
,
Dolcetti
L
,
Serafini
P
,
De Santo
C
,
Marigo
I
,
Colombo
MP
, et al
Tumors induce a subset of inflammatory monocytes with immunosuppressive activity on CD8+ T cells
.
J Clin Invest
2006
;
116
:
2777
90
.
44.
Ries
CH
,
Cannarile
MA
,
Hoves
S
,
Benz
J
,
Wartha
K
,
Runza
V
, et al
Targeting tumor-associated macrophages with anti-CSF-1R antibody reveals a strategy for cancer therapy
.
Cancer Cell
2014
;
25
:
846
59
.
45.
Cassier
PA
,
Italiano
A
,
Gomez-Roca
CA
,
Le Tourneau
C
,
Toulmonde
M
,
Cannarile
MA
, et al
CSF1R inhibition with emactuzumab in locally advanced diffuse-type tenosynovial giant cell tumours of the soft tissue: a dose-escalation and dose-expansion phase 1 study
.
Lancet Oncol
2015
;
16
:
949
56
.
46.
Wicks
IP
,
Roberts
AW
. 
Targeting GM-CSF in inflammatory diseases
.
Nat Rev Rheumatol
2016
;
12
:
37
48
.
47.
Krishna
SM
,
James
S
,
Balaram
P
. 
Expression of VEGF as prognosticator in primary nasopharyngeal cancer and its relation to EBV status
.
Virus Res
2006
;
115
:
85
90
.
48.
Yang
JC
,
Haworth
L
,
Sherry
RM
,
Hwu
P
,
Schwartzentruber
DJ
,
Topalian
SL
, et al
A randomized trial of bevacizumab, an anti-vascular endothelial growth factor antibody, for metastatic renal cancer
.
N Engl J Med
2003
;
349
:
427
34
.
49.
Matsuoka
A
,
Maeda
O
,
Ando
Y
. 
FOLFOXIRI and bevacizumab for metastatic colorectal cancer
.
N Engl J Med
2015
;
372
:
291
.
50.
Taylor
GS
,
Long
HM
,
Brooks
JM
,
Rickinson
AB
,
Hislop
AD
. 
The immunology of Epstein–Barr virus-induced disease
.
Annu Rev Immunol
2015
;
33
:
787
821
.

Supplementary data