Transcriptome analysis reveals the spinal expression profiles of non-coding RNAs involved in anorectal malformations in rat fetuses.

BACKGROUND: Despite improvements in anorectal malformation (ARM) therapy, patients might still experience post-operative problems such as fecal incontinence, constipation, and soiling. In particular, the dysplasia of the lumbosacral spinal cord in ARM patients is a major disorder that affects fecal function post-operation. However, the pathological mechanisms involved are still unclear. METHODS: The non-coding RNAs (ncRNAs) in the lumbosacral spinal cord of fetal rats with ethylenethiourea-induced ARM were identified using RNA sequencing (RNA-seq) and examined to determine their potential function. The lumbosacral spinal cord was isolated on embryonic day 17 for subsequent RNA extraction and RNA-seq. The transcriptome data was analyzed using bioinformatics analysis, followed by validation using quantitative reverse transcription PCR. RESULTS: Compared to the control group, 26 differentially expressed microRNAs (DEMs; 22 upregulated, 4 downregulated) and 112 differentially expressed long non-coding RNAs (63 upregulated, 49 downregulated) were identified in the ARM group. Several DEMs related to development, namely miR-200a-3p, miR-200b-3p, miR-200c-3p, miR-200a-5p, and miR-429, were selected for further analysis. Notably, compared to the control, the relative expression of miR-200 family members was highly upregulated in ARM fetal rats. Furthermore, GO and KEGG enrichment and miRNA-transcription factor-lncRNA/mRNA network analysis was explored to show molecular mechanism underlying DEMs. CONCLUSIONS: Our findings suggest the involvement of ncRNAs, especially the miR-200 family members, in the pathogenesis of lumbosacral spinal cord dysplasia in ARM fetal rats.