Multi-tissue transcriptome profiling linked the association between tissue-specific circRNAs and the heterosis for feed intake and efficiency in chicken

Heterosis has been widely utilized in chickens. The nonadditive inheritance of genes contributes to this biological phenomenon. However, the role of circRNAs played in the heterosis is poorly determined. In this study, we observed divergent heterosis for residual feed intake (RFI) between two crossbreds derived from a reciprocal cross between White Leghorns and Beijing You chickens. Then, circRNA landscape for 120 samples covering the hypothalamus, liver, duodenum mucosa and ovary were profiled to elucidate the regulatory mechanisms of heterosis. We detected that a small proportion of circRNAs (7.83%∼20.35%) were additively and non-additively expressed, in which non-additivity was a major inheritance of circRNAs in the crossbreds. Tissue-specific expression of circRNAs was prevalent across 4 tissues. Weighted gene co-expression network analysis revealed circRNA-mRNA co-expression modules associated with feed intake and RFI in the hypothalamus and liver, and the co-expressed genes were enriched in oxidative phosphorylation pathway. We further identified 8 nonadditive circRNAs highly correlated with 16 nonadditive genes regulating negative heterosis for RFI in the two tissues. Circ-ITSN2 was validated in the liver tissue for its significantly positive correlation with PGPEP1L. Moreover, the bioinformatic analysis indicated that candidate circRNAs might be functioned by binding the microRNAs and interacting with the RNA binding proteins. The integration of multi-tissue transcriptome firstly linked the association between tissue-specific circRNAs and the heterosis for feed intake and efficiency in chicken, which provide novel insights into the molecular mechanism underlying heterosis for feed efficiency. The validated circRNAs can act as potential biomarkers for predicting RFI and its heterosis.