Identification and Characterization of Circular RNAs (circRNAs) Using RNA-Seq in Two Breeds of Cashmere Goats

Abstract

Circular RNA (circRNA) is a type of non-coding RNA generated from back-splicing the reactions of linear RNA. It plays an important role in various cellular and biological processes. However, there are few studies about the regulatory effect of circRNAs on cashmere fiber traits in cashmere goats. In this study, the expression profiles of circRNAs in skin tissue were compared between Liaoning cashmere (LC) goats and Ziwuling black (ZB) goats, with a significant difference in cashmere fiber yield, cashmere fiber diameter, and cashmere fiber color, using RNA-seq. A total of 11,613 circRNAs were expressed in the caprine skin tissue, and their type, chromosomal distribution, and length distribution were characterized. A total of 115 up-regulated circRNAs and 146 down-regulated circRNAs in LC goats were screened compared to ZB goats. The authenticity of 10 differentially expressed circRNAs was validated by detecting their expression levels and the head-to-tail splice junction using RT-PCR and DNA sequencing, respectively. The parent genes of differentially expressed circRNA were mainly enriched in some Gene Ontology (GO) terms and pathways related to cashmere fiber traits, such as the canonical Wnt signaling pathway, which is involved in the regulation of cell promotion, stem cell proliferation, Wnt signaling pathway regulation, epithelial morphogenesis, MAPK signaling pathway, and cell adhesion molecules pathway. Eight differentially expressed circRNAs were further selected to construct a circRNA-miRNA network, and some miRNAs that were previously reported as related to fiber traits were found in the network. This study provides a deep understanding of the roles of circRNAs in the regulation of cashmere fiber traits in cashmere goats and the involvement of differential splicing in phenotypic expression according to breed and region.