A novel 3-miRNA network regulates tumour progression in oral squamous cell carcinoma
Biomarker Research volume 11, Article number: 64 (2023)
Late diagnosis is one of the major confounders in oral squamous cell carcinoma (OSCC). Despite recent advances in molecular diagnostics, no disease-specific biomarkers are clinically available for early risk prediction of OSCC. Therefore, it is important to identify robust biomarkers that are detectable using non-invasive liquid biopsy techniques to facilitate the early diagnosis of oral cancer. This study identified potential salivary exosome-derived miRNA biomarkers and crucial miRNA-mRNA networks/underlying mechanisms responsible for OSCC progression.
Small RNASeq (n = 23) was performed in order to identify potential miRNA biomarkers in both tissue and salivary exosomes derived from OSCC patients. Further, integrated analysis of The Cancer Genome Atlas (TCGA) datasets (n = 114), qPCR validation on larger patient cohorts (n = 70) and statistical analysis with various clinicopathological parameters was conducted to assess the effectiveness of the identified miRNA signature. miRNA-mRNA networks and pathway analysis was conducted by integrating the transcriptome sequencing and TCGA data. The OECM-1 cell line was transfected with the identified miRNA signature in order to observe its effect on various functional mechanisms such as cell proliferation, cell cycle, apoptosis, invasive as well as migratory potential and the downstream signaling pathways regulated by these miRNA-mRNA networks.
Small RNASeq and TCGA data identified 12 differentially expressed miRNAs in OSCC patients compared to controls. On validating these findings in a larger cohort of patients, miR-140-5p, miR-143-5p, and miR-145-5p were found to be significantly downregulated. This 3-miRNA signature demonstrated higher efficacy in predicting disease progression and clinically correlated with poor prognosis (p < 0.05). Transcriptome, TCGA, and miRNA-mRNA network analysis identified HIF1a, CDH1, CD44, EGFR, and CCND1 as hub genes regulated by the miRNA signature. Further, transfection-mediated upregulation of the 3-miRNA signature significantly decreased cell proliferation, induced apoptosis, resulted in G2/M phase cell cycle arrest and reduced the invasive and migratory potential by reversing the EMT process in the OECM-1 cell line.
Thus, this study identifies a 3-miRNA signature that can be utilized as a potential biomarker for predicting disease progression of OSCC and uncovers the underlying mechanisms responsible for converting a normal epithelial cell into a malignant phenotype.
Oral squamous cell carcinoma (OSCC) is one of the major causes of morbidity and mortality worldwide . For over four decades, OSCC has been associated with an extremely low five-year survival rate, mainly due to late-stage detection. Despite recent advances in molecular diagnostics, no disease-specific biomarkers are clinically available for early risk prediction of OSCC. Conventional tissue biopsies fail to access deeper tumour sites and are unable to detect tumour heterogeneity [2, 3]. Therefore, it is important to identify non-invasive biomarkers to facilitate the early diagnosis of OSCC.
Tumor-derived exosomes (40-150 nm in diameter) are useful sources of liquid biopsies as they are secreted in various bodily fluids including blood, saliva, urine, and cerebrospinal fluid. They mediate intercellular communication through transmission of active biomolecules (mRNA, microRNA (miRNA), DNA, proteins, and lipids) into the tumour microenvironment. These exosomes are an attractive target for early diagnosis and risk prediction of cancer as they represent the molecular and mutational landscape of the tumour bulk . Furthermore, the close proximity of saliva to the oral cavity and ease of collection makes it an ideal source to derive exosomes shed from oral cavity tumors for a liquid biopsy approach . Exosomal miRNAs have demonstrated immense clinical potential in early risk prediction, assessment of disease progression, and real-time therapeutic monitoring in various malignancies . Since miRNAs are highly tissue-specific, and negatively regulate 60% of protein-coding genes in humans, it is necessary to identify miRNA-based biomarkers and elucidate miRNA-mRNA networks that lead to the deregulation of various biological mechanisms and signalling pathways [7,8,9]. So far, miRNA-mRNA network studies have failed to conduct a comprehensive analysis of small RNA and transcriptome profiles using matched samples. Moreover, these studies involved smaller patient cohorts and lacked functional validation using in-vitro or in-vivo models.
Thus, in this study, we conducted a comparative analysis of miRNA and transcriptome profiles of matched tumour tissues and salivary exosomes derived from OSCC patients. Further, we identified regulatory networks and underlying mechanisms responsible for tumour progression in OSCC. To the best of our knowledge, this is the first study to explore the potential clinical utility of saliva and tumour-derived miRNA-mRNA regulatory networks in OSCC.
Resected tumour tissue and unstimulated saliva samples were collected from patients diagnosed with OSCC (excluding tongue, larynx, pharynx, and hypopharynx) at the HCG Cancer Centre, Ahmedabad. Patients were authorized for primary surgical resection by the institutional team after their malignancies were confirmed by histopathological examination. HPV, HIV, HCV and HbsAg-positive malignancies, pediatric cases and patients initially treated with other adjuvant therapeutic modalities were excluded from the study. All patients provided written informed consent for their participation in the study. This study was reviewed and approved by the HCG Cancer Centre’s Ethics Committee (ECR/92/Inst/GJ/2013/RR-16) for human subject research and the use of human tissues. This study also complied with the guidelines set forth by the Declaration of Helsinki (2008).
Tissue and saliva collection
Samples were collected during R0 resections of early and late-stage OSCC. Resected specimens were evaluated histopathologically and residual tumour samples were snap-frozen and stored at -80 °C for further processing. For normal controls, brush biopsy samples were taken from healthy individuals (age and sex-matched with patients) with no clinically detectable oral lesions.
Whole saliva was collected in sterile tubes from healthy individuals and OSCC patients prior to surgical resection, using a protocol described in a previous study . The saliva samples were kept at -80 °C until further use. A total of 23 patient-derived samples and 21 samples from age- and sex-matched healthy volunteers were sequenced. Sequencing results were validated by qRT-PCR using a validation cohort of patients. The details of the samples are given in Table 1.
Saliva from OSCC patients and healthy controls was centrifuged at 2000 × g for 10 min at 37 °C to remove cells and debris. The supernatant was transferred to a sterile tube, and exosomes were precipitated using a commercially available kit as per the manufacturer’s protocol (Invitrogen™ Total Exosome Isolation Reagent (from other body fluids), Thermo Fisher Scientific, Waltham, MA, USA). The exosome pellet was either re-suspended in 50µL of phosphate-buffered saline (PBS) for characterization or in 750µL of TRIzol LS reagent (Thermo Fisher Scientific) for storage at -80 °C for subsequent RNA extraction.
The concentration (particles/mL) and size of exosomes were measured using NanoSight LM10 Nanoparticle Tracking Analysis (NTA) software (Malvern Instruments Ltd., Malvern, UK). The exosome size was also determined by transmission electron microscopy (TEM) using 300 mesh carbon formvar grids. 5µL of 2% glutaraldehyde was added to the grids which were then wick-dried, rinsed, and stained with 2% uranyl acetate. The grids were examined using JEM1400 Plus transmission electron microscope (Jeol, Tokyo, Japan) at 100 kV. Images were obtained via a chamber scope film camera. The isolated exosomes were stained with anti-human CD9 allophycocyanin (CD9-APC), anti-human CD63 Alexa Fluor 488 (CD63-Alexa 488), and anti-human CD81 phycoerythrin (CD81-PE antibodies (Thermo Fisher Scientific). The samples were then analysed using the FACS Calibur™ Flow Cytometer (Becton Dickinson Biosciences, San Jose, CA).
Total RNA was extracted from tissue and salivary exosomes of OSCC patients using the described method . The RNA pellet was resuspended in 32µL of nuclease-free water. RNA yield and purity were quantified using Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) and Qubit 4 Fluorometer (Thermo Fisher Scientific).
miRNA and transcriptome sequencing
Small RNA and transcriptome sequencing of tumour tissues and salivary exosomes derived from OSCC patients was conducted at Genotypic Technology Pvt. Ltd., Bengaluru, India. The libraries were prepared using QIAseq® miRNA library kit protocol (QIAGEN, Frederick, MD, USA) as per the manufacturer’s instructions. Quality control and read mapping were performed as mentioned in our previous study . miRNAs were considered expressed when detected with a miRDeep2 score ≥ 5, reads > 15 nucleotides in length, quality score > 30 reads/sample, RPM normalized expression value ≥ the first quantile of expression distribution, false discovery rate (FDR) of ≤ 0.05, log2 fold-change (log2FC) of > 2, and statistically significant abundance in each sample (adj p-value < 0.05). This script produced a raw count file that was used as input for the DESeq2 R/Bioconductor package for further analysis (Additional File 1: Fig. S1).
For transcriptome profiling, library preparation was performed using NEBNext® Ultra™ II directional RNA library prep kit and sequencing was performed on Illumina HiSeq 2500 platform using 150 × 2 chemistry. Raw sequencing reads were subjected to quality assessment using FastQC, followed by cleaning with Cutadapt [11, 12]. The refined reads were compared with a reference genome (GRCh, version 38) using Subread aligner [10, 13]. Differential expression analysis was conducted using DESeq2 . Transcripts with read counts ≥ 30 were included. miRNAs and/or genes were considered differentially expressed based on the criteria: log2FC > 2 and the adjusted p-value < 0.01.
Reverse transcription and quantitative real-time polymerase chain reaction (qRT-PCR)
For small RNA expression analysis of miR-140, miR-143, miR-145, miR-30a, miR-21, miR-423, and miR-let7a, total RNA was reverse transcribed using the miScript-II RT kit (QIAGEN) according to the manufacturer’s instructions. For quantifying the expression of genes (E-cadherin and N-cadherin), reverse transcription of RNA extracted from OECM-1 cells transfected with miRNA mimics and SC control was carried out using the Quantitect Reverse Transcription Kit (QIAGEN) according to the manufacturer’s protocol. The expression of identified miRNAs and genes were validated using the miScript SYBR® Green PCR kit (QIAGEN) on the QuantStudio 5 Real-Time PCR system (Applied Biosystems, Foster City, CA, USA). U6 snRNA and β-actin were used as endogenous controls for miRNA and mRNA PCRs respectively. Relative expression of miRNAs and genes was calculated using the 2−ΔΔCt method. The sequences of the primers used are specified in Additional File 1: Table S1.
miRNA target prediction
Candidate miRNA-mRNA target relationships were predicted by at least one or more of the following target prediction algorithms extracted from miRDB, TargetScan 5.1, and miRWalk [15,16,17]. Genes were considered differentially expressed if log2FC > 2, and adj p-value was ≤ 0.05.
Interaction networks and pathway analysis
Protein-protein interaction and regulatory networks of the differentially expressed genes were generated using Cytoscape (v3.6.1) and Ingenuity Pathways Analysis software (IPA) (QIAGEN Inc., https://digitalinsights.qiagen.com/IPA) respectively [18, 19]. The statistically significant canonical pathway analysis was conducted using IPA. Z-score > ± 2 was defined as the threshold of significant activation and/or inhibition.
Cell Culture and transfection
The human oral squamous cell carcinoma cell line, OECM-1, was obtained from Sigma Aldrich (Merck, Darmstadt, Germany). Cells were maintained in Roswell Park Memorial Institute (RPMI) 1640 Medium (Thermo Fisher Scientific) containing 10% Fetal Bovine Serum (FBS) at 37 °C in a humid atmosphere with 5% CO2. Transfection was performed using miR-143-5p, miR-145-5p, and miR-140-5p mimics and corresponding scrambled control (SC) (Thermo Fisher Scientific) using Lipofectamine 2000 reagent (Invitrogen, Waltham, MA, USA) according to the manufacturer’s protocol. The total concentration of miR-143/miR-145/miR-140 mimics combination was 100nM and the transfection efficiency was monitored using qRT-PCR after 48 h.
Cell proliferation assay
OECM-1 cells were seeded in 96-well plates at a density of 3 × 103 cells/well. Post transfection for 48 h, cell viability was assessed using MTT (3-(4,5-Dimethylthiazol-2-yl)-2,5-Diphenyltetrazolium Bromide) (Sigma-Aldrich). 10 µl MTT (5 mg/ml) was added into each well followed by 100 µl dimethyl sulfoxide (DMSO). Absorbance was measured at 570 nm using a microplate reader (Synergy™ H1; Agilent BioTek Instruments, Winooski, VT, USA).
Cell cycle analysis
OECM-1 cells were seeded in a 6-well plate at the density of 0.5 × 106 cells/well. 48 h post-transfection, cells were harvested and fixed in 1ml of 70% ethanol. For flow cytometric analysis, the fixed cells were resuspended in 1ml of 50 µg/ml stock solution of propidium iodide (PI) (Becton Dickinson Biosciences). Further, the samples were incubated for 30 min at 4ºC and analyzed on a FACSCalibur™ Flow Cytometer (Becton Dickinson Biosciences).
Post transfection, the percentage of apoptotic cells was measured using Annexin V-fluorescein isothiocyanate (FITC)/PI kit using FACSCalibur™ Flow Cytometer (Becton Dickinson Biosciences) in accordance with the manufacturer’s instructions (Biolegend, London, UK).
Nuclear fragmentation using propidium iodide staining
After 48 h of transfection, the monolayered OECM-1 cells were stained with 10 µg/ml PI (BD Biosciences) (in PBS) for 5 min followed by visualization under the ZOE™ Fluorescent Cell Imager at 590 nm (Bio-Rad, Hercules, CA, USA).
Evaluation of mitochondrial membrane potential (MMP)
Post transfection, the cells were stained with JC-1 (2 mg/ml) and incubated for 15 min in dark. Changes in MMP were visualized under a fluorescent microscope (Axioskope II, Carl Zeiss, Oberkochen, Germany).
Wound healing assay
OECM-1 cells were seeded in a 24-well plate at a density of 0.5 × 105 cells/well. Post transfection, a scratch was made using a 200 µl pipette tip and wounds were photographed at time intervals of 0 h, 24 and 48 h. The migratory distance of the cells post-induction of the wound was measured using ImageJ software (NIH; http://imagej.nih.gov/ij/).
Transwell invasion assay
In order to assess the invasive capability, OECM-1 cells were seeded in 6-well plates with a density of 0.5 × 106 cells/well. After transfection, the cells were subjected to a transwell invasion using the QCM Collagen Cell Invasion Assay kit (Millipore Sigma, Burlington, MA, USA) as per the manufacturer’s instructions. After 48 h of incubation, the cells in the bottom chamber were detached, lysed, stained, and quantified using the 480/520nm filter set in Synergy™ H1 microplate reader (Agilent BioTek Instruments).
For imaging, the cells were stained with 10% Giemsa and fixed with formaldehyde after 48 h of incubation in the transwell chamber. The cells in the top chamber were removed and those on the bottom of the transwell were visualised at 20x magnification under the ZOE™ Fluorescent Cell Imager (Bio-Rad).
The results of the student’s t-test are represented as mean ± standard deviation (SD) and.
visualised using GraphPad Prism 6.01 (GraphPad Software,
www.graphpad.com). Receiver Operating Curves (ROC) were constructed by plotting relative miRNA expression levels using GraphPad Prism software. Principal Component Analysis (PCA) was performed to investigate relationships between samples, and the silhouette score was used to estimate the optimal number of clusters using the prcomp function in the R statistical package. The association of the miRNA signature with the overall patient survival was evaluated by Kaplan–Meier Survival Curve analysis with Log-rank statistic using TCGA data (n = 114). Survival analysis was performed using SPSS software (v19.0) (SPSS, Inc.). All experiments were performed in triplicates. Statistical significance was set at p-value < 0.05.
Identification and characterization of salivary exosomes derived from OSCC patients
The isolated salivary exosomes from OSCC patients were characterized using NTA, TEM, and tetraspanin expression. NTA results showed a single peak for the concentration of exosomes in the size range of 30-40 nm. Salivary exosomes derived from OSCC patients showed a 10-fold higher concentration (7.31 × 108 particles/mL) compared to controls (6.67 × 107 particles/mL) (Additional File 1: Fig. S2A). TEM results represent sphere-shaped vesicles with a mean radius of 40-50 nm, consistent with the NTA profiles (Additional File 1: Fig. S2B). Flow cytometry analysis indicated the presence of the tetraspanins CD9, CD63, and CD81 (86–88% exosomes expressed these markers) in salivary exosomes from both patients and healthy counterparts (Additional File 1: Fig. S2C-D), thus confirming that the isolated vesicles comprise of a pure exosomal population.
Differentially expressed miRNAs enriched in both salivary exosomes and primary tumours from OSCC patients identified by miRNA profiling
Small RNA sequencing was performed to identify differentially expressed miRNAs in tumours and salivary exosomes derived from OSCC patients as compared to age- and sex-matched controls. An average of 25.4 million reads per sample were obtained for all tissue and salivary exosome samples. After pre-processing, adapter sequence trimming and filtering low-quality reads, the mean Q20 values for tissue- and salivary exosome-derived libraries were 97.6% and 93.40%, respectively. Amongst these, an average of 68% (range 61–75%) of the reads could be mapped to hg38 with high confidence.
188 miRNAs were differentially expressed in tumour tissues, 181 miRNAs in leukoplakia samples, 316 miRNAs in salivary tumour exosomes and 69 miRNAs in salivary leukoplakia exosomes compared to their normal counterparts. Amongst these, 47 miRNAs were found to be differentially expressed across all sample types and stages, including leukoplakia lesions (Fig. 1A). A comparison of our results with the data from The Cancer Genome Atlas (TCGA) database (n = 114) showed that 12 of these miRNAs were commonly differentially expressed between both datasets with a stringent cut-off for adjusted p-value (p < 0.01) and log2FC > 2 (Fig. 1B). Seven miRNAs (miR-140, miR-143, miR-145, miR-30a, miR-21, let-7i, miR-423) were found to have significant differential expression in both OSCC patient tissue and salivary exosomes compared to controls using qRT-PCR (Fig. 1C). Since the differential expression of these miRNAs is significant in both tissue and salivary exosomes of OSCC patients, we tested their utility as biomarkers for detecting OSCC patients.
A three miRNA signature is able to detect OSCC patients
ROC curves were plotted using the different combinations of these 7 miRNAs to determine their utility in detecting OSCC patients. The combination of miR-140, miR-143, and miR-145 showed an AUC of 0.93 (p < 0.0001) with 93.02% sensitivity and 93.33% specificity in tissue samples compared to the various other combinations of 7 miRNAs. Additionally, salivary exosomes revealed an AUC of 0.99 (p < 0.0001) with 98% sensitivity and 99% specificity using the miR-143/miR-145/miR140 combination (Fig. 2A).
PCA models were constructed for salivary exosomes and tissues derived from OSCC patients as compared to their healthy counterparts, using the expression levels of the 3-miRNA combination. The model was able to segregate discrete clusters of OSCC patients and healthy controls, suggesting the effectiveness of 3-miRNA signature in differentiating cancerous patients from non-cancerous controls (Fig. 2B).
Three-miRNA signature predicts overall survival and clinically correlates with disease progression and relapse
Survival analysis was conducted based on the expression patterns of miR-143, miR-145, and miR-140 available in TCGA datasets. OSCC patients with low expression of miR-143, miR-145, and miR-140 had a significantly shorter overall survival as compared to those with high expression (p < 0.05) (Fig. 2C). On assessing the clinical correlation of the signature we observed a significant downregulation of miR-145/miR-140 and miR-143/miR-140 in the tissue and salivary exosomes of late-stage (stage III and IV; n = 30; p < 0.05) OSCC patients respectively (Fig. 2D,E). Additionally, miR-143, miR-145, and miR-140 were significantly downregulated (p < 0.05) in recurrent OSCC tumours (n = 8) as compared to non-recurrent cases (n = 8) (Fig. 2F). Collectively these findings suggest that the expression of the 3-miRNA signature can predict poor patient survival, disease progression, and recurrence in OSCC patients.
Transcriptome sequencing, TCGA datasets and target prediction analysis identifies regulatory miRNA-mRNA networks responsible for OSCC disease progression
To understand the gene regulatory networks modulated by the identified miRNAs, matched samples of OSCC patients were subjected to transcriptome sequencing. RNA-Seq analysis identified 10,660 and 11,632 differentially expressed mRNAs in tumour and leukoplakia samples of OSCC patients, respectively, compared to their non-cancerous counterparts. On comparing these results with the TCGA datasets, we found a total of 5868 mRNAs that were significantly expressed (log2FC > 2; adj p-value < 0.01) across these datasets irrespective of their disease progression and sequencing platforms.
Target genes of identified miRNA signature were predicted using TargetScan, miRDB and miRWalk software. A total of 315 genes were predicted as potential targets that were differentially expressed in leukoplakia and OSCC patients. On constructing a protein-protein interaction network of the identified genes, three enriched networks (Z-scores: 41, 36, and 35) were generated by IPA (Additional File 1: Table S2). On merging these networks, we identified top 16 nodes that were further substantiated using the STRING database and CytoHubba plugin in Cytoscape (Additional File 1: Fig. S3, Table S3). Amongst these, five genes viz. HIF1A, CDH1, CD44, EGFR, and CCND1 formed dense hubs which are putatively reported to be responsible for OSCC progression (Fig. 3A) [20,21,22,23,24,25,26,27].
Subsequently, the differential expression profiles and survival analysis of top 16 hub genes were assessed. A significant upregulation (-log2 (TPM + 1) = ± 2, p < 0.01) was observed in LOX, STAT1, PDGFRB, JAG1, HIF1A, GJA1, ETS1, FLT1 genes in OSCC patients (Additional File 1: Fig. S4). Moreover, Kaplan-Meier analysis suggested that higher expression of CD44 (p < 0.01), EGFR (p < 0.01), GJA1 (p < 0.03), HNF4A (p < 0.01), and CCND1 (p < 0.001) significantly correlated with poor overall survival of patients (Additional File 1: Fig. S5). Further, pathway analysis of the identified miRNA-mRNA networks using IPA, identified epithelial-mesenchymal transition (EMT) as the most statistically significant pathway involving a higher number of overlapping genes (Fig. 3B). Collectively, these findings indicate the relevance of identified 3-miRNA signature and its associated gene networks in disease progression and poor prognosis of OSCC patients.
Overexpression of miR-143/miR-145/miR-140 inhibits cell proliferation and induces G2/M cell cycle arrest
To understand the effect of 3-miRNA signature on functional mechanisms in OSCC, miR-143,miR-145, and miR-140 mimic combination was transfected in OECM-1 cells at 100nM concentration. The transfection efficacy of this mimic combination was measured by real-time PCR. Our results demonstrated a significant upregulation of miR-143, miR-145, and miR-140 by 5.13 ± 0.9 fold, 18.04 ± 2.55 fold, and 23 ± 2.46 fold respectively (Fig. 4A). Consequently, to evaluate the effect of 3-miRNA signature on cell proliferation, MTT assay was performed. A 20% reduction in cell proliferation (p < 0.01) compared to scrambled controls was observed 48 h post-transfection (Fig. 4B).
Cell cycle distribution of OECM1 cells transfected with mimics was measured using flow cytometry. Post transfection, a significant increase was observed in the G2/M (41.73 ± 5%) phase compared to the SC (20.93 ± 2.1%; p < 0.01). Simultaneously, there was a significant decrease in cells in the G0/G1 phase (41.16 ± 4.8%) compared to SC (58.1 ± 6%; p < 0.05) (Fig. 4C, D). Collectively these results suggest that miR-140, miR-143, and miR-145 restrict the proliferative capacity of tumour cells by inducing G2/M arrest.
Upregulation of 3-miRNA signature triggers programmed cell death
To elucidate the effect of miR-143/miR-145/miR-140 combination on apoptotic cell death, we performed Annexin-V/PI staining 48 h post-transfection. Upregulation of the 3-miRNA signature resulted in 17.33 ± 3.36% Annexin+/PI- cells compared to 4.33 ± 2.48% Annexin+/PI- cells (p < 0.01) in SC, indicating increased apoptosis in the mimic transfected cells (Fig. 5A, B). Subsequently, nuclear staining demonstrated a significant number of PI-stained cells with nuclear fragmentation and chromatin condensation which are prominent characteristics of apoptosis in mimic transfected cells compared to SC (Fig. 5C). Based on these findings, it was necessary to explore the mechanism by which this 3-miRNA signature induced apoptotic cell death. Considering that alteration in the mitochondrial membrane is a crucial event of intrinsic apoptosis, we examined the reduction in mitochondrial membrane potential using JC-1 staining. Our analysis indicated that the transfected cells exhibited a significant reduction in the MMP which was signified by a decrease in red fluorescence and a relative increase in green fluorescence compared to SC (Fig. 5D). Collectively, these findings suggest elevated levels of miR-143/miR-145/miR-140 affects the mitochondrial function and drives oral cancer cells towards mitochondria-mediated apoptosis.
Enhanced expression of 3-miRNA signature reduces the migratory and invasive potential of OECM-1 cells by reversing epithelial-mesenchymal transition
Given the aggressive nature of oral cancer, we examined the effect of miR-143/miR-145/miR-140 combination on the migratory and invasive potential of OECM-1 cells. On subjecting the cells to wound healing assay, a significant reduction was observed in the migratory ability of mimic transfected cells (230.85 ± 30.56 μm) as compared to SC (420.44 ± 35.01 μm; p < 0.01) (Fig. 6A, B). Additionally, qualitative and quantitative analysis of trans-well assay demonstrated a significant reduction in the invasive potential of mimic transfected OECM-1 cells (42.34 ± 18%) as compared to SC (81.46 ± 6%; p < 0.05) (Fig. 6C, D). Hence, these results along with the clinical findings and pathway analysis highlighted the role of miR-143/miR-145/miR-140 in disease progression and aggressiveness and prompted the assessment of its effect on EMT phenomena. Microscopic examination demonstrated an epithelial morphology in a majority of transfected cells as compared to an elongated, mesenchymal-like morphology in SC cells thereby suggesting a reversal of EMT in OECM1 cells (Fig. 6E). Further, E-cadherin and N-cadherin expression patterns were estimated in OECM-1 cells pre- and post-transfection using qRT-PCR. E-cadherin and N-cadherin expression decreased significantly compared to SC post-transfection with miRNA mimics (Fig. 6F). Therefore, these results show that overexpression of 3-miRNA signature has a vital role in inhibiting cell migration and invasion by inducing mesenchymal-epithelial transition.
OSCC has a dismal survival rate owing to late-stage tumour detection making it imperative to identify and develop novel non-invasive, high-screening value tests. The utility of miRNAs as biomarkers of disease progression holds a unique opportunity as they regulate gene transcripts during cellular processes and can serve as indicators of disease status in body fluids. In this study, we conducted a comprehensive analysis of small RNA and transcriptome profiles derived from matched salivary exosomes and tumour tissue samples of OSCC patients. Our small RNA sequencing analysis identified a 3-miRNA signature (miR-143, miR-145, and miR-140) that demonstrated significantly higher sensitivity and specificity in effectively differentiating cancerous patients from their control counterparts. On further assessing the clinical relevance of this signature, a significant association was observed between suppression of these miRNAs with disease progression and relapse in oral cancer patients. Thus, these findings reveal the potential clinical utility of this 3-miRNA signature in monitoring disease progression, poor prognosis and relapse in oral cancers. In corroboration with our findings, enhanced expression of miR-140 has been reported to significantly correlate with improved overall survival of head and neck cancer patients . A study by Bufalino et al. has revealed that suppression of tumour-derived miR-143/145 cluster clinically correlated with the progression of the disease, lymph node involvement, and poor survival rates in oral cancer patients . Moreover, downregulation of exosomal miR-145 has been previously identified as an important marker for early detection and risk prediction of oral cancer . Some of the major confounders of these studies are- small unmatched patient cohorts, the inability to monitor the disease in real-time and the inability to explain the underlying mechanisms responsible for regulating OSCC. To the best of our knowledge, this is one of the preliminary studies that highlight the potential clinical utility of this 3-miRNA signature in predicting disease progression using a liquid biopsy approach.
Given that 60% of protein-coding genes in humans are regulated by miRNAs, it is imperative to elucidate the effect of miRNA-mRNA networks on cellular function [8, 9]. These miRNA-mRNA interactions and their underlying mechanisms are responsible for the transformation of normal epithelial cells into a malignant phenotype; however, this interplay is completely unexplored in OSCC. Thus, in this study, we performed paired small RNA and transcriptome sequencing to explore the potential of these miRNAs in modulating vital genes responsible for the progression of OSCC. Integrated network analysis revealed that the identified 3-miRNA signature could potentially target 16 hub genes. Among these, HIF1α, CDH1, CD44, EGFR, and CCND1 were identified as the top hub genes that have also been reported as driver candidates responsible for oral cancer initiation and progression [31,32,33,34]. Interestingly, miR-143 has been reported to target CD44v3 which has a vital role in inhibiting migration and invasion of oral cancer cells . Additionally, Shao et al. demonstrated that miR-145 targeted cyclin D1 expression could effectively suppress OSCC cell growth by regulating cell cycle . Collectively, these findings substantiate that miR-143 and miR-145 target a few of the identified genes and have a vital role in OSCC progression. In terms of the miRNA-mRNA association, miR-140 has been reported to suppress invasive and angiogenic properties by altering EGFR levels in lung cancer . However, similar interactions between EGFR and miR-140 have not been explored in oral cancers.
Alterations in miRNA-mRNA interactions have an impact on signaling pathways and their controlled feedback mechanism resulting in OSCC progression [38,39,40]. Therefore, we explored the downstream signaling pathways modulated by these miRNA-mRNA networks. Our pathway analysis implied that the 3-miRNA signature can potentially indirectly regulate the EMT pathway by targeting 5 hub genes (EGFR, CCND1, CD44, HIF1A and CDH1). EGFR is frequently overexpressed in OSCC and is a primary factor in EMT regulation [41, 42]. It has been established that the internalization of EGFR present in OSCC exosomes leads to the transformation of normal epithelial cells to a mesenchymal phenotype . CDH1, one of the top hub genes, regulates the process of EMT in oral cancers via multiple pathways including the Wnt/β-catenin signalling cascade . In oral cancer, CD44 governs the invasive and metastatic potential of cells, especially the cancer stem cell sub-population by modulating the PI3K/Akt/GSK3β pathway . Interestingly, it has been reported that proteins of the NOTCH receptor family, one of the most frequently altered pathways in HNCs, increase the recruitment of HIF-1α to the LOX promoter region which in turn activates the SNAIL protein, resulting in migration, invasion, and EMT . Furthermore, CCND1 has also been reported to induce EMT by regulating its transcriptional factor SNAIL1/2 and modulating the G1/S phase of the cell cycle mechanism . Collectively, these observations reveal an inevitable role of the 3-miRNA cluster in modulating the EMT mechanism in oral cancers.
Further, we aimed to investigate the functional role of the 3-miRNA signature using an overexpression model in an oral cancer cell line. Our study demonstrated that overexpression of miR-143, miR-145, and miR-140 reduced cell proliferation, induced G2/M cell cycle arrest and triggered apoptosis via the mitochondrial pathway. Our findings corroborated with previous studies that have observed a significant reduction in cell proliferation and apoptosis induced by the overexpression of miR-143 and miR-145 in oral cancers [43, 44]. miR-140-5p overexpression was reported to suppress cell proliferation, promote apoptosis, and hamper cell-cycle progression in OSCC . We also found that overexpression of miR-143, miR-145, and miR-140 induced MET by reducing the invasive and migratory potential in OECM-1 cells. These findings in parallel with pathway analysis suggest a significant role for the 3-miRNA cluster in transcriptionally regulating genes responsible for the EMT process. The involvement of miR-143, miR-145, and miR-140 individually in regulating the process of EMT and metastasis has been reported in gastric and breast cancers [46,47,48]. However, the association of this 3-miRNA signature with disease progression, aggressiveness, and modulation of EMT has not been explored in oral cancers.
Inspite of the immense potential and plausible clinical utility of the identified miRNA signature, this study has substantiated it in a relatively smaller patient cohort. These findings need to be validated in a larger patient cohort using a multi-centric approach for further establishing its clinical utility. Moreover, further identifying the signaling pathways that are modulated by the identified 3-miRNA signature and substantiating these findings using ex-vivo and in-vivo models will strengthen the hypothesis and unravel the functional and biological mechanisms regulated by these miRNAs.
In summary, our findings reveal the diagnostic potential of this 3-miRNA signature using a liquid biopsy approach and uncover a plausible mechanism by which dysregulation of these tumour-suppressive miRNAs could effectively convert a normal epithelial phenotype into oral carcinoma.
The data generated in this study are available upon request from the corresponding authors.
Malik UU, Zarina S, Pennington SR. Oral squamous cell carcinoma: key clinical questions, biomarker discovery, and the role of proteomics. Arch Oral Biol. 2016;63:53–65.
Ge R, Tan E, Sharghi-Namini S, Asada HH. Exosomes in Cancer Microenvironment and Beyond: have we overlooked these extracellular messengers? Cancer Microenvironment. 2012;5(3):323–32.
Osaki M, Okada F. Exosomes and their role in cancer progression. Yonago acta medica. 2019;62(2):182–90.
Osaki M, Okada F. Exosomes and their role in Cancer Progression. Yonago Acta Med. 2019;62(2):182–90.
Patel A, Patel S, Patel P, Tanavde V. Saliva based Liquid Biopsies in Head and Neck Cancer: how far are we from the clinic? Front Oncol. 2022;12:828434.
Ventura TMO, Ribeiro NR, Taira EA, de Lima Leite A, Dionizio A, Rubira CMF, et al. Radiotherapy changes the salivary proteome in head and neck cancer patients: evaluation before, during, and after treatment. Clin Oral Invest. 2022;26(1):225–58.
Liu K, Yang Y. Incorporating link information in feature selection for identifying tumor biomarkers by using miRNA-mRNA paired expression data. Curr Proteomics. 2018;15(2):165–71.
Sun Z, Shi K, Yang S, Liu J, Zhou Q, Wang G, et al. Effect of exosomal miRNA on cancer biology and clinical applications. Mol Cancer. 2018;17(1):1–19.
Prendergast EN, de Souza Fonseca MA, Dezem FS, Lester J, Karlan BY, Noushmehr H, et al. Optimizing exosomal RNA isolation for RNA-Seq analyses of archival sera specimens. PLoS ONE. 2018;13(5):e0196913.
Patel A, Patel S, Patel P, Mandlik D, Patel K, Tanavde V. Salivary exosomal miRNA-1307-5p predicts Disease aggressiveness and poor prognosis in oral squamous cell carcinoma patients. Int J Mol Sci. 2022;23(18).
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnetjournal. 2011;17:10–2.
Andrews S. FastQC: a quality control tool for high throughput sequence data. Cambridge, United Kingdom: Babraham Bioinformatics, Babraham Institute; 2010.
Liao Y, Smyth GK, Shi W. The subread aligner: fast, accurate and scalable read mapping by seed-and-vote. Nucleic Acids Res. 2013;41(10):e108–e.
Anders S, Huber W. Differential expression analysis for sequence count data. Nat Precedings. 2010:1-.
Wang X. miRDB: a microRNA target prediction and functional annotation database with a wiki interface. RNA. 2008;14(6):1012–7.
Dweep H, Sticht C, Pandey P, Gretz N. miRWalk–database: prediction of possible miRNA binding sites by “walking” the genes of three genomes. J Biomed Inform. 2011;44(5):839–47.
Agarwal V, Bell GW, Nam J-W, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. elife. 2015;4:e05005.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Krämer A, Green J, Pollard J Jr, Tugendreich S. Causal analysis approaches in ingenuity pathway analysis. Bioinformatics. 2014;30(4):523–30.
Chaw S, Majeed AA, Dalley A, Chan A, Stein S, Farah C. Epithelial to mesenchymal transition (EMT) biomarkers–E-cadherin, beta-catenin, APC and vimentin–in oral squamous cell carcinogenesis and transformation. Oral Oncol. 2012;48(10):997–1006.
Huang C, Yoon C, Zhou X-H, Zhou Y-C, Zhou W-W, Liu H, et al. ERK1/2-Nanog signaling pathway enhances CD44 (+) cancer stem-like cell phenotypes and epithelial-to-mesenchymal transition in head and neck squamous cell carcinomas. Cell Death Dis. 2020;11(4):1–14.
Miyazaki H, Takahashi R-u, Prieto-Vila M, Kawamura Y, Kondo S, Shirota T, et al. CD44 exerts a functional role during EMT induction in cisplatin-resistant head and neck cancer cells. Oncotarget. 2018;9(11):10029.
Chen J, Zhou J, Lu J, Xiong H, Shi X, Gong L. Significance of CD44 expression in head and neck cancer: a systemic review and meta-analysis. BMC Cancer. 2014;14(1):1–9.
Kashyap T, Pramanik KK, Nath N, Mishra P, Singh AK, Nagini S, et al. Crosstalk between Raf-MEK-ERK and PI3K-Akt-GSK3β signaling networks promotes chemoresistance, invasion/migration and stemness via expression of CD44 variants (v4 and v6) in oral cancer. Oral Oncol. 2018;86:234–43.
Maggiolini M, Musti AM, De Francesco EM. Crosstalk between Notch, HIF-1α and GPER in breast Cancer EMT. Int J Mol Sci. 2018;19(7).
Chong ZX, Yeap SK, Ho WY. Unraveling the roles of miRNAs in regulating epithelial-to-mesenchymal transition (EMT) in osteosarcoma. Pharmacol Res. 2021;172:105818.
Fujiwara T, Eguchi T, Sogawa C, Ono K, Murakami J, Ibaragi S, et al. Anti-EGFR antibody cetuximab is secreted by oral squamous cell carcinoma and alters EGF-driven mesenchymal transition. Biochem Biophys Res Commun. 2018;503(3):1267–72.
Zheng M, Liu J, Meng C, Tang K, Liao J. Prognostic and clinicopathological importance of microRNA-140 expression in cancer patients: a meta-analysis. World J Surg Oncol. 2021;19(1):1–9.
Bufalino A, Cervigne NK, De Oliveira CE, Fonseca FP, Rodrigues PC, Macedo CCS, et al. Low miR-143/miR-145 cluster levels induce activin a overexpression in oral squamous cell carcinomas, which contributes to poor prognosis. PLoS ONE. 2015;10(8):e0136599.
Hang W, Feng Y, Sang Z, Yang Y, Zhu Y, Huang Q, et al. Downregulation of mir-145-5p in cancer cells and their derived exosomes may contribute to the development of ovarian cancer by targeting CT. Int J Mol Med. 2019;43(1):256–66.
Chien H-T, Cheng S-D, Liao C-T, Wang H-M, Huang S-F. Amplification of the EGFR and CCND1 are coordinated and play important roles in the progression of oral squamous cell carcinomas. Cancers. 2019;11(6):760.
Frohwitter G, Buerger H, Korsching E, van Diest PJ, Kleinheinz J, Fillies T. Site-specific gene expression patterns in oral cancer. Head Face Med. 2017;13(1):1–9.
Wen G, Wang H, Zhong Z. Associations of RASSF1A, RARβ, and CDH1 promoter hypermethylation with oral cancer risk: a PRISMA-compliant meta-analysis. Medicine. 2018;97(11).
Judd NP, Winkler AE, Murillo-Sauca O, Brotman JJ, Law JH, Lewis JS, et al. ERK1/2 regulation of CD44 modulates oral cancer aggressiveness. Cancer Res. 2012;72(1):365–74.
Xu P, Li Y, Yang S, Yang H, Tang J, Li M. Micro-ribonucleic acid 143 (MiR-143) inhibits oral squamous cell carcinoma (OSCC) cell migration and invasion by downregulation of phospho-c-Met through targeting CD44 v3. Oral surgery, oral medicine, oral pathology and oral radiology. 2015;120(1):43–51.
Shao Y, Qu Y, Dang S, Yao B, Ji M. MiR-145 inhibits oral squamous cell carcinoma (OSCC) cell growth by targeting c-Myc and Cdk6. Cancer Cell Int. 2013;13(1):1–9.
Flamini V, Dudley WGJ, Cui Y. Distinct mechanisms by which two forms of miR-140 suppress the malignant properties of lung cancer cells. Oncotarget. 2018;9(92):36474.
Huang S-F, Cheng S-D, Chien H-T, Liao C-T, Chen I-H, Wang H-M, et al. Relationship between epidermal growth factor receptor gene copy number and protein expression in oral cavity squamous cell carcinoma. Oral Oncol. 2012;48(1):67–72.
Lin L-H, Lin J-S, Yang C-C, Cheng H-W, Chang K-W, Liu C-J. Overexpression of platelet-derived growth factor and its receptor are correlated with oral tumorigenesis and poor prognosis in oral squamous cell carcinoma. Int J Mol Sci. 2020;21(7):2360.
Subarnbhesaj A, Miyauchi M, Chanbora C, Mikuriya A, Nguyen PT, Furusho H, et al. Roles of VEGF-Flt-1 signaling in malignant behaviors of oral squamous cell carcinoma. PLoS ONE. 2017;12(11):e0187092.
Laimer K, Spizzo G, Gastl G, Obrist P, Brunhuber T, Fong D, et al. High EGFR expression predicts poor prognosis in patients with squamous cell carcinoma of the oral cavity and oropharynx: a TMA-based immunohistochemical analysis. Oral Oncol. 2007;43(2):193–8.
Lo H-W, Hsu S-C, Xia W, Cao X, Shih J-Y, Wei Y, et al. Epidermal growth factor receptor cooperates with signal transducer and activator of transcription 3 to induce epithelial-mesenchymal transition in cancer cells via up-regulation of TWIST gene expression. Cancer Res. 2007;67(19):9066–76.
Sun X, Zhang L. MicroRNA-143 suppresses oral squamous cell carcinoma cell growth, invasion and glucose metabolism through targeting hexokinase 2. Biosci Rep. 2017;37(3).
Irani S, Shokri G. The role of miR-143, miR-145, and miR-590 in expression levels of CD44 and vascular endothelial cadherin in oral squamous cell carcinoma. Middle East Journal of Cancer. 2019;10(3):194–204.
Peng M, Pang C. MicroRNA-140‐5p inhibits the tumorigenesis of oral squamous cell carcinoma by targeting p21‐activated kinase 4. Cell Biol Int. 2020;44(1):145–54.
Lei C, Du F, Sun L, Li T, Li T, Min Y, et al. miR-143 and miR-145 inhibit gastric cancer cell migration and metastasis by suppressing MYO6. Cell Death Dis. 2017;8(10):e3101–e.
Fang Z, Yin S, Sun R, Zhang S, Fu M, Wu Y, et al. Mir-140-5p suppresses the proliferation, migration and invasion of gastric cancer by regulating YES1. Mol Cancer. 2017;16(1):1–11.
Lu Y, Qin T, Li J, Wang L, Zhang Q, Jiang Z, et al. MicroRNA-140-5p inhibits invasion and angiogenesis through targeting VEGF-A in breast cancer. Cancer Gene Ther. 2017;24(9):386–92.
The authors are grateful to Prof. Rakesh Rawal, Gujarat University for kindly providing access to the NanoSight LM10. We acknowledge Dept. of Bioinformatics, Gujarat University for providing computational resources for data analysis. This facility was funded by the Gujarat State Biotechnology Mission and Gujarat Council on Science and Technology, Government of Gujarat. This work was partially supported by the DBT/Wellcome Trust India Alliance Fellowship/Grant [grant number IA/E/20/1/505666] awarded to SP. KS was supported by the DBT-Ramalinagaswami fellowship (DBT-RLS/20-21/04_KS_03.25). This manuscript has been released as a pre-print at https://www.biorxiv.org/content/https://doi.org/10.1101/2022.05.31.494114v2.
This work was supported by The Gujarat State Biotechnology Mission- Financial Assistance Programme under Grant Number: 0P3C9K; and Ahmedabad University Start-up grant under the Grant Number AU/SUG/SAS/BLS/2018-19/04-Vivek Tanavde_03.22. Financial Assistantship to SP from the National Postdoctoral Fellowship Program by the Department of Science and Technology – Science and Engineering Research Board, Government of India (PDF/2017/001912) (2018–2020) is gratefully acknowledged. The Ahmedabad University Assistantship program is appreciated for supporting AP with a doctoral fellowship.
Ethical approval and Consent to participate
All procedures of human samples were conducted after the approval of the HCG Cancer Centre’s Ethics Committee (ECR/92/Inst/GJ/2013/RR-16) for human subject research. This study also complied with the guidelines set forth by the Declaration of Helsinki (2008). All patients provided written informed consent for their participation in the study and their identities have been anonymised.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
About this article
Cite this article
Patel, A., Patel, P., Mandlik, D. et al. A novel 3-miRNA network regulates tumour progression in oral squamous cell carcinoma. Biomark Res 11, 64 (2023). https://doi.org/10.1186/s40364-023-00505-5