- Open Access
Aberrant DNA hypermethylation-silenced SOX21-AS1 gene expression and its clinical importance in oral cancer
Clinical Epigeneticsvolume 8, Article number: 129 (2016)
Long noncoding RNAs (lncRNAs) are more than 200 nucleotides in length and lack transcriptional ability. The biological function of lncRNAs in oral squamous cell carcinoma (OSCC) remains unclear. The aim of this study was to identify the dysfunction of lncRNA in OSCC.
We analyzed the transcriptome profiles of human OSCC tissues and paired adjacent normal tissues from two patients through a next-generation sequencing approach. A total of 14 lncRNAs were upregulated (fold change ≥3) and 13 were downregulated (fold change ≤−3) in OSCC tissues compared with the adjacent normal tissues. SOX21-AS1 was subjected to further analysis, revealing that the expression levels of SOX21-AS1 significantly decreased in OSCC compared with the adjacent normal tissue. The promoter activity of SOX21-AS1 was obviously suppressed by in vitro methylation. The DNA methylation status of the SOX21-AS1 promoter was analyzed using combined bisulfite restriction analysis, revealing that the aberrant promoter hypermethylation of SOX21-AS1 was observed frequently in OSCC tissues. The effects of SOX21-AS1 on cell proliferation and invasion were examined through transient transfection. Our data showed that SOX21-AS1 could significantly suppress oral cancer cell growth and invasion. Furthermore, the low expression level of SOX21-AS1 was significantly correlated with an advanced stage (P = 0.047), large tumor size (P = 0.033), and poor disease-specific survival in OSCC patients (P = 0.002).
SOX21-AS1 was identified as susceptible dysfunction correlated with promoter hypermethylation in OSCC. Low SOX21-AS1 expression may be an adverse prognostic biomarker for OSCC.
Oral cancer is one of the most common cancers in developing countries and occurs in the anterior tongue, cheek, floor of the mouth, gingiva, and other parts of the oral cavity . Oral squamous cell carcinoma (OSCC) is the most common oral cancer and 300,400 incident cases and 145,400 deaths have been reported worldwide . Betel quid chewing, tobacco smoking, and alcohol drinking contribute to the increasing incidence and mortality rates of oral cancer .
Long noncoding RNAs (lncRNAs) are over 200 nucleotides in length and lack transcriptional ability . Some lncRNAs, such as transfer RNAs (tRNAs), ribosomal RNAs (rRNAs), and spliceosomal RNAs, are crucial for maintaining normal cellular mechanisms . Except for these functional lncRNAs, most lncRNAs have been frequently considered functionless sequences [6–8]. Recently, lncRNAs have been reported as crucial for regulating cell development, growth, cell cycles, and cancer metastasis in human cancer . Several aberrant lncRNAs were identified from high-throughput profiles and were associated with OSCC progression . Fang et al. reported that urothelial cancer associated 1 (UCA1) is an overexpression in tongue squamous cell carcinomas (TSCCs) and may play an oncogenic role in tumorigenesis . Jia et al.  reported that the loss of miR-26a expression resulted in MEG3 downregulation in TSCC tissues compared with adjacent normal tissues. Low MEG3 expression was highly correlated with poor prognosis of patients. Furthermore, MEG3 overexpression in SCC-15 and CAL27 cells inhibited cell proliferation and cell cycle progression and promoted apoptosis . Oncogenic lncRNA, HOX transcript antisense RNA (HOTAIR), was significantly upregulated in OSCC, and the upregulated expression was associated with poor survival and metastasis. HOTAIR knockdown can suppress oral cancer cell proliferation, colony formation, cell migration, and invasion . However, the role of lncRNA remains largely unknown, and the detailed mechanisms and functions of dysfunctional lncRNAs in OSCC have not yet been fully elucidated. On the basis of the analysis of the expression profile of two pairs of OSCC tissue samples, as well as the experimental studies of two OSCC cell lines and 86 pairs of OSCC tissue samples in this study, we determined that SOX21-AS1 expression was silenced with an aberrant methylation promoter; moreover, the low expression was highly correlated with poor prognosis of patients with OSCC.
Long noncoding RNAs identified in OSCC tissues through next-generation sequencing
Transcriptome profiles were comprehensively analyzed on laser capture microdissected tumors and corresponding adjacent normal tissues from two OSCC patients through a next-generation sequencing (NGS) approach. More than 30 million clean reads were identified, and >84% clean reads were mapped to human reference genes in four individual libraries. After the clean reads were mapped to the genome, over 560 expressed lncRNA genes were detected in four individual libraries (reads per kilobase per million [RPKM] > 1) (Additional file 1: Table S1). Dysregulated lncRNAs were selected from two OSCC patients and were highly expressed (RPKM of normal + RPKM of tumor > 5), and the fold change of RPKM was >3 or <−3 as those of tumor tissues compared with normal tissues. Our data revealed that 109 lncRNAs were dysregulated (fold change ≥3 or <−3) in tumor tissues compared with normal tissues from patient 1; among these lncRNAs, 71 and 38 were upregulated and downregulated, respectively. In addition, 185 lncRNAs were differentially expressed in patient 2; 96 and 89 of these lncRNAs were upregulated and downregulated, respectively. Combining the two profiles revealed that 14 lncRNAs were significantly upregulated (fold change ≥3), and 13 lncRNAs were downregulated (fold change ≤−3) in OSCC tissues compared with adjacent normal tissues (Fig. 1 and Table 1).
DNA methylation-silenced SOX21-AS1 expression in oral cancer cells
Our previous studies have mainly focused on DNA methylation regulating noncoding RNA expression in human cancer [14–17]. Therefore, we sought to determine whether DNA hypermethylation silences lncRNA expression in oral cancer. Because a CpG-rich region is located upstream of SOX21-AS1, we investigated SOX21-AS1 further. SOX21-AS1 is a 2986-bp long noncoding RNA, which shares a bidirectional promoter with SOX21 at human chromosome 13q32.1. SOX21 and SOX21-AS1 share a head-to-head promoter, which is also conserved in the mouse genome. Furthermore, a CpG-rich region is located upstream of SOX21-AS1 and SOX21 (Fig. 2a), implying that their transcriptional activity may be controlled through DNA methylation. Our profiles revealed that SOX21 and SOX21-AS1 are simultaneously downregulated in OSCC tissues compared with paired adjacent normal tissues (Fig. 2b). Mitchell et al. reported that methylation of the sense SOX21 was markedly lower in nonneoplastic colorectal tissues than in colorectal tumors and adenomas . Therefore, we investigated whether SOX21 and SOX21-AS1 were silenced by aberrant DNA methylation in OSCC cell lines. After 5-Aza-dC treatment, SOX21-AS1 and SOX21 expression was increased in the SAS and CAL27 cells (Fig. 2c). The coordinated effects of histone acetylation and DNA methylation were further examined. Our data showed synergistic effects on the increases of SOX21-AS1 and SOX21 expression in the SAS and CAL27 cells with TSA and 5-Aza-dC cotreatment (Fig. 2c). In addition, a DNA methylation assay revealed DNA methylation of the CpG-rich region upstream of SOX21-AS1 in two oral cancer cell lines (Fig. 2d). Consequently, TSA and a DNA demethylating reagent caused marked synergistic activation of SOX21 and SOX21-AS1 genes in oral cancer cells.
To further confirm whether SOX21-AS1 expression was silenced by DNA methylation, we assessed the promoter activity of SOX21-AS1 through an in vitro methylation assay. As shown in Fig. 3a, b, the activity of the PGL-SOX21-AS1 promoter significantly increased (>500-fold) compared with that of the negative control, implying that this DNA fragment exhibited promoter activity. The pGL-SOX21-AS1 promoter luciferase activity was almost entirely repressed after complete methylation of all CpG dinucleotides by M. SssI. However, the promoter constructs of SOX21-AS1 were partially methylated by M. HhaI or M. HpaII, revealing that the promoter activity of SOX21-AS1 was only partially suppressed compared with that of the control groups (Fig. 3c). The luciferase activity of two CpG dinucleotide-less promoter constructs, pGL-SRE and pGL-CRE, were not repressed after complete methylation of all CpG dinucleotides by M. SssI (Additional file 2: Figure S1). Thus, these data revealed that the transcriptional activity of the putative promoter of SOX21-AS1 can be repressed using DNA hypermethylation in oral cancer cells.
DNA hypermethylation results in SOX21-AS1 silencing in OSCC
Our data revealed that SOX21 and SOX21-AS1 expression was significantly reduced in OSCC tissues compared with adjacent normal tissues (SOX21: 62 of 86, P < 0.001; SOX21-AS1: 61 of 86, P < 0.001). Simultaneous downregulation of SOX21 and SOX21-AS1 was quite frequent in the 86 OSCC samples (60.5%; 52 of 86; Fig. 4a, b). Thus, SOX21 expression was significantly correlated with SOX21-AS1 expression in 86 OSCC tissues (R 2 = 0.591, P < 0.001, Fig. 4c). Previous studies have reported that combined bisulfite restriction analysis (COBRA) is well suited for analyzing numerous clinical specimens . Therefore, we further examined the methylation status of the SOX21-AS1 promoter region by using a COBRA assay. Our data revealed frequent DNA methylation in OSCC tissues (CpG1 75.6%, CpG2 60.5%, and CpG3 57%; Fig. 5a). Tumor-specific hypermethylation of CpG region 2 was frequently observed in oral cancer (36 of 86, 41.9%) and resulted in decreased SOX21-AS1 expression in OSCC tissues compared with adjacent normal tissues (25 of 36, 69.4%). In addition, the high methylation status of CpG region 2 also correlated with low SOX21-AS1 expression, but CpG regions 1 and 3 were not thus correlated (Additional file 3: Figure S2). Subsequent bisulfite sequences of three CpG-rich regions in the selected patient (patient 9) were consistent with those of the COBRA data (Fig. 5b). Thus, our data revealed that SOX21 and SOX21-AS1 expression can be silenced using a DNA-methylated promoter in oral cancer.
SOX21-AS1 suppresses oral cancer cell growth and invasion
We analyzed the subcellular localization of endogenous SOX21-AS1, revealing SOX21-AS1 expression accumulating in the nuclei of SAS cells (Additional file 4: Figure S3A). The predominantly nuclear localization of SOX21-AS1 might escape RNA interference-silencing machinery, which is the main action in the cytoplasms of the cells. Therefore, gain-of-function might be a more effective approach to determining the biological function of SOX21-AS1. We constructed a SOX21-AS1 expression vector consisting of two exons with a total length of approximately 2.9 kb. The putative biological function of SOX21-AS1 was examined through transient transfection of pCMV-SOX21-AS1 into SAS cell lines. After transfection, the SOX21-AS1 expression levels were significantly overexpressed (>60-fold increase) compared with the vector control group (Fig. 6a, b). Overexpression of SOX21-AS1 in SAS cells could significantly inhibit cancer cell growth, migration, and invasion (Fig. 6c–h). However, ectopic SOX21-AS1 expression did not influence endogenous SOX21 expression, implying that SOX21-AS executes its biological function through an SOX21-independent pathway (Additional file 4: Figure S3B).
We made two additional deletion constructs, pCMV-SOX21-AS1-E1 and pCMV-SOX21-AS1-E2, for determining which region is responsible for cell growth and motility inhibition. As shown in Fig. 6a, b, both SOX21-AS1-E1 and SOX21-AS1-E2 were significantly overexpressed (>80-fold increase) compared with the vector control group. The effects of SOX21-AS1-E1 and SOX21-AS1-E2 on the proliferation, migration, and invasion of SAS cells were then respectively examined after transient transfection. Our data revealed that ectopic overexpression of SOX21-AS1-E1 revealed no effect on proliferation or motility ability in SAS cells. Ectopic SOX21-AS1-E2 expression in SAS cells significantly suppressed the growth and motility ability of SAS cells (Fig. 6c–h). To understand the putative mechanism of SOX21-AS1 suppressed cell growth and invasion ability, we further examined the epithelial-mesenchymal transition (EMT) genes and growth-related genes. As shown in Additional file 4: Figure S3C, the expression level of E-cadherin was increased, and fibronectin and cyclin D1 were decreased in SAS cells with SOX21-AS1 overexpression. However, the details of the mechanisms of SOX21-AS1 suppressed cell growth and invasive ability remain unclear and require further study.
Correlation between SOX21-AS1 expression and clinicopathological characteristics
The effects of SOX21 and SOX21-AS1 expression on the clinicopathologic outcomes of patients were evaluated, revealing that the low SOX21-AS1 expression was significantly correlated with the pathological stage (P = 0.047) and T classification (P = 0.033, Table 2). However, no correlation was observed between SOX21 expression and clinicopathologic features. We further evaluated the association of SOX21-AS1 expression with survival in OSCC patients. Kaplan–Meier curves and the Cox regression analysis revealed that low SOX21-AS1 expression was significantly associated with a short disease-specific survival (log-rank P = 0.013, Additional file 5: Figure S4 and Table 3). The multivariate Cox regression model revealed that low SOX21-AS1 expression was significantly correlated with shorter disease-specific survival and disease-free survival for patients with OSCC (DSS: adjusted hazard ratio 5.66, 95% CI 1.85–17.30, P = 0.002; DFS: adjusted hazard ratio 2.96, 95% CI 1.14–7.64, P = 0.025, Table 3), after adjustment for cell differentiation, pathological stages, and the expression level of SOX21.
After the entire human genome was completely sequenced in 2003, it was revealed that only approximately 2% of the human genome is used for protein coding . High-throughput sequencing data revealed that 80–90% of the human genomic DNA is actively produced by RNA transcripts [20–22]. Most of these transcripts are nonprotein coding genes, including small RNAs and lncRNAs. In the present study, we comprehensively analyzed OSCC profiles by using an NGS approach, which revealed that the expression of 27 lncRNAs significantly differed between OSCC and the corresponding adjacent normal tissues (Table 1). Among lncRNAs, lnc00152 has been reported to play an oncogenic role in promoting gastric cancer cell growth, migration, and invasion and in suppressing cell apoptosis by sponging miR-18a-5p, miR-195-3p, miR-139-5p, and miR-31-5p expression [23, 24]. In addition, lnc00152 was significantly overexpressed in pancreatic cancer . However, in-depth biological function and the effect of dysfunctional lnc00152 on oral cancer remain unclear and require further investigation.
In this study, our data revealed that SOX21 and SOX21-AS1 share a head-to-head promoter, and DNA methylation of the promoter leads to simultaneous suppression of their expression in OSCC. Our data showed that SOX21-AS1 expression was increased in oral cancer cells treated with 5-Aza-dC. An in vitro methylation assay demonstrated that SOX21-AS1 promoter activity could be suppressed with complete DNA methylation. However, only the methylation status of CpG2 was well correlated with low SOX21-AS1 expression in OSCC, but CpG regions 1 and 3 were not thus correlated (Additional file 3: Figure S2). These conflicting results may be due to two putative reasons: (1) A COBRA assay does not completely represent DNA methylation status; (2) DNA methylation is not the only factor for transcriptional silencing.
A previous study revealed that SOX21 expression was controlled through DNA methylation and was a suitable biomarker for colorectal cancer diagnosis . In human glioma cells, SOX21 overexpression inhibits SOX2 expression and induces apoptosis because of the interaction between SOX21 and SOX2 [26, 27]. Ferletta et al. reported that the knockdown of SOX2 expression can suppress SOX21 expression at the protein and mRNA levels in glioma cells, implying that SOX2 can positively regulate the transcriptional activity of SOX21 . SOX21 in combination with SOX2 exerts a tumor-suppressive effect on carcinogenesis. Although SOX21 and SOX21-AS expression were reduced simultaneously in OSCC, SOX21 expression does not contribute to the clinical manifestations of patients with OSCC (Tables 2 and 3). Conversely, low SOX21-AS1 expression was significantly correlated with poor clinicopathologic features and a short survival time of patients with OSCC (Tables 2 and 3).
Although SOX21-AS1 can suppress oral cancer cell growth and invasion, the details of the underlying mechanism remain unknown. LncRNAs can modulate cellular function through different mechanisms such as RNA degradation, epigenetic modification, transcriptional regulation, and an miRNA decoy . We analyzed the subcellular fractionation localization of SOX21-AS1 in SAS cells (Additional file 4: Figure S3). Furthermore, a considerable increase in SOX21-AS1 expression was observed in the nuclei compared with the cytoplasms. According to its location and tumor-suppressive role, SOX21-AS1 is probably an enhancer that reinforces other tumor suppressor mRNAs or acts as a decoy removing certain transcription factors from binding to the DNA to prevent oncogene transcription. In conclusion, low SOX21-AS1 expression correlated with its sense gene SOX21 and can serve as an independent biomarker for poor prognosis of OSCC. SOX21-AS1 and SOX21 expression in tumors may be downregulated through DNA hypermethylation. However, the details of the biological functions of SOX21-AS1 must be further evaluated in OSCC.
Clinical tissue samples were collected from OSCC patients who provided signed informed consent and accepted a surgical operation at the Department of Dentistry and Department of Otorhinolaryngology, Kaohsiung Veterans General Hospital. This study protocol was approved by the Institutional Review Board of Kaohsiung Veterans General Hospital (Kaohsiung, Taiwan; IRB number VGHKS14-CT6-18). The total RNA and DNA of the tissues were extracted using the TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the instruction manual.
Next-generation sequencing and analysis
In this study, the transcriptome profiles of paired tumors and adjacent normal tissue samples from two OSCC patients were analyzed using an NGS approach; the detailed clinical characteristics of the samples are shown in Additional file 6: Table S2. High-purity normal and tumor tissues were individually separated through laser capture microdissection and were transferred to an RNA extraction kit (Qiagen, Hilden, Germany) for RNA extraction. All procedures of the RNA transcriptome were performed according to the manufacturer protocol from Illumina. Library construction of all samples was performed using TruSeq RNA Sample Prep Kits v2 for 50 single-end base pair sequencing on the Solexa platform. Finally, raw sequences were obtained from the Illumina GA Pipeline software CASAVA version 1.8 and expected to generate 30 million reads per sample. After low-quality data were filtered, qualified reads were analyzed using TopHat/Cufflinks for estimating gene expression, which was calculated as RPKM. For differential expression analyses, CummeRbund was used for statistically analyzing the gene expression profiles. The reference genome and gene annotations were retrieved from the Ensembl database.
Cell lines and demethylation treatment
Two human OSCC cell lines, SAS (human tongue SCC) and CAL27 (human tongue SCC) cell lines, were provided by Dr. Michael Hsiao (Genomics Research Center, Academia Sinica, Taipei, Taiwan). The oral cancer cell lines, SAS and CAL27, were cultured in Dulbecco's modified Eagle’s medium supplemented with 10% fetal bovine serum (Invitrogen, Carlsbad, CA, USA) and 1% penicillin/streptomycin (Sigma-Aldrich, St. Louis, MO, USA). To examine the role of methylation in the regulation of lncRNA expression, SAS and CAL27 cells were cultured in the presence or absence of 10 μM 5-aza-2′-deoxycytidine (5-Aza-dC) for 4 days and treated with 0.25 μM TSA on the third day.
Reverse transcriptase polymerase chain reaction of lncRNAs
For RT-PCR, 2 μg of total RNA was reverse transcribed using random primers and SuperScript III reverse transcriptase according to manufacturer instructions (Invitrogen, Carlsbad, CA, USA). Gene expression was detected using an SYBR Green I assay (Applied Biosystems, Foster City, CA, USA), and lncRNA expression was normalized to that of GAPDH (△Ct = target lncRNA Ct-GAPDH Ct). The individual primers used in the present study were as follows: SOX21-F: GCACAACTCGGAGATCAGCA and SOX21-R: CCGGGAAGGCGAACTTGTC; SOX21-AS1-exon1-F: CCGATGGGAAACCCCCAATC and SOX21-AS1-exon1-R: AACGCTTGCTCAAGCCTCAT; SOX21-AS1-exon2-F: TCACTTACATGCGCTGCTGA and SOX21-AS1-exon2-R: GCCGCAGCATACCAAAAAGT; GAPDH-F: TGCACCACCAACTGCTTAGC and GAPDH-R: GGCATGGACTGTGGTCATGAG.
SOX21-AS1 promoter and expression construction
The SOX21-AS1 promoter was amplified using the promoter-specific primer pairs (SOX21AS1-promoter-F: ATGAAGCTTCCATGAAGGCGTTCATGGGCCG and SOX21AS1-promoter-R: ATGAAGCTTAGAGGAAGACTCGAGAGGCAGGT). PCR products were digested with the restriction enzyme HindIII and cloned into the pGL4.21 luciferase expression vector (Promega, Madison, WI, USA). Full-length exons 1 and 2 of the SOX21-AS1 PCR products were amplified using individual primer pairs (exon1-F: ATGGAATTCTCTTCTTGGCTCCGGGCAGGGTG and exon1-R: ATGAAGCTTCTGAGCCGGTGCAGAGGGCG; exon2-F: ATGGAATTCGTTTAGGCGAGTGGAGAGTCCG, and exon2-R: ATGCTCGAGAATCTTTAGGACAAAACTGAGC). The PCR products were digested with the restriction enzyme EcoRI and cloned into the pCMV-Tag 2A vector (Stratagene, La Jolla, CA, USA).
In vitro methylation
We sequentially methylated pGL4-SV40, pGL4-SRE promoter, pGL4-CRE promoter, and pGL4-SOX21-AS1-P in vitro by using M. SssI, M. HpaII, and M. HhaI methyltransferase enzymes (Invitrogen, Grand Island, NY, USA). A luciferase assay was performed in SAS cells by using a Dual-Glo luciferase reporter assay system kit (Promega, Madison, WI, USA) 24 h after transfection.
Combined bisulfite restriction and sequencing analyses
DNA was subjected to bisulfite conversion by using the EZ DNA Methylation-Gold Kit (Zymo Research Corporation, Orange, CA, USA). The bisulfite-converted genomic DNA was used for methylation analysis of the promoter with the specific methylation primers. The methylation status of the genomic DNA of individual samples was examined using TaqI or BstUI digestion (New England Biolabs, MA, USA). Furthermore, the PCR products were cloned into the pJET1.2 vector (Promega, Madison, WI, USA), and several clones were used for sequencing. The individual primers used in the present study were as follows: CpG1-F: ATGAAATTTTTAATAAAATTGGAAAGGT and CpG1-R: CCAAATAAAA ACAAA AAAACCAAAC; CpG2-F: GGTTGTTTTTGGGATATTTTAATTTT and CpG2-R: CTAAAAACCCCCTTTAACACTTAAC; CpG3-F: GGAGGAGGTGGAGTTTAGGATT and CpG3-R: AAAACCACAACCAAAACAACTACA.
SAS cells were transfected with pCMV-SOX21-AS1, pCMV-SOX21-AS1-exon1, pCMV-SOX21-AS1-exon2, and empty vectors by using Lipofectamine® 2000 Reagent (Invitrogen, Grand Island, NY, USA). For the clonogenic assay, 2000 transfecting cells were plated onto 6-well plates for 2 weeks until substantial-sized colonies were formed. Subsequently, the colonies were analyzed using crystal violet staining.
Cell migration and invasion assays
The migration and invasion abilities of cells were assessed in vitro by using a transwell assay. In brief, cells were resuspended at a density of 4.5 × 105 in 2% fetal bovine serum and then added to the upper chamber of the transwells (Falcon, Corning Incorporated, USA) without Matrigel (BD Biosciences, MA, USA) for the migration assay or with a Matrigel coating for the invasion assay. Chambers were incubated in a CO2 incubator at 37 °C for 36–48 h; the remaining cells in the upper chamber were removed using cotton swabs, and the cells on the undersurface of the transwells were fixed in 10% formaldehyde solution. Cells were stained with crystal violet solution, and the numbers of cancer cells in three fields were counted under a phase-contrast microscope.
The chi-square test, Fisher exact test, Student t test, ANOVA, Mann–Whitney U test, and Kruskal–Wallis one-way ANOVA were used to evaluate the correlation of each lncRNA expression with different oral tissues or clinicopathological parameters. For clinicopathologic outcome and survival analysis, the RNA expression levels (-delta Ct) of SOX21-AS1 were dichotomized as low expression and high expression with the cutoff value (−9.8647) based on a receiver operating characteristic (ROC) curve analysis. The clinicopathologic outcome is usually defined as the time of initial diagnosis or surgery. Disease-specific survival was measured from the time of the initial resection of the primary tumor to the date of cancer-specific death or the final follow-up. Disease-free survival was calculated from the date of the initial resection of the primary tumor to the date of recurrence or the final follow-up. In this study, the median time of follow-up was 25.68 (range, 1.00–68.27) months. The cumulative survival curves were estimated using the Kaplan–Meier method, and the survival curves were compared using the log-rank test. A Cox proportional hazards model was used for determining independent predictors of survival by using factors significant in univariate analysis as covariates; P < 0.05 was considered statistically significant.
Ram H, Sarkar J, Kumar H, Konwar R, Bhatt ML, Mohammad S. Oral cancer: risk factors and molecular pathogenesis. J Maxillofac Oral Surg. 2011;10(2):132–7.
Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin. 2015;65(2):87–108.
Petti S, Masood M, Scully C. The magnitude of tobacco smoking-betel quid chewing-alcohol drinking interaction effect on oral cancer in South-East Asia. A meta-analysis of observational studies. PLoS One. 2013;8(11):e78999.
Lipovich L, Johnson R, Lin CY. MacroRNA underdogs in a microRNA world: evolutionary, regulatory, and biomedical significance of mammalian long non-protein-coding RNA. Biochim Biophys Acta. 2010;1799(9):597–615.
Wilusz JE, Sunwoo H, Spector DL. Long noncoding RNAs: functional surprises from the RNA world. Genes Dev. 2009;23(13):1494–504.
Wang J, Zhang J, Zheng H, Li J, Liu D, Li H, Samudrala R, Yu J, Wong GK. Mouse transcriptome: neutral evolution of ‘non-coding’ complementary DNAs. Nature. 2004;431(7010):1. following 757; discussion following 757.
Struhl K. Transcriptional noise and the fidelity of initiation by RNA polymerase II. Nat Struct Mol Biol. 2007;14(2):103–5.
Ebisuya M, Yamamoto T, Nakajima M, Nishida E. Ripples from neighbouring transcription. Nat Cell Biol. 2008;10(9):1106–13.
Fatima R, Akhade VS, Pal D, Rao SM. Long noncoding RNAs in development and cancer: potential biomarkers and therapeutic targets. Mol Cell Ther. 2015;3:5.
Zhang S, Tian L, Ma P, Sun Q, Zhang K, Guanchao W, Liu H, Xu B. Potential role of differentially expressed lncRNAs in the pathogenesis of oral squamous cell carcinoma. Arch Oral Biol. 2015;60(10):1581–7.
Fang Z, Wu L, Wang L, Yang Y, Meng Y, Yang H. Increased expression of the long non-coding RNA UCA1 in tongue squamous cell carcinomas: a possible correlation with cancer metastasis. Oral Surg Oral Med Oral Pathol Oral Radiol. 2014;117(1):89–95.
Jia LF, Wei SB, Gan YH, Guo Y, Gong K, Mitchelson K, Cheng J, Yu GY. Expression, regulation and roles of MiR-26a and MEG3 in tongue squamous cell carcinoma. Int J Cancer. 2014;135(10):2282–93.
Wu Y, Zhang L, Wang Y, Li H, Ren X, Wei F, Yu W, Liu T, Wang X, Zhou X, et al. Long non-coding RNA HOTAIR promotes tumor cell invasion and metastasis by recruiting EZH2 and repressing E-cadherin in oral squamous cell carcinoma. Int J Oncol. 2015;46(6):2586–94.
Chen WS, Leung CM, Pan HW, Hu LY, Li SC, Ho MR, Tsai KW. Silencing of miR-1-1 and miR-133a-2 cluster expression by DNA hypermethylation in colorectal cancer. Oncol Rep. 2012;28(3):1069–76.
Hou YY, You JJ, Yang CM, Pan HW, Chen HC, Lee JH, Lin YS, Liou HH, Liu PF, Chi CC, et al. Aberrant DNA hypomethylation of miR-196b contributes to migration and invasion of oral cancer. Oncol Lett. 2016;11(6):4013–21.
Tsai KW, Liao YL, Wu CW, Hu LY, Li SC, Chan WC, Ho MR, Lai CH, Kao HW, Fang WL, et al. Aberrant hypermethylation of miR-9 genes in gastric cancer. Epigenetics. 2011;6(10):1189–97.
Tsai KW, Wu CW, Hu LY, Li SC, Liao YL, Lai CH, Kao HW, Fang WL, Huang KH, Chan WC, et al. Epigenetic regulation of miR-34b and miR-129 expression in gastric cancer. Int J Cancer. 2011;129(11):2600–10.
Mitchell SM, Ross JP, Drew HR, Ho T, Brown GS, Saunders NF, Duesing KR, Buckley MJ, Dunne R, Beetson I, et al. A panel of genes methylated with high frequency in colorectal cancer. BMC Cancer. 2014;14:54.
Xiong Z, Laird PW. COBRA: a sensitive and quantitative DNA methylation assay. Nucleic Acids Res. 1997;25(12):2532–4.
Consortium EP, Birney E, Stamatoyannopoulos JA, Dutta A, Guigo R, Gingeras TR, Margulies EH, Weng Z, Snyder M, Dermitzakis ET, et al. Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project. Nature. 2007;447(7146):799–816.
Consortium EP. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489(7414):57–74.
Hangauer MJ, Vaughn IW, McManus MT. Pervasive transcription of the human genome produces thousands of previously unidentified long intergenic noncoding RNAs. PLoS Genet. 2013;9(6):e1003569.
Xia T, Liao Q, Jiang X, Shao Y, Xiao B, Xi Y, Guo J. Long noncoding RNA associated-competing endogenous RNAs in gastric cancer. Sci Rep. 2014;4:6088.
Cao WJ, Wu HL, He BS, Zhang YS, Zhang ZY. Analysis of long non-coding RNA expression profiles in gastric cancer. World J Gastroenterol. 2013;19(23):3658–64.
Muller S, Raulefs S, Bruns P, Afonso-Grunz F, Plotner A, Thermann R, Jager C, Schlitter AM, Kong B, Regel I, et al. Next-generation sequencing reveals novel differentially regulated mRNAs, lncRNAs, miRNAs, sdRNAs and a piRNA in pancreatic cancer. Mol Cancer. 2015;14:94.
Caglayan D, Lundin E, Kastemar M, Westermark B, Ferletta M. Sox21 inhibits glioma progression in vivo by forming complexes with Sox2 and stimulating aberrant differentiation. Int J cancer. 2013;133(6):1345–56.
Ferletta M, Caglayan D, Mokvist L, Jiang Y, Kastemar M, Uhrbom L, Westermark B. Forced expression of Sox21 inhibits Sox2 and induces apoptosis in human glioma cells. Int J Cancer. 2011;129(1):45–60.
Rinn JL, Chang HY. Genome regulation by long noncoding RNAs. Annu Rev Biochem. 2012;81:145–66.
The authors would like to thank Department of Medical Education and Research, Kaohsiung Veterans General Hospital, for providing the laser capture microdissection system and Genomics and Proteomics Core Libratory, Department of Medical Research, Kaohsiung Chang Gung Memorial Hospital, for the help with NGS data analysis.
This work was supported by grants from Kaohsiung Veterans General Hospital (VGHKS 104-051, VGHKS 104-092, VGHKS 104-103, and VGHKS104-G01-2).
Availability of data and materials
Next-generation sequencing data, methylation data, and analytical conditions can be provided upon request.
CMY and THW executed this study and wrote the draft of this manuscript. SCL performed analysis of NGS data. HCC and MCL performed clinical sample collection. HHL and YKT performed the statistical analysis. PFL and YLS helped to prepare manuscript. LPG and KWT supervised the study and edited the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
This study protocol was approved by the institutional review board of Kaohsiung Veterans General Hospital (Kaohsiung, Taiwan; IRB number VGHKS14-CT6-18).
The categories of sequence reads in the four libraries. (DOC 28 kb)
Promoter activity of CpG dinucleotideless was examined using in vitro methylation assay. Luciferase activity of CpG dinucleotide-less promoters was analyzed through in vitro methylation. (A, B) Schema of the luciferase constructs containing two CpG dinucleotide-less promoters of pGL-serum-response element and pGL-cAMP-response element (upper panels). Promoter constructs were methylated in vitro by using M. SssI methylase enzymes, and luciferase activity was examined using the Dual-Glo luciferase reporter assay system kit. (PPT 174 kb)
The correlation between od SOX21-AS1 expression and DNA methylation status. Correlation between SOX21-AS1 expression and DNA methylation status was examined in OSCC from 86 patients. (PPT 134 kb)
The cellular localization of SOX21-AS1 and putative mechanism were examined in SAS cell. (A) Cell-fractioned RNA was isolated from the nuclei and cytoplasms of SAS cells. Expression levels of SOX21-AS1 were detected using RT-qPCR. (B) Expression levels of SOX21 were examined using real-time PCR after SOX21-AS1 was transfected into SAS cells. (C) Expression levels of EMT- and cell cycle-relative genes were examined using a Western blot approach in SAS cells with SOX21-AS1 overexpression. (PPT 2411 kb)
Analysis of the prognostic significance of SOX21-AS1 expression in oral cancer. (A, B) DSS and DSF were compared according to SOX21-AS1 expression levels in oral cancer tissues. (PPT 239 kb)
The detailed clinical characteristics of two OSCC patients. (DOC 29 kb)