Skip to main content

Aberrant overexpression of transcription factor Forkhead box D1 predicts poor prognosis and promotes cancer progression in HNSCC



Forkhead box D1, the core transcription factor member of FOX family, has gradually seen as a key cancerous regulatory. However, its expression and carcinogenicity in head and neck squamous cell carcinoma (HNSCC) have not been reported yet. This study was to investigate its expression pattern, clinicopathological significance and biological roles in HNSCC.


HNSCC data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) was used to indicate the detailed expression pattern and outcome association of FOXD1, while Western Blot assay to detect FOXD1 level in a panel of HNSCC cell lines as well as immunocytochemistry to explore FOXD1 protein abundance and sublocation. Series of siRNA-mediated FOXD1 knock-down experiments to assess the proliferation, migration, invasion and anti- apoptosis ability after FOXD1 down-regulation. Bioinformatic analysis to find out which biological function and cancer-related pathways of FOXD1 associated genes involved in.


FOXD1 mRNA was significantly overexpressed in TCGA-HNSCC, GSE6631, GSE12452, GSE25099 and GSE30784. Besides, IHC results shown that nuclear location FOXD1 protein was significantly higher in primary HNSCC specimens from cohort involved in this study. Also, FOXD1 abundance was significantly correlated with cervical node metastasis and poor over-all/disease-free survival after combination analysis with patient pathological information. siRNA-mediated FOXD1 knock-down significantly inhibited cell proliferation, migration and invasion and induced apoptosis in HNSCC cells. Further analysis of GSEA, GO and KEGG showed that FOXD1 expression was significantly associated with oncological function and cancer-related pathways.


Taken together, our study implies that the potential oncogene, FOXD1, facilitates oncological behavior who can be identified as a brand-new HNSCC biomarker with diagnostic and prognostic significance.

Peer Review reports


Malignant tumor in head and neck mainly occurrences in oral cavity, nasal cavity, sinus, throat and pharynx, which pathological type is mainly squamous cell carcinoma (except thyroid tumor). It is the sixth most common malignant tumor in developed countries [1, 2]. According to the Global Cancer Statistics, there are over 550,000 new HNSCC patients and 300,000 HNSCC deaths worldwide every year [1, 3]. The main inducers of HNSCC include environmental risk, virus infection and gene mutation. Although clinical treatments of HNSCC patients has significantly been improved and improved the patient’s quality of life over the past decades, the prognosis of HNSCC is still unsatisfying, with the 5-year survival rate remaining around 50–60% [2, 3]. Therefore, exploring the potential molecular mechanisms of HNSCC will contribute to the development of molecular targeted therapies for this dreaded disease.

Forkhead box D1 (FOXD1), the transcription factor member of FOX family, mainly located on 5q13.2 (hg38.p13), encoding a vital DNA-binding protein containing 100 amino acids [4]. Numerous studies thus far have shown that FOXD1 as a key role in the occurrence, development and malignant regulation of multitudinous human malignance and its abnormal up-regulation was obviously related to cell proliferation, tumor metastasis and invasion as well as poor prognosis [5, 6]. However, FOXD1 function in the development of tumor biology and related molecular mechanism is still in preliminary exploration stage. Abnormal overexpression of FOXD1 gene may be as an oncogene, via facilitating cell proliferation and transformation, promoting the process of EMT and inhibiting cell apoptosis [7,8,9]. In addition, FOXD1 is further proved to regulate the cancer stem cells properties including tumor initiation, chemotherapy resistance and self-renewal capacities on mesenchymal glioma cells [10]. Moreover, FOXD1 can promote tumor progression and therapeutic resistance by inhibiting P27 mRNA expression in breast cancer [11].

All of the above results indicate that FOXD1 is probably a variety of new potential cancer-causing genes in human cancer. Although FOXD1 is thought to be a key oncogene in tumorigenesis, whether it takes an pivotal part in the carcinogenesis of HNSCC remains unknown. Therefore, in the present study, we collected relevant HNSCC materials from the clinical specimens, TCGA and GEO databases and executed the current systematic analysis to study the potential molecular regulatory mechanisms and clinical significance of FOXD1 in HNSCC. At the same time, the biological function of FOXD1 in HNSCC was further explored by combining small interference RNA and bioinformatics analysis techniques.

Materials and methods

Cell lines

A panel of HNSCC cell lines including Cal27, Fadu, SCC9, SCC25, HN4, HN6 and human immortalized oral mucosa keratinocytes (HOK) were used in this study. HOK, Cal27, Fadu, SCC4 and SCC25 cells were purchased from ATCC (USA). Specifically, two HNSCC cell lines HN4 and HN6 were gifted from Dr. Wantao Chen (SJTU, Shanghai, PRC). All cells lines were routinely tested as mycoplasma-free. Cells were grown in DMEM/F12 (Gibco, USA) supplemented with 10% fetal bovine serum (Gibco, USA) and penicillin streptomycin (1%, Beyotime, PRC), and culturaled under the condition of 37 °C (with 5% CO2). All regents were purchased from Sigma-Aldrich unless otherwise stated.

Small interference RNA (siRNA)

Two independent siRNA targeting h-FOXD1 mRNA region (NM_004472.3) were obtained from GenePharma (Shanghai, PRC). The detailed sequence as follows: siFOXD1–1:GAGCACUGAGAUGUCCGAUTT; siFOXD1–2:GGAAACAGACAUCGACGUGTT; siNC: TTCTCCGAACGTGTCACGT. 100 nM of siRNAs were transiently transfected into cells with Lipo-8000 (Beyotime, PRC), and harvested for further analysis at 48 h later.

Cell migration and invasion assay

In vitro experiments on cell migration were assessed by wound healing assay (experimental test time points: 0 h, 6 h, 12 h). Invasion experiments were carried on 24 aperture Transwell plank (8 μm, Conning, USA) with Corning Matrigel (BD, USA) (experimental test time points: 26 h). The detailed steps of two assays we all had previously reported [12].

Cell proliferation assay and apoptosis assay

Cell proliferation/viability was detected by CCK-8 cell viability assay (Beyotime, PRC). CCK-8 solution (10 μl per well) was added at specified time-points, incubating for 2 h and tested. Apoptosis cell were detected by FACs. Cells were treated with EDTA-free-trypsin (Gbico, USA), counted, incubated with Annexin V/PI Kit (BD, USA) and tested.

RNA extraction and qRT-PCR

Total RNA was extracted from cells with TRIzol reagent (Invitrogen, USA), reverse transcription to cDNA and using PrimeScriptTM RT-PCR kit (Takara, JP) for quantitative Real-time PCR. Relative mRNA expression was quantified as compared to internal control h-GAPDH using ΔCT method. The primers were listed as follows: h-FOXD1 forward, 5′-GATCTGTGAGTTCATCAGCGGC-3′ and reverse,5′-TGACGAAGCAGTCGTT GAGCGA-3′. h-N-cadherin forward,5′-CAAGATGGGTCAATGGAAATAG-3′ and reverse, 5′-CTCAGGAATACGAGCCTTCAC-3′; h-E-cadherin forward, 5′-AAGACAA AGAAGGCA AGGT-3’and reverse, 5′-AGAGAGTGTATGTGGCAATG-3′; h-Vimentin forward, 5′-CGAGGAGAGCAGGATTTCTC-3′ and reverse,5′-GGTATCAACCAGA GGGAGTGA-3′. h-SNAI1 (Snail) forward, 5′-CCTCCCTGTCAGATGAGGAC-3′ and reverse, 5′-CCAGGCTGAGGTATTCCTTG-3′. h-GAPDH forward, 5′-AGGTGAAGGTCGGA GTCAAC-3′ and reverse, 5′-AGTTGAGGTCAATGAAGGGG-3′.

Western blot analysis

The total protein of the cells was extracted by Cell lysis buffer for Western and IP (Beyotime, PRC) which contains with 1% PMSF (Beyotime, PRC). Protein samples were loaded and separated by 8–12% SDS-PAGE and transferred to 0.45 μm-PVDF membranes (Millipore, USA) and blocking by QuickBlock™ Blocking Buffer for Western Blot (Beyotime, PRC). These blots were incubated at 4 °C overnight with primary antibodies against FOXD1 (1:800 dilution, ab129324, Abcam, USA), E-cadherin (1:1000 dilution, #3195, CST), Vimentin(1:1000 dilution, #2118, CST), N-cadherin(1:1000 dilution, #13116, CST) and GAPDH (1:2000 dilution, #2118,CST) followed by incubations with the corresponding secondary antibodies.

Bioinformatics analysis of HNSCC

FOXD1 level and associated clinical data of HNSCC were downloaded from 2 publicly datasets: GEO ( and TCGA ( Specifically, R packages TCGAbiolinks was used for automatically sorting and normalizing the raw expression/clinical data of TCGA-HNSC. The overall and disease-free survival plot were downloaded from GEPIA2( FOXD1 correlated-gene chart (.tsv) was obtained from cBioPortal(, and genes with Spearman’s correlation > 0.2 was filtered for next step analysis. Median-cutoff was used to definition the patients from TCGA-HNSC cohort who has the high/low FOXD1 expression, then differential display analysis of gene expression between these two groups was performed by R packages limma. GSEA, GO and KEGG analysis was analysis and performed on R packages ClusterProfiler [13].

Clinical samples and immunohistochemical staining (IHC)

This study protocol was approved by the Research Ethic Committee of Nanjing Medical University. A retrospective cohort of a total number of 110 patients with primary HNSCC treated at our institution from Jan. 2011 to Dec. 2016 were enrolled. The set of standards for: no pre-treated (including tumor resection and neck dissection surgery, chemotherapy or radiotherapy) primary HNSCC patients, detailed information available including epidemiologic, clinical, pathological, follow-up data and written informed consent from these patients. We named this cohort as “NS20 cohort”. The archived tissue samples were retrieved and used for immunohistochemical staining. IHC score was qualified by senior pathologists dependently. 13 histologically confirmed normal oral mucosa samples from non-cancer surgeries were also enrolled. The detailed steps of IHC we had reported previously [14].

Statistical analyses

Data analysis was performed on GraphPad Prism 8.0.1 or R packages ggplot2. Normally distributed data were analyzed using an unpaired Student’s t test (two-tailed); multiple comparisons used nonparametric one-way ANOVA. Median cutoff and Log-rank (Mantel-Cox) test was used for survival Kapler-Merier plot. Statistical significance was defined as *P < 0.05 and **P < 0.01.


FOXD1 mRNA is obviously up-regulated in public HNSCC sets

Mounting study as well as pan-cancer analysis (Supplementary Figure1, downloaded from GEPIA2) have revealed that FOXD1 is overexpressed in a variety of cancers and has to do with unfavorable prognosis [11, 15]. Firstly, to investigate the FOXD1 mRNA expression in HNSCC, we mainly made use of GEO and TCGA databases to analysis the relevant information. Integration and analysis of the TCGA-HNSC cohort (502 cases) data showed that FOXD1 mRNA was obviously up-regulated in the TCGA-HNSC specimens compared with the normal counterpart (44 cases) using TCGA datasets (Fig. 1E). While as displayed in Figure1 A-D, data mining and questioning from GSE6631, GSE12452, GSE25099 and GSE30784 indicated that FOXD1 mRNA was significantly higher in HNSCC samples compared with the normal group, respectively. These data indicate that FOXD1 high expression may be related to the malignant regulation of HNSCC. The FOXD1 expression association among TCGA-HNSC subtypes was showed in Supplementary Figure2.

Fig. 1
figure 1

Overexpression of FOXD1 mRNA in five HNSCC datasets. A-E: The mRNA levels of FOXD1 (normalized and log2-transformed) were compared between HNSCC samples and normal counterparts in GSE6631 (A), GSE12452 (B), GSE25099(C), GSE30784(D)and TCGA-HNSC(E) datasets. The original data were retrieved from GEO database and TCGA, then plotted using Graphpad Prism 6.0 software. HNSCC, head neck squamous cell carcinoma; FOXD1, Forkhead box D1

Analysis of FOXD1 gene expression in HNSCC clinical samples and its clinical significance

Next, to further reveal the expression pattern of FOXD1 in HNSCC clinical specimens as well as its relationship with clinicopathological parameters, we semi-quantitatively detected FOXD1 protein expression pattern in 110 patients with primary HNSCC by immunohistochemical staining. All clinical data of these patients are shown in Table 1 (49 females/61 males, average 61.49 years). The patient’s latest follow-up time ranged from 3 to 83 months (mean 42.4 months). As displayed in Fig. 2, FOXD1 showed a positive staining mainly in nucleus in HNSCC paraffin sections, whereas weak/negative staining was identified in the normal oral mucosa and the stroma of HNSCC. Based on our IHC-scoring system, FOXD1 expression in HNSCC/normal mucosa was classified. Thus, the expression of FOXD1 protein can be classified as high (n = 61) or low expression group (n = 49) in HNSCC clinical samples and negative (n = 8), high (n = 1) or low (n = 4) expression in normal HNSCC clinical specimens. Therefore, these data confirmed that the FOXD1 protein was highly expressed in HNSCC (P < 0.01, Fisher’s exact test). The correlation between FOXD1 expression and clinicopathological parameters were displayed in Table 1 (Chi-square test or Fisher’s exact test). It is clear that there were no obviously association between FOXD1 abundance and patients age, drinking and heavy tobacco usage, tumor size and pathological grading. Notably, FOXD1 expression was positively associated with cervical lymph nodes metastasis and clinical stage with P value 0.0218 and 0.0105 (Fisher’s exact test), respectively.

Table 1 Expression of FOXD1 and its associations with clinicopathological parameters in 110 patients with HNSCC
Fig. 2
figure 2

Nuclear location FOXD1 enriched in human HNSCC samples. A, D: Representative negative staining of FOXD1 in normal oral epithelial (100× A, 200× D); B, E: Representative low expression of FOXD in primary human HNSCC sample (100× B, 200× E); C, F: Representative high expression of FOXD1 in primary human HNSCC sample (100× C, 200× F). Scale bar: 100 μm

Aberrant overexpression of FOXD1 in HNSCC patients associated with poor outcome

To explore prognostic significance of FOXD1 abundance in HNSCC patients, we analysis the relationship between its protein level and clinical prognostic. According to the latest follow-up data, 85 (77.28%) patients were still disease-freely alive, 8 (7.2%) survival with cervical nodal metastasis and/or local recurrences, whereas 25 (22.73%) patients died of post-surgical relapse, cancer metastases or other diseases. Further, as the Kaplan-Meier plots shown, patients with high FOXD1 abundance had obviously shorter overall-survival and disease-free survival than patients with low FOXD1 abundance (Log-rank, P = 0.0014, 0.0025, Figure3 A, B). Whereas, the similar conclusion form TCGA-HNSC cohort from GEPIA2 showed that the overall−/disease-free survival proportions in FOXD1high groups was also significantly lower than those in FOXD1low groups (Log-rank, P < 0.0001, 0.00018, Fig. 3C, D). To further evaluate the prognostic value of FOXD1 abundance in HNSCC, both univariate and multivariable survival analyses were executed by cox proportional hazards regression model. Univariate and multivariate cox analysis showed that high FOXD1 expression may be an independent prognostic factor for patients survival (as shown in Table 2.

Fig. 3
figure 3

High FOXD1 expression positively associates with reduced overall survival rates in HNSCC patients. A-A, B: Overall survival (A) and disease-free survival (B) analyses of HNSCC patients with high or low expression of FOXD1 were estimated by IHC. C, D: Overall survival (C) and disease-free survival (D) plot of FOXD1 in TCGA-HNSC patients downloaded from cBioPortal( Kaplan-Meier method with median cutoff and compared with Log-rank test

Table 2 Univariate and multivariate survival analyses (proportional hazards method) for patients with primary HNSCC

FOXD1 knockdown inhibited migration and invasion, proliferation and facilitated apoptosis in HNSCC cells

To explore the cancer-promoting effects in HNSCC, we first  examined the FOXD1 abundance in several HNSCC cell lines and carried out the loss-of-function assay by siRNA modulated knockdown. As displayed in Fig. 4A, the expression of FOXD1 protein in all HNSCC cell lines was obviously higher than that in immortalized non-tumorigenic cells (HOK). Cal27 and Fadu with relatively higher expression of endogenous FOXD1 protein were selected for further knockdown assays. Two independent siRNAs targeting human FOXD1(siFOXD1–1, siFOXD1–2) were brought in cal27 and Fadu cells to detect the changes of FOXD1 abundance and cell phenotype. As shown in Fig. 4B, the expression of FOXD1 protein decreased significantly after transfection with siRNA, which confirmed the effectiveness of our gene knockdown method. Cell proliferation was obviously repressed in two type of cells after FOXD1 expression inhibited as decision by results analysis from CCK-8 experiments (Fig. 4C). Moreover, results from Annexin V-PI Flow cytometry experiment showed that the apoptosis rate of siFOXD1-treated cells obviously elevated from 4.9 to 14.2% in Fadu and from 5.1 to 15.6% in Cal27, respectively (Fig. 4D, E). The migration and invasion ability of FOXD1 knockdown cells were also examined by wound healing and Transwell experiments, respectively. FOXD1 knockdown obviously decreased the ability of cell migration (Fig. 4F, G) and invasion (Fig. 4H-J) in vitro. Consistent with these observed phenotypical varies following FOXD1 knockdown, the protein and mRNA expression of EMT/metastasis-associated marker Vimentin、N-cadherin and Snail were low-expression concomitant with E-cadherin high-expression (Fig. 4K-M). These results suggested that FOXD1 is involved in the regulation of the malignant phenotype of HNSCC and may be a valuable potential therapeutic target.

Fig. 4
figure 4

FOXD1 knockdown inhibits cell migration and invasioncell proliferation, cell migration and invasion and triggers apoptosis in HNSCC cells. A: Endogenous FOXD1 protein expression was measured in a panel of HNSCC cell lines as compared to normal human oral keratinocytes (HOK). B: Endogenous FOXD1 was efficiently silenced by two siRNAs (siFOXD1–1, siFOXD1–2) in FADU and Cal27 cells. Non-targeting siRNA was utilized as negative control (siNC). C: Cell proliferation was remarkably suppressed when endogenous FOXD1 was silenced as measured by CCK8 assays. D, E: Increased percentages of cell undergoing apoptosis were evident following FOXD1 knockdown as assayed by Annexin V-PI staining. F, G: Cell migration was determined in cells with FOXD1 knockdown by wound-healing assay in Cal27 cells. H-G: Cell invasion was determined in cells with FOXD1 knockdown by transwell invasion assay in Cal27 and Fdau cells. K-M: The protein and mRNA abundance of migration/invasion-relevant marker E-cadherin, N-cadherin, snail and Vimentin was compared in cells infected siFOXD1 or control siNC. Scale bar: 100 μm. Data shown from three independent experiments, *P < .05, **P < .01, Kruskall-Wallis test

Functional enrichment analysis of the differentially expressed genes between FOXD1 low/high groups

In order to complement the in vitro loss-of-function experiment in investigating FOXD1 pro-tumorigenic functions. Firstly, we downloaded the FOXD1 correlated-gene chart from TCGA-HNSC on GEPIA2, filtered the correlation score > 0.2 geneset (Fig. 5A) and performing gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, and the data suggests that FOXD1-positive-corrlated genes enriched in numerous cancer related biofunction and pathways such as cell-substrate junction and focal adhesion(Fig. 5B, C). Then, we further used TCGA-HNSC data for bioinformatics analysis to recognize the nominee genes that might be potentially associated with FOXD1, the expression of which was divided into high-expression group and low-expression group (median-cutoff) by the method of bipartition in TCGA-HNSC data and then were in thrall to GO and KEGG and Gene Set Enrichment Analysis (GSEA) analyses (Fig. 5D-F). GO and KEGG analysis the up-regulated genes in high-risk group associated with FOXD1 were significantly enriched in Focal adhesion, Cytokine-cytokine receptor interaction, IL-17 signaling pathway, ECM-receptor interaction, TNF signaling pathway, TGF-beta signaling (Fig. 6A). On the contrary, KEGG and GO analysis for the down-regulated genes associated with FOXD1 in Low-risk group were related to apoptotic process, metabolism (Fig. 6B). GSEA analysis revealed that differentially expressed genes between high-FOXD1 and Low-FOXD1 groups were involved significant pathways included KEGG Pathways in cancer (Fig. 6C), EMT pathway (Fig. 6D), KEGG Focal adhesion (Fig. 6E) and Apoptosis (Fig. 6F). In conclusion, combined with cell assay in vitro and bioinformatics results extremely support the idea that FOXD1 is an newly oncogene in HNSCC.

Fig. 5
figure 5

The biofunction of FOXD1 correlation genes. A: Schematic diagram of the numbers of FOXD1-associated-genes with Spearman’s correlation > 0.2 subset. B, C: The top 15 involved significant biological process (GO, B) and pathways (KEGG, C) in these positive correlation genes. D-E: Scatter plot (D) displaying the FOXD1 level of each patient from TCGA-HNSC dataset. According to the median-cutoff FOXD1 expression, patients were divided into high- or low-expression groups and survival status plot (E) shows the patients with the high- or low- group. F: Schematic diagram of further bioinformatic analysis between high- or low- group

Fig. 6
figure 6

The GO, KEGG and GSEA analysis results for the differentially expressed genes between high-FOXD1 and Low-FOXD1 groups. A: GO (upper panel) and KEGG (low panel) analysis for the genes in high-risk group. B: GO (upper panel) and KEGG (low panel) analysis for the genes in low-risk group.C-F: The most involved significant pathways via GSEA


As a core member of the human FOX gene family, FOXD1 expression may be an oncogene involving cancer progression and has been critically participated in cell migration and apoptosis, cancer stem cell self-renewal as well as other physiological and pathological processes [16,17,18]. Previous studies have shown that FOXD1 may be an oncogene related to tumor development and act as a potential diagnostic biomarker and anti-cancer therapeutic target [9, 11, 19]. Herein, this study revealed the FOXD1 expression pattern in HNSCC, uncovered its prognostic and clinicopathological significance and also revealed its the oncogenic roles of FOXD1 by the loss-of function method and bio-informatics analysis. Our results together highly suggested that as a newly hypothesized oncogene FOXD1 promotes HNSCC development and also a novel biomarker with clinical translation potential.

HNSCC initiation and development are attributed through consecutive histopathological stages from normal epithelial cells to squamous-cell carcinoma cells which is mainly due to the activation of genetic predisposition, oncogenes and the inactivation of tumor suppressor genes [20,21,22]. In particular, dysregulation of the FOX family have been proved to be involved in nearly all stages of HNSCC tumorigenesis [23]. Among these FOX gene regulators, FOXD1 has been increasingly recognized as a key agent for cancer initiation, unrestricted growth, and metastasis, and has been proven to be a promising therapeutic target [17].

Here, we found that FOXD1 is significantly elevated in most HNSCC specimens as proofed by obviously overexpressed FOXD1 mRNA in GEO and TCGA datasets as well as up-regulation of FOXD1 protein in our HNSCC patients cohort. To our knowledge, this is the first study to confirm an abnormal pattern of FOXD1 overexpression in HNSCC. Previous findings have shown the clinical significance of FOXD1 abundance in various human cancers. For instance, elevated FOXD1 expression is obviously associated with lymph nodes metastasis, tumor size and advanced stages in non-small cell lung cancer [15], nasopharyngeal carcinoma and colorectal cancer [23, 24]. We further established overexpression pattern of FOXD1 and shown that its overexpression significantly associated with lymph node metastasis and clinical stage in HNSCC, while there was no significant correlation between FOXD1 expression and other clinicopathologic parameters. Moreover, up-regulation of FOXD1 was obviously associated with the decrease of survival rate via Kaplan-Meier survival data analysis, thus suggesting that FOXD1 might be a valuable prognostic biomarker for HNSCC. Therefore, Multivariate survival analysis using cox proportional hazards regression model results confirmed that FOXD1 expression was an independent prognostic factor for HNSCC. This was consistent with our previous research results in OSCC [16] and other associated OSCC research [22, 25]. Accumulating evidence has revealed that FOXD1 is vitally participated in tumorigenesis by facilitating cell migration, invasion, proliferation, angiogenesis and inhibiting apoptosis [9, 26]. Consistent with this, our in vitro loss-of-function experiment results indicate that FOXD1 knockdown leads to a diminished in proliferation, migration and invasion in HNSCC cells. Additionally, bioinformatics analysis results further indicated that FOXD1 was involved significant pathways included KEGG Pathways in cancer, EMT, Focal adhesion, Apoptosis, TGF-beta signaling, Cell cycle etc. Consistently, previous studies have revealed that FOXD1 facilitates cell proliferation and inhibits apoptosis by promoting the expression of PIK2 in colorectal cancer [24]. Fengping Pan and his colleagues revealed that FOXD1 promoted CRC cell malignant phenotypes by activation of the ERK 1/2 signaling pathway [6]. Complementary, in melanoma cancer, FOXD1 was critically involved in invasion and migration via indirect regulation of RAC1b and MMP-9 alternative splicing in melanoma cells [8]. Of note, the inhibition of FOXD1 expression affected cell related malignant phenotypes such as cell proliferation, cell migration and invasiveness, cell apoptosis and ultimately leads to tumor regression and reduction of tumor metastasis. Taken together, combined with bio-informatics analysis and the results of our findings and others strongly favors the notion that FOXD1 probably functions as a hypothetical oncogene by facilitating cancer cell proliferation, migration and invasion. Suppression of FOXD1 expression via genetic method might hold translational promise in the treatment of HNSCC.


In conclusion, our all results indicate that a subgroup of HNSCC patients who FOXD1 aberrantly up-regulation in an important and reveals its pro-oncogenic role in promoting the occurrence and development of HNSCC. However, there is still a large gap between the role of FOXD1 in the development of HNSCC and the effective treatment of FOXD1.

Availability of data and materials

The dataset used during the present study is available from the corresponding author upon a reasonable request.


  1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.

    Article  Google Scholar 

  2. Johnson DE, Burtness B, Leemans CR, Lui VWY, Bauman JE, Grandis JR. Head and neck squamous cell carcinoma. Nat Rev Dis Primers. 2020;6(1):92.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2019. Ca-Cancer J Clin. 2019;69(1):7–34.

    Article  PubMed  Google Scholar 

  4. Katoh M, Katoh M. Human FOX gene family (review). Int J Oncol. 2004;25(5):1495–500.

    Article  PubMed  CAS  Google Scholar 

  5. Gao YF, Zhu T, Mao XY, Mao CX, Li L, Yin JY, et al. Silencing of Forkhead box D1 inhibits proliferation and migration in glioma cells. Oncol Rep. 2017;37(2):1196–202.

    Article  PubMed  CAS  Google Scholar 

  6. Pan F, Li M, Chen W. FOXD1 predicts prognosis of colorectal cancer patients and promotes colorectal cancer progression via the ERK 1/2 pathway. Am J Transl Res. 2018;10(5):1522–30.

    PubMed  PubMed Central  CAS  Google Scholar 

  7. Cheng P, Wang J, Waghmare I, Sartini S, Coviello V, Zhang Z, et al. FOXD1-ALDH1A3 signaling is a determinant for the self-renewal and Tumorigenicity of mesenchymal glioma stem cells. Cancer Res. 2016;76(24):7219–30.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. Wu H, Larribere L, Sun Q, Novak D, Sachindra S, Granados K, et al. Loss of neural crest-associated gene FOXD1 impairs melanoma invasion and migration via RAC1B downregulation. Int J Cancer. 2018;143(11):2962–72.

    Article  PubMed  CAS  Google Scholar 

  9. Chen C, Xu ZQ, Zong YP, Ou BC, Shen XH, Feng H, et al. CXCL5 induces tumor angiogenesis via enhancing the expression of FOXD1 mediated by the AKT/NF-kappaB pathway in colorectal cancer. Cell Death Dis. 2019;10(3):178.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  10. Koga M, Matsuda M, Kawamura T, Sogo T, Shigeno A, Nishida E, et al. Foxd1 is a mediator and indicator of the cell reprogramming process. Nat Commun. 2014;5(1):3197.

    Article  PubMed  CAS  Google Scholar 

  11. Zhao YF, Zhao JY, Yue H, Hu KS, Shen H, Guo ZG, et al. FOXD1 promotes breast cancer proliferation and chemotherapeutic drug resistance by targeting p27. Biochem Biophys Res Commun. 2015;456(1):232–7.

    Article  PubMed  CAS  Google Scholar 

  12. Li Z, Wang Y, Zhu Y, Yuan C, Wang D, Zhang W, et al. The hippo transducer TAZ promotes epithelial to mesenchymal transition and cancer stem cell maintenance in oral cancer. Mol Oncol. 2015;9(6):1091–105.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  14. Li J, Li Z, Wu Y, Wang Y, Wang D, Zhang W, et al. The hippo effector TAZ promotes cancer stemness by transcriptional activation of SOX2 in head neck squamous cell carcinoma. Cell Death Dis. 2019;10(8):603.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Nakayama S, Soejima K, Yasuda H, Yoda S, Satomi R, Ikemura S, et al. FOXD1 expression is associated with poor prognosis in non-small cell lung Cancer. Anticancer Res. 2015;35(1):261–8.

    PubMed  CAS  Google Scholar 

  16. Li Z, Yan T, Wu X, Zhang W, Li J, Wang L, et al. Increased expression of FOXD1 is associated with cervical node metastasis and unfavorable prognosis in oral squamous cell carcinoma. J Oral Pathol Med. 2020;49(10):1030–6.

    Article  PubMed  CAS  Google Scholar 

  17. Jin Y, Liang Z, Lou H. The Emerging Roles of Fox Family Transcription Factors in Chromosome Replication, Organization, and Genome Stability. Cells. 2020;9(1).

  18. Gong Z, Yu J, Yang S, Lai PBS, Chen GG. FOX transcription factor family in hepatocellular carcinoma. Biochim Biophys Acta Rev Cancer. 1874;2020(1):188376.

    Article  CAS  Google Scholar 

  19. Li CH, Chang YC, Hsiao M, Liang SM. FOXD1 and Gal-3 Form a Positive Regulatory Loop to Regulate Lung Cancer Aggressiveness. Cancers (Basel). 2019;11(12).

  20. Leemans CR, Snijders PJF, Brakenhoff RH. The molecular landscape of head and neck cancer. Nat Rev Cancer. 2018;18(5):269–82.

    Article  PubMed  CAS  Google Scholar 

  21. Cramer JD, Burtness B, Le QT, Ferris RL. The changing therapeutic landscape of head and neck cancer. Nat Rev Clin Oncol. 2019;16(11):669–83.

    Article  PubMed  Google Scholar 

  22. Zhang D, Zhang Y, Sun X. LINC01133 promotes the progression of cervical cancer via regulating miR-30a-5p/FOXD1. Asia Pac J Clin Oncol. 2020;17(3):253–63.

    Article  PubMed  Google Scholar 

  23. Zhang Y, Zhang W. FOXD1, negatively regulated by miR-186, promotes the proliferation, metastasis and radioresistance of nasopharyngeal carcinoma cells. Cancer Biomark. 2020;28(4):511–21.

    Article  PubMed  CAS  Google Scholar 

  24. Han T, Lin J, Wang Y, Fan Q, Sun H, Tao Y, et al. Forkhead box D1 promotes proliferation and suppresses apoptosis via regulating polo-like kinase 2 in colorectal cancer. Biomed Pharmacother. 2018;103:1369–75.

    Article  PubMed  CAS  Google Scholar 

  25. Lin CH, Lee HH, Chang WM, Lee FP, Chen LC, Lu LS, et al. FOXD1 Repression Potentiates Radiation Effectiveness by Downregulating G3BP2 Expression and Promoting the Activation of TXNIP-Related Pathways in Oral Cancer. Cancers (Basel). 2020;12(9).

  26. Tao J, Cong H, Wang H, Zhang D, Liu C, Chu H, et al. MiR-30a-5p inhibits osteosarcoma cell proliferation and migration by targeting FOXD1. Biochem Biophys Res Commun. 2018;503(2):1092–7.

    Article  PubMed  CAS  Google Scholar 

Download references


We would like to thank the National Natural Science Foundation of China (81602378), The Natural Scientific Foundation of the Jiangsu Higher Education Institutions of China (16KJB320002). whose financial support has made this study possible.

Code availability

Not applicable.

Consent to participate

Not applicable.


This work is financially supported, in whole or in part, This work is financially supported, in whole or in part, by National Natural Science Foundation of China (81602378), Jiangsu Higher Education Institutions of China (16KJB320002).

Author information

Authors and Affiliations



JL and TY performed the experimental study, data collection and analysis and manuscript writing. XW, XL and XK carried out the most experiments. YZ and JY performed histological and statistical analyses. ZL conceived and supervised the whole project. All authors have made an intellectual contribution to the manuscript and approved the submission. JL and TY are co-frst authors.

Corresponding author

Correspondence to Zhongwu Li.

Ethics declarations

Ethics approval and consent to participate

We confirm that we have passed the ethics committee of Nanjing Medical University and all methods were carried out in accordance with the guidelines and regulations. All biological samples were obtained with patients’ written informed consent in the study.

Consent for publication

Not applicable.

Competing interests

All authors declare no conflict of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1 Supplementary Figure1. FOXD1 mRNA expression pattern in pan-cancer analysis


Additional file 2 Supplementary Figure2. FOXD1 mRNA expression pattern in TCGA-HNSCC subtypes. A-F:

The mRNA levels of FOXD1 (normalized and log2-transformed) from TCGA-HNSCC datasets were compared between clinical subtypes: Clinical stage(I- IV, A), Pathological grade(I-III, B), HPV infection status(C), Tumor size(T1-T4, D), Cervical nodal metastasis(E), Distant metastasis(F). Student’s t test or Kruskall-Wallis test.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, J., Yan, T., Wu, X. et al. Aberrant overexpression of transcription factor Forkhead box D1 predicts poor prognosis and promotes cancer progression in HNSCC. BMC Cancer 21, 1205 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: