- Research article
- Open Access
Salivary extracellular vesicle-associated miRNAs as potential biomarkers in oral squamous cell carcinoma
BMC Cancer volume 18, Article number: 439 (2018)
Several studies in the past have investigated the expression of micro RNAs (miRNAs) in saliva as potential biomarkers. Since miRNAs associated with extracellular vesicles (EVs) are known to be protected from enzymatic degradation, we evaluated whether salivary EVs from patients with oral squamous cell carcinoma (OSCC) were enriched with specific subsets of miRNAs.
OSCC patients and controls were matched with regards to age, gender and risk factors. Total RNA was extracted from salivary EVs and the differential expression of miRNAs was evaluated by qRT-PCR array and qRT-PCR. The discrimination power of up-regulated miRNAs as biomarkers in OSCC patients versus controls was evaluated by the Receiver Operating Characteristic (ROC) curves.
A preliminary qRT-PCR array was performed on samples from 5 OSCC patients and 5 healthy controls whereby a subset of miRNAs were identified that were differentially expressed. On the basis of these results, a cohort of additional 16 patients and 6 controls were analyzed to further confirm the miRNAs that were up-regulated or selectively expressed in the previous pilot study. The following miRNAs: miR-302b-3p and miR-517b-3p were expressed only in EVs from OSCC patients and miR-512-3p and miR-412-3p were up-regulated in salivary EVs from OSCC patients compared to controls with the ROC curve showing a good discrimination power for OSCC diagnosis. The Kyoto Encyclopedia of Gene and Genomes (KEGG) pathway analysis suggested the possible involvement of the miRNAs identified in pathways activated in OSCC.
In this work, we suggest that salivary EVs isolated by a simple charge-based precipitation technique can be exploited as a non-invasive source of miRNAs for OSCC diagnosis. Moreover, we have identified a subset of miRNAs selectively enriched in EVs of OSCC patients that could be potential biomarkers.
Oral squamous cell carcinoma (OSCC) is the most frequent cancer of the head and neck [1, 2]. Despite outstanding diagnostic and therapeutic improvements in oncology, OSCC still holds a poor prognosis with an estimated 5-year overall survival rate of 56% both in the United States and Western Europe [1, 2]. Specifically, in northern Italy, the 3-year and 5-year overall survival rate has been estimated to be 57% and 49%  with the latter decreasing dramatically when considering advanced or metastatic cases. Although, in recent years different biological and molecular factors have been described for the prognosis of OSCC, none of them have had a real impact on routine clinical care. Histopathological staging still remains the gold standard for post-operative management and prognosis of the disease . Hence, more reliable and time saving diagnostic tools are needed.
In OSCC, metastasis spreads predominantly via a lymphatic route with cervical lymph nodes (LN) as the first location, whereas metastasis to distant sites is relatively uncommon . Efficient detection and removal of LN metastasis is therefore crucial in the treatment and survival of patients with this form of carcinoma. Therefore, great expectations lie on the identification of specific predictive factors that could be used clinically for the diagnosis, prognosis and monitoring of the therapeutic response.
In the last decade, extracellular vesicles (EVs) have gained significant attention as a conceivable source of biomarkers. These small membrane-bound vesicles are categorized into 3 different types: exosomes, microvesicles or ectosomes, and apoptotic bodies . EVs are secreted under different physiological and pathophysiological conditions into the extracellular milieu by a variety of cell types, including tumor cells. Tumor-derived EVs have been identified to influence the tumor environment by promoting cancer progression, survival, invasion, and angiogenesis . However, as EVs carry biologically active proteins and nucleic acids from the parent cells, tumor-derived EVs could therefore act as molecular signatures of cancer cells from which they are derived [7, 8].
MicroRNAs (miRNAs) are small non-coding RNA molecules approximately 22 nucleotides in length, which act as regulatory gatekeepers of coding genes. MiRNAs are expressed in a tissue-specific manner, and changes in their expression within a tissue can be correlated with a disease status . Furthermore, they can also modulate gene expression by regulating mRNA translation and/or degradation depending on complementarity between the miRNA and the mRNA . MiRNAs can be secreted either through EVs and/or by forming protein-miRNA complexes with molecules such as high-density lipoproteins and AGO2, which are part of the RISC complex. However, miRNAs carried in EVs are more stable once released as the encapsulation provides protection from enzyme degradation. This therefore makes them more promising as next-generation biomarkers for cancer diagnosis and prognosis . To date, cancer biomarkers carried by EVs have been studied in several types of tumors [9, 11,12,13], including head and neck cancer . However, no work has been published on the expression of miRNAs in EVs from saliva of patients with OSCC.
The aim of the present study was therefore to evaluate whether salivary EVs of OSCC patients and healthy controls express a different pattern of miRNAs and whether the differentially expressed miRNAs could be applied as potential biomarkers for OSCC.
Selection of patients
The enrolled subjects were attending the Oral Medicine Section of the Department of Surgical Sciences, University of Turin, CIR-Dental School, between January and June 2015. Patients with biopsy-proven OSCC were involved in the study with the exclusion of patients with the following criteria: 1) < 18 years of age, 2) pregnant or breast feeding, 3) psychiatrically or mentally unstable. Local ethical committee approval (n° 310/2015, “A.O.U. Città della Salute e della Scienza di Torino”, Turin, Italy) was obtained and all patients provided written informed consent. Demographic information, age at the time of diagnosis and gender, smoking, tumor site, and TNM classification  were recorded at baseline (Table 1 and Additional file 1: Table S1). Healthy subjects presenting no clinically detectable oral lesions matched for age, gender, and risk factors were recruited as controls (Table 1).
As previously reported , all subjects were asked to refrain from: eating, drinking, or oral hygiene for at least one hour prior to collection (usually between 9 and 11 a.m.); they then rinsed their mouths with water and then waited for at least 5 min before spitting into a 50 ml Falcon tube. Participants were instructed not to cough or forcefully expectorate in order to collect unstimulated saliva samples.
HPV-16 in situ hybridization
In situ hybridization (ISH) for HPV was performed on hematoxylin and eosin sections using the Bond TM Ready-to-Use ISH HPV Probe (Leica Biosystems, Newcastle, UK) which targets the following subtypes: 16, 18, 31, 33, and 51. ISH was carried out following the manufacturer’s instructions on the automated Leica BOND system (BOND-MAX, Leica Biosystems).
The sample of saliva from patients with OSCC and healthy controls was diluted 1:1 with PBS (phosphate buffered saline) and centrifuged at 3000 g for 15 min at room temperature to remove cells, debris and bacteria. The supernatant was filtered with 0.2 μm filters and transferred to a sterile tube after which, a precipitation solution (65 μL per 250 μL of saliva) was added and the mixture incubated at 4 °C overnight. The following day, samples were centrifuged at 3000 g for 30 min to precipitate EVs. The resulting supernatant was removed and samples re-centrifuged at 1500 g for 5 min to remove any remaining supernatant . The pellet was resuspended in either: 100 μL of PBS for NanoSight analysis, or 100 μl of RIPA lysis buffer (Sigma Aldrich, Milan, IT) for protein extraction, or 600 μL of Lysis/Binding Buffer (mirVana Isolation Kit, Thermo Fisher Scientific, Waltham, MA, USA) and stored at − 80 °C for subsequent RNA extraction.
The EV samples isolated from saliva were diluted 1:200 in physiologic solution and analyzed by NanoSight LM-10 (Malvern Instruments Ltd., Malvern, UK). The average number and size of EVs were measured by Nanoparticle Tracking Analysis (NTA) software (Malvern Instruments Ltd). Transmission electron microscopy (TEM) of negatively stained EVs (NanoVan, Nanoprobes, Yaphank, NK, USA) was also performed as described previously  and the images were obtained using the Joel JEM 1010 electron microscope (Jeol, Tokyo, Japan).
Western blot analysis
Protein concentration in EV samples was measured by Bradford assay. Protein samples were loaded on polyacrylamide gel at the concentration of 30 μg/well and separated by SDS/PAGE, using 4–15% precast gel (Mini-PROTEAN® Precast Gels, Bio-Rad, Hercules, CA, USA). Proteins were transferred on nitrocellulose membranes by liquid electrophoresis. Membranes were immunoblotted by polyclonal antibodies anti-CD9, CD63, TSG101, and Alix (Santa Cruz Biotechnologies, Dalls, TX, USA). Protein-bands were detected by chemiluminescent Clarity™ ECL Western Blotting Substrate (Bio-Rad) and analyzed by ChemiDoc™ XRS + System (Bio-Rad).
RNA extraction and quantification
miRNAs were extracted from purified EVs by mirVana Isolation Kit (Thermo Fisher Scientific), according to the manufacturer’s instruction. RNA concentration was measured by Nanodrop ND-1000 (Thermo Fisher Scientific), and the ratio 260/280 and 260/230 showed no contaminations. RNA integrity was assessed by a Bioanalyzer (Agilent, Santa Clara, CA) using the RNA 6000 Pico Kit (Agilent).
miRNA expression analysis by qRT-PCR array
To select differentially expressed miRNAs, qRT-PCR array analysis was performed on EVs isolated from five patients with OSCC and five healthy controls. The concentration of selected RNA samples was up to 20 ng/μl and 50 ng of total RNA were retro-transcribed to cDNA with TaqMan® MicroRNA Reverse Transcription Kit (Thermo Fisher Scientific). cDNA was pre-amplified with Megaplex™ RT Primers, Human Pool Set v3.0 and TaqMan® PreAmp Master Mix (Thermo Fisher Scientific) using a Veriti Thermal Cycler (Thermo Fisher Scientific). The expression profile of a panel of 754 human microRNAs was evaluated by a TaqMan® Array Human MicroRNA Card Set v3.0 (Thermo Fisher Scientific) using the real-time thermal cycler 7900HT (Thermo Fisher Scientific).
On the basis of results obtained from the array, we studied miRNAs up-regulated or selectively expressed by patients in a cohort of additional 16 OSCC patients and 6 controls. 500 ng of total RNA was retro-transcribed to cDNA with miScript II RT Kit (Qiagen, Hilden, D) and evaluated for the expression of five miRNAs up-regulated in patients compared to healthy controls (miR-412-3p, miR-489-3p, miR-512-3p, miR-597-5p, and miR-603). Furthermore, eight miRNAs exclusively expressed only in OSCC patients (miR-27a-3p, miR-302b-3p, miR-337-5p, miR-373-3p, miR-494-3p, miR-517b, and miR-520d-3p, miR-645) were also evaluated. Each sample was run in triplicate and each miRNA-specific primer was run in a separate reaction. SnoRNA RNU6B and miR-191 were used as endogenous control as previously described [18,19,20] due to their stable expression in saliva samples which was also confirmed in the current study by qRT-PCR in tested salivary EV samples.
The qRT-PCR reaction mix was composed of 2 ng of cDNA, 100 nM miScript Universal Primer (Qiagen), 100 nM miRNA-specific primer (Eurofins Genomics, Ebersberg, D), 5 μl QuantiTect SYBR Green PCR Master Mix (Qiagen), and nuclease free water (Qiagen) to reach a final reaction volume of 10 μl. The Real-Time Thermal Cycler Quant Studio 12 k (Thermo Fisher Scientific) was used for analysis.
Enrichment analysis for biological pathways was performed for miRNAs that were found to be up-regulated (p < 0.09) or only expressed by OSCC patients. Kyoto Encyclopedia of Gene and Genomes (KEGG) pathway analysis was performed through the DIANA-mirPath v.3.0  online software and miRNA targets were searched on microT-CDS . Results were merged by a “pathway-union” criterion. The p value was calculated by DIANA software online with: False Discovery Rate (FDR) correction, p value threshold at 0.05 and MicroT threshold at 0.8. Fisher’s exact test was used as the statistical method for the enrichment analysis.
Discrimination power analysis
The discrimination power of the up-regulated miRNAs as biomarkers for OSCC diagnosis was evaluated by the Receiver Operating Characteristic (ROC) curves . The ROC curves were constructed using the relative quantification (RQ) of the expression levels of controls and OSCC patients by the demo version of GraphPad Prism 6.01 software. Sensitivity, specificity, area under curve (AUC) and p value were calculated by the software. The optimal threshold value was decided using Youden’s index (sensitivity + specificity-1) .
For array data analysis, SDS Software v.2.3 (Thermo Fisher Scientific) was used to calculate Raw Ct values, with an automatic baseline and threshold. ExpressionSuite Software 1.1 (Thermo Fisher Scientific) was used to calculate RQ (2-∆∆Ct) values. Data were normalized using global normalization, an algorithm that finds the assays common to every sample and then uses the median Ct of those assays as the normalization factor, on a per sample basis . Ct values > 35 or with Amp score < 0.7 were excluded from the analysis. To identify candidate miRNAs differentially expressed between patients and controls, we selected miRNAs with low p values (p ≤ 0.05). P values were calculated by the ExpressionSuite software using Student’s t-test for sample group comparisons, without multiple test correction. Further to this, we screened the group for miRNAs expressed in every sample and with low variability among the same group.
For qRT-PCR data analysis, Excel software (Microsoft Office 365 ProPlus) was used to calculate ∆Ct, −∆∆Ct, and RQ for patients and controls. Statistical analysis was performed on RQ values through the demo version of GraphPad Prism 6.01 software using an unpaired non-parametric two-sided Mann-Whitney test. Confidence level was set at 95% (p value ≤0.05).
A total of 21 patients with OSCC (12 men and 9 women) were analyzed (Table 1). The TNM staging system identified the following lesion categories: T1 (n = 7), T2 (n = 8), T3 (n = 3), T4 (n = 3); according to the histology, of biopsy specimens, three patients were identified as well differentiated, twelve as moderately differentiated and six as poorly differentiated. The lateral border of the tongue was the most commonly affected site (24%), followed by the floor of the mouth and gingiva (19% respectively), the palate (14%), the pelvis (9.5%) and lastly other sites (19%) (Additional file 1: Table S1). Five, out of 21 patients, were positive for HPV (23.8%). A total of 11 controls (seven men and four women) were also analyzed. The subjects did not show oral lesions, infections, or tumor history. Controls and OSCC patients were matched based on gender, age and risk factor, as shown in Table 1.
According to NanoSight results, EVs isolated from saliva samples through charge-based precipitation appeared as a heterogeneous population with a size ranging from 100 to 300 nm (Fig. 1a). TEM analysis confirmed the characteristic shape, aspects, and dimensions of EVs (Fig. 1b). Furthermore, we observed that the size and concentration of salivary EVs from OSCC patients were slightly increased compared to healthy controls, however, the differences were not statistically significant. Western Blot analysis demonstrated the expression of the typical exosome markers: CD63, CD9, TSG 101, and Alix (Fig. 1c). Additionally, Western blot analysis performed in duplicate on samples from 10 OSCC patients and 8 controls was also positive for the exosome markers: CD63 and TSG 101 marker, confirming the reproducibility of the isolation method (Additional file 2: Figure S1). The expression of the exosome markers was similar for both OSCC patients and controls. According to Bioanalyzer results, RNA cargo is mainly constituted of RNAs ranging from 20 to 200 nucleotides, whereas the ribosomal RNAs 18 s and 28 s were absent (Figure 1d). No differences were observed between RNA profiles of OSCC patients and controls (Additional file 3: Figure S2).
miRNA expression analysis by qRT-PCR array
On analyzing the expression of miRNAs in salivary EVs, we identified five miRNAs to be up-regulated (miR-412-3p, miR-489-3p, miR-512-3p, miR-597-5p, and miR-603), and six miRNAs down-regulated (miR-193b-3p, miR-30e-3p, miR-376c-3p, miR-484, miR-720, and miR-93-3p) in tumor EVs compared to controls (Table 2). Moreover, eight miRNAs were exclusively detected in EVs from OSCC patients, while 14 miRNAs were specific only to EVs from controls (Table 3). The complete qRT-PCR array results are shown in Additional file 4: Table S2, reporting all miRNAs expressed in EVs of both groups of patients.
qRT-PCR expression analysis
After an initial screening by qRT-PCR array, we selected 11 miRNAs for subsequent analysis. We chose the five miRNAs up-regulated in OSCC patients (miR-412-3p, miR-489-3p, miR-512-3p, miR-597-5p, and miR-603) and eight miRNAs expressed only by OSCC patients (miR-27a-3p, miR-302b-3p, miR-337-5p, miR-373-3p, miR-494-3p, miR-517b, and miR-520d-3p, miR-645). The analysis of the up-regulated miRNAs showed a significant up-regulation of miR-412-3p and miR-512-3p in OSCC patients with respect to controls (Fig. 2a). The qRT-PCR analysis of miRNAs detected only in OSCC patients through qRT-PCR array showed that miR-27a-3p, miR-337-5p, miR-373-3p, miR-494-3p, and miR-520d-3p were overexpressed in patients however still present in controls (Fig. 2b). MiR-27a-3p, miR-373-3p and miR-494-3p showed a p value lower than 0.1, indicating a conserved but not statistically significant trend of up-regulation in OSCC patients (Fig. 2b). MiR-302b-3p and miR-517b-3p were confirmed to be expressed only in patients (data not shown), while miR-645 expression level was comparable to controls (data not shown).
Discrimination power of miRNAs as OSCC biomarkers
ROC curves were constructed to evaluate the discrimination power of the two up-regulated miRNAs as potential biomarkers for OSCC diagnosis. For each miRNA, the ROC curves express the sensitivity (true positive rate) versus 1-specificity (false positive rate) at various cut-off values, the AUC, indicating the discrimination power of the biomarker, and the p value (Fig. 2c, d). The optimal threshold value was set as the maximum Youden’s index (sensitivity + specificity-1) represented as a black circle. MiR-512-3p (Fig. 2c) and miR-412-3p (Fig. 2d) showed high sensitivity and specificity, with high AUC values of 0.847 and 0.871 respectively, and p values lower than 0.02.
KEGG pathway enrichment analysis
The four miRNAs (miR-512-3p, miR-412-3p, miR-27a-3p, and miR-494-3p) confirmed to be up-regulated and the two miRNAs (miR-302b-3p and miR-517b-3p) confirmed to be expressed only in OSCC patients by qRT-PCR were selected for KEGG pathway enrichment analysis. Eight pathways were found to be significantly enriched for at least two of the tested miRNAs (Fig. 3a). Furthermore, Fig. 3b shows the number of predicted target genes involved in each pathway and Fig. 3c shows, for each miRNA, the number of predicted target genes in each pathway and the respective p value.
The use of saliva as a diagnostic biofluid has been widely recognized, and it has many advantages over other specimens like blood, exfoliated cells and urine [26, 27]. Salivary biomarkers have the potential to serve as non-invasive, widely accessible screening tools. In fact, the collection is inexpensive and can be easily performed. Identifying the proper salivary biomarker profile could contribute to the current screening method of oral cancer, which is limited to physical exam and biopsy of suspicious lesions .
Several works describe the possibility to detect RNA biomarkers of numerous diseases in saliva [26,27,28] and, more specifically, miRNAs associated with oral cancer [26,27,28,29,30,31]. Evidence demonstrates that it is possible to isolate EV-associated RNA from saliva and oral samples [28, 32, 33]. However, to date, miRNA expression analysis in EVs from OSCC have never been reported.
In this work, according to previous evidence [17, 32,33,34,35], we successfully isolated EVs from saliva. A previous study  showed that most of the salivary RNA was associated with EVs. Through Bioanalyzer RNA profiles, we observed that EVs were enriched with RNAs ranging from 20 to 200 nucleotides whereas ribosomal RNAs were nearly absent. Our results are in accordance with other published data of EV RNA cargo [33, 36, 37]. Molecular analysis of miRNAs revealed an up-regulation of miR-412-3p, miR-512-3p, miR-27a-3p, miR-373-3p, miR-494-3p in salivary EVs from OSCC patients. Furthermore, we found that miR-302b-3p and miR-517b-3p were expressed specifically only in samples from the OSCC group. KEGG pathway enrichment analysis based on predicted miRNA targets provides speculative information of miRNA functions. Eight pathways showed a statistically significant enrichment with each pathway predicted to involve two or more miRNAs. For instance, miR-512-3p and miR-27a-3p could target respectively 7 and 20 genes involved in the ErbB signaling pathway. The pathway is known to promote cell proliferation and survival in several solid tumors  and has been shown to be activated in OSCC as well [39,40,41]. MiR-512-3p, miR-27a-3p, and miR-302b-3p could target respectively 14, 30, and 14 genes regulating proteoglycan in cancer pathways. Evidence has shown that CD44, which can be targeted by both miR-512-3p and miR-302b-3p, and the downstream pathway promote cell invasion and migration upon c-Fos stimulation in OSCC . Increased CD44 expression has been associated with ERK1/2 phosphorylation, and increased tumor aggressiveness . Moreover, high CD44 levels have been described as a characteristic feature of cancer stem-like cells in OSCC . In addition, miR-512-3p, miR-412-3p, miR-27a-3p, and miR-302b-3p could target several genes of the TGFβ signaling pathway, including TGFβR2 gene. Interestingly, it has been previously reported that TGFβR2 is commonly reduced in oral epithelium and stroma in OSCC patients . In line with cancer stem like cells, 46 genes involved in signaling pathways regulating pluripotency of stem cells can also be targeted by miR-512-3p, miR-27a-3p, miR-494-3p, miR-517b-3p, and miR-412-3p. The overexpression of Bmi1, which can be targeted by miR-494 and miR-27a-3p, has been shown to promote formation, growth, migration, and metastasis in a subpopulation of tumor cells of the HNSCC . Basing on these observations, we speculate that the increase of miRNA that target genes involved in tumor progression in salivary EVs might represent a defenses mechanism of tumor cells to eliminate anti-tumor miRNAs. However, due to the speculative nature of this analysis, the relation between miRNAs and target genes and pathways should be further proven experimentally.
To better evaluate the discrimination power of the up-regulated miRNAs as OSCC biomarkers, we constructed ROC curves. MiR-512-3p and miR-412-3p resulted to be either sensitive and specific, as shown by high AUC values (0.847 and 0.871 respectively, with p values < 0.02) and maximum Youden’s Index. This indicates that the two miRNAs are good predictors and can be suggested as new candidate biomarkers for OSCC, which can be evaluated through further studies on a larger population. On the other hand, the up-regulation of miR-27a-3p and miR-494-3p can be used as indicators, but are not sufficient as diagnostic biomarkers. Nevertheless, the involvement of miR-494-3p and miR-27a-3p in OSCC is supported by the literature. MiR-494 has been previously isolated from blood of OSCC patients and proposed as a biomarker . MiR-27a-3p is involved in the progression of OSCC by targeting YAP1 and therefore inhibiting epithelial to mesenchymal transition processes . Furthermore, miR-27a-3p can also target MCPH1, which acts as an onco-suppressor gene . It has also been proposed that high levels of miR-27a increase heat sensitivity in OSCC cells, enhancing hyperthermia-induced-cell death . Moreover, miR-27a-3p seems to play a role in progression and metastasis of nasopharyngeal carcinoma , gastric cancer , esophageal cancer , and has been proposed as an EV-associated biomarker for colorectal cancer . No evidence has been reported confirming the involvement of miR-412-3p, miR-512-3p, miR-302b-3p, and miR-517b-3p in OSCC, thus our work provides new insights about the dysregulation of miRNAs in the tumor environment. It is worth mentioning that miR-512-3p has been reported to be up-regulated in metastatic prostate cancer , and conversely shown anti-tumor activity in non-small cell lung cancer  and hepatocellular carcinoma .
On the other hand, our molecular analysis revealed the expression, although not statistically significant, of several miRNAs that have been reported in the literature as up-regulated in whole saliva or plasma of OSCC patients and have been proposed as biomarkers. For example, we observed the expression of both miR-31-5p and miR-31-3p in OSCC patients, however the levels were comparable with controls (Additional file 4: Table S2). Several studies have reported that miR-31 is overexpressed and/or involved in OSCC [58,59,60,61,62,63]. MiR-184, expressed in OSCC patients and controls, has been reported to have a good diagnostic value . MiR-708 is up-regulated in progressing oral premalignant lesions  and was increased but without any statistically significant up-regulation in our group of patients (Additional file 4: Table S2). The discrepancies among the miRNAs up-regulated in our study and in other studies in the literature may be due to the relatively small number of recruited patients and controls and to the different geographical origin of the group of cases enrolled in each study. Moreover, most of the studies detected miRNAs in whole saliva [18, 20, 29, 30] or plasma , instead of salivary EVs. This may lead to different results since miRNAs can be differentially represented in whole saliva and salivary EVs, as it has been previously described for total plasma or plasma-derived EVs [64, 65].
In this work, we demonstrated the possibility to use salivary EVs as a non-invasive source of miRNAs for OSCC diagnosis and we identified two miRNAs (miR-412-3p and miR-512-3p) overexpressed in OSCC patients and two miRNAs (miR-302b-3p and miR-517b-3p) selectively enriched in EVs from OSCC patients. These four miRNAs have the potential to be used as biomarkers.
Head and neck squamous-cell carcinoma
Kyoto Encyclopedia of Gene and Genomes
Oral squamous cell carcinoma
RNA-induced silencing complex
Receiver Operating Characteristic
Transmission electron microscopy
Moore S, Johnson N, Pierce A. The epidemiology of mouth cancer: a review of global incidence. Oral Dis. 2000;6(2):65–74.
Bagan J, Sarrion G, Jimenez Y. Oral cancer: clinical features. Oral Oncol. 2010;46(6):414–7.
Arduino PG, Carrozzo M, Chiecchio A, Broccoletti R, Tirone F, Borra E, et al. Clinical and histopathologic independent prognostic factors in oral squamous cell carcinoma: a retrospective study of 334 cases. J Oral Maxillofac Surg. 2008;66(8):1570–9.
Woolgar JA. Histopathological prognosticators in oral and oropharyngeal squamous cell carcinoma. Oral Oncol. 2006;42(3):229–39.
Zanaruddin SN, Saleh A, Yang YH, Hamid S, Mustafa WM, Khairul Bariah AA, et al. Four proteins signature accurately predicts lymph node metastasis and survival in oral squamous cell carcinoma. Hum Pathol. 2013;44(3):417–26.
Raposo G, Stoorvogel W. Extracellular vesicles: exosomes, microvesicles, and friends. J Cell Biol. 2013;200(4):373–83.
Javeed N, Mukhopadhyay D. Exosomes and Their role in the micro−/macro-environment: a comprehensive review. J Biomed Res 2016 Mar 15;30. https://doi.org/10.7555/JBR.30.20150162. [Epub ahead of print] Review. PubMed PMID: 2829018210.
Kinoshita T, Yip KW, Spence T, Liu FF. MicroRNAs in extracellular vesicles: potential cancer biomarkers. J Hum Genet 2017 Jan;62(1):67–74. https://doi.org/10.1038/jhg.2016.87. Epub 2016 Jul 7. Review. PubMed PMID: 27383658.
Hutvagner G, Zamore PD. A microRNA in a multiple-turnover RNAi enzyme complex. Science. 2002;297(5589):2056–60.
He Y, Lin J, Kong D, Huang M, Xu C, Kim TK, Etheridge A, Luo Y, Ding Y, Wang K. Current state of circulating MicroRNAs as Cancer biomarkers. Clin Chem 2015 Sep;61(9):1138–1155. https://doi.org/10.1373/clinchem.2015.241190. Epub 2015 Aug 3. Review.
Nawaz M, Camussi G, Valadi H, Nazarenko I, Ekström K, Wang X, et al. The emerging role of extracellular vesicles as biomarkers for urogenital cancers. Nat Rev Urol. 2014;11(12):688–701.
Sadovska L, Eglītis J, Linē A. Extracellular vesicles as biomarkers and therapeutic targets in breast Cancer. Anticancer Res. 2015;35(12):6379–90.
Melo SA, Luecke LB, Kahlert C, Fernandez AF, Gammon ST, Kaye J, et al. Glypican-1 identifies cancer exosomes and detects early pancreatic cancer. Nature. 2015;523(7559):177–82.
Bergmann C, Strauss L, Wieckowski E, Czystowska M, Albers A, Wang Y, et al. Tumor-derived microvesicles in sera of patients with head and neck cancer and their role in tumor progression. Head Neck. 2009;31(3):371–80.
Hermanek P, Sobin L. Classification of malignant tumors. Forth. Berlin: Springer-Verlag; 1988.
Arduino PG, Menegatti E, Cappello N, Martina E, Gardino N, Tanteri C, et al. Possible role for interleukins as biomarkers for mortality and recurrence in oral cancer. Int J Biol Markers. 2015;30(2):262–6.
Deregibus MC, Figliolini F, D'Antico S, Manzini PM, Pasquino C, Lena De, et al. Charge-based precipitation of extracellular vesicles. Int J Mol Med. 2016;38(5):1359–66.
Park NJ, Zhou H, Elashoff D, Henson BS, Kastratovic DA, Abemayor E, et al. Salivary microRNA: discovery, characterization, and clinical utility for oral cancer detection. Clin Cancer Res. 2009;15(17):5473–7.
Sauer E, Madea B, Courts C. An evidence based strategy for normalization of quantitative PCR data from miRNA expression analysis in forensically relevant body fluids. Forensic Sci Int Genet. 2014;11:174–81.
Momen-Heravi F, Trachtenberg AJ, Kuo WP, Cheng YS. Genomewide study of salivary MicroRNAs for detection of oral Cancer. J Dent Res. 2014;93(7 Suppl):86S–93S.
Vlachos IS, Zagganas K, Paraskevopoulou MD, Georgakilas G, Karagkouni D, Vergoulis T, et al. DIANA-miRPath v3.0: deciphering microRNA function with experimental support. Nucleic Acids Res. 2015;43(1):460–6.
Paraskevopoulou MD, Georgakilas G, Kostoulas N, Vlachos IS, Vergoulis T, Reczko M, Filippidis C, Dalamagas T, Hatzigeorgiou AG. DIANA-microT web server v5.0: service integration into miRNA functional analysis workflows. Nucleic Acids Res. 2013 Jul;41(Web Server issue):W169–73. https://doi.org/10.1093/nar/gkt393.
Hajian-Tilaki K. Receiver operating characteristic (ROC) curve analysis for medical diagnostic test evaluation. Caspian J Intern Med. 2013;4(2):627–35.
Ruopp MD, Perkins NJ, Whitcomb BW, Schisterman EF. Youden index and optimal cut-point estimated from observations affected by a lower limit of detection. Biom J. 2008;50(3):419–30.
Mestdagh P, Van Vlierberghe P, De Weer A, Muth D, Westermann F, Speleman F, Vandesompele J. A novel and universal method for microRNA RT-qPCR data normalization. Genome Biol. 2009;10(6):R64. https://doi.org/10.1186/gb-2009-10-6-r64.
Kaczor-Urbanowicz KE, Martin Carreras-Presas C, Aro K, Tu M, Garcia-Godoy F, Wong DT. Saliva diagnostics - current views and directions. Exp Biol Med (Maywood). 2017;242(5):459–72.
Li Y, St John MA, Zhou X, Kim Y, Sinha U, Jordan RC, et al. Salivary transcriptome diagnostics for oral cancer detection. Clin Cancer Res. 2004;10:8442–50.
Michael A, Bajracharya SD, Yuen PS, Zhou H, Star RA, Illei GG, et al. Exosomes from human saliva as a source of microRNA biomarkers. Oral Dis. 2010;16(1):34–8.
Yang Y, Li YX, Yang X, Jiang L, Zhou ZJ, Zhu YQ. Progress risk assessment of oral premalignant lesions with saliva miRNA analysis. BMC Cancer. 2013;13:129.
Zahran F, Ghalwash D, Shaker O, Al-Johani K, Scully C. Salivary microRNAs in oral cancer. Oral Dis. 2015;21(6):739–47.
Zhou Y, Kolokythas A, Schwartz JL, Epstein JB, Adami GR. microRNA from brush biopsy to characterize oral squamous cell carcinoma epithelium. Cancer Med. 2017;6(1):67–78.
Yap T, Vella LJ, Seers C, Nastri A, Reynolds E, Cirillo N, et al. Oral swirl samples - a robust source of microRNA protected by extracellular vesicles. Oral Dis. 2017;23(3):312–7.
Lässer C, Alikhani VS, Ekström K, Eldh M, Paredes PT, Bossios A, Sjöstrand M, Gabrielsson S, Lötvall J, Valadi H. Human saliva, plasma and breast milk exosomes contain RNA: uptake by macrophages. J Transl Med. 2011;9:9. https://doi.org/10.1186/1479-5876-9-9.
Marzesco AM, Janich P, Wilsch-Bräuninger M, Dubreuil V, Langenfeld K, Corbeil D, et al. Release of extracellular membrane particles carrying the stem cell marker prominin-1 (CD133) from neural progenitors and other epithelial cells. J Cell Sci. 2005;118(13):2849–58.
Ogawa Y, Kanai-Azuma M, Akimoto Y, Kawakami H, Yanoshita R. Exosome-like vesicles with dipeptidyl peptidase IV in human saliva. Biol Pharm Bull. 2008;31(6):1059–62.
Willms E, Johansson HJ, Mäger I, Lee Y, Blomberg KE, Sadik M, Alaarg A, Smith CI, Lehtiö J, El Andaloussi S, Wood MJ, Vader P. Cells release subpopulations of exosomes with distinct molecular and biological properties. Sci Rep. 2016;6:22519. https://doi.org/10.1038/srep22519.
Santonocito M, Vento M, Guglielmino MR, Battaglia R, Wahlgren J, Ragusa M, Barbagallo D, Borzì P, Rizzari S, Maugeri M, Scollo P, Tatone C, Valadi H, Purrello M, Di Pietro C. Molecular characterization of exosomes and their microRNA cargo in human follicular fluid: bioinformatic analysis reveals that exosomal microRNAs control pathways involved in follicular maturation. Fertil Steril. 2014;102(6):1751–61.e1. https://doi.org/10.1016/j.fertnstert.2014.08.005.
Appert-Collin A, Hubert P, Crémel G, Bennasroune A. Role of ErbB receptors in Cancer cell migration and invasion. Front Pharmacol. 2015;6:283.
Li X, Sun R, Geng X, Wang S, Zen D, Pei J, et al. A comprehensive analysis of candidate gene signatures in oral squamous cell carcinoma. Neoplasma. 2017;64(2):167–74.
Ohnishi Y, Yasui H, Kakudo K, Nozaki M. Lapatinib-resistant cancer cells possessing epithelial cancer stem cell properties develop sensitivity during sphere formation by activation of the ErbB/AKT/cyclin D2 pathway. Oncol Rep. 2016;36(5):3058–64.
Li SX, Yang YQ, Jin LJ, Cai ZG, Sun Z. Detection of survivin, carcinoembryonic antigen and ErbB2 level in oral squamous cell carcinoma patients. Cancer Biomark. 2016;17(4):377–82.
Dong C, Ye DX, Zhang WB, Pan HY, Zhang ZY, Zhang L. Overexpression of c-fos promotes cell invasion and migration via CD44 pathway in oral squamous cell carcinoma. J Oral Pathol Med. 2015;44(5):353–60.
Judd NP, Winkler AE, Murillo-Sauca O, Brotman JJ, Law JH, Lewis JS Jr, et al. ERK1/2 regulation of CD44 modulates oral cancer aggressiveness. Cancer Res 2012;72(1):365–374.
Ghuwalewala S, Ghatak D, Das P, Dey S, Sarkar S, Alam N, et al. CD44(high)CD24(low) molecular signature determines the Cancer stem cell and EMT phenotype in oral squamous cell carcinoma. Stem Cell Res. 2016;16(2):405–17.
Meng W, Xia Q, Wu L, Chen S, He X, Zhang L, et al. Downregulation of TGF-beta receptor types II and III in oral squamous cell carcinoma and oral carcinoma-associated fibroblasts. BMC Cancer. 2011;11:88.
Liu W, Feng JQ, Shen XM, et al. Two stem cell markers, ATP-binding cassette, G2 subfamily (ABCG2) and BMI-1, predict the transformation of oral leukoplakia to cancer: a long-term follow-up study. Cancer. 2012;118:1693–700.
Ries J, Vairaktaris E, Agaimy A, Kintopp R, Baran C, Neukam FW, et al. miR-186, miR-3651 and miR-494: potential biomarkers for oral squamous cell carcinoma extracted from whole blood. Oncol Rep. 2014;31(3):1429–36.
Zeng G, Xun W, Wei K, Yang Y, Shen H. MicroRNA-27a-3p regulates epithelial to mesenchymal transition via targeting YAP1 in oral squamous cell carcinoma cells. Oncol Rep. 36(3):1475–82.
Venkatesh T, Nagashri MN, Swamy SS, Mohiyuddin SM, Gopinath KS, Kumar A. Primary microcephaly gene MCPH1 shows signatures of tumor suppressors and is regulated by miR-27a in oral squamous cell carcinoma. PLoS One. 2013;8(3):e54643. https://doi.org/10.1371/journal.pone.0054643.
Kariya A, Furusawa Y, Yunoki T, Kondo T, Tabuchi Y. A microRNA-27a mimic sensitizes human oral squamous cell carcinoma HSC-4 cells to hyperthermia through downregulation of Hsp110 and Hsp90. Int J Mol Med. 2014;34(1):334–40.
Li L, Luo Z. Dysregulated miR-27a-3p promotes nasopharyngeal carcinoma cell proliferation and migration by targeting Mapk10. Oncol Rep. 2017; https://doi.org/10.3892/or.2017.5544.
Zhou L, Liang X, Zhang L, Yang L, Nagao N, Wu H, et al. MiR-27a-3p functions as an oncogene in gastric cancer by targeting BTG2. Oncotarget. 2016;7(32):51943–54.
Wu XZ, Wang KP, Song HJ, Xia JH, Jiang Y, Wang YL. MiR-27a-3p promotes esophageal cancer cell proliferation via F-box and WD repeat domain-containing 7 (FBXW7) suppression. Int J Clin Exp Med. 2015;8(9):15556–62.
Ostenfeld MS, Jensen SG, Jeppesen DK, Christensen LL, Thorsen SB, Stenvang J, et al. miRNA profiling of circulating EpCAM(+) extracellular vesicles: promising biomarkers of colorectal cancer. J Extracell Vesicles. 2016;5:31488.
Sadeghi M, Ranjbar B, Ganjalikhany MR, M Khan F, Schmitz U, Wolkenhauer O, Gupta SK. MicroRNA and transcription factor gene regulatory network analysis reveals key regulatory elements associated with prostate Cancer progression. PLoS One. 2016 Dec 22;11(12):e0168760. https://doi.org/10.1371/journal.pone.0168760.
Chen F, Zhu HH, Zhou LF, Wu SS, Wang J, Chen Z. Inhibition of c-FLIP expression by miR-512-3p contributes to taxol-induced apoptosis in hepatocellular carcinoma cells. Oncol Rep. 2010;23(5):1457–62.
Zhu X, Gao G, Chu K, Yang X, Ren S, Li Y, et al. Inhibition of RAC1-GEF DOCK3 by miR-512-3p contributes to suppression of metastasis in non-small cell lung cancer. Int J Biochem Cell Biol. 2015;61:103–14.
Liu CJ, Kao SY, Tu HF, Tsai MM, Chang KW, Lin SC. Increase of microRNA miR-31 level in plasma could be a potential marker of oral cancer. Oral Dis. 2010;16(4):360–4.
Lu WC, Kao SY, Yang CC, Tu HF, Wu CH, Chang KW, et al. EGF up-regulates miR-31 through the C/EBPβ signal cascade in oral carcinoma. PLoS One. 2014;9(9):e108049.
Siow MY, Ng LP, Vincent-Chong VK, Jamaludin M, Abraham MT. Abdul Rahman Zaet al. Dysregulation of miR-31 and miR-375 expression is associated with clinical outcomes in oral carcinoma. Oral Dis. 2014;20(4):345.
Hung PS, Tu HF, Kao SY, Yang CC, Liu CJ, Huang TY, et al. miR-31 is up-regulated in oral premalignant epithelium and contributes to the immortalization of normal oral keratinocytes. Carcinogenesis. 2014;35(5):1162–71.
Yan ZY, Luo ZQ, Zhang LJ, Li J, Liu JQ. Integrated analysis and MicroRNA expression profiling identified seven miRNAs associated with progression of oral squamous cell carcinoma. J Cell Physiol. 2016; https://doi.org/10.1002/jcp.25728.
Hung KF, Liu CJ, Chiu PC, Lin JS, Chang KW, Shih WY, et al. MicroRNA-31 up-regulation predicts increased risk of progression of oral potentially malignant disorder. Oral Oncol. 2016;53:42–7.
Cheng L, Sharples RA, Scicluna BJ, Hill AF. Exosomes provide a protective and enriched source of miRNA for biomarker profiling compared to intracellular and cell-free blood. J Extracell Vesicles. 2014;3 https://doi.org/10.3402/jev.v3.23743.
Endzeliņš E, Berger A, Melne V, Bajo-Santos C, Soboļevska K, Ābols A, et al. Detection of circulating miRNAs: comparative analysis of extracellular vesicle-incorporated miRNAs and cell-free miRNAs in whole plasma of prostate cancer patients. BMC Cancer. 2017;17(1):730. https://doi.org/10.1186/s12885-017-3737-z.
The authors thank Unicyte for providing the precipitation solution to purify EVs.
The authors thanks Dr. Sharad Kholia for the revision of English.
This work was supported by the Associazione Italiana Ricerca Cancro (grant number AIRC IG 2015.16973) to GC and by “Fondi Universitari ex-60% 2014” to PGA. Funding agencies had not any influence in the design of the study, collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and materials
All data generated or analyzed during this study are included in this published article in Additional file 1: Table S1, Additional file 2: Figure S1, Additional file 3: Figure S2 and Additional 4: Table S2.
Ethics approval and consent to participate
The study was approved by local ethical committee “A.O.U. Città della Salute e della Scienza di Torino”, Turin, Italy (n° 310/2015) and all patients gave written informed consent to participate.
GC is named as inventors in a related patent application. All other authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Characteristics and tumor staging of OSCC patients enrolled in the study. (DOCX 16 kb)
Figure S1. Protein expression analysis of exosome markers CD63 and TSG101. For each protein, the picture shows two western blot experiments performed on salivary EVs from 10 OSCC patients (left) and 8 controls (right). (TIF 373 kb)
Figure S2. Bioanalyzer RNA profiles. Profiles of RNA isolated from salivary EVs of four controls (left) and four patients (right) in duplicate (a, b). The graphs show fluorescence intensity [FU]/nucleotide length [nt] and were obtained by bioanalyzer analysis. (TIF 1968 kb)
Table S2. miRNAs expressed in both controls and OSCC patients. (DOCX 44 kb)
About this article
Cite this article
Gai, C., Camussi, F., Broccoletti, R. et al. Salivary extracellular vesicle-associated miRNAs as potential biomarkers in oral squamous cell carcinoma. BMC Cancer 18, 439 (2018). https://doi.org/10.1186/s12885-018-4364-z
- Oral squamous cell carcinoma (OSCC)
- Extracellular vesicles (EVs)
- microRNA (miRNA)
- miRNA-512-3p (miR-512)
- miRNA-412-3p (miR-412)
- miR-302b-3p (miR-302b)
- miR-517b-3p (miR-517b)
- miR-27a-3p (miR-27a)