Analysis of the Long Non-Coding and Messenger RNA Expression Profiles in the Skin Tissue of Super Merino and Small-Tailed Han Sheep

Abstract

Wool quality and yield are two important economic livestock traits. However, there are relatively few molecular studies on lncRNA for improving sheep wool, so these require further exploration. In this study, we examined skin tissue from the upper scapula of Super Merino (SM) and Small-Tailed Han (STH) sheep during the growing period. The apparent difference was verified via histological examination. High-throughput RNA sequencing identified differentially expressed (DE) long non-coding (lncRNAs) and messenger RNAs (mRNAs). The target gene of DE lncRNA and DE genes were enrichment analyzed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). A Reverse Transcription quantitative Polymerase Chain Reaction (RT-qPCR) was used to verify randomly selected DE lncRNAs and mRNAs. Finally, the DE, RAC2, WNT11, and FZD2 genes, which were enriched in the Wnt signaling pathway, were detected via immunohistochemistry. The results showed that a total of 20,888 lncRNAs and 31,579 mRNAs were identified in the skin tissues of the two sheep species. Among these, 56 lncRNAs and 616 mRNAs were differentially expressed. Through qRT-PCR, the trends in the randomly selected DE genes’ expression were confirmed to be aligned with the RNA-seq results. GO and KEGG enrichment analysis showed that DE lncRNA target genes were enriched in GO terms as represented by epidermal and skin development and keratin filature and in KEGG terms as represented by PI3K-Akt, Ras, MAPK, and Wnt signaling pathways, which were related to hair follicle growth and development. Finally, immunohistochemistry staining results indicated that RAC2, WNT11, and FZD2 were expressed in dermal papilla (DP). The lncRNAs MSTRG.9225.1 and MSTRG.98769.1 may indirectly participate in the regulation of hair follicle growth, development, and fiber traits by regulating their respective target genes, LOC114113396(KRTAP15-1), FGF1, and IGF1. In addition, MSTRG.84658.1 may regulate the Wnt signaling pathway involved in the development of sheep hair follicles by targeting RAC2. This study provides a theoretical reference for improving sheep breeding in the future and lays a foundation for further research on the effects of MSTRG.84658.1 and the target gene RAC2 on dermal papilla cells (DPC).

Keywords

• super merino sheep
• small-tailed han sheep
• lncRNA
• mRNA
• RNA sequencing
• immunohistochemistry