Article Text
Abstract
Objective Exhausted hepatitis B virus (HBV)-specific CD8 T cells in chronic HBV infection are broadly heterogeneous. Characterisation of their functional impairment may allow to distinguish patients with different capacity to control infection and reconstitute antiviral function.
Design HBV dextramer+CD8 T cells were analysed ex vivo for coexpression of checkpoint/differentiation markers, transcription factors and cytokines in 35 patients with HLA-A2+chronic hepatitis B (CHB) and in 29 control HBsAg negative CHB patients who seroconverted after NUC treatment or spontaneously. Cytokine production was also evaluated in HBV peptide-stimulated T cell cultures, in the presence or absence of antioxidant, polyphenolic, PD-1/PD-L1 inhibitor and TLR-8 agonist compounds and the effect on HBV-specific responses was further validated on additional 24 HLA-A2 negative CHB patients.
Results Severely exhausted HBV-specific CD8 T cell subsets with high expression of inhibitory receptors, such as PD-1, TOX and CD39, were detected only in a subgroup of chronic viraemic patients. Conversely, a large predominance of functionally more efficient HBV-specific CD8 T cell subsets with lower expression of coinhibitory molecules and better response to in vitro immune modulation, typically detected after resolution of infection, was also observed in a proportion of chronic viraemic HBV patients. Importantly, the same subset of patients who responded more efficiently to in vitro immune modulation identified by HBV-specific CD8 T cell analysis were also identified by staining total CD8 T cells with PD-1, TOX, CD127 and Bcl-2.
Conclusions The possibility to distinguish patient cohorts with different capacity to respond to immune modulatory compounds in vitro by a simple analysis of the phenotypic CD8 T cell exhaustion profile deserves evaluation of its clinical applicability.
- chronic viral hepatitis
- immune response
- cellular immunity
- T lymphocytes
Data availability statement
All data relevant to the study are included in the article or uploaded as online supplemental information.
This is an open access article distributed in accordance with the Creative Commons Attribution 4.0 Unported (CC BY 4.0) license, which permits others to copy, redistribute, remix, transform and build upon this work for any purpose, provided the original work is properly cited, a link to the licence is given, and indication of whether changes were made. See: https://creativecommons.org/licenses/by/4.0/.
Statistics from Altmetric.com
WHAT IS ALREADY KNOWN ON THIS TOPIC
Hepatitis B virus (HBV)-specific CD8 T cells play a crucial role in antiviral protection and failure to control HBV is associated with a severe impairment of virus-specific CD8 T cell responses.
Exhausted HBV-specific CD8 T cells are heterogeneous in chronic HBV infection.
WHAT THIS STUDY ADDS
Simultaneous staining of HBV-specific CD8 T cells with exhaustion and differentiation/memory markers allows to define a phenotypic score for the classification of chronic hepatitis B patients according to the degree of HBV-specific CD8 T cell exhaustion.
Calculation of this score in individual patients can distinguish different cohorts of patients with distinct probability of response to immune modulatory compounds in vitro.
Costaining with PD-1, TOX, CD127 and Bcl-2 gives similar results when applied to the total CD8 T cell population providing a simplified approach to predict individual HBV-specific CD8 T cell functionality without requiring direct analysis of virus-specific CD8 T cells.
HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY
In vitro development of an exhaustion score based on total instead of HBV-specific CD8 T cell phenotyping represents a promising experimental basis in the perspective of novel T cell-based diagnostic predictors.
Introduction
In hepatitis B, chronic evolution of infection is associated with impaired hepatitis B virus (HBV)-specific CD8 T cell responses which are believed to play a key role in the pathogenesis of virus control and persistence. For this reason, a possible therapeutic strategy for chronic HBV infection is to reconstitute an efficent HBV-specific CD8 T cell function.1 2
Available data indicate that exhausted virus-specific CD8 T cells are heterogeneous in their antiviral activity and different studies in mouse chronic lymphocytic choriomeningitis virus infection and in human chronic HCV infection indicate that simultaneous detection of different transcription factors, coinhibitory receptors and differentiation molecules can allow to distinguish T cell subsets with different degree of functional impairement and different propensity to reconstitute their antiviral activity.3–11
Recently, novel information about the heterogeneity of exhausted virus-specific CD8 T cells was also reported in chronic HBV infection.12 13 Even though the HBV-specific T cell population appears globally dysfunctional, different levels of impairment have been reported for T cell subsets targeting different HBV antigen specificities with higher expression of exhaustion markers and lower expansion capacity of polymerase-specific compared with core-specific CD8 T cells.14–17
Based on this background, the aim of our study was to characterise in better detail the functional heterogeneity of HBV-specific and total CD8 T cells in chronic hepatitis B (CHB) patients, to assess whether this can affect individual responses to immune modulatory strategies and to identify CD8 T cell functional and phenotypic parameters to quantify HBV-specific T cell impairment in individual patients in the perspective of novel diagnostic tools based on T cell analysis.
Materials and methods
Patient populations
A total of 376 CHB patients were enrolled and screened for HLA-A2 expression. A total of 160 resulted HLA-A2+ and 60 of them showed detectable frequencies of core18-27 dextramer+CD8 T cells (figure 1A). Seventy-one HLA-A2+patients were tested also with a polymerase-specific dextramer containing the pol sequence 455–463 and 23 of them came out to be positive. All patients who showed the presence of at least one of the two tested CD8 T cell epitopes (64 patients) represent the final HLA-A2+study population, which comprises the following categories of patients (see also table 1 and online supplemental material):
Supplemental material
Thirty-five treatment-naïve patients (CHB).
Sixteen immune subjects spontaneously recovered from chronic HBV infection (spontaneous seroconversion, SS).
Thirteen NUC resolved patients, as shown by HBsAg clearance following NUC therapy (NUC-RES).
Phenotypical analysis of HBV-specific CD8 T cells. (A) Flow chart of the patient cohorts enrolled in the study. Pie plots indicate core18-27-dextramer+ CD8 T cell frequency. (B) PD-1, TOX, CD39, CD127, Bcl-2 and TCF1 MdFI of HBV-core18-27 and influenza-specific CD8 T cells in the indicated study groups and in the negative CD8 population calculated for each specific marker within the total CD8+ T cells (negative control). Box-whisker plots show median values and 5th/95th percentiles; each dot represents a single patient. Statistics by the Kruskal-Wallis with Dunn’s correction test. (C) Representative histogram plots for each parameter. (D) Correlation of exhaustion and differentiation/memory markers by core18-27-specific CD8 T cells in CHB patients. Statistics by the Spearman’s correlation test. CHB, chronic hepatitis B; HBV, hepatitis B virus; MdFI, median fluorescence intensity; NUC RES, NUC resolved patients; SS, spontaneous seroconversion.
Demographic and clinical details of different HLA-A2+ patient categories
Seven healthy subjects (H) were enrolled as controls.
Additional 24 HLA-A2 negative treatment-naïve CHB patients were enrolled and studied to further validate the association between phenotype of total CD8 T cells and T cell responsiveness to modulatory compounds (table 2).
Demographic and clinical details of HLA-A2- CHB patients
Phenotypic analysis of HBV-specific CD8 T cells
To define the stage of CD8+T cell differentiation, dextramer+CD8+ cells were stained with antibodies to cytokine receptors (CD127), anti-apoptosis/cell survival (Bcl-2 and TCF-1) and T cell exhaustion (PD-1, TOX and CD39) markers, as reported in online supplemental methods.
Ex vivo functional assessment of HBV-specific CD8 T cells
Intracellular cytokine staining (ICS) for IFN-γ and TNF-α was performed on dextramer+CD8+ T cells following stimulation with phorbol-12-myristate-13-acetate (PMA) (100 ng/mL) and ionomycin (1 µg/mL) for 4 hours (see online supplemental methods).
‘z-Score’ and ‘Exhaustion index’
Median fluorescence intensity (MdFI) values of exhaustion and memory markers (PD-1, TOX, CD39 and CD127, Bcl-2, TCF-1, respectively) expressed by HBV-specific CD8 T cells were standardised by calculation of the z-Score (details in online supplemental materials) to overcome the problem raised by the different order of magnitude of MdFI values detected for individual parameters. All z-Scores of each patient’s MdFI were calculated relative to reference populations (original set of data from which m and s were derived) represented by NUC-RES and SS. A new variable, named Exhaustion Index (EI), was then created by averaging PD-1, TOX, CD39, CD127, Bcl-2 z-Scores: EI=mean (z-PD-1 + z TOX + z-CD39 – z-CD127 – z-Bcl-2). This new variable takes into account the contribution of each MdFI value irrespective of its original scale. A threshold of 2 (like in the gaussian case) was choosen as cut-off to separate patients with high vs low EI.
In vitro T-cell expansion and treatment with immunomodulators
T cell cytokines (IFN-γ, TNF-α) were tested by intracellular cytokine staining (ICS), as previously described18 on short-term T-cell lines treated with the following compounds (see online supplemental materials):
Trans-Resveratrol.
MitoTEMPO.
Anti-programmed cell death ligand 1 (PD-L1).
Small molecule PD-L1 inhibitor.
Direct sequencing analysis
The HBV core genome was sequenced in 29 patients with chronic HBV infection (see details in online supplemental materials).
Statistical methods
Two-tailed Mann-Whitney U test and Wilcoxon matched-pairs test were used for statistical comparisons. Differences between multiple patient groups were evaluated by the Kruskal-Wallis non-parametric test corrected for pairwise multiple comparisons. The Spearman’s rank-correlation test was applied for correlation analysis. Further details in online supplemental materials.
Results
Phenotypic heterogeneity of HBV-specific CD8 T cells in CHB
HLA-A02-restricted core18-27-specific CD8+T cells were analysed ex vivo in 35 CHB patients before starting antiviral therapy, all with stably elevated or fluctuating viraemia and ALT levels (as typically observed in HBeAg- CHB patients), in 13 NUC treated patients after anti-HBs seroconversion and therapy suspension (NUC-RES) and in 16 subjects who became anti-HBs+following chronic HBV carriage (‘SS’) (table 1 and figure 1A).
Flow cytometry analysis of the exhaustion and memory markers PD-1, TOX, CD39, CD127, Bcl-2 and TCF14 6 7 19 showed higher PD-1, TOX, CD39 and lower Bcl-2, CD127 expression on core18-27-specific CD8 T cells from untreated CHB patients compared with resolved subjects (NUC-RES and SS) and to influenza-specific CD8 T cells from healthy controls. Instead, TCF1 expression was similar in all patient categories (figure 1B,C). A direct correlation was observed in CHB patients between individual exhaustion markers (PD-1 vs TOX, CD39 vs TOX and CD39 vs PD-1) on one hand and between individual differentiation/memory markers (Bcl-2 vs TCF1 and CD127 vs Bcl-2) on the other (figure 1D). Instead, exhaustion and differentiation/memory markers were inversely correlated to each other (figure 1D).
Quantification of the overall level of CD8 T cell exhaustion in each patient by an exhaustion score (EI) that we generated from z-standardised MdFI values of PD-1, TOX, CD39, CD127, Bcl-2, allowed to distinguish two cohorts of CHB patients with EI values lower or greater than 2, comprising 72% (23 patients; figure 2A) and 28% (9 patients; figure 2A, black bars within the yellow square) of the overall CHB cohort, respectively. Combined analysis of the different exhaustion and memory markers is well recapitulated by the radar plot in figure 2B, where blue and orange lines correspond to each z-normalised phenotypic parameter in CHB patients with high and low EI values, respectively. The high EI CHB patient group was also characterised by the presence of a unique CD8 T cell subset composed of PD-1highCD127low/- cells that was poorly represented or totally negative in the remaining group of CHB patients and completely absent in the resolver groups (NUC-RES and SS) (figure 2C,D). As depicted also in figure 2B, a significant difference in the expression of PD-1, TOX, CD39, Bcl-2 and CD127 was detected in core18-27-specific CD8 T cells of CHB patients with high and low EI (figure 2E,F). A statistically significant positive correlation was detected in naive CHB patients by comparing EI with ALT levels but not with viral load or HBsAg values (online supplemental figure S1).
Supplemental material
Classification of CHB patients according to CD8 T cells expression of exhaustion and memory/differentiation markers. (A) EI values were defined by PD-1, CD39, TOX, CD127, Bcl-2 and TCF1 expression on core18-27-specific CD8 T cells for each study cohort. Each bar represents an individual patient. The yellow area indicates patients with EI values above 2. (B) Radar plots depict the mean z-Score values of the indicated HBV-specific CD8 T cell exhaustion and memory/differentiation markers. Blue and orange lines indicate z-Score values of patients with high or low EI values, respectively; the grey line NUC-induced and spontaneous anti-HBs seroconverters. (C) Frequencies of PD-1highCD127low/- and CD127+/PD-1+ cells among the core18-27-specific CD8 T cell population (red and grey bars, respectively). Each bar represents an individual patient. (D) Representative FACS plots from two chronic patients with or without the PD-1highCD127low/- subset. (E) PD-1, TOX, CD39, CD127, Bcl-2 and TCF1 expression on core18-27-specific CD8 T cells of CHB patients with high and low EI values. Box-whisker plots as in figure 1; each dot represents a single patient. Statistics by the Mann-Whitney U test. (F) Representative histogram plots for each parameter in each patient cohort. CHB, chronic hepatitis B; EI, Exhaustion Index; HBV, hepatitis B virus.
In order to extend T cell analysis to another HBV antigen, CD8 T cells from patients with active (n=45) and resolved (n=26) chronic infections were also stained with the pol455-463 dextramer. Pol455-463-specific CD8 T cells were found only in a limited proportion of the HBV-infected chronic patients (11 of 45 tested, 24%) (table 1 and figure 3A). In all patient categories they were significantly less represented than core18-27-specific CD8 T cells (figure 3B,C, p<0.0001 and p=0.0012 in CHB and resolved patients, respectively). In line with previous studies,14–16 in chronic patients pol455-463-specific CD8+T cells expressed higher levels of TOX and lower levels of PD-1 (figure 3D, p=0.0156 and p=0.07, respectively) compared with core18-27-specific CD8+T cells, while no differences in CD39, Bcl-2, TCF1 and CD127 expression were observed between CD8 cells of different specificity (figure 3D). No significant differences were also detected in the expression of exhaustion and differentiation/memory markers between pol455-463-specific CD8+T cells from HBV-infected chronic patients and individuals with spontaneous or NUC-induced control of infection (figure 3E). Overall, our results suggest that pol455-463 and core18-27-specific CD8 T cells in CHB significantly differ with respect to peripheral blood frequencies and phenotypic parameters highlighting two distinct exhaustion profiles.
Different detection rate and phenotype of core18-27-specific and pol455-463-specific CD8 T cells. (A, B) Frequencies of core18-27- and pol455-463-specific CD8 T cells within the total CD8+T cell population from 31 CHB and 18 resolved patients who showed a detectable frequency for at least one HBV dextramer (paired core and polymerase data are illustrated in panel B). (C) Representative FACS plots from two CHB and three resolved patients. (D) PD-1, TOX, CD39, CD127, Bcl-2 and TCF1 expression within core18-27- and pol455-463-specific CD8 T cells from 7 CHB patients; statistics by Wilcoxon matched-pairs test. (E) PD-1, TOX, CD39, CD127, Bcl-2 and TCF1 expression on pol455-463-specific CD8 T cells in CHB and resolved patients (n=10 for both groups). CHB, chronic hepatitis B; HBV, hepatitis B virus.
Different CD8 T cell functionality of CHB patients with different EI values
Cytokine production by core18-27-specific CD8 T cells was then analysed ex vivo (figure 4A). The same analysis on pol-specific CD8 T cells was precluded by their extremely low frequency. Functionality of core-specific CD8 T cell was very variable, because some CHB patients showed cytokine production levels comparable to NUC and spontaneous anti-HBs seroconverters, while other CHB patients produced cytokines less efficiently (figure 4B, left graph).
Ex vivo HBV-specific T cell functional analysis. Cytokine production by HBV-specific and influenza-specific CD8 T cells: (A) representative plots; (B) different patient categories assessed after PBMC stimulation with PMA and Ionomycin (box whisker plots, on the left); CHB patients split according to EI values (right graph). (C) Correlation between cytokine production by the core18-27-specific CD8 T cells and EI in chronic patients with high (orange dots) and low (blue dots) EI, respectively, and in HBV resolved patients (black dots). (D, E) Cytokine production by PD-1highCD127low/- and PD-1+CD127+ HBV-specific CD8 T cells from CHB patients and from PD-1+CD127+ CD8 T cells from the different patient categories. (F) Correlation between cytokine production and Bcl-2 (left graph) or TOX (right graph) expression by core18-27-specific CD8 T cells. Statistics by the Kruskal-Wallis with Dunn’s correction test (B, E) and the Spearman’s correlation test (C, F). CHB, chronic hepatitis B; EI, Exhaustion Index; HBV, hepatitis B virus; NUC RES, NUC resolved patients; PBMC, peripheral blood mononuclear cell; SS, spontaneous seroconversion.
HBV-specific CD8 T cells from CHB patients with lower EI were more efficient in cytokine production (figure 4B, right graph) which was inversely correlated to the EI (figure 4C). PD-1highCD127low/- CD8 T cells were functionally less efficient than PD-1+CD127+ cells (p=0.0113, figure 4D). A functional comparison of the PD-1highCD127low/- CD8 T cell subset in the different patient cohorts was not possible because this subset was selectively present only in CHB patients. Instead, PD-1+CD127+ CD8 T cells were detectable in all patient cohorts and their cytokine production was variable in chronic patients, some of whom were less efficient than NUC and spontaneous resolvers (figure 4E).
Total cytokine production by the overall HBV-specific CD8 T cell population showed a positive and negative correlation with Bcl-2 and TOX levels, respectively (figure 4F). The same trend (but without significant levels of correlation) was observed with the other exhaustion (PD-1 and CD39) and memory (CD127, TCF-1) markers.
CD8 T cell phenotypic profile and response to immune modulatory interventions in vitro
We then assessed whether HBV-specific CD8 T cells of CHB patients with different EI levels, were more or less sensitive to the effect in vitro of antioxidant, polyphenolic and PD-1/PD-L1 targeting compounds18 20–26 and a selective toll-like receptor 8 (TLR8) agonist.27 28
In the absence of these immune modulatory compounds, T-cell lines produced in CHB patients with low EI showed better expansion and cytokine production compared with T cell lines from patients with high EI (figure 5A,B), in keeping with ex vivo results.
Effect of immune modulatory interventions on the HBV-specific CD8 T cell function. Expansion capacity calculated as the ratio (fold-change) between in vitro and ex vivo frequencies of core18-27 dextramer+ CD8 T cells (A) and percentage of cytokine-positive CD8 T cells (B) in short-term T-cell lines generated by core18-27 peptide stimulation of PBMC from CHB patients with high (n=7) and low (n=13) EI; statistics by Mann-Whitney U test. (C) Percentage of double-positive IFNγ+TNFα+ CD8 T cells in paired short-term T-cell lines generated as in (B) in the presence (treated) or absence (untreated) of Resveratrol (RSV), MitoTempo (MT), anti-PD-L1, a small PD-L1 inhibitor molecule (PD-L1 SM) and a TLR8 agonist (TLR8a) from chronic patients with high (n=7) and low (n=13) EI; statistics by the Wilcoxon-matched-paired test. (D) Delta values of double-positive IFNγ+TNFα+ CD8 T cells derived by subtracting CD8 T cell frequencies of untreated from treated short-term T-cell lines (black lines indicate the median values; statistics by Mann-Whitney U test). (E) Representative examples of cytokine production in two CHB patients. (F) Hierarchical-clustering of HBV-specific CD8 T-cell responses induced by in vitro core18-27 stimulation in presence of the different treatments in chronic patients with high (n=7 orange) and low (n=13 blue) EI. CHB, chronic hepatitis B; EI, Exhaustion Index; HBV, hepatitis B virus; N/A, not available; PBMC, peripheral blood mononuclear cell.
All tested compounds had significantly greater effect on HBV-specific CD8 T cells of patients with low EI (figure 5C–E, online supplemental figures S2, S3). Functional changes induced by these immune modulatory compounds are well recapitulated by hierarchical-clustering analysis, which demonstrates a clear segregation of upregulated and downregulated functional parameters in chronic patients with low and high EI (figure 5F). Interestingly, the anti-PD-L1 monoclonal antibody and the PD-L1 small molecule inhibitor24 displayed a common clustering indicating similar quality and intensity of CD8 T cell modulation. Likewise, the effects of Resveratrol and the TLR8-agonist clustered together, suggesting a comparable pattern of modulatory activity (figure 5F). By direct comparison of the different compounds, the TLR8-agonist was the most effective in the reconstitution of the CD8 T cell function (online supplemental figure S4).
Supplemental material
Supplemental material
Supplemental material
Prediction of immune T cell reconstitution by phenotypic T cell analysis of the whole CD8 T cell population
Evidence of CD8 T cell rejuvenation by immune modulatory treatment in vitro and prediction of response to immune modulation are based on the analysis of a single antigenic specificity and may thus be partially relevant with respect to the multispecific T cell response primed in vivo by infection. Interestingly, we noted a significant positive correlation between core18-27-specific and total CD8 T cells relative to the frequency of PD-1highCD127low/- cells (figure 6A) and to the expression of differentiation/exhaustion markers, such as Bcl-2 and CD39, and transcription factors, such as TCF-1 and TOX (figure 6B). This evidence prompted us to assess whether the identification of patients most likely to respond to immunomodulatory compounds in vitro can be accomplished also by phenotypic analysis of total CD8 T cells.
Phenotypic analysis of the total CD8 T cell population. (A) Correlation between logarithmic frequency values of PD-1highCD127low/- T cells among core18-27-specific and total CD8 T cells. (B) Correlations of TCF1, Bcl-2, TOX and CD39 expression in total and HBV-specific CD8 T cells. (C) Representative plots of PD-1high/CD127low/- and TOX+/PD-1high CD8 T cell subsets (yellow boxes) in two CHB patients with high and low EI. (D) Top: phenotypic profiles of total CD8 T cells of individual patients (grey bars) in relation to the corresponding EI values derived from HBV-specific CD8 T cell analysis (black bars, on the left). Red areas indicate the significance threshold for each phenotypic CD8 T cell subpopulation calculated by using a ROC curve analysis based on exhaustion index classification criteria (details in online supplemental materials). In the internal squares, each dot represents the frequency of total CD8 T cells expressing the indicated phenotypic profile in chronic patients with high or low EI (Mann-Whitney test). Bottom: correlation between frequency of cytokine positive core18-27-specific CD8 T cells ex vivo and frequency of PD-1high/TOXhigh, PD-1high/TOXhigh/ Bcl-2-, PD-1high/TOXhigh/CD127- or PD-1high/TOXhigh/TCF1- among total CD8 T cells in CHB patients. Statistics by the Spearman’s correlation test (A, B, D). CHB, chronic hepatitis B; EI, Exhaustion Index; HBV, hepatitis B virus; ROC, receiver operating characteristic.
Remarkably, the same distribution of patients derived from the calculation of the EI in HBV-specific CD8 T cells, was also found by analysing the expression of different combinations of PD-1, CD127, Bcl-2 and TOX on total CD8 T cells (figure 6C,D). In particular, the frequency of PD-1highTOXhigh, PD-1highTOXhighBcl-2-, PD-1highTOXhighTCF1- and PD-1highTOXhighCD127- subsets in total CD8 T cells was significantly greater in CHB patients with high EI (details in online supplemental results and figure S5). These observations indicate that total CD8 T cell phenotyping can reflect HBV-specific CD8 T cell analysis for the identification of patient cohorts with different antiviral function and different in vitro responsiveness to modulatory compounds. This conclusion is further supported by the significant inverse correlation between cytokine production by HBV-specific T cells and the phenotypic profile of total CD8 T cells (figure 6D, bottom charts) and by the significant inverse correlation between HBV-specific T cell responses to the distinct immune modulations and the frequency of PD-1highTOXhigh, PD-1highTOXhighBcl-2-, PD-1highTOXhighCD127- and PD1highTOXhighTCF1- CD8 T cell subsets (online supplemental figure S6). Instead, significant correlations between phenotypes of the bulk CD8+T cell population and clinical/virological parameters (HBV-DNA, ALT, HBsAg) were not observed in treatment naïve CHB patients, with the exception of a positive trend observed only for ALT levels (online supplemental figure S7).
Supplemental material
Supplemental material
Supplemental material
CD8 T cell phenotype can predict recovery of multi-specific T cell responses to immune modulation. Delta values of IFNγ+, TNFα+ and double-positive IFNγ+TNFα+ CD8 T cells derived by subtracting CD8 T cell frequencies of untreated from treated short-term T-cell lines generated by stimulation of PBMC from CHB patients (n=32) with 15-mer overlapping HBV core and polymerase peptide pools. Data are shown as sum of T cell responses against the two antigens. (A) Segregation of individual CD8 T cell responses to core and polymerase peptides based on the PD-1highTOXhigh CD8 T cell frequency threshold (0.4% value) previously obtained using an ROC curve analysis as illustrated in figure 6D. Statistics by the Mann-Whitney U test. (B–D) Correlation between PD-1highTOXhigh CD8 T cell subset frequency and percentage of IFNγ+, TNFα+ and double-positive IFNγ+TNFα+ CD8 T cells in short-term T-cell lines generated as described above, illustrated for each immune modulation (B), as cumulative data of all treatments (C) and as numbers of immune modulatory agents able to induce a positive response (D). Statistics by the Spearman’s correlation test. (E) Representative examples of cytokine production in a CHB patient with low PD-1highTOXhigh CD8 T cell frequency. CHB, chronic hepatitis B; HBV, hepatitis B virus; PBMC, peripheral blood mononuclear cell.
We then attempted to validate the T cell functional restoration effect of the immune modulatory compounds and the predictive value of the CD8 T cell-based phenotypic score in the context of a multispecific CD8 T cell response stimulated with broader peptide panels covering the entire core and polymerase sequences, which obviously reflect more closely the complexity of the immune response induced by infection.
The validation cohort was composed of additional 32 newly enrolled treatment naïve CHB patients (24 and 6 HLA-A2 negative and positive, respectively, details in tables 1 and 2). Levels of cytokine production were assessed upon PBMC stimulation with 15-mer overlapping peptide pools covering HBV core and polymerase antigens in the presence or absence of the different immune modulatory agents. A frequency threshold of 0.4% PD-1highTOXhigh CD8 T cells among the overall CD8 T cell population, defined by an ROC curve analysis of phenotypic data generated in HLA-A2 positive patients on core18-27-specific and total CD8 T cells and shown to distinguish patients with different levels of EI values, was applied to analyse the new data set. Results confirm that CHB patients with a PD-1highTOXhigh CD8 T cell frequency below 0.4% show a better response to immune modulatory compounds also when CD8 T cells are stimulated with peptide pools corresponding to the overall core and polymerase sequences able to stimulate a multispecific T cell response (figure 7A). In line with this observation, a statistically significant inverse correlation was observed between cytokine production by core and polymerase stimulated CD8 T cells and PD-1highTOXhigh CD8 T cell subset frequency (figure 7B–E). Although both core-specific and polymerase-specific CD8 T cell responses were improved by in vitro modulation, a slightly greater contribution to the overall HBV-specific T-cell restoration was given by core-specific responses (online supplemental figure S8).
Supplemental material
Impact of mutations on core18-27-specific CD8 T cell responses. (A) Ex vivo frequency of wild-type-specific or variant-specific core18-27 CD8 T cells by staining with the corresponding wild-type of mutated peptide HLA class I dextramers from seven CHB patients; next to each dot, the CHB patient ID number is indicated; statistics by Wilcoxon matched-pairs test. (B) Ratio (fold-change) between the MdFI values of the indicated phenotypic markers expressed by CD8 T cells targeting the mutated or the wild-type core18-27 epitopes (left graph) and EI values defined by PD-1, CD39, TOX, CD127, Bcl-2 and TCF1 expression (right graph). (C) Percentage of cytokine-positive CD8 T cells in short-term T-cell lines generated by WT or mutated core18-27 peptide stimulation of PBMCs from seven CHB patients. The column chart represents median values (upper), while the dot graph shows individual patient responses to wild-type and mutated peptides (lower); statistics by the Wilcoxon-signed-rank test. (D) Dot-plots of ex vivo frequencies of wild-type- or mutated-specific core18-27 CD8 T cells in CHB patients 2, 5, 1. On the right, the corresponding graphs illustrating frequency, phenotype (EI) and function are shown. CHB, chronic hepatitis B; EI, Exhaustion Index; MdFI, median fluorescence intensity; PBMC, peripheral blood mononuclear cell.
These experiments confirm the results obtained with HBV-specific CD8 T cells of a single specificity and validate in the context of a multispecific CD8 T cell response the feasibility of distinguishing patients more prone to respond in vitro to immune modulation by a simple phenotypic assessment of CD8 exhaustion parameters.
Impact of viral sequence variation on phenotype and function of core18-27-specific CD8 T cells
Since HBV is a highly variable virus we investigated whether mutations within the core18-27 sequence can contribute to the functional heterogeneity of HBV-specific CD8 T cell responses. The HBV core genome was sequenced in 29 patients with chronic HBV infection (table 1). Notably, a wide sequence variation with a high mutation rate was detected in the HLA-A2 restricted core18-27 epitope in CHB viraemic patients, because mutated core sequences were detected in 73% of the analysed epitopes (21 patients out of 29). The effect of mutations was quite variable, although frequency, phenotype and function of core18-27-specific CD8 T cells did not change significantly in most of the patients. For example, only CHB 2 and 5 showed more evident frequency variations in opposite directions (figure 8A). However, their phenotypic profile did not change significantly because the EI increase in patient 5 did not reach a value greater than 2 (which is the threshold discriminating high and low levels of exhaustion) (figure 8B) and the decline in IFN-γ production was not associated with a parallel and comparable decline in TNF-α production (figure 8C,D). Similarly, in patient 2 the slight decline in SI (too modest to induce a decline below the threshold of 2) (figure 8B) was associated with a very modest increase in IFN-γ and TNF-α production (figure 8C,D).
Discussion
In chronic HBV infection HBV-specific CD8 T cells appear to be functionally heterogeneous.2 12 13 16 29 To elucidate further the phenotypic and functional features of this CD8 T cell heterogeneity and to identify CD8 T cell-based predictors of response to immune reconstitution therapies we first selected a homogeneous chronic active hepatitis patient population with persistent viraemia and liver inflammation, eligible for therapy and expected to harbour terminally exhausted HBV-specific CD8 T cells. This is particularly relevant because part of the data so far generated in this area are derived from heterogeneous groups of patients, frequently including also patients with an inactive infection. Particular attention was also paid to identify the rare patient populations able to achieve functional cure of chronic infection either spontaneously or after NUC therapy, which are essential to understand whether reacquisition of a fully protective antiviral activity is achievable after long-term exposure to high antigen and viral loads. Despite low frequency of HBV-specific CD8 T cells, an in-depth phenotypic and functional analysis was performed directly ex vivo.
Second, we studied the effect in vitro of antioxidant, polyphenolic, TLR8 agonist and PD-1/PD-L1 blocking compounds on phenotype and function of HBV-specific CD8 T cells to select optimal candidate targets for immune modulatory therapies. Finally, we searched for predictors of in vitro response to immune modulation, in the perspective of a possible clinical translation aimed at personalising immune therapies for CHB.
By combining the analysis of different exhaustion (PD-1, CD39, TOX) and memory (CD127, Bcl-2, TCF1) markers we identified a specific CD8 T cell phenotypic score, that we named EI, suitable to distinguish cohorts of CHB patients with different CD8 T cell exhaustion levels. By focusing on PD-1 and CD127 staining, we found that a subset of PD-1+CD127+ HBV-specific CD8 T cells was present in all patient categories, irrespective of the clinical condition and level of virus control. This is consistent with the detection of TCF1+CD127+PD-1+ memory like CD8 T cells in chronic HCV infection, that persist after therapy-induced HCV clearance and show higher levels of PD-1 and Eomes expression, as well as reduced cytokine production compared with conventional memory HCV-specific T cells generated after spontaneous resolution of HCV infection.4 30 31 Instead, PD-1hiCD127low/- HBV-specific CD8 T cells were present in a more limited proportion of untreated chronic viraemic patients with higher EI levels.
Higher cytokine production was detected in CHB patients with a less exhausted phenotype, marked by low EI levels and the total lack of PD-1hiCD127low/- HBV-specific CD8 T cells. Instead, CHB patients with higher EI levels and a greater representation of PD-1hiCD127low/- CD8 T cells appeared to be significantly less efficient in cytokine production. All together, our observations suggest that the phenotype-based EI score actually reflects the CD8 T cell function and may thus represent a composite multiparametric score for quantification of CD8 T cell exhaustion.
Interestingly, a functional exhaustion score (FES) was also developed by Bengsch et al with the same predictive purpose in HIV and human cancer settings, based on the assessment of exhaustion and effector/memory functional parameters. Although Bengsch’s FES score and our EI score are not necessarily expected to be concordant because of the different clinical settings of analysis and application, their agreement would strongly support the robustness our phenotype-based EI. Preliminary analysis of IFN-γ and TNF-α production profiles (the only components of the FES score available in our dataset) does not seem, however, to show a different behaviour in patients with high or low EI, as predictable if the two scores were concordant.32
To the interindividual CD8 T cell heterogeneity may contribute the presence of mutations within the core18-27 epitope. Thus, the core 18–27 region was sequenced in 29 CHB patients and the effect of mutations on CD8 responses was assessed in 7 of them who harboured variant core18-27 epitopes. HBV mutations caused variable effects on HBV-specific CD8 T cells. Most of them were neutral with no effects on CD8 frequency, phenotype and function, but a few mutations caused increase or decline in CD8 T cell frequencies. These changes, however, had only a marginal impact on phenotypic and functional CD8 T cell profiles because distribution of the patients into high or low EI classes was not modified by the results obtained with mutated instead of prototype sequences. Additional studies are however warranted to drive definitive conclusions because only seven patients with mutated epitopes were analysed in this study.
The possibility to rank the level of CD8 T cell exhaustion of single chronic HBV patients may be extremely helpful for the identification of patients with immunological conditions still amenable of functional reconstitution, for example, with antioxidant and polyphenolic compounds, such as MitoQ and MitoTempo or Resveratrol and Oleuropein, that have recently been proposed by us as promising candidates for CD8 T cell-based therapeutic strategies.18 20 To investigate this issue, we first tested the capacity of different metabolic and immune modulatory compounds to improve HBV-specific CD8 T cell responses in our cohort of chronic viraemic patients with well characterised levels of CD8 T cell exhaustion, and, second, we looked for specific predictive parameters of response to immune modulation. In addition to mitochondrial antioxidant and polyphenolic compounds, HBV-specific CD8 T cells were treated in culture with the TLR8 agonist Selgantolimod and a novel small PD-L1 inhibitor molecule, which represents a promising alternative to monoclonal antibodies.24 All these compounds induced a significant improvement of CD8 T cell functions, with a slightly better effect of the TLR8 agonist over the other compounds. The effect of the different modulatory agents was greater in low EI patients, namely those with predominance of less exhausted T cell phenotypes.
Although very promising, our results suffer, however, from a major flaw, since our analysis was based on a single epitope specificity. To overcome this drawback, first we tried to select and characterise also polymerase-specific CD8 T cells though they were detected only at very low frequencies and in a limited number of viraemic CHB patients. A simultaneous phenotypic profiling of pol455-463- and core18-27-specific CD8 T cells was performed in some CHB patients and, in line with previous studies,14–16 we observed higher levels of TOX and lower levels of PD-1 in pol455-463-specific CD8+T cells. These findings highlight intrinsic differences in exhaustion profiles between pol455-463 and core18-27-specific CD8 T cells in patients with chronic active hepatitis, making the EI score not applicable to pol455-463-specific CD8 T cells.
As a second strategy to overcome the single specificity issue, we then investigated whether the subset distribution detected among core-specific CD8 T cells can actually reflect the features of the overall HBV multispecific CD8 T cell population. To do this, we analysed the phenotypic profile of the total CD8 T cell population, assuming that expression of exhaustion markers in the overall CD8 T cell population should be mostly contributed by the multispecific repertoire of the different HBV antigenic specificities. Remarkably, the same stratification of chronic patients with different levels of functional T cell impairement, as defined by the calculation of the EI score for HBV core-specific CD8 T cells, was reproduced by staining total CD8 T cells with PD-1, CD127, TOX and Bcl-2. This means that a simple and easy to be performed phenotypic assessment of the overall CD8 T cell population may allow to predict individual patient’s responsiveness to in vitro modulation without the need of virus-specific T cell analysis.
In order to validate the predictive value of the total CD8 T cell phenotypic profile on the response to immune modulation, we then used a different cohort of patients and a different experimental approach looking at the effect of immune modulatory compounds on multispecific CD8 responses induced in vitro by PBMC stimulation with overlapping peptides covering the entire core and polymerase proteins. Results indicate that the phenotypic cut-off value derived from the initial series of experiments with core18-27-specific CD8 T cells can be applied to total CD8 T cells from a randomly selected cohort of CHB patients to identify infected hosts who are more or less likely to respond in vitro to immune modulatory compounds. Moreover, these experiments show that CD8 T cell functional recovery by immune modulation is effective also in the context of multispecific T cell responses stimulated by entire HBV proteins. Indeed, functional improvement was observed among both pol-specific and core-specific T cells with a major contribution of the latter to the overall recovery of HBV-specific T-cell responses.
In conclusion, our study demonstrates that co-expression profiling of specific exhaustion and memory markers by HBV core-specific and total CD8 T cells can allow to generate a predictive score suitable for the identification of distinct CHB patient cohorts with various degrees of CD8 T cell dysfunction and different susceptibility to functional T cell restoration in vitro by immune modulatory treatments. A crucial limitation, however, is that the predictive value of the EI score has been tested only in vitro and that clinical studies have not been performed so far to understand whether this simple phenotypic score can have an application in the clinical setting. This is particularly relevant because without a clinical assessment we cannot know whether the in vitro effect of the tested compounds can be reproduced in vivo and whether the level of restoration detected in vitro can actually translate into a clinically meaningful functional improvement. In favour of our prediction score based on a phenotypic profiling of total CD8 T cells is its simplicity of application which makes it a promising basis for the development of clinically useful T cell-based predictors. Indeed, a clinical trial with mitochondrial antioxidant compounds is being organised in CHB patients to characterise their effect in vivo and to validate the clinical value of our phenotypic CD8 T cell score.
Data availability statement
All data relevant to the study are included in the article or uploaded as online supplemental information.
Ethics statements
Patient consent for publication
Ethics approval
The study was approved by the ethics committee of the coordinating centre, Comitato Etico per Parma (0047243_20Dec2017), and the local ethics committees of all collaborating sites accepted the favorable opinion.
References
Supplementary materials
Supplementary Data
This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.
Footnotes
MR and AV contributed equally.
Contributors MR and AV: execution of experiments, acquisition of data, statistical analysis, analysis and interpretation of data. CT, VB, PF, AP and IM: execution of experiments, acquisition and interpretation of data. AL, AA, FB, MM and GA: recruitment and characterisation of the patients. DL, RS, MB and CB: administrative support. GP: statistical analysis. TP and CM: sequencing analysis. GM and PL: interpretation of data, critical revision of the manuscript. SD, SPF, AG, JM, MG and LL: reagents provision, critical revision of the manuscript. CF: study concept and design, critical revision and editing of the manuscript, obtained funding, study supervision and interpretation of data. CB: study concept and design, execution of experiments, acquisition of data, analysis and interpretation of data, writing of the manuscript and guarantor.
Funding This work was supported by a grant from Emilia-Romagna Region, Italy (Programma di Ricerca Regione-Università 2010-2012; PRUa1RI-2012-006), by a grant from the italian Ministry of Health (Ricerca Finalizzata RF 2013-02359333), by a grant from Gilead Sciences Europe (grant ID number 01553 to CF), by a PRIN project from the italian Ministry of the University and Research (protocol code n. 2017MPCWPY), and by a grant from the European Union's Horizon 2020 research and innovation programme (grant agreement no. 848223).
Competing interests CF, Grant: Gilead, Abbvie. Consultant: Gilead, Abbvie, Vir Biotechnology, Arrowhead, Transgene, BMS; PL: advisor and speaker bureau for Gilead, Roche, BMS, GSK, MSD, Arrowhead, Alnylam, Spring Bank, Janssen, EIGER, Myr Pharma; AL: Consultant: MyrPharma, Gilead; MM: advisory board for Abbvie. SD, SPF, AG, JM, MG and LL are employees of and stock-holders in Gilead Sciences. The remaining authors disclose no conflicts.
Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.
Provenance and peer review Not commissioned; externally peer reviewed.
Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.