Description of studies
In our meta-analysis, we searched for 350 relevant studies in the databases. 35 duplicated data from the same population were removed. Due to lack of statistical analysis and clinical data (N=265), reviews (N=7), case reports (N=5), editorial (N=15) and abstract (N=11), 303 articles were excluded. Finally, 12 articles were included [15-26] (Fig. 1). We listed the 12 articles in Table 1. Overall, the 12 articles included 1,583 GC patients were published between 2013 and 2021. 2 studies came from Korea, 1 study was from Japan, 1 study was from Brazil, 8 studies were from China.
Table 1. Clincal data included in the literature
No
|
First
author
|
Year
|
No
|
Gender (M/F)
|
pNstage
(N+/N-)
|
During
|
Country
|
Method
|
NOS
|
1
|
Wang
|
2021
|
115
|
90/25
|
79/36
|
2016-2018
|
China
|
IHC
|
7
|
2
|
Ouyang
|
2019
|
86
|
48/38
|
51/35
|
2014-2017
|
China
|
IHC
|
6
|
3
|
Neves
|
2018
|
33
|
20/13
|
-
|
-
|
Brazil
|
IHC
|
6
|
4
|
Deng
|
2013
|
90
|
64/26
|
-
|
2006-2007
|
China
|
IHC
|
6
|
5
|
Suh
|
2015
|
117
|
-
|
69/48
|
2006-2009
|
Korea
|
IHC
|
7
|
6
|
Yoshii
|
2012
|
92
|
44/48
|
34/58
|
2000-2010
|
Japan
|
IHC
|
7
|
7
|
Guo
|
2020
|
268
|
210/58
|
183/85
|
2009-2011
|
China
|
IHC
|
7
|
8
|
Geng
|
2014
|
100
|
61/39
|
68/32
|
2005-2006
|
China
|
IHC
|
6
|
9
|
Wang
|
2020
|
106
|
67/39
|
21/85
|
2005-2015
|
China
|
IHC
|
8
|
10
|
Li
|
2017
|
197
|
128/69
|
102/95
|
2007-2010
|
China
|
IHC
|
8
|
11
|
Ma
|
2014
|
197
|
130/67
|
86/101
|
-
|
China
|
INC
|
6
|
12
|
Won
|
2017
|
182
|
-
|
82/100
|
2006-2012
|
Korea
|
IHC
|
7
|
M,male;F,female;-,no data.
Foxp3 expression and clinicopathological data
To study the association between Foxp3 expression and clinicopathological parameters, we observed the association between Foxp3 expression and age, sex, tumor lymph node metastasis, tumor differentiation, TNM stage, and tumor size by meta-analysis. In our results, we found that Foxp3 expression was not related to age (OR = 0.899, 95% CI = 0.549-1.472, p = 0.672), gender (OR = 0.836, 95% CI = 0.647-1.080, p = 0.171), differentiation (OR = 1.083, 95% CI = 0.616-1.904, p = 0.781), and TNM stage (OR = 0.967, 95% CI = 0.583-1.605, p = 0.898) in GC patients. We observed that Foxp3 expression was related to tumor size (OR = 0.489, 95% CI = 0.2857-0.839, p = 0.009) and lymph node metastasis (OR = 0.514, 95% CI = 0.280-0.944, p = 0.032) in GC patients (Table 2 and Fig. 2).
Table 2. Results for the meta-analysis between Foxp3 expression and clinicopathological parameters of GC patients
Clinical parameters
|
No. of studies
|
Overall
OR (95%CI)
|
Heterogeneity test (Q, I2, P)
|
Age
|
14,15,16,18,19,20,22
|
0.899(0.549-1.472)
|
16.62,63%,0.2632 (random-effect)
|
Sex
|
13,14,15,16,18,19,20,22,23
|
0.836(0.647-1.080)
|
6.98,0.00%,0.538 (random-effect)
|
Lymph node metastasis
(N0 vs N1-N3)
|
12,13,14,15,16,18,20,21,22,23
|
0.514(0.280-0.944)
|
63.20,85.8%,0.000 (random-effect)
|
Differentiation (Well/moderate vs poor)
|
12,14,15,17,18,19,20,21
|
1.083(0.616-1.904)
|
32.97,78.8%,0.000 (random-effect)
|
TNM stage (I II vs III IV)
|
12,14,15,16,17,18,22,23
|
0.967(0.583-1.605)
|
35.17,77.73%,0.000 (random-effect)
|
tumor size (≥5cm vs<5cm)
|
14,15,18,20
|
0.489(0.285-0.839)
|
6.06,50.05%,0.109 (random-effect)
|
Abbrebiation:OR,odds rates.
Publication bias between Foxp3 expression and clinicopathological parameters
We did not find any obvious publication bias in age, sex, tumor differentiation, tumor lymph node metastasis, TNM stage, and tumor size. p value was more than 0.05 in Begg’s test (p = 0.634, p = 0.054, p = 0.954, p = 0.954, p = 0.098, and p = 0.780, respectively) (Fig. 3).
Foxp3 expression and OS
To study the association between Foxp3 expression and OS of GC patients, we found that the association between Foxp3 expression and the OS of GC patients in Table 3 and Fig. 4. Compared with lower Foxp3 expression, high Foxp3 expression had a better 1-year OS (OR = 0.606, 95% CI = 0.399-0.921, p = 0.019), 3-year OS (OR = 0.629, 95% CI = 0.436-0.908, p = 0.013), and 5-year OS (OR = 0.412, 95% CI = 0.231-0.734, p = 0.003) of GC patients.
Publication bias between Foxp3 expression and OS
As shown in Fig. 5, there is no publication bias between Foxp3 expression and 1-year, 3-year, and 5-year OS (Bgger’s test, p = 0.637, p = 0.275, and p = 0.236) of GC patients.
Table3. Patient OS in relationship with expression of Foxp3
Author(year)
|
|
1-year OS
|
|
3-year OS
|
|
5-year OS
|
|
OR (95%CI)
|
Weight
|
OR (95%CI)
|
Weight
|
OR (95%CI)
|
|
Weight
|
Wang (2021)
|
|
0.39(0.04,3.82)
|
3.32%
|
|
0.47(0.09,2.43)
|
4.30%
|
|
-
|
|
-
|
Guo (2020)
|
|
1.07(0.43,2.67)
|
20.88%
|
|
0.92(0.52,1.64)
|
18.50%
|
|
0.66(0.37,1.16)
|
|
38.07%
|
Suh (2015)
|
|
0.45(0.17,1.16)
|
19.11%
|
|
0.56(0.24,1.29)
|
12.24%
|
|
-
|
|
-
|
Li (2017)
|
|
0.14(0.03,0.73)
|
6.56%
|
|
1.02(0.48,2.17)
|
13.84%
|
|
-
|
|
-
|
WON (2015)
|
|
0.43(0.13,1.45)
|
11.82%
|
|
0.43(0.18,1.00)
|
11.94%
|
|
0.23(0.10,0.52)
|
|
27.50%
|
Ouyang (2019)
|
|
0.48(0.09,2.63)
|
6.04%
|
|
-
|
-
|
|
-
|
|
-
|
Yoshii (2012)
|
|
1.52(0.34,6.75)
|
7.82%
|
|
1.01(0.43,2.38)
|
11.83%
|
|
-
|
|
-
|
Ma (2014)
|
|
0.69(0.30,1.62)
|
24.45%
|
|
0.65(0.34,1.24)
|
16.50%
|
|
0.39(0.23,0.73)
|
|
34.43%
|
Overall
|
|
0.61(0.40,0.92)
|
100%
|
|
0.63(0.44,0.91)
|
100%
|
|
0.41(0.23,0.73)
|
|
100%
|
Relative risks refer to no relationship in survival for patients with expression of Foxp3. One year: heterogeneity χ2 = 6.97; p = 0.432; I2 = 0.0%; z-test = 2.35; p = 0.019; Three years: heterogeneity χ2 = 11.30; p = 0.126; I2 = 38.1%; z-test = 2.48; p = 0.013; Five years: heterogeneity χ2 = 4.48; p = 0.107; I2 = 55.3%; z-test = 3.00; p = 0.003.
TCGA
By TCGA data, we explored the association between Foxp3 expression and prognosis in GC patients. The Gene Expression Profiling Interactive Analysis (GEPIA) included 384 patients who were equally divided into low Foxp3 group and high Foxp3 group, with 192 patients in each group (Fig.6). We found that the OS of the high Foxp3 expression group was higher than the low Foxp3 expression group, but it was not statistically significant (p=0.34).