Abstract
Anti-PD-1 therapy has shown promising outcomes in the treatment of different types of cancer. It is of fundamental interest to analyze the efficacy of anti-PD-1 therapy in cancer patients infected with hepatitis B virus (HBV) since the comorbidity of HBV and cancer is widely documented. We designed a multicenter retrospective study to evaluate the efficacy of anti-PD-1 therapy on non-liver cancer patients infected with HBV. We found anti-PD-1 therapy achieved much better outcomes in HBV+ non-liver cancer patients than their HBV– counterparts. We performed single-cell RNA sequencing (scRNA-seq) on peripheral blood mononuclear cells (PBMCs) from esophageal squamous cell carcinoma (ESCC) patients. We found both cytotoxicity score of T cells and MHC score of B cells significantly increased after anti-PD-1 therapy in HBV+ ESCC patients. We also identified CX3CR1high TEFF, a subset of CD8+ TEFF, associated with better clinical outcome in HBV+ ESCC patients. Lastly, we found CD8+ TEFF from HBV+ ESCC patients showing higher fraction of Exhaustionhi T than their HBV– counterpart. In summary, anti-PD-1 therapy on HBV+ non-liver cancer patients is safe and achieves better outcomes than that on HBV– non-liver cancer patients, potentially because HBV+ patients had higher fraction of Exhaustionhi T, which made them more efficiently respond to anti-PD-1 therapy.
Similar content being viewed by others
Introduction
Immune checkpoint blockade (ICB) therapy is a promising immunotherapy approach, among which anti-programmed cell death-1 (anti-PD-1) therapy yielded promising outcomes in treatment of dozen of types of cancer, including melanoma, hepatocellular carcinoma (HCC), gastric cancer, and esophageal cancer [1,2,3,4,5,6]. However, only a fraction of cancer patients responded to anti-PD-1 therapy. Of note, the efficacy of anti-PD-1 therapy is affected by many factors including tumor microenvironment, commensal microbiota, antibiotics, steroids, and viral infections [3, 7,8,9]. In particular, about 54% of HCC cases were associated with hepatitis B virus (HBV) infection [10], and HCC patients with HBV infection were excluded in early anti-PD-1 therapy trials due to concerns about the reactivation of HBV and uncertain about immune microenvironment. However, enrollment of HCC patients with HBV infection showed that low and intermediate HBV-DNA level (<500 IU/mL) did not impact the efficacy and safety of anti-PD-1 therapy [1, 11], even HCC patients with high HBV-DNA level did not show increased incidence of HBV-associated hepatitis [12, 13].
Simultaneously, the success of anti-PD-1 in cancer therapy potentially indicates that anti-PD-1 might be effective for treating infectious diseases since chronic infection also showed high expression of PD1 and T-cell exhaustion [14]. Over 350 million people in the world and 70 million people in China were infected with HBV, which was one of the most common comorbidities with non-liver cancer [15]. Recently, many studies analyzed anti-PD-1 therapy on patients with comorbidity of non-liver cancer and HBV infection [16,17,18,19,20,21,22,23,24,25,26]. Some of these studies only focused on the safety of anti-PD-1 therapy due to concerns about HBV reactivation [16, 21, 25, 26]. Other studies analyzed the efficacy of anti-PD-1 therapy [17,18,19,20, 22,23,24], which only enrolled a few non-liver cancer patients with HBV virus load (<7 samples). Since these studies did not design HBV+non-liver cancer patients and their HBV– counterparts matched case-control study [17,18,19, 23], it is still unknown whether patients with comorbidity of non-liver cancer and HBV infection could achieve similar efficiency compared with their HBV– counterparts.
Here, we performed a multicenter retrospective study on 84 non-liver cancer patients, including 35 HBV+ patients and 49 HBV– patients. We further performed single-cell RNA sequencing (scRNA-seq) on six peripheral blood mononuclear cells (PBMCs) samples from three esophageal squamous cell carcinoma (ESCC) patients (pre- and post-anti-PD-1 therapy), to analyze their response to anti-PD-1 therapy. We further analyzed the differences in the immune microenvironment between HBV+ ESCC patients and HBV– ESCC patients to explore the mechanism underlying their different response to anti-PD-1 therapy.
Results
HBV– patients and HBV+ patients are well matched
In order to analyze the effect of HBV infection on anti-PD-1 therapy in non-liver cancer patients, a total of 35 HBV+ non-liver cancer patients were screened from 7,231 cancer patients who visited during 2018–2021 (Fig. 1A, Supplementary Fig. S1A, Supplementary Table 1). We selected 49 HBV– non-liver cancer patients that matched HBV+ non-liver cancer patients in terms of cancer types, age, and gender (Fig. 1A, Supplementary Fig. S1A, Supplementary Table 2). The median age of HBV+ and HBV– non-liver cancer patients was 55-year-old (range: 33–74 years, with 7 patients ≥65) and 55-year-old (range: 37–79 years, with 9 patients ≥65), respectively. Both HBV+ and HBV– non-liver cancer patients were mainly males, with 27 males (77.1%) and 34 males (69.4%), respectively (Supplementary Table 3). In our study, 42.9% of HBV+ patients and 34.7% of HBV– patients received anti-PD-1 therapy as the first-line therapy, with no significant difference between the two groups (p-value = 0.499) (Supplementary Table 3).
Anti-PD-1 therapy achieves favorable outcomes in HBV+ non-liver cancer patients
Surprisingly, anti-PD-1 therapy achieved significantly better outcomes in HBV+ non-liver cancer patients than in HBV– non-liver cancer patients (p-value = 0.0052), with 15 (42.86%) HBV+ non-liver cancer patients showing partial response (PR) while only 7 (14.29%) HBV– patients showing PR (Fig. 1B, Supplementary Tables 1–3). Moreover, 22 (44.9%) HBV– non-liver cancer patients developed progressive disease (PD), whereas only 6 (17.1%) HBV+ non-liver cancer patients developed PD.
Based on the objective response rate (ORR), the odds ratio (OR) of HBV+ non-liver cancer patients compared with their HBV– counterparts was 4.41 (95% CI; 1.426–14.985). HBV+ patients showed much longer overall survival (OS), further indicating that HBV+ patients have much better outcome than their HBV– counterparts (Fig. 1B, Supplementary Fig. S1B). Although these patients received different PD-1 antibody treatments, there was no significant difference in the outcome of different PD-1 antibodies (Fig. 1B). The outcomes of combination anti-PD-1 therapy showed higher ORR events (30.19%) than monotherapy (19.35%), although the difference is not significant (p-value = 0.315) (Fig. 1B) potentially due to the small sample size in this study. Furthermore, OR for first-line treatment compared with second- and beyond-line is 2.49 (95%CI; 0.83–7.679), which is not significant difference and was essentially consistent with a recent study by Hughes et al. [27].
Change of population abundances in PBMCs after anti-PD-1 therapy
Three ESCC patients, with one HBV-negative patient (HBV–#1) and two HBV-positive patients (HBV+#1 and HBV+#2), were enrolled from the Southwest Hospital for exploring their response to anti-PD-1 therapy. HBV+#1 and HBV+#2, both with hepatitis B surface antigen (HBsAg)-positive and hepatitis B core antibody (HBcAb)-positive, had 1910 IU/mL HBV-DNA and 48,600 IU/mL HBV-DNA in blood, respectively (Supplementary Table 4). We performed 10x genomics scRNA-seq on the PBMCs of the three patients before anti-PD-1 therapy. All three patients received sintilimab (anti-PD-1, IBI 308) and paclitaxel liposome as first-line therapy. After one cycle of therapy (4 weeks), blood samples were collected from patients for routine blood test and scRNA-seq library preparation. HBV–#1 achieved stable disease (SD) to the anti-PD-1 therapy after one cycle of therapy. HBV+#1 and HBV+#2 achieved SD and PR, respectively. In total, 44,253 PBMCs from the 3 ESCC patients passed quality control and were clustered into 17 subsets that were annotated according to their specific expression of classic markers (Fig. 2A, B, Supplementary Fig. S2). The fractions of CD4+ TN, CD4+ TM, CD8+ TN, BN, and BM consistently decreased across patients after anti-PD-1 therapy, while the fraction of cMo increased in both HBV–#1 and HBV+#2 (Fig. 2B). We developed population change index (PopIndex) to describe change of immune cell subsets pre- and post-anti-PD-1 therapy. We found the PopIndex of CD4+ TEFF, CD8+ TEFF and plasma cells had pronounced increase in HBV+#2, which may contribute to partial response of HBV+#2 (Fig. 2C). These results were consistent with previous reports that increase of TEFF or increase of plasma cells is correlated with tumor regression [28,29,30].
Change of T-cell features and a T-cell subset after anti-PD-1 therapy
To analyze the effect of anti-PD-1 therapy on cell state of T cells, we identified the significantly differentially expressed genes (DEGs) between samples pre- and post-anti-PD-1 therapy in each of the 6T cell subsets, namely CD4+ TN, CD4+ TM, CD4+ TEFF, CD8+ TN, CD8+ TM, and CD8+ TEFF (Supplementary Fig. S3A, Supplementary Table 5). GO enrichment analyses showed that cytokine signaling and T-cell activation signaling was significantly enriched in HBV+ ESCC patients, particularly in CD8+ TEFF (Supplementary Fig. S3B, Supplementary Table 5) [31]. We further found that cytotoxicity scores of CD8+ TEFF increased significantly after anti-PD-1 therapy in HBV+ ESCC patients, with HBV+#2 showing the most pronounced increase, while cytotoxicity scores of CD8+ TEFF significantly decreased in HBV–#1 (Fig. 3A).
We further analyzed CD8+ TEFF since it has the highest increase of population size and the highest increase of cytotoxicity score in HBV+#2. We clustered CD8+ TEFF and identified two subsets, namely CX3 CR1hi TEFF and CX3 CR1low TEFF (Fig. 3B). Particularly, CX3CR1, a marker of T-cell differentiation, was identified as a predictive biomarker in response to ICB therapy [32, 33]. CX3CR1hi TEFF highly expressed cytotoxicity-associated genes such as GZMA and GZMB (Fig. 3C, Supplementary Table 6), consistent with recent reports that CX3CR1+ TEFF lowly expressed co-inhibitor and highly expressed cytotoxicity-associated genes [32, 33]. In contrast, CX3CR1low TEFF highly expressed activator protein-1 (AP-1) family genes (FOS, FOSB, JUNB, JUN and JUND), which were cancer-related transcription factors [34]. CX3CR1low TEFF highly expressed regulators of G protein signaling (RGS1 and RGS2) that was related to lower infiltration to cancer [35]. It also highly expressed CXCR4 (Fig. 3C) that was correlated with T cells migration and promoted tumorigenesis [36]. These results revealed opposite functions of CX3CR1low TEFF and CX3CR1hi TEFF. CX3CR1low TEFF was mainly enriched in HBV–#1 after anit-PD1 therapy, while CX3CR1hi TEFF were mainly in HBV+patient after anit-PD1 therapy (Fig. 3D), thus the enrichment of CX3CR1hi TEFF is associated with response to anti-PD-1 therapy. We performed survival analysis in esophagus carcinoma (ESCA) patients in The Cancer Genome Atlas (TCGA) [37]. We found the ESCA patients highly expressed CX3CR1low TEFF-specific genes, namely CX3CR1low TEFF signature high patients, had lower survival than CX3CR1low TEFF signature low patients (Fig. 3E). Furthermore, CX3CR1low TEFF signature high patients displayed lower survival than CX3CR1low TEFF signature low patients in BLCA, LUSC, OV, STAD and COAD, supporting that CX3CR1low and CX3CR1hi TEFF have different function in patients (Supplementary Fig. S4).
Change of B/plasma cell features and subsets after anti-PD-1 therapy
Recent studies demonstrated that B cells played an important role in antitumor responses, which presented tumor-derived antigens to T cells and provoke T cells to regress tumor [28]. Therefore, we identified DEGs of each B-cell subset pre- and post-anti-PD-1 therapy in each patient. We found that major histocompatibility complex (MHC)-associated genes were significantly increased in BN and BM in HBV+ patients after anti-PD-1 therapy, while it significantly decreased in HBV–#1 (Fig. 3F and Supplementary Table 7). Furthermore, we found MHC I score and MHC II score of BN and BM in HBV+ patients significantly increased after anti-PD-1 therapy, while there is no change or decrease in HBV–#1 (Fig. 3G). In addition, many immunoglobulin genes such as IGHG3, IGHG4, IGHG1, IGHM, IGLC2, IGLC1 and IGKC were up-regulated in HBV+ ESCC patients and down-regulated in HBV–#1 after anti-PD-1 therapy (Supplementary Fig. S5, Supplementary Table 8). Therefore, B cells showed a higher MHC score and higher expression of immunoglobulin genes in the 2 HBV+ ESCC patients after anti-PD-1 therapy.
HBV+ non-liver cancer patients have higher fraction of Exhaustionhi T before treatment
Our analyses showed that immune function of HBV+ ESCC patients was significantly enhanced after anti-PD-1 therapy, while HBV– patient had little change, which might be caused by prevailing differences in the immune microenvironment between HBV+ patients and HBV– patients. We conducted a systemic comparison of the immune microenvironment between HBV+ ESCC patients and HBV– ESCC patients by integrating scRNA-seq data from Dinh et al. [38]. A total of 40,434 cells from 6 untreated ESCC patients were projected into UMAP plot (Fig. 4A, B). Interestingly, the distribution of exhaustion scores of CD8+ TEFF showed there are two CD8+ TEFF subsets (Fig. 4C), namely Exhaustionhi T and Exhaustionlow T based on whether the normalized exhaustion score is >0 or not. HBV+ ESCC patients had much higher fraction of Exhaustionhi T than HBV– ESCC patients (Fig. 4D). In addition, the cytotoxicity scores in CD8+ TEFF in HBV+ ESCC patients were significantly higher than that in HBV– ESCC patients (Fig. 4E). We further analyzed 16 HBV infected patients from Zhang et al. [39], and found CD8T_c03-CX3CR1 subset reported in their study has the highest similarity to CD8+ TEFF in this study (Supplementary Fig. S6A). We found the proportion of Exhaustionhi T in CD8T_c03-CX3CR1 in HBV+ patients was much higher than HBV– individuals (Supplementary Fig. S6B–D). These results are not affected by anti-PD-1 therapy since we only analyzed the samples before therapy. Since anti-PD-1 therapy is mainly targeted at exhausted T cells and other immunosuppressed immune cells, thus a higher fraction of Exhaustionhi T might provide more therapeutic targets, which could explain anti-PD-1 therapy is more effective in HBV+ non-liver cancer patients.
Discussion
Different from these studies focusing on HBV infection and HCC comorbidity [1, 11, 13, 18, 40], this study focuses on HBV infection and non-liver cancer comorbidity. Our analyses showed patients with HBV+ non-liver cancer did not show increased irAEs than their HBV– counterparts during anti-PD-1 therapy. Intriguingly, HBV+ non-liver patients responded to anti-PD-1 therapy much better than their HBV– counterparts. Similar phenomena have been reported in other virus infections and cancer comorbidities. For example, either patients with Epstein–Barr virus infection and gastric cancer comorbidity, or patients with Epstein–Barr virus infection and non-Hodgkin lymphoma comorbidity, responded better to anti-PD-1 therapy than their Epstein–Barr virus negative counterparts [3, 41]. Finally, our results suggested that anti-PD-1 therapy should be encouraged to treat HBV+ non-liver cancer patients since they could achieve favorable outcome.
To explore the potential mechanism underlying the enhanced response of HBV+ non-liver cancer patients to anti-PD-1 therapy, we performed scRNA-seq on PBMCs and compared the features of each cell subset pre- and post-therapy in the 3 ESCC patients. We found cytotoxicity scores and MHC scores were significantly increased in CD8+ TEFF and B cells in HBV+ ESCC patients after anti-PD-1 therapy, respectively. The enhanced immune function of PBMCs in HBV+ ESCC patients after anti-PD-1 therapy suggested that there is a stronger immunosuppressive microenvironment in HBV+ non-liver cancer patients before anti-PD-1 therapy. We indeed found that HBV+ ESCC patients had a higher fraction of Exhaustionhi T compared with HBV– patients by integrating public data. We further found the HBV+patients had a higher fraction of Exhaustionhi T than HBV– individuals with another 16 samples. Our observation that HBV+ non-liver cancer patients are subject to stronger immunosuppression are reasonable since it is well-known that both HBV infection [14, 42] and cancer [43] can lead immunosuppression and T-cell exhaustion. A higher fraction of Exhaustionhi T in HBV+ non-liver cancer patients may lead to an enhanced response to anti-PD-1 therapy.
Despite the limited sample size of this study, our clinical data provided solid evidence to support the notion that anti-PD-1 therapy is a fairly good treatment for non-liver cancers. Our scRNA-seq data generated from the 3 ESCC patients provided important clues about the changes of immune cells. We compared the exhaustion scores of CD8+ TEFF between HBV+ ESCC patients and HBV– ESCC patients and provided biological insight into their differential. In conclusion, we found anti-PD-1 therapy achieved better outcomes in non-liver cancer patients with HBV infection. We also found HBV+ non-liver patients had a higher fraction of Exhaustionhi T, which could explain the differential response to anti-PD-1 between HBV+ and HBV– non-liver cancer patients.
Methods
Study design and patients
We designed a multicenter retrospective study to evaluate the efficacy of anti-PD-1 therapy on HBV+ non-liver cancer patients by comparing them with their HBV– counterparts (Supplementary Fig. S1A). This study was approved by IRB at the Southwest Hospital, AMU (ID: KY2021112). A total of 84 eligible non-liver cancer patients (35 HBV+ patients and 49 matched HBV– patients) were screened from 7231 cancer patients in clinical database of the Southwest Hospital, Xinqiao Hospital and the Third Affiliated Hospital, CQMU from 2018 to 2021 (Supplementary Fig. S1A).
The 35 HBV+ patients are HBsAg-positive (35; 100%) and HBcAb-positive (33; 94.3%). Furthermore, HBV-DNA was detected in all the 35 HBV+ patients, indicating HBV was active in those patients (Supplementary Table 1). Clinical data was collected for each patient, including age, gender, cancer types, metastasis, medical history, and therapy history. HBV+ non-liver cancer patients and HBV– non-liver cancer patients had matched for tumor types, including lung carcinoma, digestive system carcinoma, nasopharyngeal carcinoma, urogenital neoplasms, melanoma, and esophageal carcinoma. The 49 HBV– non-liver cancer patients were selected to match the 35 HBV+ non-liver cancer patients in regard to cancer types, age, and gender, thus the features of HBV– patients were similar to that of HBV+ patients (Supplementary Fig. S1A, Supplementary Tables 1–3).
Anti-PD-1 therapy
Overall, 26 out of 35 HBV+ non-liver cancer patients received anti-PD-1 combination therapy (anti-PD-1 plus chemotherapy), while the other 9 patients received anti-PD-1 monotherapy. The anti-PD-1 antibodies used in this study include Toripalimab, Sintilimab, Pembrolizumab, Camrelizumab, Nivolumab, and Tislelizumab. Similarly, 39 out of the 49 HBV– non-liver cancer patients received anti-PD-1 combination therapy, while the other 10 patients received anti-PD-1 monotherapy. In addition, some patients were followed up for obtaining survival information after treatment (Supplementary Tables 1–3).
Patients enrolled in this study did not show any immune-related adverse event (irAE) or low grade (1 grade or 2 grade) irAEs during and after anti-PD-1 therapy, except three 3-grade irAEs, namely maculopapular rash in a HBV+ patient, which may be caused by drug treatment [44], cardiotoxicity in a HBV– patient and interstitial pneumonia in a HBV– patient. None of the patients died of irAEs or developed HBV-associated disease (Supplementary Tables 1–3). Overall, irAEs in HBV+ non-liver cancer patients were not significantly different from that in their HBV– counterparts, consistent with recent reports that anti-PD-1 therapy was safe for treatment of HBV+ non-liver cancer patients [17, 18, 40].
scRNA-seq library preparation and sequencing
We performed scRNA-seq of 6 PBMCs samples derived from the three ESCC patients. This study was approved by IRB at the Southwest Hospital, AMU (ID: KY2021112). All individuals signed an informed consent form approved by the IRBs at the time of enrollment. Patients received the treatment according to clinical routine. The scRNA-seq library preparation and data analyses following our previous studies [45,46,47].
Approximately 5 mL blood sample was phlebotomized from each of the three ESCC patients pre- and post-treatment and was transported to the laboratory in an EDTA- containing anticoagulant tube (WEGO, Blood Collection Tube with Vacuum for Single Use, EDTA-K2). Blood samples were centrifuged at 1750 rpm for 8 min to collect the white blood cells. Lymphocyte Separation Medium (Human) (TBD, LTS1077) was added to substratum cell precipitation to obtain PBMCs after centrifugating at 1950rpm for 20 min and washed in 0.04% BSA of PBS (Gibco PBS pH7.4 basic 1x, 8121456) for twice with centrifuging at 1000 rpm for 5 min. Trypan blue (Beyotime Trypan blue 10 mL, C360I-2) staining was used to check for cell viability. The single-cell RNA sequencing libraries of PBMCs were generated using Chromium Next GEM Single-Cell 3’ Reagent Kit v3.1 from 10x genomics. Briefly, 20,000 cells were loaded into each channel of the 10x Chromium controller to generate single-cell GEMs. Single-cell RNA-seq libraries were prepared using the Chromium Single-Cell 3’ Gel Bead and Library Kit (10x Genomics), according to manufacturers’ instructions. Sequencing libraries were loaded at 2.4 pM on an Illumina NovaSeq 6000 with 2× 75 paired-end kits.
Single-cell RNA-seq data pre-process
Briefly, we mapped the reads to the human reference genome GRCh38 and aggregated the gene expression matrixes using Cell Ranger (V.5.0.0). We constituted a Seurat object using the filter expression matrixes from Cell Ranger. Seurat package (V.3.2.3) [48] and Harmony package (V.1.0) [49] were used for the following analyses. Genes expressed in <20 cells were filtered out. Low-quality cells were filtered out if they met the criteria: (1) <200 features number; (2) <800 UMI number; or (3) >15% mitochondrial genome UMIs. We further normalized the gene expression matrix using the NormalizeData function. We identified the top 3,000 highly variable genes using the FindVariableFeatures function and further scaled the gene expression matrix using the ScaleData function. We conducted principal component analysis (PCA) for linear dimension reduction using the RunPCA function and the top 30 PCs were used for further analyses. The RunHarmony function in Harmony package was used to remove the batch effect in our data. Finally, Uniform Manifold Approximation and Projection for Dimension Reduction (UMAP) were conducted using the RunUMAP function.
Clustering and identification of cell subsets
The cell clusters were inferred using FindNeighbors and FindClusters functions based on the normalized gene expression matrix. We identified 17 clusters that were annotated according to their specific expression of classic markers, including 9T cell subsets (CD3D), 2 natural killer cell (NK) subsets (KLRF1), 2 monocyte subsets (FCN1), 2 B-cell subsets (CD19), mDC (CLEC10A) and plasma cell (MZB1) (Supplementary Fig. S2B). The 2 NK cell subsets are NCAM1bright NK and NCAM1dim NK. The 2 B-cell subsets are naive B cells (BN) (TCL1A+) and memory B cells (BM) (CD27+). The 2 monocyte subsets are classic monocyte (cMo) (CD14+) and non-classical Monocyte (ncMo) (FCGR3A+). The 9T cell subsets are CD4+ Naive T cells (CD4+ TN) (CD4+ CCR7+), CD8+ Naive T cells (CD8+ TN) (CD8+ CCR7+), CD4+ memory T cells (CD4+ TM) (CD4+ SELLlo CCR7lo), CD8+ memory T cells (CD8+ TM) (CD8+ SELLlo CCR7lo), CD4+ effector T cells (CD4+ TEFF) (CD4+ GZMA+ GZMB+), CD8+ effector T cells (CD8+ TEFF) (CD8+ GZMA+ GZMB+), CD4+ regulatory T cells (CD4+ Treg) (CD4+ FOXP3+), γδT (TRDC+ TRGC1+) and mucosal-associated invariant T-cell (MAIT) (ZBTB16+ SLC4A10+) (Fig. 2A, B).
Comparison of population size and DEGs
We calculated the PopIndex of each cell subset to measure the relative change of population size after anti-PD-1 therapy. The PopIndex was defined as PopIndex = (Fpost – Fpre)/Fpre, where Fpre is the fraction of a cell subset in PBMCs before anti-PD-1 therapy, and Fpost is the fraction of the same cell subset in PBMCs after anti-PD-1 therapy. A positive value of PopIndex indicates the increase of population size after anti-PD-1 therapy, and a negative value indicates population size decrease after anti-PD-1 therapy.
The FindMarkers function was utilized to find DEGs between the cell subsets. DEGs with bonferroni-adjusted p-values < 0.05 and average Log2 (fold-change) (avg_log2FC) > 0.2 were kept for further analyses.
Comparison of functional features for each cell subset pre- and post-anti-PD-1 therapy
In order to evaluate the changes of functional features of each cell subset pre- and post-anti-PD-1 therapy, we calculated the cytotoxicity score, exhaustion score, major histocompatibility complex (MHC)-I score and MHC-II score of each cell using the AddModuleScore function. All the scores were measured according to expression of specific functional genes. For instance, the cytotoxicity score was calculated based on 12 cytotoxicity-associated genes: PRF1, IFNG, GNLY, NKG7, GZMB, GZMA, GZMH, KLRK1, KLRB1, KLRD1, CTSW, and CST7. The exhaustion score was calculated based on 6 exhaustion-associated genes: LAG3, TIGIT, PDCD1, CTLA4, HAVCR2 and TOX. The MHC-I score was calculated based on 9 MHC-I-associated genes: HLA-A, HLA-B, HLA-C PSMB8, PSMB9, TAP1, TAP2, TAPBP and B2M. The MHC-II score was calculated based on 11 MHC-II-associated genes: HLA-DRA, HLA-DRB1, HLA-DRB5, HLA-DPA1, HLA-DPB1, HLA-DQA1, HLA-DQA2, HLA-DQB1, HLA-DQB2, HLA-DMA, and HLA-DMB.
Statistical analyses
The OR value and p-value of clinical statistical analyses were calculated using Fisher’s exact test. The statistical significance of the score of gene sets was determined by Student’s t-test. Log-rank-test was utilized to calculate the p-value of the survival analyses. Kruskal–Wallis test was used to determine whether the proportion of Exhaustionhi T were significantly different between two samples.
Data availability
The sequencing data have been deposited in Genome Sequence Archive in BIG Data Center under the accession number HRA002492.
Code availability
Code used for analyses is available upon request.
References
El-Khoueiry AB, Sangro B, Yau T, Crocenzi TS, Kudo M, Hsu C, et al. Nivolumab in patients with advanced hepatocellular carcinoma (CheckMate 040): an open-label, non-comparative, phase 1/2 dose escalation and expansion trial. Lancet. 2017;389:2492–502.
Schachter J, Ribas A, Long GV, Arance A, Grob JJ, Mortier L, et al. Pembrolizumab versus ipilimumab for advanced melanoma: final overall survival results of a multicentre, randomised, open-label phase 3 study (KEYNOTE-006). Lancet. 2017;390:1853–62.
Kim ST, Cristescu R, Bass AJ, Kim KM, Odegaard JI, Kim K, et al. Comprehensive molecular characterization of clinical responses to PD-1 inhibition in metastatic gastric cancer. Nat Med. 2018;24:1449–58.
Weber JS, D’Angelo SP, Minor D, Hodi FS, Gutzmer R, Neyns B, et al. Nivolumab versus chemotherapy in patients with advanced melanoma who progressed after anti-CTLA-4 treatment (CheckMate 037): a randomised, controlled, open-label, phase 3 trial. Lancet Oncol. 2015;16:375–84.
Reck M, Rodriguez-Abreu D, Robinson AG, Hui R, Csoszi T, Fulop A, et al. Pembrolizumab versus chemotherapy for PD-L1-positive non-small-cell lung cancer. N Engl J Med. 2016;375:1823–33.
Kudo T, Hamamoto Y, Kato K, Ura T, Kojima T, Tsushima T, et al. Nivolumab treatment for oesophageal squamous-cell carcinoma: an open-label, multicentre, phase 2 trial. Lancet Oncol. 2017;18:631–9.
Matson V, Chervin CS, Gajewski TF. Cancer and the microbiome-influence of the commensal microbiota on cancer, immune responses, and immunotherapy. Gastroenterology. 2021;160:600–13.
Elkrief A, Derosa L, Kroemer G, Zitvogel L, Routy B. The negative impact of antibiotics on outcomes in cancer patients treated with immunotherapy: a new independent prognostic factor? Ann Oncol. 2019;30:1572–9.
Arbour KC, Mezquita L, Long N, Rizvi H, Auclin E, Ni A, et al. Impact of baseline steroids on efficacy of programmed cell death-1 and programmed death-ligand 1 blockade in patients with non-small-cell lung cancer. J Clin Oncol. 2018;36:2872–8.
Park JW, Chen M, Colombo M, Roberts LR, Schwartz M, Chen PJ, et al. Global patterns of hepatocellular carcinoma management from diagnosis to death: the BRIDGE Study. Liver Int. 2015;35:2155–66.
Zhu AX, Finn RS, Edeline J, Cattan S, Ogasawara S, Palmer D, et al. Pembrolizumab in patients with advanced hepatocellular carcinoma previously treated with sorafenib (KEYNOTE-224): a non-randomised, open-label phase 2 trial. Lancet Oncol. 2018;19:940–52.
He MK, Peng C, Zhao Y, Liang RB, Lai ZC, Kan A, et al. Comparison of HBV reactivation between patients with high HBV-DNA and low HBV-DNA loads undergoing PD-1 inhibitor and concurrent antiviral prophylaxis. Cancer Immunol Immunother. 2021;70:3207–16.
Finn RS, Ryoo BY, Merle P, Kudo M, Bouattour M, Lim HY, et al. Pembrolizumab as second-line therapy in patients with advanced hepatocellular carcinoma in KEYNOTE-240: a randomized, double-blind, phase III trial. J Clin Oncol. 2020;38:193–202.
Wykes MN, Lewin SR. Immune checkpoint blockade in infectious diseases. Nat Rev Immunol. 2018;18:91–104.
Ramsey SD, Unger JM, Baker LH, Little RF, Loomba R, Hwang JP, et al. Prevalence of hepatitis B virus, hepatitis C virus, and HIV infection among patients with newly diagnosed cancer from academic and community oncology practices. JAMA Oncol. 2019;5:497–505.
Wen X, Yao W, Ding Y, Li D, Zhang XS. Safety of immune checkpoint inhibitors in Chinese patients with melanoma. Melanoma Res. 2016;26:284.
Tio M, Rai R, Ezeoke OM, Mcquade JL, Zimmer L, Khoo C, et al. Anti-PD-1/PD-L1 immunotherapy in patients with solid organ transplant, HIV or hepatitis B/C infection. Eur J Cancer. 2018;104:137–44.
Kothapalli A, Khattak MA. Safety and efficacy of anti-PD-1 therapy for metastatic melanoma and non-small-cell lung cancer in patients with viral hepatitis: a case series. Melanoma Res. 2018;28:155–8.
Xu F, Zeng Z, Yan B, Fu Y, Sun Y, Yang G, et al. Safety and efficacy of anti-PD-1 inhibitors in Chinese patients with advanced lung cancer and hepatitis B virus infection: a retrospective single-center study. Transl Lung Cancer Res. 2021;10:1819–28.
Zhang X, Tian D, Chen Y, Chen C, He LN, Zhou Y, et al. Association of hepatitis B virus infection status with outcomes of non-small cell lung cancer patients undergoing anti-PD-1/PD-L1 therapy. Transl Lung Cancer Res. 2021;10:3191–202.
Lin Z, Zhang X, Zhou Y, Chen C, He LN, Li H, et al. Hepatotoxicity associated with PD-1 blockade antibodies in cancer patients co-infected with hepatitis B virus. Cancer Immunol Immunother. 2022;71:1247–55.
Zhong L, Zhong P, Liu H, Li Z, Nie Q, Peng W. Hepatitis B virus infection does not affect the clinical outcome of anti-programmed death receptor-1 therapy in advanced solid malignancies: real-world evidence from a retrospective study using propensity score matching. Medicine (Baltimore). 2021;100:e28113.
Zhao Z, Hu R, Chen Y, Zhou G, Yu S, Feng J. Efficacy and safety of PD-1 immune checkpoint inhibitors in locally advanced and advanced non-small-cell lung cancer patients with chronic infection. Oncol Res Treat. 2022;45:366–74.
Cheng YK, Chen P, Chen DW, Lin ZS, Ye SB, Lan P. Comparative safety, efficacy and survival outcome of anti-PD-1 immunotherapy in colorectal cancer patients with vs without hepatitis B virus infection: a multicenter cohort study. Clin Transl Gastroenterol. 2022;13:e00475.
Zhang X, Zhou Y, Chen C, Fang W, Cai X, Zhang X, et al. Hepatitis B virus reactivation in cancer patients with positive Hepatitis B surface antigen undergoing PD-1 inhibition. J Immunother Cancer. 2019;7:322.
Yoo S, Lee D, Shim JH, Kim KM, Lim Y-S, Lee HC, et al. Risk of hepatitis B virus reactivation in patients treated with immunotherapy for anti-cancer treatment. Clin Gastroenterol Hepatol. 2022;20:898–907.
Hughes BGM, Munoz-Couselo E, Mortier L, Bratland Å, Gutzmer R, Roshdy O, et al. Pembrolizumab for locally advanced and recurrent/metastatic cutaneous squamous cell carcinoma (KEYNOTE-629 study): an open-label, nonrandomized, multicenter, phase II trial. Ann Oncol. 2021;32:1276–85.
Sharonov GV, Serebrovskaya EO, Yuzhakova DV, Britanova OV, Chudakov DM. B cells, plasma cells and antibody repertoires in the tumour microenvironment. Nat Rev Immunol. 2020;20:294–307.
Lei X, Lei Y, Li JK, Du WX, Li RG, Yang J, et al. Immune cells within the tumor microenvironment: Biological functions and roles in cancer immunotherapy. Cancer Lett. 2020;470:126–33.
Gajewski TF, Schreiber H, Fu YX. Innate and adaptive immune cells in the tumor microenvironment. Nat Immunol. 2013;14:1014–22.
Zhou Y, Zhou B, Pache L, Chang M, Chanda SK. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10:1523.
Yamauchi T, Hoki T, Oba T, Jain V, Chen H, Attwood K, et al. T-cell CX3CR1 expression as a dynamic blood-based biomarker of response to immune checkpoint inhibitors. Nat Commun. 2021;12:1402.
Yamauchi T, Hoki T, Oba T, Saito H, Attwood K, Sabel MS, et al. CX3CR1-CD8+ T cells are critical in antitumor efficacy but functionally suppressed in the tumor microenvironment. JCI Insight. 2020;5:e133920.
Eferl R, Wagner EF. AP-1: a double-edged sword in tumorigenesis. Nat Rev Cancer. 2003;3:859–68.
Huang D, Chen X, Zeng X, Lao L, Li J, Xing Y, et al. Targeting regulator of G protein signaling 1 in tumor-specific T cells enhances their trafficking to breast cancer. Nat Immunol. 2021;22:865–79.
Seo YD, Jiang X, Sullivan KM, Jalikis FG, Smythe KS, Abbasi A, et al. Mobilization of CD8(+) T cells via CXCR4 blockade facilitates PD-1 checkpoint therapy in human pancreatic cancer. Clin Cancer Res. 2019;25:3934–45.
Tang Z, Kang B, Li C, Chen T, Zhang Z. GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis. Nucleic Acids Res. 2019;47:W556–w60.
Dinh HQ, Pan F, Wang G, Huang QF, Olingy CE, Wu ZY, et al. Integrated single-cell transcriptome analysis reveals heterogeneity of esophageal squamous cell carcinoma microenvironment. Nat Commun. 2021;12:7335.
Zhang C, Li J, Cheng Y, Meng F, Song JW, Fan X, et al. Single-cell RNA sequencing reveals intrahepatic and peripheral immune characteristics related to disease phases in HBV-infected patients. Gut 2023;72:153–67.
Wen X, Wang Y, Ding Y, Li D, Li J, Guo Y, et al. Safety of immune checkpoint inhibitors in Chinese patients with melanoma. Melanoma Res. 2016;26:284–9.
Kim SJ, Hyeon J, Cho I, Ko YH, Kim WS. Comparison of efficacy of pembrolizumab between epstein-barr virus‒positive and ‒negative relapsed or refractory non-Hodgkin lymphomas. Cancer Res Treat. 2019;51:611–22.
Tang LSY, Covert E, Wilson E, Kottilil S. Chronic hepatitis B infection: a review. JAMA 2018;319:1802–13.
Zheng Y, Chen Z, Han Y, Han L, Zou X, Zhou B, et al. Immune suppressive landscape in the human esophageal squamous cell carcinoma microenvironment. Nat Commun. 2020;11:6268.
Chu T, Zhong R, Zhong H, Zhang B, Zhang W, Shi C, et al. Phase 1b study of sintilimab plus anlotinib as first-line therapy in patients with advanced NSCLC. J Thorac Oncol. 2021;16:643–52.
Zhou B, Jin W. Visualization of single cell RNA-Seq data using t-SNE in R. Methods Mol Biol. 2020;2117:159–67.
Qin P, Pang Y, Hou W, Fu R, Zhang Y, Wang X, et al. Integrated decoding hematopoiesis and leukemogenesis using single-cell sequencing and its medical implication. Cell Discov. 2021;7:2.
Wang X, Shen X, Chen S, Liu H, Hong N, Zhong H, et al. Reinvestigation of classic T cell subsets and identification of novel cell subpopulations by single-cell RNA sequencing. J Immunol. 2022;208:396–406.
Stuart T, Butler A, Hoffman P, Hafemeister C, Satija R. Comprehensive integration of single-cell data. Cell. 2019;177:1888–902.e21.
Korsunsky I, Millard N, Fan J, Slowikowski K, Raychaudhuri S. Fast, sensitive and accurate integration of single-cell data with Harmony. Nat Methods. 2019;16:1–8.
Acknowledgements
We thank all members from the Jin lab for the helpful discussion of the project. We acknowledge the assistance of Core Facilities of SUSTech. The computational work was supported by Center for Computational Science and Engineering at SUSTech.
Funding
This study was supported by Chongqing International Institute for Immunology (2020YJC08), National Key R&D Program of China (2021YFF1200900), National Natural Science Foundation of China (81802783), Guangdong Basic and Applied Basic Research Foundation (2023A1515011908), Shenzhen Innovation Committee of Science and Technology (JCYJ20220818100401003, JCYJ20180504170158430), Shenzhen Science and Technology Program (KQTD20180411143432337), Natural Science Foundation of Chongqing Science & Technology Commission (CSTB2022NSCQ-MSX0220).
Author information
Authors and Affiliations
Contributions
YZ and WJ conceived the project. JZ, ST, XX, XH, YZ, and YP collected the clinical samples. JW, XS, JW, and WS performed experiments with help from TZ and YW. GC carried out the computational analysis with help from BZ, HZ, WJ, ZR, YZ, and NH supervised the project. GC and WJ prepared the manuscript, with all authors’ contribution.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Ethics approval and consent to participate
This study was approved by IRB at the Southwest Hospital, AMU (ID: KY2021112). All individuals signed an informed consent form approved by the IRBs at the time of enrollment.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Zhou, J., Chen, G., Wang, J. et al. Anti-PD-1 therapy achieves favorable outcomes in HBV-positive non-liver cancer. Oncogenesis 12, 22 (2023). https://doi.org/10.1038/s41389-023-00468-0
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41389-023-00468-0