Abstract
BACKGROUND: Cashmere has long been used as the raw material for wool textiles. The diameter of the cashmere fibre determines its quality and economic value. However, the regulatory role of noncoding RNAs (ncRNAs) in cashmere fineness remains unclear, especially regarding the interaction between ncRNAs and coding RNAs. RESULTS: Transcriptome sequencing was used to identify the expression profiles of long noncoding RNAs (lncRNAs), circular RNAs (circRNAs) and microRNAs (miRNAs) in the skin tissues of Jiangnan cashmere goats with different cashmere fineness levels. Integration analysis of ncRNA and coding RNA was performed in combination with previous research results. The results showed that 16,437 lncRNAs, 2234 circRNAs, and 1322 miRNAs were identified in 8 skin samples of cashmere goats. A total of 403 differentially expressed (DE) lncRNAs, 62 DE circRNAs and 30 DE miRNAs were identified in the skin tissues of the fine groups (Fe) and coarse groups (Ce). We predicted the target gene of DE lncRNA, the target gene of DE miRNA and the host gene of DE circRNA. Based on functional annotation and enrichment analysis of target genes, we found that DE lncRNAs could be involved in regulating the fineness traits of cashmere. The most potential lncRNAs were MSTRG.42054.1, MSTRG.18602.3, and MSTRG.2199.13. CONCLUSIONS: The data from this study enriched the cashmere goat noncoding RNA database and helped to supplement the annotation of the goat genome. The results provided a new direction for the breeding of cashmere characters.