Identification of characteristic genes and construction of regulatory network in gallbladder carcinoma.

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. 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). 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. 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.

Shao H ，Zhu J ，Zhu Y ，Liu L ，Zhao S ，Kang Q ，Liu Y ，Zou H ... -  《BMC Medical Genomics》

Identification of cell cycle as the critical pathway modulated by exosome-derived microRNAs in gallbladder carcinoma.

Gallbladder cancer (GBC), the most common malignancy in the biliary tract, is highly lethal malignant due to seldomly specific symptoms in the early stage of GBC. This study aimed to identify exosome-derived miRNAs mediated competing endogenous RNAs (ceRNA) participant in GBC tumorigenesis. A total of 159 differentially expressed miRNAs (DEMs) was identified as exosome-derived miRNAs, contains 34 upregulated exo-DEMs and 125 downregulated exo-DEMs based on the expression profiles in GBC clinical samples downloaded from the Gene Expression Omnibus database with the R package. Among them, 2 up-regulated exo-DEMs, hsa-miR-125a-3p and hsa-miR-4647, and 5 down-regulated exo-DEMs, including hsa-miR-29c-5p, hsa-miR-145a-5p, hsa-miR-192-5p, hsa-miR-194-5p, and hsa-miR-338-3p, were associated with the survival of GBC patients. Results of the gene set enrichment analysis showed that the cell cycle-related pathways were activated in GBC tumor tissues, mainly including cell cycle, M phase, and cell cycle checkpoints. Furthermore, the dysregulated ceRNA network was constructed based on the lncRNA-miRNA-mRNA interactions using miRDB, TargetScan, miRTarBase, miRcode, and starBase v2.0., consisting of 27 lncRNAs, 6 prognostic exo-DEMs, and 176 mRNAs. Together with prognostic exo-DEMs, the STEAP3-AS1/hsa-miR-192-5p/MAD2L1 axis was identified, suggesting lncRNA STEAP3-AS1, might as a sponge of exosome-derived hsa-miR-192-5p, modulates cell cycle progression via affecting MAD2L1 expression in GBC tumorigenesis. In addition, the biological functions of genes in the ceRNA network were also annotated by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes. Our study promotes exploration of the molecular mechanisms associated with tumorigenesis and provide potential targets for GBC diagnosis and treatment.

Su L ，Zhang J ，Zhang X ，Zheng L ，Zhu Z ... -  《-》

Identification of Serum Exosome-Derived circRNA-miRNA-TF-mRNA Regulatory Network in Postmenopausal Osteoporosis Using Bioinformatics Analysis and Validation in Peripheral Blood-Derived Mononuclear Cells.

Osteoporosis is one of the most common systemic metabolic bone diseases, especially in postmenopausal women. Circular RNA (circRNA) has been implicated in various human diseases. However, the potential role of circRNAs in postmenopausal osteoporosis (PMOP) remains largely unknown. The study aims to identify potential biomarkers and further understand the mechanism of PMOP by constructing a circRNA-associated ceRNA network. The PMOP-related datasets GSE161361, GSE64433, and GSE56116 were downloaded from the Gene Expression Omnibus (GEO) database and were used to obtain differentially expressed genes (DEGs). Gene ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were applied to determine possible relevant functions of differentially expressed messenger RNAs (mRNAs). The TRRUST database was used to predict differential transcription factor (TF)-mRNA regulatory pairs. Afterwards, combined CircBank and miRTarBase, circRNA-miRNA as well as miRNA-TF pairs were constructed. Then, a circRNA-miRNA-TF-mRNA network was established. Next, the correlation of mRNAs, TFs, and PMOP was verified by the Comparative Toxicogenomics Database. And expression levels of key genes, including circRNAs, miRNAs, TFs, and mRNAs in the ceRNA network were further validated by quantitative real-time PCR (qRT-PCR). Furthermore, to screen out signaling pathways related to key mRNAs of the ceRNA network, Gene Set Enrichment Analysis (GSEA) was performed. A total of 1201 DE mRNAs, 44 DE miRNAs, and 1613 DE circRNAs associated with PMOP were obtained. GO function annotation showed DE mRNAs were mainly related to inflammatory responses. KEGG analysis revealed DE mRNAs were mainly enriched in osteoclast differentiation, rheumatoid arthritis, hematopoietic cell lineage, and cytokine-cytokine receptor interaction pathways. We first identified 26 TFs and their target mRNAs. Combining DE miRNAs, miRNA-TF/mRNA pairs were obtained. Combining DE circRNAs, we constructed the ceRNA network contained 6 circRNAs, 4 miRNAs, 4 TFs, and 12 mRNAs. The expression levels of most genes detected by qRT-PCR were generally consistent with the microarray results. Combined with the qRT-PCR validation results, we eventually identified the ceRNA network that contained 4 circRNAs, 3 miRNAs, 3 TFs, and 9 mRNAs. The GSEA revealed that 9 mRNAs participate in many important signaling pathways, such as "olfactory transduction", "T cell receptor signaling pathway", and "neuroactive ligand-receptor interaction". These pathways have been reported to the occurrence and development of PMOP. To sum up, key mRNAs in the ceRNA network may participate in the development of osteoporosis by regulating related signal pathways. A circRNA-associated ceRNA network containing TFs was established for PMOP. The study may help further explore the molecular mechanisms and may serve as potential biomarkers or therapeutic targets for PMOP.

Dong Q ，Han Z ，Tian L 《Frontiers in Endocrinology》

Complex integrated analysis of lncRNAs-miRNAs-mRNAs in oral squamous cell carcinoma.

This study aims to reveal regulatory network of lncRNAs-miRNAs-mRNAs in oral squamous cell carcinoma (OSCC) through gene expression data. Differentially expressed lncRNAs, miRNAs and mRNAs (cut-off: False discovery rate (FDR)<0.05 and |fold change|>1.5) were unveiled by package edgeR of R. Cox regression analysis was performed to screen prognostic factors in OSCC related with overall survival (OS) and relapse-free survival (RFS). Protein-protein interaction (PPI) network was constructed for differentially expressed mRNAs using BioGRID, HPRD and DIP. Key hub genes were identified from top 100 differentially expressed mRNAs ranked by betweenness centrality using recursive feature elimination. LncRNA-miRNA and miRNA-mRNA regulatory network were constructed and combined into ceRNAs regulatory network. Gene ontology biological terms and Kyoto Encyclopedia of Genes and Genomes pathways were identified using Fisher's exact test. A total of 929 differentially expressed mRNAs, 23 differentially expressed lncRNAs and 29 differentially expressed miRNAs were identified. 59 mRNAs, 6 miRNAs (hsa-mir-133a-1, hsa-mir-1-2, hsa-mir-486, hsa-mir-135b, hsa-mir-196b, hsa-mir-193b) and 6 lncRNAs (C10orf91, C2orf48, SFTA1P, FLJ41941,PART1,TTTY14) were related with OS; and 52 mRNAs, 4 miRNAs (hsa-mir-133a-1, hsa-mir-135b, hsa-mir-196b, hsa-mir-193b) and 2 lncRNAs (PART1, TTTY14) were associated with RFS. A support vector machine (SVM) classifier containing 37 key hub genes was obtained. A ceRNA regulatory network containing 417 nodes and 696 edges was constructed. ECM-receptor interaction, cytokine-cytokine receptor interaction, focal adhesion, arachidonic acid metabolism, and p53 signaling pathway were significantly enriched in the network. These findings uncover the pathogenesis of OSCC and might provide potential therapeutic targets.

Li S ，Chen X ，Liu X ，Yu Y ，Pan H ，Haak R ，Schmidt J ，Ziebolz D ，Schmalz G ... -  《-》

Microarray data analysis on gene and miRNA expression to identify biomarkers in non-small cell lung cancer.

Jin X ，Guan Y ，Zhang Z ，Wang H ... -  《BMC CANCER》