Oncotarget

Research Papers:

Oncofetal gene SALL4 and prognosis in cancer: A systematic review with meta-analysis

PDF |  HTML  |  Supplementary Files  |  How to cite  |  Order a Reprint

Oncotarget. 2017; 8:22968-22979. https://doi.org/10.18632/oncotarget.14952

Metrics: PDF 1234 views  |   HTML 1473 views  |   ?  

Lorenzo Nicole, Tiziana Sanavia, Nicola Veronese, Rocco Cappellesso, Claudio Luchini, Paolo Dabrilli and Ambrogio Fassina _

Abstract

Lorenzo Nicolè1,*, Tiziana Sanavia2,*, Nicola Veronese3, Rocco Cappellesso1,7, Claudio Luchini4,5,6, Paolo Dabrilli1, Ambrogio Fassina1

1Department of Medicine, Surgical Pathology & Cytopathology Unit, University of Padova, Padova, Italy

2Department of Biomedical Informatics, Harvard Medical School, Boston, MA, USA

3National Research Council, Neuroscience Institute, Aging Branch, Padova, Italy

4Department of Pathology and Diagnostics, University and Hospital Trust of Verona, Verona, Italy

5ARC-NET Research Center, University and Hospital Trust of Verona, Verona, Italy

6Department of Pathology, Santa Chiara Hospital, Trento, Italy

7Pathological Anatomy and Histology Unit, Veneto Institute of Oncology, Padova, Italy

*These authors have contributed equally to this work

Correspondence to:

Ambrogio Fassina, email: ambrogio.fassina@unipd.it

Keywords: SALL4, cancer, prognosis, meta-analysis, targeted therapy

Received: April 05, 2016    Accepted: January 06, 2017    Published: February 01, 2017

ABSTRACT

The Spalt-Like Transcription Factor 4 (SALL4) oncogene plays a central function in embryo-fetal development and is absent in differentiated tissues. Evidence suggests that it can be reactivated in several cancers worsening the prognosis. We aimed at investigating the risk associated with SALL4 reactivation for all-cause mortality and recurrence in cancer using the current literature. A PubMed and SCOPUS search until 1st September 2016 was performed, focusing on perspective studies reporting prognostic parameters in cancer data. In addition, 17 datasets of different cancer types from The Cancer Genome Atlas were considered. A total of 9,947 participants across 40 cohorts, followed-up for about 5 years on average, were analyzed comparing patients showing SALL4 presence (SALL4+, n = 1,811) or absence (SALL4-, n = 8,136). All data were summarised using risk ratios (RRs) for the number of deaths/recurrences and hazard ratios (HRs) for the time-dependent risk related to SALL4+, adjusted for potential confounders. SALL4+ significantly increased overall mortality (RR = 1.34, 95% confidence intervals (CI)=1.21-1.48, p<0.0001, I2=66%; HR=1.4; 95%CI: 1.19-1.65; p<0.0001; I2=63%) and recurrence of disease (RR = 1.25, 95% CI = 1.1-1.42, p=0.0006, I2=62%); HR=1.52; 95% CI: 1.22-1.89, p=0.0002; I2=69%) compared to SALL4-. Moreover, SALL4 remained significantly associated with poor prognosis even using HRs adjusted for potential confounders (overall mortality: HR=1.4; 95%CI: 1.19-1.65; p<0.0001; I2=63%; recurrence of disease: HR=1.52; 95% CI: 1.22-1.89, p=0.0002; I2=69%). These results suggest that SALL4 expression increases both mortality and recurrence of cancer, confirming this gene as an important prognostic marker and a potential target for personalized medicine.


INTRODUCTION

The stem-like phenotype in cancer is the result of epigenetic and genetic alterations leading to the expression of genes involved in cell migration, invasion, angiogenesis, self-renewal, anti-apoptosis, and immune-escape, which are fundamental for the embryo-fetal development. The expression of a stem-like phenotype seems to play a central role in defining the malignant potential of different cancers. During the last decades, several stemness-related genes have been proposed as diagnostic markers for cancer, sometimes with prognostic significance. In particular, the fetal oncogene Spalt-Like Transcription Factor 4 (SALL4) has recently emerged as a potential prognostic marker in many tumors. SALL4 encodes for a zinc-finger transcription factor that plays an essential role during embryo-fetal development forming a regulatory network with other stemness-related genes, such as the Octamer-Binding Transcription Factor 4 (OCT-4), the Nanog Homeobox (NANOG), and the Sex-Determining Region Y-Box 2 (SOX2), [13] and then gradually disappears until it remains strongly silenced in fully differentiated tissues (except for the germline cells and the hematopoietic stem/progenitor cells) [4, 5].

Analyses of SALL4 expression and its epigenetic status as well as studies on cellular models have shown its oncogenic role in several tumors, such as precursor B-cell lymphoblastic lymphoma, acute and chronic myeloid leukemia, gastrointestinal, breast, and lung cancers. SALL4 expression is generally assessed by immunohistochemistry (IHC) on whole section or tissue microarray (TMA), or by molecular testing, such as real time PCR and methylation analysis of the promoter region. Recently, several studies have proposed SALL4 as possible prognostic marker for cancer [6]. However, due to the lack of a comprehensive investigation of its prognostic value and to the different assessment techniques and protocols, the reliability of SALL4 as prognostic marker in cancer is still doubtful. In this work, we presented a systematic review and meta-analysis in order to investigate the prognostic role of SALL4 presence (SALL4+) in cancer patients by considering all-cause mortality and recurrence of cancer, evaluating whether SALL4+ can be associated with a poorer prognosis with respect to the absence of SALL4 (SALL4-).

RESULTS

Study and patient characteristics

The bibliographic search included 175 not-redundant articles. After excluding 137 articles based on title/abstract review, 38 articles were retrieved for full text review. 22 studies published after 2012 [729] were collected, reporting 23 independent cohorts (Supplementary Figure 1). The quality of the studies, assessed through NOS score [30], was generally good without any study with possible high risk of bias (Supplementary Table 1). Most of these 22 studies were carried out in Asia (n=18), three in USA, and one in Europe, mainly focusing on hepatocellular carcinoma (HCC) (n=11) and SALL4 was mainly assayed by IHC (n=18) rather than by molecular tests (n=6). In addition, 17 datasets of The Cancer Genome Atlas (TCGA) [31] were selected (see Methods section for additional details). All these studies were carried out in USA and SALL4 transcriptional activity was investigated by RNA sequencing technology. The selected TCGA datasets represent the following cancer types: Bladder Urothelial Carcinoma (BLCA), Breast invasive carcinoma (BRCA), Cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC), Colon and Rectum adenocarcinoma (COADREAD), Esophageal carcinoma (ESCA), Glioblastoma multiforme (GBM), Head and Neck squamous cell carcinoma (HNSC), Kidney renal clear cell carcinoma (KIRC), Liver hepatocellular carcinoma (LIHC), Lung adenocarcinoma (LUAD), Lung squamous cell carcinoma (LUSC), Ovarian serous cystadenocarcinoma (OV), Sarcoma (SARC), Skin Cutaneous Melanoma (SKCM), Stomach adenocarcinoma (STAD), Thyroid carcinoma (THCA) and Uterine Corpus Endometrial Carcinoma (UCEC).

40 different cohorts were considered for the meta-analysis, screening 9,947 patients (1,811 SALL4+ vs. 8,136 SALL4-) with a median follow-up period of 63.75 months (range: 7-250).

After dividing the participants by SALL4 status, no significant differences were found in terms of age (p=0.54), gender (p=0.9), tumor node metastasis (TNM) stage I-II percentage (p=0.21), grade (p=0.22), nodal (p=0.4) or vascular (p=0.31) invasion. Full descriptive details of all data are reported in Supplementary Tables 2a-2b.

Risk ratios and adjusted hazard ratios for overall and disease-free survival

The presence of SALL4 resulted significantly associated with higher death rate (55.2% SALL4+ vs. 39.8% SALL4- on average, RR=1.34; 95%CI: 1.21-1.48, p<0.0001; I2=66%) (Figure 1a). Similar findings were obtained also considering the recurrence of cancer as outcome (26 cohorts [12, 14, 16, 17, 23, 25, 28 and the TCGA cohorts]; recurrence rate: 51.5% SALL4+ vs. 42.2% SALL4-; RR=1.25; 95% CI: 1.1-1.42, p=0.0006; I2=62%) (Figure 1b). We then investigated the time-dependent risk of death/recurrence in terms of hazard ratios (adjusted for the maximum number of the covariates available in each study). 25 cohorts [14, 16, 18, 19, 23, 25, 26, 28 and the 17 TCGA cohorts] reported a survival analysis adjusted for a number of covariates ranging from 1 to 14 (Supplementary Table 3).

Figure 1A: Forest plot of the association between associating the presence/absence of SALL4 with and all-cause mortality.

Figure 1A: Forest plot of the association between associating the presence/absence of SALL4 with and all-cause mortality. Abbreviations: CI: Confidence Interval, SALL4: Spalt-Like Transcription Factor 4, RE: Random-Effects.

Figure 1B: Forest plot of the association between associating the presence/absence of SALL4 with recurrence of cancer.

Figure 1B: Forest plot of the association between associating the presence/absence of SALL4 with recurrence of cancer. Abbreviations: CI: Confidence Interval, SALL4: Spalt-Like Transcription Factor 4, RE: Random-Effects.

About all-causes mortality, the pooled adjusted HR showed a risk of death in SALL4+ patients significantly higher compared to SALL4- patients (HR=1.4; 95% CI: 1.19-1.65; p<0.0001; I2=63%) (Figure 2a). For the cancer recurrence, considering the available 4 studies [16, 23, 25, 28] and the 17 TCGA cohorts, a similar result was obtained confirming a worse prognostic role of SALL4+ (HR=1.52; 95% CI: 1.22-1.89, p=0.0002; I2=69%) (Figure 2b).

Figure 2A: Forest plot of the association associating the presence/absence of between SALL4 with and all-cause mortality using the hazard ratios adjusted for the maximum number of potential confounders available in each paper.

Figure 2A: Forest plot of the association associating the presence/absence of between SALL4 with and all-cause mortality using the hazard ratios adjusted for the maximum number of potential confounders available in each paper. Abbreviations: CI: Confidence Interval, SALL4: Spalt-Like Transcription Factor 4, RE: Random-Effects.

Figure 2B: Forest plot of the association associating the presence/absence of between SALL4 with recurrence of cancer using the hazard ratios adjusted for the maximum number of potential confounders available in each paper.

Figure 2B: Forest plot of the association associating the presence/absence of between SALL4 with recurrence of cancer using the hazard ratios adjusted for the maximum number of potential confounders available in each paper. Abbreviations: CI: Confidence Interval, SALL4: Spalt-Like Transcription Factor 4, RE: Random-Effects.

Meta-regression, sensitivity analyses and publication bias

The outcomes showed high heterogeneity (as I2≥50%, p<0.001), therefore the effects of possible moderators were also considered in the meta-analysis models in order to check whether they can explain to some extent the observed heterogeneity.

In Table 1, the sensitivity analysis for pre-planned moderators (i.e. continent, type of cancer, method of assessment of SALL4) seems to poorly explain this heterogeneity. It is worth to stress the paucity of the studies for some strata (e.g. for the hazard ratio of cancer recurrence), which does not allow a robust assessment for the effect of some moderators. Table 2 reports the meta-regression analyses for other potential moderators: difference of mean age in SALL4+ vs. SALL4- patients, and differences of their percentages in terms of female number, low stage/high grade tumors, node metastases, vascular invasions. Follow-up period and number of adjustments in the survival analysis were considered as well. The results showed that a possible moderator for both relative risk and adjusted hazard ratio in all-cause mortality can be tumor stage (p<0.05). However, the residual heterogeneity test demonstrated that there should be also other moderators, not available for the observed data, able to influence the overall survival (p<0.005).

Table 1: Stratification for some potential categorical moderator variables describing the association between SALL4 and all-cause mortality or recurrence of cancer

Moderator

Strata

Analysis details

RR for all-cause mortality

RR for recurrence of cancer

HR for all-cause mortality

HR for recurrence of cancer

Continent

Asia




Others

Pooled estimate, RR/HR (95%CI)
P-value for RR/HR
Heterogeneity, I2 (P-value)
Number of studies
Pooled estimate, RR/HR (95%CI)
P-value for RR/HR
Heterogeneity, I2 (P-value)
Number of studies
P-value *

1.42
(1.24-1.63)
<0.0001

70 (<0.0001)
20
1.26
(1.06-1.48)
0.0075

63 (<0.0001)
20
0.19

1.25
(0.98-1.58)
0.07

69 (0.004)
7
1.25
(1.07-1.47)
0.005

60 (0.0004)
19
0.96

1.71
(1.3-2.24)
<0.0001

41 (0.12)
7
1.28
(1.07-1.53)
0.007

57 (0.002)
18
0.08

2.46
(1.64-3.7)
<0.0001

38 (0.2)
3
1.37
(1.11-1.68)
0.003

57 (0.001)
18
0.02

Type of cancer

HCC



Others

Pooled estimate, RR/HR (95%CI)
P-value for RR/HR
Heterogeneity, I2 (P-value)
Number of studies
Pooled estimate, RR/HR (95%CI)
P-value for RR/HR
Heterogeneity, I2 (P-value)
Number of studies
P-value*

1.35
(1.14-1.61)
0.000569

(0.0004)
11
1.34
(1.17-1.53)
<0.0001

67 (<0.0001)
29
0.85

1.32
(1.03-1.69)
0.03

78 (0.0001)
7
1.2
(1.04-1.39)
0.01

44 (0.02)
19
0.51

1.61
(1.23-2.11)
0.0006

39 (0.15)
6
1.32
(1.1-1.58)
0.003

58 (0.0008)
19
0.21

2.09
(1.47-2.97)
<0.0001

12 (0.29)
2
1.45
(1.16-1.8)
0.001

63 (0.001)
19
0.24

SALL4 assessment

IHC



Others

Pooled estimate, RR/HR (95%CI)
P-value for RR/HR
Heterogeneity, I2 (P-value)
Number of studies
Pooled estimate, RR/HR (95%CI)
P-value for RR/HR
Heterogeneity, I2 (P-value)
Number of studies
P-value *

1.5
(1.26-1.78)
<0.0001

77 (<0.0001)
18
1.24
(1.09-1.4)
0.0008

52 (0.003)
22
0.12

1.5
(1.18-1.89)
0.0008

70 (0.001)
8
1.14
(0.99-1.32)
0.06

49 (0.01)
18
0.05

1.54
(1.24-1.92)
0.0001

50 (0.05)
8
1.29
(1.05-1.6)
0.02

59 (0.001)
17
0.19

1.99
(1.14-3.47)
0.02

86 (<0.0001)
4
1.4
(1.1-1.79)
0.006

56 (0.003)
17
0.24

* P-value represents the interaction between strata.

Abbreviations: CI: confidence interval, RR: Relative Risk, HCC: hepatocellular carcinoma, SALL4: Spalt-Like Transcription Factor 4, IHC: immunohistochemistry.

Table 2: Meta-regression analysis of potential continuous moderator variables the association between SALL4 and all-cause mortality or recurrence of cancer

RR for all-cause mortality

RR for Recurrence of cancer

HR for all-cause mortality

HR for Recurrence of cancer

Difference in mean age

-0.003 (-0.02 to 0.01)
(23 studies, p=0.73)

-0.01 (-0.04 to 0.02)
(19 studies, p=0.63)

-0.02 (-0.08 to 0.04)
(18 studies, p=0.49)

-0.003 (-0.06 to 0.06)
(18 studies, p=0.91)

Difference in women (%)

-0.001 (-0.01 to 0.01)
(36 studies, p=0.77)

-0.003 (-0.01 to 0.01)
(26 studies, p=0.5)

0.004 (-0.01 to 0.02)
(25 studies, p=0.59)

0.01 (-0.01 to 0.03)
(21 studies, p=0.33)

Difference in TNM 1-2 (%)

-0.01 (-0.02 to 0.001)
(30 studies, p=0.03)

-0.005 (-0.01 to 0.004)
(19 studies, p=0.3)

-0.02 (-0.03 to -0.003)
(20 studies, p=0.01)

-0.01 (-0.03 to 0.01)
(16 studies, p=0.21)

Difference in G3 (%)

-0.002 (-0.01 to 0.01)
(21 studies, p=0.68)

0.01 (-0.002 to 0.01)
(13 studies, p=0.16)

0.01 (-0.002 to -0.03)
(13 studies, p=0.1)

0.01 (-0.01 to 0.03)
(11 studies, p=0.35)

Difference in nodes metastasis (%)

0.001 (-0.02 to 0.02)
(7 studies, p=0.96)

Few studies

Few studies

Few studies

Difference in vascular invasion (%)

-0.002 (-0.01 to 0.01)
(10 studies, p=0.73)

0.001 (-0.03 to 0.03)
(5 studies, p=0.93)

-0.01 (-0.04 to 0.02)
(5 studies, p=0.65)

Few studies

Follow-up duration (months)

0.002 (-4e05 to 0.004)
(39 studies, p=0.05)

-0.0001 (-0.002 to 0.002)
(24 studies, p=0.93)

-0.003 (-0.0001 to 0.01)
(25 studies, p=0.06)

0.0002 (-0.004 to 0.01)
(21 studies, p=0.5)

Number of adjustments

Not applied

Not applied

0.02 (-0.04 to 0.07)
(25 studies, p=0.56)

0.03 (-0.03 to 0.1)
(21 studies, p=0.34)

All differences are calculated as values in SALL4+ vs. SALL4-, except follow-up duration and number of adjustments.

The data in cells are expressed as coefficients of the meta-regression model (with 95% confidence intervals).

Bold: significant results (p<0.05)

Abbreviations: SALL4: Spalt-Like Transcription Factor 4, RR: relative risk, HR: hazard ratio. TNM: Tumor Staging system, G: tumor grading.

Funnel plots inspection (Figure 3) and Egger’s test results (Table 3) showed that the publication bias was unlikely across all the outcomes, except for the relative risk in all-cause mortality. However, the trim-and-fill analysis suggested that there is a low risk of publication bias. Moreover, the high fail-safe number for each outcome confirmed the significance of our findings.

Funnel plots indicating that the publication bias was unlikely across all the outcomes, except for the relative risk in all-cause mortality.

Figure 3: Funnel plots indicating that the publication bias was unlikely across all the outcomes, except for the relative risk in all-cause mortality.

Table 3: Publication bias analysis

Outcome

Publication bias

Classic fail safe N

Egger bias (z-score)
p value

Kendall’s tau
>p-value

Trim and fill (95% CI)
[trimmed]

RR for all cause mortality

4.15
p < 0.0001

0.31
p=0.004

1.22 (1.09-1.37)
[8]

1278

RR for recurrence of cancer

0.52
p=0.6

0.02
p=0.9

1.25 (1.1-1.42)
[0]

279

HR for all cause mortality

1.98
p=0.05

0.23
p=0.1

1.2 (0.98-1.46)
[6]

361

HR for recurrence of cancer

0.56
p=0.57

0.07
0.7

1.52 (1.22-1.89)
[0]

285

Bold: significant results (p<0.05)

Abbreviations: RR: relative risk, HR: hazard ratio, CI: confidence interval.

DISCUSSION

SALL4 gene exerts its physiological role during embryo-fetal development, then gradually disappears and remains silenced in fully differentiated tissue, with the exception of germ cells and the hematopoietic stem/progenitor cells [4, 5]. Homozygous loss-of function mutation of SALL4 are lethal for the embryo, while heterozygous mutations cause renal-ocular syndrome and the Okihiro syndrome, associated with multi limbs defects, deficient eye movements, renal malformations, and deafness [32]. As for OCT4, evidence suggests that SALL4 is a major stemness factor and both SALL4 and OCT4 are expressed from the 2-cells stage embryo [33]. Moreover, SALL4 and OCT4 form a transcriptional auto-regulating core network, interacting with NANOG and SOX2 to form a multi-protein complex able to directly regulate both their own expression and the expression of hundreds of downstream target genes involved in pluripotency maintenance, such as Estrogen Related Receptor Beta (ESSRB), REST Corepressor 2 (RCOR2), Replication Timing Regulatory Factor 1 (RIF1) [34]. In addition, SALL4 complex can control gene accessibility recruiting the epigenetic repressor complex Mi-2/Nucleosome Remodelling and Deacetylase (NuRD) involved in Phosphatase and TENsin homologue (PTEN) and other SALL-family genes down-regulation [35]. The regulation of SALL4 expression is largely unknown. Putative upstream regulators of SALL4 include several promoter activators such as the Signal Transducer and Activator of Transcription 3 (STAT3), the Wnt/β-catenin pathway, and the Caudal type homeobox 1 (CDX1) [6]. Epigenetic regulation of SALL4 has also been proposed, but it needs further confirmation. Specifically, in induced pluripotent stem (iPS) cells, in embryonal stem cells (ESCs), and in SALL4+ cancers, the expression of SALL4 has been related to the hypo-methylation of its promoter [36]. Moreover, post-transcriptional regulation of SALL4 has been reported and, in particular, the inverse relationship between miR-107 and its expression on human glioma [37].

In cancers where aberrantly reactivated, SALL4 has been associated with the expression of many stemness-related genes (i.e. OCT4, NANOG, c-Myc, and SOX2) conferring to the cancer cells self-renewal pluripotency abilities [3]. Moreover, SALL4 expression allows cancer cells to acquire a stem-like phenotype, including: 1) increased mobility, invasion, and neoangiogenetic functions through the expression of epithelial-to-mesenchymal transition-related genes (e.g. SNAI1, CXCR4, TWIST1, CDH1, Vimentin, and ZEB1) [38]; 2) silencing of pro-apoptotic genes (e.g PTEN); 3) inducing the expression of chemo-resistance-related genes (e.g. ATP-Binding Cassette Multidrug Transporter) [39]; 4) acquisition of immune evasion abilities. All the above-mentioned features are related with the biological behavior of the SALL4+ cancers, which are more aggressive and associated with a worse prognosis compared with SALL4- ones. Recently, a meta-analysis showed a possible correlation between SALL4 expression and poor prognosis in cancer patients [40]. However, a large systematic evaluation about the role of SALL4 expression as prognostic marker including also the large amount of data available from public databases is still lacking. In this study, 22 perspective studies retrieved by the literature and 17 different TCGA cohorts were considered for the meta-analysis, involving 9,947 participants during a median follow-up period of about 5 years. Results showed that the expression of SALL4 was significantly associated with increased cancer mortality and recurrence, also after adjusting for potential confounders in the survival analyses. These findings suggest that the expression of this gene should be early assessed in patients. Therefore, the impact of SALL4 expression analysis should not be considered as a diagnostic and/or prognostic tool only, but it should be also investigated as possible target for new personalized treatments.

As a HCC prognostic factor, SALL4 seems to be particularly useful, as suggested by most of these 22 studies (11 HCC focused) and by the large TCGA dataset of HCC included in our analysis. However, these studies were mostly performed on patients with resectable tumors, thus with early stage of disease, for which the survival is longer as demonstrated by the mean follow-up period of 5 years; and so far a large investigation of SALL4 expression in advanced HCC is still lacking.

Interestingly, although SALL4 was discovered in 1994 [41, 42], all the studies included in this study were published only after 2012, highlighting that its role as prognostic marker and possible molecular target has been investigated only recently. SALL4 was firstly studied in germline cells from solid tumors and it resulted a valid diagnostic marker with a good sensitivity [43], but only later it was associated with poor prognosis in digestive system cancers [18]. A main point of our study was the long follow-up period (5 years on average), which was proper to evaluate the outcomes.

From a methodological point of view, most of the perspective studies used IHC, the most affordable technique to test indirectly genes expression, especially for large case studies. An advantage of IHC is the possibility to evaluate the expression of a marker in the whole tumor section allowing the detection of even a focal positivity in heterogeneous tumors (as it could be for SALL4, a stemness-related marker). However, two studies applied IHC on tissue microarrays instead of whole tumor sections and the low number of tissue cores obtained from each tumor sample could represent a limitation, increasing the risk of false negative cases [43]. In addition, some inconsistencies emerged by comparing the results from IHC studies, probably related to the different antibodies, IHC protocols and assessment criteria. Consensus standardization on SALL4 IHC would improve the reproducibility of the results. Real time PCR (RT-PCR) allows a precise quantification of SALL4 expression, although this method can be affected by tumor heterogeneity and it does not provide any information about the subcellular localization of SALL4. Looking at the TCGA datasets, SALL4 activity was investigated by RNA sequencing, which potentially allows the study of its expression in the entire transcriptome landscape, revealing possible regulative networks of SALL4 and related genes. However, as underlined for RT-PCR, gene expression could be affected by tumor heterogeneity, especially if the sample is not representative of the entire tumor mass. IHC, instead, can differently localize nuclear and cytoplasmic regions, allowing the detection of the position where a transcription factor is supposed to be active. Therefore, in our opinion, the best method to assess SALL4 presence in cancer is standardized IHC on whole tumor section.

The findings of our meta-analysis should be interpreted within their limitations and the most important of them is the inclusion of a limited number of studies. Secondly, through sensitivity and meta-regression analyses we were not able not explain the high heterogeneity found for the investigated outcomes. Finally, in the adjusted survival analyses, several important confounders (like cardiovascular diseases, disability, or other co-morbidities) were not available, further limiting the results. Despite these limitations, our meta-analysis suggests that SALL4 expression shortens overall survival as well as increases the rate of recurrences, even taking in consideration potential confounders. Since many epithelial cancers are characterized by SALL4 reactivation, this gene should be considered for developing future targeted therapeutic strategies. Due to its prognostic value, SALL4 expression should be considered as potential marker in the next-generation histopathological diagnosis, hopefully by standardized IHC protocols, integrating cancer morphological features and molecular targets information [4452].

MATERIALS AND METHODS

Literature search strategy

Two investigators independently conducted a literature search using PubMed and SCOPUS with no language restriction, from database inception to 1st September 2016, for perspective studies comparing relevant outcomes (i.e., all-cause mortality, cancer mortality and recurrence of disease/cancer) in patients with a diagnosis of cancer with loss vs. presence of expression of SALL4. In PubMed, the following search strategy was used: (“SALL4” OR “NM_020436”) AND (mortality OR mortalities OR fatality OR fatalities OR death* OR survival OR prognosis OR “hazard ratio” OR HR” OR “relative risk” OR “RR” OR “prognosis” OR “progression” OR “disease free survival” OR “DFS”). Conference abstracts and reference lists of included articles and those relevant to the topic were also hand-searched for identification of additional relevant articles. Any inconsistencies were resolved by consensus.

Study selection and quality assessment

Newcastle-Ottawa Scale (NOS) was used to evaluate study quality [30]. This systematic review was performed following the Meta-Analysis Of Observational Studies in Epidemiology (MOOSE) [53] guidelines and Preferred Reporting Items for Systematic reviews and Meta-Analyses (PRISMA) [54] statement. Inclusion criteria for this meta-analysis were: 1) perspective, longitudinal cohort study, 2) immuno-histochemical or molecular investigation of SALL4 expression, 3) diagnosis of cancer, 4) data about mortality or cancer recurrence. Exclusion criteria were: 1) no presence of cancer, 2) no data about relevant outcomes in the title/abstract, 3) no comparison between patients with SALL4+ vs. SALL4-, and 4) in vitro or animal research. To avoid overlaps between cohorts, in the presence of two or more studies from the same patient cohort, only the more recent study was taken into account and included in the meta-analysis.

Data extraction

Two investigators extracted key data from the included articles and a third independent investigator checked these data. For each study, information about authors, publication year, continent, histotype, SALL4 assessment methods, other genes analyzed, participant characteristics according by SALL4 expression data (e.g., age, gender, tumor stage and grading, percentage of participants with nodal metastasis and vascular invasion), number and type of adjustments in survival analysis, and the period of follow-up. When some information was missing, first and/or corresponding authors of the original article were contacted at least four times to obtain unpublished data.

In addition to these studies, one investigator retrieved molecular and clinical data of different cancer types from cBioPortal for Cancer Genomics (http://www.cbioportal.org) using the R package cgdsr version 1.2.5 (http://cran.r-project.org/web/packages/cgdsr/index.html). This package provides R functions for querying the Cancer Genomic Data Server (CGDS) hosted by the Computational Biology Center (cBio) at the Memorial Sloan-Kettering Cancer Center. For each dataset, SALL4 expression levels from RNAseq data were downloaded as z-scores, representing the number of standard deviations from the mean of expression using tumors diploid as the reference population. Datasets reporting a low number of subjects per class for the survival analysis (i.e. less than 10 individuals) were not considered.

Outcomes

The primary outcomes were number of deaths independent of the cause (all-cause mortality) and number of cancer recurrences during the follow-up period, in relation with the SALL4 presence or absence. The number of deaths due to cancer was preliminary considered as primary outcome, but no perspective study reported this information. As secondary outcomes, we considered hazard ratios (HRs) for all-cause mortality and recurrence, adjusted for the maximum number of confounders present for each study. Survival analysis was performed on TCGA data using the R packages survival version 2.39-5 (https://cran.r-project.org/web/packages/survival/index.html). SALL4 expression was classified in a binary way: one class of patients was identified by the presence of SALL4 (SALL4+) according to an over-expression quantified by a z-score>=2 and the other class was characterized by samples reporting a z-score<2 (SALL4-). For each TCGA dataset, HR for all-cause mortality and recurrence was estimated. Possible confounders of the final model were selected with a step-down procedure: the decision to remove confounders was based on a chi-square test for goodness of fit. All the available confounders were first included in the full model considering main effects only, and then they were sequentially removed if their removal did not result in a significant change of the estimates, using a threshold of 0.05 on the resulting p-values.

Data synthesis and statistical analysis

All statistical analyses were performed using R language version 3.3.0. R package metaphor version 1.9-9 (https://cran.r-project.org/web/packages/metafor/index.html) was used for the meta-analysis.

Descriptive characteristics of the patients, divided according to the presence or absence of SALL4, were compared using Wilcoxon rank-sum test. In particular, gender, tumor stage and tumor grading were represented in terms of percentages of females, patients with low-stage tumors (i.e. I-II) and patients with high-grade tumors (i.e. G3, G4), respectively. Vascular invasion and node metastasis were represented as percentages as well.

Pooled risk ratios (RRs) and adjusted hazard ratios (HRs) with 95% CIs were calculated for all-cause mortality and recurrence of cancer in patients with SALL4+ vs. SALL4- patients using DerSimonian-Laird random-effects models [55]. Heterogeneity across studies was assessed by the Cochrane I2 metric and chi square statistics. Given significant heterogeneity (p<0.05), meta-regression and sensitivity analyses were performed considering potential moderators and according to SALL4 status [56]. The following moderators were tested independently: continent (categorized as Asia vs. other continents), type of cancer (HCC or others), assessment methods (IHC or others), period of follow-up, number of adjustments applied in the survival model, and differences between SALL4+ and SALL4- participants in mean age, gender, tumor stage (divided in TNM stage 1-2, indicating low stage in percentage of SALL4+ vs. SALL4-patients), tumor grading (G3/G4, indicating high grade in percentage of SALL4+ vs. SALL4 patients), node metastasis and vascular invasions as percentage in SALL4+ vs. SALL4- patients.

Publication bias was assessed by visually inspecting funnel plots and using Egger’s bias test [56] and Begg-Mazumdar Kendall tau [57]. Then, to account for publication bias, we used the trim-and-fill method, based on the assumption that the effect sizes of all the studies are normally distributed around the center of a funnel plot; in the event of asymmetries, it adjusts for the potential effect of unpublished (imputed) studies [58].

For all the analyses, a p-value less than 0.05 was considered statistically significant.

CONFLICTS OF INTEREST

The authors declare that they have no conflicts of interest.

REFERENCES

1. Lim CY, Tam WL, Zhang J, Ang HS, Jia H, Lipovich L, Ng HH, Wei CL, Sung WK, Robson P, Yang H, Lim B. Sall4 regulates distinct transcription circuitries in different blastocyst-derived stem cell lineages. Cell Stem Cell. 2008; 3:543–54.

2. Zhang J, Tam WL, Tong GQ, Wu Q, Chan HY, Soh BS, Lou Y, Yang J, Ma Y, Chai L, Ng HH, Lufkin T, Robson P, Lim B. Sall4 modulates embryonic stem cell pluripotency and early embryonic development by the transcriptional regulation of Pou5f1. Nat Cell Biol. 2006; 8:1114–23.

3. Wang J, Rao S, Chu J, Shen X, Levasseur DN, Theunissen TW, Orkin SH. A protein interaction network for pluripotency of embryonic stem cells. Nature. 2006; 444:364–68.

4. Miettinen M, Wang Z, McCue PA, Sarlomo-Rikala M, Rys J, Biernat W, Lasota J, Lee YS. SALL4 expression in germ cell and non-germ cell tumors: a systematic immunohistochemical study of 3215 cases. Am J Surg Pathol. 2014; 38:410–20.

5. Gao C, Kong NR, Li A, Tatetu H, Ueno S, Yang Y, He J, Yang J, Ma Y, Kao GS, Tenen DG, Chai L. SALL4 is a key transcription regulator in normal human hematopoiesis. Transfusion. 2013; 53:1037–49.

6. Tatetsu H, Kong NR, Chong G, Amabile G, Tenen DG, Chai L. SALL4, the missing link between stem cells, development and cancer. Gene. 2016; 584:111–19.

7. Li A, Jiao Y, Yong KJ, Wang F, Gao C, Yan B, Srivastava S, Lim GS, Tang P, Yang H, Tenen DG, Chai L. SALL4 is a new target in endometrial cancer. Oncogene. 2015; 34:63–72.

8. Deng G, Zhu L, Huang F, Nie W, Huang W, Xu H, Zheng S, Yi Z, Wan T. SALL4 is a novel therapeutic target in intrahepatic cholangiocarcinoma. Oncotarget. 2015; 6:27416–26. doi: 10.18632/oncotarget.4862.

9. Han SX, Wang JL, Guo XJ, He CC, Ying X, Ma JL, Zhang YY, Zhao Q, Zhu Q. Serum SALL4 is a novel prognosis biomarker with tumor recurrence and poor survival of patients in hepatocellular carcinoma. J Immunol Res. 2014; 2014:262385.

10. Ma JC, Qian J, Lin J, Qian W, Yang J, Wang CZ, Chai HY, Li Y, Chen Q, Qian Z. Aberrant hypomethylation of SALL4 gene is associated with intermediate and poor karyotypes in acute myeloid leukemia. Clin Biochem. 2013; 46:304–07.

11. Liu L, Zhang J, Yang X, Fang C, Xu H, Xi X. SALL4 as an epithelial-mesenchymal transition and drug resistance inducer through the regulation of c-Myc in endometrial cancer. PLoS One. 2015; 10:e0138515.

12. Liu TC, Vachharajani N, Chapman WC, Brunt EM. SALL4 immunoreactivity predicts prognosis in Western hepatocellular carcinoma patients but is a rare event: a study of 236 cases. Am J Surg Pathol. 2014; 38:966–72.

13. Osada M, Aishima S, Hirahashi M, Takizawa N, Takahashi S, Nakamura K, Tanaka M, Maehara Y, Takayanagi R, Oda Y. Combination of hepatocellular markers is useful for prognostication in gastric hepatoid adenocarcinoma. Hum Pathol. 2014; 45:1243–50.

14. Park H, Lee H, Seo AN, Cho JY, Choi YR, Yoon YS, Han HS, Park YN, Kim H. SALL4 Expression in Hepatocellular Carcinomas Is Associated with EpCAM-Positivity and a Poor Prognosis. J Pathol Transl Med. 2015; 49:373–81.

15. Zeng SS, Yamashita T, Kondo M, Nio K, Hayashi T, Hara Y, Nomura Y, Yoshida M, Hayashi T, Oishi N, Ikeda H, Honda M, Kaneko S. The transcription factor SALL4 regulates stemness of EpCAM-positive hepatocellular carcinoma. J Hepatol. 2014; 60:127–34.

16. Shibahara J, Ando S, Hayashi A, Sakamoto Y, Hesegawa K, Kokudo N, Fukayama M. Clinicopathologic characteristics of SALL4-immunopositive hepatocellular carcinoma. Springerplus. 2014; 3:721.

17. Tanaka Y, Aishima S, Kohashi K, Okumura Y, Wang H, Hida T, Kotoh K, Shirabe K, Maehara Y, Takayanagi R, Oda Y. Spalt-like transcription factor 4 immunopositivity is associated with epithelial cell adhesion molecule expression in combined hepatocellular carcinoma and cholangiocarcinoma. Histopathology. 2016; 68:693–701.

18. Yong KJ, Gao C, Lim JS, Yan B, Yang H, Dimitrov T, Kawasaki A, Ong CW, Wong KF, Lee S, Ravikumar S, Srivastava S, Tian X, et al. Oncofetal gene SALL4 in aggressive hepatocellular carcinoma. N Engl J Med. 2013; 368:2266–76.

19. Zhang L, Yan Y, Jiang Y, Cui Y, Zou Y, Qian J, Luo C, Lu Y, Wu X. The expression of SALL4 in patients with gliomas: high level of SALL4 expression is correlated with poor outcome. J Neurooncol. 2015; 121:261–68.

20. Jung YK, Jang K, Paik SS, Kwon YJ, Kim HJ, Lee KG, Park HK, Choi D. Positive immunostaining of Sal-like protein 4 is associated with poor patient survival outcome in the large and undifferentiated Korean hepatocellular carcinoma. Ann Surg Treat Res. 2016; 91:23–28.

21. Kilic E, Tennstedt P, Högner A, Lebok P, Sauter G, Bokemeyer C, Izbicki JR, Wilczak W. The zinc-finger transcription factor SALL4 is frequently expressed in human cancers: association with clinical outcome in squamous cell carcinoma but not in adenocarcinoma of the esophagus. Virchows Arch. 2016; 468:483–92.

22. Oikawa T, Kamiya A, Zeniya M, Chikada H, Hyuck AD, Yamazaki Y, Wauthier E, Tajiri H, Miller LD, Wang XW, Reid LM, Nakauchi H. Sal-like protein 4 (SALL4), a stem cell biomarker in liver cancers. Hepatology. 2013; 57:1469–83.

23. Zhou S, Venkatramani R, Gomulia E, Shillingford N, Wang L. The diagnostic and prognostic value of SALL4 in hepatoblastoma. Histopathology. 2016; 69:822–30.

24. He J, Zhou M, Chen X, Yue D, Yang L, Qin G, Zhang Z, Gao Q, Wang D, Zhang C, Huang L, Wang L, Zhang B, et al. Inhibition of SALL4 reduces tumorigenicity involving epithelial-mesenchymal transition via Wnt/β-catenin pathway in esophageal squamous cell carcinoma. J Exp Clin Cancer Res. 2016; 35:98.

25. Yang M, Xie X, Ding Y. SALL4 is a marker of poor prognosis in serous ovarian carcinoma promoting invasion and metastasis. Oncol Rep. 2016; 35:1796–806.

26. Yin F, Han X, Yao SK, Wang XL, Yang HC. Importance of SALL4 in the development and prognosis of hepatocellular carcinoma. World J Gastroenterol. 2016; 22:2837–43.

27. Singhi AD, Seethala RR, Nason K, Foxwell TJ, Roche RL, McGrath KM, Levy RM, Luketich JD, Davison JM. Undifferentiated carcinoma of the esophagus: a clinicopathological study of 16 cases. Hum Pathol. 2015; 46:366–75.

28. He L, Ding H, Wang JH, Zhou Y, Li L, Yu YH, Huang L, Jia WH, Zeng M, Yun JP, Luo RZ, Zheng M. Overexpression of karyopherin 2 in human ovarian malignant germ cell tumor correlates with poor prognosis. PLoS One. 2012; 7:e42992.

29. Wang F, Guo Y, Chen Q, Yang Z, Ning N, Zhang Y, Xu Y, Xu X, Tong C, Chai L, Cui W. Stem cell factor SALL4, a potential prognostic marker for myelodysplastic syndromes. J Hematol Oncol. 2013; 6:73.

30. Wells G, Shea B, O’Connell D, Peterson J, Welch V, Losos M, Tugwell P. The Newcastle-Ottawa Scale (NOS) for assessing the quality if nonrandomized studies in meta-analyses. (Available from: URL: http://www.ohri.ca/programs/clinical_epidemiology/oxford.asp).

31. Weinstein JN, Collisson EA, Mills GB, Shaw KR, Ozenberger BA, Ellrott K, Shmulevich I, Sander C, Stuart JM, and Cancer Genome Atlas Research Network. The Cancer Genome Atlas Pan-Cancer analysis project. Nat Genet. 2013; 45:1113–20.

32. Al-Baradie R, Yamada K, St Hilaire C, Chan WM, Andrews C, McIntosh N, Nakano M, Martonyi EJ, Raymond WR, Okumura S, Okihiro MM, Engle EC. Duane radial ray syndrome (Okihiro syndrome) maps to 20q13 and results from mutations in SALL4, a new member of the SAL family. Am J Hum Genet. 2002; 71:1195–99.

33. Yang J, Gao C, Chai L, Ma Y. A novel SALL4/OCT4 transcriptional feedback network for pluripotency of embryonic stem cells. PLoS One. 2010; 5:e10766.

34. Loh YH, Wu Q, Chew JL, Vega VB, Zhang W, Chen X, Bourque G, George J, Leong B, Liu J, Wong KY, Sung KW, Lee CW, et al. The Oct4 and Nanog transcription network regulates pluripotency in mouse embryonic stem cells. Nat Genet. 2006; 38:431–40.

35. Lu J, Jeong HW, Kong N, Yang Y, Carroll J, Luo HR, Silberstein LE, Yupoma, Chai L. Stem cell factor SALL4 represses the transcriptions of PTEN and SALL1 through an epigenetic repressor complex. PLoS One. 2009; 4:e5577.

36. Amabile G, Di Ruscio A, Müller F, Welner RS, Yang H, Ebralidze AK, Zhang H, Levantini E, Qi L, Martinelli G, Brummelkamp T, Le Beau MM, Figueroa ME, et al. Dissecting the role of aberrant DNA methylation in human leukaemia. Nat Commun. 2015; 6:7091.

37. He J, Zhang W, Zhou Q, Zhao T, Song Y, Chai L, Li Y. Low-expression of microRNA-107 inhibits cell apoptosis in glioma by upregulation of SALL4. Int J Biochem Cell Biol. 2013; 45:1962–73.

38. Fassina A, Cappellesso R, Guzzardo V, Dalla Via L, Piccolo S, Ventura L, Fassan M. Epithelial-mesenchymal transition in malignant mesothelioma. Mod Pathol. 2012; 25:86–99.

39. Zhang X, Yuan X, Zhu W, Qian H, Xu W. SALL4: an emerging cancer biomarker and target. Cancer Lett. 2015; 357:55–62.

40. Cheng J, Gao J, Shuai X, Tao K. Oncogenic protein SALL4 and ZNF217 as prognostic indicators in solid cancers: a meta-analysis of individual studies. Oncotarget. 2016; 7:24314–25. doi: 10.18632/oncotarget.8237.

41. Kohlhase J, Schuh R, Dowe G, Kühnlein RP, Jäckle H, Schroeder B, Schulz-Schaeffer W, Kretzschmar HA, Köhler A, Müller U, Raab-Vetter M, Burkhardt E, Engel W, Stick R. Isolation, characterization, and organ-specific expression of two novel human zinc finger genes related to the Drosophila gene spalt. Genomics. 1996; 38:291–98.

42. Kühnlein RP, Frommer G, Friedrich M, Gonzalez-Gaitan M, Weber A, Wagner-Bernholz JF, Gehring WJ, Jäckle H, Schuh R. spalt encodes an evolutionarily conserved zinc finger protein of novel structure which provides homeotic gene function in the head and tail region of the Drosophila embryo. EMBO J. 1994; 13:168–79.

43. Mei K, Liu A, Allan RW, Wang P, Lane Z, Abel TW, Wei L, Cheng H, Guo S, Peng Y, Rakheja D, Wang M, Ma J, et al. Diagnostic utility of SALL4 in primary germ cell tumors of the central nervous system: a study of 77 cases. Mod Pathol. 2009; 22:1628–36.

44. Fassan M, Simbolo M, Bria E, Mafficini A, Pilotto S, Capelli P, Bencivenga M, Pecori S, Luchini C, Neves D, Turri G, Vicentini C, Montagna L, et al. High-throughput mutation profiling identifies novel molecular dysregulation in high-grade intraepithelial neoplasia and early gastric cancers. Gastric Cancer. 2014; 17:442–49.

45. Luchini C, Veronese N, Solmi M, Cho H, Kim JH, Chou A, Gill AJ, Faraj SF, Chaux A, Netto GJ, Nakayama K, Kyo S, Lee SY, et al. Prognostic role and implications of mutation status of tumor suppressor gene ARID1A in cancer: a systematic review and meta-analysis. Oncotarget. 2015; 6:39088–97. doi: 10.18632/oncotarget.5142.

46. Luchini C, Capelli P, Fassan M, Simbolo M, Mafficini A, Pedica F, Ruzzenente A, Guglielmi A, Corbo V, Scarpa A. Next-generation histopathologic diagnosis: a lesson from a hepatic carcinosarcoma. J Clin Oncol. 2014; 32:e63–66.

47. Luchini C, Nottegar A, Solmi M, Sergi G, Manzato E, Capelli P, Scarpa A, Veronese N. Prognostic implications of extranodal extension in node-positive squamous cell carcinoma of the vulva: A systematic review and meta-analysis. Surg Oncol. 2016; 25:60–65.

48. Luchini C, Veronese N, Pea A, Sergi G, Manzato E, Nottegar A, Solmi M, Capelli P, Scarpa A. Extranodal extension in N1-adenocarcinoma of the pancreas and papilla of Vater: a systematic review and meta-analysis of its prognostic significance. Eur J Gastroenterol Hepatol. 2016; 28:205–09.

49. Veronese N, Nottegar A, Pea A, Solmi M, Stubbs B, Capelli P, Sergi G, Manzato E, Fassan M, Wood LD, Scarpa A, Luchini C. Prognostic impact and implications of extracapsular lymph node involvement in colorectal cancer: a systematic review with meta-analysis. Ann Oncol. 2016; 27:42–48.

50. Veronese N, Luchini C, Nottegar A, Kaneko T, Sergi G, Manzato E, Solmi M, Scarpa A. Prognostic impact of extra-nodal extension in thyroid cancer: A meta-analysis. J Surg Oncol. 2015; 112:828–33.

51. Luchini C, Nottegar A, Pea A, Solmi M, Stubbs B, Capelli P, Sergi G, Manzato E, Fassan M, Wood LD, Scarpa A, Veronese N. Extranodal extension is an important prognostic parameter for both colonic and rectal cancer. Ann Oncol. 2016; 27:955–56.

52. Nottegar A, Veronese N, Senthil M, Roumen RM, Stubbs B, Choi AH, Verheuvel NC, Solmi M, Pea A, Capelli P, Fassan M, Sergi G, Manzato E, et al. Extra-nodal extension of sentinel lymph node metastasis is a marker of poor prognosis in breast cancer patients: A systematic review and an exploratory meta-analysis. Eur J Surg Oncol. 2016; 42:919–25.

53. Stroup DF, Berlin JA, Morton SC, Olkin I, Williamson GD, Rennie D, Moher D, Becker BJ, Sipe TA, Thacker SB. Meta-analysis of observational studies in epidemiology: a proposal for reporting. Meta-analysis Of Observational Studies in Epidemiology (MOOSE) group. JAMA. 2000; 283:2008–12.

54. Liberati A, Altman DG, Tetzlaff J, Mulrow C, Gøtzsche PC, Ioannidis JP, Clarke M, Devereaux PJ, Kleijnen J, Moher D. The PRISMA statement for reporting systematic reviews and meta-analyses of studies that evaluate health care interventions: explanation and elaboration. PLoS Med. 2009; 6:e1000100.

55. DerSimonian R, Laird N. Meta-analysis in clinical trials revisited. Contemp Clin Trials. 2015; 45:139–45.

56. Egger M, Davey Smith G, Schneider M, Minder C. Bias in meta-analysis detected by a simple, graphical test. BMJ. 1997; 315:629–34.

57. Begg CB, Mazumdar M. Operating characteristics of a rank correlation test for publication bias. Biometrics. 1994; 50:1088–101.

58. Duval S, Tweedie R. Trim and fill: A simple funnel-plot-based method of testing and adjusting for publication bias in meta-analysis. Biometrics. 2000; 56:455–63.


Creative Commons License All site content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 License.
PII: 14952