miRNA signature identification of retinoblastoma and the correlations between differentially expressed miRNAs during retinoblastoma progression.

PURPOSE: Retinoblastoma (RB) is a common pediatric cancer. The study aimed to uncover the mechanisms of RB progression and identify novel therapeutic biomarkers.

METHODS: The miRNA expression profile GSE7072, which includes three RB samples and three healthy retina samples, was used. After data normalization using the preprocessCore package, differentially expressed miRNAs (DE-miRs) were selected by the limma package. The targets of the DE-miRs were predicted based on two databases, followed by construction of the miRNA-target network. Pathway enrichment analysis was conducted for the targets of the DE-miRNAs using DAVID. The CTD database was used to predict RB-related genes, followed by clustering analysis using the pvclust package. The correlation network of DE-miRs was established. MiRNA expression was validated in another data set, GSE41321.

RESULTS: In total, 24 DE-miRs were identified whose targets were correlated with the cell cycle pathway. Among them, hsa-miR-373, hsa-miR-125b, and hsa-miR-181a were highlighted in the miRNA-target regulatory network; 14 DE-miRs, including hsa-miR-373, hsa-miR-125b, hsa-miR-18a, hsa-miR-25, hsa-miR-20a, and hsa-let-7 (a, b, c), were shown to distinguish RB from healthy tissue. In addition, hsa-miR-25, hsa-miR-18a, and hsa-miR-20a shared the common target BCL2L11; hsa-let-7b and hsa-miR-125b targeted the genes CDC25A, CDK6, and LIN28A. Expression of three miRNAs in GSE41321 was consistent with that in GSE7072.

CONCLUSIONS: Several critical miRNAs were identified in RB progression. Hsa-miR-373 might regulate RB invasion and metastasis, hsa-miR-181a might involve in the CDKN1B-mediated cell cycle pathway, and hsa-miR-125b and hsa-let-7b might serve as tumor suppressors by coregulating CDK6, CDC25A, and LIN28A. The miRNAs hsa-miR-25, hsa-miR-18a, and hsa-miR-20a might exert their function by coregulating BCL2L1.