Abstract
AIM: This study aimed to identify key genes, non-coding RNAs, and their possible regulatory interactions during gallbladder cancer (GBC). BACKGROUND: The early detection of GBC, i.e. before metastasis, is restricted by our limited knowledge of molecular markers and mechanism(s) involved during carcinogenesis. Therefore, identifying important disease-associated transcriptome-level alterations can be of clinical importance. METHODS: In this study, six NCBI-GEO microarray dataseries of GBC and control tissue samples were analyzed to identify differentially expressed genes (DEGs) and non-coding RNAs {microRNAs (DEmiRNAs) and long non-coding RNAs (DElncRNAs)} with a computational meta-analysis approach. A series of bioinformatic methods were applied to enrich functional pathways, create protein-protein interaction networks, identify hub genes, and screen potential targets of DEmiRNAs and DElncRNAs. Expression and interaction data were consolidated to reveal putative DElncRNAs:DEmiRNAs:DEGs interactions. RESULTS: In total, 351 DEGs (185 downregulated, 166 upregulated), 787 DEmiRNAs (299 downregulated, 488 upregulated), and 7436 DElncRNAs (3127 downregulated, 4309 upregulated) were identified. Eight genes (FGF, CDK1, RPN2, SEC61A1, SOX2, CALR, NGFR, and NCAM) were identified as hub genes. Genes associated with ubiquitin ligase activity, N-linked glycosylation, and blood coagulation were upregulated, while those for cell-cell adhesion, cell differentiation, and surface receptor-linked signaling were downregulated. DEGs-DEmiRNAs-DElncRNAs interaction network identified 46 DElncRNAs to be associated with 28 DEmiRNAs, consecutively regulating 27 DEGs. DEmiRNAs-hsa-miR-26b-5p and hsa-miR-335-5p; and DElnRNAs-LINC00657 and CTB-89H12.4 regulated the highest number of DEGs and DEmiRNAs, respectively. CONCLUSION: The current study has identified meaningful transcriptome-level changes and gene-miRNA-lncRNA interactions during GBC and laid a platform for future studies on novel prognostic and diagnostic markers in GBC.