High-throughput sequencing of small RNA transcriptomes reveals critical biological features targeted by microRNAs in cell models used for squamous cell cancer research

  • Patricia Severino5Email author,

    Affiliated with

    • Liliane Santana Oliveira5, 6,

      Affiliated with

      • Natalia Torres5,

        Affiliated with

        • Flavia Maziero Andreghetto5,

          Affiliated with

          • Maria de Fatima Guarizo Klingbeil7,

            Affiliated with

            • Raquel Moyses4,

              Affiliated with

              • Victor Wünsch-Filho5,

                Affiliated with

                • Fabio Daumas Nunes6,

                  Affiliated with

                  • Monica Beatriz Mathor7,

                    Affiliated with

                    • Alexandre Rossi Paschoal6, 7 and

                      Affiliated with

                      • Alan Mitchell Durham6

                        Affiliated with

                        BMC Genomics201314:735

                        DOI: 10.1186/1471-2164-14-735

                        Received: 29 October 2012

                        Accepted: 17 October 2013

                        Published: 26 October 2013

                        Abstract

                        Background

                        The implication of post-transcriptional regulation by microRNAs in molecular mechanisms underlying cancer disease is well documented. However, their interference at the cellular level is not fully explored. Functional in vitro studies are fundamental for the comprehension of their role; nevertheless results are highly dependable on the adopted cellular model. Next generation small RNA transcriptomic sequencing data of a tumor cell line and keratinocytes derived from primary culture was generated in order to characterize the microRNA content of these systems, thus helping in their understanding. Both constitute cell models for functional studies of microRNAs in head and neck squamous cell carcinoma (HNSCC), a smoking-related cancer. Known microRNAs were quantified and analyzed in the context of gene regulation. New microRNAs were investigated using similarity and structural search, ab initio classification, and prediction of the location of mature microRNAs within would-be precursor sequences. Results were compared with small RNA transcriptomic sequences from HNSCC samples in order to access the applicability of these cell models for cancer phenotype comprehension and for novel molecule discovery.

                        Results

                        Ten miRNAs represented over 70% of the mature molecules present in each of the cell types. The most expressed molecules were miR-21, miR-24 and miR-205, Accordingly; miR-21 and miR-205 have been previously shown to play a role in epithelial cell biology. Although miR-21 has been implicated in cancer development, and evaluated as a biomarker in HNSCC progression, no significant expression differences were seen between cell types. We demonstrate that differentially expressed mature miRNAs target cell differentiation and apoptosis related biological processes, indicating that they might represent, with acceptable accuracy, the genetic context from which they derive. Most miRNAs identified in the cancer cell line and in keratinocytes were present in tumor samples and cancer-free samples, respectively, with miR-21, miR-24 and miR-205 still among the most prevalent molecules at all instances. Thirteen miRNA-like structures, containing reads identified by the deep sequencing, were predicted from putative miRNA precursor sequences. Strong evidences suggest that one of them could be a new miRNA. This molecule was mostly expressed in the tumor cell line and HNSCC samples indicating a possible biological function in cancer.

                        Conclusions

                        Critical biological features of cells must be fully understood before they can be chosen as models for functional studies. Expression levels of miRNAs relate to cell type and tissue context. This study provides insights on miRNA content of two cell models used for cancer research. Pathways commonly deregulated in HNSCC might be targeted by most expressed and also by differentially expressed miRNAs. Results indicate that the use of cell models for cancer research demands careful assessment of underlying molecular characteristics for proper data interpretation. Additionally, one new miRNA-like molecule with a potential role in cancer was identified in the cell lines and clinical samples.

                        Background

                        MicroRNAs (miRNAs) are small RNA molecules, typically between 19 and 22 nucleotides in length, that regulate protein-coding genes through sequence-specific binding to messenger RNAs (mRNAs). These molecules were first implicated in Caenorhabditis elegans development in the early 90s [1], and have been associated with a variety of biological processes since then [2]. They are believed to have important roles in cancer aethiology and progression, and are currently being evaluated for cancer classification and prognosis [3].

                        In order to identify cell processes that are affected by miRNAs, overexpression and inhibition of miRNA genes are routinely performed for in vitro functional studies. Well-established cell lines are generally used for this purpose, since they are readily available from certified sources that guarantee the genetic identity. Primary cultures, on the other hand, have the advantage of not presenting genetic changes associated with the process of obtaining immortalized cell lines. It is clear, however, that experimental data will be biased by the adopted cell model and should be interpreted with caution.

                        Head and neck squamous cell carcinoma (HNSCC) is a smoking-related cancer for which, despite being one of the most common malignancies worldwide, reliable diagnostic and prognostic markers are not available [4]. Recent studies have addressed deregulation of microRNAs in the context of HNSCC, suggesting that these molecules could be used to improve diagnosis and the outcome of this disease [5].

                        Cell models are being broadly used in order to comprehend the function of specific miRNAs in this disease, but there is no current information on the miRNA content of these cells. Due to the complexity of miRNA regulatory networks, where one miRNA may target multiple genes, and where a single gene might be targeted by several miRNAs, miRNA expression levels (the miRNome) of a cell line may lead to a significant bias in functional studies results. Aiming to understand the miRNA background of cell models used for functional studies in HNSCC, we report the high-throughput sequencing analysis of the small RNA transcriptome of an oral squamous cell carcinoma cell line (SCC25) and of normal oral keratinocytes obtained from primary cultures. We believe that a deep understanding of the molecular background of cell models should greatly improve knowledge on mechanisms targeted by miRNAs and other gene regulators.

                        Results

                        MiRNAs expressed in the carcinoma cell line and in normal keratinocytes

                        We aimed to identify miRNAs that are expressed in a human carcinoma cell line and in a cell type representing its normal counterpart for miRNA functional studies in the context of HNSCC. With this purpose we chose SCC25, a tongue squamous cell carcinoma cell line, and normal oral keratinocytes derived from primary cultures. Three small RNA libraries were constructed for each cell type and sequenced on a SOLiD sequencer. Table  1 shows that, after filtering out rRNA, tRNA and other contaminants, 2.62% and 3.10% of reads accounted for miRNAs previously annotated in miRBase (release 18) for SCC25 and normal keratinocytes, respectively. When we considered only reads matching mature miRNAs, these numbers dropped to about 2% in both cases. Three hundred and thirty one mature miRNAs were identified in keratinocytes and 288 in SCC25, with 202 common miRNAs. The complete set of detected mature miRNAs for both cell types can be found in Additional file 1.
                        Table 1

                        Mapping of reads to genome and miRBase v18 after rRNA, tRNA, repeated DNA and adaptor sequences filtering

                         

                        SCC25 (million)

                        %

                        Kerat (million)

                        %

                        Total of reads

                        3.4

                        100

                        4.4

                        100

                        Reads matching miRBase V18

                        0.09

                        2.64

                        0.14

                        3.1

                        Reads matching mature miRNAs

                        0.07

                        2.00

                        0.9

                        1.94

                        Reads mapped to genome and not matching v18

                        0.13

                        3.97

                        0.07

                        4.53

                        Figure  1 depicts the 10 most represented miRNAs in SCC25 and keratinocytes. These miRNAs represent 72% of reads mapping mature miRNAs in keratinocytes and 80% in SCC25. The three most expressed miRNAs were common to both datasets: miR-21, miR-205 and miR-24. Both miR-21 and miR-205 have been reported as HNSCC biomarkers [6], but here there was no differential expression between the normal cell and the cancer cell line. Noteworthy is the fact that miR-21 and miR-205 are involved in keratinocyte proliferation and migration [2, 7], a common characteristic that might justify their ubiquous presence in these epithelium-derived cells. MiR-24, on the other hand, has not been implicated in HNSCC. A few reports mention a tumor suppressor function for miR-24 through the regulation of apoptosis [8]. We assessed the expression levels of miR-21 and miR-24 by real-time PCR, confirming that there was no significant difference in the expression of the two miRNAs between the two cell types (data not shown).
                        http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-14-735/MediaObjects/12864_2012_5450_Fig1_HTML.jpg
                        Figure 1

                        Histogram indicating the expression levels of 10 most represented miRNAs in both datasets. (A) 10 most expressed miRNAs in keratinocytes; (B) 10 most expressed miRNAs in the cell line (SCC25).

                        Together miR-21, miR-24 and miR-205 may target 66 genes (Additional file 2A). The task of identifying miRNA gene targets is an essential step for the understanding of their role in a given cellular environment. However, it is also a challenging undertaking due to the fact that miRNAs are usually imperfectly complementary to the 3′UTR region of their mRNA targets. Considering these difficulties, in order to access cellular processes targeted by the three most and commonly expressed miRNAs in the cell line and the normal keratinocyte, we selected only gene targets experimentally verified, as reported by miRecords and TarBase databases. Most of the 66 validated targets are reportedly involved in cell cycle regulation and programmed cell death (Additional file 2B). We emphasize the involvement of these miRNAs in the regulation of cyclin-dependent kinases (CDKs), key components in cell cycle regulation and traditional targets for cancer therapy development, in the regulation of SMAD and of type II receptor of TGFB (TGFBR2), components of the transforming growth factor beta signaling pathway, as well as in the regulation of NOTCH1, PI3K and PTEN. A recent review addresses the relationship between Smad/TGFBR2 and NOTCH1, PTEN and PI3K, considered relatively new players in HNSCC [9].

                        Only seven miRNAs were detected uniquely in keratinocytes, when a read count of at least ten reads was considered: miR-4500, miR-154, miR-337, miR-493, miR-326, miR-369, miR-381 and miR-627. Only one miRNA was uniquely expressed in SCC25 with this cut-off: miR-615. To our knowledge, no targets have been validated for these molecules up to now.

                        Differential regulation of miRNAs between the cancer cell line and normal keratinocytes

                        The relative abundance of miRNAs was then compared between the carcinoma cell line and normal keratinocytes. Among the 202 miRNAs expressed in both cell types, 65 miRNAs were overexpressed in keratinocytes and 42 in SCC25, when a 2-fold change in expression was considered. This result is in agreement with previous findings reporting overall lower levels of miRNAs in cancer as compared to normal tissues. Twenty miRNAs presenting higher fold-changes are depicted in Figure  2, where blue bars correspond to miRNAs mostly expressed in keratinocytes and white bars correspond to miRNAs mostly expressed in the cell line. Additional file 3 summarizes differences in gene expression for 107 miRNAs expressed with a 2-fold change between keratinocytes and the cell line.
                        http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-14-735/MediaObjects/12864_2012_5450_Fig2_HTML.jpg
                        Figure 2

                        Twenty most differentially expressed miRNAs between keratinocytes and SCC25. The expression level of miRNAs mostly expressed in keratinocytes is represented in blue and in white is the expression of those more expressed in the cell line.

                        From the 107 miRNAs, we found experimentally validated targets for only 6 molecules – miR-1, miR-125b and miR-133a up-regulated in keratinocytes and miR-196a, miR-511 and miR-7 up-regulated in the cell line. Among these, we chose to confirm by real time PCR the deregulation of the miRNAs presenting the lowest level of deregulation – miR-1 and miR-133a, up-regulated ~ 2-fold in keratinocytes, and of miR-196a, up-regulated ~ 2 fold in the cell line. Results between real time PCR and sequencing were consistent (data not shown). These 6 miRNAs may target 315 gene targets (Additional file 4A). Programmed cell death-related genes were the most represented among targeted genes, mostly grouped under the GO term Regulation of Apoptosis (GO:0042981, FDR corrected p value < 0.05) (Additional file 4B). This result indicates differences in apoptosis regulation between the cancer cell line and the normal counterpart. Illustrating the importance of this scenario for HNSCC, p53, a key tumor suppressor frequently mutated in HNSCC [10] is a validated target of miR-125b [11], and this miRNA has also been recently implicated in the carcinogenesis of squamous cell carcinomas [12].

                        Cell differentiation-related genes, grouped under GO term Positive Regulation of Cell Differentiation (GO:0045598, FDR corrected p value < 0.05), were also enriched in the dataset (Additional file 4B), an expected result considering intrinsic differences between a cell line and a keratinocyte derived from primary cultures, the later retaining its full capacity for cell differentiation.

                        Top over-expressed miRNAs in keratinocytes included miR-409, reportedly a tumor suppressor in gastric cancer [13], miR-3545, a recently described miRNA corresponding to the antisense strand of miR-203, which is a miRNA preferentially expressed in the skin, and which has been described as a tumor suppressor molecule silenced in different malignancies [14, 15], and miR-376c, shown to induce apoptosis and described as down regulated in HNSCC [16].

                        MiR-495 and miR-379 have not been previously associated with HNSCC biology or epithelial differentiation, but miR-495 has been reported to be a tumor suppressor [17] and miR-379 has been shown to be down-regulated in HPV positive HNSCC samples [18].

                        Among top over-expressed miRNAs in SCC25, miR-210 and miR-181c have been previously associated with HNSCC. MiR-210 has been described as a marker for tumor hypoxia and a prognostic factor in HNSCC [19]. MiR-181c, reported as up-regulated in tongue squamous cell carcinoma (TSCC), was the second most up-regulated miRNA in SCC25 [20].

                        Taken together these results highlight important differences between a cancer cell line derived from the oral cavity, and its normal counterpart, a keratinocyte derived from a primary culture. These differences are consistent with the genetic background they represent and should be taken into consideration when cell models are chosen.

                        mRNA expression patterns in the cell line and keratinocytes

                        DNA microarrays were used in order to evaluate gene expression differences between keratinocytes and the cancer cell line that might support processes targeted by the differential miRNA expression. Even though miRNAs are mostly posttranscriptional regulators and their effect on mRNA abundance would be small in our cell models, we aimed to access differentially regulated biological processes and, when possible, associate these differences with miRNA regulation potential. A total of 978 genes were over-expressed in SCC25 when a 2-fold difference in gene expression was considered, while 523 were down-regulated (FDR adjusted p value of 0.05) (Additional file 5).

                        Gene Ontology term enrichment analysis of the differentially expressed genes dataset identified processes associated with cell cycle regulation (cell cycle progression, chromosome segregation and DNA replication) enriched when genes up-regulated in SCC25 were considered (Additional file 6A). Genes up-regulated in keratinocytes were mostly involved with epithelium development (Additional file 6B). These results are mostly in agreement with expected differences between the cell types, since up-regulation of proliferation-related processes should be a characteristic of a cancer cell line and epithelium development of a normal keratinocyte.

                        We then looked for gene targets for the 7 miRNAs that possess experimentally validated targets addressed in the previous section – miR-1, miR-133a and miR-125b up-regulated in keratinocytes and miR-196a, miR-511 and miR-7 up-regulated in the cell line – in the dataset of deregulated mRNAs. A total of 21 targets were found among deregulated genes (Table  2). However, only eight of these targets presented expression levels inversely correlated with the related miRNA, an indication of regulation. Although the correlation of expression levels does not confirm the regulation of targets by a miRNA, nor does the lack of gene expression differences indicate that the targets are not being regulated - since regulation might not be at the gene expression level and we did not check protein levels - we chose to present only targets previously shown to be regulated at the gene expression level.
                        Table 2

                        Deregulated miRNA targets in dataset comparing mRNA expression between the cell line and keratinocytes

                        Gene target

                        Up-regulated in cell line

                        Up-regulated in Keratinocyte

                        Related miRNA

                        ADAR

                        3.6

                         

                        miR-1

                        ANKRD29

                         

                        4.1

                        miR-1

                        ANP32B

                        3.4

                         

                        miR-1

                        ATP6V1B2

                         

                        2

                        miR-1

                        CDCP1

                         

                        2.3

                        miR-1

                        EPB41L4B

                         

                        4.3

                        miR-1

                        FSTL1

                         

                        4.7

                        miR-1

                        KCNQ1

                        8

                         

                        miR-1

                        POLA2

                        2.9

                         

                        miR-1

                        PREX1

                        12

                         

                        miR-1

                        SSNA1

                        2.6

                         

                        miR-1

                        UST

                        4.6

                         

                        miR-1

                        CDK9

                        2.3

                         

                        miR-1

                        B3GALT4

                        2.7

                         

                        miR-125

                        BMF

                        31

                         

                        miR-125

                        CDKN2A

                         

                        21

                        miR-125

                        IGFBP3

                        700

                         

                        miR-125

                        RHOA

                         

                        2.3

                        miR-133a

                        SET

                        2.5

                         

                        miR-196a

                        SPRR2C

                         

                        71

                        miR-196a

                        SNCA

                         

                        5

                        miR-7

                        The table reports 21 validated gene targets presenting deregulation between SCC25 and keratinocytes. Fold-change in bold-type indicates that the fold-change is in agreement with the expression level of the regulator miRNA in the cell.

                        Identification of miRNA-like molecules

                        For the identification of novel miRNA-like molecules, all reads shorter that 35 nt, having more than 10 copies in our libraries and no match to known miRNAs were tested. Reads were mapped to the genome and, at each genomic locus of a read to be tested, one longer sequence covering the read was extracted for secondary structure analysis. This sequence extended 100 nt upstream and 100 nt downstream from the read and this length was established following the average length of known human precursor miRNAs (Additional file 7). A set of 448 sequences obtained this way was submitted to a four-step annotation procedure: (i) sequence similarity search with the BLAST program against a locally curated ncRNA database; (ii) structural search against all RFAM families using INFERNAL, (iii) ab initio characterization using RNAfold and HHMMiR; (iv) manual curation of the results (Figure  3). Candidates that mapped to known coding regions were discarded. Of the 448 sequences, 13 presented some evidence in at least one of the procedures. The results are shown in Table  3.
                        http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-14-735/MediaObjects/12864_2012_5450_Fig3_HTML.jpg
                        Figure 3

                        Workflow for the analysis of small RNA transcriptome focusing on miRNA identification and discovery. Each library was converted to FASTA format, subjected to RNA2MAP tool and reads that did not match the human genome were discarded. After the filtering protocol, that matched miRBase were submitted to differential expression analysis and those that did not match were mapped to the genome. At each genomic locus the read was extended by100 nt up and downstream. Resulting sequences were subjected to the following pipeline: secondary structure prediction using RNAfold and Infernal against RFAM, blast against a local non-coding RNA database and, ab intio characterization using HHMMiR and RNAFold (HHMMiR in fact is performed using the data produced by RNAfold).

                        Table 3

                        Putative miRNAs identified using similarity search, structural search and ab initio prediction

                          

                        Cell line count(normalized)

                        Blast search against miRNAs in local DB1

                        HHMMR

                        RNAfold

                        Rfam/INFERNAL2

                        Candidate number

                        Genomic location

                        Kerat

                        SCC5

                         

                        score

                        Original read overlaps Helix?

                        Keal/mol

                        Family

                        ID

                        Strand

                        Bit score

                        Start

                        End

                        Original reads overlaps Helix?

                        1

                        Intergenic

                        3.9

                        9.6

                        -

                        0.67

                        No

                        -72,10

                        -

                        -

                        -

                        -

                        -

                        -

                        -

                        2

                        Intergenic

                        46

                        82.5

                        -

                        0.69

                        No

                        -64,49

                        -

                        -

                        -

                        -

                        -

                        -

                        -

                        3

                        Intron

                        4.4

                        -

                         

                        0.69

                        Partially

                        -51,30

                        -

                        -

                        -

                        -

                        -

                        -

                        -

                        4

                        Intergenic

                        3.9

                        51.6

                        -

                        0.66

                        Partially

                        -44,30

                        MIR245

                        RF00816

                        +

                        29.76

                        97

                        182

                        Total

                        5

                        Intron

                        -

                        8

                        -

                        0.7

                        No

                        53,52

                        -

                        -

                        -

                        -

                        -

                        -

                        -

                        6

                        Intergenic

                        -

                        4.8

                        -

                        0.66

                        No

                        -48,90

                        -

                        -

                        -

                        -

                        -

                        -

                        -

                        7

                        Intergenic

                        -

                        78.9

                        Original DB: smiRNAdb; Type: unclassified RNA QueryCov: 13%; SbjCov: 97% Id: 100%

                        0.66

                        No

                        -57,34

                        -

                        -

                        -

                        -

                        -

                        -

                        -

                        8

                        Intron

                        -

                        5.6

                        Origianl DB: smiRNAdb; Type: miRNA; QueryCov: 11%; SbjCov: 96% Id: 100%

                        -

                        -

                        57,10

                        -

                        -

                        -

                        -

                        -

                        -

                        -

                        9

                        Intergenic

                        -

                        11.2

                        Original DB: smiRNAdb; Type: unclassified RNA; QueryCov: 13%; SubjectCov: 94%; Id: 100%;

                        0.67

                        Total

                        -35,35

                        mir-450

                        RF00708

                        +

                        30.4

                        32

                        112

                        Parcial

                        mir-245

                        RF00816

                        +

                        29.76

                        35

                        120

                        Parcial

                        10

                        Intergenic

                        -

                        4.8

                        -

                        -

                        -

                        -122,82

                        mir-689

                        RF00871

                        +

                        66.15

                        17

                        91

                        No

                        11

                        Intron

                        -

                        4.8

                        -

                        0.66

                        Total

                        -83,69

                        -

                        -

                        -

                        -

                        -

                        -

                        -

                        12

                        Intron

                        -

                        4.0

                        Original DB: smiRNAdb; Type: unclassified RNA; QueryCov: 12%; SubjectCov: 100%: Id: 100%

                        0.68

                        Total

                        -50,30

                        -

                        -

                        -

                        -

                        -

                        -

                        -

                        13

                        Intron

                        -

                        15.2

                        -

                        0.67

                        No

                        -94,14

                        -

                        -

                        -

                        -

                        -

                        -

                        -

                        Notation: Original DB: Database from which sequence was downloaded into local database: Type: RNA type (misRNA if not specified); Query Cov: percentage of database sequence covered in the alignment; Subject Cov: percentage of candidate covered by the alignment; Id: alignment identity; 1: local miRNA database with sequences from databases designated by the Non-coding RNA Databases Resource (NRDR).

                        Four candidates presented good BLAST alignments with sequences in ncRNA databases (Cand7, Cand8, Cand9 and Cand12) and, of these one aligned with a ncRNA classified as miRNA (Cand8). Three candidates (Cand4, Cand9 and Cand10) presented good structural alignments to RFAM families, but of these only two alignments included the corresponding original reads (candidates Cand4 and Cand9) and only two candidates aligned to RFAM families that included mammalian sequences (Cand9 and Cand10). All candidates but for Cand8 and Cand10 were predicted by HHMMiR as including miRNAs, however only candidates Cand3, Cand4, Cand9, Cand11 and Cand12 had the original reads included in predicted stems (Additional file 8).

                        Five of the thirteen candidates showed evidence of the original read mapping (totally or partially) in a predicted stem in the precursor sequence considered (Additional file 8). These full sequences (reads plus extensions constituting putative precursors) were analyzed with the computational tool MatureBayes for the identification of possible mature sequences. As depicted in the Additional file 9, the algorithm did not find a potential functional part in the proposed precursor sequences (putative mature miRNAs) that contained the original read. This analysis underlines difficulties in novel miRNA discovery, showing that different approaches are needed for complementarity.

                        Comparison of results in cells with clinical samples

                        In order to compare miRNAs expressed in cell cultures with clinical samples, small RNA libraries from eight HNSCC patients were constructed and sequenced, following the same procedures described for the cells (results in Additional file 10). In each case we sequenced a tumor sample and surgical margins, the later representing cancer-free tissue. A PCA plot illustrates global characteristics of these samples: despite heterogeneity within the two groups, the expression level of shared molecules is able to differentiate cancer from cancer-free samples (Figure  4). Table  4 shows that similarities in miRNA population, when cancer samples were compared to their cell model counterpart, varied from 46% to about 70%. Since these percentages did not correlate with the sequencing coverage obtained in this study, we concluded that sequencing depth for each clinical sample was adequate.
                        http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-14-735/MediaObjects/12864_2012_5450_Fig4_HTML.jpg
                        Figure 4

                        Unsupervised classification by principal components analysis of cancer and cancer-free samples. Principal components analysis (PCA) was used to classify 15 samples (8 cancer samples and 7 cancer-free samples) based on the expression profile of 193 miRNAs expressed in all samples. Sample 306 M was not included in this analysis due to its high heterogeneity. The PCA plot depicts 70% of variability in the dataset.

                        Table 4

                        Common miRNAs identified in cells and clinical samples.

                        Sample

                        Total of miRNAs

                        Common miRNAs

                        % of common miRNAs with related cell type

                        196 M

                        313

                        219

                        69.97

                        196 T

                        279

                        178

                        63.80

                        240 M

                        321

                        223

                        69.47

                        240 T

                        298

                        184

                        61.74

                        277 M

                        558

                        294

                        52.69

                        277 T

                        335

                        187

                        55.82

                        296 M

                        642

                        309

                        48.13

                        296 T

                        393

                        212

                        53.94

                        306 M

                        589

                        286

                        48.56

                        306 T

                        641

                        269

                        41.97

                        321 M

                        495

                        273

                        55.15

                        321 T

                        515

                        243

                        47.18

                        333 M

                        350

                        240

                        68.57

                        333 T

                        503

                        246

                        48.91

                        349 M

                        585

                        292

                        49.91

                        349 T

                        480

                        221

                        46.04

                        Tumor samples were compared to the cancer cell line and tumor-free samples with normal keratinocytes.

                        Despite differences in absolute quantities, miRNAs highlighted as most expressed in cells are frequently present among the most expressed miRNAs in clinical samples (Figure  5). Noteworthy is again the widespread expression of miR-21, miR-24 and miR-205, already addressed in cells. This result corroborates the assumption that these abundant molecules are involved in processes common to epithelial cells and/or maintenance, and do not seem to be affected by the cancer phenotype.
                        http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-14-735/MediaObjects/12864_2012_5450_Fig5_HTML.jpg
                        Figure 5

                        Most expressed miRNAs in keratinocytes, cell line and in clinical samples. A: Most expressed miRNAs in keratinocytes (Krt) and in tumor-free samples; B: Most expressed miRNAs in the cancer cell line (SCC25) and in tumor samples. MiRNAs are reported in alphabetical and numerical order. The presence of a given miRNA in the dataset is indicated by gray color.

                        On the other hand, the scenario describing the differential expression between cancer samples vs cancer-free tissue is more complex and not directly comparable with results obtained for cell lines. This is expected due to the complexity of the addressed tumor type. For instance, miR-409, presenting the highest fold-change between the keratinocytes and the cancer cell line was over-expressed in cancer-free tissue in only four patients, while miR-210, highly over-expressed in the cell line when compared to keratinocytes, was not differentially expressed in our clinical samples. In Additional file 11 we show the comparisons between common miRNAs in cancer and cancer-free tissue and their cell counterpart.

                        Finally, we considered the expression levels in the genomic region containing the miRNA candidate proposed in this study. Table  5 depicts the expression of the two reads that constitute the candidate in each clinical sample. Samples 196, 306, 312, and 333 showed similar levels of expression found in cells. Results suggest that the molecule might be functional and have a role in HNSCC.
                        Table 5

                        Expression of Cand4 and Cand9 (reads from predicted miRNA) in clinical samples

                        Sample 196

                        Read count (Raw)

                        Normalized read count

                         

                        196 M

                        196 T

                        196 M

                        196 T

                        Cand99

                        159

                        129

                        41.19

                        63.66

                        Cand209

                        4

                        1

                        1.04

                        1.60

                        Sample 240

                        Read Count (Raw)

                        Normalized Read Count

                         

                        240 M

                        240 T

                        240 M

                        240 T

                        Cand99

                        1076

                        26

                        278.77

                        10.41

                        Cand209

                        40

                        -

                        10.36

                        -

                        Sample 277

                        Read Count (Raw)

                        Normalized Read Count

                         

                        277 M

                        277 T

                        277 M

                        277 T

                        Cand99

                        391

                        16

                        101.30

                        6.41

                        Cand209

                        8

                        15

                        2.07

                        6.01

                        Sample 296

                        Read Count (Raw)

                        Normalized Read Count

                         

                        296 M

                        296 T

                        296 M

                        296 T

                        Cand99

                        1566

                        335

                        405.72

                        134.13

                        Cand209

                        9

                        10

                        2.33

                        4.00

                        Sample 306

                        Read Count (Raw)

                        Normalized Read Count

                         

                        306 M

                        306 T

                        306 M

                        306 T

                        Cand99

                        222

                        166

                        57.52

                        66.46

                        Cand209

                        2

                        3

                        0.52

                        1.20

                        Sample 321

                        Read Count (Raw)

                        Normalized Read Count

                         

                        321 M

                        321 T

                        321 M

                        321 T

                        Cand99

                        86

                        170

                        22.28

                        68.06

                        Cand209

                        -

                        1

                        -

                        0.40

                        Sample 333

                        Read Count (Raw)

                        Normalized Read Count

                         

                        333 M

                        333 T

                        333 M

                        333 T

                        Cand99

                        164

                        985

                        42.49

                        394.38

                        Cand209

                        -

                        3

                        -

                        1.20

                        Sample 349

                        Read Count (Raw)

                        Normalized Read Count

                         

                        349 M

                        349 T

                        349 M

                        349 T

                        Cand99

                        177

                        50

                        45.86

                        20.02

                        Cand209

                        -

                        1

                        -

                        0.40

                        Samples 196, 306, 321 and 333 presented expression levels similar to findings in cells.

                        Discussion

                        Complex regulatory networks connect genes within cellular processes. MiRNAs have been recently added to this scenario, constituting an additional layer in gene regulation. The involvement of miRNAs in fundamental biological processes such as cell differentiation and programmed cell death, as well as their implication in innumerous human diseases, are currently known. The comprehension of their roles at a systems-level, however, is far from complete.

                        Once a miRNA is identified as a player in a pathological condition, functional studies using cell models are commonly used to address its involvement in gene regulation and consequent phenotypes. Due to the complex interactions between miRNAs and gene targets, the knowledge regarding the expression of miRNAs in available cell models should be carefully considered before tackling miRNA functional assays. Additionally, the response of cell lines to different stimuli varies due to, among other issues, their particular genetic background. For squamous cell carcinoma cell lines, distinct patterns of proliferation and survival upon treatment with drugs or environmental stress have been reported in the literature [2123]. Regulation of miRNAs certainly impacts these conclusions, even though it was not addressed in the above-mentioned studies.

                        In this study we addressed the miRNA constitution of a cancer cell line and normal oral keratinocytes, both used as models for head and neck squamous cell carcinoma and its cancer-free counterpart, respectively. Our results comparing the expression of miRNAs expressed in both cell types corroborate literature findings. For instance, the expression of miR-205 was comparable in both cell types, in agreement with its description as a marker of squamous epithelia [6], and higher expression levels of miR-125b in keratinocytes corroborates recently published data implicating the loss of miR-125 and HNSCC carcinogenesis [12]. On the other hand, unexpected results such as the finding of miR-21 equally expressed between the cell models, a miRNA commonly deregulated in cancer [24], was confirmed by differences in the expression level of this miRNA found between our clinical samples.

                        Through the relationship between miRNA expression levels and the targets they could be regulating we show that important cancer hallmarks i.e. cell death and cell proliferation regulation, could be intensely affected by the expression of these molecules. On the other hand, we also demonstrate that the expression of certain miRNAs, and, consequently, processes they target, might be independent of the cancer phenotype. For this result we addressed miRNAs that did not show differences in expression between the cell types and their experimentally validated gene targets reported in the literature. As expected, expression levels of miRNAs in clinical samples varied and did not always correlate with findings in the cell models.

                        We have also performed an in silico analysis to search for new miRNA candidates. It is currently known that high throughput sequencing technologies allow for the discovery of novel molecules, due to their inherent sensibility and accuracy. The analysis used structural alignment against known miRNA families, ab initio prediction using HHMMiR and RNAfold, and similarity search against a curated miRNA local database. Since we were sequencing small RNAs, we postulated that any miRNA sequenced in the process was a mature miRNA, and therefore had to be part of a stem in the precursor miRNA secondary structure. Of the 13 candidates with any miRNA evidence, in only five there was any evidence that the original read mapped totally or partially in a predicted stem (Cand03, Cand04, Cand09, Cand11, Cand12), these will be discussed in detail. Cand3 had only a HHMMiR prediction, which partially included the original sequence in a stem, and was expressed only in keratinocytes. Cand11 presented a good HHMMiR prediction, and was expressed only in SCC25 cells. Cand12, had evidence coming from a HHMMIR positive prediction and an alignment against an uncharacterized small ncRNA in smiRNAdb. Cand4 had a positive HHMMiR prediction and an alignment to a known miRNA family (RF00816). Cand9 also mapped to the same RFAM family (RF00816), with the original read located in a predicted stem. In fact, the genomic mapping of the candidates showed that they Cand4 and Cand9 correspond to adjacent genomic locations (Figure  6). The original reads of each candidate covered both sides of the same predicted stem in the miRNA family. Additionally both candidates were more expressed in SCC25 cells when compared to keratinocytes (51.6 vs 3.9 for Cand4, 11.2 vs. 0 for Cand9). A close investigation of the RF00816 family in RFAM shows miRNAs have been characterized for 14 different species (6 species of nematodes and 8 species or arthropods). Finally Cand9 had a hit against an undefined small ncRNA deposited in the smiRNAdb database. Cand9 also had a structural alignment to the RFAM family RF00708, but this alignment included only 12 bases of the original read, and in spite of the fact that the family has been characterized in human, the alignment did not cover the whole consensus sequence.
                        http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-14-735/MediaObjects/12864_2012_5450_Fig6_HTML.jpg
                        Figure 6

                        Structural alignment of candidates Cand4 and Cand9 against RFAM family RF00816. A: Secondary structure of family RF00816 and regions corresponding to the original reads of candidates Cand4 and Cand9; B: INFERNAL alignment of candidates Cand4 and Cand9 against family RF00816 - the original reads correspond to positions 100–130 of each candidate.

                        All candidates were also submitted to the software MatureBayes, which identifies mature sequences within putative precursors, but in none of the candidates the original read was considered as a mature miRNA. However, in this software, the definition of the precursor sequence could significantly impact the results. Our putative precursor sequences were based on the average length of precursors in humans. The extension of 100nt on both sides of the read derives from the fact that we did not have enough information to the position the read within its precursor. This may have affected the quality of the software’s predictions.

                        This conjunction of evidences point to the genomic region comprising candidates Cand4 and Cand9 as a strong candidate for a new miRNA gene. Other two candidates, Cand11 and Cand12, also presented consistent evidence, but based only on HMMIR, which is heavily dependent on the folding provided by the RNAfold algorithm, and suffers from the same dependence of the putative precursor sequence as MatureBayes.

                        The genomic region comprising the strongest candidate for a new miRNA was evaluated in clinical samples. Expression levels of the sequences that originated Cand4 and Cand9 were up-regulated in several tumor samples, when compared to cancer-free surgical margins, suggesting that this molecule might, indeed, have a role in HNSCC. Further experiments are necessary in order to confirm if the molecule is truly a miRNA.

                        Conclusions

                        MiRNA content of two cell models used for HNSCC cancer research was characterized by deep sequencing. Several miRNAs were equally expressed between a cancer cell line and keratinocytes, suggesting that the regulation of processes targeted by these molecules may be independent from the cancer phenotype. On the other hand, we provide evidences that pathways commonly deregulated in HNSCC, such as apoptosis and cell differentiation, may be targeted by miRNAs differentially expressed between cell types. Our results also imply that the use of cell models for microRNA functional studies in cancer research demands careful assessment of underlying molecular characteristics for proper data interpretation. The characterization of putative novel microRNA molecules carried out here revealed one strong new miRNA gene to be experimentally validated. This candidate was mostly expressed in the cancer cell line and its expression was validated in clinical samples, indicating a possible role in cancer.

                        Methods

                        Cell culture

                        We used the HNSCC cell line SCC25, derived from a SCC of the tongue for small RNA transcriptome sequencing. It was obtained from the American Type Culture Collection (ATCC catalog number CRL-1628). For the purpose of this study, the cell line was 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.

                        Human oral epithelial tissue was obtained from healthy volunteers undergoing dental surgeries, following previously published procedures [6], after approval by the Research Ethics Committee of IPEN under License Number 087/CEP-IPEN/SP and with informed consent signature. Keratinocytes used for small RNA transcriptome sequencing and gain-of-function experiments were grown on a fibroblast feeder-layer in a Dulbecco’s Modified Eagle Medium (DMEM; Gibco, New York, NY, USA) F-12 Nutrient Mixture (HAM, Gibco, New York, NY, USA) (2:1), with 10% Bovine Serum Product Fetal Clone III (Hyclone, Logan, Utah, USA), penicillin (100 U/ml), streptomycin (100 lg/ml), gentamicin (50 lg/ml), and amphotericin B (2.5 lg/ml) glutamine (4 mM), adenine (0.18 mM) (Sigma–Aldrich, St. Louis, MO, USA), insulin (5 lg/ml) (Sigma–Aldrich, St Louis, MO, USA), hydrocortisone (0.4 lg/ml) (Sigma–Aldrich, St Louis, MO, USA), cholera toxin (0.1 nM) (Sigma–Aldrich, St Louis, MO, USA), triiodotyronine (20 pM) (Sigma–Aldrich, St Louis, MO, USA) and epidermal growth factor (10 ng/ml) (R&D Systems, Minneapolis, MN, USA).

                        Clinical samples

                        Eight patients with oral squamous cell carcinoma (tongue and floor of the mouth) were selected for this study. The clinical and pathological profile of patients is shown in Table  6. 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 were snap-frozen in liquid nitrogen immediately after surgery. Analysis of hematoxylin and eosin-stained sections confirmed that >75% tumor cells in all HNSCC samples and that surgical margins were tumor-free.
                        Table 6

                        Clinical data of patients in this study

                        Patient

                        Tumor site

                        Gender

                        Age (yr)

                        196

                        OC-T

                        Male

                        49

                        240

                        OC-FOM

                        Male

                        75

                        277

                        OC-T

                        Male

                        66

                        296

                        OC-FOM

                        Male

                        53

                        306

                        OC-FOM

                        Female

                        82

                        321

                        OC-FOM

                        Male

                        69

                        333

                        OC-T

                        Male

                        52

                        349

                        OC-FOM

                        Male

                        59

                        OC-T: Oral Cavity – Tongue, OC-FOM: Oral Cavity - Floor of the Mouth.

                        Small RNA library construction and sequencing

                        Total RNA was obtained from confluent cell cultures or clinical samples using the mirVana Isolation Kit (Ambion Inc., USA). The concentration and quality were determined using a Nanovue spectrophotometer (GE Healthcare). Library construction followed, strictly, the SOLiD Total RNA-Seq Kit for Small RNA Libraries protocols (Ambion Inc., USA). Three libraries were constructed for each cell type. Eight hundred ng of total RNA were used as a template to obtain the small RNA library and we used the SOLiD RNA Barcoding System (Ambion Inc., USA) for library multiplexing. The SOLiD3 sequencing system (Life Technologies, CA, USA) was used to generate reads that were 35 bp long. Default parameters were used at all instances during sequencing.

                        Sequencing data analysis

                        Sequence analysis was performed using the Small RNA Analysis Tool (RNA2MAP) [25] using the following parameters: three color-space mismatches within the 'seed sequence’ (first 18 bases of the reads), and six color-space mismatches on the following positions of the 35 bp reads. All sequences that matched tRNA, rRNA, DNA repeats and adaptor molecules were filtered out and then the remaining reads were matched against miRNA precursor sequences in miRBase release 18. Then, using a Perl script we selected only reads containing mature sequences.

                        The RNA2MAP mapping tool generated two types of alignment: reads uniquely mapped to miRNAs and reads generating multiple hits. In order to select for molecules most likely to represent a mature miRNA we restricted multiple alignments to 5 hits and such hits should be within variations of a single miRNA family. For instance, reads mapping identically to hsa-mir-103a-1 and hsa-mir-103a-2 counted as "hsa-mir-103a". Reads generating multiple hits that did not conform to these parameters were discarded.

                        Following the identification of miRNAs present in each sample, we also used Perl scripts to compare normal keratinocytes and SCC25, and cancer and cancer-free patients.

                        To compare the expression data levels, the expression of each mature miRNA was normalized using the highest expression value in the dataset [26] and miRNAs presenting a difference in expression level of at least 2-fold were considered differentially expressed between datasets.

                        For data visualization of clinical samples, expression patterns of miRNA detected in every sample were clustered by Principal Components Analysis (PCA) using Partek Genomics Suite (v6.6).

                        The miRNA targets were searched using Ingenuity Pathway Analysis ( http://​www.​ingenuity.​com), through the integration of Tarbase and miRecords databases. For pathway mapping and Gene Ontology term enrichment analysis we used DAVID Bioinformatics Resources ( http://​david.​abcc.​ncifcrf.​gov).

                        Relative quantitation of gene expression by real-time PCR

                        For miRNA expression analysis, the cDNA was synthesized from 100 ng of total RNA using sequence-specific stem-loop primers for hsa-miRNA-1, hsa-miRNA-7, hsa-miRNA-21, hsa-miRNA-24, hsa-miRNA-133, hsa-miRNA-196 and for the endogenous control RNU48 (TaqMan miRNA RT kit, Life Technologies, Carlsbad, CA, USA). The relative quantitation of miRNA was carried out using the Taqman Universal PCR Master Mix (Life Technologies, Carlsbad, CA, USA), according to the manufacturer’s instructions. Data was normalized to the expression of RNU48 and analyzed using the delta delta Ct method.

                        Microarray analysis

                        Microarray analysis was performed using Agilent Whole Human Genome Microarray 4 × 44 K arrays and labeled using the One Color Quick Amp Labeling Kit (Agilent Technologies). A total of four samples were analyzed: two biological replicates of SCC25 and two biological replicates of normal keratinocytes. Hybridization and washing followed protocols described by the manufacturer (Agilent Technologies). The one-color arrays were scanned by GenePix 4000B Scanner (Axon), and analyzed using the Agilent Feature extraction software (version 9.5). The quality control of the microarrays was assessed using the standard Agilent controls to verify that the arrays met the expected criteria. The gProcessedSignal from each array was loaded into Partek Genomics Suite (v6.6) and normalized between arrays using quantile normalization. For subsequent statistical analysis we used the ANOVA implementation of Partek. Differences between cell lines and between clinical samples are presented in fold-changes. The raw data can be assessed at Gene Expression Omnibus under the accession number GSE41436.

                        Discovery of novel miRNA

                        All reads that did not match miRBase v18 were subjected to a bioinformatics pipeline represented in Figure  3. Reads were initially clustered based on their genomic positions and extended by 100 nucleotides upstream and downstream from their genomic coordinates. The size of the extension was based on an analysis of the sizes of known precursor human miRNAs deposited in miRBase v 16, which revealed that 98% of these miRNAs had size smaller than 135nt (Additional file 7). Including most of the hypothetical precursor sequence was important for the structural characterization of the candidates (INFERNAL, RNAFold and HHMMiR). The 448 sequences extracted in this manner formed a miRNA candidate set.

                        All 448 candidate miRNA sequences obtained as described above were then submitted to a 5-step analysis: similarity search against miRNA sequences, false positive detection, structural search against RFAM, ab initio classification using HHMMiR, folding using RNAfold (version 1.8) and tabulation of the results.

                        Similarity search against known miRNA sequences was performed against a local miRNA database with sequences from databases designated by NRDR [27] as containing miRNAs. Sequences from 15 different databases were included [28], miRBase [29], miRNAMap [30], microRNA.org [31], Argonaute, ASRP [32], CSRDB [33], fRNAdb [34], ncRNAdb [35], NONCODE [36], RFAM [37], RNAdb [38], smiRNAdb [39], TarBase [40], UCSC Genome Browser human miRNAs [41]. Similarity search was performed using WU-BLAST v.2.2.6 with the DUST filter turned on, wordsize 6, and minimum e-value of 0.0001. Results were filtered to the following cutoff values: minimum 80% coverage of the query or of the subject sequence, 80% minimum identity. False positive detection adopted the conservative approach of excluding all candidates mapping on known exons, all candidates that matched other ncRNA types on our local database and all candidates postulated as tRNAs by tRNAScan-SE version 1.23 [42].

                        Structural search against RFAM version 9.1 [37] was performed using INFERNAL version 0.81 with the recommended cutoff value 25. Ab initio prediction with HHMMiR version 1.2 was performed using the recommended cutoff value of 0.71 MLE.

                        Identification of mature miRNAs within putative precursors

                        In order to identify a mature miRNA sequence within putative miRNA precursors, following the analysis presented in Figure  3 candidates that showed the original read mapping to a predicted stem (at least partially) were analyzed using the MatureBayes algorithm [43]. This computational tool incorporates a Naive Bayes classifier to identify mature miRNA candidates based on sequence and secondary structure information of their miRNA putative precursors.

                        Declarations

                        Acknowledgements

                        The authors acknowledge the contribution of GENCAPO (Brazilian Head and Neck Genome Project) for clinical samples and for clinical and pathological data collection (complete list of members and affiliations presented at http://​www.​gencapo.​famerp.​br). This work was supported by FAPESP (grants 09/04166-5 and 10/51168-0) and by Hospital Israelita Albert Einstein.

                        Authors’ Affiliations

                        (1)
                        Albert Einstein Research and Education Institute, Hospital Israelita Albert Einstein
                        (2)
                        Instituto de Matemática e Estatística, University of Sao Paulo
                        (3)
                        Nuclear and Energy Research Institute IPEN/CNEN
                        (4)
                        Division of Head and Neck Surgery, Department of Surgery, School of Medicine, University of Sao Paulo
                        (5)
                        Department of Epidemiology, Faculty of Public Health, University of Sao Paulo
                        (6)
                        Department of Stomatology, Faculty of Dentistry, University of Sao Paulo
                        (7)
                        Federal University of Technology

                        References

                        1. Lee RC, Feinbaum RL, Ambros V: The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell 1993,75(5):843–854.PubMedView Article
                        2. Gkirtzou K, Tsamardinos I, Tsakalides P, Poirazi P: MatureBayes: a probabilistic algorithm for identifying the mature miRNA within novel precursors. PLoS One 2010,5(8):e11843.PubMedView Article
                        3. Tanzer A, Amemiya CT, Kim CB, Stadler PF: Evolution of microRNAs located within Hox gene clusters. J Exp Zool Part B 2005,304B(1):75–85.View Article
                        4. 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.PubMedView Article
                        5. 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.PubMedView Article
                        6. Kimura S, Naganuma S, Susuki D, Hirono Y, Yamaguchi A, Fujieda S, Sano K, Itoh H: Expression of microRNAs in squamous cell carcinoma of human head and neck and the esophagus: miR-205 and miR-21 are specific markers for HNSCC and ESCC. Oncol Rep 2010,23(6):1625–1633.PubMed
                        7. Yu J, Peng H, Ruan Q, Fatima A, Getsios S, Lavker RM: MicroRNA-205 promotes keratinocyte migration via the lipid phosphatase SHIP2. FASEB J 2010,24(10):3950–3959.PubMedView Article
                        8. Walker JC, Harland RM: microRNA-24a is required to repress apoptosis in the developing neural retina. Gene Dev 2009,23(9):1046–1051.PubMedView Article
                        9. Du L, Shen J, Weems A, Lu SL: Role of phosphatidylinositol-3-kinase pathway in head and neck squamous cell carcinoma. J Oncol 2012, 2012:450179.PubMedView Article
                        10. Nemunaitis J: Head and neck cancer: response to p53-based therapeutics. Head Neck 2011,33(1):131–134.PubMedView Article
                        11. Le MTN, Teh C, Shyh-Chang N, Xie HM, Zhou BY, Korzh V, Lodish HF, Lim B: MicroRNA-125b is a novel negative regulator of p53. Gene Dev 2009,23(7):862–876.PubMedView Article
                        12. Nakanishi H, Taccioli C, Palatini J, Fernandez-Cymering C, Cui R, Kim T, Volinia S, Croce CM: Loss of miR-125b-1 contributes to head and neck cancer development by dysregulating TACSTD2 and MAPK pathway. Oncogene 2013.
                        13. Zheng B, Liang L, Huang S, Zha R, Liu L, Jia D, Tian Q, Wang Q, Wang C, Long Z, et al.: MicroRNA-409 suppresses tumour cell invasion and metastasis by directly targeting radixin in gastric cancers. Oncogene 2012,31(42):4509–4516.PubMedView Article
                        14. Viticchie G, Lena AM, Latina A, Formosa A, Gregersen LH, Lund AH, Bernardini S, Mauriello A, Miano R, Spagnoli LG, et al.: MiR-203 controls proliferation, migration and invasive potential of prostate cancer cell lines. Cell Cycle 2011,10(7):1121–1131.PubMedView Article
                        15. Sonkoly E, Love J, Meisgen F, Wei T, Brodin P, Jaks V, Kasper M, Shimokawa T, Harada M, Heilborn J, Hedblad M-A, Hippe A, Grander D, Homey B, Zaphiropoulos PG, Arsenian-Henriksson M, Stahe M, Pivarcsi A: MicroRNA-203 functions as a tumor suppressor in basal cell carcinoma. Oncog 2012, 1–9.
                        16. Kozaki K, Imoto I, Mogi S, Omura K, Inazawa J: Exploration of tumor-suppressive microRNAs silenced by DNA hypermethylation in oral cancer. Cancer Res 2008,68(7):2094–2105.PubMedView Article
                        17. Jiang X, Huang H, Li ZJ, He CJ, Li YY, Chen P, Gurbuxani S, Arnovitz S, Hong GM, Price C, et al.: miR-495 is a tumor-suppressor microRNA down-regulated in MLL-rearranged leukemia. Proc Natl Acad Sci USA 2012,109(47):19397–19402.PubMedView Article
                        18. Lajer CB, Nielsen FC, Friis-Hansen L, Norrild B, Borup R, Garnaes E, Rossing M, Specht L, Therkildsen MH, Nauntofte B, et al.: Different miRNA signatures of oral and pharyngeal squamous cell carcinomas: a prospective translational study. Br J Cancer 2011,104(5):830–840.PubMedView Article
                        19. Gee HE, Camps C, Buffa FM, Patiar S, Winter SC, Betts G, Homer J, Corbridge R, Cox G, West CM, et al.: hsa-mir-210 is a marker of tumor hypoxia and a prognostic factor in head and neck cancer. Cancer 2010,116(9):2148–2158.PubMed
                        20. Wong TS, Liu XB, Wong BYH, Ng RWM, Yuen APW, Wei WI: Mature miR-184 as potential oncogenic microRNA of squamous cell carcinoma of tongue. Clin Cancer Res 2008,14(9):2588–2592.PubMedView Article
                        21. Chou SC, Azuma Y, Varia MA, Raleigh JA: Evidence that involucrin, a marker for differentiation, is oxygen regulated in human squamous cell carcinomas. Br J Cancer 2004,90(3):728–735.PubMedView Article
                        22. Kim YY, Lee EJ, Kim YK, Kim SM, Park JY, Myoung H, Kim MJ: Anti-cancer effects of celecoxib in head and neck carcinoma. Mol Cells 2010,29(2):185–194.PubMedView Article
                        23. Liang CH, Wang GH, Hung WJ, Lin RJ, Cheng DL, Chou TH: Apoptosis effect of Sinularia leptoclados, S-depressan and S-inflate extracts in human oral squamous cell carcinomas. J Taiwan Inst Chem E 2010,41(1):86–91.View Article
                        24. Fu X, Han Y, Wu Y, Zhu X, Lu X, Mao F, Wang X, He X, Zhao Y: Prognostic role of microRNA-21 in various carcinomas: a systematic review and meta-analysis. Eur J Clin Invest 2011,41(11):1245–1253.PubMedView Article
                        25. SOLiDTM 3 system application documentation small RNA analysis tool. Applied boisystems 2009.
                        26. Ribeiro-dos-Santos A, Khayat AS, Silva A, Alencar DO, Lobato J, Luz L, Pinheiro DG, Varuzza L, Assumpcao M, Assumpcao P, et al.: Ultra-deep sequencing reveals the microRNA expression pattern of the human stomach. PLoS One 2010,5(10):e13205.PubMedView Article
                        27. Paschoal AR, Maracaja-Coutinho V, Setubal JC, Simoes ZL, Verjovski-Almeida S, Durham AM: Non-coding transcription characterization and annotation: a guide and web resource for non-coding RNA databases. RNA Biol 2012.,9(3):
                        28. Maselli V, Di Bernardo D, Banfi S: CoGemiR: a comparative genomics microRNA database. BMC Genomics 2008, 9:457.PubMedView Article
                        29. Kozomara A, Griffiths-Jones S: miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res 2011, 39:D152-D157. Database issuePubMedView Article
                        30. Hsu SD, Chu CH, Tsou AP, Chen SJ, Chen HC, Hsu PW, Wong YH, Chen YH, Chen GH, Huang HD: miRNAMap 2.0: genomic maps of microRNAs in metazoan genomes. Nucleic Acids Res 2008, 36:D165-D169. Database issuePubMedView Article
                        31. Betel D, Wilson M, Gabow A, Marks DS, Sander C: The microRNA.org resource: targets and expression. Nucleic Acids Res 2008, 36:D149-D153. Database issuePubMedView Article
                        32. Backman TW, Sullivan CM, Cumbie JS, Miller ZA, Chapman EJ, Fahlgren N, Givan SA, Carrington JC, Kasschau KD: Update of ASRP: the Arabidopsis Small RNA Project database. Nucleic Acids Res 2008, 36:D982-D985. Database issuePubMedView Article
                        33. Johnson C, Bowman L, Adai AT, Vance V, Sundaresan V: CSRDB: a small RNA integrated database and browser resource for cereals. Nucleic Acids Res 2007, 35:D829-D833. Database issuePubMedView Article
                        34. Kin T, Yamada K, Terai G, Okida H, Yoshinari Y, Ono Y, Kojima A, Kimura Y, Komori T, Asai K: fRNAdb: a platform for mining/annotating functional RNA candidates from non-coding RNA sequences. Nucleic Acids Res 2007, 35:D145-D148. Database issuePubMedView Article
                        35. Szymanski M, Erdmann VA, Barciszewski J: Noncoding RNAs database (ncRNAdb). Nucleic Acids Res 2007, 35:D162-D164. Database issuePubMedView Article
                        36. Bu D, Yu K, Sun S, Xie C, Skogerbo G, Miao R, Xiao H, Liao Q, Luo H, Zhao G, et al.: NONCODE v3.0: integrative annotation of long noncoding RNAs. Nucleic Acids Res 2012, 40:D210-D215. Database issuePubMedView Article
                        37. Griffiths-Jones S, et al.: Annotating non-coding RNAs with Rfam. In Current protocols in bioinformatics Edited by: Baxevanis AD. 2005. Chapter 12:Unit 12 15
                        38. Pang KC, Stephen S, Dinger ME, Engstrom PG, Lenhard B, Mattick JS: RNAdb 2.0--an expanded database of mammalian non-coding RNAs. Nucleic Acids Res 2007, 35:D178-D182. Database issuePubMedView Article
                        39. 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.PubMedView Article
                        40. Vergoulis T, Vlachos IS, Alexiou P, Georgakilas G, Maragkakis M, Reczko M, Gerangelos S, Koziris N, Dalamagas T, Hatzigeorgiou AG: TarBase 6.0: capturing the exponential growth of miRNA targets with experimental support. Nucleic Acids Res 2012, 40:D222-D229. Database issuePubMedView Article
                        41. Fujita PA, Rhead B, Zweig AS, Hinrichs AS, Karolchik D, Cline MS, Goldman M, Barber GP, Clawson H, Coelho A, et al.: The UCSC Genome Browser database: update 2011. Nucleic Acids Res 2011, 39:D876-D882. Database issuePubMedView Article
                        42. Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res 1997,25(5):955–964.PubMed
                        43. Jun J, Oh KM: Asian and hispanic Americans’ cancer fatalism and colon cancer screening. Am J Health Behav 2013,37(2):145–154.PubMedView Article

                        Copyright

                        © Severino et al.; licensee BioMed Central Ltd. 2013

                        This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License ( http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

                        Advertisement