In the present study, we overexpressed miR-449a-5p in a rat model of SCI and used BBB scores to test the recovery and movement of rat hind limbs, as well as bioinformatics to analyze the potential targets of miR-449a-5p on axon regeneration after SCI. The 76 potential gene targets of miR-449a-5p were obtained from the databases, and were performed with GO enrichment analysis, and KEGG signaling pathways analysis, and PPI network. In GO enrichment analysis, regulation of ion transmembrane transport is the most abundant BP, and synaptic membrane is the most important chemical component in CC. Phospholipase D Signaling Pathway is the most important pathway in KEGG signaling pathway analysis, with the most significant differential expression, and we believe that miR-449a-5p may regulate the rehabilitation of axonal regeneration after SCI through this signaling pathway. PPI network analysis showed that CNTN2, ANK3 and CNTNAP1 were hug genes of miR-499a-5p.
SCI is a relatively common motor system injury in clinical practice, and its incidence is increasing year by year (Zhou, Wang et al. 2016, Eckert and Martin 2017). SCI results in varying degrees of paralysis or loss of work capacity in patients and are characterized by high morbidity, disability, and cost. SCI include primary and secondary injuries. Primary SCI is irreversible, while secondary SCI develops gradually after primary injuries and is often accompanied by axon demyelination, inflammation, apoptosis and glial scar formation, changes in gene expression and cellular metabolism, thus secondary SCI causes more severe damages to patients than primary SCI (Ratan and Noble 2009, Oyinbo 2011, Anwar, Al Shehabi et al. 2016, Cizkova, Murgoci et al. 2020).
As a kind of short non-coding single-stranded small molecule RNA in organisms, miRNA can regulate the expression of at least 30% of human protein-coding genes. More than 1000 miRNAs related to human beings have been identified, providing a new direction for the diagnosis, treatment and mechanism research of various clinical diseases (Shi, Zhou et al. 2017). MiR-449 family is located in the intron of cell division cycle gene 20B (CDC20B) on chromosome 5. This region has been identified as a strong cancer susceptibility site in earlier genome-wide association studies (Sucheston-Campbell, Clay-Gilmour et al. 2018). MiR-449a, located in the first intron of CDC20B on chromosome 5Q11, acts as a tumor suppressor by regulating cell proliferation, apoptosis, migration and invasion in a variety of cancers (Kumar, Sharad et al. 2016).
It was shown that TPD52 is a direct target of miR-449, and the expression of miR-449 was down-regulated and the expression of TPD52 was up-regulated in breast cancer cell line MDAMB-231, and it was speculated that the down-regulation of miR-449 may promote the migration and invasion of breast cancer cells by eliminating the inhibition of TPD52 (Zhang, Wang et al. 2017). Additionally, miR-449a acts as a tumor suppressor gene by suppressing the expression of CAPN6 or POU2F1 in hepatoma carcinomacells, inhibiting the proliferation of cancer cells, and promoting the apoptosis of HCC cells (Liu, Wang et al. 2016). Down-regulation of miR-449a-5p promotes the proliferation of esophageal squamous cell carcinoma cells, and silencing Bcl-2, a target gene of miR-449a-5p, reversed the effect of miR-449a-5p inhibitors on the proliferation of ECA-109 cells (Jiang, Zhao et al. 2020). It was shown that inhibition of miR-449a promotes cartilage regeneration and prevents osteoarthritis in an acute defect modelin vivo,intra-articular injection of LNA-anti-miR-449a prevents disease progression in osteoarthritis models (Baek, Lee et al. 2018). MiR-449a promotes neuron-like differentiation by targeting histone deacetylase 1 (HDAC1) expression (Liu, Hou et al. 2015). In this study, the rats exhibited the paralysis of hind limb and lost mobility after SCI modeling. Additionally, the comparison of BBB scores at 1d, 7d, and 21d after SCI surgery suggested that the overexpression of miR-499a-5p improved the recovery of motor function of rats with SCI. It is indicated that miR-449a-5p could induce or inhibit the repair process after SCI and change the rehabilitation process of SCI.
Biological function involved by gene targets of miR-499a-5p
In this study, the target genes of miR-449a-5p were predicted and GO analysis of target genes was carried out by bioinformatics. Protein localization to axonal and synaptic tissues were significantly enriched in BP related to axon regeneration. The axons of developing neurons produce an active structure called growth cone that guides their normal physiological processes. The growth cone regulates the growth, steering and localization of axons through continuous polymerization and depolymerization of microtubules and microfilaments. Tau is a microtubule-associated protein that is thought to be localized to axons. Moreover, Tau axon localization does not require stable microtubule binding. Tau lacks the microtubule binding domain (MTBD) and shows high diffusivity, but correctly localizes to axons (Iwata, Watanabe et al. 2019). The results of this study suggests that miR-449a may act mainly by interfering in the process of protein axonal localization.
Hub genes of miR-499a-5p
In this study, PPI network analysis of the predicted intersected gene targets of miR-499a-5p showed that HDAC1, Notch1 and Kitlg were centrally located and connected to many other nodes, and the hub genes CNTN2, ANK3 and CNTNAP1 were obtained from PPI network and GO enrichment analyses. CNTN2 is expressed in CNS and oligodendrocytes. CNTN2 in neurons is highly regulated during development and plays an important role in axon growth, guidance and neuronal migration. On the other hand, CNTN2 expressed in oligodendrocytes interferes with myelination, and its ablation leads to hypomyelination. It was found that CNTN2 in neurons may function through inhibition of Akt signaling pathway (Savvaki, Kafetzis et al. 2021). Also, the biological effect of miR-3075 on the phenotype of Schwann cells may be realized through the negative regulation of CNTN2, which regulates cell migration and myelination (Wang, He et al. 2018). The role of the ankyrin-G protein encoded by the ANK3 gene is to bind the skeleton protein to the cell membrane, and it plays an important role in many biological activities as an important linker protein. In recent years, it has been found that ankyrin-G proteins may be involved in signal transduction and apoptosis, and are involved in many key processes of signal transduction. The ANK3 gene,encoding multiple ankyrin-G subtypes, is widely expressed in various tissues, was first identified in the initial segment of axons and Ranferian node (Kordeli, Lambert et al. 1995), and is speculated to be involved in the process of axon reconstruction. CNTNAP1, also known as Caspr1, is a member of the Caspr family, and is involved in neural excitation and conduction and releases neurotransmitters in myelinated axons (Gollan, Salomon et al. 2003). Caspr1 is not only involved in the formation of myelinated axons, but also involved in maintenance the stability of adjacent connections, the formation, differentiation and proliferation of neurons and astrocytes, as well as motor control and cognitive functions (Sherman, Tait et al. 2005, Buttermore, Dupree et al. 2011). These results suggests that CNTN2, ANK3 and CNTNAP1 play an important role in neural repair during axonal regeneration and development after SCI.