High-throughput transcriptome analysis reveals potentially important relationships between lncRNAs and genes in broilers affected by Valgus-varus Deformity (Gallus gallus).

Valgus-varus Deformity (VVD) is an outward or inward deviation of the tibiotarsus or tarsometatarsus, which results in physical distress of chickens and economic loss in poultry industry. While the etiology and pathogenesis of VVD at the molecular level are still not fully understood so far. Here, based on a case/control design with VVD birds and normal birds, we identified genes and lncRNAs which associated with VVD using RNA sequencing. Transcriptome analysis revealed 231 differentially expressed mRNAs and 23 differentially expressed lncRNAs between case and control of leg cartilage. We identified the cis- and trans-regulatory targets of the differentially expressed lncRNAs, and we constructed a functional lncRNA-mRNA co-expression network. Analysis of the network showed that the differentially expressed mRNAs and the target genes of the differentially expressed lncRNAs were enriched in the signaling pathways associated with bone development, including p53, MAPK, Toll-like receptor, Jak-STAT, Hedgehog, and PPAR. The expression levels of DENND4A, FGF10, FGF12 and BMP3 were also determined in cartilage and other six tissues. Overall, our study predicted the mRNAs and lncRNAs related with leg diseases by transcriptome analyses, which might contribute to understand the etiology and pathogenesis of VVD. It established the foundation for the further research on the function of -mRNAs and lncRNAs in skeleton development.