Comprehensive analysis of differently expression mRNA and non-coding RNAs, and their regulatory mechanisms on relationship in thiram-induced tibial dyschondroplasia in chicken.

Thiram pollution is one of the main causes of tibial dyschondroplasia (TD) induced by feed sources. Several studies have speculated that miRNA, circRNA and lncRNA may have significant impact on the development of TD, however, the specific mRNAs and noncoding RNAs and their respective regulatory mechanisms and functions in the development of TD have not been explored. Therefore, in this present study, we screened the differentially expressed mRNA, miRNA, circRNA and lncRNA by whole-transcriptome sequencing (RNA-seq) and differentially expressed genes (DEGs) enrichment, as well as constructed the interaction network among the mRNA-miRNA, mRNA-lncRNA and mRNA-miRNA-circRNA. The sequencing results were verified by fluorescence real-time quantitative PCR (RT-qPCR). The results obtained in this study, revealed that the cells were atrophied and disordered in the TD group, and the expression of BMP6, TGF-β and VEGF were significantly reduced. A total of 141 mRNAs, 10 miRNAs, 23 lncRNAs and 35 circRNAs of DEGs were obtained (p<0.05) Theses DEGs were enriched in the adhere junction and insulin signaling pathways. In addition, the mRNA-miRNA-circRNA network suggested that several pivotal ceRNA showed a regulatory relationship between the transcripts with miRNA, circRNA or lncRNA. Taken together, the results in the present study, represent an insight for further functional research on the ceRNA regulatory mechanism of TD in broilers.