As one of the most obvious phenotypic traits, the coat color of sheep is an ideal model to study the genetic mechanisms underlying coat color varieties of mammals. One distinguishable coat color is the black-headed type, such as the famous black-headed Dorper sheep from Africa and Bayinbuluke sheep from Asia. In this study, we compared the genome sequences of black-headed and all-white sheep to identify causative genes for the black-headed sheep, including black-headed Dorper versus white-headed Dorper, as well as Bayinbuluke (black-headed) versus Small-tailed Han (all-white). The most differentiating region between black-headed sheep and all-white sheep was found to harbor a haplotype covering melanocortin receptor 1 (MC1R) gene. The share of this haplotype by the black-headed sheep from Africa and Asia suggested that the convergent change in the MC1R region is likely to determine this unique coat color. Two missense mutations (g. 14251947T>A and g. 14252090G>A) within this haplotype of MC1R gene were found. We further analyzed whole genome sequence data of 460 world-wide sheep with diverse coat colors and confirmed the association between the MC1R haplotype with pigmentation variations. Our study provides novel insights into coat color genetics in sheep and expands our knowledge of the link between MC1R gene and varying pigmentation patterns in sheep.

Despite recent advances in generating high-quality reference genome assemblies, the genome sequences for most livestock species, including goats, are still poorly annotated. Single-molecule long-read sequencing has greatly facilitated gene annotation by obtaining full-length transcripts. In this study, we generated full-length transcriptome data for samples from abomasum (n = 2) and testicle (n = 1), using PacBio Iso-Seq technology. We further combined these data with published data from abomasum (5ZY, SRR8618141) to evaluate and improve the gene annotation of the goat genome. We identified 14.5-16.3% of novel genes per sample from the four Iso-Seq datasets. At the transcript level, 40.6% of them were novel, including 29.7% novel transcripts from known genes and 10.9% from novel genes. We further verified the expression of novel genes in four additional RNA-seq data and found that the expression level of novel genes was significantly lower than that of known genes, indicating that the lowly expressed genes tend to be missed in the current genome annotation. This study shows the superiority of full-length transcriptome data in gene annotation, and more such data are required to improve the gene annotation for goat genome and other species.

Marker-assisted selection is an important method for livestock breeding. In recent years, this technology has been gradually applied to livestock breeding to improve the body conformation traits. In this study, the LRRC8B (Leucine Rich Repeat Containing 8 VRAC Subunit B) gene was selected to evaluate the association between its genetic variations and the body conformation traits in two native sheep breeds in China. Four body conformation traits, including withers height, body length, chest circumference, and body weight, were collected from 269 Chaka sheep. We also collected the body length, chest width, withers height, chest depth, chest circumference, cannon bone circumference, and height at hip cross of 149 Small-Tailed Han sheep. Two different genotypes, ID and DD, were detected in all sheep. Our data showed that the polymorphism of the LRRC8B gene was significantly associated with chest depth (p < 0.05) in Small-Tailed Han sheep, and it is greater in sheep with DD than those with ID. In conclusion, our data suggested that the LRRC8B gene could serve as a candidate gene for marker-assisted selection in Small-Tailed Han sheep.