Deep integrative analysis of microRNA-mRNA regulatory networks for biomarker and target discovery in chondrosarcoma.

Chondrosarcoma (CHS) is a common malignant bone sarcoma and its occurrence increases with age. microRNAs (miRNAs) are a class of noncoding RNAs that participate in various biological processes and disease pathogenesis by targeting functional messenger RNA (mRNA). However, the modulation of miRNAs in CHS remains largely unknown. In this study, we performed integrative analysis to explore the expression profiles of miRNAs and mRNAs, together with their interaction networks in human CHS tissues and cell lines by RNA-seq (miRNA and mRNA). A total of 133 and 796 differentially expressed miRNAs and mRNAs were identified (|Fold change| ≥ 2 and P-value ≤ 0.5). miRNA-mRNA regulatory interactions between 55 miRNAs and 242 mRNAs were screened by the Pearson correlation analysis and target prediction. mRNAs in the network were enriched to 145 Gene Ontology terms and 35 Kyoto Encyclopedia of Genes and Genomes pathways. Specifically, some key regulators (hub-miRNAs) in the network (miR-622, miR-4539, miR-145, miR-25, and miR-96) were suggested to play important regulatory roles in the pathogenesis of CHS. In addition, functional experiments validated that miR-622 regulated CHS cell proliferation by targeting bone morphogenetic protein 1 (BMP1).