Skip to main content
Download PDF

Systematic analysis of genes involved in oral cancer metastasis to lymph nodes

Cellular & Molecular Biology Letters volume 23, Article number: 53 (2018) Cite this article

Abstract

Oral cancer remains a deadly disease worldwide. Lymph node metastasis and invasion is one of the causes of death from oral cancer. Elucidating the mechanism of oral cancer lymph node metastasis and identifying critical regulatory genes are important for the treatment of this disease. This study aimed to identify differentially expressed genes (gene signature) and pathways that contribute to oral cancer metastasis to lymph nodes. The GSE70604-associated study compared gene profiles in lymph nodes with metastasis of oral cancer to those of normal lymph nodes. The GSE2280-associated study compared gene profiles in primary tumor of oral cancer with lymph node metastasis to those in tumors without lymph node metastasis. There are 28 common differentially expressed genes (DEGs) showing consistent changes in both datasets in overlapping analysis. GO biological process and KEGG pathway analysis of these 28 DEGs identified the gene signature CCND1, JUN and SPP1, which are categorized as key regulatory genes involved in the focal adhesion pathway. Silencing expression of CCND1, JUN and SPP1 in the human oral cancer cell line OECM-1 confirmed that those genes play essential roles in oral cancer cell invasion. Analysis of clinical samples of oral cancer found a strong correlation of these genes with short survival, especially JUN expression associated with metastasis. Our study identified a unique gene signature – CCND1, JUN and SPP1 – which may be involved in oral cancer lymph node metastasis.

Introduction

Oral cancer is a major devastating head and neck cancer subtype and is any cancerous tissue growth located in the oral cavity or oropharynx [1]. There are several types of oral cancers, but around 90% of oral cancers are squamous cell carcinomas. Oral cancer is the sixth most common malignancy in humans. Its incidence and mortality have also increased over the past decades. In 2013, oral cancer resulted in 135,000 deaths, up from 84,000 deaths in 1990 [2]. Oral cancer is characterized by poor prognosis and a low survival rate despite sophisticated surgical and radiotherapeutic modalities. The five-year survival rate is only 63% in the United States.

Metastasis is one of the main causes of oral cancer patients’ death. The process of oral cancer metastasis is a series of sequential and interdependent events involving detachment of cells from tumor tissue, increased cell motility and local invasion, angiogenesis, intravasation of invading cells into the vasculature or lymphatic systems, extravasation and subsequent deposition and proliferation at a second site. Oral cancer tends to spread primarily to the regional lymph nodes of the neck first before it spreads to remote sites. Lymph node metastasis is called locoregional metastasis and lung (or other organs) metastasis is called distant metastasis [3]. Early diagnosis of lymph node metastasis is important for improving clinical outcomes of oral cancer patients [4].

The clinical diagnosis of lymph node metastasis of oral cancer is currently based on imaging techniques and sentinel lymph node biopsy (SLNB) [5]. However, current imaging tests have been proven to be unreliable, especially in the detection of early nodal diseases. Localization of SLNB in patients with certain cancers such as floor-of-mouth carcinoma is difficult. The difficulties of detecting micrometastasis in frozen sections also limit SLNB as an adequate guide to clinical decision-making [6]. Therefore, better understanding of the lymphatic metastasis of oral cancer and developing new diagnostic strategies to predict the clinical behavior of the disease are desired [4].

Identifying reliable gene signatures or molecular biomarkers for oral cancer lymph node metastasis is highly valuable to use as a potential diagnostic or prognostic tool for routine clinical practice of oral cancer. Microarray analysis has been extensively used to examine the global gene expression changes in non-metastatic versus metastatic human tumor samples. Comparing the gene expression profiles of primary tumors from non-metastatic and metastatic tumor generated a list of genes whose expression was significantly different in these two groups. The list of genes is called “metastasis signatures”. Several “metastasis signatures” have been generated in oral cancer aimed at predicting lymph node metastasis [7,8,9,10,11,12,13,14,15,16,17]. However, defined “metastasis signatures” from different studies are highly variable, which hinders their translation to clinical applications. To date, none of the identified signatures or molecular biomarkers has been successfully implemented as a diagnostic or prognostic tool applicable to routine clinical practice of oral cancer.

In this study, in order to identify differentially expressed genes (gene signature) and pathways that contribute to oral cancer metastasis to lymph nodes, we designed a new bioinformatic analysis strategy. The GSE70604-associated study is aimed at identifying lymph node metastasis-associated genes and pathways in oral cancer by comparing gene profiles in lymph nodes with metastasis to those of normal lymph nodes. The GSE2280-associated study tried to achieve this goal by comparing gene profiles in squamous cell carcinomas of the oral cavity with lymph node metastasis to those in tumors without lymph node metastasis. Then overlapping analysis was performed to identify common DEGs (gene signature) and pathways in these two studies. We focused on those 28 common DEGs that show consistent changes in both datasets and defined a “metastasis signature”: CCND1, JUN and SPP1 for oral cancer lymph node metastasis. Both Jun and SPP1 are known molecules involved in cancer cell invasiveness and tumor metastasis. CCND1 (cyclinD1) is a known critical factor in cell cycle progression, but the function of CCND1 in human cancer cell migration/invasion is not sufficiently understood despite a recent study showing that it directly regulates the focal adhesion pathway and promotes R3327 rat prostatic tumor cell migration/invasion and tumor metastasis [18] and affects cell migration and invasion in breast cancer [18,19,20,21]. Our study identified a novel “focal adhesion”-related gene signature (CCND1, JUN and SPP1) that might be applicable for diagnosis of oral cancer metastasis to lymph nodes.

Materials and methods

Cell lines

OECM-1 is a human oral cancer cell line derived from gingival epidermoid carcinoma of a patient [22]. The cells were cultured in RPMI1640 supplemented with 10% fetal bovine serum (FBS), 2 mmol/ l glutamine, 100 units/mL penicillin and 100 μg/mL streptomycin in a humidified atmosphere of 95% air and 5% CO2 at 37 °C.

siRNA interference

SPP1, JUN and CCND1 siRNA smart pool (Cat. L-009152-00-0005) was purchased from Dharmacon/Thermo Fisher Scientific. Transfection of the siRNA oligonucleotide duplexes was performed in a 6-well plate (1 × 10 5 cells per well) with Lipofectamine 2000 (Invitrogen, Inc.), using the methods recommended by the manufacturer. Knockdown of SPP1, JUN and CCND1 with siRNA was examined 72 h after siRNA transfection through western blotting.

Western blot

Cells were lysed in a lysis buffer containing 50 mmol/L TRIS-HCl, pH 7.4, 150 mmol/L NaCl, 0.5% NP40, 50 mmol/L NaF, 1 mmol/L Na3VO4, 1 mmol/L phenyl-methylsulfonyl fluoride, 25 μg/mL leupeptin, and 25 μg/mL aprotinin and clarified by centrifugation (14,000 g for 30 min at 4 °C). The protein concentration of the cell lysates was determined using the Bradford Coomassie blue method (Pierce Chemical Corp.). Whole-cell lysates were separated by sodium dodecyl sulfate (SDS)-PAGE and transferred onto nitrocellulose membrane. The membranes were blocked with PBS containing 5% (w/v) skim milk at 4 °C for 2 h, washed with PBST (PBS with 0.05% Tween-20), and then incubated overnight with primary antibody. After washing with PBST, the membrane was incubated with a second antibody at room temperature for 2 h, washed with PBST and then developed with the ECL system. The results of Western blot were analyzed with Odyssey software version 3.0.

Cell invasion assay

Cell invasion potential was measured with a Boyden transwell chamber consisting of upper inserts with 8-μm-pore-size filter membranes at the bottom of the inserts and lower wells in 24-well cell culture plates (Corning Life Sciences). 20 μl of 1:6 diluted Matrigel (2–3 mg/ml protein) was added to the center of each cell well insert. Coated inserts were placed in an incubator to allow the Matrigel to solidify for 20–30 min. Cells (3.5 × 105 cells in 0.2 mL) suspended in serum-free medium with 0.1% bovine serum albumin were seeded into the inserts of the chambers. The inserts were then placed over the wells filled with 0.5 mL of 10% FBS culture medium and incubated in a 37 °C incubator for 24 h. Cells that had not penetrated the filter membrane in the inserts were wiped off with cotton swabs, and the cells on the underside of the filter membrane were fixed and stained with the HEMA-3 kit (Fisher Diagnostics). Invaded cells were counted in a total of 10 fields for each sample under a microscope with the 10X objective and the stained cell number per field was calculated.

Differentially expressed genes (DEG) analysis

Differentially expressed genes of the GSE2280-associated study and GSE2280-associated study were analyzed. The expression abundance (FPKM) value of each gene was estimated by running cufflinks [23] and the differentially expressed genes were assessed by cuffdiff. Statistically differentially expressed genes between two groups were those genes with logFC (fold change) > 0.45 and adjusted p value < 0.05. The adjusted p value was obtained through applying Benjamini and Hochberg’s (BH) false discovery rate correction on the original p value, and the fold change threshold was selected based on our purpose of focusing on significantly differentially expressed genes.

Hierarchical clustering

Hierarchical clustering was conducted [24] to classify analyzed samples based on gene expression profiles. Hierarchical clustering using differentially expressed genes (DEGs) demonstrated the global gene expression patterns in the samples. In addition, the DEGs were further extracted and classified in specific biological processes (Gene Ontology terms) and KEGG pathways. The expression pattern of those DEGs was characterized and heat maps of the DEGs were classified in targeted biological processes or KEGG pathways using the R package.

GO and KEGG pathway analysis

We used the R packages GO.db, KEGG.db and KEGGREST to detect Gene Ontology categories and KEGG pathways with significant enrichment in DEGs for comparison across all measured genes. The significantly enriched biological processes were identified by p value less than the threshold value 0.05. For the KEGG pathway, the p value was also set to less than 0.05.

Results

Identification of potential genes related to oral cancer metastasis to lymph nodes through screening GEO database

In order to find the key genes regulating oral cancer metastasis to lymph nodes, we screened the GEO (gene expression omnibus) database for GSE70604 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE70604) [25] and GSE2280 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE2280) [12] as shown in Additional file 1: Table S1 and Additional file 2: Table S3.

In GSE70604, the comparisons of gene expression profiles were made between lymph nodes with metastasis of oral squamous cell carcinoma (OSCC) and normal lymph nodes (comparison 1).

In GSE2280, the comparisons of gene expression profiles were made between primary tumors of OSCC which has lymph node metastasis and nonmetastatic primary OSCC without lymph node metastasis (comparison 2).

Differentially expressed genes (DEGs) of both comparisons were obtained where lymph nodes with OSCC metastasis were compared to normal lymph nodes in comparison 1 and metastatic OSCC primary tumors were compared to non-metastatic OSCC primary tumors in comparison 2. Both comparisons had the |log(fold change)|(logFC) > 0.45 and p value < 0.05, indicating the overall changes as statistically significant. In comparison 1, gene expression of 7 lymph nodes with metastasis of OSCC was compared to that of a normal lymph node. Figure 1 shows the distribution of DEGs in comparison 1; we found that 1392 genes had expression changes (Additional file 1: Table S1). Among those genes, we identified 723 down-regulated genes (Additional file 3: Table S2) and 699 up-regulated genes (Additional file 2: Table S3). In comparison 2, gene expression of 5 primary tumors of OSCC with lymph node metastasis was compared to that of 8 non-metastatic primary tumors of OSCC. Figure 2 shows the distribution of DEGs in comparison 2; we found 890 genes that had expression changes (Additional file 4: Table S4). Among those genes, we identified 477 down-regulated genes (Additional file 5: Table S5) and 413 up-regulated genes (Additional file 6: Table S6). To identify potential genes related to oral cancer metastasis to lymph nodes, we then analyzed the overlapping DEGs in comparison 1 and comparison 2 and found 114 overlapping DEGs. Among those DEGs, 28 of 114 DEGs had a consistent change trend (consistent increase or decrease in both comparisons) (Additional file 7: Table S7) and the remaining 86 had an opposite change trend (Additional file 8: Table S8). The 28 genes with consistent changes in comparison 1 and comparison 2 are the potential genes related to oral cancer lymph node metastasis.

Fig. 1
figure 1

Heat map of distribution of DEGs in GSE70604. The comparisons of gene expression profiles were made between lymph nodes with metastasis of oral squamous cell carcinoma (OSCC) and normal lymph nodes

Full size image
Fig. 2
figure 2

Heat map of distribution of DEGs in GSE2280. The comparisons of gene expression profiles were made between primary tumors of squamous cell carcinomas of the oral cavity (OSCC) which has lymph node metastasis and nonmetastatic primary OSCC without lymph node metastasis

Full size image

GO-KEGG analysis of consistent overlapping DEGs

We further analyzed the 28 consistent overlapping DEGs by GO-KEGG analysis. GO biological process analysis indicated that these 28 genes were involved in 448 biological processes; Additional file 9: Table S9 shows the top 10 biological processes. KEGG pathway analysis showed that those genes mainly participated in 12 signaling pathways (Additional file 10: Table S10) and the top 3 KEGG signaling pathways are colorectal cancer, Toll-like receptor signaling pathway and Chagas disease. We also identified 7 interacting gene pairs among these 28 overlapping genes (Additional file 11: Table S11). Among those genes, JUN, ATF3 and FOS were presented as the key connecting nodes in the network of signaling pathways according to connectivity. JUN and FOS cooperatively participated in 9 common KEGG signaling pathway among the 12 KEGG pathways we identified. JUN, ATF3 and FOS were involved in 257, 115 and 183 GO-biological processes respectively.

Definition of key regulatory genes involved in oral cancer metastasis to lymph nodes

Among 12 KEGG pathways we identified by analyzing the 28 consistent overlapping DEGs (Additional file 10: Table S10), we mainly focus on the focal adhesion pathway, which is important for cancer cell migration/invasion and tumor metastasis.

Three DEGs associated with the focal adhesion pathway in the 28 consistent overlapping DEGs are CCND1, JUN and SPP1. Both JUN and SPP1 are genes known to promote cancer cell invasiveness and tumor metastasis. CCND1 (cyclinD1), a critical gene regulating cell cycle progress, recently has been reported to directly regulate the focal adhesion pathway and promote R3327 rat prostatic tumor cell migration/invasion and tumor metastasis [18]. The direct involvement of CCND1 in human cancer cell migration/invasion is virtually unknown. Our study identified a novel focal adhesion-related gene signature (CCND1, JUN and SPP1) that might be important for oral cancer metastasis to lymph nodes.

Gene signature CCND1, JUN and SPP1 are required for oral cancer cell invasion

To confirm the function of CCND1, JUN and SPP1 in oral cancer cell invasion, we knocked down expression of CCND1, JUN and SPP1 through siRNA in the human oral cancer cell line OECM-1 and examined the alteration of cell invasion ability. Deletion of CCND1, JUN and SPP1 expression in OECM-1 via siRNA interference significantly impaired cell invasion although in varying degrees (Fig. 3a-c). The result further supports the hypothesis that the focal adhesion-related gene signature (CCND1, JUN and SPP1) is important for oral cancer invasion.

Fig. 3
figure 3

CCND1, JUN and SPP1 are involved in oral cancer cell invasion. a Human oral cancer cell line OECM-1 was transfected with SPP1 siRNA, SPP1 protein expression was examined through western blot, and alteration of cell invasion ability was examined and quantified as a histogram. b Human oral cancer cell line OECM-1 was transfected with JUN siRNA, JUN protein expression was examined through western blot, and alteration of cell invasion ability was examined and quantified as a histogram. c Human oral cancer cell line OECM-1 was transfected with CCND1 siRNA, CCND1 protein expression was examined through western blot, and alteration of cell invasion ability was examined and quantified as a histogram

Full size image

C- JUN expression is correlated with oral cancer metastasis in clinical practice

We further tried to explore the correlation of expression of CCND1, JUN and SPP1 with oral cancer metastasis through statistically analyzing the clinic data in the literature. Based on the original data of Rickman’s work published in 2008 in Oncogene [26], we compared JUN expression in oropharyngeal squamous cell carcinoma in patients without metastatic events at 5 years with that of patients with metastatic events at 5 years [26]. Metastatic oropharyngeal squamous cell carcinoma had 1.69-fold higher JUN gene expression than non-metastatic oropharyngeal squamous cell carcinoma (p = 0.025). Figure 4 compares JUN expression in oral cavity squamous cell carcinoma in patients without metastatic events at 3 years with that of patients with metastatic events at 3 years [26]. Metastatic oral cavity squamous cell carcinoma had 2.69-fold higher JUN gene expression than non-metastatic oral cavity squamous cell carcinoma (p = 0.012) (Fig. 4). Figure 5 compared JUN expression in oral cavity squamous cell carcinoma in patients without metastatic events at 5 years than that of patients with metastatic events at 5 years [26]. Metastatic oral cavity squamous cell carcinoma had 2.39-fold higher JUN gene expression than non-metastatic oral cavity squamous cell carcinoma (p = 0.025) (Fig. 5). Figure 6 compares JUN expression in oropharyngeal squamous cell carcinoma in patients without metastatic events at 5 years with that of patients with metastatic events at 5 years [26]. Metastatic oropharyngeal squamous cell carcinoma had 1.69-fold higher JUN gene expression than non-metastatic oropharyngeal squamous cell carcinoma (p = 0.025) (Fig. 6). The results suggest that c-JUN is strongly correlated with oral cancer metastasis. We have not found similar gene expression data of CCND1 and SPP1 in oral cancer metastasis. However, we observed that higher expression of Jun or CCND1 or SPP1 was associated with short survival of oral squamous cell carcinoma patients (Figs. 7, 8 and 9).

Fig. 4
figure 4

Comparison of JUN expression in oral cavity squamous cell carcinoma without metastatic events at 3 years and with metastatic events at 3 years

Full size image
Fig. 5
figure 5

Comparison of JUN expression in oral cavity squamous cell carcinoma without metastatic events at 5 years and with metastatic events at 5 years

Full size image
Fig. 6
figure 6

Comparison of JUN expression in oropharyngeal squamous cell carcinoma without metastatic events at 5 years and with metastatic events at 5 years

Full size image
Fig. 7
figure 7

Comparison of JUN expression in oropharyngeal squamous cell carcinoma

Full size image
Fig. 8
figure 8

Comparison of CCND1 expression in hypopharyngeal squamous cell carcinoma

Full size image
Fig. 9
figure 9

Comparison of SPP1 expression in oropharyngeal squamous cell carcinoma

Full size image

Discussion

Through bioinformatics analysis of the gene expression profile of oral cancer metastatic and non-metastatic lymph nodes and primary tumors, we identified a new oral cancer metastatic gene signature: CCND1, JUN and SPP1. The gene knockdown experiment on the oral cancer cell line and clinical data correlation analysis indicated strong association of this gene signature expression with oral cancer invasion and metastasis.

c-Jun is a subunit of transcription factor AP-1. Activated AP-1 increases the transcription of target genes and plays roles in cell division, proliferation, differentiation, apoptosis and so on. Although overexpression of c-Jun promotes invasion and metastasis of various tumors [27,28,29], it has not been reported in oral cancer. c-Jun plays a role during the initiation and progression of OSCC. High expression of c-Jun is associated with poor prognosis of OSCC [30]. Through a gene knockdown experiment in oral cancer cell line and clinical data analysis, our study confirmed that it is also required for oral cancer invasion and metastasis.

SPP1 (secreted phosphoprotein 1) encodes osteopontin (OPN). Osteopontin is a secreted non-collagenous, sialic acid rich, chemokine-like, matricellular phosphoglycoprotein that facilitates cell–matrix interactions and promotes tumor progression [31]. Osteopontin is a multifunctional cytokine regulating cell proliferation, survival, drug resistance, invasion, and stem-like behavior. Its aberrant expression and/or splicing is functionally responsible for many disease pathologies including cancer [31]. The expression of OPN was elevated in 95% of OSCC and can be used as a diagnostic marker for oral cancer [32]. It is also a poor prognostic factor in OSCC treated with cisplatin-based IC followed by CCRT [33]. However, its correlation with OSCC invasion and metastasis has not been studied anywhere else. The gene knockdown experiment in an oral cancer cell line and clinical data analysis showed that SPP1 is important for oral cancer invasion and metastasis in this study. Besides aberrant expression, the pattern of isoform expression (gene splicing) and post-translational modification are other SPP1 regulation methods. This regulation is cell-type specific and may influence the potential role of OPN in malignancy as a cancer biomarker [34] For example, invasive breast tumor cells generate three splice variants of OPN, while non-invasive breast cells express only the unspliced form or no OPN at all [35]. We do not know the splice variants of SPP1 in non-metastatic oral cancer and metastatic oral cancer with lymph node metastasis beside the expression difference, which deserves further investigation.

A recent study showed that CCND1 (cyclin D1) together with its binding partner CDK4 does not only simply act as a transcriptional regulator to control cell proliferation, but also controls cell adhesion, migration and metastasis under normal and pathological conditions. The focal adhesion component paxillin is a cytoplasmic substrate of CCND1·CDK4. This complex phosphorylates a fraction of paxillin specifically associated with the cell membrane, and promotes Rac1 activation, thereby triggering membrane ruffling and cell invasion in both normal fibroblasts and tumor cells [18]. Cytoplasmic CCND1 controls the migration and invasiveness of mantle lymphoma cells [36]. However, more evidence is needed to support the function of CCND1 in cell invasion and metastasis, especially in clinic samples. CCND1 is one of the DEGS with elevated expression in both metastatic lymph nodes and metastatic primary tumor of oral cancer in our analysis. The results suggested the important function of CCND1 in promoting oral cancer lymph node metastasis although the detained mechanism still needs to be explored.

In sum, this study identifies a unique gene signature – CCND1, JUN and SPP1 – which could be a new early biomarker for diagnosing oral cancer lymph node invasion and metastasis.

Abbreviations

CCND1:

Cyclin D1

DEGs:

Differentially expressed genes

GEO:

Gene expression omnibus

OPN:

Osteopontin

OSCC:

Oral squamous cell carcinoma

SLNB:

Sentinel lymph node biopsy

SPP1:

Secreted phosphoprotein 1

References

  1. Werning JW. Oral Cancer: Diagnosis, Management, and Rehabilitation2007 May 16; 2007.

    Google Scholar 

  2. Collaborators GMaCoD. Global, regional, and national age-sex specific all-cause and cause-specific mortality for 240 causes of death, 1990-2013: a systematic analysis for the global burden of disease study 2013. Lancet. 2015;385(9963):117–71.

    Article  Google Scholar 

  3. Noguti J, De Moura CF, De Jesus GP, Da Silva VH, Hossaka TA, Oshima CT, et al. Metastasis from oral cancer: an overview. Cancer Genomics Proteomics. 2012;9(5):329–35.

    CAS  PubMed  Google Scholar 

  4. Liu X, Kolokythas A, Wang J, Huang H, Zhou X. Gene expression signatures of lymph node metastasis in Oral Cancer: molecular characteristics and clinical significances. Curr Cancer Ther Rev. 2010;6(4):294–307.

    Article  CAS  Google Scholar 

  5. Civantos F, Zitsch R, Bared A. Sentinel node biopsy in oral squamous cell carcinoma. J Surg Oncol. 2007;96(4):330–6.

    Article  Google Scholar 

  6. Civantos F Jr, Zitsch R, Bared A, Amin A. Sentinel node biopsy for squamous cell carcinoma of the head and neck. J Surg Oncol. 2008;97(8):683–90.

    Article  Google Scholar 

  7. Kashiwazaki H, Hassan NM, Hamada J, Moriuchi T, Yamazaki Y, Tei K, et al. Gene expression profile changes correlated with lymph node metastasis in oral squamous cell carcinoma. Odontology. 2008;96(1):38–43.

    Article  CAS  Google Scholar 

  8. Kato Y, Uzawa K, Saito K, Nakashima D, Kato M, Nimura Y, et al. Gene expression pattern in oral cancer cervical lymph node metastasis. Oncol Rep. 2006;16(5):1009–14.

    CAS  PubMed  Google Scholar 

  9. Kondoh N, Ishikawa T, Ohkura S, Arai M, Hada A, Yamazaki Y, et al. Gene expression signatures that classify the mode of invasion of primary oral squamous cell carcinomas. Mol Carcinog. 2008;47(10):744–56.

    Article  CAS  Google Scholar 

  10. Nagata M, Fujita H, Ida H, Hoshina H, Inoue T, Seki Y, et al. Identification of potential biomarkers of lymph node metastasis in oral squamous cell carcinoma by cDNA microarray analysis. Int J Cancer. 2003;106(5):683–9.

    Article  CAS  Google Scholar 

  11. Nguyen ST, Hasegawa S, Tsuda H, Tomioka H, Ushijima M, Noda M, et al. Identification of a predictive gene expression signature of cervical lymph node metastasis in oral squamous cell carcinoma. Cancer Sci. 2007;98(5):740–6.

    Article  CAS  Google Scholar 

  12. O'Donnell RK, Kupferman M, Wei SJ, Singhal S, Weber R, O'Malley B, et al. Gene expression signature predicts lymphatic metastasis in squamous cell carcinoma of the oral cavity. Oncogene. 2005;24(7):1244–51.

    Article  CAS  Google Scholar 

  13. Roepman P, Wessels LF, Kettelarij N, Kemmeren P, Miles AJ, Lijnzaad P, et al. An expression profile for diagnosis of lymph node metastases from primary head and neck squamous cell carcinomas. Nat Genet. 2005;37(2):182–6.

    Article  CAS  Google Scholar 

  14. Schmalbach CE, Chepeha DB, Giordano TJ, Rubin MA, Teknos TN, Bradford CR, et al. Molecular profiling and the identification of genes associated with metastatic oral cavity/pharynx squamous cell carcinoma. Arch Otolaryngol Head Neck Surg. 2004;130(3):295–302.

    Article  Google Scholar 

  15. Warner GC, Reis PP, Jurisica I, Sultan M, Arora S, Macmillan C, et al. Molecular classification of oral cancer by cDNA microarrays identifies overexpressed genes correlated with nodal metastasis. Int J Cancer. 2004;110(6):857–68.

    Article  CAS  Google Scholar 

  16. Watanabe H, Mogushi K, Miura M, Yoshimura R, Kurabayashi T, Shibuya H, et al. Prediction of lymphatic metastasis based on gene expression profile analysis after brachytherapy for early-stage oral tongue carcinoma. Radiother Oncol. 2008;87(2):237–42.

    Article  CAS  Google Scholar 

  17. Zhou X, Temam S, Oh M, Pungpravat N, Huang BL, Mao L, et al. Global expression-based classification of lymph node metastasis and extracapsular spread of oral tongue squamous cell carcinoma. Neoplasia. 2006;8(11):925–32.

    Article  CAS  Google Scholar 

  18. Fuste NP, Fernandez-Hernandez R, Cemeli T, Mirantes C, Pedraza N, Rafel M, et al. Cytoplasmic cyclin D1 regulates cell invasion and metastasis through the phosphorylation of paxillin. Nat Commun. 2016;7:11581.

    Article  CAS  Google Scholar 

  19. Zhong Z, Yeow WS, Zou C, Wassell R, Wang C, Pestell RG, et al. Cyclin D1/cyclin-dependent kinase 4 interacts with filamin a and affects the migration and invasion potential of breast cancer cells. Cancer Res. 2010;70(5):2105–14.

    Article  CAS  Google Scholar 

  20. Li Z, Jiao X, Wang C, Ju X, Lu Y, Yuan L, et al. Cyclin D1 induction of cellular migration requires p27(KIP1). Cancer Res. 2006;66(20):9986–94.

    Article  CAS  Google Scholar 

  21. Meng H, Tian L, Zhou J, Li Z, Jiao X, Li WW, et al. PACSIN 2 represses cellular migration through direct association with cyclin D1 but not its alternate splice form cyclin D1b. Cell Cycle. 2011;10(1):73–81.

    Article  CAS  Google Scholar 

  22. Chen YW, Huang HS, Shieh YS, Ma KH, Huang SH, Hueng DY, et al. A novel compound NSC745885 exerts an anti-tumor effect on tongue cancer SAS cells in vitro and in vivo. PLoS One. 2014;9(8):e104703.

    Article  Google Scholar 

  23. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 28(5):511–5.

    Article  CAS  Google Scholar 

  24. Tavazoie S, Hughes JD, Campbell MJ, Cho RJ, Church GM. Systematic determination of genetic network architecture. Nat Genet. 1999;22(3):281–5.

    Article  CAS  Google Scholar 

  25. Oka R, Nakashiro K, Goda H, Iwamoto K, Tokuzen N, Hamakawa H. Annexin A8 is a novel molecular marker for detecting lymph node metastasis in oral squamous cell carcinoma. Oncotarget. 2016;7(4):4882–9.

    Article  Google Scholar 

  26. Rickman DS, Millon R, De Reynies A, Thomas E, Wasylyk C, Muller D, et al. Prediction of future metastasis and molecular characterization of head and neck squamous-cell carcinoma based on transcriptome and genome analysis by microarrays. Oncogene. 2008;27(51):6607–22.

    Article  CAS  Google Scholar 

  27. Zenz R, Scheuch H, Martin P, Frank C, Eferl R, Kenner L, et al. c-Jun regulates eyelid closure and skin tumor development through EGFR signaling. Dev Cell. 2003;4(6):879–89.

    Article  CAS  Google Scholar 

  28. Zhao C, Qiao Y, Jonsson P, Wang J, Xu L, Rouhi P, et al. Genome-wide profiling of AP-1-regulated transcription provides insights into the invasiveness of triple-negative breast cancer. Cancer Res. 2014;74(14):3983–94.

    Article  CAS  Google Scholar 

  29. Eferl R, Ricci R, Kenner L, Zenz R, David JP, Rath M, et al. Liver tumor development. C-Jun antagonizes the proapoptotic activity of p53. Cell. 2003;112(2):181–92.

    Article  CAS  Google Scholar 

  30. Xu H, Jin X, Yuan Y, Deng P, Jiang L, Zeng X, et al. Prognostic value from integrative analysis of transcription factors c-Jun and Fra-1 in oral squamous cell carcinoma: a multicenter cohort study. Sci Rep. 2017;7(1):7522.

    Article  Google Scholar 

  31. Shevde LA, Samant RS. Role of osteopontin in the pathophysiology of cancer. Matrix Biol. 2014;37:131–41.

    Article  CAS  Google Scholar 

  32. Subramani VN, Narasimhan M, Thiyagarajan M, Munuswamy BD, Jayamani L. Expression of Osteopontin in Oral squamous cell carcinoma and its surgical margins-an Immunohistochemical study. J Clin Diagn Res. 2015;9(11):ZC66–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  33. Luo SD, Chen YJ, Liu CT, Rau KM, Chen YC, Tsai HT, et al. Osteopontin involves cisplatin resistance and poor prognosis in Oral squamous cell carcinoma. Biomed Res Int. 2015;2015:508587.

    PubMed  PubMed Central  Google Scholar 

  34. Anborgh PH, Mutrie JC, Tuck AB, Chambers AF. Pre- and post-translational regulation of osteopontin in cancer. J Cell Commun Signal. 2011;5(2):111–22.

    Article  Google Scholar 

  35. Shi Z, Mirza M, Wang B, Kennedy MA, Weber GF. Osteopontin-a alters glucose homeostasis in anchorage-independent breast cancer cells. Cancer Lett. 2014;344(1):47–53.

    Article  CAS  Google Scholar 

  36. Body S, Esteve-Arenys A, Miloudi H, Recasens-Zorzo C, Tchakarska G, Moros A, et al. Cytoplasmic cyclin D1 controls the migration and invasiveness of mantle lymphoma cells. Sci Rep. 2017;7(1):13946.

    Article  Google Scholar 

Download references

Acknowledgements

This work was supported by Bureau of Science & Technology and Intellectual Property Nanchong City, China (16YFZJ0026).

Funding

This work was supported by the Bureau of Science & Technology and Intellectual Property Nanchong City, China (16YFZJ0026).

Availability of data and materials

Not applicable.

Author information

Authors and Affiliations

  1. Department of Stomatology, Nanchong Central Hospital, The Second Clinical Medical College of North Sichuan Medical College, Nanchong, Sichuan, 637000, People’s Republic of China

    Xing’an Zhang, Xiaoyao Tan, Xinsheng Han, Huadong Wang & Huawei Ming

  2. Institute of Tissue Engineering and Stem Cells, Nanchong Central Hospital, The Second Clinical Medical College of North Sichuan Medical College, No. 95, People’s south Road, Shunqing District, Nanchong, Sichuan, 637000, People’s Republic of China

    Xing’an Zhang, Qiujiang Li, Kang Liu & Gang Feng

  3. Department of Burn and Plastic Surgery, Nanchong Central Hospital, The Second Clinical Medical College of North Sichuan Medical College, Nanchong, Sichuan, 637000, People’s Republic of China

    Lanfang Zhang

  4. Department of Science and Education, Nanchong Central Hospital, The Second Clinical Medical College of North Sichuan Medical College, Nanchong, Sichuan, 637000, People’s Republic of China

    Ying Lin

Authors
  1. Xing’an Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Lanfang Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Xiaoyao Tan
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Ying Lin
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Xinsheng Han
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Huadong Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Huawei Ming
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Qiujiang Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Kang Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Gang Feng
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

LY and LQ carried out the molecular genetic studies, participated in the sequence alignment. LK and ZX drafted the manuscript. HX and WH carried out the immunoassays. ZL and TX participated in the sequence alignment. ZX and MH participated in the design of the study and performed the statistical analysis. LK and FG conceived of the study and participated in its design and coordination. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Kang Liu or Gang Feng.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional files

Additional file 1:

Table S1. DEGs of GSE70604 Associated Study. (XLSX 106 kb)

Additional file 2:

Table S3. Upregulated DEGs of GSE70604 Associated Study. (XLSX 57 kb)

Additional file 3:

Table S2. Downregulated DEGs of GSE70604 Associated Study. (XLSX 59 kb)

Additional file 4:

Table S4. DEGs of GSE2280 Associated Study. (XLSX 70 kb)

Additional file 5:

Table S5. Downregulated DEGs of GSE2280 Associated Study. (XLSX 42 kb)

Additional file 6:

Table S6. Upregulated DEGs of GSE2280 Associated Study. (XLSX 39 kb)

Additional file 7:

Table S7. Consistent Overlapping DEGs. (XLSX 11 kb)

Additional file 8:

Table S8. Opposite Overlapping DEGs. (XLSX 13 kb)

Additional file 9:

Table S9. Top 10 GO Biological Processes of Consistent Overlapping DEGs. (XLSX 12 kb)

Additional file 10:

Table S10. KEGG Pathways of Consistent Overlapping DEGs. (XLSX 10 kb)

Additional file 11:

Table S11. Interacting Gene Pairs of Consistent Overlapping DEGs. (XLSX 9 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, X., Zhang, L., Tan, X. et al. Systematic analysis of genes involved in oral cancer metastasis to lymph nodes. Cell Mol Biol Lett 23, 53 (2018). https://doi.org/10.1186/s11658-018-0120-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s11658-018-0120-2

Keywords

Cellular & Molecular Biology Letters

ISSN: 1689-1392

Contact us