B cells in esophago-gastric adenocarcinoma are highly differentiated, organize in tertiary lymphoid structures and produce tumor-specific antibodies

ABSTRACT

Tumor-infiltrating lymphocytes (TILs) are correlated to prognosis of several kinds of cancer. Most studies focused on T cells, while the role of tumor-associated B cells (TABs) has only recently gained more attention. TABs contain subpopulations with distinct functions, potentially promoting or inhibiting immune responses. This study provides a detailed analysis of TABs in gastro-esophageal adenocarcinoma (EAC). Flow cytometric analyses of single cell suspensions of tumor samples, mucosa, lymph nodes and peripheral blood mononuclear cells (PBMC) of EAC patients and healthy controls revealed a distinct B cell compartment in cancer patients. B cells were increased in tumor samples and subset-analyses of TILs showed increased proportions of differentiated and activated B cells and an enrichment for follicular T helper cells. Confocal microscopy demonstrated that TABs were mainly organized in tertiary lymphoid structures (TLS), which resemble lymphoid follicles in secondary lymphoid organs. A panel of 34 tumor-associated antigens (TAAs) expressed in EAC was identified based on public databases and TCGA data to analyze tumor-specific B cell responses using a LUMINEX^TM bead assay and flow cytometry. Structural analyses of TLS and the detection of tumor-specific antibodies against one or more TAAs in 48.1% of analyzed serum samples underline presence of anti-tumor B cell responses in EAC. Interestingly, B cells were decreased in tumors with expression of Programmed Death Ligand 1 or impaired HLA-I expression. These data demonstrate that anti-tumor B cell responses are an additional and underestimated aspect of EAC. Our results are of immediate translational relevance to emerging immunotherapies.

KEYWORDS:

• Plasma cells
• tumor associated antigen
• antibody
• esophageal cancer

Introduction

Adenocarcinomas of the stomach and the esophago-gastric junction are among the most frequent causes of cancer-related deaths worldwide.¹ Modern multimodal treatments combining surgery with neoadjuvant radiochemotherapy or perioperative chemotherapy significantly improved survival of patients with locally advanced disease.^2,3 Nevertheless, the overall prognosis remains poor and more than 50% of patients will experience disease recurrence after neoadjuvant radiochemotherapy and curatively intended resection.² As therapeutic options in metastatic or recurrent disease are highly limited, there is a high medical need in this complex disease.⁴

Cancer immunotherapy has extended therapeutic options across different tumor entities. Especially immune checkpoint inhibition represents a major breakthrough for cancer therapy,^5,6 including treatment of gastric cancer: The programmed death 1 (PD-1) inhibitor pembrolizumab has recently been approved for the treatment of programmed death ligand 1 (PD-L1) positive gastric and esophageal adenocarcinoma. In contrast to molecular targeted therapies, response to immune checkpoint inhibition can only partially be predicted by expression of the targeted molecule on cancer cells and susceptibility to immunotherapy seems to depend on a large panel of immune related factors (e.g. endogenous immune response, neoantigen burden, composition of the lymphocytic immune infiltrate or site-specific features of the tumor microenvironment).⁷

The lymphocytic composition of the tumor microenvironment is a possible correlate of tumor immunogenicity. A standardized quantification of tumor infiltrating lymphocytes has been proposed based on the observation that tumor infiltrating effector memory T cells are associated with superior survival in colorectal cancer.^8,9 In line with these results, Erdag et al. described a positive prognostic impact of a high infiltration by CD8⁺ T cells in melanoma. Interestingly, a high B cell infiltration was also associated with superior prognosis in this cohort of melanoma patients.¹⁰ The impact of tumor associated B cells on the prognosis of cancer has only been addressed by few studies. In pancreatic adenocarcinoma or hepatocellular carcinoma, a high B cell infiltrate is associated with superior prognosis and a recent publication identified T and B cells as major lymphocytic subsets of prognostic relevance also for gastro-esophageal adenocarcinoma.^11–13 In contrast, patients with B cells or plasma cells in the tumor microenvironment of melanoma or lung adenocarcinoma showed an inferior prognosis.^14,15 These opposing results can be explained by different B cell functions. Similar to T cells, specific phenotypes define functionally distinct B cell subsets, which can promote or inhibit immune responses. The distribution of functionally different B cell subsets within the whole B cell infiltrate in cancer has only been addressed by very few studies and is still widely unknown in gastro-esophageal adenocarcinoma.¹² The relevance of B cells for immune responses to cancer is further highlighted by recent publications focusing on the spatial distribution of tumor associated B cells in cancer. B cells often form clusters in the tumor microenvironment, which are similar to lymphoid follicles in secondary lymphoid organs and have therefore been termed tertiary or ectopic lymphoid structures. Tertiary lymphoid structures (TLS) seem to contribute to anti-tumor immune responses in several kinds of cancer and a formation of these peritumoral B cell clusters in the microenvironment is often associated with superior prognosis.^13,16 Functionally, these anti-tumor effects could be mediated by antigen presentation or antibody production. First evidence described a prognostic relevance of B cells in the tumor microenvironment (based on retrospective immunohistochemical analyses) and we decided to further investigate the composition of tumor associated B cells in gastro-esophageal adenocarcinoma.¹¹ This study provides a comprehensive prospective analysis of tumor associated B cell subsets in gastro-esophageal adenocarcinoma including tumor-specific B cell response. Furthermore, we examined the impact of factors enhancing or inhibiting anti-tumor immune responses on tumor associated B cells.

Results

Tumor infiltrating T and B cells are increased in primary tumor samples and mainly show an activated and differentiated phenotype

Tumor infiltrating lymphocytes (TILs) could be isolated from tumor tissue of 44 patients. Primary tumor samples contained significantly more CD45⁺ lymphocytes (8.5% ± 2.5) than samples obtained from patients following neoadjuvant radiochemotherapy or chemotherapy (3.0% ± 0.96, p < 0.05) or mucosa samples (0.29% ± 0.1, p < 0.01; Figure 1A). Overall lymphocyte counts were low in most pretreated tumor samples and normal gastric mucosa (only 5/19 single cell suspensions of pretreated tumor samples and 4/23 analyzed single cell suspensions of mucosa samples contained relevant B cell infiltrates). As these proportions are too low to allow quantification of B cell subsets, we only present data showing B cell infiltrates in PBMCs, TDLNs and TILs from patients that did not receive any preoperative treatment. As additional analyses, the 19 patients treated with neoadjuvant Radiochemotherapy according to the CROSS protocol were stratified in minor (> 10% vital tumor cells, n = 9) and major responders (< 10% vital tumor cells, n = 10). Comparison of B cell subsets in PBMC of these patients did not show significant differences (data not shown).

Figure 1. Comprehensive flow cytometric analyses of lymphocytic subsets in gastro esophageal adenocarcinoma. Single cell suspensions of primary tumor samples (n = 28), tumors of patients following neoadjuvant treatment (n = 19) or normal mucosa (n = 23) were analyzed for the percentage of CD45⁺ lymphocytes by flow cytometry (A). Lymphocytic subsets in PBMC of healthy controls (PBMC HC, n = 20), PBMC of untreated tumor patients (PBMC AC, n = 46), tumor samples (n = 28) and tumor-draining lymph nodes (TDLN, n = 23) were analyzed for the percentages of total B cells (B), CD86⁺ activated B cells (C), IgD⁻CD27⁺ memory B cells (D), CD3⁺PD-1⁺ activated T cells (E) and CCR7⁻CD45RA⁻ effector memory T cells (F). Additional plots compare B cell subsets in low (UICC I and II) and high (UICC III and IV) tumor stages (B-F). Graphs show mean values ± SEM, p = ANOVA, *p < 0.05.

B cells (CD19⁺CD20⁺ in % of CD45⁺ cells) were enriched in tumor samples (11.0% ± 3.5) compared to PBMC of tumor patients (3.9% ± 2.4, p < 0.01; PBMC AC) or PBMC of healthy controls (5.1% ± 3.5, p < 0.01; PBMC HC; Figure 1B). Single cells suspensions of TDLN contained the highest fraction of B cells (17.9% ± 1.8, p < 0.01; Figure 1B). Analyses of tumor associated B cell subsets revealed a higher percentage of activated (CD86⁺ in % of CD19⁺CD20⁺ lymphocytes, Figure 1C) and effector memory B cell subsets (IgD⁻CD27⁺ in % of CD19⁺CD20⁺ lymphocytes, Figure 1D) in the tumor microenvironment (14.4% ± 2.1 and 34.0% ± 3.6, respectively) compared to PBMC AC (6.9% and 22.0%, respectively, p < 0.01). In line with increased differentiation of B cells, T cells in tumor samples were also mainly of an activated (CD3⁺PD-1⁺ in % of CD3⁺ cells, Figure 1E) and effector-memory phenotype (CCR7⁻CD45RA⁻ in % of CD3⁺ cells, Figure 1F). The total percentage of B or T cells and the relative proportions of the aforementioned T and B cell subsets did not differ between PBMC of tumor patients and PBMC of healthy controls (Figure 1B-F). The distribution of B cell subsets did not show significant differences, if samples were stratified according to tumor grading (G1/G2 vs. G3) and localization of the tumor (esophago-gastric junction vs. stomach).

B cells in the tumor microenvironment are organized in tertiary lymphoid structures

As we observed partial similarities of B cell subsets in tumor samples and TDLN, we hypothesized structural similarities. The spatial distribution of tumor associated B cells (TABs) was analyzed by immunohistochemical staining of CD20 in 72 available paraffin embedded tumor samples. The decreased lymphocytic infiltrate of pretreated samples obtained by flow cytometric analyses was confirmed by immunohistochemical analyses of CD20 in tumor specimens. Whereas 51.6% (15/31) of pretreated tumor samples did not contain B cells in the tumor microenvironment, 92.7% (38/41) of primary tumors showed a relevant B cell infiltrate (p < 0.01). TABs in these primary samples were mainly organized in B cell clusters along the invasive margin (Figure 2A). In line with a formation of tertiary lymphoid structures, follicular T helper cells (TFH; CD3⁺CD4⁺CCR7⁻CD45RA⁻CXCR5⁺ in % of CD3⁺ lymphocytes) were elevated in tumor samples (11.3% ± 2.0) and TDLN (9.7% ± 3.3) compared to PBMC AC (3.4% ± 0.6; Figure 2B, E; p < 0.05). Plasmablasts (CD20⁻CD27⁺CD38⁺⁺ in % of CD19⁺ lymphocytes) were also enriched in tumor samples (12.9% ± 2.1) compared to PBMC AC (5.2% ± 1.1) and TDLN (0.8% ± 0.2; p < 0.05; Figure 2C). Whereas the percentage of plasmablasts was not increased in PBMC AC compared to PBMC HC, our analyses revealed a systemic plasmacytosis (CD20⁻CD27⁺CD38⁺⁺CD138⁺ in % of CD19⁺ lymphocytes) in PBMC of tumor patients (3.5% ± 0.9 vs. 0.3% ± 0.1; p < 0.01). Plasmacells were also detectable in tumor samples (3.4% ± 0.7), but not in TDLN (0.3% ± 0.1; p < 0.01; Figure 2D). Presence of follicular dendritic cells and lymphatic vessels as additional important structural components of TLS could be visualized by multicolor confocal immunofluorescence microscopy (Figure 2E). CD19⁺CD38⁺ plasmablasts were mainly located in the surrounding of TLS along the invasive margin of tumors (Figure 2F).

Figure 2. B cells in gastro-esophageal adenocarcinoma are mainly localized in tertiary lymphoid structures (TLS). B cells in tumor sections of gastro-esophageal adenocarcinoma were stained by immunohistochemistry for CD20 (A). Lymphocytic subsets in PBMC of healthy controls (PBMC HC, n = 20), PBMC of untreated tumor patients (PBMC AC, n = 46), tumor samples (n = 28) and tumor-draining lymph nodes (TDLN, n = 23) were analyzed by flow cytometry for the percentages of CD4⁺CXCR5⁺CCR7⁻CD45RA⁻ follicular T helper cells (B), CD20⁻CD27⁺CD38⁺⁺ plasmablasts (C) and CD20⁻CD27⁺CD38⁺⁺CD138⁺ plasmacells (D). Presence of lymphatic vessels (PNAd) and follicular dendritic cells (FDC-M1) as additional components of TLS is demonstrated by confocal microscopy (E). CD38⁺ B cells were visualized by confocal microscopy of CD19, CD38 and CD3 (F). Graphs show mean values ± SEM, p = ANOVA, *p < 0.05.

Figure 3 summarizes results revealed for B and T cell subsets in PBMC AC, tumor and TDLN normalized for results obtained in PBMC HC. The pattern of differential expression was similar to TDLN for B and T cell subsets involved in antigen presentation (antigen-presenting B cells (BAPC, CD21⁻CD86⁺ in % of CD19⁺CD20⁺) and T follicular helper cells (TFH, CD3⁺CD4⁺CCR7⁻CD45RA⁻CXCR5⁺ in % of CD3⁺), whereas other subpopulations (e.g. CD3⁺PD-1⁺) were only increased in tumor samples or tumor samples and PBMC (e.g. plasmablasts (CD20⁻CD27⁺CD38⁺⁺ in % of CD19⁺, Figure 3).

Figure 3. Heatmap showing an overview of analyzed lymphocytic subsets normalized to results obtained from PBMC of healthy controls (PBMC HC = 1).

Tumor specific B cell response

We searched databases (CTDatabase, Proteinatlas) to identify cancer testis antigens and other tumor associated antigens frequently expressed in esophageal adenocarcinoma. Expression of identified antigens was analyzed using a dataset combining data from The Cancer Genome Atlas (TCGA) with data from the Genotype Tissue Expression program (GTEx) and ArrayExpress (n = 185 EAC samples, 319 normal esophageal samples and 207 normal testicle samples). All analyzed antigens were highly expressed in normal human testicle and several antigens showed higher expression in tumor samples compared to normal tissue (Supplementary Figure 2 and Table 1).

Table 1. Mean fluorescence intensities for 34 tumor associated antigens analyzed in serum samples of gastroesophageal adenocarcinoma patients (positive = MFI> 250, negative = MFI< 250) and healthy controls. Expression levels of analyzed antigens in external cohorts of esophageal adenocarcinoma samples (TCGA n = 185) and normal esophageal mucosa samples (GTEx n = 300, ArrayExpress n = 6, TCGA n = 13).

	Luminex pos.		Luminex Neg.		Healthy controls		Gene expression esophagus tumor/esophagus normal
Antigen	mean MFI	n	mean MFI	n	mean	n	Log 2 fold change	Adjusted p-value
Anti- IgG	13619,6	34	n.a.		13708,8	5	n.a.	n.a.
BSA-myc-DDK	n.a.		20,3	34	17,0	5	n.a.	n.a.
CCDC62	535,9	2	35,1	32	17,2	5	1.667282	4.95E-47
CCNA1	434,3	2	25,8	32	18,8	5	1.94532	5.51E-08
CDKN2A	990,8	3	22,8	31	24,8	5	−0.5388	0.615508
CLP1	330,6	2	23,5	32	14,1	5	−0.52529	2.73E-23
CT45A1	n.a.		39,9	34	19,3	5	4.287943	7.70E-47
CTAG1A	11765,1	5	57,9	29	50,8	5	2.82866	3.74E-30
G250/CAIX	883,4	2	28,3	32	14,6	5	−0.02061	0.293723
GAGE13	876,0	2	43,3	32	36,3	5	0.240787	2.19E-12
GKAP1	384,8	3	30,2	31	17,0	5	−0.77163	1.49E-13
HS Albumin	n.a.		15,1	34	13,0	5	n.a.	n.a.
LAGE-1	12855,8	4	48,8	30	38,5	5	3.477191	2.04E-25
LUZP4	n.a.		28,3	34	27,5	5	2.419656	8.37E-63
MAGEA1	593,8	4	24,7	30	21,5	5	3.374806	0.37873
MAGEA3	2019,7	4	13,7	30	11,6	5	5.92568	1.70E-14
MAGEA4	4169,2	4	37,9	30	22,0	5	4.409419	4.34E-07
MAGEB1	1229,0	2	24,4	32	30,1	5	1.080169	9.37E-28
MAGEB2	1164,8	3	20,5	31	22,6	5	2.662383	1.25E-16
MAGEC2	3493,3	2	39,8	32	23,8	5	1.104916	8.20E-13
MAGED2	290,0	1	25,7	33	13,8	5	−1.31776	8.10E-54
MAGEF1	823,1	2	25,3	32	26,0	5	−0.43488	6.82E-12
MAGEH1	690,2	3	16,4	31	18,2	5	−1.07827	1.13E-22
NXF2	476,4	6	75,3	28	44,5	5	2.765306	1.53E-13
NY-ESO-1	11702,1	4	33,6	30	23,4	5	3.551951	7.61E-33
OIP5	432,8	1	19,5	33	15,2	5	1.604937	1.05E-26
PRAME	307,3	2	34,1	32	29,8	5	3.228187	2.52E-23
RAGE1	556,0	1	48,9	33	20,8	5	0.321387	5.14E-05
SCP1	858,8	3	29,9	31	43,2	5	−0.0383	0.384819
SPAG1	805,5	1	27,4	33	21,1	5	1.106646	1.86E-27
SSX1	782,3	2	17,9	32	16,1	5	2.604618	0.019442
SSX2	1499,9	2	34,3	32	28,0	5	−0.95589	0.001309
SSX4	1464,5	2	30,3	32	25,8	5	3.691177	2.11E-15
TP53	1694,1	5	34,5	29	24,4	5	−0.24766	0.401955
TSGA10	898,7	6	75,3	28	86,3	5	1.029976	2.79E-30
XAGE2	1471,0	1	30,3	33	21,2	5	0.362781	8.11E-35
XAGE3	528,0	2	36,8	32	65,2	5	0.302256	2.06E-26

Based on these results, we designed a custom LUMINEX^TM bead array and could detect antibody responses against 31/34 tested tumor associated antigens (TAAs) in serum samples of tumor patients. Antibodies against at least one of the analyzed TAAs could be detected in 16 out of 33 analyzed serum samples of untreated patients (Figure 4A and Table 1). We did not detect antibody responses against CT45A and LUZP4. In contrast, antibodies targeting CTAG1A and MAGEA4 could be detected in 5/33 analyzed serum samples of previously untreated patients (Figure 4B left and 4B right, respectively). MFI of positive results showed a large variety. Using MFI > 250 as cutoff, our analyses did not reveal any positive results in serum samples of healthy controls, whereas 87 out of 1156 tests in samples of tumor patients revealed an MFI > 250 (p < 0.001; Figure 4C). In addition, we analyzed antibody responses in serum samples of patients which received neoadjuvant radiochemotherapy. CTA-specific antibodies were detected in 2/7 samples. Comparison of antibodies against TAA prior and post neoadjuvant treatment (approximately 2–3 months between the two samples) revealed similar MFIs for most analyzed TAAs (Patient 32–34, Figure 4B). Indirect staining of the B cell receptor using biotinylated protein and PE-conjugated Streptavidin tetramers can be used as an alternative method to identify antigen-specific B cell response. Using this approach, we identified NY-ESO specific B cells in tumor draining lymph nodes of 1/3 patients with esophageal adenocarcinoma (Figure 4D).

Figure 4. Analysis of the tumor antigen-specific B cell response in gastro-esophageal adenocarcinoma. LUMINEX^TM analyses of 34 TAAs were performed in serum samples of untreated gastro-esophageal adenocarcinoma patients (n = 34), serum samples following neoadjuvant chemoradiotherapy (n = 7) and healthy controls (n = 5) was analyzed by LUMINEX^TM (A). Individual MFIs for CTAG1A and MAGEA4 (B). Summary of antibody responses detected in serum samples of tumor patients and healthy controls (C) Exemplary flow cytometry analyses of b cells specific for NY-ESO-1 (CTAG1A) in PBMC and TDLN of a gastro-esophageal adenocarcinoma patient using biotinylated NY-ESO-1 and a streptavidin tetramer (D). Heatmap and bar graphs of LUMINEX^TM data show mean MFIs of duplicates, p = Chi² test.

B cell infiltrates are decreased in tumors containing factors conferring low immunogenicity (HLA-loss and PD-L1 expression)

According to the “cancer immunoediting” hypothesis, reciprocal shaping of tumors and their immune microenvironment frequently induces immune escape.¹⁷ To assess the impact of immune escape on tumor associated B cell subsets, we analyzed loss of HLA-I (HLA-ABC) and expression of PD-L1 as two mechanisms of immune evasion described in gastroesophageal adenocarcinoma.^18,19 Altered HLA-expression with either partial or complete loss was detected in 6/27 (22.2%) of primary tumor samples (Figure 5A, C). Immunohistochemistry showed positive staining for PD-L1 on tumor cells or tumor infiltrating lymphocytes in 16/27 (59.3%) of analyzed samples (Figure 5B, C). Expression of PD-L1 and HLA-loss were not correlated (p = 0.68, Figure 5C). B cells (CD19⁺CD20⁺ in % of CD45⁺ lymphocytes) were reduced in the tumor microenvironment of patients with expression of PD-L1 (7.8% ± 1.4) or loss of HLA-ABC (6.5% ± 2.1) compared to patients without expression of PD-L1 (14.6% ± 3.2, p < 0.01; Figure 5E) or unaltered HLA-ABC expression (14.2% ± 2.0; p < 0.01; Figure 5E), respectively. Comparison of B cell subsets revealed similar percentages of activated B cells, plasmablasts and memory B cells in patients with and without HLA-loss or PD-L1 expression in the tumor microenvironment (Figure 5 D, E). Microsatellite instability can increase anti-tumor immune response and is correlated to improved response to immune checkpoint inhibition.^20,21 As only 3/28 primary tumor samples were identified as microsatellite instable, our cohort was too small to allow conclusions regarding the impact of MSI on TABs.

Figure 5. Influence of immune escape mechanisms on tumor associated B cells. HLA-loss (A, C) and expression of PD-L1 (B, C) was assessed by immunohistochemistry in 27 primary tumor samples. PBMC of these patients (PBMC AC) or single cell suspensions of tumor samples were analyzed by flow cytometry for the presence of different B cell subsets according to their expression of PD-L1 (D) or HLA-loss (E). Bar graphs show mean percentages of flow cytometric data, p = ANOVA, p < 0.05.

Immunosuppressive B cell subsets, which inhibit immune responses by production of Interleukin-10 (IL-10), Granzyme B or expression of PD-L1 have been described.^22–25 However, in our study regulatory B cell subsets (CD24⁺⁺CD38⁺⁺ or CD24^highCD27⁺ in % of CD 19⁺CD20⁺) were not increased in tumor samples (1.6% ± 0.3 and 10.1% ± 2.0) compared to TDLN (0.5% ± 0.1 and 5.3% ± 1.1, p = 0.4 and 0.1), PBMC of EAC patients (3.1% ± 0.5 and 18.5% ± 2.1, p < 0.05) and PBMC of healthy controls (3.1% ± 0.4 and 20.0% ± 2.2, p < 0.05 and 0.2; Supplementary Figure 3A, B). In addition, PD-L1⁺ B cells were not enriched in the tumor microenvironment (1.9% ± 0.4) compared to PBMC AC (1.8% ± 0.3, p = 1.0) and PBMC HC (1.0% ± 0.1, respectively, p = 0.9, Supplemetary Figure 3C). In contrast, regulatory T cells (CD3⁺CD4⁺CD25⁺CD127^low in % of CD3⁺) were elevated in tumor samples (8.2% ± 0.8) compared to PBMC AC (4.4% ± 0.3) and PBMC HC (2.0% ± 0.4, p < 0.01; Supplementary Figure 3D).

Discussion

Type, density and localization of tumor infiltrating T cells is correlated to prognosis in several kinds of cancer.^20,26–28 Only a small subset of these studies included quantification of TABs.²⁸ A high density of TABs (analyzed based on detection of CD19 or CD20 by immunohistochemistry) is associated with superior survival in the majority of studies in gastro-esophageal adenocarcinoma.^11,29–31 Svensson et al. analyzed tumor associated B and T cells in a cohort of esophageal adenocarcinoma patients. Patients with a combination of a high B and T cell infiltrate had the best prognosis, suggesting that tumor associated B and T cells synergistically inhibit tumor growth.²⁹

Distinct B cell subsets have been identified, which potentially promote or inhibit anti-tumor immune response.^32–34 The data described above highlight the prognostic relevance of B cells in solid cancer, but only very few studies analyzed the composition of B cell subsets in the tumor microenvironment.^35–37 Our current study demonstrates, that B cells are regularly detectable in the microenvironment of gastro- esophageal adenocarcinoma and predominantly consist of activated and differentiated subsets. Plasmablasts were particularly elevated in PBMC and tumor samples of gastro-esophageal adenocarcinoma patients. A similar systemic plasmacytosis and elevation of plasmablasts in tumor samples has been described in melanoma and other solid tumors, most likely reflecting anti-tumor B cell responses in these patients.^38–41

Garnelo et al. found a similar increase of plasmablasts using multicolor immunofluorescence microscopy to identify B cell subpopulations in hepatocellular carcinoma. They demonstrated a positive prognostic impact of B cells and CD38⁺ plasmacells in this disease. As costimulatory molecules and inflammatory cytokines were elevated in tumor samples with a high B cell infiltrate, they concluded that B and T cells act synergistically in the tumor microenvironment.¹² Furthermore, TABs in hepatocellular carcinoma were localized in tertiary lymphoid structures (TLS), whose presence is correlated to superior prognosis in several other types of solid cancer.^{13,16,34,42,43} In line with an increase of differentiated B cell subsets in our flow cytometric data, TLS were detectable in the vast majority of our primary tumor samples (40/43). In addition to a clustering of T and B cells we could confirm the presence of high endothelial vessels and follicular dendritic cells as additional important structural components of TLS (Figure 2E). CD19⁺CD38⁺ plasmablasts were mainly localized along the invasive margin and outside of TLS (Figure 2F), which has been described similarly for other tumor entities.^12,34,41 One recent publication suggested that immunosuppressive effects of circulating regulatory B cell subsets promote tumor growth in esophageal adenocarcinoma.⁴⁴ Our results do not confirm these data, since the two immunosuppressive B cell populations (CD24^highCD38^high and CD24^highCD27⁺ B cells) were not enriched in the tumor microenvironment.

The majority of clinical studies described a positive prognostic impact of tumor associated B cells and our data describing highly differentiated and organized B cell infiltrates further supports this observation.⁴⁵ However, the exact mode of action of anti-tumor B cell effects is still widely unknown and this aspect needs further investigation. Antibody responses against TAAs could be described in cancers of different origins, mainly focusing on one specific antigen. Frequencies of antibodies against TAAs varied according to analyzed antigens and cancer sites. Some of these publications demonstrated a prognostic relevance of antibody-mediated anti-tumor immune responses.^46,47 Two larger studies analyzed antibody responses against a panel of TAAs in gastro-esophageal adenocarcinoma.^48,49 Antibodies could be detected in 0–39.1% of patients, depending on the analyzed antigen. Chen et al. described coexpression of more than four out of eight cancer testis antigens in 31% of esophageal adenocarcinoma patients.⁴⁸ In the present study we identified a larger panel of 34 TAAs and could detect antibody response in at least one patient for 31 of these antigens. Although our expanded panel increased the frequency of detectable antibodies compared to previous studies, we only detected antibodies in 48.5% of serum samples. Using a similarly large panel of antigens, Germain et al. described tumor-specific antibodies in supernatants of isolated and in vitro stimulated tumor infiltrating B cells of lung adenocarcinoma patients in 41% of analyzed samples. In line with our results, some patients showed response to a broad spectrum of antigens, whereas the majority of patients had detectable antibodies against only one or two different TAAs.³⁴

Finally, our study provides first evidence for impaired anti-tumor B cell responses in patients with expression of PD-L1 or HLA-loss as two important mechanisms of immune escape. The total B cell infiltrate was decreased in both conditions, but the composition of infiltrating B cell subsets did not show differences in tumors with and without evidence of immune escape. Decreased lymphocytic infiltrates have been described for patients with evidence of HLA-loss in the tumor microenvironment of other cancers,^50,51 whereas increased expression of PD-L1 was mostly associated with increased lymphocytic infiltration.^52,53 With few exceptions, these studies did not include B cells into their analyses. Nevertheless, our data provides first evidence for a correlation of the analyzed immune escape mechanisms and B cell infiltrates and the impact of immune escape mechanisms on TAB should be further analyzed. The fraction of MSI-positive tumor samples was too low, as our cohort contained only 6 tumor samples with microsatellite instability of which 2 had received neoadjuvant treatment. In line with previous data showing a link between HLA-loss and MSI, 5 out of 6 MSI⁺ tumors showed impaired HLA-expression.⁵⁴

Taken together our results are highly suggestive for tumor-specific B cell responses in gastro-esophageal adenocarcinoma. The present study is of immediate translational relevance since B cells could be an important additional aspect influencing therapeutic efficacy in ongoing clinical trials evaluating the role of immunotherapy in gastro-esophageal adenocarcinoma.

Materials and methods

Patients

89 patients were included into this prospective study. Peripheral blood mononuclear cells (PBMC) were isolated from peripheral blood using density gradient centrifugation (Pancoll, PAN Biotech). Resection specimens were transferred to our Institute of Pathology immediately after surgery (29 gastrectomies, 53 Ivor-Lewis esophagectomies). Samples with very small or undefined tumors had to be excluded to ensure state of the art pathological assessment according to the 8^th edition of the Union for International Cancer Control (UICC) staging system. Tumor samples, normal mucosa and tumor draining lymph nodes (TDLN) were provided from 35 primary tumor samples and 20 cases following neoadjuvant chemotherapy or chemoradiotherapy. Processing of 28/35 untreated and 15/20 pretreated tumor samples revealed sufficient cell numbers for further analyses. Patient characteristics are summarized in supplementary Table 1. Written informed consent was obtained from all patients and this study was approved by our institutional ethics committee (No. 11–116).

Production of single cell suspensions and ex vivo 10-color flow cytometry

Tissue specimens were processed to single cell suspensions using mechanical (“GentleMacs” Tissue dissociator, Miltenyi Biotech) and enzymatic digestion (320U/ml Collagenase IV, Worthington, and 100U/ml DNAse I, Applichem). Subsequently, isolated PBMC and single cell suspensions were stained for 10-color flow cytometry (see supplementary table 2 for used clones and providers). Data were acquired on a Gallios flow cytometer (Beckman Coulter) and analyzed using the Kaluza software (Version 1.5a, Beckman Coulter). Dead cells were excluded using aqua dead cell staining (Life Technologies).

Immunohistochemistry (IHC) and immunofluorescence (IF) microscopy

Whole slides of tumor specimens were stained for CD20 and HLA-ABC using standard IHC (AEC EnVision+, Dako). Antigen retrieval for IHC and IF was performed on a LabVision PTmodule (Thermo Fisher) using pH6.0 citrate buffer. Stainings for PD-L1 were performed using the Bond-Max System by Leica. Ready to use anti-MLH-1 and anti-MSH-6 antibodies (Roche/Ventana) were used with DAB on a Ventana Benchmark according to the manufacturer’s instructions. Slides were analyzed by an experienced GI-pathologist (Quaas A.). Expression of HLA-ABC was categorized in overall positivity, partial loss and complete loss. PD-L1 was considered positive, if > 1% of tumor cells or tumor infiltrating lymphocytes showed membranous staining. Patients with loss of MLH-1 or MSH-6 were categorized as microsatellite instable (MSI). Confocal microscopy was performed on a Leica TCS SP8 gSTED superresolution microscope. FIJI ImageJ (Version 2.0.0) was used to adjust for brightness and autofluorescence. See supplementary table 2 for detailed information on used antibodies and dilutions.

Gene expression analyses

A set of 34 tumor associated antigens was designed based on available data (Pubmed, Proteinatlas, CTAdatabase). Expression of these antigens in esophageal adenocarcinoma was analyzed in pooled data of 185 esophageal adenocarcinoma patients downloaded from the cancer genome atlas (TCGA),⁵⁰ and 319 normal controls from GTEx⁵⁵ and ArrayExpress.⁵⁶ RNA-Seq data were mapped against the human genome version hg19 with TopHat2-2.0.12. R-3.4.1 and Bioconductor 3.6 were used for the RNA-Seq analysis. Reads were counted using the R package GenomicAlignments (mode = ’Union’, inter.feature = FALSE), only primary read alignments were retained. In order to normalize the read counts we used the equation below.

$normalizedreadcounts=\frac{readcounts}{DESeq2NormalizationFactor\ast \left(genelength/1000\right)}$

DESeq2 normalization factors were obtained using DESeq2. Log 2 fold change values were obtained by comparing the normalized read counts. P-values were calculated using the Wilcoxon Rank Sum Test and corrected using the Benjamini & Hochberg method. Supplementary Fig. 2 was created using ggplot2_1.0.0 on the log 2 values of the normalized read counts. The dendogram of the heatmap was generated using hierarchical clustering with complete linkage.

Analyses of tumor-specific antibodies

A customized TruePlex^TM (Origene) antibody profiling array was used to analyze the predicted 34 tumor associated antigens. Median fluorescence intensities of anti-IgG antibodies were analyzed in 38 serum samples of esophageal adenocarcinoma patients on a LUMINEX 200^TM platform according to the manufacturer’s instructions (Origene). Briefly, HEK293T-expressed recombinant human proteins with C-terminal myc-DDK tag or E. coli-overexpressed and purified human proteins with N-terminal His tag were coupled to LUMINEX^TM beads and incubated with serum samples. Median fluorescence intensity (MFI) for each of the analyzed antigens was measured using fluorescence-labeled anti-human IgG antibodies. In addition to included negative (human and bovine serum albumin) and positive controls (IgG), serum samples of five healthy controls (living-kidney donors) were analyzed.

Antigen-specific B cell response

NY-ESO-1 specific B cells in TDLN were identified based on our previously established protocol using antigen-biotin-tetramers.⁵⁷ Briefly, NY-ESO-1 protein (Origene) was biotinylated (EZ-Link NHS-Biotin reagent, Thermo Fisher) and coupled to PE-conjugated streptavidin (Biolegend) in a ratio of 1:3. Cells were incubated with tetramer for 15 min in addition to B cell antibodies and analyzed by flow cytometry.

Statistical analyses

Statistical analyses were performed using Prism 7 (Graphpad). Means of lymphocytic subsets were compared using ANOVA or Kruskal-Wallis, depending on normal distribution as assessed by D’Agostino-Pearson omnibus normality test. Data was corrected for multiple comparisons (either Dunn’s or Tukey’s Test) and p values < 0.05 were considered significant. Results are presented as mean values ± standard error of the mean (SEM).

Conflicts of Interest

S.R.: Honoraria for advisory boards from BMS and MSD. Honoraria for invited talks from BMS. Financial support for research projects from AstraZeneca; M.B.,: Honoraria for advisory boards, for invited talks from BMS and financial support for research projects from Astellas, Roche and MSD. All other authors declare no conflicts of interest. TZ Honoraria for invited talks from MSD, Merck, Roche, Lilly, Novartis, BMS.;

Additional information

Funding

The “Sander Stiftung” (Nr. 2014.001.01, to H.S., M.B., A.H.), a “Cologne Fortune” and “Gerok” local research grant (to H.S.) and the German Research Foundation (DFG 325827080 to H.S.) supported our research;Sander Foundation [2014.001.01]; Cologne Fortune and Gerok (local research grants).