Search and selection results
The initial search yielded a total of 9470 studies. After removal of duplicates, a total of 7659 studies were screened by analysis of title and abstract, and 7435 studies were removed because they met one or more exclusion criteria. The leaving 224 studies were then assessed for eligibility by full-text examination, and a further 199 were excluded for ineligibility. Reasons for exclusion were recorded. Finally, 25 studies[3, 8, 15-18, 25-43] were included in the final analysis (Fig. 1).
Study characteristics
Among the 25 included studies, 7 had a multicentre setting and 18 had a single centre design. Five of the included studies were prospectively performed; the remaining twenty were retrospective. 22 studies were considered of high quality (N-O score ≥ 6), and 3 studies were considered of low quality. Complete characteristics of the included studies are available in Table 1.
Factors not included in the quantitative synthesis
Six of clinicopathological and molecular factors could not be included in quantitative synthesis because they had only a single study of their subgroup, or their methodology did not permit pooling data. The six factors were serum CEA[25], serum CA125[32], CTGF (connective tissue growth factor)[40], DDR2 (discoidindomain receptor 2)[30], VIM (vimentin)[42], and TP53[28] respectively. We included these factors in table 1 for completeness, but not in the final quantitative synthesis through meta-analysis.
Finally, 21 studies about 13 factors were included in the quantitative synthesis through meta-analysis, 7 studies on gender, 4 studies on tumor invasion depth, 3 studies on lymph node metastasis, 5 studies on differentiation, 6 studies on primary tumor site, 7 studies on histology, 2 studies on serum CA19-9, 2 studies on PROK1/PROKR2, 9 studies on BRAF, 6 studies on KRAS, 2 studies on NRAS, 2 studies on PIK3CA and 4 studies on MSI-H/dMMR status.
Gender
Seven studies[8, 15, 17, 18, 25, 33, 41], including data on 160679 patients (30366 synchronous pmCRC, 130313 non-pmCRC) regarding gender, were included for eligibility in the meta-analysis. The pooled analysis indicated that female was associated positively with synchronous CPM compared with male (OR 1.299; 95% CI, 1.118 to 1.509; P = 0.001) (Fig. 2a). There was significant heterogeneity (Cochran Q, P < 0.001; I² = 76.9 percent). In order to explore possible sources of heterogeneity, sensibility analysis was performed by omitting each study sequentially to test the influence of each individual study on the pooled result. When one study[17] was omitted, there was no significant heterogeneity (Cochran Q, P = 0.099; I² = 46.0 percent), with no noticeable influence on the pooled OR and confidence interval. It’s noteworthy that the rate of female in the synchronous CPM group was > 50 percent in that one study, but the others were < 50 percent.
Tumor invasion depth
Four studies[3, 15, 25, 41], including data on 19432 patients (809 synchronous pmCRC, 18623 non-pmCRC) regarding tumor invasion depth, were included for eligibility in the meta-analysis. The pooled analysis indicated that T4 was associated positively with synchronous CPM compared with T1-3 (OR 12.331; 95% CI, 7.734 to 19.660; P < 0.001) (Fig. 2b). There was significant heterogeneity (Cochran Q, P = 0.009; I² = 74.2 percent). When one study[41] was omitted through sensibility analysis, there was no significant heterogeneity (Cochran Q, P = 0.593; I² = 0 percent), with no noticeable influence on the pooled result.
Lymph node metastasis
Three studies[3, 15, 41], including data on 16097 patients (702 synchronous pmCRC, 15395 non-pmCRC) comparing lymph node metastasis, were included for eligibility in the meta-analysis. The pooled analysis indicated that N1-2 was associated positively with synchronous PM compared with N0 (OR 5.665; 95% CI, 3.628 to 8.848; P < 0.001) (Fig. 2c). There was significant heterogeneity (Cochran Q, P = 0.068; I² = 62.7 percent). When one study[3] was omitted through sensibility analysis, there was no significant heterogeneity (Cochran Q, P = 0.765; I² = 0 percent), with no noticeable influence on the pooled result.
Differentiation
Five studies[15, 17, 25, 41, 43], including data on 108360 patients (21986 synchronous pmCRC, 86374 non-pmCRC) comparing differentiation, were included for eligibility in the meta-analysis. The pooled analysis indicated that poorly differentiated grade was associated positively with synchronous CPM compared with well/moderately differentiated grade (OR 2.560; 95% CI, 1.537 to 4.265; P < 0.001) (Fig. 2d). There was significant heterogeneity (Cochran Q, P < 0.001; I² = 94.5 percent). When one study[17] was omitted through sensibility analysis, there was no significant heterogeneity (Cochran Q, P = 0.636; I² = 0 percent), with no noticeable influence on the pooled result.
Location of primary tumor
Six studies[3, 8, 15, 18, 25, 41] including data on 24331 patients (1610 synchronous pmCRC, 22721 non-pmCRC) regarding right-sided colon cancer were included for eligibility in the meta-analysis. Synchronous CPM was associated positively with right-sided colon cancer location (OR 2.468; 95% CI, 2.050 to 2.970; P < 0.001) (Fig. 3a). There was no significant heterogeneity (Cochran Q, P = 0.119; I² = 42.9 percent).
Six studies[3, 8, 15, 18, 25, 41] including data on 24331 patients (1610 synchronous pmCRC, 22721 non-pmCRC) regarding left-sided colon cancer were included for eligibility in the meta-analysis. Synchronous CPM was not associated with left-sided colon cancer location (OR 1.000; 95% CI, 0.761 to 1.314; P = 0.998) (Fig. 3b). There was significant heterogeneity (Cochran Q, P = 0.004; I² = 71.4 percent). When one study[15] was omitted through sensibility analysis, there was less significant heterogeneity (Cochran Q, P = 0.049; I² = 58.0 percent), with no noticeable influence on the pooled result.
Five studies[3, 8, 18, 25, 41] including data on 23278 patients (1519 synchronous pmCRC, 21759 non-pmCRC) regarding rectal cancer were included for eligibility in the meta-analysis. Synchronous CPM was associated negatively with rectal cancer location (OR 0.323; 95% CI, 0.284 to 0.368; P < 0.001) (Fig. 3c). There was no significant heterogeneity (Cochran Q, P = 0.969; I² = 0 percent).
Histology
Six studies[3, 26, 36, 38, 41, 43] including data on 24252 patients (1600 synchronous pmCRC, 22652 non-pmCRC) regarding non-mucinous adenocarcinoma (NMC) were included for eligibility in the meta-analysis. Synchronous CPM was associated negatively with NMC (OR 0.319; 95% CI, 0.237 to 0.429; P < 0.001) (Fig. 4a). There was significant heterogeneity (Cochran Q, P = 0.005; I² = 70.4 percent). When one study[36] was omitted through sensibility analysis, there was no significant heterogeneity (Cochran Q, P = 0.106; I² = 47.5 percent), with no noticeable influence on the pooled OR and confidence interval.
Seven studies[3, 17, 26, 36, 38, 41, 43] including data on 154377 patients (29448 synchronous pmCRC, 124929 non-pmCRC) regarding mucinous adenocarcinoma (MC) were included for eligibility in the meta-analysis. Synchronous CPM was associated positively with MC (OR 3.565; 95% CI, 2.095 to 6.064; P < 0.001) (Fig. 4b). There was significant heterogeneity (Cochran Q, P < 0.001; I² = 97.1 percent). In order to explore possible sources of heterogeneity, subgroup analysis was performed. Two studies[17, 36] were divided into the subgroup one that had no significant heterogeneity (Cochran Q, P = 0.228; I² = 31.2 percent), with no noticeable influence on the pooled result, and the others[3, 26, 38, 41, 43] were divided into the subgroup two that also had no significant heterogeneity (Cochran Q, P = 0.174; I² = 37.0 percent), with no noticeable influence on the pooled result. The subgroup one had much higher OR value in each study than subgroup two.
Three studies[3, 26, 38] including data on 5741 patients (673 synchronous pmCRC, 5068 non-pmCRC) regarding signet-ring cell carcinoma (SRCC) were included for eligibility in the meta-analysis. Synchronous CPM was associated positively with SRCC (OR 4.480; 95% CI, 1.836 to 10.933; P = 0.001) (Fig. 4c). There was significant heterogeneity (Cochran Q, P = 0.036; I² = 69.7 percent). When one study[3] was omitted through sensibility analysis, there was no significant heterogeneity (Cochran Q, P = 0.656; I² = 0 percent), with no noticeable influence on the pooled result. It’s noteworthy that the study had a much higher OR value.
Serum CA19-9
Two studies[25, 39], including data on 728 patients (24 synchronous pmCRC, 704 non-pmCRC) regarding serum CA19-9 status, were included for eligibility in the meta-analysis. Serum CA19-9 level of up to 37.0 u/ml was taken as upper cut-off values for reference ranges. Synchronous CPM was associated positively with elevated serum CA19-9 (OR 12.868; 95% CI, 5.196 to 31.867; P < 0.001) (Fig. 5a). There was no significant heterogeneity (Cochran Q, P = 0.710; I² = 0 percent).
PROK1/PROKR2
Two studies[34, 37], including data on 944 patients (29 synchronous pmCRC, 915 non- pmCRC) regarding PROK1/PROKR2 status, were included for eligibility in the meta-analysis. Synchronous CPM was associated positively with PROK1/PROKR2-positive (OR 2.244; 95% CI, 1.031 to 4.884; P = 0.042) (Fig. 5b). There was no significant heterogeneity (Cochran Q, P = 0.344; I² = 0 percent).
BRAF status
Nine studies[8, 16, 18, 27-29, 31, 33, 35], including data on 4979 patients (704 synchronous pmCRC, 4275 non-pmCRC) regarding BRAF status, were included for eligibility in the meta-analysis. Synchronous PM was associated positively with BRAF mutations (OR 2.586; 95% CI, 1.674 to 3.994; P < 0.001) (Fig. 5c). There was significant heterogeneity (Cochran Q, P = 0.019; I² = 56.3 percent). When one study[28] was omitted through sensibility analysis, there was no significant heterogeneity (Cochran Q, P = 0.073; I² = 45.9 percent), with no noticeable influence on the pooled OR and confidence interval. It’s clearly seen that the study has a smaller sample size.
KRAS status
Six studies[8, 16-18, 28, 33], including data on 134197 patients (28362 synchronous pmCRC, 105835 non-pmCRC) regarding KRAS status, were included for eligibility in the meta-analysis. Synchronous CPM was not associated with KRAS mutations (OR 0.972; 95% CI, 0.576 to 1.638; P = 0.914) (Fig. 6a). There was significant heterogeneity (Cochran Q, P < 0.001; I² = 92.4 percent). When one study[17] was omitted through sensibility analysis, there was no significant heterogeneity (Cochran Q, P = 0.774; I² = 0 percent), with no noticeable influence on the pooled result. It’s noteworthy that the rate of KRAS mutations in the synchronous PM group was much lower in the study.
NRAS status
Two studies[16, 28], including data on 731 patients (43 synchronous pmCRC, 688 non- pmCRC) regarding NRAS status, were included for eligibility in the meta-analysis. Synchronous CPM was not associated with NRAS mutations (OR 1.140; 95% CI, 0.133 to 9.748; P = 0.905) (Fig. 6b). There was no significant heterogeneity (Cochran Q, P = 0.373; I² = 0 percent).
PIK3CA status
Two studies[16, 33], including data on 897 patients (93 synchronous pmCRC, 804 non-pmCRC) regarding PIK3CA status, were included for eligibility in the meta-analysis. Synchronous CPM was not associated with PIK3CA mutations (OR 0.667; 95% CI, 0.289 to 1.540; P = 0.343) (Fig. 6c). There was no significant heterogeneity (Cochran Q, P = 0.415; I² = 0 percent).
MSI-H/dMMR status
Four studies[16-18, 31], including data on 131015 patients (27922 synchronous pmCRC, 103093 non-pmCRC) regarding MSI-H/dMMR status, were included for eligibility in the meta-analysis. Synchronous CPM was not associated with MSI-H/dMMR (OR 1.087; 95% CI, 0.351 to 3.367; P = 0.885) (Fig. 6d). There was significant heterogeneity (Cochran Q, P = 0.097; I² = 52.5 percent). When one study[18] was omitted through sensibility analysis, there was no significant heterogeneity (Cochran Q, P = 0.153; I² = 46.6 percent), with no noticeable influence on the pooled result.
Publication bias
No significant publication bias was found, according to visual inspection of funnel plot and to Begg’s test (supplementary Fig. S1-S5).