Identification of long noncoding RNAs involved in adaptability to chronic hypoxic by whole transcriptome sequencing

Hypoxia affects the physiology of cells and organisms; however, the mechanisms associated with hypoxia adaptation remain unknown in Tibetan chickens. In this study, we aimed to identify long noncoding RNAs (lncRNAs) involved in hypoxia adaptation in Tibetan chickens and Daheng broilers, to provide insights into the mechanisms underlying hypoxia induction. RNA sequencing results revealed that a total of 5504 lncRNAs and 16,779 microRNAs were differentially expressed in four Tibetan chickens and four Daheng broilers; 70 lncRNAs were up-regulated and 113 lncRNAs were down-regulated in the Tibetan chickens compared to the expression levels in the Daheng broilers. The differentially expressed lncRNAs (DElncRNAs) were enriched in the following Gene ontology terms: protein complex localization, small-molecule metabolic process, and RNA splicing. Kyoto Encyclopedia of Genes and Genomes analyses revealed that the DElncRNAs were mainly enriched in pathways that regulate cell junctions and intercellular spaces and oxygen or energy metabolism, mainly involved in hypoxic adaption. Moreover, a predicted ceRNA network with five DElncRNAs interacted with three miRNAs that acted on 42 pathways through 19 target genes. Quantitative real-time polymerase chain reaction was used to verify that the expression levels of ENSGALG00000008047, ENSGALG00000050044, and ENSGALG00000053982 were significantly lower in Tibetan chickens than in the Daheng broilers, consistent with the RNA sequencing results. We obtained lncRNA expression profiles for the heart tissue of Tibetan chickens for the first time and have provided novel data that may aid research on biological adaptation to hypoxic stress.