Article Text

Original research
Helicobacter pylori infection disturbs the tumor immune microenvironment and is associated with a discrepant prognosis in gastric de novo diffuse large B-cell lymphoma
  1. Yuwei Deng1,
  2. Wenjia Su2,
  3. Junwen Zhu1,
  4. Hongfei Ji3,
  5. Xiaoping Zhou1,
  6. Jingshu Geng4,
  7. Jiayu Zhu1 and
  8. Qingyuan Zhang1
  1. 1Department of Medical Oncology, Harbin Medical University Cancer Hospital, Harbin, Heilongjiang, People's Republic of China
  2. 2Department of Hematology, First Affiliated Hospital of Harbin Medical University, Harbin, Heilongjiang, China
  3. 3Institute of Cancer Prevention and Treatment, Heilongjiang Academy of Medical Sciences, Harbin Medical University, Harbin, Heilongjiang, People's Republic of China
  4. 4Department of Pathology, Harbin Medical University Cancer Hospital, Harbin, Heilongjiang, People's Republic of China
  1. Correspondence to Dr Qingyuan Zhang; 0566{at}hrbmu.edu.cn

Abstract

Background Gastric diffuse large B-cell lymphoma (gDLBCL) related to Helicobacter pylori infection exhibits a wide spectrum of prognosis, and the tumor immune microenvironment (TIME) affects tumor progression. However, there are few studies on the correlation between prognosis and changes of TIME induced by H. pylori infection in de novo gDLBCL.

Methods A retrospective study was performed to determine the prognostic value of TIME related to H. pylori infection in de novo gDLBCL. A total of 252 patients were included and have been treated with standard rituximab to cyclophosphamide, doxorubicin, vincristine, and prednisone chemotherapy or other similar regimens in addition to H. pylori eradication (HPE). All patients were stratified by H. pylori infection, HPE efficacy, and preliminary TIME evaluation using conventional criteria. Statistical analyses were conducted. To assess the mechanism, 30 subjects were assessed for H. pylori infection. The components and spatial distributions of TIME were analyzed.

Results The median follow-up of the 252 patients was 66.6 months (range 0.7–119.2), and the 5-year overall survival (OS) was 78.0%. A total of 109 H. pylori-positive cases with pathological complete remission and high tumor-infiltrating T lymphocytes (cohort 1) had significantly higher 5-year progression-free survival (88.1% vs 70.5%, p<0.001) and OS (89.2% vs 76.6%, p<0.001) than the other 143 patients (cohort 2). Among 30 patients, 19 were cytotoxin-associated gene A-marked as the cohort 1 subset. Compared with cohort 2, cohort 1 exhibited increased inflammatory factors (tumor necrosis factor-α, interferon γ, etc) and decreased immunosuppressive components (PD-L1, PD-1, IL-10, etc). There was reduced NF-kB activation. Cancer-promoting immune cells (PD-1hiTim-3+ CTL, Tregs, M2-like macrophages, etc) occupied a minor spatial distribution, while the antitumor subtypes increased, corresponding to favorable survival.

Conclusion H. pylori-evoked inflammatory responses disturb the TIME, causing a differential prognosis in de novo gDLBCL, which can be used to identify patients who could benefit from HPE and immunochemotherapy.

  • immune reconstitution
  • hematologic neoplasms
  • lymphocytes
  • tumor-infiltrating
  • tumor microenvironment
  • inflammation

Data availability statement

Data are available on reasonable request.

http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

Introduction

Diffuse large B cell lymphoma (DLBCL), which is the leading type of aggressive non-Hodgkin’s lymphoma from activated B-cell (ABC) and germinal center B-cell (GCB) origins, arises from extranodal organs in about 30%–45% of cases. The stomach is the most common site, whereby DLBCL accounts for <5% of all gastric malignancies.1 2 Primary gastric DLBCL (de novo gDLBCL) is a heterogeneous disease lacking mucosa-associated lymphoid tissue (MALT) characteristics, which constitutes approximately 30% of all gastric lymphomas.3 In China, the overall prevalence of Helicobacter pylori infection in the gastrointestinal tract is 54%, which is much higher than in other countries in which the figure stands at a maximum of 31%.4 Several studies have established that H. pylori affects biological processes and the prognosis of gastric MALT lymphoma. H. pylori eradication (HPE) applied as the standard first-line treatment of MALT, results in 65% pathological complete remission (pCR), with a median of 7.7 years.5 Compared with Westerners, Asians exhibit a more favorable response rate (84% vs 74%).6

At present, HPE combined with rituximab plus anthracycline-based immunochemotherapy (R-cyclophosphamide, doxorubicin, vincristine, and prednisone, CHOP) has reached clinical consensus due to its marked pCR and favorable overall survival (OS) in gDLBCL, with 50% pCR in early-stage patients.7 However, the initial failure, defined as local relapse, was found to be the most common pattern with systemic therapy (R-CHOP) alone, accounting for 44%–63% of all relapse cases in Asians.8 Although several studies have confirmed the correlation between H. pylori infection and prognosis in de novo gDLBCL,9 a worrying fluctuation of 27%–87% pCR has also been observed among the H. pylori-positive subgroup adopting HPE alone in the Chinese population,10 causing a delay of available immunotherapy due to a lack of pre-evaluation.

Latest studies have emphasized the reliable prognostic stratification based on immune cytotypes related to the tumor immune microenvironment (TIME), demonstrating that ‘Stromal-1’ and ‘Stromal-2’, respectively, predicted favorable and inferior outcomes in DLBCL patients with R-CHOP, regardless of the cell-of-origin.11 Numerous studies have suggested that tumor programmed death ligand-1 (PD-L1) expansion and nuclear factor kappa beta (NF-κB) growth factor signaling modulates programmed death protein-1 (PD-1)/PD-L1 immunosuppression and tumor progression, leading to an inferior prognosis in DLBCL.9 12 A parallel study has shown that H. pylori-positive de novo gDLBCL presenting pCR after HPE indicates a favorable prognosis, especially in the limited-stage subgroup.13 H. pylori infiltration can trigger effective inflammatory immune responses involving the recruitment and activation of dendritic cells, CD4+ T cells (Th1, Th17, and Tregs), neutrophils, and macrophages.14 However, there are few reports on the effect of TIME on the prognosis.

Here, we investigated the intriguing prognostic correlation by evaluating H. pylori infection with pCR after HPE treatment and the corresponding initial TIME status in de novo gDLBCL patients (especially early stage) in a major cohort. The results highlight that H. pylori-induced inflammatory disturbances in the TIME mapping to a wide distribution of phenotypically diverse immune cells may accurately reflect the prognostic differences in H. pylori-related de novo gDLBCL, which was defined by H. pylori-inherent cytotoxin-associated gene A (CagA) toxic factor.15

Methods

Patients

A total of 252 paraffin-embedded cancer specimens collected between 2005 and 2019 were retrieved from the Harbin Medical University Cancer Hospital and the First Affiliated Hospital of Harbin Medical University. All participants provided signed informed consent. Endoscopic biopsy and surgical resection were performed for the acquisition of pathological specimens, which were re-examined by experienced histopathologists to confirm the diagnosis according to the WHO hematological malignant tumor classification system. De novo gDLBCL was defined as excluding histological features of MALT lymphoma. Tumors of MALT-origin proof and secondary gastric lymphoma were excluded. Six cycles of standardized CHOP (cyclophosphamide, doxorubicin, vincristine, and prednisone) or CHOP-like chemotherapy accompanied by rituximab (R) were administered to all the patients. Second-line chemotherapy or other treatments were also available based on indications. The total baseline clinical features of the patients according to medical records are shown in table 1.

Table 1

Clinicopathological features of gastric de novo diffuse large B-cell lymphoma with and without Helicobacter pylori infection

Helicobacter pylori infection

At least one of the following tests was used to confirm the clinical presence of H. pylori infection: histological reports (Warthin–Starry staining was conducted in all cases), bacterial culture, and serological examination. HPE treatment was administered to all H. pylori-positive patients, comprizing bismuth compounds, proton pump inhibitors (PPI), and two selective antibiotics out of clarithromycin, amoxicillin, and metronidazole.

Preliminary TIME assessment

The stored slides, pathologically diagnosed using immunohistochemistry were retrieved to analyze tumor-infiltrating T lymphocytes (TIL-Ts), including all 76 surgical specimens and 176 biopsies. CD3 (routine clinicopathological marker) was used to define TIL-Ts. All CD3+TIL-Ts counts were evaluated by two experienced pathologists, and were divided into high-TIL-Ts and low-TIL-Ts groups based on the cut-off value defined as 14%, which was based on Receiver Operating Characteristic (ROC) curves according to a previous study.16

Clinical parameters

All parameters were retrieved from individual medical records. Physical examinations included peripheral lymph nodes, Waldeyer ring, specific hematological indexes, CT, gastroendoscopy, ultrasonography, positron emission tomography-CT, and bone marrow biopsy. For patients undergoing surgery, further consideration was given to intraoperative exploration and pathological authentication. The routine clinical pathological staining results for CD10, BCL6 and MUM-1 were used to classify all tumors into GCB or non-GCB subgroups by pathologists as described previously.9 The Lugano staging system was employed for staging and classification of the tumors. In addition, the performance status was evaluated according to the Eastern Cooperative Oncology Group (ECOG) criteria and International Prognostic Index (IPI). The response of each patient to HPE or chemotherapy was assessed using periodical endoscopic biopsy and other imaging examinations. Tumors that exhibited pCR were classified toward therapeutic effectiveness clinically.

Immunohistochemical staining

Thirty paraffin-embedded tissue slides were subjected to immunohistochemical staining using the following antibodies: PD-L1 (AG27557, Proteintech), PD-1 (AG12470, Proteintech), p65 (AG1199, Proteintech), TRAF6 (AG24296, Proteintech), CagA (A-10, Santa Cruz Biotechnology), tumor necrosis factor-α (TNF-α) (A11534, ABclonal), interferon (IFN)-γ (AG8321, Proteintech), and IL-10 (AG14870, Proteintech). Staining of preprocessed endoscopic biopsy samples and surgical specimens was carried out on embedded slices using the indirect peroxidase method (pH 6.0 in sodium citrate buffer). Two expert pathologists scored the slides independently, bypassing the clinical data. An IScope microscope (IS.1159EPLi, Euromex) and Image Focus V.4.0 (Euromex) software were applied to scan and capture the images. The H-score to quantify the P65, TRAF6, TNF-α, IFN-γ, IL-10, PD-1, and PD-L1 staining intensities was calculated as follows: ∑ (3×percentage of strong staining +2× percentage of moderate staining +percentage of weak staining) using ImageJ with 0–300 fluctuations. The graphic data were displayed as mean±SE of the mean (SEM). For CagA markers, the positive expression exceeded 10% in cells with moderate or strong immunostaining.15

Multiplexed ion beam imaging by time of flight (MIBI-TOF) tissue image

Antibodies and panels

Lymphocyte types as well as activation, cytokines, vascular, and spatial structure of other tissue cells was based on all antibody panels. Staining with the primary antibodies (online supplemental table 1) with metal-conjugation was conducted using the Maxpar labeling kit (Fluidigm) using concentrations based on BioTek (Berten Instruments) at 0.5 g/L and stored long-term at 4°C.

Supplemental material

Tissue antibody labeling

Formalin-fixed and paraffin-embedded tissue slides were stained as listed in online supplemental table 1. After dewaxing, antigen retrieval (R&D systems) and blocking were performed. The samples were incubated with the antibody cocktail overnight at 4°C and then stained with Intercalator-Ir (Fluidigm, 201 192A) for the detection of DNA.

Imaging mass cytometry

According to H&E staining, appropriate 500×500 µm sites were selected for scanning. IMC Images were captured using a multiplexed ion beam imaging by time of flight Imaging System (Fluidigm) at 200 Hz with interspersed capture of isotope polymers (Fluidigm) to monitor the stability. MCDViewer, CellProfiler, and HistoCAT were used for the processing of successful image captures.

Single-cell identification

The original MCD file saved in TIFF format was transformed using imctools (https://github.com/BodenmillerGroup/imctools). The data were segmented into single cells using CellProfiler V.3.1.8.

Data transformation and normalization

The single-cell markers were measured using histoCat v1.75, multiplied by a factor of 107 to yield values larger than 1 and then log-transformed. The single-cell data were deleted at the 95th percentile to exclude the outliers. The heatmaps were presented with z-scored means of clusters. Then, the data were standardized using Harmony in t-Distributed Stochastic Neighbor Embedding (t-SNE) and PhenoGraph.

Clustering

The aggregation of subpopulation groups according to their markers was performed using PhenoGraph (V.2.0) with 15 nearest neighbors clustering. Twenty defined clusters were assigned into a larger aggregation as the hierarchical clustering (Euclidean distance and Ward’s linkage) according to their mean pertinence of related markers. The uncertainty of each subtree was assessed by multiscale bootstrap resampling (R package pvclust, V.2.0).

Barnes-Hut t-SNE

Single-cell data with a high dimension was degraded into two dimensions by the nonlinear dimensionality reduction algorithm (t-SNE). The Barnes-Hut of t-SNE was harmoniously standardized to exclude the default parameters (perplexity=25).

Neighborhood analysis

The pairwise neighbor interactions for cell type significance involved in enrichment or depletion were detected using the CellProfiler Measure Object Neighbors module and neighborhood (https://github.com/BodenmillerGroup/neighbouRhood) functions accompanied by permutation-test-based analysis of the spatial distribution. Four pixels (4 µm) showed the boundary of neighboring cells. Differences with p<0.01 were identified as significant.

Tandem mass tag peptide labeling proteomics

Tandem mass tag peptide (TMT) Mass Tagging Kits and Reagents (Thermo Fisher Scientific) were applied for TMT labeling of the peptides with lysed (Sigma) samples. The tissues were fractionated and shotgun proteomics analyses were performed using an EASY-nLCTM 1200 ultrahigh-pressure liquid chromatography system (Thermo Fisher) interfaced with a Q Exactive HF (X) mass spectrometer (Thermo Fisher) operated in the data-dependent acquisition (DDA) mode. The final spectra of each run were searched solely according to the homo_sapiens_uniprot_2020_7_2.fasta (192 320 sequences) database by using the search software Proteome Discoverer V.2.2 (Thermo), followed by filtering of the retrieved data. The proteins that exhibited significantly different abundance between two groups were classified using the t-test, and the differentially expressed proteins were identified according to fold change (FC) ≥1.2 and p≤0.05 for the upregulated proteins and FC ≤0.83 and p≤0.05 for the downregulated proteins. Gene Ontology (GO), Clusters of Orthologous Groups, Kyoto Encyclopedia of Genes and Genomes (KEGG) and InterPro (IPR) functional analyses were carried out with the non-redundant protein databases (Pfam, PRINTS, ProDom, SMART, ProSite, PANTHER). Volcano plots and cluster heatmaps as well as GO, IPR, and KEGG enrichments were marked based on DPEs. The STRING-db server (http://string.embl.de/) was used to predict the probable protein–protein interactions.

Statistical analyses

Pearson’s χ2 test and Fisher’s exact test were applied to determine the correlation between H. pylori status and clinicopathological characteristics. The primary endpoints were progression-free survival (PFS) and OS. PFS denotes the length of time from initial diagnosis to disease progression, relapse, or death due to random causes. OS refers to the length of time between the initial diagnosis and death or final records during the follow-up of survivors. The Kaplan-Meier (K-M) method was used to draw the survival curves, and the log-rank test was employed to estimate the significance of differences between the K-M curves. The effective variables in univariate analysis were included in the subsequent multivariate Cox hazards regression analysis for estimating the underlying independent risk factors for PFS and OS. The accurate 5-year PFS and OS were described in the statistical analysis of the 10-year follow-up. The statistical analyses were conducted using SPSS V.25.0 (IBM) and GraphPad Prism V.7 (GraphPad), and the results were presented as means±SEM or percentages. Student’s t-test was used for comparisons of test subjects from two samples with different conditions. One-way analysis of variance and Tukey’s multiple comparisons test were used for multiple comparisons. For all tests, a p<0.05 (two pairs) was considered to indicate significance (*p<0.05, **p<0.01, and ***p<0.001; not significant (NS) p>0.05).

Results

Patient characteristics

The baseline clinicopathological parameters of the 252 patients with de novo gDLBCL are listed in table 1. They were segregated into H. pylori positive (128, 50.8%) and negative (124, 49.2%) groups. Major gastric ulcer lesions (68.0% vs 46.8%, p=0.001), massive lodging in the proximal stomach (cardia, stomach fundus) (24.2% vs 6.5%, p<0.001), younger age (<60%, 69.5% vs 54.0%, p=0.011), less B symptoms (89.8% vs 54.8%, p<0.001), higher proportion of normal lactate dehydrogenase (LDH) level (82%.0 vs 62.1%, p<0.001), more GCB origin (39.1% vs 23.4%, p=0.016), and minor IPI score (0–1, 66.4% vs 31.5%, p<0.001) were correlated with the better prognosis of the H. pylori-positive group. Specially, patients with H. pylori infection tended to remain in the IE stage while patients in the H. pylori-negative group had a higher proportion of tumors at IIE stage (IE, 72.7% vs 54.0%, p<0.001). There were no significant differences of ECOG, gender, and chemotherapy regimens between the two groups (p>0.05). Intriguingly, more favorable pCR (39.1% vs 19.4%, p=0.001) and higher TIL-Ts (≥14%, 72.7% vs 46.0%, p=0.012) were noted in the H. pylori-positive group.

Over a median follow-up of 66.6 months (range 0.7–119.2), the 5-year OS was 78.0% for all patients. During the entire study, the median PFS and OS were 49.7 (95% CI 31.5 to 67.9) and 56.2 months (95% CI 26.6 to 85.8), respectively. Significantly favorable 5-year PFS and OS were present in the H. pylori-infected patients with pCR and high-TIL-Ts (termed as ‘cohort 1) than in the others (named as ‘cohort 2’) (5-year PFS, 88.1% vs 70.5%, p<0.001; 5 year OS, 89.2% vs 76.6%, p<0.001), as shown in tables 2 and 3.

Table 2

Univariate and multivariate analysis for PFS in 252 patients

Table 3

Univariate and multivariate analysis for overall survival (OS) in 252 patients

H. pylori-positive patients with pCR and benign TIME had a favorable prognosis of de novo gDLBCL

On follow-up, the median PFS was 120.2 months (95% CI 93.2 to 147.2) for cohort 1, while it was 37.8 months (95% CI 19.1 to 56.5) for cohort 2 (p<0.001). The median OS for cohort 1 and cohort 2 was 138.9 months (95% CI 111.8 to 166.0) and 28.1 months (95% CI 4.9 to 51.3) (p<0.001), respectively. Univariate analysis demonstrated that the absence of B symptoms (p<0.001 for PFS, p=0.001 for OS), earlier Lugano stage (p=0.017 for PFS, p=0.046 for OS), good ECOG state (p=0.037 for PFS, p=0.023 for OS), low IPI scores (p<0.001 for PFS and OS), normal LDH level (p<0.001 for PFS and OS), and H. pylori infection with pCR and high-TIL-Ts (p<0.001 for PFS and OS) were remarkably linked to favorable PFS and OS of de novo gDLBCL. Younger age (p=0.038) contributed to a lower risk for OS. Multivariate analysis revealed that aberrant LDH levels (HR=31.721, p<0.001 for PFS; HR=45.390, p<0.001 for OS), weak ECOG performance (HR=2.698, p=0.047 for PFS; HR=3.137, p=0.043 for OS), and emergent B symptoms (HR=0.406, p=0.021 for PFS; HR=0.392, p=0.035 for OS) were independent predictors of shorter PFS and OS. Importantly, H. pylori-positivity, pCR and high-TIL-Ts state were marked as independent positive prognostic factors for PFS (HR=6.255; p=0.002) and OS (HR=26.305; p<0.001) (tables 2 and 3).

Benign TIME and high response to HPE improved the survival of patients with H. pylori-positive de novo gDLBCL

K-M survival curves indicated that H. pylori-positivity with pCR and high-TIL-Ts (p<0.001), young age (p=0.022), absence of B symptoms (p<0.001), early Lugano stage (p<0.001), good ECOG performance status (p<0.001), regular LDH level (p<0.001), low IPI risk (p=0.024), and proximal location in the stomach (p=0.032) were significantly correlated with better OS of patients with de novo gDLBCL (figure 1A–H).

Figure 1

Effects of Helicobacter pylori with pCR and TIL-Ts status, presence of B symptoms, age, ECOG performance, IPI risk group, LDH, Lugano stage and tumor sites on the overall survival of patients with de novo gastric diffuse large B-cell lymphoma (gDLBCL) according to Kaplan-Meier analysis. (A) H. pylori-positive with pCR and high TIL-Ts (p˂0.001), (B) age <60 (p=0.022), (C) absence of B symptoms (p˂0.001), (D) Lugano stage IE (p˂0.001), (E) ECOG <2 (p<0.001), (F) normal LDH (p=0.001), (G) low and intermediate risk (p=0.024), (H) proximal or ≥2 components (p=0.032) indicated better os in de novo gDLBCL. Distal, antrum, angle or lower body; ECOG, Eastern Cooperative Oncology Group; IPI, International Prognostic Index; LDH, lactate dehydrogenase; pCR, complete pathological remission; Proximal: middle body, upper body, fundus or cardia; TIL-Ts, tumor-infiltrating T-lymphocytes.

Decrease of immunosuppressive components in de novo gDLBCL of H. pylori positive patients with pCR and high TIL-Ts

Among 30 subjects, 26 H. pylori-positive tissues were confirmed by CagA staining (figure 2A). A total of 19 cases were further defined as the cluster 1 subset according to medical records. They were marked by increased inflammatory cytokines such as TNF-α and IFN-γ and decreased immunosuppressive factors such as PD-1, PD-L1, and IL-10 in contrast to the cluster 2 subset (figure 2A,C), and the difference was statistically significant (figure 2B,D). Based on high-quality control of expressed proteins (online supplemental figure 1A–E), a total of 2193 multiple functional proteins covering diverse biological processes were identified by searching against common databases (online supplemental figure 2A–F). Principal component analysis revealed significant differences of protein abundance between groups, with favorable coefficients of variance (online supplemental figure 1F, G), which were grouped in six expression clusters (figure 3A). Compared with the cluster 2 subset, the cluster 1 subset had reduced immunosuppressive markers such as PD-1, PD-L1, CD206, CD163, IL-10, and TGF-β as well as related chemokines such as CCL17, CCL22, CCL5, and CXL5. Conversely, proinflammatory cytokines such as TNF-α, IFN-γ, and IL-12 were comparatively augmented (figure 3B,C). Multidimensional analysis including function, pathways, and structural domains supported the different functional components (figure 4). Specifically, proteins related to the immune response were upregulated (p=0.0029) (figure 4A), while the NF-kB (P ˂ 0.001) and PI3K-AKT-mTOR (p=0.012) growth-factor signaling pathways (figure 4B) were downregulated in the cohort 1 subset compared with the cohort 2 subset. Furthermore, intriguing proteins within TIME exhibited characteristic localization and interactions with intergroups differences (online supplemental figure 3). These data strongly imply that H. pylori-elicited inflammatory responses affect the expression of immunosuppressive components and reshape the TIME in de novo gDLBCL.

Supplemental material

Supplemental material

Figure 2

Higher expression of proinflammatory cytokines combined with lower expression of immunosuppressive factors, NF-κB-related proteins and PD-L1 in de novo gDLBCL of Helicobacter pylori-positve patients with pCR and high TIL-Ts. (A, C) Representative immunohistochemistry (IHC) detection of cagA, p65, TRAF6, TNF-α, IFN-γ, PD-L1, PD-1 and IL-10 of cohort 1 subset and cohort 2 subset. Magnification ×20; scale bar: 100 µm. (B, D) H-score of each protein staining degree. Data in graphs are presented as means±SEM Statistical analyses by Student’s t-test. *P<0.05, **p<0.01; ***p<0.001; ****p<0.0001. Cohort 1, H. pylori-positive with PCR and high TIL-Ts; cohort 2, others. gDLBCL, gastric diffuse large B-cell lymphoma; IFN-γ, interferon γ; IL-10, interleukin; pCR, complete pathological remission; SEM, SE of the mean; TIL-Ts, tumor-infiltrating T-lymphocytes; TNF-α, tumor necrosis factor-α.

Figure 3

Decreased expression of immunosuppressive proteins in the cohort 1 subset compared with the cohort 2 subset. (A) Immune cells were classified as six clusters in all samples. (B) The volcano plot illustrating the differentially abundant proteins in cohort 1 subset compared with the cohort 2 subset. The x axis specifies the fold-changes (FC) and the y axis specifies the negative logarithm to the base 10 of the t-test p values. Vertical and horizontal lines reflect the filtering criteria (Fc ≥1.2 and p≤0.05 for upregulated proteins, Fc ≤0.83 and p≤0.05 for downregulated proteins). Red and green dots represent probe sets for protein expressed at significantly higher (32) and lower (67) levels in the cohort 1 subset compared with the cohort 2 subset, respectively. (C) Heatmap showing the specific upregulated and downregulated immuno-related proteins between the cohort 1 and cohort 2 subsets. Red and blue strips represent significantly higher or lower protein abundance with Z scores ranging from −1 to 1. The x axis llists the samples and the y axis lists gene identification section (IDS) with corresponding protein names.

Figure 4

GO, KEGG and IPR enrichment analyses of the differentially abundant proteins showing an enhanced immune response and downregulated NF-kB and PI3K-AKT-mTOR growth-factor signaling pathways in the cohort 1 subset compared with the cohort 2 subset. (A) GO enrichment of all, downregulated and upregulated proteins categorized into BP, CC and MF suggested a higher enrichment of immune response-related proteins (p=0.0029). (B) KEGG enrichment of all, downregulated and upregulated signaling pathways suggested the downregulation of NF-kB (p˂0.001) and PI3K-AKT-mTOR growth-factor signaling pathways (p=0.012) in the cohort 1 subset compared with the cohort 2 subset. (C) IPR enrichment of all, downregulated and upregulated IPR categories. BP, biological process; CC, cellular component; GO, Gene Ontology; IPR, InterPro (domain); KEGG, Kyoto Encyclopedia of Genes and Genomes; MF, molecular function. MHC; major histocompatibility complex. EVE; novel pseudouridine synthase and archaeosine transglycosylase (PUA) domain-like domin.

Immune crosstalk in de novo gDLBCL of H. pylori-positive tumors with pCR and high TIL-Ts affects the prognosis

Based on the prognostic value of TIL-Ts, a holistic and detailed scenario of TIME was developed. The cohort 1 subset had large amounts of aggregated CD45+ immune cells with an equal distribution of myeloid and lymphocytic lineages. Adjacent CD4+ and CD8+ T cells with similar proportions were located close to the tumor center. Large amounts of effective CD8+ T cells partially displayed co-expression of Tim-3 and PD-1, and bound to PD-L1 niches. FOXP3+ Tregs, Bcl-6+ follicular T helper cells, and non-regulatory CD4+ CD25+ Th cells were in an activated state. CD3-CD56+ natural killer cells were scattered throughout the tissue sections. CD11c+ CD16low macrophages were in a comparatively immunosuppressed state. Non-tumor CD19+Bcl-6- Ki67high B cells implied low malignancy. The striking contrasts in the cohort 2 subset were characterized by myeloid-dominated clusters with sporadic lymphocytes. Less scattered T cells maintained the CD8+PD-1- phenotype indicative of a static status, restricted by PD-1/PD-L1 modulation due to upregulated tumor PD-L1, and were far away from a few individual CD4+ cells. CD19+ Bcl-6+ Ki67high tumor B cells suggested high malignancy (figure 5 and online supplemental figure 4).

Supplemental material

Figure 5

The different spatial distributions of TIME represented by typical immune cell infiltrates between the cohort 1 subset and cohort 2 subset. Upward side, 500×500 µm region, scale bar: 100 µm; downward side, 100×100 µm, scale bar: 20 µm. Cohort 1, Helicobacter pylori-positve with pCR and high TIL-Ts; cohort 2, others; TIME, tumor immune microenvironment. pCR, complete pathological remission; TIL-Ts, tumor-infiltrating T-lymphocytes.

We quantified the above characteristics and found differences between the two groups (figure 6A,B). After the comparison of density and proximity in each related cluster, the cohort 1 subset had higher proportions of functional immune cells and a closer intracellular spatial distance than the cohort 2 subset, such as S9 for CTL and S10 for Th, which is associated with better antitumor immunity (figure 6C–F). Taken together, these data demonstrate that H. pylori-triggered inflammatory reactions perturbed the TIME in favor of an improved prognosis in de novo gDLBCL.

Figure 6

(A) Tumor cells and infiltrating immune cells in the TIME from all regions were clustered according to protein expression markers (total 21). The expression values of each protein are scaled from −2 to 2. (B) The different expression values of 21 resulting clusters (y axis) between the cohort 1 subset and cohort 2 subset (x axis) is shown, scaled from −1 to 1. (C) tSNE embedding for tumor cells and infiltrating immune cells in the time. Cells are colored by either the sample (left) or cluster (right). Heatmaps showing the spatial proximity z scores between pairs of clusters in regions of all neighbors (D), in selected regions of R1 and R2 of the cohort 1 subset (E) as well as those of the cohort 2 subset (F). TIME, tumor immune microenvironment.

Discussion

Despite of pronounced pCR after HPE,17 18 there is a wide spectrum of prognosis among H. pylori-positive de novo gDLBCL patients.19 A study initially proposed that H. pylori infection affects the masking of de novo gDLBCL involved in the immune cross-reactions between H. pylori and the B tumor cells.20 We found a statistically significant higher proportion of high TIL-Ts counts in H. pylori-positive patients than in H. pylori-negative ones, preliminarily indicating a benign TIME. The H. pylori-positive patients with pCR and high-TIL-Ts had significantly prolonged survival, emphasizing the impact of the TIME associated with H. pylori infection. Previous evidence indicates that PD-1hi T cells promote malignant progression in MALT lymphoma.21 We further dissected the precise changes of immune function induced by the H. pylori-elicited inflammatory response and TIME by retaining tissues features. Meaningfully, the reduced immunosuppressive molecules and augmented inflammatory components in the cluster 1 subset possibly hint at partial retrieval of antitumor immunity, indicating a relatively benign TIME in the cluster 1.

We proposed that the H. pylori-elicited inflammatory response endowed autoimmune Th cells with a long-lasting activation status, forming a benign TIME.22 When the malignant B-cell-imposed suppressive TIME comes in contact with an active H. pylori-induced antitumor inflammatory environment, PD-1/PD-L1-restricted suppressive T cells or other negative regulators would be impeded, contributing to an ameliorative ‘cancer-immune set point’ and a favorable prognosis.23 Conversely, in de novo gDLBCL with an initial refractory or inferior prognosis,24 the pace may be such that the formation of tumor-derived TIME lagged behind H. pylori-triggered autoimmunity. The tumor effectively performed immune escape, outstripping the diminishing inflammatory activities. In addition, a clinical situation has emerged in which the tumor recurred in spite of a seemingly durable remission at initial treatment.20 25 One of the possible interpretations is an additional break of cancer-immune homeostasis, ascribed to accelerated time-dependent tumor-imposed immunosuppression, outcompeting the diminishing H. pylori-elicited inflammatory immunity over time. Another novel viewpoint focuses on CagA-impacted activation and polarization of T cells. Indeed, CagA elicits Th1-mediated inflammatory immunity and minor Tregs in the earlier stages of gastric carcinogenesis, while the advanced stages are dominated by Th2-mediated inhibitory immunity.26 Accordingly, our data showed that the CagA-marked cluster 1 subset displayed expansive CD4+CD25+ Th1 cell populations. This observation deepened our understanding of the time-dependent changes of the TIME, which directly governs cancer progression.

We intuitively illustrated the prognostic determinants based on immune-crosstalk insights.14 27 Antitumor immunity was subdivided into the immune-desert phenotype with deficiency of immune activation, the immune–exclusion phenotype with specific immune barriers, and the inflamed phenotype with massive immunocyte infiltration,28 reflecting the divergent clinical outcomes.29 30 Based on an analysis of the spatial distributions of the TIME at initial diagnosis without clinical intervention, we reasoned that cohort 1 represented inflamed tumors with more infiltration by functional immune cells. We accordingly defined cohort 2 as representing immune-desert tumors due to only a few sporadic T cells. More CD3+ T cells got closer to the tumor center and were adjacent to other immune cells in the cohort 1 subset. This was associated with better survival, similar to a precious report for DLBCL.31 This spatial distinction of the TIME was in agreement with the idea of a disturbance due to H. pylori infection and its pivotal role in the ameliorative prognosis.

Given the above evidence that the underlying determinants of survival in de novo gDLBCL are related to the crosstalk within paradoxical microenvironments derived from tumor-mediated and H. pylori-induced states, we felt that this provides novel key points for treatment strategies based on an evaluation of H. pylori and TIME status. H. pylori-elicited spatiotemporal diversification of the TIME poses a challenge in ascertaining the optimal therapeutic opportunity. Faced with relapse in a subset of H. pylori-positive patients following conventional treatment,32 it pointed to the caveat that the effect of TIME is non-negligible and could be used to stratify patients for personalized treatment. At present, the Food and Drug Administration (FDA) has approved two chimeric antigen receptor-T medicines targeting relapsed and refractory B cell lymphoma.33 It is further highlighted that not only H. pylori status but also the diverse TIME should be taken into consideration for better clinical outcomes in de novo gDLBCL. Given the limited patient samples and other unknown adaptive intrinsic differences derived from H. pylori-evoked biological behaviors,34 the influence of H. pylori on the TIME and its correlation with the prognosis should be further verified in a larger cohort. By all accounts, this study is the first to illustrate the mechanism of H. pylori-related prognosis in de novo gDLBCL, providing valuable insights for novel clinical treatment directions.

Conclusions

In summary, H. pylori-elicited inflammatory responses may improve the prognosis via crosstalk with the TIME in de novo gDLBCL. H. pylori infection with pCR and high-TIL-Ts can hence be termed as an independent prognostic factor for de novo gDLBCL and help identify patients that will benefit from precision treatment.

Supplemental material

Supplemental material

Data availability statement

Data are available on reasonable request.

Ethics statements

Patient consent for publication

Ethics approval

All study procedures were approved by the institutional review boards of the Harbin Medical University Cancer Hospital and First-affiliated Hospital of Harbin Medical University. This study was performed in accordance with the Declaration of Helsinki.

Acknowledgments

Harbin Xingyun Medical Inspection Institute and Beijing Jinkorui Medical Laboratory provided technical assistance. The clinical data collection was supported by the Harbin Medical University Cancer Hospital and the First-affiliated Hospital of Harbin Medical University.

References

Supplementary materials

Footnotes

  • YD and WS contributed equally.

  • Collaborators Not avaliable.

  • Contributors YD and WS contribute equally in this study.

  • Funding This study was financed by a major project of the National Natural Science Foundation of China (No. 81730074) and N10 project in the Cancer Hospital of Harbin Medical University.

  • Competing interests No, there are no competing interests.

  • 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.

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.