Abstract
Background:
Radiofrequency ablation (RFA) is the primary curative treatment for hepatocellular carcinoma (HCC) patients who are not eligible for surgery. However, the effects of RFA on the global tumor immune response remain unclear.
Method:
In this study, we examined the phenotypic and functional changes in peripheral blood mononuclear cells (PBMCs) from recurrent HCC patients who had undergone two RFA treatments using mass cytometry and high-throughput mRNA assays.
Results:
We observed significant increase in monocytes and decrease in T cell subpopulations three days after the first RFA treatment and three days after the second RFA treatment. The down-regulation of GZMB, GZMH, GZMK, and CD8A, which are involved in the cytotoxic function of T cells, was observed following RFA. Furthermore, the population of CD8 effector and memory T cells (CD8 Teff and CD8 Tem) significantly decreased after RFA. The expression of CD5 and CD161 in various T cell subpopulations also showed significant reductions. Additionally, elevated secretion of VEGF was observed in monocytes, B cells, regulatory T cells (Tregs), and CD4 naive T cells.
Conclusion:
In recurrent HCC patients, serum components derived from radiofrequency therapy can enhance the antigen-presenting capacity of monocytes. However, they also inhibit the anti-cancer immune response by reducing the population of CD8 effector and memory T cells and suppressing the activation of T cells, as well as down-regulating the expression of CD161 and CD5 in various T cell subpopulations. These tumor-derived components also contribute to an immunosuppressive microenvironment by promoting the secretion of VEGF in monocytes, Tregs, B cells, and CD4 naive T cells.
Introduction
Hepatocellular carcinoma is one of the most common malignancies, ranking as the third most lethal malignancy worldwide (, ). Following primary treatment, a large proportion of patients suffered cancer recurrence. For HCC patients ineligible for surgical resection or liver transplantation, radiofrequency ablation (RFA) has emerged as a frequently employed curative approach. Other than direct damage to tumor tissue, RFA also has the potential to influence the overall immune response (), which may lead to a durable anti-cancer immune response and improved tumor control. Nevertheless, this alone does not guarantee complete remission. Approximately 13.3% patients receiving RFA experienced local recurrence, while 45% suffered from distant recurrence within 3 years (). Notably, inflammation induced by incomplete radiofrequency ablation accelerates tumor progression and hinders PD-1 immunotherapy. Therefore, the detailed understanding of the effect of RFA on immune cells in patients with recurrent HCC is urgently needed for further clinical trials.
Peripheral blood mononuclear cells (PBMC) are believed to be the origin of the tumor infiltrating immune cells. Given the challenges in tumor sample acquisition in RFA, PBMCs are considered to be the optimal target to analyze the overall immune status following RFA. Notably, PBMC components were reported to change remarkably after RFA and can serve as biomarkers for the effectiveness of RFA treatment (, ). However, these clinical studies could study only a few subtypes of the immune cells because of the limited number of signal channels of flow cytometry. Mass cytometry, also known as cytometry by time-of-flight (CyTOF), uniquely combines metal-labeling of antibodies with mass spectrometry to enable high-dimensional measurement of the characteristics of individual cells (). This technology enables us to analyze immune cell subsets and the expression of functional molecules at the single-cell level.
In this study, we used mass cytometry to comprehensively characterize the phenotypic and functional alterations of PBMCs from the patients with recurrent HCC who received RFA treatment, aiming to unveil the dynamic landscape of the anti-cancer immune response following both the primary and secondary RFA. Also, the immune–related mRNA profiles of the PBMC were analyzed using a high throughput RNA assay system, identifying the down-regulated expression of markers for a set of cytotoxic function-related genes after RFA. These findings were further validated by mass cytometry analysis on health donor-derived PBMC treated with serum from HCC patients before and after RFA. This research offers valuable insights for future investigations on the potential combination of immunopotentiating therapy with RFA in the treatment of HCC.
Materials and methods
Patients
In this study, blood samples were obtained from 21 recurrent HCC patients who received primary and secondary RFA therapy in YouAn Hospital from January 2011 to January 2015. The primary eligibility criteria included patients with primary HCC, who would undergo RFA with curative intent. The diagnosis of HCC was preferably histologically confirmed. In cases when tumor biopsy results were unavailable, diagnosis was established by contrast-enhanced CT or MRI. Treatment selection in each patient was decided by the clinical doctors guided by the NCCN and Chinese HCC treatment guidelines. Blood samples were collected before RFA treatment and three days after RFA and stored at. -80°C. Follow-up assessments and imaging examinations were conducted periodically for a minimum of 3 years. Patients who had HCC recurrence diagnosed by imaging examinations and received second RFA were finally enrolled. Patient characteristics and disease classification are shown in Table 1. All patients provided written informed consent and this study was approved by the Medical Ethics Committee of YouAn Hospital.
Table 1
| ID | Age | Gender | Hepatitis | Pathological results | |||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Pathological type | Degree of differentiation; | HBS Ag | HBC Ag | Hepa | CK19 | CD34 | P53 | GPC3 | Ki67 | ||||
| 1 | 48 | Male | HBV | HCC | NA | + | – | + | – | + | – | + | + |
| 2 | 46 | Male | HBV | HCC | Moderate | + | – | + | – | + | – | + | + |
| 3 | 72 | Male | HBV | HCC | Moderate | +++ | – | +++ | – | + | + | +++ | + |
| 4 | 63 | Male | HBV | HCC | Moderate | + | – | + | – | – | + | – | + |
| 5 | 52 | Male | HBV | HCC | Low | + | – | + | – | + | + | + | + |
| 6 | 55 | Male | HBV | HCC | High | + | – | + | – | + | – | – | – |
| 7 | 68 | Male | HBV | HCC | High | + | – | +++ | – | ++ | – | ++ | + |
| 8 | 61 | Male | HBV | HCC | Moderate | ++ | – | +++ | – | – | + | – | + |
| 9 | 67 | Male | No | HCC | Moderate | – | – | +++ | – | – | + | + | ++ |
| 10 | 47 | Male | No | HCC | Moderate | – | – | + | – | – | + | + | + |
| 11 | 56 | Male | No | HCC | Moderate | – | – | + | – | + | +++ | – | ++ |
| 12 | 40 | Male | No | HCC | Moderate | – | – | + | – | + | – | – | + |
| 13 | 68 | Male | HBV | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 14 | 68 | Male | No | HCC | Moderate | – | – | +++ | – | + | ++ | ± | + |
| 15 | 44 | Male | No | HCC | Low | – | – | + | – | – | +++ | +++ | + |
| 16 | 55 | Male | HBV | HCC | High | + | – | + | – | + | – | – | – |
| 17 | 56 | Male | HBV | HCC | Moderate | – | – | + | – | + | ++ | +++ | + |
| 18 | 72 | Male | No | HCC | Moderate-Low | – | – | + | – | + | + | ± | + |
| 19 | 47 | Male | No | HCC | Moderate | – | – | + | ± | + | – | + | + |
| 20 | 53 | Male | No | HCC | High | + | – | +++ | – | ++ | – | ++ | + |
| 21 | 56 | Male | HBV | HCC | High | + | – | + | – | + | – | – | – |
The Baseline Clinical and Pathological characters of the patients.
Patient 13 was diagnosed by CT and ultrasonography examinations, while other patients were diagnosed by pathological examinations. HCC, Hepatocellular carcinom; HBV, Hepatitis B virus; HBSAg, Hepatitis B surface antigen; HBC Ag, Hepatitis C surface antigen; CK19, Cytokeratin 19; CD34, Cluster of differentiation 34; GPC3, glypican; NA, not available.
-: Negative;+: weakly positive ++ :moderate positive +++: strongly positive.
PBMC isolation and storage
Venous blood samples were obtained on the day before the first RFA treatment, as well as 3 days after the first treatment, the day before the second RFA treatment, and 3 days after the second treatment. PBMCs were separated from whole blood using Human Peripheral Blood Lymphocytes Separation Medium (#P8610, Solarbio, Beijing, China) by means of density gradient centrifugation following the instruction of the kit. The isolated PBMCs were stored at -80°C until the mass cytometry or mRNA examination.
In vitro PBMC stimulation
Serum samples was collected from patients who received RFA before and after RFA treatment, inactivated at 56°C for 30mins, and then used for PBMC culture. For removal of the big molecules, the supernatant of serum was transferred to the centrifuge tube, centrifuged at 4500rpm for 5 min, the serum in the centrifuge tube containing big molecules was removed and the serum in the bottom was collected and filtered by 0.22um filter mesh, then the serum was used for PBMC culture. PBMCs from healthy donor were cultured overnight in 6-well plates at a concentration of 1 × 106cells/well, followed by exposure to prepared HCC patient’s serum (preparation was described above). After 48 hours incubation, PBMCs were stimulated with leukocyte activation cocktail (plus protein transport inhibitors, BD Biosciences, USA) for 4 hours before mass cytometry examination.
Mass cytometry
Pre-conjugated antibodies panels (Maxpar human T cell phenotyping panel kit and expansion panel kit) were purchased from Fluidigm Company. Samples were rapidly rewarmed, then washed and stained with cell surface antibodies for 30 minutes on ice. After that, PBMCs were permeabilized at 4°C overnight and stained with intracellular antibodies for 30 minutes on ice. The antibody-labeled PBMCs were washed and incubated in 0.125nM Ir intercalator (Fluidigm) diluted in PBS containing 2% formaldehyde, and stored at 4°C until the mass cytometry examination.Before acquisition, samples were washed with de-ionized water and then resuspended at a concentration of 1×106 cells/ml with deionized water containing a 1/20 dilution of EQ 4 Element Beads. PBMCs were then examined by a CyTOF2 mass cytometry (Fluidigm DVS Sciences Inc., Toronto, ON, Canada).
High-throughput mRNA analysis
Total RNA was extracted from PBMCs using the QIAamp RNA Mini Kit following the manufacturer’s instructions. The random primer was used to synthesis cDNA following the standard protocol of SuperScript III RT kit (Invitrogen). The synthesized cDNA was added to 1.5ml tube added with 190ul MILIQ H2O. The diluted cDNA was added with 200ul saturated phenolic chloroform mixture in the ice for 10min and centrifuged at 12000g for 10min. The upper aqueous phase was collected into a new 1.5ml tube added with 2ug glycogen and 500ul ethanol. After stored at -80 degree refrigerator for 8h, the sample was centrifuged at 14000g for 30min. The pellet was wash with 1ml 70% ethanol and centrifuged at 10000g for 5min. The cDNA was dissolved in 10ul H2O. cDNA from PBMCs was screened for expression of housekeeping genes ACTB and GAPDH and then were analyzed with TaqMan Universal PCR Master Mix (Applied Biosystems, Foster City, CA, USA), 20X GE Sample Loading Reagent (Fluidigm, South San Francisco, CA, USA), 2X Assay Loading Reagent (Fluidigm, South San Francisco, CA, USA) and individual TaqMan Gene Expression Assay using 96.96 Dynamic Array™ IFC (Fluidigm, South San Francisco, CA, USA) on a BioMark HD System (Fluidigm, South San Francisco, CA, USA). Ct (threshold cycle) values were processed and analyzed using the BioMark Real-Time PCR Analysis software (Fluidigm South San Francisco, CA, USA).
Statistical analysis
Mass Cytometry data were firstly normalized using EQ Four Element Calibration Beads (EQ Beads, 201078, Fluidigm) according to manufacturer’s instructions, then the cell debris were removed according the 191Ir and 193Ir channel. Doublets were removed according to the Even Length. Furthermore, CD45 gate were used to isolate all the blood leukocyte, then a clustering panel including CD3, CD4, CD8A, CCR4, CCR5, CCR7, HLA-DR, CD11A, CD16, CD25, CD44, CD45RA, CD45RO and CD127 were used to manually differentiate subpopulations of the blood leukocytes. Meanwhile, blood leukocytes were also clustered into 200 nodes using the SPADE analysis with a 20% down-sampling rate. Then repeated measurements ANOVA were performed to analyze the difference of cell count and functional marker expression of the subpopulations between the four time points. Paired t test were used to analyze the difference of mRNA expression, cell count and protein expression of the subpopulations if only samples of two time points are available. To further explore the alteration in subpopulations and functions of T cells, a viSNE analysis was performed using the panel including CD4, CD8A, CCR4, CCR5, CCR7, CD45RA, CD45RO, PD-1, LAG3, TIM3, CTLA4, OX40, 4-1BB, ICOS, CD28, FAS, and CXCR3. To validate the findings of high-throughput mRNA assay, a additional antibody panel were used in the mass cytometry analysis. In this panel, TNFα, IFNγ, TGFβ, pNFkB_P65, CD86, GM-CSF, IL-10, CD284, IL-2, MHCII, CD206 antibodies were added, while functional markers which have been already detected in the first antibody panel, including LAG3, TIM3, CTLA4, OX40,4-1BB, ICOS, CD28, FAS, were removed.
Bead-based normalization was performed using the CyTOF software. Manual gating, SPADE analysis and vi-SNE analysis were performed using the Cytobank software (https://premium.cytobank.org). Repeated measurements ANOVA were performed using R software and its CAR package. Bubble plot, box plot and heat maps were prepared using R software (version 3.3), pheatmap and ggplot2 package and Graphpad Prism 5 software. KEGG pathway enrichment were performed using the DAVID software (version 6.8, https://david.ncifcrf.gov/). Protein-protein interaction were analyzed using the STRING software (version 10.5, https://string-db.org/) and cytoscape software (version 3.4.0).
Results
Patients
In this study, 24 recurrent HCC patients were enrolled, and all of them received second RFA. Three of them were excluded because of loss of follow-up after the second RFA. The number of isolated PBMC of four patient was insufficient to complete mass cytometry. The clinical characteristics of the enrolled patients were shown in Table 1. Blood samples were collected at four time points: the time before the first RFA, 3 days after the first RFA, before the second RFA, and 3 days after the second RFA. Then mass cytometry was performed on the PBMC samples from 17 patients. High-throughput mRNA analysis was performed on the pre-RFA and post-RFA serum collected from 4 patients.
Manual gating and SPADE analysis
To explore the multi-dimension mass cytometry data, we analyzed CD45+ blood leukocytes through the implementation of manual gating and SPADE analysis. The gating strategy employed were displayed in Figure 1A. The significant altered PBMC subpopulations identified by manual gating were shown in Supplementary Figure S1.
Figure 1
In the SPADE analysis, numbers of clusters were set to 200 to avoid omission of rare and small cell subpopulation. These clusters were then annotated according to the median expression of cell markers. The analysis revealed a diverse range of clusters, including 43 CD4 T cell clusters, 50 CD8 T cell clusters, 10 Treg clusters, 6 NK clusters, 2 NKT clusters, 26 Monocyte clusters, 3 DPT (double positive T cells) clusters, 11 DNT (double negative T cell) clusters, 18 B cell clusters and one DC (Dendritic Cell) clusters (Figure 1B). Eight HLA-DR+CD3-CD16-CD45RA-CD44lowCD11Alow antigen-presenting cells (APC) clusters and 22 CD3- HLA-DR-CD16- clusters were not recognized by our panel. The detailed marker expression of every SPADE node was shown in Supplementary Figure S2.
The alteration of APCs after RFA treatment
APCs are major players in initiation of anti-tumor immune response. Among CD45+ cells, HLA-DR+ APCs constituted a substantial proportion, accounting for 25.11% (25 and 75 percentiles: 17.88-39.73%). Within the HLA-DR+ APC population, monocytes constituted 13.23% (8.18-25.8%), DCs accounted for 0.09% (0.03-0.17%), and B cells represented 4.02% (2.8-7.84%), and a big proportion of APCs cannot be recognized by our panel, which took 4.18% (2.08-7.03%) of all the CD45+ cells.
Following RFA treatment, there was a rapid and significant increase in the percentage of monocytes, which subsequently returned to normal levels prior to the second RFA (P=0.026, Supplementary Figure S1A). A sustained reduction of B cell count was observed after the first RFA (P=0.021, Supplementary Figure S1B), and the B cell count remained at a low level before and after the secondary RFA. The B cell count was also reduced in 13 clusters defined by SPADE analysis.
Populational and functional landscape of CD3+ T cells
Following the first RFA treatment, there was a significant reduction in the percentage of CD4+ T cells within the CD45+ cells, and this reduction was observed again after the second RFA (P=0.024, Supplementary Figure S1C). CD8+ T cells showed the same alteration as that of CD4+ cells (P=0.03, Supplementary Figure S1D). We then studied the subpopulation of CD4+ and CD8+ T cells using SPADE analysis. Primary SPADE analysis yielded 45 CD4 T cell clusters, 49 CD8 T cell clusters, 3 DPT clusters and 11 DNT clusters. Among the CD4+ T cell clusters, 10 clusters represented CD45RAhighCD45ROlowCCR7- CD4 Teff (T effector) clusters, 4 clusters represented CD45RAhighCD45ROlowCCR7+ CD4 T naïve cells, 15 clusters represented CD45RAlow CD45ROhighCCR7- CD4 Tem (T effector memory) cells and 16 clusters represented CD45RAlow CD45ROhighCCR7+ CD4 Tcm (T central memory) cells. Among the CD8 T cell clusters, 14 clusters represented CD45RAhighCD45ROlowCCR7- CD8 Teff cells, 20 clusters represented CD45RAhighCD45ROlowCCR7+ CD8 T naïve cells, 9 clusters represented CD45RAlow CD45ROhighCCR7- CD8 Tem cells and 6 clusters represented CD45RAlow CD45ROhighCCR7+ CD8 Tcm cells.
Furthermore, we annotated the T effector clusters according to the expression of the surface marker CCR4, CCR5, and CCR7. Within the CD4 T effector cell clusters, five clusters were identified as CCR4-CCR5-CCR7- Th0 cells, one cluster was identified as CCR4-CCR5+CCR7- Th1 cell and four clusters were identified as CCR4+CCR5-CCR7- Th2 cells. Within the CD8 T effector cell clusters, 11 clusters were identified as CCR4-CCR5-CCR7- Tc0 cells and 3 clusters were identified as CCR4+CCR5-CCR7- Tc2 clusters. (Figures 2A, B, Supplementary Table S1).
Figure 2
We then analyzed the change in cell count of T cells at the four designated time points. Six CD4 Tem clusters showed a reduced population after the first RFA. After the second RFA, the number of Tem clusters showed no change or moderate increase (P<0.05 for all the six Tem clusters) (Figure 3A). Two CD4 Tcm clusters showed reduced populations after the first RFA, which were then restored to a lower level before the second RFA started (P=0.002, 0.021, respectively, Figure 3A). Two CD4 T naive clusters showed increased cell counts at time points 3 and 4 (P=0.023, 0.042, respectively, Figure 3A). With respect to CD4 effector T cells, 3 Th2 clusters showed significant population alteration across the 4 time points. All the Th2 clusters showed a more than 3-fold decrease after the first RFA. Cluster 38 and 58 had a slight decrease while Node 73 showed a significant increase after the second RFA (P=0.017, 0.026, 0.024, respectively, Supplementary Table S1). Three Tc0 clusters showed major or moderate decrease after the first and second RFA (P=0.022, 0.015, 0.006, respectively) (Supplementary Table S1).
Figure 3
Next, we focused on the expression of co-stimulatory molecule of T cells. Significant alterations of CD28 expression were found in 7 T cells clusters, including two CD4 Tem clusters, two CD8 T naïve clusters, one DPT cluster, and one DNT cluster (P<0.05 Figure 3B). However, following the first RFA, these seven clusters displayed decreased CD28 expression. Interestingly, 2 CD8 T Naïve clusters showed increased CD28 expression after the second RFA, while the other 5 clusters showed decreased CD28 expression again after the second RFA (Figure 3B). Cluster 167, identified as a CD4 Tnaive cluster, showed significant ICOS expression after the first RFA, and remained in a high ICOS expression level before and after the second RFA (P=0.049). A reduction in PD-1 expression of cluster 167, a CD4 Tnaive cluster, and cluster 193, a CD4 Tcm cluster, was found after the first and second RFA (P=0.022). No significant alterations of 4-IBB, OX-40 were found in CD3+ T cell clusters (Supplementary Table S2).
Furthermore, we investigated the expression of inhibitory molecules, and expression of PD-1 were altered in Node 19, 56, 167,190 and 193, which were identified as CD4Tcm, CD8 Tem, CD4 Tnaive, Treg and CD4 Tcm, respectively (Supplementary Table S2). An elevated PD-1 expression after the first and second RFA were found in Node 19 and 56. While no significant alterations were found in the expressions of TIM3, CTLA-4 and LAG-3.
Furthermore, a viSNE analysis was performed in 5 patients with more than 10 thousand CD3+ T cells at all four time points to explore the functional subpopulations of T cells (Supplementary Figure S3), We found that CD3+ T cells of patient 9 showed remarkable CTLA-4 expression and the expression of CTLA-4 was greatly reduced, accompanied with an increase of CD28 expression at time points 3 and 4 (Supplementary Figure S4). In the other four patients, CTLA-4 was mainly expressed in DPT cells.
NK cells showed down-regulated ICOS and OX40 expression after the first RFA
NKs cells are crucial components in the anti-tumor immune response after RFA, which could eradicate the tumor cells directly or induce the apoptosis of tumor cells by secreting IFN-α and TNF-α (). Notably,we examined the cell count of NK cells and NKT cells and found no significant alterations across the four time points.
High-throughput mRNA assays revealed immunosuppressant mRNA profile after RFA
To extend our knowledge of impact of RFA on the functional alteration of PBMC, we proceeded to examine the mRNA profile of PBMC using a high-throughput gene assay system. The expression of a panel of 148 immune-related genes was examined in samples from four HCC patients. Overall, we discovered 21 differently expressed genes in PBMC collected after RFA, compared with the PBMC collected before RFA (Figure 4A). Interestingly, CD2 and CD8A, which is known as T cell-specific surface markers, and GZMB, GZMH, GZMK, and NKG7, which is closely related to the cytotoxic function of T cells, showed decreased mRNA expression in PBMC after RFA. Also, the expression of chemokine CCL4 and CCL5 significant decreased (Figure 5A). Furthermore, the protein-protein interaction analysis again revealed an inhibited PPI network relating to the T cell function (Figure 4B). Further functional pathway analysis using DAVID software demonstrated that the Natural killer cell mediated cytotoxicity, and the JAK-STAT signaling pathway were significantly inhibited after RFA (Figure 4C, Supplementary Table S3). These data again showed an immunosuppressant effect of RFA on PBMC, especially on T cells.
Figure 4
Figure 5
In vitro stimulating of PBMC by post-RFA serum of HCC patients
We have demonstrated above that RFA could cause the population and gene signature of PBMC to change towards an immunosuppressant status. Given that the major mechanism of RFA involves the direct destruction of tumor tissue and cause the necrosis in the remaining cancer cells, we hypothesized that the components released by the necrosis tumor tissue is the main cause of the immunosuppression effect of RFA. We then treated the PBMC from a health donor with before-RFA and post-RFA serum from four HCC patients to validate the immunosuppression effect of post-RFA serum. We found that after treated by post-RFA serum for 48 hours, the percentage of CD8 Teff and CD8 Tem in CD45+ cells significantly decreased (Figure 5A; P= 0.002 and 0.040, respectively).
In terms of functional markers, the expression of MHCII and TLR4 of monocytes significantly elevated compared to the PBMC treated with post-RFA serum, indicating an enhanced antigen-presenting function of monocytes (Figure 5C; P= 0.033 and 0.022, respectively). However, no significant changes were found in GM-CSF, TNFα, IFNγ and p-NFκB expression of monocytes. Instead, the secretion of VEGF, which is increased in tumor-associated macrophages, significantly increased in monocytes (Figure 5G, P<0.001). Increased VEGF expression was also found in B cell, Treg, and CD4 T Naïve (Figure 5G P= 0.024, 0.002 and 0.038, respectively). Meanwhile, a significant reduction of CD5 expression in a wide span of T cells was found, including CD8 Tcm, CD8 Teff, CD8 Tem, CD8 T Naïve, CD4 Teff and CD4 Tem (Figure 5E; P=0.041, 0.028, 0.026 0.023, 0.046 and 0.042, respectively), as well as a significant reduction of CD161 expression in CD4 Teff, CD8 Tem, and CD8 Teff (Figure 5F; P=0.025, 0.025 and 0.038, respectively), were also found. These data further support the immunosuppressive effect of post-RFA serum.
To further locate the key component with immunosuppression effect, and find out whether the immunosuppression effect of post-RFA serum is caused by the components released by necrotic tumor tissue, we removed serum components bigger than 10KD, mostly of which are proteins, in the before-RFA and post-RFA serum with protein chromatography, and treated the PBMC from the same donor for 48h, followed by mass cytometry assay. We detected no population alteration in the post-RFA serum treated PBMC (Figure 5B). Meanwhile, we found no alteration in VEGF or CD161 expression, and the expression of MHCII and TLR4 of monocytes showed no difference (Figure 5D).
Discussion
As a method for a precise elimination of tumor loci with minimal invasion and fast recovery, RFA is widely used in the treatment of HCC. Apart from its ability to induce coagulate necrosis directly, RFA also exerts an impact on the anti-tumor immune response (). The immune phenotype of PBMC, which is the direct reflection of the global immune status, was reported to remarkably correlate with cancer recurrence (, ). In this article, we studied the influence of RFA on anti-tumor immune response in the PBMC of recurrent HCC patients using a combination of mass cytometry and high-throughput RNA assay.
Enhanced antigen-presenting capacity of monocyte after RFA
RFA could activate the antigen-presenting cells by increasing exposure to tumor antigens (). Following RFA treatment, necrotic HCC cells release tumor specific antigens and other heat stressed cells release acute phase proteins, pro-inflammatory cytokines and heat shock proteins (HSPs) (, ). HSPs function as intracellular molecular chaperones facilitating the binding of tumor peptide antigens. The HSP-tumor antigen complex can be effectively recognized by APCs, promoting their maturation (), In our study, we found that the number of monocytes in peripheral blood of recurrent HCC patients significantly increased after the RFA, and culturing PBMC with post-RFA serum didn’t directly increase the population of monocytes, these findings suggest that RFA could enhance the recruitment of monocytes from the lymphatic organs. Additionally, we discovered that post-RFA serum of recurrent HCC patients induced increased expression of MHCII and TLR4 in monocytes. However, after removing components bigger than 10KD from the post-RFA serum, it demonstrated no influence on the phenotype of monocytes. MHC II and TLR4 are constitutionally expressed by antigen presenting cells and present antigenic peptides to T lymphocytes. MHC II and TLR4 expression correlated significantly with the activity of intratumoral inflammation and the infiltration of CD4+ T lymphocytes (, ). These findings suggest that cancer-derived components in the serum of recurrent HCC patients can also enhance the antigen-presenting capacity of monocytes.
Alteration in peripheral T cell subpopulation and recruitment of T cells
CD4+T cells and CD8+ T cells play pivotal roles in the anti-cancer immune response. Upon recognition of tumor specific antigen presented by APCs, CD4 and CD8 naïve T cells undergo differentiation and become effector T cells. CD4 Teff could induce humeral immune response and enhance cellular immune response by secreting pro-inflammatory cytokines, and CD8 Teff could directly eliminate cancer cells by secreting granzymes, perforin and granulysin. T cell exhaustion has been thought to be an important cause of the failure of anti-cancer immunity. Chronic infection and cancer are believed to be common factors inducing T cell exhaustion (–).
An increase in peripheral T cells following RFA has been reported in several studies (–), and the increased level of peripheral T cells, especially the number of CD8 T cells, could improve the prognosis of the HCC patients (). However, our current study found a temporary decrease in peripheral CD4 and CD8 T cells after RFA in recurrent HCC patients, which further declined to lower levels compared to pre-RFA PBMCs. Previous studies suggested that incomplete RFA induces inflammation and accelerates tumor progression (). In this context, HCC recurrence after RFA, similar to incomplete radiofrequency ablation, may affect inflammation.
We then examined the expression of co-stimulatory molecules in peripheral T cells and found that only a small proportion of CD8 Tnaive showed increased CD28 expression after RFA, and only one out of four CD4 Tnaive clusters showed increased ICOS expression. Also, no alterations in OX40 and 4-IBB expression were found. Meanwhile, we found that the post-RFA serum of recurrent HCC patients could reduce population of CD8 Tem and CD8 Teff, and the CD161 expression in CD4 Teff, CD4 Tnaive, CD8 Tem, CD8Teff, and CD8 Tnaive. CD161 is a marker of all human IL17 producing T cell subsets, and it is also a pan-cancer prognostic factor with favoring outcome (, ), and these inhibitory function of RFA over T cells could be eliminated by removal of serum components which is bigger than 10KD. We also found the down-regulation of mRNA expression of GZMB (granzyme B), GZMK (granzyme K), GZMH (granzyme H), which are key enzymes in the cytotoxic function of CD8+T cells (, ), and the down-regulation of ISG15 (Interferon-Stimulated Gene 15), NKG7 (natural killer cell granule protein 7), which were reported to be related to T cell and NK cell related immune response (). RFA was reported to enhance host immune response by increasing tumor-associated antigen (TAA)-specific T cell responses (, ). But on the contrary, we observed a significant reduction in the population and cytotoxic function of CD8 Teff and CD8 Tem cells after incomplete RFA, despite these cells may receive more tumor antigens because of the enhancement of the antigen-presenting capacity of monocytes in incomplete RFA. These findings suggest that the inhibitory effect of proteins secreted by necrotic tumor tissues may be an important reason for inhibition of T cell- related anti-tumor response and tumor recurrence. In recent years targeting immune checkpoint therapy showed promising therapeutic effect in HCC. Wen et al. () reported that combining RFA with anti-PD-1 showed improved RFS and was deemed safe for patients with recurrent HCC who had previously undergone RFA treatment alone. But some studies also suggest that incomplete radiofrequency ablation accelerates tumor progression and hinders PD-1 immunotherapy by inducing the accumulation of tumor-associated macrophages. Overall, the combination of RFA and anti-PD-1 therapy represents a promising treatment strategy for recurrent HCC patients, but further studies are warranted to enhance the efficacy of this treatment approach.
Up-regulation of VEGF in post-RFA serum treated lymphocytes
In this study, we observed a significant increase in VEGF production in monocytes, Tregs, B cells and CD4 Tnaive cells showed significant increase following treated by post-RFA serum, and this manifest could also be eliminated by removal of serum components which is bigger than 10KD. VEGF is best known for its role in the tumor angiogenesis, and it is known to be extensively secreted by tumor associated macrophages and Tregs (). It has also been reported to selectively kill CD8+T cells and produce a immunosuppressant microenvironment for cancer cells (). These findings suggest that cancer-derived components in the serum may contribute to enhanced immunosuppression and angiogenesis by promoting the secretion of VEGF.
In conclusion, by using a combination of mass cytometry and high-throughput mRNA assay, our study demonstrated that in recurrent HCC patients, RFA therapy can enhance the antigen-presenting capacity of monocytes, However, it is insufficient to induce a robust anti-tumor immune response. This is because RFA therapy can inhibit the anti-tumor immune response by declining the population of T cells, inhibiting the activation of T cells and downregulating the expression of CD161 and CD5 in various T cell subpopulations. Additionally, tumor-derived components contribute to the creation of an immunosuppressive microenvironment by promoting the secretion of VEGF in Monocytes, Tregs, B cells and CD4 T naïve cells. Therefore, the immunosuppression caused by tumor-derived big molecules may be an important reason of HCC recurrence after RFA therapy, and application of immunopotentiating therapy in recurrent HCC patients should be further studied.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by Medical Ethics Committee of YouAn Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
YZ: Data curation, Formal Analysis, Writing – original draft, Writing – review & editing. TY: Methodology, Investigation, Writing – original draft, Writing – review & editing. YO: Methodology, Investigation, Writing – original draft, Writing – review & editing. FW: Methodology, Investigation, Writing – original draft, Writing – review & editing. DL: Methodology, Investigation, Writing – original draft, Writing – review & editing. LQ: Methodology, Investigation, Writing – original draft, Writing – review & editing. ZX: Methodology, Investigation, Writing – original draft, Writing – review & editing. WR: Resources, Writing – original draft, Writing – review & editing. KL: Resources, Writing – original draft, Writing – review & editing. JZ: Data curation, Writing – original draft, Writing – review & editing. FL: Data curation, Writing – original draft, Writing – review & editing. YS: Data curation, Writing – original draft, Writing – review & editing. ZY: Formal Analysis, Methodology, Project administration, Supervision, Writing – review & editing. DC: Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Supervision, Writing – review & editing. WW: Conceptualization, Data curation, Investigation, Methodology, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Beijing Natural Science Foundation project (7234384, Z230008), Peking Union Medical College Hospital Special Clinical Research Project (LY22B0100010288), and the National Natural Science Foundation of China (81672026, 81601796).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2024.1339213/full#supplementary-material
Supplementary Figure 1Significant altered subpopulation of PBMC identified by manual gating (A) Population of B cells significantly reduced after the first RFA (P=0.021), and the B cell count remained at a low level before and after the secondary RFA. (B) The percentage of monocytes showed a fast increase after RFA, and was reduced to the normal level before the second RFA (P=0.026). (C) and (D). The percentage of CD4+ T cells and CD8 Tcells in CD45+ cells was significantly reduced after both the first and the second RFA (P=0.024 and 0.03, respectively).
Supplementary Figure 2Expression of clustering markers of SPADE nodes in patients one. Median marker expression were marked using the color of the nodes and cell counts were marked with the cell count of the nodes. The cell type of every nodes were identified according to its clustering marker expression. The expression of nodes with small cell counts may vary slightly and these nodes are annotated according to the expression of these nodes in other fcs files with larger cell count.
Supplementary Figure 3Clustering markers and functional markers of T cells. A viSNE analysis was performed on CD3+ T cells. 10 T cell related cells markers were included in the analysis to generate two variables (tSNE1 and tSNE2) to differentiate the T cells populations. Plots are colored according to 10 populational and functional markers. Also, CXCR3 and CCR4+ showed similar expression pattern in all these five patients.
Supplementary Figure 4CTLA4 and CD28 expression of T cells of patients 9. Compared with other patients, patient 9 showed remarkable CTLA-4 expression in CD3+ T cells and the expression of CTLA-4 greatly reduced, accompanied with increase of CD28 expression in time point 3 and 4, while in the other four patients, the CTLA-4 were mainly expressed in DPT cells.
References
1
El-SeragHBRudolphKL. Hepatocellular carcinoma: epidemiology and molecular carcinogenesis. Gastroenterology (2007) 132:2557–76. doi: 10.1053/j.gastro.2007.04.061
2
SiegelRLMillerKDJemalA. Cancer statistics, 2017. CA Cancer J Clin (2017) 67:7–30. doi: 10.3322/caac.21387
3
MizukoshiEYamashitaTAraiKSunagozakaHUedaTAriharaFet al. Enhancement of tumor-associated antigen-specific T cell responses by radiofrequency ablation of hepatocellular carcinoma. Hepatology (2013) 57:1448–57. doi: 10.1002/hep.26153
4
TampakiMDoumbaPPDeutschMKoskinasJ. Circulating biomarkers of hepatocellular carcinoma response after locoregional treatments: New insights. World J Hepatol (2015) 7:1834–42. doi: 10.4254/wjh.v7.i14.1834
5
RozenblumNZeiraEBulvikBGourevitchSYotvatHGalunEet al. Radiofrequency ablation: inflammatory changes in the periablative zone can induce global organ effects, including liver regeneration. Radiology (2015) 276:416–25. doi: 10.1148/radiol.15141918
6
ChenHLauMCWongMTNewellEWPoidingerMChenJ. Cytofkit: A bioconductor package for an integrated mass cytometry data analysis pipeline. PloS Comput Biol (2016) 12:e1005112. doi: 10.1371/journal.pcbi.1005112
7
BironCANguyenKBPienGCCousensLPSalazar-MatherTP. Natural killer cells in antiviral defense: function and regulation by innate cytokines. Annu Rev Immunol (1999) 17:189–220. doi: 10.1146/annurev.immunol.17.1.189
8
WuF. Heat-based tumor ablation: role of the immune response. Adv Exp Med Biol (2016) 880:131–53. doi: 10.1007/978-3-319-22536-4_8
9
AriharaFMizukoshiEKitaharaMTakataYAraiKYamashitaTet al. Increase in CD14+HLA-DR -/low myeloid-derived suppressor cells in hepatocellular carcinoma patients and its impact on prognosis. Cancer Immunol Immunother (2013) 62:1421–30. doi: 10.1007/s00262-013-1447-1
10
ZerbiniAPilliMFagnoniFPelosiGPizziMGSchivazappaSet al. Increased immunostimulatory activity conferred to antigen-presenting cells by exposure to antigen extract from hepatocellular carcinoma after radiofrequency thermal ablation. J Immunother (2008) 31:271–82. doi: 10.1097/CJI.0b013e318160ff1c
11
SchuellerGKettenbachJSedivyRStiftAFriedlJGnantMet al. Heat shock protein expression induced by percutaneous radiofrequency ablation of hepatocellular carcinoma in vivo. Int J Oncol (2004) 24:609–13. doi: 10.3892/ijo.24.3.609
12
AliMYGrimmCFRitterMMohrLAllgaierHPWethRet al. Activation of dendritic cells by local ablation of hepatocellular carcinoma. J Hepatol (2005) 43:817–22. doi: 10.1016/j.jhep.2005.04.016
13
GaidaMMWelschTHerpelETschaharganehDFFischerLSchirmacherPet al. MHC class II expression in pancreatic tumors: a link to intratumoral inflammation. Virchows Arch (2012) 460:47–60. doi: 10.1007/s00428-011-1175-x
14
ChaudhuriAWilsonNSYangBPaler MartinezALiuJZhuCet al. Host genetic background impacts modulation of the TLR4 pathway by RON in tissue-associated macrophages. Immunol Cell Biol (2013) 91:451–60. doi: 10.1038/icb.2013.27
15
SenDRKaminskiJBarnitzRAKurachiMGerdemannUYatesKBet al. The epigenetic landscape of T cell exhaustion. Science (2016) 354:1165–9. doi: 10.1126/science.aae0491
16
SakuishiKNgiowSFSullivanJMTengMWKuchrooVKSmythMJet al. TIM3+FOXP3+ regulatory T cells are tissue-specific promoters of T-cell dysfunction in cancer. Oncoimmunology (2013) 2:e23849. doi: 10.4161/onci.23849
17
YuLHLiNShiJGuoWXWuMCChengSQ. Does anti-HBV therapy benefit the prognosis of HBV-related hepatocellular carcinoma following hepatectomy? Ann Surg Oncol (2014) 21:1010–5. doi: 10.1245/s10434-013-3320-z
18
WidenmeyerMShebzukhovYHaenSPSchmidtDClasenSBossAet al. Analysis of tumor antigen-specific T cells and antibodies in cancer patients treated with radiofrequency ablation. Int J Cancer (2011) 128:2653–62. doi: 10.1002/ijc.25601
19
HanslerJWissniowskiTTSchuppanDWitteABernatikTHahnEGet al. Activation and dramatically increased cytolytic activity of tumor specific T lymphocytes after radio-frequency ablation in patients with hepatocellular carcinoma and colorectal liver metastases. World J Gastroenterol (2006) 12:3716–21. doi: 10.3748/wjg.v12.i23.3716
20
NapoletanoCTaurinoFBiffoniMDe MajoACoscarellaGBellatiFet al. RFA strongly modulates the immune system and anti-tumor immune responses in metastatic liver patients. Int J Oncol (2008) 32:481–90. doi: 10.3892/ijo.32.2.481
21
HiroishiKEguchiJBabaTShimazakiTIshiiSHiraideAet al. Strong CD8(+) T-cell responses against tumor-associated antigens prolong the recurrence-free interval after tumor treatment in patients with hepatocellular carcinoma. J Gastroenterol (2010) 45:451–8. doi: 10.1007/s00535-009-0155-2
22
ShiLWangJDingNZhangYZhuYDongSet al. Inflammation induced by incomplete radiofrequency ablation accelerates tumor progression and hinders PD-1 immunotherapy. Nat Commun (2019) 10:5421. doi: 10.1038/s41467-019-13204-3
23
GentlesAJNewmanAMLiuCLBratmanSVFengWKimDet al. The prognostic landscape of genes and infiltrating immune cells across human cancers. Nat Med (2015) 21:938–45. doi: 10.1038/nm.3909
24
MaggiLSantarlasciVCaponeMPeiredAFrosaliFCromeSQet al. CD161 is a marker of all human IL-17-producing T-cell subsets and is induced by RORC. Eur J Immunol (2010) 40:2174–81. doi: 10.1002/eji.200940257
25
MouchaccaPChassonLFrickMForayCSchmitt-VerhulstAMBoyerC. Visualization of granzyme B-expressing CD8 T cells during primary and secondary immune responses to Listeria monocytogenes. Immunology (2015) 145:24–33. doi: 10.1111/imm.12420
26
BadeBBoettcherHELohrmannJHink-SchauerCBratkeKJenneDEet al. Differential expression of the granzymes A, K and M and perforin in human peripheral blood lymphocytes. Int Immunol (2005) 17:1419–28. doi: 10.1093/intimm/dxh320
27
BurksJReedREDesaiSD. Free ISG15 triggers an antitumor immune response against breast cancer: a new perspective. Oncotarget (2015) 6:7221–31. doi: 10.18632/oncotarget.3372
28
NobuokaDMotomuraYShirakawaHYoshikawaTKuronumaTTakahashiMet al. Radiofrequency ablation for hepatocellular carcinoma induces glypican-3 peptide-specific cytotoxic T lymphocytes. Int J Oncol (2012) 40:63–70. doi: 10.3892/ijo.2011.1202
29
WenZWangJTuBLiuYYangYHouLet al. Radiofrequency ablation combined with toripalimab for recurrent hepatocellular carcinoma: A prospective controlled trial. Cancer Med (2023) 12:20311–20. doi: 10.1002/cam4.6602
30
ColegioORChuNQSzaboALChuTRhebergenAMJairamVet al. Functional polarization of tumour-associated macrophages by tumour-derived lactic acid. Nature (2014) 513(7519):559–63. doi: 10.1038/nature13490
31
MotzGTSantoroSPWangLPGarrabrantTLastraRRHagemannISet al. Tumor endothelium FasL establishes a selective immune barrier promoting tolerance in tumors. Nat Med (2014) 20:607–15. doi: 10.1038/nm.3541
Summary
Keywords
mass cytometry, radiofrequency ablation, recurrent hepatocellular carcinoma, PBMC, immunosuppresion
Citation
Zhao Y, Yang T, Ouyang Y, Rao W, Liu K, Zheng J, Lv F, Shi Y, Wang F, Liu D, Qiao L, Xia Z, Zhang Y, Chen D and Wang W (2024) Radiofrequency ablation plays double role in immunosuppression and activation of PBMCs in recurrent hepatocellular carcinoma. Front. Immunol. 15:1339213. doi: 10.3389/fimmu.2024.1339213
Received
06 December 2023
Accepted
11 January 2024
Published
29 January 2024
Volume
15 - 2024
Edited by
Ming Yi, Zhejiang University, China
Reviewed by
Wang Zehua, Shandong University, China
Jia-quan Zhou, Hainan General Hospital, China
Updates
Copyright
© 2024 Zhao, Yang, Ouyang, Rao, Liu, Zheng, Lv, Shi, Wang, Liu, Qiao, Xia, Zhang, Chen and Wang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wenjing Wang, wangwenjing85321@ccmu.edu.cn; Dexi Chen, dexichen@ccmu.edu.cn; Yushi Zhang, zhangyushipumch@163.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.