Identification of Competing Endogenous RNA Regulatory Networks in Vitamin A Deficiency-Induced Congenital Scoliosis by Transcriptome Sequencing Analysis

Abstract

Background/Aims: Congenital scoliosis (CS) is a result of anomalous development of vertebrae and is frequently associated with somitogenesis malformation. Although noncoding RNAs (ncRNAs) have been recently determined to be involved in the pathogenesis of CS, the competing endogenous RNA (ceRNA) regulatory networks in CS remain largely unknown. Methods: Sequencing was conducted to explore the ncRNA expression profiles in rat embryos (gestation day 9) following vitamin A deficiency (VAD) (n = 9 for the vitamin A deficiency-induced congenital scoliosis (VAD-CS) group and n = 4 for the control group). Real-time reverse transcriptase polymerase chain reaction (RT-PCR) was conducted to verify the expression levels of selected mRNAs, long noncoding RNAs (lncRNAs), circular RNAs (circRNAs), and microRNAs (miRNAs). Bioinformatics analysis was used to discover the possible relationships and functions of the ceRNAs. Results: A total of 749 mRNAs, 56 miRNAs, 685 lncRNAs, and 70 circRNAs were identified to have significantly different expression levels in the two groups. Wnt, PI3K-ATK, FoxO, EGFR, and mTOR were found to be the most significant pathways involved in VAD-CS pathogenesis. The circRNA/miRNA/mRNA and lncRNA/miRNA/mRNA networks of CS were built, and the gene expression mechanisms regulated by ncRNAs were unveiled via the ceRNA regulatory networks. Conclusion: We comprehensively identified ceRNA regulatory networks of embryonic somite development in VAD-CS as well as revealed the contribution of different ncRNA expression profiles. Our data demonstrate the association between mRNAs and ncRNAs in the pathogenic mechanism of CS.