Prognostic role of tumor-infiltrating lymphocytes in gastric cancer: a meta-analysis and experimental validation.

INTRODUCTION: We performed a meta-analysis and an experimental validation to investigate the association between tumor infiltrating lymphocytes (TILs) and the outcome of gastric cancer (GC) patients to provide prognostic indicators for clinical practice.
MATERIAL AND METHODS: The relative literature of TILs in tumor tissue from patients with gastric cancer was searched from PubMed, Embase, NIH databases, from April 2000 to 31 December 2016. Studies on the prognostic value of TILs as CD3+, CD4+, CD8+, GrB+, and FOXP3+ lymphocytes for GC were retrieved, and also the related references were traced as supplements. Independent screening documents, extracting information and evaluating quality were implemented independently by 2 evaluators according to the inclusion and exclusion criteria, which were then analyzed by meta-analysis using STATA version 12.0 software.
RESULTS: The results indicated that high levels of intratumoral CD8+, CD3+ and CD4+ T cell infiltration were associated with better overall survival(OS) in gastric cancer patients, while high density of intratumoral FOXP3+ T cells was not closely associated with a worse outcome. Additionally, in our study, higher density of granzyme B+ (GrB+) T cell infiltration indicated an optimistic prognosis, and infiltration of a larger number of general TILs also suggested a favorable prognosis by log-rank test analysis.
CONCLUSIONS: This meta-analysis clarified that high levels of CD8+, CD3+, and CD4+ T cell infiltration in tumor tissue showed better OS in GC patients, whereas high density of FOXP3+ T cell infiltration may not be recognized as a negative prognostic factor. These results may provide some useful prognostic indicators for clinical application in gastric cancer. Copyright:

Introduction

Early in 1922, McCarty [1] had already brought up the concept of TILs, and considered the infiltration of lymphocytes into tumor tissue as an antitumor activity of the immune system. However, increasing evidence has indicated that the interactions between immune cells and cancers are amazingly complicated [2]. Recently, a bidirectional role of TILs within the tumor microenvironment in tumor progression has been confirmed. Namely, the TILs can not only suppress tumor growth by destroying cancer cells or inhibiting their outgrowth, but also promote tumor progression either by selecting for tumor cells that are more fit to survive in an immunocompetent host or by establishing conditions within the tumor microenvironment which could facilitate tumor outgrowth [3, 4], which may result from the fact that TILs are a heterogeneous population including T cells with different functional phenotypes (occupying the main parts), B cells, and natural killer (NK) cells [5], which could interrupt the immune balance during cancer progression, because of their different physiological and pathological effects in the tumor microenvironment. Therefore, it is important to discuss the prognostic role of the main parts of TIL subsets separately. The correlation between TILs and clinical outcome have been investigated in many solid cancers [6-8]. Plenty of TILs related studies showed conflicting results in gastric cancer field, the prognostic value of TILs and/or TILs subsets was not yet determined [9, 10]. As a result, the meta-analysis was conducted to identify the prognostic value of TILs and TILs subsets in GC patient stratified by infiltration sites. According to the lymphocyte functional phenotypes, we chose the following subtypes of T lymphocytes as our analysis objectives. CD3, a biomarker of T lymphocytes, is expressed in almost all T lymphocytes; CD8+ cytotoxic T lymphocytes (CTLs) are directly capable of killing tumor cells; granzyme B (GrB) is the main effector molecule of activated cytotoxic T lymphocytes (CTLs) and natural killer (NK) cells, and also is their specific activation marker; CD4+ T helper lymphocytes (Th) are a heterogeneous cytokine secreting class of T lymphocytes; T helper type 1 lymphocytes (Th1) have a crucial role in activating CTLs; T helper type 2 lymphocytes stimulate humoral immunity and activate eosinophils; FOXP3+ regulatory cells (Tregs) play a critical role in immunotolerance and deficiency of anti-tumor immunity and some immune diseases [11-13].

Material and methods

Literature search

We searched the literature using the PubMed database, Embase, and NIH databases from April 2000 to 31 December 2016. The retrieval approach was adopted by combining the keywords with free words. Search terms included: immune cells, immune responses, tumor infiltrating lymphocytes, gastric, stomach, cardia or pyloric, and cancer, carcinoma, tumor or malignancy, and prognosis, prognostic or survival outcome. There were no other limitations in the database searching process. Additionally, all the reference lists of identified articles were also reviewed in order to find potential studies. Two authors were responsible for searching the comprehensive database and evaluating availability independently.

Eligibility criteria

Papers that were eligible for inclusion in this meta-analysis should meet the following criteria: (1) patients were diagnosed with gastric cancer; (2) prognostic roles of the general TILs or T lymphocyte subsets (including CD3+, CD8+, GrB+, CD4+, and FoxP3+ lymphocytes) were evaluated; (3) the TIL measurement must be detected in situ of tumor tissue by application of HE or immunohistochemistry; (4) the lymphocyte infiltration site of the study should be within the tumor tissue, such as tumor parenchyma, tumor stroma; (5) adequate data (HR, and its associated 95% CI) were provided for the further analysis; (6) original article; when two or more studies were published by the same authors, the better quality paper was included. Exclusion criteria were as follows: (1) the studies were of insufficient data; (2) as small sample size would induce publication bias, studies with sample size < 50 were excluded; (3) some styles of literature such as case reports, letters, reviews, conference abstracts and animal trials were excluded.

Data extraction and quality assessment

Two authors extracted the data independently from all eligible studies. Uncertainties were resolved through discussion. The following information was extracted from all included publications: author, country of study, publication year, number of patients, lymph node metastasis, tumor staging, TILs detection method, TILs subsets and distribution site, the cut-off definition, follow-up period and outcome measures. The outcome measure of OS (overall survival), DFS (disease-free survival) or RFS (relapse-free survival) and odds ratio (OR) or hazard ratio (HR) as well as their associated 95% confidence intervals (95% CI) were extracted for each incorporated study.

Patients and tissue samples

To validate the prognostic role of the TIL subset in gastric cancer, a cohort of 147 patients were diagnosed with gastric adenocarcinoma , and with complete 10 years follow-up data were recruited from the First Affiliated Hospital of Xi’an Jiao Tong University at the Department of Pathology. The median survival was 66 months, ranging from 3 to 144 months. The median age was 60 years; 116 patients were male, and 31 were female. The basic information of 147 cancer samples included in this study is shown in Table I. This study was approved by the local ethics committee.

Table I

Patient characteristics (n = 147)

Feature	Frequency, N (%)
Age [years]:
≤ 45	20 (13.6)
45–60	52 (35.4)
≥ 60	75 (51.0)
Sex:
Female	31 (21.1)
Male	116 (78.9)
Tumor location:
Cardia	38 (25.9)
Antrum	66 (44.9)
Corpus	43 (29.2)
Lymph node involvement:
Negative	62 (42.2)
Positive	85 (57.8)
TNM stage:
IA–IB	39 (26.5)
IIA–IIB	40 (27.2)
IIIA–IIIC	68 (46.3)
IV	0 (0)
Grade:
G1	3 (2.0)
G2	58 (39.5)
G3	70 (47.6)
G4	16 (10.9)

Immunohistochemistry

Immunohistochemical staining was performed on paraffin-embedded cancer tissue sections that had been serially sectioned at 4-mm thickness. The sections were deparaffinized in xylene and rehydrated in descending concentrations of ethanol. Antigen retrieval was performed by heat-induced epitope retrieval methods in citrate buffer (pH 6.0). Endogenous peroxidase activity was blocked by incubating in 3% hydrogen peroxide in methanol for 15 min. The sections were incubated for 2 h at room temperature (RT) with primary monoclonal antibodies: CD8 (DakoCytomation, Glostrup, Denmark; 1 : 100 dilution); granzyme B (Novocastra, Newcastle, UK; 1 : 100); OX40 (Novocastra; 1 : 30); FOXP3 (Abcam, Cambridge, UK; 1 : 50). Incubation with Dako EnVision (DakoCytomation, Denmark) for 30 min at RT, and followed by color development with DAB and counterstaining with hematoxylin. Surgical resected chronic inflammatory human tonsil tissue was used as the positive control. The negative control for immunostaining was prepared by incubating tissue sections without primary antibody. The positive TILs were counted by the eyepiece micrometer under the light microscope at 400× magnification. The density of TILs within the tumor specimens were evaluated as follows: (1) immune cells were counted in ten of the highest densities of TILs in different fields of each section; (2) the cells were counted in the intratumoral compartment and the necrotic areas were excluded; (3) two observers counted the cells at the same time, in the same field, using a multiple-lens microscope at 400×; (4) the results were expressed as the mean ± standard error of the mean; (5) the densities of the positive TILs were grouped into TILs-high and TILs-low groups based on the median for further analysis.

Statistical analysis

Reading all the incorporated literature carefully; extracting and collating data according to the requirements of meta-analysis; creating data tables; and reviewing data to ensure its accuracy were performed. STATA version 12.0 (Stata Corporation, College Station, TX, USA) was used to perform the further analysis, with significance defined as a p-value less than 0.05. The fixed-effects model or the random-effects model was considered according to heterogeneity, and the I2 statistic was used to evaluate heterogeneity. If each independent study was homogeneous (I2 < 50%), the fixed-effect model was selected for analysis; if significant heterogeneity existed (I2 > 50%), the random effects model would be used to deal with it. The death risk of gastric cancer patients with different levels of TILs was evaluated by hazard ratios (HRs) and 95% CIs. An observed HR > 1 indicated worse prognosis in patients with high density TILs and a HR < 1 suggested better prognosis. The publication bias was evaluated by visually assessing the asymmetry of an inverted funnel plot, and was quantified by Begg’s and Egger’s tests. Survival rates of our verification samples were analyzed with the Kaplan-Meier method, and differences were compared with the log-rank test. Statistical analyses were performed with SPSS 13.0 for Windows (SPSS, Chicago, IL) and two-sided p < 0.05 was considered statistically significant.

Results

Using the described searching strategy, a total of 995 records were identified in a primary literature search. After screening the titles, abstracts, publication types and full text of each publication, 955 studies were excluded, including reviews, small sample size articles, duplicate articles, animal studies and irrelevant studies. Only 40 studies investigated the correlation between TILs and patients’ outcome in gastric cancer. However, 22 articles were excluded, due to lacking essential results concerning data extraction. Finally, 18 articles were include in the meta-analysis involving 3020 cases. The flow diagram for study selection is shown in Figure 1.

Figure 1

Flow diagram of the study selection process

Study characteristics

Table II summarized the characteristics of all included studies [14-31]. The summarizes should represented by summarized incorporated studies were published from 2003 to 2016 and conducted in China (n = 6), Germany (n = 1), Korea (n = 6), Japan (n = 3), America (n = 1), and Italy (n = 1). Among these, only 3 articles that reported the prognostic value of generalized TILs and TIL subsets were evaluated in 15 studies. Sample size of each study ranged from 52 to 487 patients; H & E (hematoxylin-eosin) staining and immunohistochemical (IHC) staining were applied for the detection of generalized TILs and specific TIL subsets respectively in the primary tumor tissue. The cut-off values contained the median level (n = 14), mean level (n = 1), and other semiquantitative methods (n = 3). HRs and 95% CI for overall survival (OS), disease-free survival (DFS), or recurrence-free survival (RFS) were extracted directly from the included studies. The assessment of publication bias for individual studies is presented in Table III, which showed that 18 studies were of high quality.

Table II

Main characteristics of the included studies

Author, year [ref.]	Country	Sample number (male/female)	Lymph node metastasis (positive/negative)	Tumor stage (I/II/III/IV)	Detection method	TILs phenotype	TILs location	Cut-off for overexpression	Follow-up [months] Median (range)	Outcome
Grogg 2003 [14]	America	110 (72/38)	29/81	15/15/43/37	HE	TIL	IT	Median	1–131	OS
Lee 2008 [15]	Korea	220 (156/164)	165/55	67/53/55/45	IHC	CD3, CD8	IT	75th percentile	66.4 (1–108)	OS
Perrone 2008 [16]	Italy	110 (53/57)	79/31	NA	IHC	FOXP3	IT	Median	36–108	OS, RFS
Haas 2009 [17]	Germany	52 (40/12)	28/24	20/19/10/3	IHC	FOXP3, CD3, CD8, GrB	TS, TP	Median	61 (36–132)	DFS
Shen 2010 [18]	China	135 (90/45)	84/51	40/28/60/7	IHC	FOXP3, CD4, CD8	IT	Median	57 (5–120)	OS
Kashimura 2011 [19]	Japan	123 (89/34)	42/81	80/10/13/20	IHC	FOXP3	IT	Median	1–108	OS, DFS
Kim HI 2011 [20]	Korea	180 (126/54)	107/73	0/81/99/0	IHC	CD3, CD4, CD8 FOXP3, GrB	IT	Median	45	OS
Wang 2011 [21]	China	107 (69/38)	69/38	NA	IHC	FOXP3	IT, PT	Median	62 (2–120)	OS
Chen 2012 [22]	China	152 (117/35)	98/35	10/31/93/18	HE	TIL	TN	19.05/HPF	12–144	OS, DFS
Kim KJ 2013 [23]	Korea	99 (55/44)	57/42	Not clear	IHC	CD8	IT	Median and 60th percentile	59 (1–96)	OS
Kim JW 2014 [24]	Korea	243 (152/91)	124/119	120/41/82/0	IHC	CD3, CD4, CD8	IT	Median	74 (0–123)	OS, DFS
Arigami 2014 [25]	Japan	120 (74/46)	66/54	46/16/36/22	IHC	CD3	IT	Median	36 (1–112)	OS
Geng 2014 [26]	China	100 (61/39)	68/32	NA	ICH	FOXP3	IT	25.0 % field	> 60	OS
Liu 2015 [27]	China	166 (125/41)	118/48	23/41/80/22	IHC	CD3, CD4, CD8, FOXP3	IT, PT	Median	79.59 (36–96)	OS
Kang 2016 [28]	Korea	120 (96/24)	33/87	74/26/19/1	HE	TIL	TS	Median	22.2 (2.1–50.8)	RFS, DFS
Kawazoe 2016 [29]	Japan	487 (327/160)	449/38	0/0/358/129	IHC	CD3, CD4, CD8, FOXP3	IT	Median	NA	OS
Kim EK 2016 [30]	Korea	396 (270/126)	246/156	165/77/143/11	ICH	CD3, FOXP3	IT	Mean value	53.9 (0–84.5)	OS
Dong 2016 [31]	China	100 (72/28)	NA	0/36/34/30	ICH	CD8	IT	Median	36.5 (2–88)	OS, DFS

TIL – tumor-infiltrating lymphocyte, PT – peritumoral, IT – intratumoral, TS – tumor stroma, TN – tumor nest, TP – tumor parenchyma, OS – overall survival, DFS – disease-free survival, RFS – relapse-free survival, NA – not available, IHC – immunohistochemistry, HE – hematoxylin-eosin, DFS – disease-free survival.

Table III

Pooled associations between intratumoral TIL subsets and overall survial of patients with gastric cancer

Subsets	No. of studies	No. of patients	Analysis model	HR (95% CI)	P-value	Heterogeneity I2, p	Publication bias Begg’s p, Egger’s p
CD8+ TIL	7	1630	Fixed	0.699 (0.625–0.782)	< 0.001	4.3%, 0.394	0.368, 0.825
CD4+ TIL	5	1211	Fixed	0.770 (0.673–0.881)	< 0.001	47.5%, 0.107	0.806, 0.578
CD3+ TIL	7	1812	Fixed	0.610 (0.526–0.708)	< 0.001	0.0%, 0.926	0.230, 0.391
FOXP3+ TIL	9	1793	Random	1.382 (0.944–2.025)	0.096	77.4%, < 0.001	0.466, 0.080

The prognostic value of whole TILs in general

Only three articles [14, 22, 27] described the relation between the density of whole TILs and patients’ survival. One of the articles evaluated the association of TILs in tumor nest with OS (HR = 0.55, 95% CI: 0.30–0.99, p = 0.047) and DFS (HR = 0.56, 95% CI: 0.32–0.98, p = 0.043), and one article assessed the relation between TILs in tumor stroma with RFS (HR = 12.193, 95% CI: 1.039–143.055, p = 0.047) as well as DFS (HR = 4.836, 95% CI: 0.917–25.525, p = 0.063). One article described the connection between the intratumoral TILs and OS (HR = 0.21, 95% CI: 0.145–0.286, p = 0.013). As the above results showed, we can approximately draw a conclusion that high density of whole TILs was related to a favorable prognosis in GC.

Pooled analysis of TIL subsets

Fifteen studies were assessed for the prognostic value of TIL subsets in GC patients.

CD8+ lymphocytes

There were 7 articles [15, 18, 20, 24, 27, 29, 31] investigating the prognostic value of intratumoral CD8+ T cells for OS. Our analysis results showed I 2 = 4.3% and p = 0.394, with a good homogeneity. The fixed-effects model was chosen to make a further analysis, and the results indicated that a higher density of CD8+ T lymphocytes demonstrated a significantly longer OS (HR = 0.699, 95% CI: 0.625–0.782, p < 0.001, Figure 2 A). Only 2 works [24, 31] studied the contact of CD8+T cells with DFS in intratumoral tissue, so we did not make the further analysis, due to its unreliable credibility. For the activated CTL (GrB+) cells, only one work [17] separately described the relationship between GrB+ TILs and DFS in tumor stromal tissue (HR = 0.87, 95% CI: 0.33–2.29, p = 0.78) and tumor parenchyma (HR = 0.72, 95% CI: 0.21–2.53, p = 0.61). One work [20] discussed the effect of GrB + TILs in the intratumoral site on OS (HR = 0.828, 95% CI: 0.500–1.369, p = 0.461). It seems that GrB + T cells have no contribution to the prognosis of patients with gastric cancer.

Figure 2

Forest plots of studies evaluating hazard ratios of high level of intratumoral TIL subsets in gastric cancer. Forest plots of the prognostic effect on OS of CD8+ T cells (A), CD3+ T cells (B). Forest plots of studies evaluating hazard ratios of high level of intratumoral TIL subsets in gastric cancer. Forest plots of the prognostic effect on OS of CD4+ T cells (C), FOXP3+ T cells (D)

CD3+ T lymphocytes

Eight papers [15, 17, 20, 24, 25, 27, 29, 30] investigated the prognostic value of CD3+ TILs in gastric cancer patients. Seven studies [15, 20, 24, 25, 27, 29, 30] researched the connection of CD3+ TILs in the intratumoral site with OS, with a pooled HR of 0.610 for TILs-high cases (95% CI: 0.526–0.708, p < 0.001), without any heterogeneity (Figure 2 B). Only one study [17] reported the CD3+ TILs in tumor stroma (TS) and tumor parenchyma (TP), and the results showed that HRs for DFS were 0.86 (95% CI: 0.33–2.23, p = 0.76) and 0.76 (95% CI: 0.29–1.98, p = 0.58) respectively. Thus, it is reliable that a high level of intratumoral CD3+ TILs is associated with better outcome in GC.

CD4+ T lymphocyte

For CD4+ TILs, we found a pooled HR of 0.770 (95% CI: 0.673–0.881, p < 0.001, Figure 2 C) using 5 incorporated studies [18, 20, 24, 27, 29], indicating that a higher CD4+ TILs level in the intratumoral region significantly decreased risk of death for GC patients. Heterogeneity between studies remained not obvious with I2 = 47.5%.

FOXP3+ T lymphocyte

Nine studies [16, 18, 19, 20, 21, 26, 27, 29, 30] investigated the effect of intratumoral FOXP3+ TILs on OS, with a pooled HR of 1.382 with a 95% CI of 0.944–2.025, not reaching statistical significance in this meta-analysis (p = 0.096). A random effects model was used due to the high heterogeneity noted between studies (I2 = 77.4%, Figure 2 D).

Publication bias

Begg’s and Egger’s test were conducted to assess the publication bias of all included studies. The p values for Begg’s test were 0.230–0.806, and for Egger’s test were 0.080–0.825. Additionally, a funnel plot was applied to detect publication bias for the FOXP3+ TILs subgroup, with the largest number of studies. The funnel plot of the intratumoral FOXP3+ T cell infiltration was substantially symmetric (Figure 3). Thus, the results above suggest that publication bias was not evident in this meta-analysis (Table III).

Figure 3

Funnel plots of association of intratumoral FOXP3+ T cell infiltration with overall survival in select studies

Validation of the prognostic value of intratumoral TILs in primary gastric tumor tissue

To validate the prognostic role of the TIL subset in gastric cancer, we did a similar experiment to validate the meta-analysis results. The immune markers we chose were not completely the same as it was mentioned in the meta-analysis. CD3 molecules are expressed on all T cells and functionally T cells are classified as cytotoxic T cells (CD8+) and helper T cells (CD4+). CD8+ T cells are well known to be designated CTLs, as for the CD4+ T helper cells, OX40 are used for identifying the primed CD4+ T cells. The OX40 has induced expression after activation, and is mainly expressed on the activated CD4+ T cells, inducing cytokine production and maintaining a normal immune response. So in our experiment, in addition to considering the quantity of T cells infiltrating in the tumor microenvironment, the functional status was taken into consideration too. Consequently, we selected OX40 (expressing on the primed CD4+ T cells), granzyme B (expressing on the activated CD8+ T cells), and FOXP3 as well as CD8, the same marker with the meta-analysis. CD8+ and OX40+ cells were browned in the cell membrane; the brownish positive particles of GrB were located in the cytoplasm; and the FOXP3+ T cells are colored brown in nuclei (Figure 4). Our experimental results showed that there was no obvious correlation between the intratumoral CD8+, OX40+ as well as FOXP3+ TILs and the overall survival rate respectively, while a higher density of GrB+T cell infiltration indicated an optimistic prognosis (p = 0.007, Figure 5 A). Regardless of the subsets of infiltrating lymphocytes, the density of whole TIL suggested a favorable prognosis (p = 0.009, Figure 5 B).

Figure 4

Images (400×) were the positive T cells of FOXP3 (A), OX40 (B), CD8 (C), and GrB (D), infiltrating intratumoral Tissue

Figure 5

Kaplan-Meier curves. A, B – the survival rate with respect to the density of GrB+ T cells (p = 0.007), and the total number of CD8+, FOXP3+, OX40+ T cells (p = 0.009)

Discussion

The gastric tumor microenvironment is frequently filled with a wide range of immune cells, which have been reported to impact on cancer development, progression, and cancer-related immune reactions, emerging as the hotspot of cancer research [32, 33]. A large number of studies have shown that tumor infiltrating immune cells are closely related to the patients’ prognosis, tumor recurrence and metastasis [34-36]. Due to variance in study design, sample size, and detection site, the direct impact of the TIL status on patient survival remains inconclusive and contradictory. Therefore, the TILs could not be used for grading or staging in clinics. This meta-analysis incorporated 18 studies involving 3020 patients aiming to systematically review the prognostic value of TILs in gastric cancer. Additionally, we also performed a similar study to validate the analysis results. To our limited knowledge, this is the first experimental demonstration with meta-analysis to discuss the prognostic values of various kinds of immune cells in GC. The results from the 3 included studies and our investigation confirmed that high density of TILs was associated with a favorable prognosis. Simultaneously, analysis of the TIL subsets including CD8+, GrB+, CD3+, CD4+ and FOXP3+ T cells was performed too. CD8+ lymphocytes are the main effective cells in the immune response, and a better OS in patients with high level of intratumoral CD8+ T cell infiltration was found in the results of this meta-analysis. However, CD8+ T cells did not reveal a positive prognostic value in our validation study, which may, on the one hand, be due to the limited samples, and on the other hand, be due to the complexity and multiplicity of CD8+ T cell function. CD8+ Treg cells have drawn more attention recently, such as the CD8+ CD28– T cells infiltrating human cancer which could inhibit both T cell proliferation and cytotoxicity [37]. Tumor-infiltrating CD8+ T cells with high CD39 expression exhibited features of exhaustion, and the exhausted CD39+CD8+T cells were present within tumors and invaded or metastatic lymph nodes, but were barely detectable within noninvaded lymph nodes, and inhibited IFN-γ production by responder CD8+ T cells [38]. Additionally, there are many distinct states of CD8+ T cell dysfunction in tumors, such as T cell tolerance to self-antigens and tumor-induced T cell dysfunction [39]. Genome-wide transcriptional analysis has revealed self-tolerant CD8+T cells to harbor a tolerance-specific gene program, with hundreds of genes differentially expressed compared to their naive or memory counterparts [40]. Comprehensive analysis of CD8+ T cell status are needed imperatively. Research for activated CD8+ T cells with the presence of GrB is limited, and the results from other studies did not reveal a positive prognostic value. Due to GrB+ T cells playing a main role of the cytotoxic effect, we added the GrB+ lymphocytes to our experiment, and the result showed a better outcome with higher density of intratumoral GrB+ T cell infiltration (p = 0.007, Figure 5 A). As for the CD4 cells, a high density of intatumoral CD4+ T cell infiltration would decrease the death risk in our meta-analysis results. In our validating experiment, we took OX40 phenotype to denote the CD4 T cells, but it did not indicate better prognostic value as CD4+ T cells. It is known that CD4+ T cells are essential organizers of cell-mediated immunity, participating in each stage of the immune response. Naive CD4+ T cells can be induced to differentiate to specific lineages according to the local cytokine milieu, towards Th1, Th2, Th17 and Treg [41]. And Th17 cells play a potent pro-inflammatory role in certain infections, tumor and autoimmune disorders [42, 43]. The complicated role of Th17 cells in tumors tissue was well reviewed by Qi et al. [44], and Th17 cells provided functions of both antitumoral and protumoral activity in different cancers with different stage. The author pointed out that the tumor microenvironment could affect Th17 cells’ developmental plasticity, which can convert into other subsets of T cells, including Th1 cells and Tregs. Thus, future studies need to focus on defining the molecular mechanisms and identifying the tumor microenvironment factors responsible for the function of Th17 cells in tumor immunity. For the much more heterogeneous CD4+ T cell population, more research is needed.

A high density of FOXP3+ T cell infiltration may increase the death risk in GC [16, 18]. However, one study [17] found that FOXP3+ T cells indicated better progression-free survival, while another study [20] found they had no prognostic value, and our results of the analysis revealed that FOXP3+ lymphocytes may not indicate the patients’ prognosis. These results suggested a complex role of FOXP3+ lymphocytes for GC, which may result from the intricate function of FOXP3+ T cells with respect to their surface phenotypes, cell subsets, and functional heterogeneity [45]. The FOXP3+ T cells are known as a subpopulation of T cells which modulate the immune system, with immunosuppressive activities by producing anti-inflammatory cytokines, such as IL-10 and TGF-β [46], and related to poor prognosis in many cancers [47]. However, in some types of cancer, a high level of FOXP3+ T cells was associated with a positive prognosis, such as colorectal carcinoma and follicular lymphoma. This could be due to Tregs’ ability to suppress general inflammation, which is known to trigger cell proliferation and metastasis [48, 49]. These results suggested that FOXP3+ lymphocytes’ role in the development of cancer is highly dependent on both type and location of the tumor. Consequently, the molecular mechanisms governing Treg recruitment and function in cancer are urgently required to be understood in more detail in individual tumor types. T cells with CD3 positivity, a pan T cell marker, were found associated with better OS by the integrated analysis of the seven studies. Our data showing that the density of whole TILs indicated a favorable prognosis were in line with the above-mentioned results. A study in colon cancer revealed that type, density, and location of immune cells within human colorectal tumors could predict clinical outcome [50], and many other reports favor the better indicator of CD3+ T cell infiltrates in digestive tracts [51-53]. A meta-analysis of T cell infiltrates in ovarian cancer [54] found that higher CD3+ lymphocyte levels were a worse prognostic marker, for which further validation is needed because of insufficient studies included.

The present meta-analysis results revealed that various intratumoral TIL subsets were indeed associated with outcome of patients with GC, and the pooled HRs showed dramatic statistical significance. The publication bias detected by Begg’s and Egger’s test showed less significance and the conclusion was reliable. However, the results from our study should be interpreted with caution because we might have failed to identify some published and unpublished studies with negative results or with limited data that would have affected our pooled estimates. We presumed that the potential sources of bias were as follows: (1) Different cut-off values could result in bias. In our meta-analysis, the median cutoff values were mainly used, and there were still 4 studies that used other methods. (2) Not all of the incorporated studies provided a median follow-up time, so the results were uncertain. (3) The number of included documents was too small; and the research objects were of different race and different region. (4) Although all studies have used immunohistochemical methods, the discrepancy in tissue fixation, antibodies used for T cell detection and other operational inconsistencies may lead to inevitable bias. Although there are some limitations in this meta-analysis, it is possible that the density of lymphocyte infiltration may be a reliable predictor for prognosis in gastric cancer patients. Global multi-center collaborating investigation [55] have been carried out on colorectal carcinoma by Immunoscore by which the infiltrating CD3+ and CD8+ T cells within the cancer center and cancer margin were extensively calculated by digital pathology technique. This initiative may pave the way to utilize tumor infiltrating lymphocytes as a predictive marker for oncology clinical practice in the future.

In conclusion, this meta-analysis suggests that TILs were prognostic markers for OS in GC patients to a large degree. In addition, a high density of intratumoral CD8+, CD3+, and CD4+ lymphocyte indicated good prognosis in GC, while FOXP3+ lymphocytes demonstrated no obvious effect on survival outcomes. However, due to the limitations of this study, TILs cannot fully explain their impact on prognosis. Therefore, further studies should focus on high-quality prospective studies, including a comprehensive clinical pathology, evaluation of information, follow-up strategy, and standardized cut-off values, similar treatment strategies, and multivariate analyses of clinicopathological variables of the patients, which will make the study more standardized. Furthermore, quantitative studies of TILs alone are far from explaining the complex effects of the tumor microenvironment, requiring more rigorous design, larger sample size, more standardized survival analyses, and longer follow-up studies to produce more credible statistics.

Conflict of interest

The authors declare no conflict of interest.