- Open Access
SOX9 drives the epithelial–mesenchymal transition in non-small-cell lung cancer through the Wnt/β-catenin pathway
Journal of Translational Medicine volume 17, Article number: 143 (2019)
The distant metastasis of cancer cells is a risk factor for tumor lethality and poor prognosis in non-small-cell lung carcinoma (NSCLC). Increased SOX9 expression has been associated with clinical stage and poor prognosis in NSCLC, but the molecular mechanisms by which SOX9 promotes metastasis in NSCLC are still unknown.
The relationship between SOX9 expression and T, N, M classification was assessed using the χ2 test and Spearman’s analysis in 142 immunohistochemically diagnosed specimens of NSCLC. We also generated SOX9-overexpression and SOX9-knockdown cells lines and their corresponding control cell lines by transfection with lentiviral constructs. In vivo assay, SOX9-overexpressing and SOX9-knockdown NSCLC cells were injected in zebrafish to examine distance metastasis. Gene set enrichment analysis (GSEA) was applied to analysis the correlation between SOX9 overexpression and Wnt/β-catenin pathway. Luciferase assay was used to check transcriptional activity of TCF/LEF and western blot and immunofluorescence was employed to detect β-catenin translocation in SOX9-overexpression, SOX9-knockdown and their corresponding control cell lines.
We found that SOX9 overexpression correlates with the T, N and M stage significantly (p = 0.03, 0.000, and 0.032 respectively) in 142 immunohistochemically diagnosed specimens of NSCLC. SOX9 overexpression was found to decrease the expression of the epithelial cell markers E-cadherin and γ-catenin and increase the expression of the mesenchymal cell markers N-cadherin and vimentin. An in vivo assay showed distant metastasis of the SOX9-overexpressing cells, which was not observed in the SOX9-knockdown cells. These findings indicate that SOX9 promotes distant metastasis by promoting EMT in NSCLC cells. GSEA showed that SOX9 overexpression was significantly correlated with the Wnt/β-catenin pathway which was corroborated by the expression of EMT-associated proteins in this pathway and its downstream target genes. SOX9 overexpression was also found to enhance the transcriptional activity of TCF/LEF, promote the nuclear translocation of β-catenin and increase the phosphorylation of GSK3β at Ser9. Further, inhibition of β-catenin suppressed the metastasis-promoting effects of SOX9 overexpression.
This study is the first to report that SOX9 is associated with clinical TNM stage and indicates that SOX9 promotes migration, invasion and the EMT process through the Wnt/β-catenin pathway.
According to the Cancer Registration Report released by the National Central Cancer Registry of China in 2015, there are ~ 3.1 million newly diagnosed cancer patients and ~ 2 million cancer-associated deaths in China each year . Lung cancer has the highest incidence and mortality rate among all tumors, and non-small cell lung cancer (NSCLC) is the most common lung cancer type, accounting for 80% of lung cancer-associated deaths . Local cancer recurrence and distant metastases are the major causes of 90% deaths, and they are also the main prognostic factors in patients with NSCLC. Unfortunately, the majority of NSCLC patients are diagnosed with terminal-stage cancer upon hospital admission and are insensitive to chemotherapy. Due to the heterogeneous nature of the cancer, the causes of distant metastasis remain unclear. From the viewpoint of early diagnosis and timely treatment of NSCLC, it is important to understand the causes and mechanisms of distant metastasis, as well as the molecules involved.
The transcription factor SOX9 (SRY-related high mobility group-box 9) belongs to the SOX protein family, which includes SOX8, SOX9, SOX10 and SOX E . SOX9 is evolutionarily conserved across many species and has essential functions in mediating embryonic development and sex determination . Many basic and clinical studies have associated SOX9 overexpression with tumor malignancy and tumorigenesis in lung, colorectal, prostate and breast cancer [5,6,7,8,9]. For example, SOX9 knockout was found to decrease lung cancer tumorigenesis, whereas exogenous SOX9 overexpression enhanced tumor volume in mice . Further, SOX9 upregulation enhances breast cancer cell proliferation in highly metastatic breast cancers  and is an independent marker of poor prognosis and a low 5-year overall survival rate [11, 12]. Our previous research found that the overexpression of SOX9 is associated with clinical stage and poor survival in NSCLC patients . With regard to the underlying mechanism, it was found that activation of cGK-II (a cGMP-dependent protein kinase) led to the phosphorylation and inhibition of the SOX9 protein, which in turn suppressed AKT phosphorylation and arrested cultured glioma cells in the G1 phase; siRNA-induced suppression of SOX9 was also found to reduce glioma cell proliferation . These data indicate that SOX9 is a potential biomarker for early stage diagnosis and may be a key point of clinical intervention in NSCLC. However, the biological functions and mechanism of action of SOX9 in promoting invasion and metastasis in NSCLC remain unclear.
The overexpression of SOX9 has also been detected in the lymph nodes of prostate cancer xenografts, in which SOX9 was found to enhance cellular growth, angiogenesis and invasion by promoting epithelial cells to expand into mesenchyme cells during fetal prostate development . Epithelial–mesenchymal transition (EMT) is known as a key indicator of the early stage in the NSCLC metastasis process. Thus, the role of SOX9 in promoting EMT might be a good starting point to investigate the mechanisms via which it may influence metastasis in NSCLC, and this was the aim of the present study.
Here, we report that SOX9 overexpression decreased the protein level of epithelial cell markers and increased the protein level of mesenchymal cell markers in NSCLC cell lines, thus confirming the influence of SOX9 on EMT. The findings were corroborated by in vivo experiments in zebrafish models, in which SOX9 overexpression was found to promote distant metastasis. Further investigation into the underlying mechanism revealed that the Wnt/β-catenin pathway was involved in the metastasis-promoting effects of SOX9. In our previous study, we found that SOX9 was a predictor of poor prognosis and over expression of SOX9 had a significant correlation with the early histological stage subgroup (stages I and II) and the late histological stage subgroup (stages III and IV) in 142 NSCLC tissue sections . To our knowledge, this study is the first to report that SOX9 is associated with clinical TNM stage in NSCLC, and to describe the mechanism by which SOX9 enhances Wnt/β-catenin pathway via promoting GSK3β phosphorylation at serine 9. Moreover, our findings support a significant role of SOX9 in facilitating migration, invasion, and the EMT process through the Wnt/β-catenin pathway and provide a molecular understanding of the involvement of SOX9 in NSCLC metastasis which indicates that SOX9 may be an important biomarker in NSCLC progression.
Patients and tissue specimens
This study was conducted on a total of 142 paraffin-embedded lung cancer specimens, which were obtained from patients who were diagnosed histopathologically at the Second Affiliated Hospital of Guangzhou Medical College between 2006 and 2010 . Prior patient consent and approval from the Institutional Research Ethics Committee were obtained for the use of the clinical materials for research. The tumors were staged according to the 7th edition of the Cancer Stage Manual of the American Joint Committee on Cancer .
Cell lines, culture conditions and reagents
Lung cancer cell lines A549 and NCI-H460 were purchased from the American Type Culture Collection (Manassas, VA, USA) and cultured in the Dulbecco’s modified Eagle medium (DMEM; Invitrogen, Carlsbad, USA) with 10% fetal bovine serum (FBS, Invitrogen) and 1% penicillin/streptomycin (Gibco, USA) at 37 °C in a 5% CO2 atmosphere, and 0.25% trypsin (Gibco, USA) was used to passage the cell lines till they reached 90% confluence. The β-catenin inhibitor XAV-939 was purchased from Selleck (S1180; Houston, TX).
Plasmids, retroviral infection, and transfection
Human SOX9 cDNA was amplified by PCR and cloned into the pSin-EF2 lentiviral vector. Retroviral production and infection were performed as described previously . Stable cell lines that expressed SOX9 were selected for 10 days with 0.5 mg/ml puromycin. We constructed stable SOX9-overexpressing cell lines A549-SOX9 and NCI-H460-SOX9, and their corresponding stable control cell lines A549-pSin-Vector and NCI-H460-pSin-Vector respectively. shRNA-mediated interference was used to create SOX9-knockdown A549 and NCI-H460 cells by transfection with lentiviral constructs expressing shRNAs against SOX9 (Ruibo, China) or control shRNA for 24 h. Cells that were positive for SOX9 knockdown were selected with 0.5 mg/ml puromycin for 10 days. We then constructed stable SOX9-knockdown cell lines A549-SOX9 sh1# and NCI-H460-SOX9 sh1# and their corresponding stable control cell lines A549-pSuper-Vector and NCI-H460-pSuper-Vector respectively.
RNA extraction and real-time reverse transcription-polymerase chain reaction
Total RNA was extracted from cultured cells with Trizol reagent (Invitrogen), and the RNA was purified using a PureLink RNA mini kit (Invitrogen) according to the manufacturer’s instructions. Real-time reverse transcription-polymerase chain reaction (RT-PCR) was used to quantify Wnt/β-catenin pathway downstream target gene mRNA levels in A549-SOX9, NCI-H460-SOX9, A549-SOX9 sh1# and NCI-H460-SOX9 sh1# cells, as well as their corresponding control cells. Primers and probes for SOX9 and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) were designed using the Primer Express software, version 5.0 (Applied Biosystems, Foster City, CA, USA). Relative fold expression was calculated with the comparative threshold cycle (2−ΔΔCt) method.
The sequences of the primers used were as follows:
SOX9 Forward: 5ʹ-CAAGAAGGACCACCCGGATT-3ʹ
SOX9 Reverse: 5ʹ-AAGATGGCGTTGGGGGAGAT-3ʹ
Axin2 Forward: 5ʹ-CAACACCAGGCGGAACGAA-3ʹ
Axin2 Reverse: 5ʹ-GCCCAATAAGGAGTGTAAGGACT-3ʹ
Slug Forward: 5ʹ-CGAACTGGACACACATACAGTG-3ʹ
Slug Reverse: 5ʹ-CTGAGGATCTCTGGTTGTGGT-3ʹ
Twist Forward: 5ʹ-GCCTAGAGTTGCCGACTTATG-3ʹ
Twist Reverse: 5ʹ-TGCGTTTCCTGTTAAGGTAGC-3ʹ
Snail Forward: 5ʹ-TCGGAAGCCTAACTACAGCGA-3ʹ
Snail Reverse: 5ʹ-AGATGAGCATTGGCAGCGAG-3ʹ
CDH1 Forward: 5ʹ-ATTTTTCCCTCGACACCCGAT-3ʹ
CDH1 Reverse: 5ʹ- TCCCAGGCGTAGACCAAGA-3ʹ
MMP9 Forward: 5ʹ-TGTACCGCTATGGTTACACTCG-3ʹ
MMP9 Reverse: 5ʹ-GGCAGGGACAGTTGCTTCT-3ʹ
GAPDH Forward: 5ʹ-GACTCATGACCACAGTCCATGC-3ʹ
GAPDH Reverse: 5ʹ-AGAGGCAGGGATGATGTTCTG-3ʹ
Western blotting was performed as previously described . Briefly, protein concentrations were determined by the BCA assay (Pierce, Rockford, USA), and equal amounts of protein were separated on 10% sodium dodecyl sulfate–polyacrylamide gels and then transferred onto polyvinylidene fluoride membranes (Millipore, Bedford, USA). The membranes were probed with a primary antibody and a secondary antibody. Expression was determined by using the Super Signal West Pico Chemiluminescent Substrate (#34580; Thermo Fisher Scientific, USA) according to the manufacturer’s instructions. The membranes were stripped and re-probed with anti-β-tubulin antibody as a loading control. The antibodies used included: anti-SOX9 (1:5000, #ab182579; Abcam), anti-E-cadherin (1:1000, #610181; BD Biosciences), anti-γ-catenin (1:1000, #610254; BD Biosciences), anti-N-cadherin (1:1000, #610920; BD Biosciences), anti-vimentin (1:1000, #550513; BD Biosciences), anti-β-catenin (1:1000, #610154; BD Biosciences), anti-GAPDH (1:1000, sc-47724; Santa Cruz Biotechnology), anti-GSK-3β (1:1000, #12456; CST), anti-GSK-3β pSer9 (1:1000, #5558; Cell Signaling Technology), anti-β-tubulin (1:5000, #2146; Cell Signaling Technology), anti-H2A (1:1000, #12349; Cell Signaling Technology), goat anti-mouse secondary antibody (1:5000, #DC02L; Merck), and goat anti-rabbit secondary antibody (1:5000, #31460; Thermo Fisher Scientific, USA).
Cell migration was analyzed in vitro using 8.0 μm Transwell Permeable Supports (Corning, New York, USA). Briefly, 5 × 104 A549-pSin-Vector, A549-SOX9, NCI-H460-pSin-Vector, NCI-H460-SOX9, A549-pSuper-Vector, A549-SOX9 sh1#, NCI-H460-pSuper-Vector and NCI-H460-SOX9 sh1# cells in serum-free medium (100 μl) were seeded into the upper part of a Transwell chamber. The lower chamber was supplemented with 650 μl DMEM medium containing 20% FBS, and the chamber was incubated at 37 °C for 24 h. Cells on the upper side of the membrane were removed with cotton swabs, and cells on the lower side of the membrane were fixed with methanol for 10 min and then stained with crystal violet for 5 min. After the membrane was washed with distilled water, the cells that had invaded the membrane were counted in three randomly chosen fields of view per Transwell, and images were captured under an Olympus® CKX53 microscope. A similar procedure using a Transwell insert pre-coated with 80 μl Matrigel was performed to measure cell invasion ability in three independent assays. Comparisons between each two corresponding groups were made by two-tailed paired Student’s t test. Differences were considered statistically significant at p < 0.05 and three independent experiments were performed.
Wound healing assay
A549-pSin-Vector, A549-SOX9, NCI-H460-pSin-Vector, NCI-H460-SOX9, A549-pSuper-Vector, A549-SOX9 sh1#, NCI-H460-pSuper-Vector and NCI-H460-SOX9 sh1# cells at a density of 1 × 106 were seeded into six-well plates and cultured in DMEM containing 10% FBS to 90% confluence. The confluent cell monolayer was wounded using a sterile 200-μl pipette tip, and the cells in suspension were washed in normal growth medium. Images of the monolayer wound were captured after 0 h, 24 h and 48 h under an Olympus® CKX53 microscope in three randomly chosen fields of view. The migratory ability of the cells was calculated as the ratio of the open area after 24 h and 48 h to the open area at 0 h. Comparisons between each two corresponding groups were made by two-tailed paired Student’s t-test. Three independent assays were performed. Differences were considered statistically significant at p < 0.05.
Zebra fish assay
Transgenic Tg (flil: EGFP) zebrafish (provided by the Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, China) were maintained at 28 °C in aquaria under a 14 h/10 h day/night cycle. A549-pSin-Vector, A549-SOX9, NCI-H460-pSin-Vector, NCI-H460-SOX9, A549-pSuper-Vector, A549-SOX9 sh1#, NCI-H460-pSuper-Vector and NCI-H460-SOX9 sh1# cells were labeled with DiI stain. A total of 500 cells of each cell line were injected into the perivitelline space of 48-h post-fertilization embryos using a micro-injector (World Precision Instruments Inc., Sarasota, FL, USA). The embryos were incubated at 28 °C for 3 days before live images were captured under an Olympus® BX51 microscope. The number of disseminated foci from the tumor mass in the zebrafish embryo tails was counted. Comparisons between each two corresponding groups were made by two-tailed paired Student’s t-test. Differences were considered statistically significant at p < 0.05. Five zebrafish were used in each experiment and three independent experiments were performed.
Cells (A549-pSin-Vector, A549-SOX9, NCI-H460-pSin-Vector, NCI-H460-SOX9, A549-pSuper-Vector, A549-SOX9 sh1#, NCI-H460-pSuper-Vector and NCI-H460-SOX9 sh1#) for immunofluorescence staining were grown and treated in chamber slides, fixed in 4% formaldehyde in phosphate-buffered saline (PBS) for 10 min, permeabilized for 10 min with 0.2% Triton X-100 in PBS, and blocked with 2% bovine serum albumin (BSA) for 1 h. Primary antibodies against SOX9 (#ab182579, Abcam), E-cadherin (#610181, BD Biosciences), γ-catenin (#610254, BD Biosciences), N-cadherin (#610920, BD Biosciences), vimentin (#550513, BD Biosciences), and β-catenin (#610154, BD Biosciences) were diluted to 1:400 in PBS containing 1% BSA and incubated for 1 h at room temperature. Secondary antibody was purchased from Life Technologies® (Grand Island, NY), diluted to 1:250 in 1% BSA in PBS and incubated for 1 h. Images were captured with the Nikon® TS2 microscope.
Dual-luciferase reporter assay
Cells (A549-pSin-Vector, A549-SOX9, NCI-H460-pSin-Vector, NCI-H460-SOX9, A549-pSuper-Vector, A549-SOX9 sh1#, NCI-H460-pSuper-Vector and NCI-H460-SOX9 sh1#) were plated in 100-mm cell culture dishes, until they proliferated to 60–80% confluence after 24 h of culture. TOP flash or FOP flash and Renilla pRL-TK plasmids were transfected into the cells with Lipofectamine 3000 (Life Technologies) according to the manufacturer’s protocol. After 48 h of incubation, the transfection medium was replaced, and the cells were harvested, washed with PBS and lysed with passive lysis buffer (Promega). The cell lysates were analyzed immediately using a 96-well plate luminometer (Biotech, Germany). Luciferase and Renilla luciferase activity were measured using a dual-luciferase reporter assay system (Promega) according to the manufacturer’s instructions. The luciferase activity of each lysate was normalized to Renilla luciferase activity. The relative transcriptional activity was converted into fold induction above the vehicle control value .
The gene set enrichment analysis was performed according to a standard protocol as described previously [16, 17]. We first prepared four date files, an expression dataset file, phenotype labels file, gene sets file, and chip annotations file. Next, Input data to GSEA software and run the gene set enrichment analysis. To explore the potential pathways that may be involved in the metastatic effect of SOX9, the gene set enrichment analysis (GSEA, http://software.broadinstitute.org/gsea/index.jsp) software program was used to analyze the GEO database (GEO/GSE42127, http://software.broadinstitute.org/gsea/msigdb/index.jsp), which is the transcription profiling by array of NSCLC patients according to a standard protocol.
All statistical analyses were carried out using the statistical software package SPSS, version 22.0 (IBM SPSS, Chicago, USA). Comparisons between groups were performed using the two-tailed paired Student’s t-test. The relationship between SOX9 expression and TNM classification was assessed using the χ2 test. Bi-variate correlations between variables were calculated using Spearman’s rank correlation coefficients. A p-value of less than 0.05 was considered to indicate statistical significance in all cases.
Correlation of SOX9 overexpression with NSCLC TNM stage
We previously performed immunohistochemical staining with an antibody against human SOX9 and reported that SOX9 overexpression is significantly correlated to NSCLC histological stage (p = 0.017) and poor prognosis (p < 0.001) . Here, we investigated whether the SOX9 protein expression levels also correlate with NSCLC TNM stages. We classified the 142 paraffin-embedded, archived NSCLC clinical specimens according to TNM stage (Table 1 and Fig. 1a) and found that SOX9 expression levels were positively associated with the T (p = 0.030), N (p = 0.000), and M (p = 0.032) stages (Table 2 and Fig. 1a). Spearman rank analysis confirmed that there was a significant correlation between increasing SOX9 expression levels and progression from the T to the N and then the M stage (Table 3 and Fig. 1b). These results indicate that NSCLC metastasis is significantly associated with SOX9 overexpression.
Effect of SOX9 on the migration and invasion ability of NSCLC cell lines
To investigate the role of SOX9 played in NSCLC metastasis, we transduced A549 and NCI-H460 NSCLC cells with a pSin/SOX9-derived retroviral construct to establish stable SOX9-over-expressing cell lines (Fig. 2a) or pSuper/SOX9 shRNA-derived retroviral construct to establish SOX9-knock down cell lines (Fig. 2b).
A Transwell migration assay was used to determine the invasive and migratory ability of the SOX9-overexpressing cell lines A549-SOX9 and NCI-H460-SOX9 and the SOX9-knockdown cell lines A549-SOX9 sh1# and NCI-H460-SOX9 sh1#, as well as their corresponding control cells. The assay results showed that overexpression of SOX9 resulted in an increase in the number of cells that passed through the bottom of the Transwell chamber with or without the Matrigel layer, as compared to the control groups. Downregulation of SOX9 had converse effects. These data indicate that the invasive and migratory capacity of SOX9-overexpressing cells was enhanced, while the capacity of SOX9-knockdown cells was reduced in comparison to the control groups (Fig. 2c–h).
The wound healing assay was employed to detect the cell migration ability, and the results also showed that SOX9 overexpression could accelerate the migration rate as compared to the control cells (Fig. 2i–l). Contrasting results were obtained in the SOX9-knockdown cell lines as compared to the control cells.
Taken together, our data indicate that overexpression of SOX9 could improve the invasion ability of NSCLC cells and thereby promote their transformation into the malignant phenotype.
In vivo and in vitro effect of SOX9 on inducing EMT and promoting distant metastasis
Western blotting and immunofluorescence analyses revealed that expression of the epithelial cell markers (E-cadherin and γ-catenin) was downregulated and expression of the mesenchymal cell markers (vimentin and N-cadherin) was upregulated at the protein level in the SOX9-overexpressing cell lines as compared to the corresponding control cell lines (Fig. 3a, b). In contrast, expression of the epithelial cell markers (E-cadherin and γ-catenin) was upregulated and expression of the mesenchymal cell markers (vimentin and N-cadherin) was downregulated at the protein level in the SOX9-knockdown cells as compared to the corresponding control cells (Fig. 3c, d).
Next, we examined the in vivo effect of SOX9 over-expression on metastasis in a zebrafish model. We microinjected each constructed cell line and after 1 week found increased numbers of SOX9-over-expressing disseminated tumor foci (A549-SOX9 and NCI-H460-SOX9) in the zebrafish tail compared to the control groups (A549-pSin-Vector and NCI-H460-pSin-Vector). Conversely, we found decreased numbers of SOX9-knockdown disseminated tumor foci (A549-SOX9 sh1# and NCI-H460-SOX9 sh1#) in the zebrafish tail compared to the control group (Fig. 3e, f). These results indicate that SOX9 over-expression can promote in vivo distant metastasis of NSCLC.
Promotion of EMT by SOX9 via activation of the Wnt/β-catenin pathway in NSCLC cell lines
To investigate the mechanism of SOX9-induced EMT, we performed gene enrichment analysis with the GSEA software. We found that SOX9 was significantly associated with the Wnt/β-catenin signaling pathway (Fig. 4a). Therefore, RT-PCR analysis was used to quantify the mRNA levels of genes involved in the Wnt/β-catenin signaling pathway that are associated with metastasis, invasion, migration and EMT. Further, the RNA expression levels of downstream target genes of the Wnt/β-catenin pathway were detected by real-time PCR. The data showed that CDH1 (E-cadherin) expression was significantly decreased, and that Twist, snail, MMP9, Axin2, and Slug expression was significantly increased in SOX9-overexpressing cells as compared to the corresponding control cells. We also found that the expression of Twist, snail, MMP9, Axin2 and Slug was significantly decreased and the expression of CDH1 (E-cadherin) was significantly increased in the SOX9-knockdown cells as compared to the corresponding control cells (Fig. 4b).
In the classical Wnt/β-catenin pathway, β-catenin accumulates in the nucleus and binds to TCF1/LEF1 to promote downstream target gene expression. We therefore performed a dual luciferase reporter assay to detect TCF1/LEF1 transcriptional activity. The luciferase fluorescence signal versus Renilla fluorescence signal was increased in the SOX9-overexpressing A549-SOX9 and NCI-H460-SOX9 cells but was decreased in the SOX9-knockdown A549-SOX9 sh1# and NCI-H460-SOX9 sh1# cells, as compared to the corresponding control cells. Thus, upregulation of SOX9 expression enhanced TCF1/LEF1 transcriptional activity (Fig. 4c). This means that SOX9 might activate the Wnt/β-catenin pathway by promoting β-catenin translocation.
Next, we analyzed β-catenin protein expression in nuclear extracts from SOX9-overexpressing cells and SOX9-knockdown cells. The β-catenin expression in the nucleus was increased in SOX9-overexpressing cells and decreased in SOX9-knockdown cells (Fig. 4d). Further, immunofluorescence assay showed enhanced β-catenin translocation from the cytoplasm to the nucleus in the SOX9-overexpressing cells but not in the SOX9-knockdown cells (Fig. 4e). This finding indicates that SOX9 promotes the localization of β-catenin to the nucleus, where it acts as a co-activator for enhancing the expression of downstream target genes of the Wnt/β-catenin signaling pathway.
GSK3β is a key regulatory molecule in both the Wnt/β-catenin and the phosphatidylinositol 3-kinase/AKT pathways. Akt activated by PI3 K is known to have the capacity to inactivate GSK3β by phosphorylating it at residue Ser9. Once GSK3β is phosphorylated at residue Ser9, the complex formed by GSK3β, Axin and APC is destroyed; this is the mechanism responsible for the translocation of β-catenin to the nucleus. In this study, we found that the expression of GSK-3β was not significantly changed by SOX9 overexpression; however, the phosphorylation level of GSK-3β at Ser9 was markedly increased (Fig. 4f). Thus, activation of the Wnt/β-catenin pathway by SOX9 might be initiated by the phosphorylation of GSK3β at Ser9 and the subsequent nuclear translocation of β-catenin.
Effect of the β-catenin inhibitor XAV-939 on EMT induced by SOX9
Finally, we performed an assay to determine whether SOX9-induced EMT could be reversed by inhibiting β-catenin function with the β-catenin inhibitor XAV-939. We cultured A549-SOX9 and NCI-H460-SOX9 (SOX9-overexpression) cells in the presence of XAV-939 (100 μM) and evaluated EMT marker expression by western blotting. Inhibition of β-catenin expression was found to increase E-cadherin expression and decrease N-cadherin expression; thus, the EMT process might be reversed by the β-catenin inhibitor XAV-939 (Fig. 5a, b). The wound healing assay and Transwell assay also showed that XAV-939 application inhibits the migratory capacity of SOX9-overexpressin cells (Fig. 5c–e).
In the present study, we were able to confirm the role of SOX9 in the progression of NSCLC. The pathological data obtained from 142 NSCLC patients confirmed that SOX9 overexpression is significantly associated with TNM stage. We thus focused on the role of SOX9 in distant metastasis in NSCLC and its underlying molecular mechanism.
Distant metastasis is the last step in solid tumor progression and the highest risk factor for host death . Distant metastasis involves the detachment of tumor cells from their primary location or tumor, their invasion into circulation or the lymphatic system, and their proliferation at a location that is distant from the primary tumor . In thyroid cancer cells, knockdown of SOX9 has been found to inhibit proliferation, invasion, and EMT . We thus modeled this process by microinjecting fluorescently labeled SOX9-over-expressing or SOX9-knockdown cells into zebrafish larvae. We found that zebrafish injected with SOX9-over-expressing cells developed more disseminated tumor foci and zebrafish injected with SOX9-knockdown cells developed less disseminated tumor foci than control zebrafish. These data support our hypothesis that SOX9 overexpression promotes distant metastasis in NSCLC.
EMT is an early event in tumor metastasis. It is characterized by downregulation of epithelial cell markers (E-cadherin, γ-catenin and claudin-1) and upregulation of mesenchymal cell markers (vimentin, fibronectin, and N-cadherin) [21, 22]. In a human chondrocyte cell line, miR-19a was found to promote cell viability and migration of chondrocytes via positive regulation of SOX9 expression . Further, overexpression of miR-216b suppressed NSCLC cell proliferation and invasion by targeting SOX9. Interestingly, the effects of miR-216b can be mimicked by siRNA against SOX9 . Here, we found that SOX9 overexpression results in an increase in N-cadherin and vimentin expression and a decrease in E-cadherin and γ-catenin expression. The opposite results are obtained in SOX9 knockdown cells. These data indicate that SOX9 overexpression may promote NSCLC metastasis by inducing EMT.
Many signaling pathways are involved in EMT-mediated metastasis, including the TGF-β, MAPK, AKT, PI3 K, Notch and WNT pathways [25,26,27,28]. It was reported recently that SOX9 participates in the feedback loop of MALAT1 (metastasis associated lung adenocarcinoma transcript 1)-miR-101-SOX9, which modulates the chemoresistance of lung cancer cells to Cisplatin via the Wnt signaling pathway . In the present study, we found that SOX9 was significantly associated with the Wnt/β-catenin pathway, based on the results of GSEA with the GEO database (GEO/GSE42127). This was confirmed by the increased transcription of downstream metastasis-associated genes of the Wnt/β-catenin signaling pathway in the SOX9-overexpressing cells.
Our nuclear extraction assay and immunofluorescence analysis showed that SOX9 overexpression promoted β-catenin translocation from the cytoplasm to nucleus, but this was not observed in the SOX9-knockdown cells. This finding combined with the dual luciferase reporter assay finding that SOX9 overexpression indirectly enhanced TCF1/LEF1 transcriptional activity indicates that SOX9 promoted the translocation of β-catenin to the nucleus. Similarly, Blache et al. reported that SOX9 overexpression is the precursor of β-catenin accumulation in the nucleus, and that it acts on TCF4 to form a β-catenin–TCF4 complex in colon cancer . To further verify the role of SOX9 in inducing EMT by β-catenin nuclear translocation, we used the β-catenin inhibitor XAV-939. Our data showed that SOX9-induced EMT could be reversed upon exogenous application of XAV-939. Taken together, these results demonstrate that SOX9 promotes EMT by activating the Wnt/β-catenin pathway through the nuclear translocation of β-catenin.
A study on chondroplasia showed that the N-terminal of SOX9 binds to β-catenin and induces β-catenin translocate to the nucleus. Once SOX9 induced the nuclear translocation of β-catenin, the phosphorylation and degradation of β-catenin was accelerated by GSK-3β . In the NSCLC cells used in the present study, we found that the GSK-3β levels were not altered as a result of SOX9 overexpression; however, GSK-3β phosphorylation at serine 9 was significantly increased. Thus, SOX9 might promote metastasis in NSCLC by activating the Wnt/β-catenin pathway through enhancing the phosphorylation of GSK-3β at Ser9. The interaction between SOX9 and GSK-3β needs to be further investigated, especially with regard to whether SOX9 is associated closely with the reduction of β-catenin phosphorylation.
Our study confirms that SOX9 promotes distant metastasis by driving the EMT in NSCLC. We show that the Wnt/β-catenin pathway is activated as a result of SOX9 promoting β-catenin translocation from the cytoplasm to the nucleus to enhance TCF1/LEF1 transcriptional activity. This effect may be initiated by SOX9-mediated phosphorylation of GSK-3β at Ser9. Further studies on the relevance of SOX9-mediated GSK3β phosphorylation at Ser9 are now required to better understand the mechanisms of SOX9. Collectively, these findings suggest that manipulation of SOX9 to impede NSCLC metastasis and malignancy may be clinically useful and provide a new insight for the intervention of distal metastasis in cancer therapy.
protein kinase B
American Type Culture Collection
bovine serum albumin
complimentary deoxyribonucleic acid
cyclic guanosine monophosphate-dependent protein kinase II
Dulbecco Modified Eagle Medium
epidermal growth factor receptor
fetal bovine serum
glyceraldehyde 3-phosphate dehydrogenase
glycogen synthase kinase-3β
gene set enrichment analysis
high mobility group
lymphoid enhancer-binding factor
non-small cell lung cancer
phosphate buffered saline
polymerase chain reaction
real-time reverse transcription-polymerase chain reaction
small cell lung cancer
SRY-related High Mobility Group-box 9
sex-determining region Y
tris buffered saline
T cell transcription factor
tumor node metastasis
Gong L, Song J, Lin X, Wei F, Zhang C, Wang Z, Zhu J, Wu S, Chen Y, Liang J, et al. Serine-arginine protein kinase 1 promotes a cancer stem cell-like phenotype through activation of Wnt/beta-catenin signalling in NSCLC. J Pathol. 2016;240:184–96.
Travis LB, Curtis RE, Bennett WP, Hankey BF, Travis WD, Boice JD. Lung cancer after Hodgkin’s disease. J Natl Cancer Inst. 1995;87:1324–7.
Wagner T, Wirth J, Meyer J, Zabel B, Held M, Zimmer J, Pasantes J, Bricarelli FD, Keutel J, Hustert E, et al. Autosomal sex reversal and campomelic dysplasia are caused by mutations in and around the SRY-related gene SOX9. Cell. 1994;79:1111–20.
Gruber HE, Norton HJ, Ingram JA, Hanley EN. The SOX9 transcription factor in the human disc: decreased immunolocalization with age and disc degeneration. Spine (Phila Pa 1976). 2005;30:625–30.
Muller P, Crofts JD, Newman BS, Bridgewater LC, Lin CY, Gustafsson JA, Strom A. SOX9 mediates the retinoic acid-induced HES-1 gene expression in human breast cancer cells. Breast Cancer Res Treat. 2010;120:317–26.
Jiang SS, Fang WT, Hou YH, Huang SF, Yen BL, Chang JL, Li SM, Liu HP, Liu YL, Huang CT, et al. Upregulation of SOX9 in lung adenocarcinoma and its involvement in the regulation of cell growth and tumorigenicity. Clin Cancer Res. 2010;16:4363–73.
Wang H, Leav I, Ibaragi S, Wegner M, Hu GF, Lu ML, Balk SP, Yuan X. SOX9 is expressed in human fetal prostate epithelium and enhances prostate cancer invasion. Can Res. 2008;68:1625–30.
Lu B, Fang Y, Xu J, Wang L, Xu F, Xu E, Huang Q, Lai M. Analysis of SOX9 expression in colorectal cancer. Am J Clin Pathol. 2008;130:897–904.
Cortese R, Hartmann O, Berlin K, Eckhardt F. Correlative gene expression and DNA methylation profiling in lung development nominate new biomarkers in lung cancer. Int J Biochem Cell Biol. 2008;40:1494–508.
Chakravarty G, Moroz K, Makridakis NM, Lloyd SA, Galvez SE, Canavello PR, Lacey MR, Agrawal K, Mondal D. Prognostic significance of cytoplasmic SOX9 in invasive ductal carcinoma and metastatic breast cancer. Exp Biol Med. 2011;236:145–55.
Zhou CH, Ye LP, Ye SX, Li Y, Zhang XY, Xu XY, Gong LY. Clinical significance of SOX9 in human non-small cell lung cancer progression and overall patient survival. J Exp Clin Cancer Res. 2012;31:18.
Lefebvre V, Dumitriu B, Penzo-Mendez A, Han Y, Pallavi B. Control of cell fate and differentiation by Sry-related high-mobility-group box (Sox) transcription factors. Int J Biochem Cell Biol. 2007;39:2195–214.
Swartling FJ, Ferletta M, Kastemar M, Weiss WA, Westermark B. Cyclic GMP-dependent protein kinase II inhibits cell proliferation, Sox9 expression and Akt phosphorylation in human glioma cell lines. Oncogene. 2009;28:3121–31.
Edge SB, Compton CC. The American Joint Committee on Cancer: the 7th Edition of the AJCC cancer staging manual and the future of TNM. Ann Surg Oncol. 2010;17:1471–4.
Gong LY, Song JW, Lin X, Wei FK, Zhang CC, Wang ZM, Zhu JR, Wu S, Chen Y, Liang J, et al. Serine-arginine protein kinase 1 promotes a cancer stem cell-like phenotype through activation of Wnt/beta-catenin signalling in NSCLC. J Pathol. 2016;240:184–96.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005;102:15545–50.
Reimand J, Isserlin R, Voisin V, Kucera M, Tannus-Lopes C, Rostamianfar A, Wadi L, Meyer M, Wong J, Xu CJ, et al. Pathway enrichment analysis and visualization of omics data using g:profiler, GSEA. Cytoscape and EnrichmentMap. Nat Protocols. 2019;14:482–517.
Weigelt B, Peterse JL, Veer LJ. Breast cancer metastasis: Markers and models. Nat Rev Cancer. 2005;5:591–602.
Mehes G, Witt A, Kubista E, Ambros PF. Circulating breast cancer cells are frequently apoptotic. Am J Pathol. 2001;159:17–20.
Huang J, Guo L. Knockdown of SOX9 inhibits the proliferation, invasion, and EMT in thyroid cancer cells. Oncol Res. 2017;25:167–76.
Polyak K, Weinberg RA. Transitions between epithelial and mesenchymal states: acquisition of malignant and stem cell traits. Nat Rev Cancer. 2009;9:265–73.
Thiery JP. Epithelial–mesenchymal transitions in tumour progression. Nat Rev Cancer. 2002;2:442–54.
Yu C, Wang Y. MicroRNA-19a promotes cell viability and migration of chondrocytes via up-regulating SOX9 through NF-kappaB pathway. Biomed Pharmacother. 2018;98:746–53.
Liu S, Dong H, Dai H, Liu D, Wang Z. MicroRNA-216b regulated proliferation and invasion of non-small cell lung cancer by targeting SOX9. Oncol Lett. 2018;15:10077–83.
Bakin AV, Tomlinson AK, Bhowmick NA, Moses HL, Arteaga CL. Phosphatidylinositol 3-kinase function is required for transforming growth factor beta-mediated epithelial to mesenchymal transition and cell migration. J Biol Chem. 2000;275:36803–10.
Bhowmick NA, Ghiassi M, Bakin A, Aakre M, Lundquist CA, Engel ME, Arteaga CL, Moses HL, Bhowmick NA, Ghiassi M, Bakin A, Aakre M, Lundquist CA, Engel ME, et al. Transforming growth factor-b1 mediates epithelial to mesenchymal transdifferentiation through a RhoA-dependent mechanism. Mol Biol Cell. 2001;12:27–36.
Smith BN, Burton LJ, Henderson V, Randle DD, Morton DJ, Smith BA, Taliaferro-Smith L, Nagappan P, Yates C, Zayzafoon M, et al. Snail promotes epithelial mesenchymal transition in breast cancer cells in part via activation of nuclear ERK2. PLoS ONE. 2014;9:1.
Zhu Z, Dai J, Liao Y, Wang T. Sox9 protects against human lung fibroblast cell apoptosis induced by LPS through activation of the AKT/GSK3beta pathway. Biochemistry (Mosc). 2017;82:606–12.
Chen W, Zhao W, Zhang L, Wang L, Wang J, Wan Z, Hong Y, Yu L. MALAT1-miR-101-SOX9 feedback loop modulates the chemo-resistance of lung cancer cell to DDP via Wnt signaling pathway. Oncotarget. 2017;8:94317–29.
Blache P, van de Wetering M, Duluc I, Domon C, Berta P, Freund JN, Clevers H, Jay P. SOX9 is an intestine crypt transcription factor, is regulated by the Wnt pathway, and represses the CDX2 and MUC2 genes. J Cell Biol. 2004;166:37–47.
Topol LL, Chen W, Song H, Day TF, Yang YZ. Sox9 inhibits Wnt signaling by promoting beta-catenin phosphorylation in the nucleus. J Biol Chem. 2009;284:3323–33.
JQH, WFK and XLX performed the experiments. SXY and JWS analyzed the data. PKD, JZ, HFL and XPL contributed reagents/materials/analysis tools. HG and LS helped in data collection. LYG and LY conceived of the study, participated in its design and wrote the paper. All authors read and approved the final manuscript.
The authors would like to thank Dr. Jessica Tamanini (ETediting and Shenzhen University) for her critical review and editing of the manuscript prior to submission.
The authors declare that they have no competing interests.
Availability of data and materials
All data are fully available without restriction. The datasets generated and analysed during the current study are available in the [GEO/GSE42127] repository, (http://software.broadinstitute.org/gsea/msigdb/index.jsp).
Consent for publication
Ethics approval and consent to participate
All subjects have provided written informed consent before participation in the study which has been approved by the Ethics Committee of the Guangzhou Medical University.
This work was supported by the Natural Science Foundation of China (81101501, 81872448, and 31761133012), the science and technology program of Guangdong province in China (2017B030301016), the Science and Technology Bureau of Shenzhen City (JCYJ 20130401144744187, 20170302144320949, 2017041211200974 and 20180305163414905), the Science and Technology Foundation of Shenzhen Nanshan district (2018001), Discipline Construction Funding of Shenzhen (2016-1452) and Medicine Foundation of Guangdong Province (A2017149 and A2018265).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.