MicroRNA-21 , 204 and 125 b Play Potential Roles in Tumorigenesis of Melanoma

Aim: The high mortality rate of melanoma is due to partially the lack of good diagnostic markers and treatment strategies. Over the past several years, several microRNA (miRNA) profiling studies have been performed on melanoma tissues, but with extremely inconsistency, the diagnostic value of miRNA candidates in melanoma remains under debate. Thus, this study aims to systematically evaluate the consistency of miRNAs tissue in multiple independent studies in melanoma. Method: Eligible studies were screened and selected from the PubMed, EMBASE, and Web of Science. A systematic analysis of published miRNA expression studies that compared the miRNA expression profiles between melanoma tissues and normal skin tissue was conducted. A vote-counting strategy was followed with the collection of information. Real time PCRs were employed to validate miRNA candidates with high consistency. Targets of consistent miRNAs were predicted by online programs (like miRTarBase, microRNA.org and TargetScanHuman 6.2). Enrichment analyses for gene ontology (GO) terms and Kyoto encyclopedia of genes and genomes (KEGG) pathways were carried out with Database for Annotation, Visualization, and Integrated Discovery (DAVID). Results: A total of 303 differentially expressed miRNAs were reported in the 10 miRNA-profiling studies during comparison of melanoma tissues with normal tissues; 132 were up-regulated in melanoma, and 171 were down-regulated. However, in the group of consistently reported miRNAs (cutoff > 3 times), only moderate numbers of consistent and differentially expressed miRNAs were selected. miRNA-21 was found increased in 5 different studies, miRNA-146b, miRNA-17 and miRNA-18a were reported up-regulated in 4 profiling studies. Meanwhile, miRNA-204 and miRNA125b were found down-regulated in 5 studies, miRNA-141, miRNA-149, miRNA-224, miRNA-200b, miRNA-200c were consistently decreased in just 4 out of 10 profiling studies in total. The directions of differential expression of these miRNA candidates were confirmed by real time PCRs. Enrichment analyses demonstrated that programmed cell death and transcription regulation Chengtan Li and Xiayu Wang contribute to the manuscript equally. Corresponding author.


Introduction
Cutaneous melanoma, as currently the sixth most common cancer in white men and women in the United States, is the most aggressive form of skin cancer characterized by poor prognosis [1], with the median no more than 1 year of survival time and is less than 5% of 5-year survival rate [2] [3].The incidence and mortality of melanoma continues to increase faster than any other cancer worldwide and there is no curative therapy for advanced stages of melanoma [4] [5]; as an example, almost 70 thousand new cases of melanoma and around 9 thousand deaths were predicted in 2009 [6].Unfortunately, current treatment options are limited for metastatic melanoma.Therefore, new indicators of prognosis and therapeutic targets are pressing required.Genetic factors and exposure to ultraviolet radiation are risk factors for melanoma pathogenesis [7].
Since discovered in 1993, microRNAs (miRNAs), as small, evolutionarily conserved, single-stranded and non-coding RNA molecules, have been found to be important to regulate the expression of up to 30% encoding genes by binding to specific mRNA targets and promoting their degradation and/or translational inhibition [8], implying their important roles in physiological function and pathological condition.Growing evidence suggests that miRNAs, besides for a variety of essential cellular processes including development, cell differentiation, inflammation, proliferation, and apoptosis [9], play very important roles in initiating tumorigenesis [10], functioning as tumor suppressors or oncogenes (oncomir), evidenced by the miRNA deregulation in various cancer types like lung cancer [11].Thereafter, miRNAs are implied as biomarkers and intervention targets against various forms of cancer.Further, miRNAs can be a superior method for cancer subtype classification and prognostication [12] [13].An important feature of miRNAs is their remarkable stability; they can be well preserved in tissue samples even after years of formalin-fixation and paraffin embedding, and can be efficiently extracted from such specimens and variety of cell lines [14].Several studies have been conducted to search for biomarkers or therapeutic agents by identifying the differential expression of miRNAs between melanoma tissue samples and corresponding non-tumor skin tissue from the same patient [15]- [25].These studies generated hundreds of differentially expressed miRNAs, however, many of these are likely to be false positives and finally only a small number of them may have clinical significance and act as diagnostic and prognostic biomarkers.Further, different profiling studies showed quite inconsistency in the identified differentially expressed miRNAs.Third, in miRNA-based therapeutics, miRNAs are tissue specific; their role depends on cellular context.For instance, overexpression of miR-143 in pancreatic xenografts reduces tumor growth, and anti-miR-143 is required to suppress skin tumor formation [26].It is therefore crucial to identify skin tissue specific and consistently differentially expressed miRNAs to investigate the potential as biomarkers or therapeutic agents.Thus, a logical approach to distinguish important miRNAs from a large number of candidate miRNA lists is to search for the intersection of miRNAs identified in multiple independent studies [27].Although this kind of methods has become increasingly popular, no published study has identified the intersections of melanoma related miRNAs based on a large number of miRNA expression profiling studies thus far.Abnormal activation and inactivation of oncogenes and tumor suppressor genes are important factors leading to malignancy (including melanoma).Clarifying the molecular mechanisms of melanoma could provide a basis for melanoma risk assessment, early diagnosis, effective treatment and intervention.
We conducted this systematic analysis to identify the most important differentially expressed miRNAs that had been consistently reported in a series of independent miRNA expression profiling studies in melanoma patients.Moreover, we further validated some of the miRNAs that were most up-or down-regulated using realtime PCR in 6 pairs of melanoma and matched adjacent non-tumor tissue samples.

Search Strategy
Potential studies published in English were collected from PubMed, EMBASE, Web of Science published from June, 2006 to Feb 2015 (last accessed on April 15, 2015) using the following medical subject headings terms: "miRNA" OR "microRNA" OR "miR", "skin cancer" OR "melanoma", "profiling" OR "microarray".

Inclusion Criteria of the Literature
For a study to be included in this systematic review, several criteria had to be met: 1) studies had to be miRNA profiling studies in melanoma patients; 2) studies had to use melanoma tissues and their corresponding adjacent non-tumor tissues for comparison; 3) miRNA microarray and sequencing methods were comprised; they reported cut-off criteria of differentially expressed miRNAs, and 4) they reported the validation method set.Therefore, miRNA profiling studies that used serum samples of melanoma patients or skin relevant cell lines, those that compared melanoma biopsies from tumors with different stages of disease, or those using different miRNA techniques were excluded.Review articles were also excluded.

Data Extraction
Two investigators (X.W. and Y.C.) independently evaluated and extracted the data using standard protocols, and all discrepancies were resolved by the corresponding investigator (Y.Y.).From the full text and corresponding supplementary information, the following eligibility items were collected and recorded for each study: author, journal and year of publication, location of study, methods and characteristics of animal modeling, platform of miRNA expression profiling, author defined cut-off criteria of statistically significant differentially expressed miRNAs.

Ranking
Each included profiling/microarray study [15]- [25] provided a list of differentially expressed miRNAs related to melanoma.For example, Poliseno et al. [24] carried out comprehensive analysis of melanoma-related miRNAs based on a large number of miRNA profiling studies.Griffith et al. devised a method to rank potential molecular biomarkers for comparison groups [28], which had been used for miRNA profiling like meta-analysis studies.Briefly, miRNAs were ranked to the criteria in the following order of importance: (i) the miRNA was consistently reported as differentially expressed in a consistent direction of change; (ii) the frequency of the miRNAs was reported in the microarray studies; (iii) the total sample size for each consistent reported miRNAs.

Validation of the miRNAs Using Real Time PCR
To validate the profiling results, 6 fresh melanoma tissues and their paired non-tumor skin tissues were obtained from the Second Hospital of Hangzhou City, affiliated to Hangzhou Normal University School of Medicine.Total RNA was extracted from 6 pairs of matched human melanoma specimens (including cancer and adjacent noncancerous tissues) using TRIzol reagent (Life technology, Casbad, CA, USA) according to the manufacturer's instructions.The differentially expressed amount of the miRNAs was validated in triplicate by real time PCR.Briefly, Reverse transcription from 3 ug RNA was done using SuperScript III First-Strand Synthesis System (Life technology, Casbad, CA, USA) according to the manufacturer's protocol.Real-time PCR was performed using iQ SYBR Green Supermix kit (Bio-Rad, Hercules, CA) with the iCycler sequence detection system (Bio-Rad) with specific primers (Table 1).Small nucleolar RNA U6 was used as an internal standard for normalization.The cycle threshold (C T ) value was calculated.The 2 −ΔCT (ΔC T = C TmiR − C TU6 RNA ) method was used to quantify the relative amount of miRNAs.

Identify the Experimentally Validated microRNA Target Genes
miRTarBase (mirtarbase.mbc.nctu.edu.tw),microRNA.organd TargetScanHuman 6.2 were used to explore all the target genes of miRNAs, and highlighted the genes that are experimentally validated by luciferase reporter assay, western blot, or microarray experiments with over expression or knock-down of miRNAs.

Enrichment Analysis
Enrichment analyses for gene ontology (GO) terms and Kyoto encyclopedia of genes and genomes (KEGG) pathways were carried out with Database for Annotation, Visualization, and Integrated Discovery (DAVID) [29].DAVID is a web-accessible program that integrates functional genomics annotations with intuitive graphical summaries.List of gene or protein identifiers are rapidly annotated and summarized according to shared categorical data for Gene Ontology, protein domain, and biochemical pathway.For this, we listed all the target genes of all the differentially expressed miRNAs, including up-regulated and down-regulated miRNAs.The top 17 or 18 GO terms and KEGG pathways showing association with target genes were listed with GO terms, KEGG pathway, number of genes in the GO term, number of genes in the KEGG pathways.

Statistical Analysis
Student's t-test was used to compare values between two independent groups.

Included Independent Studies
A total of 68 studies were recorded using PubMed, EMBASE and web of science.38 of which were excluded after screening the titles and abstracts, and 20 studies were excluded after reading the full text based on the inclusion and exclusion criteria, only 10 independent studies were included in this systematic analysis.The detailed workflow used in our analysis was shown in Figure 1.The detailed characteristics of each study are given in Table 2 in descending order of their year of publication.

Differentially Expressed miRNAs
These 10 studies from 10 different groups with different platforms, different number of samples, and various statistical analyses had been employed for microRNA profiling analysis to compare melanoma tissue with corresponding noncancerous skin tissue.The number of differentially expressed microRNAs ranges from 4 to 53.A total of 303 differentially expressed miRNAs were reported in the 10 profiling studies.171 microRNAs are downregulated in melanoma, and 132 microRNAs are up-regulated.Among the 171 down-regulated miRNAs, two of them were reported in 5 microarray studies (miRNA-204 and miRNA-125b), five miRNAs were reported in four profiling studies (miRNA-141, miRNA-149, miRNA-224, miRNA-200b, miRNA-200c).Among the 132 upregulated miRNAs, there was only one different miRNA in 5 studies (miRNA-21), three studies were reported  in 4 profiling studies (miRNA-146b, miRNA-17 and miRNA-18a).Here, consistently expressed miRNAs, the corresponding microarray study and the total number of tissue samples were shown in Table 3.We also noticed that these selected miRNAs differentially expressed in a very inconsistent manner in 10 miRNA profiling studies.For example, in the groups of inconsistently reported miRNAs, miRNA-211 was reported down-regulated in three studies but up-regulated in two studies; miRNA-199a was increased in one study but down-regulated in 1 study (Table 4), respectively.

Primary miRNAs Analysis
The identified and selected miRNAs all play specific roles in skin homeostasis, or pathogenesis of melanoma.Particularly, all these miRNAs are either oncogenic miRNAs or tumor suppressors.Specifically, miRNA-21, miRNA-146b and miRNA-18a are oncogenic miRNAs (oncomir), and miRNA-204, miRNA-125b, miRNA-141, miRNA-200b and miRNA-200care tumor suppressors.We need to point out that these selected oncomires and tumor suppressors are differentially expressed only in moderate number of profiling studies (less than 5 or 4 times), implying that miRNA plays marginal roles in the tumorigeneisis and development of melanoma.

Target Genes of Differentially Expressed microRNAs
After we identified the most consistently differentially expressed microRNAs that were either down-regulated  or up-regulated in melanoma, we screened the their potential target genes with programs such as miRTarBase (mirtarbase.mbc.nctu.edu.tw),microRNA.organd TargetScanHuman 6.2, focusing on the genes which were confirmed either by real time PCR, western blot, microarray or luciferase assay.As a result, we identified 211 and 318 target genes (519 in total, some target genes are shared by two different kinds of miRNAs, reasons are to be defined) corresponding to those down-regulated and up-regulated miRNAs in melanoma (Table 5 and Table 6), respectively.

Discussion
In the current study, we established the concept, by conducting a systematic analysis, which only a few of consistently differentially expressed miRNAs may be of significance for tumorigenesis and development of melanoma.The potential targets of these selected miRNAs conduct the signaling pathways critical for programmed cell death and the regulation of transcription.Together, this research shed some light on the potential of certain miRNAs as biomarkers and therapeutic target for melanoma.Many conventional tumor biomarkers, such as ovarian cancer biomarker CA-125 and melanoma biomarker S-100, have encountered some limitation like poor sensitivity of specificity, for example, only 50% of patients with early stage ovarian cancer have elevated levels of CA125; specificity of S-100 for melanoma is about 75% [30].Thus, there is a critical need for reliable biomarkers that allow diagnosis with high sensitivity and precise monitoring the pathogenic process and therapeutic alternatives of melanoma.Many miRNAs are deregulated significantly even in early stages of tumor development and involved in tumorigenesis and metastasis of many caner categories, exhibiting their diagnostic and therapeutic potential utility.However, as many miRNA profiling studies on melanoma were performed on different platforms and very different processing methods by different investigators, there lacks the common concert among these studies, consequently, which ended within consistent miRNAs expression profile [15]- [25].
Based on our systematic analysis, these oncomirs and tumor suppressors were consistently and differentially expressed in relative less investigations (no more than 5 profiling studies), suggesting that miRNAs may only    As an oncomir, miRNA-21 negatively regulates the expression of many different tumor suppressor genes, such as PTEN [31], Bcl-2 [32] and Il-12p35 [33] etc, by which is involved in the tumorigenesis of different cancers like colon [34], liver [31] and lung [35].The involvement of miR-17 in pathogenesis of cancer is controversial and open to debate; miR-17 is up-regulated as oncomir in breast cancer in some groups [36]- [38], but act as tumor suppressor in breast cancer in another group [39], specifically, miR-17 inhibits tumor growth by stimulating T-cell mediated host immune response in melanoma [40].In our systematic study, miR-125bperforms tumor-suppressor functions in several  cancers including melanoma, by suppressing melanoma cell proliferation and metastasis.miR-125b is downregulated in metastatic melanoma and controls melanoma progression by direct regulation of c-Jun protein expression [18] [41].As a tumor suppressor, miR-204 was decreased in melanoma [42] and can inhibit cell proliferation by interfering with target gene IGFBP5 or RAB22A [43] [44].miRNAs function as oncogenes or tumor suppressors by involving multiple molecular mechanisms.By enrichment analysis, we found these oncomirs and tumor suppressors can targets variety of substrate genes, it is of great significance to point out that these targets concentrate on the cancer-related cellular process, mainly including apoptosis, which is in concert with the fact that apoptosis plays checkpoint roles in tumorigenesis of any cancers [45] [46], which further highlight the importance of these miRNAs as potential biomarkers or therapeutic alternatives.KEGG pathways analysis by DAVID demonstrated that the target genes of these miRNAs Figure 9. KEGG pathway (pathways in cancer) analysis with predicated target genes from all selected differentially expressed miRNAs.These pathways in cancer are a combination of many different but relevant pathways, such as MAPK signaling pathway, cell cycle signaling pathway and PI3K-Akt signaling pathway etc. Green boxes with red star represent the potential genes targeted by selected miRNAs.built up a very rich networks, combined with different pathways and molecule-molecule interactions.Wnt signaling pathway is very important for the integrity of endothelial cells and for melanoma tissue invasion and metastasis, inhibition of mTOR signaling pathway can reduce cellular apoptosis; p53 signaling pathway contributes to inhibit apoptosis; MAPK signaling pathway are crucial for cell proliferation; All these pathways work synergistically to evade apoptosis and increase proliferation.Some other pathways related to block of differentiation, resistance to chemotherapy, failed repair of genes, insensitivity to anti-growth signals are also implied.All these miRNA/target genes related pathways mimic exactly the pathways involved in the tumorigenesis of melanoma, which make these miRNAs even more important as potential biomarkers and therapeutic candidates for clinical application.Extensive studies should focus these miRNAs to establish them as biomarkers for cancer detection and progression.

Limitations
Several limitations of this research should be considered when interpreting the results due to some unsolvable reasons.Firstly, our literature searching was depended on English databases only, as a result, language bias may present.Secondly, our study did not include all the populations (only Chinese, British, American, German, Canadian and Luxemburg), so the result may not be able to apply to other populations such as Latin American and African.

Conclusion
We performed comprehensive literature search in multiple databases by limiting publication language and date.By systematic analysis, we filtered the most consistently expressed 4 up-regulated and 7 down-regulated miR-NAs, and further confirmed by real time PCR, specially, highlight would be given to miRNA-21, miRNA-204 and miR-125b since they consistently differentially expressed the most (5 studies respectively).By DAVID analysis, we showed that the target genes were mainly involved in the regulation of programmed cell death mediated by combined signaling pathways such as cell cycle signaling pathway.In short, these selected miRNAs

Figure 1 .
Figure 1.Flowchart for this systematic analysis.

Figure 3 .
Figure 3. Real-time PCR analysis of some most consistent up-regulated expressed miRNAs.Expression levels of miR-17, miR-146b, miR-21 and miR-18a in melanoma and adjacent noncancerous tissue samples were analyzed by real time PCR, U6 here functions as an internal control.The expressions of miR-21, miR-146b and miR-18a were significantly higher in melanoma tissues, while the expression of miR-17 was significantly reduced.Here, values represent means ± S.E. of three determinations.*, p < 0.05; **, p < 0.01.

Figure 4 .
Figure 4. GO analysis of target genes from most up-regulated miRNAs in selected profiling studies.Predicted target genes are mainly involved in cancer related cellular process including programmed cell death and metabolic process; further, these involvement are probably through the regulation of transcription regulation.

Figure 5 .
Figure5.KEGG pathway (pathways in cancer) analysis with predicated target genes from up-regulated miRNAs.These pathways in cancer are a combination of many different but relevant pathways, such as MAPK signaling pathway and cell cycle signaling pathway, etc. Green boxes with red star represent the potential genes targeted by selected miRNAs.

Figure 6 .
Figure 6.GO analysis of target genes from most down-regulated miRNAs in selected profiling studies.Predicted target genes are mainly involved in cancer related cellular process including programmed cell death, metabolic process; further, these involvements are probably through the regulation of transcription regulation.play marginal roles in the tumorigenesis and development of melanoma, which is concert with the fact that environmental factors like UV plays important roles in the onset of melanoma.Here, we found that three of the selected miRNAs (miRNA-21, miRNA-17, miRNA-146b and miRNA-18a) are enhanced, displaying oncogenic characteristics (oncomir), and seven of screened miRNAs (miRNA-204, miRNA-125b, miRNA-141, miRNA-200b and miRNA-200c) decreased markedly, functioning as tumor suppressors.As an oncomir, miRNA-21 negatively regulates the expression of many different tumor suppressor genes, such as PTEN[31], Bcl-2[32] and Il-12p35[33] etc, by which is involved in the tumorigenesis of different cancers like colon[34], liver[31] and lung[35].The involvement of miR-17 in pathogenesis of cancer is controversial and open to debate; miR-17 is up-regulated as oncomir in breast cancer in some groups[36]-[38], but act as tumor suppressor in breast cancer in another group[39], specifically, miR-17 inhibits tumor growth by stimulating T-cell mediated host immune response in melanoma[40].In our systematic study, miR-125bperforms tumor-suppressor functions in several

Figure 7 .
Figure 7. KEGG pathway (pathways in cancer) analysis with predicated target genes from down-regulated miRNAs.These pathways in cancer are a combination of many different but relevant pathways, such as MAPK signaling pathway, cell cycle signaling pathway, PI3K-Akt signaling pathway etc. Green boxes with red star represent the potential genes targeted by selected miRNAs.

Figure 8 .
Figure 8. GO analysis of target genes from all the differentially expressed miRNAs in selected profiling studies.Predicted target genes are mainly involved in cancer related cellular process including programmed cell death, which is performed probably though cell transcription regulation.

Table 1 .
Primers used for real time PCR.

Table 2 .
Characteristics of miRNAs datasets in hepatocellular carcinoma.

Table 5 .
Selected up-regulated miRNAs and their predicated target genes.

Table 6 .
Selected down-regulated miRNAs and their predicated target genes.