Abstract:[Objective] Tibetan sheep is one of three coarse wool sheep breeds in China,distributing mainly in Tibet and its adjacent alpine pastoral areas,such as Qinghai,Gansu,Sichuan,Yunnan.The wool of Tibetan sheep,felt smooth and warm duo to the long and thick fibers,is good raw material for manufacture of Tibetan carpet.Coat color is an important economic traits.However,at present,the mechanism of molecular regulation of sheep coat color is not clear.In this study,the transcriptome of skin tissues of different color (black and white) was sequenced with a aim to explore the role of miRNAs at post-transcriptional levels in different coat color skin tissues and possible regulatory pathways.[Method] Healthy sheep with black and white coat color was sacrificed to provide skin tissues for RNA extraction and miRNA analysis.Through miRNA sequencing and bioinformatics analysis,miRNA expression profiles of skin tissues of black and white wool color were obtained.Then the differentially expressed miRNAs were screened and the related target genes were predicted.Through the target gene analysis,the signal pathway which is related to Tibetan coat color skin traits were proposed.[Result] A total of 85.76 million original reads and 85.08 million clean reads were obtained from the analyzed tissues.Among the clean reads were 334 known miRNAs and 59 newly identified miRNAs and 23 of them were differentially expressed between white and black color tissues.In the 23 differentially expressed miRNAs,14 and 9 were up-regulated and down-regulated in white skin tissues respectively.miR-2284b and miR-744 were the type that is expressed only in white Tibetan sheep skin,and miR-23b,miR-41 la-5p,miR-30c,miR-423-3p,miR-324-5p were the miRNAs that expressed only in black skin.Among them,miR-10a can participate in the regulation of a variety of signaling pathways,involving proliferation,differentiation,apoptosis,cell adhesion.MiR-23b is associated with Wnt,Notch and other signal pathways,which may related to melanin production.MiR-41 la-5p,miR103 and miR-200b may also closely relate to regulation of hair color.A total of 981 target genes were predicted to be involved in coat color control,most of which were related to WNT signaling pathway,MAPK pathway,EDNRB signaling pathway and Melanogenesis pathway.The melanogenesis pathway was found to be associated with the WNT signal pathway,the KIT signal pathway,and the EDNRB signal pathway.The synthesis of melanin may be regulated by downstream TYR gene (tyrosinase gene) through MITF gene.Quantitative analysis showed that five miRNAs were up-regulated in black skin tissue,and two miRNAs were up-regulated in white tissues.The quantitative results were consistent with the sequencing results.[Conclusion] The miRNA expression profile of the skin of Tibetan sheep were obtained and the miRNAs and signal pathways which may be related to control of coat color were obtained by bioinformatics analysis.These results indicated that the color traits could be regulated by many miRNAs and involved many signal pathways.This may help understand the level of post-transcriptional regulation of the process,and lay the foundation for further functional verification.

Profiles of Mirnas and Target Gene Analysis with White and Black Skin Tissues of the Tibetan Sheep