Research article | Open | Open Peer Review | Published:
MicroRNA expression profile in head and neck cancer: HOX-cluster embedded microRNA-196a and microRNA-10b dysregulation implicated in cell proliferation
BMC Cancervolume 13, Article number: 533 (2013)
Current evidence implicates aberrant microRNA expression patterns in human malignancies; measurement of microRNA expression may have diagnostic and prognostic applications. Roles for microRNAs in head and neck squamous cell carcinomas (HNSCC) are largely unknown. HNSCC, a smoking-related cancer, is one of the most common malignancies worldwide but reliable diagnostic and prognostic markers have not been discovered so far. Some studies have evaluated the potential use of microRNA as biomarkers with clinical application in HNSCC.
MicroRNA expression profile of oral squamous cell carcinoma samples was determined by means of DNA microarrays. We also performed gain-of-function assays for two differentially expressed microRNA using two squamous cell carcinoma cell lines and normal oral keratinocytes. The effect of the over-expression of these molecules was evaluated by means of global gene expression profiling and cell proliferation assessment.
Altered microRNA expression was detected for a total of 72 microRNAs. Among these we found well studied molecules, such as the miR-17-92 cluster, comprising potent oncogenic microRNA, and miR-34, recently found to interact with p53. HOX-cluster embedded miR-196a/b and miR-10b were up- and down-regulated, respectively, in tumor samples. Since validated HOX gene targets for these microRNAs are not consistently deregulated in HNSCC, we performed gain-of-function experiments, in an attempt to outline their possible role. Our results suggest that both molecules interfere in cell proliferation through distinct processes, possibly targeting a small set of genes involved in cell cycle progression.
Functional data on miRNAs in HNSCC is still scarce. Our data corroborate current literature and brings new insights into the role of microRNAs in HNSCC. We also show that miR-196a and miR-10b, not previously associated with HNSCC, may play an oncogenic role in this disease through the deregulation of cell proliferation. The study of microRNA alterations in HNSCC is an essential step to the mechanistic understanding of tumor formation and could lead to the discovery of clinically relevant biomarkers.
MicroRNAs (miRNAs) are ~22 nt non-coding RNA molecules that negatively regulate gene expression by degrading or destabilizing the messenger RNA (mRNA) or by inhibiting protein translation ; some reports demonstrate that they may also function as positive regulators [2, 3]. MiRNAs have been shown to contribute to cancer development and progression, and are differentially expressed between normal tissues and cancers . Although the function of most of the miRNAs identified to date has yet to be determined, their use as potential biomarkers or therapeutic targets has been considered in several human diseases and cancers [5, 6].
Head and neck squamous cell carcinoma (HNSCC) is a significant public health entity, representing the sixth leading cancer by incidence worldwide [7, 8]. Genetic changes that bring about HNSCC are usually a consequence of continued exposure to carcinogens associated with tobacco. Despite advances in medical and surgical treatment, the overall 5-year survival rate for patients with HNSCC remains around 50% . A recent work by Liu et al., 2009  analyzed data compiled by the American Cancer Society and points out that new cases of HNSCC increased 25% during the past 5 years, highlighting the need for a better understanding of the molecular events leading to the development of this disease.
The number of studies addressing the contribution of miRNA deregulation in the context of HNSCC is, however, limited [10, 11]. Some of these studies have evaluated the potential use of miRNAs as biomarkers with clinical application, associating the expression levels of some of these miRNAs with survival rates or metastatic potential [12–16]. Overall, results are promising, but still preliminary and lacking corroboration.
In our study we determined the miRNA expression profile of oral squamous cell carcinoma (OSCC) samples, a type of cancer that represents 90% of all HNSCC . We also performed functional assays for two differentially expressed miRNAs, miR-196 and miR-10b, since neither have been previously associated with HNSCC. Despite studies linking these miRNAs to the regulation of homeobox (HOX) genes [17, 18] our data suggest that, in the case of HNSCC, they have little impact in HOX gene expression. We demonstrate that miR-10b and miR-196a interfere in cell proliferation through distinct processes and in a cell-type dependent manner.
Fifteen patients with OSCC (tongue and floor of the mouth) were selected for the microarray experiments. In order to validate the microarray results, 35 additional patients with HNSCC (oral cavity, oropharynx and larynx) were selected. The clinical and pathological profile of patients is shown in Table 1. The average age of patients was 55.5 years (SD 9.8, range 38–82 years), and the male/female ratio was 24:1. Most patients were smokers or former smokers and had a history of chronic alcohol abuse. Tumor and corresponding cancer free surgical margins containing the corresponding epithelium were collected from patients submitted to surgical resection of primary tumor at Hospital das Clinicas, Hospital Heliopolis and Arnaldo Vieira de Carvalho Cancer Institute, in Sao Paulo, Brazil. All patients provided written informed consent, and the research protocol was approved by review boards of all institutions involved and by the National Committee of Ethics in Research (CONEP 1763/05). Samples corresponding to the oral cavity, base of the tongue and larynx were snap-frozen in liquid nitrogen immediately after surgery and stored in liquid nitrogen until RNA preparation. Frozen samples were sectioned using a cryostat, and tissue sections were stained with RNAse-free reagents. Analysis of hematoxylin and eosin-stained sections by the study pathologists confirmed >75% tumor cells in all HNSCC samples and that surgical margins were tumor-free. The diagnosis of HNSCC was confirmed, and tumors were histologically examined for perineural invasion (tumor cells in the perineural space or epineurium), tumor differentiation (well, moderated or poorly differentiated, according to the WHO guidelines), lymphatic-vascular invasion, surgical margins, and lymph node metastasis. Tumors were staged according to the TNM clinical staging system, as proposed by the International Union Against Cancer.
Total RNA was prepared from tissue using mirVana miRNA Isolation Kit (Ambion, Austin, TX) in compliance with the manufacturer’s protocol. RNA integrity and concentration were assessed using the RNA 6000 Nano Assay kit with Agilent 2100 Bioanalyzer according to the manufacturer’s instructions (Agilent Technologies, Palo Alto, CA).
miRNA microarray expression profiling
MiRNA expression profiling was performed using the Illumina miRNA arrays version 1.0. Sample preparation and hybridization followed the manufacturer’s instructions. Briefly, 200 ng of total RNA was first polyadenylated and converted to biotinylated cDNA, which was attached to a solid phase and hybridised with a pool of miRNA-specific oligonucleotides (MSO). Each single MSO is used to assay one miRNA on the panel. Universal PCR amplification was then performed, creating fluorescently labeled products identifiable by their unique MSO sequence. These products were hybridized on the Illumina miRNA array, with the address sequence from each MSO enabling the hybridization of specific miRNA products to specific locations on the BeadArray substrate. Hybridization signals were detected and quantified using Illumina scanner and BeadStudio version 3.2.7. Average signals were quantile normalized and then filtered according their detection p-value: only miRNAs for which the detection p-value was consistently equal or lower than 0.05 across at least one group of samples (marginal and/or tumor samples) were considered expressed and further analyzed. All data is MIAME compliant and the raw data has been deposited in a MIAME compliant database (Gene Expression Omnibus database) under accession Number GSE31277. Differentially expressed miRNAs were determined by using the Rank Product, non-parametric statistical method based on ranks of fold-changes . We used the RankProd package on R for this analysis; the percentage of false positives was calculated and p-values were accordingly corrected for multiple comparisons.
Relative quantification of miRNA levels using real time-PCR
To validate the microarray expression data, miRNAs were subjected to quantitative Real Time-PCR using the TaqMan miRNA assay system (Applied Biosystems, Foster City, CA). Briefly, about 100 ng of total RNA was subjected to a reverse transcription reaction using miRNA-specific looped primers according to the manufacturer’s protocol to obtain the cDNA. Subsequent PCRs used miRNA specific forward and reverse primers along with appropriate cDNA product and TaqMan universal mix. PCR was carried out in AB7500 (Applied Biosystems, Foster City, CA) in a 20 ul volume reaction following thermal cycling parameters suggested by the manufacturer: 50°C for 2 min, 95°C for 10 min and 45 cycles of 95°C for 15 s and 60°C for 1 min.
The expression data was normalized to the RNU48 expression. RNU48 was chosen as a suitable endogenous control gene following analysis of gene expression stability of three candidate genes across our samples. For a detailed description of this step refer to the next Methods’ section. Expression levels were determined using the comparative ∆Ct method .
For miRNAs individually studied in independent sets of samples by quantitative real-time PCR, the nonparametric test Wilcoxon Signed Rank Test was used to detect the statistically significant differences between paired normal tissue (N) and tumor (T) samples obtained from the same individual. This test was performed using SPSS for Windows® Software. The same software was used to calculate the mean and standard deviation of all variables.
Identification of suitable endogenous control gene for microRNA gene expression analysis by real-time PCR
The expression of three snoRNAs (RNU6B, RNU44 and RNU48) was measured by quantitative real-time PCR with TaqMan miRNA assays (Applied Biosystems, Foster City, CA), as previously described for all samples assayed by miRNA microarrays. This data was analyzed using the SLqPCR package in R  to determine the expression stability of these snoRNAs across samples. The stability factor M was calculated for each snoRNA (M (for RNU48) = 0.69; M (for RNU44) = 0.78; M (for RNU6B) = 0.75). Since high expression stability is associated to low M values, RNU48 appeared to be the snoRNA with most stable expression across the set of samples analyzed, hence was chosen as control for normalisation.
Prediction of miRNA targets and their functional analysis
Potential miRNA targets were identified using Ingenuity Pathway Analysis (IPA Ingenuity Systems, http://www.ingenuity.com). Only experimentally validated targets were selected, using miRecords (http://mirecords.biolead.org/), Tarbase (http://microrna.gr/tarbase) or TargetScan (http://www.targetscan.org/). For fuctional annotation of potential targets we used KEGG pathways term enrichment analysis using the computational tool Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.7 (http://david.abcc.ncifcrf.gov/home.jsp).
HNSCC cell line and keratinocyte cell culture
The HNSCC cell lines SCC25 and SCC9, derived from a SCC of the tongue, and FaDu, derived from a SCC of the hypopharynx were used in this study. They were obtained from American Type Culture Collection (SCC25 catalog number CRL-1628, SCC9 catalog number CRL-1629, and FaDu catalog number HTB-43). The cell lines were grown in a Dulbecco’s Modified Eagle’s medium/Nutrient Mixture F-12 Ham (DMEM/F12) supplemented with 10% fetal bovine serum in a humidified atmosphere of 5% CO2 and 95% air at 37°C. Oral keratinocytes were obtained from primary cultures of the buccal mucosa, from voluntary donor patients undergoing surgery performed in out-patient clinics in the Dentistry School of USP. The patients were informed and signed the required Informed Consent. This study was approved by the Research Ethics Committee of the Instituto de Pesquisas Energéticas e Nucleares (IPEN/CNEN-SP) [Institute of Energy and Nuclear Research] (approval number 087/CEP-IPEN/SP). Keratinocytes were plated on a support layer, called feeder-layer, composed of murine fibroblasts of the type 3T3-Swiss albino (ATCC, catalog number CCL-92), which were irradiated (60 Gy), and maintained in an incubator at 37°C, in a humidified atmosphere containing 5% CO2 and grown as previously described .
Transfection of cultured cells for up-regulation of miRNAs
The siPORT NeoFx reagent (Ambion) was used for transfection following the manufacturer’s protocol. For up-regulation, the Ambion Pre-miR™ miRNA Precursor Molecule (hsa-miR-10b and hsa-miR-196a) was used, with Ambion’s Pre-miR negative control #1. Successful up-regulation was achieved with 50 nM of final Pre-miR miRNA Precursor concentration.
Immunofluorescence assay for proliferation analysis
Normal keratinocytes transfected with the miRNA precursor and the negative control were cultured in Lab-Tek Chamber Slides (Nalge Nunc International, Rochester, NY, USA) for the immunofluorescence assay. Cells were fixed with methanol, blocked with 3% bovine serum in PBS, and incubated for 1 h with antihuman Ki67 (Monoclonal Mouse, clone MIB-1, DAKO, Denmark A/5), diluted 1:400. Cells were washed with PBS and incubated at room temperature for 45 minutes with secondary antibody conjugated with fluorescein (1:50) (Novocastra Laboratories, UK), in a dark chamber. Following washing, chambers containing the cells were mounted with VECTASHIELD Mounting Medium with DAPI (Vector Laboratories, Ind. Buelingame, CA 94010). Results were analyzed by fluorescence microscopy (Zeiss Axiophot II, Carl Zeiss, Oberköchen, Germany). The percentage of cells displaying Ki67 labeling was determined by counting the number of positive Ki67 stained cells as a proportion of the total number of cells counted. Cells were counted manually in the whole chamber area.
Proliferation assay by flow cytometry
Cell lines SCC9, SCC25 and FaDu were stained with Cell Trace Violet (Molecular Probes®), according to the manufacturer protocol. Briefly, the cells were incubated with 5 μM Cell Trace Violet for 20 minutes at 37°C, washed twice with fresh and warmed medium and cultured under regular conditions. The cells were run on BD LSR Fortessa flow cytometer with 405 nm laser at day zero and after 72 hours of cell culture for cell proliferation rate assessment. Proliferation rate was determined by fluorescence decay. Analysis was performed using Flow Jo software (Tree Star™). For cell proliferation rates after transfection, cell lines SCC25 and FaDu were stained 24 h after transfection (at the time of medium exchange). Proliferation rates were compared between scramble (negative control) and cells overexpressing miR-10b.
mRNA microarray expression profiling and analysis
Following the transfection assays, the global gene expression analysis was conducted using the Agilent Human Whole Genome Oligonucleotide Microarray (44K; Agilent Technologies) following the manufacturer’s protocols. Oligonucleotide microarrays were scanned using the GenePix 4000B Microarray Scanner (Molecular Devices) and features were automatically extracted and analyzed for quality control using Agilent Feature Extraction Software. Raw data was deposited in a MIAME compliant database (Gene Expression Omnibus database) under the accession Number GSE31277. Partek Genomics Suite 6.6 (Partek Incorporated) was used for normalization of gene expression levels and for fold-change in gene expression calculation. To gain insights into the potential mechanisms affected by the overexpression of the miR-10b and miR-196a in cells, deregulated genes were mapped to regulatory networks using Ingenuity Pathway Analysis (IPA Ingenuity Systems, http://www.ingenuity.com).
Western blotting was performed using a specific antibody against annexin 1 (1:1000 dilution) (Zymed Laboratories - 713400), and β-Actin (1:12.000 dilution) (Cell Signaling Technology, Danvers, MA, USA). Briefly, 72 hours after transfection cells were lysed in RIPA buffer (150 mM NaCl, 10 mM Tris/HCl, pH 7.4, 0.5% Triton X-100 and protease and phosphatase inhibitors (Biogene). Protein concentration was estimated using the BCA Protein Assay Kit (BioAgency, London, UK). 20 ug of protein lysate was separated in 15% SDS gel (GE Healthcare, Piscataway, NJ, USA) and subsequently transferred to nitrocellulose membrane of 0,45 μm (GE Healthcare, Piscataway, NJ, USA). The membranes were blocked using 3% non-fat dry milk, and incubated with primary antibodies overnight at 4°C. The membranes were washed in 1x TBS eith 0.1% Tween-20 (TBS/T), incubated for 1 h with anti-rabbit secondary antibodies conjugated to horseradish peroxidase (Abcam - ab102779) and visualized with a chemiluminescence reagent (ECL) system (Amersham Biosciences, Arlington Heights, IL).
Results and discussion
MiRNA deregulation in OSCC samples: implication in tumor progression
HNSCC can involve multiple anatomical sites, each with individual molecular characteristics, and highly affected by the drinking and smoking habits of patients [13, 23, 24]. In an attempt to limit data variability due to HNSCC subsites and environmental factors, we assessed miRNA expression levels in 15 OSCC samples limited to tongue and floor of the mouth, from patients possessing similar demographic and clinico-pathological characteristics (Table 1, detailed in Methods). Samples were paired with tumor-free surgical margins. The expression profiles of tumor samples revealed significant differential expression for 72 miRNAs compared to their corresponding tumor-free margins (Table 2). Several studies have analysed the miRNA expression profile of OSCC cell lines and tumor samples, with little overlap among results [25, 26]. This inconsistency in results justifies additional studies.
In order to access biological processes possibly targeted by deregulated miRNAs we performed a functional analysis of validated targets through KEGG term enrichment analysis using the computational tool DAVID. Thirty-eight of the 72 deregulated miRNAs possessed mRNA targets that have been experimentally observed; in total 609 genes are potentially regulated (Additional file 1). These genes were mapped to KEGG pathways and were shown to be broadly involved in cancer development (Additional file 2).
Specifically, members of the miR-17-92 cluster were deregulated in our dataset: miR-19a and miR-19b were strongly up-regulated, in addition to moderate up-regulation of miR-17-3p/miR-17-5p and miR-92b. These results are in line with the observation that the miR-17-92 cluster is up-regulated in many cancer types, including lung cancer and lymphoma [27, 28]. Accordingly, miR-17-92 cluster members have been shown to take part in feedback loops determining the role of c-MYC as tumor suppressor and/or oncogene [29, 30]. Specifically, c-MYC apparently possesses a tumorigenic role in HNSCC, constituting a current candidate for anticancer strategies . Recently, the miR-17-92 cluster has been also shown to regulate multiple components of the TGF-β pathway in neuroblastoma . Other cancer-related miRNAs up-regulated in our OSCC samples are members of the miR-34 family: miR-34b and miR-34c. To our knowledge this is the first report of their altered expression profile in HNSCC, although the deregulation of miR-34a has been recently addressed in HNSCC . These results are interesting in light of the finding that miR-34 is a direct target of p53, functioning downstream of the p53 pathway as a tumor suppressor [34, 35]. Similar to other types of cancer, inactivation of p53 is an extremely common event in head and neck cancers, with mutant p53 status found in nearly 50% of the cases and commonly associated with poor prognosis . However, the role of miR-34b/c in the context of p53 regulation has not been addressed in HNSCC.
In agreement with most miRNA profiles in HNSCC samples and tumor cell lines, miR-133a was also down-regulated in our cancer set as compared to tumor-free samples. Its tumor suppressor activity, for instance by controlling the target genes actin-related protein 2/3 complex subunit 5 (ARPC5) and moesin (MSN), has been already demonstrated in squamous cell carcinoma of the tongue [37–39]. Since this seems to be a robust characteristic in HNSCC, its function should be further investigated as well as its possible use as a biomarker for early cancer detection.
Deregulation of homeobox cluster-encoded miRNAs miR-196a/b and miR-10b
MiR-196a/b was over-expressed and miR-10b was down-regulated in the OSCC samples compared with tumor-free surgical margins (Table 2). Both miRNAs are dysregulated in a variety of cancers [40, 41], but have not been previously associated with OSCC. We validated our microarray results in an additional subset of OSCC samples as well as in samples belonging to other HNSCC subsites (Figure 1, and Table 1 for sample characteristics). Both miRNAs clearly presented differential expression between tumor and tumor-free samples, suggesting a role in HNSCC.
MiR-196a/b and miR-10b are embedded within homeobox (HOX) clusters of developmental regulators . Schimanski and collaborators  demonstrated that HOX genes are targeted by miR-196, and HOX transcripts were also experimentally validated as miR-10 targets [44, 45]. Given that molecular events involved in carcinogenesis interfere in the regulation of cell identity, it has been proposed that HOX proteins could be oncogenic regulators . HOX genes have not been implicated in the development of HNSCC, as judged from reviewing the available literature, including HNSCC gene expression profiles. This suggests that the homeobox-cluster embedded miRNA could have a different role in HNSCC. Thus, we performed gain-of-function experiments aiming to outline a possible role for these molecules.
Gain-of-function of miR-10b and miR-196a lead to impaired cell proliferation
Precursor molecules of miR-10b were transfected into SCC25 and SCC9 (tongue squamous cell carcinoma-derived cell lines) and FaDu (a cell line derived from hypopharyngeal squamous cell carcinoma), while miR-196a precursor molecules were transfected into human keratinocytes derived from normal oral epithelium. We chose SCC cell lines and oral keratinocytes as models for the investigation of miRNA function in a HNSCC genetic background, emulating cancer and tumor-free cellular systems, respectively.
Two SCC cell lines were initially chosen for the gain-of-function experiments due to differences in their proliferation rates, as reported in Figure 2 and Table 3, a characteristic that could affect results.
As expected, in untreated SCC cell lines, miR-196 was up-regulated and miR-10b was downregulated when expression levels were compared to untreated keratinocytes (data not shown).
After transfection, we assessed the over-expresssion of the respective mature miRNAs in each cell line (Figure 3). Despite the successful overexpression of miR-10b in SCC9, these cells were very sensitive to the transfection, with intense effects in cell proliferation seen also for the negative control. Thus this cell line was not used in subsequent experiments.
Since tumor cells evade programmed cell death and sustain proliferative status , we tested whether miR-10b and/or miR-196a could play a role in this scenario. Assessment of Ki-67 antigen expression, a cell proliferation marker, revealed that keratinocytes over-expressing miR-196a were mostly quiescent, as defined by a lack of Ki-67 antigen expression. When compared to transfection controls, cell proliferation was reduced approximately 5 fold in keratinocytes over-expressing miR-196a (Figure 4).
Despite inhibiting cell proliferation, specific occurrences at the cell surface level, such as surface blebbing, considered as a marker for apoptosis, were absent at 72 h after transfection (data not shown). Changes in nuclear morphology indicating late apoptosis were also absent.
Inhibition of cell proliferation upon over-expression of miR-10b in SCC25 and FaDu was assessed by flow cytometry (Figure 5). This is consistent with the miRNA expression data: miR-10b was detected in low levels in HNSCC in this study, a downregulation that would thwart the inhibitory effect of miR-10b on cell proliferation. The effect of the overexpression of miR-10b was clearly more intense in FaDu as compared to SCC25.
Given the fact that miR-196a was found to be overexpressed in HNSCC samples, the inhibitory effect of miR-196a overexpression on proliferation of normal keratinocyes cannot easily put into context. One hypothesis could be that miR-196a overexpression in HNSCC could also be the consequence of uncontrolled proliferation, as a means to counteract it, rather than the cause of its perturbation.
Efforts to understand the global effects of miR-196a, which might be cell-type dependent, are essential considering that it has been recently addressed as a potential therapeutic target .
Gene expression profiles upon overexpression of miR-10b and miR-196a do not show regulation of known targets
The identification of miRNA target genes is critical in order to understand their roles. However, this task is challenging. MiRNAs are usually imperfectly complementary to the 3′UTR region of their mRNA targets in mammalian cells, with target effect hardly detected at the gene expression level. Additionally, the cellular environment is key in determining miRNA functions, which will vary depending on the cell type.
While keeping these drawbacks in mind, global gene expression profiling using DNA microarrays was performed in this study to identify deregulated cellular processes upon the overexpression of miR-10b and miR-196a in each cell model. Broad consequences of the overexpression should be detected by global gene expression profiling, even when the deregulation of specific target genes might not be detected by this kind of experiment.
Experimentally validated mRNA targets were searched in Tarbase and miRecord databases. None of the miR-10b targets HOXA1, HOXD10 and KLF4[40, 48, 49] were affected at the mRNA level by the overexpression of miR-10b in SCC25 or FaDu (Additional files 3 and 4, respectively). The same was true for the miR-196a gene targets ANXA1, HOXA7, HOXB8, HOXC8, HOXD8, KRT5 and S100A9[42, 50, 51] (Additional file 5). These results suggest that, at least at the mRNA levels these genes are not targeted by miR-196a in the cells used here.
Among the above mentioned gene targets, only ANXA1 down-regulation has been previously reported in HNSCC [52, 53]. For this reason, we checked for alterations of this target also at the protein level. Our results demostrate that ANXA1 is not targeted by miR-196a under the conditions studied here (Figure 6).
MiR10b and miR196a lead to cell cycle arrest through different mechanisms
We performed a functional analysis of deregulated genes aiming to pinpoint alterations that could explain impaired proliferation. A total of 353 annotated genes were downregulated (at 1.5 fold-change) following miR-196a over-expression in keratinocytes (Additional file 5). The relationships among these genes were assessed using Ingenuity Pathway Analysis (IPA), while considering only experimentally proven connections between genes or proteins. The most significant interaction network consisted of genes associated with DNA replication, recombination and repair, cell cycle and, consequently, cancer. Figure 7 depicts this network and genes involved in cell cycle arrest are highlighted. This network contains 8 deregulated genes from our dataset: CDK2, SYNM (DMN), TP73, AKT1, NFATC4, HOXA9, HSPB3 and CD40LG. Of particular interest is the downregulation of CDK2 and AKT1 and the upregulation of TP73. CDK2 is a subunit of the cyclin-dependent protein kinase complex, expressed in G1-S phase, and essential for cell cycle G1/S phase transition. TP73, up-regulated in cells overexpressing miR-196a, transcriptionally activates target genes leading to apoptosis and growth arrest . The activation of PI3K/AKT pathway in HNSCC is well known; the pathway regulates cell proliferation and has been addressed as a therapeutical target. Thus, the expression patterns of these three genes, following over-expression of miR-196a, would be in agreement with the observed arrest of the cell cycle. However, none of them are direct targets of this miRNA and further studies are needed in order to comprehend the observed effect.
Overexpression of miR-10b in SCC25 and in FaDu provided relatively similar results. Two hundred and ten annotated genes were downregulated and 169 were upregulated when SCC25 cells overexpressing miR-10b were compared to controls (Additional file 3) while 161 genes were downregulated and 169 upregulated in FaDu overexpressing miR-10b (Additional file 4), when at least a 2-fold difference was considered. Sixteen common genes were downregulated in both cell lines, but none of these genes were miR-10b predicted targets.
Regulatory networks provided by IPA did not contain a significant number of genes directly implicated in cell proliferation or cell cycle arrest for SCC25 cell line. This analysis, however, highlighted enrichment of terms belonging to the G-protein-coupled-receptor signaling pathway - with 9 molecules regulated in our dataset (ADRB3, AVPR2 and GRM2, upregulated 2-fold in SCC25 overexpressing miR-10b; and CNR1 (16-fold), DRD3 (4-fold), HCAR2 (2-fold) and OPRK1 (14-fold), downregulated in these cells). A recent review addresses mechanisms by which G-protein-coupled-receptors participate in the regulation of cell cycle  and, in the context of HNSCC, G-protein-coupled-receptors have been associated with EGFR signaling and cell survival .
A significant regulatory network built with deregulated genes upon overexpression of miR-10b in FaDu includes genes involved in the regulation of cell cycle progression and arrest (Figure 8). Although none of these genes have been implicated in HNSCC or heavily studied in the context of cancer, it is noteworthy the fact that they relate to cell cycle regulation through key players in HNSCC: TP53, NOTCH1, MYC and HRAS .
From this analysis it became clear that the effect of the overexpression of miR-10b in SCC25 and FaDu, and miR-196a in keratinocytes do not act upon a large number of cellular processes but may rather target a small set of genes, some of which directly or indirectly involved in the progression of cell cycle.
Data on miRNA effects in tumorigenesis and cancer progression is still controversial and should vary with cell and cancer types. While individual miRNAs might possess numerous and distinct targets, they should be able to contribute to the same tumorigenic processes through complex, and still mostly unknown, networks. In HNSCC little is known about the contribution of miRNA to tumor development and progression, with several studies lacking corroboration. Besides presenting data matching to current knowledge, in this study we show that two miRNAs, miR-196a and miR-10b, play distinct roles in the regulation of cell proliferation within a HNSCC background.
Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5.
Jopling CL, Yi M, Lancaster AM, Lemon SM, Sarnow P: Modulation of hepatitis C virus RNA abundance by a liver-specific MicroRNA. Science. 2005, 309 (5740): 1577-1581. 10.1126/science.1113329.
Vasudevan S, Tong YC, Steitz JA: Switching from repression to activation: MicroRNAs can up-regulate translation. Science. 2007, 318 (5858): 1931-1934. 10.1126/science.1149460.
Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, Sweet-Cordero A, Ebet BL, Mak RH, Ferrando AA, et al: MicroRNA expression profiles classify human cancers. Nature. 2005, 435 (7043): 834-838. 10.1038/nature03702.
Michael MZ, O’Connor SM, Pellekaan NGV, Young GP, James RJ: Reduced accumulation of specific microRNAs in colorectal neoplasia. Mol Cancer Res. 2003, 1 (12): 882-891.
Bartels CL, Tsongalis GJ: MicroRNAs: novel biomarkers for human cancer. Clin Chem. 2009, 55 (4): 623-631. 10.1373/clinchem.2008.112805.
Forastiere A, Koch W, Trotti A, Sidransky D: Head and neck cancer. N Engl J Med. 2001, 345 (26): 1890-1900. 10.1056/NEJMra001375.
Mao L, Hong WK, Papadimitrakopoulou VA: Focus on head and neck cancer. Cancer Cell. 2004, 5 (4): 311-316. 10.1016/S1535-6108(04)00090-X.
Liu X, Chen Z, Yu J, Xia J, Zhou X: MicroRNA profiling and head and neck cancer. Comp Funct Genet. 2009, 837514-
Babu JM, Prathibha R, Jijith VS, Hariharan R, Pillai MR: A miR-centric view of head and neck cancers. Biochim Biophys Acta. 2011, 1816 (1): 67-72.
Wua BH, Xiong XP, Jia J, Zhang WF: MicroRNAs: New actors in the oral cancer scene. Oral Oncol. 2011, 47 (5): 314-319. 10.1016/j.oraloncology.2011.03.019.
Childs G, Fazzari M, Kung G, Kawachi N, Brandwein-Gensler M, McLemore M, Chen Q, Burk RD, Smith RV, Prystowsky MB, et al: Low-level expression of microRNAs let-7d and miR-205 are prognostic markers of head and neck squamous cell carcinoma. Am J Pathol. 2009, 174 (3): 736-745. 10.2353/ajpath.2009.080731.
Avissar M, McClean MD, Kelsey KT, Marsit CJ: MicroRNA expression in head and neck cancer associates with alcohol consumption and survival. Carcinogenesis. 2009, 30 (12): 2059-2063. 10.1093/carcin/bgp277.
Christensen BC, Moyer BJ, Avissar M, Ouellet LG, Plaza SL, McClean MD, Marsit CJ, Kelsey KT: A let-7 microRNA-binding site polymorphism in the KRAS 3’ UTR is associated with reduced survival in oral cancers. Carcinogenesis. 2009, 30 (6): 1003-1007. 10.1093/carcin/bgp099.
Avissar M, Christensen BC, Kelsey KT, Marsit CJ: MicroRNA expression ratio is predictive of head and neck squamous cell carcinoma. Clin Cancer Res. 2009, 15 (8): 2850-2855. 10.1158/1078-0432.CCR-08-3131.
Fletcher AM, Heaford AC, Trask DK: Detection of metastatic head and neck squamous cell carcinoma using the relative expression of tissue-specific mir-205. Transl Oncol. 2008, 1 (4): 202-208.
Tanzer A, Amemiya CT, Kim CB, Stadler PF: Evolution of microRNAs located within Hox gene clusters. J Exp Zool B Mol Dev Evol. 2005, 304B (1): 75-85. 10.1002/jez.b.21021.
Yekta S, Tabin CJ, Bartel DP: MicroRNAs in the Hox network: an apparent link to posterior prevalence. Nat Rev Genet. 2008, 9 (10): 789-796. 10.1038/nrg2400.
Breitling R, Armengaud P, Amtmann A, Herzyk P: Rank products: a simple, yet powerful, new method to detect differentially regulated genes in replicated microarray experiments. FEBS Lett. 2004, 573 (1–3): 83-92.
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(T)(−Delta Delta C) method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3 (7): @@-Epub 2002 Jun 18
Klingbeil MFG, Herson MR, Cristo EB, Pinto DD, Yoshito D, Mathor MB: Comparison of two cellular harvesting methods for primary human oral culture of keratinocytes. Cell Tissue Bank. 2009, 10 (3): 197-204. 10.1007/s10561-009-9122-7.
Ziober AF, Patel KR, Alawi F, Gimotty P, Weber RS, Feldman MM, Chalian AA, Weinstein GS, Hunt J, Ziober BL: Identification of a gene signature for rapid screening of oral squamous cell carcinoma. Clin Cancer Res. 2006, 12 (20 Pt 1): 5960-5971.
Severino P, Alvares AM, Michaluart P, Okamoto OK, Nunes FD, Moreira-Filho CA, Tajara EH: Global gene expression profiling of oral cavity cancers suggests molecular heterogeneity within anatomic subsites. BMC Res Notes. 2008, 1: 113-10.1186/1756-0500-1-113.
Wu BH, Xiong XP, Jia J, Zhang WF: MicroRNAs: new actors in the oral cancer scene. Oral Oncol. 2011, 47 (5): 314-319. 10.1016/j.oraloncology.2011.03.019.
Tran N, O’Brien CJ, Clark J, Rose B: Potential role of micro-RNAs in head and neck tumorigenesis. Head Neck. 2010, 32 (8): 1099-1111. 10.1002/hed.21356.
Hayashita Y, Osada H, Tatematsu Y, Yamada H, Yanagisawa K, Tomida S, Yatabe Y, Kawahara K, Sekido Y, Takahashi T: A polycistronic microRNA cluster, miR-17-92, is overexpressed in human lung cancers and enhances cell proliferation. Cancer Res. 2005, 65 (21): 9628-9632. 10.1158/0008-5472.CAN-05-2352.
Olive V, Jiang I, He L: mir-17-92, a cluster of miRNAs in the midst of the cancer network. Int J Biochem Cell Biol. 2010, 42 (8): 1348-1354. 10.1016/j.biocel.2010.03.004.
Woods K, Thomson JM, Hammond SM: Direct regulation of an oncogenic micro-RNA cluster by E2F transcription factors. J Biol Chem. 2007, 282 (4): 2130-2134.
Aguda BD, Kim Y, Piper-Hunter MG, Friedman A, Marsh CB: MicroRNA regulation of a cancer network: consequences of the feedback loops involving miR-17-92, E2F, and Myc. Proc Natl Acad Sci U S A. 2008, 105 (50): 19678-19683. 10.1073/pnas.0811166106.
Xu B, Liu P, Li J, Lu H: c-MYC depletion potentiates cisplatin-induced apoptosis in head and neck squamous cell carcinoma: involvement of TSP-1 up-regulation. Ann Oncol. 2010, 21 (3): 670-672.
Mestdagh P, Bostrom AK, Impens F, Fredlund E, Van Peer G, De Antonellis P, von Stedingk K, Ghesquiere B, Schulte S, Dews M, et al: The miR-17-92 microRNA cluster regulates multiple components of the TGF-beta pathway in neuroblastoma. Mol Cell. 2010, 40 (5): 762-773. 10.1016/j.molcel.2010.11.038.
Kumar B, Yadav A, Lang J, Teknos TN, Kumar P: Dysregulation of microRNA-34a expression in head and neck squamous cell carcinoma promotes tumor growth and tumor angiogenesis. PLoS One. 2012, 7 (5): e37601-10.1371/journal.pone.0037601.
Corney DC, Flesken-Nikitin A, Godwin AK, Wang W, Nikitin AY: MicroRNA-34b and MicroRNA-34c are targets of p53 and cooperate in control of cell proliferation and adhesion-independent growth. Cancer Res. 2007, 67 (18): 8433-8438. 10.1158/0008-5472.CAN-07-1585.
Hermeking H: The miR-34 family in cancer and apoptosis. Cell Death Differ. 2010, 17 (2): 193-199. 10.1038/cdd.2009.56.
Koch WM, Brennan JA, Zahurak M, Goodman SN, Westra WH, Schwab D, Yoo GH, Lee DJ, Forastiere AA, Sidransky D: p53 mutation and locoregional treatment failure in head and neck squamous cell carcinoma. J Natl Cancer Inst. 1996, 88 (21): 1580-1586. 10.1093/jnci/88.21.1580.
Wong TS, Liu XB, Ho ACW, Yuen APW, Ng RWM, Wei WI: Identification of pyruvate kinase type M2 as potential oncoprotein in squamous cell carcinoma of tongue through microRNA profiling. Int J Cancer. 2008, 123 (2): 251-257. 10.1002/ijc.23583.
Kinoshita T, Nohata N, Watanabe-Takano H, Yoshino H, Hidaka H, Fujimura L, Fuse M, Yamasaki T, Enokida H, Nakagawa M, et al: Actin-related protein 2/3 complex subunit 5 (ARPC5) contributes to cell migration and invasion and is directly regulated by tumor-suppressive microRNA-133a in head and neck squamous cell carcinoma. Int J Oncol. 2012, 40 (6): 1770-1778.
Kinoshita T, Nohata N, Fuse M, Hanazawa T, Kikkawa N, Fujimura L, Watanabe-Takano H, Yamada Y, Yoshino H, Enokida H, et al: Tumor suppressive microRNA-133a regulates novel targets: Moesin contributes to cancer cell proliferation and invasion in head and neck squamous cell carcinoma. Biochem Biophys Res Commun. 2012, 418 (2): 378-383. 10.1016/j.bbrc.2012.01.030.
Lund AH: miR-10 in development and cancer. Cell Death Differ. 2010, 17 (2): 209-214. 10.1038/cdd.2009.58.
Chen C, Zhang Y, Zhang L, Weakley SM, Yao Q: MicroRNA-196: critical roles and clinical applications in development and cancer. J Cell Mol Med. 2011, 15 (1): 14-23. 10.1111/j.1582-4934.2010.01219.x.
Yekta S, Shih IH, Bartel DP: MicroRNA-directed cleavage of HOXB8 mRNA. Science. 2004, 304 (5670): 594-596. 10.1126/science.1097434.
Schimanski CC, Frerichs K, Rahman F, Berger M, Lang H, Galle PR, Moehler M, Gockel I: High miR-196a levels promote the oncogenic phenotype of colorectal cancer cells. World J Gastroenterol. 2009, 15 (17): 2089-2096. 10.3748/wjg.15.2089.
Baffa R, Fassan M, Volinia S, O’Hara B, Liu CG, Palazzo JP, Gardiman M, Rugge M, Gomella LG, Croce CM, et al: MicroRNA expression profiling of human metastatic cancers identifies cancer gene targets. J Pathol. 2009, 219 (2): 214-221. 10.1002/path.2586.
Ma L, Teruya-Feldstein J, Weinberg RA: Tumour invasion and metastasis initiated by microRNA 10b in breast cancer. Nature. 2007, 449 (7163): 682-U682. 10.1038/nature06174.
Wynter CVA: The dialectics of cancer: a theory of the initiation and development of cancer through errors in RNAi. Med Hypotheses. 2006, 66 (3): 612-635. 10.1016/j.mehy.2005.10.024.
Oganesian L, Karlseder J: Mammalian 5’ C-rich telomeric overhangs are a mark of recombination-dependent telomere maintenance. Mol Cell. 2011, 42 (2): 224-236. 10.1016/j.molcel.2011.03.015.
Garzon R, Pichiorri F, Palumbo T, Iuliano R, Cimmino A, Aqeilan R, Volinia S, Bhatt D, Alder H, Marcucci G, et al: MicroRNA fingerprints during human megakaryocytopoiesis. Proc Natl Acad Sci U S A. 2006, 103 (13): 5078-5083. 10.1073/pnas.0600587103.
Tian YY, Luo AP, Cai YR, Su Q, Ding F, Chen HY, Liu ZH: MicroRNA-10b promotes migration and invasion through KLF4 in human esophageal cancer cell lines. J Biol Chem. 2010, 285 (11): 7986-7994. 10.1074/jbc.M109.062877.
Luthra R, Singh RR, Luthra MG, Li YX, Hannah C, Romans AM, Barkoh BA, Chen SS, Ensor J, Maru DM, et al: MicroRNA-196a targets annexin A1: a microRNA-mediated mechanism of annexin A1 downregulation in cancers. Oncogene. 2008, 27 (52): 6667-6678. 10.1038/onc.2008.256.
Maru DM, Singh RR, Hannah C, Albarracin CT, Li YX, Abraham R, Romans AM, Yao H, Luthra MG, Anandasabapathy S, et al: MicroRNA-196a is a potential marker of progression during Barrett’s metaplasia-dysplasia-invasive adenocarcinoma sequence in esophagus. Am J Pathol. 2009, 174 (5): 1940-1948. 10.2353/ajpath.2009.080718.
Pedrero JMG, Fernandez MP, Morgan RO, Zapatero AH, Gonzalez MV, Nieto CS, Rodrigo JP: Annexin A1 down-regulation in head and neck cancer is associated with epithelial differentiation status. Am J Pathol. 2004, 164 (1): 73-79. 10.1016/S0002-9440(10)63098-2.
Alves VA, Nonogaki S, Cury PM, Wunsch-Filho V, de Carvalho MB, Michaluart-Junior P, Moyses RA, Curioni OA, Figueiredo DL, Scapulatempo-Neto C, et al: Annexin A1 subcellular expression in laryngeal squamous cell carcinoma. Histopathology. 2008, 53 (6): 715-727. 10.1111/j.1365-2559.2008.03186.x.
Melino G, De Laurenzi V, Vousden KH: p73: Friend or foe in tumorigenesis. Nat Rev Cancer. 2002, 2 (8): 605-615. 10.1038/nrc861.
New DC, Wong YH: Molecular mechanisms mediating the G protein-coupled receptor regulation of cell cycle progression. J Mol Signal. 2007, 2: 2-10.1186/1750-2187-2-2.
Cassell A, Grandis JR: Investigational EGFR-targeted therapy in head and neck squamous cell carcinoma. Expert Opin Investig Drugs. 2010, 19 (6): 709-722. 10.1517/13543781003769844.
Agrawal N, Frederick MJ, Pickering CR, Bettegowda C, Chang K, Li RJ, Fakhry C, Xie TX, Zhang J, Wang J, et al: Exome sequencing of head and neck squamous cell carcinoma reveals inactivating mutations in NOTCH1. Science. 2011, 333 (6046): 1154-1157. 10.1126/science.1206923.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/13/533/prepub
Financial support was provided by Fundaçao de Amparo a Pesquisa do Estado de Sao Paulo - FAPESP (05/51467-0 and 09/04166-5) and Albert Einstein Research and Education Institute – IIEP (IEP.PE.08-0125). CC and JR were funded by The Wellcome Trust. The authors acknowledge Dr. Flavio Borges for statistical analysis of miRNA expression and clinical data. The authors acknowledge the contribution of the GENCAPO (Brazilian Head and Neck Genome Project) for sample collection, clinical and pathological data collection and interpretation, and fruitful discussions (complete list of members and affiliations presented at http://www.gencapo.famerp.br).
The authors declare there are no competing interests.
Conceived and designed the experiments: PS JR. Performed the experiments: PS CC FMA WOP RMS MFGK. Interpreted experimental data: PS CC HB. Clinical data analysis and sample selection: RM VWF EHT FDN. Wrote the paper: PS HB. All authors read and approved the final manuscript.
Electronic supplementary material
Authors’ original submitted files for images
About this article
- Oral Squamous Cell Carcinoma
- Oral Squamous Cell Carcinoma
- Oral Squamous Cell Carcinoma Cell
- Oral Keratinocytes
- Oral Squamous Cell Carcinoma Cell Line