Genome-Wide Analysis Reveals Extensive Changes in LncRNAs during Skeletal Muscle Development in Hu Sheep

As an important type of noncoding RNA molecules, long non-coding RNAs (lncRNAs) act as versatile players in various biological processes. However, little is known about lncRNA regulators during sheep muscle growth. To explore functional lncRNAs during sheep muscle growth, we systematically investigated lncRNAs using strand-specific Ribo-Zero RNA sequencing at three key developmental stages in Hu sheep. A total of 6924 lncRNAs were obtained, and the differentially expressed lncRNAs and genes were screened from (control vs. experiment) fetus vs. lamb, lamb vs. adult, and fetus vs. adult comparisons, respectively. The quantitative real-time polymerase chain reaction (qRT-PCR) analysis results correlated well with the sequencing data. Moreover, functional annotation analysis based on the Gene Ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) databases showed that the target genes of the differentially expressed lncRNAs were significantly enriched in organ morphogenesis, skeletal system development as well as response to stimulus and some other terms related to muscle. Furthermore, a co-expression network of the differentially expressed target genes and lncRNAs was constructed and well-known muscle growth regulators such as retrotransposon-like 1 and Junctophilin-2 were included. Finally, we investigated the expression profiles of seven lncRNAs and their target genes, and found that they played vital roles in muscle growth. This study extends the sheep muscle lncRNA database and provides novel candidate regulators for future genetic and molecular studies on sheep muscle growth, which is helpful for optimizing the production of mutton.