Identification of Potential miRNA-mRNA Regulatory Network in Denervated Muscular Atrophy by Bioinformatic Analysis

Muscle atrophy caused by long-term denervation leads to the loss of skeletal muscle mass and strength, resulting in a poor recovery of functional muscles and decreasing quality of life. Increasing differentially expressed microRNAs (DEMs) have been reported to be involved in the pathogenesis of denervated muscle atrophy. However, there is still insufficient evidence to explain the role of miRNAs and their target genes in skeletal muscle atrophy. Therefore, an integrative exploration of the miRNA-mRNA regulatory network in denervated muscle atrophy is necessary. A total of 21 (16 upregulated and 5 downregulated) DEMs were screened out in the GSE81914 dataset. Med1, Myod1, Nfkb1, Rela, and Camta1 were predicted and verified to be significantly upregulated in denervated muscle atrophy, from which 6 key TF-miRNA relationship pairs, including Med1-mir-1949, Med1-mir-146b, Myod1-mir-29b, Nfkb1-mir-21, Rela-mir-21, and Camta1-mir-132, were obtained. 60 target genes were then predicted by submitting candidate DEMs to the miRNet database. GO and KEGG pathway enrichment analysis showed that target genes of DEMs were mainly enriched in the apoptotic process and PI3K/Akt signaling pathway. Through the PPI network construction, key modules and hub genes were obtained and potentially modulated by mir-29b, mir-132, and mir-133a. According to the qRT-PCR results, the expression of COL1A1 and Ctgf is opposite to their related miRNAs in denervated muscle atrophy. In the study, a potential miRNA-mRNA regulatory network was firstly constructed in denervated muscle atrophy, in which the mir-29b-COL1A1 and mir-133a-Ctgf pathways may provide new insights into the pathogenesis and treatment.