Identification of Characteristic Genes and Construction of Regulatory Network in Gallbladder Carcinoma

Background: Gallbladder carcinoma (GBC) is a highly malignant tumor with a poor overall prognosis. This study aimed to identify the characteristic microRNAs (miRNAs) of GBC and the competing endogenous RNA (ceRNA) regulatory mechanisms.

Methods: The microarray data of GBC tissue samples and normal gallbladder (NGB) tissue samples from the Gene Expression Omnibus (GEO) database was downloaded. GBC-related differentially expressed miRNAs (DE-miRNAs) were identified by inter-group differential expression analysis and weighted gene co-expression network analysis (WGCNA). Machine learning algorithms were used to screen the characteristic miRNA based on the intersect between least absolute shrinkage and selection operator (LASSO) and Support vector machine-recursive feature elimination (SVM-RFE). Based on the differential expression analysis of GEO database, the ceRNA network of characteristic miRNA was predicted and constructed. The biological functions of the ceRNA network were revealed by carrying out the gene enrichment analysis was implemented. We further screened the key genes of ceRNA network and constructed a protein-protein interaction (PPI) network, and predicted and generated the transcription factors (TFs) network of signature miRNAs. The expression of characteristic miRNA in clinical samples was verified by quantitative real-time polymerase chain reaction (qRT-PCR).

Results: A total of 131 GBC-related DE-miRNAs were obtained. The hsa-miR-4770 was defined as characteristic miRNA for GBC. The ceRNA network containing 211 mRNAs, one miRNA, two lncRNAs, and 48 circRNAs was created. Gene enrichment analysis suggested that the downstream genes were mainly involved in actin filament organization, cell-substrate adhesion, cell-matrix adhesion, reactive oxygen species metabolic process, glutamine metabolic process and extracellular matrix (ECM)-receptor interaction pathway. 10 key genes in the network were found to be most correlated with disease, and involved in cell cycle-related processes, p53, and extrinsic apoptotic signaling pathways. The qRT-PCR result demonstrated that hsa-miR-4770 is down-regulated in GBC, and the expression trend is consistent with the public database.

Conclusions: We identified hsa-miR-4770 as the characteristic miRNA for GBC. The ceRNA network of hsa-miR-4770 may play key roles in GBC. This study provided some basis for potential pathogenesis of GBC.