Prognostic value of immune cells in the tumor microenvironment of early-stage lung cancer: a meta-analysis

Abstract

Background: Early-stage non-small cell lung cancer (NSCLC) patients carry significant risk of recurrence post-surgery. In-depth characterization of the immune tumor microenvironment (TME) can have prognostic value. This study aimed to evaluate the association of individual immune cell types in the TME with clinical outcomes in surgically resected, early-stage NSCLC.

Methods: We performed a systematic literature search of the National Library of Medicine database through November 2019, investigating predefined biomarkers (CD3+ T cells, CD4+ T helper cells, CD8+ cytotoxic T cells, CD20+ B cells, CD56+ & CD57+ Natural Killer (NK) cells, CD68+ Tissue Associated Macrophages (TAMS), FoxP3+ T regulatory cells, and Mast Cells (MC)), and their association with survival following PRISMA guidelines.

Results: Studies that adjusted for important clinical covariates (such as stage and age) showed that higher levels of CD8+ cytotoxic T cells were associated with improved OS (HR = 0.68; 95% CI, 0.50–0.93) and DFS (HR = 0.60; 95% CI, 0.41–0.87), while increased CD20+ B cells (HR = 0.16; 95% CI, 0.04–0.64) and CD 56/57+ NK cells (HR = 0.50; 95% CI, 0.26–0.95) were associated with improved OS; lung cancers with increased FoxP3+ T regulatory cells (HR = 2.22; 95% CI, 1.47–3.34) had worse OS.

Conclusions: Immune cell components of the TME have prognostic value in early-stage, surgically resected NSCLC, and may reveal which patients are more likely to need additional systemic treatment, including immunotherapy. Clinical covariates need to be considered when evaluating the prognostic value of immune cells in the TME.

Introduction

Immune cells within the tumor microenvironment (TME) play an important role in the development, progression and outcomes of non-small cell lung cancer (NSCLC). Innate and adaptive immune cells are able to detect and eliminate malignant transformed cells through the process of immunosurveillance [1]. However, lung cancers that become clinically apparent acquire resistance mechanisms to escape the anti-tumor immune response [1–3]. The relative balance of antagonistic effector (i.e., CD8+ cytotoxic T cells) and regulatory (i.e., FoxP3+ T regulatory cells [Tregs]) immune cell subpopulations may tilt the TME to be either detrimental or supportive of tumorigenesis, and will have a profound impact on the tumor’s eventual destiny [4]. Therefore in-depth characterization of the immune cell composition of the TME is critical to understanding cancer outcomes and may help guide treatment decisions for lung cancer patients.

Resection of early stage NSCLC represents the best opportunity for meaningful long-term survival and cure. However, despite complete removal of all detectable disease, there remains a significant risk of lung cancer recurrence [5]. Predicting which patients are most likely to have recurrence following surgery is of great clinical importance. Patients within the same TNM stage exhibit wide variations in recurrence rates [5]. Efforts to improve post-surgical outcomes using adjuvant chemotherapy have only provided marginal benefit in a subset of patients, while exposing all patients to significant toxicity. New prognostic biomarkers based on immune cell signatures that predict survival outcomes of early stage NSCLC can help identify patients that are most likely to receive benefit from additional systemic treatment, including neo-adjuvant or adjuvant immunotherapies which are now being studied.

A number of retrospective studies have shown individual immune cells within the TME are associated with survival in various malignancies. In colorectal cancers, the type and density of tumor infiltrating lymphocytes (TILs) were found to be more powerful prognostic factors than standard anatomical staging criteria (i.e., TNM) [6]. While CD8+ TILs have been the most studied, mounting data is demonstrating many other immune cell types in the TME also play an important role in cancer outcomes [7]. However, findings from these studies are often not consistent and are limited by small sample size. Three previous meta-analyses evaluating the association between immune cells in lung cancer outcomes are also limited because they included studies on patients with both local and advanced disease and did not account for potentially important clinical confounders [8–10].

Here, we have conducted a meta-analysis of studies evaluating the association of individual immune cell types in the TME with clinical outcomes of surgically resected, early stage NSCLC in order to determine novel prognostic biomarkers for this subset of patients.

Materials and Methods

Search strategy

We conducted a systematic literature search of the National Library of Medicine database to search for all original, retrospective observational studies reporting postoperative survival outcomes of surgically resected, stage I-III NSCLC patients according to immune cell biomarkers (measured by immunohistochemistry [IHC]). There were no date restrictions and the search was finalized in November 2019. Additionally, the cited references of each study (including published reviews and other meta-analyses) were reviewed and evaluated for eligibility. The specific immune cell types studied were selected based on prior literature showing their role in lung cancer prognosis and included: CD3+ T cells, CD4+ T helper cells, CD8+ cytotoxic T cells, CD20+ B cells, CD56+ & CD57+ Natural Killer (NK) cells, CD68+ Tissue Associated Macrophages (TAMS), FoxP3+ T regulatory cells, and Mast Cells (MC). For each of these biomarkers the search terms included: the relevant biomarker, “lung cancer”, and “survival”. A separate Medline search was performed for each biomarker.

Selection criteria

Articles were first screened for relevancy within the scope of the project by independent review of titles and abstracts by two reviewers (ST and WLC). Articles that met qualification underwent further scrutiny through full text review to assure they met all the selection criteria. Disagreements in screening and selection were adjudicated by group consensus involving a third reviewer (ET). Studies were considered eligible for inclusion in this systematic review if they reported on: 1) human subjects, 2) stage I-III NSCLC patients who underwent surgical resection, 3) at least 10 patients, 4) survival data for at least one of the predefined biomarkers or provided sufficient data to estimate survival, and 5) follow-up of at least 4 years. Outcomes of interest were overall survival (OS) and disease-free survival (DFS), either as unadjusted estimates or after adjustment for clinical covariates.

Data extraction

All relevant descriptive information was extracted from each study to create a standardized tabular summary, including author, year of publication, biomarker data, tumor histology, tumor stage, duration of follow-up, number of patients included, gender (when reported), smoking status (when reported), and hazard ratio [HR] and 95% confidence interval [CI]) for OS and/or DFS. When studies reported separate survival for immune cell biomarkers in the tumor compartment and the stromal compartment, tumor HR was preferentially chosen over the stromal HR as tumor infiltrating cells are likely more clinically relevant.

Statistical methods

A meta-analysis was conducted using a linear mixed-effects model to determine the meta-estimate of the average effect, either OS or DFS [11]. The presence of heterogeneity across studies was tested with the Q and I² statistics [12]. The results of the meta-analyses were graphically summarized using forest plots created with the metafor package in R Studio (version 3.2.2; R Foundation for Statistical Computing, Vienna, Austria) [13].

Results

We identified and screened for relevance 2,650 studies; 2,139 of them were irrelevant to the aims of the study. The full text of the remaining 511 articles were reviewed for inclusion, and 45 unique studies were found to be eligible, accounting for 8,471 patients (see Figure 1 for Preferred Reporting Items for Systematic Reviews and Meta-analyses [PRISMA] guidelines). Reasons for ineligibility included having reported survival less than 4 years after surgical resection, inclusion of stage IV NSCLC cases, and not having reported a HR and 95% CI. Of the included studies, 22 reported survival measures for multiple biomarkers of interest, and 16 studies reported both unadjusted and adjusted HRs (Table 1).

Figure 1: PRISMA. * Eligibility Criteria 1) human subjects, 2) stage I-III NSCLC patients who underwent surgical resection, 3) at least 10 patients, 4) survival data for at least one of the predefined biomarkers or provided sufficient data to estimate survival, and 5) follow-up of at least 4 years.

Table 1: Characteristics of included studies

Study, Year	Biomarker	Location	Outcome	HR	Covariates in adjustment model	n of Patients	Male (%)	Stage (%)	Histology (%)	Smoking Status	Cutoff Value
Takanami, 2000 [25]	MCs	Tumor	OS	Adj	gender, T, N, differentiation, microvascular density	180	56	I: 55	ADC: 100	NR	MCD>21
								II: 14
								III: 31
Pelletier, 2001 [26]	CD20	Peritumor	OS	Adj	stage, histology, gender	113	60	I: 58	ADC: 50	Current: 92%	>10%
								II: 18	SCC: 42
								III: 24	LCC: 3
									Other: 5
Kojima, 2002 [27]	CD68 MCs	Intratumor	OS	Unadj and Adj	VEGF, microvessel density	132	67	I: 100	ADC: 71	NR	Mean
									SCC: 29
Villegas, 2002 [28]	CD57	Intratumor	OS	Unadj and Adj	stage, age, endoscopy localization	50	98	I: 82	SCC: 100	Current: 74%	Median
								II: 12		Former: 22%
								III: 6		Never: 4%
Pelosi, 2004 [29]	MCs	Tumor	OS & DFS	Unadj and Adj	age, diameter, pT, Ki-67 labelling index, HER-2, tumor grade, symptoms, blc-2, synaptophysin	201	91	I: 100	ADC: 44	NR	≥5%
									SCC: 56
Kojima, 2005 [30]	CD68 MCs	Tumor		Unadj and Adj (adenocarcinoma only)	tumor VEGF, tumor VEGF-C, tumor VEGFR-3, microvessel density	129	48	I: 100	ADC: 100	Smoker: 40%	Staining score >3
										Non-smoker: 60%
Petersen, 2006 [31]	CD3	Tumor	DFS	Unadj	n/a	64	53	I: 100	ADC: 46	Mean pack-year 51 ± 33	Score ≥2 (median)
									SCC: 34
									Other: 19
Kikuchi, 2007 [32]	CD8	Cancer-Nest	OS	Unadj	n/a	161	68	I: 59*	ADC: 52	Never: 24%	1-cell increment
	CD56							II/III/IV: 41	SCC: 42	Ever: 67%
									Other: 6	Unknown: 9%
Dieu-Nosjean, 2008 [33]	CD3	Tumor	OS & DFS	Unadj	n/a	74	81	I: 84	ADC: 62	Current: 91%	1.5 mean score
	CD20							II: 16	SCC: 38	Never: 9%
Shimizu, 2010 [34]	FoxP3	Intratumor	DFS	Adj	nodal involvement, COX-2 expression	100	60	I: 68	ADC: 69	NR	≥3%
								II: 14	SCC: 31
								III: 18
da Costa Souza, 2012 [35]	CD4		OS	Adj	stage and histological type	65	69	I: 31	ADC: 58	Pack-years median 41(0–120)	CD4: ≥16.1%
	CD8							II: 51	SCC: 31		CD8: ≥1.8%
	CD68							III: 18	LCC: 11		CD68: ≥4.5%
Hanagiri, 2013 [36]	FoxP3	Regional Lymph Nodes	OS	Unadj and Adj	gender, pack-year index, T factor, N factor	158	65	I: 72	ADC: 59	NR	>0.5% PBL
								II: 11	SCC: 22		>1.1% RLNL
								III: 19	Other: 9
Suzuki, 2013 [37]	FoxP3	Stroma	DFS	Unadj and Adj	gender, stage, lymphatic invasion, IL-12R, IL-7R	956 (478)	36	I: 100	ADC: 100	Current: 16%	Score ≥2
										Former: 69%
										Never: 15%
Feng, 2014 [38]	CD68	Islet and Stroma	OS & DFS	Adj	mediastinal down staging, islet/stromal machrophage ratio	28	54	III: 100	ADC: 70	Smoker: 39%	median
									SCC: 21	Never: 61%
Germain, 2014 [39]	CD20	Tumor	OS	Unadj	n/a	74	81	I: 84	ADC: 62	Current: 91%	CD20: 0.0255 mm2/tumor IPF
								II: 16	SCC: 38	Never: 9%
Hanagiri, 2014 [40]	FoxP3	Regional Lymph Nodes	OS	Unadj and Adj	gender, age, histology	131	57	I: 100	ADC: 76	NR	Relative expression > 0.06
									SCC: 16
									Other: 8
Lee, 2014 [41]	CD68		DFS	Adj	procedure, stage, CD68, FoxP3/CD3	151	40	I: 100	ADC: 77	Current: 7%	Score ≥2
									SCC: 23	Former: 72%
										Never: 21%
Djenidi, 2015 [42]	CD3	Total (tumor + stroma)	OS & DFS	Unadj	n/a	101	68	I: 100	ADC: 45	Ever: 89%	Continuous variables
	CD8								SCC: 42
									Other: 13
Hernandez-Prieto, 2015 [43]	CD3	Stroma	DFS	Unadj	n/a	84	86	I: 71	ADC: 48	Current: 43%	Moderate 10%
	CD4							II: 29	SCC: 46	Former: 51%	Strong 20%
	CD8								LCC: 2	Never: 6%
	CD20								Other: 4
	CD57
Kadota, 2015 [44]	CD3	Tumor	OS	Unadj	n/a	331	60	I: 57	SCC: 100	≤90 pack-year: 81%	CD3: High ≥50
	CD4							II: 28		>90 pack-year: 19%	CD4: High ≥20
	CD8							III: 15			CD8: High ≥50
	CD20										CD20: High ≥20
	CD68										CD68: High ≥50
	FoxP3										FoxP3: High ≥20
	Neutrophil-(CD10)										CD10: High ≥10
Kim, 2015 [45]	CD8	Tumor	OS & DFS	Unadj	n/a	331	96	I: 40	SCC: 100	Smoker: 91%	Median
								II: 36		Never: 5%
								III: 24
Donnem, 2015 [46]	CD8	Stroma	OS & DFS	Adj	stage, differentiation, histology, age	797	64	I: 51	ADC: 48	NR	High density >50%
								II: 35	SCC: 45
								III: 15	LCC: 7
Li, 2015 [47]	CD68		OS & DFS	Adj	TNM, N stage, subcarinal lymph node, n number of nodal stations involved, number of involved nodes	159	69	I: 23	ADC: 26	Smoker: 31%
								II: 36	SCC: 52	Never: 69%
								III: 24
O’Callaghan, 2015 [48]	CD3	TI/S	OS	Unadj and Adj	tumor size, stage, lymph node stage, WHO overall stage, positive resection margins, TI/S biomarker ratios	197 (186)	62	I: 55	ADC: 45	Current & Former: 94%	Median
	CD8							II: 24	SCC: 46	Never: 6%
	FoxP3							III: 21	LCC: 3
									Other: 6
Paulsen, 2015 [49]	CD8	Tumor	OS	Unadj	n/a	536	68	I: 48	ADC: 38	Current: 33%	score of ≥2
								II: 36	SCC: 54	Former: 64%
								III: 16		Never: 3%
Tian, 2015 [50]	CD3	Tumor	OS	Adj	histology, pT, pN, pTNM stage, CD3, CD8, IL-2	129	71	I: 37	ADC: 37	Smoker: 37%	H-score 0–4 = low expression
	CD8							II: 22	SCC: 47	Never: 63%	H-score 5–12 = high expression
								III: 40	Other: 16
Ameratunga, 2016 [51]	CD8	Stroma	OS & DFS	Unadj and Adj	age, sex, pneumonectomy status, nodal stage, histology, smoking status, CD8+, FOXP3, PD-L1	527 (509)	69	NR	ADC: 55	Heavy: 73%	Score of % of positively stained cells > 3
									SCC: 34	Light: 17%
									Other: 11	Never: 7%
										Unknown: 4%
Kinoshita, 2016 [52]	CD4	Tumor	OS & DFS	Unadj and Adj (for CD20 only)	age, diameter, nodal metastasis, pleural invasion, FoxP3/CD4	218	65	I: 71	ADC: 72	Ever: 61%	CD4: > 251 cells/mm2
	CD8							II: 11	SCC: 18	Never: 39%	CD8: > 110 cells/mm2
	CD20							III: 18	Other: 10		FoxP3: > 163 cells/mm2
	FoxP3
Parra, 2016 [53]	CD57	Peri and intra tumor	OS (CD57 only) & DFS	Adj	tumor stage, adjuvant therapy	254	55	I: 50	ADC: 57	Current: 45%	Median
	CD68							II: 30	SCC: 43	Former: 30%
								III: 20		Never: 17%
Teng, 2016 [54]	CD8	Tumor	OS & DFS	Unadj and Adj (for CD8 only)	high risk, FoxP3/CD8+TILS	126	67	I: 100	ADC: 45	NR	CD3: ≥30%
	FoxP3								SCC: 33		FoxP3: ≥45/HPF
									LCC: 22
Uso, 2016 [55]	CD8	Tumor and Stroma	OS & DFS (CD8 only)	Adj	Tumoral CD8+ cells, Stromal FOXP3+ cells, FOXP3+ stroma/CD4+ tumor, FOXP3+ stroma/CD8+ tumor	122	85	I: 59	ADC: 42	Current: 48%	CD8: Median
	FoxP3							II: 21	SCC: 47	Former: 38%	FoxP3: >10%
								III: 20	Other: 11	Never: 14%
Yang, 2016 [56]	CD4	Tumor and Stroma	OS & DFS	Adj	age, sex, stage, PD-L1 expression, Stromal CD4+ T cell, Stromal regulatory T cell, Epithelial CD8+ T cell	105	85	I: 100	SCC: 100	Smoker: 75%	CD4: ≥5% (tumor) & ≥25% (stroma)
	CD8									Non Smoker: 25%	CD8: ≥5% (tumor) & ≥50% (stroma)
	FoxP3										FoxP3: ≥20%
Yazdi, 2016 [57]	CD8	Stroma	OS	Unadj and Adj	stage, sex, HLA-E	197	50	I: 31	ADC: 100	NR	mean
								II: 38
								III: 18
Huang, 2017 [58]	CD8	Tumor	OS & DFS	Adj	age, sex, differentiation, histology, tumor, node, CD133, OCT-4, CD8, CD56, HLA I, PD-L1	172	65	NR	SCC: 42	Smoker: 55%	Median
	CD56								Other: 58	Non-Smoker: 45%
Kinoshita, 2017 [59]	CD4	Tumor	OS & DFS	Unadj	n/a	164	5	I: 68	ADC: 100	Non-smoker: 100%	CD4: > 585 cells/mm2
	CD8							II, III: 32			CD8: > 900 cells/mm2
	CD20										CD20: > 1070 cells/mm2
	FoxP3										FoxP3: > 81 cells/mm2
Koh, 2017 [60]	CD8 [60]	Tumor	DFS	Adj	age, lymph node metastasis, CD103+ TILs, CD8+ TILs	378	95	I, II: 79	SCC: 100	Smoker: 90%	mean
								III: 21		Never: 5%
Sepesi, 2017 [61]	CD3	Intratumor	OS	Unadj and Adj	age, histology, thoracotomy, adjuvant therapy, DLCO %, Zubrod score, Lobectomy, PD-L1 tumor H score, PD-L1 tumor % expression, PD-L1 macrophages H score, PD-L1 macrophages % expression, CD3, CD4, CD8, CD45RO, CD57, CD68, FoxP3, PD-1	113	49	I: 100	ADC: 70	Current and Former: 89%	CD3: ≥ 827.3
	CD4								SCC: 30	Never: 11%	CD4 ≥ 852.1
	CD8										CD8: ≥ 292.3
	CD57										CD57: ≥ 1408
	CD68										CD68: ≥ 515
	FoxP3										FoxP3: ≥ 460.3
Ye, 2017 [62]	CD8		OS	Unadj and Adj	gender, age, pathology grade, T, N, clinical stage	102	54	I: 20	ADC: 100	NR	final score>6
								II: 59
								III: 21
Barua, 2018 [63]	CD4	Tumor	OS	Unadj	n/a	120	43	I: 4	ADC: 60	Current & Former: 83%	G-cross signatures
	CD8							II: 11	SCC: 29	Never: 17%
	CD68							III: 85	Other: 11
	FoxP3
Jackute, 2018 [64]	CD68	Total	OS	Unadj	n/a	80	80	I: 29	ADC: 48	Smoker: 83%	Median
								II: 33	SCC: 45	Non-Smoker: 17%
								III: 39	Other: 7
Mazzaschi, 2018 [65]	CD3		OS & DFS	Unadj	n/a	100	4	I: 35	ADC: 42	Current: 37%	Continuous
	CD8							II: 41	SCC: 58	Former: 52%
								III: 24		Never: 9%
Su, 2018 [66]	CD8	Stroma	OS & DFS	Adj	age, sex, smoking, micropapillary pattern, CEA, PD-1 expression	223	54	I: 67	ADC: 100	Smoker: 30%	≥25%
								II: 14		Non-Smoker: 70%
								III: 19
Matsubara, 2019 [67]	CD3	Tumor Margins	OS	Unadj and Adj	age, stage, lymphatic invasion, PD-L1, PD-L2	211	83	I: 54	SCC: 100	Smokers ≥30 pack-years: 85%	Median
	CD4							II: 35
	CD8							III: 11
Meng, 2019 [68]	CD4^	Tumor (FoxP3 Stroma)	OS	Adj	T, N, CD4, CD8, FoxP3, PD-L1	197	65	I: 57	SCC: 43	Smoker: 58%	CD4: 5%
	CD8*							II: 45	Non-SCC: 57	Non-smoker: 42%	CD8: 5%
	FoxP3*							III: 40			FoxP3: 20%
Cao, 2019 [69]	CD68	Tumor	OS & DFS	Unadj and Adj	T, N	137	56	I: 34	ADC: 68	Former or Current: 48%	Median
								II: 34	SCC: 32	Never: 53%
								III: 33

^*HR reported for stage I only. ^{^}HR reported for squamous only.

^#ADC = adenocarcinoma, SCC = squamous cell carcinoma, LCC = large cell carcinoma, NR = not reported, n/a = not applicable, OS = overall survival, DFS = disease free survival, Adj = adjusted for covariates, Unadj = not adjusted for covariates.

CD3

The CD3 complex is a defining feature of T cell lineage. T lymphocytes are an important aspect of adaptive immunity, subsets of which accomplish various immune system functions [14]. Ten eligible studies (n = 1,385 patients) reported on CD3+T lymphocytes as a prognostic biomarker. Of these, 7 articles reported unadjusted OS HRs (n = 1,108), 5 articles unadjusted DFS HR information (n = 415), and 3 articles reported adjusted OS HRs (n = 428); no articles reported an adjusted DFS HR estimate.

CD4

CD4+ T helper cells function as part of the adaptive immune system, where they are involved in priming the immune response through their interaction with major histocompatibility complex (MHC) class II proteins on antigen presenting cells (APCs). Upon activation, they differentiate and release a variety of cytokines to help B cells make antibodies, induce macrophage activity, and recruit neutrophils, eosinophils and basophils [15]. There were 10 studies (n = 1,496 patients) on CD4+ T helper cells as a prognostic maker in NSCLC patients. Five articles reported unadjusted OS HR (n = 1,037) and 3 articles reported unadjusted DFS (n = 466). Six articles reported an adjusted OS HR (n = 699), while only 1 article gave an adjusted DFS estimate (n = 105).

CD8

CD8+ cytotoxic T cells recognize antigen peptides presented on MHC class I molecules and are critical for immune defense against intracellular pathogens (i.e., viruses, bacteria) and for tumor surveillance [16]. There were 26 studies (n = 5,624 patients) that provided HR estimates for the association between CD8+ cytotoxic T cells and survival in stage I-III NSCLC patients. Fifteen of these articles gave unadjusted OS estimates (n = 3,317) and 8 reported on unadjusted DFS (n = 1,603). Thirteen articles gave an adjusted OS HR estimate (n = 2,456), while 9 reported an adjusted DFS HR (n = 2,596).

CD20

CD20 is expressed on the surface of B lymphocytes, which provide humoral immunity through the production of antibodies. The specific function of tumor infiltrating B cells is still being explored, though it’s suggested that these cells play a role in anti-tumor immunity through either directly presenting antigens to T cells or by generating tumor antigen-specific antibodies that form immune complexes with tumor antigens that are presented by professional APCs [17]. Seven eligible studies (n = 1,058) provided HRs on CD20+ B cells as a prognostic biomarker in NSCLC; 5 reported unadjusted OS (n = 897) and 4 reported unadjusted DFS (n = 540). Only one article reported an estimate for adjusted OS (n = 113) or adjusted DFS (n = 218).

FoxP3

T regulatory cells expressing FoxP3 are important for maintaining immune homeostasis by suppressing the proliferation and activation of cytotoxic T cell response [18]. There were 14 eligible (n = 2,464 patients) articles that reported on the predictive value of FoxP3+ T regulatory cells for survival; 9 articles included unadjusted OS (n = 1,547), 3 reported unadjusted DFS (n = 5080), 7 reported adjusted OS (n = 934) and 3 reported unadjusted DFS (n = 683).

CD56 and 57 (NK Cells)

CD56/57 expressing NK cells, a part of the innate immune system, are capable of quickly recognizing and killing both infected and malignant cells that have lost their MHC I receptors, without priming or prior activation [19]. Six articles reported on survival measures for CD56+ & CD57+ NK cells (n = 622). Three of these reported unadjusted OS (n = 258) and 1 reported unadjusted DFS (n = 84). Four articles reported adjusted OS (n = 443) with 2 reporting adjusted DFS (n = 280).

CD68

TAMS (as identified by CD68 antibodies) are thought to be driven by immunosuppressive cytokines such as IL-10 and TGF-beta and have been associated with suppressing T cell tumor response and promoting tumor growth and spread [20]. There were 12 eligible articles (n = 1,699); 6 reported unadjusted OS (n = 922), 1 study reported on unadjusted DFS (n = 137). Seven articles reported adjusted OS (n = 716), and 6 reported adjusted DFS (n = 861).

MCs

When activated, MCs release inflammatory mediators, increasing vascular permeability and recruiting other immune cells [21]. Their role in cancer development, however, is largely unknown. There were 4 unique eligible articles (n = 771 patients) that investigated OS according to MC infiltration. Three unique articles reported an unadjusted OS (n = 591), with only one study reporting unadjusted DFS (n = 201). Four articles reported adjusted OS (n = 733), while 1 study reported adjusted DFS (n = 201).

Immuno-infiltrates and survival

There were no significant differences between high vs. low CD3+ T cells or CD4+ T Helper cell infiltration in terms of OS or DFS (Table 2; Supplementary Figures 1 and 2). Adjusted data showed that increased levels of CD8+ T cytotoxic cells were predictive of improved OS (HR_{OS (adj)} = 0.68; 95% CI, 0.50–0.93) and DFS (HR_{DFS (adj)} = 0.60; 95% CI, 0.41–0.87; Table 2; Supplementary Figure 3A). However, CD8+ T cells were not significantly associated with improved survival (OS and DFS) in unadjusted analysis. Both unadjusted OS and DFS estimates indicated an association between CD20+ B cells and better survival (HR_{OS (unadj)} = 0.45; 95% CI, 0.22–0.93 and HR_{DFS (unadj)} = 0.57; 95% CI, 0.33–1.00, respectively). There was only one study reporting adjusted OS, which confirmed the finding (HR_{OS (adj)} = 0.16; 95% CI, 0.04–0.64; Table 2; Supplementary Figure 4A). High levels of FoxP3+ Treg lymphocytes were associated with worse OS when looking at both unadjusted (HR_{OS (unadj)} = 1.78; 95% CI, 1.20–2.64) and adjusted (HR_{OS (adj)} = 2.22; 95% CI, 1.47–3.34) data, and with worse DFS (HR_{DFS (adj)} = 2.07; 95% CI, 1.10–3.90) when adjusting for covariates (Table 2; Supplementary Figure 5A). Patients with a higher level of NK cells (CD56/CD57+) had significantly better adjusted OS than those with lower levels (HR_{OS (adj)} = 0.50; 95% CI, 0.26–0.95), although this association was not significant in the meta-estimate of unadjusted studies (Table 2; Supplementary Figure 6A). Increased presence of macrophages (CD68+) did not show a statistically significant prognostic value in terms of OS, although CD68+ macrophages were associated with worse adjusted DFS (HR_{DFS (adj)} = 1.84; 95% CI, 1.02–3.34), (Table 2; Supplementary Figure 7A). Higher levels of MCs were indicative of worse OS in the analysis of adjusted studies (HR_{OS (adj)} = 2.13; 95% CI, 1.14–3.96; Table 2; Supplementary Figure 8A).

Table 2: Unadjusted and adjusted pooled survival estimates

Overall Survival (OS)
Biomarker	Unadjusted HR (95% CI)	Number of Articles	n	Heterogeneity	Adjusted HR (95% CI)	Number of Articles	n	Heterogeneity
CD3	0.99 (0.87, 1.12)	7	1,108	I2 = 0.00%	0.71 (0.37, 1.37)	3	428	I2 = 0.00%
				Q = 1.64 (p = 0.90)				Q = 1.76 (p = 0.42)
CD4	0.76 (0.44, 1.32)	5	1,037	I2 = 34.65%	1.00 (0.55, 1.81)	6	699	I2 = 33.79%
				Q = 3.47 (p = 0.19)				Q = 7.55 (p = 0.18)
CD8	0.99 (0.95, 1.04)	15	3,317	I2 = 0.00%	0.68 (0.50, 0.93)	13	2,456	I2 = 23.78%
				Q = 10.96 (p = 0.69)				Q = 15.74 (p = 0.20)
CD20	0.45 (0.22, 0.93)	5	861	I2 = 52.29%	0.16 (0.04, 0.64)	1	113	–
				Q = 8.38 (p = 0.08)
FoxP3	1.78 (1.20, 2.64)	9	1,547	I2 = 14.93%	2.22 (1.47, 3.34)	7	934	I2 = 0.00%
				Q = 9.40 (p = 0.31)				Q = 4.74 (p = 0.58)
CD56/CD57	0.66 (0.35, 1.25)	3	258	I2 = 13.63%	0.50 (0.26, 0.95)	4	443	I2 = 0.00%
				Q = 2.32 (p = 0.31)				Q = 0.22 (p = 0.97)
CD68	1.36 (0.75, 2.45)	6#	922	I2 = 40.21%	1. 13 (0.77, 1.65)	7	716	I2 = 0.00%
				Q = 10.03 (p = 0.12)				Q = 4.18 (p = 0.65)
MCs	1.81 (1.01, 3.15)	3^	462	I2 = 0.00%	2.01 (0.88–3.92)	3	604	I2 = 43.99%
				Q = 0.09 (p = 0.90)				Q = 5.36 (p = 0.15)

Disease Free Survival (DFS)
Biomarker	Unadjusted HR (95% CI)	Number of Articles	n	Heterogeneity	Adjusted HR (95% CI)	Number of Articles	n	Heterogeneity
CD3	0.98 (0.86, 1.12)	5	415	I2 = 0.00%	–	–	–	–
				Q = 2.63 (p = 0.62)
CD4	0.72 (0.39, 1.32)	3	466	I2 = 0.00%	0.59 (0.11, 3.12)	1	105	–
				Q = 0.60 (p = 0.74)
CD8	0.94 (0.77, 1.16)	8	1,630	I2 = 16.65%	0.60 (0.41, 0.87)	9	2,596	I2 = 20.43%
				Q = 8.39 (p = 0.30)				Q = 10.05 (p = 0.26)
CD20	0.57 (0.33, 1.00)	4	540	I2 = 0.00%	0.51 (0.20, 1.32)	1	218	–
				Q = 0.27 (p = 0.97)
FoxP3	1.45 (0.81, 2.59)	3	508	I2 = 0.00%	2.07 (1.10, 3.90)	3	683	I2 = 0.00%
				Q = 0.07 (p = 0.96)				Q = 0.02 (p = 0.99)
CD56/CD57	1.35 (0.39, 4.66)	1	84	–	0.59 (0.27, 1.28)	2	280	I2 = 0.00%
								Q = 0.34 (p = 0.56)
CD68	1.81 (0.66, 4.92	1	137	–	1.84 (1.02, 3.34)	6#	861	I2 = 29.96%
								Q = 7.14 (p = 0.21)
MCs	2.30 (1.20–4.70)	1	201	–	1.50 (0.60, 3.60)	1	201	–

^#6 unique articles but 7 estimates. For Parra, 2016 adenocarcinoma and squamous histologies were evaluated separately. ^{^}3 unique articles but 4 estimates.

HR >1 refers to high level of immune cell subtype associated with worse survival.

DISCUSSION

To our knowledge, this is the first meta-analysis demonstrating immune cell subpopulations within the TME having prognostic value in early-stage NSCLC patients undergoing surgical resection. Furthermore, this meta-analysis is also the first to present findings according to whether or not studies adjusted for clinical covariates. In the analysis of unadjusted studies, only CD20+ B cells were associated with improved OS and DFS, while FOXP3+ Tregs were associated with worse OS. However, when evaluating studies that adjusted for important clinical covariates (such as stage and age), higher levels of CD8+ cytotoxic T cells were associated with improved OS and DFS, and increased CD20+ B cells and CD 56/57+ NK cells were associated with improved OS. Lung cancers with increased FoxP3+ T regulatory cells or increased MCs had worse OS, and cancers with increased CD68+ macrophages had worse DFS. Our results are in keeping with what is known about the function of these immune cells; we now show these cells appear to have an important role in clinical outcomes of early-stage lung cancer.

Data evaluating the role of immune cells of the TME and lung cancer outcomes mostly come from small case series that are limited by sample size. Three prior meta-analyses attempting to combine the data from these studies have reported improved outcomes with increased CD8+ T cells and worse outcomes associated with FOXP3+ T regulatory cells [8–10]. However, the impact of CD20+ B cells, CD3+ T cells, CD4+ T cells and TAMs on survival was not consistent across studies. These meta-analyses presented significant heterogeneity by including data from both early- and advanced-stage NSCLC patients, possibly conferring different stages of the immunoediting process. Furthermore, patients with metastatic disease are subject to different treatments strategies (i.e., systemic chemotherapies, immunotherapies, targeted treatments) that can affect both the immune composition of the TME and survival outcomes. Additionally, tissue samples from metastatic disease are typically biopsies (i.e., core or fine needle aspirates), and it is unclear whether immunophenotyping of these samples completely characterizes the TME. Our study overcomes several of these limitations by focusing on the immune cell composition of resected early-stage lung cancers, enabling us to evaluate studies using 1) a more homogenous NSCLC population that is likely treatment naïve (except possibly the small percentage receiving neo-adjuvant chemotherapy) and 2) larger surgical samples, to study the immune composition of the TME in a more comprehensive way. Additionally, given that surgical resection of early-stage lung cancer represents the best chance for cure and long-term survival, the analysis of this very important subgroup is an important and innovative addition to the existing literature.

This meta-analysis is also the first to assess immune cell biomarkers in lung cancer while stratifying studies based on their adjustment for clinically relevant covariates. Patient factors such as sex and stage may influence both the immune composition of the TME and post-surgical survival, and therefore may distort the true relationship between the TME and survival. Our study in fact showed differences in the adjusted and unadjusted estimates of CD8+ cytotoxic T cells, NK cells, CD68+ macrophages, MCs and survival. Being that different studies were used in the unadjusted and adjusted analyses it is possible that this finding could be the result of differing methodological approaches instead of confounding. However, when we assessed direct comparisons between unadjusted and adjusted findings for a single study and biomarker, adjustment for clinical covariates impacted statistical significance in 41% of instances (Supplementary Table 1). Our ability to demonstrate significantly variable findings based on the adjustment of important clinical confounders highlights the need for future research to account for several clinical factors in order to determine the independent association of immune cell biomarkers and lung cancer survival, both overall and on a per patient basis. A consensus on standard protocol for these studies is especially needed as more immune cell subtypes in relation to cancer continues to be explored.

Our study has several limitations worth discussing. As with any meta-analysis, an important limitation is the possibility of publication bias, although the funnel plots were not suggestive of publication bias, with the exception of CD8+ T cells (Supplementary Figures 3B–3E, 4B–4C, 5B–5E, 6B–6C, 7B–7D, and 8B–8C). By including only stage I-IIIA NSCLC that underwent surgical resection we minimized heterogeneity. We were also unable to directly adjust for important clinical confounders of the relationships between immune cell biomarkers and survival. While early-stage patients undergoing surgical resection are typically treatment näive, in this meta-analysis we were not able to verify this. We did stratify the results according to the presence of adjustment for important clinical confounders, but this was still a limited approach as the included studies adjusted for a variety of different clinical covariates. For instance, only 27% of eligible articles adjusted for stage, 22% for sex, 16% for histology and just 7% for smoking. Available studies also used different cutoffs to separate high versus low biomarker infiltration, introducing additional inter-study variability. Despite these possible sources of variability, there was not statistically significant heterogeneity among studies as determined by the I² statistic. While we attempted to investigate multiple relevant immune cell types, for some immune cell biomarkers, such as CD3+ T cells, CD20+ B cells, NK cells, TAMs and MCs, there was limited literature, thus comparisons between unadjusted and adjusted OS and DFS estimates was not possible.

Despite these limitations, by summarizing the results from several studies, we were able to overcome the limitation of small sample size observed in individual studies evaluating multiple immune cell biomarkers, and in doing so come to a more accurate estimate of their prognostic value. This meta-analysis is also the first to investigate if these immune cells are biomarkers for survival in surgically-resectable NSCLC patients, an important subgroup of patients that are likely growing as the result of screening guidelines. It is also the first to take clinical covariates into account.

Our results suggest that there are immune cells infiltrating the TME that can be considered biomarkers of survival in early-stage NSCLC. Specifically, we demonstrate the significant association of CD8+ cytotoxic T cells, CD20+B cells, NK cells and FoxP3+ T regulatory cells with survival in cases of early, resectable disease. The presence or absence of these immune cells within the TME of resected lung cancers may be used to stratify patients according to risk of recurrence and survival, with implications for who may be more likely to benefit from neo-adjuvant or adjuvant therapies [22]. Currently, there is ongoing debate about the effectiveness of immunotherapy treatment for early-stage lung cancers, with a concern that the delay in surgery to allow for immunotherapy neo-adjuvant treatment is more detrimental than beneficial to overall outcomes. However, if the immune infiltration profiles of those most likely to have disease recurrence after surgery can be identified, a patient’s immune infiltration can be used to determine which patients will receive the most benefit from either adjuvant or neo-adjuvant immunotherapy.

In conclusion, future randomized clinical trials should verify the prognostic value of CD8+ cytotoxic T cells, CD20+B cells, NK cells, CD68+ macrophages and FoxP3+ T regulatory cells and other immune components of the TME. Our meta-analysis was limited by what is published in the existing literature, but future studies should also attempt to determine the prognostic value of other immune cells that were beyond the scope of this work, including dendritic cells, Th1 cells, Th2 cells, Th17 cells and eosinophils, all of which may have an important impact on tumor immune escape and tumorigenesis, but at present are understudied [23, 24]. This will allow for clinicians to consider multi-immune cell panels when evaluating potential cancer outcomes and the need for additional treatment of patients with early-stage, surgically resected NSCLC.

Author contributions

ET, RV and RF conceptualized study design, ST and WLC performed data abstraction and statistical analysis, SG, PW and FP contributed additional expertise. All authors participated in manuscript preparation and have approved this final version.

CONFLICTS OF INTEREST

RV has received consulting honorarium from Onconova Pharmaceuticals. SG has received research support from Immune Design Corp.

FUNDING

RV was partially supported by internal grant KL02539874. FP and PW are supported by grant U24 CA210993, from the National Cancer Institute Clinical Proteomic Tumor Analysis Consortium (CPTAC).