Salivary MicroRNAs for Early Detection of Head and Neck Squamous Cell Carcinoma: A Case-Control Study in the High Altitude Mestizo Ecuadorian Population

Salazar-Ruales, Carolina; Arguello, Jessica-Viviana; López-Cortés, Andrés; Cabrera-Andrade, Alejandro; García-Cárdenas, Jennyfer M.; Guevara-Ramírez, Patricia; Peralta, Patricio; Leone, Paola E.; Paz-y-Miño, César

doi:https://doi.org/10.1155/2018/9792730

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Data Availability Conflicts of Interest Authors’ Contributions References Copyright Related Articles

Research Article | Open Access

Volume 2018 | Article ID 9792730 | https://doi.org/10.1155/2018/9792730

Salivary MicroRNAs for Early Detection of Head and Neck Squamous Cell Carcinoma: A Case-Control Study in the High Altitude Mestizo Ecuadorian Population

Carolina Salazar-Ruales,¹Jessica-Viviana Arguello,²Andrés López-Cortés,¹Alejandro Cabrera-Andrade,¹Jennyfer M. García-Cárdenas,¹Patricia Guevara-Ramírez,¹Patricio Peralta,³Paola E. Leone,¹and César Paz-y-Miño¹

Academic Editor: Xin-yuan Guan

Received30 Jul 2018

Revised06 Oct 2018

Accepted23 Oct 2018

Published21 Nov 2018

Abstract

Head and neck squamous cell carcinoma (HNSCC) is the sixth most common cancer with the highest incidence worldwide. HNSCC is often diagnosed at advanced stages, incurring significant high mortality and morbidity. The use of saliva, as a noninvasive tool for the diagnosis of cancer, has recently increased. Salivary microRNAs (miRNAs) have emerged as a promising molecular tool for early diagnosis of HNSCC. The aim was to identify the differential expression of salivary miRNAs associated with HNSCC in the high altitude mestizo Ecuadorian population. Using PCR Arrays, miR-122-5p, miR-92a-3p, miR-124-3p, miR-205-5p, and miR-146a-5p were found as the most representative ones. Subsequently, miRNAs expression was confirmed in saliva samples from 108 cases and 108 controls. miR-122-5p, miR-92a-3p, miR-124-3p, and miR-146a-5p showed significant statistical difference between cases and controls with areas under the curve (AUC) of 0.73 (p < 0.001), 0.70 (p < 0.001), 0.71 (p = 0.002), and 0.66 (p = 0.008), respectively. miRNAs were also deregulated in between HNSCC localizations. A differentiated expression of miR-122-5p between oral cancer and oropharynx cancer (AUC of 0.96 p = 0.01) was found: miR-124-3p between larynx and pharynx (AUC = 0.97, p < 0.01) and miR-146a-5p between larynx, oropharynx, and oral cavity (AUC = 0.96, p = 0.01). Moreover, miR-122-5p, miR-124-3p, miR-205-5p, and miR-146a-5p could differentiate between HPV+ and HPV- (p=0.004). Finally, the expression profiles of the five miRNAs were evaluated to discriminate HNSCC patient’s tumor stages (TNM 2-4). miR-122-5p differentiates TNM 2 and 3 (p = 0.002, AUC = 0.92), miR-124-3p TNM 2, 3, and 4 (p < 0.001, AUC = 98), miR-146a-5p TNM 2 and 3 (p < 0.001, AUC = 0.97), and miR-92a-3p TNM 3 (p < 0.001, AUC = 0.99). Taken together, these findings show that altered expression of miRNAs could be used as biomarkers for HNSCC diagnosis in the high altitude mestizo Ecuadorian population.

1. Introduction

Head and neck squamous cell carcinoma (HNSCC) is the sixth most common cancer with the highest incidence in the world, with 600,000 new cases per year and a mortality rate between 223,000 and 300,000 deaths [1–5]. The overall incidence of HNSCC is 8.8 and 5.1 cases per 100,000 inhabitants in men and women, respectively [1, 5]. The tumors originate in the epithelial cells of the mucosal linings of the upper airway and food passages (oral cavity, oropharynx, larynx, and pharynx) [6]. Oral cancer is the most common HNSCC. The primary screening test for oral cancer is a systematic examination, including inspection and palpitation of the oral cavity [7]. The highest incidence of oral cavity cancer is found in Southern and Central Asia, Melanesia, Southwestern Europe, and Southern Africa [4]. On the other hand, serum Epstein Barr Virus- (EBV-) associated antibodies and circulating cell-free EBV DNA testing have been used for nasopharyngeal cancer diagnosis and screening [8]. This type of cancer has a higher incidence in Southern China [4]. Additionally, there is no screening test to find early laryngeal cancer at the moment. The only way to screen this type of cancer would be a test called flexible endoscopy, taking a biopsy from the lining of the larynx. Laryngeal cancer has a higher incidence in Southern and Eastern Europe, South America, and Western Asia [4, 9]. In Ecuador, HNSCC is the third most frequent cancer for women and sixth for men, with 1637 new cases per year, being larynx cancer the most common one [10].

The risk factors for developing HNSCC include tobacco consumption, either smoking or chewing, and alcohol consumption. Both of these factors account for nearly 90% of the cases and are associated with age, sex, and ethnicity [11, 12]. Additionally, the Human Papillomavirus (HPV) infection is a risk factor for oropharyngeal cancer and EBV infection for nasopharyngeal cancer. The relative frequency of these risk factors contribute to the variation in the observed distribution of HNSCC worldwide [3, 9].

HNSCC is often diagnosed at advanced stages (T3 or T4) when survival rate is reduced to 20%, incurring significantly high mortality and morbidity [3, 13, 14]. On the other hand, the 5-year survival rate is 80%, if detected at early stages (T1 or T2). Therefore, the availability of biomarkers of different racial/ethnic populations with a potential application for HNSCC screening, early diagnosis, and monitoring of response to therapy are indispensable in clinical practice nowadays [12, 15].

The use of ‘liquid biopsies,’ such as saliva, for the analysis of biomarkers might predict relapse at the earliest stage [6]. Saliva analysis has recently increased as a noninvasive and inexpensive tool for cancer diagnosis. Saliva has advantages at the time of collection, transport, and processing, compared to blood or tissue samples [16, 17]. In addition, salivary transcriptome has been reported with more than 3,000 species of RNAs and microRNAs (miRNAs) that can be used as biomarkers for the noninvasive diagnosis of HNSCC [3, 18–20]. According to Koshizuka et al., 2016, there are several studies of downregulated and upregulated miRNA expression profiles in HNSCC specially located in hypopharyngeal, larynx, oropharynx, maxillary sinus, oral cavity, and nasopharyngeal [21–31]. Hence, miRNAs have emerged as a promising molecular tool for diagnosis, prognosis, and therapy of HNSCC. miRNAs are currently defined as a small noncoding RNA type of 19 to 25 nucleotides in length, which regulate gene expression at the posttranscriptional level by blocking translation or inducing degradation of messenger RNA [32]. Multiple genes are regulated by a single miRNA, and multiple miRNAs can regulate a gene [6]. The specificity and biological properties of miRNAs make them a potential biomarker in the pathogenesis of cancer. miRNAs are expressed intracellular and extracellular in many tissues such as saliva, an easily accessible and noninvasive fluid with high stability [3]. Specific miRNAs and their expression pattern have been established in different primary tumors in HNSCC [33]. Zahran et al., 2015, identified three salivary miRNAs as potential markers of malignant transformation in oral mucosal lesions [34]. The aim of this study was to identify the differential expression of salivary miRNAs in order to develop a noninvasive and inexpensive tool for the diagnosis of HNSCC in the high altitude mestizo Ecuadorian population.

2. Materials and Methods

2.1. Study Subjects

All 216 individuals belong to the high altitude mestizo Ecuadorian population (2,800 meters above sea level). Unstimulated whole mouth resting saliva was collected from 108 HNSCC patients (stages 1 to 4), from the Hospital Oncológico Solón Espinosa Ayala, Quito, Ecuador, between the years 2015 and 2016. Demographic and clinicopathological characteristics of patients are listed in Table 1. Likewise, as for the control group, saliva samples were collected from 108 healthy individuals who were nonsmokers, with good oral hygiene and with no history of the disease. All samples were immediately processed in the laboratory. The Bioethics Committee at the Universidad de las Americas previously approved this study (No. 2013-1101). All participants provided written informed consent before sample collection.

2.2. Total RNA and miRNA Extraction

Total RNA from 200 μL of saliva of 10 healthy controls and 10 patients was extracted using TRIzol, following manufacturer’s instructions (Invitrogen, CA, USA). Total RNA was used to perform profiler PCR Arrays (miScript miRNA PCR Array) (Qiagen, Valencia, CA, USA), selecting the most representative miRNA in the Ecuadorian HNSCC population. In order to validate the selected miRNAs from the PCR Arrays, miRNAs were extracted from 300 μL of saliva from all 216 participants using mirVana™ miRNA Isolation kit (ThermoFisher Scientific, CA, USA), according to the manufacturer’s instructions. miRNAs were eluted in 30 μL of RNase free water. Total RNA and miRNA concentrations were measured using a NanoDrop 2000 (ThermoScientific, Waltham, MA, USA) and samples with OD 260:280 ratios ≥1.8 were included in the study.

2.3. Total RNA and miRNA cDNA Synthesis

The complementary DNA (cDNA) synthesis from total RNA was carried out using the RT² First Standard Kit (Cat. No. 330401) (Qiagen Sciences, Maryland 20874, USA). Prior to cDNA synthesis, a genomic elimination reaction was performed, according to manufactures instructions. Briefly, 2 μL of GE buffer were added to 400 ng of RNA in a final reaction volume of 10 μL. The reactions were incubated 5 min at 42°C in a SureCycler 8800 (Agilent, Santa Clara, CA, USA) and then kept on ice for 5 min. For the reverse transcription, all the genomic DNA elimination mix was converted into cDNA adding 10 μL of the reverse transcription mix and incubated for 15 minutes at 42°C and 5 minutes at 95°C. Then 91 μL of RNase free water were added to the reactions and mixed by pipetting. The reactions were kept at -20°C until the development of the real-time PCR Arrays. miRNA cDNA synthesis was performed with miScript II RT cDNA synthesis kit using Hispec buffer (Qiagen, Valencia, CA, USA), following the manufacturer instructions.

2.4. miRNA Profiler PCR Arrays

To select the most representative miRNAs in the high altitude mestizo Ecuadorian population affected with HNSCC, miScript miRNA PCR Array (Cat. No. MIHS-102Z) (Qiagen, Valencia, CA, USA) was used, according to manufacturer’s instructions. First, cDNA from the total RNA, containing miRNAs, of 10 patients and 10 controls were diluted 1:10 in RNase free water. The miRNAs expressions were run in the Mx3005p qPCR System (Agilent, CA, United States) with the following cycling conditions: 95°C for 15 min, followed by 40 cycles of 94°C for 15 s, 55°C for 30 s, and 70°C for 30 s.

2.5. Relative miRNA Expression and the Comparative Livak Method (2^-ΔΔCT)

Once the miRNA PCR Array was done, we calculated the relative miRNA expression using SNORD96A as the housekeeping normalizer through the Livak method () [35]. Relative quantification relates the PCR signal of the target transcript in a treatment group to that of another sample such as an untreated control [35]. Ct values >35 were excluded and miRNAs were selected based on fold-change >2 and p < 0.05, and p values were calculated based on a Student’s t-test of the replicate values for each miRNA in the case and control groups. Primers used for the selected miRNAs are the following: hsa-miR-92a-3p (UAUUGCACUUGUCCCGGCCUGU), hsa-miR-205-5p (UCCUUCAUUCCACCGGAGUCUG), hsa-miR-124-3p (UAAGGCACGCGGUGAAUGCC), hsa-miR-122-5p (UGGAGUGUGACAAUGGUGUUUG), hsa-miR-146a-5p (UGAGAACUGAAUUCCAUGGGUU), and SNORD96A (CCUGGUGAUGACAGAUGGCAUUGUCAGCCAAUCCCCAAGUGGGAGUGAGGACAUGUCCUGCAAUUCUGAAGG).

2.6. miRNA Quantitative Real-Time PCR

miRNA expression was confirmed by RT-qPCR. Saliva of 108 patients and 108 controls was analyzed with miScript SYBR Green PCR kit (Qiagen, Valencia, CA, USA). RT-qPCR was run in triplicates in a Mx3005p qPCR System (Agilent, CA, United States) with the following conditions: 95°C for 15 min, 40 cycles of 94°C for 15 s, 55°C for 30s,and 70°C for 30s. Ct values >35 were excluded, and method was used to determine fold-changes.

2.7. Statistical Analysis

Relative quantification method was used to compare differences between patients and controls. The data were analyzed with miScript miRNA PCR Array Data Analysis Excel web portal (https://dataanalysis.sabiosciences.com/pcr/arrayanalysis). Additionally, Mann–Whitney U test for nonparametric analysis was applied to compare miRNA expression levels after Livak method [35]. Finally, in order to determinate the utility of the miRNAs as diagnostic biomarkers, receiver-operating characteristics (ROC) curves were generated. Ct values were normalized using SNORD96A as housekeeping gene. All statistical analyses were performed using the Graph Pad Prism 7 software version 7.02 (Graph Pad Software Inc., USA).

3. Results

Out of the 108 patients, 56 were females and 52 males; the median age was 53.5 (26-81) years old. 28.57% of patients were regular smokers (≥ 1 cigarette per day); out of those 3.9% were females and 26.7% males. Alcohol consumption was higher in males (14.3%) than females (6.7%). Moreover, 12.7% of patients were HPV positive (HPV+), 54.5% had tumor stage (TNM) 2, 30.9% TNM 3 and 14.5% TNM 4 (Table 1).

From the PCR Arrays, miR-122-5p, miR-92a-3p, miR-124-3p, miR-205-5p, and miR-146a-5p were selected as the most representative ones in the Ecuadorian HNSCC population based on fold-change >2 (p < 0.05). miR-122-5p and miR-146a-5p were upregulated while miR-92a-3p, miR-124-3p, and miR-205-5p were downregulated (fold-changes 11.5, 10.79, 11.24, 25.51, and 55.7, respectively). Subsequently, expression of the five selected miRNAs was confirmed in saliva samples from 108 patients and 108 controls. Interestingly, all the miRNAs studied showed regulation in the fold-change values, miR-122-5p: 4.47 (p < 0.05), miR-92a-3p: 3.50 (p < 0.05), miR-124-3p: 3.18 (p < 0.01), and miR-146a-5p: 1.12 (p < 0.01) (Figure 1), showing a statistical difference between the miRNAs in patients and controls with the exception of miR-205-5p (p > 0.05).

In order to determinate the discriminatory differentiation of miR-122-5p, miR-92a-3p, miR-124-3p, and miR-146a-5p, receiver operator characteristic (ROC) curves were generated. The four miRNAs showed significant statistical difference with AUC of 0.73 (p < 0.001), 0.70 (p < 0.001), 0.71 (p < 0.05), and 0.66 (p < 0.01), respectively.

The ability of miRNAs to differentiate among HNSCC localizations was also evaluated (Figure 2(a)). miR-122-5p was upregulated in patients with HNSCC in larynx (fold-change of 8.27, p < 0.05) and downregulated in oropharynx HNSCC patients (fold-change of 3.83, p < 0.01), AUC of 0.96 p < 0.01. Moreover, an upregulation was found in miR-124-3p and patients with larynx and pharynx cancer (fold-changes of 5.91 and 2.77, respectively, p < 0.001). On the other hand, a downregulation was found in oropharynx and oral cavity patients (fold-changes of 2.03 and 3.9, respectively, p < 0.01). A differentiation among these localization sites was done with AUC of 0.97. miR-146a-5p could differentiate patients with HNSCC in the larynx, oropharynx, pharynx, and oral cavity (AUC of 0.96 p < 0.01) with upregulation in the larynx (fold-change of 5.98) and downregulation in the oropharynx, pharynx, and oral cavity (fold-changes of 2.21, 7.88, 2.2, respectively, p < 0.01). miR-92a-3p showed no significant difference (p > 0.05) to distinguish HNSCC localizations.

(a)

(b)

(c)

Additionally, we also evaluated the differentiation power of the miRNAs between HPV+ and HPV- patients (Figure 2(b)). Interestingly, in HPV+ patients all miRNAs were downregulated and in HPV- patients all miRNAs were upregulated. The relative expression between HPV+ and HPV- was significant in miR-205-5p, miR-122-5p, miR-124-3p, and miR-146a-5p (p < 0.05). However, miR-92a-3p showed no significant difference (p > 0.05).

Finally, the expression profiles of the five miRNAs were evaluated to discriminate HNSCC patient’s tumor stages (TNM 2-4) (Figure 2(c)). miR-122-5p could differentiate TNM 2 and 3 (p < 0.01) with AUC of 0.92. miR-124-3p was able to differentiate stages 2, 3, and 4 (p < 0.001) and AUC of 0.98. miR-146a-5p differentiates between TNM 2 and 3 (p < 0.001), AUC of 0.97. Finally, miR-92a-3p could differentiate HNSCC patients on TNM 3 (p < 0.001, AUC of 0.99).

4. Discussion

The increased incidence of HNSCC patients contrasts with the lack of early diagnosis biomarkers that reduce cancer-related morbidity and mortality rates [36]. miRNA profiles have been shown as potentially useful biomarkers in cancer diagnosis [36–40]. Different miRNA expression has been identified in several studies of different racial/ethnic populations [3, 19, 34, 36, 37, 41], and the association of miRNA profiles with the population genetic background could increase the understanding of the molecular mechanisms of HNSCC [42]. Several studies worldwide have identified downregulated and upregulated miRNAs as biomarkers for different tumor locations of HNSCC. Kikkawa et al., 2010, identified miRNAs in hypopharyngeal [21]. Hui et al., 2010, identified miRNAs in larynx, oropharynx, and hypopharynx [22]. Liu et al., 2010, identified miRNAs in head and neck [24]. Nohata et al., 2011, identified miRNAs in maxillary sinus [25]. Lajer et al., 2011, identified miRNAs in oral cavity [26]. Severino et al., 2013, identified miRNAs in oral cavity [27]. Fukumoto et al., 2014 and 2015, identified miRNAs in hypopharyngeal [28, 43]. Zhang et al., 2014, identified miRNAs in larynx [29]. Martinez et al., 2015, identified miRNAs in head and neck [30]. Wang et al., 2016, identified miRNAs in nasopharyngeal [31]. Finally, Lubov et al., 2017, have demonstrated that several studies in Caucasian and Asian populations showed elevated expression of miRNAs associated with poor prognosis in HNSCC and decreased expression of miRNAs correlated with lower survival and metastasis [13].

Analysis of ‘liquid biopsies’ is being evaluated as a tool in the care of patients with cancer in order to predict relapse at the earliest stage [6, 44]. For instance, plasma EBV DNA is a biomarker for nasopharyngeal carcinoma [8]. Over the last decade, saliva has emerged as an alternative diagnostic medium to detect both local and systemic events. Saliva medium carries a number of advantages compared with other tissue testing. First, it is produced and collected locally near HNSCC malignancies. Second, saliva sample is noninvasive and easy to perform. Third, the readily available nature of saliva collection also permits multiple samples to be collected for validation and follow-up [45]. The detection of salivary miRNAs is a promising noninvasive method to determine the risk of premalignant head and neck lesions [46]. Park et al., 2009, identified miR-125a and miR-200a as downregulated miRNAs in oral squamous cell carcinoma [19]. Liu et al., 2012, identified miR-31 as upregulated miRNA in oral squamous cell carcinoma [32, 47]. Salazar et al., 2014, developed a panel of saliva-based diagnostic biomarkers for the detection of HNSCC independently validated, using miRNA expression data from The Cancer Genome Atlas (TCGA). They concluded that the saliva-derived miRNAs miR-9, miR-134, and miR-191 may serve as novel biomarkers to reliably detect HNSCC [3].

Regarding our study, we have studied the expression of five different miRNAs (miR-122-5p; miR-92a-3p; miR-124-3p; miR-205-5p; and miR-146a-5p), selected from a miRNA Profiler PCR Arrays, in 108 saliva samples of HNSCC patients and 108 controls, and examined the association with clinicopathological features to determine whether these miRNAs are useful as diagnostic biomarkers.

Previous studies describe miR-92a as a microRNA belonging to the miR-17 ~92 group, which was described as oncogenic. However, subsequent studies have analyzed the miRNAs belonging to this group, indicating that not all act as oncogenes [48]. A possible link between the differentiated expression of miR-92a and tumorigenesis in plasma of patients with different types of malignant neoplasia has been demonstrated. In leukemia this miRNA has been identified as an effective biomarker to distinguish between healthy patients and individuals [49]. Different expression patterns have also been identified in patients with colorectal cancer [50]. In addition, in patients with hepatocellular cancer it was observed that miR-92a had a downregulation compared to healthy controls [51]. A similar result was obtained in the Brazilian population where it was verified that the miR17 ~92 group had different expression (downregulation) in patients with HNSCC, different from our study where this miRNA was found upregulated. However, MacLellan et al., 2012, found an upregulation of miR-92a-3p in serum of HNSCC patients [52]. Moreover, several experiments have found miR-92a acting on the proliferation of cancer cells and tumor growth, inhibiting cell apoptosis and metastasis [53].

On the other hand, it has been shown that miR-124 appears as a tumor suppressor, and the low expression of this miRNA can trigger tumor initiation or progression [54]. Moreover, miR-124 was found reduced in patients with HNSCC, and it was seen that miR-124 binds to the 3′-UTR region of transcription to negatively regulate ESX levels (Epithelial Restriction with Serine) responsible for regulating the epidermal growth factor receptor in HNSCC [55]. However, according to Liu et al., 2010, miR-124 is downregulated in HNSCC patients [24], contrary to Wong et al., 2008, who found that miR-124 was upregulated in HNSCC patients [56], similar to our results. Additionally, Sethi et al., 2014, found that miR-124-3p works as a potential tumor-suppressor miRNA [32].

Regarding miR-122-5p, a study by Ma et al., 2010, showed that overexpression of miR-122 could induce apoptosis and arrest of cancerous cell development through reduced expression of the Bcl-W and/or CCNG1 genes. These findings suggest that miR-122 may behave as a tumor suppressor in different types of cancer [57]. In a previous study, it was also shown that miR-122 plays an important role in the inhibition of tumorigenesis through regulation of the PI3K/Akt/mTOR/p70S6K pathway [58]. However, studies on this particular miRNA are limited, so this finding could define miR-122 as a specific biomarker for the diagnosis of HNSCC, empowering pharmacogenomics and precision medicine in the high altitude Ecuadorian mestizo population [12, 59–64].

miR-146a has been extensively characterized in head and neck cancer from several studies. In thyroid cancer miR-146a has been found to have a 3.9 fold-change lower in patients suffering from the disease compared to healthy individuals [65]. Another study showed that there is a relationship between miR-146a overexpression and progression to metastatic tumors in patients with nasopharyngeal cancer [56]. On the other hand, the involvement of miR-146a in carcinogenesis is demonstrated by its altered expression in several types of tumor, such as thyroid, prostate, pancreatic, gastric, and oral cancer [66]. miR-146a acts by suppressing the expression of the kinase-1 associated with the interleukin-1 receptor, implicated in the NF-κB pathway, which has been shown to participate in human carcinogenesis [67]. Orsos et al., 2013, found that miR-146-5p increases the risk of developing HNSCC when it is upregulated [66], similar to our results.

It has been described that miR-205 is a tissue-specific miRNA highly expressed in the epidermis being essential both in epithelial biogenesis and in the maintenance of the epithelium; a low regulation can lead to alterations, such as the proliferation of cancer cells [68, 69]. In several studies miR205 has been shown to be a promising biomarker in HNSCC, with high sensitivity and specificity [70]. This is also confirmed in a study by Childs et al., 2009, who concluded that miR-205 appears to act as a tumor suppressor in HNSCC, and the reduced expression of this miRNA should result in increased expression of target oncogenes [71]. However, we did not find a statistical difference with this miRNA between patients and controls.

In this study, we have identified a panel of four biomarkers for the diagnosis of HNSCC (hsa-miR-92a-3p, has-miR-124-3p, hsa-miR-122-5p, and hsa-miR-146a-5p) in the high altitude mestizo Ecuadorian population. However, different populations with HNSCC show a differential expression of miRNAs. For example, Li et al., 2009, identify, in Chinese patients with HNSCC, that an increased expression of miR-21 leads to a decrease in the expression of phosphatidylinositol-3,4,5-trisphosphate 3-phosphatase (PTEN). In addition, when miR-31 is deregulated in different types of cancer, including head and neck cancer, the hypoxia-inducible factor (HIF) is activated, promoting tumor angiogenesis [72]. According to Xie et al., 2013, miR-144, miR-21, miR-451, miR-634, miR-486-5p, miR-10b, and several regulatory miRNAs were found in HNSCC patients; miR-10b, miR-144, and miR-451 were found in whole saliva and miR-10b, miR-144, miR-21, and miR-451 in the saliva supernatant. With this study it was possible to establish the presence of miRNAs in all salivary fractions for the detection of esophageal cancer [73].

Moreover, the deregulation of miRNAs in different TNM in HNSCC has been found. According to Chang et al., 2013, both miR-17 and miR-20a are downregulated in advanced TNM in oral squamous cell carcinoma [74]. We found that miR-205 was downregulated in 4 TNM patients, mir-122 upregulated in TNM 2 and 3, miR-124 upregulated in TNM 2 and 3 and downregulated in TNM4, miR-146 upregulated in TNM 2 and 3, and finally miR-92a upregulated in TNM 2 and 3. Furthermore, miR-146, -149, -196, and -499 were found to increase the risk of HNSCC associated to HPV16 [32, 75]. This is similar to our findings, where miR-146 was upregulated in HPV+ patients.

In conclusion, the increased incidence of HNSCC patients contrasts with the lack of early diagnosis biomarkers that could reduce cancer-related morbidity and mortality rates worldwide. miRNA profiles have been shown as potentially useful biomarkers in cancer diagnosis. Our findings have shown that altered expression of miRNAs is different from that one shown in other populations worldwide. That is, hsa-miR-92a-3p, has-miR-124-3p, hsa-miR-122-5p, and hsa-miR-146a-5p miRNAs could be promising biomarkers for HNSCC diagnosis in the high altitude mestizo Ecuadorian population who lives 2,800 meters above sea level.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Carolina Salazar-Ruales designed the study; Carolina Salazar-Ruales and Andrés López-Cortés wrote the manuscript; Jessica-Viviana Arguello and Patricio Peralta collected the samples; Carolina Salazar-Ruales, Jessica-Viviana Arguello, Jennyfer M. García-Cárdenas, Andrés López-Cortés, and Patricia Guevara-Ramírez performed the molecular and data analyses; Paola E. Leone and César Paz-y-Miño supervised the project.

References

Global Burden of Disease Cancer Collaboration, “Global, regional, and national cancer incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 32 cancer groups, 1990 to 2015,” Journal of the American Medical Association (JAMA) Oncology, vol. 3, no. 4, 2017.
View at: Google Scholar
W. Blaszczak, W. Barczak, A. Wegner, W. Golusinski, and W. M. Suchorska, “Clinical value of monoclonal antibodies and tyrosine kinase inhibitors in the treatment of head and neck squamous cell carcinoma,” Medical Oncology, vol. 34, no. 4, 2017.
View at: Google Scholar
C. Salazar, R. Nagadia, P. Pandit et al., “A novel saliva-based microRNA biomarker panel to detect head and neck cancers,” Cellular Oncology, vol. 37, no. 5, pp. 331–338, 2014.
View at: Publisher Site | Google Scholar
R. Lambert, C. Sauvaget, M. De Camargo Cancela, and R. Sankaranarayanan, “Epidemiology of cancer from the oral cavity and oropharynx,” European Journal of Gastroenterology & Hepatology, vol. 23, no. 8, pp. 633–641, 2011.
View at: Publisher Site | Google Scholar
A. K. Chaturvedi, W. F. Anderson, J. Lortet-Tieulent et al., “Worldwide trends in incidence rates for oral cavity and oropharyngeal cancers,” Journal of Clinical Oncology, vol. 31, no. 36, pp. 4550–4559, 2013.
View at: Publisher Site | Google Scholar
C. R. Leemans, P. J. Snijders, and R. H. Brakenhoff, “The molecular landscape of head and neck cancer,” Nature Reviews Cancer, vol. 18, no. 5, pp. 269–282, 2018.
View at: Google Scholar
V. A. Moyer, “Screening for prostate cancer: U.S. preventive services task force recommendation statement,” Annals of Internal Medicine, vol. 157, no. 8, pp. 567–570, 2012.
View at: Publisher Site | Google Scholar
K. A. Chan, J. K. Woo, and A. King, “Analysis of Plasma Epstein–Barr Virus DNA to Screen for Nasopharyngeal Cancer,” The New England Journal of Medicine, vol. 377, no. 6, pp. 513–522, 2017.
View at: Publisher Site | Google Scholar
F. Bray, J.-S. Ren, E. Masuyer, and J. Ferlay, “Global estimates of cancer prevalence for 27 sites in the adult population in 2008,” International Journal of Cancer, vol. 132, no. 5, pp. 1133–1145, 2013.
View at: Publisher Site | Google Scholar
J. Ferlay, I. Soerjomataram, R. Dikshit, and S. Eser, “Mortality World Wide: IARC CancerBase, E. M. Cancer Incidence and Mortality Worldwide: IARC CancerBase No. 11,” Globocan 2012 v1.0, 2013, accessed on 12/08/2014, Available from: http://globocan.iarc.fr.
View at: Google Scholar
S. Núñez-González, J. A. Delgado-Ron, C. Gault, and D. Simancas-Racines, “Trends and spatial patterns of oral cancer mortality in Ecuador, 2001–2016,” International Journal of Dentistry, vol. 2018, Article ID 6086595, 7 pages, 2018.
View at: Publisher Site | Google Scholar
S. Guerrero, A. López-Cortés, A. Indacochea et al., “Analysis of racial/ethnic representation in select basic and applied cancer research studies,” Scientific Reports, vol. 8, no. 1, 2018.
View at: Publisher Site | Google Scholar
J. Lubov, M. Maschietto, I. Ibrahim et al., “Meta-analysis of micrornas expression in head and neck cancer: Uncovering association with outcome and mechanisms,” Oncotarget , vol. 8, no. 33, pp. 55511–55524, 2017.
View at: Google Scholar
J. M. Yoshizawa and D. T. W. Wong, “Salivary microRNAs and oral cancer detection,” Methods in Molecular Biology, vol. 936, pp. 313–324, 2013.
View at: Publisher Site | Google Scholar
A. López-Cortés, S. Guerrero, M. A. Redal, A. T. Alvarado, and L. A. Quiñones, “State of art of cancer pharmacogenomics in Latin American populations,” International Journal of Molecular Sciences, vol. 18, no. 6, 2017.
View at: Google Scholar
C. Punyadeera and P. D. Slowey, “Saliva as an emerging biofluid for clinical diagnosis and applications of MEMS/NEMS in salivary diagnostics,” Nanobiomaterials in Clinical Dentistry, pp. 453–473, 2012.
View at: Google Scholar
B. L. Schulz, J. Cooper-White, and C. K. Punyadeera, “Saliva proteome research: current status and future outlook,” Critical Reviews in Biotechnology, vol. 33, no. 3, pp. 246–259, 2013.
View at: Publisher Site | Google Scholar
Y. Li, X. Zhou, M. A. R. St. John, and D. T. W. Wong, “RNA profiling of cell-free saliva using microarray technology,” Journal of Dental Research, vol. 83, no. 3, pp. 199–203, 2004.
View at: Publisher Site | Google Scholar
N. J. Park, H. Zhou, D. Elashoff et al., “Salivary microRNA: discovery, characterization, and clinical utility for oral cancer detection,” Clinical Cancer Research, vol. 15, no. 17, pp. 5473–5477, 2009.
View at: Publisher Site | Google Scholar
R. S. Patel, A. Jakymiw, B. Yao et al., “High resolution of microRNA signatures in human whole saliva,” Archives of Oral Biolog, vol. 56, no. 12, pp. 1506–1513, 2011.
View at: Publisher Site | Google Scholar
N. Kikkawa, T. Hanazawa, L. Fujimura et al., “MiR-489 is a tumour-suppressive miRNA target PTPN11 in hypopharyngeal squamous cell carcinoma (HSCC),” British Journal of Cancer, vol. 103, no. 6, pp. 877–884, 2010.
View at: Publisher Site | Google Scholar
A. B. Y. Hui, M. Lenarduzzi, T. Krushel et al., “Comprehensive MicroRNA profiling for head and neck squamous cell carcinomas,” Clinical Cancer Research, vol. 16, no. 4, pp. 1129–1139, 2010.
View at: Publisher Site | Google Scholar
K. Koshizuka, T. Hanazawa, I. Fukumoto, N. Kikkawa, Y. Okamoto, and N. Seki, “The microRNA signatures: Aberrantly expressed microRNAs in head and neck squamous cell carcinoma,” Journal of Human Genetics, vol. 62, no. 1, pp. 3–13, 2017.
View at: Publisher Site | Google Scholar
C. J. Liu, M. M. Tsai, P. S. Hung et al., “miR-31 ablates expression of the HIF regulatory factor FIH to activate the HIF pathway in head and neck carcinoma,” Cancer Research, vol. 70, no. 4, pp. 1635–1644, 2010.
View at: Publisher Site | Google Scholar
N. Nohata, T. Hanazawa, N. Kikkawa et al., “Tumour suppressive microRNA-874 regulates novel cancer networks in maxillary sinus squamous cell carcinoma,” British Journal of Cancer, vol. 105, no. 6, pp. 833–841, 2011.
View at: Publisher Site | Google Scholar
C. B. Lajer, F. C. Nielsen, L. Friis-Hansen et al., “Different miRNA signatures of oral and pharyngeal squamous cell carcinomas: a prospective translational study,” British Journal of Cancer, vol. 104, no. 5, pp. 830–840, 2011.
View at: Publisher Site | Google Scholar
P. Severino, H. Brüggemann, F. M. Andreghetto et al., “MicroRNA expression profile in head and neck cancer: HOX-cluster embedded microRNA-196a and microRNA-10b dysregulation implicated in cell proliferation,” BMC Cancer, vol. 13, pp. 1–15, 2013.
View at: Google Scholar
I. Fukumoto, T. Kinoshita, T. Hanazawa et al., “Identification of tumour suppressive microRNA-451a in hypopharyngeal squamous cell carcinoma based on microRNA expression signature,” British Journal of Cancer, vol. 111, no. 2, pp. 386–394, 2014.
View at: Publisher Site | Google Scholar
Y. Zhang, Y. Chen, J. Yu, G. Liu, and Z. Huang, “Integrated transcriptome analysis reveals miRNA-mRNA crosstalk in laryngeal squamous cell carcinoma,” Genomics, vol. 104, no. 4, pp. 249–256, 2014.
View at: Publisher Site | Google Scholar
B. V. Martinez, J. M. Dhahbi, Y. O. Nunez Lopez et al., “Circulating small non coding RNA signature in head and neck squamous cell carcinoma,” Oncotarget , vol. 6, no. 22, pp. 19246–19263, 2015.
View at: Google Scholar
F. Wang, J. Lu, X. Peng et al., “Integrated analysis of microRNA regulatory network in nasopharyngeal carcinoma with deep sequencing,” Journal of Experimental & Clinical Cancer Research, vol. 35, no. 1, article 17, 2016.
View at: Publisher Site | Google Scholar
N. Sethi, A. Wright, H. Wood, and P. Rabbitts, “MicroRNAs and head and neck cancer: reviewing the first decade of research,” European Journal of Cancer, vol. 50, no. 15, pp. 2619–2635, 2014.
View at: Publisher Site | Google Scholar
H. Polanska, M. Raudenska, J. Gumulec et al., “Clinical significance of head and neck squamous cell cancer biomarkers,” Oral Oncology, vol. 50, no. 3, pp. 168–177, 2014.
View at: Publisher Site | Google Scholar
F. Zahran, D. Ghalwash, O. Shaker, K. Al-Johani, and C. Scully, “Salivary microRNAs in oral cancer,” Oral Diseases, vol. 21, no. 6, pp. 739–747, 2015.
View at: Publisher Site | Google Scholar
K. J. Livak and T. D. Schmittgen, “Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method,” Methods, vol. 25, no. 4, pp. 402–408, 2001.
View at: Google Scholar
M. Avissar, M. D. McClean, K. T. Kelsey, and C. J. Marsit, “MicroRNA expression in head and neck cancer associates with alcohol consumption and survival,” Carcinogenesis, vol. 30, no. 12, pp. 2059–2063, 2009.
View at: Publisher Site | Google Scholar
G. A. Calin, M. Ferracin, A. Cimmino et al., “A microRNA signature associated with prognosis and progression in chronic lymphocytic leukemia,” The New England Journal of Medicine, vol. 353, no. 17, pp. 1793–1801, 2005.
View at: Publisher Site | Google Scholar
Y. Chen and R. L. Stallings, “Differential patterns of microRNA expression in neuroblastoma are correlated with prognosis, differentiation, and apoptosis,” Cancer Research, vol. 67, no. 3, pp. 976–983, 2007.
View at: Publisher Site | Google Scholar
O. Slaby, M. Svoboda, and P. Fabian, “Altered expression of miR-21, miR-31, miR-143 and miR-145 is related to clinicopathologic features of colorectal cancer,” Oncology, vol. 72, no. 5-6, pp. 397–402, 2008.
View at: Publisher Site | Google Scholar
G. Heller, M. Weinzierl, C. Noll et al., “Genome-wide miRNA expression profiling identifies miR-9-3 and miR-193a as targets for DNA methylation in non-small cell lung cancers,” Clinical Cancer Research, vol. 18, no. 6, pp. 1619–1629, 2012.
View at: Publisher Site | Google Scholar
R. S. Huang, E. R. Gamazon, D. Ziliak et al., “Population differences in microRNA expression and biological implications,” RNA Biology, vol. 8, no. 4, pp. 692–701, 2011.
View at: Publisher Site | Google Scholar
S. Irani, “miRNAs signature in head and neck squamous cell carcinoma metastasis: a literature review,” Journal of Dentistry, vol. 17, no. 172, pp. 71–83, 2016.
View at: Google Scholar
I. Fukumoto, T. Hanazawa, T. Kinoshita et al., “MicroRNA expression signature of oral squamous cell carcinoma: functional role of microRNA-26a/b in the modulation of novel cancer pathways,” British Journal of Cancer, vol. 112, no. 5, pp. 891–900, 2015.
View at: Publisher Site | Google Scholar
X. Cai, F. Janku, Q. Zhan, and J.-B. Fan, “Accessing genetic information with liquid biopsies,” Trends in Genetics, vol. 31, no. 10, pp. 564–575, 2015.
View at: Publisher Site | Google Scholar
D. A. Ovchinnikov, M. A. Cooper, P. Pandit et al., “Tumor-suppressor gene promoter hypermethylation in saliva of head and neck cancer patients,” Translational Oncology, vol. 5, no. 5, pp. 321–326, 2012.
View at: Publisher Site | Google Scholar
Y. Yang, Y.-X. Li, X. Yang, L. Jiang, Z.-J. Zhou, and Y.-Q. Zhu, “Progress risk assessment of oral premalignant lesions with saliva miRNA analysis,” BMC Cancer, vol. 13, article 129, 2013.
View at: Publisher Site | Google Scholar
C.-J. Liu, S.-C. Lin, C.-C. Yang, H.-W. Cheng, and K.-W. Chang, “Exploiting salivary miR-31 as a clinical biomarker of oral squamous cell carcinoma,” Head & Neck, vol. 34, no. 2, pp. 219–224, 2012.
View at: Publisher Site | Google Scholar
E. Mogilyansky and I. Rigoutsos, “The miR-17/92 cluster: a comprehensive update on its genomics, genetics, functions and increasingly important and numerous roles in health and disease,” Cell Death & Differentiation, vol. 20, no. 12, pp. 1603–1614, 2013.
View at: Publisher Site | Google Scholar
M. Tanaka, K. Oikawa, M. Takanashi et al., “Down-regulation of miR-92 in human plasma is a novel marker for acute leukemia patients,” PLoS ONE, vol. 4, no. 5, Article ID e5532, 2009.
View at: Publisher Site | Google Scholar
Z. Huang, D. Huang, S. Ni, Z. Peng, W. Sheng, and X. Du, “Plasma microRNAs are promising novel biomarkers for early detection of colorectal cancer,” International Journal of Cancer, vol. 127, no. 1, pp. 118–126, 2010.
View at: Publisher Site | Google Scholar
M. Shigoka, A. Tsuchida, T. Matsudo et al., “Deregulation of miR-92a expression is implicated in hepatocellular carcinoma development,” Pathology International, vol. 60, no. 5, pp. 351–357, 2010.
View at: Publisher Site | Google Scholar
S. A. MacLellan, J. Lawson, J. Baik, M. Guillaud, C. F. Poh, and C. Garnis, “Differential expression of miRNAs in the serum of patients with high-risk oral lesions,” Cancer Medicine, vol. 1, no. 2, pp. 268–274, 2012.
View at: Publisher Site | Google Scholar
M. Li, X. Guan, Y. Sun et al., “MiR-92a family and their target genes in tumorigenesis and metastasis,” Experimental Cell Research, vol. 323, no. 1, pp. 1–6, 2014.
View at: Publisher Site | Google Scholar
Y. Zhao, Z. Ling, Y. Hao et al., “MiR-124 acts as a tumor suppressor by inhibiting the expression of sphingosine kinase 1 and its downstream signaling in head and neck squamous cell carcinoma,” Oncotarget , vol. 8, no. 15, pp. 25005–25020, 2017.
View at: Google Scholar
P. Wang, L. Chen, J. Zhang et al., “Methylation-mediated silencing of the miR-124 genes facilitates pancreatic cancer progression and metastasis by targeting Rac1,” Oncogene, vol. 33, no. 4, pp. 514–524, 2014.
View at: Publisher Site | Google Scholar
T.-S. Wong, X.-B. Liu, B. Y.-H. Wong, R. W.-M. Ng, A. P.-W. Yuen, and W. I. Wei, “Mature miR-184 as potential oncogenic microRNA of squamous cell carcinoma of tongue,” Clinical Cancer Research, vol. 14, no. 9, pp. 2588–2592, 2008.
View at: Publisher Site | Google Scholar
L. Ma, J. Liu, J. Shen et al., “Expression of miR-122 mediated by adenoviral vector induces apoptosis and cell cycle arrest of cancer cells,” Cancer Biology & Therapy, vol. 9, no. 7, pp. 554–561, 2014.
View at: Publisher Site | Google Scholar
B. Qian, S. A. Nag, Y. Su et al., “miRNAs in cancer prevention and treatment and as molecular targets for natural product anticancer agents,” Current Cancer Drug Targets, vol. 13, no. 5, pp. 519–541, 2013.
View at: Publisher Site | Google Scholar
A. López-Cortés, S. Guerrero, M. Redal, A. Alvarado, and L. Quiñones, “State of Art of Cancer Pharmacogenomics in Latin American Populations,” International Journal of Molecular Sciences, vol. 18, no. 6, p. 639, 2017.
View at: Publisher Site | Google Scholar
L. Quinones, M. Lavanderos, J. Cayun et al., “Perception of the usefulness of drug/gene Pairs and barriers for pharmacogenomics in Latin America,” Current Drug Metabolism, vol. 15, no. 2, pp. 202–208, 2014.
View at: Publisher Site | Google Scholar
A. López-Cortés, G. Jaramillo-Koupermann, and M. J. Muñoz, “Genetic polymorphisms in MTHFR (C677T, A1298C), MTR (A2756G) and MTRR (A66G) genes associated with pathological characteristics of prostate cancer in the ecuadorian population,” The American Journal of the Medical Sciences, vol. 346, no. 6, pp. 447–454, 2013.
View at: Publisher Site | Google Scholar
C. Paz-y-Miño, M. J. Muñoz, A. López-Cortés et al., “Frequency of polymorphisms pro198leu in GPX-1 gene and ile58thr in MnSOD gene in the altitude ecuadorian population with bladder cancer,” Oncology Research : Featuring Preclinical and Clinical Cancer Therapeutics, vol. 18, no. 8, pp. 395–400, 2010.
View at: Publisher Site | Google Scholar
A. López-Cortés, A. Cabrera-Andrade, F. Oña-Cisneros et al., “Breast cancer risk associated with genotype polymorphisms of the aurora kinase a gene (aurka): a case-control study in a high altitude ecuadorian mestizo population,” Pathology & Oncology Research, pp. 1–9, 2017.
View at: Google Scholar
A. López-Cortés, P. E. Leone, B. Freire-Paspuel et al., “Mutational Analysis of Oncogenic AKT1 Gene Associated with Breast Cancer Risk in the High Altitude Ecuadorian Mestizo Population,” BioMed Research International, vol. 2018, Article ID 7463832, 10 pages, 2018.
View at: Publisher Site | Google Scholar
K. Jazdzewski, E. L. Murray, K. Franssila, B. Jarzab, D. R. Schoenberg, and A. De La Chapelle, “Common SNP in pre-miR-146a decreases mature miR expression and predisposes to papillary thyroid carcinoma,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 105, no. 20, pp. 7269–7274, 2008.
View at: Publisher Site | Google Scholar
Z. Orsós, I. Szanyi, A. Csejtei, I. Gerlinger, I. Ember, and I. Kiss, “Association of pre-miR-146a rs2910164 polymorphism with the risk of head and neck cancer,” Anticancer Reseach, vol. 33, no. 1, pp. 341–346, 2013.
View at: Google Scholar
M. Masroor, R. Mir, J. Javid et al., “Cell free EGFR mRNA expression and implications for survival and metastasis in non-small cell lung cancer cases,” Asian Pacific Journal of Cancer Prevention, vol. 16, no. 15, pp. 6445–6449, 2015.
View at: Publisher Site | Google Scholar
M. Teta, Y. S. Choi, T. Okegbe et al., “Inducible deletion of epidermal Dicer and Drosha reveals multiple functions for miRNAs in postnatal skin,” Development, vol. 139, no. 8, pp. 1405–1416, 2012.
View at: Publisher Site | Google Scholar
L. F. Sempere, M. Christensen, A. Silahtaroglu et al., “Altered microRNA expression confined to specific epithelial cell subpopulations in breast cancer,” Cancer Research, vol. 67, no. 24, pp. 11612–11620, 2007.
View at: Publisher Site | Google Scholar
N. Tran, T. McLean, X. Zhang et al., “MicroRNA expression profiles in head and neck cancer cell lines,” Biochemical and Biophysical Research Communications, vol. 358, no. 1, pp. 12–17, 2007.
View at: Publisher Site | Google Scholar
G. Childs, M. Fazzari, G. Kung et al., “Low-level expression of microRNAs let-7d and miR-205 are prognostic markers of head and neck squamous cell carcinoma,” The American Journal of Pathology, vol. 174, no. 3, pp. 736–745, 2009.
View at: Publisher Site | Google Scholar
J. Li, H. Huang, L. Sun et al., “MiR-21 indicates poor prognosis in tongue squamous cell carcinomas as an apoptosis inhibitor,” Clinical Cancer Research, vol. 15, no. 12, pp. 3998–4008, 2009.
View at: Publisher Site | Google Scholar
Z. Xie, G. Chen, X. Zhang et al., “Salivary microRNAs as promising biomarkers for detection of esophageal cancer,” PLoS ONE, vol. 8, no. 4, Article ID e57502, 2013.
View at: Publisher Site | Google Scholar
C.-C. Chang, Y.-J. Yang, Y.-J. Li et al., “MicroRNA-17/20a functions to inhibit cell migration and can be used a prognostic marker in oral squamous cell carcinoma,” Oral Oncology, vol. 49, no. 9, pp. 923–931, 2013.
View at: Publisher Site | Google Scholar
C. Salazar, D. Calvopiña, and C. Punyadeera, “MiRNAs in human papilloma virus associated oral and oropharyngeal squamous cell carcinomas,” Expert Review of Molecular Diagnostics, vol. 14, no. 8, pp. 1033–1040, 2014.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2018 Carolina Salazar-Ruales et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

3440

Downloads

1918

Citations

BioMed Research International

Salivary MicroRNAs for Early Detection of Head and Neck Squamous Cell Carcinoma: A Case-Control Study in the High Altitude Mestizo Ecuadorian Population

Abstract

1. Introduction

2. Materials and Methods

2.1. Study Subjects

2.2. Total RNA and miRNA Extraction

2.3. Total RNA and miRNA cDNA Synthesis

2.4. miRNA Profiler PCR Arrays

2.5. Relative miRNA Expression and the Comparative Livak Method (2-ΔΔCT)

2.6. miRNA Quantitative Real-Time PCR

2.7. Statistical Analysis

3. Results

4. Discussion

Data Availability

Conflicts of Interest

Authors’ Contributions

References

Copyright

2.5. Relative miRNA Expression and the Comparative Livak Method (2^-ΔΔCT)