Integration Analysis of circRNA–miRNA–mRNA and Identification of Critical Networks in Valgus-Varus Deformity (Gallus gallus)

Valgus-valgus deformity (VVD) is a common leg deformity in broilers with inward or outward deviation of the tibiotarsus and tarsometatarsus. The competing endogenous RNA (ceRNA) network plays an essential role in the study of leg disease. However, its role in the etiology and pathogenesis of VVD remains unclear. Here, based on case (VVD) and control (normal) group design, we performed analyses of differentially expressed circRNAs (DEcircRNAs), differentially expressed miRNAs (DEmiRNAs) and differentially expressed mRNAs (DEmRNAs). Transcriptome data derived 86 DEcircRNAs, 13 DEmiRNAs and 410 DEmRNAs. Functional analysis showed that DEmRNAs were significantly enriched in cell cycle, apoptosis, ECM-receptor interaction, FoxO signaling pathway and protein processing synthesis. DEcirc/miRNA-associated DEmRNAs were associated with skeletal and muscle growth and development pathways, including mTOR, Wnt, and VEGF signaling pathways. Subsequently, a circRNA–miRNA–mRNA regulatory network was constructed based on the ceRNA hypothesis, including 8 circRNAs, 6 miRNAs, and 31 mRNAs, which were significantly enriched in the skeletal developmental pathway. Finally, two key mRNAs (CDC20 and CTNNB1) and their regulatory axes were screened by the PPI network and cytohubba. The expression levels of CDC20 and CTNNB1 in cartilage and seven other tissues were also quantified by qPCR. In conclusion, we analyzed the functions of DEmRNA, DEcircRNA and DEmiRNA and constructed the hub ceRNA regulatory axis, and obtained two hub genes, CDC20 and CTNNB1. The study more deeply explored the etiology and pathogenesis of VVD and lays the foundation for further study of the role of the ceRNA network on skeletal development.