Dysregulation of microRNA in cholangiocarcinoma identified through a meta-analysis of microRNA profiling.

BACKGROUND: In the past decades, the potential of microRNA (miRNA) in cancer diagnostics and prognostics has gained a lot of interests. In this study, a meta-analysis was conducted upon the pooled miRNA microarray data of cholangiocarcinoma (CCA). AIM: To identify differentially expressed (DE) miRNAs and perform functional analyses in order to gain insights to understanding miRNA-target interactions involved in tumorigenesis pathways of CCA.
METHODS: Raw data from 8 CCA miRNA microarray datasets, consisting of 443 samples in total, were integrated and statistically analyzed to identify DE miRNAs via comparison of levels of miRNA expression between CCA and normal bile duct samples using t-tests (P < 0.001). The 10-fold cross validation was performed in order to increase the robustness of the t-test results.
RESULTS: Our data showed 70 up-regulated and 48 down-regulated miRNAs in CCA. Gene Ontology and pathway enrichment analyses revealed that mRNA targets of DE miRNAs were significantly involved in several biological processes. The most prominent dysregulated pathways included phosphatidylinositol-3 kinases/Akt, mitogen-activated protein kinase and Ras signaling pathways.
CONCLUSION: DE miRNAs found in our meta-analysis revealed dysregulation in major cancer pathways involved in the development of CCA. These results indicated the necessity of understanding the miRNA-target interactions and the significance of dysregulated miRNAs in terms of diagnostics and prognostics of cancers. ©The Author(s) 2020. Published by Baishideng Publishing Group Inc. All rights reserved.

Core tip: At present, there is an accumulating mass of cholangiocarcinoma microRNA (miRNA) profiling data, however, it is challenging to gain the maximal information from these data because the experimental designs in each study tend to focus on only a few specific research questions. This work therefore integrates and inter-validates the cholangiocarcinoma miRNA expression profiles from multiple independent datasets to identify the differential dysregulation of miRNA and their corresponding downstream pathways underlying mechanism of pathogenesis. The significant merit of our findings offers a valuable reference for future studies and further investigation of these miRNA/genes and their interactions will eventually lead to the identification of genes and pathways important to the overall mechanism of the dysregulated processes in cholangiocarcinoma development.

INTRODUCTION

Cholangiocarcinoma (CCA), first described by Durand-Fardel in 1840, is a form of malignant tumor that originate from biliary epithelial cells in the liver and/or extrahepatic bile ducts[1,2]. CCA accounts for 10%-15% of hepato-biliary neoplasm. Thus, it is the second most common primary hepato-biliary malignancy after hepatocellular carcinoma (HCC)[3]. Moreover, the incidence and mortality rate of CCA have been reportedly increasing worldwide over the past three decades[2,4-6]. However, the prevalence of CCA vary greatly among different geographical regions of the world. Incidence of CCA in most Western countries ranges from 2 to 6 cases per 100000 people per year[7]. There is a higher prevalence of CCA in Asia and in people of Asian descent, which has been attributed to endemic chronic parasitic infestation[6,7]. Primary sclerosing cholangitis, inflammation that causes scars within the bile ducts, is the most common known predisposing factor for CCA[6]. In East and Southeast Asia, where the disease is common, CCA has been pathogenically associated with liver fluke infestation, particularly the endemic Clonorchis sinensis[8] and Opisthorcis viverrini[9]. In addition, hepatitis C virus infection and liver cirrhosis have been suggested as potential risk factors for CCA[6].

MicroRNAs (miRNAs) are a family of endogenous, non-coding RNAs found in plants, animals, and some viruses[10-13]. miRNA genes are highly-conserved and may be located either within the introns or exons of protein-coding genes (70%) or in intergenic areas (30%)[10,14]. A miRNA in its single-stranded functional form is usually 21-22 nucleotides long (though it can vary from 19-25 nucleotides)[10,13]. In present, the latest release of miRBase (version 22) (http://www.mirbase.org/) contains 38589 hairpin precursors and 48860 mature miRNA from 271 organisms[15]. Several hundreds of miRNA genes in the human genome have been discovered in the last decade[10,13,16]. However, it is estimated that the human genome may encode over 1000 miRNAs in total[17,18]. The primary function of miRNAs is to control gene expression at post-transcriptional level. miRNAs suppress the target mRNA expression, mostly through interaction with the 3’ untranslated region[13,19], resulting in inhibition of target mRNA translation activity and, to a lesser extent, targeting mRNA cleavage[11,16,20]. Each miRNA may be responsible for regulation of the expression of hundreds of gene targets[19].

Although the functions of dysregulated miRNAs in human cancers remain largely a mystery, multiple miRNAs and their corresponding target genes have been reported to be associated with tumor initiation and progression[18,21,22]. Many transcriptional profiling data demonstrated that miRNA expression profiling efficiently classified different tumor types more reliably than did mRNA profiling[23,24]. Systematic expression analyses using miRNA microarray technology have been performed in several types of cancer; for example, hepatic[25], colorectal[26], lung[27], and breast cancers[28].

However, due to insufficient control of false positives and the small sample sizes relative to the large sets of microarray probes, individual microarray-based studies are often deficient in terms of statistical robustness[29,30]. Meta-analysis is the use of statistical techniques to combine results from independent but related studies, hence it is one of the most preferable ways to increase the statistical power of the readily available microarray data. Moreover, it is relatively inexpensive in terms of financial and time investments[31].

Our meta-analysis of miRNA microarray datasets was aimed to identify the differentially expressed (DE) miRNA in various CCA samples compared to the non-cancerous counterparts. The results from the robust statistical tests would provide insights to understanding the regulatory potential of miRNAs in tumorigenesis pathways of CCA.

MATERIALS AND METHODS

miRNA microarray data collection

The miRNA microarray datasets were retrieved from the public repository database Gene Expression Omnibus (http: //ncbi.nlm.nih.gov/geo) via the computerized search using combinations of relevant keywords, including (miRNA OR microRNA) AND (cholangiocarcinoma OR CCA OR CCC). All dataset search hits were initially checked whether raw data were provided. Datasets without raw data were promptly omitted. A matrix series tables and raw data package(s) of each available dataset were downloaded and extracted. Eight miRNA microarray datasets from independent research studies were employed in our meta-analysis. Six out of these 8 datasets (GSE32958, GSE47764, GSE50894, GSE51429, GSE53870, GSE53992) were conducted using biopsied tissue samples. One study (GSE59856) used serum miRNA and another (GSE47396) used miRNA samples from cell lines.

Data processing

After extraction and background-correction, the GPR or TXT raw data files of each dataset were converted to comma-separated values file format and then imported as a data frame into R software version 3.1.2[32]. Before pooling samples of all datasets, the identification name of every miRNA probe was checked against identification entries registered in miRBase database (www.mirbase.org/)[15] and was subsequently renamed in accordance with miRBase identification, in order to avoid confusion from naming system of different microarray platforms. Log2 transformation was applied to all intensity values. The transformed data were then normalized following 2 steps, i.e., within-array normalization using a median centering method, followed by between-array normalization using a quantile normalization method. Hence normalized data would exhibit normal distribution with the same standard deviation across datasets which is an essential assumption of parametric statistical tests.

Statistical analysis

In order to identify DE miRNA in CCA, the pooled samples were grouped into 2 groups: “CCA” and “Normal”. The normalized intensity values of each miRNA were compared between these two groups using t-tests conducted in MultiViewer Experiment version 4.6 in TM4 software suite[33]. P value threshold was set at below 0.001 (P < 0.001). As clinical validation in our study is restricted, k-fold cross validation was applied to verify the integrity of significant DE miRNAs from t-test analysis. The 10-fold cross validation was performed by randomly assigning samples into 10 different sets and repeating t-test statistical analysis with the same parameter (P < 0.001) for 10 rounds. In each round, one set was chosen as the validation set whereas the rest were test sets. DE miRNA which showed statistical significance among test sets and validation sets from all 10 rounds were collected as validated DE miRNAs for further analysis.

Bioinformatics analysis and visualization of miRNA-target interactions

In order to assess the biological functions of the gene targets of DE miRNAs, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were employed using DIANA miRPath version 3.0[34]. The lists of up- and down-regulated miRNAs were categorized and separately uploaded as inputs to DIANA miRPath. The human KEGG and GO analyses were selected with P-value threshold at 0.01. The KEGG pathways, which were the most enriched by target genes of DE miRNAs, were selected to analyze the miRNA-target interactions. Unique targets of each DE miRNA in the selected pathways were filtered into the list, which was then used in the network visualizer software. Visualized miRNA-target interaction networks were constructed using Cytoscape version 3.3.0[35] with CyTargetLinker plugin version 2.1[36].

RESULTS

Differentially-expressed miRNA

A framework of our meta-analysis approach is depicted in Figure 1. In total, there were eligible 246 CCA and 197 normal samples from 8 independent miRNA microarray datasets, which included tissue samples, sera and cell lines as detailed in Table 1. Initial t-tests results identified 224 DE miRNAs consisting 114 up-regulated and 110 down-regulated miRNAs. Following statistical validation of these results using k-fold cross validation, the numbers of up-regulated and down-regulated miRNAs were confined to 70 (Table 2) and 48 (Table 3), respectively. The overall expression profiling of DE miRNAs is presented in Figure 2.

Figure 1

Overview of a meta-analysis approach in this study. CCA: Cholangiocarcinoma; DE: Differentially expressed; miRNA: MicroRNA.

Table 1

Summary of microRNA microarray datasets used in this study

Dataset	Sample type	Number of samples		Microarray platform
		CCA	Normal
GSE32958	Tissue	19	5	Nanostring nCounter human miRNA expression
GSE47396	Cell lines	2		miRCURY LNA miRNA Array, 5th generation
GSE47764	Tissue	3	3	Agilent-021827 human miRNA microarray
GSE50894	Tissue	16		Exiqon miRCURY LNA miRNA array version 10.0 extended and version 11.0 (condensed version)
GSE51429	Tissue	29		Applied Biosystems human miRNA array version 2.0
GSE53870	Tissue	63	9	State Key Laboratory human miRNA array 1104
GSE53992	Tissue	16	30	Agilent-031181 unrestricted human miRNA version 16.0
GSE59856	Serum	98	150	3D-Gene human miRNA version 20_1.0.0
Total		246	197

CCA: Cholangiocarcinoma; miRNA: MicroRNA.

Table 2

List of 10-fold-cross-validated differentially-expressed microRNAs which were up-regulated in cholangiocarcinoma

miRNA	FC	Abs. t value	DF	Adj. P value	FDR
hsa-let-7b-3p	1.482	6.603	74	0.000	0.000
hsa-miR-1193	1.579	7.409	112	0.000	0.000
hsa-miR-1207-3p	1.516	10.533	120	0.000	0.000
hsa-miR-1224-3p	5.164	7.024	66	0.000	0.000
hsa-miR-1227-3p	2.039	7.963	117	0.000	0.000
hsa-miR-1237-3p	1.594	5.7	355	0.000	0.000
hsa-miR-1247-5p	1.781	7.406	126	0.000	0.000
hsa-miR-1249-3p	3.877	9.138	137	0.000	0.000
hsa-miR-1267	2.433	7.758	119	0.000	0.000
hsa-miR-1269a	1.647	7.209	72	0.000	0.000
hsa-miR-1273d	1.665	7.159	80	0.000	0.000
hsa-miR-1284	1.606	7.019	127	0.000	0.000
hsa-miR-1304-5p	1.651	6.863	89	0.000	0.000
hsa-miR-1322	1.731	9.986	132	0.000	0.000
hsa-miR-147b	2.727	8.195	52	0.000	0.000
hsa-miR-149-3p	3.022	8.11	88	0.000	0.000
hsa-miR-181a-2-3p	1.476	7.004	90	0.000	0.000
hsa-miR-185-3p	5.044	9.684	113	0.000	0.000
hsa-miR-1908-5p	9.384	7.977	150	0.000	0.000
hsa-miR-1909-3p	2.19	6.502	129	0.000	0.000
hsa-miR-1913	2.222	6.861	137	0.000	0.000
hsa-miR-1914-5p	1.575	6.615	140	0.000	0.000
hsa-miR-1972	2.043	8.605	111	0.000	0.000
hsa-miR-2113	1.634	7.74	358	0.000	0.000
hsa-miR-2116-3p	1.656	9.782	115	0.000	0.000
hsa-miR-219a-2-3p	1.6	6.436	125	0.000	0.000
hsa-miR-2355-5p	2.762	9.902	113	0.000	0.000
hsa-miR-30c-1-3p	1.976	7.144	114	0.000	0.000
hsa-miR-3131	2.446	6.211	114	0.000	0.000
hsa-miR-3147	1.371	7.617	115	0.000	0.000
hsa-miR-3150a-3p	1.894	7.004	112	0.000	0.000
hsa-miR-3151-5p	1.387	6.831	115	0.000	0.000
hsa-miR-3153	1.994	9.741	115	0.000	0.000
hsa-miR-3178	4.752	8.999	97	0.000	0.000
hsa-miR-3180	4.161	9.232	92	0.000	0.000
hsa-miR-3184-5p	4.809	6.949	81	0.000	0.000
hsa-miR-3185	3.193	8.217	96	0.000	0.000
hsa-miR-3186-3p	1.446	6.611	104	0.000	0.000
hsa-miR-3197	4.454	13.341	95	0.000	0.000
hsa-miR-325	1.498	5.994	168	0.000	0.000
hsa-miR-330-5p	1.808	8.061	161	0.000	0.000
hsa-miR-412-3p	1.933	7.414	144	0.000	0.000
hsa-miR-423-3p	1.784	8.826	81	0.000	0.000
hsa-miR-4258	3.644	8.72	99	0.000	0.000
hsa-miR-4270	3.705	6.216	113	0.000	0.000
hsa-miR-4288	2.393	7.991	84	0.000	0.000
hsa-miR-4292	2.93	7.318	83	0.000	0.000
hsa-miR-4301	2.558	8.404	93	0.000	0.000
hsa-miR-4310	2.22	8.196	99	0.000	0.000
hsa-miR-4312	1.462	7.004	115	0.000	0.000
hsa-miR-4323	1.915	9.85	109	0.000	0.000
hsa-miR-512-5p	1.484	6.579	172	0.000	0.000
hsa-miR-515-3p	1.833	7.931	97	0.000	0.000
hsa-miR-517c-3p	1.635	6.22	186	0.000	0.000
hsa-miR-519d-3p	1.833	7.499	164	0.000	0.000
hsa-miR-520a-5p	1.627	6.71	151	0.000	0.000
hsa-miR-520d-5p	5.392	7.514	342	0.000	0.000
hsa-miR-520g-3p	2.047	6.687	119	0.000	0.000
hsa-miR-548a-5p	1.79	6.719	131	0.000	0.000
hsa-miR-548d-3p	2.938	6.573	62	0.000	0.000
hsa-miR-612	1.79	6.952	376	0.000	0.000
hsa-miR-615-5p	1.691	7.028	116	0.000	0.000
hsa-miR-625-3p	3.763	8.618	104	0.000	0.000
hsa-miR-637	2.081	7.404	159	0.000	0.000
hsa-miR-668-3p	1.437	6.341	153	0.000	0.000
hsa-miR-765	2.014	7.266	338	0.000	0.000
hsa-miR-891a-5p	2.297	13.332	105	0.000	0.000
hsa-miR-891b	2.052	7.342	117	0.000	0.000
hsa-miR-92b-5p	2.802	6.499	106	0.000	0.000
hsa-miR-938	1.556	8.445	155	0.000	0.000

The P value was set to lower than 0.001 (P < 0.001). miRNA: MicroRNA; FC: Fold change (ratio of mean signal intensities of cholangiocarcinoma to those of normal samples); DF: Degree of freedom; Abs. t value: Absolute t value; Adj. P value: Adjusted P value; FDR: False discovery rate.

Table 3

List of 10-fold-cross-validated differentially-expressed microRNAs which were down-regulated in cholangiocarcinoma

miRNA	FC	Abs. t value	DF	Adj. P value	FDR
hsa-let-7b-5p	0.336	6.655	186	0.000	0.000
hsa-let-7c-3p	0.686	9.519	138	0.000	0.000
hsa-let-7f-5p	0.226	6.65	89	0.000	0.000
hsa-miR-100-5p	0.424	7.984	144	0.000	0.000
hsa-miR-10a-5p	0.458	5.943	196	0.000	0.000
hsa-miR-10b-3p	0.352	8.834	72	0.000	0.000
hsa-miR-1225-5p	0.097	6.448	88	0.000	0.000
hsa-miR-127-3p	0.62	6.167	84	0.000	0.000
hsa-miR-1288-3p	0.394	7.17	61	0.000	0.000
hsa-miR-1305	0.311	7.255	66	0.000	0.000
hsa-miR-130a-3p	0.459	7.775	101	0.000	0.000
hsa-miR-136-5p	0.51	7.184	191	0.000	0.000
hsa-miR-139-3p	0.626	7.227	95	0.000	0.000
hsa-miR-145-5p	0.359	7.026	129	0.000	0.000
hsa-miR-150-3p	0.51	7.109	139	0.000	0.000
hsa-miR-17-3p	0.638	7.108	115	0.000	0.000
hsa-miR-181c-3p	0.602	6.651	133	0.000	0.000
hsa-miR-181d-5p	0.598	7.741	138	0.000	0.000
hsa-miR-192-5p	0.267	7.582	69	0.000	0.000
hsa-miR-194-5p	0.283	8.454	90	0.000	0.000
hsa-miR-195-5p	0.449	7.441	101	0.000	0.000
hsa-miR-20a-5p	0.386	6.929	99	0.000	0.000
hsa-miR-215-5p	0.271	8.404	70	0.000	0.000
hsa-miR-26b-5p	0.305	7.657	121	0.000	0.000
hsa-miR-27b-3p	0.317	8.02	146	0.000	0.000
hsa-miR-29b-3p	0.432	6.216	110	0.000	0.000
hsa-miR-29c-3p	0.25	9.578	93	0.000	0.000
hsa-miR-29c-5p	0.755	6.521	122	0.000	0.000
hsa-miR-3120-3p	0.673	6.341	90	0.000	0.000
hsa-miR-330-3p	0.691	6.605	99	0.000	0.000
hsa-miR-339-3p	0.746	7.662	90	0.000	0.000
hsa-miR-342-3p	0.516	8.228	193	0.000	0.000
hsa-miR-345-5p	0.531	8.154	69	0.000	0.000
hsa-miR-374a-5p	0.439	8.25	99	0.000	0.000
hsa-miR-374b-5p	0.568	6.771	117	0.000	0.000
hsa-miR-378b	0.489	7.181	101	0.000	0.000
hsa-miR-4294	0.471	7.279	96	0.000	0.000
hsa-miR-4326	0.69	6.513	115	0.000	0.000
hsa-miR-451a	0.143	7.636	85	0.000	0.000
hsa-miR-483-5p	0.412	7.272	166	0.000	0.000
hsa-miR-505-3p	0.625	6.612	105	0.000	0.000
hsa-miR-575	0.311	7.35	62	0.000	0.000
hsa-miR-614	0.341	7.8	412	0.000	0.000
hsa-miR-627-5p	0.601	6.504	190	0.000	0.000
hsa-miR-650	0.252	9.369	148	0.000	0.000
hsa-miR-671-5p	0.594	6.372	402	0.000	0.000
hsa-miR-922	0.323	8.522	148	0.000	0.000
hsa-miR-99a-5p	0.328	9.611	117	0.000	0.000

The P value was set to lower than 0.001 (P < 0.001). miRNA: MicroRNA; FC: Fold change (ratio of mean signal intensities of cholangiocarcinoma to those of normal samples); DF: Degree of freedom; Abs. t value: Absolute t value; Adj. P value: Adjusted P value; FDR: False discovery rate.

Figure 2

Heatmap of microRNA expression in cholangiocarcinoma and normal samples. The color gradient of each cell represents the log2 of normalized intensity value of microRNA microarray spot. CCA: Cholangiocarcinoma; miRNA: MicroRNA.

Enrichment analyses

To explore the biological functions of DE miRNAs found in the meta-analysis results, GO and KEGG pathway enrichment analyses were performed using DIANA miRPath v3.0 with P value threshold at 0.01 and MicroT score threshold at 0.8. The miRPath v.3.0 results indicated that there were 4407 and 7236 predicted target genes of up-regulated and down-regulated miRNAs, respectively, involved in GO biological processes. GO enrichment analysis showed 72 biological processes associated with up-regulated miRNAs (Supplementary Table 1), whereas 95 biological processes were identified to be associated with down-regulated miRNAs (Supplementary Table 2).

KEGG pathway enrichment analysis showed that gene targets of up-regulated miRNAs were significantly involved (P < 0.01) in 48 molecular biological processes (Supplementary Table 3) while gene targets of down-regulated miRNAs were significantly involved (P < 0.01) in 32 processes (Supplementary Table 4). Figure 3 represents the most prominent dysregulated pathways of up-regulated miRNAs including phosphatidylinositol-3 kinases/Akt (PI3K/Akt) signaling pathway (215 gene targets of 57 miRNAs, P = 0.00424), mitogen-activated protein kinase (MAPK) signaling pathway (169 gene targets of 55 miRNAs, P = 0.00034), and Ras signaling pathway (159 gene targets of 58 miRNAs, P = 2.34E-06). These three pathways were also the most prominent dysregulated pathways of down-regulated miRNAs, which were PI3K/Akt signaling pathway (209 gene targets of 42 miRNAs, P = 1.91E-05), MAPK signaling pathway (147 gene targets of 39 miRNAs, P = 0.0079), and Ras signaling pathway (147 gene targets of 39 miRNAs, P = 2.03E-06).

Figure 3

Pathway enrichment analyses of predicted target genes of differentially expressed microRNAs. The differentially expressed microRNAs obtained from meta-analysis were input to DIANA miRPath version 3.0. P value thresholds were set at 0.01. PI3K: Phosphatidylinositol-3 kinases; MAPK: Mitogen-activated protein kinase.

Visualization of miRNA-target interaction networks

According to the results of KEGG pathway analysis, PI3K/Akt, MAPK and Ras signaling pathways were the most enriched biological processes with which target genes of DE miRNAs were associated. For up-regulated miRNAs associated with these three pathways, 7 miRNAs including miR-330-5p, miR-519d-3p, miR-548a-5p, miR-548d-3p, miR-1207-3p, miR-1304-5p, and miR-2113 were chosen as representative miRNAs of this group. For down-regulated miRNAs, 9 miRNAs including let-7b-5p, let-7c-3p, let-7f-5p, miR-195-5p, miR-20a-5p, miR-26b-5p, miR-27b-3p, miR-29b-3p, and miR-330-3p were chosen as representatives. Using regulatory interaction networks data from TargetScan and miRTarBase, 35 and 46 unique target genes of up- and down-regulated miRNAs, respectively, were identified to be involved in these three pathways (Figure 4).

Figure 4

MicroRNA-target interaction networks. A: Up-regulated; B: Down-regulated. MicroRNA (miRNA)-target interaction networks of 7 up-regulated and 9 down-regulated miRNAs associating in phosphatidylinositol-3 kinases/Akt, mitogen-activated protein kinase, and Ras signaling pathways based on pathway enrichment analysis via DIANA miRPath version 3.0. Blue and red lines indicate the miRNA-target prediction or information based on TargetScan and miRTarBase databases, respectively. PI3K: Phosphatidylinositol-3 kinases; MAPK: Mitogen-activated protein kinase.

DISCUSSION

At present, a global view of miRNA roles in the development of cancers remains incomplete. Due to widespread usage of microarray technology, there has been an enormous expansion of publicly available datasets[37], which could be integrated and analyzed with a statistically robust meta-analysis approach. In our meta-analysis, miRNA microarray datasets from multiple independent studies were analyzed with highly stringent statistics and cross validation, leading to identification of DE miRNAs in CCA compared to non-cancerous cells across the pooled samples. The lack information on whether the samples were from primary or metastatic sites would pose as one of the limitations of this study.

Many DE miRNAs observed in our study are significantly related to 3 major cancer signaling pathways, namely the MAPK signaling pathway, PI3K/Akt signaling pathway, and Ras signaling pathway. Among the up-regulated miRNAs found in our meta-analysis, miR-519d and miR-330 are of particular interest as dysregulated expression of these miRNAs have been reported in several types of cancer. miR-519d belongs to the chromosome 19 miRNA cluster, which is the largest human miRNA cluster described so far[17]. Although there is no direct evidence on the relationship of miR-519d and CCA development, miR-519d has been shown to be up-regulated in HCC patient’s tissues, exerting oncogenic activity by inhibiting the tumor suppressor proteins such as CDKN1A/p21, PTEN, AKT3 and TIMP2[38]. In contrast, overexpression of miR-519d in a human HCC cell line QGY-7703 has been shown to block cell proliferation. On the contrary, down-regulation of miR-519d has been reported to promote cell proliferation in many cancers, including cervical cancer, breast cancer, and ovarian cancer[39-42].

Another up-regulated miRNA identified in our study includes miR-330. This miRNA has been previously reported to be up-regulated in glioblastoma[43], colorectal cancer[44], non-small cell lung cancer[45], and esophageal cancer[46] whereas its down-regulated expression has been demonstrated in prostate cancer[47-49]. miR-330 has been shown to promote cancer cell proliferation via suppression of CDC42, a Rho GTPase-associated with MAPK signaling pathway[44]. Besides, hypoxia-induced upregulation of integrin-alpha 5, a predicted target of miR-330 and a critical receptor in PI3K/Akt signaling pathway, has been shown to enhance cell proliferation, metastasis and apoptosis resistance of CCA cell lines[50-52]. Paradoxically, miR-330 has been reported to induce apoptosis in prostate cancer cells through E2F1-mediated suppression of Akt phosphorylation in prostate cancer cells[47].

One of the down-regulated miRNAs identified in our study was miR-20a, one of the mature miRNA products of the miR-17-92 cluster pri-miRNA. Down-regulation of miR-20a has been shown to mediate cellular differentiation and growth arrest induced by HIF-1 in acute myeloid leukemia cells by targeting p21 and STAT3[52], demonstrating an oncogenic role of miR-20a. In addition, the miR-17-92 cluster miRNAs have been shown to be up-regulated and play oncogenic roles in many types of cancer including CCA[53]. In contrast, down-regulation of miR-20a has been shown to promote HCC cell proliferation via upregulation of Mcl-1, an antiapoptotic member of the Bcl-2 family, suggesting a tumor-suppressor role[54].

Besides, our meta-analysis of CCA miRNA microarrays has revealed down-regulation of several members of the let-7 miRNA family (let-7b, -7c, -7f), whose members are estimated to comprise 1%-5% of the mammalian genome[10,17]. This miRNA family has been shown to generally play a tumor suppressor role[55], where ectopic expression of the let-7 family inhibits cell proliferation through down-regulation of c-Myc in nasopharyngeal carcinoma cells, whereas down-regulation of let-7 promotes cancer cell growth by increasing the activity of Ras protein, in lung cancer[56] and liver cancer[57] cells.

Altogether, the meta-analysis of miRNA microarray datasets with highly stringent statistical methodology provides new insights into the role of miRNA and its dysregulations in CCA. Our findings of miRNA dysregulations in the cancer signaling pathways including PI3K/Akt pathway, MAPK pathway, and Ras pathway give clues into underlying miRNA-mRNA interplays of CCA. However, the analyses reported herein were based on the different origins of miRNAs, validations of such findings are warranted. Of clinical relevance, since many miRNAs have been reported as highly specific biomarkers for several types of cancer[58-60], the identified miRNA in this study may have predictive values for CCA cases. Also, there are conflicting results in adjuvant settings for CCA[61], the detection of a specific miRNA may be associated with an increased risk of recurrence after surgery. Further investigation of the miRNAs reported herein will bring about the novel knowledge of the dysregulated processes in CCA development at post-transcriptional level which could offer novel diagnostic and therapeutic approaches in the future.

ARTICLE HIGHLIGHTS

Research background

The incidence of cholangiocarcinoma (CCA) is alarmingly elevating in many countries. Patients with CCA usually have poor prognosis as there is still no effective screening and treatment available. Therefore, it is essential to identify biomarkers for CCA.

Research motivation

Differential expression profiles of microRNA (miRNA) have been reported for many different types of cancer. Thus, a growing number of miRNA microarray data can be a valuable resource for the discovery of biomarkers to tackle challenges in the clinical management of CCA.

Research objectives

This work integrates and intervalidates the CCA miRNA expression profiles from multiple independent datasets to identify the differential dysregulation of miRNA and their corresponding downstream pathways underlying mechanism of pathogenesis.

Research methods

Eight independent CCA miRNA profiling microarray datasets, including 246 CCA and 197 normal samples were assimilated into a meta-analysis and cross-validation to identify a cohort of miRNA that were significantly dysregulated in CCA.

Research results

Of 118 dysregulated miRNA identified in our study, 70 were up-regulated and 48 were down-regulated miRNAs in CCA. Bioinformatic analyses revealed that mRNA targets of differentially expressed miRNAs were significantly distributed across various biological processes. The most prominent dysregulated pathways included phosphatidylinositol-3 kinases/Akt, mitogen-activated protein kinase and Ras signaling pathways.

Research conclusions

This current study represents the meta-analysis of miRNA microarray datasets with highly stringent statistical methodology and provides new insights into the role of miRNA and its dysregulations in CCA.

Research perspectives

The merit of our findings offers a valuable reference for future studies and further investigation of these miRNA/genes and their interactions will eventually lead to the identification of genes and pathways important to the overall mechanism of the dysregulated processes in CCA development.

ACKNOWLEDGEMENTS

We appreciate a scholarship from the Development and Promotion of Science and Technology Talented Project awarded to Likhitrattanapisal S.