LncRNAs are essential for regulating skeletal muscle. However, the expression profile and function of lncRNAs in goat muscle remains unclear. Here, an average of ~14.58 Gb high-quality reads were obtained from longissimus dorsi tissues of 1-month-old (n = 3) and 9-month-old (n = 3) Wu'an black goats using RNA sequencing. Of a total of 3441 lncRNAs, 1281 were lincRNAs, 805 were antisense lncRNAs, and 1355 were sense_overlapping lncRNAs. These lncRNAs shared some properties with goats, such as fewer exons, shorter transcript, and open reading frames (ORFs) length. Among them, 36 differentially expressed lncRNAs (DE lncRNA) were identified, and then 10 random lncRNAs were validated by RT-qPCR. Furthermore, 30 DE lncRNAs were neighboring 71 mRNAs and several genes were functionally enriched in muscle development-related pathways, such as APC, IFRD1, NKX2-5, and others. Additionally, 36 DE lncRNAs and 2684 mRNAs were included in co-expression interactions. A lncRNA-miRNA-mRNA network containing 4 lncRNAs, 3 miRNAs, and 8 mRNAs was finally constructed, of which XR_001296113.2 might regulate PDLIM7 expression by interaction with chi-miR-1296 to affect skeletal muscle development. This study revealed the expression profile of goat lncRNAs for further investigative studies and provides a fuller understanding of skeletal muscle development.
Keywords: RNA sequencing; goat; lncRNAs; regulation networks; skeletal muscle.