Comprehensive analysis of coding and non-coding RNA transcriptomes related to hypoxic adaptation in Tibetan chickens

Abstract Background Tibetan chickens, a unique native breed in the Qinghai-Tibet Plateau of China, possess a suite of adaptive features that enable them to tolerate the high-altitude hypoxic environment. Increasing evidence suggests that long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) play roles in the hypoxic adaptation of high-altitude animals, although their exact involvement remains unclear. Results This study aimed to elucidate the global landscape of mRNAs, lncRNAs, and miRNAs using transcriptome sequencing to construct a regulatory network of competing endogenous RNAs (ceRNAs) and thus provide insights into the hypoxic adaptation of Tibetan chicken embryos. In total, 354 differentially expressed genes (DE genes), 389 differentially expressed lncRNAs (DE lncRNAs), and 73 differentially expressed miRNAs (DE miRNAs) were identified between Tibetan chickens (TC) and control Chahua chickens (CH). GO and KEGG enrichment analysis revealed that several important DE miRNAs and their target DE lncRNAs and DE genes are involved in angiogenesis (including blood vessel development and blood circulation) and energy metabolism (including glucose, carbohydrate, and lipid metabolism). The ceRNA network was then constructed with the predicted DE gene-DE miRNA-DE lncRNA interactions, which further revealed the regulatory roles of these differentially expressed RNAs during hypoxic adaptation of Tibetan chickens. Conclusions Analysis of transcriptomic data revealed several key candidate ceRNAs that may play high-priority roles in the hypoxic adaptation of Tibetan chickens by regulating angiogenesis and energy metabolism. These results provide insights into the molecular mechanisms of hypoxic adaptation regulatory networks from the perspective of coding and non-coding RNAs.